我试图在每个节点或每个集成点处提取应力。我会感到压力,但我不知道它在哪里价值观。如何添加节点号或节点坐标?下面是我提取所有应力分量的代码
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()