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

    基于熵值比-DTW度量奇異值相似度指標(biāo)的滾動軸承退化歷程辨識

    2023-01-09 09:48:24周建清朱文昌王恒
    軸承 2023年1期
    關(guān)鍵詞:歷程軸承模板

    周建清,朱文昌,王恒

    (1.常州市高級職業(yè)技術(shù)學(xué)校 電氣工程學(xué)院,江蘇 常州 213161; 2.南通大學(xué) 機械工程學(xué)院,江蘇 南通 226019)

    有效監(jiān)測滾動軸承性能退化歷程,制定針對性的維護(hù)計劃,可避免設(shè)備損壞造成的財產(chǎn)損失及人員傷亡,具有重要意義[1]。基于數(shù)據(jù)分析的滾動軸承退化歷程辨識得到了廣泛研究,其關(guān)鍵在于如何構(gòu)造單調(diào)性好且能準(zhǔn)確刻畫不同狀態(tài)的性能退化指標(biāo)。傳統(tǒng)方法采用軸承振動信號的時域(如峭度、均方根)、頻域(如重心頻率、均方頻率)及時頻域(如小波包分解、經(jīng)驗?zāi)B(tài)分解)等特征進(jìn)行分析[2]。文獻(xiàn)[3]將 Kullback-Leibler散度作為軸承健康退化指標(biāo),可有效量化滾動軸承的不同退化階段;文獻(xiàn)[4]改進(jìn)了卷積神經(jīng)網(wǎng)絡(luò)并用于深層次振動特征提取,結(jié)合遷移學(xué)習(xí)實現(xiàn)了不同工況下軸承運行狀態(tài)的識別;文獻(xiàn)[5]使用經(jīng)驗?zāi)B(tài)分解算法提取軸承故障特征,對特征向量進(jìn)行 k-medoids 聚類以構(gòu)建模型,可更精確地檢測早期退化;文獻(xiàn)[6]利用應(yīng)用累積并對JRD(Jensen Renyi Divergence)進(jìn)行改進(jìn),提高了退化指標(biāo)的穩(wěn)定性及單調(diào)性;文獻(xiàn)[7]將排列熵(Permutation Entropy,PE)作為軸承退化特征,通過改進(jìn)構(gòu)造出熵能比(Entropy Energy Rate, EER)并作為滾動軸承性能退化指標(biāo),其對早期故障更敏感且與軸承故障發(fā)展趨勢更加一致。

    隨著人工智能、機器學(xué)習(xí)的發(fā)展,以卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Network, CNN)為代表的深度學(xué)習(xí)算法在軸承特征提取方面也得到了廣泛的應(yīng)用[8],與傳統(tǒng)方法相比,其能夠挖掘軸承數(shù)據(jù)更深層次的信息,但由于池化層、隱含層的層數(shù)等超參數(shù)需要人為設(shè)定,缺乏一定的科學(xué)性。近年來,許多學(xué)者將矩陣論應(yīng)用于滾動軸承狀態(tài)監(jiān)測領(lǐng)域,文獻(xiàn)[9]提出一種非負(fù)矩陣分解,用于軸承特征提取并實現(xiàn)了軸承故障診斷;文獻(xiàn)[10]基于隨機矩陣?yán)碚搶S承特征矩陣進(jìn)行分解,將矩陣特征值用于軸承退化指標(biāo)構(gòu)建與檢測閾值設(shè)定,實現(xiàn)了滾動軸承的早期異常檢測。

    目前,常利用滾動軸承某個時刻退化指標(biāo)的數(shù)值大小及變化趨勢監(jiān)測軸承當(dāng)前時刻的狀態(tài),然而軸承性能退化是一個連續(xù)變化的過程,即當(dāng)前狀態(tài)與歷史狀態(tài)相關(guān),如何將軸承當(dāng)前狀態(tài)及歷史信息相結(jié)合,構(gòu)建能準(zhǔn)確反映軸承當(dāng)前狀態(tài)與正常狀態(tài)之間差異的性能指標(biāo)并提高指標(biāo)對不同階段的敏感性,避免數(shù)據(jù)波動干擾軸承退化狀態(tài)的判定是值得深入研究的問題。動態(tài)時間規(guī)整(Dynamic Time Warping,DTW)是一種計算時間序列相似性的算法,將時間規(guī)整和距離測度計算相結(jié)合,利用對2個時序信號的拉伸、對齊等操作獲得時間校準(zhǔn)匹配路徑,通過計算路徑間的最短距離描述時間序列間的相似度與差異性[11]。文獻(xiàn)[12]利用DTW算法計算時間序列間的相似性并進(jìn)行狀態(tài)匹配,與歐幾里得距離相比,其對時間序列相似度的計算更為準(zhǔn)確;文獻(xiàn)[13]利用DTW算法分析手語運動軌跡間的相似性并判斷其是否屬于同一類別,提高了手語軌跡識別的準(zhǔn)確率。傳統(tǒng)DTW算法假設(shè)序列中的各個點有著相同的作用,從而計算序列間的相似度,而在軸承的實際退化過程中,不同時間段采集到的數(shù)據(jù)結(jié)構(gòu)及組成可能存在較大的變化,因此利用DTW算法計算軸承信號相似度時不能僅考慮數(shù)據(jù)序列之間的距離,也應(yīng)考慮不同狀態(tài)數(shù)據(jù)間的差異性。

    綜上所述,本文基于奇異值分解(Singular Value Decomposition,SVD)算法對軸承數(shù)據(jù)進(jìn)行特征提取,結(jié)合信息熵和DTW算法計算不同狀態(tài)軸承退化數(shù)據(jù)間的相似度,并用其表征軸承全壽命過程中的不同退化歷程。

    1 基于奇異值分解的軸承特征提取

    根據(jù)滾動軸承健康監(jiān)測數(shù)據(jù)的采集特點,在時間T內(nèi)對軸承進(jìn)行N次采樣,采樣時間記為ti(i=1,2,…,N),采樣點數(shù)為M,各數(shù)據(jù)點分別記為x(ti,j)(j=1,2,…,M),在ti時刻采集的數(shù)據(jù)可構(gòu)成向量x(ti),即

    x(ti)=[x(ti,1),x(ti,2),…,x(ti,M)] 。

    (1)

    向量x(ti)中涵蓋ti時刻軸承不同位置的運行信息,可以看作矩陣x(ti)∈R1×M。對于一個實矩陣B∈RK×K,可引入特征值λ及特征向量y將其替代,λ可用于B中的有用信息,則

    By=λy。

    (2)

    將采集的數(shù)據(jù)構(gòu)建矩陣并進(jìn)行特征分解,將特征值構(gòu)造退化指標(biāo)以識別軸承的不同退化階段[14]。在矩陣論中,矩陣奇異值分解與特征值分解的作用較為相似,不同的是奇異值可適用于矩陣不滿秩的情況,具有更大的應(yīng)用范圍。而在矩陣論中,對于一個實矩陣A∈RH×J,可對其進(jìn)行奇異值分解,則A可寫為

    (3)

    Λ=diag(σ1,σ2,…,σl),

    式中:V與VT為正交特征向量矩陣;Λ為對角矩陣;σi為采樣時刻ti對軸承數(shù)據(jù)進(jìn)行奇異值分解后所對應(yīng)的奇異值,且l=min(H,J)。

    對ti時刻采集到的軸承數(shù)據(jù)x(ti)∈R1×M進(jìn)行分解后,可從對角矩陣Λ中提取到一個奇異值,故在采樣時間T內(nèi),對不同時刻的軸承數(shù)據(jù)進(jìn)行分解并提取不同時刻的奇異值,可構(gòu)造出退化時間序列G∈R1×N,即

    G=[σ1,σ2,σ3,…,σi,…,σN]。

    (4)

    2 基于熵值比-動態(tài)時間規(guī)整構(gòu)造奇異值相似度指標(biāo)

    滾動軸承由正常運行至完全失效是連續(xù)變化的過程,當(dāng)軸承受到損傷進(jìn)入退化狀態(tài),其損傷的衍變規(guī)律應(yīng)與前n個歷史狀態(tài)有關(guān),單獨取某一時刻的特征去判斷軸承狀態(tài)是不合理的,因此本文將相互鄰近的多個特征點構(gòu)成時間序列,利用DTW算法分析軸承的退化歷程。

    確定一段軸承正常狀態(tài)信號序列作為參考模版P=[p1,p2,…,pa],構(gòu)造其他時間段的退化序列作為測試模板Q=[q1,q2,…,qb],其中pa和qb分別表示參考模板的第a個和測試模板的第b個特征矢量值。對齊路徑li用于描述P與Q之間數(shù)據(jù)點的對齊關(guān)系,其被定義為包含s個二元組的集合,每個二元組包含2個分別來自時間序列P和Q的數(shù)據(jù)點,l可表示為

    (5)

    序列P與Q之間所有對齊路徑的集合記為AP,Q。DTW的目標(biāo)是最小化兩段時序數(shù)據(jù)中所有對應(yīng)數(shù)據(jù)點的局部距離值之和,其定義為

    (6)

    d(pa,qb)=∣pa-qb∣2,

    (7)

    式中:d(pa,qb)為矢量間的距離。

    DTW算法采用動態(tài)規(guī)劃思想,利用遞歸公式將以上問題轉(zhuǎn)換為對P和Q中特征矢量距離的求解問題,即

    (8)

    將規(guī)整后路徑間的距離作為量化指標(biāo),可對不同序列的相似性進(jìn)行有效度量,規(guī)整后的距離越小則序列間相似度越大。在軸承的實際退化過程中,不同狀態(tài)對應(yīng)的數(shù)據(jù)結(jié)構(gòu)及組成可能產(chǎn)生較大的變化,所蘊含的信息不同,在計算序列間相似性時應(yīng)考慮不同狀態(tài)數(shù)據(jù)結(jié)構(gòu)間的差異性。熵值可反映信息可靠程度,系統(tǒng)混亂程度越高則熵值越大,本文將信息熵之比引入DTW算法中用于優(yōu)化奇異值指標(biāo),其步驟如下:

    1)從退化時間序列G∈R1×N選取a個連續(xù)的正常數(shù)據(jù)點作為參考模板P并保持不變。采用滑動時間窗口從時間序列G中提取不同時刻的數(shù)據(jù)作為測試模板,設(shè)置時間窗長度為b,則ti時刻構(gòu)造的測試模板Q(ti)為

    Q(ti)=[σi-b,…,σi-2,σi-1,σi] ;b≥i≥N,

    (9)

    從σb開始,滑動時間窗鎖定的第1個測試模板為Q(t1)=[σ1,σ2,σ3,…,σb],隨后時間窗每次后移1個單位,共移動N-b+1次。利用DTW算法計算2個模板間的相似度,在ti時刻測試模板Q(ti)與參考模板P的距離(即相似度)記為di。

    2)分別計算參考模板和測試模版中奇異值對應(yīng)相似信號的信息熵值及熵值比,即

    (10)

    (11)

    式中:a,b分別為參考模板和測試模板中相似信號對應(yīng)奇異值的個數(shù);xk(ti,j)為將ti時刻所采集數(shù)據(jù)x(ti,j)分解k次后相似信號中第j個數(shù)據(jù)的值;wi為ti時刻測試模板與參考模板距離間的權(quán)值。

    3)將ti時刻對應(yīng)的權(quán)值wi與DTW相似度距離di點乘得到加權(quán)后的相似度dw,i。將不同時間段的相似度通過權(quán)熵值進(jìn)行優(yōu)化可得到加權(quán)后的距離矩陣D=(dw,1,dw,2,…,dw,i,…,dw,N-b+1),將其歸一化處理后用于對軸承退化歷程進(jìn)行辨識。

    綜上所述,本文所提基于熵值比-DTW度量奇異值相似度指標(biāo)的構(gòu)造流程如圖1所示。

    圖1 基于熵值比-DTW度量奇異值相似度指標(biāo)構(gòu)建流程圖Fig.1 Constructed flowchart of singular value similarity index based on DTW optimized by entropy ratio

    3 應(yīng)用研究

    3.1 數(shù)據(jù)來源

    本文采用辛辛那提大學(xué)智能維護(hù)系統(tǒng)(IMS)中心的滾動軸承全壽命試驗中軸承1的數(shù)據(jù)進(jìn)行應(yīng)用研究,該試驗采用加速度傳感器每 10 min采集一次軸承振動信號,采樣頻率為20 kHz,當(dāng)軸承1外圈發(fā)生故障時結(jié)束試驗,整個試驗共經(jīng)歷9 830 min。取奇異值序列G中正常狀態(tài)下的序列段[σ1,σ2,σ3,σ4,σ5]為參考模板,測試模板長度b=5,利用DTW計算其相似度并結(jié)合熵值比進(jìn)行優(yōu)化,采用優(yōu)化后的相似度指標(biāo)辨識軸承1的退化歷程。

    3.2 基于SVD算法的軸承奇異值特征提取

    利用SVD算法對軸承各個時刻的數(shù)據(jù)進(jìn)行分解,提取各個時刻對應(yīng)的奇異值并進(jìn)行歸一化處理后,構(gòu)建的軸承全壽命歷程退化指標(biāo)如圖2所示,通過3σ準(zhǔn)則可在5 370 min檢測出軸承早期異常的發(fā)生[15],與文獻(xiàn)[16-17]分別為5 350,5 330 min的檢測結(jié)果比較接近,與峭度指標(biāo)(圖3)相比可提前1 122 min檢測出軸承早期異常的發(fā)生,且奇異值指標(biāo)在退化時期整體呈單調(diào)向上的趨勢,數(shù)據(jù)曲線波動小,穩(wěn)定性高。然而,在軸承“自愈現(xiàn)象”發(fā)生(7 000~9 000 min)時,奇異值指標(biāo)存在明顯“上升—下降—再上升”的起伏波動[18],給軸承狀態(tài)判斷帶來了干擾,難以判斷圖3中A,B范圍內(nèi)具有相同指標(biāo)值數(shù)據(jù)點所對應(yīng)的退化狀態(tài)。

    圖2 基于奇異值指標(biāo)軸承全壽命歷程Fig.2 Bearing full life history based on singular value index

    圖3 基于峭度指標(biāo)的軸承全壽命歷程Fig.3 Bearing full life history based on kurtosis index

    3.3 基于熵值比-DTW構(gòu)造奇異值相似度指標(biāo)

    為進(jìn)一步提高退化指標(biāo)整體的單調(diào)性及對早期異常的敏感性,將動態(tài)時間規(guī)整算法與熵值法相結(jié)合,從時間序列及數(shù)據(jù)結(jié)構(gòu)兩方面共同對軸承奇異值指標(biāo)進(jìn)行處理。以10~50 min軸承正常狀態(tài)奇異值序列為參考模版,提取1 210~1 250 min與8 010~8 050 min奇異值序列作為測試模板分別計算相似度,計算過程如圖4、圖5所示。在參考模板不變時,由圖4a可發(fā)現(xiàn)測試模板為1 210~1 250 min時,對于2個原始序列變化趨勢不同的曲線,可通過DTW算法對曲線進(jìn)行伸長、對齊操作后將其轉(zhuǎn)變?yōu)榫哂邢嗤兓厔莸那€,計算對齊后各對應(yīng)點間的最小距離并最終得到規(guī)整后的最短路徑(從數(shù)值較小的深色區(qū)經(jīng)過),如圖4b所示,累積距離大小為0.867,這是由于軸承在1 210~1 250 min期間處于正常狀態(tài),參考模板也處于正常狀態(tài),2個序列的相似度較高,規(guī)整得到累積距離較??;當(dāng)測試模板為8 010~8 050 min 時,序列間的差異較大,DTW算法難以將其規(guī)整為相同退化趨勢的曲線,且規(guī)劃出的路徑經(jīng)過淺色區(qū)域,累計距離較大(11.415,約正常狀態(tài)下相似度的13倍),如圖5所示。因此, DTW算法可放大不同時間序列間的差異性。

    (a) 時間序列規(guī)整

    (a) 時間序列規(guī)整

    利用DTW處理軸承數(shù)據(jù)時只考慮了不同時間段時間序列的相似性,未考慮不同時間段數(shù)據(jù)結(jié)構(gòu)之間的差異性。因此,本文采用熵值比分析軸承的全壽命歷程,結(jié)果如圖6所示:當(dāng)軸承處于正常狀態(tài),熵值比近似為一條直線;當(dāng)軸承進(jìn)入異常狀態(tài),熵值比有明顯的上升趨勢;隨著軸承故障加劇,熵值比不斷增大,表明熵值比表征數(shù)據(jù)結(jié)構(gòu)間的差異是可行的。

    圖6 滾動軸承全壽命歷程的熵值比變化曲線Fig.6 Change curve based on entropy ratio of rolling bearing's full life history

    以(10)、(11)式計算測試模板中奇異值所對應(yīng)相似信號的熵值比作為權(quán)值,用其優(yōu)化動態(tài)時間規(guī)整算法所測的相似度,最終構(gòu)造出軸承1的性能退化曲線,結(jié)果如圖7所示:利用3σ準(zhǔn)則檢測出軸承1的早期異常點在5 330 min,優(yōu)化后的相似度指標(biāo)可更早檢測出軸承的早期異常點,且指標(biāo)的整體單調(diào)性及對異常狀態(tài)的敏感性均得到了較大的提高。

    圖7 基于奇異值相似度指標(biāo)的軸承退化曲線Fig.7 Bearing degradation curve based on singular value similarity index

    3.4 軸承退化歷程辨識

    基于奇異值相似度指標(biāo)對軸承的退化歷程辨識結(jié)果如圖8所示:

    1)在5 330 min前指標(biāo)變化平緩,波動較小,可劃分為正常階段。

    2)5 330~6 910 min期間退化指標(biāo)較正常狀態(tài)有較為明顯的上升趨勢,可定義為軸承的早期退化階段。

    3)7 020~8 640 min期間指標(biāo)幅值比上一階段更大,但上升并不明顯,可定義為軸承的中期退化階段。這是由于軸承1發(fā)生了“自愈現(xiàn)象”,即當(dāng)前階段軸承表面由于長時間運行出現(xiàn)了裂紋及小的剝落,滾動體不斷運轉(zhuǎn)將此缺陷撫平使軸承仍可保持較穩(wěn)定的運行狀態(tài),但軸承的疲勞磨損仍在發(fā)生,指標(biāo)仍有上升趨勢。與處于相同階段的指標(biāo)相比(圖2、圖3),基于熵值比-DTW構(gòu)建的奇異值相似度指標(biāo)的單調(diào)性較好,避免了數(shù)據(jù)波動對狀態(tài)判斷帶來的干擾。

    4)軸承的嚴(yán)重退化發(fā)生在8 640~9 560 min期間,在此之前經(jīng)歷了“損傷—愈合—再損傷”的階段,在8 640 min后,軸承損傷加劇,運行狀態(tài)極不穩(wěn)定,退化加速;直到9 560 min后,退化曲線上升劇烈,軸承失效并將完全損壞。

    經(jīng)過熵值比-DTW構(gòu)造的奇異值相似度指標(biāo)檢測的早期異常點、嚴(yán)重故障的檢測結(jié)果與文獻(xiàn)[19-20]基本一致,證明了基于相似度退化曲線對軸承退化歷程劃分的有效性。

    圖8 基于奇異值相似度指標(biāo)的軸承全壽命狀態(tài)識別Fig.8 Identification of bearing full life state based on singular value similarity index

    3.5 不同退化指標(biāo)性能對比

    為進(jìn)一步研究奇異值相似度指標(biāo)的有效性,從2個方面對指標(biāo)進(jìn)行量化分析:1)單調(diào)性,滾動軸承退化具有不可逆性,構(gòu)建的退化指標(biāo)是否隨軸承運行時間共同增長;2)魯棒性,構(gòu)造出的退化指標(biāo)是否具有抵抗數(shù)據(jù)頻繁波動干擾的能力。單調(diào)性、魯棒性分別定義為

    m(D)=

    (12)

    (13)

    提取軸承1 在5 330 min(早期異常點)后的退化序列進(jìn)行單調(diào)性及魯棒性分析,并與基于隨機矩陣?yán)碚?RMT)和主成分分析(PCA)結(jié)合構(gòu)造的融合特征指標(biāo)[14],基于優(yōu)化經(jīng)驗小波變換及卷積神經(jīng)網(wǎng)絡(luò)的EWT-CNN指標(biāo)[19],傳統(tǒng)峭度和均方根(RMS)指標(biāo)進(jìn)行對比,結(jié)果見表1。

    表1 不同退化指標(biāo)的單調(diào)性及魯棒性Tab.1 Monotonicity and robustness of different degradation indexes

    由表1可知:與其他退化指標(biāo)相比,基于熵值比-DTW算法的奇異值相似度指標(biāo)的單調(diào)性及魯棒性均有明顯提升;這是由于DTW算法通過計算軸承不同時間序列的相似度放大了時間段的差異,結(jié)合信息熵比進(jìn)一步凸顯了正常與異常狀態(tài)數(shù)據(jù)結(jié)構(gòu)內(nèi)部的差異性,顯著抑制了軸承全壽命歷程中退化指標(biāo)的“上升—下降”波動現(xiàn)象,可更準(zhǔn)確地表征軸承全壽命歷程中的損傷衍變過程。另外,奇異值相似度指標(biāo)劃分確定的軸承早期異常點及嚴(yán)重故障點與EWT-CNN指標(biāo)的劃分結(jié)果較為接近,可更好地辨識軸承的不同退化階段。

    4 結(jié)束語

    本文借鑒矩陣特征值在軸承健康監(jiān)測領(lǐng)域中的良好應(yīng)用效果,對軸承信號進(jìn)行奇異值分解并構(gòu)造奇異值退化序列,同時考慮到軸承退化的連續(xù)性及所采集軸承數(shù)據(jù)間的差異性,將DTW及熵值比優(yōu)化后構(gòu)建出的奇異值相似度指標(biāo)用于表征軸承異常的衍變過程,為軸承健康狀態(tài)監(jiān)測的研究提供了一種新思路。

    對退化指標(biāo)單調(diào)性及魯棒性的量化分析表明本文奇異值相似度指標(biāo)可有效克服軸承退化過程中由于數(shù)據(jù)頻繁波動對不同退化階段識別帶來的干擾,且與其他退化指標(biāo)相比具有更好的敏感性,整體性能較好;但本文所提優(yōu)化算法的整體流程較復(fù)雜,需進(jìn)一步研究算法復(fù)雜度的降低問題。

    猜你喜歡
    歷程軸承模板
    鋁模板在高層建筑施工中的應(yīng)用
    鋁模板在高層建筑施工中的應(yīng)用
    百年大黨壯闊歷程
    黨員文摘(2022年14期)2022-08-06 07:19:22
    軸承知識
    哈爾濱軸承(2022年2期)2022-07-22 06:39:32
    軸承知識
    哈爾濱軸承(2022年1期)2022-05-23 13:13:24
    百年大黨 壯闊歷程
    黨員文摘(2022年7期)2022-04-28 08:57:58
    百年大黨 壯闊歷程
    黨員文摘(2022年1期)2022-02-14 06:09:00
    百年大黨壯闊歷程 《百色起義》
    黨員文摘(2022年3期)2022-02-12 10:38:14
    軸承知識
    哈爾濱軸承(2021年2期)2021-08-12 06:11:46
    軸承知識
    哈爾濱軸承(2021年1期)2021-07-21 05:43:16
    一二三四社区在线视频社区8| 最近视频中文字幕2019在线8| 国产精品影院久久| 国产免费男女视频| 久久精品国产99精品国产亚洲性色| 我要搜黄色片| 最好的美女福利视频网| 亚洲国产欧美网| 女人被狂操c到高潮| 一进一出好大好爽视频| 欧美大码av| 色吧在线观看| 精品久久久久久久久久久久久| 国产91精品成人一区二区三区| 国内久久婷婷六月综合欲色啪| 欧美午夜高清在线| 最近最新中文字幕大全电影3| 国产精品久久久人人做人人爽| 男女那种视频在线观看| 99热只有精品国产| 伊人久久大香线蕉亚洲五| 亚洲一区二区三区不卡视频| 欧美国产日韩亚洲一区| 欧美黑人巨大hd| 午夜福利在线观看吧| 天堂动漫精品| 日本黄色视频三级网站网址| 国产探花在线观看一区二区| 国产69精品久久久久777片| 观看美女的网站| 国产精品免费一区二区三区在线| 大型黄色视频在线免费观看| 午夜日韩欧美国产| 国内久久婷婷六月综合欲色啪| 一本久久中文字幕| 久久久精品欧美日韩精品| 男女下面进入的视频免费午夜| 午夜福利在线观看免费完整高清在 | 99国产综合亚洲精品| av欧美777| 久久精品91蜜桃| 真人做人爱边吃奶动态| 日本撒尿小便嘘嘘汇集6| 中亚洲国语对白在线视频| 九色成人免费人妻av| 亚洲国产精品999在线| 999久久久精品免费观看国产| 特级一级黄色大片| 亚洲在线观看片| 性欧美人与动物交配| 免费观看的影片在线观看| 90打野战视频偷拍视频| 久久久成人免费电影| 最好的美女福利视频网| 欧美三级亚洲精品| 少妇的逼水好多| 欧美色欧美亚洲另类二区| 草草在线视频免费看| 国产激情欧美一区二区| 最近最新中文字幕大全免费视频| 精品国产超薄肉色丝袜足j| 国产中年淑女户外野战色| 午夜影院日韩av| 亚洲av电影不卡..在线观看| 91麻豆精品激情在线观看国产| 日韩欧美在线乱码| 免费搜索国产男女视频| 又粗又爽又猛毛片免费看| 蜜桃亚洲精品一区二区三区| 波多野结衣高清作品| 国产不卡一卡二| 成人国产一区最新在线观看| 国产精品1区2区在线观看.| 久久精品国产综合久久久| 精品久久久久久久久久免费视频| 美女高潮的动态| 日本一本二区三区精品| 午夜精品一区二区三区免费看| 亚洲欧美日韩东京热| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产精品合色在线| 国产探花在线观看一区二区| 亚洲国产精品久久男人天堂| 一级黄片播放器| 男插女下体视频免费在线播放| 女人十人毛片免费观看3o分钟| 国产亚洲欧美在线一区二区| 午夜影院日韩av| 日韩人妻高清精品专区| 国产精品一区二区三区四区免费观看 | 久久香蕉国产精品| 黄片小视频在线播放| 99久久精品热视频| 亚洲国产精品999在线| 国产精品影院久久| 伊人久久大香线蕉亚洲五| 国产免费一级a男人的天堂| 免费av不卡在线播放| 国产免费男女视频| 亚洲国产欧美网| 一区二区三区国产精品乱码| 国产精品一区二区免费欧美| www.熟女人妻精品国产| 成熟少妇高潮喷水视频| 一级黄片播放器| 丰满人妻熟妇乱又伦精品不卡| 国产成人啪精品午夜网站| 成人鲁丝片一二三区免费| 无人区码免费观看不卡| 亚洲欧美日韩高清专用| 十八禁人妻一区二区| 丝袜美腿在线中文| 精品久久久久久久毛片微露脸| АⅤ资源中文在线天堂| 精品乱码久久久久久99久播| 国产免费av片在线观看野外av| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利欧美成人| 啦啦啦免费观看视频1| 国产高潮美女av| 国内精品久久久久精免费| 我的老师免费观看完整版| 欧美又色又爽又黄视频| 日韩欧美三级三区| 亚洲av电影在线进入| 亚洲av五月六月丁香网| 操出白浆在线播放| 精品国产美女av久久久久小说| 国产伦一二天堂av在线观看| 最近在线观看免费完整版| 日韩欧美一区二区三区在线观看| 内射极品少妇av片p| 亚洲成av人片在线播放无| 国产成人a区在线观看| 国产精品永久免费网站| 黄色丝袜av网址大全| 亚洲不卡免费看| 亚洲第一欧美日韩一区二区三区| 精品熟女少妇八av免费久了| 免费观看的影片在线观看| 禁无遮挡网站| 一本久久中文字幕| 亚洲最大成人手机在线| 久久久久国内视频| 午夜免费激情av| 18美女黄网站色大片免费观看| 啪啪无遮挡十八禁网站| 国产高清有码在线观看视频| 制服丝袜大香蕉在线| 天美传媒精品一区二区| 成人午夜高清在线视频| av视频在线观看入口| 成人午夜高清在线视频| 99国产精品一区二区蜜桃av| 日日干狠狠操夜夜爽| 亚洲欧美日韩卡通动漫| 两个人视频免费观看高清| 欧美激情在线99| 91久久精品国产一区二区成人 | 国产精品嫩草影院av在线观看 | 精品无人区乱码1区二区| 在线观看日韩欧美| 丁香欧美五月| 2021天堂中文幕一二区在线观| 免费在线观看影片大全网站| 99久久久亚洲精品蜜臀av| 一进一出抽搐gif免费好疼| 成人国产一区最新在线观看| 免费看a级黄色片| 丁香欧美五月| 国产主播在线观看一区二区| 色综合亚洲欧美另类图片| 岛国视频午夜一区免费看| 五月玫瑰六月丁香| 成人特级av手机在线观看| 亚洲av电影在线进入| 国产一区二区三区在线臀色熟女| 麻豆国产97在线/欧美| 久久久国产精品麻豆| 欧美一区二区国产精品久久精品| 久久久久久大精品| 亚洲av日韩精品久久久久久密| 俄罗斯特黄特色一大片| av在线天堂中文字幕| 国产精品女同一区二区软件 | 欧洲精品卡2卡3卡4卡5卡区| 色综合欧美亚洲国产小说| 亚洲av电影不卡..在线观看| 99久久成人亚洲精品观看| 日韩欧美国产在线观看| 国产探花在线观看一区二区| 国产精品野战在线观看| 麻豆一二三区av精品| 18+在线观看网站| 国产伦一二天堂av在线观看| 成人特级黄色片久久久久久久| 亚洲成人中文字幕在线播放| 美女免费视频网站| 他把我摸到了高潮在线观看| 男女床上黄色一级片免费看| 琪琪午夜伦伦电影理论片6080| 国产精品女同一区二区软件 | 丁香六月欧美| 国内精品一区二区在线观看| 99视频精品全部免费 在线| 日韩精品青青久久久久久| 3wmmmm亚洲av在线观看| 日韩欧美一区二区三区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 69人妻影院| 麻豆国产av国片精品| 中文字幕精品亚洲无线码一区| 欧美成人免费av一区二区三区| 精品久久久久久久久久久久久| 国产激情欧美一区二区| 在线十欧美十亚洲十日本专区| 久久国产精品影院| 国产爱豆传媒在线观看| 亚洲片人在线观看| 九色成人免费人妻av| 狠狠狠狠99中文字幕| 亚洲精品色激情综合| 88av欧美| 少妇丰满av| 亚洲av成人av| 岛国在线免费视频观看| 给我免费播放毛片高清在线观看| 在线看三级毛片| 九色国产91popny在线| 亚洲va日本ⅴa欧美va伊人久久| 动漫黄色视频在线观看| av天堂在线播放| av欧美777| 亚洲18禁久久av| 久久久久久国产a免费观看| 日日干狠狠操夜夜爽| av女优亚洲男人天堂| 国产亚洲欧美98| 亚洲天堂国产精品一区在线| 脱女人内裤的视频| 亚洲成av人片在线播放无| 免费看a级黄色片| 少妇裸体淫交视频免费看高清| 最新在线观看一区二区三区| 国产精品国产高清国产av| 国产一区二区三区视频了| 日韩欧美国产一区二区入口| 久久香蕉精品热| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美国产一区二区入口| 欧美bdsm另类| 午夜a级毛片| 亚洲国产欧美人成| 欧美一区二区亚洲| 久久天躁狠狠躁夜夜2o2o| 男人和女人高潮做爰伦理| 又爽又黄无遮挡网站| 亚洲久久久久久中文字幕| 人妻丰满熟妇av一区二区三区| 特级一级黄色大片| 国产99白浆流出| 亚洲自拍偷在线| 九九久久精品国产亚洲av麻豆| 九九热线精品视视频播放| 免费人成视频x8x8入口观看| 久久久色成人| 色综合欧美亚洲国产小说| 少妇高潮的动态图| 成人高潮视频无遮挡免费网站| 亚洲国产高清在线一区二区三| 嫩草影视91久久| 香蕉av资源在线| av在线天堂中文字幕| 成人无遮挡网站| 麻豆国产av国片精品| 日本黄色片子视频| 国产 一区 欧美 日韩| 一区二区三区国产精品乱码| 免费无遮挡裸体视频| 日韩精品青青久久久久久| 国产91精品成人一区二区三区| 90打野战视频偷拍视频| 亚洲七黄色美女视频| 禁无遮挡网站| 成熟少妇高潮喷水视频| 免费一级毛片在线播放高清视频| 欧美日本视频| 亚洲成人免费电影在线观看| 亚洲精品亚洲一区二区| 免费av毛片视频| 欧美性感艳星| 男插女下体视频免费在线播放| 欧美另类亚洲清纯唯美| 中文字幕高清在线视频| 亚洲电影在线观看av| 久久精品国产自在天天线| or卡值多少钱| 日韩亚洲欧美综合| 亚洲av不卡在线观看| 波多野结衣巨乳人妻| 国产精品久久久久久久电影 | 久久久色成人| 内地一区二区视频在线| 老司机深夜福利视频在线观看| 国产三级中文精品| 美女黄网站色视频| 久久天躁狠狠躁夜夜2o2o| 女人被狂操c到高潮| 天堂网av新在线| 国产精品国产高清国产av| 此物有八面人人有两片| 免费高清视频大片| 一边摸一边抽搐一进一小说| 精品午夜福利视频在线观看一区| 国产欧美日韩精品一区二区| 亚洲电影在线观看av| 热99在线观看视频| 在线观看美女被高潮喷水网站 | 在线视频色国产色| 久久精品综合一区二区三区| 午夜福利在线观看吧| 亚洲精品456在线播放app | 亚洲美女视频黄频| 一个人免费在线观看的高清视频| 欧美中文日本在线观看视频| 国产精品国产高清国产av| 一区福利在线观看| 久久6这里有精品| 在线观看午夜福利视频| 色综合欧美亚洲国产小说| 校园春色视频在线观看| 欧美大码av| 亚洲av五月六月丁香网| 日韩欧美国产一区二区入口| 一进一出好大好爽视频| 国产伦人伦偷精品视频| 在线看三级毛片| 真人做人爱边吃奶动态| 国产视频一区二区在线看| 国产男靠女视频免费网站| 在线免费观看的www视频| 久久亚洲真实| 人妻久久中文字幕网| 国产亚洲精品综合一区在线观看| 国产高潮美女av| 亚洲18禁久久av| 成人无遮挡网站| 成人亚洲精品av一区二区| 国产av麻豆久久久久久久| 中文亚洲av片在线观看爽| 搡女人真爽免费视频火全软件 | 国产伦在线观看视频一区| 国产亚洲精品久久久久久毛片| 五月伊人婷婷丁香| 天美传媒精品一区二区| 色视频www国产| 香蕉久久夜色| 国产真实伦视频高清在线观看 | 十八禁网站免费在线| 欧美日韩国产亚洲二区| 欧美国产日韩亚洲一区| 亚洲成人久久爱视频| 美女cb高潮喷水在线观看| 欧美一级毛片孕妇| 亚洲成av人片在线播放无| av在线蜜桃| 19禁男女啪啪无遮挡网站| 男插女下体视频免费在线播放| 亚洲人成网站在线播| 老司机深夜福利视频在线观看| 亚洲精品亚洲一区二区| 色吧在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲国产中文字幕在线视频| 色综合婷婷激情| 每晚都被弄得嗷嗷叫到高潮| 久久精品影院6| 色综合站精品国产| av视频在线观看入口| av片东京热男人的天堂| 小蜜桃在线观看免费完整版高清| 在线天堂最新版资源| 久久久成人免费电影| 制服丝袜大香蕉在线| 无人区码免费观看不卡| 成人精品一区二区免费| 一区二区三区高清视频在线| av欧美777| 99久久无色码亚洲精品果冻| 欧美日韩亚洲国产一区二区在线观看| 欧美乱码精品一区二区三区| 成人av一区二区三区在线看| 19禁男女啪啪无遮挡网站| 亚洲人与动物交配视频| 真人做人爱边吃奶动态| 精品人妻偷拍中文字幕| 久久精品综合一区二区三区| 可以在线观看的亚洲视频| 一级毛片高清免费大全| 亚洲国产精品999在线| 99国产精品一区二区蜜桃av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 嫁个100分男人电影在线观看| 亚洲第一电影网av| 亚洲国产高清在线一区二区三| 国产高清视频在线观看网站| 在线免费观看不下载黄p国产 | 88av欧美| 久久午夜亚洲精品久久| 看免费av毛片| 日韩欧美国产一区二区入口| 在线观看美女被高潮喷水网站 | 变态另类丝袜制服| 99riav亚洲国产免费| 亚洲av成人精品一区久久| 国产高潮美女av| 欧美激情在线99| 欧美最黄视频在线播放免费| 九九在线视频观看精品| 岛国视频午夜一区免费看| 18禁黄网站禁片午夜丰满| 内射极品少妇av片p| 搡老岳熟女国产| 亚洲av熟女| 亚洲在线自拍视频| 久久久久亚洲av毛片大全| 久久6这里有精品| 久久久久久久亚洲中文字幕 | 全区人妻精品视频| 亚洲av免费在线观看| 欧美日韩一级在线毛片| 亚洲无线观看免费| 免费av毛片视频| 午夜激情欧美在线| 久久久久国产精品人妻aⅴ院| 国产精品av视频在线免费观看| 桃红色精品国产亚洲av| 成人欧美大片| avwww免费| 高清毛片免费观看视频网站| 日韩大尺度精品在线看网址| 女人高潮潮喷娇喘18禁视频| 男女视频在线观看网站免费| 日韩大尺度精品在线看网址| 欧美日本亚洲视频在线播放| 免费看十八禁软件| 精品久久久久久,| 18禁在线播放成人免费| 久久草成人影院| 国产一级毛片七仙女欲春2| 18禁在线播放成人免费| 色尼玛亚洲综合影院| 欧美一区二区精品小视频在线| 香蕉丝袜av| 国产 一区 欧美 日韩| 日韩人妻高清精品专区| 91在线精品国自产拍蜜月 | 91av网一区二区| 性色avwww在线观看| 99精品在免费线老司机午夜| 国产真实伦视频高清在线观看 | 日韩国内少妇激情av| 欧美黑人欧美精品刺激| 欧美3d第一页| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区亚洲| 日韩大尺度精品在线看网址| 亚洲真实伦在线观看| 日本 欧美在线| 中文字幕av在线有码专区| 亚洲成人中文字幕在线播放| 高清毛片免费观看视频网站| 色综合婷婷激情| 久久这里只有精品中国| 欧美日韩一级在线毛片| 日韩精品中文字幕看吧| 亚洲18禁久久av| 国产精品亚洲av一区麻豆| 变态另类成人亚洲欧美熟女| 美女高潮喷水抽搐中文字幕| 桃色一区二区三区在线观看| 女同久久另类99精品国产91| 欧美性猛交黑人性爽| 一卡2卡三卡四卡精品乱码亚洲| 久久精品91蜜桃| 欧美日韩乱码在线| 免费在线观看日本一区| 欧美在线黄色| 熟女人妻精品中文字幕| 老司机在亚洲福利影院| 国产国拍精品亚洲av在线观看 | 亚洲人成网站高清观看| 久久欧美精品欧美久久欧美| 精品久久久久久久末码| 免费av观看视频| 国产精品久久久人人做人人爽| 久久久久性生活片| 欧美三级亚洲精品| 国产伦人伦偷精品视频| 日日摸夜夜添夜夜添小说| 日日干狠狠操夜夜爽| 久久久久免费精品人妻一区二区| 欧美一区二区国产精品久久精品| 真人做人爱边吃奶动态| 亚洲熟妇中文字幕五十中出| 免费看日本二区| av在线蜜桃| 嫩草影院入口| 亚洲中文日韩欧美视频| 精品国产美女av久久久久小说| 热99re8久久精品国产| h日本视频在线播放| 亚洲av电影在线进入| 国产在视频线在精品| 亚洲欧美日韩无卡精品| 欧美日韩福利视频一区二区| 国产黄片美女视频| 国产精品精品国产色婷婷| 九色国产91popny在线| 中文字幕人成人乱码亚洲影| 很黄的视频免费| 国产免费一级a男人的天堂| 亚洲五月天丁香| 久久精品亚洲精品国产色婷小说| 国产欧美日韩精品亚洲av| 色尼玛亚洲综合影院| 无遮挡黄片免费观看| 热99re8久久精品国产| 夜夜看夜夜爽夜夜摸| 亚洲av第一区精品v没综合| 国产激情偷乱视频一区二区| 91麻豆精品激情在线观看国产| 蜜桃亚洲精品一区二区三区| 精品电影一区二区在线| 国产成人影院久久av| 白带黄色成豆腐渣| 午夜亚洲福利在线播放| 午夜免费男女啪啪视频观看 | 亚洲国产精品sss在线观看| 国产成人av教育| 12—13女人毛片做爰片一| 最后的刺客免费高清国语| 久久精品国产亚洲av涩爱 | 国产伦精品一区二区三区四那| 香蕉av资源在线| 精品久久久久久成人av| av天堂在线播放| 在线观看日韩欧美| 日本黄大片高清| 性色avwww在线观看| 欧美日韩精品网址| 性色av乱码一区二区三区2| 久久精品亚洲精品国产色婷小说| 国产一区二区三区在线臀色熟女| 美女被艹到高潮喷水动态| 88av欧美| 高清在线国产一区| 在线观看一区二区三区| 国产视频内射| 午夜福利在线观看免费完整高清在 | 午夜免费激情av| 婷婷亚洲欧美| 听说在线观看完整版免费高清| 99热精品在线国产| 97超视频在线观看视频| 国内久久婷婷六月综合欲色啪| 国产野战对白在线观看| 18禁黄网站禁片午夜丰满| 欧美三级亚洲精品| 91在线精品国自产拍蜜月 | 在线免费观看的www视频| 午夜福利成人在线免费观看| 日本与韩国留学比较| 三级毛片av免费| 中文字幕av成人在线电影| 久9热在线精品视频| 欧美日韩黄片免| 日本一二三区视频观看| 三级男女做爰猛烈吃奶摸视频| 久久久久久久精品吃奶| 国产探花极品一区二区| 俄罗斯特黄特色一大片| 麻豆一二三区av精品| 哪里可以看免费的av片| 两个人看的免费小视频| 国产精品一及| 亚洲不卡免费看| 国产高清激情床上av| 首页视频小说图片口味搜索| 国产三级中文精品| 好男人在线观看高清免费视频| 在线天堂最新版资源| 亚洲成av人片免费观看| 日本黄色视频三级网站网址| 禁无遮挡网站| 午夜两性在线视频| 亚洲精品久久国产高清桃花| 一个人看的www免费观看视频| av视频在线观看入口| 国产精品国产高清国产av| 欧美成人性av电影在线观看| 成人高潮视频无遮挡免费网站| 国产高清有码在线观看视频| 国产精品一区二区三区四区免费观看 | 国产精品影院久久| 女人被狂操c到高潮| 国产高清激情床上av| 高潮久久久久久久久久久不卡| 久久久久久久午夜电影| 国产成人av激情在线播放| 亚洲激情在线av| 国产色爽女视频免费观看| 久久久久久久久久黄片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产又黄又爽又无遮挡在线|