An elasto-plastic self-consistent model for damaged polycrystalline materials: Theoretical formulation and numerical implementation

Elasto-plastic multiscale approaches are known to be suitable to model the mechanical behavior of metallic materials during forming processes. These approaches are classically adopted to explicitly link relevant microstructural effects to the macroscopic behavior. This paper presents a finite strain elasto-plastic self-consistent model for damaged polycrystalline aggregates and its implementation into ABAQUS/Standard finite element (FE) code. Material degradation is modeled by the introduction of a scalar damage variable at each crystallographic slip system for each individual grain. The single crystal plastic flow is described by both the classical and a regularized version of the Schmid criterion. To integrate the single crystal constitutive equations, two new numerical algorithms are developed (one for each plastic flow rule). Then, the proposed single crystal modeling is embedded into the self-consistent scheme to predict the mechanical behavior of elasto-plastic polycrystalline aggregates in the finite strain range. This strategy is implemented into ABAQUS/Standard FE code through a user-defined material (UMAT) subroutine. Special attention is paid to the satisfaction of the incremental objectivity and the efficiency of the convergence of the global resolution scheme, related to the computation of the consistent tangent modulus. The capability of the new constitutive modeling to capture the interaction between the damage evolution and the microstructural properties is highlighted through several simulations at both single crystal and polycrystalline scales. It appears from the numerical tests that the use of the classical Schmid criterion leads to a poor numerical convergence of the self-consistent scheme (due to the abrupt changes in the activity of the slip systems), which sometimes causes the computations to be prematurely stopped. By contrast, the use of the regularized version of the Schmid law allows a better convergence of the self-consistent approach, but induces an important increase in the computation time devoted to the integration of the single crystal constitutive equations (because of the high value of the power-law exponent used to regularize the Schmid yield function). To avoid these difficulties, a numerical strategy is built to combine the benefits of the two approaches: the classical Schmid criterion is used to integrate the single crystal constitutive equations, while its regularized version is used to compute the microscopic tange