Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
1/10
Organization (S): EDF/AMA, IPSN
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
V6.04.150 document
SSNV150 - Triaxial Traction with the law
of behavior BETON_DOUBLE_DP
Summary
This case of validation is intended to check the model of behavior 3D BETON_DOUBLE_DP formulated in
tally of the thermo plasticity, for the description of the nonlinear behavior of the concrete, in traction, and in
compression, with the taking into account of the irreversible variations of the thermal characteristics and
mechanics of the concrete, particularly sensitive at high temperature.
The description of cracking is treated within the framework of plasticity, using an energy equivalence,
by identifying the density of energy of cracking in mode I, with the plastic work of a homogeneous medium
equivalent, where the plastic deformation is uniformly distributed, in an “elementary” zone. This approach
preserve the continuity of the formulation of the model, on the whole of its behavior, and contributes to avoid
possible numerical difficulties during the change of state of material.
Pathological sensitivity of the numerical solution to the space discretization (grid), generated by
the introduction of a softening behavior of the concrete in traction and compression, is partially solved
by introducing an energy of cracking or rupture, dependant a characteristic length LLC, related to
cut elements.
The resolution of the equations constitutive of the model is carried out by an implicit scheme.
It is about a cube with 8 nodes subjected to a triaxial traction, in imposed displacement. This led loading
with the particular case of a hydrostatic state of stress, solved by projection at the top of the cone of traction,
when one places oneself in a hydrostatic diagram forced equivalent/forced. It is about a case test
with analytical solution.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
2/10
1
Problem of reference
1.1 Geometry
It is about a cube with 8 nodes, whose three faces have a normal displacement no one and the three faces
opposed an imposed and identical normal displacement.
The cube makes 1 mm on side. In modeling A, the cube is directed according to the Oxyz reference mark.
Modeling A
U2
Face1xy
Face1xz
Face1yz
U1
U
Faceyz
3
Ux = 0
Z
y
Uz = 0
N1
N2
X
Facexy
U
Facexz
2 = U1= U3
1.2
Material properties
To test the irreversible evolution of the mechanical characteristics with the temperature, one applies
a field of temperature decreasing. Certain variables depend on the temperature, others of
drying. Lastly, one applies a coefficient of withdrawal of desiccation not no one, equal to the coefficient of
thermal dilation, to test “data-processing” operation. The thermal deformations will be
thus equal and opposed to the deformations of withdrawal of desiccation. These dependences
intervene only for purely data-processing checks, the mechanical characteristics
can be regarded as constants.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
3/10
For the usual linear mechanical characteristics:
Young modulus:
E = 32.000 MPa
of
0°C with 20°C
E = 15.000 MPa
with
400°C (linear decrease)
E = 5.000 MPa
with
800°C (linear decrease)
Poisson's ratio:
N = 0.18
Thermal dilation coefficient:
has = 105/°C
Coefficient of withdrawal of desiccation: K = 105
For the nonlinear mechanical characteristics of model BETON_DOUBLE_DP:
Resistance in uniaxial pressing f' C = 40 NR/mm ²
of
0°C with 400°C
:
f' C = 15 NR/mm ²
with
800°C (linear decrease)
Resistance in uniaxial traction:
f' T = 4 NR/mm ²
of
0°C with 400°C
f' T = 1.5 NR/mm ²
with
800°C (linear decrease)
Report/ratio of resistances in compression
biaxial/uniaxial pressing:
B = 1.16
Energy of rupture in compression:
Gc =10 Nm/mm ²
Energy of rupture in traction:
WP =0.1 Nm/mm ²
Report/ratio of the limit elastic to resistance
in uniaxial pressing:
30%
1.3
Boundary conditions and loadings mechanical
Field of temperature decreasing of 20°C with 0°C.
Lower face of the cube (facexy):
blocked according to OZ.
Higher face of the cube (face1xy):
Uz displacement = 0,15 mm
Left face of the cube (faceyz):
blocked according to OX.
Right face of the cube (face1yz):
Ux displacement = 0,15 mm
Front face of the cube (facexz):
blocked according to OY.
Face postpones cube (face1xz):
Uy displacement = 0,15 mm
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
4/10
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution is calculated in an analytical way, knowing that in traction, only the criterion of
traction is activated, and that in the case of a hydrostatic loading, one is projected at the top of
cone of traction. It is thus necessary to solve a linear system of an equation to an unknown factor, which allows
to obtain the plastic deformation cumulated in traction. This one makes it possible to calculate then
strains and stresses.
2.2
Calculation of the reference solution of reference
For more detail on the notations and the setting in equation, one will refer to the reference document
[R7.01.03]. Only, the principal equations are pointed out here.
One notes “has”, imposed displacement following directions X, y and Z. The tensor of deformation is
form (has, has, has, 0., 0., 0.) by taking the usual notations of Code_Aster (three components
principal, three components of shearing).
The tensor of constraint is form (, 0., 0., 0.), in modeling A.
General equations of the model:
The equations constitutive of the model are written by distinguishing the isotropic part of the part
deviatoric of the tensors of constraints and deformations.
1
1
1
1
=
()
= -
=
~
= - tr
H
tr
S
tr () I
()
() I
3
3
H
tr
3
3
= S + I
= ~ +
H
and
H I
3
The equivalent constraint is written then: eq =
tr (s2)
2
In the case of an incremental formulation, and of a variable law of behavior, while noting with one
exhibitor “E” the elastic components of the constraint and the deformation, one obtains:
+
µ
K +
=
S +
+
2µ
~
E =
- +
+
-
3K
µ
and
H
-
H
H
K
The criteria in compression (fcomp) and traction (ftrac) are expressed in the following way:
+ A.
2
has
F
Oct.
Oct.
=
- F
eq
=
+ -
C (c)
F
comp
H
C (c)
B
B
3
B
+ C.
2
C
F
Oct.
Oct.
=
- F
eq
=
+ -
T (T)
F
trac
H
T (T)
D
D
3
D
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
5/10
tr (s2)
with Oct. =
3
tr ()
Oct. =
3
C: plastic multiplier in compression
T: plastic multiplier in traction
and has, B, C, D the coefficients of the model
The plastic deformations in traction and compression are expressed:
S
has
~
p
C
p
C =
H =
2 eq
B
C
C
B
3
S
C
~
p
T
p
T =
H =
2 eq
D
T
T
D
3
One obtains for the constraint:
S = -
+
2µ (~ p
~ p
E
+
p
p
C +
T)
=
- 3K (
C +
H
H
H
H T)
1
+
has
C
S
+
C
T
1 2µ
=
-
+
E
=
- 3K
+
H
H
C
T
B
2
2D eeq
B
3
D
3
for the equivalent constraint:
eq = eeq - 2µ +
C
T
+
B
2
2D
The two criteria lead then to a system of two equations to two unknown factors C and T with
to solve:
2
has
2 +
+ 2
µ
K has
2µ+ K +ac
E eq
E
+
-
+
-
+
- F - +
=
H
C
2
2
T
C (C
c)
0
3b
B
3b
B
3bd
data base
2
C
2µ + K +ac
2 +
+ 2
µ
K C
E eq
E
+
-
+
-
+
- F - +
=
H
C
T
2
2
T (T
T)
0
3D
D
3bd
data base
3D
D
In a similar way, in the case of the only criterion of traction activated, configuration of the case test, one
obtains a system of an equation to an unknown factor T to be solved:
2
C
2 +
+ 2
µ
K C
E eq
E
+
-
+
- F -
+
=
H
T
2
2
T (T
T)
0
3D
D
3D
D
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
6/10
Resolution with projection at the top of the cone of traction:
One thus seeks to solve this system, by using the particular shape of the tensors of constraints
and of deformations, uniforms on the structure.
On the basis of = (has, has, has, 0., 0., 0.) and of = (, 0., 0., 0.). one obtains:
= 3 +
X
has (2µ)
The elastic tensor of constraint
= 3 +
y
has (2µ)
= 3 +
Z
has (2µ)
S =
X
0
The elastic diverter of constraint S =
y
0
S =
Z
0
The hydrostatic constraint elastic E H = (
3 + 2µ) (A) = has
3 K
The elastic equivalent constraint eeq = 0
In the case of a curve of linear work hardening post-peak in traction, the expression of the parameter
of work hardening is as follows:
p
2.G ()
F
F
p
T
,
=, = 1
=
T (
T
) () F ()
T
() with ()
U
U
L. F ()
C
T
where the maximum of temperature during the history of loading indicates, F ()
T
resistance in traction.
~
p
C
p
+
C
T = 0
E
H
=
T
T
=
- 3K
D
3
H
H
T
D
3
The equation characterizing projection at the top of the cone of traction is as follows:
C
K +c2
L. F
E
-
- F
C
T
1 -
= 0
éq
2.2-1
D H
T
D 2
T
T
2.G
T
WP being the energy of rupture in traction (characteristic of material).
What makes it possible to obtain the plastic multiplier:
C E
C
- F
(aK
3
) -
H
T
ft
D
D
=
=
T
2
2
K +c2
L.
+ 2
.
C (F T)
K C
LLC (ft)
-
-
D 2
2.G
2
T
D
2.Gt
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
7/10
+
C
C
Then the constraint:
= E - 3K
K 3.a
H
H
T
T
D
3 =
-
D
Knowing has, imposed displacement, one obtains all the unknown factors of the problem.
2.3
Uncertainty on the solution
The solution being analytical, uncertainty is negligible, about the precision of the machine.
2.4 References
bibliographical
The model was defined starting from the following theses and is described in the report/ratio of specification:
[1]
G. Heinfling, at the time of its thesis “Contribution with the numerical modeling of the behavior of
concrete and of the concrete structures reinforced under mechanical thermo stresses with high
temperature ",
[2]
J.F. Georgin, at the time of its thesis “Contribution with the modeling of the concrete under stress of
fast dynamics. The taking into account of the effect speed by viscoplasticity ".
[3]
SCSA/128IQ1/RAP/00.034 Version 1.2, Développement of a model of behavior 3D
concrete with double criterion of plasticity in Code_Aster - Spécifications “.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
8/10
3 Modeling
With
3.1
Characteristics of modeling
3D (HEXA8)
1 element, stress field and uniform deformation.
U2
Face1xy
Face1xz
Face1yz
U1
U
Faceyz
3
Ux = 0
Z
y
Uz = 0
N1
N2
X
Facexy
Facexz
3.2
Characteristics of the grid
A number of nodes: 8
A number of meshs and type: 1 HEXA8
3.3
Functionalities tested
Commands Options
“MECHANICAL” AFFE_MODELE
“3D”
DEFI_MATERIAU “BETON_DOUBLE_DP”
DEFI_MATERIAU “ELAS_FO”
“K_DESSIC”
AFFE_CHAR_MECA “SECH_CALCULEE”
STAT_NON_LINE “BETON_DOUBLE_DP”
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
9/10
4
Results of modeling A
4.1 Values
tested
Components xx and zz of the stress field SIEF_ELNO_ELGA were tested, and
plastic deformation cumulated in traction (second variable internal, second component of
field VARI_ELNO_ELGA). Displacement being imposed, field EPSI_ELNO_DEPL is not
tested.
The three moments correspond to a displacement of 0.005, 0.01 and 0.015 Misters.
Component field SIEF_ELNO_ELGA SIXX
Identification Reference
Aster %
difference
For an imposed displacement
1.9182065 1.9182066 7.106
in load U1= U2= U3= 0.005
For an imposed displacement
1.161
1.1616769 3.106
in load U1= U2= U3= 0.010
6770
For an imposed displacement
0.4051470 0.4051473 7.105
in load U1= U2= U3= 0.015
Component field SIEF_ELNO_ELGA SIZZ
Identification Reference
Aster %
difference
For an imposed displacement
1.9182065 1.9182066 7.106
in load U1= U2= U3= 0.005
For an imposed displacement
1.1616770 1.1616769 3.106
in load U1= U2= U3= 0.010
For an imposed displacement
0.4051470 0.4051473 7.105
in load U1= U2= U3= 0.015
Component field VARI_ELNO_ELGA V2 (plastic deformation cumulated in traction)
Identification Reference
Aster %
difference
For an imposed displacement
0.0099232717 0.0099232717 4.107
in load U1= U2= U3= 0.005
For an imposed displacement
0.0199535329 0.0199535329 1.107
in load U1= U2= U3= 0.010
For an imposed displacement
0.0299837941 0.0299837941 3.108
in load U1= U2= U3= 0.015
4.2 Parameters
of execution
Version: 5.04.22
Machine: Claster
Obstruction memory: 32 Mo
Time CPU To use: 17.66 seconds
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
SSNV150 - Tractiontriaxiale with law BETON_DOUBLE_DP
Date:
22/01/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.150-A Page:
10/10
5
Summary of the results
This case test offers very satisfactory results compared to the analytical solution, lower than
7.10-5 % with a low iteration count (1 or 2 iterations). The solution is obtained from one
linear equation in the case of a linear curve of work hardening in traction, but the resolution uses
an algorithm of Newton within a more general framework.
One can note the work hardening of the criterion of traction which takes place during the loading, involving one
reduction in the constraint (component xx, yy and zz) in addition equalizes with the hydrostatic constraint.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Outline document