Script python pour extraire les déplacements à partir d'Abaqus

Bonjour,

Dans le cadre d'un projet, je dois extraire la déformée modale ( et donc les déplacements) d'un modèle abaqus pour chaque mode. Pour cela, j'ai créé un 'Path' dans abaqus de manière à sélectionner les déplacements selon une ligne de ma poutre. Cependant, j'ai déjà essayé différents codes Python afin d'extraire ses données mais je n'y arrive pas... Est ce que quelqu'un aurait une idée svp ?

Merci d'avance pour ceux qui prendront le temps d'y regarder

J'avais dans un premier temps fait le code suivant mais il semble y avoir plusieurs erreurs.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

from odbAccess import * from abaqus import * from abaqusConstants import * import types import visualization from math import * from odbMaterial import * from odbSection import * frame_k=[] odb=openOdb(path='demoplaque.odb') # opens database with specified name historyRegion=('Path-1') historyOutput=['U'] for k in range (len(odb.steps['modal'].frames)): frame_k.append(odb.steps['modal'].frames[k]) frame_k.fieldOutputs['U'] fichier = open('results.rpt','w') i=0 for data in frame_k: if str(data) != '[]': fichier.write(str(data)) i=i+1 fichier.write('\n') fichier.close()

---

Salut,

Quelle différence entre ceci:

Code:

---

1 2 3

for k in range (len(odb.steps['modal'].frames)): frame_k.append(odb.steps['modal'].frames[k])

---

et ceci:

Code:

---

1 2	frame_k = odb.steps['modal'].frames

---

Essaye pour voir ...

Utilise des print pour t'assurer de la validité des étapes de ton code:

Code:

---

1 2	print('Nombre de frames: %s' % len(frame_k))

---

Qu'est-ce que ceci:

Code:

---

1 2	frame_k.fieldOutputs['U']

---

frame_k est une liste, les listes n'ont pas d'attribut fieldOutputs et que penser du "['U']" ?


Pour finir:

Code:

---

1 2 3 4 5 6 7 8 9

fichier = open('results.rpt','w') i=0 # Inutile ! for data in frame_k: if str(data) != '[]': # ?? fichier.write(str(data)) i=i+1 # Inutile ! fichier.write('\n') fichier.close()

---

Est-ce que ce sont les listes vides que tu veux supprimer ?

Il est souvent plus simple de préparer ses objets avant de les utiliser, avec l'avantage d'un code plus clair pour tout le monde.

Par exemple:

Code:

---

1 2 3 4 5 6 7

# Supprimer les listes vides frames = [str(f) for f in frame_k if f] # ! Supprimera aussi les 0 ! # Lier les données avec un newline text = '\n'.join(frames) with open('results.rpt','w') as outf: outf.write(text)

---

pas de close() avec with open()

Demo:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

>>> l = [25, [], 'hfgtr', [64, 72], 'kjhg'] >>> f = [str(i) for i in l if i] >>> f ['25', 'hfgtr', '[64, 72]', 'kjhg'] >>> txt = '\n'.join(f) >>> print(txt) 25 hfgtr [64, 72] kjhg

---

D'autre part, si ton code lève une erreur, copie colle le traceback complet avec ton message.

Tout d'abord un très grand merci déjà pour ta réponse. Pour le 1er changement, ce que tu fais fonctionne donc j'ai gardé cet syntaxe , merci .
Ensuite, pour la ligne

Code:

---

frame_k.fieldOutputs['U']

---

J'essayais d'afficher les déplacements selon y d'un Path-1 que j'ai défini dans Abaqus, aurais tu une idée de comment faire ca stp ? Cela peut sembler facile mais je n'y arrive pas.
Et également , aurais tu une idée de comment je peux extraire le coefficient d'amortissement?

Hélas, je ne connais rien à Abaqus.

J'ai regardé sur leur site, c'est pas open source, il y a bien une version student téléchargeable (c'est peut-être celle-là que tu utilises) mais je n'ai vu aucune mention d'OS or moi je suis sous Linux, pas Windows.

Tu as vu qu'il existe ce PDF particulièrement complet (en anglais):
http://www.maths.cam.ac.uk/computing...CRIPT_USER.pdf