曾翔龍， 王奇志， 蘇 飛

(北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院， 北京100191)

?

含缺陷C/SiC平紋機(jī)織復(fù)合材料拉伸力學(xué)行為數(shù)值模擬

曾翔龍， 王奇志， 蘇 飛

(北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院， 北京100191)

采用APDL語言實現(xiàn)ANSYS的二次開發(fā)，建立含預(yù)制缺陷的纖維束截面卵圓形多尺度單胞模型。首先計算纖維束單胞的初始模量，強(qiáng)度以及最大應(yīng)變；隨后利用掃描電鏡圖中的缺陷建立單胞模型，并引入周期性邊界條件，預(yù)測材料的初始各向材料常數(shù)。同時利用Linde提出的逐漸損傷準(zhǔn)則，進(jìn)行單軸拉伸力學(xué)行為的數(shù)值模擬，并闡述該平紋機(jī)織復(fù)合材料單胞模型在經(jīng)向拉伸載荷作用下其纖維束的損傷及演化過程。該模型計算得到的最大拉應(yīng)度為0.65%，強(qiáng)度為256.46 MPa。結(jié)果表明，該模型給出的數(shù)值模擬結(jié)果與實驗數(shù)據(jù)吻合較好，證明了模型的有效性，為該類材料的優(yōu)化設(shè)計及其力學(xué)性能分析提供了一種有效方法。

平紋機(jī)織C/SiC復(fù)合材料;漸進(jìn)損傷模型;失效準(zhǔn)則;強(qiáng)度預(yù)測

C/SiC復(fù)合材料是一種兼有金屬材料、陶瓷材料和碳材料性能優(yōu)點(diǎn)的熱結(jié)構(gòu)/功能一體化新型材料，具有耐高溫、低密度、高比強(qiáng)、高比模、抗氧化、抗燒蝕，對裂紋不敏感，不發(fā)生災(zāi)難性損毀等特點(diǎn)，自1998年實現(xiàn)國產(chǎn)化以來，在航空航天領(lǐng)域有非常廣泛的應(yīng)用。

國內(nèi)外對C/SiC復(fù)合材料的研究主要集中在實驗上[1-4]。在機(jī)織樹脂基復(fù)合材料數(shù)值模擬上，王新峰[5]、盧子興等[6]基于平紋機(jī)織復(fù)合材料的細(xì)觀結(jié)構(gòu)單胞模型，考慮其制備過程中產(chǎn)生的孔隙缺陷為隨機(jī)分布的特征，基于漸進(jìn)損傷演化方程建立了平紋機(jī)織復(fù)合材料的單軸拉伸模型，利用有限元數(shù)值方法模擬了其拉伸應(yīng)力-應(yīng)變行為。左中鵝等[7]利用ANSYS軟件建立了平紋織物復(fù)合材料代表體積單元(RVE)模型，并根據(jù)復(fù)合材料內(nèi)纖維與基體的失效準(zhǔn)則預(yù)測了模型的拉伸強(qiáng)度。然而，這些數(shù)值方法大都基于樹脂基復(fù)合材料，未能體現(xiàn)陶瓷基體的破壞特點(diǎn)，且大部分基于基體填滿的模型計算，并未考慮孔洞的影響[8-9]。文獻(xiàn)[6]描述基體空隙時采用的是概率分布的方法，一方面，假設(shè)的這些孔隙分布規(guī)律缺乏必要的統(tǒng)計數(shù)據(jù)支持；另一方面，SiC基體中的孔洞與樹脂基體中呈現(xiàn)明顯的不同，由于制作工藝的原因，中間部分的孔隙明顯更大。此外，多數(shù)研究者常以Hashin失效準(zhǔn)則作為基體或纖維的損傷判據(jù)[5, 10-11]，并采用“剛度折減”方法來研究材料的損傷失效過程。由于剛度折減系數(shù)的大小直接影響材料的損傷演化過程，且各個研究者對剛度折減系數(shù)的選取有所不同。因此，采用傳統(tǒng)的“突變”剛度折減方法研究材料的逐漸損傷失效并不十分合理，難以揭示實際材料的逐漸損傷失效的演化過程。

本工作基于考慮纖維束相互擠壓的卵圓形纖維束截面單胞模型，引入周期性位移邊界條件，采用細(xì)觀多尺度有限元方法，基于實際材料的掃描電鏡圖，建立含缺陷C/SiC平紋機(jī)織復(fù)合材料力學(xué)模型，模擬其拉伸力學(xué)行為，并對其強(qiáng)度進(jìn)行預(yù)測。該模型在兩個尺度上對其單胞的力學(xué)性能進(jìn)行考慮。在第一尺度上建立了由SiC基體包裹的T300單絲碳纖維束模型以模擬纖維束的性能退化。在第二尺度上依據(jù)第一尺度均勻化方法得到的力學(xué)性能，利用掃描電鏡圖建立含缺陷平紋編織C/SiC復(fù)合材料單胞模型。結(jié)合周期性邊界條件，利用Linde等[12]提出的逐步失效準(zhǔn)則定義材料的典型損傷類型，通過ANSYS有限元軟件的二次開發(fā)程序APDL，對含缺陷平紋編織復(fù)合材料的單軸載荷下非線性拉伸性能進(jìn)行數(shù)值模擬。該模型基于纖維束和純基體相應(yīng)損傷類型所造成的材料性能退化，模擬單調(diào)拉伸狀態(tài)下?lián)p傷產(chǎn)生、發(fā)展的過程。通過數(shù)值結(jié)果與實驗結(jié)果的對比，驗證模型的適用性。

1 單胞彈性常數(shù)的推斷

(1)

其單胞的平均力學(xué)性能即為整體結(jié)構(gòu)的性能，則有：

(2)

其中[S]為柔度矩陣，若僅對單元體進(jìn)行模量分析，則施加六組線性無關(guān)的位移邊界條件，式(2)可擴(kuò)展為：

(3)

因此，可以求得單胞柔度矩陣[S]：

(4)

根據(jù)單元體柔度矩陣可以得到材料的全部彈性常數(shù)。

2 纖維束單胞模型的建立

2.1纖維束單胞幾何模型

針對陶瓷基復(fù)合材料，這里的纖維束是指經(jīng)過CVI高溫環(huán)境沉積過熱解碳界面和內(nèi)部滲入碳化硅基體后的單向纖維增強(qiáng)復(fù)合材料，是位于單胞下一個層次的力學(xué)分析模型。在纖維束理論分析中，假設(shè)纖維束為纖維在基體中均勻規(guī)則分布，因此可按照圖1中矩形區(qū)域作為纖維束單胞，模型如圖2所示，其中a,b,r的比例由實際制備時C纖維的體積分?jǐn)?shù)決定，本研究針對纖維束中C纖維體積分?jǐn)?shù)為70%的試件進(jìn)行計算。可以看出，該單胞滿足周期性要求，且為規(guī)則的長方體結(jié)構(gòu)，便于周期性邊界條件的施加。

圖1 纖維束代表性體積單元Fig.1 Representative volume unit of fiber bundle

圖2 紗線浸漬體單胞模型Fig.2 Unit cell model of impregnated yarns

2.2纖維束單胞強(qiáng)度準(zhǔn)則

在常規(guī)使用環(huán)境下，一般認(rèn)為基體是均質(zhì)各向同性線彈性材料，各種適合于各向同性材料的強(qiáng)度準(zhǔn)則均可用于判斷基體的損傷。對于SiC陶瓷，拉伸和壓縮破壞的抵抗能力存在很大差別，抗壓強(qiáng)度遠(yuǎn)遠(yuǎn)地大于抗拉強(qiáng)度，因此宜采用莫爾強(qiáng)度理論作為其強(qiáng)度準(zhǔn)則，另一方向，C纖維為強(qiáng)度較高的各向異性材料，在此采用最大拉應(yīng)力準(zhǔn)則判定其破壞。

本研究中通過ANSYS用戶二次開發(fā)程序APDL，在材料本構(gòu)關(guān)系中嵌入損傷準(zhǔn)則和材料性能退化方案。通過不同損傷準(zhǔn)則對組分材料單元進(jìn)行判斷，由此在損傷過程中退化相應(yīng)的材料性能；利用Newton-Raphson方法計算各增量步單胞中單元應(yīng)力應(yīng)變，根據(jù)損傷準(zhǔn)則更新單元本構(gòu)矩陣，非線性求解達(dá)到收斂，再繼續(xù)增加載荷步，如此循環(huán)完成分析。

3 復(fù)合材料單胞模型的建立

3.1單胞幾何模型

在完成纖維束尺度的幾何建模后，由圖3中C/SiC平紋機(jī)織復(fù)合材料板的橫截面掃描電鏡圖可以看出，由于C/SiC平紋機(jī)織復(fù)合材料的CVI制備工藝，板材的中部基體分布較少，且分布形式與通常采用的以Weibull分布布置在基體中的微小缺陷有明顯不同，孔洞顯然更大。本研究為考慮板材中部缺陷較大的極端不利情況，參照文獻(xiàn)[14]給出的二維t模型，按照圖4選取代表性體積單元。這一模型考慮了編織陶瓷基復(fù)合材料纖維束之間的間隙；考慮了陶瓷基復(fù)合材料中的大孔洞，體現(xiàn)了和樹脂基復(fù)合材料建模的極大不同。

圖3 2D C/SiC復(fù)合材料內(nèi)部缺陷Fig.3 Internal defects of 2D C/SiC composites

圖4 2D C/SiC單胞代表性體積單元Fig.4 RVE of 2 D C/SiC unit cell

3.2組分材料失效準(zhǔn)則與損傷演化

C/SiC平紋機(jī)織復(fù)合材料在單向拉伸載荷作用下，材料中損傷逐步產(chǎn)生、發(fā)展直至破壞。在各種碳纖維增強(qiáng)復(fù)合材料力學(xué)行為模擬模型中，大多采用的是達(dá)到某一各向異性材料破壞準(zhǔn)則后殺死單元的方法[5, 15]。顯然，材料的破壞是一個逐步演進(jìn)的過程，突變的方式難以完全揭示實際材料的連續(xù)損傷演化的失效過程；因此，采用Linde等[12]提出的以失效應(yīng)變作為判斷依據(jù)的失效準(zhǔn)則。該準(zhǔn)則假設(shè)材料是連續(xù)損傷失效的，其損傷演化是通過纖維、基體中各自的應(yīng)變能來控制。該準(zhǔn)則的具體形式如下。

(1)纖維的損傷失效[12]

(5)

(6)

式中：Lc為材料單元的特征長度；Gf為纖維的斷裂能密度。

(2)基體的損傷失效

(7)

(8)

4 數(shù)值結(jié)果與討論

4.1紗線浸漬體有限元分析

利用2.1中的幾何模型，采用通用有限元軟件ANSYS對C/SiC纖維束單胞的剛度、強(qiáng)度進(jìn)行分析?；w單元和纖維單元損傷判斷和剛度折減通過自行編寫的ANSYS參數(shù)化編程語言APDL實現(xiàn)。外載荷采用應(yīng)變逐級加載形式，按第1節(jié)中公式計算出各載荷步單胞的平均應(yīng)力，并繪出應(yīng)力應(yīng)變的關(guān)系曲線。所求得的纖維浸漬體單胞各方向材料常數(shù)和各方向強(qiáng)度由表1和表2[16]給出。

表1 纖維束各方向初始彈性模量

表2 纖維束各方向強(qiáng)度Table 2 Strength of fiber bundle in different directions

在單向復(fù)合材料中，縱向拉、壓性能主要是由纖維決定，纖維破壞將導(dǎo)致整體材料的最終失效。圖5給出了纖維束單胞在70%纖維體積含量時X/Z方向的拉伸應(yīng)力-應(yīng)變曲線，其他方向的材料常數(shù)和模量同理可得。其中，Y方向指沿單絲纖維軸向方向，X,Z方向則位于與Y方向垂直的平面中。

圖5 纖維束單胞應(yīng)力-應(yīng)變曲線Fig.5 Stress-strain curve of fiber bundle unit cell

由圖5明顯看出，當(dāng)X/Z方向應(yīng)變達(dá)到0.55%時，應(yīng)力-應(yīng)變曲線存在明顯偏折，這是由于SiC基體剛度相對較大而強(qiáng)度較低，因此首先破壞，并使得整體平均應(yīng)力下降，這一點(diǎn)和樹脂基復(fù)合材料的纖維束單胞破壞主要由C纖維體現(xiàn)[5]有極大不同。

4.2復(fù)合材料單胞有限元模型

在完成纖維束尺度的力學(xué)性能分析后，按照3.1節(jié)圖4中的方法取復(fù)合材料中典型機(jī)織結(jié)構(gòu)為一個單胞，如圖6所示為機(jī)織結(jié)構(gòu)剖面圖[17]，假設(shè)機(jī)織結(jié)構(gòu)中，纖維束截面形狀為如圖示經(jīng)紗所示的壓扁形狀，由兩段相交圓弧組成，圓弧半徑為R1。纖維束由兩段同心圓弧段組成，兩段圓弧半徑分別為R1,R1+h/2。模型中相關(guān)參數(shù)可由電鏡照片中經(jīng)紗緯紗圖形測量得到，并進(jìn)一步得到組成經(jīng)紗與緯紗形狀的各圓弧方程。利用如圖7所示的C/SiC復(fù)合材料局部電鏡掃描圖，依照圖6中的模型截面[17]，可量取其中的具體尺寸。

圖6 二維C/SiC機(jī)織復(fù)合材料單胞模型截面圖Fig.6 Unit-cell sectional drawing of 2D C/SiC composite

圖7 二維C/SiC機(jī)織復(fù)合材料電鏡掃描圖Fig.7 SEM photograph of 2D C/SiC composite

依據(jù)3.1中的幾何模型和量得的尺寸建立如圖8所示的二維C/SiC復(fù)合材料單胞有限元模型，其中的纖維束單元為在第一步計算中得到其初始模量及漸進(jìn)損傷演化規(guī)律的平均化模型。由于纖維束的分布特點(diǎn)，四根交疊纖維束由八個有不同方向材料參數(shù)的部分組成。為了順利施加周期性邊界條件[14]，在有限元網(wǎng)格劃分時，必須保證主從面上相應(yīng)網(wǎng)格節(jié)點(diǎn)一一對應(yīng)，采用網(wǎng)格掃掠的方法來實現(xiàn)，模型主要采用二十節(jié)點(diǎn)六面體solid187單元離散，纖維束與纖維束分界面及纖維束與基體分界面均采用共節(jié)點(diǎn)方式連接，以保證界面處變形的連續(xù)性。此外，由于纖維束單元的形成方式，有限元模型中纖維束單元坐標(biāo)系方向如圖9所示沿軸線旋轉(zhuǎn)。

圖8 二維C/SiC復(fù)合材料有限元模型 (a)基體;(b)紗線Fig.8 Finite element model of C/SiC composites (a)matrix;(b)yarn

圖9 二維C/SiC復(fù)合材料纖維束單元坐標(biāo)方向Fig.9 Orientation of element coordinate system for 2D C/SiC composites

4.3復(fù)合材料單胞細(xì)觀應(yīng)力分布

利用4.2節(jié)中建立的有限元模型，同時代入第4.1節(jié)中求得的纖維束材料常數(shù)，可得到圖10給出的單胞內(nèi)各組分材料在經(jīng)向拉伸載荷作用下的第一主應(yīng)力應(yīng)力分布云圖。由圖10(a)可見，根據(jù)掃描電鏡圖建立的含缺陷單胞中間位置有較大的集中應(yīng)力，極易成為外部基體起裂位置。由圖10(b)可見，經(jīng)向纖維束承受大部分的拉伸載荷，與基體交接的部位由于承受來自彈性模量較大的基體的作用力，因此產(chǎn)生較大的應(yīng)力集中，這是與樹脂基復(fù)合材料最大的區(qū)別。待基體破壞較大后這一部分的損傷會趨于緩解。利用1節(jié)中平均方法求得的初始加載方向拉伸模量為118.827 GPa，比文獻(xiàn)[18]中142.06±13.69 GPa略小，這是由于本研究依照空隙較大的危險狀態(tài)建模，實際情況中更多單胞有較多的基體存在。求得的其他材料常數(shù)如表3所示，其中σu為最大拉伸強(qiáng)度，εf為極限應(yīng)變。

4.4應(yīng)力-應(yīng)變曲線的模擬

利用表1，表2給出的C/SiC各組分材料的性能參數(shù)，其中根據(jù)Total Materia數(shù)據(jù)庫查得的數(shù)據(jù)纖維和基體的斷裂能密度分別取為Gm=1 N/mm2,Gf=17.8 N/mm2。圖11同時給出了本研究中預(yù)測拉伸應(yīng)力-應(yīng)變曲線的結(jié)果與文獻(xiàn)[19]中的實驗曲線及有限元結(jié)果的比較?？梢?，復(fù)合材料的本構(gòu)曲線呈現(xiàn)偽塑性，在加載的初始階段，SiC基體作為主要的承載對象，表現(xiàn)為較高的初始模量，隨著應(yīng)力的增加，纖維束外基體由于失效破壞，導(dǎo)致應(yīng)力-應(yīng)變曲線上出現(xiàn)轉(zhuǎn)折點(diǎn)，隨后纖維束內(nèi)纖維因具有較高的抗拉強(qiáng)度而繼續(xù)承載，應(yīng)力-應(yīng)變又近似呈線性關(guān)系。當(dāng)應(yīng)力水平接近纖維束的抗拉強(qiáng)度時，材料的模量下降較為迅速，并最終達(dá)到失效破壞，表現(xiàn)出明顯的脆性破壞特征。應(yīng)力-應(yīng)變曲線中預(yù)估的C/SiC復(fù)合材料抗拉強(qiáng)度為256.46 MPa，與文獻(xiàn)[19]中265.28±15.95 MPa一致。

圖10 單胞經(jīng)向拉伸應(yīng)力分布 (a)基體;(b)紗線Fig.10 Stress field distribution along wrap direction (a)matrix;(b)yarn

表3 實驗數(shù)據(jù)與計算數(shù)據(jù)對比Table 3 Experimental data compared to calculated data

圖11 二維機(jī)織復(fù)合材料應(yīng)力-應(yīng)變曲線實驗和模擬結(jié)果對比Fig.11 2D woven composites stress-strain experiment curve compared to simulation

4.5損傷演化分析

二維機(jī)織C/SiC復(fù)合材料的力學(xué)行為主要由作為增強(qiáng)相的纖維束決定。在ε=0.035%之前，單胞未發(fā)生損傷，應(yīng)力-應(yīng)變曲線呈現(xiàn)線性；在ε=0.035%之后，基體單元開始發(fā)生損傷，隨著基體大片損傷的開始，應(yīng)力-應(yīng)變曲線開始呈現(xiàn)“偽塑性”；在ε=0.3%時，繼續(xù)加載，纖維束各種損傷相對穩(wěn)定擴(kuò)展，并保持一定的承載能力，材料具有一定的延展性；在ε=0.6%時，材料達(dá)到應(yīng)變極限，隨后材料屈服，承受載荷的能力消失。如圖12(a),(b)分別為C/SiC平紋機(jī)織復(fù)合材料紗線和基體的損傷擴(kuò)展圖。對比兩圖的損傷擴(kuò)展趨勢可以看出，基體由于存在模量大、強(qiáng)度低的特點(diǎn)，因此在較低的應(yīng)變水平下即發(fā)生破壞，而纖維束中由于存在著高強(qiáng)碳纖維使得其強(qiáng)度大為提升，可以看出在應(yīng)變率為ε=0.06%時，纖維束和基體的損傷帶從中部缺陷處產(chǎn)生，隨著應(yīng)變的增加，損傷帶逐漸向周圍擴(kuò)展，且由于建模特點(diǎn)，其損傷場呈現(xiàn)中心對稱特點(diǎn)，當(dāng)應(yīng)變達(dá)到屈服時應(yīng)變時，基體已呈現(xiàn)大范圍破壞趨勢，在加載的末期ε=1.2%時，纖維束中心位置發(fā)生斷裂，材料完全失效。

圖12 C/SiC平紋機(jī)織復(fù)合材料損傷擴(kuò)展圖 (a-1)紗線,ε=0.06% ;(a-2)紗線,ε=0.6% ; (a-3)紗線,ε=1.2%(b-1)基體,ε=0.06%;(b-2)基體,ε=0.2%;(b-3)基體,ε=0.6%Fig.12 Evolution of plain weave damage (a-1)yarn,ε=0.06%;(a-2)yarn,ε=0.6%; (a-3)yarn,ε=1.2%(b-1)matrix,ε=0.06%; (b-2) matrix, ε=0.2%;(b-3) matrix,ε=0.6%

5 結(jié)論

(1) 采取跨尺度分析方法，在多個層次上建立了含缺陷平紋機(jī)織復(fù)合材料的細(xì)觀結(jié)構(gòu)單胞模型，并應(yīng)用APDL語言實現(xiàn)了ANSYS的二次開發(fā)，模擬了其拉伸應(yīng)力-應(yīng)變行為。預(yù)估的最大拉應(yīng)變?yōu)?.65%，強(qiáng)度為256.46 MPa，與文獻(xiàn)中的實驗曲線吻合較好，證明了模型的有效性。

(2)基于該模型可以看出，預(yù)制缺陷后，由于應(yīng)力集中，纖維束強(qiáng)度和極限應(yīng)變與理想填滿狀況相比較小。并從細(xì)觀角度演示了該平紋機(jī)織復(fù)合材料單胞模型在經(jīng)向拉伸載荷作用下其纖維束的損傷模式及其演化過程：對于纖維束而言，破壞從中間應(yīng)力較為集中處開始并隨后向兩側(cè)沿較窄范圍擴(kuò)展；對于基體而言，首先在與受載方向垂直方向產(chǎn)生破壞，并隨即擴(kuò)散到四周基體，發(fā)生大片破壞。

[1] 李龍彪. 纖維增強(qiáng)陶瓷基復(fù)合材料疲勞遲滯回線模型研究[J]. 力學(xué)學(xué)報, 2014，46(5):710-729.

(LI L B. Investigation on fatigue hysteresis loops models of fibre-reinforced ceramic-matrix composites[J]. Chinese Journal of Theoretical and Applied Mechanics, 2014, 46(5): 710-729.)

[2] 葉龍海, 鄧瓊, 湯忠斌. 2D-C/SiC復(fù)合材料低溫動態(tài)Z向壓縮性能[J]. 復(fù)合材料學(xué)報, 2014,31(5):1337-1342.

(YE L H,DENG Q,TANG Z B.Dynamic compressive properties of 2D-C/SiC composites Z-direction in low temperature[J].Acta Materiae Compositae Sinica,2014,31(5):1337-1342.

[3] 張青, 成來飛, 王芙愿, 等. 氧化對2D C/SiC復(fù)合材料阻尼行為的影響[J]. 新型炭材料, 2014，29(6):522-528.

(ZHANG Q, CHENG L F, WANG F Y,etal. Effect of oxidation on the damping behavior of 2D C/ SiC composites[J]. New Carbon Materials, 2014, 29(6): 522-528.)

[4] 馮炎建, 馮祖德, 劉永勝, 等. 2D C/SiC復(fù)合材料高溫蠕變損傷的顯微CT分析[J]. 金屬熱處理, 2011(增刊1):482-485.

(FENG Y J,FENG Z D,LIU Y S,etal.Micro-CT analysis of high temperature creep damage of 2D C/SiC composites[J].Heat Treatment of Metals, 2011(Suppl 1):482-485.)

[5] 王新峰. 機(jī)織復(fù)合材料多尺度漸進(jìn)損傷研究[D].南京：南京航空航天大學(xué), 2007.

(WANG X F. Multi-scale analyses of damage evolution plain woven composites[D]. Nanjing:Nanjing University of Aeronautics and Astronautics,2007.)

[6] 盧子興, 徐強(qiáng), 王伯平, 等. 含缺陷平紋機(jī)織復(fù)合材料拉伸力學(xué)行為數(shù)值模擬[J]. 復(fù)合材料學(xué)報, 2011,28(6):200-207.

(LU Z X,XU Q,WANG B P,etal.Numerical simulation of plain weave composites with defects under unidirectional tension[J].Acta Materiae Compositae Sinica,2011,28(6):200-207.)

[7] 左中鵝, 王瑞, 徐磊. 基于有限單元法的平紋織物復(fù)合材料強(qiáng)度預(yù)測:1.RVE的有限元模型[J]. 紡織學(xué)報, 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.)

[8] 李京菁，趙美英．平紋機(jī)織復(fù)合材料靜力學(xué)性能分析[J]. 機(jī)械科學(xué)與技術(shù), 2015，34(8)：1294-1297．

(LI J J，ZHAO M Y. Mechanical analysis of plain woven composite materials[J]．Mechanical Science and Technology for Aerospace Engineering，2015，34(8)：1294-1297．)

[9] ROMEL T P, CUNNINGHAM P R. A multi-scale finite element approach for modelling damage progression in woven composite structures[J]. Composite Structures, 2012,94(3):977-986.

[10] 徐焜, 許希武. 三維編織復(fù)合材料彈性性能數(shù)值預(yù)測及細(xì)觀應(yīng)力分析[J]. 復(fù)合材料學(xué)報, 2007,24(3):178-185.

(XU K,XU X W. Prediction of elastic constants and simulation of stress field of 3D braided composites based on the finite element method[J]. Acta Materiae Compositae Sinica, 2007,24(3):178-185.)

[11] HASHIN Z. Failure criteria for unidirectional fiber composites[J]. Journal of Applied Mechanics, 1980,47(2):329-334.

[12] LINDE P, PLEITNER J, BOER H D,etal. Modelling and simulation of fibre metal laminates:ABAQUS User’s Conference[Z]. Boston,USA:Dassault Systemes Company,2004.

[13] XU Y, ZHANG W, BASSIR D. Stress analysis of multi-phase and multi-layer plain weave composite structure using global/local approach[J]. Composite Structures, 2010,92(5):1143-1154.

[14] 張增光. 2D-C/SiC復(fù)合材料的彈性常數(shù)預(yù)測及失效模式研究[D].西安:西北工業(yè)大學(xué), 2005.

(ZHANG Z G. Uniaxial macro-mechanical property and failure mode of a 2D-woven C/SiC composite[D].Xi′an:Northwestern Polytechnical University,2005.)

[15] 王丹勇. 層合板接頭損傷失效與疲勞壽命研究[D]. 南京：南京航空航天大學(xué), 2006.

(WANG D Y.Damage failure and fatigue life prediction in composite laminates[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2006.)

[16] 曾翔龍,王奇志. 二維機(jī)織C/SiC復(fù)合材料非線性力學(xué)行為數(shù)值模擬[J]. 宇航材料工藝, 2017(1):29-36.

(ZENG X L,WANG Q Z.Numerical simulation of nonlinear mechanics behavior for 2D weave composites[J]. Aerospace Materials & Technology,2017(1):29-36.)

[17] 張超, 許希武, 嚴(yán)雪. 紡織復(fù)合材料細(xì)觀力學(xué)分析的一般性周期性邊界條件及其有限元實現(xiàn)[J]. 航空學(xué)報, 2013,34(7):1636-1645.

(ZHANG C,XU X W,YAN X.General periodic boundary conditions and their application to micromechanical finite element analysis of textile composites[J].Acta Aeronautica et Astronautica Sinica,2013,34(7):1636-1645.)

[18] 李俊. 二維C/SiC復(fù)合材料的非線性本構(gòu)關(guān)系研究[D].西安：西北工業(yè)大學(xué), 2014.

(LI J. Research on the nonlinear constitutive relationship of 2D C/SiC composites[D].Xi′an:Northwestern Polytechnical University,2014.)

[19] 張立同. 纖維增韌碳化硅陶瓷復(fù)合材料[M].北京：化學(xué)工業(yè)出版社, 2009.

(ZHANG L T.Fiber-reinforced silicon carbide ceramic composites[M].Beijing:Chemical Industry Press,2009.)

(責(zé)任編輯：張 崢)

NumericalSimulationofC/SiCPlainWeaveCompositeswithDefectsunderUnidirectionalTension

ZENG Xianglong, WANG Qizhi, SU Fei

(School of Aeronautical Science and Engineering, Beihang University, Beijing 100191, China)

The second time develop of the ANSYS was realized by using APDL language, and an ovoid sectional multi-scale unit cell model for fiber bundle section with pre-made defects was established. At first the initial modulus, the strength and the ultimate strain of fiber bundle were calculated. Then the unit cell model of C/SiC woven composites was established according to the SEM photographs. By introducing the periodic boundary conditions, the initial anisotropic constants of C/SiC woven composites were predicted. The failure criteria proposed by Linde was utilized to set up the progressive damage mode,and the C/SiC plain weave composites with defects under unidirectional tension was simulated. So,the damage evolution process of the fiber bundle could be interpreted, when the unit cell of C/SiC plain weave composites was under the unidirection load. The numerical results show that the proposed model accurately captures the data from the experiments, which demonstrates the validity of the present analytical model. Furthermore, the numerical model provides an alternate way to design and predict the mechanical properties of plain weave composites.

woven C/SiC composites; progressive damage model; failure criteria; strength prediction

2016-05-17；

2016-07-11

王奇志(1964—)，女，博士，教授，主要從事復(fù)合材料力學(xué)、斷裂力學(xué)和損傷力學(xué)方面的研究，(E-mail)wqz@buaa.edu.cn。

10.11868/j.issn.1005-5053.2016.000086

TB332；O341

： A

： 1005-5053(2017)04-0061-08