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

    基于D-CA和R-EEMD的液壓系統(tǒng)故障識別

    2015-12-28 06:40:24張梅軍
    噪聲與振動控制 2015年1期
    關(guān)鍵詞:端點(diǎn)頻域時(shí)域

    柴 凱,張梅軍,黃 杰,馮 霞

    (解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院,南京 210007)

    基于D-CA和R-EEMD的液壓系統(tǒng)故障識別

    柴 凱,張梅軍,黃 杰,馮 霞

    (解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院,南京 210007)

    剖析液壓系統(tǒng)故障特征,采用了一種雙相關(guān)分析(D-CA)和改進(jìn)的集合經(jīng)驗(yàn)?zāi)B(tài)分解(R-EEMD)相結(jié)合的液壓系統(tǒng)故障提取新方法。該方法首先對原信號進(jìn)行自相關(guān)分析,突出信號中的周期成分和去噪,利用支持向量回歸機(jī)(SVR)延拓來改進(jìn)的EEMD對原信號的自相關(guān)函數(shù)進(jìn)行分解;得到理論意義上的固有模態(tài)函數(shù)(IMF)。再通過求取IMF分量與自相關(guān)處理的信號的頻域而非傳統(tǒng)時(shí)域上的互相關(guān)系數(shù),去除虛假IMF分量。最后對去除虛假分量重構(gòu)信號進(jìn)行Hilbert譜分析提取信號的故障特征。該方法提高了信噪比,減少了IMF的數(shù)量,抑制了端點(diǎn)效應(yīng),成功地提取了液壓系統(tǒng)故障特征頻率。

    振動與波;故障診斷;雙相關(guān)分析;集合經(jīng)驗(yàn)?zāi)B(tài)分解

    液壓系統(tǒng)故障信號由于測試儀器、測試環(huán)境和人為因素等影響,常被淹沒在強(qiáng)大的背景噪聲中,采集的信號是典型的非線性、非平穩(wěn)信號,很難被發(fā)現(xiàn)和提取[1]。EMD(empirical mode decomposition)方法從信號本身尺度特征出發(fā)對信號進(jìn)行分解,具有良好的自適應(yīng)性,能夠得到信號的時(shí)間—頻率—幅度分布特征[2]。但是EMD分解中存在IMF的虛假分量、端點(diǎn)效應(yīng)、模態(tài)混疊等缺陷影響其實(shí)際中的應(yīng)用[3]。為減小EMD分解過程中的頻譜混疊現(xiàn)象,Wu與Huang等提出了EEMD(Ensemble Empirical Mode Decomposition),但是EEMD并未有效地解決EMD的虛假分量和端點(diǎn)效應(yīng)[4]。文獻(xiàn)[5,6]分別利用相似波形延拓、RBF神經(jīng)網(wǎng)絡(luò)來抑制端點(diǎn)效應(yīng),但延拓結(jié)果很難保證分解的統(tǒng)一,難以實(shí)現(xiàn)信號較準(zhǔn)確地預(yù)測,不可避免的引入誤差。目前利用相關(guān)分析與EMD(EEMD)結(jié)合主要用于以下兩個(gè)方面:一是通過時(shí)域互相關(guān)系數(shù)剔除分解的虛假模態(tài)分量[7,8],但相關(guān)程度大的分量較多,相關(guān)系數(shù)數(shù)值差別較小,難以確定閾值;二是將EMD與自相關(guān)分析結(jié)合,通過分解后的各分量作自相關(guān)分析,判別每層IMF是否含有微弱周期信號[9],但是自相關(guān)次數(shù)較多。在時(shí)域內(nèi)噪聲對信號幅值影響較大,從而使EEMD分解的互相關(guān)主分量增多;而在頻域內(nèi),信號的功率譜能有效突出信號的峰值頻率,受噪聲影響小。因此,本文提出了基于雙相關(guān)分析和SVR延拓改進(jìn)的EEMD故障特征提取方法,在EEMD分解前先對原信號求自相關(guān)函數(shù),對原信號的自相關(guān)函數(shù)進(jìn)行EEMD分解,利用自相關(guān)突出信號的周期成分和去噪。利用SVR(Support Vector Regression)來保證抑制端點(diǎn)效應(yīng)使包絡(luò)線失真度達(dá)到最小。利用頻域互相關(guān)系數(shù)消除虛假IMF分量,并將該方法引入液壓系統(tǒng)故障特征提取中,有效解決了影響EMD應(yīng)用中的主要問題,取得了良好的效果。

    1 雙相關(guān)分析和改進(jìn)的EEMD基本理論

    雙相關(guān)分析是指對原始信號的自相關(guān)函數(shù)和兩個(gè)信號頻域內(nèi)的互相關(guān)系數(shù)。利用信號時(shí)域內(nèi)的自相關(guān)函數(shù)不僅可以去噪,而且能夠突出故障信號的周期分量,去除不相關(guān)的非周期趨勢。借助SVR延拓來改進(jìn)EEMD,抑制了端點(diǎn)效應(yīng)。經(jīng)IMF與自相關(guān)處理后信號的頻域互相關(guān)系數(shù)確定信號中代表故障特征信息的主要IMF分量,剔除虛假分量。最后根據(jù)主要分量重構(gòu)信號,進(jìn)行Hilbert譜分析來提取相應(yīng)的故障特征頻率。

    1.1 自相關(guān)函數(shù)去噪

    自相關(guān)函數(shù)定義為動態(tài)信號在某一時(shí)刻的數(shù)據(jù)x(t)與另一時(shí)刻的數(shù)據(jù)之間的相互依賴關(guān)系。

    同理,互相關(guān)函數(shù)是反映兩個(gè)動態(tài)信號x(t)和之間的相互依賴關(guān)系[10]。

    其中τ為延遲時(shí)間,T為采樣周期。

    當(dāng)信號x(t)被噪聲信號n(t)污染后變?yōu)閟(t),則

    式(3)中,當(dāng)τ不為0且較大時(shí),Rxn(τ),Rnx(τ)及Rn(τ)均趨于0,所以染噪信號的自相關(guān)函數(shù)具有周期性,且周期與原信號相同,因此噪聲經(jīng)自相關(guān)后得到消除和抑制[11]。

    1.2 SVR延拓改進(jìn)的EEMD

    EEMD是在EMD“篩分”的基礎(chǔ)上,通過在原數(shù)據(jù)中加入適當(dāng)大小的白噪聲,給EMD提供一個(gè)相對一致的參照尺度分布,經(jīng)多次計(jì)算后求取平均值使所加的白噪聲互相抵消,最后得到尺度一致的分量,一定程度上解決了模態(tài)混疊的問題[12]。本文通過SVR延拓來改進(jìn)EEMD,抑制分解過程中存在的端點(diǎn)效應(yīng)。具體步驟是:

    (1)對自相關(guān)處理后的信號Rx(t)進(jìn)行SVR延拓;

    (2)在延拓后的信號中加入一組均值為零、方差相等的隨機(jī)白噪聲序列;

    (3)對延拓和加噪處理后的復(fù)合信號Rx1(t)進(jìn)行EMD分解;

    (4)復(fù)重m次1~3過程,得到m組經(jīng)過SVR延拓和加噪預(yù)處理的復(fù)合信號,經(jīng)EMD分解后得到m組IMF分量和m組余量;

    (5)求出m組IMF分量相應(yīng)的均值Ci和m組余量的均值R;

    (6)截?cái)郈i和R中SVR延拓部分,作為信號Rx(t)改進(jìn)的EEMD分解結(jié)果。

    1.3 頻域互相關(guān)系數(shù)去除虛假分量

    與傳統(tǒng)的時(shí)域相關(guān)系數(shù)去除偽分量不同,本文采用改進(jìn)的EEMD分解出的各個(gè)IMF分量功率譜與自相關(guān)處理后信號功率譜的互相關(guān)系數(shù)大小來去除虛假IMF分量。這是考慮在時(shí)域內(nèi)噪聲對信號影響較大,從而影響時(shí)域相關(guān)系數(shù),而在頻域內(nèi),噪聲雖然覆蓋整個(gè)頻段,但是噪聲功率譜密度小,信號的功率譜頻率峰值突出,信號在頻域內(nèi)功率譜的互相關(guān)系數(shù)受噪聲干擾小。兩個(gè)時(shí)間序列xi,yi的相關(guān)系數(shù)ρxy定義如下

    推廣到頻域內(nèi)的互相關(guān)系數(shù),xi,yi分別用對應(yīng)的功率譜Gx,Gy替代,fa為分析頻率,即

    1.4 D-CA和R-EEMD的算法流程

    D-CA和R-EEMD的算法流程如圖1。

    2 仿真信號分析

    采用含有三個(gè)頻率成分:5 Hz,10 Hz和30 Hz的仿真信號進(jìn)行實(shí)驗(yàn),采樣頻率為400 Hz,采樣點(diǎn)數(shù)為800,原始信號

    加入高斯白噪聲,時(shí)域波形如圖2—圖3,圖4是相應(yīng)的功率譜。對比發(fā)現(xiàn),在時(shí)域內(nèi)噪聲對信號影響較大,而在頻域內(nèi)雖然噪聲覆蓋整個(gè)頻段,但是并未對信號的功率譜主要成分產(chǎn)生明顯污染。

    圖1 D-CA和R-EEMD算法流程圖

    圖2 原始信號時(shí)域波形

    圖3 染噪信號時(shí)域波形

    圖4 仿真信號的功率譜

    2.1 仿真信號的EEMD分析

    對染噪信號直接進(jìn)行EEMD分析得到8個(gè)IMF分量和一個(gè)余量(圖5),并對主要IMF進(jìn)行Hilbert變換得到時(shí)頻譜(圖6)。從圖6可知原始信號的三個(gè)頻率成分淹沒在噪聲中,EEMD分解未能將信號的特征提取出來,頻譜混疊嚴(yán)重,主要的IMF分量無法表達(dá)出應(yīng)有的物理意義,而且產(chǎn)生較多的虛假分量,存在明顯的端點(diǎn)效應(yīng)。

    圖5 染噪信號直接EEMD分解圖

    圖6 染噪信號EEMD分解的Hilbert譜

    2.2 仿真信號D-CA和R-EEMD分析

    運(yùn)用D-CA和R-EEMD算法對仿真信號進(jìn)行處理,信號自相關(guān)處理后,進(jìn)行改進(jìn)的EEMD分解,SVR延拓選擇線性核函數(shù),左右延拓100個(gè)數(shù)據(jù)點(diǎn),分解后截?cái)嘌油夭糠?,將得到的IMF分量與自相關(guān)處理的信號進(jìn)行頻域互相關(guān),分析結(jié)果如表1。相關(guān)程度較大三個(gè)IMF分量:IMF1,IMF2,IMF3。比較表2的時(shí)域互相關(guān)系數(shù),有五個(gè)IMF分量相關(guān)程度較大:IMF1-IMF5。

    頻域互相關(guān)和時(shí)域相關(guān)得到的主要IMF分量不一致的原因是噪聲對時(shí)域信號影響嚴(yán)重,而頻域功率譜信號峰值頻率突出,抑制了噪聲的干擾。頻域相關(guān)程度較大三個(gè)IMF分量:IMF1,IMF2,IMF3正好,分別對應(yīng)仿真信號的三個(gè)特征頻率成分。因此,頻域互相關(guān)系數(shù)比時(shí)域互相關(guān)更能篩選出真實(shí)的IMF分量。具有更大的優(yōu)越性。

    表1 自相關(guān)處理的仿真信號與各IMF分量頻域相關(guān)系數(shù)

    表2 自相關(guān)處理的仿真信號與各IMF分量時(shí)域相關(guān)系數(shù)

    取IMF主分量:IMF1,IMF2,IMF3作為信號EEMD的分解結(jié)果(圖7),然后對信號進(jìn)行Hilbert譜分析(圖8)。從圖7看出經(jīng)D-CA和R-EEMD算法處理可以有效解決虛假IMF分量,模態(tài)混疊得到有效解決,從Hilbert譜可以看出端點(diǎn)效應(yīng)已基本消除,提高了信號的信噪比,三個(gè)頻率成分:5 Hz,10 Hz和30 Hz,特征明顯。圖9是未經(jīng)SVR延拓的EEMD分解的Hilbert譜,端點(diǎn)效應(yīng)明顯,影響EEMD對模態(tài)混疊的抑制,低頻混疊嚴(yán)重。圖10是未用DCA處理的EEMD分解的Hilbert時(shí)頻譜,虛假IMF增多,時(shí)頻譜混亂,過多的IMF分量使端點(diǎn)效應(yīng)的抑制效果不明顯。因此,圖9和圖10的Hilbert時(shí)頻譜都不能有效提取仿真信號的特征頻率。

    圖7 染噪信號D-CA和R-EEMD算法處理后的結(jié)果

    圖8 D-CA和R-EEMD處理后的Hilbert譜

    圖9 D-CA處理后的Hilbert譜

    圖10 R-EEMD處理后的Hilbert譜

    3 液壓故障信號分析

    實(shí)驗(yàn)對象為自主開發(fā)的液壓故障綜合實(shí)驗(yàn)平臺,技術(shù)規(guī)格及主要參數(shù)為:額定壓力8 MPa,額定流量3.7 L/min,電機(jī)功率1.5 KW,額定轉(zhuǎn)速3 000 r/ min。通過EMT690D采集數(shù)據(jù),它具有14位AD、4通道振動數(shù)據(jù)并行處理。采樣頻率為5 kHz,采樣點(diǎn)數(shù)為2 048,選取1 500個(gè)點(diǎn)作為分析對象。圖11是液壓系統(tǒng)某故障的原始波形,故障信號淹沒在振動信號中,無法直接識別。

    圖11 液壓故障信號原始波形

    3.1 實(shí)驗(yàn)信號的EEMD分析

    對信號直接進(jìn)行EEMD分解(圖12)得到9個(gè)IMF分量和一個(gè)殘余項(xiàng)。對分解得到的IMF分量進(jìn)行Hilbert變換得到信號的時(shí)頻譜(圖13)。時(shí)頻譜混疊嚴(yán)重、特征頻率不明顯,只能勉強(qiáng)提取100 Hz、240 Hz兩個(gè)頻率成分,其它故障信息不能有效提取。

    圖12 液壓故障信號直接EEMD分解圖

    圖13 液壓故障信號的Hilbert譜

    3.2 實(shí)驗(yàn)信號D-CA和R-EEMD分析

    運(yùn)用本文方法對故障信號進(jìn)行分析和處理,頻域互相關(guān)系數(shù)如表3,分解的主要IMF分量和時(shí)頻譜分別如圖14和圖15。與直接進(jìn)行EEMD分解比較,主要IMF數(shù)量大幅下降,利用頻率互相關(guān)剔除了與故障無關(guān)的IMF干擾,利用自相關(guān)函數(shù)去除了高頻噪聲,強(qiáng)化了故障信息,時(shí)頻圖清晰、不混亂,可以看到明顯的四個(gè)特征頻率:100 Hz、240 Hz、468 Hz和1 168 Hz,分別對應(yīng)故障信號功率譜(圖16)中四個(gè)主要頻率成分。經(jīng)過D-CA和R-EEMD處理后,液壓主要故障頻率全部順利提取,清晰表達(dá)了IMF的物理意義。

    圖14 故障信號D-CA和R-EEMD算法處理后的結(jié)果

    圖15 故障信號D-CA和R-EEMD處理后的Hilbert譜

    圖16 故障信號的功率譜

    表3 自相關(guān)處理后的故障信號與各IMF分量頻域相關(guān)系數(shù)

    4 結(jié)語

    針對液壓故障特征非線性、非平穩(wěn)、多噪聲的特點(diǎn),提出了基于D-CA和R-EEMD的信號處理方法,有效地提取了故障特征頻率。

    (1)利用對原始信號的自相關(guān)處理能有效地降噪和突出故障特征的周期分量。EEMD分解的IMF分量與自相關(guān)處理后的原信號作頻域互相關(guān)能減少了噪聲對相關(guān)系數(shù)的影響,比傳統(tǒng)時(shí)域互相關(guān)系數(shù)更有效地剔除虛假IMF分量,從而減少EEMD分解的IMF分量。利用SVR延拓改進(jìn)EEMD能抑制端點(diǎn)效應(yīng)和頻譜混疊。有效解決了影響EMD應(yīng)用中虛假分量、端點(diǎn)效應(yīng)和模態(tài)混疊問題;

    (2)通過與直接EEMD分解比較,仿真和實(shí)驗(yàn)結(jié)果表明,該方法能提高信號的信噪比,時(shí)頻譜更加清晰。虛假IMF分量明顯減少,兩側(cè)端點(diǎn)基本無波動,保證了IMF的物理意義,有效提取了信號中的特征頻率,揭示了液壓系統(tǒng)故障特征信息。

    [1]陳彥龍,張培林,徐超,等.基于DCT和EMD的滾動軸承故障診斷[J].電子測量技術(shù),2012,35(2):121-124.

    [2]雷亞國.基于改進(jìn)Hilbert-Huang變換的機(jī)械故障診斷[J].機(jī)械工程學(xué)報(bào),2011,47(5):71-77.

    [3]齊天,裘焱,吳亞峰.利用聚合經(jīng)驗(yàn)?zāi)B(tài)分解抑制振動信號中的模態(tài)混疊[J].噪聲與振動控制,2010,(2):103-106.

    [4]彭暢,柏林,謝小亮.基于EEMD、度量因子和快速譜峭度的滾動軸承故障診斷方法[J].振動與沖擊,2012,31 (20):143-146.

    [5]高強(qiáng),段晨東,趙艷青,等.基于最大波形相關(guān)延拓的經(jīng)驗(yàn)?zāi)J椒纸舛它c(diǎn)效應(yīng)抑制方法[J].振動與沖擊,2013,32 (2):62-66.

    [6]周龍寶,馬景輝,李強(qiáng).GA和RBF神經(jīng)網(wǎng)絡(luò)的Hilbert.Huang變換端點(diǎn)問題研究[J].計(jì)算機(jī)測量與控制,2012,20(2):447-450.

    [7]佟雨燕,陸森林.信號相關(guān)性和EEMD-Hilbert包絡(luò)在滾動軸承故障診斷中的應(yīng)用[J].噪聲與振動控制,2013,(5):144-149.

    [8]黃迪山.經(jīng)驗(yàn)?zāi)B(tài)分解中虛假模態(tài)分量消除法[J].振動、測試與診斷,2011,31(3):381-384.

    [9]黃姣英,袁海文,崔勇,等.基于EMD和自相關(guān)分析的輪胎音頻信號處理[J].電子測量與儀器學(xué)報(bào),2009,23(9):33-37.

    [10]劉樹林,趙海峰,喬波,等.基于經(jīng)驗(yàn)?zāi)B(tài)分解及自相關(guān)分析的微弱信號提取方法[J].大慶石油學(xué)院學(xué)報(bào),2007,31(5):80-84.

    [11]明安波,儲福磊,張煒.滾動軸承復(fù)合故障特征分離的小波頻譜自相關(guān)系數(shù)[J].機(jī)械工程學(xué)報(bào),2013,49(3):80-86.

    [12]張梅軍,陳灝,曹勤,等.基于SVM信號延拓的改進(jìn)的EEMD方法[J].振動、測試與診斷,2013,33(1):93-97.

    Fault Feature Extraction of Hydraulic Systems Based on D-CAand R-EEMD Methods

    CHAIKai,ZHANG Mei-jun,HUANG Jie,Feng Xia
    (College of Field Engineering,PLAUniv.of Sci.&Tech.,Nanjing 210007,China)

    A new method for hydraulic-system fault-feature extraction was proposed based on double-correlation analysis(D-CA)and refined ensemble empirical-mode-decomposition(R-EEMD).Firstly,the original signal was processed by adaptive correlation analysis to extrude the periodic components and eliminate the noise.Secondly,the R-EEMD based on support vector regression(SVR)was used to analyze the adaptive correlation processing signal to obtain the theoretical intrinsic mode function(IMF).Thirdly,the false IMF components were removed by extracting the correlation coefficient of IMFs and the adaptive correlation processing signal in the frequency domain instead of traditional time domain.Finally,the reconstructed signal was analyzed by the Hilbert spectrum to extract the fault features.Simulation and experimental results show that this method can increase signal-to-noise ratio,reduce the number of IMFs,depress the end effect,and effectively extract the faults feature frequencies of the hydraulic system.

    vibration and wave;fault diagnosis;double-correlation analysis;refined ensemble empirical mode decomposition

    TB53;TN911.7;TH165+.3

    :A

    :10.3969/j.issn.1006-1335.2015.01.042

    1006-1355(2014)05-0204-05+224

    2014-06-13

    2011年國家自然科學(xué)基金項(xiàng)目(51175511)

    柴凱(1989-),男,湖南益陽人,碩士研究生,主要研究方向:軍用裝備狀態(tài)監(jiān)測與故障診斷。

    張梅軍,女,碩士生導(dǎo)師。E-mail:dingdf0365@sina.com

    猜你喜歡
    端點(diǎn)頻域時(shí)域
    非特征端點(diǎn)條件下PM函數(shù)的迭代根
    不等式求解過程中端點(diǎn)的確定
    基于時(shí)域信號的三電平逆變器復(fù)合故障診斷
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設(shè)計(jì)
    參數(shù)型Marcinkiewicz積分算子及其交換子的加權(quán)端點(diǎn)估計(jì)
    基于極大似然準(zhǔn)則與滾動時(shí)域估計(jì)的自適應(yīng)UKF算法
    基丁能雖匹配延拓法LMD端點(diǎn)效應(yīng)處理
    基于改進(jìn)Radon-Wigner變換的目標(biāo)和拖曳式誘餌頻域分離
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    一種基于頻域的QPSK窄帶干擾抑制算法
    日韩 亚洲 欧美在线| 亚洲国产精品专区欧美| 色综合色国产| 日韩国内少妇激情av| 国产免费福利视频在线观看| 少妇人妻一区二区三区视频| 成人黄色视频免费在线看| 国产爱豆传媒在线观看| 老师上课跳d突然被开到最大视频| 一本久久精品| 有码 亚洲区| 亚洲av国产av综合av卡| 欧美人与善性xxx| 欧美日韩亚洲高清精品| 欧美日韩国产mv在线观看视频 | 99久久九九国产精品国产免费| 国产在线一区二区三区精| 午夜福利视频1000在线观看| 国产黄a三级三级三级人| 只有这里有精品99| 久久久久久久久久久免费av| 老师上课跳d突然被开到最大视频| 黄色视频在线播放观看不卡| 一级二级三级毛片免费看| 久久精品国产自在天天线| 精品国产一区二区三区久久久樱花 | 岛国毛片在线播放| 乱码一卡2卡4卡精品| 简卡轻食公司| 日本爱情动作片www.在线观看| 久久久久久伊人网av| av网站免费在线观看视频| 日韩国内少妇激情av| 免费黄频网站在线观看国产| 精品人妻熟女av久视频| 少妇裸体淫交视频免费看高清| 国产精品av视频在线免费观看| 寂寞人妻少妇视频99o| 免费av观看视频| 亚洲久久久久久中文字幕| 日本欧美国产在线视频| 舔av片在线| 菩萨蛮人人尽说江南好唐韦庄| 两个人的视频大全免费| 高清视频免费观看一区二区| 久久久久久久久久久丰满| 国产片特级美女逼逼视频| 欧美一区二区亚洲| 18禁裸乳无遮挡动漫免费视频 | 国产亚洲精品久久久com| 欧美性猛交╳xxx乱大交人| 一本色道久久久久久精品综合| 99久久精品一区二区三区| 国产精品伦人一区二区| 春色校园在线视频观看| 欧美国产精品一级二级三级 | 日本黄色片子视频| 能在线免费看毛片的网站| 亚洲精品国产成人久久av| 搡女人真爽免费视频火全软件| 边亲边吃奶的免费视频| 天天一区二区日本电影三级| 精品久久久噜噜| 如何舔出高潮| 在线 av 中文字幕| 免费观看性生交大片5| 亚洲自偷自拍三级| 国产真实伦视频高清在线观看| 大话2 男鬼变身卡| 亚洲自拍偷在线| 亚洲成人av在线免费| 国产精品国产三级国产专区5o| 欧美zozozo另类| 精品久久国产蜜桃| 日韩精品有码人妻一区| 午夜日本视频在线| 精品国产一区二区三区久久久樱花 | 99久久中文字幕三级久久日本| 午夜日本视频在线| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品一区蜜桃| 人妻系列 视频| 禁无遮挡网站| 啦啦啦中文免费视频观看日本| 一本色道久久久久久精品综合| 免费高清在线观看视频在线观看| 国产色婷婷99| 久久久色成人| 日韩强制内射视频| 18禁裸乳无遮挡动漫免费视频 | 亚洲精品第二区| 亚洲最大成人手机在线| 免费观看无遮挡的男女| 一边亲一边摸免费视频| av免费观看日本| 亚洲一区二区三区欧美精品 | 日日摸夜夜添夜夜爱| 一级毛片电影观看| 国产成人免费观看mmmm| 亚洲成人精品中文字幕电影| 大话2 男鬼变身卡| 国产精品一区二区在线观看99| kizo精华| 精品视频人人做人人爽| 久久精品国产a三级三级三级| 色视频www国产| 久久久久精品性色| 午夜老司机福利剧场| 精品少妇久久久久久888优播| 九色成人免费人妻av| 亚洲怡红院男人天堂| 精品人妻视频免费看| 久久精品人妻少妇| 亚洲电影在线观看av| 三级经典国产精品| 久久精品国产亚洲av涩爱| 亚洲人成网站在线观看播放| 国产成人a区在线观看| 亚洲人成网站在线播| 天堂俺去俺来也www色官网| 亚洲av免费高清在线观看| 久久韩国三级中文字幕| 在线观看美女被高潮喷水网站| 日本猛色少妇xxxxx猛交久久| 日日撸夜夜添| 看非洲黑人一级黄片| 免费高清在线观看视频在线观看| 三级国产精品欧美在线观看| 免费av观看视频| 中文字幕亚洲精品专区| av国产精品久久久久影院| 亚洲精品国产成人久久av| 国产男女内射视频| 中文精品一卡2卡3卡4更新| 一级毛片 在线播放| 国产伦精品一区二区三区视频9| 少妇人妻一区二区三区视频| 久久精品国产a三级三级三级| 男人爽女人下面视频在线观看| 少妇熟女欧美另类| 人人妻人人看人人澡| 国产免费一级a男人的天堂| 黄片无遮挡物在线观看| 国产免费一级a男人的天堂| 在线a可以看的网站| 女人久久www免费人成看片| 亚洲国产日韩一区二区| 在线精品无人区一区二区三 | 欧美一级a爱片免费观看看| 国产日韩欧美亚洲二区| 99久久人妻综合| 成年av动漫网址| 成人二区视频| 美女高潮的动态| 日本黄色片子视频| 国产精品一二三区在线看| 国产中年淑女户外野战色| 人妻少妇偷人精品九色| 男插女下体视频免费在线播放| 亚洲欧美日韩卡通动漫| 久久精品久久精品一区二区三区| 黄色怎么调成土黄色| 亚洲av一区综合| 国产精品精品国产色婷婷| 一级毛片电影观看| 亚洲精品国产色婷婷电影| 大片电影免费在线观看免费| xxx大片免费视频| 国产成人精品婷婷| 大片电影免费在线观看免费| 亚洲av福利一区| 亚洲av福利一区| 大片电影免费在线观看免费| 亚洲自偷自拍三级| 国产亚洲91精品色在线| 欧美精品国产亚洲| 精品午夜福利在线看| av又黄又爽大尺度在线免费看| 一本色道久久久久久精品综合| 亚洲精品乱久久久久久| 色吧在线观看| 亚洲av免费在线观看| 一本色道久久久久久精品综合| 一级毛片 在线播放| 九九久久精品国产亚洲av麻豆| 亚洲伊人久久精品综合| 日韩制服骚丝袜av| 日韩制服骚丝袜av| 久久久久久伊人网av| 高清日韩中文字幕在线| 日韩大片免费观看网站| 国产精品人妻久久久影院| 午夜激情久久久久久久| 国产一区二区亚洲精品在线观看| 丰满乱子伦码专区| 18禁在线播放成人免费| 深爱激情五月婷婷| 亚洲色图综合在线观看| av卡一久久| 亚洲欧美清纯卡通| 好男人在线观看高清免费视频| 熟女人妻精品中文字幕| 狂野欧美白嫩少妇大欣赏| 天堂俺去俺来也www色官网| 干丝袜人妻中文字幕| 在线亚洲精品国产二区图片欧美 | 国产欧美日韩一区二区三区在线 | 日本三级黄在线观看| 观看免费一级毛片| av专区在线播放| 91精品一卡2卡3卡4卡| 日韩国内少妇激情av| 成人亚洲精品一区在线观看 | 国产成人一区二区在线| 亚洲精品久久久久久婷婷小说| av一本久久久久| 午夜爱爱视频在线播放| 亚洲精品第二区| 极品教师在线视频| 亚洲国产精品国产精品| av在线老鸭窝| av专区在线播放| 国产精品99久久久久久久久| 三级经典国产精品| 大话2 男鬼变身卡| 日日啪夜夜撸| 26uuu在线亚洲综合色| 丰满人妻一区二区三区视频av| 有码 亚洲区| 观看美女的网站| 91精品伊人久久大香线蕉| 大香蕉97超碰在线| 亚洲精品成人av观看孕妇| 中文字幕制服av| 国产成人a区在线观看| 美女主播在线视频| 26uuu在线亚洲综合色| 亚洲av免费在线观看| 狂野欧美激情性xxxx在线观看| 一区二区三区四区激情视频| 一本久久精品| 秋霞在线观看毛片| 免费黄网站久久成人精品| 国产色婷婷99| 久久久久久久久久成人| 久久久久久久久久人人人人人人| 人妻夜夜爽99麻豆av| 国产精品偷伦视频观看了| 日韩制服骚丝袜av| 色哟哟·www| 日本爱情动作片www.在线观看| 久久久久九九精品影院| 麻豆成人午夜福利视频| 免费看av在线观看网站| 国产大屁股一区二区在线视频| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩东京热| 国产午夜福利久久久久久| av在线观看视频网站免费| 亚洲国产av新网站| 久久久精品94久久精品| 少妇猛男粗大的猛烈进出视频 | 免费av不卡在线播放| 国产乱人视频| 一区二区av电影网| 久久久久久久午夜电影| 91精品伊人久久大香线蕉| 国产精品一区www在线观看| 午夜免费男女啪啪视频观看| 人妻制服诱惑在线中文字幕| 在线看a的网站| 3wmmmm亚洲av在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 久热这里只有精品99| 女人久久www免费人成看片| 欧美精品人与动牲交sv欧美| eeuss影院久久| 精品久久久久久久久av| 校园人妻丝袜中文字幕| 男人和女人高潮做爰伦理| 婷婷色av中文字幕| 欧美一级a爱片免费观看看| 日本熟妇午夜| 国内少妇人妻偷人精品xxx网站| 在线观看一区二区三区| 国产精品三级大全| 丝袜美腿在线中文| 狂野欧美白嫩少妇大欣赏| av线在线观看网站| 亚洲精品一二三| 日本免费在线观看一区| 街头女战士在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 在线观看美女被高潮喷水网站| 禁无遮挡网站| 免费看日本二区| 日韩欧美一区视频在线观看 | av国产精品久久久久影院| 下体分泌物呈黄色| 18+在线观看网站| 国产成人91sexporn| 99热这里只有是精品在线观看| 精品国产露脸久久av麻豆| 欧美性猛交╳xxx乱大交人| 国产久久久一区二区三区| 亚洲无线观看免费| 亚洲第一区二区三区不卡| 精华霜和精华液先用哪个| 一级片'在线观看视频| 性插视频无遮挡在线免费观看| 日韩一区二区三区影片| 天堂网av新在线| 成人免费观看视频高清| 国产黄色免费在线视频| 日韩大片免费观看网站| 偷拍熟女少妇极品色| 五月天丁香电影| av线在线观看网站| 色5月婷婷丁香| 在线天堂最新版资源| 亚洲熟女精品中文字幕| 国产亚洲91精品色在线| 黄色怎么调成土黄色| 中文乱码字字幕精品一区二区三区| 午夜精品国产一区二区电影 | 人妻夜夜爽99麻豆av| 国产黄频视频在线观看| 免费看a级黄色片| 日韩制服骚丝袜av| 精品一区二区免费观看| 天美传媒精品一区二区| 国产精品99久久99久久久不卡 | 97超视频在线观看视频| 69人妻影院| 国产黄色视频一区二区在线观看| 九九爱精品视频在线观看| 看黄色毛片网站| 永久网站在线| 久久精品国产自在天天线| 激情 狠狠 欧美| 男人爽女人下面视频在线观看| 99久久人妻综合| 国产探花极品一区二区| 亚洲精品视频女| 免费观看av网站的网址| 在线观看一区二区三区| 国产综合精华液| 人妻系列 视频| 看黄色毛片网站| av在线蜜桃| 国产精品一及| 国产成人免费无遮挡视频| 99热全是精品| 欧美日本视频| 亚洲av欧美aⅴ国产| 一本色道久久久久久精品综合| av在线播放精品| 久久精品夜色国产| 青春草视频在线免费观看| h日本视频在线播放| 国产午夜精品久久久久久一区二区三区| 日日撸夜夜添| av卡一久久| 超碰97精品在线观看| 日韩在线高清观看一区二区三区| 亚洲av欧美aⅴ国产| av又黄又爽大尺度在线免费看| 亚洲国产最新在线播放| 蜜桃久久精品国产亚洲av| 不卡视频在线观看欧美| 日本与韩国留学比较| 亚洲精品国产色婷婷电影| 国产美女午夜福利| 免费人成在线观看视频色| 人妻夜夜爽99麻豆av| 中文精品一卡2卡3卡4更新| 亚洲四区av| 国产精品麻豆人妻色哟哟久久| 熟女电影av网| 七月丁香在线播放| 一本久久精品| 国产 一区 欧美 日韩| 国产老妇女一区| 免费电影在线观看免费观看| 日韩欧美精品免费久久| 国产日韩欧美亚洲二区| 精品国产一区二区三区久久久樱花 | 久久精品久久久久久久性| 高清视频免费观看一区二区| 亚洲成人精品中文字幕电影| 少妇裸体淫交视频免费看高清| 亚洲成人一二三区av| 日韩,欧美,国产一区二区三区| 啦啦啦在线观看免费高清www| 韩国高清视频一区二区三区| 日韩av在线免费看完整版不卡| 久久久久精品性色| 国产精品一区二区性色av| 一区二区三区乱码不卡18| av天堂中文字幕网| 欧美zozozo另类| 欧美激情国产日韩精品一区| 我要看日韩黄色一级片| 午夜免费观看性视频| 国产精品国产三级国产av玫瑰| 国产黄色视频一区二区在线观看| 美女被艹到高潮喷水动态| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美在线一区| 乱系列少妇在线播放| 免费av观看视频| 亚洲欧美日韩无卡精品| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 搞女人的毛片| 少妇被粗大猛烈的视频| videossex国产| 国产成人精品一,二区| 两个人的视频大全免费| 91久久精品国产一区二区三区| 国产69精品久久久久777片| 日本熟妇午夜| 久久99热这里只有精品18| 春色校园在线视频观看| 美女被艹到高潮喷水动态| 亚洲真实伦在线观看| 日日摸夜夜添夜夜爱| 尤物成人国产欧美一区二区三区| 亚洲国产高清在线一区二区三| 一边亲一边摸免费视频| 国产视频内射| 五月开心婷婷网| 亚洲精品乱久久久久久| 黄色视频在线播放观看不卡| 色5月婷婷丁香| 国产精品av视频在线免费观看| 高清欧美精品videossex| 人妻 亚洲 视频| 亚洲天堂av无毛| 精品久久久精品久久久| 亚洲丝袜综合中文字幕| 自拍偷自拍亚洲精品老妇| 久久韩国三级中文字幕| 在线观看一区二区三区激情| 亚洲精品视频女| 国精品久久久久久国模美| 最近最新中文字幕大全电影3| 中文天堂在线官网| 一级毛片久久久久久久久女| 日本三级黄在线观看| 1000部很黄的大片| 一级二级三级毛片免费看| 久久精品国产亚洲av天美| 寂寞人妻少妇视频99o| 中国国产av一级| 亚洲av在线观看美女高潮| 国产精品一区二区性色av| av在线天堂中文字幕| 亚洲av在线观看美女高潮| 欧美精品人与动牲交sv欧美| 日本av手机在线免费观看| 乱码一卡2卡4卡精品| 国产黄色免费在线视频| av卡一久久| 久久久精品免费免费高清| 国产精品偷伦视频观看了| 97热精品久久久久久| 99久久精品一区二区三区| av国产免费在线观看| 2021少妇久久久久久久久久久| 91精品一卡2卡3卡4卡| 久久99热这里只有精品18| 国产成人aa在线观看| 欧美zozozo另类| 亚洲美女搞黄在线观看| 纵有疾风起免费观看全集完整版| 黄色怎么调成土黄色| 久久久久久久亚洲中文字幕| 久久久午夜欧美精品| 黄色怎么调成土黄色| 亚洲人成网站高清观看| 自拍偷自拍亚洲精品老妇| 精品久久久精品久久久| 欧美高清成人免费视频www| 亚洲国产av新网站| 99热6这里只有精品| 老司机影院毛片| 久久久精品欧美日韩精品| 国产乱人偷精品视频| 午夜福利在线在线| 一区二区三区乱码不卡18| 亚洲欧美成人精品一区二区| 久久久久久伊人网av| 久久99热6这里只有精品| 国产有黄有色有爽视频| 青春草亚洲视频在线观看| 丝瓜视频免费看黄片| 亚洲精品日韩在线中文字幕| 国产精品久久久久久av不卡| 国产午夜精品一二区理论片| 久久99热这里只有精品18| 人妻系列 视频| av.在线天堂| 亚洲,一卡二卡三卡| 黄色一级大片看看| 91在线精品国自产拍蜜月| 一个人看视频在线观看www免费| 激情五月婷婷亚洲| 中文字幕人妻熟人妻熟丝袜美| 免费看a级黄色片| 大香蕉97超碰在线| 久久久久国产网址| 久久久精品94久久精品| 国产日韩欧美亚洲二区| 少妇人妻一区二区三区视频| 在线 av 中文字幕| 久久女婷五月综合色啪小说 | 男人爽女人下面视频在线观看| 国产成人精品福利久久| 狂野欧美激情性xxxx在线观看| 免费黄频网站在线观看国产| 免费观看a级毛片全部| 日韩 亚洲 欧美在线| 国产精品久久久久久久电影| 免费观看性生交大片5| 看黄色毛片网站| 亚州av有码| 久久精品国产亚洲av天美| 国产精品福利在线免费观看| 久久久久网色| 欧美区成人在线视频| 久久97久久精品| 欧美高清成人免费视频www| 黑人高潮一二区| av线在线观看网站| 哪个播放器可以免费观看大片| 最近的中文字幕免费完整| 日本-黄色视频高清免费观看| 久久韩国三级中文字幕| 人妻一区二区av| 97在线视频观看| 国产大屁股一区二区在线视频| 久久精品国产自在天天线| 国产精品国产三级国产专区5o| 在线看a的网站| 午夜老司机福利剧场| 国产毛片a区久久久久| 美女xxoo啪啪120秒动态图| 美女被艹到高潮喷水动态| 我的女老师完整版在线观看| 亚洲av不卡在线观看| 国产精品偷伦视频观看了| 日韩 亚洲 欧美在线| av在线蜜桃| 禁无遮挡网站| 美女脱内裤让男人舔精品视频| 一二三四中文在线观看免费高清| 人妻一区二区av| 亚洲精品国产色婷婷电影| 一级毛片aaaaaa免费看小| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品乱码久久久v下载方式| 亚洲精品日韩av片在线观看| 青春草国产在线视频| 亚洲精品亚洲一区二区| 中文乱码字字幕精品一区二区三区| 九九在线视频观看精品| 久久韩国三级中文字幕| 成人高潮视频无遮挡免费网站| 欧美日韩视频高清一区二区三区二| 国产午夜精品一二区理论片| 精品久久国产蜜桃| 国产淫语在线视频| 国产男人的电影天堂91| 男女那种视频在线观看| 七月丁香在线播放| 高清在线视频一区二区三区| 亚洲欧洲国产日韩| 在线亚洲精品国产二区图片欧美 | 黄片无遮挡物在线观看| 建设人人有责人人尽责人人享有的 | 亚洲,一卡二卡三卡| www.av在线官网国产| 久久精品久久精品一区二区三区| 麻豆成人av视频| 18禁在线播放成人免费| 日韩成人伦理影院| 久久影院123| 尾随美女入室| 久久综合国产亚洲精品| 亚洲性久久影院| 欧美日韩视频高清一区二区三区二| 久久亚洲国产成人精品v| 久久久精品欧美日韩精品| 97超碰精品成人国产| 我要看日韩黄色一级片| 欧美日本视频| 久久人人爽av亚洲精品天堂 | 另类亚洲欧美激情| 精品一区在线观看国产| 国产又色又爽无遮挡免| 久久精品综合一区二区三区| 精品久久久久久久人妻蜜臀av| 欧美另类一区| 欧美高清成人免费视频www| 舔av片在线| 特级一级黄色大片| 69人妻影院| 亚洲色图av天堂| 成人无遮挡网站| 亚洲第一区二区三区不卡| 色婷婷久久久亚洲欧美| 91久久精品电影网| 午夜激情福利司机影院| 一区二区三区精品91|