Improved explicit integration in plasticity

Abstract

SUMMARYThe article presents a simple but efficient numerical scheme for the integration of non-linear constitutive equations, in which the principal reason for the inaccuracy of the classical explicit schemes, for example forward-Euler scheme, is effectively eliminated. In the newly developed explicit scheme, where there is no need for iteration, the implementation simplicity of the forward-Euler scheme and accuracy of the approach, which is using the backward-Euler scheme to integrate the constitutive equations, are successfully combined. Computational performance of the proposed next increment corrects error (NICE) integration scheme, particularly regarding the accuracy and the CPU time consumption, is first analysed on a case of complex loading of a material point. When comparing it to the forward-Euler, backward-Euler, trapezoidal and midpoint integration schemes, it turns out that because of its capability of a fast and relatively accurate integration of the constitutive equations, the NICE scheme is very convenient for the integration of constitutive models, where a direct solution technique is used to solve a boundary value problem. Although the deduction of the new integration scheme is general, its implementation for shell applications needs particular care. Namely, in order to satisfy the zero normal stress condition during the whole integration, a through-thickness strain increment has to be adequately chosen in each integration step. The NICE scheme, which was also implemented into ABAQUS/Explicit via User Material Subroutine (VUMAT) interface platform, has been additionally compared with the ABAQUS/Explicit default integration scheme (backward-Euler) and forward-Euler scheme. Two loading case-studies, namely the bending of a square plate and the stretching of a specimen including the onset of necking, are considered with two constitutive models-the von Mises and GTN material model being adopted. Generally, the NICE scheme has demonstrated to be advantageous in cases, where reasonable accuracy and very fast integration of the constitutive model is demanded, which is mostly the case in engineering computations with a direct solution method, for example explicit dynamics and metal forming process simulations.

Algorithmic issues for three-invariant hyperplastic Critical State models

“…The popularity of the fully implicit BE approach over explicit schemes (for example [43]) is due to its relatively high accuracy for a given numerical effort, particularly when large strain increments are applied [5,42]. Working with elastic strains as the primary unknown, we can express the return mapping as…”

NURBS plasticity: Yield surface representation and implicit stress integration for isotropic inelasticity

“…Nayak and Zienkiewicz [24] were the first to formulate a general explicit stress integration procedure for use within the finite element method and several models have been subsequently implemented using their method. However, the major drawback of these methods is that they do not enforce the consistency condition at the end of the stress-strain path [37].…”

Reuleaux plasticity: Analytical backward Euler stress integration and consistent tangent

“…The BE scheme gives a fully implicit approximation. The popularity of this approach over explicit schemes (for example [15]) is due to its high level of accuracy for a given numerical effort, particularly when large strain increments are applied [13]. The associated perfect plasticity model presented here may be thought of as a simple hybrid sitting between the Drucker-Prager (D-P) and Mohr-Coulomb (M-C) formulations.…”