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

    基于旋轉多面體錐的滾動軸承智能故障診斷

    2022-11-14 02:55:33李寶慶楊路航楊宇程軍圣
    湖南大學學報·自然科學版 2022年6期
    關鍵詞:滾動軸承模式識別故障診斷

    李寶慶 楊路航 楊宇 程軍圣

    摘要:針對多面體錐分類器形成的凸包區(qū)域分類邊界數量有限且無法進行不同尺度伸縮的問題,在范數向量化的基礎上增加旋轉因子,提出一種旋轉多面體錐分類器(Rotated and Extended Polyhedral Conic Classifier,REPCC). REPCC增加了凸包區(qū)域的分類邊界數量,并且分類邊界可以在每個維度上自適應伸縮,從而可以更好地擬合正類區(qū)域,提高了分類精度.在兩個不同的滾動軸承數據集上進行實驗驗證,實驗結果表明:REPCC具有更好的分類精度、魯棒性和泛化能力,能夠準確識別滾動軸承的工作狀態(tài)和故障類型,可以用于滾動軸承智能故障診斷.

    關鍵詞:多面體錐分類器;旋轉因子;滾動軸承;模式識別;故障診斷

    中圖分類號:TH 113;TH 165+.3文獻標志碼:A

    Rolling Bearing Intelligent Fault Diagnosis Based on Rotated and Extended Polyhedron Cone

    LI Baoqing,YANG Luhang,YANG Yu,CHENG Junsheng

    (College of Mechanical and Vehicle Engineering,Hunan University,Changsha 410082,China)

    Abstract:In order to solve the problem that the number of classification boundaries of convex hull region formed by polyhedral cone classifier is limited and it is unable to scale in different scales,a rotated and extended polyhedral conic classifier(REPCC)is proposed by adding a rotation factor based on norm vectorization. REPCC increases the number of classification boundaries of the convex hull region,and the classification boundaries can be adaptively scaled in each dimension,which can better fit the positive region and improve the classification accuracy. Experimental verification is carried out on two different rolling bearing datasets. The results show that REPCC has better classification accuracy,robustness and generalization ability,and can accurately identify the working state and fault type of rolling bearing. REPCC can be used for intelligent fault diagnosis of rolling bearing.

    Key words:polyhedral conic classifier;rotation factor;rolling bearing;pattern recognition;fault diagnosis

    滾動軸承是機械設備的關鍵部件,長期在復雜多變載荷環(huán)境下高速運行,工作環(huán)境惡劣,極易發(fā)生故障,直接影響機械系統(tǒng)的總體性能[1].因此,滾動軸承故障診斷研究具有重要意義,已得到許多學者的廣泛關注[2-3].

    滾動軸承智能故障診斷的核心是模式識別,隨著近年來人工智能技術的快速發(fā)展,滾動軸承智能故障診斷方法的研究也迅速發(fā)展[4].人工神經網絡(Artificial Neural Network,ANN)、支持向量機(Support Vector Machine,SVM)等方法被廣泛應用于滾動軸承故障識別,取得了良好的效果[5-6].特別是深度學習(Deep Learning,DL)可以利用深層神經網絡結構直接從原始信號中提取特征,近幾年已在滾動軸承智能故障診斷領域獲得廣泛應用[7-8].但是,這些智能診斷方法也存在一些問題,其中:神經網絡和深度學習需要大量的訓練樣本,而工程實際故障大多為小樣本[9];SVM方法對核的選取不敏感,確定核參數和懲罰系數需要花費較多的時間,另外,由于引入核函數計算,當處理大樣本或者樣本數據本身維度相對較高時,會遇到核函數計算緩慢的問題.

    合正類區(qū)域.本質上,REPCC等基于PCF的分類模型具有與線性SVM相似的性質和復雜性,都是基于統(tǒng)計學習理論的有監(jiān)督機器學習方法,可以使用與SVM相關聯的凸二次規(guī)劃問題進行最大邊距訓練求解.但是,REPCC等分類模型不需要引入核函數,并能夠保證較高的識別率,在參數選擇上比SVM更加便捷.

    當滾動軸承出現故障時,其振動信號的時頻分布與正常情況不同,振動信號的時域和頻域統(tǒng)計參數通常被用作滾動軸承故障診斷的特征[13].本文將REPCC應用于滾動軸承故障診斷,首先提取振動加速度信號的時域和頻域特征作為候選故障特征;然后采用距離評估技術[14](Distance Evaluation Technique,DET)篩選候選故障特征中的敏感特征,DET通過有效因子來評估特征在樣本類別變化時的敏感程度,選擇類間變化明顯和類內變化小的特征,評估因子越大特征越敏感;最后將篩選后的敏感特征作為分類器的輸入特征,輸入到REPCC進行故障識別.

    在兩種不同的滾動軸承數據集上開展實驗驗證,分析REPCC的識別精度、魯棒性和泛化能力.通過對比實驗,與PCC和EPCC相比,REPCC具有更高的識別精度、更好的抗干擾能力、更好的泛化能力.REPCC可以有效識別滾動軸承的工作狀態(tài)和故障類型,在小樣本情況下具有更優(yōu)越的故障識別性能.

    1多面體錐分類器及其擴展型

    1.1多面體錐分類器

    在n維向量空間R中,PCF可以用來分割不相交的兩個點集,對于測試樣本x∈R,分類函數定義如下:

    式(1)所示的函數是一個多面體錐,其頂點位于(a,-γ)∈R×R,PCF的每個水平截面都是凸多邊形.三維向量空間下的PCF示意圖如圖1所示.

    從圖中可以看出,在平面g=0上,多面體錐的四個錐面將平面分為內部區(qū)域和外部區(qū)域,由四條線性邊界包圍構成,這四條邊界包圍的區(qū)域可以看作是一個封閉的二維凸包,此時在凸包內部區(qū)域的點映射到超平面上的值g小于0,外部區(qū)域的點映射到超平面上的值g大于0,實現線性可分.因此,當一個二分類數據集樣本特征維度為2,且其正類分布較為緊湊,負類分布較為分散時,一個二維PCF就能夠很好地對其實現分類擬合.當樣本特征維度大于2 時,與二維類似,PCF會在特征空間形成封閉凸包,并在更高維空間形成決策項來實現分類.

    為直觀表示PCC的分類特點,隨機生成二維數據樣本點,方框點為正類,圓圈點為負類,然后使用PCC進行優(yōu)化分類,分類結果如圖2所示.從圖中可以看出,分類邊界完全線性,實現了對隨機樣本點的正確分類.

    1.2擴展型多面體錐分類器

    EPCC中的ξ是一個屬于R空間的向量,其幾何意義為:分類截面所形成的凸包頂點在n個維度上的不同伸縮因子.因此,PCC中的ξ為常數,只能實現分類截面整體尺度的縮放,而EPCC的分類邊界則更加靈活,可以調整凸包邊界實現對正類樣本的更好擬合.當ξ向量的每個值都相等時,EPCC就退化成了PCC.

    對于測試樣本x,g(x)<0屬于正類,g(x)>0屬于負類.對PCC進行最大邊界訓練,對于正類樣本強制g(x)<-1,負類樣本g(x)>+1.為了實現這種方式的訓練,對樣本特征向量和目標函數進行向量維度的擴展,對于PCC而言,將特征向量擴展為如下形式:

    權重向量擴展為:

    由于權重變化為負的,所以對于正類樣本f>0,負類樣本f<0.

    同樣地,對于EPCC而言,將特征向量擴展為如下形式:

    權重向量擴展為:

    2旋轉多面體錐分類器

    2.1旋轉多面體錐函數

    雖然EPCC可以使得分類決策的凸包邊界能夠在每個維度上實現不同權重的伸縮,改善了原始分類區(qū)域擬合正類樣本的精度,但是EPCC的分類邊界數量始終等于向量維度的2倍,以二維特征向量分類為例,形成的分類邊界是四邊形,分類邊界只有4 條,有限的分類邊界仍然無法保證對所有正類樣本進行充分地擬合,來實現更高的分類精度.

    為了擴大分類截面多邊形的邊數,從而使得分類截面能夠更好地擬合正類區(qū)域,將坐標軸進行旋轉,然后融合重新繪制在原始的坐標軸中,形成新的凸包區(qū)域,得到旋轉多面體錐函數(Rotated Polyhedral Conic Function,RPCF).三維向量空間下,RPCF截面區(qū)域的旋轉融合過程如圖3所示.從示意圖中可以看出,原始錐面形成的凸包邊界只有4條,融合后的凸包邊界增加到了8條,并且對于新的凸包區(qū)域,凸包上的8個頂點同時能夠在四個方向維度上進行不同比例的伸縮,從而更好地擬合正類區(qū)域.

    2.2旋轉因子計算

    基于RPCF定義REPCC,公式如式(9)所示.式中,A是維數為n的旋轉變換矩陣,即旋轉因子,μ是一個權重向量,ξ和μ控制四個維度方向的伸縮尺度.

    在高維空間進行坐標旋轉須保證旋轉因子A是正交矩陣[15].在二維笛卡爾坐標系中,旋轉因子A的形式如式(10)所示.坐標繞原點的旋轉可以通過原始坐標向量和坐標變換矩陣相乘得到,如式(11)所示.其中,0≤θ≤π,x表示在原坐標系中x逆時針旋轉θ弧度得到的新坐標.

    在n維空間中,在兩個特征維度的坐標系上的旋轉因子A可以采用一個n維的單位方陣表示,如式(12)所示.

    由于是在n維空間進行旋轉,需要兩兩組成子坐標空間進行旋轉,共有n(n-1)/2組子空間,因此需要將每個子空間的旋轉矩陣進行左乘.旋轉因子的計算流程如表1所示.

    選擇一個初始的旋轉角度0,根據表1計算相對應的旋轉因子A.確定旋轉因子后,與PCC和EPCC的變換策略相同,REPCC將特征向量擴展為式(13)所示的形式,將權重向量擴展為式(14)所示的形式.

    REPCC變換后的特征向量和權重向量的維度是原始的3倍,屬于R空間.分類截面的邊界數量是原始邊界數量的2倍,并且繼承了EPCC的優(yōu)點,可以實現分類截面在2n個維度上的不同尺度伸縮,可以擬合出更為合理的正類區(qū)域,從而提高分類精度.

    為直觀表示,根據式(1)給出一個二維特征的分類函數,具體參數如式16所示,該仿真信號的三維空間示意圖如圖4所示.根據REPCC,給定一個初始旋轉角度后計算旋轉因子,其三維空間示意圖如圖5所示.

    從圖4中可以看出,PCF在平面g=0形成的二維凸包區(qū)域存在4條線性邊界,分別為:0.44x+0.44x-2.76≤0,0.22x-0.22x-0.12≤0,0.44x-0.22x-1.44≤0,-0.22x+0.44x-1.44≤0.

    從圖5中可以看出,經過旋轉融合后的多面體錐函數在平面g=0形成的凸包區(qū)域存在8條線性邊界,RPCF能夠以更加合理的尺度形成新的凸包邊界來擬合正類樣本.

    2.3多分類模型

    根據決策函數,PCC、EPCC和REPCC三種分類器的±1邊距都轉化為類似SVM方法的±1邊距,采用凸二次規(guī)劃算法進行最大邊距訓練,可以得到支持向量,然后根據決策函數逆變換獲得權重向量.目標函數如式(17)所示.

    式中:y為樣本的類別標簽,對于正類y=1,負類

    采用一對多的策略將二分類問題推廣到多分類模型,對于m類問題,以當前類為正類樣本,其他類為負類樣本,構造m個分類器進行訓練[16].當進行新樣本的預測時,將樣本輸入到每個分類器中,以輸出值最大的那個分類器對應的正類類別作為新樣本的最終決策.決策函數如式(19)所示.

    3基于REPCC的滾動軸承故障診斷

    3.1基于REPCC的故障診斷步驟

    選取滾動軸承振動信號的時域和頻域統(tǒng)計參數作為故障診斷的特征.時域統(tǒng)計參數選取16個,具體有:均值、方根幅值、均方根值、峰值、標準差、偏斜度、峭度、峰度指標、裕度指標、波形指標、脈沖指標、方差、最大值、最小值、峰峰值、絕對平均值.頻域統(tǒng)計參數選取8個,具體有:頻率均值、頻率均方根值、頻率斜度、頻率峭度、頻率方差、頻率標準差、中心頻率、重心頻率.時域和頻域統(tǒng)計參數共計24個.

    不是所有的特征都對區(qū)分不同的滾動軸承故障敏感,因此在進行分類訓練時,需要從候選特征中選擇對故障敏感的顯著特征.DET方法通過有效因子來評估特征在樣本類別變化時的敏感程度,可以有效提取候選特征中的敏感特征.

    基于REPCC的滾動軸承故障診斷的主要步驟如下:

    1)采集滾動軸承不同故障狀態(tài)的振動加速度信號,并將每種故障狀態(tài)的數據分為訓練樣本和測試樣本,訓練樣本和測試樣本各占50%;

    2)提取振動加速度信號的時域和頻域統(tǒng)計參數,共計24個特征參數,進行歸一化處理后構建候選特征集;

    3)采用DET方法從候選特征集中篩選敏感特征,評估因子的閾值設置為全部特征評估因子的均值,選擇大于評估因子閾值的敏感特征作為故障特征,構建新的特征集;

    4)將經過DET篩選后的特征集輸入REPCC中進行訓練和測試,進行滾動軸承故障診斷.

    3.2實驗數據集

    為驗證基于REPCC的滾動軸承故障診斷方法的有效性和優(yōu)越性,采用數據集A和B兩個實驗數據集進行分析驗證,具體如表2所示.

    數據集A的數據來源為美國凱斯西儲大學軸承數據中心,實驗軸承為型號6205-2RS的SKF深溝球軸承,采用電火花加工設置單點損傷,本文采用的數據是驅動端軸承的振動加速度測試數據,采樣頻率為48kHz,電機轉速為1 797 r/min、1 772 r/min、1 750 r/min、1 730 r/min,以每12 000個數據點劃分為一個樣本,總共400個樣本.

    數據集B的數據來源于旋轉機械故障實驗臺的實際測試數據,實驗滾動軸承型號為6205-2RS- SKF,實驗臺如圖6所示.所有的故障均采用激光切割加工,故障深度為0.2 mm,電機轉速為986.1 r/min,采樣頻率10 kHz,以6 000個數據點作為一個樣本,總共400個樣本.

    3.3特征集構建

    為更加客觀地體現REPCC的分類性能,不對采集得到的振動加速度信號進行處理,直接提取其時域和頻域統(tǒng)計參數作為候選特征集.然后,采用DET方法從候選特征集中篩選敏感特征構建新的特征集.在進行敏感特征篩選時,選擇大于評估因子閾值的敏感特征作為故障特征,評估因子的閾值設置為全部特征評估因子的均值,有效因子的計算方法參照文獻[14].

    數據集A和數據集B經過時域和頻域特征提取并進行歸一化后,各個特征的評估因子及其均值如圖7~圖8所示.

    在圖7-8中,黑色虛線為平均值,是評估因子閾值,選擇評估因子位于黑色虛線上方的特征構建新的故障特征集.從兩圖中可以看出,篩選后重新構建的特征集的維數下降,數據集A的特征維數從24下降為14,數據集B的特征維數從24下降為11.由于數據集A和數據集B的故障尺寸設置以及電機轉速等試驗工況不同,因此經DET篩選后構建的故障特征集不同.

    3.4實驗結果及分析

    將經過DET篩選構建的新特征集輸入分類器進行滾動軸承智能故障診斷.為驗證REPCC的有效性和優(yōu)越性,在表2所示的兩個數據集上采用SVM、PCC、EPCC和REPCC進行分類識別,考察不同分類器的識別精度、魯棒性和泛化能力.

    3.4.1分類器的識別精度

    在使用不同的分類器時,需要確定分類器所需的相關參數.SVM需要確定的參數為高斯核參數δ和懲罰系數C采用網格搜索法獲取最優(yōu)參數.PCC和EPCC無須進行核計算,需要確定的參數為懲罰系數C.由于引入了旋轉因子A,REPCC需要確定的參數為懲罰系數C和旋轉弧度θ,采用控制變量法,旋轉弧度0在EPCC的懲罰系數C基礎上進行優(yōu)化選擇.

    在數據集A和數據集B上,對DET篩選后的數據特征集進行五折交叉驗證實驗,獨立重復30次,統(tǒng)計測試數據30次的平均準確率.PCC和EPCC的識別精度隨懲罰系數C的變化如圖9所示.對于REPCC,在EPCC懲罰系數C的基礎上,REPCC的識別精度隨旋轉弧度e的變化如圖10所示.利用圖9 和圖10選擇最優(yōu)化參數,各分類器確定的最優(yōu)化參數如表3所示.

    為分析REPCC的識別精度,在表3所示的最優(yōu)化參數下,隨機獨立重復30次實驗,得到平均識別準確率如表4所示.從表4中可以看出,在兩個滾動軸承實驗數據集上,REPCC的識別準確率最高,優(yōu)于改進前的PCC和EPCC的識別準確率,也優(yōu)于SVM的識別準確率.實驗結果表明,REPCC具有更高的識別精度.

    從表2中可以知道,與數據集A相比,數據集B的試驗工況更加單一.從表4中還可以看出,各分類器在數據集B上的分類精度都高于數據集A,符合實際情況.

    3.4.2分類器的魯棒性

    為分析REPCC的魯棒性,在兩個滾動軸承實驗數據集上進行添加離群點實驗,從其他類的訓練樣本中隨機選擇m個樣本作為離群點,添加到當前類的訓練樣本集中再進行訓練,然后驗證測試樣本的識別準確率.對于不同的m值,每次都隨機獨立重復30次,然后統(tǒng)計測試樣本的平均識別精度.實驗結果如圖11和圖12所示.

    從圖11和圖12中可以看出,在兩個滾動軸承實驗數據集上,隨著訓練集中混雜的離群點個數增加,測試集的識別準確率都出現下降,SVM的下降程度最嚴重.圖11和圖12對比分析可以看出,對于不同的數據集,PCC和EPCC的下降趨勢變化較大,REPCC的下降趨勢基本一致.整體來看,REPCC的下降趨勢最低,表明REPCC具有更好的魯棒性.

    3.4.3分類器的泛化能力

    為分析REPCC的泛化能力,在兩個滾動軸承實驗數據集上,以不同的數量樣本作為訓練集,分析測試樣本集的識別準確率隨訓練集樣本數量的變化情況.對于兩個實驗數據集,每種故障狀態(tài)的數據分為訓練樣本和測試樣本,隨機選取其中一半為訓練樣本,其余的為測試樣本.在訓練樣本中根據不同比例構建不同數量樣本的訓練集.在最優(yōu)參數下,不同樣本數量的訓練集每次隨機獨立重復30次,然后統(tǒng)計測試樣本的平均識別精度.實驗結果如圖13和圖14所示.

    從圖13-14中可以看出,兩個數據集的趨勢基本一致,隨著訓練樣本數量的增加,各分類器的識別準確率也增加;REPCC的識別準確率最高,而且在訓練樣本較少的情況下REPCC的優(yōu)勢更加明顯,識別準確率大于其他分類器.REPCC不僅提高了PCC和EPCC的分類精度,而且具有更好的泛化能力.

    4結論

    本文在EPCC的基礎上引入旋轉因子,提出了REPCC模式識別方法,并將其用于滾動軸承智能故障診斷.REPCC增加了分類邊界數量,并且分類邊界可以在每個維度上自適應伸縮,改進了分類精度,并且具有更好的抗干擾性能和泛化能力.在兩個滾動軸承數據集上進行實驗驗證,結果表明REPCC可以準確識別滾動軸承的工作狀態(tài)和故障類型,REPCC可以用于滾動軸承智能故障診斷.

    REPCC是一種二分類約束算法,但是其分類模型也可以作為一種單分類模型,在后續(xù)的研究中可以尋找一種在只有單一樣本分布下的模型優(yōu)化求解算法,實現單分類.

    參考文獻

    [1]王玉靜,呂海巖,康守強,等.不同型號滾動軸承故障診斷方法[J].中國電機工程學報,2021,41(1):267-276.

    WANG Y J,LYU H Y,KANG S Q,et al. Fault diagnosis method for different types of rolling bearings [J]. Proceedings of the CSEE,2021,41(1):267-276(In Chinese)

    [2]程軍圣,馬利,潘海洋,等.基于EEMD和改進VPMCD的滾動軸承故障診斷方法[J].湖南大學學報(自然科學版),2014,41(10):22-26.

    CHENG J S,MA L,PAN H Y,et al. A fault diagnosis method for rolling bearing based on EEMD and improved VPMCD [J]. Journal of Hunan University (Natural Sciences),2014,41(10):22-26.(In Chinese)

    [3] QIAO M Y,TANG X X,LIU Y X,et al. Fault diagnosis method of rolling bearings based on VMD and MDSVM [J]. Multimedia Tools and Applications,2021,80(10):14521-14544.

    [4]ALSHORMAN O,IRFAN M,SAAD N,et al. A review of artificial intelligence methods for condition monitoring and fault diagnosis of rolling element bearings for induction motor [J]. Shock and Vibration,2020,2020:8843759.

    [5] LI B,CHOW M Y,TIPSUWAN Y,et al. Neural-network-based motor rolling bearing fault diagnosis [J]. IEEE Transactions on Industrial Electronics,2000,47(5):1060-1069.

    [6]YANG Y,YU D J,CHENG J S. A fault diagnosis approach for roller bearing based on IMF envelope spectrum and SVM[J].Measurement,2007,40(9/10):943-950.

    [7]LI X,ZHANG W,DING Q. Understanding and improving deep learning-based rolling bearing fault diagnosis with attention mechanism[J].Signal Processing,2019,161:136-154.

    [8]JIA F,LEI Y G,LIN J,et al. Deep neural networks:a promising tool for fault characteristic mining and intelligent diagnosis of rotating machinery with massive data [J]. Mechanical Systems and Signal Processing,2016,72/73:303-315.

    [9]SHAO H D,JIANG H K,ZHANG X,et al. Rolling bearing fault diagnosis using an optimization deep belief network [J]. Measurement Science and Technology,2015,26(11):115002.

    [10] OZTURK G,TAHIR CIFTCI M.Clustering based polyhedralconic functions algorithm in classification[J].Journal of Industrial & Management Optimization,2015,11(3):921-932.

    [11] CEVIKALP H,YAVUZ H S.Fast and accurate face recognition with image sets[C]//2017 IEEE International Conference on Computer Vision Workshops.October 22-29,2017,Venice,Italy.IEEE,2017:1564-1572.

    [12] CEVIKALP H,TRIGGS B.Polyhedral conic classifiers for visual object detection and classification [C]//2017 IEEE Conference on Computer Vision and Pattern Recognition.July 21-26,2017,Honolulu,HI,USA.IEEE,2017:4114-4122.

    [13]程軍圣,鄭近德,楊宇,等.基于部分集成局部特征尺度分解與拉普拉斯分值的滾動軸承故障診斷模型[J].振動工程學報,2014,27(6):942-950.

    CHENG J S,ZHENG J D,YANG Y,et al.Fault diagnosis model for rolling bearing based on partly ensemble local characteristicscale decomposition and Laplacian score[J].Journal of Vibration Engineering,2014,27(6):942-950.(In Chinese)

    [14] GEORGOULAS G,LOUTAS T,STYLIOS C D,et al.Bearing fault detection based on hybrid ensemble detector and empirical mode decomposition[J].Mechanical Systems and Signal Processing,2013,41(1/2):510-525.

    [15] REDI M H. Oceanic isopycnal mixing by coordinate rotation [J]. Journal of Physical Oceanography,1982,12(10):1154-1158.

    [16] VAPNIK V N.An overview of statistical learning theory[J].IEEE Transactions on Neural Networks,1999,10(5):988-999.

    猜你喜歡
    滾動軸承模式識別故障診斷
    淺談模式識別在圖像識別中的應用
    電子測試(2017年23期)2017-04-04 05:06:50
    第四屆亞洲模式識別會議
    滾動軸承裝配方法及注意事項
    滾動軸承安裝使用中的幾個問題
    基于小波包絡分析的滾動軸承典型故障診斷技術研究
    商情(2016年11期)2016-04-15 07:39:54
    因果圖定性分析法及其在故障診斷中的應用
    基于包絡分析法的軸承故障分析
    第3屆亞洲模式識別會議
    基于LCD和排列熵的滾動軸承故障診斷
    基于WPD-HHT的滾動軸承故障診斷
    機械與電子(2014年1期)2014-02-28 02:07:31
    制服人妻中文乱码| 99久久99久久久精品蜜桃| 日本一本二区三区精品| 精品人妻一区二区三区麻豆 | 国产野战对白在线观看| 神马国产精品三级电影在线观看| 每晚都被弄得嗷嗷叫到高潮| 精品人妻偷拍中文字幕| 国产主播在线观看一区二区| 最近最新免费中文字幕在线| 法律面前人人平等表现在哪些方面| 国产伦一二天堂av在线观看| 久久久久亚洲av毛片大全| 亚洲av免费在线观看| 一边摸一边抽搐一进一小说| 午夜久久久久精精品| 欧美不卡视频在线免费观看| 国产一区在线观看成人免费| 免费人成视频x8x8入口观看| 色尼玛亚洲综合影院| 激情在线观看视频在线高清| 国产高清激情床上av| 欧美一区二区国产精品久久精品| 性色av乱码一区二区三区2| 免费大片18禁| 制服丝袜大香蕉在线| 69人妻影院| 舔av片在线| 精品免费久久久久久久清纯| 中出人妻视频一区二区| 老司机深夜福利视频在线观看| 神马国产精品三级电影在线观看| 欧美性感艳星| 午夜激情福利司机影院| 国产熟女xx| 97碰自拍视频| 哪里可以看免费的av片| 三级毛片av免费| 精品久久久久久久人妻蜜臀av| 1000部很黄的大片| 欧美性猛交黑人性爽| 欧美日韩福利视频一区二区| 欧美高清成人免费视频www| 日韩欧美 国产精品| 岛国视频午夜一区免费看| 99riav亚洲国产免费| www.熟女人妻精品国产| 中文亚洲av片在线观看爽| 国产精品久久久人人做人人爽| 日日干狠狠操夜夜爽| 女人被狂操c到高潮| 亚洲一区二区三区不卡视频| 村上凉子中文字幕在线| 露出奶头的视频| 亚洲av日韩精品久久久久久密| 国产乱人伦免费视频| 美女被艹到高潮喷水动态| 国产高清三级在线| 国产又黄又爽又无遮挡在线| 成人一区二区视频在线观看| 国产免费男女视频| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 99热6这里只有精品| 一级毛片高清免费大全| 欧美日本亚洲视频在线播放| 欧美成人a在线观看| 欧洲精品卡2卡3卡4卡5卡区| 丁香六月欧美| 国产探花极品一区二区| 亚洲av五月六月丁香网| 欧美三级亚洲精品| 伊人久久大香线蕉亚洲五| 国产午夜精品久久久久久一区二区三区 | 天天一区二区日本电影三级| 精品国产美女av久久久久小说| 18美女黄网站色大片免费观看| 久久久久精品国产欧美久久久| 小说图片视频综合网站| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 国产亚洲精品av在线| 亚洲真实伦在线观看| 超碰av人人做人人爽久久 | 日本黄大片高清| 亚洲18禁久久av| 久久久国产成人免费| 99精品在免费线老司机午夜| 午夜影院日韩av| 男女之事视频高清在线观看| 久久久久精品国产欧美久久久| 日韩中文字幕欧美一区二区| 国产成年人精品一区二区| 99在线视频只有这里精品首页| 日日夜夜操网爽| 又粗又爽又猛毛片免费看| 免费电影在线观看免费观看| 最近最新中文字幕大全免费视频| 女警被强在线播放| 精品久久久久久成人av| 韩国av一区二区三区四区| 国产极品精品免费视频能看的| 美女黄网站色视频| 成人无遮挡网站| 午夜福利成人在线免费观看| 国产探花极品一区二区| 亚洲自拍偷在线| www.999成人在线观看| av视频在线观看入口| 美女 人体艺术 gogo| 亚洲精品乱码久久久v下载方式 | 一二三四社区在线视频社区8| 欧美成人性av电影在线观看| 给我免费播放毛片高清在线观看| 亚洲av电影不卡..在线观看| 人妻夜夜爽99麻豆av| 欧美性猛交黑人性爽| 欧美3d第一页| 欧美最新免费一区二区三区 | 亚洲性夜色夜夜综合| 蜜桃亚洲精品一区二区三区| 又黄又粗又硬又大视频| 日韩国内少妇激情av| 精品99又大又爽又粗少妇毛片 | av国产免费在线观看| 中文字幕人妻丝袜一区二区| or卡值多少钱| 日韩欧美精品免费久久 | 久久精品国产99精品国产亚洲性色| 中文在线观看免费www的网站| 欧美日韩黄片免| 中文字幕熟女人妻在线| 黄色丝袜av网址大全| 日本黄大片高清| 淫秽高清视频在线观看| 久久久久精品国产欧美久久久| 三级国产精品欧美在线观看| 18禁美女被吸乳视频| 久久久久精品国产欧美久久久| 国产免费一级a男人的天堂| 国产精品爽爽va在线观看网站| 亚洲精品456在线播放app | 国产精品99久久久久久久久| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 在线天堂最新版资源| 长腿黑丝高跟| 国产单亲对白刺激| 色哟哟哟哟哟哟| av在线天堂中文字幕| 欧美另类亚洲清纯唯美| 精品电影一区二区在线| 久久天躁狠狠躁夜夜2o2o| 亚洲在线自拍视频| 日韩免费av在线播放| 18禁国产床啪视频网站| 可以在线观看毛片的网站| 国内精品一区二区在线观看| 日本五十路高清| 成年女人永久免费观看视频| 成年人黄色毛片网站| 可以在线观看毛片的网站| 国产精品综合久久久久久久免费| 国产免费av片在线观看野外av| 日韩欧美免费精品| 午夜精品久久久久久毛片777| 国产精品久久电影中文字幕| 亚洲午夜理论影院| 夜夜看夜夜爽夜夜摸| 搡女人真爽免费视频火全软件 | 国产高清videossex| 在线观看av片永久免费下载| 一夜夜www| 日韩 欧美 亚洲 中文字幕| 日韩欧美精品免费久久 | 女人高潮潮喷娇喘18禁视频| 青草久久国产| 久久精品亚洲精品国产色婷小说| 高清日韩中文字幕在线| 最新在线观看一区二区三区| 国产一级毛片七仙女欲春2| 极品教师在线免费播放| 高潮久久久久久久久久久不卡| 女人高潮潮喷娇喘18禁视频| 日韩亚洲欧美综合| 韩国av一区二区三区四区| 国产精品久久久人人做人人爽| 国产精品美女特级片免费视频播放器| 日韩精品中文字幕看吧| 女生性感内裤真人,穿戴方法视频| 精品午夜福利视频在线观看一区| 国产午夜福利久久久久久| 叶爱在线成人免费视频播放| 啦啦啦韩国在线观看视频| 无遮挡黄片免费观看| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 18禁黄网站禁片免费观看直播| 18美女黄网站色大片免费观看| 丁香六月欧美| 黄片小视频在线播放| 无人区码免费观看不卡| 欧美区成人在线视频| 欧美性猛交黑人性爽| 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 国模一区二区三区四区视频| 日韩国内少妇激情av| 亚洲国产欧洲综合997久久,| 国产久久久一区二区三区| 中文字幕av成人在线电影| 亚洲国产精品成人综合色| 夜夜躁狠狠躁天天躁| 亚洲国产精品sss在线观看| 禁无遮挡网站| 天天一区二区日本电影三级| 久久香蕉国产精品| 国产高清videossex| 99国产极品粉嫩在线观看| 久久精品91无色码中文字幕| 欧美精品啪啪一区二区三区| 极品教师在线免费播放| xxx96com| 国产伦精品一区二区三区视频9 | 又粗又爽又猛毛片免费看| 色精品久久人妻99蜜桃| 欧美+日韩+精品| 亚洲激情在线av| 日韩欧美精品免费久久 | 国产野战对白在线观看| 国产亚洲欧美在线一区二区| 俺也久久电影网| 校园春色视频在线观看| 亚洲午夜理论影院| 久久久久久人人人人人| 国产国拍精品亚洲av在线观看 | 男女那种视频在线观看| 999久久久精品免费观看国产| 亚洲性夜色夜夜综合| 亚洲人成网站高清观看| 久久99热这里只有精品18| x7x7x7水蜜桃| 欧美日韩黄片免| 久久久久久久精品吃奶| 亚洲成av人片在线播放无| 91麻豆精品激情在线观看国产| 欧美黄色淫秽网站| 久久久久亚洲av毛片大全| 国产亚洲av嫩草精品影院| 老司机深夜福利视频在线观看| 哪里可以看免费的av片| 男女那种视频在线观看| 国内揄拍国产精品人妻在线| 在线观看日韩欧美| 欧美高清成人免费视频www| 国产三级在线视频| 欧美日韩乱码在线| 亚洲国产精品久久男人天堂| 午夜两性在线视频| 黄色女人牲交| 免费看十八禁软件| 真实男女啪啪啪动态图| 亚洲中文字幕一区二区三区有码在线看| 18禁国产床啪视频网站| 禁无遮挡网站| 美女被艹到高潮喷水动态| 亚洲美女视频黄频| 色综合欧美亚洲国产小说| 欧美一区二区精品小视频在线| 亚洲精华国产精华精| 在线播放无遮挡| 99热只有精品国产| 午夜激情福利司机影院| 黑人欧美特级aaaaaa片| 性欧美人与动物交配| 成人永久免费在线观看视频| 人人妻人人澡欧美一区二区| 天天躁日日操中文字幕| 亚洲国产精品久久男人天堂| 一级黄色大片毛片| 他把我摸到了高潮在线观看| 母亲3免费完整高清在线观看| 99热这里只有精品一区| 麻豆成人av在线观看| 亚洲av成人av| 97超级碰碰碰精品色视频在线观看| 999久久久精品免费观看国产| 老鸭窝网址在线观看| 成人永久免费在线观看视频| 久久中文看片网| 三级毛片av免费| 亚洲va日本ⅴa欧美va伊人久久| 热99在线观看视频| 国产欧美日韩一区二区精品| 精品久久久久久久久久久久久| 全区人妻精品视频| 丰满乱子伦码专区| 亚洲精品乱码久久久v下载方式 | 嫩草影院精品99| 色视频www国产| 色老头精品视频在线观看| 午夜福利欧美成人| 色播亚洲综合网| 欧美大码av| 我的老师免费观看完整版| 少妇的逼好多水| 国产成人欧美在线观看| 精品人妻1区二区| 日韩有码中文字幕| 精品电影一区二区在线| 一本一本综合久久| 日本黄大片高清| 午夜老司机福利剧场| www日本在线高清视频| av在线天堂中文字幕| 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| 国产av麻豆久久久久久久| 精品乱码久久久久久99久播| 两个人看的免费小视频| 深夜精品福利| 亚洲av中文字字幕乱码综合| 成年女人永久免费观看视频| 国产激情偷乱视频一区二区| 亚洲av日韩精品久久久久久密| 免费观看精品视频网站| 亚洲内射少妇av| 成人18禁在线播放| 日韩免费av在线播放| 亚洲成人免费电影在线观看| 久久欧美精品欧美久久欧美| 51午夜福利影视在线观看| 国产高清三级在线| 午夜福利18| 成人午夜高清在线视频| 日日摸夜夜添夜夜添小说| 最新中文字幕久久久久| 亚洲人成电影免费在线| 成人高潮视频无遮挡免费网站| 美女黄网站色视频| 亚洲一区高清亚洲精品| 午夜两性在线视频| 免费观看人在逋| 日本一本二区三区精品| 亚洲成人中文字幕在线播放| 国产精品嫩草影院av在线观看 | 欧美日韩中文字幕国产精品一区二区三区| 色av中文字幕| 亚洲 国产 在线| 欧美成人一区二区免费高清观看| 日韩欧美在线乱码| 丰满人妻一区二区三区视频av | www日本在线高清视频| www.999成人在线观看| 亚洲人成网站在线播| 亚洲精品亚洲一区二区| 69av精品久久久久久| 成人精品一区二区免费| 国产高清激情床上av| 一级作爱视频免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品自产拍在线观看55亚洲| 国产aⅴ精品一区二区三区波| 午夜影院日韩av| 国产精品 国内视频| 国产激情偷乱视频一区二区| 日韩欧美三级三区| 国产欧美日韩一区二区精品| 夜夜躁狠狠躁天天躁| 成人av一区二区三区在线看| АⅤ资源中文在线天堂| 真人一进一出gif抽搐免费| 亚洲欧美日韩东京热| 久久精品影院6| 欧美日韩亚洲国产一区二区在线观看| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 中文字幕精品亚洲无线码一区| 中文资源天堂在线| 91麻豆精品激情在线观看国产| 一边摸一边抽搐一进一小说| 午夜视频国产福利| 搡女人真爽免费视频火全软件 | 十八禁网站免费在线| xxx96com| 丰满的人妻完整版| 女人被狂操c到高潮| 国产av在哪里看| e午夜精品久久久久久久| 久久久久亚洲av毛片大全| 97超级碰碰碰精品色视频在线观看| 国产野战对白在线观看| 一卡2卡三卡四卡精品乱码亚洲| 精品国产三级普通话版| 欧美大码av| 欧美黑人欧美精品刺激| 国产伦精品一区二区三区四那| 国产精品日韩av在线免费观看| 人人妻人人看人人澡| 最近视频中文字幕2019在线8| 国产亚洲精品一区二区www| 免费大片18禁| 国产精品一区二区三区四区免费观看 | 免费在线观看日本一区| 久久久久久国产a免费观看| 国产欧美日韩精品一区二区| 天天添夜夜摸| 又黄又爽又免费观看的视频| 最近最新免费中文字幕在线| 不卡一级毛片| 国产精品美女特级片免费视频播放器| 免费看a级黄色片| www国产在线视频色| 国产精品香港三级国产av潘金莲| 最新美女视频免费是黄的| 亚洲熟妇中文字幕五十中出| 欧美高清成人免费视频www| 免费大片18禁| 亚洲精品一卡2卡三卡4卡5卡| 俺也久久电影网| 我的老师免费观看完整版| 国产成人av教育| 亚洲精品一区av在线观看| 亚洲精品456在线播放app | 国产av在哪里看| 免费看日本二区| 在线观看一区二区三区| 观看免费一级毛片| 欧美一区二区国产精品久久精品| 香蕉久久夜色| 91久久精品国产一区二区成人 | 欧美在线一区亚洲| 手机成人av网站| 亚洲国产精品999在线| 国产成人av教育| 成人特级黄色片久久久久久久| 香蕉丝袜av| 日韩欧美在线乱码| 国产精品一区二区三区四区久久| 婷婷丁香在线五月| av在线蜜桃| 久久久久国内视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品日韩av在线免费观看| 黑人欧美特级aaaaaa片| 欧美bdsm另类| 国产亚洲欧美98| 男人舔奶头视频| 51国产日韩欧美| 亚洲乱码一区二区免费版| 校园春色视频在线观看| 美女黄网站色视频| 国产成人a区在线观看| www.色视频.com| 亚洲国产精品久久男人天堂| 午夜a级毛片| 亚洲欧美精品综合久久99| 欧美av亚洲av综合av国产av| 在线观看免费午夜福利视频| 色哟哟哟哟哟哟| 别揉我奶头~嗯~啊~动态视频| 国产免费一级a男人的天堂| 成年人黄色毛片网站| 亚洲中文字幕日韩| 一本久久中文字幕| 黄色成人免费大全| 色精品久久人妻99蜜桃| av福利片在线观看| av黄色大香蕉| 九色成人免费人妻av| 国产午夜精品论理片| 国产精品久久久久久人妻精品电影| 全区人妻精品视频| 免费在线观看成人毛片| 欧美bdsm另类| 中文字幕久久专区| 深爱激情五月婷婷| av福利片在线观看| 18禁黄网站禁片午夜丰满| 国产成人系列免费观看| 国产综合懂色| 天天添夜夜摸| 变态另类成人亚洲欧美熟女| 18禁黄网站禁片午夜丰满| 18禁黄网站禁片免费观看直播| 中文在线观看免费www的网站| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 九九久久精品国产亚洲av麻豆| 亚洲欧美精品综合久久99| 亚洲美女黄片视频| 欧美黄色淫秽网站| 亚洲精品在线观看二区| 国产成人影院久久av| 老司机在亚洲福利影院| 日日干狠狠操夜夜爽| 欧美最新免费一区二区三区 | 黄色片一级片一级黄色片| 精华霜和精华液先用哪个| 99热这里只有精品一区| 亚洲熟妇熟女久久| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人影院久久av| 波多野结衣高清无吗| 亚洲av不卡在线观看| 国产高清视频在线观看网站| 国产精品亚洲美女久久久| 最后的刺客免费高清国语| 日韩免费av在线播放| 欧美极品一区二区三区四区| 国产主播在线观看一区二区| 亚洲不卡免费看| 国产主播在线观看一区二区| 日本在线视频免费播放| 久久婷婷人人爽人人干人人爱| svipshipincom国产片| 国内揄拍国产精品人妻在线| 日本三级黄在线观看| 他把我摸到了高潮在线观看| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美精品v在线| av片东京热男人的天堂| 午夜免费成人在线视频| 欧美黄色淫秽网站| 亚洲国产精品999在线| 色综合欧美亚洲国产小说| 国产一区在线观看成人免费| 欧美极品一区二区三区四区| 欧美av亚洲av综合av国产av| 欧美中文日本在线观看视频| 久久人人精品亚洲av| 久久久国产成人精品二区| 村上凉子中文字幕在线| 老司机午夜福利在线观看视频| 日韩欧美精品v在线| 蜜桃亚洲精品一区二区三区| 中文字幕人成人乱码亚洲影| 亚洲国产精品久久男人天堂| 久久久成人免费电影| 99久久成人亚洲精品观看| 91麻豆精品激情在线观看国产| 国产 一区 欧美 日韩| 高清日韩中文字幕在线| 国产在线精品亚洲第一网站| 18禁国产床啪视频网站| 国产高清三级在线| АⅤ资源中文在线天堂| 精品国内亚洲2022精品成人| 亚洲av成人av| 日日摸夜夜添夜夜添小说| 九色成人免费人妻av| 小蜜桃在线观看免费完整版高清| 嫩草影院精品99| 搞女人的毛片| 久久亚洲真实| 亚洲中文字幕日韩| 一进一出抽搐动态| 国产蜜桃级精品一区二区三区| 黑人欧美特级aaaaaa片| 床上黄色一级片| 亚洲av美国av| 国产成人av教育| 亚洲国产精品999在线| 国产成人啪精品午夜网站| 一进一出好大好爽视频| 天堂网av新在线| 久久久国产精品麻豆| 欧美绝顶高潮抽搐喷水| www.www免费av| 又粗又爽又猛毛片免费看| 欧美中文日本在线观看视频| 国产高清激情床上av| 精品久久久久久久久久免费视频| 久久久久久久午夜电影| 熟妇人妻久久中文字幕3abv| 亚洲七黄色美女视频| 手机成人av网站| 男插女下体视频免费在线播放| 色哟哟哟哟哟哟| 国产精品美女特级片免费视频播放器| eeuss影院久久| 真实男女啪啪啪动态图| 久久久精品欧美日韩精品| 黄色成人免费大全| 99riav亚洲国产免费| 悠悠久久av| 国产成人aa在线观看| 精华霜和精华液先用哪个| 日本黄色片子视频| 国产精品,欧美在线| 中文字幕高清在线视频| 老鸭窝网址在线观看| 欧美一区二区精品小视频在线| 久久精品国产99精品国产亚洲性色| 久久中文看片网| 91九色精品人成在线观看| 免费观看人在逋| 亚洲色图av天堂| 俺也久久电影网| 亚洲真实伦在线观看| 日韩欧美国产在线观看| 淫秽高清视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 女警被强在线播放| 国产一区在线观看成人免费| 久久人妻av系列| 老汉色av国产亚洲站长工具| 又黄又粗又硬又大视频| 日韩高清综合在线| 99国产极品粉嫩在线观看| 色尼玛亚洲综合影院| 日本精品一区二区三区蜜桃| 免费一级毛片在线播放高清视频| 欧美性猛交╳xxx乱大交人| 老司机深夜福利视频在线观看|