A finite element formulation for piezoelectric shell structures considering geometrical and material non‐linearities

In this paper, we present a non‐linear finite element formulation for piezoelectric shell structures. Based on a mixed multi‐field variational formulation, an electro‐mechanical coupled shell element is developed considering geometrically and materially non‐linear behavior of ferroelectric ceramics. The mixed formulation includes the independent fields of displacements, electric potential, strains, electric field, stresses, and dielectric displacements. Besides the mechanical degrees of freedom, the shell counts only one electrical degree of freedom. This is the difference in the electric potential in the thickness direction of the shell. Incorporating non‐linear kinematic assumptions, structures with large deformations and stability problems can be analyzed. According to a Reissner–Mindlin theory, the shell element accounts for constant transversal shear strains. The formulation incorporates a three‐dimensional transversal isotropic material law, thus the kinematic in the thickness direction of the shell is considered. The normal zero stress condition and the normal zero dielectric displacement condition of shells are enforced by the independent resultant stress and the resultant dielectric displacement fields. Accounting for material non‐linearities, the ferroelectric hysteresis phenomena are considered using the Preisach model. As a special aspect, the formulation includes temperature‐dependent effects and thus the change of the piezoelectric material parameters due to the temperature. This enables the element to describe temperature‐dependent hysteresis curves. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Hybrid‐stabilized solid‐shell model of laminated composite piezoelectric structures under non‐linear distribution of electric potential through thickness

Eighteen‐node solid‐shell finite element models have been developed for the analysis of laminated composite plate/shell structures embedded with piezoelectric actuators and sensors. The explicit hybrid stabilization method is employed to formulate stabilization vectors for the uniformly reduced integrated 18‐node three‐dimensional composite solid element. Unlike conventional piezoelectric elements, the concept of the electric nodes introduced in this paper can effectively eliminate the burden of constraining the equality of the electric potential for the nodes lying on the same electrode. Furthermore, the non‐linear distribution of electric potential in the piezoelectric layer is expressed by introducing internal electric potential, which not only can simplify modelling but also obtains the same as the exact solution. Copyright © 2003 John Wiley & Sons, Ltd.     

A robust non‐linear mixed hybrid quadrilateral shell element

In the paper a non‐linear quadrilateral shell element for the analysis of thin structures is presented. The variational formulation is based on a Hu–Washizu functional with independent displacement, stress and strain fields. The interpolation matrices for the mid‐surface displacements and rotations as well as for the stress resultants and strains are specified. Restrictions on the interpolation functions concerning fulfillment of the patch test and stability are derived. The developed mixed hybrid shell element possesses the correct rank and fulfills the in‐plane and bending patch test. Using Newton's method the finite element approximation of the stationary condition is iteratively solved. Our formulation can accommodate arbitrary non‐linear material models for finite deformations. In the examples we present results for isotropic plasticity at finite rotations and small strains as well as bifurcation problems and post‐buckling response. The essential feature of the new element is the robustness in the equilibrium iterations. It allows very large load steps in comparison to other element formulations. Copyright © 2005 John Wiley & Sons, Ltd.     

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 new finite‐element formulation for electromechanical boundary value problems

In this paper, a new finite‐element formulation for the solution of electromechanical boundary value problems is presented. As opposed to the standard formulation that uses scalar electric potential as nodal variables, this new formulation implements a vector potential from which components of electric displacement are derived. For linear piezoelectric materials with positive definite material moduli, the resulting finite‐element stiffness matrix from the vector potential formulation is also positive definite. If the material is non‐linear in a fashion characteristic of ferroelectric materials, it is demonstrated that a straightforward iterative solution procedure is unstable for the standard scalar potential formulation, but stable for the new vector potential formulation. Finally, the method is used to compute fields around a crack tip in an idealized non‐linear ferroelectric material, and results are compared to an analytical solution. Copyright © 2002 John Wiley & Sons, Ltd.     

A mixed shell formulation accounting for thickness strains and finite strain 3d material models

A non‐linear quadrilateral shell element for the analysis of thin structures is presented. The Reissner–Mindlin theory with inextensible director vector is used to develop a three‐field variational formulation with independent displacements, stress resultants and shell strains. The interpolation of the independent shell strains consists of two parts. The first part corresponds to the interpolation of the stress resultants. Within the second part independent thickness strains are considered. This allows incorporation of arbitrary non‐linear 3d constitutive equations without further modifications. The developed mixed hybrid shell element possesses the correct rank and fulfills the in‐plane and bending patch test. The essential feature of the new element is the robustness in the equilibrium iterations. It allows very large load steps in comparison with other element formulations. We present results for finite strain elasticity, inelasticity, bifurcation and post‐buckling problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Non‐linear exact geometry 12‐node solid‐shell element with three translational degrees of freedom per node

This paper presents the finite rotation exact geometry (EG) 12‐node solid‐shell element with 36 displacement degrees of freedom. The term ‘EG’ reflects the fact that coefficients of the first and second fundamental forms of the reference surface and Christoffel symbols are taken exactly at each element node. The finite element formulation developed is based on the 9‐parameter shell model by employing a new concept of sampling surfaces (S‐surfaces) inside the shell body. We introduce three S‐surfaces, namely, bottom, middle and top, and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows one to represent the finite rotation higher order EG solid‐shell element formulation in a very compact form and to derive the strain–displacement relationships, which are objective, that is, invariant under arbitrarily large rigid‐body shell motions in convected curvilinear coordinates. The tangent stiffness matrix is evaluated by using 3D analytical integration and the explicit presentation of this matrix is given. The latter is unusual for the non‐linear EG shell element formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Four‐node semi‐EAS element in six‐field nonlinear theory of shells

We propose a new four‐node C⁰ finite element for shell structures undergoing unlimited translations and rotations. The considerations concern the general six‐field theory of shells with asymmetric strain measures in geometrically nonlinear static problems. The shell kinematics is of the two‐dimensional Cosserat continuum type and is described by two independent fields: the vector field for translations and the proper orthogonal tensor field for rotations. All three rotational parameters are treated here as independent. Hence, as a consequence of the shell theory, the proposed element has naturally six engineering degrees of freedom at each node, with the so‐called drilling rotation. This property makes the element suitable for analysis of shell structures containing folds, branches or intersections. To avoid locking phenomena we use the enhanced assumed strain (EAS) concept. We derive and linearize the modified Hu–Washizu principle for six‐field theory of shells. What makes the present approach original is the combination of EAS method with asymmetric membrane strain measures. Based on literature, we propose new enhancing field and specify the transformation matrix that accounts for the lack of symmetry. To gain knowledge about the suitability of this field for asymmetric strain measures and to assess the performance of the element, we solve typical benchmark examples with smooth geometry and examples involving orthogonal intersections of shell branches. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical integration of the stiff dynamics of geometrically exact shells: an energy‐dissipative momentum‐conserving scheme

This paper presents a new family of time‐stepping algorithms for the integration of the dynamics of non‐linear shells. We consider the geometrically exact shell theory involving an inextensible director field (the so‐called five‐parameter shell model). The main characteristic of this model is the presence of the group of finite rotations in the configuration manifold describing the deformation of the solid. In this context, we develop time‐stepping algorithms whose discrete solutions exhibit the same conservation laws of linear and angular momenta as the underlying physical system, and allow the introduction of a controllable non‐negative energy dissipation to handle the high numerical stiffness characteristic of these problems. A series of algorithmic parameters for the different components of the deformation of the shell (i.e. membrane, bending and transverse shear) fully control this numerical dissipation, recovering existing energy‐momentum schemes as a particular choice of these algorithmic parameters. We present rigorous proofs of the numerical properties of the resulting algorithms in the full non‐linear range. Furthermore, it is argued that the numerical dissipation is introduced in the high‐frequency range by considering the proposed algorithm in the context of a linear problem. The finite element implementation of the resulting methods is described in detail as well as considered in the final arguments proving the aforementioned conservation/dissipation properties. We present several representative numerical simulations illustrating the performance of the newly proposed methods. The robustness gained over existing methods in these stiff problems is confirmed in particular. Copyright © 2002 John Wiley & Sons, Ltd.     

Efficient and accurate multilayer solid‐shell element: non‐linear materials at finite strain

We present in this paper an efficient and accurate low‐order solid‐shell element formulation for analyses of large deformable multilayer shell structures with non‐linear materials. The element has only displacement degrees of freedom (dofs), and an optimal number of enhancing assumed strain (EAS) parameters to pass the patch tests (both membrane and out‐of‐plane bending) and to remedy volumetric locking. Based on the mixed Fraeijs de Veubeke‐Hu‐Washizu (FHW) variational principle, the in‐plane and out‐of‐plane bending behaviours are improved and the locking associated with (nearly) incompressible materials is avoided via a new efficient enhancement of strain tensor. Shear locking and curvature thickness locking are resolved effectively by using the assumed natural strain (ANS) method. Two non‐linear 3‐D constitutive models (Mooney–Rivlin material and hyperelastoplastic material at finite strain) are applied directly without requiring the enforcement of the plane‐stress assumption. In particular, we give a simple derivation for the hyperelastoplastic model using spectral representations. In addition, the present element has a well‐defined lumped mass matrix, and provides double‐side contact surfaces for shell contact problems. With the dynamics referred to a fixed inertial frame, the present element can be used to analyse multilayer shell structures undergoing large overall motion. Numerical examples involving static analyses and implicit/explicit dynamic analyses of multilayer shell structures with both material and geometric non‐linearities are presented, and compared with existing results obtained from other shell elements and from a meshless method. It is shown that elements that did not pass the out‐of‐plane bending patch test could not provide accurate results, as compared to the present element formulation, which passed the out‐of‐plane bending patch test. The present element proves to be versatile and efficient in the modelling and analyses of general non‐linear composite multilayer shell structures. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical analysis of an elasto‐piezoelectric problem with damage

In this paper, we analyze an algorithm for the quasistatic evolution of the mechanical state of an elasto‐piezoelectric body with damage. Both damage, caused by the development and the growth of internal microcracks, and piezoelectric effects, are included in the model. The mechanical problem is expressed as an elliptic system for the displacement field coupled with a non‐linear parabolic partial differential equation for the damage field and a linear partial differential equation for the electric potential. The variational formulation leads to a coupled system composed of two linear variational equations for the displacement field and the electric potential, and a non‐linear parabolic variational equation for the damage field. The existence of a unique weak solution is stated. Then, a fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. Error estimates are derived on the approximate solutions, from which the linear convergence of the algorithm is deduced under suitable regularity conditions. Finally, some numerical simulations are performed, in one, two and three dimensions, to demonstrate the accuracy of the scheme and the behaviour of the solution. Copyright © 2008 John Wiley & Sons, Ltd.     

Optimal solid shell element for large deformable composite structures with piezoelectric layers and active vibration control

In this paper, we present an optimal low‐order accurate piezoelectric solid‐shell element formulation to model active composite shell structures that can undergo large deformation and large overall motion. This element has only displacement and electric degrees of freedom (dofs), with no rotational dofs, and an optimal number of enhancing assumed strain (EAS) parameters to pass the patch tests (both membrane and out‐of‐plane bending). The combination of the present optimal piezoelectric solid‐shell element and the optimal solid‐shell element previously developed allows for efficient and accurate analyses of large deformable composite multilayer shell structures with piezoelectric layers. To make the 3‐D analysis of active composite shells containing discrete piezoelectric sensors and actuators even more efficient, the composite solid‐shell element is further developed here. Based on the mixed Fraeijs de Veubeke–Hu–Washizu (FHW) variational principle, the in‐plane and out‐of‐plane bending behaviours are improved via a new and efficient enhancement of the strain tensor. Shear‐locking and curvature thickness locking are resolved effectively by using the assumed natural strain (ANS) method. We also present an optimal‐control design for vibration suppression of a large deformable structure based on the general finite element approach. The linear‐quadratic regulator control scheme with output feedback is used as a control law on the basis of the state space model of the system. Numerical examples involving static analyses and dynamic analyses of active shell structures having a large range of element aspect ratios are presented. Active vibration control of a composite multilayer shell with distributed piezoelectric sensors and actuators is performed to test the present element and the control design procedure. Copyright © 2005 John Wiley & Sons, Ltd.     

Two refined non‐conforming quadrilateral flat shell elements

Two refined quadrilateral flat shell elements named RSQ20 and RSQ24 are constructed in this paper based on the refined non‐conforming element method, and the elements can satisfy the displacement compatibility requirement at the interelement of the non‐planar elements by introducing the common displacements suggested by Chen and Cheung. A refined quadrilateral plate element RPQ4 and a plane quadrilateral isoparametric element are combined to obtain the refined quadrilateral flat shell element RSQ20, and a refined quadrilateral flat shell element RSQ24 is constructed on the basis of a RPQ4 element and a quadrilateral isoparametric element with drilling degrees of freedom. The numerical examples show that the present method can improve the accuracy of shell analysis and that the two new refined quadrilateral flat shell elements are efficient and accurate in the linear analysis of some shell structures. Copyright © 2000 John Wiley & Sons, Ltd.     

Geometrically non-linear analysis of composite structures with integrated piezoelectric sensors and actuators

This paper deals with the geometrically non-linear analysis of thin plate/shell laminated structures with embedded integrated piezoelectric actuators or sensors layers and/or patches. The motivation for the present developments is the lack of studies in the behavior of adaptive structures using geometrically non-linear models, where only very few published works were found in the open literature.

The model is based on the Kirchhoff classical laminated theory and can be applied to plate and shell adaptive structures with arbitrary shape, general mechanical and electrical loadings.

The finite element model is a non-conforming single layer triangular plate/shell element with 18 degrees of freedom for the generalized displacements and one electrical potential degree of freedom for each piezoelectric layer or patch.

An updated Lagrangian formulation associated to Newton–Raphson technique is used to solve incrementally and iteratively the equilibrium equations.

The model is applied in the solution of four illustrative cases, and the results are compared and discussed with alternative solutions when available.     

Energy and momentum conserving elastodynamics of a non‐linear brick‐type mixed finite shell element

The paper presents aspects of the finite element formulation of momentum and energy conserving algorithms for the non‐linear dynamic analysis of shell‐like structures. The key contribution is a detailed analysis of the implementation of a Simó–Tarnow‐type conservation scheme in a recently developed new mixed finite shell element. This continuum‐based shell element provides a well‐defined interface to strain‐driven constitutive stress updates algorithms. It is based on the classic brick‐type trilinear displacement element and is equipped with specific gradient‐type enhanced strain modes and shell‐typical assumed strain modifications. The excellent performance of the proposed dynamic shell formulation with respect to conservation properties and numerical stability behaviour is demonstrated by means of three representative numerical examples of elastodynamics which exhibit complex free motions of flexible structures undergoing large strains and large rigid‐body motions. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite element model for thermomechanical analysis of sheet metal forming

A thermal model based on explicit time integration is developed and implemented into the explicit finite element code DYNA3D to model simultaneous forming and quenching of thin‐walled structures. A staggered approach is used for coupling the thermal and mechanical analysis, wherein each analysis is performed with different time step sizes. The implementation includes a thermal shell element with linear temperature approximation in the plane and quadratic in the thickness direction, and contact heat transfer. The material behaviour is described by a temperature‐dependent elastic–plastic model with a non‐linear isotropic hardening law. Transformation plasticity is included in the model. Examples are presented to validate and evaluate the proposed model. The model is evaluated by comparison with a one‐sided forming and quenching experiment. Copyright © 2004 John Wiley & Sons, Ltd.     

A simplified co‐rotational method for quadrilateral shell elements in geometrically nonlinear analysis

This paper presents a simplified co‐rotational formulation for quadrilateral shell elements inheriting the merit of element‐independence from the traditional co‐rotational approach in literature. With the objective of application to nonlinear analysis of civil engineering structures, the authors further simplify the formulation of the geometrical stiffness using the small strain assumption, which is valid in the co‐rotational approach, with the warping effects considered as eccentricities. Compared with the traditional element‐independent co‐rotational method, the projector is neglected both in the tangent stiffness matrix and in the internal force vector for simplicity in formulation. Meanwhile, a quadrilateral flat shell element allowing for drilling rotations is adopted and incorporated into this simplified co‐rotational algorithm for geometrically nonlinear analysis involved with large displacements and large rotations. Several benchmark problems are presented to confirm the efficiency and accuracy of the proposed method for practical applications.