A four-node solid shell element formulation with assumed strain

A set of four-node shell element models based on the assumed strain formulation is considered here. The formulation allows for changes in the shell thickness. As a result, the kinematics of deformation are described by purely vectorial variables, without using rotational angles. The present study investigates the use of bubble function displacements and the assumed strain field. Careful selection of the assumed strain terms generates an element whose order of numerical integration does not increase even when the bubble function displacements are added. Results for the four-node element without any bubble function terms show sensitivity to element distortion. Use of the bubble functions with a carefully chosen assumed strain field greatly improves element performance. © 1998 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.     

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.     

An assumed strain triangular solid shell element with bubble function displacements for analysis of plates and shells

Bubble function displacements are used in conjunction with the assumed strain formulation to construct efficient triangular solid shell elements tailored for shell analysis. Two versions of 36‐DOF triangular elements are presented with different bubble function displacements and the corresponding assumed strain fields. In the first version the bubble function displacement is independent of the thickness while in the second version the bubble function varies linearly in the thickness direction. The assumed strain fields are carefully selected to alleviate locking effect. The first version models curved shells with flat triangular elements. The second version is effective in alleviating the membrane locking and thus allows more accurate modelling of curved shells. Various numerical examples demonstrate the validity and effectiveness of the present assumed strain formulation elements with the bubble function displacements. Copyright © 2001 John Wiley & Sons, Ltd.     

A co-rotational 8-node assumed strain shell element for postbuckling analysis of laminated composite plates and shells

 The formulation of a nonlinear composite shell element is presented for the solution of stability problems of composite plates and shells. The formulation of the geometrical stiffness presented here is exactly defined on the midsurface and is efficient for analyzing stability problems of thin and thick laminated plates and shells by incorporating bending moment and transverse shear resultant forces. The composite element is free of both membrane and shear locking behaviour by using the assumed natural strain method such that the element performs very well as thin shells. The transverse shear stiffness is defined by an equilibrium approach instead of using the shear correction factor. The proposed formulation is computationally efficient and the test results showed good agreement. In addition the effect of the viscoelastic material is investigated on the postbuckling behaviour of laminated composite shells. Received: 6 February 2002 / Accecpted: 6 January 2003 ID=" Present address: School of Civil Engineering, Asian Institute of Technology     

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.     

On using degenerated solid shell elements in adaptive refinement analysis

Three different degenerated shell elements are studied in an adaptive refinement procedure for the solution of shell problems. The stress recovery procedure expressed in a convective patch co‐ordinate system is used for the construction of continuous smoothed stress fields for the a posteriori error estimation. The performance of the stress recovery procedure, the error estimator and the adaptive refinement strategy are tested by solving three benchmark shell problems. It is found that when adaptive refinement is used, the adverse effects of boundary layers and stress singularities are eliminated and all the elements tested are able to achieve their optimal convergence rates. It is also found that the accuracy of the shell elements increases with the number of polynomial terms included in the stress and strain approximations. In addition, if complete Lagrangian polynomial terms are used, the element will be less sensitive to shape distortion than the one in which only complete polynomial terms are employed. Copyright © 1999 John Wiley & Sons, Ltd.     

Automatic adaptive refinement for plate bending problems using Reissner–Mindlin plate bending elements

The influence of the presence of singular points and boundary layers associated with the edge effects in a Reissner–Mindlin (RM) plate in the design of an optimal mesh for a finite element solution is studied, and methods for controlling the discretization error of the solution are suggested. An effective adaptive refinement strategy for the solution of plate bending problems based on the RM plate bending model is developed. This two-stage adaptive strategy is designed to control both the total and the shear error norms of a plate in which both singular points and boundary layers are present. A series of three different order assumed strain RM plate bending elements has been used in the adaptive refinement procedure. The locations of optimal sampling points and the effect of element shape distortions on the theoretical convergence rate of these elements are given and discussed. Numerical experiments show that the suggested refinement procedure is effective and that optimally refined meshes can be generated. It is also found that all the plate bending elements used can attain their full convergence rates regardless of the presence of singular points and boundary layers inside the problem domain. Boundary layer effects are well captured in all the examples tested and the use of a second stage of refinement to control the shear error is justified. In addition, tests on the Zienkiewicz–Zhu error estimator show that their performances are satisfactory. Finally, tests of the relative effectiveness of the plate bending elements used have also been made and it is found that while the higher order cubic element is the most accurate element tested, the quadratic element tested is the most efficient one in terms of CPU time used. © 1998 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.     

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.     

Free and forced vibration analysis of laminated composite plates and shells using a 9-node assumed strain shell element

The natural frequencies of isotropic and composite laminates are presented. The forced vibration analysis of laminated composite plates and shells subjected to arbitrary loading is investigated. In order to overcome membrane and shear locking phenomena, the assumed natural strain method is used. To develop a laminated shell element for free and forced vibration analysis, the equivalent constitutive equation that makes the computation of composite structures efficient was applied. The Mindlin-Reissner theory which allows the shear deformation and rotary inertia effect to be considered is adopted for development of nine-node assumed strain shell element. The present shell element offers significant advantages since it consistently uses the natural co-ordinate system. Results of the present theory show good agreement with the 3-D elasticity and analytical solutions. In addition the effect of damping is investigated on the forced vibration analysis of laminated composite plates and shells.     

A new one‐point quadrature enhanced assumed strain (EAS) solid‐shell element with multiple integration points along thickness—part II: nonlinear applications

In this work the recently proposed Reduced Enhanced Solid‐Shell (RESS) finite element, based on the enhanced assumed strain (EAS) method and a one‐point quadrature integration scheme, is extended in order to account for large deformation elastoplastic thin‐shell problems. One of the main features of this finite element consists in its minimal number of enhancing parameters (one), sufficient to circumvent the well‐known Poisson and volumetric locking phenomena, leading to a computationally efficient performance when compared to other 3D or solid‐shell enhanced strain elements. Furthermore, the employed numerical integration accounts for an arbitrary number of integration points through the thickness direction within a single layer of elements. The EAS formulation comprises an additive split of the Green–Lagrange material strain tensor, making the inclusion of nonlinear kinematics a straightforward task. A corotational coordinate system is used to integrate the constitutive law and to ensure incremental objectivity. A physical stabilization procedure is implemented in order to correct the element's rank deficiencies. A variety of shell‐type numerical benchmarks including plasticity, large deformations and contact are carried out, and good results are obtained when compared to well‐established formulations in the literature. Copyright © 2006 John Wiley & Sons, Ltd.     

Dual adaptive finite element refinement for multiple local quantities in linear elastostatics

In this paper, we summarize how dual analysis techniques can be used to determine upper bounds of the discretization error, both in terms of global and local outputs. We present formulas for the bounds of the error in local outputs, based on the approach proposed by Greenberg in 1948 and we show that the resulting intervals are the same as those previously presented, based on the approach proposed by Washizu in 1953. We then explain how the elemental contributions to these bounds can be used to define an adaptive strategy that considers multiple quantities and we present its application to a simple plane stress problem. Copyright © 2010 John Wiley & Sons, Ltd.     

Patch‐averaged assumed strain finite elements for stress analysis

A finite element model for linear‐elastic small deformation problems is presented. The formulation is based on a weighted residual that requires a priori the satisfaction of the kinematic equation. In this approach, an averaged strain‐displacement matrix is constructed for each node of the mesh by defining an appropriate patch of elements, yielding a smooth representation of strain and stress fields. Connections with traditional and similar procedure are explored. Linear quadrilateral four‐node and linear hexahedral eight‐node elements are derived. Various numerical tests show the accuracy and convergence properties of the proposed elements in comparison with extant finite elements and analytic solutions. Specific examples are also included to illustrate the ability to resist numerical locking in the incompressible limit and insensitive response in the presence of shape distortion. Furthermore, the numerical inf‐sup test is applied to a selection of problems to show the stability of the present formulation. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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.     

Two‐ and three‐dimensional transition element families for adaptive refinement analysis of elasticity problems

In this paper, new enhanced assumed strain (EAS) and hybrid stress transition element families are developed for 2D and 3D adaptive refinement analysis of elasticity problems. The EAS element families are based on some existing incompatible transition element families. By using the EAS method and the previous incompatible modes, the B ‐matrix columns associated with the EAS modes can be directly designed such that their domain integrals vanish automatically and they can be computed more efficiently. For 2D hybrid stress transition element families, it is possible to derive different stress fields that lead to rank‐sufficient transition elements. However, the task becomes intractable for 3D hybrid stress transition elements in which many combinations of mid‐side and mid‐face nodes are possible. This paper proposes to use hybrid stress transition element families in which the assumed stress fields are linearly complete. The new 2D element family is more accurate than the 2D rank‐sufficient element family. The new 3D element family is more accurate than the one with additional bilinear stress modes. Numerical examples reveal that the most accurate transition element families are the newly developed hybrid stress families followed by the EAS families, the incompatible families and then the compatible families. Copyright © 2008 John Wiley & Sons, Ltd.     

Three-node macro triangular shell element based on the assumed natural strains

 Finite element models with simple triangular geometry facilitate pre-processing for the finite element analysis. In the present study, a robust shear deformable triangular shell element formulation is presented for general plate and shell analysis, using three-node mesh discretization. The present formulation is developed on the basis of the assumed natural strain (ANS) formulation to attenuate the shear locking effect. Furthermore, this study proposes macro triangular element scheme by dividing a three-node triangular mesh into three parts to produce three individual ANS triangular elements and then merges these by condensing out the virtual center node. The performance of the macro ANS element is highly enhanced by reducing the number of sampling locations of three sub triangular element, but still utilizes three-node mesh discretization same as does the mesh used for three-node triangular elements. The macro ANS element has invariant stiffness, possesses no commutable zero energy mechanism. Numerical tests are presented to illustrate the high performance nature of the macro ANS element in general shell analysis. In particular, the numerical study demonstrates that the macro ANS element completely removes the shear locking effect that has been detrimental in shear deformable linear triangular elements so far. Received: 30 November 2001 / Accepted: 12 July 2002 This research was supported by the Ministry of Science and Technology through National Research Laboratory Programs under contract number 00-N-NL-01-C-026. The authors also thank for the financial support from the Brain Korea 21 Project in 2001.