高濤， 漆文凱,2， 沈承,2,*

1.南京航空航天大學(xué) 能源與動(dòng)力學(xué)院， 南京 210016 2.江蘇省航空動(dòng)力系統(tǒng)重點(diǎn)實(shí)驗(yàn)室， 南京 210016

基于一種新的均勻化實(shí)施方法的FRP剛度預(yù)測(cè)

高濤1， 漆文凱1,2， 沈承1,2,*

1.南京航空航天大學(xué) 能源與動(dòng)力學(xué)院， 南京 210016 2.江蘇省航空動(dòng)力系統(tǒng)重點(diǎn)實(shí)驗(yàn)室， 南京 210016

基于一種新的漸近均勻化(AH)實(shí)施方法，預(yù)測(cè)并討論了單向纖維增強(qiáng)復(fù)合材料(FRP)的宏觀等效彈性性能及纖維排列方式對(duì)等效力學(xué)性能的影響。該方法方便地將有限元分析(FEA)軟件作為一個(gè)工具箱使用，只需在單胞上施加簡(jiǎn)單位移周期邊界條件開展靜力學(xué)分析，即可經(jīng)過(guò)簡(jiǎn)單計(jì)算得到等效彈性常數(shù)，相比傳統(tǒng)均勻化實(shí)施方法顯著降低了實(shí)施難度并簡(jiǎn)化了計(jì)算過(guò)程。通過(guò)對(duì)比不同數(shù)值方法的結(jié)果驗(yàn)證了該方法的有效性和精確性。數(shù)值結(jié)果表明：六邊形排列下單向纖維增強(qiáng)復(fù)合材料呈現(xiàn)橫觀各向同性，而正方形排列下則呈現(xiàn)宏觀正交各向異性，經(jīng)過(guò)剛度平均化過(guò)程可得到橫觀各向異性材料性質(zhì)，纖維體積含量對(duì)兩種排列方式下材料等效彈性模量影響顯著但有所差別。

均勻化方法； 等效彈性性能； 周期邊界條件； 有限元； 纖維排列； 復(fù)合材料

隨著復(fù)合材料力學(xué)性能的優(yōu)越性和可設(shè)計(jì)性得到了廣泛的認(rèn)可，越來(lái)越多的輕質(zhì)高強(qiáng)復(fù)合材料被大量應(yīng)用于航空、航天和船舶領(lǐng)域，如：飛機(jī)的機(jī)翼、機(jī)身和螺旋槳，火箭發(fā)動(dòng)機(jī)殼體及船舶殼體等[1]。鑒于復(fù)合材料大多具有復(fù)雜的微細(xì)觀結(jié)構(gòu)，其宏觀力學(xué)性質(zhì)較難得到準(zhǔn)確的預(yù)測(cè)。若能在設(shè)計(jì)階段獲知其宏觀力學(xué)性能，將對(duì)復(fù)合材料的工程應(yīng)用以及設(shè)計(jì)開發(fā)具有重要的價(jià)值，因此復(fù)合材料的等效宏觀性能預(yù)測(cè)一直是材料及力學(xué)研究者關(guān)注的重要課題。

迄今為止，學(xué)者們提出了預(yù)測(cè)復(fù)合材料有效性能的許多模型和計(jì)算方法，詳見(jiàn)綜述文獻(xiàn)[2-3]。其中，自洽(SCS)方法[4]、廣義自洽(GSCS)方法[5]及Mori-Tanaka(M-T)方法[6]多被應(yīng)用于推導(dǎo)具有簡(jiǎn)單微結(jié)構(gòu)復(fù)合材料的數(shù)值解析公式[7-8]。而代表體元(Representative Volume Element，RVE)法[9]和漸近均勻化(Asymptotic Homogenization，AH)方法是近些年發(fā)展起來(lái)并得到廣泛應(yīng)用的數(shù)值計(jì)算方法，兩者均可處理具有復(fù)雜周期性微結(jié)構(gòu)的復(fù)合材料，且借助于有限元法便于實(shí)施和應(yīng)用。具體來(lái)說(shuō)，RVE法由于具有清晰的力學(xué)概念且應(yīng)用過(guò)程較為簡(jiǎn)單，常被應(yīng)用于復(fù)合材料以及點(diǎn)陣結(jié)構(gòu)等效彈性性能的預(yù)測(cè)。例如，Sun和Vaidya[10]基于三維代表體單元和周期邊界條件預(yù)測(cè)了復(fù)合材料的等效性能。Kanit等[11]提出了隨機(jī)分布復(fù)合材料代表體積單元的定量選擇模型。值得注意的是，RVE法并非基于嚴(yán)格的數(shù)學(xué)推導(dǎo)，因此一般僅適用于等效彈性性能的近似估計(jì)。

AH方法是一種基于攝動(dòng)理論的多尺度數(shù)值等效方法，它通過(guò)近似求解定義在單胞上的偏微分方程進(jìn)而得到材料的等效彈性性能[12]，該方法完整的推導(dǎo)過(guò)程可參考文獻(xiàn)[13-14]。邢譽(yù)峰和陳磊[15]研究了攝動(dòng)階次和單元階次對(duì)均勻化方法計(jì)算精度的影響，強(qiáng)調(diào)了二階攝動(dòng)對(duì)精度的作用。AH方法被廣泛應(yīng)用于求解各種材料和結(jié)構(gòu)的等效彈性性能及宏細(xì)觀應(yīng)力分析[16-18]。例如，Hassan等[19]基于AH方法得到了三維編織復(fù)合材料、層合板、薄壁殼以及蜂窩夾芯三明治殼等結(jié)構(gòu)的等效力學(xué)性能[20-23]。具體實(shí)施過(guò)程中，AH方法通常借助于有限元理論，需針對(duì)不同問(wèn)題建立相應(yīng)的有限元結(jié)構(gòu)單元，進(jìn)行大量編程和積分運(yùn)算，操作起來(lái)比較復(fù)雜且求解效率比較低。為此，Cheng等[12]提出了一種基于漸近均勻化理論的新型有限元實(shí)施(Novel Implementation of Asymptotic Homogenization， NIAH)方法，已經(jīng)被成功應(yīng)用于蜂窩板的等效剛度預(yù)測(cè)[24]。該方法可以避免大量編程和積分運(yùn)算，但實(shí)施過(guò)程中求解每個(gè)參數(shù)需對(duì)所有節(jié)點(diǎn)分3次加載等效初始位移場(chǎng)和周期性邊界條件，加載分析過(guò)程較繁瑣且運(yùn)算量大。因此，本文改進(jìn)了傳統(tǒng)實(shí)施方法，可以更方便地借助于商業(yè)有限元軟件作為工具箱來(lái)實(shí)施AH方法。通過(guò)引入一個(gè)新的特征函數(shù)(在Wang等的研究[25]中已有詳細(xì)證明)，將原AH實(shí)施方法中的初始應(yīng)變場(chǎng)和周期性邊界條件簡(jiǎn)化成簡(jiǎn)單的位移型周期邊界。新的實(shí)施方法借助于商業(yè)有限元軟件作為工具箱，只需對(duì)單胞模型的邊界施加簡(jiǎn)單位移型邊界條件，進(jìn)行一次靜力分析，即可直接得到等效彈性常數(shù)。該方法具有廣泛的適用性，可應(yīng)用于具有周期性微結(jié)構(gòu)的材料彈性性能預(yù)測(cè)，如纖維增強(qiáng)復(fù)合材料，格柵結(jié)構(gòu)材料和點(diǎn)陣結(jié)構(gòu)等。

眾所周知，隨機(jī)分布的長(zhǎng)纖維增強(qiáng)復(fù)合材料(FRP)具有橫觀各向同性的材料性質(zhì)。在以往的解析和數(shù)值研究中，通常利用正方形或六邊形排列模型來(lái)預(yù)測(cè)其等效力學(xué)性能，其中纖維六邊形排列模型可完全等效為橫觀各向同性，但正方形排列下的模型卻具有正交各向異性的材料性質(zhì)，這在以往的研究中通常被忽視。因此，深入分析排列方式對(duì)單向纖維增強(qiáng)復(fù)合材料等效力學(xué)性能的影響具有重要的理論價(jià)值。本文基于新的漸近均勻化實(shí)施方法，通過(guò)建立可考慮不同纖維排列方式的統(tǒng)一模型，系統(tǒng)研究了排列方式以及纖維體積含量對(duì)單向纖維增強(qiáng)復(fù)合材料等效彈性參數(shù)的影響。

1 新的漸近均勻化方法

(1)

根據(jù)彈性力學(xué)基本理論，應(yīng)變與位移、應(yīng)力與應(yīng)變的關(guān)系可以表示為

(2)

式中：εkl和σij為應(yīng)變和應(yīng)力張量；Eijkl為彈性常數(shù)張量。

根據(jù)虛位移原理，虛位移方程可以描述為

(3)

圖1 周期結(jié)構(gòu)及其單胞示意圖 Fig.1 Periodic stucture and corresponding unit cell

將式(1)代入式(3)，線彈性單胞問(wèn)題可以描述為[13]

(4)

對(duì)于如圖2所示的長(zhǎng)方體單胞，周期邊界條件可以描述為

(5)

式中：U?代表單胞沿y1正方向上的兩個(gè)面；V?代表單胞沿y2正方向上的兩個(gè)面；W?代表單胞沿y3正方向上的兩個(gè)面。

等效彈性常數(shù)可由式(6)得到：

(6)

或?qū)懗删仃囆问綖?/p>

(7)

圖2 長(zhǎng)方體單胞 Fig.2 A rectangular parallelepiped unit cell

(8)

對(duì)于三維單胞問(wèn)題，有

(9)

(10)

式中：δkn和δln為Kronecker函數(shù)，n取1，2，3。

通過(guò)引入這個(gè)特征函數(shù)，等效彈性常數(shù)式(6)可表示為

(11)

通過(guò)引入新的特征函數(shù)張量，原來(lái)的對(duì)單胞加載初始應(yīng)變和周期邊界條件的問(wèn)題轉(zhuǎn)化成簡(jiǎn)單的位移周期邊界問(wèn)題。周期性邊界條件式(5)可以進(jìn)一步簡(jiǎn)化為

(12)

結(jié)合式(5)和式(12)，簡(jiǎn)單位移型周期邊界條件可以描述為6種初始位移邊界加載，即

(13)

(14)

(15)

(16)

(17)

(18)

可以看出，這6種位移邊界條件格式統(tǒng)一，只與單胞的尺寸有關(guān)，故加載更加簡(jiǎn)便。

根據(jù)文獻(xiàn)[26]，式(11)可進(jìn)一步表示為

(19)

矩陣形式為

(20)

(21)

2 新方法實(shí)施流程簡(jiǎn)述

圖3 傳統(tǒng)AH方法和新實(shí)施方法流程圖 Fig.3 Flowcharts of traditional AH and new implementation methods

為了更清楚地展現(xiàn)新的方法，以及新實(shí)施方法和傳統(tǒng)AH方法的不同和聯(lián)系，繪制兩者的流程圖，如圖3所示. 從對(duì)比流程圖中可以看到，傳統(tǒng)AH求解方法中，需要在每個(gè)單元上積分，這就要求必須清楚有限元單元的每一個(gè)細(xì)節(jié)，如形函數(shù)、本構(gòu)關(guān)系和應(yīng)變-位移矩陣。對(duì)于不同的單元類型，這些矩陣也都不同。因此必須針對(duì)不同單胞和單元自編不同的有限元模型、列式和程序并進(jìn)行多次求解。在新的實(shí)施方法中，可以用商業(yè)有限元軟件來(lái)完成建模和一次直接求解，比如ANSYS。新的方法實(shí)施起來(lái)更為簡(jiǎn)單方便，只需編寫一個(gè)小的腳本，針對(duì)不同單胞不同單元，都可以直接計(jì)算，不用考慮具體的細(xì)節(jié)。商業(yè)軟件作為一個(gè)工具箱，僅需要一次靜力加載即可直接輸出結(jié)果，大大簡(jiǎn)化漸進(jìn)均勻化方法的實(shí)施流程，同時(shí)，無(wú)需自編程序和積分運(yùn)算，提高了計(jì)算效率。在實(shí)施求解時(shí)也應(yīng)注意：① 為了保證周期性邊界條件的準(zhǔn)確施加，有限元單胞模型在建模時(shí)需要在平行的兩個(gè)邊界面對(duì)應(yīng)節(jié)點(diǎn)的保持對(duì)稱性；② 若單胞模型采用二維單元或者一維單元(如殼或梁?jiǎn)卧?，施加位移周期邊界時(shí)需要將施加的節(jié)點(diǎn)平動(dòng)自由度轉(zhuǎn)換成轉(zhuǎn)角自由度。

根據(jù)復(fù)合材料力學(xué)，橫觀各向同性材料的剛度系數(shù)與彈性模量的關(guān)系有

S=D-1=

(22)

式中:下標(biāo)a和t分別代表縱向(即沿纖維方向)和橫向(即垂直于纖維方向)；E、G和ν分別為彈性模量、剪切模量和泊松比。

本文基于商業(yè)有限元軟件ANSYS建立單胞模型，為便于施加周期邊界條件，單胞相對(duì)邊界網(wǎng)格需保證完全一致，然后運(yùn)用APDL語(yǔ)言可方便地施加初始位移場(chǎng)。等效剛度系數(shù)矩陣可由ANSYS結(jié)果直接輸出，通過(guò)橫觀各向同性材料的剛度系數(shù)矩陣與彈性模量的關(guān)系式(22)最終得到材料等效彈性模量。

3 算例分析

3.1 有效性驗(yàn)證

圖5為該單胞模型的有限元網(wǎng)格，基體和纖維結(jié)合部經(jīng)適當(dāng)加密。模型采用網(wǎng)格數(shù)為24 420的Solid185實(shí)體單元進(jìn)行有限元離散。本例假設(shè)為硼/鋁復(fù)合材料，纖維和基體材料參數(shù)為：Ef=379.3 GPa，νf=0.1；Em=68.3 GPa，νm=0.3,Vm=0.47。

圖4 單向纖維增強(qiáng)復(fù)合材料及對(duì)應(yīng)單胞示意圖 Fig.4 Illustration of unidirectional fibre reinforced composite and corresponding unit cell

圖5 單胞的有限元網(wǎng)格 Fig.5 Finite element mesh of unit cell

本文模型和其他模型對(duì)比結(jié)果見(jiàn)表1，其中，F(xiàn)EM[10]實(shí)際是一種基于有限元理論的數(shù)值逼近方法。GMC方法[27]、HFGMC方法[27]和ECM方法[28]是MOC方法[29]的變形，本質(zhì)是利用Legendre多項(xiàng)式的不同階形式對(duì)不同尺度下的位移進(jìn)行展開。VAMUCH方法[30]是另外一種漸近展開法，它通過(guò)對(duì)單胞能量泛函漸進(jìn)展開來(lái)預(yù)測(cè)有效彈性常數(shù)。從表1中可以看出，MOC和GMC預(yù)測(cè)的剪切模量G12偏小，而FEM方法預(yù)測(cè)的G12偏大，本文和另外一種漸進(jìn)均勻化方法VAMUCH方法以及HFGMC方法計(jì)算的結(jié)果較為接近。諸多方法中，本文方法實(shí)施最為簡(jiǎn)單有效，且結(jié)果較精確，充分說(shuō)明了該方法的有效性。

表1不同方法預(yù)測(cè)的硼/鋁復(fù)合材料等效材料性能

Table1Equivalentmaterialpropertiesofbron/aluminumcompositesobtainedusingdifferentmethods

MethodE11/GPaE22/GPaG12/GPaG23/GPaν12ν23VAMUCH[30]215．3144．154．345．90．1950．255FEM[10]21514457．245．90．1950．29MOC[29]215142．651．343．70．200．25GMC[27]21514151．243．700．1950．261HFGMC[27]215．414454．345．830．1950．255ECM[28]215143．454．345．10．1950．26Present215．2143．954．345．810．1950．255

Notes：Ef=379.3 GPa,νf=0.1,Em=68.3 GPa,νm=0.3,Vm=0.47.

3.2 不同單胞模型對(duì)結(jié)果的影響

圖6給出了同一纖維排列下選取不同單胞模型的示意圖，表2列出了C1～C4相同纖維體積分?jǐn)?shù)下不同單胞模型的彈性模量及泊松比結(jié)果，C4是經(jīng)過(guò)坐標(biāo)變換后的結(jié)果。結(jié)果顯示：C1、C2和C3的結(jié)果完全一致，只有C4彈性模量有細(xì)小差別，這種細(xì)小的差異是由有限元網(wǎng)格本身的差異造成的。這說(shuō)明同一纖維排列方式下，彈性模量的結(jié)果與單胞模型的選取無(wú)關(guān)。

圖6 不同的單胞模型 Fig.6 Different unit cell models

表2 不同單胞模型的結(jié)果Table 2 Results using different unit cell models

ModelofunitcellEa/GPaEt/GPaGt/GPaGa/GPaνaνtC137．9410．932．902．190．0870．410C237．9410．932．902．190．0870．410C337．9410．932．902．190．0870．410C437．9310．922．902．190．0870．410

Notes：Ef=73 GPa,νf=0.22,Em=2.9 GPa,νm=0.4,Vm=0.5.

3.3 正方形排列與六邊形排列

圖7 菱形纖維排列下復(fù)合材料的橫截面和單胞 Fig.7 Cross-section and unit cell of composites with rhombic fiber arrangement

圖8 單胞模型的有限元網(wǎng)格 Fig.8 Finite element meshes of unit cell model

圖9 正方形纖維排列下的坐標(biāo)變換 Fig.9 Coordinate transformation of square fiber arrangements

D′=T(θ)DTT(θ)

(23)

式中：T(θ)為坐標(biāo)轉(zhuǎn)換矩陣，且

(24)

其中：li,mi,ni(i=1,2,3)為方向余弦。

材料參數(shù)同表2。計(jì)算可得到等效剛度矩陣，再根據(jù)式(23)的坐標(biāo)變換，得到兩種纖維排列方式下等效剛度系數(shù)隨θ的變化規(guī)律(如圖10所示)。結(jié)果表明，兩種纖維排列下任意θ(-45°≤θ≤45°)均滿足如下關(guān)系：D11=D22,D13=D23,D44=D55,D16=-D26。纖維排列方式為六邊形時(shí)，所有剛度系數(shù)都不隨θ變化，說(shuō)明該纖維排列下模型呈現(xiàn)為橫觀各向同性材料屬性。纖維呈正方形排列時(shí)，模型呈現(xiàn)正交各向異性，只有D13(=D23),D33,D44(=D55)不隨θ變化，并且其拉伸-彎曲耦合剛度系數(shù)D16(=-D26)在某些角度θ下并不為零。兩種排列方式下剛度系數(shù)D13(=D23),D33,D44(=D55)非常接近。

為了消除轉(zhuǎn)角對(duì)有效彈性常數(shù)的影響，令剛度系數(shù)在所有可能轉(zhuǎn)角θ下取平均值，此時(shí)平均剛度系數(shù)可以表示為

(25)

式中：Dsquare(θ)為角度θ下的剛度矩陣。此時(shí)，正方形排列下平均化后的剛度矩陣為

(26)

平均化后的剛度矩陣不再隨轉(zhuǎn)角θ變化，進(jìn)而可得到橫觀各向同性材料性質(zhì)下的剛度矩陣，最終由式(22)獲得等效彈性模量。

圖10 兩種纖維排列方式下的等效剛度系數(shù)隨θ的變化 Fig.10 Variations of effective stiffness coefficients with θ for two fiber arrangements

3.4 纖維體積分?jǐn)?shù)對(duì)彈性模量的影響

圖11給出了正方形和六邊形纖維排列下的拉伸和剪切模量隨纖維體積分?jǐn)?shù)變化曲線，其中，正方形排列的等效模量是經(jīng)過(guò)剛度平均化后的結(jié)果。從圖中可以看出：兩種排列方式在相同體積含量下的縱向拉伸模量(Ea)相同，且都隨纖維體積分?jǐn)?shù)的增加呈線性增加。說(shuō)明等效的縱向拉伸模量主要受纖維縱向拉伸模量的影響，與纖維體積分?jǐn)?shù)有關(guān)，與纖維排列方式無(wú)關(guān)。在較低纖維體積含量(Vf≤0.35)下，正方形和六邊形纖維排列方式下的橫向拉伸模量(Et)、橫向剪切模量(Gt)和縱向剪切模量(Ga)隨纖維體積含量的增大緩慢增加，且兩種排列下的值都非常接近，而在較高纖維體積含量(Vf>0.35)下，這些模量隨纖維體積分?jǐn)?shù)增大的速率不斷增大，正方形排列下的模量都比六邊形排列要大，且這種差距隨纖維體積分?jǐn)?shù)的增加而逐漸增大。總的來(lái)看，相同體積分?jǐn)?shù)條件下，除了縱向拉伸模量外，正方形排列材料的拉伸和剪切彈性模量均大于六邊形排列材料的。說(shuō)明在較低纖維體積含量時(shí)，其橫向拉伸模量、橫向剪切模量和縱向剪切模量受基體材料的模量影響較大，在較高纖維體積含量時(shí)，這些模量受纖維材料的模量影響較大。物理機(jī)理[31]在于正方形排列纖維之間的最小間距比六邊形排列時(shí)要小，即纖維排列更緊密，因此其等效橫向拉伸模量和剪切模量更大。

圖11 不同體積分?jǐn)?shù)下的拉伸模量和剪切模量 Fig.11 Tensile modulus and shear modulus under different volume fraction

4 結(jié) 論

1) 基于一種新型漸近均勻化實(shí)施方法，建立了單向纖維增強(qiáng)復(fù)合材料的等效彈性性能預(yù)測(cè)模型。通過(guò)和其他理論對(duì)比，驗(yàn)證了該方法的準(zhǔn)確性和簡(jiǎn)便性。

2) 數(shù)值結(jié)果表明：?jiǎn)伟倪x取方式對(duì)求解結(jié)果沒(méi)有影響；通過(guò)建立不同纖維排列方式下的統(tǒng)一單胞模型，可討論排列方式對(duì)材料等效參數(shù)的影響，六邊形排列下的單胞模型具有橫觀各向同性的材料性質(zhì)，而正方形排列下的單胞模型則具有正交各向異性的性質(zhì)，此時(shí)，可通過(guò)坐標(biāo)變換和剛度平均化由正方形排列單胞模型得到橫觀各向同性的材料性質(zhì)；此外，纖維體積含量對(duì)不同排列方式下材料的拉伸及剪切模量影響顯著但有所差別。

[1] 杜善義, 章繼峰, 張博明. 先進(jìn)復(fù)合材料格柵結(jié)構(gòu)(AGS)應(yīng)用與研究進(jìn)展[J]. 航空學(xué)報(bào), 2007, 28(2): 419-424.

DU S Y, ZHANG J F, ZHANG B M. Overview of application and research on advanced composite grid structures[J]. Acta Aeronautica et Astronautica Sinica, 2007, 28(2): 419-424 (in Chinese).

[2] DIXIT A, MALI H S. Modeling techniques for predicting the mechanical properties of woven-fabric textile composites: A review[J]. Mechanics of Composite Materials, 2013, 49(1): 1-20.

[3] YOUNES R, HALLAL A, FARDOUN F, et al. Comparative review study on elastic properties modeling for unidirectional composite materials[M]//Composites and Their Properties. Rijeka: InTech, 2012: 261-290.

[4] HILL R. A self-consistent mechanics of composite materials[J]. Journal of the Mechanics and Physics of Solids, 1965, 13(4): 213-222.

[5] CHRISTENSEN R M, LO K H. Solutions for effective shear properties in three phase sphere and cylinder models[J]. Journal of the Mechanics and Physics of Solids, 1979, 27(4): 315-330.

[6] MORI T, TANAKA K. Average stress in matrix and average elastic energy of materials with misfitting inclusions[J]. Acta Metallurgica, 1973, 21(5): 571-574.

[7] LIU X N, HU G K. A continuum micromechanical theory of overall plasticity for particulate composites including particle size effect[J]. International Journal of Plasticity, 2005, 21(4): 777-799.

[8] HU N, WANG B, TAN G W, et al. Effective elastic properties of 2-D solids with circular holes: Numerical simulations[J]. Composites Science and Technology, 2000, 60(9): 1811-1823.

[9] LIU T, DENG Z C, LU T J. Design optimization of truss-cored sandwiches with homogenization[J]. International Journal of Solids & Structures, 2006, 43(25-26): 7891-7918.

[10] SUN C T, VAIDYA R S. Prediction of composite properties from a representative volume element[J]. Composites Science and Technology, 1996, 56(2): 171-179.

[11] KANIT T, FOREST S, GALLIET I, et al. Determination of the size of the representative volume element for random composites: Statistical and numerical approach[J]. International Journal of Solids and Structures, 2003, 40(13-14): 3647-3679.

[12] CHENG G D, CAI Y W, XU L. Novel implementation of homogenization method to predict effective properties of periodic materials[J]. Acta Mechanica Sinica, 2013, 29(4): 550-556.

[13] BENSSOUSAN A, LIONS J L, PAPANICOLAOU G. Asymptotic analysis for periodic structures[M]. New York: North-Holland Publishing Co., 1978: 102-121.

[14] HASSANIA B, HINTON E. A review of homogenization and topology optimization I—Homogenization theory for media with periodic structure[J]. Computers and Structures, 1998, 69(6): 707-717.

[15] 邢譽(yù)峰, 陳磊. 若干周期性復(fù)合材料結(jié)構(gòu)數(shù)學(xué)均勻化方法的計(jì)算精度[J]. 航空學(xué)報(bào), 2015, 36(5): 1520-1529.

XING Y F, CHEN L. Accuracy analysis of mathematical homogenization method for several periodical composite structure[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(5): 1520-1529 (in Chinese).

[16] ANDRIANOV I V, DANISHEVS’KYY V V, WEICHERT D. Asymptotic determination of effective elastic properties of composite materials with fibrous square-shaped inclusions[J]. European Journal of Mechanics—A/Solids, 2002, 21(6): 1019-1036.

[17] ASKARI D, VEEDU V P, GHASEMI-NEJHAD M N. Chirality dependence of carbon single-walled nanotube material properties: Axial coefficient of thermal expansion[J]. Journal of Nanoscience & Nanotechnology, 2006, 6(7): 2167-2174.

[18] NASUTION M R E, WATANABE N, KONDO A, et al. Thermomechanical properties and stress analysis of 3-D textile composites by asymptotic expansion homogenization method[J]. Composites Part B: Engineering, 2014, 60: 378-391.

[19] HASSAN E M, GEORGIADES A V, SAVI M A, et al. Analytical and numerical analysis of 3D grid-reinforced orthotropic composite structures[J]. International Journal of Engineering Science, 2011, 49(7): 589-605.

[20] KALAMKAROV A L, GEORGIADES A V. Asymptotic homogenization models for smart composite plates with rapidly varying thickness: Part I—Theory[J]. International Journal for Multiscale Computational Engineering, 2004, 2(1): 133-148.

[21] CHALLAGULLA K S, GEORGIADES A V, KALAMKAROV A L. Asymptotic homogenization modeling of smart composite generally orthotropic grid-reinforced shells: Part I—Theory[J]. European Journal of Mechanics—A/Solids, 2010, 29(4): 530-540.

[22] GEORGIADES A V, CHALLAGULLA K S, KALAMKAROV A L. Asymptotic homogenization modeling of smart composite generally orthotropic grid-reinforced shells: Part II—Applications[J]. European Journal of Mechanics—A/Solids, 2010, 29(4): 541-556.

[23] KALAMKAROV A L, ANDRIANOV I V, DANISHEVS’KYY V V. Asymptotic homogenization of composite materials and structures[J]. Applied Mechanics Reviews, 2009, 62(3): 669-676.

[24] CAI Y W, XU L, CHENG G D. Novel numerical implementation of asymptotic homogenization method for periodic plate structures[J]. International Journal of Solids and Structures, 2014, 51(1): 284-292.

[25] WANG W X, LUO D M, TAKAO Y, et al. New solution method for homogenization analysis and its application to the prediction of macroscopic elastic constants of materials with periodic microstructures[J]. Computers and Structures, 2006, 84(15-16): 991-1001.

[26] SIGMUND O, TORQUATO S. Design of smart composite materials using topology optimization[J]. Smart Materials and Structures, 1999, 8(3): 365-379.

[27] ABOUDI J, PINDERA M J, ARNOLD S M. Linear thermoelastic higher-order theory for periodic multiphase materials[J]. Journal of Applied Mechanics, 2001, 68(5): 697-707.

[28] WILLIAMS T O. A two-dimensional, higher-order, elasticity-based micromechanics model[J]. International Journal of Solids and Structures, 2005, 42(3-4): 1009-1038.

[29] ABOUDI J. A continuum theory for fiber-reinforced elastic-viscoplastic composites[J]. International Journal of Engineering Science, 1982, 20(5): 605-621.

[30] YU W B, TANG T. Variational asymptotic method for unit cell homogenization of periodically heterogeneous materials[J]. International Journal of Solids and Structures, 2007, 44(11-12): 3738-3755.

[31] 許震宇， 張若京. 纖維排列方式對(duì)單向纖維加強(qiáng)復(fù)合材料彈性常數(shù)的影響[J]. 力學(xué)季刊, 2002, 23(1): 59-63.

XU Z Y, ZHANG R J. Influence of fibers’ packing modes on elastic constants of fiber reinforced polymer (FRP) material[J]. Chinese Quarterly of Mechanics, 2002, 23(1): 59-63 (in Chinese).

(責(zé)任編輯： 徐曉)

URL：www.cnki.net/kcms/detail/11.1929.V.20170217.1614.012.html

PredictionofFRPstiffnessbasedonanewimplementationofhomogenizationmethod

GAOTao1,QIWenkai1,2,SHENCheng1,2,*

1.CollegeofEnergyandPowerEngineering,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China2.JiangsuProvinceKeyLaboratoryofAerospacePowerSystem,Nanjing210016,China

Basedonanewimplementationofasymptotichomogenization(AH)method,themacroscopicequivalentelasticpropertiesofunidirectionalfiberreinforcedpolymer(FRP)andtheeffectoffiberarrangementontheequivalentmechanicalpropertieswerepredictedanddiscussed.Thismethodcanbeeasilyimplementedtoobtaintheeffectiveelasticconstantsusingcommercialfiniteelementanalysis(FEA)softwareasatoolbox,andnodaldisplacementfieldcorrespondingtotheunitstrainfieldwasappliedwithconsiderationofperiodicboundaryconditions.Comparedtotraditionalhomogenizationmethodsthismethodsignificantlyreducedthedifficultyofimplementationandsimplifiedthecalculationprocess.Acomparisonofseveralmethodsrevealsthatthisimplementationmethodissimpleandtheresultobtainedisaccurate.ThenumericalresultsshowfiberreinforcedpolymerwithhexagonalarrangementhastransverselyisotropicmaterialbehaviorbutFRPwithsquarearrangementhasorthotropicmaterialbehavior;astiffnessaveragenessprocedureisrequiredtoobtainthetransverselyisotropicstiffnessmatrixofFRP;thefibervolumefractionhasasignificantbutdifferenteffectontheequivalentelasticmodulusofthetwokindsofarrangement.

homogenizationmethod;equivalentelasticproperties;periodicboundarycondition;finiteelement;fiberarrangement;composites

2016-06-30；Revised2017-01-18；Accepted2017-01-23；Publishedonline2017-02-171614

s：AeronauticalScienceFoundationofChina(2013ZB52019);NationalNaturalScienceFoundationofChina(11502110);NaturalScienceFoundationofJiangsuProvinceofChina(BK20150737);theFundamentalResearchFundsforCentralUniversities(NJ20150005)

.E-mailcshen@nuaa.edu.cn

2016-06-30；退修日期2017-01-18；錄用日期2017-01-23； < class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間

時(shí)間：2017-02-171614

www.cnki.net/kcms/detail/11.1929.V.20170217.1614.012.html

航空科學(xué)基金 (2013ZB52019)； 國(guó)家自然科學(xué)基金 (11502110)； 江蘇省自然科學(xué)基金 (BK20150737)； 中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金 (NJ20150005)

.E-mailcshen@nuaa.edu.cn

高濤， 漆文凱， 沈承. 基于一種新的均勻化實(shí)施方法的FRP剛度預(yù)測(cè)J． 航空學(xué)報(bào),2017,38(5):220579.GAOT,QIWK,SHENC.PredictionofFRPstiffnessbasedonanewimplementationofhomogenizationmethodJ.ActaAeronauticaetAstronauticaSinica,2017,38(5):220579.

http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

10.7527/S1000-6893.2017.220579

V258+.3

A

1000-6893(2017)05-220579-11