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

    基于Kalman濾波和譜峭度的滾動(dòng)軸承故障診斷

    2014-07-22 01:07:06劉志川唐力偉曹立軍
    軸承 2014年8期
    關(guān)鍵詞:峭度內(nèi)圈外圈

    劉志川,唐力偉,曹立軍

    (軍械工程學(xué)院 火炮工程系,石家莊 050003)

    滾動(dòng)軸承廣泛應(yīng)用于旋轉(zhuǎn)機(jī)械中,是常見的易損件。軸承故障會(huì)導(dǎo)致機(jī)器出現(xiàn)振動(dòng)和噪聲,嚴(yán)重時(shí)會(huì)影響整臺(tái)設(shè)備的正常運(yùn)行,因此軸承的狀態(tài)檢測(cè)和故障診斷具有重要意義。

    共振解調(diào)技術(shù)是目前軸承故障診斷中普遍使用的方法,但是存在帶通濾波器參數(shù)難以確定的問(wèn)題,傳統(tǒng)的小波和小波包分析技術(shù)均不能自適應(yīng)選擇最優(yōu)帶通,導(dǎo)致應(yīng)用受限。譜峭度[1]的基本原理是通過(guò)計(jì)算頻域中每根譜線的峭度值,從而在頻域中定位瞬態(tài)沖擊。文獻(xiàn)[2]提出將譜峭度值作為短時(shí)Fourier窗口函數(shù),通過(guò)譜峭度圖選擇帶通濾波器的參數(shù)。針對(duì)譜峭度圖計(jì)算耗時(shí)的問(wèn)題,文獻(xiàn)[3]提出了快速譜峭度圖的概念。文獻(xiàn)[4]則指出在信噪比較高的情況下,譜峭度能夠較好地識(shí)別瞬態(tài)沖擊信號(hào)。軸承工作環(huán)境相對(duì)復(fù)雜,環(huán)境噪聲較大,而復(fù)雜工況對(duì)譜峭度計(jì)算值會(huì)產(chǎn)生影響,從而給最優(yōu)帶通濾波的帶寬和中心頻率的計(jì)算帶來(lái)誤差,影響診斷效果,因此有必要對(duì)信號(hào)進(jìn)行降噪預(yù)處理。

    Kalman濾波是一種動(dòng)態(tài)數(shù)據(jù)處理方法,屬于遞推式濾波方法[5],可以克服傳統(tǒng)時(shí)頻分析方法不適用于處理非平穩(wěn)信號(hào)的問(wèn)題,對(duì)測(cè)量數(shù)據(jù)逐一處理,由遞推方程隨時(shí)給出某一時(shí)刻新的狀態(tài)估計(jì)值,并與該時(shí)刻后的數(shù)據(jù)一起做出下一狀態(tài)的最優(yōu)估計(jì),達(dá)到信號(hào)降噪目的。在此,針對(duì)軸承內(nèi)、外圈裂紋故障,將Kalman濾波與快速譜峭度圖相結(jié)合,以實(shí)現(xiàn)軸承的故障診斷。

    1 AR模型

    1.1 參數(shù)估計(jì)

    對(duì)于時(shí)間序列{xt},t=1,2,…,N,AR模型為

    xt=φ1xt-1+φ2xt-2+…+φnxt-n+at,

    (1)

    式中:φ2,…,φn及σa是AR模型[6]的參數(shù);n為模型階次。參數(shù)估計(jì)即估計(jì)出這n+1個(gè)數(shù)值,由于

    at=xt-φ1xt-1-φ2xt-2-…-φnxt-n,

    (2)

    (3)

    1.2 模型階數(shù)

    AR模型的階次n未知,需要在遞推過(guò)程中確定。在遞推過(guò)程中,可以給出低階到高階的每一組參數(shù),模型的最小預(yù)測(cè)誤差功率P是遞減的,當(dāng)Pn達(dá)到最小值或者不再變化時(shí),此時(shí)的階次就是AR模型的正確階次。

    定階是建立AR模型的關(guān)鍵過(guò)程,其中AIC準(zhǔn)則(Akaike Information Criterion)應(yīng)用較為廣泛。它由2項(xiàng)構(gòu)成:第1項(xiàng)隨階數(shù)增大而變小,體現(xiàn)模型擬合好壞;第2項(xiàng)隨階數(shù)增大而變大,體現(xiàn)模型參數(shù)多少。AIC準(zhǔn)則表示為

    AIC(n)=Nln(Pn)+2n。

    (4)

    從n=0開始增加模型階數(shù),AIC(n)先減小,后增大。當(dāng)階數(shù)n達(dá)到某一值n0時(shí),AIC(n0)最小,則n0即為AR模型的最合適階次。

    2 Kalman濾波

    Kalman濾波是工程應(yīng)用較為成熟的一種遞推線性最小方差估計(jì),能夠通過(guò)前一時(shí)刻的最優(yōu)值與當(dāng)前時(shí)刻的測(cè)量值一起給出當(dāng)前時(shí)刻的最優(yōu)估計(jì),能夠?qū)y(cè)量數(shù)據(jù)進(jìn)行實(shí)時(shí)更新。通過(guò)對(duì)原始數(shù)據(jù)建立的AR模型可以構(gòu)建Kalman濾波的5個(gè)基本方程。

    設(shè)齒輪箱振動(dòng)真實(shí)信號(hào)為x(t),t=1,2,3,…,N,觀測(cè)噪聲為R,觀測(cè)值為y(t),則觀測(cè)方程可表示為

    y(t)=Hx(t)+R,

    (5)

    其中H=[1 0 … 0]1×n。設(shè)系統(tǒng)狀態(tài)為s(t)=[x(t),x(t-1),…,x(t-n+1)]T,則狀態(tài)轉(zhuǎn)移方程可表示為

    s(t)=As(t-1)+a(t),

    (6)

    其中In-1為n-1維單位陣。由此可建立Kalman濾波的5個(gè)基本方程:

    (1)假設(shè)系統(tǒng)現(xiàn)在狀態(tài)為t,根據(jù)上一時(shí)刻狀態(tài)預(yù)測(cè)現(xiàn)在的狀態(tài)為

    s(t|t-1)=As(t-1|t-1)+a(t);

    (7)

    (2)s(t|t-1)時(shí)刻協(xié)方差更新為

    P(t|t-1)=AP(t-1|t-1)A′+Q;

    (8)

    (3)更新Kalman增益為

    K(t)=P(t|t-1)H′/[HP(t|t-1)H′+R];

    (9)

    (4)新觀測(cè)得到的狀態(tài)最優(yōu)估計(jì)值為

    s(t|t)=X(t|t-1)+K(t)[y(t)-

    Hs(t|t-1)];

    (10)

    (5)s(t|t)時(shí)刻協(xié)方差更新為

    P(t|t)=[I-K(t)H]P(t|t-1),

    (11)

    式中:P為協(xié)方差;Q為系統(tǒng)過(guò)程噪聲;R為測(cè)量噪聲;K(t)為系統(tǒng)增益;I為單位矩陣。

    3 譜峭度

    信號(hào)的Wold-Cramer分解在頻域的表達(dá)式為

    (12)

    式中:H(t,f)為t時(shí)刻沖擊響應(yīng)h(t,f)的頻域表示,是一個(gè)復(fù)包絡(luò)函數(shù)。

    定義Y(t)過(guò)程的4階譜累積量為

    (13)

    2n階的譜瞬時(shí)矩為

    S2nY(f)E{|H(t,f)dX(f)|2n}/df。

    (14)

    則譜峭度定義為

    譜峭度[7-8]是關(guān)于中心頻率和帶寬的函數(shù),當(dāng)帶寬無(wú)限小時(shí),譜峭度為0,當(dāng)帶寬過(guò)大時(shí),譜峭度無(wú)法反映頻帶范圍中的瞬態(tài)沖擊現(xiàn)象。快速譜峭度圖的原理是找到一個(gè)中心頻率fs和帶寬Bw,使峭度達(dá)到最大值,這時(shí)對(duì)應(yīng)的中心頻率和帶寬即最優(yōu)濾波器的2個(gè)參數(shù)??焖僮V度圖是平面上的二維圖,圖中顏色的深淺代表不同中心頻率和帶寬組合下的峭度值,顏色最深處峭度值最大,對(duì)應(yīng)獲得最優(yōu)濾波器的中心頻率和帶寬。

    4 基于Kalman濾波和譜峭度的故障診斷方法

    軸承工作環(huán)境常常有較大噪聲,為降低振動(dòng)信號(hào)噪聲,提高故障特征,先采用Kalman濾波技術(shù)對(duì)振動(dòng)信號(hào)進(jìn)行降噪預(yù)處理,然后利用快速譜峭度圖確定最優(yōu)濾波器的中心頻率和帶寬。應(yīng)用能量算子解調(diào)法解決傳統(tǒng)Hilbert包絡(luò)解調(diào)過(guò)程中存在的信號(hào)兩端調(diào)制和計(jì)算量大的問(wèn)題,然后對(duì)包絡(luò)進(jìn)行FFT變換得到包絡(luò)譜,將包絡(luò)譜中峰值較大的頻率成分與軸承故障特征頻率相比較,實(shí)現(xiàn)診斷,其流程如圖1所示。

    圖1 軸承故障診斷方法流程圖

    5 試驗(yàn)驗(yàn)證

    軸承內(nèi)、外圈故障實(shí)測(cè)信號(hào)來(lái)自二級(jí)減速器試驗(yàn)平臺(tái),試驗(yàn)軸承為6206型深溝球軸承,鋼球個(gè)數(shù)Z=9,鋼球直徑Dw=9.5 mm,球組節(jié)圓直徑Dpw=46.5 mm。采用線切割分別在軸承的內(nèi)圈和外圈上加工寬度0.5 mm、深度1 mm的細(xì)小裂紋作為故障。采集信號(hào)時(shí),振動(dòng)加速度傳感器B&K4508固定在故障軸承所在的軸承座上,采樣頻率為10 kHz;轉(zhuǎn)速由JN338型轉(zhuǎn)速轉(zhuǎn)矩傳感器測(cè)量。

    5.1 內(nèi)圈故障

    輸入軸轉(zhuǎn)速為610.20 r/min,故障頻率f=27.55 Hz,中間軸轉(zhuǎn)動(dòng)頻率fr=5.08 Hz,采樣時(shí)間為2.5 s。

    首先通過(guò)原始信號(hào)建立系統(tǒng)AR模型,利用AIC準(zhǔn)則計(jì)算得到的AIC曲線如圖2所示,從曲線中可以看出,階數(shù)取50時(shí),AIC有極小值,因此確定模型階數(shù)為50。

    圖2 AR模型的AIC曲線

    Kalman濾波過(guò)程中,過(guò)程噪聲和觀測(cè)噪聲是2個(gè)重要參數(shù)[9]。觀測(cè)噪聲使用小波變換降噪的方法提取,對(duì)振動(dòng)信號(hào)進(jìn)行小波多尺度分解,得到小波系數(shù),設(shè)置合適的閾值,大于閾值的小波系數(shù)置零,對(duì)保留的小波系數(shù)進(jìn)行逆變換即為噪聲信號(hào)。過(guò)程噪聲在濾波過(guò)程中由大量試驗(yàn)數(shù)據(jù)確定。軸承內(nèi)圈故障信號(hào)濾波前、后的時(shí)域圖如圖3所示,從圖中可以明顯看出,濾波后信號(hào)幅值減小,沖擊更加明顯。

    圖3 軸承內(nèi)圈故障信號(hào)濾波前、后時(shí)域圖

    對(duì)濾波后的信號(hào)進(jìn)行譜峭度計(jì)算,繪制快速譜峭度圖,如圖4所示。

    圖4 濾波后內(nèi)圈故障信號(hào)的快速譜峭度圖

    由圖4可知,在中心頻率為4 791.67 Hz、帶寬為416.67 Hz的頻帶內(nèi),峭度達(dá)到最大值35(在3.6層內(nèi)),其他頻帶信噪比較低。以此中心頻率和帶寬為帶通濾波器的最優(yōu)參數(shù)對(duì)信號(hào)進(jìn)行帶通濾波,并對(duì)帶通濾波后信號(hào)進(jìn)行能量算子包絡(luò)解調(diào),結(jié)果如圖5所示。圖中可以看出與內(nèi)圈裂紋故障特征頻率接近的27.48 Hz及其倍頻等,而且故障軸承所在軸的轉(zhuǎn)動(dòng)頻率(5.09 Hz)非常明顯,從結(jié)果可以判斷是內(nèi)圈故障。

    圖5 帶通濾波后內(nèi)圈故障信號(hào)時(shí)頻圖

    5.2 外圈故障

    輸入軸轉(zhuǎn)速為1 488.5 r/min,故障頻率f=44.41 Hz,中間軸轉(zhuǎn)動(dòng)頻率fr=12.4 Hz,采樣時(shí)間為1 s。對(duì)采集到的信號(hào)進(jìn)行Kalman濾波及譜峭度計(jì)算,繪制的快速譜峭度圖如圖6所示。由圖可知,在中心頻率為4 375 Hz、帶寬為1 250 Hz的頻帶內(nèi)峭度達(dá)到最大值13.6(在2層內(nèi))。對(duì)信號(hào)進(jìn)行帶通濾波和能量算子包絡(luò)解調(diào),結(jié)果如圖7所示,可以看出與外圈裂紋故障特征頻率接近的44.26 Hz及其倍頻等,而且故障軸承所在軸的轉(zhuǎn)動(dòng)頻率(12.52 Hz)非常明顯,從結(jié)果可以判斷是外圈故障。

    圖6 濾波后外圈故障信號(hào)的快速譜峭度圖

    圖7 帶通濾波后外圈故障信號(hào)時(shí)頻圖

    6 結(jié)束語(yǔ)

    為解決傳統(tǒng)包絡(luò)分析方法中心頻率和帶寬參數(shù)難以確定的問(wèn)題,提出了基于Kalman濾波和譜峭度圖的軸承故障診斷方法。首先建立系統(tǒng)AR模型,對(duì)原始信號(hào)進(jìn)行Kalman濾波降噪預(yù)處理,提高信噪比,然后利用譜峭度自適應(yīng)確定最優(yōu)濾波器的中心頻率和帶寬,最后利用能量算子解調(diào)法對(duì)帶通濾波后的信號(hào)進(jìn)行包絡(luò)解調(diào)分析,結(jié)合軸承故障特征頻率實(shí)現(xiàn)故障診斷。

    對(duì)工程信號(hào)的分析處理驗(yàn)證了本方法的可行性,而且在轉(zhuǎn)速相對(duì)較低的情況下也可以成功診斷出軸承內(nèi)、外圈故障,為低速無(wú)載荷軸承故障診斷提供了新的研究方向。

    猜你喜歡
    峭度內(nèi)圈外圈
    基于MCKD和峭度的液壓泵故障特征提取
    深溝球軸承外圈表面凹坑缺陷分析
    哈爾濱軸承(2022年1期)2022-05-23 13:13:16
    特種復(fù)合軸承內(nèi)圈推力滾道磨削用工裝設(shè)計(jì)
    哈爾濱軸承(2021年4期)2021-03-08 01:00:48
    角接觸球軸承外圈鎖口高度自動(dòng)檢測(cè)規(guī)改進(jìn)
    哈爾濱軸承(2020年2期)2020-11-06 09:22:34
    聯(lián)合快速峭度圖與變帶寬包絡(luò)譜峭度圖的輪對(duì)軸承復(fù)合故障檢測(cè)研究
    主軸軸承內(nèi)圈鎖緊用臺(tái)階套的裝配
    基于峭度分析的聲發(fā)射故障檢測(cè)
    電子世界(2018年12期)2018-07-04 06:34:38
    基于魯棒性小波包峭度圖的滾動(dòng)軸承故障診斷*
    內(nèi)圈帶缺陷中介軸承的動(dòng)力學(xué)建模與振動(dòng)響應(yīng)分析
    軸承內(nèi)圈與軸的配合過(guò)盈量分析
    軸承(2011年10期)2011-07-25 01:36:22
    亚洲性久久影院| 久久久a久久爽久久v久久| 九九在线视频观看精品| 精品无人区乱码1区二区| 亚洲综合色惰| 人人妻人人澡欧美一区二区| 天天躁夜夜躁狠狠久久av| 国内精品美女久久久久久| 国产精品国产高清国产av| 乱码一卡2卡4卡精品| 中出人妻视频一区二区| 国产69精品久久久久777片| 久久99热6这里只有精品| 国产 一区精品| 最好的美女福利视频网| 神马国产精品三级电影在线观看| ponron亚洲| 国产一级毛片七仙女欲春2| 日本黄色视频三级网站网址| 欧美日本视频| 麻豆一二三区av精品| 亚洲图色成人| 啦啦啦啦在线视频资源| 麻豆乱淫一区二区| 非洲黑人性xxxx精品又粗又长| 久久韩国三级中文字幕| 久久久久久久午夜电影| 长腿黑丝高跟| 亚洲av美国av| 丰满的人妻完整版| 久久99热这里只有精品18| www日本黄色视频网| 亚洲精品亚洲一区二区| 午夜福利在线观看免费完整高清在 | 久久精品人妻少妇| 免费观看在线日韩| 国产黄片美女视频| 看免费成人av毛片| 搞女人的毛片| 中国美白少妇内射xxxbb| 中文字幕av在线有码专区| 免费观看精品视频网站| 精品一区二区三区av网在线观看| 美女 人体艺术 gogo| 丝袜美腿在线中文| 亚洲成人av在线免费| 久久午夜亚洲精品久久| 国产精品福利在线免费观看| 国产淫片久久久久久久久| 黄色配什么色好看| 国产淫片久久久久久久久| 女生性感内裤真人,穿戴方法视频| 亚洲中文字幕日韩| 色综合站精品国产| 看十八女毛片水多多多| 天堂av国产一区二区熟女人妻| 成人二区视频| 午夜老司机福利剧场| 欧美性猛交黑人性爽| 国产高清激情床上av| 99久久久亚洲精品蜜臀av| 白带黄色成豆腐渣| 国产精品亚洲美女久久久| 淫妇啪啪啪对白视频| 亚洲精品乱码久久久v下载方式| 国产亚洲欧美98| 十八禁网站免费在线| 国产老妇女一区| 亚洲七黄色美女视频| 日韩中字成人| 亚洲人成网站在线观看播放| 啦啦啦观看免费观看视频高清| 日日啪夜夜撸| 看十八女毛片水多多多| 亚洲av中文字字幕乱码综合| 免费观看人在逋| 久久人人爽人人爽人人片va| 欧美日本亚洲视频在线播放| 一级毛片aaaaaa免费看小| 91久久精品国产一区二区成人| 免费不卡的大黄色大毛片视频在线观看 | 国产真实乱freesex| 日韩欧美国产在线观看| 精品一区二区三区视频在线观看免费| 国产精品一区二区三区四区久久| 精品欧美国产一区二区三| 一本精品99久久精品77| 精品一区二区三区人妻视频| 国产av不卡久久| 俺也久久电影网| 色尼玛亚洲综合影院| 国内精品一区二区在线观看| avwww免费| 两个人视频免费观看高清| 成年女人毛片免费观看观看9| 97超碰精品成人国产| 成人国产麻豆网| 观看美女的网站| 99热全是精品| 欧美人与善性xxx| 免费看a级黄色片| 露出奶头的视频| 一本一本综合久久| 欧美潮喷喷水| av在线观看视频网站免费| 亚洲国产欧美人成| 成人鲁丝片一二三区免费| 少妇的逼水好多| 一区二区三区高清视频在线| 成人漫画全彩无遮挡| 亚洲欧美清纯卡通| 啦啦啦韩国在线观看视频| 国产精品一区二区性色av| 亚州av有码| 国产精华一区二区三区| 天堂av国产一区二区熟女人妻| 精品人妻偷拍中文字幕| 日韩欧美免费精品| 内地一区二区视频在线| 黄色日韩在线| 天堂动漫精品| 非洲黑人性xxxx精品又粗又长| 在线国产一区二区在线| 亚洲欧美日韩无卡精品| 18禁在线无遮挡免费观看视频 | 色综合亚洲欧美另类图片| 97热精品久久久久久| 天堂影院成人在线观看| 日韩欧美 国产精品| 女生性感内裤真人,穿戴方法视频| 狂野欧美白嫩少妇大欣赏| 波多野结衣巨乳人妻| 波多野结衣高清无吗| 听说在线观看完整版免费高清| 在线看三级毛片| 国产成人一区二区在线| 日本 av在线| 美女内射精品一级片tv| 黄色欧美视频在线观看| 婷婷精品国产亚洲av在线| 亚洲性夜色夜夜综合| 麻豆精品久久久久久蜜桃| 日韩成人伦理影院| 久久这里只有精品中国| 别揉我奶头~嗯~啊~动态视频| 亚洲人成网站在线观看播放| 国产不卡一卡二| 一个人观看的视频www高清免费观看| 亚洲人成网站在线播| 黄色日韩在线| 亚洲精品粉嫩美女一区| 久久久国产成人精品二区| 亚洲性久久影院| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久com| 一级毛片aaaaaa免费看小| 久久欧美精品欧美久久欧美| 天天躁日日操中文字幕| 国产精品电影一区二区三区| 一级a爱片免费观看的视频| 久久6这里有精品| 亚洲国产精品久久男人天堂| 12—13女人毛片做爰片一| 丰满乱子伦码专区| 亚洲人成网站在线播| 国产视频一区二区在线看| 国产欧美日韩精品亚洲av| 国产精品一区www在线观看| 成人精品一区二区免费| 久久精品国产自在天天线| 亚洲精品久久国产高清桃花| 日韩欧美在线乱码| 热99在线观看视频| 人妻少妇偷人精品九色| 久久久久国产网址| 欧美精品国产亚洲| 美女高潮的动态| 人妻制服诱惑在线中文字幕| 亚洲18禁久久av| 国产美女午夜福利| 欧美不卡视频在线免费观看| 精品久久久久久久久av| 特大巨黑吊av在线直播| 欧美潮喷喷水| 人妻制服诱惑在线中文字幕| 国产精品亚洲一级av第二区| 免费观看精品视频网站| 蜜桃久久精品国产亚洲av| 乱码一卡2卡4卡精品| 亚洲精华国产精华液的使用体验 | 国产av麻豆久久久久久久| 特级一级黄色大片| 人人妻人人澡人人爽人人夜夜 | 少妇熟女aⅴ在线视频| 国产成人精品久久久久久| 最新中文字幕久久久久| 国产精品精品国产色婷婷| 成人无遮挡网站| 成人无遮挡网站| 男女做爰动态图高潮gif福利片| 99久久精品热视频| 久久精品国产亚洲av香蕉五月| 久久久精品94久久精品| 久久午夜亚洲精品久久| 精品人妻一区二区三区麻豆 | 国产免费一级a男人的天堂| 婷婷色综合大香蕉| 男女之事视频高清在线观看| 欧美中文日本在线观看视频| 亚洲中文日韩欧美视频| 永久网站在线| 午夜福利18| 国产av不卡久久| 午夜亚洲福利在线播放| 熟女人妻精品中文字幕| 国产精品久久电影中文字幕| 欧美最黄视频在线播放免费| 国内精品一区二区在线观看| 欧美又色又爽又黄视频| 丰满乱子伦码专区| 麻豆国产av国片精品| 亚洲内射少妇av| 亚洲七黄色美女视频| 禁无遮挡网站| 少妇人妻精品综合一区二区 | 菩萨蛮人人尽说江南好唐韦庄 | 精品久久久久久久末码| 中文字幕免费在线视频6| 久久亚洲国产成人精品v| 国产午夜福利久久久久久| av国产免费在线观看| 国产一区二区在线av高清观看| а√天堂www在线а√下载| 最近2019中文字幕mv第一页| 国产精品久久电影中文字幕| 春色校园在线视频观看| 老司机福利观看| 欧美最黄视频在线播放免费| .国产精品久久| 一a级毛片在线观看| 亚洲国产欧洲综合997久久,| 日本a在线网址| 女人被狂操c到高潮| 久久久国产成人精品二区| 久久久久精品国产欧美久久久| 午夜视频国产福利| 国产免费男女视频| 亚洲精品亚洲一区二区| 色噜噜av男人的天堂激情| 免费大片18禁| 久久久久九九精品影院| 午夜福利成人在线免费观看| 日韩大尺度精品在线看网址| 大又大粗又爽又黄少妇毛片口| 小说图片视频综合网站| 日韩制服骚丝袜av| 亚洲熟妇熟女久久| 我要搜黄色片| 麻豆av噜噜一区二区三区| 久久精品国产亚洲av天美| 91狼人影院| 亚洲经典国产精华液单| 午夜a级毛片| 卡戴珊不雅视频在线播放| 国产精品一区www在线观看| 九九爱精品视频在线观看| 久久久久国产精品人妻aⅴ院| 一区二区三区高清视频在线| 国产精品三级大全| 日韩成人av中文字幕在线观看 | 搡老妇女老女人老熟妇| 黄色视频,在线免费观看| 一级黄片播放器| 一个人看视频在线观看www免费| 男女做爰动态图高潮gif福利片| 国内精品宾馆在线| 国产亚洲精品av在线| 一a级毛片在线观看| 欧美成人a在线观看| 亚洲国产精品国产精品| 精品一区二区三区视频在线观看免费| 变态另类丝袜制服| 免费在线观看成人毛片| 不卡视频在线观看欧美| 少妇丰满av| 最近手机中文字幕大全| 欧美一级a爱片免费观看看| 不卡一级毛片| 免费电影在线观看免费观看| 噜噜噜噜噜久久久久久91| 中文字幕精品亚洲无线码一区| 精品国内亚洲2022精品成人| 91午夜精品亚洲一区二区三区| 欧美+亚洲+日韩+国产| 丝袜喷水一区| 天堂av国产一区二区熟女人妻| 精品一区二区三区视频在线观看免费| 搡老岳熟女国产| 黑人高潮一二区| 国产精品国产高清国产av| 自拍偷自拍亚洲精品老妇| 日韩制服骚丝袜av| 人人妻,人人澡人人爽秒播| 欧美xxxx性猛交bbbb| 中文字幕久久专区| 欧美日韩一区二区视频在线观看视频在线 | av中文乱码字幕在线| 欧美日韩乱码在线| 内地一区二区视频在线| 国产乱人视频| 亚洲精品国产成人久久av| 91精品国产九色| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av不卡在线观看| 91狼人影院| 男女下面进入的视频免费午夜| 国产成人影院久久av| 国产精品不卡视频一区二区| 最近的中文字幕免费完整| 晚上一个人看的免费电影| 亚洲av中文字字幕乱码综合| 亚洲精品在线观看二区| 日日摸夜夜添夜夜爱| 国产在视频线在精品| 欧美一区二区精品小视频在线| 成年女人永久免费观看视频| 成年女人毛片免费观看观看9| 国产伦在线观看视频一区| 久久久久久九九精品二区国产| 男女边吃奶边做爰视频| 最新在线观看一区二区三区| 少妇熟女aⅴ在线视频| 亚洲综合色惰| 99久久精品国产国产毛片| 啦啦啦韩国在线观看视频| 欧美日本亚洲视频在线播放| 日韩欧美精品免费久久| 一级毛片aaaaaa免费看小| 日本色播在线视频| 国产精华一区二区三区| 国产亚洲精品av在线| 免费电影在线观看免费观看| 国产精品一二三区在线看| 国产黄片美女视频| 露出奶头的视频| 内地一区二区视频在线| 人人妻人人澡人人爽人人夜夜 | 欧美一区二区精品小视频在线| 免费观看人在逋| 日韩欧美国产在线观看| 免费在线观看影片大全网站| 黄色视频,在线免费观看| 亚洲人成网站在线播| 不卡视频在线观看欧美| 日韩人妻高清精品专区| 两个人视频免费观看高清| 又黄又爽又刺激的免费视频.| 日日啪夜夜撸| 国产麻豆成人av免费视频| 永久网站在线| 在线观看午夜福利视频| 亚洲av二区三区四区| 成人漫画全彩无遮挡| 午夜日韩欧美国产| 国产黄a三级三级三级人| 亚洲图色成人| 一a级毛片在线观看| 成人二区视频| 国产 一区精品| 欧美3d第一页| 亚洲成人中文字幕在线播放| 精品久久久噜噜| 国产蜜桃级精品一区二区三区| 特大巨黑吊av在线直播| 午夜a级毛片| 亚洲av电影不卡..在线观看| av国产免费在线观看| 免费看av在线观看网站| 国产免费一级a男人的天堂| 免费看美女性在线毛片视频| 级片在线观看| 国产午夜福利久久久久久| 有码 亚洲区| 99久久精品国产国产毛片| 欧美+日韩+精品| 免费电影在线观看免费观看| 国产美女午夜福利| 欧美国产日韩亚洲一区| 日韩欧美 国产精品| 丝袜美腿在线中文| 精品国内亚洲2022精品成人| 欧美日韩精品成人综合77777| 日本欧美国产在线视频| 中文字幕av在线有码专区| 亚洲国产日韩欧美精品在线观看| www日本黄色视频网| 麻豆国产av国片精品| 黄色一级大片看看| 亚洲熟妇中文字幕五十中出| 亚洲一区高清亚洲精品| 身体一侧抽搐| 色播亚洲综合网| 日韩一区二区视频免费看| 日韩欧美一区二区三区在线观看| 成人鲁丝片一二三区免费| av在线蜜桃| 天天躁日日操中文字幕| 成人永久免费在线观看视频| 97超视频在线观看视频| 亚洲va在线va天堂va国产| 村上凉子中文字幕在线| 午夜日韩欧美国产| 男女下面进入的视频免费午夜| 久久综合国产亚洲精品| 永久网站在线| 亚洲18禁久久av| 在线免费十八禁| 亚洲,欧美,日韩| 国内揄拍国产精品人妻在线| 99久国产av精品国产电影| 欧美人与善性xxx| 国产国拍精品亚洲av在线观看| 亚洲av熟女| 在线播放国产精品三级| 麻豆久久精品国产亚洲av| 天堂√8在线中文| 国产视频内射| 三级国产精品欧美在线观看| 少妇猛男粗大的猛烈进出视频 | 国产在线男女| 亚洲av.av天堂| 直男gayav资源| 亚洲av.av天堂| 久久久a久久爽久久v久久| 97超视频在线观看视频| 成人欧美大片| 中国国产av一级| 国产探花极品一区二区| 日韩欧美免费精品| 日韩精品中文字幕看吧| 久久久久九九精品影院| 国产女主播在线喷水免费视频网站 | 久久久精品欧美日韩精品| 精品欧美国产一区二区三| 亚洲丝袜综合中文字幕| 男女啪啪激烈高潮av片| av卡一久久| 一级毛片久久久久久久久女| 久久久a久久爽久久v久久| 九九久久精品国产亚洲av麻豆| 又爽又黄无遮挡网站| 午夜福利在线在线| 国产成人aa在线观看| 精品一区二区三区人妻视频| 性色avwww在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产国拍精品亚洲av在线观看| 欧美高清性xxxxhd video| 此物有八面人人有两片| 亚洲经典国产精华液单| 欧美zozozo另类| 黄色日韩在线| 国内精品宾馆在线| 欧美又色又爽又黄视频| 国产精品伦人一区二区| 国产成人91sexporn| 国产在视频线在精品| 亚州av有码| 精品久久久久久久人妻蜜臀av| 国产精品精品国产色婷婷| 久久鲁丝午夜福利片| 男女视频在线观看网站免费| 亚洲中文字幕一区二区三区有码在线看| 日本欧美国产在线视频| 日韩精品青青久久久久久| 最新在线观看一区二区三区| 国内揄拍国产精品人妻在线| 成人美女网站在线观看视频| 久久久久久大精品| 免费看光身美女| 91久久精品国产一区二区成人| 成人永久免费在线观看视频| 色综合站精品国产| av黄色大香蕉| 看片在线看免费视频| 精品不卡国产一区二区三区| 在线观看午夜福利视频| 看十八女毛片水多多多| 在线观看美女被高潮喷水网站| 丝袜喷水一区| 91麻豆精品激情在线观看国产| 国内久久婷婷六月综合欲色啪| 亚洲美女视频黄频| 老司机福利观看| 日产精品乱码卡一卡2卡三| 午夜爱爱视频在线播放| 亚洲成人av在线免费| 久久韩国三级中文字幕| 久久精品国产亚洲av涩爱 | 亚洲精品456在线播放app| 一级黄色大片毛片| 综合色av麻豆| 成人一区二区视频在线观看| 日韩精品青青久久久久久| 黄色一级大片看看| 日本免费a在线| 99久久精品一区二区三区| 国产一区二区亚洲精品在线观看| 国产精品野战在线观看| 亚洲av.av天堂| 内地一区二区视频在线| avwww免费| 人妻久久中文字幕网| 又粗又爽又猛毛片免费看| 99riav亚洲国产免费| 日本成人三级电影网站| 欧美成人a在线观看| 婷婷亚洲欧美| 夜夜爽天天搞| 黄色一级大片看看| 丝袜美腿在线中文| 亚洲av五月六月丁香网| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| 亚洲真实伦在线观看| 色哟哟·www| 亚洲天堂国产精品一区在线| 婷婷色综合大香蕉| 99热这里只有是精品在线观看| 成人av一区二区三区在线看| 一区二区三区免费毛片| 不卡视频在线观看欧美| 露出奶头的视频| 国产精品人妻久久久久久| 色哟哟·www| 又黄又爽又刺激的免费视频.| 免费电影在线观看免费观看| 日本三级黄在线观看| 男女下面进入的视频免费午夜| 在线看三级毛片| 欧美又色又爽又黄视频| 国产欧美日韩一区二区精品| 在线观看66精品国产| 观看美女的网站| 尾随美女入室| 亚洲精品粉嫩美女一区| 十八禁国产超污无遮挡网站| eeuss影院久久| a级毛色黄片| 国产免费一级a男人的天堂| 97热精品久久久久久| 又黄又爽又免费观看的视频| 免费大片18禁| 欧美高清成人免费视频www| 精品久久国产蜜桃| 精品一区二区三区av网在线观看| 精品国产三级普通话版| 国产v大片淫在线免费观看| 亚洲av中文av极速乱| 好男人在线观看高清免费视频| h日本视频在线播放| 午夜亚洲福利在线播放| 小蜜桃在线观看免费完整版高清| 国产av不卡久久| 丝袜喷水一区| 丰满的人妻完整版| 三级毛片av免费| 大香蕉久久网| 欧美日韩乱码在线| 人妻丰满熟妇av一区二区三区| 婷婷精品国产亚洲av在线| 亚洲国产精品成人久久小说 | 久久精品国产亚洲av天美| eeuss影院久久| 中文字幕免费在线视频6| 日本色播在线视频| 国产精品伦人一区二区| 午夜福利18| 麻豆av噜噜一区二区三区| 日韩一本色道免费dvd| 欧美色视频一区免费| av.在线天堂| 久久久久国内视频| .国产精品久久| 亚洲精品影视一区二区三区av| 看免费成人av毛片| 99国产精品一区二区蜜桃av| 日韩精品中文字幕看吧| 女人十人毛片免费观看3o分钟| 欧美激情国产日韩精品一区| 一级毛片久久久久久久久女| 欧美绝顶高潮抽搐喷水| 久久精品国产自在天天线| 波多野结衣高清作品| 国产aⅴ精品一区二区三区波| 99热只有精品国产| 最近最新中文字幕大全电影3| 国产精品久久久久久精品电影| 在线播放国产精品三级| 免费电影在线观看免费观看| 亚洲自偷自拍三级| 大又大粗又爽又黄少妇毛片口| 国产av不卡久久| 亚洲三级黄色毛片| АⅤ资源中文在线天堂| 日韩一区二区视频免费看| 国产黄a三级三级三级人| 天堂av国产一区二区熟女人妻| 在线观看美女被高潮喷水网站| 尾随美女入室| 亚洲五月天丁香| 久久久久久九九精品二区国产| 三级毛片av免费| 亚洲乱码一区二区免费版| 91精品国产九色|