py源代码|平均应力应变位移输出至Excel

当提取模型文件计算结果数据量较大时，依次设置历程输出任务繁重。实际上结合公众号(有限元与力学)提供的脚本文件(output_constants.py和output_main.py)，只需几行代码即可高效提取有用数据。

话不多说实现方式如下：（不需要下载任何软件）

（1）以记事本形式打开output_constants.py文件，修改odbName和stepName

odbName='RS1'        #odb文件名
stepName='Step-1'    #step名

（2）在odb文件同目录下新建一个txt文件，将以下内容作为表头固定输入

#!/usr/bin/python
#-*-coding:utf-8-*-
import csv
import output_main

（3）以提取模型的分析时间、S11应力和E11应变数据为例，输入下列代码

time=[]     #储存时间数据
time=output_main.output_time(time)

stress=[]          #储存应力数据,其中'CONCRETE-1'部件名，'SET-1'单元集合名
stress=output_main.output_S11('CONCRETE-1','SET-1',stress)

strain=[]          #储存应力数据,其中'CONCRETE'部件名，'SET-1'单元集合名
stress=output_main.output_E11('CONCRETE-1','SET-1',strain)

（4）输出数据至Excel

New=open('Output.csv','wb')   #新建一个Excel文件储存数据
New.write('time,stress,strain\n')  #输入Excel表头分别为time，stress，strain
for i in range(len(stress)):   #循环输出数据
    New.write('%s,%s,%s\n'%(time[i],stress[i],strain[i]))
    New.close()

（5）保存新建txt文件，将.txt后缀修改为.py，在Abaqus中以脚本形式运行该文件

output_main.output_S11中的S11用于控制输出结果

结果控制参数如下

S11—x方向应力；S22—y方向应力；S33—z方向应力；mises—mises应力；Smax—最大主应力；Smin—最小主应力；E11—x方向应变；E22—y方向应变；E33—z方向应变；Emax—最大主应变；U11—x方向位移；U22—y方向位移；U33—z方向位移；RF1—x方向支反力；RF2—y方向支反力；RF3—z方向支反力

output_constants.py、output_main.py获取方式如下

公众号(有限元与力学)任意两篇文章转发朋友圈12h后截图加客服(QQ:416140933)索取