4. Modélisation B#

4.1. Caractéristiques de la modélisation#

Dans cette modélisation, on valide l’utilisation de la macro-commande MACR_ECREVISSE en présence de 2 fissures verticales. La géométrie est légèrement différente de celle de la modélisation A. Elle est décrite avec le maillage dans la . Les caractéristiques de l’écoulement et du matériau ont étés décrites au § 1.

Comme pour la modélisation A, on travaille en déformations planes pour la mécanique (D_PLAN), et avec une matrice de masse diagonalisée pour la thermique (PLAN_DIAG).

4.2. Maillage et conditions aux limites#

Le maillage utilisé est présenté :

Figure 4.2-a : Maillage complet

Les dimensions des chacun des trois blocs sont de \(0.50m\) de haut et \(0.25m\) de large.

Nombre de nœuds : 2275

Nombre de mailles : 2108 QUAD4

Les paramètres de croissance des tailles de mailles utilisés pour réaliser le maillage sont les suivants : de \(0.002m\) (intrados) à \(0.05m\) (extrados), et idem de BFISG1/BFISD2 à ENCASTRE ; de \(0.025m\) au MILIEU à \(0.002m\) sur BFISM1/BFISM2.

Les fissures sont dénommées avec les nom des lèvres : BFISG1-BFISM1 (FISSURE 1), (FISSURE 2).

Les conditions aux limites :

Les cotés nommés ENCASTRE sont encastrés (\(\mathit{DX}=0\) , \(\mathit{DY}=0\) ).
Sur le coté MILIEU (ligne de symétrie) on bloc le déplacement en direction \(X\) (\(\mathit{DX}=0\) ).
Une condition de contact est définie entre le lèvres BFISG1-BFISM1 et BFISM2-BFISD2 (DEFI_CONTACT).
La pression et la température sur l’intrados INTRADOS sont plus grandes que sur l’extrados EXTRADOS (\(1.E6\mathrm{Pa}\) et \(140°C\) versus \(1.E6\mathrm{Pa}\) et \(20°C\) ), ce qui cause l’écoulement de l’intrados à l’extrados.
Le matériau est considéré, au début du calcul, à la température ambiance \(20°C\) , qui est la température de référence.

4.3. Grandeurs et résultats#

On teste à 2 instants différents, les déplacements et la température en entrée de la fissure de gauche et en sortie de la fissure de droite. Enfin, on teste le débit de la fissure de droite à \(t=500s\) et de la fissure de gauche à \(t=10000s\) .

A \(t=500s\) :

Identification	Tolérance
Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{MD}}\) \([m]\)	1.0E-06
Température du matériau au nœud \({E}_{\mathit{MD}}\) \([°C]\)	1.0E-06
Déplacement \(\mathit{DX}\) nœud \({I}_{\mathit{GD}}\) \([m]\)	1.0E-06
Température du matériau au nœud \({I}_{\mathit{GD}}\) \([°C]\)	1.0E-06
Fissure BFISG1-BFISM1, débit	1.0E-06

A \(t=10000s\) :

Identification	Tolérance
Déplacement \(\mathit{DX}\) au nœud \({E}_{\mathit{MD}}\) \([m]\)	1.0E-06
Température du matériau au nœud \({E}_{\mathit{MD}}\) \([°C]\)	1.0E-06
Déplacement \(\mathit{DX}\) au nœud \({I}_{\mathit{GD}}\) \([m]\)	1.0E-06
Température du matériau au nœud \({I}_{\mathit{GD}}\) \([°C]\)	1.0E-06
Fissure BFISM2-BFISD2, débit	1.0E-06

4.4. Remarques#

L’ouverture de la fissure à \(t=10000s\) est de l’ordre de \({10}^{-9}\) , donc est presque refermé à cause de la dilatation thermique due à la haute température à l’intrados. Dans ce cas aussi, c’est l’ouverture rémanente qui est retenue pour le calcul de l’écoulement.