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

    基于CEEMD與SVM的離心泵轉(zhuǎn)子不對(duì)中故障診斷方法研究

    2022-08-25 14:04:22肖幸鑫宋禮威張翊勛張宇航
    流體機(jī)械 2022年7期
    關(guān)鍵詞:離心泵特征參數(shù)時(shí)域

    肖幸鑫,宋禮威,張翊勛,董 亮,張宇航

    (1.江蘇大學(xué) 流體機(jī)械工程技術(shù)研究中心,江蘇鎮(zhèn)江 212000;2.中廣核工程有限公司 核電監(jiān)控技術(shù)與裝備國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣東深圳 518172)

    0 引言

    離心泵作為在各類(lèi)工程中常用的機(jī)械設(shè)備,廣泛應(yīng)用于例如石油,煤礦等能源的開(kāi)采;排水灌溉等農(nóng)業(yè)生產(chǎn)方向;航空航天航海等領(lǐng)域,占到泵產(chǎn)品總量的70%[1]。在實(shí)際工程應(yīng)用中,離心泵的轉(zhuǎn)子故障成為了困擾國(guó)內(nèi)外相關(guān)研究人員的難題[2];判斷離心泵是否發(fā)生故障,主要依靠離心泵在運(yùn)行過(guò)程中出現(xiàn)異常的噪聲與振動(dòng)[3]。

    離心泵轉(zhuǎn)子不對(duì)中故障主要由轉(zhuǎn)子間存在高度差或轉(zhuǎn)子間存在距離引起的,如圖1所示。

    圖1 典型的不對(duì)中故障Fig.1 Typical misalignment fault

    2018年李峰等[5]提出了基于電機(jī)電流時(shí)頻特征的不對(duì)中故障診斷方法,有效地提高了識(shí)別轉(zhuǎn)子系統(tǒng)的軸系不對(duì)中故障的正確率,因此通過(guò)提取信號(hào)時(shí)頻特征參數(shù),用于故障診斷是有效的;唐貴基等[6]提出一種譜峭度與變分模態(tài)分解的轉(zhuǎn)子故障診斷方法,解決了在噪聲環(huán)境下轉(zhuǎn)子微弱不對(duì)中故障特征難以提取的問(wèn)題,劉忠等[7]提出一種IVMD算法,并將其應(yīng)用到離心泵空化聲發(fā)射信號(hào)特征提取,因此模態(tài)分解算法廣泛用于特征提取中,降低故障特征參數(shù)提取的難度;XIE等[8]將不對(duì)中故障的特征參數(shù)離散和模糊化構(gòu)造故障決策表,能夠有效地診斷汽輪發(fā)電機(jī)組轉(zhuǎn)子不對(duì)中故障;YAO等[9]提出了基于分段閾值小波去噪、經(jīng)驗(yàn)?zāi)B(tài)分解、遺傳算法和最小二乘支持向量機(jī)的混合鐵路滾動(dòng)軸承故障診斷方法;肖軍等[10]針對(duì)往復(fù)壓縮機(jī)氣閥故障,將氣閥振動(dòng)信號(hào)不同頻段的能量作為故障特征,提出了結(jié)合模糊C均值聚類(lèi)和支持向量機(jī)的二叉樹(shù)多分類(lèi)診斷方法。

    小波分析與經(jīng)驗(yàn)?zāi)B(tài)分解(empirical moded ecomposition,簡(jiǎn)稱(chēng)EMD)系列方法作為信號(hào)處理的常見(jiàn)方法,小波分析在提取故障特征信息時(shí)需要選擇合適的小波基,不同的小波基對(duì)于故障特征信息的提取影響巨大,在使用小波分析時(shí)必須進(jìn)行多次試驗(yàn)來(lái)確定最合適的小波基,并且確定后不能更改,對(duì)信號(hào)不具備自適應(yīng)性;而EMD克服了基函數(shù)無(wú)法自適應(yīng)的問(wèn)題,可以將一段信號(hào)進(jìn)行固定模式開(kāi)始分解,從高頻至低頻拆分成若干個(gè)IMF分量。轉(zhuǎn)子不對(duì)中故障主要故障頻率位于低頻,因此使用EMD作為轉(zhuǎn)子不對(duì)中故障信號(hào)處理方法能夠更有效地提取故障特征[11]。SVM作為有監(jiān)督算法的典型算法之一,利用給定的訓(xùn)練數(shù)據(jù)集進(jìn)行訓(xùn)練學(xué)習(xí)從而得到函數(shù)模型,并根據(jù)函數(shù)模型預(yù)測(cè)新數(shù)據(jù)類(lèi)型,王金東等[12]將EMD信息熵和支持向量機(jī)(SVM)結(jié)合應(yīng)用于往復(fù)式壓縮機(jī)的軸承故障診斷,發(fā)現(xiàn)該方法能夠準(zhǔn)確識(shí)別壓縮機(jī)軸承故障;杜星洲[13]提出了基于振動(dòng)信號(hào)與SVM的旋轉(zhuǎn)機(jī)械故障診斷研究方法,利用SVM對(duì)旋轉(zhuǎn)機(jī)械的常見(jiàn)故障進(jìn)行分類(lèi)和識(shí)別;而無(wú)監(jiān)督算法對(duì)于樣本數(shù)據(jù)類(lèi)型未知,主要依靠數(shù)據(jù)間的相似性對(duì)數(shù)據(jù)進(jìn)行分類(lèi),黃惠中[14]利用聚類(lèi)分析為基礎(chǔ),構(gòu)建了包含信號(hào)處理與特征提取及特征選擇技術(shù)在內(nèi)的智能故障診斷模型;李志國(guó)等[15]針對(duì)潛水磨碎泵流道堵塞的故障診斷,提出了一種基于奇異值分解和固有時(shí)間尺度分解的電流信號(hào)分析和構(gòu)建支持向量機(jī)故障識(shí)別模型相結(jié)合的方法。

    本文主要是針對(duì)離心泵轉(zhuǎn)子不對(duì)中故障與正常狀態(tài)區(qū)分,是一個(gè)二分類(lèi)問(wèn)題,而SVM作為一個(gè)二分類(lèi)模型相較于其他神經(jīng)網(wǎng)絡(luò)模型具有較大的優(yōu)勢(shì)。本文利用互補(bǔ)經(jīng)驗(yàn)?zāi)B(tài)分解(CEEMD)對(duì)電渦流位移傳感器采集的振動(dòng)位移信號(hào)進(jìn)行分解,通過(guò)相關(guān)系數(shù)法與閾值選取有效IMF分量,通過(guò)重構(gòu)有效IMF分量,通過(guò)提取重構(gòu)后信號(hào)中的時(shí)域故障特征參數(shù)組成特征向量,再將該特征向量輸入SVM中進(jìn)行轉(zhuǎn)子不對(duì)中故障判斷。

    1 互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解

    EMD在處理非平穩(wěn)及非線性數(shù)據(jù)上,具有非常明顯的優(yōu)勢(shì),適合于分析非線性、非平穩(wěn)信號(hào)序列[16],具有很高的信噪比。

    CEEMD算法相較于EMD在故障特征提取中存在模態(tài)混疊問(wèn)題和集中經(jīng)驗(yàn)?zāi)B(tài)分解(ensemble empirical mode decomposition,簡(jiǎn)稱(chēng)EEMD)存在運(yùn)算量大,計(jì)算效率低的問(wèn)題,CEEMD算法具備迭代次數(shù)少,計(jì)算成本低,診斷精度高等優(yōu)點(diǎn)[17-20]。

    1.1 CCEMD基本原理

    CEEMD主要計(jì)算步驟如下:

    (1)向原始信號(hào)s(t)中不同的正負(fù)白噪聲x(t):

    xna(t)、xnb(t)分別為加入正負(fù)白噪聲后的信號(hào),其中n為正整數(shù),由此可以得到2n組IMF集合。

    (2)將得到的2n組IMF集合中的每個(gè)IMF分量使用EMD進(jìn)行分解,每個(gè)信號(hào)都可以得到一組IMF,將第i個(gè)信號(hào)的第j個(gè)IMF分量記為Cij。

    (3)將得到的2n組IMF進(jìn)行平均,從而得到最終的IMF分量。

    該算法需要相較于EMD需要添加2個(gè)參數(shù):輔助白噪聲幅值k和對(duì)數(shù)N,一般當(dāng)N取100時(shí),k 取 0.01~0.10[21]。

    1.2 相關(guān)系數(shù)法分析

    基于CEEMD的分解流程,可以知道從原信號(hào)分解的IMF中提取有效的IMF,是實(shí)現(xiàn)轉(zhuǎn)子不對(duì)中故障診斷的重要步驟。采用相關(guān)系數(shù)法,提取出與原信號(hào)高于閾值的有效IMF分量來(lái)重構(gòu)信號(hào),對(duì)重構(gòu)后的信號(hào)進(jìn)行時(shí)域特征提取。

    設(shè)信號(hào)采樣信號(hào)數(shù)為M,原始信號(hào)為S,第k個(gè)IMF分量為IMFk,則IMF分量與原信號(hào)相關(guān)系數(shù)的計(jì)算式為:

    式中 Rk——相關(guān)系數(shù);

    E ——平均。

    閾值TH定義為相關(guān)系數(shù)的標(biāo)準(zhǔn)差,計(jì)算式為:

    式中 n ——分解出來(lái)的IMF的個(gè)數(shù)。

    2 時(shí)域故障特征提取

    時(shí)域信息是以時(shí)間為變量描繪出信號(hào)的波形,時(shí)域信號(hào)包括有量綱特征參數(shù)與無(wú)量綱特征參數(shù),本文選用平均值和均方根、峭度、波形因子、峰值因子、脈動(dòng)因子、裕度因子7個(gè)時(shí)域特征參數(shù)進(jìn)行分析[17]。時(shí)域特征參數(shù)公式見(jiàn)表1。表1中,Ns表示采樣點(diǎn)數(shù),i=1~Ns;x(i)表示采集點(diǎn)數(shù)中第i個(gè)信號(hào)。

    表1 時(shí)域特征參數(shù)公式Tab.1 Time domain formula

    3 SVM原理

    SVM是一類(lèi)按監(jiān)督學(xué)習(xí)方式對(duì)數(shù)據(jù)進(jìn)行二元分類(lèi)的廣義線性分類(lèi)器,主要功能是解決模式識(shí)別區(qū)域的數(shù)據(jù)分類(lèi)等問(wèn)題,在解決小樣本、非線性及高維模式識(shí)別問(wèn)題中具有特有的優(yōu)勢(shì),基本原理在文獻(xiàn)[22-23]中有詳細(xì)描述。本文中不再贅述。

    4 轉(zhuǎn)子角不對(duì)中故障試驗(yàn)

    為了模擬離心泵角不對(duì)中故障,搭建離心泵故障模擬實(shí)驗(yàn)臺(tái),其主要設(shè)備包括臥式離心泵,三相異步電動(dòng)機(jī),電磁流量計(jì),電渦流振動(dòng)位移傳感器等。在試驗(yàn)中離心泵轉(zhuǎn)速通過(guò)變頻電機(jī)調(diào)節(jié)控制,系統(tǒng)流量通過(guò)出口閥門(mén)進(jìn)行調(diào)節(jié),離心泵轉(zhuǎn)子不對(duì)中故障通過(guò)在聯(lián)軸器放置墊片來(lái)模擬實(shí)現(xiàn)。

    試驗(yàn)使用的電渦流振動(dòng)位移傳感器參數(shù)見(jiàn)表2,利用電渦流振動(dòng)位移傳感器采集轉(zhuǎn)子振動(dòng)位移信號(hào),采樣頻率為25 600 Hz,離心泵的額定轉(zhuǎn)速為2 900 r/min,額定流量為10.6 m3/h。

    表2 電渦流傳感器參數(shù)Tab.2 Eddy current sensor parameters

    本次實(shí)驗(yàn)采集離心泵在同一流量下轉(zhuǎn)速比分別為0.7n,0.85n,1.0n的工況,正常運(yùn)轉(zhuǎn)狀態(tài)下的振動(dòng)位移信號(hào)和轉(zhuǎn)子角不對(duì)中故障振動(dòng)位移信號(hào);采集在同一轉(zhuǎn)速下,流量比分別為0.7q,0.85q,1.0q工況下的正常狀態(tài)下的振動(dòng)信號(hào)和角不對(duì)中故障振動(dòng)信號(hào)。

    試驗(yàn)主要針對(duì)轉(zhuǎn)子不對(duì)中故障中的角不對(duì)中故障進(jìn)行研究,如圖2,3所示,在電機(jī)端添加墊片來(lái)模擬實(shí)際工程中轉(zhuǎn)子角不對(duì)中故障的發(fā)生,利用百分表測(cè)量聯(lián)軸器位置的偏移距離來(lái)判斷是否故障為轉(zhuǎn)子角不對(duì)中。進(jìn)行多次測(cè)量后,求得平均值記錄見(jiàn)表3。

    表3 聯(lián)軸器偏移距離Tab.3 Coupling offset distance

    圖2 垂直方向偏移距離Fig.2 Vertical offset distance

    圖3 水平方向偏移距離Fig.3 Horizontal offset distance

    5 試驗(yàn)數(shù)據(jù)分析以及驗(yàn)證

    本文試驗(yàn)使用16通道采集卡,采集頻率為25 600 Hz,采集時(shí)間為1 s,每組試驗(yàn)采集樣本個(gè)數(shù)為 25 600個(gè);選取2 900 r/min,10.6 m3/h工況下其中10 240個(gè)樣本點(diǎn)繪制原始數(shù)據(jù)時(shí)域波形,如圖4,5所示。

    圖4 正常工況時(shí)域Fig.4 Time domain diagram of normal working conditions

    圖5 故障工況時(shí)域Fig.5 Time domain diagram of fault working conditions

    圖6示出了正常工況與故障工況頻譜對(duì)比,其中f'/f為倍頻,f'為不同工況的實(shí)際頻率,f是額定轉(zhuǎn)速下對(duì)應(yīng)的基頻(48.33 Hz),實(shí)驗(yàn)中通過(guò)對(duì)工常工況以及故障工況在額定轉(zhuǎn)速,額定流量下的頻譜圖進(jìn)行分析可以發(fā)現(xiàn),正常工況下的主頻出現(xiàn)在1APF,故障工況下的1APF為主頻,但是2APF的振幅明顯增大,符合轉(zhuǎn)子角不對(duì)中故障發(fā)生時(shí)的特征。

    圖6 正常工況與故障工況頻譜對(duì)比Fig.7 Spectrum comparison between normal working conditions and fault working conditions

    5.1 CEEMD分解以及信號(hào)重構(gòu)

    對(duì)正常工況下的額定轉(zhuǎn)速,額定流量下的信號(hào)進(jìn)行CEEMD分解,我們選取一組數(shù)據(jù),數(shù)據(jù)點(diǎn)個(gè)數(shù)為2 560個(gè),結(jié)果如圖7所示。

    圖7 CEEMD分解Fig.7 CEEMD decomposition diagram

    各IMF分量與原信號(hào)的相關(guān)系數(shù)見(jiàn)表4,為了方便表示,各相關(guān)系數(shù)均乘以100%。得到相關(guān)系數(shù)后,利用式(4)可以計(jì)算出閾值TH=18.18%,因此選取有效 IMF 分量為 IMF1,IMF2,IMF6,IMF7,進(jìn)行信號(hào)重構(gòu),從而得到一個(gè)新的信號(hào)NewIMF,NewIMF的時(shí)域波形如圖8所示。

    表4 相關(guān)系數(shù)Tab.4 Correlation coefficient

    圖8 重構(gòu)后NewIMF信號(hào)Fig.8 Reconstructed NewIMF signal

    5.2 故障特征提取

    選取離心泵正常運(yùn)行狀態(tài)信號(hào),離心泵泵轉(zhuǎn)子角不對(duì)中故障狀態(tài)信號(hào)各50組,一共100組信號(hào),每組2 560個(gè)數(shù)據(jù)點(diǎn),對(duì)每組信號(hào)都進(jìn)行上述方法分析,通過(guò)重構(gòu)后信號(hào),得到平均值和均方根2個(gè)有量綱時(shí)域特征參數(shù)以及峭度、波形因子、峰值因子、脈動(dòng)因子、裕度因子5個(gè)無(wú)量綱時(shí)域特征參數(shù),并對(duì)特征參數(shù)進(jìn)行min-max歸一化處理使數(shù)據(jù)映射至[-1,1]區(qū)間內(nèi):

    式中 y1——?dú)w一化后數(shù)據(jù);

    y ——原始樣本數(shù)據(jù);

    min ——樣本數(shù)據(jù)最小值;

    max ——樣本數(shù)據(jù)最大值。

    歸一化處理后組成特征向量,由于篇幅原因表中僅列出正常狀態(tài)下和故障狀態(tài)下各10組特征向量見(jiàn)表5。

    表5 部分正常與故障特征向量Tab.5 Part of the normal and fault characteristic vectors

    5.3 利用SVM進(jìn)行故障識(shí)別

    支持向量機(jī)選用RBF作為核函數(shù),利用交叉驗(yàn)證法選用最合適的懲罰系數(shù)c和參數(shù)g,在數(shù)據(jù)集中隨機(jī)選擇80組數(shù)據(jù),作為訓(xùn)練集,20組數(shù)據(jù)作為測(cè)試集。由于是隨機(jī)選擇數(shù)據(jù)作為訓(xùn)練集以及測(cè)試集,為了避免識(shí)別誤差,因此進(jìn)行20次識(shí)別,并取其平均值,由表可知平均值為93%。圖9示出第一次識(shí)別結(jié)果,將識(shí)別結(jié)果記錄見(jiàn)表6。

    圖9 第一次SVM故障識(shí)別結(jié)果Fig.9 Identification results of the first SVM fault

    表6 SVM識(shí)別結(jié)果Tab.6 Identification results of SVM

    6 結(jié)論

    (1)離心泵轉(zhuǎn)子不對(duì)中故障發(fā)生時(shí),轉(zhuǎn)子振動(dòng)信號(hào)的時(shí)域特征參數(shù)會(huì)發(fā)生變化,且在頻譜中二倍頻會(huì)明顯增大。

    (2)使用CEEMD算法進(jìn)行信號(hào)分解,采用相關(guān)系數(shù)法和閾值能夠有效的選擇出包含故障特征的IMF分量。

    (3)對(duì)重構(gòu)后的故障信號(hào)提取時(shí)域特征參數(shù),組成特征向量輸入SVM中,可有效的識(shí)別出離心泵轉(zhuǎn)子不對(duì)中故障,并具有較高的準(zhǔn)確率。

    猜你喜歡
    離心泵特征參數(shù)時(shí)域
    一種改善離心泵運(yùn)行狀態(tài)的方法
    大型立式單級(jí)引黃離心泵軸向力平衡的研究
    故障診斷中信號(hào)特征參數(shù)擇取方法
    離心泵流場(chǎng)外特性及激勵(lì)數(shù)值計(jì)算
    高比速離心泵的高效設(shè)計(jì)研究
    基于特征參數(shù)化的木工CAD/CAM系統(tǒng)
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    基于PSO-VMD的齒輪特征參數(shù)提取方法研究
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    两性午夜刺激爽爽歪歪视频在线观看 | 久久人人97超碰香蕉20202| 国产精品久久久av美女十八| 午夜福利乱码中文字幕| 亚洲国产精品一区二区三区在线| bbb黄色大片| 9色porny在线观看| 又黄又粗又硬又大视频| 好男人电影高清在线观看| 热99久久久久精品小说推荐| 亚洲精品日韩在线中文字幕| 国产成人免费观看mmmm| 国产在视频线精品| 欧美+亚洲+日韩+国产| 99精品欧美一区二区三区四区| 国产精品久久久久久人妻精品电影 | 在线观看免费午夜福利视频| 久久久久久久精品精品| 国产精品久久久久久精品古装| 国产精品久久久人人做人人爽| 国产免费福利视频在线观看| 亚洲综合色网址| 国产亚洲精品一区二区www | 中文字幕制服av| 午夜福利,免费看| 最新的欧美精品一区二区| 国产亚洲欧美精品永久| a 毛片基地| 无遮挡黄片免费观看| 久久精品人人爽人人爽视色| 一区二区av电影网| 天天躁日日躁夜夜躁夜夜| 亚洲专区中文字幕在线| 人人妻人人澡人人爽人人夜夜| 亚洲天堂av无毛| av天堂久久9| 老司机影院成人| 在线观看人妻少妇| 成年美女黄网站色视频大全免费| 熟女少妇亚洲综合色aaa.| 亚洲欧美清纯卡通| 免费少妇av软件| 俄罗斯特黄特色一大片| 在线av久久热| 9热在线视频观看99| 午夜影院在线不卡| 天堂中文最新版在线下载| 首页视频小说图片口味搜索| 曰老女人黄片| 免费黄频网站在线观看国产| 亚洲欧美精品自产自拍| 熟女少妇亚洲综合色aaa.| a级毛片黄视频| 午夜影院在线不卡| 电影成人av| 大片电影免费在线观看免费| 色婷婷久久久亚洲欧美| 美女午夜性视频免费| 亚洲精品中文字幕一二三四区 | 无限看片的www在线观看| 99久久人妻综合| 国产一区二区激情短视频 | 免费黄频网站在线观看国产| 久久人妻福利社区极品人妻图片| 成人三级做爰电影| 亚洲国产精品成人久久小说| 欧美午夜高清在线| 国产欧美日韩一区二区精品| 99久久综合免费| 伊人久久大香线蕉亚洲五| www.av在线官网国产| tocl精华| 啦啦啦视频在线资源免费观看| 少妇粗大呻吟视频| 欧美激情极品国产一区二区三区| 精品一区二区三区四区五区乱码| 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 精品福利永久在线观看| 2018国产大陆天天弄谢| 欧美日韩精品网址| 午夜福利,免费看| 中文欧美无线码| 久久综合国产亚洲精品| 超碰97精品在线观看| 亚洲成av片中文字幕在线观看| 中文欧美无线码| 精品国内亚洲2022精品成人 | av国产精品久久久久影院| 国产成人欧美| 欧美日韩中文字幕国产精品一区二区三区 | 天天操日日干夜夜撸| 色综合欧美亚洲国产小说| www日本在线高清视频| 巨乳人妻的诱惑在线观看| 男女之事视频高清在线观看| 国产片内射在线| 亚洲精品第二区| 国产免费av片在线观看野外av| 欧美亚洲 丝袜 人妻 在线| 丰满迷人的少妇在线观看| 久久久精品区二区三区| 色94色欧美一区二区| 亚洲欧美成人综合另类久久久| 丝袜美足系列| 波多野结衣av一区二区av| 欧美xxⅹ黑人| 欧美另类一区| 亚洲国产成人一精品久久久| 中文字幕另类日韩欧美亚洲嫩草| 欧美中文综合在线视频| 亚洲精品国产区一区二| 国产精品二区激情视频| 热99re8久久精品国产| 精品国内亚洲2022精品成人 | 久久精品国产综合久久久| av在线app专区| 免费观看人在逋| 91麻豆av在线| 亚洲精品在线美女| 视频区欧美日本亚洲| 黑人巨大精品欧美一区二区蜜桃| 精品久久久久久久毛片微露脸 | 在线亚洲精品国产二区图片欧美| 中文字幕另类日韩欧美亚洲嫩草| 2018国产大陆天天弄谢| 黄网站色视频无遮挡免费观看| 一本一本久久a久久精品综合妖精| 日本av免费视频播放| 精品高清国产在线一区| 免费少妇av软件| 视频在线观看一区二区三区| 99久久99久久久精品蜜桃| 涩涩av久久男人的天堂| 这个男人来自地球电影免费观看| 丝袜美足系列| 他把我摸到了高潮在线观看 | 高潮久久久久久久久久久不卡| 久久久久久人人人人人| 美女大奶头黄色视频| 亚洲精品国产色婷婷电影| 精品国产超薄肉色丝袜足j| www.av在线官网国产| 日本一区二区免费在线视频| 欧美久久黑人一区二区| a级毛片黄视频| 国产真人三级小视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | av天堂久久9| 久久国产精品影院| av线在线观看网站| 在线天堂中文资源库| 国产色视频综合| 91大片在线观看| 热re99久久精品国产66热6| 国产xxxxx性猛交| 狠狠狠狠99中文字幕| 极品人妻少妇av视频| 精品少妇内射三级| 97在线人人人人妻| 亚洲人成77777在线视频| 黄色片一级片一级黄色片| 欧美久久黑人一区二区| 下体分泌物呈黄色| 新久久久久国产一级毛片| 国产精品二区激情视频| 国产成人精品久久二区二区91| 亚洲人成电影免费在线| 中国国产av一级| 老司机影院成人| 亚洲少妇的诱惑av| 国产在线视频一区二区| 超色免费av| 美女扒开内裤让男人捅视频| av国产精品久久久久影院| 下体分泌物呈黄色| 少妇精品久久久久久久| 精品久久久久久久毛片微露脸 | 黄片播放在线免费| 免费少妇av软件| 国产精品成人在线| 不卡av一区二区三区| 欧美变态另类bdsm刘玥| 人人妻人人澡人人看| 一本大道久久a久久精品| 久热这里只有精品99| 亚洲天堂av无毛| 老司机深夜福利视频在线观看 | 国产成人a∨麻豆精品| 亚洲国产看品久久| 日本一区二区免费在线视频| 曰老女人黄片| 天天添夜夜摸| 亚洲情色 制服丝袜| 日韩欧美国产一区二区入口| 亚洲一区二区三区欧美精品| 99国产精品99久久久久| 久久av网站| 一区二区av电影网| 亚洲精品第二区| 久久女婷五月综合色啪小说| 欧美精品av麻豆av| 丝袜喷水一区| 99香蕉大伊视频| 欧美 亚洲 国产 日韩一| 人妻一区二区av| 99久久精品国产亚洲精品| 亚洲欧美精品综合一区二区三区| 91成年电影在线观看| av网站在线播放免费| 飞空精品影院首页| 99re6热这里在线精品视频| 久久天躁狠狠躁夜夜2o2o| 多毛熟女@视频| 久久久久久久国产电影| 这个男人来自地球电影免费观看| 国产精品熟女久久久久浪| 在线精品无人区一区二区三| 亚洲第一欧美日韩一区二区三区 | 亚洲av成人不卡在线观看播放网 | 亚洲国产精品成人久久小说| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲专区国产一区二区| 久久久国产精品麻豆| 国产亚洲欧美精品永久| 日本五十路高清| 精品福利永久在线观看| √禁漫天堂资源中文www| 在线观看www视频免费| 久久人妻福利社区极品人妻图片| 国产免费视频播放在线视频| 99精品欧美一区二区三区四区| 精品久久蜜臀av无| 欧美+亚洲+日韩+国产| 麻豆av在线久日| 中文字幕av电影在线播放| 日本91视频免费播放| 在线观看免费日韩欧美大片| av免费在线观看网站| 国产一区二区三区在线臀色熟女 | 免费高清在线观看日韩| 色视频在线一区二区三区| tube8黄色片| 香蕉国产在线看| 99国产精品一区二区三区| 中文字幕人妻丝袜一区二区| 国产精品 欧美亚洲| 日本五十路高清| 国产精品成人在线| 久久国产精品影院| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| a级片在线免费高清观看视频| 免费观看a级毛片全部| 亚洲欧美日韩高清在线视频 | 国产成人精品在线电影| 一区二区三区四区激情视频| 久久亚洲国产成人精品v| 777米奇影视久久| 少妇粗大呻吟视频| 亚洲美女黄色视频免费看| 9热在线视频观看99| 国产精品自产拍在线观看55亚洲 | 国产成人啪精品午夜网站| 成人国产一区最新在线观看| 最黄视频免费看| 国产日韩欧美视频二区| 亚洲av日韩在线播放| 制服人妻中文乱码| 国产免费av片在线观看野外av| 成人国语在线视频| 精品人妻1区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩 亚洲 欧美在线| 国产亚洲精品久久久久5区| 看免费av毛片| 一二三四社区在线视频社区8| 欧美少妇被猛烈插入视频| 免费黄频网站在线观看国产| 色婷婷av一区二区三区视频| 在线亚洲精品国产二区图片欧美| 久9热在线精品视频| 777久久人妻少妇嫩草av网站| 国产精品自产拍在线观看55亚洲 | 高清在线国产一区| 人妻 亚洲 视频| 成人国语在线视频| 91成人精品电影| 免费在线观看黄色视频的| 韩国高清视频一区二区三区| 视频在线观看一区二区三区| 一二三四在线观看免费中文在| 精品免费久久久久久久清纯 | 精品国产乱子伦一区二区三区 | 午夜福利视频精品| 欧美日本中文国产一区发布| 久久综合国产亚洲精品| 久久天堂一区二区三区四区| 国产亚洲午夜精品一区二区久久| 日本欧美视频一区| 亚洲精品乱久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人免费观看mmmm| 好男人电影高清在线观看| 如日韩欧美国产精品一区二区三区| 一边摸一边抽搐一进一出视频| 黑丝袜美女国产一区| 亚洲 欧美一区二区三区| 欧美老熟妇乱子伦牲交| 亚洲精品国产色婷婷电影| 精品视频人人做人人爽| 免费久久久久久久精品成人欧美视频| 十八禁高潮呻吟视频| 女警被强在线播放| 成人国产一区最新在线观看| 精品国产一区二区久久| 国产在线一区二区三区精| 成人黄色视频免费在线看| 又紧又爽又黄一区二区| 久久久久久久久久久久大奶| 曰老女人黄片| 国产免费福利视频在线观看| 国产老妇伦熟女老妇高清| 久久女婷五月综合色啪小说| netflix在线观看网站| 亚洲久久久国产精品| 丰满人妻熟妇乱又伦精品不卡| 美女午夜性视频免费| 十分钟在线观看高清视频www| 日韩欧美国产一区二区入口| 日韩一卡2卡3卡4卡2021年| 999久久久国产精品视频| 久久久国产成人免费| 另类精品久久| 成年人黄色毛片网站| 亚洲中文字幕日韩| 国产高清国产精品国产三级| 亚洲成人国产一区在线观看| 久久人人爽av亚洲精品天堂| 亚洲精品美女久久久久99蜜臀| 欧美久久黑人一区二区| 国产不卡av网站在线观看| 大片电影免费在线观看免费| 在线观看一区二区三区激情| 桃花免费在线播放| 新久久久久国产一级毛片| 亚洲国产av新网站| 人人妻,人人澡人人爽秒播| 一区二区av电影网| av在线app专区| 日韩一卡2卡3卡4卡2021年| 亚洲av成人不卡在线观看播放网 | 欧美97在线视频| 91麻豆av在线| 一区二区日韩欧美中文字幕| 中文字幕最新亚洲高清| 99精品欧美一区二区三区四区| 欧美久久黑人一区二区| 国产欧美日韩精品亚洲av| 一级毛片电影观看| 伊人久久大香线蕉亚洲五| 电影成人av| 久久精品国产综合久久久| 国产免费现黄频在线看| 国精品久久久久久国模美| 19禁男女啪啪无遮挡网站| 50天的宝宝边吃奶边哭怎么回事| www.av在线官网国产| 亚洲精品av麻豆狂野| 在线精品无人区一区二区三| 黄色毛片三级朝国网站| 日本黄色日本黄色录像| 亚洲精品av麻豆狂野| 五月开心婷婷网| 淫妇啪啪啪对白视频 | 男女国产视频网站| 日韩欧美一区二区三区在线观看 | 一级片免费观看大全| 成年人午夜在线观看视频| 人成视频在线观看免费观看| 国产精品1区2区在线观看. | 一区二区日韩欧美中文字幕| 中文字幕最新亚洲高清| 99国产精品免费福利视频| 看免费av毛片| 精品久久久久久久毛片微露脸 | 啦啦啦 在线观看视频| 日本av手机在线免费观看| 亚洲男人天堂网一区| 国产一区二区 视频在线| 91精品三级在线观看| 50天的宝宝边吃奶边哭怎么回事| av有码第一页| 午夜福利影视在线免费观看| 男女下面插进去视频免费观看| 精品福利观看| 国产黄频视频在线观看| 亚洲国产欧美日韩在线播放| 国产一区二区在线观看av| 欧美+亚洲+日韩+国产| 欧美老熟妇乱子伦牲交| 亚洲情色 制服丝袜| 十八禁网站免费在线| 日本av手机在线免费观看| 19禁男女啪啪无遮挡网站| 国产成人精品久久二区二区免费| 少妇裸体淫交视频免费看高清 | 成人免费观看视频高清| 亚洲欧美一区二区三区久久| 亚洲一码二码三码区别大吗| 日韩 亚洲 欧美在线| 欧美日韩成人在线一区二区| 在线亚洲精品国产二区图片欧美| 国产成人免费观看mmmm| 日韩精品免费视频一区二区三区| 制服人妻中文乱码| 美女高潮到喷水免费观看| 黑人操中国人逼视频| 久久中文字幕一级| 九色亚洲精品在线播放| 桃红色精品国产亚洲av| 黄色 视频免费看| 欧美日韩福利视频一区二区| 国产在线一区二区三区精| 18禁裸乳无遮挡动漫免费视频| 精品一区二区三卡| 精品一区在线观看国产| 日本黄色日本黄色录像| 美女福利国产在线| 亚洲人成77777在线视频| av有码第一页| 精品国产乱子伦一区二区三区 | 下体分泌物呈黄色| 国产色视频综合| 久久久久久久国产电影| 少妇的丰满在线观看| 精品国产一区二区三区四区第35| 中文字幕最新亚洲高清| 1024香蕉在线观看| 热99re8久久精品国产| 乱人伦中国视频| 亚洲成国产人片在线观看| 男人添女人高潮全过程视频| 欧美久久黑人一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 视频区欧美日本亚洲| 丰满少妇做爰视频| 满18在线观看网站| 一边摸一边做爽爽视频免费| 日本一区二区免费在线视频| 日本欧美视频一区| 青春草亚洲视频在线观看| 岛国毛片在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 免费不卡黄色视频| tube8黄色片| 欧美少妇被猛烈插入视频| 精品少妇黑人巨大在线播放| 亚洲国产欧美网| 国产精品av久久久久免费| 老熟妇乱子伦视频在线观看 | 9热在线视频观看99| 桃红色精品国产亚洲av| 久久精品国产a三级三级三级| 91麻豆av在线| 深夜精品福利| 91精品伊人久久大香线蕉| 欧美xxⅹ黑人| 麻豆乱淫一区二区| 亚洲国产av新网站| 久久亚洲精品不卡| av国产精品久久久久影院| 国产有黄有色有爽视频| 久久人人爽人人片av| 国产欧美日韩一区二区精品| 久久久国产欧美日韩av| 在线av久久热| tube8黄色片| 一本—道久久a久久精品蜜桃钙片| 黑人巨大精品欧美一区二区蜜桃| 18禁黄网站禁片午夜丰满| 久久久久久久国产电影| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久人人做人人爽| 欧美国产精品一级二级三级| 国产三级黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美 亚洲 国产 日韩一| 亚洲一区二区三区欧美精品| 嫩草影视91久久| 高清av免费在线| 欧美黄色片欧美黄色片| 91精品伊人久久大香线蕉| 国产福利在线免费观看视频| 国产一级毛片在线| 欧美黄色淫秽网站| 国产av一区二区精品久久| 精品少妇久久久久久888优播| 久久中文字幕一级| 亚洲熟女毛片儿| 精品久久久久久电影网| 久久毛片免费看一区二区三区| 国产成人免费无遮挡视频| 最近最新中文字幕大全免费视频| 老汉色∧v一级毛片| 一本色道久久久久久精品综合| 国产av一区二区精品久久| 一级毛片精品| 男女免费视频国产| 免费在线观看影片大全网站| 一本色道久久久久久精品综合| 欧美老熟妇乱子伦牲交| 人妻人人澡人人爽人人| 一区二区三区精品91| 手机成人av网站| 如日韩欧美国产精品一区二区三区| 精品人妻1区二区| 亚洲欧洲精品一区二区精品久久久| 国产免费现黄频在线看| 亚洲精品一卡2卡三卡4卡5卡 | 动漫黄色视频在线观看| 欧美黄色淫秽网站| 人妻一区二区av| 一本久久精品| 国产成人影院久久av| 伊人久久大香线蕉亚洲五| 成人亚洲精品一区在线观看| 9热在线视频观看99| 欧美日韩福利视频一区二区| 在线观看免费视频网站a站| 久久人人爽av亚洲精品天堂| 免费女性裸体啪啪无遮挡网站| 久久国产精品影院| 操美女的视频在线观看| 久久av网站| 中文字幕精品免费在线观看视频| 成年人黄色毛片网站| 国产成人系列免费观看| 久久综合国产亚洲精品| 51午夜福利影视在线观看| 久久久精品区二区三区| 成人国产av品久久久| 香蕉丝袜av| 国产亚洲一区二区精品| 夜夜夜夜夜久久久久| 久久久欧美国产精品| 在线观看舔阴道视频| 亚洲人成电影观看| 亚洲国产毛片av蜜桃av| 黑人欧美特级aaaaaa片| 叶爱在线成人免费视频播放| 另类精品久久| av网站在线播放免费| 99国产精品免费福利视频| 五月天丁香电影| 99热国产这里只有精品6| 亚洲国产看品久久| 91精品伊人久久大香线蕉| 亚洲精品自拍成人| 十八禁高潮呻吟视频| 两个人看的免费小视频| 亚洲欧美一区二区三区黑人| 亚洲精品第二区| 桃花免费在线播放| 十八禁网站免费在线| 国产精品久久久久久精品电影小说| 18禁黄网站禁片午夜丰满| 久久毛片免费看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品自拍成人| 亚洲成人免费av在线播放| 我的亚洲天堂| 在线看a的网站| 好男人电影高清在线观看| 亚洲欧洲精品一区二区精品久久久| 丝袜美腿诱惑在线| 老司机靠b影院| 欧美日韩亚洲国产一区二区在线观看 | 9191精品国产免费久久| 最近最新免费中文字幕在线| 国产日韩欧美亚洲二区| 国产日韩一区二区三区精品不卡| 久久人妻福利社区极品人妻图片| 丝袜在线中文字幕| 欧美少妇被猛烈插入视频| 亚洲九九香蕉| 免费在线观看影片大全网站| 欧美国产精品一级二级三级| 国产精品av久久久久免费| 老汉色∧v一级毛片| 国产精品欧美亚洲77777| 午夜福利在线免费观看网站| 麻豆乱淫一区二区| 777米奇影视久久| 成人亚洲精品一区在线观看| 韩国高清视频一区二区三区| 高清av免费在线| 天天添夜夜摸| 丝袜在线中文字幕| 久久久久精品国产欧美久久久 | 男男h啪啪无遮挡| 欧美日韩成人在线一区二区| 日韩 欧美 亚洲 中文字幕| 国产免费视频播放在线视频| 免费高清在线观看视频在线观看| 国产xxxxx性猛交| 亚洲国产欧美网| 波多野结衣一区麻豆| avwww免费| 国产亚洲精品一区二区www | 国产黄色免费在线视频| 午夜影院在线不卡| 丝袜美足系列| 91国产中文字幕|