周 霞, 李 凱, 陳成杭, 陳星馳
(大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室,運載工程與力學(xué)學(xué)部工程力學(xué)系,遼寧 大連 116024)
纖維金屬層合板(Fiber metal laminates, FMLs)是將纖維增強樹脂層與金屬薄板交替鋪排,并在一定溫度和壓力下制備的新型復(fù)合材料,它能充分發(fā)揮纖維復(fù)合材料與金屬材料的力學(xué)特性,具有高比強度、比剛度、優(yōu)良的疲勞性能及抗沖擊損傷性能。以纖維增強鋁合金和鎂合金為代表的該類混雜復(fù)合材料憑借其獨特的力學(xué)性能,在汽車、航空航天和軍工等領(lǐng)域具有廣泛的應(yīng)用前景。隨著現(xiàn)代化進(jìn)程的不斷推進(jìn),在汽車、航空航天和軍工領(lǐng)域中的運載工具和飛行器遭受外物撞擊以及防護工具發(fā)生非貫穿性損傷的現(xiàn)象日益普遍,因此纖維增強金屬層合板的低速沖擊性能研究很有必要。
目前,面向航空航天領(lǐng)域的玻璃纖維增強鋁板已有較多的研究和應(yīng)用[1- 2],相比鋁合金,鎂合金具有更低的密度,更高的比強度及比剛度,用碳纖維或碳纖維與玻璃纖維混雜增強的鎂合金層合板有望成為下一代的航空航天替代材料[3]。近年來,已有學(xué)者對纖維增強鎂合金層合板展開了研究,Cortés等[4]分別對玻璃纖維和碳纖維增強的AZ31鎂合金層合板進(jìn)行了落錘沖擊試驗,研究了層合板的分層和剪切破壞情況,發(fā)現(xiàn)鎂基纖維層合板比鋁基纖維層合板具有更優(yōu)越的沖擊能吸收能力。P?rn?nen等[5]和Múgica等[6]對纖維增強鎂基層合板和纖維增強鋁基層合板的低速沖擊性能進(jìn)行了對比研究,發(fā)現(xiàn)纖維鎂合金層合板雖然在較低的沖擊能時就會開裂,但它與纖維鋁合金層合板的侵徹極限在同等水平。盡管已有國內(nèi)外學(xué)者對纖維鎂合金層板的基本力學(xué)性能和沖擊性能進(jìn)行了一些研究[7-8],但對其在沖擊載荷下的動態(tài)響應(yīng)和損傷變化過程的研究探索較少,對混雜纖維增強的鎂合金層合板的研究還未見報道。
本文通過ABAQUS/Explicit分析算法,分別對玻璃纖維/環(huán)氧(Glass Fiber Reinforced Epoxy Resin, GF)、碳纖維/環(huán)氧(Carbon Fiber Reinforced Epoxy Resin, CF)和二者不同混雜方式(簡記為:H C-G1、H C-G2)增強的四種鎂合金層合板在低速沖擊下的動態(tài)響應(yīng)及損傷破壞過程進(jìn)行了數(shù)值模擬。基于鎂合金各向異性屈服本構(gòu)、界面指數(shù)形式的內(nèi)聚力本構(gòu)以及纖維復(fù)合材料的三維Hashin失效準(zhǔn)則,采用 FORTRAN 語言編寫了纖維/樹脂復(fù)合材料層及粘結(jié)界面層本構(gòu)及損傷模型的用戶材料子程序(VUMAT),進(jìn)而采用ABAQUS 自帶的材料模型和本文建立的復(fù)合材料三維損傷力學(xué)模型,研究了在不同落錘沖擊能下,纖維增強鎂合金層合板的能量吸收、變形損傷以及錘頭接觸力的變化。在模擬與試驗結(jié)果吻合良好的基礎(chǔ)上,還對不同鋪層形式的CF /GF /AZ31鎂合金混雜復(fù)合材料的抗沖擊性能進(jìn)行了模擬研究。本研究為準(zhǔn)確評價新型纖維增強鎂合金層合板復(fù)合材料的沖擊行為及其性能優(yōu)化提供了合理的依據(jù)。
為了驗證FMLs低速沖擊響應(yīng)預(yù)測方法的正確性,本文參考P?rn?nen等研究中的試驗樣品建立了同樣幾何尺寸的纖維鎂層合板分析模型,該層合板由3層AZ31B鎂板和2層纖維復(fù)合材料層(4個單向纖維層構(gòu)成1層)交錯堆疊而成,其結(jié)構(gòu)如圖1(a)所示。層合板為直徑80 mm的圓盤,其中鎂合金層厚度為0.5 mm,每2層鎂合金板之間的纖維復(fù)合材料層由4層單層板按鋪層方式鋪設(shè)而成,纖維單層板厚度為0.1 mm,總厚度為0.4 mm,纖維/樹脂層鋪設(shè)方向與鎂合金薄板軋制方向相同。根據(jù)不同纖維層材料,上述纖維鎂層合板分析模型包括玻璃纖維增強鎂板(GF/Mg)、碳纖維增強鎂板(CF/Mg)和二者混雜增強鎂板(H C-G1/Mg, H C-G2/Mg)四種結(jié)構(gòu),不同纖維增強鎂合金層合板的纖維鋪層方式見表1,其中玻璃纖維材料為FM94/S2-玻璃纖維/樹脂復(fù)合材料,碳纖維材料為HTS40 12K碳纖維/樹脂復(fù)合材料。
表1 纖維增強鎂合金FMLs的纖維鋪層方式
纖維/鎂層合板有限元模型如圖1(b)所示,纖維層與鎂合金層均采用三維實體單元建模,單元類型為8節(jié)點實體單元(C3D8R),單層板之間則插入一層零厚度的界面內(nèi)聚力單元(COH3D8)來模擬層合板中粘結(jié)層(Cohesive Layer)的粘結(jié)及層間分層作用。為了減小計算量及同時保證計算結(jié)果精度,對與沖擊彈頭直接接觸的層合板區(qū)域的網(wǎng)格進(jìn)行了細(xì)化處理,層合板每層最小網(wǎng)格尺寸為0.83 mm×0.83 mm。層合板模型包含了3層鎂合金層、6層單向鋪設(shè)的纖維層(相鄰的兩個90°鋪設(shè)層用一層單元模擬)及8層連接界面cohesive層,層合板中每個單層的網(wǎng)格分布一致,表2為層合板模型中的單元和節(jié)點數(shù)。
(a) 鎂基纖維層板橫截面
(b) 鎂基纖維層板有限元模型圖1 纖維/鎂合金層合板截面及其有限元網(wǎng)格模型Fig.1 Cross-section and FE mesh model of fiber/Mg FMLs
本文模擬中,約束板外圈各節(jié)點的所有自由度,沖頭垂直于板的中心撞擊。模型中剛性半球形沖頭的質(zhì)量為2.78 kg,直徑為15.9 mm,通過改變沖頭的速度來獲得層合板在不同沖擊能下的動態(tài)力學(xué)響應(yīng),沖頭下落速度分別為3.793 m/s、4.646 m/s及5.998 m/s,對應(yīng)的動能為20 J、30 J及50 J。此外,正碰撞問題中摩擦的影響較小[9],且在P?rn?nen等提及的試驗中在沖擊位置放置了塑料薄膜來進(jìn)行防摩擦處理,因此,模擬中不考慮沖擊摩擦的影響。
表2 模型的節(jié)點數(shù)和單元數(shù)
FMLs在低速沖擊過程中,其主要的損傷失效模式可分為鎂合金層和纖維/樹脂層的層內(nèi)損傷及其層間損傷。因此,各單層的層內(nèi)損傷與層間損傷本構(gòu)模型和失效準(zhǔn)則是能否準(zhǔn)確地預(yù)測纖維增強鎂合金層合板低速沖擊響應(yīng)的關(guān)鍵。
鎂合金薄板在軋制方向與橫向存在較為明顯的各向異性,本文采用被廣泛使用的Hill48正交各向異性屈服[10]準(zhǔn)則,其勢函數(shù)表示為:
(1)
式中:σij為應(yīng)力分量;F,G,H,L,M,N分別為材料的各向異性參數(shù),用lankford’s系數(shù)r0,r45,r90表示,r0,r45,r90分別表示板材在拉伸實驗中拉伸方向與板材軋制方向成0°、45°、90°時的值。
(2)
對于薄板材料,各向異性參數(shù)L,M不能由單向拉伸試驗測得,本文取L=M=N。當(dāng)F=G=H=L/3=M/3=N/3=1/2時,Hill48各向異性屈服準(zhǔn)則退化為各向同性的Mises屈服準(zhǔn)則。表3為本文所用的AZ31B鎂合金板各向異性參數(shù)。
表3 AZ31B鎂合金各向異性參數(shù)[11]
對鎂合金損傷破壞過程的模擬,同時采用了Johnson-Cook (J-C) 失效本構(gòu)模型和延性損傷準(zhǔn)則。J-C失效本構(gòu)模型:
(3)
式中:D為材料的損傷因子。D初始值為0,當(dāng)材料發(fā)生損傷時,D>0,隨著損傷破壞的加劇,D逐漸增大到1,當(dāng)D=1時,材料破壞單元刪除;Δε為當(dāng)前時間步的等效塑性應(yīng)變增量;εf為當(dāng)前時間步的失效應(yīng)變,其表達(dá)式為:
(4)
延性損傷準(zhǔn)則可表示為:
(5)
表4 AZ31B鎂合金J-C損傷參數(shù)[12]
三維 Hashin 準(zhǔn)則是常用的纖維/樹脂復(fù)合材料損傷失效判據(jù)[13],本文纖維層采用3D-Hashin[14]損傷失效模型,該損傷模型主要考慮四種損傷模式:纖維拉伸、纖維壓縮、基體拉伸、基體壓縮。
纖維拉伸:σ11≥0
(6)
纖維壓縮:σ11<0
(7)
基體拉伸:σ22+σ33≥0
(8)
基體壓縮:σ22+σ33<0
(9)
損傷演化:
當(dāng)材料滿足以上四種損傷模式的任意一種時,有實驗表明材料并未完全失效,仍具有承載能力,其承載能力的下降主要體現(xiàn)為宏觀彈性模量的減小,本文按照表5所示的準(zhǔn)則進(jìn)行剛度折減[15]。材料中的單元達(dá)到某一種損傷模式后仍繼續(xù)承載,當(dāng)纖維達(dá)到最大應(yīng)變時發(fā)生纖維破壞,基體達(dá)到最大應(yīng)變時發(fā)生基體破壞,單元滿足任意一種破壞時,表明該單元已徹底失去承載能力,應(yīng)從模型中刪除該單元?;谝陨侠w維增強復(fù)合材料的損傷失效模型,編寫用戶材料子程序VUMAT來實現(xiàn)纖維樹脂層破壞的模擬。纖維/環(huán)氧樹脂復(fù)合材料性能參數(shù),分別見表6和表7。
表5 本文所用的剛度折減方案[15]
表6 不同纖維/環(huán)氧樹脂復(fù)合材料的彈性參數(shù)[19-20]
表7 不同纖維/環(huán)氧樹脂復(fù)合材料的強度參數(shù)[19-20]
單元刪除判據(jù):
(10)
式中:ψ為單元破壞狀態(tài)變量,等于0時單元破壞并刪除,等于1時單元未刪除。
本文采用指數(shù)形式的粘結(jié)表面內(nèi)聚力-相對位移函數(shù)關(guān)系(Transction-Sepration Law, T-S Law),對鎂合金層和纖維樹脂層之間脫層的初始和演化進(jìn)行模擬。指數(shù)T-S關(guān)系的內(nèi)聚力(cohesive)模型最先由Xu等[16]提出,之后Bosch等[17]研究發(fā)現(xiàn)其中存在一些不足后對其改進(jìn),提出以下T-S關(guān)系:
(11)
(12)
式中:δn,φn,Δn分別為法向開裂初始位移、破壞斷裂能、當(dāng)前位移;δt,φt,Δt分別為切向開裂初始位移、破壞斷裂能、當(dāng)前位移。
判斷單元是否刪除的依據(jù)是當(dāng)前開裂能是否達(dá)到最大開裂能,對單元當(dāng)前所吸收開裂能的計算,可通過開裂應(yīng)力對相應(yīng)開裂位移積分求得:
(13)
(14)
(15)
(16)
式中:ψ*為單元破壞狀態(tài)變量,當(dāng)其小于1時,單元未破壞,當(dāng)其大于等于1時,單元破壞并刪除;系數(shù)0.95為斷裂能截斷系數(shù)。由于本文采用指數(shù)形式的cohesive模型,界面吸收的斷裂能會在破壞過程中逐步的無限接近材料開裂斷裂能,并不會等于開裂斷裂能,因此需要在某個位置進(jìn)行截斷,本文取0.95。指數(shù)關(guān)系的T-S界面本構(gòu)關(guān)系通過編寫VUMAT子程序?qū)崿F(xiàn),并結(jié)合纖維復(fù)合材料的子程序在ABAQUS/Explicit中同時用于層合板的沖擊模擬。脫粘的初始和演化的損傷參數(shù),見表8。
表8 粘接層參數(shù)[21, 22]
沖擊能為20 J時(圖2),從模擬結(jié)果中可以看出,第一層(沖擊面)、第二層鎂合金出現(xiàn)了較大的變形、鼓包現(xiàn)象,僅第三層(非沖擊面)鎂合金出現(xiàn)輕微開裂,與試驗中僅在非沖擊面出現(xiàn)裂縫結(jié)果一致。沖擊能為50 J時(圖3),在沖擊面鎂合金呈花型開裂,在非沖擊面,一條較長的主裂縫貫穿沖擊作用點,并在主裂縫中有一條橫向的小裂縫。從鎂合金層的破壞形貌來看,模擬結(jié)果與試驗基本吻合。
圖2 沖擊能為20 J 時,模擬預(yù)測的GF/Mg形貌及其與P?rn?nen等實驗結(jié)果的比較 Fig.2 Comparison of the damage morphology between the experiment in P?rn?nen et al and numerical simulations for GF/Mg under a impact energy of 20 J
圖3 沖擊能為50 J時的GF/Mg損傷形貌數(shù)值預(yù)測及其與P?rn?nen等實驗結(jié)果的比較Fig.3 The predicted damage morphology and comparison with the experiments by P?rn?nen et al. for GF/Mg under 50 J
圖4為層合板沖擊過程中的載荷-時程對比曲線。從曲線中可以看出,沖擊能為20 J時,模擬與試驗曲線兩者吻合良好,沖擊力的峰值基本一致,誤差僅為3.1%,試驗曲線走勢較為光滑,模擬曲線在沖擊開始后不久由于材料第一次由彈性進(jìn)入塑性階段而出現(xiàn)了短暫的平動及下降走勢,在峰值處出現(xiàn)了較大的下降,模擬曲線與P?rn?nen等[5]的試驗結(jié)果走勢一致。在30 J沖擊能下,載荷-時程曲線仍能保持較高的吻合度,沖擊力峰值誤差為6.6%。當(dāng)沖擊能增大到50 J時,曲線走勢、沖擊過程持續(xù)時間仍能保持一致,僅在沖擊的后半段,模擬的沖擊力與試驗相比有一定的偏差,這可能與沖擊實驗過程中鎂合金表面并非十分光滑,以及鎂合金在高能沖擊加載下的變形和斷裂機理比較復(fù)雜有關(guān)。沖擊加載下鎂合金的變形失效不僅與溫度、應(yīng)變量和應(yīng)變速率有關(guān),還與加載方向、組織演化如變形孿晶、絕熱剪切帶形成有關(guān),這些因素的綜合作用可能導(dǎo)致試樣在沖擊過程中失效較慢,而模擬并未考慮上述因素的綜合作用。
圖4 不同沖擊能(20 J、30 J及50 J) 時, GF/ Mg載荷-時程曲線的模擬預(yù)測與P?rn?nen等[5]實驗結(jié)果的比較Fig.4 Comparison of force vs. time relations between the experiments by P?rn?nen et al. and simulation predictions for GF/ Mg under different impact energies of 20 J, 30 J and 50 J
圖5為載荷-撓度曲線的模擬預(yù)測與試驗結(jié)果的對比。沖擊能為20 J時,模擬與試驗曲線整體走勢基本一致,試驗中板的中心最大撓度為5.72 mm,模擬得到的最大撓度為5.62 mm,誤差僅為1.85%。當(dāng)沖擊能為30 J時,層合板的最大撓度試驗值為7.41 mm,模擬值為7.08 mm,誤差為4.36%。沖擊能達(dá)到50 J時,兩者曲線吻合度稍差,試驗中最大撓度為10.29 mm,模擬結(jié)果為9.36 mm,誤差為9.35%。
圖5 不同沖擊能(20 J、30 J及50 J)時,GF/ Mg載荷-撓度曲線的模擬預(yù)測與P?rn?nen等實驗結(jié)果的比較Fig.5 Comparison of force vs. central deflection relations between the experiments by P?rn?nen et al. and simulation predictions for GF/ Mg under different impact energies of 20 J, 30 J and 50 J
圖6為20 J沖擊能下三種鋪層形式的纖維鎂合金層板的損傷形貌圖。結(jié)合圖2中GF/Mg的損傷形貌,可以發(fā)現(xiàn),在單一碳纖維增強鎂合金復(fù)合材料中添加玻璃纖維層,且玻璃纖維的鋪設(shè)位置靠近下表面時,對應(yīng)的混雜層合板H C-G1/Mg中的三層鎂合金層損傷面積均較小,說明玻璃纖維的存在可以提高CF/Mg復(fù)合材料的抗沖擊能力。圖7 為50 J沖擊能量下單一碳纖維增強鎂合金復(fù)合材料CF/Mg和混雜層合板H C-G1/Mg中各纖維鋪層沖擊后的破壞形貌??梢钥闯?,對于CF/Mg和H C-G1/Mg,它們的第1層碳纖維受到?jīng)_擊載荷后破壞均較嚴(yán)重,較多的單元被刪除。比較而言,H C-G1/Mg中下層玻璃纖維的破壞較上層碳纖維層較輕,最低層玻璃纖維只出現(xiàn)較小破壞,從而限制了沖頭完全穿透層合板,提高了材料耐沖擊強度。
圖8為沖擊能為50 J時層合板破壞后的橫截面形貌,可以看出,添加了玻璃纖維的混雜層合板雖然也產(chǎn)生了纖維斷裂,但纖維斷裂伸長較小,且填加玻璃纖維層的H C-G1/Mg混雜復(fù)合材料完整性較好,而全部為玻璃纖維的層合板出現(xiàn)了穿透現(xiàn)象,說明單一玻璃纖維增強鎂合金層合板在較高能沖擊時較單一碳纖維增強層合板抗沖擊能力差,但二者混雜纖維增強的H C-G1/Mg層合板在50 J時仍具有較好的抗沖擊性能。
圖9表示了單一碳纖維鎂合金層合板CF/Mg與混雜纖維增強的H C-G1/Mg層合板在50 J沖擊能作用下各個粘結(jié)層的脫粘情況對比??梢钥闯鲈诙N層合板的鎂合金板中間的粘結(jié)層即第2、3粘結(jié)層及第6、7粘結(jié)層的脫粘情況均較為嚴(yán)重。二種層合板比較而言,雖然在沖擊能較高時粘結(jié)層均遭受不同程度的破壞,但在靠近下表層鋪設(shè)玻璃纖維的H C-G1/Mg混雜層合板,通過混雜纖維的相互作用能夠減少靠近下表面的粘結(jié)層(如第7粘結(jié)層)的脫粘面積,因而有一定的增韌效果和抗沖擊效果。
圖6 20 J沖擊能下鎂合金層的損傷形貌Fig.6 The damage morphology of magnesium alloy layers under impact energy of 20 J
圖7 50 J沖擊能下纖維層的損傷形貌Fig.7 Damage morphology of fiber layers under impact energy of 50 J
圖8 50 J沖擊能下不同纖維/鎂層合板破壞后的斷口形貌Fig.8 Damage morphology of different fiber/Mg alloy composites laminates
圖9 50 J沖擊能下界面脫粘損傷Fig.9 Interface debonding damage under impact energy of 50 J
圖10為四種不同纖維鎂合金FMLs在沖擊過程中的載荷-時程曲線。從曲線中可以看出,四種FMLs的曲線走勢比較接近,沖擊過程中的最大沖擊載荷相差較小。CF/ Mg沖擊力曲線上升的斜率最大,兩種混雜增強的層合板次之, GF/Mg層合板沖擊力曲線上升最慢,這是由于碳纖維復(fù)合材料層比玻璃纖維層的剛度更大造成的,混雜增強的層合板的沖擊力曲線走勢介于兩者之間,但更接近與碳纖維增強的鎂合金層合板。圖11為四種纖維鎂合金FMLs在沖擊過程中的載荷-位移曲線。在20 J、30 J、50 J沖擊能下,四種FMLs的載荷-位移曲線走勢基本一致,但相比之下,GF/Mg在沖擊過程中的變形及沖擊后的殘余撓度均最大,CF/Mg層合板的殘余撓度則最小,兩種混雜纖維鎂合金FMLs的曲線基本介于兩者之間,且沖擊能較大時,H C-G1/Mg復(fù)合材料具有較大的峰值沖擊載荷和較小的殘余變形。
圖10 不同沖擊能(20 J、30 J及50 J)時,四種不同纖維鎂合金FMLs的載荷-時程曲線模擬結(jié)果Fig.10 Simulation of force vs. time relations for four fiber/Mg FMLs under different impact energies of 20 J, 30 J and 50 J
圖11 不同沖擊能(20 J、30 J及50 J)時,四種不同纖維鎂合金FMLs的載荷-撓度曲線模擬結(jié)果Fig.11 Simulation of force vs. central deflection relations for four fiber/Mg FMLs under different impact energies of 20 J, 30 J and 50 J
圖12表示了纖維鎂合金層板在不同沖擊能量下的吸能-時程曲線??梢钥闯?,在低速沖擊中,當(dāng)沖擊能分別為20 J、30 J、50 J時,層合板吸收能量的比例隨沖擊能增加而增大,沖擊結(jié)束時,沖頭回彈速度隨沖擊能增加而減小。對比四種不同纖維鎂合金FMLs的低速沖擊模擬結(jié)果還發(fā)現(xiàn),隨著沖擊能增加,當(dāng)沖擊能大于30 J時,單一玻璃纖維增強的層合板最終吸收能要大于單一碳纖維增強的層合板,吸能現(xiàn)象能夠增加纖維和基體的損傷及鎂合金層的損傷失效,但產(chǎn)生最高吸能的時刻滯后。由于玻璃纖維的韌性要遠(yuǎn)高于碳纖維,因此采用混雜碳纖維與玻璃纖維增強的鎂合金層合板(H C-G1/Mg),在研究的沖擊能量范圍內(nèi),能在保持沖擊撓度較小的同時,也有較好的吸能特性,因而能改善碳纖維增強鎂合金復(fù)合材料的抗沖擊性。
圖12 四種不同纖維鎂合金FMLs的吸能-時程曲線 Fig.12 Absorbed energy vs. time curves of four fiber/Mg FMLs
(1) 根據(jù)鎂合金塑性各向異性本構(gòu)、纖維層的三維Hashin失效準(zhǔn)則以及粘接界面的指數(shù)內(nèi)聚力模型編寫的復(fù)合材料損傷子程序,能較好模擬再現(xiàn)試驗中纖維增強鎂合金層合板的沖擊響應(yīng)和損傷過程。纖維增強鎂合金層合板在落錘低速沖擊中,隨著沖擊能量的增加,非沖擊面首先出現(xiàn)裂縫,說明纖維鎂合金層合板中鎂合金層性能對層合板低速沖擊性能影響較大。
(2) 在沖擊能為20 J、30 J、50 J時,不同纖維鎂合金層合板的金屬損傷、纖維破壞面積、脫粘面積及板的中心最大撓度均隨沖擊能的增加而變大,但回彈能量卻隨初始沖擊能能量的增加而減小。當(dāng)沖擊能量大于30 J后,單一纖維鎂合金層板的沖擊能量閥值要大于單一玻璃纖維鎂合金層合板。
(3) 由于玻璃纖維增強鎂合金復(fù)合材料的高沖擊韌性,在單一碳纖維增強的鎂合金層合板中添加玻璃纖維以后,混雜層合板的抗低速沖擊能力提高。在沖擊能量為20~50 J范圍內(nèi),當(dāng)玻璃纖維位于混雜復(fù)合材料第2 層時,玻璃纖維/碳纖維/鎂合金復(fù)合材料的沖擊韌性最好,材料整體既有較高的抗沖擊損傷能力。