Variational projection methods for gradient crystal plasticity using Lie algebras

Abstract

Summary A computational method is developed for evaluating the plastic strain gradient hardening term within a crystal plasticity formulation. While such gradient terms reproduce the size effects exhibited in experiments, incorporating derivatives of the plastic strain yields a nonlocal constitutive model. Rather than applying mixed methods, we propose an alternative method whereby the plastic deformation gradient is variationally projected from the elemental integration points onto a smoothed nodal field. Crucially, the projection utilizes the mapping between Lie groups and algebras in order to preserve essential physical properties, such as orthogonality of the plastic rotation tensor. Following the projection, the plastic strain field is directly differentiated to yield the Nye tensor. Additionally, an augmentation scheme is introduced within the global Newton iteration loop such that the computed Nye tensor field is fed back into the stress update procedure. Effectively, this method results in a fully implicit evolution of the constitutive model within a traditional displacement‐based formulation. An elemental projection method with explicit time integration of the plastic rotation tensor is compared as a reference. A series of numerical tests are performed for several element types in order to assess the robustness of the method, with emphasis placed upon polycrystalline domains and multi‐axis loading. Copyright © 2016 John Wiley & Sons, Ltd.

FY17 Status Report on the Micromechanical Finite Element Modeling of Creep Fracture of Grade 91 Steel

“…However, in the traditional finite element implementation of crystal plasticity, the elastic deformation gradient is treated as a local quantity that is calculated solely at the integration points inside elements of the mesh. Therefore, two techniques for projecting the elastic rotation tensor and evaluating its derivative were developed and compared through as a series of benchmarks [37]. As a candidate crystal plasticity model, the Mechanical Threshold Strength (MTS) model that is implemented in WARP3D was chosen.…”

“…As an alternative, a staggered procedure was developed for simultaneously solving for the stress tensor, hardening variables, and kinematic variables [37]. Observe that the calculation of the Nye tensor…”

A configurational force for adaptive re-meshing of gradient-enhanced poromechanics problems with history-dependent variables

“…In this work, we employ a property-preserving Lie algebra mapping to extrapolate scalar and tensorial state variables. While there are previous works that employ similar strategy for the recovery of stress [72] and the internal variables in finite strain regime [73,74], our works actually find that the Lie algebra mapping can also be crucial in the small-deformation regime, especially when the residuals of the equilibrium equations are highly sensitive to perturbation of internal variables (e.g., strain softening). We compute the spectral decomposition of the tensor and obtain the logarithms of the rotation and stretch component separately.…”

Combined crystal plasticity and grain boundary modeling of creep in ferritic-martensitic steels: I. Theory and implementation

“…The CP model implementation uses an implicit backward Euler time integration scheme and solves the resulting local nonlinear equations using the Newton-Raphson method. The details of its implementation are given elsewhere [85,86,97,101]. The key equations are summarized in the appendix C, with the inclusion of the diffusional creep contribution.…”