v1.01.340 ZZZZ340 – Validation du mot-clé TEMP_CONTINUE pour AFFE_CHAR_THER#
Résumé
Ce test valide le mot clé TEMP_CONTINUE du mot-clé facteur ECHANGE_PAROI des opérateurs AFFE_CHAR_THER et AFFE_CHAR_THER_F [U4.44.02]. Ce mot clé ne concerne que les modèles X-FEM, et conduit à annuler tous les degrés de liberté enrichis («Heaviside» et «crack-tip») correspondant à la présence des fissures données par l’utilisateur sous le mot-clé FISSURE. La mise à zéro de ces degrés de liberté revient à résoudre l’équation de la chaleur sur un domaine ne comportant pas ces fissures.
On valide cette fonctionnalité sur un problème simple de thermique linéaire stationnaire en dimension 2. On considère une plaque carrée à température imposée, comportant une fissure droite débouchante disposée dans la direction orthogonale au gradient de température. L’utilisation du mot-clé TEMP_CONTINUE permet de ne pas prendre en compte la discontinuité du champ de température à travers cette fissure et de se ramener à la résolution du problème «sain», ce dernier admettant une solution analytique.
On valide également dans ce test le calcul du gradient de la température aux points de Gauss des éléments X-FEM dans l’opérateur THER_LINEAIRE (champ TEMP_ELGA).
On ne considère qu’une seule modélisation:
modélisation \(A\) : \(\text{X-FEM}\mathrm{2D}\) (fissure au milieu des éléments)
Solution de référence#
Méthode de calcul#
Il s’agit d’une solution analytique. Dans cette configuration le problème est unidimensionnel et l’équation de la chaleur se réduit à l’équation différentielle suivante:
\(\frac{{d}^{2}T}{{\mathit{dy}}^{2}}=0\) avec les conditions \(T(y=0)=\text{}\stackrel{ˉ}{T}{\rbrace ^{\inf}\) et \(T(y=\mathit{LY})=\text{}\stackrel{ˉ}{T}{\rbrace ^{\sup}\)
Cette équation admet la solution suivante:
\(T(y)=(\text{}\stackrel{ˉ}{T}{\rbrace ^{\sup}-\text{}\stackrel{ˉ}{T}{\rbrace ^{\inf})\frac{y}{\mathit{LY}}+\text{}\stackrel{ˉ}{T}{\rbrace ^{\inf}\)
Grandeurs et résultats de référence#
On teste la température aux points \({P}^{+}\) , \({P}^{-}\) et \(Q\) (voir Figure ) :
\(T(\frac{\mathit{LY}}{2})=\frac{(\text{}\stackrel{̄}{T}{\rbrace ^{\sup}+\text{}\stackrel{̄}{T}{\rbrace ^{\inf})}{2}\) , application numérique: \(T(\frac{\mathit{LY}}{2})=15°C\)
On teste également les deux composantes (constantes en espace) du gradient de température:
\({\partial}_{x}T=0\) , application numérique: \({\partial}_{x}T=0°{\mathit{C.m}}^{-1}\)
\({\partial}_{y}T=\frac{(\text{}\stackrel{̄}{T}{\rbrace ^{\sup}-\text{}\stackrel{̄}{T}{\rbrace ^{\inf})}{\mathit{LY}}\) , application numérique: \({\partial}_{y}T=10°{\mathit{C.m}}^{-1}\)
Identification |
Type de référence |
Valeur de référence |
Points \({P}^{+}\) , \({P}^{-}\) et \(Q\) - \(\mathit{TEMP}\) |
“ANALYTIQUE” |
\(15°C\) |
En tout point, \(\mathit{DTX}\) |
“ANALYTIQUE” |
\(0°{\mathit{C.m}}^{-1}\) |
En tout point, \(\mathit{DTY}\) |
“ANALYTIQUE” |
\(10°{\mathit{C.m}}^{-1}\) |
Modélisation A : fissure non-maillée en dimension 2#
La fissure n’étant pas maillée, la condition d’échange entre les lèvres de la fissure est appliquée à l’aide du mot-clé FISSURE du mot-clé facteur ECHANGE_PAROI de l’opérateur AFFE_CHAR_THER [U4.44.02].
Caractéristiques de la modélisation#
On utilise la modélisation PLAN du phénomène THERMIQUE. La méthode des éléments finis étendue (X-FEM) est utilisée.
Caractéristiques du maillage#
La structure est modélisée par un maillage régulier composé de \(101\times 101\) QUAD4, respectivement suivant les axes \(x\) , \(y\) . La fissure n’est pas maillée.
Grandeurs testées et résultats#
Les grandeurs ci-dessous sont d’abord testées dans le cas où le chargement permettant d’imposer la continuité du champ température (mot-clé TEMP_CONTINUE) résulte de l’opérateur AFFE_CHAR_THER, elles sont ensuite testées lorsqu’il résulte l’opérateur AFFE_CHAR_THER_F.
On teste la température aux points \({P}^{+}\) , \({P}^{-}\) et \(Q\) (voir Figure ). Pour cela on teste le champ de température après appel aux opérateurs POST_MAIL_XFEM et POST_CHAM_XFEM.
On teste ensuite les deux composantes du gradient de la température stockées dans le champ TEMP_ELGA dans le résultat produit par l’opérateur THER_LINEAIRE. Ces composantes sont testées:
en un point de Gauss d’un élément de type «Heaviside»
en un point de Gauss d’un élément de type «crack-tip»
en un point de Gauss d’un élément de type «Heaviside/crack-tip»
Identification |
Type de référence |
Valeur de référence |
Tolérance |
Points \({P}^{+}\) , \({P}^{-}\) et \(Q\) - \(\mathit{TEMP}\) |
“ANALYTIQUE” |
15 |
0.1% |
En 3 points de Gauss \(\mathit{DTX}\) |
“ANALYTIQUE” |
0 |
1.E-9 |
En 3 points de Gauss \(\mathit{DTY}\) |
“ANALYTIQUE” |
10 |
0.1% |
Synthèses des résultats#
Les objectifs de ce test sont atteints:
valider pour les modèles X-FEM le mot clé TEMP_CONTINUE du mot-clé facteur ECHANGE_PAROI des opérateurs AFFE_CHAR_THER et AFFE_CHAR_THER_F.
valider le calcul du gradient de la température aux points de Gauss des éléments X-FEM dans l’opérateur THER_LINEAIRE