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

    基于FFT-SVD 算法的轉(zhuǎn)子軸心軌跡提純與故障檢測(cè)

    2022-02-07 13:43:04肖鑫龍楊洪濤郭曉軍
    惠州學(xué)院學(xué)報(bào) 2022年6期
    關(guān)鍵詞:特征頻率軸心差分

    肖鑫龍,楊洪濤.2*,陳 賀,郭曉軍,李 莉,沈 梅

    (1.安徽理工大學(xué) 機(jī)械工程學(xué)院;2.安徽理工大學(xué) 礦山智能裝備與技術(shù)安徽省重點(diǎn)實(shí)驗(yàn)室,安徽 淮南 232001;3.無(wú)錫市新久揚(yáng)機(jī)械制造有限公司,江蘇 無(wú)錫 214181)

    轉(zhuǎn)子作為重要的旋轉(zhuǎn)機(jī)械零件,工作環(huán)境一般極為惡劣。在轉(zhuǎn)子工作的過(guò)程最容易出現(xiàn)轉(zhuǎn)子不對(duì)中故障[1]。轉(zhuǎn)子不對(duì)中的原因主要是安裝誤差、承載后的變形和環(huán)境溫度的變化[2]。目前,一種方法通過(guò)直接對(duì)轉(zhuǎn)子振動(dòng)信號(hào)的相位、幅值進(jìn)行分析,提取故障特征[3-5],另一種方法是通過(guò)轉(zhuǎn)子軸向和徑向振動(dòng)信號(hào)合成軸心軌跡來(lái)識(shí)別轉(zhuǎn)子故障[6]。軸心軌跡是轉(zhuǎn)子旋轉(zhuǎn)中心軸相對(duì)于軸承座振動(dòng)位移軌跡,它可以實(shí)時(shí)反映出轉(zhuǎn)子的工作狀態(tài),相比于直接分析振動(dòng)信號(hào)具有更強(qiáng)的分辨性,因此被廣泛應(yīng)用于轉(zhuǎn)子的故障診斷[7]。但在采集轉(zhuǎn)子振動(dòng)信號(hào)過(guò)程中往往因?yàn)榘罅康脑肼暥构收闲盘?hào)特征淹沒(méi),很難通過(guò)軸心軌跡來(lái)識(shí)別出轉(zhuǎn)子產(chǎn)生的故障。因此把噪聲和干擾信號(hào)從原始信號(hào)中分離出來(lái)對(duì)轉(zhuǎn)子故障軸心軌跡的識(shí)別至關(guān)重要。

    現(xiàn)在對(duì)軸心軌跡信號(hào)進(jìn)行提純的算法很多,例如主成分分析(Principal Component Analysis,PCA)[8]、小波 閾值(Wavelet Threshold,WT)[9]、變分模態(tài) 分解(Variational Modal Decomposition,VMD)[10]、經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Modal Decomposition,EMD)[11]、數(shù)學(xué)形態(tài)學(xué)法(Mathematical Morphology,MM)[12]等。PCA 算法方差小的非主成分可能含有對(duì)樣本差異的重要信息;WT算法存在各頻帶間能量交疊的現(xiàn)象;VMD算法會(huì)隨著時(shí)間的推移在全局范圍內(nèi)造成信號(hào)重疊;EMD算法會(huì)造成信號(hào)的模態(tài)混疊;MM 算法存在結(jié)構(gòu)元素形狀和大小難以選擇等問(wèn)題[13]。

    奇異值分解(Singular Value Decomposition,SVD)是一種非線性信號(hào)濾波處理方法,非常適合應(yīng)用于非平穩(wěn)信號(hào)和非線性信號(hào)的特征提取。錢征文等[14]通過(guò)仿真驗(yàn)證奇異值分解對(duì)白噪聲和色噪聲有較好的降噪效果;張景潤(rùn)等[15]通過(guò)奇異差分譜理論清晰的對(duì)軸心軌跡進(jìn)行提純;郭明軍等[16]從能量損失角度提出有效奇異值的選取方法,進(jìn)一步完善了差分譜理論并對(duì)軸心軌跡降噪提純成功。但是在研究中發(fā)現(xiàn)當(dāng)奇異值差分譜的第二個(gè)或者第三個(gè)峰值同樣很大的時(shí),只通過(guò)差分譜選取峰值進(jìn)行信號(hào)重構(gòu)會(huì)造成信號(hào)失真,故障特征會(huì)和噪聲一起被過(guò)濾掉。

    因此,本文針對(duì)轉(zhuǎn)子軸心軌跡存在噪聲、干擾等問(wèn)題,研究了一種FFT-SVD 利用特征頻率選取奇異值,對(duì)軸心軌跡進(jìn)行降噪提純的算法。經(jīng)實(shí)驗(yàn)分析,F(xiàn)FTSVD 算法能有效提純軸心軌跡,并成功檢測(cè)出轉(zhuǎn)子不對(duì)中故障。和奇異值差分譜、EMD、小波閾值降噪算法相比較,本文研究算法對(duì)軸心軌跡的降噪提純更清晰光滑有效。

    1 FFT-SVD降噪原理

    圖1 為FFT-SVD 降噪原理流程圖,首先將采集到的轉(zhuǎn)子軸向信號(hào)X和徑向信號(hào)Y轉(zhuǎn)換為Hankel 矩陣,再進(jìn)行奇異值分解,對(duì)原始信號(hào)進(jìn)行FFT 變換得到頻譜圖,根據(jù)特征頻率選取有效奇異值來(lái)重構(gòu)信號(hào),最后將重構(gòu)后的信號(hào)合成軸心軌跡。

    圖1 FFT-SVD降噪流程圖

    1.1 Hankel矩陣的構(gòu)造

    假設(shè)原始信號(hào)為x= [x1,x2, … ,xl],則構(gòu)成的m行n列的Hankel矩陣X為:

    其中l(wèi)=m+n- 1,研究表明Hankel 矩陣的行數(shù)和列數(shù)越相近信號(hào)提純的效果越好,因此,當(dāng)l為偶數(shù)時(shí),m=l/2,n=l/2 + 1;當(dāng)l為奇數(shù)時(shí),m= (l+ 1)/2,n= (l+ 1)/2。

    1.2 SVD分解

    任意實(shí)矩陣X Rm n必存在正交矩陣U Rm m和V Rn n使下式成立

    式(2)稱為X的奇異值分解,其中左奇異矩陣

    為X的奇異值個(gè)數(shù),12r0。

    式(2)可展開為:

    或?qū)⒕仃嘦、V寫成分量組合的形式

    1.3 有效奇異值選取

    上述將原始信號(hào)x構(gòu)造成Hankel 矩陣X,假設(shè)有r個(gè)奇異值,則經(jīng)奇異值分解得到各分量矩陣如式(5)所示:

    本文通過(guò)FFT-SVD 算法利用特征頻率選取有效奇異值,假設(shè)有效奇異值個(gè)數(shù)為k,即

    1.4 矩陣重構(gòu)

    由以上選取前k個(gè)有效分量矩陣進(jìn)行信號(hào)重構(gòu),則重構(gòu)后的矩陣X*為:

    其中左奇異矩陣[U1,U2, ,Uk],右奇異矩陣,奇異值S[diag(1,2, ,k)],則重構(gòu)后的矩陣為:

    2 實(shí)驗(yàn)分析

    實(shí)驗(yàn)流程如圖2 所示,通過(guò)搭建轉(zhuǎn)子故障試驗(yàn)臺(tái)采集轉(zhuǎn)子信號(hào)來(lái)合成軸心軌跡。

    圖2 實(shí)驗(yàn)流程圖

    2.1 實(shí)驗(yàn)裝置

    本文研究對(duì)象為如圖3 所示的轉(zhuǎn)子實(shí)驗(yàn)臺(tái)裝置,主要由800 mm*150 mm的實(shí)驗(yàn)臺(tái)體、400 W伺服電機(jī)、10 mm 直軸鋼、6200 深溝球軸承、可調(diào)偏心軸承座、聯(lián)軸器和負(fù)載盤組成。

    圖3 轉(zhuǎn)子實(shí)驗(yàn)臺(tái)裝置圖

    2.2 軸心軌跡數(shù)據(jù)采集系統(tǒng)

    轉(zhuǎn)子軸心軌跡信號(hào)采集系統(tǒng)主要由HD9200 數(shù)據(jù)采集器、2 個(gè)WT0150 電渦流傳感器、HTX-V3 信號(hào)調(diào)理器、轉(zhuǎn)速控制箱和上位機(jī)組成,如圖4 所示。2 個(gè)相互垂直的WT0150電渦流傳感器安裝位置如圖5所示,WT0150 電渦流傳感器的頻率響應(yīng)為0~10 kHz,靈敏度為8 V/mm,通過(guò)采集轉(zhuǎn)子的軸向位移信號(hào)X和徑向位移信號(hào)Y合成軸心軌跡。

    圖4 HD9200數(shù)據(jù)采集系統(tǒng)

    圖5 電渦流傳感器安裝位置

    2.3 模擬轉(zhuǎn)子不對(duì)中及信號(hào)采集

    本實(shí)驗(yàn)主要模擬轉(zhuǎn)子不對(duì)中故障,通過(guò)電渦流傳感器采集的轉(zhuǎn)子軸向和徑向位移信號(hào),經(jīng)FFT-SVD降噪提純后擬合軸心軌跡。其HD9200數(shù)據(jù)采集系統(tǒng)中信號(hào)的采樣率為10 240 Hz,采樣時(shí)間為0.5 s,采樣點(diǎn)數(shù)為5 120個(gè),電機(jī)轉(zhuǎn)速設(shè)置為3 000 r/min。

    通過(guò)調(diào)節(jié)如圖6所示可調(diào)節(jié)偏心軸承座的不對(duì)中距離,來(lái)模擬轉(zhuǎn)子不對(duì)中,用HD9200數(shù)據(jù)采集設(shè)備采集轉(zhuǎn)子不對(duì)中軸向信號(hào)X和徑向信號(hào)Y,如表1為部分信號(hào)數(shù)據(jù)(前15個(gè))。

    表1 信號(hào)X和Y的部分?jǐn)?shù)據(jù)

    圖6 可調(diào)偏心軸承座

    2.4 合成軸心軌跡

    對(duì)實(shí)驗(yàn)所采集的轉(zhuǎn)子軸向信號(hào)X和徑向信號(hào)Y,首先畫出如圖7所示的信號(hào)X時(shí)域圖和Y的時(shí)域圖,再將信號(hào)X和Y合成如圖8所示的軸心軌跡圖,可見(jiàn)直接通過(guò)X和Y合成的原始軸心軌跡雜亂無(wú)章,無(wú)法通過(guò)軸心軌跡來(lái)檢測(cè)轉(zhuǎn)子故障。

    圖7 信號(hào)X和Y的時(shí)域圖

    圖8 原始軸心軌跡圖

    3 FFT-SVD軸心軌跡降噪

    根據(jù)轉(zhuǎn)子不對(duì)中實(shí)驗(yàn)采集轉(zhuǎn)子的橫向信號(hào)X和徑向信號(hào)Y按照?qǐng)D1流程進(jìn)行降噪,步驟如下:

    (1)對(duì)采集的轉(zhuǎn)子軸向信號(hào)X和徑向信號(hào)Y按照式(1)構(gòu)造Hankel矩陣;

    (2)將矩陣按照式(2)進(jìn)行SVD 降噪,如表2 所示為信號(hào)X和信號(hào)Y的奇異值(取前10個(gè)奇異值);

    (3)對(duì)信號(hào)X和信號(hào)Y進(jìn)行FFT變換得到如圖9所示的頻域圖(取前500 Hz);

    (4)根據(jù)頻域圖中主要頻率個(gè)數(shù)選取奇異值個(gè)數(shù);

    (5)將確定的奇異值個(gè)數(shù)按照式(7)進(jìn)行矩陣重構(gòu),得到重構(gòu)信號(hào)X*和Y*;

    (6)把重構(gòu)后的信號(hào)X*和Y*合成軸心軌跡,完成軸心軌跡的降噪提純;

    圖9 信號(hào)X和Y的頻域圖

    表2 為信號(hào)X和信號(hào)Y前10 個(gè)奇異值的具體數(shù)據(jù),表3為信號(hào)X和信號(hào)Y頻域圖中前5個(gè)幅值對(duì)應(yīng)的頻率數(shù)據(jù)。

    表2 信號(hào)X和Y的奇異值

    表3 信號(hào)X和Y頻率幅值及對(duì)應(yīng)頻率

    隨機(jī)取表2 中信號(hào)X和信號(hào)Y第1 對(duì)奇異值(22.520,22.499)和(18.769,18.764)、第2 對(duì)奇異值(11.388,11.367)和(8.636,8.630)、信號(hào)X的第3 對(duì)奇異值(1.386,1.385)、信號(hào)Y的第5 對(duì)奇異值(1.308,1.307)來(lái)分別重構(gòu)信號(hào),經(jīng)FFT變換后得到如表4所示的每對(duì)奇異值對(duì)應(yīng)的頻率和原始頻率對(duì)比圖。

    實(shí)驗(yàn)研究發(fā)現(xiàn),原始信號(hào)幅值頻率和每對(duì)奇異值按照從大到小排列,每2個(gè)奇異值對(duì)應(yīng)一個(gè)頻率幅值,因此可以用特征頻率個(gè)數(shù)選取有效奇異值來(lái)重構(gòu)信號(hào)。

    表4 奇異值重構(gòu)后頻率和幅值

    由圖9中信號(hào)X和信號(hào)Y頻域圖可知,特征頻率主要有50 Hz 和100 Hz,分別為信號(hào)的基頻和二倍頻,和理論上轉(zhuǎn)子發(fā)生不對(duì)中特征頻率(主要是基頻常伴隨著二倍頻)一致,根據(jù)頻率和奇異值的對(duì)應(yīng)關(guān)系,選取前4個(gè)奇異值按照式(7)來(lái)重構(gòu)信號(hào),則重構(gòu)后的信號(hào)分別為X*和Y*如圖10所示,把信號(hào)X*和Y*合成軸心軌跡,如圖11為FFT-SVD算法降噪提純后的軸心軌跡。

    圖10 信號(hào)X*和Y*的時(shí)域圖

    圖11 FFT-SVD降噪后軸心軌跡

    文獻(xiàn)[13]表明轉(zhuǎn)子系統(tǒng)出現(xiàn)不對(duì)中故障時(shí)理論軸心軌跡為香蕉形,嚴(yán)重為8字形。

    圖11 中降噪提純后的軸心軌跡清晰光滑為香蕉形,與理論相符。因此,F(xiàn)FT-SVD 算法根據(jù)特征頻率選取奇異值,能有效對(duì)轉(zhuǎn)子軸心軌跡進(jìn)行提純,并成功檢測(cè)出轉(zhuǎn)子不對(duì)中故障。

    4 其他算法軸心軌跡降噪

    奇異值差分譜算法是通過(guò)上一個(gè)奇異值和下一個(gè)奇異值做差,按照差分譜峰值對(duì)應(yīng)的序列來(lái)確定有效奇異值,再根據(jù)信號(hào)重構(gòu)實(shí)現(xiàn)信號(hào)降噪提純。如圖12所示為信號(hào)X和Y奇異值和奇異值差分譜,圖中差分譜的峰值均對(duì)應(yīng)序列2,因此選取前兩個(gè)奇異值進(jìn)行信號(hào)重構(gòu),將重構(gòu)后的信號(hào)合成軸心軌跡如圖13所示。

    圖12 信號(hào)X和Y的奇異值和差分譜

    圖13 奇異值差分譜降噪后軸心軌跡

    按照EMD算法和小波閾值算法對(duì)信號(hào)X和信號(hào)Y進(jìn)行降噪處理,將降噪后的信號(hào)合成的軸心軌跡如圖14所示。

    如圖12 信號(hào)X和Y的奇異值差分譜所示,當(dāng)差分譜的第1 個(gè)峰值和第2 個(gè)峰值相差不大時(shí)(信號(hào)X為89.8%,信號(hào)Y為70.6%),通過(guò)第1 個(gè)峰值選取奇異值重構(gòu)信號(hào)會(huì)造成信號(hào)失真,合成的軸心軌跡畸形,無(wú)法通過(guò)軸心軌跡檢測(cè)轉(zhuǎn)子故障。

    實(shí)驗(yàn)結(jié)果表明:和奇異值差分譜、EMD 和小波閾值算法降噪提純相比較,F(xiàn)FT-SVD 算法通過(guò)特征頻率選取奇異值,能有效對(duì)軸心軌跡進(jìn)行降噪提純,并成功檢測(cè)出轉(zhuǎn)子不對(duì)中故障。

    圖14 EMD和小波閾值降噪后的軸心軌跡

    5 結(jié)論

    本文將采集的轉(zhuǎn)子軸向信號(hào)X和徑向信號(hào)Y進(jìn)行FFT 變換和SVD 分解,發(fā)現(xiàn)頻率和奇異值存在對(duì)應(yīng)關(guān)系,因此,根據(jù)特征頻率個(gè)數(shù)確定有效奇異值,進(jìn)行信號(hào)重構(gòu)并合成軸心軌跡完成降噪提純,最后和奇異值差分譜、EMD 和小波閾值降噪算法相比較,得到以下結(jié)論:

    (1)實(shí)驗(yàn)發(fā)現(xiàn),信號(hào)頻域圖中頻率和SVD分解后的奇異值存在對(duì)應(yīng)關(guān)系,按照從大到小的排列順序,每?jī)蓚€(gè)奇異值重構(gòu)后的信號(hào)頻率和原始信號(hào)頻率對(duì)應(yīng)。由信號(hào)X和Y頻譜圖中可見(jiàn)轉(zhuǎn)子主要頻率為2個(gè)(信號(hào)X和信號(hào)Y均為50 Hz,100 Hz),因此,選取前4個(gè)奇異值進(jìn)行重構(gòu),將重構(gòu)后的信號(hào)合成軸心軌跡。重構(gòu)后的軸心軌跡清晰光滑,有效去除了干擾信號(hào),并成功檢測(cè)出轉(zhuǎn)子不對(duì)中故障。

    (2)使用奇異值差分譜算法對(duì)轉(zhuǎn)子軸心軌跡降噪提純發(fā)現(xiàn):當(dāng)差分譜中出現(xiàn)第二個(gè)峰值且和第一個(gè)峰值相差不大(信號(hào)X是89.8%,信號(hào)Y是70.6%)時(shí),按照差分譜使用第一個(gè)峰值選取奇異值進(jìn)行重構(gòu)信號(hào),會(huì)導(dǎo)致信號(hào)失真、丟失重要特征、合成的軸心軌跡畸形,很難去檢測(cè)轉(zhuǎn)子故障。

    (3)和EMD、奇異值差分譜、小波閾值降噪算法相比較,F(xiàn)FT-SVD 算法對(duì)軸心軌跡的提純最干凈有效,小波閾值算法降噪次之,EMD分解算法降噪殘留噪聲最多。因此,本文研究算法能有效能用于轉(zhuǎn)子軸心軌跡提純,并成功檢測(cè)出轉(zhuǎn)子不對(duì)中故障。

    猜你喜歡
    特征頻率軸心差分
    數(shù)列與差分
    瓷磚檢測(cè)機(jī)器人的聲音信號(hào)處理
    鋼結(jié)構(gòu)軸心受壓構(gòu)件穩(wěn)定性分析
    光學(xué)波前參數(shù)的分析評(píng)價(jià)方法研究
    基于振動(dòng)信號(hào)特征頻率的數(shù)控車床故障辨識(shí)方法
    基于小波去噪和EMD算法在齒輪故障檢測(cè)中的應(yīng)用
    CFRP和角鋼復(fù)合加固混凝土矩形柱軸心受壓承載力
    以門靜脈-腸系膜上靜脈為軸心的腹腔鏡胰十二指腸切除術(shù)16例報(bào)道
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    相對(duì)差分單項(xiàng)測(cè)距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    亚洲av在线观看美女高潮| 女人十人毛片免费观看3o分钟| 能在线免费看毛片的网站| 秋霞在线观看毛片| 观看免费一级毛片| av福利片在线观看| 少妇人妻一区二区三区视频| 国产高清国产精品国产三级 | 国产成年人精品一区二区| 亚洲国产精品成人综合色| 少妇人妻精品综合一区二区| 国产欧美日韩精品一区二区| 久久久亚洲精品成人影院| 男女边摸边吃奶| 日日啪夜夜撸| 18禁在线播放成人免费| 欧美3d第一页| 伊人久久精品亚洲午夜| 国产一区二区在线观看日韩| 国产成人a∨麻豆精品| 男人狂女人下面高潮的视频| 美女被艹到高潮喷水动态| 亚洲第一区二区三区不卡| 久久精品国产自在天天线| 18禁在线播放成人免费| 91精品国产九色| 亚洲无线观看免费| 精品人妻视频免费看| 淫秽高清视频在线观看| 亚洲av免费高清在线观看| 成年免费大片在线观看| 久久久久网色| 免费观看a级毛片全部| 欧美精品一区二区大全| 久久久久精品久久久久真实原创| 乱人视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲婷婷狠狠爱综合网| 最近2019中文字幕mv第一页| 老女人水多毛片| 精品人妻一区二区三区麻豆| 国产成人一区二区在线| 亚洲av男天堂| 欧美区成人在线视频| 中文字幕久久专区| 成年人午夜在线观看视频 | 亚洲激情五月婷婷啪啪| 免费无遮挡裸体视频| 中文资源天堂在线| 国产亚洲精品av在线| 国产在线男女| 日韩视频在线欧美| 国产 一区 欧美 日韩| 免费看日本二区| 国产麻豆成人av免费视频| 人妻少妇偷人精品九色| 18禁裸乳无遮挡免费网站照片| 嫩草影院入口| 嫩草影院入口| 精品久久久久久久末码| 国产高清有码在线观看视频| 成人午夜高清在线视频| 欧美激情在线99| 亚洲精品影视一区二区三区av| 女人十人毛片免费观看3o分钟| 日韩成人av中文字幕在线观看| 最近视频中文字幕2019在线8| 亚洲最大成人av| 搡老乐熟女国产| 91久久精品国产一区二区三区| or卡值多少钱| 国产人妻一区二区三区在| 欧美日韩国产mv在线观看视频 | 伦精品一区二区三区| 特级一级黄色大片| 国内精品一区二区在线观看| 亚洲,欧美,日韩| 久久久久久久久久成人| 中文字幕av在线有码专区| 亚洲国产精品国产精品| 午夜老司机福利剧场| 美女国产视频在线观看| 亚洲最大成人中文| 午夜福利在线观看吧| 国产免费一级a男人的天堂| 久久久久久久大尺度免费视频| 久久精品夜夜夜夜夜久久蜜豆| 秋霞伦理黄片| 乱人视频在线观看| 高清午夜精品一区二区三区| 欧美xxxx性猛交bbbb| 免费播放大片免费观看视频在线观看| 你懂的网址亚洲精品在线观看| 国产成人精品福利久久| 亚洲国产日韩欧美精品在线观看| 97超碰精品成人国产| 99久久人妻综合| 天堂俺去俺来也www色官网 | 男人狂女人下面高潮的视频| 国产精品美女特级片免费视频播放器| 日本一二三区视频观看| 美女cb高潮喷水在线观看| 天美传媒精品一区二区| 高清av免费在线| 国产黄a三级三级三级人| 男女国产视频网站| 色哟哟·www| 欧美日韩视频高清一区二区三区二| 天美传媒精品一区二区| av在线天堂中文字幕| 国产av不卡久久| 床上黄色一级片| 久久久久精品久久久久真实原创| av播播在线观看一区| 男女边摸边吃奶| 日韩av不卡免费在线播放| 永久网站在线| 国产女主播在线喷水免费视频网站 | 久久精品久久精品一区二区三区| 夫妻午夜视频| 中文精品一卡2卡3卡4更新| 久久久久久伊人网av| 亚洲美女视频黄频| 热99在线观看视频| 看非洲黑人一级黄片| 亚洲av国产av综合av卡| 国产av在哪里看| 一级毛片黄色毛片免费观看视频| av.在线天堂| 黄色日韩在线| 国产高清不卡午夜福利| 亚洲第一区二区三区不卡| 全区人妻精品视频| 免费大片黄手机在线观看| 天堂av国产一区二区熟女人妻| 久久精品国产鲁丝片午夜精品| 国产美女午夜福利| 成年人午夜在线观看视频 | 久久亚洲国产成人精品v| 午夜激情福利司机影院| 国产精品熟女久久久久浪| 嫩草影院新地址| 欧美性猛交╳xxx乱大交人| 天堂√8在线中文| 精品一区二区免费观看| 亚洲自拍偷在线| 三级国产精品欧美在线观看| 肉色欧美久久久久久久蜜桃 | 真实男女啪啪啪动态图| 美女高潮的动态| 麻豆乱淫一区二区| 三级国产精品片| 三级毛片av免费| 最近2019中文字幕mv第一页| 久久精品夜色国产| 能在线免费观看的黄片| 国产视频内射| 超碰97精品在线观看| 欧美丝袜亚洲另类| 免费观看无遮挡的男女| 国产黄色小视频在线观看| or卡值多少钱| 综合色丁香网| 校园人妻丝袜中文字幕| 国产精品精品国产色婷婷| 欧美精品国产亚洲| 亚洲一级一片aⅴ在线观看| 亚洲av成人精品一二三区| 亚洲精品亚洲一区二区| 亚洲精品,欧美精品| 日本wwww免费看| 国产 亚洲一区二区三区 | 亚洲精品日韩在线中文字幕| 亚洲av国产av综合av卡| 国产精品无大码| 日本与韩国留学比较| 国产一区二区在线观看日韩| 日本与韩国留学比较| 3wmmmm亚洲av在线观看| 久久99热这里只有精品18| 欧美一区二区亚洲| 成人av在线播放网站| 一级毛片aaaaaa免费看小| 精品熟女少妇av免费看| 亚洲成人中文字幕在线播放| 久久久久久国产a免费观看| 男插女下体视频免费在线播放| 在线观看人妻少妇| 国内精品美女久久久久久| 亚洲精品乱久久久久久| 好男人视频免费观看在线| 亚洲一区高清亚洲精品| 国产伦一二天堂av在线观看| 午夜日本视频在线| 久久久久久国产a免费观看| 美女cb高潮喷水在线观看| 国产在视频线在精品| 80岁老熟妇乱子伦牲交| 国产成人精品福利久久| 日韩精品青青久久久久久| 久久精品综合一区二区三区| 亚洲欧美清纯卡通| 99久久人妻综合| 人人妻人人澡人人爽人人夜夜 | 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| 久久精品熟女亚洲av麻豆精品 | 国产亚洲最大av| 99热全是精品| 99热这里只有是精品50| 天堂中文最新版在线下载 | 免费观看a级毛片全部| 亚洲aⅴ乱码一区二区在线播放| 嘟嘟电影网在线观看| 午夜免费男女啪啪视频观看| 国产乱人视频| 精品一区在线观看国产| 能在线免费观看的黄片| 日韩电影二区| 午夜福利网站1000一区二区三区| 亚洲内射少妇av| 亚洲精品日本国产第一区| 国产精品麻豆人妻色哟哟久久 | 偷拍熟女少妇极品色| 少妇熟女欧美另类| 亚洲av中文字字幕乱码综合| 欧美三级亚洲精品| 国产伦精品一区二区三区视频9| 精品久久久精品久久久| 蜜桃久久精品国产亚洲av| av福利片在线观看| 高清日韩中文字幕在线| 一本久久精品| 亚洲国产精品国产精品| 欧美日韩一区二区视频在线观看视频在线 | 高清毛片免费看| av福利片在线观看| 久久久国产一区二区| 精品人妻熟女av久视频| 久久久欧美国产精品| 国产精品福利在线免费观看| 青青草视频在线视频观看| 免费在线观看成人毛片| 亚洲精华国产精华液的使用体验| 久久99蜜桃精品久久| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 在线a可以看的网站| 国产亚洲91精品色在线| 成年人午夜在线观看视频 | 亚洲真实伦在线观看| 日本爱情动作片www.在线观看| 天堂俺去俺来也www色官网 | 婷婷色综合大香蕉| 一个人观看的视频www高清免费观看| 精品亚洲乱码少妇综合久久| 免费av不卡在线播放| 国产成年人精品一区二区| 亚洲欧美中文字幕日韩二区| 日韩视频在线欧美| 男女啪啪激烈高潮av片| 成人漫画全彩无遮挡| 中文字幕亚洲精品专区| 91狼人影院| 久久国产乱子免费精品| 在线观看免费高清a一片| 女的被弄到高潮叫床怎么办| 大陆偷拍与自拍| 国产欧美另类精品又又久久亚洲欧美| 久久久久网色| 国产成人午夜福利电影在线观看| 九九在线视频观看精品| 亚洲怡红院男人天堂| 午夜老司机福利剧场| 免费少妇av软件| 国产女主播在线喷水免费视频网站 | 在线免费十八禁| 中文资源天堂在线| 亚洲成人中文字幕在线播放| 能在线免费看毛片的网站| 小蜜桃在线观看免费完整版高清| 国产精品无大码| 18禁裸乳无遮挡免费网站照片| 国产免费又黄又爽又色| 日韩av在线免费看完整版不卡| 国产亚洲91精品色在线| 麻豆乱淫一区二区| 久久久久久久久久成人| 日本wwww免费看| 亚洲国产色片| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆| 久久6这里有精品| 国产精品麻豆人妻色哟哟久久 | 美女高潮的动态| 久久久久久久久久人人人人人人| 美女黄网站色视频| 97精品久久久久久久久久精品| 国产成人freesex在线| 亚州av有码| 免费在线观看成人毛片| 亚洲av电影不卡..在线观看| 欧美一级a爱片免费观看看| av天堂中文字幕网| 午夜激情欧美在线| 免费少妇av软件| 色网站视频免费| 免费看a级黄色片| 超碰av人人做人人爽久久| 成年免费大片在线观看| 大香蕉97超碰在线| 又爽又黄a免费视频| 成人毛片a级毛片在线播放| 97在线视频观看| 久久精品夜色国产| 国产一区有黄有色的免费视频 | 边亲边吃奶的免费视频| 国产成人精品久久久久久| 亚洲av福利一区| 午夜福利视频1000在线观看| 欧美变态另类bdsm刘玥| 成人亚洲欧美一区二区av| 日本色播在线视频| 日韩av在线大香蕉| 亚洲久久久久久中文字幕| 国产免费福利视频在线观看| 国产爱豆传媒在线观看| 色综合站精品国产| 少妇熟女欧美另类| 午夜免费男女啪啪视频观看| 大香蕉久久网| 日韩不卡一区二区三区视频在线| 午夜老司机福利剧场| 久久久欧美国产精品| 两个人的视频大全免费| 99久久中文字幕三级久久日本| 午夜福利视频1000在线观看| 国产免费又黄又爽又色| 亚洲精品亚洲一区二区| 在线观看免费高清a一片| 国产成人精品一,二区| 亚洲电影在线观看av| 精品久久久久久久久av| 国产单亲对白刺激| 免费观看av网站的网址| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影 | 欧美高清成人免费视频www| 久久精品久久久久久噜噜老黄| 高清欧美精品videossex| 欧美变态另类bdsm刘玥| 黄片无遮挡物在线观看| 免费大片18禁| videos熟女内射| 插逼视频在线观看| 边亲边吃奶的免费视频| 国产精品国产三级国产av玫瑰| 韩国av在线不卡| 欧美激情国产日韩精品一区| 国产亚洲91精品色在线| 国产v大片淫在线免费观看| 人妻夜夜爽99麻豆av| 久久6这里有精品| 亚洲av电影不卡..在线观看| 男的添女的下面高潮视频| 一级片'在线观看视频| 亚洲一区高清亚洲精品| 国产高清不卡午夜福利| 成人特级av手机在线观看| 日本三级黄在线观看| 美女脱内裤让男人舔精品视频| 久久精品久久久久久噜噜老黄| 成人午夜高清在线视频| 少妇熟女欧美另类| 观看免费一级毛片| 国产亚洲5aaaaa淫片| 黄色日韩在线| 波多野结衣巨乳人妻| 看非洲黑人一级黄片| 欧美另类一区| 亚洲av在线观看美女高潮| 久久精品人妻少妇| av国产久精品久网站免费入址| 看黄色毛片网站| 婷婷色综合大香蕉| 国产精品.久久久| 国产黄片美女视频| 久久人人爽人人爽人人片va| 成人午夜高清在线视频| 波多野结衣巨乳人妻| 国产 亚洲一区二区三区 | 欧美zozozo另类| 80岁老熟妇乱子伦牲交| 99久久精品热视频| 国产亚洲一区二区精品| 精品久久久久久成人av| 亚洲18禁久久av| 亚洲欧美一区二区三区国产| 精品酒店卫生间| 噜噜噜噜噜久久久久久91| 免费看日本二区| 大话2 男鬼变身卡| 日日摸夜夜添夜夜爱| 日韩欧美一区视频在线观看 | 搡老乐熟女国产| 十八禁国产超污无遮挡网站| 蜜臀久久99精品久久宅男| 亚洲欧美中文字幕日韩二区| 成人亚洲精品一区在线观看 | 久久久久久久久中文| 久久久成人免费电影| 三级毛片av免费| 老司机影院毛片| 亚洲伊人久久精品综合| 精品国产三级普通话版| 亚洲国产欧美在线一区| 久久久久久国产a免费观看| 成人特级av手机在线观看| 五月伊人婷婷丁香| 国产 一区精品| 久久99精品国语久久久| 嫩草影院入口| 高清视频免费观看一区二区 | 久久精品国产亚洲av天美| 日日啪夜夜撸| 免费观看无遮挡的男女| 18+在线观看网站| av又黄又爽大尺度在线免费看| 欧美极品一区二区三区四区| 神马国产精品三级电影在线观看| 久久国内精品自在自线图片| 久久精品国产亚洲av涩爱| 国产精品一区二区三区四区久久| 精品99又大又爽又粗少妇毛片| 日韩一本色道免费dvd| 97热精品久久久久久| 国产伦精品一区二区三区四那| 国产一区亚洲一区在线观看| 久久精品国产亚洲av天美| 亚洲精品日韩av片在线观看| 国产精品嫩草影院av在线观看| 成人亚洲精品一区在线观看 | 日本一本二区三区精品| 五月伊人婷婷丁香| 建设人人有责人人尽责人人享有的 | 日韩av在线大香蕉| 亚洲,欧美,日韩| 精品久久国产蜜桃| 国产 一区 欧美 日韩| 三级国产精品欧美在线观看| 插阴视频在线观看视频| 久久久欧美国产精品| 亚洲av.av天堂| 丝瓜视频免费看黄片| 亚洲成人中文字幕在线播放| 亚洲精品成人av观看孕妇| 1000部很黄的大片| 亚洲精品国产av成人精品| 女人久久www免费人成看片| 亚洲av一区综合| 亚洲人成网站高清观看| 精品人妻视频免费看| 麻豆乱淫一区二区| 日韩精品青青久久久久久| 国产老妇伦熟女老妇高清| 网址你懂的国产日韩在线| 成人一区二区视频在线观看| 国产白丝娇喘喷水9色精品| 欧美性猛交╳xxx乱大交人| 日韩一区二区视频免费看| 国产成年人精品一区二区| 高清在线视频一区二区三区| eeuss影院久久| 日本黄色片子视频| 2022亚洲国产成人精品| 2021少妇久久久久久久久久久| 波野结衣二区三区在线| 国产av国产精品国产| 麻豆成人av视频| 舔av片在线| 午夜免费男女啪啪视频观看| 国产黄a三级三级三级人| 日韩电影二区| 精品一区二区免费观看| 亚洲精品国产av蜜桃| 五月伊人婷婷丁香| 国产精品伦人一区二区| 久久久午夜欧美精品| 日韩 亚洲 欧美在线| 久久久久久九九精品二区国产| 亚洲精品国产av成人精品| 欧美区成人在线视频| 好男人视频免费观看在线| 成年版毛片免费区| 91在线精品国自产拍蜜月| 久久久久久久久久成人| 嫩草影院精品99| 亚洲人成网站在线观看播放| 麻豆国产97在线/欧美| 成年女人在线观看亚洲视频 | 一级片'在线观看视频| 久久久久久久亚洲中文字幕| 99热网站在线观看| 蜜桃亚洲精品一区二区三区| 亚洲av中文av极速乱| 久久精品夜色国产| 国产成人精品久久久久久| 乱人视频在线观看| av专区在线播放| 中文精品一卡2卡3卡4更新| 中文在线观看免费www的网站| 99久久精品热视频| 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 晚上一个人看的免费电影| 最新中文字幕久久久久| 国产精品久久久久久久久免| 国产av在哪里看| 天堂俺去俺来也www色官网 | 最近中文字幕高清免费大全6| 欧美激情国产日韩精品一区| 一级毛片久久久久久久久女| 欧美人与善性xxx| 日本av手机在线免费观看| 日本熟妇午夜| 欧美性感艳星| 国产伦精品一区二区三区四那| 国国产精品蜜臀av免费| 中文资源天堂在线| 欧美性感艳星| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 精品酒店卫生间| 久久精品久久精品一区二区三区| 国产毛片a区久久久久| 久久97久久精品| 免费无遮挡裸体视频| videossex国产| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 三级国产精品欧美在线观看| 久久久久精品性色| 91午夜精品亚洲一区二区三区| 国产视频内射| 亚洲国产精品sss在线观看| 韩国高清视频一区二区三区| 神马国产精品三级电影在线观看| 建设人人有责人人尽责人人享有的 | 日本wwww免费看| 高清在线视频一区二区三区| 在线免费观看的www视频| 联通29元200g的流量卡| 国国产精品蜜臀av免费| 国产精品人妻久久久久久| 青青草视频在线视频观看| 伊人久久精品亚洲午夜| av在线天堂中文字幕| 日韩,欧美,国产一区二区三区| 久久6这里有精品| 色网站视频免费| 亚洲av免费高清在线观看| av在线亚洲专区| 国产av国产精品国产| 欧美日韩综合久久久久久| 亚洲国产最新在线播放| 在线天堂最新版资源| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 能在线免费看毛片的网站| 人妻夜夜爽99麻豆av| 国产中年淑女户外野战色| 亚洲熟女精品中文字幕| 干丝袜人妻中文字幕| 久久精品国产自在天天线| 如何舔出高潮| 亚洲欧美成人综合另类久久久| 国产大屁股一区二区在线视频| 日韩欧美三级三区| 亚洲国产精品国产精品| 亚洲伊人久久精品综合| 亚洲精品国产av蜜桃| 久久99蜜桃精品久久| 久久精品夜色国产| 身体一侧抽搐| 午夜老司机福利剧场| 黄色一级大片看看| 青青草视频在线视频观看| 99久久精品一区二区三区| 国产精品国产三级国产av玫瑰| 特大巨黑吊av在线直播| 日韩精品青青久久久久久| 午夜爱爱视频在线播放| 精品熟女少妇av免费看| 亚洲高清免费不卡视频| 不卡视频在线观看欧美| 噜噜噜噜噜久久久久久91| 欧美一级a爱片免费观看看| 久久久久久久久大av| 成年av动漫网址| 欧美日韩精品成人综合77777| 日韩一区二区三区影片| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区 | 欧美bdsm另类| 欧美3d第一页| 婷婷色麻豆天堂久久| av女优亚洲男人天堂| 三级男女做爰猛烈吃奶摸视频| 69人妻影院| 亚洲国产日韩欧美精品在线观看| 日本猛色少妇xxxxx猛交久久| 十八禁国产超污无遮挡网站| 搡女人真爽免费视频火全软件| 特大巨黑吊av在线直播| 日本与韩国留学比较|