• <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久久网站| 亚洲第一av免费看| 亚洲一区中文字幕在线| 人妻久久中文字幕网| 女性生殖器流出的白浆| 好男人电影高清在线观看| 我的亚洲天堂| 夜夜夜夜夜久久久久| 久久免费观看电影| 99国产精品免费福利视频| 999久久久精品免费观看国产| 国产精品一区二区免费欧美| 国产成人av教育| 性色av乱码一区二区三区2| 国产精品久久久人人做人人爽| 国产亚洲av高清不卡| 真人做人爱边吃奶动态| 国产成人av教育| 纵有疾风起免费观看全集完整版| 一区福利在线观看| 美女视频免费永久观看网站| 亚洲色图综合在线观看| 一本色道久久久久久精品综合| 嫁个100分男人电影在线观看| 日韩视频一区二区在线观看| 999久久久国产精品视频| 日韩免费高清中文字幕av| av又黄又爽大尺度在线免费看| 最近最新中文字幕大全免费视频| 久久精品国产亚洲av香蕉五月 | 精品久久久久久电影网| 精品国产乱码久久久久久男人| 精品视频人人做人人爽| 亚洲欧美激情在线| 黑人操中国人逼视频| 熟女少妇亚洲综合色aaa.| 丝袜人妻中文字幕| 亚洲男人天堂网一区| 日韩欧美一区视频在线观看| 久久狼人影院| 国产不卡一卡二| 日韩视频在线欧美| www.精华液| 最新美女视频免费是黄的| 黄色怎么调成土黄色| 午夜日韩欧美国产| 欧美乱码精品一区二区三区| av不卡在线播放| 一级黄色大片毛片| 精品亚洲乱码少妇综合久久| 国产精品久久久av美女十八| 男女床上黄色一级片免费看| 精品久久久久久久毛片微露脸| 国产一区二区激情短视频| 亚洲黑人精品在线| 久久久久久久精品吃奶| 中亚洲国语对白在线视频| 亚洲av欧美aⅴ国产| 亚洲国产毛片av蜜桃av| 国产精品秋霞免费鲁丝片| 国产精品av久久久久免费| 成人国产一区最新在线观看| 亚洲一区中文字幕在线| 国产精品熟女久久久久浪| 欧美日韩av久久| 久久青草综合色| 国产精品一区二区免费欧美| 精品第一国产精品| 精品久久久久久电影网| 久久午夜亚洲精品久久| 亚洲av日韩在线播放| 久久精品人人爽人人爽视色| 久久国产精品男人的天堂亚洲| 美女视频免费永久观看网站| 777久久人妻少妇嫩草av网站| 99国产精品99久久久久| 国产片内射在线| 亚洲国产av影院在线观看| 亚洲国产欧美一区二区综合| 女性生殖器流出的白浆| 日韩视频在线欧美| 国产精品影院久久| 亚洲精品中文字幕在线视频| 中文字幕高清在线视频| 高清av免费在线| 成年人黄色毛片网站| 亚洲黑人精品在线| 成在线人永久免费视频| 精品久久久久久电影网| 大型av网站在线播放| 久久中文字幕一级| 国产av精品麻豆| 在线天堂中文资源库| 免费少妇av软件| 少妇粗大呻吟视频| 中文亚洲av片在线观看爽 | 男女之事视频高清在线观看| 国产一卡二卡三卡精品| 亚洲视频免费观看视频| 午夜福利,免费看| 午夜福利,免费看| 久久av网站| 99国产精品免费福利视频| 99国产精品免费福利视频| 午夜福利在线免费观看网站| 欧美在线黄色| 少妇 在线观看| 免费日韩欧美在线观看| 亚洲一区中文字幕在线| 色播在线永久视频| 亚洲中文日韩欧美视频| 黑人操中国人逼视频| 国产1区2区3区精品| 亚洲国产中文字幕在线视频| 后天国语完整版免费观看| 制服诱惑二区| 波多野结衣一区麻豆| a级毛片在线看网站| 国产精品av久久久久免费| 国产一区二区激情短视频| 99在线人妻在线中文字幕 | 色老头精品视频在线观看| 午夜91福利影院| 啦啦啦在线免费观看视频4| 成人18禁高潮啪啪吃奶动态图| 国产精品 欧美亚洲| 大码成人一级视频| 亚洲伊人久久精品综合| 制服诱惑二区| 亚洲精品自拍成人| 一级片免费观看大全| av超薄肉色丝袜交足视频| 久久人妻熟女aⅴ| 99精国产麻豆久久婷婷| 黄色视频不卡| 天天添夜夜摸| 国产精品av久久久久免费| 欧美午夜高清在线| 人人妻,人人澡人人爽秒播| 亚洲熟女毛片儿| 激情在线观看视频在线高清 | www.熟女人妻精品国产| 大片电影免费在线观看免费| 精品久久蜜臀av无| 久久亚洲真实| 91av网站免费观看| 一夜夜www| 麻豆av在线久日| 国产一区二区在线观看av| 啪啪无遮挡十八禁网站| 亚洲精品粉嫩美女一区| 久久av网站| 99re6热这里在线精品视频| 99精品欧美一区二区三区四区| 精品少妇黑人巨大在线播放| 国产亚洲欧美在线一区二区| 亚洲精品国产精品久久久不卡| 国产伦理片在线播放av一区| 亚洲男人天堂网一区| 欧美老熟妇乱子伦牲交| 国产精品国产av在线观看| 搡老岳熟女国产| 在线观看一区二区三区激情| 少妇精品久久久久久久| 日韩欧美国产一区二区入口| 好男人电影高清在线观看| 啦啦啦在线免费观看视频4| 91老司机精品| 一级毛片精品| 久久影院123| 一级毛片电影观看| 成人三级做爰电影| 欧美日韩亚洲国产一区二区在线观看 | 欧美人与性动交α欧美精品济南到| tocl精华| 首页视频小说图片口味搜索| 最新美女视频免费是黄的| 午夜免费鲁丝| 在线十欧美十亚洲十日本专区| 两个人免费观看高清视频| 十八禁网站免费在线| 在线天堂中文资源库| 日韩有码中文字幕| 午夜福利视频精品| 精品福利观看| 久久精品国产a三级三级三级| 国产伦理片在线播放av一区| 亚洲av欧美aⅴ国产| 国产单亲对白刺激| 亚洲午夜理论影院| 欧美日韩黄片免| 日韩制服丝袜自拍偷拍| 日韩欧美一区二区三区在线观看 | 中文字幕精品免费在线观看视频| www日本在线高清视频| 国产在线一区二区三区精| 一进一出好大好爽视频| 女人久久www免费人成看片| 最近最新中文字幕大全免费视频| 欧美乱妇无乱码| 91大片在线观看| 亚洲精品成人av观看孕妇| 亚洲国产毛片av蜜桃av| 天天操日日干夜夜撸| 老司机午夜十八禁免费视频| 一级毛片女人18水好多| 国产亚洲一区二区精品| 巨乳人妻的诱惑在线观看| 国产亚洲精品第一综合不卡| 精品亚洲成a人片在线观看| 99国产精品免费福利视频| 69av精品久久久久久 | 在线观看66精品国产| 久久精品国产综合久久久| 人人妻人人爽人人添夜夜欢视频| 成年人午夜在线观看视频| 欧美国产精品一级二级三级| 精品欧美一区二区三区在线| 亚洲综合色网址| 亚洲欧洲日产国产| 成人18禁高潮啪啪吃奶动态图| 熟女少妇亚洲综合色aaa.| 国产免费现黄频在线看| 18禁美女被吸乳视频| 电影成人av| 极品人妻少妇av视频| 高清av免费在线| 一级a爱视频在线免费观看| 天天躁夜夜躁狠狠躁躁| 欧美激情 高清一区二区三区| 高清毛片免费观看视频网站 | 成年版毛片免费区| 午夜福利视频在线观看免费| 国产伦理片在线播放av一区| 亚洲av电影在线进入| 久久久国产成人免费| 午夜成年电影在线免费观看| 亚洲自偷自拍图片 自拍| 狠狠狠狠99中文字幕| 欧美日韩福利视频一区二区| 成人18禁在线播放| 色播在线永久视频| 成人特级黄色片久久久久久久 | 两性夫妻黄色片| 午夜日韩欧美国产| 国产又色又爽无遮挡免费看| 交换朋友夫妻互换小说| 在线播放国产精品三级| 麻豆乱淫一区二区| 国产精品.久久久| av福利片在线| www日本在线高清视频| 久久精品熟女亚洲av麻豆精品| 免费人妻精品一区二区三区视频| 久久久久久久久免费视频了| 午夜精品国产一区二区电影| 精品亚洲成a人片在线观看| 国产一区二区激情短视频| 波多野结衣一区麻豆| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 精品人妻熟女毛片av久久网站| tocl精华| 亚洲 欧美一区二区三区| 新久久久久国产一级毛片| 精品福利永久在线观看| 人人澡人人妻人| 狠狠精品人妻久久久久久综合| 免费高清在线观看日韩| 自拍欧美九色日韩亚洲蝌蚪91| 国产视频一区二区在线看| 久久精品人人爽人人爽视色| 午夜激情久久久久久久| 日本黄色视频三级网站网址 | 91麻豆精品激情在线观看国产 | 亚洲精品中文字幕在线视频| a级片在线免费高清观看视频| 久久香蕉激情| 大片免费播放器 马上看| 欧美日韩av久久| 日韩精品免费视频一区二区三区| 91国产中文字幕| 在线观看舔阴道视频| 另类亚洲欧美激情| 国产欧美日韩综合在线一区二区| 欧美大码av| 十八禁高潮呻吟视频| 麻豆av在线久日| 午夜福利视频精品| 久久精品熟女亚洲av麻豆精品| 精品国产乱码久久久久久小说| www日本在线高清视频| 精品亚洲乱码少妇综合久久| 91av网站免费观看| 大片电影免费在线观看免费| 不卡一级毛片| 成人亚洲精品一区在线观看| 国产男女内射视频| 每晚都被弄得嗷嗷叫到高潮| 精品国产一区二区三区久久久樱花| 高清在线国产一区| 久久久久久久大尺度免费视频| 久久精品人人爽人人爽视色| 日本撒尿小便嘘嘘汇集6| 脱女人内裤的视频| 久久久久久久国产电影| 久久精品人人爽人人爽视色| 久久精品亚洲av国产电影网| 丰满迷人的少妇在线观看| 天堂俺去俺来也www色官网| 欧美日本中文国产一区发布| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩亚洲高清精品| cao死你这个sao货| 最新美女视频免费是黄的| 在线观看人妻少妇| 久久精品国产亚洲av高清一级| 国产在线精品亚洲第一网站| 国产精品影院久久| 精品久久久久久久毛片微露脸| 中文字幕人妻熟女乱码| 色老头精品视频在线观看| videosex国产| 国产aⅴ精品一区二区三区波| 97在线人人人人妻| 国产黄频视频在线观看| 夜夜骑夜夜射夜夜干| 亚洲成人国产一区在线观看| 国产极品粉嫩免费观看在线| 又黄又粗又硬又大视频| 国产精品麻豆人妻色哟哟久久| 制服诱惑二区| 激情在线观看视频在线高清 | 欧美精品一区二区大全| 国产1区2区3区精品| 一进一出好大好爽视频| 可以免费在线观看a视频的电影网站| 久久久久久久大尺度免费视频| 久久中文字幕人妻熟女| 大片电影免费在线观看免费| 韩国精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 热re99久久国产66热| 在线观看一区二区三区激情| 人成视频在线观看免费观看| 亚洲av片天天在线观看| 国产91精品成人一区二区三区 | 久久精品熟女亚洲av麻豆精品| 一区二区av电影网| 麻豆成人av在线观看| 久久人妻熟女aⅴ| 人人澡人人妻人| 国产又色又爽无遮挡免费看| 亚洲伊人久久精品综合| 黄色 视频免费看| 母亲3免费完整高清在线观看| 99九九在线精品视频| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 十八禁高潮呻吟视频| 精品国产乱码久久久久久小说| 黄网站色视频无遮挡免费观看| 老汉色av国产亚洲站长工具| 国产区一区二久久| 国产精品国产高清国产av | cao死你这个sao货| 天天躁日日躁夜夜躁夜夜| 欧美黑人欧美精品刺激| 午夜激情久久久久久久| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 九色亚洲精品在线播放| 免费观看人在逋| 法律面前人人平等表现在哪些方面| 动漫黄色视频在线观看| 啦啦啦 在线观看视频| 亚洲精品国产区一区二| 人成视频在线观看免费观看| 无限看片的www在线观看| 十八禁网站网址无遮挡| 伦理电影免费视频| 午夜福利乱码中文字幕| 精品欧美一区二区三区在线| 免费看a级黄色片| 国产精品秋霞免费鲁丝片| 午夜视频精品福利| 欧美乱妇无乱码| 亚洲专区国产一区二区| 亚洲九九香蕉| 亚洲午夜理论影院| 十八禁网站网址无遮挡| 国产一区二区 视频在线| 乱人伦中国视频| 麻豆av在线久日| 国产亚洲精品一区二区www | 99re6热这里在线精品视频| 色尼玛亚洲综合影院| 国产亚洲精品第一综合不卡| 日本一区二区免费在线视频| 老司机靠b影院| 黄片大片在线免费观看| 午夜激情av网站| 女人久久www免费人成看片| 精品欧美一区二区三区在线| 日韩一区二区三区影片| 黄色视频不卡| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 久热爱精品视频在线9| 亚洲国产欧美一区二区综合| 99国产精品一区二区三区| 国产男女超爽视频在线观看| 91精品国产国语对白视频| 男女床上黄色一级片免费看| 欧美精品av麻豆av| av视频免费观看在线观看| 欧美另类亚洲清纯唯美| 精品国产一区二区三区四区第35| 国产欧美日韩一区二区精品| www日本在线高清视频| 91成年电影在线观看| 老司机午夜福利在线观看视频 | cao死你这个sao货| 老司机亚洲免费影院| 国产精品av久久久久免费| 男男h啪啪无遮挡| 国产欧美亚洲国产| 色视频在线一区二区三区| 91麻豆精品激情在线观看国产 | 两性夫妻黄色片| 精品国产超薄肉色丝袜足j| 久久精品91无色码中文字幕| 亚洲欧美精品综合一区二区三区| 大型av网站在线播放| 在线十欧美十亚洲十日本专区| 精品免费久久久久久久清纯 | 久久99热这里只频精品6学生| 欧美老熟妇乱子伦牲交| 色尼玛亚洲综合影院| 午夜福利一区二区在线看| 国产成人啪精品午夜网站| 99re在线观看精品视频| 欧美一级毛片孕妇| 欧美av亚洲av综合av国产av| 99久久精品国产亚洲精品| 欧美中文综合在线视频| 嫩草影视91久久| 菩萨蛮人人尽说江南好唐韦庄| 狠狠精品人妻久久久久久综合| 女警被强在线播放| 国产一区二区三区视频了| 久久精品成人免费网站| 欧美日韩av久久| 男女边摸边吃奶| 国产亚洲av高清不卡| 中国美女看黄片| 国产精品亚洲一级av第二区| 久久人妻熟女aⅴ| 亚洲精品av麻豆狂野| 捣出白浆h1v1| 国产1区2区3区精品| 亚洲av电影在线进入| 搡老熟女国产l中国老女人| 日韩欧美免费精品| 亚洲中文日韩欧美视频| 美女高潮到喷水免费观看| 啦啦啦在线免费观看视频4| 五月开心婷婷网| 99精品在免费线老司机午夜| 欧美日本中文国产一区发布| 久久国产精品人妻蜜桃| kizo精华| 久久人人爽av亚洲精品天堂| 日本精品一区二区三区蜜桃| 免费在线观看日本一区| 捣出白浆h1v1| 精品高清国产在线一区| 国产av精品麻豆| 一本大道久久a久久精品| 肉色欧美久久久久久久蜜桃| 热99国产精品久久久久久7| 国产成人精品在线电影| 欧美大码av| 国产欧美日韩精品亚洲av| 黄色怎么调成土黄色| 国产精品九九99| 亚洲av片天天在线观看| 国产一区二区三区在线臀色熟女 | 18禁国产床啪视频网站| 69精品国产乱码久久久| 欧美日韩av久久| 露出奶头的视频| 9色porny在线观看| 在线天堂中文资源库| 精品卡一卡二卡四卡免费| 日韩视频一区二区在线观看| 久久 成人 亚洲| 美女高潮喷水抽搐中文字幕| 国产aⅴ精品一区二区三区波| 男女床上黄色一级片免费看| 中文字幕av电影在线播放| 欧美日本中文国产一区发布| 一级毛片精品| 黄色视频不卡| 日日摸夜夜添夜夜添小说| 亚洲av日韩精品久久久久久密| 女人爽到高潮嗷嗷叫在线视频| 亚洲七黄色美女视频| 国产精品熟女久久久久浪| 狠狠婷婷综合久久久久久88av| 曰老女人黄片| 欧美日韩亚洲综合一区二区三区_| 一区二区三区乱码不卡18| 天天添夜夜摸| av电影中文网址| 高清黄色对白视频在线免费看| 久热爱精品视频在线9| 嫁个100分男人电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久中文字幕一级| 一本久久精品| 精品久久久久久久毛片微露脸| 天天躁日日躁夜夜躁夜夜| 成人18禁高潮啪啪吃奶动态图| 国产日韩欧美亚洲二区| 18禁黄网站禁片午夜丰满| 丝袜在线中文字幕| 在线观看免费午夜福利视频| 国产成人欧美| 久久性视频一级片| 日韩中文字幕欧美一区二区| 成年人黄色毛片网站| 丝袜在线中文字幕| 99热网站在线观看| 亚洲第一av免费看| 国产高清国产精品国产三级| 天堂动漫精品| 精品国内亚洲2022精品成人 | 国产精品免费一区二区三区在线 | 成年女人毛片免费观看观看9 | 黄色视频在线播放观看不卡| 国产免费福利视频在线观看| 岛国在线观看网站| 99国产精品一区二区蜜桃av | 无限看片的www在线观看| 香蕉丝袜av| 一本大道久久a久久精品| 啪啪无遮挡十八禁网站| 亚洲中文字幕日韩| 精品国产一区二区久久| 无遮挡黄片免费观看| 国产高清videossex| 大香蕉久久网| 亚洲va日本ⅴa欧美va伊人久久| 电影成人av| 精品少妇一区二区三区视频日本电影| 丝袜人妻中文字幕| 欧美成人午夜精品| 欧美亚洲 丝袜 人妻 在线| 欧美日韩亚洲高清精品| 操美女的视频在线观看| 国产精品久久久av美女十八| 国产成人免费观看mmmm| 国产日韩欧美在线精品| 久久午夜综合久久蜜桃| 成人亚洲精品一区在线观看| 色老头精品视频在线观看| 两性夫妻黄色片| 男女午夜视频在线观看| 国产男女内射视频| 日韩欧美一区视频在线观看| 亚洲精品国产色婷婷电影| 一本一本久久a久久精品综合妖精| 精品少妇久久久久久888优播| av有码第一页| 91九色精品人成在线观看| 国产亚洲精品久久久久5区| 国产精品1区2区在线观看. | 久久性视频一级片| 免费看a级黄色片| 国产高清videossex| videosex国产| 黄色视频在线播放观看不卡| 亚洲天堂av无毛| 国产亚洲精品第一综合不卡| 在线播放国产精品三级| 成人免费观看视频高清| 国产精品久久久久久精品电影小说| 亚洲熟女精品中文字幕| 男女午夜视频在线观看| 国产人伦9x9x在线观看| 亚洲精品国产精品久久久不卡| 大片免费播放器 马上看| 久久人妻熟女aⅴ| 欧美成人午夜精品| 日韩一卡2卡3卡4卡2021年| 久久精品熟女亚洲av麻豆精品| 一级毛片女人18水好多| 肉色欧美久久久久久久蜜桃| 一区二区三区激情视频| 精品久久久久久电影网| 他把我摸到了高潮在线观看 | 91av网站免费观看| 成人特级黄色片久久久久久久 | 在线观看免费视频日本深夜| 精品熟女少妇八av免费久了| 人人妻,人人澡人人爽秒播| 国产成人啪精品午夜网站| 国产在线免费精品| 国产精品久久久久久精品电影小说| 香蕉国产在线看| 极品少妇高潮喷水抽搐| 国产精品欧美亚洲77777|