5. Modélisation C#
5.1. Caractéristiques de la modélisation#
Dans cette modélisation, on valide l’utilisation de la macro-commande MACR_ECREVISSE par rapport a une rotation et à une translation. Il ne s’agit pas donc d’un test de non-régression.
En effet, les résultats d” Écrevisse ne dépendent que des caractéristiques de la fissure, y compris la direction. Celle-ci est influencée par la gravité. Ainsi, des fissures avec les même caractéristiques et ayant une direction symétrique par rapport à la verticale doivent donner les mêmes écoulements.
On travaille avec les mêmes paramètres (du matériau, de l’écoulement de la fissure) que dans la modélisation A (§ 1.1). Le maillage est aussi le même, sauf que il est tourné et translaté. Le détail du maillage est donné au §5.2.
5.2. Maillage et conditions aux limites.#
Les deux maillages utilisés ont été obtenus par rotation et translation du maillage de la modélisation A, donc ils en conservent les mêmes dimensions. Ils sont présentés dans la et .
Pour le premier maillage, la première abscisse curviligne de la fissure est le point d’origine des axes \((0;0)\) .
Le second maillage résulte d’une rotation de 120 degrés du premier maillage autour de l’origine et d’une translation de vecteur \((-10;20)\) .
Les conditions aux limites pour les deux maillages sont analogues à la modélisation A :
Les cotés nommés FIXE1 (ou FIXE2) sont encastrés (\(\mathit{DX}=0\) , \(\mathit{DY}=0\) ).
Une condition de contact est définie entre le lèvres BFISH1-BFISB1 (BFISH2-BFISB2) (DEFI_CONTACT).
La pression et la température sur l’intrados INTRA1 (INTRA2) sont plus grandes que sur l’extrados EXTRA1 (EXTRA2) : (\(1.E6\mathrm{Pa}\) et \(140°C\) versus \({P}_{\mathrm{atm}}\) et \(20°C\) ). Cela cause l’écoulement de l’intrados à l’extrados.
Le matériau est considéré, au début du calcul, à la température ambiante \(20°C\) , qui est la température de référence.
Figure 5.2-a : Maillage 1 pour la direction \(theta=120°\) (par rapport à la verticale).
Figure 5.2-b : Maillage 2 pour la direction \(theta=-120°\) (par rapport à la verticale).
5.3. Grandeurs testées et résultats#
On teste les résultats à l’instant \(t=10000s\) . On extrait les valeurs du résultat du maillage 2, et on les compare à celles correspondant au maillage 1.
Maillage 1 |
Maillage 2 |
Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{HD1}}\) |
par symétrie : (-1) * Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{BD2}}\) |
Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{BD1}}\) |
par symétrie : (-1) * Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{HD2}}\) |
Déplacement \(\mathit{DY}\) au nœud \({E}_{\mathit{HD1}}\) |
par symétrie : Déplacement \(\mathit{DY}\) au nœud \({E}_{\mathit{BD2}}\) |
Déplacement \(\mathit{DY}\) au nœud \({E}_{\mathit{BD1}}\) |
par symétrie : Déplacement \(\mathit{DY}\) au nœud \({E}_{\mathit{HD2}}\) |
Température du matériau au nœud \({I}_{\mathit{HD1}}\) |
Température du matériau au nœud \({I}_{\mathit{HD2}}\) |
Débit total |
Débit total |
Coefficient Convection première abs. curv. |
Coefficient Convection première abs. curv. |
Pression première abs. curv. |
Pression première abs. curv. |
Température du fluide première abs. curv. |
Température du fluide première abs. curv. |
Flux de chaleur du fluide première abs. curv. |
Flux de chaleur du fluide première abs. curv. |
Grandeur |
Tolérance |
Déplacement \(\mathit{DX}\) aux nœud \({E}_{\mathit{HD1}}\) |
1.0E-05 |
Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{BD1}}\) |
1.0E-05 |
Déplacement \(\mathit{DY}\) au nœud \({E}_{\mathit{HD1}}\) |
1.0E-05 |
Déplacement \(\mathit{DY}\) au nœud \({E}_{\mathit{BD1}}\) |
1.0E-05 |
Température du matériau au nœud \({I}_{\mathit{HD1}}\) |
1.0E-04 |
Débit total |
1.0E-03 |
Coefficient Convection première abs. curv. |
1.0E-04 |
Pression première abs. curv. |
1.0E-04 |
Température du fluide première abs. curv. |
1.0E-04 |
Flux de chaleur du fluide première abs. curv. |
1.0E-03 |
5.4. Remarques#
Le lèvres sont en contact, la commande DEFI_CONTACT est activée. Le calcul d’écoulement est mené avec l’ouverture rémanente.