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

    基于改進(jìn)經(jīng)驗(yàn)小波變換的機(jī)車(chē)軸承故障診斷

    2019-04-09 11:29:14段晨東
    中國(guó)機(jī)械工程 2019年6期
    關(guān)鍵詞:峭度頻帶時(shí)頻

    段晨東 張 榮

    長(zhǎng)安大學(xué)電子與控制工程學(xué)院,西安,710064

    0 引言

    機(jī)車(chē)軸承是鐵路機(jī)車(chē)行走部的關(guān)鍵部件,在重載工況下工作,振動(dòng)信號(hào)噪聲大, 其損傷特征提取難[1]。王曉冬等[2]針對(duì)機(jī)車(chē)軸承信號(hào)特征,構(gòu)造自適應(yīng)多小波,以峭度最大為目標(biāo)函數(shù),采用基于遺傳算法的多小波譜峭度來(lái)確定最佳頻帶,有效地提取了軸承信號(hào)的特征分量。LEI等[3]提出了基于小波包變換的譜峭度方法,以小波包分量的峭度最大為指標(biāo),從多層小波包分量中尋找最佳小波包分量,以此分離出了機(jī)車(chē)包含軸承損傷特征的有效分量。ZHAO等[4]提出了基于集成經(jīng)驗(yàn)?zāi)B(tài)分解的最佳模式分量的包絡(luò)階次跟蹤分析方法,在軸承轉(zhuǎn)速不恒定時(shí),采用該方法可以撲捉軸承故障的特征。LI等[5]提出了一種基于相關(guān)分析的獨(dú)立變分模態(tài)分解方法,以提取有效分量用于包絡(luò)分析,發(fā)現(xiàn)了機(jī)車(chē)軸承的微弱故障和混合故障。上述這些方法中,小波變換方法對(duì)頻帶的劃分無(wú)法做到自適應(yīng),經(jīng)驗(yàn)?zāi)B(tài)分解方法缺乏完備的理論基礎(chǔ),因此,采用高效的信號(hào)處理方法是提取軸承故障特征的關(guān)鍵。經(jīng)驗(yàn)小波變換(empirical wavelet transform, EWT)[6]結(jié)合了小波分析的完備理論性和經(jīng)驗(yàn)?zāi)J椒纸獾淖赃m應(yīng)性,通過(guò)信號(hào)的頻譜局部極大值的頻率自適應(yīng)劃分其分析頻帶,并基于劃分的子頻帶構(gòu)造一系列合適的正交小波濾波器組對(duì)信號(hào)進(jìn)行分解。然而,信噪比較低的工程信號(hào)的噪聲頻譜較寬,頻譜上噪聲分量的譜峰會(huì)影響信號(hào)的峰值分布,干擾子頻帶劃分的合理性,使EWT提取特征的有效性受到影響。為了克服EWT基于頻譜局部極大值劃分頻帶方法的不足,針對(duì)機(jī)車(chē)軸承的故障特征提取,筆者提出一種采用信號(hào)時(shí)頻峭度譜局部極小值劃分頻帶的方法,在時(shí)頻變換的基礎(chǔ)上求取信號(hào)的時(shí)頻峭度譜,由其峭度譜局部極小值對(duì)應(yīng)的頻率來(lái)確定子頻帶的邊界,構(gòu)造正交小波濾波器組對(duì)信號(hào)進(jìn)行EWT。仿真實(shí)驗(yàn)和工程應(yīng)用表明,改進(jìn)后的EWT能夠克服噪聲的干擾,有效地提取軸承的損傷故障特征。

    1 基于時(shí)頻峭度譜的經(jīng)驗(yàn)小波變換

    1.1 經(jīng)驗(yàn)小波變換(EWT)

    經(jīng)驗(yàn)小波變換是經(jīng)驗(yàn)?zāi)B(tài)分解原理與小波框架理論相結(jié)合的一種信號(hào)分析方法[6],即把任意的一個(gè)時(shí)域信號(hào)看成是一組調(diào)頻調(diào)幅信號(hào)分量的組合,應(yīng)用小波框架理論,基于信號(hào)的頻域分布特征構(gòu)造一系列小波濾波器組,以此提取不同頻帶的調(diào)頻調(diào)幅分量。經(jīng)驗(yàn)小波變換原理如下:

    (1)對(duì)信號(hào)進(jìn)行傅里葉變換,求取其頻譜。設(shè)信號(hào)f(t)的傅里葉變換為

    (1)

    (2)根據(jù)信號(hào)的頻譜特征,分割頻帶。在頻譜f(ω)上尋找M個(gè)局部極大值,并按降序排列,然后將它們對(duì)應(yīng)的頻率按升序排列后記為ωm(m=1,2,…,M)。將分析頻帶[0, π]劃分成N(N≤M)個(gè)子頻帶。設(shè)以ωn為中心的調(diào)幅調(diào)頻分量所在的子頻帶為[Ωn-1,Ωn][7],其中,Ωn-1=(ωn-1+ωn)/2,n=2,3,…,N。令Ω0=0,ΩN=π,這樣可把信號(hào)的分析頻帶 [0, π]劃分成N個(gè)連續(xù)的、互不交疊的子頻帶。

    (2)

    (3)

    對(duì)于信號(hào)f(t),在N個(gè)子頻帶上,利用式(2)、式(3)可以設(shè)計(jì)出1個(gè)低通濾波器(尺度函數(shù))和N-1個(gè)帶通濾波器(小波函數(shù))組成的自適應(yīng)小波濾波器組。

    (4)

    逼近系數(shù)為

    (5)

    *φ1(t)

    (6)

    (7)

    其中,*表示函數(shù)的卷積運(yùn)算。

    經(jīng)驗(yàn)小波變換是完全重構(gòu)的,重構(gòu)表達(dá)式為[9]

    (8)

    設(shè)試驗(yàn)信號(hào)為

    y(t)=sin(240πt)+sin(80πt)+cos(20πt)

    (9)

    圖1為信號(hào)時(shí)域波形和頻譜,采樣頻率為2 000 Hz,數(shù)據(jù)長(zhǎng)度為1 024,圖1b中的短虛線為采用上述方法得到的頻帶劃分邊界線。加入白噪聲的信號(hào)y(t)的時(shí)域波形如圖1c所示,信噪比為-1.5 dB,圖1d為其頻帶分割圖,可以看到,噪聲分量產(chǎn)生的譜峰(局部極大值)影響了頻帶分割的有效性,信號(hào)y(t)的3個(gè)分量的峰值并沒(méi)有處在分割子頻帶的中心位置。當(dāng)信號(hào)被噪聲污染時(shí),噪聲信號(hào)分量的無(wú)效局部峰值會(huì)影響經(jīng)驗(yàn)小波變換的頻帶劃分的合理性[9],使得EWT分解得到的信號(hào)分量的有效性降低。

    1.2 基于時(shí)頻峭度譜的經(jīng)驗(yàn)小波變換改進(jìn)算法

    (a)無(wú)噪聲時(shí)的波形圖

    (b)無(wú)噪聲時(shí)的頻譜分割圖

    (c)有噪聲時(shí)的波形圖

    (d)有噪聲時(shí)的頻譜分割圖圖1 信號(hào)的波形和頻譜分割Fig.1 Waveform and spectral segmentation of a signal

    為了克服信號(hào)噪聲分量的譜峰對(duì)頻帶劃分的影響,本文用一種信號(hào)時(shí)頻峭度譜代替信號(hào)的頻譜,通過(guò)峭度峰值對(duì)應(yīng)的頻率來(lái)確定子頻帶的邊界,以此為基礎(chǔ)構(gòu)造小波濾波器組,提取信號(hào)中的幅值調(diào)制分量。改進(jìn)的經(jīng)驗(yàn)小波變換原理如下。

    (10)

    (11)

    (2)以時(shí)頻分解結(jié)果為基礎(chǔ),求信號(hào)的時(shí)頻峭度譜。在信號(hào)f(t)的時(shí)頻空間上,定義頻帶寬度為ΔW,步長(zhǎng)為fp,此處的頻帶寬度ΔW和步長(zhǎng)fp均取為3倍的特征頻率,則在時(shí)頻子空間[0,tN-1, (k-1)fp, (k-1)fp+ΔW]上的重構(gòu)信號(hào)分量為

    (12)

    k=1,2,…,kf

    式中,kf為fs/(2fp)的整數(shù)部分。

    則yk(t)的峭度為

    (13)

    信號(hào)f(t)在頻帶[0,fs/2]的峭度為KR=(Kr(1),Kr(2),…,Kr(kf))。令時(shí)頻子空間[0,tN-1, (k-1)fp, (k-1)fp+ΔW]的中心頻率為

    則各個(gè)時(shí)頻子空間中心頻率的所組成的向量Fc=(fc(1),fc(2),…,fc(kf)),從而得到信號(hào)f(t)的一種時(shí)頻峭度曲線。

    (3)找出信號(hào)時(shí)頻峭度譜最大值相鄰的局部極小值及其對(duì)應(yīng)的頻率,確定子頻帶邊界。在時(shí)頻峭度譜上找出L-1個(gè)局部極小值對(duì)應(yīng)的頻率{ωm|m=1,2,…,L-1}和最大值對(duì)應(yīng)的頻率fmax。從這L-1個(gè)局部極小值中找出與最大值相鄰的2個(gè)局部極小值。 2個(gè)相鄰的局部極小值之間必定存在一個(gè)局部最大峰值,該峰值對(duì)應(yīng)的頻率為調(diào)頻調(diào)制分量的中心頻率,本文就只劃分了3個(gè)頻帶,因此有子頻帶分界點(diǎn)Ω0、Ω1、Ω2、Ω3,其中,Ω0=0,Ω3=π,可把信號(hào)的分析頻帶[0,π]劃分成3個(gè)連續(xù)的、互不交疊子頻帶。

    (4)在子頻帶為[Ωn-1,Ωn](n=1,2,3)上按照式(2)、式(3)構(gòu)造濾波器組,再由式(4)、式(5)求出小波系數(shù)和逼近系數(shù),最終由式(6)、式(7)求出信號(hào)的各個(gè)子頻帶的分解分量。

    2 仿真實(shí)驗(yàn)

    設(shè)軸承損傷的沖擊響應(yīng)分量s(t)=exp(-2πf2t)sin(2πf1t),其中,固有頻率f1=2 000 Hz,f2=70 Hz;損傷沖擊重復(fù)頻率為120 Hz,疊加白噪聲的試驗(yàn)信號(hào)如圖2所示,加噪后信號(hào)的信噪比為-11.76 dB,采樣頻率為12 000 Hz,數(shù)據(jù)長(zhǎng)度為4 096。

    圖2 試驗(yàn)信號(hào)Fig.2 Test signal

    圖3為試驗(yàn)信號(hào)的頻譜和包絡(luò)譜。由圖3a可以隱約觀察到試驗(yàn)信號(hào)以2 000 Hz為中心頻率的共振頻帶。圖3b的包絡(luò)譜上顯示出120 Hz沖擊頻率及其倍頻所對(duì)應(yīng)的峰值。

    (a)信號(hào)頻譜

    (b)信號(hào)包絡(luò)譜圖3 試驗(yàn)信號(hào)的頻譜和包絡(luò)譜Fig.3 Spectrum and envelope spectrum of the test signal

    為了簡(jiǎn)化說(shuō)明,試驗(yàn)僅使用頻譜中的前3個(gè)較大峰值進(jìn)行子頻帶邊界運(yùn)算。圖4所示為采用傳統(tǒng)方法的邊界劃分結(jié)果,可以看到中心頻率2 000 Hz及其邊頻帶被劃分到2個(gè)相鄰子頻帶,因此,EWT分解后,圖5中的子頻帶Ⅱ和子頻帶Ⅲ的信號(hào)分量都會(huì)含有重復(fù)性沖擊響應(yīng)分量,兩個(gè)頻帶分量的包絡(luò)譜都包含120 Hz及其倍頻的譜峰。顯然,噪聲分量對(duì)頻譜的影響改變了頻譜的局部最大值分布,使子頻帶邊界的計(jì)算產(chǎn)生了偏差,出現(xiàn)了不合理的頻帶分割。

    圖6所示為采用改進(jìn)方法的邊界劃分結(jié)果,可以看到,以中心頻率為2 000 Hz及其邊頻帶被劃在子頻帶Ⅱ中。

    圖7為峭度最大峰值所在的子頻帶Ⅱ的分量包絡(luò)譜,3個(gè)顯著的峰值對(duì)應(yīng)的頻率分別為重復(fù)沖擊頻率及其2倍頻、3倍頻。

    圖4 EWT頻譜分割(試驗(yàn)信號(hào))Fig.4 EWT spectrum segmentation(test signal)

    (a)子頻帶Ⅱ的分量包絡(luò)譜

    (b)子頻帶Ⅲ的分量包絡(luò)譜圖5 EWT分量的包絡(luò)譜(試驗(yàn)信號(hào))Fig.5 Envelope spectrum of EWT component(test signal)

    (a)時(shí)頻峭度譜的劃分

    (b)對(duì)應(yīng)的頻譜劃分圖6 改進(jìn)的EWT頻譜分割(試驗(yàn)信號(hào))Fig.6 Improved EWT spectrum segmentation(test signal)

    圖7 改進(jìn)的EWT子頻帶Ⅱ的分量包絡(luò)譜(試驗(yàn)信號(hào))Fig.7 Envelope spectrum of improved EWT subbandⅡ(test signal)

    3 工程應(yīng)用

    東風(fēng)43418型機(jī)車(chē)行走部的軸承型號(hào)為52732QT,滾子直徑為34 mm,軸承內(nèi)徑為160 mm,外徑為290 mm,滾子個(gè)數(shù)為17。實(shí)驗(yàn)數(shù)據(jù)來(lái)自于機(jī)車(chē)軸承試驗(yàn)臺(tái),實(shí)驗(yàn)?zāi)M了不同轉(zhuǎn)速和載荷條件下軸承的工作狀態(tài),采用加速度傳感器采集軸承軸瓦的振動(dòng)信號(hào),采樣頻率為12.8 kHz。分別采用傳統(tǒng)EWT算法和改進(jìn)EWT算法分析軸承的振動(dòng)信號(hào),信號(hào)長(zhǎng)度為8 192。

    圖8為一組無(wú)損傷軸承的振動(dòng)信號(hào),實(shí)驗(yàn)載荷為2 002 kg,轉(zhuǎn)速為378 r/min(6.3 Hz),此時(shí),內(nèi)圈、外圈和滾動(dòng)體的損傷特征頻率分別為61.64 Hz、45.46 Hz 和40.74 Hz。

    圖8 無(wú)損傷時(shí)振動(dòng)信號(hào)Fig.8 Vibration signal without damage

    由圖9可以看出,子頻帶Ⅱ和子頻帶Ⅲ的幅值較大,因此分別對(duì)該子頻帶的分量進(jìn)行包絡(luò)譜分析。

    圖10所示為子頻帶Ⅱ和子頻帶Ⅲ的分量包絡(luò)譜,這兩個(gè)頻帶的分量包絡(luò)譜在6.25 Hz、20.31 Hz處有較大的峰值,它們分別對(duì)應(yīng)轉(zhuǎn)頻及其3倍頻,但由于噪聲的干擾,還存在其他干擾峰值,轉(zhuǎn)頻并不明顯。

    圖11所示為采用改進(jìn)EWT算法時(shí),該信號(hào)以時(shí)頻峭度譜劃的分子頻帶邊界,可以看出,子頻帶Ⅲ的峭度峰值最大,對(duì)該子頻帶的分量進(jìn)行包絡(luò)譜分析。

    圖9 EWT頻譜邊界劃分(無(wú)損傷)Fig.9 Boundary division of EWT spectrum(no damage)

    圖10 EWT分量的包絡(luò)譜(無(wú)損傷)Fig.10 Envelope spectrum of EWT component(no damage)

    圖11 時(shí)頻峭度譜的邊界劃分(無(wú)損傷)Fig.11 Boundary division of time-frequency kurtosis spectrum(no damage)

    圖12中,包絡(luò)譜在6.25 Hz、20.31 Hz處有較大的峰值,它們分別對(duì)應(yīng)轉(zhuǎn)頻及其3倍頻,20.31 Hz是基頻的3倍頻,沒(méi)有發(fā)現(xiàn)其他的損傷特征頻率。

    圖12 改進(jìn)的EWT子頻帶Ⅲ的分量包絡(luò)譜(無(wú)損傷)Fig.12 Envelope spectrum of improved EWT subbandⅢ(no damage)

    圖13為外圈存在損傷的軸承振動(dòng)信號(hào),此時(shí),載荷為1 992 kg,轉(zhuǎn)速為611 r/min(10.18 Hz),計(jì)算可知其內(nèi)圈、外圈和滾動(dòng)體的損傷特征頻率分別為99.64 Hz、73.48 Hz和65.85 Hz。

    圖13 外圈損傷的軸承振動(dòng)信號(hào)Fig.13 Bearing vibration signal of outer ring damage

    圖14所示為采用傳統(tǒng)EWT算法的邊界劃分結(jié)果,可以看出子頻帶Ⅱ的幅值較大,因此對(duì)該子頻帶的分量進(jìn)行包絡(luò)譜分析。

    圖14 EWT頻譜邊界劃分(外圈)Fig.14 Boundary division of EWT spectrum(outer ring)

    圖15中,該頻帶的分量包絡(luò)譜在53.13 Hz、73.44 Hz、79.69 Hz處有較大的峰值,但由于噪聲的干擾,還存在其他峰值,使得故障判斷結(jié)果存在較大誤差。圖16所示為采用改進(jìn)EWT算法時(shí),該信號(hào)基于時(shí)頻峭度譜劃的分子頻帶邊界。

    圖15 EWT子頻帶Ⅱ的分量包絡(luò)譜(外圈)Fig.15 Envelope spectrum of EWT subbandⅡ(outer ring)

    圖16 時(shí)頻峭度譜的邊界劃分(外圈)Fig.16 Boundary division of time-frequency 6kurtosis spectrum(outer ring)

    圖17是最大峭度所在的子頻帶Ⅲ的分量的包絡(luò)譜,可以觀察到包絡(luò)譜上2個(gè)較大峰值對(duì)應(yīng)的頻率分別26.56 Hz、73.44 Hz,它們分別為外圈損傷的特征頻率及其1/3倍頻。

    圖17 子頻帶Ⅲ分量的包絡(luò)譜(外圈)Fig.17 Envelope spectrum of subbandⅢ(outer ring)

    圖18所示為滾動(dòng)體損傷的振動(dòng)信號(hào),此時(shí),載荷為2 001 kg,轉(zhuǎn)速為531 r/min(8.85 Hz),計(jì)算可知,其內(nèi)圈、外圈和滾動(dòng)體的損傷特征頻率分別為86.59 Hz、63.86 Hz和57.23 Hz。

    圖18 滾動(dòng)體損傷的振動(dòng)信號(hào)Fig.18 Vibration signal of rolling element damage

    圖19所示為采用傳統(tǒng)EWT算法的邊界劃分結(jié)果,可以看出子頻帶Ⅱ的幅值較大,因此對(duì)該子頻帶的分量進(jìn)行包絡(luò)譜分析,如圖20所示。該頻帶的分量包絡(luò)譜在26.56 Hz、46.88 Hz、57.81 Hz處有較大的峰值,但由于噪聲的干擾,滾子損傷特征頻率并不明顯,且還存在其他干擾峰值。

    圖20 EWT子頻帶Ⅱ的分量包絡(luò)譜(滾動(dòng)體)Fig.20 Envelope spectrum of EWT subbandⅡ(rolling element)

    改進(jìn)EWT算法分析時(shí),僅以較大的3個(gè)峰值為主進(jìn)行子頻帶劃分,如圖21所示。

    圖21 時(shí)頻峭度譜的邊界劃分(滾動(dòng)體)Fig.21 Boundary division of time-frequency kurtosis spectrum (rolling element)

    同樣,對(duì)最大峰值所在的子頻帶Ⅱ的分量做包絡(luò)譜分析(見(jiàn)圖22),57.81 Hz處有顯著的峰值,該頻率正是滾子損傷的特征頻率,另外,另一個(gè)峰值對(duì)應(yīng)的頻率為26.56 Hz,它是轉(zhuǎn)頻的3倍頻。

    圖22 子頻帶Ⅱ的分量包絡(luò)譜(滾動(dòng)體)Fig.22 Envelope spectrum of subbandⅡ(rolling element)

    4 結(jié)論

    在信號(hào)時(shí)頻分解的基礎(chǔ)上,通過(guò)時(shí)頻區(qū)間分量重構(gòu)的方法來(lái)構(gòu)建時(shí)頻峭度譜,該時(shí)頻峭度譜能夠反映信號(hào)沖擊特性的頻域分布特征;采用時(shí)頻峭度譜代替頻譜確定子頻帶的邊界,克服了噪聲分量對(duì)子頻帶邊界劃分的影響,獲得了合理的子頻帶劃分;以此為基礎(chǔ)構(gòu)造的正交濾波器組能夠有效地分離出信號(hào)的特征分量;時(shí)頻峭度譜的峰值可用于檢測(cè)軸承振動(dòng)信號(hào)子頻帶中的最佳共振頻帶,該最佳共振頻帶可作為選用EWT的分解分量提取損傷特征頻率的依據(jù)。改進(jìn)EWT方法有效地提取了重載工況下機(jī)車(chē)軸承的損傷特征。

    猜你喜歡
    峭度頻帶時(shí)頻
    基于MCKD和峭度的液壓泵故障特征提取
    Wi-Fi網(wǎng)絡(luò)中5G和2.4G是什么?有何區(qū)別?
    聯(lián)合快速峭度圖與變帶寬包絡(luò)譜峭度圖的輪對(duì)軸承復(fù)合故障檢測(cè)研究
    單音及部分頻帶干擾下DSSS系統(tǒng)性能分析
    基于峭度分析的聲發(fā)射故障檢測(cè)
    電子世界(2018年12期)2018-07-04 06:34:38
    雙頻帶隔板極化器
    基于魯棒性小波包峭度圖的滾動(dòng)軸承故障診斷*
    調(diào)諧放大器通頻帶的計(jì)算及應(yīng)用
    基于時(shí)頻分析的逆合成孔徑雷達(dá)成像技術(shù)
    對(duì)采樣數(shù)據(jù)序列進(jìn)行時(shí)頻分解法的改進(jìn)
    久久人人97超碰香蕉20202| videossex国产| 在现免费观看毛片| 国产在线视频一区二区| 成年av动漫网址| 日日啪夜夜爽| 秋霞在线观看毛片| 超碰97精品在线观看| 亚洲少妇的诱惑av| 波野结衣二区三区在线| 极品人妻少妇av视频| 最近中文字幕高清免费大全6| 亚洲欧美日韩另类电影网站| 秋霞伦理黄片| 午夜福利网站1000一区二区三区| 精品久久国产蜜桃| 熟女人妻精品中文字幕| 精品视频人人做人人爽| 欧美xxⅹ黑人| 精品一区二区三区四区五区乱码 | videossex国产| 蜜桃国产av成人99| 日韩,欧美,国产一区二区三区| 国产极品天堂在线| 91久久精品国产一区二区三区| 老司机影院成人| 午夜福利视频在线观看免费| 国产男女超爽视频在线观看| 亚洲成国产人片在线观看| 99久久综合免费| 久久久久久久国产电影| 夫妻性生交免费视频一级片| 国产精品秋霞免费鲁丝片| 在线观看美女被高潮喷水网站| 在线免费观看不下载黄p国产| 国产xxxxx性猛交| 国产熟女欧美一区二区| xxx大片免费视频| 夜夜骑夜夜射夜夜干| 97在线人人人人妻| 亚洲国产精品一区二区三区在线| 国产精品国产三级国产av玫瑰| 9191精品国产免费久久| 久久精品国产亚洲av天美| 九草在线视频观看| 亚洲综合色惰| www日本在线高清视频| 精品亚洲成国产av| 国产福利在线免费观看视频| 免费av中文字幕在线| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 婷婷色麻豆天堂久久| 激情五月婷婷亚洲| 51国产日韩欧美| 日韩欧美一区视频在线观看| 婷婷色麻豆天堂久久| 国产高清国产精品国产三级| 人妻系列 视频| 免费看不卡的av| 22中文网久久字幕| 91久久精品国产一区二区三区| 在线免费观看不下载黄p国产| 18+在线观看网站| 亚洲av中文av极速乱| 丰满乱子伦码专区| 99热这里只有是精品在线观看| 亚洲精华国产精华液的使用体验| 美女主播在线视频| 久久精品国产综合久久久 | 精品午夜福利在线看| 大香蕉97超碰在线| 欧美97在线视频| 亚洲一码二码三码区别大吗| 激情五月婷婷亚洲| 在线亚洲精品国产二区图片欧美| 国产精品久久久久久久电影| 男女边吃奶边做爰视频| 亚洲欧美一区二区三区国产| 人人妻人人爽人人添夜夜欢视频| 视频区图区小说| 高清不卡的av网站| 欧美成人精品欧美一级黄| 在线天堂中文资源库| 久久久久精品久久久久真实原创| 精品卡一卡二卡四卡免费| 国产乱人偷精品视频| 黑丝袜美女国产一区| 亚洲欧美成人综合另类久久久| 国产精品久久久久久av不卡| 新久久久久国产一级毛片| 国产精品久久久久久久电影| 亚洲精品国产色婷婷电影| 最近2019中文字幕mv第一页| 色5月婷婷丁香| 国产精品无大码| 在线观看美女被高潮喷水网站| 欧美亚洲日本最大视频资源| 欧美最新免费一区二区三区| 一级毛片 在线播放| av电影中文网址| 老女人水多毛片| 制服丝袜香蕉在线| 男女国产视频网站| 天堂俺去俺来也www色官网| 国产伦理片在线播放av一区| av免费观看日本| 美女福利国产在线| 国产日韩欧美亚洲二区| av网站免费在线观看视频| 伦理电影免费视频| 久久99热6这里只有精品| 伦精品一区二区三区| 如日韩欧美国产精品一区二区三区| 视频中文字幕在线观看| 男人舔女人的私密视频| 日本黄大片高清| √禁漫天堂资源中文www| 秋霞在线观看毛片| 久久精品人人爽人人爽视色| 高清av免费在线| 国产精品国产三级国产专区5o| 精品熟女少妇av免费看| 久久综合国产亚洲精品| 精品亚洲成a人片在线观看| 一级毛片黄色毛片免费观看视频| 我的女老师完整版在线观看| 亚洲精品日本国产第一区| 欧美变态另类bdsm刘玥| 欧美bdsm另类| 国产精品熟女久久久久浪| 日韩视频在线欧美| 日韩av不卡免费在线播放| 久久久国产一区二区| 日韩 亚洲 欧美在线| 国产精品不卡视频一区二区| 七月丁香在线播放| 高清欧美精品videossex| h视频一区二区三区| 久久久久久久久久久久大奶| 色视频在线一区二区三区| 午夜精品国产一区二区电影| 你懂的网址亚洲精品在线观看| 大片电影免费在线观看免费| 国产精品久久久久久久久免| 男女下面插进去视频免费观看 | 亚洲国产日韩一区二区| 最近最新中文字幕免费大全7| 国产精品熟女久久久久浪| 欧美变态另类bdsm刘玥| 日本wwww免费看| 亚洲精品国产av成人精品| 在线观看免费视频网站a站| 国产黄色免费在线视频| 久久午夜福利片| 日本欧美视频一区| 如日韩欧美国产精品一区二区三区| 18禁在线无遮挡免费观看视频| 久久久a久久爽久久v久久| 交换朋友夫妻互换小说| 视频中文字幕在线观看| 日本欧美视频一区| 欧美人与性动交α欧美精品济南到 | 欧美bdsm另类| 在线 av 中文字幕| 日韩,欧美,国产一区二区三区| 国产精品国产三级国产av玫瑰| 午夜激情av网站| 久久av网站| 人人妻人人澡人人看| 欧美人与性动交α欧美软件 | 国精品久久久久久国模美| 新久久久久国产一级毛片| 国产精品 国内视频| 9色porny在线观看| 在线观看人妻少妇| 国产精品一国产av| 国产精品久久久久久精品电影小说| 色5月婷婷丁香| 亚洲美女视频黄频| 秋霞在线观看毛片| 老司机影院毛片| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩另类电影网站| 插逼视频在线观看| 亚洲情色 制服丝袜| 亚洲图色成人| 久久久久久久久久久久大奶| 久久午夜福利片| 久久这里有精品视频免费| 欧美日韩视频精品一区| 亚洲精品国产av成人精品| 热99久久久久精品小说推荐| 久久精品久久久久久噜噜老黄| 欧美国产精品va在线观看不卡| 国产午夜精品一二区理论片| 亚洲综合色惰| 婷婷色av中文字幕| 中文字幕人妻丝袜制服| 久久精品久久精品一区二区三区| 精品亚洲乱码少妇综合久久| 亚洲精品自拍成人| 99久国产av精品国产电影| 国产伦理片在线播放av一区| 国产欧美另类精品又又久久亚洲欧美| 黄色 视频免费看| 免费观看在线日韩| 人人澡人人妻人| 99久久人妻综合| 久久久亚洲精品成人影院| av在线播放精品| 精品一区二区三卡| 搡女人真爽免费视频火全软件| 国产精品嫩草影院av在线观看| 欧美+日韩+精品| 日韩大片免费观看网站| 国产精品一二三区在线看| 欧美最新免费一区二区三区| 欧美日韩视频高清一区二区三区二| 精品人妻一区二区三区麻豆| 男女边摸边吃奶| 91精品伊人久久大香线蕉| 两个人免费观看高清视频| 观看av在线不卡| 日本欧美国产在线视频| 色94色欧美一区二区| 国产亚洲午夜精品一区二区久久| 高清在线视频一区二区三区| 91国产中文字幕| 99九九在线精品视频| 午夜福利乱码中文字幕| 成人国产av品久久久| 中文字幕最新亚洲高清| videos熟女内射| 在线观看免费日韩欧美大片| 水蜜桃什么品种好| 国产男女超爽视频在线观看| 欧美日韩视频高清一区二区三区二| 91成人精品电影| 日产精品乱码卡一卡2卡三| 国产免费一区二区三区四区乱码| av视频免费观看在线观看| 看非洲黑人一级黄片| 国产亚洲精品久久久com| 18禁国产床啪视频网站| 在线观看三级黄色| 精品亚洲乱码少妇综合久久| 久久精品久久久久久噜噜老黄| 水蜜桃什么品种好| 国产av码专区亚洲av| 欧美精品av麻豆av| 久久久久久久亚洲中文字幕| 日韩伦理黄色片| 国产精品99久久99久久久不卡 | 人妻一区二区av| 成人国语在线视频| 丰满饥渴人妻一区二区三| 最近中文字幕2019免费版| 国产视频首页在线观看| 91成人精品电影| 91国产中文字幕| 欧美国产精品va在线观看不卡| 日本黄色日本黄色录像| 精品第一国产精品| h视频一区二区三区| 日韩伦理黄色片| 国产 精品1| 亚洲经典国产精华液单| 久久青草综合色| 成年女人在线观看亚洲视频| 丝袜脚勾引网站| 人人妻人人澡人人爽人人夜夜| 亚洲情色 制服丝袜| 国产淫语在线视频| 日韩三级伦理在线观看| 狠狠精品人妻久久久久久综合| 色网站视频免费| 啦啦啦在线观看免费高清www| 久久久精品区二区三区| 人体艺术视频欧美日本| 9色porny在线观看| 国产成人午夜福利电影在线观看| 国产激情久久老熟女| 丝瓜视频免费看黄片| 制服丝袜香蕉在线| 男女高潮啪啪啪动态图| 三级国产精品片| 亚洲,欧美精品.| 男的添女的下面高潮视频| 日日撸夜夜添| 最后的刺客免费高清国语| 国产淫语在线视频| 男女边摸边吃奶| 欧美精品国产亚洲| 色婷婷久久久亚洲欧美| h视频一区二区三区| 日本免费在线观看一区| 丰满迷人的少妇在线观看| 狠狠精品人妻久久久久久综合| 久久精品国产亚洲av涩爱| 波多野结衣一区麻豆| 9色porny在线观看| 免费黄频网站在线观看国产| 男女免费视频国产| √禁漫天堂资源中文www| 曰老女人黄片| 日韩精品免费视频一区二区三区 | 2018国产大陆天天弄谢| 国产精品99久久99久久久不卡 | 亚洲av电影在线观看一区二区三区| 欧美成人午夜精品| 亚洲在久久综合| 蜜臀久久99精品久久宅男| 久久午夜综合久久蜜桃| 国产69精品久久久久777片| 国产在线免费精品| 国产成人91sexporn| 热99国产精品久久久久久7| 成人综合一区亚洲| 日韩av不卡免费在线播放| 免费大片18禁| 国产爽快片一区二区三区| 欧美日韩综合久久久久久| 国产精品人妻久久久影院| 亚洲欧美日韩另类电影网站| 91精品三级在线观看| 亚洲欧洲国产日韩| 亚洲,一卡二卡三卡| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| 男女国产视频网站| 如何舔出高潮| 美女脱内裤让男人舔精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美人与性动交α欧美软件 | 午夜久久久在线观看| 在线免费观看不下载黄p国产| 久久精品国产亚洲av涩爱| 波多野结衣一区麻豆| 亚洲高清免费不卡视频| 日本欧美视频一区| 国产女主播在线喷水免费视频网站| 欧美精品av麻豆av| 99热全是精品| 精品国产露脸久久av麻豆| 欧美少妇被猛烈插入视频| 国产亚洲一区二区精品| 精品酒店卫生间| 日韩精品有码人妻一区| 日韩 亚洲 欧美在线| 午夜福利,免费看| 免费观看a级毛片全部| 久久精品久久精品一区二区三区| 在线观看免费视频网站a站| 午夜日本视频在线| 国精品久久久久久国模美| 久久久a久久爽久久v久久| 国产在线视频一区二区| 一本—道久久a久久精品蜜桃钙片| 亚洲国产精品成人久久小说| 蜜桃在线观看..| 国产一区二区在线观看日韩| 最后的刺客免费高清国语| 在线观看免费视频网站a站| 天天躁夜夜躁狠狠久久av| 久久韩国三级中文字幕| 在线观看人妻少妇| 寂寞人妻少妇视频99o| 国产麻豆69| 丝袜在线中文字幕| 新久久久久国产一级毛片| 这个男人来自地球电影免费观看 | 久久狼人影院| 日韩免费高清中文字幕av| 高清欧美精品videossex| av在线老鸭窝| 午夜av观看不卡| 大香蕉97超碰在线| 国产白丝娇喘喷水9色精品| 国产亚洲最大av| 国产色爽女视频免费观看| 最黄视频免费看| 韩国av在线不卡| 伊人亚洲综合成人网| 国产精品嫩草影院av在线观看| 亚洲经典国产精华液单| 毛片一级片免费看久久久久| 国产白丝娇喘喷水9色精品| av卡一久久| 九草在线视频观看| 国产男女内射视频| 精品一品国产午夜福利视频| 九色成人免费人妻av| 国产色爽女视频免费观看| 亚洲一级一片aⅴ在线观看| 乱人伦中国视频| 欧美bdsm另类| 自线自在国产av| 永久免费av网站大全| 欧美国产精品va在线观看不卡| 亚洲av电影在线进入| 亚洲 欧美一区二区三区| 欧美激情国产日韩精品一区| 午夜免费鲁丝| 精品人妻熟女毛片av久久网站| 成人18禁高潮啪啪吃奶动态图| 亚洲精品久久成人aⅴ小说| 国产一区亚洲一区在线观看| 亚洲欧美色中文字幕在线| 十分钟在线观看高清视频www| 精品亚洲乱码少妇综合久久| 亚洲一级一片aⅴ在线观看| tube8黄色片| 美女国产视频在线观看| 丁香六月天网| 色吧在线观看| 久久久国产欧美日韩av| 亚洲一级一片aⅴ在线观看| 国产亚洲精品第一综合不卡 | 久久久a久久爽久久v久久| 视频区图区小说| 丝袜在线中文字幕| 免费观看a级毛片全部| 性色av一级| 亚洲精品日韩在线中文字幕| 99热全是精品| 超色免费av| 99热这里只有是精品在线观看| 91午夜精品亚洲一区二区三区| 狂野欧美激情性bbbbbb| 亚洲精品久久久久久婷婷小说| 最近2019中文字幕mv第一页| 1024视频免费在线观看| 一级黄片播放器| freevideosex欧美| 黄网站色视频无遮挡免费观看| 午夜日本视频在线| 侵犯人妻中文字幕一二三四区| 制服丝袜香蕉在线| 成年人午夜在线观看视频| 寂寞人妻少妇视频99o| www日本在线高清视频| 国产免费一级a男人的天堂| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 久久久久网色| 日韩av不卡免费在线播放| 免费黄色在线免费观看| 国产黄频视频在线观看| 日韩制服丝袜自拍偷拍| 我的女老师完整版在线观看| 夜夜爽夜夜爽视频| 国产男女内射视频| kizo精华| 欧美性感艳星| 精品亚洲成国产av| 美女国产视频在线观看| 久久狼人影院| 久久热在线av| www.av在线官网国产| 亚洲欧美日韩卡通动漫| 欧美激情极品国产一区二区三区 | 国产成人精品无人区| 国产老妇伦熟女老妇高清| 男女下面插进去视频免费观看 | 婷婷色麻豆天堂久久| 久久人人爽人人片av| 国产不卡av网站在线观看| 日韩三级伦理在线观看| 18禁在线无遮挡免费观看视频| videosex国产| 亚洲综合精品二区| 日本黄大片高清| 久久精品国产综合久久久 | 高清欧美精品videossex| 亚洲国产精品一区二区三区在线| 下体分泌物呈黄色| 欧美日韩成人在线一区二区| 国产精品蜜桃在线观看| 亚洲伊人久久精品综合| 久久久久久久亚洲中文字幕| 日韩av免费高清视频| 国产av一区二区精品久久| 国产片内射在线| 欧美+日韩+精品| 人人妻人人添人人爽欧美一区卜| 久久久久国产网址| 一区二区三区四区激情视频| 少妇高潮的动态图| 丝袜美足系列| 狠狠婷婷综合久久久久久88av| 国产精品久久久久成人av| 国产免费福利视频在线观看| 国产 精品1| 最新的欧美精品一区二区| 国产69精品久久久久777片| 精品国产乱码久久久久久小说| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产看品久久| 国产片内射在线| 91精品三级在线观看| 日日摸夜夜添夜夜爱| 久久国产精品大桥未久av| 高清欧美精品videossex| 22中文网久久字幕| 国产亚洲一区二区精品| 日韩伦理黄色片| 精品亚洲成国产av| 国国产精品蜜臀av免费| 精品少妇黑人巨大在线播放| 国产精品久久久久久av不卡| 欧美日韩一区二区视频在线观看视频在线| 精品一区二区三区视频在线| 丝袜人妻中文字幕| 激情视频va一区二区三区| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆| 亚洲精品乱久久久久久| 国产又色又爽无遮挡免| 男男h啪啪无遮挡| 一级毛片黄色毛片免费观看视频| av女优亚洲男人天堂| 久久99蜜桃精品久久| 午夜精品国产一区二区电影| 色婷婷av一区二区三区视频| av黄色大香蕉| 欧美精品国产亚洲| 国产精品人妻久久久久久| 黄片播放在线免费| 伦理电影大哥的女人| 在线 av 中文字幕| 成人手机av| 久久人人97超碰香蕉20202| 亚洲av.av天堂| 日韩欧美一区视频在线观看| 亚洲av电影在线观看一区二区三区| 国产一区二区激情短视频 | 午夜福利视频精品| 久久久欧美国产精品| 99久久中文字幕三级久久日本| 日日摸夜夜添夜夜爱| 亚洲精品中文字幕在线视频| 国产精品一区二区在线观看99| 亚洲成人手机| 欧美激情国产日韩精品一区| 久久热在线av| 亚洲性久久影院| 精品第一国产精品| 色5月婷婷丁香| 9色porny在线观看| 日韩一区二区视频免费看| 女人精品久久久久毛片| 自线自在国产av| 看非洲黑人一级黄片| 999精品在线视频| av播播在线观看一区| 国产无遮挡羞羞视频在线观看| 哪个播放器可以免费观看大片| 亚洲精品456在线播放app| 免费高清在线观看视频在线观看| 三级国产精品片| 1024视频免费在线观看| 在现免费观看毛片| 中文字幕另类日韩欧美亚洲嫩草| 少妇猛男粗大的猛烈进出视频| 国产成人一区二区在线| 多毛熟女@视频| 国产精品久久久av美女十八| 亚洲av成人精品一二三区| 在线观看一区二区三区激情| 制服人妻中文乱码| 夫妻性生交免费视频一级片| 久久久久国产精品人妻一区二区| 看非洲黑人一级黄片| 国产又爽黄色视频| 亚洲av电影在线观看一区二区三区| 色94色欧美一区二区| av电影中文网址| 2022亚洲国产成人精品| 99久久人妻综合| 多毛熟女@视频| videosex国产| 最近2019中文字幕mv第一页| 少妇人妻久久综合中文| 少妇的逼水好多| av国产久精品久网站免费入址| 久久久亚洲精品成人影院| 飞空精品影院首页| 五月玫瑰六月丁香| 赤兔流量卡办理| 亚洲欧美日韩卡通动漫| 国产av精品麻豆| 一二三四在线观看免费中文在 | 亚洲 欧美一区二区三区| 午夜福利在线观看免费完整高清在| 熟女人妻精品中文字幕| 日韩伦理黄色片| 日韩精品免费视频一区二区三区 | 人妻 亚洲 视频| 久久精品久久精品一区二区三区| 国产精品久久久久久久电影| 99热6这里只有精品| 久久这里只有精品19| 日韩欧美一区视频在线观看| av卡一久久| 日产精品乱码卡一卡2卡三| 国产乱人偷精品视频| 国产视频首页在线观看| 老司机影院毛片| 亚洲五月色婷婷综合| 欧美bdsm另类| 一区二区日韩欧美中文字幕 | 另类亚洲欧美激情| 99九九在线精品视频| 国产一区二区三区av在线|