Nonlinear finite elements/Homework 10/Solutions



Problem 1: Kinematics and Stress Rates

Given:

Figure 1 shows a linear three-noded triangular element in the reference configuration.

File:NFE HW10Prob1.png Figure 1. Three-noded triangular element.

The motion of the nodes is given by:

$\begin{array}{cccccc}\text{Node}1:& & x_1(t)& =0& y_1(t)& =0\\ \text{Node}2:& & x_2(t)& =2(1+at)\cos \frac{\pi t}{2}& y_2(t)& =2(1+at)\sin \frac{\pi t}{2}\\ \text{Node}3:& & x_3(t)& =-(1+bt)\sin \frac{\pi t}{2}& y_3(t)& =(1+bt)\cos \frac{\pi t}{2}\end{array}$

The configuration ($x,y$) of the element at time $t$ is given by

$x(\mathit{\xi },t)=\sum _{i=1}^3{x}_i(t)N_i(\mathit{\xi });y(\mathit{\xi },t)=\sum _{i=1}^3{y}_i(t)N_i(\mathit{\xi })$

Solutions

Part 1

In the initial configuration $t=0$. Therefore,

$X_1=x_1(0)=0;Y_1=y_1(0)=0;X_2=x_2(0)=2;Y_2=y_2(0)=0;X_3=x_3(0)=0;Y_3=y_3(0)=1.$

Therefore, the initial configuration is given by

$X(\mathit{\xi })=\sum _{i=1}^3{X}_i{N}_i(\mathit{\xi });Y(\mathit{\xi })=\sum _{i=1}^3{Y}_i{N}_i(\mathit{\xi }).$

Substituting the values of $X_i$ and $Y_i$, we get

$X(\mathit{\xi })=2N_2(\mathit{\xi });Y(\mathit{\xi })=N_3(\mathit{\xi }).$

Hence

$N_2(\mathit{\xi })=\frac{X(\mathit{\xi })}{2};N_3(\mathit{\xi })=Y(\mathit{\xi }).$

Also, the shape functions must satisfy the partition of unity condition

$\sum _{i=1}^3{N}_i(\mathit{\xi })=1.$

Therefore,

$N_1(\mathit{\xi })=1-N_2(\mathit{\xi })-N_3(\mathit{\xi })=1-\frac{X(\mathit{\xi })}{2}-Y(\mathit{\xi }).$

The required expressions are

$N_1(\mathit{\xi })=1-\frac{X(\mathit{\xi })}{2}-Y(\mathit{\xi });N_2(\mathit{\xi })=\frac{X(\mathit{\xi })}{2};N_3(\mathit{\xi })=Y(\mathit{\xi }).$

Part 2

The deformation gradient is given by

$𝐅=\left[\begin{array}{ccccc}\frac{\partial x}{\partial X}& & \frac{\partial x}{\partial Y}& & \frac{\partial x}{\partial Z}\\ \frac{\partial y}{\partial X}& & \frac{\partial y}{\partial Y}& & \frac{\partial y}{\partial Z}\\ \frac{\partial z}{\partial X}& & \frac{\partial z}{\partial Y}& & \frac{\partial z}{\partial Z}\end{array}\right].$

Before computing the derivatives, let us express $x,y,z$ in terms of $X,Y,Z$. Recall

$x(\mathit{\xi },t)=\sum _{i=1}^3{x}_i(t)N_i(\mathit{\xi });y(\mathit{\xi },t)=\sum _{i=1}^3{y}_i(t)N_i(\mathit{\xi })$

Therefore,

$\begin{array}{cc}x(\mathit{\xi },t)& =x_1(t)[1-\frac{X(\mathit{\xi })}{2}-Y(\mathit{\xi })]+x_2(t)\frac{X(\mathit{\xi })}{2}+x_3(t)Y(\mathit{\xi })\\ y(\mathit{\xi },t)& =y_1(t)[1-\frac{X(\mathit{\xi })}{2}-Y(\mathit{\xi })]+y_2(t)\frac{X(\mathit{\xi })}{2}+y_3(t)Y(\mathit{\xi })\\ z(\mathit{\xi },t)& =Z(\mathit{\xi })\end{array}$

Substituting in the expressions for $x_i$ and $y_i$, we get

$\begin{array}{cc}x(\mathit{\xi },t)& =[2(1+at)\cos \frac{\pi t}{2}]\frac{X(\mathit{\xi })}{2}-[(1+bt)\sin \frac{\pi t}{2}]Y(\mathit{\xi })\\ y(\mathit{\xi },t)& =[2(1+at)\sin \frac{\pi t}{2}]\frac{X(\mathit{\xi })}{2}+[(1+bt)\cos \frac{\pi t}{2}]Y(\mathit{\xi })\\ z(\mathit{\xi },t)& =Z(\mathit{\xi }).\end{array}$

In the above expression, the parent coordinates $\mathit{\xi }$ are no longer useful. Therefore, we write

$\begin{array}{cc}x(𝐗,t)& =X(1+at)\cos \frac{\pi t}{2}-Y(1+bt)\sin \frac{\pi t}{2}\\ y(𝐗,t)& =X(1+at)\sin \frac{\pi t}{2}+Y(1+bt)\cos \frac{\pi t}{2}\\ z(𝐗,t)& =Z\end{array}$

where $𝐗=(X,Y,Z)$. Taking derivatives, we get

$\begin{array}{cccccc}\frac{\partial x}{\partial X}& =(1+at)\cos \frac{\pi t}{2}& \frac{\partial x}{\partial Y}& =-(1+bt)\sin \frac{\pi t}{2}& \frac{\partial x}{\partial Z}& =0\\ \frac{\partial y}{\partial X}& =(1+at)\sin \frac{\pi t}{2}& \frac{\partial y}{\partial Y}& =(1+bt)\cos \frac{\pi t}{2}& \frac{\partial y}{\partial Z}& =0\\ \frac{\partial z}{\partial X}& =0& \frac{\partial z}{\partial Y}& =0& \frac{\partial z}{\partial Z}& =1\end{array}$

Therefore,

$𝐅(t)=\left[\begin{array}{ccccc}(1+at)\cos \frac{\pi t}{2}& & -(1+bt)\sin \frac{\pi t}{2}& & 0\\ (1+at)\sin \frac{\pi t}{2}& & (1+bt)\cos \frac{\pi t}{2}& & 0\\ 0& & 0& & 1\end{array}\right].$

The Jacobian determinant is

$J(t)=\det 𝐅(t)=[(1+at)\cos \frac{\pi t}{2}][(1+bt)\cos \frac{\pi t}{2}]+[(1+bt)\sin \frac{\pi t}{2}][(1+at)\sin \frac{\pi t}{2}].$

or

$J(t)=(1+at)(1+bt).$

Part 3

For isochoric motion, $J=1$. Therefore,

$(1+at)(1+bt)=1.$

One possibility is

$a=b=0$

This is a pure rotation.

Another possibility is that

$b=-\frac{a}{1+at}$

This is a combination of shear and rotation where the volume remains constant.

Part 4

We get invalid motions when $J\le 0$. Let us consider the case where $J=0$. Then

$(1+at)(1+bt)=0$

This is possible when

$at=-1\text{or}bt=-1.$

If $J<0$ then

$1+at<0\text{and}1+bt>0\text{or}1+at>0\text{and}1+bt<0$

That is,

$at<-1\text{and}bt>-1\text{or}at>-1\text{and}bt<-1$

Therefore the values at which we get invalid motions are

$a\le -\frac{1}{t}\text{or}b\le -\frac{1}{t}.$

Part 5

Derive the expression for the Green (Lagrangian) strain tensor

for the element as a function of time.

The Green strain tensor is given by

$𝐄=\frac{1}{2}({𝐅}^T𝐅-𝐈)$

Recall

$𝐅(t)=\left[\begin{array}{ccccc}(1+at)\cos \frac{\pi t}{2}& & -(1+bt)\sin \frac{\pi t}{2}& & 0\\ (1+at)\sin \frac{\pi t}{2}& & (1+bt)\cos \frac{\pi t}{2}& & 0\\ 0& & 0& & 1\end{array}\right].$

Let us make the following substitutions

$A:=(1+at);B:=(1+bt);c:=\cos \frac{\pi t}{2};s:=\sin \frac{\pi t}{2}.$

Then

$𝐅=\left[\begin{array}{ccc}Ac& -Bs& 0\\ As& Bc& 0\\ 0& 0& 1\end{array}\right]\text{and}{𝐅}^T=\left[\begin{array}{ccc}Ac& As& 0\\ -Bs& Bc& 0\\ 0& 0& 1\end{array}\right].$

Therefore,

${𝐅}^T𝐅=\left[\begin{array}{ccc}{A}^2& 0& 0\\ 0& B^2& 0\\ 0& 0& 1\end{array}\right]=\left[\begin{array}{ccc}(1+at{)}^2& 0& 0\\ 0& (1+bt{)}^2& 0\\ 0& 0& 1\end{array}\right].$

Hence,

$𝐄=\frac{1}{2}(\left[\begin{array}{ccc}(1+at{)}^2& 0& 0\\ 0& (1+bt{)}^2& 0\\ 0& 0& 1\end{array}\right]-\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right])$

The Green strain is

$𝐄(t)=\frac{1}{2}\left[\begin{array}{ccc}2at+a^2{t}^2& 0& 0\\ 0& 2bt+b^2{t}^2& 0\\ 0& 0& 0\end{array}\right]$

Part 6

Derive an expression for the velocity gradient as a function of

time.

The velocity is the material time derivative of the motion.

Recall that the motion is given by

$\begin{array}{cc}x(𝐗,t)& =X(1+at)\cos \frac{\pi t}{2}-Y(1+bt)\sin \frac{\pi t}{2}\\ y(𝐗,t)& =X(1+at)\sin \frac{\pi t}{2}+Y(1+bt)\cos \frac{\pi t}{2}\\ z(𝐗,t)& =Z\end{array}$

Therefore,

$\begin{array}{cc}{v}_x(𝐗,t)=\frac{\partial x}{\partial t}& =X[a\cos \frac{\pi t}{2}-(1+at)\frac{\pi }{2}\sin \frac{\pi t}{2}]-Y[b\sin \frac{\pi t}{2}+(1+bt)\frac{\pi }{2}\cos \frac{\pi t}{2}]\\ v_y(𝐗,t)=\frac{\partial y}{\partial t}& =X[a\sin \frac{\pi t}{2}+(1+at)\frac{\pi }{2}\cos \frac{\pi t}{2}]+Y[b\cos \frac{\pi t}{2}-(1+bt)\frac{\pi }{2}\sin \frac{\pi t}{2}]\\ v_z𝐗,t)=\frac{\partial z}{\partial t}& =0\end{array}$

We could compute the velocity gradient using

$𝐥=\mathrm{\nabla }𝐯=\left[\begin{array}{ccccc}\frac{\partial v_x}{\partial x}& & \frac{\partial v_x}{\partial y}& & \frac{\partial v_x}{\partial z}\\ \frac{\partial v_y}{\partial x}& & \frac{\partial v_y}{\partial y}& & \frac{\partial v_y}{\partial z}\\ \frac{\partial v_z}{\partial x}& & \frac{\partial v_z}{\partial y}& & \frac{\partial v_z}{\partial z}\end{array}\right]$

after expressing $𝐯$ in terms of $𝐱$. However, that makes the expression quite complicated. Instead, we will use the relation

$𝐥=\dot{𝐅}{𝐅}^{-1}.$

The time derivative of the deformation gradient is

$\dot{𝐅}=\left[\begin{array}{ccccc}a\cos \frac{\pi t}{2}-(1+at)\frac{\pi }{2}\sin \frac{\pi t}{2}& & -b\sin \frac{\pi t}{2}-(1+bt)\frac{\pi }{2}\cos \frac{\pi t}{2}& & 0\\ a\sin \frac{\pi t}{2}+(1+at)\frac{\pi }{2}\cos \frac{\pi t}{2}& & b\cos \frac{\pi t}{2}-(1+bt)\frac{\pi }{2}\sin \frac{\pi t}{2}& & 0\\ 0& & 0& & 0\end{array}\right].$

The inverse of the deformation gradient is

${𝐅}^{-1}=\frac{1}{(1+at)(1+bt)}\left[\begin{array}{ccccc}(1+bt)\cos \frac{\pi t}{2}& & (1+bt)\sin \frac{\pi t}{2}& & 0\\ -(1+at)\sin \frac{\pi t}{2}& & (1+at)\cos \frac{\pi t}{2}& & 0\\ 0& & 0& & (1+at)(1+bt)\end{array}\right].$

Using the substitutions

$A:=(1+at);B:=(1+bt);c:=\cos \frac{\pi t}{2};s:=\sin \frac{\pi t}{2}$

we get

$\dot{𝐅}=\left[\begin{array}{ccccc}ac-As\frac{\pi }{2}& & -bs-Bc\frac{\pi }{2}& & 0\\ as+Ac\frac{\pi }{2}& & bc-Bs\frac{\pi }{2}& & 0\\ 0& & 0& & 0\end{array}\right]=\left[\begin{array}{ccccc}ac& & -bs& & 0\\ as& & bc& & 0\\ 0& & 0& & 0\end{array}\right]+\frac{\pi }{2}\left[\begin{array}{ccccc}-As& & -Bc& & 0\\ Ac& & -Bs& & 0\\ 0& & 0& & 0\end{array}\right]$

and

${𝐅}^{-1}=\frac{1}{AB}\left[\begin{array}{ccccc}Bc& & Bs& & 0\\ -As& & Ac& & 0\\ 0& & 0& & AB\end{array}\right].$

The product is

$𝐥=\frac{1}{AB}\left[\begin{array}{ccccc}Ba{c}^2+Ab{s}^2& & Bacs-Abcs& & 0\\ Bacs-Abcs& & Ba{s}^2+Ab{c}^2& & 0\\ 0& & 0& & 0\end{array}\right]+\frac{\pi }{2}\left[\begin{array}{ccccc}0& & -1& & 0\\ 1& & 0& & 0\\ 0& & 0& & 0\end{array}\right]$

Note that the first matrix is symmetric while the second is skew-symmetric.

Therefore, the velocity gradient is

$𝐥=\left[\begin{array}{ccccc}\frac{a}{1+at}{\cos }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\sin }^2\frac{\pi t}{2}& & [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}-\frac{\pi }{2}& & 0\\ [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}+\frac{\pi }{2}& & \frac{a}{1+at}{\sin }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\cos }^2\frac{\pi t}{2}& & 0\\ 0& & 0& & 0\end{array}\right]$

Part 7

The rate of deformation is the symmetric part of the velocity gradient:

$𝐝=\left[\begin{array}{ccccc}\frac{a}{1+at}{\cos }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\sin }^2\frac{\pi t}{2}& & [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}& & 0\\ [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}& & \frac{a}{1+at}{\sin }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\cos }^2\frac{\pi t}{2}& & 0\\ 0& & 0& & 0\end{array}\right]$

The rate of deformation is the skew-symmetric part of the velocity gradient:

$𝐰=\frac{\pi }{2}\left[\begin{array}{ccccc}0& & -1& & 0\\ 1& & 0& & 0\\ 0& & 0& & 0\end{array}\right]$

Part 8

Assume that $a=1$ and $b=2$. Sketch the undeformed configuration and the deformed configuration at $t=1$ and $t=2$.  Draw both the deformed and undeformed configurations on the same plot

and label.

Recall that in the initial configuration

$X_1=x_1(0)=0;Y_1=y_1(0)=0;X_2=x_2(0)=2;Y_2=y_2(0)=0;X_3=x_3(0)=0;Y_3=y_3(0)=1.$

Also, the motion is

$\begin{array}{cc}x(𝐗,t)& =X(1+at)\cos \frac{\pi t}{2}-Y(1+bt)\sin \frac{\pi t}{2}\\ y(𝐗,t)& =X(1+at)\sin \frac{\pi t}{2}+Y(1+bt)\cos \frac{\pi t}{2}\\ z(𝐗,t)& =Z\end{array}$

Plugging in the values of $a$ and $b$, we get

$\begin{array}{cc}x(𝐗,t)& =X(1+t)\cos \frac{\pi t}{2}-Y(1+2t)\sin \frac{\pi t}{2}\\ y(𝐗,t)& =X(1+t)\sin \frac{\pi t}{2}+Y(1+2t)\cos \frac{\pi t}{2}\\ z(𝐗,t)& =Z\end{array}$

At $t=1$,

$\begin{array}{cc}x(𝐗,1)& =2X\cos \frac{\pi }{2}-3Y\sin \frac{\pi }{2}=-3Y\\ y(𝐗,1)& =2X\sin \frac{\pi }{2}+3Y\cos \frac{\pi }{2}=2X\\ z(𝐗,1)& =Z\end{array}$

At $t=2$,

$\begin{array}{cc}x(𝐗,2)& =3X\cos \pi -5Y\sin \pi =-3X\\ y(𝐗,2)& =3X\sin \pi +5Y\cos \pi =-5Y\\ z(𝐗,2)& =Z\end{array}$

The deformed and undeformed configurations are shown below.

File:NFE HW10Prob1 1.png Deformed and undeformed configurations.

Part 9

The deformation gradient is

$𝐅(t)=\left[\begin{array}{ccccc}(1+t)\cos \frac{\pi t}{2}& & -(1+2t)\sin \frac{\pi t}{2}& & 0\\ (1+t)\sin \frac{\pi t}{2}& & (1+2t)\cos \frac{\pi t}{2}& & 0\\ 0& & 0& & 1\end{array}\right].$

The right Cauchy-Green deformation tensor is

$𝐂={𝐅}^T𝐅=\left[\begin{array}{ccc}(1+t{)}^2& 0& 0\\ 0& (1+2t{)}^2& 0\\ 0& 0& 1\end{array}\right]$

The eigenvalue problem is

$(𝐂-{\lambda }^2𝐈)𝐍=0.$

This problem has a solution if

$\det (𝐂-{\lambda }^2𝐈)=0$

i.e.,

$[{\lambda }^4-(1+t{)}^2{\lambda }^2-{\lambda }^2(1+2t{)}^2+(1+t{)}^2(1+2t{)}^2](1-{\lambda }^2)=0.$

The eigenvalues are (as expected)

${\lambda }_1^2=(1+t{)}^2;{\lambda }_2^2=(1+2t{)}^2;{\lambda }_3^2=1.$

The principal stretches are

${\lambda }_1=(1+t);{\lambda }_2=(1+2t);{\lambda }_3=1.$

The principal directions are (by inspection)

${𝐍}_1=[100{]}^T;{𝐍}_2=[010{]}^T;{𝐍}_3=[001{]}^T.$

Now, the right stretch tensor is given by

$𝑼=\sum _{i=1}^3{\lambda }_i{𝑵}_i\otimes {𝑵}_i$

Therefore,

$𝐔={\lambda }_1{𝐍}_1{𝐍}_1^T+{\lambda }_2{𝐍}_2{𝐍}_2^T+{\lambda }_3{𝐍}_3{𝐍}_3^T.$

Hence the right stretch is

$𝐔=\left[\begin{array}{ccc}1+t& 0& 0\\ 0& 1+2t& 0\\ 0& 0& 1\end{array}\right]$

At $t=0,1,2$, we have

$𝐔(0)=\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right];𝐔(1)=\left[\begin{array}{ccc}2& 0& 0\\ 0& 3& 0\\ 0& 0& 1\end{array}\right];𝐔(2)=\left[\begin{array}{ccc}3& 0& 0\\ 0& 5& 0\\ 0& 0& 1\end{array}\right].$

Now

$𝐑=𝐅{𝐔}^{-1}$

and

${𝐔}^{-1}=\left[\begin{array}{ccc}\frac{1}{1+t}& 0& 0\\ 0& \frac{1}{1+2t}& 0\\ 0& 0& 1\end{array}\right]$

Therefore, the rotation is

$𝐑=\left[\begin{array}{ccccc}\cos \frac{\pi t}{2}& & -\sin \frac{\pi t}{2}& & 0\\ \sin \frac{\pi t}{2}& & \cos \frac{\pi t}{2}& & 0\\ 0& & 0& & 1\end{array}\right].$

At $t=0,1,2$, we have

$𝐑(0)=\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right];𝐑(1)=\left[\begin{array}{ccc}0& -1& 0\\ 1& 0& 0\\ 0& 0& 1\end{array}\right];𝐑(2)=\left[\begin{array}{ccc}-1& 0& 0\\ 0& -1& 0\\ 0& 0& 1\end{array}\right].$

Part 10

A hypoelastic material behaves according to the relation

$\stackrel{\mathrm{△}}{\mathit{\sigma }}=𝖢:𝐝.$

For an isotropic material

$𝖢=\lambda 1\otimes 1+2\mu 𝑰$

Therefore,

$\stackrel{\mathrm{△}}{\mathit{\sigma }}=\lambda 1\otimes 1:𝐝+2\mu 𝑰:𝐝=\lambda \text{tr}(𝐝)1+2\mu 𝐝$

Recall that

$𝐝=\left[\begin{array}{ccccc}\frac{a}{1+at}{\cos }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\sin }^2\frac{\pi t}{2}& & [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}& & 0\\ [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}& & \frac{a}{1+at}a{\sin }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\cos }^2\frac{\pi t}{2}& & 0\\ 0& & 0& & 0\end{array}\right]$

Using the values of $a$ and $b$ from the previous part, at $t=1$,

$𝐝=\left[\begin{array}{ccccc}\frac{1}{2}{\cos }^2\frac{\pi }{2}+\frac{2}{3}{\sin }^2\frac{\pi }{2}& & [\frac{1}{2}-\frac{2}{3}]\cos \frac{\pi }{2}\sin \frac{\pi }{2}& & 0\\ [\frac{1}{2}-\frac{2}{3}]\cos \frac{\pi }{2}\sin \frac{\pi }{2}& & \frac{1}{2}{\sin }^2\frac{\pi }{2}+\frac{2}{3}{\cos }^2\frac{\pi }{2}& & 0\\ 0& & 0& & 0\end{array}\right]=\left[\begin{array}{ccc}\frac{2}{3}& 0& 0\\ 0& \frac{1}{2}& 0\\ 0& 0& 0\end{array}\right]$

Therefore, the trace of the rate of defromation is

$\text{tr}(𝐝)=\frac{7}{6}$

Therefore,

$\stackrel{\mathrm{△}}{\mathit{\sigma }}=\frac{7\lambda }{6}\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right]+2\mu \left[\begin{array}{ccc}\frac{2}{3}& 0& 0\\ 0& \frac{1}{2}& 0\\ 0& 0& 0\end{array}\right]=\left[\begin{array}{ccc}\frac{7\lambda }{6}+\frac{4\mu }{3}& 0& 0\\ 0& \frac{7\lambda }{6}+\mu & 0\\ 0& 0& \frac{7\lambda }{6}\end{array}\right]$

For the Jaumann rate

$\stackrel{\mathrm{△}}{\mathit{\sigma }}=\dot{\mathit{\sigma }}+\mathit{\sigma }\bullet 𝐰-𝐰\bullet \mathit{\sigma }$

where the spin is

$𝐰=\frac{\pi }{2}\left[\begin{array}{ccccc}0& & -1& & 0\\ 1& & 0& & 0\\ 0& & 0& & 0\end{array}\right]$

Therefore,

$\dot{\mathit{\sigma }}=\left[\begin{array}{ccc}\frac{7\lambda }{6}+\frac{4\mu }{3}& 0& 0\\ 0& \frac{7\lambda }{6}+\mu & 0\\ 0& 0& \frac{7\lambda }{6}\end{array}\right]-\frac{\pi }{2}\left[\begin{array}{ccc}{\sigma }_{xx}& {\sigma }_{xy}& {\sigma }_{xz}\\ {\sigma }_{xy}& {\sigma }_{yy}& {\sigma }_{yz}\\ {\sigma }_{xz}& {\sigma }_{yz}& {\sigma }_{zz}\end{array}\right]\left[\begin{array}{ccccc}0& & -1& & 0\\ 1& & 0& & 0\\ 0& & 0& & 0\end{array}\right]+\frac{\pi }{2}\left[\begin{array}{ccccc}0& & -1& & 0\\ 1& & 0& & 0\\ 0& & 0& & 0\end{array}\right]\left[\begin{array}{ccc}{\sigma }_{xx}& {\sigma }_{xy}& {\sigma }_{xz}\\ {\sigma }_{xy}& {\sigma }_{yy}& {\sigma }_{yz}\\ {\sigma }_{xz}& {\sigma }_{yz}& {\sigma }_{zz}\end{array}\right]$

or,

$\dot{\mathit{\sigma }}=\left[\begin{array}{ccc}\frac{7\lambda }{6}+\frac{4\mu }{3}& 0& 0\\ 0& \frac{7\lambda }{6}+\mu & 0\\ 0& 0& \frac{7\lambda }{6}\end{array}\right]+\frac{\pi }{2}\left[\begin{array}{ccc}-{\sigma }_{xy}& {\sigma }_{xx}-{\sigma }_{yy}& -{\sigma }_{yz}\\ {\sigma }_{xx}-{\sigma }_{yy}& 2{\sigma }_{xy}& {\sigma }_{xz}\\ -{\sigma }_{yz}& {\sigma }_{xz}& 0\end{array}\right]$

For the Truesdell rate

$\stackrel{\mathrm{△}}{\mathit{\sigma }}=\stackrel{\circ }{\mathit{\sigma }}=\dot{\mathit{\sigma }}-\mathit{\sigma }\bullet {𝐥}^T-𝐥\bullet \mathit{\sigma }+\text{tr}(𝐥)\mathit{\sigma }.$

Therefore,

$\dot{\mathit{\sigma }}=\stackrel{\circ }{\mathit{\sigma }}+\mathit{\sigma }\bullet {𝐥}^T+𝐥\bullet \mathit{\sigma }-\text{tr}(𝐥)\mathit{\sigma }.$

Recall,

$𝐥=\left[\begin{array}{ccccc}\frac{a}{1+at}{\cos }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\sin }^2\frac{\pi t}{2}& & [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}-\frac{\pi }{2}& & 0\\ [\frac{a}{1+at}-\frac{b}{1+bt}]\cos \frac{\pi t}{2}\sin \frac{\pi t}{2}+\frac{\pi }{2}& & \frac{a}{1+at}{\sin }^2\frac{\pi t}{2}+\frac{b}{1+bt}{\cos }^2\frac{\pi t}{2}& & 0\\ 0& & 0& & 0\end{array}\right]$

For $a=1$, $b=2$, $t=1$, we have

$𝐥=\left[\begin{array}{ccccc}\frac{1}{2}{\cos }^2\frac{\pi }{2}+\frac{2}{3}{\sin }^2\frac{\pi }{2}& & [\frac{1}{2}-\frac{2}{3}]\cos \frac{\pi }{2}\sin \frac{\pi }{2}-\frac{\pi }{2}& & 0\\ [\frac{1}{2}-\frac{2}{3}]\cos \frac{\pi }{2}\sin \frac{\pi }{2}+\frac{\pi }{2}& & \frac{1}{2}{\sin }^2\frac{\pi }{2}+\frac{2}{3}{\cos }^2\frac{\pi }{2}& & 0\\ 0& & 0& & 0\end{array}\right]=\left[\begin{array}{ccc}\frac{2}{3}& -\frac{\pi }{2}& 0\\ \frac{\pi }{2}& \frac{1}{2}& 0\\ 0& 0& 0\end{array}\right].$

Therefore,

$\text{tr}(𝐥)=\frac{7}{6}=\text{tr}(𝐝)$

and

$\dot{\mathit{\sigma }}=\left[\begin{array}{ccc}\frac{7\lambda }{6}+\frac{4\mu }{3}& 0& 0\\ 0& \frac{7\lambda }{6}+\mu & 0\\ 0& 0& \frac{7\lambda }{6}\end{array}\right]+\left[\begin{array}{ccc}{\sigma }_{xx}& {\sigma }_{xy}& {\sigma }_{xz}\\ {\sigma }_{xy}& {\sigma }_{yy}& {\sigma }_{yz}\\ {\sigma }_{xz}& {\sigma }_{yz}& {\sigma }_{zz}\end{array}\right]\left[\begin{array}{ccc}\frac{2}{3}& \frac{\pi }{2}& 0\\ -\frac{\pi }{2}& \frac{1}{2}& 0\\ 0& 0& 0\end{array}\right]+\left[\begin{array}{ccc}\frac{2}{3}& -\frac{\pi }{2}& 0\\ \frac{\pi }{2}& \frac{1}{2}& 0\\ 0& 0& 0\end{array}\right]\left[\begin{array}{ccc}{\sigma }_{xx}& {\sigma }_{xy}& {\sigma }_{xz}\\ {\sigma }_{xy}& {\sigma }_{yy}& {\sigma }_{yz}\\ {\sigma }_{xz}& {\sigma }_{yz}& {\sigma }_{zz}\end{array}\right]-\frac{7}{6}\left[\begin{array}{ccc}{\sigma }_{xx}& {\sigma }_{xy}& {\sigma }_{xz}\\ {\sigma }_{xy}& {\sigma }_{yy}& {\sigma }_{yz}\\ {\sigma }_{xz}& {\sigma }_{yz}& {\sigma }_{zz}\end{array}\right]$

Hence,

$\dot{\mathit{\sigma }}=\left[\begin{array}{ccc}\frac{7\lambda }{6}+\frac{4\mu }{3}& 0& 0\\ 0& \frac{7\lambda }{6}+\mu & 0\\ 0& 0& \frac{7\lambda }{6}\end{array}\right]+\left[\begin{array}{ccc}\frac{1}{6}{\sigma }_{xx}+\pi {\sigma }_{xy}& \frac{\pi }{2}(-{\sigma }_{xx}+{\sigma }_{yy})& \frac{1}{2}(-{\sigma }_{xz}+\pi {\sigma }_{yz})\\ \frac{\pi }{2}(-{\sigma }_{xx}+{\sigma }_{yy})& -\frac{1}{6}{\sigma }_{yy}-\pi {\sigma }_{xy}& -\frac{2}{3}{\sigma }_{yz}-\frac{\pi }{2}{\sigma }_{xz}\\ \frac{1}{2}(-{\sigma }_{xz}+\pi {\sigma }_{yz})& -\frac{2}{3}{\sigma }_{yz}-\frac{\pi }{2}{\sigma }_{xz}& -\frac{7}{6}{\sigma }_{zz}\end{array}\right]$

Problem 2: Hyperelastic Pinched Cylinder Problem

Read the following paper on shells:

Buchter, N., Ramm, E., and Roehl, D., 1994, "Three-dimensional extension of non-linear shell formulation based on the enhanced assumed strain concept," Int. J. Numer. Meth. Engng., 37, pp. 2551-2568.

Answer the following questions:

Solution

Part 1

See Wikipedia article on [w:Enhanced assumed strain|enhanced assumed strain]] .

Part 2

The following material properties are used:

$E=16800$ kN/cm$2^{}$, $\nu =0.4$, $\mu =6000$ kN/cm$2^{}$, and $d=2/K=7.14\times 10^{-5}$ cm$2^{}$/kN. Symmetry is used and only half of the model is meshed. At the support, the model is constraint in all directions. The load of 36 kN is applied on the top of the cylinder (Fig~2. ANSYS input listing is shown Fig~6 of Appendix.

File:NFE HW10 p2 model.png Figure 2. Meshed model

Part 3

The plot of force vs. edge displacement (vertical) and the deformed model are shown in Fig 3. From the plot, one sees that a load of 35.1 kN is required to deform the edge by 16 cm. This is less than 1% difference compared to the result given in the paper using 7-parameter shell theory.

File:NFE HW10 p2 fd.png Figure 3. Left: Force vs. edge displacement. Right: Deformed model.

Problem 3: Elastic-Plastic Punch Indentation

Read the following paper on elastic-viscoplastic FEA:

Rouainia, M. and Peric, D., 1998, "A computational model for elasto-viscoplastic solids at finite strain with reference to thin shell applications," Int. J. Numer. Meth. Engng., 42, pp. 289-311.

Answer the following questions:

Solution

Part 1

The following data are used for the thin plate:$E=70$ GPa, $\nu =0.3$, ${\sigma }_Y=240$ Mpa, ${\sigma }_t=7800$ Mpa, see Fig 4 for the stress-strain data. Symmetry is used and only half of the model is meshed (Fig 4). At the support, the model is constrained in all directions. A plastic finite strained shell (SHELL43) is chosen for this simulation.

The following material properties are used for the punch and die:$E=100\times 10^6$ GPa, $\nu =0.3$. The value of Young's modulus is arbitrary chosen so long as it is high enough to remain rigid during the simulation.

The meshed model is illustrated in Fig 4. A load of 10 kN is applied on the top of the punch. Load steps are split into two steps. First load step the punch is moved close to the plate to activate the contact elements (CONTAC49). The second load step, 30 kN is applied on top of the punch. ANSYS input listing is shown Fig 7 of the Appendix.

File:NFE HW10 p3 model.png Figure 4. Left: Bi-linear stress-strain data. Right: Meshed model.

Part 2

The punch force vs. punch travel and the plot of the deformation at the final load step is shown in Fig 5. The curve displays similar pattern as those shown in the paper.

File:NFE HW10 p3 fd.png Figure 5. Left: Punch force vs. punch travel. Right: Sketch of deformation at final load step.

Appendix

File:NFE HW10 p2 input.pdf Figure 6. ANSYS input listing for Problem 3.

File:NFE HW10 p3 input.pdf Figure 7. ANSYS input listing for Problem 3.