Discussion :

[manumanu12]

Bonjour à tous,

je souhaite créer un script python permettant d'extraire des infos d'un odb (fichier résultat abaqus).

Concrètement, je veux sortir par face géométrique la somme des flux radiatifs (RADFLA) présents sur chaque facette d'élément de cette face géométrique.

Je peux lister les surfaces géométriques (odb.rootassembly.surfaces) et les valeurs de RADFLA (odb.steps.frames.fieldoutputs(radfla).values.data)
Mais je suis incapable de sortir la liste des faces d'éléments correspondant à la liste des RADFLA et de sortir la liste des faces d'éléments associés à chaque surface géométrique.

Quelqu'un a-t-il déjà fait du script sous Abaqus et peut m'aider??

Merci d'avance!!

[Julien N]

Bonjour,

Je peux lister les surfaces géométriques (odb.rootassembly.surfaces) et les valeurs de RADFLA (odb.steps.frames.fieldoutputs(radfla).values.data)

Est-ce que cela veut dire qu'il existe une API python pour Abaqus? Si oui, c'est du côté de sa doc qu'il faut chercher. Sinon, est-ce tu pourais montrer un exemple representatif (mais court) d'un fichier de sortie odb?

J

[manumanu12]

Bonjour Julien,

Effectivement, il existe une doc mais pas aussi détaillée que je le souhaite (je suis peut être aussi trop mauvais !!)

[Julien N]

Salut,

Si la doc est incomplete, dans un cas similaire j'essaye de "lire" le code de L'API. Sans chercher pour autant à décripter en detail tout le code, je commence par determiner comment celui-ci est contruit, puis je repère les elements que je connnais déjà (ici pour odb.rootassembly.surfaces et odb.steps.frames.fieldoutputs(radfla).values.data, je rechercherais l'objet odb, puis la methode rootassembly, etc). J'utilise aussi abondamment la methode __dict__ qui liste les attributs d'un objet

Comme ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

print odb.rootassembly.__dict__

Bref, si tu veux résoudre ton probleme, tu vas être obligé de fouiller. Il y a très peu de chance que quelqu'un ici utilise aussi cette API (mais on ne sait jamais!!)

[EDIT]
La fonction dir() liste les methodes:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

print dir(odb.rootassembly)

J

[manumanu12]

Bonsoir Julien,

C'est effectivement ce que j'ai tenté de faire en passant par l'éditeur présent dans Abq viewer, où je passais les différents attributs.
Je ne connaissais pas _dict_ par contre, merci pour le conseil.

Voici l'ébauche du scripttest.py

[Julien N]

Salut Manu,

Malheureusement je n'ai pas Abaqus, ni sur mon PC ni au boulot. Donc je ne peux pas vraiment tester ton script. Cependant je peux te conseiller sur certains morceaux du code:

Par exemple ici:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
liste_surfaces=[]
for surface in myassembly.surfaces.keys():
  liste_surfaces.append(surface)
 
nb_surfaces=len(liste_surfaces)

i=0
while i < nb_surfaces:
  print liste_surfaces[i] 
  i+=1

Je suppose que tu souhaites récupérer les surfaces de ton modèle et les afficher. Tu peux tout simplement faire ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
liste_surfaces = myassembly.surfaces.keys()
print liste_surfaces

La méthode keys() du dictionnaire par défaut sort une liste des clés. Tu peux directement les afficher. Si tu souhaites uniquement les afficher tu peux tout aussi bien faire:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

print myassembly.surfaces.keys()

Les compter est tout aussi simple:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

print len(myassembly.surfaces.keys())

A ce que je comprends, les surfaces sont constituées de faces. Et tu aimerais également les récupérer. Est-ce que ces faces sont les "valeurs" associées aux "clés" de myassembly.surfaces? Est-ce que tu peux tester ceci et poster le résultat:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
for key, value in myassembly.surfaces.iteritems():
    print "clé: ", key, type(key)
    print "valeur: ", value, type(value)

Ju

[manumanu12]

Salut Julien,

Effectivement mon code n'est pas forcément optimal!

Je n'ai pas Abaqus ici, je testerai seulement lundi...

[manumanu12]

Salut,

Je n'ai pas pu faire fonctionner ton bout de code, tant pis!

Je me suis aussi rapproché du support abaqus, qui m'a conseillé de revoir la méthode car il n'est pas possible de lier surfaces géométriques et faces d'éléments...
Je vais essayer de m'en sortir par le JDD, je vais regarder ça dans la semaine.

Merci en tout cas pour ton aide et tes conseils!

[Julien N]

De rien, et bon courage.