Finite element analysis of laminated composite plates using a higher-order displacement model

A C° continuous displacement finite element formulation of a higher-order theory for flexure of thick arbitrary laminated composite plates under transverse loads is presented. The displacement model accounts for non-linear and constant variation of in-plane and transverse displacement model eliminates the use of shear correction coefficients. The discrete element chosen is a nine-noded quadrilateral with nine degrees-of-freedom per node. Results for plate deformations, internal stress-resultants and stresses for selected examples are shown to compare well with the closed-form, the theory of elasticity and the finite element solutions with another higher-order displacement model by the same authors. A computer program has been developed which incorporates the realistic prediction of interlaminar stresses from equilibrium equations.     

Interlaminar stresses in laminated composite plates, cylindrical/spherical shell panels damaged by low-velocity impact

Prediction of damage caused by low-velocity impact in laminated composite plate cylindrical/spherical shell panels is an important problem faced by designers using composites. Not only the in-plane stresses but also the interlaminar normal and shear stresses play a role in estimating the damage caused. The work reported here is an effort in getting better predictions of damage in composite plate cylindrical/spherical shell panels subjected to low-velocity impact.

The low-velocity impact problem is treated as a quasi-static problem. First, the in-plane stresses are calculated by 2-D nonlinear finite element analysis using a 48 degrees of freedom laminated composite shell element. The damage analysis is then carried out using a Tsai-Wu quadratic failure criterion and a maximum stress criteria. Interlaminar normal and shear stresses are predicted after taking into account the in-plane damage caused by low-velocity impact. The interlaminar stresses are obtained by integrating the 3-D equations of equilibrium through the thickness. The deformed geometry is taken into account in the third equation of equilibrium (in the thickness direction). After evaluating the formulation and the computer program developed for correctness, the interlaminar stresses are predicted for composite plates/shell panels which are damaged by low-velocity impact.     

A simple finite element formulation of a higher-order theory for unsymmetrically laminated composite plates

A higher-order theory which satisfies zero transverse shear stress conditions on the bounding planes of a generally laminated fibre-reinforced composite plate subjected to transverse loads is developed. The displacement model accounts for non-linear distribution of inplane displacement components through the plate thickness and the theory requires no shear correction coefficients. A C∘ continuous displacement finite element formulation is presented and the coupled membrane-flexure behaviour of laminated plates is investigated. The nodal unknowns are the three displacements, two rotations and two higher-order functions as the generalized degrees of freedom. The simple isoparametric formulation developed here is capable of evaluating transverse shears and transverse normal stress accurately by using the equilibrium equations. The accuracy of the nine-noded Lagrangian quadrilateral element is then established by comparing the present results with the closed-form, three-dimensional elasticity and other finite element available solutions.     

A global-local higher order theory for multilayered shells and the analysis of laminated cylindrical shell panels

Based on the global-local superposition technique proposed by Li and Liu [Li XY, Liu D. Generalized laminate theories based on double superposition hypothesis. Int J Numer Meth Eng 1997;40:1197–212.], a global-local higher order laminated shell model is proposed for predicting both displacement and stress distributions through the thickness of laminated shells. This shell model satisfies transverse shear stress continuity conditions at interfaces as well as free surface conditions of transverse shear stresses. The merit of this model is that transverse shear stresses can be accurately predicted directly from constitutive equations without smoothing techniques. Cylindrical bending of laminated and sandwich shell panels is chosen to assess the present model wherein the results from several 2D laminated shell models and three-dimensional elasticity solution are available for comparison. In addition, thermal bending and thermal expansion of laminated cylindrical shell panels are also considered in this paper.     

层板弯曲与振动的近似分析

横向剪切变形对复合材料层合板弯曲与振动的影响甚大。在本文的近似分析中,假定板在弯曲时横向位移沿整个板厚为常量。横向剪切应变沿各层厚度方向也分别为常量,但各层不同。文中以特殊正交各向异性层合板为例,采用两种不同的方法建立了各层剪切应变间的关系,推演了层合板横向弯曲与振动的微分方程组及边界条件。算例表明,即使层合板的跨——厚比很小,用本文两种分析方案计算位移、应力及固有频率,都仍具有较高的精度。      

Reddy-type zig-zag model for multilayered composite plate based on the HW variational theorem

Reddy's higher-order theory is quite attractive, but it could not describe a zig-zag shape distribution of in-plane displacement through the thickness direction and violates the continuity of transverse shear stresses at interfaces. This is due to neglect of the zig-zag function in the in-plane displacement field. Thus, a Reddy-type higher-order zig-zag theory is developed for analysis of multilayered composite plates. The developed model differs from existing ones by two features. First, a Reddy-type zig-zag function (RZZF) satisfying the bounding surface free traction condition is constructed. By introducing the RZZF into Reddy's model, a Reddy-type higher-order zig-zag model can be obtained. Second, a functional suitable for composite plate has been presented to obtain improved transverse shear stresses by employing the three-field Hu–Washizu (HW) variational principle. It is significant that the higher-order derivatives of displacement parameters in expression of transverse shear stresses have been eliminated, which is convenient for the model's finite element implementation. Equilibrium equations and analytical solution can be also presented by means of the HW variational principle. The performance of the proposed model is tested with different numerical examples, and numerical results show its accuracy and range of applicability.     

Thermal buckling of cross-ply laminated composite and sandwich plates according to a global higher-order deformation theory

A two-dimensional global higher-order deformation theory is presented for thermal buckling of cross-ply laminated composite and sandwich plates. By using the method of power series expansion of continuous displacement components, a set of fundamental governing equations which can take into account the effects of both transverse shear and normal stresses is derived through the principle of virtual work. Several sets of truncated Mth-order approximate theories are applied to solve the eigenvalue problems of a simply supported multilayered plate. Modal transverse shear and normal stresses can be calculated by integrating the three-dimensional equations of equilibrium in the thickness direction, and satisfying the continuity conditions at the interface between layers and stress boundary conditions at the external surfaces. Numerical results are compared with those of the published three-dimensional layerwise theory in which both in-plane and normal displacements are assumed to be C⁰ continuous in the continuity conditions at the interface between layers. Effects of the difference of displacement continuity conditions between the three-dimensional layerwise theory and the global higher-order theory are clarified in thermal buckling problems of multilayered composite plates.     

A quasi-conforming triangular laminated composite shell element based on a refined first-order theory

A quasi-conforming triangular laminated shell element based on a refined first-order shear deformation theory is presented. The Hu-Washizu variational principle, involving strain and displacement fields as variables, with stresses being considered as Lagrange multipliers, is used to develop the laminate composite shell element. Both strains and displacements are discretized in the element, while displacements alone are discretized at the boundary. The inter-element C ¹ continuity is satisfied a posteriori in a weak form. Due to the importance of rotations and shear deformation in the geometrically non-linear analyses of shells, 7 degrees of freedom per node are chosen, viz. three displacements, two first-derivatives in the in-plane directions of the out-of-plane displacement, and two transverse shear strains at each node. To consider the effect of transverse shear deformation on the global behavior of the laminated composite shell, the Reissner-Mindlin first-order theory, with shear correction factors of Chow and Whitney, is adopted. The transverse shear stresses are obtained through the integration of the 3-D equilibrium equations; and the warping induced by transverse shear is considered in the calculation of the in-plane stresses to improve their accuracy. Numerical examples show that the element has good convergence properties and leads to highly accurate stresses.     

A postprocess method for laminated shells with a doubly curved nine-noded finite element

A numerical method that predicts through-the-thickness stresses accurately by using in-plane displacement of Efficient Higher Order Shell Theory (EHOST) as a postprocessor is implemented in nine-noded doubly curved shell element. In the present study, an efficient postprocess method is developed in the framework of shell finite element without losing the accuracy of solutions. This method consists of two steps. First is to obtain the relationship between shear angles of First Order Shear Deformation Theory (FSDT) and EHOST. Second is to construct accurate displacement and stress fields from the FSDT solution by using EHOST displacement fields as a postprocessor. To obtain accurate transverse shear stresses, integration of equilibrium equation approach is used. In the course of calculating transverse shear stresses, the computation of third derivatives of transverse deflection is required. Simply supported curved panels and finite cylinder problems demonstrate economical and accurate solution of laminate composite shells provided by the present method. The present postprocess method should work as an efficient tool in the stress analysis of multilayered thick shells.     

An assessment of five modeling approaches for thermo-mechanical stress analysis of laminated composite panels

 A study is made of the effects of variation in the lamination and geometric parameters, and boundary conditions of multi-layered composite panels on the accuracy of the detailed response characteristics obtained by five different modeling approaches. The modeling approaches considered include four two-dimensional models, each with five parameters to characterize the deformation in the thickness direction, and a predictor-corrector approach with twelve displacement parameters. The two-dimensional models are first-order shear deformation theory, third-order theory; a theory based on trigonometric variation of the transverse shear stresses through the thickness, and a discrete layer theory. The combination of the following four key elements distinguishes the present study from previous studies reported in the literature: (1) the standard of comparison is taken to be the solutions obtained by using three-dimensional continuum models for each of the individual layers; (2) both mechanical and thermal loadings are considered; (3) boundary conditions other than simply supported edges are considered; and (4) quantities compared include detailed through-the-thickness distributions of transverse shear and transverse normal stresses. Based on the numerical studies conducted, the predictor-corrector approach appears to be the most effective technique for obtaining accurate transverse stresses, and for thermal loading, none of the two-dimensional models is adequate for calculating transverse normal stresses, even when used in conjunction with three-dimensional equilibrium equations. Received 14 September 1999     

Analysis of thick orthotropic laminated composite plates based on higher order shear deformation theory

A two-dimensional finite element model is presented to perform the linear static analysis of laminated orthotropic composite plates based on a refined higher order shear deformation theory. The theory accounts for parabolic distributions of transverse shear stresses and requires no shear correction factors. A finite element program is developed using serendipity element with seven degrees of freedom per node. The present solutions are compared with those obtained using three-dimensional elasticity theory and those obtained by other researchers. The theory accurately predicts displacements and transverse shear stresses compared to previously developed theories for thick plates and are very close to three-dimensional elasticity solutions. The effects of transverse shear deformation, material anisotropy, aspect ratio, fiber orientation and lamination sequence on transverse shear stresses are investigated. The error in values of transverse shear stresses decreases as the number of lamina increases, for a plate of same thickness. An increase in degree of anisotropy results in lower values of deflection in the plate. For cross-ply plate an increase in anisotropy results in an increase in effective stress whereas for angle-ply plate the effect is almost negligible. Through thickness variation of transverse shear stresses are independent of anisotropy. The maximum effective stress increases exponentially at lower values of anisotropy and reaches to an asymptotic value at higher values. The stacking sequence has a significant effect on the transverse deflections and shear stress. Rectangular plates experience less effective, in-plane and transverse shear stresses compared to square plates.     

Sandwich shell finite element for dynamic explicit analysis

The contribution of this paper consists of new development of transverse shear stresses through the thickness and finding an expression for the critical time step for explicit time integration of layered shells. This work presents the finite element (FE) formulation and implementation of a higher‐order shear deformable shell element for dynamic explicit analysis of composite and sandwich shells. The formulation is developed using a displacement‐based third‐order shear deformation shell theory. Using the differential equilibrium equations and the interlayer requirements, special treatment is developed for the transverse shear, resulting in a continuous, piecewise quartic distribution of the transverse shear stresses through the shell thickness. Expressions are developed for the critical time step of the explicit time integration for orthotropic homogeneous and layered shells based on the developed third‐order formulation. To assess the performance of the present shell element, it is implemented in the general non‐linear explicit dynamic FE code DYNA3D. Several problems are solved and results are presented and compared to other theoretical and numerical results. Copyright © 2002 John Wiley & Sons, Ltd.     

Efficient 2D thermo-mechanical analysis of composites and sandwich laminates

The static behavior of composites and sandwich plates in thermo-mechanical environment is investigated by a two dimensional (2D) FE model. An efficient higher-order zig-zag theory (HOZT) considering actual through-thickness temperature profile and a least square error (LSE) method to accurately predict the inter-laminar stresses is implemented in this model. The in-plane displacement field is obtained by superposing a cubically varying global displacement field on a zig-zag displacement field having different slopes at each layer. This plate theory represents parabolic through thickness variation of transverse shear stresses, which satisfy the inter-laminar continuity at the layer interfaces and zero transverse shear stress conditions at the top and bottom of the plate. In the present 2D finite element (FE) model, the first derivatives of transverse displacement have been treated as independent variables to circumvent the problem of C¹ continuity associated with the above plate theory (HOZT). The accurate through-thickness distribution of temperature is obtained by using a linear zig-zag thermal lamination theory proposed by the authors by using the thermal conduction properties of different constituent layers in the thickness direction. The LSE method is applied at the postprocessing stage to accurately calculate the inter-laminar stresses by the 3D equilibrium equations of the plate problem, after in-plane stresses are calculated. The proposed combined FE model (HOZT+LSE) is implemented to analyze the static behavior of laminated composites and sandwich plates subjected to thermo-mechanical loadings. Many new results are also presented that should be useful for the future reference.     

Accurate calculation of transverse shear stresses for soft-core sandwich laminates

Accurate evaluation of transverse stresses in soft-core sandwich laminates using the existing 2D finite element (FE) models involves cumbersome post-processing techniques. In this paper, a simple and robust method is proposed for accurate evaluation of through-the-thickness distribution of transverse stresses in soft-core sandwich laminates by using a displacement-based C⁰ continuous 2D FE model derived from refined higher-order shear deformation theory (RHSDT) and a least square error (LSE) method. In this refined higher-order shear deformation theory (RHSDT), the in-plane displacement field for the face sheets and the core is obtained by superposing a global cubically varying displacement field on a zigzag linearly early varying displacement field. The transverse displacement is assumed to have a quadratic variation within the core, and it remains constant in the faces beyond the core. The proposed C⁰ FE model satisfies the condition of transverse shear stress continuity at the layer interfaces and the zero transverse shear stress condition at the top and bottom of the sandwich plate. The nodal field variables are chosen in an efficient manner to circumvent the problem of C¹ continuity requirement of the transverse displacements associated with the RHSDT. The LSE method is applied to the 3D equilibrium equations of the plate problem at the post-processing stage, after in-plane stresses are calculated by using the above FE model based on RHSDT. Thus, the proposed method is quite simple and elegant compared to the usual method of integrating the 3D equilibrium equations at the post-processing stage for the calculation of transverse stresses in a sandwich laminates. The accuracy of the proposed method is demonstrated in the numerical examples through the comparison of the present results with those obtained from different models based on HSDT and 3D elasticity solutions.     

An improved zig-zag model for the bending of laminated composite shells

A simple layerwise higher-order zig-zag model is proposed for the bending of laminated composite shells. The model provides a cubic variation of both the in-plane displacements and the transverse shear stresses within each layer. As the displacement model satisfies the zero transverse shear stress conditions at the free surfaces, there is no need for the use of shear correction factors. By imposing the continuity of the in-plane displacements and the transverse shear stresses at the interfaces, the number of variables is shown to be the same as that given by the first-order shear deformation shell theory, irrespective of the number of layers considered. For the sake of consistency, all terms of the order of the thickness coordinate-to-radius ratio have been retained in the derivation of the governing equations. Numerical results for the cylindrical bending of thick, symmetric homogeneous orthotropic and three-layer laminated shells under sinusoidal loading show that the maximum transverse deflections and in-plane stresses are in good agreement with available exact elasticity solutions for radius-to-thickness ratios greater than or equal to four.     

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 improved isoparametric quadratic element based on refined zigzag theory to compute interlaminar stresses of multilayered anisotropic plates

An improved eight-noded isoparametric quadratic plate bending element based on refined higher-order zigzag theory (RHZT) has been developed in the present study to determine the interlaminar stresses of multilayered composite laminates. The C⁰ continuous element has been formulated by considering warping function in the displacement field based on the RHZT. Shear locking phenomenon is avoided by considering substitute shear strain field. The continuity of transverse shear stresses cannot be ensured by the proposed zigzag formulation directly, and hence, the continuity conditions of transverse shear stresses have been established by using the three-dimensional (3D) stress equilibrium equations in the present study. The transverse shear stresses are computed in a simplified manner using the differential equations of stress equilibrium. A finite element code is developed by using MATLAB software package. The performance of the present finite element model is validated by comparing the results with 3D elasticity solutions. The superiority of the proposed element in view of computational efficiency, simplicity, and accuracy has been examined by comparing the present solutions with those available in published literature using other elements.     

INTERLAMINAR STRESSES IN SHELLS OF REVOLUTION

SUMMARY

A comparative study of the interlaminar stresses in shells of revolution has been made between first order shear deformation theory (FSDT), higher order shear deformation theory with thickness stretch (HSDT7), higher order shear deformation theory with higher order inplane displacement terms (HSDT9) and three-dimensional finite element (3D) models. A semi-analytical approach is used for all the models. Interlaminar stresses are evaluated using equilibrium equations in the cases of FSDT, HSDT7 and HSDT9 models as interlaminar stresses obtained from constitutive equations are not correct. For the 3D model an eight-noded quadratic quadrilateral semi-analytic solid element is used whereas for equivalent single layer (ESL) theories a three-noded isoparametric curved element is used. Crossply parabolic and hyperbolic caps subjected to uniform external pressure and a simply supported cylindrical shell subjected to an internal sinusoidal pressure are considered in the present study.     

Modeling of functionally graded sandwich shells accounting for variation in transverse displacement

In this study, a simple C₀ isoparametric finite element formulation based on higher-order shear deformation theory is presented for static analysis of functionally graded material sandwich shells (FGMSS). To characterize the membrane-flexure behavior observed in a functionally graded shell, a displacement field involving higher-order terms in in-plane and transverse fields is considered. The proposed kinematics field incorporates for transverse normal deformation, transverse shear deformation, and nonlinear variation of the in-plane displacement field through the thickness to predict the overall response of the shell in an accurate sense. To develop the efficient C₀ formulation, the derivatives of transverse displacement are treated as independent field variables (nodal unknowns). Voigt's rule of mixture is employed to ascertain the mechanical properties of each layer's constituents along the thickness direction. A wide range of numerical problems are solved assuming various parameters: side-thickness ratio, curvature-side ratio, and gradation parameter, and their interactions with regard to static analysis of FGMSS are discussed in brief. Deflection and stresses incorporating different thickness schemes of sandwich shells are presented in the form of figures. To validate the results, a functionally graded shell without sandwich arrangement is considered. Since no results are available on static analysis of FGMSS, the present 2D model based on the finite element method might be helpful in assessing the applicability of other analytical and numerical models in this area in the future.     

High-order free vibrations of doubly-curved sandwich panels with flexible core based on a refined three-layered theory

There are few reports on the free vibration of soft core doubly-curved sandwich shells. Previous studies are largely based on the equivalent single layer theories in which the natural frequencies are grossly overestimated. This study deals with the analytical free vibrations of doubly curved sandwich shell with flexible core based on a refined general-purpose sandwich panel theory. In this theory, equations of motion are formulated based on displacements and transverse stresses at the interfaces of the core. The first order shear deformation theory and assumptions of linear distribution of transverse normal stress and uniform shear stresses over the thickness of core (based on 3D-elasticity solution of weak core) are used in the present theory. In this model, the in-plane displacements take cubic polynomial distributions and the transverse displacement has a quadratic one thorough the core thickness. Hamilton’s principle is used to obtain the equations of motion. The obtained results are validated by the analytical and numerical results published in the literatures. Parametric study is also included to investigate the effects of radius of curvature, thickness and flexibility of core.