Jeje3103
New member
- Oct 22, 2015
- 8
I m trying to extract stresses at each nodes or each integration points. I can get stresses but i have nothing to know where it is situated the values. How can i add the node number or node coordinate? Below it is my code to extract just all stresses components
from odbAccess import *
import numpy as np
odb = openOdb(path='C:/Users/Solal/Job-1.odb')
lastFrame = odb.steps['Step-1'].frames[-1]
location = lastFrame.fieldOutputs['COORD']
stress=lastFrame.fieldOutputs['S']
sortie = open('Resu_PE_maillage.txt', 'w')
sortie.write('\t S11 \t\t S22 \t\t S33 \t\t S12 \t\t S23 \t\t S13 \n')
for S in stress.values:
sortie.write('%.10e\t\t %.10e\t\t %.10e\t\t %.10e\t\t %.10e\t\t %.10e\n' % (S.data[0], S.data[1], S.data[2], S.data[3], S.data[4], S.data[5]))
sortie.close()
from odbAccess import *
import numpy as np
odb = openOdb(path='C:/Users/Solal/Job-1.odb')
lastFrame = odb.steps['Step-1'].frames[-1]
location = lastFrame.fieldOutputs['COORD']
stress=lastFrame.fieldOutputs['S']
sortie = open('Resu_PE_maillage.txt', 'w')
sortie.write('\t S11 \t\t S22 \t\t S33 \t\t S12 \t\t S23 \t\t S13 \n')
for S in stress.values:
sortie.write('%.10e\t\t %.10e\t\t %.10e\t\t %.10e\t\t %.10e\t\t %.10e\n' % (S.data[0], S.data[1], S.data[2], S.data[3], S.data[4], S.data[5]))
sortie.close()