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

被引:47
作者
Safaei, Mohsen [1 ]
Lee, Myoung-Gyu [2 ]
De Waele, Wim [3 ]
机构
[1] MIT, Impact & Crashworthiness Lab, Cambridge, MA 02139 USA
[2] Korea Univ, Dept Mat Sci & Engn, Seoul 136701, South Korea
[3] Univ Ghent, Dept Mech Construct & Prod, Fac Engn & Architecture, BE-9052 Zwijnaarde, Belgium
基金
新加坡国家研究基金会;
关键词
Stress integration schemes; Elastoplasticity; Explict; Implicit; Flow rule; UMAT VUMAT subroutine; FINITE-ELEMENT FORMULATION; YIELD FUNCTION; SHEET METALS; PREDICTION; STRAIN; ELASTOPLASTICITY;
D O I
10.1016/j.cma.2015.07.014
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
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. (C) 2015 Elsevier B.V. All rights reserved.
引用
收藏
页码:414 / 445
页数:32
相关论文
共 40 条
[1]   Plane stress yield function for aluminum alloy sheets - part 1: theory [J].
Barlat, F ;
Brem, JC ;
Yoon, JW ;
Chung, K ;
Dick, RE ;
Lege, DJ ;
Pourgoghrat, F ;
Choi, SH ;
Chu, E .
INTERNATIONAL JOURNAL OF PLASTICITY, 2003, 19 (09) :1297-1319
[2]  
Belytschko T., 2000, Nonlinear Finite Elements for Continua and Structures
[3]  
BISHOP JFW, 1951, PHILOS MAG, V42, P414
[4]   Stress integration method for a nonlinear kinematic/isotropic hardening model and its characterization based on polycrystal plasticity [J].
Cardoso, Rui P. R. ;
Yoon, Jeong Whan .
INTERNATIONAL JOURNAL OF PLASTICITY, 2009, 25 (09) :1684-1710
[5]   Modeling the Bauschinger effect for sheet metals, part I: theory [J].
Chun, BK ;
Jinn, JT ;
Lee, JK .
INTERNATIONAL JOURNAL OF PLASTICITY, 2002, 18 (5-6) :571-595
[6]   Finite element simulation of sheet forming based on a planar anisotropic strain-rate potential [J].
Chung, K ;
Lee, SY ;
Barlat, F ;
Keum, YT ;
Park, JM .
INTERNATIONAL JOURNAL OF PLASTICITY, 1996, 12 (01) :93-115
[7]  
CRISFIELD M.A., 2000, Non-linear Finite Element Analysis of Solids and Structures
[8]   A finite element formulation based on non-associated plasticity for sheet metal forming [J].
Cvitanic, Vedrana ;
Vlak, Frane ;
Lozina, Zeljan .
INTERNATIONAL JOURNAL OF PLASTICITY, 2008, 24 (04) :646-687
[9]  
Dunne F, 2005, INTRO COMPUTATIONAL
[10]  
Halilovic M., 2011, ENG COMPUT, P1