4.3.9 Deformation plasticity

Product: ABAQUS/Standard  

Deformation plasticity theory is provided in ABAQUS/Standard to allow fully plastic analysis of ductile metals, usually under small-displacement conditions, for fracture mechanics applications. The model is based on the Ramberg-Osgood relationship. In this section the detailed constitutive model is defined. The procedure for obtaining fully plastic solutions generally consists of incremental loading until the response is fully plastic.

The model is termed deformation plasticity because the stress is defined by the total mechanical strain with no history dependence. There is no “unloading” criterion (to allow recovery of the initial elastic stiffness immediately after a strain reversal), so that the model is only useful as a plasticity model in cases of continuous flow. It is, in fact, a nonlinear elastic model; but at a limit state when all of a specimen or structure is responding plastically, this model is a useful equivalent representation of the plastic response because it has such a simple form.

One-dimensional model

The basic one-dimensional model is

where is the stress, is the mechanical strain, E is Young's modulus (defined as the slope of the stress-strain curve at zero stress), is the “yield” offset (in the sense that, when , ), and n is the hardening exponent for the “plastic” (nonlinear) term: .

The material behavior described by this model is nonlinear at all stress levels, but for commonly used values of the hardening exponent ( or more) the nonlinearity becomes significant only at stress magnitudes approaching or exceeding .

Multiaxial generalization

A linear “elastic” relation is used to generalize the first term of Equation 4.3.9–1; the nonlinear term is generalized to multiaxial stress states through the use of the Mises stress potential and associated flow law, giving the multiaxial model

where

is the strain tensor,

is the stress tensor,

is the equivalent hydrostatic stress,

is the Mises equivalent stress,

is the stress deviator, and

is Poisson's ratio.

The linear part of the behavior may be compressible or incompressible depending on the value of Poisson's ratio, but the nonlinear part of the behavior is incompressible (because the flow is normal to the Mises stress potential). Since the model will generally be used for cases when the deformation is dominated by plastic flow, the use of selectively reduced integration elements or “hybrid” (mixed formulation) elements is recommended with this material model (except in plane stress).

Strain energy density

The model is often used to obtain fully plastic solutions for fracture mechanics when the J-integral is needed. For evaluation of the J-integral the strain energy density is required. This is

From Equation 4.3.9–2 this may be obtained as

Stress solution

During an analysis at each integration point the latest estimate of the kinematic solution is provided to the constitutive routines, which must provide the corresponding stress tensor calculated for the material model being used. Since this material model is nonlinear, we solve for the stresses by using the methods described below.

The uniaxial case (the only nonzero stress component is one direct stress)

In this case

where now .

We solve Equation 4.3.9–3 for using Newton's method. Writing as the correction to , the Newton equations for Equation 4.3.9–3 are

As an initial guess we use if and (with the sign chosen as the same sign as ) if .

In this case the material stiffness matrix is

All strain components defined by the kinematic solution

When all strain components are defined kinematically (that is, all cases except uniaxial and plane stress), projecting Equation 4.3.9–2 onto gives

a linear relationship for the volumetric behavior. Defining the deviatoric strain as

and using Equation 4.3.9–2, we find

Defining the equivalent deviatoric strain as

and using Equation 4.3.9–4, we obtain the scalar equation

This equation is solved for q using Newton's method:

As for the uniaxial case we will use the starting guesses:

and

Equation 4.3.9–6 can also be used to define

so that Equation 4.3.9–4 becomes

Thus, once q is known, is defined; and, hence, is known as

The material stiffness is defined as follows. From Equation 4.3.9–7 we have

and Equation 4.3.9–6 gives

and from Equation 4.3.9–5,

Using these results, Equation 4.3.9–8 becomes

Now

and

Combining these results, we obtain the material stiffness as

Plane stress

For this case we obtain a solution for assuming that the material is fully incompressible. We then use this solution as a starting guess for a Newton loop to find . Finally, we compute the corresponding material stiffness matrix.

Incompressible approximation:

In this case and , so that is known and Equation 4.3.9–6 can be solved as before for q. Equation 4.3.9–7 then defines , and the plane stress constraint requires that , so

and, hence, we obtain the initial estimate of the solution

and

Newton solution for the actual stresses:

Equation 4.3.9–2 defines the stresses, where, for this case,

and

Therefore, Equation 4.3.9–2 becomes

These equations are solved using Newton's method:

The material stiffness is directly available from Equation 4.3.9–9 as

Reference