• <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
    亚洲成人中文字幕在线播放| aaaaa片日本免费| 我的女老师完整版在线观看| 国产老妇女一区| 91麻豆av在线| 日韩中文字幕欧美一区二区| 波多野结衣巨乳人妻| 国产精品人妻久久久久久| 日日摸夜夜添夜夜添av毛片 | 亚洲精品456在线播放app | 日本撒尿小便嘘嘘汇集6| 免费在线观看亚洲国产| 美女高潮的动态| 精品日产1卡2卡| 人妻久久中文字幕网| 亚洲天堂国产精品一区在线| 免费在线观看日本一区| 欧美绝顶高潮抽搐喷水| 国产高清三级在线| 一个人看的www免费观看视频| 少妇人妻一区二区三区视频| 欧美激情在线99| 国产亚洲欧美98| 午夜免费男女啪啪视频观看 | 国产精品野战在线观看| 亚洲成人久久性| 久久99热这里只有精品18| 观看美女的网站| 99热6这里只有精品| 中文字幕久久专区| 91午夜精品亚洲一区二区三区 | 久9热在线精品视频| 国产色爽女视频免费观看| .国产精品久久| 在线观看一区二区三区| 一区二区三区免费毛片| 国产色爽女视频免费观看| 国产精品亚洲一级av第二区| 亚洲av.av天堂| 久久午夜福利片| 欧美高清性xxxxhd video| 久久久久国产精品人妻aⅴ院| 91av网一区二区| 国产精品精品国产色婷婷| 亚洲人与动物交配视频| 亚洲中文日韩欧美视频| 深夜a级毛片| 我要搜黄色片| 精品久久久久久成人av| 99久久九九国产精品国产免费| 国产真实伦视频高清在线观看 | 国产单亲对白刺激| 国产探花极品一区二区| 丰满的人妻完整版| 国产淫片久久久久久久久 | 国产欧美日韩一区二区三| 女人十人毛片免费观看3o分钟| 成人鲁丝片一二三区免费| 色综合婷婷激情| 国产免费一级a男人的天堂| 国产精品一及| 精品熟女少妇八av免费久了| 国产午夜精品论理片| 最好的美女福利视频网| 欧美一区二区亚洲| 日本在线视频免费播放| 不卡一级毛片| 午夜福利免费观看在线| 男人狂女人下面高潮的视频| 一个人看的www免费观看视频| 国产伦精品一区二区三区视频9| av视频在线观看入口| 极品教师在线视频| 伦理电影大哥的女人| 色噜噜av男人的天堂激情| 最近中文字幕高清免费大全6 | 99久久精品国产亚洲精品| 亚洲午夜理论影院| 精品久久久久久,| 制服丝袜大香蕉在线| 男人狂女人下面高潮的视频| 看十八女毛片水多多多| 欧美精品国产亚洲| 18禁黄网站禁片午夜丰满| 天堂av国产一区二区熟女人妻| 亚洲专区中文字幕在线| 岛国在线免费视频观看| 中文字幕人妻熟人妻熟丝袜美| 欧美区成人在线视频| 国产精品,欧美在线| 18美女黄网站色大片免费观看| 99久久成人亚洲精品观看| 日本成人三级电影网站| 欧美日韩中文字幕国产精品一区二区三区| 无人区码免费观看不卡| 日本熟妇午夜| 高清日韩中文字幕在线| 精品无人区乱码1区二区| 久久久久久久久中文| 免费看a级黄色片| 亚洲av一区综合| 99热这里只有是精品在线观看 | 国产亚洲av嫩草精品影院| 脱女人内裤的视频| 亚洲真实伦在线观看| 一级毛片久久久久久久久女| 成人鲁丝片一二三区免费| 久久久久国产精品人妻aⅴ院| 九九久久精品国产亚洲av麻豆| 在线观看午夜福利视频| 欧美潮喷喷水| 男人舔女人下体高潮全视频| 毛片女人毛片| 国产亚洲精品av在线| 国产午夜精品久久久久久一区二区三区 | 无人区码免费观看不卡| 亚洲欧美日韩东京热| 国产一区二区在线观看日韩| 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 日韩欧美精品v在线| 岛国在线免费视频观看| 日韩欧美国产在线观看| 人人妻,人人澡人人爽秒播| 欧美日韩乱码在线| 久久99热6这里只有精品| 日日摸夜夜添夜夜添av毛片 | 日本免费a在线| 精品日产1卡2卡| 国产视频一区二区在线看| 国产精品久久久久久精品电影| 久久久色成人| 熟妇人妻久久中文字幕3abv| 国产精品免费一区二区三区在线| 美女高潮的动态| 白带黄色成豆腐渣| 精品久久久久久久久av| 免费看日本二区| 伦理电影大哥的女人| 美女黄网站色视频| 99热这里只有是精品在线观看 | 成人av在线播放网站| 99久久99久久久精品蜜桃| 亚洲成人中文字幕在线播放| 中亚洲国语对白在线视频| 久久精品国产亚洲av天美| 亚洲一区二区三区色噜噜| 成人特级av手机在线观看| 国产精品伦人一区二区| 麻豆成人av在线观看| 久久精品影院6| 在线看三级毛片| 白带黄色成豆腐渣| 久久精品综合一区二区三区| 日本三级黄在线观看| 精品国产亚洲在线| 乱人视频在线观看| 99久久成人亚洲精品观看| 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 亚洲男人的天堂狠狠| 亚洲 欧美 日韩 在线 免费| 中文字幕熟女人妻在线| 国产 一区 欧美 日韩| 日韩人妻高清精品专区| 毛片一级片免费看久久久久 | 国产精品电影一区二区三区| 欧美bdsm另类| 美女 人体艺术 gogo| 欧美绝顶高潮抽搐喷水| 日本精品一区二区三区蜜桃| 免费在线观看亚洲国产| 成人特级黄色片久久久久久久| 色哟哟·www| 麻豆av噜噜一区二区三区| 亚洲色图av天堂| 亚洲av成人av| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久亚洲av鲁大| 老鸭窝网址在线观看| 亚洲人成网站高清观看| 在线观看一区二区三区| 免费在线观看亚洲国产| 不卡一级毛片| 毛片一级片免费看久久久久 | 美女被艹到高潮喷水动态| 亚洲专区中文字幕在线| 午夜亚洲福利在线播放| 中文字幕久久专区| 国产真实伦视频高清在线观看 | 精品人妻一区二区三区麻豆 | 国产在线男女| 少妇被粗大猛烈的视频| 免费在线观看成人毛片| 精品福利观看| 欧美性猛交黑人性爽| 亚洲男人的天堂狠狠| 色哟哟哟哟哟哟| 一级毛片久久久久久久久女| 看免费av毛片| 亚洲美女视频黄频| 99久久精品国产亚洲精品| 国产精品98久久久久久宅男小说| 精品久久国产蜜桃| 蜜桃久久精品国产亚洲av| 亚洲第一区二区三区不卡| 亚洲性夜色夜夜综合| 国内久久婷婷六月综合欲色啪| 成人av一区二区三区在线看| 精品免费久久久久久久清纯| 九色国产91popny在线| 99国产极品粉嫩在线观看| 99久久无色码亚洲精品果冻| 亚洲午夜理论影院| 国产欧美日韩精品亚洲av| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 一个人观看的视频www高清免费观看| 啦啦啦观看免费观看视频高清| 久久国产乱子免费精品| 午夜精品久久久久久毛片777| 日韩国内少妇激情av| 亚洲自偷自拍三级| 欧美日韩亚洲国产一区二区在线观看| 国产av在哪里看| 久久天躁狠狠躁夜夜2o2o| 国产精品野战在线观看| 美女高潮喷水抽搐中文字幕| 精品99又大又爽又粗少妇毛片 | 噜噜噜噜噜久久久久久91| 精品久久久久久成人av| 国产午夜精品久久久久久一区二区三区 | 两个人的视频大全免费| 毛片一级片免费看久久久久 | 久久久久国产精品人妻aⅴ院| 十八禁网站免费在线| 久久精品国产99精品国产亚洲性色| 一区二区三区激情视频| 亚洲黑人精品在线| 在线观看舔阴道视频| 亚洲av美国av| 日韩欧美在线二视频| 日韩人妻高清精品专区| 禁无遮挡网站| 看十八女毛片水多多多| 嫩草影院新地址| 亚洲三级黄色毛片| 美女黄网站色视频| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| 亚洲成av人片免费观看| 国产乱人伦免费视频| 中文字幕熟女人妻在线| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 在线看三级毛片| 亚洲 欧美 日韩 在线 免费| 一区二区三区免费毛片| 天天一区二区日本电影三级| a级毛片免费高清观看在线播放| 欧美国产日韩亚洲一区| 少妇丰满av| 可以在线观看的亚洲视频| 波野结衣二区三区在线| 狠狠狠狠99中文字幕| 成人欧美大片| 亚洲自偷自拍三级| 欧美高清性xxxxhd video| 亚洲精品在线观看二区| 亚洲成人久久爱视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av电影在线进入| 亚洲在线观看片| 男人舔女人下体高潮全视频| 亚洲欧美精品综合久久99| 欧美日韩中文字幕国产精品一区二区三区| 日韩大尺度精品在线看网址| www.色视频.com| 女人被狂操c到高潮| 我的老师免费观看完整版| 91久久精品国产一区二区成人| 真人做人爱边吃奶动态| 天堂√8在线中文| 内地一区二区视频在线| 成人毛片a级毛片在线播放| 嫩草影院入口| 国语自产精品视频在线第100页| 综合色av麻豆| 99久久精品一区二区三区| 亚洲美女视频黄频| 亚洲av不卡在线观看| 神马国产精品三级电影在线观看| 日本一二三区视频观看| 日本成人三级电影网站| 18美女黄网站色大片免费观看| 国产成+人综合+亚洲专区| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av天美| 国产野战对白在线观看| 一级黄片播放器| 不卡一级毛片| 欧美成狂野欧美在线观看| 国产白丝娇喘喷水9色精品| 天美传媒精品一区二区| 国产av麻豆久久久久久久| 国产高清激情床上av| 久久精品国产99精品国产亚洲性色| 久久久久久九九精品二区国产| 久久午夜福利片| 亚洲国产精品合色在线| 免费看美女性在线毛片视频| 欧美日韩国产亚洲二区| 亚洲狠狠婷婷综合久久图片| 丁香六月欧美| 床上黄色一级片| 国产 一区 欧美 日韩| 精品国产三级普通话版| 极品教师在线视频| 精品人妻视频免费看| 99久久精品国产亚洲精品| 91久久精品电影网| 国产人妻一区二区三区在| 午夜两性在线视频| 一本久久中文字幕| 黄色日韩在线| 日韩av在线大香蕉| 欧美不卡视频在线免费观看| 亚洲av熟女| 国产精品99久久久久久久久| 国产精品女同一区二区软件 | 嫩草影视91久久| 1024手机看黄色片| 麻豆久久精品国产亚洲av| 日韩欧美精品免费久久 | 亚洲最大成人手机在线| 精品国产亚洲在线| 一本一本综合久久| 身体一侧抽搐| 成人性生交大片免费视频hd| 两个人的视频大全免费| 麻豆国产av国片精品| 亚洲成人免费电影在线观看| 久久久久九九精品影院| 亚洲中文字幕日韩| www.www免费av| 亚洲av二区三区四区| 少妇裸体淫交视频免费看高清| 精品久久久久久久人妻蜜臀av| 变态另类丝袜制服| 精品国内亚洲2022精品成人| 国产一区二区在线av高清观看| 又黄又爽又刺激的免费视频.| 久久这里只有精品中国| 欧美成人免费av一区二区三区| 午夜福利18| 美女高潮的动态| 神马国产精品三级电影在线观看| 亚洲成a人片在线一区二区| 日韩欧美 国产精品| 午夜精品在线福利| ponron亚洲| 精品一区二区三区视频在线| 日韩免费av在线播放| 1024手机看黄色片| 91久久精品电影网| 在线免费观看不下载黄p国产 | 亚洲最大成人av| 免费一级毛片在线播放高清视频| 露出奶头的视频| av专区在线播放| 亚洲熟妇熟女久久| 国产欧美日韩精品一区二区| 露出奶头的视频| 波多野结衣巨乳人妻| 午夜福利免费观看在线| 国产精品不卡视频一区二区 | 国产黄片美女视频| 亚洲精华国产精华精| 国产乱人伦免费视频| 69av精品久久久久久| 免费无遮挡裸体视频| 欧美成人性av电影在线观看| 黄色日韩在线| 高清毛片免费观看视频网站| 精品一区二区三区人妻视频| 1024手机看黄色片| 亚洲精品亚洲一区二区| 精品无人区乱码1区二区| 日本五十路高清| 热99re8久久精品国产| 神马国产精品三级电影在线观看| 99热这里只有是精品在线观看 | 欧美国产日韩亚洲一区| 国产私拍福利视频在线观看| 成人美女网站在线观看视频| 精品久久久久久久久av| 91在线观看av| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 免费搜索国产男女视频| 国产老妇女一区| 亚洲五月天丁香| 女同久久另类99精品国产91| 日韩亚洲欧美综合| 欧美+亚洲+日韩+国产| 国产日本99.免费观看| 无遮挡黄片免费观看| 特大巨黑吊av在线直播| 日本 av在线| 色尼玛亚洲综合影院| 宅男免费午夜| 亚洲三级黄色毛片| 我要看日韩黄色一级片| 特级一级黄色大片| 中亚洲国语对白在线视频| 小说图片视频综合网站| 在现免费观看毛片| 午夜福利高清视频| 夜夜夜夜夜久久久久| 少妇人妻一区二区三区视频| 日本在线视频免费播放| 88av欧美| 18+在线观看网站| 18禁在线播放成人免费| 亚洲av五月六月丁香网| 18禁裸乳无遮挡免费网站照片| 精品国产亚洲在线| 亚洲成人久久爱视频| av中文乱码字幕在线| 性欧美人与动物交配| 午夜福利欧美成人| 精品久久久久久久久久久久久| 深夜精品福利| 国产在线男女| 最近中文字幕高清免费大全6 | 免费在线观看影片大全网站| 一级黄片播放器| 美女大奶头视频| 欧美成人免费av一区二区三区| 九九久久精品国产亚洲av麻豆| 天堂影院成人在线观看| 免费高清视频大片| 搡老熟女国产l中国老女人| 一夜夜www| 欧美最新免费一区二区三区 | 亚洲精品一卡2卡三卡4卡5卡| 最后的刺客免费高清国语| 国产一级毛片七仙女欲春2| 别揉我奶头~嗯~啊~动态视频| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久久久久久久| 亚洲熟妇中文字幕五十中出| 国产亚洲精品av在线| 成人三级黄色视频| 欧美3d第一页| 亚洲综合色惰| 免费一级毛片在线播放高清视频| 亚洲国产欧洲综合997久久,| 亚洲无线观看免费| 99久久无色码亚洲精品果冻| 国产精品一区二区性色av| 国内精品久久久久久久电影| 欧美在线黄色| 淫妇啪啪啪对白视频| 欧美一区二区国产精品久久精品| 欧美精品啪啪一区二区三区| 桃红色精品国产亚洲av| 久久精品人妻少妇| 欧美极品一区二区三区四区| 99热只有精品国产| 看片在线看免费视频| 精品免费久久久久久久清纯| 夜夜看夜夜爽夜夜摸| 熟女人妻精品中文字幕| 一级a爱片免费观看的视频| 色尼玛亚洲综合影院| 99在线人妻在线中文字幕| 亚洲国产日韩欧美精品在线观看| 久久6这里有精品| 亚洲第一区二区三区不卡| 亚洲人与动物交配视频| 国产不卡一卡二| av在线老鸭窝| 久久久精品欧美日韩精品| 亚洲av五月六月丁香网| 久久久国产成人免费| 人人妻,人人澡人人爽秒播| 亚洲avbb在线观看| 国产精品一区二区三区四区久久| 国产成+人综合+亚洲专区| 成人无遮挡网站| 日本黄大片高清| xxxwww97欧美| 99久久精品热视频| 啦啦啦韩国在线观看视频| 一个人免费在线观看电影| av天堂中文字幕网| 国产成人av教育| 五月伊人婷婷丁香| 国产精品一区二区免费欧美| 听说在线观看完整版免费高清| 少妇高潮的动态图| 别揉我奶头~嗯~啊~动态视频| 亚洲乱码一区二区免费版| 直男gayav资源| 亚洲av成人av| 少妇熟女aⅴ在线视频| 午夜免费激情av| 亚洲真实伦在线观看| 国产av一区在线观看免费| 久久人人爽人人爽人人片va | 麻豆国产97在线/欧美| 欧美激情国产日韩精品一区| 性色avwww在线观看| 国产精品久久久久久亚洲av鲁大| 国产精品久久电影中文字幕| 国产一区二区亚洲精品在线观看| 99国产精品一区二区三区| 给我免费播放毛片高清在线观看| 亚洲真实伦在线观看| 精品人妻一区二区三区麻豆 | 久久人人精品亚洲av| 级片在线观看| 噜噜噜噜噜久久久久久91| .国产精品久久| 三级毛片av免费| 最新中文字幕久久久久| 精品国内亚洲2022精品成人| 51午夜福利影视在线观看| 精品一区二区三区视频在线观看免费| 亚洲欧美激情综合另类| 香蕉av资源在线| 国产亚洲精品av在线| 欧美xxxx黑人xx丫x性爽| 此物有八面人人有两片| 别揉我奶头 嗯啊视频| 久久久久久国产a免费观看| 又黄又爽又免费观看的视频| 亚洲国产色片| 中文字幕人成人乱码亚洲影| 国产探花在线观看一区二区| 久久久久久久精品吃奶| 一区二区三区高清视频在线| 校园春色视频在线观看| 免费高清视频大片| 国产毛片a区久久久久| 麻豆成人午夜福利视频| 欧美在线黄色| 亚洲av中文字字幕乱码综合| 欧美黑人巨大hd| 日韩欧美国产一区二区入口| 国产久久久一区二区三区| 国产精品,欧美在线| 久久久精品大字幕| av在线天堂中文字幕| 亚洲人成网站在线播放欧美日韩| 久99久视频精品免费| 国产三级在线视频| 91麻豆av在线| 国产成人福利小说| 国产精品98久久久久久宅男小说| 免费在线观看亚洲国产| 黄片小视频在线播放| 伊人久久精品亚洲午夜| 99在线视频只有这里精品首页| 色噜噜av男人的天堂激情| 成人欧美大片| 一区二区三区免费毛片| 精品人妻1区二区| 日本一本二区三区精品| 一区二区三区四区激情视频 | 欧美另类亚洲清纯唯美| 色哟哟·www| 国产黄片美女视频| xxxwww97欧美| 亚州av有码| 国产毛片a区久久久久| 亚洲人与动物交配视频| 国产免费av片在线观看野外av| 好看av亚洲va欧美ⅴa在| 国产精华一区二区三区| 久久国产乱子免费精品| av专区在线播放| 亚洲精品影视一区二区三区av| 午夜免费成人在线视频| 九色国产91popny在线| 日日干狠狠操夜夜爽| 成人国产一区最新在线观看| 又爽又黄a免费视频| 日本免费a在线| 免费在线观看影片大全网站| 色尼玛亚洲综合影院| 免费观看人在逋| 每晚都被弄得嗷嗷叫到高潮| 国产精品永久免费网站| 欧美性感艳星| 我的老师免费观看完整版| 亚洲国产精品久久男人天堂| 久久久久久大精品| 国产亚洲精品久久久久久毛片| 久久久国产成人精品二区| 色哟哟·www| 熟女电影av网| 国产v大片淫在线免费观看| 欧美性猛交黑人性爽| 在线天堂最新版资源| 中出人妻视频一区二区| 亚洲无线在线观看| 国产精品亚洲一级av第二区| 特大巨黑吊av在线直播| 中文字幕熟女人妻在线| 成年女人看的毛片在线观看| 免费黄网站久久成人精品 |