Abaqus利用region读取结果数据

Abaqus计算完成后,为输出某部分模型的结果数据,需要用到region,并且配套使用的是getSubset函数。

示例一创建一个变量center,将一个叫做PUNCH的节点set赋值给这个变量,通过整个模型的位移场displacement得到该节点set的输出数据,代码如下:

center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']

centerDisplacement = displacement.getSubset(region=center)

centerValues = centerDisplacement.values

for v in centerValues:

   print v.nodeLabel, v.data



其中的displacement是已经创建好的位移场数据,centerDisplacement变量得到该节点集的所有数据,centerValues变量得到该节点集的数值数据,通过循环输出节点号和当前数据。

示例二如下:

#创建变量topCenter,得到单元集CENT

topCenter =odb.rootAssembly.instances['PART-1-1'].elementSets['CENT']

#创建变量stressField,得到位移场数据

stressField = odb.steps['Step-2'].frames[3].fieldOutputs['S']

#创建变量field,得到该单元集的积分点的数据

field = stressField.getSubset(region=topCenter,

   position=INTEGRATION_POINT, elementType = 'CAX4')

#创建变量fieldValues,得到field的数值数据

fieldValues = field.values

#循环输出单元号,单元上的积分点号以及积分点的应力分量

for v in fieldValues:

   print 'Element label = ', v.elementLabel,    

   if v.integrationPoint:

       print 'Integration Point = ', v.integrationPoint

   else:

       print

# For each tensor component.

   for component in v.data:

# Print using a format. The comma at the end of the

# print statement suppresses the carriage return.

       print '%-10.5f' % component,

# After each tuple has printed, print a carriage return.

   Print



最终的输出结果:

Element label =  1 Integration Point =  1

S : 0.01230    -0.05658   0.00892    -0.00015

Element label =  1 Integration Point =  2

S : 0.01313    -0.05659   0.00892    -0.00106

Element label =  1 Integration Point =  3

S : 0.00619    -0.05642   0.00892    -0.00023

Element label =  1 Integration Point =  4

S : 0.00697    -0.05642   0.00892    -0.00108

Element label =  11 Integration Point =  1

S : 0.01281    -0.05660   0.00897    -0.00146

Element label =  11 Integration Point =  2

S : 0.01183    -0.05651   0.00897    -0.00257

Element label =  11 Integration Point =  3 ...

另:对于getSubset的position有以下几个可选:

INTEGRATION_POINT   !积分点

NODAL                !节点(平均)

ELEMENT_NODAL      !单元节点(未平均)

CENTROID             !单元质心点

免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空