background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
1/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
Organisme(s) :
EDF-R&D/AMA, CNEPE















Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
Document : V7.31.100



WTNV100 - Essai triaxial non drainé avec le modèle
CJS (niveau 1)



Résumé

Ce test permet de valider le niveau 1 du modèle CJS. Il s'agit d'un essai triaxial en condition non drainée.
Dans les deux premières modélisations, les calculs sont effectués uniquement sur la partie solide du sol,
l'aspect non drainé étant modélisé par une déformation volumique nulle du squelette, ce sont des modélisations
3D qui diffèrent l'une de l'autre uniquement par le maillage.
Dans la troisième modélisation, le couplage hydraulique est pris en compte, l'échantillon est totalement saturé,
le squelette et le fluide sont supposés incompressibles.
Par raison de symétrie, on ne s'intéresse qu'au huitième d'un échantillon soumis à un essai triaxial.
Le niveau de confinement est de 100 kPa.
Les résultats obtenus avec le modèle CJS1 sont comparés à une solution analytique.
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
2/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
1
Problème de référence
1.1 Géométrie
z
x
y
e
l
h
C
A
B
hauteur :
h = 1 m
largeur :
l = 1 m
épaisseur : e = 1 m
Coordonnées des points (en mètres) :
A
B
C
x 0. 0.
0.5
y 0. 1.
0.5
z 0. 0.
0.5
1.2
Propriété de matériaux
E = 22,4 10
3
kPa
= 0,3
Coefficient de biot b : 1
L'eau est supposée incompressible : UN_SUR_K = 0
Paramètres CJS1 :
= - 0,03
= 0,82
R
m
= 0,289
P
a
= -100 kPa
1.3
Conditions initiales, conditions aux limites, et chargement
1.3.1 Modélisation mécanique pure
Phase 1 :
On amène l'échantillon à un état homogène :
xx
yy
zz
0
0
0
=
=
, en imposant la pression de
confinement correspondante sur les faces avant, latérale droite et supérieure. Les déplacements sont
bloqués sur les faces arrière (
u
x
=
0
), latérale gauche (
u
y
=
0
) et inférieure (
u
z
=
0
).
Phase 2 :
On maintient les déplacements bloqués sur les faces arrière (
u
x
=
0
), latérale gauche (
u
y
=
0
) et
inférieure (
u
z
=
0
). On applique un déplacement imposé sur la face supérieure :
( )
u t
z
, de façon à
obtenir une déformation
zz
= -
20 %
(comptée à partir du début de la phase 2). Sur les faces avant
et latérale droite, on impose respectivement les déplacements
( )
u t
x
et
( )
u t
y
, de façon à avoir une
déformation volumique nulle pour l'échantillon, c'est à dire finalement que l'on impose
xx
yy
zz
=
= -
2
. C'est la manière de reproduire le comportement de la phase solide lors d'un essai
triaxial non drainé.
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
3/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
1.3.2 Modélisation couplée avec hydraulique
Phase 1 :
On amène l'échantillon à un état homogène de contraintes effectives:
xx
yy
zz
0
0
0
=
=
, en imposant
la pression totale correspondante sur les faces avant, latérale droite et supérieure et en imposant
partout des pressions d'eau nulles. Les déplacements sont bloqués sur les faces arrière (
u
x
=
0
),
latérale gauche (
u
y
=
0
) et inférieure (
u
z
=
0
).
Phase 2 :
On maintient les déplacements bloqués sur les faces arrière (
u
x
=
0
), latérale gauche (
u
y
=
0
) et
inférieure (
u
z
=
0
).
Sur toutes les faces, les flux hydrauliques sont nuls.
On applique un déplacement imposé sur la face supérieure de façon à obtenir une déformation
zz
= -
20 %
(comptée à partir du début de la phase 2). Sur les faces avant et latérale droite, on
impose des conditions aux limites en contrainte totale :
.
(
)
n
kPa
=
=
0
100
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
4/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
2
Solution de référence
2.1
Solution de référence pour la pression d'eau en linéaire
0
,
0
p
0
désignant les contraintes ,déformations et pressions d'eau obtenues dans la phase on a :
(
)
(
) (
)
(
)
µ
-
=
-
+
-
-
-
= -
+
-
0
0
0
0
0
0
2
tr
b p p
m
p p
M
btr
fl
Dans cette écriture,
M
désigne le module de biot et
N
M
=
1
.
Les conditions aux limites de flux nul et la conservation de la masse d'eau donnent
m
=
0
Les conditions aux limites sur les parois latérales et le fait que l'état de contrainte est homogène
donnent :
xx
xx
-
=
0
0
On a donc finalement à résoudre les deux équations :
(
)
(
)
µ
2
2
2
xx
zz
xx
xx
zz
bp
b
p
M
Np
+
+
=
+
= -
= -



Et on obtient :
(
)
(
)
µ
µ
µ
xx
zz
zz
b
N
b
N
p
b
b
N
= -
+
+
+
= -
+
+


2
2
2
2
Dans notre cas,
µ
xx
zz
zz
p
= -
= -
2
;
2.2
Développement de la solution analytique CJS
On a en permanence :
pour les déformations :
xx
yy
zz
=
= -
2
pour les contraintes :
xx
yy
=
Phase élastique :
En écrivant simplement la loi élastique, il vient :
µ
xx
xx
zz
=
-
0
µ
zz
zz
zz
=
+
0
2
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
5/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
Par ailleurs, on sait aussi que pendant cette phase
I
1
(= trace(
)) reste constant car
v
=
0
. On en
déduit pour les composantes du déviateur :
s
I
I
xx
xx
xx
zz
zz
=
-
=
-
-
= -
µ
µ
1
0
10
3
3
et
s
I
I
zz
zz
zz
zz
zz
=
-
=
-
+
=
µ
µ
1
0
10
3
3
2
2
soit :
s
II
zz
= -
6
µ
et
( )
det s
zz
=
2
3
3
µ
Par conséquent :
( )
( )
h
s
= -
1
1 6
/
Ainsi lorsqu'on atteint le critère
f
d
=
0
, on a :
(
)
(
)
s
R I
R I
II
m
zz
m
1
6
1
0
1 6
10
1 6
10
-
+
= -
-
+
=
µ
/
/
C'est-à-dire que la transition entre les états élastique et parfaitement plastique se fait pour une
déformation axiale égale à :
(
)
µ
zz
trans
m
R I
=
-
10
1 6
6
1
/
L'état de contraintes correspondant est noté :
(
)
µ
µ
xx
trans
xx
m
R I
=
-
-
0
10
1 6
6
1
/
et
(
)
µ
µ
zz
trans
zz
m
R I
=
+
-
0
10
1 6
2
6
1
/
Phase plastique :
On note
s
d
-
le déviateur de l'inverse du tenseur s
D'une façon générale, on a les grandeurs suivantes :
(
)
s
s
xx
zz
xx
yy
= -
-
=
1
3
s
xx
zz
xx
-
=
-
-
1
3
(
)
s
xxd
zz
xx
-
=
-
-
3
2
(
)
s
zz
zz
xx
=
-
2
3
(
)
s
zz
zz
xx
-
=
-
1
3
2
s
zzd
zz
xx
-
=
-
3
soit :
(
)
s
II
zz
xx
= -
-
2
3
et
( )
(
)
det s
zz
xx
=
-
2
3
3
3
Par conséquent :
( )
(
)
h
s
= -
1
1 6
/
On en déduit :
( )
Q
xx
=
-
1
6 1
1 6
/
et
( )
Q
zz
= -
-
2
3 1
1 6
/
en outre :
(
)
f
R
d
xx
m
=
-
+
1
6 1
1 6
/
et
(
)

f
R
d
zz
m
= -
-
+
2
3 1
1 6
/
Comme on a :
( )
=
-




=
-




=
signe s
s
s
R
R
ij
ij
II
II
c
m
c
&
1
1
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
6/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
alors le tenseur
n
s'écrit :
n
xx
=
+
+




1
3
1
6
1
2
et
n
zz
=
+
-
+




1
3
2
3
1
2
Il vient alors pour
G
d
:
( )
( )
G
R
R
xx
d
m
m
=
-
+
-
-
+
+
+




1
6 1
1
3
3
1
6
1
1 6
1 6
2
/
/
( )
( )
G
R
R
zz
d
m
m
= -
-
+
-
-
+
+
-
+




2
3 1
1
3
3
2
3
1
1 6
1 6
2
/
/
On a également d'après la loi élastique :
xx
xx
trans
xx
=
+
et
zz
zz
trans
zz
=
+
où :
(
)
( )
(
)
(
)
µ
µ
µ
xx
xx
d
xx
d
v
d
d
zz
d
xx
d
d
xx
d
zz
d
G
tr G
G
G
G
=
-
+
-
= -
-
-
+
2
2
2
(
)
( )
(
)
(
)
µ
µ
µ
zz
zz
d
zz
d
v
d
d
zz
d
zz
d
d
xx
d
zz
d
G
tr G
G
G
G
=
-
+
-
=
-
-
+
2
2
2
2
et avec :
xx
xx
xx
trans
=
-
et
zz
zz
zz
trans
=
-
soit, d'après ce qui précède, on a pour
s
II
:
(
)
(
)
(
)
[
]
(
)
(
)
[
]
s
G
G
s
G
G
II
zz
trans
xx
trans
zz
zz
trans
d
zz
d
xx
d
II
trans
zz
zz
trans
d
zz
d
xx
d
= -
-
+
-
-
-
=
-
-
-
-
2
3
3
2
2
3 3
2
µ
µ
µ
µ
et pour
I
1
:
(
)
(
)
I
I
G
G
trans
d
xx
d
zz
d
1
1
3
2
2
=
-
+
+
µ
On en déduit que la fonction de charge déviatoire s'écrit :
(
)
(
)
(
)
[
]
(
)
(
)
(
)
f
s
G
G
R I
R
G
G
d
II
trans
zz
zz
trans
d
zz
d
xx
d
m
trans
m
d
xx
d
zz
d
=
-
-
-
-
-
-
+
-
+
+
1
2
3 3
2
1
3
2
2
1 6
1 6
1
µ
µ
µ
/
/
En tenant compte du fait que
( )
f
d
trans
=
0
, on trouve alors pour le multiplicateur plastique :
(
)
(
)
(
)
(
) (
)
µ
µ
µ
d
zz
d
xx
d
m
xx
d
zz
d
zz
zz
trans
G
G
R
G
G
=
-
-
-
+
+
-
3
1
2
3
2
3
2
2
1 6
/
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
7/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
ce qui donne avec les formules de
G
xx
d
et
G
zz
d
précédentes :
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
µ
µ
µ
d
m
m
zz
zz
trans
R
R
=
-
+
- -
-
-
+
-
2
3
1
3
1
2
1
3
2
1 6
2
1 6
1 6
/
/
/
On en conclut finalement l'expression analytique des contraintes :
En posant :
( )
(
)
a
b
= -
=
+
1
3
1 6
2
/
;
On a :
(
)
(
)
(
)
(
)
(
)
µ
µ
µ
µ
µ
xx
xx
trans
m
m
m
m
m
zz
zz
trans
a b
a R
a
R
b
R
a
b
R
a
a
R
-
=
-
+
+
-
+
+








+
-




-
-
+








-
2
3
2
1
6
3
1
6
1
3
2
3
2
(
)
(
)
(
)
(
)
(
)
µ
µ
µ
µ
µ
zz
zz
trans
m
m
m
m
m
zz
zz
trans
a b
a R
a
R
b
R
a
b
R
a
a
R
-
=
-
-
+
-
+
-
+










+
-




-
-
+










-
2
2
3
2
2
3
3
2
3
1
3
2
3
2

2.3
Résultats de référence
Contraintes
xx
,
yy
et
zz
aux points A, B et C.

2.4
Incertitude sur la solution
Solution analytique exacte pour CJS1.
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
8/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation
3D :
z
x
y
B
Découpage : 1 en hauteur, en largeur et en épaisseur.
Chargement de la phase 1 :
Pression de confinement :
xx
yy
zz
0
0
0
=
=
: ­100 kPa.
Niveau 1 du modèle CJS

3.2
Caractéristique du maillage
Nombre de noeuds : 8
Nombre de mailles et types : 1
HEXA8
et 6
QUA4

3.3 Fonctionnalités
testées
Commandes
DEFI_MATERIAU CJS
STAT_NON_LINE COMP_INCR
RELATION
'CJS'
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
9/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
4
Résultats de la modélisation A
4.1 Valeurs
testées
Localisation Numéro
d'ordre
Déformation
axiale
zz
(%)
Contrainte
(kPa)
Référence
Aster %
différence
Point A et B
1
­0.25
xx
­78.461538 ­78.461538 <
10­5
2
­0.50
xx
­56.923077 ­56.923077 <
10­5
3
­0.75
xx
­53.606 ­53.606 <
10
­5
4
­1.0
xx
­54.480 ­54.480 <
10
­5
8
­5.0
xx
­68.467 ­68.467 <
10
­5
23
­20.0
xx
­120.918 ­120.918
<
10
­5
1
­0.25
yy
­78.461538 ­78.461538 <
10
­5
2
­0.50
yy
­56.923077 ­56.923077 <
10
­5
3
­0.75
yy
­53.606 ­53.606 <
10
­5
4
­1.0
yy
­54.480 ­54.480 <
10
­5
8
­5.0
yy
­68.467 ­68.467 <
10
­5
23
­20.0
yy
­120.918 ­120.918
<
10
­5
1
­0.25
zz
­143,07692 ­143,07692
<
10
­5
2
­0.50
zz
­186.153846 ­186.153846
<
10
­5
3
­0.75
zz
­196.818 ­196.818
<
10
­5
4
­1.0
zz
­200.028 ­200.028
<
10
­5
8
­5.0
zz
­251.383 ­251.383
<
10
­5
23
­20.0
zz
­443.961 ­443.961
<
10
­5


background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
10/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
5 Modélisation
B
5.1
Caractéristiques de la modélisation
Cette modélisation diffère de la précédente par la finesse du maillage
3D :
z
x
y
B
Découpage : 2 en hauteur, en largeur et en épaisseur.
Chargement de la phase 1 :
Pression de confinement :
xx
yy
zz
0
0
0
=
=
: ­100 kPa.
Niveau 1 du modèle CJS

5.2
Caractéristique du maillage
Nombre de noeuds : 27
Nombre de mailles et types : 8
HEXA8
et 24
QUA4

5.3 Fonctionnalités
testées
Commandes
DEFI_MATERIAU CJS
STAT_NON_LINE COMP_INCR
RELATION
'CJS'
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
11/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
6
Résultats de la modélisation B
6.1 Valeurs
testées
Localisation Numéro
d'ordre
Déformation
axiale
zz
(%)
Contrainte
(kPa)
Référence
Aster %
différence
Point A, B et C 5
­0.2
xx
­82.76923 ­82.76923 <
10
­5
10
­0.4
xx
­65.53846 ­65.53846 <
10
­5
20
­0.8
xx
­53.78079 ­53.78079 <
10
­5
40
­1.6
xx
­56.578176 ­56.578176 <
10
­5
60
­5.6
xx
­70.565109 ­70.565109 <
10
­5
100
­20.0
xx
­120.918065 ­120.918065
<
10
­5
5
­0.2
yy
­82.76923 ­82.76923 <
10
­5
10
­0.4
yy
­65.53846 ­65.53846 <
10
­5
20
­0.8
yy
­53.78079 ­53.78079 <
10
­5
40
­1.6
yy
­56.578176 ­56.578176 <
10
­5
60
­5.6
yy
­70.565109 ­70.565109 <
10
­5
100
­20.0
yy
­120.918065 ­120.918065
<
10
­5
5
­0.2
zz
­134.46154 ­134.46154
<
10
­5
10
­0.4
zz
­168.92308 ­168.92308
<
10
­5
20
­0.8
zz
­197.460849 ­197.460849
<
10
­5
40
­1.6
zz
­207.731697 ­207.731697
<
10
­5
60
­5.6
zz
­259.085935 ­259.085935
<
10
­5
100
­20.0
zz
­443.961194 ­443.961194
<
10
­5


background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
12/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
7 Modélisation
C
7.1
Caractéristiques de la modélisation
3D_HM :
z
x
y
B
Découpage : 1 en hauteur, en largeur et en épaisseur.
Chargement de la phase 1 :
Pression de confinement :
xx
yy
zz
0
0
0
=
=
: ­100 kPa.
Niveau 1 du modèle CJS
Coefficient de biot : 1
UN_SUR_K de l'eau : 0

7.2
Caractéristique du maillage
Nombre de noeuds : 20
Nombre de mailles et types : 1
HEXA20
et 6
QUA8

7.3 Fonctionnalités
testées
Commandes
DEFI_MATERIAU CJS
STAT_NON_LINE COMP_INCR RELATION 'KIT_HM'
`RELATION_KIT' :
`CJS' `LIQU_SATU' `HYDR_UTIL'
background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
13/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
8
Résultats de la modélisation C
8.1 Valeurs
testées
Localisation Numéro
d'ordre
Déformation
axiale
zz
(%)
Contrainte
(kPa)
Référence
Aster %
différence
Point A et B
1
­0.25
xx
­78.461538 ­7,8462E+01
<
10
­5
2
­0.50
xx
­56.923077 ­5,6923E+01
<
10
­5
3
­0.75
xx
­53.606 ­5,3606E+01 <
10
­5
4
­1.0
xx
­54.480 ­5,4480E+01 <
10
­5
8
­5.0
xx
­68.467 ­6,8467E+01 <
10
­5
23
­20.0
xx
­120.918 ­1,2092E+02 <
10
­5
1
­0.25
yy
­78.461538 ­7,8462E+01
<
10
­5
2
­0.50
yy
­56.923077 ­5,6923E+01
<
10
­5
3
­0.75
yy
­53.606 ­5,3606E+01 <
10
­5
4
­1.0
yy
­54.480 ­5,4480E+01 <
10
­5
8
­5.0
yy
­68.467 ­6,8467E+01 <
10
­5
23
­20.0
yy
­120.918 ­1,2092E+02 <
10
­5
1
­0.25
zz
­143,07692 ­1,4308E+02
<
10
­5
2
­0.50
zz
­186.153846 ­1,8615E+02
<
10
­5
3
­0.75
zz
­196.818 ­1,9682E+02 <
10
­5
4
­1.0
zz
­200.028 ­2,0003E+02 <
10
­5
8
­5.0
zz
­251.383 ­2,5138E+02 <
10
­5
23
­20.0
zz
­443.961 ­4,4396E+02 <
10
­5
1
­0.25
pression eau 2,1538E+04
2.15385E+04
< 10
­5
2
­0.50
pression eau 4,3077E+04
4.30769E+04
< 10
­5

Pour la pression d'eau, on a la référence tant que le comportement est élastique linéaire

background image
Code_Aster
®
Version
7.4
Titre :
WTNV100 - Essai triaxial non drainé avec le modèle CJS (niveau 1)
Date :
05/07/05
Auteur(s) :
C. CHAVANT, Ph. AUBERT
Clé
:
V7.31.100-B
Page :
14/14
Manuel de Validation
Fascicule V7.31 : Thermo hydro mécanique en milieu poreux de structures volumiques
HT-66/05/005/A
9
Synthèse des résultats
Les valeurs du Code_Aster sont en parfait accord avec les valeurs de référence. Concernant le
couplage avec l'hydraulique, ce test prouve qu'informatiquement, le couplage CJS/THM fonctionne et
que les équations d'hydraulique sont au moins capables de redonner la variation de volume nulle
quand l'eau est incompressible.