Total Lagrangian formulation for the large displacement analysis of rectangular plates

In this paper, a method for the non-linear dynamic analysis of rectangular plates that undergo large rigid body motions and small elastic deformations is presented. The large rigid body displacement of the plate is defined by the translation and rotation of a selected plate reference. The small elastic deformation of the midplane is defined in the plate co-ordinate system using the assumptions of the classical theories of plates. Non-linear terms that represent the dynamic coupling between the rigid body displacement and the elastic deformation are presented in a closed form in terms of a set of time-invariant scalars and matrices that depend on the assumed displacement field of the plate. In this paper, the case of simple two-parameter screw displacement, where the rigid body translation and rotation of the plate reference are, respectively, along and about an axis fixed in space, is first considered. The non-linear dynamic equations that govern the most general and arbitrary motion of the plate are also presented and both lumped and consistent mass formulations are discussed. The non-linear dynamic formulation presented in this paper can be used to develop a total Lagrangian finite element formulation for plates in multibody systems consisting of interconnected structural elements.     

Combined effect of shear and large deformation on stress intensity factors

An attempt has been made to study the influence of large deformation on the stress intensity factor in a cracked plate subjected to bending including shear deformation. It is assumed that singular terms for stress resultants and strains in the case of large deformation have the same angular distribution and order of singularity as in the case of a linear problem. With this in view the small deformation singular element has been used at the crack tip region surrounded by large deformation plate bending elements. The finite element analysis, based on total Lagrangian formulation combined with the modified Newton–Raphson technique, has been used to get numerical results. Several examples connected with large deformation of cracked plates subjected to bending are studied. Using the above technique stress intensity factors for linear and non-linear cases have been compared.     

A new shear flexible cubic spline plate element for vibration analysis

Here, a new cubic B‐spline plate element is developed using field consistency principle, for vibration analysis. The formulation includes anisotropy, transverse shear deformation, in‐plane and rotary inertia effects. The element is based on a laminated refined plate theory, which satisfies the interface transverse shear stress and displacement continuity, and has a vanishing shear stress on the top and bottom surfaces of the plates. The lack of consistency in the shear strain field interpolations in its constrained physical limits produces poor convergence and results in unacceptable solutions due to locking phenomenon. Hence, numerical experimentation for the evaluation of natural frequencies of plates is carried out to check this deficiency with a series of assumed shear strain functions, redistributed in a field consistent manner. Copyright © 2002 John Wiley & Sons, Ltd.     

Three-dimensional absolute nodal co-ordinate formulation: plate problem

In this investigation, an absolute nodal co-ordinate dynamic formulation is developed for the large deformations and rotations of three-dimensional plate elements. In this formulation, no infinitesimal or finite rotations are used as nodal co-ordinates, instead global displacements and slopes are used as the plate coordinates. Using this interpretation of the plate coordinates the new method does not require the use of co-ordinate transformation to define the global inertia properties of the plates. The resulting mass matrix is the same constant matrix that appears in linear structural dynamics. The stiffness matrix, on the other hand, is a non-linear function of the nodal co-ordinates of the plate even in the case of a linear elastic problem. It is demonstrated in this paper that, unlike the incremental finite element formulations, the proposed method leads to an exact modelling of the rigid body inertia when the plate element moves as a rigid body. It is also demonstrated that by using the proposed method the conventional plate element shape function has a complete set of rigid body modes that can describe an exact arbitrary rigid body displacement. Using this fact, plate elements in the proposed new formulation can be considered as isoparametric elements. As a consequence, an arbitrary rigid body motion of the element results in zero strain. © 1997 John Wiley & Sons, Ltd.     

Dynamic analysis of shear deformable plates using the Dual Reciprocity Method

The Dual Reciprocity Method is a popular mathematical technique to treat domain integrals in the boundary element method (BEM). This technique has been used to treat inertial integrals in the dynamic thin plate bending analysis using a direct formulation of the BEM based on the elastostatic fundamental solution of the problem. In this work, this approach was applied for the dynamic analysis of shear deformable plates based on the Reissner plate bending theory, considering the rotary inertia of the plate. Three kinds of problems: modal, harmonic and transient dynamic analysis, were analyzed. Numerical examples are presented to demonstrate the efficiency and accuracy of the proposed formulation.     

Geometrically non-linear analysis of symmetrically laminated composite and sandwich shells with a higher-order theory and C finite elements

A C⁰ continuous displacement based finite element formulation of a higher order theory for linear and geometrically non-linear analysis which accounts for large displacements in the sense of von Karman of symmetrically laminated composite and sandwich shells under transverse loads is presented. The displacement model accounts for non-linear and constant variation of tangential and transverse displacement components, respectively, through the shell thickness. The assumed displacement model climinates the use of shear correction coefficients. The discrete element chosen is a nine-node quadrilateral element with nine degress of freedom per node. The accuracy of the present formulation is then established by comparing the present results with the available analytical. closed-form two-dimensional solutions, three-dimensional elasticity solutions and other finite element solutions. Some new results are generated for future comparisons to and evaluations of sandwich shells.     

A solid-shell layerwise finite element for non-linear geometric and material analysis

The application of layerwise theories to correctly model the displacement field of sandwich structures or laminates with high modulus ratios usually employs plate or facet-shell finite element formulations to compute the element stiffness and mass matrices for each layer. In this work an alternative approach is proposed, using a high performance hexahedral finite element to represent the individual layer mass and stiffness. This eight-node hexahedral finite element is formulated based on the application of the enhanced assumed strain method (EAS) to solve several locking pathologies coming from the high aspect ratio of the finite element and the usual incompressibility condition of the core materials. The solid-shell finite element formulation is introduced in the layerwise theory through the definition of a projection operator, based on the finite element variables transformation matrix. The non-linear geometric and material capabilities are introduced into the finite element formulation, allowing for the representation of large displacements, large deformation and material non-linear behaviors. The developed formulation is numerically tested and benchmarked, being validated by using published experimental results obtained from sandwich specimens.     

A BEM formulation based on Reissner's hypothesis for analysing the coupled stretching–bending problem of building floor structures

In this work, a plate bending formulation of the boundary element method (BEM) based on the Reissner's hypothesis to perform linear analysis of plates reinforced by rectangular beams is extended to consider the beams not displayed over their middle surface. Therefore eccentricity effects are taken into account. The building floor structure is modelled as a stiffened plate which is treated as a single body without dividing it into beam and plate elements. Moreover the equilibrium and compatibility conditions are automatically imposed by the integral equations. In the proposed model the final system of equation is obtained by coupling the bending problem to the stretching problem. Besides, in order to reduce the number of degrees of freedom, both the displacements and tractions are approximated along the beam width, leading to a model where the values are defined on the beams axis. In order to validate the proposed formulation, the numerical results are compared to a well know finite element code.     

Time dependent response of low velocity impact induced composite conical shells under multiple delamination

This paper presents the time dependent response of multiple delaminated angle-ply composite pretwisted conical shells subjected to low velocity normal impact. The finite element formulation is based on Mindlin’s theory incorporating rotary inertia and effects of transverse shear deformation. An eight-noded isoparametric plate bending element is employed to satisfy the compatibility of deformation and equilibrium of resultant forces and moments at the delamination crack front. A multipoint constraint algorithm is incorporated which leads to asymmetric stiffness matrices. The modified Hertzian contact law which accounts for permanent indentation is utilized to compute the contact force, and the time dependent equations are solved by Newmark’s time integration algorithm. Parametric studies are conducted with respect to triggering parameters like laminate configuration, location of delamination, angle of twist, velocity of impactor, and impactor’s displacement for centrally impacted shells.     

Large reference displacement analysis of composite plates. Part II: Computer implementation

This investigation concerns itself with the computer implementation of the dynamic formulation of thin laminated composite plates consisting of layers of orthotropic laminae that undergo large arbitrary rigid body displacements and small elastic deformations. A finite element preprocessor computer program is developed to automatically generate the invariants of the laminae, which may have arbitrary orientations. The laminae invariants are then used to obtain the invariants of the elements and the composite laminated plate. The consistent and lumped mass formulations of the invariants of motion of composite plates are compared and it is concluded that the two methods are comparable, if a fine enough finite element mesh is used. The structure of the dynamic equations of motion, based on the formulation presented in Part I of this paper, is examined. Non-linear centrifugal and Coriolis forces arising as the result of the finite rotations of the laminae are defined, and the solution schemes of the resulting non-linear differential equations of motion are discussed. Numerical examples illustrating the differences between homogeneous isotropic and laminated composite plates are presented. An RSSR (Revolute-Spherical-Spherical-Revolute) mechanism is used in the numerical examples, with the coupler modelled as a laminated plate flexible body. It is found that the inertia of the plate contributed greatly to the transverse deformation. The effects of laminae orientation is also investigated.     

An enhanced cell‐based smoothed finite element method for the analysis of Reissner–Mindlin plate bending problems involving distorted mesh

In this work, an enhanced cell‐based smoothed finite element method (FEM) is presented for the Reissner–Mindlin plate bending analysis. The smoothed curvature computed by a boundary integral along the boundaries of smoothing cells in original smoothed FEM is reformulated, and the relationship between the original approach and the present method in curvature smoothing is established. To improve the accuracy of shear strain in a distorted mesh, we span the shear strain space over the adjacent element. This is performed by employing an edge‐based smoothing technique through a simple area‐weighted smoothing procedure on MITC4 assumed shear strain field. A three‐field variational principle is utilized to develop the mixed formulation. The resultant element formulation is further reduced to a displacement‐based formulation via an assumed strain method defined by the edge‐smoothing technique. As the result, a new formulation consisting of smoothed curvature and smoothed shear strain interpolated by the standard transverse displacement/rotation fields and smoothing operators can be shown to improve the solution accuracy in cell‐based smoothed FEM for Reissner–Mindlin plate bending analysis. Several numerical examples are presented to demonstrate the accuracy of the proposed formulation.Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

基于平面弹性-板弯曲比拟理论的Wilson型平板壳单元

平面弹性—板弯曲比拟关系为充分利用现有的平面弹性单元构造新型板弯曲单元提供了有效的途径。按膜、板两部分平行列式的观点,新方法还可以利用平面弹性单元构造性能良好的新型平板壳单元。按此观点根据著名的平面弹性Wilson元(QM6)的列式方式构造出一个新的十六自由度平板壳单元。单元构造简单,无虚假自由度,数值结果表明具有很好的收敛性和精度。     

Finite element method for non-linear forced vibrations of circular plates

Geometric non-linearities for large amplitude free and forced vibrations of circular plates are investigated. In-plane displacement and in-plane inertia are included in the formulation. The finite element method is used. An harmonic force matrix for non-linear forced vibration analysis is introduced and derived. Various out-of-plane and in-plane boundary conditions are considered. The relations of amplitude and frequency ratio for different boundary conditions and various load conditions are presented.     

A quadrilateral hybrid-Trefftz plate bending element for the inclusion of warping based on a three-dimensional plate formulation

A plate formulation, for the inclusion of warping and transverse shear deformations, is considered. From a complete thick and thin plate formulation, which was derived without ad hoc assumptions from the three-dimensional equations of elasticity for isotropic materials, the bending solution, involving powers of the thickness co-ordinate z, is used for constructing a quadrilateral finite plate bending element. The constructed element trial functions, for the displacements and stresses, satisfy, a priori, the three-dimensional Navier equations and equilibrium equations, respectively. For the coupling of the elements, independently assumed functions on the boundary are used. High accuracy for both displacements and stresses (including transverse shear stresses) can be achieved with rather coarse meshes for thick and thin plates.     

An assumed strain triangular curved solid shell element formulation for analysis of plates and shells undergoing finite rotations

A formulation for 36‐DOF assumed strain triangular solid shell element is developed for efficient analysis of plates and shells undergoing finite rotations. Higher order deformation modes described by the bubble function displacements are added to the assumed displacement field. The assumed strain field is carefully selected to alleviate locking effect. The resulting element shows little effect of membrane locking as well as shear locking, hence, it allows modelling of curved shell structures with curved elements. The kinematics of the present formulation is purely vectorial with only three translational degrees of freedom per node. Accordingly, the present element is free of small angle assumptions, and thus it allows large load increments in the geometrically non‐linear analysis. Various numerical examples demonstrate the validity and effectiveness of the present formulation. Copyright © 2001 John Wiley & Sons, Ltd.     

A simple cubic displacement element for plate bending

Early attempts to construct a triangular finite element for plate bending problems from a compatible cubic displacement field are not entirely satisfactory. The present paper shows how an accurate plate element can be achieved using independent cubic polynomial assumptions for the internal and boundary displacements in conjunction with a modified potential energy principle. This approach yields a simple algebraic formulation with favourable connection quantities at the element vertices which will appeal to practical users of the conventional finite element displacement method. Moreover, in Appendix I it is shown that the cubic element is identical to a previous hybrid stress element with linear internal bending and twisting moments and cubic boundary displacements. The stresses obtained from the former hybrid finite element solution therefore satisfy the strain compatibility conditions exactly. This remarkable result has an important significance in the theory of hybrid finite elements.     

Finite element analysis of stiffened laminated plates under transverse loading

A finite element model is developed to study the behavior of stiffened laminated plates under transverse loadings. Transverse shear flexibility is incorporated in both beam and plate displacement fields. A laminated plate element with 45 degrees of freedom is used in conjunction with a laminated beam element having 12 degrees of freedom for the bending analysis of eccentrically-stiffened laminated plates. The validity of the formulation is demonstrated by comparing with the available solutions in the literature. The numerical results are presented for eccentrically-stiffened layered plates having various boundary conditions and with stiffeners varying in number.     

LOW-VELOCITY IMPACT RESPONSE OF COMPOSITE LAMINATES CONSIDERING HIGHER-ORDER SHEAR DEFORMATION AND LARGE DEFLECTION

SUMMARY

Low-velocity impact responses of composite laminates are investigated analytically and experimentally. In analytical research, the finite element analyses based on various plate theories and three-dimensional thoeory are performed. For experimental research, a drop weight type impact test system is used. In geometrical non-linear analysis, a displacement field considering higher-order shear deformation and large deflection of the laminate is assumed and the finite element formulation is derived. The modified Hertzian contact law is incorporated into the finite element program to evaluate contact force. Numerical results including impact force histories, deflections, dynamic strains in the laminate from the impact response analysis are presented and compared with the experimental results from impact test. The results of the investigation indicate that higher-order shear deformation and large deflection effects should be considered to accurately describe the low-velocity impact response including interlaminar shear stress of the laminate.