Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite element analysis of 3-ply laminated conical shell for flutter

The results of the finite element analysis of 3-ply laminated conical shells with light core for linear panel flutter are presented and certain advantages of such shells discussed.     

Free and forced vibration analysis of laminated composite plates and shells using a 9-node assumed strain shell element

The natural frequencies of isotropic and composite laminates are presented. The forced vibration analysis of laminated composite plates and shells subjected to arbitrary loading is investigated. In order to overcome membrane and shear locking phenomena, the assumed natural strain method is used. To develop a laminated shell element for free and forced vibration analysis, the equivalent constitutive equation that makes the computation of composite structures efficient was applied. The Mindlin-Reissner theory which allows the shear deformation and rotary inertia effect to be considered is adopted for development of nine-node assumed strain shell element. The present shell element offers significant advantages since it consistently uses the natural co-ordinate system. Results of the present theory show good agreement with the 3-D elasticity and analytical solutions. In addition the effect of damping is investigated on the forced vibration analysis of laminated composite plates and shells.     

Corotational non‐linear analysis of thin plates and shells using a new shell element

A new three‐node triangular shell element is developed by combining the optimal membrane element and discrete Kirchhoff triangle (DKT) plate bending element, and is modified for laminated composite plates and shells so as to include the membrane‐bending coupling effect. Using appropriate shape functions for the bending and membrane modes of the element, the ‘inconsistent’ stress stiffness matrix is formulated and the tangent stiffness matrix is determined. Non‐linear analysis of thin‐walled structures with geometric non‐linearity is conducted using the corotational method. The new element is thoroughly tested by solving few popular benchmark problems. The results of the analysis are compared with those obtained using existing membrane elements. Copyright © 2006 John Wiley & Sons, Ltd.     

复合材料层板的拟协调罚单元*

本文在文献提出的多变量拟协调元和罚函数相结合的方法基础上构造了考慮耦合和剪切效应的十五自由度三角形层合板单元。此单元适用于薄和中厚层合板的强度分析,具有收敛和应用范围广的优点。文中以各向同性和复合材料层合板多组数例证明了该单元的有效性,同时还讨论了跨厚比L/h和弹性模量比E₁/E₂与剪切效应的关系,得出了一些有意义的结论。     

A computational model for analysing interactive buckling and delamination growth in composite structures

In this paper, a unified method is presented: (i) to model delaminated stiffened laminated composite shells; (ii) for synthesising accurate multiple post-buckling solution paths under compressive loading; and (iii) for predicting delamination growth. A multi-domain modelling technique is used for modelling the delaminated stiffened shell structures. Error-free geometrically nonlinear element formulations — a 2-noded curved stiffener element (BEAM2) and a 3-noded shell element (SHELL3) — are used for the finite element analysis. An accurate and simple automated solution strategy based on Newton type iterations is used for predicting the general geometrically nonlinear and postbuckling behaviour of structures. A simple method derived from the 3-dimensionalJ-integral is used for computing the pointwise energy release rate at the delamination front in the plate/shell models. Finally, the influence of post-buckling structural behaviour and the delamination growth on each other has been demonstrated.     

Structural analysis of composite laminates using a mixed hybrid shell element

The structural analysis of thin composite structures requires robust and effective shell elements. In this paper the variational formulation is based on a Hu–Washizu functional with independent displacements, stress resultants and shell strains. For the independent shell strains an additional interpolation part is introduced. This yields an improved convergence behaviour especially for laminated shells with coupled membrane and bending stiffness. The developed mixed hybrid shell element possesses the correct rank and fulfills the in–plane and bending patch test. The formulation is tested by several nonlinear examples including bifurcation and post–buckling response. The essential feature of the new element is the robustness in nonlinear computations with large rigid body motions. It allows very large load steps in comparison to standard displacement models.     

Impact analysis of laminated composite plates and shells by super finite elements

A super finite element method that exhibits coarse-mesh accuracy is used to predict the transient response of laminated composite plates and cylindrical shells subjected to non-penetrating impact by projectiles. The governing equations are based on the classical theories of thin laminated plates and shells taking into account the von Karman kinematics assumptions for moderately large deflections. A non-linear Hertzian-type contact law accounting for curvatures of the colliding bodies is adopted to calculate the impact force . The theoretical basis of the present finite element model is verified by analysing impact-loaded laminated composite plate and shell structures that have previously been studied through analytical or other numerical procedures. The predictive capability of the present numerical approach is successfully demonstrated through comparisons between experimentally-measured and computed force-time histories for impact of carbon fibre-reinforced plastic (CFRP) plates. The current computational model offers a relatively simple and efficient means of predicting the structural impact response of laminated composite plates and shells.     

A co-rotational 8-node assumed strain shell element for postbuckling analysis of laminated composite plates and shells

 The formulation of a nonlinear composite shell element is presented for the solution of stability problems of composite plates and shells. The formulation of the geometrical stiffness presented here is exactly defined on the midsurface and is efficient for analyzing stability problems of thin and thick laminated plates and shells by incorporating bending moment and transverse shear resultant forces. The composite element is free of both membrane and shear locking behaviour by using the assumed natural strain method such that the element performs very well as thin shells. The transverse shear stiffness is defined by an equilibrium approach instead of using the shear correction factor. The proposed formulation is computationally efficient and the test results showed good agreement. In addition the effect of the viscoelastic material is investigated on the postbuckling behaviour of laminated composite shells. Received: 6 February 2002 / Accecpted: 6 January 2003 ID=" Present address: School of Civil Engineering, Asian Institute of Technology     

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.     

Nonlinear free vibration analysis of laminated composite shells in hygrothermal environments

The nonlinear free vibration behaviour of laminated composite shells subjected to hygrothermal environments is investigated using the finite element method. The present finite element formulation considers doubly curved shells, and the Green–Lagrange type nonlinear strains are incorporated into the first-order shear deformation theory. The analysis is carried out using quadratic eight-noded isoparametric elements. The validity of the model is demonstrated by comparing the present results with the solutions available in the literature. A parametric study is carried out varying the curvature ratios and side to thickness ratios of composite cylindrical shell, spherical shell and hyperbolic paraboloid shell panels with simply supported boundary conditions.     

Static,free vibration and thermal analysis of composite plates and shells using a flat triangular shell element

Finite element static, free vibration and thermal analysis of thin laminated plates and shells using a three noded triangular flat shell element is presented. The flat shell element is a combination of the Discrete Kirchhoff Theory (DKT) plate bending element and a membrane element derived from the Linear Strain Triangular (LST) element with a total of 18 degrees of freedom (3 translations and 3 rotations per node). Explicit formulations are used for the membrane, bending and membrane-bending coupling stiffness matrices and the thermal load vector. Due to a strong analogy between the induced strain caused by the thermal field and the strain induced in a structure due to an electric field the present formulation is readily applicable for the analysis of structures excited by surface bonded or embedded piezoelectric actuators. The results are presented for (i) static analysis of (a) simply supported square plates under doubly sinusoidal load and uniformly distributed load (b) simply supported spherical shells under a uniformly distributed load, (ii) free vibration analysis of (a) square cantilever plates, (b) skew cantilever plates and (c) simply supported spherical shells; (iii) Thermal deformation analysis of (a) simply supported square plates, (b) simply supported-clamped square plate and (c) simply supported spherical shells. A numerical example is also presented demonstrating the application of the present formulation to analyse a symmetrically laminated graphite/epoxy laminate excited by a layer of piezoelectric polyvinylidene flouride (PVDF). The results presented are in good agreement with those available in the literature.The work was partly sponsored by a grant (DAAHO4-95-1-0175) from the army research office with Dr. Gary Anderson as the grant monitor.     

Transient analysis of FGM and laminated composite structures using a refined 8-node ANS shell element

In this paper, we investigate the vibration analysis of functionally graded material (FGM) and laminated composite structures, using a refined 8-node shell element that allows for the effects of transverse shear deformation and rotary inertia. The properties of FGM vary continuously through the thickness direction according to the volume fraction of constituents defined by sigmoid function, but in this method, their Poisson’s ratios of the FGM plates and shells are assumed to be constant. The finite element, based on a first-order shear deformation theory, is further improved by the combined use of assumed natural strains and different sets of collocation points for interpolation the different strain components. We analyze the influence of the shell element with the various location and number of enhanced membrane and shear interpolation. Using the assumed natural strain method with proper interpolation functions the present shell element generates neither membrane nor shear locking behavior even when full integration is used in the formulation. The natural frequencies of plates and shells are presented, and the forced vibration analysis of FGM and laminated composite plates and shells subjected to arbitrary loading is carried out. In order to overcome membrane and shear locking phenomena, the assumed natural strain method is used. To validate and compare the finite element numerical solutions, the reference solutions of plates based on the Navier’s method, the series solutions of sigmoid FGM (S-FGM) plates are obtained. Results of the present theory show good agreement with the reference solutions. In addition the effect of damping is investigated on the forced vibration analysis of FGM plates and shells.     

Finite element buckling analysis of rotationally periodic laminated composite shells

A finite element (FE) buckling analysis of rotationally periodic laminated composite shells is performed in this paper. Because the buckling mode of such structures is characterized as rotationally periodic, a corresponding FE buckling analysis scheme is proposed to reduce the computational expenses. Moreover, a new kind of relative degrees‐of‐freedom element is developed, which can be connected to other solid elements with ease and can yield satisfactory results with a relatively coarse FE mesh. Numerical results of two laminated cylindrical shells subjected to lateral pressure are compared with theoretical ones. The good agreement of them shows the validity of this new computational strategy. Finally, a practical structure is analysed to demonstrate the advantage of this method. Copyright © 2001 John Wiley & Sons, Ltd.     

Delamination modeling in doubly curved laminated shells for free vibration analysis using zigzag theory-based facet shell element and hybrid continuity method

We present a finite element (FE) formulation for the free vibration analysis of doubly curved laminated composite and sandwich shells having multiple delaminations, employing a facet shell element based on the efficient third-order zigzag theory and the region approach of modeling delaminations. The methodology, hitherto not attempted, is general for delaminations occurring at multiple interfacial and spatial locations. A recently developed hybrid method is used for satisfying the continuity of the nonlinear layerwise displacement field at the delamination fronts. The formulation is shown to yield very accurate results with reference to full-field three-dimensional FE solutions, for the natural frequencies and mode shapes of delaminated shallow and deep, composite and highly inhomogeneous soft-core sandwich shells of different geometries and boundary conditions, with a significant computational advantage. The accuracy is sensitive to the continuity method used at the delamination fronts, the usual point continuity method yielding rather poor accuracy, and the proposed hybrid method giving the best accuracy. Such efficient modeling of laminated shells with delamination damage will be of immense use for model-based techniques for structural health monitoring of laminated shell-type structures.     

旋转周期复合材料层合结构的有限元屈曲分析

针对空间结构中常见的复合材料层合壳体结构发展了一种多层相对自由度层合壳元。这种实体型壳元既可以用较粗的网格很好地模拟层合壳, 又易与三维实体单元相连接, 使变厚度、带有补强的复合材料层合壳体等复杂结构得以正确建模。同时运用旋转周期有限元技术对大规模的空间复合材料层合结构成功实施了屈曲分析。数值算例验证了本文计算策略的有效性。      

一类双壳结构振动特性的半解析法

以Hamilton正则方程的半解析法为基础,为一类双壳结构的振动特性提出了一种新的数学模型。基本步骤:(1)独立地建立内外壳和连接筋的线性方程组;(2)考虑到内外壳和连接筋的界面上的应力和位移的连续性,联立内外壳和连接筋的方程,从而得到全结构的方程组。主要优点是:采用了同一种Hamiltonian等参元离散壳和连接筋,结构的转动惯性、剪切变形等因素都得到了考虑,而且不限制壳的厚度和筋的高度;该方法象一般的有限元法一样适应复杂的边界条件和由多种材料构成的结构。本文的方法可推广用来研究加筋复合材料或加筋压电材料层合壳及相应的双壳结构的动力学问题。     

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 four‐noded quadrilateral element for composite laminated plates/shells using the refined zigzag theory

A new bilinear four‐noded quadrilateral element (called quadrilateral linear refined zigzag) for the analysis of composite laminated and sandwich plates/shells based on the refined zigzag theory is presented. The element has seven kinematic variables per node. Shear locking is avoided by introducing an assumed linear shear strain field. The performance of the element is studied in several examples where the reference solution is the 3D finite element analysis using 20‐noded hexahedral elements. Copyright © 2013 John Wiley & Sons, Ltd.