• <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免费视频播放| 欧美xxⅹ黑人| 男女午夜视频在线观看 | 欧美亚洲 丝袜 人妻 在线| 成年人免费黄色播放视频| 精品亚洲乱码少妇综合久久| 男男h啪啪无遮挡| 好男人视频免费观看在线| 啦啦啦视频在线资源免费观看| 亚洲色图综合在线观看| 黑人欧美特级aaaaaa片| 久久97久久精品| 久久久久久久精品精品| 日日爽夜夜爽网站| 最近中文字幕高清免费大全6| 久久韩国三级中文字幕| 亚洲成人手机| 99久久综合免费| 欧美 亚洲 国产 日韩一| 久久精品国产综合久久久 | 十八禁高潮呻吟视频| 极品少妇高潮喷水抽搐| 香蕉精品网在线| 男女午夜视频在线观看 | 国产亚洲午夜精品一区二区久久| 日韩中文字幕视频在线看片| 18禁裸乳无遮挡动漫免费视频| 午夜免费鲁丝| 黄色毛片三级朝国网站| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av天美| 又黄又爽又刺激的免费视频.| 看免费成人av毛片| 2022亚洲国产成人精品| 国产日韩欧美在线精品| 亚洲国产精品999| 天天躁夜夜躁狠狠久久av| 2022亚洲国产成人精品| 成人黄色视频免费在线看| 中文精品一卡2卡3卡4更新| 人人妻人人添人人爽欧美一区卜| 欧美精品av麻豆av| 国产片内射在线| av卡一久久| 亚洲五月色婷婷综合| 在线观看www视频免费| 久久久久国产精品人妻一区二区| 热re99久久精品国产66热6| 国产一级毛片在线| 伊人久久国产一区二区| 精品视频人人做人人爽| 国产伦理片在线播放av一区| 成人亚洲精品一区在线观看| 中文字幕另类日韩欧美亚洲嫩草| 香蕉丝袜av| 美女大奶头黄色视频| 18在线观看网站| 美女主播在线视频| 制服诱惑二区| 一本—道久久a久久精品蜜桃钙片| 黑人猛操日本美女一级片| 国产成人a∨麻豆精品| 99热全是精品| 亚洲久久久国产精品| 国产精品偷伦视频观看了| 青春草亚洲视频在线观看| 一级黄片播放器| 亚洲精品一二三| 男的添女的下面高潮视频| 一级a做视频免费观看| 曰老女人黄片| kizo精华| 国产精品一国产av| 99久久综合免费| 亚洲精品一区蜜桃| 国产女主播在线喷水免费视频网站| 久久久国产一区二区| 国产老妇伦熟女老妇高清| 国产xxxxx性猛交| 边亲边吃奶的免费视频| 欧美精品av麻豆av| 精品99又大又爽又粗少妇毛片| 五月开心婷婷网| 国产一区二区在线观看日韩| 亚洲精品国产av成人精品| 国产一区有黄有色的免费视频| av.在线天堂| xxx大片免费视频| 成年人午夜在线观看视频| 男女无遮挡免费网站观看| 肉色欧美久久久久久久蜜桃| 日韩三级伦理在线观看| 卡戴珊不雅视频在线播放| 另类精品久久| 亚洲中文av在线| 最新中文字幕久久久久| 日本欧美国产在线视频| 久久久精品免费免费高清| 日韩av在线免费看完整版不卡| 狂野欧美激情性xxxx在线观看| 三级国产精品片| 亚洲成av片中文字幕在线观看 | 亚洲综合色惰| 亚洲av免费高清在线观看| 男女无遮挡免费网站观看| 1024视频免费在线观看| 人人妻人人添人人爽欧美一区卜| 午夜福利视频在线观看免费| 亚洲精品日韩在线中文字幕| 男男h啪啪无遮挡| 老女人水多毛片| 亚洲av日韩在线播放| 夫妻午夜视频| 如何舔出高潮| 精品国产露脸久久av麻豆| 黑人高潮一二区| 国产片特级美女逼逼视频| 考比视频在线观看| 欧美bdsm另类| 咕卡用的链子| 日韩欧美一区视频在线观看| 国产一区二区在线观看日韩| 2021少妇久久久久久久久久久| 成人综合一区亚洲| 欧美日韩精品成人综合77777| 久久国产亚洲av麻豆专区| 国产在线一区二区三区精| 亚洲欧美精品自产自拍| 国国产精品蜜臀av免费| 亚洲国产成人一精品久久久| 最近最新中文字幕大全免费视频 | www.av在线官网国产| 日韩制服骚丝袜av| 国产成人aa在线观看| 尾随美女入室| 18禁裸乳无遮挡动漫免费视频| 国产精品99久久99久久久不卡 | 国产av一区二区精品久久| 黄片无遮挡物在线观看| 咕卡用的链子| 欧美激情 高清一区二区三区| 久久精品久久久久久噜噜老黄| 亚洲欧美清纯卡通| 亚洲精华国产精华液的使用体验| 精品久久久久久电影网| 51国产日韩欧美| 久久久久久久久久成人| 久久久久视频综合| 成人二区视频| 寂寞人妻少妇视频99o| 一级毛片电影观看| 国产亚洲精品第一综合不卡 | 晚上一个人看的免费电影| 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 国产精品国产三级国产专区5o| 日本黄色日本黄色录像| 午夜久久久在线观看| 久久精品久久久久久久性| 日韩,欧美,国产一区二区三区| 色吧在线观看| 欧美亚洲日本最大视频资源| 欧美日韩亚洲高清精品| 亚洲av欧美aⅴ国产| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 欧美xxⅹ黑人| 久久 成人 亚洲| 国产精品久久久av美女十八| 免费高清在线观看视频在线观看| 亚洲精品国产av成人精品| 午夜免费男女啪啪视频观看| 亚洲综合色惰| 亚洲精品第二区| 欧美激情 高清一区二区三区| 亚洲欧洲国产日韩| 亚洲国产色片| 在线天堂中文资源库| 日韩大片免费观看网站| 欧美日韩亚洲高清精品| 亚洲精品视频女| 成人漫画全彩无遮挡| 熟妇人妻不卡中文字幕| 视频在线观看一区二区三区| 嫩草影院入口| 高清视频免费观看一区二区| 日韩av免费高清视频| 免费黄色在线免费观看| 国产av码专区亚洲av| 精品一区在线观看国产| 91成人精品电影| 日韩一本色道免费dvd| av在线播放精品| 日韩视频在线欧美| 观看av在线不卡| 在现免费观看毛片| 成人亚洲精品一区在线观看| 在线观看国产h片| 丝瓜视频免费看黄片| 久久99一区二区三区| 黄片播放在线免费| 日韩一区二区视频免费看| 美女主播在线视频| 麻豆乱淫一区二区| 寂寞人妻少妇视频99o| 天美传媒精品一区二区| 午夜视频国产福利| 国产精品秋霞免费鲁丝片| 熟女人妻精品中文字幕| 亚洲经典国产精华液单| 视频中文字幕在线观看| 精品视频人人做人人爽| 桃花免费在线播放| 久久精品久久久久久久性| 中文字幕人妻丝袜制服| 天堂8中文在线网| 男女午夜视频在线观看 | 日日啪夜夜爽| 哪个播放器可以免费观看大片| 春色校园在线视频观看| 国产精品.久久久| 男女国产视频网站| 夜夜爽夜夜爽视频| 日韩制服骚丝袜av| 久热久热在线精品观看| freevideosex欧美| 亚洲精品日本国产第一区| 捣出白浆h1v1| 国产片内射在线| 精品亚洲成a人片在线观看| 日产精品乱码卡一卡2卡三| 99九九在线精品视频| 熟女av电影| 18禁裸乳无遮挡动漫免费视频| 婷婷色av中文字幕| 日韩av不卡免费在线播放| 亚洲成人手机| 亚洲国产日韩一区二区| 免费久久久久久久精品成人欧美视频 | 免费黄网站久久成人精品| 天天躁夜夜躁狠狠久久av| 精品国产国语对白av| 尾随美女入室| 男女高潮啪啪啪动态图| 午夜影院在线不卡| 七月丁香在线播放| 日韩制服丝袜自拍偷拍| 两性夫妻黄色片 | 夫妻性生交免费视频一级片| 久久青草综合色| 午夜免费观看性视频| 日韩制服骚丝袜av| 啦啦啦啦在线视频资源| 亚洲精品乱码久久久久久按摩| 欧美日韩国产mv在线观看视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品久久国产蜜桃| 亚洲欧美色中文字幕在线| 看非洲黑人一级黄片| 最近2019中文字幕mv第一页| 最近最新中文字幕大全免费视频 | a级毛片在线看网站| 久久精品国产综合久久久 | 日韩免费高清中文字幕av| 久久精品国产亚洲av涩爱| 国产精品欧美亚洲77777| av线在线观看网站| 99热国产这里只有精品6| 少妇熟女欧美另类| 伊人亚洲综合成人网| 亚洲婷婷狠狠爱综合网| 女人久久www免费人成看片| 欧美精品av麻豆av| 一级爰片在线观看| av片东京热男人的天堂| 亚洲av综合色区一区| 人成视频在线观看免费观看| 夜夜骑夜夜射夜夜干| 国产在线免费精品| 黑人欧美特级aaaaaa片| 熟女av电影| 国产精品免费大片| 亚洲欧美色中文字幕在线| 成人国产麻豆网| 青春草国产在线视频| 男女免费视频国产| 99久国产av精品国产电影| 亚洲中文av在线| 欧美日韩视频高清一区二区三区二| 午夜福利乱码中文字幕| 亚洲内射少妇av| 伊人久久国产一区二区| 人人妻人人澡人人看| 欧美3d第一页| 免费人妻精品一区二区三区视频| 少妇人妻精品综合一区二区| 国产永久视频网站| 水蜜桃什么品种好| 国产免费一级a男人的天堂| av国产久精品久网站免费入址| 欧美人与善性xxx| 激情五月婷婷亚洲| 久久这里有精品视频免费| 亚洲成国产人片在线观看| 狠狠精品人妻久久久久久综合| 国产精品国产三级专区第一集| 下体分泌物呈黄色| 男的添女的下面高潮视频| 国产亚洲精品久久久com| 成人亚洲精品一区在线观看| 最近中文字幕高清免费大全6| 97超碰精品成人国产| 国产精品三级大全| 人妻 亚洲 视频| 高清av免费在线| 久久久国产欧美日韩av| 成人国语在线视频| 国产高清三级在线| 五月玫瑰六月丁香| 久久国产精品男人的天堂亚洲 | 成人午夜精彩视频在线观看| 国产激情久久老熟女| 久久久国产欧美日韩av| 久久久久视频综合| 日韩视频在线欧美| 2018国产大陆天天弄谢| 美女中出高潮动态图| 国产精品国产三级国产av玫瑰| 国产爽快片一区二区三区| 亚洲天堂av无毛| 久久久精品区二区三区| av网站免费在线观看视频| 女性被躁到高潮视频| 免费女性裸体啪啪无遮挡网站| 午夜91福利影院| 日本黄色日本黄色录像| 日日爽夜夜爽网站| 色94色欧美一区二区| 中文字幕av电影在线播放| 在线精品无人区一区二区三| 亚洲 欧美一区二区三区| 国产伦理片在线播放av一区| 综合色丁香网| 国产欧美亚洲国产| av卡一久久| 看十八女毛片水多多多| 亚洲成人手机| 99视频精品全部免费 在线| 毛片一级片免费看久久久久| 久久女婷五月综合色啪小说| 9191精品国产免费久久| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 肉色欧美久久久久久久蜜桃| 在线 av 中文字幕| 一级,二级,三级黄色视频| 国产精品一区www在线观看| 天天操日日干夜夜撸| 成人无遮挡网站| 国产成人精品福利久久| 9191精品国产免费久久| av.在线天堂| 春色校园在线视频观看| 男女下面插进去视频免费观看 | 一区在线观看完整版| 欧美亚洲日本最大视频资源| 久久久久精品久久久久真实原创| 九九爱精品视频在线观看| 最新中文字幕久久久久| av在线app专区| 七月丁香在线播放| 国产色婷婷99| 国产亚洲欧美精品永久| 日韩一区二区视频免费看| 色哟哟·www| 成人午夜精彩视频在线观看| 欧美性感艳星| www.av在线官网国产| 少妇猛男粗大的猛烈进出视频| 免费大片黄手机在线观看| 日日摸夜夜添夜夜爱| 国产一区二区激情短视频 | 91成人精品电影| 成人国语在线视频| 成年动漫av网址| 日本猛色少妇xxxxx猛交久久| 久久精品aⅴ一区二区三区四区 | 精品国产国语对白av| 99国产综合亚洲精品| 成人无遮挡网站| 中文字幕人妻丝袜制服| 人妻少妇偷人精品九色| 免费人妻精品一区二区三区视频| 飞空精品影院首页| 日本vs欧美在线观看视频| 亚洲少妇的诱惑av| 亚洲av成人精品一二三区| 热re99久久精品国产66热6| 伊人亚洲综合成人网| 久久精品久久久久久噜噜老黄| 如日韩欧美国产精品一区二区三区| 曰老女人黄片| 99精国产麻豆久久婷婷| 黑人欧美特级aaaaaa片| 美女中出高潮动态图| 色婷婷av一区二区三区视频| 久热久热在线精品观看| 国产极品天堂在线| 欧美 亚洲 国产 日韩一| 女人久久www免费人成看片| 日本猛色少妇xxxxx猛交久久| 少妇精品久久久久久久| 国产一区二区在线观看日韩| 久久99热6这里只有精品| 亚洲人成77777在线视频| 在线观看美女被高潮喷水网站| 精品一区二区三卡| 少妇人妻久久综合中文| 国产精品久久久久久久久免| 美女大奶头黄色视频| 18禁在线无遮挡免费观看视频| 久久亚洲国产成人精品v| 各种免费的搞黄视频| 亚洲国产成人一精品久久久| 天天躁夜夜躁狠狠躁躁| 母亲3免费完整高清在线观看 | 黄片无遮挡物在线观看| 水蜜桃什么品种好| 免费观看av网站的网址| 毛片一级片免费看久久久久| 99re6热这里在线精品视频| 视频区图区小说| 高清欧美精品videossex| 久久精品人人爽人人爽视色| 中国国产av一级| 青春草视频在线免费观看| 最黄视频免费看| 国产免费一区二区三区四区乱码| 热99久久久久精品小说推荐| 日日撸夜夜添| 国产精品一区www在线观看| 满18在线观看网站| 一边摸一边做爽爽视频免费| 观看av在线不卡| 精品人妻一区二区三区麻豆| 国产男女内射视频| 在线精品无人区一区二区三| av国产久精品久网站免费入址| 亚洲欧美一区二区三区黑人 | 日韩在线高清观看一区二区三区| 欧美成人精品欧美一级黄| 久久精品国产自在天天线| 97在线人人人人妻| 一边亲一边摸免费视频| 亚洲,欧美,日韩| 久久国产亚洲av麻豆专区| 全区人妻精品视频| 在线免费观看不下载黄p国产| 黄片无遮挡物在线观看| 日日啪夜夜爽| 欧美另类一区| 在线天堂中文资源库| 美国免费a级毛片| 男女下面插进去视频免费观看 | 26uuu在线亚洲综合色| 伊人久久国产一区二区| 多毛熟女@视频| 十八禁网站网址无遮挡| 亚洲综合精品二区| a级毛片在线看网站| 国产成人欧美| 一本—道久久a久久精品蜜桃钙片| 日韩欧美一区视频在线观看| av国产精品久久久久影院| 国产有黄有色有爽视频| 51国产日韩欧美| 五月伊人婷婷丁香| 乱人伦中国视频| 大片电影免费在线观看免费| 午夜福利,免费看| 9191精品国产免费久久| 久久久亚洲精品成人影院| 免费看av在线观看网站| 99久久中文字幕三级久久日本| 久热这里只有精品99| 观看美女的网站| 久久久精品区二区三区| 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 亚洲精品久久成人aⅴ小说| 久久久久久久精品精品| 黑丝袜美女国产一区| 水蜜桃什么品种好| 26uuu在线亚洲综合色| 亚洲第一av免费看| 天堂8中文在线网| 人体艺术视频欧美日本| 少妇 在线观看| 精品国产乱码久久久久久小说| 午夜影院在线不卡| 观看美女的网站| 在线天堂最新版资源| 一二三四在线观看免费中文在 | 一级黄片播放器| 成人国产麻豆网| 亚洲丝袜综合中文字幕| 久久青草综合色| 亚洲欧美精品自产自拍| 超碰97精品在线观看| 成人国产av品久久久| 国产无遮挡羞羞视频在线观看| 亚洲国产av新网站| 国产午夜精品一二区理论片| 丝袜美足系列| 亚洲精品视频女| 久久精品国产亚洲av天美| 精品久久国产蜜桃| 国产精品.久久久| 曰老女人黄片| 免费播放大片免费观看视频在线观看| 国产淫语在线视频| 高清不卡的av网站| 精品福利永久在线观看| 女的被弄到高潮叫床怎么办| 2018国产大陆天天弄谢| 黄色一级大片看看| 黄色视频在线播放观看不卡| 亚洲欧洲国产日韩| 亚洲人成网站在线观看播放| 亚洲婷婷狠狠爱综合网| 这个男人来自地球电影免费观看 | 亚洲成人一二三区av| 国产精品国产三级国产av玫瑰| 亚洲国产精品999| 日韩免费高清中文字幕av| av在线观看视频网站免费| av在线app专区| 国产日韩一区二区三区精品不卡| 国产有黄有色有爽视频| 妹子高潮喷水视频| 五月玫瑰六月丁香| 在线天堂中文资源库| 国产精品久久久久久精品电影小说| 少妇被粗大猛烈的视频| 久久婷婷青草| 男男h啪啪无遮挡| 亚洲精品视频女| 亚洲少妇的诱惑av| 9色porny在线观看| 亚洲欧美一区二区三区黑人 | 国产在线免费精品| 九色亚洲精品在线播放| 少妇人妻 视频| 国产又色又爽无遮挡免| av国产久精品久网站免费入址| 日产精品乱码卡一卡2卡三| 国产亚洲av片在线观看秒播厂| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 免费看av在线观看网站| 啦啦啦在线观看免费高清www| 精品人妻熟女毛片av久久网站| av不卡在线播放| 99香蕉大伊视频| 男女无遮挡免费网站观看| 久久韩国三级中文字幕| freevideosex欧美| 少妇精品久久久久久久| 久久青草综合色| 精品人妻偷拍中文字幕| 久久精品国产a三级三级三级| 久久午夜综合久久蜜桃| 永久网站在线| 香蕉精品网在线| √禁漫天堂资源中文www| 伦理电影大哥的女人| 在线免费观看不下载黄p国产| 桃花免费在线播放| 十八禁网站网址无遮挡| 黄色视频在线播放观看不卡| 色婷婷av一区二区三区视频| 欧美人与善性xxx| av女优亚洲男人天堂| 亚洲国产av影院在线观看| 亚洲欧美成人精品一区二区| 精品久久久精品久久久| 亚洲天堂av无毛| 成人午夜精彩视频在线观看| 最后的刺客免费高清国语| 欧美日韩视频精品一区| 久久久精品免费免费高清| 最后的刺客免费高清国语| 男人操女人黄网站| 日韩精品免费视频一区二区三区 | 不卡视频在线观看欧美| 曰老女人黄片| 亚洲综合色惰| 免费看光身美女| 午夜激情av网站| 亚洲国产精品专区欧美| 国产成人一区二区在线| 下体分泌物呈黄色| 亚洲av男天堂| 国产在线视频一区二区| 97精品久久久久久久久久精品| 国产精品国产三级国产av玫瑰| www.av在线官网国产| 麻豆乱淫一区二区| 亚洲精品乱久久久久久| 女人久久www免费人成看片| 99热这里只有是精品在线观看| 高清av免费在线| 国产男女超爽视频在线观看|