大胆免费视频_欧美另类人妖_快播av在线_91欧洲亚洲中日韩在线观看_亚?成人一区在线?看_亚洲免费视频免在线观看

 
齒輪工業(yè)領(lǐng)域服務(wù)平臺,行業(yè)技術(shù)的領(lǐng)航者;
把齒輪傳動之脈搏,譜信息時代之新篇!
當(dāng)前位置: 首頁 » 技術(shù)中心 » 最新技術(shù)資訊

參數(shù)化的短時傅里葉變換及齒輪箱故障診斷

發(fā)布時間:2023-10-16 | 來源:振動、測試與診斷 | 作者:包文杰等
   提出了一種基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換時頻分析方法,并將該方法用于行星齒輪箱的故障診斷。該時頻分析方法通過使用快速路徑優(yōu)化獲得瞬時頻率變化規(guī)律,在短時傅里葉變換過程中自適應(yīng)的改變時窗長度,從而獲得更恰當(dāng)?shù)臅r頻分辨率。針對行星齒輪箱運(yùn)行狀態(tài)不穩(wěn)定的特點(diǎn),通過使用筆者提出的時頻分析方法可以有效地提取出行星齒輪箱的轉(zhuǎn)速信息,利用參考轉(zhuǎn)速對故障信號角度域重采樣和階次分析,從而實(shí)現(xiàn)變轉(zhuǎn)速情況下的行星齒輪箱故障診斷。仿真分析表明,與傳統(tǒng)短時傅里葉變換相比基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換得到的時頻分布能量更加集中;試驗(yàn)分析證明了基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換方法在行星齒輪箱故障診斷中的有效性。

  行星齒輪變速箱具有傳動比大、承載能力強(qiáng)等優(yōu)點(diǎn),在現(xiàn)代工業(yè)中得到廣泛應(yīng)用。它們一般在重載和惡劣的工作環(huán)境下工作,其關(guān)鍵部件包括太陽齒輪、行星齒輪及齒圈等,都會出現(xiàn)嚴(yán)重的點(diǎn)蝕和疲勞裂紋。

  近年來,許多研究者對行星齒輪變速箱的故障診斷做出了大量貢獻(xiàn)。其中,一批學(xué)者通過動力學(xué)建模與分析研究出齒輪故障、制造誤差和負(fù)載對齒輪箱振動響應(yīng)的影響。Lei等提出了統(tǒng)計指標(biāo)為常數(shù)情況下的行星齒輪箱狀態(tài)監(jiān)測方法。Bartel- mus等提出了非平穩(wěn)運(yùn)行狀態(tài)下的行星齒輪箱狀態(tài)監(jiān)測。Feng等總結(jié)出行星齒輪箱振動信號的特征頻譜。然而,行星齒輪變速箱一般用于類似風(fēng)力發(fā)電機(jī)組這類運(yùn)行狀態(tài)不穩(wěn)定的情況,此類情況下,行星齒輪箱振動信號特征為非平穩(wěn)特征。在非平穩(wěn)條件下提取行星齒輪變速箱的時變故障信息較為困難。時頻分析方法因?yàn)槟苡行У亟沂緯r變特征,而常被用于提取非平穩(wěn)信號的故障信息。時頻分析提供了一種可用并且強(qiáng)大的方法用于提取包含重要物理參數(shù)的非穩(wěn)態(tài)信號的瞬時頻率特征。關(guān)于瞬時頻率估計的方法,已經(jīng)有了大量的研究,其中基于時頻分布的瞬時頻率估計方法在處理嚴(yán)重噪聲信號方面具有出色的特性和穩(wěn)定性,因此近幾年獲得了大量的關(guān)注。

  為了對瞬時頻率的特征做精確的估計,能否獲得高的時頻分辨率被認(rèn)為是衡量時域分析方法的關(guān)鍵。有3種傳統(tǒng)的時頻分析方法,即短時傅里葉變換、連續(xù)小波變換和魏格納分布。短時傅里葉變換和連續(xù)小波變換應(yīng)用于線性時頻分析領(lǐng)域,這兩種方法通過內(nèi)積來處理信號,并且可以定位時間和頻率。然而,根據(jù)海森伯格不確定原理,這兩種方法都不能同時獲得任意高的時頻分辨率。因此可以認(rèn)為,這些方法在保持良好時間分辨率時僅能保證較低的頻率分辨率,反之亦然。筆者提出的基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換(adaptiveshort timefouriertransform,簡稱 ASTFT)方法,希望通過快速路徑優(yōu)化方法獲得瞬時頻率變化規(guī)律來改變時窗長度,從而獲得更恰當(dāng)?shù)臅r頻分辨率。在此基礎(chǔ)上,再次通過使用快速路徑優(yōu)化方法來提高瞬時頻率估計的準(zhǔn)確性。

  一、基本原理

  ASTFT方法:傳統(tǒng)短時傅里葉變換的時頻窗時寬是固定不變的,但大多數(shù)工程信號的瞬時頻率變化規(guī)律并不恒定,使用時窗固定的短時傅里葉變換處理此類信號效果必然不好。

  自適應(yīng)短時傅里葉變換的時窗長度隨瞬時頻率的斜率(即瞬時頻率的一階導(dǎo)數(shù))變化而改變,在瞬時頻率較大時使用短時窗的時頻核函數(shù),以取得較高的時間分辨率;在瞬時頻率的斜率較小時使用長時窗的視頻核函數(shù),以取得較高的頻率分辨率。

  對于一個純頻率調(diào)制信號,例如:s(t)=eif(t) ,Co- hen給出的最優(yōu)時窗長度和瞬時頻率斜率的關(guān)系如下


  其中:f′(t)為瞬時頻率的一階導(dǎo)數(shù)。

  當(dāng)對信號s(t)采用變窗高斯核函數(shù)的短時傅里葉表達(dá)式時,時窗長度 Tt 和時變的標(biāo)準(zhǔn)差σ(t)之間的關(guān)系為


  結(jié)合式(1)和式(2)可以推出


  因此,當(dāng)時頻核函數(shù)的標(biāo)準(zhǔn)差參數(shù)σ(t)按照式(3)取值時,變窗高斯核函數(shù)的短時傅里葉變換將會得到最優(yōu)分辨率。舉例說明時頻分辨率與瞬時頻率變化關(guān)系如圖1所示。圖1中:紅線為頻率隨時間變化的時頻脊線;橫坐標(biāo)為時間t;縱坐標(biāo)為頻率 f;圖中矩形為時頻塊。從圖中可以看出時頻塊的形狀隨瞬時頻率斜率變化過程,當(dāng)瞬時頻率的斜率變大時,時頻塊的高度δf 變大,寬度δt 變小,時間分辨率變高而頻率分辨率變低。


  基于快速路徑優(yōu)化的脊線提取: 參考文獻(xiàn)給出了脊線的定義。如 圖 2 所示,使用短時傅里葉變換分析了一個二次線性調(diào)頻信號,左側(cè)為信號時域波形,右側(cè)為信號時頻分布圖,圖中紅色部分是時頻矩陣的極大值序列,其包含了大部分的振動能量,通常被認(rèn)為是時頻分布的脊線。信號由式(4)掃頻得到


  其中:瞬時頻率為f(t)=10+300t2


  文獻(xiàn)中給出了一些脊線提取的方法,然而大多數(shù)方法無法準(zhǔn)確恢復(fù)能和信號真實(shí)瞬時頻率相匹配的脊線。例如,基于瞬時頻率估計的多項(xiàng)式線性調(diào)頻小波變換方法在處理嚴(yán)重噪聲信號時效果并不好。最明顯的原因是它只考慮時頻分布中每一離散點(diǎn)的最大幅值。所以,筆者采用了一種包含快速路徑優(yōu)化(afastpathoptimization,簡稱 FPO)的瞬時頻率提取方案。該脊線提取方案的具體步驟與原理如下:

  1)使用時頻變換方法計算信號的時頻分布矩陣,得到每一時刻的幅值極大值 En(t),及其對應(yīng)的頻率 Pn(t)


  2)通過使用Pn(t)的相鄰節(jié)點(diǎn) Pn(t)+Δf 和 Pn(t)-Δf 的二次插值對Pn(t)進(jìn)行修正。優(yōu)化方案如下


  其中:Δf 為時頻矩陣的頻率分辨率;Pn由式(5)計算得出。

  3)通過對相鄰時刻極大值點(diǎn)對應(yīng)頻率的差值進(jìn)行懲罰,構(gòu)造極大值點(diǎn)Pn(tm)的權(quán)重函數(shù)


  其中:w 為保證Pn(tm)與 Pk(tm-1)之間連續(xù)性的懲罰因子


  其中:f為采樣頻率;σ為高斯窗函數(shù)的參數(shù)。

  因此,使用式(8)的權(quán)重函數(shù)構(gòu)造路徑優(yōu)化函數(shù)如下


  4)通過對上述構(gòu)造的路徑優(yōu)化函數(shù)進(jìn)行求解,得出t1~tm 時刻的最優(yōu)組合{P(t1),P(t2),…, P(tM )},即為所提脊線。

  相對于一些瞬時頻率估計方法,此方案有3個優(yōu)點(diǎn):①使用二次插值定位幅值極大值的頻率更加的精確;②用來提取脊線的路徑優(yōu)化函數(shù)不僅考慮到脊線當(dāng)前點(diǎn)的幅值大小,還考慮到脊線的連續(xù)性;③快速路徑優(yōu)化的脊線提取方法可以有效地避開隨機(jī)噪點(diǎn)對所提脊線的干擾。

  基于快速路徑優(yōu)化脊線的 ASTFT:自適應(yīng)短時傅里葉變換能夠根據(jù)瞬時頻率的斜率自動調(diào)整時窗寬度,從而在瞬時頻率斜率小的時 候保證較高的的頻率分辨率,在瞬時頻率斜率大的時候保證較高的時間分辨率。然而,在對原始信號做 ASTFT 之前,要先進(jìn)行瞬時頻率估計,以瞬時頻率的斜率作為變化時窗的參考。本研究提出了使用時頻分布的脊線來計算瞬時頻率斜率,以此計算自適應(yīng)短時傅里葉變換時所使用的變化時窗。其具體實(shí)現(xiàn)步驟為:

  1)對原始信號進(jìn)行短時傅里葉變換得到初步的時頻分布S(t,f);

  2)對S(t,f)進(jìn)行基于快速路徑優(yōu)化的脊線提取得到f(t);

  3)確定自適應(yīng)短時傅里葉變換的標(biāo)準(zhǔn)差參數(shù)

  4)以σ(t)為標(biāo)準(zhǔn)差參數(shù)對s(t)進(jìn)行自適應(yīng)短時傅里葉變換得到 ASTFT(t,f)。

  二、仿真信號分析

  仿真信號構(gòu)造:為驗(yàn)證所提方法在時頻分析方面的有效性,構(gòu)造仿真信號s(t)如下


  對信號進(jìn)行時頻分析前加入信噪比為 -1dB 的高斯白噪聲。采樣頻率為2048Hz,時域波形如圖3所示,其理論瞬時頻率f(t)如下



  仿真信號分析:通過對s(t)做固定時間窗的短時傅里葉變換,可以看出,在瞬時頻率斜率較大的時候(圖4(a)中紅圈部分),長時窗的短時傅里葉變換(高斯窗函數(shù)的標(biāo)準(zhǔn)差參數(shù)σ取0.004,σ越小時窗長度越長)效果較差,能量較分散;在瞬時頻率斜率較小的時候(圖4(b)中紅圈部分),短時窗的短時傅里葉變換(高斯窗函數(shù)的標(biāo)準(zhǔn)差參數(shù)σ取0.04,σ越大時窗長 度越短)效果較差,能量不夠集中。


  同時對s(t)做自適應(yīng)短時傅里葉變換(高斯窗函數(shù)的標(biāo)準(zhǔn)差參數(shù),如式(3)所示,f(t)為文中所提兩種脊線方法得到的脊線), 結(jié)果如圖5所示,相對圖4中時頻分布自適應(yīng)短時傅里葉變換獲得的時頻分布能量更加集中。圖 5(a)為基于最大能量峰值(maximumenergypeaks, 簡稱 MEP)的瞬時頻率估計方法做自適應(yīng)短時傅里葉變換所得時頻分布,由于隨機(jī)噪聲干擾,瞬時頻率出現(xiàn)很大的波動;圖5(b)為筆者提出的基于快速路徑優(yōu)化脊線的自適應(yīng)短時傅里葉變換所得時頻分布,由于該方法可以有效降低隨機(jī)噪聲得干擾,得到的時頻分布較為均勻、集中。為從圖5中提取的瞬時頻率與理論瞬時頻率的比較,從圖6可以看出,基于 FPO的瞬時頻率估計結(jié)果更接近于理論瞬時頻率。


  文中通過計算范數(shù)來衡量脊線的精確程度


 

  其中:fTI為真實(shí)脊線頻率;fEI為所提取脊線頻率值。

  經(jīng)計算可得本算法所得脊線誤差為0.0067,而使用 MEP脊線誤差為0.0361,故本算法所提脊線遠(yuǎn)比 MEP算法脊線準(zhǔn)確。

  三、行星齒輪箱故障診斷

  行星齒輪箱故障數(shù)據(jù):本試驗(yàn)的振動數(shù)據(jù)采集自行星齒輪箱故障模擬試驗(yàn)臺。如圖7所示,試驗(yàn)臺包括電機(jī)驅(qū)動、交流電機(jī)、行星齒輪箱、固定軸齒輪箱和制動器。轉(zhuǎn)速傳感器和加速度傳感器分別采集電機(jī)軸轉(zhuǎn)速和行星齒輪箱振動信號,電機(jī)軸轉(zhuǎn)速可由電機(jī)驅(qū)動自由調(diào)節(jié),本試驗(yàn)中電機(jī)為變轉(zhuǎn)速,包括一段升速和一段減速過程。信號采樣頻率為12800Hz,采集時長16s。


  表1給出了行星齒輪箱的齒輪參數(shù)。由參考文獻(xiàn)提供的星齒輪箱故障階次計算公式如下


  其中:Or,Os 和 Op 分別為齒圈、太陽輪、行星輪故障階次;Zr,Zs 和Zp 分別為齒圈、太陽輪、行星輪齒數(shù);np 為行星輪個數(shù)。


  表2為將表1中參數(shù)帶入式(14)~(16)計算得到的行星齒輪箱內(nèi)各零件故障特征階次表。


  試驗(yàn)信號分析:通過對含有行星輪裂紋故障的振動信號進(jìn)行階次分析,從而驗(yàn)證該方法的有效性。行星齒輪箱故障數(shù)據(jù)的振動時域圖和頻譜圖如圖8所示。圖8(a)為試驗(yàn)中直接采集的振動信息時域圖,圖8(b)為時域信號經(jīng)傅里葉變換得到的頻譜圖,在頻譜中出現(xiàn)了很多峰值,由于齒輪轉(zhuǎn)速波動,很難將頻譜圖中的峰值與故障特征頻率或者故障特征階次聯(lián)系起來。


  首先,對原始故障信號使用短時傅里葉變換獲得信號的時頻分布,并提取出信號中電機(jī)轉(zhuǎn)速對應(yīng)的瞬時頻率,如圖9所示,圖中紅線為使用基于快速路徑優(yōu)化的脊線提取方案獲得的瞬時頻率。


  然后,再次對原始故障信號使用筆者提出的基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換獲得新的時頻分布,并提取出信號中電機(jī)轉(zhuǎn)速對應(yīng)的瞬時頻率,如圖10所示,圖中紅線為從此時頻分布中提取的瞬時頻率。通過對比,可以看出本方法獲得的時頻分布(見圖10)優(yōu)于傳統(tǒng)短時傅里葉變換得到的時頻分布(見圖9)。將圖9和圖10中提出的瞬時轉(zhuǎn)速信息與實(shí)際轉(zhuǎn)速進(jìn)行對比,如圖11所示,圖中曲線1為使用STFT 得到的轉(zhuǎn)速;曲線2為本方法得到的轉(zhuǎn)速;曲線3為實(shí)際轉(zhuǎn)速。可以看出通過本方法獲得的轉(zhuǎn)速更加接近于真實(shí)轉(zhuǎn)速,可以作為實(shí)際轉(zhuǎn)速未知情況下行星齒輪箱故障信號角度域重采樣的參考轉(zhuǎn)速。


  最后,利用獲得的瞬時頻率對故障數(shù)據(jù)重新進(jìn)行角度采樣,即將時域信號對轉(zhuǎn)速進(jìn)行積分從而得到角度域信號,如圖12所示。然后對故障數(shù)據(jù)的角度域信號做階次分析得到振動信號的階次圖,如圖13所示。圖中幅值較大的1階為轉(zhuǎn)頻、2階為2倍轉(zhuǎn)頻。對應(yīng)表2可以在圖13階次中找到故障階次,故障階次為0.83階,幅 值 為0.012g,故障階次非常明顯,從而驗(yàn)證基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換方法在實(shí)際應(yīng)用中的可行性。同時添 加如下對比,以通過 MEP 方法獲得的參考轉(zhuǎn)速對信號做階次分析,得到如圖14所示階次圖,從圖中很難判斷出行星齒輪故障階次。


  四、結(jié)束語

  基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換方法依據(jù)快速路徑優(yōu)化方法獲得的瞬時頻率規(guī)律變換時頻窗長度,使自適應(yīng)短時傅里葉變換能夠獲得更加優(yōu)化的時頻分布,從而在下一次的脊線提取過程中獲得更加準(zhǔn)確的瞬時頻率估計結(jié)果。通過對比仿真信號的處理結(jié)果,可以清晰地看出在處理瞬時頻率斜率變化較快的信號時,筆者提出的基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換要優(yōu)于傳統(tǒng)的短時傅里葉變換方法。通過使用此方法處理取自行星齒輪變速箱的非平穩(wěn)信號,根據(jù)快速路徑優(yōu)化脊線提取方案獲得的齒輪箱瞬時轉(zhuǎn)速對原始故障信號重新進(jìn)行角度采樣,對角度域故障信號做階次分析,得到明顯的特征階次。通過對試驗(yàn)信號的分析證明了基于快速路徑優(yōu)化的自適應(yīng)短時傅里葉變換方法在行星齒輪變速箱故障診斷上的有效性。

  參考文獻(xiàn)略.

[ 技術(shù)中心搜索 ]  [ ]  [ 告訴好友 ]  [ 打印本文 ]  [ 關(guān)閉窗口 ]  [ 返回頂部 ]
主站蜘蛛池模板: 国产一级做a_一个色综合国产色综合_九九九99品牌的特色产品_被男人吃奶跟添下面特舒服_啄木鸟在线_99啪啪_zjzjzjzjzjzjzj视频免费播放_jvid福利在线一区二区 | 色噜噜久久综合伊人一本_伊人久久丁香色婷婷啪啪_久久狠狠爱亚洲综合影院_日韩一级免费观看_日本内射FREERAPE视频_国产欧美一区二区精品秋霞影院_亚洲一区国产二区_中文字幕一二三综合a | 拔插拔插8X8X海外华人免费视频_无码专区亚洲avl_日本黄色片免费播放_亚洲成人精品一区二区三区_国产精品一级在线观看_国产精品无码无在线观看_国语自产偷拍精品视频偷蜜芽_中文字幕亚洲乱码熟女在线萌芽 | 无码高潮喷吹在线播放亚洲_日本午夜在线视频_国产免费牲交视频_久久精品区_日日夜夜精_亚洲国产精品激情在线观看_亚洲sm另类一区二区三区_精品久久影视 | 日韩经典AV在线观看_中文无码日韩欧免费视频_捆绑白丝粉色jk震动捧喷白浆_日本慰安所一级毛片在线播放_亚洲精品久久久乳夜夜欧美_97浪潮性色91久久久美川_亚洲蜜桃精久久久久久久久久久久_最新中文字幕av | xxxx漂亮少妇hdxxxx_久久香蕉国产线看观看手机_青青草手机在线观看_人妻激情另类乱人伦人妻_狠狠色狠狠色88综合日日91_有码在线_日韩av综合_好爽好痛好湿好硬视频免费 | 亚洲三区在线视频_亚洲色图片区_久久久999精品视频_日韩一级黄色影片_黄色片免费在线观看_成人a视频片在线观看免费_青青青伊人_狠狠一区 | 久草高清_中文字幕无码中文字幕有码在线_日韩成年视频_91成人在线免费观看_日本爽爽爽爽爽爽在线观看免_特大黑人娇小亚洲女_91亚洲国产成人_精品少妇亚洲乱码 | 国产亚洲视频系列_中文字幕亚洲欧美在线不卡_国产精品成人A片在线果冻_婷婷超碰_日本a级片免费看_XXX国产精品XXX_欧美一级片在线播放_亚洲综合清纯唯美 | 一级肉体全黄裸片免费观看_免费观看成年人网站_WWW插插插无码免费视频网站_国产精品黄色小视频_999九九九久久久_51自拍视频_国产午精品午夜福利757视频播放_日本久久久www成人免费毛片 | 免费www_他一边曰一边吃我奶小说免看_日产精品卡1卡2免费_国产精品呻吟久久av凹凸_免费网站在线观看黄色_亚洲欧美日韩在线一区二区_白袜奴gay网站_欧美黄网站在线观看 午夜精品中文字幕_贪婪欲望之岛在线播放_爱情岛亚洲AV永久入口首页_免费看黄视频在线_亚洲国产A∨无码中文777_中文字字幕码一二区_国产精品亚洲а∨天堂网_国产明星裸体无码XXXX视频 | 在线观看啊啊啊啊_精品亚洲国产成人A片在线观看_国产精品永久久久久久久www_黄色国产视频_欧洲美女粗暴牲交免费观看_18国产一二三精品国产_一区二区三区日本在线观看_在线观看视频99 | 欧美精品毛片久久久久久久_日韩人妻中文无码一区二区七区_久草最新在线_5566先锋影音夜色资源站在线观看_亚洲乱码国产乱码精品精可以看_国产精品呻吟高潮_成人看片黄A免费看那个网址_久久欧美AⅤ无码精品色午夜麻 国产成人一区二区三区影院动漫_国产成人在线视频免费观看_久久99女女久久99久久_www.久久爱.com狼人_国产又爽又刺激的视频_日韩中文在线观看_天天草天天爱_狠狠爱天天操 | 日韩精选在线观看_a级毛片在线免费_久久精品国产亚洲AV日韩_国产一区二区三区影视_亚洲欧美天堂_伊人久久大香线蕉影院_久久久久成人精品亚洲国产_永久免费精品精品永久 | 又大又粗进去爽A片免费_青青青手机在线_强被迫伦轩高潮无BD_韩国V欧美V亚洲V日本V_精品国一区二区三区_人妻少妇久久中文字幕_蜜臀AV在线无码国产_天天摸天天干 | 国产黄在线观看免费观看不卡_av福利免费_国产麻豆天美果冻无码视频_午夜宅男久久久_xxxxxhd日本xxxx47_日韩高清在线一区二区三区_男人天堂视频网_少妇高潮浪潮av99网站 | 一本精品999爽爽久久久_久久av.com_最近最新中文字幕h_国产精品国产三级国产专播一_最近免费中文字幕大全高清_久久毛片一区二区三区_三年片免费观看大全国语_欧美精品无码久久久 | 狂躁美女bbbbbb视频_免费两性的视频网站_在线色_国产jjzz_久久久久久久久久久久网站_日日做夜夜爽毛片麻豆_久久少妇视频_91在线国产最新地址 | 欧美jizzhd精品欧美另类_粉嫩一区二区三区四区公司1_美女视频一级片_中文字幕无线码入口_毛片免费小视频_国产日韩欧美一区二区久久精品_免费国产黄网站在线观看动图_yellow在线观看免费观看大全 | 香蕉久久夜色精品国产_亚洲女人被黑人巨大进入_6969视频_综合色久_国产乱对白刺激视频_中文无码一区二区视频在线播放量_性爱视频在线免费_欧美特黄色片 | 午夜dj免费视频在线观看_日韩在线观看视频免费_神马久久春色视频_国产精华7777777_亚洲一区影院_九九久久网_国产日韩一区二区在线_国产精品亚欧美一区二区 | 兔子先生节目在线观看_gogogo高清在线播放韩国_林ゆな中文字幕一区二区_91精品国产91久久久久久不卡_一二三四影视在线观看免费视频_91亚洲日本aⅴ精品一区二区_成人在线午夜视频_国产在线一区二区三区在线观看 | 精品午夜久久_色惰日本视频网站_夜色综合网_亚洲中文有码字幕青青_98成人网_亚洲精选av_国产一区影视_免费观看电视在线高清视频大全 | japanese在线观看_亚洲精品黄网站_精品在线一区二区三区_双乳奶水饱满少妇视频_亚洲日本韩国欧美_欧美不卡一区二区_日韩免费大片_国产精品不卡一区二区三区 | 伊人思思_亚洲精品无码久久_中文字幕成人在线视频_国产精品免费视频软件_在线观看91免费视频_婷婷激情综合色五月久久竹菊影视_手机看日韩片_日韩视频成人 | www视频免费观看_麻豆freehdxxxⅹ传媒_abab456成人免费网址_国产激情自拍视频_青楼妓女禁脔道具调教SM_国产一区二区三区成人欧美日韩在线观看_久久久久av_精品欧美乱码久久久久久1区2区 | 天天干视频网_国产精品久久久久免费a∨大胸_在线亚洲日本_成人αv_精品人妻AV一区二区三区_最新理论三级中文在线观看_久久久免费观看_最好看的中文2018 | 91视频综合_亚洲自拍偷拍区_偷国产乱人伦偷精品视频_91九色视频网_黄色片在线观看免费_日韩精品一区二区在线观看_亚洲一级精品毛片_成年动漫av网免费 | 伊人亚洲视频_完美世界在线观看免费完整_黑人黄色一级片_久久天天躁夜夜躁狠狠ds005_操女人视频网站_51嘿嘿嘿国产精品伦理_高清av免费_多人乱p欧美在线观看 | 无遮挡拍拍拍免费观看_99日韩_亚洲国产欧美日韩另类_毛片免费全部播放无码_狂野欧美激情性xxxx喷水_视频在线国产_超碰pron人人_免费h动漫无码网站 | 嫩草精品_91精品一区二区三区四区_欧洲美女粗暴牲交_日韩精品无码一区二区三区av_伊人青草_内射爽无广熟女亚洲_欧美黄片免费看_久久本道综合色狠狠五月 | 福利片免费在线观看_久久久久久九九九九九九_欧美性xxxx69_亚洲精品aaa_亚州国产精品久久久_亚洲乱码无码永久不卡在线_成人毛毛片_久久精品一区二区免费播放 | 亚洲一区AV在线观看_欧美熟妇zozozo_热久久精品国产_日本丰满少妇XXXX_97超碰香蕉_国产日韩欧美视频在线观看_草草视频在线播放_欧美成人午夜免费全部完 | 一区在线播放视频_天堂网色_日韩国产网站_国产国语一级A毛片高清视频_色777狠狠综合秋免鲁丝_免费特级黄色片_久久亚洲第一_caoporn国产免费 | 益日韩欧群交P片内射中文_1313午夜精品理论片_综合爱爱网_纯肉无遮挡h肉动漫在线观看3d_男人的天堂久久精品_99久久精品午夜一区二区_亚洲七七久久桃花影院_国产91色在线亚洲 | 无毛一级片_9I精品福利一区二区三区蜜桃_精品免费观看_99热这里只有精品地址_亚洲国产精品一区二区成人_久久久精品一品道一区_国产的精品一区二区在线观看_国产亚洲精品久久yy5099 | 男啪女视频免费观看网站_国产成人啪精品视频网站午夜_欧亚久久日韩av久久综合_9a蜜桃久久久久久免费_国产综合18久久久久久_国产www_国产精品久久久久久久久亚洲美女_日日爱网站 | 美女91_一本大道久久a久久精品综合_青草久操_亚洲伦理99热久久_狠狠色噜噜狠狠狠888米奇_久久美女色视频_777影音_狠狠色综合7777久夜色撩人 | 黄色片特级片_亚洲性激情_在线一区二区三区四区_成人看片黄a在线观看_我的逼好痒_av软件大全_国产免费福利在线观看_粉嫩极品国产在线观看 | 日本亚洲9999aⅴ_亚洲欧美99_欧美嫩交一区二区三区_亚洲成人免费网站_日本三级视频在线观看_免费黄色一级片_国产无人区卡一卡二卡三乱码_亚洲精品成人久久一区二区 | sxeoquenetv另类重口_日本视频网址_日韩成人在线一区_综合狠狠开心_91久久高清国语自产拍_久草在线免费资源站_国产精品亚洲专区无码不卡_欧洲美熟女乱又伦av影片 |