abaqus批命令使用心得

这段时间因为一个项目工作，对abaqus的批命令使用较多，也有了一点点的心得，借此平台记录，以备将来借鉴。
首先应该明确的是，利用abaqus自带的对python的兼容性，可以十分方便快捷的完成大批量的不同工况计算，比如这次我所计算的工况有480种，之前和师兄花了近一个月时间才手动计算完毕。而这次我利用python，加上学习python程序以及编写程序，总共花费近5天解决。极大的提高了效率。
我所用到的批量处理计算包括两部分，一个是批量的提交inp文件计算，一个是批量对计算结果odb文件提取数据。

批量提交inp文件计算

首先，我批量提交inp文件的前提是，对之前已经生成的480种工况下的inp文件批量修改，比如这次的初速度发生了改变，利用python程序对inp文件中的载荷这部分内容修改完，得到新的inp文件，我们需要利用abaqus的批处理来帮助我们提交计算。
我所知有两种批处理方式：

利用bat文件
编写.py程序，利用abaqus的command窗口运行
这次我们使用的是第一种方法：新建一个bat文件，并写入对应的批处理命令：

import os
workpath=[]
for root,dirs,files in os.walk('E:/ABAQUS_work/First-impact'):if len(root)>50:workpath.append(root)
with open('finalcompute.bat','w+') as sqlfile:sqlfile.write('@echo off\n')for i in range(0,len(workpath)):sqlfile.write('pushd '+workpath[i]+'\n')sqlfile.write('call abaqus job=PA01add'+' cpus=10 int'+'\n')sqlfile.write('pause')

因为我们不同工况的inp文件放在不同文件夹里，首先需要遍历所有inp文件所在的工作路径，这也是第一个片段的作用：建立包含所有inp文件路径的列表。第二个片段是将批命令代码写入到bat文件中。其中，最开始写入的“@echo off\n”是标准开头；写入的“’pushd ‘+workpath[i]+’\n‘”作用是设定工作路径，保证生成的结果文件都是在对应inp文件所在文件夹中；写入的“‘call abaqus job=PA01add’+’ cpus=10 int’+’\n’”是调用abaqus软件计算文件名为“PA01add.inp”的文件，分配10核cpu进行计算，'int’是表示一个文件计算完之后再进行下一个的计算；最后写入的“pause”作用是让计算完成之后的cmd窗口不会自动关闭，而是保持在最后完成的状态，方便我们查看在批量计算过程中是否出现问题。

批量提取数据

上一步我们已经对480种工况的inp文件进行了计算，生成了相应的odb结果文件，现在从odb文件中提取某一个点的载荷-时间曲线，所使用的方法为此前提到过的第二种方法，运行.py文件，文件内代码如下：


import os
resultpath=[]
for root,dirs,files in os.walk('E:/ABAQUS_work/First-impact'):if len(root)>50:print(root)resultpath.append(root)
from abaqus import *
from abaqusConstants import *
from viewerModules import *
from driverUtils import executeOnCaeStartup
executeOnCaeStartup()
resultname=['PA01add.odb','PA02add.odb','PA03add.odb']
houzhui=['PA01.txt','PA02.txt','PA03.txt']
for m in range(0,80):for i in range(0,3):o2 = session.openOdb(name=resultname[i],path=resultpath[m]+'/'+resultname[i])session.viewports['Viewport: 1'].setValues(displayedObject=o2)odb = session.odbs[resultpath[m]+'/'+resultname[i]]xyList = xyPlot.xyDataListFromField(odb=odb, outputPosition=NODAL, variable=(('RF', NODAL, ((COMPONENT, 'RF2'), )), ), nodePick=(('ASSEMBLY', 1, ('[#1 ]', )), ), )name=resultpath[m].split('\\')name2=name[2].split('-')writefile=name2[0]+'-'+name2[1]+'-'+houzhui[i]session.writeXYReport(writefile,xyList)

其功能是在生成的odb结果文件中，读取节点RF的支反力RF2，并写入到文件中。值得注意的是，该代码的.py文件必须要在abaqus的command窗口读取，因为最开始调用的库都是abaqus内的库，如果直接用sublime之类的编辑器无法调用。
在abaqus的command窗口中，输入的内容为：

abaqus cae nogui=getdata1.py

getdata1.py为之前代码的文件。
要注意，上述代码，仅仅是将所有数据写入到了我们定义的txt文件中，且存放的位置无法调整，只能存放在getdata1.py所在的文件夹里面同时还有一些格式需要调整，还要在其中提取我们所需要的有效数据，于是我们用以下代码来实现：

import os
import openpyxl
writefilepath=[]
for root,dirs,files in os.walk('E:/C919_fuselage_stress/一次触地F-T曲线(new)/'):if len(root)>60:writefilepath.append(root)
resultpath=[]
for root,dirs,files in os.walk('E:/ABAQUS_work/First-impact'):if len(root)>50:#print(root)resultpath.append(root)
houzhui=['PA01.txt','PA02.txt','PA03.txt']for m in range(0,80):for i in range(0,3):name=resultpath[m].split('\\')name2=name[2].split('-')writefile=name2[0]+'-'+name2[1]+'-'+houzhui[i]excelfile=name2[0]+'-'+name2[1]+'-PA0'+str(i+1)+'.xlsx'with open(writefile,'r')as readfile:book=openpyxl.Workbook()sheet=book.create_sheet('sheet1',0)block=readfile.readlines()n=1label=[]for x in range(0,len(block)):if '             0.         \n' in block[x]:label.append(x)if len(label)>1:for y in range(label[0],label[1]+1):a=list(block[y])num1=[]num2=[]for z in range(0,23):if a[z]!=' ':num1.append(a[z])for z in range(23,46):if a[z]!=' ':num2.append(a[z])finalnum1=''.join(num1)finalnum2=''.join(num2)sheet.cell(n,1,float(finalnum1))sheet.cell(n,2,float(finalnum2))n=n+1if len(label)==1:for y in range(label[0],len(block)-4):a=list(block[y])num1=[]num2=[]for z in range(0,24):
#                       print(z)if a[z]!=' ':num1.append(a[z])for z in range(24,46):if a[z]!=' ':num2.append(a[z])finalnum1=''.join(num1)finalnum2=''.join(num2)sheet.cell(n,1,float(finalnum1))sheet.cell(n,2,float(finalnum2))n=n+1
#           print(writefilepath[m]+'/'+excelfile)book.save(writefilepath[m]+'/'+excelfile)

以上代码是利用第三方库openpyxl，创建一个excel，并把之前生成的数据文件中的有效数据内容，按对应的时间、载荷依次放入到Excel中。提取有效数据的方式是判断0出现的行数，仅仅提取两次出现中间的数据，然后将一行中的数据，先变成列表，然后将不为空格的部分组装成字符串，最后再将字符串转变为浮点数。（此部分的判断依据，是根据我自己的结果文件进行的判断，不同的结果判断标准各不一样）

以上内容，能够实现abaqus的批量提交计算与批量提取数据，代码很简陋很臃肿，仅作为工作经验留存，如有不足还请见谅。

参考内容：https://wenku.baidu.com/view/ee78ab75f242336c1eb95e52.html