A generalized displacement method for the finite element analysis of thin shells

The shear locking problem for the bilinear degenerated thick shell elements, when used in the context of thin shell structures, can be overcome by a generalized displacement method presented in this paper. The transverse shear energy in the degenerated thick shell elements is totally suppressed by introducing discrete Kirchhoff constraints in each element. The constrained variational problem based on the nodal displacement space is transformed into an unconstrained one based on a so-called generalized displacement subspace. It is shown that shear looking phenomena completely disappear and no degradation of results is observed as the ratio of thickness to span approaches zero.     

Refined discrete quadrilateral degenerated shell element by using Timoshenko's beam function

A refined discrete degenerated 20‐DOF quadrilateral shell element RQS20 is proposed. The exact displacement function of the Timoshenko's beam is used as the displacement on the element boundary. The re‐constitute method for shear strain matrix is adopted. The proposed element can be used for the analysis of both moderately thick and thin plates/shells, and the convergence for the very thin case can be ensured theoretically. Numerical examples presented show that the new model indeed possesses higher accuracy in the analysis of thin and thick plates/shells, and that it can pass the patch test required for the Kirchhoff thin plate elements. Most important of all, it is free from the membrane and shear locking phenomena for extremely thin plates/shells, on the one hand, and it can also avoid the phenomenon of oscillatory solutions for thick plates/shells case on the other. Copyright © 2005 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 triangular finite element for thin plates and shells

A finite element modelling technique which utilizes a triangular element with 45 degrees-of-freedom and seven-point integration has been tested for analysis of thin plate and shell structures. The element is based on the degenerate solid shell concept and the mixed formulation with assumed independent inplane and transverse shear strains. Numerical result indicates effectiveness of the present modelling technique which features combined use of elements with kinematic modes and those without kinematic modes in an attempt to eliminate both locking and spurious kinematic modes at global structural level.     

A 3‐node flat triangular shell element with corner drilling freedoms and transverse shear correction

The formulation, implementation and testing of simple, efficient and robust shell finite elements have challenged investigators over the past four decades. A new 3‐node flat triangular shell element is developed by combination of a membrane component and a plate bending component. The ANDES‐based membrane component includes rotational degrees of freedom, and the refined nonconforming element method‐based bending component involves a transverse shear correction. Numerical examples are carried out for benchmark tests. The results show that compared with some popular shell elements, the present one is simple but exhibits excellent all‐around properties (for both membrane‐and bending‐dominated situations), such as free of aspect ratio locking, passing the patch test, free of shear locking, good convergence and high suitability for thin to moderately thick plates. The developed element has already been adopted in a warpage simulation package for injection molding. Copyright © 2011 John Wiley & Sons, Ltd.     

Refined triangular discrete Mindlin flat shell elements

In this paper flat shell elements are formed by the assemblage of discrete Mindlin plate elements RDKTM and either the constant strain membrane element CST or the Allmans membrane element with drilling degrees of freedom LST. The element RDKTM is a robust Mindlin plate element, which can perform uniformly thick and thin plate bending analysis. It also passes the patch test for thin plate bending, and its convergence for very thin plates can be ensured theoretically. The singularity of the stiffness matrix and membrane locking are studied for the present elements. Numerical examples are presented to show that the present models indeed possess properties of simple formulations, high accuracy for thin and thick shells, and it is free from shear locking for thin plate/shell analysis.     

THE FOUR-NODE C0 SHELL ELEMENT REFORMULATED

A reformulated four-node shell element, based on the analysis of the moment redistribution mechanism development by C⁰ plate bending and shell elements, is presented. The moment redistribution mechanism of a finite shell element model is shown to be predominantly activated by the membrane flexural action of the shell. This action is triggered through the membrane strain components which participate in the moment equilibrium equations of the finite element assembly system. An equivalent elastic foundation action, along with the activation of the in-plane twisting stiffness of the shell, may also contribute to the moment redistribution mechanism of the finite shell element model. The proposed shell element formulation aims at retaining the non-spurious contribution of the transverse shear/membrane strain energy to the flexural behaviour of the shell, through the activation of the moment redistribution mechanism. Yet, any potentially spurious, whether locking or kinematic, mechanism is rejected. In warped configurations, the element activates appropriate coupling mechanisms of the bending terms to nodal translations. The so-obtained reformulated four-node shell element exhibits an excellent behaviour without experiencing any locking phenomena or zero-energy modes, while its formulation is kept simple, based on physical considerations. The proposed formulation performs equally well in flat as well as in warped shell element applications.     

A new one‐point quadrature enhanced assumed strain (EAS) solid‐shell element with multiple integration points along thickness: Part I—geometrically linear applications

Accuracy and efficiency are the main features expected in finite element method. In the field of low‐order formulations, the treatment of locking phenomena is crucial to prevent poor results. For three‐dimensional analysis, the development of efficient and accurate eight‐node solid‐shell finite elements has been the principal goal of a number of recent published works. When modelling thin‐ and thick‐walled applications, the well‐known transverse shear and volumetric locking phenomena should be conveniently circumvented. In this work, the enhanced assumed strain method and a reduced in‐plane integration scheme are combined to produce a new eight‐node solid‐shell element, accommodating the use of any number of integration points along thickness direction. Furthermore, a physical stabilization procedure is employed in order to correct the element's rank deficiency. Several factors contribute to the high computational efficiency of the formulation, namely: (i) the use of only one internal variable per element for the enhanced part of the strain field; (ii) the reduced integration scheme; (iii) the prevention of using multiple elements' layers along thickness, which can be simply replaced by any number of integration points within a single element layer. Implementation guidelines and numerical results confirm the robustness and efficiency of the proposed approach when compared to conventional elements well‐established in the literature. 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.     

A new axi-symmetric element for thin walled structures

A new axi-symmetric finite element for thin walled structures is presented in this work. It uses the solid-shell element’s concept with only a single element and multiple integration points along the thickness direction. The cross-section of the element is composed of four nodes with two degrees of freedom each. The proposed formulation overcomes many locking pathologies including transverse shear locking, Poisson’s locking and volumetric locking. For transverse shear locking, the formulation uses the selective reduced integration technique, for Poisson’s locking it uses the enhanced assumed strain (EAS) method with only one enhancing variable. The B-bar approach is used to eliminate the isochoric deformations in the hourglass field while the EAS method is used to alleviate the volumetric locking in the constant part of the deformation tensor. Several examples are shown to demonstrate the performance and accuracy of the proposed element with special focus on the numerical simulations for the beverage can industry.     

Extended rotation-free shell triangles with transverse shear deformation effects

The paper extends recent work of the authors to include transverse shear effects on rotation-free triangular element for plates (O?ate and Zárate in Int J Numer Methods Eng 83(2):196–227, 2010). Two new shell triangular elements are presented, the EBST+ and the EBST+1. Transverse shear deformation effects are important for thick shells, as well when the shell is laminated or formed by composite material. The ingredients for the element formulation are: a Hu-Washizu type mixed functional and linear interpolation for the displacement field. In both elements presented a finite volume approach is used for computing the bending moments and the curvatures over a patch of elements. The nodal translational degrees of freedom of the original enhanced basic shell triangle (EBST) are extended with the two shear deformation angles via two different approaches. The first one uses a linear interpolation of the rotation angles inside the element (EBST+) and the second one assumes a constant field for the rotation angles (EBST+1). For the thin shell case the shear angles vanish and the new elements reproduce the good behaviour of the original thin EBST element. As a consequence the elements can reproduce the solutions for thick to thin shells situations without exhibiting shear locking. The numerical solution for the thick shell case can be found iteratively starting from the deflection values for the Kirchhoff theory using the original thin EBST element. Examples of the good performance of the new rotation-free shell triangles are given.     

A new nine node degenerated shell element with enhanced membrane and shear interpolation

A new nine node degenerated shell element is presented in this paper. In the formulation of the new element, an enhanced interpolation of the transverse shear strains in the natural co-ordinate system is used to overcome the shear locking problem, and an enhanced interpolation of the membrane strains in the local Cartesian co-ordinate system is applied to avoid membrane locking behaviour. It is shown that the resulting element has the requisite number of zero eigenvalues and associated rigid body modes. The element does not exhibit membrane or shear locking for large span/thickness ratios. To illustrate the good performance of the new element some examples are presented including comparisons with the behaviour of the selectively integrated Lagrangian degenerated shell elements.     

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.     

Locking of thin plate/shell elements

Often, finite element solutions of thin plate/shell elements become very stiff and the displacement field solutions diverge from those predicted by Kirchhoff's theory. This phenomenon is known as the locking phenomenon. A theoretical fomulation demonstrating its existence is developed, and results of finite element analysis of a single element and mesh are discussed. This leads to a sufficient and necessary criterion which must be satisfied to avoid the locking phenomenon.     

拟协调九结点四边形退化壳单元

本文根据胡－鹫津变分原理和退化壳有限元概念，构造了一个拟协调（或称杂交／混合）九结点四边形退化壳单元ＣＳＨ９．该单元克服了剪切闭锁和膜闭锁现象，整体上无零能模式，通过了分片检验。数值算例表明该单元是相当准确的和有效的。     

Reduced integration and the shear-flexible beam element

A clearer insight into the ‘shear locking’ phenomenon, which appears in the development of C₀ continuous element using shear-flexible or penalty type formulations, is obtained by a careful study of the Timoshenko beam element. When a penalty type argument is used to degenerate thick elements to thin elements, the various approximations of the shear related energy terms act as different types of constraints and, depending on the formulation, two types of constraints which are classified as true or spurious may emerge. The spurious constraints, where they exist, are responsible for the ‘shear locking’ phenomenon, and its manifestation and elimination is demonstrated in a very simple example. The source of difficulty is shown to be the mathematical operations involved in the various shape function definitions and subsequent integration of functionals. It is seen that formulations that ensure only true constraints in the extreme penalty limit cases display far superior performance in the thick element situation as well, and thus guidelines for the development of efficient elements are drawn. A similar type of behaviour is observed in a shallow curved beam element and here ‘inplane locking’ can be eliminated by selective integration to obtain an improved curved beam element. However, ‘inplane locking’ does not cause a spurious constraint as the error quickly vanishes with the reduction of element size for a reasonable radius of curvature conforming with shallow shell theory.     

An improved EAS brick element for finite deformation

A new enhanced assumed strain brick element for finite deformations in finite elasticity and plasticity is presented. The element is based on an expansion of shape function derivatives using Taylor series and an extended set of orthogonality conditions that have to be satisfied for an hourglassing free EAS formulation. Such approach has not been applied so far in the context of large deformation three-dimensional problems. It leads to a surprisingly well-behaved locking and hourglassing free element formulation. Major advantage of the new element is its shear locking free performance in the limit of very thin elements, thus it is applicable to shell type problems. Crucial for the derivation of the residual and consistent tangent matrix of the element is the automation of the implementation by automatic code generation.     

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.