A quadrilateral shell element with degree of freedom to represent thickness–stretch

This paper presents a quadrilateral shell element incorporating thickness–stretch, and demonstrates its performance in small and large deformation analyses for hyperelastic material and elastoplastic models. In terms of geometry, the proposed shell element is based on the formulation of the MITC4 shell element, with additional degrees of freedom to represent thickness–stretch. To consider the change in thickness, we introduce a displacement variation to the MITC4 shell element, in the thickness direction. After the thickness direction is expressed in terms of the director vectors that are defined at each midsurface node, additional nodes are placed along the thickness direction from the bottom surface to the top surface. The thickness–stretch is described by the movement of these additional nodes. The additional degrees of freedom are used to compute the transverse normal strain without assuming the plane stress condition. Hence, the three dimensional constitutive equation can be employed in the proposed formulation without any modification. By virtue of not imposing the plane stress condition, the surface traction is evaluated at the surface where the traction is applied, whereas it is assessed at the midsurface for conventional shell elements. Several numerical examples are presented to examine the fundamental performance of the proposed shell element. In particular, the proposed approach is capable of evaluating the change in thickness and the stress distribution when the effect of the surface traction is included. The behavior of the proposed shell element is compared with that of solid elements.

[1]  Eduardo N. Dvorkin,et al.  A formulation of the MITC4 shell element for finite strain elasto-plastic analysis , 1995 .

[2]  P. Flory,et al.  Thermodynamic relations for high elastic materials , 1961 .

[3]  M. Harnau,et al.  About linear and quadratic Solid-Shell elements at large deformations , 2002 .

[4]  Peter Wriggers,et al.  A note on enhanced strain methods for large deformations , 1996 .

[5]  Erasmo Carrera,et al.  Analysis of thickness locking in classical, refined and mixed multilayered plate theories , 2008 .

[6]  Erasmo Carrera,et al.  Analysis of thickness locking in classical, refined and mixed theories for layered shells , 2008 .

[7]  Michael Ortiz,et al.  A unified approach to finite deformation elastoplastic analysis based on the use of hyperelastic constitutive equations , 1985 .

[8]  E. Stein,et al.  A 4-node finite shell element for the implementation of general hyperelastic 3D-elasticity at finite strains , 1996 .

[9]  Marcelo Epstein,et al.  A finite element formulation for multilayered and thick plates , 1983 .

[10]  Anton Matzenmiller,et al.  Consistent discretization of thickness strains in thin shells including 3D‐material models , 1999 .

[11]  Thomas J. R. Hughes,et al.  Nonlinear finite element analysis of shells: Part I. three-dimensional shells , 1981 .

[12]  Eduardo N. Dvorkin,et al.  Implementation and stability analysis of the QMITC-TLH elasto-plastic finite strain (2D) element formulation , 2000 .

[13]  O. C. Zienkiewicz,et al.  Analysis of thick and thin shell structures by curved finite elements , 1970 .

[14]  R. Hauptmann,et al.  Extension of the ‘solid‐shell’ concept for application to large elastic and large elastoplastic deformations , 2000 .

[15]  J. C. Simo,et al.  A CLASS OF MIXED ASSUMED STRAIN METHODS AND THE METHOD OF INCOMPATIBLE MODES , 1990 .

[16]  J. C. Simo,et al.  Variational and projection methods for the volume constraint in finite deformation elasto-plasticity , 1985 .

[17]  Carlo Sansour,et al.  A theory and finite element formulation of shells at finite deformations involving thickness change: Circumventing the use of a rotation tensor , 1995, Archive of Applied Mechanics.

[18]  Klaus-Jürgen Bathe,et al.  3D-shell elements for structures in large strains , 2013 .

[19]  Dr.-Ing. C. Sansour.  A Theory and finite element formulation of shells at finite deformations involving thickness change , 1995 .

[20]  L. Vu-Quoc,et al.  Efficient and accurate multilayer solid‐shell element: non‐linear materials at finite strain , 2005 .

[21]  E. H. Lee,et al.  Finite‐Strain Elastic—Plastic Theory with Application to Plane‐Wave Analysis , 1967 .

[22]  David R. Owen,et al.  A refined analysis of laminated plates by finite element displacement methods—I. Fundamentals and static analysis , 1987 .

[23]  J. Degrieck,et al.  A novel versatile multilayer hybrid stress solid-shell element , 2013 .

[24]  R. M. Natal Jorge,et al.  An enhanced strain 3D element for large deformation elastoplastic thin-shell applications , 2004 .

[25]  H. Parisch.  A continuum‐based shell theory for non‐linear applications , 1995 .

[26]  Klaus-Jürgen Bathe,et al.  A 4-node 3D-shell element to model shell surface tractions and incompressible behavior , 2008 .

[27]  Ekkehard Ramm,et al.  Nonlinear shell formulations for complete three-dimensional constitutive laws including composites and laminates , 1994 .

[28]  Dominique Chapelle,et al.  3D-SHELL ELEMENTS AND THEIR UNDERLYING MATHEMATICAL MODEL , 2004 .

[29]  M. Ortiz,et al.  Formulation of implicit finite element methods for multiplicative finite deformation plasticity , 1990 .

[30]  Christian Miehe,et al.  A theoretical and computational model for isotropic elastoplastic stress analysis in shells at large strains , 1998 .

[31]  R. Hauptmann,et al.  `Solid-shell' elements with linear and quadratic shape functions at large deformations with nearly incompressible materials , 2001 .

[32]  E. Dvorkin.  Nonlinear analysis of shells using the MITC formulation , 1995 .

[33]  Klaus-Jürgen Bathe,et al.  A hyperelastic‐based large strain elasto‐plastic constitutive formulation with combined isotropic‐kinematic hardening using the logarithmic stress and strain measures , 1990 .

[34]  Stefanie Reese,et al.  A large deformation solid‐shell concept based on reduced integration with hourglass stabilization , 2007 .

[35]  Sven Klinkel,et al.  A continuum based three-dimensional shell element for laminated structures , 1999 .

[36]  N. El-Abbasi,et al.  A Continuum Based Thick Shell Element for Large Deformation Analysis of Layered Composites , 2005 .

[37]  T. Hughes.  Generalization of selective integration procedures to anisotropic and nonlinear media , 1980 .

[38]  Hashem M. Mourad,et al.  A global–local discontinuous Galerkin shell finite element for small-deformation analysis of multi-layered composites , 2014 .

[39]  Sven Klinkel,et al.  A robust non-linear solid shell element based on a mixed variational formulation , 2006 .

[40]  Dong Yol Yang,et al.  Springback prediction for sheet metal forming process using a 3D hybrid membrane/shell method , 2002 .

[41]  Peter M. Pinsky,et al.  A multi-director formulation for nonlinear elastic-viscoelastic layered shells , 1986 .

[42]  Eduardo N. Dvorkin,et al.  A finite element formulation for finite strain elasto-plastic analysis based on mixed interpolation of tensorial components , 1994 .

[43]  K. Bathe,et al.  A continuum mechanics based four‐node shell element for general non‐linear analysis , 1984 .

[44]  Peter Wriggers,et al.  Finite element concepts for finite elastoplastic strains and isotropic stress response in shells: theoretical and computational analysis , 1999 .

[45]  Ted Belytschko,et al.  Resultant-stress degenerated-shell element , 1986 .

[46]  Sven Klinkel,et al.  A mixed shell formulation accounting for thickness strains and finite strain 3d material models , 2008 .

[47]  T. Hughes,et al.  Nonlinear finite element shell formulation accounting for large membrane strains , 1983 .

[48]  F. Kikuchi,et al.  A 9-node mixed shell element based on the Hu-Washizu principle , 1991 .

[49]  Noritoshi Iwata,et al.  Numerical Prediction of Springback Shape of Severely Bent Sheet Metal , 2007 .

[50]  Peter Betsch,et al.  Numerical implementation of multiplicative elasto-plasticity into assumed strain elements with application to shells at large strains , 1999 .

[51]  E. Ramm,et al.  Three‐dimensional extension of non‐linear shell formulation based on the enhanced assumed strain concept , 1994 .

[52]  E. Stein,et al.  An assumed strain approach avoiding artificial thickness straining for a non‐linear 4‐node shell element , 1995 .

[53]  R. Hauptmann,et al.  A SYSTEMATIC DEVELOPMENT OF 'SOLID-SHELL' ELEMENT FORMULATIONS FOR LINEAR AND NON-LINEAR ANALYSES EMPLOYING ONLY DISPLACEMENT DEGREES OF FREEDOM , 1998 .

[54]  K. Bathe.  Finite Element Procedures , 1995 .

[55]  L. Vu-Quoc,et al.  Geometrically exact sandwich shells: The dynamic case , 2000 .

[56]  Erasmo Carrera,et al.  Reissner’s mixed variational theorem toward MITC finite elements for multilayered plates , 2013 .

[57]  J. C. Simo,et al.  On a stress resultant geometrically exact shell model , 1990 .

[58]  Multilayer shells: Geometrically-exact formulation of equations of motion , 2000 .

[59]  S. A. Meguid,et al.  A new shell element accounting for through-thickness deformation , 2000 .

[60]  Ekkehard Ramm,et al.  EAS‐elements for two‐dimensional, three‐dimensional, plate and shell structures and their equivalence to HR‐elements , 1993 .

[61]  P. Wriggers,et al.  A fully nonlinear multi-parameter shell model with thickness variation and a triangular shell finite element , 2004 .

[62]  J. C. Simo,et al.  On stress resultant geometrically exact shell model. Part I: formulation and optimal parametrization , 1989 .

[63]  Results on best theories for metallic and laminated shells including Layer-Wise models , 2015 .