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

    基于BFD和MSCNN的風(fēng)電滾動軸承智能故障診斷

    2021-06-10 00:25:50鄧敏強(qiáng)鄧艾東史曜煒馬天霆
    關(guān)鍵詞:階次準(zhǔn)確率卷積

    鄧敏強(qiáng) 鄧艾東 朱 靜 史曜煒 馬天霆,3

    (1東南大學(xué)能源與環(huán)境學(xué)院, 南京 210096)(2東南大學(xué)火電機(jī)組振動國家工程研究中心, 南京 210096)(3 國華太倉發(fā)電有限公司, 蘇州 215000)

    滾動軸承作為風(fēng)電傳動鏈的關(guān)鍵零部件[1],其健康狀態(tài)的準(zhǔn)確識別對提高風(fēng)電機(jī)組運行的安全性與經(jīng)濟(jì)性意義重大[2].受時變風(fēng)況的影響,風(fēng)電機(jī)組常處于變速變載荷條件下運行[3].針對變工況下滾動軸承的故障診斷問題,國內(nèi)外學(xué)者開展了大量的研究工作[4].其中,以振動信號分析為基礎(chǔ)的診斷方法受到了廣泛關(guān)注,其主要根據(jù)包絡(luò)譜中的故障特征頻率實現(xiàn)滾動軸承健康狀態(tài)的有效識別[5].然而,該方法一定程度上依賴專家經(jīng)驗,難以實現(xiàn)設(shè)備故障的自動診斷[6].

    隨著機(jī)器學(xué)習(xí)和人工智能技術(shù)的發(fā)展,數(shù)據(jù)驅(qū)動的旋轉(zhuǎn)機(jī)械智能故障診斷受到了越來越多的重視.唐貴基等[7]建立了一種基于快速譜相關(guān)和粒子群優(yōu)化支持向量機(jī)的故障診斷模型,趙小強(qiáng)等[8]提出一種基于改進(jìn)AlexNet的滾動軸承變工況故障診斷方法,施杰等[9]提出了一種變分模態(tài)分解與深度遷移想結(jié)合的診斷方法.已有的研究成果主要以特征參數(shù)、原始時域波形和傅里葉譜等為特征量,其易受轉(zhuǎn)速、載荷等運行條件的影響,一定程度上降低了模型在不同工況下的泛化能力.因此,已有的智能診斷方法在模型構(gòu)建過程中都需要有標(biāo)簽或無標(biāo)簽訓(xùn)練樣本覆蓋完整的運行工況.然而,受現(xiàn)場采集條件的限制,難以獲得覆蓋全工況的訓(xùn)練樣本.因此,僅根據(jù)部分工況下的數(shù)據(jù)構(gòu)建適用于全工況的智能診斷模型有待深入研究.

    針對該問題,本文從包絡(luò)階次譜(EOS)智能識別的角度出發(fā),提出了一種基于帶寬傅里葉分解(BFD)和多尺度卷積神經(jīng)網(wǎng)絡(luò)(MSCNN)的風(fēng)電滾動軸承故障診斷方法.BFD根據(jù)傅里葉譜能量分布將原始信號分解為窄帶的帶寬模態(tài)函數(shù)(BMFs),能有效地從振動信號中分離出包含故障特征的窄帶分量[10].然后,提出了一種希爾伯特階次變換(HOT)算法并將其應(yīng)用于BMF的包絡(luò)階次譜計算.最后,根據(jù)特征階次比選擇分解結(jié)果中包含故障信息最多的有效分量,并通過MSCNN學(xué)習(xí)其包絡(luò)階次譜與故障類別之間的映射關(guān)系以實現(xiàn)滾動軸承健康狀態(tài)的自動識別.實驗結(jié)果表明,本文所提方法以BFD分解結(jié)果有效分量的包絡(luò)階次譜作為故障識別的特征量,可有效提高模型的泛化能力,降低了對樣本量的要求,能準(zhǔn)確識別出變工況下滾動軸承的健康狀態(tài).

    1 算法原理

    1.1 帶寬傅里葉分解

    假設(shè)時域信號x(t)是定義在區(qū)間[a,b]上的能量有限信號,BFD的目標(biāo)是將時域信號x(t)通過傅里葉變換轉(zhuǎn)換成頻域信號X(ω)后,在頻域內(nèi)將其自適應(yīng)分解為K個稀疏的以{ωk}為中心的窄帶子信號{Uk(ω)},然后通過傅里葉逆變換將{Uk(ω)}轉(zhuǎn)換為帶寬模態(tài)函數(shù){uk(t)}.在BFD算法中,子信號的稀疏性和窄帶性通過對總帶寬的優(yōu)化實現(xiàn),其構(gòu)造的約束性優(yōu)化問題為[11]

    L0({Uk(ω)},{ωk})=

    (1)

    采用拉格朗日乘數(shù)法將式(1)中的約束性優(yōu)化問題轉(zhuǎn)換為如下無約束性求極值問題:

    L1({Uk(ω)},{ωk})=

    (2)

    通過對式(2)求鞍點即可得到Uk(ω)的解析式:

    (3)

    將式(3)代入式(2),原始優(yōu)化問題L0({Uk(ω)},{ωk})可簡化為中心頻率ωk的搜索問題L′1({ωk}):

    (4)

    其中

    (5)

    通過對L1({Uk(ω)},{ωk})求鞍點可得中心頻率ωk滿足的隱式方程:

    (6)

    式(6)迭代計算至其收斂后,對{Uk(ω)}求傅里葉逆變換即可得到時域信號x(t)的BFD分解結(jié)果{uk(t)}.

    由式(1)中構(gòu)造的優(yōu)化問題可知,中心頻率的數(shù)量K對BFD的分解結(jié)果有較大的影響.K過小會導(dǎo)致信號的欠分解,從而出現(xiàn)模態(tài)混疊現(xiàn)象,K過大可能會導(dǎo)致特征信息被分解在2個或以上的BMF之中,從而降低算法的特征提取能力.因此,BFD算法中心頻率數(shù)量K的確定方法如下:

    ①令K=2.

    ②對原始信號x(t)進(jìn)行BFD運算.

    ④令K=K+1,執(zhí)行步驟②.

    其中,步驟③通過BMF中心頻率的分布判斷BFD是否出現(xiàn)過度分解.在BFD算法中μ取0.1,即2個相鄰的BMF疊加為一個分量后仍滿足窄帶信號的條件時,算法出現(xiàn)過度分解.

    1.2 希爾伯特階次變換

    傳統(tǒng)的階次分析以角域重采樣為基礎(chǔ),然而該方法需要大量的插值運算,一定程度上增加了計算難度.本文在希爾伯特變換的基礎(chǔ)上提出了希爾伯特階次變換算法,可通過積分運算直接計算出窄帶分量的包絡(luò)階次譜,有效降低了計算的復(fù)雜性.

    假設(shè)有一調(diào)頻信號y(t),其解析式為

    (7)

    式中,A為幅值;ω1(t)為角頻率;fr(t)為時變轉(zhuǎn)頻;o1為特征階次.其階次譜計算式為

    (8)

    式中,T為時間長度.文獻(xiàn)[12]的研究驗證了下式的正確性:

    (9)

    則式(8)可進(jìn)一步簡化為

    (10)

    式(10)表明,根據(jù)式(8)計算得到的階次譜在特征階次o1處存在一條突出的譜線,即該算法能有效應(yīng)用于特征階次的檢測.因此,對于任意窄帶信號x(t),其包絡(luò)階次譜F(o)的計算式可寫為

    (11)

    式中,xe(t)為信號x(t)的包絡(luò)函數(shù),根據(jù)希爾伯特變換得

    (12)

    1.3 多尺度卷積神經(jīng)網(wǎng)絡(luò)

    與傳統(tǒng)機(jī)器學(xué)習(xí)算法相比,深度學(xué)習(xí)模型能更有效地從高維復(fù)雜輸入中自動提取特征.卷積神經(jīng)網(wǎng)絡(luò)(CNN)是應(yīng)用最廣泛的深度學(xué)習(xí)模型之一,通過卷積、池化等操作提取原始數(shù)據(jù)的特征,并通過權(quán)連接層輸出模型的計算結(jié)果.其中,卷積核的大小在一定程度上影響著特征提取的效果和模型的故障識別能力.MSCNN是一種改進(jìn)的卷積神經(jīng)網(wǎng)絡(luò),通過不同大小的卷積核從多尺度挖掘特征信息,有效解決了傳統(tǒng)CNN模型卷積核的自適應(yīng)選擇問題[13].在本文所提方法中,模型輸入為包絡(luò)階次譜,其階次范圍為0~50,階次分辨率為0.05,輸入特征的維度為1 000,輸出的類別數(shù)為4.圖1給出了本文構(gòu)建的MSCNN模型網(wǎng)絡(luò)結(jié)構(gòu),其參數(shù)設(shè)置如表1所示.

    圖1 MSCNN網(wǎng)絡(luò)結(jié)構(gòu)圖

    表1 MSCNN網(wǎng)絡(luò)參數(shù)

    首先,考慮到轉(zhuǎn)速測量誤差導(dǎo)致的階次偏差,MSCNN通過池化層提取池化范圍內(nèi)的最大值以降低階次偏差對故障識別的影響;然后,利用不同大小的卷積核從多尺度提取池化輸出的特征,并對其進(jìn)行線性加權(quán);最后,通過全連接層和Softmax激活函數(shù)計算模型輸入屬于不同類別的概率,并根據(jù)最大概率對設(shè)備健康狀態(tài)進(jìn)行分類.該模型的正向傳播過程如下:

    (13)

    式中,x為模型輸入;pc為分類器輸出的類別概率;pool(·)表示池化操作,其輸出為y;?為卷積運算,其通過ci線性加權(quán)后的輸出為z;Wi和bi分別為不同尺度的卷積核及其偏置;W為權(quán)連接矩陣;f(·)為偏置是bi的Softmax激活函數(shù).

    2 故障診斷方法

    變工況下風(fēng)電滾動軸承的故障特征不僅受強(qiáng)背景噪聲的干擾,而且受到轉(zhuǎn)速、扭矩等時變運行參數(shù)的影響.本文通過信號分解技術(shù)從原始數(shù)據(jù)中分離出包含故障信息的單分量信號,并通過包絡(luò)階次分析降低變工況對故障特征的影響,最后通過MSCNN自動識別設(shè)備的健康狀態(tài),具體的方法步驟如下:

    ①提取包絡(luò)階次譜.通過BFD算法將原始振動信號分解為窄帶的BMF,并通過希爾伯特階次變換計算各BMF的包絡(luò)階次譜.

    ②有效分量選擇.信號分解后僅部分BMF中包含相對較多的故障信息.參考文獻(xiàn)[10]的研究結(jié)果,本文根據(jù)特征階次比(COR)挑選出分解結(jié)果中包含故障信息最多的有效分量,其計算式為

    (14)

    式中,F(xiàn)為BMF的包絡(luò)階次譜,其階次范圍為0~oend;M為特征階次的最大倍數(shù),本文取3.由式(14)可知,COR反映了階次譜中最突出階次的比例,本文取COR最大的BMF作為分解結(jié)果的有效分量,并將其包絡(luò)階次譜作為MSCNN模型的輸入.

    ③模型訓(xùn)練.將部分工況下的數(shù)據(jù)作為訓(xùn)練樣本,構(gòu)建智能診斷模型.

    ④模型測試.對未知工況下的測試樣本,通過BFD算法和階次分析提取包絡(luò)階次譜,并通過訓(xùn)練好的MSCNN自動識別設(shè)備健康狀態(tài).

    3 實驗驗證

    3.1 數(shù)據(jù)集

    實驗數(shù)據(jù)采用國外公開數(shù)據(jù)集,以模擬風(fēng)電傳動鏈高速端滾動軸承的運行工況.

    數(shù)據(jù)集Ⅰ是美國凱斯西儲大學(xué)軸承數(shù)據(jù)中心提供的滾動軸承振動加速度數(shù)據(jù),其實驗臺如圖2(a)示,主要由驅(qū)動電機(jī)、扭矩傳感器、編碼器和測力計等部件組成.本文使用的實驗數(shù)據(jù)采集自電機(jī)驅(qū)動端的6205-2RS-JEM SKF型深溝球滾動軸承,信號采樣頻率為12 kHz,其故障特征階次如表2所示.數(shù)據(jù)集Ⅰ包括4種健康狀態(tài),分別為正常、內(nèi)圈故障、外圈故障和滾動體故障,其信號時域波形如圖3(a)所示.為模擬不同的運行工況,實驗分別在4種工況下進(jìn)行:①轉(zhuǎn)速1 797 r/min,負(fù)載0;②轉(zhuǎn)速1 772 r/min,負(fù)載0.75 kW;③轉(zhuǎn)速1 750 r/min,負(fù)載 1.49 kW;④轉(zhuǎn)速1 730 r/min,負(fù)載2.24 kW.表3給出了數(shù)據(jù)集Ⅰ的樣本分布,其中信號的時間長度為0.5 s,每種工況下各類別的樣本數(shù)為117,總樣本數(shù)為1 872.

    數(shù)據(jù)集Ⅱ是加拿大渥太華大學(xué)提供的滾動軸承振動加速度數(shù)據(jù),其實驗臺如圖2(b)所示.測試軸承的型號為ER16K,數(shù)據(jù)集主要包括正常、內(nèi)圈故障和外圈故障3種類別的振動信號,其故障特征階次如表2所示.實驗在勻變速工況下進(jìn)行,轉(zhuǎn)速的變化范圍為900~1 620 r/min,信號的采樣頻率為200 kHz.數(shù)據(jù)集Ⅱ的樣本分布如表3所示,其中信號的時間長度為0.5 s,每種工況下各類別的樣本數(shù)為60,總樣本數(shù)為720.

    (a) 數(shù)據(jù)集Ⅰ實驗臺

    (b) 數(shù)據(jù)集Ⅱ?qū)嶒炁_

    表2 滾動軸承故障特征階次

    表3 數(shù)據(jù)集樣本

    3.2 對比方法

    用于對比的故障診斷方法主要包括以下5種:

    1)本文所提方法,即將BFD分解結(jié)果中有效分量的包絡(luò)階次譜作為模型輸入,以MSCNN為分類器;將該方法記為BFD-MSCNN.

    2)將原始信號的包絡(luò)階次譜作為MSCNN的輸入,以驗證BFD算法提取原始信號中有效分量的必要性;將該方法記為EOS-MSCNN.

    3)將原始振動信號作為模型輸入,通過多核殘差卷積神經(jīng)網(wǎng)絡(luò)(ResCNN)構(gòu)建診斷模型[14];將該方法記為waveform-ResCNN.

    4)將振動信號的傅里葉譜作為模型輸入,通過多注意力一維卷積神經(jīng)網(wǎng)絡(luò)(MACNN) 構(gòu)建診斷模型[15];將該方法記為spectrum-MACNN.

    5)將振動信號的特征參數(shù)作為模型輸入,以支持向量機(jī)(SVM)為分類器.將該方法記為featureSVM.

    其中,本文方法和EOS-MSCNN方法中無超參數(shù)優(yōu)化;waveform-ResCNN、spectrum-MACNN和feature-SVM方法分別以原始時域信號、傅里葉譜及特征參數(shù)為模型輸入,以驗證本文所提方法中以BFD分解結(jié)果的包絡(luò)階次譜作為故障識別特征量的優(yōu)越性.waveform-ResCNN和spectrum-MACNN方法中振動信號的時間長度取0.5 s,神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)可參閱文獻(xiàn)[14-15];feature-SVM方法中提取的特征參數(shù)主要包括6個時域特征(均方根、峭度、偏度、波形因子、峰值因子和標(biāo)準(zhǔn)差)、5個頻域特征(中心頻率、均方頻率、均方根頻率、頻率方差、頻譜散度)和4個時頻特征(包絡(luò)階次譜中旋轉(zhuǎn)階次和故障特征階次處的幅值),SVM的超參數(shù)通過貝葉斯優(yōu)化算法搜索最優(yōu)值.

    3.3 實驗結(jié)果

    3.3.1 故障信號分析

    以數(shù)據(jù)集Ⅱ中內(nèi)圈故障為例,說明本文所提方法中BFD算法和HOT算法在故障特征提取方面的有效性.圖3(a)和(b)分別給了BFD分解結(jié)果的波形圖和包絡(luò)階次譜.由圖3(b)可知,BMF5~BMF9的包絡(luò)階次譜中能清楚地觀察到內(nèi)圈故障特征階次(5.45)及其倍數(shù),即本文所采用的信號處理方法能有效地提取出原始信號中的故障特征.表4給出了各BMF特征階次比的計算結(jié)果.由表可知,BMF9的特征階次比高于其他分量,即BMF9的故障特征相對較突出.因此,取BMF9為有效分量,并將其包絡(luò)階次譜作為MSCNN的輸入.

    (a) 時域波形

    表4 分解結(jié)果的特征階次比

    3.3.2 全工況交叉驗證

    在智能故障診斷研究中,通常將不同工況下的歷史數(shù)據(jù)隨機(jī)劃分為訓(xùn)練集和測試集以驗證模型的有效性.圖4給出了各診斷方法10次三折交叉驗證的平均準(zhǔn)確率.由圖可知,本文所提出的BFD-MSCNN方法、EOS-MSCNN、waveform-ResCNN和spectrum-MACNN的測試準(zhǔn)確率都超過99%,即深度學(xué)習(xí)模型能有效挖掘并學(xué)習(xí)故障信息,實現(xiàn)設(shè)備健康狀態(tài)的準(zhǔn)確識別.以特征參數(shù)為模型輸入的feature-SVM方法準(zhǔn)確率約為95%,略低于其他方法,表明該方法對故障信息的表征與識別能力略低于其他4種方法.圖4中結(jié)果還表明,訓(xùn)練樣本覆蓋全工況的情況下,有監(jiān)督深度學(xué)習(xí)和機(jī)器學(xué)習(xí)模型的故障識別效果均比較理想.

    圖4 三折交叉驗證準(zhǔn)確率

    3.3.3 變工況故障識別

    實際應(yīng)用中往往難以獲得覆蓋完整運行工況的故障樣本,有限的訓(xùn)練數(shù)據(jù)對模型在變工況場景下的泛化能力提出了更高的要求.如表5所示,本文依據(jù)訓(xùn)練樣本所屬工況的不同,將數(shù)據(jù)集Ⅰ和數(shù)據(jù)集Ⅱ分別劃分為4個診斷任務(wù),每個診斷任務(wù)中訓(xùn)練樣本與測試樣本均采集自不同的工況.

    表5 診斷任務(wù)設(shè)置

    表6和表7分別給出了數(shù)據(jù)集Ⅰ和數(shù)據(jù)集Ⅱ的測試準(zhǔn)確率.由表可知,在訓(xùn)練樣本所屬工況與測試樣本完全不同的情況下,BFD-MSCNN仍然有97%以上的平均準(zhǔn)確率,即本文所提方法僅根據(jù)部分工況下的監(jiān)測數(shù)據(jù)即可準(zhǔn)確地識別出設(shè)備在其他工況下的健康狀態(tài).EOS-MSCNN的測試準(zhǔn)確率明顯低于BFD-MSCNN,驗證了基于BFD的信號分解與有效分量選擇對故障識別的有效性.waveform-ResCNN、spectrum-MACNN和feature-SVM在數(shù)據(jù)集Ⅰ和數(shù)據(jù)集Ⅱ上的平均準(zhǔn)確率低于EOS-MSCNN,驗證了本文所提方法中以包絡(luò)階次譜為故障識別特征量的優(yōu)越性.

    表6 數(shù)據(jù)集Ⅰ測試準(zhǔn)確率

    表7 數(shù)據(jù)集Ⅱ測試準(zhǔn)確率

    3.3.4 小樣本分析

    圖5給出了本文所提方法在變工況故障識別中的平均準(zhǔn)確率隨各類別樣本數(shù)的變化趨勢.由圖可知,BFD-MSCNN僅需10個以上類別的訓(xùn)練樣本即可獲得95%以上的測試準(zhǔn)確率,即該方法能有效地應(yīng)用于小樣本情況下的故障診斷.

    圖5 BFD-MSCNN測試準(zhǔn)確率隨樣本數(shù)變化趨勢

    3.3.5 特征可視化分析

    為進(jìn)一步分析本文所提方法的有效性和優(yōu)越性,取數(shù)據(jù)集Ⅰ變工況故障識別中的診斷任務(wù)A進(jìn)行特征可視化分析,其結(jié)果如圖6所示.其中,前4種方法采用文獻(xiàn)[14]的研究思路,取神經(jīng)網(wǎng)絡(luò)最后一層的輸入,通過T-SNE進(jìn)行特征分布分析.由圖6(a)可知,在本文所提BFD-MSCNN方法中,同類別樣本的特征聚集性較好,不同類別的樣本特征混淆較小,具有較強(qiáng)的可識別性;此外,訓(xùn)練樣本與測試數(shù)據(jù)的分布較接近,即該方法具有較強(qiáng)的跨工況泛化能力.由圖6(b)~(e)可知,EOS-MSCNN和feature-SVM方法出現(xiàn)了一定程度的特征混淆,waveform-ResCNN和spectrum-MACNN方法中訓(xùn)練樣本和測試數(shù)據(jù)的特征分布存在較大的差異,一定程度上降低了故障識別能力.圖6(f)給出了小樣本分析的特征可視化結(jié)果,其中各類別的樣本數(shù)為10.圖6(f)中同類別特征的聚集性和不同類別樣本的可識別性進(jìn)一步驗證了本文所提方法在小樣本識別中的有效性.

    (a) BFD-MSCNN

    (d) spectrum-MACNN

    4 結(jié)論

    1)本文所提方法僅根據(jù)部分工況下的訓(xùn)練數(shù)據(jù)即可準(zhǔn)確識別變轉(zhuǎn)速滾動軸承的健康狀態(tài),其平均準(zhǔn)確率達(dá)到97%以上.

    2)本文所提方法對訓(xùn)練樣本的依賴性較小,在小樣本情況下也有較高的識別率,僅需10個以上類別訓(xùn)練樣本即可獲得95%以上的測試準(zhǔn)確率.

    3)通過BFD信號分解提取振動信號的有效分量,并從包絡(luò)階次譜自動識別的角度構(gòu)建智能診斷模型,有效降低了變工況對故障特征的影響,也降低了模型對樣本的依賴.

    猜你喜歡
    階次準(zhǔn)確率卷積
    基于3D-Winograd的快速卷積算法設(shè)計及FPGA實現(xiàn)
    乳腺超聲檢查診斷乳腺腫瘤的特異度及準(zhǔn)確率分析
    健康之家(2021年19期)2021-05-23 11:17:39
    不同序列磁共振成像診斷脊柱損傷的臨床準(zhǔn)確率比較探討
    2015—2017 年寧夏各天氣預(yù)報參考產(chǎn)品質(zhì)量檢驗分析
    階次分析在驅(qū)動橋異響中的應(yīng)用
    基于Vold-Kalman濾波的階次分析系統(tǒng)設(shè)計與實現(xiàn)*
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    高速公路車牌識別標(biāo)識站準(zhǔn)確率驗證法
    基于傅里葉域卷積表示的目標(biāo)跟蹤算法
    基于齒輪階次密度優(yōu)化的變速器降噪研究
    價值工程(2017年28期)2018-01-23 20:48:29
    国精品久久久久久国模美| 亚洲最大成人中文| 中文乱码字字幕精品一区二区三区| xxx大片免费视频| 街头女战士在线观看网站| 黄色日韩在线| 国产在线男女| 国产免费又黄又爽又色| 麻豆乱淫一区二区| 特大巨黑吊av在线直播| 赤兔流量卡办理| 人妻 亚洲 视频| 你懂的网址亚洲精品在线观看| 婷婷色麻豆天堂久久| 26uuu在线亚洲综合色| 国产欧美另类精品又又久久亚洲欧美| 日韩欧美精品免费久久| 超碰97精品在线观看| 女人被狂操c到高潮| 国产一区二区在线观看日韩| 免费人成在线观看视频色| 国产色爽女视频免费观看| 亚洲人成网站高清观看| 国产精品99久久99久久久不卡 | 嘟嘟电影网在线观看| 麻豆乱淫一区二区| 波多野结衣巨乳人妻| 大话2 男鬼变身卡| 亚洲自偷自拍三级| 美女内射精品一级片tv| 久久久久久久久久久丰满| 人妻 亚洲 视频| 人人妻人人澡人人爽人人夜夜| 男人舔奶头视频| 国产精品国产三级国产av玫瑰| 国产一区亚洲一区在线观看| 18禁在线无遮挡免费观看视频| 国产美女午夜福利| 欧美一级a爱片免费观看看| 国产探花极品一区二区| 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区国产| 亚洲va在线va天堂va国产| 啦啦啦中文免费视频观看日本| 成人特级av手机在线观看| 久久久久久国产a免费观看| 久久久国产一区二区| 美女国产视频在线观看| 少妇丰满av| 最近最新中文字幕免费大全7| 超碰av人人做人人爽久久| 国产亚洲av嫩草精品影院| 久久久久精品性色| 欧美97在线视频| 国产一区二区三区综合在线观看 | 亚洲av二区三区四区| 日韩一本色道免费dvd| 久久久成人免费电影| 亚洲精品国产成人久久av| 国产黄色视频一区二区在线观看| 国产淫语在线视频| 国产熟女欧美一区二区| 永久网站在线| 中文资源天堂在线| 久久久久久久精品精品| 亚洲成人精品中文字幕电影| 99精国产麻豆久久婷婷| 色网站视频免费| 夫妻性生交免费视频一级片| www.av在线官网国产| 午夜福利网站1000一区二区三区| 嫩草影院精品99| 好男人在线观看高清免费视频| 久久精品国产亚洲av天美| 国产欧美亚洲国产| 日本午夜av视频| 一级a做视频免费观看| 精品久久久精品久久久| 国产免费又黄又爽又色| 国产爱豆传媒在线观看| 大片电影免费在线观看免费| 免费在线观看成人毛片| 小蜜桃在线观看免费完整版高清| 男女啪啪激烈高潮av片| 有码 亚洲区| 国产黄色视频一区二区在线观看| 波多野结衣巨乳人妻| 亚洲欧美一区二区三区黑人 | 少妇高潮的动态图| 色哟哟·www| 少妇人妻一区二区三区视频| 成年女人看的毛片在线观看| 又爽又黄a免费视频| a级毛色黄片| 三级国产精品片| 在线观看免费高清a一片| 99九九线精品视频在线观看视频| 一级毛片久久久久久久久女| 国产综合精华液| 最近手机中文字幕大全| 国产在线一区二区三区精| 国产探花在线观看一区二区| 91久久精品电影网| 国产男人的电影天堂91| 久久久色成人| 久久久久久久亚洲中文字幕| 国产精品一区二区性色av| 99热国产这里只有精品6| 在线观看免费高清a一片| 亚洲国产成人一精品久久久| 欧美最新免费一区二区三区| 2018国产大陆天天弄谢| 午夜爱爱视频在线播放| 欧美变态另类bdsm刘玥| 黄片wwwwww| 免费观看性生交大片5| 麻豆精品久久久久久蜜桃| 小蜜桃在线观看免费完整版高清| 欧美bdsm另类| 亚洲,一卡二卡三卡| 欧美3d第一页| 伦理电影大哥的女人| 免费观看a级毛片全部| 亚洲av电影在线观看一区二区三区 | 午夜福利视频精品| 国产亚洲av片在线观看秒播厂| av在线老鸭窝| 高清午夜精品一区二区三区| 97在线视频观看| 婷婷色综合www| 亚洲人与动物交配视频| 香蕉精品网在线| 赤兔流量卡办理| 天天一区二区日本电影三级| 1000部很黄的大片| 日韩三级伦理在线观看| 国产精品嫩草影院av在线观看| 亚洲欧美一区二区三区国产| 欧美+日韩+精品| 国产精品偷伦视频观看了| 久久久久久伊人网av| 欧美性感艳星| 韩国高清视频一区二区三区| 亚洲最大成人中文| 日本与韩国留学比较| 亚洲精品色激情综合| 国产精品久久久久久av不卡| 亚洲av在线观看美女高潮| 三级经典国产精品| 99久久精品热视频| 啦啦啦在线观看免费高清www| 春色校园在线视频观看| 秋霞伦理黄片| 国产黄色免费在线视频| 免费av毛片视频| 一级毛片久久久久久久久女| 丝瓜视频免费看黄片| av福利片在线观看| 高清日韩中文字幕在线| 国内精品宾馆在线| 国产日韩欧美在线精品| 久久久午夜欧美精品| 国产成人福利小说| 高清av免费在线| 午夜免费鲁丝| 久久久久性生活片| 久久久精品免费免费高清| 国产大屁股一区二区在线视频| 三级男女做爰猛烈吃奶摸视频| 九九久久精品国产亚洲av麻豆| videossex国产| 免费黄色在线免费观看| av在线亚洲专区| av播播在线观看一区| 一级二级三级毛片免费看| 亚洲精品456在线播放app| 熟妇人妻不卡中文字幕| 久久久a久久爽久久v久久| 黄色欧美视频在线观看| 国产精品女同一区二区软件| 中文精品一卡2卡3卡4更新| 我的女老师完整版在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人人妻人人爽人人添夜夜欢视频 | 国产美女午夜福利| 亚洲国产精品成人综合色| 国产国拍精品亚洲av在线观看| av在线观看视频网站免费| 午夜福利在线观看免费完整高清在| 亚洲欧美成人精品一区二区| 伊人久久国产一区二区| 建设人人有责人人尽责人人享有的 | 亚洲欧美日韩卡通动漫| 国产亚洲5aaaaa淫片| 久久国产乱子免费精品| 国产国拍精品亚洲av在线观看| 国产精品秋霞免费鲁丝片| 男女无遮挡免费网站观看| 男女国产视频网站| 免费播放大片免费观看视频在线观看| 午夜日本视频在线| 在线观看av片永久免费下载| 91午夜精品亚洲一区二区三区| 国产精品福利在线免费观看| 22中文网久久字幕| av免费观看日本| 黄片wwwwww| 99热网站在线观看| 精品国产三级普通话版| 一级毛片我不卡| 国产亚洲av片在线观看秒播厂| 欧美日韩一区二区视频在线观看视频在线 | 日韩,欧美,国产一区二区三区| 2022亚洲国产成人精品| 一级二级三级毛片免费看| 久久99热6这里只有精品| 久久亚洲国产成人精品v| 免费高清在线观看视频在线观看| 精品久久久久久电影网| 国产淫片久久久久久久久| 久久精品熟女亚洲av麻豆精品| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| a级毛色黄片| 免费看光身美女| 我的老师免费观看完整版| 亚洲精品视频女| av在线天堂中文字幕| 成人二区视频| 欧美日韩一区二区视频在线观看视频在线 | 神马国产精品三级电影在线观看| 视频区图区小说| videossex国产| 国产 一区 欧美 日韩| 一个人看的www免费观看视频| 99热网站在线观看| 亚洲国产精品成人久久小说| 人人妻人人看人人澡| 深夜a级毛片| 麻豆国产97在线/欧美| 狂野欧美激情性bbbbbb| 国产又色又爽无遮挡免| 亚洲精品成人久久久久久| 男男h啪啪无遮挡| 嘟嘟电影网在线观看| 夫妻性生交免费视频一级片| 国产视频首页在线观看| 亚洲av在线观看美女高潮| 亚洲在久久综合| 日韩中字成人| 日本欧美国产在线视频| 亚洲国产精品999| 久久久久久九九精品二区国产| 国产男人的电影天堂91| 日韩强制内射视频| 女的被弄到高潮叫床怎么办| av福利片在线观看| 国产成人精品婷婷| 51国产日韩欧美| 国产探花在线观看一区二区| 只有这里有精品99| 18禁在线无遮挡免费观看视频| 欧美一级a爱片免费观看看| 国产永久视频网站| av免费在线看不卡| 少妇人妻精品综合一区二区| videos熟女内射| 国产亚洲av片在线观看秒播厂| 99热这里只有是精品在线观看| 久久精品久久精品一区二区三区| 听说在线观看完整版免费高清| 91久久精品国产一区二区三区| 日本猛色少妇xxxxx猛交久久| 精品国产三级普通话版| 国产 一区精品| 国产精品偷伦视频观看了| 欧美精品国产亚洲| 视频中文字幕在线观看| 亚洲欧美日韩另类电影网站 | 搡女人真爽免费视频火全软件| 校园人妻丝袜中文字幕| 99久久精品国产国产毛片| 97在线人人人人妻| 麻豆久久精品国产亚洲av| 亚洲av.av天堂| 插阴视频在线观看视频| 亚洲av成人精品一区久久| 1000部很黄的大片| 欧美日韩视频高清一区二区三区二| 国产亚洲av嫩草精品影院| 女的被弄到高潮叫床怎么办| 免费大片黄手机在线观看| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽人人片av| 精品少妇黑人巨大在线播放| 国产精品一及| 国产69精品久久久久777片| 亚洲丝袜综合中文字幕| 三级国产精品欧美在线观看| 一二三四中文在线观看免费高清| 联通29元200g的流量卡| 国产黄频视频在线观看| 免费观看无遮挡的男女| 国产成人午夜福利电影在线观看| 大码成人一级视频| 国产日韩欧美亚洲二区| 国产精品成人在线| 嫩草影院入口| 成人欧美大片| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 欧美日韩视频高清一区二区三区二| 天堂俺去俺来也www色官网| 国产精品一区二区在线观看99| 国产av不卡久久| 蜜桃久久精品国产亚洲av| 亚洲经典国产精华液单| 美女高潮的动态| 亚洲国产欧美人成| 久久99热这里只频精品6学生| 久久影院123| 国产精品熟女久久久久浪| 又黄又爽又刺激的免费视频.| 少妇的逼水好多| 国产伦在线观看视频一区| 亚洲国产精品999| 伊人久久精品亚洲午夜| 成人毛片a级毛片在线播放| 99热网站在线观看| 视频中文字幕在线观看| 成年免费大片在线观看| 色综合色国产| 日韩强制内射视频| 成人一区二区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 三级国产精品片| 日日摸夜夜添夜夜添av毛片| 日韩三级伦理在线观看| 狠狠精品人妻久久久久久综合| 中文字幕av成人在线电影| 丰满少妇做爰视频| 女人久久www免费人成看片| 啦啦啦在线观看免费高清www| 国产 一区 欧美 日韩| av国产免费在线观看| 亚洲欧美中文字幕日韩二区| 在线天堂最新版资源| 日本爱情动作片www.在线观看| 久久久精品免费免费高清| 国产成人免费观看mmmm| 国产69精品久久久久777片| 2022亚洲国产成人精品| 色婷婷久久久亚洲欧美| 全区人妻精品视频| 国产成人a∨麻豆精品| 精品酒店卫生间| 久久ye,这里只有精品| 丝袜喷水一区| 3wmmmm亚洲av在线观看| 韩国高清视频一区二区三区| 精品熟女少妇av免费看| 成人特级av手机在线观看| 又黄又爽又刺激的免费视频.| 亚洲精品自拍成人| 午夜免费鲁丝| 国产又色又爽无遮挡免| 亚洲天堂av无毛| 亚洲欧美一区二区三区国产| 永久网站在线| 欧美日韩视频高清一区二区三区二| 国产在线男女| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 午夜福利在线在线| 插阴视频在线观看视频| 黄色日韩在线| 久久6这里有精品| 久久久欧美国产精品| 91aial.com中文字幕在线观看| 成年人午夜在线观看视频| 国产男女超爽视频在线观看| 有码 亚洲区| 国产精品久久久久久精品电影小说 | 97人妻精品一区二区三区麻豆| 国产成人aa在线观看| kizo精华| 97人妻精品一区二区三区麻豆| 最近的中文字幕免费完整| 大陆偷拍与自拍| 午夜免费鲁丝| 人妻 亚洲 视频| 大码成人一级视频| 大香蕉97超碰在线| 日日摸夜夜添夜夜添av毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产在视频线精品| 亚洲经典国产精华液单| av国产精品久久久久影院| 成人综合一区亚洲| 七月丁香在线播放| 三级国产精品片| 白带黄色成豆腐渣| 久久6这里有精品| 一级av片app| 国产女主播在线喷水免费视频网站| 亚洲欧美日韩卡通动漫| 丝袜喷水一区| 一级毛片 在线播放| av线在线观看网站| av天堂中文字幕网| 国国产精品蜜臀av免费| 国产视频内射| 亚洲成色77777| 岛国毛片在线播放| 色综合色国产| 神马国产精品三级电影在线观看| 国产成人免费无遮挡视频| 久久99热这里只频精品6学生| 国产乱来视频区| 内地一区二区视频在线| 国产伦在线观看视频一区| 全区人妻精品视频| 各种免费的搞黄视频| 麻豆精品久久久久久蜜桃| 色网站视频免费| 日韩三级伦理在线观看| 啦啦啦啦在线视频资源| 欧美变态另类bdsm刘玥| 波野结衣二区三区在线| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区国产| 亚洲av成人精品一区久久| 亚洲精品视频女| 女人十人毛片免费观看3o分钟| 国产爽快片一区二区三区| 国产精品不卡视频一区二区| 国产精品福利在线免费观看| 69人妻影院| 日本一本二区三区精品| 免费看光身美女| 一级毛片 在线播放| 亚洲精品日韩av片在线观看| 97人妻精品一区二区三区麻豆| 狂野欧美激情性xxxx在线观看| 精品99又大又爽又粗少妇毛片| 国产精品国产三级专区第一集| 亚洲电影在线观看av| 国产男女内射视频| 精品一区二区三区视频在线| 亚洲天堂国产精品一区在线| 国产精品一区二区性色av| 国产综合精华液| 久久久久精品性色| 久久精品国产a三级三级三级| 99精国产麻豆久久婷婷| 久久影院123| 最近最新中文字幕大全电影3| 国产极品天堂在线| 在线观看免费高清a一片| 赤兔流量卡办理| 精品国产乱码久久久久久小说| av在线播放精品| 狂野欧美白嫩少妇大欣赏| 白带黄色成豆腐渣| 亚洲精品国产色婷婷电影| 中文字幕免费在线视频6| 国产女主播在线喷水免费视频网站| 超碰av人人做人人爽久久| 国产一区有黄有色的免费视频| 欧美日韩在线观看h| 久久久色成人| 联通29元200g的流量卡| 日产精品乱码卡一卡2卡三| 国产高清国产精品国产三级 | 欧美xxxx性猛交bbbb| 波野结衣二区三区在线| 高清av免费在线| 婷婷色av中文字幕| 久久久久国产网址| 国产老妇女一区| 精品久久久久久电影网| 免费大片18禁| 久久久色成人| 91狼人影院| 久久精品久久久久久噜噜老黄| 日本一二三区视频观看| 亚洲在久久综合| 一区二区三区精品91| 日本欧美国产在线视频| 成人高潮视频无遮挡免费网站| 欧美另类一区| 超碰97精品在线观看| 午夜福利高清视频| 日韩欧美精品免费久久| 99re6热这里在线精品视频| 国产av码专区亚洲av| 亚洲av一区综合| 麻豆乱淫一区二区| 国产亚洲av片在线观看秒播厂| 自拍偷自拍亚洲精品老妇| 中国三级夫妇交换| 狂野欧美白嫩少妇大欣赏| 色网站视频免费| 日韩不卡一区二区三区视频在线| 搡老乐熟女国产| 直男gayav资源| 免费av不卡在线播放| 久久久成人免费电影| 丝袜美腿在线中文| 青春草视频在线免费观看| 搡女人真爽免费视频火全软件| 51国产日韩欧美| 精品国产一区二区三区久久久樱花 | 欧美国产精品一级二级三级 | 久热这里只有精品99| 男女无遮挡免费网站观看| 欧美精品人与动牲交sv欧美| 爱豆传媒免费全集在线观看| 免费看日本二区| 久久久久久久久久久免费av| 欧美另类一区| 精品人妻一区二区三区麻豆| av天堂中文字幕网| 欧美日韩一区二区视频在线观看视频在线 | 久久午夜福利片| 可以在线观看毛片的网站| av.在线天堂| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| 久久久成人免费电影| 亚洲精华国产精华液的使用体验| 日本欧美国产在线视频| 久久久欧美国产精品| 成人亚洲欧美一区二区av| 欧美一区二区亚洲| 男女无遮挡免费网站观看| 久久久久精品久久久久真实原创| 中文乱码字字幕精品一区二区三区| 人人妻人人看人人澡| 丰满人妻一区二区三区视频av| 狂野欧美白嫩少妇大欣赏| 网址你懂的国产日韩在线| 精品久久久久久久久亚洲| 99热全是精品| 男的添女的下面高潮视频| 亚洲四区av| 熟女电影av网| 国产视频首页在线观看| 久久精品国产a三级三级三级| 两个人的视频大全免费| 欧美xxxx黑人xx丫x性爽| 在线精品无人区一区二区三 | 18+在线观看网站| 久久人人爽人人爽人人片va| 国产又色又爽无遮挡免| 日日摸夜夜添夜夜爱| 精品久久久久久久人妻蜜臀av| 国产精品三级大全| 新久久久久国产一级毛片| 男的添女的下面高潮视频| 免费观看的影片在线观看| 十八禁网站网址无遮挡 | 亚洲精品乱久久久久久| 亚洲久久久久久中文字幕| 久久久成人免费电影| 亚洲人成网站高清观看| 国产亚洲91精品色在线| 亚洲人成网站在线观看播放| 中国国产av一级| 欧美激情久久久久久爽电影| av国产久精品久网站免费入址| 肉色欧美久久久久久久蜜桃 | 成年女人看的毛片在线观看| 亚洲国产精品成人久久小说| 2022亚洲国产成人精品| 欧美成人精品欧美一级黄| 亚洲伊人久久精品综合| 在线观看一区二区三区激情| 欧美另类一区| 十八禁网站网址无遮挡 | 亚洲成人av在线免费| 狂野欧美激情性xxxx在线观看| 秋霞伦理黄片| 亚洲内射少妇av| 夜夜爽夜夜爽视频| 亚洲av免费在线观看| 直男gayav资源| 亚洲国产av新网站| 内地一区二区视频在线| 草草在线视频免费看| 亚洲三级黄色毛片| 日本av手机在线免费观看| 你懂的网址亚洲精品在线观看| 欧美成人一区二区免费高清观看| 国产高潮美女av| 18禁动态无遮挡网站| 97超视频在线观看视频| av在线亚洲专区| 一本久久精品| 国产午夜精品久久久久久一区二区三区| 99久久人妻综合| 有码 亚洲区| 日韩国内少妇激情av| 热99国产精品久久久久久7| 亚洲精品乱码久久久久久按摩| 韩国av在线不卡| 欧美日韩视频精品一区| 又黄又爽又刺激的免费视频.| 搡老乐熟女国产| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区国产| 免费电影在线观看免费观看| 看黄色毛片网站| 嫩草影院精品99| 国产精品一区二区在线观看99|