2. Benchmark solution#

2.1. Calculation method used for the reference solution for modeling A#

It is a flat problem. We can look for the solution in the form of a rigid rotation followed by an expansion by a factor \(a\) in one direction and \(b\) in the other:

_images/Object_4.svg

The transformation gradient and the Green-Lagrange deformation are then:

_images/Object_5.svg

The behavioral relationship leads to a diagonal Lagrangian stress tensor (with \(\lambda\) and \(\mu\) the Lamé coefficients):

_images/Object_7.svg

We deduce the Cauchy stress tensor, which is also diagonal:

_images/Object_8.svg

Finally, the boundary conditions are written as:

_images/Object_9.svg

It is also possible to calculate the forces exerted on the faces:

_images/Object_10.svg

where

_images/Object_11.svg

represent the initial surfaces of the faces.

2.2. Model B reference#

The reference results for this modeling are provided by the same calculation to which the rotation step was not applied.

2.3. Benchmark results#

As reference results, we adopt the displacements, the Grenn-Lagrange deformations, the Cauchy stresses and the forces exerted on faces \([\mathrm{1,}3]\), \([\mathrm{3,}4]\) and \([\mathrm{1,}\mathrm{2,}\mathrm{3,}4]\) at the end of loading (\(t=2s\)).

We’re looking for

_images/Object_12.svg

such as dilation

_images/Object_13.svg

be \(p=–26610.3\mathrm{MPa}\).

Dilation \(b\) and displacements are then:

_images/Object_16.svg

The Cauchy constraints are as follows:

_images/Object_17.svg

Finally, the forces exerted are:

_images/Object_18.svg

2.4. Uncertainty about the solution#

Analytical solution.

2.5. Bibliographical references#

  1. Eric LORENTZ « A non-linear hyperelastic behavior relationship » Internal Note EDF/DER HI‑74/95/011/0