最近使用Python提取节点集的应力,编辑了一段脚本,出现的问题是,运行脚本后,文本里没有输出任何结果;另外,还有个问题不是很明白:应力的提取是针对积分点上的,那么怎么才能提取积分点上的应力,怎么才能将积分点上的应力映射到节点上,即提取节点上的应力。以下是我编辑的脚本,请各位指导一下问题所在,多谢。
from odbAccess import*
from abaqusConstants import*
odb=openOdb(path='D:\Program Files (x86)\ABQUS2 Abaqustemp\shijianpython-0214\Job-1.odb')
cpFile=open('yingliF1.txt','w')
topCenter=odb.rootAssembly.nodeSets['SET-1']
stressField=odb.steps['Step-1'].frames[-1].fieldOutputs['S']
field =stressField.getSubset(region=topCenter,position=INTEGRATION_POINT, elementType ='C3D8R')
fieldValues= field.values
for v in fieldValues:
cpFile.write(v.elementLabel,)
cpFile.write('Integration Point = ',v.integrationPoint)
cpFile.write('%d %7.4f %7.4f\n' % (v.data[0],v.data[1],v.data[2]))
cpFile.write('\n')
else:
cpFile.close() |