笔者为科研界最后的摆烂王,目前利用python代码对Abaqus进行二次开发尚在学习中。欢迎各位摆烂的仁人志士们和我一起摆烂!ps:搞什么科研,如果不是被逼无奈,谁要搞科研!
该代码是学习过程中,对前人已有工作所做的稍加修改。为什么是稍加修改,是因为原代码跑不出来!!笔者在提取场输出的位移数据时,渴望偷懒,打算百度一下草草了事,奈何发现网上代码多半驴头不对马嘴,笔者明明是想提取位移,而不是节点和单元的集合!!所以被逼无奈之下,只好硬着头皮修改!欢迎各位大佬们把小弟代码更优化,然后也发给小弟,让小弟尝尝被带飞的滋味!万分感谢,给您老们拜个大年,祝各位新年新气象,文章多多,money多多!!
首先,利用python脚本对Abaqus进行数据提取时,要对python脚本做前步的导入处理。
第一个是一定要给出python脚本遍历查询的路径提示,该代码如下:
#!/usr/bin/python# -*- coding:UTF-8-*-
第二个是在调用odb数据文件时,要导入Abaqus内置的模块,代码如下:
其中 from odbAccess import* 是调用odb数据文件必须要用的。
from odbAccess import*from abaqusConstants import *from odbMaterial import *from odbSection import *import csvimport string
随后的代码编写则按照以下思路进行:
(1)打开odb对象文件
(2)创建变量表示第一个分析步
(3)创建变量表示第一个分析步的最后一帧
(4)创建变量表示节点集
(5)创建变量表示2,3步中输出的位移
其代码如下:
通过该代码便完成了对已生成的odb文件中的特定节点集合进行的位移提取
笔者为了验证代码的正确性,做了简单的弹性材料的梁受压模拟,利用python脚本遍历循环输出的结果如下:
若想对输出的数据再进一步优化,可以利用open()与write()等函数进行二次处理,将输出的数据直接写入到文档中方便使用。