Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
1/8
Organization (S): EDF-R & D/AMA
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of structures
voluminal
Document: V7.31.124
WTNV124 - Essai of désaturation-consolidation
with the model of Barcelona
Summary:
This test makes it possible to validate the model of Barcelona, which integrates an elastoplastic mechanical law coupled to
hydraulics (and possibly with thermics) in condition of nonsaturation of the liquid phase. This law
integrate an elastoplastic hydrostatic mechanism (of which the elastic part is non-linear and the threshold
of flow corresponds to a pressure of variable consolidation with suction) coupled to a mechanism
deviatoric elastoplastic. The characteristics of these mechanisms depend on suction (i.e.
difference between gas pressure and pressure of liquid). There are in particular two mechanisms
of work hardening in completely coupled pressure and suction. The surface of load of the model of Barcelona
present (in the diagram pressure hydrostatic-diverter and for a given suction) in the form
of an ellipse cutting the hydrostatic axis in two points: the value of the pressure of consolidation and
cohesion of material proportional to suction. In condition of complete saturation, this criterion is reduced to
that of the Cam_Clay model specific to the saturated normally consolidated grounds.
This test carried out in hydro-mechanical coupling (modeling HHM) includes/understands two ways of
loading:
1) a way of desaturation while making increase the capillary pressure beyond the threshold of plasticity
hydrous.
2) a hydrostatic way of compression on the désaturé sample.
Modeling is carried out in 3D.
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
2/8
1
Problem of reference
1.1 Geometry
Z
E
H
y
X
L
height: H = 1m
width: L = 1 m
thickness: E = 1 m
1.2
Properties of material
Elastic thermo properties:
E = 22.4E7 Pa
= 0.3
3
= 2500 kg/m
Parameters specific to the model of CAM_CLAY:
·
Initial porosity PORO = 0.14
·
Modulate plastic compressibility with of saturated = 0.25,
·
Modulate elastic compressibility = 0.05,
·
Slope of the straight line criticizes M = 0.9,
·
Critical pressure equalizes with half of the pressure of consolidation to saturation
CLOSE _ CRIT = E
.
3 7Pa,
·
Pressure of reference Pa = 1.E5 Pa
Additional parameters specific to BARCELONA:
·
Parameters allowing to calculate the module of compressibility according to suction
(p) = (
C
[)
0 1
(- R) exp (- p
) + R,
C
]
R =
,
75
.
0
=
5
.
12th - 6
·
Slope of cohesion K =
6
.
0
C
·
Initial threshold of suction PC0_INIT p
()
0 = 6 7
c0
E
·
Modulate elastic compressibility of suction =
01
.
0
S
·
Modulate plastic compressibility of suction =
05
.
0
S
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
3/8
Hydraulic properties: the hydraulic properties of material which are independent of the model
from Barcelona but nevertheless necessary to carry out coupled calculation are presented in
table below:
Liquid water
Density (kg.m-3)
1.103
Heat with constant pressure (J.K-1)
4180
Opposite of the coefficient of compressibility (Pa-1)
0.510-9
Intrinsic permeability (m2)
1.10-18
viscosity
1.10-3
Initial State
Porosity
0,14
Temperature
293°K
Capillary pressure
1.51 107
Gas pressure
1.51 105
Initial saturation in liquid
0,99
Constants
Constant of perfect gases
8,315
Coefficients
Homogenized density
2400
homogenized
Capillary curve
S (
9
-
C
P) = 0 99
. (1 - 6 10 PC)
Coefficient of Biot
1
1.3
Boundary conditions and loadings
The first way of loading consists in carrying out a desaturation, with a pressure
hydrostatic constant P = E
.
1 7 Pa (the gas pressure is maintained constant
with PRE2 = E
.
1 5 Pa during all the test). The capillary pressure varies from zero until
PRE1 = E
.
7 7 Pa, beyond the threshold of plasticization equal initially to PRE1 = E
.
6 7 Pa.
Work hardening due to the increase in suction causes an increase in the threshold of consolidation,
since two work hardenings hydrous and mechanical are coupled in the model of Barcelona. One
checks while following the second way consisting in exerting a hydrostatic pressure exceeding it
initial threshold (P
= E
.
6
Pa
idiot
7
) without causing plasticization.
1.4 Conditions
initial
The initial constraint (forced effective of Bishop) is selected in such way that the constraint used
in the behavior (
D
= T + pgz1) is inside the surface of load.
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
4/8
2
Reference solution
In the absence of loading deviatoric, an exact solution is available for the deformations and
thresholds of work hardening at all the stages of the loading:
P
Reversible voluminal deformation in mechanical loading =
Ln
v
1+ e0
0
P
S
PC + p
Reversible voluminal deformation in hydrous loading
atm
v =
Ln
1+ E
p
0
atm
Total voluminal deformation in hydrous loading, after crossing of the threshold
:
+
S
PC + patm
=
Ln
if p > p
v
-
C
0
1
(
C
+ e0)
PC + patm
Total voluminal deformation in mechanical loading, after crossing of the threshold of
(p) dP
consolidation: D
C
v =
1+ E
P
0
dp 0
- dP
Coupling of the thresholds:
C
Cr
=
PC + p
0
atm
S - S Cr
P
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
5/8
3 Modeling
With
3.1
Characteristics of modeling
Modeling 3D
Z
NO3
NO5
NO8
NO2
NO7
y
NO1
NO4
NO6
NO12
NO10
NO9
NO11
X
NO15
NO17
NO20
NO14
NO19
NO13
NO18
NO16
3.2
Characteristics of the grid
A number of nodes:
20
A number of meshs:
1 of type HEXA 20
6 of type QUAD 8
The following meshs are defined:
DROITE
NO3 NO5 NO8 NO10 NO12 NO15 NO17 NO20
GAUCHE
NO1 NO4 NO6 NO9 NO11 NO13 NO16 NO18
DEVANT
NO6 NO7 NO8 NO11 NO12 NO18 NO19 NO20
DERRIERE
NO1 NO2 NO3 NO9 NO10 NO13 NO14 NO15
BAS
NO13 NO14 NO15 NO16 NO17 NO18 NO19 NO20
HAUT
NO1 NO2 NO3 NO4 NO5 NO6 NO7 NO8
To represent the 1/8th structure, the boundary conditions in displacement imposed are:
On face BAS: DZ = 0
On face GAUCHE: DY = 0
On face DERRIERE: DX = 0
The loading is consisted of the same pressure divided into compression on the 3 meshs: `HAUT,
“DROITE” and “DEVANT” to simulate a hydrostatic test. All the nodes are compelled with one
constant pressure of gas and with a suction varying of 0 with 7.107Pa
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
6/8
3.3 Functionalities
tested
Commands
DEFI_MATERIAU BARCELONA
STAT_NON_LINE COMP_INCR
RELATION= `KIT_HHM'
`RELATION_KIT'=
“BARCELONA”
“LIQU_SATU”
“HYDR_UTIL”
NEWTON
STAMP = “TANGENT”
3.4
Sizes tested and results
It is about a homogeneous test, the place of observation of the fields is indifferent. It will be tested
uz displacement with node 8 at moment 1 (end of the hydrostatic way), at moment 6.0 (crossing of
hydrous threshold) as well as the variables intern of indicator of plasticity and pressure criticizes with same
node. One tests then same the fields at moment 10 (loading purely hydrostatic in
crossing the old threshold thus without plasticization) and finally at the moment 20 where the new one is crossed
mechanical threshold (plasticization).
C
p
6
10
20
1
C
p
C
p 0
C
L 0
C
L
1
C
L
2
1
P
Values of uz:
Moment
Reference
Aster
1st loading
1. - 5.607E-02
- 5.593E-02
2nd loading
6.0 - 7.565E-02
- 7.647E-02
3rd loading
10.0 - 1.3097E-01
- 1.3097E-01
4th loading
20. - 2.11952E-01
- 2.11951E-01
Plastic indicator (mechanical threshold):
Moment
Reference
Aster
1st loading
1. 0
0
2nd loading
6.0 0
0
3rd loading
10.0 0
0
4th loading
20. 1
1
Hydrous indicator of irreversibility:
Moment
Reference
Aster
1st loading
1. 0
0
2nd loading
6.0 1
1
3rd loading
10.0 0
0
4th loading
20. 0
0
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
7/8
Value of the hydrous threshold:
Moment
Reference
Aster
1st loading
1. 6e+07
6e+07
2nd loading
6.0 7e+07
7e+07
3rd loading
10.0 7e+07
7e+07
4th loading
20. 1.108E+10
1.108E+10
Value of the mechanical threshold:
Moment
Reference
Aster
1st loading
1. 3e+07
3e+07
2nd loading
6.0 5.6661E+08
5.7459E+08
3rd loading
10.0 5.74596E+08
5.74596E+08
4th loading
20. 25.0631E+08
25.0630E+08
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Code_Aster ®
Version
7.3
Titrate:
WTNV124 Essai of désaturation-consolidation: model of BARCELONA Date:
27/09/04
Author (S):
G. DEBRUYNE, J. EL GHARIB Key
:
V7.31.124-A Page:
8/8
4
Summary of the results
One obtains results on displacements which present a variation compared to the solution
analytical of about 0.2% in the reversible field and of 1% in the plastic range.
even report can be made on the values of the thresholds of consolidation.
Handbook of Validation
V7.31 booklet: Thermo hydro-mechanical in porous environment of voluminal structures
HT-66/04/005/A
Outline document