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

    強(qiáng)背景噪聲振動信號中滾動軸承故障沖擊特征提取

    2021-06-06 23:33:38劉湘楠趙學(xué)智上官文斌
    振動工程學(xué)報 2021年1期
    關(guān)鍵詞:奇異值分解滾動軸承故障診斷

    劉湘楠 趙學(xué)智 上官文斌

    摘要: 針對機(jī)械早期故障引起的沖擊特征微弱,易受強(qiáng)背景信號和噪聲的干擾而難以提取的問題,提出一種奇異值分解(Singular Value Decomposition, SVD)差分譜與S變換相結(jié)合的微弱沖擊特征提取方法。將原始信號構(gòu)造成Hankel矩陣,采用SVD對重構(gòu)矩陣進(jìn)行分解;利用奇異值差分譜確定降噪階次進(jìn)行降噪;采用S變換對降噪后的信號進(jìn)行時頻分析,提取信號中的微弱沖擊特征信息。通過數(shù)值仿真和實際軸承故障數(shù)據(jù)的對比,表明該方法可有效辨別軸承振動信號中故障引起的早期微弱沖擊特征,為軸承故障診斷提供先驗信息。

    關(guān)鍵詞: 故障診斷; 滾動軸承; 沖擊特征; 奇異值分解; S變換

    中圖分類號: TH165+.3; TH133.33; TN911.7??? 文獻(xiàn)標(biāo)志碼: A??? 文章編號: 1004-4523(2021)01-0202-09

    DOI:10.16385/j.cnki.issn.1004-4523.2021.01.023

    引 言

    滾動軸承是旋轉(zhuǎn)機(jī)械設(shè)備中廣泛使用的零部件之一,也是旋轉(zhuǎn)機(jī)械最主要的故障來源之一[1]。在實際運行過程中由于工況復(fù)雜、過載、安裝精度差及潤滑不良等原因,滾動軸承外圈、滾動體及內(nèi)圈等部件容易發(fā)生故障,進(jìn)而影響機(jī)械系統(tǒng)整體運行的安全性和可靠性[2]。滾動軸承產(chǎn)生故障時,工作過程中會產(chǎn)生突變的沖擊力,該沖擊信號的頻率即為損傷點撞擊軸承元件的頻率[3]。滾動軸承振動信號蘊(yùn)含了大量的運行狀態(tài)信息,表現(xiàn)為非線性和非平穩(wěn)性的調(diào)制信號,對振動信號進(jìn)行分析,可有效獲取振動信號中所包含的軸承故障引起的沖擊特征信息。但由于旋轉(zhuǎn)機(jī)械設(shè)備結(jié)構(gòu)復(fù)雜以及工作條件的多樣性,各種激勵源產(chǎn)生的信號相互耦合,導(dǎo)致滾動軸承故障引起的沖擊特征常常淹沒在強(qiáng)背景信號和噪聲中,比較難以識別,特別是早期故障信號,在強(qiáng)噪聲背景下軸承振動信號故障沖擊特征信息微弱,更加難以提取[4?5]。因此,如何實現(xiàn)強(qiáng)噪聲背景下軸承振動信號故障沖擊特征信息的有效提取,對旋轉(zhuǎn)機(jī)械設(shè)備正常運行具有重要意義。

    傳統(tǒng)的信號特征信息提取方法中,傅里葉變換由于缺乏時間、頻率的定位功能,在非平穩(wěn)信號分析中有很大的局限性。隨著短時傅里葉變換(STFT)、連續(xù)小波變換(CWT)等時頻分析理論的逐漸發(fā)展, 為快速精確提取非平穩(wěn)信號時頻特征提供了新的途徑。但現(xiàn)有的時頻分析方法都具有一定的局限性,如:STFT其時間窗口的大小、形狀都是固定的,與頻率無關(guān),無法同時獲得精準(zhǔn)的時刻和頻率[6]。CWT結(jié)果受Heisenberg不確定性原理和小波基函數(shù)的影響,在實際處理中難以找到合適的小波基函數(shù),而且CWT得到的小波系數(shù)譜會發(fā)生能量泄漏,使得瞬時頻率能量分布被模糊化[7]。S變換是一種新的時頻分析方法,克服了STFT窗函數(shù)以及CWT基函數(shù)固定不變的缺點,是一種高效的自適應(yīng)信號時頻分析方法,適合分析非線性和非平穩(wěn)信號,通過調(diào)整參數(shù)可以改變信號的時頻分辨率,以便得到?jīng)_擊信號在時頻域中良好的能量集中性,并且其逆變換完全無損[8?9]。Pinnegar等[10]研究表明S變換在無噪聲干擾下能夠顯著提高時頻分辨率,但在強(qiáng)背景噪聲時會出現(xiàn)時頻模糊的特點,難以有效識別信號中的沖擊特征信息。由于現(xiàn)場采集的軸承故障振動信號中往往存在著較強(qiáng)的環(huán)境噪聲及其他振動源干擾,會嚴(yán)重影響S變換的時頻分析精度,因此,在利用S變換對軸承振動信號故障沖擊特征信息進(jìn)行提取前,對信號進(jìn)行降噪處理是十分必要的。

    奇異值分解作為一種非線性信號降噪方法,在軸承振動信號[11?12]、語音信號[13]、電荷放電信號[14]等不同性質(zhì)信號降噪方面有著廣泛的應(yīng)用。利用SVD對信號進(jìn)行降噪的關(guān)鍵在于如何利用信號構(gòu)造合適的矩陣以及確定有效奇異值的個數(shù)。目前SVD降噪中應(yīng)用最廣泛的矩陣形式是Hankel矩陣,采用Hankel矩陣的突出特點就是能夠消除信號中的噪聲[15?17]。另外,信號非零奇異值中有效奇異值個數(shù)的選取決定著降噪效果的好壞,奇異值個數(shù)選取過少容易導(dǎo)致信號有用信息丟失,有時甚至?xí)?dǎo)致信號波形出現(xiàn)畸變,而選取過多又會造成信號降噪效果較差[18]。趙學(xué)智等[19]研究表明,根據(jù)奇異值差分譜的峰值位置可以準(zhǔn)確地確定有效奇異值的個數(shù)。

    基于以上原因,本文提出一種基于SVD差分譜和S變換相結(jié)合的微弱故障沖擊特征信息提取方法。該方法將SVD差分譜作為S變換的前置濾波單元,對信號濾波降噪后再進(jìn)行S變換時頻分析,根據(jù)信號時頻譜圖獲取故障沖擊特征。最后,將該方法應(yīng)用于軸承故障信號分析,成功提取出其故障沖擊特征,獲得了較為滿意的結(jié)果。

    1 理論分析

    3.1 內(nèi)圈故障分析

    結(jié)合旋轉(zhuǎn)機(jī)械設(shè)備的實際工況,滾動軸承故障振動信號通常易受強(qiáng)背景噪聲干擾,因此,本文在軸承內(nèi)圈故障振動信號中添加信噪比為0的高斯白噪聲,以模擬工程實際。圖9為加噪前后軸承內(nèi)圈故障振動信號波形圖。

    由圖9(b)可知,加噪后軸承內(nèi)圈故障引起的沖擊特征微弱,從時域信號波形圖中已無法獲取沖擊特征。本文采用STFT,CWT,S變換分別對加噪前后軸承內(nèi)圈故障振動信號進(jìn)行時頻分析,其中CWT所用的小波為Morlet小波,尺度為64。

    圖10為內(nèi)圈故障振動信號加噪前后時頻譜圖。由圖10(b)可知,當(dāng)軸承內(nèi)圈故障引起的沖擊特征呈現(xiàn)出強(qiáng)噪聲背景下特征微弱時,采用STFT, CWT和S變換3種時頻分析方法均難以有效獲取內(nèi)圈故障引起的沖擊特征。為準(zhǔn)確識別強(qiáng)噪聲背景下軸承內(nèi)圈故障引起的沖擊特征,需對信號進(jìn)行降噪預(yù)處理。

    將加噪的軸承內(nèi)圈故障振動信號構(gòu)造成512×513階Hankel矩陣,采用SVD對矩陣進(jìn)行分解,根據(jù)奇異值差分譜選取有效奇異值。 圖11為信號的奇異值差分譜圖,將28個有效奇異值對應(yīng)的信號組合重構(gòu),獲得降噪后的軸承內(nèi)圈故障振動信號。

    采用STFT,CWT,SST三種方法對SVD差分譜降噪后的軸承內(nèi)圈故障振動信號進(jìn)行時頻分析,圖12為內(nèi)圈故障振動信號降噪后時頻譜圖。

    由圖12可知,利用本文所提出的方法對降噪后的內(nèi)圈故障振動信號進(jìn)行時頻分析,可有效獲取微弱沖擊特征。整個時間軸上,頻率軸約3000 Hz處出現(xiàn)明顯的周期性沖擊特性,其周期約為0.006253 s,對應(yīng)的頻率為159.92 Hz,與內(nèi)圈故障頻率157.94 Hz基本一致。

    對采用SVD差分譜降噪后的內(nèi)圈故障振動信號時頻譜進(jìn)行S逆變換,提取軸承內(nèi)圈故障振動信號中的時域沖擊特征,如圖13所示。

    對比圖9(a)和13可知,本方法所提取的內(nèi)圈故障引起的沖擊特征不可避免地出現(xiàn)一定的變形和失真,但作為最重要信息的沖擊特征出現(xiàn)頻率,則可以完全有效地提取出來。

    3.2 外圈故障分析

    將本文所提出的方法應(yīng)用于滾動軸承外圈單點故障振動信號特征提取,圖14為軸承外圈故障振動信號波形圖。

    由圖14可知,軸承外圈故障引起的沖擊特征在強(qiáng)噪聲背景下特征微弱,由信號時域圖難以獲取信號沖擊特征。

    將外圈故障振動信號構(gòu)造成512×513階Hankel矩陣,采用SVD對矩陣進(jìn)行分解,根據(jù)奇異值差分譜選取有效奇異值。圖15為信號的奇異值差分譜圖,將14個有效奇異值對應(yīng)的信號組合重構(gòu),獲得降噪后的外圈故障振動信號。

    采用S變換對降噪后外圈故障振動信號進(jìn)行時頻分析。圖16為降噪后滾動體故障振動信號的S變換時頻譜圖。

    由圖16可知,采用S變換對降噪后的外圈故障振動信號進(jìn)行時頻分析,可有效獲取微弱沖擊特征。整個時間軸上,頻率軸約3500 Hz處出現(xiàn)明顯的周期性沖擊特性,其周期約為0.009667 s,對應(yīng)的頻率為103.44 Hz,與外圈故障頻率104.57 Hz基本一致。

    對采用SVD差分譜降噪后的時頻譜進(jìn)行S逆變換,提取軸承外圈故障振動信號中的時域沖擊特征,如圖17所示。

    由圖17可知,采用本文方法能實現(xiàn)強(qiáng)噪聲背景下軸承外圈故障引起的沖擊特征信息。

    綜上所述,實驗?zāi)M的軸承內(nèi)圈分析結(jié)果驗證了本文方法對強(qiáng)噪聲背景下故障微弱沖擊特征提取的有效性。在此基礎(chǔ)上,將本文所提出的方法應(yīng)用于軸承外圈故障振動信號進(jìn)行分析,結(jié)果表明,該方法能夠有效提取外圈故障引起的沖擊特征,實現(xiàn)了對滾動軸承的故障診斷。

    4 結(jié) 論

    (1) 提出了一種基于SVD差分譜和SST相結(jié)合的信號弱沖擊特征提取方法。首先將原始信號構(gòu)造成Hankel矩陣,利用SVD對矩陣進(jìn)行分解;再利用差分譜理論確定降噪階次進(jìn)行降噪,最后,采用S變換對降噪后的信號進(jìn)行時頻分析,實現(xiàn)了信號弱沖擊特征的有效提取;

    (2) 對比分析了STFT,CWT,S變換3種方法對信號沖擊特征的識別能力,結(jié)果表明: STFT無法識別信號弱沖擊特征;CWT能識別信號弱沖擊特征,但其分辨率及時頻聚集性不如S變換;S變換提高了時頻聚集性,能夠有效識別信號中的沖擊特征。

    參考文獻(xiàn):

    [1]??????? El-Thalji I, Jantunen E. A summary of fault modelling and predictive health monitoring of rolling element bearings[J]. Mechanical Systems & Signal Processing, 2015,60-61(1): 252-272.

    [2]??????? Dyba?a J. Rolling bearing diagnosing method based on empirical mode decomposition of machine vibration signal[J]. Applied Acoustics, 2014, 77(3):195-203.

    [3]??????? 唐貴基, 蔡 偉. 應(yīng)用小波包和包絡(luò)分析的滾動軸承故障 診斷[J]. 振動、測試與診斷, 2009, 29(2): 201-204.

    Tang Guiji, Cai Wei. Rolling bearings fault diagnosis by using wavelet packet and envelope analysis[J]. Journal of Vibration, Measurement & Diagnosis, 2009, 29(2): 201-204.

    [4]??????? 陳志新, 徐金梧, 楊德斌. 基于復(fù)小波塊閾值的降噪方 法及其在機(jī)械故障診斷中的應(yīng)用[J]. 機(jī)械工程學(xué)報, 2007, 43(6):200-204.

    Chen Zhixin, Xu Jinwu, Yang Debin. Denoising method of block thresholding based on DT-CWT and its applocation in mechanical fault diagnosis[J]. Journal of Mechanical Engineering, 2007, 43(6):200-204.

    [5]??????? 李富才, 何正嘉, 陳 進(jìn). 小波域相關(guān)濾波法及其早期故 障預(yù)示應(yīng)用[J]. 振動工程學(xué)報, 2005,18(2):145-148.

    Li Fucai, He Zhengjia, Chen Jin. Wavelet transform doma domain correlation filter and its application in incipient fault prognosis [J]. Journal of Vibration Engineering, 2005,18(2):145-148.

    [6]??????? 劉義亞, 李 可, 陳 鵬. 基于同步壓縮小波變換的滾動軸承故障診斷[J]. 中國機(jī)械工程, 2018,29(5): 585-590.

    Liu Yiya, Li Ke, Chen Peng. Fault diagnosis for rolling bearings based on synchrosqueezing wavelet transform[J]. China Mechanical Engineering, 2018, 29(5): 585-590.

    [7]??????? Sinha S, Routh P S, Anno P D, et al. Spectral decomposition of seismic data with continuous-wavelet transform[J]. Geophysics, 2008, 70(6):19-25.

    [8]??????? Pinnegar C R , Mansinha L . The S-transform with windows of arbitrary and varying shape[J]. Geophysics, 2003, 68(1):381-385.

    [9]??????? Kazemi K, Amirian M, Dehghani M J. The S-transform using a new window to improve frequency and time resolutions[J]. Signal, Image and Video Processing, 2014, 8(3):533-541.

    [10]????? Pinnegar C R, Mansinha L. Time-local spectral analysis for non-stationary time series: The S-transform for noisy signals[J]. Fluctuation & Noise Letters, 2012,3(3):L357-L364.

    [11]????? 王建國, 李 健, 萬旭東. 基于奇異值分解和局域均值 分解的滾動軸承故障特征提取方法[J]. 機(jī)械工程學(xué)報,2015,(3):104-110.

    Wang Jianguo, Li Jian, Wang Xudong. Fault feature extraction method of rolling bearings based on singular value decomposition and local mean decomposition [J]. Journal of Mechanical Engineering, 2015,(3):104-110.

    [12]????? Golafshan R, Sanliturk K Y. SVD and Hankel matrix based de-noising approach for ball bearing fault detection and its assessment using artificial faults[J]. Mechanical Systems & Signal Processing, 2016, (70?71):36-50.

    [13]????? Zheng C, Peng R, Li J, et al. A constrained MMSE LP residual estimator for speech dereverberation in noisy environments[J]. IEEE Signal Processing Letters, 2014, 21(12):1462-1466.

    [14]????? Ashtiani M B, Shahrtash S M. Partial discharge de-noising employing adaptive singular value decomposition[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2014, 21(2):775-782.

    [15]????? 李 華, 劉 韜, 伍 星,等. 基于SVD和熵優(yōu)化頻帶熵的滾 動軸承故障診斷研究[J]. 振動工程學(xué)報, 2018, 31 (2):180-186.

    Li Hua, Liu Tao,Wu Xing, et al.Research on fault diagnosis of rolling bearing based on SVD and optimized frequency band entropy by entropy[J]. Journal of Vibration Engineering,2018,31(2):180-186.

    [16]????? 胥永剛,張志新,馬朝永,等. 改進(jìn)奇異譜分解及其在軸承故障診斷中的應(yīng)用[J]. 振動工程學(xué)報, 2019, 32 (3): 540-547.

    Xu Yonggang, Zhang Zhixin,Ma Chaoyong, et al. Improved singular spectrum decomposition and its application in bearing fault diagnosis[J],Journal of Vibration Engineering,2019, 32 (3): 540-547.

    [17]????? 王 超, 孔凡讓, 黃偉國,等. 改進(jìn)的奇異值分解在軸承故障診斷中的應(yīng)用[J]. 振動工程學(xué)報, 2014, 27(2): 296-303.

    Wang Chao, Kong Fanrang, Huang Weiguo,et al. Application of improved singular value decomposition in bearing fault diagnosis[J].Journal of Vibration Engineering,2014, 27(2): 296-303.

    [18]????? 趙學(xué)智, 葉邦彥, 陳統(tǒng)堅. 奇異值差分譜理論及其在 車床主軸箱故障診斷中的應(yīng)用[J]. 機(jī)械工程學(xué)報, 2010, 46(1):100-108.

    Zhao Xuezhi, Ye Bangyan, Chen Tongjian. Difference spectrum theory of singular value and its application to the fault diagnosis of headstock of lathe[J]. Journal of Mechanical Engineering, 2010, 46(1):100-108.

    [19]????? 趙學(xué)智, 葉邦彥, 陳統(tǒng)堅. 基于小波—奇異值分解差 分譜的弱故障特征提取方法[J]. 機(jī)械工程學(xué)報, 2012, 48(7):37-48.

    Zhao Xuezhi, Ye Bangyan, Chen Tongjian. Extraction method of faint fault feature based on wavelet-SVD difference spectrum[J]. Journal of Mechanical Engineering, 2012, 48(7):37-48.

    [20]????? 呂志民, 張武軍, 徐金梧, 等. 基于奇異譜的降噪方 法及其在故障診斷技術(shù)中的應(yīng)用[J]. 機(jī)械工程學(xué)報, 1999, 35(3):85-88.

    Lü Zhiming, Zhang Wujun, Xu Jinwu, et al. A noise reduction method based singular spectrum and its application in machine fault diagnosis[J]. Journal of Mechanical Engineering, 1999, 35(3):85-88.

    [21]????? 楊文獻(xiàn), 姜節(jié)勝. 機(jī)械信號奇異熵研究[J]. 機(jī)械工程 學(xué)報, 2000, 36(12):9-13.

    Yang Wenxian, Jiang Jiesheng. Research on mechanical signal singular entropy[J]. Journal of Mechanical Engineering, 2000, 36(12):9-13.

    [22]????? 何 田, 劉獻(xiàn)棟, 李其漢. 噪聲背景下檢測突變信息的 奇異值分解技術(shù)[J]. 振動工程學(xué)報, 2006, 19(3):399-403.

    He Tian, Liu Xiandong, Li Qihan. An improved method singularity value of detecting abrupt information based on decomposition in noise background [J]. Journal of Vibration Engineering, 2006, 19(3): 399-403.

    [23]????? Stockwell R G, Mansinha L, Lowe R P. Localization of the complex spectrum: The S transform[J]. IEEE Transactions on Signal Processing, 2002, 44(4):998-1001.

    Abstract: Aiming at the problem that the impact of early mechanical failure is weak and it is difficult to extract due to strong background signal and noise interference, a weak impact feature extraction method combining singular value decomposition (SVD) differential spectrum and S?transform is proposed. The original signal is constructed into a Hankel matrix, and the reconstruction matrix is decomposed by SVD. The noise reduction order is determined by singular value difference spectrum for noise reduction. The S-transformation is used to analyze the time-frequency of the denoised signal, the weak impact characteristic information of the signal is extracted. The comparison between numerical simulation and actual bearing fault data shows that the method can effectively distinguish the early weak shock characteristics caused by faults in the bearing vibration signal. It can provide a priori information for the bearing fault diagnosis.

    Key words: fault diagnosis; rolling bearing; impact feature; singular value decomposition; S?transform

    作者簡介: 劉湘楠(1992?),男,博士研究生。電話:16607319453;E-mail:lxn920613@163.com

    猜你喜歡
    奇異值分解滾動軸承故障診斷
    結(jié)合PCA及字典學(xué)習(xí)的高光譜圖像自適應(yīng)去噪方法
    滾動軸承裝配方法及注意事項
    滾動軸承安裝使用中的幾個問題
    基于分塊DWT和SVD的魯棒性數(shù)字水印算法
    一種基于奇異值分解的魯棒水印算法
    基于小波包絡(luò)分析的滾動軸承典型故障診斷技術(shù)研究
    商情(2016年11期)2016-04-15 07:39:54
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    基于包絡(luò)分析法的軸承故障分析
    協(xié)同過濾算法改進(jìn)及研究
    基于LCD和排列熵的滾動軸承故障診斷
    久久亚洲精品不卡| 亚洲成人中文字幕在线播放| 中国美白少妇内射xxxbb| 亚洲av成人精品一区久久| 内地一区二区视频在线| 精品午夜福利视频在线观看一区| 好男人在线观看高清免费视频| 久久久久国内视频| 在线观看午夜福利视频| 欧美黑人巨大hd| 欧美性感艳星| 国产精品免费一区二区三区在线| 国产精品av视频在线免费观看| 免费看a级黄色片| 日韩一本色道免费dvd| 91久久精品电影网| 在线天堂最新版资源| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品合色在线| 久久久久久久久久成人| 国产成人a区在线观看| 黄色丝袜av网址大全| 日韩 亚洲 欧美在线| 色精品久久人妻99蜜桃| 成人午夜高清在线视频| 国产一区二区在线av高清观看| 亚洲av美国av| 国产精品一及| 男人舔女人下体高潮全视频| 亚洲精品一区av在线观看| 老熟妇仑乱视频hdxx| 精品一区二区三区人妻视频| 亚洲一区二区三区色噜噜| 国语自产精品视频在线第100页| 少妇熟女aⅴ在线视频| 成人永久免费在线观看视频| 国产成人一区二区在线| 久久久久久久久久成人| 内射极品少妇av片p| 特级一级黄色大片| 两性午夜刺激爽爽歪歪视频在线观看| 中文字幕av在线有码专区| 精品一区二区三区av网在线观看| 日韩欧美国产在线观看| 99久久成人亚洲精品观看| 校园春色视频在线观看| 免费看av在线观看网站| 成人二区视频| 又爽又黄无遮挡网站| 亚洲精品456在线播放app | 91精品国产九色| 婷婷丁香在线五月| 亚洲男人的天堂狠狠| 又爽又黄无遮挡网站| 如何舔出高潮| 精品午夜福利视频在线观看一区| 久久精品国产99精品国产亚洲性色| 国产单亲对白刺激| 女人被狂操c到高潮| 欧美色欧美亚洲另类二区| 欧美色欧美亚洲另类二区| 啦啦啦啦在线视频资源| 欧美极品一区二区三区四区| 国产精品久久电影中文字幕| 亚洲avbb在线观看| 成人亚洲精品av一区二区| 精品午夜福利视频在线观看一区| 成熟少妇高潮喷水视频| 欧美日韩亚洲国产一区二区在线观看| 最新中文字幕久久久久| 久久国产乱子免费精品| 五月伊人婷婷丁香| 少妇被粗大猛烈的视频| 99久久精品热视频| 久久精品夜夜夜夜夜久久蜜豆| 精品一区二区免费观看| 免费电影在线观看免费观看| 国产精品国产三级国产av玫瑰| 久久精品影院6| 亚洲国产欧美人成| 有码 亚洲区| 女生性感内裤真人,穿戴方法视频| 床上黄色一级片| 日本精品一区二区三区蜜桃| 国产国拍精品亚洲av在线观看| 亚洲精品粉嫩美女一区| 舔av片在线| 91在线精品国自产拍蜜月| 春色校园在线视频观看| 国产欧美日韩一区二区精品| 国内精品宾馆在线| 日本 欧美在线| 久久久国产成人免费| 韩国av在线不卡| 色视频www国产| 国产精品伦人一区二区| 精品久久久久久成人av| 国产精品一区二区免费欧美| 精品国内亚洲2022精品成人| 在线观看舔阴道视频| 日韩,欧美,国产一区二区三区 | 亚洲最大成人中文| 国产精品国产高清国产av| 又紧又爽又黄一区二区| 国产一区二区三区视频了| 亚洲av中文av极速乱 | 亚洲男人的天堂狠狠| 成人综合一区亚洲| 成人欧美大片| 欧美激情国产日韩精品一区| .国产精品久久| 三级国产精品欧美在线观看| 国产精品久久久久久亚洲av鲁大| 国产午夜精品久久久久久一区二区三区 | 露出奶头的视频| 日日干狠狠操夜夜爽| 春色校园在线视频观看| 在线免费十八禁| x7x7x7水蜜桃| 一级a爱片免费观看的视频| 亚洲真实伦在线观看| 成人三级黄色视频| 国产免费av片在线观看野外av| 天堂动漫精品| 超碰av人人做人人爽久久| 少妇猛男粗大的猛烈进出视频 | 久久精品国产鲁丝片午夜精品 | 久久精品国产亚洲av香蕉五月| 99热只有精品国产| 成人特级av手机在线观看| 俄罗斯特黄特色一大片| 亚洲国产精品成人综合色| 男女视频在线观看网站免费| a在线观看视频网站| 三级毛片av免费| 精品日产1卡2卡| 午夜精品一区二区三区免费看| 国产精品不卡视频一区二区| 婷婷精品国产亚洲av| 一进一出抽搐gif免费好疼| 人人妻,人人澡人人爽秒播| 在线观看一区二区三区| 亚洲自偷自拍三级| 深爱激情五月婷婷| 日韩高清综合在线| 两人在一起打扑克的视频| 老师上课跳d突然被开到最大视频| 91在线观看av| 亚洲av五月六月丁香网| 三级国产精品欧美在线观看| 亚洲成人精品中文字幕电影| www.色视频.com| 床上黄色一级片| 校园人妻丝袜中文字幕| 亚洲欧美日韩高清专用| 两人在一起打扑克的视频| 一卡2卡三卡四卡精品乱码亚洲| 精品一区二区三区av网在线观看| 中国美白少妇内射xxxbb| 国产高清不卡午夜福利| 亚洲专区中文字幕在线| 久久久精品大字幕| 老女人水多毛片| 美女被艹到高潮喷水动态| 免费看光身美女| 日韩中文字幕欧美一区二区| 欧美+亚洲+日韩+国产| 精华霜和精华液先用哪个| 久久中文看片网| 国产精品三级大全| 日韩av在线大香蕉| 亚洲av免费高清在线观看| 成年女人毛片免费观看观看9| 国产成人a区在线观看| 日韩在线高清观看一区二区三区 | avwww免费| 国产一区二区在线av高清观看| 中文字幕熟女人妻在线| videossex国产| 国产人妻一区二区三区在| 老司机午夜福利在线观看视频| 国产精品一区二区免费欧美| 国产免费男女视频| 久久热精品热| eeuss影院久久| 九色成人免费人妻av| 日本爱情动作片www.在线观看 | 午夜老司机福利剧场| 成人性生交大片免费视频hd| 久久草成人影院| 观看免费一级毛片| 中文字幕熟女人妻在线| 长腿黑丝高跟| 女人十人毛片免费观看3o分钟| 日韩 亚洲 欧美在线| 亚洲图色成人| 亚洲性久久影院| 一进一出抽搐gif免费好疼| 日本 av在线| 国产精品爽爽va在线观看网站| 国产精品嫩草影院av在线观看 | 天堂av国产一区二区熟女人妻| 国产黄片美女视频| 又粗又爽又猛毛片免费看| 国产精品av视频在线免费观看| 乱人视频在线观看| av女优亚洲男人天堂| 亚洲男人的天堂狠狠| 久久久久久久久久成人| 悠悠久久av| 久久亚洲精品不卡| 黄色视频,在线免费观看| 免费观看的影片在线观看| 99久国产av精品| 在线国产一区二区在线| 精品一区二区三区人妻视频| 亚洲乱码一区二区免费版| 波多野结衣高清作品| 久久午夜福利片| 91精品国产九色| 国产亚洲精品综合一区在线观看| 又黄又爽又免费观看的视频| 日日啪夜夜撸| 免费电影在线观看免费观看| 99国产精品一区二区蜜桃av| 乱系列少妇在线播放| 日本在线视频免费播放| 久久久久久久久中文| 在线观看66精品国产| av天堂在线播放| 99热网站在线观看| 国产一区二区三区在线臀色熟女| 午夜免费激情av| 三级毛片av免费| 成人综合一区亚洲| 在线观看免费视频日本深夜| 亚洲成人精品中文字幕电影| 精品久久久久久久久av| 久久久久九九精品影院| 精品午夜福利视频在线观看一区| 国产精品电影一区二区三区| 欧美区成人在线视频| 自拍偷自拍亚洲精品老妇| 天堂动漫精品| 亚洲va在线va天堂va国产| 成人亚洲精品av一区二区| 国产伦在线观看视频一区| 欧美一区二区精品小视频在线| 97超视频在线观看视频| 最后的刺客免费高清国语| 日本在线视频免费播放| 国产精品一区二区三区四区免费观看 | 成熟少妇高潮喷水视频| 天天一区二区日本电影三级| 精品久久久久久,| 欧美日韩乱码在线| 国产单亲对白刺激| 色5月婷婷丁香| 久久人人精品亚洲av| 九九热线精品视视频播放| 国产精品美女特级片免费视频播放器| 最近最新中文字幕大全电影3| 国产精品亚洲一级av第二区| 成人国产麻豆网| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 精品福利观看| 亚洲美女视频黄频| 成人无遮挡网站| 亚洲美女搞黄在线观看 | 网址你懂的国产日韩在线| 成人国产麻豆网| 午夜福利在线观看吧| 俺也久久电影网| 真实男女啪啪啪动态图| 麻豆久久精品国产亚洲av| 九色成人免费人妻av| xxxwww97欧美| 国内精品一区二区在线观看| 国产伦精品一区二区三区视频9| 美女cb高潮喷水在线观看| 亚洲三级黄色毛片| 国产男靠女视频免费网站| 免费搜索国产男女视频| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| 免费看av在线观看网站| 2021天堂中文幕一二区在线观| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 国产伦人伦偷精品视频| 18+在线观看网站| 深夜精品福利| 亚洲天堂国产精品一区在线| 久久亚洲精品不卡| 亚洲avbb在线观看| 精品人妻视频免费看| 成人性生交大片免费视频hd| 日本-黄色视频高清免费观看| 丝袜美腿在线中文| 亚洲欧美日韩东京热| 成年女人永久免费观看视频| 亚洲人成网站高清观看| 国产主播在线观看一区二区| 三级毛片av免费| 少妇高潮的动态图| 免费观看人在逋| 窝窝影院91人妻| 我的女老师完整版在线观看| 亚洲第一电影网av| 一a级毛片在线观看| 桃色一区二区三区在线观看| 一进一出抽搐gif免费好疼| 亚洲va在线va天堂va国产| 亚洲国产精品合色在线| 亚洲av第一区精品v没综合| 麻豆一二三区av精品| 国产私拍福利视频在线观看| 国产精品久久电影中文字幕| 国产不卡一卡二| 色视频www国产| 大型黄色视频在线免费观看| 久99久视频精品免费| 亚洲一区二区三区色噜噜| 久久久久国内视频| 国产精品久久久久久久电影| 禁无遮挡网站| 国内精品美女久久久久久| 18禁裸乳无遮挡免费网站照片| 国产黄片美女视频| 噜噜噜噜噜久久久久久91| 人妻制服诱惑在线中文字幕| 久久精品综合一区二区三区| 成人综合一区亚洲| 久久精品影院6| 久久国产精品人妻蜜桃| 色综合亚洲欧美另类图片| 内地一区二区视频在线| 国内揄拍国产精品人妻在线| 亚洲国产精品久久男人天堂| 精品一区二区免费观看| 国产大屁股一区二区在线视频| 久久人人精品亚洲av| 18禁裸乳无遮挡免费网站照片| 久久久久久久久中文| 国产色爽女视频免费观看| 特级一级黄色大片| 国产 一区精品| 性色avwww在线观看| 人妻久久中文字幕网| 日韩国内少妇激情av| 波多野结衣高清无吗| 国产精品免费一区二区三区在线| 午夜福利在线观看吧| 免费人成在线观看视频色| 亚洲成人久久性| 亚洲内射少妇av| 国产中年淑女户外野战色| 亚洲中文日韩欧美视频| 国内精品久久久久精免费| av在线蜜桃| 嫩草影院新地址| 国产在线精品亚洲第一网站| 亚洲 国产 在线| 日日干狠狠操夜夜爽| 真实男女啪啪啪动态图| 男人舔奶头视频| 天美传媒精品一区二区| av黄色大香蕉| 男人舔奶头视频| 性插视频无遮挡在线免费观看| 中文字幕久久专区| 色噜噜av男人的天堂激情| 亚洲四区av| 波野结衣二区三区在线| www.www免费av| 伦精品一区二区三区| 深爱激情五月婷婷| 免费看光身美女| 国产精品女同一区二区软件 | 老熟妇仑乱视频hdxx| 亚洲av免费高清在线观看| av.在线天堂| 久久天躁狠狠躁夜夜2o2o| 免费av观看视频| 九九热线精品视视频播放| 午夜视频国产福利| 成年女人永久免费观看视频| 真人做人爱边吃奶动态| 我要看日韩黄色一级片| 国产精品乱码一区二三区的特点| 亚洲aⅴ乱码一区二区在线播放| 日韩精品中文字幕看吧| 久久久久国内视频| 精品一区二区三区视频在线观看免费| 在线免费十八禁| ponron亚洲| 精品人妻熟女av久视频| 日本与韩国留学比较| xxxwww97欧美| av国产免费在线观看| 久久精品国产99精品国产亚洲性色| 国产一区二区三区av在线 | 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| 国产美女午夜福利| 精品人妻1区二区| 久久精品国产亚洲网站| 神马国产精品三级电影在线观看| 麻豆成人午夜福利视频| 国内精品美女久久久久久| 亚洲五月天丁香| 一个人看的www免费观看视频| 欧美在线一区亚洲| 日本熟妇午夜| 免费人成视频x8x8入口观看| 婷婷精品国产亚洲av| 嫩草影院入口| 夜夜看夜夜爽夜夜摸| 亚洲av免费高清在线观看| 亚洲中文字幕日韩| 91在线精品国自产拍蜜月| 天堂网av新在线| 精品人妻偷拍中文字幕| 欧美zozozo另类| 国产亚洲精品av在线| 欧美性猛交黑人性爽| 赤兔流量卡办理| 九九热线精品视视频播放| 中国美白少妇内射xxxbb| 校园人妻丝袜中文字幕| 成年版毛片免费区| 国内精品宾馆在线| 午夜福利在线在线| 老司机午夜福利在线观看视频| 精品人妻视频免费看| 人妻久久中文字幕网| 九九久久精品国产亚洲av麻豆| 九色国产91popny在线| 国产精品人妻久久久影院| 联通29元200g的流量卡| 久久久久久大精品| 国产精品不卡视频一区二区| 超碰av人人做人人爽久久| 色综合婷婷激情| 亚洲av成人av| 99久久精品一区二区三区| 噜噜噜噜噜久久久久久91| 亚洲电影在线观看av| 日韩精品有码人妻一区| 国产成人av教育| a级一级毛片免费在线观看| 亚洲成人精品中文字幕电影| 美女cb高潮喷水在线观看| 国产精品国产高清国产av| 小蜜桃在线观看免费完整版高清| 91久久精品电影网| 毛片一级片免费看久久久久 | 在线天堂最新版资源| 嫩草影院精品99| 亚洲图色成人| 97碰自拍视频| 国内毛片毛片毛片毛片毛片| 亚洲图色成人| 国产精品综合久久久久久久免费| 午夜老司机福利剧场| 亚洲性久久影院| 少妇的逼好多水| 成人亚洲精品av一区二区| 久久精品人妻少妇| 欧美成人a在线观看| 特级一级黄色大片| 变态另类丝袜制服| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美精品v在线| 午夜爱爱视频在线播放| 波野结衣二区三区在线| 观看美女的网站| 欧美高清成人免费视频www| 亚洲av美国av| 久久久久久久久久久丰满 | 国产91精品成人一区二区三区| 老司机午夜福利在线观看视频| 国产伦精品一区二区三区视频9| 亚洲欧美日韩高清在线视频| 国产精品精品国产色婷婷| 日韩精品中文字幕看吧| 1024手机看黄色片| 国产精品无大码| 精品人妻一区二区三区麻豆 | 免费人成视频x8x8入口观看| 1000部很黄的大片| 免费看日本二区| 色吧在线观看| 午夜日韩欧美国产| 嫩草影院入口| 亚洲专区中文字幕在线| 精品久久久久久成人av| 少妇的逼好多水| 国产探花极品一区二区| 国产亚洲av嫩草精品影院| 99国产极品粉嫩在线观看| 日韩欧美在线二视频| 国产成人福利小说| 久久久久国产精品人妻aⅴ院| 亚洲最大成人手机在线| 亚洲精品一区av在线观看| 亚洲国产精品久久男人天堂| 男女边吃奶边做爰视频| 国产精品乱码一区二三区的特点| 午夜精品在线福利| 成人三级黄色视频| 精品久久久噜噜| 午夜亚洲福利在线播放| 亚洲中文字幕日韩| 麻豆av噜噜一区二区三区| 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 久久天躁狠狠躁夜夜2o2o| 又粗又爽又猛毛片免费看| 欧美另类亚洲清纯唯美| 色哟哟·www| 看免费成人av毛片| av女优亚洲男人天堂| 五月伊人婷婷丁香| 成人一区二区视频在线观看| 99国产精品一区二区蜜桃av| 老熟妇仑乱视频hdxx| 亚洲av不卡在线观看| 丰满人妻一区二区三区视频av| 亚洲国产日韩欧美精品在线观看| 色播亚洲综合网| 永久网站在线| 成年女人毛片免费观看观看9| 日本精品一区二区三区蜜桃| 国产免费男女视频| 好男人在线观看高清免费视频| 97人妻精品一区二区三区麻豆| 久久久久久久亚洲中文字幕| 校园人妻丝袜中文字幕| 人妻夜夜爽99麻豆av| 亚洲av不卡在线观看| 床上黄色一级片| 白带黄色成豆腐渣| 亚洲人与动物交配视频| 国产伦精品一区二区三区视频9| 中文字幕高清在线视频| 一边摸一边抽搐一进一小说| 免费看a级黄色片| 老女人水多毛片| 制服丝袜大香蕉在线| 国产极品精品免费视频能看的| 亚洲精品456在线播放app | 内射极品少妇av片p| 亚洲四区av| 亚洲精品成人久久久久久| 午夜精品一区二区三区免费看| 岛国在线免费视频观看| 亚洲专区中文字幕在线| 黄色丝袜av网址大全| 99九九线精品视频在线观看视频| 免费观看在线日韩| 亚洲欧美精品综合久久99| 哪里可以看免费的av片| 成人国产综合亚洲| 美女 人体艺术 gogo| www.色视频.com| 亚洲性久久影院| 1024手机看黄色片| 亚洲美女视频黄频| 在线观看舔阴道视频| 国内精品宾馆在线| 最新中文字幕久久久久| 欧美高清成人免费视频www| 桃红色精品国产亚洲av| 美女免费视频网站| 国产av麻豆久久久久久久| 一区二区三区免费毛片| 日本与韩国留学比较| 午夜精品一区二区三区免费看| 欧美激情在线99| 国产三级中文精品| 国产单亲对白刺激| 草草在线视频免费看| 十八禁国产超污无遮挡网站| 91麻豆精品激情在线观看国产| 精品久久久久久久久久久久久| 高清日韩中文字幕在线| 亚洲在线自拍视频| 免费电影在线观看免费观看| 久久国产精品人妻蜜桃| 91麻豆精品激情在线观看国产| 老师上课跳d突然被开到最大视频| 亚洲一级一片aⅴ在线观看| 国产色婷婷99| 久久精品综合一区二区三区| 91在线观看av| 日本爱情动作片www.在线观看 | 在线天堂最新版资源| 午夜爱爱视频在线播放| 国产男靠女视频免费网站| 精品日产1卡2卡| 丰满人妻一区二区三区视频av| 国产精品电影一区二区三区| 舔av片在线| 99久久九九国产精品国产免费| 性色avwww在线观看| 国产精品综合久久久久久久免费| 韩国av在线不卡| 亚洲成人中文字幕在线播放| 99久久成人亚洲精品观看| 久久精品影院6| 国产伦一二天堂av在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲av二区三区四区|