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所示
第二步,需要把滿足條件的單元選出,然后標(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),直至中止。
其中腳本文件如下(詳細(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)
相關(guān)標(biāo)簽搜索:ABAQUS中Model change(生死單元)的循環(huán)操作 abaqus分析培訓(xùn) abaqus技術(shù)教程 abaqus巖土分析 鋼筋混凝土仿真 abaqus分析理論 abaqus軟件下載 abaqus umat用戶子程序編程 Abaqus代做 Abaqus基礎(chǔ)知識 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn)