• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-11-20 03:43:46高濤漆文凱沈承
    航空學(xué)報 2017年5期
    關(guān)鍵詞:單胞模量正方形

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

    1.南京航空航天大學(xué) 能源與動力學(xué)院, 南京 210016 2.江蘇省航空動力系統(tǒng)重點實驗室, 南京 210016

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

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

    1.南京航空航天大學(xué) 能源與動力學(xué)院, 南京 210016 2.江蘇省航空動力系統(tǒng)重點實驗室, 南京 210016

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

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

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

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

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

    眾所周知,隨機分布的長纖維增強復(fù)合材料(FRP)具有橫觀各向同性的材料性質(zhì)。在以往的解析和數(shù)值研究中,通常利用正方形或六邊形排列模型來預(yù)測其等效力學(xué)性能,其中纖維六邊形排列模型可完全等效為橫觀各向同性,但正方形排列下的模型卻具有正交各向異性的材料性質(zhì),這在以往的研究中通常被忽視。因此,深入分析排列方式對單向纖維增強復(fù)合材料等效力學(xué)性能的影響具有重要的理論價值。本文基于新的漸近均勻化實施方法,通過建立可考慮不同纖維排列方式的統(tǒng)一模型,系統(tǒ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),線彈性單胞問題可以描述為[13]

    (4)

    對于如圖2所示的長方體單胞,周期邊界條件可以描述為

    (5)

    式中:U?代表單胞沿y1正方向上的兩個面;V?代表單胞沿y2正方向上的兩個面;W?代表單胞沿y3正方向上的兩個面。

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

    (6)

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

    (7)

    圖2 長方體單胞 Fig.2 A rectangular parallelepiped unit cell

    (8)

    對于三維單胞問題,有

    (9)

    (10)

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

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

    (11)

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

    (12)

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

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

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

    根據(jù)文獻[26],式(11)可進一步表示為

    (19)

    矩陣形式為

    (20)

    (21)

    2 新方法實施流程簡述

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

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

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

    S=D-1=

    (22)

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

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

    3 算例分析

    3.1 有效性驗證

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

    圖4 單向纖維增強復(fù)合材料及對應(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

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

    表1不同方法預(yù)測的硼/鋁復(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 不同單胞模型對結(jié)果的影響

    圖6給出了同一纖維排列下選取不同單胞模型的示意圖,表2列出了C1~C4相同纖維體積分數(shù)下不同單胞模型的彈性模量及泊松比結(jié)果,C4是經(jīng)過坐標(biāo)變換后的結(jié)果。結(jié)果顯示:C1、C2和C3的結(jié)果完全一致,只有C4彈性模量有細小差別,這種細小的差異是由有限元網(wǎng)格本身的差異造成的。這說明同一纖維排列方式下,彈性模量的結(jié)果與單胞模型的選取無關(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ù)式(23)的坐標(biāo)變換,得到兩種纖維排列方式下等效剛度系數(shù)隨θ的變化規(guī)律(如圖10所示)。結(jié)果表明,兩種纖維排列下任意θ(-45°≤θ≤45°)均滿足如下關(guān)系:D11=D22,D13=D23,D44=D55,D16=-D26。纖維排列方式為六邊形時,所有剛度系數(shù)都不隨θ變化,說明該纖維排列下模型呈現(xiàn)為橫觀各向同性材料屬性。纖維呈正方形排列時,模型呈現(xiàn)正交各向異性,只有D13(=D23),D33,D44(=D55)不隨θ變化,并且其拉伸-彎曲耦合剛度系數(shù)D16(=-D26)在某些角度θ下并不為零。兩種排列方式下剛度系數(shù)D13(=D23),D33,D44(=D55)非常接近。

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

    (25)

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

    (26)

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

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

    3.4 纖維體積分數(shù)對彈性模量的影響

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

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

    4 結(jié) 論

    1) 基于一種新型漸近均勻化實施方法,建立了單向纖維增強復(fù)合材料的等效彈性性能預(yù)測模型。通過和其他理論對比,驗證了該方法的準(zhǔn)確性和簡便性。

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

    [1] 杜善義, 章繼峰, 張博明. 先進復(fù)合材料格柵結(jié)構(gòu)(AGS)應(yīng)用與研究進展[J]. 航空學(xué)報, 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] 邢譽峰, 陳磊. 若干周期性復(fù)合材料結(jié)構(gòu)數(shù)學(xué)均勻化方法的計算精度[J]. 航空學(xué)報, 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] 許震宇, 張若京. 纖維排列方式對單向纖維加強復(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ò)出版時間

    時間:2017-02-171614

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

    航空科學(xué)基金 (2013ZB52019); 國家自然科學(xué)基金 (11502110); 江蘇省自然科學(xué)基金 (BK20150737); 中央高校基本科研業(yè)務(wù)費專項資金 (NJ20150005)

    .E-mailcshen@nuaa.edu.cn

    高濤, 漆文凱, 沈承. 基于一種新的均勻化實施方法的FRP剛度預(yù)測J. 航空學(xué)報,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

    猜你喜歡
    單胞模量正方形
    基于NURBS的點陣材料參數(shù)化建模方法
    復(fù)合材料周期結(jié)構(gòu)數(shù)學(xué)均勻化方法的一種新型單胞邊界條件
    剪正方形
    高勁度模量瀝青混合料在京臺高速車轍維修段的應(yīng)用
    剪拼正方形
    室內(nèi)回彈模量和回彈再壓縮模量試驗參數(shù)探討
    山西建筑(2020年11期)2020-06-04 00:09:48
    拼正方形
    拼正方形
    考慮界面層影響的三維機織復(fù)合材料單胞模型研究
    關(guān)于現(xiàn)行規(guī)范路基頂面回彈模量的理解和應(yīng)用
    上海公路(2018年4期)2018-03-21 05:57:24
    成人国产av品久久久| 啦啦啦在线观看免费高清www| 高清视频免费观看一区二区| 精品一区二区三区视频在线| 国产精品人妻久久久影院| 99精国产麻豆久久婷婷| 国产精品麻豆人妻色哟哟久久| 午夜福利视频精品| 久久久久久久久大av| 亚洲成人中文字幕在线播放| 久久久精品94久久精品| 国产精品不卡视频一区二区| 午夜老司机福利剧场| 26uuu在线亚洲综合色| 久久鲁丝午夜福利片| 尤物成人国产欧美一区二区三区| 午夜激情久久久久久久| 中文字幕av成人在线电影| 国产日韩欧美亚洲二区| 亚洲精品中文字幕在线视频 | 少妇 在线观看| 大片电影免费在线观看免费| 中文字幕av成人在线电影| 亚洲精华国产精华液的使用体验| 极品教师在线视频| 日韩免费高清中文字幕av| 免费av中文字幕在线| 亚洲精品乱码久久久久久按摩| 久久精品国产亚洲av涩爱| 内地一区二区视频在线| 亚洲av免费高清在线观看| 丰满迷人的少妇在线观看| 3wmmmm亚洲av在线观看| 最近中文字幕2019免费版| av视频免费观看在线观看| 亚洲精品国产av蜜桃| 我要看黄色一级片免费的| 欧美少妇被猛烈插入视频| 久久精品久久久久久噜噜老黄| 91午夜精品亚洲一区二区三区| 色综合色国产| 久热这里只有精品99| 中文字幕精品免费在线观看视频 | 亚洲精品乱码久久久v下载方式| 久久久精品94久久精品| 亚洲激情五月婷婷啪啪| 久久婷婷青草| 国产成人aa在线观看| 美女xxoo啪啪120秒动态图| 日本-黄色视频高清免费观看| 一个人免费看片子| 春色校园在线视频观看| 99久久中文字幕三级久久日本| 国产熟女欧美一区二区| 涩涩av久久男人的天堂| 一区二区三区免费毛片| 一级毛片 在线播放| 国产永久视频网站| 人妻制服诱惑在线中文字幕| av黄色大香蕉| 久久午夜福利片| 在线观看免费日韩欧美大片 | 国产成人精品一,二区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品国产av在线观看| 亚洲精品成人av观看孕妇| 免费观看性生交大片5| 99久国产av精品国产电影| 尤物成人国产欧美一区二区三区| 青春草亚洲视频在线观看| 深夜a级毛片| 你懂的网址亚洲精品在线观看| 老熟女久久久| 亚洲国产最新在线播放| 国产熟女欧美一区二区| 欧美成人精品欧美一级黄| 少妇猛男粗大的猛烈进出视频| 免费大片18禁| 九色成人免费人妻av| 成人黄色视频免费在线看| 国产男女超爽视频在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲色图综合在线观看| 国产一区二区在线观看日韩| 一级毛片aaaaaa免费看小| 久久久久性生活片| 视频中文字幕在线观看| 免费高清在线观看视频在线观看| 大又大粗又爽又黄少妇毛片口| 久热久热在线精品观看| 久久国内精品自在自线图片| 黑人猛操日本美女一级片| 久久99精品国语久久久| 亚洲无线观看免费| 精品人妻一区二区三区麻豆| 美女福利国产在线 | videossex国产| 色婷婷av一区二区三区视频| 中文天堂在线官网| 亚洲精品一二三| h视频一区二区三区| av线在线观看网站| 欧美日韩亚洲高清精品| 久久久久性生活片| 亚洲精品日本国产第一区| 最近手机中文字幕大全| 狂野欧美激情性xxxx在线观看| 国产国拍精品亚洲av在线观看| 欧美bdsm另类| 精品国产乱码久久久久久小说| 国产亚洲5aaaaa淫片| 亚洲激情五月婷婷啪啪| 亚洲av日韩在线播放| 国产精品久久久久久久电影| 成人亚洲精品一区在线观看 | 国产精品久久久久久久电影| 国产精品久久久久久久电影| 你懂的网址亚洲精品在线观看| 欧美最新免费一区二区三区| h日本视频在线播放| 国产白丝娇喘喷水9色精品| 黄色一级大片看看| 久久鲁丝午夜福利片| 国产欧美亚洲国产| 久久鲁丝午夜福利片| 日韩,欧美,国产一区二区三区| 亚洲,欧美,日韩| 国产成人午夜福利电影在线观看| 日韩视频在线欧美| 亚洲av成人精品一区久久| 丝袜喷水一区| 大又大粗又爽又黄少妇毛片口| 99热这里只有是精品50| 国产日韩欧美在线精品| 中文字幕亚洲精品专区| 成年av动漫网址| 免费黄网站久久成人精品| 亚洲图色成人| 国产久久久一区二区三区| 亚洲图色成人| 国产亚洲av片在线观看秒播厂| 国产色爽女视频免费观看| 精品久久久久久久末码| av免费观看日本| 国产精品无大码| 91aial.com中文字幕在线观看| 国产精品福利在线免费观看| 免费人妻精品一区二区三区视频| 男女啪啪激烈高潮av片| 亚洲精品日韩在线中文字幕| 亚洲真实伦在线观看| 国产精品偷伦视频观看了| 欧美最新免费一区二区三区| a 毛片基地| 国产精品久久久久久av不卡| 美女xxoo啪啪120秒动态图| 国产白丝娇喘喷水9色精品| 免费看光身美女| 夫妻午夜视频| 久久精品国产自在天天线| 国产免费福利视频在线观看| 岛国毛片在线播放| 成人高潮视频无遮挡免费网站| 久久久久国产网址| 国产精品人妻久久久久久| 又黄又爽又刺激的免费视频.| 国产精品久久久久久久久免| 久久久久久久大尺度免费视频| 国产无遮挡羞羞视频在线观看| 国产黄色视频一区二区在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品精品国产色婷婷| 有码 亚洲区| 大码成人一级视频| 高清av免费在线| 免费久久久久久久精品成人欧美视频 | 夜夜看夜夜爽夜夜摸| 在线播放无遮挡| av在线老鸭窝| 精品亚洲乱码少妇综合久久| 热re99久久精品国产66热6| 一级av片app| 成人二区视频| 黄色视频在线播放观看不卡| 久久国内精品自在自线图片| 欧美日韩国产mv在线观看视频 | 国产色爽女视频免费观看| 国产91av在线免费观看| 高清在线视频一区二区三区| 国产亚洲av片在线观看秒播厂| 老熟女久久久| 久久精品国产鲁丝片午夜精品| 99热6这里只有精品| 少妇熟女欧美另类| 久久久久性生活片| 少妇 在线观看| 最近的中文字幕免费完整| 在线免费观看不下载黄p国产| 丝袜脚勾引网站| a 毛片基地| 麻豆成人午夜福利视频| 能在线免费看毛片的网站| 一级毛片我不卡| 99九九线精品视频在线观看视频| 激情五月婷婷亚洲| 超碰av人人做人人爽久久| 成人特级av手机在线观看| 舔av片在线| 一级二级三级毛片免费看| 亚洲精品久久午夜乱码| 99九九线精品视频在线观看视频| 深爱激情五月婷婷| 超碰av人人做人人爽久久| 精品久久久久久久末码| 亚洲综合精品二区| 中文字幕精品免费在线观看视频 | 丰满人妻一区二区三区视频av| 日日撸夜夜添| 亚州av有码| 国产亚洲91精品色在线| 一级毛片 在线播放| 少妇人妻久久综合中文| 亚洲国产欧美在线一区| 亚洲欧美成人综合另类久久久| 亚洲熟女精品中文字幕| 久久99精品国语久久久| 亚洲美女黄色视频免费看| 日本午夜av视频| 一区二区av电影网| 日日撸夜夜添| 精品午夜福利在线看| 亚洲伊人久久精品综合| 99久久中文字幕三级久久日本| 亚洲一级一片aⅴ在线观看| 久久久精品94久久精品| 久久青草综合色| 18禁在线无遮挡免费观看视频| 男男h啪啪无遮挡| 一区二区三区乱码不卡18| 亚洲综合精品二区| 春色校园在线视频观看| 成人毛片60女人毛片免费| 男人添女人高潮全过程视频| 久久精品国产a三级三级三级| 精品人妻一区二区三区麻豆| 99热这里只有是精品50| 小蜜桃在线观看免费完整版高清| 男人添女人高潮全过程视频| 久久精品熟女亚洲av麻豆精品| 久热这里只有精品99| 直男gayav资源| 大又大粗又爽又黄少妇毛片口| 一个人免费看片子| 久久精品国产亚洲av涩爱| 尾随美女入室| av天堂中文字幕网| 久久久午夜欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 97超视频在线观看视频| 欧美另类一区| 久久99热这里只频精品6学生| 男人爽女人下面视频在线观看| 欧美精品亚洲一区二区| 建设人人有责人人尽责人人享有的 | 日本av手机在线免费观看| 麻豆国产97在线/欧美| 日韩成人av中文字幕在线观看| 在现免费观看毛片| 亚洲av综合色区一区| 丝袜喷水一区| 欧美人与善性xxx| 亚洲成色77777| av播播在线观看一区| 日韩三级伦理在线观看| 一区二区三区四区激情视频| 在现免费观看毛片| 欧美极品一区二区三区四区| 国产免费又黄又爽又色| 人妻一区二区av| 人体艺术视频欧美日本| 永久网站在线| 99热国产这里只有精品6| 国产视频首页在线观看| 精品久久久久久久末码| 人妻系列 视频| 国产成人精品一,二区| 大香蕉久久网| 精品人妻偷拍中文字幕| 日韩不卡一区二区三区视频在线| 丰满人妻一区二区三区视频av| 深爱激情五月婷婷| 久久久久网色| 99热这里只有是精品50| 色综合色国产| 欧美精品一区二区免费开放| 国产男女内射视频| 性色avwww在线观看| 91精品国产九色| 在线观看人妻少妇| 国产在线免费精品| 亚洲精品乱码久久久久久按摩| 亚洲av男天堂| 亚洲色图综合在线观看| 欧美xxxx性猛交bbbb| 久久精品国产a三级三级三级| 在线精品无人区一区二区三 | 国产成人a∨麻豆精品| 麻豆精品久久久久久蜜桃| 联通29元200g的流量卡| 国产精品一区二区在线观看99| 国产午夜精品久久久久久一区二区三区| 18禁裸乳无遮挡动漫免费视频| 日本爱情动作片www.在线观看| 国产成人aa在线观看| 亚洲人成网站高清观看| 五月开心婷婷网| 中文欧美无线码| 自拍欧美九色日韩亚洲蝌蚪91 | 免费av中文字幕在线| 亚洲va在线va天堂va国产| 国产美女午夜福利| 99热这里只有是精品50| 人人妻人人澡人人爽人人夜夜| 亚洲精品自拍成人| 在线观看三级黄色| 久久国产亚洲av麻豆专区| 99热6这里只有精品| 黄色欧美视频在线观看| 国产男女超爽视频在线观看| 国产精品无大码| 有码 亚洲区| 国语对白做爰xxxⅹ性视频网站| 国产91av在线免费观看| 99国产精品免费福利视频| 精品国产一区二区三区久久久樱花 | 国产久久久一区二区三区| 麻豆成人av视频| 国产成人一区二区在线| 成人高潮视频无遮挡免费网站| 制服丝袜香蕉在线| 国产高潮美女av| 我的女老师完整版在线观看| 国产精品久久久久久av不卡| 亚洲av综合色区一区| 91狼人影院| 一区二区三区精品91| 日韩av免费高清视频| 99久久精品国产国产毛片| 成年av动漫网址| 成人免费观看视频高清| 亚洲精品日韩av片在线观看| 久久国产乱子免费精品| 噜噜噜噜噜久久久久久91| 下体分泌物呈黄色| 国产一级毛片在线| 亚洲人成网站在线观看播放| 这个男人来自地球电影免费观看 | 国产男女超爽视频在线观看| 亚洲精品乱码久久久久久按摩| 一边亲一边摸免费视频| 大香蕉97超碰在线| 久热久热在线精品观看| 日本黄色日本黄色录像| 男女边摸边吃奶| 精品久久久久久久久av| 在现免费观看毛片| 少妇精品久久久久久久| 国产熟女欧美一区二区| 街头女战士在线观看网站| 日韩三级伦理在线观看| 午夜老司机福利剧场| 亚洲美女搞黄在线观看| 亚洲色图av天堂| 搡女人真爽免费视频火全软件| 毛片一级片免费看久久久久| 直男gayav资源| 国产av精品麻豆| 久久99热这里只有精品18| 国产高清有码在线观看视频| 三级国产精品片| 一个人看视频在线观看www免费| 欧美日韩精品成人综合77777| 久久影院123| 成人二区视频| 久久鲁丝午夜福利片| 欧美精品一区二区免费开放| 人人妻人人添人人爽欧美一区卜 | 精品人妻一区二区三区麻豆| 久久女婷五月综合色啪小说| 久久久久久久久大av| 亚洲成人手机| 特大巨黑吊av在线直播| 永久免费av网站大全| 免费看光身美女| 国产在线男女| 毛片女人毛片| 天美传媒精品一区二区| 国产成人aa在线观看| 欧美日本视频| 欧美日韩在线观看h| 免费观看性生交大片5| 欧美极品一区二区三区四区| 国产成人freesex在线| 欧美bdsm另类| 成人18禁高潮啪啪吃奶动态图 | 久久久久久久大尺度免费视频| 免费看光身美女| 久久久成人免费电影| 精品久久久精品久久久| 日本欧美视频一区| 亚洲精品日韩在线中文字幕| 超碰av人人做人人爽久久| 如何舔出高潮| 久久女婷五月综合色啪小说| 亚洲精品乱久久久久久| 久久久久精品性色| 久久毛片免费看一区二区三区| 91aial.com中文字幕在线观看| 卡戴珊不雅视频在线播放| 国产日韩欧美亚洲二区| 制服丝袜香蕉在线| 国产精品一区二区在线不卡| 免费看日本二区| 国产精品一区二区三区四区免费观看| 亚洲欧美成人综合另类久久久| 日本一二三区视频观看| 国产大屁股一区二区在线视频| 成人亚洲欧美一区二区av| 日韩大片免费观看网站| 在线免费十八禁| 久久精品国产亚洲av涩爱| 日韩制服骚丝袜av| 最近2019中文字幕mv第一页| 久久av网站| 视频区图区小说| 嫩草影院入口| 青春草视频在线免费观看| 日韩av不卡免费在线播放| 久久99热这里只频精品6学生| 亚洲av.av天堂| 国产免费一区二区三区四区乱码| 黄色欧美视频在线观看| 免费高清在线观看视频在线观看| 亚洲,一卡二卡三卡| av在线老鸭窝| 91午夜精品亚洲一区二区三区| 亚洲国产精品一区三区| 成年av动漫网址| 亚洲欧美一区二区三区黑人 | 久久国内精品自在自线图片| 在线观看免费视频网站a站| 国产一区二区三区av在线| 卡戴珊不雅视频在线播放| 国产高清三级在线| 男女边摸边吃奶| 日日摸夜夜添夜夜爱| 日韩电影二区| 美女高潮的动态| 日韩,欧美,国产一区二区三区| 精品一区在线观看国产| 国产深夜福利视频在线观看| 成人毛片a级毛片在线播放| 免费观看无遮挡的男女| h视频一区二区三区| 又大又黄又爽视频免费| 国产人妻一区二区三区在| 国产成人精品福利久久| 日本猛色少妇xxxxx猛交久久| 免费黄网站久久成人精品| 天天躁夜夜躁狠狠久久av| 久久韩国三级中文字幕| 久热这里只有精品99| 少妇人妻 视频| 51国产日韩欧美| av福利片在线观看| 久久影院123| 国产一区有黄有色的免费视频| 欧美日韩在线观看h| 大香蕉97超碰在线| 国产精品一区二区三区四区免费观看| 一级a做视频免费观看| 欧美高清性xxxxhd video| 91午夜精品亚洲一区二区三区| 亚洲婷婷狠狠爱综合网| 亚洲欧美成人精品一区二区| 我的女老师完整版在线观看| 免费看av在线观看网站| 亚洲第一区二区三区不卡| 高清欧美精品videossex| 日本黄色日本黄色录像| 天美传媒精品一区二区| 久久久久久久久久成人| 校园人妻丝袜中文字幕| 高清视频免费观看一区二区| 国产精品免费大片| 新久久久久国产一级毛片| 日本与韩国留学比较| 国产精品嫩草影院av在线观看| 色综合色国产| 在线观看国产h片| 国产日韩欧美在线精品| 男女下面进入的视频免费午夜| 日本黄色片子视频| 夜夜看夜夜爽夜夜摸| 久久6这里有精品| 免费看不卡的av| 亚洲精品乱久久久久久| 纯流量卡能插随身wifi吗| 男的添女的下面高潮视频| 能在线免费看毛片的网站| 99国产精品免费福利视频| 夫妻午夜视频| 丝袜喷水一区| av播播在线观看一区| 中文字幕精品免费在线观看视频 | 午夜免费观看性视频| 精品久久久久久电影网| 啦啦啦啦在线视频资源| 成人毛片a级毛片在线播放| 国产免费一区二区三区四区乱码| 欧美激情国产日韩精品一区| 国产亚洲午夜精品一区二区久久| 直男gayav资源| 久久精品国产自在天天线| 久久久久性生活片| 日韩欧美 国产精品| 插阴视频在线观看视频| 久久久久国产网址| 成人漫画全彩无遮挡| 久久婷婷青草| 亚洲精品色激情综合| 亚洲精品一区蜜桃| 少妇丰满av| 干丝袜人妻中文字幕| 在线观看av片永久免费下载| 中文字幕精品免费在线观看视频 | 久久毛片免费看一区二区三区| 黑人高潮一二区| 久久久久久久国产电影| 少妇人妻一区二区三区视频| 亚洲精品第二区| 成人毛片a级毛片在线播放| 99久久精品国产国产毛片| 欧美日韩综合久久久久久| av国产精品久久久久影院| 蜜桃在线观看..| 麻豆国产97在线/欧美| 国产亚洲av片在线观看秒播厂| 免费播放大片免费观看视频在线观看| 国产精品国产三级国产专区5o| 91久久精品电影网| 午夜福利影视在线免费观看| 免费大片黄手机在线观看| 天美传媒精品一区二区| 熟女av电影| av一本久久久久| 午夜免费观看性视频| 大码成人一级视频| av在线app专区| 国产 一区精品| 99热6这里只有精品| 欧美精品人与动牲交sv欧美| av不卡在线播放| 成人亚洲精品一区在线观看 | av女优亚洲男人天堂| 国产av精品麻豆| 男男h啪啪无遮挡| www.av在线官网国产| 男女国产视频网站| 久热这里只有精品99| 一级a做视频免费观看| 亚洲图色成人| 看非洲黑人一级黄片| 久久99热6这里只有精品| 久久99精品国语久久久| 亚洲欧美日韩另类电影网站 | 一级毛片 在线播放| 尤物成人国产欧美一区二区三区| 少妇丰满av| 久久久久久久大尺度免费视频| 深夜a级毛片| 日韩 亚洲 欧美在线| 中文天堂在线官网| 免费观看的影片在线观看| 国产av精品麻豆| 18禁裸乳无遮挡免费网站照片| 最近中文字幕2019免费版| 97热精品久久久久久| 老司机影院毛片| 六月丁香七月| 国产熟女欧美一区二区| 久久国内精品自在自线图片| 制服丝袜香蕉在线| 男人和女人高潮做爰伦理| 永久网站在线| 国产老妇伦熟女老妇高清| 中文天堂在线官网| 黄色视频在线播放观看不卡| 久久精品人妻少妇| 99热6这里只有精品| 免费人妻精品一区二区三区视频| 日本黄色片子视频| 观看av在线不卡| 高清毛片免费看| 免费播放大片免费观看视频在线观看| 国产v大片淫在线免费观看| 久久精品国产自在天天线| 狠狠精品人妻久久久久久综合| 99久久精品热视频| 国产精品麻豆人妻色哟哟久久| 国产国拍精品亚洲av在线观看| a级一级毛片免费在线观看| 国产高清国产精品国产三级 | 一边亲一边摸免费视频| 欧美变态另类bdsm刘玥|