2. Benchmark solution#

2.1. Calculation method used for the reference solution#

The reference solution is calculated analytically, knowing that in traction, only the traction criterion is activated, and that in the case of hydrostatic loading, we project ourselves to the top of the traction cone. It is therefore necessary to solve a linear system of an equation with one unknown, which makes it possible to obtain the cumulative plastic deformation under tension. This makes it possible to then calculate deformations and stresses.

2.2. Calculation of the reference reference solution#

For more details on the notations and the equation, refer to the reference document [R7.01.03]. Only the main equations are mentioned here.

We note « \(a\) « , the displacement imposed in the directions \(x\), \(y\) and \(z\). The strain tensor is of the form \((a,a,a,0.,0.,0.)\) using the usual Code_Aster notations (three main components, three shear components).

The stress tensor is of the form \((\sigma ,\sigma ,\sigma ,0.,0.,0.)\), in modeling A.

General model equations:

The equations constituting the model are written by distinguishing the isotropic part from the deviatoric part of the stress and strain tensors.

_images/Object_1.svg _images/Object_2.svg _images/Object_3.svg _images/Object_4.svg _images/Object_5.svg

and

_images/Object_6.svg

The equivalent stress is then written as:

_images/Object_7.svg

In the case of an incremental formulation, and of a variable law of behavior, by noting with an exponent « e » the elastic components of stress and deformation, we obtain:

_images/Object_8.svg

and

_images/Object_9.svg

Compression criteria

_images/Object_10.svg

And in traction

_images/Object_11.svg

are expressed in the following way:

_images/Object_12.svg _images/Object_13.svg _images/Object_14.svg _images/Object_15.svg

: plastic multiplier in compression

_images/Object_16.svg

: plastic multiplier in traction

and

_images/Object_17.svg

The coefficients of the model

Plastic deformations in tension and compression are expressed:

_images/Object_18.svg _images/Object_19.svg
_images/Object_20.svg _images/Object_21.svg

For the constraint, we obtain:

_images/Object_22.svg _images/Object_23.svg
_images/Object_24.svg _images/Object_25.svg

for the equivalent stress:

_images/Object_26.svg

The two criteria then lead to a system of two equations with two unknowns.

_images/Object_27.svg

and

_images/Object_28.svg

to be solved:

_images/Object_29.svg

Analogously, in the case of the only activated traction criterion, configuration of the test case, we obtain a system of one equation with one unknown

_images/Object_30.svg

to be solved:

_images/Object_31.svg

Resolution with projection at the top of the traction cone:

We therefore seek to solve this system, by using the particular shape of the stress and deformation tensors, which are uniform on the structure.

Starting with \(\varepsilon \mathrm{=}(a,a,a,0.,0.,0.)\) and \(\sigma \mathrm{=}(\sigma ,\sigma ,\sigma ,0.,0.,0.)\), we get:

The elastic stress tensor

_images/Object_32.svg

The elastic stress diverter

_images/Object_33.svg

Elastic hydrostatic stress

_images/Object_34.svg

The equivalent elastic stress

_images/Object_35.svg

In the case of a post-peak linear strain hardening curve under tension, the expression of the work hardening parameter is as follows:

_images/Object_36.svg

with

_images/Object_37.svg

where

_images/Object_38.svg

refers to the maximum temperature during the loading history,

_images/Object_39.svg

tensile strength.

_images/Object_40.svg _images/Object_41.svg _images/Object_42.svg

The equation characterizing the projection at the top of the traction cone is as follows:

_images/Object_43.svg

Eq 2.2-1

_images/Object_44.svg

being the tensile rupture energy (characteristic of the material).

What makes it possible to obtain the plastic multiplier:

_images/Object_45.svg

Then the constraint:

_images/Object_46.svg

Knowing a, the imposed displacement, we get all the unknowns of the problem.

2.3. Uncertainty about the solution#

Since the solution is analytical, the uncertainty is negligible, in the order of the precision of the machine.

2.4. Bibliographical references#

The model was defined based on the following theses and is described in the specification report:

    1. Heinfling, during his thesis « Contribution to the numerical modeling of the behavior of concrete and reinforced concrete structures under thermo-mechanical stresses at high temperature »,

      1. Georgin, during his thesis « Contribution to the modeling of concrete under rapid dynamic stress. Taking into account the speed effect through viscoplasticity ».

SCSA /128 IQ1/RAP /00.034 Version 1.2, Development of a 3D concrete behavior model with double plasticity criteria in the*Code_Aster - Specifications ».