趙士成,王振清,郭建明,楊斌
(哈爾濱工程大學(xué)航天與建筑工程學(xué)院,黑龍江哈爾濱150001)
混雜復(fù)合材料是指2種或2種以上的纖維增強(qiáng)同一種樹脂基體而構(gòu)成的復(fù)合材料[1]。與傳統(tǒng)增強(qiáng)纖維相比,碳纖維具有高強(qiáng)度、高模量等特性。然而,碳纖維作為增強(qiáng)纖維的沖擊韌性相對(duì)于其他纖維(如玻璃纖維)增強(qiáng)復(fù)合材料低得多[2]。與玻璃纖維的混雜,可以優(yōu)勢(shì)互補(bǔ),從而使得到的復(fù)合材料具有較佳的綜合性能。
國內(nèi)外學(xué)者對(duì)環(huán)氧樹脂基混雜復(fù)合材料的研究開展不久[3-6]。國外也有很多學(xué)者開始對(duì)混雜復(fù)合材料層合板的沖擊性能做了研究[7-9]。但目前對(duì)碳纖維和玻璃纖維共同增強(qiáng)環(huán)氧樹脂基復(fù)合材料的低速?zèng)_擊性能的研究還很少。本文采用ABAQUS有限元軟件低速?zèng)_擊載荷下的相應(yīng)情況進(jìn)行了仿真。自20世紀(jì)90年代以后,漸進(jìn)損傷破壞分析理論發(fā)展起來,結(jié)合剛度退化理論,為預(yù)測(cè)層合板各種損傷的起始和發(fā)展提供了一個(gè)較好的方法和途徑[10-13]。因此,本文編寫了用戶材料子程序VUMAT,并且基于表面內(nèi)聚力方法,考慮纖維界面的分層破壞,建立了混雜玻璃纖維增強(qiáng)環(huán)氧樹脂低速?zèng)_擊的數(shù)值模型,模擬了混雜層合板在落錘低速?zèng)_擊下能量吸收、變形以及錘頭接觸力的變化。同時(shí),對(duì)鋪層形式對(duì)CF/GF/環(huán)氧混雜復(fù)合材料的影響進(jìn)行了分析。
采用復(fù)合材料的三維Hashin失效準(zhǔn)則編寫了用于顯式計(jì)算的VUMAT用戶子程序。該失效準(zhǔn)則綜合考慮了材料在破壞時(shí)X、Y、Z這3個(gè)方向的拉伸和壓縮失效。圖1給出了復(fù)合材料的三維受力示意圖,即是在計(jì)算時(shí)考慮單元受到的正應(yīng)力及剪切應(yīng)力。
圖1 單向?qū)雍习迨芰κ疽鈭DFig.1 Schematic diagram of stress in composite laminates
具體的失效判據(jù)如下,認(rèn)為當(dāng)3個(gè)方向的判定系數(shù)大于1時(shí)材料即失效:
X向拉壓(纖維拉壓破壞):
Y向拉壓(基體斷裂與屈曲):
Z向拉壓(分層與壓縮):
當(dāng)單元受到的應(yīng)力滿足式(1)~(6)中的材料3個(gè)方向的系數(shù)任意一個(gè)時(shí)即認(rèn)為單元失效。對(duì)于失效的單元,Hashin認(rèn)為不可直接將其刪除,而是認(rèn)為單元的各剛度參數(shù)在一定程度上衰減,本文采用表1[14]中列出的剛度衰減方案進(jìn)行計(jì)算。
表1 失效單元的剛度衰減方案Table 1 Stiffness reduction methods of failure elements
本文分別選取碳纖維和玻璃纖維增強(qiáng)環(huán)氧樹脂層合板的三維力學(xué)參數(shù)進(jìn)行計(jì)算。層合板的彈性參數(shù)其強(qiáng)度在表2、3[15]中給出,沖頭的力學(xué)參數(shù)如下:E=210 GPa,μ=0.3,ρ=8 859 kg/m3。
表2 不同纖維增強(qiáng)環(huán)氧樹脂基復(fù)合材料的彈性參數(shù)Table 2 Elastic parameters of different fiber reinforced epoxy composites GPa
表3 不同纖維增強(qiáng)環(huán)氧樹脂基復(fù)合材料的強(qiáng)度參數(shù)Table 3 Strength of different fiber reinforced epoxy composites MPa
實(shí)際沖擊過程中,沖頭的剛度與被沖擊件相比要大很多,因此,在仿真時(shí)可以將沖頭看做剛體。本文在通用有限元軟件ABAQUS/Expolicit的基礎(chǔ)上,進(jìn)行了混雜復(fù)合材料低速?zèng)_擊作用下的動(dòng)態(tài)響應(yīng)分析。層合板直徑80 mm,四周全約束,沖頭直徑12 mm,速度5 m/s,可以計(jì)算出系統(tǒng)初始動(dòng)能為20 J。為了滿足強(qiáng)度要求,混雜層合板中玻璃纖維層不大于4層。因此,仿真過程中按照表5的不同混雜方式共建立了5種有限元模型,單層厚度0.192 mm,鋪層角度為復(fù)合材料層合板在低速低能沖擊下的損傷破壞模型表現(xiàn)為基體開裂、基體擠壓、纖維斷裂、分層。而在這些破壞形式中,分層損傷使復(fù)合材料的強(qiáng)度和壽命大幅下降,嚴(yán)重影響材料的使用。本文通過在每層復(fù)合材料的層間加入基于表面內(nèi)聚力行為單元,用于模擬界面分層擴(kuò)展。基于表面內(nèi)聚力行為單元的本質(zhì)是在材料之間加入與內(nèi)聚力單元方法相似的牽引分離本構(gòu)模型。同常用的裂紋分析方法相比(如內(nèi)聚力單元方法),該方法的優(yōu)勢(shì)在于不用像內(nèi)聚力方法那樣在界面處插入內(nèi)聚力單元,從而避免了內(nèi)聚力單元的存在和單元極度扭曲,提高了計(jì)算的收斂性。為提高計(jì)算精度,將與沖頭接觸的層合板網(wǎng)格細(xì)化,所建立的有限元模型如圖2所示。
表4 GF/CF/環(huán)氧復(fù)合材料混雜鋪層形式Table 4 Layer form of GF/CF/Epoxy hybrid composites
圖2 有限元模型Fig.2 The finite element model
層合板的臨界穿透能是衡量材料抵御沖擊能力的重要指標(biāo),它定義為材料剛好被完全穿透時(shí)所需要的沖擊能量。本文通過對(duì)5種不同形式的層合板進(jìn)行計(jì)算,所得到的接觸力-變形,能量-時(shí)間曲線如圖3所示。
圖3 5種層合板的接觸力-變形,能量-時(shí)間曲線Fig.3 Contact force-displacement and energy-time curves of different specimens
通過圖3可以看出,隨著玻璃纖維布鋪設(shè)位置的不同,試件承受沖擊載荷時(shí)的接觸力-變形曲線和能量-時(shí)間曲線有很大不同。從能量-時(shí)間曲線上可以看出,沖擊能量隨著接觸時(shí)間增加而減小,當(dāng)接觸時(shí)間為2.41 ms時(shí),沖頭的沖擊動(dòng)能衰減到0。而此時(shí),通過比較圖4中所列的層合板破壞形式,第1種材料剛好開始出現(xiàn)破壞,沖頭剛好能夠穿透層合板。說明對(duì)于第1種復(fù)合材料層合板來說,20 J的沖擊動(dòng)能即為其臨界穿透能。從接觸力-變形曲線來說,層合板中沒有玻璃纖維層存在時(shí),接觸力-變形曲線幾乎閉合,說明沖擊過程中20 J的沖擊能量由材料幾乎完全吸收。然而,當(dāng)添加玻璃纖維以后,接觸力-變形曲線呈現(xiàn)出半封閉狀態(tài),沖頭沒有完全穿透層合板,說明材料還可以通過進(jìn)一步的變形來吸收較多的沖擊動(dòng)能,從而來提高材料的抗沖擊能力。
混雜復(fù)合材料臨界穿透能的提高可以歸因于玻璃纖維的存在。眾所周知,玻璃纖維的韌性要比碳纖維高很多,因此,混雜復(fù)合材料在承受沖擊載荷時(shí),材料變形增加,材料通過大變形與沖頭的接觸時(shí)間變長(zhǎng),從而接觸力與變形圍成曲線的面積增加,宏觀上表現(xiàn)為混雜復(fù)合材料抗沖擊性能的提升。
圖4 沖頭能量衰減為0時(shí)層合板形狀(t=2.41 ms)Fig.4 Shape of the laminates when the impact energy decline to zero(t=2.41 ms)
圖5為5種不同形式的層合板在沖擊能量降為0時(shí)的破壞形貌。整體來看,隨著玻璃纖維層的鋪設(shè)位置靠近上下表面,層合板的破壞面積逐漸減小,說明玻璃纖維的存在對(duì)提升材料抗沖擊能力發(fā)揮了很大的作用。圖6為層合板破壞后的橫觀形貌,可以看出,未添加玻璃纖維的層合板完全被穿透,而添加玻璃纖維層的復(fù)合材料完整性較好,沖擊韌性提升。
圖5 層合板沖擊破壞面積Fig.5 Areas of composite laminates after impact damage
圖6 不同形式的層合板破壞后的斷口形貌Fig.6 Damage morphology of different composites laminates
由圖6可以看出,沖擊性能的大小順序?yàn)樾问?>形式4>形式3>形式2>形式1。碳纖維在承受沖擊載荷時(shí)表現(xiàn)出的特點(diǎn)是斷裂伸長(zhǎng)低,在承受很高應(yīng)力時(shí)伸長(zhǎng)也很小。作為復(fù)合材料的增強(qiáng)材料,伸長(zhǎng)小,意味著沖擊性能差。當(dāng)玻璃纖維與碳纖維混雜后,2種纖維之間的相互作用對(duì)混雜復(fù)合材料的變形產(chǎn)生很大的影響。隨著玻璃纖維的加入,由于玻璃纖維增強(qiáng)復(fù)合材料的強(qiáng)度低但韌性好,從而導(dǎo)致沖擊破壞引發(fā)能的減小,混雜層合板吸收的能量增加。所以,通過與玻璃纖維混雜,碳纖維增強(qiáng)復(fù)合材料沖擊性能提高的實(shí)質(zhì)是提高了材料的沖擊擴(kuò)展能。
圖7 鋪層形式為[C/G2/C3]S的混雜復(fù)合材料各層材料破壞情況Fig.7 Damage morphology of hybrid composite laminates with layer form[C/G2/C3]S
圖7為第5種混雜層合板的各層層合板的沖擊后的破壞形貌??梢钥闯觯?層玻璃纖維受到?jīng)_擊載荷后破壞嚴(yán)重,較多的單元被刪除。比較而言,上層碳纖維的破壞較下層嚴(yán)重,最低層碳纖維只出現(xiàn)較小破壞,從而限制了沖頭完全穿透層合板,提高了材料耐沖擊強(qiáng)度。可以看到,通過靠近上下表層鋪設(shè)玻璃纖維能夠很大程度地提升材料整體韌性,這是因?yàn)椴牧显谑艿經(jīng)_擊時(shí),玻璃纖維通過斷裂而將大部分沖擊動(dòng)能吸收,這種斷裂是通過受到?jīng)_擊時(shí)上表面玻璃纖維的壓縮斷裂和下表面玻璃纖維的拉伸斷裂產(chǎn)生的。然而,當(dāng)玻璃纖維層位于遠(yuǎn)離混雜復(fù)合材料表層時(shí),玻璃纖維受到的拉壓載荷小,產(chǎn)生的應(yīng)變小,纖維不能夠通過斷裂吸收沖擊動(dòng)能,從而增韌效果不明顯。
通過以上對(duì)混雜復(fù)合材料層合板低速?zèng)_擊的有限元分析方法,可以得出以下結(jié)論:
2)不同混雜方式的復(fù)合材料在承受沖擊載荷時(shí),吸收的能量不同,玻璃纖維靠近層合板表面鋪設(shè)時(shí)吸收的能量最多。
3)沖擊后的混雜復(fù)合材料的破壞面積隨著玻璃纖維鋪層靠近混雜層合板上下表面而減小,當(dāng)玻璃纖維位于混雜復(fù)合材料第2層時(shí),玻璃纖維/碳纖維/環(huán)氧樹脂復(fù)合材料的沖擊韌性最好,材料整體破壞面積最小。
[1]FARIBORZ S J,HARLOW D G.The tensile behavior of intraply hybrid compositesⅡ:Micromechanical model[J].Journal of Composite Materials,1987,21(9):856-875.
[2]趙謙,王善元.單向混雜復(fù)合材料拉伸性能的兩級(jí)分析[J].南京理工大學(xué)學(xué)報(bào),2007,31(6):766-770.
ZHAO Qian,WANG Shanyuan.Two-step analysis of tensile behaviors of unidirectional hybrid composites[J].Journal of Nanjing University of Science and Technology,2007,31 (6):766-770.
[3]張軍,張寶慶,丁艷芬,等.聚丙烯/熱致液晶聚合物/玻璃纖維混雜復(fù)合材料Ⅰ:液晶聚合物組分的預(yù)成纖及對(duì)混雜復(fù)合材料形態(tài)的控制[J].復(fù)合材料學(xué)報(bào),2004,21 (6):1-8.
ZHANG Jun,ZHANG Baoqing,DING Yanfen,et al.PP/ TLCP/Glass fiber hybrid compositesⅠ:The pregeneration of TLCP fibrils and control of the morphology of hybrid composites[J].Acta Materiae Compositae Sinica,2004,21(6): 1-8.
[4]張寶慶,張軍,丁艷芬,等.聚丙烯/熱致液晶聚合物/玻璃纖維混雜復(fù)合材料Ⅱ:混雜復(fù)合材料的靜態(tài)及動(dòng)態(tài)力學(xué)性能[J].復(fù)合材料學(xué)報(bào),2004,21(6):9-13.
ZHANG Baoqing,ZHANG Jun,DING Yanfen,et al.PP/ TLCP/Glass fiber hybrid compositesⅡ:Tensile and dynamic mechanical properties of hybrid composites[J].Acta Materiae Compositae Sinica,2004,21(6):9-13.
[5]張用兵,史俊虎,王利.混雜方式對(duì)CF/UHMWPEF混雜復(fù)合材料力學(xué)性能的影響[J].玻璃鋼/復(fù)合材料,2005 (2):17-19.
ZHANG Yongbing,SHI Junhu,WANG Li.Effect of hybrid mode on CF/UHMWPEF composite performance[J].FRP/ CM,2005(2):17-19.
[6]張梅.單向芳綸/玻璃纖維混雜復(fù)合材料板材拉伸性能研究[J].玻璃鋼/復(fù)合材料,2007(6):35-37,52.
ZHANG Mei.Tensile property study of unidirectional Kevlar/glass fiber hybrid composites[J].FRP/CM,2007(6): 35-37,52.
[7]DEHKORDI M T,NOSRATY H,SHOKRIEH M M,et al.The influence of hybridization on impact damage behavior and residual compression strength of intraply basalt/nylon hybrid composites[J].Materials&Design,2013,43:283-290.
[8]ATIQAH A,MALEQUE M A,JAWAID M,et al.Development of kenaf-glass reinforced unsaturated polyester hybrid composite for structural applications[J].Composites Part B: Engineering,2014,56:68-73.
[9]SEVKAT E,LIAW B,DELALE F.Drop-weight impact response of hybrid composites impacted by impactor of various geometries[J].Materials&Design,2013,52:67-77.
[10]CAMANHO P P,MATTHEWS F L.A progressive damage model for mechanically fastened joints in composite laminates[J].Journal of Composite Materials,1999,33(24): 2248-2280.
[11]TAN S C.A progressive failure model for composite laminates containing openings[J].Journal of Composite Materials,1991,25(5):556-577.
[12]HASHIN Z.Failure criteria for unidirectional fiber composites[J].Journal of Applied Mechanics,1980,47(2): 329-334.
[13]BüLTEN B M,KARAKUZU R.Progressive failure analysis of pin-loaded carbon-epoxy woven composite plates[J].Composites Science and Technology,2002,62(9):1259-1271.
[14]ZHOU Song,WANG Zhenqing,ZHOU Jiansheng,et al.Experimental and numerical investigation on bolted composite joint made by vacuum assisted resin injection[J].Composites Part B:Engineering,2013,45(1):1620-1628.
[15]NAIK N K,SAEKHER Y C,MEDURI S.Polymer matrix woven fabric composites subjected to low velocity impact: Part I-damage initiation studies[J].Journal of Reinforced Plastics and Composites,2000,19(12):912-954.