5. Modélisation C#
5.1. Caractéristiques de la modélisation C#
Modélisation plane: D_PLAN_HHM
2 mailles DPTR6 de la modélisation D_PLAN_HHM : HHM_ DPTR6
5.2. Résultat de la modélisation C#
Discrétisation en temps: Plusieurs pas de temps (16) pour étudier l’évolution de la pression pendant la phase transitoire jusqu’à se stabiliser. Le schéma en temps est implicite \((\theta =1)\) .
Liste des instants de calcul en secondes:
\(1,5,10,50,100,500,{10}^{3},5.\times {10}^{3},{10}^{4},5.\times {10}^{4},{10}^{5},5.\times {10}^{5},{10}^{6},5.\times {10}^{6},{10}^{7},{10}^{10}.\)
Les inconnues nodales : pressions de fluides évaluées dans Code_Aster sont des variations par rapport aux pressions initiales, c’est pourquoi ce tableau présente des variations de pression dans notre comparaison entre le calcul Code_Aster et la solution de référence. De plus les variables de pression utilisées dans Code_Aster pour évaluer les lois de comportement sont la pression totale du gaz et la pression capillaire.
Nœud / point |
Numéro d’ordre / instant \((s)\) |
Valeur |
Pression \((\mathrm{Pa})\) |
Tolérance |
1,2 / \(A\) et \(B\) |
\(1(t=1s)\) |
\(\mathrm{PRE1}\) |
-8,565.10-3 |
\({10}^{-4}\) |
\(2(t=5s)\) |
\(\mathrm{PRE1}\) |
-4,282.10-2 |
\({10}^{-4}\) |
|
\(3(t=10s)\) |
\(\mathrm{PRE1}\) |
-8,565.10-2 |
\({10}^{-4}\) |
|
\(4(t=50s)\) |
\(\mathrm{PRE1}\) |
-4,282.10-1 |
\(1\text{\%}\) |
|
\(8(t={5.10}^{3}s)\) |
\(\mathrm{PRE1}\) |
-4,26.10+1 |
\(1\text{\%}\) |
|
\(16(t={10}^{10}s)\) |
\(\mathrm{PRE1}\) |
-4,996.10+3 |
\(1\text{\%}\) |
|
\(1(t=1s)\) |
\(\mathrm{PRE2}\) |
6,796.10-6 |
\({10}^{-4}\) |
|
\(2(t=5s)\) |
\(\mathrm{PRE2}\) |
3,398.10-5 |
\({10}^{-4}\) |
|
\(3(t=10s)\) |
\(\mathrm{PRE2}\) |
6,796.10-5 |
\({10}^{-4}\) |
|
\(4(t=50s)\) |
\(\mathrm{PRE2}\) |
3,398.10-4 |
\({10}^{-4}\) |
|
\(8(t={5.10}^{3}s)\) |
\(\mathrm{PRE2}\) |
3,384.10-2 |
\({10}^{-4}\) |
|
\(16(t={10}^{10}s)\) |
\(\mathrm{PRE2}\) |
3,964 |
\({10}^{-3}\) |
|
3,4 / \(C\) et \(D\) |
\(1(t=1s)\) |
\(\mathrm{PRE1}\) |
8,565.10-3 |
\({10}^{-4}\) |
\(2(t=5s)\) |
\(\mathrm{PRE1}\) |
4,282.10-2 |
\({10}^{-4}\) |
|
\(3(t=10s)\) |
\(\mathrm{PRE1}\) |
8,565.10-2 |
\({10}^{-4}\) |
|
\(4(t=50s)\) |
\(\mathrm{PRE1}\) |
4,282.10-1 |
\(1\text{\%}\) |
|
\(8(t={5.10}^{3}s)\) |
\(\mathrm{PRE1}\) |
4,26.10+1 |
\(1\text{\%}\) |
|
\(16(t={10}^{10}s)\) |
\(\mathrm{PRE1}\) |
4,996.10+3 |
\(1\text{\%}\) |
|
\(1(t=1s)\) |
\(\mathrm{PRE2}\) |
-6,796.10-6 |
\({10}^{-4}\) |
|
\(2(t=5s)\) |
\(\mathrm{PRE2}\) |
-3,398.10-5 |
\({10}^{-4}\) |
|
\(3(t=10s)\) |
\(\mathrm{PRE2}\) |
-6,796.10-5 |
\({10}^{-4}\) |
|
\(4(t=50s)\) |
\(\mathrm{PRE2}\) |
-3,398.10-4 |
\({10}^{-4}\) |
|
\(8(t={5.10}^{3}s)\) |
\(\mathrm{PRE2}\) |
-3,384.10-2 |
\({10}^{-4}\) |
|
\(16(t={10}^{10}s)\) |
\(\mathrm{PRE2}\) |
-3,964 |
\({10}^{-3}\) |