Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
1/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V6.04 booklet: Non-linear statics of the voluminal structures
Document: V6.04.162
SSNV162 - Calculation of creep of dessication
intrinsic with the model of Bazant
Summary:
This test makes it possible to validate the model of intrinsic creep of dessication of Bazant. This type of creep
start under the effect of the drying accompanied by a mechanical loading. The results of this test are
compared with the analytical solution for two plane forced modelings and 3D.
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
2/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
1
Problem of reference
1.1 Geometry
Z
X
y
L
H
E
height: H = 1m
width: L = 1 m
thickness: E = 1 m
1.2
Properties of material
6
36400E
E
=
AP
248
.
0
=
Here one informs also the curve sorption-desorption which connects the water content
C
with the hygroscopy
H
.
In this case one supposed that the numerical values of
C
and of
H
are the same ones.
Parameters specific to the intrinsic creep of dessication:
11
875
.
1
-
=
E
1
-
AP
1.3
Boundary conditions and loadings
In this test, one creates a homogeneous field of drying in the structure which varies with time.
At initial time 0s, the hygroscopy is worth 1 or 100% and at the final moment 10 days = 864000s is worth 0.75 or
75%.
Concerning the mechanical loading, one charges in compression of 0 à10 MPa in 1s
according to X on the right breakage (called BEFORE in the file of mesh), and one it maintains
constant loading during 10 days.
1.4 Conditions
initial
The beginning of calculation is supposed the moment 1. At this moment there is neither field of drying, nor forced
mechanics.
At moment 0, one applies a field of drying corresponding to 100% of hygroscopy.
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
3/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
2
Reference solution
2.1
Method of calculation
The equation of the model is:
H
fl
as of
&
&
=
The integration of this equation between moments 0 and T gives:
()
()
() ()
(
)
0
0
H
T
H
T
fl
as of
fl
as of
-
=
-
is a fact of the case. The loading is maintained constant, therefore
is constant with all
moments.
()
0
is null and
()
0
H
1 is worth. One can deduce
()
T
fl
as of
at every moment of calculation.
The elastic component of the field of deformation remains constant during all calculation:
E
T
E
=
)
(
and the total deflection is:
fl
as of
E
+
=
2.2
Sizes and results of reference
The test is homogeneous. One tests the deformation in an unspecified node.
2.3
Uncertainties on the solution
Exact analytical result.
2.4 References
bibliographical
[1]
L. GRANGER: Behavior differed from the concrete in the chambers of nuclear thermal power stations:
analyze and modeling. Thesis of doctorate of the ENPC (1995).
[2]
F. BENBOUDJEMA, F. MEFTAH, J.M. TORRENTI, Y. LE-PAPE: Algorithm of the model of
clean creep and of dessication. UMLV coupled to an elastic model. Note HS-DG/02/?.
EDF (to be appeared)
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
4/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
3 Modeling
With
3.1
Characteristics of modeling
Modeling forced plane.
X
y
With
B
C
N1
N3
N5
N7
N2
N4
N6
N8
D
The boundary conditions and the loading are modelized by:
FACE_IMPO and PRES_REP respectively.
One imposes a field of initial drying and a field of final drying which are homogeneous.
passage of the field of drying to the hygroscopy is done with the curve sorption-desorption (given
user).
3.2
Characteristics of the mesh
A number of nodes: 8
A number of meshs: 1 of type QUAD 8
4 of type SEG3
The following meshs are defined: AB, BC, CD and DA
To represent ¼ structure, one puts the boundary conditions following:
On AB: DY = 0
On AD: DX = 0
3.3 Functionalities
tested
Controls
Key word
DEFI_MATERIAU ELAS_FO
BAZANT_FD
FONC_DESORP
LAM_VISC
CREA_CHAMP “AFFE”
NOM_CMP=' TEMP'
“TEMP_INF”
“TEMP_SUP”
TYPE_CHAM=' NOEU_TEMP_R'
CREA_RESU AFFE
TYPE_RESU=' EVOL_THER'
CHAM_GD
AFFE_CHAR_MECA FACE_IMPO
PRES_REP
SECH_CALCULEE
STAT_NON_LINE COMP_INCR
RELATION= `BAZANT_FD'
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
5/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
3.4
Sizes tested and results
The component
xx
with the N4 node was tested
Moment Reference Aster %
difference
0. 0.
0.
-
1.
2.74725 E4
2.747253018 E4
1.E4
10.7999875 E+4
2.80584 E4
2.8058467007 E4
2.39 E4
53.9999375 E+4
3.04022 E4
3.0402214294 E4
4.70 E5
86400. E+4
3.12600 E4
3.2160024759 E4
7.70 E5
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
6/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
4 Modeling
B
4.1
Characteristics of modeling
Modeling 3D
Z
X
y
NO1
NO2
NO3
NO4
NO8
NO7
NO6
NO5
4.2
Characteristics of the mesh
A number of nodes:8
A number of meshs: 1 of type HEXA 8
6 of type QUAD 4
The following meshs are defined:
REAR NO1 NO3 NO7 NO5
FRONT
NO2 NO6 NO8 NO4
STRAIGHT LINE
NO1 NO5 NO6 NO2
LEFT NO3 NO4 NO8 NO7
LOW
NO1 NO2 NO4 NO3
HIGH
NO5 NO7 NO8 NO6
To account for the 1/8
ème
structure, the boundary conditions in displacement imposed are:
On nodes NO1, NO2, NO4 and NO3: DZ = 0
On nodes NO3, NO4, NO8 and NO7: DY = 0
On nodes NO1, NO3, NO5 and NO7: DX = 0
The loading is consisted of the same field of drying and the same pressure divided into
compression on the “FRONT” mesh.
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
7/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
4.3 Functionalities
tested
Moment Reference Aster %
difference
0. 0. 0.
-
1.
2.74725 E4
2.747253018 E4
1.E4
10.7999875 E+4
2.80584 E4
2.8058467007 E4
2.39 E4
53.9999375 E+4
3.04022 E4
3.0402214294 E4
4.70 E5
86400. E+4
3.12600 E4
3.2160024759 E4
7.70 E5
4.4
Sizes tested and results
The component
xx
with node NO6 was tested.
Moment Reference Aster %
difference
0. 0. 0.
-
1.
2.74725 E4
2.747253018 E4
1.E4
10.7999875 E+4
2.80584 E4
2.8058467007 E4
2.39 E4
53.9999375 E+4
3.04022 E4
3.0402214294 E4
4.70 E5
86400. E+4
3.12600 E4
3.2160024759 E4
7.70 E5
Code_Aster
®
Version
6.3
Titrate:
SSNV162 -
Calculation of intrinsic creep of dessication
Date:
12/12/02
Author (S):
J. EL GHARIB
Key
:
V6.04.162-B
Page:
8/8
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/02/001/A
5
Summary of the results
The values obtained with Code_Aster are in agreement with the values of the analytical solution of
reference. This same test was turned with Castem at laboratory LGCU at the University of the Marne
Valley the same results were obtained.