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.     

The Displacement Perspective During Ultimate Failure of Composite Laminates

This paper deals with the studies on the state of displacement of symmetric and anti-symmetric angle-ply and cross-ply laminated composite plates during its ultimate failure, subjected to transverse static load. First-order shear deformation theory (FSDT) is employed in conjunction with the finite element approach using eight-noded quadratic isoparametric element. The free vibration analyses of isotropic and laminated composite plates are carried out to ensure the overall validity of the present finite element formulation. The mid surface of the laminate is considered as the reference plane. The principal material directions in different laminae are oriented to produce a laminated structural element capable of resisting loads in several directions. The stiffness of a composite laminate is obtained from the properties of the constituent laminae. The affected stiffness of the failed lamina is discarded completely after the failure of weakest ply. The rigidity matrix of the laminate with remaining laminae is re-established. The re-evaluation process continues until the laminate fails completely. To investigate the displacement behaviour of laminates during the ultimate failure, parametric studies are carried out for different cases by varying the stacking sequences, fiber orientations, layer thicknesses, aspect ratios and the number of layers in the laminate. The comparison of results in terms of non-dimensional natural frequencies and ply-by-ply failure analyses obtained from the present investigation are made with those available in the reported literature.     

Finite elements for elasticity with microstructure and gradient elasticity

We present a general finite element discretization of Mindlin's elasticity with microstructure. A total of 12 isoparametric elements are developed and presented, six for plane strain conditions and six for the general case of three‐dimensional deformation. All elements interpolate both the displacement and microdeformation fields. The minimum order of integration is determined for each element, and they are all shown to pass the single‐element test and the patch test. Numerical results for the benchmark problem of one‐dimensional deformation show good convergence to the closed‐form solution. The behaviour of all elements is also examined at the limiting case of vanishing relative deformation, where elasticity with microstructure degenerates to gradient elasticity. An appropriate parameter selection that enforces this degeneration in an approximate manner is presented, and numerical results are shown to provide good approximation to the respective displacements and strains of a gradient elastic solid. Copyright © 2007 John Wiley & Sons, Ltd.     

A mechanistic model for transverse damage initiation,evolution, and stiffness reduction in laminated composites

A constitutive model to predict stiffness reduction due to transverse matrix cracking is derived for laminae with arbitrary orientation, subject to in-plane stress, embedded in laminates with symmetric but otherwise arbitrary laminate stacking sequence. The moduli of the damaged laminate are a function of the crack densities in the damaging laminae, which are analyzed one by one. The evolution of crack density in each lamina is derived in terms of the calculated strain energy release rate and predicted as function of the applied load using a fracture mechanics approach. Unlike plasticity-inspired formulations, the proposed model does not postulate damage evolution functions and thus there is no need for additional experimental data to adjust material parameters. All that it is needed are the elastic moduli and critical energy release rates for the laminae. The reduction of lamina stiffness is an integral part of the model, allowing for stress redistribution among laminae. Comparisons with experimental data and some results from the literature are presented.     

实体退化板单元及其在板的振动分析中的应用

经典板壳单元是由板壳理论构造出来的,而经典的板壳理论是在空间弹性理论的基础上考虑板壳的基本假定得来的。在空间等参数单元的基础上,直接引入板壳的基本假定,修改空间等参数单元的弹性矩阵,从而构造出适合于厚薄板壳分析的20结点实体退化板单元,并将其应用于开口圆柱薄壳的静力分析和厚薄板的固有振动分析。数值算例表明,该单元收敛快,稳定性好,具有较高的精度。此外,该单元还可以用于曲边变厚度板、壳体及层合板的振动分析。     

On an isoparametric finite-element for composite laminates with finite rotations

For composite laminates consisting of an arbitrary number of orthotropic laminae first a finite-rotation theory is presented as basis of isoparametric finite-element formulations. The derivation is achieved by a Reissner-Mindlin type kinematic assumption which allows a constant shear deformation distribution across the thickness. The constitutive equations are presented in a general form such that orthotropic material behaviour with material axes varying arbitrarily across the thickness may easily be considered in numerical implementation, also when using curvilinear coordinates. Special attention is taken to predict the force distribution in the deformed shell structure. This theory is then transformed into a four-node isoparametric assumed-strain finite element. Unlike in the degeneration approach, interpolation polynomials are introduced directly for rotation variables determining the deformed position of the unit normal vector. The capability of the finite element developed to deal with strongly nonlinear situations is demonstrated by many examples. Also numerical results are presented permitting a systematical comparison of classical and isoparametric approaches concerning the numerical efficiency.     

Two simple and efficient displacement‐based quadrilateral elements for the analysis of composite laminated plates

Two simple 4‐node 20‐DOF and 4‐node 24‐DOF displacement‐based quadrilateral elements named RDKQ‐L20 and RDKQ‐L24 are developed in this paper based on the first‐order shear deformation theory (FSDT) for linear analysis of thin to moderately thick laminates. The deflection and rotation functions of the element sides are obtained from Timoshenko's laminated composite beam functions. Linear displacement interpolation functions of the standard 4‐node quadrilateral isoparametric plane element and displacement functions of a quadrilateral plane element with drilling degrees of freedom are taken as in‐plane displacements of the proposed elements RDKQ‐L20 and RDKQ‐L24, respectively. Due to the application of Timoshenko's laminated composite beam functions, convergence can be ensured theoretically for very thin laminates. The elements are simple in formulation, and shear‐locking free for extremely thin laminates even with full integration. A hybrid‐enhanced procedure is employed to improve the accuracy of stress analysis, especially for transverse shear stresses. Numerical tests show that the new elements are convergent, not sensitive to mesh distortion, accurate and efficient for analysis of thin to moderately thick laminates. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

正交各向异性多层双曲壳体理论

本文提出由多个正交各向异性层粘合而成的扁平双曲壳体的一种分析方法.文中以层间位移分量作为基本未知函数,并将层间应力平衡条件在建立基本平衡方程式之前,先行引进.此外:在计算各单层横向剪切变形影响时,采用单层中面法线变形后维持直线,但放弃层中面法线及其正交线之间变形后保持直角的假定.算题实践说明应用本文方法,计算方便,精度高.      

Static and Dynamic Analysis of Laminated Thick and Thin Plates and Shells by a Very Simple Displacement-based 3-D Hexahedral Element with Over-Integration

A very simple displacement-based hexahedral 32-node element (denoted as DPH32), with over-integration in the thickness direction, is developed in this paper for static and dynamic analyses of laminated composite plates and shells. In contrast to higher-order or layer-wise higher-order plate and shell theories which are widely popularized in the current literature, the proposed method does not develop specific theories of plates and shells with postulated kinematic assumptions, but simply uses the theory of 3-D solid mechanics and the widely-available solid elements. Over-integration is used to evaluate the element stiffness matrices of laminated structures with an arbitrary number of laminae, while only one element is used in the thickness direction without increasing the number of degrees of freedom. A stress-recovery approach is used to compute the distribution of transverse stresses by considering the equations of 3D elasticity. Comprehensive numerical results are presented for static, free vibration, and transient analyses of different laminated plates and shells, which agree well with existing solutions in the published literature, or solutions of very-expensive 3D models by commercial FEM codes. It is clearly shown that the proposed methodology can accurately and efficiently predict the structural and dynamical behavior of laminated composite plates and shells in a very simple and cost-effective manner.     

复合材料层合锥柱结合壳的有限元强度分析

本文应用各向异性层合壳体理论的有限元位移法对层合锥柱结合壳进行了强度分析。给出了各向同性锥柱结合壳和复合材料层合维柱结合壳的数值结果。各向同性锥柱结合壳的数值结果和三维弹性理论的有限元法计算结果作了比较,两者符合得很好。复合材料层合锥柱结合壳的数值结果也和用轴对称曲壳单元程序的计算结果作了比较,两者也很一致。      

平面编织复合材料层板弹性性能的预测

本文在平面编织纤维增强树脂复合材料的单方向直波纹模型的基础上,分别采用经典层合板理论和有限元应变能等效方法,预测了平面编织复合材料层合板迭层对其弹性性能的影响。根据经典层合板理论,由对平面织物复合材料的单向直纹模型的上、下表面施加不同的约束条件,得出层板弹性常数变化范围。用有限元能量法,则预测出不同铺层数的编织复合材料层板的弹性性能。与实验结果比较,表明纵向模量的预测是可靠的。     

Aspects of the formulation and finite element implementation of large strain isotropic elasticity

The paper presents a formulation of isotropic large strain elasticity and addresses some computational aspects of its finite element implementation. On the theoretical side, an Eulerian setting of isotropic elasticity is discussed exclusively in terms of the Finger tensor as a strain measure. Noval aspects are a direct representation of the Eulerian elastic moduli in terms of the Finger tensor and their rigorous decomposition into decoupled volumetric and isochoric contributions based on a multiplicative split of the Finger tensor into spherical and unimodular parts. The isochoric stress response is formulated in terms of the eigenvalues of the unimodular part of the Finger tensor. A constitutive algorithm for the computation of the stresses and tangent moduli for plane problems is developed and applied to a model problem of rubber elasticity. On the computational side, the implementation of the constitutive model in three possible finite element formulations is discussed. After pointing out algorithmic techniques for the treatment of incompressible elasticity, several numerical simulations are presented which show the performance of the proposed constitutive algorithm and the convergence behaviour of the different finite element fomulations for compressible and incompressible elasticity.     

Elastic stability of skew laminated composite plates subjected to biaxial follower forces

The present study investigates the elastic stability of skew laminated composite plates subjected to biaxial inplane follower forces by the finite element method. The plate is assumed to follow first-order shear deformation plate theory (FSDPT). The kinetic and strain energies of skew laminated composite plate and the work done by the biaxial inplane follower forces are derived by using tensor theory. Then, by Hamilton's principle, the dynamic mathematical model to describe the free vibration of this problem is formed. The finite element method and the isoparametric element are utilized to discretize the continuous system and to obtain the characteristic equations of the present problem. Finally, natural vibration frequencies, buckling loads (also the instability types) and their corresponding mode shapes are found by solving the characteristic equations. Numerical results are presented to demonstrate the effects of those parameters, such as various inplane force combinations, skew angle and lamination scheme, on the elastic stability of skew laminated composite plates subjected to biaxial inplane follower forces.     

An accurate shear-deformable six-node triangular plate element for laminated composite structures

A six-node triangle plate/shell element is developed for the analysis of laminated composite structures. This model is formulated using Hamilton's principle along with a first-order (Reissner/Mindlin) shear deformation theory. The element is based upon an isoparametric representation along with an interdependent interpolation strategy; bicubic polynomials for the transverse displacement and biquadratic polynomials for the element geometry, in-plane displacements and rotations. The resulting element, which is evaluated using exact numerical integration, has correct rank and is free of shear ‘locking’. Numerical results are presented that validate the new element and prove its outstanding convergence capabilities in comparison to existing triangular elements using standardized test problems (elastic eigenvalue analysis, patch test, static simply supported square-plate solutions) and experimentally measured vibration data of cantilevered isotropic and composite plates.     

Avoiding negative elastic moduli when using Lagrange interpolation for material grading in finite element analysis

Polynomial interpolations, one of the most common interpolations used in finite element methods (FEMs), are a workhorse of many FEM codes. These interpolations are readily available for all kinds of elements, and using them for modeling the variation of elastic moduli in graded elements is thus both convenient and natural. Yet, like all polynomial interpolations, they can be prone to oscillations that can result in regions of negative elastic modulus in the element, even with only positive nodal values of elastic moduli. The result of these negative modulus regions, even if the region is small, can be unexpected singularities in the solution. This defeats the purpose of using polynomial interpolations for capturing material grading in the element. We demonstrate the issue using three-node quadratic Lagrange interpolations of material grading in otherwise isoparametric p-type elements and show how to avoid this problem.     

A discrete constitutive model for transverse and shear damage of symmetric laminates with arbitrary stacking sequence

A damage constitutive model in conjunction with a 2-D finite element discretization is presented for predicting onset and evolution of matrix cracking and subsequent stiffness reduction of symmetric composite laminates with arbitrary stacking sequence subjected to membrane loads. The formulation uses laminae crack densities as the only state variables, with crack growth driven by both mechanical stress and residual stress due to thermal expansion. The formulation is based on fracture mechanics in terms of basic materials properties, lamina moduli, and critical strain energy release rates G_IC and G_IIC, only. No additional adjustable parameters are needed to predict the damage evolution. Spurious strain localization and mesh size dependence are intrinsically absent in this formulation. Thus, there is no need to define a characteristic length. Comparison of model results to experimental data is presented for various laminate stacking sequences. Prediction of crack initiation, evolution, and stiffness degradation compare very well to experimental data.     

Numerical characterization of compressive response and damage evolution in laminated plates containing a cutout

A micromechanical constitutive model [Liang Z, Lee HK, Suaris W. Micromechanics-based constitutive modeling for unidirectional laminated composites. Int J Solid Struct 2006;43:5674–89], based on the concept of the ensemble-volume average for laminated composites, is implemented into a finite element program to numerically characterize the compressive response and damage evolution in laminated plates containing a cutout. Prior to the implementation of the model into the finite element program, the predicted moduli of laminated composites are compared with analytical bounds and experimental data for the validation and verification of the constitutive model. A series of numerical simulations for a uniaxial test of laminated plate specimens containing a cutout are conducted using the implemented constitutive model. The predictions are compared with experimental data [Lessard LB, Chang FK. Damage tolerance of laminated composites containing an open hole and subjected to compressive loadings: Part II-Experiment. J Compos Mater 1991;25:44–64; Chang FK, Lessard LB. Damage tolerance of laminated composites containing an open hole and subjected to compressive loadings: Part I-Analysis. J Compos Mater 1991;25:2–43] to verify the accuracy of the implemented constitutive model. A parametric study is also carried out to illustrate the influence of the geometry of the specimens on the behavior of laminated plates. It is shown that the implemented constitutive model is suitable for the analysis of the constitutive behavior of laminated plates having a dilute or moderate fiber volume fraction.     

FE Stress Analysis of Thick Composite Laminates with a Hole in Bending

In this work, the stress analysis of isotropic and ortotropic laminae both with a hole and without a hole and laminated both with a hole and without a hole composite plates have been examined using the finite element method. In order to solve the problem, a computer program has been written by using characteristics of eight-node isoparametric plane element. The calculations have been tested on various composite and steel materials by this program. The results have been shown in diagrams and tables and compared with literature. It was observed that the stress distributions in the plate with a hole was totally different from the plate without a hole. The analyses have showed that, the elastic stresses have been gradually reduced as moving from the first layer to the second.     

On a finite element model for the analysis of through cracks in laminated anisotropic cylindrical shells

A finite element model for the analysis of laminated composite cylindrical shells with through cracks is presented. The analysis takes into account anisotropic elastic behaviour, bending-extensional coupling and transverse shear deformation effects. The proposed finite element model is based on the approach of dividing a cracked configuration into triangular shaped singular elements around the crack tip with adjoining quadrilateral shaped regular elements. The parabolic isoparametric cylindrical shell elements (both singular and regular) used in this model employ independent displacement and rotation interpolation in the shell middle surface. The numerical comparisons show the evidence to the conclusion that the proposed model will yield accurate stress intensity factors from a relatively coarse mesh. Through the analysis of a pressurised fibre composite cylindrical shell with an axial crack, the effect of material orthotropy on the crack tip stress intensity factors is shown to be quite significant.