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

    基于Hilbert包絡(luò)譜和SVM的齒輪故障診斷

    2018-01-09 09:12:36崔天曉周小龍劉文浩王增鵬鄭佳超
    東北電力大學(xué)學(xué)報 2017年6期
    關(guān)鍵詞:故障診斷故障信號

    崔天曉,周小龍,劉文浩,王增鵬,馬 靜,鄭佳超

    (1.東北電力大學(xué) 工程訓(xùn)練教學(xué)中心,吉林 吉林132012;2.吉林石化工程設(shè)計有限公司,吉林 吉林 132012;3.中國石油吉林石化公司,吉林 吉林 132012;4.敦化市亞聯(lián)機械制造有限公司,吉林 敦化 133700)

    基于Hilbert包絡(luò)譜和SVM的齒輪故障診斷

    崔天曉1,周小龍1,劉文浩2,王增鵬3,馬 靜3,鄭佳超4

    (1.東北電力大學(xué) 工程訓(xùn)練教學(xué)中心,吉林 吉林132012;2.吉林石化工程設(shè)計有限公司,吉林 吉林 132012;3.中國石油吉林石化公司,吉林 吉林 132012;4.敦化市亞聯(lián)機械制造有限公司,吉林 敦化 133700)

    針對齒輪振動信號非平穩(wěn)性和故障特征難以提取的問題,提出一種基于希爾伯特(Hilbert)包絡(luò)譜和支持向量機(SVM)相結(jié)合的故障診斷方法。首先,采用基于鏡像延拓的改進(jìn)經(jīng)驗?zāi)B(tài)分解(EMD)方法,將齒輪振動信號分解成一系列含有信號特征的固有模態(tài)函數(shù)(IMF);其次,根據(jù)正交性原理,選取包含信號主要信息的模態(tài)分量進(jìn)行希爾伯特變換并求出包絡(luò)譜;最后,將包絡(luò)譜所求出的特征幅值比作為支持向量機分類器的輸入來識別齒輪的工作狀態(tài)。試驗結(jié)果表明,該方法能有效地識別和診斷出齒輪的工作狀態(tài),可應(yīng)用于該類問題的故障診斷研究。

    經(jīng)驗?zāi)B(tài)分解;鏡像延拓;包絡(luò)譜;支持向量機;齒輪;故障診斷

    當(dāng)機械系統(tǒng)中的齒輪發(fā)生故障時,常以其振動信號作為故障診斷對象。然而,由于環(huán)境噪聲等外界因素和齒輪系統(tǒng)的特性,齒輪振動信號往往表現(xiàn)出非線性和非平穩(wěn)性的特點[3],這在很大程度上增加了故障診斷的難度。而常用的故障診斷方法,如:傅里葉變換和小波變換等,在分析該類信號時都具有自身的局限性[4],難以取得良好的診斷結(jié)果。由Huang[5]等人提出的EMD方法由于具有良好的自適應(yīng)性,因此在處理此類信號的故障診斷方面得到廣泛應(yīng)用[6-8]。但EMD自身的端點效應(yīng)問題限制了HHT方法在機械故障特征提取及診斷領(lǐng)域的后續(xù)發(fā)展。

    鑒于上述分析,本文提出一種基于改進(jìn)EMD的Hilbert包絡(luò)譜和SVM相結(jié)合的齒輪故障診斷方法。該方法首先利用鏡像延拓抑制EMD分解過程中所產(chǎn)生的端點效應(yīng)問題,根據(jù)各IMF分量間應(yīng)完全正交的原理,選取包含主要故障信息的IMF并求出相應(yīng)的包絡(luò)譜;最后,以不同故障所對應(yīng)包絡(luò)譜的特征幅值比作為SVM的輸入進(jìn)行齒輪的故障診斷。試驗表明,該方法能有效應(yīng)用于齒輪故障診斷。

    1 改進(jìn)的EMD

    1.1 基于鏡像延拓的EMD

    由于EMD算法忽略了對信號端點處的約束問題,導(dǎo)致在EMD分解信號的過程中,采用三次樣條函數(shù)對信號極值點進(jìn)行插值時,會出現(xiàn)擬合誤差,隨著分解的不斷進(jìn)行,誤差量逐漸增大,最終導(dǎo)致各IMF分量在端點處出現(xiàn)飛翼現(xiàn)象,失去原有物理意義[9]。

    在此,采用鏡像延拓法解決EMD的端點效應(yīng)問題。改進(jìn)EMD能夠根據(jù)信號的特性自適應(yīng)地將其分解成一系列IMF,并降低端點效應(yīng)對于EMD分解準(zhǔn)確性的影響?;阽R像延拓EMD方法的具體算法如下:

    找出信號x(t)的所有局部極值點,將鏡面放置在信號兩側(cè)極值點的位置,并將數(shù)據(jù)序列即為K1、K2。當(dāng)鏡面置于極大值或極小值位置時,分別有:K1=Km(1)、K2=Km(2)以及K1=Kn(1)、K2=Kn(2)。此時,延拓后信號序列為(t,x):

    (1)

    (2)

    經(jīng)鏡像延拓后,原數(shù)據(jù)序列和延拓后數(shù)據(jù)序列形成一封閉的環(huán)狀曲線,數(shù)據(jù)的上、下包絡(luò)線完全由內(nèi)部數(shù)據(jù)確定,從根本上避免了端點效應(yīng)問題,最終輸出的時候,只輸出兩鏡面內(nèi)的數(shù)據(jù)。

    鏡像延拓方法從端點效應(yīng)產(chǎn)生的原因入手,使信號序列的包絡(luò)線完全由端點的數(shù)據(jù)確定。同時,對于環(huán)形數(shù)據(jù)列而言,僅一次延拓就可較好地解決端點效應(yīng)問題,有效提高了EMD的分解效率。

    2.2 基于正交性的敏感IMF選擇算法

    在理想狀態(tài)下,機械振動信號所包含的各個IMF分量是完全正交的[10],但由于端點效應(yīng)的影響,各IMF分量間不可能完全正交。因此,包含信號主要信息的各IMF分量應(yīng)與原信號有較高的相關(guān)系數(shù),可以此作為判別敏感IMF分量的依據(jù)。

    信號x(t)和其各IMF分量ci(t)之間的相關(guān)系數(shù)為

    (3)

    通過計算各階IMF分量同原信號間的相關(guān)系數(shù),并由所分析信號的特點選取門限閾值,從而確定對故障信息敏感的IMF分量。

    2 Hilbert包絡(luò)譜及其特征幅值比

    2.1 Hilbert包絡(luò)譜

    由于改進(jìn)EMD方法所獲得的每個IMF分量都是一個單分量的調(diào)制信號,因此可以對每個IMF分量進(jìn)行包絡(luò)譜分析[11]。

    設(shè)ci(t)是信號x(t)經(jīng)EMD分解得到的任意1個IMF分量,對其作Hilbert變換:

    (4)

    進(jìn)一步可以求出包絡(luò)信號:

    (5)

    對求解此包絡(luò)信號的相關(guān)譜便可獲得其包絡(luò)譜。

    2.2 Hilbert包絡(luò)譜的特征幅值比

    對于信號的包絡(luò)譜而言,能夠有效反映旋轉(zhuǎn)機械的周期沖擊情況,若對此包絡(luò)信號選取一合理的評價指標(biāo),則可有效反映旋轉(zhuǎn)機械的自身特性。當(dāng)齒輪箱中的齒輪出現(xiàn)故障時,在包絡(luò)譜中故障特征頻率處會出現(xiàn)相應(yīng)的譜峰。因此,可將不同齒輪故障在其特征頻率處的幅值與齒輪箱中的齒輪(主動齒輪或從動齒輪)在包絡(luò)譜中的故障特征頻率處的幅值比定義為判別齒輪故障的特征幅值比,即

    (6)

    式中:A(frk)為齒輪箱不同工況下的特征頻率處的幅值,其中k=1,2,3,分別代表齒輪正常狀態(tài)、主從齒輪故障和從動齒輪故障;A(fi) 為齒輪故障狀態(tài)下的特征頻率處的幅值,本文中A(fi) 為從動齒輪故障的特征頻率處的幅值;γ為特征幅值比。

    3 試驗驗證

    在QPZZ-II型齒輪故障模擬平臺上進(jìn)行了齒輪箱正常運行、主動齒輪磨損和從動齒輪斷齒等三種不同工況的試驗。實驗所用齒輪箱為一級減速器,該齒輪箱由輸入軸通過聯(lián)軸器與驅(qū)動電機相連。齒輪箱中驅(qū)動齒輪和從動齒輪的齒數(shù)分別為55和75,且都為標(biāo)準(zhǔn)直齒輪,輸入軸和輸出軸的平均轉(zhuǎn)速及轉(zhuǎn)頻分別為850 r/min、 523 r/min以及14.14 Hz和10.39 Hz,嚙合頻率為779 Hz。試驗過程中,加速度傳感器安裝在靠近從動齒輪負(fù)載一側(cè)的軸承座上,設(shè)置采樣頻率為5 120 Hz,采樣時間為1 s。圖1為不同狀態(tài)下齒輪的振動信號。

    圖1 不同狀態(tài)下齒輪時域圖

    由圖1可知,不同工況下,齒輪振動信號在時域中均表現(xiàn)為脈沖沖擊,同高斯幅值調(diào)制信號相類似。雖然不同故障的時域信號存在細(xì)微差異,但是以此為依據(jù)不能準(zhǔn)確區(qū)分故障類型。

    為驗證所提改進(jìn)EMD方法的有效性,分別采用兩種EMD方法對齒輪箱正常信號進(jìn)行分析,得到的IMF分量階數(shù)分別為10和14。改進(jìn)EMD方法由于在一定程度上抑制了端點效應(yīng)問題的產(chǎn)生,因此在其分解過程中包絡(luò)擬合誤差量有所減少,有效減少了虛假IMF分量。

    為有效對比,圖2為兩種方法分解得到的前8階IMF分量。對比圖2(a)和圖2(b)可知,原始EMD方法分解出的IMF分量靠近信號兩端點處明顯存在異動,在0 s~0.02 s和0.97 s~1 s兩個時間區(qū)域內(nèi)產(chǎn)生畸變現(xiàn)象,其中低頻IMF分量尤為明顯。而改進(jìn)EMD方法分解出的各IMF分量的端點處沒有出現(xiàn)飛翼,同時畸變現(xiàn)象得到較好抑制,所分解出的各階IMF分量的波形和實際情況更加相符。由此證明改進(jìn)EMD方法能夠更加有效地提取所分析信號的特征信息,保證后續(xù)分析的準(zhǔn)確性。

    圖2 不同EMD方法分解對比

    分別對齒輪磨損故障和斷齒故障信號進(jìn)行改進(jìn)EMD分解,都分解出10階IMF分量,其中前5階IMF分量與其原信號間的相關(guān)系數(shù),如表1所示。

    表1 不同狀態(tài)下各IMF分量與自身信號的相關(guān)系數(shù)

    由表1可知,雖然所分析的齒輪振動信號處于不同的工作狀態(tài),但各種狀態(tài)信號的前3階IMF分量集中了其主要故障信息,因此,本文將各狀態(tài)下齒輪信號的前3階IMF分量作為敏感IMF分量用于故障特征的提取。不同狀態(tài)下齒輪振動信號的高頻IMF1的包絡(luò)譜,如圖3所示。

    圖3 不同狀態(tài)下齒輪高頻IMF1包絡(luò)譜

    由圖3可以看出,齒輪磨損和斷齒故障信號的包絡(luò)譜中,其譜峰分別出現(xiàn)在頻率為43.13 Hz、85.63 Hz處,以及20.63 Hz、41.88 Hz處,而這正是驅(qū)動齒輪的3倍、6倍轉(zhuǎn)頻和從動齒輪的2倍、4倍轉(zhuǎn)頻,由此可以知,改進(jìn)EMD所獲得的包絡(luò)譜可以較好地反映出齒輪的故障轉(zhuǎn)頻,表2為不同狀態(tài)下的特征幅值比。

    表2 不同狀態(tài)下的特征幅值比

    由表2知,當(dāng)齒輪處于正常狀態(tài)齒輪傳動平穩(wěn),受到的沖擊較小,當(dāng)齒輪出現(xiàn)故障時,通常表現(xiàn)為高頻沖擊響應(yīng)特性,且在時域上不連續(xù),呈現(xiàn)沖擊與調(diào)制特性,因此,相對于磨損和斷齒狀態(tài),其特征幅值比較小。由于齒輪出現(xiàn)磨損故障,齒輪嚙合過程中高頻沖擊成分增多,所以特征幅值比較正常情況明顯增大。

    4 基于SVM的齒輪故障診斷

    SVM分類器在小樣本故障診斷中的分類性能優(yōu)于神經(jīng)網(wǎng)絡(luò)分類器[12],而目前工程實際應(yīng)用時所采用的都為小樣本的情況,因此,該方法在機械故障診斷領(lǐng)域應(yīng)用較廣[13]。

    基于Hilbert包絡(luò)譜和SVM非齒輪故障診斷方法使用基于鏡像延拓的EMD對不同工況下的齒輪振動信號進(jìn)行預(yù)處理,求出不同工況信號的前3階IMF分量的包絡(luò)譜,并將其特征幅值比γ作為SVM輸入來識別齒輪的工作狀態(tài)和故障類型。具體診斷步驟如下:

    (1)分別采集不同工況下的齒輪振動信號30組,其中采用頻率為5 120 Hz,其中任選10組作為SVM的樣本,剩余作為測試數(shù)據(jù);

    (2)采用改進(jìn)EMD方法對不同工況的齒輪振動信號進(jìn)行分析,選擇各信號的前3階IMF分量作為分析對象;

    (3)對敏感IMF分量進(jìn)行Hilbert變換,求出其包絡(luò)譜,并根據(jù)公式(6)計算每組信號的特征幅值比γ;

    (4)采用2個SVM組成故障診斷分類器。并按照步驟(2)~步驟(3)求出信號的特征向量γ,分類結(jié)果以SVM分類器的輸出為判別結(jié)果。

    將SVM1定義決策函數(shù)y=+1定義為正常狀態(tài),y=-1定義為磨損或斷齒故障。將SVM2定義決策函數(shù)y=+1為磨損故障,y=-1定義斷齒故障。

    將之前所確定的訓(xùn)練樣本數(shù)據(jù)和測試樣本數(shù)據(jù)輸入SVM模型,其結(jié)果為:磨損故障全部識別正確,而正常和斷齒故障狀態(tài)下各有1組數(shù)據(jù)識別錯誤,從而確定總體識別正確率為96.67%。引起誤判的原因可能是測量過程中的儀器誤差或操作誤差所引起的,但總體識別效果較為理想。部分識別結(jié)果如表3所示。

    表3 部分測試樣本特征向量和SVM識別結(jié)果

    6 結(jié) 論

    (1)采用敏感IMF分量所得包絡(luò)譜的特征幅值比作為特征向量,可以有效增加故障特征信息提取的準(zhǔn)確性;

    (2)基于Hilbert包絡(luò)譜和SVM相結(jié)合的方法可以有效辨識齒輪的工作狀態(tài),為齒輪的故障診斷研究提供了一種切實可行的方法。

    [1] 褚福磊,彭志科,馮志鵬,等.機械故障診斷中的現(xiàn)代信號處理方法[M].北京:科學(xué)技術(shù)出版社,2009:165.

    [2] 于德介,程軍圣,楊宇.Hilbert-Huang 變換在齒輪故障診斷中的應(yīng)用[J].機械工程學(xué)報,2005,41(6):102-107.

    [3] 謝孟龍.變工況齒輪箱振動信號處理與故障診斷方法研究[D].石家莊:河北科技大學(xué),2013.

    [4] 周小龍.希爾伯特-黃變換在故障診斷中的應(yīng)用[D].長春:長春工業(yè)大學(xué),2014.

    [5] N.E.Huang,Z.Shen,R.L.Steven,et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis[J].Proc.R.Soc.Lond A,1998,454(1971):903-995.

    [6] W.X.Yang.Interpretation of mechanical signals using an improved Hilbert-Huang transform[J].Mechanical Systems and Signal Processing,2008,22(5):1061-1071.

    [7] B.Xu,S.Yuan,M.Wang,et al.Determining impact induced damage by lamb wave mode extracted by EMD method[J].Measurement,2015,65:120-128.

    [8] X.Zhao,T.H.Patel,M.J.Zuo.Multivariate EMD and full spectrum based condition monitoring for rotating machinery[J].Mechanical Systems and Signal Processing,2012,27(1):712-728.

    [9] 王傳菲,安鋼,王凱,等.基于鏡像延拓和神經(jīng)網(wǎng)絡(luò)的EMD端點效應(yīng)改進(jìn)方法[J].裝甲兵工程學(xué)院學(xué)報,2010,24(2):62-65.

    [10] 方琨,王渝,馬利兵,等.基于RO-SBM的Hilbert-Huang變換端點效應(yīng)抑制方法[J].振動、測試與診斷,2013,33(2):319-324.

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

    [12] 于德介,程軍圣,楊宇.機械故障診斷的Hilbert-Huang變換方法[M].北京:科學(xué)技術(shù)出版社,2006:155-159.

    [13] 周小龍,楊恭勇,梁秀霞,等.基于EMD重構(gòu)和SVM的滾動軸承故障診斷方法研究[J].東北電力大學(xué)學(xué)報,2016,36(6):71-75.

    GearFaultDiagnosisBasedonHilbertEnvelopeSpectrumandSVM

    CuiTianxiao1,ZhouXiaolong1,LiuWenhao2,WangZengpeng3,MaJing3,ZhengJiachao4

    (1.The Engineering Training Teaching Center,Northeast Electric Power University,Jilin Jilin 132012;2.Jilin Petrochemical Engineering Design Co.,ltd.,Jilin Jilin 132012;3.China Petroleum Jilin Petrochemical Company,Jilin Jilin 132012;4.Yalian Machinery Manufacturing Co.,ltd.,Dunhua Jilin 133700)

    Aiming at the non-stationary characteristic of the gear vibration signal and the difficulty to get the fault features in its fault diagnosis,a gear fault diagnosis method based on Hilbert envelope spectrum and support vector machine (SVM) is proposed.Firstly,the fault signal is decomposed by empirical mode decomposition (EMD) and a set of intrinsic mode function (IMF) is obtained,which contain the feature of the gear.Then,the envelope spectrums of the sensitive IMFs that selected by the orthogonality principle of each IMFs are obtained.Finally,the characteristic amplitude ratio of the envelope spectrums as input vectors of support vector machine,and identify the gear fault pattern and condition.The experiment shows that this method can identify gear fault patterns effectively and it has a practical value.

    Empirical mode decomposition;Mirror extension;Envelope spectrum;Support vector machine;Gear;Fault diagnosis

    2017-01-06

    崔天曉(1989-),男,碩士,助理實驗師,主要研究方向:金屬材料加工.

    電子郵箱:2563202296@qq.com(崔天曉);196389679@qq.com(周小龍);liu7wenhao@163.com(劉文浩);362839032@qq.com(王增鵬);164475236@qq.com(馬靜);172424852@qq.com(鄭佳超)

    1005-2992(2017)06-0056-06

    TH17;TH133.33

    A

    猜你喜歡
    故障診斷故障信號
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    故障一點通
    基于FPGA的多功能信號發(fā)生器的設(shè)計
    電子制作(2018年11期)2018-08-04 03:25:42
    奔馳R320車ABS、ESP故障燈異常點亮
    基于LabVIEW的力加載信號采集與PID控制
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    故障一點通
    江淮車故障3例
    基于LCD和排列熵的滾動軸承故障診斷
    精品国产国语对白av| 一个人免费看片子| 欧美 亚洲 国产 日韩一| 久久久国产一区二区| 欧美亚洲 丝袜 人妻 在线| 久久天堂一区二区三区四区| 777米奇影视久久| 午夜福利乱码中文字幕| 一二三四中文在线观看免费高清| 久久精品亚洲av国产电影网| 日韩欧美一区视频在线观看| 久久99热这里只频精品6学生| 99精品久久久久人妻精品| 成年美女黄网站色视频大全免费| 欧美黄色片欧美黄色片| 欧美成人精品欧美一级黄| 如何舔出高潮| 免费日韩欧美在线观看| 极品人妻少妇av视频| 久久久国产一区二区| 午夜福利免费观看在线| 飞空精品影院首页| 国产片特级美女逼逼视频| 国产精品麻豆人妻色哟哟久久| 美女国产高潮福利片在线看| 99九九在线精品视频| 岛国毛片在线播放| 五月开心婷婷网| 日本猛色少妇xxxxx猛交久久| 亚洲,欧美,日韩| 热99国产精品久久久久久7| 无限看片的www在线观看| 欧美亚洲 丝袜 人妻 在线| 永久免费av网站大全| 好男人视频免费观看在线| 国产精品一二三区在线看| 人人妻人人澡人人爽人人夜夜| 成年人免费黄色播放视频| 国产又爽黄色视频| 一级毛片我不卡| 大片电影免费在线观看免费| 日本一区二区免费在线视频| 国产成人精品无人区| 亚洲色图综合在线观看| 人人澡人人妻人| 国产成人精品久久二区二区91 | 国产女主播在线喷水免费视频网站| 一本一本久久a久久精品综合妖精| 最近的中文字幕免费完整| 精品少妇一区二区三区视频日本电影 | 国产成人免费无遮挡视频| 久久精品久久久久久久性| 性色av一级| 久久精品熟女亚洲av麻豆精品| 男女无遮挡免费网站观看| 精品免费久久久久久久清纯 | 色吧在线观看| 丝袜美足系列| av线在线观看网站| 国产欧美日韩一区二区三区在线| 日韩人妻精品一区2区三区| 国产又色又爽无遮挡免| 亚洲国产欧美日韩在线播放| 一区二区日韩欧美中文字幕| 国产成人精品在线电影| 青青草视频在线视频观看| 天天影视国产精品| 日韩伦理黄色片| 最新在线观看一区二区三区 | 青草久久国产| 99热网站在线观看| 亚洲欧美一区二区三区黑人| 天堂俺去俺来也www色官网| 日韩av免费高清视频| 欧美日韩一区二区视频在线观看视频在线| 精品少妇黑人巨大在线播放| 亚洲国产精品999| 美国免费a级毛片| 午夜免费观看性视频| 青春草视频在线免费观看| 如何舔出高潮| www.自偷自拍.com| 一边摸一边做爽爽视频免费| 亚洲欧美日韩另类电影网站| 老司机在亚洲福利影院| 久久久久久免费高清国产稀缺| 中文字幕亚洲精品专区| 中文字幕亚洲精品专区| 免费日韩欧美在线观看| 中国国产av一级| 国产成人精品久久久久久| 丝袜人妻中文字幕| 国产亚洲一区二区精品| 一级a爱视频在线免费观看| 久久影院123| 精品一区二区免费观看| 亚洲av男天堂| 亚洲一级一片aⅴ在线观看| 国产亚洲最大av| 精品一区二区三区av网在线观看 | 欧美人与性动交α欧美精品济南到| 日本av免费视频播放| 亚洲美女黄色视频免费看| av网站免费在线观看视频| 免费女性裸体啪啪无遮挡网站| 成人午夜精彩视频在线观看| 久久天堂一区二区三区四区| 成年人免费黄色播放视频| 日韩电影二区| 国产在视频线精品| 亚洲国产精品成人久久小说| 9191精品国产免费久久| 色婷婷久久久亚洲欧美| 一本—道久久a久久精品蜜桃钙片| 久久久国产精品麻豆| 欧美 日韩 精品 国产| 国产亚洲午夜精品一区二区久久| 99热网站在线观看| 日韩伦理黄色片| 日韩一卡2卡3卡4卡2021年| 日本猛色少妇xxxxx猛交久久| 丰满迷人的少妇在线观看| 欧美日韩视频精品一区| 人人澡人人妻人| 久久性视频一级片| 日韩不卡一区二区三区视频在线| 成人黄色视频免费在线看| 在线观看免费午夜福利视频| 国产精品无大码| 午夜91福利影院| 亚洲精品自拍成人| 成人毛片60女人毛片免费| 男女边摸边吃奶| 麻豆av在线久日| 亚洲精品美女久久久久99蜜臀 | 亚洲欧美日韩另类电影网站| 2021少妇久久久久久久久久久| 少妇人妻精品综合一区二区| 伊人亚洲综合成人网| 国产免费福利视频在线观看| 国产精品国产三级专区第一集| 岛国毛片在线播放| 久久ye,这里只有精品| 天堂中文最新版在线下载| 亚洲欧洲精品一区二区精品久久久 | 国产国语露脸激情在线看| 色94色欧美一区二区| 亚洲精品国产区一区二| 母亲3免费完整高清在线观看| 亚洲熟女毛片儿| 欧美日韩福利视频一区二区| 国产在线免费精品| 国产福利在线免费观看视频| 狠狠婷婷综合久久久久久88av| 午夜福利在线免费观看网站| 日本黄色日本黄色录像| 亚洲欧美成人精品一区二区| 欧美激情 高清一区二区三区| 久久久国产精品麻豆| 制服丝袜香蕉在线| 欧美精品一区二区大全| 纵有疾风起免费观看全集完整版| 日本91视频免费播放| 亚洲欧美激情在线| 日韩av免费高清视频| 可以免费在线观看a视频的电影网站 | 国产欧美亚洲国产| 高清在线视频一区二区三区| 亚洲男人天堂网一区| 欧美精品亚洲一区二区| 亚洲欧洲精品一区二区精品久久久 | 日韩免费高清中文字幕av| 亚洲一卡2卡3卡4卡5卡精品中文| 最近中文字幕2019免费版| 欧美日韩亚洲综合一区二区三区_| 大香蕉久久网| h视频一区二区三区| 久久99精品国语久久久| 久久免费观看电影| 亚洲成人国产一区在线观看 | 午夜日本视频在线| 久久久欧美国产精品| 久久久欧美国产精品| 欧美老熟妇乱子伦牲交| 国产精品.久久久| 九九爱精品视频在线观看| 亚洲自偷自拍图片 自拍| 亚洲欧洲国产日韩| 国产免费视频播放在线视频| 你懂的网址亚洲精品在线观看| 久久精品人人爽人人爽视色| 久久狼人影院| 晚上一个人看的免费电影| 不卡视频在线观看欧美| 国产av一区二区精品久久| 五月天丁香电影| 免费不卡黄色视频| 国产在线视频一区二区| 欧美av亚洲av综合av国产av | 婷婷色综合大香蕉| 免费观看人在逋| 中文字幕人妻丝袜制服| 少妇 在线观看| 久久精品久久精品一区二区三区| 欧美最新免费一区二区三区| 大陆偷拍与自拍| 一区福利在线观看| 老司机靠b影院| 成年人午夜在线观看视频| av福利片在线| 日本av免费视频播放| 欧美av亚洲av综合av国产av | 日日撸夜夜添| 国产精品偷伦视频观看了| 国产亚洲午夜精品一区二区久久| 国产一区二区三区综合在线观看| 少妇猛男粗大的猛烈进出视频| 如何舔出高潮| 王馨瑶露胸无遮挡在线观看| 大片电影免费在线观看免费| 欧美日韩成人在线一区二区| 人人妻,人人澡人人爽秒播 | 中文字幕人妻熟女乱码| 欧美精品人与动牲交sv欧美| 亚洲精品成人av观看孕妇| 国产乱来视频区| 亚洲精品久久午夜乱码| 中文字幕人妻熟女乱码| 在线观看三级黄色| √禁漫天堂资源中文www| av卡一久久| 午夜91福利影院| 妹子高潮喷水视频| 国产成人系列免费观看| 久久天堂一区二区三区四区| 人人妻人人爽人人添夜夜欢视频| 精品少妇黑人巨大在线播放| 免费女性裸体啪啪无遮挡网站| 国产黄频视频在线观看| 电影成人av| 欧美老熟妇乱子伦牲交| 黄色视频不卡| 国产在线免费精品| 卡戴珊不雅视频在线播放| 搡老乐熟女国产| 久久人妻熟女aⅴ| 99久久人妻综合| 免费黄频网站在线观看国产| 欧美精品高潮呻吟av久久| 涩涩av久久男人的天堂| 中文字幕亚洲精品专区| 丰满少妇做爰视频| 男人添女人高潮全过程视频| 男的添女的下面高潮视频| 国产成人一区二区在线| 一区福利在线观看| 国产精品.久久久| 叶爱在线成人免费视频播放| 美女脱内裤让男人舔精品视频| 美女高潮到喷水免费观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲久久久国产精品| av.在线天堂| 国产精品av久久久久免费| 精品一品国产午夜福利视频| 精品视频人人做人人爽| 青春草亚洲视频在线观看| 国产99久久九九免费精品| 色精品久久人妻99蜜桃| 中文字幕人妻熟女乱码| 国产女主播在线喷水免费视频网站| 久久久久精品久久久久真实原创| 午夜福利网站1000一区二区三区| 国产黄频视频在线观看| 免费不卡黄色视频| 69精品国产乱码久久久| 国产精品人妻久久久影院| 精品一区二区免费观看| 国产女主播在线喷水免费视频网站| 成年av动漫网址| 亚洲国产毛片av蜜桃av| 欧美日韩亚洲高清精品| 老熟女久久久| 80岁老熟妇乱子伦牲交| 国产成人精品福利久久| 国产女主播在线喷水免费视频网站| 亚洲精品,欧美精品| 亚洲美女视频黄频| 中文精品一卡2卡3卡4更新| 国产精品三级大全| 亚洲精品久久午夜乱码| av网站免费在线观看视频| 欧美日韩视频精品一区| 天天操日日干夜夜撸| 少妇人妻 视频| 国产高清不卡午夜福利| 亚洲久久久国产精品| 国产精品国产三级国产专区5o| 男女午夜视频在线观看| 久久天堂一区二区三区四区| 国产成人精品在线电影| 大香蕉久久成人网| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 制服人妻中文乱码| 国产人伦9x9x在线观看| 亚洲综合色网址| 亚洲精品av麻豆狂野| a 毛片基地| 蜜桃国产av成人99| 人人妻人人爽人人添夜夜欢视频| 国产伦人伦偷精品视频| 日韩av不卡免费在线播放| 日韩不卡一区二区三区视频在线| 美女脱内裤让男人舔精品视频| 老鸭窝网址在线观看| 国产黄频视频在线观看| 两性夫妻黄色片| 国产精品 国内视频| 亚洲精品视频女| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久精品精品| 亚洲成人手机| 女人被躁到高潮嗷嗷叫费观| 免费观看人在逋| 久久久国产一区二区| 免费看av在线观看网站| 亚洲成av片中文字幕在线观看| 久久这里只有精品19| 午夜久久久在线观看| 久久久久国产精品人妻一区二区| 免费女性裸体啪啪无遮挡网站| 大香蕉久久成人网| 国产精品熟女久久久久浪| 亚洲成色77777| 成人手机av| 国产精品嫩草影院av在线观看| 最近的中文字幕免费完整| 中文字幕人妻丝袜制服| 成人18禁高潮啪啪吃奶动态图| 欧美成人午夜精品| 国产高清不卡午夜福利| 啦啦啦在线观看免费高清www| 青春草国产在线视频| 老司机影院毛片| 国产男女内射视频| 热99国产精品久久久久久7| 国产精品香港三级国产av潘金莲 | 亚洲精品中文字幕在线视频| 99re6热这里在线精品视频| 女人精品久久久久毛片| 国产精品久久久久久久久免| 岛国毛片在线播放| www.熟女人妻精品国产| 亚洲精品国产av蜜桃| 国产熟女欧美一区二区| 国产成人av激情在线播放| 2021少妇久久久久久久久久久| 久久久久久久久久久免费av| 亚洲国产精品一区三区| 免费在线观看视频国产中文字幕亚洲 | 亚洲情色 制服丝袜| 在线观看人妻少妇| 岛国毛片在线播放| 亚洲av福利一区| 婷婷色综合www| 亚洲男人天堂网一区| 国产黄频视频在线观看| 黄色怎么调成土黄色| 老汉色∧v一级毛片| 免费女性裸体啪啪无遮挡网站| 精品国产一区二区久久| 曰老女人黄片| 亚洲欧美成人精品一区二区| 精品一区二区三区四区五区乱码 | 永久免费av网站大全| 欧美国产精品va在线观看不卡| 老司机在亚洲福利影院| 亚洲欧美成人精品一区二区| 久久人人爽av亚洲精品天堂| 99re6热这里在线精品视频| 天堂8中文在线网| 人人妻人人爽人人添夜夜欢视频| 看免费av毛片| 性高湖久久久久久久久免费观看| 国产成人a∨麻豆精品| 亚洲av日韩精品久久久久久密 | 中文精品一卡2卡3卡4更新| 老熟女久久久| 18禁裸乳无遮挡动漫免费视频| 男女午夜视频在线观看| 777米奇影视久久| 青春草视频在线免费观看| 中文字幕色久视频| 精品国产一区二区三区四区第35| 国产亚洲欧美精品永久| 午夜免费男女啪啪视频观看| 欧美日韩视频高清一区二区三区二| 国产成人一区二区在线| 日日摸夜夜添夜夜爱| 久久久久久久久久久免费av| 一区二区三区乱码不卡18| 国产av精品麻豆| 天天躁日日躁夜夜躁夜夜| 成年人免费黄色播放视频| 亚洲av福利一区| 日韩一区二区视频免费看| 欧美另类一区| 午夜免费男女啪啪视频观看| 制服人妻中文乱码| 国产一区二区三区综合在线观看| 卡戴珊不雅视频在线播放| 色94色欧美一区二区| 人人妻人人添人人爽欧美一区卜| 欧美日本中文国产一区发布| 99久久99久久久精品蜜桃| 亚洲精品美女久久久久99蜜臀 | 免费看av在线观看网站| 国产国语露脸激情在线看| 如何舔出高潮| 色婷婷av一区二区三区视频| 一个人免费看片子| 色综合欧美亚洲国产小说| 秋霞伦理黄片| 精品一区二区三区av网在线观看 | 2018国产大陆天天弄谢| 久久久精品免费免费高清| av天堂久久9| 欧美人与善性xxx| 成人影院久久| 久久av网站| 久久青草综合色| 亚洲男人天堂网一区| 久久精品久久精品一区二区三区| www.精华液| 精品一品国产午夜福利视频| 青春草亚洲视频在线观看| 精品少妇黑人巨大在线播放| 丝袜美足系列| 国产免费福利视频在线观看| 亚洲五月色婷婷综合| 亚洲美女黄色视频免费看| 国产精品久久久久久久久免| 精品一区二区免费观看| 精品一区二区三卡| 精品第一国产精品| 国产一区二区激情短视频 | 两个人看的免费小视频| 日本黄色日本黄色录像| 国产成人a∨麻豆精品| 亚洲成色77777| 午夜福利影视在线免费观看| 日韩av不卡免费在线播放| 在线观看www视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 国产在线一区二区三区精| 狂野欧美激情性xxxx| 午夜免费观看性视频| 国产一区二区 视频在线| 19禁男女啪啪无遮挡网站| 亚洲国产最新在线播放| 亚洲av日韩在线播放| 精品少妇一区二区三区视频日本电影 | e午夜精品久久久久久久| 搡老岳熟女国产| svipshipincom国产片| 一本一本久久a久久精品综合妖精| 日日摸夜夜添夜夜爱| 国产av码专区亚洲av| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 性色av一级| 久久这里只有精品19| 天天影视国产精品| 国产男人的电影天堂91| 欧美av亚洲av综合av国产av | 国产乱人偷精品视频| 菩萨蛮人人尽说江南好唐韦庄| 国产成人欧美| 午夜av观看不卡| 秋霞在线观看毛片| 老鸭窝网址在线观看| 亚洲精品久久久久久婷婷小说| av片东京热男人的天堂| 欧美成人精品欧美一级黄| av.在线天堂| 在线 av 中文字幕| 美女大奶头黄色视频| 精品卡一卡二卡四卡免费| 日本黄色日本黄色录像| 一区二区三区四区激情视频| 久久久欧美国产精品| 婷婷色综合大香蕉| 亚洲欧美激情在线| 亚洲欧美一区二区三区黑人| 色播在线永久视频| 嫩草影视91久久| 天天躁夜夜躁狠狠久久av| 亚洲自偷自拍图片 自拍| 国产成人一区二区在线| 黄色怎么调成土黄色| 久久精品国产亚洲av涩爱| 黑人猛操日本美女一级片| 巨乳人妻的诱惑在线观看| 精品一区二区免费观看| 日日啪夜夜爽| 日韩一卡2卡3卡4卡2021年| 高清在线视频一区二区三区| 亚洲欧美成人综合另类久久久| 99久久人妻综合| 久久这里只有精品19| 悠悠久久av| 精品少妇黑人巨大在线播放| 久久99热这里只频精品6学生| 在线精品无人区一区二区三| 午夜免费鲁丝| www.av在线官网国产| 国产精品女同一区二区软件| av在线app专区| 最近最新中文字幕免费大全7| 丝袜在线中文字幕| 老鸭窝网址在线观看| 最近2019中文字幕mv第一页| 免费女性裸体啪啪无遮挡网站| 最近手机中文字幕大全| 国产av精品麻豆| 午夜日韩欧美国产| 亚洲成色77777| 日韩欧美精品免费久久| 叶爱在线成人免费视频播放| 夜夜骑夜夜射夜夜干| 国产av码专区亚洲av| 国语对白做爰xxxⅹ性视频网站| 九草在线视频观看| 看免费av毛片| 深夜精品福利| 两个人看的免费小视频| 深夜精品福利| 精品少妇久久久久久888优播| 亚洲精品美女久久av网站| 国精品久久久久久国模美| 伦理电影大哥的女人| av在线播放精品| 亚洲av日韩在线播放| 国产一级毛片在线| 色综合欧美亚洲国产小说| 日韩,欧美,国产一区二区三区| 在现免费观看毛片| 久久精品亚洲av国产电影网| av在线老鸭窝| 久久综合国产亚洲精品| 精品国产乱码久久久久久男人| 建设人人有责人人尽责人人享有的| 色吧在线观看| 中文字幕人妻丝袜一区二区 | a级片在线免费高清观看视频| 午夜日韩欧美国产| 又大又爽又粗| 高清在线视频一区二区三区| 亚洲欧洲日产国产| 久久女婷五月综合色啪小说| 少妇 在线观看| 日韩一卡2卡3卡4卡2021年| 一区福利在线观看| 男男h啪啪无遮挡| 婷婷色麻豆天堂久久| 久久国产精品男人的天堂亚洲| 日本爱情动作片www.在线观看| 51午夜福利影视在线观看| 久久久久久免费高清国产稀缺| 狂野欧美激情性xxxx| 男女床上黄色一级片免费看| 国产精品嫩草影院av在线观看| 亚洲欧美精品自产自拍| 男人爽女人下面视频在线观看| 久久精品亚洲av国产电影网| 亚洲成人免费av在线播放| 亚洲国产欧美在线一区| 久久精品aⅴ一区二区三区四区| 中国国产av一级| 你懂的网址亚洲精品在线观看| 蜜桃国产av成人99| 美女扒开内裤让男人捅视频| 人人妻人人澡人人爽人人夜夜| 九色亚洲精品在线播放| 精品一区二区免费观看| 深夜精品福利| a级毛片黄视频| 亚洲国产av新网站| 亚洲第一青青草原| 美国免费a级毛片| 国产熟女午夜一区二区三区| 国产色婷婷99| 久久97久久精品| 男的添女的下面高潮视频| 99久久人妻综合| 9色porny在线观看| 永久免费av网站大全| 亚洲av成人不卡在线观看播放网 | 亚洲免费av在线视频| 亚洲成av片中文字幕在线观看| 欧美国产精品va在线观看不卡| 久久精品久久久久久噜噜老黄| 午夜免费男女啪啪视频观看| 最近手机中文字幕大全| 一区在线观看完整版| 午夜福利乱码中文字幕| 9191精品国产免费久久| 中文欧美无线码| 日韩精品免费视频一区二区三区| 丝袜美足系列| 国产一区二区三区av在线| 秋霞伦理黄片| 亚洲精品国产av成人精品| 别揉我奶头~嗯~啊~动态视频 |