is a one-to-one (bijective) mapping [4, 12, 19]
![$$\displaystyle{ \boldsymbol{\varphi }_{t}: \mathcal{B}\rightarrow \mathcal{S}_{t} \in \mathbb{R}^{3}, }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ1.gif)
(1)
![$$\boldsymbol{X} \in \mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq2.gif)
![$$\mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq3.gif)
![$$\displaystyle{ \boldsymbol{x} =\boldsymbol{\varphi } _{t}(\boldsymbol{X}) =\boldsymbol{\varphi } (\boldsymbol{X},t), }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ2.gif)
(2)
![$$\mathcal{S}_{t} \subset \mathbb{R}^{3}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq4.gif)
![$$\displaystyle{ \boldsymbol{F}(\boldsymbol{X},t) = D\boldsymbol{\varphi }(\boldsymbol{X}) = \partial \boldsymbol{x}/\partial \boldsymbol{X}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ3.gif)
(3)
![$$\displaystyle{ J(\boldsymbol{X}) =\det [F(\boldsymbol{X})]> 0. }$$
” src=”/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ4.gif”></DIV></DIV><br />
<DIV class=EquationNumber>(4)</DIV></DIV><br />
<DIV></DIV><br />
<DIV class=Para>Deformation due to swelling or shrinking can be introduced by considering a (fictitious) stress-free intermediate configuration <SPAN id=IEq8 class=InlineEquation><IMG alt=](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq8.gif)
![$$\mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq9.gif)
![$$\mathcal{S}_{t}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq10.gif)
![$$\boldsymbol{F}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq11.gif)
![$$\boldsymbol{F}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq12.gif)
![$$\boldsymbol{F}^{s}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq13.gif)
![$$\displaystyle{ \boldsymbol{F} =\boldsymbol{ F}^{e}\boldsymbol{F}^{s}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ5.gif)
(5)
![A321864_1_En_4_Fig1_HTML.gif](/wp-content/uploads/2016/10/A321864_1_En_4_Fig1_HTML.gif)
Fig. 1
Multiplicative decomposition of the deformation gradient
for the motion of two adjacent material points
to the spatial positions
![$$\boldsymbol{F} =\boldsymbol{ F}^{e}\boldsymbol{F}^{s}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq5.gif)
![$$\boldsymbol{X}_{1},\boldsymbol{X}_{2} \in \mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq6.gif)
![$$\boldsymbol{x}_{1},\boldsymbol{x}_{2} \in \mathcal{S}_{t}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq7.gif)
It should be understood that
and
are not proper gradients and that the intermediate configuration is incompatible in a global sense as indicated by the overlapping neighbourhoods
and
of the points
(Fig. 1) [13, 19]. Conceptually, the multiplicative decomposition can be thought of as the disassembly of a finite element mesh (the initial undeformed configuration) and the subsequent application of the local volumetric deformation due to swelling on the individual elements. Unless the swelling deformation is homogeneous, this intermediate configuration will be incompatible at the boundaries between elements and the elements will no longer “fit together” (more precisely, the intermediate configuration is not a proper configuration because—except for the special case of homogeneous swelling—a bijective mapping between the material particles and
does not exist [6]). Compatibility of the final (deformed) configuration is enforced by reassembling the mesh using the nodal connectivity of the elements. For a formal exposition of the geometrical details, the reader is referred to [6, 12, 13].
![$$\boldsymbol{F}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq14.gif)
![$$\boldsymbol{F}^{s}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq15.gif)
![$$\mathcal{O}_{\boldsymbol{\xi }_{1}}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq16.gif)
![$$\mathcal{O}_{\boldsymbol{\xi }_{2}}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq17.gif)
![$$\boldsymbol{\xi }_{1},\boldsymbol{\xi }_{2} \in \bar{\mathcal{B}}_{st}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq18.gif)
![$$\mathbb{R}^{3}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq19.gif)
For isotropic swelling the swelling deformation gradient can be written as
![$$\displaystyle{ \boldsymbol{F}^{s} =\lambda _{ s}\boldsymbol{I}, }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ6.gif)
where λ s is the isotropic swelling stretch and
is the identity tensor. The elastic deformation gradient can then be expressed simply as
![$$\displaystyle{ \boldsymbol{F}^{e} =\boldsymbol{ F}\left (\boldsymbol{F}^{s}\right )^{-1} = (\lambda _{ s})^{-1}\boldsymbol{F}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ7.gif)
![$$\displaystyle{ \boldsymbol{F}^{s} =\lambda _{ s}\boldsymbol{I}, }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ6.gif)
(6)
![$$\boldsymbol{I}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq20.gif)
![$$\displaystyle{ \boldsymbol{F}^{e} =\boldsymbol{ F}\left (\boldsymbol{F}^{s}\right )^{-1} = (\lambda _{ s})^{-1}\boldsymbol{F}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ7.gif)
(7)
The total Lagrangian (TL) formulation of the finite element method uses rotation invariant measures of strain (Green–Lagrange strain
) and stress (second Piola–Kirchhoff stress
) that are calculated with respect to the initial (undeformed) reference configuration
[3, 4]. The deformations
and
satisfy the invertibility and impenetrability of matter requirements so that the usual push-forward and pull-back operations can be performed to obtain the stress measures with respect to the configurations
,
and
[13]. For convenience, the (effective) elastic second Piola–Kirchhoff stress
is introduced by the push-forward of
onto the relaxed intermediate configuration
scaled by the volume ratio J s . The (total) second Piola–Kirchhoff stress with respect to the initial reference configuration
can then be expressed as
![$$\displaystyle{ \boldsymbol{S} = J^{s}(\boldsymbol{F}^{s})^{-1}\boldsymbol{S}^{e}(\boldsymbol{F}^{s})^{-T}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ8.gif)
![$$\boldsymbol{E}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq21.gif)
![$$\boldsymbol{S}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq22.gif)
![$$\mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq23.gif)
![$$\boldsymbol{F}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq24.gif)
![$$\boldsymbol{F}^{s}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq25.gif)
![$$\mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq26.gif)
![$$\bar{\mathcal{B}}_{st}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq27.gif)
![$$\mathcal{S}_{t}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq28.gif)
![$$\boldsymbol{S}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq29.gif)
![$$\boldsymbol{S}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq30.gif)
![$$\bar{\mathcal{B}}_{st}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq31.gif)
![$$\mathcal{B}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq32.gif)
![$$\displaystyle{ \boldsymbol{S} = J^{s}(\boldsymbol{F}^{s})^{-1}\boldsymbol{S}^{e}(\boldsymbol{F}^{s})^{-T}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ8.gif)
(8)
In practice, when using hyperelastic materials with the TL formulation of the FE method,
is computed using invariants of the elastic deformation gradient
, whereas the nodal forces are usually computed using
and the total deformation gradient
with respect to the initial reference configuration [3, 4]. This enables the swelling behaviour to be included entirely within the material constitutive model. Standard element formulations can then be used to calculate the element nodal forces and displacements in the usual manner. The Cauchy stress can be obtained using the inverse Piola transformation
![$$\displaystyle{ \boldsymbol{\sigma }= J^{-1}\boldsymbol{FSF}^{T}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ9.gif)
![$$\boldsymbol{S}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq33.gif)
![$$\boldsymbol{F}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq34.gif)
![$$\boldsymbol{S}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq35.gif)
![$$\boldsymbol{F}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq36.gif)
![$$\displaystyle{ \boldsymbol{\sigma }= J^{-1}\boldsymbol{FSF}^{T}. }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ9.gif)
(9)
2.2 Constitutive Material Model
At low strain rates the mechanical behaviour of soft tissues can be characterised using a hyperelastic constitutive law [22]. Viscoelastic effects are ignored due to the relatively slow loading speed involved with soft tissue swelling (on the order of a few hours).
Hyperelastic materials are characterised by the existence of a strain energy function
that relates the deformation to the second Piola–Kirchhoff stress
![$$\displaystyle{ \boldsymbol{S}^{e} = 2\frac{\partial W^{e}(\boldsymbol{C}^{e})} {\partial \boldsymbol{C}^{e}}, }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ10.gif)
where
is the elastic right Cauchy–Green deformation tensor [4]. For hyperelastic materials that are isotropic with respect to the initial, unstressed configuration the strain energy density
can be expressed in terms of the principal invariants
,
and
[4]. The classical (incompressible) neo-Hookean model with strain energy density
is based on the assumption that the deformation is isochoric (J e = 1). To account for (slight) compressibility of the material we use the modified strain energy density [24]
![$$W^{e}(\boldsymbol{C}^{e})$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq37.gif)
![$$\displaystyle{ \boldsymbol{S}^{e} = 2\frac{\partial W^{e}(\boldsymbol{C}^{e})} {\partial \boldsymbol{C}^{e}}, }$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_Equ10.gif)
(10)
![$$\boldsymbol{C}^{e} = (\boldsymbol{F}^{e})^{T}\boldsymbol{F}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq38.gif)
![$$W^{e}(I_{1}^{e},I_{2}^{e},I_{3}^{e})$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq39.gif)
![$$I_{1}^{e} =\mathrm{ trace}\ \boldsymbol{C}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq40.gif)
![$$I_{2}^{e} = \tfrac{1} {2}\{(\mathrm{trace}\ \boldsymbol{C}^{e})^{2} -\mathrm{ trace}\ (\boldsymbol{C}^{e})^{2}\}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq41.gif)
![$$I_{3}^{e} =\det \boldsymbol{ C}^{e}$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq42.gif)
![$$W^{e}(I_{1}^{e}) =\mu (I_{1}^{e} - 3)$$](/wp-content/uploads/2016/10/A321864_1_En_4_Chapter_IEq43.gif)
![](https://clinicalpub.com/wp-content/uploads/2023/09/banner1.png)