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.     

Non‐linear analysis of shells with arbitrary evolving cracks using XFEM

A new formulation and numerical procedures are developed for the analysis of arbitrary crack propagation in shells using the extended finite element method. The method is valid for completely non‐linear problems. Through‐the‐thickness cracks in sandwich shells are considered. An exact shell kinematics is presented, and a new enrichment of the rotation field is proposed which satisfies the director inextensibility condition. To avoid locking, an enhanced strain formulation is proposed for the 4‐node cracked shell element. A finite strain plane stress constitutive model based on the logarithmic corotational rate is employed. A cohesive zone model is introduced which embodies the special characteristics of the shell kinematics. Stress intensity factors are calculated for selected problems and crack propagation problems are solved. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Anisotropic elasto-plastic finite element analysis of thick and thin plates and shells

An elasto-plastic analysis of anisotropic plates and shells is undertaken by means of the finite element displacement method. A thick shell formulation accounting for shear deformation is considered, which is based on a degenerate three-dimensional continuum element. The accommodation of variable material properties, not only along the surface of the structure but also through the thickness, is made possible by a discrete layered approach. Although isoparametric elements of the Serendipity family give satisfactory solutions for thick and moderately thin shells the results exhibit ‘locking’ for an increasing ratio of span to thickness. To develop a numerical model which is applicable to thick or thin plates and shells, the nine-node Lagrangian element and the Heterosis element are also introduced into the present model. Plastic yielding is based on the Huber-Mises yield surface extended by Hill for anisotropic materials. The yield function is generalized by introducing anisotropic parameters of plasticity which are updated during the material strain hardening history. Numerical examples are presented and compared with available solutions. The effects of anisotropy on these solutions are also discussed.     

Simple and efficient shear flexible two-node arch/beam and four-node cylindrical shell/plate finite elements

Several simple and accurate C° two-node arch/beam and four-node cylindrical shell/plate finite elements are presented in this paper. The formulation used here is based on the refined theory of thick cylindrical shells and the quasi-conforming element technique. Unlike most C° elements, the element stiffness matrix presented here is given explicitly. In spite of their simplicity, these C° finite elements posseses linear bending strains and are free from the deficiencies existing in curved C° elements such as shear and membrane locking, spurious kinematic modes and numerical ill-conditioning. These finite elements are valid not only for thick/thin beams and plates, but also for arches/straight beams and cylindrical shells/plates. Furthermore, these C° elements can automatically reduce to the corresponding C¹ beam and plate elements and give the C° beam element obtained by the reduced integration as a special case. Several numerical examples indicate that the simple two-node arch/beam and four-node cylindrical shell/plate elements given in this paper are superior to the existing C° elements with the same element degrees of freedom. Only the formulation of the rectangular cylindrical shell and plate element is presented in this paper. The formulation of an arbitrarily quadrilateral plate element will be presented in a follow-up paper³².     

An assumed strain formulation of efficient solid triangular element for general shell analysis

An efficient assumed strain triangular solid element is developed for the analysis of plate and shell structures. The finite element formulation is based on the two‐field assumed strain formulation with two independent fields of assumed displacement and assumed strain. The assumed strain field is carefully selected to alleviate the shear locking effect without triggering undesirable spurious kinematic modes. The curvilinear surface of shell structures is modelled with flat facet elements to obviate the membrane locking effect. The patch tests are successfully passed, and numerical test involving various example problems demonstrates the validity and efficiency of the present element. Copyright © 2000 John Wiley & Sons, Ltd.     

A nine node finite element for analysis of geometrically non-linear shells

A nine node finite element is presented for the analysis of thin shell structures undergoing large deflection. The finite element formulation is based on the concept of degenerate solid shell element and the Hellinger-Reissner principle with independent strain. Three versions of assumed independent strain are selected to suppress spurious kinematic modes. One version leads to a finite element model which is kinematically stable at element level while the other two give globally stable models. Numerical tests indicate that the finite element model which is stable at element level may reveal the locking effect in certain cases. However, the other two models are free of locking.     

A resultant 8-node solid-shell element for geometrically nonlinear analysis

A new resultant force formulation of 8-node solid element is presented for the linear and nonlinear analysis of thin-walled structures. The global, local and natural coordinate systems were used to accurately model the shell geometry. The assumed natural strain methods with plane stress concept were implemented to remove the various locking problems appearing in thin plates and shells. The correct warping behavior in the very thin twisted beam test was obtained by using an improved Jacobian transformation matrix. The 2 × 2 Gauss integration scheme was used for the calculation of the element stiffness matrix. From the computational viewpoint, the present solid element is very efficient for a large scale of nonlinear modeling. A lot of numerical tests were carried out for the validation of the present 8-node solid-shell element and the results are in good agreement with references.An erratum to this article can be found at     

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 mixed solid‐shell element for the analysis of laminated composites

This paper presents a versatile low order locking‐free mixed solid‐shell element that can be readily employed for a wide range of linear elastic structural analyses, that is, from thick isotropic structures to multilayer anisotropic composites. This solid‐shell element has eight nodes with only displacement degrees of freedom and few assumed stress parameters that provide very accurate interlaminar stress calculations through the element thickness. These elements can be stacked on top of each other to model multilayer structures, fulfilling the interlaminar stress continuity at the interlayer surfaces and zero traction conditions on the top and bottom surfaces of the laminate. The element formulation is based on the well‐known Fraeijs de Veubeke–Hu–Washizu mixed variational principle with enhanced assumed strains formulation and assumed natural strains formulation to alleviate the different types of locking phenomena in solid‐shell elements. The distinct feature of the present formulation is its ability to accurately calculate the interlaminar stress field in multilayer structures, which is achieved by the introduction of a constraint equation on the interlaminar stresses in the Fraeijs de Veubeke–Hu–Washizu principle‐based enhanced assumed strains formulation. The intelligent computer coding of the present formulation makes the present element appropriate for a wide range of structural analyses. To assess the present formulation's accuracy, a variety of popular numerical benchmark examples related to element convergence, mesh distortion, and shell and laminated composite analyses are investigated and the results are compared with those available in the literature. These benchmark examples reveal that the proposed formulation provides very good results for the structural analysis of shells and multilayer composites. Copyright © 2011 John Wiley & Sons, Ltd.     

A nine-node assumed strain finite element for large-deformation analysis of laminated shells

An application of the element-based Lagrangian formulation is described for large-deformation analysis of both single-layered and laminated shells. Natural co-ordinate-based stresses, strains and constitutive equations are used throughout the formulation of the present shell element which offers significant implementation advantages compared with the traditional Lagrangian formulation. In order to avoid locking phenomena, an assumed strain method has been employed with judicious selection of the sampling points. Three strictly successive finite rotations are used to represent the current orientation of the shell normal. The equivalent natural constitutive equation is derived using an explicit transformation scheme to consider the multi-layer effect of laminated structures. The arc-length control method is used to trace complex load-displacement paths. Several numerical analyses are presented and discussed in order to investigate the capabilities of the present shell element. © 1998 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.     

An efficient assumed strain element model with six DOF per node for geometrically non-linear shells

The present paper describes an assumed strain finite element model with six degrees of freedom per node designed for geometrically non-linear shell analysis. An important feature of the present paper is the discussion on the spurious kinematic modes and the assumed strain field in the geometrically non-linear setting. The kinematics of deformation is described by using vector components in contrast to the conventional formulation which requires the use of trigonometric functions of rotational angles. Accordingly, converged solutions can be obtained for load or displacement increments that are much larger than possible with the conventional formulation with rotational angles. In addition, a detailed study of the spurious kinematic modes and the choice of assumed strain field reveals that the same assumed strain field can be used for both geometrically linear and non-linear cases to alleviate element locking while maintaining kinematic stability. It is strongly recommended that the element models, described in the present paper, be used instead of the conventional shell element models that employ rotational angles.     

Constitutive model and finite element formulation for large strain elasto-plastic analysis of shells

This contribution presents a refined constitutive and finite element formulation for arbitrary shell structures undergoing large elasto-plastic deformations. An elasto-plastic material model is developed by using the multiplicative decomposition of the deformation gradient and by considering isotropic as well as kinematic hardening phenomena in general form. A plastic anisotropy induced by kinematic hardening is taken into account by modifying the flow direction. The elastic part of deformations is considered by the neo-Hookean type of a material model able to deal with large strains. For an accurate prediction of complex through-thickness stress distributions a multi-layer shell kinematics is used built on the basis of a six-parametric shell theory capable to deal with large strains as well as finite rotations. To avoid membrane locking in bending dominated cases as well as volume locking caused by material incompressibility in the full plastic range the displacement based finite element formulation is improved by means of the enhanced assumed strain concept. The capability of the algorithms proposed is demonstrated by various numerical examples involving large elasto-plastic strains, finite rotations and complex through-thickness stress distributions.     

An eighteen-node solid element for thin shell analysis

An eighteen-node, three-dimensional, solid element with 54 degrees of freedom is presented for the finite element analysis of thin plates and shells. The element is based on the Hellinger-Reissner principle with independent strain. The assumed independent strain is divided into higher and lower terms. The stiffness matrix associated with the higher order independent strain plays the role of stabilization matrix. A modified stress-strain relation decoupling inplane and normal strain is used to incorporate thin shell behaviour. Numerical results demonstrate that, with a properly chosen set of assumed strain, this element is effectively free of locking even for very thin plates and shells.     

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

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

A geometrically non-linear tensorial formulation of a skewed quadrilateral thin shell finite element

A 48-degree-of-freedom (d.o.f.) skewed quadrilateral thin shell finite element, including the effect of geometrical non-linearity, is formulated and appropriate numerical procedures are adopted for the development of an efficient approach for the static and dynamic analysis of general thin shell structures. The element surface is described by a variable-order polynomial in curvilinear co-ordinates. The displacement functions are described by bicubic Hermitian polynomials in curvilinear co-ordinates. The directions of the curvilinear co-ordinates at each nodal point are uniquely defined to coincide with the directions of the boundaries of the element. In the present case of a skewed quadrilateral with non-orthogonal curvilinear coordinates, the coupling terms of the metric tensor and curvature tensor of the surface no longer vanish, such as in the case of orthogonal co-ordinates. The tensor form is used in the setup of the shape functions, geometric derivatives, stiffness matrix and computer code. This allows for the treatment of shells with irregular shapes and variable curvatures. To evaluate the efficiency and accuracy of this formulation, a systematic list of examples is chosen: (i) linear and non-linear static analysis of square and rhombic plates, cylindrical and spherical shells; (ii) linear vibrations of trapezoidal flat and curved plates; (iii) large amplitude vibrations of a rhombic plate. For the square plate and cylindrical and spherical shell, shewed element meshes with various distortion angles are used to study the effect of the distortion angles on the accuracy of the results and to demonstrate the versatility of the present element. All results are compared with alternative available solutions including those obtained using regular rectangular meshes. Pinched thin cylindrical and spherical shells are studied using different skewed meshes and various Gauss integration meshes, and no membrane locking phenomenon is observed.     

Assessment of nonlinear exact geometry sampling surfaces solid-shell elements and ANSYS solid elements for 3D stress analysis of piezoelectric shell structures

In this work, the finite rotation exact geometry four-node solid-shell element using the sampling surfaces (SaS) method is developed for the analysis of the second Piola-Kirchhoff stresses in laminated piezoelectric shells. The SaS method is based on choosing inside the layers the arbitrary number of SaS parallel to the middle surface and located at Chebyshev polynomial nodes in order to introduce the displacements and electric potentials of these surfaces as fundamental shell unknowns. The outer surfaces and interfaces are also included into a set of SaS. To circumvent shear and membrane locking, the hybrid-mixed solid-shell element on the basis of the Hu-Washizu variational principle is proposed. The tangent stiffness matrix is evaluated by 3D analytical integration throughout the finite element. This novelty provides a superior performance in the case of coarse meshes. A comparison with the SOLID226 element showed that the developed exact geometry SaS solid-shell element allows the use of load increments, which are much larger than possible with existing displacement-based finite elements. Thus, it can be recommended for the 3D stress analysis of doubly-curved laminated piezoelectric shells because the SaS formulation gives the opportunity to obtain the 3D solutions of electroelasticity with a prescribed accuracy.     

Finite Rotation Piezoelectric Exact Geometry Solid-Shell Element with Nine Degrees of Freedom per Node

This paper presents a robust non-linear piezoelectric exact geometry (EG) four-node solid-shell element based on the higher-order 9-parameter equivalent single-layer (ESL) theory, which permits one to utilize 3D constitutive equations. The term EG reflects the fact that coefficients of the first and second fundamental forms of the reference surface are taken exactly at each element node. The finite element formulation developed is based on a new concept of interpolation surfaces (I-surfaces) inside the shell body. We introduce three I-surfaces and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows us to represent the finite rotation piezoelectric higher-order EG solid-shell element formulation in a very compact form and to utilize in curvilinear reference surface coordinates the strain-displacement relationships, which are objective, that is, invariant under arbitrarily large rigid-body shell motions. To avoid shear and membrane locking and have no spurious zero energy modes, the assumed displacement-independent strain and stress resultant fields are introduced. In this connection, the Hu-Washizu variational equation is invoked. To implement the analytical integration throughout the element, the modified ANS method is applied. As a result, the present finite rotation piezoelectric EG solid-shell element formulation permits the use of coarse meshes and very large load increments.