ABAQUS中Model change(生死單元)的循環(huán)操作

2017-03-18 by:CAE仿真在線來源:互聯(lián)網(wǎng)

在ABAQUS進(jìn)行瞬態(tài)分析時,需要對每一步運(yùn)算結(jié)果進(jìn)行監(jiān)測,當(dāng)有單元應(yīng)力(或者其他條件)超過設(shè)定應(yīng)力時,需要將單元?dú)⑺?剩余單元進(jìn)行下一步計算分析,直到設(shè)定的終止?fàn)顟B(tài)。

第一步,在ABAQUS中使用Model change將符合條件的單元?dú)⑺赖牟僮魅鐖D1所示

ABAQUS中Model change(生死單元)的循環(huán)操作abaqus有限元培訓(xùn)資料圖片1

ABAQUS中Model change(生死單元)的循環(huán)操作abaqus有限元仿真圖片2

第二步,需要把滿足條件的單元選出,然后標(biāo)記,并且進(jìn)行循環(huán)設(shè)定。主要采用Python腳本語言進(jìn)行設(shè)置,現(xiàn)在采用《Python_Scripts_for_Abaqus.pdf》中10章和13章的內(nèi)容作為基礎(chǔ),懸臂板odb文件輸入應(yīng)力和位移的例子進(jìn)行進(jìn)一步分析,板一端固定一端受到兩個集中力作用,現(xiàn)在不斷調(diào)整集中力大小,分析板中應(yīng)力,當(dāng)板中某些單元的應(yīng)力超過某個選定的值時,該板單元設(shè)置為“死亡”,其余板繼續(xù)分析。

核心步驟為:

(1)建立基本懸臂板模型;

(2)建立分析步,施加外荷載,建立Job文件進(jìn)行計算,提取對應(yīng)ODB文件中的每個單元應(yīng)力;

(3)把上一步驟中應(yīng)力超過設(shè)定值的單元標(biāo)記,設(shè)定失效單元集合;

(4)建立一個新的分析步,對上一步驟的中的失效單元集合進(jìn)行Model change,然后生成一個新的Job文件重新計算,提取新的ODB文件中每個有效單元的應(yīng)力;

(5)重復(fù)(3)(4),直至中止。

ABAQUS中Model change(生死單元)的循環(huán)操作abaqus有限元仿真圖片3

ABAQUS中Model change(生死單元)的循環(huán)操作abaqus有限元仿真圖片4

ABAQUS中Model change(生死單元)的循環(huán)操作abaqus有限元仿真圖片5

其中腳本文件如下(詳細(xì)文件下載地址為:http://forum.simwe.com/thread-1241669-1-1.html)

核心代碼如下(實際代碼中無編號):

(1)elemSet =[platePart.elements[int(i)-1:int(i)] for i in highlight_element_list]

(2)platePart.Set(elements=elemSet,name='Set-elements'+repr(k))

(3)region_elementset=plateInstance.sets['Set-elements'+repr(k)]

(4)plateModel.ModelChange(name='Invalid'+repr(k+1),createStepName='Load Step'+repr(k+1),region=region_elementset,regionType=ELEMENTS, activeInStep=False, includeStrain=False)

開放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才