Thermoelastic Stress

Problem description

We solve a thermo-mechanical homogeneous model in cube domain. The dimensions of this cube model are 1 m in all directions. The boundary conditions and temperature loadings, as well as the material can refer Chapter 14 in Kolditz et al. for detailed problem description:

O. Kolditz et al. (2012): OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical (THM/C) processes in porous media. Environmental Earth Sciences, No. 2, vol. 67, p. 589-599, DOI:10.1007/s12665-012-1546-x

Results and evaluation

Result showing temperature and stresses development with time in the centre node of the model:

The analytical solution of stresses after heating is:

$$\begin{equation} \sigma_{xx} = \sigma_{yy} = \sigma_{zz} = - \frac{\alpha \Delta T E}{1 - 2 \nu} = - 3.260869\, \textrm{MPa} \end{equation}$$

The relative error between the numerical simulation and the analytical solution is 9.2⋅10-13.

References

O. Kolditz et al. (2012): OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical (THM/C) processes in porous media. Environmental Earth Sciences, No. 2, vol. 67, p. 589-599, DOI:10.1007/s12665-012-1546-x

This article was written by Xing-Yuan Miao. If you are missing something or you find an error please let us know.
Generated with Hugo 0.150.1 in CI job 743172 | Last revision: February 5, 2026
Commit: [py,cmake] Allow Python 3.14. 4d949dd99c  | Edit this page on