Chapter 3

Constitutive Models $–$ Relations between Stress and Strain

3.7 Small Strain, Rate Independent Plasticity: Metals loaded beyond yield

For many design calculations, the elastic constitutive equations outlined in Section 3.1 are sufficient, since large plastic strains are by and large undesirable and will lead to failure.  There are some applications, however, where it is of interest to predict the behavior of solids subjected to large loads, sufficient to cause permanent plastic strains.  Examples include:

Modeling metal forming, machining or other manufacturing processes

Designing crash resistant vehicles

Plastic design of structures

Plasticity theory was developed to predict the behavior of metals under loads exceeding the plastic range, but the general framework of plasticity theory has since been adapted to other materials, including polymers and some types of soil (clay).  Some concepts from metal plasticity are also used in modeling concrete and other brittle materials such as polycrystalline ceramics.

3.7.1 Features of the inelastic response of metals.

We begin by reviewing the results of a typical tension/compression test on an annealed, ductile, polycrystalline metal specimen (e.g. copper or Al).  Assume that the test is conducted at moderate temperature (less than $½$ the melting point of the solid $–$ e.g. room temperature) and at modest strains (less than 10%), at modest strain rates ( ${10}^{1}-{10}^{-2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\text{s}}^{\text{-1}}$ ).

The results of such a test are

For modest stresses (and strains) the solid responds elastically.  This means the stress is proportional to the strain, and the deformation is reversible.

If the stress exceeds a critical magnitude, the stress$—$strain curve ceases to be linear.  It is often difficult to identify the critical stress accurately, because the stress strain curve starts to curve rather gradually.

If the critical stress is exceeded, the specimen is permanently changed in length on unloading.

If the stress is removed from the specimen during a test, the stress$—$strain curve during unloading has a slope equal to that of the elastic part of the stress$—$strain curve.  If the specimen is re-loaded, it will initially follow the same curve, until the stress approaches its maximum value during prior loading.  At this point, the stress$—$strain curve once again ceases to be linear, and the specimen is permanently deformed further.

If the test is interrupted and the specimen is held at constant strain for a period of time, the stress will relax slowly.  If the straining is resumed, the specimen will behave as though the solid were unloaded elastically.  Similarly, if the specimen is subjected to a constant stress, it will generally continue to deform plastically, although the plastic strain increases very slowly.  This phenomenon is known as creep.’

If the specimen is deformed in compression, the stress$—$strain curve is a mirror image of the tensile stress$—$strain curve (of course, this is only true for modest strains.  For large strains, geometry changes will cause differences between the tension and compression tests).

If the specimen is first deformed in compression, then loaded in tension, it will generally start to deform plastically at a lower tensile stress than an annealed specimen.  This phenomenon is known as the Bauschinger effect.’

Material response to cyclic loading can be extremely complex.  One example is shown in the picture above $–$ in this case, the material hardens cyclically.  Other materials may soften.

The detailed shape of the plastic stress$—$strain curve depends on the rate of loading, and also on temperature.

We also need to characterize the multi-axial response of an inelastic solid.  This is a much more difficult experiment to do.  Some of the nicest experiments were done by G.I. Taylor and collaborators in the early part of the last century.  Their approach was to measure the response of thin-walled tubes under combined torsion, axial loading and hydrostatic pressure.

The main conclusions of these tests were

The shape of the uniaxial stress-strain curve is insensitive to hydrostatic pressure.  However, the ductility (strain to failure) can be increased by adding hydrostatic pressure, particularly under torsional loading.

Plastic strains are volume preserving, i.e. the plastic strain rate must satisfy  ${\stackrel{˙}{\epsilon }}_{kk}=0$

During plastic loading, the principal components of the plastic strain rate tensor are parallel to the components of stress acting on the solid.  This sounds obvious until you think about it…  To understand what this means, imagine that you take a cylindrical shaft and pull it until it starts to deform plastically.  Then, holding the axial stress fixed, apply a torque to the shaft.  Experiments show that the shaft will initially stretch, rather than rotate.  The plastic strain increment is proportional to the stress acting on the shaft, not the stress increment.  This is totally unlike elastic deformation.

Under multi-axial loading, most annealed polycrystalline solids obey the Levy-Mises flow rule, which relates the principal components of strain rate during plastic loading to the principal stresses as follows

$\frac{{\stackrel{˙}{\epsilon }}_{1}-{\stackrel{˙}{\epsilon }}_{2}}{{\sigma }_{1}-{\sigma }_{2}}=\frac{{\stackrel{˙}{\epsilon }}_{1}-{\stackrel{˙}{\epsilon }}_{3}}{{\sigma }_{1}-{\sigma }_{3}}=\frac{{\stackrel{˙}{\epsilon }}_{2}-{\stackrel{˙}{\epsilon }}_{3}}{{\sigma }_{2}-{\sigma }_{3}}$

In this section, we will outline the simplest plastic constitutive equations that capture the most important features of metal plasticity.  There are many different plastic constitutive laws, which are intended to be used in different applications.   There are two broad classes:

1.      Rate independent plasticity $–$ which is used to model metals deformed at low temperatures (less than half the material’s melting point) and modest strain rates (of order 0.01-10/s).  This is the focus of this section.

2.      Rate dependent plasticity, or viscoplasticity $–$ used to model high temperature creep (steady accumulation of plastic strain at contstant stress) and also to model metals deformed at high strain rates (100/s or greater), where flow strength is sensitive to deformation rate.    Viscoplasticity will be discussed in Section 3.3.

There are also various different models within these two broad categories.   The models generally differ in two respects (i) the yield criterion; (ii) the strain hardening law.  There is no completely general model that describes all the features that were just listed, so in any application, you will need to decide which aspect of material behavior is most important, and then choose a model that accurately characterizes this behavior.

Key ideas in modeling metal plasticity

Five key concepts form the basis of almost all classical theories of plasticity.  They are

1.      The decomposition of strain into elastic and plastic parts;

2.      Yield criteria, which predict whether the solid responds elastically or plastically;

3.      Strain hardening rules, which control the way in which resistance to plastic flow increases with plastic straining;

4.      The plastic flow rule, which determines the relationship between stress and plastic strain under multi-axial loading;

5.      The elastic unloading criterion, which models the irreversible behavior of the solid.

These concepts will be described in more detail in the sections below.

For simplicity, we will at this stage restrict attention to infinitesimal deformations.

Consequently, we adopt the infinitesimal strain tensor

${\epsilon }_{ij}=\frac{1}{2}\left(\frac{\partial {u}_{i}}{\partial {x}_{j}}+\frac{\partial {u}_{j}}{\partial {x}_{i}}\right)$

as our deformation measure.  We have no need to distinguish between the various stress measures and will use ${\sigma }_{ij}$ to denote stress.

It is also important to note that the plastic strains in a solid depend on the load history.  This means that the stress-strain laws are not just simple equations relating stress to strain.  Instead, plastic strain laws must either relate the strain rate in the solid to the stress and stress rate, or else specify the relationship between a small increment of plastic strain $d{\epsilon }_{ij}^{p}$ in terms of strain, stress and stress increment $d{\sigma }_{ij}$.  In addition, plasticity problems are almost always solved using the finite element method.  Consequently, numerical methods are used to integrate the plastic stress-strain equations.

3.7.2. Decomposition of strain into elastic and plastic parts

Experiments show that under uniaxial loading, the strain at a given stress has two parts: a small recoverable elastic strain, and a large, irreversible plastic strain, as shown in the picture.  In uniaxial tension, we would write

$\epsilon ={\epsilon }^{e}+{\epsilon }^{p}$

Experiments suggest that the reversible part is related to the stress through the usual linear elastic equations.   Plasticity theory is concerned with characterizing the irreversible part.

For multiaxial loading, we generalize this by decomposing a general strain increment $d{\epsilon }_{ij}$ into elastic and plastic parts, as

$d{\epsilon }_{ij}=d{\epsilon }_{ij}^{e}+d{\epsilon }_{ij}^{p}$

The elastic part of the strain is related to stress using the linear elastic equations (discussed in detail in 3.1)

${C}_{ijkl}d{\epsilon }_{kl}^{e}=d{\sigma }_{ij}$

3.7.3 Yield Criteria

The yield criterion is used to determine the critical stress required to cause permanent deformation in a material.  There are many different yield criteria $–$ here we will just list the simplest ones.  Let ${\sigma }_{ij}$ be the stress acting on a solid, and let ${\sigma }_{1},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma }_{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma }_{3}$ denote the principal values of stress.  In addition, let $Y$ denote the yield stress of the material in uniaxial tension.  Then, define

Von$—$Mises yield criterion

$f\left({\sigma }_{ij},{\overline{\epsilon }}^{p}\right)=\sqrt{\frac{1}{2}\left[{\left({\sigma }_{1}-{\sigma }_{2}\right)}^{2}+{\left({\sigma }_{1}-{\sigma }_{3}\right)}^{2}+{\left({\sigma }_{2}-{\sigma }_{3}\right)}^{2}\right]}-Y\left({\overline{\epsilon }}^{p}\right)=0$

Tresca yield criterion

$f\left({\sigma }_{ij},{\overline{\epsilon }}^{p}\right)=\mathrm{max}\left\{|{\sigma }_{1}-{\sigma }_{2}|,|{\sigma }_{1}-{\sigma }_{3}|,|{\sigma }_{2}-{\sigma }_{3}|\right\}-Y\left({\overline{\epsilon }}^{p}\right)=0$

In both cases, the criteria are defined so that the material deforms elastically for $f\left({\sigma }_{ij}\right)<0$, and plastically for $f\left({\sigma }_{ij}\right)=0$.  The yield stress $Y$ may increase during plastic straining, so we have shown that Y is a function of a measure of total plastic strain ${\overline{\epsilon }}^{p}$, to be defined in Section 3.2.5

An alternative form of Von$—$Mises criterion. For a general stress state, it is a nuisance having to compute the principal stresses in order to apply von Mises yield criterion.  Fortunately, the criterion can be expressed directly in terms of the stress tensor

$f\left({\sigma }_{ij},{\overline{\epsilon }}^{p}\right)={\sigma }_{e}-Y\left({\overline{\epsilon }}^{p}\right)$

where

${\sigma }_{e}=\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{ij}={\sigma }_{ij}-\frac{1}{3}{\sigma }_{kk}{\delta }_{ij}$

are the components of the von Mises effective stress’ and  deviatoric stress tensor’ respectively.

These yield criteria are based largely on the following experimental observations:

(1)  A hydrostatic stress (all principal stresses equal) will never cause yield, no matter how large the stress;

(2) Most polycrystalline metals are isotropic.  Since the yield criterion depends only on the magnitudes of the principal stresses, and not their directions, the yield criteria predict isotropic behavior.

Tests suggest that von Mises yield criterion provides a slightly better fit to experiment than Tresca, but the difference between them is very small.   Sometimes it simplifies calculations to use Tresca’s criterion instead of von Mises.

3.7.4 Graphical representation of the yield surface.

Any arbitrary stress state can be plotted in  principal stress space,’ with the three principal stresses as axes.

The Von$—$Mises yield criterion is plotted in this way in the picture to the right.  The yield criterion is a cylinder, radius $Y/\sqrt{3}$, with its axis parallel to the line

${\sigma }_{1}={\sigma }_{2}={\sigma }_{3}$

If the state of stress falls within the cylinder, the material is below yield and responds elastically.  If the state of stress lies on the surface of the cylinder, the material yields and deforms plastically.  If the plastic deformation causes the material to strain harden, the radius of the cylinder increases.  The stress state cannot lie outside the cylinder $–$ this would lead to an infinite plastic strain.

Because the yield criterion $f\left({\sigma }_{ij}\right)=0$ defines a surface in stress space, it is often referred to as a yield surface. The yield surface is often drawn as it would appear when viewed down the axis of the cylinder, as shown below. The Tresca yield criterion can also be plotted in this way.  It looks like a cylinder with a hexagonal cross section, as shown.

3.7.5. Strain hardening laws

Experiments show that if you plastically deform a solid, then unload it, and then try to re-load it so as to induce further plastic flow, its resistance to plastic flow will have increased.  This is known as strain hardening.

Obviously, we can model strain hardening by relating the size and shape of the yield surface to plastic strain in some appropriate way.  There are many ways to do this.  Here we describe the two simplest approaches.

Isotropic hardening

Rather obviously, the easiest way to model strain hardening is to make the yield surface increase in size, but remain the same shape, as a result of plastic straining.

This means we must devise some appropriate relationship between Y and the plastic strain.  To get a suitable scalar measure of plastic strain we define the accumulated plastic strain magnitude

${\overline{\epsilon }}^{p}=\int \sqrt{\frac{2}{3}d{\epsilon }_{ij}^{p}d{\epsilon }_{ij}^{p}}$

(the factor of 2/3 is introduced so that ${\overline{\epsilon }}^{p}={\epsilon }_{11}^{p}$ in a uniaxial tensile test in which the specimen is stretched parallel to the ${e}_{1}$ direction.  To see this, note that plastic strains do not change volume, so that $d{\epsilon }_{22}=d{\epsilon }_{33}=-d{\epsilon }_{11}/2$ and substitute into the formula.)

Then we make Y a function of ${\overline{\epsilon }}^{p}$.  People often use power laws or piecewise linear approximations in practice. A few of the more common forms of hardening functions are

Perfectly plastic solid:  $Y=\text{constant}$

Linear strain hardening solid: $Y\left({\overline{\epsilon }}^{p}\right)={Y}_{0}+h{\overline{\epsilon }}^{p}$

Power$—$law hardening material: $Y={Y}_{0}+h{\left({\overline{\epsilon }}^{p}\right)}^{1/m}$

In these formulas, ${Y}_{0}$, h and m are material properties.  These functions are illustrated in the figures below

 Perfectly plastic solid Linear strain hardening solid Power-law hardening solid

Kinematic hardening

An isotropic hardening law is generally not useful in situations where components are subjected to cyclic loading.  It does not account for the Bauschinger effect, and so predicts that after a few cycles the solid will just harden until it responds elastically.

To fix this, an alternative hardening law allows the yield surface to translate, without changing its shape.  The idea is illustrated graphically in the picture.  As you deform the material in tension, you drag the yield surface in the direction of increasing stress, thus modeling strain hardening.  This softens the material in compression, however.  So, this constitutive law can model cyclic plastic deformation.  The stress-strain curves for isotropic and kinematic hardening materials are contrasted in the figure below.

To account for the fact that the center of the yield locus is at a position ${\alpha }_{ij}$ in stress space, the Von-Mises yield criterion needs to be modified as follows

$f\left({\sigma }_{ij},{\alpha }_{ij}\right)=\sqrt{\frac{3}{2}\left({S}_{ij}-{\alpha }_{ij}\right)\left({S}_{ij}-{\alpha }_{ij}\right)}-Y=0$

Here, Y is now a constant, and hardening is modeled by the motion of the yield surface. To do so, we need to relate ${\alpha }_{ij}$ to the plastic strain history somehow.  There are many ways to do this, which can model subtle features of the plastic response of solids under cyclic and nonproportional loading. The simplest approach is to set

$d{\alpha }_{ij}=\frac{2}{3}cd{\epsilon }_{ij}^{p}$.

This hardening law predicts that the stress-plastic strain curve is a straight line with slope c.  This is known as linear kinematic hardening

A more sophisticated approach is to set

$d{\alpha }_{ij}=\frac{2}{3}cd{\epsilon }_{ij}^{p}-\gamma {\alpha }_{ij}d{\overline{\epsilon }}^{p}$

where c and $\gamma$ are material constants.  It’s not so easy to visualize what this does $–$ it turns out that that this relation can model cyclic creep $–$ the tendency of a material to accumulate strain in the direction of mean stress under cyclic loading, as illustrated in the figure on the right.  It is known as the Armstrong-Frederick hardening law.

There are many other kinematic type hardening laws.  New ones are still being developed.

3.7.6 The plastic flow law

To complete the plastic stress-strain relations, we need a way to predict the plastic strains induced by stressing the material beyond the yield point.  Specifically, given

1.      The current stress ${\sigma }_{ij}$ applied to the material

2.      The current yield stress (characterized by $Y\left({\overline{\epsilon }}^{p}\right)$ for isotropic hardening, or ${\alpha }_{ij}$ for kinematic hardening)

3.      A small increase in stress $d{\sigma }_{ij}$ applied to the solid

we wish to determine the small change in plastic strain $d{\epsilon }_{ij}^{p}$.

The formulas are given below, for isotropic and kinematic hardening.  These are just fits to experiment (specifically, to the Levy-Mises flow rule).  The physical significance and reason for the structure of the equations will be discussed later.

The plastic strains are usually derived from the yield criterion f defined in 3.6.3, and so are slightly different for isotropic and kinematic hardening.  A material that has its plastic flow law derived from f is said to have an associated’ flow law $–$ i.e.the flow law is associated with f.

Isotropic Hardening (Von-Mises yield criterion)

$d{\epsilon }_{ij}^{p}=d{\overline{\epsilon }}^{p}\frac{\partial f}{\partial {\sigma }_{ij}}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}}{Y}$

where

$f\left({\sigma }_{ij},{\overline{\epsilon }}^{p}\right)=\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{ij}={\sigma }_{ij}-\frac{1}{3}{\sigma }_{kk}{\delta }_{ij}$

denotes the Von-Mises yield criterion, and $d{\overline{\epsilon }}^{p}$ is determined from the condition that the yield criterion must be satisfied at all times during plastic straining.  This shows that

$\begin{array}{l}f\left({\sigma }_{ij}+d{\sigma }_{ij},{\overline{\epsilon }}^{p}+d{\overline{\epsilon }}^{p}\right)=f\left({\sigma }_{ij},{\overline{\epsilon }}^{p}\right)+\frac{\partial f}{\partial {\sigma }_{ij}}d{\sigma }_{ij}+\frac{\partial f}{\partial {\overline{\epsilon }}^{p}}d{\overline{\epsilon }}^{p}=0\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⇒\frac{\partial f}{\partial {\sigma }_{ij}}d{\sigma }_{ij}-\frac{\partial Y}{\partial {\overline{\epsilon }}^{p}}d{\overline{\epsilon }}^{p}=0\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⇒d{\overline{\epsilon }}^{p}=\frac{1}{h}\frac{\partial f}{\partial {\sigma }_{ij}}d{\sigma }_{ij}=\frac{1}{h}\frac{3}{2}\frac{{S}_{ij}d{\sigma }_{ij}}{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}h=\frac{dY}{d{\overline{\epsilon }}^{p}}\end{array}$

Here, $h=\partial Y/\partial {\overline{\epsilon }}^{p}$ is the slope of the plastic stress-strain curve.   The algebra involved in differentiating f with respect to stress is outlined below.

Linear Kinematic Hardening (Von-Mises yield criterion)

$d{\epsilon }_{ij}^{p}=d{\overline{\epsilon }}^{p}\frac{\partial f}{\partial {\sigma }_{ij}}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{\left({S}_{ij}-{\alpha }_{ij}\right)}{Y}$

where the yield criterion is now

$f\left({\sigma }_{ij},{\alpha }_{ij}\right)=\sqrt{\frac{3}{2}\left({S}_{ij}-{\alpha }_{ij}\right)\left({S}_{ij}-{\alpha }_{ij}\right)}-Y\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{ij}={\sigma }_{ij}-\frac{1}{3}{\sigma }_{kk}{\delta }_{ij}$

and as before $d{\overline{\epsilon }}^{p}$ is determined from the condition that the yield criterion must be satisfied at all times during plastic straining.  This shows that

$f\left({\sigma }_{ij}+d{\sigma }_{ij},{\alpha }_{ij}+d{\alpha }_{ij}\right)=f\left({\sigma }_{ij},{\alpha }_{ij}\right)+\frac{\partial f}{\partial {\sigma }_{ij}}d{\sigma }_{ij}+\frac{\partial f}{\partial {\alpha }_{ij}}d{\alpha }_{ij}=0$

Recall that for linear kinematic hardening the hardening law is

$d{\alpha }_{ij}=\frac{2}{3}cd{\epsilon }_{ij}^{p}=cd{\overline{\epsilon }}^{p}\frac{\left({S}_{ij}-{\alpha }_{ij}\right)}{Y}$

Substituting into the Taylor expansion of the yield criterion and simplifying shows that

$d{\overline{\epsilon }}^{p}=\frac{1}{c}\frac{3}{2}\frac{\left({S}_{ij}-{\alpha }_{ij}\right)d{\sigma }_{ij}}{Y}$

Comparison of flow law formulas with the Levy-Mises flow rule

The Levy-Mises flow law (based on experimental observations) states that principal values of the plastic strain increment $d{\epsilon }_{1},d{\epsilon }_{2},d{\epsilon }_{3}$ induced by a stress increment are related to the principal stresses ${\sigma }_{1},{\sigma }_{2},{\sigma }_{3}$ by

$\frac{d{\epsilon }_{1}-d{\epsilon }_{2}}{{\sigma }_{1}-{\sigma }_{2}}=\frac{d{\epsilon }_{1}-d{\epsilon }_{3}}{{\sigma }_{1}-{\sigma }_{3}}=\frac{d{\epsilon }_{2}-d{\epsilon }_{3}}{{\sigma }_{2}-{\sigma }_{3}}$

It is straightforward to show that this observation is consistent with the predictions of the flow law formulas given in this section. To see this, suppose that the principal axes of stress are parallel to the $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ directions.  In this case the only nonzero components of deviatoric stress are

${S}_{11}={\sigma }_{1}-\left({\sigma }_{1}+{\sigma }_{2}+{\sigma }_{3}\right)/3\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{22}={\sigma }_{2}-\left({\sigma }_{1}+{\sigma }_{2}+{\sigma }_{3}\right)/3\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{33}={\sigma }_{1}-\left({\sigma }_{1}+{\sigma }_{2}+{\sigma }_{3}\right)/3$

The flow law

$d{\epsilon }_{ij}^{p}=d{\overline{\epsilon }}^{p}\frac{\partial f}{\partial {\sigma }_{ij}}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}}{Y}$

gives

$d{\epsilon }_{1}^{p}=d{\epsilon }_{11}^{p}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{11}}{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}d{\epsilon }_{2}^{p}=d{\epsilon }_{22}^{p}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{22}}{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}d{\epsilon }_{3}^{p}=d{\epsilon }_{33}^{p}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{33}}{Y}$

Thus, we see that

$\begin{array}{l}d{\epsilon }_{1}^{p}-d{\epsilon }_{2}^{p}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{11}-{S}_{22}}{Y}=\text{\hspace{0.17em}}\text{\hspace{0.17em}}d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{\sigma }_{1}-{\sigma }_{2}}{Y}\\ d{\epsilon }_{1}^{p}-d{\epsilon }_{3}^{p}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{11}-{S}_{33}}{Y}=\text{\hspace{0.17em}}\text{\hspace{0.17em}}d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{\sigma }_{1}-{\sigma }_{3}}{Y}\end{array}$

with similar expressions for other components.  Some trivial algebra then yields the Levy-Mises flow law.

Differentiating the yield criterion

Differentiating the yield criterion requires some sneaky index notation manipulations. Note that

$\frac{\partial f}{\partial {\sigma }_{ij}}=\frac{3}{2}\frac{1}{2}\frac{1}{\sqrt{\frac{3}{2}{S}_{kl}{S}_{kl}}}2{S}_{pq}\frac{\partial {S}_{pq}}{\partial {\sigma }_{ij}}\text{\hspace{0.17em}}=\frac{3}{2}\frac{1}{\sqrt{\frac{3}{2}{S}_{kl}{S}_{kl}}}{S}_{pq}\frac{\partial {S}_{pq}}{\partial {\sigma }_{ij}}$

Now, recall that

${S}_{ij}={\sigma }_{ij}-\frac{1}{3}{\sigma }_{kk}{\delta }_{ij}$

and further that

$\frac{\partial {\sigma }_{ij}}{\partial {\sigma }_{kl}}={\delta }_{ik}{\delta }_{jl}$

Hence

$\frac{\partial {S}_{pq}}{\partial {\sigma }_{ij}}={\delta }_{ip}{\delta }_{jq}-\frac{1}{3}{\delta }_{ik}{\delta }_{jk}{\delta }_{pq}$

and

${S}_{pq}\frac{\partial {S}_{pq}}{\partial {\sigma }_{ij}}={S}_{pq}\left({\delta }_{ip}{\delta }_{jq}-\frac{1}{3}{\delta }_{ik}{\delta }_{jk}{\delta }_{pq}\right)={S}_{ij}-{S}_{pp}{\delta }_{ij}$

However, observe that

${S}_{pp}={\sigma }_{pp}-\frac{1}{3}{\sigma }_{kk}{\delta }_{pp}=0$

so that

${S}_{pq}\frac{\partial {S}_{pq}}{\partial {\sigma }_{ij}}={S}_{ij}$

and finally

$\frac{\partial f}{\partial {\sigma }_{ij}}=\frac{3}{2}\frac{{S}_{ij}}{\sqrt{\frac{3}{2}{S}_{kl}{S}_{kl}}}=\frac{3}{2}\frac{{S}_{ij}}{Y}$

There is one final issue to consider.  Experiments show that plastic flow is irreversible, and always dissipates energy.  If the increment in stress $d{\sigma }_{ij}$ is tangent to the yield surface, or brings the stress below yield, as shown in the picture then there is no plastic strain.

For an isotropically hardening solid, this unloading condition may be expressed as

${S}_{ij}d{\sigma }_{ij}^{}<0$

For kinematic hardening,

$\left({S}_{ij}-{\alpha }_{ij}\right)d{\sigma }_{ij}^{}<0$

In both cases, the solid deforms elastically (no plastic strain) if the condition is satisfied.

3.7.8 Complete incremental stress-strain relations for a rate independent elastic-plastic solid

We conclude by summarizing the complete elastic-plastic stress strain relations for an isotropic  solid with Von-Mises yield surface.

Isotropically hardening elastic-plastic solid

The solid is characterized by its elastic constants $E,\nu$ and by the yield stress $Y\left({\overline{\epsilon }}^{p}\right)$ as a function of accumulated plastic strain ${\overline{\epsilon }}^{p}$ and its slope $h=\frac{dY}{d{\overline{\epsilon }}^{p}}$

In this case we have that

$d{\epsilon }_{ij}=d{\epsilon }_{ij}^{e}+d{\epsilon }_{ij}^{p}$

with

$d{\epsilon }_{ij}^{e}=\frac{1+\nu }{E}\left(d{\sigma }_{ij}-\frac{\nu }{1+\nu }d{\sigma }_{kk}{\delta }_{ij}\right)$

$d{\epsilon }_{ij}^{p}=\left\{\begin{array}{c}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)<0\\ \frac{1}{h}\frac{3}{2}\frac{〈{S}_{kl}d{\sigma }_{kl}〉}{Y}\frac{3}{2}\frac{{S}_{ij}}{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)=0\end{array}$

where $〈x〉=\left\{\begin{array}{c}x\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\ge 0\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\le 0\end{array}$

These may be combined to

$d{\epsilon }_{ij}^{}=\left\{\begin{array}{c}\frac{1+\nu }{E}\left(d{\sigma }_{ij}-\frac{\nu }{1+\nu }d{\sigma }_{kk}{\delta }_{ij}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)<0\\ \frac{1+\nu }{E}\left(d{\sigma }_{ij}-\frac{\nu }{1+\nu }d{\sigma }_{kk}{\delta }_{ij}\right)+\frac{1}{h}\frac{3}{2}\frac{〈{S}_{kl}d{\sigma }_{kl}〉}{Y}\frac{3}{2}\frac{{S}_{ij}}{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)=0\end{array}$

It is sometimes necessary to invert these expressions.  A straightforward but tedious series of index notation manipulations shows that

$d{\sigma }_{ij}=\left\{\begin{array}{c}\frac{E}{1+\nu }\left\{d{\epsilon }_{ij}+\frac{\nu }{1-2\nu }d{\epsilon }_{kk}{\delta }_{ij}\right\}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)<0\\ \frac{E}{1+\nu }\left\{d{\epsilon }_{ij}+\frac{\nu }{1-2\nu }d{\epsilon }_{kk}{\delta }_{ij}-\frac{3E}{3E+2\left(1+\nu \right)h}\frac{3}{2}\frac{〈{S}_{kl}d{\epsilon }_{kl}〉}{Y}\frac{3}{2}\frac{{S}_{ij}}{Y}\right\}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}{S}_{ij}{S}_{ij}}-Y\left({\overline{\epsilon }}^{p}\right)=0\end{array}\text{\hspace{0.17em}}$

This constitutive law is the most commonly used model of inelastic deformation.  It has the following properties:

It will correctly predict the conditions necessary to initiate yield under multiaxial loading

It will correctly predict the plastic strain rate under an arbitrary multiaxial stress state

It can model accurately any uniaxial stress$—$strain curve

It has the following limitations:

It is valid only for modest plastic strains (<10%)

It will not predict creep behavior or strain rate sensitivity

It will not predict plastic strains accurately if the principal axes of stress rotate significantly (more than about 30 degrees) during inelastic deformation

Linear Kinematically hardening solid

The solid is characterized by its elastic constants $E,\nu$ and by the initial yield stress $Y$ and the strain hardening rate c. Then,

$d{\epsilon }_{ij}=d{\epsilon }_{ij}^{e}+d{\epsilon }_{ij}^{p}$

with

$d{\epsilon }_{ij}^{e}=\frac{1+\nu }{E}\left(d{\sigma }_{ij}-\frac{\nu }{1+\nu }d{\sigma }_{kk}{\delta }_{ij}\right)$

$d{\epsilon }_{ij}^{p}=\left\{\begin{array}{c}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}\left({S}_{ij}-{\alpha }_{ij}\right)\left({S}_{ij}-{\alpha }_{ij}\right)}-Y<0\\ \frac{1}{c}\frac{3}{2}\frac{〈\left({S}_{kl}-{\alpha }_{kl}\right)d{\sigma }_{kl}〉}{Y}\frac{3}{2}\frac{\left({S}_{ij}-{\alpha }_{ij}\right)}{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sqrt{\frac{3}{2}\left({S}_{ij}-{\alpha }_{ij}\right)\left({S}_{ij}-{\alpha }_{ij}\right)}-Y=0\end{array}$

where $〈x〉=\left\{\begin{array}{c}x\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\ge 0\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\le 0\end{array}$

Finally, the evolution equation for ${\alpha }_{ij}$ is

$d{\alpha }_{ij}=\frac{3}{2}cd{\epsilon }_{ij}^{p}$

This constitutive equation is used primarily to model cyclic plastic deformation, or plastic flow under nonproportional loading (where principal axes of stress rotate significantly during plastic flow). It has the following limitations:

It is valid only for modest plastic strains (<10%)

It will not predict creep behavior or strain rate sensitivity

It does not predict the shape of the stress-strain curve accurately

3.7.9 Typical values for yield stress of polycrystalline metals

Unlike elastic constants, the plastic properties of metals are highly variable, and are also very sensitive to alloying composition and microstructure (which can be influenced by heat treatment and mechanical working).  Consequently, it is impossible to give accurate values for yield stresses or hardening rates for materials.  The table below (again, taken from `Engineering Materials,’ by M.F. Ashby and D.R.H. Jones, Pergamon Press) lists rough values for yield stresses of common materials $–$ these may provide a useful guide in preliminary calculations.  If you need accurate data you will have to measure the properties of the materials you plan to use yourself.

 Material Yield Stress ${\sigma }_{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}/\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}MN{m}^{-2}$ Material Yield Stress ${\sigma }_{Y}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}/\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}MN{m}^{-2}$ Tungsten Carbide 6000 Mild steel 220 Silicon Carbide 10 000 Copper 60 Tungsten 2000 Titanium 180 - 1320 Alumina 5000 Silica glass 7200 Titanium Carbide 4000 Aluminum & alloys 40-200 Silicon Nitride 8000 Polyimides 52 - 90 Nickel 70 Nylon 49 - 87 Iron 50 PMMA 60 - 110 Low alloy steels 500-1980 Polycarbonate 55 Stainless steel 286-500 PVC 45-48

3.7.10 Perspectives on plastic constitutive equations - The Principle of Maximum Plastic Resistance

The constitutive law outlined in the preceding section has an important property, known as the principle of maximum plastic resistance.

Statement of the principle: Let ${\sigma }_{ij}$ be a stress state which causes plastic deformation, let $d{\sigma }_{ij}$ be a small change in ${\sigma }_{ij}$ and let $d{\epsilon }_{ij}^{p}$ be the resulting strain increment.  Now, let ${\sigma }_{ij}^{\ast }$ be any other stress that can be imposed on the specimen that either does not reach yield, or else just satisfies the yield criterion, i.e. $\sqrt{3{S}_{ij}^{*}{S}_{ij}^{*}}\le Y$ with ${S}_{ij}^{*}={\sigma }_{ij}^{*}-{\sigma }_{kk}^{*}{\delta }_{ij}$.

Then

$\left({\sigma }_{ij}-{\sigma }_{ij}^{\ast }\right)\frac{d{\epsilon }_{ij}}{dt}\ge 0$

Interpretation: The Principle of Maximum Plastic Resistance is a mathematical statement of the following ideas:

(i) The Mises yield surface is convex

(ii) The plastic strain rate is normal to the yield surface.

It is best to illustrate these ideas graphically. In principal stress space, the product ${\sigma }_{ij}d{\epsilon }_{ij}^{p}$ is represented by the dot product of the stress and plastic strain rate vectors.  The statement

$\left({\sigma }_{ij}-{\sigma }_{ij}^{\ast }\right)d{\epsilon }_{ij}^{p}\ge 0$

is equivalent to the requirement that the angle between the vectors formed by ${\sigma }_{ij}-{\sigma }_{ij}^{\ast }$ and $d{\epsilon }_{ij}^{p}$ must be greater than 90 degrees for all stresses and strain rates.  This is only possible if the yield stress is convex and the strain rate is normal to the yield surface.

The Principle of Maximum Plastic Resistance is important because it is the basis for a number of very important theorems concerning plastic deformation in solids.  For example, it can be shown that the stress field in a material that obeys the Principle is always unique.  In addition, the principle leads to clever techniques to estimate collapse loads for elastic-plastic solids and structures.

Proof of the principle of maximum plastic resistance

Our goal is to prove that ${\sigma }_{ij}d{\epsilon }_{ij}^{p}\ge {\sigma }_{ij}^{*}d{\epsilon }_{ij}^{p}$.  The simplest way to do so is to show that ${\sigma }_{ij}d{\epsilon }_{ij}^{p}=Yd{\overline{\epsilon }}^{p}$, while ${\sigma }_{ij}^{*}d{\epsilon }_{ij}^{p}\le Yd{\overline{\epsilon }}^{p}$, where $d{\overline{\epsilon }}^{p}=\sqrt{2d{\epsilon }_{ij}^{p}d{\epsilon }_{ij}^{p}/3}$ is the plastic strain magnitude, and Y is the yield stress.  To this end:

1.      Recall the plastic flow rule  $d{\epsilon }_{ij}^{p}=d{\overline{\epsilon }}^{p}\frac{\partial f}{\partial {\sigma }_{ij}}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}}{Y}$

2.      Multiply both sides by stress

${\sigma }_{ij}d{\epsilon }_{ij}^{p}={\sigma }_{ij}d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}}{Y}=\left({S}_{ij}+{\sigma }_{kk}{\delta }_{ij}\right)d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}}{Y}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}{S}_{ij}}{Y}$

where we have noted that ${S}_{kk}=0$.

3.      Recall that ${S}_{ij}$ causes yield, and so must satisfy the yield condition $\sqrt{3{S}_{ij}{S}_{ij}/2}-Y=0$. This shows that ${\sigma }_{ij}d{\epsilon }_{ij}^{p}=Yd{\overline{\epsilon }}^{p}$

4.      Now consider

${\sigma }_{ij}^{*}d{\epsilon }_{ij}^{p}={\sigma }_{ij}^{*}d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}}{Y}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}^{*}{S}_{ij}}{Y}$

5.      Note that

$\frac{3}{2}\left({S}_{ij}-{S}_{ij}^{*}\right)\left({S}_{ij}-{S}_{ij}^{*}\right)\ge 0⇒\frac{3}{2}\left({S}_{ij}{S}_{ij}+{S}_{ij}^{*}{S}_{ij}^{*}-2{S}_{ij}{S}_{ij}^{*}\right)\ge 0$

6.      Now, recall again that ${S}_{ij}$ causes yield, while  ${S}_{ij}^{*}$ can be at or below yield.  The yield criterion therefore requires that $3{S}_{ij}{S}_{ij}/2={Y}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}3{S}_{ij}^{*}{S}_{ij}^{*}/2\le {Y}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}$. Substituting these inequalities into (5) shows that ${Y}^{2}\ge 3{S}_{ij}{S}_{ij}^{*}/2$.  Finally, this shows

${\sigma }_{ij}^{*}d{\epsilon }_{ij}^{p}=d{\overline{\epsilon }}^{p}\frac{3}{2}\frac{{S}_{ij}^{*}{S}_{ij}}{Y}\le Yd{\overline{\epsilon }}^{p}$

Thus ${\sigma }_{ij}d{\epsilon }_{ij}^{p}\ge {\sigma }_{ij}^{\ast }d{\epsilon }_{ij}^{p}$, proving the principle.

3.7.11 Perspectives on plastic constitutive equations - Drucker’s Postulate

Constitutive models of inelastic behavior are based largely on experimental observations of plastic flow in laboratory specimens.  Similar constitutive laws are used to describe very different materials, including metals, ceramics, glasses, soils and polymers.  The mechanisms of deformation in these materials are very different, so it is surprising that their response is similar.

One perspective on the structure of constitutive laws for inelastic solids was developed by Drucker in the 1950s.  Drucker introduced the idea of a stable plastic material, as follows. Consider a deformable solid, subjected to boundary tractions ${t}_{i}$, which induce some displacement field ${u}_{i}$.  Suppose that the tractions are increased to ${t}_{i}+\Delta {t}_{i}$, resulting in an additional displacement $\Delta {u}_{i}$. The material is said to be stable in the sense of Drucker if the work done by the tractions $\Delta {t}_{i}$ through the displacements $\Delta {u}_{i}$ is positive or zero for all $\Delta {t}_{i}$:

$\Delta W=\int \left\{\underset{A}{\int }\Delta {t}_{i}\frac{d\Delta {u}_{i}}{dt}\right\}dt\ge 0$

It can be shown that, for a plastic material to be stable in this sense, it must satisfy the following conditions:

The yield surface $f\left({\sigma }_{ij}\right)$ must be convex

The plastic strain rate must be normal to the yield surface $d{\epsilon }_{ij}^{p}=d{\overline{\epsilon }}^{p}\frac{\partial f}{\partial {\sigma }_{ij}}$

The rate of strain hardening must be positive or zero $\frac{dY}{d{\overline{\epsilon }}^{p}}\ge 0$

Furthermore, a material that is stable in the sense of Drucker must satisfy the principle of maximum plastic resistance.

This does not really explain why the constitutive law should have this structure, but materials that do not satisfy the Drucker stability criterion tend to be difficult to work with in calculations, so there is a strong incentive for choosing a constitutive law that meets the condition.  It is not surprising, then, that the people developing constitutive laws ended up with a form that satisfies Drucker stability.

3.7.12 Microscopic Perspectives on Plastic Flow in Metals

It is possible to obtain some insight into the structure of the constitutive laws for metals by considering the micromechanisms responsible for plastic flow.

Plastic flow in metals is caused by dislocation motion.  Dislocations are line defects in crystalline solids $–$ you can think of a dislocation as an extra plane of atoms inserted within a perfect crystal, as shown in the picture on the right.

When the crystal is subjected to stress, these defects move through the solid and rearrange the crystal lattice.  For example, if the model crystal shown above is subjected to a shear stress, the atoms rearrange so that the top part of the crystal is shifted to the right relative to the bottom part.

Because the crystal lattice is distorted near a dislocation, only a modest shear stress is required to drive the dislocation through the solid, causing permanent plastic deformation.

Experiments and atomistic simulations suggest that dislocation motion obeys Schmidt’s law: a dislocation moves through a crystal if the shear stress on its glide plane exceeds a critical magnitude ${\tau }_{c}$.

It can be shown that a material which deforms by dislocation glide, and which obeys Schmidt’s law, will satisfy the principle of maximum plastic resistance.  This, in turn, implies that the yield surface for the solid must be convex and the plastic strain rate must be normal to the yield surface.

The notion of a yield surface and convexity for a material which deforms by dislocation glide can be illustrated with a simple double-slip model.  Consider a single crystal, which contains two dislocation glide planes oriented at ${45}^{o}$ to the principal axes of stress as shown in the picture on the right.  Assume plane stress conditions, for simplicity.

The picture below shows side views of the two slip planes. As an exercise, you should verify that the shear stresses (tangential component of traction) acting on the two slip planes are

$\begin{array}{l}\tau =±\frac{1}{2}\left({\sigma }_{1}-{\sigma }_{2}\right)\\ \tau =±\frac{1}{2}{\sigma }_{2}\end{array}$

The solid reaches yield if $\tau ={\tau }_{c}$.  The resulting yield surface is sketched below: the planes of the yield surface have been color coded to indicate which slip plane is active.

Observe that the yield surface is convex: this is a consequence of Schmidt’s law.

Now, suppose that slip is activated on one of the glide planes.  Let ${t}_{i}$ denote the tangent to the slip plane, and let ${n}_{i}$ denote the normal. To compute the strain produced by slip on a single slip system, consider the deformation of an infinitesimal line element dx under shear $\gamma$.

The deformed line element is given by

$dw=dx+\gamma \left(dx\cdot n\right)t$

or, in index notation

$d{w}_{i}=\left({\delta }_{ij}+\gamma \text{\hspace{0.17em}}{t}_{i}{n}_{j}\right)d{x}_{j}$

${F}_{ij}={\delta }_{ij}+\gamma \text{\hspace{0.17em}}{t}_{i}{n}_{j}$

and the Lagrange strain tensor is

${E}_{ij}=\frac{1}{2}\left({F}_{ki}{F}_{jk}-{\delta }_{ij}\right)=\frac{1}{2}\gamma \left({t}_{i}{n}_{j}+{t}_{j}{n}_{i}\right)+\frac{1}{2}{\gamma }^{2}{n}_{i}{n}_{j}$

For small $\gamma$, we can approximate the Lagrange strain tensor by the infinitesimal strain tensor

${\epsilon }_{ij}=\frac{1}{2}\gamma \left({t}_{i}{n}_{j}+{t}_{j}{n}_{i}\right)$

Now, suppose that the stress satisfies

${\sigma }_{2}-{\sigma }_{1}={\tau }_{c}$

as marked on the yield locus shown below.  This activates slip as shown in the picture.  The normal and tangent to the appropriate slip plane are

$t=\left[\frac{-1}{\sqrt{2}},\frac{1}{\sqrt{2}},0\right]\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=\left[\frac{1}{\sqrt{2}},\frac{1}{\sqrt{2}},0\right]$

The strain rate therefore follows as

$\left[\frac{d{\epsilon }_{ij}}{dt}\right]=\frac{d\gamma }{dt}\left[\begin{array}{ccc}-\frac{1}{2}& 0& 0\\ 0& \frac{1}{2}& 0\\ 0& 0& 0\end{array}\right]$

Thus, $\frac{d{\epsilon }_{1}}{dt}=-\frac{d{\epsilon }_{2}}{dt}$, showing that the plastic strain rate is normal to the yield locus (see the figure).You could verify as an exercise that that if the stress reaches the other limiting surfaces of the yield locus, the resulting strain rate will be normal to the yield locus.