Nonlinear finite elements/Homework 8

Problem 1: Plate Bending
Consider the flat plate shown in Figure 1.

Transverse normals are defined as the straight lines perpendicular to the midsurface. If we assume that the plate is thin, the Kirchhoff-Love hypothesis holds, i.e.,
 * Transverse normals remain straight.
 * Transverse normals remain normal to the midsurface.
 * Transverse normals are inextensible.

Answer the following questions.

Part 1
The displacement field that satisfies the Kirchhoff-Love hypothesis is given by
 * $$\begin{align}

u_1(x_1,x_2,x_3) & = u^0_1(x_1,x_2) - x_3~u^0_{3,1} \\ u_2(x_1,x_2,x_3) & = u^0_2(x_1,x_2) - x_3~u^0_{3,2} \\ u_3(x_1,x_2,x_3) & = u^0_3(x_1,x_2) \end{align}$$ where $$u^0_i$$ are the displacements of the midsurface.

Write the above equations (using the standard partial derivative notation with $$\partial$$~s) after making the following substitutions.
 * $$\text{(2)} \qquad

\begin{align} & x := x_1 ~; y := x_2 ~; z := x_3 \\ & u := u_1 ~; v := u_2 ~; w := u_3 \\ & u_0 := u^0_1 ~; v_0 := u^0_2 ~; w_0 := u^0_3 ~. \end{align} $$

Part 2
Assume that the plate undergoes infinitesimal strains and rotations. Then the strain-displacement relations are

\boldsymbol{\varepsilon} = \frac{1}{2}\left[\boldsymbol{\nabla} \mathbf{u} + (\boldsymbol{\nabla} \mathbf{u})^T\right]~. $$ Show that, in terms of the midsurface displacements, the strain-displacement relations can be written in index notation as
 * $$\begin{align}

\varepsilon_{\alpha\beta} & = \frac{1}{2} \left(u^0_{\alpha,\beta} + u^0_{\beta,\alpha}\right) - x_3~u^0_{3,\alpha\beta} \\ \varepsilon_{\alpha 3} & = 0 \\ \varepsilon_{33} & = 0 \end{align}$$ where $$\alpha = 1,2$$ and $$\beta = 1,2$$.

Part 3
The internal virtual work in the plate is given by (using index notation)

\delta W_{\text{int}} = \int_{\Omega_0}\left(\int_{-h}^{h} \sigma_{\alpha\beta}~\delta \varepsilon_{\alpha\beta}~dx_3\right) d\Omega_0 ~. $$ Note that if a component of the strain tensor is zero, then the variation of that strain component is also zero.

Show that the internal virtual work can be written as
 * $$\text{(3)} \qquad

\delta W_{\text{int}} = \int_{\Omega_0}\left( N_{\alpha\beta}\delta \varepsilon_{\alpha\beta}^0 + M_{\alpha\beta}\delta \varepsilon_{\alpha\beta}^1\right)~d\Omega_0 $$ where
 * $$\begin{align}

N_{\alpha\beta} & = \int_{-h}^h \sigma_{\alpha\beta}~dx_3 \\ M_{\alpha\beta} & = \int_{-h}^h x_3~\sigma_{\alpha\beta}~dx_3 \\ \varepsilon_{\alpha\beta}^0 & = \frac{1}{2} \left(u^0_{\alpha,\beta} + u^0_{\beta,\alpha}\right) \\ \varepsilon_{\alpha\beta}^1 & = - u^0_{3,\alpha\beta} ~. \end{align}$$ Write down the expanded form of equation (3).

Part 4
Reddy's book uses the notation of equation (2) to derive the following  equilibrium equations for the plate.
 * $$\begin{align}

\frac{\partial N_{xx}}{\partial x} + \frac{\partial N_{xy}}{\partial y} & = 0 \\ \frac{\partial N_{xy}}{\partial x} + \frac{\partial N_{yy}}{\partial y} & = 0 \\ \frac{\partial^2 M_{xx}}{\partial x^2} + 2\frac{\partial^2 M_{xy}}{\partial x\partial y} + \frac{\partial^2 M_{yy}}{\partial y^2} + \mathcal{N}(u_0,v_0,w_0) + q & = 0 \end{align}$$ where

\mathcal{N}(u_0,v_0,w_0) = \frac{\partial }{\partial x} \left(N_{xx}\frac{\partial w_0}{\partial x} + N_{xy}\frac{\partial w_0}{\partial y}\right) + \frac{\partial }{\partial y}\left(N_{xy}\frac{\partial w_0}{\partial x} + N_{yy}\frac{\partial w_0}{\partial y} \right) ~. $$ Write the above equations in index notation. Then, show that in index-free  Gibbs notation, the equilibrium equations for the plate can be written as
 * $$\begin{align}

\boldsymbol{\nabla} \bullet \boldsymbol{N} & = \mathbf{0} \text{(4)} \qquad \\ \boldsymbol{\nabla} \bullet (\boldsymbol{\nabla}\bullet\boldsymbol{M}) + \boldsymbol{\nabla} \bullet (\boldsymbol{N}\bullet\boldsymbol{\nabla u^0_3)} + q & = 0 \text{(5)} \qquad \end{align}$$ where

\begin{align} \boldsymbol{N} & = N_{\alpha\beta}~\mathbf{e}_{\alpha}\otimes\mathbf{e}_{\beta} ~;& \boldsymbol{M} & = M_{\alpha\beta}~\mathbf{e}_{\alpha}\otimes\mathbf{e}_{\beta}\\ \boldsymbol{\nabla} \varphi & = \varphi_{,\alpha} ~\mathbf{e}_{\alpha} ~;& \boldsymbol{\nabla} \mathbf{v} & = v_{\alpha,\beta} ~\mathbf{e}_{\alpha}\otimes\mathbf{e}_{\beta} ~;& \boldsymbol{\nabla} \boldsymbol{S} & = S_{\alpha\beta,\gamma} ~\mathbf{e}_{\alpha}\otimes\mathbf{e}_{\beta}\otimes\mathbf{e}_{\gamma} \\ \boldsymbol{\nabla} \bullet \mathbf{v} & = v_{\alpha,\alpha} ~;& \boldsymbol{\nabla} \bullet \boldsymbol{S} & = S_{\alpha\beta,\beta}~\mathbf{e}_{\alpha} ~;& \boldsymbol{\nabla} \bullet \boldsymbol{T} & = T_{\alpha\beta\gamma,\gamma} ~\mathbf{e}_{\alpha}\otimes\mathbf{e}_{\beta} \end{align} $$ and $$\alpha=1,2$$, $$\beta=1,2$$, and $$\gamma=1,2$$.$$\varphi$$ is a scalar, $$\mathbf{v}$$ is a vector, $$\boldsymbol{S}$$ is a second-order tensor, and $$\boldsymbol{T}$$ is a third-order tensor.

Part 5
Derive the symmetric weak forms of the equilibrium equations starting with equations (4) and (5). Use the coordinate-free Gibbs notation in your derivation.

You will find the following identities useful in your derivation
 * $$\begin{align}

\boldsymbol{\nabla} \bullet (\varphi\mathbf{v}) & = \varphi (\boldsymbol{\nabla} \bullet \mathbf{v}) + \mathbf{v}\bullet(\boldsymbol{\nabla} \varphi)\\ \boldsymbol{\nabla} \bullet (\boldsymbol{S}^T\bullet\mathbf{v}) & = \mathbf{v}\bullet(\boldsymbol{\nabla} \bullet \boldsymbol{S}) + \boldsymbol{S}:(\boldsymbol{\nabla} \mathbf{v}) \end{align}$$

Part 6

 * Write down the natural boundary conditions in Gibbs notation.
 * What are the essential boundary conditions that can be applied to the plate?

Problem 2: Transient Heat Conduction
Consider the model problem shown in Figure 2.

A 4340 steel casing is used to contain a solid high energy (HE) material. The HE material reacts on the inner (star-shaped) surface and generates gases.Since the gases are confined, the temperature of the gases increases according to the ideal gas law ($$PV = nRT$$).

We will try to simulate the transient heat transfer and the associated thermal stresses in the HE material and the steel casing.

Since the actual problem is quite complex and involves chemical reactions, equations of state, and fluid-structure interaction, we will simplify the problem.

The assumptions we will make are:


 * 1) The burn-surface does not change position.
 * 2) The temperature is the same at all points on the burn surface.
 * 3) The temperature at the burn surface increases exponentially with time according to the relation

T = T_0 + A~\exp{(B~t)} $$
 * where $$T_0$$ is the initial temperature (300 K), $$T$$ is the temperature at time $$t$$, $$A = 2$$ K, and $$B = 3$$ K/s.

The inner radius of the steel casing is 5 cm. The outer radius of the steel casing is 6 cm.

The inner and outer radii of the star shaped region are 4 cm and 2 cm, respectively. You can make your own assumptions about the shape of the curves in the star shaped region.

Use standard values for the density, Young's modulus, Poisson's ratio, coefficient of thermal expansion, and thermal conductivity for the steel. You can use Matweb to find the properties you need.

Assume that the high energy material has the following properties:
 * Density = 1200 kg/m$$^3$$~.
 * Young's modulus = 1 GPa~.
 * Poisson's ratio = 0.49~.
 * Coefficient of thermal expansion = 1.2$$\times$$10$$^{-4}$$ /K~.
 * Thermal Conductivity = 5 W/m-K (this is 10 times higher than theactual value).
 * Specific Heat Capacity = 1000 J/kg-K.


 * 1) Assume that the initial temperature of the steel and the HE material is 300 K.Also assume that there is no heat flux on the external surface of the steel casing. Plot the temperature contours and stresses in the two materials at $$t = 1$$ s and $$t = 2$$ s.
 * 2) Assume that the initial temperature of the steel and the HE material is 300 K.Also assume that the external surface of the steel casing is maintained at 300 K. Plot the temperature contours and stresses in the two materials at $$t = 1$$ s and $$t = 2$$ s.

You should look at examples given in VM 33 and VM 174 in the ANSYS Verification manual for ideas on how to solve this problem.

Problem 3: Explicit Finite Elements
Figure 3(a) shows a beam that is built-in at one and loaded at the other. Figure 3(b) shows a plate that is built-in at one end and loaded at the other. The load that is applied is transient and has the form shown in Figure 3(c).

Assume that the beam and the plate are made of 4340 steel. The beam has dimensions 12 in. $$\times$$ 0.25 in. $$\times$$ 0.25 in. The plate has dimensions 12 in. $$\times$$ 0.25 in. $$\times$$ 6 in.

Simulate the problem for 10 secs. using LS-DYNA or some other explicit FE code. Use beam elements for the beam and shell elements for the plate.Use an element size of 2 in. for the beam and 2 in. $$\times$$ 2 in. for the shell.


 * 1) Plot the displacement of the end of the beam as a function of time.
 * 2) Plot the axial stress at two points on the top of the beam as a function of time.
 * 3) Plot the displacement of the end of the plate as a function of time.
 * 4) Plot the axial stress at two points on the top of the plate as a function of time.

After you run the LS-DYNA, the program will create a binary file 'd3plot' and ASCII file 'nodout'. To view the result:


 * 1) Run LS-Prepost and open the binary file 'd3plot'.
 * 2) To view the deflection at free end: click ASCII button.
 * 3) In the window below it, select 'nodout', then click LOAD.
 * 4) Available nodes will be shown in the next window below.
 * 5) Select the node.  For shell element, those are 2 nodes at the free end.
 * 6) Select 'y-displacement' for beam, and 'z-displacement' for shell, then click PLOT.

You should also try the animation tool to get a feel for how the beam and shell deform over time.