Abstract:Numerically robust algorithmic formulations suitable for rate-independent crystal plasticity are presented. They cover classic local models as well as gradient-enhanced theories in which the gradients of the plastic slips are incorporated by means of the micromorphic approach. The elaborated algorithmic formulations rely on the underlying variational structure of (associative) crystal plasticity. To be more precise and in line with so-called variational constitutive updates or incremental energy minimization principles, an incrementally defined energy derived from the underlying time-continuous constitutive model represents the starting point of the novel numerically robust algorithmic formulations. This incrementally defined potential allows to compute all variables jointly as minimizers of this energy. While such discrete variational constitutive updates are not new in general, they are considered here in order to employ powerful techniques from non-linear constrained optimization theory in order to compute robustly the aforementioned minimizers. The analyzed prototype models are based on (1) nonlinear complementarity problem (NCP) functions as well as on (2) the augmented Lagrangian formulation. Numerical experiments show the numerical robustness of the resulting algorithmic formulations. Furthermore, it is shown that the novel algorithmic ideas can also be integrated into classic, non-variational, return-mapping schemes.

What problem does this paper attempt to address?

The paper primarily aims to address the robustness and non-uniqueness issues of numerical algorithms in crystal plasticity theory, particularly for rate-independent crystal plasticity theory. Specifically, the paper addresses the following key issues: 1. **Non-uniqueness Issue**: In rate-independent local crystal plasticity theory, the active state calculation of slip systems is not unique, which can lead to severe numerical problems. This is because the spatial dimension of the stress deviator tensor (5 dimensions) is smaller than the number of potentially active slip systems (for example, for face-centered cubic lattice structures, there are usually 12 slip systems). This non-uniqueness makes it very difficult to determine which slip systems are active. 2. **Proposing Robust Algorithms**: To solve the above problem, the authors propose two robust algorithms based on nonlinear constrained optimization theory. These two methods are the Nonlinear Complementarity Problem (NCP) function-based method and the Augmented Lagrangian framework-based method. These methods allow the simultaneous calculation of all relevant variables as minimizers of an energy functional, thereby ensuring the robustness of the algorithms. 3. **Extension to Gradient-Enhanced Models**: In addition to traditional local models, the paper also considers gradient-enhanced crystal plasticity theory, which incorporates the gradient of slip systems into the model through micromorphic approaches. By introducing incrementally defined energy, these non-local models can also be incorporated into the aforementioned variational framework and solved using the same algorithms. 4. **Algorithm Improvements**: The paper also proposes a new scaling method for the NCP function to improve the numerical performance of the algorithm. Additionally, by using the regularized Fischer-Burmeister NCP function, the local non-uniqueness problem in classical rate-independent crystal plasticity theory can be eliminated. In summary, the main contribution of this paper is the proposal of a series of robust algorithms to handle the non-uniqueness and computational challenges in crystal plasticity theory. These methods are not only applicable to traditional local models but can also be applied to more complex gradient-enhanced models.

Rate-Independent Gradient Crystal Plasticity Theory -- Robust Algorithmic Formulations based on Incremental Energy Minimization

On the Numerical Implementation of Non-Smooth Domains. Part I: Rate-dependent Formulations

A multiplicative finite strain crystal plasticity formulation based on additive elastic corrector rates: Theory and numerical implementation

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Quantized plastic deformation

A canonical rate-independent model of geometrically linear isotropic gradient plasticity with isotropic hardening and plastic spin accounting for the Burgers vector

On the Selection of Active Slip Systems in Crystal Plasticity

Vectorized MATLAB Implementation of the Incremental Minimization Principle for Rate-Independent Dissipative Solids Using FEM: A Constitutive Model of Shape Memory Alloys

An Effective Semi-Implicit Integration Scheme for Rate Dependent Crystal Plasticity Using Explicit Finite Element Codes

Finite Element Implementation of A Generalised Non-Local Rate-Dependent Crystallographic Formulation for Finite Strains

Modelling of Micro-Plasticity Evolution in Crystalline Materials

A trust region algorithm for finite-strain plasticity with strongly coupled hardening

A Micro‐mechanical Damage Model Based on Gradient Plasticity: Algorithms and Applications

Dissipation-consistent modelling and classification of extended plasticity formulations

On boundary conditions and plastic strain-gradient discontinuity in lower-order gradient plasticity

Calibration and surrogate model-based sensitivity analysis of crystal plasticity finite element models

Elastic-gap free strain gradient crystal plasticity model that effectively account for plastic slip gradient and grain boundary dissipation

Efficient numerical strategies for an implicit volume fraction transfer scheme for single crystal plasticity including twinning and secondary plasticity on the example of magnesium

A physics-based crystal plasticity model for the prediction of the dislocation densities in micropillar compression

Mechanism-based strain gradient crystal plasticity—I. Theory

Identifying material parameters in crystal plasticity by Bayesian optimization