On the use of a reduced enhanced solid-shell (RESS) element for sheet forming simulations

Abstract A recently proposed reduced enhanced solid-shell (RESS) element [Alves de Sousa, R.J., Cardoso, R.P.R., Fontes Valente, R.A., Yoon, J.W., Gracio, J.J., Natal Jorge, R.M., 2005. A new one-point quadrature enhanced assumed strain (EAS) solid-shell element with multiple integration points along thickness: Part I – Geometrically Linear Applications. International Journal for Numerical Methods in Engineering 62, 952–977; Alves de Sousa, R.J., Cardoso, R.P.R., Fontes Valente, R.A., Yoon, J.W., Gracio, J.J., Natal Jorge, R.M., 2006. A new one-point quadrature enhanced assumed strain (EAS) solid-shell element with multiple integration points along thickness: Part II – Nonlinear Applications. International Journal for Numerical Methods in Engineering, 67, 160–188.] is based on the enhanced assumed strain (EAS) method with a one-point quadrature numerical integration scheme. In this work, the RESS element is applied to large-deformation elasto-plastic thin-shell applications, including contact and plastic anisotropy. One of the main advantages of the RESS is its minimum number of enhancing parameters (only one), which when associated with an in-plane reduced integration scheme, circumvents efficiently well-known locking phenomena, leading to a computationally efficient performance when compared to conventional 3D solid elements. It is also worth noting that the element accounts for an arbitrary number of integration points through thickness direction within a single element layer. This capability has proven to be efficient, for instance, for accurately describing springback phenomenon in sheet forming simulations. A physical stabilization procedure is employed in order to correct the element’s rank deficiency. A general elasto-plastic model is also incorporated for the constitutive modelling of sheet forming operations with plastic anisotropy. Several examples including contact, anisotropic plasticity and springback effects are carried out and the results are compared with experimental data.

[1]  A. P. Karafillis,et al.  A general anisotropic yield criterion using bounds and a transformation weighting tensor , 1993 .

[2]  O. C. Zienkiewicz,et al.  Reduced integration technique in general analysis of plates and shells , 1971 .

[3]  Michael A. Puso,et al.  A highly efficient enhanced assumed strain physically stabilized hexahedral element , 2000 .

[4]  R. H. Wagoner,et al.  A multiplicative finite elasto-plastic formulation with anisotropic yield functions , 2003 .

[5]  Stefanie Reese,et al.  A new stabilization technique for finite elements in non-linear elasticity , 1999 .

[6]  Jeong Whan Yoon,et al.  A new one‐point quadrature enhanced assumed strain (EAS) solid‐shell element with multiple integration points along thickness—part II: nonlinear applications , 2006 .

[7]  Christian Miehe,et al.  Anisotropic elastic–plastic analysis of shells at large strains. A comparison of multiplicative and additive approaches to enhanced finite element design and constitutive modelling , 2004 .

[8]  J. C. Simo,et al.  Geometrically non‐linear enhanced strain mixed methods and the method of incompatible modes , 1992 .

[9]  Peter Wriggers,et al.  A stabilization technique to avoid hourglassing in finite elasticity , 2000 .

[10]  Frédéric Barlat,et al.  Generalization of Drucker's Yield Criterion to Orthotropy , 2001 .

[11]  W. Hosford.  A Generalized Isotropic Yield Criterion , 1972 .

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

[13]  Renato Natal Jorge,et al.  Development of shear locking‐free shell elements using an enhanced assumed strain formulation , 2002 .

[14]  Dong-Yol Yang,et al.  Elasto-plastic finite element method based on incremental deformation theory and continuum based shell elements for planar anisotropic sheet materials , 1999 .

[15]  K. Chung,et al.  Finite element simulation of sheet metal forming for planar anisotropic metals , 1992 .

[16]  S. Reese.  On the Equivalent of Mixed Element Formulations and the Concept of Reduced Integration in Large Deformation Problems , 2002 .

[17]  Frédéric Barlat,et al.  Development of a one point quadrature shell element for nonlinear applications with contact and anisotropy , 2002 .

[18]  Arif Masud,et al.  A stabilized 3-D co-rotational formulation for geometrically nonlinear analysis of multi-layered composite shells , 2000 .

[19]  Jacques Besson,et al.  A yield function for anisotropic materials Application to aluminum alloys , 2004 .

[20]  Frédéric Barlat,et al.  Linear transfomation-based anisotropic yield functions , 2005 .

[21]  Manabu Gotoh,et al.  A finite element analysis of rigid-plastic deformation of the flange in a deep-drawing process based on a fourth-degree yield function , 1978 .

[22]  Antoine Legay,et al.  Elastoplastic stability analysis of shells using the physically stabilized finite element SHB8PS , 2003 .

[23]  F. Barlat,et al.  Anisotropic strain hardening behavior in simple shear for cube textured aluminum alloy sheets , 2005 .

[24]  S. Reese.  On a physically stabilized one point finite element formulation for three-dimensional finite elasto-plasticity , 2005 .

[25]  Frédéric Barlat,et al.  A general elasto-plastic finite element formulation based on incremental deformation theory for planar anisotropy and its application to sheet metal forming , 1999 .

[26]  Rainald Löhner,et al.  A general advancing front technique for filling space with arbitrary objects , 2004 .

[27]  Medhat A. Haroun,et al.  Reduced and selective integration techniques in the finite element analysis of plates , 1978 .

[28]  R. M. Natal Jorge,et al.  A new volumetric and shear locking‐free 3D enhanced strain element , 2003 .

[29]  K. Chung,et al.  A deformation theory of plasticity based on minimum work paths , 1993 .

[30]  Jeong Whan Yoon,et al.  A new one‐point quadrature enhanced assumed strain (EAS) solid‐shell element with multiple integration points along thickness: Part I—geometrically linear applications , 2005 .

[31]  R. E. Dick,et al.  Plane stress yield functions for aluminum alloy sheets , 2002 .

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

[33]  F. Barlat,et al.  Plane stress yield function for aluminum alloy sheets—part 1: theory , 2003 .

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

[35]  Dong-Yol Yang,et al.  Earing predictions based on asymmetric nonquadratic yield function , 2000 .

[36]  Frédéric Barlat,et al.  Plastic behavior and stretchability of sheet metals. Part I: A yield function for orthotropic sheets under plane stress conditions , 1989 .

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

[38]  S. Reese,et al.  A new locking-free brick element technique for large deformation problems in elasticity ☆ , 2000 .

[39]  Yeong-Maw Hwang,et al.  Analysis of tube hydroforming in a square cross-sectional die , 2005 .

[40]  Wing Kam Liu,et al.  A multiple-quadrature eight-node hexahedral finite element for large deformation elastoplastic analysis , 1998 .

[41]  Sheet Metal Forming Simulation for Aluminum Alloy Sheets , 2000 .

[42]  Frédéric Barlat,et al.  Prediction of six or eight ears in a drawn cup based on a new anisotropic yield function , 2006 .