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

    基于多尺度形態(tài)學(xué)分解譜熵的電機(jī)軸承性能退化特征提取

    2013-07-21 05:18:56王冰李洪儒許葆華
    軸承 2013年8期
    關(guān)鍵詞:形態(tài)學(xué)特征提取尺度

    王冰,李洪儒,許葆華

    (軍械工程學(xué)院,石家莊 050003)

    滾動(dòng)軸承是電動(dòng)機(jī)“轉(zhuǎn)子-軸承”系統(tǒng)中最為關(guān)鍵的機(jī)械零件,其性能往往會(huì)從完好逐漸經(jīng)歷一系列不同的退化狀態(tài)直至完全失效,軸承運(yùn)行狀態(tài)直接影響整臺(tái)機(jī)器的精度、可靠性及壽命[1]。當(dāng)前,針對(duì)電機(jī)軸承的研究重點(diǎn)一般集中于故障類型的診斷和故障位置的確定,相關(guān)技術(shù)方法也是以判斷不同故障類型的特征頻率為基礎(chǔ)。

    近年來,隨著維修理論及其相關(guān)技術(shù)的發(fā)展,基于狀態(tài)的維修越來越得到關(guān)注,故障預(yù)測(cè)技術(shù)則是實(shí)現(xiàn)基于狀態(tài)維修的核心。廣義上講,故障預(yù)測(cè)包括3個(gè)步驟:(1)早期故障檢測(cè)和預(yù)測(cè)特征提取;(2)退化狀態(tài)識(shí)別(估計(jì)當(dāng)前故障的嚴(yán)重程度);(3)預(yù)測(cè)故障發(fā)生時(shí)間或部件剩余使用壽命[2-3]。其中,預(yù)測(cè)特征提取和退化狀態(tài)識(shí)別是實(shí)現(xiàn)故障預(yù)測(cè)的基礎(chǔ)。

    預(yù)測(cè)特征提取又稱為性能退化特征提取,是設(shè)備退化狀態(tài)識(shí)別與故障預(yù)測(cè)的瓶頸,直接關(guān)系到退化狀態(tài)識(shí)別的準(zhǔn)確性和故障預(yù)測(cè)的可靠性[4]。傳統(tǒng)的特征提取技術(shù)側(cè)重于對(duì)不同故障類型的區(qū)分,而預(yù)測(cè)特征信息提取則側(cè)重于分析該信息能否反映設(shè)備從完好逐漸經(jīng)歷一系列不同退化狀態(tài)直至完全失效的性能退化過程[4-8]。文獻(xiàn)[4-5]提取小波相關(guān)特征尺度熵作為性能退化的描述信息;文獻(xiàn)[6-7]則將小波包分解的節(jié)點(diǎn)能量構(gòu)成特征向量,定量評(píng)估樣本的退化程度;文獻(xiàn)[8]則將通過循環(huán)平穩(wěn)分析得到的組合切片累積能量作為預(yù)測(cè)特征值,取得了較好的效果。

    數(shù)學(xué)形態(tài)分析是基于積分幾何和隨機(jī)集的不同于時(shí)、頻域分析的非線性方法。數(shù)學(xué)形態(tài)學(xué)在考察信號(hào)時(shí)使用結(jié)構(gòu)元素探針,通過結(jié)構(gòu)元素探針在信號(hào)中不斷移動(dòng)來提取有用信息從而進(jìn)行特征分析和描述[9-10]。多尺度形態(tài)學(xué)以數(shù)學(xué)形態(tài)分析為基礎(chǔ),直接從波形的幾何形態(tài)出發(fā),采用不同尺度的結(jié)構(gòu)元素對(duì)波形進(jìn)行變換,通過形態(tài)學(xué)膨脹、腐蝕和開、閉運(yùn)算,得到各尺度下的形態(tài)譜,進(jìn)而獲得對(duì)研究對(duì)象直觀的理解[11]。

    在此,分析了傳統(tǒng)無量綱時(shí)域特征參數(shù)在反映性能退化過程中的不足,并將多尺度形態(tài)學(xué)與信息熵理論相結(jié)合,提出一種基于多尺度形態(tài)分解譜熵的電機(jī)軸承性能退化特征提取方法。該方法以多尺度形態(tài)分解為基礎(chǔ),以分解后的譜熵(能譜熵和奇異譜熵)作為表征性能退化過程的特征向量,并通過仿真與試驗(yàn)數(shù)據(jù)對(duì)此方法的有效性進(jìn)行了驗(yàn)證。

    1 多尺度形態(tài)分解

    數(shù)學(xué)形態(tài)學(xué)的基本思想是利用具有一定形態(tài)的結(jié)構(gòu)元素去度量和提取信號(hào)中的對(duì)應(yīng)形態(tài),以達(dá)到對(duì)信號(hào)進(jìn)行分析和識(shí)別的目的。腐蝕和膨脹是數(shù)學(xué)形態(tài)學(xué)的基本運(yùn)算[12],分別等價(jià)于離散函數(shù)在滑動(dòng)濾波窗(相當(dāng)于結(jié)構(gòu)元素)內(nèi)的最小值和最大值濾波。

    在基于數(shù)學(xué)形態(tài)學(xué)的一維信號(hào)分析與處理中,結(jié)構(gòu)元素g(n)的選擇十分關(guān)鍵。單尺度形態(tài)學(xué)用固定類型和尺度的結(jié)構(gòu)元素度量和提取信號(hào)中的對(duì)應(yīng)形態(tài),以達(dá)到分析信號(hào)的目的。同理,如果能夠給出一個(gè)結(jié)構(gòu)元素序列,該序列中的所有結(jié)構(gòu)元素類型相同而尺度不同,運(yùn)用該序列對(duì)一維信號(hào)進(jìn)行分析,就可以用不同尺度的結(jié)構(gòu)元素提取信號(hào)中的對(duì)應(yīng)形態(tài),從而得到在不同尺度和“分辨率”下的一組處理結(jié)果。這種利用不同尺度的結(jié)構(gòu)元素進(jìn)行形態(tài)學(xué)變換的算法稱為多尺度形態(tài)學(xué)變換[13]。

    在多尺度形態(tài)學(xué)(Multiscale Morphology)中,多尺度開、閉是2種基本運(yùn)算,則對(duì)于原始信號(hào)f(n)有

    fng=fΘng⊕ng,

    (1)

    f?ng=f⊕ngΘng,

    (2)

    如果g為凸函數(shù),則ng可以通過基本結(jié)構(gòu)元素g對(duì)其自身進(jìn)行n-1次膨脹得到[14]

    (3)

    多尺度形態(tài)學(xué)分解建立在多尺度形態(tài)變換的基礎(chǔ)上。設(shè)f(x)為已知信號(hào),n為分解的層數(shù)。di(x)(1≤i≤n-2)為f(x)在第i尺度上的形態(tài)分解信號(hào),則有[15]

    (4)

    hi=0.5×[(fig?ig)(x)+(f?igig)·

    (x)]。

    可以看出,hi即為不同尺度結(jié)構(gòu)元素下的開-閉和閉-開組合形態(tài)濾波器[16]。而(4)式中的di(x)即為在i尺度上由濾波器hi所濾除的信號(hào)。在基于數(shù)學(xué)形態(tài)學(xué)的信號(hào)濾波中,其效果與所采用的結(jié)構(gòu)元素有著密切關(guān)系,只有與結(jié)構(gòu)元素的尺寸和形狀相匹配的信號(hào)基元才能被保留[17],因此,形態(tài)分解信號(hào)di(x)在小尺度下體現(xiàn)信號(hào)細(xì)節(jié),大尺度下體現(xiàn)信號(hào)輪廓,尺度與頻率相對(duì)應(yīng),從本質(zhì)上實(shí)現(xiàn)了信號(hào)在不同頻率分辨率下的多尺度劃分。

    多尺度形態(tài)學(xué)分解本質(zhì)是利用數(shù)學(xué)形態(tài)學(xué)變換,將復(fù)雜信號(hào)X分割成一系列互不相交的簡(jiǎn)單子集X1,X2,…,Xn,并對(duì)其求并集得到原始信號(hào)

    (5)

    由此可見,多尺度形態(tài)學(xué)分解的運(yùn)算結(jié)果是互相獨(dú)立的,運(yùn)算過程是完備的。

    2 多尺度形態(tài)分解譜熵的特征提取

    如何選取特征參數(shù)來表征電機(jī)軸承性能退化程度是進(jìn)行退化狀態(tài)識(shí)別的基礎(chǔ)。分析信號(hào)變化的本質(zhì),發(fā)現(xiàn)電機(jī)軸承性能退化的過程就是隨機(jī)成分不斷減小,特定頻率能量所占比例不斷增強(qiáng)的過程,因此,用信息熵理論定量刻畫性能退化過程是順其自然的。

    信息熵是從平均意義上表征信源總體信息測(cè)度的一個(gè)量,同時(shí)又是對(duì)信源輸出信息的不確定性和事件發(fā)生的隨機(jī)性的度量。當(dāng)信源中各個(gè)變量的概率分布越均勻,信息熵的值越大,信源包含的信息量也就越大;反之,當(dāng)信源中各個(gè)變量的概率分布越不均勻時(shí),信息熵的值就越小,信源包含的信息量也就越小[18]。從概率分布的角度來考慮各種譜分析,并運(yùn)用熵概念對(duì)譜分析結(jié)果進(jìn)行度量,多種譜熵的概念相繼提出,如小波能譜熵和Hilbert多尺度譜熵[19-20]等。

    多尺度形態(tài)分解將原信號(hào)在不同尺度和頻率下進(jìn)行劃分,同時(shí)也實(shí)現(xiàn)了信號(hào)能量在不同尺度上的分布,將該能量分布特征從概率角度分析,即可提出基于多尺度形態(tài)分解的譜熵的概念,即多尺度形態(tài)分解能譜熵和多尺度形態(tài)分解奇異譜熵。

    2.1 多尺度形態(tài)分解能譜熵

    對(duì)原信號(hào)進(jìn)行多尺度形態(tài)分解,可得到不同尺度下的信號(hào)分量di(i=1,2,…,n),計(jì)算各尺度下的能量值Ei,可得到多尺度形態(tài)分解能譜E=[E1,E2,…,En],形成了對(duì)信號(hào)在頻域上的一種劃分,則信號(hào)的多尺度形態(tài)分解能譜熵為

    (6)

    式中:pi=Ei/E,表示信號(hào)經(jīng)多尺度形態(tài)分解后,第i層的能量Ei占總能量E的比重。

    2.2 多尺度形態(tài)分解奇異譜熵

    將多尺度形態(tài)分解理論與奇異譜分析理論[21]相結(jié)合,可得到多尺度形態(tài)分解奇異譜熵。對(duì)原信號(hào)進(jìn)行多尺度形態(tài)分解,得到不同尺度下的信號(hào)分量di(i=1,2,…,n),將分量組成初始特征向量矩陣A,A=[d1,d2,…,dn]T。

    (7)

    多尺度形態(tài)分解奇異譜熵HMMQS反映了經(jīng)由多尺度形態(tài)分解后奇異成分能量分布的復(fù)雜程度。當(dāng)信號(hào)能量分布不均勻時(shí),奇異值分布結(jié)果分布于少數(shù)幾個(gè)分量,HMMQS就??;反之,當(dāng)信號(hào)能量分布更分散時(shí),奇異值分解結(jié)果越分散,HMMQS就越大。

    綜上所述,當(dāng)電機(jī)軸承出現(xiàn)故障,其振動(dòng)信號(hào)必然出現(xiàn)奇異現(xiàn)象,這樣,在其某一分量下的局部能量必然增大(或變小),相應(yīng)的概率也隨之變大(或變小),從而引起各分量能量分布的變化,從而必然導(dǎo)致HMMES和HMMQS值的變化,故障程度越深,其值越小。

    2.3 仿真驗(yàn)證

    為驗(yàn)證HMMES和HMMQS反映軸承性能退化過程的有效性,采用仿真信號(hào)[22]對(duì)其進(jìn)行分析。設(shè)仿真信號(hào)為

    x(t)=cos(2π×50t)+0.1t2cos(2π×10t+2),

    (8)

    其中,0.1t2cos(2π×10t+2)為故障模擬信號(hào),cos(2π×50t)為常規(guī)振動(dòng)信號(hào),幅值0.1t2則反映故障變化過程。 取信號(hào)采樣點(diǎn)數(shù)為N=10 240,采樣頻率為1 024 Hz,信號(hào)的時(shí)域波形如圖1所示。

    圖1 仿真信號(hào)時(shí)域圖形

    將信號(hào)等分為10段并順序標(biāo)記,用連續(xù)的10組數(shù)據(jù)描述仿真信號(hào)故障程度不斷加深的性能退化過程。為了說明HMMES和HMMQS對(duì)于描述性能退化過程的有效性,首先選取常見的無量綱時(shí)域特征參數(shù)對(duì)10組數(shù)據(jù)進(jìn)行對(duì)比分析,根據(jù)計(jì)算結(jié)果繪制的變化曲線如圖2所示。

    圖2 無量綱時(shí)域特征參數(shù)變化曲線

    由圖2可以看出,隨著故障程度的不斷加深,5條曲線均是先不斷增大,在第4組數(shù)據(jù)時(shí)達(dá)到峰值,之后則不斷減小。相比而言,波形指標(biāo)變化最為平緩。因此,傳統(tǒng)的無量綱時(shí)域特征參數(shù)的變化趨勢(shì)與性能退化的趨勢(shì)是不一致的,均不能作為表征其退化程度的特征參數(shù)。

    依次對(duì)每組數(shù)據(jù)進(jìn)行多尺度形態(tài)分解,分解層數(shù)n=3,基本結(jié)構(gòu)元素g=[0 0 0],以第1組數(shù)據(jù)為例,其分解結(jié)果如圖3所示。圖3中由上至下分別為原始信號(hào)及n=1,2,3時(shí)的分解信號(hào)。

    圖3 多尺度形態(tài)分解圖

    根據(jù)分解結(jié)果計(jì)算每組數(shù)據(jù)的HMMES和HMMQS值,并繪制變化曲線如圖4所示。

    圖4 HMMES與HMMQS變化曲線

    由圖4可以看出,2個(gè)參數(shù)隨故障程度的不斷加深而逐漸變小,與故障程度的變化趨勢(shì)一致,說明其對(duì)故障退化趨勢(shì)有很好的線性反映能力。另外,2條曲線在初期變化均比較平緩,這說明其對(duì)于初期微弱故障的反映能力較弱,整體而言,選用這2個(gè)參數(shù)來描述性能退化過程是可行的。

    3 實(shí)例分析

    為了驗(yàn)證多尺度形態(tài)分解譜熵的有效性,應(yīng)用某軸承試驗(yàn)數(shù)據(jù)[23]進(jìn)行實(shí)例分析。測(cè)試軸承為SKF6205-2RS深溝球軸承,采樣頻率12 kHz。用電火花機(jī)在軸承內(nèi)圈上人工加工局部損傷。損傷直徑為0.18 ,0.36 和0.54 mm,分別用于模擬軸承內(nèi)圈輕度、中度和重度損傷。試驗(yàn)數(shù)據(jù)在4種工況下取得,分別為工況1(轉(zhuǎn)速1 797 r/min,功率0 kW)、工況2(轉(zhuǎn)速1 772 r/min,功率0.75 kW)、工況3(轉(zhuǎn)速1 750 r/min,功率1.5 kW)、工況4(轉(zhuǎn)速1 730 r/min,功率2.25 kW)。

    在不同的工況下,分別選取一組正常、輕度、中度和重度故障軸承的狀態(tài)數(shù)據(jù)作為該工況下的標(biāo)準(zhǔn)退化試驗(yàn)數(shù)據(jù),數(shù)據(jù)長(zhǎng)度為6 000。

    以工況1下的標(biāo)準(zhǔn)退化試驗(yàn)數(shù)據(jù)為例進(jìn)行說明。該組數(shù)據(jù)包括4種不同狀態(tài)的軸承振動(dòng)信號(hào)。其時(shí)域波形如圖5所示。顯然,當(dāng)軸承內(nèi)圈處于正常狀態(tài)時(shí),信號(hào)呈現(xiàn)出幾乎隨機(jī)性的分布,不確定因素最高。當(dāng)內(nèi)圈出現(xiàn)故障后,信號(hào)的規(guī)律性逐漸增強(qiáng),幅值也不斷增大。

    圖5 不同狀態(tài)軸承時(shí)域波形

    利用多尺度形態(tài)分解方法分別對(duì)4種工況下的標(biāo)準(zhǔn)退化試驗(yàn)數(shù)據(jù)進(jìn)行處理,分解層數(shù)n=3,基本結(jié)構(gòu)元素g=[0 0 0]。計(jì)算得到不同損傷程度下的HMMES和HMMQS參數(shù)值,結(jié)果見表1。從表中數(shù)據(jù)分析可得:HMMES和HMMQS隨著故障程度的加深而呈現(xiàn)出下降趨勢(shì),且二者受工況差異的影響甚小,因此,用HMMES和HMMQS作為電機(jī)軸承性能退化特征是可行的。

    表1 4種工況下的多尺度形態(tài)分解譜熵

    4 結(jié)束語(yǔ)

    基于多尺度形態(tài)分解的能譜熵與奇異譜熵從不同角度描述了信號(hào)的復(fù)雜程度,且其變化趨勢(shì)與性能退化過程趨勢(shì)是一致的,因此,可以作為特征向量對(duì)電機(jī)軸承性能退化的過程進(jìn)行描述。

    所用的多尺度形態(tài)分解方法均采用扁平型結(jié)構(gòu)元素進(jìn)行處理。結(jié)構(gòu)元素類型、長(zhǎng)度對(duì)分解結(jié)果的影響均未作詳細(xì)分析。下一步可通過詳細(xì)對(duì)比分析,進(jìn)一步增強(qiáng)多尺度形態(tài)分解方法的科學(xué)性與指導(dǎo)性。

    猜你喜歡
    形態(tài)學(xué)特征提取尺度
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    基于Daubechies(dbN)的飛行器音頻特征提取
    電子制作(2018年19期)2018-11-14 02:37:08
    Bagging RCSP腦電特征提取算法
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    醫(yī)學(xué)微觀形態(tài)學(xué)在教學(xué)改革中的應(yīng)用分析
    基于MED和循環(huán)域解調(diào)的多故障特征提取
    9
    數(shù)學(xué)形態(tài)學(xué)濾波器在轉(zhuǎn)子失衡識(shí)別中的應(yīng)用
    巨核細(xì)胞數(shù)量及形態(tài)學(xué)改變?cè)谒姆N類型MPN中的診斷價(jià)值
    一May-Hegglin異常家系細(xì)胞形態(tài)學(xué)觀察
    久久精品亚洲熟妇少妇任你| 热re99久久国产66热| 久久中文字幕一级| 亚洲欧美精品综合一区二区三区| 男人的好看免费观看在线视频 | 国产xxxxx性猛交| 久久精品国产亚洲av香蕉五月| 三上悠亚av全集在线观看| 无限看片的www在线观看| 国产av一区二区精品久久| 女性被躁到高潮视频| 久久久久国产一级毛片高清牌| 这个男人来自地球电影免费观看| 欧美精品一区二区免费开放| 我的亚洲天堂| 97人妻天天添夜夜摸| 国产一区二区三区视频了| 伊人久久大香线蕉亚洲五| 亚洲专区字幕在线| 精品国产乱码久久久久久男人| 伊人久久大香线蕉亚洲五| 侵犯人妻中文字幕一二三四区| 在线观看免费高清a一片| 亚洲一区二区三区不卡视频| 久久青草综合色| 国产精品美女特级片免费视频播放器 | 欧美老熟妇乱子伦牲交| 色综合欧美亚洲国产小说| 俄罗斯特黄特色一大片| 俄罗斯特黄特色一大片| av天堂在线播放| 极品人妻少妇av视频| 久久国产精品影院| 中亚洲国语对白在线视频| 岛国视频午夜一区免费看| 天堂√8在线中文| 日韩欧美一区二区三区在线观看| av视频免费观看在线观看| 男女下面进入的视频免费午夜 | 午夜福利影视在线免费观看| 成人国语在线视频| 日韩高清综合在线| 在线观看午夜福利视频| 久久九九热精品免费| 久久香蕉激情| 久久精品成人免费网站| 免费av毛片视频| 黑人巨大精品欧美一区二区mp4| 精品国产亚洲在线| 亚洲欧美日韩另类电影网站| 悠悠久久av| 宅男免费午夜| 五月开心婷婷网| 国产一区二区三区在线臀色熟女 | 午夜免费激情av| 欧美色视频一区免费| 色婷婷av一区二区三区视频| 午夜老司机福利片| 久久精品91无色码中文字幕| 亚洲精品粉嫩美女一区| 午夜精品在线福利| 母亲3免费完整高清在线观看| 久久精品91无色码中文字幕| 免费看十八禁软件| 午夜福利在线观看吧| 久久亚洲真实| 中文字幕最新亚洲高清| 欧美日韩黄片免| 亚洲精品美女久久久久99蜜臀| 一级片免费观看大全| 亚洲色图 男人天堂 中文字幕| 大陆偷拍与自拍| avwww免费| 亚洲欧美精品综合一区二区三区| 又黄又爽又免费观看的视频| 美国免费a级毛片| 亚洲 欧美 日韩 在线 免费| 搡老岳熟女国产| 国产91精品成人一区二区三区| 欧美+亚洲+日韩+国产| 美女国产高潮福利片在线看| 99久久99久久久精品蜜桃| 国产精品乱码一区二三区的特点 | 精品无人区乱码1区二区| www.熟女人妻精品国产| 一级毛片精品| 精品国产乱子伦一区二区三区| bbb黄色大片| a在线观看视频网站| 一个人观看的视频www高清免费观看 | 90打野战视频偷拍视频| 国产精品野战在线观看 | 国产精品影院久久| 老司机靠b影院| 国产无遮挡羞羞视频在线观看| 天天影视国产精品| 日韩高清综合在线| 看黄色毛片网站| 99国产极品粉嫩在线观看| 人人妻人人爽人人添夜夜欢视频| 长腿黑丝高跟| 国产蜜桃级精品一区二区三区| 精品电影一区二区在线| 757午夜福利合集在线观看| av免费在线观看网站| 777久久人妻少妇嫩草av网站| 国产精品99久久99久久久不卡| 久久久久久久精品吃奶| av网站免费在线观看视频| 首页视频小说图片口味搜索| 少妇 在线观看| 18禁国产床啪视频网站| 99精品久久久久人妻精品| 在线观看免费高清a一片| 亚洲精品在线美女| 精品人妻1区二区| 波多野结衣av一区二区av| 亚洲 欧美 日韩 在线 免费| 在线观看免费视频日本深夜| 成人特级黄色片久久久久久久| 精品一区二区三区视频在线观看免费 | 精品国产乱子伦一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 一级作爱视频免费观看| 久久天躁狠狠躁夜夜2o2o| 国产亚洲欧美在线一区二区| 国产成人精品在线电影| 两人在一起打扑克的视频| 色综合婷婷激情| 亚洲精品国产区一区二| 天天躁夜夜躁狠狠躁躁| 新久久久久国产一级毛片| 一边摸一边抽搐一进一出视频| 欧美日韩国产mv在线观看视频| 99国产精品99久久久久| 久久青草综合色| 啦啦啦在线免费观看视频4| 欧美日韩瑟瑟在线播放| 少妇被粗大的猛进出69影院| 69av精品久久久久久| 国产色视频综合| 国产精品av久久久久免费| 狠狠狠狠99中文字幕| 桃红色精品国产亚洲av| 国产极品粉嫩免费观看在线| 欧美日韩一级在线毛片| e午夜精品久久久久久久| 18禁黄网站禁片午夜丰满| 夜夜夜夜夜久久久久| 天堂影院成人在线观看| 婷婷精品国产亚洲av在线| 久久精品国产亚洲av高清一级| 久久久久精品国产欧美久久久| 国产精品爽爽va在线观看网站 | 最好的美女福利视频网| 在线天堂中文资源库| 女生性感内裤真人,穿戴方法视频| 女生性感内裤真人,穿戴方法视频| 宅男免费午夜| 在线天堂中文资源库| 天堂中文最新版在线下载| 亚洲国产精品一区二区三区在线| 国产成人影院久久av| 国内久久婷婷六月综合欲色啪| 在线观看www视频免费| 亚洲av日韩精品久久久久久密| 大码成人一级视频| av免费在线观看网站| 国产成+人综合+亚洲专区| 日日夜夜操网爽| 欧美成人午夜精品| 亚洲成人精品中文字幕电影 | 日韩免费av在线播放| 老司机在亚洲福利影院| 一级a爱视频在线免费观看| 亚洲精品国产区一区二| 美女午夜性视频免费| 男女之事视频高清在线观看| 精品第一国产精品| 亚洲熟女毛片儿| 波多野结衣av一区二区av| 国产av在哪里看| 黑人巨大精品欧美一区二区mp4| 午夜精品国产一区二区电影| 两个人免费观看高清视频| 中文字幕另类日韩欧美亚洲嫩草| 视频在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| xxxhd国产人妻xxx| 亚洲国产毛片av蜜桃av| 久久久久九九精品影院| 不卡av一区二区三区| 精品国产一区二区三区四区第35| 丝袜美腿诱惑在线| 一区二区三区精品91| 亚洲黑人精品在线| 精品免费久久久久久久清纯| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利在线观看吧| 可以免费在线观看a视频的电影网站| 欧美成人午夜精品| 亚洲精品中文字幕一二三四区| 国产精品一区二区免费欧美| 自线自在国产av| 黄片小视频在线播放| 少妇裸体淫交视频免费看高清 | av中文乱码字幕在线| 国产一区二区三区视频了| 亚洲国产精品sss在线观看 | 亚洲午夜理论影院| 精品少妇一区二区三区视频日本电影| 一边摸一边做爽爽视频免费| 日韩中文字幕欧美一区二区| 亚洲一区二区三区色噜噜 | 久久亚洲精品不卡| 国产精品一区二区三区四区久久 | 亚洲欧美日韩高清在线视频| 多毛熟女@视频| 国产欧美日韩一区二区精品| 欧美成人性av电影在线观看| av在线天堂中文字幕 | x7x7x7水蜜桃| 亚洲午夜理论影院| 人人妻人人爽人人添夜夜欢视频| 在线观看午夜福利视频| 国内久久婷婷六月综合欲色啪| 久久久久久亚洲精品国产蜜桃av| 很黄的视频免费| 久久伊人香网站| 日日摸夜夜添夜夜添小说| 国产乱人伦免费视频| 日本欧美视频一区| 久久中文字幕一级| 男人操女人黄网站| 视频区图区小说| av网站在线播放免费| √禁漫天堂资源中文www| 色综合婷婷激情| 欧美激情久久久久久爽电影 | 亚洲三区欧美一区| 亚洲av熟女| 国产成人av激情在线播放| 国产亚洲精品第一综合不卡| 69av精品久久久久久| 人妻久久中文字幕网| xxx96com| 岛国在线观看网站| 欧美激情久久久久久爽电影 | 欧美人与性动交α欧美软件| 久热爱精品视频在线9| ponron亚洲| 久热这里只有精品99| 淫秽高清视频在线观看| 日本wwww免费看| 国产黄a三级三级三级人| 精品久久久久久久毛片微露脸| 99国产综合亚洲精品| 亚洲黑人精品在线| 国产精品偷伦视频观看了| 最新在线观看一区二区三区| 亚洲精品久久午夜乱码| 午夜两性在线视频| 黄色女人牲交| 亚洲欧美一区二区三区黑人| 一区在线观看完整版| 首页视频小说图片口味搜索| 久久人人爽av亚洲精品天堂| 亚洲三区欧美一区| 久久国产乱子伦精品免费另类| 亚洲av日韩精品久久久久久密| 日本撒尿小便嘘嘘汇集6| 一本综合久久免费| 免费在线观看完整版高清| 亚洲人成电影观看| 久久精品国产99精品国产亚洲性色 | 亚洲色图av天堂| 久久精品91蜜桃| 99精品欧美一区二区三区四区| 男女下面插进去视频免费观看| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 国产精品自产拍在线观看55亚洲| 可以在线观看毛片的网站| 亚洲国产看品久久| 欧美中文日本在线观看视频| 亚洲三区欧美一区| 啦啦啦 在线观看视频| 亚洲人成电影免费在线| 国产激情久久老熟女| 91老司机精品| 岛国在线观看网站| 村上凉子中文字幕在线| 亚洲av成人av| 国产精品成人在线| 正在播放国产对白刺激| 在线视频色国产色| 99香蕉大伊视频| 黑人操中国人逼视频| 嫩草影视91久久| 国产精品久久久av美女十八| 美女高潮喷水抽搐中文字幕| 麻豆av在线久日| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文日韩欧美视频| 精品人妻1区二区| 欧美成人性av电影在线观看| 美国免费a级毛片| 长腿黑丝高跟| 久久久久国产一级毛片高清牌| 老熟妇乱子伦视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 1024香蕉在线观看| 欧美激情高清一区二区三区| 久久这里只有精品19| 久久久久久久久中文| 久久久久久人人人人人| 国产午夜精品久久久久久| 极品人妻少妇av视频| 亚洲成人久久性| 黄片大片在线免费观看| 色综合婷婷激情| 国产又色又爽无遮挡免费看| 国产亚洲精品久久久久5区| 久久精品亚洲av国产电影网| 欧美中文综合在线视频| 一本综合久久免费| 免费少妇av软件| 国产精品偷伦视频观看了| 50天的宝宝边吃奶边哭怎么回事| 18禁观看日本| 超色免费av| 国产精品免费视频内射| 满18在线观看网站| 老司机深夜福利视频在线观看| 国产免费av片在线观看野外av| 亚洲五月婷婷丁香| 午夜福利影视在线免费观看| 一区二区日韩欧美中文字幕| 中国美女看黄片| xxx96com| 精品国产超薄肉色丝袜足j| 亚洲av美国av| 欧美日韩一级在线毛片| 91精品国产国语对白视频| 午夜影院日韩av| 视频在线观看一区二区三区| 窝窝影院91人妻| 久久狼人影院| 波多野结衣一区麻豆| 免费av中文字幕在线| 成人永久免费在线观看视频| 男女午夜视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 久久久国产成人精品二区 | 婷婷丁香在线五月| 亚洲第一av免费看| av网站在线播放免费| 啦啦啦免费观看视频1| 亚洲人成伊人成综合网2020| 黑人巨大精品欧美一区二区mp4| 久久精品影院6| 成人黄色视频免费在线看| xxxhd国产人妻xxx| 97碰自拍视频| 99国产综合亚洲精品| 欧美另类亚洲清纯唯美| 色老头精品视频在线观看| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美98| 黑人巨大精品欧美一区二区蜜桃| 国产成人影院久久av| 欧美乱码精品一区二区三区| 久久精品国产亚洲av香蕉五月| 国产精华一区二区三区| 久久热在线av| 久久人人97超碰香蕉20202| 日本黄色日本黄色录像| 欧美精品一区二区免费开放| 黑人欧美特级aaaaaa片| 欧美精品亚洲一区二区| 久久精品国产清高在天天线| 69精品国产乱码久久久| 99久久99久久久精品蜜桃| 人人妻,人人澡人人爽秒播| 桃红色精品国产亚洲av| 欧美日韩黄片免| 欧美色视频一区免费| 欧美日韩亚洲高清精品| 岛国视频午夜一区免费看| 久久影院123| 亚洲男人天堂网一区| 黄色a级毛片大全视频| 日本黄色日本黄色录像| 色综合婷婷激情| 人人妻人人爽人人添夜夜欢视频| 亚洲色图综合在线观看| 嫩草影院精品99| 欧美最黄视频在线播放免费 | 亚洲精品在线观看二区| 精品一区二区三卡| 久久国产精品男人的天堂亚洲| 成人手机av| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 国产激情欧美一区二区| 精品国产美女av久久久久小说| 日本撒尿小便嘘嘘汇集6| 老熟妇乱子伦视频在线观看| 久久中文看片网| 欧美人与性动交α欧美精品济南到| 午夜91福利影院| 久久精品aⅴ一区二区三区四区| 欧美精品一区二区免费开放| 久久精品亚洲熟妇少妇任你| 国产欧美日韩精品亚洲av| 大码成人一级视频| 久久欧美精品欧美久久欧美| 欧美中文日本在线观看视频| 男女做爰动态图高潮gif福利片 | 久久精品91无色码中文字幕| 亚洲专区字幕在线| 在线免费观看的www视频| 精品无人区乱码1区二区| 国产视频一区二区在线看| 午夜福利在线免费观看网站| 桃色一区二区三区在线观看| 淫妇啪啪啪对白视频| 999久久久国产精品视频| 欧美成人性av电影在线观看| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女 | 精品久久久久久,| 高清黄色对白视频在线免费看| 可以在线观看毛片的网站| 国产主播在线观看一区二区| 女警被强在线播放| 久久精品91蜜桃| 精品一品国产午夜福利视频| 一二三四社区在线视频社区8| 这个男人来自地球电影免费观看| 成人三级做爰电影| 亚洲国产欧美日韩在线播放| 日本a在线网址| 麻豆一二三区av精品| 91麻豆精品激情在线观看国产 | 免费在线观看视频国产中文字幕亚洲| 国产精品久久久久久人妻精品电影| 久久香蕉精品热| 国产欧美日韩一区二区精品| 91成人精品电影| 校园春色视频在线观看| 美女高潮喷水抽搐中文字幕| 精品国内亚洲2022精品成人| 亚洲精品粉嫩美女一区| 超碰成人久久| 超色免费av| 男人舔女人的私密视频| 99久久综合精品五月天人人| 91麻豆精品激情在线观看国产 | 性色av乱码一区二区三区2| 亚洲精品一卡2卡三卡4卡5卡| 如日韩欧美国产精品一区二区三区| 97碰自拍视频| 精品久久久久久成人av| 免费高清视频大片| 无限看片的www在线观看| av超薄肉色丝袜交足视频| 国产99白浆流出| 极品人妻少妇av视频| 国产不卡一卡二| 亚洲情色 制服丝袜| 午夜福利免费观看在线| 欧美日韩av久久| 欧美成狂野欧美在线观看| 99国产极品粉嫩在线观看| 熟女少妇亚洲综合色aaa.| 亚洲av美国av| 精品久久久精品久久久| 一边摸一边抽搐一进一小说| 三上悠亚av全集在线观看| 国产单亲对白刺激| 国产精品免费一区二区三区在线| 美女福利国产在线| 日韩av在线大香蕉| 超碰97精品在线观看| aaaaa片日本免费| 国产主播在线观看一区二区| 国产激情欧美一区二区| 亚洲国产精品一区二区三区在线| 国产99白浆流出| 夜夜夜夜夜久久久久| 日韩欧美一区二区三区在线观看| 一级毛片高清免费大全| 夜夜爽天天搞| 天天躁夜夜躁狠狠躁躁| 国产精品久久电影中文字幕| 亚洲精品粉嫩美女一区| 岛国在线观看网站| 久久久国产成人免费| 亚洲国产毛片av蜜桃av| 久久久久久久久免费视频了| 久久中文字幕一级| 日韩中文字幕欧美一区二区| 夜夜夜夜夜久久久久| 久久天堂一区二区三区四区| 99久久99久久久精品蜜桃| 最近最新免费中文字幕在线| 自线自在国产av| www.精华液| 国产一区二区激情短视频| 女性被躁到高潮视频| 欧美日韩视频精品一区| 丁香欧美五月| 免费观看精品视频网站| 亚洲av熟女| 亚洲成人久久性| 91麻豆av在线| 亚洲av成人av| 日韩精品中文字幕看吧| 在线观看日韩欧美| 精品国产一区二区三区四区第35| av国产精品久久久久影院| 欧美成人性av电影在线观看| 成人av一区二区三区在线看| 高清av免费在线| 19禁男女啪啪无遮挡网站| www日本在线高清视频| 男人操女人黄网站| 巨乳人妻的诱惑在线观看| 午夜免费鲁丝| 亚洲av熟女| 女人被躁到高潮嗷嗷叫费观| 大码成人一级视频| 黄色视频不卡| 777久久人妻少妇嫩草av网站| 18禁黄网站禁片午夜丰满| 久久精品aⅴ一区二区三区四区| 免费一级毛片在线播放高清视频 | 精品欧美一区二区三区在线| 亚洲av成人不卡在线观看播放网| 亚洲成a人片在线一区二区| www日本在线高清视频| 99国产极品粉嫩在线观看| 午夜精品在线福利| 亚洲一区二区三区色噜噜 | 亚洲avbb在线观看| 国产成人系列免费观看| 国产单亲对白刺激| 午夜福利欧美成人| 黄色丝袜av网址大全| 欧美黄色片欧美黄色片| 精品卡一卡二卡四卡免费| 在线观看午夜福利视频| 欧美日本亚洲视频在线播放| 亚洲情色 制服丝袜| 亚洲av第一区精品v没综合| 欧美日韩av久久| 中国美女看黄片| 亚洲色图综合在线观看| 老熟妇仑乱视频hdxx| 成人黄色视频免费在线看| 亚洲专区中文字幕在线| 精品国产一区二区久久| 神马国产精品三级电影在线观看 | 一区二区三区激情视频| 中文欧美无线码| 天堂俺去俺来也www色官网| 久久久国产成人免费| 久久欧美精品欧美久久欧美| 女人被狂操c到高潮| 午夜a级毛片| 高清在线国产一区| 精品乱码久久久久久99久播| 亚洲欧美日韩高清在线视频| 五月开心婷婷网| 久久精品国产综合久久久| 中文字幕另类日韩欧美亚洲嫩草| 国产精品国产高清国产av| 免费看十八禁软件| 国产男靠女视频免费网站| 精品一区二区三区av网在线观看| 男女下面插进去视频免费观看| 嫩草影视91久久| 丰满迷人的少妇在线观看| 亚洲七黄色美女视频| 91字幕亚洲| 国产日韩一区二区三区精品不卡| а√天堂www在线а√下载| 亚洲全国av大片| 制服人妻中文乱码| 日日干狠狠操夜夜爽| 成人国产一区最新在线观看| 久久人妻av系列| 国产成人精品无人区| 日韩欧美国产一区二区入口| 成人18禁高潮啪啪吃奶动态图| 狠狠狠狠99中文字幕| 欧美+亚洲+日韩+国产| 久久热在线av| 69av精品久久久久久| av中文乱码字幕在线| 一级毛片女人18水好多| 久久精品91无色码中文字幕| 在线视频色国产色| 国产1区2区3区精品| 视频在线观看一区二区三区| 国产成人系列免费观看| 在线观看免费高清a一片| 男人的好看免费观看在线视频 | 一级毛片女人18水好多| 色精品久久人妻99蜜桃| 久久久国产欧美日韩av| 看免费av毛片| 国产黄色免费在线视频| 国产真人三级小视频在线观看|