Remerciez-le!

Remerciez @Admin pour avoir partagé cet document gratuitement, de la manière la plus simple, en partageant sur les réseaux sociaux.

TABLE DES MATIÈRES 1 TABLE DES ÉVOLUTIONS Tutorial 1 : Post-traitements des rés

TABLE DES MATIÈRES 1 TABLE DES ÉVOLUTIONS Tutorial 1 : Post-traitements des résultats Utilisation des outils connexes à Code Aster R ⃝ Version 1 28 janvier 2007 Table des matières 1 Table des évolutions 1 2 Avant propos 2 3 Introduction 2 4 Types de résultats exploitables en post-traitement 2 5 Post-traitements sous GMSH 2 5.1 Valeurs moyennes aux éléments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 5.2 Extrapolation aux noeuds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 5.2.1 Mise en données sous EFICAS R ⃝ . . . . . . . . . . . . . . . . . . . . . . . . 3 5.2.2 Calcul thermique linéaires - stationnaire . . . . . . . . . . . . . . . . . . . . . 3 5.2.3 Calcul mécanique - statique . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 5.3 Valeurs aux Points d’intégration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 6 Post-traitements sous GIBI R ⃝ 8 6.1 Lancement/arrêt de GIBI R ⃝. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 6.2 Valeurs moyennes aux éléments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 6.3 Extrapolation aux noeuds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 6.3.1 Mise en données sous EFICAS R ⃝ . . . . . . . . . . . . . . . . . . . . . . . . 11 6.3.2 Calcul thermique linéaires - stationnaire . . . . . . . . . . . . . . . . . . . . . 11 6.3.3 Calcul mécanique - statique . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 6.3.4 Fichier de données GIBI R ⃝pour le calcul thermique . . . . . . . . . . . . . . . 14 6.3.5 Fichier de données GIBI R ⃝pour le calcul mécanique . . . . . . . . . . . . . . . 19 7 Post-traitement sous SALOME R ⃝ 21 8 Post-traitement sous GRACE R ⃝ 21 9 Conclusion, remerciements, auteur(s) 22 1 Table des évolutions Version Date Modiﬁcations Auteur(s) 01 28/01/2007 Création PC Page 1 1 28 janvier 2007 5 POST-TRAITEMENTS SOUS GMSH Mots clefs : Post-traitement graphique, GMSH, GIBI R ⃝, SALOME R ⃝, GRACE R ⃝. 2 Avant propos J’ai pris le parti de séparer les calculs à proprement parlé des considérations de pré & post-traitement. En effet, nombres d’astuces pourront être partagées par la communauté et il m’apparaît plus simple (et plus pertinent) de modiﬁer un seul et unique document plutôt que de mettre à jours x tutoriaux pour la même idée ... et cela évite de rendre les documents indigestes par une abondance de détails sans rapport forcement avec le thème du document. Enﬁn, les exemples traités ici reprennent les résultats des tutoriaux présents sur le site de CAELI- NUX, ces derniers faisant eux-même référence à cette note pour la partie post-traitement. 3 Introduction Code aster R ⃝est un solveur c’est-à-dire qu’il ne traite que de la partie calcul. Concrètement cela signiﬁe que le maillage et la visualisation graphique des résultats requièrent l’utilisation d’outils spé- ciﬁques (libres ou propriétaires). Fort heureusement Code Aster R ⃝est livré avec certains d’entre eux que nous nous proposons d’utili- ser ici : 1. GMSH et GIBI R ⃝pour la partie visualisation des isovaleurs, 2. GRACE R ⃝pour la réalisation de courbes, Cette liste n’est pas exhaustive ; il est possible d’utiliser d’autres outils comme PARAVIEW R ⃝puisque Code Aster R ⃝exporte les résultats au format ENSIGHT R ⃝, ou CGX. Nous nous intéresserons également à SALOME R ⃝(pré & post-traitement au format MED), disponible sous licence GPL auquel EDF c ⃝participe activement. 4 Types de résultats exploitables en post-traitement Le post-traitement graphique des résultats peut se faire : ▷Par élément (valeur moyenne par élément), ▷Par extrapolation aux noeuds, ▷Directement aux points d’intégration, 5 Post-traitements sous GMSH GMSH est un outil assez complet et probablement le plus simple à utiliser. Il permet en effet de vi- sualiser champs scalaires (contraintes équivalentes, isovaleurs de déplacements, etc ...) et champs vectoriels ; ils disposent de plugins intéressants (coupes, visualisations de groupes de mailles, etc. ...). Il est surtout d’une prise en main facile où tout se fait à la souris (translations, rotations, grossisse- ments, réductions, etc. ....). Il présente néanmoins l’inconvénient de ne pas prendre en compte les éléments du second ordre (SEG3, QUAD8, TETRA10 et HEXA20) ; les développeurs de Code Aster ont contourné le problème en créant un nouveau maillage tel que les noeuds deviennent tous des noeuds sommet. Le principal Page 2 2 28 janvier 2007 5.1 Valeurs moyennes aux éléments 5 POST-TRAITEMENTS SOUS GMSH inconvénient se situe au niveau de la taille des ﬁchiers de post-traitement qui peuvent "‘exploser"’ pour de gros maillages quadratiques (particulièrement en transitoire ou en non-linéaires si plusieurs pas de temps sont archivés). NOTA : GMSH est assez intuitif et on ne saurait trop conseiller le lecteur de consacrer un peu de son temps au test des différentes fonctionnalités ; nous traiterons donc uniquement de la façon d’exporter les résultats depuis Code Aster. 5.1 Valeurs moyennes aux éléments Sauf erreur, GMSH ne permet pas ce type d’afﬁchage et il faudra lui préférer un autre outil comme GIBI R ⃝par exemple. 5.2 Extrapolation aux noeuds 5.2.1 Mise en données sous EFICAS R ⃝ Nous pensons qu’il est plus pertinent de lancer les simulations en 2 temps : ▷Le calcul à proprement parlé, ▷Le post-traitement, Il est en effet dommage de relancer tout un calcul (particulièrement si celui-ci est long) pour une problématique de post-traitement (recherche d’un instant particulier, calcul d’une grandeur supplé- mentaire ... ou tout simplement pour une erreur de mise en données). Sous ASTK R ⃝, le ﬁchier est par exemple de la forme : F comm /symetrie_cyclique_post.comm D 1 R base /base_symetrie_cyclique D 0 F mess /symetrie_cyclique_post.mess R 6 F erre /symetrie_cyclique_post.erre R 9 F resu /symetrie_cyclique_post.resu R 8 F pos /symetrie_cyclique_vmis.pos R 51 F pos /symetrie_cyclique_depl.pos R 52 RAPPEL : ▷Les unités 1,6,9,8 sont spéciﬁques et correspondent respectivement au .comm, .mess, .erre et .resu ... ils sont attribués par ASTK R ⃝et le lecteur n’a pas à les modiﬁer : 51 & 52 correspondent aux unités des ﬁchiers GMSH qui ont été affectées dans le .comm., ▷Les lettres R et D correspondantes signiﬁent R(esult) & D(ata), ▷Tout à droite les lettres F & R désignent F(ile) & R(epertory). Par la suite, nous distingueront les calculs linéaires (mécanique & thermiques) des calculs non- linéaires & transitoires. 5.2.2 Calcul thermique linéaires - stationnaire L ’idée est ici de visualiser les champs de températures dans la pièce (champ scalaire). Le ﬁchier de post-traitement sous Code Aster R ⃝est de la forme : POURSUITE(); RESOL_TH=CALC_ELEM(reuse =RESOL_TH, RESULTAT=RESOL_TH, OPTION=’FLUX_ELNO_TEMP’,); RESOL_TH=CALC_NO(reuse =RESOL_TH, Page 3 3 28 janvier 2007 5.2 Extrapolation aux noeuds 5 POST-TRAITEMENTS SOUS GMSH RESULTAT=RESOL_TH, OPTION=’FLUX_NOEU_TEMP’,); IMPR_RESU(FORMAT=’GMSH’, UNITE=51, RESU=_F(RESULTAT=RESOL_TH, NOM_CHAM=’TEMP’,),); FIN(); CALC_ELEM et CALC_NO permettent de calculer les températures aux noeuds par extrapolation. L ’information est contenue dans FLUX_NOEU_TEMP (ﬂux thermique). Les ﬁgures 1 et 2 présentent le résultat graphique. 5.2.3 Calcul mécanique - statique En premier lieu, il est intéressante de mettre visuellement en évidence la transformation d’un maillage quadratique sous gmsh ; la ﬁgure 3 montre le maillage avant calcul alors que la ﬁgure 4 le montre après calcul1. En statique linéaire, les grandeurs physiques recherchées sont généralement les isovaleurs de dé- placements et les contraintes équivalentes (Von Mises) ... des grandeurs scalaires dans les 2 cas. Le ﬁchier de données de Code Aster est de la forme : POURSUITE(); RESOL_ME=CALC_ELEM(reuse =RESOL_ME, RESULTAT=RESOL_ME, OPTION=(’SIGM_ELNO_DEPL’,’EQUI_ELNO_SIGM’,),); RESOL_ME=CALC_NO(reuse =RESOL_ME, RESULTAT=RESOL_ME, OPTION=(’SIGM_NOEU_DEPL’,’EQUI_NOEU_SIGM’,),); IMPR_RESU(FORMAT=’GMSH’, UNITE=51, RESU=_F(RESULTAT=RESOL_ME, NOM_CHAM=’EQUI_NOEU_SIGM’, NOM_CMP=’VMIS’,),); IMPR_RESU(FORMAT=’GMSH’, UNITE=52, RESU=_F(RESULTAT=RESOL_ME, NOM_CHAM=’DEPL’,),); FIN(); A l’instar du cas précédent, les valeurs de Von Mises extrapolées sont contenues dans EQUI_NOEU_DEPL et VMIS. Les résultats sont générés dans 2 ﬁchiers distincts (un pour les contraintes uploads/Sante/ caelinux-post-traitement-pdf.pdf