Report errors

Intended learning outcomes

  • Explain how the deformation gradient F\boldsymbol{ F} represents the deformation of a volume element

  • Show that the small strain tensor, ϵ\boldsymbol{ \epsilon}, is a reasonable strain measure for small rotations

Strain

The strain is a local measure of the shape-changes of a structure. As it is a local measure, its connection with the stress is not related to the structure's shape. Therefore, it is suitable to model the material behavior as the stress as a function of the strain: This material behavior can be characterized independently from the structure. For example, we can perform a laboratory test on a test sample to measure the material's stress-strain relationship. We can then apply this relationship to determine the local stress states when simulating an entire structure. But the question then is, how to define such a local strain measure? We have two key requirements

  1. It should describe how an infitesimal volume element (think cube) changes shape

  2. It should be independent of rigid body motion (translation and rotation)

In addition, for convenience, we would like to have a strain measure such that if it is zero there is no deformation

  1. It should be zero if there are no deformations.

The deformation gradient and the small strain tensor

The deformation gradient, F\boldsymbol{ F}, measures the gradient of the deformed coordinates, x=X+u\underline{\boldsymbol{ x}}=\underline{\boldsymbol{ X}}+\underline{\boldsymbol{ u}} (the name displacement gradient is sometimes used for FI=u/X\boldsymbol{ F}-\boldsymbol{ I}=\partial \underline{\boldsymbol{ u}}/\partial \underline{\boldsymbol{ X}})

F=xX=xX=uX+I\begin{aligned} \boldsymbol{ F} = \underline{\boldsymbol{ x}}\otimes\underline{\boldsymbol{ \nabla}}_{\underline{\boldsymbol{ X}}} = \frac{\partial \underline{\boldsymbol{ x}}}{\partial \underline{\boldsymbol{ X}}} = \frac{\partial \underline{\boldsymbol{ u}}}{\partial \underline{\boldsymbol{ X}}} + \boldsymbol{ I} \end{aligned}

Here, x\underline{\boldsymbol{ x}} are the deformed coordinates and X\underline{\boldsymbol{ X}} the undeformed coordinates. That is, if we have the displacements u\underline{\boldsymbol{ u}}, then x=X+u\underline{\boldsymbol{ x}}=\underline{\boldsymbol{ X}}+\underline{\boldsymbol{ u}}. In 2 dimensions, the effect of the deformation gradient on an infitesimal area, which is initially square, can be illustrated as

Based on this figure, it is clear that the deformation gradient describes how an infitesimal volume element changes shape. But it is not zero when the strain is zero - it is the identity tensor. Defining the strain as FI\boldsymbol{ F}-\boldsymbol{ I} solves this problem. However, consider a pure rotation of the square:

For the small rotation, β=5o\beta=5^\mathrm{o}, we see that we obtain a change in the deformation gradient's off-diagonal (shear) terms. However, these are of equal magnitude but with opposite signs. So in the symmetric part of the deformation gradient is negligibly affected by the small rotation. So for small rotations, we can thus define a strain measure as

ϵ=[xX]symI=[xX]symI=[uX]sym\begin{aligned} \boldsymbol{ \epsilon} = \left[ \underline{\boldsymbol{ x}}\otimes\underline{\boldsymbol{ \nabla}}_{\underline{\boldsymbol{ X}}}\right]^{ \mathrm{sym} } - \boldsymbol{ I} = \left[\frac{\partial \underline{\boldsymbol{ x}}}{\partial \underline{\boldsymbol{ X}}}\right]^{ \mathrm{sym} } - \boldsymbol{ I} = \left[\frac{\partial \underline{\boldsymbol{ u}}}{\partial \underline{\boldsymbol{ X}}}\right]^{ \mathrm{sym} } \end{aligned}

This strain measure is the small strain tensor.

Considering the larger rotation, β=30o\beta=30^\mathrm{o}, even the normal terms (on the diagonal) are affected. Therefore, ϵ\boldsymbol{ \epsilon} is not suitable for large rotations.

When considering nonlinear solid mechanics, so-called finite strain measures are introduced. One example is the Green-Lagrange strain, E\boldsymbol{ E}, defined as

E=0.5[FTFI]\begin{aligned} \boldsymbol{ E} = 0.5\left[\boldsymbol{ F}^{\mathrm{T}}\boldsymbol{ F} - \boldsymbol{ I}\right] \end{aligned}

The current coordinates after a pure rotation are given by x=RX\underline{\boldsymbol{ x}} = \boldsymbol{ R}\underline{\boldsymbol{ X}}, where R\boldsymbol{ R} is a proper orthogonal rotation tensor. For this case, we obtain F=R\boldsymbol{ F} = \boldsymbol{ R} and our two strain measures become

ϵ=0.5[R+RT]I0E=0.5[RTRI]=0\begin{aligned} \boldsymbol{ \epsilon} &= 0.5 \left[\boldsymbol{ R} + \boldsymbol{ R}^{\mathrm{T}}\right] - \boldsymbol{ I} \neq \boldsymbol{ 0} \\ \boldsymbol{ E} &= 0.5\left[\boldsymbol{ R}^{\mathrm{T}}\boldsymbol{ R} - \boldsymbol{ I}\right] = \boldsymbol{ 0} \end{aligned}

So why do we not use this strain tensor instead? It involves square terms of the deformation gradient, and is therefore not linear. When solving structural problems, the linearity of the problem makes it much easier to solve. Let's us consider the linearization of the Green-Lagrange strain tensor: We would like to see how it is affected by a small deformation, i.e. F=I+δF\boldsymbol{ F} = \boldsymbol{ I} + \delta\boldsymbol{ F}, where δFI\boldsymbol{ \delta F}\ll\boldsymbol{ I} we have

E=0.5[[IT+δFT][I+δF]I]=0.5[I+δF+δFT+δFTδFI]δFsym=[I+δF]symI=ϵ\begin{aligned} \boldsymbol{ E} &= 0.5\left[\left[\boldsymbol{ I}^{\mathrm{T}}+\boldsymbol{ \delta F}^{\mathrm{T}}\right] \left[\boldsymbol{ I}+\boldsymbol{ \delta F}\right] - \boldsymbol{ I} \right] = 0.5\left[\boldsymbol{ I} + \boldsymbol{ \delta F} + \boldsymbol{ \delta F}^{\mathrm{T}} + \boldsymbol{ \delta F}^{\mathrm{T}}\boldsymbol{ \delta F} - \boldsymbol{ I}\right] \\ &\approx \boldsymbol{ \delta F}^{ \mathrm{sym} } = \left[\boldsymbol{ I} + \boldsymbol{ \delta F}\right]^{ \mathrm{sym} } - \boldsymbol{ I} = \boldsymbol{ \epsilon} \end{aligned}

where the fact that δFI\boldsymbol{ \delta F}\ll \boldsymbol{ I} causes the term δFTδF\boldsymbol{ \delta F}^{\mathrm{T}}\boldsymbol{ \delta F} to be negligible compared to δF\boldsymbol{ \delta F}.