碳纖維增強樹脂基斷裂的數(shù)值模擬研究

程黨根，朱路紅，文 韜

(長沙航空職業(yè)技術(shù)學(xué)院，湖南 長沙 410124)

碳纖維增強樹脂基復(fù)合材料已經(jīng)廣泛地應(yīng)用于航空、航天等各個領(lǐng)域。然而，目前碳纖維增強樹脂基復(fù)合材料作為結(jié)構(gòu)應(yīng)用的最大障礙之一是其斷裂韌性比金屬材料差很多。斷裂韌性是應(yīng)用線彈性斷裂力學(xué)在各向同性材料中建立的描述材料抵抗裂紋擴展能力的力學(xué)性能參數(shù)。復(fù)合材料是非均質(zhì)的材料，含有基體、纖維和界面三相 ，這幾種相結(jié)構(gòu)對于裂紋擴展的阻力不同;另外復(fù)合材料的斷裂過程往往伴隨著界面脫粘、纖維橋接、纖維撥出等能量消耗過程，這便造成復(fù)合材料內(nèi)部的裂紋不遵循自相似擴展的規(guī)律，即裂紋體在受到垂直于裂紋面的拉伸載荷時，裂紋沿著原裂紋面擴展的規(guī)律。一般認(rèn)為線彈性斷裂力學(xué)不能適用于復(fù)合材料，因此有必要建立適合纖維增強復(fù)合材料斷裂測試的方法［1］，并選擇合適參數(shù)來表征其斷裂韌性。對于金屬材料，一般選用KIC和JIC來表征其斷裂韌性，如彈性材料即裂紋尖端無塑性區(qū)或者塑性區(qū)很小的自相似擴展I型裂紋，用KIC來表征其斷裂韌性;延性材料即裂紋塑性區(qū)較大，并且滿足J主導(dǎo)區(qū)條件時，用JIC來表征其斷裂韌性。對于碳纖維增強樹脂基復(fù)合材料，由于斷裂機理的復(fù)雜性，其破壞模式隨著纖維、基體及界面性能的不同，可能是I型、II型甚至復(fù)合型，因此不能直接將金屬材料斷裂韌性的表征參數(shù)用于復(fù)合材料中。為此，雖然有許多科學(xué)工作者在這方面做了大量研究工作，但是復(fù)合材料斷裂韌性的測試和表征仍沒有完全解決，仍有很多問題還需要進一步的研究。

1 碳纖維增強樹脂基斷裂模型

根據(jù)試驗情況，在Abaqus平臺上建立碳纖維增強樹脂基的J積分有限元計算的裂紋模型。試件幾何尺寸為 60mm×10mm×5mm，跨距為40mm，均勻鋪層單層厚度為0.125mm，共40層，層合板的鋪層參數(shù)為［45°/0°/90°/－45°］s，復(fù)合材料單層板的材料參數(shù)采用表1的參數(shù)。最終在ABAQUS建立的復(fù)合材料有限元模型［2］如圖1所示，模型選用三維二次減縮殼單元(S8R5)，本文中的模擬對網(wǎng)格并不敏感，因此只在裂紋尖端處細化網(wǎng)格，并進行奇異性處理。

表1 試驗件的材料特性

E1、E2、E3分別是1，2，3 方向楊氏(壓縮、拉伸)彈性模量;G12G13G23分別是12、13、23方向切變彈性模量;V12材料12方向斷裂韌性參數(shù)。

圖1 有限元模型網(wǎng)格

1.1 裂紋尖端奇異性處理

在圍線積分區(qū)域內(nèi)，對裂紋尖端處應(yīng)用1/4奇異性單元來模擬裂紋尖端的奇異性能夠更真實地模擬裂紋尖端的奇異性，另外也能更準(zhǔn)確地計算圍線的J積分。奇異單元是通過調(diào)整高階單元中間節(jié)點的位置而來的，如把二階單元的中間節(jié)點放在1/4位置，如圖2所示，就會在裂紋尖端處產(chǎn)生平方根的奇異性［3］。

圖2 1/4節(jié)點奇異單元

在Abaqus網(wǎng)格劃分過程中，目的是為了得到在裂紋區(qū)域很密的單元，而在其他區(qū)域較疏松，同時在裂紋的周圍n(取決于自己定義圍線積分的計算層數(shù))層單元必需用六面體單元。通常的做法是用“蜘蛛網(wǎng)”形式來實現(xiàn)這種網(wǎng)格的疏密過渡，如圖3所示，特別在Abaqus里面網(wǎng)格的過渡過程參數(shù)設(shè)置不合適，造成過程區(qū)域網(wǎng)格的扭曲經(jīng)常會導(dǎo)致網(wǎng)格的劃分失敗［4］。

圖3 裂紋尖端網(wǎng)格

1.2 載荷邊界條件

在數(shù)值仿真中，載荷邊界條件的設(shè)置原則是盡可能的按實驗的實際情況加載。如圖所示模型的邊界條件為一端固定X、Y方向的自由度，另一端固定Y方向的自由度。載荷加在模型中間裂紋面對應(yīng)的節(jié)點上，具體的載荷邊界條件如圖4所示。

圖4 模型載荷邊界條件

將建立好的模型，在標(biāo)準(zhǔn)模式(Standard)下提交計算，并在分析步中設(shè)置好需要在后處理中得到的斷裂參數(shù)J積分，于是就能方便的在后處理模塊中查看結(jié)果。在歷史輸出中設(shè)置15道圍線積分，Abaqus后處理輸出的數(shù)據(jù)文件查看到J積分，由于J積分與路徑無關(guān)，如圖5所示，因此用15個積分值得出其平均值，結(jié)果比較合理。

圖5 J積分路徑

2 實驗結(jié)果

2.1 仿真結(jié)果及分析

初始裂紋長度為5mm時，載荷為起裂載荷P=1490N，有限元模擬計算得到的應(yīng)力應(yīng)變云圖和變形后X方向的等值線圖如圖6a－c所示:

圖6 應(yīng)力應(yīng)變云圖和變形后X方向的等值線圖

從圖7可以看出有限元模擬得到的J積分隨圍線位置基本不變，但是當(dāng)積分回路非常靠近裂紋尖端時所得到的結(jié)果并不總是嚴(yán)格一致的。因此，一般取裂紋尖端后面幾個穩(wěn)定的積分值比較的合理。可以看出在有效積分區(qū)域內(nèi)J積分的值基本穩(wěn)定，這也驗證了J積分與路徑無關(guān)的理論。

圖7 J積分隨圍線位置的分布

圖8 J積分仿真值與試驗值的對比

從Abaqus的模擬仿真結(jié)果與實驗值對比可以看出，如圖8，這兩種結(jié)果有比較好的一致性。這說明了采用的圍道積分法計算復(fù)合材料的J積分是可行的，本文建立的模型能夠準(zhǔn)確地得到裂紋尖端周圍的J積分，與此同時數(shù)值仿真也進一步地證明了J積分與路徑無關(guān)，此外也驗證了理論上采用J積分來表征CFRP層合板的斷裂韌性特征，得出比較合理的數(shù)據(jù)。

2.2 虛擬裂紋閉合法模擬裂紋擴展

在采用圍道積分方法模擬得到了CFRP層合板起裂點的斷裂韌性值，并且由試驗值的對比可知圍道積分法能比較準(zhǔn)確地計算J積分。但是圍道積分法在模擬斷裂問題時有一個不足的地方就是不能模擬裂紋的擴展［5］。而虛擬裂紋閉合技術(shù)(VCCT)很好的解決了這一個問題，并且Abaqus推出了專門的VCCT模塊。因此我們采用Abaqus自帶的VCCT模塊來模擬CFRP層合板的彎曲斷裂。

虛擬裂紋閉合技術(shù)(VCCT)最先由Rybicki和Kanninen在他們的著名論文［6］里第一次使用。近幾年隨著計算機性能的提升和有限元軟件的應(yīng)用，再加上VCCT的簡潔性和有效性，虛擬裂紋閉合法廣泛各工程領(lǐng)域。VCCT的基本假設(shè)是擴展一定長度列裂紋所需要的能量與閉合相同長度裂紋所需的能量相等，只需計算裂紋閉合和裂紋張開位移就可得到應(yīng)變能量釋放率。

2.3 碳纖維增強樹脂基層合板應(yīng)變能釋放率的計算

VCCT模型如圖9所示，圖中的模型的初始裂紋長度取為5mm，外形尺寸為長度60mm，寬度10mm，加載跨距為40mm。單層板的材料參數(shù)如表1，均勻鋪層單層厚度為0.125mm，共40層，層合板的鋪層參數(shù)為［45°/0°/90°/－45°］s。模型的載荷邊界條件為一端固定x和y方向的位移，另一端固定y方向的位移。該模型的建立主要有一下幾個關(guān)鍵的部分:①定義主、從接觸面;②設(shè)置分裂節(jié)點;③收斂性問題的解決。斷裂準(zhǔn)則采用VCCT的BK準(zhǔn)則，假設(shè)整個模型為平面應(yīng)變問題，本文采用平面應(yīng)變四邊形單元對試件進行對稱模型有限元網(wǎng)格劃分，圖10a－e所示為相對縫深a0/W=0.5時的試件有限元網(wǎng)格圖，單元總數(shù)為4315，節(jié)點總數(shù)為8839。

圖9 有限元模型

圖10 試件有限元網(wǎng)格圖

3 結(jié)論

建立了復(fù)合材料層合板的斷裂模型，在網(wǎng)格劃分時對裂紋尖端處的網(wǎng)格進行了奇異化的處理，在裂紋尖端處建立奇異單元模擬裂紋尖端的奇異性。采用路徑積分法計算了該模型的斷裂韌性J積分，并將模擬得到的J積分值與理論計算得到的J積分對比分析，發(fā)現(xiàn)這兩個值吻合的很好。結(jié)果表明該模型能比較好地模擬復(fù)合材料的斷裂，并能精確的計算J積分。用面狀裂紋的虛擬裂紋閉合法模擬了斷裂過程，雖然該方法能很好地模擬裂紋的擴展過程，但是由于分析前需指定初始裂紋，且只能沿著事先確定的擴展面擴展，不能模擬纖維的拔出等過程。

［1］解德，錢勤，李長安.斷裂力學(xué)中的數(shù)值計算方法及工程應(yīng)用［M］.北京:科學(xué)出版社，2009.

［2］莊茁.Abaqus非線形有限元分析與實例［M］.北京:科學(xué)出版社，2005.

［3］燕偉.鋼筋混凝土梁裂縫寬度的影響因素分析［D］.重慶:重慶交通大學(xué)，2010.

［4］裂縫擴展［EB/OL］.http://blog.sina.com.cn/s/blog_6e0aa6420100rpkx.html.2011 －03 －19.

［5］張正藝，解德，趙耀，等.結(jié)構(gòu)尺度上三維線狀裂紋擴展的數(shù)值模擬［J］.結(jié)構(gòu)強度研究(增刊)，2008，(2).

［6］E.F.Rybiciki，M.F.Kanninen.A finite element calculation of stress intensity factors by a modified crack－closure integral［J］.Engineering Fracture Mechanics，1977，(1).