1、引言
傳統(tǒng)的汽車通常是全鋼車身,所用鋼材延展性較好,碰撞后車體結構主要發(fā)生潰縮和彎折變形,只有極少量結構件可能斷裂失效。所以在進行碰撞仿真時只要求材料卡片能夠準確描述材料彈塑性行為,一般無需對材料失效進行數(shù)值模擬,少量部件失效的影響可以基于工程經(jīng)驗進行估計。
近年來,輕量化材料,如鎂鋁合金、高強熱成型鋼和非金屬復合材料在汽車行業(yè)得到了廣泛的應用,這些材料的延展性遠低于普通鋼材,在車輛發(fā)生碰撞時極有可能發(fā)生斷裂,如圖1。所以有必要考慮其失效行為并在仿真分析中應用合適的材料失效模型。
圖1 鋁合金車體碰撞時失效開裂
2、應力三軸度與失效應變
在碰撞仿真中,我們可以基于單向拉伸試驗的失效塑性應變來判斷材料是否失效,這就是所謂的常應變失效準則。常應變失效準則按照所采用應變指標的不同,又可細分為多種,其中最常見的有最大主應變準則和等效塑性應變準則,分別如公式(1)和公式(2)所示。
其中εmax為單拉工況的失效應變,εpmax為單拉工況的失效塑性應變,對于大多數(shù)材料,發(fā)生失效時塑性應變遠大于彈性應變,所以我們可以認為失效應變和失效塑性應變近似相等。εp1、εp2和εp3是三個塑性主應變。εp是等效塑性應變,也叫Von Mises塑性應變,其表達式與Von Mises應力類似,但前面多了一個系數(shù)2/3。
此處需要強調一點,通常認為塑性應變不會導致體積變化(即塑性應變的泊松比為0.5),所以三個塑性主應變之和恒定為0,
常應變失效準則最簡單最易實現(xiàn),但它的表達形式中沒有考慮到材料失效中的諸多因素,因而計算結果誤差較大。
對現(xiàn)有金屬材料的研究發(fā)現(xiàn),應力狀態(tài)對于失效等效應變數(shù)值起著決定性作用,材料所受應力狀態(tài)不同時,材料內(nèi)產(chǎn)生的塑性變形與應力集中程度將不同,材料失效應變數(shù)值也將發(fā)生變化。
在大多數(shù)的失效準則中,結構一點處的應力狀態(tài)通常采用應力三軸度η來表示,其表達式為
從應力三軸度的定義式可以看出,應力三軸度為靜水應力與Von Mises應力的比值。靜水應力導致體積變化,Von Mises應力反映了形狀改變。應力三軸度作為描述應力狀態(tài)的參數(shù),同時反映了體積變化和形狀改變。
需要注意的是,有些文獻中定義應力三軸度時,在公式(6)表達式前加了一個負號,材料受拉時應力三軸度為負,受壓時為正,與通常的定義方式恰好相反。
圖2展示了幾種典型加載工況下的應力三軸度數(shù)值,都是比例加載情況,應力三軸度在整個加載過程中保持恒定。
圖2 不同加載工況下的應力三軸度
將材料試件設計為不同的缺口形式,進行單雙軸拉伸以及剪切等形式的加載,以對應不同的應力三軸度,通過對其失效應變進行測量,就可以建立起應力三軸度和材料失效時的等效塑性應變之間的對應關系,如圖3所示。從圖中可見,應力三軸度和等效塑性應變之間并不是簡單的單調函數(shù)關系,而是一條形狀比較復雜的曲線。
圖3 試驗測定應力三軸度和失效等效塑性應變關系曲線
3、為什么需要復雜的失效模型?
得到應力三軸度η與失效等效塑性應變εf的關系曲線后,似乎在碰撞仿真中就很容易判斷材料是否失效了。我們發(fā)現(xiàn)某個單元的等效塑性應變數(shù)值達到了當前應力三軸度所對應的失效塑性應變,就認為該單元失效,在下一時間步將該單元刪除即可。這樣看來我們并不需要再去研究復雜的失效模型,只要用常應變失效準則同時考慮應力三軸度η即可。
可惜的是,上述方案只適用于比例加載(η為恒定值)的情況,對于非比例加載工況并不適用。對于非比例加載情況,失效等效塑性應變不僅取決于加載結束時的η值,加載過程中的η值也有影響。
圖4展示了某種材料的η-εf曲線和4條加載路徑。兩條灰色加載路徑的應力三軸度η是恒定的,所以失效時的(η,εf)必然位于曲線上,我們根據(jù)η值就可以確定加載到何時會發(fā)生斷裂失效。而紅色和綠色路徑是非比例加載,應力三軸度η在加載過程中不斷變化,導致這兩條加載路徑的失效點不在η-εf曲線上,我們無法直接判斷出加載到何時發(fā)生失效。
圖4 比例加載和非比例加載示例
為了應對非比例加載情況,我們需要一種合適的失效模型,能夠考慮整個加載過程中材料的損傷演化,進而判斷出失效點。
LS-Dyna軟件中可以采用GISSMO模型、MMC模型和Johnson-Cook模型等多種材料失效模型,它們都引入了應力三軸度作為應力狀態(tài)參數(shù),都考慮加載過程中的損傷累積效應。其中GISSMO模型的卡片輸入比較簡潔,所需參數(shù)便于試驗測定,計算結果也比較符合實際,因此獲得了廣泛的工程應用。
4、GISSMO模型對損傷的累積計算
金屬材料的韌性斷裂是指金屬材料經(jīng)過塑性變形后發(fā)生的宏觀斷裂。一般金屬材料都是由多種成分組成,包括第二相粒子及夾雜物。金屬材料受力后,隨著宏觀變形的發(fā)展,金屬材料中局部產(chǎn)生微孔洞。而第二相粒子和夾雜物也會在材料中成核。微孔洞成核之后在變形過程中不斷的生長,與周圍的孔洞相互聚合最終形成宏觀的裂紋。
既然金屬韌性斷裂的微觀機理表現(xiàn)為微孔洞的成核、長大和聚合等演變過程,我們就可以利用內(nèi)部損傷變量來表征細觀結構缺陷引起的材料性能劣化的程度,以損傷變量的演化函數(shù)來描述材料從變形到破壞的整個過程。
GISSMO模型以及其他考慮損傷積累的模型,都是定義了一個損傷變量D∈{0,1},在加載過程中,材料的等效塑性應變εp逐步增加,材料的損傷D也在不斷累積,當損傷值D=1時,材料發(fā)生失效。
GISSMO模型定義了損傷值變化率和塑性應變率之間的關聯(lián),
其中εf為斷裂等效塑性應變,是當前的應力三軸度η的函數(shù),其數(shù)值根據(jù)試驗εf -η曲線確定;n為材料的損傷累積指數(shù)。
在有限元顯式分析中,實際使用的是公式(7)的增量形式,
每一個時間步計算一次損傷增量ΔD,當某個時間步的損傷值D≥1 ,即認為材料發(fā)生失效。
對于比例加載情況,應力三軸度η在加載過程中保持不變,所以εf也保持為常量。這樣公式(1)可以積分為更簡潔的形式,
顯然,對于比例加載情況,當εp=εf時D=1,此時材料發(fā)生失效。
5、GISSMO模型對材料本構曲線的修正
上節(jié)提到金屬韌性斷裂的過程實質是微孔洞的成核、長大和聚合的過程,這種損傷累積會對材料的本構關系產(chǎn)生影響。如圖5所示,材料的截面積為S,材料中的微孔洞導致材料的有效承載面積?