[计算机]abaqusPython实例-操作excel文件.doc

《[计算机]abaqusPython实例-操作excel文件.doc》由会员分享，可在线阅读，更多相关《[计算机]abaqusPython实例-操作excel文件.doc（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、计算机abaqusPython实例-操作excel文件abaqus Python实例-操作excel文件目前处理数据离不开excel，所以pythoner必须学会用python操作excel表格.Python与excel交互方法也比较多，我一开始就接触的xlrd/xlwt包。直到现在也没有发现什么bug或者缺点，而且上次从ujs505那里知道Win32Com的方法在64bit系统下无效了,但是xlrd/xlwt方法依然好用，我欣喜万分.这里给出我一个从材料单轴拉伸数据TRA文件中处理得到拉伸过程中各种应力应变量并记录到excel中的例子,其实这个例子比较乱，但是实在没有时间给弄个明了的版本了。

2、大家将就看看from math import *import sysimport reimport xlrdfrom xlwt import Workbookbook=Workbook（)sheet=book。add_sheet(test_NT,cell_overwrite_ok=True）sheet1=book.add_sheet（engineer，cell_overwrite_ok=True)sheet2=book。add_sheet(True,cell_overwrite_ok=True)sheet3=book。add_sheet（Plastic,cell_overwrite_ok=Tr

3、ue)sheet4=book。add_sheet(Plastic_modify,cell_overwrite_ok=True）total_nubmer=46nclown=0nclown_e=0nclown_t=0nclown_p=0nclown_pm=0for i in range(total_nubmer）:myfile=tensile-+str（i+1)f=open（myfile+。TRA,r）engineer_strain=true_strain=engineer_stress=true_stress=plastic_strain=test_force=s1=f。readline（）s2

4、=f.readline（）s3=f。readline（)s4=f.readline（）s1=re.sub（|t;|n, ,s1)s1=re。split( +，s1）name=s12little=name-2：s2=re.sub(”t|；n, ，s2）s2=re。split（ +,s2)times=s22print name+run+str（times）bk=xlrd。open_workbook(Static report。xls)sh=bk。sheet_by_name(”Stat”)diameter=sh。cell_value(i+2，4)area=pi（diameter)*2）/4E_mod

5、ulus=sh.cell_value(i+2,7）state=Trues0=f.readlines（）s0=s02s0=re.sub（;, ，s0)s0=re.split( +,s0)f。seek(0)s1=f。readline(）s2=f。readline（）s3=f。readline（）s4=f.readline(）while state:s=f。readline（）s=re.sub(;， ，s）s=re.split( +，s)print s#print s0if （s=”EOF)（s=s0)：state=Falseelse：force=float（s1）eee=float(s0)/100

6、.0sss=force/areatest_force。append（force）engineer_strain.append(eee)engineer_stress。append(sss）true_strain。append（log(1.0+eee）)true_stress.append（sss*（1。0+eee)plastic_strain。append（log(1.0+eee)-sss（1。0+eee)/E_modulus)f.close(）#3simu_strain=simu_P_strain=simu_triax=simu_force=simu_E=204323。0simu_littl

7、e=str(int（little)+3)simu_name=”New_pass+simu_little+test.dat”print refer to+simu_namesimu_f=open(simu_name，r）simu_state=Trues0=simu_f.readlines()s0=s01s0=re.sub（；， ,s0）s0=re。split（ +，s0）print s0simu_f。seek(0）simu_s1=simu_f。readline()while simu_state：simu_s2=simu_f。readline（）simu_s2=re.sub（;， ，simu_s

8、2)simu_s2=re.split（ +,simu_s2)if （simu_s2=EOF”)|(simu_s2=s0):simu_state=Falseelse:print simu_s2s_force=float(simu_s25)s_eee=float(simu_s22）s_triax=float(simu_s23)s_sss=s_force/areasimu_force.append(s_force）simu_triax。append(s_triax）simu_strain。append(s_eee）simu_P_strain。append（log（1。0+s_eee）-s_sss（1

9、.0+s_eee)/simu_E）simu_f。close（）#3nrows=3sheet。write（0，nclown,sample)sheet.write（0,nclown+2,name)sheet.write(1，nclown，”area”）sheet.write（1，nclown+2，area)sheet.write（1,nclown+1,diameter)sheet。write（1，nclown+3,”E_modulus)sheet.write（1,nclown+4，E_modulus）sheet.write(2，nclown,E_strain）sheet.write(2，nclow

10、n+1,E_stress”)sheet.write(2，nclown+2，”T_strain）sheet。write（2，nclown+3，T_stress）sheet.write(2,nclown+4,PL_strain)for i in range（len（engineer_strain）：sheet.write(nrows,nclown,engineer_straini）sheet。write(nrows,nclown+1,engineer_stressi）sheet。write(nrows,nclown+2,true_straini)sheet.write（nrows，nclown+3

11、，true_stressi)sheet。write（nrows,nclown+4，plastic_straini）nrows+=1nclown+=5output engineer strain and stressnrows=3sheet1.write(0,nclown_e，file)sheet1。write(0，nclown_e+1，myfile）sheet1.write(0,nclown_e+2，E_modulus”）sheet1.write(0，nclown_e+3,E_modulus）sheet1.write(1，nclown_e，sample)sheet1。write（1,nclow

12、n_e+1,name）sheet1。write（1,nclown_e+2，run）sheet1。write(1,nclown_e+3,times）sheet1。write（2，nclown_e,E_strain”）sheet1.write(2，nclown_e+1，E_stress”）for i in range（len（engineer_strain）:sheet1.write(nrows,nclown_e，engineer_straini)sheet1.write(nrows,nclown_e+1,engineer_stressi)nrows+=1nclown_e+=4#output tr

13、ue strain and stressnrows=3sheet2.write（0，nclown_t,file)sheet2。write（0,nclown_t+1,myfile）sheet2。write（0，nclown_t+2，E_modulus)sheet2。write（0,nclown_t+3，E_modulus）sheet2。write(1,nclown_t，sample)sheet2。write(1，nclown_t+1，name）sheet2。write(1,nclown_t+2,run)sheet2。write（1，nclown_t+3,times)sheet2.write(2,

14、nclown_t,”T_strain）sheet2。write(2,nclown_t+1，”T_stress”）for i in range(len（engineer_strain）)：sheet2.write(nrows，nclown_t，true_straini）sheet2。write（nrows,nclown_t+1，true_stressi)nrows+=1nclown_t+=4output plastic strain and stressnrows=3sheet3.write（0，nclown_p，file）sheet3.write（0，nclown_p+1，myfile)she

15、et3.write（0,nclown_p+2,E_modulus)sheet3。write（0,nclown_p+3,E_modulus）sheet3。write（1,nclown_p,sample)sheet3.write（1,nclown_p+1，name)sheet3。write（1,nclown_p+2，run）sheet3。write(1,nclown_p+3，times)sheet3.write(2,nclown_p,”PL_strain）sheet3.write（2,nclown_p+1，”T_stress”）for i in range（len(engineer_strain）

16、)：sheet3。write（nrows，nclown_p,plastic_straini)sheet3.write（nrows,nclown_p+1,true_stressi)nrows+=1nclown_p+=4output modified plastic strain and true stressnrows=3i_temp=1sheet4.write(0,nclown_pm,file）sheet4。write(0，nclown_pm+1,myfile）sheet4。write(0,nclown_pm+2，E_modulus）sheet4。write(0,nclown_pm+3,E_m

17、odulus)sheet4.write(1,nclown_pm，sample）sheet4.write(1，nclown_pm+1,name）sheet4。write（1，nclown_pm+2，run）sheet4.write（1,nclown_pm+3，times）sheet4。write（2,nclown_pm,”PL_strain）sheet4.write(2，nclown_pm+1，”Force”）running=Truefor i in range（len(engineer_strain)）:if running:mincr=150temp_e=（test_forceitest_f

18、orcei+mincr）/(plastic_strainiplastic_straini+mincr)if temp_e400000.0:i_temp=irunning=Falseelse:modifed_strain=engineer_straini-engineer_straini_tempmodifed_plastic_strain=log（1。0+modifed_strain)-engineer_stressi*(1。0+modifed_strain）/E_modulussheet4。write(nrows，nclown_pm，modifed_plastic_strain)sheet4

19、.write（nrows,nclown_pm+1,test_forcei)nrows+=1output simulation resultnrows=3sheet4.write(2，nclown_pm+2，”PL_strain）sheet4。write（2,nclown_pm+3,”Force）sheet4.write（2,nclown_pm+4,”TRIAX)for j in range(len(simu_P_strain)）:sheet4。write(nrows，nclown_pm+2，simu_P_strainj)sheet4.write（nrows,nclown_pm+3，simu_f

20、orcej）sheet4。write（nrows，nclown_pm+4，simu_triaxj)nrows+=1#next file(Test result)nclown_pm+=5book。save(material_treat。xls）另一例子import csvfrom odbAccess import from abaqusConstants import filename=getInput（Please input the ODB file name)下面这样都是我定义的字典或list，用来存保存提取数据的elementArea= elementConn=nodeArea =tim

21、eTP =#下面可以看做一般的从odb文件中提取结果的步骤打开指定的odb文件odb=openOdb（path=filename）得到assembly中所有的instanceinst=odb.rootAssembly。instances #或取第一个instance中的所有单元,因为我这个odb里面只有一个instance#.keys（）方法可以获得一个list有所有的instance名字elments=instinst。keys（)0。elements下面一段对单元循环，得到每个单元的几个结点,然后记录下来for el in elments： label=el。label nodes=el。

22、connectivity elementConnlabel=nodes#节点nodes=instinst。keys()0。nodesfor nd in nodes： label=nd。label nodeArealabel=0打开指定的stepst=odb。stepsodb.steps.keys(）0对指定step的某个特定场变量做循环for v in st。frames1。fieldOutputsEVOL.values: label=v.elementLabel data =v.data elementArealabel=datafor k，v in elementArea.iteritem

23、s（)： nds=elementConnk for nd in nds: nodeAreand+=0。25v 下面是提取每个frame中每个节点的NT11值frames=st。frames for fr in frames： nt11=fr.fieldOutputsNT11.values time=fr.frameValue sumTp=0 sumAr=0 for va in nt11: label=va.nodeLabel data=va。data if data0： area=nodeArealabel sumTp+=data*area sumAr+=area try： mean=float(sumTp)/float(sumAr) timeTP。append( （time，mean)） except: print fr.frameIdfilename=getInput（pleae input the csv file n name to save)wr=csv.writer(file(filename ,wb)#这个把结果写进csv文件，只要一句,很简单，也可以在这一句之前写个表头wr.writerows(timeTP）