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

2017-03-18 by:CAE仿真在線來源:互聯網

在ABAQUS進行瞬態分析時,需要對每一步運算結果進行監測,當有單元應力(或者其他條件)超過設定應力時,需要將單元殺死,剩余單元進行下一步計算分析,直到設定的終止狀態。

第一步,在ABAQUS中使用Model change將符合條件的單元殺死的操作如圖1所示

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

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

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

核心步驟為:

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

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

(3)把上一步驟中應力超過設定值的單元標記,設定失效單元集合;

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

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

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

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

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

其中腳本文件如下(詳細文件下載地址為: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)

開放分享：優質有限元技術文章,助你自學成才