Hybrid‐stress solid elements for shell structures based upon a modified variational functional

In this paper, we start with a modified generalized laminate stiffness matrix that serves as a remedy to resolve the thickness locking and some abnormalities encountered by solid‐shell elements in laminate analyses. A modified Hellinger–Reissner functional having displacement and a set of generalized stresses as independent fields is devised. Based upon the functional, eight‐node and 18‐node hybrid‐stress solid‐shell elements are proposed. A number of benchmark tests on homogenous and laminated plates/shells are conducted. The accuracy of the elements is promising. Copyright © 2002 John Wiley & Sons, Ltd.     

A hybrid stress ANS solid‐shell element and its generalization for smart structure modelling. Part I—solid‐shell element formulation

In the recent years, solid‐shell finite element models which possess no rotational degrees of freedom and applicable to thin plate/shell analyses have attracted considerable attention. Development of these elements are not straightforward. Shear, membrane, trapezoidal, thickness and dilatational lockings must been visioned. In this part of this paper, a novel eight‐node solid‐shell element is proposed. To resolve the shear and trapezoidal lockings, the assumed natural strain (ANS) method is resorted to. The hybrid‐stress formulation is employed to rectify the thickness and dilatational locking. The element is computationally more efficient than the conventional hybrid elements by adopting orthogonal‐assumed stress modes and enforcing admissible sparsity in the flexibility matrix. Popular benchmark tests are exercised to illustrate the efficacy of the elements. In Part II of the paper, the element will be generalized for smart structure modelling by including the piezoelectric effect. Copyright © 2000 John Wiley & Sons, Ltd.     

A hybrid stress ANS solid‐shell element and its generalization for smart structure modelling. Part II—smart structure modelling

In Part I of the paper, a hybrid‐stress‐assumed natural strain eight‐node solid‐shell element immune to shear, membrane, trapezoidal, thickness and dilatational lockings has been developed. Moreover, the element computational cost is reduced by enforcing admissible sparsity in the flexibility matrix. In this part of the paper, the solid‐shell element is generalized to a piezoelectric solid‐shell element. Using the two solid‐shell elements, smart structures with segmented piezoelectric sensors and actuators can be conveniently modelled. A number of problems are studied and comparisons with other ad hoc element models for smart structure modelling are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Curved quadratic triangular degenerated‐ and solid‐shell elements for geometric non‐linear analysis

Compared to the large number of curved quadrilateral degenerated‐ and solid‐shell elements, there are only a very few curved triangular degenerated‐ and solid‐shell elements. Based on the assumed natural strain sampling scheme previously developed for a quadratic degenerated‐shell element for linear analysis, this paper devises geometric non‐linear six‐node degenerated‐shell and twelve‐node solid‐shell elements. Both elements can be curved and are only equipped with the standard nodal d.o.f.s. Careful consideration has been exercised to circumvent various locking phenomena that plague degenerated‐ and solid‐shell elements. Numerical examples are presented to illustrate their efficacy. Copyright © 2003 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‐stress six‐node prismatic elements

This paper presents three novel hybrid‐stress six‐node prismatic elements. Starting from the element displacement interpolation, the equilibrating non‐constant stress modes for the first element are identified and orthogonalized with respect to the constant stress modes for higher computational efficiency. For the second element, the non‐constant stress modes are non‐equilibrating and chosen for the sake of stabilizing the reduced‐integrated element. The first two elements are intended for three‐dimensional continuum analysis with both passing the patch test for three‐dimensional continuum elements. The third element is primarily intended for plate/shell analysis. Shear locking is alleviated by a new assumed strain scheme which preserves the element accuracy with respect to the twisting load. Furthermore, the Poisson's locking along the in‐plane and out‐of‐plane directions is overcome by using the hybrid‐stress modes of the first element. The third element passes the patch test for plate/shell elements. Unless the element assumes the right prismatic geometry, it fails the patch test for three‐dimensional continuum elements. It will be seen that all the proposed elements are markedly more accurate than the conventional fully integrated element. Copyright © 2004 John Wiley & Sons, Ltd.     

A new reduced integration solid‐shell element based on EAS and ANS with hourglass stabilization

In this paper, a novel reduced integration eight‐node solid‐shell finite element formulation with hourglass stabilization is proposed. The enhanced assumed strain method is adopted to eliminate the well‐known volumetric and Poisson thickness locking phenomena with only one internal variable required. In order to alleviate the transverse shear and trapezoidal locking and correct rank deficiency simultaneously, the assumed natural strain method is implemented in conjunction with the Taylor expansion of the inverse Jacobian matrix. The projection of the hourglass strain‐displacement matrix and reconstruction of its transverse shear components are further employed to avoid excessive hourglass stiffness. The proposed solid‐shell element formulation successfully passes both the membrane and bending patch tests. Several typical examples are presented to demonstrate the excellent performance and extensive applicability of the proposed element. Copyright © 2015 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.     

An optimal versatile partial hybrid stress solid‐shell element for the analysis of multilayer composites

In the present contribution we propose an optimal low‐order versatile partial hybrid stress solid‐shell element that can be readily employed for a wide range of geometrically linear elastic structural analyses, that is, from shell‐like isotropic structures to multilayer anisotropic composites. This solid‐shell element has eight nodes with only displacement degrees of freedom and only a few internal parameters that provide the locking‐free behavior and accurate interlaminar shear stress resolution through the element thickness. These elements can be stacked on top of each other to model multilayer composite structures, fulfilling the interlaminar shear stress continuity at the interlayer surfaces and zero traction conditions on the top and bottom surfaces of composite laminates. The element formulation is based on the modified form of the well‐known Fraeijs de Veubeke–Hu–Washizu multifield 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 shear stress field in multilayer structures, which is achieved by the introduction of the assumed interlaminar shear stress field in a standard enhanced assumed strains formulation based on the Fraeijs de Veubeke–Hu–Washizu principle. The numerical testing of the present formulation, employing a variety of popular numerical benchmark examples related to element patch test, convergence, mesh distortion, shell and laminated composite analyses, proves its accuracy for a wide range of structural analyses.Copyright © 2012 John Wiley & Sons, Ltd.     

Development of geometrically exact new shell elements based on general curvilinear co‐ordinates

In the present study first‐order shear deformable shell finite elements based on general curvilinear co‐ordinates are proposed. For the development of the present shell elements, a partial mixed variational functional with independently assumed strains is provided in order to avoid the severe locking troubles known as transverse shear and membrane lockings. Bubble functions are included in the shape function of displacement to improve the performance of the developed element. The proposed assumed strain four‐ and nine‐node elements based on the general tensor shell theory provide an efficient linkage framework for shell surface modelling and finite element analysis. In the several benchmark problems, the present shell elements with exact geometric representations demonstrate their performance compared to previously reported results. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Quadratic NMS Mindlin‐plate‐bending element

The development of a robust and efficient quadratic Mindlin‐plate‐bending elements mainly by the use of non‐conforming displacement modes is presented in this paper. A brief review on the previous efforts to develop efficient non‐conforming Mindlin plate bending elements is also given. The behaviour of the newly proposed plate element is further improved by the combined use of nonconforming displacement modes, the selectively reduced integration scheme, and the assumed shear strain fields. Thus, the newly developed element has been designated as ‘NMS‐8P’. The improvement achieved may be attributable to the fact that the merits of these improvement techniques are merged in the formation of the new element in a complementary manner. The proposed 8‐node element passes the patch tests, does not show spurious mechanism, and does not produce shear locking phenomena even with distorted meshes. It is also shown that the element produces reliable solutions through numerical tests for standard benchmark problems. Copypright © 1999 John Wiley & Sons, Ltd.     

A co‐rotational grid‐based model for fabric drapes

Fabric drapes are typical large displacement, large rotation and small strain problems. Compared to conventional geometric non‐linear shell analyses, computational fabric drape analysis is particularly challenging due to the extremely weak bending rigidities of fabrics. Compared to continuum shell finite element methods, grid‐ or particle‐based methods appear to be more successful in high drapeability problems. The latter methods often resort to simple particle mechanics and formulate the elastic energy in terms of the inter‐particle distances and trigonometrical functions of the angles between the straight lines joining adjacent particles. In this paper, the co‐rotational approach and commonly employed assumptions for small strain problems in finite element analysis will be adopted to formulate the elastic energy. It will be seen that the internal force vector and the stiffness matrix are considerably simpler than other grid‐based models, yet the sparsity of the tangential stiffness matrix remains unchanged. A number of examples are considered and the predicted results are promising. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

High‐order layerwise finite element for the damped free‐vibration response of thick composite and sandwich composite plates

A high‐order layerwise finite element methodology is presented, which enables prediction of the damped dynamic characteristics of thick composite and sandwich composite plates. The through‐thickness displacement field in each discrete layer of the laminate includes quadratic and cubic polynomial distributions of the in‐plane displacements, in addition to the linear approximations assumed by linear layerwise theories. Stiffness, mass and damping matrices are formulated from ply to structural level. Interlaminar shear stress compatibility conditions are imposed on the discrete layer matrices, leading to both size reduction and prediction of interlaminar shear stresses at the laminate interfaces. The C¹ continuous finite element implemented yields an element damping matrix in addition to element stiffness and mass matrices. Application cases include thick [0/90/0], [±θ]_S and [±θ] composite plates with interlaminar damping layers and sandwich plates with composite faces and foam core. In the latter case, modal frequencies and damping were also experimentally determined and compared with the finite element predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

Exact static element stiffness matrices of non‐symmetric thin‐walled curved beams

An evaluation procedure of exact static stiffness matrices for curved beams with non‐symmetric thin‐walled cross section are rigorously presented for the static analysis. Higher‐order differential equations for a uniform curved beam element are first transformed into a set of the first‐order simultaneous ordinary differential equations by introducing 14 displacement parameters where displacement modes corresponding to zero eigenvalues are suitably taken into account. This numerical technique is then accomplished via a generalized linear eigenvalue problem with non‐symmetric matrices. Next, the displacement functions of displacement parameters are exactly calculated by determining general solutions of simultaneous non‐homogeneous differential equations. Finally an exact stiffness matrix is evaluated using force–deformation relationships. In order to demonstrate the validity and effectiveness of this method, displacements and normal stresses of cantilever thin‐walled curved beams subjected to tip loads are evaluated and compared with those by thin‐walled curved beam elements as well as shell elements. Copyright © 2004 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.     

An efficient co‐rotational formulation for curved triangular shell element

A 6‐node curved triangular shell element formulation based on a co‐rotational framework is proposed to solve large‐displacement and large‐rotation problems, in which part of the rigid‐body translations and all rigid‐body rotations in the global co‐ordinate system are excluded in calculating the element strain energy. Thus, an element‐independent formulation is achieved. Besides three translational displacement variables, two components of the mid‐surface normal vector at each node are defined as vectorial rotational variables; these two additional variables render all nodal variables additive in an incremental solution procedure. To alleviate the membrane and shear locking phenomena, the membrane strains and the out‐of‐plane shear strains are replaced with assumed strains in calculating the element strain energy. The strategy used in the mixed interpolation of tensorial components approach is employed in defining the assumed strains. The internal force vector and the element tangent stiffness matrix are obtained from calculating directly the first derivative and second derivative of the element strain energy with respect to the nodal variables, respectively. Different from most other existing co‐rotational element formulations, all nodal variables in the present curved triangular shell formulation are commutative in calculating the second derivative of the strain energy; as a result, the element tangent stiffness matrix is symmetric and is updated by using the total values of the nodal variables in an incremental solution procedure. Such update procedure is advantageous in solving dynamic problems. Finally, several elastic plate and shell problems are solved to demonstrate the reliability, efficiency, and convergence of the present formulation. Copyright © 2007 John Wiley & Sons, Ltd.     

Collapse of thin shell structures—stress resultant plasticity modelling within a co‐rotated ANDES finite element formulation

Due to the very non‐linear behaviour of thin shells under collapse, numerical simulations are subject to challenges. Shell finite elements are attractive in these simulations. Rotational degrees of freedom do, however, complicate the solution. In the present study a co‐rotated formulation is employed. The deformation of the shell is decomposed in to a contribution from large rigid body rotation and a strain producing term. A triangular assumed strain shell finite element is used. Hence, a high performance elastic element is combined with the co‐rotated formulation. In the co‐rotated co‐ordinate system the plasticity is accounted for by a simplifyed Ilyushin stress resultant yield surface. The stress update is determined from the backward Euler difference, and a consistent geometrical and material tangent stiffness is derived. Comparison with other published analysis results show that the present formulation gives acceptable accuracy. Copyright © 1999 John Wiley & Sons, Ltd.