An important consider in nonlinear finite element analysis is
material frame indifference or objectivity of the material response.
The idea is that the position of the observer frame should not affect the
constitutive relations of a material. You can find more details and a history of
the idea in Truesdell and Noll (1992) - sections 17, 18, 19, and 19A.
[1]
We have already talked about the objectivity of kinematic quantities and stress rates.
Let us now discuss the same ideas with a particular constitutive model in mind.
A detailed description of thermoelastic materials can be found in Continuum mechanics/Thermoelasticity. In this discussion we will avoid the complications induced by
including the temperature.
In the material configuration, a hyperelastic material satisfies two requirements:
- a stored energy function (
) exists for the material.
- the stored energy function depends locally only on the deformation gradient.
Given these requirements, if
is the nominal stress (
is
the first Piola-Kirchhoff stress tensor), then
![{\displaystyle {\boldsymbol {N}}^{T}(\mathbf {X} ,t)={\boldsymbol {P}}(\mathbf {X} ,t)=\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {F}}}}\left[\mathbf {X} ,{\boldsymbol {F}}(\mathbf {X} ,t)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8d1edd1eea6840b8113617de431230fc4e8355e)
The stored energy function
is said to be objective or frame indifferent
if
![{\displaystyle W(\mathbf {X} ,{\boldsymbol {Q}}\cdot {\boldsymbol {F}})=W(\mathbf {X} ,{\boldsymbol {F}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cafcdede7c61bb79c6bddaf36f773e0c3feec98a)
where
is an orthogonal tensor with
.
This objectivity condition can be achieved only if (in the material configuration)
![{\displaystyle W(\mathbf {X} ,{\boldsymbol {F}})={\hat {W}}(\mathbf {X} ,{\boldsymbol {C}})~;~~{\boldsymbol {C}}={\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/120811fb6e94495a17d89ce5774218fd859e21e7)
since
.
We can show that
Constitutive relations for hyperelastic materials
![{\displaystyle {\begin{aligned}{\boldsymbol {N}}&=2~\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\cdot {\boldsymbol {F}}^{T}\\{\boldsymbol {P}}&=2~\rho _{0}~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\\{\boldsymbol {S}}&=2~\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\\{\boldsymbol {\sigma }}&=2~\rho ~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\cdot {\boldsymbol {F}}^{T}\\{\boldsymbol {\tau }}&=2~\rho _{0}~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\cdot {\boldsymbol {F}}^{T}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b290f9a180ae6bb21fdf8b5808655d8f8fd8b46a)
|
Proof:
The stress strain relation for a hyperelastic material is
![{\displaystyle {\boldsymbol {P}}=\rho _{0}~{\frac {\partial W}{\partial {\boldsymbol {F}}}}\left[\mathbf {X} ,{\boldsymbol {F}}(\mathbf {X} ,t)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3c292d13628035f797ff1b69a23c463e90451e35)
The chain rule then implies that
![{\displaystyle {\boldsymbol {P}}:{\boldsymbol {T}}=\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}:\left({\frac {\partial {\boldsymbol {C}}}{\partial {\boldsymbol {F}}}}:{\boldsymbol {T}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/671af172f899d8c4f87882768d3d921326de9e75)
for any second order tensor
.
Now, using the product rule of differentiation,
![{\displaystyle {\frac {\partial {\boldsymbol {C}}}{\partial {\boldsymbol {F}}}}:{\boldsymbol {T}}={\frac {\partial }{\partial {\boldsymbol {F}}}}\left({\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}\right):{\boldsymbol {T}}=\left({\frac {\partial {\boldsymbol {F}}^{T}}{\partial {\boldsymbol {F}}}}:{\boldsymbol {T}}\right)\cdot {\boldsymbol {F}}+{\boldsymbol {F}}^{T}\cdot \left({\frac {\partial {\boldsymbol {F}}}{\partial {\boldsymbol {F}}}}:{\boldsymbol {T}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab6f1e7ed641ffbae8b1497634ceb5ab8cea6747)
or,
![{\displaystyle {\frac {\partial {\boldsymbol {C}}}{\partial {\boldsymbol {F}}}}:{\boldsymbol {T}}=\left({\boldsymbol {\mathsf {I}}}^{T}:{\boldsymbol {T}}\right)\cdot {\boldsymbol {F}}+{\boldsymbol {F}}^{T}\cdot \left({\boldsymbol {\mathsf {I}}}:{\boldsymbol {T}}\right)={\boldsymbol {T}}^{T}\cdot {\boldsymbol {F}}+{\boldsymbol {F}}^{T}\cdot {\boldsymbol {T}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59f91728ba17c3fe1436a5a55a01c656941f7a78)
where
is the fourth order identity tensor. Therefore,
![{\displaystyle {\boldsymbol {P}}:{\boldsymbol {T}}=\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}:\left({\boldsymbol {T}}^{T}\cdot {\boldsymbol {F}}+{\boldsymbol {F}}^{T}\cdot {\boldsymbol {T}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f2317468f800981184508c7a31b0e9d02ce4d761)
Using the identity
![{\displaystyle {\boldsymbol {A}}:({\boldsymbol {B}}^{T}\cdot {\boldsymbol {C}})=({\boldsymbol {C}}\cdot {\boldsymbol {A}}^{T}):{\boldsymbol {B}}=({\boldsymbol {B}}\cdot {\boldsymbol {A}}):{\boldsymbol {C}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8536e1be998e510d27fa5dbce4f41bad84f182da)
we have
![{\displaystyle {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}:({\boldsymbol {T}}^{T}\cdot {\boldsymbol {F}})=\left[{\boldsymbol {F}}\cdot \left({\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\right)^{T}\right]:{\boldsymbol {T}}\quad {\text{and}}\quad {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}:({\boldsymbol {F}}^{T}\cdot {\boldsymbol {T}})=\left[{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\right]:{\boldsymbol {T}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f08a2da1577dc1038a1316b6fa2a8280533dcb3)
Therefore, invoking the arbitrariness of
, we have
![{\displaystyle {\boldsymbol {P}}=\rho _{0}~{\boldsymbol {F}}\cdot \left({\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\right)^{T}+\rho _{0}~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/faef9839eeb0a7199a424a369b9a8368aeb23c94)
Since
we have
![{\displaystyle \left({\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\right)^{T}={\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9f48c0bef3e8f32aea6261ab699feaeb3655439)
which implies that
![{\displaystyle {\boldsymbol {P}}=2~\rho _{0}~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\quad {\text{or}}\quad {\boldsymbol {N}}=2~\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\cdot {\boldsymbol {F}}^{T}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2693f85ab0dc2dbb3767d1f377cc118a97c4dda9)
Recall the relations between the 2nd Piola-Kirchhoff stress tensor and the
first Piola-Kirchhoff stress tensor (and the nominal stress tensor)
![{\displaystyle {\boldsymbol {S}}={\boldsymbol {F}}^{-1}\cdot {\boldsymbol {P}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a4ca42a23ee9f34987b281703aeab140e08341c)
Therefore, we have
![{\displaystyle {\boldsymbol {S}}=2~\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/72f1e89c497654606b351c94868dddda2eff3aa7)
Also from the relation between the Cauchy stress and the 2nd Piola-Kirchhoff stress
tensor
![{\displaystyle {\boldsymbol {\sigma }}=J^{-1}~{\boldsymbol {F}}\cdot {\boldsymbol {S}}\cdot {\boldsymbol {F}}^{T}={\cfrac {\rho }{\rho _{0}}}~{\boldsymbol {F}}\cdot {\boldsymbol {S}}\cdot {\boldsymbol {F}}^{T}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d193bd1ac9b1aa6933bf7461acb34a3ceeb69d8)
we have
![{\displaystyle {\boldsymbol {\sigma }}=2~\rho ~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\cdot {\boldsymbol {F}}^{T}\quad {\text{or}}\quad {\boldsymbol {\tau }}=2~\rho _{0}{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {C}}}}\cdot {\boldsymbol {F}}^{T}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3526d662a2f9b6d42972305294cc796da1443252)
We may also express these relations in terms of the Lagrangian Green strain
![{\displaystyle {\boldsymbol {E}}={\frac {1}{2}}({\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}-{\boldsymbol {\mathit {1}}})={\frac {1}{2}}({\boldsymbol {C}}-{\boldsymbol {\mathit {1}}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1650abda3243db200eaab1217e69e394bd2bfdd7)
Then we have
![{\displaystyle {\frac {\partial W}{\partial {\boldsymbol {C}}}}={\frac {\partial W}{\partial {\boldsymbol {E}}}}:{\frac {\partial {\boldsymbol {E}}}{\partial {\boldsymbol {C}}}}={\frac {1}{2}}~{\frac {\partial W}{\partial {\boldsymbol {E}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/17c54b10c60d9a243c37bede0bc58b2856c992ba)
Hence, we can write
![{\displaystyle {\begin{aligned}{\boldsymbol {N}}&=\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {E}}}}\cdot {\boldsymbol {F}}^{T}\\{\boldsymbol {P}}&=\rho _{0}~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {E}}}}\\{\boldsymbol {S}}&=\rho _{0}~{\frac {\partial {\hat {W}}}{\partial {\boldsymbol {E}}}}\\{\boldsymbol {\sigma }}&=\rho ~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {E}}}}\cdot {\boldsymbol {F}}^{T}\\{\boldsymbol {\tau }}&=\rho _{0}~{\boldsymbol {F}}\cdot {\frac {\partial {\hat {W}}}{\partial {\boldsymbol {E}}}}\cdot {\boldsymbol {F}}^{T}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a41783d45d2ff3938b228fffe5731e5841c4d5d)
|
The stored energy function
is objective if and only if the Cauchy stress tensor is symmetric, i.e., if the balance of angular momentum holds. Show this.
- ↑ C. Truesdell and W. Noll, 1992,
The Nonlinear Field Theories of Mechanics:2nd ed., Springer-Verlag, Berlin