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

    基于EMD和PCA 的滾動軸承故障信號特征提取研究*

    2015-11-18 12:27:54楊魯偉楊俊玲
    機(jī)電工程 2015年10期
    關(guān)鍵詞:外圈特征值維度

    張 穎,馬 波*,張 明,楊魯偉,楊俊玲

    (1.北京化工大學(xué) 診斷與自愈工程研究中心,北京 100029;2.中國科學(xué)院 低溫工程學(xué)重點(diǎn)實(shí)驗(yàn)室(理化技術(shù)研究所),北京 100190)

    0 引言

    滾動軸承是離心泵重要部件,在持續(xù)不斷的交變載荷作用下,極易發(fā)生故障,一旦發(fā)生故障,輕則引起停機(jī)影響生產(chǎn),重則發(fā)生爆炸事故。因此,實(shí)現(xiàn)離心泵滾動軸承故障診斷具有重要意義。

    滾動軸承發(fā)生故障后,經(jīng)過故障點(diǎn)時產(chǎn)生瞬時沖擊,使信號呈現(xiàn)非平穩(wěn)性[1]。針對非平穩(wěn)信號特征提取的問題,國內(nèi)外進(jìn)行了大量的研究。文獻(xiàn)[2]運(yùn)用小波包變換提取滾動軸承故障信號特征,并和神經(jīng)網(wǎng)絡(luò)方法相結(jié)合,實(shí)現(xiàn)滾動軸承故障診斷。文獻(xiàn)[3]將小波變換和能量算子結(jié)合提取氣閥故障特征,實(shí)驗(yàn)證明該方法能較好地區(qū)分氣閥不同故障,并表現(xiàn)出抑噪特性。文獻(xiàn)[4]將連續(xù)小波變換和決策樹分類方法結(jié)合實(shí)現(xiàn)了電動離心泵的故障診斷。文獻(xiàn)[5]通過分析對比多種時頻分析方法,驗(yàn)證Hilbert 變換提取非平穩(wěn)信號上的優(yōu)勢?,F(xiàn)有的大多數(shù)非平穩(wěn)信號處理方法,如小波變換、小波包變換等,在處理非平穩(wěn)信號時,由于受到海森伯測不準(zhǔn)原理約束,具有局限性。

    1998年,Hilbert-Huang 提出的一種新的時頻分析方法—經(jīng)驗(yàn)?zāi)B(tài)分解,從原理上分析該方法不受測不準(zhǔn)原理約束,且該方法具有自適應(yīng)性,使非平穩(wěn)信號的特征提取不受限于窗函數(shù)的選取,已經(jīng)在很多方面已經(jīng)得到廣泛的應(yīng)用。文獻(xiàn)[6]運(yùn)用經(jīng)驗(yàn)?zāi)B(tài)分解方法(empirical mode decomposition,EMD)和AR 模型提取滾動軸承故障信號特征,并通過試驗(yàn)驗(yàn)證了該方法的有效性。文獻(xiàn)[7]運(yùn)用EMD 提取故障信號特征,神經(jīng)網(wǎng)絡(luò)對故障進(jìn)行診斷,試驗(yàn)驗(yàn)證了該方法的可有效實(shí)現(xiàn)滾動軸承智能故障診斷。文獻(xiàn)[8-10]提出基于EMD和峭度的Hilbert 包絡(luò)解調(diào)方法,用于診斷滾動軸承故障,并驗(yàn)證了該方法的有效性。但是,這些方法在提取故障特征時,得到的特征值維度較高,其中包含很多無效的特征值,為了去除冗余的特征值,降低計算復(fù)雜度,提取有效特征,提高診斷效率,可選用主成分分析方法(principal component analysis,PCA)方法解決上述問題[11-12]。

    本研究擬將EMD和PCA 方法相結(jié)合,運(yùn)用EMD方法去除信號中的噪聲,提取出有用的故障信號特征,PCA 方法去除冗雜的特征值,并提取出主成分,以提高診斷效率,降低故障診斷計算復(fù)雜度。

    1 EMD 方法

    經(jīng)驗(yàn)?zāi)J椒纸夥椒ㄕJ(rèn)為,任何信號由不同的固有模式、簡單的振蕩信號組成。線性和非平穩(wěn)模式的零交叉點(diǎn)和極值點(diǎn)的數(shù)量都是相同的。EMD 分解的目的就是解決瞬時頻率的問題。

    任何一個獲取的特征信號都可以分解為多個固有模態(tài)函數(shù)(IMF),同時每個IMF 也有自己的要求,必須滿足以下幾個定義[13]:

    (1)在整個實(shí)數(shù)域或復(fù)數(shù)域中,極值與零交叉點(diǎn)的數(shù)量頂多相差一個。

    (2)在函數(shù)上的任一點(diǎn),信號的上、下包絡(luò)線所構(gòu)成的平均包絡(luò)值必須為零。

    上、下包絡(luò)平均值為m1,信號x(t)和中間的差異m1是第一部分:

    理論上,h1是x(t)的第一部分,如果h1是固有模態(tài)函數(shù)。

    如果h1不是,那么有:

    停止篩選直到hk是本征模量時,有:

    并可表示為:

    從x(t)信號中分離c1,可得:

    再將r1作為原始信號,重復(fù)上述過程,得出第2個本征量。依此類推,可得:

    這里,假設(shè)原始信號為x(t)。通過上式綜合得出:

    式中:imfi(t)—x(t)經(jīng)驗(yàn)?zāi)B(tài)分解得到的第i個固有模態(tài)函數(shù)(IMF);rn(t)—去除n個固有模態(tài)函數(shù)(IMF)后所剩的剩余信號。

    直到r(t)極限值小于兩個或者|r(t)|很小時可結(jié)束分解。

    式(7,8)以上過程體現(xiàn)了EMD 分解的完備性。

    2 PCA 方法

    主成分分析(PCA)主要用于縮減數(shù)據(jù)的維度。一般情況下,如果系統(tǒng)有m個變量,就可以用k個不相關(guān)主成分(k <m)來表示樣本之間的差異,然后就可以把原m個變量的變化用k個主成分來代替。這就將原數(shù)據(jù)空間分成了兩個完全不相關(guān)的空間,一個是由k個主成分構(gòu)成的主元子空間,另一個是由其他剩余變量構(gòu)成的殘差子空間,主元子空間用來描述整個系統(tǒng)變化在理想的情況下,而隨機(jī)噪聲包含在殘差子空間[14]。

    分解樣本數(shù)據(jù)矩陣X:

    式中:UUT=I,VVT=I;σi(i=1,2,…,m)—矩陣X 的奇異值,且σ1>σ2>…σm。從而可得:

    本研究進(jìn)行維度縮減的主要步驟如下:

    (1)從非平穩(wěn)信號中提取n個特征值Xi(i=1,2,3,…,n),特征值Xi維度為N{Xi1,Xi2,…,Xij,…,XiN}(j=1,2,…,N)

    (2)為特征值建立樣本矩陣M:

    (3)根據(jù)樣本矩陣M,求協(xié)方差矩陣;

    (4)計算協(xié)方差矩陣M'的特征值和特征向量。

    (5)計算協(xié)方差矩陣M'特征值的絕對值并按值得大小進(jìn)行降序排列,同時調(diào)整相應(yīng)特征向量的順序,得到矩陣P。

    (6)計算特征矩陣:)

    (7)計算特征值得累計貢獻(xiàn)率,根據(jù)累計貢獻(xiàn)率判定特征矩陣K 的前n 行即為對多特征參數(shù)處理后最能代表故障征兆的n個特征參數(shù)。

    3 基于EMD和PCA 的特征提取方法

    本研究采取EMD和PCA 相結(jié)合的方法,提取滾動軸承典型的故障特征的流程如圖1 所示。

    圖1 基于EMD和PCA 的特征提取方法流程

    (1)用傳感器采集滾動軸承和軸系故障的加速度信號,隨機(jī)選取各狀態(tài)下的樣本信號作為分析信號。

    (2)對分析信號分別進(jìn)行EMD 分解,得到多層IMF 函數(shù)和一層剩余函數(shù)r(c1,c2,c3,c4…cn),不同的IMF 分量包含著不同的原始信號的部分信息并且頻率分布由高到低,因此,信號能夠在不同的分辨率下準(zhǔn)確地顯示出來。

    (3)對每一個IMF 分量進(jìn)行傅里葉變換,得到每一層IMF 函數(shù)的頻譜圖,頻譜圖能夠反映信號中的不同頻率成分的能量。

    (4)計算每一層IMF 的能量值,根據(jù)計算選擇大部分能量所在的前n 層IMF 函數(shù),根據(jù)頻率成分分布的不同把前n 層IMF 函數(shù)對應(yīng)的各層頻譜圖分成M段,對每一段的能量值進(jìn)行計算,得到一個N×M 的高維度特征向量。能量值計算,如下式所示:

    式中:E—能量值,n—IMF 的頻率值數(shù),A—幅值。

    (5)為了減小計算負(fù)擔(dān),找出最能代表故障特征的特征值,本研究運(yùn)用PCA 方法對上述特征向量進(jìn)行維度縮減,針對高維度的故障特征的特征值建立樣本矩陣,對樣本矩陣的特征值和特征向量分別進(jìn)行計算,如果所得到的前n個特征的累計貢獻(xiàn)率達(dá)到85%,即可認(rèn)為該組特征的前n個特征具有代表性,可代表整組的特征。

    4 實(shí)驗(yàn)及結(jié)果分析

    本次試驗(yàn)所用的試驗(yàn)臺如圖2 所示。其包含了兩個支撐座,雙支承轉(zhuǎn)子試驗(yàn)臺,左端為電機(jī)。傳感器布置如圖3 所示。

    圖2 滾動軸承故障試驗(yàn)臺

    圖3 傳感器布置

    該次試驗(yàn)所用軸承型號為NU205EM,軸承轉(zhuǎn)速為1 200 r/min,采樣頻率為25.6 kHz,每組故障的采樣點(diǎn)數(shù)為450 000個左右,該實(shí)驗(yàn)截取其中的5 000個點(diǎn)進(jìn)行分析,轉(zhuǎn)速為1 200 r/min。該實(shí)驗(yàn)所設(shè)計的試驗(yàn)故障為滾動軸承外圈故障、滾動軸承內(nèi)圈故障、滾動軸承滾動體故障。

    筆者以滾動軸承為外圈故障為例,說明EMD 分解和特征值計算過程。運(yùn)用EMD 方法將滾動軸承外圈故障信號分解為11個IMF 分量和一個剩余函數(shù)如圖4 所示。

    本研究對上述11 層IMF 函數(shù)進(jìn)行傅里葉變換之后得到頻譜如圖5 所示。計算11個IMF 分量的能量值及其占總能量的百分比結(jié)果如表1 所示。

    表1 前5 層IMF 能量及占總能量的百分比

    圖4 滾動軸承外圈故障原始波形及EMD 分解結(jié)果

    通過表1 可以看出前5 層的IMF 函數(shù)能量占總能量的97.2%,能量主要集中在前5 層IMF 函數(shù)中,因此筆者選取經(jīng)過傅里葉變換的前5 層IMF 函數(shù)計算特征值。從頻譜圖中可以看出能量主要分布在0~12 000 Hz頻段,將每層IMF 數(shù)據(jù)0~12 000 Hz 頻段均分成6 段,即0~2 000 Hz,2 000 Hz~4 000 Hz,4 000 Hz~6 000 Hz,6 000 Hz~8 000 Hz,8 000 Hz~10 000 Hz,10 000 Hz~12 000 Hz,求每一段的能量值,每層IMF提取6個特征值,每組故障提取6×5個特征值,滾動軸承外圈故障特征值如圖6 所示。

    圖5 IMF 分量及剩余函數(shù)的傅變頻譜

    本研究通過計算得到30個特征值,計算每組故障得到的每個特征值的累計貢獻(xiàn)率,根據(jù)累計貢獻(xiàn)率,選擇最能代表該組故障特征的特征值進(jìn)行后序計算。滾動軸承外圈故障得到的特征值及其相應(yīng)的累計貢獻(xiàn)率計算結(jié)果如表2 所示。

    圖6 滾動軸承外圈故障特征值

    表2 滾動軸承外圈故障的特征值及累計貢獻(xiàn)率

    一般認(rèn)為,累計貢獻(xiàn)率達(dá)到85%左右的前n個特征值即可作為代表本組故障的特征值,從表2 可以看出,該組特征值的前8個特征的累計貢獻(xiàn)率達(dá)到85%,即前8個特征值可作為代表滾動軸承外圈故障縮減后的特征值,PCA 方法有效地實(shí)現(xiàn)了特征值維度縮減。

    本研究嘗試將每組故障的前3 維特征表示在3 維坐標(biāo)系中,結(jié)果如圖7 所示。從圖7 中可以看出,同一種故障的特征值點(diǎn)聚集在一起,不同故障的前三維特征值被清晰分開。由此可以得出結(jié)論,前三維特征數(shù)據(jù)即可將不同的故障類型進(jìn)行區(qū)分,PCA 方法能去除冗余特征值,有效實(shí)現(xiàn)故障特征維度縮減。

    圖7 4 種工況特征值維度縮減后在三維坐標(biāo)中的表示

    5 結(jié)束語

    本研究針對滾動軸承典型故障信號的非平穩(wěn)性,運(yùn)用EMD 方法提取信號特征,然后運(yùn)用PCA 方法去除冗余的特征值。主要結(jié)論如下;

    (1)EMD 方法分解非平穩(wěn)信號,可以得到多個頻率由高至低分布的IMF 分量,選擇能量集中的IMF,可有效去除低頻噪聲干擾,提取故障信號的高頻分量,計算得到信號特征值。

    (2)PCA 方法可依據(jù)累計貢獻(xiàn)率,選擇少數(shù)特征值,來代替高維的特征值,該方法能夠保留較完整的故障特征信息,同時PCA 對滾動軸承典型故障特征分類有較好的效果。

    [1]GELMAN L,MURRAY B,PATEL T H,et al.Vibration diagnostics of rolling bearings by novel nonlinear non-stationary wavelet bicoherence technology[J].Engineering Structures,2014(80):514-520.

    [2]劉樂平,林鳳濤.基于小波包特征向量與神經(jīng)網(wǎng)絡(luò)的滾動軸承故障診斷[J].軸承,2008(4):46-48.

    [3]魏中青,馬 波,么子云,等.運(yùn)用小波包變換與能量算子的氣閥故障特征提?。跩].振動.測試與診斷,2011,31(1):50-54.

    [4]MURALIDHARAN,SUGUMARAN V.Feature extraction using wavelets and classification through decision tree algorithm for fault diagnosis of mono-block centrifugalpump[J].Measurement,2013,46(1):353-359

    [5]向 玲,唐貴基,胡愛軍.旋轉(zhuǎn)機(jī)械非平穩(wěn)振動信號的時頻分析比較[J].振動與沖擊,2010(2):42-45.

    [6]JUNSHENG C,DEJIE Y,YU Y.A fault diagnosis approach for roller bearings based on EMD method and AR model[J].Mechanical Systems and Signal Processing,2006,20(2):350-362.

    [7]ALI J B,F(xiàn)NAIECH N,SAIDI L,et al.Application of empirical mode decomposition and artificial neural network for automatic bearing fault diagnosis based on vibration signals[J].Applied Acoustics,2015(89):16-27.

    [8]周 浩,賈民平.基于EMD和峭度的Hilbert 包絡(luò)解調(diào)在滾動軸承故障診斷中的應(yīng)用分析[J].機(jī)電工程,2014,31(9):1136-1139.

    [9]趙月剛.4N6Z 凝結(jié)水泵軸承故障分析及處理[J].流體機(jī)械,2011,39(10):55-56.

    [10]王逍夢,任 章.基于小波變換的二階傳熱感器故障診斷研究[J].兵工自動化,2013(5):78-82.

    [11]SHAO R,HU W,WANG Y.The fault feature extraction and classification of gear using principal component analysis and kernel principal component analysis based on the wavelet packet transform[J].Measurement,2014,54(10):118-132.

    [12]ALVAREZ D G,F(xiàn)UENTE M J,SAINZ G I.Faultdetection and isolation in transient states using principal component analysis[J].Journal of Process Control,2012,22(3):551-563.

    [13]YU D,CHENG J,YANG Y.Application of EMD method and Hilbert spectrum to the fault diagnosis of roller bearings[J].Mechanical systems and signal processing,2005,19(2):259-270.

    [14]劉永斌,何清波,孔凡讓,等.基于PCA和SVM 的內(nèi)燃機(jī)故障診斷[J].振動.測試與診斷,2012,32(2):250-255.

    猜你喜歡
    外圈特征值維度
    深溝球軸承外圈表面凹坑缺陷分析
    哈爾濱軸承(2022年1期)2022-05-23 13:13:16
    一類帶強(qiáng)制位勢的p-Laplace特征值問題
    單圈圖關(guān)聯(lián)矩陣的特征值
    角接觸球軸承外圈鎖口高度自動檢測規(guī)改進(jìn)
    哈爾濱軸承(2020年2期)2020-11-06 09:22:34
    淺論詩中“史”識的四個維度
    中華詩詞(2019年7期)2019-11-25 01:43:00
    光的維度
    燈與照明(2016年4期)2016-06-05 09:01:45
    “五個維度”解有機(jī)化學(xué)推斷題
    基于商奇異值分解的一類二次特征值反問題
    關(guān)于兩個M-矩陣Hadamard積的特征值的新估計
    人生三維度
    吐魯番(2014年2期)2014-02-28 16:54:43
    中文欧美无线码| 国产精品国产三级国产av玫瑰| 国产亚洲最大av| 18禁在线无遮挡免费观看视频| 国产中年淑女户外野战色| 老师上课跳d突然被开到最大视频| 91狼人影院| 在线观看av片永久免费下载| 亚洲av.av天堂| 一二三四中文在线观看免费高清| 欧美丝袜亚洲另类| 精品久久久久久久人妻蜜臀av| 全区人妻精品视频| 欧美bdsm另类| 亚洲婷婷狠狠爱综合网| 别揉我奶头 嗯啊视频| 国产免费男女视频| 久久国产乱子免费精品| av视频在线观看入口| 麻豆精品久久久久久蜜桃| 成人欧美大片| 永久网站在线| 少妇的逼好多水| 男插女下体视频免费在线播放| 欧美高清性xxxxhd video| 亚洲在久久综合| 久久久午夜欧美精品| 久久精品影院6| 午夜福利视频1000在线观看| 高清在线视频一区二区三区 | 国产三级在线视频| 九九热线精品视视频播放| 精品一区二区三区视频在线| 99热这里只有是精品50| 午夜福利网站1000一区二区三区| 黑人高潮一二区| 一边摸一边抽搐一进一小说| 成年av动漫网址| 日韩,欧美,国产一区二区三区 | 直男gayav资源| 黄色配什么色好看| 亚洲伊人久久精品综合 | 天堂av国产一区二区熟女人妻| 国产午夜福利久久久久久| 99热网站在线观看| 国产精品嫩草影院av在线观看| videos熟女内射| 在线a可以看的网站| 国内揄拍国产精品人妻在线| 国产黄片视频在线免费观看| 国产精品伦人一区二区| 国产精品.久久久| 日韩中字成人| 好男人视频免费观看在线| 亚洲久久久久久中文字幕| www.av在线官网国产| 美女大奶头视频| 青青草视频在线视频观看| 桃色一区二区三区在线观看| 久久久欧美国产精品| 日韩av不卡免费在线播放| 欧美一区二区国产精品久久精品| 简卡轻食公司| 国产探花在线观看一区二区| 国产免费一级a男人的天堂| 亚州av有码| 最后的刺客免费高清国语| videos熟女内射| 在线天堂最新版资源| 99久久中文字幕三级久久日本| 国产高清有码在线观看视频| 久久草成人影院| av在线观看视频网站免费| 国产精品一区www在线观看| 精品久久久噜噜| 免费电影在线观看免费观看| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 美女脱内裤让男人舔精品视频| 久久精品国产鲁丝片午夜精品| 91精品伊人久久大香线蕉| 搡女人真爽免费视频火全软件| 亚洲中文字幕日韩| 精品国产三级普通话版| 亚洲aⅴ乱码一区二区在线播放| 国产美女午夜福利| 精品99又大又爽又粗少妇毛片| 国产免费男女视频| 日韩精品有码人妻一区| 国模一区二区三区四区视频| 亚洲欧美精品综合久久99| av在线亚洲专区| 亚洲精品乱码久久久久久按摩| 九九在线视频观看精品| 蜜桃亚洲精品一区二区三区| 深爱激情五月婷婷| 欧美日本亚洲视频在线播放| 亚洲av中文字字幕乱码综合| 女人十人毛片免费观看3o分钟| 国产伦一二天堂av在线观看| 欧美性猛交╳xxx乱大交人| 午夜福利网站1000一区二区三区| 亚洲国产精品久久男人天堂| 青春草亚洲视频在线观看| 永久网站在线| 国产高清三级在线| 亚洲国产精品成人久久小说| 欧美一级a爱片免费观看看| 在线天堂最新版资源| 午夜a级毛片| av在线播放精品| 国产一级毛片在线| av福利片在线观看| 亚洲av免费在线观看| 日本一二三区视频观看| 国产精品野战在线观看| 久久99蜜桃精品久久| 欧美一区二区亚洲| 精品久久久久久成人av| 99视频精品全部免费 在线| 婷婷六月久久综合丁香| 99久久精品国产国产毛片| 亚洲av成人精品一区久久| 亚洲中文字幕日韩| 久久久精品94久久精品| 男人狂女人下面高潮的视频| 伊人久久精品亚洲午夜| 国产成人91sexporn| 亚洲五月天丁香| 久久久久久久久久久丰满| 在线观看66精品国产| 免费播放大片免费观看视频在线观看 | 欧美xxxx黑人xx丫x性爽| 老司机福利观看| 欧美激情久久久久久爽电影| 国内揄拍国产精品人妻在线| 天天躁日日操中文字幕| 成人美女网站在线观看视频| 午夜福利网站1000一区二区三区| 久久久久久久午夜电影| 嘟嘟电影网在线观看| 日本熟妇午夜| 日韩大片免费观看网站 | 国产单亲对白刺激| av在线亚洲专区| 熟女电影av网| 日韩强制内射视频| 亚洲av电影不卡..在线观看| 亚洲成色77777| 亚洲一级一片aⅴ在线观看| 美女cb高潮喷水在线观看| 日韩一区二区视频免费看| 蜜臀久久99精品久久宅男| 免费观看人在逋| 免费黄网站久久成人精品| 欧美区成人在线视频| 欧美精品国产亚洲| 亚洲精品456在线播放app| 国产白丝娇喘喷水9色精品| 国产精品一及| a级毛色黄片| 久久久久国产网址| 国产亚洲91精品色在线| 亚洲精品456在线播放app| 在线a可以看的网站| 精品午夜福利在线看| 少妇熟女欧美另类| 欧美区成人在线视频| 国产黄a三级三级三级人| 国产乱来视频区| 九草在线视频观看| 高清av免费在线| 99热这里只有是精品在线观看| 亚洲国产精品合色在线| 亚洲内射少妇av| 中文字幕精品亚洲无线码一区| 久久久久久久久久成人| 亚洲,欧美,日韩| 两个人视频免费观看高清| 成人亚洲精品av一区二区| 全区人妻精品视频| 免费观看人在逋| 级片在线观看| 精品人妻视频免费看| 午夜福利在线观看吧| 中文字幕久久专区| 少妇被粗大猛烈的视频| 午夜爱爱视频在线播放| 老司机影院毛片| 久久草成人影院| 欧美激情国产日韩精品一区| 亚洲av男天堂| 女人被狂操c到高潮| 一级黄色大片毛片| 日本猛色少妇xxxxx猛交久久| 久久精品国产亚洲网站| 免费搜索国产男女视频| 看片在线看免费视频| 久久国产乱子免费精品| 午夜爱爱视频在线播放| 欧美日韩综合久久久久久| 一边亲一边摸免费视频| 九九爱精品视频在线观看| 三级国产精品欧美在线观看| 99九九线精品视频在线观看视频| 在线播放无遮挡| 美女大奶头视频| 春色校园在线视频观看| 国产伦一二天堂av在线观看| 99久久中文字幕三级久久日本| 国产美女午夜福利| 欧美日韩一区二区视频在线观看视频在线 | 亚洲成av人片在线播放无| 精品国产三级普通话版| 韩国高清视频一区二区三区| 精品不卡国产一区二区三区| 日韩强制内射视频| 成人综合一区亚洲| 亚洲美女视频黄频| 亚洲一区高清亚洲精品| 国产女主播在线喷水免费视频网站 | 99久久无色码亚洲精品果冻| 国产色婷婷99| 日本黄色视频三级网站网址| 精品久久久久久电影网 | 插逼视频在线观看| 日韩制服骚丝袜av| 丰满少妇做爰视频| 直男gayav资源| 国产精品一区二区三区四区免费观看| 村上凉子中文字幕在线| 麻豆国产97在线/欧美| 国产真实乱freesex| 又爽又黄a免费视频| 日韩,欧美,国产一区二区三区 | 久久午夜福利片| 最近的中文字幕免费完整| 又黄又爽又刺激的免费视频.| 欧美最新免费一区二区三区| 人妻夜夜爽99麻豆av| 三级国产精品片| 成人一区二区视频在线观看| 波野结衣二区三区在线| 精品久久久久久电影网 | 在线观看66精品国产| 亚洲欧洲国产日韩| 女人十人毛片免费观看3o分钟| 亚洲国产精品久久男人天堂| 国产亚洲5aaaaa淫片| 97热精品久久久久久| 中国美白少妇内射xxxbb| 日韩大片免费观看网站 | 波多野结衣高清无吗| 精品久久久久久久末码| 欧美激情国产日韩精品一区| 高清av免费在线| 成年av动漫网址| 91午夜精品亚洲一区二区三区| 国产精品久久视频播放| 国产精品乱码一区二三区的特点| 亚洲av成人av| 精品人妻一区二区三区麻豆| 久久99蜜桃精品久久| 欧美色视频一区免费| 国产精品一区二区三区四区免费观看| 免费看美女性在线毛片视频| 亚洲三级黄色毛片| 黑人高潮一二区| 欧美成人免费av一区二区三区| 亚洲精品aⅴ在线观看| 嫩草影院入口| 国产精品一区二区三区四区免费观看| 午夜亚洲福利在线播放| 中文字幕av在线有码专区| 国产成人91sexporn| 啦啦啦啦在线视频资源| 日韩一本色道免费dvd| 日韩欧美精品免费久久| 亚洲av二区三区四区| 久久久精品欧美日韩精品| 国产成人a区在线观看| 少妇丰满av| 少妇人妻精品综合一区二区| 国产午夜福利久久久久久| 岛国在线免费视频观看| 日日摸夜夜添夜夜添av毛片| 亚洲18禁久久av| 免费黄网站久久成人精品| 国产av在哪里看| 久久久色成人| 亚洲第一区二区三区不卡| 91精品伊人久久大香线蕉| 色5月婷婷丁香| 99在线人妻在线中文字幕| 别揉我奶头 嗯啊视频| 亚洲精品,欧美精品| 国产精品嫩草影院av在线观看| 岛国毛片在线播放| 搡老妇女老女人老熟妇| 中文资源天堂在线| 一级毛片电影观看 | 26uuu在线亚洲综合色| 麻豆乱淫一区二区| 久久久亚洲精品成人影院| 一本久久精品| 国产av码专区亚洲av| 91在线精品国自产拍蜜月| 看黄色毛片网站| 亚洲自偷自拍三级| 国产日韩欧美在线精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久6这里有精品| 国产精品人妻久久久久久| 九色成人免费人妻av| 久久久a久久爽久久v久久| 国产精品乱码一区二三区的特点| 麻豆国产97在线/欧美| 国产高清三级在线| 免费不卡的大黄色大毛片视频在线观看 | 最近最新中文字幕大全电影3| 婷婷色综合大香蕉| 女的被弄到高潮叫床怎么办| 国产三级在线视频| 成年免费大片在线观看| 亚洲真实伦在线观看| 国产片特级美女逼逼视频| 人体艺术视频欧美日本| 日韩欧美精品免费久久| 精品国产露脸久久av麻豆 | 国产在视频线精品| 国语自产精品视频在线第100页| 超碰av人人做人人爽久久| 久久久欧美国产精品| АⅤ资源中文在线天堂| 联通29元200g的流量卡| 国产成人精品婷婷| 国产免费又黄又爽又色| 久久精品影院6| 老女人水多毛片| 国内揄拍国产精品人妻在线| 午夜福利成人在线免费观看| 精华霜和精华液先用哪个| 建设人人有责人人尽责人人享有的 | 国产黄a三级三级三级人| 国产成人a区在线观看| 国产精品一及| 亚洲人成网站高清观看| 亚洲熟妇中文字幕五十中出| 国产综合懂色| 丰满人妻一区二区三区视频av| 蜜桃久久精品国产亚洲av| 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 久久99精品国语久久久| a级毛色黄片| 波多野结衣高清无吗| 国产黄色视频一区二区在线观看 | 99久国产av精品| 亚洲国产精品sss在线观看| 三级国产精品片| 汤姆久久久久久久影院中文字幕 | 国产久久久一区二区三区| 国产一区二区三区av在线| 久久精品久久久久久久性| 国产三级中文精品| 少妇熟女欧美另类| 国产午夜精品久久久久久一区二区三区| 最近的中文字幕免费完整| 两个人的视频大全免费| 亚洲成人久久爱视频| 麻豆久久精品国产亚洲av| 男人和女人高潮做爰伦理| 国产精品久久久久久久电影| 大香蕉97超碰在线| 三级国产精品欧美在线观看| 国产淫片久久久久久久久| 身体一侧抽搐| 简卡轻食公司| 又粗又爽又猛毛片免费看| 国产成人aa在线观看| 久久精品人妻少妇| 久久精品久久久久久久性| 黄片无遮挡物在线观看| 国产一级毛片在线| 搡女人真爽免费视频火全软件| 亚洲不卡免费看| 国内揄拍国产精品人妻在线| 99久久人妻综合| 国产精品av视频在线免费观看| 久久久久久久国产电影| 国产精品久久久久久精品电影| 久久精品国产亚洲av天美| 国产精品av视频在线免费观看| 国产精品永久免费网站| 黑人高潮一二区| 啦啦啦观看免费观看视频高清| 亚洲熟妇中文字幕五十中出| 看片在线看免费视频| 久久久精品94久久精品| 国产免费男女视频| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久成人| 日本免费a在线| 国产毛片a区久久久久| 99久久无色码亚洲精品果冻| 九九在线视频观看精品| 亚洲四区av| 午夜福利在线观看免费完整高清在| 免费观看的影片在线观看| 国产av码专区亚洲av| 国产精品熟女久久久久浪| 欧美一区二区精品小视频在线| 国产精品美女特级片免费视频播放器| 神马国产精品三级电影在线观看| 干丝袜人妻中文字幕| 3wmmmm亚洲av在线观看| 直男gayav资源| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩高清专用| 神马国产精品三级电影在线观看| 国产单亲对白刺激| 亚洲伊人久久精品综合 | 久久99热这里只有精品18| 级片在线观看| 麻豆国产97在线/欧美| 日韩在线高清观看一区二区三区| 最近中文字幕高清免费大全6| 国产毛片a区久久久久| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 午夜爱爱视频在线播放| 91aial.com中文字幕在线观看| 特大巨黑吊av在线直播| 精品一区二区免费观看| 啦啦啦韩国在线观看视频| 嫩草影院精品99| 久久精品熟女亚洲av麻豆精品 | 久久午夜福利片| 国产午夜福利久久久久久| 午夜爱爱视频在线播放| 国产一区二区三区av在线| 成年av动漫网址| 最近最新中文字幕免费大全7| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| kizo精华| 国产成人a∨麻豆精品| 国产伦一二天堂av在线观看| 中文字幕熟女人妻在线| 久久久久久久久久成人| 成人无遮挡网站| 少妇熟女aⅴ在线视频| 乱人视频在线观看| 久久精品国产亚洲av涩爱| 亚洲aⅴ乱码一区二区在线播放| 色网站视频免费| 国产精品熟女久久久久浪| 日韩精品青青久久久久久| 国产精品久久久久久久久免| 国产成人91sexporn| 大又大粗又爽又黄少妇毛片口| 狂野欧美激情性xxxx在线观看| 特级一级黄色大片| 日韩欧美精品免费久久| 日韩av不卡免费在线播放| 久久久国产成人精品二区| 亚洲美女视频黄频| 免费人成在线观看视频色| 中文字幕熟女人妻在线| h日本视频在线播放| 一级毛片我不卡| 欧美精品一区二区大全| 91精品国产九色| 国产精品一及| 免费av观看视频| 欧美激情在线99| 插逼视频在线观看| 热99re8久久精品国产| 亚洲精品乱久久久久久| a级一级毛片免费在线观看| 我要搜黄色片| 全区人妻精品视频| 麻豆乱淫一区二区| 精品久久久久久电影网 | 岛国在线免费视频观看| 午夜福利在线观看免费完整高清在| 非洲黑人性xxxx精品又粗又长| 黄片wwwwww| 国产在视频线在精品| 久久精品影院6| 伦精品一区二区三区| 卡戴珊不雅视频在线播放| 国产老妇女一区| 精品人妻熟女av久视频| 欧美日韩精品成人综合77777| 级片在线观看| or卡值多少钱| 三级男女做爰猛烈吃奶摸视频| 国产精品美女特级片免费视频播放器| 国产一区二区在线观看日韩| 中文字幕久久专区| 国产黄片视频在线免费观看| 国产一级毛片七仙女欲春2| 国产精品一二三区在线看| 久久久成人免费电影| 国产色婷婷99| 免费观看性生交大片5| 国产极品精品免费视频能看的| 色视频www国产| 搡女人真爽免费视频火全软件| 一级毛片aaaaaa免费看小| 日本爱情动作片www.在线观看| 日本免费在线观看一区| 美女高潮的动态| 草草在线视频免费看| 夜夜爽夜夜爽视频| 国产精品精品国产色婷婷| 欧美最新免费一区二区三区| 国产爱豆传媒在线观看| 男人和女人高潮做爰伦理| 亚洲欧美日韩无卡精品| 中文字幕免费在线视频6| 七月丁香在线播放| 国产精品人妻久久久久久| 综合色丁香网| 国产精品.久久久| 亚洲av熟女| 三级男女做爰猛烈吃奶摸视频| 禁无遮挡网站| 22中文网久久字幕| 女人被狂操c到高潮| 热99在线观看视频| 国产又黄又爽又无遮挡在线| 天堂√8在线中文| 久久久欧美国产精品| 久久精品久久精品一区二区三区| 成年版毛片免费区| 亚洲综合精品二区| 人体艺术视频欧美日本| 在现免费观看毛片| 精品熟女少妇av免费看| 毛片女人毛片| 国产高潮美女av| 人妻系列 视频| 97人妻精品一区二区三区麻豆| 99热这里只有精品一区| 99热网站在线观看| 蜜桃久久精品国产亚洲av| 高清日韩中文字幕在线| 乱码一卡2卡4卡精品| 国产精品三级大全| 在线播放国产精品三级| 亚州av有码| 波多野结衣巨乳人妻| 91在线精品国自产拍蜜月| 人体艺术视频欧美日本| 久99久视频精品免费| 国产高潮美女av| 久久精品久久久久久久性| 男人舔女人下体高潮全视频| videos熟女内射| 国产精品三级大全| 久久鲁丝午夜福利片| 春色校园在线视频观看| 亚洲精品乱码久久久久久按摩| 欧美成人精品欧美一级黄| 高清日韩中文字幕在线| 久久精品熟女亚洲av麻豆精品 | 2021天堂中文幕一二区在线观| 人妻少妇偷人精品九色| 中文字幕av成人在线电影| 99热全是精品| 一个人看的www免费观看视频| 免费播放大片免费观看视频在线观看 | 91在线精品国自产拍蜜月| 精品国产露脸久久av麻豆 | 天堂√8在线中文| 97超碰精品成人国产| 热99在线观看视频| 啦啦啦观看免费观看视频高清| 99九九线精品视频在线观看视频| 国产精品.久久久| 久久久国产成人免费| 91aial.com中文字幕在线观看| 成人亚洲精品av一区二区| 精品国产三级普通话版| 国产色爽女视频免费观看| 欧美一区二区国产精品久久精品| 欧美高清成人免费视频www| 天天躁日日操中文字幕| 1000部很黄的大片| 免费看av在线观看网站| 亚洲av男天堂| 看黄色毛片网站| 一个人观看的视频www高清免费观看| 在线观看66精品国产| 精品人妻熟女av久视频| 少妇被粗大猛烈的视频| 久久久精品大字幕| 色综合站精品国产| 国产白丝娇喘喷水9色精品| 简卡轻食公司| 3wmmmm亚洲av在线观看| 精品久久久久久久久久久久久| 变态另类丝袜制服| 91狼人影院| 最后的刺客免费高清国语| 中文字幕av成人在线电影| 女人被狂操c到高潮| 亚洲精品色激情综合| 伦精品一区二区三区| 美女cb高潮喷水在线观看| 只有这里有精品99| 久久草成人影院| 久久久久精品久久久久真实原创|