Objective stress rate

Under rigid body rotations (${\displaystyle {\boldsymbol {Q}}}$), the Cauchy stress tensor ${\displaystyle {\boldsymbol {\sigma }}}$ transforms as

$${\displaystyle {\boldsymbol {\sigma }}_{r}={\boldsymbol {Q}}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {Q}}^{T}~;~~{\boldsymbol {Q}}\cdot {\boldsymbol {Q}}^{T}={\boldsymbol {\mathit {1}}}}$$

Since ${\displaystyle {\boldsymbol {\sigma }}}$ is a spatial quantity and the transformation follows the rules of tensor transformations, ${\displaystyle {\boldsymbol {\sigma }}}$ is objective. However,

$${\displaystyle {\cfrac {d}{dt}}({\boldsymbol {\sigma }}_{r})={\dot {\boldsymbol {\sigma }}}_{r}={\dot {\boldsymbol {Q}}}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {Q}}^{T}+{\boldsymbol {Q}}\cdot {\dot {\boldsymbol {\sigma }}}\cdot {\boldsymbol {Q}}^{T}+{\boldsymbol {Q}}\cdot {\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {Q}}}^{T}\neq {\boldsymbol {Q}}\cdot {\dot {\boldsymbol {\sigma }}}\cdot {\boldsymbol {Q}}^{T}\,.}$$

Therefore, the stress rate is not objective unless the rate of rotation is zero, i.e. ${\displaystyle {\boldsymbol {Q}}}$ is constant.

Figure 1. Undeformed and deformed material element, and an elemental cube cut out from the deformed element.

For a physical understanding of the above, consider the situation shown in Figure 1. In the figure the components of the Cauchy (or true) stress tensor are denoted by the symbols ${\displaystyle S_{ij}}$. This tensor, which describes the forces on a small material element imagined to be cut out from the material as currently deformed, is not objective at large deformations because it varies with rigid body rotations of the material. The material points must be characterized by their initial Lagrangian coordinates ${\displaystyle X_{i}}$. Consequently, it is necessary to introduce the so-called objective stress rate ${\displaystyle {\overset {\circ }{S}}_{ij}}$, or the corresponding increment ${\displaystyle \Delta S_{ij}={\overset {\circ }{S}}_{ij}\Delta t}$. The objectivity is necessary for ${\displaystyle {\overset {\circ }{S}}_{ij}}$ to be functionally related to the element deformation. It means that ${\displaystyle {\overset {\circ }{S}}_{ij}}$ must be invariant with respect to coordinate transformations, particularly the rigid-body rotations, and must characterize the state of the same material element as it deforms.

The objective stress rate can be derived in two ways:

• by tensorial coordinate transformations,[5] which is the standard way in finite element textbooks[6]
• variationally, from strain energy density in the material expressed in terms of the strain tensor (which is objective by definition)[7][8]

While the former way is instructive and provides useful geometric insight, the latter way is mathematically shorter and has the additional advantage of automatically ensuring energy conservation, i.e., guaranteeing that the second-order work of the stress increment tensor on the strain increment tensor be correct (work conjugacy requirement).

The relation between the Cauchy stress and the 2nd P-K stress is called the Piola transformation. This transformation can be written in terms of the pull-back of ${\displaystyle {\boldsymbol {\sigma }}}$ or the push-forward of ${\displaystyle {\boldsymbol {S}}}$ as

$${\displaystyle {\boldsymbol {S}}=J~\phi ^{*}[{\boldsymbol {\sigma }}]~;~~{\boldsymbol {\sigma }}=J^{-1}~\phi _{*}[{\boldsymbol {S}}]}$$

The Truesdell rate of the Cauchy stress is the Piola transformation of the material time derivative of the 2nd P-K stress. We thus define

$${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~\phi _{*}[{\dot {\boldsymbol {S}}}]}$$

Expanded out, this means that

$${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~{\boldsymbol {F}}\cdot {\dot {\boldsymbol {S}}}\cdot {\boldsymbol {F}}^{T}=J^{-1}~{\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left(J~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}=J^{-1}~{\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]}$$

where the Kirchhoff stress ${\displaystyle {\boldsymbol {\tau }}=J~{\boldsymbol {\sigma }}}$ and the Lie derivative of the Kirchhoff stress is

$${\displaystyle {\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]={\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\tau }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}~.}$$

This expression can be simplified to the well known expression for the Truesdell rate of the Cauchy stress

Truesdell rate of the Cauchy stress $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}-{\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}-{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}+{\text{tr}}({\boldsymbol {l}})~{\boldsymbol {\sigma }}}$$ where ${\displaystyle {\boldsymbol {l}}}$ is the velocity gradient: ${\displaystyle {\boldsymbol {l}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {F}}^{-1}}$. Proof We start with $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~{\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left(J~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}~.}$$ Expanding the derivative inside the square brackets, we get $${\displaystyle {\begin{aligned}{\overset {\circ }{\boldsymbol {\sigma }}}&=J^{-1}~{\boldsymbol {F}}\cdot ({\dot {J}}~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T})\cdot {\boldsymbol {F}}^{T}+J^{-1}~{\boldsymbol {F}}\cdot (J~{\dot {{\boldsymbol {F}}^{-1}}}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T})\cdot {\boldsymbol {F}}^{T}\\&+J^{-1}~{\boldsymbol {F}}\cdot (J~{\boldsymbol {F}}^{-1}\cdot {\dot {\boldsymbol {\sigma }}}\cdot {\boldsymbol {F}}^{-T})\cdot {\boldsymbol {F}}^{T}+J^{-1}~{\boldsymbol {F}}\cdot (J~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\dot {{\boldsymbol {F}}^{-T}}})\cdot {\boldsymbol {F}}^{T}\end{aligned}}}$$ or, $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~{\dot {J}}~{\boldsymbol {\sigma }}+{\boldsymbol {F}}\cdot {\dot {{\boldsymbol {F}}^{-1}}}\cdot {\boldsymbol {\sigma }}+{\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\dot {{\boldsymbol {F}}^{-T}}}\cdot {\boldsymbol {F}}^{T}}$$ Now, $${\displaystyle {\boldsymbol {F}}\cdot {\boldsymbol {F}}^{-1}={\boldsymbol {\mathit {1}}}}$$ Therefore, $${\displaystyle {\frac {d}{dt}}\left({\boldsymbol {F}}\cdot {\boldsymbol {F}}^{-1}\right)=0\quad \implies \quad {\dot {\boldsymbol {F}}}\cdot {\boldsymbol {F}}^{-1}+{\boldsymbol {F}}\cdot {\dot {{\boldsymbol {F}}^{-1}}}=0}$$ or, $${\displaystyle {\dot {{\boldsymbol {F}}^{-1}}}=-{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {l}}\quad \implies \quad {\dot {{\boldsymbol {F}}^{-T}}}=-{\boldsymbol {l}}^{T}\cdot {\boldsymbol {F}}^{-T}}$$ where the velocity gradient ${\displaystyle {\boldsymbol {l}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {F}}^{-1}}$. Also, the rate of change of volume is given by $${\displaystyle {\dot {J}}=J~{\text{tr}}({\boldsymbol {d}})=J~{\text{tr}}({\boldsymbol {l}})}$$ where ${\displaystyle {\boldsymbol {d}}}$ is the rate of deformation tensor. Therefore, $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~J~{\text{tr}}({\boldsymbol {l}})~{\boldsymbol {\sigma }}-{\boldsymbol {F}}\cdot {\boldsymbol {F}}^{-1}\cdot {\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}+{\dot {\boldsymbol {\sigma }}}-{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}\cdot {\boldsymbol {F}}^{-T}\cdot {\boldsymbol {F}}^{T}}$$ or, $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}-{\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}-{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}+{\text{tr}}({\boldsymbol {l}})~{\boldsymbol {\sigma }}}$$

It can be shown that the Truesdell rate is objective.

The Truesdell rate of the Kirchhoff stress can be obtained by noting that

$${\displaystyle {\boldsymbol {S}}=\phi ^{*}[{\boldsymbol {\tau }}]~;~~{\boldsymbol {\tau }}=\phi _{*}[{\boldsymbol {S}}]}$$

and defining

$${\displaystyle {\overset {\circ }{\boldsymbol {\tau }}}=\phi _{*}[{\dot {\boldsymbol {S}}}]}$$

Expanded out, this means that

$${\displaystyle {\overset {\circ }{\boldsymbol {\tau }}}={\boldsymbol {F}}\cdot {\dot {\boldsymbol {S}}}\cdot {\boldsymbol {F}}^{T}={\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\tau }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}={\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]}$$

Therefore, the Lie derivative of ${\displaystyle {\boldsymbol {\tau }}}$ is the same as the Truesdell rate of the Kirchhoff stress.

Following the same process as for the Cauchy stress above, we can show that

Truesdell rate of the Kirchhoff stress $${\displaystyle {\overset {\circ }{\boldsymbol {\tau }}}={\dot {\boldsymbol {\tau }}}-{\boldsymbol {l}}\cdot {\boldsymbol {\tau }}-{\boldsymbol {\tau }}\cdot {\boldsymbol {l}}^{T}}$$

This is a special form of the Lie derivative (or the Truesdell rate of the Cauchy stress). Recall that the Truesdell rate of the Cauchy stress is given by

$${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~{\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left(J~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}~.}$$

From the polar decomposition theorem we have

$${\displaystyle {\boldsymbol {F}}={\boldsymbol {R}}\cdot {\boldsymbol {U}}}$$

where ${\displaystyle {\boldsymbol {R}}}$ is the orthogonal rotation tensor (${\displaystyle {\boldsymbol {R}}^{-1}={\boldsymbol {R}}^{T}}$) and ${\displaystyle {\boldsymbol {U}}}$ is the symmetric, positive definite, right stretch.

If we assume that ${\displaystyle {\boldsymbol {U}}={\boldsymbol {\mathit {1}}}}$ we get ${\displaystyle {\boldsymbol {F}}={\boldsymbol {R}}}$. Also since there is no stretch ${\displaystyle J=1}$ and we have ${\displaystyle {\boldsymbol {\tau }}={\boldsymbol {\sigma }}}$. Note that this doesn't mean that there is not stretch in the actual body - this simplification is just for the purposes of defining an objective stress rate. Therefore,

$${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\boldsymbol {R}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {R}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {R}}^{-T}\right)\right]\cdot {\boldsymbol {R}}^{T}={\boldsymbol {R}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {R}}^{T}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {R}}\right)\right]\cdot {\boldsymbol {R}}^{T}}$$

We can show that this expression can be simplified to the commonly used form of the Green-Naghdi rate

Green-Naghdi rate of the Cauchy stress $${\displaystyle {\overset {\square }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {\Omega }}-{\boldsymbol {\Omega }}\cdot {\boldsymbol {\sigma }}}$$ where ${\displaystyle {\boldsymbol {\Omega }}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}}$. Proof Expanding out the derivative $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\boldsymbol {R}}\cdot {\dot {{\boldsymbol {R}}^{T}}}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {R}}\cdot {\boldsymbol {R}}^{T}+{\boldsymbol {R}}\cdot {\boldsymbol {R}}^{T}\cdot {\dot {\boldsymbol {\sigma }}}\cdot {\boldsymbol {R}}\cdot {\boldsymbol {R}}^{T}+{\boldsymbol {R}}\cdot {\boldsymbol {R}}^{T}\cdot {\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}}$$ or, $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\boldsymbol {R}}\cdot {\dot {{\boldsymbol {R}}^{T}}}\cdot {\boldsymbol {\sigma }}+{\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}}$$ Now, $${\displaystyle {\boldsymbol {R}}\cdot {\boldsymbol {R}}^{T}={\boldsymbol {\mathit {1}}}\quad \implies \quad {\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}=-{\boldsymbol {R}}\cdot {\dot {{\boldsymbol {R}}^{T}}}}$$ Therefore, $${\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}-{\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}\cdot {\boldsymbol {\sigma }}}$$ If we define the angular velocity as $${\displaystyle {\boldsymbol {\Omega }}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}}$$ we get the commonly used form of the Green–Naghdi rate $${\displaystyle {\overset {\square }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {\Omega }}-{\boldsymbol {\Omega }}\cdot {\boldsymbol {\sigma }}}$$

The Green–Naghdi rate of the Kirchhoff stress also has the form since the stretch is not taken into consideration, i.e.,

$${\displaystyle {\overset {\square }{\boldsymbol {\tau }}}={\dot {\boldsymbol {\tau }}}+{\boldsymbol {\tau }}\cdot {\boldsymbol {\Omega }}-{\boldsymbol {\Omega }}\cdot {\boldsymbol {\tau }}}$$

The Zaremba-Jaumann rate of the Cauchy stress is a further specialization of the Lie derivative (Truesdell rate). This rate has the form

Zaremba-Jaumann rate of the Cauchy stress $${\displaystyle {\overset {\triangle }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {w}}-{\boldsymbol {w}}\cdot {\boldsymbol {\sigma }}}$$ where ${\displaystyle {\boldsymbol {w}}}$ is the spin tensor.

The Zaremba-Jaumann rate is used widely in computations primarily for two reasons

1. it is relatively easy to implement.
2. it leads to symmetric tangent moduli.

Recall that the spin tensor ${\displaystyle {\boldsymbol {w}}}$ (the skew part of the velocity gradient) can be expressed as

$${\displaystyle {\boldsymbol {w}}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}+{\frac {1}{2}}~{\boldsymbol {R}}\cdot ({\dot {\boldsymbol {U}}}\cdot {\boldsymbol {U}}^{-1}-{\boldsymbol {U}}^{-1}\cdot {\dot {\boldsymbol {U}}})\cdot {\boldsymbol {R}}^{T}}$$

Thus for pure rigid body motion

$${\displaystyle {\boldsymbol {w}}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}={\boldsymbol {\Omega }}}$$

Alternatively, we can consider the case of proportional loading when the principal directions of strain remain constant. An example of this situation is the axial loading of a cylindrical bar. In that situation, since

$${\displaystyle {\boldsymbol {U}}={\begin{bmatrix}\lambda _{X}\\&\lambda _{Y}\\&&\lambda _{Z}\end{bmatrix}}}$$

we have

$${\displaystyle {\dot {\boldsymbol {U}}}={\begin{bmatrix}{\dot {\lambda }}_{X}\\&{\dot {\lambda }}_{Y}\\&&{\dot {\lambda }}_{Z}\end{bmatrix}}}$$

Also,

$${\displaystyle {\boldsymbol {U}}^{-1}={\begin{bmatrix}1/\lambda _{X}\\&1/\lambda _{Y}\\&&1/\lambda _{Z}\end{bmatrix}}}$$

of the Cauchy stress. Therefore,

$${\displaystyle {\dot {\boldsymbol {U}}}\cdot {\boldsymbol {U}}^{-1}={\begin{bmatrix}{\dot {\lambda }}_{X}/\lambda _{X}\\&{\dot {\lambda }}_{Y}/\lambda _{Y}\\&&{\dot {\lambda }}_{Z}/\lambda _{Z}\end{bmatrix}}=U^{-1}{\dot {U}}}$$

This once again gives

$${\displaystyle {\boldsymbol {w}}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}={\boldsymbol {\Omega }}}$$

In general, if we approximate

$${\displaystyle {\boldsymbol {w}}\approx {\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}}$$

the Green–Naghdi rate becomes the Zaremba-Jaumann rate of the Cauchy stress

$${\displaystyle {\overset {\triangle }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {w}}-{\boldsymbol {w}}\cdot {\boldsymbol {\sigma }}}$$

There can be an infinite variety of objective stress rates. One of these is the Oldroyd stress rate

$${\displaystyle {\overset {\triangledown }{\boldsymbol {\sigma }}}={\mathcal {L}}_{\varphi }[{\boldsymbol {\sigma }}]={\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}}$$

In simpler form, the Oldroyd rate is given by

$${\displaystyle {\overset {\triangledown }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}-{\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}-{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}}$$

If the current configuration is assumed to be the reference configuration then the pull back and push forward operations can be conducted using ${\displaystyle {\boldsymbol {F}}^{T}}$ and ${\displaystyle {\boldsymbol {F}}^{-T}}$ respectively. The Lie derivative of the Cauchy stress is then called the convective stress rate

$${\displaystyle {\overset {\diamond }{\boldsymbol {\sigma }}}={\boldsymbol {F}}^{-T}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{T}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}\right)\right]\cdot {\boldsymbol {F}}^{-1}}$$

In simpler form, the convective rate is given by

$${\displaystyle {\overset {\diamond }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}}$$

Many materials undergo inelastic deformations caused by plasticity and damage. These material behaviors cannot be described in terms of a potential. It is also often the case that no memory of the initial virgin state exists, particularly when large deformations are involved.[9] The constitutive relation is typically defined in incremental form in such cases to make the computation of stresses and deformations easier.[10]

The incremental loading procedure

For a small enough load step, the material deformation can be characterized by the small (or linearized) strain increment tensor[11]

$${\displaystyle {\boldsymbol {e}}={\tfrac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}\right]\quad \equiv \quad e_{ij}={\tfrac {1}{2}}(u_{i,j}+u_{j,i})}$$

where ${\displaystyle \mathbf {u} }$ is the displacement increment of the continuum points. The time derivative

$${\displaystyle {\frac {\partial {\boldsymbol {e}}}{\partial t}}={\dot {\boldsymbol {e}}}={\tfrac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {v} +({\boldsymbol {\nabla }}\mathbf {v} )^{T}\right]\quad \equiv \quad {\dot {e}}_{ij}={\tfrac {1}{2}}(v_{i,j}+v_{j,i})}$$

is the strain rate tensor (also called the velocity strain) and ${\displaystyle \mathbf {v} ={\dot {\mathbf {u} }}}$ is the material point velocity or displacement rate. For finite strains, measures from the Seth–Hill family (also called Doyle–Ericksen tensors) can be used:

$${\displaystyle \mathbf {E} _{(m)}={\frac {1}{2m}}(\mathbf {U} ^{2m}-\mathbf {I} )}$$

where ${\displaystyle \mathbf {U} }$ is the right stretch. A second-order approximation of these tensors is

$${\displaystyle \mathbf {E} _{(m)}\approx {\boldsymbol {e}}+{\tfrac {1}{2}}(\nabla \mathbf {u} )^{T}\cdot \nabla \mathbf {u} -(1-m){\boldsymbol {e}}\cdot {\boldsymbol {e}}}$$

Variation of work done

Then the variation in work done can be expressed as

$${\displaystyle \delta W={\boldsymbol {S}}_{(m)}:\delta {\boldsymbol {E}}_{(m)}={\boldsymbol {P}}:\delta \nabla \mathbf {u} }$$

where the finite strain measure ${\displaystyle {\boldsymbol {E}}_{(m)}}$ is energy conjugate to the stress measure ${\displaystyle {\boldsymbol {\sigma }}^{(m)}}$. Expanded out,

$${\displaystyle \delta W=\left({\boldsymbol {S}}+{\boldsymbol {\sigma }}_{0}\right):\delta {\boldsymbol {E}}_{(m)}=\left({\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0}\right):\delta \nabla \mathbf {u} \,.}$$

The objectivity of stress tensor ${\displaystyle {\boldsymbol {S}}_{(m)}}$ is ensured by its transformation as a second-order tensor under coordinate rotations (which causes the principal stresses to be independent from coordinate rotations) and by the correctness of ${\displaystyle {\boldsymbol {S}}_{(m)}:\delta {\boldsymbol {E}}_{(m)}}$ as a second-order energy expression.

From the symmetry of the Cauchy stress, we have

$${\displaystyle {\boldsymbol {\sigma }}_{0}:\delta \nabla \mathbf {u} ={\boldsymbol {\sigma }}_{0}:\delta {\boldsymbol {e}}\,.}$$

For small variations in strain, using the approximation

$${\displaystyle {\boldsymbol {S}}:\delta {\boldsymbol {E}}_{(m)}\approx {\boldsymbol {S}}:\delta \nabla \mathbf {u} }$$

and the expansions

$${\displaystyle {\boldsymbol {\sigma }}_{0}:\delta {\boldsymbol {E}}_{(m)}={\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]~,~~{\boldsymbol {\sigma }}_{0}:\delta {\boldsymbol {e}}={\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]}$$

we get the equation

$${\displaystyle {\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]+{\boldsymbol {S}}:\delta \nabla \mathbf {u} ={\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]+{\boldsymbol {T}}:\delta \nabla \mathbf {u} \,.}$$

Imposing the variational condition that the resulting equation must be valid for any strain gradient ${\displaystyle \delta \nabla \mathbf {u} }$, we have [7]

$${\displaystyle {\boldsymbol {S}}={\boldsymbol {T}}-{\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}-{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}\right]}$$

(1)

We can also write the above equation as

$${\displaystyle {\boldsymbol {S}}_{(m)}={\boldsymbol {P}}-{\boldsymbol {\sigma }}_{0}:{\frac {\partial }{\partial \nabla \mathbf {u} }}\left[{\boldsymbol {E}}_{(m)}-{\boldsymbol {e}}\right]\,.}$$

(2)

Time derivatives

The Cauchy stress and the first Piola-Kirchhoff stress are related by (see Stress measures)

$${\displaystyle {\boldsymbol {\sigma }}={\boldsymbol {P}}\cdot {\boldsymbol {F}}^{T}J^{-1}=({\boldsymbol {P}}+{\boldsymbol {P}}\cdot \nabla \mathbf {u} ^{T})J^{-1}\,.}$$

For small incremental deformations,

$${\displaystyle J^{-1}\approx 1-\nabla \cdot \mathbf {u} \,.}$$

Therefore,

$${\displaystyle \Delta {\boldsymbol {\sigma }}={\boldsymbol {\sigma }}-{\boldsymbol {\sigma }}_{0}\approx ({\boldsymbol {P}}+{\boldsymbol {P}}\cdot \nabla \mathbf {u} ^{T})(1-\nabla \cdot \mathbf {u} )-{\boldsymbol {\sigma }}_{0}\,.}$$

Substituting ${\displaystyle {\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0}={\boldsymbol {P}}}$,

$${\displaystyle \Delta {\boldsymbol {\sigma }}\approx [{\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0}+({\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0})\cdot \nabla \mathbf {u} ^{T}](1-\nabla \cdot \mathbf {u} )-{\boldsymbol {\sigma }}_{0}\,.}$$

For small increments of stress ${\displaystyle {\boldsymbol {T}}}$ relative to the initial stress ${\displaystyle {\boldsymbol {\sigma }}_{0}}$, the above reduces to

$${\displaystyle \Delta {\boldsymbol {\sigma }}\approx {\boldsymbol {T}}-{\boldsymbol {\sigma }}_{0}(\nabla \cdot \mathbf {u} )+{\boldsymbol {\sigma }}_{0}\cdot \nabla \mathbf {u} ^{T}\,.}$$

(3)

From equations (1) and (3) we have

$${\displaystyle {\boldsymbol {S}}=\Delta {\boldsymbol {\sigma }}+{\boldsymbol {\sigma }}_{0}(\nabla \cdot \mathbf {u} )-{\boldsymbol {\sigma }}_{0}\cdot \nabla \mathbf {u} ^{T}-{\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}-{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}\right]}$$

(4)

Recall that ${\displaystyle {\boldsymbol {S}}}$ is an increment of the stress tensor measure ${\displaystyle {\boldsymbol {S}}_{(m)}}$. Defining the stress rate

$${\displaystyle {\boldsymbol {S}}=:{\overset {\circ }{\boldsymbol {S}}}_{(m)}\Delta t}$$

and noting that

$${\displaystyle \Delta {\boldsymbol {\sigma }}={\dot {\boldsymbol {\sigma }}}\Delta t}$$

we can write equation (4) as

$${\displaystyle {\overset {\circ }{\boldsymbol {S}}}_{(m)}\Delta t={\dot {\boldsymbol {\sigma }}}\Delta t+{\boldsymbol {\sigma }}_{0}(\nabla \cdot \mathbf {v} )\Delta t-{\boldsymbol {\sigma }}_{0}\cdot \nabla \mathbf {v} ^{T}\Delta t-{\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}-{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}\right]}$$

(5)

Taking the limit at ${\displaystyle \Delta t\rightarrow 0}$, and noting that ${\displaystyle {\boldsymbol {\sigma }}_{0}={\boldsymbol {\sigma }}}$ at this limit, one gets the following expression for the objective stress rate associated with the strain measure ${\displaystyle {\boldsymbol {E}}_{(m)}}$:

$${\displaystyle {\overset {\circ }{\boldsymbol {S}}}_{(m)}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}(\nabla \cdot \mathbf {v} )-{\boldsymbol {\sigma }}\cdot \nabla \mathbf {v} ^{T}-{\boldsymbol {\sigma }}:{\frac {\partial }{\partial t}}\left[{\frac {\partial }{\partial \nabla \mathbf {u} }}\left({\boldsymbol {E}}_{(m)}-{\boldsymbol {e}}\right)\right]\,.}$$

(6)

Here ${\displaystyle {\dot {\sigma }}_{ij}=\partial \sigma _{ij}/\partial t}$ = material rate of Cauchy stress (i.e., the rate in Lagrangian coordinates of the initial stressed state).

Work-conjugate stress rates

A rate for which there exists no legitimate finite strain tensor ${\displaystyle {\boldsymbol {E}}_{(m)}}$ associated according to Eq. (6) is energetically inconsistent, i.e., its use violates energy balance (i.e., the first law of thermodynamics).

Evaluating Eq. (6) for general ${\displaystyle m}$ and for ${\displaystyle m=2}$, one gets a general expression for the objective stress rate:[7][8]

$${\displaystyle {\overset {\circ }{\boldsymbol {S}}}_{(m)}={\overset {\circ }{\boldsymbol {S}}}_{(2)}+{\tfrac {1}{2}}(2-m)[{\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {e}}}+({\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {e}}})^{T}]}$$

(7)

where ${\displaystyle {\overset {\circ }{\boldsymbol {S}}}_{(2)}}$ is the objective stress rate associated with the Green-Lagrangian strain (${\displaystyle m=2}$).

In particular,

• ${\displaystyle m=2}$ gives the Truesdell stress rate
• ${\displaystyle m=0}$ gives the Zaremba-Jaumann rate of Kirchhoff stress
• ${\displaystyle m=1}$ gives the Biot stress rate

(Note that m = 2 leads to Engesser's formula for critical load in shear buckling, while m = -2 leads to Haringx's formula which can give critical loads differing by >100%).

The tangential stress-strain relation has generally the form

$${\displaystyle {\dot {S}}_{ij}^{(m)}=C_{ijkl}^{(m)}{\dot {e}}_{kl}}$$

(6)

where ${\displaystyle C_{ijkl}^{(m)}}$ are the tangential moduli (components of a 4th-order tensor) associated with strain tensor ${\displaystyle \epsilon _{ij}^{(m)}}$. They are different for different choices of ${\displaystyle m}$, and are related as follows:

$${\displaystyle \left[C_{ijkl}^{(m)}-C_{ijkl}^{(2)}-{\tfrac {1}{4}}(2-m)(S_{ik}\delta _{jl}+S_{jk}\delta _{il}+S_{il}\delta _{jk}+S_{jl}\delta _{ik})\right]v_{k,l}=0}$$

(7)

From the fact that Eq. (7) must hold true for any velocity gradient ${\displaystyle v_{k,l}}$, it follows that:[7]

$${\displaystyle C_{ijkl}^{(m)}=C_{ijkl}^{(2)}+(2-m)[S_{ik}\delta _{jl}]_{\mathrm {sym} },~~[S_{ik}\delta _{jl}]_{\mathrm {sym} }={\tfrac {1}{4}}(S_{ik}\delta _{jl}+S_{jk}\delta _{il}+S_{il}\delta _{jk}+S_{jl}\delta _{ik})}$$

(8)

where ${\displaystyle C_{ijkl}^{(2)}}$ are the tangential moduli associated with the Green–Lagrangian strain (${\displaystyle m=2}$), taken as a reference, ${\displaystyle S_{ij}}$ = current Cauchy stress, and ${\displaystyle \delta _{ij}}$ = Kronecker delta (or unit tensor).

Eq. (8) can be used to convert one objective stress rate to another. Since ${\displaystyle S_{ij}{\dot {e}}_{kk}=(S_{ij}\delta _{kl})\delta e_{kl}}$, the transformation[7][8]

$${\displaystyle C_{ijkl}^{\mathrm {conj} }=C_{ijkl}^{\mathrm {nonconj} }+S_{ij}\delta _{kl}}$$

(9)

can further correct for the absence of the term ${\displaystyle S_{ij}v_{k,k}}$ (note that the term ${\displaystyle S_{ij}\delta _{km}}$ does not allow interchanging subscripts ${\displaystyle ij}$ with ${\displaystyle kl}$, which means that its absence breaks the major symmetry of the tangential moduli tensor ${\displaystyle C_{ijkl}^{\mathrm {nonconj} }}$).

Large strain often develops when the material behavior becomes nonlinear, due to plasticity or damage. Then the primary cause of stress dependence of the tangential moduli is the physical behavior of material. What Eq. (8) means that the nonlinear dependence of ${\displaystyle C_{ijkl}}$ on the stress must be different for different objective stress rates. Yet none of them is fundamentally preferable, except if there exists one stress rate, one ${\displaystyle m}$, for which the moduli can be considered constant.

• Cauchy stress tensor
• Stress measures
• Principle of material objectivity
• Hypoelastic material

• Objectivity in Classical Mechanics