Discussion :

[amandine2402]

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 :

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

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()

[VinsS]

Salut,

Quelle différence entre ceci:

Code :

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

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

et ceci:

Code :

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

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 :

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

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

Qu'est-ce que ceci:

Code :

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

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 :

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

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 :

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

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 :

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

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.

[amandine2402]

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 :

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

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?

[VinsS]

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