3. Modélisation A#
3.1. Caractéristiques de la modélisation#
Le but de cette modélisation est d’obtenir la solution au problème thermique à partir du calcul réalisé sur 1/8ème de la structure.
Le cas test se déroule la façon suivante :
lecture du maillage, 1/8ème de structure, à partir du fichier au format “MED”, commande LIRE_MAILLAGE.
lecture du champ de température à partir du fichier au format “MED”, commande LIRE_CHAMP.
création d’un résultat à partir du champ précédemment lu, commande CREA_RESU.
création d’un modèle à partir du maillage précédemment lu, commande AFFE_MODELE.
création d’un modèle pour les groupes de mailles qui appartiennent au plan de symétrie, commande AFFE_MODELE (Cf remarque n°1).
lecture du maillage, 1/8ème de structure, à partir du fichier au format “MED” et modification du maillage par symétrie par rapport au plan (\({\pi}_{\mathrm{sym}}\) ) défini par :
AXE_1= \((\mathrm{1.0,}1.0,0.0)\) , AXE_2= \((0.0,0.0,-1.0)\) , POINT= \((0.0,0.0,0.0)\)
commandes LIRE_MAILLAGE et MODI_MAILLAGE.
lecture du champ de température à partir du fichier au format “MED”, commande LIRE_CHAMP.
création d’un résultat à partir du champ précédemment lu, commande CREA_RESU.
création d’un modèle à partir du maillage symétrique, commande AFFE_MODELE.
lecture du maillage représentant 1/4 de la structure, commande LIRE_MAILLAGE.
création d’un modèle à partir du maillage précédemment lu, commande AFFE_MODELE.
projection des 3 résultats créés à partir de la solution calculée sur 1/8ème de structure, commande PROJ_CHAMP avec le mot clé DISTANCE_MAX.
extraction des champs de température des 3 résultats issus de la projection, commande CREA_CHAMP.
création d’un champ nul sur le modèle construit sur 1/4 de la structure, commande CREA_CHAMP (Cf remarque n°2).
combinaison de tous les champs, commande CREA_CHAMP/OPERATION=’ASSE’.
création d’un résultat à partir de la combinaison des champs, CREA_RESU.
Remarque n°1 :
Dans ce cas test, des nœuds du maillage, [Figure ], appartiennent au plan de symétrie \({\pi}_{\mathrm{sym}}\) *. Pour ces nœuds, la projection du champ de température sera donc comptabilisé 2 fois. Pour éviter de cumuler le champ de température une fois de trop, une solution est de créer un modèle contenant uniquement ces nœuds et de réaliser la projection sur le maillage complet. Ce champ projeté sera ensuite soustrait à l’aide de la commande* CREA_CHAMP/OPERATION=’ASSE’.
Remarque n°2 :
La projection des champs à l’aide de la commande PROJ_CHAMP et du mot clé DISTANCE_MAXpermet de ne pas créer de champ sur les nœuds qui ne sont pas contenus dans un des éléments du maillage initial et qui sont à une distance supérieure à DISTANCE_MAXde l’élément le plus proche. Lorsque l’on combine des champs à l’aide de la commande CREA_CHAMP/OPERATION=’ASSE’, c’est le premier champ aux nœuds qui est pris comme référence. Si le champ est incomplet, comme dans notre cas, la combinaison ne va pas donner le résultat escompté. La solution est donc de créer un champ nul sur tout le modèle et de s’en servir comme champ de référence.
3.2. Caractéristiques du maillage#
La projection du champ de température calculé sur 1/8ème de structure (maillage de la [Figure1.1-a]) est projeté sur le modèle construit à partir du maillage du 1/4 de structure [Figure3.2-a]).
|
|
Figure 3.2-a: Maillage du 1/4 de la structure.
3.3. Grandeurs testées et résultats#
Le test est réalisé sur le champ de température.
Du résultat obtenu par projection et combinaison sur 1/4 de structure est extrait le champ de température TEMP1. Du résultat calculé sur 1/4 de structure est extrait le 2nd champ de température TEMP2. Ces 2 champs sont soustrait à l’aide de la commande CREA_CHAMP/OPERATION=”ASSE”, et le test est réalisé sur la valeur maximale et minimale du champ résultant. En tous points de la structure, on doit avoir TEMP1=TEMP2, le champ résultant de la commande CREA_CHAMP/OPERATION=”ASSE” doit donc être nul en tous points.
Valeurs testées |
Référence |
Code_Aster |
Précision |
Maximum |
0.0 |
2.694E-05 |
1.0E-4 |
Minimum |
0.0 |
-1.808E-05 |
1.0E-4 |
La [Figure] donne la carte des différences entre les deux champs de température obtenus lors de ce cas test.
Figure 3.3-a: Carte des différences de température.