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

    基于小波包模糊熵的往復(fù)壓縮機(jī)軸承故障特征提取方法研究

    2018-08-30 16:10:04趙海洋王金東
    噪聲與振動(dòng)控制 2018年4期
    關(guān)鍵詞:波包壓縮機(jī)尺度

    趙海洋,韓 輝,王金東,李 穎

    (東北石油大學(xué) 機(jī)械科學(xué)與工程學(xué)院,黑龍江 大慶 163318)

    往復(fù)壓縮機(jī)是用于壓縮和輸送氣體的機(jī)械設(shè)備,因其壓力適用范圍廣和適用性強(qiáng)等特點(diǎn),廣泛應(yīng)用于石油、化工行業(yè)。傳動(dòng)機(jī)構(gòu)是往復(fù)壓縮機(jī)動(dòng)力傳遞和運(yùn)動(dòng)形式轉(zhuǎn)換的重要部件,其連桿和各部件之間通常使用滑動(dòng)軸承連接,運(yùn)動(dòng)時(shí)間一久,滑動(dòng)軸承常因磨損而出現(xiàn)間隙過大故障[1]。因此,對(duì)其進(jìn)行故障診斷方法研究意義重大。往復(fù)壓縮機(jī)運(yùn)行時(shí),因其結(jié)構(gòu)復(fù)雜,激勵(lì)源眾多,運(yùn)動(dòng)形式多樣,其振動(dòng)信號(hào)呈現(xiàn)出強(qiáng)非平穩(wěn)性、非線性和特征耦合特性。這就導(dǎo)致用傳統(tǒng)的以線性系統(tǒng)為前提的時(shí)域和頻域信號(hào)處理技術(shù)對(duì)往復(fù)壓縮機(jī)軸承間隙狀態(tài)進(jìn)行故障診斷和評(píng)估存在一定的局限性。

    近年來,隨著非線性科學(xué)理論的發(fā)展,分形維數(shù)、近似熵、樣本熵和多尺度熵等非線性動(dòng)力學(xué)方法的提出為識(shí)別和預(yù)測(cè)復(fù)雜非線性非平穩(wěn)問題提供了解決思路。例如,基于熵的參量的非線性分析方法已成為故障特征提取領(lǐng)域的研究熱點(diǎn),并應(yīng)用于機(jī)械故障診斷領(lǐng)域。胥永剛等將近似熵應(yīng)用在機(jī)械設(shè)備的故障診斷中,結(jié)果表明近似熵在表征信號(hào)復(fù)雜性方面具有很強(qiáng)的能力[2];鄭近德等將樣本熵和多尺度熵應(yīng)用于滾動(dòng)軸承和轉(zhuǎn)子系統(tǒng)故障診斷,取得了良好的效果[3–4]。

    Pincus為統(tǒng)計(jì)量化非線性時(shí)間序列復(fù)雜性最早提出近似熵,并將其應(yīng)用于生理信號(hào)分析[5]。隨后,Richman等提出另一種度量序列復(fù)雜度的新方法—樣本熵,相比于近似熵等非線性動(dòng)力學(xué)指標(biāo),其具有所需數(shù)據(jù)長(zhǎng)度少、對(duì)噪聲不敏感、在大取值范圍一致性好等特點(diǎn)[6]。近年來,陳偉婷等提出了模糊熵的概念,采用指數(shù)模糊函數(shù)定義相似性克服了近似熵和樣本熵基于階躍函數(shù)定義相似性與實(shí)際樣本類邊緣不符的缺陷,避免了二分類性質(zhì)[7]。然而,近似熵、樣本熵和模糊熵只能從單一尺度上衡量時(shí)間序列復(fù)雜性,與復(fù)雜度沒有直接的關(guān)系,Costa等通過對(duì)時(shí)間序列進(jìn)行粗?;?,提出了多尺度熵的概念[8–9]。鄭近德將粗粒分割與模糊熵結(jié)合提出多尺度模糊熵,并將其應(yīng)用于不同軸承不同類型故障診斷中[10]。黃曉琳等指出實(shí)現(xiàn)多尺度分析的粗?;椒ū举|(zhì)是對(duì)時(shí)間序列做低通濾波[11]。因此,通過粗粒分割的多尺度模糊熵存在只分析時(shí)間序列的低頻成分而忽略高頻成分的局限。

    針對(duì)上述問題,本文在模糊熵的基礎(chǔ)上,利用小波包將頻帶多層次劃分優(yōu)勢(shì),提出小波包模糊熵的概念,用來衡量時(shí)間序列在不同頻率下的復(fù)雜性。與多尺度模糊熵相比,小波包模糊熵既可反映信號(hào)高頻分量的復(fù)雜度又可反映低頻分量的復(fù)雜度。為了能夠更加全面分析非線性、非平穩(wěn)信號(hào)在各個(gè)頻段的特征,本文將小波包模糊熵作為故障特征提取的工具應(yīng)用到往復(fù)壓縮機(jī)故障診斷中,利用小波包模糊熵提取的故障特征,結(jié)合支持向量機(jī)提出一種新的往復(fù)圧縮機(jī)軸承診斷方法。在此基礎(chǔ)上,將該方法與多尺度模糊熵進(jìn)行比較,實(shí)驗(yàn)結(jié)果表明該方法可有效完成軸承不同故障類型的準(zhǔn)確識(shí)別診斷。

    1 小波包模糊熵方法

    1.1 模糊熵

    (1)對(duì)一個(gè)N點(diǎn)的時(shí)間序列u(1),u(2),…u(N),按順序構(gòu)造一個(gè)m維向量

    代表連續(xù)m個(gè)u的值去掉均值u0(i)。

    (4)定義函數(shù)

    (5)類似地再對(duì)維數(shù)m+1重復(fù)上述步驟(1)至(4),得

    (6)定義模糊熵為

    當(dāng)N為有限數(shù)時(shí),上式表示成

    模糊熵同樣本熵一樣都是時(shí)間序列復(fù)雜性的度量,時(shí)間序列越復(fù)雜,熵值越大。然而,模糊熵不僅具備了樣本熵所需數(shù)據(jù)量小、保持一致性的特點(diǎn),而且在對(duì)數(shù)據(jù)處理過程中在細(xì)節(jié)和性能方面更為優(yōu)越。其一,模糊熵所采用的指數(shù)函數(shù)模糊化取代樣本熵中基于單位階躍函數(shù)定義向量間相似性的方式保證了模糊熵值隨參數(shù)變化且平穩(wěn)連續(xù)變化。其二,模糊熵中的均值運(yùn)算消除了樣本熵中通過絕對(duì)幅值差計(jì)算向量相似性所帶來的基線漂移效應(yīng),從而使相似性度量模糊化。

    1.2 多尺度模糊熵

    多尺度模糊熵定義為時(shí)間序列不同尺度的模糊熵,計(jì)算過程如下:

    (1)對(duì)長(zhǎng)度為N的原始時(shí)間序列{u(i):1≤i≤N},預(yù)先設(shè)定嵌入維數(shù)m和相似容限r(nóng),通過粗?;瘶?gòu)造粗粒序列

    其中τ為尺度因子,{yτ(j)}就是先通過把時(shí)間序列u(i)分成長(zhǎng)度為τ的窗口,然后計(jì)算窗口里數(shù)據(jù)的均值得到的。當(dāng)τ=1時(shí),{yτ(1)}就是原始序列,每個(gè)粗粒序列的長(zhǎng)度等于原時(shí)間序列的長(zhǎng)度除以尺度因子τ。序列粗粒化過程如圖1所示。

    (2)計(jì)算不同尺度因子τ個(gè)粗粒序列的模糊熵,觀察原時(shí)間序列不同尺度下模糊熵變化,稱為多尺度模糊熵分析。

    圖1 時(shí)間序列粗?;^程示意圖

    多尺度模糊熵通過對(duì)時(shí)間序列粗?;?,計(jì)算其在不同尺度下的模糊熵值,充分反映了時(shí)間序列在多個(gè)尺度上的復(fù)雜程度。同多尺度熵類似,如果一個(gè)時(shí)間序列模糊熵值在大部分尺度上都比另一個(gè)時(shí)間序列熵值大,則認(rèn)為前者比后者復(fù)雜性高。

    1.3 小波包模糊熵

    小波包模糊熵計(jì)算步驟如下:

    (1)對(duì)一個(gè)N點(diǎn)的時(shí)間序列{u(i):1≤i≤N}進(jìn)行小波包k層分解,分別提取第k層從低頻到高頻2k個(gè)頻率成分的節(jié)點(diǎn)信號(hào)特征,分解樹結(jié)構(gòu)如圖2所示。

    圖2 時(shí)間序列u(i)小波包分解樹結(jié)構(gòu)示意圖(k=3)

    (2)計(jì)算每個(gè)節(jié)點(diǎn)信號(hào)的模糊熵值,并將其表示為節(jié)點(diǎn)的函數(shù),觀察原時(shí)間序列在低頻和高頻不同頻率段的模糊熵變化,稱為小波包模糊熵分析。

    小波包分解和粗?;际嵌喑叨确治鰧?shí)現(xiàn)方式。然而,用通過粗?;姆绞接?jì)算得到的多尺度熵衡量時(shí)間序列時(shí)只考慮了時(shí)間序列的低頻信息,忽略了高頻信息。這樣對(duì)于故障信息主要頻率分布較豐富的時(shí)間序列,多尺度模糊熵不能滿足要求。而小波包分析[12]將時(shí)間序列分解成低頻和高頻,對(duì)頻帶進(jìn)行多層次劃分,并能根據(jù)信號(hào)特征自適應(yīng)選擇頻帶,與粗?;亩喑叨确椒ㄏ啾?,既考慮時(shí)間序列的高頻信息又考慮了時(shí)間的低頻信息,避免了遺漏在高頻成分中的重要信息。所以,用小波模糊熵分析往復(fù)壓縮機(jī)軸承故障振動(dòng)信號(hào)能夠提取出比多尺度熵更豐富的故障特征信息。

    1.4 參數(shù)選取

    根據(jù)小波包模糊熵的定義,在進(jìn)行小波包模糊熵計(jì)算前需要設(shè)定一個(gè)小波基函數(shù)和4個(gè)參數(shù),分別為嵌入維數(shù)m、模糊函數(shù)的相似容限r(nóng)、模糊函數(shù)的邊界梯度n、小波包分解的層數(shù)k。

    (1)基于往復(fù)壓縮機(jī)軸承間隙振動(dòng)信號(hào)特性,采用db10小波基函數(shù)。

    (2)嵌入維數(shù)m的大小影響序列的聯(lián)合概率動(dòng)態(tài)重構(gòu)時(shí)信息量的多少,m過小則信息會(huì)丟失,m過大則會(huì)有很多詳細(xì)信息,因此綜合考慮設(shè)定m=2。

    (3)相似容限r(nóng)的大小影響統(tǒng)計(jì)特性,r過大則統(tǒng)計(jì)信息丟失,r過小則統(tǒng)計(jì)特性不準(zhǔn)確,因此r一般取0.1~0.25SD(SD為原始序列的標(biāo)準(zhǔn)差)。

    (4)模糊函數(shù)邊界梯度n在向量間的相似性計(jì)算中起著權(quán)重的作用,n>1時(shí),更多計(jì)入較近向量對(duì)其相似度的貢獻(xiàn),而更少計(jì)入較遠(yuǎn)向量對(duì)其相似性的貢獻(xiàn);n過大會(huì)導(dǎo)致細(xì)節(jié)信息丟失,而n<1時(shí)則相反。為獲取更多細(xì)節(jié)信息,參考文獻(xiàn)[7],取較小整數(shù)值n=2。

    (5)分解層數(shù)k值過大影響計(jì)算效率并且會(huì)導(dǎo)致節(jié)點(diǎn)分量計(jì)算的點(diǎn)減少,k值過小則信號(hào)頻帶劃分不夠詳細(xì),無法獲取足夠的從低頻到高頻的節(jié)點(diǎn)分量,綜合考慮取k=3。

    2 基于小波包模糊熵的故障特征提取方法

    基于小波包模糊熵分析的故障特征提取過程如下:

    (1)信號(hào)預(yù)處理,實(shí)現(xiàn)對(duì)往復(fù)壓縮機(jī)軸承振動(dòng)信號(hào)的降噪處理。

    (2)對(duì)處理后的信號(hào)進(jìn)行小波包分解,得到8個(gè)分解節(jié)點(diǎn)信號(hào)。

    (3)將各節(jié)點(diǎn)信號(hào)的模糊熵FuzzyEn(u3,0)、FuzzyEn(u3,1)、…、FuzzyEn(u3,7)作為故障特征向量。

    3 往復(fù)壓縮機(jī)軸承振動(dòng)信號(hào)特征提取實(shí)驗(yàn)

    3.1 故障試驗(yàn)研究

    本試驗(yàn)以2D12型雙作用對(duì)動(dòng)式往復(fù)壓縮機(jī)為研究對(duì)象。其主要參數(shù)如下:軸功率為500 kW、排氣量為活塞行程為240 mm、電機(jī)轉(zhuǎn)速為根據(jù)壓縮機(jī)結(jié)構(gòu)特點(diǎn),設(shè)置測(cè)點(diǎn)在曲軸箱正面、上端和右端、各級(jí)氣缸蓋側(cè)、各級(jí)十字頭滑道下端7個(gè)故障敏感部位,利用加速度傳感器測(cè)試壓縮機(jī)振動(dòng)信號(hào),測(cè)點(diǎn)布局如圖3所示。

    圖3 往復(fù)壓縮機(jī)測(cè)點(diǎn)布局

    鑒于往復(fù)壓縮機(jī)激勵(lì)源眾多,運(yùn)行時(shí)會(huì)產(chǎn)生高頻振動(dòng)成分,采樣頻率設(shè)置為50kHz。數(shù)據(jù)采集系統(tǒng)選用北京東方所的INV306U-6660多通道智能數(shù)據(jù)采集儀,采樣時(shí)間為4 s。試驗(yàn)共模擬了往復(fù)壓縮機(jī)傳動(dòng)機(jī)構(gòu)軸承正常間隙狀態(tài)以及兩級(jí)連桿的小頭軸承、大頭軸承共四個(gè)位置的間隙大故障狀態(tài)。

    3.2 實(shí)驗(yàn)信號(hào)分析

    分別提取兩個(gè)周期的5種往復(fù)壓縮機(jī)軸承狀態(tài)的振動(dòng)加速度信號(hào)數(shù)據(jù)進(jìn)行分析,由采樣時(shí)間和電機(jī)轉(zhuǎn)速可知,其數(shù)據(jù)長(zhǎng)度為12 028點(diǎn)。眾所周知,噪聲會(huì)影響信號(hào)的復(fù)雜性,從而影響模糊熵的大小,使其無法準(zhǔn)確反映系統(tǒng)的復(fù)雜程度。為了降低噪聲干擾,應(yīng)用局部均值分解方法[13]對(duì)往復(fù)壓縮機(jī)軸承振動(dòng)數(shù)據(jù)進(jìn)行分解預(yù)處理,利用原始信號(hào)與各分解分量之間的互相關(guān)系數(shù)選擇包含主要故障信息的分量進(jìn)行重構(gòu),從而實(shí)現(xiàn)降噪處理。圖4至圖8所示分別為往復(fù)壓縮機(jī)軸承五種不同狀態(tài)數(shù)據(jù)降噪前后時(shí)域波形圖。

    由上述信號(hào)預(yù)處理結(jié)果可知,無論降噪前后,都很難從軸承故障信號(hào)的時(shí)域波形中判斷出軸承屬于何種狀態(tài),因此,需要對(duì)軸承振動(dòng)加速度信號(hào)做進(jìn)一步處理以實(shí)現(xiàn)故障特征提取。

    圖4 正常機(jī)組軸承數(shù)據(jù)降噪前后時(shí)域波形

    圖9為5種軸承狀態(tài)信號(hào)在不同尺度下的模糊熵。從圖中可以明顯看出正常間隙狀態(tài)信號(hào)的模糊熵隨著尺度的增大而呈現(xiàn)明顯的上升趨勢(shì),說明正常間隙狀態(tài)信號(hào)尺度越大,自相似程度越低,復(fù)雜度越高。

    而4種間隙大故障狀態(tài)的模糊熵隨著尺度的變化無趨勢(shì)上的明顯變化,在各個(gè)尺度上都比較接近,可分性差。因此,提出采用小波包模糊熵分析的方法。

    圖5 一級(jí)連桿小頭軸瓦間隙大故障數(shù)據(jù)降噪前后時(shí)域波形

    圖6 一級(jí)連桿大頭軸瓦間隙大故障數(shù)據(jù)降噪前后時(shí)域波形

    圖8 二級(jí)連桿大頭軸瓦間隙大故障數(shù)據(jù)降噪前后時(shí)域波形

    圖9 不同軸承狀態(tài)的多尺度模糊熵

    圖10為5種軸承狀態(tài)信號(hào)在n=3時(shí)小波包分解各節(jié)點(diǎn)模糊熵分布情況。

    圖10 不同軸承狀態(tài)的小波包模糊熵

    可以看出正常間隙狀態(tài)信號(hào)的模糊熵值總體趨勢(shì)是從低頻到高頻先下降后平穩(wěn),說明正常間隙狀態(tài)振動(dòng)信息主要分布在低頻部分。其余四種間隙大故障狀態(tài)模糊熵值隨頻率的變化無趨勢(shì)上的明顯變化,說明故障狀態(tài)的特征信息在高低頻都有分布,其次特征向量曲線具有無交叉、間隔大、可分性強(qiáng)特點(diǎn)。

    3.3 基于小波包模糊熵和支持向量機(jī)的故障診斷

    選擇適合小樣本分類的支持向量機(jī)(Support Vector Machine,SVM),由于不同位置軸承間隙故障分類屬于多類分類,本文采用一對(duì)一的多分類方法。對(duì)于SVM,核參數(shù)和誤差懲罰參數(shù)C對(duì)SVM的性能優(yōu)劣至關(guān)重要。本文核函數(shù)選用徑向基函數(shù),核函數(shù)參數(shù)為σ。利用LibSVM的遺傳算法對(duì)參數(shù)進(jìn)行優(yōu)化尋找最優(yōu)參數(shù)C和核函數(shù)參數(shù)σ,并建立支持向量機(jī)。

    利用小波包模糊熵分析考慮信號(hào)高頻成分的優(yōu)點(diǎn),采用小波包3層分解時(shí)各節(jié)點(diǎn)的模糊熵值作為故障特征向量。選取上述五種軸承間隙狀態(tài)的測(cè)試數(shù)據(jù)各120組,分別以各個(gè)狀態(tài)的80組特征向量作為訓(xùn)練樣本,40組特征向量作為測(cè)試樣本。把訓(xùn)練樣本數(shù)據(jù)輸入SVM,訓(xùn)練優(yōu)化結(jié)果為誤差懲罰參數(shù)C取2.14,核參數(shù)σ取3.21。對(duì)各個(gè)狀態(tài)測(cè)試樣本進(jìn)行識(shí)別,結(jié)果如表1所示。

    為了檢驗(yàn)小波包模糊熵方法的有效性,同樣選取上述5種軸承間隙狀態(tài)數(shù)據(jù)樣本各120組。首先,計(jì)算每組樣本中多尺度模糊熵的前8個(gè)尺度上的模糊熵,作為故障特征向量,分別以各個(gè)狀態(tài)的80組特征向量作為訓(xùn)練樣本,另外40組特征向量作為測(cè)試樣本,利用相同的方法訓(xùn)練SVM并進(jìn)行分類識(shí)別。多尺度模糊熵診斷的識(shí)別率如表1所示。

    由表1可知,兩種方法都可以實(shí)現(xiàn)往復(fù)壓縮機(jī)軸承不同位置故障的分類識(shí)別,但由于小波包分解具有自身將頻帶多層次劃分的優(yōu)勢(shì),基于小波包分解和模糊熵分析方法不但在總體識(shí)別率上高于多尺度模糊熵分析,而且各工況的識(shí)別率均為最高,識(shí)別效果最好,驗(yàn)證了該方法的優(yōu)越性。

    4 結(jié)語

    在模糊熵的基礎(chǔ)上,結(jié)合小波包可對(duì)頻帶多層次劃分的優(yōu)勢(shì)提出一種新的特征提取方法—小波包模糊熵方法。與多尺度模糊熵相比,小波包模糊熵同時(shí)分析了信號(hào)中低頻分量和高頻分量的故障信息,對(duì)故障特征描述更加全面;通過對(duì)往復(fù)壓縮機(jī)軸承振動(dòng)信號(hào)作3層小波包分解,取分解節(jié)點(diǎn)上的模糊熵值作為故障特征向量,輸入支持向量機(jī)進(jìn)行分類,識(shí)別效果更好,實(shí)現(xiàn)了故障的準(zhǔn)確診斷。

    表1 小波包模糊熵和多尺度模糊熵故障診斷準(zhǔn)確率比較

    猜你喜歡
    波包壓縮機(jī)尺度
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    BOG壓縮機(jī)在小型LNG船舶上的應(yīng)用
    淺析特大型往復(fù)壓縮機(jī)組的安裝
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    壓縮機(jī)組智能控制與節(jié)能增效
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    基于小波包變換的電力系統(tǒng)諧波分析
    小波包理論與圖像小波包分解
    9
    基于小波包的全信息解調(diào)方法及其應(yīng)用
    久久久久性生活片| 亚洲综合色惰| 国产真实伦视频高清在线观看| av免费观看日本| 看非洲黑人一级黄片| 观看av在线不卡| 亚洲美女视频黄频| 日韩精品有码人妻一区| 99久久精品热视频| 欧美xxⅹ黑人| 日本免费在线观看一区| 肉色欧美久久久久久久蜜桃| 亚洲国产av新网站| 亚洲国产最新在线播放| 久久精品国产a三级三级三级| 久久鲁丝午夜福利片| 久久久午夜欧美精品| 亚洲av福利一区| videossex国产| 国产爽快片一区二区三区| 国产毛片在线视频| 日日啪夜夜撸| 伦理电影免费视频| 亚洲精品成人av观看孕妇| 国产色婷婷99| 99精国产麻豆久久婷婷| 18禁裸乳无遮挡动漫免费视频| 在线观看国产h片| xxx大片免费视频| 日本欧美视频一区| 秋霞在线观看毛片| 国产精品伦人一区二区| 久久人人爽人人片av| 狂野欧美激情性bbbbbb| tube8黄色片| 精品少妇黑人巨大在线播放| 成人毛片a级毛片在线播放| 永久网站在线| 日日啪夜夜撸| 国产精品一及| 欧美国产精品一级二级三级 | 亚洲欧美精品自产自拍| 亚洲精品久久久久久婷婷小说| 欧美日韩在线观看h| 久久久久视频综合| 久久综合国产亚洲精品| 久久久久人妻精品一区果冻| 最近手机中文字幕大全| 国产无遮挡羞羞视频在线观看| 在线观看免费视频网站a站| 精品一区二区三卡| 黄色一级大片看看| 九九在线视频观看精品| 久久久久久九九精品二区国产| 国产亚洲5aaaaa淫片| 成人免费观看视频高清| 精品酒店卫生间| 亚洲欧美精品自产自拍| 成人黄色视频免费在线看| 国产精品免费大片| 一级毛片我不卡| 精品99又大又爽又粗少妇毛片| 久久这里有精品视频免费| 色5月婷婷丁香| 欧美精品人与动牲交sv欧美| 91aial.com中文字幕在线观看| 亚洲国产精品专区欧美| 久久ye,这里只有精品| 大片免费播放器 马上看| 少妇人妻精品综合一区二区| 亚洲人成网站在线播| av黄色大香蕉| 另类亚洲欧美激情| 免费观看av网站的网址| .国产精品久久| 国产一区二区三区av在线| 亚洲一级一片aⅴ在线观看| 日本猛色少妇xxxxx猛交久久| 各种免费的搞黄视频| 国产精品三级大全| 免费观看性生交大片5| 美女福利国产在线 | 老师上课跳d突然被开到最大视频| 精品一区二区三卡| 亚洲久久久国产精品| 欧美日韩一区二区视频在线观看视频在线| 交换朋友夫妻互换小说| 一区二区av电影网| 男女国产视频网站| 婷婷色综合www| 高清毛片免费看| 欧美变态另类bdsm刘玥| 国产免费又黄又爽又色| 午夜福利视频精品| 成年人午夜在线观看视频| 国产无遮挡羞羞视频在线观看| 久久久a久久爽久久v久久| 亚洲无线观看免费| 天堂中文最新版在线下载| 日本av免费视频播放| 国产成人精品婷婷| 91久久精品国产一区二区成人| 夜夜骑夜夜射夜夜干| 国产精品伦人一区二区| av视频免费观看在线观看| 午夜福利在线观看免费完整高清在| 国产男人的电影天堂91| 少妇熟女欧美另类| av天堂中文字幕网| 在线观看三级黄色| 免费观看无遮挡的男女| 国产成人aa在线观看| 嘟嘟电影网在线观看| 亚洲精品国产av蜜桃| 久久久精品免费免费高清| 边亲边吃奶的免费视频| 深爱激情五月婷婷| 国产欧美日韩一区二区三区在线 | 国产欧美日韩一区二区三区在线 | 国产欧美亚洲国产| 18禁裸乳无遮挡动漫免费视频| 国产精品熟女久久久久浪| h视频一区二区三区| 亚洲精品亚洲一区二区| 亚洲国产精品一区三区| av黄色大香蕉| 国产爱豆传媒在线观看| 亚洲aⅴ乱码一区二区在线播放| 直男gayav资源| 欧美三级亚洲精品| 欧美性感艳星| 久久精品国产a三级三级三级| 肉色欧美久久久久久久蜜桃| 色视频在线一区二区三区| 内地一区二区视频在线| 人妻夜夜爽99麻豆av| 男女国产视频网站| 午夜精品国产一区二区电影| 边亲边吃奶的免费视频| 校园人妻丝袜中文字幕| 日本欧美视频一区| 欧美亚洲 丝袜 人妻 在线| 午夜免费鲁丝| videos熟女内射| 国产女主播在线喷水免费视频网站| 亚洲精品国产成人久久av| 嘟嘟电影网在线观看| 成人亚洲欧美一区二区av| xxx大片免费视频| 久久久久性生活片| 久久久久久久久久久免费av| 一级片'在线观看视频| 建设人人有责人人尽责人人享有的 | 身体一侧抽搐| 多毛熟女@视频| 国产精品伦人一区二区| 永久免费av网站大全| 搡老乐熟女国产| 在线亚洲精品国产二区图片欧美 | 国产黄片美女视频| 成人影院久久| 亚洲精品456在线播放app| 一个人看的www免费观看视频| 各种免费的搞黄视频| 99热全是精品| 一个人看的www免费观看视频| 亚洲国产欧美在线一区| 婷婷色av中文字幕| 寂寞人妻少妇视频99o| 我要看日韩黄色一级片| videos熟女内射| 中文字幕av成人在线电影| 国产高清有码在线观看视频| 精品久久久噜噜| 在线观看免费日韩欧美大片 | 欧美日韩视频高清一区二区三区二| 免费黄网站久久成人精品| 成年美女黄网站色视频大全免费 | 国产亚洲欧美精品永久| 欧美日韩亚洲高清精品| 久久热精品热| 亚洲色图综合在线观看| 乱码一卡2卡4卡精品| 日韩av在线免费看完整版不卡| 亚洲欧美日韩卡通动漫| h视频一区二区三区| av免费观看日本| 中文字幕精品免费在线观看视频 | 午夜福利视频精品| 亚洲欧美成人精品一区二区| 国产欧美另类精品又又久久亚洲欧美| 美女视频免费永久观看网站| 交换朋友夫妻互换小说| 欧美成人a在线观看| 日韩电影二区| 少妇被粗大猛烈的视频| 在线免费十八禁| 日韩,欧美,国产一区二区三区| 国产一区亚洲一区在线观看| 内地一区二区视频在线| 亚洲精品视频女| 欧美一区二区亚洲| 欧美日韩国产mv在线观看视频 | 国语对白做爰xxxⅹ性视频网站| 亚洲国产最新在线播放| 老熟女久久久| 久久这里有精品视频免费| 久久久久久久久久久免费av| 一区在线观看完整版| 久久人人爽av亚洲精品天堂 | 男女边摸边吃奶| a级毛片免费高清观看在线播放| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 免费观看在线日韩| 日韩中文字幕视频在线看片 | 亚洲欧美精品自产自拍| 亚洲自偷自拍三级| 久久久久久久久久久丰满| 色网站视频免费| 一级毛片 在线播放| 亚洲av免费高清在线观看| 一级爰片在线观看| 永久免费av网站大全| 九九久久精品国产亚洲av麻豆| 丝袜脚勾引网站| 欧美高清成人免费视频www| 最新中文字幕久久久久| 51国产日韩欧美| 成人亚洲精品一区在线观看 | 美女高潮的动态| 成人免费观看视频高清| 激情 狠狠 欧美| 亚洲精品国产av蜜桃| 欧美精品一区二区大全| 卡戴珊不雅视频在线播放| 久热这里只有精品99| 一区二区av电影网| 特大巨黑吊av在线直播| 久久99蜜桃精品久久| 美女cb高潮喷水在线观看| 国产伦在线观看视频一区| 亚洲av男天堂| 三级国产精品片| 精品国产露脸久久av麻豆| 激情 狠狠 欧美| 国产亚洲av片在线观看秒播厂| 汤姆久久久久久久影院中文字幕| 亚洲av免费高清在线观看| 欧美日本视频| 国产精品人妻久久久影院| 嘟嘟电影网在线观看| 中文天堂在线官网| 性色avwww在线观看| 最近最新中文字幕大全电影3| 你懂的网址亚洲精品在线观看| 噜噜噜噜噜久久久久久91| 秋霞在线观看毛片| 国产熟女欧美一区二区| 国产成人免费无遮挡视频| 大又大粗又爽又黄少妇毛片口| 亚洲欧美日韩无卡精品| 亚洲精品成人av观看孕妇| 欧美三级亚洲精品| 一区二区三区乱码不卡18| 精品少妇黑人巨大在线播放| 亚洲国产精品一区三区| 国产精品国产av在线观看| 丰满人妻一区二区三区视频av| 极品教师在线视频| 亚洲精品久久久久久婷婷小说| av不卡在线播放| 亚洲色图av天堂| 丰满迷人的少妇在线观看| 大又大粗又爽又黄少妇毛片口| 99久久精品热视频| 成年av动漫网址| 免费人妻精品一区二区三区视频| 国产成人aa在线观看| 免费观看av网站的网址| 国语对白做爰xxxⅹ性视频网站| 国产在线视频一区二区| 视频中文字幕在线观看| 另类亚洲欧美激情| 成人特级av手机在线观看| 国产精品麻豆人妻色哟哟久久| 美女主播在线视频| 亚洲精品456在线播放app| 亚洲精品自拍成人| 精品一区二区三区视频在线| 夜夜看夜夜爽夜夜摸| 国产av国产精品国产| 人人妻人人看人人澡| 日韩中文字幕视频在线看片 | 亚洲自偷自拍三级| 免费黄频网站在线观看国产| 免费播放大片免费观看视频在线观看| 中文字幕亚洲精品专区| 欧美亚洲 丝袜 人妻 在线| 在线免费十八禁| 一级毛片黄色毛片免费观看视频| 欧美成人一区二区免费高清观看| 成人亚洲精品一区在线观看 | 蜜桃在线观看..| xxx大片免费视频| 97在线人人人人妻| 国产精品av视频在线免费观看| 超碰97精品在线观看| 男女边摸边吃奶| 精品一区二区三区视频在线| 成人毛片60女人毛片免费| 最近手机中文字幕大全| 肉色欧美久久久久久久蜜桃| 春色校园在线视频观看| 免费观看av网站的网址| av不卡在线播放| 欧美区成人在线视频| 91精品国产国语对白视频| 欧美人与善性xxx| 我要看日韩黄色一级片| 亚洲国产精品999| 黑人猛操日本美女一级片| 免费少妇av软件| 成人国产av品久久久| 简卡轻食公司| 高清黄色对白视频在线免费看 | 欧美一级a爱片免费观看看| 亚洲av男天堂| av在线app专区| 久久精品久久久久久噜噜老黄| 久久精品久久久久久久性| 亚洲精品乱码久久久久久按摩| 色婷婷久久久亚洲欧美| 日本黄大片高清| 一本久久精品| 亚洲电影在线观看av| 六月丁香七月| 人人妻人人澡人人爽人人夜夜| 天美传媒精品一区二区| 免费看光身美女| 免费大片18禁| 麻豆成人av视频| 日韩强制内射视频| av专区在线播放| 少妇的逼好多水| 日本欧美视频一区| 免费久久久久久久精品成人欧美视频 | 最近中文字幕2019免费版| 国产成人a∨麻豆精品| www.av在线官网国产| 免费观看性生交大片5| 亚洲在久久综合| 午夜免费鲁丝| 久久久久国产网址| 一二三四中文在线观看免费高清| 伦理电影大哥的女人| 91在线精品国自产拍蜜月| 久久精品国产亚洲av涩爱| 国产综合精华液| 麻豆成人av视频| 国产综合精华液| 亚洲,欧美,日韩| 青春草亚洲视频在线观看| 免费看不卡的av| 狂野欧美激情性xxxx在线观看| 免费看不卡的av| 亚洲精品国产成人久久av| 人人妻人人看人人澡| 亚洲精品国产成人久久av| 国产在线男女| 舔av片在线| 国产伦理片在线播放av一区| 亚洲精品久久午夜乱码| 三级经典国产精品| 最近手机中文字幕大全| 1000部很黄的大片| 熟女人妻精品中文字幕| 日韩一本色道免费dvd| 免费黄色在线免费观看| 日本与韩国留学比较| 热re99久久精品国产66热6| 日日啪夜夜撸| 国产精品伦人一区二区| 亚洲欧美成人综合另类久久久| av线在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站高清观看| 一边亲一边摸免费视频| 蜜桃在线观看..| 一级二级三级毛片免费看| 欧美日本视频| 啦啦啦中文免费视频观看日本| 尾随美女入室| av在线播放精品| 三级经典国产精品| 久久国内精品自在自线图片| 久久久久久九九精品二区国产| 狂野欧美白嫩少妇大欣赏| 日韩一本色道免费dvd| 秋霞在线观看毛片| 亚洲成人中文字幕在线播放| 内地一区二区视频在线| 一本色道久久久久久精品综合| 一个人免费看片子| 夫妻午夜视频| 亚洲性久久影院| 久久精品久久精品一区二区三区| 精品酒店卫生间| av国产久精品久网站免费入址| 高清午夜精品一区二区三区| 男人爽女人下面视频在线观看| 日韩一本色道免费dvd| 免费在线观看成人毛片| 观看av在线不卡| 丝瓜视频免费看黄片| 我的女老师完整版在线观看| 日韩电影二区| 美女xxoo啪啪120秒动态图| 国产成人freesex在线| 国产免费视频播放在线视频| 亚洲国产高清在线一区二区三| 99久国产av精品国产电影| 男的添女的下面高潮视频| 免费在线观看成人毛片| 少妇被粗大猛烈的视频| 在线观看人妻少妇| 精品久久国产蜜桃| 国产精品99久久久久久久久| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 黄色视频在线播放观看不卡| 国产一区二区三区av在线| 国产91av在线免费观看| 一本—道久久a久久精品蜜桃钙片| 成人综合一区亚洲| 舔av片在线| 日韩,欧美,国产一区二区三区| av卡一久久| 亚洲电影在线观看av| 国产人妻一区二区三区在| 天天躁夜夜躁狠狠久久av| 欧美丝袜亚洲另类| 亚洲欧美成人综合另类久久久| 免费黄频网站在线观看国产| a级毛片免费高清观看在线播放| 高清视频免费观看一区二区| 91精品伊人久久大香线蕉| 夜夜看夜夜爽夜夜摸| 欧美精品亚洲一区二区| 国产午夜精品久久久久久一区二区三区| 国产爽快片一区二区三区| 日韩大片免费观看网站| 另类亚洲欧美激情| 国产成人免费无遮挡视频| 一级毛片黄色毛片免费观看视频| 97精品久久久久久久久久精品| 日韩,欧美,国产一区二区三区| 国产精品久久久久久av不卡| 身体一侧抽搐| 国产69精品久久久久777片| 日韩伦理黄色片| 免费黄频网站在线观看国产| 亚洲精品国产成人久久av| 日本黄色日本黄色录像| 一本久久精品| 国产黄色免费在线视频| 久久久久久人妻| 亚洲精品乱久久久久久| 国产精品不卡视频一区二区| av在线播放精品| 全区人妻精品视频| 午夜激情福利司机影院| 又大又黄又爽视频免费| 亚洲人成网站高清观看| 99九九线精品视频在线观看视频| 久久女婷五月综合色啪小说| 久久久久网色| 久久国产乱子免费精品| 在线观看人妻少妇| 我要看日韩黄色一级片| 免费av不卡在线播放| 国产精品嫩草影院av在线观看| 午夜日本视频在线| 嘟嘟电影网在线观看| 99久久中文字幕三级久久日本| 亚洲av中文字字幕乱码综合| 啦啦啦视频在线资源免费观看| 99视频精品全部免费 在线| 最近2019中文字幕mv第一页| 少妇丰满av| 国产色婷婷99| 久久久久久久精品精品| 精品熟女少妇av免费看| 久久精品久久久久久久性| 国产成人91sexporn| 午夜免费观看性视频| av在线播放精品| 久久人妻熟女aⅴ| 99热这里只有是精品在线观看| 国产伦精品一区二区三区视频9| 国产色婷婷99| 久久99热这里只频精品6学生| 一级毛片黄色毛片免费观看视频| a级毛色黄片| 少妇人妻久久综合中文| 赤兔流量卡办理| 午夜福利视频精品| 天天躁夜夜躁狠狠久久av| 亚洲av不卡在线观看| 亚洲av福利一区| 在线精品无人区一区二区三 | 国产黄色免费在线视频| 国产在线一区二区三区精| 国国产精品蜜臀av免费| 亚洲在久久综合| 国产亚洲精品久久久com| 久久99蜜桃精品久久| videos熟女内射| 青春草亚洲视频在线观看| 黑丝袜美女国产一区| 久久精品夜色国产| 日韩制服骚丝袜av| 97超碰精品成人国产| 欧美日韩国产mv在线观看视频 | 国产久久久一区二区三区| 免费久久久久久久精品成人欧美视频 | 国产一区二区在线观看日韩| 联通29元200g的流量卡| 久久6这里有精品| 身体一侧抽搐| 亚洲国产成人一精品久久久| 久久人人爽人人爽人人片va| 国产成人精品福利久久| 乱系列少妇在线播放| 成人国产麻豆网| 国产在线男女| 亚洲国产最新在线播放| 国产黄频视频在线观看| 99热这里只有是精品50| 免费观看性生交大片5| 国产免费一级a男人的天堂| 女的被弄到高潮叫床怎么办| 2018国产大陆天天弄谢| 国产精品.久久久| 狂野欧美白嫩少妇大欣赏| 久久99热这里只有精品18| 18+在线观看网站| 97热精品久久久久久| 日日啪夜夜撸| 丝袜喷水一区| 免费观看性生交大片5| 永久免费av网站大全| 亚洲人成网站在线播| 欧美国产精品一级二级三级 | 丰满人妻一区二区三区视频av| 亚洲性久久影院| 亚洲av免费高清在线观看| 精品久久久久久久久亚洲| 久久这里有精品视频免费| 99久久综合免费| 亚洲av中文av极速乱| 老师上课跳d突然被开到最大视频| 国产精品女同一区二区软件| 亚洲欧美精品自产自拍| 中文字幕制服av| 夫妻午夜视频| 美女脱内裤让男人舔精品视频| 国产69精品久久久久777片| 在线观看三级黄色| 女性生殖器流出的白浆| 韩国av在线不卡| 97超视频在线观看视频| 久久鲁丝午夜福利片| 亚洲人与动物交配视频| 精品一区二区免费观看| 伦精品一区二区三区| 国产片特级美女逼逼视频| 一区二区av电影网| 国产毛片在线视频| 麻豆成人午夜福利视频| 一区二区av电影网| 色5月婷婷丁香| 国产老妇伦熟女老妇高清| 亚洲欧洲国产日韩| 熟女av电影| 亚洲av免费高清在线观看| 又大又黄又爽视频免费| 久久亚洲国产成人精品v| 久久久色成人| 草草在线视频免费看| 熟女av电影| 成年人午夜在线观看视频| 99久久精品热视频| 久久久久久久久久成人| 黄色视频在线播放观看不卡| 七月丁香在线播放| 亚洲久久久国产精品| 国产精品99久久99久久久不卡 | 黄片wwwwww| 蜜桃在线观看..| 色婷婷久久久亚洲欧美| 亚洲怡红院男人天堂| 国产精品不卡视频一区二区| 国产精品嫩草影院av在线观看| 少妇人妻 视频| 欧美+日韩+精品| 欧美精品亚洲一区二区| 丰满人妻一区二区三区视频av| 免费大片黄手机在线观看| 免费黄色在线免费观看| 国产精品人妻久久久影院| 亚洲经典国产精华液单| 亚洲av综合色区一区| 夫妻午夜视频| 成年美女黄网站色视频大全免费 |