• <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)合 科技共振
    基于魯棒性小波包峭度圖的滾動軸承故障診斷*
    成人特级黄色片久久久久久久| 看免费av毛片| 中文字幕人妻熟女乱码| 亚洲精华国产精华精| 日韩大码丰满熟妇| 欧美乱妇无乱码| 后天国语完整版免费观看| 欧美中文日本在线观看视频| 精品久久久久久久久久久久久 | 人人澡人人妻人| 久久午夜综合久久蜜桃| 亚洲在线自拍视频| 亚洲 欧美一区二区三区| 男女那种视频在线观看| 日本一本二区三区精品| 中文字幕久久专区| 久久久久久久久中文| 亚洲欧美精品综合一区二区三区| 一级毛片女人18水好多| 国产成人av激情在线播放| 国产片内射在线| 亚洲精品国产一区二区精华液| 亚洲一区高清亚洲精品| e午夜精品久久久久久久| 欧美成人一区二区免费高清观看 | 欧美精品亚洲一区二区| 黄色毛片三级朝国网站| 久久青草综合色| 99久久99久久久精品蜜桃| 国产视频一区二区在线看| 黄色 视频免费看| 日韩精品青青久久久久久| 久久人人精品亚洲av| 精品国产美女av久久久久小说| 淫秽高清视频在线观看| 极品教师在线免费播放| 欧美黑人巨大hd| 日韩有码中文字幕| 精品日产1卡2卡| 欧美另类亚洲清纯唯美| 级片在线观看| 日韩大尺度精品在线看网址| 亚洲专区字幕在线| 国产精品乱码一区二三区的特点| 一级毛片女人18水好多| 黄片播放在线免费| 制服人妻中文乱码| 观看免费一级毛片| 国产高清videossex| 九色国产91popny在线| xxx96com| 国产午夜精品久久久久久| 亚洲天堂国产精品一区在线| 特大巨黑吊av在线直播 | 99热6这里只有精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲在线自拍视频| 中文字幕av电影在线播放| 国产精品永久免费网站| 国产熟女xx| 亚洲欧美激情综合另类| 99热6这里只有精品| 精品欧美一区二区三区在线| 亚洲一区高清亚洲精品| 午夜成年电影在线免费观看| 欧美一级a爱片免费观看看 | 91大片在线观看| 窝窝影院91人妻| 一级毛片精品| 国产精品电影一区二区三区| 久久青草综合色| 50天的宝宝边吃奶边哭怎么回事| 国产精品野战在线观看| 一本综合久久免费| 色老头精品视频在线观看| 身体一侧抽搐| 女人高潮潮喷娇喘18禁视频| 日本黄色视频三级网站网址| 老司机靠b影院| 国产又爽黄色视频| av中文乱码字幕在线| 国产成人欧美在线观看| 桃色一区二区三区在线观看| 高清在线国产一区| 国产精品野战在线观看| 岛国在线观看网站| 亚洲国产高清在线一区二区三 | 日韩欧美一区二区三区在线观看| 国产av不卡久久| 51午夜福利影视在线观看| 青草久久国产| 欧美丝袜亚洲另类 | 大香蕉久久成人网| 丁香六月欧美| 女人高潮潮喷娇喘18禁视频| 久久亚洲真实| 我的亚洲天堂| 人人妻人人澡欧美一区二区| 男女之事视频高清在线观看| 亚洲天堂国产精品一区在线| 欧美日韩黄片免| 哪里可以看免费的av片| 少妇的丰满在线观看| 在线十欧美十亚洲十日本专区| 男女下面进入的视频免费午夜 | 欧美日韩瑟瑟在线播放| av欧美777| 国产精品 国内视频| 成人三级做爰电影| 国产精品永久免费网站| 国产亚洲欧美在线一区二区| 日韩视频一区二区在线观看| 99久久无色码亚洲精品果冻| 两性夫妻黄色片| 1024手机看黄色片| 免费看日本二区| 50天的宝宝边吃奶边哭怎么回事| 757午夜福利合集在线观看| 一区二区三区高清视频在线| 国产主播在线观看一区二区| 一个人免费在线观看的高清视频| 桃红色精品国产亚洲av| 国产视频内射| 人人澡人人妻人| 亚洲午夜精品一区,二区,三区| 自线自在国产av| 国产欧美日韩精品亚洲av| 色精品久久人妻99蜜桃| 在线av久久热| 精品卡一卡二卡四卡免费| 精品一区二区三区视频在线观看免费| 日韩成人在线观看一区二区三区| 久久亚洲精品不卡| 深夜精品福利| 一个人观看的视频www高清免费观看 | 亚洲中文av在线| 久久精品91无色码中文字幕| 看黄色毛片网站| 波多野结衣av一区二区av| 99久久无色码亚洲精品果冻| 久久 成人 亚洲| 国产野战对白在线观看| 免费看十八禁软件| 久久精品国产综合久久久| 啦啦啦免费观看视频1| 久久精品成人免费网站| 欧美乱码精品一区二区三区| 亚洲第一青青草原| xxxwww97欧美| 久久精品91蜜桃| 久久久久久大精品| 亚洲欧美精品综合久久99| 欧美精品亚洲一区二区| 天堂√8在线中文| 日韩一卡2卡3卡4卡2021年| 美女国产高潮福利片在线看| 国产一区二区激情短视频| 久久久久久亚洲精品国产蜜桃av| 夜夜看夜夜爽夜夜摸| av有码第一页| 亚洲欧美精品综合久久99| 中文亚洲av片在线观看爽| 国产午夜精品久久久久久| www国产在线视频色| 变态另类成人亚洲欧美熟女| 亚洲av电影在线进入| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产高清在线一区二区三 | 91老司机精品| 精品久久久久久,| 日本 欧美在线| 精品一区二区三区av网在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲美女久久久| 日本 av在线| 国产高清有码在线观看视频 | 国产亚洲精品综合一区在线观看 | 满18在线观看网站| 成人18禁在线播放| 午夜影院日韩av| 亚洲五月天丁香| 亚洲精品av麻豆狂野| 亚洲国产精品999在线| www.www免费av| 看片在线看免费视频| 搡老岳熟女国产| 久久性视频一级片| 欧美精品啪啪一区二区三区| 美女高潮到喷水免费观看| 老司机福利观看| 亚洲 欧美 日韩 在线 免费| 搡老妇女老女人老熟妇| 免费在线观看日本一区| 久久精品夜夜夜夜夜久久蜜豆 | 99久久精品国产亚洲精品| 激情在线观看视频在线高清| 天天添夜夜摸| xxx96com| 亚洲欧美日韩无卡精品| 99精品久久久久人妻精品| 日韩国内少妇激情av| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 亚洲欧美精品综合一区二区三区| 国产精品九九99| 国产区一区二久久| 国产成人欧美| 国产高清有码在线观看视频 | 亚洲一区二区三区不卡视频| 亚洲国产欧洲综合997久久, | 1024香蕉在线观看| 99热6这里只有精品| 在线观看66精品国产| 久久久久久久午夜电影| 欧美久久黑人一区二区| 男女下面进入的视频免费午夜 | 日本黄色视频三级网站网址| 69av精品久久久久久| 中文字幕av电影在线播放| 午夜a级毛片| 中文字幕人妻熟女乱码| 国产乱人伦免费视频| 色综合站精品国产| 精品一区二区三区av网在线观看| 草草在线视频免费看| 美女 人体艺术 gogo| 一级a爱视频在线免费观看| 国产精品日韩av在线免费观看| 校园春色视频在线观看| 国产亚洲精品第一综合不卡| 久9热在线精品视频| 亚洲国产精品合色在线| 叶爱在线成人免费视频播放| 真人做人爱边吃奶动态| 久久久久久久久久黄片| 亚洲精品国产一区二区精华液| 大型黄色视频在线免费观看| 91国产中文字幕| 日本熟妇午夜| 久久久精品欧美日韩精品| xxx96com| 欧美午夜高清在线| 91麻豆av在线| av天堂在线播放| www.www免费av| 欧美乱色亚洲激情| 91成年电影在线观看| 国产又黄又爽又无遮挡在线| 亚洲精品久久成人aⅴ小说| 国产麻豆成人av免费视频| 女同久久另类99精品国产91| 国产成人精品久久二区二区91| 中文字幕av电影在线播放| 在线永久观看黄色视频| 欧美另类亚洲清纯唯美| 国产精品久久久av美女十八| 老司机午夜十八禁免费视频| 午夜福利高清视频| 一本久久中文字幕| 久热爱精品视频在线9| av在线天堂中文字幕| 国内精品久久久久精免费| 国产精品香港三级国产av潘金莲| 欧美丝袜亚洲另类 | 成年版毛片免费区| 国产亚洲精品久久久久5区| 99精品在免费线老司机午夜| 国产精品免费一区二区三区在线| 亚洲国产精品久久男人天堂| 一a级毛片在线观看| 国内精品久久久久久久电影| 身体一侧抽搐| 精品久久久久久久末码| 色婷婷久久久亚洲欧美| 91av网站免费观看| 亚洲av中文字字幕乱码综合 | 亚洲av片天天在线观看| 久久久久久国产a免费观看| 日韩欧美国产在线观看| 久久婷婷人人爽人人干人人爱| 国产视频内射| 精品人妻1区二区| 亚洲成av人片免费观看| 国产视频一区二区在线看| 成人午夜高清在线视频 | 男女视频在线观看网站免费 | 波多野结衣av一区二区av| 日韩有码中文字幕| 午夜视频精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区三区视频了| a级毛片在线看网站| 亚洲一区二区三区不卡视频| 99国产综合亚洲精品| 精品国内亚洲2022精品成人| 久久天堂一区二区三区四区| 国产av一区在线观看免费| 日韩成人在线观看一区二区三区| bbb黄色大片| 欧美三级亚洲精品| 午夜免费激情av| 白带黄色成豆腐渣| 国产成人精品无人区| 在线观看舔阴道视频| 在线天堂中文资源库| xxx96com| 99精品欧美一区二区三区四区| 午夜久久久在线观看| 国产三级在线视频| 免费在线观看成人毛片| 久久欧美精品欧美久久欧美| 三级毛片av免费| 丝袜在线中文字幕| 男人操女人黄网站| 妹子高潮喷水视频| 久久久久九九精品影院| 夜夜看夜夜爽夜夜摸| 麻豆av在线久日| 久久久久国产精品人妻aⅴ院| 欧美日韩亚洲综合一区二区三区_| 久久 成人 亚洲| 免费高清视频大片| 国内精品久久久久精免费| 亚洲,欧美精品.| 在线播放国产精品三级| 久久亚洲精品不卡| 特大巨黑吊av在线直播 | 天堂影院成人在线观看| 一进一出好大好爽视频| 国产又黄又爽又无遮挡在线| 深夜精品福利| 啦啦啦观看免费观看视频高清| 1024香蕉在线观看| 法律面前人人平等表现在哪些方面| 欧美一级a爱片免费观看看 | 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 制服丝袜大香蕉在线| 国产成人精品久久二区二区免费| 波多野结衣高清作品| 日本一区二区免费在线视频| 国产午夜福利久久久久久| 欧美成人一区二区免费高清观看 | 亚洲成国产人片在线观看| 韩国精品一区二区三区| 亚洲欧美激情综合另类| 久久久国产欧美日韩av| 欧美一级毛片孕妇| 欧美在线黄色| 亚洲国产精品久久男人天堂| 亚洲激情在线av| √禁漫天堂资源中文www| 一区福利在线观看| 国产av一区在线观看免费| 久久亚洲精品不卡| 伦理电影免费视频| 久久久久久久久久黄片| 久久青草综合色| 一二三四社区在线视频社区8| 亚洲在线自拍视频| 国产黄a三级三级三级人| 精品国内亚洲2022精品成人| 久久久国产精品麻豆| 精品久久久久久久人妻蜜臀av| 天堂√8在线中文| 日韩大尺度精品在线看网址| 啦啦啦 在线观看视频| 精品国内亚洲2022精品成人| 亚洲av日韩精品久久久久久密| 国产激情欧美一区二区| 久久久久久大精品| av电影中文网址| 亚洲性夜色夜夜综合| 精品卡一卡二卡四卡免费| or卡值多少钱| 一本大道久久a久久精品| 国产一区二区三区在线臀色熟女| 欧美乱码精品一区二区三区| 一夜夜www| 757午夜福利合集在线观看| 欧美zozozo另类| 免费无遮挡裸体视频| 国产精品久久久久久人妻精品电影| 在线十欧美十亚洲十日本专区| 欧美成人性av电影在线观看| 老司机在亚洲福利影院| 久久久久九九精品影院| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆成人av免费视频| 亚洲av成人不卡在线观看播放网| 成人国语在线视频| 日韩精品免费视频一区二区三区| 日韩欧美国产在线观看| 一二三四社区在线视频社区8| 欧美最黄视频在线播放免费| 色综合站精品国产| 一本大道久久a久久精品| 亚洲欧美一区二区三区黑人| 757午夜福利合集在线观看| 国产精品久久久久久人妻精品电影| 国语自产精品视频在线第100页| 啪啪无遮挡十八禁网站| 亚洲五月色婷婷综合| 最近最新免费中文字幕在线| 国产午夜福利久久久久久| 国产欧美日韩一区二区精品| 国产精品影院久久| 99久久综合精品五月天人人| 日韩欧美一区视频在线观看| 国产aⅴ精品一区二区三区波| 免费av毛片视频| 午夜久久久久精精品| 黄片小视频在线播放| 国产成人av激情在线播放| 免费高清视频大片| 黑人操中国人逼视频| 久久中文字幕人妻熟女| 曰老女人黄片| 禁无遮挡网站| 自线自在国产av| 久久久久久亚洲精品国产蜜桃av| 久久精品夜夜夜夜夜久久蜜豆 | 看黄色毛片网站| 99国产综合亚洲精品| 亚洲国产精品合色在线| 高清毛片免费观看视频网站| 12—13女人毛片做爰片一| 无人区码免费观看不卡| 国产精品久久久av美女十八| 黄色丝袜av网址大全| 国产精品二区激情视频| 人人妻人人看人人澡| 免费观看人在逋| 99国产精品一区二区三区| 国产成人精品无人区| 亚洲精品一卡2卡三卡4卡5卡| 这个男人来自地球电影免费观看| 久久久久久免费高清国产稀缺| 最近在线观看免费完整版| 日本免费一区二区三区高清不卡| 热99re8久久精品国产| 97碰自拍视频| 可以免费在线观看a视频的电影网站| 欧美av亚洲av综合av国产av| 国产精品久久电影中文字幕| 欧美激情 高清一区二区三区| 高潮久久久久久久久久久不卡| 老汉色av国产亚洲站长工具| 一级毛片高清免费大全| 免费高清在线观看日韩| 美女国产高潮福利片在线看| 一本久久中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 成人国产一区最新在线观看| 首页视频小说图片口味搜索| 国产av一区二区精品久久| 制服人妻中文乱码| 黄色视频不卡| 亚洲午夜理论影院| 成人午夜高清在线视频 | 制服诱惑二区| 久久欧美精品欧美久久欧美| 女人被狂操c到高潮| av在线播放免费不卡| 国产成人啪精品午夜网站| 国产精品久久久人人做人人爽| 男女那种视频在线观看| 露出奶头的视频| 又黄又粗又硬又大视频| 免费搜索国产男女视频| 身体一侧抽搐| 99热6这里只有精品| 日韩成人在线观看一区二区三区| 亚洲男人的天堂狠狠| 国产精品免费一区二区三区在线| 国产精品野战在线观看| 黄频高清免费视频| 侵犯人妻中文字幕一二三四区| tocl精华| 人人妻人人澡欧美一区二区| 亚洲精品国产一区二区精华液| 午夜免费观看网址| 丁香六月欧美| 性色av乱码一区二区三区2| av免费在线观看网站| 大香蕉久久成人网| 成人国产一区最新在线观看| 亚洲欧美精品综合久久99| 欧美午夜高清在线| 免费一级毛片在线播放高清视频| 国产伦人伦偷精品视频| 亚洲专区国产一区二区| 一级作爱视频免费观看| 在线视频色国产色| www.熟女人妻精品国产| 国产午夜精品久久久久久| 成人永久免费在线观看视频| 亚洲五月天丁香| 妹子高潮喷水视频| av片东京热男人的天堂| 久久婷婷人人爽人人干人人爱| 午夜免费鲁丝| 国产精品av久久久久免费| 757午夜福利合集在线观看| 中国美女看黄片| 欧美一区二区精品小视频在线| 欧美丝袜亚洲另类 | 午夜久久久久精精品| 变态另类丝袜制服| 欧美 亚洲 国产 日韩一| 国产精品 国内视频| 国产精品综合久久久久久久免费| 亚洲免费av在线视频| 久久精品aⅴ一区二区三区四区| 亚洲在线自拍视频| 亚洲九九香蕉| 国产视频一区二区在线看| 热99re8久久精品国产| 精品国产一区二区三区四区第35| 亚洲欧美日韩高清在线视频| 国内毛片毛片毛片毛片毛片| 女同久久另类99精品国产91| 日本一区二区免费在线视频| 久久狼人影院| 这个男人来自地球电影免费观看| 婷婷精品国产亚洲av在线| 国产av一区在线观看免费| 90打野战视频偷拍视频| 亚洲专区国产一区二区| 老司机午夜十八禁免费视频| 1024香蕉在线观看| 美国免费a级毛片| 这个男人来自地球电影免费观看| 日本五十路高清| 中文在线观看免费www的网站 | 性欧美人与动物交配| 18禁黄网站禁片免费观看直播| 国产精品久久视频播放| 午夜久久久久精精品| 精品国产美女av久久久久小说| 久久久国产成人精品二区| 哪里可以看免费的av片| 免费观看精品视频网站| 亚洲无线在线观看| 女人爽到高潮嗷嗷叫在线视频| 免费高清在线观看日韩| 少妇熟女aⅴ在线视频| 久久精品91无色码中文字幕| 黄网站色视频无遮挡免费观看| 99在线视频只有这里精品首页| 老司机靠b影院| 在线观看日韩欧美| 日韩欧美国产一区二区入口| 欧美中文日本在线观看视频| 男女做爰动态图高潮gif福利片| www.熟女人妻精品国产| 99热只有精品国产| 50天的宝宝边吃奶边哭怎么回事| 19禁男女啪啪无遮挡网站| 脱女人内裤的视频| 国产一区二区激情短视频| 91九色精品人成在线观看| 亚洲精品中文字幕一二三四区| 又黄又爽又免费观看的视频| 久久亚洲真实| 一级黄色大片毛片| 欧美成狂野欧美在线观看| 欧美一区二区精品小视频在线| 亚洲欧美激情综合另类| 精品久久久久久,| 欧美另类亚洲清纯唯美| 亚洲激情在线av| 啦啦啦韩国在线观看视频| 成人亚洲精品一区在线观看| 欧美乱色亚洲激情| av超薄肉色丝袜交足视频| 亚洲人成网站在线播放欧美日韩| 久久香蕉激情| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色| 久热这里只有精品99| 亚洲午夜理论影院| 香蕉国产在线看| 久久精品夜夜夜夜夜久久蜜豆 | 成熟少妇高潮喷水视频| 亚洲av成人av| 成年版毛片免费区| 日韩一卡2卡3卡4卡2021年| 久久久久久久久久黄片| 女性被躁到高潮视频| 三级毛片av免费| 精品乱码久久久久久99久播| 女性被躁到高潮视频| 午夜激情福利司机影院| 亚洲av美国av| 成年版毛片免费区| 97超级碰碰碰精品色视频在线观看| 亚洲,欧美精品.| 亚洲一区中文字幕在线| 免费在线观看成人毛片| www.999成人在线观看| 亚洲色图av天堂| 国产区一区二久久| 变态另类丝袜制服| 在线观看免费视频日本深夜| 在线av久久热| 国产亚洲欧美在线一区二区| 久久香蕉精品热| 成熟少妇高潮喷水视频| 嫩草影院精品99| 国语自产精品视频在线第100页| 精品欧美国产一区二区三| 国产精品久久久久久亚洲av鲁大| 成熟少妇高潮喷水视频|