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 :

  1. Les cotés nommés FIXE1 (ou FIXE2) sont encastrés (\(\mathit{DX}=0\) , \(\mathit{DY}=0\) ).

  2. Une condition de contact est définie entre le lèvres BFISH1-BFISB1 (BFISH2-BFISB2) (DEFI_CONTACT).

  3. 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.

  4. 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.

_images/10000000000003810000025E03E35BE38FC02D12.png

Figure 5.2-a : Maillage 1 pour la direction \(theta=120°\) (par rapport à la verticale).

_images/10000000000003AA00000260983095712DC68E77.png

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.