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

    纖維增強復(fù)合材料層壓板的疲勞壽命預(yù)測方法

    2024-03-28 14:40:08徐蓉霞高建雄朱鵬年吳志峯
    上海交通大學(xué)學(xué)報 2024年3期
    關(guān)鍵詞:層壓板鋪層單向

    徐蓉霞, 高建雄, 朱鵬年, 吳志峯

    (新疆大學(xué) 機械工程學(xué)院,烏魯木齊 830017)

    纖維增強復(fù)合材料(Fiber Reinforced Polymer, FRP)具有比強度高、比剛度大、壽命長、結(jié)構(gòu)可設(shè)計[1-3]等諸多優(yōu)點,因此其在近幾十年來應(yīng)用較為廣泛[4]. 在工程實際應(yīng)用中,導(dǎo)致其結(jié)構(gòu)破壞的重要原因包括疲勞損傷造成的斷裂,但FRP有多種形式的損傷,且疲勞破壞機理十分復(fù)雜[5-6]. 在疲勞過程中,各種形式的損傷互相影響,不斷累積,FRP的機械性能也隨之發(fā)生退化[7]. FRP的疲勞壽命模型一般是根據(jù)應(yīng)力-壽命曲線(S-N曲線)或Goodman疲勞壽命曲線來預(yù)測特定載荷下的疲勞壽命. 這種方法不考慮疲勞損傷的微觀機理、材料的性能退化、損傷之間的耦合和損傷的累積,而且往往需要大量的實驗數(shù)據(jù)來建立疲勞壽命曲線. 或者通過漸進疲勞損傷分析模型來預(yù)測其疲勞壽命,該方法涵蓋疲勞加載過程中應(yīng)力分析、疲勞失效準(zhǔn)則和材料性能退化模型3部分. 因此,FRP層壓板的疲勞壽命,可以基于經(jīng)典層壓板理論分析單層板的應(yīng)力來預(yù)測,但該計算過程繁瑣,且認(rèn)為同一單層板的受力均勻,因此預(yù)測誤差較大. 基于以上問題,本文提出了單向板的載荷水平-疲勞壽命曲線(q-N曲線),與層壓板的有限元分析緊密結(jié)合來預(yù)測FRP層壓板的疲勞壽命方法.

    目前用于表征復(fù)合材料疲勞損傷擴展的參量有疲勞模量、剩余強度、應(yīng)變等物理量,基于這些物理量,可以構(gòu)造出多種疲勞損傷和疲勞壽命模型[8-9]. 本文用疲勞模量來描述碳纖維增強復(fù)合材料的損傷演化規(guī)律,從而建立疲勞壽命的函數(shù)表達式. 以單向板在單軸循環(huán)應(yīng)力加載下載荷比與疲勞壽命的試驗數(shù)據(jù)為基礎(chǔ),利用灰色系統(tǒng)等間距GM(1, 1)模型[10],構(gòu)造單向板的疲勞壽命和損傷函數(shù)模型.

    復(fù)合材料的疲勞過程從本質(zhì)上講是損傷不斷累積,材料性質(zhì)不斷退化,應(yīng)力重新分布的動態(tài)過程[11-12]. 本文將FRP層壓板的疲勞失效過程簡化為:在應(yīng)力加載下,單層板的危險單元失效,導(dǎo)致該單層的材料性質(zhì)發(fā)生改變,從而對該層材料進行剛度退化,隨之應(yīng)力分布改變的循環(huán)過程. 利用Abaqus對層壓板建模分析,獲取每單層的應(yīng)力分布,進而提取危險點的應(yīng)力,通過Hashin損傷失效準(zhǔn)則[13],判斷該層的失效模式,然后對失效單層進行剛度退化,循環(huán)該過程,確定單層的失效強度以及層壓板的拉伸強度. 通過該有限元分析過程中所得的應(yīng)力分布,與灰色系統(tǒng)建立的單向板的疲勞壽命計算公式結(jié)合,以此來預(yù)測FRP層壓板的疲勞壽命. 通過與實驗數(shù)據(jù)以及其他文獻方法的預(yù)測結(jié)果對比,來驗證該方法的有效性.

    1 基于灰色系統(tǒng)GM(1, 1)模型建立 單向板的q-N函數(shù)

    灰色系統(tǒng)理論是一種研究“小樣本”“貧信息”不確定性問題的方法,主要提取“部分”已知信息,實現(xiàn)對系統(tǒng)演化規(guī)律的正確描述[14]. 灰色預(yù)測建模方法種類繁多,其中等間距GM(1, 1)模型具有樣本需求小、建模過程簡單等優(yōu)點. 當(dāng)建立單向板的疲勞壽命預(yù)測模型時,在實驗數(shù)據(jù)較少的情況下,若用最小二乘法擬合函數(shù)曲線,將導(dǎo)致誤差較大. 因此,借助灰色系統(tǒng)建立單向板的疲勞壽命函數(shù),首先引入疲勞模量衰減的雙參數(shù)模型.

    1.1 疲勞模量的概念

    Hwang等[15]首先引入疲勞模量的概念,提出通過疲勞模量來預(yù)測疲勞壽命. 文獻[16]中以碳纖維為研究對象,給出了循環(huán)應(yīng)力最大值與應(yīng)力循環(huán)次數(shù)n有關(guān)的疲勞模量的定義:

    (1)

    式中:σmax為外加循環(huán)應(yīng)力的最大值;En為第n個循環(huán)時的疲勞模量;載荷比q表征應(yīng)力水平,q=σmax/σu,σu為極限強度;εn為第n個循環(huán)后的累積應(yīng)變,包含了第n次循環(huán)時的彈性應(yīng)變.

    En的初值、終值條件為

    (2)

    式中:E0為材料的靜態(tài)彈性模量,在工程上,假定n=0時的疲勞模量等于E0;EN為經(jīng)過N個循環(huán)后,材料發(fā)生疲勞破壞時的疲勞模量,N即為疲勞壽命.因此,E0是常數(shù),EN取決于外載荷. Hwang等[17]提出了疲勞模量衰減三參數(shù)模型以及可以簡化的雙參數(shù)模型. 其中,如下式所示的雙參數(shù)模型的疲勞壽命預(yù)測值與實驗數(shù)據(jù)吻合較好,而且簡單,因此本文選用此模型預(yù)測疲勞壽命:

    (3)

    式中:A、B均為材料常數(shù).

    1.2 FRP單向板的疲勞壽命分析

    以單向板受縱向拉-拉循環(huán)載荷為例,來說明單軸循環(huán)應(yīng)力下單向板疲勞壽命的計算方法. 定義如下物理量[18]:

    (4)

    式中:Xt為單向板的縱向拉伸強度;εut為縱向靜態(tài)極限應(yīng)變;E11為縱向靜態(tài)彈性模量;εN表示疲勞破壞時的極限應(yīng)變.將式(3)從第n1個循環(huán)到第n2個循環(huán)進行積分:

    (5)

    令nn=n,n1=1,式(5)可以簡化為

    (6)

    (7)

    假設(shè)單向板在承受N次循環(huán)載荷后發(fā)生疲勞破壞,那么EN便可作為疲勞模量的臨界值.令k1=εN/εut,則N可以表示為

    (8)

    (9)

    式中:a1、k1均為材料參數(shù).當(dāng)B=1時,式(9)簡化為N=exp[a1(k1-q)].

    1.3 損傷函數(shù)的構(gòu)成

    假設(shè)復(fù)合材料無初始損傷,疲勞損傷函數(shù)的定義Dn=(εn-ε0)/(εN-ε0)應(yīng)滿足下列初值、終值條件:

    (10)

    以單向板受縱向單軸循環(huán)應(yīng)力作用為例,以應(yīng)變和疲勞模量為中間量,建立損傷函數(shù),其過程如下:

    (11)

    同樣的方法,可以建立單向板在其他兩個主軸方向上(橫向和面內(nèi)剪切)形如式(11)所示的疲勞損傷函數(shù). 因此三參數(shù)和雙參數(shù)疲勞損傷累積公式分別如下所示:

    (12)

    (13)

    根據(jù)式(12)和式(13)所示的損傷累積公式,計算單向板在3個主軸方向上循環(huán)應(yīng)力加載造成的損傷,分別表示為:D1、D2、D12.考慮多向應(yīng)力的綜合作用效果,定義復(fù)合型損傷:

    (14)

    式中:D11表示多軸循環(huán)應(yīng)力加載下纖維的損傷程度;D22則反映基體的損傷程度. 只要其中任何一個首先達到1,單向板就發(fā)生了疲勞失效. 但疲勞極限本身存在模糊性,因此當(dāng)復(fù)合型損傷很接近1時單向板也有可能發(fā)生了疲勞失效.

    (15)

    式中:Yt為橫向拉伸強度;Sin為面內(nèi)剪切強度;σ1,max、σ2,max表示縱向和橫向兩個方向上的最大主應(yīng)力;τ12,max為剪切方向最大應(yīng)力;νij表示應(yīng)力作用在i方向時,在j方向引起的橫向變形的泊松比,i,j=1, 2.將等效載荷比代入建立的q-N函數(shù)表達式中,計算對應(yīng)的疲勞壽命.

    1.4 灰色系統(tǒng)等間距GM(1, 1)建模

    灰色系統(tǒng)模型通過序列的累加生成,揭示系統(tǒng)的發(fā)展趨勢,累加生成處理后的序列呈現(xiàn)單調(diào)遞增的規(guī)律. 設(shè)X(0)為原始序列,X(1)為一次累加生成序列,Z(1)為均值序列,表達式為X(0)=[x(0)(1)x(0)(2) …x(0)(n)],X(1)=[x(1)(1)x(1)(2) …x(1)(n)],Z(1)=[z(1)(1)z(1)(2) …z(1)(n)].

    其中,累加序列當(dāng)序列數(shù)k=1時,x(1)(1)=x(0)(0),當(dāng)k為其他值時,如下式所示:

    x(1)(k)=(x(0)(k)+x(1)(k-1))

    (16)

    k=2,3,…,n

    均值序列為

    z(1)(k)=0.5(x(1)(k)+x(1)(k-1)),

    k=2,3,…,n

    (17)

    則x(1)上的灰微分方程為x(0)(k)+az(1)(k)=b,其白化方程為dx(1)/dt+ax(1)=b.其中,參數(shù)a和b根據(jù)下式計算:

    (18)

    (19)

    參數(shù)C、D、E、F可以根據(jù)原始序列、一次累加序列以及均值序列得到:

    (20)

    (21)

    (22)

    (23)

    白化方程式的解為x(1)(t)=(x(1)(1)-b/a)exp(-at)+b/a,其中t表示時間,則灰色系統(tǒng)GM(1, 1)模型的時間響應(yīng)序列為x(1)(k+1)=(x(0)(1)-b/a)exp(-ak)+b/a,其中k=1,2,…,n.因此,該模型的預(yù)測結(jié)果計算式如下:

    x(0)(k+1)=x(1)(k+1)-x(1)(k)=

    (24)

    k=1,2,…,n

    設(shè)單層板在單軸循環(huán)應(yīng)力作用下,不同載荷比下的疲勞壽命視為一個灰色系統(tǒng)中,疲勞壽命構(gòu)成一個原始序列,根據(jù)灰色等間距GM(1, 1)建模,擴展實驗數(shù)據(jù),然后擬合q-N函數(shù)曲線.

    1.5 基于灰色系統(tǒng)建立FRP單層板在單軸循環(huán)應(yīng) 力作用下的q-N曲線和損傷函數(shù)

    對于T300/QY8911復(fù)合材料單向板[0]16, 文獻[19]中給出了應(yīng)力比R=0.1時的縱向拉-拉疲勞試驗數(shù)據(jù),如表1所示.

    表1 [0]16單向板在縱向應(yīng)力作用下疲勞試驗數(shù)據(jù)

    首先,用最小二乘法擬合雙參數(shù)和三參數(shù)q-N曲線的函數(shù)表達式如下所示:

    q=-0.035lnN+1.189

    (25)

    q=(1.778×10-5lnN+2.132)-6.597

    (26)

    建立灰色等間距模型時,需要等間距載荷比下的疲勞壽命,根據(jù)插值法,計算q=0.75時的N=2.71×105. 因此灰色等間距GM(1, 1)模型原始序列可以設(shè)為X(0)=[2.50×1032.40×1041.43×1052.71×1054.57×105],根據(jù)式(16)~(23),計算一次累加序列和均值序列,C=1.086 5×106,D=8.95×105,E=4.027 5×1011,F=5.504×1011,a=-0.625 4,b=53 882.根據(jù)式(24)計算當(dāng)q=0.65,0.6,0.55,0.5時,N=9.40×105,1.76×106,3.28×106,6.14×106.在灰色系統(tǒng)建模時,發(fā)現(xiàn)q=0.85時曲線突變,因此將該曲線分段,擬合的雙參數(shù)和三參數(shù)q-N函數(shù),分別如下所示:

    (27)

    (28)

    將4種方法進行誤差對比分析,結(jié)果如表2所示.

    表2 擬合[0]16單向板q-N函數(shù)不同方法的精度對比Tab.2 Precision comparison of different methods for fitting q-N function of [0]16 uni-directional plates

    由表2可知,最小二乘法擬合三參數(shù)q-N函數(shù)的方法誤差太大, 因此不再做函數(shù)曲線圖.其他3種方法的q-lgN函數(shù)曲線如圖1所示.

    圖1 3種建模方式擬合[0]16單向板的q-lg N曲線Fig.1 Fitting of q-lg N curves of [0]16 uni-directional plates of three modeling methods

    通過表2和圖1可以發(fā)現(xiàn),通過灰色系統(tǒng)建立的q-N函數(shù),精度明顯得到了提高. 本文在計算縱向損傷量時,選用擬合曲線與實驗數(shù)據(jù)相對誤差最小的灰色系統(tǒng)擬合的三參數(shù)q-N函數(shù). 因此損傷函數(shù)為

    (29)

    橫向拉-拉建模數(shù)據(jù)根據(jù)[90]20單向板拉-拉疲勞實驗獲得[19],如表3所示.

    表3 [90]20單向板在橫向應(yīng)力作用下的疲勞實驗數(shù)據(jù)(R=0.1)

    首先,用最小二乘法擬合q-N函數(shù). 擬合雙參數(shù)q-N曲線的函數(shù)表達式如下所示:

    q=-0.036 4lnN+1.128

    (30)

    擬合三參數(shù)q-N函數(shù)的發(fā)展趨勢不符合復(fù)合材料q-N函數(shù)發(fā)展趨勢,因此暫不考慮該方法. 根據(jù)插值法,計算出q=0.95,0.85,0.75時,N=550.5, 2.245×103, 8.124×104. 灰色等間距模型原始序列設(shè)為X(0)=[550.5 1.01×1032.245×1033.48×1038.124×1041.59×105]. 根據(jù)式(16)~(23)計算一次累加序列和均值序列,C=225 215,D=246 975,E=3.063 4×1010,F=3.056 7×1010,a=-0.955 3.根據(jù)式(24),計算當(dāng)q=0.65,0.6,0.55時,N=1.369 1×106,3.559×106,9.251 4×106.將這9組數(shù)據(jù),用最小二乘法擬合雙參數(shù)和三參數(shù)q-N函數(shù),結(jié)果如下所示:

    q=-0.034 1lnN+1.111

    (31)

    q=(0.27lnN-0.519)-0.391

    (32)

    將基于橫向拉-拉疲勞試驗數(shù)據(jù)的3種擬合方法的誤差進行對比分析,結(jié)果如表4所示.

    表4 擬合[90]20單向板q-N函數(shù)不同方法的精度對比Tab.4 Precision comparison of different methods for fitting q-N function of [90]20 uni-directional plates

    3種方法擬合的q-lgN函數(shù)曲線,如圖2所示.

    灰色系統(tǒng)擬合的q-N三參數(shù)方法的誤差最小,根據(jù)損傷函數(shù),計算3個參數(shù)時發(fā)現(xiàn),參數(shù)k不是實數(shù),不能計算本文運用的損傷函數(shù). 因此,選用灰色系統(tǒng)擬合的雙參數(shù)q-N函數(shù). 所以損傷函數(shù)為

    (33)

    面內(nèi)剪切疲勞實驗數(shù)據(jù),根據(jù)[0/90]4s剪-剪疲勞實驗獲得[19],如表5所示.

    根據(jù)這5組數(shù)據(jù),用最小二乘法擬合雙參數(shù)和三參數(shù)q-N函數(shù),結(jié)果如下所示:

    圖2 3種建模方式擬合[90]20單向板的q-lg N曲線Fig.2 Fitting of [90]20 uni-directional plates q-lg N of three modeling methods

    表5 [0/90]4s層壓板在面內(nèi)剪切應(yīng)力作用下的疲勞實驗數(shù)據(jù)(R=0.1)

    q=-0.051 1lnN+1.117

    (34)

    q=(-0.042 0lnN+1.108)1.354

    (35)

    根據(jù)插值法,計算q=0.75,0.65,0.55時,N=1 874, 14 075, 101 400. 設(shè)原始序列為X(0)=[2293981 8743.35×1031.407 5×1042.48×1041.014×1051.59×105]. 根據(jù)式(16)~(23),計算一次累加序列和均值序列,C=381 930,D=323 897,E=5.252 9×1010,F=6.562 2×1010,a=-0.778 3.根據(jù)式(24)計算q=0.45,0.4,0.35時,N=1.399 7×106,3.048 4×106,6.638 8×106.基于這9組數(shù)據(jù),用最小二乘法擬合雙參數(shù)和三參數(shù)q-N函數(shù),結(jié)果如下所示:

    q=-0.046 2lnN+1.081

    (36)

    q=(-0.030 8lnN+1.082)1.956 9

    (37)

    將4種方法進行誤差對比分析,結(jié)果如表6所示.

    4種方法擬合q-lgN函數(shù)曲線,擬合結(jié)果如圖3所示.

    通過表6和圖3可以發(fā)現(xiàn),雙參數(shù)q-N曲線的誤差較小. 因此,本文在計算剪切損傷量時,選用雙參數(shù)q-N曲線. 所以損傷函數(shù)為

    (38)

    圖3 4種建模方式擬合[0/90]4s層壓板q-lg N曲線Fig.3 Fitting of q-lg N curves of the [0/90]4s laminate of four modeling methods

    表6 擬合[0/90]4 s層壓板q-N函數(shù)不同方法的精度對比Tab.6 Precision comparison of different methods for fitting [0/90]4 s laminate q-N function

    2 FRP多向?qū)訅喊宓钠趬勖?/h2>

    本文FRP層壓板疲勞壽命預(yù)測的步驟為兩大部分,第1部分是通過應(yīng)力分析和失效分析,計算該層壓板的拉伸強度. 第2部分是基于FRP層壓板的拉伸強度和典型FRP單向板的疲勞壽命預(yù)測函數(shù)以及損傷函數(shù)來計算層壓板的疲勞壽命. 對FRP層壓板的疲勞壽命預(yù)測過程作如下假設(shè):

    (1) 多向?qū)訅喊逯械母鲉蜗虬寰哂邢嗤钠谛阅?且層壓板沒有初始損傷[18].

    (2) 把層壓板的疲勞失效看成是不同單層逐次發(fā)生失效的間斷過程,把材料性質(zhì)的連續(xù)退化過程簡化為單層的危險部位失效后,材料性質(zhì)退化,應(yīng)力分布改變.

    2.1 FRP多向?qū)訅喊宓膽?yīng)力分析和失效分析

    假設(shè)FRP層壓板處于平面應(yīng)力狀態(tài),利用Abaqus建模對層壓板進行應(yīng)力分析. 具體步驟如下.

    (1) 通過Abaqus仿真拉伸實驗,在加載外力P時,設(shè)定循環(huán)周期增量P,僅對Δn、Δ2n、Δ3n循環(huán)周期下的層壓板進行應(yīng)力分析和單元失效判斷. 在Δn內(nèi)只進行材料性能漸降,忽略性能突降. 選擇合適的面內(nèi)強度準(zhǔn)則,判斷是否失效. 本文選用Tsai-Hill強度準(zhǔn)則[20]判定:

    (39)

    式中:σ11、σ22、τ12分別表示縱向、橫向及剪切方向的應(yīng)力.

    (2) 計算單向板的失效強度.首先計算各單層的強度比RL,

    (40)

    找出最小強度比RLmin,強度比最小的單層為最先失效層. 然后計算該層的失效強度FS,

    FS=PRLmin

    (41)

    根據(jù)Hashin失效準(zhǔn)則,判斷失效模式:

    (42)

    (43)

    式中:S12為xy面內(nèi)切應(yīng)力.

    若滿足式(42)則為纖維拉伸失效,若滿足式(43)則為基體拉伸失效. 若根據(jù)式(40)判斷出該層已失效,判斷失效模式時,兩個都未達到失效,則根據(jù)相對大小判斷失效模式,造成這種現(xiàn)象的原因是本文沒考慮層間應(yīng)力的影響.

    (3) 對失效的單層進行剛度退化.本文采用Tsai提出的0.4Em剛度退化準(zhǔn)則[21]:纖維斷裂后,該層便不起作用了,其正軸模量矩陣Q退化為零矩陣,即Q11=Q22=Q12=Q66=0;基體的破壞退化規(guī)則為:縱向拉伸模量E11,泊松比ν12保持不變,橫向拉伸模量E22下降至0.56E22,面內(nèi)切變模量G12下降至0.

    (4) 重復(fù)該過程,直到層壓板失效,以單向板的失效強度最大值作為層壓板的拉伸強度.

    應(yīng)力分析和失效分析的具體過程如圖4所示.

    2.2 FRP多向?qū)訅喊宓钠趬勖?/h3>

    (1) 根據(jù)第1部分的結(jié)果得知層壓板的極限拉伸強度SE,確定循環(huán)載荷的最大外載荷Pb:

    Pb=SEq

    (44)

    (2) 復(fù)合材料層壓板在載荷加載時,失效層的破壞可能是靜強度破壞或者疲勞破壞. 如果Pb>FS,屬于靜力破壞,發(fā)生靜力破壞的單層也消耗了一定數(shù)量的疲勞壽命,計算等效的疲勞累積損傷. 如果發(fā)生纖維斷裂,則將該層的累積損傷量D1、D2、D12均置為1,如果發(fā)生基體破壞, 將面內(nèi)橫向應(yīng)力引起的累積損傷量D2置為1[18]. 若Pb

    (4) 分別根據(jù)q-N函數(shù)表達式(28)、(31)、(36)計算其壽命n1、n2、n12,選擇三者中最小的記為nmin,再根據(jù)損傷函數(shù)式(29)、(33)、(38)以nmin為初值進行迭代.

    (5) 根據(jù)式(14)計算復(fù)合損傷D11、D22,判斷D11、D22其中之一是否滿足大于0.999,滿足時迭代結(jié)束,若不滿足,進行循環(huán)迭代,具體迭代方法如圖5 所示. 迭代終值N為該單向板的疲勞壽命. 隨后對該單向板根據(jù)其失效模式,進行剛度退化,即對材料屬性進行改變. 因為材料的疲勞極限存在模糊性,累積損傷也不一定為1,所以當(dāng)循環(huán)迭代時,n增加的緩慢時迭代結(jié)束,根據(jù)本文算例結(jié)果,將本文迭代終值復(fù)合損傷定義為0.999.

    圖5 層壓板疲勞壽命計算Fig.5 Calculation of fatigue life of laminate

    (6) 重復(fù)上述步驟,直到層壓板失效.各單向板的累積壽命和為層壓板的疲勞壽命.

    3 實例驗證與對比分析

    材料T300/QY8911的彈性常數(shù)[22]縱向拉伸模量E11=135 GPa,橫向拉伸模量E22=8.8 GPa,法向拉伸模量E33=8.8 GPa;泊松比ν12=0.33,ν23=0.48,ν13=0.33;Gij表示i-j面內(nèi)切變模量,G12=4.47 GPa,G23=3.2 GPa,G13=4.47 GPa;強度參數(shù)[23]Xt=1 627.5 MPa,縱向壓縮強度Xc=1 226 MPa,Yt=68.4 MPa,橫向壓縮強度Yc=218 MPa,Sin=89.9 MPa.

    3.1 算例1

    層壓板的鋪層方式為[0/90]4s,長40 mm,寬20 mm,單層厚度為0.125 mm[23].根據(jù)材料的彈性參數(shù)和強度參數(shù),采用Abaqus建立分析計算模型,R=0.1,q=0.85.為了提高計算精度,網(wǎng)格沿寬度方向從邊界向中間由密變疏,對網(wǎng)格劃分進行收斂性分析以及網(wǎng)格無關(guān)性檢查,最終確定有限元模型的網(wǎng)格大小是:兩端是邊長為0.5 mm的正方體,中間是邊長為1 mm的正方體. 單元所用類型為Continum shell: SC8R,單元數(shù)量為 9 600,所建有限元模型如圖6所示,鋪層方式如圖7所示.圖中:Angle=0表示該單向板鋪層角為0°,即單向板鋪層的縱向與層壓板參考坐標(biāo)X軸之間的夾角為0°,Angle=90表示該單向板鋪層角為90°;Ply-1—Ply-16分別表示該層壓板的鋪層,共有16層.

    圖6 [0/90]4s層壓板的模型Fig.6 Model of [0/90]4s laminate

    圖7 [0/90]4s的鋪層方式Fig.7 Laying method of [0/90]4s laminate

    在層壓板左側(cè)端面施加完全約束,在層壓板右側(cè)端面施加均布力載荷. 然后施加載荷,當(dāng)P逐漸增加到560 MPa時:90°鋪層危險點的應(yīng)力分布為σ11=-6.271 MPa,σ22=68.34 MPa,τ12=4.644 MPa. 根據(jù)式(39)判斷該層是否失效:

    1.001 1>1

    (45)

    根據(jù)式(45)結(jié)果大于1,所以90° 鋪層失效,該層的失效強度比RL(90°)=0.999 5,失效強度FS(90°)=559.72 MPa. 根據(jù)失效準(zhǔn)則判斷該層的失效模式為基體失效. 0° 鋪層危險點的應(yīng)力分布為σ11=1 096 MPa,σ22=19.93 MPa,τ12=13.93 MPa. 判斷該層是否失效:

    0.554 2<1

    (46)

    根據(jù)式(46)結(jié)果小于1,因此0°鋪層未失效. 所以對90°鋪層根據(jù)本文選定的剛度退化準(zhǔn)則進行剛度退化. 然后施加載荷,當(dāng)載荷增加到738 MPa時:0°鋪層危險點的應(yīng)力分布為σ11=1 464 MPa,σ22=29.8 MPa,τ12=12.27 MPa. 判斷該層是否失效:

    1.001 1>1

    (47)

    由式(47)計算結(jié)果大于1,因此0°鋪層失效,該層的失效強度比RL(0°)=0.999 4,失效強度FS(0°)=737 MPa,失效模式為纖維斷裂. 極限強度SE=737 MPa,根據(jù)式(44),計算循環(huán)載荷的最大外載荷,Pb=SEq=626.5 MPa. 因為FS(90°)Pb,因此,0°鋪層因為疲勞破壞而失效.

    (48)

    D11>1,則該單向板的壽命為2.330 6×105. 因此,該層壓板的疲勞壽命為2.331 06×105. 根據(jù)文獻[23],該層壓板的疲勞壽命實驗數(shù)據(jù)為 2.33×105,將本文預(yù)測結(jié)果、多個文獻方法與實驗數(shù)據(jù)計算相對誤差,結(jié)果進行對比,如表7所示.

    表7 不同方法預(yù)測[0/90]4s層壓板的疲勞壽命結(jié)果對比

    3.2 算例2

    層壓板的鋪層方式為[-60/0/60]3s,其在Abaqus建模時,材料屬性、模型尺寸、網(wǎng)格劃分、單元類型以及加載方式均與算例1一致,q=0.8. 施加載荷,當(dāng)P逐漸增加到320 MPa時:-60° 鋪層危險點的應(yīng)力分布為σ11=-160.3 MPa,σ22=64.17 MPa,τ12=29.98 MPa. 判斷該層是否失效:

    1.004 9>1

    (49)

    根據(jù)式(49)的計算結(jié)果大于1,因此-60°鋪層失效,失效強度比RL(-60°)=0.997 5,失效強度FS(-60°)=319.2 MPa,失效模式是基體破壞. 60°鋪層危險點的應(yīng)力分布為σ11=-160.9 MPa,σ22=64.14 MPa,τ12=29.92 MPa. 判斷該層是否失效:

    1.003 8>1

    (50)

    根據(jù)式(50)的計算結(jié)果大于1,因此60°鋪層失效,失效強度比RL(60°)=0.998 1,失效強度FS(60°)=319.4 MPa,失效模式是基體破壞. 0°鋪層危險點的應(yīng)力分布為σ11=1 235 MPa,σ22=8.705 MPa,τ12=20.93 MPa. 判斷該層是否失效:

    0.642 2<1

    (51)

    根據(jù)式(51)的計算結(jié)果小于1,因此0° 鋪層未失效. 對 -60° 和60° 鋪層,根據(jù)本文選定的剛度退化準(zhǔn)則進行剛度退化,然后施加載荷,當(dāng)載荷增加到350 MPa時:0° 鋪層危險點的應(yīng)力分布為σ11=1 558 MPa,σ22=10.06 MPa,τ12=26.25 MPa. 判斷該層是否失效:

    1.017 4>1

    (52)

    表8 不同方法預(yù)測[-60/0/60]3s層壓板的疲勞壽命結(jié)果對比

    4 結(jié)論

    復(fù)合材料層壓板內(nèi)部損傷會引起剛度退化,而損傷不斷累積會導(dǎo)致其疲勞破壞. 以疲勞模量為參量,借助灰色系統(tǒng)構(gòu)造FRP單向板的損傷函數(shù)和疲勞壽命函數(shù),通過建模對FRP層壓板有限元分析,對其疲勞破壞過程模擬. 這種理論與仿真結(jié)合預(yù)測FRP多向?qū)訅喊迤趬勖姆椒梢缘玫揭韵陆Y(jié)論:

    (1) 通過灰色系統(tǒng)GM(1, 1)模型建模,擬合了q-N函數(shù)表達式以及損傷函數(shù)表達式,且擬合結(jié)果與最小二乘法擬合的結(jié)果相比提高了預(yù)測精度,擬合的疲勞壽命和損傷函數(shù)可直接用于同種材料任意鋪層的多向?qū)訅喊宓膲勖A(yù)測.

    (2) 基于Abaqus模擬仿真,對多層板進行應(yīng)力分析和失效分析,然后與建立的q-N函數(shù)結(jié)合,以此來預(yù)測FRP多層板的疲勞壽命. 該方法與傳統(tǒng)的經(jīng)典層壓板理論做應(yīng)力分析的方法相比減小了誤差. 因為經(jīng)典層壓板理論認(rèn)為,單層板的應(yīng)力分布是均勻的,但事實上,受力點與受力方向都會使單層板存在危險點,而通過Abaqus建模,可以解決這個問題.

    (3) 本文理論與模擬仿真緊密結(jié)合的方法突破了以往完全依賴試件疲勞試驗或者只有仿真的預(yù)測方法. 通過算例驗證,證明了該方法的有效性,具有一定的工程應(yīng)用價值.

    附錄

    圖1 PANI/CF電極紅外光譜圖Fig.1 FTIR spectra of PANI/CF electrodes

    表1 PANI/CF電極紅外特征峰位置Tab.1 FTIR characteristic peaks of PANI/CF electrodes

    圖2 PANI/CF電極XPS測試結(jié)果Fig.2 XPS test results of PANI/CF electrodes

    表2 PANI/CF電極實測幅值Tab.2 Actual measurement amplitude of PANI/CF electrodes mV

    猜你喜歡
    層壓板鋪層單向
    碳纖維/PPS熱塑性單向預(yù)浸帶進入市場
    用“單向?qū)m排除法”解四宮數(shù)獨
    單向截止閥密封失效分析
    含表面裂紋的復(fù)合材料層壓板剩余強度的工程算法
    不同工藝參數(shù)對自動鋪帶碳纖維層壓板的性能影響
    準(zhǔn)靜態(tài)壓痕力作用下復(fù)合材料層壓板損傷分析方法
    衛(wèi)星天線復(fù)合材料框架的鋪層優(yōu)化設(shè)計
    復(fù)合材料軸結(jié)構(gòu)力學(xué)性能預(yù)測及鋪層方案設(shè)計
    CPD鋪層自動創(chuàng)建技術(shù)
    FIBERSIM軟件在工字梁鋪層設(shè)計上的應(yīng)用
    教練機(2015年2期)2015-04-03 06:51:48
    成人一区二区视频在线观看| 日韩一本色道免费dvd| 亚洲欧美成人综合另类久久久 | 亚洲欧美成人综合另类久久久 | 淫秽高清视频在线观看| 亚洲国产精品成人久久小说 | 九九在线视频观看精品| 日韩欧美一区二区三区在线观看| videossex国产| 亚洲欧美精品综合久久99| 中文字幕av在线有码专区| 精品一区二区三区人妻视频| 久久精品影院6| 欧美成人a在线观看| 成年女人毛片免费观看观看9| 免费一级毛片在线播放高清视频| 麻豆久久精品国产亚洲av| 一进一出好大好爽视频| 亚洲高清免费不卡视频| 观看免费一级毛片| 亚洲中文字幕一区二区三区有码在线看| 久久国产乱子免费精品| 欧美日韩乱码在线| 久久精品夜色国产| 国产伦在线观看视频一区| 亚洲激情五月婷婷啪啪| or卡值多少钱| 啦啦啦啦在线视频资源| 国产成人一区二区在线| 2021天堂中文幕一二区在线观| 女生性感内裤真人,穿戴方法视频| 男插女下体视频免费在线播放| 亚洲美女搞黄在线观看 | 久久九九热精品免费| 在线播放国产精品三级| 婷婷亚洲欧美| 丰满的人妻完整版| 免费看美女性在线毛片视频| 久久久久国内视频| 亚洲av中文字字幕乱码综合| 国产麻豆成人av免费视频| 久久午夜福利片| 最近2019中文字幕mv第一页| 一卡2卡三卡四卡精品乱码亚洲| 国产精品人妻久久久影院| 少妇丰满av| av中文乱码字幕在线| 国产女主播在线喷水免费视频网站 | 一进一出抽搐动态| 午夜激情欧美在线| 日本黄色片子视频| 成年av动漫网址| 国内精品久久久久精免费| 亚洲成av人片在线播放无| 又黄又爽又免费观看的视频| 免费av不卡在线播放| 我的老师免费观看完整版| av中文乱码字幕在线| 欧美一区二区精品小视频在线| 欧美区成人在线视频| 亚洲精品一区av在线观看| 97人妻精品一区二区三区麻豆| 深夜精品福利| 成人午夜高清在线视频| 国产成人a∨麻豆精品| 可以在线观看的亚洲视频| 干丝袜人妻中文字幕| 欧美潮喷喷水| 亚洲av二区三区四区| 久久久久久九九精品二区国产| 99久久久亚洲精品蜜臀av| 九九在线视频观看精品| 久久久久久久久久黄片| 国产真实乱freesex| 精品久久久噜噜| 99久久久亚洲精品蜜臀av| 亚洲精品粉嫩美女一区| 女人被狂操c到高潮| 国产伦在线观看视频一区| 国产真实伦视频高清在线观看| 在线观看美女被高潮喷水网站| a级毛色黄片| 欧美激情国产日韩精品一区| 亚洲精品日韩av片在线观看| 一级毛片我不卡| 国产高清三级在线| 悠悠久久av| 亚洲内射少妇av| 最近2019中文字幕mv第一页| 人妻制服诱惑在线中文字幕| avwww免费| 男女视频在线观看网站免费| 国产色爽女视频免费观看| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| 99热全是精品| 人妻制服诱惑在线中文字幕| 在线免费观看的www视频| 1024手机看黄色片| 亚洲人成网站在线观看播放| 在线观看一区二区三区| 免费观看的影片在线观看| 亚洲丝袜综合中文字幕| 久久国产乱子免费精品| 欧美xxxx黑人xx丫x性爽| 中文字幕人妻熟人妻熟丝袜美| 亚洲一级一片aⅴ在线观看| 成人无遮挡网站| 久久久久精品国产欧美久久久| av天堂中文字幕网| 天堂动漫精品| 欧美中文日本在线观看视频| 一级黄色大片毛片| 久久久久久久久久黄片| 中文字幕精品亚洲无线码一区| 亚洲成av人片在线播放无| 97人妻精品一区二区三区麻豆| 国产精品久久久久久av不卡| 美女被艹到高潮喷水动态| 日日啪夜夜撸| 日韩成人av中文字幕在线观看 | 18禁黄网站禁片免费观看直播| av卡一久久| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费| 国产国拍精品亚洲av在线观看| 最近2019中文字幕mv第一页| 亚洲色图av天堂| 成人国产麻豆网| 亚洲欧美清纯卡通| 欧美日韩国产亚洲二区| 亚洲精品成人久久久久久| 国产成人a∨麻豆精品| 亚洲图色成人| 两个人的视频大全免费| 欧美日韩乱码在线| 淫妇啪啪啪对白视频| 一区二区三区高清视频在线| 露出奶头的视频| 免费看av在线观看网站| 亚洲人成网站在线播放欧美日韩| 激情 狠狠 欧美| 午夜福利成人在线免费观看| 身体一侧抽搐| 中文字幕熟女人妻在线| 亚洲国产精品国产精品| 又粗又爽又猛毛片免费看| 免费大片18禁| 精品久久久久久久久久免费视频| 国产真实伦视频高清在线观看| 性插视频无遮挡在线免费观看| 淫妇啪啪啪对白视频| 午夜福利在线在线| 亚洲国产精品成人综合色| 国产成人a区在线观看| 又粗又爽又猛毛片免费看| av在线播放精品| 亚洲国产精品国产精品| 国产男靠女视频免费网站| 久久6这里有精品| 国产成人a区在线观看| 日韩欧美 国产精品| 亚洲18禁久久av| 久久久欧美国产精品| 极品教师在线视频| a级毛色黄片| 禁无遮挡网站| 一进一出抽搐动态| 亚洲欧美日韩高清专用| 久久草成人影院| 久久人妻av系列| 国产成人freesex在线 | 亚洲精品色激情综合| 两个人视频免费观看高清| 国产精品美女特级片免费视频播放器| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精华国产精华液的使用体验 | 国产午夜福利久久久久久| 男人舔女人下体高潮全视频| 日韩欧美免费精品| 少妇熟女欧美另类| 亚洲,欧美,日韩| 久久中文看片网| 在线看三级毛片| 网址你懂的国产日韩在线| 天堂动漫精品| 色尼玛亚洲综合影院| 12—13女人毛片做爰片一| 免费看美女性在线毛片视频| 日韩大尺度精品在线看网址| 精品人妻视频免费看| 亚洲婷婷狠狠爱综合网| 欧美激情国产日韩精品一区| 国产精品三级大全| 天堂影院成人在线观看| 亚洲欧美日韩无卡精品| 欧美色视频一区免费| 搞女人的毛片| 久久久精品94久久精品| 一本精品99久久精品77| 日韩中字成人| 少妇人妻精品综合一区二区 | 国产精品一区二区性色av| 听说在线观看完整版免费高清| 九九爱精品视频在线观看| 天天躁夜夜躁狠狠久久av| 校园人妻丝袜中文字幕| 久久久久久久亚洲中文字幕| 久久亚洲精品不卡| 成年免费大片在线观看| 国产老妇女一区| 久久久久久大精品| 成年av动漫网址| 中文字幕av成人在线电影| 国产高清有码在线观看视频| av专区在线播放| 国产不卡一卡二| 国产淫片久久久久久久久| 国产探花极品一区二区| 国产一区二区亚洲精品在线观看| 精品久久久久久久久久久久久| 亚洲精品日韩在线中文字幕 | 中文在线观看免费www的网站| 亚洲精品亚洲一区二区| 好男人在线观看高清免费视频| 精品少妇黑人巨大在线播放 | 日韩欧美三级三区| 色av中文字幕| 欧美日韩在线观看h| 久久久久九九精品影院| 国产白丝娇喘喷水9色精品| 精品久久国产蜜桃| 最近最新中文字幕大全电影3| 国产av在哪里看| 久久国内精品自在自线图片| av天堂中文字幕网| 亚洲欧美日韩卡通动漫| 国产精品亚洲美女久久久| 亚洲熟妇中文字幕五十中出| av黄色大香蕉| 3wmmmm亚洲av在线观看| 在线看三级毛片| 精品熟女少妇av免费看| 免费搜索国产男女视频| 国产伦一二天堂av在线观看| 国产精品爽爽va在线观看网站| 在线播放国产精品三级| 婷婷精品国产亚洲av在线| 乱人视频在线观看| 一级毛片aaaaaa免费看小| 少妇丰满av| 少妇熟女欧美另类| 一区二区三区高清视频在线| 亚洲最大成人手机在线| 黄色日韩在线| 狠狠狠狠99中文字幕| 成人一区二区视频在线观看| 香蕉av资源在线| 国产人妻一区二区三区在| 色综合色国产| 免费观看人在逋| 国产av麻豆久久久久久久| 国产精品一及| 亚洲精品456在线播放app| 天天躁夜夜躁狠狠久久av| 最后的刺客免费高清国语| 国产精品亚洲美女久久久| 国产成人福利小说| 久久久久久久久久久丰满| 男人舔奶头视频| 久久久久国内视频| 免费av不卡在线播放| 欧美日韩乱码在线| 日韩欧美在线乱码| 国产一区二区激情短视频| 老女人水多毛片| 搡老岳熟女国产| 精品人妻熟女av久视频| 国产一级毛片七仙女欲春2| 一区二区三区四区激情视频 | 国产三级在线视频| 国产真实乱freesex| 亚洲人成网站在线播放欧美日韩| 中国美女看黄片| 欧美+日韩+精品| 丰满乱子伦码专区| 成人一区二区视频在线观看| 免费高清视频大片| av在线老鸭窝| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av熟女| 久久久欧美国产精品| 十八禁网站免费在线| 天堂影院成人在线观看| 日韩欧美一区二区三区在线观看| 日本三级黄在线观看| 亚洲乱码一区二区免费版| 国产淫片久久久久久久久| 亚洲成人av在线免费| 成人午夜高清在线视频| 黄色配什么色好看| 狠狠狠狠99中文字幕| 国产亚洲精品av在线| 国产亚洲91精品色在线| 国产大屁股一区二区在线视频| 看非洲黑人一级黄片| 韩国av在线不卡| 久久久久久九九精品二区国产| 午夜福利视频1000在线观看| 插逼视频在线观看| 欧美激情在线99| 亚洲,欧美,日韩| 综合色av麻豆| 亚洲精华国产精华液的使用体验 | 少妇人妻精品综合一区二区 | 3wmmmm亚洲av在线观看| 欧美成人免费av一区二区三区| 最好的美女福利视频网| 熟女人妻精品中文字幕| 悠悠久久av| 国内精品宾馆在线| 久久精品国产亚洲av涩爱 | 99国产精品一区二区蜜桃av| 国产伦精品一区二区三区四那| 自拍偷自拍亚洲精品老妇| 久99久视频精品免费| 亚洲第一区二区三区不卡| 97热精品久久久久久| 白带黄色成豆腐渣| 日韩av在线大香蕉| 99久久九九国产精品国产免费| 男女之事视频高清在线观看| 欧美+日韩+精品| 欧美性猛交╳xxx乱大交人| 中文字幕av在线有码专区| 女生性感内裤真人,穿戴方法视频| 国产色婷婷99| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 久久精品国产清高在天天线| 国产在视频线在精品| 九九久久精品国产亚洲av麻豆| 国产91av在线免费观看| 男女做爰动态图高潮gif福利片| 91精品国产九色| 精品人妻视频免费看| 国产69精品久久久久777片| av天堂在线播放| 可以在线观看的亚洲视频| 国产高清视频在线播放一区| 好男人在线观看高清免费视频| 午夜福利在线观看免费完整高清在 | 男人的好看免费观看在线视频| 日本成人三级电影网站| 最后的刺客免费高清国语| 在线观看66精品国产| 亚洲五月天丁香| 国语自产精品视频在线第100页| 丝袜美腿在线中文| av福利片在线观看| 男插女下体视频免费在线播放| 99视频精品全部免费 在线| 国产精品久久视频播放| 欧美日韩一区二区视频在线观看视频在线 | 18禁在线播放成人免费| 日日摸夜夜添夜夜添av毛片| 久久综合国产亚洲精品| 久久久久国产网址| 国产精品久久久久久av不卡| 欧美日韩综合久久久久久| 嫩草影院入口| 亚洲经典国产精华液单| 欧美成人一区二区免费高清观看| 国产精品亚洲一级av第二区| 99久国产av精品国产电影| 不卡一级毛片| 简卡轻食公司| 99国产极品粉嫩在线观看| 国产日本99.免费观看| 人人妻人人澡欧美一区二区| 亚洲欧美精品综合久久99| 亚洲精品成人久久久久久| 人妻久久中文字幕网| 久久精品夜夜夜夜夜久久蜜豆| 高清毛片免费看| www.色视频.com| 精品久久久久久久久久久久久| 久久久久久伊人网av| 国产 一区 欧美 日韩| 深夜精品福利| 亚洲av中文字字幕乱码综合| 深爱激情五月婷婷| 日本熟妇午夜| 亚洲av熟女| 国产黄片美女视频| 禁无遮挡网站| 国产乱人视频| 国产69精品久久久久777片| 91av网一区二区| 亚洲人成网站在线播| 菩萨蛮人人尽说江南好唐韦庄 | 欧美色视频一区免费| 91久久精品电影网| 99热6这里只有精品| 最好的美女福利视频网| 午夜福利成人在线免费观看| 一进一出好大好爽视频| 成人综合一区亚洲| 国产一区二区三区在线臀色熟女| 欧美人与善性xxx| 亚洲欧美日韩卡通动漫| 亚洲人成网站在线播放欧美日韩| 国产真实乱freesex| 亚洲精品日韩在线中文字幕 | 久久鲁丝午夜福利片| 69av精品久久久久久| 亚洲av五月六月丁香网| 久久久久国内视频| 一区二区三区高清视频在线| 久久99热6这里只有精品| 搡老岳熟女国产| 亚洲不卡免费看| 99热网站在线观看| 国产精品久久久久久av不卡| 精品熟女少妇av免费看| 午夜激情欧美在线| av国产免费在线观看| 别揉我奶头 嗯啊视频| 香蕉av资源在线| 天堂网av新在线| 国产黄a三级三级三级人| 午夜日韩欧美国产| 亚洲av熟女| 在线观看免费视频日本深夜| 亚洲欧美中文字幕日韩二区| 亚洲激情五月婷婷啪啪| 日本黄色片子视频| 九九久久精品国产亚洲av麻豆| 男女那种视频在线观看| 亚洲国产欧美人成| 欧美+日韩+精品| 男女视频在线观看网站免费| 日本黄大片高清| 午夜激情福利司机影院| 美女大奶头视频| 日本在线视频免费播放| 亚洲熟妇熟女久久| 亚洲性久久影院| 成人午夜高清在线视频| 欧美三级亚洲精品| 狂野欧美白嫩少妇大欣赏| 国产三级在线视频| 久久久久久久久久久丰满| 波多野结衣高清作品| 天堂av国产一区二区熟女人妻| 久久久久性生活片| 97人妻精品一区二区三区麻豆| 别揉我奶头 嗯啊视频| 波多野结衣高清作品| 欧美高清成人免费视频www| 午夜亚洲福利在线播放| 国产精品1区2区在线观看.| 日本欧美国产在线视频| 人人妻人人澡欧美一区二区| 欧美国产日韩亚洲一区| 亚洲av成人av| 久久亚洲国产成人精品v| 69人妻影院| 国产高清不卡午夜福利| 国产美女午夜福利| 日本色播在线视频| 日本一本二区三区精品| 淫秽高清视频在线观看| 国产白丝娇喘喷水9色精品| 久久午夜亚洲精品久久| 日日啪夜夜撸| 日本三级黄在线观看| 久久人人爽人人爽人人片va| 搞女人的毛片| 日本一二三区视频观看| 人人妻人人澡欧美一区二区| 尤物成人国产欧美一区二区三区| 色综合站精品国产| 成年版毛片免费区| 久久午夜亚洲精品久久| 日韩欧美精品v在线| 日本在线视频免费播放| 国产精品日韩av在线免费观看| 国产成人freesex在线 | 久久人人精品亚洲av| 久久精品国产鲁丝片午夜精品| 精品午夜福利在线看| 尾随美女入室| 一级毛片我不卡| 亚洲激情五月婷婷啪啪| 亚洲不卡免费看| 床上黄色一级片| 免费看光身美女| 99久久精品国产国产毛片| 国产av在哪里看| 在线播放国产精品三级| 丰满乱子伦码专区| 国语自产精品视频在线第100页| 99久久精品国产国产毛片| 99国产极品粉嫩在线观看| 亚洲精品乱码久久久v下载方式| 日本精品一区二区三区蜜桃| 尾随美女入室| 插阴视频在线观看视频| 日韩精品青青久久久久久| 女的被弄到高潮叫床怎么办| 免费高清视频大片| 人妻制服诱惑在线中文字幕| 日韩亚洲欧美综合| 精华霜和精华液先用哪个| 久久久精品大字幕| 嫩草影院新地址| 亚洲中文字幕日韩| 欧美3d第一页| 五月玫瑰六月丁香| 国产成人一区二区在线| 春色校园在线视频观看| 亚洲真实伦在线观看| 久久精品国产亚洲av天美| 少妇的逼好多水| 又黄又爽又免费观看的视频| 在线免费十八禁| 97人妻精品一区二区三区麻豆| 99热这里只有是精品50| 综合色丁香网| 久久精品91蜜桃| 欧美3d第一页| 国产精品一区二区三区四区久久| 欧美+日韩+精品| 毛片女人毛片| av国产免费在线观看| 大又大粗又爽又黄少妇毛片口| 一级毛片久久久久久久久女| 日本-黄色视频高清免费观看| a级毛色黄片| 美女黄网站色视频| 国产免费男女视频| 国产片特级美女逼逼视频| 国产 一区精品| 人妻丰满熟妇av一区二区三区| 亚洲七黄色美女视频| 国产女主播在线喷水免费视频网站 | 此物有八面人人有两片| 尤物成人国产欧美一区二区三区| 乱人视频在线观看| 亚洲电影在线观看av| 久久精品国产99精品国产亚洲性色| 国产精品免费一区二区三区在线| 熟妇人妻久久中文字幕3abv| 99在线视频只有这里精品首页| 日本色播在线视频| 成人一区二区视频在线观看| 国产精品久久久久久av不卡| 国产精品久久久久久亚洲av鲁大| 国产aⅴ精品一区二区三区波| 我要搜黄色片| 亚洲成人久久性| 精品久久久久久成人av| 日韩强制内射视频| 小说图片视频综合网站| 性欧美人与动物交配| 久久人人爽人人爽人人片va| 三级经典国产精品| 免费在线观看成人毛片| 99热只有精品国产| 精品少妇黑人巨大在线播放 | 春色校园在线视频观看| 美女免费视频网站| 99热精品在线国产| 精品久久久久久久末码| 亚洲精品在线观看二区| 国产又黄又爽又无遮挡在线| 一夜夜www| 久久热精品热| 尤物成人国产欧美一区二区三区| 久久综合国产亚洲精品| 99久久久亚洲精品蜜臀av| 日本三级黄在线观看| 成人漫画全彩无遮挡| 日本黄大片高清| 蜜桃久久精品国产亚洲av| 欧美性猛交黑人性爽| 你懂的网址亚洲精品在线观看 | 久久99热6这里只有精品| 国产精品一区二区三区四区久久| 欧美不卡视频在线免费观看| 韩国av在线不卡| 国产高清视频在线观看网站| 国产白丝娇喘喷水9色精品| 麻豆精品久久久久久蜜桃| 日本精品一区二区三区蜜桃| 亚洲色图av天堂| 久久99热6这里只有精品| 久久人人爽人人爽人人片va| 国产精品1区2区在线观看.| 夜夜夜夜夜久久久久| 变态另类丝袜制服| 国产老妇女一区| 欧美成人精品欧美一级黄| 欧美性猛交黑人性爽| 成年女人看的毛片在线观看| 免费黄网站久久成人精品| or卡值多少钱| 国产精品一区二区性色av| 直男gayav资源| 国产黄色小视频在线观看| 91午夜精品亚洲一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产高清有码在线观看视频| 欧美又色又爽又黄视频|