Python在ABAQUS二次開發(fā)中的應(yīng)用實例2課件

Python在ABAQUS二次開發(fā)中的應(yīng)用實例Python在ABAQUS二次開發(fā)中主要內(nèi)容一、ABAQUS腳本概述二、Python語言簡介三、ABAQUS腳本編寫主要內(nèi)容一、ABAQUS腳本概述一、ABAQUS腳本概述一、ABAQUS腳本概述ABAQUS腳本概述ABAQUS軟件簡介ABAQUS是國際著名的CAE軟件，它以解決實際工業(yè)問題能力和強大的非線性功能贏得廣泛聲譽。航空工業(yè)是ABAQUS最重要的應(yīng)用領(lǐng)域之一，波音、空中客車、洛克希德?馬丁等是其長期合作的用戶。ABAQUS分析過程由三個步驟組成：前處理（ABAQUS/CAE）；模擬計算(ABAQUS/Standard

或者ABAQUS/Explicit）；后處理（ABAQUS/CAE）。建模過程主要集中在前處理步驟中。ABAQUS腳本概述ABAQUS軟件簡介ABAQUS分析過ABAQUS腳本概述ABAQUS建模方法常用的ABAQUS建模方法主要有三種：①在ABAQUS/CAE界面建模優(yōu)點：界面直觀，簡單易懂；缺點：操作步驟繁瑣，不利于修改，缺乏通用性。②通過inp文件建模優(yōu)點：可以直接在ABAQUS中提交；缺點：只包含模型節(jié)點信息，代碼較多，修改不便；③使用Python語言建模優(yōu)點：代碼少，語言友好，可移

植性好，可進(jìn)行參數(shù)化修改；缺點：不夠直觀。ABAQUS腳本概述ABAQUS建模方法ABAQUS腳本概述ABAQUS腳本接口ABAQUS有限元分析軟件二次開發(fā)環(huán)境提供的腳本接口，是基于Python語言進(jìn)行的定制開發(fā)。在ABAQUS/CAE中進(jìn)行建模和進(jìn)行后處理時，對話框中做的所有設(shè)置都由ABAQUS/CAE從內(nèi)部發(fā)出與之對應(yīng)的命令（command）。ABAQUS腳本接口直接與內(nèi)核進(jìn)行通信（communicate），而與ABAQUS/CAE的圖形用戶界面（GUI）無關(guān)。如果將所有的腳本接口命令存儲于文件中，該文件則稱為腳本（script）。腳本由一系列純ASCII格式的Python語句組成，擴展名一般為.py。ABAQUS腳本概述ABAQUS腳本接口在ABAQUS/CABAQUS腳本概述ABAQUS腳本接口編寫腳本可以實現(xiàn)以下功能：自動執(zhí)行重復(fù)任務(wù)；進(jìn)行參數(shù)分析；創(chuàng)建和修改模型；訪問輸出數(shù)據(jù)庫（ODB文件）；定制ABAQUS環(huán)境文件；創(chuàng)建ABAQUS插件程序。編寫完成的腳本文件，可以直接在ABAQUS運行:File→RunScript，打開相應(yīng)的.py文件，即可執(zhí)行；剛打開ABAQUS界面，直接選擇RunScript也可；可以在GUI界面下部的命令行接口輸入腳本語句執(zhí)行。ABAQUS腳本概述ABAQUS腳本接口ABAQUS腳本概述創(chuàng)建ABAQUS腳本創(chuàng)建腳本有3種常用方法：①在GUI建模時錄制宏（macro）File→MacroManager可以按需錄制任何操作對應(yīng)的命令②直接編輯abaqus.rpy文件建模時工作目錄下自動生成abaqus.rpy文件使用軟件EditPlus可以直接打開編輯保存.cae文件并退出后才能生成.rpy文件③借助PythonReader.exe軟件中國石油大學(xué)焦中良開發(fā)可以實時獲取CAE每步操作對應(yīng)的腳本命令A(yù)BAQUS腳本概述創(chuàng)建ABAQUS腳本ABAQUS腳本概述調(diào)試ABAQUS腳本調(diào)試腳本有多種方法，這里介紹Abaqus自帶的PDE中的調(diào)試。PDE（PythonDevelopmentEnvironment，Python開發(fā)環(huán)境）有以下功能：編輯Python文件、調(diào)試Python腳本或插件；創(chuàng)建guiLog腳本，錄制所有的GUI動作；通過GUI（guiLog）運行Python腳本或在ABAQUS/CAE之外運行Python腳本；設(shè)置延遲（delay）和斷點（breakpoint），并查看代碼的執(zhí)行情況；在watch窗口中查看變量值等。打開方式File→AbaqusPDE，然后打開相應(yīng)的.py文件，即可編輯調(diào)試。ABAQUS腳本概述調(diào)試ABAQUS腳本ABAQUS腳本概述調(diào)試ABAQUS腳本點按鈕可以調(diào)試整個腳本文件；點

按鈕可以逐行調(diào)試腳本文件，調(diào)試行顯示為藍(lán)色；點

按鈕可以停止調(diào)試文件；點

按鈕可以設(shè)置斷點，調(diào)試時將調(diào)到此行為止；點

按鈕，然后點按鈕，可以錄制CAE中的操作，點可以停止錄制。ABAQUS腳本概述調(diào)試ABAQUS腳本二、Python語言簡介二、Python語言簡介Python語言簡介基本規(guī)則代碼按級別從大到小解讀，之間用操作符“.”連接。如：

mdb.models[‘Model-1’].materials[‘BPCOMPOSITE’].Depvar(n=5)

表示將模型‘model-1’的材料‘BPCOMPOSITE’里的子程序相關(guān)參數(shù)設(shè)為5。語句前加上#號，表示是注釋內(nèi)容，這個語句將不被執(zhí)行；使用中文注釋，要在文件最開始加上：#-*-coding:UTF-8-*-以換行來識別一個邏輯語句的結(jié)束。如果需要在一個物理行中編寫多個邏輯語句，需要用分號分開。如：x=1;y=2;z=3。如果一個邏輯語句太長需要分行寫，需要用反斜杠“\”進(jìn)行連接（有括號括起來的內(nèi)容就不用）。不同代碼塊之間要注意用冒號和縮進(jìn)來區(qū)分代碼塊之間的層次關(guān)系，如：

foriinrange(0,b3):

if(iAnglelowfg[i]=='/'):

forlinrange(0,i-j-1):

c3[k]=c3[k]+iAnglelowfg[j+l+1]

j=i

elif(iAnglelowfg[i]==']'):Python語言簡介基本規(guī)則Python語言簡介保留字段定義變量名、函數(shù)名的時候，避免使用保留字段。如：

保留字段說明and\not\or邏輯“與”“非”“或”break中止循環(huán)語句continue繼續(xù)執(zhí)行下一次循環(huán)class定義類def定義函數(shù)或者方法if\else\elif條件判斷語句print輸出語句while循環(huán)語句from\import導(dǎo)入模塊return返回函數(shù)計算結(jié)果Python語言簡介保留字段保留字段說明and\nPython語言簡介數(shù)據(jù)類型數(shù)據(jù)類型主要有數(shù)字類型、字符串類型；數(shù)字類型主要有四種：整型int、長整型long、浮點型float、復(fù)數(shù)型complex；Python語言可自動判斷數(shù)據(jù)類型，不需要編寫者定義說明類型，如：i=10.5，即可。不需要用：floati=10.5

；不同數(shù)據(jù)類型之間可通過函數(shù)轉(zhuǎn)換，如int（）、float（）?？梢杂锰娲柋硎径ㄖ谱址?，如

stress=1234.56

print‘Thestressis%fMPa’%stress

即可輸出：Thestressis1234.56Mpa

符號描述符號描述%c定制字符及其ASCII碼%o定制無符號八進(jìn)制數(shù)%s定制字符串%x定制無符號十六進(jìn)制數(shù)%d定制整數(shù)%f定制浮點型數(shù)據(jù)%u定制無符號整數(shù)%e用科學(xué)計數(shù)法定義浮點型數(shù)據(jù)Python語言簡介數(shù)據(jù)類型符號描述符號描述%c定制字符Python語言簡介運算符主要有賦值運算符、算術(shù)運算符、邏輯運算符、關(guān)系運算符；下表按優(yōu)先級排列，從上往下優(yōu)先級增加。

運算符描述or邏輯或and邏輯與not邏輯非<,<=,>,>=,!=,==小于，小于等于，大于，大于等于，不等于，等于+，-加、減運算*,/,%乘、除、取余數(shù)運算**指數(shù)運算不支持自加或自減運算，如：i++、i--；但類似i+=1這樣的語句是可以的。

Python語言簡介運算符運算符描述or邏輯或and邏輯與三、ABAQUS腳本編寫三、ABAQUS腳本編寫ABAQUS腳本編寫定義和導(dǎo)入模塊在每個腳本文件的最開始處，首先要導(dǎo)入相應(yīng)的ABAQUS模塊；每個步驟均有對應(yīng)的模塊，如part模塊，material模塊、assembly模塊等，通過以下語句可以直接導(dǎo)入cae模塊的所有對象:

fromcaeModulesimport*如果需要用到各種符號常數(shù)，如單元類型C3D8R，還要導(dǎo)入符號常數(shù)模塊：

fromabaqusConstantsimport*定義和導(dǎo)入對象還需要導(dǎo)入相應(yīng)的ABAQUS對象，包括Session對象（包括視口、遠(yuǎn)程隊列、視圖對象）Mdb對象（包括Model對象和Job對象）Odb對象通過類似如下的語句導(dǎo)入對象：

fromabaqusimport*ABAQUS腳本編寫定義和導(dǎo)入模塊ABAQUS腳本編寫創(chuàng)建新實體myModel=mdb.models['Model-1']myViewport=session.Viewport(name='Regionsyntax',origin=(20,20),width=200,height=100)mySketch=myModel.Sketch(name='SketchA',sheetSize=200.0)mySketch.rectangle(point1=(-40.0,30.0),point2=(-10.0,0.0))mySketch.rectangle(point1=(10.0,30.0),point2=(40.0,0.0))door=myModel.Part(name='Door',dimensionality=THREE_D,type=DEFORMABLE_BODY)door.BaseSolidExtrude(sketch=mySketch,depth=20.0)通過以上語句，在Part模塊中生成兩個三維實體。#建立新模型Model-1#創(chuàng)建新的視口#創(chuàng)建一個草圖，并繪制兩個矩形#對上面創(chuàng)建的兩個矩形進(jìn)行拉伸操作，生成三維部件

ABAQUS腳本編寫創(chuàng)建新實體#建立新模型Model-1ABAQUS腳本編寫裝配實體、導(dǎo)入邊界條件myAssembly=myModel.rootAssemblydoorInstance=myAssembly.Instance(name='Door-1',part=door,dependent=OFF)pillarVertices=doorInstance.vertices.findAt(((-40,30,0),),((40,0,0),))myModel.StaticStep(name='impact',previous='Initial',initialInc=1,timePeriod=1)myPillarLoad=myModel.ConcentratedForce(name='pillarForce',createStepName='impact',region=(pillarVertices,),cf1=12.50E4)通過以上語句，裝配實體，并施加邊界條件。此處使用了findAt語句來查找頂點。#導(dǎo)入實體，完成裝配#選擇兩個頂點#創(chuàng)建靜力分析步（static）#在選擇的頂點上施加集中力

ABAQUS腳本編寫裝配實體、導(dǎo)入邊界條件#導(dǎo)入實體，完成裝ABAQUS腳本編寫findAt語句簡介在Abaqus建模操作生成的rpy文件中，Abaqus將自動為模型的特征進(jìn)行編號，每個節(jié)點(node)、單元(element)、頂點(vertex)、邊(edge)、面(face)、體(cell)，都有對應(yīng)的ID。實例：對方形實體切割

p=mdb.models['Model-1'].parts['Door']

c=p.cells

pickedCells=c.getSequenceFromMask(mask=('[#1]',),)

v,e,d=p.vertices,p.edges,p.datums

p.PartitionCellByPlaneThreePoints(point1=v[0],

point2=v[2],point3=v[6],cells=pickedCells)如果后續(xù)對模型進(jìn)行編輯（增加、刪除），特征

ID會改變，這時再用原來的語句就報錯。解決方法：用該特征的空間坐標(biāo)來表示；對于點，坐標(biāo)是唯一的；對于線、面、體，上面的任一坐標(biāo)都可以，但不要用兩個特征的交點坐標(biāo)。p.PartitionCellByPlaneThreePoints(point1=(-10,30,0),point2=(-40,30,20),point3=(-40,0,20),cells=pickedCells)c.findAt(((-25,30,20),))ABAQUS腳本編寫findAt語句簡介p.PartitioABAQUS腳本編寫定義材料屬性mdb.models['Model-1'].Material(name='Steel')mdb.models['Model-1'].materials['Steel'].Elastic(table

=((220000.0,0.35),))mdb.models['Model-1'].HomogeneousSolidSection(name=

'Section-1',material='Steel',thickness=None)pickedCells=c[:]p.Set(cells=pickedCells,name='Set-all')p=mdb.models['Model-1'].parts['Door']c=p.cellscells=c[:]region=p.sets['Set-all']p=mdb.models['Model-1'].parts['Door']p.SectionAssignment(region=region,sectionName='Section-1',offset=0.0)#定義材料“steel”。#定義材料“steel”的工程彈性常數(shù)。#定義截面屬性“steel”，將材料“steel”賦予截面#c[:]可以選取全部的實體#定義集合“Set-all”，將全部實體定義成一個集合#將截面屬性賦予給集合“Set-all”

ABAQUS腳本編寫定義材料屬性#定義材料“steel”。ABAQUS腳本編寫布置種子a=mdb.models['Model-1'].rootAssemblye1=a.instances['Door-1'].edgespickedEdges=e1[:]#e1[:]可以選取全部的邊a.seedEdgeByNumber(edges=pickedEdges,number=4)通過以上語句布置種子。a=mdb.models['Model-1'].rootAssemblypartInstances=(a.instances['Door-1'],)a.generateMesh(regions=partInstances)通過以上語句，把實體劃分網(wǎng)格。提交工作mdb.Job(name='Job-1',model='Model-1',description='',type=ANALYSIS,atTime=None,waitMinutes=0,waitHours=0,queue=None,memory=90,memoryUnits=PERCENTAGE,getMemoryFromAnalysis=True,

explicitPrecision=SINGLE,nodalOutputPrecision=SINGLE,echoPrint=OFF,