• <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
    女人被躁到高潮嗷嗷叫费观| 欧美日韩av久久| 国产精品.久久久| 亚洲第一区二区三区不卡| 欧美精品一区二区大全| 色网站视频免费| 精品国产露脸久久av麻豆| 国产成人午夜福利电影在线观看| 国产精品一区二区在线观看99| 亚洲中文av在线| 国产在线免费精品| 波野结衣二区三区在线| 亚洲一区二区三区欧美精品| 日韩伦理黄色片| 免费在线观看黄色视频的| 十分钟在线观看高清视频www| 久久久久久伊人网av| 欧美日韩亚洲国产一区二区在线观看 | 9热在线视频观看99| 九九爱精品视频在线观看| 国产伦理片在线播放av一区| 这个男人来自地球电影免费观看 | 久久国内精品自在自线图片| 婷婷成人精品国产| 考比视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 色吧在线观看| 美国免费a级毛片| 18禁观看日本| 国产精品.久久久| 精品国产乱码久久久久久小说| 精品一区二区三区四区五区乱码 | 亚洲色图综合在线观看| 日韩制服骚丝袜av| 人人妻人人澡人人爽人人夜夜| 一边摸一边做爽爽视频免费| 各种免费的搞黄视频| 人成视频在线观看免费观看| 一区二区日韩欧美中文字幕| 午夜激情av网站| 两个人看的免费小视频| 一本大道久久a久久精品| 一级黄片播放器| 久久久久久久国产电影| 精品99又大又爽又粗少妇毛片| 久久久久久久亚洲中文字幕| 91精品三级在线观看| 国产有黄有色有爽视频| 婷婷色麻豆天堂久久| 国产欧美日韩综合在线一区二区| 国产精品三级大全| 久久精品国产综合久久久| 亚洲国产精品成人久久小说| 美女午夜性视频免费| 国产免费现黄频在线看| 精品久久蜜臀av无| 国产人伦9x9x在线观看 | 精品久久蜜臀av无| 日韩大片免费观看网站| 一本色道久久久久久精品综合| 高清av免费在线| xxxhd国产人妻xxx| 亚洲精品日本国产第一区| 黄片播放在线免费| 国产av一区二区精品久久| av一本久久久久| 一级,二级,三级黄色视频| 日本黄色日本黄色录像| 91精品三级在线观看| 香蕉丝袜av| 日韩人妻精品一区2区三区| 黄片无遮挡物在线观看| 亚洲一级一片aⅴ在线观看| 精品一区在线观看国产| 9191精品国产免费久久| 亚洲三区欧美一区| 精品亚洲成a人片在线观看| 男人添女人高潮全过程视频| 1024香蕉在线观看| 国产探花极品一区二区| 亚洲一码二码三码区别大吗| 精品国产一区二区久久| 久久 成人 亚洲| 97精品久久久久久久久久精品| 国产成人aa在线观看| 久久久a久久爽久久v久久| 欧美人与善性xxx| 日本-黄色视频高清免费观看| 欧美精品一区二区免费开放| 男女国产视频网站| 成年人午夜在线观看视频| 香蕉精品网在线| 欧美精品高潮呻吟av久久| av视频免费观看在线观看| www.av在线官网国产| 成人毛片a级毛片在线播放| 久久精品国产自在天天线| 一级片'在线观看视频| 国产av国产精品国产| 国产精品久久久久成人av| 丝袜在线中文字幕| 在线观看免费视频网站a站| 国产欧美亚洲国产| 欧美日韩成人在线一区二区| 99热网站在线观看| 看免费av毛片| 国产成人aa在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲一区中文字幕在线| 色视频在线一区二区三区| 韩国高清视频一区二区三区| 看免费av毛片| 午夜激情av网站| 丝瓜视频免费看黄片| 国产又色又爽无遮挡免| 香蕉丝袜av| 免费在线观看黄色视频的| 亚洲欧洲精品一区二区精品久久久 | 国产又爽黄色视频| 老汉色av国产亚洲站长工具| 男的添女的下面高潮视频| 欧美bdsm另类| 国产精品三级大全| 成人影院久久| 人人澡人人妻人| 最近手机中文字幕大全| 婷婷成人精品国产| 亚洲视频免费观看视频| 制服丝袜香蕉在线| 又大又黄又爽视频免费| 国产高清国产精品国产三级| 亚洲精品国产av成人精品| 亚洲国产欧美日韩在线播放| 欧美精品高潮呻吟av久久| 五月开心婷婷网| 男女啪啪激烈高潮av片| 日本vs欧美在线观看视频| 亚洲国产精品一区三区| 国产片特级美女逼逼视频| 中文字幕制服av| 宅男免费午夜| 中文字幕人妻熟女乱码| 国产在视频线精品| h视频一区二区三区| 超色免费av| av在线播放精品| 国产深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 一边亲一边摸免费视频| 超碰97精品在线观看| 老汉色av国产亚洲站长工具| 亚洲精品国产一区二区精华液| 国产精品国产三级国产专区5o| 美女高潮到喷水免费观看| 日本vs欧美在线观看视频| 成年人午夜在线观看视频| 亚洲国产精品999| 日本wwww免费看| 欧美国产精品va在线观看不卡| h视频一区二区三区| 日本av免费视频播放| av国产久精品久网站免费入址| 国产又爽黄色视频| 五月伊人婷婷丁香| 久久久久网色| 久久精品人人爽人人爽视色| 性色av一级| 性高湖久久久久久久久免费观看| 国产精品三级大全| 久久久久久免费高清国产稀缺| 97人妻天天添夜夜摸| 亚洲av男天堂| 国产精品亚洲av一区麻豆 | 亚洲美女视频黄频| 国产日韩一区二区三区精品不卡| 亚洲国产精品一区三区| av天堂久久9| av在线app专区| 精品亚洲乱码少妇综合久久| 国产xxxxx性猛交| 丝袜喷水一区| 久久精品aⅴ一区二区三区四区 | 午夜福利在线观看免费完整高清在| 少妇人妻 视频| 在线观看国产h片| 亚洲国产最新在线播放| 在线观看www视频免费| 国产成人精品一,二区| 日韩熟女老妇一区二区性免费视频| 国产精品一国产av| www.精华液| 99久国产av精品国产电影| 黄色 视频免费看| 赤兔流量卡办理| 国产免费福利视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 老司机影院成人| 国产成人一区二区在线| 午夜久久久在线观看| 天堂俺去俺来也www色官网| 美女脱内裤让男人舔精品视频| 免费av中文字幕在线| 中文字幕亚洲精品专区| 亚洲精品久久久久久婷婷小说| 毛片一级片免费看久久久久| 精品一品国产午夜福利视频| 99热网站在线观看| 激情视频va一区二区三区| 国产精品一区二区在线观看99| www.熟女人妻精品国产| 老女人水多毛片| 最新的欧美精品一区二区| 99九九在线精品视频| 国产不卡av网站在线观看| 少妇人妻 视频| 五月天丁香电影| 在线精品无人区一区二区三| 欧美中文综合在线视频| 极品少妇高潮喷水抽搐| 青草久久国产| 国产成人aa在线观看| 国产成人精品在线电影| www.熟女人妻精品国产| 亚洲第一av免费看| 国产麻豆69| 午夜激情av网站| 亚洲成色77777| 久久久精品94久久精品| 老女人水多毛片| 两性夫妻黄色片| 超碰成人久久| 久久久久久久久久人人人人人人| 欧美精品av麻豆av| 亚洲精品久久成人aⅴ小说| 免费在线观看黄色视频的| 国产成人aa在线观看| 久久久久久久久久人人人人人人| 亚洲精品av麻豆狂野| 日韩电影二区| 亚洲国产av影院在线观看| 亚洲三区欧美一区| 久久久亚洲精品成人影院| 啦啦啦视频在线资源免费观看| 日韩精品免费视频一区二区三区| 一本色道久久久久久精品综合| 欧美成人午夜精品| 大话2 男鬼变身卡| 国产精品久久久久久av不卡| 中文字幕人妻丝袜一区二区 | 欧美日韩国产mv在线观看视频| h视频一区二区三区| 最黄视频免费看| 欧美日韩av久久| 亚洲精品美女久久久久99蜜臀 | 精品国产乱码久久久久久小说| 国产精品免费视频内射| 男女免费视频国产| 99九九在线精品视频| 亚洲在久久综合| 欧美国产精品一级二级三级| 国产日韩欧美在线精品| 欧美日韩精品成人综合77777| 免费观看无遮挡的男女| 久久国内精品自在自线图片| 久久久久国产一级毛片高清牌| 亚洲综合精品二区| 午夜福利乱码中文字幕| av在线观看视频网站免费| 久久久久久久精品精品| 免费黄频网站在线观看国产| 黄频高清免费视频| 久久影院123| 老汉色∧v一级毛片| 美女高潮到喷水免费观看| 2018国产大陆天天弄谢| 亚洲欧美中文字幕日韩二区| 精品一品国产午夜福利视频| a级毛片在线看网站| 一级毛片电影观看| 美女国产高潮福利片在线看| 国产色婷婷99| 男女免费视频国产| 男女国产视频网站| 国产欧美日韩一区二区三区在线| 国产精品不卡视频一区二区| 制服诱惑二区| 国精品久久久久久国模美| 女人被躁到高潮嗷嗷叫费观| 亚洲人成77777在线视频| 成人国语在线视频| 亚洲精品日本国产第一区| 欧美变态另类bdsm刘玥| 国产黄色免费在线视频| 青春草亚洲视频在线观看| 国产野战对白在线观看| 国产麻豆69| 欧美成人精品欧美一级黄| 美女大奶头黄色视频| 夜夜骑夜夜射夜夜干| 九草在线视频观看| xxxhd国产人妻xxx| 国产成人一区二区在线| 精品国产乱码久久久久久男人| 天堂8中文在线网| 精品一品国产午夜福利视频| 国产精品免费大片| 国产xxxxx性猛交| 搡女人真爽免费视频火全软件| 欧美国产精品va在线观看不卡| 精品国产一区二区三区四区第35| 90打野战视频偷拍视频| 狂野欧美激情性bbbbbb| 国产亚洲一区二区精品| 男女高潮啪啪啪动态图| 亚洲婷婷狠狠爱综合网| 国产爽快片一区二区三区| 尾随美女入室| 18在线观看网站| 亚洲国产精品一区三区| 如何舔出高潮| 亚洲人成77777在线视频| 高清av免费在线| 一边摸一边做爽爽视频免费| 亚洲成av片中文字幕在线观看 | 夫妻午夜视频| 国产精品久久久久久久久免| 大陆偷拍与自拍| 欧美日韩亚洲国产一区二区在线观看 | 国产精品嫩草影院av在线观看| 亚洲经典国产精华液单| 精品国产国语对白av| 九草在线视频观看| 久久韩国三级中文字幕| 午夜福利在线免费观看网站| 欧美日韩国产mv在线观看视频| videosex国产| 午夜91福利影院| 久久鲁丝午夜福利片| 黑人欧美特级aaaaaa片| 亚洲男人天堂网一区| 搡老乐熟女国产| 国产午夜精品一二区理论片| 熟女少妇亚洲综合色aaa.| 天天躁狠狠躁夜夜躁狠狠躁| 各种免费的搞黄视频| av在线app专区| 极品人妻少妇av视频| 韩国精品一区二区三区| 亚洲国产av新网站| 国产乱人偷精品视频| 又大又黄又爽视频免费| 亚洲精品一区蜜桃| 精品一区二区免费观看| 免费黄频网站在线观看国产| xxx大片免费视频| 国产亚洲av片在线观看秒播厂| 搡女人真爽免费视频火全软件| 男人舔女人的私密视频| 五月伊人婷婷丁香| 韩国精品一区二区三区| 一级毛片 在线播放| 亚洲av在线观看美女高潮| 久久精品久久久久久噜噜老黄| 18禁动态无遮挡网站| 国产精品无大码| 另类精品久久| 欧美人与性动交α欧美精品济南到 | 亚洲精品一二三| 国产精品一区二区在线不卡| 深夜精品福利| 看免费av毛片| 伦理电影大哥的女人| av在线播放精品| 国产一区二区激情短视频 | 纵有疾风起免费观看全集完整版| 亚洲欧美一区二区三区国产| 免费高清在线观看视频在线观看| 欧美亚洲日本最大视频资源| 一区二区三区精品91| 精品国产超薄肉色丝袜足j| 国产精品国产三级专区第一集| 亚洲三区欧美一区| 五月伊人婷婷丁香| 亚洲精品一区蜜桃| 精品午夜福利在线看| 精品一区二区三卡| 国产视频首页在线观看| 亚洲精品美女久久久久99蜜臀 | 9色porny在线观看| www.av在线官网国产| 美女福利国产在线| 日本免费在线观看一区| 亚洲精品日本国产第一区| 交换朋友夫妻互换小说| 一级a爱视频在线免费观看| 婷婷色麻豆天堂久久| 黄网站色视频无遮挡免费观看| freevideosex欧美| 妹子高潮喷水视频| 成人国产麻豆网| 十八禁高潮呻吟视频| 久久久久久伊人网av| 少妇精品久久久久久久| 欧美国产精品一级二级三级| 欧美+日韩+精品| 啦啦啦在线免费观看视频4| 人人妻人人澡人人爽人人夜夜| 咕卡用的链子| 天天操日日干夜夜撸| 国产精品免费视频内射| 男人操女人黄网站| 久久久久久人妻| 免费播放大片免费观看视频在线观看| 精品人妻一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看 | 观看美女的网站| 黄色一级大片看看| 国产成人av激情在线播放| 不卡视频在线观看欧美| av卡一久久| 欧美激情高清一区二区三区 | 久热这里只有精品99| 国产精品一区二区在线不卡| 欧美激情极品国产一区二区三区| 嫩草影院入口| 欧美xxⅹ黑人| 麻豆乱淫一区二区| 久久久久精品久久久久真实原创| 韩国精品一区二区三区| 大码成人一级视频| 亚洲人成网站在线观看播放| 亚洲 欧美一区二区三区| 国产白丝娇喘喷水9色精品| videossex国产| 亚洲欧美一区二区三区黑人 | 欧美日韩一级在线毛片| 老女人水多毛片| 国产精品免费大片| 中文字幕人妻丝袜一区二区 | 精品人妻一区二区三区麻豆| 一边亲一边摸免费视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美精品高潮呻吟av久久| 国产精品免费大片| 亚洲在久久综合| 午夜日韩欧美国产| 欧美精品一区二区大全| 一区二区三区激情视频| 久久久久精品久久久久真实原创| 在线观看一区二区三区激情| 中国国产av一级| av有码第一页| 五月天丁香电影| 亚洲精品久久久久久婷婷小说| 如何舔出高潮| 久久久久视频综合| 亚洲欧美成人精品一区二区| 国产成人a∨麻豆精品| 国产精品av久久久久免费| 丝袜脚勾引网站| 欧美精品一区二区免费开放| 丝瓜视频免费看黄片| 在线观看三级黄色| 国产成人一区二区在线| 午夜91福利影院| 在线看a的网站| 建设人人有责人人尽责人人享有的| 91精品三级在线观看| 女人高潮潮喷娇喘18禁视频| 国产麻豆69| 黄色 视频免费看| 1024视频免费在线观看| 久久人人爽av亚洲精品天堂| 日韩在线高清观看一区二区三区| 欧美人与性动交α欧美精品济南到 | 大片电影免费在线观看免费| 三级国产精品片| 交换朋友夫妻互换小说| 亚洲国产欧美网| 亚洲,欧美,日韩| 99久久中文字幕三级久久日本| 亚洲少妇的诱惑av| 国产成人91sexporn| 人妻系列 视频| 亚洲av日韩在线播放| 在线 av 中文字幕| 免费人妻精品一区二区三区视频| 伊人久久国产一区二区| 男女高潮啪啪啪动态图| 久久久久久久久久久免费av| 女性被躁到高潮视频| 2022亚洲国产成人精品| 久久人人97超碰香蕉20202| 中国三级夫妇交换| 欧美激情 高清一区二区三区| 大片免费播放器 马上看| 久久人妻熟女aⅴ| 麻豆精品久久久久久蜜桃| 99热国产这里只有精品6| 欧美av亚洲av综合av国产av | 免费看av在线观看网站| 亚洲伊人久久精品综合| 久久毛片免费看一区二区三区| 亚洲色图综合在线观看| 高清黄色对白视频在线免费看| 国产精品国产三级国产专区5o| 国产精品秋霞免费鲁丝片| 国产xxxxx性猛交| 观看美女的网站| 国产精品av久久久久免费| 亚洲精品成人av观看孕妇| 亚洲欧美成人综合另类久久久| 天天影视国产精品| 精品少妇久久久久久888优播| 国产av精品麻豆| 秋霞在线观看毛片| av片东京热男人的天堂| 久久毛片免费看一区二区三区| 男男h啪啪无遮挡| 精品人妻偷拍中文字幕| 三上悠亚av全集在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美一区二区三区久久| 天天躁夜夜躁狠狠久久av| 18禁动态无遮挡网站| 又粗又硬又长又爽又黄的视频| 人人妻人人爽人人添夜夜欢视频| 我要看黄色一级片免费的| 国产精品成人在线| 亚洲欧美色中文字幕在线| 亚洲av电影在线观看一区二区三区| 国产免费现黄频在线看| 少妇 在线观看| 日韩精品有码人妻一区| 久久久精品免费免费高清| 男男h啪啪无遮挡| 91午夜精品亚洲一区二区三区| 狠狠精品人妻久久久久久综合| 欧美精品国产亚洲| 黄色一级大片看看| 最近2019中文字幕mv第一页| 久久久精品国产亚洲av高清涩受| 男女无遮挡免费网站观看| xxxhd国产人妻xxx| 欧美精品一区二区大全| 欧美成人午夜免费资源| 亚洲欧洲精品一区二区精品久久久 | 亚洲精华国产精华液的使用体验| 毛片一级片免费看久久久久| 午夜日本视频在线| 国产乱人偷精品视频| 国产日韩一区二区三区精品不卡| 美女脱内裤让男人舔精品视频| 少妇熟女欧美另类| 18禁裸乳无遮挡动漫免费视频| 亚洲久久久国产精品| 免费在线观看视频国产中文字幕亚洲 | 老司机亚洲免费影院| 午夜av观看不卡| 亚洲三级黄色毛片| 亚洲,一卡二卡三卡| 午夜福利网站1000一区二区三区| 亚洲国产成人一精品久久久| 97在线视频观看| 色94色欧美一区二区| 一边摸一边做爽爽视频免费| 国产成人91sexporn| 人人妻人人爽人人添夜夜欢视频| 99久国产av精品国产电影| 哪个播放器可以免费观看大片| 久久婷婷青草| 久久精品国产综合久久久| 免费高清在线观看视频在线观看| 国产精品麻豆人妻色哟哟久久| 欧美日韩视频高清一区二区三区二| 久久99一区二区三区| 自线自在国产av| 亚洲av电影在线观看一区二区三区| 日本-黄色视频高清免费观看| 国产精品蜜桃在线观看| 久久久精品94久久精品| 一级毛片黄色毛片免费观看视频| 妹子高潮喷水视频| 免费久久久久久久精品成人欧美视频| 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| 在线观看三级黄色| 久久精品人人爽人人爽视色| 日韩中文字幕欧美一区二区 | 精品久久久久久电影网| a级片在线免费高清观看视频| 夜夜骑夜夜射夜夜干| 中文字幕人妻熟女乱码| 国产又色又爽无遮挡免| av不卡在线播放| 五月伊人婷婷丁香| 久久毛片免费看一区二区三区| 中文字幕精品免费在线观看视频| 久久精品国产自在天天线| 考比视频在线观看| 国产精品久久久久成人av| 日韩伦理黄色片| 久久亚洲国产成人精品v| 久久久久网色| 亚洲综合色惰| 欧美日本中文国产一区发布| 免费黄色在线免费观看| av又黄又爽大尺度在线免费看| 岛国毛片在线播放| 国产精品国产三级国产专区5o| 午夜激情久久久久久久| 久久精品人人爽人人爽视色| 亚洲精品在线美女| 免费av中文字幕在线| 久久人妻熟女aⅴ|