Elasticity/Principle of minimum potential energy

Strain Energy Density

The strain energy density is defined as

$U({\boldsymbol {\varepsilon }})=\int _{0}^{\varepsilon }{\boldsymbol {\sigma }}:d{\boldsymbol {\varepsilon }}$

If the strain energy density is path independent, then it acts as a potential for stress, i.e.,

$\sigma _{ij}={\frac {\partial U(\varepsilon )}{\partial \varepsilon _{ij}}}$

For adiabatic processes, $U$ is equal to the change in internal energy per unit volume.

For isothermal processes, $U$ is equal to the Helmholtz free energy per unit volume.

The natural state of a body is defined as the state in which the body is in stable thermal equilibrium with no external loads and zero stress and strain.

When we apply energy methods in linear elasticity, we implicitly assume that a body returns to its natural state after loads are removed. This implies that the Gibb's condition is satisfied :

$U({\boldsymbol {\varepsilon }})\geq 0~~{\text{with}}~~U({\boldsymbol {\varepsilon }})=0~~{\text{iff}}~~{\boldsymbol {\varepsilon }}=0$

Principle of Minimum Potential Energy

This principle states that

• If the prescribed traction and body force fields are independent of the deformation
• then the actual displacement field makes the potential energy functional an absolute minimum.

In other words, the principle of minimum potential energy states that the potential energy functional

$\Pi [\mathbf {u} ]={\frac {1}{2}}\int _{\mathcal {B}}{\boldsymbol {\sigma }}:{\boldsymbol {\varepsilon }}~dV-\int _{\mathcal {B}}\mathbf {f} \bullet \mathbf {u} ~dV-\int _{\partial {\mathcal {B}}^{t}}\mathbf {t} \bullet \mathbf {u} ~dA$

is minimized by the actual displacement field.

Proof

The first step in the proof is to show that the actual displacements make the function $\Pi \,$ stationary. The second step is to show that the stationary point is actually the minimum.

Proof of the Principle of Stationary Potential Energy

The first variation of the potential energy functional $\Pi \,$ is

$\delta \Pi [\mathbf {u} ,\delta \mathbf {u} ]=\int _{B}C_{ijkl}~u_{i,j}~\delta u_{k,l}~dV-\int _{B}{\tilde {f}}_{i}~\delta u_{i}~dV-\int _{\partial B^{t}}{\tilde {t}}_{i}~\delta u_{i}~dA$

or,

${\begin{aligned}\delta \Pi [\mathbf {u} ,\delta \mathbf {u} ]&=\int _{B}\sigma _{ji}~\delta u_{i,j}~dV-\int _{B}{\tilde {f}}_{i}~\delta u_{i}~dV-\int _{\partial B^{t}}{\tilde {t}}_{i}~\delta u_{i}~dA\\&=\int _{B}(\sigma _{ji}~\delta u_{i})_{,j}~dV-\int _{B}\sigma _{ji,j}~\delta u_{i}~dV-\int _{B}{\tilde {f}}_{i}~\delta u_{i}~dV-\int _{\partial B^{t}}{\tilde {t}}_{i}~\delta u_{i}~dA\\&=\int _{\partial B}\sigma _{ji}~\delta u_{i}~n_{j}~dA-\int _{B}\sigma _{ji,j}~\delta u_{i}~dV-\int _{B}{\tilde {f}}_{i}~\delta u_{i}~dV-\int _{\partial B^{t}}{\tilde {t}}_{i}~\delta u_{i}~dA\\&=\int _{\partial B^{t}}\sigma _{ji}~n_{j}~\delta u_{i}~dA-\int _{\partial B^{t}}{\tilde {t}}_{i}~\delta u_{i}~dA-\int _{B}\sigma _{ji,j}~\delta u_{i}~dV-\int _{B}{\tilde {f}}_{i}~\delta u_{i}~dV\\&=\int _{\partial B^{t}}\left[\sigma _{ji}~n_{j}-{\tilde {t}}_{i}\right]~\delta u_{i}~dA-\int _{B}\left[\sigma _{ji,j}+{\tilde {f}}_{i}\right]~\delta u_{i}~dV\end{aligned}}$

or,

$\delta \Pi [\mathbf {u} ,\delta \mathbf {u} ]=\int _{\partial B^{t}}\left[{\widehat {\mathbf {n} }}{}\bullet {\boldsymbol {\sigma }}-{\tilde {\mathbf {t} }}\right]\bullet \delta \mathbf {u} ~dA-\int _{B}\left[{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}+{\tilde {\mathbf {f} }}\right]\bullet \delta \mathbf {u} ~dV$

Therefore, $\delta \Pi [\mathbf {u} ,\delta \mathbf {u} ]=0$ (i.e. $\Pi$ is stationary) only if

${\begin{aligned}{\widehat {\mathbf {n} }}{}\bullet {\boldsymbol {\sigma }}&={\tilde {\mathbf {t} }}\\{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}+{\tilde {\mathbf {f} }}&=0\end{aligned}}$

which are the conditions that only the actual displacement field satisfies.

Proof of the Principle of Minimum Potential Energy

To prove that $\Pi [\mathbf {u} ,\delta \mathbf {u} ]$ is not only stationary, but also the global minimum all we now need to show is that

$\delta ^{2}\Pi [\mathbf {u} ,\delta \mathbf {u} ]>0$

Now,

$\delta ^{2}\Pi [\mathbf {u} ,\delta \mathbf {u} ]=\int _{B}C_{ijkl}~\delta u_{i,j}~\delta u_{k,l}~dV$

If the displacement field is a pure rigid body motion, then the strain energy density

$U={\frac {1}{2}}C_{ijkl}w_{ij}w_{kl}=0$

where $\mathbf {w}$ is the spin tensor given by

$\mathbf {w} ={\frac {1}{2}}\left({\boldsymbol {\nabla }}u-{\boldsymbol {\nabla }}u^{T}\right)~~;~~w_{ij}={\frac {1}{2}}\left(u_{i,j}-u_{j,i}\right)$

If the displacement field does not contain any rigid body motion, then the strain energy density is given by

$U={\frac {1}{2}}C_{ijkl}\varepsilon _{ij}\varepsilon _{kl}>0$

where ${\boldsymbol {\varepsilon }}$ is the strain tensor given by

${\boldsymbol {\varepsilon }}={\frac {1}{2}}\left({\boldsymbol {\nabla }}u+{\boldsymbol {\nabla }}u^{T}\right)~~;~~\varepsilon _{ij}={\frac {1}{2}}\left(u_{i,j}+u_{j,i}\right)$

Therefore, for a displacement field containing both spin and strain

$U={\frac {1}{2}}C_{ijkl}~\delta u_{i,j}~\delta u_{k,l}\geq 0$

or,

$\delta ^{2}\Pi [\mathbf {u} ,\delta \mathbf {u} ]=\int _{B}C_{ijkl}~\delta u_{i,j}~\delta u_{k,l}~dV\geq 0$

This means that $\delta ^{2}\Pi >0$ for all values of $\delta \mathbf {u}$ other than rigid body motion, in which case $\delta ^{2}\Pi =0$.

Hence, for mixed boundary value problems $\delta ^{2}\Pi >0$ for all $\delta \mathbf {u}$, as long as the displacement BCs prevent rigid body motions. Therefore, the potential energy functional is minimized by the actual displacement field.

Elasticity