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

[后处理] abaqus后处理提取应力应变曲线

[复制链接]
发表于 2017-3-14 15:33:41 | 显示全部楼层 |阅读模式 来自 陕西西安
悬赏25仿真币未解决
模拟再生混凝土裂缝开展问题,选用的cdp本构模型,加位移荷载,在后处理中,请问如何获取全应力应变关系曲线?(五相复合材料,本人只会提取单元或者节点的应力应变)

发表于 2017-11-7 21:42:36 | 显示全部楼层 来自 陕西宝鸡
Simdroid开发平台
from odbAccess  import*

#
# get field
#

odb = openOdb(path='Job-shikeli.odb')
framesnumber=len(odb.steps['Step-1'].frames)
print "framesnumber=",framesnumber
i=0
while i<framesnumber:
        endSet = odb.rootAssembly.elementSets['WHOLE']   #zhe ge dan yuan ji he shi step li mian de shu chu de dan yuan ji he
        field1=odb.steps['Step-1'].frames[i].fieldOutputs['S']
        field2=odb.steps['Step-1'].frames[i].fieldOutputs['EVOL']
        subField1=field1.getSubset(region=endSet)
        subField2=field2.getSubset(region=endSet)
        sum_mises=0
        vol=0
        for val_2 in subField2.values:
            for val_1 in subField1.values:
                if val_2.elementLabel == val_1.elementLabel:
                       
                       mises = val_2.data * val_1.mises
                       break
            sum_mises = sum_mises + mises
            vol = vol + val_2.data
            ave_mises= sum_mises/vol
        
        print ave_mises
        i+=1
提取应力的,希望能帮到你
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-4-20 22:24 , Processed in 0.026624 second(s), 9 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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