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

    基于增強層次注意熵和極限學習機的軸承故障檢測方法*

    2023-12-20 12:38:48朱文軒張書維
    機電工程 2023年12期
    關鍵詞:特征提取準確率軸承

    朱文軒,張書維,王 琳

    (江蘇安全技術職業(yè)學院 機械工程學院,江蘇 徐州 221011)

    0 引 言

    滾動軸承是旋轉機械中最為常見的零部件,大約有40%的機械故障是由滾動軸承故障引發(fā)的[1]。滾動軸承發(fā)生故障后,輕則會加劇設備的振動噪聲,影響工作人員的操作舒適性和加工精度,重則造成設備無法工作并引發(fā)安全事故。因此,研究滾動軸承的故障診斷方法對確保安全穩(wěn)定生產具有重大意義。

    滾動軸承在運行時受到負載、摩擦和阻尼等非線性因素的影響,導致采集的振動信號表現(xiàn)出強烈的非線性和非平穩(wěn)性,如何從非線性的時間序列中提取出高質量的故障特征,成為滾動軸承故障診斷的關鍵[2-3]。

    目前,由于理論解釋性好,特征提取性能優(yōu)異,基于熵理論的特征提取方法被廣泛應用于滾動軸承的故障診斷領域[4]。HAN Ming-hong等人[5]將局部均值分解和樣本熵相結合,進行了故障特征提取,有效分辨了軸承的不同故障狀態(tài);但是,樣本熵在分析長時間序列時,故障特征提取效率較為低下。隨后,楊云等人[6]提出了一種基于變分模態(tài)分解和排列熵的滾動軸承故障診斷模型,借此模型有效識別了軸承的故障;然而排列熵忽視了數(shù)據(jù)振幅之間的區(qū)別,導致了軸承故障信息的丟失。為此,陳劍等人[7]提出了基于改進固有時間尺度分解和模糊熵的滾動軸承故障診斷模型,借此模型對滾動軸承不同故障進行了準確分類;然而模糊熵的特征提取效率非常低下。為緩解上述熵的缺陷,ROSTAGHI M等人[8]提出了一種新的測量時間序列復雜度的方法——散布熵,其不僅能準確區(qū)分信號的幅值差異,而且可以顯著提高故障分析的效率。葉震等人[9]將模糊信息?;蜕⒉检叵嘟Y合,用于滾動軸承的故障診斷,有效識別了軸承的不同故障類型。

    無論是排列熵、模糊熵還是散布熵,在特征提取階段均需要設置2個及以上的超參數(shù),削弱了算法的泛化性。為此,YANG Jia-wei等人[10]提出了一種新的測量時間序列復雜性的工具——注意熵(attention entropy,ATE)。與傳統(tǒng)熵值關注信號中全部數(shù)據(jù)的頻率分布不同,ATE僅關注了信號峰值點間隔的頻率分布情況,因此ATE具有參數(shù)少、運行效率高和對數(shù)據(jù)長度不敏感的優(yōu)點。

    然而,ATE只提取了信號的單個尺度特征,滾動軸承的振動包含多個模態(tài)特征,僅開展單一尺度分析,無法全面提取滾動軸承的故障特征頻率信息[11]。為此,陳飛等人[12]359-360提出了多尺度注意熵(multiscale attention entropy,MATE)和時移多尺度注意熵(time-shifted multiscale attention entropy,TSMATE),并將其用于水電機組的故障診斷,有效提取了多個尺度的故障特征,從而準確區(qū)分了不同故障狀態(tài);然而MATE的粗?;幚泶嬖谳^大的缺陷,丟失了大量故障信息,且造成了較大的偏差。CHEN F等人[13]對粗?;M行了改進,提出了改進多尺度注意熵(improved multiscale attention entropy,IMATE),其有效區(qū)分了旋轉機械的故障狀態(tài);但是,上述基于粗?;幚淼淖⒁忪刂笜酥荒芴崛≌駝有盘栔械牡皖l故障特征,卻忽略了信號中同樣重要的高頻故障特征[14]。

    為提取滾動軸承振動信號的低頻和高頻特征信息,朱可恒等人[15]通過層次分析進行了多尺度計算,提出了基于層次樣本熵的滾動軸承故障診斷方法,對滾動軸承故障進行了準確識別;然而傳統(tǒng)的層次分析面臨著和傳統(tǒng)粗?;幚眍愃频娜毕?即隨著分解層數(shù)的增加,層次序列的長度會顯著減少,使得分析誤差增大且丟失特征信息[16]。為此,柏世兵等人[17]對傳統(tǒng)的層次分析方法進行了改進,提出了一種改進的層次分析方法,并采用軸承的診斷結果驗證了改進方法的有效性;但其采用的熵值指標需要設置較多的超參數(shù)。

    在損傷識別方面,極限學習機(ELM)具有較高效率,同時還具有較強的分類性能,因此被廣泛應用于滾動軸承的故障識別領域[18]。但是ELM的輸入層權重和隱含層偏置需要人為設定,為此,筆者引入基于灰狼算法(GWO)優(yōu)化的極限學習機(GWO-ELM)代替原始的ELM,進行特征樣本的識別。

    綜上所述,為了實現(xiàn)在較少參數(shù)下提取非線性、非平穩(wěn)振動信號的故障特征,筆者將增強的層次分析和注意熵相結合,提出增強層次注意熵(EHATE);并在此基礎上,提出一種基于EHATE和GWO-ELM的滾動軸承故障診斷方法(其中,EHATE方法用于提取滾動軸承振動信號的低頻和高頻特征信息,而GWO-ELM用于進行滾動軸承的故障識別);利用3組滾動軸承數(shù)據(jù)集對該方法進行實驗驗證和分析,并開展多個方面的對比研究。

    1 增強層次注意熵

    目前,排列熵、模糊熵、散布熵等指標被廣泛用于軸承的故障診斷領域,但這些方法的多個參數(shù)需要預先設置,算法的泛化性較差。同時,為了實現(xiàn)信號的多維度分析,以粗?;幚頌榇淼亩喑叨确椒捌涓倪M方法應用廣泛,但這些方法無法提取信號的高頻特征,遺漏了大量故障信息。

    為此,筆者充分考慮了這些方法的缺陷,提出了增強層次注意熵(EHATE),其不僅避免了算法在特征提取過程中需要設置大量參數(shù)的缺陷,而且能夠提取信號的高頻特征。

    1.1 注意熵

    和傳統(tǒng)的熵值方法相比,注意熵(ATE)對分析數(shù)據(jù)的長度不敏感,且不需要設置任何的超參數(shù),優(yōu)勢較為顯著,其實現(xiàn)過程如圖1所示。

    圖1 注意熵的示意圖Fig.1 The diagram of attention entropy

    ATE的計算過程可總結如下:

    1)若振動信號中的每個數(shù)據(jù)點都被假定為一個系統(tǒng),其狀態(tài)的變化可以視為系統(tǒng)對環(huán)境的調整。峰值點可以準確反映局部狀態(tài)的上下邊界的變化情況,因而將局部峰值點定義為核心點;

    2)按照{min-min}、{min-max}、{max-min}和{max-max}4種不同策略給定核心點,并獲取相鄰核心點的間隔;

    3)求解相鄰核心點間隔香農熵,公式如下:

    (1)

    式中:p(x)為x出現(xiàn)的概率;b為間隔點種類個數(shù);

    4)將4種不同策略求解得到的香農熵均值定義為該振動信號的注意熵。

    1.2 增強層次注意熵

    ATE只是從單個尺度描述信號的復雜性,無法全面地表征信號的固有特性。為此,COSTA M等人[19]提出了多尺度熵的概念,將時間序列分割為多個長度較短的子序列,實現(xiàn)了從多個尺度分析信號的固有特性目的。然而,基于粗?;幚淼亩喑叨扔嬎惴椒ㄖ皇翘崛⌒盘柕牡皖l特征,忽略了信號的高頻特征。為此,JIANG Ying等人[20]對時間序列進行層次化處理,提出了平均算術符和差分算術符的概念,對信號進行多尺度分析的同時還提取了信號的低頻和高頻特征。但隨著分解層數(shù)的增加,信號的長度會縮短,導致其可靠性下降。此外,傳統(tǒng)層次分析對信號長度有特殊要求,限制了算法的廣泛應用。

    為此,筆者將層次熵和注意熵相結合,使用移動平均和移動差分過程來取代原始的層次分解過程,提出增強層次注意熵(EHATE),以緩解傳統(tǒng)層次熵層次分解不徹底和參數(shù)難以準確設置等問題。

    EHATE的原理如下:

    1)對于信號{x(i),i=1,2,…,L},定義平均算術符Q0和差分算術符Q1如下:

    (2)

    式中:Q0(x)為信號x(i)的低頻分量;Q1(x)為信號x(i)的高頻分量;

    (3)

    3)構造向量{γ1,γ2,…,γk},則整數(shù)e可以定義如下:

    (4)

    式中:γm為第m層的平均或差分運算符;

    4)基于向量{γ1,γ2,…,γk},則信號x(i)的層次分量定義為:

    (5)

    5)求解每個層次分量的注意熵,則增強層次注意熵可以定義為:

    EHATE(x,k)=ATE(Xk,e)

    (6)

    1.3 參數(shù)分析和仿真對比

    在EHATE中,需要預先輸入的參數(shù)是層次分解層數(shù)k,同時信號x(i)的長度L也需要合理的設置。通常來說,參數(shù)k設定為3或4,當小于3時,所提取的特征數(shù)量太少,無法體現(xiàn)信號在多個尺度上的復雜性,導致信號固有特性表征的不夠全面;反之當k大于4時,會導致特征的數(shù)量過多,如k=5時,特征的數(shù)量達到了32個,這不僅增加了特征提取時間,而且增加了特征冗余度[21]。

    因此,筆者將k設置為3,后續(xù)將就k的取值對算法性能的影響做進一步分析。

    信號長度L對EHATE性能的影響較小,但長度通常不小于1 024,否則樣本中包含的峰值點較少,無法準確表征信號的趨勢和固有特性。

    在陳飛等人[12]361-362的研究中,對不同長度1/f噪聲的TSMATE和MATE性能進行了研究,結果表明TSMATE對數(shù)據(jù)長度不敏感。

    為此,筆者在前人工作的基礎上,對不同長度白噪聲(white gaussian nosie,WGN)和1/f噪聲的EHATE、層次注意熵(hierarchical attention entropy,HATE)和MATE進行了分析。

    3種方法的噪聲分析結果隨長度的變化如圖2所示。

    圖2 EHATE、HATE和MATE的分布情況Fig.2 Distribution of EHATE, HATE and MATE

    在圖2中,對比2種不同長度噪聲的EHATE、HATE和MATE的結果發(fā)現(xiàn),隨著數(shù)據(jù)長度的減小,HATE和MATE的波動逐漸變大,特別是在節(jié)點較大的時候熵值發(fā)生了突變(以數(shù)據(jù)長度L=256的WGN噪聲信號為例,在層次節(jié)點為5、6和7時,HATE熵值變化非常劇烈)。

    此外,不同長度噪聲信號的HATE和MATE值分布的差異過大,對于同一節(jié)點,WGN和1/f噪聲的HATE和MATE波動的最大幅值分別達到了0.626、0.481、0.790和0.445。

    以上結果證明,在分析同一種信號時,短信號的分析結果和長信號的HATE和MATE分析結果存在較大的差異。

    而不同長度噪聲信號的EHATE的分布差異較小,在各個節(jié)點上僅存在較小的波動,WGN和1/f噪聲的EHATE的最大波動幅值僅為0.393和0.471,說明EHATE在分析不同長度的信號時能產生一致的熵值分布;并且,隨著節(jié)點的增加,EHATE值重疊得更加緊密,表明EHATE對信號長度有著較強的魯棒性。

    為了進一步評估EHATE、HATE和MATE的分析穩(wěn)定性,筆者首先構造了50組長度為2 048的WGN和1/f噪聲信號,隨后利用3種方法計算這50組噪聲信號的熵值,并提取其標準差,結果如圖3所示。

    圖3 WGN和1/f信號的EHATE,HATE和MATE分析標準差Fig.3 Analysis standard deviations of EHATE, HATE and MATE for WGN and 1/f signals

    由圖3可知:在尺度大于2后,EHATE的標準差就小于HATE和MATE方法的標準差,趨于一個較小的水平。對于1/f信號而言,HATE和MATE的標準差隨著節(jié)點的增加而逐漸增大,特別是在節(jié)點為7時,HATE和MATE的標準差達到了最大,表明這2種方法的分析誤差在節(jié)點數(shù)較大時會比較大,結果不夠穩(wěn)定。

    雖然,在節(jié)點數(shù)小于2時,EHATE的標準差大于MATE方法的標準差,但總體上,通過構造平均和差分算術符,使得EHATE具有提取信號低頻和高頻特征的能力,因而優(yōu)于MATE。

    2 基于EHATE和GWO-ELM的診斷方法

    2.1 灰狼算法優(yōu)化的極限學習機

    2.1.1 極限學習機

    極限學習機(ELM)為單隱層的神經網絡,若有N個隨機樣本(xi,yi)|xi∈Rm,yi∈Rn,則包含L個隱藏層節(jié)點的ELM模型可表示如下:

    (7)

    式中:βj為鏈接第j個隱藏層節(jié)點與模型輸出節(jié)點的輸出權重矢量;g(x)為隱藏層神經元的激活函數(shù);ωj為鏈接模型輸出節(jié)點與第j個隱藏層節(jié)點的輸入權重矢量;bj為第j個隱藏層節(jié)點的偏置;yi為樣本的模型輸出,yi=[yi1,yi2,…,yin]T。

    (8)

    上式變?yōu)榫仃囆问?如下:

    (9)

    此時,輸出權重矩陣β=[β1,β2,…,βL]T可定義如下:

    β=H+T

    (10)

    式中:H+為隱藏層輸出矩陣H的廣義逆矩陣。

    2.1.2 優(yōu)化流程

    極限學習機的分類性能優(yōu)異,并且效率高,但是其2個核心參數(shù)需要人為設置,這導致算法容易過擬合或欠擬合,降低算法的性能。為此,筆者引入了基于灰狼算法優(yōu)化的極限學習機,進行特征的識別。

    其優(yōu)化流程如下:

    1)對原始數(shù)據(jù)進行歸一化預處理,并將樣本隨機地劃分出訓練樣本和測試樣本;

    2)設置GWO算法的種群規(guī)模和迭代次數(shù),將ELM的輸入權重和隱藏層閾值設置為優(yōu)化目標,適應度函數(shù)設置為故障識別的平均錯誤率,通過不斷地迭代,以搜索錯誤率最小時所對應的參數(shù);

    3)計算灰狼個體的適應度函數(shù),并不斷更新參數(shù),如果迭代過程中符合終止條件,則輸出對應的最優(yōu)參數(shù),否則繼續(xù)執(zhí)行迭代;

    4)將迭代得到的最優(yōu)參數(shù)代入至ELM模型中,并利用樣本對ELM進行訓練和測試,輸出測試樣本的故障識別結果。

    GWO-ELM的詳細理論以及參數(shù)設置可參考欒孝馳等人[22]的研究。

    2.2 故障診斷流程

    基于前面的分析,EHATE只需要輸入2個參數(shù),分別是信號數(shù)據(jù)長度L和分解層數(shù)k。對于信號數(shù)據(jù)長度,分析結果表明EHATE的魯棒性較強,對于數(shù)據(jù)長度不敏感,為此筆者設置數(shù)據(jù)長度L=2 048。在該數(shù)據(jù)長度下,2種噪聲信號的EHATE、HATE和MATE分析結果表明,EHATE的分析穩(wěn)定性優(yōu)于另外2種方法。

    因此,鑒于EHATE所需參數(shù)少、性能優(yōu)異的優(yōu)點,結合已被證明具有優(yōu)異性能的GWO-ELM分類器,筆者提出一種基于EHATE和GWO-ELM的滾動軸承故障診斷方法。

    具體診斷流程如下所示:

    1)利用振動加速度計采集滾動軸承的不同工況下的振動信號,將其作為原始的數(shù)據(jù)樣本;

    2)將每個工況的振動信號樣本平均分割成長度為2 048的50組樣本,其中25組樣本作為訓練樣本,剩余25組樣本作為測試樣本;

    3)利用EHATE提取訓練和測試樣本的前8個熵值,構造故障特征;

    4)將訓練樣本的故障特征輸入至GWO-ELM模型進行訓練,待訓練完畢后輸入測試樣本,以進行故障識別;

    5)根據(jù)GWO-ELM的輸出標簽對滾動軸承樣本的故障類型和故障程度進行診斷,并進行模型有效性的評估。

    3 故障診斷實驗和結果分析

    為了驗證EHATE+GWO-ELM方法的有效性,并評估該模型對2個參數(shù)的敏感性,筆者分別利用江南大學[23]、東南大學[24]和印度某大學的滾動軸承數(shù)據(jù)集[25]進行實驗驗證。

    3.1 江南大學軸承數(shù)據(jù)

    江南大學的滾動軸承數(shù)據(jù)集包含4種不同的工況,分別是健康(Nor)、內圈故障(IR)、外圈故障(OR)和滾動體故障(BA);數(shù)據(jù)的采樣頻率為50 kHz;軸承的轉速包含600 r/min、800 r/min和1 000 r/min 3種,筆者選擇800 r/min下的數(shù)據(jù)進行實驗。訓練樣本為100組,測試樣本為100組。

    4種軸承樣本的振動信號的時域波形如圖4所示。

    圖4 滾動軸承振動信號的波形Fig.4 Waveform of vibration signal of rolling bearing

    由圖4可以發(fā)現(xiàn):不同振動信號中所包含的沖擊成分不一致,具有較大的差異性,但無法直接根據(jù)波形確定區(qū)分準則。

    筆者利用EHATE提取滾動軸承振動信號的故障特征,隨后將故障特征輸入至GWO-ELM模型中,進行故障識別,得到滾動軸承故障診斷結果,如圖5所示。

    從圖5(a)可以發(fā)現(xiàn):EHATE特征曲線在各個節(jié)點都具有相同的趨勢,無法直接根據(jù)特征曲線來分辨樣本的故障類型。

    圖5 EHATE故障特征和對應的故障識別結果Fig.5 EHATE fault characteristics and corresponding fault identification results

    而從圖5(b)可以發(fā)現(xiàn):EHATE+GWO-ELM方法取得了100%的識別準確率,證明該方法可以有效地識別滾動軸承的不同故障類型。

    隨后,為了驗證EHATE方法的特征提取性能對數(shù)據(jù)長度L和分解層數(shù)k2個不同參數(shù)的敏感性,筆者對不同參數(shù)組合下EHATE+GWO-ELM方法的分類結果進行了評估,結果如圖6所示。

    由圖6可知:隨著數(shù)據(jù)長度和分解層數(shù)的增加,EHATE+GWO-ELM方法的分類準確率也隨之增加,表明這2個參數(shù)對算法的性能影響是正向的。當數(shù)據(jù)長度大于2 048后,即使分解層數(shù)為1,準確率也能夠達到90%以上;當分解層數(shù)為2時,準確率也達到了98%;當分解層數(shù)為3時,準確率達到了100%,證明該方法僅需要常規(guī)的參數(shù)設置(L=2 048,k=3)即可獲得較為可靠的故障識別結果。

    為了評估EHATE+GWO-ELM方法的優(yōu)越性,筆者利用HATE、MATE、TSMATE、精細復合多尺度樣本熵(refined composite multiscale sample entropy,RCMSE)、精細復合廣義多尺度模糊熵(refined composite generalized multiscale fuzzy entropy,RCGMFE)、層次排列熵(hierarchical permutation entropy,HPE)、增強層次散布熵(enhanced hierarchical dispersion entropy,EHDE)進行對比,將上述方法提取的故障特征輸入至GWO-ELM中進行故障識別,并重復10次,同時統(tǒng)計每種方法的特征提取時間和平均分類準確率,結果如表1所示。

    表1 不同故障特征提取方法的診斷結果對比

    Table 1 Comparison of diagnosis results of different fault feature extraction methods

    特征提取方法參數(shù)設置平均準確率/%標準差特征提取時間/sEHATEk=398.201.1493.52HATEk=372.102.9231.17MATEs=896.002.113.92TSMATEs=896.700.9512.71RCMSEm=2;r=0.2;t=1;s=885.803.9966.10RCGMFEm=2;r=0.2;n=2;t=1;s=872.802.20122.55HPEm=5;t=1;k=368.402.6857.13EHDEm=3;c=6;t=1;k=398.501.08108.88

    由表1可知:在準確率方面,EHDE的平均準確率最高,達到了98.50%,EHATE的準確率次之,為98.20%,相差較小;在特征提取時間方面,MATE的效率最高,為3.92 s,對應的平均準確率為96%;在標準差方面,TSMATE的標準差最小,為0.95,對應的準確率為96.70%。雖然EHDE的準確率高于EHATE,但其特征提取效率低于EHATE,同時其需要額外設置3個參數(shù),該算法對參數(shù)的敏感性較強。而MATE和TSMATE方法不僅具有較高的準確率,同時特征提取的效率也較高。

    綜合上述對比結果,雖然EHATE的特征提取時間增加了,但在準確率和穩(wěn)定性等多個指標上都優(yōu)于其他方法。此外,最關鍵的是EHATE設置的參數(shù)少,這提高了該方法的泛化性。雖然MATE、TSMATE、HATE方法都只需要設置一個參數(shù),但EHATE方法的準確率明顯高于這幾種方法。

    綜合來看,該方法的性能是比較優(yōu)秀的。

    3.2 東南大學軸承數(shù)據(jù)

    為了進一步驗證EHATE+GWO-ELM方法的有效性,筆者利用東南大學提供的滾動軸承數(shù)據(jù)進行實驗。實驗平臺如圖7所示。

    圖7 東南大學滾動軸承故障模擬實驗臺Fig.7 Rolling bearing fault simulation test bench in Southeast University

    在該實驗平臺中,振動傳感器布置在行星齒輪箱表面,以5 120 Hz的頻率采集滾動軸承的振動信號,筆者利用模擬器將系統(tǒng)的轉速和負載調整至20 Hz和0 V,將健康軸承(NOR)分別替換為存在滾珠故障(BA)、內圈故障(IR)、外圈故障(OR)、內圈&外圈故障(IO)的軸承,以此來模擬系統(tǒng)不同狀態(tài)下的運行情況。

    訓練樣本為125組,測試樣本為125組,5種軸承樣本的振動信號的時域波形如圖8所示。

    圖8 滾動軸承振動信號的波形Fig.8 Waveform of vibration signal of rolling bearing

    由圖8可以發(fā)現(xiàn):振動信號波形比較復雜,缺乏顯著的規(guī)律,有必要采用EHATE+GWO-ELM方法進行樣本的故障識別。

    筆者利用EHATE提取滾動軸承振動信號的故障特征,隨后,將故障特征輸入至GWO-ELM模型中進行故障識別,診斷結果如圖9所示。

    由圖9可知:EHATE熵值特征在各個節(jié)點上具有比較好的區(qū)分度,證明利用EHATE提取振動信號的低頻和高頻故障特征,能夠較好地區(qū)分滾動軸承的不同故障狀態(tài)。

    圖9(b)的GWO-ELM故障識別結果也驗證了EHATE+GWO-ELM方法的良好性能,僅IO樣本出現(xiàn)了一個被錯誤識別的樣本,其他樣本都被準確地識別出來,總的準確率為99.2%,證明該方法能夠準確地識別滾動軸承樣本的故障類型。

    為了進一步驗證EHATE方法的性能對數(shù)據(jù)長度L和分解層數(shù)k的參數(shù)敏感性,筆者對不同參數(shù)組合下EHATE+GWO-ELM方法的性能進行了評估,如圖10所示。

    圖10 不同參數(shù)組合下EHATE+GWO-ELM方法的診斷準確率Fig.10 The diagnostic accuracy of the EHATE+GWO-ELM method under different parameter combinations

    由圖10可知:在數(shù)據(jù)長度較小時,EHATE+GWO-ELM方法的準確率較低,即使特征數(shù)量增加到32個,準確率也僅為85.6%。

    而當數(shù)據(jù)長度為4 096時,即使特征的數(shù)量只有2個,準確率也達到了84%,這證明數(shù)據(jù)長度是影響算法性能的主要因素,若同時增加數(shù)據(jù)長度和特征的數(shù)量,則準確率會得到顯著提高。

    當數(shù)據(jù)長度為2 048,分解層數(shù)為3時,該方法的準確率已經可以達到99.2%,證明此時該方法能夠有效地識別滾動軸承的故障類型,同時進一步證明了該方法在常規(guī)的參數(shù)設置下能夠獲得較為可靠的故障識別結果。

    為了進一步評估EHATE+GWO-ELM方法的優(yōu)越性,筆者利用前面的7種方法進行同樣的對比,結果如表2所示。

    表2 不同故障特征提取方法的診斷結果對比

    由表2可知:TSMATE的平均準確率最高,達到了99.28%,同時標準差也只有0.79,優(yōu)于EHATE方法,證明時移粗?;幚硪材軌蛴行У靥崛L動軸承振動信號中的故障信息。除了TSMATE方法之外,EHATE方法的準確率最高,同時穩(wěn)定性也較高,證明了該方法的有效性。

    綜合來看,EHATE+GWO-ELM方法的準確率比較高,雖然效率不高,特征提取的時間僅少于EHTE和RCGMFE,但由于其對參數(shù)不敏感,因此具有相對均衡的性能。

    3.3 印度大學軸承數(shù)據(jù)

    為更進一步驗證EHATE+GWO-ELM方法在診斷滾動軸承不同故障類型和不同嚴重程度中的有效性,筆者利用印度Sant Longowal工程技術學院提供的滾動軸承數(shù)據(jù)集進行實驗。

    滾動軸承實驗平臺如圖11所示。

    圖11 滾動軸承實驗平臺Fig.11 Rolling bearing test platform

    實驗中,滾動軸承的型號為NU205E的滾子軸承,筆者利用加速度計,以70 kHz的頻率采集其在運行時的振動信號(其中,滾動軸承的轉速為2 050 r/min,負載為200 N)。

    為了模擬存在不同軸承故障時系統(tǒng)的運行狀態(tài),筆者將正常軸承替換為不同故障類型和不同故障程度的軸承,進行振動測量實驗。

    實驗中,共有4種類型的軸承,分別是正常軸承、滾珠故障軸承、內圈故障軸承和外圈故障軸承。對于3種故障類型的軸承,存在4種故障程度,分別是輕度、中度、重度和超重度,因此,總共有13種故障類型的軸承樣本。

    樣本的詳細信息如表3所示。

    表3 滾動軸承樣本的故障信息

    13種軸承樣本的振動信號的時域波形如圖12所示。

    由圖12可以發(fā)現(xiàn):由于故障狀態(tài)較多,且波形成分復雜,缺乏顯著的規(guī)律,無法通過人為觀察來判斷樣本的狀態(tài),有必要采用EHATE+GWO-ELM方法進行樣本的故障識別。

    圖12 滾動軸承振動信號的波形Fig.12 Waveform of vibration signal of rolling bearing

    筆者利用EHATE提取滾動軸承樣本的故障特征,隨后將故障特征輸入至GWO-ELM中,進行故障識別,結果如圖13所示。

    圖13 EHATE故障特征和對應的故障識別結果Fig.13 EHATE fault characteristics and corresponding fault identification results

    由圖13(a)可知:不同樣本的熵值特征曲線具有類似的趨勢,并且由于樣本類型較多,曲線互相重疊,難以直接根據(jù)波形判斷樣本的狀態(tài)。

    由圖13(b)可知:EHATE方法的GWO-ELM診斷準確率為96.92%,出現(xiàn)了10個被錯誤識別的樣本,有5種類型的樣本出現(xiàn)了錯誤樣本(即:4個BA2樣本被錯誤識別為BA4,2個BA3樣本被錯誤識別為BA2,1個BA4樣本被錯誤識別為BA2,1個IR3樣本被錯誤識別為IR4,1個OR3樣本被錯誤識別為IR3,1個OR3樣本被錯誤識別為IR4),剩余8個類型的樣本都被準確地分類。

    根據(jù)該診斷結果,可以驗證EHATE+GWO-ELM方法能夠較為有效地識別滾動軸承不同故障類型和故障程度,具有一定的應用潛力。

    同樣地,筆者對EHATE方法在不同參數(shù)設置下的特征提取性能進行了研究,結果如圖14所示。

    圖14 不同參數(shù)組合下所提方法的診斷準確率Fig.14 The diagnostic accuracy of the proposed method under different parameter combinations

    由圖14可知:當數(shù)據(jù)長度不小于1 024,分解層數(shù)不小于3層時,EHATE+GWO-ELM方法能夠取得了不低于90%的分類準確率;同時隨著參數(shù)的增加,準確率隨之增加,但特征提取的時間也隨之顯著提高,因此需要在效率和性能之間進行綜合考慮。

    當數(shù)據(jù)長度為2 048,分解層數(shù)為3時,該方法的準確率為96.92%;而當數(shù)據(jù)長度為4 096,分解層數(shù)為4時,準確率達到了100%,然而此時特征提取時間是前者的7倍。

    因此,該方法只需要采用常規(guī)的參數(shù)組合,即可以獲得較高的性能和效率。

    同樣的,筆者對前述7種方法進行了多次實驗,取平均準確率和特征提取時間進行對比,結果如表4所示。

    表4 不同故障特征提取方法的診斷結果對比

    由表4可知:EHATE的平均分類準確率最高,達到了95.26%,其次是TSMATE方法,為89.23%;而EHATE方法的標準差也最小,為0.88;特征提取效率最高的是MATE方法,為15.39 s,而EHATE方法的效率為370.07 s。

    綜合準確率、穩(wěn)定性和效率3個方面而言,EHATE+GWO-ELM方法的綜合性能最優(yōu)。因此,該方法能夠穩(wěn)定且有效地診斷滾動軸承的不同故障類型和故障程度。

    此外,EHATE方法由于僅需設置較少的參數(shù),更適合用于故障特征提取階段。

    4 結束語

    鑒于常規(guī)滾動軸承故障診斷方法面臨準確率不高、穩(wěn)定性差和參數(shù)設置較多的缺陷,筆者提出了一種基于EHATE和GWO-ELM的滾動軸承故障診斷方法,并利用3種滾動軸承故障數(shù)據(jù)集對該方法進行了實驗和分析。

    研究結論如下:

    1)EHATE方法由于僅需要設置數(shù)據(jù)長度和分解層數(shù)2個參數(shù),其相較于其他方法對參數(shù)不敏感,采用常規(guī)的參數(shù)設置即可以獲得良好的故障識別結果;當樣本的數(shù)據(jù)長度為N=2 048,而分解層數(shù)k=3時,準確率能夠達到95%以上;

    2)EHATE的性能優(yōu)異,能夠有效地反映信號的復雜度特性,仿真信號的分析結果驗證了EHATE優(yōu)于HATE和MATE;

    3)基于EHATE和GWO-ELM的故障診斷方法能夠有效地診斷滾動軸承的不同故障類型,識別準確率分別達到了100%、99.2%和96.92%,而多次實驗的平均識別準確率分別為98.2%、99.12%和95.26%,在識別穩(wěn)定性和識別準確率方面優(yōu)于其他7種對比方法,更適合用于滾動軸承的故障診斷。

    雖然EHATE方法取得了不錯的故障識別結果,但是其特征提取效率有待提高。后續(xù)筆者將對粗粒化方式進行優(yōu)化,以進一步提高特征提取的效率。

    猜你喜歡
    特征提取準確率軸承
    軸承知識
    哈爾濱軸承(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
    乳腺超聲檢查診斷乳腺腫瘤的特異度及準確率分析
    健康之家(2021年19期)2021-05-23 11:17:39
    不同序列磁共振成像診斷脊柱損傷的臨床準確率比較探討
    2015—2017 年寧夏各天氣預報參考產品質量檢驗分析
    基于Daubechies(dbN)的飛行器音頻特征提取
    電子制作(2018年19期)2018-11-14 02:37:08
    高速公路車牌識別標識站準確率驗證法
    Bagging RCSP腦電特征提取算法
    黄色日韩在线| 深夜精品福利| 国产伦一二天堂av在线观看| 日本一二三区视频观看| 免费观看人在逋| 极品教师在线免费播放| 成年女人永久免费观看视频| 欧美黑人欧美精品刺激| 国产亚洲91精品色在线| 国内揄拍国产精品人妻在线| 午夜福利高清视频| 啦啦啦韩国在线观看视频| 久久精品国产99精品国产亚洲性色| 色视频www国产| 一级黄片播放器| 国产成人av教育| 能在线免费观看的黄片| 五月玫瑰六月丁香| 亚洲av日韩精品久久久久久密| 欧美在线一区亚洲| 最近中文字幕高清免费大全6 | 国产亚洲欧美98| 久久久国产成人精品二区| 国产成人aa在线观看| 简卡轻食公司| 日本欧美国产在线视频| 久99久视频精品免费| 两人在一起打扑克的视频| 亚洲人与动物交配视频| 男人的好看免费观看在线视频| 一个人观看的视频www高清免费观看| 男女视频在线观看网站免费| 性插视频无遮挡在线免费观看| 久久久久国内视频| 欧美色视频一区免费| 18禁黄网站禁片午夜丰满| aaaaa片日本免费| 欧美中文日本在线观看视频| 可以在线观看的亚洲视频| 欧美xxxx性猛交bbbb| 亚洲成人久久爱视频| 久久久国产成人免费| 国产美女午夜福利| 中文字幕精品亚洲无线码一区| 琪琪午夜伦伦电影理论片6080| 国产成人av教育| 蜜桃亚洲精品一区二区三区| 高清日韩中文字幕在线| 制服丝袜大香蕉在线| 免费电影在线观看免费观看| 国产精品人妻久久久影院| 精品久久国产蜜桃| 亚洲性夜色夜夜综合| 免费看美女性在线毛片视频| 我要看日韩黄色一级片| 国产伦人伦偷精品视频| 男女视频在线观看网站免费| 国产亚洲精品久久久com| 在线免费观看不下载黄p国产 | 村上凉子中文字幕在线| 婷婷六月久久综合丁香| 亚洲18禁久久av| 露出奶头的视频| 18禁在线播放成人免费| 一边摸一边抽搐一进一小说| 可以在线观看的亚洲视频| 国内久久婷婷六月综合欲色啪| 精品久久久噜噜| 国产精品亚洲美女久久久| 国产精品亚洲一级av第二区| 国产精品福利在线免费观看| 国产91精品成人一区二区三区| 国产精品日韩av在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 成人av在线播放网站| 狂野欧美白嫩少妇大欣赏| 久久精品国产99精品国产亚洲性色| 别揉我奶头 嗯啊视频| av在线老鸭窝| 国产在线男女| 又紧又爽又黄一区二区| 亚洲精品在线观看二区| 日韩精品有码人妻一区| 十八禁网站免费在线| 听说在线观看完整版免费高清| 搡老妇女老女人老熟妇| 嫁个100分男人电影在线观看| 久久久久性生活片| av黄色大香蕉| 国产探花极品一区二区| 免费在线观看影片大全网站| 国产午夜精品久久久久久一区二区三区 | 久久久久久国产a免费观看| 欧美一区二区亚洲| 久久久久九九精品影院| 亚洲欧美日韩卡通动漫| bbb黄色大片| 日本在线视频免费播放| 人人妻人人看人人澡| 免费无遮挡裸体视频| 不卡一级毛片| 一区福利在线观看| 超碰av人人做人人爽久久| 国产精品久久视频播放| 少妇人妻精品综合一区二区 | 亚洲美女黄片视频| 天堂动漫精品| 99久久无色码亚洲精品果冻| 欧美日韩中文字幕国产精品一区二区三区| 午夜福利欧美成人| 亚洲熟妇熟女久久| 国产精品亚洲一级av第二区| 又粗又爽又猛毛片免费看| 国产真实伦视频高清在线观看 | 性欧美人与动物交配| 亚洲av不卡在线观看| 联通29元200g的流量卡| 白带黄色成豆腐渣| 麻豆久久精品国产亚洲av| 中文字幕av成人在线电影| 午夜影院日韩av| 99久久精品热视频| 国国产精品蜜臀av免费| 久久6这里有精品| 国产精品美女特级片免费视频播放器| 亚洲人成伊人成综合网2020| av黄色大香蕉| 日韩国内少妇激情av| 中文亚洲av片在线观看爽| 日韩欧美国产一区二区入口| 免费不卡的大黄色大毛片视频在线观看 | 极品教师在线视频| 国内精品久久久久精免费| 亚洲av中文字字幕乱码综合| 色5月婷婷丁香| 夜夜看夜夜爽夜夜摸| 波多野结衣高清无吗| 成人国产麻豆网| 亚洲av一区综合| 1024手机看黄色片| 久久99热6这里只有精品| 欧美国产日韩亚洲一区| 日韩一区二区视频免费看| 亚洲色图av天堂| 日韩欧美国产一区二区入口| 88av欧美| 欧美性猛交黑人性爽| a在线观看视频网站| 中出人妻视频一区二区| 啦啦啦韩国在线观看视频| 亚洲av五月六月丁香网| 日韩欧美国产一区二区入口| 国产欧美日韩精品亚洲av| 欧美zozozo另类| 国产麻豆成人av免费视频| 亚洲第一电影网av| 搡老熟女国产l中国老女人| 99久久久亚洲精品蜜臀av| 精华霜和精华液先用哪个| 国产男人的电影天堂91| 69av精品久久久久久| 最新在线观看一区二区三区| 国产色爽女视频免费观看| 国内精品一区二区在线观看| 身体一侧抽搐| 他把我摸到了高潮在线观看| 欧美性猛交╳xxx乱大交人| videossex国产| 亚洲一区高清亚洲精品| 日日夜夜操网爽| 免费观看在线日韩| 国产精品av视频在线免费观看| 精品久久久久久久久av| 一本久久中文字幕| 日本黄色片子视频| 看片在线看免费视频| 熟女电影av网| www日本黄色视频网| 国产一区二区激情短视频| 亚洲精华国产精华精| 亚洲av不卡在线观看| 国产伦在线观看视频一区| 日韩国内少妇激情av| 12—13女人毛片做爰片一| 99热网站在线观看| 天天躁日日操中文字幕| 国内精品久久久久久久电影| 精品午夜福利视频在线观看一区| 小蜜桃在线观看免费完整版高清| 我的老师免费观看完整版| 日日撸夜夜添| 精品免费久久久久久久清纯| 99精品久久久久人妻精品| 黄色欧美视频在线观看| 在线天堂最新版资源| 国产不卡一卡二| 美女xxoo啪啪120秒动态图| 国产成人aa在线观看| 欧美一级a爱片免费观看看| 在线天堂最新版资源| 国产精品久久久久久精品电影| 亚洲性夜色夜夜综合| 无人区码免费观看不卡| 亚洲在线自拍视频| 国产亚洲精品久久久久久毛片| 欧美+亚洲+日韩+国产| 麻豆国产97在线/欧美| 午夜福利视频1000在线观看| 国产精品伦人一区二区| 亚洲三级黄色毛片| av.在线天堂| 99热只有精品国产| 久久久久九九精品影院| 如何舔出高潮| 欧美激情国产日韩精品一区| 身体一侧抽搐| 国产精品人妻久久久影院| 性色avwww在线观看| 亚洲男人的天堂狠狠| 欧美最黄视频在线播放免费| 网址你懂的国产日韩在线| 婷婷六月久久综合丁香| 国内毛片毛片毛片毛片毛片| 春色校园在线视频观看| 亚洲久久久久久中文字幕| 不卡视频在线观看欧美| 免费在线观看日本一区| 不卡一级毛片| 久久久久久久久久黄片| 国产高清有码在线观看视频| 少妇人妻精品综合一区二区 | 午夜福利欧美成人| 制服丝袜大香蕉在线| 国产精品一区www在线观看 | 最新中文字幕久久久久| 岛国在线免费视频观看| 日本撒尿小便嘘嘘汇集6| 午夜免费男女啪啪视频观看 | 免费看a级黄色片| 亚洲性夜色夜夜综合| 韩国av在线不卡| 成人精品一区二区免费| 免费在线观看影片大全网站| 哪里可以看免费的av片| 日韩欧美精品免费久久| 欧洲精品卡2卡3卡4卡5卡区| 中国美白少妇内射xxxbb| 亚洲精品在线观看二区| 中文字幕熟女人妻在线| 日韩亚洲欧美综合| 欧美一级a爱片免费观看看| 亚洲中文字幕日韩| 免费av毛片视频| 99国产精品一区二区蜜桃av| 女人被狂操c到高潮| 熟妇人妻久久中文字幕3abv| 日韩一本色道免费dvd| 亚洲无线在线观看| 亚洲人与动物交配视频| 男女那种视频在线观看| 欧美区成人在线视频| 88av欧美| 亚洲国产色片| 国产激情偷乱视频一区二区| 91久久精品电影网| 狂野欧美激情性xxxx在线观看| 午夜激情福利司机影院| 啦啦啦韩国在线观看视频| 欧美三级亚洲精品| 精品午夜福利在线看| 国产精品嫩草影院av在线观看 | 亚洲精品在线观看二区| 日韩强制内射视频| 久久久国产成人精品二区| 久久久久免费精品人妻一区二区| 午夜福利18| 嫩草影院精品99| 久久久久性生活片| av在线亚洲专区| 黄片wwwwww| 99久久精品国产国产毛片| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 永久网站在线| 国产一区二区三区在线臀色熟女| 日本黄色片子视频| 国产精品1区2区在线观看.| 国产精品免费一区二区三区在线| 美女cb高潮喷水在线观看| 午夜免费男女啪啪视频观看 | 国产蜜桃级精品一区二区三区| 午夜福利视频1000在线观看| 啦啦啦观看免费观看视频高清| 国产免费一级a男人的天堂| 国内揄拍国产精品人妻在线| videossex国产| 国产激情偷乱视频一区二区| 午夜免费成人在线视频| 亚洲国产精品成人综合色| 黄色丝袜av网址大全| 午夜福利在线观看免费完整高清在 | 极品教师在线视频| 全区人妻精品视频| 亚洲熟妇熟女久久| 国产v大片淫在线免费观看| 一个人看的www免费观看视频| 长腿黑丝高跟| 在线看三级毛片| 色综合色国产| 十八禁网站免费在线| 桃色一区二区三区在线观看| 一个人看视频在线观看www免费| 亚洲综合色惰| 最近在线观看免费完整版| 日日摸夜夜添夜夜添小说| 3wmmmm亚洲av在线观看| 日本一二三区视频观看| 亚洲av成人精品一区久久| 国产乱人视频| 岛国在线免费视频观看| 欧美国产日韩亚洲一区| 亚洲中文字幕一区二区三区有码在线看| 免费高清视频大片| 一进一出抽搐动态| 国产午夜精品论理片| 欧美日韩精品成人综合77777| 美女cb高潮喷水在线观看| 亚洲中文日韩欧美视频| 3wmmmm亚洲av在线观看| 给我免费播放毛片高清在线观看| 免费一级毛片在线播放高清视频| 精品久久久久久,| 成人无遮挡网站| 真人做人爱边吃奶动态| 欧美日韩精品成人综合77777| 日韩亚洲欧美综合| av在线亚洲专区| 亚洲一级一片aⅴ在线观看| 亚洲人与动物交配视频| 搡老熟女国产l中国老女人| 久久久成人免费电影| 午夜激情欧美在线| 国产黄片美女视频| 亚洲一级一片aⅴ在线观看| 午夜激情福利司机影院| 亚洲va日本ⅴa欧美va伊人久久| 久久中文看片网| 欧美丝袜亚洲另类 | 99九九线精品视频在线观看视频| 亚洲不卡免费看| 人妻夜夜爽99麻豆av| 亚洲综合色惰| 国产精品一区www在线观看 | 91在线精品国自产拍蜜月| 国产在线精品亚洲第一网站| 亚洲成人久久性| 丰满乱子伦码专区| 亚洲人成网站在线播| 日韩欧美 国产精品| 国产精品精品国产色婷婷| 如何舔出高潮| 日本-黄色视频高清免费观看| 琪琪午夜伦伦电影理论片6080| 久久精品国产99精品国产亚洲性色| h日本视频在线播放| 欧美日韩中文字幕国产精品一区二区三区| 免费看日本二区| 国产极品精品免费视频能看的| 日本爱情动作片www.在线观看 | a级一级毛片免费在线观看| 午夜a级毛片| 最近在线观看免费完整版| 亚洲精华国产精华液的使用体验 | 欧美最黄视频在线播放免费| 国产真实伦视频高清在线观看 | 永久网站在线| 嫩草影院新地址| 亚洲av中文字字幕乱码综合| 欧美日韩精品成人综合77777| 欧美色视频一区免费| 国内精品一区二区在线观看| 99久久精品热视频| 狂野欧美白嫩少妇大欣赏| 亚洲色图av天堂| 噜噜噜噜噜久久久久久91| 又爽又黄a免费视频| 此物有八面人人有两片| 色播亚洲综合网| 日日摸夜夜添夜夜添小说| 国内揄拍国产精品人妻在线| 99精品久久久久人妻精品| 在线a可以看的网站| 国产精品爽爽va在线观看网站| 国产久久久一区二区三区| 欧美一区二区精品小视频在线| av视频在线观看入口| 日韩大尺度精品在线看网址| 日本与韩国留学比较| 国产精华一区二区三区| 搞女人的毛片| 变态另类丝袜制服| 午夜福利成人在线免费观看| 精品一区二区三区视频在线观看免费| 超碰av人人做人人爽久久| 美女被艹到高潮喷水动态| 高清日韩中文字幕在线| 蜜桃亚洲精品一区二区三区| 国产午夜福利久久久久久| 日韩欧美免费精品| 99热精品在线国产| 香蕉av资源在线| xxxwww97欧美| 不卡一级毛片| 国产精品久久久久久av不卡| 国产午夜精品论理片| 亚洲一级一片aⅴ在线观看| 99久国产av精品| 最好的美女福利视频网| 久久午夜亚洲精品久久| 免费无遮挡裸体视频| 综合色av麻豆| 亚洲在线自拍视频| 国产精品av视频在线免费观看| 国产美女午夜福利| 欧美三级亚洲精品| 亚洲七黄色美女视频| 窝窝影院91人妻| 亚州av有码| 国产精品永久免费网站| www日本黄色视频网| 午夜爱爱视频在线播放| av福利片在线观看| 欧美性感艳星| 可以在线观看毛片的网站| 精品一区二区免费观看| 久久草成人影院| 婷婷精品国产亚洲av在线| 久久九九热精品免费| 免费无遮挡裸体视频| 亚洲精品国产成人久久av| 老司机午夜福利在线观看视频| 精品人妻1区二区| 毛片女人毛片| 少妇高潮的动态图| 美女xxoo啪啪120秒动态图| 亚洲中文日韩欧美视频| 色哟哟·www| 久久久久久久亚洲中文字幕| 日韩精品有码人妻一区| 欧美丝袜亚洲另类 | 美女cb高潮喷水在线观看| 看免费成人av毛片| 国产精品女同一区二区软件 | 俄罗斯特黄特色一大片| 国产精品电影一区二区三区| av在线蜜桃| 欧美成人性av电影在线观看| 国产探花在线观看一区二区| 久久亚洲精品不卡| 性欧美人与动物交配| 欧美日韩黄片免| 俺也久久电影网| 不卡视频在线观看欧美| 欧美性猛交黑人性爽| 精品一区二区三区av网在线观看| 97人妻精品一区二区三区麻豆| 欧美色视频一区免费| 国产视频一区二区在线看| 69av精品久久久久久| 成人永久免费在线观看视频| 国产精品爽爽va在线观看网站| 国产色爽女视频免费观看| 久久久久久伊人网av| 亚洲avbb在线观看| 日本色播在线视频| 久久久久性生活片| av在线蜜桃| 亚洲欧美日韩高清专用| 亚洲经典国产精华液单| 亚洲七黄色美女视频| av在线亚洲专区| 听说在线观看完整版免费高清| 大又大粗又爽又黄少妇毛片口| 久久中文看片网| 国产又黄又爽又无遮挡在线| 少妇人妻一区二区三区视频| 国产精品一区www在线观看 | 久久久色成人| 国产极品精品免费视频能看的| 男人的好看免费观看在线视频| 国产高清不卡午夜福利| 嫩草影院入口| 少妇高潮的动态图| 欧美成人性av电影在线观看| 亚洲午夜理论影院| 成人特级av手机在线观看| 国产激情偷乱视频一区二区| 深夜精品福利| 国产又黄又爽又无遮挡在线| av中文乱码字幕在线| 很黄的视频免费| 又粗又爽又猛毛片免费看| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| 99riav亚洲国产免费| 国产在线精品亚洲第一网站| 91麻豆精品激情在线观看国产| 老司机福利观看| 午夜影院日韩av| 欧美中文日本在线观看视频| 午夜免费成人在线视频| 免费在线观看成人毛片| 我的老师免费观看完整版| 日韩欧美精品免费久久| 日本成人三级电影网站| 一个人观看的视频www高清免费观看| 婷婷亚洲欧美| 国产黄色小视频在线观看| 国产亚洲精品av在线| 欧美性猛交黑人性爽| 日韩一区二区视频免费看| 久久精品综合一区二区三区| 国产探花在线观看一区二区| 亚洲av免费在线观看| 亚洲三级黄色毛片| 国产一区二区在线av高清观看| 国产又黄又爽又无遮挡在线| 香蕉av资源在线| 黄色丝袜av网址大全| 色吧在线观看| 成人综合一区亚洲| 欧美一区二区国产精品久久精品| 亚洲中文日韩欧美视频| 免费看美女性在线毛片视频| 成人永久免费在线观看视频| 丝袜美腿在线中文| 久久热精品热| 欧美国产日韩亚洲一区| 亚洲国产欧美人成| 国产精品国产高清国产av| 国内精品宾馆在线| www日本黄色视频网| 亚洲人成网站高清观看| 舔av片在线| 色5月婷婷丁香| 男插女下体视频免费在线播放| 神马国产精品三级电影在线观看| 国产在视频线在精品| 99九九线精品视频在线观看视频| 免费观看在线日韩| 99久国产av精品| 国产黄a三级三级三级人| 免费看光身美女| 午夜精品久久久久久毛片777| 午夜福利欧美成人| 赤兔流量卡办理| 最后的刺客免费高清国语| 亚洲乱码一区二区免费版| 中文字幕人妻熟人妻熟丝袜美| 99九九线精品视频在线观看视频| 午夜福利在线观看免费完整高清在 | 九色成人免费人妻av| 少妇熟女aⅴ在线视频| 午夜精品在线福利| 男人舔女人下体高潮全视频| 淫妇啪啪啪对白视频| 99国产精品一区二区蜜桃av| 日韩欧美 国产精品| 亚洲欧美精品综合久久99| 波多野结衣高清作品| 尾随美女入室| 日韩亚洲欧美综合| 日本 av在线| 成人一区二区视频在线观看| 69av精品久久久久久| 嫩草影视91久久| av黄色大香蕉| 日本爱情动作片www.在线观看 | 日韩强制内射视频| 日韩中文字幕欧美一区二区| 免费看日本二区| 国产亚洲91精品色在线| h日本视频在线播放| 国产日本99.免费观看| 91麻豆精品激情在线观看国产| netflix在线观看网站| 成年人黄色毛片网站| 日韩精品有码人妻一区| 欧美日本视频| 午夜免费激情av| 无遮挡黄片免费观看| 成人欧美大片| 少妇猛男粗大的猛烈进出视频 | 在线观看av片永久免费下载| 天堂影院成人在线观看| 国产一区二区三区在线臀色熟女| 欧美日韩综合久久久久久 | 男人舔奶头视频| 久久久久久久久久久丰满 | 看黄色毛片网站| 免费观看的影片在线观看| 男人舔女人下体高潮全视频| 免费一级毛片在线播放高清视频| 性插视频无遮挡在线免费观看| 国产三级在线视频| 亚洲av成人精品一区久久| а√天堂www在线а√下载| 国产女主播在线喷水免费视频网站 | 啦啦啦观看免费观看视频高清| 亚洲av熟女| 国产男靠女视频免费网站| 啦啦啦韩国在线观看视频| 毛片一级片免费看久久久久 |