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

[二次开发] 单元积分点的主应力、应力分量可以提取;节点的可以吗?

[复制链接]
发表于 2022-10-12 16:07:55 | 显示全部楼层 |阅读模式 来自 四川
如题

请问各位前辈、老师:

总所周知,abaqus可以用python提取单元积分点的三个主应力以及6个应力分量。程序如下:
  1. # -*- coding: mbcs -*-
  2. from odbAccess import *
  3. from abaqusConstants import *
  4. from textRepr import *  # 为了使用prettyPrint()函数,这是abaqus内置的函数

  5. o = session.openOdb(name='myJob.odb', readOnly=False)
  6. session.viewports['Viewport: 1'].setValues(displayedObject=o)
  7. steps = o.steps
  8. print steps  # 除了初始分析步,只有一个分析步,是字典形式
  9. step1 = steps[steps.keys()[0]]
  10. frame = step1.frames  # 该分析步的所有帧,是个list
  11. print frame, '\n', len(frame)

  12. f_last_step1 = frame[-1]  # 最后一帧
  13. prettyPrint(f_last_step1)  # 漂亮的打印
  14. f_last_outputs = f_last_step1.fieldOutputs  # 获取最后一帧的结果
  15. prettyPrint(f_last_outputs)  # {'S': 'FieldOutput object', 'U': 'FieldOutput object'}
  16. # 最后一帧应力
  17. fop_S = f_last_outputs['S']
  18. print '======================================================='
  19. # 最后一帧位移
  20. fop_U = f_last_outputs['U']
  21. prettyPrint(fop_S.locations[0])  # locations为记录数据依附点的FieldLocation对象序列,就一个值,告诉你在积分点
  22. print '======================第一个点的数据==========================='
  23. prettyPrint(fop_S.values[0])

  24. print len(fop_S.values)  # 3300个单元
  25. prettyPrint(fop_S.values[0].data)  # 6个应力分量 S11/S22/S33/S12/S13/S23
  26. prettyPrint(fop_S.values[0].maxPrincipal)  # 第一主应力
  27. prettyPrint(fop_S.values[0].instance)
复制代码



问题:那么,单元节点可以通过类似于编程提取积分点数据的方式,一次性提取6个分量与3个主应力吗?
           可以通过python程序进行提取,保存在程序内存中吗?因为我需要节点的应力数据进行后续计算,并不想产生额外的文件。

还是说节点的主应力与应力分量只能通过可视化模块->查询->导出csv文件查看?如图:


自我思考:

我知道可以通过以set建立path(路径)的方式提取,但是set的建立又有另外的问题(请见我上一个提问帖,如能解答,感激不尽)而且一次只能提取一个量(一个主应力或一个应力分量),而且程序需要多次遍历,时间成本过大;
也能通过建立场变量xy数据的方式提取,但是这样会产生大量xydata,因为abaqus会将每一个点的每一个数据都新建一个xydata,而不是形成大列表,这样真的太蠢了,而且我不需要真的创建xydata,我只是编程计算中需要数据。
如图:



求前辈、老师们给予指导与解答,不胜感激!




本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

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

本版积分规则

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

GMT+8, 2024-4-28 21:09 , Processed in 0.026008 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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