Methodology
============

It is a double simulation, the first in thermo-mechanics, the second in pure mechanics. The first will be validated in comparison with the second, assuming of course that the tested behavior provides a correct solution in pure mechanics.


Simulation 1
------------

The first simulation (thermo-mechanical solution that we are trying to validate) consists in applying a temperature variation to a material point, by blocking the deformations in the direction :math:`x` :math:`{\varepsilon }_{\mathrm{xx}}=0`. The imposed temperature increases linearly as a function of time. The temperature varies from :math:`{T}_{0}\mathrm{=}20°C` to :math:`{T}_{\mathit{max}}\mathrm{=}500°C`. The transient consists of NCAL steps. The reference temperature is :math:`{T}_{\mathit{réf}}\mathrm{=}20°C`.


Simulation 2
------------

The second simulation (which must be equivalent to the first) consists in applying to the same material point an imposed deformation following :math:`x`: :math:`{\varepsilon }_{\mathit{xx}}\mathrm{=}\mathrm{-}{\varepsilon }^{\mathit{th}}\mathrm{=}\mathrm{-}\alpha (T)(T\mathrm{-}{T}_{\mathit{ref}})`, in pure mechanics. For each calculation :math:`i`, the imposed load consists of the thermal deformation :math:`{\varepsilon }_{\mathit{xx}}\mathrm{=}\mathrm{-}{\varepsilon }_{\mathit{th}}\mathrm{=}\mathrm{-}\alpha (T)({T}_{i}\mathrm{-}{T}_{\mathit{réf}})`. The initial load consists of the deformations, stresses and internal variables from the previous mechanical calculation.

Indeed, for any behavior (assuming the additive decomposition of deformations):

:math:`{\sigma }_{\mathrm{xx}}=E(T)({\varepsilon }_{\mathrm{xx}}-{\varepsilon }^{\mathrm{th}}-{\varepsilon }_{\mathrm{xx}}^{p})`

in the first case, :math:`{\sigma }_{\mathrm{xx}}=E(T)(0-{\varepsilon }^{\mathrm{th}}-{\varepsilon }_{\mathrm{xx}}^{p})`, and in the second: :math:`{\sigma }_{\mathrm{xx}}=E(T)(\varepsilon -{\varepsilon }_{\mathrm{xx}}^{p})`.

It is therefore sufficient, at any moment, to apply :math:`{\varepsilon }_{\mathrm{xx}}=-{\varepsilon }^{\mathrm{th}}=-\alpha (T)(T-{T}_{\mathrm{ref}})` for mechanical calculation.

Moreover, in order to obtain the same results in both cases, it is necessary, at each time step of the second simulation, to perform the pure mechanical calculation with coefficients whose values are interpolated as a function of the temperature at the current moment. This interpolation is performed in the test command file, in a loop outside of STAT_NON_LINE.