# Nonlinear finite elements/Rate form of hyperelastic laws
## Rate equations

Let us now derive rate equations for a hyperelastic material.

### First elasticity tensor

We start off with the relation
    ${\boldsymbol {P}}=\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {F}}}}$
Then the material time derivative of ${\boldsymbol {P}}$ is given by
where the fourth order tensor ${\boldsymbol {\mathsf {A}}}$ is call the first elasticity tensor. This tensor has major symmetries but not minor symmetries. In index notation with respect to an orthonormal basis
    ${\dot {P}}_{iJ}={\mathsf {A}}_{iJkL}~{\dot {F}}_{kL}=\rho _{0}~{\frac {\partial ^{2}W}{\partial F_{iJ}\partial F_{kL}}}~{\dot {F}}_{kL}$
Proof:
   We have
       ${\dot {\boldsymbol {P}}}=\rho _{0}~{\frac {\partial }{\partial t}}\left({\frac {\partial W}{\partial {\boldsymbol {F}}}}\right)=\rho _{0}~{\frac {\partial }{\partial {\boldsymbol {F}}}}\left({\frac {\partial W}{\partial t}}\right)=\rho _{0}~{\frac {\partial }{\partial {\boldsymbol {F}}}}\left({\frac {\partial W}{\partial {\boldsymbol {F}}}}:{\dot {\boldsymbol {F}}}\right)$
   Using the product rule, we have
       ${\dot {\boldsymbol {P}}}=\rho _{0}~{\frac {\partial ^{2}W}{\partial {\boldsymbol {F}}\partial {\boldsymbol {F}}}}:{\dot {\boldsymbol {F}}}+\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {F}}}}:{\frac {\partial ^{2}{\boldsymbol {F}}}{\partial {\boldsymbol {F}}\partial t}}=\rho _{0}~{\frac {\partial ^{2}W}{\partial {\boldsymbol {F}}\partial {\boldsymbol {F}}}}:{\dot {\boldsymbol {F}}}+\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {F}}}}:{\frac {\partial }{\partial t}}\left({\frac {\partial {\boldsymbol {F}}}{\partial {\boldsymbol {F}}}}\right)$
   Therefore,
       ${\dot {\boldsymbol {P}}}=\rho _{0}~{\frac {\partial ^{2}W}{\partial {\boldsymbol {F}}\partial {\boldsymbol {F}}}}:{\dot {\boldsymbol {F}}}={\boldsymbol {\mathsf {A}}}:{\dot {\boldsymbol {F}}}\qquad \square$

### Second elasticity tensor

Similarly, if we start off with the relation
    ${\boldsymbol {S}}=2~\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {C}}}}$
the material time derivative of ${\boldsymbol {S}}$ can be expressed as
where the fourth order tensor ${\boldsymbol {\mathsf {C}}}$ is called the material elasticity tensor or the second elasticity tensor. Since this tensor relates symmetric second order tensors it has minor symmetries. It also has major symmetries because the two partial derivatives are with the same quantity and an interchange does not change things. In index notation with respect to an orthonormal basis
    ${\dot {S}}_{IJ}={\frac {1}{2}}~{\mathsf {C}}_{IJKL}~{\dot {C}}_{KL}={\frac {1}{2}}~\left(4~\rho _{0}~{\frac {\partial ^{2}W}{\partial C_{IJ}\partial C_{KL}}}\right)~{\dot {C}}_{KL}$
Proof:
   We have
       ${\dot {\boldsymbol {S}}}=2~\rho _{0}~{\frac {\partial }{\partial t}}\left({\frac {\partial W}{\partial {\boldsymbol {C}}}}\right)=2~\rho _{0}~{\frac {\partial }{\partial {\boldsymbol {C}}}}\left({\frac {\partial W}{\partial t}}\right)=2~\rho _{0}~{\frac {\partial }{\partial {\boldsymbol {C}}}}\left({\frac {\partial W}{\partial {\boldsymbol {C}}}}:{\dot {\boldsymbol {C}}}\right)$
   Again using the product rule, we have
       ${\dot {\boldsymbol {S}}}=2~\rho _{0}~{\frac {\partial ^{2}W}{\partial {\boldsymbol {C}}\partial {\boldsymbol {C}}}}:{\dot {\boldsymbol {C}}}+2~\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {C}}}}:{\frac {\partial ^{2}{\boldsymbol {C}}}{\partial {\boldsymbol {C}}\partial t}}=2~\rho _{0}~{\frac {\partial ^{2}W}{\partial {\boldsymbol {C}}\partial {\boldsymbol {C}}}}:{\dot {\boldsymbol {C}}}$
   Therefore,
       ${\dot {\boldsymbol {S}}}=4~\rho _{0}~{\frac {\partial ^{2}W}{\partial {\boldsymbol {C}}\partial {\boldsymbol {C}}}}:{\frac {1}{2}}~{\dot {\boldsymbol {C}}}={\boldsymbol {\mathsf {C}}}:{\frac {1}{2}}~{\dot {\boldsymbol {C}}}\qquad \square$
Proof:
   Recall that the first and second Piola-Kirchhoff stresses are related by
       ${\boldsymbol {P}}={\boldsymbol {F}}\cdot {\boldsymbol {S}}$
   Taking the material time derivative of both sides gives
       ${\dot {\boldsymbol {P}}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {S}}+{\boldsymbol {F}}\cdot {\dot {\boldsymbol {S}}}$
   Using the expression for ${\dot {\boldsymbol {S}}}$ above, we get
       ${\dot {\boldsymbol {P}}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {S}}+{\frac {1}{2}}~{\boldsymbol {F}}\cdot ({\boldsymbol {\mathsf {C}}}:{\dot {\boldsymbol {C}}})$
   Now
       ${\boldsymbol {C}}={\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}\quad \implies \quad {\dot {\boldsymbol {C}}}={\dot {{\boldsymbol {F}}^{T}}}\cdot {\boldsymbol {F}}+{\boldsymbol {F}}^{T}\cdot {\dot {\boldsymbol {F}}}$
   Therefore,
       ${\dot {\boldsymbol {P}}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {S}}+{\frac {1}{2}}~{\boldsymbol {F}}\cdot [{\boldsymbol {\mathsf {C}}}:({\dot {{\boldsymbol {F}}^{T}}}\cdot {\boldsymbol {F}}+{\boldsymbol {F}}^{T}\cdot {\dot {\boldsymbol {F}}})]$
   Now
       ${\boldsymbol {\mathsf {C}}}:({\dot {{\boldsymbol {F}}^{T}}}\cdot {\boldsymbol {F}})\equiv {\mathsf {C}}_{IJKL}~{\dot {F}}_{mK}~F_{mL}\quad {\text{and}}\qquad {\boldsymbol {\mathsf {C}}}:({\boldsymbol {F}}^{T}\cdot {\dot {\boldsymbol {F}}})\equiv ={\mathsf {C}}_{IJLK}~F_{mL}~{\dot {F}}_{mK}={\mathsf {C}}_{IJKL}~F_{mL}~{\dot {F}}_{mK}$
   That means
       ${\boldsymbol {\mathsf {C}}}:({\dot {{\boldsymbol {F}}^{T}}}\cdot {\boldsymbol {F}})={\boldsymbol {\mathsf {C}}}:({\boldsymbol {F}}^{T}\cdot {\dot {\boldsymbol {F}}})$
   which gives us
       ${\dot {\boldsymbol {P}}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {S}}+{\boldsymbol {F}}\cdot [{\boldsymbol {\mathsf {C}}}:({\boldsymbol {F}}^{T}\cdot {\dot {\boldsymbol {F}}})]$
   In index notation,
       ${\begin{aligned}{\dot {P}}_{iJ}&={\dot {F}}_{iL}~S_{LJ}+F_{iN}~[{\mathsf {C}}_{NJML}~F_{kM}~{\dot {F}}_{kL}]\\&=S_{LJ}~\delta _{ik}{\dot {F}}_{kL}+{\mathsf {C}}_{NJML}~F_{iN}~F_{kM}~{\dot {F}}_{kL}\\&=(\delta _{ik}~S_{JL}+F_{iN}~{\mathsf {C}}_{NJML}~F_{kM})~{\dot {F}}_{kL}&={\mathsf {A}}_{iJkL}~{\dot {F}}_{kL}\end{aligned}}$
   Therefore,
       ${\mathsf {A}}_{iJkL}=\delta _{ik}~S_{JL}+F_{iN}~{\mathsf {C}}_{NJML}~F_{kM}\qquad \square$

### Fourth elasticity tensor

Now we will compute the spatial elasticity tensor for the rate constitutive equation for a hyperelastic material. This tensor relates an objective rate of stress (Cauchy or Kirchhoff) to the rate of deformation tensor. We can show that
The fourth order tensor ${\boldsymbol {\mathsf {c}}}$ is called the spatial elasticity tensor or the fourth elasticity tensor. Clearly, ${\mathsf {c}}_{ijkl}$ cannot be derived from the store energy function $W$ because of the dependence on the deformation gradient.
Proof:
   Recall that the Lie derivative of the Kirchhoff stress is defined as
       ${\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]={\overset {\circ }{\boldsymbol {\tau }}}={\boldsymbol {F}}\cdot {\dot {\boldsymbol {S}}}\cdot {\boldsymbol {F}}^{T}$
   We have found that
       ${\dot {\boldsymbol {S}}}={\boldsymbol {\mathsf {C}}}:{\frac {1}{2}}~{\dot {\boldsymbol {C}}}$
   We also know from Continuum_mechanics/Time_derivatives_and_rates#Time_derivative_of_strain that
       ${\dot {\boldsymbol {C}}}=2~{\boldsymbol {F}}^{T}\cdot {\boldsymbol {d}}\cdot {\boldsymbol {F}}$
   where ${\boldsymbol {d}}$ is the spatial rate of deformation tensor. Therefore,
       ${\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]={\boldsymbol {F}}\cdot [{\boldsymbol {\mathsf {C}}}:({\boldsymbol {F}}^{T}\cdot {\boldsymbol {d}}\cdot {\boldsymbol {F}})]\cdot {\boldsymbol {F}}^{T}$
   In index notation,
       $({\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}])_{ij}=F_{iK}~{\mathsf {C}}_{KLMN}~F_{kM}~d_{kl}~F_{lN}~F_{jL}=F_{iK}~F_{jL}~F_{kM}~F_{lN}~{\mathsf {C}}_{KLMN}~d_{kl}={\mathsf {c}}_{ijkl}~d_{kl}$
   or,
       ${\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]={\boldsymbol {\mathsf {c}}}:{\boldsymbol {d}}$
   where
       ${\mathsf {c}}_{ijkl}=F_{iK}~F_{jL}~F_{kM}~F_{lN}~{\mathsf {C}}_{KLMN}\qquad \square$
Alternatively, we may define ${\boldsymbol {\mathsf {c}}}$ in terms of the Cauchy stress ${\boldsymbol {\sigma }}$, in which case the constitutive relation is written as
The proof of this relation between the spatial and material elasticity tensors is very similar to that for the rate of Kirchhoff stress. Many authors define this quantity ${\mathsf {c}}_{ijkl}$ as the spatial elasticity tensor. Note the factor of $J^{-1}$. This form of the spatial elasticity tensor is crucial for some of the calculations that follow.
Proof:
   Recall that
       ${\dot {\boldsymbol {P}}}={\boldsymbol {\mathsf {A}}}:{\dot {\boldsymbol {F}}}$
   Therefore,
       ${\dot {\boldsymbol {P}}}\cdot {\boldsymbol {F}}^{T}=({\boldsymbol {\mathsf {A}}}:{\dot {\boldsymbol {F}}})\cdot {\boldsymbol {F}}^{T}$
   Also recall that
       ${\mathsf {A}}_{iJkL}=\delta _{ik}~S_{JL}+F_{iN}~{\mathsf {C}}_{NJML}~F_{kM}$
   Therefore, using index notation,
       ${\dot {P}}_{iJ}~F_{lJ}=\delta _{ik}~S_{JL}~F_{lJ}~{\dot {F}}_{kL}+F_{iN}~{\mathsf {C}}_{NJML}~F_{kM}~F_{lJ})~{\dot {F}}_{kL}$
   Now,
       ${\boldsymbol {l}}={\dot {\boldsymbol {F}}}\cdot {\boldsymbol {F}}^{-1}\qquad {\text{or}}\qquad {\dot {\boldsymbol {F}}}={\boldsymbol {l}}\cdot {\boldsymbol {F}}$
   In index notation
       ${\dot {F}}_{kL}=l_{km}~F_{mL}$
   Using this we get
       ${\dot {P}}_{iJ}~F_{lJ}=\delta _{ik}~S_{JL}~F_{lJ}~l_{km}~F_{mL}+F_{iN}~{\mathsf {C}}_{NJML}~F_{kM}~F_{lJ})~l_{km}~F_{mL}$
   or,
       ${\dot {P}}_{iJ}~F_{lJ}=(\delta _{ik}~F_{lJ}~S_{JL}~F_{mL}+F_{iN}~F_{lJ}~F_{kM}~F_{mL}~{\mathsf {C}}_{NJML})~l_{km}$
   Now,
       ${\boldsymbol {\tau }}={\boldsymbol {F}}\cdot {\boldsymbol {S}}\cdot {\boldsymbol {F}}^{T}\qquad {\text{or}}\qquad \tau _{lm}=F_{lJ}~S_{JL}~F_{mL}$
   Therefore,
       ${\dot {P}}_{iJ}~F_{lJ}=(\delta _{ik}~\tau _{lm}+F_{iN}~F_{lJ}~F_{kM}~F_{mL}~{\mathsf {C}}_{NJML})~l_{km}$
   Also,
       ${\mathsf {c}}_{ilkm}=F_{iN}~F_{lJ}~F_{kM}~F_{mL}~{\mathsf {C}}_{NJML}$
   So we have
       ${\dot {P}}_{iJ}~F_{lJ}=(\delta _{ik}~\tau _{lm}+{\mathsf {c}}_{ilkm})~l_{km}={\mathsf {a}}_{ilkm}~l_{km}$
   Note:
      The fourth order tensor
          ${\mathsf {a}}_{ijkl}=\delta _{ik}~\tau _{jl}+{\mathsf {c}}_{ijkl}$
      which depends on the symmetry of ${\boldsymbol {\tau }}$ is called the third elasticity tensor, i.e.,
          ${\dot {\boldsymbol {P}}}\cdot {\boldsymbol {F}}^{T}={\boldsymbol {\mathsf {a}}}:{\boldsymbol {l}}$
   Therefore, the relation between the first and third elasticity tensors is
       ${\dot {\boldsymbol {P}}}\cdot {\boldsymbol {F}}^{T}=({\boldsymbol {\mathsf {A}}}:{\dot {\boldsymbol {F}}})\cdot {\boldsymbol {F}}^{T}={\boldsymbol {\mathsf {a}}}:{\boldsymbol {l}}={\boldsymbol {\mathsf {a}}}:({\dot {\boldsymbol {F}}}\cdot {\boldsymbol {F}}^{-1})$
   or,
       ${\boldsymbol {\mathsf {A}}}:{\dot {\boldsymbol {F}}}=[{\boldsymbol {\mathsf {a}}}:({\dot {\boldsymbol {F}}}\cdot {\boldsymbol {F}}^{-1})]\cdot {\boldsymbol {F}}^{-T}$
   In index notation
       ${\mathsf {A}}_{iJkL}~{\dot {F}}_{kL}={\mathsf {a}}_{ijkl}~{\dot {F}}_{kL}~F_{Ll}^{-1}~F_{jJ}^{-1}$
   Therefore,
       ${\mathsf {A}}_{iJkL}={\mathsf {a}}_{ijkl}~F_{Ll}^{-1}~F_{jJ}^{-1}=(\delta _{ik}~\tau _{jl}+{\mathsf {c}}_{ijkl})~F_{Ll}^{-1}~F_{jJ}^{-1}\qquad \square$

## An isotropic spatial elasticity tensor?

An isotropic spatial elasticity tensor cannot be derived from a stored energy function if the constitutive relation is of the form
    ${\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]={\boldsymbol {\mathsf {c}}}:{\boldsymbol {d}}$
where
    ${\boldsymbol {\mathsf {c}}}=\lambda ~{\boldsymbol {\mathit {1}}}\otimes {\boldsymbol {\mathit {1}}}+2~\mu ~{\boldsymbol {\mathsf {I}}}^{s}~.$
Since a significant number of finite element codes use such a constitutive equation, (also called the equation of a hypoelastic material of grade 0) it is worth examining why such a model is incompatible with elasticity.

### Start with a constant and isotropic material elasticity tensor

Let us start of with an isotropic elastic material model in the reference configuration. The simplest such model is the St. Venant-Kirchhoff hyperelastic model
    ${\boldsymbol {S}}=\lambda ~{\text{tr}}({\boldsymbol {E}})~{\boldsymbol {\mathit {1}}}+2~\mu ~{\boldsymbol {E}}$
where ${\boldsymbol {S}}$ is the second Piola-Kirchhoff stress, ${\boldsymbol {E}}$ is the Lagrangian Green strain, and $\lambda ,\mu$ are material constants. We can show that this equation can be derived from a stored energy function.
Taking the material time derivative of this equation, we get
    ${\dot {\boldsymbol {S}}}=\lambda ~{\text{tr}}({\dot {\boldsymbol {E}}})~{\boldsymbol {\mathit {1}}}+2~\mu ~{\dot {\boldsymbol {E}}}$
Now,
    ${\dot {\boldsymbol {S}}}={\boldsymbol {\mathsf {C}}}:{\frac {1}{2}}~{\dot {\boldsymbol {C}}}={\boldsymbol {\mathsf {C}}}:{\dot {\boldsymbol {E}}}$
where ${\boldsymbol {\mathsf {C}}}$ is the second (material) elasticity tensor.
Therefore,
    ${\boldsymbol {\mathsf {C}}}:{\dot {\boldsymbol {E}}}=\lambda ~{\text{tr}}({\dot {\boldsymbol {E}}})~{\boldsymbol {\mathit {1}}}+2~\mu ~{\dot {\boldsymbol {E}}}$
which implies that
    ${\boldsymbol {\mathsf {C}}}=\lambda ~{\boldsymbol {\mathit {1}}}\otimes {\boldsymbol {\mathit {1}}}+2~\mu ~{\boldsymbol {\mathsf {I}}}^{s}~.$
In index notation,
    ${\mathsf {C}}_{IJKL}=\lambda ~\delta _{IJ}~\delta _{KL}+\mu ~(\delta _{IK}~\delta _{JL}+\delta _{JK}~\delta _{IL})$
Now, from the relations between the second elasticity tensor and the fourth (spatial) elasticity tensor, we have
    ${\mathsf {c}}_{ijkl}=F_{iI}~F_{jJ}~F_{kK}~F_{lL}~{\mathsf {C}}_{IJKL}$
Therefore, in this case,
    ${\mathsf {c}}_{ijkl}=F_{iI}~F_{jJ}~F_{kK}~F_{lL}~\left[\lambda ~\delta _{IJ}~\delta _{KL}+\mu ~(\delta _{IK}~\delta _{JL}+\delta _{JK}~\delta _{IL})\right]$
or,
    ${{\mathsf {c}}_{ijkl}=\lambda ~b_{ij}~b_{kl}+\mu ~(b_{ik}~b_{jl}+b_{il}~b_{jk})}$
where ${\boldsymbol {b}}={\boldsymbol {F}}\cdot {\boldsymbol {F}}^{T}$. So we see that the spatial elasticity tensor ${\boldsymbol {\mathsf {c}}}$ cannot be a constant tensor unless ${\boldsymbol {b}}={\boldsymbol {\mathit {1}}}$.
Alternatively, if we define
    ${\mathsf {c}}_{ijkl}=J^{-1}~F_{iI}~F_{jJ}~F_{kK}~F_{lL}~{\mathsf {C}}_{IJKL}$
we get
    ${{\mathsf {c}}_{ijkl}=J^{-1}~\lambda ~b_{ij}~b_{kl}+J^{-1}~\mu ~(b_{ik}~b_{jl}+b_{il}~b_{jk})}$

### Start with a constant and isotropic spatial elasticity tensor

Let us now look at the situation where we start off with a constant and isotropic spatial elasticity tensor, i.e.,
    ${\boldsymbol {\mathsf {c}}}=\lambda ~{\boldsymbol {\mathit {1}}}\otimes {\boldsymbol {\mathit {1}}}+2~\mu ~{\boldsymbol {\mathsf {I}}}^{s}$
In index notation,
    ${\mathsf {c}}_{ijkl}=\lambda ~\delta _{ij}~\delta _{kl}+\mu (\delta _{ik}~\delta _{jl}+\delta _{jk}~\delta _{il})$
Since
    ${\mathsf {c}}_{ijkl}=F_{iI}~F_{jJ}~F_{kK}~F_{lL}~{\mathsf {C}}_{IJKL}$
multiplying both sides by $F_{Ai}^{-1}~F_{Bj}^{-1}~F_{Ck}^{-1}~F_{Dl}^{-1}$ we have,
    $F_{Ai}^{-1}~F_{Bj}^{-1}~F_{Ck}^{-1}~F_{Dl}^{-1}~{\mathsf {c}}_{ijkl}={\mathsf {C}}_{ABCD}$
Therefore, substituting in the expression for a constant and isotropic ${\boldsymbol {\mathsf {c}}}$, we have
    ${\mathsf {C}}_{ABCD}=F_{Ai}^{-1}~F_{Bj}^{-1}~F_{Ck}^{-1}~F_{Dl}^{-1}~\left[\lambda ~\delta _{ij}~\delta _{kl}+\mu (\delta _{ik}~\delta _{jl}+\delta _{jk}~\delta _{il})\right]$
or,
    ${\mathsf {C}}_{ABCD}=\lambda ~F_{Aj}^{-1}~F_{Bj}^{-1}~F_{Cl}^{-1}~F_{Dl}^{-1}+\mu \left(F_{Ak}^{-1}~F_{Bl}^{-1}~F_{Ck}^{-1}~F_{Dl}^{-1}+F_{Al}^{-1}~F_{Bk}^{-1}~F_{Ck}^{-1}~F_{Dl}^{-1}\right)$
or,
    ${\mathsf {C}}_{ABCD}=\lambda ~({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {F}}^{-T})_{AB}~({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {F}}^{-T})_{CD}+\mu \left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {F}}^{-T})_{AC}~({\boldsymbol {F}}^{-1}\cdot ~{\boldsymbol {F}}^{-T})_{BD}+{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {F}}^{-T})_{AD}~({\boldsymbol {F}}^{-1}\cdot ~{\boldsymbol {F}}^{-T})_{BC}\right)$
Since ${\boldsymbol {C}}={\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}$ which gives us ${\boldsymbol {C}}^{-1}={\boldsymbol {F}}^{-1}\cdot {\boldsymbol {F}}^{-T}$, we can write
    ${{\mathsf {C}}_{ABCD}=\lambda ~{\boldsymbol {C}}_{AB}^{-1}~{\boldsymbol {C}}_{CD}^{-1}+\mu \left({\boldsymbol {C}}_{AC}^{-1}~{\boldsymbol {C}}_{BD}^{-1}+{\boldsymbol {C}}_{AD}^{-1}~{\boldsymbol {C}}_{BC}^{-1}\right)}$
Alternatively, if we define
    ${\mathsf {c}}_{ijkl}=J^{-1}~F_{iI}~F_{jJ}~F_{kK}~F_{lL}~{\mathsf {C}}_{IJKL}$
we get
    ${\mathsf {C}}_{ABCD}=J~F_{Ai}^{-1}~F_{Bj}^{-1}~F_{Ck}^{-1}~F_{Dl}^{-1}~{\mathsf {c}}_{ijkl}$
and therefore,
    ${{\mathsf {C}}_{ABCD}=J~\lambda ~{\boldsymbol {C}}_{AB}^{-1}~{\boldsymbol {C}}_{CD}^{-1}+J~\mu \left({\boldsymbol {C}}_{AC}^{-1}~{\boldsymbol {C}}_{BD}^{-1}+{\boldsymbol {C}}_{AD}^{-1}~{\boldsymbol {C}}_{BC}^{-1}\right)}$

## Hypoelastic material of grade 0

A hypoelastic material of grade zero is one for which the stress-strain relation in rate form can be expressed as
    ${\mathcal {L}}_{\varphi }[{\boldsymbol {\sigma }}]={\boldsymbol {\mathsf {c}}}:{\boldsymbol {d}}$
where ${\boldsymbol {\mathsf {c}}}$ is constant. When the material is isotropic we have
    ${\boldsymbol {\mathsf {c}}}=\lambda ~{\boldsymbol {\mathit {1}}}\otimes {\boldsymbol {\mathit {1}}}+2~\mu ~{\boldsymbol {\mathsf {I}}}^{s}$
We want to show that hypoelastic material models of grade 0 cannot be derived from a stored energy function. To do that, recall that
    ${\mathsf {C}}_{IJKL}=4~\rho _{0}~{\frac {\partial ^{2}W}{\partial C_{IJ}\partial C_{KL}}}$
and
    $S_{IJ}=2~\rho _{0}~{\frac {\partial W}{\partial C_{IJ}}}$
For a material elasticity tensor ${\boldsymbol {\mathsf {C}}}$ to be derivable from a stored energy function it has to satisfy the Bernstein integrability conditions. We have
    ${\frac {\partial S_{IJ}}{\partial C_{KL}}}=2~\rho _{0}~{\frac {\partial ^{2}W}{\partial C_{IJ}\partial C_{KL}}}={\frac {1}{2}}~{\mathsf {C}}_{IJKL}$
Also, due to the interchangeability of derivatives,
    ${\frac {\partial }{\partial C_{MN}}}\left({\frac {\partial S_{IJ}}{\partial C_{KL}}}\right)={\frac {\partial }{\partial C_{KL}}}\left({\frac {\partial S_{IJ}}{\partial C_{MN}}}\right)$
Therefore,
    ${{\frac {\partial {\mathsf {C}}_{IJKL}}{\partial C_{MN}}}={\frac {\partial {\mathsf {C}}_{IJMN}}{\partial C_{KL}}}}$
These integrability conditions have to be satisfied by any material elasticity tensor.
At this stage we will use the relation
    ${\mathsf {C}}_{ABCD}=J~\lambda ~{\boldsymbol {C}}_{AB}^{-1}~{\boldsymbol {C}}_{CD}^{-1}+J~\mu \left({\boldsymbol {C}}_{AC}^{-1}~{\boldsymbol {C}}_{BD}^{-1}+{\boldsymbol {C}}_{AD}^{-1}~{\boldsymbol {C}}_{BC}^{-1}\right)$
If we plug this into the integrability condition we will see that
    ${(\lambda +\mu )~C_{KL}^{-1}(C_{IM}^{-1}~C_{JN}^{-1}+C_{IN}^{-1}~C_{JM}^{-1})=(\lambda +\mu )~C_{MN}^{-1}(C_{IK}^{-1}~C_{JL}^{-1}+C_{IL}^{-1}~C_{JK}^{-1})}$
If we multiply both sides by $C_{MN}~C_{IK}~C_{JL}$ we are left with
    ${\lambda +\mu =0}$
This is an unphysical situation and hence shows that a hypoelastic material of grade zero requires that $\lambda =-\mu$ for it to be derivable from a stored energy function.
Proof:
   Let us simplify the notation by writing ${\boldsymbol {A}}:={\boldsymbol {C}}^{-1}$. Then,
       ${\mathsf {C}}_{IJKL}=J~\lambda ~A_{IJ}~A_{KL}+J~\mu \left(A_{IK}~A_{JL}+A_{JK}~A_{IL}\right)$
   Then,
       ${\begin{aligned}{\frac {\partial {\mathsf {C}}_{IJKL}}{\partial C_{MN}}}=&\lambda ~{\frac {\partial J}{\partial C_{MN}}}~A_{IJ}~A_{KL}+\\&J~\lambda ~\left[{\frac {\partial A_{IJ}}{\partial C_{MN}}}~A_{KL}+A_{IJ}~{\frac {\partial A_{KL}}{\partial C_{MN}}}\right]+\\&\mu ~{\frac {\partial J}{\partial C_{MN}}}~(A_{IK}~A_{JL}+A_{JK}~A_{IL})+\\&J~\mu ~\left[{\frac {\partial A_{IK}}{\partial C_{MN}}}~A_{JL}+A_{IK}~{\frac {\partial A_{JL}}{\partial C_{MN}}}+{\frac {\partial A_{JK}}{\partial C_{MN}}}~A_{IL}+A_{JK}~{\frac {\partial A_{IL}}{\partial C_{MN}}}\right]\end{aligned}}$
   and
       ${\begin{aligned}{\frac {\partial {\mathsf {C}}_{IJMN}}{\partial C_{KL}}}=&\lambda ~{\frac {\partial J}{\partial C_{KL}}}~A_{IJ}~A_{MN}+\\&J~\lambda ~\left[{\frac {\partial A_{IJ}}{\partial C_{KL}}}~A_{MN}+A_{IJ}~{\frac {\partial A_{MN}}{\partial C_{KL}}}\right]+\\&\mu ~{\frac {\partial J}{\partial C_{KL}}}~(A_{IM}~A_{JN}+A_{JM}~A_{IN})+\\&J~\mu ~\left[{\frac {\partial A_{IM}}{\partial C_{KL}}}~A_{JN}+A_{IM}~{\frac {\partial A_{JN}}{\partial C_{KL}}}+{\frac {\partial A_{JM}}{\partial C_{KL}}}~A_{IN}+A_{JM}~{\frac {\partial A_{IN}}{\partial C_{KL}}}\right]\end{aligned}}$
   As this stage we use the identities (see Nonlinear finite elements/Kinematics#Some_useful_results for proofs)
       ${\frac {\partial J}{\partial C_{IJ}}}={\cfrac {J}{2}}~C_{IJ}^{-1}={\cfrac {J}{2}}~A_{IJ}$
   and
       ${\frac {\partial C_{IJ}^{-1}}{\partial C_{KL}}}=-{\frac {1}{2}}~(C_{IK}^{-1}~C_{JL}^{-1}+C_{JK}^{-1}~C_{IL}^{-1})\qquad {\text{or}}\qquad {\frac {\partial A_{IJ}}{\partial C_{KL}}}=-{\frac {1}{2}}~(A_{IK}~A_{JL}+A_{JK}~A_{IL})$
   Therefore we have
       ${\begin{aligned}{\frac {\partial {\mathsf {C}}_{IJKL}}{\partial C_{MN}}}=&{\cfrac {J~\lambda }{2}}~A_{MN}~A_{IJ}~A_{KL}-\\&{\cfrac {J~\lambda }{2}}~\left[(A_{IM}~A_{JN}+A_{JM}~A_{IN})~A_{KL}+(A_{KM}~A_{LN}+A_{LM}~A_{KN})~A_{IJ}\right]+\\&{\cfrac {J~\mu }{2}}~A_{MN}~(A_{IK}~A_{JL}+A_{JK}~A_{IL})-\\&{\cfrac {J~\mu }{2}}~\left[(A_{IM}~A_{KN}+A_{KM}~A_{IN})~A_{JL}+(A_{JM}~A_{LN}+A_{LM}~A_{JN})~A_{IK}\right.+\\&\qquad \quad \left.(A_{JM}~A_{KN}+A_{KM}~A_{JN})~A_{IL}+(A_{IM}~A_{LN}+A_{LM}~A_{IN})A_{JK}\right]\end{aligned}}$
   and
       ${\begin{aligned}{\frac {\partial {\mathsf {C}}_{IJMN}}{\partial C_{KL}}}=&{\cfrac {J~\lambda }{2}}~A_{KL}~A_{IJ}~A_{MN}-\\&{\cfrac {J~\lambda }{2}}~\left[(A_{IK}~A_{JL}+A_{JK}~A_{IL})~A_{MN}+(A_{KM}~A_{LN}+A_{LM}~A_{KN})A_{IJ}\right]+\\&{\cfrac {J~\mu }{2}}~A_{KL}~(A_{IM}~A_{JN}+A_{JM}~A_{IN})-\\&{\cfrac {J~\mu }{2}}~\left[(A_{IK}~A_{LM}+A_{KM}~A_{IL})~A_{JN}+(A_{JK}~A_{LN}+A_{JL}~A_{KN})~A_{IM}+\right.\\&\qquad \qquad \left.(A_{JK}~A_{LM}+A_{JL}~A_{KM})~A_{IN}+(A_{IK}~A_{LN}+A_{IL}~A_{KN})~A_{JM}\right]\end{aligned}}$
   Equating the two, we see that the terms that cancel out are
       ${\cfrac {J~\lambda }{2}}~A_{KL}~A_{IJ}~A_{MN}~;~~{\cfrac {J~\lambda }{2}}~(A_{KM}~A_{LN}+A_{LM}~A_{KN})A_{IJ}~;~~$
   and
       ${\begin{aligned}{\cfrac {J~\mu }{2}}~&\left[(A_{IM}~A_{KN}+A_{KM}~A_{IN})~A_{JL}+(A_{JM}~A_{LN}+A_{LM}~A_{JN})~A_{IK}\right.+\\&\left.(A_{JM}~A_{KN}+A_{KM}~A_{JN})~A_{IL}+(A_{IM}~A_{LN}+A_{LM}~A_{IN})A_{JK}\right]\end{aligned}}$
   Therefore,
       ${\frac {\partial {\mathsf {C}}_{IJKL}}{\partial C_{MN}}}={\frac {\partial {\mathsf {C}}_{IJMN}}{\partial C_{KL}}}$
   implies that
       ${\begin{aligned}-{\cfrac {J~\lambda }{2}}~&(A_{IM}~A_{JN}+A_{JM}~A_{IN})~A_{KL}+{\cfrac {J~\mu }{2}}~A_{MN}~(A_{IK}~A_{JL}+A_{JK}~A_{IL})=\\&-{\cfrac {J~\lambda }{2}}~(A_{IK}~A_{JL}+A_{JK}~A_{IL})~A_{MN}+{\cfrac {J~\mu }{2}}~A_{KL}~(A_{IM}~A_{JN}+A_{JM}~A_{IN})\end{aligned}}$
   or,
       $(\lambda +\mu )~A_{MN}~(A_{IK}~A_{JL}+A_{JK}~A_{IL})=(\lambda +\mu )~A_{KL}~(A_{IM}~A_{JN}+A_{JM}~A_{IN})$
   In other words,
       $(\lambda +\mu )~C_{KL}^{-1}(C_{IM}^{-1}~C_{JN}^{-1}+C_{IN}^{-1}~C_{JM}^{-1})=(\lambda +\mu )~C_{MN}^{-1}(C_{IK}^{-1}~C_{JL}^{-1}+C_{IL}^{-1}~C_{JK}^{-1})$
   Now, if we multiply both sides by $C_{MN}$ we get
       $(\lambda +\mu )~C_{KL}^{-1}(\delta _{IN}~C_{JN}^{-1}+\delta _{IM}~C_{JM}^{-1})=(\lambda +\mu )~\delta _{MM}(C_{IK}^{-1}~C_{JL}^{-1}+C_{IL}^{-1}~C_{JK}^{-1})$
   or,
       $2~(\lambda +\mu )~C_{KL}^{-1}~C_{JI}^{-1}=3~(\lambda +\mu )~(C_{IK}^{-1}~C_{JL}^{-1}+C_{IL}^{-1}~C_{JK}^{-1})$
   Next, multiplying both sides by $C_{IK}$ gives
       $2~(\lambda +\mu )~\delta _{IL}~C_{JI}^{-1}=3~(\lambda +\mu )~(\delta _{KK}~C_{JL}^{-1}+\delta _{KL}~C_{JK}^{-1})$
   or,
       $2~(\lambda +\mu )~C_{JL}^{-1}=12~(\lambda +\mu )~C_{JL}^{-1}$
   Finally, multiplying both sides by $C_{JL}$ gives
       $(\lambda +\mu )~\delta _{JJ}=6~(\lambda +\mu )~\delta _{JJ}$
   Therefore,
       $\lambda +\mu =0$