Exact solutions for dynamic analysis of composite plates with distributed piezoelectric layers

Exact solutions are presented for analyzing dynamics of composite plates with piezoelectric layers bonded at the top and the bottom surfaces. The expressions for mechanical displacements, stresses, electric displacements and potential are derived from constitutive relations and field equations for the piezoelectric medium under applied surface traction and electric potential. The procedure is illustrated with a simply supported symmetric cross-ply (0°/90°/0°) graphite–epoxy composite plate covered with piezoelectric material polyvinylidene fluoride (PVDF). Results are in good agreement with those obtained from finite element model.     

A class of C finite elements for the static and dynamic analysis of laminated plates

A general higher-order deformation theory is developed to analyse the behaviour of an arbitrary laminated fibre-reinforced composite plate. Three-dimensional effects such as the warping of sections and the presence of interlaminar stress field components are taken into account assuming a power series expansion of displacements along the thickness. A class of C⁰ finite element models based on this theory is then developed for mono- and bi-dimensional elements. Applications of the models to bending and vibration of laminated plates are then discussed. The present solutions are compared with those obtained using the three-dimensional elasticity theory, classical laminate theory and other higher-order theories.     

A discrete layer beam finite element for the dynamic analysis of composite sandwich beams with integral damping layers

A discrete layer finite element is presented for the dynamic analysis of laminated beams. The element uses C⁰ continuous linear and quadratic polynominals to interpolate the in-plane and transverse displacement field, respectively, and is free from the effects of shear locking. Modal frequencies and damping are estimated using both the modal strain energy method and the complex modulus method. A forced response version of the model is also presented. The model predictions are compared with experimental data for composite sandwich beams with integral damping layers. Four damping configurations are considered, a constrained layer treatment, a segmented constrained layer treatment and two internal treatments.     

Characterization of delamination effect on composite laminates using a new generalized layerwise approach

A dynamic analysis method has been developed to investigate and characterize the effect due to presence of discrete single and multiple embedded delaminations on the dynamic response of composite laminated structures with balanced/unbalanced and arbitrary stacking sequences in terms of number, placement, mode shapes and natural frequencies. A new generalized layerwise finite element model is developed to model the presence of multiple finite delamination in laminated composites. The new theory accurately predicts the interlaminar shear stresses while maintaining computational efficiency.     

Mixed least-squares finite element model for the static analysis of laminated composite plates

This paper presents a mixed finite element model for the static analysis of laminated composite plates. The formulation is based on the least-squares variational principle, which is an alternative approach to the mixed weak form finite element models. The mixed least-squares finite element model considers the first-order shear deformation theory with generalized displacements and stress resultants as independent variables. Specifically, the mixed model is developed using equal-order C⁰ Lagrange interpolation functions of high p-levels along with full integration. This mixed least-squares-based discrete model yields a symmetric and positive-definite system of algebraic equations. The predictive capability of the proposed model is demonstrated by numerical examples of the static analysis of four laminated composite plates, with different boundary conditions and various side-to-thickness ratios. Particularly, the mixed least-squares model with high-order interpolation functions is shown to be insensitive to shear-locking.     

Coordinative optimization method of composite laminated structures based on system reliability

This paper develops the coordinative optimization method based on system reliability for laminated structures. The proposed method improves the rough RBO based on first layer failure (FLF) criterion for composite laminates, and the coupling optimization method of thickness and sequence in traditional RBO strategy based on last layer failure criterion (LLF) is improved. In this paper, the finite element analysis is used to obtain the response for the failure based on two-dimension Hashin failure criterion (the limit function). Obviously, the stiffness of composite materials will decline due to destruction of elements. Therefore, stiffness degradation is considered to describe the process of damage evolution. Subsequently, combining with the branch-bound method (B&B), we can complete the search of main failure sequences and calculate the system reliability with the help of the second-order upper bound theory. In order to guarantee the efficiency and accuracy of optimization, the adaptive GA algorithm is introduced in the whole optimization procedure. After the proposed optimization policy is given in detail, two laminated structures are presented and the results are compared with the traditional optimal method based on safety factor, which demonstrates the validity and reasonability of the developed methodology.

     

基于假定应变法的协同转动四边形复合材料壳单元

为对复合材料层合板壳结构进行精确的大变形数值模拟,提出一种采用假定应变法的能分析层合结构大转动问题的协同转动四边形壳单元.该方法在建立有限元公式时引入假定应变法以克服膜闭锁和剪切闭锁的不利影响.与其他能分析大转动问题的复合材料壳单元相比,在新的协同转动框架中采用矢量型转动变量,可大大降低在非线性增量求解过程中更新转动变量的难度,且能得到对称的单元切线刚度矩阵,提高单元的计算效率.分析两个典型算例,并与其他学者的结果进行对比,结果表明在计算层合结构大转角问题时拥有较好的精度和收敛性.     

Finite element analysis of bimodulus composite stiffened thin shells of revolution

This paper is a sequel to the work published by the first and third authors[l] on stiffened laminated shells of revolution made of unimodular materials (materials having identical properties in tension and compression). A finite element analysis of laminated bimodulus composite thin shells of revolution, reinforced by laminated bimodulus composite stiffeners is reported herein. A 48 dot doubly curved quadrilateral laminated anisotropic shell of revolution finite element and it's two compatible 16 dof stiffener finite elements namely: (i) a laminated anisotropic parallel circle stiffener element (PCSE) and (ii) a laminated anisotropic meridional stiffener element (MSE) have been used iteratively.The constitutive relationship of each layer is assumed to depend on whether the fiberdirection strain is tensile or compressive. The true state of strain or stress is realized when the locations of the neutral surfaces in the shell and the stiffeners remain unaltered (to a specified accuracy) between two successive iterations. The solutions for static loading of a stiffened plate, a stiffened cylindrical shell. and a stiffened spherical shell, all made of bimodulus composite materials, have been presented.     

Three-dimensional finite-element analysis of laminated composites

A three-dimensional finite-element analysis treating the mechanical response of thick laminated composite plates in bending is presented. An isoparametric solid element with a cubic displacement expansion in planform and a linear variation through the thickness is used to model each layer of the laminate. The degrees-of-freedom of the element are retained at its boundaries so that interconnections between lamina with different fiber orientations can be made at their interfaces. An incore version of the conjugate gradient technique, which does not have bandwidth restrictions, is used to minimize the total potential energy of the system with the number of iterations to convergence being about one-fifth the total global degrees-of-freedom. Because a three-dimensional analysis is used, the effects of thickness-stretching, transverse shear, extension, and bending deformations are obtained. Comparisons with three-dimensional elasticity solutions are in excellent agreement and show the necessity of having individual elements for each layer when they have different fiber orientations and when the plates are thick.     

A three-dimensional hybrid stress isoparametric element for the analysis of laminated composite plates

In view of the increasing interest in using composite materials for aerospace structures, the analysis of laminated composite plates becomes essential. A three-dimensional eight-node hybrid stress finite element method is developed for the analysis of laminated plates. The hybrid stress model is based on the modified complementary energy principle and takes into account the transverse shear deformation effects. The displacement field is interpolated through shape functions and nodal displacements. All three displacement components are assumed to vary linearly through the thickness of each lamina. The stress field is interpolated through assumed stress polynomials with 55 stress parameters for each lamina. All six stresses are included and satisfy the homogeneous equilibrium equations. The validity of the hybrid stress finite element model is determined by comparing the predicted numerical results with the existing three-dimensional elasticity solutions. Excellent accuracy and fast convergence are observed in the numerical results.     

Sensitivity analysis and optimization of thin laminated structures with a nonsmooth eigenvalue based criterion

This paper presents a discrete model for the design sensitivity analysis of thin laminated angle-ply composite structures using a plate shell element based on a Kirchhoff discrete theory for the bending effects. To overcome the nondifferentiability of multiple eigenvalues, which may occur during a structural optimization involving free vibrations or buckling design situations, a nonsmooth eigenvalue based criterion is implemented. Angle-ply design variables and vectorial distances from the laminated midle surface to the upper surface of each layer are considered as design variables. The design sensitivities and the directional derivatives are evaluated analytically. The efficiency and accuracy of the model developed is discussed with two illustrative cases which show the need to compute sensitivities of multiple eigenvalues as directional derivatives for laminated composite structures.     

A new finite element for free edge effect analysis in laminated composites

This paper presents a finite element model based on the Multiparticle Model of Multilayered Materials (M4) developed in Institut Navier-LAMI during the last years. The laminated plate is considered as a superposition of Reissner plates coupled by interfacial stresses. An eight-node multiparticle element is developed here. This element has 5n degrees of freedom per node (n is the layer’s number of the laminate). A finite element program called MPFEAP has been developed for the implementation of the multiparticle element. The proposed finite element model is capable of computing interlaminar stresses and other localized effect which is impossible with classical 2D finite element model. Some examples of free edge problem are provided to illustrate the accuracy of the present element in predicting interlaminar stresses.     

Estimation of low-velocity impact damage in laminated composite circular plates using nonlinear finite element analysis

Nonlinear finite element analysis is used for the estimation of damage due to low-velocity impact loading of laminated composite circular plates. The impact loading is treated as an equivalent static loading by assuming the impactor to be spherical and the contact to obey Hertzian law. The stresses in the laminate are calculated using a 48 d.o.f. laminated composite sector element. Subsequently, the Tsai-Wu criterion is used to detect the zones of failure and the maximum stress criterion is used to identify the mode of failure. Then the material properties of the laminate are degraded in the failed regions. The stress analysis is performed again using the degraded properties of the plies. The iterative process is repeated until no more failure is detected in the laminate. The problem of a typical T300/N5208 composite [45 °/0 °/ − 45 °/90 °]_s circular plate being impacted by a spherical impactor is solved and the results are compared with experimental and analytical results available in the literature. The method proposed and the computer code developed can handle symmetric, as well as unsymmetric, laminates. It can be easily extended to cover the impact of composite rectangular plates, shell panels and shells.     

Damping optimization of viscoelastic laminated sandwich composite structures

Recent developments on the optimization of passive damping for vibration reduction in sandwich structures are presented in this paper, showing the importance of appropriate finite element models associated with gradient based optimizers for computationally efficient damping maximization programs. A new finite element model for anisotropic laminated plate structures with viscoelastic core and laminated anisotropic face layers has been formulated, using a mixed layerwise approach. The complex modulus approach is used for the viscoelastic material behavior, and the dynamic problem is solved in the frequency domain. Constrained optimization is conducted for the maximization of modal loss factors, using gradient based optimization associated with the developed model, and single and multiobjective optimization based on genetic algorithms using an alternative ABAQUS finite element model. The model has been applied successfully and comparative optimal design applications in sandwich structures are presented and discussed.     

Three-dimensional finite element analysis of interlaminar stresses in thick composite laminates

The three-dimensional finite element computer program has been developed to investigate interlaminar stresses in thick composite laminates. The finite element analysis is based on displacement formulation employing curved isoparametric 16-node elements. By using substructure technique, the program developed is capable of handling any composite laminates which consist of any number of orthotropic laminae and any orientations. In this paper, solid laminates and laminates with a circular hole were taken to study interlaminar stresses at the straight edge and the curved edge, respectively. Various solid laminates such as [45_n/0_n − 45_n/90_n]_s, [45/0/ − 45/90]_ns, and [45/0/ − 45/90]_sn (n = 1˜4) were analyzed. Also, [45/0/ − 45/90]_sn laminates with a circular hole were studied for n = 1 ˜ 20. The effect of laminate thickness and stacking sequence on the interlaminar stresses near the free edge was investigated. Interlaminar stresses were governed by stacking sequence rather than laminate thickness. The boundary layer width did not increase with laminate thickness but with the number of plies in the repeating unit.     

Design of laminated composite plates for optimal dynamic characteristics using a constrained global optimization technique

The lamination arrangements of moderately thick laminated composite plates for optimal dynamic characteristics are studied via a constrained multi-start global optimization technique. In the optimization process, the dynamical analysis of laminated composite plates is accomplished by utilizing a shear deformable laminated composite finite element, in which the exact expressions for determining shear correction factors were adopted and the modal damping model constructed based on an energy concept. The optimal layups of laminated composite plates with maximum fundamental frequency or modal damping are then designed by maximizing the frequency or modal damping capacity of the plate via the multi-start global optimization technique. The effects of length-to-thickness ratio, aspect ratio and number of layer groups upon the optimum fiber orientations or layer group thicknesses are investigated by means of a number of examples of the design of symmetrically laminated composite plates.     

Geometrically non-linear transient analysis of laminated composite and sandwich shells with a refined theory and C finite elements

A C⁰ continuous finite element formulation of a higher order shear deformation theory is presented for predicting the linear and geometrically non-linear, in the sense of von Karman, transient responses of composite and sandwich laminated shells. The displacement model accounts for the non-linear cubic variation of the tangential displacement components through the thickness of the shell and the theory requires no shear correction coefficients. In the time domain, the explicit central difference integrator is used in conjunction with the special mass matrix diagonalization scheme which conserves the total mass of the element and includes effects due to rotary inertia terms. Numerical results for central transverse deflection and stresses are presented for composite and sandwich laminated shells with various boundary conditions subjected to different types of loads and are compared with the results from other sources. Some new results are also included for future reference.     

A penalty model for the analysis of curved composite beams

A simple one-dimensional mechanical model for curved laminated beams is presented. The laminae composing the beam are modelled as Timoshenko beams, perfectly bonded at the interfaces. Because the laminae can rotate differently from one to the other, the cross-sections of the composite beam can warp. The elasto-static problem of the beam is formulated through the principle of stationary potential energy, imposing constraint conditions between the displacements of adjacent laminae by a penalty technique. This approach produces an approximation of radial and tangential interactions between adjacent laminae. By using four-node isoparametric finite elements, numerical values of interlaminar stresses in straight and curved laminated beams are given. They are compared with the results obtained by other authors under different conditions.