许可优化
许可优化
产品
产品
解决方案
解决方案
服务支持
服务支持
关于
关于
软件库
当前位置:服务支持 >  软件文章 >  提取Abaqus中节点集合的应力应变

提取Abaqus中节点集合的应力应变

阅读数 3
点赞 0
article_banner
节点是没有应力应变的,只有积分点才有,所以直接在odb文件里面用FieldOutput是整不出来的(结果是个空列表),只能用session这种让abaqus自己帮你处理。具体的代码需要根据自己操作的时候abaqus.rpy这个文件中代码的变化来更改
import numpy
from abaqus import *
from abaqusConstants import *
from caeModules import *
import visualization
from driverUtils import executeOnCaeStartup
import numpy as np
from driverUtils import executeOnCaeStartup
from odbAccess import *
import numpy as np

job_name = "v1100-2"
o1 = session.openOdb(name='https://www.gofarlic.com/Abaqus/file/'+job_name+'.odb')
odb = session.odbs['https://www.gofarlic.com/Abaqus/file/'+job_name+'.odb']

# F
session.xyDataListFromField(odb=odb, outputPosition=NODAL, variable=(('CF',
                                                                      NODAL, ((COMPONENT, 'CF1'), (COMPONENT, 'CF2'),
                                                                              (COMPONENT, 'CF3'),)),),
                            nodeSets=("RP_F",))
# c0 = session.Curve(xyData=xy0)
xy1 = session.XYDataFromHistory(name='RF2 PI: rootAssembly N: 1 NSET RP_F-1',
                                odb=odb,
                                outputVariableName='Reaction force: RF2 PI: rootAssembly Node 1 in NSET RP_F',
                                steps=('Step-1',), __linkedVpName__='Viewport: 1')
# c1 = session.Curve(xyData=xy1)
xy2 = session.XYDataFromHistory(name='RF3 PI: rootAssembly N: 1 NSET RP_F-1',
                                odb=odb,
                                outputVariableName='Reaction force: RF3 PI: rootAssembly Node 1 in NSET RP_F',
                                steps=('Step-1',), __linkedVpName__='Viewport: 1')

# c1
session.xyDataListFromField(odb=odb, outputPosition=NODAL,
                            variable=(('LE', INTEGRATION_POINT, ((COMPONENT, 'LE22'),)),),
                            nodeSets=("C1_LE22_P8",))
xy1 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 575']
xy2 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2501']
xy3 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2506']
xy4 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2511']
xy5 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2516']
xy6 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2521']
xy7 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2526']
xy8 = session.xyDataObjects['LE:LE22 (Avg: 75%) PI: V11-1 N: 2531']
xy9 = sum(xy1, xy2, xy3, xy4, xy5, xy6, xy7, xy8) /
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删


相关文章
QR Code
微信扫一扫,欢迎咨询~
customer

online

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

* 公司名称:

姓名不为空

姓名不为空

姓名不为空
手机不正确

手机不正确

手机不正确
公司不为空

公司不为空

公司不为空