A primal–dual interior point method to implicitly update Gurson–Tvergaard–Needleman model

Abstract

This study proposes an implicit algorithm applying the primal–dual interior point method (PDIP method) to stabilize the stress update when using a class of the Gurson–Tvergaard–Needleman model (GTN model). The GTN model is widely used to realize the change in void volume fraction that governs ductile fracture in metals, but numerical instabilities arise due to shrinkage of the yield surface and the accelerated void growth. In fact, such shrinkage can lead to misjudgment of yield conditions when using conventional return mapping algorithms, since trial elastic stresses are computed assuming zero incremental plastic strain. In addition, the change in void volume fraction is often approximated in bilinear form to represent the acceleration of void growth, but should be smooth to apply nonlinear solution methods such as the Newton’s method. To avoid such inconvenience in the implicit stress update for the GTN model and ensure numerical stability, we propose an algorithm that replaces the constitutive equations with inequality constraints with an equivalent constrained optimization problem by applying the PDIP method. After verifying the numerical accuracy and convergence of the proposed implicit algorithm using iso-error maps, we demonstrate its capability through several numerical examples that cannot be solved by the conventional return mapping algorithm or the PDIP method applied only to the inequality constraint corresponding to the yield condition.