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

    基于等距映射和最小二乘支持向量機(jī)的轉(zhuǎn)盤軸承故障識(shí)別方法

    2019-07-22 01:58:24趙祥龍陳捷洪榮晶潘裕斌
    軸承 2019年6期
    關(guān)鍵詞:降維分形軸承

    趙祥龍,陳捷,洪榮晶,潘裕斌

    (南京工業(yè)大學(xué) 機(jī)械與動(dòng)力工程學(xué)院,南京 211816)

    轉(zhuǎn)盤軸承作為旋轉(zhuǎn)機(jī)械的核心部件,工作環(huán)境復(fù)雜且惡劣,在實(shí)際使用過程中易產(chǎn)生故障并帶來巨大損失,因此,實(shí)現(xiàn)轉(zhuǎn)盤軸承早期故障的識(shí)別具有十分重要的意義[1-2],而分析軸承的振動(dòng)信號(hào)并提取相關(guān)特征信息是故障診斷研究的關(guān)鍵技術(shù)[3-5]。

    目前,基于時(shí)域、頻域以及時(shí)頻域的特征提取方法已經(jīng)得到廣泛應(yīng)用并取得了較好的效果[6-7]。然而,由于振動(dòng)信號(hào)具有很強(qiáng)的非線性、非平穩(wěn)性特征,極易被其他噪聲所掩蓋,為提高診斷效果,可結(jié)合數(shù)據(jù)約簡(jiǎn)技術(shù)實(shí)現(xiàn)對(duì)振動(dòng)信號(hào)特征信息的有效提取。依據(jù)振動(dòng)信號(hào)在一定尺度范圍內(nèi)包含分形特征的特性,文獻(xiàn)[8]利用振動(dòng)信號(hào)多分形特征實(shí)現(xiàn)了機(jī)械設(shè)備的狀態(tài)監(jiān)測(cè)和故障診斷;文獻(xiàn)[9]利用主成分分析方法提取時(shí)域及頻域特征中的敏感特征,實(shí)現(xiàn)降維的目的,解決了特征提取時(shí)維數(shù)過多而導(dǎo)致的信息冗余問題;文獻(xiàn)[10]則利用多維尺度分析方法對(duì)提取的軸承信號(hào)統(tǒng)計(jì)指標(biāo)進(jìn)行降維處理。

    上述維數(shù)約簡(jiǎn)方法針對(duì)線性問題時(shí)效果較好,對(duì)于非線性問題的使用效果并不明顯。而等距映射方法(Isometric Mapping,ISOMAP)是一種典型的流行學(xué)習(xí)算法,能夠很好的處理非線性問題,其通過建立源數(shù)據(jù)與降維數(shù)據(jù)之間的對(duì)等關(guān)系實(shí)現(xiàn)數(shù)據(jù)約簡(jiǎn)目的,且具有降維前后距離保持不變的特點(diǎn),能確保降維后的特征矩陣保留降維前的幾何信息。綜上所述,提出基于ISOMAP的降維方法,并利用最小二乘支持向量機(jī)(LSSVM)模型對(duì)獲得的低維特征進(jìn)行討論分析,以驗(yàn)證基于ISOMAP特征提取方法的識(shí)別效果。

    1 基于wavelet leader的多分形特征提取

    機(jī)械設(shè)備的振動(dòng)信號(hào)具有一定的分形特性,多分形分析方法不僅能從整體方面分析振動(dòng)信號(hào)的特性,而且能夠精細(xì)表征振動(dòng)信號(hào)的局部特性[11-12]。wavelet leader算法[13-15]具有堅(jiān)實(shí)的數(shù)學(xué)理論基礎(chǔ),計(jì)算簡(jiǎn)便,能夠很好的獲得振動(dòng)信號(hào)的多分形特征?;趙avelet leader方法得到的多分形譜D(q)和奇異指數(shù)h(q)的經(jīng)驗(yàn)公式為

    式中:wj為權(quán)重;q為多分量矩的階數(shù);L為上確界;U,V為統(tǒng)計(jì)量。

    多分形特征的D(q)-h(huán)(q)圖如圖1所示,可以看出其具有顯著的幾何結(jié)構(gòu)特征,可利用wavelet leader計(jì)算轉(zhuǎn)盤軸承振動(dòng)信號(hào)的多分形特征,并構(gòu)造出特征向量組 Ti={h(q)i1,D(q)i1,h(q)i6,D(q)i6,h(q)i11,D(q)i11},其中 i=1,2,…,n;n為樣本總數(shù)。由于多分形特征能夠反映振動(dòng)信號(hào)的幾何信息,因此,提取的多分形特征信息可用于表征轉(zhuǎn)盤軸承不同故障狀態(tài)下的信息。

    圖1 D(q)-h(huán)(q)圖Fig.1 Diagram of D(q)-h(huán)(q)

    2 基于ISOMAP與LSSVM分析的轉(zhuǎn)盤軸承故障類型識(shí)別

    2.1 等距映射流行學(xué)習(xí)方法

    流行學(xué)習(xí)是一種新的非線性降維算法,可用于機(jī)器學(xué)習(xí)及模式識(shí)別[16],在數(shù)據(jù)降維方面具有廣泛的應(yīng)用,其主要思想是將高維中的數(shù)據(jù)通過某種方式映射到低維空間中,并且依然能夠反映高維空間數(shù)據(jù)的結(jié)構(gòu)特征。

    高維數(shù)據(jù)空間中數(shù)據(jù)點(diǎn)之間的距離在低維空間中直接采用直線距離計(jì)算具有一定的誤差[16],因此,等距映射算法引進(jìn)鄰域圖的概念,采用鄰近點(diǎn)距離方法逼近測(cè)地線距離,實(shí)現(xiàn)以局部鄰域距離對(duì)樣本點(diǎn)間的真實(shí)測(cè)地線距離進(jìn)行評(píng)估的目的。源數(shù)據(jù)樣本點(diǎn)僅與其相鄰的樣本點(diǎn)相連,樣本點(diǎn)間的距離可直接通過歐式距離獲得,較遠(yuǎn)點(diǎn)間的距離通過最小路徑獲得,最終確立高維空間數(shù)據(jù)間距離與低維數(shù)據(jù)空間的對(duì)等關(guān)系,達(dá)到數(shù)據(jù)降維的效果。等距映射算法的詳細(xì)流程如下:

    1)構(gòu)造鄰域圖G。假設(shè)觀測(cè)空間數(shù)據(jù)集為X,xi為 X中的樣本點(diǎn),其中 i=1,2,…,n且 xi∈RS,S為高維空間維數(shù)。計(jì)算出任意2個(gè)樣本點(diǎn)xi與xj之間的歐式距離 dx(i,j),選取離各數(shù)據(jù)點(diǎn)距離最近的k個(gè)點(diǎn)或者選取以該樣本點(diǎn)為圓心,ε為半徑的圓內(nèi)所有點(diǎn)作為該樣本點(diǎn)的近鄰點(diǎn),記為以樣本點(diǎn)為節(jié)點(diǎn),歐式距離為邊即可構(gòu)造出樣本點(diǎn)間的鄰域關(guān)系圖G。

    2)計(jì)算測(cè)地距離矩陣D。利用Floyd算法或Dijkstra算法計(jì)算出最短路徑 dG(i,j),以逼近高維空間測(cè)地線距離矩陣 DG,DG={dG(i,j)}。

    3)利用MDS算法實(shí)現(xiàn)對(duì)高維數(shù)據(jù)集X的降維處理。假設(shè)數(shù)據(jù)集X中xi與xj之間的距離為距離矩陣DG中的元素dij,空間對(duì)應(yīng)表示為Z,Z={z1,z2,...,zt,...,zn},zt∈Rs,其中 s為低維空間的維數(shù),n為距離向量個(gè)數(shù)。同時(shí)任意2個(gè)樣本點(diǎn)在低維、高維空間的歐式距離相等,即‖zi-zj‖=dij。構(gòu)造矩陣B=ZTZ,其中B為降維后的內(nèi)積矩陣則有

    2.2 最小二乘支持向量機(jī)

    最小二乘支持向量機(jī)[17]是通過最小二乘法借助誤差平方和選擇超平面的方法構(gòu)造出平方損失函數(shù),將支持向量機(jī)的不等式約束問題轉(zhuǎn)化為等式約束,實(shí)現(xiàn)將二次規(guī)劃問題轉(zhuǎn)化為線性問題的方法。

    LSSVM建立的目標(biāo)函數(shù)和約束條件為

    式中:δi為松弛變量;λ為懲罰系數(shù);ω為權(quán)向量;b為偏差值。構(gòu)建Lagrange函數(shù),實(shí)現(xiàn)對(duì)最優(yōu)問題的求解

    對(duì)每個(gè)變量求偏導(dǎo)可得

    得到的線性方程組為

    式中:Y=[y1,y2,…,yn],α=[α1,α2,...,αn]T,Ω為 n×n的對(duì)稱矩陣,且 Ωij=φ(xi)Tφ(xj)=K(xi,xj);j=1,2,...,n。K(xi,xj)是滿足 Mercer條件的核函數(shù)。常用的核函數(shù)主要有徑向基核函數(shù)、多項(xiàng)式核函數(shù)、Sigmoid核函數(shù)。由于多項(xiàng)式核函數(shù)和Sigmoid核函數(shù)所需考慮的參數(shù)較多,相比單一參數(shù)的徑向基核函數(shù)存在參數(shù)設(shè)定工作量大,花費(fèi)時(shí)間較長(zhǎng)等問題。因此,引入徑向基函數(shù)作為核函數(shù),即 K(xi,xj)=exp(-γ‖xi-xj‖2),γ為核參數(shù)。最終得到的決策函數(shù)為

    由于LSSVM模型具有良好的數(shù)理統(tǒng)計(jì)性能,能夠有效識(shí)別出轉(zhuǎn)盤軸承不同的故障類型,因此將特征集的訓(xùn)練樣本對(duì)模型進(jìn)行訓(xùn)練,測(cè)試樣本檢測(cè)模型性能。轉(zhuǎn)盤軸承故障類型識(shí)別技術(shù)路線如圖2所示。

    圖2 轉(zhuǎn)盤軸承故障類型識(shí)別技術(shù)路線Fig.2 Technical route for recognition of fault type for slewing bearing

    3 實(shí)例分析

    3.1 轉(zhuǎn)盤軸承全壽命加速試驗(yàn)

    為驗(yàn)證所提方案的可行性與有效性,采用如圖3所示的試驗(yàn)臺(tái)進(jìn)行轉(zhuǎn)盤軸承全壽命加速試驗(yàn),試驗(yàn)臺(tái)主要由機(jī)械、液壓和測(cè)控3部分組成。由液壓馬達(dá)驅(qū)動(dòng)轉(zhuǎn)盤軸承旋轉(zhuǎn),3個(gè)液壓缸聯(lián)合作用實(shí)現(xiàn)轉(zhuǎn)盤軸承工作受力模擬;測(cè)控系統(tǒng)采用LabVIEW軟件平臺(tái)以及配套的硬件設(shè)備,通過OPC協(xié)議與S7-200控制系統(tǒng)完成通信,實(shí)現(xiàn)測(cè)控一體化的功能。

    圖3 轉(zhuǎn)盤軸承試驗(yàn)臺(tái)Fig.3 Test rig for slewing bearing

    試驗(yàn)選用QNA-730-22型內(nèi)齒式轉(zhuǎn)盤軸承(具體參數(shù)見表1),以4 r/min的轉(zhuǎn)速持續(xù)運(yùn)行12 d直至基本失效。拆解發(fā)現(xiàn),溝道及鋼球出現(xiàn)了嚴(yán)重?fù)p傷,如圖4所示。考慮到轉(zhuǎn)盤軸承工況的特殊性,為準(zhǔn)確獲取轉(zhuǎn)盤軸承不同部位的受力情況,采用4個(gè)加速度傳感器沿圓周方向間隔90°均勻分布在轉(zhuǎn)盤軸承的定圈端面,傳感器采樣頻率設(shè)置為2 048 Hz。試驗(yàn)結(jié)束后,對(duì)比分析發(fā)現(xiàn)4號(hào)加速度傳感器所采集信號(hào)的變化最為明顯(圖5),因此采用4號(hào)傳感器所采集信號(hào)進(jìn)行后續(xù)分析。

    表1 轉(zhuǎn)盤軸承的基本參數(shù)Tab.1 Basic parameters for slewing bearing

    圖4 試驗(yàn)前后轉(zhuǎn)盤軸承部件對(duì)比圖Fig.4 Contrast diagram of slewing bearing parts before and after experiment

    圖5 轉(zhuǎn)盤軸承的振動(dòng)信號(hào)Fig.5 Vibration signal of slewing bearing

    3.2 數(shù)據(jù)處理分析

    采用EEMD-相關(guān)系數(shù)對(duì)圖4所示信號(hào)進(jìn)行降噪處理,通過wavelet leader算法獲取降噪后信號(hào)多分形特征信息并重構(gòu)出特征向量組Ti。任選Ti中某3個(gè)特征作為x,y,z軸的數(shù)據(jù)進(jìn)行繪圖,結(jié)果如圖6所示,圖中 f1,f2,f3分別表示特征參數(shù)D(q)1,D(q)6,D(q)11。從圖中可以看出 3種故障狀態(tài)的區(qū)分效果并不明顯。將Ti特征集直接作為特征向量導(dǎo)入模型進(jìn)行訓(xùn)練,必定會(huì)對(duì)模型精度產(chǎn)生一定影響,同時(shí)多維特征向量也會(huì)影響模型計(jì)算速度,因此,需要對(duì)特征集進(jìn)行降維操作,在確保不丟失原特征信息前提下消除干擾信息,以達(dá)到約簡(jiǎn)維數(shù)的目的。

    圖6 轉(zhuǎn)盤軸承原始特征集Fig.6 Original feature set of slewing bearing

    利用ISOMAP算法對(duì)提取的轉(zhuǎn)盤軸承振動(dòng)信號(hào)數(shù)據(jù)幾何特征信息進(jìn)行降維處理。選取近鄰參數(shù)k=7,利用低維嵌入剩余方差R確定特征信息的低維嵌入空間d。轉(zhuǎn)盤軸承特征信息剩余方差與低維嵌入維數(shù)之間的關(guān)系如圖7所示。從圖中可以看出,剩余方差隨低維嵌入維數(shù)的增加而降低,當(dāng)?shù)途S嵌入維數(shù)增加到3以后,剩余方差基本保持在0.023附近,由此確定低維嵌入維數(shù)取值為3。選取特征信息低維嵌入特征矩陣前3維特征并繪制出3種狀態(tài)識(shí)別圖,結(jié)果如圖8所示。圖中f1,f2,f3分別表示經(jīng)ISOMAP降維處理后前3維量綱一的數(shù)據(jù),與圖6對(duì)比可以看出,經(jīng)等距映射方法降維處理后的特征信息具有很好的識(shí)別度。

    圖7 剩余方差與維度的關(guān)系Fig.7 Relationship between residual variance and dimension

    圖8 經(jīng)ISOMAP降維處理后的轉(zhuǎn)盤軸承數(shù)據(jù)集Fig.8 Data set of slewing bearing after dimension reduction of ISOMAP

    為進(jìn)一步詳細(xì)判斷轉(zhuǎn)盤軸承幾何特征信息降維前后的識(shí)別效果,通過LSSVM數(shù)據(jù)統(tǒng)計(jì)模型進(jìn)行量化說明。將經(jīng)過等距映射方法處理前、后的特征向量組作為樣本集并分為訓(xùn)練樣本和測(cè)試樣本。訓(xùn)練樣本的輸入為提取的某一狀態(tài)特征量,輸出為該狀態(tài)對(duì)應(yīng)的故障類型函數(shù)值f。在進(jìn)行LSSVM訓(xùn)練建模之前,首先確立輸入樣本與輸出樣本之間的函數(shù)對(duì)應(yīng)關(guān)系,即函數(shù)值f與故障類型之間的關(guān)系為f={1—正常狀態(tài),2—螺栓破壞,3—外圈破壞}。測(cè)試樣本的輸入為未知狀態(tài)對(duì)應(yīng)的特征量,輸出為L(zhǎng)SSVM模型的預(yù)測(cè)f=selection{1,2,3}。ISOMAP降維前后各選取 450組樣本集,其中訓(xùn)練樣本集為300組,測(cè)試樣本集為150組。利用訓(xùn)練樣本對(duì)LSSVM模型進(jìn)行訓(xùn)練,然后用測(cè)試樣本檢驗(yàn)?zāi)P蛯?duì)故障類型識(shí)別的效果,最終分類結(jié)果如圖9所示,模型預(yù)測(cè)精度見表2。

    圖9 ISOMAP降維前后軸承故障類型識(shí)別結(jié)果Fig.9 Recognition results of fault type for bearingbefore and after dimension reduction of ISOMAP

    表2 降維前后識(shí)別結(jié)果Tab.2 Recognition results before and after dimension reduction

    綜合分析可知,未經(jīng)降維處理的故障類型整體識(shí)別率為84.67%,經(jīng)過ISOMAP降維處理的總體識(shí)別率為96.67%,經(jīng)過ISOMAP降維處理后的故障類型識(shí)別結(jié)果明顯優(yōu)于未經(jīng)降維處理的結(jié)果。另外,降維前、后正常狀態(tài)識(shí)別精度均為90%,螺栓破壞的識(shí)別精度均為100%,降維前對(duì)外圈破壞的識(shí)別準(zhǔn)確率為64%,而經(jīng)過ISOMAP降維后對(duì)外圈破壞的識(shí)別準(zhǔn)確率達(dá)到了100%,這是由于ISOMAP降維處理后的特征向量剔除了重復(fù)、干擾信息,保留了有效特征信息,進(jìn)而改進(jìn)了模型整體的故障識(shí)別效果。正常狀態(tài)的識(shí)別結(jié)果相比降維前識(shí)別精度未進(jìn)一步得到改進(jìn),其原因是由于構(gòu)建原始特征向量時(shí)僅通過提取特殊位置點(diǎn)的信息以表征信號(hào)的幾何結(jié)構(gòu)信息,雖然減少了工作量,同時(shí)可能造成原始特征信息不完整的問題。因此,如何深度挖掘振動(dòng)信號(hào)有效的幾何結(jié)構(gòu)特征需要進(jìn)一步研究。

    4 結(jié)束語

    以轉(zhuǎn)盤軸承為研究背景,利用wavelet leader方法計(jì)算振動(dòng)信號(hào)多分形特征,提取多分形特征幾何信息進(jìn)行特征重構(gòu)。針對(duì)高維特征包含干擾信息易對(duì)模型訓(xùn)練產(chǎn)生影響的問題,提出了等距映射特征降維方法,既保留了振動(dòng)信號(hào)數(shù)據(jù)中的幾何特征信息,又降低了特征維數(shù),減少了計(jì)算量。結(jié)合最小二乘支持向量機(jī)在數(shù)據(jù)統(tǒng)計(jì)分析方面的優(yōu)勢(shì),利用降維后的特征樣本集對(duì)LSSVM模型進(jìn)行訓(xùn)練和測(cè)試,診斷精度得到了有效提高,證明了等距映射特征降維方法在故障識(shí)別方面的可行性及優(yōu)勢(shì)。

    猜你喜歡
    降維分形軸承
    Three-Body’s epic scale and fiercely guarded fanbase present challenges to adaptations
    軸承知識(shí)
    哈爾濱軸承(2022年2期)2022-07-22 06:39:32
    軸承知識(shí)
    哈爾濱軸承(2022年1期)2022-05-23 13:13:24
    軸承知識(shí)
    哈爾濱軸承(2021年2期)2021-08-12 06:11:46
    軸承知識(shí)
    哈爾濱軸承(2021年1期)2021-07-21 05:43:16
    感受分形
    降維打擊
    海峽姐妹(2019年12期)2020-01-14 03:24:40
    分形之美
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應(yīng)用
    拋物化Navier-Stokes方程的降維仿真模型
    欧美成人午夜免费资源| 在线精品无人区一区二区三| 91精品国产国语对白视频| 99久久精品热视频| 人人妻人人澡人人爽人人夜夜| 国产免费又黄又爽又色| 建设人人有责人人尽责人人享有的| 男人爽女人下面视频在线观看| 国内精品宾馆在线| 久久久久网色| 少妇熟女欧美另类| 国产视频内射| 啦啦啦啦在线视频资源| 91精品国产九色| 久久久久久久国产电影| 91精品国产九色| 久久久久久人妻| 极品少妇高潮喷水抽搐| 国产爽快片一区二区三区| 色婷婷av一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片 在线播放| 亚洲av.av天堂| 国产极品粉嫩免费观看在线 | 亚洲怡红院男人天堂| a级毛片在线看网站| 男女免费视频国产| 三级经典国产精品| 少妇的逼水好多| 乱人伦中国视频| 人妻系列 视频| 亚洲综合精品二区| 久久97久久精品| 搡老乐熟女国产| 亚洲美女视频黄频| 大话2 男鬼变身卡| .国产精品久久| 熟女电影av网| 边亲边吃奶的免费视频| 精品人妻偷拍中文字幕| 精华霜和精华液先用哪个| 久久精品久久久久久噜噜老黄| 只有这里有精品99| 亚洲精品日韩在线中文字幕| 黄色日韩在线| 久久女婷五月综合色啪小说| 精品熟女少妇av免费看| 只有这里有精品99| 又大又黄又爽视频免费| 久久国产乱子免费精品| 午夜日本视频在线| 久久女婷五月综合色啪小说| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 国产熟女午夜一区二区三区 | 国产精品麻豆人妻色哟哟久久| 极品少妇高潮喷水抽搐| 国产色婷婷99| av线在线观看网站| av播播在线观看一区| 久久97久久精品| 免费av中文字幕在线| 美女xxoo啪啪120秒动态图| 啦啦啦在线观看免费高清www| 久久久久视频综合| 久久国产乱子免费精品| 欧美日韩亚洲高清精品| 日韩中文字幕视频在线看片| 男女啪啪激烈高潮av片| 国产熟女欧美一区二区| 亚洲av在线观看美女高潮| 热99国产精品久久久久久7| 丰满少妇做爰视频| 国产男女内射视频| 最近的中文字幕免费完整| 22中文网久久字幕| 亚洲三级黄色毛片| 成人二区视频| 亚洲高清免费不卡视频| 观看美女的网站| 精品少妇内射三级| 国产 精品1| 国产日韩欧美视频二区| a级毛色黄片| a级毛片免费高清观看在线播放| 国产精品国产三级专区第一集| 肉色欧美久久久久久久蜜桃| 免费人成在线观看视频色| 99久久综合免费| 日韩精品有码人妻一区| 国产精品国产av在线观看| 久久 成人 亚洲| 国产亚洲午夜精品一区二区久久| 久久久国产一区二区| 久久久久视频综合| 熟女电影av网| 一级毛片黄色毛片免费观看视频| 亚洲国产成人一精品久久久| 免费观看性生交大片5| 一级a做视频免费观看| 亚洲国产av新网站| 精品午夜福利在线看| 在线观看三级黄色| 一边亲一边摸免费视频| 亚洲av福利一区| 伦理电影免费视频| 妹子高潮喷水视频| 久久久精品免费免费高清| 亚洲精品aⅴ在线观看| 亚洲伊人久久精品综合| 有码 亚洲区| 亚洲人与动物交配视频| 日本91视频免费播放| 国产精品不卡视频一区二区| 日本av免费视频播放| 男女边吃奶边做爰视频| 国产亚洲一区二区精品| 男女无遮挡免费网站观看| 男女免费视频国产| 大码成人一级视频| 国产高清国产精品国产三级| 99久久综合免费| 少妇的逼好多水| 极品少妇高潮喷水抽搐| 下体分泌物呈黄色| 99九九在线精品视频 | 天堂俺去俺来也www色官网| 啦啦啦啦在线视频资源| 国产精品麻豆人妻色哟哟久久| 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 免费观看性生交大片5| 亚洲人与动物交配视频| 亚洲无线观看免费| 亚洲av中文av极速乱| 日韩三级伦理在线观看| 极品人妻少妇av视频| 天堂8中文在线网| 欧美精品国产亚洲| 亚洲av成人精品一二三区| 国产男人的电影天堂91| 午夜免费鲁丝| 性色avwww在线观看| 久久精品久久久久久久性| 亚洲av在线观看美女高潮| 成人免费观看视频高清| 午夜av观看不卡| a级毛片免费高清观看在线播放| 韩国高清视频一区二区三区| 黄色一级大片看看| 少妇人妻一区二区三区视频| av国产久精品久网站免费入址| 亚洲国产成人一精品久久久| 精品熟女少妇av免费看| 精品久久久久久电影网| 久久久久久久国产电影| 99久久精品一区二区三区| 边亲边吃奶的免费视频| 青青草视频在线视频观看| av女优亚洲男人天堂| 大又大粗又爽又黄少妇毛片口| 九草在线视频观看| 精品午夜福利在线看| 亚洲欧美清纯卡通| 黑人巨大精品欧美一区二区蜜桃 | 亚洲天堂av无毛| 色视频在线一区二区三区| 国产美女午夜福利| 99国产精品免费福利视频| 亚洲精品,欧美精品| 91在线精品国自产拍蜜月| 丰满饥渴人妻一区二区三| 欧美日本中文国产一区发布| av福利片在线| 2021少妇久久久久久久久久久| 五月天丁香电影| 国产日韩欧美在线精品| 国产午夜精品久久久久久一区二区三区| 97在线视频观看| 啦啦啦啦在线视频资源| 精品人妻熟女毛片av久久网站| av免费观看日本| 国产精品欧美亚洲77777| 妹子高潮喷水视频| 亚洲在久久综合| 交换朋友夫妻互换小说| 又粗又硬又长又爽又黄的视频| 一本久久精品| 午夜福利,免费看| 日本av免费视频播放| 好男人视频免费观看在线| 亚洲精品国产av蜜桃| 免费观看av网站的网址| 自拍欧美九色日韩亚洲蝌蚪91 | 成年人免费黄色播放视频 | 国产黄片视频在线免费观看| 国产视频内射| 国产精品久久久久久久电影| 97超视频在线观看视频| 国产精品久久久久久精品古装| 熟女av电影| 精品国产一区二区久久| 18禁在线无遮挡免费观看视频| 日本欧美视频一区| 内射极品少妇av片p| 久久久久久久久大av| 大片免费播放器 马上看| 99热这里只有是精品50| 久久热精品热| 日本av手机在线免费观看| 国产精品久久久久久久久免| 成年人免费黄色播放视频 | 一区在线观看完整版| 精品国产国语对白av| 国产深夜福利视频在线观看| 亚洲精品自拍成人| 色视频在线一区二区三区| 观看免费一级毛片| 永久网站在线| 久久6这里有精品| 中文字幕制服av| 亚洲欧美精品自产自拍| 久久久精品94久久精品| 久久精品国产鲁丝片午夜精品| 欧美激情国产日韩精品一区| 色吧在线观看| 国产老妇伦熟女老妇高清| 不卡视频在线观看欧美| 一区二区av电影网| 中文乱码字字幕精品一区二区三区| 日日啪夜夜爽| 男的添女的下面高潮视频| 夫妻午夜视频| 偷拍熟女少妇极品色| 人妻少妇偷人精品九色| 丁香六月天网| 夜夜看夜夜爽夜夜摸| 最新的欧美精品一区二区| 一级片'在线观看视频| 少妇裸体淫交视频免费看高清| 最近手机中文字幕大全| 特大巨黑吊av在线直播| 伊人久久精品亚洲午夜| 这个男人来自地球电影免费观看 | av福利片在线| 午夜av观看不卡| 国产综合精华液| 国内精品宾馆在线| 成年女人在线观看亚洲视频| 看免费成人av毛片| 婷婷色综合大香蕉| 大片免费播放器 马上看| 在线观看www视频免费| 日本免费在线观看一区| 校园人妻丝袜中文字幕| 国产日韩欧美亚洲二区| 午夜视频国产福利| 日韩一本色道免费dvd| 国产真实伦视频高清在线观看| 少妇被粗大猛烈的视频| 免费av中文字幕在线| 亚洲成人手机| 国产黄色视频一区二区在线观看| 久久久久久久大尺度免费视频| 亚洲丝袜综合中文字幕| 久久人妻熟女aⅴ| 欧美bdsm另类| 9色porny在线观看| 久久99热6这里只有精品| 秋霞伦理黄片| 日本免费在线观看一区| 十分钟在线观看高清视频www | 如日韩欧美国产精品一区二区三区 | 18+在线观看网站| 亚洲av日韩在线播放| 久久 成人 亚洲| 观看免费一级毛片| 观看av在线不卡| 99九九在线精品视频 | 桃花免费在线播放| 欧美xxxx性猛交bbbb| a级片在线免费高清观看视频| 女人精品久久久久毛片| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费 | 人人妻人人添人人爽欧美一区卜| 午夜免费观看性视频| 日日撸夜夜添| 国产日韩欧美视频二区| 纵有疾风起免费观看全集完整版| √禁漫天堂资源中文www| 久久人人爽人人片av| 日本免费在线观看一区| 国产综合精华液| 秋霞伦理黄片| 成人国产麻豆网| 丰满饥渴人妻一区二区三| 国产精品一区二区在线观看99| 中国国产av一级| 另类精品久久| av在线观看视频网站免费| 亚洲高清免费不卡视频| 黄色一级大片看看| 国产精品偷伦视频观看了| 高清不卡的av网站| 亚洲国产精品国产精品| 亚洲av欧美aⅴ国产| 日韩强制内射视频| 极品教师在线视频| 色94色欧美一区二区| 丰满乱子伦码专区| 亚洲欧美清纯卡通| 五月开心婷婷网| 26uuu在线亚洲综合色| 国产中年淑女户外野战色| av在线老鸭窝| 色5月婷婷丁香| 丝袜脚勾引网站| 狂野欧美激情性bbbbbb| 中文字幕制服av| av有码第一页| 汤姆久久久久久久影院中文字幕| 极品人妻少妇av视频| av卡一久久| 建设人人有责人人尽责人人享有的| 少妇精品久久久久久久| 日本欧美视频一区| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| 丰满人妻一区二区三区视频av| 午夜免费鲁丝| 国产日韩欧美视频二区| 老女人水多毛片| 制服丝袜香蕉在线| 草草在线视频免费看| 18禁在线无遮挡免费观看视频| 妹子高潮喷水视频| 欧美人与善性xxx| 男女免费视频国产| 国产成人免费观看mmmm| 人妻人人澡人人爽人人| 日韩电影二区| 欧美3d第一页| 一个人免费看片子| 在线精品无人区一区二区三| 亚洲天堂av无毛| 亚州av有码| 成年人午夜在线观看视频| 久久毛片免费看一区二区三区| 国产精品成人在线| av天堂久久9| 另类精品久久| 欧美激情极品国产一区二区三区 | 黄色怎么调成土黄色| 久热久热在线精品观看| 国产精品国产av在线观看| 国产真实伦视频高清在线观看| 最近中文字幕高清免费大全6| av国产久精品久网站免费入址| 亚洲精品日韩av片在线观看| 人人澡人人妻人| 熟女人妻精品中文字幕| 丝袜喷水一区| 亚洲精品456在线播放app| 人人妻人人爽人人添夜夜欢视频 | 观看免费一级毛片| 国产av码专区亚洲av| 三级经典国产精品| 99久久精品一区二区三区| a级毛片在线看网站| 亚洲色图综合在线观看| 黑人巨大精品欧美一区二区蜜桃 | 中文字幕久久专区| 国产精品偷伦视频观看了| 一级a做视频免费观看| 少妇熟女欧美另类| 内射极品少妇av片p| 日韩一本色道免费dvd| 国产精品秋霞免费鲁丝片| 亚洲av免费高清在线观看| 一级爰片在线观看| 国产精品无大码| 少妇猛男粗大的猛烈进出视频| 精品人妻一区二区三区麻豆| 男人舔奶头视频| 日韩强制内射视频| 高清在线视频一区二区三区| 妹子高潮喷水视频| av在线老鸭窝| 极品少妇高潮喷水抽搐| 岛国毛片在线播放| 日本与韩国留学比较| 国产午夜精品久久久久久一区二区三区| 成人国产麻豆网| 色婷婷久久久亚洲欧美| 亚洲图色成人| 人人妻人人看人人澡| 午夜av观看不卡| 亚洲精品国产av成人精品| 久久久久久久久久久丰满| 九草在线视频观看| 亚州av有码| 一本一本综合久久| 最近手机中文字幕大全| 国产无遮挡羞羞视频在线观看| 国产精品不卡视频一区二区| 一区二区av电影网| 中文乱码字字幕精品一区二区三区| 亚洲在久久综合| 男的添女的下面高潮视频| 久久久久久久久大av| 99久国产av精品国产电影| 在线亚洲精品国产二区图片欧美 | 国产精品欧美亚洲77777| 亚洲精品视频女| 欧美国产精品一级二级三级 | 亚洲四区av| 大片电影免费在线观看免费| 最后的刺客免费高清国语| 99久久人妻综合| 精品视频人人做人人爽| 欧美日韩在线观看h| 久久韩国三级中文字幕| av视频免费观看在线观看| 亚洲伊人久久精品综合| 亚洲av综合色区一区| 国产白丝娇喘喷水9色精品| 亚洲精品乱久久久久久| 欧美日韩亚洲高清精品| 欧美激情极品国产一区二区三区 | 91久久精品国产一区二区三区| 精品少妇内射三级| 亚洲欧美一区二区三区国产| 国精品久久久久久国模美| 99九九在线精品视频 | 黑人巨大精品欧美一区二区蜜桃 | 国产一区二区在线观看av| 日韩强制内射视频| 国精品久久久久久国模美| 亚洲欧美一区二区三区国产| 亚洲三级黄色毛片| 国产精品一区二区性色av| 久久青草综合色| 亚洲av综合色区一区| 国产精品99久久久久久久久| 欧美bdsm另类| 婷婷色麻豆天堂久久| 噜噜噜噜噜久久久久久91| 国产无遮挡羞羞视频在线观看| 亚洲av二区三区四区| 久久久久人妻精品一区果冻| 3wmmmm亚洲av在线观看| 久久99热这里只频精品6学生| 最近手机中文字幕大全| 欧美变态另类bdsm刘玥| 91aial.com中文字幕在线观看| 国产在视频线精品| 成年人免费黄色播放视频 | 最近最新中文字幕免费大全7| 韩国av在线不卡| 成人毛片a级毛片在线播放| 黑人高潮一二区| 国产中年淑女户外野战色| 欧美精品一区二区大全| 日日爽夜夜爽网站| 欧美变态另类bdsm刘玥| 一级毛片黄色毛片免费观看视频| 久久精品久久久久久噜噜老黄| 免费看日本二区| 人妻制服诱惑在线中文字幕| 少妇人妻一区二区三区视频| 国精品久久久久久国模美| 97在线人人人人妻| 亚洲欧美精品专区久久| 91在线精品国自产拍蜜月| 亚洲成人一二三区av| 午夜精品国产一区二区电影| 日韩av不卡免费在线播放| 在线观看三级黄色| 一本大道久久a久久精品| 国产爽快片一区二区三区| 大话2 男鬼变身卡| 国产日韩一区二区三区精品不卡 | 天天躁夜夜躁狠狠久久av| 91在线精品国自产拍蜜月| 麻豆精品久久久久久蜜桃| 久久久久久久久久成人| 日韩不卡一区二区三区视频在线| 国产片特级美女逼逼视频| 亚洲,欧美,日韩| 男女啪啪激烈高潮av片| 黄色怎么调成土黄色| 下体分泌物呈黄色| 国产成人精品福利久久| 99视频精品全部免费 在线| 亚洲精品乱久久久久久| 欧美日韩国产mv在线观看视频| 日韩欧美一区视频在线观看 | 日韩不卡一区二区三区视频在线| 国产视频内射| 久久久久久久大尺度免费视频| 国产精品蜜桃在线观看| 岛国毛片在线播放| 国产成人精品婷婷| 在线观看人妻少妇| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 精品99又大又爽又粗少妇毛片| 亚洲精品日韩av片在线观看| 亚洲精品国产av成人精品| 亚洲第一av免费看| 九九久久精品国产亚洲av麻豆| 在线观看www视频免费| 99精国产麻豆久久婷婷| 久久影院123| 97精品久久久久久久久久精品| 一级毛片黄色毛片免费观看视频| 日韩中文字幕视频在线看片| 曰老女人黄片| 久久狼人影院| 蜜臀久久99精品久久宅男| 久久久久久久精品精品| 妹子高潮喷水视频| 男的添女的下面高潮视频| 国产一区二区在线观看日韩| 曰老女人黄片| 高清黄色对白视频在线免费看 | www.色视频.com| 免费观看的影片在线观看| tube8黄色片| 不卡视频在线观看欧美| 一区二区av电影网| 最近的中文字幕免费完整| 国产精品麻豆人妻色哟哟久久| 我的老师免费观看完整版| 人人澡人人妻人| videos熟女内射| 天堂中文最新版在线下载| 啦啦啦视频在线资源免费观看| 日本欧美视频一区| 日韩不卡一区二区三区视频在线| 亚洲精品乱码久久久v下载方式| 久久国产精品大桥未久av | 久久久精品94久久精品| 极品少妇高潮喷水抽搐| 久久国产乱子免费精品| 欧美亚洲 丝袜 人妻 在线| 婷婷色综合大香蕉| 极品教师在线视频| 黄色日韩在线| 纵有疾风起免费观看全集完整版| 在线观看三级黄色| 能在线免费看毛片的网站| 亚洲国产精品999| 国产伦理片在线播放av一区| 久久久久久久久久久久大奶| 日本欧美视频一区| 国产av码专区亚洲av| 爱豆传媒免费全集在线观看| 国产中年淑女户外野战色| 女性生殖器流出的白浆| av又黄又爽大尺度在线免费看| 成年av动漫网址| 国产永久视频网站| 国产成人免费无遮挡视频| 精品少妇黑人巨大在线播放| 亚洲国产精品专区欧美| 一二三四中文在线观看免费高清| 国产乱来视频区| 波野结衣二区三区在线| 黄色怎么调成土黄色| 十八禁网站网址无遮挡 | 午夜激情久久久久久久| 精品人妻一区二区三区麻豆| 十分钟在线观看高清视频www | 老熟女久久久| 亚洲av日韩在线播放| 国产精品久久久久成人av| 日韩在线高清观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产爽快片一区二区三区| 在线观看av片永久免费下载| av一本久久久久| 性色av一级| 国产精品无大码| 久久人人爽人人片av| 亚洲欧美精品自产自拍| a级毛片在线看网站| 大香蕉97超碰在线| 在线天堂最新版资源| 国产在视频线精品| 国语对白做爰xxxⅹ性视频网站| h视频一区二区三区| 国内揄拍国产精品人妻在线| 久久精品久久久久久久性| 最新中文字幕久久久久| 成人二区视频| 欧美日韩一区二区视频在线观看视频在线| 嫩草影院新地址| 黄色欧美视频在线观看| 日韩不卡一区二区三区视频在线| 黑人巨大精品欧美一区二区蜜桃 | 七月丁香在线播放| 色视频www国产| 欧美精品一区二区大全| av女优亚洲男人天堂| 欧美老熟妇乱子伦牲交| 街头女战士在线观看网站| videossex国产| 在线天堂最新版资源| 亚洲欧洲精品一区二区精品久久久 | 精品少妇黑人巨大在线播放| 亚洲情色 制服丝袜| 精品酒店卫生间| 色婷婷av一区二区三区视频|