Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法

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

1目的

在有限元計(jì)算中,可能會(huì)碰到材料參數(shù)(這里以線彈性材料彈性模量E為例)在模型中并不為常值的情況,例如,E與某一坐標(biāo)相關(guān)。大家知道,abaqus材料中材料是在*section對(duì)截面屬性定義的時(shí)候?qū)卧M進(jìn)行賦值的。并且GUI中也沒(méi)有提供變參數(shù)的接口(這里順路提一下,載荷是可以通過(guò)函數(shù)給出變載荷的,GUI中有接口)。這就需要我們采取一些特殊的措施。本教程給出2種行之有效的方法,它們分別是:(1)編輯inp文件和(2)使用Umat。

2方法

方法	要求	難度
編輯inp	會(huì)一門(mén)編程語(yǔ)言,能夠完成對(duì)文本文件的讀取和寫(xiě)入	4星
使用Umat	會(huì)使用abaqus的Umat用戶子程序,計(jì)算文件中不需要使用其它的用戶子程序	3星(如果對(duì)abaqus用戶子程序不熟悉,則為5星)

上表中給出兩種方法的一些基本要求,編輯inp文件的方法要求更低,它通過(guò)直接編輯inp文件的方法來(lái)實(shí)現(xiàn)。而使用Umat的方法要求更高,需要用戶對(duì)子程序有一定的了解,并且也需要會(huì)一門(mén)編程語(yǔ)言。

同樣是根據(jù)坐標(biāo)確定材料參數(shù),使用Umat時(shí),可以直接通過(guò)abaqus用戶子程序自帶的子函數(shù)直接獲得單元的節(jié)點(diǎn)坐標(biāo),進(jìn)而通過(guò)坐標(biāo)確定材料參數(shù);而使用編輯inp的方法則需要首先通過(guò)編寫(xiě)程序?qū)δ阋呀?jīng)生成的inp文件(包含節(jié)點(diǎn)、單元信息)進(jìn)行讀取,之后對(duì)每一個(gè)單元進(jìn)行材料賦值。從這方面講,對(duì)熟悉Umat的用戶而言,使用Umat解決這個(gè)問(wèn)題更為方便。

3編輯inp方法的簡(jiǎn)要實(shí)施過(guò)程

以Fortran+VS為例,在編輯inp文件時(shí),涉及到了單元節(jié)點(diǎn)信息的讀取、計(jì)算單元(質(zhì)心)坐標(biāo)、計(jì)算彈性模量E、和輸出這幾個(gè)子程序。如圖1所示,prereadnew為讀取節(jié)點(diǎn)、單元信息子程序;processmat為通過(guò)單元坐標(biāo)計(jì)算彈性模量子程序;writeabaqus為輸出子程序。各個(gè)子程序的詳細(xì)信息這里不再贅述,有興趣或者有需求的同學(xué)可以關(guān)注文章最后的公眾號(hào)或者和作者私下聯(lián)系。

Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法abaqus有限元資料圖片1

圖1 項(xiàng)目及子程序

Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法abaqus有限元資料圖片2

圖2 inp文件讀取子程序(部分代碼)

Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法abaqus有限元分析案例圖片3

圖3 inp文件輸出-材料部分(部分代碼)

通過(guò)以上敘述不難發(fā)現(xiàn),通過(guò)編輯inp文件可以較好的解決材料參數(shù)不為常數(shù)的情況。但是,這需要用戶具有一定的編程基礎(chǔ),對(duì)abaqus的inp計(jì)算文件有一定的了解。同理,在其它軟件進(jìn)行有限元計(jì)算時(shí),也可以通過(guò)編輯計(jì)算文件的方法來(lái)解決材料參數(shù)不為常值得情況。

4 Umat方法的簡(jiǎn)要實(shí)施過(guò)程

作為Abaqus提供的用戶自定義子程序中的一種,Umat為用戶自己定義一種新的材料提供了接口。UMAT子程序具有強(qiáng)大的功能,使用UMAT子程序:

(1)可以定義材料的本構(gòu)關(guān)系,對(duì)ABAQUS 材料庫(kù)進(jìn)行擴(kuò)充;

(2)幾乎可以將任何本構(gòu)關(guān)系運(yùn)用到ABAQUS中的任何單元中;

但是,用戶必須在UMAT中提供材料本構(gòu)模型的雅可比(Jacobian)矩陣,即應(yīng)力增量對(duì)應(yīng)變?cè)隽康淖兓?而實(shí)際上,這也是Umat子程序最主要的任務(wù)。

承前所述,由于假定彈性模量與坐標(biāo)相關(guān),使用Umat子程序的首要任務(wù)就是找到彈性模量E和坐標(biāo)之間的關(guān)系;之后,才是根據(jù)坐標(biāo)得到彈性模量和泊松比之后對(duì)DDEDDS的編寫(xiě)。這里,依據(jù)Abaqus的幫助文檔,給出了Umat子程序中一些變量的定義:

Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法abaqus有限元分析案例圖片4

圖4 Abaqus提供的Umat函數(shù)接口

Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法abaqus有限元分析案例圖片5

圖5 部分Umat變量定義

詳細(xì)的參數(shù)介紹大家可以查看abaqus幫助。不難發(fā)現(xiàn)其中并不包含有單元或者節(jié)點(diǎn)的坐標(biāo)信息,但是卻包含有單元的單元編號(hào)信息(NOEL)。而通過(guò)單元編號(hào)又可以關(guān)聯(lián)起來(lái)單元的坐標(biāo)。

Abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法abaqus有限元分析案例圖片6

圖6 線彈性材料的DDEDDS

圖6中EMOD=PROPS(1)*A(NOEL)/0.125 就是根據(jù)單元NOEL的坐標(biāo)確定其彈性模量的語(yǔ)句,不難看出這個(gè)彈性模量是INP文件中給定的值得A(NOEL)/0.125倍,其中A數(shù)組并不是自定義的數(shù)組,為單元的坐標(biāo)值數(shù)組。

5 算例

對(duì)于如圖8所示的4階魔方結(jié)構(gòu),x軸左端面約束x方向自由度,x軸右端面加載位移。材料彈性模量由約束面到加載面分別為10、30、50、70GPa,而其對(duì)應(yīng)的單元質(zhì)心x軸坐標(biāo)分別為0.125、0.375、0.625、0.875。圖8中所示云圖為單元應(yīng)變E11。