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

    基于譜估計與核模糊聚類的往復壓縮機軸承故障評估方法

    2022-09-15 06:26:34劉巖王金東趙海洋王斌武韓興國
    機床與液壓 2022年12期
    關鍵詞:分形軸承聚類

    劉巖,王金東,趙海洋,王斌武,韓興國

    (1.桂林航天工業(yè)學院能源與建筑環(huán)境學院,廣西桂林 541004;2.東北石油大學機械科學與工程學院,黑龍江大慶 163318 )

    0 前言

    在能源動力行業(yè)中,往復式壓縮機是廣泛用于壓縮和輸送高危險介質的核心設備,其滑動軸承與組件既能進行旋轉運動又能進行往復運動,故障的隱蔽性、危害性極大,且表現(xiàn)出軸瓦(包括十字頭)磨損的多特征耦合模式?;谥鲃泳S護的評估與預示技術研究一直是故障診斷技術的延伸方向和技術攻關的重點。特別是自美國智能維護系統(tǒng)中心(IMS Center)提出設備性能退化評估的概念以來,壓縮機運動副間隙或摩擦性能的劣化所表現(xiàn)出的性能衰退問題已成為設備壽命評估研究的熱點與難點,諸多滾動軸承的性能退化評估模型被相繼提出。WANG等針對大規(guī)模的拉普拉斯矩陣分解問題,引入隨機相空間重構策略,提出了一種改進的圖譜重構方法來增強信號的特征,提取出滾動軸承退化數(shù)據(jù)的本征退化流形,并進行了驗證。LI等利用分段隱馬爾可夫模型(HMMs),結合多元指數(shù)加權滑動平均(MEWMA)控制圖來識別劣化的軸承退化多元信號,實現(xiàn)區(qū)分穩(wěn)定信號和劣化信號的有效性能退化評估。柏林等人采用相似近鄰傳播(AP)聚類方法,結合自組織神經(jīng)網(wǎng)絡(SOM)和自適應混沌粒子群(ACPSO)算法對多維特征集進行聚類和篩選,減少多維特征集之間相關冗余信息對壽命預測的影響。以上研究和方法為往復機的滑動軸承性能退化研究提供了有益參考。

    然而,往復機軸承性能衰退表現(xiàn)出的強烈多特征耦合的強沖擊時變特性,使評估模型與指標的建立更具模糊性和不確定性。結合軸承間隙故障振動信號傳遞路徑復雜性并考慮統(tǒng)計參數(shù)的一致性和敏感性要求,本文作者提出VMD分解與多重分形奇異譜特征參數(shù)結合的精細多重分形分析思想。利用正交化奇異值分解(Singular Value Decomposition,SVD)方法處理奇異譜特征參數(shù),建立核模糊C均值聚類(Kernel Fuzzy C-means clustering,KFCM)與二叉樹SVM結合的評估模型,實現(xiàn)不同滑動軸承磨損故障特征的有效聚類,并為預測和故障分類提供歸一化指標。

    1 基于模態(tài)分解的譜估計算法

    由于多重分形譜參數(shù)具有明確的物理意義,是反映系統(tǒng)動力學特征的有效參數(shù),奇異譜估計也成為非線性動力學研究中最活躍的分支之一。特別是基于自適應分解算法和分頻技術,融合的數(shù)據(jù)處理技術,是展現(xiàn)非平穩(wěn)時變時間序列局部分形尺度、反映系統(tǒng)內(nèi)部運行狀態(tài)的有力工具。

    1.1 變分模態(tài)分解

    對任意信號(),分解可概括為兩個步驟:構造變分模型、求解變分模型。

    1.1.1 構造變分模型

    (1)定義限定帶寬本征模態(tài)函數(shù)(Band-Limited Intrinsic Mode Function,BLIMF)為

    ()=()cos[()]

    (1)

    (2)解析信號獲取單邊譜的構建與混頻:

    (2)

    (3)帶寬估計。通過2范數(shù)梯度的平方根對信號解調,可得模態(tài)函數(shù)帶寬;

    (4)模型優(yōu)化。引入約束條件,構造優(yōu)化變分模型如下:

    (3)

    其中,為BLIMF分量{}={,,…,}的個數(shù);{}={,,…,}是()中心頻率。

    1.1.2 求解變分模型

    (1)為得到以上變分模型的最優(yōu)解,VMD通過引入二次罰因子和Lagrange乘子()構建增廣拉格朗日函數(shù)({},{},),將式(3)變?yōu)闊o約束問題。其中,二次懲罰因子用于保證噪聲背景下重構信號的準確性,Lagrange乘子用于保證模型約束剛度,即:

    ({},{},)=

    (4)

    (2)用交替方向乘子法更新、{}和{},并尋求增廣拉格朗日函數(shù)的鞍點,變分模型的最優(yōu)值由方程(3)求得,因此輸入信號()分解成個分量。以上VMD算法流程如圖1所示。

    圖 1 VMD算法流程

    1.2 多重分形譜特征值估計

    分形譜直接算法是以尺度為的盒子覆蓋所研究的多重分形集,通過計算點落在第個盒子的概率()而構造一個測度族,即:

    (5)

    該多重分形集的豪斯道夫維數(shù)為

    (6)

    該分形集整體奇異性均值為

    (7)

    以文中研究的壓縮機某測點振動信號奇異譜特征為例,得到振動信號多重分形譜表現(xiàn)出復雜分形體特征的連續(xù)單峰圖像,即()是關于的凸函數(shù)分析4個形態(tài)參數(shù)、Δ、和Δ=()-(),結果如圖2所示。

    圖2 多重分形奇異譜特征參數(shù)

    各參數(shù)的定義與物理意義:定義為中心值,是信號的長程相關性的表征,值越大相關性越強;Δ為譜寬度,反映振動信號的波動程度,值越大波動越劇烈;為附近的譜曲線擬合值,稱其為對稱度,當它大于0時,曲線形狀左傾,這時對應奇異性弱,反之曲線奇異性強;Δ為峰值差,反映振動信號峰值大小所占的比例,其值小于0,則概率最大子集數(shù)目大于概率最小子集數(shù)目,反之亦然。以上參數(shù)無量綱。

    VMD與多重分形融合可構成精細多重分形譜估計算法,是系統(tǒng)隱含信息提取、有效表征和識別的有效手段,具體表現(xiàn)為

    (1)系統(tǒng)特征成分常表現(xiàn)為幅值和頻率隨時間變化的調頻調幅(AM-FM)信號,有堅實理論基礎的VMD算法能將復雜信號中具有不同中心頻率的模態(tài)有效篩分;

    (2)多重分形理論是刻畫非線性信號多層次信息的有效語言,以維數(shù)譜、信息熵和奇異譜等算法構成較完備混沌與分形理論體系;

    (3)VMD較強的抗噪能力為精細多重分形譜分析的準確性提供有效的算法支撐,同時,基于奇異譜參數(shù)的精細多重分形與其他非線性分類和識別算法的有機結合,豐富和延伸了其理論應用領域。

    2 基于SVD與VMD的譜參數(shù)優(yōu)選

    2.1 VMD分解與奇異譜參數(shù)優(yōu)化

    2.1.1 數(shù)據(jù)采集與譜特征指標分析

    以大慶天然氣分公司2D12型雙作用對動式往復壓縮機為研究對象,其現(xiàn)場監(jiān)測圖如圖3所示,主要參數(shù)為軸功率500 kW、排氣量70 m/min、活塞行程240 mm、電機轉速496 r/min。

    圖3 2D12往復式壓縮機現(xiàn)場監(jiān)測

    結合長期維修實踐經(jīng)驗并根據(jù)API-618標準和活塞式壓縮機使用技術手冊,以壓縮機傳動機構中曲軸與二級連桿大頭軸瓦的間隙狀態(tài)為研究對象,結合標準和手冊要求將軸瓦間隙分別調至0.15、0.25、0.32和 0.40 mm,即劃分為正常狀態(tài)、輕微磨損、中度磨損和重度磨損4種狀態(tài)。針對采集的30周期(采樣頻率50 kHZ)二級連桿軸瓦測點振動信號,提取多重分形奇異譜(Multifractal Singular Spectrum,MSS),特征值如圖4所示。

    圖4 不同條件下原始數(shù)據(jù)MSS特征值

    特征參數(shù)對狀態(tài)的敏感性是評價狀態(tài)特征提取方法的重要指標,MSS參數(shù)可全面、直觀地刻畫非線性系統(tǒng)內(nèi)部通過振動信號表現(xiàn)的奇異性。由圖4可知:譜中心值和譜寬Δ相對其余參數(shù)而言,穩(wěn)定與可分性表現(xiàn)較高;由于缺乏有效的信號特征增強和歸一化手段,各周期的MSS譜形態(tài)特征值波動劇烈,可分性較差,無法作為有效的性能評估參數(shù)。

    2.1.2 VMD分解與SVD歸一化

    將奇異譜算法用于衰退性能評估,要求長時間采樣分析具有一致性,因此SVD法的周期信號探測與基于信噪分離的降維運算能有效解決該問題。

    (1) SVD重構矩陣

    對于包含噪聲或突變信息的振動時間序列重構,其吸引子軌跡矩陣可以用矩陣加法表示:

    =++

    (8)

    其中:為特征值矩陣;為突變特征信息矩陣;為噪聲矩陣。若已知、未知,可以通過研究矩陣的奇異值而得到并去除,即實現(xiàn)特征增強。

    本文作者以長時采樣信號整周期截斷法,形成固定列寬的連續(xù)截斷型矩陣,采用SVD降噪以形成穩(wěn)定性更好的特征指標。

    (2)參數(shù)選擇與優(yōu)化

    為實現(xiàn)SVD與VMD結合的特征增強算法,需要解決兩個問題:

    ①固定列寬的長序列矩陣不同周期VMD分解偏差對結果的影響

    考慮到特征成分往往包含在高頻成分中,采用能量歸一化的方法,消除周期差異保證零偏移特性,歸一化公式為

    (9)

    式中:()表示第個分量的原始信號

    ②確定SVD分解特征增強算法逆運算個數(shù)

    提出MRmR法與中心差商法融合互驗證,實現(xiàn)VMD預設尺度下的降噪與降維運算。結合測點的鍵相信號,每行為一周期,對30個采樣周期構造連續(xù)截斷型矩陣,得到中心差商如圖5所示。

    圖5 中心差商法奇異值分布曲線

    由圖5可知:最優(yōu)SVD分解個數(shù),即最優(yōu)模態(tài)個數(shù)(=4)處出現(xiàn)局部峰值??梢?,基于預設尺度的VMD構造整周期截斷型矩陣,因事先明確了各個特征BLIMF模態(tài),使得VMD與SVD結合的降噪目的明確,結果的普適性更利于評估指標建立。

    利用上述算法優(yōu)化,重新提取圖4數(shù)據(jù)的MSS特征值,如圖6所示。

    圖6 圖4中數(shù)據(jù)經(jīng)算法處理后的MMS特征值

    由圖6可知:所提的特征增強算法提高了譜參數(shù)的平穩(wěn)性和可分性,特別是譜中心值和譜寬兩參數(shù),其不同故障程度間的特征辨識度明顯增強,為后續(xù)的性能評估提供特征識別的方法支撐。

    3 KFCM性能衰減模型

    3.1 核模糊C均值聚類算法

    KFCM算法在約束函數(shù)中引入模糊數(shù)學中的隸屬度,并通過核函數(shù),以隱式映射實現(xiàn)數(shù)據(jù)樣本在高維空間可分,具有較好的普適性。該算法原理如式(10)所示:

    (10)

    式中:={,…,}∈,表示樣本數(shù);={,…,},表示將其劃分為類;=1,2,…,;為第類的聚類中心,所有構成隸屬度矩陣。

    由Lagrange乘子法,目標函數(shù)值最小的條件為

    (11)

    (12)

    令高維空間中距離表示為(,),則核空間的歐氏距離為

    ()-2()()+()

    (13)

    KFCM算法可簡單描述為4個步驟:

    步驟1,初始化矩陣,設定迭代次數(shù)和閾值;

    步驟2,計算距離(,);

    步驟3,更新隸屬度矩陣;

    3.2 算法步驟和流程

    綜上,對軸承性能評估過程可歸結為:首先,通過長期周期采樣構造連續(xù)截斷型矩陣,利用VMD分解模態(tài)分量;然后,經(jīng)歸一化SVD分解求逆過程,再通過多重分形理論提取譜形態(tài)參數(shù);最后,以KFCM算法優(yōu)選特征值,形成穩(wěn)定且可分性良好的特征向量,為特征衰減指標的建立與評估構建了基本算法流程,如圖7所示。

    圖7 基于KFCM與譜指標的軸承性能退化評估流程

    (1) 以軸承測點數(shù)據(jù)構造×的連續(xù)截斷型矩陣,確定預分解尺度并利用VMD方法分解各行向量,得BLIMF分量矩陣;

    (2)經(jīng)SVD分解后,以min(,,…,)為奇異值的重構數(shù)降噪,實現(xiàn)模態(tài)特征增強;

    (3)計算多重分形奇異譜,優(yōu)選奇異譜參數(shù)特征值;

    (4)基于KFCM算法,找到聚類中心以形成該狀態(tài)譜特征指標;對各性能衰減工況按步驟(1)~(4)計算,形成不同特征譜聚類中心;

    (5)對軸承間隙故障按正常、輕度磨損、中度磨損和重度磨損的實測數(shù)據(jù)建立全壽命周期的狀態(tài)類別,建立基于不同譜參數(shù)的閾值指標;

    (6)長期采集監(jiān)測數(shù)據(jù),利用聚類中心的歐氏距離平均值,建立設備性能退化全壽命指標衰減評估模型。

    4 往復壓縮機軸承間隙故障性能評估實例

    4.1 譜估計指標聚類分析

    對第2.1節(jié)中的現(xiàn)場采集數(shù)據(jù),根據(jù)評估流程計算4種工況下(正常狀態(tài)、輕微磨損、中度磨損和重度磨損)的最佳VMD分解值分別為=2、=3、=4、=4。按照公式(9)進行能量歸一化,經(jīng)SVD分解重構得到特征增強矩陣,提取MSS的4種譜特征向量。在KFCM算法中,核函數(shù)取為高斯核函數(shù),形式為

    (,)=exp[-‖-‖2]

    式中:為核參數(shù),取1。聚類數(shù)目取4,加權指數(shù)=2,當相鄰迭代間隸屬度差小于1×10時算法終止。同時,為分析和比較MSS各特征參數(shù)的可分性,選擇最具穩(wěn)定性的譜中心值分別與其余2個參數(shù)構成二維聚類以進行對比。各訓練樣本特征向量聚類結果如圖8所示。

    圖8 算法處理前后譜參數(shù)聚類分析對比

    由圖8可知:處理后各參數(shù)的團聚性明顯增加,考慮到譜中心值和譜寬Δ的穩(wěn)定性較高,聚類中譜寬Δ和對稱比兩個參數(shù)的可分性好,宜優(yōu)選三者作為軸承性能衰退的評估指標。經(jīng)計算,聚類指標和各狀態(tài)的聚類中心如表1所示。

    表1 各工況下譜特征參數(shù)聚類中心

    4.2 基于FSMBTC的狀態(tài)評估

    SVM通過高維映射可分性將特征向量樣本集作概率統(tǒng)計分類,KFCM 是類故障模式分類問題,將二者結合的算法稱為模糊二叉樹SVM。結合第4.1節(jié)特征向量聚類的二分特性,可建立往復壓縮機故障狀態(tài)特征衰減指標評估模型。

    4.2.1 評估流程

    結合第3節(jié)所提出的評估流程,分別采集4種狀態(tài)數(shù)據(jù)各30組MSS特征值作為訓練樣本;以相同測點實測4種間隙程度各30組數(shù)據(jù),經(jīng)VMD分解重構后計算得到的MSS特征向量為測試樣本,即30組訓練樣本,30組測試樣本。具體評估算法步驟如下:

    (1)首先,計算學習樣本模糊聚類中心={,,,},得到4個聚類得中心值、、、,則每類對應一個聚類中心;利用模糊聚類將聚成兩類,設聚成一類,記類,?,、、記為類,?;將和對應的學習樣本分別置為正類與負類,∩??,∪?,即可構造二叉樹分類器SVM1;

    (2)將正樣本聚類成兩類和,同理把兩類聚類中心對應的樣本置為正類和負類,∩??,∪?,?,?,構造二叉樹分類器SVM2;

    (3)以此類推,構造子分類器SVM3,直至每類只含一個聚類中心,所有分類器構成了一個基于模糊聚類的二叉樹結構。

    具體算法流程如圖9所示。明顯看出對于類問題,僅需構建-1個二叉樹分類器,在一對剩余(OVR)算法下,需構建個二類分類,而用一對一(OVO)算法,則需(-l)/2個分類器,基于模糊聚類的二叉樹算法分類訓練和識別速度較高。

    圖9 基于FSMBTC的分類流程

    4.2.2 算法評估

    為對比和驗證譜估計聚類算法有效性并判斷FSMBTC分類效果,在相同模擬數(shù)據(jù)和分類問題的情況下,采用基于EMD法(其中EMD模態(tài)重構個數(shù)按文獻[8]提出的互相關法確定)和基于VMD法聚類;采用直接計算特征值與聚類中心的歐氏距離法(ED)和FSMBTC法進行辨識。不同算法識別結果比較和聚類效果指標分別如表2和表3所示。

    表2 不同算法識別結果比較

    表3 不同算法聚類效果指標比較

    5 結論

    (1)融合VMD模態(tài)分解與多重分形理論,提出精細多重分形譜算法模型,結合SVD分解重構技術歸一化并突出模態(tài)特征成分,獲取復雜非線性系統(tǒng)狀態(tài)空間特征,并提高狀態(tài)參數(shù)的可分性和穩(wěn)定性;

    (2)為提升穩(wěn)定的狀態(tài)特征指標,將多重分形譜特征向量引入核模糊聚類算法,建立基于譜參量的軸承間隙故障衰減性能指標,給出完整算法流程;

    (3)引入模糊二叉樹SVM分類器,結合壓縮機軸承磨損全壽命狀態(tài)模擬,提高了整體算法模型對軸承磨損間隙衰退特征指標的有效率識別,為壽命預測提供參考。

    猜你喜歡
    分形軸承聚類
    軸承知識
    哈爾濱軸承(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
    感受分形
    分形之美
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應用
    基于DBSACN聚類算法的XML文檔聚類
    電子測試(2017年15期)2017-12-18 07:19:27
    基于改進的遺傳算法的模糊聚類算法
    一種層次初始的聚類個數(shù)自適應的聚類方法研究
    两个人的视频大全免费| 国产免费av片在线观看野外av| 日日啪夜夜撸| 男女下面进入的视频免费午夜| 国产主播在线观看一区二区| 少妇丰满av| 久久精品久久久久久噜噜老黄 | 两个人的视频大全免费| 国产中年淑女户外野战色| 亚洲熟妇熟女久久| 无遮挡黄片免费观看| 日韩人妻高清精品专区| 国内精品一区二区在线观看| 国产精品亚洲美女久久久| 精品人妻一区二区三区麻豆 | 色哟哟·www| 国产一区二区三区av在线 | 日韩中文字幕欧美一区二区| 内射极品少妇av片p| 国产成人a区在线观看| 亚洲一区高清亚洲精品| 免费大片18禁| 欧美日韩中文字幕国产精品一区二区三区| 久久人人爽人人爽人人片va| 国产蜜桃级精品一区二区三区| 国产精品久久久久久精品电影| 久久99热这里只有精品18| 热99在线观看视频| 国产不卡一卡二| 国产精品av视频在线免费观看| 国产精品无大码| 村上凉子中文字幕在线| 午夜福利18| 久久人人精品亚洲av| 嫩草影视91久久| 国产免费一级a男人的天堂| 麻豆久久精品国产亚洲av| 嫩草影院新地址| 久久久久久大精品| 亚洲精品粉嫩美女一区| 亚洲国产精品sss在线观看| 日韩欧美 国产精品| av天堂中文字幕网| 日韩一区二区视频免费看| 中文字幕免费在线视频6| 中亚洲国语对白在线视频| 极品教师在线免费播放| 99久久精品热视频| 亚洲欧美日韩无卡精品| 日韩一本色道免费dvd| 欧美区成人在线视频| 欧美+日韩+精品| netflix在线观看网站| 成人亚洲精品av一区二区| 国产精品野战在线观看| 精品一区二区三区av网在线观看| 在线看三级毛片| 18禁黄网站禁片午夜丰满| а√天堂www在线а√下载| 精品久久久久久久久av| 日本精品一区二区三区蜜桃| 久久这里只有精品中国| 国产av在哪里看| 欧美一区二区精品小视频在线| .国产精品久久| 亚洲,欧美,日韩| 成人美女网站在线观看视频| 成熟少妇高潮喷水视频| 麻豆国产av国片精品| 麻豆国产av国片精品| a在线观看视频网站| 成年人黄色毛片网站| 色av中文字幕| 亚洲色图av天堂| 日日啪夜夜撸| 国产黄a三级三级三级人| 最近中文字幕高清免费大全6 | 在现免费观看毛片| 极品教师在线视频| 国产免费男女视频| 国产av在哪里看| 精品午夜福利在线看| 欧美激情在线99| 亚洲人与动物交配视频| 最近最新中文字幕大全电影3| 免费无遮挡裸体视频| 波多野结衣高清无吗| 草草在线视频免费看| 亚洲精品国产成人久久av| 久久久久国内视频| 国产久久久一区二区三区| 看免费成人av毛片| 日本与韩国留学比较| 舔av片在线| 麻豆国产av国片精品| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久免费视频| 伊人久久精品亚洲午夜| 婷婷精品国产亚洲av在线| 亚洲18禁久久av| 我要搜黄色片| 欧美丝袜亚洲另类 | 日韩av在线大香蕉| 天天躁日日操中文字幕| 麻豆精品久久久久久蜜桃| 国产色爽女视频免费观看| 内射极品少妇av片p| 1000部很黄的大片| 五月玫瑰六月丁香| 国产欧美日韩一区二区精品| 又爽又黄a免费视频| 看免费成人av毛片| 91av网一区二区| 久久婷婷人人爽人人干人人爱| 免费人成视频x8x8入口观看| 一级黄片播放器| 欧美极品一区二区三区四区| 人人妻人人看人人澡| 成年版毛片免费区| 欧美bdsm另类| 中文字幕熟女人妻在线| 成年人黄色毛片网站| 久久精品国产亚洲av天美| 亚洲欧美日韩高清在线视频| 国产免费男女视频| 国产午夜福利久久久久久| 国产视频一区二区在线看| 婷婷精品国产亚洲av在线| 国产国拍精品亚洲av在线观看| 午夜福利高清视频| 亚洲自偷自拍三级| 成人无遮挡网站| 高清毛片免费观看视频网站| 精品久久久久久久久av| 日韩 亚洲 欧美在线| 成年免费大片在线观看| 美女 人体艺术 gogo| 久久午夜福利片| 日本精品一区二区三区蜜桃| 五月玫瑰六月丁香| 亚洲av成人av| 久久精品国产自在天天线| 91av网一区二区| 成人国产麻豆网| 男女做爰动态图高潮gif福利片| 久久婷婷人人爽人人干人人爱| 免费电影在线观看免费观看| 欧美bdsm另类| 国产伦精品一区二区三区视频9| 他把我摸到了高潮在线观看| 丰满人妻一区二区三区视频av| 99九九线精品视频在线观看视频| 国产麻豆成人av免费视频| 国产久久久一区二区三区| 欧美日韩瑟瑟在线播放| 午夜免费成人在线视频| 天堂√8在线中文| 亚洲精华国产精华精| 欧美丝袜亚洲另类 | 久久午夜福利片| 一卡2卡三卡四卡精品乱码亚洲| 国产在线男女| 日韩欧美 国产精品| 男人舔女人下体高潮全视频| 欧美又色又爽又黄视频| 久久久久久大精品| 欧美xxxx性猛交bbbb| 老师上课跳d突然被开到最大视频| 男女那种视频在线观看| 我的老师免费观看完整版| 国产精品亚洲美女久久久| 少妇被粗大猛烈的视频| 69人妻影院| aaaaa片日本免费| 日日夜夜操网爽| 亚洲成人免费电影在线观看| 18禁黄网站禁片免费观看直播| 国产一区二区在线av高清观看| 嫁个100分男人电影在线观看| 成熟少妇高潮喷水视频| 欧美日韩综合久久久久久 | 国产av不卡久久| 中亚洲国语对白在线视频| 精品久久久久久成人av| 又紧又爽又黄一区二区| 日本五十路高清| 美女高潮喷水抽搐中文字幕| 久久中文看片网| 午夜老司机福利剧场| 女的被弄到高潮叫床怎么办 | 少妇熟女aⅴ在线视频| 国产午夜精品久久久久久一区二区三区 | 最近中文字幕高清免费大全6 | 欧美在线一区亚洲| 人妻久久中文字幕网| 波多野结衣巨乳人妻| 97碰自拍视频| а√天堂www在线а√下载| 日本撒尿小便嘘嘘汇集6| 亚洲三级黄色毛片| av.在线天堂| 精品乱码久久久久久99久播| 欧美黑人巨大hd| 国产久久久一区二区三区| 波野结衣二区三区在线| 国产高清激情床上av| 精品免费久久久久久久清纯| bbb黄色大片| 小蜜桃在线观看免费完整版高清| 一进一出抽搐gif免费好疼| 欧美黑人欧美精品刺激| .国产精品久久| 国内揄拍国产精品人妻在线| 嫩草影院入口| 桃色一区二区三区在线观看| 免费观看精品视频网站| 国产爱豆传媒在线观看| 国产精品人妻久久久影院| 一本一本综合久久| 国产精品不卡视频一区二区| av国产免费在线观看| 18禁在线播放成人免费| 美女高潮的动态| 国产午夜精品论理片| 精品人妻视频免费看| 国产高清有码在线观看视频| 婷婷六月久久综合丁香| 国产一区二区三区在线臀色熟女| 日韩av在线大香蕉| 男插女下体视频免费在线播放| 欧美色欧美亚洲另类二区| 成人精品一区二区免费| 日韩一本色道免费dvd| 日本与韩国留学比较| 国产黄a三级三级三级人| 在线观看美女被高潮喷水网站| 精品无人区乱码1区二区| 搡老熟女国产l中国老女人| 午夜免费成人在线视频| 狂野欧美激情性xxxx在线观看| xxxwww97欧美| 亚洲av二区三区四区| 亚洲无线观看免费| 国产一区二区三区视频了| 国产蜜桃级精品一区二区三区| 国产亚洲精品久久久com| 乱码一卡2卡4卡精品| 一区二区三区激情视频| 久久久久久久午夜电影| 成年免费大片在线观看| av在线老鸭窝| 九色国产91popny在线| 国产高清三级在线| 在线播放无遮挡| 婷婷亚洲欧美| 老师上课跳d突然被开到最大视频| 看片在线看免费视频| 成人永久免费在线观看视频| 免费看光身美女| 亚洲国产精品合色在线| 在线观看免费视频日本深夜| 两个人的视频大全免费| 国内少妇人妻偷人精品xxx网站| 十八禁国产超污无遮挡网站| 亚洲成人久久性| 狂野欧美激情性xxxx在线观看| 国产精品嫩草影院av在线观看 | www.www免费av| 成人鲁丝片一二三区免费| 别揉我奶头 嗯啊视频| 免费黄网站久久成人精品| 国产成人aa在线观看| 一卡2卡三卡四卡精品乱码亚洲| 人妻久久中文字幕网| 午夜日韩欧美国产| 久久久久免费精品人妻一区二区| 日本撒尿小便嘘嘘汇集6| 真实男女啪啪啪动态图| 一个人看视频在线观看www免费| 天堂av国产一区二区熟女人妻| 免费人成在线观看视频色| 久久久精品大字幕| 婷婷精品国产亚洲av在线| 搡老妇女老女人老熟妇| 欧美一级a爱片免费观看看| 日本熟妇午夜| 神马国产精品三级电影在线观看| 国产蜜桃级精品一区二区三区| 久久香蕉精品热| 国内精品一区二区在线观看| 91午夜精品亚洲一区二区三区 | 狠狠狠狠99中文字幕| 日本一二三区视频观看| 精品久久久久久成人av| 国产精品亚洲美女久久久| 国产熟女欧美一区二区| 99国产极品粉嫩在线观看| a级毛片a级免费在线| 亚洲国产精品成人综合色| 久99久视频精品免费| av在线老鸭窝| 亚洲最大成人中文| av在线蜜桃| 亚洲国产欧洲综合997久久,| 成人午夜高清在线视频| 午夜激情福利司机影院| 人人妻人人看人人澡| 欧美区成人在线视频| 日本五十路高清| 深爱激情五月婷婷| www.www免费av| 久99久视频精品免费| av在线老鸭窝| 亚洲无线观看免费| 最近最新中文字幕大全电影3| 成人国产综合亚洲| 亚洲国产精品sss在线观看| 我要看日韩黄色一级片| 桃红色精品国产亚洲av| 亚洲成a人片在线一区二区| 淫妇啪啪啪对白视频| 黄色一级大片看看| av在线蜜桃| 久久精品国产清高在天天线| 嫩草影院新地址| 免费观看的影片在线观看| 老司机福利观看| 精品久久久久久成人av| 免费在线观看日本一区| 婷婷亚洲欧美| 国产三级在线视频| 成人三级黄色视频| 国产欧美日韩精品亚洲av| 久久精品影院6| 嫩草影院入口| 桃色一区二区三区在线观看| 小蜜桃在线观看免费完整版高清| 成人毛片a级毛片在线播放| 亚洲精品日韩av片在线观看| 亚洲美女黄片视频| 久久久久久大精品| 久久国内精品自在自线图片| 亚洲久久久久久中文字幕| 天堂影院成人在线观看| 天堂av国产一区二区熟女人妻| 中文字幕免费在线视频6| 免费在线观看日本一区| 日本黄大片高清| 婷婷六月久久综合丁香| 色综合亚洲欧美另类图片| 在线观看午夜福利视频| 亚洲国产精品成人综合色| 久久午夜亚洲精品久久| 久久久久久久久久久丰满 | 女生性感内裤真人,穿戴方法视频| 久久久久久久精品吃奶| 国产日本99.免费观看| 日本 欧美在线| 亚洲欧美精品综合久久99| 日本熟妇午夜| 精品人妻偷拍中文字幕| 18+在线观看网站| 久久精品久久久久久噜噜老黄 | 成人午夜高清在线视频| 黄色丝袜av网址大全| 999久久久精品免费观看国产| 久久人人爽人人爽人人片va| 一个人免费在线观看电影| 亚洲自偷自拍三级| 亚洲av第一区精品v没综合| 麻豆av噜噜一区二区三区| 人妻制服诱惑在线中文字幕| 久久久久久伊人网av| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 不卡一级毛片| 少妇裸体淫交视频免费看高清| 午夜日韩欧美国产| 狂野欧美激情性xxxx在线观看| 国产 一区 欧美 日韩| 老熟妇仑乱视频hdxx| 日本 欧美在线| 日本在线视频免费播放| 12—13女人毛片做爰片一| 国产高清视频在线观看网站| 毛片一级片免费看久久久久 | 搡老熟女国产l中国老女人| 女同久久另类99精品国产91| 国产91精品成人一区二区三区| 亚洲精品久久国产高清桃花| 日本在线视频免费播放| 国产视频一区二区在线看| 日韩欧美三级三区| av女优亚洲男人天堂| 日韩一区二区视频免费看| 日本与韩国留学比较| 日本a在线网址| 久久精品国产99精品国产亚洲性色| 国产一级毛片七仙女欲春2| 99riav亚洲国产免费| 欧美最新免费一区二区三区| 91久久精品国产一区二区成人| 国产在线男女| 淫秽高清视频在线观看| 狂野欧美白嫩少妇大欣赏| 大型黄色视频在线免费观看| 欧美xxxx性猛交bbbb| 国产乱人伦免费视频| 日本免费a在线| 九九爱精品视频在线观看| 无遮挡黄片免费观看| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看 | 日韩国内少妇激情av| 色哟哟·www| 午夜久久久久精精品| 18禁黄网站禁片免费观看直播| 成人无遮挡网站| 亚洲三级黄色毛片| 亚洲一级一片aⅴ在线观看| 老司机午夜福利在线观看视频| ponron亚洲| 一本久久中文字幕| 欧美三级亚洲精品| 国模一区二区三区四区视频| 免费av观看视频| 99久久中文字幕三级久久日本| 久久草成人影院| 成人国产一区最新在线观看| 一区二区三区免费毛片| 午夜影院日韩av| 久久久久性生活片| 亚洲自拍偷在线| 蜜桃亚洲精品一区二区三区| 午夜a级毛片| 久久精品国产鲁丝片午夜精品 | 日本三级黄在线观看| av中文乱码字幕在线| 午夜精品久久久久久毛片777| 亚洲第一区二区三区不卡| 少妇猛男粗大的猛烈进出视频 | 亚洲美女视频黄频| 在现免费观看毛片| 看黄色毛片网站| 亚洲av成人av| 久久久久久久亚洲中文字幕| 国产国拍精品亚洲av在线观看| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 两人在一起打扑克的视频| 亚洲第一区二区三区不卡| 天天躁日日操中文字幕| 亚洲最大成人中文| 国产黄a三级三级三级人| 久久久久久久久中文| 免费高清视频大片| 亚洲精品色激情综合| 91精品国产九色| 欧美又色又爽又黄视频| 能在线免费观看的黄片| 三级国产精品欧美在线观看| 男人的好看免费观看在线视频| 一个人看视频在线观看www免费| 日韩强制内射视频| 99精品在免费线老司机午夜| 国产男人的电影天堂91| 国内久久婷婷六月综合欲色啪| 亚洲国产精品sss在线观看| 免费在线观看成人毛片| 日本成人三级电影网站| 亚洲av一区综合| 日韩在线高清观看一区二区三区 | 久久天躁狠狠躁夜夜2o2o| 久久久久免费精品人妻一区二区| 欧美一区二区国产精品久久精品| 不卡视频在线观看欧美| 一个人免费在线观看电影| 亚洲经典国产精华液单| 久久九九热精品免费| 久久精品91蜜桃| 亚洲精品久久国产高清桃花| 在线a可以看的网站| 毛片女人毛片| 欧美绝顶高潮抽搐喷水| 午夜亚洲福利在线播放| 日本一本二区三区精品| 嫩草影视91久久| 深爱激情五月婷婷| 国产av一区在线观看免费| 国产老妇女一区| 九九热线精品视视频播放| 91麻豆av在线| 亚洲乱码一区二区免费版| 在线观看午夜福利视频| 女的被弄到高潮叫床怎么办 | 一级毛片久久久久久久久女| 日日夜夜操网爽| 亚洲 国产 在线| 俄罗斯特黄特色一大片| 午夜精品一区二区三区免费看| 波多野结衣高清无吗| 99精品久久久久人妻精品| 欧美激情久久久久久爽电影| 日韩欧美免费精品| 欧美最黄视频在线播放免费| av在线天堂中文字幕| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 日韩在线高清观看一区二区三区 | 97热精品久久久久久| 亚洲第一区二区三区不卡| 给我免费播放毛片高清在线观看| 亚洲电影在线观看av| 精品无人区乱码1区二区| 此物有八面人人有两片| 内地一区二区视频在线| 少妇被粗大猛烈的视频| 老司机福利观看| 亚洲av美国av| 国内精品久久久久精免费| 欧美性感艳星| 在线观看66精品国产| 97热精品久久久久久| 日韩一本色道免费dvd| 中亚洲国语对白在线视频| 国产女主播在线喷水免费视频网站 | 国产三级在线视频| 亚洲一区高清亚洲精品| 亚洲av.av天堂| 成人欧美大片| or卡值多少钱| 一级黄色大片毛片| 免费搜索国产男女视频| 成人永久免费在线观看视频| 亚洲av免费在线观看| 亚洲精品456在线播放app | 少妇裸体淫交视频免费看高清| 99精品在免费线老司机午夜| 亚洲av成人av| 亚洲精品乱码久久久v下载方式| 国产一区二区三区av在线 | 九九爱精品视频在线观看| 97碰自拍视频| 桃色一区二区三区在线观看| 国产高潮美女av| 最近在线观看免费完整版| 亚洲专区国产一区二区| 91久久精品国产一区二区成人| 欧美日韩国产亚洲二区| or卡值多少钱| 中文字幕久久专区| 哪里可以看免费的av片| 高清在线国产一区| 欧美+亚洲+日韩+国产| 国产精品一区www在线观看 | 亚洲av美国av| 两人在一起打扑克的视频| 国产三级在线视频| 国产精品嫩草影院av在线观看 | 在线天堂最新版资源| 99riav亚洲国产免费| 黄色一级大片看看| 校园人妻丝袜中文字幕| 欧美最黄视频在线播放免费| 国产蜜桃级精品一区二区三区| 成人亚洲精品av一区二区| 最好的美女福利视频网| 别揉我奶头~嗯~啊~动态视频| 免费在线观看成人毛片| 日韩,欧美,国产一区二区三区 | 精品久久国产蜜桃| 久久国产精品人妻蜜桃| av国产免费在线观看| 日韩 亚洲 欧美在线| 亚洲成人免费电影在线观看| 亚洲精品粉嫩美女一区| 免费黄网站久久成人精品| 欧美又色又爽又黄视频| 亚洲欧美日韩高清专用| av天堂中文字幕网| www.www免费av| 亚洲色图av天堂| 最近最新免费中文字幕在线| 少妇熟女aⅴ在线视频| 精品免费久久久久久久清纯| 在线观看舔阴道视频| 免费看美女性在线毛片视频| 中文字幕高清在线视频| 乱人视频在线观看| 啦啦啦观看免费观看视频高清| 国产欧美日韩精品亚洲av| 亚洲天堂国产精品一区在线| 久久中文看片网| 婷婷丁香在线五月| 国产真实伦视频高清在线观看 | 免费看日本二区| 亚洲精品成人久久久久久| 麻豆国产av国片精品| 日本免费一区二区三区高清不卡| 午夜老司机福利剧场| 一区二区三区激情视频| 欧美成人a在线观看| 九九在线视频观看精品| 欧美成人一区二区免费高清观看| 小说图片视频综合网站| 欧美xxxx黑人xx丫x性爽| 欧美日韩综合久久久久久 | 亚洲精品国产成人久久av| 可以在线观看毛片的网站| 琪琪午夜伦伦电影理论片6080| 熟妇人妻久久中文字幕3abv| 国产欧美日韩一区二区精品| 97碰自拍视频| 又黄又爽又刺激的免费视频.| 少妇裸体淫交视频免费看高清|