Code_Aster ®
Version
7.4
Titre :
WTNP110 - Ecoulement orthotrope saturé 2D
Date :
20/10/04
Auteur(s) :
S. GRANET Clé
:
V7.32.110 -A Page :
1/4
Organisme(s) : EDF-R&D/AMA
Manuel de Validation
Fascicule V7.32 : Thermo-hydro-mécanique en milieu poreux non saturé
Document : V7.32.110
WTNP110 - Ecoulement orthotrope saturé 2D
Résumé :
Le test présenté ici permet de vérifier le bon fonctionnement des opérateurs utilisés pour la résolution des
équations d'un écoulement en milieu saturé orthotrope. Ce test correspond au test 1.2 du plan de qualification
du projet ALLLIANCES [bib1].
La solution de référence est une solution analytique.
Manuel de Validation
Fascicule V7.32 : Thermo-hydro-mécanique en milieu poreux non saturé
HT-66/04/005/A
Code_Aster ®
Version
7.4
Titre :
WTNP110 - Ecoulement orthotrope saturé 2D
Date :
20/10/04
Auteur(s) :
S. GRANET Clé
:
V7.32.110 -A Page :
2/4
1
Problème de référence
1.1 Géométrie
D C
(0 ;0)
A B
Coordonnées des points (m) :
A -0,1
-0,1
C 0,1 0,1
B 0,1 -0,1
D -0,1 0,1
1.2
Propriétés du matériau
On ne donne ici que les propriétés dont la solution dépend, sachant que le fichier de commandes
contient d'autres données de matériau (modules d'élasticité, ...) qui finalement ne jouent aucun rôle
dans la solution du problème traité.
Eau liquide
Masse volumique (kg.m-3)
1
Viscosité
1
Paramètres
Perméabilité intrinsèque (composante en X)
1
homogénéisés
Perméabilité intrinsèque (composante en Y)
3/4
Perméabilité intrinsèque (composante en Z)
1
Etat initial
Porosité
1
Pression de liquide
0
1.3
Conditions aux limites et chargements
1.3.1 Modélisation
A
Sur AB P = P(X) = -45X+30.5
Sur BC P = P(X) = -80Y+18
Sur CD P = P(X) = -45X+14.5
Sur DA P = P(X) = -80Y+27
Manuel de Validation
Fascicule V7.32 : Thermo-hydro-mécanique en milieu poreux non saturé
HT-66/04/005/A
Code_Aster ®
Version
7.4
Titre :
WTNP110 - Ecoulement orthotrope saturé 2D
Date :
20/10/04
Auteur(s) :
S. GRANET Clé
:
V7.32.110 -A Page :
3/4
2
Solution de référence
La solution analytique en pression est un polynôme de degré en x et y , la vitesse est constante et
horizontale :
P(x, y) = -45x - 80y +
5
,
22
V (x, y) = (45K 80
; K ) = (
)
60
;
45
x
y
2.1
Résultats de référence
On donne la valeur de la pression et de la vitesse en 3 points :
x -0,05 0,00 +0,05
y -0,05 0.00 0.05
P 28,8 22.5 16.3
Vx 45 45 45
Vy 60 60 60
3 Modélisation
A
3.1
Caractéristiques de la modélisation A
Modélisation en déformations planes. 20*20 éléments Q8
3.2 Fonctionnalités
testées
Commande Option
AFFE_MODELE
D_PLAN_HMD
DEFI_MATERIAU
THM_LIQU
THM_DIFFU
PERMIN_X
PERMIN_Y
PERMIN_Z
THM_INIT
ELAS
AFFE_CHAR_MECA DDL_IMPO
PRE1
DX
DY
STAT_NON_LINE COMP_INCR RELATION
KIT_HM
RELATION_KIT
ELAS
LIQU_SATU
HYDR_UTIL
Discrétisation en temps : un seul pas de temps : 1 s.
Manuel de Validation
Fascicule V7.32 : Thermo-hydro-mécanique en milieu poreux non saturé
HT-66/04/005/A
Code_Aster ®
Version
7.4
Titre :
WTNP110 - Ecoulement orthotrope saturé 2D
Date :
20/10/04
Auteur(s) :
S. GRANET Clé
:
V7.32.110 -A Page :
4/4
3.3 Valeurs
testées
Noeud Type
de
Instant (s)
Référence
Aster Différence
valeur
(analytique)
(%)
N441
PRE1 1 16,3
16,25 0,3
%
(0.05 ;0.05)
N689
PRE1 1 28,8
28,5 0.96
%
(-0.05 ;-0.05)
N341
PRE1 1 22,5
22,5 0.
%
(0 ;0)
Mailles Type
de
Instant (s)
Référence
Aster Différence
valeur
(analytique)
(%)
DROIT FH11X
1 45
45 0
%
DROIT FH11Y
1 60
60 0
%
4
Synthèse des résultats
Très bon accord avec la solution de référence.
5 Bibliographie
[1]
Projet Alliances plan de qualification, note ANDRA CNT-ASCS 02-075B
Manuel de Validation
Fascicule V7.32 : Thermo-hydro-mécanique en milieu poreux non saturé
HT-66/04/005/A
Document Outline