Consistent arbitrary Lagrangian Eulerian formulation for large deformation thermo-mechanical analysis

Arbitrary Lagrangian Eulerian (ALE) method is widely used for simulation of large deformation problems, such as metal forming. However, in many such applications, modeling of the heat generation and transfer in conjunction with the stress analysis is necessary. In this work, a fully coupled dynamic ALE formulation is developed. The ALE form of energy balance equation is derived, and is coupled with the dynamic, rate dependent ALE stress analysis. The proposed formulation is used for simulation of a few thermo-mechanical problems. The effectiveness and efficiency of the ALE method is verified by comparing the results of this simulation with available experimental and numerical results.     

An extended BEM formulation for plates reinforced by rectangular beams

This article presents a BEM formulation developed particularly for analysis of plates reinforced by rectangular beams. This is an extended version of a previous paper that only took into account bending effects. The problem is now re-formulated to consider bending and membrane force effects. The effects of the reinforcements are taken into account by using a simplified scheme that requires application of an initial stress field to locally correct the bending and stretching stiffness of the reinforcement regions. The domain integrals due to the presence of the reinforcements are then transformed to the reinforcement/plate interface. To reduce the number of degrees of freedom related to the presence of the reinforcement, the proposed model was simplified to consider only bending and stretching rigidities in the direction of the beams. The complete model can be recovered by applying all six internal force correctors, corresponding to six degrees of freedom per node. Examples are presented to confirm the accuracy of the formulation and to illustrate the level of simplification introduced by this strong reduction in the number of degrees of freedom.     

An enhanced co‐rotational approach for large displacement analysis of plates

This paper presents a new co‐rotational approach for the large displacement analysis of plates employing 4‐noded quadrilateral flat shell elements. The proposed approach benefits from (i) a simple local co‐rotational system invariant to the element nodal ordering, (ii) the choice of the two smallest components of the nodal normal vector as global rotational degrees of freedom, and (iii) the use of hierarchic freedoms, that are unaffected by the co‐rotational transformations, for higher‐order accuracy. Important additional benefits that arise from the aforementioned features include symmetry of the tangent stiffness matrix and complete insensitivity of the large displacement transformations to the size of the incremental step. The applicability of the new approach to moderately thick as well as thin plates is illustrated by considering two alternative local formulations based on the Reissner–Mindlin and discrete Kirchhoff hypotheses. Several examples are finally presented which demonstrate the accuracy, step‐insensitivity and computational benefits of the proposed co‐rotational approach for large displacement analysis of plate structures. Copyright © 2005 John Wiley & Sons, Ltd.     

On the displacement formulation of torsion of shafts with rectangular cross-sections

The use of warping displacement functions for the torsional stiffness representation of beams with rectangular cross-sections is studied. These functions can directly be employed in the displacement-based formulation of Hermitian and isoparametric three-dimensional beam elements for linear, elastic-plastic or large displacement/large rotation analysis. the results of some studies are given to demonstrate the applicability and use of the proposed functions.     

Large reference displacement analysis of composite plates part I: Finite element formulation

This investigation concerns itself with the dynamic analysis of thin, laminated composite plates consisting of layers of orthotropic laminae that undergo large arbitrary rigid body displacements and small elastic deformations. A non-linear finite element formulation is developed which utilizes the assumption that the bonds between the laminae are infinitesimally thin and shear non-deformable. Using the expressions for the kinetic and strain energies, the lamina mass and stiffness matrices are identified. The non-linear mass matrix of the lamina is expressed in terms of a set of invariants that depend on the assumed displacement field. By summing the kinetic and strain energies of the laminae of an element, the element mass and stiffness matrix can be defined in terms of the set of element invariants. It is shown that the element invariants can be expressed explicitly in terms of the invariants of its laminae. By assembling the finite elements of the deformable body, the body invariants can be identified and expressed explicitly in terms of the invariants of the laminae of its elements. In the dynamic formulation presented in this paper, the shape functions of the laminae are assumed to have rigid body modes that need to describe only large rigid body translations. The computer implementation and the use of the formulation developed in this investigation in multibody dynamics are discussed in the second part of this paper.     

Effect of the coupling between stretching and bending in the large displacement analysis of plates

The effect of the rotary inertia on the non-linear dynamics of plates that undergo a large reference displacement is examined in this paper. An assumed displacement field that accounts for the coupling between the stretching and bending of the plate as the result of considering the effect of the rotary inertia is used to identify the configuration of the plate. Furthermore, the coupling between the stretching and bending of the plate as the result of finite rotation is also considered in this investigation. Based on the assumed displacement field that accounts for the effect of the rotary inertia, a non-linear finite element formulation is developed for the large displacement analysis of plates. The element equations of motion are expressed in terms of a set of element invariants that depend on the assumed displacement field as well as the rotary inertia. The use of the formulation presented in this paper is demonstrated using numerical examples.     

A mortar formulation for 3D large deformation contact using NURBS-based isogeometric analysis and the augmented Lagrangian method

NURBS-based isogeometric analysis is applied to 3D frictionless large deformation contact problems. The contact constraints are treated with a mortar-based approach combined with a simplified integration method avoiding segmentation of the contact surfaces, and the discretization of the continuum is performed with arbitrary order NURBS and Lagrange polynomial elements. The contact constraints are satisfied exactly with the augmented Lagrangian formulation proposed by Alart and Curnier, whereby a Newton-like solution scheme is applied to solve the saddle point problem simultaneously for displacements and Lagrange multipliers. The numerical examples show that the proposed contact formulation in conjunction with the NURBS discretization delivers accurate and robust predictions. In both small and large deformation cases, the quality of the contact pressures is shown to improve significantly over that achieved with Lagrange discretizations. In large deformation and large sliding examples, the NURBS discretization provides an improved smoothness of the traction history curves. In both cases, increasing the order of the discretization is either beneficial or not influential when using isogeometric analysis, whereas it affects results negatively for Lagrange discretizations.     

The large deflection dynamic plastic response of rectangular plates

By assuming a kinematically admissible, time-dependent velocity field and considering the global equilibrium of all the forces acting on each rigid segment during large deflection of a plate, a complete theoretical investigation is undertaken herein to trace the large deflection dynamic plastic response (including the transient phase) of simply-supported or fully-clamped rectangular plates. This procedure deduces two modifying factors ƒ₁ and ƒ₂ which reflect the effect of the membrane forces and are employed to formulate the governing equations in the case of large deflections. The present prediction of the final plate deflection coincides excellently with the corresponding experimental results for deflections up to 5–10 times the plate thickness. This theory greatly improves the estimates based on the bending-only theory and provides a new way to trace the transient phase of dynamically loaded plates when the effect of membrane forces is significant.     

Elastic large deflection analysis of plates subjected to uniaxial thrust using meshfree Mindlin-Reissner formulation

A meshfree approach for plate buckling/post-buckling problems in the case of uniaxial thrust is presented. A geometrical nonlinear formulation is employed using reproducing kernel approximation and stabilized conforming nodal integration. The bending components are represented by Mindlin–Reissner plate theory. The formulation has a locking-free property in imposing the Kirchhoff mode reproducing condition. In addition, in-plane deformation components are approximated by reproducing kernels. The deformation components are coupled to solve the general plate bending problem with geometrical non-linearity. In buckling/post-buckling analysis of plates, the in-plane displacement of the edges in their perpendicular directions is assumed to be uniform by considering the continuity of plating, and periodic boundary conditions are considered in assuming the periodicity of structures. In such boundary condition enforcements, some node displacements/rotations should be synchronized with others. However, the enforcements introduce difficulties in the meshfree approach because the reproducing kernel function does not have the so-called Kronecker delta property. In this paper, the multiple point constraint technique is introduced to treat such boundary conditions as well as the essential boundary conditions. Numerical studies are performed to examine the accuracy of the multiple point constraint enforcements. As numerical examples, buckling/post-buckling analyses of a rectangular plate and stiffened plate structure are presented to validate the proposed approach.     

Lagrangian formulation of the Maxwell-Cattaneo hydrodynamics

The Maxwell-Cattaneo modification of the Fourier constitutive relation allows to formulate hydrodynamics as dynamics of a material and a caloric point. The Maxwell-Cattaneo hydrodynamics is compatible with equilibrium thermodynamics, it represents an approximation of a kinetic equation, and it reduces to the standard hydrodynamics if the new transport coefficient (a relaxation time) tends to zero.     

Boundary element formulation for the coupled stretching-bending analysis of thin laminated plates

The boundary integral equations for the coupled stretching-bending analysis of thin laminated plates involve an integral which will be singular when the field point approaches the source point. To avoid the singular problem occurring in the numerical programming, the boundary integral equations are modified in which the integrals of singular part are integrated analytically. The analytical solutions for the free term coefficients and singular integrals are obtained in explicit closed-form. By dividing the boundary into elements and using suitable interpolation polynomials for basic functions, the set of equations necessary for boundary element programming are written explicitly for regular nodes and corner nodes. The equations for the determination of displacements and stresses at internal points are also presented in this paper.     

BEM formulation for reinforced plates

This article presents a BEM formulation developed to analyse reinforced plate bending. The reinforcements are formulated using a simplified scheme based on applying an initial moment field adopted to locally correct the stiffness of the reinforcement regions. The domain integrals due to the presence of the reinforcements are then transformed to the reinforcement/plate interface. The increase in system stiffness due to the reinforcements can be taken into account independently for each coefficient. Thus, one can conveniently reduce the number of degrees of freedom required in considering the reinforcement. Only one degree-of-freedom is required at each internal node when taking into account only the flexural stiffness of beams. Examples are presented to confirm the accuracy of the formulation.     

Vibration analysis of orthotropic rectangular plates on elastic foundations

The vibration responses of orthotropic plates on nonlinear elastic foundations are numerically modeled using the differential quadrature method. The differential quadrature technique is utilized to transform partial differential equations into a discrete eigenvalue problem. Numerical results and those from literature closely correspond to each other. Numerical results demonstrate that elastically restrained stiffness, plate aspect ratio and foundation stiffness significantly impact the dynamic behavior of orthotropic plates.     

Boundary element frequency domain formulation for dynamic analysis of Mindlin plates

A new boundary element formulation for analysis of shear deformable plates subjected to dynamic loading is presented. Fundamental solutions for the Mindlin plate theory are derived in the Laplace transform domain. The characteristics of the three flextural waves are studied in the time domain. It is shown that the new fundamental solutions exhibit the same strong singularity as in the static case. Two numerical examples are presented to demonstrate the accuracy of the boundary element method and comparisons are made with the finite element method. Copyright © 2006 John Wiley & Sons, Ltd.     

A BEM formulation based on Reissner’s theory to perform simple bending analysis of plates reinforced by rectangular beams

In this work, the plate bending formulation of the boundary element method (BEM), based on the Reissner’s hypothesis, is extended to the analysis of plates reinforced by rectangular beams. This composed structure is modelled by a zoned plate, being the beams represented by narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each sub-region, and summing them to get the equation for the whole plate. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. In order to decrease the number of degrees of freedom, some approximations are considered for both displacements and tractions along the beam width. The accuracy of the proposed model is illustrated by simple examples whose exact solution are known as well as by more complex examples whose numerical results are compared with a well-known finite element code.     

Asymptotic post-buckling analysis of rectangular plates by HC finite elements

The present paper extends the finite element perturbation approach already presented for pin-jointed and framed structures¹⁵ to rectangular thin plates. Koiter's asymptotic strategy² is coupled with a High-Continuity finite element discretization of the plate.²² The consistency of the discrete model is discussed from the kinematical and numerical points of view and several numerical tests are reported. It appears that use of the HC elements makes the perturbation algorithm insensitive to the locking phenomenon occurring in the evaluation of the postbuckling behaviour. It also allows the use of very fine discretization meshes at low computational cost.     

Vibration analysis for piezoceramic rectangular plates using Ritz's method with equivalent constants

The transverse vibration of piezoceramic rectangular thin plates is investigated theoretically and experimentally using the Ritz's method incorporated with the defined equivalent constants. The equivalent constants are derived by comparing the characteristic equations of transverse resonant frequencies between isotropic and piezoceramic disks. By replacing the Poisson's ratio and flexural rigidity with the equivalent constants, the well-known Ritz's method can be used to investigate the transverse vibration of piezoceramic rectangular plates. Two different types of boundary conditions - clamped-free-free-free (CFFF) and clamped-free-clamped-free (CFCF) - are analyzed in this paper. For the experimental measurement, two optical techniques - amplitude-fluctuation electronic speckle pattern interferometry (AF-ESPI) and laser Doppler vibrometer (LDV) - are used to validate the analytical results. Both the transverse vibration modes and resonant frequencies of piezoceramic rectangular plates are obtained by the AF-ESPI method. Numerical calculations using the finite-element method (FEM) are performed, and the results are compared with the theoretical analysis and experimental measurements. Excellent agreements are obtained for results of both resonant frequencies and mode shapes. According to the theoretical calculations with different equivalent Poisson's ratios, resonant frequency variations versus aspect ratios ranging from 0.1 to 10 also are discussed for the first several modes in the work.     

A boundary element formulation for analysis of elastoplastic plates with geometrical nonlinearity

In this paper a new boundary element method formulation for elastoplastic analysis of plates with geometrical nonlinearities is presented. The von Mises criterion with linear isotropic hardening is considered to evaluate the plastic zone. Large deflections are assumed but within the context of small strain. To derive the boundary integral equations the von Kármán’s hypothesis is taken into account. An initial stress field is applied to correct the true stresses according to the adopted criterion. Isoparametric linear elements are used to approximate the boundary unknown values while triangular internal cells with linear shape function are adopted to evaluate the domain value influences. The nonlinear system of equations is solved by using an implicit scheme together with the consistent tangent operator derived along the paper. Numerical examples are presented to demonstrate the accuracy and the validity of the proposed formulation.     

Elastic-plastic buckling analysis of rectangular plates subjected to biaxial loads

In order to calculate the buckling load of a rectangular plate, the analytical approach is used in this study. The plate is assumed to be simply supported on four edges and loaded by uniform stresses along the edges. If the plate is slender, the buckling is elastic. However, if the plate is sturdy, it buckles in the plastic range. Then, the instantaneous moduli in the constitutive equations depend on the external loading. In this study, the elastic and plastic buckling equations are derived for rectangular plates under biaxial loading, and the corresponding interaction curves are presented. The influences of aspect ratios, load ratios and hardening factors on the buckling stresses are investigated for rectangular plates. From the plastic buckling analysis, the optimal combination of loads is given for the buckling strength.