基于“離位”增韌技術(shù)Z向注射RTM成型的浸潤(rùn)研究

, , , , ,

(1 山東理工大學(xué) 材料科學(xué)與工程學(xué)院，山東 淄博 255049;2 中航復(fù)合材料有限責(zé)任公司 復(fù)合材料中心，北京 100095;3 濟(jì)南大學(xué) 材料科學(xué)與工程學(xué)院 建筑材料制備與測(cè)試技術(shù)重點(diǎn)實(shí)驗(yàn)室，濟(jì)南 250022;4 山東大學(xué) 材料液固結(jié)構(gòu)演變與加工教育部重點(diǎn)實(shí)驗(yàn)室，濟(jì)南 250061)

基于“離位”增韌技術(shù)Z向注射RTM成型的浸潤(rùn)研究

董抒華1,李偉東2,丁妍羽3,賈玉璽4,劉剛2,魏春城1

(1山東理工大學(xué)材料科學(xué)與工程學(xué)院，山東淄博255049;2中航復(fù)合材料有限責(zé)任公司復(fù)合材料中心，北京100095;3濟(jì)南大學(xué)材料科學(xué)與工程學(xué)院建筑材料制備與測(cè)試技術(shù)重點(diǎn)實(shí)驗(yàn)室，濟(jì)南250022;4山東大學(xué)材料液固結(jié)構(gòu)演變與加工教育部重點(diǎn)實(shí)驗(yàn)室，濟(jì)南250061)

針對(duì)“離位”增韌技術(shù)和Z-RTM成型技術(shù)，引入飽和度參數(shù)修正Darcy定律，建立描述樹(shù)脂在纖維預(yù)制件中非穩(wěn)態(tài)流動(dòng)的偏微分方程，研究恒流注射過(guò)程中體積流量、樹(shù)脂黏度和纖維預(yù)制件滲透率等工藝參數(shù)對(duì)非穩(wěn)態(tài)浸潤(rùn)過(guò)程注入壓力的影響，模擬樹(shù)脂在層間未增韌和增韌纖維預(yù)制件束內(nèi)和束間的流動(dòng)。結(jié)果表明：數(shù)值模擬結(jié)果具有可靠性；隨著注射時(shí)間的增加，纖維預(yù)制件內(nèi)部各點(diǎn)的壓力增加；隨著體積流量、樹(shù)脂黏度的增加，注入壓力線性增加，而隨著纖維滲透率的增加，注入壓力減少，符合Darcy定律；實(shí)現(xiàn)了樹(shù)脂在纖維預(yù)制件細(xì)微觀層次浸潤(rùn)的可視化，這種可視化結(jié)果為預(yù)測(cè)樹(shù)脂在預(yù)制件中的宏觀流動(dòng)提供了重要補(bǔ)充，并為實(shí)際工藝提供了一定指導(dǎo)作用。

“離位”增韌；非穩(wěn)態(tài)浸潤(rùn)；有限元模擬；RTM成型

樹(shù)脂傳遞模塑(Resin Transfer Molding, RTM)因其制件具有整體化、低成本、凈尺寸和較高性能的特點(diǎn)而被廣泛應(yīng)用[1]。近年來(lái)，我國(guó)的航空航天事業(yè)正在快速發(fā)展，對(duì)這一先進(jìn)的低成本復(fù)合材料制備技術(shù)有著迫切的需求[2-4]。RTM工藝要求樹(shù)脂具有極低的初始黏度，然而低黏度樹(shù)脂固化后較脆，這導(dǎo)致復(fù)合材料的抗沖擊損傷性能較差。北京航空材料研究院先進(jìn)復(fù)合材料重點(diǎn)實(shí)驗(yàn)室提出了具有自主知識(shí)產(chǎn)權(quán)的“離位”增韌技術(shù)和Z向流動(dòng)RTM成型技術(shù)，從而成功解決了RTM成型樹(shù)脂低黏度和復(fù)合材料高韌性之間的矛盾[5,6]。對(duì)于“離位”增韌復(fù)合材料，一些學(xué)者已經(jīng)研究了增韌層對(duì)制品力學(xué)性能的影響[7]。然而，“離位”增韌技術(shù)的成功與其成型工藝密切相關(guān)，增韌層的引入不僅改變了樹(shù)脂的充模流動(dòng)速率，而且改變了樹(shù)脂的流動(dòng)模式，因此，“離位”增韌RTM成型中樹(shù)脂流動(dòng)行為較之傳統(tǒng)工藝更為復(fù)雜。

RTM成型中樹(shù)脂對(duì)纖維增強(qiáng)體的充分浸潤(rùn)是影響產(chǎn)品質(zhì)量的一個(gè)重要影響因素。如果成型過(guò)程中樹(shù)脂對(duì)纖維不能完全浸潤(rùn)，就會(huì)使得成型后復(fù)合材料產(chǎn)品中產(chǎn)生干斑、氣孔等缺陷，從而降低了復(fù)合材料的性能[8]。由于RTM工藝大多采用非透明模具，實(shí)驗(yàn)過(guò)程中樹(shù)脂浸潤(rùn)纖維的表觀過(guò)程和內(nèi)部過(guò)程都不可見(jiàn)，難以預(yù)測(cè)和控制，因而，常采用數(shù)值模擬方法預(yù)測(cè)充模時(shí)間、模腔內(nèi)壓力場(chǎng)分布、任意時(shí)刻樹(shù)脂流動(dòng)前鋒位置等，進(jìn)而揭示樹(shù)脂的流動(dòng)規(guī)律，預(yù)測(cè)主要工藝缺陷等[9]。目前，國(guó)內(nèi)外研究人員在RTM工藝充模過(guò)程的模擬方面已開(kāi)展了卓有成效的工作。秦偉等[10]研究了RTM工藝樹(shù)脂流動(dòng)過(guò)程的特點(diǎn)，采用貼體坐標(biāo)/有限差分法模擬了樹(shù)脂滲流過(guò)程；戴福洪等[11]模擬了復(fù)雜形狀三維薄壁構(gòu)件的RTM制造工藝注模過(guò)程，得到任意時(shí)刻的樹(shù)脂流動(dòng)前峰曲線、壓力場(chǎng)分布；Laurenzi等[12]采用有限元方法對(duì)大型復(fù)合材料飛機(jī)橫梁在成型時(shí)充填過(guò)程進(jìn)行了模擬；Tan等[13]提出了在織物中非飽和流動(dòng)的多尺度耦合新方法，采用兩套網(wǎng)格分別計(jì)算束間和束內(nèi)流動(dòng)并進(jìn)行耦合；流體體積(VOF)方法是采用運(yùn)動(dòng)界面技術(shù)跟蹤流動(dòng)前沿的一種方法[14]，而有限元/控制體(FE/CV)方法進(jìn)一步簡(jiǎn)化了算法并被廣泛應(yīng)用[15,16]。在這些研究工作中，樹(shù)脂在模腔中的流動(dòng)被視為流體在多孔介質(zhì)中的流動(dòng)，Darcy定律是RTM成型數(shù)值模擬中的基礎(chǔ)，它在宏觀尺度上描述了流動(dòng)方向上壓力梯度和流體平均速度之間的關(guān)系。在“離位”增韌RTM成型中，樹(shù)脂非穩(wěn)態(tài)浸潤(rùn)過(guò)程中會(huì)同時(shí)發(fā)生在纖維束間和束內(nèi)的細(xì)微觀層次的耦合流動(dòng)，樹(shù)脂的微觀浸潤(rùn)行為對(duì)RTM工藝中樹(shù)脂的宏觀流動(dòng)有重要影響，特別是在纖維含量較高的先進(jìn)復(fù)合材料成型過(guò)程中，樹(shù)脂在纖維預(yù)制件中的細(xì)微觀浸潤(rùn)是決定產(chǎn)品質(zhì)量和性能的關(guān)鍵影響因素[17]。到目前為止，針對(duì)“離位”增韌RTM成型的浸潤(rùn)過(guò)程的數(shù)值模擬鮮有報(bào)道。本工作針對(duì)“離位”增韌技術(shù)和Z-RTM成型技術(shù)，引入了飽和度參數(shù)修正Darcy定律描述樹(shù)脂在纖維預(yù)制件中非穩(wěn)態(tài)流動(dòng)的偏微分方程，這能顯著降低CPU計(jì)算運(yùn)行時(shí)間，從而提高工作效率[18]；研究了無(wú)紡布“離位”增韌RTM成型非穩(wěn)態(tài)浸潤(rùn)過(guò)程注入壓力與時(shí)間關(guān)系，模擬了樹(shù)脂在層間未增韌和“離位”增韌纖維預(yù)制件束內(nèi)和束間的流動(dòng)，實(shí)現(xiàn)了樹(shù)脂在纖維預(yù)制件細(xì)微觀層次浸潤(rùn)的可視化，這對(duì)提高RTM工藝產(chǎn)品質(zhì)量有著非常重要的意義。

1 有限元模型及其數(shù)學(xué)公式

1.1浸潤(rùn)過(guò)程控制方程的建立

樹(shù)脂充填纖維預(yù)制件的過(guò)程可視為流體通過(guò)多孔介質(zhì)過(guò)程。樹(shù)脂在纖維預(yù)制件中的流動(dòng)遵循牛頓流體在多孔介質(zhì)中的流動(dòng)定律Darcy定律[19,20]：

(1)

(2)

式中：φ為纖維預(yù)制件孔隙率，由土壤力學(xué)中飽和度與壓力的關(guān)系可得[18]：

(3)

式中：c(P)是持水量，在流體填充多孔介質(zhì)過(guò)程中：

(4)

式中：α是數(shù)值形狀因子，則可得：

(5)

將式(5)帶入式(3)可得：

(6)

由式(1)，(2)和式(6)可得：

(7)

根據(jù)上述描述樹(shù)脂在纖維預(yù)制件中非穩(wěn)態(tài)流動(dòng)的偏微分方程，采用有限元/控制體法進(jìn)行計(jì)算。Z向注射RTM成型模腔尺寸為200mm×150mm×10mm。模具頂部有9個(gè)出口，底部中心為注射入口。模具的結(jié)構(gòu)示意圖如圖1所示，邊界條件[11]：注射口處：P=P(t)；流動(dòng)前沿P=0；模具邊界：?P/?n=0。

圖1 Z向注射模具結(jié)構(gòu)圖Fig.1 Mould structural diagram of Z-direction injection

1.2流動(dòng)前沿與充填時(shí)間關(guān)系式

當(dāng)樹(shù)脂單向流動(dòng)時(shí)，可得：

(8)

聯(lián)合式(1)和(8)可得：

(9)

假設(shè)流動(dòng)前沿的壓力為0，則可得樹(shù)脂在預(yù)制件內(nèi)的壓力場(chǎng)解析式：

(10)

式中：Pinj是注射口的壓力值；P為樹(shù)脂在纖維預(yù)制件內(nèi)離注射口距離為X處的壓力；Xf是流動(dòng)前沿的位置，結(jié)合式(1)可得模腔內(nèi)樹(shù)脂宏觀浸潤(rùn)速率：

(11)

在恒壓注射條件下，由式(11)可得樹(shù)脂流動(dòng)前沿位置與浸潤(rùn)時(shí)間的關(guān)系：

(12)

由式(12)可得流動(dòng)前沿與時(shí)間的關(guān)系式：

(13)

2 有限元模型的驗(yàn)證

2.1有限元模型與解析解的比較驗(yàn)證

為驗(yàn)證RTM工藝的浸潤(rùn)控制方程及有限元求解的正確性與可靠性，將一維流動(dòng)數(shù)值模擬結(jié)果與解析解進(jìn)行比較。圖2是分別采用數(shù)值模擬和解析解方法在恒壓注射時(shí)流動(dòng)前沿與充填時(shí)間的關(guān)系曲線。從圖中可看出，數(shù)值模擬結(jié)果與解析解結(jié)果吻合一致，采用數(shù)值模擬獲得的填充時(shí)77.5s，而解析解獲得的填充時(shí)間為80s，兩者相對(duì)誤差為3%，證明了數(shù)值模擬的合理性和正確性。

圖2 流動(dòng)前沿與充填時(shí)間的關(guān)系Fig.2 Flow front location vs filling time

2.2有限元模型與實(shí)驗(yàn)結(jié)果的比較驗(yàn)證

未增韌與增韌纖維預(yù)制件的注射壓力隨時(shí)間變化曲線如圖3所示，纖維體積分?jǐn)?shù)為57.5%，未增韌預(yù)制件Z向滲透率為2.9×10-14m2，面內(nèi)滲透率為2×10-13m2；尼龍無(wú)紡布增韌預(yù)制件Z向滲透率為3.6×10-15m2，面內(nèi)滲透率為1.6×10-14m2；測(cè)試液體麥芽糖漿水溶液黏度為0.05Pa·s。采用恒流注射方式，以5mL/min流量沿下模面中心位置注入模腔，流體經(jīng)過(guò)樹(shù)脂分配器快速均勻注入。從圖3中可看出，注入口壓力隨著浸潤(rùn)時(shí)間增加而呈線性增加，在1560s時(shí)刻未增韌與增韌纖維預(yù)制件的模擬計(jì)算的注入壓力與實(shí)驗(yàn)結(jié)果相比，相對(duì)誤差分別為16%和5%。結(jié)果表明在一定范圍內(nèi)該模型能夠精確預(yù)測(cè)樹(shù)脂在纖維預(yù)制件中的非穩(wěn)態(tài)流動(dòng)。對(duì)于剛性高孔隙率的樹(shù)脂分配器，在樹(shù)脂注入時(shí)內(nèi)部的壓力損失很小，樹(shù)脂迅速充滿分配器并同時(shí)以近一維線性流動(dòng)方式進(jìn)行流動(dòng)。隨著流動(dòng)時(shí)間增加，尼龍無(wú)紡布增韌的纖維預(yù)制件注射壓力遠(yuǎn)大于非增韌的纖維預(yù)制件注射壓力，這說(shuō)明在纖維預(yù)制件層間引入增韌層后，使得模具內(nèi)纖維體積分?jǐn)?shù)升高，滲透率降低，因而“離位”增韌的預(yù)制件注射壓力升高且遠(yuǎn)大于非增韌的纖維預(yù)制件的注入壓力。

圖3 預(yù)制件在Z向注射過(guò)程中注入壓力變化曲線Fig.3 Inlet pressure variation of preform during Z-direction injection

3 工藝參數(shù)對(duì)樹(shù)脂浸潤(rùn)過(guò)程的影響

3.1樹(shù)脂在充填時(shí)沿Z向流動(dòng)的壓力分布

圖4是恒流條件下樹(shù)脂在增韌預(yù)制件面內(nèi)中心沿Z向流動(dòng)的壓力分布情況(No.1～5分別是離注射口距離為0，2，4，6，8mm的位置)。從圖中可以看出，模腔內(nèi)各點(diǎn)壓力隨著注射時(shí)間的增加而增加，這主要是因?yàn)殡S著注射時(shí)間的延長(zhǎng)，流動(dòng)前沿增加，由Darcy定律知，為了保證流動(dòng)速度和壓力梯度不變，預(yù)制件內(nèi)部點(diǎn)壓力與前沿壓力差增加，因此，該點(diǎn)壓力隨著流動(dòng)前沿推進(jìn)而增加。圖5是樹(shù)脂完全浸潤(rùn)預(yù)制件后其內(nèi)部的壓力場(chǎng)分布等值面圖。從圖中可看出，從入口到出口壓力逐漸減小，在出口之間模具內(nèi)部接近上表面處由于壓力差變小使得樹(shù)脂流動(dòng)困難，從而易出現(xiàn)干斑等缺陷，因此需要根據(jù)實(shí)際情況關(guān)閉某一個(gè)或幾個(gè)出口閥門(mén)進(jìn)行多個(gè)循環(huán)充填以避免缺陷。

圖4 增韌預(yù)制件在Z向注射過(guò)程中內(nèi)部各點(diǎn)的壓力變化Fig.4 Inner pressure variation of toughened preform during Z-direction injection

3.2體積流量對(duì)注入壓力的影響

圖6是“離位”增韌預(yù)制件在恒流注射時(shí)注入壓力與樹(shù)脂體積流量的關(guān)系。從圖中看出，隨著體積流量增加樹(shù)脂注入壓力線性增加，為了保證樹(shù)脂在預(yù)制件中的恒流速，隨著樹(shù)脂的浸潤(rùn)注入壓力增加，控制體積流量可調(diào)控注入壓力，體積流量過(guò)大則壓力過(guò)大，過(guò)高的壓力會(huì)使流體快速?zèng)_擊到纖維預(yù)制件中，使其局部變形，孔隙率降低，使得滲透率分布不均勻，從而使得產(chǎn)品易出現(xiàn)缺陷，壓力過(guò)高也會(huì)損壞模具，降低設(shè)備使用壽命，而過(guò)低的體積流速會(huì)使得壓力過(guò)低，不能保證纖維預(yù)制件完全被浸潤(rùn)。

圖5 預(yù)制件被完全浸潤(rùn)后其內(nèi)部壓力等值面圖Fig.5 Inner pressure isosurface of preform after fully infiltrated

圖6 預(yù)制件注入壓力與體積流量的關(guān)系Fig.6 Injection pressure of preform vs volume flow rate

3.3樹(shù)脂黏度對(duì)注入壓力的影響

圖7是“離位”增韌預(yù)制件在恒流注射時(shí)注入壓力與樹(shù)脂黏度的關(guān)系。從圖中看出，隨著樹(shù)脂黏度的增加樹(shù)脂注入壓力線性增加。黏度是影響浸漬速度和浸漬效果的重要因素之一，黏度過(guò)高，則需要較大的泵壓力，使得模具厚度必須增加，而模內(nèi)局部纖維有可能被沖走或移位，導(dǎo)致產(chǎn)品出現(xiàn)缺陷，另外，浸漬速度慢，不適宜長(zhǎng)程浸漬；黏度太低，浸漬速度太快，容易夾帶空氣，使制品出現(xiàn)針孔缺陷。因此，要求樹(shù)脂體系具有適宜的黏度特性，有利于樹(shù)脂流動(dòng)浸漬增強(qiáng)材料，使得大尺寸纖維預(yù)成型體的樹(shù)脂能夠長(zhǎng)程滲流充分浸漬，從而實(shí)現(xiàn)大型復(fù)合材料構(gòu)件的成型。

圖7 預(yù)制件注射壓力與黏度的關(guān)系Fig.7 Injection pressure of preform vs viscosity

3.4纖維預(yù)制件滲透率對(duì)注入壓力的影響

圖8 預(yù)制件注射壓力與滲透率的關(guān)系Fig.8 Injection pressure preform vs permeability

圖8是“離位”增韌預(yù)制件在恒流注射時(shí)注入壓力與預(yù)制件滲透率的關(guān)系。從圖中看出，隨著纖維預(yù)制件滲透率的增加樹(shù)脂注入壓力減小。滲透率是影響樹(shù)脂在纖維復(fù)合材料預(yù)制件中流動(dòng)的一個(gè)重要參數(shù)，它描述了預(yù)制件對(duì)樹(shù)脂流動(dòng)的阻礙作用。在恒流注射時(shí)它是影響注入壓力大小的重要因素之一，隨著纖維預(yù)制件滲透率的增加，纖維預(yù)制件孔隙率增大，保持恒速時(shí)注入壓力則隨之降低。

3.5樹(shù)脂在纖維預(yù)制件內(nèi)的流動(dòng)行為

為了研究“離位”增韌Z-RTM成型過(guò)程中樹(shù)脂在纖維預(yù)制件纖維束間和束內(nèi)的浸潤(rùn)流動(dòng)規(guī)律，按[0/45/90/-45]鋪層方式建立幾何模型，增韌層的尼龍無(wú)紡布的纖維是面內(nèi)無(wú)序隨機(jī)分布，因此本工作中尼龍無(wú)紡布被視為面內(nèi)各向同性介質(zhì)。圖9是鋪層纖維預(yù)制件的內(nèi)部浸潤(rùn)圖，其中圖9(a)是在沒(méi)有增韌層時(shí)樹(shù)脂在纖維預(yù)制件內(nèi)的浸潤(rùn)行為圖。從圖中可看出，在一定壓力梯度的驅(qū)動(dòng)下，樹(shù)脂沿預(yù)制件Z向流入纖維束間和束內(nèi)孔隙，樹(shù)脂優(yōu)先通過(guò)纖維束之間較大的流道空隙，同時(shí)浸入纖維束內(nèi)部，因?yàn)槔w維束的滲透率遠(yuǎn)小于束間流道滲透率，因而纖維束內(nèi)部的流動(dòng)前沿滯后，造成束內(nèi)和束間的流動(dòng)前沿參差不齊，在纖維束內(nèi)未浸潤(rùn)完成情況下有可能流動(dòng)前沿已經(jīng)在大的空隙前匯集，因此容易形成空隙缺陷。圖9(b)是在有增韌層時(shí)樹(shù)脂在纖維預(yù)制件內(nèi)的浸潤(rùn)行為圖。從圖中可看出，當(dāng)樹(shù)脂沿預(yù)制件纖維層束間優(yōu)先流過(guò)后，前沿遇到增韌層，增韌層對(duì)樹(shù)脂流動(dòng)產(chǎn)生了一定阻礙，從而使得樹(shù)脂在前沿處橫向匯合，這在一定程度上延緩了樹(shù)脂前沿向前擴(kuò)展，避免了與束內(nèi)流動(dòng)前沿造成大的差距并產(chǎn)生空隙，因此，增韌層對(duì)流動(dòng)前沿起到了二次分配作用。圖9(b)與9(a)相比較，還可看出，樹(shù)脂在“離位”增韌預(yù)制件中流動(dòng)的前沿相對(duì)平緩和平滑，樹(shù)脂在纖維束間和束內(nèi)的流動(dòng)前沿相對(duì)一致，避免了氣體的裹入和浸潤(rùn)不完全，從而減少了孔隙和干斑缺陷的形成，提高了產(chǎn)品的質(zhì)量。

圖10(a)，(b)分別是未增韌的纖維預(yù)制件與增韌預(yù)制件的內(nèi)部形貌照片。從圖中可看出，引入層間增韌層后，缺陷顯著減少。樹(shù)脂在纖維預(yù)制件的纖維束內(nèi)和束間浸潤(rùn)過(guò)程的可視化為預(yù)測(cè)樹(shù)脂在預(yù)制件中的宏觀流動(dòng)提供了重要補(bǔ)充并為實(shí)際工藝提供了一定指導(dǎo)作用。

4 結(jié)論

(1) 建立了描述樹(shù)脂在纖維預(yù)制件中非穩(wěn)態(tài)流動(dòng)的偏微分方程，對(duì)比了解析解和實(shí)驗(yàn)結(jié)果，證明了數(shù)值模擬結(jié)果具有合理性與可靠性。

(2) 隨著注射時(shí)間的增加，纖維預(yù)制件內(nèi)部的各點(diǎn)壓力增加；隨著體積流量、樹(shù)脂黏度的增加注入壓力線性增加，而隨著纖維滲透率的增加，注入壓力減少，符合Darcy定律。

圖9 鋪層纖維預(yù)制件的內(nèi)部浸潤(rùn)圖 (a)未增韌；(b)層間增韌Fig.9 Inner infiltration simulation of fiber preform (a)untoughened；(b)interlamination toughened

圖10 纖維預(yù)制件的內(nèi)部形貌 (a)未增韌；(b)層間增韌Fig.10 Morphologies in fiber preform (a)untoughened；(b)interlamination toughened

(3) 實(shí)現(xiàn)了樹(shù)脂在纖維預(yù)制件細(xì)微觀層次浸潤(rùn)的可視化，這種可視化結(jié)果為預(yù)測(cè)樹(shù)脂在預(yù)制件中的宏觀流動(dòng)提供了重要補(bǔ)充，并為實(shí)際工藝提供了一定指導(dǎo)作用。

[1]BODAGHIM,CRISTVOC,GOMESR,etal.ExperimentalcharacterizationofvoidsinhighfibrevolumefractioncompositesprocessedbyhighinjectionpressureRTM[J].CompositesPartA:AppliedScienceandManufacturing,2016,82(3):88-99.

[2] 益小蘇, 許亞洪, 程群峰, 等. 航空樹(shù)脂基復(fù)合材料的高韌性化研究進(jìn)展[J]. 科技導(dǎo)報(bào),2008,26(6):84-92.

YIXS,XUYH,CHENGQF,etal.Developmentofstudiesonpolymermatrixaircraftcompositematerialshighlytoughened[J].Science&TechnologyReview,2008,26(6):84-92.

[3] 劉偉, 曹臘梅, 王嶺, 等.RTM成型工藝對(duì)Cf/SiBCN陶瓷基復(fù)合材料性能的影響[J]. 材料工程,2015,43(6):1-6.

LIUW,CAOLM,WANGL,etal.EffectofRTMprocessonthepropertiesofCf/SiBCNceramicmatrixcomposites[J].JournalofMaterialsEngineering,2015,43(6):1-6.

[4]NGUYENVH,DELéGLISE-LAGARDéREM,PARKCH.Modelingofresinflowinnaturalfiberreinforcementforliquidcompositemoldingprocesses[J].CompositesScienceandTechnology,2015,113:38-45.

[5] 李偉東, 劉剛, 安學(xué)鋒, 等.Z向流動(dòng)RTM工藝樹(shù)脂的流動(dòng)浸潤(rùn)行為[J]. 復(fù)合材料學(xué)報(bào),2013,30(6):82-89.

LIWD,LIUG,ANXF,etal.InvestigationofresinflowingandinfiltrationbehaviorduringZdirectionflowingRTMprocess[J].ActaMateriaeCompositaeSinica,2013,30(6):82-89.

[6] 劉剛, 張朋, 李偉東, 等. 結(jié)構(gòu)化增韌層增韌RTM復(fù)合材料預(yù)成型體的滲透特性[J]. 復(fù)合材料學(xué)報(bào),2015,32(2):586-593.

LIUG,ZHANGP,LIWD,etal.PermeabilityoftoughenedRTMcompositepreformsbystructuraltougheninglayer[J].ActaMateriaeCompositaeSinica,2015,32(2):586-593.

[7] 張朋, 劉剛, 胡曉蘭, 等. 結(jié)構(gòu)化增韌層增韌RTM復(fù)合材料性能[J]. 復(fù)合材料學(xué)報(bào),2012,29(4):1-9.

ZHANGP,LIUG,HUXL,etal.PropertiesoftoughenedRTMcompositesbystructuraltougheninglayer[J].ActaMateriaeCompositaeSinica,2012,29(4):1-9.

[8] 董抒華, 王成國(guó), 賈玉璽, 等. 纖維復(fù)合材料預(yù)制件滲透率與其結(jié)構(gòu)相關(guān)性的研究進(jìn)展[J]. 材料工程,2013(5):94-100.

DONGSH,WANGCG,JIAYX,etal.Researchprogressonthepermeabilityoffibercompositepreformswithstructuraldependence[J].JournalofMaterialsEngineering,2013(5):94-100.

[9] 齊文, 劉東, 趙俊利, 等.RTM工藝充模過(guò)程模擬研究進(jìn)展[J]. 玻璃鋼/復(fù)合材料,2015(12):105-109.

QIW,LIUD,ZHAOJL,etal.ProgressinnumericalsimulationofmoldfillinginRTM[J].FiberReinforcedPlastics/Composites,2015(12):105-109.

[10] 秦偉, 李海晨, 張志謙, 等.RTM工藝樹(shù)脂流動(dòng)過(guò)程數(shù)值模擬及實(shí)驗(yàn)比較[J]. 復(fù)合材料學(xué)報(bào),2003,20(4):77-80.

QINW,LIHC,ZHANGZQ,etal.ComparisonbetweennumericalsimulationandexperimentalresultofresinflowinRTM[J].ActaMateriaeCompositaeSinica,2003,20(4):77-80.

[11] 戴福洪, 張博明, 杜善義, 等. 復(fù)雜形狀三維薄壁構(gòu)件RTM制造工藝注模過(guò)程模擬[J]. 復(fù)合材料學(xué)報(bào),2004,21(2):87-91.

DAIFH,ZHANGBM,DUSY,etal.Simulationofmould-fillinginRTMprocessfor3Dcomplexshapethinshellparts[J].ActaMateriaeCompositaeSinica,2004,21(2):87-91.

[12]LAURENZIS,GRILLIA,PINNAM,etal.Processsimulationforalargecompositeaeronauticbeambyresintransfermolding[J].CompositesPartB:Engineering,2014,57:47-55.

[13]TANH,PILLAIKM.Multiscalemodelingofunsaturatedflowindual-scalefiberpreformsofliquidcompositemoldingI:Isothermalflows[J].CompositesPartA:AppliedScienceandManufacturing,2012,43(1):1-13.

[14]LUOMAJA,VOLLERVR.Anexplicitschemefortrackingthefillingfrontduringpolymermoldfilling[J].AppliedMathematicalModelling,2000,24(8/9):575-590.

[15]LIMST,LEEWI.Ananalysisofthethree-dimensionalresin-transfermoldfillingprocess[J].CompositesScienceandTechnology,2000,60(7):961-975.

[16]SHOJAEIA.Anumericalstudyoffillingprocessthroughmultilayerpreformsinresininjection/compressionmolding[J].CompositesScienceandTechnology,2006,66(11/12):1546-1557.

[17] 李永靜, 晏石林, 嚴(yán)飛, 等. 注射條件對(duì)LCM工藝非飽和流動(dòng)特性影響[J]. 復(fù)合材料學(xué)報(bào),2016,33(11):2688-2697.

LIYJ,YANSL,YANF,etal.TheinfluencesofinjectionconditionsontheunsaturatedflowofLCMprocess[J].ActaMateriaeCompositaeSinica,2016,33(11):2688-2697.

[18]KLUNKERF,ELSENHANSC,ARANDAS,etal.Modellingtheresininfusionprocess,partI:flowmodellingandnumericalinvestigationforconstantgeometries[J].JournalofPlasticsTechnology,2011,7(5):179-200.

[19]NGOND,TAMMAKK.Microscalepermeabilitypredictionsofporousfibrousmedia[J].InternationalJournalofHeatandMassTransfer,2001,44(16):3135-3145.

[20] 楊波, 王時(shí)龍, 畢鳳陽(yáng). 基于混合網(wǎng)格方法的VARTM工藝充模仿真與實(shí)驗(yàn)驗(yàn)證[J/OL]. 復(fù)合材料學(xué)報(bào),http://www.cnki.net/kcms/detail/11.1801.TB.20161115.1550.016.html.

YANGB,WANGSL,BIFY.Simulationandexperimentalvalidationforthemold-fillingprocessofVARTMbasedonmixedgridapproach[J/OL].MateriaeCompositaeSinica,http://www.cnki.net/kcms/detail/11.1801.TB.20161115.1550.016.html.

[21]LINM,HAHNHT,HUHH.Afiniteelementsimulationofresintransfermoldingbasedonpartialnodalsaturationandimplicittimeintegration[J].CompositesPartA:AppliedScienceandManufacturing,1998,29(5/6):541-550.

(本文責(zé)編：解 宏)

InfiltrationofZ-directionInjectionRTMProcessBasedonEx-situTougheningTechnology

DONGShu-hua1,LIWei-dong2,DINGYan-yu3,JIAYu-xi4,LIUGang2,WEIChun-cheng1

(1SchoolofMaterialsScienceandEngineering,ShandongUniversityofTechnology,Zibo255049,Shandong,China;2CompositeCenter，AVICCompositeCorporationLtd.,Beijing100095,China;3ShandongProvincialKeyLaboratoryofPreparationandMeasurementofBuildingMaterials,SchoolofMaterialsScience&Engineering,UniversityofJinan,Jinan250022,China;4KeyLaboratoryforLiquid-SolidStructuralEvolution&ProcessingofMaterials(MinistryofEducation),ShandongUniversity,Jinan250061,China)

Aimed atex-situtoughening technology andZ-direction RTM process，Darcy’s law was modified by introducing the saturation parameter. The partial differential equation describing the unsteady flow of the resin in the fiber preform was established. The effect of process parameters such as volume flow rate, resin viscosity and fiber preform’s permeability during the constant flow process on the injection pressure was investigated. The resin flow between intra-tow and inter-tow of the preform with untoughened layers and toughened layers was simulated. The results show that the numerical simulation results are reliable. The inner pressure in the fiber performs increases with the increase of injection time. The injection pressure increases linearly with the increase of volume flow rate and resin viscosity, while decreases with the increase of fiber preform’s permeability, which accords with Darcy’s law. The infiltration visualization of resin flow through meso-scale and micro-scale fiber preform is realized, which provides an important supplement for prediction of the macro-flow in fiber preforms and provides guidance for actual process.

ex-situtoughening;unsteady infiltration;finite element simulation;RTM process

10.11868/j.issn.1001-4381.2016.001418

TB332

： A

： 1001-4381(2017)09-0052-07

國(guó)家自然科學(xué)基金(51173100,51373090)；山東省自然科學(xué)基金(ZR2015QZ05,ZR2014EMQ014,ZR2014JL032)

2016-11-28；

：2017-04-05

董抒華(1975-)，女，副教授，博士，主要研究方向?yàn)楦叻肿訌?fù)合材料的制備與仿真,聯(lián)系地址：山東省淄博市山東理工大學(xué)材料科學(xué)與工程學(xué)院(255049)，E-mail:dongshuhua@sdut.edu.cn