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

    基于自適應(yīng)TQWT與小波包奇異譜熵的滾動(dòng)軸承早期故障診斷

    2023-03-29 02:55:54謝鋒云劉慧胡旺姜永奇
    關(guān)鍵詞:峭度層數(shù)波包

    謝鋒云,劉慧,胡旺,姜永奇

    (華東交通大學(xué) 機(jī)電與車輛工程學(xué)院,江西 南昌 330013)

    滾動(dòng)軸承是機(jī)械傳動(dòng)設(shè)備的重要組成部分,廣泛應(yīng)用于各類交通運(yùn)輸設(shè)備、工業(yè)生產(chǎn)設(shè)備,其健康狀況與設(shè)備正常工作有著密不可分的關(guān)系,若在滾動(dòng)軸承早期故障階段就敏感地檢測(cè)到其故障的發(fā)生,可以達(dá)到減少經(jīng)濟(jì)損失以及保護(hù)人身安全的效果。在機(jī)械設(shè)備生產(chǎn)工作過(guò)程中,由于變轉(zhuǎn)速和變工況的影響,傳感器采集到的信號(hào)大多表現(xiàn)為非線性和非平穩(wěn)性。短時(shí)傅里葉變換(STFT)、經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)和小波分析等是常用的非線性、非平穩(wěn)信號(hào)處理方法。STFT 的窗函數(shù)一經(jīng)選定,在變換過(guò)程中便不能進(jìn)行更改,不能反映頻率隨時(shí)間的變化[1];EMD的理論基礎(chǔ)不夠完善,且在信號(hào)分解過(guò)程中存在端點(diǎn)效應(yīng)和模態(tài)混疊等問(wèn)題[2];小波分析具有優(yōu)良的時(shí)頻分析能力,受到了業(yè)內(nèi)人士的熱愛[3]。但運(yùn)用小波分析對(duì)信號(hào)進(jìn)行處理的效果依賴小波基的選擇,且在處理過(guò)程中小波基一旦選定便不能更改,選擇不當(dāng)對(duì)信號(hào)的處理效果會(huì)有很大的影響?;诖?,SELESNICK[4]提出一種新的品質(zhì)因子可調(diào)的時(shí)頻分析方法,即可調(diào)品質(zhì)因子小波變換(TQWT),它能夠通過(guò)調(diào)節(jié)參數(shù)品質(zhì)因子Q和冗余因子r,為信號(hào)匹配符合其振蕩特性的小波基函數(shù),達(dá)到更好的分解效果。龍瑩等[5]采用改進(jìn)譜峭度(Improved spectral kurtosis, ISK)的TQWT,對(duì)齒輪箱軸承故障進(jìn)行了有效診斷。ZHANG 等[6]采用TQWT 對(duì)信號(hào)進(jìn)行處理,并利用頻域定義的1.5D 譜峭度分離出有用的故障特征頻帶和干擾頻帶。KONG 等[7]提出一種新的基于自適應(yīng)TQWT 濾波器的特征提取方法來(lái)檢測(cè)重復(fù)瞬態(tài),將該方法應(yīng)用于故障軸承振動(dòng)信號(hào),驗(yàn)證了其有效性。小波包奇異譜熵是將小波包分解與奇異譜熵相結(jié)合的一種故障特征提取方法。小波包分解不僅可以對(duì)信號(hào)的全頻率成分進(jìn)行較為完整的分解[8],還可根據(jù)信號(hào)特性和分析要求自適應(yīng)選擇相匹配的頻帶與頻譜,從而從信號(hào)中更有效地提取信息[9]。奇異譜熵能夠反映信號(hào)的復(fù)雜性與時(shí)序性。周建民等[10]通過(guò)將小波包奇異譜熵和支持向量數(shù)據(jù)描述(SVDD)相結(jié)合對(duì)滾動(dòng)軸承性能退化進(jìn)行了評(píng)估。針對(duì)滾動(dòng)軸承早期故障的特點(diǎn),本文提出一種將TQWT 與小波包奇異譜熵相結(jié)合提取故障特征,SVM 模式識(shí)別的滾動(dòng)軸承故障診斷方法。利用網(wǎng)格搜索法對(duì)TQWT的最佳品質(zhì)因子進(jìn)行搜索,并以峭譜積最大化為評(píng)價(jià)標(biāo)準(zhǔn)確定最佳品質(zhì)因子,利用中心頻率比確定合適的分解層數(shù),將得到的參數(shù)對(duì)信號(hào)進(jìn)行TQWT 分解并單支重構(gòu),將峭譜積最大的重構(gòu)分量的小波包奇異譜熵值作為故障特征向量,運(yùn)用SVM 進(jìn)行故障狀態(tài)的判斷,并與未進(jìn)行TQWT 處理、以峭度和包絡(luò)譜峰值因子單獨(dú)作為TQWT 優(yōu)化指標(biāo)的信號(hào)的SVM 識(shí)別結(jié)果比較,對(duì)比結(jié)果顯示,本文所提方法處理后的信號(hào)的識(shí)別率得到了一定程度的提升。

    1 基本原理

    1.1 TQWT

    TQWT 是一種從頻域?yàn)V波角度設(shè)計(jì)而成的離散小波變換方法,與傳統(tǒng)小波變換品質(zhì)因子(信號(hào)的中心頻率fc 和帶寬BW 的比值)恒定相比,它根據(jù)待分析信號(hào)的振蕩特性來(lái)選擇品質(zhì)因子,通過(guò)迭代雙通道濾波器組和離散傅里葉變換(DFT)以迭代的方式實(shí)現(xiàn)信號(hào)的分解、重構(gòu)[11]。

    TQWT 由一個(gè)個(gè)小波基函數(shù)組成,每個(gè)基函數(shù)由3 個(gè)可調(diào)的參數(shù)(品質(zhì)因子Q,冗余度r和分解層數(shù)J)控制,參數(shù)不同,輸出的小波基函數(shù)的形狀和性質(zhì)也不同。其中,Q和r與低通尺度因子α以及高通尺度因子β的關(guān)系滿足:

    TQWT 的Q和r確定了,可以通過(guò)式(2)確定TQWT理論允許的最大分解層數(shù)Jmax。

    其中:[·]為向0取整;N為待分析信號(hào)的長(zhǎng)度。

    品質(zhì)因子Q定義為濾波器中心頻率與帶寬的比值。它用來(lái)描述信號(hào)的振蕩屬性(小波振蕩持續(xù)的時(shí)間),可以反映信號(hào)的振蕩程度,Q越大,振蕩越強(qiáng)烈。品質(zhì)因子Q的選擇決定著最終濾波頻帶的選擇,合適的頻帶評(píng)價(jià)指標(biāo)是保證最終濾波效果的關(guān)鍵。本文使用能夠表征信號(hào)沖擊性和周期性的時(shí)頻域指標(biāo)——包絡(luò)譜峰值因子(Crest factor of envelope spectrum, Ec)[12]與信號(hào)的峭度值的乘積——峭譜積,對(duì)最佳的Q進(jìn)行挑選。根據(jù)文獻(xiàn)[11, 13],本文Q的遍歷范圍選擇為[1,5]。

    r是計(jì)算無(wú)窮多級(jí)時(shí)TQWT 的冗余度,通常r≥3。冗余因子決定了小波頻響間的重合度,r越大,則頻響重合度越高。根據(jù)文獻(xiàn)[14-15],本文設(shè)定r=3。

    分解層數(shù)J表示雙通道濾波器組的數(shù)目。分解層數(shù)過(guò)大不僅會(huì)降低運(yùn)算速度還會(huì)導(dǎo)致冗余信號(hào)分量過(guò)多,這些分量由于帶寬太窄而無(wú)法囊括足夠的特征信息,對(duì)于故障特征提取來(lái)說(shuō)意義不大,并且會(huì)造成計(jì)算資源的浪費(fèi)。本文使用中心頻率比來(lái)確定合適的分解層數(shù)[13]。

    1.2 峭譜積

    包絡(luò)譜峰值因子考慮了共振頻段能量較大且濾波信號(hào)周期性較強(qiáng)的特點(diǎn) 。定義為:

    式中:Amax為包絡(luò)譜在[n*fr,fs/2]范圍內(nèi)的最大幅值;Arms為均方根值。Ec值越大,信號(hào)包含的周期性和沖擊性越強(qiáng),反之亦然。

    峭度是表征信號(hào)的沖擊性強(qiáng)弱的時(shí)域指標(biāo),其計(jì)算公式為:

    其中:μ為信號(hào)的均值;σ為信號(hào)的標(biāo)準(zhǔn)差。

    峭譜積的表達(dá)式為:

    它在包絡(luò)譜峰值因子的基礎(chǔ)上,加強(qiáng)了對(duì)信號(hào)沖擊性的度量,KEc越大,信號(hào)的周期性和沖擊性越明顯,反之亦然。

    1.3 中心頻率比

    中心頻率比的定義為:

    合理分解層數(shù)停止準(zhǔn)則表示為:

    其中:T為合理的中心頻率比閾值,推薦選擇值為0.02~0.05[13],本文設(shè)置閾值T=0.05。因此,合理的分解層數(shù)Ja選擇為:

    1.4 小波包奇異譜熵

    1.4.1 小波包奇異譜熵原理

    小波包奇異譜熵是將小波包分解、奇異值分解與信息熵融為一體的特征提取方法。求取信號(hào)的小波包奇異譜熵的步驟為:對(duì)信號(hào)進(jìn)行小波包分解,對(duì)分解后的系數(shù)進(jìn)行重構(gòu),并將重構(gòu)后的系數(shù)矩陣進(jìn)行奇異值分解,求歸一化后的奇異值的信息熵,從而對(duì)信號(hào)的復(fù)雜程度給出一個(gè)確定的量度[10]。

    1.4.2 小波包分解層數(shù)及小波基的選擇

    由文獻(xiàn)[16]可知,如果只是檢測(cè)信號(hào)奇變的突變點(diǎn),按照規(guī)則性系數(shù)相似方法,選擇規(guī)則性系數(shù)較小的小波基波會(huì)得到較好的結(jié)果,所以本文選擇db10小波基進(jìn)行小波包分解。

    王婉秋等[17]指出,對(duì)于滾動(dòng)軸承振動(dòng)信號(hào),小波包的分解層數(shù)應(yīng)該滿足:

    式中:J為最大分解層數(shù);fs為采樣頻率;f為信號(hào)的故障頻率。經(jīng)過(guò)計(jì)算,本文所用信號(hào)的小波包分解的最大分解層數(shù)為6層,本文小波包分解層數(shù)選擇3層。

    2 故障診斷流程

    為提取滾動(dòng)軸承早期故障的微弱故障特征,對(duì)滾動(dòng)軸承的早期故障狀態(tài)進(jìn)行有效識(shí)別,本文提出一種自適應(yīng)可調(diào)品質(zhì)因子與小波包奇異譜熵相結(jié)合提取故障特征的方法,運(yùn)用支持向量機(jī)對(duì)提取的故障特征進(jìn)行模式識(shí)別,診斷流程如圖1所示。

    圖1 故障診斷流程圖Fig. 1 Fault diagnosis flow chart

    具體的步驟為:

    1) 確定TQWT最佳的品質(zhì)因子Q。

    依據(jù)文獻(xiàn)[11, 13],設(shè)定Q的尋優(yōu)范圍為[1, 5],遍歷步長(zhǎng)為0.1。分別對(duì)41種情況下的[Q,r]組合進(jìn)行TQWT 分解,并單支重構(gòu),計(jì)算單支重構(gòu)后各重構(gòu)分量的峭譜積KEc,最大的KEc對(duì)應(yīng)的Q值即為最佳的Q值。

    2) 利用中心頻率比確定TQWT 合適的分解層數(shù)J。

    3) 利用得到的Q和J對(duì)待分析的信號(hào)進(jìn)行TQWT 分解,并單支重構(gòu),計(jì)算各重構(gòu)分量的KEc,挑選出每組信號(hào)KEc最大對(duì)應(yīng)的重構(gòu)分量為新的待分析的信號(hào)。

    4) 對(duì)得到的新的待分析信號(hào)提取小波包奇異譜熵值,作為故障特征向量。

    5) 運(yùn)用支持向量機(jī)對(duì)得到的故障特征向量進(jìn)行模式識(shí)別。

    3 實(shí)驗(yàn)平臺(tái)與實(shí)驗(yàn)分析

    運(yùn)用XJTU-SY 滾動(dòng)軸承加速壽命試驗(yàn)數(shù)據(jù)集驗(yàn)證本文所提方法的有效性,試驗(yàn)所用的軸承平臺(tái)如圖2 所示[18-19]。該平臺(tái)由數(shù)字式力顯示器、電機(jī)速度控制器、支撐軸承、測(cè)試軸承、交流電機(jī)、液壓加載系統(tǒng)以及加速度傳感器等設(shè)備搭建而成,可用來(lái)測(cè)試在不同工況下的軸承的加速度振動(dòng)信號(hào)。實(shí)驗(yàn)驗(yàn)證所用數(shù)據(jù)的具體信息見表1。選取外圈、內(nèi)圈和保持架故障產(chǎn)生初期的信號(hào)作為早期故障信號(hào),加上軸承正常運(yùn)行時(shí)的數(shù)據(jù),共4種運(yùn)行狀態(tài),運(yùn)用本文所提方法進(jìn)行分析驗(yàn)證。

    表1 實(shí)驗(yàn)相關(guān)參數(shù)信息Table 1 Experimental related parameter information

    圖2 試驗(yàn)平臺(tái)Fig. 2 Test platform

    由采樣頻率和轉(zhuǎn)速等信息可計(jì)算出軸承轉(zhuǎn)一圈大概會(huì)產(chǎn)生683個(gè)采樣點(diǎn),為保證每組樣本包含足夠的故障信息,設(shè)置每組樣本長(zhǎng)度為2 048 個(gè)采樣點(diǎn),每種狀態(tài)200 組樣本,共800 組樣本。任選保持架故障的一組樣本進(jìn)行驗(yàn)證展示,圖3為該組保持架故障原始信號(hào)的時(shí)域波形圖。

    圖3 原始信號(hào)時(shí)域圖Fig. 3 Time domain diagram of original signal

    對(duì)41 種[Q,r]組合進(jìn)行TQWT 分解與單支重構(gòu),計(jì)算重構(gòu)后分量的KEc值,圖4 為41 種[Q,r]組合下得到的該組樣本的KEc最大值的曲線圖。從圖中可以得到,最大KEc值對(duì)應(yīng)的Q值為Q=1+利用得到的Q值,結(jié)合式(2)計(jì)算出該組樣本TQWT 分解時(shí)的最大分解層數(shù)為35,利用中心頻率比計(jì)算合適的分解層數(shù)為24。

    圖4 不同[Q,r]組合下的最大KEc值曲線Fig. 4 Maximum KEc curves under different [Q,r] combinations

    對(duì)待處理信號(hào)運(yùn)用得到的[Q,r,J]值進(jìn)行TQWT分解并單支重構(gòu),計(jì)算重構(gòu)后各分量的KEc值,挑選最大的KEc值對(duì)應(yīng)的重構(gòu)分量為新的待處理的信號(hào)。圖5為TQWT處理后的信號(hào)的時(shí)域波形圖。可以看到TQWT 處理后信號(hào)的幅值有所下降,但信號(hào)的周期性得到了增強(qiáng)。

    圖5 TQWT處理后信號(hào)時(shí)域圖Fig. 5 Time domain diagram of TQWT processed signal

    對(duì)TQWT 處理后的信號(hào)提取小波包奇異譜熵作為故障特征向量。圖6 為4 種狀態(tài)TQWT 處理后一組信號(hào)的小波包奇異譜熵值對(duì)比,圖7 為4 種狀態(tài)原始信號(hào)的小波包奇異譜熵值對(duì)比。從2幅圖的對(duì)比可以看出,TQWT 處理后4種狀態(tài)的小波包奇異譜熵值的區(qū)分度更高。

    圖6 TQWT處理后4種狀態(tài)的小波包奇異譜熵值對(duì)比Fig. 6 Comparison of singular spectral entropy values of wavelet packets in four states after TQWT processing

    圖7 4種狀態(tài)原始信號(hào)的小波包奇異譜熵值對(duì)比Fig. 7 Comparison of the singular spectral entropy values of wavelet packets of four original signals

    對(duì)得到的故障特征向量運(yùn)用SVM 進(jìn)行故障狀態(tài)的識(shí)別,其中50組用于模型的訓(xùn)練,150組用于測(cè)試。圖8 為TQWT 處理后的信號(hào)的SVM 識(shí)別結(jié)果。圖9 為原始信號(hào)的SVM 識(shí)別結(jié)果,圖10 和圖11 分別為TQWT 優(yōu)化指標(biāo)為峭度和包絡(luò)譜峰值因子的識(shí)別結(jié)果。

    圖8 TQWT處理后信號(hào)SVM識(shí)別結(jié)果Fig. 8 SVM recognition result of TQWT processed signal

    圖9 原始信號(hào)識(shí)別結(jié)果Fig. 9 Original signal recognition result

    圖10 TQWT優(yōu)化指標(biāo)為峭度識(shí)別結(jié)果Fig. 10 TQWT optimization index is kurtosis recognition result

    圖11 TQWT優(yōu)化指標(biāo)為Ec識(shí)別結(jié)果Fig. 11 TQWT optimization index is Ec identification result

    為使上述識(shí)別結(jié)果更加直觀,引入了混沌矩陣,詳細(xì)分析了本文所提方法、未進(jìn)行TQWT 處理、TQWT 優(yōu)化指標(biāo)為峭度和包絡(luò)譜峰值因子的識(shí)別結(jié)果,如圖12所示。

    圖12 4種模型識(shí)別結(jié)果混沌矩陣Fig. 12 Confusion matrix of identification results of four models

    從以上SVM 的處理結(jié)果對(duì)比發(fā)現(xiàn),利用本文所提峭譜積為優(yōu)化指標(biāo)對(duì)信號(hào)TQWT 處理后,識(shí)別率較未進(jìn)行TQWT 處理的提高了約7%,與TQWT 的優(yōu)化指標(biāo)為峭度和Ec 的識(shí)別結(jié)果相比識(shí)別率也有提高,但以峭度和Ec 為優(yōu)化指標(biāo)時(shí),雖然識(shí)別率與未進(jìn)行TQWT 處理時(shí)相比有了提高,但出現(xiàn)了將內(nèi)圈和外圈故障識(shí)別為正常狀態(tài)的樣本,這在故障診斷領(lǐng)域是不被允許的。

    4 結(jié)論

    1) 提出一種基于自適應(yīng)TQWT 與小波包奇異譜熵相結(jié)合提取特征的滾動(dòng)軸承早期故障診斷方法,可將原始信號(hào)不易被檢測(cè)的微弱信息進(jìn)行放大,有效診斷出滾動(dòng)軸承所處的工作狀態(tài)。

    2) 與確定TQWT 的最佳品質(zhì)因子并篩選最佳分量的指標(biāo)為峭度、包絡(luò)譜峰值因子相比,指標(biāo)為本文所提的KEc時(shí)的最終故障診斷識(shí)別率是最高的,且沒有將故障樣本識(shí)別為正常樣本的情況出現(xiàn)。

    猜你喜歡
    峭度層數(shù)波包
    基于MCKD和峭度的液壓泵故障特征提取
    填筑層數(shù)對(duì)土石壩應(yīng)力變形的影響研究
    上海發(fā)布藥品包裝物減量指南
    康復(fù)(2022年31期)2022-03-23 20:39:56
    聯(lián)合快速峭度圖與變帶寬包絡(luò)譜峭度圖的輪對(duì)軸承復(fù)合故障檢測(cè)研究
    MoS2薄膜電子性質(zhì)隨層數(shù)變化的理論研究
    電子制作(2019年11期)2019-07-04 00:34:50
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    基于峭度分析的聲發(fā)射故障檢測(cè)
    電子世界(2018年12期)2018-07-04 06:34:38
    基于魯棒性小波包峭度圖的滾動(dòng)軸承故障診斷*
    基于小波包變換的電力系統(tǒng)諧波分析
    小波包理論與圖像小波包分解
    成人午夜精彩视频在线观看| 18禁裸乳无遮挡免费网站照片| 精品人妻熟女av久视频| 男女那种视频在线观看| 插逼视频在线观看| 亚洲国产av新网站| 18禁在线无遮挡免费观看视频| 成人二区视频| 久久久久国产网址| 亚洲熟女精品中文字幕| 成人午夜高清在线视频| 中文字幕亚洲精品专区| 天天一区二区日本电影三级| 只有这里有精品99| 日本爱情动作片www.在线观看| 乱系列少妇在线播放| 亚洲在线自拍视频| 91久久精品电影网| 国产黄a三级三级三级人| 大香蕉97超碰在线| av在线老鸭窝| 色网站视频免费| 免费播放大片免费观看视频在线观看| 国产免费视频播放在线视频 | 亚洲精品久久久久久婷婷小说| 又粗又硬又长又爽又黄的视频| 精品久久久精品久久久| 国产一区有黄有色的免费视频 | 美女xxoo啪啪120秒动态图| 精品人妻熟女av久视频| 日日摸夜夜添夜夜添av毛片| 日本色播在线视频| 国产黄片视频在线免费观看| 午夜精品国产一区二区电影 | 老师上课跳d突然被开到最大视频| 丝瓜视频免费看黄片| www.av在线官网国产| 日韩欧美 国产精品| 精华霜和精华液先用哪个| 久久久国产一区二区| 一夜夜www| 最近2019中文字幕mv第一页| 久久久久久久亚洲中文字幕| 波野结衣二区三区在线| 自拍偷自拍亚洲精品老妇| 99热网站在线观看| 超碰av人人做人人爽久久| 亚洲av免费高清在线观看| 中文字幕制服av| 人人妻人人看人人澡| 国产一级毛片七仙女欲春2| 久久久久久久久大av| 久久精品久久精品一区二区三区| 少妇熟女欧美另类| 亚洲av不卡在线观看| 午夜视频国产福利| 在线a可以看的网站| 国产乱人视频| 3wmmmm亚洲av在线观看| 国产精品一及| 久久99热这里只频精品6学生| 超碰av人人做人人爽久久| av卡一久久| 日本爱情动作片www.在线观看| av国产久精品久网站免费入址| 建设人人有责人人尽责人人享有的 | 一级片'在线观看视频| 国产成人免费观看mmmm| 亚洲内射少妇av| 国产综合精华液| 2021少妇久久久久久久久久久| 麻豆乱淫一区二区| 尤物成人国产欧美一区二区三区| 久久久精品免费免费高清| 国产精品1区2区在线观看.| 26uuu在线亚洲综合色| 少妇猛男粗大的猛烈进出视频 | 黄色日韩在线| 免费av毛片视频| 最近手机中文字幕大全| 青春草视频在线免费观看| 国产免费一级a男人的天堂| 国产国拍精品亚洲av在线观看| av女优亚洲男人天堂| 成年女人在线观看亚洲视频 | 综合色av麻豆| 成年av动漫网址| av天堂中文字幕网| 有码 亚洲区| 国产中年淑女户外野战色| 欧美高清成人免费视频www| 只有这里有精品99| 一区二区三区免费毛片| 亚洲丝袜综合中文字幕| 亚洲精品成人久久久久久| 啦啦啦中文免费视频观看日本| 国产91av在线免费观看| 国产熟女欧美一区二区| 中国国产av一级| 色视频www国产| 内射极品少妇av片p| 看非洲黑人一级黄片| 欧美成人a在线观看| av黄色大香蕉| 亚洲精品自拍成人| 欧美激情在线99| 亚洲av在线观看美女高潮| 欧美日韩综合久久久久久| 国产 亚洲一区二区三区 | 亚洲精品久久午夜乱码| 18禁裸乳无遮挡免费网站照片| 国产在视频线在精品| 你懂的网址亚洲精品在线观看| 天堂√8在线中文| 亚洲人与动物交配视频| 99re6热这里在线精品视频| 日韩伦理黄色片| 最后的刺客免费高清国语| 国产精品麻豆人妻色哟哟久久 | 国产精品伦人一区二区| 午夜福利在线观看免费完整高清在| 精品久久久久久久久av| 亚洲精品中文字幕在线视频 | 日本av手机在线免费观看| 亚洲国产精品成人久久小说| 国产精品美女特级片免费视频播放器| 草草在线视频免费看| 激情 狠狠 欧美| 91精品伊人久久大香线蕉| av国产久精品久网站免费入址| 简卡轻食公司| 亚洲美女搞黄在线观看| 人人妻人人看人人澡| 人妻夜夜爽99麻豆av| 国产精品av视频在线免费观看| 有码 亚洲区| 亚洲婷婷狠狠爱综合网| 国产精品国产三级国产专区5o| 亚洲美女搞黄在线观看| 精品少妇黑人巨大在线播放| av网站免费在线观看视频 | av一本久久久久| 一本久久精品| 内地一区二区视频在线| 色5月婷婷丁香| 日韩av在线免费看完整版不卡| 最近最新中文字幕免费大全7| 国精品久久久久久国模美| 人妻系列 视频| 午夜免费观看性视频| 亚洲欧美日韩无卡精品| av卡一久久| 欧美97在线视频| 国产三级在线视频| 午夜精品在线福利| 特级一级黄色大片| 久久久国产一区二区| 麻豆av噜噜一区二区三区| 国产白丝娇喘喷水9色精品| 亚洲精品乱久久久久久| 97精品久久久久久久久久精品| 十八禁网站网址无遮挡 | 中文字幕av在线有码专区| 美女大奶头视频| 日韩欧美精品免费久久| 国产精品人妻久久久久久| 少妇被粗大猛烈的视频| 亚洲成人久久爱视频| 伊人久久精品亚洲午夜| 高清av免费在线| 久久精品久久久久久久性| 国产白丝娇喘喷水9色精品| 欧美精品国产亚洲| 亚洲精品中文字幕在线视频 | 18禁在线无遮挡免费观看视频| 国产探花在线观看一区二区| 最新中文字幕久久久久| 久久精品国产亚洲网站| 亚洲欧美精品自产自拍| 少妇裸体淫交视频免费看高清| 亚洲av在线观看美女高潮| 22中文网久久字幕| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线观看播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 女的被弄到高潮叫床怎么办| 男女啪啪激烈高潮av片| 欧美xxxx黑人xx丫x性爽| 午夜视频国产福利| 亚洲av.av天堂| 在线免费观看不下载黄p国产| 免费无遮挡裸体视频| freevideosex欧美| 18禁裸乳无遮挡免费网站照片| 人人妻人人看人人澡| 草草在线视频免费看| 亚洲人与动物交配视频| 欧美三级亚洲精品| 午夜福利视频精品| 亚洲,欧美,日韩| 色吧在线观看| 99热全是精品| av天堂中文字幕网| 日日干狠狠操夜夜爽| 黄色欧美视频在线观看| 日韩不卡一区二区三区视频在线| 国产黄色视频一区二区在线观看| 欧美区成人在线视频| 99九九线精品视频在线观看视频| 一级毛片aaaaaa免费看小| 久久国产乱子免费精品| 秋霞在线观看毛片| 1000部很黄的大片| 能在线免费看毛片的网站| 欧美变态另类bdsm刘玥| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品人妻一区二区三区麻豆| 亚洲熟女精品中文字幕| 午夜激情福利司机影院| 三级国产精品欧美在线观看| 丰满少妇做爰视频| 黄色一级大片看看| 97在线视频观看| 直男gayav资源| 汤姆久久久久久久影院中文字幕 | 黄色一级大片看看| 免费看美女性在线毛片视频| 秋霞伦理黄片| 亚洲精品色激情综合| 亚洲精品中文字幕在线视频 | 亚洲婷婷狠狠爱综合网| 欧美日本视频| 午夜激情久久久久久久| 欧美极品一区二区三区四区| 免费看a级黄色片| 禁无遮挡网站| 夫妻性生交免费视频一级片| 欧美一区二区亚洲| a级一级毛片免费在线观看| 日韩制服骚丝袜av| 亚洲精品456在线播放app| av在线蜜桃| 亚洲aⅴ乱码一区二区在线播放| 成人美女网站在线观看视频| 国产精品久久久久久精品电影| 99热这里只有是精品在线观看| 天堂中文最新版在线下载 | 禁无遮挡网站| 成人亚洲欧美一区二区av| 国产精品女同一区二区软件| 免费观看精品视频网站| 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 成人国产麻豆网| 久久人人爽人人片av| 熟女电影av网| 在线观看免费高清a一片| 午夜久久久久精精品| 亚洲精品乱码久久久久久按摩| 一个人观看的视频www高清免费观看| 精品久久久噜噜| 久久这里只有精品中国| 国产视频内射| 亚洲精品亚洲一区二区| 免费观看的影片在线观看| 精品酒店卫生间| 亚洲精品国产成人久久av| 日本一本二区三区精品| 高清日韩中文字幕在线| 日本三级黄在线观看| 蜜臀久久99精品久久宅男| 老司机影院成人| av天堂中文字幕网| 久久97久久精品| 在线观看人妻少妇| 边亲边吃奶的免费视频| 又大又黄又爽视频免费| 成人一区二区视频在线观看| 日韩人妻高清精品专区| 国产黄色小视频在线观看| 国产麻豆成人av免费视频| 熟妇人妻久久中文字幕3abv| 国产精品无大码| 一二三四中文在线观看免费高清| 又爽又黄a免费视频| 日韩av不卡免费在线播放| av在线播放精品| 久久久久性生活片| 熟女电影av网| 女人十人毛片免费观看3o分钟| 可以在线观看毛片的网站| 久久精品夜夜夜夜夜久久蜜豆| 久久综合国产亚洲精品| 精品久久久久久久久久久久久| 亚洲精品自拍成人| 日本免费a在线| 国产 一区 欧美 日韩| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人综合色| 久久精品国产鲁丝片午夜精品| 国产黄频视频在线观看| 精品欧美国产一区二区三| 成人高潮视频无遮挡免费网站| 一级片'在线观看视频| 亚洲久久久久久中文字幕| 国产精品一区二区性色av| 综合色丁香网| 乱码一卡2卡4卡精品| 国产亚洲5aaaaa淫片| 亚洲国产最新在线播放| av在线亚洲专区| av又黄又爽大尺度在线免费看| 青春草国产在线视频| 亚洲第一区二区三区不卡| 一个人免费在线观看电影| 午夜精品一区二区三区免费看| 看非洲黑人一级黄片| 亚洲国产色片| 一个人看的www免费观看视频| 久久久久久国产a免费观看| 国产在线一区二区三区精| 日韩精品青青久久久久久| 一级毛片我不卡| 成人欧美大片| 一级毛片 在线播放| 韩国av在线不卡| 亚洲成色77777| 男女国产视频网站| av线在线观看网站| 亚洲精品,欧美精品| 美女cb高潮喷水在线观看| 天堂av国产一区二区熟女人妻| 国产不卡一卡二| 简卡轻食公司| 国产精品99久久久久久久久| www.色视频.com| 99热这里只有是精品在线观看| 国产高清不卡午夜福利| 亚洲精品乱久久久久久| 成年免费大片在线观看| 国产精品一区二区三区四区免费观看| 久久人人爽人人爽人人片va| 国产视频内射| 老司机影院毛片| 国产男人的电影天堂91| 国产单亲对白刺激| 午夜福利在线观看免费完整高清在| 淫秽高清视频在线观看| 熟妇人妻不卡中文字幕| 国产精品嫩草影院av在线观看| 伦精品一区二区三区| 又粗又硬又长又爽又黄的视频| 18禁动态无遮挡网站| av福利片在线观看| 噜噜噜噜噜久久久久久91| 乱码一卡2卡4卡精品| freevideosex欧美| 日日干狠狠操夜夜爽| 成人综合一区亚洲| 男女下面进入的视频免费午夜| 51国产日韩欧美| 国产成人freesex在线| 亚洲av成人av| 亚洲欧美中文字幕日韩二区| 欧美潮喷喷水| 国产精品久久久久久精品电影| 国产单亲对白刺激| 亚洲电影在线观看av| 日韩欧美精品v在线| 在线 av 中文字幕| 高清视频免费观看一区二区 | 亚洲18禁久久av| 91精品国产九色| 精品熟女少妇av免费看| 成人午夜高清在线视频| 成人一区二区视频在线观看| 夜夜看夜夜爽夜夜摸| 青春草国产在线视频| 色哟哟·www| 免费大片18禁| 亚洲久久久久久中文字幕| 99视频精品全部免费 在线| 亚洲成人中文字幕在线播放| 少妇裸体淫交视频免费看高清| 午夜福利成人在线免费观看| 熟女电影av网| 亚洲av成人av| 欧美丝袜亚洲另类| 亚洲国产成人一精品久久久| 欧美xxⅹ黑人| 久久99热这里只频精品6学生| 成年av动漫网址| 黄片无遮挡物在线观看| 免费大片18禁| 亚洲人成网站高清观看| 日韩一本色道免费dvd| 嫩草影院新地址| 蜜桃久久精品国产亚洲av| 全区人妻精品视频| 成人午夜高清在线视频| 中文在线观看免费www的网站| 免费高清在线观看视频在线观看| 亚洲欧美日韩卡通动漫| 观看免费一级毛片| 精品熟女少妇av免费看| 丰满人妻一区二区三区视频av| 国产69精品久久久久777片| 国产老妇女一区| 午夜亚洲福利在线播放| 亚洲在线观看片| 日韩欧美精品v在线| 亚洲四区av| 久久久久久久久久久免费av| 淫秽高清视频在线观看| 中文在线观看免费www的网站| 国产女主播在线喷水免费视频网站 | 听说在线观看完整版免费高清| 国产 一区 欧美 日韩| 少妇被粗大猛烈的视频| 精品久久久久久成人av| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻一区二区三区麻豆| 高清日韩中文字幕在线| 联通29元200g的流量卡| 国产亚洲最大av| 欧美日本视频| 久久97久久精品| 日韩,欧美,国产一区二区三区| 91在线精品国自产拍蜜月| 欧美不卡视频在线免费观看| 69av精品久久久久久| 国产男女超爽视频在线观看| 精品欧美国产一区二区三| 永久免费av网站大全| 亚洲国产精品成人久久小说| 自拍偷自拍亚洲精品老妇| 韩国高清视频一区二区三区| 精品久久久久久久久亚洲| 亚洲久久久久久中文字幕| 亚洲av男天堂| 国产淫语在线视频| 成人漫画全彩无遮挡| 一个人看视频在线观看www免费| 啦啦啦中文免费视频观看日本| 人妻一区二区av| av在线播放精品| 精品酒店卫生间| 91精品伊人久久大香线蕉| 国产有黄有色有爽视频| 能在线免费观看的黄片| 国产真实伦视频高清在线观看| 女的被弄到高潮叫床怎么办| 色综合亚洲欧美另类图片| 十八禁国产超污无遮挡网站| 麻豆久久精品国产亚洲av| 一个人免费在线观看电影| av在线观看视频网站免费| 99热全是精品| 久久久久精品久久久久真实原创| 极品少妇高潮喷水抽搐| 床上黄色一级片| kizo精华| 亚洲aⅴ乱码一区二区在线播放| 三级男女做爰猛烈吃奶摸视频| 国产亚洲av嫩草精品影院| 午夜福利网站1000一区二区三区| 青春草国产在线视频| 亚洲成人中文字幕在线播放| 观看免费一级毛片| or卡值多少钱| 精品酒店卫生间| 国产亚洲精品久久久com| 一区二区三区高清视频在线| 成年免费大片在线观看| 成人国产麻豆网| 十八禁国产超污无遮挡网站| 噜噜噜噜噜久久久久久91| 久久鲁丝午夜福利片| videossex国产| 波野结衣二区三区在线| 免费av毛片视频| 国产视频内射| 天堂中文最新版在线下载 | 国产欧美另类精品又又久久亚洲欧美| 国产午夜精品一二区理论片| 久久精品久久精品一区二区三区| 大陆偷拍与自拍| 国产精品一区www在线观看| 亚洲精品久久久久久婷婷小说| 欧美 日韩 精品 国产| 国产精品一及| 亚洲精华国产精华液的使用体验| 久久韩国三级中文字幕| 插逼视频在线观看| 最近的中文字幕免费完整| av免费观看日本| 老师上课跳d突然被开到最大视频| 麻豆国产97在线/欧美| 久久久久久久久久成人| 性色avwww在线观看| 午夜福利高清视频| 久热久热在线精品观看| 丰满人妻一区二区三区视频av| 久久久欧美国产精品| 欧美zozozo另类| 久久鲁丝午夜福利片| 大香蕉97超碰在线| 久久99精品国语久久久| 久久精品国产鲁丝片午夜精品| 国国产精品蜜臀av免费| 日韩在线高清观看一区二区三区| 18禁在线无遮挡免费观看视频| 纵有疾风起免费观看全集完整版 | 一区二区三区高清视频在线| 99久国产av精品| 一级毛片 在线播放| 国产激情偷乱视频一区二区| 69人妻影院| 精品不卡国产一区二区三区| 好男人视频免费观看在线| av网站免费在线观看视频 | 插逼视频在线观看| 麻豆乱淫一区二区| 九九久久精品国产亚洲av麻豆| 大又大粗又爽又黄少妇毛片口| 美女xxoo啪啪120秒动态图| .国产精品久久| 亚洲成人久久爱视频| 美女主播在线视频| 午夜精品一区二区三区免费看| 国产一区二区亚洲精品在线观看| 久久久久国产网址| 精品一区二区三区视频在线| 日韩伦理黄色片| 一级毛片aaaaaa免费看小| 最近中文字幕高清免费大全6| 夜夜看夜夜爽夜夜摸| 国产69精品久久久久777片| 美女黄网站色视频| 人妻一区二区av| 午夜福利视频1000在线观看| 久久99蜜桃精品久久| videos熟女内射| 青春草亚洲视频在线观看| 久久久久久久久久成人| 男女视频在线观看网站免费| 欧美xxⅹ黑人| 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 国产成人一区二区在线| 高清欧美精品videossex| 国产一区二区三区av在线| 最近最新中文字幕大全电影3| 嫩草影院精品99| 99视频精品全部免费 在线| 久久99热6这里只有精品| 久久久久精品久久久久真实原创| 麻豆精品久久久久久蜜桃| 人妻夜夜爽99麻豆av| 在现免费观看毛片| 国产精品久久久久久久久免| 国产成人一区二区在线| 日韩成人av中文字幕在线观看| 国产伦一二天堂av在线观看| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 亚洲自偷自拍三级| 亚洲精品一区蜜桃| 男人舔女人下体高潮全视频| 精品欧美国产一区二区三| 边亲边吃奶的免费视频| 大香蕉97超碰在线| 美女脱内裤让男人舔精品视频| 最近中文字幕高清免费大全6| 国产午夜精品一二区理论片| 亚洲欧洲国产日韩| 少妇人妻精品综合一区二区| 高清午夜精品一区二区三区| av在线播放精品| 久久精品国产亚洲av涩爱| 亚洲av一区综合| 91狼人影院| 国产av国产精品国产| 国产精品99久久久久久久久| 久久6这里有精品| 久久久久久九九精品二区国产| 别揉我奶头 嗯啊视频| 搡老乐熟女国产| 夜夜看夜夜爽夜夜摸| 欧美 日韩 精品 国产| 五月玫瑰六月丁香| 成年版毛片免费区| 日韩一区二区三区影片| 精品午夜福利在线看| 亚洲色图av天堂| 免费在线观看成人毛片| 色吧在线观看| 色5月婷婷丁香| 在现免费观看毛片| 搡女人真爽免费视频火全软件| 成人一区二区视频在线观看| 久久99热这里只频精品6学生| 天天躁夜夜躁狠狠久久av| av黄色大香蕉| 久久99热这里只频精品6学生| 日韩一本色道免费dvd| 国产 一区 欧美 日韩| 国产av国产精品国产| 神马国产精品三级电影在线观看| 国产色爽女视频免费观看| 欧美zozozo另类| 久久99精品国语久久久| 日韩电影二区| 性色avwww在线观看| 少妇被粗大猛烈的视频|