我试图在每个节点或每个集成点处提取应力。 我可以得到压力,但我不知道它的价值在哪里。 如何添加节点号或节点坐标?下面是我的代码,只提取所有应力分量
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()
答案 0 :(得分:2)
按数字索引在节点上循环:
for i in range (len(....nodes)):
N=....nodes[i]
找到每个节点的压力