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

    基于時間序列差分特性的LMD端點效應(yīng)處理方法

    2016-06-04 08:24:36杜冬梅
    中國機械工程 2016年10期

    杜冬梅 張 昭 何 青 李 紅

    華北電力大學(xué),北京,102206

    ?

    基于時間序列差分特性的LMD端點效應(yīng)處理方法

    杜冬梅張昭何青李紅

    華北電力大學(xué),北京,102206

    摘要:局部均值分解方法在計算過程中會使信號在端部產(chǎn)生端點效應(yīng),為減弱LMD中存在的端點效應(yīng),提出了一種基于信號時間序列差分特性的端點延拓方法。將信號曲線看作一段離散函數(shù),利用該函數(shù)段末端的一階差分和二階差分的符號判斷曲線的特征和發(fā)展趨勢,根據(jù)具體情況判定所延拓極值點的值。使用該方法對仿真信號和實驗軸承外圈故障和內(nèi)圈故障信號進行了分解,并將其結(jié)果與鏡像延拓改進的LMD方法分解的結(jié)果進行了對比。該方法在分解中的良好表現(xiàn)驗證了其在信號分析中的有效性。

    關(guān)鍵詞:局部均值分解;端點效應(yīng);時間序列差分;延拓

    0引言

    振動是反映機械設(shè)備運行狀態(tài)的重要參數(shù),通過對機械設(shè)備振動信號的監(jiān)測分析,可以對設(shè)備的運行狀態(tài)進行監(jiān)視,及時發(fā)現(xiàn)設(shè)備出現(xiàn)的異常狀態(tài),并對其原因作出分析判斷。時頻分析可以更好地判斷設(shè)備發(fā)生異常時的振動頻率變化,傳統(tǒng)的時頻分析方法有短時Fourier變換[1]、Gabor變換[2]、 Wignar-Ville分布[3]和小波變換[4-5]。短時Fourier變換、Gabor變換、Wignar-Ville分布和小波變換是在Fourier變換的基礎(chǔ)上利用窗函數(shù)將信號劃分為若干段分別求子段頻率的頻譜分析方法,它們可以確定頻率發(fā)生變化的時刻。但是利用窗函數(shù)對信號進行機械的劃分,不能自適應(yīng)地將信號中包含的固有模態(tài)成分分離出來,在信號分析中有一定的局限性。

    經(jīng)驗?zāi)B(tài)分解(empirical mode decomposition, EMD)是Huang等[6]提出的一種自適應(yīng)的分解方法,能夠自適應(yīng)地將信號分解為一組頻率從高到低自動排列的本征模態(tài)分量和一個殘余分量,但EMD本身存在的端點效應(yīng)、模態(tài)混疊等問題較嚴(yán)重。Simth[7]提出了一種新的信號分解算法——局部均值分解(local mean decomposition, LMD),采用移動平均的方式求取信號包絡(luò)線和均值,減小了EMD因樣條插值造成的嚴(yán)重的端點效應(yīng),同時瞬時頻率具有了明確的物理意義。學(xué)者們對LMD作了許多研究,Yang等[8]利用高斯白噪聲具有頻率均勻分布的統(tǒng)計特性,在目標(biāo)信號中加入白噪聲以減弱LMD算法的模態(tài)混疊現(xiàn)象。Park等[9]將LMD引入復(fù)數(shù)域,從而使LMD算法可以進行多變量的信號分析。任達千等[10]將LMD的端點效應(yīng)與EMD作了比較,并提出了一種基于能量的端點效應(yīng)評價方法和鏡像延拓的端點效應(yīng)抑制方法,但鏡像延拓的極值點確定缺乏具體的理論依據(jù)。孟宗等[11]提出了基于對稱差分能量算子解調(diào)的LMD端點效應(yīng)抑制方法,實現(xiàn)了旋轉(zhuǎn)機械故障的有效診斷,但LMD分解過程中產(chǎn)生的各分量的端點效應(yīng)難以避免。

    本文針對LMD算法中存在的端點效應(yīng)問題,提出了基于信號時間序列差分特性的方法進行端點延拓,根據(jù)信號兩端時間序列的一階差分和二階差分的符號判斷信號的發(fā)展趨勢,對信號進行端點的延拓,從而大大減弱了LMD算法計算過程中的端點效應(yīng),通過仿真信號和實驗軸承數(shù)據(jù)的分析,驗證了改進方法的有效性。

    1LMD原理

    LMD可以將原始信號自適應(yīng)地分解為一組頻率由高到低自動排列的乘積函數(shù)(product function, PF)分量,每個分量由一個純調(diào)頻信號和一個包絡(luò)信號相乘得到。LMD原理如下[12]。

    (1)確定原始信號x(t)所有的局部極值點ni,j,計算相鄰兩個極值點ni,j和ni,j+1的平均值mk,i,j:

    mk,i,j=(ni,j+ni,j+1)/2

    (1)

    將所有平均值mk,i,j用直線連接起來,然后利用滑動平均方法進行平滑處理,得到局部均值函數(shù)mk,i(t)。計算局部包絡(luò)估計值:

    ak,i,j=|ni,j-ni,j+1|/2

    (2)

    將所有包絡(luò)估計值ak,i,j用直線連接起來,同樣利用滑動平均方法進行平滑處理,得到包絡(luò)函數(shù)ak,i(t)。其中,k為PF分量的個數(shù),k=1,2,…,l;i為求取PF分量過程中的迭代次數(shù),i=1,2,…,n;j為極值點下標(biāo)。

    (2)將局部均值函數(shù)mk,i(t)從原始信號x(t)中分離出來:

    hk,i(t)=x(t)-mk,i(t)

    (3)

    (3)用hk,i(t)除以包絡(luò)估計函數(shù)ak,i(t),從而對hk,i(t)進行解調(diào),得

    sk,i(t)=hk,i(t)/ak,i(t)

    (4)

    理想地,sk,i(t)是一個純調(diào)頻函數(shù),其包絡(luò)值為1。如果sk,i(t)不滿足此條件,則再將其作為原始時間序列,重復(fù)步驟(1)~(3),直到某一調(diào)頻函數(shù)sk(t)的包絡(luò)值為1時停止。但是在實際應(yīng)用中,由于種種原因,純調(diào)頻信號的包絡(luò)值不能精確地為1,可設(shè)置一個閾值Δ,當(dāng)1-Δ≤a1n(t)≤1+Δ時,結(jié)束迭代。迭代得到的純調(diào)頻信號sk(t)求反余弦即可得到分量瞬時相位角:

    φk=arccos(sk(t))/(2π)

    (5)

    由瞬時相位角對時間求導(dǎo)可得到瞬時頻率ωk:

    ωk=dφk/dt

    (6)

    (4)將上述迭代過程中得到的所有包絡(luò)估計函數(shù)相乘,從而得到包絡(luò)信號ak(t):

    (7)

    (5)將包絡(luò)估計函數(shù)ak(t)與純調(diào)頻函數(shù)sk(t)相乘,得到原始信號的第k個攜帶有用信息的分量CPFk(t):

    CPFk(t)=ak(t)sk(t)

    (8)

    (6) 將第一個分量CPF1(t)從原始時間序列x(t)中分離出來后,得到一個新的信號u1(t),再將u1(t)視為原始信號,重復(fù)步驟(1)~(5),循環(huán)l次,直到ul為一個單調(diào)函數(shù)或常數(shù)函數(shù)為止。最終可以得到一組攜帶有用信息的分量和一個殘余分量ul,即

    (9)

    循環(huán)結(jié)束后得到l個PF分量及其對應(yīng)的瞬時頻率,通過時頻圖能很容易地獲得各分量頻率隨時間變化的軌跡信息,從而從頻域?qū)π盘柛鹘M分進行更準(zhǔn)確的分析。本文涉及的LMD方法的迭代停止準(zhǔn)則都應(yīng)用了文獻[12]提出的停止準(zhǔn)則,減少了迭代過程中的計算量。

    2LMD的改進

    由于離散時間序列求取包絡(luò)估計值和平均值后,序列長度較原始時間序列會變短,故需要在端點延拓值才能進行后續(xù)的計算。但是所延拓點的值是一個未知數(shù),無論采用何種延拓方法,分解得到的PF分量較信號原始組分都會有一定誤差。因此,可設(shè)法延拓序列使其更接近信號發(fā)展趨勢,盡可能地減弱端點效應(yīng),從而使分析更真實可靠。

    表1 基于序列差分特性的端點處理方法

    (a)(t)>0,(t)>0,xe<0

    (b)(t)>0,(t)>0,xe>0

    (c)(t)>0,(t)<0,xe<0

    圖1 信號右端曲線函數(shù)特征及延拓示意圖

    x(s)-x(s-1)

    (10)

    其中,s∈t且s=k,k+1,…,e,k是離端點最近的拐點的下標(biāo)。函數(shù)x(t)外延拓點幅值為xe+Δx1,這樣可以避免如鏡像延拓方法中延拓至該點幅值為Ma而造成幅值過大的問題;如圖1d所示,若xe>0時延拓至xe+Δx1也是為了防止幅值過大。其他情況可作同理分析。

    3仿真分析

    為驗證新的極值點延拓方法的有效性,對仿真信號

    x1=10sin(10πt)+30cos(100πt)+20cos(50πt)

    進行分析,采樣頻率為1 kHz,采樣點數(shù)為1000。對比LMD方法改進前后的分解結(jié)果如圖2和圖3所示,圖中,CPF1、CPF2、CPF3、CPF4及C分別為各PF分量及殘余分量幅值。圖2所示為原始LMD方法的分解結(jié)果,原LMD方法單純將信號端點值視為極值點進行計算,將仿真信號分解為4個分量和1個殘余分量,第4分量為因誤差產(chǎn)生的虛假分量。從圖2中可以很容易看出在分解到第3分量時,組分兩端出現(xiàn)了嚴(yán)重偏差,并影響到信號內(nèi)部區(qū)域。

    圖2 原始LMD方法分解結(jié)果一

    圖3 改進LMD方法分解結(jié)果一

    圖3所示為改進的LMD分解結(jié)果,包含3個組分和1個殘余分量,符合原始信號有3個不同頻率組分的事實,且到第3分量時端點效應(yīng)明顯減弱。端點效應(yīng)累積到第3分量時使信號產(chǎn)生較大的偏差,導(dǎo)致第3分量波形失真,第3分量分解結(jié)果與原組分的絕對誤差如圖4所示,從而驗證了改進方法的有效性。

    圖4 第3分量絕對誤差對比一

    用同樣的采樣頻率和采樣點數(shù)對仿真信號

    x2=10sin(10πt)+30cos(120πt)+20sin(50πt)

    進行分解得到圖5和圖6所示結(jié)果。對比圖5和圖6結(jié)果中的前3個分量可以發(fā)現(xiàn),原始LMD結(jié)果在分解到第3分量時產(chǎn)生嚴(yán)重的端點效應(yīng),并且影響到了信號內(nèi)部,而圖6中的端點效應(yīng)則有了明顯的改善。

    圖5 原始LMD方法分解結(jié)果二

    圖6 改進LMD方法分解結(jié)果二

    圖7所示為LMD方法改進前后,分解到第3分量時產(chǎn)生的與原分量的絕對誤差對比。從圖7中可以很容易看出改進前產(chǎn)生了很大的誤差,并且嚴(yán)重影響到信號內(nèi)部;而應(yīng)用改進方法后,信號的分解結(jié)果有了很大的改善,使分解結(jié)果更可靠。

    圖7 第3分量絕對誤差對比二

    4實驗信號分析

    實驗信號的驗證選取美國凱斯西儲大學(xué)軸承數(shù)據(jù)中心公開的電機驅(qū)動端SKF6205-2RS深溝球軸承外圈和內(nèi)圈故障數(shù)據(jù),軸承參數(shù)見表2,故障特征頻率分別為108 Hz和162 Hz。采樣頻率為12 kHz,分析點數(shù)為6000。實驗信號的分解采用了鏡像延拓[13]改進的LMD方法和本文改進的LMD方法兩種方法進行對比。

    表2 SKF6205-2RS軸承結(jié)構(gòu)尺寸參數(shù) mm

    對于軸承外圈故障信號,鏡像延拓LMD分解結(jié)果如圖8所示,得到2個分量;新方法改進的LMD分解如圖9所示,得到3個分量。圖8中第1分量對應(yīng)的頻譜中不能找到外圈故障特征頻率108 Hz,該分量含有較寬頻帶的噪聲成分和因故障產(chǎn)生的共振成分,第2分量對應(yīng)的頻譜中可以找到故障特征頻率,但其仍含有大量的噪聲成分。圖9中第1分量對應(yīng)頻譜中噪聲頻率能量明顯提高,并且改進后的算法在第2和第3分量中都能找到外圈故障特征頻率,說明改進的LMD方法對各頻率段的成分分解得更徹底,使故障分量信號特征更明顯,驗證了改進算法的優(yōu)越性。

    圖8 外圈信號鏡像延拓LMD分解及頻譜

    圖9 外圈信號新改進LMD分解及頻譜

    對于軸承內(nèi)圈故障信號,鏡像延拓LMD分解結(jié)果如圖10所示,分解得到3個分量;新方法改進的LMD分解得到4個分量,如圖11所示。圖10、圖11在第2和第3分量對應(yīng)的頻譜中都能找到內(nèi)圈故障特征頻率162 Hz。從頻譜來看,圖10中第3分量包含的低頻成分比圖11中第3分量的低頻成分要多,而圖11分出4個分量,將圖10中第3分量進行進一步分解,使不同頻率段的信號得到了更為徹底的分解,由此驗證了新方法改進LMD的有效性。

    圖10 內(nèi)圈信號鏡像延拓LMD分解及頻譜

    圖11 內(nèi)圈信號新改進LMD分解及頻譜

    5結(jié)語

    本文針對LMD方法中存在的端點效應(yīng)問題,結(jié)合曲線函數(shù)的一階導(dǎo)數(shù)和二階導(dǎo)數(shù)判斷函數(shù)趨勢的原理,提出了基于信號序列差分的端點延拓方法。該方法兼顧了延拓點的誤差和算法的復(fù)雜程度兩個方面,在信號端點延拓一個盡可能接近下一個極值點的點以減弱端點效應(yīng)。利用該方法對仿真信號和實驗軸承外圈和內(nèi)圈故障信號進行了分析,驗證了該方法的實用性。

    參考文獻:

    [1]Lee J Y. Variable Short-time Fourier Transform for Vibration Signals with Transients[J]. Journal of Vibration and Control, 2015, 21(7): 1383-1397.

    [2]申永軍, 張光明, 楊紹普, 等. 基于Gabor變換的欠定盲信號分離新方法[J]. 振動、測試與診斷, 2011, 31(3): 309-313.

    Shen Yongjun, Zhang Guangming, Yang Shaopu, et al. New Method for Blind Source Separation in Under-determined Mixtures Based on Gabor Transform[J]. Journal of Vibration, Measurement & Diagnosis, 2011, 31(3): 309-313.

    [3]石林鎖, 張亞洲, 米文鵬. 基于WVD的譜峭度法在軸承故障診斷中的應(yīng)用[J]. 振動、測試與診斷, 2011, 31(1): 27-31.

    Shi Linsuo, Zhang Yazhou, Mi Wenpeng. Application of Wigner-Ville-distribution Based on Spectral Kurtosis Algorithm to Fault Diagnosis of Rolling Bearing[J].Journal of Vibration, Measurement & Diagnosis, 2011, 31(1): 27-31.

    [4]向東陽, 吳正國, 胡文彪, 等. 改進的小波變換系數(shù)相關(guān)去噪方法[J]. 振動、測試與診斷, 2010, 30(5): 561-565.

    Xiang Dongyang, Wu Zhengguo, Hu Wenbiao, et al. Improved Denoising Algorithm Using Correlation of Wavelet Coefficient[J]. Journal of Vibration, Measurement & Diagnosis, 2010, 30(5): 561-565.

    [5]吳一全, 陶飛翔, 曹照清. 利用雙樹復(fù)小波變換和SURF的圖像配準(zhǔn)算法[J].系統(tǒng)工程與電子技術(shù), 2014, 36(5): 997-1003.

    Wu Yiquan, Tao Feixiang, Cao Zhaoqing. Image Registration Algorithm Based on Dual Tree Complex Wavelet Transform and SURF[J].Systems Engineering and Electronics, 2014, 36(5): 997-1003.

    [6]Huang N E, Shen Z, Long S R,et al. The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Non-stationary Time Series Analysis[J].Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 1998, 454(1971): 903-995.

    [7]Smith J S. The Local Mean Decomposition and Its Application to EEG Perception Data[J]. Journal of the Royal Society Interface, 2005, 2(5): 443-454.

    [8]Yang Yu, Cheng Junsheng, Zhang Kang. An Ensemble Local Means Decomposition Method and Its Application to Local Rub-impact Fault Diagnosis of the Rotor Systems[J]. Measurement, 2012, 45(3): 561-570.

    [9]Park C, Looney D, Hulle M M V, et al. The Complex Local Mean Decomposition[J].Neurocomputing, 2011, 74(6): 867-875.

    [10]任達千, 楊世錫, 吳昭同, 等. LMD時頻分析方法的端點效應(yīng)在旋轉(zhuǎn)機械故障診斷中的影響[J]. 中國機械工程, 2012, 23(8): 951-956.

    Ren Daqian, Yang Shixi, Wu Zhaotong, et al. Research on End Effect of LMD Based Time-frequency Analysis in Rotating Machinery Fault Diagnosis[J].China Mechanical Engineering, 2012, 23(8): 951-956.

    [11]孟宗, 李珊珊, 季艷. 基于對稱差分能量算子解調(diào)的局部均值分解端點效應(yīng)抑制方法[J].機械工程學(xué)報, 2014, 50(13): 80-87.

    Meng Zong, Li Shanshan, Ji Yan. Restraining Method for End Effect of Local Mean Decomposition Based on Energy Operator Demodulation of Symmetrical Differencing[J]. Journal of Mechanical Engineering, 2014, 50(13): 80-87.

    [12]唐貴基, 王曉龍. 基于局部均值分解和切片雙譜的滾動軸承故障診斷研究[J].振動與沖擊, 2013, 32(24): 83-88.

    Tang Guiji, Wang Xiaolong. Fault Diagnosis of Roller Bearings Based on Local Mean Decomposition and Slice Bi-spectrum[J]. Journal of Vibration and Shock, 2013, 32(24): 83-88.

    [13]韓建平, 錢炯, 董小軍. 采用鏡像延拓和RBF神經(jīng)網(wǎng)絡(luò)處理EMD中端點效應(yīng)[J].振動、測試與診斷, 2010, 30(4):414-417.

    Han Jianping, Qian Jiong, Dong Xiaojun. Suppression of End-effect in Empirical Mode Decomposition by Mirror Extension and Radial Basis Function Neural Network Prediction[J]. Journal of Vibration, Measurement and Diagnosis, 2010, 30(4):414-417.

    (編輯陳勇)

    A Processing Method to Weaken End Effects of LMD Based on Difference Performance of Time Series

    Du DongmeiZhang ZhaoHe QingLi Hong

    North China Electric Power University,Beijing,102206

    Abstract:End effects occured at the ends of signals after calculating their envelope and mean value estimation with LMD algorithm. In order to weaken the influences, a continuation approach was proposed based on the difference performance of time series. It considered the signals as a discrete function, and judged the signal tendency by the signs of their first-order difference and second-order difference, determining the value of extreme points of continuation according to the specific situation.This new method was used to decompose the synthetic signals and test signals of rolling bearings with outer raceway faults and inner raceway faults, their results were compared with those decomposed by the LMD improved by mirror extension. Its good performance in decomposition validates the effectiveness of this new approach in signal analysis.

    Key words:local mean decomposition(LMD); end effect; difference of time series; continuation

    收稿日期:2015-06-25

    基金項目:中央高校基本科研業(yè)務(wù)費專項資金資助項目(2014MS17)

    中圖分類號:TH133.3; TH17

    DOI:10.3969/j.issn.1004-132X.2016.10.007

    作者簡介:杜冬梅,女,1964年生。華北電力大學(xué)能源動力與機械工程學(xué)院教授。主要研究方向為狀態(tài)監(jiān)測與故障診斷。張昭,男,1990年生。華北電力大學(xué)能源動力與機械工程學(xué)院碩士研究生。何青,男,1962年生。華北電力大學(xué)能源動力與機械工程學(xué)院教授。李紅,女,1977年生。華北電力大學(xué)能源動力與機械工程學(xué)院講師。

    国产亚洲5aaaaa淫片| 日韩人妻高清精品专区| 六月丁香七月| 又爽又黄无遮挡网站| 欧美丝袜亚洲另类| 欧美成人a在线观看| 99久久无色码亚洲精品果冻| 丰满乱子伦码专区| 久久久国产成人免费| 国产精品国产三级国产专区5o | 麻豆一二三区av精品| 女人久久www免费人成看片 | 老女人水多毛片| 波多野结衣巨乳人妻| 精品国产露脸久久av麻豆 | 久久这里只有精品中国| 毛片一级片免费看久久久久| 亚洲av电影不卡..在线观看| 国产成人freesex在线| 精品国产露脸久久av麻豆 | 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区在线av高清观看| av卡一久久| 人妻少妇偷人精品九色| 国产欧美日韩精品一区二区| 边亲边吃奶的免费视频| 日韩制服骚丝袜av| 亚洲av熟女| 午夜激情福利司机影院| 欧美变态另类bdsm刘玥| 日韩av不卡免费在线播放| 亚洲av免费在线观看| 超碰av人人做人人爽久久| 男人和女人高潮做爰伦理| 亚洲av电影不卡..在线观看| 免费无遮挡裸体视频| 亚洲av中文av极速乱| 丰满人妻一区二区三区视频av| 久久99热这里只频精品6学生 | 天堂中文最新版在线下载 | 人人妻人人澡人人爽人人夜夜 | 久久久久久久久大av| 亚洲av熟女| 免费看日本二区| 中文字幕av在线有码专区| av福利片在线观看| 在线天堂最新版资源| 91久久精品电影网| av播播在线观看一区| 精华霜和精华液先用哪个| 日本色播在线视频| 国产一区二区三区av在线| 亚洲精品乱码久久久久久按摩| 黄色日韩在线| 成人亚洲欧美一区二区av| 色5月婷婷丁香| 亚洲久久久久久中文字幕| 国产老妇伦熟女老妇高清| 久久精品久久久久久噜噜老黄 | 国产高清有码在线观看视频| 久久久久久久午夜电影| 高清日韩中文字幕在线| 99视频精品全部免费 在线| 国产又色又爽无遮挡免| 国产亚洲一区二区精品| 免费搜索国产男女视频| 免费搜索国产男女视频| 午夜免费激情av| 一级毛片我不卡| 日韩成人av中文字幕在线观看| 丰满乱子伦码专区| 亚洲伊人久久精品综合 | 久久99精品国语久久久| 国产精品一区二区三区四区免费观看| 亚洲精品一区蜜桃| 乱码一卡2卡4卡精品| 狂野欧美白嫩少妇大欣赏| 久久久亚洲精品成人影院| 欧美一区二区精品小视频在线| 久久久午夜欧美精品| 特级一级黄色大片| 美女国产视频在线观看| 波多野结衣高清无吗| 黄色一级大片看看| 国产精品一区二区三区四区免费观看| 亚洲精品亚洲一区二区| 99久久中文字幕三级久久日本| 一级毛片电影观看 | 91av网一区二区| 欧美日本亚洲视频在线播放| 久久久午夜欧美精品| 国产免费一级a男人的天堂| 精品少妇黑人巨大在线播放 | 国产 一区精品| 午夜福利在线观看免费完整高清在| 国产精品久久久久久久电影| 美女脱内裤让男人舔精品视频| 看黄色毛片网站| 蜜桃亚洲精品一区二区三区| 99久久精品国产国产毛片| 亚洲自偷自拍三级| 又粗又硬又长又爽又黄的视频| 男女下面进入的视频免费午夜| 十分钟在线观看高清视频www| 亚洲天堂av无毛| 国产欧美日韩综合在线一区二区| 街头女战士在线观看网站| 亚洲精品,欧美精品| 亚洲av欧美aⅴ国产| 亚洲精品美女久久av网站| 精品99又大又爽又粗少妇毛片| 国产综合精华液| 成年人午夜在线观看视频| 在线观看美女被高潮喷水网站| 一二三四中文在线观看免费高清| 飞空精品影院首页| 老司机影院毛片| 亚洲国产av影院在线观看| 好男人视频免费观看在线| 亚洲av在线观看美女高潮| 中文字幕另类日韩欧美亚洲嫩草| 亚洲一级一片aⅴ在线观看| 9色porny在线观看| 日本av手机在线免费观看| 两个人看的免费小视频| 男女边吃奶边做爰视频| 久久久国产欧美日韩av| 亚洲,欧美,日韩| 最近最新中文字幕大全免费视频 | 99久久中文字幕三级久久日本| av福利片在线| 亚洲精品aⅴ在线观看| av.在线天堂| 最近中文字幕2019免费版| av.在线天堂| 国产有黄有色有爽视频| 国产精品久久久久久久电影| 亚洲欧美精品自产自拍| 免费人妻精品一区二区三区视频| 久久精品国产自在天天线| 亚洲熟女精品中文字幕| 一本色道久久久久久精品综合| 观看美女的网站| 丰满少妇做爰视频| 99热6这里只有精品| 999精品在线视频| 曰老女人黄片| 一本大道久久a久久精品| 十分钟在线观看高清视频www| 制服丝袜香蕉在线| 波野结衣二区三区在线| 亚洲一码二码三码区别大吗| av.在线天堂| 男女边摸边吃奶| 欧美bdsm另类| 亚洲成国产人片在线观看| 五月玫瑰六月丁香| av免费在线看不卡| 黑人高潮一二区| 黄色 视频免费看| 啦啦啦啦在线视频资源| 日韩制服丝袜自拍偷拍| 日韩欧美一区视频在线观看| 性高湖久久久久久久久免费观看| 大香蕉久久成人网| 九色成人免费人妻av| 免费观看无遮挡的男女| 美国免费a级毛片| 欧美 日韩 精品 国产| 午夜福利乱码中文字幕| 你懂的网址亚洲精品在线观看| 午夜av观看不卡| 有码 亚洲区| 欧美精品av麻豆av| 少妇被粗大的猛进出69影院 | 高清在线视频一区二区三区| 一级爰片在线观看| 国产成人精品无人区| 精品久久国产蜜桃| 丰满饥渴人妻一区二区三| a级片在线免费高清观看视频| 国产黄色免费在线视频| xxxhd国产人妻xxx| 亚洲精品中文字幕在线视频| 国产成人午夜福利电影在线观看| 亚洲国产精品一区二区三区在线| 国国产精品蜜臀av免费| 下体分泌物呈黄色| 夜夜爽夜夜爽视频| 国产欧美亚洲国产| 久久99一区二区三区| 黑丝袜美女国产一区| 99精国产麻豆久久婷婷| 欧美精品一区二区大全| 亚洲av中文av极速乱| 欧美日韩精品成人综合77777| 免费观看在线日韩| 国产精品秋霞免费鲁丝片| videossex国产| 日韩制服丝袜自拍偷拍| 亚洲,欧美,日韩| 日本午夜av视频| 国产精品久久久久久av不卡| 美女国产高潮福利片在线看| 91在线精品国自产拍蜜月| 欧美成人午夜精品| 国产色婷婷99| 99国产综合亚洲精品| 精品一区二区免费观看| 在现免费观看毛片| 国产精品久久久久久久久免| 久久久久久久久久成人| 久久99一区二区三区| 七月丁香在线播放| 日本爱情动作片www.在线观看| 2021少妇久久久久久久久久久| 国产69精品久久久久777片| 日韩成人伦理影院| 91aial.com中文字幕在线观看| 制服人妻中文乱码| 999精品在线视频| 五月开心婷婷网| 美女脱内裤让男人舔精品视频| 国产一区二区在线观看日韩| 久久97久久精品| 国产乱来视频区| 亚洲精品一二三| 制服丝袜香蕉在线| 777米奇影视久久| 中国国产av一级| 99热这里只有是精品在线观看| 亚洲成av片中文字幕在线观看 | 国产一区有黄有色的免费视频| 亚洲av日韩在线播放| 亚洲精品一二三| 国产精品国产三级国产专区5o| 一区二区三区乱码不卡18| videossex国产| 午夜影院在线不卡| 一级,二级,三级黄色视频| 久久影院123| 伊人久久国产一区二区| 美女大奶头黄色视频| 男女边吃奶边做爰视频| 国产精品蜜桃在线观看| 亚洲精品美女久久av网站| 爱豆传媒免费全集在线观看| 中文乱码字字幕精品一区二区三区| 国产av一区二区精品久久| 大香蕉久久成人网| 综合色丁香网| 免费黄色在线免费观看| 国产在线一区二区三区精| 99热全是精品| 成人午夜精彩视频在线观看| 欧美精品一区二区大全| 国产精品一区二区在线不卡| 99热这里只有是精品在线观看| 考比视频在线观看| 日本av免费视频播放| 黄色毛片三级朝国网站| 18禁观看日本| 伦理电影大哥的女人| 国产熟女午夜一区二区三区| 青春草亚洲视频在线观看| 国产乱来视频区| 亚洲,欧美,日韩| 亚洲美女视频黄频| 黄色毛片三级朝国网站| 国产成人欧美| 国产欧美亚洲国产| 国产国语露脸激情在线看| 99九九在线精品视频| 高清毛片免费看| 在线天堂最新版资源| 男人操女人黄网站| 人人妻人人添人人爽欧美一区卜| 人妻系列 视频| 久久狼人影院| 交换朋友夫妻互换小说| 久久精品国产亚洲av天美| videos熟女内射| 久久亚洲国产成人精品v| 黄色配什么色好看| 久久久久久久久久久久大奶| 黄片无遮挡物在线观看| 精品一区二区免费观看| 大片免费播放器 马上看| 乱人伦中国视频| 国产精品不卡视频一区二区| 欧美精品人与动牲交sv欧美| 精品午夜福利在线看| 欧美 亚洲 国产 日韩一| 久久精品aⅴ一区二区三区四区 | 麻豆精品久久久久久蜜桃| 性高湖久久久久久久久免费观看| 极品人妻少妇av视频| 欧美精品亚洲一区二区| 亚洲精品美女久久av网站| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 亚洲欧美精品自产自拍| 亚洲精品国产色婷婷电影| 在线亚洲精品国产二区图片欧美| 80岁老熟妇乱子伦牲交| 侵犯人妻中文字幕一二三四区| 日韩欧美一区视频在线观看| 五月伊人婷婷丁香| 在线观看免费高清a一片| 亚洲欧美色中文字幕在线| 精品国产露脸久久av麻豆| 五月玫瑰六月丁香| 亚洲欧美一区二区三区国产| 99久久人妻综合| 90打野战视频偷拍视频| 制服诱惑二区| 国内精品宾馆在线| 丁香六月天网| 国产伦理片在线播放av一区| 日本av免费视频播放| 大香蕉97超碰在线| 久久久精品免费免费高清| 亚洲精品视频女| 国产极品天堂在线| 综合色丁香网| 国产无遮挡羞羞视频在线观看| 国产男女超爽视频在线观看| 国产男女内射视频| 国产无遮挡羞羞视频在线观看| 久久人妻熟女aⅴ| 欧美日韩精品成人综合77777| 久久鲁丝午夜福利片| 欧美最新免费一区二区三区| 波多野结衣一区麻豆| 日本av手机在线免费观看| 狠狠精品人妻久久久久久综合| 免费女性裸体啪啪无遮挡网站| 亚洲精品久久午夜乱码| 久久婷婷青草| 亚洲国产精品专区欧美| 成人二区视频| 美女大奶头黄色视频| 国产成人免费无遮挡视频| 日韩欧美一区视频在线观看| 久久久久久伊人网av| 色婷婷久久久亚洲欧美| 久久精品人人爽人人爽视色| 国产男女内射视频| 十八禁高潮呻吟视频| 国产黄频视频在线观看| 丰满少妇做爰视频| 国产不卡av网站在线观看| 黄片无遮挡物在线观看| 成人漫画全彩无遮挡| 精品亚洲乱码少妇综合久久| 久久久久久久久久人人人人人人| 国产精品一国产av| 色视频在线一区二区三区| tube8黄色片| 国产永久视频网站| 精品久久蜜臀av无| 91精品国产国语对白视频| 在线观看人妻少妇| 天堂中文最新版在线下载| www.色视频.com| av国产久精品久网站免费入址| 高清毛片免费看| 老司机影院毛片| 一级片'在线观看视频| 国产乱来视频区| 日本黄色日本黄色录像| 成人黄色视频免费在线看| 看免费av毛片| 自线自在国产av| 男的添女的下面高潮视频| 校园人妻丝袜中文字幕| 十八禁高潮呻吟视频| 22中文网久久字幕| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 九色亚洲精品在线播放| 伊人久久国产一区二区| 亚洲国产精品一区二区三区在线| 亚洲国产精品专区欧美| 蜜桃国产av成人99| 男女下面插进去视频免费观看 | 成年动漫av网址| 黄片无遮挡物在线观看| 亚洲婷婷狠狠爱综合网| 久久狼人影院| 一区二区av电影网| 美女视频免费永久观看网站| 久久久国产精品麻豆| 又粗又硬又长又爽又黄的视频| 啦啦啦啦在线视频资源| 成人毛片60女人毛片免费| 高清不卡的av网站| 丝袜喷水一区| 丝袜人妻中文字幕| av国产久精品久网站免费入址| 国产一区有黄有色的免费视频| 亚洲精品乱码久久久久久按摩| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美日韩在线播放| 人体艺术视频欧美日本| 成人国产av品久久久| 久久狼人影院| 在线观看免费日韩欧美大片| 亚洲欧美日韩卡通动漫| 啦啦啦啦在线视频资源| 国产男女超爽视频在线观看| 又粗又硬又长又爽又黄的视频| 另类精品久久| 涩涩av久久男人的天堂| 欧美人与性动交α欧美精品济南到 | 日日摸夜夜添夜夜爱| 国产精品不卡视频一区二区| 国产男女内射视频| 久久精品国产a三级三级三级| 亚洲欧美一区二区三区黑人 | 亚洲高清免费不卡视频| 熟女av电影| 色哟哟·www| 老女人水多毛片| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 久久久久国产网址| 亚洲精品一二三| 五月开心婷婷网| 久久ye,这里只有精品| av黄色大香蕉| 中文欧美无线码| 午夜福利影视在线免费观看| 精品少妇黑人巨大在线播放| 久久人人爽av亚洲精品天堂| 蜜桃在线观看..| 欧美精品国产亚洲| 精品国产一区二区三区四区第35| 色吧在线观看| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 亚洲精品第二区| 国产免费一级a男人的天堂| 美女福利国产在线| 久久久久国产网址| xxxhd国产人妻xxx| 精品少妇黑人巨大在线播放| 黄片播放在线免费| 午夜福利乱码中文字幕| 国产精品人妻久久久影院| 蜜桃国产av成人99| 亚洲精品久久久久久婷婷小说| 亚洲欧美中文字幕日韩二区| 制服丝袜香蕉在线| 一区二区三区四区激情视频| 99久久精品国产国产毛片| 黑人欧美特级aaaaaa片| 97在线人人人人妻| 成人黄色视频免费在线看| 狂野欧美激情性bbbbbb| 亚洲av电影在线进入| 国产精品三级大全| 乱人伦中国视频| 久久99精品国语久久久| 99九九在线精品视频| 精品人妻熟女毛片av久久网站| 亚洲综合色网址| 午夜福利视频在线观看免费| 99热国产这里只有精品6| 人妻系列 视频| 狠狠婷婷综合久久久久久88av| 日韩成人av中文字幕在线观看| 久久精品熟女亚洲av麻豆精品| 免费播放大片免费观看视频在线观看| 香蕉精品网在线| 欧美日韩成人在线一区二区| 国产精品国产三级国产专区5o| 亚洲精品日韩在线中文字幕| 亚洲欧美日韩另类电影网站| 一级,二级,三级黄色视频| 2021少妇久久久久久久久久久| 亚洲国产日韩一区二区| 色哟哟·www| 成年女人在线观看亚洲视频| 18禁裸乳无遮挡动漫免费视频| 欧美日韩国产mv在线观看视频| 大香蕉久久网| 久久久久人妻精品一区果冻| 亚洲丝袜综合中文字幕| 午夜激情久久久久久久| 男人添女人高潮全过程视频| 日韩三级伦理在线观看| 亚洲精品国产av蜜桃| 欧美日韩精品成人综合77777| 国产精品一二三区在线看| 国产女主播在线喷水免费视频网站| 好男人视频免费观看在线| 亚洲精品美女久久久久99蜜臀 | 有码 亚洲区| 久久精品aⅴ一区二区三区四区 | 亚洲经典国产精华液单| 欧美激情国产日韩精品一区| 免费人妻精品一区二区三区视频| 人妻少妇偷人精品九色| 日韩av不卡免费在线播放| 老司机亚洲免费影院| 国产男女超爽视频在线观看| 中文字幕av电影在线播放| 亚洲久久久国产精品| 亚洲情色 制服丝袜| 久久久欧美国产精品| 高清av免费在线| 日韩精品免费视频一区二区三区 | 国产视频首页在线观看| 久久精品熟女亚洲av麻豆精品| 大码成人一级视频| 乱码一卡2卡4卡精品| 国产精品成人在线| a级毛片黄视频| 亚洲欧美精品自产自拍| 夜夜骑夜夜射夜夜干| 90打野战视频偷拍视频| 毛片一级片免费看久久久久| 成人亚洲精品一区在线观看| av国产久精品久网站免费入址| 亚洲人成77777在线视频| av国产久精品久网站免费入址| 性高湖久久久久久久久免费观看| videossex国产| 狠狠婷婷综合久久久久久88av| av不卡在线播放| 免费女性裸体啪啪无遮挡网站| 中文乱码字字幕精品一区二区三区| 亚洲精品aⅴ在线观看| 亚洲精品第二区| 91国产中文字幕| 欧美日韩亚洲高清精品| 亚洲欧美成人综合另类久久久| 99九九在线精品视频| 亚洲高清免费不卡视频| 丁香六月天网| 七月丁香在线播放| 精品人妻一区二区三区麻豆| 日本av免费视频播放| 欧美性感艳星| 精品一品国产午夜福利视频| 欧美国产精品一级二级三级| 亚洲精品美女久久久久99蜜臀 | 自线自在国产av| 精品久久蜜臀av无| av视频免费观看在线观看| 97在线视频观看| 久久av网站| 人人澡人人妻人| 久久精品国产亚洲av涩爱| 视频中文字幕在线观看| 成年人免费黄色播放视频| 午夜免费鲁丝| 亚洲成人一二三区av| 蜜臀久久99精品久久宅男| 丝袜在线中文字幕| 亚洲国产av新网站| 日日啪夜夜爽| 日本猛色少妇xxxxx猛交久久| 亚洲性久久影院| 欧美日韩视频高清一区二区三区二| 免费大片黄手机在线观看| 国语对白做爰xxxⅹ性视频网站| 久久久国产精品麻豆| 国产精品嫩草影院av在线观看| 国产免费一级a男人的天堂| 九色亚洲精品在线播放| 汤姆久久久久久久影院中文字幕| 另类精品久久| 亚洲美女视频黄频| 国产一区有黄有色的免费视频| 2018国产大陆天天弄谢| 少妇被粗大的猛进出69影院 | 免费黄色在线免费观看| 街头女战士在线观看网站| 久久久久久人人人人人| 亚洲国产av影院在线观看| av在线app专区| 大片电影免费在线观看免费| 成人午夜精彩视频在线观看| 制服人妻中文乱码| 国产成人精品福利久久| 一区二区三区乱码不卡18| 成人毛片a级毛片在线播放| 黄网站色视频无遮挡免费观看| 如日韩欧美国产精品一区二区三区| 国产精品蜜桃在线观看| 亚洲一级一片aⅴ在线观看| 一区二区日韩欧美中文字幕 | 中文字幕亚洲精品专区| 久久精品久久久久久噜噜老黄| 欧美精品亚洲一区二区| 激情视频va一区二区三区| 亚洲av国产av综合av卡| 国产成人欧美| 一区二区av电影网| 欧美精品高潮呻吟av久久| 一级a做视频免费观看| 日韩精品免费视频一区二区三区 | 国产精品一二三区在线看| 亚洲伊人久久精品综合| 国产免费又黄又爽又色| 久久久久人妻精品一区果冻| 精品人妻熟女毛片av久久网站| 国产亚洲精品久久久com| 日韩一本色道免费dvd| 在线观看www视频免费| 日韩制服骚丝袜av| 国产又爽黄色视频| 丰满乱子伦码专区|