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

    譜峭度在軸承故障振動信號共振頻帶優(yōu)選中的應用

    2017-06-05 09:35:10于明奇夏均忠陳成法汪治安劉鯤鵬
    軍事交通學院學報 2017年5期
    關鍵詞:峭度頻帶共振

    于明奇,夏均忠,陳成法,汪治安,劉鯤鵬

    (1.軍事交通學院 研究生管理大隊,天津 300161; 2.軍事交通學院 軍用車輛系,天津 300161)

    ?

    譜峭度在軸承故障振動信號共振頻帶優(yōu)選中的應用

    于明奇1,夏均忠2,陳成法2,汪治安1,劉鯤鵬1

    (1.軍事交通學院 研究生管理大隊,天津 300161; 2.軍事交通學院 軍用車輛系,天津 300161)

    在滾動軸承故障振動信號分析中,為減少采集信號共振頻帶對解調效果的影響,利用譜峭度方法確定共振頻帶的帶寬和中心頻率。首先解釋譜峭度的含義,分析基于短時傅里葉變換的譜峭度算法的不足;然后提出快速峭度圖和突起度圖兩種改進方法,并分析其特點;最后通過仿真和試驗手段,對比分析它們在振動信號共振頻帶優(yōu)選中的應用效果。結果表明,突起度圖在信噪比較低的情況下,對共振頻帶的辨別能力優(yōu)于快速峭度圖。

    滾動軸承;故障診斷;譜峭度;快速峭度圖;突起度圖

    滾動軸承廣泛用于旋轉機械中,其運行狀態(tài)是否正常往往直接影響整臺設備的性能,因此對滾動軸承技術狀態(tài)監(jiān)測和早期故障診斷具有重要意義。滾動軸承通常由內圈、外圈、滾動體和保持架組成,當其中某一部分表面發(fā)生局部故障時,產生的脈沖會激起軸承及其相鄰部件的共振,產生調制現(xiàn)象,這使得采集到的振動信號中包含多個共振頻帶[1]。為獲取更好的解調效果,在解調前往往需要選取最優(yōu)共振頻帶,以增強故障特征提取效果。

    譜峭度(Spectral Kurtosis)及其基于短時傅里葉變換(short-time fourier transform,STFT)的算法提出以來[2],常用來確定最優(yōu)共振頻帶的中心頻率和帶寬,但是基于STFT的譜峭度計算過程中,只使用一個固定窗,確定的中心頻率和帶寬不一定最優(yōu)[3]。Antoni引入峭度圖(Kurtogram)概念,提出了快速峭度圖(Fast Kurtogram),使譜峭度得到迅速應用[4],但在信號信噪比較低或者含有隨機脈沖噪聲等情況下容易失效。為此,增強或改進方法不斷提出,其中,Tomasz Barszcz等首先嘗試轉變譜峭度的計算對象,提出了突起度圖(Protrugram)方法,取得了更為理想的效果[5]。

    1 譜峭度

    Antoni從Wold-Cramer表達式角度提出非平穩(wěn)信號譜峭度的定義,并被廣泛認可。Wold-Cramer表達式可以將任何零均值非平穩(wěn)隨機過程x(n)分解為

    (1)

    式中:dZx(f)為正交譜增量;H(n,f)為x(n)在頻率f上的復包絡。

    譜峭度定義為四階歸一化累積量:

    (2)

    基于短時傅里葉變換的譜峭度計算方法為

    (3)

    式中:KY(f)為譜峭度的估計量;fs為采樣頻率;fd為脈沖重復速率;KX為脈沖波動強度;Nw為STFT窗口長度;γ4w為分析窗平方的時間—帶寬積;ρ(f)為信號的信噪比。

    該算法的提出使譜峭度得到實際應用,但利用STFT分解信號時,只使用一個固定窗,因此只分析了部分信號,計算出的譜峭度不是真正意義上的“全域”譜峭度,確定的中心頻率和帶寬不一定最優(yōu)。

    針對該算法存在的不足,引入了峭度圖的概念,并利用準解析濾波器組代替STFT,提出了快速峭度圖。

    2 快速峭度圖

    2.1 快速峭度圖算法

    對于非平穩(wěn)過程,譜峭度的計算依賴估計量頻率分辨率的選擇,譜峭度可以看作是頻率和頻率分辨率(f,Δf)的函數(shù),因此計算譜峭度的過程就轉化為如何在既定頻率上選取頻率分辨率的問題。為便于表達,用峭度圖表示(f,Δf))平面上的譜峭度?;跇錉疃嗨俾蕿V波器組結構的快速峭度圖算法為:

    (1)設置一個低通標準濾波器h(n),截止頻率為fc=1/8+ε(ε≥0),由h(n)構造頻帶低通濾波器h0(n)和高通濾波器h1(n),表達式為

    h0(n)=h(n)ejπn/4
    h1(n)=h(n)ej3πn/4

    (4)

    頻帶位置分別為[0∶1/4]和[1/4∶1/2]。

    (5)

    式中:i=0,1,…,2k-1;k=0,1,…,K-1。

    依據(jù)式(5),可以得出(f,Δf)平面上頻段(fi,Δfk)的峭度圖。

    通過上述算法,峭度圖可以得到快速計算,但在一些窄帶瞬態(tài)采樣時明顯過寬,進而確定的共振頻帶不是十分精確。為此,額外添加了3個準解析帶通濾波器,提出了濾波器組的1/3-二叉樹結構,從而在(f,Δf)平面上的采樣范圍得到明顯改善,并且未影響計算速度。

    2.2 快速峭度圖算法改進

    快速峭度圖提出后得到迅速應用,但在信號存在下列情形時能力受到限制:

    (1)信號的信噪比較低,尤其是含有高峰值的非高斯噪聲;

    (2)信號存在隨機脈沖噪聲;

    (3)故障脈沖重復速率較高。

    通過以下3種方式對快速峭度圖進行改進,增強快速峭度圖的應用效果。

    (1)與其他處理方法結合,增強快速峭度圖應用效果[6-9]。

    (2)改進分解方法。如利用小波包變換、經驗模式分解等替代STFT、準解析濾波器組進行分解,從而達到改進效果[10-13]。

    (3)通過改變譜峭度計算對象來改進算法。如由計算包絡信號的譜峭度得到計算包絡譜、功率譜、平方包絡譜等的譜峭度[14-15]。Tomasz Barszcz和Adam Jabonski首先在這方面做出了嘗試,由計算時域信號的譜峭度轉變?yōu)橛嬎憬浗庹{后的信號包絡譜的譜峭度,提出了突起度圖方法。

    3 突起度圖

    譜峭度定義為四階標準矩,即為:K=μk/σk,其中,μk是關于數(shù)列均值的四階矩,σk是標準差。對于矢量X=[X1,X2,…,XN],譜峭度值計算方法為

    (6)

    突起度圖中譜峭度的意義在于表征信號中脈沖突起強度,在信號中“突起”數(shù)目一定時,其幅值離散程度越大,峭度值也就越大,從而可以確定最大峭度值下共振頻帶的中心頻率,進而確定最優(yōu)共振頻帶。

    突起度圖算法流程如圖1所示。首先確定最優(yōu)帶寬Bw,帶寬選擇滿足兩個約束條件:①盡可能“窄”,要滿足信號中只包含尋找的故障特征信息,避免引入過多噪聲;②盡可能“寬”,以支持多組故障軸承同時診斷,提高運算效率,這在工業(yè)應用中至關重要。然后設置迭代步長step(步長越小,精度越高,但運算時間增加),確定迭代范圍[Bw/2,fs/2](fs為采樣頻率),以Bw/2為初始中心頻率fc,利用帶通濾波器選取頻帶[fc-Bw/2,fc+fs/2],計算該頻帶包絡譜上的譜峭度值,按此循環(huán),直至中心頻率為fs/2。最后取其最大譜峭度值所對應的中心頻率,即為最優(yōu)共振頻帶中心頻率。

    圖1 突起度圖算法流程

    4 仿真分析

    構造周期性振動脈沖仿真信號。單周期脈沖仿真信號x(t)表達式為

    (7)

    式中:t為采樣時刻;fn為載波頻率,fn=1 700Hz;x0為位移常數(shù),x0= 3;θ為阻尼系數(shù),θ=0.1。

    設置重復周期T= 0.01s,即故障特征頻率fm=100Hz,生成周期性振動脈沖信號,如圖2(a)所示。添加高斯白噪聲(信噪比SNR=-10dB),其加噪仿真信號如圖2(b)所示。

    利用快速峭度圖、突起度圖分別對加噪仿真信號進行處理,結果如圖3、圖4所示。圖3中顏色深度代表譜峭度值大小,在中心頻率fc=2 906.25Hz、帶寬Bw=187.5Hz處峭度值最大,其值為0.1,因此,快速峭度圖確定的最優(yōu)共振頻帶中心頻率為2 906.25Hz,帶寬為187.5Hz。

    在突起度計算過程中,為便于之后的解調分析,選擇包含故障特征頻率3至4次諧波的帶寬,設置帶寬Bw= 400Hz>3fm,同時為兼顧計算速度,并不影響識別效果,選擇迭代步長step=100Hz。圖4中,曲線值代表不同中心頻率下的峭度值,在1 600Hz中心頻率處譜峭度值最大,其值為0.016 5,因此,突起度圖所確定的最優(yōu)共振頻帶中心頻率為1 600Hz,帶寬為400Hz。

    (a)仿真信號

    (b)加噪仿真信號圖2 仿真信號時域波形

    圖3 加噪仿真信號快速峭度圖

    圖4 加噪仿真信號突起度圖

    分別應用快速峭度圖、突起度圖選取共振頻帶,濾波后信號的平方包絡譜如圖5所示。圖5(a)中虛線代表故障特征頻率及其倍頻,圖中未能識別與虛線重合的諧波,顯然快速峭度圖選定的共振頻帶不正確;在圖5(b)中能清晰識別設置的故障特征頻率(100 HZ)及其倍頻,說明突起度圖在信號信噪比較低的情況下對共振頻帶辨別能力強于快速峭度圖。

    (a)應用快速峭度圖優(yōu)選共振頻帶

    (b)應用突起度圖優(yōu)選共振頻帶圖5 仿真信號平方包絡譜

    預置帶寬的大小直接影響突起度圖的辨別能力,分別設置帶寬為100 Hz、500 Hz、1 500 Hz、2 500 Hz和4 000 Hz,利用突起度圖分別處理仿真信號、加噪仿真信號,確定其對應中心頻率(見表1)。仿真信號中由于沒有噪聲干擾,其中心頻率收斂于1 650 Hz。加噪仿真信號下,隨著預置帶寬的增大,選定共振頻帶的中心頻率先減小后增大。為直觀表達帶寬選擇對突起度圖的影響,各自共振頻帶的包絡譜如圖6所示,可見帶寬設置過小(100 Hz)或者設置過大(4 000 Hz)均使突起度圖失效(其包絡譜圖中難以識別故障特征頻率),當帶寬設置為500 Hz時,包絡譜中的故障特征頻率及其倍頻最為清晰,這說明帶寬的預置對于突起度圖非常重要??赏ㄟ^支持向量機等方法自適應確定最優(yōu)預置帶寬來改進突起度圖方法。

    表1 預置帶寬及其相應中心頻率 Hz

    圖6 突起度圖優(yōu)選共振頻帶包絡譜

    5 試驗驗證

    選取美國Case Western Reserve University公開的軸承故障數(shù)據(jù)進行故障診斷[6]。試驗裝置由驅動電機、振動加速度傳感器、轉矩解碼/編碼器、聯(lián)軸器和功率計等組成(如圖7所示)。試驗軸承為SKF 6205-2RS深溝球軸承,支撐驅動電機驅動端,其技術參數(shù)見表2。

    圖7 試驗裝置

    滾動體直徑d/mm節(jié)圓直徑D/mm內徑dm/mm外徑Do/mm滾動體數(shù)Z接觸角α/(°)839255290

    對采集的軸承內圈故障振動信號進行共振頻帶優(yōu)選,分別得到快速峭度圖(如圖8所示)、突起度圖(如圖9所示)。

    快速峭度圖利用準解析濾波器組計算譜峭度。從圖8中可以確定最優(yōu)共振頻帶的中心頻率fc=2 625Hz、帶寬Bw=250Hz,其信號的平方包絡譜如圖10(a)所示。從圖中較難識別軸承內圈故障特征頻率(158Hz)及其倍頻。

    突起度圖方法中,設置帶寬Bw=500Hz,選擇迭代步長step= 100Hz。從圖9中可以確定在中心頻率fc=3 550Hz時譜峭度值最大,即最優(yōu)共振頻帶中心頻率fc=3 550Hz、帶寬Bw=500Hz,其信號的平方包絡譜如圖10(b)所示。從圖中能夠清晰識別軸承內圈故障特征頻率(158 Hz)及其倍頻(與理論值157.94 Hz不完全一致的原因是由于滾動體隨機滑動的影響)。

    圖8 軸承內圈故障振動信號快速峭度圖

    圖9 軸承內圈故障振動信號突起度圖

    (a)應用快速峭度圖優(yōu)選共振頻帶

    (b)應用突起度圖優(yōu)選共振頻帶圖10 軸承內圈故障振動信號平方包絡譜

    通過軸承內圈故障振動信號分析可知,由突起度圖確定的共振頻帶,其平方包絡譜中故障特征更為突出,這說明在實際應用中突起度圖相比于快速峭度圖,在振動信號共振頻帶優(yōu)選中較優(yōu)。

    6 結 語

    本文首先系統(tǒng)研究了譜峭度的含義,指出基于STFT的譜峭度算法只使用一個固定窗,確定的共振頻帶不一定最優(yōu);針對該算法存在的不足,提出了快速峭度圖方法,但其在信號信噪比較低或者含有隨機脈沖噪聲等情況下容易失效;然后提出改進思路,重點分析了突起度圖方法;最后通過仿真和試驗手段分析驗證了快速峭度圖和突起度圖在低信噪比以及實際應用中的效果。突起度圖方法在信噪比較低的情況下,對共振頻帶的優(yōu)選能力較之快速峭度圖更為突出,并且在實際應用中取得了更為理想的效果。

    應用譜峭度方法可以迅速準確確定最優(yōu)共振頻帶,一方面可以增強滾動軸承故障診斷效果,另一方面可以極大提高運算速度,為實現(xiàn)在線技術狀態(tài)監(jiān)測提供了可能。如何改進或增強其應用效果依舊是未來研究重點:一是融合多種改進方式;二是改進完善譜峭度指標,比如,引入熱力學中熵的概念代替?zhèn)鹘y(tǒng)意義上的譜峭度等。

    [1] RANDALL R B, ANTONI J. Rolling element bearing diagnostics:a tutorial[J]. Mechanical Systems and Signal Processing, 2011, 25(2):485-520.

    [2] ANTONI J. The spectral kurtosis: a useful tool for characterizing nonstationarysignals[J]. Mechanical Systems and Signal Processing, 2006, 20(2): 282-307.

    [3] ANTONI J,RANDALL R B. The spectral kurtosis: application to the vibratory surveillance and diagnostics of rotating machines[J]. Mechanical Systems and Signal Processing, 2006, 20(2): 308-331.

    [4] ANTONI J. Fast computation of the Kurtogram for the detection of transient faults[J]. Mechanical Systems and Signal Processing, 2007, 21(1):108-124.

    [5] BARSZCZ T,JABLONSKi A. A novel method for the optimal band selection for vibration signal demodulation and comparison with the kurtogram[J]. Mechanical Systems and Signal Processing, 2011, 25(1): 431-451.

    [6] GUO W, TSE P W, A Djordjevich. Faulty bearing signal recovery from large noise using a hybrid method based on spectral kurtosis and ensemble empirical mode decomposition[J]. Measurement, 2012, 45(5): 1308-1322.

    [7] 彭暢,柏林,謝小亮.基于EEMD、度量因子和快速峭度圖的滾動軸承故障診斷[J].振動與沖擊,2012,31(20):143-146.

    [8] HUANG W,LIU H,ZHU Z,et al. Fault feature extracting for rotating machinery vibration based on blind deconvolution and spectral kurtosis[J]. Mechanical Systems and Signal Processing, 2012, 43(2): 261-266.

    [9] 周智,朱永生,張優(yōu)云,等.基于MMSE和譜峭度的滾動軸承故障診斷方法[J].振動與沖擊,2013,32(6):73-77.

    [10] LEI Y G, LIN J, HE Z J, et al. Application of an improved kurtogram method for fault diagnosis of rolling element bearings[J]. Mechanical Systems and Signal Processing, 2011, 25(5): 1738-1749.

    [11] 程軍圣,楊怡,楊宇.基于LMD的譜峭度方法在齒輪故障診斷中的應用[J].振動與沖擊,2012,31(18):20-23.

    [12] 代士超,郭瑜,伍星,等.基于子頻帶譜峭度平均的快速譜峭度圖算法改進[J].振動與沖擊,2015,34(7):98-102.

    [13] 彭暢,柏林,劉小峰.基于魯棒性小波包峭度圖的滾動軸承故障診斷[J].振動、測試與診斷,2016,36(1):11-16.

    [14] WANG D,TSE P W,TSUI K L. An enhanced Kurtogram method for fault diagnosis of rolling element bearings[J].Mechanical Systems and Signal Processing, 2013, 35(1-2): 176-199.

    [15] ANTONI J. The infogram: entropic evidence of the signature of repetitive transients[J]. Mechanical Systems and Signal Processing, 2016, 74(1): 73-94.

    [16] SMITH W A,RANDALL R B.Rolling element bearing diagnostics using the Case Western Reserve University data:a benchmark study[J]. Mechanical Systems and Signal Processing, 2015, 64(2): 100-131.

    (編輯:張峰)

    Application of Spectral Kurtosis in Optimal Selection of Resonance Frequency Bands for Bearing Fault Vibration Signal

    YU Mingqi1, XIA Junzhong2, CHEN Chengfa2, WANG Zhian1, LIU Kunpeng1

    (1.Postgraduate Training Brigade, Military Transportation University, Tianjin 300161, China; 2.Military Vehicle Department, Military Transportation University, Tianjin 300161, China)

    To reduce the impact of signal resonance frequency bands on demodulation effect while analyzing rolling bearing fault vibration signal, the paper determines bandwidth and center frequency of resonance frequency bands with spectral kurtosis. It firstly interprets the meaning of spectral kurtosis and analyzes the deficiency of spectral kurtosis algorithm based on short-time Fourier transform. Then, it puts forward two improved methods: fast kurtogram and protrugram, and analyzes their characteristics. Finally, it compares their application effect in optimal selection of vibration signal resonance frequency bands. The result shows that protrugram has better distinguishing ability than fast kurtogram on resonance frequency bands in low SNR.

    rolling bearing; fault diagnosis; spectral kurtosis; fast kurtogram; protrugram

    2016-12-30;

    2017-01-13.

    于明奇(1991—),男,碩士研究生; 夏均忠(1967—),男,博士,教授,碩士研究生導師.

    10.16807/j.cnki.12-1372/e.2017.05.010

    TN911.23;TB53

    A

    1674-2192(2017)05- 0042- 06

    ● 車輛工程 Vehicle Engineering

    猜你喜歡
    峭度頻帶共振
    基于MCKD和峭度的液壓泵故障特征提取
    機床與液壓(2023年1期)2023-02-03 10:14:18
    Wi-Fi網絡中5G和2.4G是什么?有何區(qū)別?
    聯(lián)合快速峭度圖與變帶寬包絡譜峭度圖的輪對軸承復合故障檢測研究
    單音及部分頻帶干擾下DSSS系統(tǒng)性能分析
    安然 與時代同頻共振
    選硬人打硬仗——紫陽縣黨建與脫貧同頻共振
    當代陜西(2018年12期)2018-08-04 05:49:22
    基于峭度分析的聲發(fā)射故障檢測
    電子世界(2018年12期)2018-07-04 06:34:38
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    CTA 中紡院+ 化纖聯(lián)盟 強強聯(lián)合 科技共振
    基于魯棒性小波包峭度圖的滾動軸承故障診斷*
    国产精品久久电影中文字幕 | 久久久久精品国产欧美久久久| 久久久久久久精品吃奶| a级毛片黄视频| 可以免费在线观看a视频的电影网站| 宅男免费午夜| 高清av免费在线| 无人区码免费观看不卡 | 欧美 日韩 精品 国产| 999久久久国产精品视频| 妹子高潮喷水视频| 欧美乱码精品一区二区三区| 两个人看的免费小视频| 高清av免费在线| 蜜桃在线观看..| 高清欧美精品videossex| 日本撒尿小便嘘嘘汇集6| 色老头精品视频在线观看| 色视频在线一区二区三区| 国产一区二区 视频在线| 中文字幕另类日韩欧美亚洲嫩草| 欧美日本中文国产一区发布| 国产精品自产拍在线观看55亚洲 | 99re在线观看精品视频| 丁香欧美五月| h视频一区二区三区| 80岁老熟妇乱子伦牲交| 一本一本久久a久久精品综合妖精| 五月开心婷婷网| 精品国产一区二区三区久久久樱花| 国产精品偷伦视频观看了| 色尼玛亚洲综合影院| 国产精品国产高清国产av | 大陆偷拍与自拍| 免费在线观看黄色视频的| 国产一区二区 视频在线| 一边摸一边做爽爽视频免费| 国产在线视频一区二区| 无人区码免费观看不卡 | 久久精品国产亚洲av高清一级| 亚洲人成伊人成综合网2020| 男男h啪啪无遮挡| 好男人电影高清在线观看| 精品少妇黑人巨大在线播放| 午夜老司机福利片| 欧美日韩国产mv在线观看视频| 男人舔女人的私密视频| 91av网站免费观看| 大片电影免费在线观看免费| 一级a爱视频在线免费观看| 久久精品熟女亚洲av麻豆精品| 精品免费久久久久久久清纯 | 女人爽到高潮嗷嗷叫在线视频| 久久久久久人人人人人| 成年人免费黄色播放视频| 国产熟女午夜一区二区三区| 亚洲av成人一区二区三| 亚洲国产看品久久| 99国产精品一区二区蜜桃av | 热99国产精品久久久久久7| 日本五十路高清| 视频区欧美日本亚洲| 免费在线观看黄色视频的| 在线观看免费高清a一片| 国产成人欧美| 国产精品九九99| 操出白浆在线播放| 午夜福利影视在线免费观看| 夜夜骑夜夜射夜夜干| 在线观看免费日韩欧美大片| 欧美人与性动交α欧美软件| 法律面前人人平等表现在哪些方面| 操美女的视频在线观看| www.熟女人妻精品国产| 亚洲成人手机| 亚洲成人手机| 国产1区2区3区精品| 2018国产大陆天天弄谢| 精品国内亚洲2022精品成人 | 国产不卡一卡二| 伊人久久大香线蕉亚洲五| 国产真人三级小视频在线观看| 国产无遮挡羞羞视频在线观看| 欧美精品啪啪一区二区三区| 热re99久久国产66热| 一本综合久久免费| 欧美国产精品va在线观看不卡| av免费在线观看网站| 久久婷婷成人综合色麻豆| 国产精品久久久久久精品古装| 亚洲av成人不卡在线观看播放网| 日韩精品免费视频一区二区三区| 免费少妇av软件| 大型黄色视频在线免费观看| 麻豆乱淫一区二区| 水蜜桃什么品种好| 动漫黄色视频在线观看| 建设人人有责人人尽责人人享有的| 国产精品亚洲一级av第二区| 一边摸一边抽搐一进一出视频| 丰满少妇做爰视频| 大片电影免费在线观看免费| 国产在线精品亚洲第一网站| 国产高清videossex| 久久久久国产一级毛片高清牌| 麻豆av在线久日| 满18在线观看网站| 久久中文字幕一级| 露出奶头的视频| 国内毛片毛片毛片毛片毛片| 丝袜美腿诱惑在线| 99国产综合亚洲精品| 黄片大片在线免费观看| 国产成人欧美在线观看 | 国产精品一区二区在线不卡| 男女下面插进去视频免费观看| 日韩欧美一区二区三区在线观看 | 丰满迷人的少妇在线观看| av电影中文网址| 一级黄色大片毛片| 成年人免费黄色播放视频| 中文字幕人妻丝袜制服| 9热在线视频观看99| 午夜激情久久久久久久| 亚洲国产精品一区二区三区在线| 啦啦啦免费观看视频1| 夫妻午夜视频| 后天国语完整版免费观看| 免费观看人在逋| 天堂动漫精品| 欧美+亚洲+日韩+国产| 亚洲午夜精品一区,二区,三区| 啦啦啦中文免费视频观看日本| 精品欧美一区二区三区在线| 欧美乱妇无乱码| 国产成人一区二区三区免费视频网站| 超碰成人久久| 欧美变态另类bdsm刘玥| 免费观看a级毛片全部| 亚洲 欧美一区二区三区| 中文字幕制服av| 国产免费现黄频在线看| 桃花免费在线播放| 一二三四在线观看免费中文在| 日本av手机在线免费观看| 香蕉国产在线看| 91精品三级在线观看| 午夜福利视频在线观看免费| 三上悠亚av全集在线观看| 9热在线视频观看99| 香蕉久久夜色| 国产1区2区3区精品| 久久热在线av| 国产亚洲午夜精品一区二区久久| 国产97色在线日韩免费| 免费在线观看黄色视频的| 十八禁人妻一区二区| cao死你这个sao货| 国产伦人伦偷精品视频| 俄罗斯特黄特色一大片| 国产成人精品在线电影| 在线天堂中文资源库| 精品国产一区二区三区四区第35| 91成人精品电影| 伊人久久大香线蕉亚洲五| 妹子高潮喷水视频| 精品国产亚洲在线| 亚洲三区欧美一区| 久久国产精品影院| 9热在线视频观看99| 欧美精品啪啪一区二区三区| 99久久99久久久精品蜜桃| 99精品在免费线老司机午夜| 黄片播放在线免费| 久久影院123| 亚洲国产av影院在线观看| 免费观看a级毛片全部| 一区二区三区精品91| 一级,二级,三级黄色视频| 人成视频在线观看免费观看| 久久人妻熟女aⅴ| 叶爱在线成人免费视频播放| 人人妻人人澡人人爽人人夜夜| 亚洲熟妇熟女久久| 捣出白浆h1v1| 精品人妻熟女毛片av久久网站| 亚洲熟女精品中文字幕| 欧美日韩成人在线一区二区| 一边摸一边抽搐一进一小说 | 正在播放国产对白刺激| 国产成人啪精品午夜网站| 一本色道久久久久久精品综合| 日本撒尿小便嘘嘘汇集6| 三级毛片av免费| 亚洲一卡2卡3卡4卡5卡精品中文| 制服诱惑二区| 欧美人与性动交α欧美精品济南到| 青草久久国产| 久久精品亚洲av国产电影网| 在线观看免费日韩欧美大片| 一边摸一边抽搐一进一小说 | 精品免费久久久久久久清纯 | 性高湖久久久久久久久免费观看| 大香蕉久久网| 性高湖久久久久久久久免费观看| 午夜福利免费观看在线| 露出奶头的视频| 亚洲av日韩精品久久久久久密| 国产在线视频一区二区| 国产欧美亚洲国产| 亚洲av日韩精品久久久久久密| 人人妻人人澡人人爽人人夜夜| 老司机影院毛片| 99精国产麻豆久久婷婷| 久久人妻福利社区极品人妻图片| 久久人妻福利社区极品人妻图片| 亚洲av第一区精品v没综合| 51午夜福利影视在线观看| 一夜夜www| 大型av网站在线播放| 一级片免费观看大全| 男女边摸边吃奶| 国产成人影院久久av| 91国产中文字幕| 国产欧美日韩一区二区精品| 在线观看免费午夜福利视频| 人人妻人人澡人人看| 桃花免费在线播放| 一级片'在线观看视频| 国产成人啪精品午夜网站| 男女高潮啪啪啪动态图| 久久久久精品国产欧美久久久| 久久ye,这里只有精品| 高潮久久久久久久久久久不卡| 午夜免费鲁丝| 在线观看免费高清a一片| netflix在线观看网站| 久久精品亚洲精品国产色婷小说| 色视频在线一区二区三区| 伦理电影免费视频| 人人妻人人添人人爽欧美一区卜| 国产av国产精品国产| 亚洲午夜理论影院| 亚洲va日本ⅴa欧美va伊人久久| 97人妻天天添夜夜摸| 欧美成人午夜精品| 动漫黄色视频在线观看| 三上悠亚av全集在线观看| 丁香六月天网| 一区福利在线观看| 黄片小视频在线播放| 亚洲熟妇熟女久久| 中文字幕人妻丝袜制服| 一边摸一边抽搐一进一小说 | 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到| 国产成人一区二区三区免费视频网站| 首页视频小说图片口味搜索| 免费av中文字幕在线| 脱女人内裤的视频| 香蕉久久夜色| 亚洲精品成人av观看孕妇| 老司机靠b影院| 无遮挡黄片免费观看| 久热这里只有精品99| 中文字幕人妻丝袜一区二区| 精品久久久久久久毛片微露脸| 高清黄色对白视频在线免费看| 亚洲一区二区三区欧美精品| 成人永久免费在线观看视频 | 1024香蕉在线观看| 91成年电影在线观看| 在线看a的网站| 亚洲人成电影观看| 热99re8久久精品国产| 亚洲国产欧美在线一区| 久久人妻av系列| 亚洲国产欧美网| 手机成人av网站| 国产精品久久久久久精品古装| 亚洲国产毛片av蜜桃av| 天天操日日干夜夜撸| 高清在线国产一区| 另类精品久久| 露出奶头的视频| 成人三级做爰电影| 亚洲欧美日韩另类电影网站| 男女床上黄色一级片免费看| 99精品久久久久人妻精品| 精品久久久精品久久久| 久久精品熟女亚洲av麻豆精品| 亚洲一区二区三区欧美精品| xxxhd国产人妻xxx| 新久久久久国产一级毛片| 交换朋友夫妻互换小说| 欧美久久黑人一区二区| 日韩免费av在线播放| 黄色视频,在线免费观看| 91麻豆av在线| 亚洲成人免费电影在线观看| 99在线人妻在线中文字幕 | 美女扒开内裤让男人捅视频| 国产视频一区二区在线看| 国产精品二区激情视频| 国产成人一区二区三区免费视频网站| 欧美 亚洲 国产 日韩一| 久热这里只有精品99| 欧美人与性动交α欧美精品济南到| 亚洲av片天天在线观看| 亚洲国产毛片av蜜桃av| 日韩大码丰满熟妇| 窝窝影院91人妻| 国产成人系列免费观看| 男女边摸边吃奶| 啦啦啦在线免费观看视频4| 免费看a级黄色片| 精品福利永久在线观看| 国产aⅴ精品一区二区三区波| 成人av一区二区三区在线看| 亚洲九九香蕉| 男女高潮啪啪啪动态图| 18禁观看日本| 亚洲精品中文字幕一二三四区 | 国产欧美亚洲国产| 国精品久久久久久国模美| 国产精品免费一区二区三区在线 | 精品少妇久久久久久888优播| 国产91精品成人一区二区三区 | 欧美激情极品国产一区二区三区| 制服人妻中文乱码| 亚洲精品一二三| 国产精品久久电影中文字幕 | 国产精品免费一区二区三区在线 | 在线观看免费高清a一片| 99在线人妻在线中文字幕 | 国产成人欧美| 成人黄色视频免费在线看| 天堂俺去俺来也www色官网| av超薄肉色丝袜交足视频| 一本大道久久a久久精品| 91成年电影在线观看| 国产免费视频播放在线视频| 亚洲中文av在线| 精品国产乱子伦一区二区三区| 欧美午夜高清在线| 亚洲国产中文字幕在线视频| av欧美777| 高清av免费在线| 亚洲欧美色中文字幕在线| 桃花免费在线播放| 丁香六月天网| 老司机午夜福利在线观看视频 | bbb黄色大片| 悠悠久久av| 中文欧美无线码| 国产日韩欧美视频二区| 高清毛片免费观看视频网站 | 男女高潮啪啪啪动态图| 国产91精品成人一区二区三区 | 久久精品人人爽人人爽视色| 黑人操中国人逼视频| 久久久水蜜桃国产精品网| 老司机福利观看| 成人特级黄色片久久久久久久 | 国产黄色免费在线视频| 咕卡用的链子| 99精品欧美一区二区三区四区| 欧美黄色淫秽网站| 一边摸一边做爽爽视频免费| 久久国产精品大桥未久av| 成人18禁高潮啪啪吃奶动态图| 丝袜喷水一区| 制服诱惑二区| 下体分泌物呈黄色| 五月开心婷婷网| √禁漫天堂资源中文www| 男人舔女人的私密视频| 久久亚洲真实| 久久国产亚洲av麻豆专区| 另类亚洲欧美激情| 亚洲免费av在线视频| 日韩中文字幕欧美一区二区| 欧美激情高清一区二区三区| 久久人人爽av亚洲精品天堂| 久久狼人影院| 黄片大片在线免费观看| 久久精品亚洲熟妇少妇任你| 国产在线免费精品| 国产人伦9x9x在线观看| 国产成人系列免费观看| 免费在线观看黄色视频的| 精品少妇一区二区三区视频日本电影| 丝袜美腿诱惑在线| 色视频在线一区二区三区| 久久精品熟女亚洲av麻豆精品| 精品国产一区二区久久| 国产视频一区二区在线看| 性色av乱码一区二区三区2| 一级毛片女人18水好多| 日韩精品免费视频一区二区三区| 国产精品美女特级片免费视频播放器 | 高清黄色对白视频在线免费看| 久久青草综合色| 精品一品国产午夜福利视频| 久久影院123| 日韩欧美国产一区二区入口| 久久九九热精品免费| 9热在线视频观看99| 黑丝袜美女国产一区| 丰满迷人的少妇在线观看| 国产又爽黄色视频| 黄色片一级片一级黄色片| 又紧又爽又黄一区二区| 亚洲性夜色夜夜综合| 50天的宝宝边吃奶边哭怎么回事| 国产淫语在线视频| 国产成人一区二区三区免费视频网站| 日本黄色视频三级网站网址 | 国产精品成人在线| 9热在线视频观看99| 丝袜美腿诱惑在线| 国产免费视频播放在线视频| 日韩欧美三级三区| √禁漫天堂资源中文www| 国产精品麻豆人妻色哟哟久久| 99精品久久久久人妻精品| 中文亚洲av片在线观看爽 | 女同久久另类99精品国产91| 69av精品久久久久久 | 99精国产麻豆久久婷婷| 看免费av毛片| 国产精品免费一区二区三区在线 | 国产成人一区二区三区免费视频网站| 国产亚洲一区二区精品| 电影成人av| 欧美人与性动交α欧美软件| 99精品在免费线老司机午夜| 久久久国产欧美日韩av| 免费在线观看黄色视频的| 久久久精品区二区三区| 天天操日日干夜夜撸| 又黄又粗又硬又大视频| 成人亚洲精品一区在线观看| 久久久久久久久免费视频了| 国产av一区二区精品久久| 中文字幕色久视频| 久久精品国产a三级三级三级| 日本av免费视频播放| 九色亚洲精品在线播放| 免费观看av网站的网址| 久久精品国产综合久久久| 成年人免费黄色播放视频| 国产在线精品亚洲第一网站| 欧美黑人精品巨大| av又黄又爽大尺度在线免费看| 午夜免费成人在线视频| 别揉我奶头~嗯~啊~动态视频| 国产视频一区二区在线看| 母亲3免费完整高清在线观看| 欧美国产精品va在线观看不卡| 叶爱在线成人免费视频播放| 亚洲成人手机| 国产欧美亚洲国产| 日本av手机在线免费观看| 亚洲第一青青草原| 欧美精品av麻豆av| av在线播放免费不卡| 18禁裸乳无遮挡动漫免费视频| 自线自在国产av| 日韩有码中文字幕| 欧美亚洲日本最大视频资源| 最新的欧美精品一区二区| 一级a爱视频在线免费观看| 亚洲国产精品一区二区三区在线| 日韩免费高清中文字幕av| 国产免费现黄频在线看| 免费一级毛片在线播放高清视频 | 69精品国产乱码久久久| 成人特级黄色片久久久久久久 | 欧美日韩亚洲国产一区二区在线观看 | 久久午夜综合久久蜜桃| 80岁老熟妇乱子伦牲交| 老司机福利观看| 日本欧美视频一区| 精品欧美一区二区三区在线| 波多野结衣av一区二区av| 日本wwww免费看| 亚洲av成人不卡在线观看播放网| 久久香蕉激情| 亚洲精品中文字幕在线视频| 亚洲人成电影观看| 午夜两性在线视频| 亚洲精品自拍成人| 久久国产精品影院| 亚洲一区二区三区欧美精品| 一级a爱视频在线免费观看| 日韩人妻精品一区2区三区| 757午夜福利合集在线观看| 精品免费久久久久久久清纯 | 久久久精品区二区三区| 国产黄频视频在线观看| 国产视频一区二区在线看| 色精品久久人妻99蜜桃| 变态另类成人亚洲欧美熟女 | 久久影院123| 国产伦理片在线播放av一区| 国产高清videossex| 考比视频在线观看| 多毛熟女@视频| 久久精品aⅴ一区二区三区四区| 亚洲欧洲日产国产| 一区福利在线观看| √禁漫天堂资源中文www| 国产av一区二区精品久久| 国产片内射在线| 欧美老熟妇乱子伦牲交| 国产高清国产精品国产三级| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 91老司机精品| 亚洲熟妇熟女久久| 国产成人免费无遮挡视频| 午夜福利免费观看在线| 黑人操中国人逼视频| 国产精品av久久久久免费| 一本大道久久a久久精品| 久久人妻熟女aⅴ| 亚洲一区中文字幕在线| 人成视频在线观看免费观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品成人免费网站| 国产亚洲av高清不卡| 日韩大码丰满熟妇| 妹子高潮喷水视频| 最黄视频免费看| www.自偷自拍.com| 精品久久久久久电影网| 日本五十路高清| 中文字幕色久视频| 亚洲精品久久成人aⅴ小说| 亚洲精品国产色婷婷电影| 国产91精品成人一区二区三区 | 捣出白浆h1v1| 亚洲中文日韩欧美视频| 亚洲av电影在线进入| 91九色精品人成在线观看| 国产伦人伦偷精品视频| 国产精品国产高清国产av | 亚洲人成电影免费在线| 天天躁夜夜躁狠狠躁躁| 99热网站在线观看| 欧美av亚洲av综合av国产av| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区日韩欧美中文字幕| av一本久久久久| 日韩欧美三级三区| 亚洲精品av麻豆狂野| 视频区欧美日本亚洲| 日韩免费高清中文字幕av| 久久久久久久精品吃奶| 国产亚洲av高清不卡| 久9热在线精品视频| 亚洲全国av大片| 国产精品熟女久久久久浪| cao死你这个sao货| 国产亚洲av高清不卡| 久久久国产成人免费| 精品亚洲乱码少妇综合久久| 少妇的丰满在线观看| 久久精品国产亚洲av香蕉五月 | 757午夜福利合集在线观看| av福利片在线| 国产精品久久久人人做人人爽| 咕卡用的链子| 亚洲色图 男人天堂 中文字幕| 午夜福利欧美成人| 免费一级毛片在线播放高清视频 | 91字幕亚洲| 黑人猛操日本美女一级片| 国产精品av久久久久免费| 国产精品久久久久久精品古装| 国产男女内射视频| 亚洲av日韩精品久久久久久密| 午夜激情av网站| 欧美日韩一级在线毛片| 精品人妻在线不人妻| 精品久久久精品久久久| 免费观看a级毛片全部| 国产高清激情床上av| 91九色精品人成在线观看| tocl精华| 国产淫语在线视频| 亚洲欧美一区二区三区久久| 久久午夜亚洲精品久久| 午夜福利乱码中文字幕| 一区二区三区乱码不卡18| 欧美乱妇无乱码| 精品一品国产午夜福利视频| 久久精品国产综合久久久| 欧美日韩视频精品一区| 国产高清激情床上av| 国产亚洲精品一区二区www | 日本av手机在线免费观看| 亚洲国产欧美一区二区综合| 午夜免费成人在线视频| 精品国产乱子伦一区二区三区| 亚洲国产欧美一区二区综合| 国产野战对白在线观看| 国产av国产精品国产| 汤姆久久久久久久影院中文字幕| 亚洲国产毛片av蜜桃av| 中文字幕高清在线视频| 桃红色精品国产亚洲av| 日本欧美视频一区|