Overview of Elasticity of Materials/Example Problems

Problem Statement:
Given a point on an elastic body with stress state $\boldsymbol{\sigma}$, determine the principal stresses and the angle between the original orientation and the orientation with maximum shear stress. Solve this problem both using the analytical equations and using Mohr's circle construction. Stresses are given in MPa.

\boldsymbol{\sigma}= \begin{pmatrix} 72 & 21 & 0 \\ 21 & 32 & 0 \\ 0 & 0  & 0 \end{pmatrix} $$

Solution:
We're going to solve this by two methods. First using the analytic solution then using Mohr's circle.

Analytic Solution
We have the equations
 * $$\begin{align}

\sigma'_{11}=\frac{\sigma_{11}+\sigma_{22}}{2}+\frac{\sigma_{11}-\sigma_{22}}{2}cos\left ( 2\theta \right )+\sigma_{12}sin\left ( 2\theta \right ) \\

\sigma'_{22}=\frac{\sigma_{11}+\sigma_{22}}{2}-\frac{\sigma_{11}-\sigma_{22}}{2}cos\left ( 2\theta \right )-\sigma_{12}sin\left ( 2\theta \right ) \end{align}$$ and
 * $$\begin{align}

\sigma'_{12}=\frac{\sigma_{22}-\sigma_{11}}{2}sin\left ( 2\theta \right )+\sigma_{12}cos\left ( 2\theta \right ) \end{align}$$

where the prime stresses are in the rotated reference frame. When $$\sigma'_{12}=0$$ the system is in its principal orientation and $$\sigma'_{11}=\sigma_{p1}$$ and $$\sigma'_{22}=\sigma_{p2}$$. Setting $$\sigma'_{12}=0$$ and solving yields

tan\left ( 2\theta \right )=\frac{2\sigma_{12}}{\sigma_{22}-\sigma_{11}} $$ and upon substituting the stresses given in this problem we find $$\theta=23.2$$°. Substituting this value of $$\theta$$ and the given values for $$\boldsymbol{\sigma}$$ yields $$\sigma'_{11}=\sigma_{p1}=81$$ MPa $$\sigma'_{22}=\sigma_{p2}=23$$ MPa and $$\sigma'_{12}=0$$ MPa, which verifies that this is indeed the principal axis.

If we are given a stress state that is the principal orientation then in our rotated reference frame above $$\boldsymbol{\sigma'}$$ the $$\sigma_{12}$$ term is zero. We find the direction with maximum shear by taking the derivative of $$\sigma'_{12}$$, setting it to zero, and solving for $$\theta$$.

\frac{d \sigma'_{12}}{d\theta}=\frac{d}{d\theta}\left(\frac{\sigma_{22}-\sigma_{11}}{2}sin\left ( 2\theta \right )\right)=\left ( \sigma_{22}-\sigma_{11} \right ) cos\left ( 2\theta \right ) $$ and the maximum occurs when $$\theta=45$$°. This means the maximum or minimum occurs at $$23.2 + 45 = 68.2$$°. Substituting this into our equations along with $$\boldsymbol{\sigma}$$ yields $$\sigma'_{11}=\sigma_{p1}=52$$ MPa, $$\sigma'_{22}=\sigma_{p2}=52$$ MPa and $$\sigma'_{12}=-29$$ MPa, the minimum. By the symmetry of the system rotating an additional $$90$$° yields the maximum, $$\sigma'_{11}=\sigma_{p1}=52$$ MPa, $$\sigma'_{22}=\sigma_{p2}=52$$ MPa, and $$\sigma'_{12}=29$$ MPa.

Mohr's Circle Solution
Following the diagram given here,

Radius=AB=BC=\sqrt{BD^2+CD^2} $$ and

tan\left ( 2\theta \right )=\frac{CD}{BD} $$ Given $$CD=21$$, $$OD=72$$, and $$OE=32$$, we concluded that $$AD=40$$ and $$EB=BD=20$$. Therefore, $$Radius=29$$ and $$\theta=23.2$$°.From the principal orientation$$P1=OB+Radius=OE+EB+Radius=81$$ MPa and $$P2=P1-2 Radius=23$$ MPa. The minimum and maximum shear stress is $$\pm Radius=\pm 29$$ MPa and are found by rotating the system $$23.2+45=68.2$$° and $$23.2 + 45 + 90 = 158.2 = -21.8^\circ$$.

Problem Statement:
Given a point on an elastic body with stress state $\boldsymbol{\sigma}$, determine the principal stresses. Stresses are given in MPa. [Hint: There is only one way to draw the Mohr's circle representation. Use this to simplify your work. You'll find that there is almost no math once you draw the picture.]

\boldsymbol{\sigma}= \begin{pmatrix} 0 & 6  & 8 \\ 6 &  0  & 10 \\ 8 &  10 & 0 \end{pmatrix} $$

Problem Statement:
Given a point on an elastic body with stress state $\boldsymbol{\sigma}$, and that rotating about $x_{3}$ yields stress state $\boldsymbol{\sigma'}$ , fully determine both stress states and the unknown parameter $q$. Stresses are given in MPa.

\boldsymbol{\sigma}= \begin{pmatrix} 25q & \sigma_{12} & 0 \\ \sigma_{21} & 5q & 0 \\ 0 & 0  & 0 \end{pmatrix} ~\boldsymbol{\sigma'}= \begin{pmatrix} 50 & 5 & 0 \\ 5 & \sigma'_{22} & 0 \\ 0 & 0  & 0 \end{pmatrix} $$

Solution:
We have the equations
 * $$\begin{align}

\sigma'_{11}=\frac{\sigma_{11}+\sigma_{22}}{2}+\frac{\sigma_{11}-\sigma_{22}}{2}cos\left ( 2\theta \right )+\sigma_{12}sin\left ( 2\theta \right ) \\

\sigma'_{22}=\frac{\sigma_{11}+\sigma_{22}}{2}-\frac{\sigma_{11}-\sigma_{22}}{2}cos\left ( 2\theta \right )-\sigma_{12}sin\left ( 2\theta \right ) \end{align}$$ and
 * $$\begin{align}

\sigma'_{12}=\frac{\sigma_{22}-\sigma_{11}}{2}sin\left ( 2\theta \right )+\sigma_{12}cos\left ( 2\theta \right ) \end{align}$$ where the prime stresses are in the rotated reference frame.

This immediately allows us to substitute for $$\sigma'_{12}$$ and $$\theta$$ to find $$q=-1/2$$. This give $$\sigma_{11}=-12.5$$ MPa and $$\sigma_{22}=-2.5$$. The invariant relation

\sigma_{11}+\sigma_{22}+\sigma_{33}=I_1 $$ allows

\sigma_{11}+\sigma_{22}+\sigma_{33}=\sigma'_{11}+\sigma'_{22}+\sigma'_{33} $$ with substitutions determines $$\sigma'_{22}=-65$$ MPa. At this point all the parameters are determined except $$\sigma_{12}$$ and simply substituting into either the $$\sigma'_{11}$$ or $$\sigma'_{22}$$ allows us to find $$\sigma_{12}=57.5$$ MPa. The resulting stress tensors, in units of MPa, are

\boldsymbol{\sigma}= \begin{pmatrix} -12.5 & 57.5 & 0 \\ 57.5 & -2.5 & 0 \\ 0 & 0  & 0 \end{pmatrix} ~\boldsymbol{\sigma'}= \begin{pmatrix} 50 & 5 & 0 \\ 5 & -65 & 0 \\ 0 & 0  & 0 \end{pmatrix} $$

The Mohr's circle representation of this solution is shown here. Note that the shaded triangles are similar and therefore can be used to simplify the solution if sought graphically.

Problem Statement:
Given a point on an elastic body with stress state $\boldsymbol{\sigma}$, determine the principal stresses and the angle between the original orientation and the orientation with maximum shear stress. Stresses are given in MPa.

\boldsymbol{\sigma}= \begin{pmatrix} 18 & 12 & 9 \\ 12 & 12 & -6 \\ 9 & -6  & 6 \end{pmatrix} $$

Problem Statement:
Both mathematically and with words explain the impact of applying the transformation tensors $\boldsymbol{T_1}$, $\boldsymbol{T_2}$ ,$\boldsymbol{T_3}$ , and $\boldsymbol{T_4}$ on vector $\boldsymbol{a}$. Apply the same transformation to rank two tensor $\boldsymbol{Z}$.

\boldsymbol{T_1}= \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0  & 1 \end{pmatrix} ~ \boldsymbol{T_2}= \begin{pmatrix} -1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0  & -1 \end{pmatrix} ~ \boldsymbol{T_3}= \begin{pmatrix} -1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0  & -1 \end{pmatrix} ~ \boldsymbol{T_4}= \begin{pmatrix} cos\left(\theta\right) & 0 & sin\left(\theta\right) \\ 0 & 1 & 0 \\ -sin\left(\theta\right) & 0  & cos\left(\theta\right) \end{pmatrix} $$



\boldsymbol{a}= \begin{pmatrix} a_1 \\ a_2 \\ a_3 \end{pmatrix} ~ \boldsymbol{Z}= \begin{pmatrix} Z_{11} & Z_{12} & Z_{13} \\ Z_{21} & Z_{22} & Z_{23} \\ Z_{31} & Z_{32} & Z_{33} \end{pmatrix} $$

Solution:
Remembering the definition of tensor multiplication and the implicit summations used in Einstein notation, we know that the rank two $$\boldsymbol{T}$$ tensor acting on the vector $$\boldsymbol{a}$$ will necessarily result in terms that look like

a'_2=a_1T_{11}+a_2T_{12}+a_3T_{13} $$ which is similar to "normal" matrix multiplication as you've seen before. In contrast applying the rank two transformation tensor on the rank two $$\boldsymbol{Z}$$ tensor requires a double sum which results in terms that look similar to

Z'_{12}=T_{11} T_{21} Z_{11} + T_{11} T_{22} Z_{12} + T_{11} T_{23} Z_{13} + T_{12} T_{21} Z_{21} + T_{12} T_{22} Z_{22} + T_{12} T_{23} Z_{23} + T_{13} T_{21} Z_{31} + T_{13} T_{22} Z_{32} + T_{13} T_{23} Z_{33} $$ When working with these type of 9-term sums it is usually advisable to use a software package to simply the work. The results are given here.

Transformation tensor 1 results in

\boldsymbol{a'}= \begin{pmatrix} a_1 \\ a_2 \\ a_3 \end{pmatrix} ~ \boldsymbol{Z'}= \begin{pmatrix} Z_{11} & Z_{12} & Z_{13} \\ Z_{21} & Z_{22} & Z_{23} \\ Z_{31} & Z_{32} & Z_{33} \end{pmatrix} $$ which corresponds to the identity transformation, i.e., leave the tensors unmodified.

Transformation tensor 2 results in

\boldsymbol{a'}= \begin{pmatrix} -a_1 \\ -a_2 \\ -a_3 \end{pmatrix} ~ \boldsymbol{Z'}= \begin{pmatrix} Z_{11} & Z_{12} & Z_{13} \\ Z_{21} & Z_{22} & Z_{23} \\ Z_{31} & Z_{32} & Z_{33} \end{pmatrix} $$ This transformation is the inversion transformation. This can be seen in the behavior of the vector. Interestingly enough it leaves $$\boldsymbol{Z}$$ unmodified.

Transformation tensor 3 results in

\boldsymbol{a'}= \begin{pmatrix} -a_1 \\ a_2 \\ -a_3 \end{pmatrix} ~ \boldsymbol{Z'}= \begin{pmatrix} Z_{11} & -Z_{12} & Z_{13} \\ -Z_{21} & Z_{22} & -Z_{23} \\ Z_{31} & -Z_{32} & Z_{33} \end{pmatrix} $$ This transformation involves mirroring across the $$x_1$$ and $$x_3$$ directions. This is equivalent to rotation around the $$x_2$$-axis by $$180$$°.

Transformation tensor 4 results in $$ \boldsymbol{a'}= \begin{pmatrix} a_1 cos\left(\theta\right) + a_3 sin\left(\theta\right) \\ a_2 \\ a3 cos\left(\theta\right) - a_1 sin\left(\theta\right) \end{pmatrix} ~ \boldsymbol{Z'}= \begin{pmatrix} Z_{11} cos\left(\theta\right)^2 + Z_{13} cos\left(\theta\right) sin\left(\theta\right) + Z_{31} cos\left(\theta\right) sin\left(\theta\right) + Z_{33} sin\left(\theta\right)^2 & Z_{12} cos\left(\theta\right) + Z_{32} sin\left(\theta\right) & Z_{13} cos\left(\theta\right)^2 - Z_{11} cos\left(\theta\right) sin\left(\theta\right) + Z_{33} cos\left(\theta\right) sin\left(\theta\right) - Z_{31} sin\left(\theta\right)^2 \\ Z_{21} cos\left(\theta\right) + Z_{23} sin\left(\theta\right) & Z_{22} & Z_{23} cos\left(\theta\right) - Z_{21} sin\left(\theta\right) \\ Z_{31} cos\left(\theta\right)^2 - Z_{11} cos\left(\theta\right) sin\left(\theta\right) + Z_{33} cos\left(\theta\right) sin\left(\theta\right) - Z_{13} sin\left(\theta\right)^2 & Z_{32} cos\left(\theta\right) - Z_{12} sin\left(\theta\right) & Z_{33} cos\left(\theta\right)^2 - Z_{13} cos\left(\theta\right) sin\left(\theta\right) - Z_{31} cos\left(\theta\right) sin\left(\theta\right) + Z_{11} sin\left(\theta\right)^2 \end{pmatrix} $$ This transformation involves rotation around the $$x_2$$-axis. Substituting $$180$$° for $$\theta$$° yields identical results to transformation tensor 3.

Problem Statement:
Consider the stress state $\boldsymbol{\sigma}$. Write a transformation tensor that rotates the reference frame into the principal directions. Stresses are given in MPa.

\boldsymbol{\sigma}= \begin{pmatrix} 72 & 21 & 0 \\ 21 & 32 & 0 \\ 0 & 0  & 0 \end{pmatrix} $$

Solution:
We know from Example 1 that the solution is to rotate by $$\theta=23.2$$ around the $$x_3$$ axis and we know the form of the rotation transformation matrix from Example 5, therefore the solution is

\boldsymbol{\sigma}= \begin{pmatrix} cos\left(\theta\right) & sin\left(\theta\right) & 0 \\ sin\left(\theta\right) & cos\left(\theta\right) & 0 \\ 0 & 0  & 0 \end{pmatrix} $$ where $$\theta=23.2$$.

Problem Statement:
Given the displacement tensor $\boldsymbol{e}$, identify the rotation tensor and strain tensor.

\boldsymbol{e}= \begin{pmatrix} 0.001 & 0.001 & 0.002 \\ -0.001 & 0.002 & 0.002 \\ 0.001 & 0.003  & 0.002 \end{pmatrix} $$

Solution:
We can break any tensor into a fully symmetric and fully anti-symmetric tensor resulting in an anti-symmetric rotation tensor

\boldsymbol{\omega}= \begin{pmatrix} 0 & \frac{e_{12}-e_{21}}{2} & \frac{e_{13}-e_{31}}{2} \\ \frac{e_{21}-e_{12}}{2} & 0 & \frac{e_{23}-e_{32}}{2} \\ \frac{e_{31}-e_{13}}{2} & \frac{e_{32}-e_{23}}{2}  & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0.001 & 0.0005 \\ -0.001 & 0 & -0.0005 \\ -0.0005 & 0.0005  & 0 \end{pmatrix} $$ and symmetric strain tensor

\boldsymbol{\varepsilon}= \begin{pmatrix} e_{11} & \frac{e_{12}+e_{21}}{2} & \frac{e_{13}+e_{31}}{2} \\ \frac{e_{21}+e_{12}}{2} & e_{22} & \frac{e_{23}+e_{32}}{2} \\ \frac{e_{31}+e_{13}}{2} & \frac{e_{32}+e_{23}}{2}  & e_{33} \end{pmatrix} = \begin{pmatrix} 0.001 & 0 & 0.0015 \\ 0 & 0 & 0.0025 \\ 0.0015 & 0.0025  & 0 \end{pmatrix} $$ resulting ultimately in $$\boldsymbol{e}=\boldsymbol{\omega}+\boldsymbol{\varepsilon}$$.

Problem Statement:
Within linear, homogeneous, isotropic elasticity theory, for a given stress state (units in MPa) determine the strain state given the Poisson ratio is 0.40 and the shear modulus is 50 GPa. Identify the hydrostatic stress, deviatoric stresses, strain dilatation, and strain deviator.

\boldsymbol{\sigma}= \begin{pmatrix} 10 & 2 & 2 \\ 2 & 5 & -3 \\ 2 & -3  & -5 \end{pmatrix} $$

Solution:
We're given the stress, so returning a solution requires determining the strain, which can be determined by

\varepsilon_{ij}=\frac{1+\nu}{E}\sigma_{ij}-\frac{\nu}{E}\sigma_{kk}\delta_{ij} $$ We're given the Poisson ratio, $$\nu$$ and shear modulus $$G$$ but to use this equation we need the elastic modulus $$E$$. This can be found from

G=\frac{E}{2\left(1 - \nu \right)} $$ which can be inverted to

E=2 G \left(1 - \nu \right) $$ Upon substitution this yields $$E=60$$ GPa. The resulting strain tensor is

\boldsymbol{\varepsilon}= \begin{pmatrix} \frac{1+\nu}{E}\sigma_{11}-\frac{\nu}{E}\left(\sigma_{22}+\sigma_{33}\right) & \frac{\sigma_{12}}{2G} & \frac{\sigma_{13}}{2G} \\ \frac{\sigma_{21}}{2G} & \frac{1+\nu}{E}\sigma_{22}-\frac{\nu}{E}\left(\sigma_{11}+\sigma_{33}\right) & \frac{\sigma_{23}}{2G} \\ \frac{\sigma_{31}}{2G} & \frac{\sigma_{32}}{2G}  & \frac{1+\nu}{E}\sigma_{33}-\frac{\nu}{E}\left(\sigma_{11}+\sigma_{22}\right) \end{pmatrix} = \begin{pmatrix} 0.1667 & 0.0467 & 0.0467 \\ 0.0467 & 0.0500 & -0.0700 \\ 0.0467 & -0.0700  & -0.1833 \end{pmatrix}. $$

The hydrostatic stress is

\sigma_m=\frac{\sigma_{kk}}{3}=\frac{\sigma_{11}+\sigma_{22}+\sigma_{33}}{3}=\frac{10+5-5}{3}=3.33 \text{ GPa} $$ The strain dilatation is

\Delta=\left(1+\varepsilon_{11}\right)\left(1+\varepsilon_{22}\right)\left(1+\varepsilon_{33}\right)-1=\left(1+0.1667\right)\left(1+0.0500\right)\left(1+0.1833\right)-1=0.4496 $$ which results in a mean strain of

\varepsilon_m=\frac{\Delta}{3}=\frac{0.4496}{3}=0.1499 $$ (Note that when strain is small $$\Delta \approx \varepsilon_{11}+\varepsilon_{22}+\varepsilon_{33}=0.400$$, which yields $$\varepsilon_m \approx 0.1333$$. Unfortunately, in this example the strain is relatively large.) The deviatoric stress and strain are then determined by subtracting the mean stress and strain from their respective tensors diagonal.

\boldsymbol{\sigma'}= \begin{pmatrix} 10-3.33 & 2 & 2 \\ 2 & 5-3.33 & -3 \\ 2 & -3  & -5-3.33 \end{pmatrix} = \begin{pmatrix} 6.67 & 2 & 2 \\ 2 & 1.67 & -3 \\ 2 & -3  & -8.33 \end{pmatrix} $$ and

\boldsymbol{\varepsilon'}= \begin{pmatrix} 0.1667-0.1499 & 0.0467 & 0.0467 \\ 0.0467 & 0.0500-0.1499 & -0.0700 \\ 0.0467 & -0.0700  & -0.1833-0.1499 \end{pmatrix} = \begin{pmatrix} 0.0168 & 0.0467 & 0.0467 \\ 0.0467 & 0.0999 & -0.0700 \\ 0.0467 & -0.0700  & -0.3332 \end{pmatrix} $$

Problem Statement:
Within linear, homogeneous, isotropic elasticity theory, for a given strain state determine the stress state given the bulk modulus of 100 GPa and Lam$$\acute{e}$$ parameter of 50 GPa. Identify the hydrostatic stress, deviatoric stresses, strain dilatation, and strain deviator.

\boldsymbol{\varepsilon}= \begin{pmatrix} 0.002 & 0.002 & 0.001 \\ 0.002 & 0.002 & 0.001 \\ 0.001 & 0.001  & -0.001 \end{pmatrix} $$

Solution:
We are given the strain, so returning a solution requires determining the stress, which can be determined by

\sigma_{ij}=\frac{E}{1+\nu}\varepsilon_{ij}-\frac{\nu E}{\left(1 + \nu \right)\left(1 - 2 \nu \right)}\varepsilon_{kk}\delta_{ij} $$ where $$\frac{\nu E}{\left(1 + \nu \right)\left(1 - 2 \nu \right)}$$ is the Lamé parameter $$\lambda$$. Since

G=\frac{E}{2\left(1+\nu\right)} $$ and

G=\frac{3\left(K-\lambda\right)}{2} $$ we know

\frac{E}{\left(1+\nu\right)}=3\left(K-\lambda\right) $$ The stress expression becomes

\sigma_{ij}=3\left(K-\lambda\right)\varepsilon_{ij}+\lambda\varepsilon_{kk}\delta_{ij} $$ Substituting and solving yields

\boldsymbol{\sigma}= \begin{pmatrix} 0.45 & 0.3 & 0.15 \\ 0.3 & 0.45 & 0.15 \\ 0.15 & 0.15  & 0.0 \end{pmatrix} $$GPa

The hydrostatic stress is

\sigma_m=\frac{\sigma_{kk}}{3}=\frac{\sigma_{11}+\sigma_{22}+\sigma_{33}}{3}=\frac{0.45+0.45+0.0}{3}=0.30 \text{ GPa} $$ The strain dilatation is

\Delta=\left(1+\varepsilon_{11}\right)\left(1+\varepsilon_{22}\right)\left(1+\varepsilon_{33}\right)-1=\left(1+0.002\right)\left(1+0.002\right)\left(1-0.001\right)-1=0.003 $$ which results in a mean strain of

\varepsilon_m=\frac{\Delta}{3}=\frac{0.003}{3}=0.001 $$ An interesting observation is that the strain dilation and hydrostatic (mean) stress are related by the bulk modulus $$\sigma_m=K\Delta$$, which in this case is $$0.3=100\times0.003$$. When performing calculations it is important to use known checkpoints such as these to validate your work.

The deviatoric stress and strain are then determined by subtracting the mean stress and strain from their respective tensors diagonal.

\boldsymbol{\sigma'}= \begin{pmatrix} 0.45-0.30 & 0.30 & 0.15 \\ 0.30 & 0.45-0.30 & 0.15 \\ 0.15 & 0.15  & 0.00-0.30 \end{pmatrix} = \begin{pmatrix} 0.15 & 0.30 & 0.15 \\ 0.30 & 0.15 & 0.15 \\ 0.15 & 0.15  & -0.30 \end{pmatrix} $$ and

\boldsymbol{\varepsilon'}= \begin{pmatrix} 0.002-0.001 & 0.002 & 0.001 \\ 0.002 & 0.002-0.001 & 0.001 \\ 0.001 & 0.001  & -0.001-0.001 \end{pmatrix} = \begin{pmatrix} 0.001 & 0.002 & 0.001 \\ 0.002 & 0.001 & 0.001 \\ 0.001 & 0.001  & -0.002 \end{pmatrix} $$

Problem Statement:
For a single crystal of cubic zirconia, which has a elastic constants approximately $$c_{11}=575$$, $$c_{12}=115$$, and $$c_{44}=75$$ GPa, determine the elastic energy required to apply a uniaxial strain of 0.001 in the $$\left[ 100 \right]$$ direction. Determine the elastic energy required to apply a uniaxial strain of 0.001 in the $$\left[ 110 \right]$$ direction. Compute the Zener anisotropy ratio. Using isotropic elasticity theory, the elastic modulus of 200 GPa, and Poisson ratio of 0.3 the elastic modulus of $$200 \ GPa $$ and the Poisson's ratio of $$0.3$$ to compute the elastic energy to apply a uniaxial strain of 0.001.

Solution:
Say that the $$\left[100\right]$$ is in the $$x_1$$ direction. The elastic energy, $$U_o$$, is

U_o=\frac{1}{2}\sigma_{ij}\varepsilon_{ij} $$ so

dU_o=\frac{1}{2}\sigma_{ij}d\varepsilon_{ij} $$ In the case of anisotropic elasticity theory

\sigma_{ij}=c_{ijkl}\varepsilon_{kl}. $$ In this case all $$\varepsilon_{ij}=0$$ except $$\varepsilon_{11}=0.001$$. Therefore all $$\sigma_{ij}=0$$ except $$\sigma_{11}=c_{1111}\varepsilon_{11}$$, $$\sigma_{22}=c_{1122}\varepsilon_{22}$$, and $$\sigma_{33}=c_{1122}\varepsilon_{33}$$.

Substituting into the above equationsː
 * $$\begin{align}

dU_0 &= \varepsilon_{11} c_{1111} d\varepsilon_{11} + \underbrace{\varepsilon_{11} c_{1122} \ d\varepsilon_{22} + \varepsilon_{11} c_{1122} \ d\varepsilon_{33}}_{= \ 0 \ upon \ integration} \\ &= \varepsilon_{11} c_{1111} d\varepsilon_{11} \\ \\ U_0 &= \int_0^{\varepsilon_{11} \ = \ 0.001} \varepsilon_{11} c_{11} d\varepsilon_{11} \\ &= \left. c_{11} {\varepsilon_{11}^2 \over 2} \ \right |_{0}^{0.001} \\ &= (575 *10^9 \ Pa){0.001^2 \over 2} \\ &= 288 \ kJ \end{align}$$

Problem Statement:
For a single crystal of cubic zirconia, which has a elastic constants approximately $$c_{11}=575$$, $$c_{12}=115$$, and $$c_{44}=75$$ GPa, determine the elastic energy required to apply a uniaxial strain of 0.001 in the $$\left[ 100 \right]$$ direction followed by a shear strain of 0.001 on the $$\left( 001 \right)$$ face in the $$\left[ 010 \right]$$ direction.

Solution:
Obeying the given order of operation, first we apply a uniaxial load of $$0.001$$ in the $$[100]$$ direction. (Also covered in Example 10)

Note that our general equations are $dU = \sigma_{ij} \ d\varepsilon_{ij}$, and $\sigma_{ij} = c_{ijkl} \varepsilon_{kl}$.

Restating the strain tensorː

\varepsilon = \left ( \begin{matrix} \varepsilon_{11} & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{matrix} \right )$$

From this we get the non-zero stress tensorː

\sigma = \left ( \begin{matrix} \varepsilon_{11} c_{11} & 0 & 0 \\ 0 & \varepsilon_{11} c_{12} & 0 \\ 0 & 0 & \varepsilon_{11} c_{12} \end{matrix} \right )$$

Thus the only non-zero $$d\varepsilon_{ij}$$ is $$d\varepsilon_{11}$$, and by integrating our equation for $dU$ we getː $$U_1 = 288 \ kJ$$

Now let's apply the shear $$0.001$$ on the $$(100)$$ face in the $$[010]$$ direction.

Our new stress tensor isː

\sigma = \left ( \begin{matrix} \varepsilon_{11} c_{11} & \varepsilon_{12} c_{44} & 0 \\ \varepsilon_{12} c_{44} & \varepsilon_{11} c_{12} & 0 \\ 0 & 0 & \varepsilon_{11} c_{12} \end{matrix} \right )$$

Here, the only non-zero terms are $$d\varepsilon_{12}$$, and $$d\varepsilon_{21}$$. As these are equivalent terms, due to symmetry, we can solve for one and multiply the answer by two.

Once again utilizing our basic energy equation we getː

\begin{align} dU &= \sigma_{ij} \ d\varepsilon_{ij} = \varepsilon_{12} c_{44} \ d\varepsilon_{12} \\ \\ U_2 &= 2\int_0^{\varepsilon_{12}} \varepsilon_{12} c_{44} \ d\varepsilon_{12} \\ &= 2 \left [c_{44} {\varepsilon_{12}^2 \over 2} \right ]_0^{\varepsilon_{12}} \\ &= (75*10^9 \ Pa) 0.001^2 \\ U_2 &= 75 \ kJ \end{align}$$

Finally, adding $$U_1$$ and $$U_2$$ together to get the total energy of this combined transformation gives us a final answer of $$363 \ kJ$$.

Problem Statement:
For a polycrystal specimen of cubic zirconia, which has a elastic constants approximately $$c_{11}=575$$, $$c_{12}=115$$, and $$c_{44}=75$$ GPa, use isotropic elasticity theory and the elastic modulus of $$200 \ GPa $$ and the Poisson's ratio of $$0.3$$ to compute the elastic energy to apply a strain state

\boldsymbol{\varepsilon}= \begin{pmatrix} 0.001 & 0.001 & 0.000 \\ 0.001 & 0.000 & 0.000 \\ 0.000 & 0.000 & 0.000 \end{pmatrix} $$ If the polycrystal material has a porosity of 2% approximately how much will this change the elastic modulus? Approximately how much will this change the elastic energy for this applied strain?

Solution:
Note that our general equations are $dU = \sigma_{ij} \ d\varepsilon_{ij}$, and $\sigma_{ij} = c_{ijkl} \varepsilon_{kl}$.

Here, we can put the stress in terms of the elastic modulus and the Poisson's ratioː


 * $$\sigma_{ij} = {E \over 1 + \nu} \varepsilon_{ij} +

{\nu E \over (1 + \nu) (1 - 2 \nu)} \varepsilon_{kk} \partial_{ij} $$

Where the latter constant is equivalent to the Lamé Constant ($$\lambda $$)ː $$\lambda = {\nu E \over (1 + \nu)(1 - 2\nu)} $$

Briefly solving for the Lamé Constant yields usː

\begin{align} \lambda &= {\nu E \over (1 + \nu)(1 - 2\nu)} \\ &= {0.3 * 200 \ GPa \over (1 + 0.3) (1 - 2(0.3))} = 115.3846 \ GPa \end{align} $$

Keeping in mind that $$\varepsilon_{12}$$ equals $$\varepsilon_{21}$$ thanks to symmetry, the non-zero stresses for this problem areː

\sigma = \left ( \begin{matrix} {E \over 1+ \nu} \varepsilon_{11} + \lambda \varepsilon_{11} & {E \over 1 + \nu} \varepsilon_{12} & 0 \\ {E \over 1 + \nu} \varepsilon_{12} & \lambda \varepsilon_{11} & 0 \\ 0 & 0 & \lambda \varepsilon_{11} \end{matrix} \right )$$

Therefore, we can write out the non-zero energy terms asː


 * $$\begin{align}

dU &= \sigma_{ij} \ d\varepsilon_{ij} \\ &= \left ( {E \over 1 + \nu} +\lambda \right ) \varepsilon_{11} \operatorname{d} \! \varepsilon_{11} + \left ( {E \over 1 + \nu} +\lambda \right ) \varepsilon_{12} \operatorname{d} \! \varepsilon_{12} + \left ( {E \over 1 + \nu} +\lambda \right ) \varepsilon_{21} \operatorname{d} \! \varepsilon_{21} + \lambda \varepsilon_{11} \operatorname{d} \varepsilon_{11} + \lambda \varepsilon_{11} \operatorname{d} \varepsilon_{11} \\ &= \left ( {E \over 1 + \nu} +\lambda \right ) \varepsilon_{11} \operatorname{d} \! \varepsilon_{11} + 2 \left ( {E \over 1 + \nu} +\lambda \right ) \varepsilon_{12} \operatorname{d} \! \varepsilon_{12} \\ \\ U &= \left. \left (    {E \over 1+\nu} + \lambda \right ) {1 \over 2} \varepsilon_{11}^2 \right |_{0}^{\varepsilon_{11}} + \left. \left ( {E \over 1 + \nu} \right ) \varepsilon_{12}^2 \right |_0^{\varepsilon_{12}} \\ &= 288 \ kJ \end{align}$$