Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
1/14
Organization (S): EDF/AMA, CS IF
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
V6.04.143 document
SSNV143 - Biaxial 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 case test includes/understands two modelings 3D, the loading consists of a load followed by a discharge.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
2/14
1
Problem of reference
1.1 Geometry
It is about a cube with 8 nodes, whose two faces have a normal displacement no one, and the two faces
opposed have an imposed normal displacement, different one from the other of a coefficient 2.
The cube makes 1 mm on side. The cases tests are composed of a load, followed by a discharge. In
modeling A, the cube is directed according to the Oxyz reference mark. In modeling B, it is turned of 30°
by around axis OY.
U2
Face1xy
Modeling A
Face1yz
U1
Faceyz
Ux = 0
Z
y Uz = 0
N1
X
N2
Facexy
Modeling B
U2
Face1yz
U1
Face1xy
One = 0
Faceyz
Facexy
Z
NR
y
One = 0
2
X
NR
1
U2 = 2.U1
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
3/14
1.2
Material properties
To test the evolution of the mechanical characteristics in an irreversible way with the temperature, one
apply 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
thermal dilation coefficient, to test “data-processing” operation. Deformations
thermics thus equal and will be opposed to the deformations of withdrawal of desiccation. These
dependences intervene only for checks purely data-processing, the characteristics
mechanics can be regarded as constants.
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:
= 0.18
Thermal dilation coefficient:
= 10-5
Coefficient of withdrawal of desiccation: = 10-5
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 = 1.16
biaxial/uniaxial pressing:
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 30%
in uniaxial pressing:
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):
displacement 0.30 mm imposed followed by a discharge of
0.1 mm
Left face of the cube (faceyz):
blocked according to OX.
Right face of the cube (face1yz):
displacement 0.15 mm imposed followed by a discharge of
0.05 mm
Lower nodes front face (N1, N2):
blocked according to OY (Suppression of the movements of
solid body).
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
4/14
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution is calculated in an semi-analytical way, knowing that in traction, only the criterion
of traction is activated. It is thus necessary to solve a system of an equation to an unknown factor, which allows
to obtain by dichotomy for example, plastic deformation cumulated in traction. This one allows
to calculate strains and stresses then. This is possible, knowing displacement, and thus
deformation in the two imposed directions. Displacement in the third direction is then
an unknown factor of the problem.
The reference solution is calculated only in traction. The solution is determined by one
program resolution by dichotomy in independent FORTRAN. In compression, discharge,
exact solution was not recomputed, and constitutes a solution of nonregression of the code, been dependant on
version 5.02.14.
For modeling B, the results result by rotation from the tensor from constraint from
modeling A, of the intrinsic reference mark of the cube to the reference mark user, the stress field of both
configurations being identical in the intrinsic reference mark of the cube.
2.2
Calculation of the reference solution of reference
For more details on the notations and the setting in equation, one will refer to the document of
reference. Only, the principal equations are pointed out here.
One notes “has”, imposed displacement following direction X, and “2.a” following imposed displacement
direction Z. The tensor of deformation is form (has, y, 2.a, 0., 0., 0.) by taking the notations
usual of Code_Aster (three principal components, three components of shearing).
The tensor of constraint is form (X, 0., Z, 0., 0., 0.), in modeling A.
The criterion of traction is expressed in the form:
+ C.
2
C
F
Oct.
Oct.
=
- F
eq
=
+ -
T (T)
F
trac
H
T (T)
D
D
3
D
The constitutive equations are written by distinguishing the isotropic part of the deviatoric part of
tensors of constraints and deformations.
1
1
1
1
=
= -
=
~
= - tr
H
tr
H
tr () S
tr () I
()
() I
3
3
3
3
= S + I
H
= ~ + H I
3
The equivalent constraint is written then: eq =
tr ()
S
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
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
5/14
The criteria in compression and traction 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
The plastic deformations in traction and compression are expressed:
S
has
p
~
C
p
C =
H =
C
2 eq
B
C
B
3
S
C
p
~
T
p
T =
H =
T
2 eq
D
T
D
3
One obtains for the constraint:
S = -
+
2µ (~ p
~ p
E
+
p
p
C +
T)
=
- 3
+
H
H
K (H C
H T)
1
+
has
C
S
+
C
T
1 2µ
= -
+
E
=
- 3
+
H
H
K
C
T
B
2
2D eeq
B
3
D
3
eq
E eq
+
C
for the equivalent constraint:
=
- 2µ
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.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
6/14
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.
2 A
y
) .0
.,
0
.,
0
,
and of = (X.,
0
, Z
).
0
.,
0
.,
0
,
. one obtains:
X = has (+ 2µ) + y + 2 A
.
The elastic tensor of constraint
y = has
. + y (+ 2µ) + 2a
.
Z = has
. + y + 2a (. + 2µ)
2
S = -
µ
X
. y
3
4
The elastic diverter of constraint S = - 2.µ.a +
µ
y
. y
3
2
S = 2.µ.a -
µ
Z
. y
3
1
The hydrostatic constraint elastic E H = (
3 + 2µ) has + .y
3
The equivalent constraint elastic E
2
2
eq = µ 4
- 12. +
y
y 12.a has
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
,) = (,) = F
T
1
=
U
T
T
T ()
with ()
U ()
LLC. ft ()
One thus seeks to solve the equation:
2
C
2 +
+ 2
µ
K C
L. F
E eq
E
+
-
+
- F 1
C
T
-
= 0
éq
2.2-1
3
D H
T3 2
2
D
D
D
T
T
2.G
T
Knowing that the constraint in the direction is null there, one second equation is obtained:
+
1
T
E
E
=
+
= 0 = 1 - 2
+
y
sy H
µ
S y H
D eeq
1
+
.
T
4
C
E
T
= 0 = 1 - 2
-
y
µ
µ
2.µ.a
K
D eeq
y
H
3
+
-
D
4 µ
E
-
y
2.µ.a
H
3
+
From where: =
T
that one can substitute in the expression of the criterion
2µ 4 µ
- 2.µ
C
.a
K
D. eeq
y
3
+ D
[éq 2.2-1].
Knowing has, imposed displacement, one obtains a nonlinear equation with an unknown factor, that one
can solve simply by dichotomy, and which makes it possible to calculate the deformation y, then the unit
unknown factors of the system.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
7/14
2.3
Uncertainty on the solution
It is negligible, about the precision machine.
2.4 References
bibliographical
The model was defined starting from the following theses:
[1]
J.F. GEORGIN, at the time of its thesis “Contribution with the numerical modeling of the behavior
concrete concrete and structures reinforced under mechanical thermo stresses with high
temperature ",
[2]
G. HEINFLING, 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 ", and is described
in the report/ratio of specification:
[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.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
8/14
3 Modeling
With
3.1
Characteristics of modeling
3D (HEXA8)
1 element, stress field and uniform deformation.
U2
U1
Ux = 0
Z
y
Uz = 0
X
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
AFFE_MODELE
“MECANIQUE”
“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.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
9/14
4
Results of modeling A
4.1 Values
tested
The nonnull components of the stress field SIEF_ELNO_ELGA (component were tested
xx and zz), the component yy of the field of deformation EPSI_ELNO_DEPL, which constitutes an unknown factor
system (deformations in the two other directions being imposed), deformation
figure cumulated in traction (second variable internal, second component of the field
VARI_ELNO_ELGA), and finally, only for the fourth case of loading (discharge),
plastic deformation cumulated in compression, (first internal variable, first component of
field VARI_ELNO_ELGA).
The first three loadings correspond to the load, and have results of reference.
The fourth loading corresponds to the discharge, and constitutes a result of nonregression of
code.
Component field SIEF_ELNO_ELGA SIXX
Identification Reference
Aster
% difference
For a displacement imposed in
0.1235611 0.1235380 0.019
charge U1=0.1 and U2= 0.05
For a displacement imposed in
6.882374.102 6.878218.102 0.060
charge U1=0.2 and U2= 0.10
For a displacement imposed in
1.408764.102 1.402639.102 0.435
charge U1=0.3 and U2= 0.15
For a displacement imposed in
(*) 4.195092.105 -
discharge U1=0.1 and U2= 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Component field SIEF_ELNO_ELGA SIZZ
Identification Reference
Aster
% difference
For a displacement imposed in
0.239212 0.239174 0.016
charge U1=0.1 and U2 = 0.05
For a displacement imposed in
0.133243 0.133165 0.059
charge U1=0.2 and U2 = 0.10
For a displacement imposed in
2.727403.102 2.725569.102 0.434
charge U1=0.3 and U2 = 0.15
For a displacement imposed in
(*) 4.959258.105 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
10/14
Component field EPSI_ELNO_DEPL EPYY
Identification Reference
Aster
% difference
For a displacement imposed in
3.419463.103 3.419464.103 2.107
charge U1=0.1 and U2 = 0.05
For a displacement imposed in
6.835813.103 6.835815.103 2.107
charge U1=0.2 and U2 = 0.10
For a displacement imposed in
1.025216.102 1.025216.102 2.107
charge U1=0.3 and U2 = 0.15
For a displacement imposed in
(*) 4.357498.101 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Component field VARI_ELNO_ELGA VARI_2 (plastic deformation cumulated in traction)
Identification Reference
Aster
% difference
For an imposed displacement
1.085728.102 1.085728.102 5.109
U1=0.1 and U2 = 0.05
For an imposed displacement
2.171556.102 2.171556.102 5.109
U1=0.2 and U2 = 0.10
For an imposed displacement
3.257385.102 3.257385.102 4.109
U1=0.3 and U2 = 0.15
For an imposed displacement
3.257385.102 3.257385.102 4.109
U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Component field VARI_ELNO_ELGA VARI_1 (plastic deformation cumulated in
compression)
Identification Reference
Aster
% difference
For a displacement imposed in
(*) 3.528401.101 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
11/14
5 Modeling
B
5.1
Characteristics of modeling
3D (HEXA8)
1 element, stress field and uniform deformation.
U2
U
1
One = 0
Z
One = 0
y
X
5.2
Characteristics of the grid
A number of nodes: 8
A number of meshs and type: 1 HEXA8
5.3
Functionalities tested
Commands Options
AFFE_MODELE
“MECANIQUE”
“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.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
12/14
6
Results of modeling B
The nonnull components of the stress field SIEF_ELNO_ELGA (component were tested
xx, zz and xz), plastic deformation cumulated in traction (second variable internal, second
component of field VARI_ELNO_ELGA), and finally, only for the fourth case of loading
(discharge), plastic deformation cumulated in compression, (first internal variable, first
component of field VARI_ELNO_ELGA).
The first three loadings correspond to the load, and have results of reference.
The fourth loading corresponds to the discharge, and constitutes a result of nonregression of
code.
6.1 Values
tested
Component field SIEF_ELNO_ELGA SIXX
Identification Reference
Aster
% difference
For a displacement imposed in
0.152474 0.1524472 0.018
charge U1=0.1 and U2 = 0.05
For a displacement imposed in
8.492877.102 8.487797.102 0.060
charge U1=0.2 and U2 = 0.10
For a displacement imposed in
1.732484.102 1.730871.102 0.434
charge U1=0.3 and U2 = 0.15
For a displacement imposed in
(*) 4.386134.105 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Component field SIEF_ELNO_ELGA SIZZ
Identification Reference
Aster
% difference
For a displacement imposed in
0.210300 0.210265 0.016
charge U1=0.1 and U2 = 0.05
For a displacement imposed in
0.117138 0.117069 0.059
charge U1=0.2 and U2 = 0.10
For a displacement imposed in
2.397743.102 2.387336.102 0.434
charge U1=0.3 and U2 = 0.15
For a displacement imposed in
(*) 4.768217.105 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
13/14
Component field SIEF_ELNO_ELGA SIXZ
Identification Reference
Aster
% difference
For a displacement imposed in
5.007871.102 5.007226.102 0.013
charge U1=0.1 and U2 = 0.05
For a displacement imposed in
2.789472.102 2.787871.102 0.057
charge U1=0.2 and U2 = 0.10
For a displacement imposed in
5.709873.103 5.685155.103 0.433
charge U1=0.3 and U2 = 0.15
For a displacement imposed in
(*) 3.308936.106 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no analytical solution.
Component field VARI_ELNO_ELGA VARI_2 (plastic deformation cumulated in traction)
Identification Reference
Aster
% difference
For a displacement imposed in
1.085728.102 1.085728.102 5.109
charge U1=0.1 and U2 = 0.05
For a displacement imposed in
2.171556.102 2.171556.102 5.109
charge U1=0.2 and U2 = 0.10
For a displacement imposed in
3.257385.102 3.257385.102 4.109
charge U1=0.3 and U2 = 0.15
For a displacement imposed in
3.257385.102 3.257385.102 4.109
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Component field VARI_ELNO_ELGA VARI_1 (plastic deformation cumulated in
compression)
Identification Reference
Aster
% difference
For a displacement imposed in
(*) 3.528401.101 -
discharge U1=0.1 and U2 = 0.05
(*) discharges some, one carries out a test of nonregression. There is no calculated analytical solution.
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Code_Aster ®
Version
5.2
Titrate:
SSNV143 - Biaxial Traction with law BETON_DOUBLE_DP
Date
:
23/09/02
Author (S):
C. CHAVANT, B. CIREE Key
:
V6.04.143-A Page:
14/14
7
Summary of the results
This case test offers satisfactory results compared to the results of reference, lower than 0.06%
for the first two cases of loading, more important for the third, which is explained by one
relatively low level of constraint (one reaches the end of the curve of work hardening in traction).
The test discharges some (fourth loading) makes it possible to check nonthe regression of the code.
The iteration count is relatively important with the first step of calculation, about 13, then
drop to 7, 4 and 1, which is explained by the passage of the plastic threshold to the first step of calculation, for
to reach a quasi linear behavior thereafter (curved post-peak linear).
One obtains also a more significant number of iterations to step 31 (beginning of the fourth case of
loading), then an iteration count dropping up to 1, because of the passage in discharge, with one
change of behavior, follow-up of a quasi linear behavior thereafter (curved post-peak
linear).
Handbook of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Outline document