Python编写abaqus后处理脚本（学习笔记）

时间：2024-04-15 17:33:15

本节内容参考自书籍《Python语言在Abaqus中的应用》，注意：以下代码为伪代码，仅供参考

1、导入必要的模块，加载后处理odb文件

from abaqus import *
from abaqusConstants import *
from odbAccess import*
import os
odb=openOdb(path=\'job-1.odb\')

2、常用获取模型数据的api

#输出所有部件实例名称，在cae界面能看到结果,该部分代码只能在abaqus command环境下运行
for instanceName in odb.rootAssembly.instances.keys():
    print instanceName

#调用部件实例所有节点集名
print \'Node set=\',odb.rootAssembly.instances[\'PART-1-1\'].nodeSets.keys()

#所有单元集名
print \'Element set=\',odb.rootAssembly.instances[\'PART-1-1\'].elementSets.keys()

#输出所有材料的名称
allMaterials=odb.materials
for materialName in allMaterials.keys():
    print \'Material Name:\',materialName

#输出分析步的关键字
for stepName in odb.steps.keys():
    print stepName

#帧赋值
lastFrame=odb.steps[\'Step-1\'].frames[-1]

#读取场输出数据,列出分析步最后一帧的所有变量
for fieldName in lastFrame.fieldOutputs.keys():
    print fieldName
#输出最后一帧场输出的变量名、描述和成员类型
for f in lastFrame.fieldOutputs.values():
    print f.name,\':\',f.description
    print \'Type:\',f.type
    #对于每个计算值，输出其位置
    for loc in f.locations:
        print \'Position:\',loc.position
    print

3、输出odb文件中的位移信息

#获取WALL部分位移---节点集合
wall=odb.rootAssembly.instances[\'PART-1-1\'].nodeSets[\'WALL\']
print(odb.steps[\'water5\'])
lastFrame=odb.steps[\'water5\'].frames[1]
displacement=lastFrame.fieldOutputs[\'U\']
wall_Disp=displacement.getSubset(region=wall)
wall_U=wall_Disp.values
print(len(wall_U))

#位移输出为txt文件
full_path=\'D:\\abaqus1\\\'+\'1111111111111.txt\'
f=open(full_path,\'a+\')
for i in [400,800,1200]:
    f.write(\'%7.4f %7.4f %7.4f \' % (wall_U[i].data[0],wall_U[i].data[1],wall_U[i].data[2]))
f.write(\'\n\')
f.close()

4、输出odb文件中的应力信息

#获取防渗墙wall部分的节点集合,输出应力值
#获取防渗墙的单元集合
wall_1=odb.rootAssembly.instances[\'PART-1-1\'].elementSets[\'HMPROP_WALL1\']
stress=lastFrame.fieldOutputs[\'S\']
wall_Stress=stress.getSubset(region=wall_1)
wall_S=wall_Stress.values
print(\'wall_S number:\',len(wall_S))

# 输出应力
full_path=\'C:\\abaqus2\\\'+\'wall_S.txt\'
f=open(full_path,\'w+\')
for i in wall_S:
    f.write(\'%d %7.4f %7.4f %7.4f \n\' % (i.elementLabel,i.mises,i.maxPrincipal,i.minPrincipal))
else:
    f.close()

相关文章