[U2.51.02] Tracé de courbes avec Code_Aster
Contents
1. Le document U1 03 02 d taille les m thodes Python disponibles sur les objets de type Table Pour l extraction des valeurs on dispose notamment d un m thode values qui retourne un dictionnaire dont les cl s sont les noms de param tres ex DY et les valeurs les listes des valeurs de la table Attention Les listes Python sont index es de O0 n 1 pour n l ments l quivalent de tabl Dy 2 estdonc tabpy DY 1 Pour le trac de courbes partir de listes de r els Python voir 4 Manipuler les valeurs en Python On donne ici quelques exemples de manipulation des valeurs obtenues pr c demment sous forme de listes ou de tableaux Numeric Numeric est un module Python optionnel c est dire non inclus la distribution de Python fournie sur www python org mais indispensable pour utiliser Code Aster on peut donc faire import Numeric sur toutes installations de Code Aster Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster ou Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 5 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 3 1 Avec des listes Python Les listes Python sont facilement manipulables en utilisant les boucles Prenons l exemple de la notice U2 51 01 pour 10 lt x lt 10 en 100 pas y 1 5 s2. NOM PARA INST VMIS UNITE unit FORMAT XMGRACE IMPR TABLE TABLE tab mfac L int r t tant bien entendu de pouvoir construire le dictionnaire comme on le souhaite Remarque Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster oc default Titre Trac de courbes avec Code_Aster Date 01 03 2013 Page 10 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 Les mots cl s facteurs ici FILTRE peuvent tre construits en utilisant F mot cl valeur mais il est plus souple de les voir comme une liste de dictionnaires Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html
3. 02 R vision 10626 2 Transformer les valeurs d une table ou d une fonction en objets Python Remarque Les traitements en Python qui sont utilis s partir de maintenant imposent d tre en PAR LOT NON dans POURSUITE ou DEBUT L objet est ici de r cup rer les valeurs d une table ou d une fonction dans un objet Python pour les manipuler ensuite Notons qu il est parfois pratique de produire une fonction partir des donn es d une table en filtrant ventuellement certaines lignes de la table c est une utilisation de RECU FONCTION U4 32 03 que nous n aborderons pas ici Sur les objets de type fonction on dispose des m thodes 1 Valeurs pour r cup rer les abscisses et les ordonn es dans 2 listes Python de r els Avec les donn es du paragraphe pr c dent gt gt gt 1x ly fDY2 Valeurs On obtient gt gt gt print 1x ML Os Ale L gt gt gt print ly 0 2 07127 0 Zy 1 38 2451 414928 0 792728 0 688227 0 343029 1 Absc et Ordo permettent de r cup rer les abscisses et les ordonn es s par ment gt gt gt 1x fDY2 Absc gt gt gt 1x fDY2 0rdo On peut acc der au contenu d une cellule d une table avec nom param tre num ro _lignel gt gt gt print tab DY 2 8 14498E 2 On peut galement transformer l objet JEVEUX table en une instance de la classe Python Table gt gt gt tabpy tab EXTR TABLE
4. X X R F IMPR_FONCTION Filtrage de la table pour en extraire Dy f t INTITULE Umoyen QUANTITE MOMENT_0 Filtrage de la table pour en extraire Fy f t INTITULE Fresultante Trac des ordonn es de Fy en fonction de Dy UNITE 29 FORMAT XMGRACE COURBE _F FONC X Dy FONC Y Fy TITRE Plaque trou e en traction LEGENDE X D placement moyen LEGENDE Y Effort r sultant Ce qui nous donne la courbe suivante Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster default Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 8 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 Plaque trou e en traction 13125 10500 7875 Effort r sultant 5250 2625 0 048199 0 096398 0 1446 0 1928 0 24099 D placement moyen 4 3 Trac d un grand nombre de courbes Dans certaines applications on est amen tracer de nombreuses courbes Supposons que lon souhaite comparer nos r sultats 50 points de mesure obtenus par ailleurs Il serait alors fastidieux de d finir 50 fichiers dans astk et d utiliser 50 unit s logiques diff rentes dans le fichier de commandes Il suffit alors d utiliser l
5. 000E 00 1 00000E 01 0 00000E 00 7 97278E 2 7 92728E 01 1 80000E 00 1 00000E 01 0 00000E 00 7 86827E 2 6 88227E 01 2 45000E 00 1 00000E 01 0 00000E 00 8 48309E 2 3 43029E 01 On a plusieurs blocs qui n ont pas forc ment le m me nombre de colonnes Supposons que l on veuille comparer le d placement DY obtenus avec les deux valeurs de PARA on utilisera par exemple fDY1 LIRE_ FONCTION TYPE FONCTION INDIC_PARA 1 1 les abscisses sont prises dans le bloc 1 colonne 1 INDIC RESU 1 3 les ordonn es sont prises dans le bloc 1 colonne 3 UNITE 38 NOM PARA INST NOM RESU DY fDY2 LIRE FONCTION TYPE FONCTION INDIC PARA 2 1 les abscisses sont prises dans le bloc 2 colonne 1 INDIC RESU 2 5 les abscisses sont prises dans le bloc 2 colonne 5 UNITE 38 NOM PARA INST NOM RESU DY trac classique de deux fonctions avec IMPR FONCTION IMPR FONCTION FORMAT XMGRACE UNITE 29 COURBE _F FONCTION fDYI LEGENDE PARA 1 23 _F FONCTION FDY2 LEGENDE PARA 1 98 TITRE DY t Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster eul Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 4 10 Responsable Mathieu COURTOIS Cl U2 51
6. 1 04 00 ou STANLEY U4 81 31 pour simplifier le post traitement on peut utiliser les possibilit s offertes par Python qui n cessitent d tre en PAR LOT NON dans POURSUITE ce qui emp che l utilisation d eficas pour ce type de post traitement alors qu il est plus simple d utiliser eficas pour construire le jeu de donn es principal 1 Extraire les donn es l aide des commandes Aster On suppose que l utilisateur dispose d un r sultat de calcul obtenu partir par exemple de la commande MODE_IT statique non lin aire ou DYNA NON LINE pour un calcul dynamique non lin aire ER SIMULT pour un calcul de modes propres STAT NON LINE pour un calcul On dispose dans tous les cas d un concept r sultat qui sera par exemple de type mode meca mode flamb dyna trans tran gene evol elas evol noli evol ther etc selon la commande utilis e et qui contient des champs de valeurs que l on souhaite repr senter sous forme de courbes 1 1 Produire une fonction ou une table Rappel Une r els Pour Une fonction est compos e de deux listes de valeurs abscisses et ordonn es les abscisses sont n cessairement monotones table est un agglom rat de valeurs non n cessairement de m me type auxquelles on acc de via un param tre nom de colonne Dans l utilisation des tables qui nous int resse ici on produira g n ralement des colonnes de nombres leur variation est quel
7. 2 DEPL 1 1 00000E 00 2 38745E 02 MOMENT _0 Umoyen SOLNL2 DEPL 1 1 00000E 00 3 70291E 04 MOMENT_1 Umoyen SOLNL2 DEPL 1 1 00000E 00 2 36709E 02 MINIMUM Umoyen SOLNL2 DEPL 1 1 00000E 00 2 40291E 02 MAXIMUM Umoyen SOLNL2 DEPL 1 1 00000E 00 2 40597E 02 MOYE_INT Umoyen SOLNL2 DEPL 1 1 00000E 00 2 36894E 02 MOYE_EXT Umoyen a SOLNL2 DEPL 2 1 20000E 00 2 86494E 02 MOMENT_0 Umoyen SOLNL2 DEPL 2 1 20000E 00 4 44349E 04 MOMENT_1 Umoyen SOLNL2 DEPL 2 1 20000E 00 2 84050E 02 MINIMUM Umoyen a SOLNL2 DEPL 2 1 20000E 00 2 88350E 02 MAXIMUM Umoyen SOLNL2 DEPL 2 1 20000E 00 2 88716E 02 MOYE_INT Umoyen SOLNL2 DEPL 2 1 20000E 00 2 84273E 02 MOYE EXT Umoyen SOLNL2 DEPL 3 1 40000E 00 3 34244E 02 MOMENT _O Umoyen SOLNL2 DEPL 3 1 40000E 00 5 18504E 04 MOMENT 1 Umoyen ai SOLNL2 DEPL 3 1 40000E 00 3 31393E 02 MINIMUM Umoyen cs SOLNL2 DEPL 3 1 40000E 00 3 36410E 02 MAXIMUM Umoyen SOLNL2 DEPL 3 1 40000E 00 3 36837E 02 MOYE_INT Umoyen SOLNL2 DEPL 3 1 40000E 00 3 31652E 02 MOYE_EXT Fresultante a SOLNL2 FORC_NODA 1 1 00000E 00 2 50000E 03 Fresultante SOLNL2 FORC_NODA 2 1 20000E 00 3 00000E 03 S Fresultante SOLNL2 FORC_NODA 3 1 40000E 00 3 50000E 03 1 Dy RECU FONCTION TABLE tab PARA X INST PARA Y DY FILTRE _F NOM PARA VALE K _F NOM PARA VALEK Fy RECU FONCTION TABLE tab PARA CINSTE DY PARA E FILT NOM PARA VALE K
8. 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster defaki Titre Trac de courbes avec Code_Aster Date 01 03 2013 Page 3 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 1 2 Traiter des blocs de donn es Un bloc de donn es est simplement un tableau de valeurs N lignes P colonnes Dans certains cas on dispose d un ou plusieurs fichiers compos s chacun d un ou plusieurs blocs de donn es s par s par des lignes de texte Pour construire des fonctions partir de tels fichiers on peut utiliser LIRE FONCTION U4 32 02 Par exemple on effectue une tude param trique chaque calcul produit un tableau qui vient enrichir un fichier de r sultat le fichier peut aussi avoir t construit par un calcul unique ou encore manuellement peu importe Un tel fichier pourrait ressembler cela AVEC PARA 1 23 INST COOR X DY 1 00000E 00 1 00000E 01 3 02717E 2 1 20000E 00 1 00000E 01 8 14498E 2 1 40000E 00 1 00000E 01 7 97278E 2 1 60000E 00 1 00000E 01 3 86827E 2 1 80000E 00 1 00000E 01 8 48309E 2 2 00000E 00 1 00000E 01 9 37561E 2 2 20000E 00 1 00000E 01 7 18293E 2 2 40000E 00 1 00000E 01 6 05322E 2 AVEC PARA 1 98 INST COOR X COOR Y DX DY 1 00000E 00 1 00000E 01 0 00000E 00 3 02717E 2 2 07127E 01 1 10000E 00 1 00000E 01 0 00000E 00 8 14498E 2 4 14928E 01 1 20
9. Code Aster Po Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 1 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 Trac de courbes avec Code Aster R sum Ce document explique comment partir d un r sultat de calcul on peut produire des tables ou des fonctions extraire les valeurs de ces tables ou fonctions les manipuler et enfin en dernier lieu utiliser les commandes d impression pour tracer des courbes Si vous disposez d une fonction ou d une table et que vous souhaitez simplement la repr senter sous forme d une courbe allez directement au paragraphe 4 et consultez les documentations des commandes IMPR FONCTION U4 33 01 et IMPR TABLE U4 91 03 Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster s Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 2 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 Remarque pr liminaire Le 1 2 3 post traitement doit tre effectu en POURSUITE et non pas la suite du calcul Plusieurs raisons cela en cas d erreur on ne perd pas les heures de calcul qui ont permis d atteindre le r sultat on peut effectuer de nombreux post traitements directement en lan ant Stanley sur la base r sultat du calcul voir astk U
10. NCTION en prenant soin de convertir les tableaux Numeric en liste par exemple cf 4 ABSCISSE lx tolist Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version default Titre Trac de courbes avec Code Aster Responsable Mathieu COURTOIS 4 Date 01 03 2013 Page 6 10 Cl U2 51 02 R vision 10626 Exemples d utilisation de IMPR FONCTION IMPR TABLE P 2 2 4 1 Fonction reconstruite partir de deux listes Python Exemple o l on recale les r sultats obtenus pour pouvoir les comparer une solution de r f rence les abscisses sont invers es il faut comparer la valeur absolue SSLS501a on utilise les m thodes permettant d extraire les abscisses et les ordonn es d une fonction IMPR FONCTION FORMAT XMGRACE UNITE 53 COURBE F ABSCISSE 57766 1 x for x in MTAST Absc ORDONNEE abs y for y in MTAST Ordo TITRE Courbe recal e GENDE X P2 GENDE Y MT 42 Trac d un r sultat en fonction d un autre Cet exemple est extrait en partie du cas test FORMAO3a1 il s agit dune plaque trou e en traction Apr s un calcul r alis avec STAT NON LINE on souhaite tracer l eff
11. conque plus de d tails sur ce qu est une table on pourra consulter la documentation de IMPR TABLE U4 91 03 Les valeurs peuvent tre extraites en utilisant les commandes suivantes e RI ECU FONCTION U4 32 03 produit une fonction partir d un r sultat d un champ d une table Exemple volution temporelle d une composante d un champ en un point particulier POST_RELEVE_T U4 81 21 produit une table partir d un r sultat ou d un champ On peut extraire une quantit associ e aux composantes d un champ une composante un invariant en certains points particuliers ou le long d un chemin non n cessairement rectiligne MACR LIGN COUPE U4 81 13 produit une table partir d un r sultat ou d un champ le long d une ligne de coupe ligne droite compos e de segments r guliers RI d peut aussi extraire une table de quelques structures de donn es particuli res C ECU TABII E U4 71 02 produit une table partir des valeurs d un ou plusieurs param tres un r sultat Par exemple l volution du param tre de pilotage au cours d un calcul On REA CHAM P U4 72 04 permet d extraire un champ d une structure de donn es r sultat Ceci peut tre utile quand une commande ne sait pas traiter certains r sultats On peut par exemple ensuite r cup rer une fonction via RECU FONCTION CHAM GD Manuel d utilisation Fascicule u2
12. e type repe en r sultat dans astk cf U1 04 00 On proc de ensuite ainsi dans le fichier de commandes en PAR LOT NON dans POURSUIT repe Local 4 Resultats courbes OJE PJ d finition des 50 groupes de n uds de d pouillement Lgrno point01 point0O2 points50 pour chaque n ud de d pouillement for point in lgrno unit 29 les fichiers DEPL point0i dat seront recopi s dans Resultats courbes EFI FICHIER UNITE unit FICHIER REPE OUT DEPL point dat jsl kd tab POST_RELEVE T ACTION F INTITULE VomMises RESULTAT resM NOM CHA SIEQ ELNO NOM CMP VMIS TOUT ORDRE OUI GROUP NO point OPERATION EXTRACTION IMPR_TABLE Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Ast Pr oae _ASier default Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 9 10 Responsable Mathieu COURTOIS Cl U2 51 02 R vision 10626 UNIT unit TABLI tab FORMAT XMGRACE NOM PARA INST VMIS EYE on lib re l unit pour la r utiliser pour la courbe suivante DEFI FICHIER UNITE unit ACTION LIBERER 5 Quelques a
13. i sin x lt 0 y 5 sinon On cherche donc construire deux listes de r els x et ly Comme toujours en Python il y a plusieurs moyens de faire plus ou moins simplement plus ou moins l gamment x0 10 pas 20 100 1x ou bien lx x0 i pas for i in range 101 for i in range 101 l1x append x0 i pas def f x if sin x lt O return 1 5 else return 5 ly map f 1x qui applique la fonction f tous les l ments de lx On peut tracer cette courbe en utilisant les mots cl s ABSCISSE et ORDONNEE de IMPR_FONCTION cf 4 3 2 Avec des tableaux Numeric La manipulation des donn es sous forme de tableaux Numeric est simplifi e par l utilisation d op rations tr s performantes appel es ufunc sur le tableau entier dans l exemple suivant on utilise sin Numeric g re galement les tableaux plusieurs dimensions En reprenant l exemple pr c dent from Numeric import 1x arrayrange 10 10 0 2 0 2 Float ly array map f 1x ou bien sans utiliser f ly 1 5 less sin lx 0 5 1 less sin lx 0 Notons que map retourne une liste et non un tableau Numeric La deuxi me expression est entre 10 et 20 fois plus rapide sur de tr s gros tableaux 10 10 termes ce qui est cependant assez peu souvent le cas des fonctions ou tables issues d Aster On peut tracer cette courbe en utilisant les mots cl s ABSCISSE et ORDONNEE de IMPR FO
14. ort r sultant de traction en fonction du d placement vertical moyen de la partie sup rieure de l prouvette D tails du post traitement POURSUITE SOLNL2 CALC CHAMP Calcul des forces nodales reuse SOLNL2 RESULTAT SOLNL2 OPTION FORC NODA M DEFI GROUP D finition du groupe de n uds de reuse M d pouillement la ligne sup rieure de la plaque AILLAGE M CREA GROUP NO F GROUP MA LFG NOM LIGNE tab POST RELEVE T ACTION _F INTITULE Umoyen Relev du d placement moyen tous les RESULTAT SOLNL2 instants du calcul NOM CHA DEPL NOM CMP artip TOUT ORDRE OUl GROUP NO LIGNE OPERATION MOYENNE LA _F INTITULE Fresultante Relev de l effort r sultant RESULTAT SOLNI2 NOM CHA FORC NODA TOUT ORDRE OUl GROUP NO LIGNE RESULTANTE DY OPERATION EXTRACTION Manuel d utilisation Fascicule u2 51 Outils de post traitement Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Version default Titre Trac de courbes avec Code Aster Date 01 03 2013 Page 7 10 Responsable Mathieu COURTOIS IMPR TABL E TABL E tab Cl U2 51 02 R vision 10626 Juste pour rep rer le nom des param tres Contenu partiel limit aux 3 premiers instants de la table INTITULE NOEUD RESU NOM _CHAM NUME ORDRE INST DY SIYY QUANTITE Umoyen SOLNL
15. stuces utiles On propose ici quelques manipulations des donn es des tables en Python r currentes lorsque de l on veut aller plus loin dans la g n ration de courbes depuis Code Aster Extraire de la table issue de POST RELEVE T la liste des n uds de post traitement Lorsque l on post traite une grandeur sur un groupe de n uds plus d un n ud pour plusieurs instants les n uds apparaissent pour chaque instant il est donc n cessaire d extraire la liste de ces n uds sans r p tition tabpy tab EXTR TABLE Cr ation de l objet Table Python tno tabpy NOEUD values On extrait les valeurs de la colonne NOEUD lno list set i strip for i in tno M thode performante pour liminer les doublons utilisant set lno sort Tri par ordre croissant Construire un ou plusieurs mots cl s dynamiquement Ceci peut notamment tre utile pour renseigner le mot cl facteur FILTRE d IMPR TABLE en fonction du contexte on construit dans ce cas un dictionnaire qui est ensuite fourni en argument de la commande Ceci IMPR TABLE UNITE unit TABLE tab FORMAT XMGRACE FILTRE F NOM PARA NOEUD VALE K N4 _F NOM PARA NOEUD VALE K N4 NOM PARA INST VMIS est quivalent cela mfac FILTRE NOM PARA NOEUD VALE K N4 NOM PARA INTITULE VALE K exemple
Download Pdf Manuals
Related Search