A large strain plasticity model for implicit finite element analyses

Abstract

Annual:11-1-89 to 12-31-90 14.The theoretical basis and numerical implementation of a plasticity model suitable for finite strains and rotations are described. The constitutive equations governing J2 flow the0!y' are formulated using strainsstresses and their rates defined on the unrotated frame of reference. Unhke models based on the classical J aumann (or corotational) stress rate, the present model predicts physically acceptable responses for homogeneous deformations of exceedingly large magnitude. The associated numerical algorithms accommodate the large strain increments that arise in finite-element formulations employing an implicit solution of the global equilibrium equations. The resulting computational framework divorces the finite rotation effects on . Consequently, all of the numerical refinements developed previously for small-strain plasticity (radial return with subincrementation, plane stress modifications, kinematic hardening, consistent tangent operators) are utilized without modificatIOn. Details of the numerical algorithms are provided including the necessary transformation matrices and additional techniques re9.uired for finite deformations in plane stress. Several numerical examples are presented to illustrate the realIstic responses predicted by the model and the robustness of the numerical procedures. ABSTRACTThe theoretical basis and numerical implementation of a plasticity model suitable for finite strains and rotations are described. The constitutive equations governing J2 flow theory are formulated using strains-stresses and their rates defined on the unrotated frame of reference. Unlike models based on the classical Jaumann (or corotational) stress rate, the present model predicts physically acceptable responses for homogeneous deformations of exceedingly large magnitude. The associated numerical algorithms accommodate the large strain increments that arise in finite-element formulations employing an implicit solution of the global equilibrium equations. The resulting computational framework divorces the finite rotation effects on strain-stress rates from integration of the rates to update the material response over a load (time) step. Consequently, all of the numerical refinements developed previously for small-strain plasticity (radial return with subincrementation, plane stress modifications, kinematic hardening, consistent tangent operators) are utilized without modification. Details of the numerical algorithms are provided including the necessary transformation matrices and additional techniques required for finite deformations in plane stress. Several numerical examples are presented to illustrate the realistic responses predicted by the model and the robustness of the numerical procedures.-ii - ACKNOWLEDGMENTS

A transferability model for brittle fracture including constraint and ductile tearing effects: a probabilistic approach

“…The stress-update proceeds as follows (see [59] for full details): (1) using Rn + 1/2 rotate the spatial increment of the deformation tensor (D), evaluated from Bn+l/2 . fl.u e , to the unrotated configuration, d n + 1 / 2 = R~+1/2Dn+l/2Rn+1/2; (2) compute the unrotated Cauchy stress at n+ 1 (t n + 1) using a conventional small-strain, backward Euler procedure; and (3) compute the spatial Cauchy stress at n + 1 as (In+l = R n + 1 t n + 1 R;+1' The polar decompositions insure accuracy in the rotational operations independent of the displacement gradient magnitudes over n~n + 1.…”

Numerical modeling of ductile crack growth in 3-D using computational cell elements

“…The stress update proceeds as follows (see [37] for full details): (1) using Rn + 1/2 rotate the spatial increment of the deformation tensor (D), evaluated from B n + 1/2 . /J..u e , to the unrotated configuration, dn+1/2=R~+1/2Dn+1/2Rn+1/2; (2) compute the unrotated Cauchy stress at n+ 1 (t n + 1 ) using a conventional small-strain, backward Euler procedure; and (3) compute the spatial Cauchy stress at n + 1 as a n + 1 = R n + 1 t n + 1 R;+1' The polar decompositions insure accuracy in the rotational operations independent of the displacement gradient magnitudes over n ~ n + 1.…”

Line‐search methods in general return mapping algorithms with application to porous plasticity

“…A typical class of algorithms preserving incremental objectivity are rotationally neutralized algorithms that operate on the level of kinematics only, e.g. [2,[16][17][18][19]. They make use of a rotated configuration where an evolution equation for the proper orthogonal rotation tensor K at a generalized midpoint at time t n+ , ∈ [0, 1], is defined as follows:…”