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

    基于峭度與互相關(guān)的IEWT軸承故障診斷方法研究

    2019-03-24 01:25:38郭輝伍川輝劉澤潮黃衍于昊明
    關(guān)鍵詞:尺度空間峭度頻帶

    郭輝,伍川輝,劉澤潮,黃衍,于昊明

    基于峭度與互相關(guān)的IEWT軸承故障診斷方法研究

    郭輝,伍川輝,劉澤潮,黃衍,于昊明

    (西南交通大學(xué) 機(jī)械工程學(xué)院,四川 成都 610031)

    經(jīng)驗(yàn)小波變換(EWT)作為自適應(yīng)的信號(hào)分解方法,其分解結(jié)果與頻帶邊界的劃分密切相關(guān)。尺度空間表征(SSR)算法可以根據(jù)信號(hào)的頻譜自適應(yīng)地實(shí)現(xiàn)信號(hào)的頻帶劃分。但是,在對(duì)軸承信號(hào)處理時(shí)SSR容易產(chǎn)生過度劃分,這就導(dǎo)致一個(gè)共振頻帶被分解成多個(gè)分量,從而無法有效地提取故障沖擊。因此,提出基于峭度準(zhǔn)則與互相關(guān)分析的改進(jìn)EWT算法:首先,根據(jù)尺度空間表征對(duì)信號(hào)的頻譜進(jìn)行劃分;然后,計(jì)算每個(gè)分量的峭度值,對(duì)峭度值大于閾值的分量進(jìn)行相關(guān)性分析,選取相關(guān)性高的分量進(jìn)行合并;最后,對(duì)合并后的分量進(jìn)行包絡(luò)分析,提取故障沖擊。通過使用仿真與實(shí)驗(yàn)數(shù)據(jù)對(duì)所提出的IEWT算法進(jìn)行驗(yàn)證,證明該方法可以有效地對(duì)軸承早期微弱故障進(jìn)行診斷。

    尺度空間;經(jīng)驗(yàn)小波變換;峭度;互相關(guān)分析;故障診斷

    滾動(dòng)軸承作為高速列車機(jī)械傳動(dòng)系統(tǒng)中應(yīng)用最廣泛的關(guān)鍵部件之一,一旦發(fā)生故障,將造成巨大的人員與經(jīng)濟(jì)損失。因此,針對(duì)滾動(dòng)軸承的故障診斷研究具有重大現(xiàn)實(shí)意義。軸承的故障信號(hào)一般是非平穩(wěn)、非線性信號(hào)[1-2],針對(duì)這一類信號(hào),使用小波包變換(WPT)[3]、經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)[4]等方法可以有效地對(duì)軸承早期微弱故障沖擊進(jìn)行提取。WPT在進(jìn)行信號(hào)分解時(shí),其分量個(gè)數(shù)由分解層數(shù)所決定。因此,當(dāng)分解層數(shù)固定時(shí),其各個(gè)分量的頻帶也固定。EMD是一種自適應(yīng)的信號(hào)分解算法,可以根據(jù)信號(hào)的特征自適應(yīng)地實(shí)行信號(hào)的分解。但是,EMD存在模態(tài)混疊、端點(diǎn)效應(yīng)的缺點(diǎn)[5],這些缺點(diǎn)對(duì)故障信息的提取具有嚴(yán)重影響很大。因此,Gilles等[6-7]結(jié)合EMD的自適應(yīng)性和小波變換的理論框架等優(yōu)點(diǎn),提出一種新的信號(hào)處理方法即經(jīng)驗(yàn)小波變換[8-9](empirical wavelet transform)。在此基礎(chǔ)上,李志農(nóng)等[10]將該方法應(yīng)用于轉(zhuǎn)子磨碰故障診斷中,區(qū)分出碰磨故障的嚴(yán)重程度,驗(yàn)證了經(jīng)驗(yàn)小波在轉(zhuǎn)子故障診斷的有效性。向玲等[11?12]通過在頻域中找到2個(gè)連續(xù)極大值之間的極小值作為頻域劃分界限,提出IEWT和FSK聯(lián)合診斷的方法,識(shí)別出齒輪和滾動(dòng)軸承故障頻率。但是這些方法在進(jìn)行頻帶劃分時(shí)仍然存在一個(gè)共振頻帶被劃分為幾個(gè)分量的現(xiàn)象[13],文獻(xiàn)中并沒有提出改善的方法。在此,本文針對(duì)共振頻帶過度劃分的現(xiàn)象,在尺度空間表征的基礎(chǔ)上根據(jù)峭度與互相關(guān)準(zhǔn)則提出改進(jìn)的IEWT 算法。首先,根據(jù)尺度空間表征對(duì)信號(hào)的頻譜進(jìn)行劃分;然后,計(jì)算每個(gè)分量的峭度值,對(duì)分量進(jìn)行篩選;再對(duì)篩選后的分量進(jìn)行互相關(guān)性分析,選取相關(guān)性高的分量進(jìn)行合并;最后,對(duì)合并后的分量進(jìn)行包絡(luò)譜分析,識(shí)別軸承的故障類型。使用IEWT算法對(duì)仿真信號(hào)與實(shí)驗(yàn)信號(hào)進(jìn)行分析,證明IEWT算法可以有效地對(duì)軸承早期微弱故障進(jìn)行診斷。

    1 基本理論

    1.1 EWT原理

    EWT的核心是根據(jù)信號(hào)的頻譜劃分,自適應(yīng)地構(gòu)造正交小波濾波器組以提取包含故障沖擊的 AM-FM 分量。

    信號(hào)()的傅里葉變換為(),頻率范圍為?[0,π]。ω表示頻帶邊界的中心頻率,將[0,π]劃分為個(gè)頻帶,得到+1條邊界(包含[0,π])。頻譜分界示意圖如圖所示,每個(gè)頻帶可以表示為:

    過度區(qū)間可以定義為以ωn為中心,寬度為Tn=2τn的區(qū)間,如圖1陰影區(qū)。

    其中:傅里葉變?yōu)閇?],傅里葉反變換為?1[?],則信號(hào)()的經(jīng)驗(yàn)小波變換的細(xì)節(jié)系數(shù)可以由經(jīng)驗(yàn)小波函數(shù)和信號(hào)內(nèi)積產(chǎn)生:

    ()的經(jīng)驗(yàn)小波變換的近似系數(shù)可以由尺度函數(shù)和信號(hào)內(nèi)積產(chǎn)生:

    對(duì)于一個(gè)時(shí)間序列(),經(jīng)過經(jīng)驗(yàn)小波分解后可以得到

    1.2 尺度空間

    尺度空間方法是使用高斯函數(shù)對(duì)信號(hào)的頻譜進(jìn)行平滑,隨著平滑次數(shù)的增加,一些虛假的頻帶分界點(diǎn)消失,而真實(shí)的分界點(diǎn)保留。所以每平滑1次,求取1次頻帶的分界點(diǎn)并記錄下來。最后生成一張關(guān)于尺度和頻帶分界點(diǎn)的二維圖像,這個(gè)圖像記錄的是每一次平滑后的分界點(diǎn)位置。尺度空間定義為:

    其中:*表示卷積;()表示需要平滑的信號(hào)的頻譜。離散的尺度空間定義:

    其中:為高斯核函數(shù)的長(zhǎng)度,且

    圖2 尺度空間曲線

    1.3 峭度

    峭度是概率密度分布尖峭程度的度量,是對(duì)軸承振動(dòng)信號(hào)進(jìn)行時(shí)域處理最常用的無綱量參考指標(biāo),數(shù)學(xué)描述為:

    式中:為振動(dòng)信號(hào);為信號(hào)的均值;為信號(hào)的標(biāo)準(zhǔn)差。

    峭度在軸承中表現(xiàn)出很強(qiáng)規(guī)律性,軸承在剛投入使用期間,其振動(dòng)信號(hào)幅值分布接近正態(tài)分布,峭度值保持大約等于3。在軸承進(jìn)入使用后期,出現(xiàn)故障的時(shí)候,信號(hào)幅值會(huì)明顯偏離正態(tài)分布,峭度會(huì)隨之增大,此時(shí)可以認(rèn)為出現(xiàn)了故障初期現(xiàn)象,篩選出峭度值超過4的分量進(jìn)行下一步計(jì)算。

    1.4 互相關(guān)系數(shù)

    在信號(hào)處理中,用互相關(guān)來衡量2個(gè)時(shí)間序列()和()在2個(gè)不同時(shí)刻取值之間的相似程度。對(duì)于互相關(guān)函數(shù)定義為:

    互相關(guān)函數(shù)歸一化處理:

    2 故障診斷流程

    因?yàn)椋叨瓤臻g表征在對(duì)軸承故障信號(hào)的頻譜進(jìn)行劃分時(shí),容易存在過度劃分,因此,本文在峭度與互相關(guān)準(zhǔn)則的基礎(chǔ)上提出一種改進(jìn)的EWT分解——IEWT,見圖3。方法診斷流程如下:

    1) 對(duì)原始信號(hào)進(jìn)行傅里葉變換,用尺度空間自適應(yīng)地對(duì)頻帶進(jìn)行劃分。

    2) 計(jì)算各個(gè)分量的峭度值,設(shè)置閾值進(jìn)行篩選。

    3) 篩選出的分量若存在相鄰分量峭度值相近,則對(duì)相鄰分量作互相關(guān)分析。

    4) 若在計(jì)算相關(guān)分析結(jié)果中,相關(guān)系數(shù)超過0.5,則對(duì)相鄰分量進(jìn)行合并,否則不合并。

    5) 對(duì)合并后的分量進(jìn)行包絡(luò)譜分析,找出故障頻率。

    圖3 診斷流程圖

    3 仿真分析

    通過SSR對(duì)仿真信號(hào)進(jìn)行頻帶劃分,頻帶劃分如圖5。分解得到9個(gè)分量,計(jì)算各個(gè)分量的峭度,結(jié)果如圖6所示。c5和c6的峭度值超過閾值,表明c5和c6 2個(gè)分量包含了故障頻率,且分量又相鄰,可能是被過度劃分的共振頻帶,所以對(duì)其做進(jìn)一步的相關(guān)性分析,結(jié)果如表1所示。由表1可知,c5和c6高度相關(guān),說明c5和c6分量確是一個(gè)共振頻帶被過度劃分為2個(gè)子頻帶,所以對(duì)其合并重構(gòu)。最后對(duì)合并分量進(jìn)行包絡(luò)譜分析,結(jié)果如圖7所示。通過對(duì)比圖7~8,可以發(fā)現(xiàn)合并后分量信號(hào)包含更多故障頻率倍頻,進(jìn)一步說明了算法對(duì)共振頻帶過度劃分的現(xiàn)象起到了抑制作用,從而能更有效診斷出軸承早期故障。

    圖4 仿真信號(hào)時(shí)域圖

    圖5 SSR頻帶劃分

    圖6 IEWT分量峭度值計(jì)算結(jié)果(閾值為4)

    表1 IEWT分量互相關(guān)系數(shù)計(jì)算

    圖7 合并分量頻譜分析

    (a) c5分量頻譜分析;(b) c6分量頻譜分析

    4 試驗(yàn)分析

    為了進(jìn)一步驗(yàn)證算法有效性,應(yīng)用在西南交通大學(xué)牽引動(dòng)力國(guó)家重點(diǎn)實(shí)驗(yàn)室試驗(yàn)的滾動(dòng)軸承外圈故障數(shù)據(jù)進(jìn)行驗(yàn)證,并與EWT進(jìn)行對(duì)比,試驗(yàn)設(shè)備如圖9所示。試驗(yàn)采樣頻率為10 kHz,軸旋轉(zhuǎn)頻率為10.28 Hz。由式(15)可知,軸承外圈故障頻率BPFO= 41.65 Hz。

    圖9 試驗(yàn)臺(tái)

    試驗(yàn)信號(hào)時(shí)域波形如圖10所示。根據(jù)SSR得到的頻譜劃分邊界如圖11所示,然后計(jì)算每個(gè)分量包絡(luò)的峭度值,結(jié)果如圖12所示。由圖12可知,c2,c3和c4峭度值超過閾值,表明這三者都包含了軸承故障頻率,又因?yàn)楸舜讼噜彛赡苁枪舱耦l帶過度劃分導(dǎo)致,所以對(duì)其進(jìn)行互相關(guān)性分析,結(jié)果如表2所示。根據(jù)表2所示,c3和c4具有高度的相關(guān)性,說明c3和c4是一個(gè)共振頻帶被過度劃分成的2個(gè)分量。因此,對(duì)這2個(gè)分量進(jìn)行合并。合并后分量的包絡(luò)譜如圖13所示,圖14為合并前c3和c4的包絡(luò)譜。通過對(duì)比表明,圖13中故障特征頻率及其諧波更加明顯,進(jìn)一步證明了該算法對(duì)共振頻帶被過度劃分的現(xiàn)象起到了抑制作用,IEWT可以更加有效地實(shí)現(xiàn)軸承早期微弱故障的診斷。

    圖10 原始數(shù)據(jù)時(shí)域圖

    圖11 SSR頻帶劃分

    圖12 IEWT分量峭度值計(jì)算結(jié)果(閾值為4)

    表2 IEWT分量互相關(guān)系數(shù)計(jì)算

    圖13 合并分量包絡(luò)譜分析

    (a) c3分量頻譜分析;(b) c4分量頻譜分析

    5 結(jié)論

    1) 尺度空間表征可以自適應(yīng)地實(shí)現(xiàn)信號(hào)頻帶的劃分。

    2) 在對(duì)軸承故障信號(hào)進(jìn)行頻帶劃分時(shí),尺度空間表征并不能得到最優(yōu)的帶寬。

    3) 基于峭度和互相關(guān)準(zhǔn)則的IEWT可以得到包含故障沖擊更多的頻帶。

    [1] 向玲, 唐貴基, 胡愛軍. 旋轉(zhuǎn)機(jī)械非平穩(wěn)振動(dòng)信號(hào)的時(shí)頻分析比較[J]. 振動(dòng)與沖擊, 2010, 29(2): 42?45. XIANG Ling, TANG Guiji, HU Aijun. Time frequency analysis of nonstationary vibration signals in rotating machinery[J]. Journal of Vibration and Shock, 2010, 29(2): 42?45.

    [2] Yeh Poliang, LIU Peiling. Application of the wavelet transform and the enhanced Fourier spectrum in the impact echo test[J]. NDT and E International, 2008, 41(5): 382?394.

    [3] Jivet I, Tomoroga M. Ultrasound image enhancement for visualization using spatially variant morphological filters[M]. Springer International Publishing, 2014: 255? 258.

    [4] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis[J]. Proceedings of the Royal Society of London Series A, 1998, 454: 903?995.

    [5] 鐘佑明, 秦樹人, 湯寶平. 一種振動(dòng)信號(hào)新變換法的研究[J]. 振動(dòng)工程學(xué)報(bào), 2002, 15(2): 113?118. ZHONG Youming, QIN Shuren, TANG Baoping. A new method of vibration signal transformation[J]. Journal of Vibration Engineering, 2002, 15(2): 113?118.

    [6] Gilles J. Empirical wavelet transform[J]. IEEE Transactions on Signal Processing, 2013, 61(16): 3999? 4010.

    [7] Gilles J, Tran G, Osher S. 2D empirical transforms. Wavelets, ridgelets and curvelets revisited[J]. Siam Journal on Imaging Sciences, 2014, 7(7): 157?186.

    [8] GU X, YANG S, LIU Y, et al. Rolling element bearing faults diagnosis based on kurtogram and frequency domain correlated kurtosis[J]. Measurement Science and Technology, 2016, 27(12): 125?129.

    [9] GU X, YANG S, LIU Y, et al. Compound faults detection of the rolling element bearing based on the optimal complex Morlet wavelet filter[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2018, 232(10): 1786? 1801.

    [10] 李志農(nóng), 朱明, 褚福磊, 等. 基于經(jīng)驗(yàn)小波變換的機(jī)械故障診斷方法研究[J]. 儀器儀表學(xué)報(bào), 2014, 35(11): 2423?2432. LI Zhinong, ZHU Ming, CHU Fulei, et al. Research on Mechanical Fault Diagnosis Based on Empirical Wavelet Transform[J]. Journal of Instrumentation, 2014, 35(11): 2423?2432.

    [11] 向玲, 高雪媛, 張力佳, 等. IEWT和FSK在齒輪和滾動(dòng)軸承故障診斷中的應(yīng)用[J]. 振動(dòng)、測(cè)試與診斷, 2017, 37(6): 1256?1261. XIANG Ling, GAO Xueyuan, ZHANG Lijia, et al. Application of IEWT and FSK in fault diagnosis of gears and rolling bearings[J]. Vibration, Testing and Diagnosis, 2017, 37(6): 1256?1261.

    [12] 李媛媛. 基于經(jīng)驗(yàn)小波變換的機(jī)械故障診斷方法[D]. 北京: 華北電力大學(xué), 2016: 1?50. LI Yuanyuan. Mechanical fault diagnosis method based on empirical wavelet transform[D]. Beijing: North China Electric Power University, 2016: 1?50.

    [13] Antoni J, Randall R B. The spectral kurtosis: application to the vibratory surveillance and diagnostics of rotating machines[J]. Systems and Signal Processing, 2006, 20(2): 308?331.

    [14] LIANG Ming, Soltani Bozchalooi I. An energy operator approach to joint application of amplitude and frequency- demodulations for bearing fault detection[J]. Mechanical Systems and Signal Processing, 2010, 24(1): 1473?1494.

    Research on fault diagnosis method of IEWT bearing based on kurtosis and cross correlation

    GUO Hui, WU Chuanhui, LIU Zechao, HUANG Yan, YU Haoming

    (School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China)

    Empirical wavelet transform (EWT) was an adaptive signal decomposition method, and its decomposition result was closely related to the division of frequency band boundaries. The Scale Space Representation (SSR) algorithm can adaptively implement the frequency band division of the signal according to the spectrum of the signal. However, the SSR is prone to transitional division when processing the bearing signals, which causes one resonance frequency band to be decomposed into a plurality of components, so that the fault impact cannot be effectively extracted. Therefore, this paper proposes an Improved EWT algorithm based on kurtosis criteria and cross-correlation analysis. First, the spectrum of the signal was divided according to the scale space representation. Then, the kurtosis value of each component was calculated, and the component with the kurtosis value greater than the threshold was calculated. Correlation analysis was performed, the components with high correlation for merging was selected; Finally, perform envelope analysis on the combined components to extract the fault impact. The proposed IEWT algorithm was verified by using simulation and experimental data, which proves that the method can effectively diagnose the early weak faults of the bearing.

    scale space; empirical wavelet transform; kurtosis; cross correlation analysis; fault diagnosis

    TH133.33;TN911.72

    A

    1672 ? 7029(2019)07? 1774 ? 07

    10.19713/j.cnki.43?1423/u.2019.07.023

    2018?09?25

    國(guó)家重點(diǎn)研發(fā)計(jì)劃先進(jìn)軌道交通重點(diǎn)專項(xiàng)資助項(xiàng)目(2017YFB1201004)

    伍川輝(1964?),男,四川成都人,副教授,從事先進(jìn)測(cè)控技術(shù)及應(yīng)用研究;E?mail:chwu126@126.com

    (編輯 陽麗霞)

    猜你喜歡
    尺度空間峭度頻帶
    基于MCKD和峭度的液壓泵故障特征提取
    基于AHP的大尺度空間域礦山地質(zhì)環(huán)境評(píng)價(jià)研究
    Wi-Fi網(wǎng)絡(luò)中5G和2.4G是什么?有何區(qū)別?
    聯(lián)合快速峭度圖與變帶寬包絡(luò)譜峭度圖的輪對(duì)軸承復(fù)合故障檢測(cè)研究
    單音及部分頻帶干擾下DSSS系統(tǒng)性能分析
    基于峭度分析的聲發(fā)射故障檢測(cè)
    電子世界(2018年12期)2018-07-04 06:34:38
    雙頻帶隔板極化器
    居住區(qū)園林空間尺度研究
    基于魯棒性小波包峭度圖的滾動(dòng)軸承故障診斷*
    基于降采樣歸一化割的多尺度分層分割方法研究
    一级av片app| 九九在线视频观看精品| av在线观看视频网站免费| 日韩欧美三级三区| 69人妻影院| 亚洲av男天堂| 一夜夜www| 国产真实乱freesex| 欧洲精品卡2卡3卡4卡5卡区| h日本视频在线播放| 亚洲欧美精品专区久久| a级毛片a级免费在线| 国产亚洲精品久久久com| 日本三级黄在线观看| 熟妇人妻久久中文字幕3abv| 亚洲成人精品中文字幕电影| 国产一区二区亚洲精品在线观看| 高清毛片免费观看视频网站| 男女那种视频在线观看| 亚洲欧美精品自产自拍| 中文精品一卡2卡3卡4更新| 国产精品久久视频播放| 乱码一卡2卡4卡精品| 少妇高潮的动态图| 久久这里只有精品中国| 性色avwww在线观看| 欧美另类亚洲清纯唯美| 国产免费男女视频| 久久人人精品亚洲av| 丰满乱子伦码专区| 青春草国产在线视频 | 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| а√天堂www在线а√下载| 狂野欧美白嫩少妇大欣赏| 中文字幕久久专区| 两个人的视频大全免费| 国产一区亚洲一区在线观看| 亚洲最大成人手机在线| 国产精品电影一区二区三区| 2022亚洲国产成人精品| 午夜亚洲福利在线播放| 99热这里只有精品一区| 全区人妻精品视频| 国产国拍精品亚洲av在线观看| 午夜福利视频1000在线观看| 午夜精品一区二区三区免费看| 日本黄色视频三级网站网址| 色哟哟·www| 国产老妇伦熟女老妇高清| 精品人妻一区二区三区麻豆| 永久网站在线| 一边亲一边摸免费视频| 久久精品久久久久久久性| 亚洲欧美成人精品一区二区| 国产午夜精品久久久久久一区二区三区| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看 | 在线a可以看的网站| 在线观看一区二区三区| 毛片一级片免费看久久久久| 又黄又爽又刺激的免费视频.| 国产亚洲精品av在线| 亚洲精品乱码久久久久久按摩| 一级毛片电影观看 | 麻豆av噜噜一区二区三区| 18+在线观看网站| 日韩人妻高清精品专区| 成年女人永久免费观看视频| 麻豆乱淫一区二区| 午夜免费男女啪啪视频观看| 国产三级在线视频| 一个人免费在线观看电影| 又爽又黄无遮挡网站| 亚洲国产日韩欧美精品在线观看| 久久欧美精品欧美久久欧美| 国产视频首页在线观看| 国产单亲对白刺激| 看免费成人av毛片| 美女国产视频在线观看| 久久久色成人| 看黄色毛片网站| 国内少妇人妻偷人精品xxx网站| 国产一区二区在线av高清观看| 国产美女午夜福利| 国产精品不卡视频一区二区| 97人妻精品一区二区三区麻豆| 青青草视频在线视频观看| 久久久国产成人精品二区| 久久久久国产网址| 亚州av有码| 亚洲欧美日韩高清专用| 国产成人午夜福利电影在线观看| 亚洲激情五月婷婷啪啪| 69人妻影院| 一卡2卡三卡四卡精品乱码亚洲| 免费无遮挡裸体视频| 菩萨蛮人人尽说江南好唐韦庄 | 麻豆成人午夜福利视频| 少妇的逼水好多| 免费看a级黄色片| 欧美性猛交╳xxx乱大交人| 国产亚洲欧美98| 日韩在线高清观看一区二区三区| avwww免费| 91久久精品国产一区二区成人| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 亚洲三级黄色毛片| 亚洲欧美精品自产自拍| 免费观看人在逋| 欧美另类亚洲清纯唯美| 三级国产精品欧美在线观看| 久久综合国产亚洲精品| 高清在线视频一区二区三区 | 波多野结衣高清无吗| 九九热线精品视视频播放| 久久精品久久久久久久性| 亚洲在久久综合| 日本五十路高清| 国国产精品蜜臀av免费| 午夜免费激情av| 久久久国产成人精品二区| 精品久久久久久久久久免费视频| 中文字幕久久专区| 国产亚洲av片在线观看秒播厂 | 3wmmmm亚洲av在线观看| 国产高清激情床上av| 91久久精品国产一区二区三区| 免费搜索国产男女视频| 日本一二三区视频观看| 成人无遮挡网站| 日韩 亚洲 欧美在线| 中国美女看黄片| 一边亲一边摸免费视频| 校园人妻丝袜中文字幕| 久久久久网色| 日韩,欧美,国产一区二区三区 | 亚洲va在线va天堂va国产| 青春草视频在线免费观看| 男人舔奶头视频| 观看美女的网站| 一级av片app| 国产白丝娇喘喷水9色精品| 少妇被粗大猛烈的视频| 嫩草影院新地址| 中文字幕精品亚洲无线码一区| 亚洲国产精品久久男人天堂| 麻豆av噜噜一区二区三区| 国产欧美日韩精品一区二区| 岛国毛片在线播放| 成人特级黄色片久久久久久久| 成人特级黄色片久久久久久久| 欧美+亚洲+日韩+国产| 国产精品福利在线免费观看| 国产蜜桃级精品一区二区三区| 亚洲av男天堂| 国产一区二区亚洲精品在线观看| 亚洲精品久久国产高清桃花| 两个人视频免费观看高清| 能在线免费看毛片的网站| 禁无遮挡网站| 欧美色欧美亚洲另类二区| 欧美日韩国产亚洲二区| 亚洲精品色激情综合| 久久精品夜色国产| 久久鲁丝午夜福利片| 高清日韩中文字幕在线| 毛片一级片免费看久久久久| 91精品一卡2卡3卡4卡| 男人舔女人下体高潮全视频| 三级毛片av免费| 久久鲁丝午夜福利片| 亚州av有码| 免费在线观看成人毛片| 亚洲中文字幕日韩| 两个人视频免费观看高清| 国产在线精品亚洲第一网站| 我的老师免费观看完整版| 国产午夜精品久久久久久一区二区三区| 国产成人一区二区在线| 亚洲精品久久久久久婷婷小说 | 日韩欧美在线乱码| 乱码一卡2卡4卡精品| 人人妻人人澡人人爽人人夜夜 | 91久久精品国产一区二区成人| 国产91av在线免费观看| 婷婷亚洲欧美| 神马国产精品三级电影在线观看| 日韩强制内射视频| 国产精品嫩草影院av在线观看| 国产色婷婷99| 欧美一区二区精品小视频在线| 有码 亚洲区| 乱系列少妇在线播放| 别揉我奶头 嗯啊视频| 亚洲最大成人中文| 91aial.com中文字幕在线观看| 青春草亚洲视频在线观看| 蜜桃亚洲精品一区二区三区| 久久久久久久久中文| 免费观看在线日韩| 国产精品国产三级国产av玫瑰| 色播亚洲综合网| 国内久久婷婷六月综合欲色啪| 女同久久另类99精品国产91| 亚洲av.av天堂| 日韩欧美精品免费久久| 欧洲精品卡2卡3卡4卡5卡区| 天堂中文最新版在线下载 | 在线观看免费视频日本深夜| 最好的美女福利视频网| 嫩草影院精品99| 久久精品久久久久久久性| 岛国毛片在线播放| 男女做爰动态图高潮gif福利片| 成人无遮挡网站| 十八禁国产超污无遮挡网站| 寂寞人妻少妇视频99o| 九色成人免费人妻av| 精品无人区乱码1区二区| 色5月婷婷丁香| 少妇被粗大猛烈的视频| 成人无遮挡网站| 亚洲av中文字字幕乱码综合| 日本与韩国留学比较| 亚洲欧美日韩卡通动漫| 国产又黄又爽又无遮挡在线| 22中文网久久字幕| 观看美女的网站| 又爽又黄a免费视频| 日本-黄色视频高清免费观看| 99热网站在线观看| 成人性生交大片免费视频hd| 久久久久免费精品人妻一区二区| 黄片无遮挡物在线观看| 欧美日韩乱码在线| 99热这里只有精品一区| 人人妻人人澡人人爽人人夜夜 | 一本久久中文字幕| 国产精品电影一区二区三区| 22中文网久久字幕| 色哟哟·www| 亚洲成a人片在线一区二区| 国产精品人妻久久久久久| 午夜福利在线在线| 欧美性猛交黑人性爽| 亚洲欧美中文字幕日韩二区| 99久久人妻综合| 亚洲一级一片aⅴ在线观看| 听说在线观看完整版免费高清| 最近2019中文字幕mv第一页| 久久人妻av系列| 国产极品精品免费视频能看的| 丰满的人妻完整版| 亚洲av熟女| 婷婷色av中文字幕| kizo精华| 亚洲自拍偷在线| 综合色av麻豆| 中文字幕精品亚洲无线码一区| 伦精品一区二区三区| 日日啪夜夜撸| 国产精品女同一区二区软件| 精品国内亚洲2022精品成人| 久久精品国产清高在天天线| 亚洲美女搞黄在线观看| 五月伊人婷婷丁香| 久久人人爽人人片av| 寂寞人妻少妇视频99o| 国产成人a∨麻豆精品| 欧美高清成人免费视频www| 久久久久久久久久黄片| 国产av一区在线观看免费| 精品人妻视频免费看| 日韩 亚洲 欧美在线| 国产精品1区2区在线观看.| 有码 亚洲区| 成人二区视频| 天天一区二区日本电影三级| 欧美日韩精品成人综合77777| 91久久精品国产一区二区三区| 欧美3d第一页| 色5月婷婷丁香| 国产三级在线视频| 综合色丁香网| 中文字幕免费在线视频6| 久久鲁丝午夜福利片| 婷婷色综合大香蕉| 免费搜索国产男女视频| 日本三级黄在线观看| .国产精品久久| 菩萨蛮人人尽说江南好唐韦庄 | 精品久久国产蜜桃| 国产成人aa在线观看| 国产伦一二天堂av在线观看| 国产成人91sexporn| 成人二区视频| 少妇人妻一区二区三区视频| 婷婷色av中文字幕| 久久欧美精品欧美久久欧美| 午夜激情欧美在线| 此物有八面人人有两片| 热99re8久久精品国产| av在线蜜桃| 成人无遮挡网站| 99久久九九国产精品国产免费| 免费大片18禁| 伦精品一区二区三区| 综合色av麻豆| 国产欧美日韩精品一区二区| 国产精品三级大全| 国产精品电影一区二区三区| 国产成人福利小说| 成人一区二区视频在线观看| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 国产精品人妻久久久影院| 国产精品永久免费网站| 久久亚洲精品不卡| 99热这里只有精品一区| 亚洲欧洲国产日韩| 国产精品久久电影中文字幕| 久久亚洲精品不卡| 久久精品国产鲁丝片午夜精品| 色噜噜av男人的天堂激情| 中文字幕精品亚洲无线码一区| 变态另类成人亚洲欧美熟女| 中文字幕熟女人妻在线| 免费看av在线观看网站| 欧美激情国产日韩精品一区| 中文欧美无线码| 少妇猛男粗大的猛烈进出视频 | 免费av观看视频| 亚洲真实伦在线观看| 观看美女的网站| 亚洲av电影不卡..在线观看| 伦理电影大哥的女人| 久久久成人免费电影| 高清毛片免费观看视频网站| 国产在视频线在精品| 久久精品91蜜桃| 97超视频在线观看视频| 美女 人体艺术 gogo| 国产白丝娇喘喷水9色精品| 一边亲一边摸免费视频| 99久国产av精品| 寂寞人妻少妇视频99o| 日本免费a在线| 老师上课跳d突然被开到最大视频| 午夜福利高清视频| 夜夜爽天天搞| 成人一区二区视频在线观看| 日韩中字成人| 久久久久网色| 国产精品乱码一区二三区的特点| 国产av不卡久久| 成年av动漫网址| 欧美日韩精品成人综合77777| 一级av片app| 亚洲欧美日韩无卡精品| 99久国产av精品| 亚洲成人中文字幕在线播放| 亚洲av电影不卡..在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 免费人成在线观看视频色| 精品午夜福利在线看| 亚洲在线观看片| 亚洲熟妇中文字幕五十中出| 国产av麻豆久久久久久久| 国产成人精品婷婷| 国产av麻豆久久久久久久| 精品久久国产蜜桃| 免费av观看视频| av在线天堂中文字幕| 亚洲成av人片在线播放无| 综合色av麻豆| 内射极品少妇av片p| 久久综合国产亚洲精品| 久久久久国产网址| 一区二区三区免费毛片| 国产一级毛片七仙女欲春2| 国产精品嫩草影院av在线观看| 精品一区二区三区人妻视频| 狂野欧美白嫩少妇大欣赏| 久久欧美精品欧美久久欧美| 性欧美人与动物交配| 日日啪夜夜撸| 91精品国产九色| 日韩中字成人| 国产av一区在线观看免费| 精品一区二区免费观看| 男女那种视频在线观看| 少妇的逼好多水| 久久这里有精品视频免费| 爱豆传媒免费全集在线观看| 十八禁国产超污无遮挡网站| 一边亲一边摸免费视频| 在线国产一区二区在线| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 欧美在线一区亚洲| 免费av不卡在线播放| 亚洲精品久久国产高清桃花| 永久网站在线| 欧美xxxx性猛交bbbb| 日韩av不卡免费在线播放| 99热6这里只有精品| 久久久午夜欧美精品| 午夜激情福利司机影院| 少妇人妻一区二区三区视频| 内地一区二区视频在线| 蜜臀久久99精品久久宅男| 最近的中文字幕免费完整| 亚洲av中文av极速乱| 啦啦啦啦在线视频资源| 亚洲欧美成人综合另类久久久 | 成人特级黄色片久久久久久久| 国产精品嫩草影院av在线观看| 91久久精品电影网| 久久99热6这里只有精品| 久久国产乱子免费精品| 国产av一区在线观看免费| 日韩三级伦理在线观看| 少妇的逼水好多| 免费看美女性在线毛片视频| 老熟妇乱子伦视频在线观看| 男插女下体视频免费在线播放| 成人毛片60女人毛片免费| 老司机影院成人| 天天躁夜夜躁狠狠久久av| 久久久久性生活片| eeuss影院久久| 99久久成人亚洲精品观看| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 欧美性感艳星| 黄色日韩在线| av女优亚洲男人天堂| 亚洲精品久久久久久婷婷小说 | 亚洲在线自拍视频| 日韩强制内射视频| www.av在线官网国产| 亚洲精品久久国产高清桃花| 午夜福利高清视频| 一级av片app| 国产蜜桃级精品一区二区三区| 国产日韩欧美在线精品| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 国模一区二区三区四区视频| 欧美一区二区精品小视频在线| 夜夜夜夜夜久久久久| 久久精品国产自在天天线| 老师上课跳d突然被开到最大视频| 最近的中文字幕免费完整| 三级经典国产精品| 亚洲精品国产成人久久av| 国产成人精品一,二区 | 婷婷色av中文字幕| 国产精品综合久久久久久久免费| 一个人观看的视频www高清免费观看| 在线播放国产精品三级| 18禁黄网站禁片免费观看直播| 在线免费观看不下载黄p国产| 日本-黄色视频高清免费观看| 亚洲五月天丁香| 人妻制服诱惑在线中文字幕| 欧美性猛交黑人性爽| 久久久a久久爽久久v久久| 男女做爰动态图高潮gif福利片| 性欧美人与动物交配| 国产精品久久久久久亚洲av鲁大| 亚洲四区av| 亚洲在久久综合| 免费大片18禁| 亚洲av中文字字幕乱码综合| 久久久久久久久久黄片| 亚洲天堂国产精品一区在线| 午夜a级毛片| 国产v大片淫在线免费观看| 欧美日本亚洲视频在线播放| 久久国内精品自在自线图片| 色综合色国产| 久久综合国产亚洲精品| 赤兔流量卡办理| 亚洲内射少妇av| 精品久久久噜噜| 美女国产视频在线观看| 中文字幕免费在线视频6| 3wmmmm亚洲av在线观看| 亚洲精品乱码久久久v下载方式| 在线免费观看的www视频| 在线国产一区二区在线| 六月丁香七月| 久久精品人妻少妇| 欧洲精品卡2卡3卡4卡5卡区| 麻豆国产97在线/欧美| 国产免费男女视频| 免费观看的影片在线观看| 日韩精品青青久久久久久| 国产91av在线免费观看| 少妇的逼好多水| 免费看光身美女| 一本久久精品| 2021天堂中文幕一二区在线观| 免费黄网站久久成人精品| 哪里可以看免费的av片| 九草在线视频观看| 亚洲国产色片| 国产 一区精品| 国产色爽女视频免费观看| 九九热线精品视视频播放| 99久久九九国产精品国产免费| 欧美色视频一区免费| 日韩欧美精品v在线| 精品人妻一区二区三区麻豆| 中文字幕熟女人妻在线| 在线观看免费视频日本深夜| 麻豆精品久久久久久蜜桃| 九九久久精品国产亚洲av麻豆| 99热精品在线国产| 69人妻影院| 身体一侧抽搐| 国产精品久久久久久久久免| 黑人高潮一二区| 亚洲,欧美,日韩| 国产亚洲精品久久久com| 村上凉子中文字幕在线| 最近最新中文字幕大全电影3| 中国国产av一级| 小说图片视频综合网站| 亚洲成a人片在线一区二区| 久久精品国产清高在天天线| av在线亚洲专区| 日本一二三区视频观看| 国产真实伦视频高清在线观看| 观看免费一级毛片| 国产精品精品国产色婷婷| 插逼视频在线观看| 欧美3d第一页| 亚洲18禁久久av| 亚洲最大成人中文| 人妻少妇偷人精品九色| 久久精品国产自在天天线| 可以在线观看毛片的网站| 欧美在线一区亚洲| 久久精品人妻少妇| 精品人妻一区二区三区麻豆| 偷拍熟女少妇极品色| 成年女人看的毛片在线观看| 免费看av在线观看网站| 波多野结衣高清作品| 内射极品少妇av片p| 好男人在线观看高清免费视频| 久久国产乱子免费精品| 国产探花在线观看一区二区| 亚洲七黄色美女视频| 国产精品综合久久久久久久免费| 亚洲不卡免费看| 噜噜噜噜噜久久久久久91| 国产精品国产高清国产av| 在线观看午夜福利视频| 三级国产精品欧美在线观看| 欧美成人免费av一区二区三区| 哪里可以看免费的av片| 日日摸夜夜添夜夜添av毛片| 精品免费久久久久久久清纯| 九九在线视频观看精品| 国产精品伦人一区二区| 亚洲内射少妇av| 黄色日韩在线| 成年免费大片在线观看| 亚洲av成人av| 亚洲色图av天堂| 午夜激情福利司机影院| 日韩欧美三级三区| 看免费成人av毛片| 午夜激情欧美在线| 爱豆传媒免费全集在线观看| 黄色配什么色好看| 一进一出抽搐动态| 91麻豆精品激情在线观看国产| 久久久久久久久久成人| 成人av在线播放网站| 久久久成人免费电影| 成人永久免费在线观看视频| 在现免费观看毛片| 国产高清激情床上av| 尾随美女入室| 99久国产av精品国产电影| 男人舔女人下体高潮全视频| 日韩欧美在线乱码| 亚洲av免费高清在线观看| 成人国产麻豆网| kizo精华| 1024手机看黄色片| 亚洲美女视频黄频| 国产伦一二天堂av在线观看| 亚洲国产精品成人久久小说 | 精品国内亚洲2022精品成人| 嫩草影院入口| 国产真实伦视频高清在线观看| 国产av麻豆久久久久久久| 一级毛片我不卡| 欧美日韩精品成人综合77777| 亚洲人成网站在线观看播放| 99久久精品国产国产毛片| 在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看| 精品欧美国产一区二区三| 国产高清激情床上av| 国产伦在线观看视频一区| 人妻少妇偷人精品九色| 日本黄色视频三级网站网址|