Extension of the enhanced assumed strain method based on the structure of polyconvex strain-energy functions

In this work, two well-known approaches for mixed finite elements are combined to render three novel classes of elements. First, the widely used enhanced assumed strain (EAS) method is considered. Its key idea is to enhance a compatible kinematic field with an incompatible part. The second concept is a framework for mixed elements inspired by polyconvex strain-energy functions, in which the deformation gradient, its cofactor and determinant are three principal kinematic fields. The key idea for the novel elements is to treat enhancement of those three fields separately. This approach leads to a plethora of novel enhancement strategies and promising mixed finite elements. Some key properties of the newly proposed mixed approaches are that they are based on a Hu-Washizu type variational functional, fulfill the patch test, are frame-invariant, can be constructed completely locking free and show no spurious hourglassing in elasticity. Furthermore, they give additional insight into the mechanisms of standard EAS elements. Extensive numerical investigations are performed to assess the elements' behavior in elastic and elasto-plastic simulations.     

A highly efficient enhanced assumed strain physically stabilized hexahedral element

A method which combines the incompatible modes method with the physical stabilization method is developed to provide a highly efficient formulation for the single point eight‐node hexahedral element. The resulting element is compared to well‐known enhanced elements in standard benchmark type problems. It is seen that this single‐point element is nearly as coarse mesh accurate as the fully integrated EAS elements. A key feature is the novel enhanced strain fields which do not require any matrix inversions to solve for the internal element degrees of freedom. This, combined with the reduction of hourglass stresses to four hourglass forces, produces an element that is only 6.5 per cent slower than the perturbation stabilized single‐point brick element commonly used in many explicit finite element codes. Copyright © 2000 John Wiley & Sons, Ltd.     

Enhanced transverse shear strain shell formulation applied to large elasto‐plastic deformation problems

In this work, a previously proposed Enhanced Assumed Strain (EAS) finite element formulation for thin shells is revised and extended to account for isotropic and anisotropic material non‐linearities. Transverse shear and membrane‐locking patterns are successfully removed from the displacement‐based formulation. The resultant EAS shell finite element does not rely on any other mixed formulation, since the enhanced strain field is designed to fulfil the null transverse shear strain subspace coming from the classical degenerated formulation. At the same time, a minimum number of enhanced variables is achieved, when compared with previous works in the field. Non‐linear effects are treated within a local reference frame affected by the rigid‐body part of the total deformation. Additive and multiplicative update procedures for the finite rotation degrees‐of‐freedom are implemented to correctly reproduce mid‐point configurations along the incremental deformation path, improving the overall convergence rate. The stress and strain tensors update in the local frame, together with an additive treatment of the EAS terms, lead to a straightforward implementation of non‐linear geometric and material relations. Accuracy of the implemented algorithms is shown in isotropic and anisotropic elasto‐plastic problems. Copyright © 2004 John Wiley & Sons, Ltd.     

EAS concept for higher‐order finite shell elements to eliminate volumetric locking

The deficiency of volumetric locking phenomena in finite elements using higher‐order shell element formulations based on Lagrangean polynomials and a linear finite shell kinematics cannot be avoided by the existent enhanced assumed strain (EAS) concept established for low‐order elements. In this paper a consistent modification of the EAS concept is proposed to extend its applicability to higher‐order shell elements. This modification, affecting the transversal normal strain for polynomial orders p>1, eliminates pathological modes caused by volumetric locking. The efficiency of the proposed extended EAS method is demonstrated by means of eigenvalue analyses and two representative numerical examples. Copyright © 2008 John Wiley & Sons, Ltd.     

A study of stabilization and projection in the 4-node Mindlin plate element

In this paper, a Mindlin plate element is formulated based on the Hellinger–Reissner principle and the γ-technique. The stiffness consists of a constant stress (one-point quadrature) matrix and a stabilization matrix. The stabilization matrix is compared with those previously proposed. In addition, the element uses a projection to modify the nodal displacements so that the patch test is satisfied. The projection matrix is based on a mode decomposition. Several numerical cases are presented, and it is shown that the mode decomposition projection is necessary both for satisfaction of the patch test and convergence.     

土体三维弹塑性分析中的体积闭锁分析及单元选择

在材料不可压缩或胀/缩塑性流动情况下，传统低阶单元有体积闭锁问题。以摩尔－库仑模型为例，推导了塑性剪切应变和塑性体积应变的关系，揭示闭锁产生的原因。分析8节点等参元、Wilson非协调元、EAS单元和14节点单元的闭锁性态，表明8节点单元有严重闭锁性，Wilson非协调元也有闭锁性，EAS单元和采用降阶积分的14节点单元能克服闭锁。单元测试和方形基础的承载力计算两个数值算例证实了分析的结果，为土体三维分析中选择有效可靠的单元提供依据。     

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.     

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.     

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.     

Finite element formulations of strain gradient theory for microstructures and the C0–1 patch test

Based on finite element formulations for the strain gradient theory of microstructures, a convergence criterion for the C^0–1 patch test is introduced, and a new approach to devise strain gradient finite elements that can pass the C^0–1 patch test is proposed. The displacement functions of several plane triangular elements, which satisfy the C⁰ continuity and weak C¹ continuity conditions are evaluated by the C^0–1 patch test. The difference between the proposed C^0–1 patch test and the C⁰ constant stress and C¹ constant curvature patch tests is elucidated. An 18-DOF plane strain gradient triangular element (RCT9+RT9), which passes the C^0–1 patch test and has no spurious zero energy modes, is proposed. Numerical examples are employed to examine the performance of the proposed element by carrying out the C^0–1 patch test and eigenvalue test. The proposed element is found to be without spurious zero energy modes, and it possesses higher accuracy compared with other strain gradient elements. Copyright © 2004 John Wiley & Sons, Ltd.     

A novel versatile multilayer hybrid stress solid-shell element

This paper presents a versatile multilayer locking free hybrid stress solid-shell element that can be readily employed for a wide range of geometrically linear elastic structural analyses, i.e. from shell-like isotropic structures to multilayer anisotropic composites. This solid-shell element has eight nodes with only displacement degrees of freedom and a few internal parameters that provide the locking free behavior and accurate interlaminar stress resolution 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 composite laminates. The element formulation is based on the modified form of the well-known Fraeijs de Veubeke–Hu–Washizu (FHW) multifield variational principle with enhanced assumed strains (EAS formulation) and assumed natural strains (ANS 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 incorporating an assumed stress field in a standard EAS formulation based on the FHW principle. To assess the present formulation’s accuracy, a variety of popular numerical benchmark examples related to element patch tests, convergence, mesh distortion, shell and laminated composite analyses are investigated and the results are compared with those available in the literature. This assessment reveals that the proposed solid-shell formulation provides very accurate results for a wide range of structural analyses.     

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.     

IMPROVED ENHANCED STRAIN FOUR-NODE ELEMENT WITH TAYLOR EXPANSION OF THE SHAPE FUNCTIONS

A class of enhanced strain four-node elements with Taylor expansion of the shape function derivatives is presented. A new concept of enhancement using besides the ‘standard’ enhanced strain fields also two other enhanced fields is developed on the basis of the Hu–Washizu principle. For first-order Taylor expansion enhanced modes become uncoupled, thus only a negligible amount of computing effort for the static condensation of enhanced modes is needed. Furthermore, the formulation permits a symbolic integration, which leads to a closed-form solution for the element tangent matrix. Several numerical examples show that the element is stable, invariant, passes the patch test and yields good results especially in the highly distorted regime. © 1997 by John Wiley & Sons, Ltd.     

A set of pathological tests to validate new finite elements

The finite element method entails several approximations. Hence it is essential to subject all new finite elements to an adequate set of pathological tests in order to assess their performance. Many such tests have been proposed by researchers from time to time. We present an adequate set of tests, which every new finite element should pass. A thorough account of the patch test is also included in view of its significance in the validation of new elements.     

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.     

Development of shear locking‐free shell elements using an enhanced assumed strain formulation

The degenerated approach for shell elements of Ahmad and co‐workers is revisited in this paper. To avoid transverse shear locking effects in four‐node bilinear elements, an alternative formulation based on the enhanced assumed strain (EAS) method of Simo and Rifai is proposed directed towards the transverse shear terms of the strain field. In the first part of the work the analysis of the null transverse shear strain subspace for the degenerated element and also for the selective reduced integration (SRI) and assumed natural strain (ANS) formulations is carried out. Locking effects are then justified by the inability of the null transverse shear strain subspace, implicitly defined by a given finite element, to properly reproduce the required displacement patterns. Illustrating the proposed approach, a remarkably simple single‐element test is described where ANS formulation fails to converge to the correct results, being characterized by the same performance as the degenerated shell element. The adequate enhancement of the null transverse shear strain subspace is provided by the EAS method, enforcing Kirchhoff hypothesis for low thickness values and leading to a framework for the development of shear‐locking‐free shell elements. Numerical linear elastic tests show improved results obtained with the proposed formulation. Copyright © 2001 John Wiley & Sons, Ltd.     

Mindlin板弯曲和振动分析的高阶杂交应力三角形单元

现有的Mindlin板单元只能通过零剪力分片检验,而不能通过非零常剪力分片检验。该文根据Reissner- Mindlin一阶剪切变形理论,基于余能原理,提出了一种高阶杂交应力六节点三角形Mindlin板单元。该单元特点是不仅能通过零剪力分片检验,而且能通过严格的非零常剪力增强型分片检验。构造单元时特别注意了单元边界位移以及域内应力的插值函数的选取。采用任意阶Timoshenko梁函数作为边界位移插值函数,应力插值函数选取为满足平衡方程的多项式。对不同厚度不同边界条件的方板进行弯曲和自由振动分析,质量矩阵采用集中质量阵。数值结果表明无论对薄板还是中厚板,该单元均是准确有效的。     

Integration‐free Coons macroelements for the solution of 2D Poisson problems

Large isoparametric macroelements with closed‐form cardinal global shape functions under the label ‘Coons‐patch macroelements’ (CPM) have been previously proposed and used in conjunction with the finite element method and the boundary element method. This paper continues the research on the performance of CPM in conjunction with the collocation method. In contrast to the previous CPM that was based on a Galerkin/Ritz formulation, no domain integration is now required, a fact that justifies the name ‘integration‐free Coons macroelements’. Therefore, in addition to avoiding mesh generation, and saving human effort, the proposed technique has the additional advantage of further reducing the computer effort. The theory is supported by five test cases concerning Poisson and Laplace problems within 2D smooth quadrilateral domains. Copyright © 2008 John Wiley & Sons, Ltd.