利用python提取abaqus单元质心应力及单元节点坐标

# -*- coding: utf-8 -*- 
from odbAccess import *
import os

fileName = r"E:\Abaqus-model\HouZhuangDAM\OLD-ZC-TEMP-DOWN\old-zc-stress-up.odb"
txtName=r"E:\Abaqus-model\HouZhuangDAM\OLD-ZC-TEMP-DOWN\result.txt"

cpFile=open(txtName,'w')
#打开odb，获取应力场数据
odb = openOdb(path=fileName)
fieldOuts=odb.steps['bottomwater'].frames[-1].fieldOutputs['S']


#获取odb实例及单元
assembly=odb.rootAssembly
elementSet=assembly.elementSets['mulElements']
#getSubset参数，region指定区域（集合），position指定位置：INTEGRATION_POINT积分点，CENTROID质心，ELEMENT_NODAL单元结点,NODAL没有返回值？
stressSet=fieldOuts.getSubset(region=elementSet,position=CENTROID)

cpFile.write("print element stress\n")
#对于c3d8单元，因为有8个积分点，如果取积分点应力的话会循环8次，把各个积分点应力都列出来
for s in stressSet.values:
    cpFile.write("element label:"+str( s.elementLabel)