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

    基于主曲線相似度的軸承健康狀態(tài)評估方法*

    2019-06-25 06:32:40尹愛軍梁子曉王冬磊
    振動、測試與診斷 2019年3期
    關(guān)鍵詞:流形軸承距離

    尹愛軍, 梁子曉, 張 波, 王冬磊

    (1.重慶大學(xué)機(jī)械傳動國家重點實驗室 重慶,400044) (2.中國石油西南油氣田分公司重慶氣礦工藝研究所 重慶,400021) (3.中國工程物理研究院化工材料研究所 成都,621900)

    引 言

    旋轉(zhuǎn)機(jī)械設(shè)備在現(xiàn)代工業(yè)中起著不可替代的作用[1]。滾動軸承是旋轉(zhuǎn)機(jī)械中重要的部件,是最容易受損的部件之一,它直接影響著旋轉(zhuǎn)機(jī)械的性能和運(yùn)行[2]。

    目前,滾動軸承健康狀態(tài)評估方法的研究受到了學(xué)者的廣泛關(guān)注。蘇文勝等[3]將軸承振動信號進(jìn)行小波包分解后作為滾動軸承故障特征。黃海鳳等[4]利用盲源分離后軸承振動信號的峭度值作為軸承性能評估的特征。張云強(qiáng)等[5]提出一種基于中心對稱局部二值模式的時頻特征提取方法。然而,以上單一特征或?qū)S承早期故障不敏感,或穩(wěn)定性差,難以實現(xiàn)滾動軸承狀態(tài)的準(zhǔn)確評估。為此研究人員綜合運(yùn)用時域、頻域及小波包時頻域的多個特征參數(shù)進(jìn)行狀態(tài)評估[6]。Shao等[7]提出一種基于局部線性嵌入的連續(xù)深度信念網(wǎng)絡(luò)的方法實現(xiàn)滾動軸承故障檢測。Tong等[8]提出了一種基于雙樹復(fù)小波包變換、改進(jìn)時間尺度分解和極值學(xué)習(xí)機(jī)結(jié)合的滾動軸承故障識別新方法。Sohaib等[9]提出了一種基于混合特征模型和深層學(xué)習(xí)的滾動軸承狀態(tài)評估方法。李巍華等[10]將小波包熵與高斯混合模型(Gaussian mixture model ,簡稱GMM)相結(jié)合,通過計算各個運(yùn)行狀態(tài)對應(yīng)的GMM相對基準(zhǔn)GMM模型的偏離程度來判斷軸承的實際退化程度。然而,高維特征集之間往往存在比較強(qiáng)的相關(guān)性,有大量的“冗余”信息,并且維數(shù)過高,影響后續(xù)的算法效率[11]。因此,基于流形學(xué)習(xí)主曲線將原始高維特征空間轉(zhuǎn)換到低維空間再進(jìn)行評價等方法得到相關(guān)研究[12-13],這類方法在保留原始數(shù)據(jù)有效信息量的基礎(chǔ)上,消除“冗余”信息,提高效率。

    曲線的相似性評價是計算機(jī)圖形學(xué)、模式識別中的中心問題之一,目前主要的判別方法有特征值法和相似性比較法,相似性比較法在一定程度上比特征值法效果好。歐氏距離、馬氏距離、Frechet距離都是評估曲線的相似度有效方法,F(xiàn)rechet距離曲線相似性評價表現(xiàn)具有良好的性能[14]。

    基于此,筆者提出了一種基于流形空間主曲線Frechet距離相似度的狀態(tài)評估方法。實驗結(jié)果表明,該方法克服了單一特征在進(jìn)行故障診斷時準(zhǔn)確率不高的問題,對設(shè)備的早期故障更加敏感,且能有效評估滾動軸承的性能退化程度,實現(xiàn)滾動軸承健康狀態(tài)的定量評估。

    1 高維特征流形空間主曲線提取

    1.1 高維特征提取

    時域特征、頻域特征和時頻域特征是軸承故障特征的3種類型。時域特征性能不夠穩(wěn)定; 而頻域方法作為一種信號整體變換方法,處理對于統(tǒng)計特征隨時間不斷變化的非平穩(wěn)軸承故障信號效果并不理想;時頻域特征只適用于對緩變信號的分析,作為一種非動態(tài)分析方法,對非平穩(wěn)信號和突變信號的分析有局限[15-18]。

    為克服單一特征在表征振動信號特性時存在的片面性及其評估能力的不足,綜合運(yùn)用時域、頻域及小波包時頻域的高維特征參數(shù)進(jìn)行評估,對滾動軸承全壽命周期振動信號提取時域統(tǒng)計指標(biāo)、頻域統(tǒng)計指標(biāo)以及小波包能量特征指標(biāo),構(gòu)造融合38個特征量高維特征集[19-20]。

    1.2 基于LE流形學(xué)習(xí)的降維

    拉普拉斯特征映射(LE)是一種經(jīng)典的流形學(xué)習(xí)算法[21],可以在獲得原始高維數(shù)據(jù)的低維流形的基礎(chǔ)上,同時也能對離群點保持較強(qiáng)的魯棒性。在高維空間中距離較近的點,在低維空間中也應(yīng)該距離較近,是一種典型的局部非線性流形學(xué)習(xí)方法。

    給定高維樣本集合X={x1,x2,…,xM}T,其中:M為樣本點個數(shù);X∈RD。 LE的目標(biāo)是找到一個最優(yōu)的流形嵌入。設(shè)點集Y={y1,y2,…,yM}T,其中:Y∈RL是X的低維流形表示,Y=F(X),F(xiàn)則表示空間轉(zhuǎn)換映射

    (1)

    上述最優(yōu)化問題進(jìn)一步轉(zhuǎn)化為求解

    ΓF=-div▽‖F(xiàn)‖

    (2)

    其中:Γ為拉普拉斯特征算子。

    若f0,f1,…,fM是上述方程的解,則稱其為拉普拉斯特征算子Γ的廣義特征向量,且其對應(yīng)的廣義特征值按升序排列為0=λ0≤λ1≤…≤λM, 從而最優(yōu)的嵌入映射可以表示為F={f1,f2,…,fM}T,因此,高維空間數(shù)據(jù)X可以映射到低維空間數(shù)據(jù)Y

    Y=F(X)={f1(X),f2(X),…,fM(X)}

    (3)

    1.3 軟-K主曲線提取

    主曲線是由Hastie在1984年提出[22],是流形學(xué)習(xí)的一個分支,其理論基礎(chǔ)是尋找嵌入高維空間的非歐氏低維流形。主曲線是通過數(shù)據(jù)分布“中間”的曲線[23],能真實反映數(shù)據(jù)的形態(tài),即數(shù)據(jù)集合的“骨架”主曲線是一條通過數(shù)據(jù)“中心”的自相合光滑曲線。相比主成分分析(principal component analysis,簡稱PCA)在線性表達(dá)上對數(shù)據(jù)集合的描述,主曲線可以更大程度地刻畫數(shù)據(jù)集的趨勢和分布。Verbeek在多邊形主曲線的基礎(chǔ)上提出了軟-K主曲線算法[24],采用局部主成分方法來形成K條線段,并根據(jù)光滑性來連接形成主曲線。相對于其他方法,軟-K主曲線存在且唯一,曲線估計無偏差,其計算過程較為簡捷。

    軟-K主曲線算法步驟主要分為5步。

    1) 初始化。計算數(shù)據(jù)集X={x1,x2,…,xn}的第1主成分,且設(shè)s1=3σ為初始線段(σ為第1主成分標(biāo)準(zhǔn)差),求出s1的Voronoi區(qū)域可表示為V1={x1,x2,…,xn};K=1。

    2) 產(chǎn)生新線段。計算點xk

    當(dāng)mind(xi,sj)-‖xi-xj‖2>0時,則

    g(xi,xj)=mind(xi,sj)-‖xi-xj‖2

    (4)

    當(dāng)mind(xi,sj)-‖xi-xj‖2≤0時,則

    g(xi,xj)=0

    (5)

    xk的Voronoi區(qū)域可表示為

    xk={x∈X|x-xk≤mind(xi,sj)}

    (6)

    計算Vk第1主成分,則sk=3σ;K=K+1;

    3) 線段調(diào)整。計算新Voronoi區(qū)域

    4) 構(gòu)造優(yōu)化。用K條線段構(gòu)造出哈密爾頓路徑,并使用2-opt的旅行商問題算法(traveling salesman problem,簡稱TSP)來優(yōu)化該路徑,其優(yōu)化目標(biāo)函數(shù)為

    其中:l為路徑長度。

    5) 計算評判函數(shù)

    (7)

    其中:l為步驟4中的哈密爾頓路徑;Δ(xi,l)為點xi到l最短距離的平方,通過該評判函數(shù)是否收斂來控制迭代結(jié)束。

    2 基于Frechet距離的主曲線相似性狀態(tài)評價

    2.1 基于Frechet距離的曲線相似性評價

    離散Frechet距離是在連續(xù)Frechet距離的基礎(chǔ)上的一種曲線相似度評估算法,具有良好的性能[25]。

    給定兩個曲線P={u1,u2,…,up},Q={v1,v2,…,vq},P和Q之間的一個單調(diào)配對L為一系列對應(yīng)點(ua1,vb1)(ua2,vb2),…,(uam,vbm),其中a1=1,b1=1,am=p,bm=q,對所有的i=1,2,…,m,當(dāng)且僅當(dāng)下面是3個中的1個成立

    (11)

    其中:dist(uai,vbi)為uai,vbi間的歐氏距離。

    曲線P和Q間的離散Frechet距離為

    (12)

    距離數(shù)值越小,表明曲線P和Q越相似。

    圖1 滾動軸承性能退化評估方法流程圖Fig.1 A flow chart for evaluating the performance degradation of rolling bearings

    圖2 滾動軸承加速度退化實驗系統(tǒng)Fig.2 Experimental system of accelerated degradation of rolling bearing

    2.2 滾動軸承性能退化評估流程

    以開始階段的健康樣本主曲線作為基準(zhǔn)主曲線,通過計算各階段樣本主曲線和基準(zhǔn)主曲線的離散Frechet距離,實現(xiàn)對軸承退化的定量評估,流程如圖1所示。

    3 對比實驗分析

    3.1 滾動軸承健康狀態(tài)評估實驗系統(tǒng)

    實驗數(shù)據(jù)來自于美國辛辛那提大學(xué)IMS實驗室,軸承實驗裝置如圖2所示,4個ZA-2115雙列滾子軸承安裝在測試臺,由交流電機(jī)通過摩擦帶輪驅(qū)動旋轉(zhuǎn),轉(zhuǎn)速恒定為2 kr/min,在測試臺橫梁施加2 722 kg徑向載荷,以磁性塞累積磨損碎片反映軸承退化程度,當(dāng)累積碎片量超過一定閾值時,停止實驗。歷時9 840 min后,軸承1出現(xiàn)外圈故障。實驗采集軸承座的加速度振動信號,采樣頻率為20 kHz。

    3.2 滾動軸承健康狀態(tài)評估

    滾動軸承性能退化狀態(tài)一般分為初始故障狀態(tài)、中度故障狀態(tài)以及深度故障狀態(tài)[19]。 針對實驗中滾動軸承全壽命周期數(shù)據(jù),選取前600 min作為健康樣本,600 min至實驗結(jié)束均為測試樣本,評估滾動軸承性能狀態(tài)變化。

    圖3分為軸承健康樣本與測試樣本的流形空間二維主曲線。由圖3 (a)和(b)可知,運(yùn)行早期,主曲線形狀與趨勢比較相似;而后期軸承會產(chǎn)生磨損等缺陷,因此,所提取的主曲線有明顯的差異,如圖3(c)所示。

    圖4為3種不同“距離”評價下所得到的軸承狀態(tài)評估曲線。從圖4中的3條軸承狀態(tài)評估曲線變化趨勢可看出,軸承從開始運(yùn)行到5 000min左右相當(dāng)長的一段時間內(nèi),“距離”評估值很小,接近于0,且趨勢變化平穩(wěn),表明滾動軸承性能狀態(tài)良好。 評估曲線a在5 200 min左右時,有明顯上升, 偏離正常狀態(tài),軸承開始進(jìn)入初始性能退化狀態(tài)。 曲線b,c分別在6 400和6 500 min左右才檢測到軸承的初始性能退化狀態(tài)。Frechet距離對軸承早期故障更敏感。

    圖4 不同“距離”狀態(tài)評估曲線對比Fig.4 Comparison of different “distance” state assessment curves

    在7 000 min左右,3條曲線均有非常明顯的上升趨勢,顯示滾動軸承開始進(jìn)入中度性能退化狀態(tài)。在7 000~9 500 min左右,曲線呈現(xiàn)下凹趨勢,實質(zhì)為出現(xiàn)故障而故障又被消除的過程[26]。 曲線b和c的距離評估指標(biāo)在下凹階段與初始故障狀態(tài)的評估值有重疊,非常容易產(chǎn)生狀態(tài)識別錯誤,而曲線a并無重疊,表明筆者所提方法對中度性能退化狀態(tài)評估更好。

    在9 500~9 800 min左右,曲線b,c中的距離評估指標(biāo)出現(xiàn)震蕩,雖然呈現(xiàn)上升趨勢,但與中度性能退化狀態(tài)相比并無明顯區(qū)別,無法實現(xiàn)對滾動軸承健康狀態(tài)的定量評估。而曲線a有明顯上升趨勢,表明滾動軸承進(jìn)入深度性能退化狀態(tài),此時軸承已經(jīng)基本失去工作特性,直至到 9 800 min左右最終完全損壞失效。由此可知,筆者所提距離評估指標(biāo)對軸承初始性能退化狀態(tài)較敏感;在故障點平滑的下凹階段與初始退化指標(biāo)無重疊現(xiàn)象,且誤估率很低;同時,對軸承深度故障及失效狀態(tài)能夠進(jìn)行準(zhǔn)確的評估。因此,本研究方法具有更好的評估效果。

    作為非線性降維算法的代表,局部線性嵌入 (locally linear embedding ,簡稱LLE)能夠使降維后的數(shù)據(jù)較好地保持原有流形結(jié)構(gòu),而線性判別分析 (linear discriminant analysis,簡稱LDA)是一種經(jīng)典的有監(jiān)督的線性降維算法。本研究方法與LLE及LDA對比如圖5所示。對比分析3種降維方法,LLE算法隨時間變化呈現(xiàn)出明顯的波動,對軸承振動信號比較敏感,能大致反映出軸承的運(yùn)行狀態(tài)與退化趨勢,但波動過大,不能準(zhǔn)確識別故障程度;而LDA算法對滾動軸承性能退化狀態(tài)的評估曲線對軸承故障不敏感。直到9 800min左右軸承完全失效,評估曲線才出現(xiàn)上升,評估效果并不理想。

    圖5 不同流行學(xué)習(xí)算法狀態(tài)評估曲線對比Fig.5 Comparison of different characteristic spatial state assessment curves

    3.3 對比分析

    基于概率以及機(jī)器學(xué)習(xí)的算法在機(jī)械設(shè)備故障診斷、狀態(tài)評估等領(lǐng)域得到較多研究[7-9]。本研究方法與HMM及DBN對比如圖6所示。其中: HMM選取隱含狀態(tài)數(shù)為2,混合高斯函數(shù)個數(shù)為3;DBN選取2層隱含層,各層節(jié)點個數(shù)為100。

    圖6 狀態(tài)評估方法實驗對比Fig.6 Experimental comparison of state assessment methods

    由圖6可以看出,在5 000~6 000min之間,歸一化幅值無明顯變化,表明HMM與DBN模型對軸承早期故障不敏感;在7 000min之后,只有波動,沒有類似主曲線評估算法中突變-回落-波動的過程,即無法準(zhǔn)確描述軸承故障加深逐漸磨平的運(yùn)行狀態(tài)。相比較而言,這兩種方法均能在一定程度上描述滾動軸承的狀態(tài)變化趨勢,但均沒有本研究基于主曲線相似度的方法更準(zhǔn)確,且不能發(fā)現(xiàn)軸承的早期故障。

    定量分析3種方法,為此分別比較了評估曲線值變化率,表1為曲線值變化率為2%,5%時3種方法所用的時間對比。

    表1 評估曲線值變化率對比

    Tab.1 Comparison of the rate of change of the value of the evaluation curve

    /%/minDBN/minHMM/min255 3206 5306 8806 9906 8907 010

    由表1可以看出,主曲線相似度方法可以更快達(dá)到2%,5%的變化率,可知主曲線相似度對軸承振動信號以及早期故障更敏感。

    4 結(jié)束語

    針對單一特征在表征振動信號特性時存在的片面性及其評估能力的不足,融合多分析域特征,提出一種基于流形空間主曲線相似度的滾動軸承健康狀態(tài)評估方法,該方法能夠有效描述滾動軸承全壽命周期的性能退化規(guī)律。所提方法克服了單一特征監(jiān)測結(jié)果的不穩(wěn)定,能夠?qū)L動軸承健康狀態(tài)進(jìn)行定量評估,性能退化趨勢更為直觀。相對于歐式距離等距離評估方法,對早期損傷更敏感。本研究所提方法僅涉及低維流形空間中主曲線的提取,后續(xù)研究將考慮更高維度流形空間。

    猜你喜歡
    流形軸承距離
    軸承知識
    哈爾濱軸承(2022年2期)2022-07-22 06:39:32
    軸承知識
    哈爾濱軸承(2022年1期)2022-05-23 13:13:24
    軸承知識
    哈爾濱軸承(2021年2期)2021-08-12 06:11:46
    軸承知識
    哈爾濱軸承(2021年1期)2021-07-21 05:43:16
    緊流形上的Schr?dinger算子的譜間隙估計
    迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
    Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
    算距離
    每次失敗都會距離成功更近一步
    山東青年(2016年3期)2016-02-28 14:25:55
    基于多故障流形的旋轉(zhuǎn)機(jī)械故障診斷
    在线观看免费日韩欧美大片 | 久久久久精品久久久久真实原创| 少妇人妻久久综合中文| 黄色毛片三级朝国网站 | 成人国产麻豆网| 伦理电影免费视频| 免费在线观看成人毛片| 91精品一卡2卡3卡4卡| 天堂中文最新版在线下载| 精品久久久久久久久亚洲| 久久影院123| 男人和女人高潮做爰伦理| 极品少妇高潮喷水抽搐| 欧美区成人在线视频| 欧美精品一区二区免费开放| 人人妻人人添人人爽欧美一区卜| 一级av片app| 亚洲精品一二三| 日韩人妻高清精品专区| 夫妻午夜视频| 最近2019中文字幕mv第一页| 国产精品人妻久久久久久| 精品国产一区二区三区久久久樱花| 99热国产这里只有精品6| 青青草视频在线视频观看| 在线免费观看不下载黄p国产| 美女中出高潮动态图| 国产精品久久久久久久电影| xxx大片免费视频| 人妻人人澡人人爽人人| 亚洲av成人精品一区久久| 人体艺术视频欧美日本| 在线观看三级黄色| 亚洲欧美一区二区三区国产| 亚洲欧美成人精品一区二区| av播播在线观看一区| 国产高清国产精品国产三级| 热re99久久精品国产66热6| 免费看不卡的av| 在线观看av片永久免费下载| 青春草亚洲视频在线观看| 少妇人妻精品综合一区二区| 亚洲无线观看免费| 国产精品成人在线| 国产中年淑女户外野战色| 伊人久久国产一区二区| 老司机影院毛片| 日韩不卡一区二区三区视频在线| 高清欧美精品videossex| 日日啪夜夜爽| 在线看a的网站| 黑人高潮一二区| 国产免费一区二区三区四区乱码| 丝袜喷水一区| 亚洲av欧美aⅴ国产| 免费久久久久久久精品成人欧美视频 | 日韩精品免费视频一区二区三区 | 日韩电影二区| 免费高清在线观看视频在线观看| 国产高清国产精品国产三级| 亚洲国产成人一精品久久久| 黑丝袜美女国产一区| 欧美bdsm另类| av卡一久久| 丝瓜视频免费看黄片| 久久久国产精品麻豆| 免费播放大片免费观看视频在线观看| 高清欧美精品videossex| 久久久久久久亚洲中文字幕| 日韩av不卡免费在线播放| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡 | 在线观看av片永久免费下载| 人人妻人人添人人爽欧美一区卜| 亚洲av在线观看美女高潮| 老司机影院毛片| 国产成人免费观看mmmm| av女优亚洲男人天堂| 秋霞伦理黄片| 亚洲精品视频女| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 人妻制服诱惑在线中文字幕| 九九久久精品国产亚洲av麻豆| av国产久精品久网站免费入址| 涩涩av久久男人的天堂| 日韩一本色道免费dvd| 夫妻午夜视频| 国产精品蜜桃在线观看| 夜夜看夜夜爽夜夜摸| 各种免费的搞黄视频| 人妻人人澡人人爽人人| 亚洲人成网站在线观看播放| 精品久久久久久久久亚洲| 热re99久久国产66热| av黄色大香蕉| 九九在线视频观看精品| √禁漫天堂资源中文www| 99热全是精品| 黄色日韩在线| 在线观看三级黄色| 春色校园在线视频观看| 国产乱来视频区| 亚洲av成人精品一区久久| 一级毛片aaaaaa免费看小| 日韩精品免费视频一区二区三区 | 美女大奶头黄色视频| 国产精品无大码| 高清av免费在线| 国产成人a∨麻豆精品| 青春草国产在线视频| 欧美精品高潮呻吟av久久| 亚洲欧美日韩东京热| 免费在线观看成人毛片| 一级毛片我不卡| 嘟嘟电影网在线观看| 国产精品不卡视频一区二区| 欧美精品一区二区大全| 18禁裸乳无遮挡动漫免费视频| 久久国产亚洲av麻豆专区| 99久国产av精品国产电影| 狠狠精品人妻久久久久久综合| 亚洲国产精品999| 精品久久久精品久久久| 国产亚洲精品久久久com| 亚洲图色成人| 十八禁网站网址无遮挡 | 国模一区二区三区四区视频| 丝袜在线中文字幕| 在线播放无遮挡| 成人毛片a级毛片在线播放| 七月丁香在线播放| 免费观看性生交大片5| 亚洲不卡免费看| 高清午夜精品一区二区三区| 久久久久人妻精品一区果冻| 老熟女久久久| 亚洲精品乱码久久久久久按摩| 久久人人爽av亚洲精品天堂| 91成人精品电影| 22中文网久久字幕| 国产免费视频播放在线视频| 日韩一区二区视频免费看| 亚洲真实伦在线观看| 久久青草综合色| 国产在线一区二区三区精| 大陆偷拍与自拍| 少妇的逼好多水| 韩国高清视频一区二区三区| 老司机影院毛片| 久久人人爽人人爽人人片va| 一级av片app| 欧美日韩综合久久久久久| 久久人人爽人人爽人人片va| 中文字幕制服av| 九九在线视频观看精品| 18禁在线播放成人免费| av线在线观看网站| 男女边摸边吃奶| 久久久久久久精品精品| 欧美97在线视频| 蜜臀久久99精品久久宅男| 午夜福利视频精品| 人人妻人人看人人澡| 国产熟女午夜一区二区三区 | 又爽又黄a免费视频| 两个人的视频大全免费| 午夜免费鲁丝| 一级毛片aaaaaa免费看小| 少妇人妻一区二区三区视频| 女性被躁到高潮视频| a级片在线免费高清观看视频| 夜夜看夜夜爽夜夜摸| 一级毛片电影观看| 中文在线观看免费www的网站| 天天躁夜夜躁狠狠久久av| 天堂8中文在线网| 国产乱人偷精品视频| 纵有疾风起免费观看全集完整版| av卡一久久| 免费人成在线观看视频色| 久久99蜜桃精品久久| videossex国产| 成人综合一区亚洲| 又大又黄又爽视频免费| 秋霞在线观看毛片| 亚洲成人av在线免费| 国国产精品蜜臀av免费| 在线亚洲精品国产二区图片欧美 | 少妇的逼好多水| 男人狂女人下面高潮的视频| 伊人久久精品亚洲午夜| 国产精品偷伦视频观看了| 99久久精品一区二区三区| 久久久久久久久久成人| 亚洲av欧美aⅴ国产| 日韩大片免费观看网站| 最后的刺客免费高清国语| 99热这里只有精品一区| 免费看光身美女| 日韩中字成人| 人妻人人澡人人爽人人| 久久久久久久久久成人| 夫妻午夜视频| 日本91视频免费播放| 国产亚洲5aaaaa淫片| 六月丁香七月| 精品少妇黑人巨大在线播放| √禁漫天堂资源中文www| 狠狠精品人妻久久久久久综合| 在线看a的网站| 最近中文字幕高清免费大全6| 日韩视频在线欧美| 色婷婷av一区二区三区视频| 亚洲精品中文字幕在线视频 | 自拍欧美九色日韩亚洲蝌蚪91 | 一本—道久久a久久精品蜜桃钙片| 亚洲精品成人av观看孕妇| 波野结衣二区三区在线| 免费大片18禁| av在线老鸭窝| 国产精品秋霞免费鲁丝片| 久久久久国产精品人妻一区二区| 夜夜看夜夜爽夜夜摸| 观看美女的网站| 一个人免费看片子| 美女主播在线视频| 最近2019中文字幕mv第一页| 欧美xxⅹ黑人| 王馨瑶露胸无遮挡在线观看| 人妻一区二区av| 成人国产av品久久久| 久久婷婷青草| 精品国产一区二区久久| 国产精品伦人一区二区| 成人毛片60女人毛片免费| 一二三四中文在线观看免费高清| 人妻一区二区av| xxx大片免费视频| 三上悠亚av全集在线观看 | 特大巨黑吊av在线直播| 在线观看免费日韩欧美大片 | 国产淫语在线视频| 桃花免费在线播放| 纯流量卡能插随身wifi吗| 久久免费观看电影| 久久国内精品自在自线图片| 亚洲成人av在线免费| 国产精品女同一区二区软件| 国产伦精品一区二区三区四那| 草草在线视频免费看| 国产乱来视频区| 成年美女黄网站色视频大全免费 | av线在线观看网站| 少妇被粗大猛烈的视频| 亚洲人与动物交配视频| 欧美精品一区二区大全| 99久国产av精品国产电影| 国产成人精品福利久久| 成人影院久久| 成人美女网站在线观看视频| 精品一区二区免费观看| 欧美3d第一页| 午夜日本视频在线| 亚洲av免费高清在线观看| 黄色一级大片看看| 人人澡人人妻人| 午夜免费观看性视频| 日韩制服骚丝袜av| 午夜精品国产一区二区电影| 少妇人妻一区二区三区视频| 天堂俺去俺来也www色官网| 成年美女黄网站色视频大全免费 | 内地一区二区视频在线| 全区人妻精品视频| 午夜福利网站1000一区二区三区| 亚洲成色77777| 亚洲国产精品专区欧美| 国产精品久久久久成人av| 国产淫语在线视频| 亚洲熟女精品中文字幕| 男人和女人高潮做爰伦理| 18+在线观看网站| 国内少妇人妻偷人精品xxx网站| 久久久精品94久久精品| 国产亚洲av片在线观看秒播厂| 国产熟女午夜一区二区三区 | 伊人久久国产一区二区| 亚洲美女视频黄频| 高清毛片免费看| 少妇人妻一区二区三区视频| 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 老司机影院成人| 黄色日韩在线| 美女内射精品一级片tv| 国产欧美亚洲国产| 国产av国产精品国产| 国产午夜精品一二区理论片| 777米奇影视久久| 黑人巨大精品欧美一区二区蜜桃 | 欧美激情极品国产一区二区三区 | 久久久久久久久久久久大奶| 国产伦精品一区二区三区视频9| 精品卡一卡二卡四卡免费| 久久久久久久久久久丰满| 如日韩欧美国产精品一区二区三区 | 噜噜噜噜噜久久久久久91| 午夜影院在线不卡| 国产一区亚洲一区在线观看| 午夜影院在线不卡| 免费不卡的大黄色大毛片视频在线观看| 中文字幕制服av| 日韩av免费高清视频| 亚洲人成网站在线播| 自拍偷自拍亚洲精品老妇| 乱系列少妇在线播放| 亚洲欧美一区二区三区黑人 | 肉色欧美久久久久久久蜜桃| 午夜91福利影院| 精品人妻熟女毛片av久久网站| 激情五月婷婷亚洲| 免费看不卡的av| √禁漫天堂资源中文www| 久久国产精品男人的天堂亚洲 | 亚洲av中文av极速乱| 国产精品免费大片| 少妇被粗大猛烈的视频| 在线观看免费高清a一片| 亚洲av在线观看美女高潮| 三级国产精品欧美在线观看| 91久久精品国产一区二区三区| 晚上一个人看的免费电影| 欧美日韩一区二区视频在线观看视频在线| 久久人人爽人人爽人人片va| 香蕉精品网在线| 亚洲欧美一区二区三区黑人 | 夫妻性生交免费视频一级片| 成人无遮挡网站| 久久久国产精品麻豆| 日日爽夜夜爽网站| 国产精品一区二区三区四区免费观看| 久久 成人 亚洲| 99re6热这里在线精品视频| 18禁在线播放成人免费| 久久精品国产亚洲av涩爱| 伦理电影免费视频| 色5月婷婷丁香| 久久久久人妻精品一区果冻| 亚洲性久久影院| 99九九线精品视频在线观看视频| 18+在线观看网站| 91久久精品国产一区二区三区| 久热这里只有精品99| 亚洲内射少妇av| 久久精品熟女亚洲av麻豆精品| 各种免费的搞黄视频| 丝袜喷水一区| 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看 | 亚洲欧美日韩东京热| 视频中文字幕在线观看| 免费看光身美女| 日本爱情动作片www.在线观看| 人妻一区二区av| 久久综合国产亚洲精品| 国产成人freesex在线| 校园人妻丝袜中文字幕| 五月开心婷婷网| 啦啦啦视频在线资源免费观看| 欧美日本中文国产一区发布| 91精品伊人久久大香线蕉| 伦理电影大哥的女人| 国产免费视频播放在线视频| 色网站视频免费| 一区二区三区乱码不卡18| 人体艺术视频欧美日本| 一本一本综合久久| 日韩av不卡免费在线播放| 少妇人妻精品综合一区二区| 亚洲欧美精品专区久久| 久久久久久久国产电影| 亚洲精品日韩在线中文字幕| 人人澡人人妻人| 国产高清国产精品国产三级| 久久久久久久国产电影| 亚洲精品aⅴ在线观看| 一区二区三区四区激情视频| 精品久久久噜噜| 婷婷色综合大香蕉| av女优亚洲男人天堂| 大陆偷拍与自拍| 少妇裸体淫交视频免费看高清| 中文字幕久久专区| 婷婷色综合大香蕉| 国产黄色视频一区二区在线观看| 久久久国产精品麻豆| 六月丁香七月| 久久人妻熟女aⅴ| 在现免费观看毛片| 男女边摸边吃奶| 免费在线观看成人毛片| 国产熟女午夜一区二区三区 | 精品一品国产午夜福利视频| 丰满乱子伦码专区| 亚洲,欧美,日韩| 日韩成人av中文字幕在线观看| 久久鲁丝午夜福利片| 日日摸夜夜添夜夜爱| 久久久久久久久久久免费av| 国产精品久久久久成人av| 我的老师免费观看完整版| 亚洲国产精品成人久久小说| 只有这里有精品99| 国产91av在线免费观看| 97超视频在线观看视频| h视频一区二区三区| 亚洲精品国产成人久久av| 午夜福利,免费看| av国产精品久久久久影院| 精品酒店卫生间| 久久午夜福利片| 在线天堂最新版资源| 久久久久网色| 日日啪夜夜撸| 中文字幕免费在线视频6| 精品久久久久久电影网| 亚洲三级黄色毛片| 丰满迷人的少妇在线观看| 国产成人freesex在线| 美女福利国产在线| 日本-黄色视频高清免费观看| 免费观看的影片在线观看| 色94色欧美一区二区| 亚洲情色 制服丝袜| 涩涩av久久男人的天堂| 日韩人妻高清精品专区| 一二三四中文在线观看免费高清| 秋霞伦理黄片| 中国三级夫妇交换| 99热这里只有精品一区| 中文欧美无线码| 麻豆成人午夜福利视频| 国产熟女欧美一区二区| av卡一久久| 内射极品少妇av片p| 在线观看免费高清a一片| 全区人妻精品视频| 久久久久国产精品人妻一区二区| 啦啦啦啦在线视频资源| 在线观看免费日韩欧美大片 | 精品人妻一区二区三区麻豆| 亚洲一级一片aⅴ在线观看| 国产日韩欧美在线精品| 一级毛片 在线播放| 久久久久久人妻| 色吧在线观看| 80岁老熟妇乱子伦牲交| 亚洲欧美日韩卡通动漫| 国产一区二区三区综合在线观看 | 2021少妇久久久久久久久久久| 成年人免费黄色播放视频 | 岛国毛片在线播放| 少妇人妻 视频| 色视频www国产| 亚洲四区av| 亚洲成人av在线免费| 日日摸夜夜添夜夜爱| 国产一区二区三区综合在线观看 | 日韩一区二区视频免费看| 亚洲国产av新网站| 丝袜在线中文字幕| 国内少妇人妻偷人精品xxx网站| 高清欧美精品videossex| 国产高清不卡午夜福利| 国产成人免费无遮挡视频| 亚洲国产日韩一区二区| 天堂8中文在线网| 亚洲精华国产精华液的使用体验| 亚洲国产欧美在线一区| 18+在线观看网站| 一二三四中文在线观看免费高清| 亚洲av日韩在线播放| 韩国高清视频一区二区三区| 少妇人妻久久综合中文| 成人午夜精彩视频在线观看| av在线播放精品| 免费观看在线日韩| 韩国高清视频一区二区三区| 亚洲av男天堂| 80岁老熟妇乱子伦牲交| 美女福利国产在线| 男的添女的下面高潮视频| 一级av片app| 丰满人妻一区二区三区视频av| 久久毛片免费看一区二区三区| 高清不卡的av网站| 亚洲国产精品一区二区三区在线| 国产精品一二三区在线看| 日韩不卡一区二区三区视频在线| 久久婷婷青草| 免费观看无遮挡的男女| 国产女主播在线喷水免费视频网站| 亚洲成人一二三区av| 亚洲av男天堂| 中国美白少妇内射xxxbb| 国产免费一区二区三区四区乱码| 亚洲激情五月婷婷啪啪| 黑人巨大精品欧美一区二区蜜桃 | 一本—道久久a久久精品蜜桃钙片| 亚洲人成网站在线观看播放| 美女内射精品一级片tv| 亚洲精品自拍成人| 亚洲国产精品国产精品| 日韩强制内射视频| 三级经典国产精品| 日韩 亚洲 欧美在线| 少妇人妻一区二区三区视频| 国产日韩欧美在线精品| 国产黄片视频在线免费观看| 欧美xxⅹ黑人| 91精品国产九色| 18禁裸乳无遮挡动漫免费视频| 日日摸夜夜添夜夜添av毛片| 丰满少妇做爰视频| 久久久亚洲精品成人影院| 在线亚洲精品国产二区图片欧美 | 狂野欧美激情性bbbbbb| 精品国产一区二区久久| 国产又色又爽无遮挡免| 午夜91福利影院| 汤姆久久久久久久影院中文字幕| 精品人妻一区二区三区麻豆| 亚洲av电影在线观看一区二区三区| 国产男女超爽视频在线观看| 日韩一区二区三区影片| 中国三级夫妇交换| 内射极品少妇av片p| 亚洲四区av| 亚洲精品日韩av片在线观看| 亚洲三级黄色毛片| 国产成人免费无遮挡视频| 久久久久久久久久久久大奶| 国产深夜福利视频在线观看| 狂野欧美激情性xxxx在线观看| 欧美日韩综合久久久久久| 少妇熟女欧美另类| 最近最新中文字幕免费大全7| 午夜福利,免费看| 久久国内精品自在自线图片| av免费观看日本| 日韩精品免费视频一区二区三区 | 亚洲伊人久久精品综合| 欧美变态另类bdsm刘玥| 成人无遮挡网站| 日本与韩国留学比较| 一级爰片在线观看| 日本爱情动作片www.在线观看| 国产成人精品无人区| 99热这里只有精品一区| 亚洲欧洲精品一区二区精品久久久 | 赤兔流量卡办理| 99九九在线精品视频 | 97在线人人人人妻| 色哟哟·www| 国产白丝娇喘喷水9色精品| 免费黄色在线免费观看| 中国三级夫妇交换| 日本黄色日本黄色录像| a 毛片基地| 一区二区av电影网| 中文资源天堂在线| 春色校园在线视频观看| 国产熟女午夜一区二区三区 | 亚洲av中文av极速乱| 丁香六月天网| 最近中文字幕高清免费大全6| 中文天堂在线官网| 国产深夜福利视频在线观看| 丰满迷人的少妇在线观看| 欧美xxⅹ黑人| 99热这里只有精品一区| 久久午夜综合久久蜜桃| 寂寞人妻少妇视频99o| 中文字幕制服av| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美日韩东京热| 欧美精品亚洲一区二区| 精品久久久久久电影网| 欧美+日韩+精品| 女性被躁到高潮视频| 搡女人真爽免费视频火全软件| 欧美老熟妇乱子伦牲交| 五月天丁香电影| 最新中文字幕久久久久| 在线播放无遮挡| 日韩亚洲欧美综合| 人人妻人人澡人人看| 欧美高清成人免费视频www| 欧美 亚洲 国产 日韩一| 国产极品粉嫩免费观看在线 | 国产高清不卡午夜福利| 国产精品久久久久久精品古装| 亚洲精品乱久久久久久| 国产在线免费精品| 日韩一区二区三区影片| 欧美精品国产亚洲| av天堂中文字幕网| 久久精品久久久久久久性| 大片电影免费在线观看免费| 老司机影院毛片| 国产一级毛片在线| 国产亚洲午夜精品一区二区久久| 精品国产露脸久久av麻豆| 国产精品蜜桃在线观看|