• <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)軸承故障診斷*
    基于降采樣歸一化割的多尺度分層分割方法研究
    成人美女网站在线观看视频| 日本爱情动作片www.在线观看| 高清欧美精品videossex| 亚洲精品第二区| 男女无遮挡免费网站观看| 国产黄片美女视频| 亚洲av成人精品一二三区| 久久精品久久精品一区二区三区| 亚洲人成网站在线播| 另类亚洲欧美激情| 男人舔奶头视频| 又粗又硬又长又爽又黄的视频| 亚洲精品国产av成人精品| 国产亚洲一区二区精品| av专区在线播放| 日本-黄色视频高清免费观看| 国产免费福利视频在线观看| 18禁裸乳无遮挡动漫免费视频| 日韩三级伦理在线观看| 亚洲婷婷狠狠爱综合网| 女的被弄到高潮叫床怎么办| 三级经典国产精品| 日韩制服骚丝袜av| 午夜福利视频精品| 如何舔出高潮| 亚洲成色77777| 日本午夜av视频| 偷拍熟女少妇极品色| 最近的中文字幕免费完整| tube8黄色片| 99精国产麻豆久久婷婷| 各种免费的搞黄视频| videossex国产| 中文字幕制服av| 成人高潮视频无遮挡免费网站| 自拍偷自拍亚洲精品老妇| 久久久久人妻精品一区果冻| 啦啦啦啦在线视频资源| 久久毛片免费看一区二区三区| 两个人的视频大全免费| 多毛熟女@视频| 成年人午夜在线观看视频| 国产免费又黄又爽又色| av网站免费在线观看视频| 久久人人爽人人爽人人片va| 麻豆乱淫一区二区| 99国产精品免费福利视频| 欧美精品一区二区免费开放| 美女福利国产在线 | 欧美三级亚洲精品| 一级毛片 在线播放| 男男h啪啪无遮挡| 在线亚洲精品国产二区图片欧美 | 国产免费又黄又爽又色| 这个男人来自地球电影免费观看 | 亚洲精品色激情综合| 最后的刺客免费高清国语| 久久久久久久久久久丰满| 看十八女毛片水多多多| 99久久中文字幕三级久久日本| 美女内射精品一级片tv| 国产色爽女视频免费观看| 国产有黄有色有爽视频| 国产国拍精品亚洲av在线观看| 男人和女人高潮做爰伦理| 99热这里只有是精品50| 日本黄大片高清| 高清av免费在线| 亚洲国产精品成人久久小说| 久久av网站| 日本与韩国留学比较| 欧美性感艳星| 欧美成人午夜免费资源| 高清午夜精品一区二区三区| 激情五月婷婷亚洲| 韩国av在线不卡| 久久久久久久亚洲中文字幕| 蜜桃在线观看..| 免费观看av网站的网址| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 国产日韩欧美在线精品| 亚洲精品aⅴ在线观看| 国产国拍精品亚洲av在线观看| 亚洲色图综合在线观看| 午夜激情福利司机影院| 亚洲在久久综合| 亚洲精品日韩av片在线观看| 精品久久久久久久末码| 黄色欧美视频在线观看| 这个男人来自地球电影免费观看 | 欧美少妇被猛烈插入视频| 欧美精品亚洲一区二区| 在线精品无人区一区二区三 | 免费少妇av软件| 成人毛片a级毛片在线播放| 大片电影免费在线观看免费| 建设人人有责人人尽责人人享有的 | 久久久久久久久久久免费av| 日韩大片免费观看网站| 国产美女午夜福利| 亚洲精品国产av蜜桃| 欧美日本视频| 久久久久精品久久久久真实原创| 夜夜骑夜夜射夜夜干| 精品亚洲成国产av| 久久久亚洲精品成人影院| 亚洲国产欧美在线一区| 亚洲婷婷狠狠爱综合网| 日本黄色片子视频| 又黄又爽又刺激的免费视频.| 尾随美女入室| 国产在线视频一区二区| 国产免费福利视频在线观看| 22中文网久久字幕| 国产精品秋霞免费鲁丝片| 国产精品一区二区在线不卡| 熟妇人妻不卡中文字幕| 国产av精品麻豆| 亚洲av福利一区| a 毛片基地| 18+在线观看网站| 在线观看免费日韩欧美大片 | 有码 亚洲区| 秋霞在线观看毛片| av福利片在线观看| 97在线视频观看| 免费看av在线观看网站| 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 国产精品久久久久成人av| 亚洲丝袜综合中文字幕| 婷婷色综合大香蕉| 最近手机中文字幕大全| 亚洲欧美成人精品一区二区| 黄片wwwwww| 性色avwww在线观看| 成年女人在线观看亚洲视频| 久久精品久久久久久噜噜老黄| 亚洲精品国产av成人精品| 国产伦精品一区二区三区四那| 欧美精品人与动牲交sv欧美| 久久韩国三级中文字幕| 国产深夜福利视频在线观看| 少妇 在线观看| 一个人看视频在线观看www免费| av免费观看日本| 日韩欧美精品免费久久| 99久久人妻综合| 日韩欧美一区视频在线观看 | 国模一区二区三区四区视频| 亚洲四区av| 激情 狠狠 欧美| 韩国av在线不卡| 日韩电影二区| 91午夜精品亚洲一区二区三区| 丰满乱子伦码专区| 秋霞在线观看毛片| 久热这里只有精品99| 精品国产一区二区三区久久久樱花 | 看免费成人av毛片| 国产黄频视频在线观看| 色吧在线观看| 午夜福利视频精品| 国产深夜福利视频在线观看| 欧美xxⅹ黑人| 免费久久久久久久精品成人欧美视频 | 亚洲色图综合在线观看| 妹子高潮喷水视频| 国产成人一区二区在线| 亚洲精品日韩av片在线观看| 免费观看无遮挡的男女| 久久久久久久国产电影| 亚洲久久久国产精品| 2022亚洲国产成人精品| 日韩中字成人| 日韩制服骚丝袜av| 久久人人爽人人爽人人片va| 91精品伊人久久大香线蕉| 麻豆成人午夜福利视频| 国产精品久久久久久久久免| 日日啪夜夜爽| 久久婷婷青草| 亚洲最大成人中文| 成人亚洲精品一区在线观看 | 我要看日韩黄色一级片| 欧美成人a在线观看| 国产亚洲精品久久久com| 麻豆成人午夜福利视频| 极品少妇高潮喷水抽搐| 亚洲精华国产精华液的使用体验| 人人妻人人爽人人添夜夜欢视频 | 99热这里只有是精品50| 成人午夜精彩视频在线观看| 国产精品99久久久久久久久| 三级经典国产精品| 欧美xxⅹ黑人| 日本与韩国留学比较| 欧美bdsm另类| av又黄又爽大尺度在线免费看| 夫妻午夜视频| av国产精品久久久久影院| 国产在视频线精品| 最近的中文字幕免费完整| 午夜精品国产一区二区电影| 日日啪夜夜爽| 亚洲人成网站在线观看播放| 久久精品国产亚洲av涩爱| 亚洲国产高清在线一区二区三| 日韩av不卡免费在线播放| 少妇被粗大猛烈的视频| 国产一区亚洲一区在线观看| av一本久久久久| 午夜福利视频精品| 丝袜脚勾引网站| 欧美成人精品欧美一级黄| 欧美成人一区二区免费高清观看| 日韩亚洲欧美综合| 热re99久久精品国产66热6| 久久精品国产a三级三级三级| 色综合色国产| 国产精品人妻久久久影院| 精品人妻熟女av久视频| 在线播放无遮挡| 久久精品久久久久久噜噜老黄| 一区二区三区乱码不卡18| 国产日韩欧美在线精品| 女人久久www免费人成看片| 在线观看av片永久免费下载| 最新中文字幕久久久久| 一边亲一边摸免费视频| 国产黄片视频在线免费观看| 婷婷色综合www| 久久精品久久久久久久性| 国产一区有黄有色的免费视频| 国产精品.久久久| 99久久精品热视频| 国产精品熟女久久久久浪| 国产黄色视频一区二区在线观看| 韩国av在线不卡| 国产v大片淫在线免费观看| 国产深夜福利视频在线观看| 少妇精品久久久久久久| 永久网站在线| 有码 亚洲区| 日韩亚洲欧美综合| 久久精品熟女亚洲av麻豆精品| 亚洲精品,欧美精品| 欧美xxxx黑人xx丫x性爽| 一级av片app| 国产精品.久久久| 国产亚洲91精品色在线| av不卡在线播放| 久久久久久久亚洲中文字幕| 永久网站在线| 五月玫瑰六月丁香| 亚洲欧美中文字幕日韩二区| 成人午夜精彩视频在线观看| 亚洲精品乱久久久久久| 中文字幕制服av| 亚洲精品乱码久久久久久按摩| 久久午夜福利片| 欧美人与善性xxx| 麻豆国产97在线/欧美| 欧美成人精品欧美一级黄| 成人毛片a级毛片在线播放| videossex国产| 内射极品少妇av片p| 国产亚洲欧美精品永久| 在线观看人妻少妇| 亚洲综合色惰| 国产在线一区二区三区精| 久久99精品国语久久久| 99热6这里只有精品| 制服丝袜香蕉在线| 国产在视频线精品| 国产在线免费精品| 噜噜噜噜噜久久久久久91| 久久久a久久爽久久v久久| 观看av在线不卡| 国产黄色视频一区二区在线观看| 国产免费一级a男人的天堂| 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 色视频在线一区二区三区| 成年女人在线观看亚洲视频| 欧美精品国产亚洲| 日本欧美视频一区| 色吧在线观看| 国产 一区 欧美 日韩| 久久久久精品久久久久真实原创| 欧美日韩亚洲高清精品| av免费观看日本| 国产中年淑女户外野战色| 身体一侧抽搐| 春色校园在线视频观看| 一级毛片黄色毛片免费观看视频| 午夜福利视频精品| 三级经典国产精品| 国产精品久久久久久av不卡| 三级经典国产精品| 在线观看免费高清a一片| 91狼人影院| 欧美另类一区| 少妇的逼水好多| 国产欧美另类精品又又久久亚洲欧美| 国产久久久一区二区三区| 日本色播在线视频| 妹子高潮喷水视频| 国产一区二区三区av在线| 亚洲欧美日韩另类电影网站 | 天美传媒精品一区二区| 欧美日本视频| 深爱激情五月婷婷| 欧美+日韩+精品| 国产精品一区二区在线不卡| 国产av一区二区精品久久 | 蜜桃久久精品国产亚洲av| 少妇人妻一区二区三区视频| 美女高潮的动态| 国产永久视频网站| 狂野欧美白嫩少妇大欣赏| 人人妻人人澡人人爽人人夜夜| 九色成人免费人妻av| 嘟嘟电影网在线观看| 国产精品一区二区在线观看99| 中国三级夫妇交换| 狂野欧美激情性bbbbbb| 九草在线视频观看| 2018国产大陆天天弄谢| 国产v大片淫在线免费观看| 国内少妇人妻偷人精品xxx网站| 美女内射精品一级片tv| av视频免费观看在线观看| 久久99精品国语久久久| 大香蕉97超碰在线| 欧美一级a爱片免费观看看| 精品国产一区二区三区久久久樱花 | 亚洲内射少妇av| 欧美丝袜亚洲另类| 高清视频免费观看一区二区| 久久99蜜桃精品久久| 啦啦啦视频在线资源免费观看| 99久久人妻综合| 3wmmmm亚洲av在线观看| 免费黄网站久久成人精品| 3wmmmm亚洲av在线观看| 蜜桃久久精品国产亚洲av| 婷婷色综合大香蕉| 成年av动漫网址| 观看av在线不卡| 91精品国产九色| 少妇人妻一区二区三区视频| 亚洲国产成人一精品久久久| 亚洲美女视频黄频| 国产深夜福利视频在线观看| 国产成人免费观看mmmm| 亚洲av成人精品一区久久| 亚洲av男天堂| 久久久a久久爽久久v久久| 国产精品一及| av在线老鸭窝| 九色成人免费人妻av| 久久精品熟女亚洲av麻豆精品| 一级片'在线观看视频| 成人美女网站在线观看视频| 国产男女内射视频| 国产精品国产av在线观看| 国产精品爽爽va在线观看网站| 亚洲人成网站在线观看播放| 国产深夜福利视频在线观看| 亚洲欧洲国产日韩| 乱码一卡2卡4卡精品| 中国三级夫妇交换| 大又大粗又爽又黄少妇毛片口| 各种免费的搞黄视频| 最近最新中文字幕大全电影3| 最近最新中文字幕免费大全7| 乱码一卡2卡4卡精品| 99久久中文字幕三级久久日本| 乱码一卡2卡4卡精品| 亚洲一级一片aⅴ在线观看| 一个人免费看片子| 国产亚洲91精品色在线| 小蜜桃在线观看免费完整版高清| 激情 狠狠 欧美| 日日摸夜夜添夜夜添av毛片| 在线观看国产h片| 国产毛片在线视频| 啦啦啦中文免费视频观看日本| 九草在线视频观看| 国模一区二区三区四区视频| 夜夜看夜夜爽夜夜摸| 日韩免费高清中文字幕av| 熟女电影av网| 色婷婷久久久亚洲欧美| 亚洲精品视频女| 久久99热6这里只有精品| av在线观看视频网站免费| a 毛片基地| 亚洲精品亚洲一区二区| 成人二区视频| 亚洲美女视频黄频| 国产一区有黄有色的免费视频| 在线观看三级黄色| 国产伦精品一区二区三区四那| 伦理电影免费视频| 亚洲三级黄色毛片| 亚洲伊人久久精品综合| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 午夜福利影视在线免费观看| 国产女主播在线喷水免费视频网站| 乱系列少妇在线播放| 国产乱人视频| 男人狂女人下面高潮的视频| 亚洲av成人精品一二三区| 熟妇人妻不卡中文字幕| 久久99蜜桃精品久久| 欧美精品亚洲一区二区| 成年美女黄网站色视频大全免费 | 三级国产精品欧美在线观看| 一区二区三区四区激情视频| 亚洲精品一二三| 在线观看美女被高潮喷水网站| 国产欧美亚洲国产| 18禁裸乳无遮挡免费网站照片| 91狼人影院| 一级毛片电影观看| 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 一本一本综合久久| 大香蕉久久网| 搡老乐熟女国产| 国产日韩欧美在线精品| 国产伦精品一区二区三区视频9| 成人二区视频| 欧美精品人与动牲交sv欧美| 少妇猛男粗大的猛烈进出视频| 小蜜桃在线观看免费完整版高清| 国产成人免费观看mmmm| 日韩一区二区三区影片| 国产av码专区亚洲av| 久久综合国产亚洲精品| 男人添女人高潮全过程视频| 美女主播在线视频| 男人和女人高潮做爰伦理| 人妻制服诱惑在线中文字幕| 多毛熟女@视频| 伊人久久国产一区二区| 日日啪夜夜撸| 国产精品久久久久成人av| 高清视频免费观看一区二区| 2018国产大陆天天弄谢| 久久精品国产鲁丝片午夜精品| 久久热精品热| 国产精品久久久久成人av| 免费大片黄手机在线观看| 多毛熟女@视频| 日日啪夜夜爽| 免费高清在线观看视频在线观看| 性色avwww在线观看| 亚洲,欧美,日韩| 大片电影免费在线观看免费| 99久久人妻综合| 久久99蜜桃精品久久| 日韩制服骚丝袜av| 特大巨黑吊av在线直播| 少妇人妻 视频| 亚洲国产毛片av蜜桃av| 国产亚洲91精品色在线| 日韩一区二区视频免费看| 国产亚洲最大av| 国产精品国产av在线观看| 久久综合国产亚洲精品| 亚洲中文av在线| 人人妻人人添人人爽欧美一区卜 | 制服丝袜香蕉在线| 国模一区二区三区四区视频| 看十八女毛片水多多多| 亚洲,一卡二卡三卡| 国产一区二区三区综合在线观看 | 久久久久精品性色| 亚洲国产日韩一区二区| 亚洲精品久久午夜乱码| 欧美一级a爱片免费观看看| 久久久久久久久大av| 亚洲精品第二区| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 久久久久国产网址| 97在线视频观看| 免费黄网站久久成人精品| 欧美日韩在线观看h| 亚洲人与动物交配视频| 五月开心婷婷网| 欧美精品国产亚洲| 91精品国产九色| 永久网站在线| 在线观看美女被高潮喷水网站| 亚洲国产毛片av蜜桃av| 欧美日韩在线观看h| 国产 精品1| 午夜福利高清视频| 国产 精品1| av天堂中文字幕网| 日本爱情动作片www.在线观看| 赤兔流量卡办理| 自拍欧美九色日韩亚洲蝌蚪91 | 乱系列少妇在线播放| 永久网站在线| 美女内射精品一级片tv| 蜜桃在线观看..| 亚洲一级一片aⅴ在线观看| 人妻系列 视频| 在线天堂最新版资源| 丰满迷人的少妇在线观看| 中文在线观看免费www的网站| 一区二区三区精品91| 国模一区二区三区四区视频| 一个人看视频在线观看www免费| 免费久久久久久久精品成人欧美视频 | 国产精品欧美亚洲77777| 新久久久久国产一级毛片| 免费不卡的大黄色大毛片视频在线观看| 又爽又黄a免费视频| 一区二区av电影网| 内地一区二区视频在线| 日本与韩国留学比较| 少妇人妻精品综合一区二区| 久久久久久九九精品二区国产| 欧美精品亚洲一区二区| 亚洲精品国产色婷婷电影| www.av在线官网国产| 在线免费十八禁| 免费观看a级毛片全部| 亚洲第一av免费看| 日韩成人av中文字幕在线观看| 成年美女黄网站色视频大全免费 | 大香蕉97超碰在线| 亚洲美女黄色视频免费看| 最近最新中文字幕大全电影3| 久久久久人妻精品一区果冻| 熟女电影av网| 简卡轻食公司| 成人国产麻豆网| 午夜免费观看性视频| 国产精品久久久久久av不卡| 中文字幕av成人在线电影| 国产精品久久久久久久久免| 中文精品一卡2卡3卡4更新| 久久久久国产网址| 久久久久网色| 精品久久久久久久久亚洲| 国产高清三级在线| 最近中文字幕高清免费大全6| .国产精品久久| 日本一二三区视频观看| av.在线天堂| 国产精品熟女久久久久浪| 亚洲av日韩在线播放| 乱系列少妇在线播放| h日本视频在线播放| 舔av片在线| 国产伦精品一区二区三区四那| 91久久精品国产一区二区三区| 国产精品人妻久久久久久| 国产成人精品婷婷| 一级爰片在线观看| 亚洲av男天堂| 最近最新中文字幕大全电影3| 国产成人免费观看mmmm| 精品人妻一区二区三区麻豆| www.色视频.com| 蜜桃久久精品国产亚洲av| 亚洲怡红院男人天堂| 精品亚洲成a人片在线观看 | 亚洲精品456在线播放app| 3wmmmm亚洲av在线观看| 久久久久久人妻| 我的女老师完整版在线观看| 国产免费又黄又爽又色| 一级片'在线观看视频| kizo精华| 天堂8中文在线网| 99re6热这里在线精品视频| 好男人视频免费观看在线| 插逼视频在线观看| av不卡在线播放| 精品人妻偷拍中文字幕| 免费黄色在线免费观看| 国产亚洲5aaaaa淫片| 国产精品伦人一区二区| 黄片wwwwww| 乱码一卡2卡4卡精品| 国产精品免费大片| 在线播放无遮挡| 在线观看免费高清a一片| 欧美日韩综合久久久久久| 亚洲四区av| 亚洲av中文av极速乱| 韩国高清视频一区二区三区| 国产精品一及| 十八禁网站网址无遮挡 | 美女中出高潮动态图| 女人十人毛片免费观看3o分钟| 国产女主播在线喷水免费视频网站| 国产久久久一区二区三区| 亚洲国产精品一区三区| 国产精品一区二区性色av| 国产在线一区二区三区精| 一级毛片久久久久久久久女| 国产爱豆传媒在线观看| av黄色大香蕉| 亚洲,一卡二卡三卡|