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

    集成多策略改進FCM算法的旋轉機械故障數(shù)據聚類分析研究

    2019-11-27 03:21鄧林峰張愛華趙榮珍
    振動工程學報 2019年5期
    關鍵詞:聚類分析故障診斷

    鄧林峰 張愛華 趙榮珍

    摘要: 針對旋轉機械故障數(shù)據聚類分析中的初始聚類中心不確定和孤立點敏感問題,提出了一種集成多策略改進的模糊C均值(FCM)聚類方法。首先以故障數(shù)據集的決策屬性為等價關系對數(shù)據集進行劃分,得到若干個由等價關系導出的等價類;然后以每個等價類為可行域,采用均值漂移方法搜索故障數(shù)據類中心;最后以搜索到的類中心為FCM算法的初始聚類中心,通過核技術計算故障數(shù)據樣本與相應類中心在高維特征空間中的歐氏距離,從而實現(xiàn)數(shù)據樣本相似性的有效度量,并完成故障數(shù)據的模糊聚類。通過標準數(shù)據集和旋轉機械故障數(shù)據集對方法的性能進行了驗證及比較分析。結果顯示,改進FCM算法的聚類性能相比傳統(tǒng)FCM算法的聚類性能得到了明顯提升,在收斂速度和聚類準確性兩個性能指標上,改進的FCM算法比FCM算法具有顯著優(yōu)勢。

    關鍵詞: 故障診斷; 旋轉機械; 聚類分析; 模糊C均值; 聚類性能

    中圖分類號: TH165+.3; TP391.4 文獻標志碼: A 文章編號: 1004-4523(2019)05-0918-09

    DOI:10.16385/j.cnki.issn.1004-4523.2019.05.020

    引 言

    故障診斷技術為保障旋轉機械安全可靠運行發(fā)揮著極其重要的作用[1]。然而由于旋轉機械故障振動信號具有明顯的非平穩(wěn)特性,且易受到噪聲干擾,使得部分故障樣本的特征值將發(fā)生偏移,從而脫離實際故障的特征分布區(qū)間,甚至與其他故障的特征數(shù)據交疊在一起,形成分布狀態(tài)十分復雜的故障數(shù)據,導致分析診斷故障的難度急劇增大[2]。

    為了獲取準確的故障特征信息,各種數(shù)據挖掘技術被應用于故障數(shù)據的處理[3-5]。聚類分析是進行數(shù)據挖掘的主要方法之一[5],作為一種優(yōu)良的聚類算法,模糊C均值(Fuzzy C-means, FCM)聚類算法受到旋轉機械故障數(shù)據分析研究的廣泛關注。Rodríguez Ramos等[6]將基于密度的模糊C均值聚類與核模糊C均值聚類算法相結合,對攪拌反應器故障進行分類識別。鄭直等[7]提出一種EEMD形態(tài)譜和核模糊C均值聚類集成的故障診斷新方法,實現(xiàn)了滾動軸承故障的有效診斷。李懷俊等[8]研究了基于核特征模糊聚類及模糊關聯(lián)熵的故障識別方法,可顯著提高故障數(shù)據的聚類效果,并識別齒輪故障。田再克等[9]提出了一種基于改進多重分形去趨勢波動分析和半監(jiān)督馬氏距離模糊C均值的液壓泵退化狀態(tài)識別方法,并利用液壓泵實測數(shù)據驗證了方法的有效性。

    上述研究使FCM聚類算法與故障診斷技術有了很好的結合,可有效識別故障狀態(tài)。但經研究發(fā)現(xiàn),該算法依然存在初始聚類中心不確定和孤立點敏感的缺陷[10],這將導致FCM聚類算法在處理復雜數(shù)據時容易出現(xiàn)收斂速度慢、準確性偏低的問題,從而影響故障診斷效率和精度。近年來,雖然針對FCM算法的上述缺陷開展了一些研究工作[5, 10],但算法性能的提升還存在一定空間,有待進一步對其改進以滿足應用需求。

    基于以上分析,本文提出了一種集成多策略改進的FCM(Improved FCM, IFCM)聚類算法。該算法采用建立聚類中心的可行域和均值漂移方法解決故障數(shù)據集初始聚類中心不確定問題,并通過核方法計算故障數(shù)據樣本在高維特征空間中的距離,以減小孤立樣本對聚類準確性的影響。最后利用UCI標準數(shù)據集和轉子-軸承系統(tǒng)故障數(shù)據集對所提出的改進FCM算法的有效性進行驗證,通過與FCM算法進行對比分析,進一步說明改進算法具有更好的聚類性能。

    1 模糊C均值聚類算法

    模糊C均值聚類算法的出發(fā)點是基于對目標函數(shù)的優(yōu)化,通過迭代算法尋找一組中心矢量,使各樣本到其距離平方和達到最小,即通過對平方誤差函數(shù)求極值來確定最佳模糊分類矩陣和聚類中心[11]。該算法的基本流程描述如下:

    在運行聚類算法時,可行域設置的越小,相應的計算復雜度就越低,但過小的可行域可能并不包含真實的聚類中心,將使搜索過程失效;反之,可行域設置的越大,搜索范圍也相應增大,從而增加了計算復雜度。由于一個FCM聚類模型可以看成由若干個相互獨立的噪聲聚類問題的子模型構成,且每個子模型與一個單獨的類別相對應,因此對包含C個類別的數(shù)據集進行FCM聚類就可以轉化為求解這C個獨立優(yōu)化子問題各自的解[12]。也就是說,這C個子問題與C個類之間具有一一對應關系。

    另一方面,以數(shù)據集的決策屬性(數(shù)據類別)作為數(shù)據集X上的一個等價關系(表示為D)可形成對X的一個劃分,并得到X上一個由等價關系D導出的等價類集合XD,且每個等價類與數(shù)據集的一個類別相對應。因此,若數(shù)據集X由包含C個類別的數(shù)據構成,通過數(shù)據的決策屬性對其進行劃分,則可得到相應的C個等價類,而每個等價類就包含了數(shù)據集中某一類別的所有元素。

    通過比較劃分和可行域兩個不同概念的定義可以發(fā)現(xiàn),利用決策屬性D對數(shù)據集X形成的劃分將產生一個關于X的子集族,這個子集族必然滿足可行域的條件,即XD實際上是數(shù)據集X的聚類中心可行域的一種特殊情形。顯然,若X包含C個類別,則XD中將具有C個等價類,這些等價類就可作為FCM聚類算法C個獨立優(yōu)化子問題的可行域。

    建立了FCM算法優(yōu)化子問題的可行域之后,在運行FCM算法時,就可設定每個子目標函數(shù)在相應的可行域內搜索最優(yōu)值。由于XD中任意兩個等價類的交集為空集,因此以X〖〗D作為可行域能夠避免FCM算法產生重合聚類。

    2.2 均值漂移算法

    均值漂移(Mean-Shift, MS)是一種有效的統(tǒng)計迭代算法,使每一個點都“漂移”到密度函數(shù)的局部極大值點,并利用核密度估計來表示特征空間數(shù)據的分布情況[13]。MS算法的基本原理如下。

    2.3 核方法的樣本相似性度量

    數(shù)據樣本之間的相似性度量是進行聚類分析的基礎,對于最終聚類結果的形狀具有很大影響。在FCM及大多數(shù)目前常用的聚類算法中,通常都利用歐氏距離在原始特征空間對樣本間的相似性進行度量。但由于原始特征空間數(shù)據分布的復雜性,樣本的相似性度量會產生較大誤差,導致聚類結果的精度偏低。

    核方法[14]可以通過某種非線性變換將低維空間中的樣本映射到高維空間中,然后在這個高維空間中對樣本的相似性進行度量,以提高樣本相似性度量的準確度。因此,通過引入核函數(shù)將樣本映射到高維空間并設計新的相似性度量函數(shù),應該可以實現(xiàn)提升FCM算法聚類性能的目的。

    3 改進FCM算法的性能驗證及應用

    為評價本文提出算法的聚類性能,分別使用標準數(shù)據和轉子-軸承系統(tǒng)故障實驗數(shù)據進行聚類分析,并與FCM算法的聚類結果進行比較,以驗證本文算法的有效性。在進行聚類分析時,對于標準數(shù)據和轉子故障實驗數(shù)據,同一算法均采用相同的參數(shù)。其中,F(xiàn)CM算法和改進FCM算法的迭代閾值都設為0.001,模糊聚類指數(shù)m都設為2[6, 11];MS算法的迭代閾值ε1設為0.001,核函數(shù)G選取高斯核,核半徑h取2[13];改進FCM算法的核函數(shù)K也采用高斯核,核參數(shù)σ取1.34[6];FCM 和改進FCM算法的最大迭代次數(shù)都設定為100。

    3.1 標準數(shù)據

    從UCI數(shù)據庫中選擇兩個被廣泛用于聚類分析的數(shù)據集Iris[15]和Soybean[16],對FCM算法和改進FCM算法的性能進行比較,以驗證改進FCM算法的有效性。如上兩個數(shù)據集的數(shù)據特性如表1所示。利用兩種聚類算法對這兩個數(shù)據集分別進行處理,其目標函數(shù)的迭代誤差如圖1所示。

    ? 表1中的Iris數(shù)據集具有4維特征,包含3個不同類別的數(shù)據,每一類數(shù)據都有50個樣本;Soybean數(shù)據集具有35維特征值,包含4個不同類別的數(shù)據,其中3個類別每一類都有10個樣本,剩下一類有17個樣本。為便于描述,Iris的3個類別將采用符號數(shù)字1,2,3進行表示;Soybean的4個類別將采用符號數(shù)字1,2,3,4進行表示。

    從圖1可見,對于數(shù)據集Iris和Soybean,利用兩種聚類算法處理時,對應的目標函數(shù)迭代誤差存在明顯差別。顯然,F(xiàn)CM聚類算法的目標函數(shù)迭代誤差在迭代初始階段較大,隨著迭代次數(shù)的增加,迭代誤差迅速減小并趨近于零;而改進FCM聚類算法的目標函數(shù)迭代誤差在迭代初始階段就很小。這主要是因為改進的FCM聚類算法在進行模糊聚類之前先利用MS算法搜索到了數(shù)據集的初始聚類中心,通常情況下,這些初始的聚類中心與最終的聚類中心已非常接近,所以改進的FCM聚類算法只需要很少的迭代次數(shù)即可收斂。因此,與FCM聚類算法相比,改進FCM聚類算法的數(shù)據處理效率更高,收斂速度更快。

    兩種算法對兩個標準數(shù)據集的聚類準確率及運行時間如表2所示;具體的聚類結果如表3所示。

    ? ? ?表2顯示,利用兩種算法分別對Iris和Soybean兩個數(shù)據集進行聚類分析時,改進FCM算法的聚類準確率都更高。雖然改進的FCM算法比FCM算法的運行時間增長了1倍左右,這主要是由于利用MS算法搜索數(shù)據的初始聚類中心,所以增加了算法的運行時間,但都沒有超過0.3 s,具有較快的數(shù)據處理速度。

    從表3可見,與FCM算法相比,改進FCM算法的聚類精度之所以得到提高,主要是其將FCM算法中聚類錯誤的部分數(shù)據樣本劃分到了正確的類別,從而增加了這些類別的正確樣本,與實際樣本數(shù)更加接近。雖然改進的FCM算法對于Iris數(shù)據集第2類數(shù)據樣本的正確劃分減少了2個,但第3類數(shù)據樣本的正確劃分增加了7個,從而使得Iris數(shù)據集的總體聚類準確率仍然提高了3.34%;對于Soybean數(shù)據集,其第3、第4類數(shù)據樣本的正確劃分數(shù)分別增加了1個和3個,從而將改進的FCM算法的聚類準確率提高了8.51%。

    3.2 轉子-軸承系統(tǒng)故障數(shù)據

    將本文提出的改進FCM聚類算法應用于圖2所示轉子-軸承系統(tǒng)故障數(shù)據的聚類分析。故障模擬實驗設置了質量不平衡、轉子不對中、動靜碰摩、軸承松動及正常共5種運行狀態(tài)。采集轉子轉速為3000 r/min時的振動信號,采樣頻率設置為5000 Hz,采樣點數(shù)確定為12500,每一種運行狀態(tài)都采集100個數(shù)據樣本。

    ? 由文獻[17-19]的研究結果可知,均方根、均方根頻率、頻率標準差可有效刻畫不同故障狀態(tài),故選取這3個特征用于構造故障數(shù)據集。轉子實驗臺上共設置了12個振動傳感器,即采集的每一個樣本包含著12個通道,而每個通道的信號都能產生上述3個特征,因此可構成36維特征的故障數(shù)據集,如表4所示。

    表4所示的故障特征數(shù)據由12個通道的振動信號計算得到,每一行代表著一個振動數(shù)據樣本,也就是說,每一個數(shù)據樣本都包含著轉子系統(tǒng)在某種運行狀態(tài)下12個不同測點位置上的振動特性。顯然,每3個特征來自一個信號通道,任何兩個通道的特征類型都是相同的,但是由于測點不同,即使同一特征也具有不同的特征值,因此這些特征包含著能夠相互補充的故障狀態(tài)信息。

    表4雖然只給出了全部故障數(shù)據的一小部分,但是從表中數(shù)據仍然可以發(fā)現(xiàn),對于相同的故障特征,轉子系統(tǒng)不同測點處的特征值存在明顯差別。這說明轉子系統(tǒng)不同測振點的振動狀況是不一樣的,利用多個傳感器組成多通道的數(shù)據采集陣列可最大程度上獲取轉子系統(tǒng)的完整狀態(tài)信息。

    利用FCM算法和改進FCM算法對轉子故障數(shù)據集分別進行處理,其目標函數(shù)的迭代誤差如圖3所示。

    從圖3可見,當?shù)螖?shù)到達26時,兩種聚類算法目標函數(shù)的迭代誤差都趨于零。但是,很明顯,F(xiàn)CM算法目標函數(shù)的初始迭代誤差較大,而改進的FCM算法目標函數(shù)的迭代誤差在迭代初期就很小,其收斂速度明顯更快,說明MS算法搜索的初始聚類中心已經非常接近故障數(shù)據集的真實聚類中心,從而可以有效加快改進FCM算法的收斂速度。兩種算法對轉子故障數(shù)據集的聚類準確率及運行時間如表5所示;具體的聚類結果如表6所示。

    Li Huaijun, Xie Xiaopeng. Gear fault pattern recognition based on kernel feature fuzzy clustering and fuzzy association entropy[J]. Chinese Journal of Scientific Instrument, 2015, 36(4): 848-855.

    [9] 田再克, 李洪儒, 孫 健, 等. 基于改進MF-DFA和SSM-FCM的液壓泵退化狀態(tài)識別方法[J]. 儀器儀表學報, 2016, 37(8): 1851-1860.

    Tian Zaike, Li Hongru, Sun Jian, et al. Degradation state identification method of hydraulic pump based on improved MF-DFA and SSM-FCM[J]. Chinese Journal of Scientific Instrument, 2016, 37(8): 1851-1860.

    [10] 陳海鵬, 申鉉京, 龍建武, 等. 自動確定聚類個數(shù)的模糊聚類算法[J]. 電子學報, 2017, 45(3):687-694.

    Chen Haipeng, Shen Xuanjing, Long Jianwu, et al. Fuzzy clustering algorithm for automatic identification of clusters[J]. Acta Electronica Sinica, 2017, 45(3): 687-694.

    [11] 董曉峰. 基于RCM分析的智能化汽輪機組故障診斷系統(tǒng)研究[D]. 北京: 華北電力大學, 2012.

    Dong Xiaofeng. Study on intelligent fault diagnosis system for turbogenerator based on RCM analysis[D]. Beijing: North China Electric Power University, 2012.

    [12] 胡雅婷. 可能性聚類方法研究及應用[D]. 長春: 吉林大學, 2012.

    Hu Yating. Research and application of possibilistic clustering method[D]. Changchun: Jilin University, 2012.

    [13] 宋 新, 羅 軍, 王魯平, 等. 基于Mean Shift聚類的邊緣檢測方法[J]. 彈箭與制導學報, 2007, 27(1): 366-368.

    Song Xin, Luo Jun, Wang Luping, et al. Edge detection method based on Mean Shift clustering method[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2007, 27(1): 366-368.

    [14] 康家銀, 紀志成, 龔成龍. 一種核模糊C均值聚類算法及其應用[J]. 儀器儀表學報, 2010, 31(7): 1657-1663.

    Kang Jiayin, Ji Zhicheng, Gong Chenglong. Kernelized fuzzy C-means clustering algorithm and its application[J]. Chinese Journal of Scientific Instrument, 2010, 31(7): 1657-1663.

    [15] Fisher R A. Iris Data Set[DB/OL]. http://archive.ics.uci.edu/ml/datasets/Iris, 1988-07-01.

    [16] Michalski R S. Soybean (Small) Data Set[DB/OL]. http://archive.ics.uci.edu/ml/datasets/Soybean+(Small), 1987-01-01.

    [17] Lei Yaguo, He Zhengjia, Zi Yanyang. A new approach to intelligent fault diagnosis of rotating machinery[J]. Expert Systems with Applications, 2008, 35(4): 1593-1600.

    [18] 張 恒, 趙榮珍. 故障特征選擇與特征信息融合的加權KPCA方法研究[J]. 振動與沖擊, 2014, 33(9): 89-93.

    Zhang Heng, Zhao Rongzhen. Weighted KPCA based on fault feature selection and feature information fusion[J]. Journal of Vibration and Shock, 2014, 33(9): 89-93.

    [19] 程曉涵, 汪愛明, 花如祥, 等. 24種特征指標對軸承狀態(tài)識別的性能研究[J]. 振動、測試與診斷, 2016, 36(2): 351-358.

    Chen Xiaohan, Wang Aiming, Hua Ruxiang, et al. Studying on property of 24 characteristic indexes to bearing state recognition[J]. Journal of Vibration, Measurement & Diagnosis, 2016, 36(2): 351-358.

    Abstract: To solve the problems that initial clustering centers are uncertain and outliers are sensitive in the clustering analysis of rotating machine fault data, an improved fuzzy C-means (FCM) clustering algorithm based on integrated multi-strategy is proposed in this paper. Firstly, the fault data set is divided through the decision attribute of itself, and several equivalence classes derived from the equivalence relation are obtained. Then, each equivalence class is regarded as the feasible domain, and mean shift method is used to search the class center in the feasible domain. Finally, treating the searched class center as the initial clustering center of the clustering algorithm, the similarity between the fault data sample and the corresponding cluster center is measured via the kernel method in the high-dimensional feature space to achieve fuzzy clustering of fault data set. The performance of the proposed method is verified and compared via standard data sets and rotating machinery fault data sets. The results show that the clustering performance of the proposed method is significantly improved compared with the FCM algorithm and the proposed method has significant advantages over the FCM algorithm in terms of convergence speed and clustering accuracy.

    Key words: fault diagnosis; rotating machinery; clustering analysis; fuzzy C-means; clustering performance

    作者簡介: 鄧林峰(1984-),男,副教授。電話: 15193117545; E-mail: denglinfeng2002@163.com

    猜你喜歡
    聚類分析故障診斷
    風力發(fā)電機組齒輪箱軸承故障診斷分析
    基于人工神經網絡的故障診斷專利淺析
    淺談機電一體化設備的故障診斷技術研究
    基于EMD和SSAE的滾動軸承故障診斷方法
    農村居民家庭人均生活消費支出分析
    基于省會城市經濟發(fā)展程度的實證分析
    基于聚類分析的互聯(lián)網廣告投放研究
    “縣級供電企業(yè)生產經營統(tǒng)計一套”表輔助決策模式研究
    基于R8C的汽車OBD通用故障診斷儀設計
    少妇熟女欧美另类| 国产探花在线观看一区二区| 成人高潮视频无遮挡免费网站| 99视频精品全部免费 在线| 男女边摸边吃奶| 99热网站在线观看| 一夜夜www| 天美传媒精品一区二区| 两个人视频免费观看高清| 欧美区成人在线视频| 欧美日韩综合久久久久久| 国产精品久久久久久久电影| 亚洲精品乱码久久久久久按摩| 亚洲精品乱码久久久久久按摩| 免费观看a级毛片全部| 九九在线视频观看精品| 国产精品一区二区三区四区免费观看| 精品不卡国产一区二区三区| 国产av码专区亚洲av| h日本视频在线播放| 欧美日韩一区二区视频在线观看视频在线 | 99热6这里只有精品| videos熟女内射| 国产探花极品一区二区| 国产三级在线视频| 乱人视频在线观看| 少妇熟女欧美另类| 久久久欧美国产精品| 1000部很黄的大片| xxx大片免费视频| 亚洲综合精品二区| 中文字幕人妻熟人妻熟丝袜美| 亚洲色图av天堂| 五月伊人婷婷丁香| 国产av国产精品国产| 婷婷色av中文字幕| 日日摸夜夜添夜夜添av毛片| 身体一侧抽搐| 91aial.com中文字幕在线观看| 日韩制服骚丝袜av| 五月玫瑰六月丁香| 国产精品熟女久久久久浪| 国产 一区 欧美 日韩| 国产黄片美女视频| 成人国产麻豆网| 视频中文字幕在线观看| 免费看av在线观看网站| 欧美日本视频| 亚洲乱码一区二区免费版| 淫秽高清视频在线观看| 免费黄色在线免费观看| a级一级毛片免费在线观看| 1000部很黄的大片| 国产精品日韩av在线免费观看| 99久久精品国产国产毛片| 日日啪夜夜爽| 99视频精品全部免费 在线| 久久久久久伊人网av| 性插视频无遮挡在线免费观看| 亚洲av电影不卡..在线观看| 九草在线视频观看| 男人舔女人下体高潮全视频| 久久韩国三级中文字幕| 亚洲国产精品国产精品| 免费av不卡在线播放| 久久99热这里只有精品18| 美女内射精品一级片tv| 色尼玛亚洲综合影院| 国产精品1区2区在线观看.| 亚洲怡红院男人天堂| 久久午夜福利片| av线在线观看网站| 国产av码专区亚洲av| 亚洲一级一片aⅴ在线观看| 国产在视频线精品| 国产精品麻豆人妻色哟哟久久 | 亚洲精品日本国产第一区| 日日啪夜夜爽| 亚洲欧美日韩东京热| 国产精品人妻久久久久久| 精品熟女少妇av免费看| 免费av观看视频| 国产亚洲精品av在线| 国产老妇女一区| 日韩一本色道免费dvd| 日韩三级伦理在线观看| 亚洲高清免费不卡视频| 午夜福利视频1000在线观看| 男人舔女人下体高潮全视频| 色吧在线观看| 精品久久久久久久久av| 欧美zozozo另类| 久久久久网色| 久久久久久久午夜电影| 精品久久久精品久久久| 99热这里只有是精品在线观看| 国产乱人偷精品视频| 一个人看视频在线观看www免费| 国产一区二区三区综合在线观看 | 嫩草影院新地址| 国产高潮美女av| 国产视频内射| 成人性生交大片免费视频hd| 成年免费大片在线观看| 三级国产精品片| 久久精品国产鲁丝片午夜精品| 日韩欧美一区视频在线观看 | 一级毛片我不卡| 伊人久久国产一区二区| 国产熟女欧美一区二区| 亚洲av.av天堂| 在现免费观看毛片| 欧美变态另类bdsm刘玥| 国产午夜精品久久久久久一区二区三区| 国产精品女同一区二区软件| 女的被弄到高潮叫床怎么办| 亚洲色图av天堂| 最近最新中文字幕免费大全7| 日本欧美国产在线视频| 国产在线男女| 婷婷色综合大香蕉| 午夜老司机福利剧场| 欧美97在线视频| 日韩一本色道免费dvd| 美女xxoo啪啪120秒动态图| 亚洲精品456在线播放app| 国产精品一区二区三区四区免费观看| 久久久久久伊人网av| 在现免费观看毛片| 日日摸夜夜添夜夜添av毛片| 搡女人真爽免费视频火全软件| 插阴视频在线观看视频| 美女被艹到高潮喷水动态| 国产av在哪里看| 国产黄色小视频在线观看| 久久精品国产亚洲av天美| 一级毛片黄色毛片免费观看视频| 国产精品日韩av在线免费观看| 啦啦啦啦在线视频资源| av福利片在线观看| 国产亚洲午夜精品一区二区久久 | 亚洲精品一二三| av免费观看日本| 毛片一级片免费看久久久久| 国产一区有黄有色的免费视频 | 最近手机中文字幕大全| 赤兔流量卡办理| 美女大奶头视频| 嫩草影院入口| 观看免费一级毛片| 国产男女超爽视频在线观看| 婷婷色综合大香蕉| 舔av片在线| 一区二区三区乱码不卡18| 天堂俺去俺来也www色官网 | 熟女电影av网| 黄片wwwwww| or卡值多少钱| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 一级a做视频免费观看| 又爽又黄无遮挡网站| 久久久久久久久久久免费av| 中文精品一卡2卡3卡4更新| 久久精品夜色国产| 少妇高潮的动态图| 床上黄色一级片| 色视频www国产| 欧美人与善性xxx| 国产真实伦视频高清在线观看| 一级毛片久久久久久久久女| 99热6这里只有精品| 成年女人在线观看亚洲视频 | 波野结衣二区三区在线| freevideosex欧美| 国产男人的电影天堂91| 欧美成人a在线观看| 女的被弄到高潮叫床怎么办| 欧美不卡视频在线免费观看| 成人毛片a级毛片在线播放| 亚洲精品第二区| kizo精华| 久久久午夜欧美精品| 久久人人爽人人爽人人片va| 久99久视频精品免费| 午夜免费男女啪啪视频观看| 美女脱内裤让男人舔精品视频| 亚洲美女视频黄频| 日韩强制内射视频| av女优亚洲男人天堂| 免费观看精品视频网站| 国产免费又黄又爽又色| 久久精品国产亚洲av涩爱| 在线观看av片永久免费下载| 国产精品av视频在线免费观看| 亚洲最大成人av| 精品久久久久久成人av| 男女啪啪激烈高潮av片| 久久午夜福利片| 看黄色毛片网站| 在线观看一区二区三区| 乱系列少妇在线播放| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 18禁在线播放成人免费| 国产日韩欧美在线精品| 精品人妻视频免费看| 午夜免费激情av| 日日摸夜夜添夜夜添av毛片| 亚州av有码| 麻豆av噜噜一区二区三区| 国产精品久久久久久久电影| 一夜夜www| 中文字幕久久专区| 美女主播在线视频| 日韩 亚洲 欧美在线| 看十八女毛片水多多多| 久久久久久久午夜电影| 中文字幕人妻熟人妻熟丝袜美| 汤姆久久久久久久影院中文字幕 | 日韩欧美国产在线观看| 91aial.com中文字幕在线观看| 春色校园在线视频观看| 国内精品一区二区在线观看| 国产成人福利小说| av线在线观看网站| 国产午夜精品久久久久久一区二区三区| 最近手机中文字幕大全| 少妇熟女aⅴ在线视频| 精品午夜福利在线看| 欧美成人精品欧美一级黄| 亚洲av日韩在线播放| 成人午夜高清在线视频| 中文字幕av在线有码专区| 成人国产麻豆网| 亚洲精品成人av观看孕妇| 99视频精品全部免费 在线| 18禁在线播放成人免费| 成人一区二区视频在线观看| 国产成人91sexporn| 日韩一区二区三区影片| 国产免费一级a男人的天堂| 日本爱情动作片www.在线观看| 联通29元200g的流量卡| 搡女人真爽免费视频火全软件| 亚洲精品第二区| 国产有黄有色有爽视频| 日本爱情动作片www.在线观看| 久久精品熟女亚洲av麻豆精品 | 三级国产精品片| 欧美xxⅹ黑人| 久久午夜福利片| 亚洲性久久影院| 亚洲自偷自拍三级| 亚洲不卡免费看| 国产av在哪里看| 老司机影院成人| 亚洲经典国产精华液单| 日本免费a在线| 国产淫语在线视频| 久久久成人免费电影| 少妇熟女欧美另类| 国产女主播在线喷水免费视频网站 | 久久久久久久亚洲中文字幕| 少妇丰满av| 久久精品国产亚洲网站| 成人毛片60女人毛片免费| 午夜福利视频精品| 一个人看视频在线观看www免费| 国产精品三级大全| av在线老鸭窝| 99久久精品热视频| 欧美xxxx性猛交bbbb| 欧美激情国产日韩精品一区| 最新中文字幕久久久久| 精品久久久精品久久久| 男女那种视频在线观看| 熟妇人妻久久中文字幕3abv| 亚洲精品久久久久久婷婷小说| 一级a做视频免费观看| 99久久精品国产国产毛片| 亚洲三级黄色毛片| av免费在线看不卡| 亚洲av电影不卡..在线观看| 国产免费福利视频在线观看| 舔av片在线| 亚洲图色成人| 2022亚洲国产成人精品| 国产淫片久久久久久久久| 色网站视频免费| 欧美日韩视频高清一区二区三区二| 美女被艹到高潮喷水动态| 高清日韩中文字幕在线| 熟妇人妻不卡中文字幕| 亚洲精品日韩av片在线观看| 欧美xxxx性猛交bbbb| 寂寞人妻少妇视频99o| 毛片一级片免费看久久久久| 九草在线视频观看| 欧美性猛交╳xxx乱大交人| 成人一区二区视频在线观看| 美女黄网站色视频| 精品人妻一区二区三区麻豆| 久久精品久久久久久噜噜老黄| 99热这里只有是精品在线观看| 成人毛片a级毛片在线播放| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 亚洲国产精品成人综合色| 欧美高清成人免费视频www| 免费不卡的大黄色大毛片视频在线观看 | 免费人成在线观看视频色| 22中文网久久字幕| 日本av手机在线免费观看| 久久久久久久大尺度免费视频| 91狼人影院| 亚洲成人精品中文字幕电影| 久久6这里有精品| 亚洲精品456在线播放app| 国产精品99久久久久久久久| 国产一区二区三区av在线| 最后的刺客免费高清国语| 秋霞伦理黄片| 国产成人精品婷婷| 亚洲三级黄色毛片| 国产国拍精品亚洲av在线观看| 国产色爽女视频免费观看| www.色视频.com| 欧美成人精品欧美一级黄| 国产黄频视频在线观看| 天堂俺去俺来也www色官网 | 国产精品嫩草影院av在线观看| 黄色日韩在线| 国产免费视频播放在线视频 | 精品久久久精品久久久| 国产老妇女一区| 性插视频无遮挡在线免费观看| www.av在线官网国产| 日韩成人av中文字幕在线观看| 中文资源天堂在线| 深爱激情五月婷婷| 尾随美女入室| 国产精品精品国产色婷婷| 国产欧美日韩精品一区二区| 欧美bdsm另类| 亚州av有码| 久久久久性生活片| 欧美日韩综合久久久久久| 听说在线观看完整版免费高清| 嘟嘟电影网在线观看| 插逼视频在线观看| 一级片'在线观看视频| 毛片女人毛片| 赤兔流量卡办理| 国产一区二区在线观看日韩| 在现免费观看毛片| 国内精品美女久久久久久| 一区二区三区乱码不卡18| av天堂中文字幕网| 能在线免费看毛片的网站| 国产高清不卡午夜福利| 国产精品久久视频播放| 国产白丝娇喘喷水9色精品| 搡老妇女老女人老熟妇| 国产精品一区二区三区四区免费观看| 国产成人a∨麻豆精品| 伦理电影大哥的女人| 国精品久久久久久国模美| 国产成人freesex在线| 亚洲精品久久午夜乱码| 少妇高潮的动态图| videossex国产| 夫妻性生交免费视频一级片| 亚洲精品国产av成人精品| 欧美区成人在线视频| 亚洲国产最新在线播放| 亚洲国产高清在线一区二区三| 春色校园在线视频观看| 1000部很黄的大片| 51国产日韩欧美| 午夜日本视频在线| 麻豆av噜噜一区二区三区| 黄片无遮挡物在线观看| 中文字幕亚洲精品专区| 国产伦精品一区二区三区视频9| 亚洲美女搞黄在线观看| 免费电影在线观看免费观看| 亚洲在线观看片| 国产精品人妻久久久影院| 国产成人精品一,二区| 亚洲精品乱久久久久久| 美女被艹到高潮喷水动态| 国国产精品蜜臀av免费| 色5月婷婷丁香| 一个人免费在线观看电影| 亚洲精品一二三| 欧美一区二区亚洲| 色综合色国产| 日韩欧美精品免费久久| 国产成人精品久久久久久| 日本一二三区视频观看| 国产爱豆传媒在线观看| 久久99热6这里只有精品| 欧美日韩在线观看h| 女人被狂操c到高潮| 日韩欧美精品v在线| 欧美 日韩 精品 国产| 欧美bdsm另类| 中文在线观看免费www的网站| 一级毛片 在线播放| 国产成人a区在线观看| 国产精品一区二区三区四区久久| 国产伦精品一区二区三区四那| 成人特级av手机在线观看| 午夜免费男女啪啪视频观看| 天美传媒精品一区二区| 国产精品一区二区在线观看99 | 中国美白少妇内射xxxbb| 搡老妇女老女人老熟妇| 国产成人午夜福利电影在线观看| 亚洲av福利一区| 十八禁国产超污无遮挡网站| 精品久久久久久电影网| 久久精品国产鲁丝片午夜精品| 男女下面进入的视频免费午夜| 干丝袜人妻中文字幕| 免费高清在线观看视频在线观看| 午夜激情欧美在线| 国产乱人视频| 亚洲国产色片| 成年女人看的毛片在线观看| 国产精品久久视频播放| 成人亚洲精品一区在线观看 | 午夜视频国产福利| 有码 亚洲区| 国产成人精品一,二区| 亚洲一级一片aⅴ在线观看| 亚洲av电影不卡..在线观看| av播播在线观看一区| 久久久久性生活片| 久久久久久九九精品二区国产| 五月天丁香电影| 精品不卡国产一区二区三区| 国产黄片视频在线免费观看| 免费看不卡的av| 国产成人91sexporn| 少妇人妻精品综合一区二区| av.在线天堂| 久久99热这里只频精品6学生| 国产日韩欧美在线精品| 国产黄频视频在线观看| 日韩 亚洲 欧美在线| 日韩精品青青久久久久久| 美女高潮的动态| 国内揄拍国产精品人妻在线| 免费少妇av软件| 91久久精品国产一区二区三区| 国产探花极品一区二区| 亚洲综合精品二区| 亚洲精品亚洲一区二区| 晚上一个人看的免费电影| 免费av毛片视频| 国产淫语在线视频| 综合色丁香网| 九九在线视频观看精品| 天天躁夜夜躁狠狠久久av| 在线免费观看不下载黄p国产| 久久99蜜桃精品久久| 高清毛片免费看| 精品午夜福利在线看| 免费av观看视频| 欧美激情在线99| 国产有黄有色有爽视频| 一级黄片播放器| 国产不卡一卡二| 老师上课跳d突然被开到最大视频| 精品久久久久久久久av| 国产一区二区三区综合在线观看 | 日韩强制内射视频| 免费无遮挡裸体视频| 2021少妇久久久久久久久久久| 亚洲精品视频女| 欧美极品一区二区三区四区| 在线免费观看的www视频| 看黄色毛片网站| 免费观看无遮挡的男女| 国产免费福利视频在线观看| 美女高潮的动态| 久久久成人免费电影| 一级毛片aaaaaa免费看小| 一级二级三级毛片免费看| 亚洲经典国产精华液单| 乱系列少妇在线播放| 日本色播在线视频| 老师上课跳d突然被开到最大视频| 国产精品日韩av在线免费观看| 亚洲欧美日韩东京热| 美女cb高潮喷水在线观看| 午夜福利在线在线| 亚洲国产精品国产精品| 久久久精品94久久精品| 内射极品少妇av片p| 国产在视频线精品| 只有这里有精品99| 中文字幕人妻熟人妻熟丝袜美| 久久久久久国产a免费观看| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久com| av在线观看视频网站免费| 国产在线一区二区三区精| 国产一区二区亚洲精品在线观看| 国内精品一区二区在线观看| 久久久久网色| 日韩中字成人| 欧美97在线视频| 亚洲精品久久午夜乱码| 国产成人精品久久久久久| 啦啦啦啦在线视频资源| 亚洲精品乱久久久久久| 麻豆成人av视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品视频女| 亚洲图色成人| 国产视频首页在线观看| 国产精品一区二区在线观看99 | 你懂的网址亚洲精品在线观看| 韩国av在线不卡| 色综合亚洲欧美另类图片| 久久久久久久久中文| 美女内射精品一级片tv| 寂寞人妻少妇视频99o| 欧美激情在线99| 国产精品一区二区三区四区免费观看| 卡戴珊不雅视频在线播放| 国产男女超爽视频在线观看| 18禁裸乳无遮挡免费网站照片| 午夜福利在线观看免费完整高清在| 男人爽女人下面视频在线观看| 亚洲精品亚洲一区二区| 麻豆精品久久久久久蜜桃| 男女边摸边吃奶| 久久精品国产亚洲av天美| 久久精品综合一区二区三区| 中国美白少妇内射xxxbb| 一区二区三区乱码不卡18| 国产av在哪里看| 国产 一区 欧美 日韩| 久久久久久久午夜电影| 久久久久免费精品人妻一区二区| 久久99热这里只频精品6学生| 国产白丝娇喘喷水9色精品| 精品久久久精品久久久| 欧美三级亚洲精品| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| 国产成人aa在线观看| 亚洲精品久久久久久婷婷小说| 久久精品久久久久久久性| 免费观看a级毛片全部| 高清毛片免费看| 在线免费观看的www视频| 久久久久免费精品人妻一区二区| 毛片女人毛片| 少妇被粗大猛烈的视频| 天堂av国产一区二区熟女人妻| 亚洲精品影视一区二区三区av| 久久久午夜欧美精品| 最近最新中文字幕大全电影3| 综合色丁香网| 边亲边吃奶的免费视频| 99热这里只有精品一区| 久久人人爽人人爽人人片va| 女人久久www免费人成看片| 最近的中文字幕免费完整| 久久久久久伊人网av| 伦精品一区二区三区| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 小蜜桃在线观看免费完整版高清| 一级毛片我不卡| 亚洲熟女精品中文字幕| 大又大粗又爽又黄少妇毛片口| 在线观看免费高清a一片| 亚洲精品自拍成人| 丝袜喷水一区| 日日干狠狠操夜夜爽| 日韩一本色道免费dvd| 午夜久久久久精精品| 一夜夜www| 天天躁夜夜躁狠狠久久av| 极品少妇高潮喷水抽搐| 国产黄色小视频在线观看| 中文天堂在线官网| 精品久久久久久久久久久久久| 成人漫画全彩无遮挡| 人妻夜夜爽99麻豆av| 久久精品国产亚洲av天美| 国产精品女同一区二区软件| 亚洲第一区二区三区不卡| 亚洲av不卡在线观看| 国产亚洲最大av| 亚洲第一区二区三区不卡| 成人亚洲精品av一区二区| 91久久精品国产一区二区三区| 免费观看在线日韩| 国产91av在线免费观看| 2018国产大陆天天弄谢| 日日摸夜夜添夜夜添av毛片| 大片免费播放器 马上看| av在线蜜桃| 91精品国产九色| 中文精品一卡2卡3卡4更新| 中文资源天堂在线|