二維三軸編織復(fù)合材料預(yù)壓?jiǎn)伟Ｐ徒⒓捌鋸椥砸?guī)律數(shù)值預(yù)測(cè)

張芳芳， 段永川

(1. 燕山大學(xué) 機(jī)械工程學(xué)院， 河北 秦皇島 066004；2. 燕山大學(xué) 先進(jìn)鍛壓成形技術(shù)與科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室， 河北 秦皇島 066004)

二維三軸編織復(fù)合材料由3個(gè)方向的纖維束在平面內(nèi)相互編織而成，這樣在不同角度下可兼顧材料的橫向力學(xué)行為，有效克服傳統(tǒng)層合板力學(xué)性能差的缺點(diǎn)。由于編織復(fù)合材料中纖維束相互交織，因此，編織纖維束的截面形狀在不同位置會(huì)發(fā)生畸變，構(gòu)造實(shí)際的纖維束幾何模型已十分困難，如果再對(duì)纖維束幾何模型離散網(wǎng)格構(gòu)造有限元模型更加困難。

目前多數(shù)學(xué)者都是基于一定的模型簡(jiǎn)化來(lái)構(gòu)造編織復(fù)合材料的幾何和有限元網(wǎng)格模型。Tsai等[1]采用彈簧單元建立了簡(jiǎn)化的單胞網(wǎng)格模型，并對(duì)彈性力學(xué)參數(shù)進(jìn)行了預(yù)測(cè)。Miravete等[2]通過沿編織方向增加4條邊界線，建立了簡(jiǎn)化的“米”字形單胞模型，雖降低了構(gòu)造模型的難度，但也降低了復(fù)合材料彈性、損傷等力學(xué)性能的預(yù)測(cè)精度。Benzley等[3]通過模擬對(duì)比得出，六面體單元具有更多的自由度和抗畸變能力，在力學(xué)性能預(yù)測(cè)時(shí)能得到更高的精度[4-5]。Kim等[6]采用一種體素單元法建立了復(fù)合材料模型，該方法在表征復(fù)合材料兩相材料的交界面時(shí)精度降低，但這種方法隨著分辨率的提高，可消除這個(gè)問題。Gao等[7]假設(shè)紗線橫截面為矩形，建立出二維三軸編織預(yù)制件的動(dòng)態(tài)幾何模型。嚴(yán)雪等[8]基于傳統(tǒng)有限元建模方法建立了考慮纖維束彎曲扭轉(zhuǎn)和空間交錯(cuò)特征的有限元模型，并預(yù)測(cè)了其彈性性能規(guī)律。Kier等[9]假設(shè)纖維束橫截面為橢圓形，編織纖維束路徑為正弦曲線，利用CAD軟件生動(dòng)再現(xiàn)了二維三軸編織復(fù)合材料的幾何模型。張芳芳等[10-11]基于直線的纖維束路徑建立了三維編織復(fù)合材料模型，并對(duì)其進(jìn)行了損傷分析，然而當(dāng)纖維束路徑波動(dòng)較大時(shí)預(yù)測(cè)誤差會(huì)逐漸增大。

目前學(xué)者大都利用簡(jiǎn)化模型對(duì)材料的微觀形態(tài)進(jìn)行表征，但復(fù)合材料纖維束在固化時(shí)會(huì)受到一定的擠壓變形，空間形態(tài)改變會(huì)對(duì)材料的彈性及損傷等性能造成影響?；诖耍疚奶岢龌诹W(xué)原理考慮空間擠壓的模型快速建立方法。利用變形后纖維束空間分布及形態(tài)再形成體素單胞模型，基于該模型預(yù)測(cè)了二維三軸編織復(fù)合材料的彈性規(guī)律。

1 細(xì)觀幾何模型

1.1 編織結(jié)構(gòu)

二維三軸編織復(fù)合材料的編織結(jié)構(gòu)如圖1所示。復(fù)合材料由3個(gè)方向的纖維束組成，其中與y軸平行的為直纖維束，另外2個(gè)方向?yàn)榫幙椑w維束，他們沿著軸向纖維束相互交織，在z軸方向上看，軸向纖維束位于2個(gè)編織纖維束中間，在厚度方向上軸向纖維束位于材料中間。圖中白色框?yàn)橐粋€(gè)周期單元，稱為一個(gè)單胞，其寬度為W，高度為L(zhǎng)，編織角指軸向纖維束與編織向纖維束所夾的銳角(見圖1 中α)，單胞內(nèi)編織向纖維束長(zhǎng)Lb。在圖中所示坐標(biāo)系下，y軸與復(fù)合材料軸向纖維束的軸線平行，x軸與復(fù)合材料軸向纖維束的軸向垂直。z軸與復(fù)合材料的厚度方向平行。

圖1 二維三軸編織復(fù)合材料的結(jié)構(gòu)示意圖

將纖維束中心線定義為空間函數(shù)，選用三次樣條曲線描述纖維束路徑，其中第i段樣條曲線pi(t)的參數(shù)方程一般形式為

pi(t)=Ai+Bit+Cit2+Dit3

(1)

式中：t為該樣條曲線的參數(shù)， 0≤t≤1,i=0,1,…,n-2，Ai、Bi、Ci、Di分別為第i段樣條曲線的冪次項(xiàng)系數(shù)向量。該方程是在單胞內(nèi)定義，在單胞兩端面上要滿足幾何連續(xù)性，在樣條曲線的起始點(diǎn)處應(yīng)滿足的連續(xù)條件為

p″0(0)=p″n-2(1)

(2)

為驗(yàn)證樣條曲線建立的正確性，在空間取了5個(gè)坐標(biāo)點(diǎn)，每個(gè)坐標(biāo)點(diǎn)的坐標(biāo)分別為(0,0,0)，(0.5,0.7,0.4)，(1,0.6,0.0)，(1.5,0.7,-0.4)和(2,1.2,0)。在邊界約束條件下生成的周期樣條曲線如圖2所示?？梢钥闯?，樣條線光滑連續(xù)，可以滿足實(shí)際要求。

圖2 周期性空間樣條曲線

參考文獻(xiàn)[12]中的纖維束橫截面拓?fù)湫螤?，將纖維束的橫截面假設(shè)為透鏡形，如圖3所示。透鏡形纖維束橫截面由半徑分別為r1和r2的2個(gè)圓相交而得，2個(gè)圓各自偏離圖中坐標(biāo)系橫軸的距離是o1和o2。參數(shù)r1、r2、o1和o2可根據(jù)透鏡的寬度w、高度h和偏移距離d計(jì)算得到。

圖3 透鏡形纖維束截面

透鏡形橫截面的參數(shù)方程為：

(3)

(4)

(5)

(6)

(7)

1.2 編織結(jié)構(gòu)參數(shù)

根據(jù)圖1所示幾何關(guān)系，單胞內(nèi)編織向纖維長(zhǎng)Lb為

(8)

式中：L為單胞高度，mm；α為編織角，(°)。

定義軸向纖維束橫截面面積為Aa及圓弧段所對(duì)圓心角為αa，其計(jì)算公式分別為：

(9)

(10)

式中：ra為圓弧段的半徑，mm；w、h分別為透鏡形截面的寬度和高度，mm。同理，可計(jì)算出編織向纖維束橫截面面積Ab和圓心角αb。從圖1可以看出，在一個(gè)矩形單胞模型中含有2根軸向纖維束和4根編織向纖維束，則編織向纖維束和軸向纖維束的體積分別為：

Va=2AaLa

(11)

Vb=4AbLb

(12)

(13)

式中：Va、Vb分別為單胞中軸線纖維束和編織向纖維束的體積，mm3；Aa、Ab分別為軸線纖維束和編織向纖維束的橫截面面積，mm2；La、Lb分別為單胞中軸線纖維束和編織向纖維束的長(zhǎng)度，mm；Vf為單胞中的纖維體積含量，%；L、W、H分別為單胞模型的長(zhǎng)度、寬度和高度，mm；φ為紗線填充因子，本文認(rèn)為編織向纖維束與軸向纖維束的紗線填充因子相同。

2 體素法建立有限元模型

2.1 單元材料主方向計(jì)算

纖維束在空間的走向可以用其中心線在空間的走向進(jìn)行描述，由于纖維束路徑的波動(dòng)導(dǎo)致纖維束各個(gè)部位有著不同的材料主方向。應(yīng)用有限元方法對(duì)實(shí)體模型進(jìn)行離散后，為準(zhǔn)確地描述各單元的材料主方向，需要計(jì)算單元形心處的材料主方向，如圖4所示。令第i段曲線中某單元的形心點(diǎn)為P，過該點(diǎn)做垂直于纖維束中心線的平面，計(jì)算出該平面與纖維束中心線的交點(diǎn)，交點(diǎn)處的參數(shù)坐標(biāo)為tp，該交點(diǎn)沿纖維束走向的切線向量P′i(tp)計(jì)算方程為：

P′i(tp)·(P-Pi(tp))=0

(14)

圖4 纖維方向求解

依據(jù)右手法則，確定單元材料坐標(biāo)系與總體坐標(biāo)系的變換方向和角度。本文采用ANSYS軟件對(duì)有限元模型進(jìn)行建立和計(jì)算，在該軟件內(nèi)采用LOCAL命令為每個(gè)單元建立局部坐標(biāo)系，該坐標(biāo)系采用2個(gè)相對(duì)轉(zhuǎn)動(dòng)角度θxy和θxz定義，如圖5所示。假設(shè)一個(gè)單元的纖維方向?yàn)長(zhǎng)(x,y,z)，L(x,y,z)與圖中坐標(biāo)系的z軸構(gòu)成一個(gè)平面，其法線方向?yàn)镈(x,y,z)。

圖5 材料主方向轉(zhuǎn)角定義

圖5中x、y和z軸與全局笛卡爾坐標(biāo)系的三軸平行，其中θ為y軸與向量D所形成的銳角，β為向量L與z軸所成銳角，2個(gè)相對(duì)轉(zhuǎn)角θxy和θxz的計(jì)算公式為：

(15)

(16)

2.2 擠壓力學(xué)模型的建立

為考慮纖維束間的擠壓扭曲影響，首先建立二維三軸編織預(yù)制件模型如圖6(a)所示。因預(yù)制件未與基質(zhì)黏接固化，無(wú)法構(gòu)造固化后預(yù)制件的周期單胞，對(duì)預(yù)制件進(jìn)行預(yù)壓時(shí)，采用材料的真實(shí)邊界進(jìn)行約束，一個(gè)中心單胞利用其領(lǐng)域內(nèi)8個(gè)單胞進(jìn)行自然邊界的約束。在預(yù)制件有限元模型中纖維上下表面分別加入接觸單元，在上下2層分別利用平面進(jìn)行一定的預(yù)壓縮，將預(yù)制件模型厚度壓縮至目標(biāo)單胞模型厚度，從而得到纖維束呈擠壓接觸狀態(tài)下的幾何形態(tài)。開發(fā)周期單胞提取程序，該程序基于FORTRAN語(yǔ)言編寫，具有與ANSYS軟件兼容的進(jìn)出接口，首先將變形后的單元節(jié)點(diǎn)信息導(dǎo)出，給出切取單胞的邊界盒坐標(biāo)，在邊界面上程序可自動(dòng)細(xì)分單元，完成周期單胞單元提取。為便于后期處理，提取前對(duì)不同方位的纖維束進(jìn)行集合劃分，單元切割程序可保持單元的集合和單元坐標(biāo)系的繼承能力，多個(gè)單元可能由一個(gè)父單元切割再分而來(lái)，這幾個(gè)單元就要保持與父單元同樣的集合和單元坐標(biāo)系方位。切割出的增強(qiáng)相單胞網(wǎng)格模型，如圖6(b)所示。切割出的增強(qiáng)相單胞模型單元坐標(biāo)系方位顯示如圖6(c)所示。

圖6 二維三軸編織預(yù)制件網(wǎng)格提取

2.3 體素網(wǎng)格模型的建立

應(yīng)用上述切割算法已經(jīng)獲取了變形后纖維束單胞模型，本文利用體素模型法建立變形后的單胞模型，在該模型區(qū)域內(nèi)填充正方體素單元。根據(jù)已提取的預(yù)制件單胞單元來(lái)對(duì)所填充的正方體素單元進(jìn)行集合創(chuàng)建與單元坐標(biāo)系的繼承。利用預(yù)制件單胞單元在每根纖維束表面生成表面單元，然后判斷體素單元與表面單元的相對(duì)位置，當(dāng)體素單元的形心坐標(biāo)在表面單元的包絡(luò)區(qū)域內(nèi)時(shí)，將該體素單元納入該集合內(nèi)，其他纖維束計(jì)算方法相似，最后生成的體素模型如圖7所示。

圖7 采用體素法建立的復(fù)合材料單胞網(wǎng)格模型

該模型采用全正方六面體網(wǎng)格建立，邊界非常容易滿足邊界節(jié)點(diǎn)一一對(duì)應(yīng)的關(guān)系。

3 周期性邊界條件

復(fù)合材料內(nèi)部細(xì)觀結(jié)構(gòu)周期重復(fù)，在外載荷作用下其應(yīng)力應(yīng)變場(chǎng)同樣滿足周期重復(fù)性，為降低計(jì)算量，選取一個(gè)代表單元進(jìn)行力學(xué)行為預(yù)測(cè)，在單胞2個(gè)對(duì)應(yīng)的邊界面建立周期對(duì)應(yīng)的耦合約束方程[13]，約束表達(dá)式如下:

(i,j,k=x,y,z)

(17)

圖8 編織角和纖維體積含量對(duì)二維三軸編織復(fù)合材料橫向和軸向模量的影響

4 結(jié)果分析

以文獻(xiàn)[12]中所提供的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證，模型中所采用的材料組分力學(xué)性能參數(shù)如表1所示。為驗(yàn)證體素模型的收斂性，不同編織角和纖維體積含量下分別建立了3組不同網(wǎng)格密度的體素模型，進(jìn)行軸向模量的預(yù)測(cè)，并與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，結(jié)果如表2所示。隨著體素模型網(wǎng)格密度的增加，預(yù)測(cè)出的軸向彈性模量變化不大，認(rèn)為得到收斂解。二維三軸編織復(fù)合材料的力學(xué)行為主要受編織角和纖維體積含量影響，在不同編織角和不同纖維體積含量下材料宏觀力學(xué)行為的變化規(guī)律如圖8所示。

表1 組分材料性能參數(shù)

從圖8(a)可見，隨著纖維體積含量的升高，二維三軸編織復(fù)合材料橫向，即坐標(biāo)系的x軸方向的彈性模量Ex有所升高，在材料內(nèi)部主要靠纖維束承載，當(dāng)預(yù)壓縮提高纖維含量時(shí)材料的彈性性能有所升高。

表2 數(shù)值預(yù)測(cè)結(jié)果與實(shí)驗(yàn)對(duì)比

從圖8(b)可以看出，隨著纖維體積含量的升高，材料軸向即坐標(biāo)系的y軸方向的彈性模量Ey有所升高，其升高原因與橫向模量相似。隨著編織角的增大，材料的橫向彈性模量Ex不斷升高，當(dāng)編織角升高時(shí)，軸向纖維束以外的兩軸方向不斷變化，編織角越大，纖維束軸向的性能向單胞橫向貢獻(xiàn)越多。而材料的軸向彈性模量Ey與橫向彈性模量Ex變化規(guī)律正好相反。

從圖8(c)可發(fā)現(xiàn)，當(dāng)提高纖維體積含量時(shí)，材料的縱向彈性性能Ez有所升高，但對(duì)材料模量的提高程度較橫向和軸向小，隨著編織角的增加，3個(gè)軸的纖維束交織重疊區(qū)域減小，這樣造成了在厚度方向上彈性模量逐漸降低，但是降低的幅度很小。

從圖8(d)可以看出，當(dāng)預(yù)壓縮提高纖維含量時(shí),材料的剪切模量Gxy逐漸增大，隨著編織角的增加，剪切模量Gxy逐漸增大而后逐漸減小，當(dāng)編織角在40°左右時(shí)，剪切模量Gxy達(dá)到最大值。從圖8(e)可以看出，剪切模量Gxz隨著纖維體積含量和編織角的增加均呈現(xiàn)出逐漸增大的趨勢(shì)。從圖8(f)可以看出，剪切模量Gyz對(duì)纖維體積含量和編織角變化的敏感程度均較低。

從圖8(g)可以看出，隨著編織角的增加，泊松比μxy呈線性增加，然而改變纖維體積含量對(duì)泊松比μxy影響甚微。從圖8(h)發(fā)現(xiàn)：當(dāng)提高纖維體積含量時(shí),泊松比μxz逐漸減小，編織角較小時(shí)，編織角變化對(duì)泊松比影響顯著；隨著編織角的增加，泊松比近似呈線性降低，當(dāng)編織角大于40°后，編織角變化對(duì)泊松比的影響程度逐漸降低。從圖8(i)可以看出：泊松比μyz隨著編織角的增加呈現(xiàn)先減小后增大的變化趨勢(shì)，并且在編織角較小時(shí)，纖維體積含量的變化對(duì)泊松比影響很大；當(dāng)編織角小于40°，預(yù)壓縮提高纖維含量時(shí),泊松比μyz有較大幅度的增加；當(dāng)編織角大于40°時(shí)，纖維體積含量變化對(duì)泊松比的影響程度逐漸降低；當(dāng)編織角達(dá)到50°時(shí)，纖維體積含量變化對(duì)泊松比的影響程度幾乎為零。

5 結(jié) 論

1) 采用空間周期三次樣條曲線表述二維三軸編織復(fù)合材料中纖維束路徑周期波動(dòng)現(xiàn)象，并求解出任意一點(diǎn)處的單元材料主方向。

2) 考慮擠壓變形建立了預(yù)制件有限元模型，預(yù)壓縮后利用單元切割程序生成了預(yù)制件單胞模型?；陬A(yù)制件單胞網(wǎng)格，利用體素法生成了考慮纖維束間相互擠壓影響的二維三軸編織復(fù)合材料的參數(shù)化單胞模型。

3) 基于參數(shù)化的體素單胞模型，預(yù)測(cè)了二維三軸編織復(fù)合材料的彈性性能規(guī)律。該方法為二維三軸編織復(fù)合材料工藝參數(shù)優(yōu)化編織復(fù)合材料設(shè)計(jì)的軟件化程序化奠定了較好的基礎(chǔ)。