Hypoelasticity (Jaumann)

Recall that \boldsymbol{\sigma} and \mathbf{B} were shown to be work-conjugate, and this permitted us to develop a constitutive relationship between the two tensors. Now we would like to develop a constitutive relationship in hypoelastic form. Thus, we need a stress rate and a strain rate that are work-conjugate. Unfortunately, the time derivative of the Cauchy stress tensor, \boldsymbol{\dot{\sigma}}, is not work-conjugate with \mathbf{D}, where we recall that \mathbf{D} is the “rate of deformation tensor.” Consider the following:

(1)   \begin{equation*} \overset{\bigtriangledown}{\boldsymbol{\sigma}}=\frac{1}{det\mathbf{F}}\mathbf{F}\cdot\dot{\hat{\boldsymbol{\sigma}}}\cdot\mathbf{F}^{T}=\frac{1}{det\mathbf{F}}\mathbf{F}\frac{d}{dt}\left[(det\mathbf{F})\mathbf{F}^{-1} \cdot \boldsymbol{\sigma} \cdot \mathbf{F}^{-T}\right]\mathbf{F}^{T} \end{equation*}

Here, \overset{\bigtriangledown}{\boldsymbol{\sigma}} is the “Truesdell rate.” Essentially, we wrote \boldsymbol{\sigma} in terms of the Second Piola-Kirchhoff Stress, and then, within this expression, we took the time derivative of only the Second Piola-Kirchhoff Stress (recall from eq. 2 in Section 4: Alternative Measures of Stress, that the Second-Piola Kirchhoff Stress, \hat{\boldsymbol{\sigma}}, is (det\mathbf{F})\mathbf{F}^{-1} \cdot \boldsymbol{\sigma} \cdot \mathbf{F}^{-T} ).

The time derivative expressed in eq. 1 can be taken, which requires the use of the product rule for the four terms within the square brackets. An alternative derivation is given in Appendix D.2. In either case, eq. 1 becomes:

(2)   \begin{equation*} \overset{\bigtriangledown}{\boldsymbol{\sigma}}=\dot{\boldsymbol{\sigma}}-\mathbf{L} \cdot \boldsymbol{\sigma}-\boldsymbol{\sigma} \cdot \mathbf{L}^{T}+tr(\mathbf{L})\boldsymbol{\sigma} \end{equation*}

Recall the definition of the velocity gradient, \mathbf{L}, from eq. 1 in Section 3: Rate of Deformation and Spin Tensors. Additionally, note that tr\mathbf{L}=tr\mathbf{D}.

Eq. 2 is the most commonly used form of the Truesdell rate. In eq. 2, if we eliminate all of the \mathbf{D} terms except for those that appear in the constitutive relationship (we can also think of this simplification as taking \mathbf{L} \approx \dot{\mathbf{R}} \cdot \mathbf{R}^{T} \approx \mathbf{W}), then we arrive at the following result:

(3)   \begin{equation*} \mathring{\boldsymbol{\sigma}}=\dot{\boldsymbol{\sigma}}-\mathbf{W} \cdot \boldsymbol{\sigma}-\boldsymbol{\sigma} \cdot \mathbf{W}^{T} \end{equation*}

or

(4)   \begin{equation*} \mathring{\boldsymbol{\sigma}}=\dot{\boldsymbol{\sigma}}-\mathbf{W} \cdot \boldsymbol{\sigma}+\boldsymbol{\sigma} \cdot \mathbf{W} \end{equation*}


note: \mathbf{W}^{T}=-\mathbf{W}


In hypoelasticity, eq. 3 is called the Jaumann rate of Cauchy Stress. We will form our hypoelastic constitutive relationship in terms of this rate. Then, we’ll see how this constitutive relationship would be used to solve real FEA problems, through a simple three step time-stepping algorithm. But first, we need to prove that the Jaumann rate of Cauchy Stress is indeed work conjugate with \mathbf{D}i.e. does \mathring{\boldsymbol{\sigma}}^*=\mathbf{Q} \cdot \mathring{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T}?


Note that \hat{\boldsymbol{\sigma}} can represent the linear infinitesimal stress. It turns out that most FEA software use the Jaumann rate for linear infinitesimal stress as well. The reason for this will be explained at the end of this chapter, along with a brief disussion of linear infinitesimal elasticity time-stepping.


Recall \boldsymbol{\sigma}^*=\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T}

\dot{\boldsymbol{\sigma}}^*=\dot{\mathbf{Q}} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T}+\mathbf{Q} \cdot \dot{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T}+\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \dot{\mathbf{Q}}^{T}
=\dot{\mathbf{Q}} \cdot \underbrace{\mathbf{Q}^{-1} \cdot \mathbf{Q}}_{\mathbf{I}} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T}+\mathbf{Q} \cdot \dot{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T}+\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \underbrace{\mathbf{Q}^{T} \cdot \mathbf{Q}}_{\mathbf{I}} \cdot \dot{\mathbf{Q}}^{T}
=\dot{\mathbf{Q}} \cdot \mathbf{Q}^{-1} \cdot \boldsymbol{\sigma}^*+\mathbf{Q} \cdot \dot{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T}+\boldsymbol{\sigma}^* \cdot (\dot{\mathbf{Q}} \cdot \mathbf{Q}^{-1})^{T}

Since \mathbf{W}^*=\dot{\mathbf{Q}} \cdot \mathbf{Q}^{-1}+\mathbf{Q} \cdot \mathbf{W} \cdot \mathbf{Q}^{T} \longrightarrow \dot{\mathbf{Q}} \cdot \mathbf{Q}^{-1}=\mathbf{W}^*-\mathbf{Q \cdot W} \cdot \mathbf{Q}^{T}
\longrightarrow \dot{\boldsymbol{\sigma}}^*=(\mathbf{W}^*-\mathbf{Q \cdot W \cdot Q}^{T}) \cdot \boldsymbol{\sigma}^*+\mathbf{Q} \cdot \dot{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T} +\boldsymbol{\sigma}^* \cdot (\underbrace{-\mathbf{W}^*+\mathbf{Q \cdot W \cdot Q}^{T}}_{(\dot{\mathbf{Q}} \cdot \mathbf{Q}^{-1})^{T}})
Moving some of the terms to the other side:
\longrightarrow \dot{\boldsymbol{\sigma}}^*-\mathbf{W}^* \cdot \boldsymbol{\sigma}^*+\boldsymbol{\sigma}^* \cdot \mathbf{W}^*=\mathbf{Q} \cdot \dot{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T}+\boldsymbol{\sigma}^* \cdot \mathbf{Q} \cdot \mathbf{W} \cdot \mathbf{Q}^{T}-\mathbf{Q} \cdot \mathbf{W} \cdot \mathbf{Q}^{T} \cdot \boldsymbol{\sigma}^*

Since \mathbf{Q \cdot W \cdot Q}^{T} \cdot \boldsymbol{\sigma}^*=\mathbf{Q} \cdot \mathbf{W} \cdot \cancel{\mathbf{Q}^{T}} \cdot \underbrace{\cancel{\mathbf{Q}} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T}}_{\boldsymbol{\sigma}^*}=\mathbf{Q \cdot W} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T}
and \boldsymbol{\sigma}^* \cdot \mathbf{Q \cdot W \cdot Q}^{T}=\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T} \cdot \mathbf{Q \cdot W \cdot Q}^{T}=\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \mathbf{W} \cdot \mathbf{Q}^{T}

\longrightarrow \dot{\boldsymbol{\sigma}}^*-\mathbf{W}^* \cdot \boldsymbol{\sigma}^*+\boldsymbol{\sigma}^* \cdot \mathbf{W}^*=\mathbf{Q} \cdot (\dot{\boldsymbol{\sigma}}-\mathbf{W} \cdot \boldsymbol{\sigma}+\boldsymbol{\sigma} \cdot \mathbf{W}) \cdot \mathbf{Q}^{T}

Since we know that \mathring{\boldsymbol{\sigma}}=\dot{\boldsymbol{\sigma}}-\mathbf{W} \cdot \boldsymbol{\sigma} +\boldsymbol{\sigma} \cdot \mathbf{W}

(5)   \begin{equation*} \longrightarrow \mathring{\boldsymbol{\sigma}}^*=\mathbf{Q} \cdot \mathring{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T} \end{equation*}


Note the similarity to \boldsymbol{\sigma}^*=\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T} and \mathbf{D}^*=\mathbf{Q} \cdot \mathbf{D} \cdot \mathbf{Q}^{T}. This identical behavior under superimposed rigid body rotation is what we wanted to show. These special time derivatives of stress (or strain) are sometimes called Lie time derivatives. An alternative proof uses the deformation gradient instead of the arbitrary rigid body rotation, and takes advantage of the idea of the “pull-back operation” and “push-forward operation.” [Holzapfel] does a good job of explaining this latter approach in his section on Lie time derivatives.


Now, just as we have derived \boldsymbol{\sigma} \longleftrightarrow \mathbf{B} in the past, we can now derive \mathring{\boldsymbol{\sigma}} \longleftrightarrow \mathbf{D} (\mathring{\boldsymbol{\sigma}} will also depend on \boldsymbol{\sigma}, \mathring{\mathbf{B}}, and \mathbf{B}).

Truesdell et al. [Truesdell] coined the phrase “hypoelasticity.” All hyperelastic constitutive relationships can be put into rate form (and hence are hypoelastic). The reverse is not always possible [Simo]. Materials that are viscoelastic, or any material that is expected to experience plasticity, can be conveniently expressed in rate form, where history dependent variables can be tracked, naturally.


note:
\mathbf{B}^*=\mathbf{Q} \cdot \mathbf{B} \cdot \mathbf{Q}^{T}
\boldsymbol{\sigma}^*=\mathbf{Q} \cdot \boldsymbol{\sigma} \cdot \mathbf{Q}^{T}
\mathbf{D}^*=\mathbf{Q} \cdot \mathbf{D} \cdot \mathbf{Q}^{T}
\mathring{\boldsymbol{\sigma}}^*=\mathbf{Q} \cdot \mathring{\boldsymbol{\sigma}} \cdot \mathbf{Q}^{T}


Since all of the tensors above are work-conjugate, it is convenient to use them to derive our final form of \mathring{\boldsymbol{\sigma}}.

Recall \boldsymbol{\sigma}=\frac{1}{det\mathbf{B}^{1/2}}(\mathbf{B}\frac{d\phi}{d\mathbf{B}}+\frac{d\phi}{d\mathbf{B}}\mathbf{B}):

We will form our hypoelastic constitutive relationship by applying the Jaumann “operator” (\mathring{\boldsymbol{\mathcal{X}}}) to both sides of this expression. But, before we do this, we ought to find \mathring{det\mathbf{F}} and \mathring{\mathbf{B}}.

\mathring{det\mathbf{F}}=\dot{det\mathbf{F}} (i.e. presumably, the analog to -\mathbf{W \cdot B + B \cdot W} would simply be zero, since any asymmetries cancel out when taking the determinant)

From eq. 1 in Section 2: Volume and Area Change we know that dV=dV_0det\mathbf{F} and from eq. 1 in Section 3: Other Rates of Change we know that \frac{d}{dt}(dV)=D_{kk}dV.

So, \frac{d}{dt}dV=D_{kk}dV=\dot{det\mathbf{F}}dV_0. Thus,

(6)   \begin{equation*} \dot{det\mathbf{F}}=\underbrace{\frac{dV}{dV_0}}_{det\mathbf{F}}D_{kk}=(det\mathbf{F})tr\mathbf{D} \end{equation*}

In terms of \mathbf{B}: det\mathbf{F}=det\mathbf{B}^{1/2} \longrightarrow \frac{d}{dt}det\mathbf{B}^{1/2}=det\mathbf{B}^{1/2}tr\mathbf{D}

Next, we need to find \mathring{\mathbf{B}}:

\mathbf{B}=\mathbf{F} \cdot \mathbf{F}^{T}
\longrightarrow \dot{\mathbf{B}}=\dot{\mathbf{F}} \cdot \mathbf{F}^{T}+\mathbf{F} \cdot \dot{\mathbf{F}}^{T}=\underbrace{\dot{\mathbf{F}} \cdot \mathbf{F}^{-1}}_{\mathbf{L}} \cdot \underbrace{\mathbf{F} \cdot \mathbf{F}^{T}}_{\mathbf{B}} +\mathbf{F} \cdot \overbrace{\mathbf{F}^{T} \cdot \mathbf{F}^{-T}}^{\mathbf{I}} \cdot \dot{\mathbf{F}}^{T}
=\mathbf{L} \cdot \mathbf{B} + \mathbf{B} \cdot \mathbf{L}^{T}

\dot{\mathbf{B}}=\mathbf{L} \cdot \mathbf{B} + \mathbf{B} \cdot \mathbf{L}^{T}=\mathbf{D} \cdot \mathbf{B} +\mathbf{B} \cdot \mathbf{D}+\mathbf{W} \cdot \mathbf{B} -\mathbf{B} \cdot \mathbf{W}, (since \mathbf{L}=\mathbf{D}+\mathbf{W}
and \mathbf{L}^{T} = \mathbf{D}^{T} +\mathbf{W}^{T}=\mathbf{D-W})

\mathring{\mathbf{B}}=\underbrace{(\mathbf{D \cdot B} + \mathbf{B \cdot D} +\mathbf{W \cdot B} -\mathbf{B \cdot W})}_{\dot{\mathbf{B}}}-\mathbf{W \cdot B} + \mathbf{B \cdot W}=\mathbf{D \cdot B +B \cdot D}


Reminder: The direct calculation of the time derivative of \mathbf{B} (i.e. \dot{\mathbf{B}}), similar to \dot{\boldsymbol{\sigma}}, is not work-conjugate with \mathbf{D}. This is why we need to use \mathring{\mathbf{B}}, just as we are using \mathring{\boldsymbol{\sigma}}.


det\mathbf{B}^{1/2} \boldsymbol{\sigma}=\mathbf{B}\frac{d\phi}{d\mathbf{B}}+\frac{d\phi}{d\mathbf{B}}\mathbf{B} | apply the Jaumann operator (\mathring{\boldsymbol{\mathcal{X}}}) to this entire expression.

Noting that \underbrace{\mathbf{B}\frac{d\phi}{d\mathbf{B}}=\mathring{\mathbf{B}}\frac{d\phi}{d\mathbf{B}}+\mathbf{B}\frac{d^{2}\phi}{d\mathbf{B}^{2}}:\mathring{\mathbf{B}}}_{\text{product rule and chain rule}}

\longrightarrow \boldsymbol{\sigma}det\mathbf{B}^{1/2}tr\mathbf{D}+det\mathbf{B}^{1/2}\mathring{\boldsymbol{\sigma}}=\mathring{\mathbf{B}}\frac{d\phi}{d\mathbf{B}}+\mathbf{B}\frac{d^{2}\phi}{d\mathbf{B}^{2}}:\mathring{\mathbf{B}}+(\frac{d^{2}\phi}{d\mathbf{B}^{2}}:\mathring{\mathbf{B}})\cdot \mathbf{B}+\frac{d\phi}{d\mathbf{B}}\mathring{\mathbf{B}}


note: Recall from Section Section 1. Trace, Scalar Product, Eigenvalues that \mathbf{A}:\mathbf{B}=tr(\mathbf{A} \cdot \mathbf{B}^{T})=A_{ij}B_{ij}


So, \mathring{\sigma}_{ij}=-\sigma_{ij}D_{kk}+A_{ijmn}\mathring{B}_{mn}

where (\text{skipped work}):

A_{ijmn}=\frac{1}{det\mathbf{B}^{1/2}}[\delta_{im}\frac{d\phi}{dB_{nj}}+\frac{d\phi}{dB_{im}}\delta_{jn}+B_{ik}\frac{d^{2}\phi}{dB_{kj}dB_{mn}}+\frac{d^{2}\phi}{dB_{ik}dB_{mn}}B_{kj}]

Since \mathring{\mathbf{B}}=\mathbf{B \cdot D}+\mathbf{D \cdot B}, all terms in the above expression involve \mathbf{D}. Thus, we arrive at the following relationship:

(7)   \begin{equation*} \mathring{\sigma}_{ij}=\Lambda_{ijkl}D_{kl} \end{equation*}

Eq. 7 is the Jaumann rate of Cauchy Stress in terms of deformation, stress, and rate of deformation.


Note that we could’ve formed a slightly different hypoelastic relationship using the Truesdell rate. Either way, the hypoelastic expression is going to be complicated (the Truesdell rate, slightly more so). This is, however, meant to be broad enough to treat materials that are both nonlinear elastic, and subjected to large strain.


  • C. Truesdell and W. Noll, The Non-Linear Field Theories of Mechanics, 3rd ed., Berlin, DE: Pringer, 2004.
    [Bibtex]
    @BOOK{Truesdell,
        Address        = {Berlin, DE},
        Author         = {Clifford Truesdell and Walter Noll},
        Edition        = {3rd},
        Publisher      = {Pringer},
        Title          = {The Non-Linear Field Theories of Mechanics},
        Year           = {2004}
    }
  • G. A. Holzapfel, Nonlinear Solid Mechanics: A Continuum Approach for Engineering, Baffins Lane, Chichester, West Sussex PO19 1UD, England: John Wiley & Sons Ltd., 2000.
    [Bibtex]
    @book{holzapfel,
    
        address = {Baffins Lane, Chichester, West Sussex PO19 1UD, England},
    
        author = {Holzapfel, Gerhard A.},
        publisher = {John Wiley \& Sons Ltd.},
    
        title = {{Nonlinear Solid Mechanics: A Continuum Approach for Engineering}},
    
        year = {2000}
    }
  • J. C. Simo and T. J. R. Hughes, Computational inelasticity, Springer, 1998.
    [Bibtex]
    @book{hughes,
      title={Computational inelasticity},
      author={Simo, J.C. and Hughes, T.J.R.},
      Publisher = {Springer},
      year={1998}
    }
  • G. Jaumann, “Geschlossenes System physikalischer und chemischer Differentialgesetze,” Sitzungsberichte Akad. Wiss. Wien, IIa, pp. 385-530, 1911.
    [Bibtex]
    @article{Jaumann,
      title={Geschlossenes {S}ystem physikalischer und chemischer {D}ifferentialgesetze},
      author={Jaumann, Gustav},
      journal={Sitzungsberichte Akad. Wiss. Wien, IIa},
      pages={385--530},
      year={1911}
    }