Evaluation of stress integration algorithms for elastic–plastic constitutive models based on associated and non-associated flow rules

Abstract This paper presents in-depth analyses of four stress integration algorithms for finite deformation elastic–plastic constitutive relations. In particular, the four integration schemes were developed for both non-Associated Flow Rule (non-AFR) and Associated Flow Rule (AFR) in the continuum level of plasticity theory. The four integration schemes are (1) fully implicit backward Euler, (2) semi explicit convex cutting plane, (3) fully explicit classical forward Euler and (4) fully explicit forward Euler named Next Increment Corrects Error (NICE-1), which were implemented into the user material subroutines of finite element software ABAQUS. Analysis on numerical accuracy was carried out by using uniaxial tensile simulations at various orientations and time increments. The same analysis was also conducted for uniaxial tension/compression tests to investigate the effect of time increment on the calculated hardening curve. Finally, cylindrical cup deep drawing (manufacturing process which stamps a cylindrical cup) simulations were performed to compare computation time and accuracy for both AFR and non-AFR schemes for the evaluation of more realistic forming application. From the systematic comparative analyses, the following conclusions are made; (1) Caution should be exercised when finding an optimum time increment for explicit integration schemes. (2) The computation time and accuracy for both AFR and non-AFR are not much different when identical integration scheme is used. (3) For the simulation of large scale sheet metal forming applications, the explicit type stress integration algorithm can be a more practical choice considering that it does not deteriorate the computational accuracy and efficiency compared to the fully implicit algorithm.

[1]  Wim De Waele,et al.  Study on the definition of equivalent plastic strain under non-associated flow rule for finite element formulation , 2014 .

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

[3]  Marko Vrh,et al.  NICE - An explicit numerical scheme for efficient integration of nonlinear constitutive equations , 2009, Math. Comput. Simul..

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

[5]  S. Y. Lee,et al.  Finite element simulation of sheet forming based on a planar anisotropic strain-rate potential , 1996 .

[6]  R. Hill,et al.  XLVI. A theory of the plastic distortion of a polycrystalline aggregate under combined stresses. , 1951 .

[7]  Kwansoo Chung,et al.  Spring-back evaluation of automotive sheets based on isotropic-kinematic hardening laws and non-quadratic anisotropic yield functions: Part II: characterization of material properties , 2005 .

[8]  Wilfred D. Iwan,et al.  COMPUTATIONAL ASPECTS OF STRAIN SPACE PLASTICITY , 1983 .

[9]  Z. Marciniak,et al.  The mechanics of sheet metal forming , 1992 .

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

[11]  E. P. Popov,et al.  Accuracy and stability of integration algorithms for elastoplastic constitutive relations , 1985 .

[12]  R. Hill A theory of the yielding and plastic flow of anisotropic metals , 1948, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[13]  Cv Clemens Verhoosel,et al.  Non-Linear Finite Element Analysis of Solids and Structures , 1991 .

[14]  Željan Lozina,et al.  A finite element formulation based on non-associated plasticity for sheet metal forming , 2008 .

[15]  Earing Prediction in Cup Drawing Based on Non-Associated Flow Rule , 2007 .

[16]  J. C. Simo,et al.  A return mapping algorithm for plane stress elastoplasticity , 1986 .

[17]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[18]  Jeong Whan Yoon,et al.  A pressure-sensitive yield criterion under a non-associated flow rule for sheet metal forming , 2004 .

[19]  Jeong Whan Yoon,et al.  Stress integration method for a nonlinear kinematic/isotropic hardening model and its characterization based on polycrystal plasticity , 2009 .

[20]  Boris Štok,et al.  Improved explicit integration in plasticity , 2010 .

[21]  Wim De Waele,et al.  Evaluation of anisotropic constitutive models: Mixed anisotropic hardening and non-associated flow r , 2013 .

[22]  Jeong Whan Yoon,et al.  Anisotropic hardening and non-associated flow in proportional loading of sheet metals , 2009 .

[23]  Wim De Waele,et al.  An evolutionary anisotropic model for sheet metals based on non-associated flow rule approach , 2014 .

[24]  H. W. Swift Plastic instability under plane stress , 1952 .

[25]  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 .

[26]  J. C. Simo,et al.  Consistent tangent operators for rate-independent elastoplasticity☆ , 1985 .

[27]  Marko Vrh,et al.  NICEh: a higher-order explicit numerical scheme for integration of constitutive models in plasticity , 2011, Engineering with Computers.

[28]  Jeong Whan Yoon,et al.  A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming , 2010 .

[29]  R. D. Krieg,et al.  Accuracies of Numerical Solution Methods for the Elastic-Perfectly Plastic Model , 1977 .

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

[31]  R. F. Kulak,et al.  Accurate Numerical Solutions for Elastic-Plastic Models , 1979 .

[32]  Thomas B. Stoughton,et al.  A non-associated flow rule for sheet metal forming , 2002 .

[33]  Kwansoo Chung,et al.  Non-associated flow rule with symmetric stiffness modulus for isotropic-kinematic hardening and its application for earing in circular cup drawing , 2012 .

[34]  N. Petrinic,et al.  Introduction to computational plasticity , 2005 .

[35]  J. K. Lee,et al.  Modeling the Bauschinger effect for sheet metals, part I: theory , 2002 .

[36]  Michael Ortiz,et al.  An analysis of a new class of integration algorithms for elastoplastic constitutive relations , 1986 .

[37]  S. Thuillier,et al.  A model of one-surface cyclic plasticity and its application to springback prediction , 2011 .

[38]  Mohsen Safaei,et al.  Constitutive modelling of anisotropic sheet metals based on a non-associated flow rule , 2013 .