Part 1

The displacement field that satisfies the Kirchhoff-Love hypothesis is given by

${\displaystyle {\begin{aligned}u_{1}(x_{1},x_{2},x_{3})&=u_{1}^{0}(x_{1},x_{2})-x_{3}~u_{3,1}^{0}\\u_{2}(x_{1},x_{2},x_{3})&=u_{2}^{0}(x_{1},x_{2})-x_{3}~u_{3,2}^{0}\\u_{3}(x_{1},x_{2},x_{3})&=u_{3}^{0}(x_{1},x_{2})\end{aligned}}}$

where ${\displaystyle u_{i}^{0}}$ are the displacements of the midsurface.

Write the expanded form of the above equations after making the following substitutions.

${\displaystyle {\text{(2)}}\qquad {\begin{aligned}&x:=x_{1}~;~~y:=x_{2}~;~~z:=x_{3}\\&u:=u_{1}~;~~v:=u_{2}~;~~w:=u_{3}\\&u_{0}:=u_{1}^{0}~;~~v_{0}:=u_{2}^{0}~;~~w_{0}:=u_{3}^{0}~.\end{aligned}}}$

After making the substitutions we get

${\displaystyle {\begin{aligned}u(x,y,z)&=u_{0}(x,y)-z~w_{0,1}\\v(x,y,z)&=v_{0}(x,y)-z~w_{0,2}\\w(x,y,z)&=w_{0}(x,y)\end{aligned}}}$

The expanded form is

${\displaystyle {\begin{aligned}u(x,y,z)&=u_{0}(x,y)-z~{\frac {\partial w_{0}}{\partial x}}\\v(x,y,z)&=v_{0}(x,y)-z~{\frac {\partial w_{0}}{\partial y}}\\w(x,y,z)&=w_{0}(x,y)\end{aligned}}}$

Part 2

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

${\displaystyle {\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

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}\left(u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0}\right)-x_{3}~u_{3,\alpha \beta }^{0}\\\varepsilon _{\alpha 3}&=0\\\varepsilon _{33}&=0\end{aligned}}}$

where ${\displaystyle \alpha =1,2}$ and ${\displaystyle \beta =1,2}$.

The strain-displacement relations in index notation are

${\displaystyle \varepsilon _{ij}={\frac {1}{2}}\left[u_{i,j}+u_{j,i}\right]}$

where ${\displaystyle i,j=1,2,3}$. If the indices ${\displaystyle \alpha }$ and ${\displaystyle \beta }$ run from 1 to 2, we can write the above strain-displacement relations as

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}\left[u_{\alpha ,\beta }+u_{\beta ,\alpha }\right]\\\varepsilon _{\alpha 3}=\varepsilon _{3\alpha }&={\frac {1}{2}}\left[u_{\alpha ,3}+u_{3,\alpha }\right]\\\varepsilon _{33}&={\frac {1}{2}}\left[u_{3,3}+u_{3,3}\right]=u_{3,3}\end{aligned}}}$

The displacements are given by

${\displaystyle {\begin{aligned}u_{1}&=u_{1}^{0}-x_{3}~u_{3,1}^{0}\\u_{2}&=u_{2}^{0}-x_{3}~u_{3,2}^{0}\\u_{3}&=u_{3}^{0}\end{aligned}}}$

If the index ${\displaystyle \alpha }$ runs from 1 to 2, the displacements can be written as

${\displaystyle {\begin{aligned}u_{\alpha }&=u_{\alpha }^{0}-x_{3}~u_{3,\alpha }^{0}\\u_{3}&=u_{3}^{0}\end{aligned}}}$

The partial derivatives of the displacements are (since ${\displaystyle u_{1}^{0},u_{2}^{0},u_{3}^{0}}$ are functions only of ${\displaystyle x_{1}}$ and ${\displaystyle x_{2}}$)

${\displaystyle {\begin{aligned}u_{\alpha ,\beta }&=u_{\alpha ,\beta }^{0}-x_{3}~u_{3,\alpha \beta }^{0}\\u_{\alpha ,3}&=u_{\alpha ,3}^{0}-x_{3}~u_{3,\alpha 3}^{0}=0\\u_{3,\alpha }&=u_{3,\alpha }^{0}=0\\u_{3,3}&=u_{3,3}^{0}=0\end{aligned}}}$

Therefore,

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}\left[u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0}\right]-x_{3}~u_{3,\alpha \beta }^{0}\\\varepsilon _{\alpha 3}=\varepsilon _{3\alpha }&=0\\\varepsilon _{33}&=0\end{aligned}}}$

Part 3

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

${\displaystyle \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

${\displaystyle {\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

${\displaystyle {\begin{aligned}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_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0}\right)\\\varepsilon _{\alpha \beta }^{1}&=-u_{3,\alpha \beta }^{0}~.\end{aligned}}}$

Write down the expanded form of equation (3).

The strains are given by

${\displaystyle \varepsilon _{\alpha \beta }={\frac {1}{2}}\left[u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0}\right]-x_{3}~u_{3,\alpha \beta }^{0}=\varepsilon _{\alpha \beta }^{0}+x_{3}~\varepsilon _{\alpha \beta }^{1}}$

Therefore, the variations of the strains are

${\displaystyle \delta \varepsilon _{\alpha \beta }=\delta \varepsilon _{\alpha \beta }^{0}+x_{3}~\delta \varepsilon _{\alpha \beta }^{1}}$

The internal virtual work is

${\displaystyle {\begin{aligned}\delta W_{\text{int}}&=\int _{\Omega _{0}}\left[\int _{-h}^{h}\sigma _{\alpha \beta }~\left(\delta \varepsilon _{\alpha \beta }^{0}+x_{3}~\delta \varepsilon _{\alpha \beta }^{1}\right)~dx_{3}\right]d\Omega _{0}\\&=\int _{\Omega _{0}}\left[\delta \varepsilon _{\alpha \beta }^{0}\int _{-h}^{h}\sigma _{\alpha \beta }~dx_{3}+\delta \varepsilon _{\alpha \beta }^{1}\int _{-h}^{h}x_{3}~\sigma _{\alpha \beta }~dx_{3}\right]~d\Omega _{0}\\&=\int _{\Omega _{0}}\left[\delta \varepsilon _{\alpha \beta }^{0}~N_{\alpha \beta }+\delta \varepsilon _{\alpha \beta }^{1}~M_{\alpha \beta }\right]~d\Omega _{0}\end{aligned}}}$

Therefore,

${\displaystyle \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}}$

The expanded form of the above equation is

${\displaystyle {\begin{aligned}\delta W_{\text{int}}=\int _{\Omega _{0}}(&N_{11}\delta \varepsilon _{11}^{0}+N_{12}\delta \varepsilon _{12}^{0}+N_{21}\delta \varepsilon _{21}^{0}+N_{22}\delta \varepsilon _{22}^{0}+\\&M_{11}\delta \varepsilon _{11}^{1}+M_{12}\delta \varepsilon _{12}^{1}+M_{21}\delta \varepsilon _{21}^{1}+M_{22}\delta \varepsilon _{22}^{1})~d\Omega _{0}\end{aligned}}}$

or,

${\displaystyle {\begin{aligned}\delta W_{\text{int}}=\int _{\Omega _{0}}(&N_{11}\delta \varepsilon _{11}^{0}+2N_{12}\delta \varepsilon _{12}^{0}+N_{22}\delta \varepsilon _{22}^{0}+\\&M_{11}\delta \varepsilon _{11}^{1}+2M_{12}\delta \varepsilon _{12}^{1}+M_{22}\delta \varepsilon _{22}^{1})~d\Omega _{0}\end{aligned}}}$

Part 4

The equilibrium equations for the plate can be written as

${\displaystyle {\begin{aligned}{\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{aligned}}}$

where

${\displaystyle {\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.

The equilibrium equations are (written using indices ${\displaystyle 1,2,3}$):

${\displaystyle {\begin{aligned}&N_{11,1}+N_{12,2}=0\\&N_{21,1}+N_{22,2}=0\\&M_{11,11}+2M_{12,12}+M_{22,22}+\left(N_{11}~u_{3,1}^{0}+N_{12}~u_{3,2}^{0}\right)_{,1}+\left(N_{21}~u_{3,1}^{0}+N_{22}~u_{3,2}^{0}\right)_{,2}+q=0\end{aligned}}}$

Using the summation convention, we get

${\displaystyle {\begin{aligned}&N_{1\beta ,\beta }=0\\&N_{2\beta ,\beta }=0\\&M_{\alpha \beta ,\alpha \beta }+\left(N_{1\beta }~u_{3,\beta }^{0}\right)_{,1}+\left(N_{2\beta }~u_{3,\beta }^{0}\right)_{,2}+q=0\end{aligned}}}$

The first two equations can then be written as

${\displaystyle N_{\alpha \beta ,\beta }=0}$

and the third equation can be written as

${\displaystyle M_{\alpha \beta ,\alpha \beta }+\left(N_{\alpha \beta }~u_{3,\beta }^{0}\right)_{,\alpha }+q=0}$

Therefore the equilibrium equations in index notation are

${\displaystyle {\begin{aligned}&N_{\alpha \beta ,\beta }=0\\&M_{\alpha \beta ,\alpha \beta }+\left(N_{\alpha \beta }~u_{3,\beta }^{0}\right)_{,\alpha }+q=0\end{aligned}}}$

Part 5

Show that in index-free Gibbs notation the equilibrium equations for the plate can be written as

${\displaystyle {\begin{aligned}{\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_{3}^{0})}}+q&=0{\text{(5)}}\qquad \end{aligned}}}$

where

${\displaystyle {\begin{aligned}{\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{aligned}}}$

and ${\displaystyle \alpha =1,2}$, ${\displaystyle \beta =1,2}$, and ${\displaystyle \gamma =1,2}$.

Equation (4) in index notation is

${\displaystyle {\boldsymbol {\nabla }}\bullet {\boldsymbol {N}}=N_{\alpha \beta ,\beta }~\mathbf {e} _{\alpha }=\mathbf {0} \qquad \implies N_{\alpha \beta ,\beta }=0}$

The above equations have been shown to be equivalent to the first two equilibrium equations in the previous part.

Let us now look at equation (5). Let

${\displaystyle \mathbf {p} :={\boldsymbol {\nabla }}\bullet {\boldsymbol {M}}\qquad \implies p_{\alpha }\mathbf {e} _{\alpha }=M_{\alpha \beta ,\beta }\mathbf {e} _{\alpha }}$

Then

${\displaystyle {\boldsymbol {\nabla }}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})={\boldsymbol {\nabla }}\bullet \mathbf {p} =p_{\alpha ,\alpha }=M_{\alpha \beta ,\alpha \beta }}$

For the next term, let

${\displaystyle \mathbf {w} :={\boldsymbol {\nabla }}u_{3}^{0}\qquad \implies w_{\alpha }\mathbf {e} _{\alpha }=u_{3,\alpha }^{0}\mathbf {e} _{\alpha }\implies w_{\alpha }=u_{3,\alpha }^{0}}$

Then,

${\displaystyle \mathbf {q} :={\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0}={\boldsymbol {N}}\bullet \mathbf {w} \qquad \implies q_{\alpha }~\mathbf {e} _{\alpha }=N_{\alpha \beta }~w_{\beta }~\mathbf {e} _{\alpha }=N_{\alpha \beta }~u_{3,\beta }^{0}~\mathbf {e} _{\alpha }}$

Hence

${\displaystyle {\boldsymbol {\nabla }}\bullet ({\boldsymbol {N}}\bullet {\boldsymbol {\nabla u_{3}^{0})}}={\boldsymbol {\nabla }}\bullet \mathbf {q} =q_{\alpha ,\alpha }=(N_{\alpha \beta }~u_{3,\beta }^{0})_{,\alpha }}$

Therefore, equation (5) can be written as

${\displaystyle M_{\alpha \beta ,\alpha \beta }+(N_{\alpha \beta }~u_{3,\beta }^{0})_{,\alpha }+q=0}$

The above equation has been shown to be equivalent to the third equilibrium equations in the previous part.

Part 6

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

For equation 4, we multiply the equation by a vector weighting function ${\displaystyle \mathbf {w} _{1}}$ and integrate over the midsurface ${\displaystyle \Omega _{0}}$. Then,

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}\mathbf {w} _{1}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {N}})~d\Omega _{0}&=0\end{aligned}}}$

Using the identity

${\displaystyle {\boldsymbol {\nabla }}\bullet ({\boldsymbol {S}}^{T}\bullet \mathbf {v} )=\mathbf {v} \bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {S}})+{\boldsymbol {S}}:({\boldsymbol {\nabla }}\mathbf {v} )}$

we get

${\displaystyle \mathbf {w} _{1}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {N}})={\boldsymbol {\nabla }}\bullet ({\boldsymbol {N}}^{T}\bullet \mathbf {w} _{1})-{\boldsymbol {N}}:({\boldsymbol {\nabla }}\mathbf {w} _{1})}$

From the symmetry to the stress tensor ${\displaystyle {\boldsymbol {N}}={\boldsymbol {N}}^{T}}$. Hence,

${\displaystyle \mathbf {w} _{1}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {N}})={\boldsymbol {\nabla }}\bullet ({\boldsymbol {N}}\bullet \mathbf {w} _{1})-{\boldsymbol {N}}:({\boldsymbol {\nabla }}\mathbf {w} _{1})}$

Therefore

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}\mathbf {w} _{1}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {N}})~d\Omega _{0}&=\int _{\Omega _{0}}{\boldsymbol {\nabla }}\bullet ({\boldsymbol {N}}\bullet \mathbf {w} _{1})~d\Omega _{0}-\int _{\Omega _{0}}{\boldsymbol {N}}:({\boldsymbol {\nabla }}\mathbf {w} _{1})~d\Omega _{0}\\&=\int _{\Gamma _{0}}\mathbf {n} \bullet ({\boldsymbol {N}}\bullet \mathbf {w} _{1})~d\Gamma _{0}-\int _{\Omega _{0}}{\boldsymbol {N}}:({\boldsymbol {\nabla }}\mathbf {w} _{1})~d\Omega _{0}\end{aligned}}}$

where ${\displaystyle \mathbf {n} }$ is the normal to the boundary ${\displaystyle \Gamma _{0}}$ of the midsurface.

Hence the symmetric weak form of the first equation is

${\displaystyle \int _{\Omega _{0}}{\boldsymbol {N}}:{\boldsymbol {\nabla }}\mathbf {w} _{1}~d\Omega _{0}=\int _{\Gamma _{0}}\mathbf {n} \bullet ({\boldsymbol {N}}\bullet \mathbf {w} _{1})~d\Gamma _{0}}$

For equation 5, we multiply the equation by a scalar weighting function ${\displaystyle w_{2}}$ and integrate over the midsurface ${\displaystyle \Omega _{0}}$. Then,

${\displaystyle \int _{\Omega _{0}}w_{2}\left[{\boldsymbol {\nabla }}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})+{\boldsymbol {\nabla }}\bullet ({\boldsymbol {N}}\bullet {\boldsymbol {\nabla u_{3}^{0})}}+q\right]~d\Omega _{0}=0}$

Using the identity,

${\displaystyle {\boldsymbol {\nabla }}\bullet (\varphi \mathbf {v} )=\varphi ({\boldsymbol {\nabla }}\bullet \mathbf {v} )+\mathbf {v} \bullet ({\boldsymbol {\nabla }}\varphi )}$

we get

${\displaystyle w_{2}[{\boldsymbol {\nabla }}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})]={\boldsymbol {\nabla }}\bullet [w_{2}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})]-({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})\bullet ({\boldsymbol {\nabla }}w_{2})}$

Therefore,

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}w_{2}[{\boldsymbol {\nabla }}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})]~d\Omega _{0}&=\int _{\Omega _{0}}{\boldsymbol {\nabla }}\bullet [w_{2}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})]~d\Omega _{0}-\int _{\Omega _{0}}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}\\&=\int _{\Gamma _{0}}\mathbf {n} \bullet [w_{2}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})]~d\Gamma _{0}-\int _{\Omega _{0}}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}\\&=\int _{\Gamma _{0}}w_{2}~[\mathbf {n} \bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})]~d\Gamma _{0}-\int _{\Omega _{0}}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}\end{aligned}}}$

In order to have natural BCs directly in terms of moments, we have to perform a similar set of operations the last term above. Thus,

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}&=\int _{\Omega _{0}}{\boldsymbol {\nabla }}\bullet [{\boldsymbol {M}}\bullet ({\boldsymbol {\nabla w_{2})]}}~d\Omega _{0}-\int _{\Omega _{0}}{\boldsymbol {M}}:[{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}{w_{2})}]~d\Omega _{0}\\&=\int _{\Gamma _{0}}\mathbf {n} \bullet [{\boldsymbol {M}}\bullet {({\boldsymbol {\nabla }}w_{2})}]~d\Gamma _{0}-\int _{\Omega _{0}}{\boldsymbol {M}}:[{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}w_{2})]~d\Omega _{0}\end{aligned}}}$

Hence,

${\displaystyle \int _{\Omega _{0}}w_{2}~{\boldsymbol {\nabla }}\bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})~d\Omega _{0}=\int _{\Gamma _{0}}\left\{w_{2}~\mathbf {n} \bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})-\mathbf {n} \bullet {\boldsymbol {M}}\bullet ({\boldsymbol {\nabla }}w_{2})\right\}d\Gamma _{0}+\int _{\Omega _{0}}{\boldsymbol {M}}:\left({\boldsymbol {\nabla }}{\boldsymbol {\nabla }}{w_{2}}\right)~d\Omega _{0}}$

Similarly,

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}w_{2}\left[{\boldsymbol {\nabla }}\bullet ({\boldsymbol {N}}\bullet {\boldsymbol {\nabla u_{3}^{0})}}\right]~d\Omega _{0}&=\int _{\Omega _{0}}{\boldsymbol {\nabla }}\bullet [w_{2}({\boldsymbol {N}}\bullet {\boldsymbol {\nabla u_{3}^{0})]}}~d\Omega _{0}-\int _{\Omega _{0}}({\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}\\&=\int _{\Gamma _{0}}\mathbf {n} \bullet [w_{2}(\mathbf {d} ot{{\boldsymbol {N}}{{\boldsymbol {\nabla }}u_{3}^{0}})]}~d\Gamma _{0}-\int _{\Omega _{0}}({\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}\\&=\int _{\Gamma _{0}}w_{2}~[\mathbf {n} \bullet ({\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0})]~d\Gamma _{0}-\int _{\Omega _{0}}({\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0})\bullet ({\boldsymbol {\nabla }}w_{2})~d\Omega _{0}\end{aligned}}}$

Therefore, the symmetric weak form of the second equation is

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}&\left[({\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0})\bullet {\boldsymbol {\nabla }}w_{2}-{\boldsymbol {M}}:\left({\boldsymbol {\nabla }}{\boldsymbol {\nabla }}w_{2}\right)\right]~d\Omega _{0}=\\&\int _{\Omega _{0}}w_{2}~q~d\Omega _{0}+\int _{\Gamma _{0}}\left[w_{2}\left\{\mathbf {n} \bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}})+\mathbf {n} \bullet {\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0}\right\}-\mathbf {n} \bullet ({\boldsymbol {M}}\bullet {\boldsymbol {\nabla }}w_{2})\right]~d\Gamma _{0}\end{aligned}}}$

Using the distributive law for vectors

${\displaystyle (\mathbf {a} +\mathbf {b} )\bullet \mathbf {c} =\mathbf {a} \bullet \mathbf {c} +\mathbf {b} \bullet \mathbf {c} }$

we can also write the symmetric weak form for the second equation as:

${\displaystyle {\begin{aligned}\int _{\Omega _{0}}&\left[({\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0})\bullet {{\boldsymbol {\nabla }}w_{2}}-{\boldsymbol {M}}:\left({\boldsymbol {\nabla }}{\boldsymbol {\nabla }}w_{2}\right)\right]~d\Omega _{0}=\\&\int _{\Omega _{0}}w_{2}~q~d\Omega _{0}+\int _{\Gamma _{0}}\mathbf {n} \bullet \left({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}}+{\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0}\right)~w_{2}~d\Gamma _{0}-\int _{\Gamma _{0}}\mathbf {n} \bullet ({\boldsymbol {M}}\bullet {\boldsymbol {\nabla }}w_{2})~d\Gamma _{0}\end{aligned}}}$

Part 7

Write down the natural boundary conditions in Gibbs notation.

From the expressions for the weak form, the weighted natural boundary conditions are

${\displaystyle {\begin{aligned}&\int _{\Gamma _{0}}\mathbf {n} \bullet (\mathbf {d} ot{{\boldsymbol {N}}{\mathbf {w} _{1}})}~d\Gamma _{0}\\&\int _{\Gamma _{0}}\mathbf {n} \bullet \left({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}}+{\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0}\right)~w_{2}~d\Gamma _{0}\\&-\int _{\Gamma _{0}}\mathbf {n} \bullet (\mathbf {d} ot{{\boldsymbol {M}}{{\boldsymbol {\nabla }}w_{2}})}~d\Gamma _{0}\end{aligned}}}$

Therefore, the natural boundary conditions are:

${\displaystyle {\begin{aligned}{\text{Boundary Membrane Forces}}~&&\mathbf {n} \bullet {\boldsymbol {N}}\\{\text{Boundary Transverse Forces}}~&&\mathbf {n} \bullet \left({\boldsymbol {\nabla }}\bullet {\boldsymbol {M}}+{\boldsymbol {N}}\bullet {\boldsymbol {\nabla }}u_{3}^{0}\right)\\{\text{Boundary Moments}}~&&\mathbf {n} \bullet {\boldsymbol {M}}\end{aligned}}}$

Part 8

What are the essential boundary conditions that can be applied to the plate?

The essential boundary conditions that can be applied are

${\displaystyle {\begin{aligned}{\text{Boundary Membrane Displacements}}~&&u_{1}^{0},u_{2}^{0}\\{\text{Boundary Transverse Displacements}}~&&u_{3}^{0}\\{\text{Boundary Out-of-plane Rotations}}~&&{\boldsymbol {\nabla }}u_{3}^{0}\end{aligned}}}$

Part 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 ${\displaystyle t=1}$ s and ${\displaystyle t=2}$ s.

The temperature contours for the two materials are shown in Figure 7. The low thermal conductivity of the HE material leads to a slow rate of diffusion of heat into the material. Therefore, the steel casing sees no change in temperature. However, the casing does deform due to the deformation of the HE material.

Figure 7(a). Temperature contours for no flux BC at t = 1 sec. Units are degrees K.

Figure 7(b). Temperature contours for no flux BC at t = 2 sec. Units are degrees K.

The stress contours for the two materials are shown in Figures 8, 9, 10, and 11. The principal stress ${\displaystyle \sigma _{1}}$ is tensile in the steel casing and reaches a values of around 200 MPa (close to its yield point). In the HE material, ${\displaystyle \sigma _{1}}$ is compressive in some regions and tensile in others and reaches values well beyond the ulimate tensile strength of the material inside some elements.

Figure 8(a). First Principal Stress contours for no flux BC at t = 1 sec. Units are Pa.

Figure 8(b). First Principal Stress contours for no flux BC at t = 2 sec. Units are Pa.

Figure 9(a). Second Principal Stress contours for no flux BC at t = 1 sec. Units are Pa.

Figure 9(b). Second Principal Stress contours for no flux BC at t = 2 sec. Units are Pa.

Figure 10(a). Third Principal Stress contours for no flux BC at t = 1 sec. Units are Pa.

Figure 10(b). Third Principal Stress contours for no flux BC at t = 2 sec. Units are Pa.

Figure 11(a). Shear Stress contours for no flux BC at t = 1 sec. Units are Pa.

Figure 11(b). Shear Stress contours for no flux BC at t = 2 sec. Units are Pa.

The ANSYS input file used for this problem is given below. The preprocessing stage is shown in the PDF file NFE_HW8Prob2Ans.pdf.

Figure 12. Preprocessing ANSYS input.

The solution stage is shown below:

!
! Solve
!
/solu
antype, trans
timint, off, struc
cnvtol, heat
cnvtol, f
autots, on
solcontrol, on
!
! Initial condition
!
ic, all, temp, 300
!
! Time stepping
!
time, 2
deltim, 0.1
outres, all, 10
solve
finish

Part 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 ${\displaystyle t=1}$ s and ${\displaystyle t=2}$ s.

The temperature distribution from the previous problem suggests that the heat flux at the external surface is zero if the outside is maintained at 300 K (because of the low thermal conductivity of the HE material.) We don't expect to see any significant difference in the temperature and stress distributions.

Part 1

Plot the displacement of the end of the beam as a function of time.

Figure 14 shows the displacement of the end of the beam as a function of time. Some numerical damping appears to occur if Hughes-Liu beam elements are used with 2${\displaystyle \times }$2 quadrature.

Figure 14. Displacement (in) of the end of the beam as a function of time (secs).

Part 2

Plot the axial stress at two points on the top of the beam as a function of time.

Figure 15 shows the axial stress at a Gauss point near the top of the beam for elements 2 and 5.

Figure 15 (a). Axial stress (psi) at element 2 at the top of the beam as a function of time (secs).

Figure 15 (b). Axial stress (psi) at element 5 at the top of the beam as a function of time (secs).

The LS-Dyna keyword file for the beam is shown below.

*keyword
*title
HW8Prob3_1.k - cantilever beam
*control_termination
10
*database_binary_d3plot
0.1
*database_history_node
1, 2, 3, 4, 5, 6, 7
*database_nodout
0.1
*database_history_beam
1, 2, 3, 4, 5, 6
*database_elout
0.1
*part
steel beam
1,1,1
*section_beam
1, 1, 0.833, 2.0, 0.0
0.25, 0.25, 0.25, 0.25, 0.0, 0.0
*mat_elastic
1, 0.284, 29.7e6, 0.29
*node
1, 0, 0, 0, 7, 7
2, 2.0, 0, 0, 0, 0
3, 4.0, 0, 0, 0, 0
4, 6.0, 0, 0, 0, 0
5, 8.0, 0, 0, 0, 0
6, 10.0, 0, 0, 0, 0
7, 12.0, 0, 0, 0, 0
8, 5.0, 1.0, 0, 0, 0
*element_beam
1, 1, 1, 2, 8, 0, 0, 0, 0, 2
2, 1, 2, 3, 8, 0, 0, 0, 0, 2
3, 1, 3, 4, 8, 0, 0, 0, 0, 2
4, 1, 4, 5, 8, 0, 0, 0, 0, 2
5, 1, 5, 6, 8, 0, 0, 0, 0, 2
6, 1, 6, 7, 8, 0, 0, 0, 0, 2
*load_node_point
7, 2, 1, 1.0
*define_curve
1
0.0, 0.0
0.1, 200.0
0.2, 200.0
0.25, 0.0
10.0, 0.0
*end

Part 3

Plot the displacement of the end of the plate as a function of time.

Figure 16 shows the displacement of the end of the plate as a function of time. Hughes-Liu shell elements with 3 integration points are used.

Figure 16. Displacement of the end of the plate as a function of time.

Part 4

Plot the axial stress at two points on the top of the plate as a function of time.

Figure 17 shows the axial stress at a Gauss point near the top of the plate for elements 8 and 11.

Figure 17(a). Axial stress at element 8 at the top of the plate as a function of time.

Figure 17(b). Axial stress at element 11 at the top of the plate as a function of time.

The LS-Dyna keyword file for the beam is shown below.

*keyword
*title
HW8Prob3_2 - bending of shell element
*control_termination
10
*database_binary_d3plot
0.1
*database_history_node_set
2
*database_nodout
0.1
*database_history_shell_set
1
*database_elout
0.1
*part
steel shell
1,1,1
*section_shell
1, 1, 0.833, 3.0, 0.0
0.25, 0.25, 0.25, 0.25, 0.0, 0.0
*mat_elastic
1, 0.284, 29.7e6, 0.29
*node
1,0,0,0
2,2,0,0
3,4,0,0
4,6,0,0
5,8,0,0
6,10,0,0
7,12,0,0
8,0,2,0
9,2,2,0
10,4,2,0
11,6,2,0
12,8,2,0
13,10,2,0
14,12,2,0
15,0,4,0
16,2,4,0
17,4,4,0
18,6,4,0
19,8,4,0
20,10,4,0
21,12,4,0
22,0,6,0
23,2,6,0
24,4,6,0
25,6,6,0
26,8,6,0
27,10,6,0
28,12,6,0
*set_node_list
1
1, 8, 15, 22
*set_node_list
2
7, 14, 21, 28
*element_shell
1,1,1,2,9,8
2,1,2,3,10,9
3,1,3,4,11,10
4,1,4,5,12,11
5,1,5,6,13,12
6,1,6,7,14,13
7,1,8,9,16,15
8,1,9,10,17,16
9,1,10,11,18,17
10,1,11,12,19,18
11,1,12,13,20,19
12,1,13,14,21,20
13,1,15,16,23,22
14,1,16,17,24,23
15,1,17,18,25,24
16,1,18,19,26,25
17,1,19,20,27,26
18,1,20,21,28,27
*set_shell_list
1
8, 11
*boundary_spc_set
1, 0, 1, 1, 1, 1, 1, 1
*load_node_set
2, 3, 1, 1.0
*define_curve
1
0.0, 0.0
0.1, 200.0
0.2, 200.0
0.25, 0.0
10.0, 0.0
*end