胡 妞，李布楠，高本征

（清華大學(xué) 深圳研究生院 新材料研究所， 廣東 深圳 518055）

炭纖維/樹脂復(fù)合材料導(dǎo)熱性能的數(shù)值模擬

胡 妞，李布楠，高本征

（清華大學(xué) 深圳研究生院 新材料研究所， 廣東 深圳 518055）

采用有限元方法對炭纖維/樹脂復(fù)合材料的導(dǎo)熱性能進(jìn)行了數(shù)值模擬，分別建立一維結(jié)構(gòu)和二維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料計(jì)算分析模型，研究炭纖維含量、界面接觸熱阻、以及炭纖維直徑對復(fù)合材料有效熱導(dǎo)率的影響。研究結(jié)果表明炭纖維作為復(fù)合材料增強(qiáng)相，其含量越高復(fù)合材料的熱導(dǎo)率越高；界面的接觸熱阻在10-3～10-5(m2K)/W范圍內(nèi)對復(fù)合材料有效熱導(dǎo)率有較大的影響，超出范圍之后改變接觸熱阻對材料熱導(dǎo)率的影響可以忽略；接觸熱阻比較大時，炭纖維的直徑對復(fù)合材料的熱導(dǎo)率有較大的的影響，當(dāng)接觸熱阻比較小時，炭纖維的直徑對于復(fù)合材料熱導(dǎo)率的影響非常小。

有限元；炭纖維/樹脂復(fù)合材料；熱導(dǎo)率

隨著化工、能源、電子信息、電氣工程、航空航天領(lǐng)域的發(fā)展[1,2]，對材料的導(dǎo)熱性能、耐腐蝕性能、抗氧化性能以及力學(xué)性能等提出了更高的要求。由于高分子樹脂材料具有優(yōu)良的耐腐蝕性能和優(yōu)異的力學(xué)性能，但大部分的高分子材料導(dǎo)熱性能比較差。中間相瀝青基炭纖維具有高度晶體取向的有序結(jié)構(gòu)而呈現(xiàn)出高導(dǎo)熱及輕質(zhì)高強(qiáng)度等特點(diǎn)，中間相瀝青基炭纖維的熱導(dǎo)率可以高達(dá)1 000 W/(m?K)[3]，但其很難單獨(dú)應(yīng)用于工程領(lǐng)域，因此其常作為導(dǎo)熱增強(qiáng)相用于制備 C/C復(fù)合材料或炭/樹脂復(fù)合材料[4-7]。目前炭纖維/樹脂復(fù)合材料的研究和應(yīng)用比較廣泛，炭纖維復(fù)合材料通過調(diào)整炭纖維直徑、含量、復(fù)合材料界面的接觸熱阻，可以對其熱物理性能和力學(xué)性能進(jìn)行設(shè)計(jì)，從而滿足工程領(lǐng)域的應(yīng)用要求。對于單相材料，其熱導(dǎo)率與材料的晶體結(jié)構(gòu)和純度密切相關(guān)。而復(fù)合材料的熱導(dǎo)率取決于各組元的導(dǎo)熱能力[8]、體積分?jǐn)?shù)[9]、增強(qiáng)相的尺寸[10]以及界面接觸熱阻[11]等因素。

有限元法的基本思想是將求解區(qū)域離散為有限個按一定的規(guī)則相互連接的單元，然后組合為系統(tǒng)方程組求解，使用計(jì)算機(jī)可提高分析的效率。在計(jì)算機(jī)上數(shù)值模擬固體復(fù)合材料有效熱導(dǎo)率是一種高效的手段，其具有突出的特點(diǎn)，可以模擬各種增強(qiáng)相狀態(tài)下的有效熱導(dǎo)率。Jose[12]運(yùn)用Mathematica軟件模擬計(jì)算二維結(jié)構(gòu)和三維結(jié)構(gòu)復(fù)合材料的有效熱導(dǎo)率，并研究了填充料的形狀和分布對于材料熱導(dǎo)率的影響。D. Marcos-Gómez[13]利用有限元離散化方法研究界面接觸熱阻和填料排布對復(fù)合材料有效熱導(dǎo)率的影響，其研究表明當(dāng)接觸熱阻較大時，增強(qiáng)相的排布對復(fù)合材料有效熱阻有較大的影響，當(dāng)接觸熱阻較小時，增強(qiáng)相的排布對復(fù)合材料熱導(dǎo)率的影響非常小。Siddiqui, MOR[14]采用有限元法建立單元體模型研究織物材料的有效熱導(dǎo)率，其通過模擬跟實(shí)驗(yàn)的對比得到材料的界面接觸熱阻，研究了纖維含量和其熱導(dǎo)性能對織物有效熱導(dǎo)率的影響。Chmielewski, M[15]采用有限元方法去分析模擬研究Cu-AlN/陶瓷復(fù)合材料中孔隙率對復(fù)合材料導(dǎo)熱性能的影響。本文采用有限元方法對炭纖維/樹脂復(fù)合材料的導(dǎo)熱性能進(jìn)行了數(shù)值模擬，分別建立一維結(jié)構(gòu)和二維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料計(jì)算分析模型，研究炭纖維含量、界面接觸熱阻、以及炭纖維直徑對復(fù)合材料有效熱導(dǎo)率的影響。

1 材料及性能

采用熱壓成型制備炭纖維/樹脂復(fù)合材料，以中間相瀝青纖維作為復(fù)合材料的增強(qiáng)材料，環(huán)氧樹脂作為基體材料。首先將環(huán)氧樹脂樹脂溶入適量的丙酮中，然后將石墨化的炭纖維按照一定的要求排布在磨具內(nèi)，然后再把環(huán)氧樹脂丙酮溶液倒入磨具，待丙酮揮發(fā)后，將磨具置入熱壓機(jī)中熱壓，然后制的一維排布炭纖維/樹脂復(fù)合材料和二維炭纖維/樹脂復(fù)合材料。一維結(jié)構(gòu)炭纖維復(fù)合材料炭纖維沿一個方向定向排布，二維結(jié)構(gòu)炭纖維復(fù)合材料是炭纖維交替相互垂直排布。圖1和圖2分別是炭纖維體積分?jǐn)?shù)為30%一維結(jié)構(gòu)和二維結(jié)構(gòu)復(fù)合材料軸向截面和徑向截面的偏光圖（PLM）。從圖中可以看出炭纖維排布比較規(guī)整，纖維與樹脂之間結(jié)合緊密。

有限元中所采用的炭纖維的軸向熱導(dǎo)率約為780 W/(m K)，炭纖維徑向熱導(dǎo)率為25 W/(m K)，炭纖維的密度約為1.86 g/cm3，比熱為700 J/(kg K)，環(huán)氧樹脂的密度為1.1 g/cm3，比熱為1 300 J/(kg K)，界面接觸熱阻約為10-5(m2K)/W。采用LFA447熱導(dǎo)儀對部分的復(fù)合材料熱導(dǎo)率進(jìn)行測試，通過對比數(shù)值計(jì)算得到的背面升溫曲線和 LFA測試的紅外曲線確定數(shù)值計(jì)算模型。

圖1 一維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料的偏光圖Fig.1 Polarized light microscopy photographs of the one dimension carbon fiber reinforced composites

2 有限元模擬

激光閃光法（LFA）是一種有效測試熱導(dǎo)率的瞬態(tài)測試方法，激光閃光法的基本原理在一個絕熱薄試樣的正面照射一個垂直于試樣正面的均勻脈沖激光。試樣背面的溫度隨時間的變化可用方程（1）表示[16]：

紅外探測器探測到背面的溫升變化轉(zhuǎn)化成溫度變化的信號曲線，通過對比數(shù)值計(jì)算得到試樣背面的溫升曲線和實(shí)驗(yàn)紅外探測器探測到的溫升信號曲線確定數(shù)值計(jì)算模型，然后對模型施加穩(wěn)態(tài)的邊界條件，根據(jù)穩(wěn)態(tài)模型研究復(fù)合材料微觀性能對復(fù)合材料有效熱導(dǎo)率的影響。

圖2 二維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料的偏光圖Fig.2 Polarized light microscopy photographs of the two dimension carbon fiber reinforced composites

從樣品的 PLM圖中可以看出炭纖維分布比較均勻，但是纖維有出現(xiàn)了劈裂，為了簡化數(shù)值計(jì)算，將炭纖維簡化成圓柱狀，由于樣品四周絕熱，且樣品具有二維周期性，則一維和二維的數(shù)值計(jì)算模型如圖3和圖4所示（注：從PLM圖，二維結(jié)構(gòu)的炭纖維復(fù)合材料中炭纖維不是嚴(yán)格按照模型單層交替排布，而是多層交替排布，選擇炭纖維20層橫向排布，然后再20層縱向垂直排布計(jì)算得到的背面溫升曲線與實(shí)驗(yàn)測得背面溫升曲線有較好的一致性，其計(jì)算單元體如圖4所示。）

數(shù)值計(jì)算模型圖根據(jù)實(shí)驗(yàn)測得的碳纖維幾何尺寸和碳纖維體積百分比確定，本次研究中所采用的是直徑約為30μm的中間相瀝青基石墨化纖維。根據(jù)圖3、4所示的穩(wěn)態(tài)的邊界條件，則一維、二維炭纖維/樹脂復(fù)合材料中的熱傳導(dǎo)方程如方程（2）、(3)所示。

圖3 一維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料計(jì)算模型Fig.3 The numerical simulation model of one dimension carbon-fiber/resin composite

圖4 二維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料計(jì)算模型Fig.4 The numerical simulation model of two dimension carbon-fiber/resin composite

式中:△T=Thot-Tcold；L表示樣品的厚度；q是任意垂直截面的平均熱流密度。

3 結(jié)果與分析

3.1 炭纖維的含量對復(fù)合材料導(dǎo)熱性能的影響

圖5 復(fù)合材料有效熱導(dǎo)率隨纖維質(zhì)量分?jǐn)?shù)變化曲線Fig.5 The changes of effective thermal conductivities of composites with mass fractions of fibers

通過數(shù)值計(jì)算得到不同溫度條件下，其他參數(shù)不變，改變復(fù)合材料中的炭纖維質(zhì)量分?jǐn)?shù)，復(fù)合材料的有效熱導(dǎo)率。圖5（a）所示的是一維結(jié)構(gòu)下復(fù)合材料有效熱導(dǎo)率隨炭纖維含量的變化曲線，從圖中可以看出在不同的溫度條件下復(fù)合材料的有效熱導(dǎo)率都隨炭纖維含量的增加而增加，圖5（b）所示的是二維結(jié)構(gòu)下復(fù)合材料有效熱導(dǎo)率隨炭纖維含量的變化曲線，該結(jié)構(gòu)的炭纖維復(fù)合材料的有效熱導(dǎo)率也隨炭纖維含量的增加而增加。一維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料的有效熱導(dǎo)率隨炭纖維含量增加而增加的速率大于二維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料。炭纖維/樹脂基復(fù)合材料有效熱導(dǎo)率隨著炭纖維含量的增加而增加主要是因?yàn)樘坷w維的軸向熱導(dǎo)率遠(yuǎn)遠(yuǎn)的高于樹脂的熱導(dǎo)率，在熱傳遞過程中，從炭纖維上傳遞的熱流量遠(yuǎn)遠(yuǎn)的大于樹脂上的熱流量，隨著炭纖維含量的增加，單位截面內(nèi)炭纖維的截面積增大，傳遞熱量的能力增大，材料的導(dǎo)熱性能增加。

3.2 接觸熱阻對復(fù)合材料導(dǎo)熱性能的影響

圖6 接觸熱阻與復(fù)合材料熱導(dǎo)率之間的關(guān)系圖Fig.6 The relationship between the thermal conductivity of composites and the interface thermal resistance.

采用數(shù)值模擬計(jì)算出來的接觸熱阻與復(fù)合材料導(dǎo)熱性能之間的關(guān)系曲線如圖6所示。圖6(a)所示的是采用數(shù)值模擬計(jì)算出的不同界面接觸熱阻條件下一維炭纖維/樹脂基復(fù)合材料的半升溫時間 t1/2（注：一維碳纖維/樹脂復(fù)合材料穩(wěn)態(tài)傳熱過程中熱量沒有穿過界面接觸處，因?yàn)闉榱吮碚鹘缑鏌嶙璧膶υ摻Y(jié)構(gòu)復(fù)合材料熱傳遞性能的影響就采用了瞬態(tài)傳熱過程研究），半升溫時間與材料的熱擴(kuò)散系數(shù)成反比，在一定程度上反映材料傳遞溫度的能力，半升溫時間越短代表材料傳遞溫度的能力越強(qiáng)，其熱擴(kuò)散系數(shù)越大。從圖 6(a)中可以看出，當(dāng)接觸熱阻從10-1(m2K)/W減小到10-3(m2K)/W，和接觸熱阻從10-5(m2K)/W減小到10-8(m2K)/W時，半升溫時間t1/2基本上沒有變化，及材料的熱擴(kuò)散系數(shù)未變化，當(dāng)接觸熱阻從10-3(m2K)/W減小到10-5(m2K)/W，半升溫時間t1/2明顯變小，即表明熱傳導(dǎo)性能增加，這表明接觸熱阻在此范圍，減小接觸熱阻能有效的增加復(fù)合材料的熱導(dǎo)率。圖 6(b)所示的數(shù)值計(jì)算求解的二維結(jié)構(gòu)炭纖維復(fù)合材料的熱導(dǎo)率隨著接觸熱阻變化而變化的圖。從圖中可以看出，界面接觸熱阻在10-3(m2K)/W～10-5(m2K)/W之間，復(fù)合材料的有效熱導(dǎo)率會隨著接觸熱阻的減小而增大，超出此范圍界面接觸熱阻，增大或減小界面的接觸熱阻對復(fù)合材料的有效熱導(dǎo)率的影響非常小。

3.3 炭纖維直徑對復(fù)合材料熱性能影響

圖7 不同接觸熱阻下，復(fù)合材料有效熱導(dǎo)率隨炭纖維直徑變化曲線Fig.7 The changes of effective thermal conductivities of composites with mass diameter of fibers under different interface thermal resistance

圖7是不同接觸熱阻下纖維直徑對二維結(jié)構(gòu)炭纖維/樹脂復(fù)合材料有效熱導(dǎo)率的影響，從圖7中可以看出接觸熱阻比較大時，其他條件不變的情況下，復(fù)合材料的有效熱導(dǎo)率會隨著炭纖維直徑的增加而增加，隨著接觸熱阻的減小，這種變化規(guī)律越不明顯。這種變化的主要原因是，當(dāng)界面接觸熱阻存在的條件下，炭纖維越粗，相對比表面就越少，界面接觸熱阻就相對比較小，復(fù)合材料的有效熱導(dǎo)率就會相應(yīng)的越大，當(dāng)接觸熱阻越大，這種影響效果就越明顯。從圖 7中可以看出當(dāng)界面接觸熱阻為10-1(m2K)/W時，復(fù)合材料有效熱導(dǎo)率隨著炭纖維直徑增大而增大的趨勢明顯，而當(dāng)界面接觸熱阻為10-10(m2K)/W時，改變炭纖維的直徑，復(fù)合材料的有效熱導(dǎo)率基本上沒有變化。該研究表明接觸熱阻比較大時，炭纖維的直徑對復(fù)合材料的熱導(dǎo)率有較大的的影響，當(dāng)接觸熱阻比較小時，炭纖維的直徑對于復(fù)合材料熱導(dǎo)導(dǎo)熱的影響非常小。

4 結(jié) 論

（1）研究表明炭纖樹脂基復(fù)合材料的熱導(dǎo)率隨著炭纖維增強(qiáng)相的增加而增加。

（2）接觸熱阻的存在對復(fù)合材料的熱導(dǎo)率有一定的影響，通過一維結(jié)構(gòu)和二維結(jié)構(gòu)的炭纖維數(shù)值復(fù)合材料的模擬研究表明在一定范圍 10-5(m2K)/W～10-3(m2K)/W內(nèi)，復(fù)合材料的熱導(dǎo)率隨著接觸熱阻的增加而減小，超出該范圍以后接觸熱阻對復(fù)合材料的傳熱性能的影響非常小。

（3）由于接觸熱阻的存在炭纖維直徑對于復(fù)合材料的熱導(dǎo)率有一定的影響，當(dāng)接觸熱阻越大時，復(fù)合材料的熱導(dǎo)率隨著炭纖維的直徑增大而增大的趨勢越明顯。當(dāng)接觸熱阻比較小時，這種影響就減弱，接觸熱阻小到一定程度后，炭纖維直徑變化對接觸熱阻的影響可以忽略不計(jì)。

［1］Ferrer-Anglada N, Gomis V, EI-Hachemi Z, et al. Carbon nanotubes based composites for electronic application: CNT-conducting polymers, CNT-Cu[J]. Physical status Solidi Application and Materials Science, 2006, 203(6): 1082-1087.

［2］儲九榮，張曉輝，徐傳驤．導(dǎo)熱高分子材料的研究與應(yīng)用［J］．高分子材料科學(xué)與工程，2000，16（4）：17-21．

［3］高曉晴, 郭全貴, 劉朗. 高導(dǎo)熱炭材料的研究進(jìn)展[J]. 功能材料, 2006, (02): 173-177.

［4］Bonal J P,Wu C H. Neutron irradiation effects on the thermal conductivity and dimensional stability of carbon fiber composites at divertor conditions [J]. Journal of Nuclear Materials, 1996, 228(2): 155-161.

［5］Lewis J S, Lackey W J, Vaidyaraman S. Model for prediction of matrix microstructure for carbon/carbon composites prepared by forced flow-thermal gradient CVI [J]. Carbon, 1997, 35(1): 103-112.

［6］Chen Y M ,Ting J M. Ultra high thermal conductivity polymer composites [J]. Carbon, 2002, 40(3): 359-362.

［7］Soo-Jin Park, Min-seok Cho, jac-Rock Lee, et al. Influence of molybdenum disilicide filler on carbon carbon composites [J]. Carbon, 1999, 37(11): 1685-1689.

［8］Molina J M，Rheme M，Carron J. Thermal conductivity of aluminum matrix compositess reinforced with mixtures of diamond and SiC particles［J］.Scripta Material，2008，58: 39－396.

［9］龍樂. 電子封裝中的鋁碳化硅及其應(yīng)用［J］.電子與封裝，2006，(6) :16－20.

［10］Molina J M，Narciso J，Weber L. Thermal conductivity of Al-SiC compositess with monomodal and bimodal particle size distribution［J］. Materials Science and Engineering，2008，480:483－488.

［11］Nan C-W, Birringer R, Clarke D R. Effective thermal conductivity of particulate composites with interfacial thermal resistance[J]. Journal of Applied Physics, 1997, 81(10): 6692.

［12］Gómez-Mu?oz J L, Bravo-Castillero J. Calculation of effective conductivity of 2D and 3D composite materials with anisotropic constituents and different inclusion shapes in Mathematica[J]. Computer Physics Communications, 2008, 179(4): 275-287.

［13］Marcos-Gómez D, Ching-Lloyd J, Elizalde M R. Predicting the thermal conductivity of composite materials with imperfect interfaces[J]. Composites Science and Technology, 2010, 70(16): 2276-2283.

［14］Siddiqui M O R, Sun D M. Finite element analysis of thermal conductivity and thermal resistance behaviour of woven fabric[J]. Computational Materials Science, 2013, 75: 45-51.

［15］Chmielewski M, Weglewski W. Comparison of experimental and modelling results of thermal properties in Cu-AlN composite materials[J]. Bulletin of the Polish Academy of Sciences-Technical Sciences, 2013, 61(2): 507-514.

［16］Parker W J, Jenkins R J, Abbott G L, et al. Flash method of determining thermal diffusivity, heat capacity, and thermal conductivity[J]. Journal of Applied Physics, 1961, 32(9): 1679-1684.

硅藻泥搶奪壁紙涂料市場蛋糕

從硅藻土中提取的硅藻泥壁材正以環(huán)保、美化等功能迅速進(jìn)入公眾視線，且謀切著壁紙、涂料市場蛋糕。在近期于四川成都召開的全國硅藻泥產(chǎn)業(yè)發(fā)展研討會上，業(yè)內(nèi)專家預(yù)測，十年后硅藻泥將與壁紙、乳膠漆在壁材市場形成三分天下的格局。 在此次全國硅藻泥產(chǎn)業(yè)發(fā)展研討會上，中國建筑材料聯(lián)合會生態(tài)環(huán)境建材分會成立了“硅藻泥專業(yè)技術(shù)委員會”，以加強(qiáng)對行業(yè)的服務(wù)和引領(lǐng)，助推硅藻泥行業(yè)健康發(fā)展。同時，生態(tài)環(huán)境建材分會為硅藻泥企業(yè)設(shè)立、頒發(fā)了“中國硅藻泥示范單位”等三個獎項(xiàng)，并首發(fā)了硅藻泥行業(yè)第一本專業(yè)技術(shù)書籍《硅藻泥裝飾壁材》，對行業(yè)進(jìn)行科普。

有“硅藻泥壁材創(chuàng)始人”之稱的中國建筑材料聯(lián)合會生態(tài)環(huán)境建材分會秘書長冀志江透露，制造硅藻泥的原料之一硅藻土在我國的儲量居世界第二和亞洲首位，遠(yuǎn)景儲量超過20億噸。

Numerical Simulation of Thermal Conductivity of Carbon-fiber /Resin Composites

HU Niu，LI Bu-nan，GAO Ben-zheng
（New Materials Institute , Graduate School at Shenzhen,Tsinghua University, Guangdong Shenzhen 508055, China）

Thermal conductivity of carbon-fiber /resin composites was simulated by finite element method. Simulation models of one-dimensional structure and two-dimensional structure were respectively built. Effects of carbon fiber content, interface contact resistance and diameter of carbon fiber on effective thermal conductivity of carbon-fiber/resin composite were studied. The results show that thermal conductivity of carbon-fiber /resin composites increases with increasing of carbon fiber volume fraction; the interface contact resistance in the range of 10-3～ 10-5m2K/W has obvious influence on the effective thermal conductivity, while out of range, the influence can be ignored; the diameter of carbon fiber has a remarkable influence on the thermal conductivity when the interface contact resistance is larger, but the influence is very small when the interface contact resistance is lesser.

Finite element method; Carbon-fiber/resin; Thermal conductivity

TQ 018

A

1671-0460（2014）12-2636-04

2014-06-02

胡妞（1987-），女，湖北孝感人，碩士研究生，研究方向：從事炭纖維/樹脂復(fù)合材料的數(shù)值模擬計(jì)算研究。E-mail：huniu0603@126.com。