An improved solid‐shell element based on ANS and EAS concepts

This paper presents an eight‐node nonlinear solid‐shell element for static problems. The main goal of this work is to develop a solid‐shell formulation with improved membrane response compared with the previous solid‐shell element (MOS2013), presented in 1 . Assumed natural strain concept is implemented to account for the transverse shear and thickness strains to circumvent the curvature thickness and transverse shear locking problems. The enhanced assumed strain approach based on the Hu–Washizu variational principle with six enhanced assumed strain degrees of freedom is applied. Five extra degrees of freedom are applied on the in‐plane strains to improve the membrane response and one on the thickness strain to alleviate the volumetric and Poisson's thickness locking problems. The ensuing element performs well in both in‐plane and out‐of‐plane responses, besides the simplicity of implementation. The element formulation yields exact solutions for both the membrane and bending patch tests. The formulation is extended to the geometrically nonlinear regime using the corotational approach, explained in 2 . Numerical results from benchmarks show the robustness of the formulation in geometrically linear and nonlinear problems. Copyright © 2016 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.     

Mixed shell element for seven‐parameter formulation

A new mixed shell element is developed for a seven‐parameter formulation in this paper. The mixed shell element is constructed by assuming stress field and displacement field together. Assumed stress field and assumed displacement field can be combined by stress–strain relationship with Hu‐Washizu functional. The developed mixed shell element can provide more flexible stiffness than other commercial softwares. Additionally, seven‐parameter shell formulation is used instead of Reissner/Mindlin formulation, since it can provide the thickness change. Even though some commercial engineering software are not proper for very thick shell structure, the developed mixed shell element for seven‐parameter formulation can be used without distinction of thick shell and thin shell. An example of shell models with different thickness is provided with solid model. Static and modal analyses are also performed for verification. Copyright © 2005 John Wiley & Sons, Ltd.     

A rotation‐free shell formulation using nodal integration for static and dynamic analyses of structures

This paper proposed a rotation‐free thin shell formulation with nodal integration for elastic–static, free vibration, and explicit dynamic analyses of structures using three‐node triangular cells and linear interpolation functions. The formulation is based on the classic Kirchhoff plate theory, in which only three translational displacements are treated as the filed variables. Based on each node, the integration domains are further formed, where the generalized gradient smoothing technique and Green divergence theorem that can relax the continuity requirement for trial function are used to construct the curvature filed. With the aid of strain smoothing operation and tensor transformation rule, the smoothed strains in the integration domain can be finally expressed by constants. The principle of virtual work is then used to establish the discretized system equations. The translational boundary conditions are imposed same as the practice of standard finite element method, while the rotational boundary conditions are constrained in the process of constructing the smoothed curvature filed. To test the performance of the present formulation, several numerical examples, including both benchmark problems and practical engineering cases, are studied. The results demonstrate that the present method possesses better accuracy and higher efficiency for both static and dynamic problems. Copyright © 2015 John Wiley & Sons, Ltd.     

A 4‐node assumed strain quasi‐conforming shell element with 6 degrees of freedom

Quasi‐conforming formulations of 4‐node stress‐resultant shell elements are presented. The element formulations use interrelated displacement–rotation interpolations. The formulation also includes drilling degrees of freedom, which improves membrane behavior and allows the modeling of stiffened plates and shells. The proposed treatment for bending provides very good results in the 4‐node shell element. The stiffness matrices for the present elements are explicitly expressed and the stresses are taken accurately at the nodal points. Compared to elements using Gauss integration, where the stresses are most accurate at the integration points, the extrapolation procedure needed for post‐processing is eliminated in the present shell element. A lot of numerical tests were carried out for the validation of the present 4‐node shell element and the results are in good agreement with references. Copyright © 2003 John Wiley & Sons, Ltd.     

A large deformation solid‐shell concept based on reduced integration with hourglass stabilization

In this paper a new eight‐node (brick) solid‐shell finite element formulation based on the concept of reduced integration with hourglass stabilization is presented. The work focuses on static problems. The starting point of the derivation is the three‐field variational functional upon which meanwhile established 3D enhanced strain concepts are based. Important additional assumptions are made to transfer the approach into a powerful solid‐shell. First of all, a Taylor expansion of the first Piola–Kirchhoff stress tensor with respect to the normal through the centre of the element is carried out. In this way the stress becomes a linear function of the shell surface co‐ordinates whereas the dependence on the thickness co‐ordinate remains non‐linear. Secondly, the Jacobian matrix is replaced by its value in the centre of the element. These two assumptions lead to a computationally efficient shell element which requires only two Gauss points in the thickness direction (and one Gauss point in the plane of the shell element). Additionally three internal element degrees‐of‐freedom have to be determined to avoid thickness locking. One important advantage of the element is the fact that a fully three‐dimensional stress state can be modelled without any modification of the constitutive law. The formulation has only displacement degrees‐of‐freedom and the geometry in the thickness direction is correctly displayed. Copyright © 2006 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.     

A finite‐strain quadrilateral shell element based on discrete Kirchhoff–Love constraints

This paper improves the 16 degrees‐of‐freedom quadrilateral shell element based on pointwise Kirchhoff–Love constraints and introduces a consistent large strain formulation for this element. The model is based on classical shell kinematics combined with continuum constitutive laws. The resulting element is valid for large rotations and displacements. The degrees‐of‐freedom are the displacements at the corner nodes and one rotation at each mid‐side node. The formulation is free of enhancements, it is almost fully integrated and is found to be immune to locking or unstable modes. The patch test is satisfied. In addition, the formulation is simple and amenable to efficient incorporation in large‐scale codes as no internal degrees‐of‐freedom are employed, and the overall calculations are very efficient. Results are presented for linear and non‐linear problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Sandwich shell finite element for dynamic explicit analysis

The contribution of this paper consists of new development of transverse shear stresses through the thickness and finding an expression for the critical time step for explicit time integration of layered shells. This work presents the finite element (FE) formulation and implementation of a higher‐order shear deformable shell element for dynamic explicit analysis of composite and sandwich shells. The formulation is developed using a displacement‐based third‐order shear deformation shell theory. Using the differential equilibrium equations and the interlayer requirements, special treatment is developed for the transverse shear, resulting in a continuous, piecewise quartic distribution of the transverse shear stresses through the shell thickness. Expressions are developed for the critical time step of the explicit time integration for orthotropic homogeneous and layered shells based on the developed third‐order formulation. To assess the performance of the present shell element, it is implemented in the general non‐linear explicit dynamic FE code DYNA3D. Several problems are solved and results are presented and compared to other theoretical and numerical results. Copyright © 2002 John Wiley & Sons, Ltd.     

A four‐node,shear‐deformable shell element developed via explicit Kirchhoff constraints

An efficient, four‐node quadrilateral shell element is formulated using a linear, first‐order shear deformation theory. The bending part of the formulation is constructed from a cross‐diagonal assembly of four three‐node anisoparametric triangular plate elements, referred to as MIN3. Closed‐form constraint equations, which arise from the Kirchhoff constraints in the thin‐plate limit, are derived and used to eliminate the degrees‐of‐freedom associated with the ‘internal’ node of the cross‐diagonal assembly. The membrane displacement field employs an Allman‐type, drilling degrees‐of‐freedom formulation. The result is a displacement‐based, fully integrated, four‐node quadrilateral element, MIN4T, possessing six degrees‐of‐freedom at each node. Results for a set of validation plate problems demonstrate that the four‐node MIN4T has similar robustness and accuracy characteristics as the original cross‐diagonal assembly of MIN3 elements involving five nodes. The element performs well in both moderately thick and thin regimes, and it is free of shear locking. Shell validation results demonstrate superior performance of MIN4T over MIN3, possibly as a result of its higher‐order interpolation of the membrane displacements. It is also noted that the bending formulation of MIN4T is kinematically compatible with the existing anisoparametric elements of the same order of approximation, which include a two‐node Timoshenko beam element and a three‐node plate element, MIN3. Copyright © 2000 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.     

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.     

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.     

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.     

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.     

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.     

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 continuum eight‐parameter shell finite element for large deformation analysis

Abstract

In this paper, a finite element formulation, using eight independent parameters and high‐order spectral/hp functions, for nonlinear analysis is presented. This formulation allows the use of a third‐order thickness stretch kinematics, which also avoids Poisson's locking. Full nonlinear terms up to quadratic in the Green–Lagrange strain tensorare retained. Several nontrivial problems are solved using the presented formulation. A comparison between this formulation and others found in the literature,and with shell and solid elements in commercial codes ABAQUS and ANSYS are presented and the differences are brought out.