找回密码
 注册
Simdroid-非首页
查看: 658|回复: 10

[后处理] python提取节点集应力

[复制链接]
发表于 2020-2-15 16:09:12 | 显示全部楼层 |阅读模式 来自 四川资阳
悬赏1仿真币未解决
  最近使用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()

发表于 2020-2-17 04:34:45 | 显示全部楼层 来自 英国
Simdroid开发平台
正常情况下你只能获得单元积分点上的应力值和节点上的位移值。
但是在ABAQUS CAE 中的确是可以通过设立Path的办法来获得节点上的应力值的。
回复

使用道具 举报

发表于 2020-2-21 15:14:59 | 显示全部楼层 来自 台湾
field =stressField.getSubset(region=topCenter,position=INTEGRATION_POINT, elementType ='C3D8R')
上面這行的position我記得可以換成ELEMENT_NODAL或是NODAL?
回复

使用道具 举报

 楼主| 发表于 2020-2-22 18:20:03 | 显示全部楼层 来自 四川资阳
duncanusn 发表于 2020-2-21 15:14
field =stressField.getSubset(region=topCenter,position=INTEGRATION_POINT, elementType ='C3D8R')
上面 ...

嗯嗯 好的 就是这个问题 现在有结果了 多谢
回复

使用道具 举报

发表于 2021-2-10 12:04:19 | 显示全部楼层 来自 中国
打扰啦~ 请问您怎么建的节点集合Set-1吗?
回复

使用道具 举报

发表于 2021-5-6 17:00:00 | 显示全部楼层 来自 江苏南京
yr1998 发表于 2021-2-10 12:04
打扰啦~ 请问您怎么建的节点集合Set-1吗?

可以用getByBoundingBox或Cylinder
回复

使用道具 举报

发表于 2021-5-6 18:00:15 | 显示全部楼层 来自 安徽淮南
duoxielouzhufenxiang
回复

使用道具 举报

发表于 2021-5-12 22:10:25 | 显示全部楼层 来自 甘肃兰州
猴子oooooo 发表于 2020-2-22 18:20
嗯嗯 好的 就是这个问题 现在有结果了 多谢

可以加一下您联系方式请教一下吗?qq:1519505084
回复

使用道具 举报

发表于 2021-5-18 09:06:33 | 显示全部楼层 来自 甘肃兰州
你好!打扰了,请问提取这个方向的值V,V2,代码该怎么写。
回复

使用道具 举报

发表于 2021-5-18 09:15:50 | 显示全部楼层 来自 甘肃兰州
v2Field=lifangtiodb.steps['Step-2'].frames[-1].fieldOutputs['V'],请问怎么提V2方向的?
回复

使用道具 举报

发表于 2023-8-19 09:34:57 | 显示全部楼层 来自 浙江宁波
谢谢吧友 救大命
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-20 01:00 , Processed in 0.036900 second(s), 10 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表