Nonlinear finite elements/Kinematics - spectral decomposition

Spectral decompositions

Many numerical algorithms use spectral decompositions to compute material behavior.

Spectral decompositions of stretch tensors

Infinitesimal line segments in the material and spatial configurations are related by

d\mathbf {x} ={\boldsymbol {F}}\cdot d{\boldsymbol {X}}={\boldsymbol {R}}\cdot ({\boldsymbol {U}}\cdot d{\boldsymbol {X}})={\boldsymbol {V}}\cdot ({\boldsymbol {R}}\cdot d{\boldsymbol {X}})~.

So the sequence of operations may be either considered as a stretch of in the material configuration followed by a rotation or a rotation followed by a stretch.

Also note that

{\boldsymbol {V}}={\boldsymbol {R}}\cdot {\boldsymbol {U}}\cdot {\boldsymbol {R}}^{T}~.

Let the spectral decomposition of ${\boldsymbol {U}}$ be

{\boldsymbol {U}}=\sum _{i=1}^{3}\lambda _{i}~{\boldsymbol {N}}_{i}\otimes {\boldsymbol {N}}_{i}

and the spectral decomposition of ${\boldsymbol {V}}$ be

{\boldsymbol {V}}=\sum _{i=1}^{3}{\hat {\lambda }}_{i}~\mathbf {n} _{i}\otimes \mathbf {n} _{i}~.

Then

{\boldsymbol {V}}={\boldsymbol {R}}\cdot {\boldsymbol {U}}\cdot {\boldsymbol {R}}^{T}=\sum _{i=1}^{3}\lambda _{i}~{\boldsymbol {R}}\cdot ({\boldsymbol {N}}_{i}\otimes {\boldsymbol {N}}_{i})\cdot {\boldsymbol {R}}^{T}=\sum _{i=1}^{3}\lambda _{i}~({\boldsymbol {R}}\cdot {\boldsymbol {N}}_{i})\otimes ({\boldsymbol {R}}\cdot {\boldsymbol {N}}_{i})

Therefore the uniqueness of the spectral decomposition implies that

\lambda _{i}={\hat {\lambda }}_{i}\quad {\text{and}}\quad \mathbf {n} _{i}={\boldsymbol {R}}\cdot {\boldsymbol {N}}_{i}

The left stretch ( ${\boldsymbol {V}}$ ) is also called the spatial stretch tensor while the right stretch ( ${\boldsymbol {U}}$ ) is called the material stretch tensor.

Spectral decompositions of deformation gradient

The deformation gradient is given by

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

In terms of the spectral decomposition of ${\boldsymbol {U}}$ we have

{\boldsymbol {F}}=\sum _{i=1}^{3}\lambda _{i}~{\boldsymbol {R}}\cdot ({\boldsymbol {N}}_{i}\otimes {\boldsymbol {N}}_{i})=\sum _{i=1}^{3}\lambda _{i}~({\boldsymbol {R}}\cdot {\boldsymbol {N}}_{i})\otimes {\boldsymbol {N}}_{i}=\sum _{i=1}^{3}\lambda _{i}~\mathbf {n} _{i}\otimes {\boldsymbol {N}}_{i}

Therefore the spectral decomposition of ${\boldsymbol {F}}$ can be written as

{\boldsymbol {F}}=\sum _{i=1}^{3}\lambda _{i}~\mathbf {n} _{i}\otimes {\boldsymbol {N}}_{i}

Let us now see what effect the deformation gradient has when it is applied to the eigenvector ${\boldsymbol {N}}_{i}$ .

We have

{\boldsymbol {F}}\cdot {\boldsymbol {N}}_{i}={\boldsymbol {R}}\cdot {\boldsymbol {U}}\cdot {\boldsymbol {N}}_{i}={\boldsymbol {R}}\cdot \left(\sum _{j=1}^{3}\lambda _{j}~{\boldsymbol {N}}_{j}\otimes {\boldsymbol {N}}_{j}\right)\cdot {\boldsymbol {N}}_{i}

From the definition of the dyadic product

(\mathbf {u} \otimes \mathbf {v} )\cdot \mathbf {w} =(\mathbf {w} \cdot \mathbf {v} )~\mathbf {u}

Since the eigenvectors are orthonormal, we have

({\boldsymbol {N}}_{j}\otimes {\boldsymbol {N}}_{j})\cdot {\boldsymbol {N}}_{i}={\begin{cases}0&{\mbox{if}}~i\neq j\\{\boldsymbol {N}}_{i}&{\mbox{if}}~i=j\end{cases}}

Therefore,

\left(\sum _{j=1}^{3}\lambda _{j}~{\boldsymbol {N}}_{j}\otimes {\boldsymbol {N}}_{j}\right)\cdot {\boldsymbol {N}}_{i}=\lambda _{i}~{\boldsymbol {N}}_{i}{\text{no sum on}}~i

That leads to

{\boldsymbol {F}}\cdot {\boldsymbol {N}}_{i}=\lambda _{i}~({\boldsymbol {R}}\cdot {\boldsymbol {N}}_{i})=\lambda _{i}~\mathbf {n} _{i}

So the effect of ${\boldsymbol {F}}$ on ${\boldsymbol {N}}_{i}$ is to stretch the vector by $\lambda _{i}$ and to rotate it to the new orientation $\mathbf {n} _{i}$ .