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

    基于VMD和SVDD的滾動軸承早期微弱故障檢測和性能退化評估研究

    2019-12-02 06:05:48房立清趙玉龍齊子元
    振動與沖擊 2019年22期
    關鍵詞:壽命軸承模態(tài)

    王 斐, 房立清, 趙玉龍, 齊子元

    (1.武警士官學校 軍械系,杭州 310023; 2.陸軍工程大學 石家莊校區(qū)火炮工程系,石家莊 050003)

    滾動軸承是大型機械系統(tǒng)的重要零部件,一旦軸承失效將會導致整個機械系統(tǒng)的災難性故障,造成巨大的經(jīng)濟損失和嚴重事故,目前針對滾動軸承的狀態(tài)監(jiān)測和維修策略一直是研究的熱點?;跔顟B(tài)的維修CBM(Condition Based Maintenance),或稱視情維修,是對機械系統(tǒng)的健康狀態(tài)進行實時監(jiān)測,管理和維護,提高整個系統(tǒng)的運行可靠性[1]。軸承從正常運轉到完全失效往往要經(jīng)歷一系列不同的性能退化過程,如果能夠在軸承運行過程中監(jiān)測到早期微弱故障及軸承性能退化的規(guī)律,跟蹤其故障的發(fā)展,那么就可以使設備維護計劃的制定具有針對性,因此滾動軸承的故障退化程度的評估是開展視情維修的基礎和重點[2]。故障退化程度評估要求狀態(tài)特征要對性能退化程度要有量化反應能力,進而實現(xiàn)對全壽命周期中任意時刻的退化程度進行定量的度量,而非單純描述不同故障之間的差異;這是故障退化程度評估區(qū)別于傳統(tǒng)故障診斷的方面。

    故障退化程度評估包含兩個方面:一是性能退化特征的提取;二是退化特征指標的評估。針對這兩方面,國內外學者進行了大量的研究,文獻[3]提出了一種基于模糊C均值聚類(Fuzzy C-means Clustering,FCM)的設備性能退化評估方法,以正常狀態(tài)和失效時刻的數(shù)據(jù)為基礎,建立智能評估模型,在對滾動軸承性能退化評估中取得了良好的效果;文獻[4]提出了基于流行學習特征提取和高斯混合模型的退化評估方法,效果優(yōu)于時域特征提取方法;文獻[5]提出了基于隱馬爾可夫的多通道融合性能退化評估模型,通過滾動軸承加速壽命數(shù)據(jù)驗證了方法的有效性;文獻[6]提出了結合EMD(Empirical Mode Decomposition)分解特征提取和K均值聚類的故障退化評估方法,可以較好的發(fā)現(xiàn)早期故障退化階段;這些方法很好的實現(xiàn)了特征指標的提取與故障程度的評估,但對于早期微弱故障的檢測效果并不是很明顯。

    變分模態(tài)分解(Variational Mode Decomposition,VMD)是文獻[7]提出的一種非遞歸的自適應信號分解方法,VMD通過求解變分模態(tài)最優(yōu)解實現(xiàn)模態(tài)分解,可將一個非平穩(wěn)振動信號自適應地分解為若干個不同中心頻率的本征模態(tài)函數(shù)(Intrinsic Mode Function, IMF)之和,具有精度高、收斂速度快、對噪聲有較好的魯棒性等特點,與EMD和EEMD方法相比,有效克服了端點效應與模態(tài)混疊現(xiàn)象。目前基于VMD方法的機械故障診斷領域取得了很好的效果[8-13]。

    支持向量數(shù)據(jù)描述(Support Vector Data Description, SVDD)[14]是一種有效的單值(One-class)分類方法,其核心思想是通過非線性映射將樣本映射到高維的內積空間,在特征空間中尋找一個包含全部或大部分樣本且體積最小的超球體(最優(yōu)超球體)。近年來,SVDD方法在故障檢測、臉部識別、語音識別、圖像恢復等領域取得了良好的效果。利用SVDD的單值分類特性,區(qū)分正常狀態(tài)和其他退化狀態(tài),可有效實現(xiàn)機械零件性能退化程度的評估[15-18]。

    基于上述描述,本文提出了結合VMD分解和SVDD的滾動軸承性能退化評估方法,利用VMD方法對振動信號進行分解,選取對故障較為敏感的模態(tài)分量,提取分量的奇異值并結合時域特征、復雜度特征組成滾動軸承綜合特征指標;利用SVDD在單值分類問題上的優(yōu)勢,僅使用正常樣本的綜合特征指標進行SVDD評估模型的構建,通過滾動軸承全壽命數(shù)據(jù)進行模型的驗證,驗證結果表明本文所提出的方法不僅能有效的描述滾動軸承在全壽命周期內的性能退化規(guī)律,同是對軸承早期微弱故障的發(fā)生具有良好的預警能力。

    1 變分模態(tài)分解原理介紹

    VMD算法的核心思想是在變分框架中對信號分解,每個模態(tài)分量的頻率中心及帶寬在迭代求解變分模型的過程中不斷更新,最后將信號自適應分解為K個IMF函數(shù)。分解過程的數(shù)學模型為

    (1)

    式中:f為輸入信號;uk為分解得到的K個IMF分量;ωk為各IMF對應的頻率中心。

    (2)

    (3)

    VMD方法相比EMD方法,在含噪脈沖信號分解,故障特征頻率提取,克服模態(tài)混疊方面有很大的優(yōu)越性。因此,本文將其應用于滾動軸承信號的分解,并進一步提取故障特征作為軸承性能退化評估的特征指標。

    2 支持向量數(shù)據(jù)描述

    SVDD可以有效解決單值分類問題,其核心思想是將給定的訓練數(shù)據(jù)映射到高維空間中,在高維空間中尋找一個體積最小的超球體,使得訓練數(shù)據(jù)盡可能多地包含于在該超球體中,不屬于該類的數(shù)據(jù)則位于超球體之外。其數(shù)學模型構建如下:

    給定一個目標類樣本集X={x1,x2,…,xn},定義一個包含所有或幾乎所有目標類樣本的超球體,超球體球心為a,半徑為R??紤]到目標樣本中可能存在少數(shù)野點,為了提高算法對訓練樣本中野點的魯棒性,引入松弛因子ξi,允許部分樣本點分布于超球體之外,此時SVDD優(yōu)化問題描述為

    (4)

    其中,C>0為指定的懲罰參數(shù),用以調節(jié)錯分訓練樣本數(shù)(球外樣本數(shù))和r的大小。為了解決上述優(yōu)化問題,構造以下拉格朗日方程

    式中:αi,γi為拉格朗日乘子。對R,a,ξi求偏導數(shù),令其為0,得到

    (6)

    將式(6)代入式(2)中得到優(yōu)化函數(shù)

    (7)

    當樣本點位于超球體內是,αi=0當樣本點位于超球體的邊界上,則0<αi

    R2=‖φ(xsv)-φ(a)‖2=K(xsv,xsv)-

    (8)

    對于新樣本z,它與球心的距離可表示為

    d=‖φ(z)-φ(a)‖=K(z,z)-

    (9)

    核函數(shù)K的使用完成由原空間向高維空間的轉換,把原來線性不可分的問題轉換成高維空間中線性可分的問題,同是并沒有增加太多的計算量。

    在軸承退化性能評估過程中,將提取到的正常狀態(tài)下的樣本特征向量作為訓練樣本輸入SVDD進行訓練,得到包含正常狀態(tài)樣本的超球體,球心為O,半徑為R;新樣本與超球體球心的廣義距離d可由式(9)求出,當d≤R時,則樣本為正常樣本;當d>R樣本屬于故障樣本,故障樣本d的取值越大則故障程度越深。

    3 滾動軸承性能退化程度評估

    3.1 滾動軸承全壽命實驗數(shù)據(jù)

    本文所使用的滾動軸承全壽命數(shù)據(jù)來自于辛辛那提智能維護中心軸承疲勞壽命試驗臺,如圖1所示。軸承振動信號每隔10 min采集一次,采樣頻率20 kHz,采樣時間1 s,試驗臺運行163 h后軸承2外圈出現(xiàn)嚴重剝落故障,壽命終止,共采集到984組樣本,此數(shù)據(jù)即為軸承的全壽命周期試驗數(shù)據(jù),經(jīng)計算軸承2外圈故障頻率約為234 Hz,文中采用每組樣本的第一列數(shù)據(jù)進行性能退化評估研究。

    圖1 軸承疲勞壽命試驗臺Fig.1 Rolling bearing fatigue life test platform

    3.2 特征指標的選擇

    滾動軸承全壽命周期一般分為正常狀態(tài),輕微退化階段,嚴重退化階段和失效階段,一般壽命評估研究的重點為輕微退化階段和嚴重退化階段。良好的性能退化特征指標一是要在全壽命周期的各個退化階段要保持相對穩(wěn)定,另一方面需要在初始退化階段有較強的敏感性。目前,有多種特征指標被應用于性能退化指標的提取,如時域特征指標(均方根值,峰—峰值,裕度,峭度等),復雜度指標(關聯(lián)維數(shù)、近似熵、樣本熵、模糊熵、譜熵等)。本文經(jīng)過大量分析研究表明,VMD奇異值,均方根值,樣本熵值可以更好的反映滾動軸承全壽命周期指標的變化趨勢。

    3.2.1 VMD奇異值指標

    在進行VMD分解之前,需要確定分解層數(shù)K的取值,以避免出現(xiàn)過分解可能帶來的模態(tài)混疊現(xiàn)象。分解層數(shù)的選取通過比較各IMF分量的中心頻率的方法確定,VMD 分解得到的K個 IMF 分量的中心頻率是從小到大依次排列的,隨著K的增大,各個分量的中心頻率之間的差值會減小,最終出現(xiàn)過分解現(xiàn)象,通過比較各IMF中心頻率的差值,來判斷是否過分解,進而便可以確定模態(tài)個數(shù)K。這里以1 000 min時的振動信號為例進行分析,表1為不同模態(tài)個數(shù)下的中心頻率的值,從中可知,當K=7時,分解得到的IMF4和IMF5分量的中心頻率相差較小,因此,當K≥7時,信號有過分解的風險,而當模態(tài)個數(shù)過小時,各模態(tài)分量之間中心頻率的差值過大,又不足以表征原信號所包含的不同頻率成分;考慮到降低特征提取維數(shù)的目的,因此最終選取K=5作為VMD分解的模態(tài)個數(shù)。其它時刻模態(tài)個數(shù)的選取方法相同,這里不做過多說明。

    表1 不同模態(tài)個數(shù)下IMF分量的中心頻率Tab.1 Center frequency of the IMF under different K

    將滾動軸承振動信號進行VMD分解,得到5個IMF分量,結果如圖2所示,進一步計算每個分量的奇異值,奇異值可依據(jù)下面的奇異值分解公式進行求解

    A=UQVT

    (10)

    式中:A∈Rm×n為待分解的矩陣;U∈Rm×m和V∈Rn×n為正交矩陣;Q∈Rm×n為對角矩陣。

    (11)

    式中:0為零矩陣;q=min(m,n);且α1≥α2≥…≥αq>0,αi(i=1,2,…,q)為矩陣A的奇異值。由于IMF分量為一維時域信號,因此每個IMF分量經(jīng)奇異值分解后得到一個奇異值。經(jīng)計算得到5個分量在全壽命周期內奇異值的變化情況,如圖3所示;圖3反映出這5個IMF分量的奇異值在全壽命周期內變化趨勢基本一致,但IMF3~IMF5分量的奇異值在故障早期階段相比IMF1和IMF2變化更為敏感,其包含了更多故障信息,更能表征全壽命特征指標的變化,因此選取IMF3~IMF5奇異值作為部分故障特征向量。

    圖2 軸承振動信號VMD分解Fig.2 Bearing vibration signal VMD decomposition

    圖3 IMF信號奇異值在全壽命周期內的變化值Fig.3 The variation of singular value of IMF signal in the whole life period

    3.2.2均方根值指標與樣本熵值指標

    圖4 均方根值特征Fig.4 RMS feature

    樣本熵是常用復雜度特征指標之一,是近似熵的一種改進形式,它通過衡量時間序列復雜性度量信號中產(chǎn)生新模式的概率大小。熵值越大時間序列復雜度越大。樣本熵克服了數(shù)據(jù)偏差,具有較強的抗噪能力和優(yōu)異的一致性,采用較少的數(shù)據(jù)段即可得到穩(wěn)定的熵值。對樣本熵的具體計算過程可參考文獻[19]。滾動軸承全壽命周期內樣本熵值的變化情況如圖5所示。

    分析奇異值,均方根值和樣本熵指標的變化情況,奇異值和均方根值指標在各個不同階段數(shù)值較為穩(wěn)定,對輕微退化階段D1和嚴重退化階段D2都具有較好的表征能力,樣本熵指標雖然在數(shù)值穩(wěn)定性上不如均方根值,但在D1,D2階段當中值的變化更大一些,對故障狀態(tài)的變化更為敏感。因此,經(jīng)綜合考慮將這三種指標結合起來作為滾動軸承性能退化綜合特征指標,記為V=[SV3,SV4,SV5,RMS,SampEn]。

    圖5 樣本熵特征Fig.5 Sample entropy feature

    3.3 滾動軸承性能退化評估模型

    通過上述分析,本文提出了基于VMD分解和SVDD的滾動軸承性能退化評估模型,如圖6所示。其具體步驟為:

    步驟1特征提取。對信號進行5層VMD分解,提取IMF3~IMF5分量的奇異值,與信號均方根值和樣本熵值組成特征向量V作為軸承綜合特征指標。

    步驟2利用正常樣本信號的綜合特征指標作為訓練樣本構建SVDD分類器,得到包含正常樣本超球面O,并得到超球體半徑R。

    步驟3計算測試樣本綜合特征指標Vtest,輸入SVDD分類器中求得Vtest與超球體球心之間的距離d,并計算得到性能退化指標(Degradation Indicator,DI),以DI值作為性能退化評估的依據(jù)

    (12)

    圖6 性能退化評估模型Fig.6 Performance degradation assessment model

    3.4 基于VMD和SVDD的滾動軸承性能退化評估

    從圖3~圖5可知,滾動軸承在0~4 000 min運行期間內,各特征值保持穩(wěn)定,可認為期間軸承處于正常工作狀態(tài),本文選取運行時間在1 000~3 000 min的200組振動信號作為樣本數(shù)據(jù)用于SVDD的訓練樣本。其中,懲罰參數(shù)C取值為1, 核函數(shù)參數(shù)σ=90。將滾動軸承全壽命周期的984組數(shù)據(jù)作為測試樣本輸入退化評估模型最終得到全壽命周期的DI值,以此作為全壽命周期的評估結果,如圖7(a)所示,圖7(b)為500~700組樣本的DI變化情況,圖7(c)為只采用樣本熵與均方根值作為特征指標的評估結果,與文中所提特征提取方法相比,監(jiān)測到初始故障發(fā)生的時間晚了約260 min,這說明本文提出的SV3,SV4,SV5特征指標能很好的反映軸承運行狀態(tài)信息。

    圖7 基于SVDD的滾動軸承性能退化評估結果Fig.7 Performance degradation assessment results by SVDD

    圖7(a)和圖7(b)中,在運行的0~5 340 min期間樣本的DI=0,此階段的滾動軸承處于正常工作狀態(tài),在5 340 min樣本的DI>0,說明此時開始出現(xiàn)極其輕微的早期故障,從5 550 min開始故障開始逐漸加深,直到7 050 minDI值急劇變化,說明滾動軸承性能退化加深,故障程度較為嚴重,之后DI值時而增大時而減小,規(guī)律性較弱,從9 100 min開始DI值大幅增大,軸承失效。從對評估結果的分析可以得到以下推論:①5 340~5 540 min,軸承出現(xiàn)極其微弱的早期故障,幾乎難以發(fā)現(xiàn);②在5 500~7 040 min,故障程度呈逐漸加深態(tài)勢;③7 050~9 840 min,DI變化情況復雜,無明顯規(guī)律,可以認為從7 050 min以后軸承已經(jīng)處于嚴重故障狀態(tài)。

    為了驗證本文的推論,分別對1 000 min,5 330 min,5 340 min,5 400 min,5 550 min,7 050 min的振動信號作VMD分解,模態(tài)個數(shù)K=5,取IMF3~IMF5分量重構信號,并對重構信號做包絡譜分析,如圖8所示。5 340 min時的包絡譜分析可以明顯檢測到234.4 Hz的頻率成分,非常接近外圈故障特征頻率的理論計算結果234 Hz,5 340 min時刻之前的5 330 min及1 000 min時的正常狀態(tài)則沒有此頻率成分,可以認為5 340 min第一次發(fā)生了外圈故障;此時的包絡譜還存在很多其他頻率成分,說明此時的外圈故障非常微弱,5 340 min之后的5 440 min還能較為明顯的觀察到234.4 Hz及其倍頻成分459 Hz,但頻率的幅值較5 340 min時更小,這是因為出現(xiàn)微弱故障后故障點被“磨平”使得故障程度減弱,相比5 340 min,5 340~5 550 min的故障更弱,更加難以檢測,可以認為5 340~5 550 min滾動軸承處于早期微弱故障緩慢退化階段;從5 550 min開始,外圈故障特征頻率及其倍頻成分愈發(fā)明顯,直到7 050 min時包絡譜中幾乎只存在故障特征頻率及其倍頻成分。由此可認為5 550~7 050 min軸承處于性能加速退化階段;7 050 min以后外圈故障更為明顯,本文認為7 050 min以后軸承其實已經(jīng)處于性能嚴重退化狀態(tài)。與文獻[20]的評估結果相比,本文方法所檢測到的初始故障退化階段提前了約1 690 min,相比文獻[21]提前了約690 min,這充分說明本文所提出的退化評估模型對滾動軸承性能退化有較強的敏感性,而盡早檢測到滾動軸承故障情況的發(fā)生對機械系統(tǒng)的視情維修又具有重要意義。

    圖8 不同時間的振動信號包絡譜分析Fig.8 Envelope spectral analysis results under different times

    為了進一步說明文中所提方法的優(yōu)越性,同樣利用模糊聚類分析算法進行滾動軸承性能退化評估,使用正常樣本和故障樣本訓練FCM評估模型,依然選取V為綜合特征指標,將全壽命的特征指標作為測試樣本輸入評估模型,將測試樣本與正常樣本特征的隸屬度值作為性能退化指標DI,DI∈[0,1],進而得到滾動軸承全壽命周期的評估結果,結果如圖9(a)所示,圖9(b)為局部放大情況。由圖9可知,F(xiàn)CM評估模型的DI值的趨勢與圖7變化趨勢相似,在0~5 000 min期間內DI值波動較大,且FCM模型在6 100 min左右才出現(xiàn)較為明顯的DI值減小趨勢,相比本文所提方法晚了約550 min,本文所提方法對數(shù)據(jù)異常點具有較強的魯棒性,并且在早期微弱故障監(jiān)測方面明顯優(yōu)于FCM評估模型。同時,本文所提方法只使用了正常狀態(tài)下的樣本作為訓練樣本,一方面可以減少計算量,簡化評估模型;另一方面可以克服實際工業(yè)生產(chǎn)中故障樣本較難獲取的問題,對于進行機械系統(tǒng)性能監(jiān)測具有重要意義。

    圖9 基于FCM的滾動軸承全壽命性能退化評估結果Fig.9 FCM-based evaluation of the life-cycle performance degradation of rolling bearings

    4 結 論

    本文針對滾動軸承微弱故障特征較難提取問題及性能退化程度監(jiān)測問題,提出了一種基于VMD分解和SVDD的滾動軸承性能退化評估模型,利用軸承全壽命數(shù)據(jù)對所提出的評估模型進行了驗證,主要結論如下:

    (1)對振動信號進行VMD分解,通過提取IMF分量的奇異值與信號均方根值、樣本熵值組成綜合特征指標,可以較好的表征軸承性能退化情況。

    (2)以正常狀態(tài)的樣本作為訓練樣本進行評估模型的構建,所得到性能退化評估結果可以準確檢測到軸承早期微弱故障的發(fā)生,并可以對軸承性能退化情況進行很好的評估,對實際生產(chǎn)中軸承微弱故障的發(fā)生和性能退化評估具有很好的指導意義。

    猜你喜歡
    壽命軸承模態(tài)
    軸承知識
    哈爾濱軸承(2022年2期)2022-07-22 06:39:32
    軸承知識
    哈爾濱軸承(2022年1期)2022-05-23 13:13:24
    人類壽命極限應在120~150歲之間
    中老年保健(2021年8期)2021-12-02 23:55:49
    軸承知識
    哈爾濱軸承(2021年2期)2021-08-12 06:11:46
    軸承知識
    哈爾濱軸承(2021年1期)2021-07-21 05:43:16
    倉鼠的壽命知多少
    馬烈光養(yǎng)生之悟 自靜其心延壽命
    華人時刊(2018年17期)2018-12-07 01:02:20
    人類正常壽命為175歲
    奧秘(2017年12期)2017-07-04 11:37:14
    國內多模態(tài)教學研究回顧與展望
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    欧美老熟妇乱子伦牲交| 嫩草影视91久久| 啦啦啦在线免费观看视频4| 19禁男女啪啪无遮挡网站| 乱人伦中国视频| 亚洲色图av天堂| 国产成人精品久久二区二区免费| 很黄的视频免费| 黄色视频,在线免费观看| 亚洲国产看品久久| 亚洲专区中文字幕在线| 侵犯人妻中文字幕一二三四区| 精品国内亚洲2022精品成人| 亚洲成国产人片在线观看| 色精品久久人妻99蜜桃| 国产亚洲精品久久久久5区| 欧美午夜高清在线| 丰满饥渴人妻一区二区三| 日日摸夜夜添夜夜添小说| 狂野欧美激情性xxxx| 中亚洲国语对白在线视频| 免费看十八禁软件| 日韩一卡2卡3卡4卡2021年| 可以免费在线观看a视频的电影网站| 高清在线国产一区| 国产熟女xx| 亚洲色图av天堂| av在线播放免费不卡| 国产精品一区二区三区四区久久 | 亚洲一区二区三区欧美精品| 桃红色精品国产亚洲av| 午夜免费激情av| 18禁黄网站禁片午夜丰满| 校园春色视频在线观看| 在线观看免费高清a一片| 欧美性长视频在线观看| 日本免费一区二区三区高清不卡 | 每晚都被弄得嗷嗷叫到高潮| 日本欧美视频一区| 久久中文看片网| 午夜福利免费观看在线| av网站免费在线观看视频| 欧美成狂野欧美在线观看| 十分钟在线观看高清视频www| 久久伊人香网站| 窝窝影院91人妻| 岛国在线观看网站| 亚洲自拍偷在线| 少妇被粗大的猛进出69影院| 久久久国产精品麻豆| 一二三四社区在线视频社区8| 国产有黄有色有爽视频| 麻豆国产av国片精品| 久久精品国产99精品国产亚洲性色 | 淫妇啪啪啪对白视频| 国产免费男女视频| 午夜精品国产一区二区电影| 久久狼人影院| 黄色视频,在线免费观看| 一区二区三区精品91| 精品久久久久久电影网| 制服人妻中文乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久久免费视频了| 热99国产精品久久久久久7| 亚洲一区二区三区色噜噜 | 大型黄色视频在线免费观看| 亚洲国产毛片av蜜桃av| av欧美777| 亚洲中文av在线| 少妇裸体淫交视频免费看高清 | 性少妇av在线| 香蕉丝袜av| 亚洲自偷自拍图片 自拍| 成人国语在线视频| 亚洲成av片中文字幕在线观看| 精品一区二区三区av网在线观看| 免费在线观看黄色视频的| 亚洲中文av在线| 国产欧美日韩一区二区三区在线| 91大片在线观看| 国产麻豆69| 国产欧美日韩精品亚洲av| 国产亚洲av高清不卡| 97超级碰碰碰精品色视频在线观看| 久久久久久大精品| 波多野结衣av一区二区av| 变态另类成人亚洲欧美熟女 | 亚洲色图 男人天堂 中文字幕| 欧美精品一区二区免费开放| 欧美日韩福利视频一区二区| 满18在线观看网站| 51午夜福利影视在线观看| 久久人妻福利社区极品人妻图片| 极品教师在线免费播放| 成在线人永久免费视频| 色老头精品视频在线观看| 日韩国内少妇激情av| 黑人欧美特级aaaaaa片| www.www免费av| av有码第一页| 精品电影一区二区在线| 日本黄色视频三级网站网址| 一二三四在线观看免费中文在| 真人做人爱边吃奶动态| www.精华液| 狂野欧美激情性xxxx| 99riav亚洲国产免费| 国产免费现黄频在线看| av在线播放免费不卡| 亚洲国产看品久久| 午夜影院日韩av| av中文乱码字幕在线| 精品国产国语对白av| 亚洲精品国产色婷婷电影| 人妻久久中文字幕网| 黄片大片在线免费观看| 18美女黄网站色大片免费观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲狠狠婷婷综合久久图片| 一区福利在线观看| 午夜a级毛片| 巨乳人妻的诱惑在线观看| 国产精品综合久久久久久久免费 | 亚洲激情在线av| 麻豆久久精品国产亚洲av | 国产亚洲av高清不卡| 亚洲熟女毛片儿| 国产视频一区二区在线看| 精品国内亚洲2022精品成人| 国产麻豆69| 国产熟女午夜一区二区三区| 69av精品久久久久久| 久久久国产精品麻豆| 美女大奶头视频| 国产精品综合久久久久久久免费 | 黑丝袜美女国产一区| 亚洲视频免费观看视频| 99精品欧美一区二区三区四区| 久久国产精品男人的天堂亚洲| 一区福利在线观看| 在线免费观看的www视频| 国产男靠女视频免费网站| 国产精品影院久久| 激情视频va一区二区三区| 免费在线观看视频国产中文字幕亚洲| 欧美最黄视频在线播放免费 | 夜夜爽天天搞| 大陆偷拍与自拍| 18禁国产床啪视频网站| 精品人妻在线不人妻| 99精品欧美一区二区三区四区| 在线天堂中文资源库| 国产一区在线观看成人免费| 精品人妻1区二区| 日韩精品青青久久久久久| 久久久久久免费高清国产稀缺| 精品久久久精品久久久| 免费一级毛片在线播放高清视频 | 亚洲欧美精品综合一区二区三区| 国产精品免费视频内射| 精品少妇一区二区三区视频日本电影| 又黄又粗又硬又大视频| 搡老熟女国产l中国老女人| 色综合婷婷激情| 成人影院久久| 精品久久久久久久久久免费视频 | 久久午夜亚洲精品久久| 成人黄色视频免费在线看| a级毛片在线看网站| 中文字幕人妻丝袜一区二区| 亚洲视频免费观看视频| 久久精品人人爽人人爽视色| 国产免费现黄频在线看| 99国产综合亚洲精品| 19禁男女啪啪无遮挡网站| 中文字幕人妻丝袜一区二区| 免费高清在线观看日韩| 亚洲精品国产精品久久久不卡| 美女高潮到喷水免费观看| 99在线人妻在线中文字幕| 亚洲视频免费观看视频| 黄色女人牲交| 欧美在线黄色| 免费av毛片视频| 女人高潮潮喷娇喘18禁视频| 丰满饥渴人妻一区二区三| 一区二区三区激情视频| 淫妇啪啪啪对白视频| 18禁黄网站禁片午夜丰满| 十分钟在线观看高清视频www| 日本撒尿小便嘘嘘汇集6| 波多野结衣一区麻豆| 桃红色精品国产亚洲av| 欧美成人免费av一区二区三区| 黑人欧美特级aaaaaa片| 亚洲一区中文字幕在线| 新久久久久国产一级毛片| 国产精品一区二区三区四区久久 | 国产在线精品亚洲第一网站| 国产精品永久免费网站| 国产片内射在线| 国产麻豆69| 夜夜夜夜夜久久久久| 十分钟在线观看高清视频www| 午夜成年电影在线免费观看| 天堂影院成人在线观看| 久久国产精品人妻蜜桃| 一级毛片高清免费大全| 久久精品亚洲精品国产色婷小说| 高潮久久久久久久久久久不卡| 亚洲国产欧美网| 水蜜桃什么品种好| 国产欧美日韩一区二区三区在线| 国产高清视频在线播放一区| 国产av在哪里看| 无遮挡黄片免费观看| 国产精品香港三级国产av潘金莲| 国产三级黄色录像| 久久精品人人爽人人爽视色| 黄频高清免费视频| 国产成人影院久久av| 国产色视频综合| 亚洲七黄色美女视频| 在线天堂中文资源库| 午夜福利在线免费观看网站| 最近最新免费中文字幕在线| 高清黄色对白视频在线免费看| 久久精品亚洲av国产电影网| 女警被强在线播放| 国产亚洲欧美98| 夜夜夜夜夜久久久久| www国产在线视频色| 亚洲人成伊人成综合网2020| 免费久久久久久久精品成人欧美视频| 精品一区二区三区四区五区乱码| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 1024视频免费在线观看| 俄罗斯特黄特色一大片| 亚洲aⅴ乱码一区二区在线播放 | 老司机午夜福利在线观看视频| 99久久国产精品久久久| 精品欧美一区二区三区在线| 国产在线精品亚洲第一网站| 久久久久久久久免费视频了| 欧美久久黑人一区二区| 免费在线观看日本一区| 可以在线观看毛片的网站| 搡老熟女国产l中国老女人| 欧美日韩精品网址| 国产亚洲精品久久久久久毛片| 欧美日韩瑟瑟在线播放| 麻豆成人av在线观看| 久久久久久免费高清国产稀缺| 黑人巨大精品欧美一区二区mp4| 色哟哟哟哟哟哟| 亚洲少妇的诱惑av| 一级a爱片免费观看的视频| 天天影视国产精品| 另类亚洲欧美激情| 一进一出好大好爽视频| 男女做爰动态图高潮gif福利片 | 中文字幕av电影在线播放| 99国产精品99久久久久| 91精品国产国语对白视频| 日韩精品中文字幕看吧| 国产av在哪里看| 久久久久久久午夜电影 | 亚洲人成电影观看| 精品免费久久久久久久清纯| 一a级毛片在线观看| 天天影视国产精品| 黄频高清免费视频| 夜夜爽天天搞| 午夜精品在线福利| xxx96com| 欧美 亚洲 国产 日韩一| 久久国产精品影院| 身体一侧抽搐| 久久精品aⅴ一区二区三区四区| 成在线人永久免费视频| 久久性视频一级片| 成人亚洲精品av一区二区 | 欧美不卡视频在线免费观看 | 国产欧美日韩一区二区三区在线| 一区福利在线观看| 国产伦人伦偷精品视频| 国产欧美日韩综合在线一区二区| svipshipincom国产片| 看黄色毛片网站| 777久久人妻少妇嫩草av网站| 亚洲片人在线观看| 亚洲国产精品一区二区三区在线| x7x7x7水蜜桃| 久久久久久久久中文| 悠悠久久av| 久久精品国产清高在天天线| 亚洲黑人精品在线| 日本vs欧美在线观看视频| 黄色怎么调成土黄色| 女同久久另类99精品国产91| 80岁老熟妇乱子伦牲交| 免费女性裸体啪啪无遮挡网站| 美女 人体艺术 gogo| 精品卡一卡二卡四卡免费| 一本综合久久免费| 国产欧美日韩一区二区三区在线| 亚洲成av片中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放 | 69精品国产乱码久久久| 老司机午夜福利在线观看视频| 中文字幕高清在线视频| 超碰成人久久| 亚洲精品在线观看二区| 国产午夜精品久久久久久| 一区在线观看完整版| 久久热在线av| 九色亚洲精品在线播放| 美女午夜性视频免费| 高清黄色对白视频在线免费看| 99久久精品国产亚洲精品| 久久影院123| 国产精品98久久久久久宅男小说| 国产又爽黄色视频| 国产有黄有色有爽视频| 三级毛片av免费| 我的亚洲天堂| 亚洲精品中文字幕在线视频| 99在线人妻在线中文字幕| 成年人黄色毛片网站| 欧洲精品卡2卡3卡4卡5卡区| 久久伊人香网站| 两性午夜刺激爽爽歪歪视频在线观看 | 黄频高清免费视频| 免费在线观看亚洲国产| 欧美最黄视频在线播放免费 | 亚洲一区二区三区色噜噜 | 免费女性裸体啪啪无遮挡网站| 18禁裸乳无遮挡免费网站照片 | 男人的好看免费观看在线视频 | 两人在一起打扑克的视频| 美女福利国产在线| 欧美一区二区精品小视频在线| 国产精品二区激情视频| а√天堂www在线а√下载| 性少妇av在线| 国产亚洲精品综合一区在线观看 | 神马国产精品三级电影在线观看 | 国产一区二区三区在线臀色熟女 | 国产伦一二天堂av在线观看| 岛国在线观看网站| 日韩欧美免费精品| 免费高清在线观看日韩| 女性被躁到高潮视频| 一本大道久久a久久精品| 两人在一起打扑克的视频| 久久天躁狠狠躁夜夜2o2o| 19禁男女啪啪无遮挡网站| 淫秽高清视频在线观看| 国产免费现黄频在线看| 中国美女看黄片| 黄片大片在线免费观看| 欧美精品一区二区免费开放| 精品欧美一区二区三区在线| 欧美久久黑人一区二区| 啦啦啦在线免费观看视频4| 亚洲av片天天在线观看| 久久午夜综合久久蜜桃| 50天的宝宝边吃奶边哭怎么回事| av电影中文网址| 国产三级在线视频| 一a级毛片在线观看| a级片在线免费高清观看视频| 欧美老熟妇乱子伦牲交| 日日干狠狠操夜夜爽| 欧美精品一区二区免费开放| 国产麻豆69| 婷婷丁香在线五月| 好看av亚洲va欧美ⅴa在| 在线观看舔阴道视频| 亚洲成人免费电影在线观看| 新久久久久国产一级毛片| 99re在线观看精品视频| 男女床上黄色一级片免费看| 亚洲久久久国产精品| 久久天堂一区二区三区四区| 亚洲专区中文字幕在线| 一边摸一边抽搐一进一出视频| 热99国产精品久久久久久7| 国产精品野战在线观看 | 亚洲全国av大片| 亚洲成人免费电影在线观看| 丁香欧美五月| 99热只有精品国产| 亚洲七黄色美女视频| 国产av又大| 亚洲成国产人片在线观看| 极品人妻少妇av视频| 国产成人av教育| 久久人妻av系列| 亚洲精品一卡2卡三卡4卡5卡| 韩国精品一区二区三区| 搡老岳熟女国产| 桃色一区二区三区在线观看| 国产野战对白在线观看| x7x7x7水蜜桃| 日韩精品青青久久久久久| 免费在线观看黄色视频的| 国产一区二区三区视频了| 欧美色视频一区免费| 身体一侧抽搐| av福利片在线| 纯流量卡能插随身wifi吗| 欧美人与性动交α欧美精品济南到| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 国产精品1区2区在线观看.| 久久99一区二区三区| 成人永久免费在线观看视频| 国产伦人伦偷精品视频| 香蕉久久夜色| 91av网站免费观看| a级毛片黄视频| 免费人成视频x8x8入口观看| 又黄又爽又免费观看的视频| 最近最新中文字幕大全免费视频| 曰老女人黄片| 久久亚洲真实| 国产精品1区2区在线观看.| av超薄肉色丝袜交足视频| 精品福利观看| 电影成人av| 两性午夜刺激爽爽歪歪视频在线观看 | 视频区图区小说| 侵犯人妻中文字幕一二三四区| 亚洲国产精品一区二区三区在线| 九色亚洲精品在线播放| 久久亚洲真实| 久久精品91蜜桃| 亚洲五月婷婷丁香| 999精品在线视频| 欧美乱妇无乱码| 欧美亚洲日本最大视频资源| 99在线视频只有这里精品首页| 悠悠久久av| 99在线视频只有这里精品首页| 大型黄色视频在线免费观看| 精品一区二区三区四区五区乱码| 村上凉子中文字幕在线| 国产成人欧美在线观看| 男女下面进入的视频免费午夜 | 亚洲黑人精品在线| 一二三四在线观看免费中文在| av超薄肉色丝袜交足视频| 亚洲成人免费电影在线观看| 国产熟女午夜一区二区三区| ponron亚洲| 一级片免费观看大全| 欧美乱码精品一区二区三区| 免费观看人在逋| 成人永久免费在线观看视频| cao死你这个sao货| 国产视频一区二区在线看| 国产精品av久久久久免费| 久久影院123| 99国产精品免费福利视频| 黄色怎么调成土黄色| 男人舔女人下体高潮全视频| 日本vs欧美在线观看视频| 亚洲第一av免费看| av天堂久久9| 岛国视频午夜一区免费看| 日日爽夜夜爽网站| 亚洲狠狠婷婷综合久久图片| 免费在线观看视频国产中文字幕亚洲| 欧美成人性av电影在线观看| 少妇 在线观看| 嫁个100分男人电影在线观看| 国产一区二区三区视频了| 99久久国产精品久久久| 国产三级在线视频| 久久伊人香网站| 亚洲视频免费观看视频| 免费在线观看日本一区| 视频区欧美日本亚洲| 亚洲一区二区三区不卡视频| 国产亚洲欧美在线一区二区| 好看av亚洲va欧美ⅴa在| 亚洲一码二码三码区别大吗| 亚洲中文字幕日韩| 午夜精品久久久久久毛片777| 日韩精品青青久久久久久| 久久精品国产99精品国产亚洲性色 | 亚洲国产欧美日韩在线播放| 国产麻豆69| 在线观看免费视频网站a站| 欧美久久黑人一区二区| 国产av精品麻豆| 国产精品免费视频内射| 欧美日韩视频精品一区| 男女午夜视频在线观看| 国产精品久久电影中文字幕| 亚洲精品美女久久久久99蜜臀| 脱女人内裤的视频| 免费观看人在逋| 一边摸一边做爽爽视频免费| 99久久国产精品久久久| 亚洲精品久久成人aⅴ小说| 久久青草综合色| 日韩视频一区二区在线观看| 97人妻天天添夜夜摸| 超色免费av| 97碰自拍视频| 亚洲欧美精品综合久久99| 在线观看一区二区三区| 午夜免费鲁丝| 黑丝袜美女国产一区| 午夜激情av网站| 两性夫妻黄色片| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 亚洲熟妇中文字幕五十中出 | 亚洲精华国产精华精| 欧美日韩黄片免| 夜夜夜夜夜久久久久| 精品一品国产午夜福利视频| 亚洲精品久久午夜乱码| 俄罗斯特黄特色一大片| 国产精品99久久99久久久不卡| 欧美色视频一区免费| 精品久久久久久,| 伦理电影免费视频| 国产一区二区在线av高清观看| 高清欧美精品videossex| 丝袜人妻中文字幕| 老汉色∧v一级毛片| 久久久久久亚洲精品国产蜜桃av| 久久狼人影院| 久久 成人 亚洲| 亚洲aⅴ乱码一区二区在线播放 | 精品无人区乱码1区二区| 日本欧美视频一区| 亚洲色图综合在线观看| 少妇裸体淫交视频免费看高清 | 天堂中文最新版在线下载| 少妇 在线观看| 这个男人来自地球电影免费观看| 精品久久久久久成人av| 免费高清视频大片| 不卡av一区二区三区| 国产色视频综合| 美女大奶头视频| 久久热在线av| 亚洲九九香蕉| 中文亚洲av片在线观看爽| 99在线视频只有这里精品首页| 午夜影院日韩av| 人人澡人人妻人| 成在线人永久免费视频| 麻豆成人av在线观看| 亚洲成a人片在线一区二区| 久久精品国产清高在天天线| 国产亚洲精品综合一区在线观看 | 不卡一级毛片| 亚洲一区高清亚洲精品| 宅男免费午夜| 国产熟女xx| 嫁个100分男人电影在线观看| 91字幕亚洲| 久久久久国产一级毛片高清牌| 亚洲成人国产一区在线观看| 又黄又粗又硬又大视频| 两性夫妻黄色片| 丁香欧美五月| 亚洲精品av麻豆狂野| 精品福利永久在线观看| 亚洲国产欧美日韩在线播放| 欧洲精品卡2卡3卡4卡5卡区| 19禁男女啪啪无遮挡网站| 国产精品1区2区在线观看.| aaaaa片日本免费| 丰满的人妻完整版| 搡老熟女国产l中国老女人| 黄片大片在线免费观看| 欧美av亚洲av综合av国产av| 日日夜夜操网爽| 日韩av在线大香蕉| 黄色 视频免费看| 成人亚洲精品av一区二区 | 老司机亚洲免费影院| 真人一进一出gif抽搐免费| 亚洲伊人色综图| 男女高潮啪啪啪动态图| 黑丝袜美女国产一区| 久久精品91蜜桃| 亚洲avbb在线观看| 视频区欧美日本亚洲| 一级a爱视频在线免费观看| 女人精品久久久久毛片| 亚洲男人天堂网一区| 欧美人与性动交α欧美软件| 欧美日韩一级在线毛片| 久久精品国产亚洲av香蕉五月| 精品国产国语对白av| 在线观看午夜福利视频| 女人高潮潮喷娇喘18禁视频| 人人妻人人添人人爽欧美一区卜| 最近最新中文字幕大全免费视频| 午夜精品国产一区二区电影| 日本三级黄在线观看| 90打野战视频偷拍视频| 欧美成人午夜精品| 亚洲欧美激情综合另类| 91精品三级在线观看|