3. Modélisation B#
3.1. Géométrie#
3.2. Propriétés du matériau#
Le matériau est celui défini pour le cas-test de THM wtnl100a
3.3. Conditions aux limites et chargements#
Le calcul est en mécanique non linéaire hydro-mécanique saturée. Après chaque adaptation, le calcul est initialisé par les résultats du calcul précédent, interpolés sur le nouveau maillage. On regardera l’évolution du déplacement sur un nœud de la surface supérieure.
Face |
Mécanique |
Hydraulique |
Supérieure |
Contrainte imposée |
Flux nul |
Inférieure |
Déplacement nul |
Flux nul |
Latérale |
Contrainte nulle |
Pression imposée |
Problème mécanique :
La pièce est bloquée sur la face inférieure :
Face Z_MIN : DX = DY = DZ = 0
On applique une pression sur la face supérieure :
Face Z_MAX: PRES = 1.0 105
Les autres bords sont à contrainte nulle.
Problème hydraulique :
On applique une pression sur les faces latérales :
Faces COTE_0 , COTE_1 , COTE_2 , COTE_3 : PRE1 = 1.0 105
Les autres bords sont à flux nul.
3.4. Caractéristiques du maillage#
Le maillage initial avant raffinement.
Nœuds : 622
TRIA6 : 148
TET10 : 339
3.5. Résultats de référence#
Déplacement DZ pour le groupe de nœud A, constitué d’un seul nœud, après la 3ème adaptation :
DZ = -6,23819060503x10-2
3.6. Maillages adaptés#
La boucle python de raffinement de maillage comporte 3 itérations à partir du saut du champ de déplacement mécanique d’un nœud à son voisin. Pour chaque itération, on décrit les caractéristiques de chaque maillage produit par la macro-commande MACR_ADAP_MAIL.
Nœuds : 1611
TRIA6 : 362
TET10 : 901
3.7. Remarques#
On peut constater que les nœuds issus des découpages des segments sur la frontière vont se placer sur la description analytique de la frontière.
On regardera avec attention le mécanisme employé pour relire les champs aux points de Gauss.