2. Modélisation A

2.1. Caractéristiques de la modélisation

Modélisation 3D _HH2MS sur 1 maille HEXA20.

1 seconde est simulée par 500 pas de temps.

2.2. Résultat de la modélisation A

La Figure 2.2-a montre l’évolution de la contrainte totale en fonction de la pression capillaire homogène en tout point. Le post traitement est ici fait au nœud \(N8\) de coordonnées \((0,5;0,5;0,5)\) . Dans la partie saturée (\(\mathit{Pc}\le 0\) ), la diminution de la pression capillaire correspond à une augmentation de pression d’eau et la contrainte totale croît linéairement. On constate que la pente de la courbe est continue.

Les paramètres \(A\) et \({\beta}_{m}\) ont été calculés de manière à retrouver une pression de gonflement de \(7\mathit{MPa}\) . Lorsque la saturation atteint 1 (ou la pression capillaire 0), la pression de gonflement est donnée par la formule suivante:

\(\frac{{P}_{\mathit{gf}}}{A}=\frac{\sqrt{\pi}}{2\sqrt{{\beta}_{m}}}+\frac{1}{2{\beta}_{m}}\)

On retrouve donc bien l’allure classique de la contrainte de gonflement et on vérifie que la courbe coupe bien l’axe des ordonnées (\(\mathit{PC}=0\) ) avec une valeur de \(7\mathit{Mpa}\) .

Figure 2.2-a test de gonflement

On rappelle sur la figure l’évolution de la pression capillaire en fonction du tempscorrespondant au chargement du problème:

Figure 2.2-b : pression capillaire (N3)

2.3. Grandeurs testées et résultats

Ce cas-test n’a pas de valeur de référence, on en fait donc un cas de non régression.

On effectue des tests sur deux valeurs :

\(N\)	Temps ( \(s\) )	\(\mathrm{SIXX}\) Aster	Erreur relative autorisée ( \(\text{\%}\) )
\(\mathrm{N8}\)	0,6	-4,56.104	0.1%
\(\mathrm{N8}\)	0.8163	-5,67.106	0.1%