Python提取Abaqus结果数据

View Post

1. Python提取Abaqus结果数据
参考资料

1. Python提取Abaqus结果数据

1.1 Python选择节点和单元

getByBoundingBox()
getByBoundingCylinder()
getByBoundingSphere()
可用tab键查找函数的变量

通过位置选择节点

p = mdb.models[\'Thermal2D-L\'].parts[\'Part-2\']
n = p.nodes
nlist = n.getByBoundingBox(-100,-100,0,100,100,0)#通过矩形来框选节点
p.Set(nodes=nlist, name=\'Set-44\')#创建set

通过位置选择单元

p = mdb.models[\'Thermal2D-L\'].parts[\'Part-2\']
e = p.elements
elist = e.getByBoundingBox(-0.2,-12.6,-0.1,0.2,12.6,0.1)
p.Set(elements=elist,name=\'Set-5\')

注：选中节点的单元都会被选中

获取节点的label

p = mdb.models[\'Thermal2D-L\'].parts[\'Part-2\']
e = p.elements
elist = e.getByBoundingBox(-0.2,-12.6,-0.1,0.2,12.6,0.1)
e1 = elist[1]
elabel1 = e1.label

1.2 Python提取结果

结构都存储在*.odb文件中，可以采用python直接进行读取

提取位移

from odbAccess import*
from abaqusConstants import*
import os

myodb=openOdb(path=r\'H:\Abaqus\MultiScaleMeshSize\Job-2D-SubThermal.odb\')
RS=myodb.steps[\'Step-1\'].frames[1].fieldOutputs[\'U\'].values
cpFile=open(\'distance.txt\',\'w\')
for i in RS:
    cpFile.write(\'%d %7.4f %7.4f\n\' % (i.nodeLabel,i.data[0],i.data[1]))
    
else:
    cpFile.close()
    
# frames为第几个增量
# fieldOutputs[\'U\']结果类型，位移

提取mises应力

RF=myodb.steps[\'Step-1\'].frames[1].fieldOutputs[\'S\'].values
RF0 = RF[0]
mises = RF0.mises #mises应力
data = RF0.data #S11,S22,S33,S12
eleType = RF0.baseElementType #单元类型
eleLabel = RF0.elementLabel
maxPrincipal = RF0.maxInPlanePrincipal #最大平面主应力，可参考具体

获取节点集合上的位移

from odbAccess import*
from abaqusConstants import*
import sys
import os

odb=openOdb(path=r\'H:\Abaqus\MultiScaleMeshSize\Job-2D-Thermal-Displace.odb\')

lastFrame= odb.steps[\'Step-1\'].frames[-1]#创建变量lastFrame，得到载荷步Step-1的最后一帧
displacement= lastFrame.fieldOutputs[\'U\']#创建变量displacement，得到最后一帧的位移场数据

#创建变量center，得到节点集PUNCH
#其中PART-2-1为装配体中的几何体，SET-1为创建的集合
center =odb.rootAssembly.instances[\'PART-2-1\'].nodeSets[\'SET-1\']

#创建变量centerDisplacement，得到region center的位移场数据
centerDisplacement= displacement.getSubset(region=center)

#输出各种信息，该节点集只有一个节点
for v in centerDisplacement.values:
    print \'Position = \', v.position,\'Type =\',v.type
    print \'Node label = \', v.nodeLabel
    print \'X displacement = \', v.data[0]
    print \'Y displacement = \', v.data[1]
    print \'Displacement magnitude =\',v.magnitude

odb.close()

获取单元的应力值

from odbAccess import*
from abaqusConstants import*
import sys
import os

# cpFile=open(\'artlcF1.txt\',\'w\')
odb=openOdb(path=r\'H:\Abaqus\MultiScaleMeshSize\Job-2D-Thermal-Displace.odb\')

# topCenter =odb.rootAssembly.instances[\'PART-2-1\'].nodeSets[\'SET-1\']#创建变量topCenter

topCenter =odb.rootAssembly.instances[\'PART-2-1\'].elementSets[\'SET-4\']
#创建变量 stressField，得到位移场数据
stressField= odb.steps[\'Step-1\'].frames[-1].fieldOutputs[\'S\']

field =stressField.getSubset(region=topCenter)

#创建变量 fieldValues，得到field的数值数据
fieldValues= field.values

#输出label值
for v in fieldValues:
    print(v.elementLabel)

# 输出应力值
for v in fieldValues:
    print(v.mises)

# cpFile.write (\'%-15.5f\' % component)
# cpFile.close()

获取节点的温度值

from odbAccess import*
from abaqusConstants import*
import sys
import os

odb=openOdb(path=r\'H:\Abaqus\MultiScaleMeshSize\Job-2D-Thermal-Displace.odb\')

lastFrame= odb.steps[\'Step-1\'].frames[-1]#创建变量lastFrame，得到载荷步Step-1的最后一帧
displacement= lastFrame.fieldOutputs[\'NT11\']#创建变量displacement，得到最后一帧的位移场数据

#创建变量center，得到节点集PUNCH
#其中PART-2-1为装配体中的几何体，SET-1为创建的集合
center =odb.rootAssembly.instances[\'PART-2-1\'].nodeSets[\'SET-1\']

#创建变量centerDisplacement，得到region center的位移场数据
centerDisplacement= displacement.getSubset(region=center)

#输出各种信息，该节点集只有一个节点
for v in centerDisplacement.values:
    print \'Node label = \', v.nodeLabel
    print \'templabel= \', v.data

1.3 Python结果提取方法

可以直接提取所有的单元或者几何节点的值
首先获取节点或者单元的label值，然后再读取响应的结果
在求解先，先建立要提取结果的set集合，然后在结果中按照set进行提取

1.4 单元面积计算方法

首先获取组成单元的节点
然后获取节点的坐标
通过节点坐标来计算面积

p = mdb.models[\'Thermal2D-L\'].parts[\'Part-2\']
n = p.nodes
e = p.elements

e1 = e[1]
nodeIndexList = e1.connectivity

# 获取节点的坐标
PointList = []
for Index in nodeIndexList:
    temp = n[Index].coordinates
    PointList.append(temp)
    
# 计算多边形的面积

参考资料

【1】Abaqus Python getByBoundingBox command
【2】Python提取Abaqus结果数据
【3】Abaqus 6.14 帮助文档