Levenberg–Marquardt vs Powell’s dogleg method for Gurson–Tvergaard–Needleman plasticity model

Abstract The GTN continuous damage model is very popular in academia and industry for structural integrity assessment and ductile fracture simulation. Following Aravas’ influential 1987 paper, Newton’s method has been used widely to solve the GTN equations. However, if the starting point is far from the solution, then Newton’s method can fail to converge. Hybrid methods are preferred in such cases. In this work we translate the GTN equations into a non-linear minimization problem and then apply the Levenberg–Marquardt and Powell’s ‘dogleg’ hybrid methods to solve it. The methods are tested for accuracy and robustness on two simple single finite element models and two 3D models with complex deformation paths. In total nearly 137,000 different GTN problems were solved. We show that the Levenberg–Marquardt method is more robust than Powell’s method. Our results are verified against the Abaqus’ own solver. The superior accuracy of the Levenberg–Marquardt method allows for larger time increments in implicit time integration schemes.

[1]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[2]  B. K. Dutta,et al.  A phenomenological form of the q2 parameter in the Gurson model , 2008 .

[3]  N. Aravas.  On the numerical integration of a class of pressure-dependent plasticity models , 1987 .

[4]  Philip E. Gill,et al.  Practical optimization , 1981 .

[5]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[6]  J. Devaux,et al.  A Methodology for Ductile Fracture Analysis Based on Damage Mechanics: An Illustration of a Local Approach of Fracture , 1988 .

[7]  A. Needleman,et al.  Void Nucleation Effects in Biaxially Stretched Sheets , 1980 .

[8]  Miroslav Zivkovic,et al.  Implicit stress integration procedure for small and large strains of the Gurson material model , 2002 .

[9]  Yonggang Huang,et al.  Effects of superimposed hydrostatic pressure on sheet metal formability , 2009 .

[10]  Anne Habraken,et al.  On the numerical integration of an advanced Gurson model , 2011 .

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

[12]  Zhiliang Zhang,et al.  On the accuracies of numerical integration algorithms for Gurson-based pressure-dependent elastoplastic constitutive models , 1995 .

[13]  L. E. Scales,et al.  Introduction to Non-Linear Optimization , 1985 .

[14]  Jorge J. Moré,et al.  The Levenberg-Marquardt algo-rithm: Implementation and theory , 1977 .

[15]  Mjd Powell,et al.  A Fortran subroutine for solving systems of non-linear algebraic equations , 1968 .

[16]  E. Roos,et al.  Investigation of failure behavior of ferritic–austenitic type of dissimilar steel welded joints , 2011 .

[17]  Joonmo Choung,et al.  Impact test simulations of stiffened plates using the micromechanical porous plasticity model , 2010 .

[18]  A. Needleman,et al.  Analysis of the cup-cone fracture in a round tensile bar , 1984 .

[19]  D. W. Beardsmore,et al.  An Implementation of the Gurson-Tvergaard-Needleman Plasticity Model for ABAQUS Standard Using a Trust Region Method , 2006 .

[20]  A. Saxena,et al.  Nonlinear Fracture Mechanics: Volume II Elastic-Plastic Fracture , 1988 .

[21]  P. Manach,et al.  Aluminum alloy damage evolution for different strain paths – Application to hemming process , 2009 .