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

最新推荐文章于 2024-06-11 20:29:50 发布

CSDNJBG

最新推荐文章于 2024-06-11 20:29:50 发布

阅读量8.9k

点赞数 4

分类专栏： abaqus开发文章标签： python

本文链接：https://blog.csdn.net/CSDNJBG/article/details/111636565

版权

本文介绍了如何利用Python脚本在ABAQUS中提取单元的质心应力和节点坐标信息，详细阐述了相关步骤和代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

# -*- 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)