2. Benchmark solution#

The solution to the problem posed can be determined analytically.

We note:

  • \({\varepsilon }_{0}\) the deformation applied in the \(z\) direction,

  • \({\varepsilon }_{1}\), \({\varepsilon }_{2}\), and \({\varepsilon }_{3}\) the main deformations

2.1. First loading step: simple compression#

  • The deformation tensor is equal to:

    _images/Object_6.svg

with \({\varepsilon }_{0}<0\)

  • The equivalent deformation is therefore equal to:

    _images/Object_7.svg
  • As soon as

    _images/Object_8.svg

, there is an evolution of the damage which is equal to:

_images/Object_10.svg
  • Finally, constraint \({\sigma }_{\mathit{zz}}\) is equal to:

_images/Object_11.svg

2.2. Second loading stage: thermal expansion in plane deformations#

  • The tensor of the total deformations is equal to:

_images/Object_12.svg

with fixed \({\varepsilon }_{0}<0\)

  • The elastic deformation equal

    _images/Object_13.svg

, the equivalent deformation is equal to:

_images/Object_14.svg
  • The damage is worth:

_images/Object_16.svg
  • Finally, constraint \({\sigma }_{\mathit{zz}}\) is equal to:

_images/Object_17.svg

Note:

    • In a given state, the material parameters used are those defined at the maximum temperature seen by the material and not at the current temperature.

    • The assessment of damage \(D\) involves the concept of maximum reached during the history of the load; the solution is therefore not completely analytical but involves discretization. In the case where there is no thermal influence, simply take \(\stackrel{~}{\mathrm{\epsilon }}\) equivalent to the maximum equivalent deformation reached. When the thermal aspect is taken into account, heating can contribute to « reducing » or « delaying » damage at a given deformation; this is the case with the evolution of \({B}_{c}\) retained. In this case, it is in fact necessary to discretize the load quite finely to have the right damage value \(D\) (which in fact has a maximum in our case) .