玄碳混雜增強樹脂基復(fù)合材料混雜比優(yōu)化設(shè)計

徐 虹， 盧 巖， 張 可， 劉 栓， 王祥輝， 竇艷麗

(吉林大學 材料科學與工程學院，長春130025)

?

玄碳混雜增強樹脂基復(fù)合材料混雜比優(yōu)化設(shè)計

徐 虹， 盧 巖， 張 可， 劉 栓， 王祥輝， 竇艷麗

(吉林大學 材料科學與工程學院，長春130025)

對碳纖維-玄武巖纖維混雜增強樹脂基復(fù)合材料最優(yōu)混雜比范圍進行研究。以碳纖維與玄武巖纖維平紋織物為增強體，制備9種具有不同混雜比的混雜纖維復(fù)合材料(Hybrid Fiber Reinforced Polymer，HFRP)試樣，并進行拉伸實驗。依據(jù)平紋織物結(jié)構(gòu)特點，計算得出平紋織物單胞性能參數(shù)，在ANSYS中，以SHELL181殼單元體建立HFRP有限元模型。該模型對試樣剛度的模擬值與實驗值近似。分析模型受力時的應(yīng)力云圖發(fā)現(xiàn)，存在將HFRP破壞形式分為一次破壞與二次破壞的臨界混雜比。有限元模擬研究樹脂含量為45%時，不同混雜比的HFRP剛度、強度和拉伸極限應(yīng)變。當混雜比為60%時，可保證HFRP強度無折減的情況下，較玄武巖復(fù)合材料(Basalt Fiber Reinforced Polymer，BFRP)剛度提高93.4%，較碳纖維復(fù)合材料(Carbon Fiber Reinforced Polymer，CFRP)拉伸極限應(yīng)變提高11.3%。

混雜纖維復(fù)合材料；碳纖維；玄武巖纖維；平紋織物；混雜比

輕質(zhì)高強型復(fù)合材料是航空航天、汽車、船舶、化工、電子和建筑等領(lǐng)域都極其需要的。纖維增強樹脂基復(fù)合材料是其中應(yīng)用范圍最廣、最具研究意義的一種。美國波音公司、歐洲空客集團已大力開展CFRP在民用飛機上的應(yīng)用，波音787的CFRP用量高達50%，波音B767油耗降低17%[1-3]。

混雜纖維復(fù)合材料(HFRP)最早被日本的Hayashi[4]提出，認為其可綜合兩種、甚至多種纖維性能上的優(yōu)點。有關(guān)HFRP最為典型的研究為摻雜玻璃纖維來改善碳纖維的斷裂韌性。但近期發(fā)現(xiàn)，玻璃纖維對人體危害極大，不宜再大量生產(chǎn)應(yīng)用。玄武巖纖維力學性能優(yōu)良，適宜與碳纖維混雜，但相關(guān)研究較少[5-6]。

近年來，國內(nèi)外學者對HFRP力學性能的改進主要通過對混雜比與鋪層結(jié)構(gòu)進行設(shè)計。為得出HFRP的強度與剛度，除測試方法外，還有很多學者依靠混合定律、斷裂力學模型、微觀力學模型、隨機臨界核模型等來計算。除此之外，有限元模擬可以更加準確的預(yù)測材料的力學性能，減少實驗工作量[7-13]。

目前，對于HFRP的研究多集中于以單向纖維為增強體的復(fù)合材料，以平紋織物為增強體的HFRP剛度與強度預(yù)測模型以及有限元模擬還很少見。事實上，單向纖維復(fù)合材料易于發(fā)生纖維間的劈裂，不能實際投入生產(chǎn)應(yīng)用，平紋織物復(fù)合材料才更具優(yōu)良的力學性能，大量應(yīng)用于高強度結(jié)構(gòu)制品中[14-17]。

本工作以力學性能優(yōu)良的碳纖維與玄武巖纖維平紋織物為增強體，制備碳纖維-玄武巖纖維層間混雜增強樹脂基復(fù)合材料。在ANSYS中，通過計算平紋織物的單胞參數(shù)，提供一種較為簡便而準確的以平紋織物為增強體的HFRP建模方式。將模型與平紋織物復(fù)合材料細觀結(jié)構(gòu)模型比較，將其模擬所得值與實驗值比較；憑借應(yīng)力云圖明確混雜比對HFRP拉伸時的破壞機制的影響，最終給出HFRP的最優(yōu)混雜比范圍。

1 實驗設(shè)計

1.1材料

碳纖維平紋織物，吉林豐帆特種纖維有限公司；玄武巖纖維平紋織物，浙江石金玄武巖纖維有限公司；鳳凰牌環(huán)氧樹脂E51，南通星辰合成材料有限公司；2-乙基-4-甲基咪唑，天津市化學試劑廠。碳纖維及玄武巖纖維編織參數(shù)見表1。

表1 碳纖維及玄武巖纖維編織參數(shù)Table 1 Weaving parameters of carbon fiber and basalt fiber

1.2試件制備與拉伸實驗

按照GB/T 3354—2014定向纖維增強聚合物基復(fù)合材料拉伸性能試驗方法制備9組具有不同混雜比的碳纖維-玄武巖纖維層間混雜增強樹脂基復(fù)合材料層合板，并且以織物經(jīng)向為試樣長度方向進行切割，每組5件試件，均采用對稱嵌層式鋪層結(jié)構(gòu)(見表2)。試件采用XLB50X2型平板硫化機熱壓法成型、WDW-200型電子萬能試驗機進行拉伸實驗。試樣沿長度方向拉伸，在試樣居中位置安裝50 mm標距引伸計，實驗加載速率控制為3 mm/min，直至試件斷裂。

表2 碳纖維-玄武巖纖維混雜增強樹脂基復(fù)合材料分組Table 2 Groups of carbon/basalt hybrid fiber reinforced resin matrix composites

2 HFRP有限元模型

2.1平紋織物細觀結(jié)構(gòu)分析

以往的單向纖維增強復(fù)合材料有限元建模方法是以ANSYS中專為復(fù)合材料所設(shè)計的殼單元體劃分實體模型，輸入樹脂與纖維紗線性能參數(shù)并設(shè)定層厚、層數(shù)以及鋪層角度,這顯然不適用于平紋織物復(fù)合材料。平紋織物包含經(jīng)、緯兩向紗線，在受拉方向上，經(jīng)向紗線承受主要載荷，緯向紗線起穩(wěn)固連接作用。

Vandeurzen等[18]提出的紗線軸線為正弦曲線，截面為凸透鏡的雙凸模型，可以非常直觀地表現(xiàn)出平紋織物的細觀結(jié)構(gòu)，圖1所示為平紋織物的單胞形態(tài)。

圖1 平紋織物的單胞模型Fig.1 Unit cell of plain weave

如果將一層平紋織物簡單化為殼單元體中經(jīng)、緯兩層，并賦予不同性能參數(shù)，則會喪失織物原本結(jié)構(gòu)特征，大大降低模擬精準度。因為這種建模方式將原本起穩(wěn)固連接作用的緯向纖維全部抽離出來，轉(zhuǎn)化為承載部分。緯向纖維力學參數(shù)在受拉方向上低于經(jīng)向纖維兩個數(shù)量級， HFRP在受拉過程中各部分等位移延伸，這將導致模擬過程中材料過早失效。因此，必須準確計算出平紋織物單胞各方向上的性能參數(shù)，將其輸入為殼單元體平紋織物層的材料屬性，完成有限元建模。

2.2模型建立

根據(jù)混合定律，以經(jīng)、緯紗線不同方向的彈性模量乘以相應(yīng)體積含量求得平紋織物單胞各方向上的彈性模量：

(1)

式中：表示彈性模量；表示纖維含量，可根據(jù)織物編制參數(shù)求得；下角標c，b分別表示碳纖維和玄武巖纖維;i,j分別表示平紋織物經(jīng)向方向與緯向方向；l,t分別表示纖維紗線軸線方向和徑向方向。其他性能參數(shù)，如剪切模量G、泊松比可以用同樣方法求得。

確定平紋織物單胞參數(shù)后，創(chuàng)建有限元模型。根據(jù)模型中的纖維與實際纖維的體積相等，得到SHELL單元中碳纖維的片層厚度和玄武巖纖維的片層厚度：

(2)

式中：δ表示織物的面密度；ρ表示纖維的密度。樹脂層厚度Hm為：

(3)

在ANSYS中選定四節(jié)點三維殼單元SHELL181定義單元類型，按照上述計算，賦予纖維體正交各向異性屬性，樹脂為各向同性屬性，按表2設(shè)置層數(shù)與層厚，鋪層角為0°。在GUI界面中完成試件實體模型建立并劃分網(wǎng)格，至此完成HFRP有限元模型的建立。

2.3模型結(jié)構(gòu)精準度檢測

事實上，以殼單元體建立有限元模型仍然簡化了平紋織物復(fù)合材料的真實結(jié)構(gòu)。目前，有學者提出了單層平紋織物復(fù)合材料細觀模型的建立方法。但HFRP結(jié)構(gòu)較為復(fù)雜，其中包含鋪層數(shù)目，鋪層角度，混雜比例等多個變量。如考慮細觀模型建立方法則需建立多層復(fù)合材料，為保證正弦曲線結(jié)構(gòu)不至于失真，當中所創(chuàng)建的節(jié)點數(shù)目之多，將使得計算機運算極為困難，運算周期較長。

通過圖1可以看出，平紋織物復(fù)合材料是極為扁平的結(jié)構(gòu)，其內(nèi)部極微小的形狀波動對最終模擬結(jié)果影響不大。為驗證建立的模型對復(fù)合材料結(jié)構(gòu)的還原程度，建立了相同尺寸的單層碳纖維平紋織物復(fù)合材料的細觀結(jié)構(gòu)模型與殼單元體模型，細觀結(jié)構(gòu)模型如圖2所示。分別固定模型單側(cè)端面，對另一側(cè)端面進行等伸長拉伸，施加位移載荷。通過對節(jié)點反力的面內(nèi)積分，求出受力面的等效拉伸應(yīng)力。結(jié)果顯示,細觀結(jié)構(gòu)模型等效拉伸應(yīng)力為123.302 MPa，根據(jù)E=δ/ε求得剛度為61.651 GPa。殼單元體模型等效拉伸應(yīng)力為121.898 MPa，換算為剛度為60.949 GPa。證明本工作建立的模型適用于一般情況下對以平紋織物為增強體的HFRP力學性能的預(yù)測。這保證了模擬的精準度，又大大減輕了計算負擔。

圖2 細觀結(jié)構(gòu)模型Fig.2 Model of mesostructure

2.4模型數(shù)值精準度檢測

為進一步檢測本工作所建模型的準確性，將有限元模擬值與實驗值進行了比對，對模擬精準度進行檢測。按照表2所示信息對試件進行有限元建模，同樣方式施加位移載荷。根據(jù)E=σ/ε求得復(fù)合材料的拉伸剛度。HFRP剛度的實驗值與模擬值對比見圖3。

圖3 碳纖維-玄武巖纖維混雜增強樹脂基復(fù)合材料彈性模量實驗值與模擬值對比Fig.3 Comparison of experimental values and simulated values of elasticity modulus for carbon/basalt hybrid fiber reinforced resin matrix composites

由圖3可以看出，實驗值與模擬值基本近似，實驗所測得數(shù)值略小于有限元模擬數(shù)值是因為試件制備難以避免孔隙等缺陷。

3 HFRP層合板破壞機制分析

3.1HFRP的強度分析

為確定材料的最優(yōu)混雜比范圍，必須明確材料的破壞機制。實驗所測HFRP強度如圖4所示，碳纖維含量低于25%左右時，HFRP的抗拉強度隨著碳纖維含量的增加而減小，碳纖維含量高于25%左右時，趨勢相反。強度值在碳纖維含量為25%左右時出現(xiàn)了最小值，此值被稱為臨界混雜比。產(chǎn)生此現(xiàn)象的原因是，碳纖維的伸長率小，在HFRP受拉時，其先承受大部分載荷并斷裂。這時，若玄武巖纖維含量較多，碳纖維斷裂后，其仍具備一定承載能力，HFRP可以繼續(xù)承載，直至玄武巖纖維失效HFRP才徹底失效。最終HFRP發(fā)生二次破壞，強度僅為玄武巖纖維這部分所能承擔的極限強度，如式(4)，其中Vc表示混雜比。因此，隨著碳纖維含量的增加，HFRP的抗拉強度隨之降低。若碳纖維斷裂時，玄武巖纖維的含量較少，其不能單獨承載，便與碳纖維一起斷裂失效，最終HFRP為一次破壞，強度為碳纖維與玄武巖纖維共同承載可經(jīng)受的極限強度，如式(5)。因此，隨著碳纖維含量的增加，HFRP強度也隨之增加。

σH=(1-Vc)σb

(4)

(5)

圖4 碳纖維-玄武巖纖維混雜增強樹脂基復(fù)合材料抗拉強度與混雜比的關(guān)系Fig.4 Relationship between tensile strength and hybrid ratio of carbon/basalt hybrid fiber reinforced resin matrix composites

3.2HFRP層合板破壞機制分析

為驗證以上結(jié)論，基于ANSYS對HFRP進行有限元分析。實際上，以平紋織物為增強體的復(fù)合材料中，纖維層的抗拉強度為經(jīng)向紗線軸向強度與緯向紗線徑向強度的疊加，但紗線徑向強度要比軸向強度小兩個數(shù)量級，所以可以忽略緯向紗線的作用。纖維層的抗拉強度可以近似考慮為：

σ≈λi×El

(6)

已知碳纖維紗線抗拉強度為2600 MPa，玄武巖纖維紗線抗拉強度為2200 MPa，根據(jù)式(6)可以得到HFRP中碳纖維層和玄武巖纖維層的抗拉強度分別為1266 MPa和905 MPa。

圖5顯示臨界混雜比發(fā)生在試樣5(混雜比26.94%左右)。實驗測得試樣5最大可經(jīng)受335.31 MPa的應(yīng)力，在其有限元模型上一端施加全約束，一端施加最大應(yīng)力可得到圖5(a)所示的應(yīng)力云圖?？梢娞祭w維所受應(yīng)力剛好達到其抗拉強度，將碳纖維層退化，使其作用與樹脂一樣，可以得到如圖5(b)的應(yīng)力云圖，此時玄武巖纖維所受應(yīng)力也超過了其抗拉強度?？梢娫谔祭w維破壞之后玄武巖纖維也無法繼續(xù)承載，HFRP發(fā)生一次斷裂。

圖5 試樣5應(yīng)力云圖 (a)碳纖維未失效；(b)碳纖維失效Fig.5 Stress nephogram of P05 (a)carbon fiber does not lose efficacy;(b) carbon fiber loses efficacy

再來考察試樣4(見圖6)。實驗測得試樣4最大可經(jīng)受387.19 MPa的應(yīng)力。對其加載可以得到圖6(a)所示的應(yīng)力云圖?？梢娞祭w維層的應(yīng)力已達到1280 MPa，超過了抗拉極限，此時碳纖維應(yīng)已失效，將碳纖維層退化可以得到如圖6(b)所示的應(yīng)力云圖，此時玄武巖層所受應(yīng)力剛好達到抗拉強度。顯然，試樣4斷裂時是碳纖維先發(fā)生斷裂，載荷繼續(xù)由玄武巖纖維承擔，最終HFRP發(fā)生二次破壞。

圖6 試樣4應(yīng)力云圖 (a)碳纖維未失效;(b)碳纖維失效Fig.6 Stress nephogram of P04 (a)carbon fiber does not lose efficacy;(b)carbon fiber loses efficacy

4 混雜比的優(yōu)化設(shè)計

根據(jù)以上分析，確實存在一個臨界混雜比值，HFRP的強度與極限拉伸應(yīng)變必將會在此值處存在轉(zhuǎn)折點。此臨界混雜比區(qū)分HFRP層合板發(fā)生一次破壞還是二次破壞。二次破壞是實際應(yīng)用中不被允許的，最優(yōu)混雜比范圍不得將其包含其中。

ANSYS模擬可以排除孔隙率，制備工藝，樹脂含量不同等影響因素所帶來的實驗誤差，在理論上給出準確結(jié)論。因此，模擬得出HFRP的剛度、強度、拉伸極限應(yīng)變與混雜比的關(guān)系曲線。模擬時將樹脂含量定為45%。

按2.3節(jié)所述即可對剛度進行模擬。強度模擬則按照Hashin破壞準則通過UPFs自行編譯,用戶自定義失效準則。在ANSYS中不斷對HFRP增加應(yīng)力載荷直至纖維層全部失效，即輸出強度模擬結(jié)果，同時可獲得材料的拉伸極限應(yīng)變。最終模擬曲線見圖7。

圖7 碳纖維與玄武巖纖維混雜增強樹脂基復(fù)合材料樹脂含量為45%時力學性能模擬曲線 (a)彈性模量與混雜比的關(guān)系；(b)抗拉強度與混雜比的關(guān)系；(c)拉伸極限應(yīng)變與混雜比的關(guān)系Fig.7 Simulation curve of mechanical property for carbon/basalt hybrid fiber reinforced resin matrix composites with 45% resin content (a) relationship between elasticity modulus and hybrid ratio；(b)relationship between tensile strength and hybrid ratio；(c)relationship between limit tensile strain and hybrid ratio

由圖7可見，在排除HFRP發(fā)生二次破壞范圍(0%

5 結(jié) 論

(1)對以平紋織物為增強體的HFRP的有限元建模方式可以準確預(yù)測HFRP的各項性能。碳纖維與玄武巖纖維混雜增強樹脂基復(fù)合材料剛度隨混雜比(碳纖維含量)增加而增加，符合混合定律；強度在混雜比為25%左右時出現(xiàn)最小值，這種現(xiàn)象與HFRP破壞機制有關(guān)。提供的ANSYS分析方法可以驗證不同混雜比對HFRP的破壞機制的影響。

(2)當混雜比大于60%且不等于100%時，可保證HFRP強度、剛度、拉伸極限應(yīng)變均有提高。當混雜比為60%時可保證HFRP強度無折減的情況下，剛度較玄武巖纖維復(fù)合材料提高93.4%，拉伸極限應(yīng)變較碳纖維復(fù)合材料提高11.3%。

[1] VALENZA A，F(xiàn)IORE V，DI BELLA G．Effect of UD carbon on the specific mechanical properties of glass mat composites for marine applications[J]．Journal of Composite Materials，2010，44(11)：1350-1364．

[2] SUN X S，CHEN Y，TAN V B C，etal．Homogenization and stress analysis of mutilayered composite offshore production risers[J]．Journal of Applied Mechanics，2014，81(3)：1-11．

[3] KOTHA S P，CUZELSU N．Tensile behavior of corti-cal bone：dependence of organic matrix material properties on bone mineral content[J]．Journal of Biomechanics，2007，40(1)：35-45．

[4] HAYASHI T．On the improvement of mechanical properties of composites by hybrid composition [C]∥ Proceedings of the 8th International Reinforced Plastics Conference.Brighten,UK:[s.n.],1972．

[5] ZHANG Y，LI Y，MA H，etal．Tensile and interfacial properties of unidirectional flax/glass fiber reinforced hybrid composites[J]．Composites Science and Technology，2013，88(14)：172-177．

[6] 張普，朱虹，孟少平，等．CFRP/BFRP混雜效應(yīng)試驗研究和理論分析[J]．工業(yè)建筑，2011，41(3)：103-107．

(ZHANG P，ZHU H，MENG S P，etal．Theoretical analysis and test research on hybrid effect of CFRP/BFRP hybrid fiber[J]．Industrial Construction，2011，41(3)：103-107．)

[7] SWOLFS Y，GORBATIKH L，VERPOEST I．Fiber hybridisation in polymer composites：a review[J]．Composites Part A：Applied Science and Manufacturing，2014，67(12)：181-200．

[8] 張小東，王耀杰，鄧宗才．混雜纖維片材拉伸性能研究[J]．國防交通工程與技術(shù)，2010，8(1)：21-25. (ZHANG X D，WANG Y J，DENG Z C．An experimental study of the tensile properties of IHFRP[J]．Traffic Engineering and Technology for National Defence，2010，8(1)：21-25．)

[9] 范賦群，曾慶敦，黃小清．復(fù)合材料的統(tǒng)計斷裂理論[J]．力學與實踐，2000,22(2)：1-5．

(FAN F Q，ZENG Q D，HUANG X Q．Probabilistic fracture theories for composites[J]．Mechanics in Engineering，2000，22(2)：1-5．)

[10] 李衛(wèi)東，曹海琳，魏斌，等．BF/CF層間混雜結(jié)構(gòu)對復(fù)合材料性能影響[J]．熱固性樹脂，2009，24(5)：39-43.

(LI W D，CAO H L，WEI B，etal．Effect of CF/BF interlaminated hybrid structure on composite properties[J]．Thermosetting Resin，2009，24(5)：39-43．)

[11] QU P，GUAN X J．JIA Y X，etal．Effective elastic properties and stress distribution of 2D biaxial nonorthogonally braided composites[J]．Journal of Composite Materials，2012，46(8)：997-1008．

[12] AYRANCI C and CAREY J P．Predicting the longitudenal elastic modulus of braided tubular composites using a curbed unit-cell geometry[J]．Composites Part B：Engineering，2010，41(3)：229-235．

[13] NASU S，OHTANI A，NAKAI A，etal．Deformation behavior and mechanical properties of braided rectangular pipes[J]．Composite Structures 2010：92(3)：752-756．

[14] PAPANICOLAOU G C，MOUZAKIS D E，KOSMIDOU T V．Static and dynamic behavior of single-edge notched glass fabric composites[J]．Polymer Composites．2006．27(2)：177-183．

[15] 魏廣恒．基于變分原理的機織復(fù)合材料細觀力學模型[D]．天津：天津工業(yè)大學，2006．

(WEI G H．The micro-mechanics model of the woven fabric composite based on variational principle[D]．Tianjin：Tianjin Ploytechnic University，2006．)

[16] SHROTRIYA P，SOTTOS N R．Viscoelastic response of woven composite substrates[J]．Composites Science and Technology，2005，65(3/4)：621-634．

[17] 左中鵝，王瑞，徐磊．基于有限單元法的平紋織物復(fù)合材料強度預(yù)測：1．RVE的有限元模型[J]．紡織學報，2009，30(12)：45-49．

(ZUO Z E，WANG R，XU L．Mechanical strength prediction of plain woven fabric composite 1．finite element model of composite RVE[J]．Journal of Textile Research，2009，30(12)：45-49．)

[18] VANDEURZEN P，IVENS J，VERPOEST I．Structure performance analysis of two-dimensional woven fabric composites[J]．Polymers & Polymer Composites，1996，4(5)：361-367．

(責任編輯：徐永祥)

OptimalDesignforHybridRatioofCarbon/BasaltHybridFiberReinforcedResinMatrixComposites

XU Hong， LU Yan， ZHANG Ke， LIU Shuan， WANG Xianghui， DOU Yanli

(School of Materials Science and Engineering, Jilin University, Changchun 130025, China)

The optimum hybrid ratio range of carbon/basalt hybrid fiber reinforced resin composites was studied. Hybrid fiber composites with nine different hybrid ratios were prepared before tensile test.According to the structural features of plain weave, the unit cell’s performance parameters were calculated. Finite element model was established by using SHELL181 in ANSYS. The simulated values of the sample stiffness in the model were approximately similar to the experimental ones. The stress nephogram shows that there is a critical hybrid ratio which divides the failure mechanism of HFRP into single failure state and multiple failure state. The tensile modulus, strength and limit tensile strain of HFRP with 45% resin are simulated by finite element method. The result shows that the tensile modulus of HFRP with 60% hybrid ratio increases by 93.4% compared with basalt fiber composites (BFRP), and the limit tensile strain increases by 11.3% compared with carbon fiber composites(CFRP).

hybrid fiber composites； carbon fiber； basalt fiber；plain weave； hybrid ratio

2017-04-10；

2017-05-24

吉林省“雙十”重大科技成果轉(zhuǎn)化項目(3D515AI92416)；吉林大學大學生創(chuàng)新創(chuàng)業(yè)訓練計劃(2016A43127)

竇艷麗(1979—)，女，博士，副教授，研究方向為混雜纖維制備與性能，(E-mail)douyl@jlu.edu.cn。

10.11868/j.issn.1005-5053.2017.000047

TB332

： A

： 1005-5053(2017)04-0007-07