The complementary strain energy density is given by
![{\displaystyle U^{c}({\boldsymbol {\sigma }})=\int _{0}^{\boldsymbol {\sigma }}\varepsilon _{ij}~d\sigma _{ij}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1abcbe90973b6588356e4b32c20bea690a0e02b6)
For linear elastic materials
![{\displaystyle U^{c}={\frac {1}{2}}{\boldsymbol {\sigma }}:{\boldsymbol {\varepsilon }}=U}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8298a4e843261747d111cd2441931c69369bc0bb)
Let
be the set of all statically admissible states of stress.
Let the complementary energy functional be
![{\displaystyle \Pi ^{c}[{\boldsymbol {\sigma }}]=\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }})~dV-\int _{\partial {\mathcal {B}}^{u}}\mathbf {t} \bullet {\widehat {\mathbf {u} }}~dA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/54a21bbc8fd996c8752aeffffd08e456fefadd33)
Then the Principle of Stationary Complementary Energy states that:
The Principle of Minimum Complementary Energy states that:
For linear elastic materials, the complementary energy functional is rendered an absolute minimum by the actual stress field.
Note that the complementary energy corresponding to the actual stress field is the negative of the potential energy corresponding to the actual displacement field.
Let
be a solution of a the mixed boundary value problem of linear elasticity.
Let
.
Define
![{\displaystyle {\boldsymbol {\sigma }}^{'}={\tilde {\boldsymbol {\sigma }}}-{\boldsymbol {\sigma }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/693dce697665cee82fe4125318916d79d8cee658)
Then
satisfies equilibrium and the traction BCs, i.e.,
![{\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\bullet {{\boldsymbol {\sigma }}^{'}}=0&~~\forall \mathbf {x} \in {\mathcal {B}}\\\mathbf {t} ^{'}={\widehat {\mathbf {n} }}{}\bullet {\boldsymbol {\sigma }}^{'}=0&~~\forall \mathbf {x} \in \partial {\mathcal {B}}^{t}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/236dc66adc60bb3aff9ca7baba30f440cfd0b294)
Since
![{\displaystyle {\boldsymbol {\varepsilon }}={\text{S}}:{\boldsymbol {\sigma }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/24fa184d9c720f2e55716b47e8aafe6e29c8d207)
and
![{\displaystyle \int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }}^{'}+{\boldsymbol {\sigma }})~dV=\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }}^{'})~dV+\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }})~dV+\int _{\mathcal {B}}{\boldsymbol {\sigma }}^{'}:({\text{S}}:{\boldsymbol {\sigma }})~dV}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cdcda43a7e4c4134775be3a37a5de4275e8db58a)
we have,
![{\displaystyle {\text{(1)}}\qquad \int _{\mathcal {B}}\left[U^{c}({\tilde {\boldsymbol {\sigma }}})-U^{c}({\boldsymbol {\sigma }})\right]~dV=\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }}^{'})~dV+\int _{\mathcal {B}}{\boldsymbol {\sigma }}^{'}:{\boldsymbol {\varepsilon }}~dV}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bed8f1be60bf2bdf14f1158983177e1f256b48e7)
We can also show that
![{\displaystyle \int _{\partial {\mathcal {B}}}\mathbf {t} \bullet \mathbf {u} ~dA=\int _{\mathcal {B}}\mathbf {u} \bullet ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})~dV+{\frac {1}{2}}\int _{\mathcal {B}}{\boldsymbol {\sigma }}:({\boldsymbol {\nabla }}\mathbf {u} +{\boldsymbol {\nabla }}\mathbf {u} ^{T})~dV}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4066432ec8a3a46b3de4a8db3a9c365c9cca7856)
Therefore,
![{\displaystyle {\text{(2)}}\qquad \int _{\partial {\mathcal {B}}^{u}}\mathbf {t} ^{'}\bullet \mathbf {u} ~dA=\int _{\mathcal {B}}{\boldsymbol {\sigma }}^{'}:{\boldsymbol {\varepsilon }}~dV}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ce1767f827a538c72201ad601d854791248d905)
Now,
![{\displaystyle \Pi ^{c}[{\boldsymbol {\sigma }}]=\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }})~dV-\int _{\partial {\mathcal {B}}^{u}}\mathbf {t} \bullet {\widehat {\mathbf {u} }}~dA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/54a21bbc8fd996c8752aeffffd08e456fefadd33)
Hence,
![{\displaystyle \Pi ^{c}[{\tilde {\boldsymbol {\sigma }}}]=\int _{\mathcal {B}}U^{c}({\tilde {\boldsymbol {\sigma }}})~dV-\int _{\partial {\mathcal {B}}^{u}}(\mathbf {t} +\mathbf {t'} )\bullet {\widehat {\mathbf {u} }}~dA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd7840300b2e6feebaa7ea7b285d0a1492c28df7)
Therefore,
![{\displaystyle {\text{(3)}}\qquad \Pi ^{c}[{\tilde {\boldsymbol {\sigma }}}]-\Pi ^{c}[{\boldsymbol {\sigma }}]=\int _{\mathcal {B}}\left[U^{c}({\tilde {\boldsymbol {\sigma }}})-U^{c}({\boldsymbol {\sigma }})\right]~dV-\int _{\partial {\mathcal {B}}^{u}}\mathbf {t} ^{'}\bullet {\widehat {\mathbf {u} }}~dA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dddfee691f92e378f2b32f2ad2ed2d243fbf87f8)
From equations (1), (2), and (3), we have,
![{\displaystyle \Pi ^{c}[{\tilde {\boldsymbol {\sigma }}}]-\Pi ^{c}[{\boldsymbol {\sigma }}]=\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }}^{'})~dV+\int _{\partial {\mathcal {B}}^{u}}\mathbf {t} ^{'}\bullet \mathbf {u} ~dA-\int _{\partial {\mathcal {B}}^{u}}\mathbf {t} ^{'}\bullet {\widehat {\mathbf {u} }}~dA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f80e1ef5f737c0f817d98587d2c6766b35621304)
Since
on
, we have,
![{\displaystyle \Pi ^{c}[{\tilde {\boldsymbol {\sigma }}}]-\Pi ^{c}[{\boldsymbol {\sigma }}]=\int _{\mathcal {B}}U^{c}({\boldsymbol {\sigma }}^{'})~dV>0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b149668b728a4082c07eb5640604cc525caaac29)
Hence, proved.