一種新的圖譜域滾動軸承早期故障檢測與識別方法

陳子旭， 朱振杰， 盧國梁

(1. 山東大學(xué) 機械工程學(xué)院，濟南 250061； 2. 山東大學(xué) 高效潔凈機械制造教育部重點實驗室，濟南 250061；3. 山東大學(xué) 機械工程國家級實驗教學(xué)示范中心，濟南 250061)

滾動軸承是旋轉(zhuǎn)機械的關(guān)鍵部件，其損壞會直接影響設(shè)備的運行過程，導(dǎo)致設(shè)備停機、甚至引發(fā)人員傷亡。對滾動軸承的健康狀態(tài)進行實時、在線診斷具有重要意義[1]。

滾動軸承頻發(fā)的故障大多為局部缺陷形式的潛在故障，類型較復(fù)雜，其機械信號特點有：①強烈的非線性、非平穩(wěn)特性；②常伴有強噪聲干擾；②調(diào)制源弱導(dǎo)致早期故障特征通常很微弱。這些都給軸承故障在早期階段的檢測與識別帶來很大困難[2]。軸承早期故障檢測與識別的關(guān)鍵是從振動信號中提取故障特征信息。傳統(tǒng)的時域方法通過提取軸承振動信號的峰度、峭度、均方根值等統(tǒng)計指標，對信號進行動態(tài)描述，但其容易受到噪聲干擾，魯棒性不足[3]。以傅里葉變換為基礎(chǔ)的傳統(tǒng)頻域分析方法在處理非線性、非平穩(wěn)信振動號時，常常無法準確描述信號頻率隨時間的動態(tài)特征，制約了其在早期故障檢測中的應(yīng)用[4]。為了彌補上述方法的不足，時頻分析方法被相繼提出，代表性的時頻分析方法包括短時傅里葉變換、小波變換等。這些方法被認為可提供信號完整的時頻分布，因此已被廣泛應(yīng)用于軸承故障診斷領(lǐng)域[5]。隨著數(shù)字信號處理技術(shù)的不斷進步，近些年來，在傳統(tǒng)時頻分析方法基礎(chǔ)上，自適應(yīng)時頻分析方法得到了較大發(fā)展。Huang等[6]提出的經(jīng)驗?zāi)B(tài)分解(empirical modal decomposition，EMD)解決了傳統(tǒng)時頻分析存在的非自適應(yīng)問題，更加適合分析多分量信號，但其也存在端點效應(yīng)、模態(tài)混淆的缺點。Smith[7]提出的局部均值分解(local mean decomposition, LMD)是在EMD的基礎(chǔ)上，對端點效應(yīng)、模態(tài)混淆問題進行了改善，能獲得原始信號更加精確的時頻分布。盡管LMD方法在軸承故障診斷領(lǐng)域獲得了較多應(yīng)用，但其仍存在兩個問題：

(1) 由LMD分解得到多尺度分量信號數(shù)據(jù)量龐大、復(fù)雜，因此仍需進一步計算以提取更低維度的數(shù)據(jù)特征。Tian等[8]提出基于乘積函數(shù)(product functions，PF)分量構(gòu)建漢克爾矩陣，利用奇異值分解壓縮數(shù)據(jù)尺度以獲得更加穩(wěn)定的特征向量值；王洪明等[9]通過計算全矢能量熵從PF分量中提取有用信息，獲取得到了低維度的數(shù)據(jù)特征。

(2) 由LMD分解可以獲取信號的強數(shù)據(jù)特征，而軸承早期故障特征微弱，因此需要對LMD進行增強處理以獲取軸承早期微弱故障特征。最大相關(guān)峰度反褶積處理能夠有效消除噪聲對PF分量的影響，用于增強故障特征信息[10]；基于PF分量構(gòu)建Teager能量譜能夠有效識別滾動軸承故障信號中的沖擊特性，從而能夠提取微弱故障特征[11]；王志堅等[12]提出了一種將掩膜信號法與LMD相結(jié)合的LMD-MS(local mean decomposition-mask signal)方法，成功提取了軸承的微弱故障特征。

隨著圖論的不斷發(fā)展，圖模型作為一種高效的數(shù)學(xué)模型，已被應(yīng)用于機械設(shè)備狀態(tài)監(jiān)測領(lǐng)域[13]。Lu等[14]采用圖模型的方式進行數(shù)據(jù)建模，通過對模型進行距離度量實現(xiàn)對旋轉(zhuǎn)機械的動態(tài)描述。Wang等[15]首先采用無向加權(quán)圖的建模方式基于信號的周期譜進行建模，通過計算模型中位圖和距離指標來描述機械設(shè)備的動態(tài)狀態(tài)。Wang等[16]在之前工作的基礎(chǔ)上進一步擴展了圖建模策略，構(gòu)建了時空圖模型并通過比較試驗驗證了該方法的有效性。受以上研究的啟發(fā)，本文提出了一種新的圖譜域滾動軸承早期故障檢測與識別方法。本方法實現(xiàn)了兩方面改進：①建立的異常度動態(tài)指標實現(xiàn)了低維度軸承健康狀態(tài)動態(tài)描述；②基于PF分量構(gòu)建圖模型進行圖譜分析，實現(xiàn)了對LMD的增強處理。通過上述改進，本文所提方法可有效提取軸承早期故障的微弱特征。

1 方法簡介

方法主要包含以下3個步驟，如圖1所示。

圖1 方法流程圖

步驟1滾動軸承振動信號動態(tài)建模。采用LMD對滾動軸承振動信號進行分解，將分解得到的多尺度分量信號作為輸入構(gòu)建圖模型，計算相鄰圖模型的相似性指標以得到能描述滾動軸承健康狀態(tài)的動態(tài)指標。

步驟2早期故障檢測?；谏鲜鰳?gòu)建得到的動態(tài)指標，基于假設(shè)檢驗的方法，采用拉依達準則對軸承早期故障進行決策。

步驟3故障識別。對于檢測到的故障異常，進一步采用模式識別技術(shù)進行故障類型的識別。

2 滾動軸承健康狀態(tài)動態(tài)建模

滾動軸承由于復(fù)雜、惡劣的工況導(dǎo)致其振動信號經(jīng)常受到噪聲等冗余信息的干擾，并且不利于微弱故障特征的提取。LMD可有效將噪聲等冗余信息同有效信息隔離，進行多尺度分析。同時，隨著圖論的發(fā)展，圖建模方式已經(jīng)被證實在提取微弱故障信息方面的重要作用。因此，針對以上方法的優(yōu)勢，本文基于LMD和圖論方法，對滾動軸承健康狀態(tài)進行動態(tài)建模。

2.1 局部均值分解

作為一種自適應(yīng)分解方法，LMD能將信號分解為一系列具有物理意義的分量信號，反映信號能量在多尺度空間下的分布情況。針對監(jiān)測獲取的滾動軸承振動信號，記為x(t)，其經(jīng)過LMD分解的步驟如下。

步驟1根據(jù)所有的局部極值點ni，求出所有局部極值平均值mi和包絡(luò)估計值ai，用滑動平均法處理后得到局部均值函數(shù)m11(t)和a11(t)。

步驟3將上述迭代過程的包絡(luò)函數(shù)相乘得到包絡(luò)信號，即a1(t)=a11(t)·a12(t)…a1n(t)。

步驟4將純調(diào)頻信號s1n(t)與包絡(luò)信號a1(t)相乘，得到第一個PF分量，即FP1(t)=s1n(t)·a1(t)。

步驟5將FP1(t)從x(t)中分離出來，得到u1(t)，將其作為振動信號重復(fù)以上計算過程得到u2(t)，直到uq(t)為單調(diào)函數(shù)時停止迭代。

步驟6滾動軸承振動信號經(jīng)過上述迭代過程，被分解為一系列PF分量和一個殘余分量。

(1)

本文首先對滾動軸承振動信號采用無重疊滑動窗的方式進行LMD分解，原始信號被分解為一系列PF分量，即x(t)→{FP1(t),FP2(t),…,uq(t)}。從LMD定義來看，殘余分量uq(t)為一單調(diào)函數(shù)，其包含的振動信號極其微弱，具有的能量接近于零；并且振動信號分解后的前5個分量幾乎包含了滾動軸承振動信號的所有頻率成分，其能量和與信號的總能量近似相等[17]。因此，本文選取前5個PF分量用于圖建模。

2.2 圖建模

圖模型通常由一系列的節(jié)點V和邊L組成，即G={V,L}。通過對每對節(jié)點vi和vj之間計算距離，比如歐幾里得距離、曼哈頓距離、赫凌格距離等，以表示兩節(jié)點所構(gòu)成的邊lij的權(quán)重。圖模型可以有效的表示相鄰節(jié)點以及互連節(jié)點之間的連接信息，并且應(yīng)用到信號的多尺度分析中。本文將分解得到的前5個PF分量作為圖模型的節(jié)點序列，并采用最常用的歐氏距離度量方式，賦予圖模型每條邊權(quán)重進行圖建模，如圖2所示。本文以t1時刻為例，給出了圖建模過程示意圖，其建模過程分為以下幾個步驟。

圖2 圖建模示意圖

步驟1將分量FP1～FP5作為圖模型的節(jié)點序列，5個分量在t1時刻的取值分別為{FP1,1,FP2,1,FP3,1,FP4,1,FP5,1}，作為圖模型的節(jié)點{v1,v2,v3,v4,v5}，每對節(jié)點vi和vj之間的連接作為圖模型的邊lij，其中i,j∈[1,5],Z。

步驟2計算每對節(jié)點vi和vj之間的歐式距離dij，作為邊lij的權(quán)重對其進行賦值。

di,j=|vi-vj|

(2)

步驟3將計算得到的歐式距離dij如式(3)所示，以鄰接矩陣的形式表示，得到圖模型的矩陣表示。

(3)

使用來自每個PF分量的多尺度信息在每個數(shù)據(jù)點上構(gòu)建圖模型。因此，每個滑動窗內(nèi)部的振動信號則可以表示為一系列的圖模型，即{G1,G2,…,Gl}，l為滑動窗的長度，本文將l設(shè)為5 000。構(gòu)建的圖模型除了可以描述時序特性，還包含信號不同尺度之間的相關(guān)信息。

在圖論中，均值圖是用于度量一組圖模型平均水平的常見方式。對于在每個滑動窗內(nèi)構(gòu)建得到的圖模型序列，計算其均值圖用以對該段數(shù)據(jù)狀態(tài)進行描述，計算公式如下

(4)

2.3 相似性度量

構(gòu)建得到的圖模型具有較高的維度，同時，為了更好地描述滾動軸承動態(tài)特性、便于故障檢測的決策，引入了相似性度量指標。本文認為LMD分解得到的各分量對于提取故障信息具有相同的貢獻度，因此，采用歐式距離度量方式計算相鄰圖模型之間相似性，得到異常度分數(shù)序列，對滾動軸承運行過程中狀態(tài)變化的可能性進行量化。計算公式如下所示

(5)

3 早期故障檢測

當滾動軸承正常運轉(zhuǎn)時，其振動信號分布接近正態(tài)分布，此狀態(tài)下圖模型結(jié)構(gòu)穩(wěn)定，異常度分數(shù)服從高斯分布；而當滾動軸承在某一時刻發(fā)生故障時，隨之而來的是故障特征頻率的出現(xiàn)和沖擊成分的增多，其振動信號將偏離正態(tài)分布?；诖耍愿咚狗植紴榧僭O(shè)，采用拉依達準則對異常度分數(shù)進行假設(shè)檢驗(本文采用6σ準則)，實現(xiàn)早期故障時刻的定位。原假設(shè)和備擇假設(shè)如下所示

(6)

式中，μh-1和σh-1分別為歷史數(shù)據(jù){s1,s2,…,sh-1}的平均值和標準差，計算公式如下

(7)

(8)

假設(shè)H0設(shè)定：當前時刻數(shù)據(jù)與歷史數(shù)據(jù)服從同一個正態(tài)分布，即滾動軸承處于正常運行狀態(tài)。備擇假設(shè)H1設(shè)定：當前時刻數(shù)據(jù)出現(xiàn)了異于正常的狀態(tài)變化，該時刻即為檢測到的由正常轉(zhuǎn)向故障的狀態(tài)變化點。實線為異常度分數(shù)曲線，虛線為上下6σ控制限，圈表示檢測到的狀態(tài)變化點，如圖3所示。當滾動軸承運行出現(xiàn)了異于正常的狀態(tài)變化時，圖模型結(jié)構(gòu)的變化會導(dǎo)致異常度分數(shù)發(fā)生較大變化，異常度分數(shù)不再服從先前分布，通過判斷當前時刻異常度分數(shù)值是否處于此前分布的置信區(qū)間，便可以確定當前時刻滾動軸承是否發(fā)生狀態(tài)變化。

圖3 早期故障檢測示意圖

4 故障識別

基于第3章早期故障檢測，一旦檢測到滾動軸承運行狀態(tài)發(fā)生異常，可利用模式識別方法(如支持向量機(support vector machines, SVM)、人工神經(jīng)網(wǎng)絡(luò)、隨機森林、深度學(xué)習(xí)等)實現(xiàn)故障類型的識別。在常見的識別方法中，K最近鄰(k-nearest neighbors, KNN) 算法被廣泛應(yīng)用。KNN算法通過測量不同特征值之間的距離，依照多數(shù)表決的決策規(guī)則進行分類，即由待測樣本的K個最近鄰的訓(xùn)練實例中的多數(shù)類別決定待測樣本的類別。本文給出了KNN算法用于分類的示意圖，如圖4所示。首先計算測試數(shù)據(jù)與各訓(xùn)練數(shù)據(jù)之間的距離，按照距離的遞增關(guān)系進行排序，如果一個樣本在特征空間中的K個最相似(即特征空間中最近鄰的前K個點)的樣本中的大多數(shù)屬于某一個類別，則該樣本也屬于這個類別。KNN是一種非參的、惰性算法模型，其簡單易用，對異常值不敏感，可用于非線性分類，同時具有較低的運算復(fù)雜度，無需大量數(shù)據(jù)進行網(wǎng)絡(luò)訓(xùn)練[18]。

圖4 KNN算法原理圖

本文采用KNN開展故障類型識別的試驗。需要特別指出的是，本文主要關(guān)注點在于提出了一種基于圖的振動信號建模與故障表征方法，因此未對其他模式識別方法進行驗證。實際上，本文方法也可與其他模式識別方法進行結(jié)合，以進一步提高其準確率。一旦檢測到滾動軸承運行狀態(tài)發(fā)生異常，則狀態(tài)變化點之后的數(shù)據(jù)為待識別數(shù)據(jù)，將待測數(shù)據(jù)的均值圖結(jié)合KNN算法對故障類型進行識別。值得提出的是，均值圖表示為鄰接矩陣的形式，在輸入KNN之前需將其轉(zhuǎn)換為行向量。

5 試驗及討論

5.1 試驗數(shù)據(jù)

分別在XJTU-SY[19]和凱斯西儲大學(xué)(Case Western Reserve University, CWRU)[20]數(shù)據(jù)集上進行了驗證。其中，XJTU-SY數(shù)據(jù)集用于驗證早期故障檢測，該數(shù)據(jù)集是由西安交通大學(xué)設(shè)計科學(xué)研究所和基礎(chǔ)構(gòu)件研究所以及長興晟陽科技有限公司提供的滾動軸承加速退化試驗數(shù)據(jù)。該試驗數(shù)據(jù)在3種工況下取得：①轉(zhuǎn)速2 100 r/min，徑向力12 kN；②轉(zhuǎn)速2 250 r/min，徑向力11 kN；③轉(zhuǎn)速2 400 r/min，徑向力10 kN，每個工況下有5個軸承。數(shù)據(jù)采樣頻率為25.6 kHz，采樣間隔為1 min，每個采樣持續(xù)時間為1.28 s，失效原因包括外圈磨損、內(nèi)圈磨損、滾動體損壞、保持架斷裂等。每一個測試軸承的詳細信息，如表1所示。有關(guān)該試驗系統(tǒng)以及試驗詳情見文獻[19]。

表1 XJTU-SY軸承數(shù)據(jù)集信息

CWRU數(shù)據(jù)集用于驗證故障類型識別，CWRU數(shù)據(jù)集來自凱斯西儲大學(xué)軸承數(shù)據(jù)中心。數(shù)據(jù)集包含正常數(shù)據(jù)、采樣率12 kHz的驅(qū)動端故障數(shù)據(jù)、采樣率48 kHz的驅(qū)動端故障數(shù)據(jù)、風(fēng)扇端故障數(shù)據(jù)。本文采用48 kHz驅(qū)動端故障數(shù)據(jù)進行試驗，試驗數(shù)據(jù)在4種工況下獲得：①負載0，轉(zhuǎn)速1 797 r/min；②負載745.69 W，轉(zhuǎn)速1 772 r/min；③負載1 491.39 W，轉(zhuǎn)速1 750 r/min；④負載2 237.07 W，轉(zhuǎn)速1 730 r/min。每種工況下包含滾動軸承4種工作狀態(tài)：①正常(Normal)；②外圈磨損(OF)；③內(nèi)圈磨損(IF)；④滾動體磨損(RF)，其中，3種故障類型均是通過放電加工產(chǎn)生，每種故障類型包括3種不同的破壞尺寸，如表2所示。數(shù)據(jù)集共包括10種不同類型的數(shù)據(jù)。有關(guān)該數(shù)據(jù)集更為詳細的介紹見文獻[20]。

表2 CWRU 驅(qū)動端故障數(shù)據(jù)類型

5.2 故障檢測試驗

XJTU-SY數(shù)據(jù)集既包含緩慢失效的故障信號，又包含突然失效的故障信號。本文方法在故障檢測階段的兩個例子，如圖5、圖6所示。圖5(a)、圖6(a)為XJTU-SY數(shù)據(jù)集中數(shù)據(jù)Bearing1_3和Bearing3_3，分別對應(yīng)這兩種不同退化形式的信號。圖5(b)、圖6(b)為本文檢測結(jié)果，其中圈為檢測得到的狀態(tài)變化點。若故障發(fā)生后，方法發(fā)出預(yù)警，認為方法成功檢測到狀態(tài)變化，并且預(yù)警時刻與真正開始發(fā)生故障時刻的延遲越小則檢測效果越好；若方法發(fā)出預(yù)警，但軸承此時實際并未發(fā)生故障，則認為方法此時為誤報。對于突然失效的故障信號，其時域波形發(fā)生較大突變，通過肉眼便能確定故障開始時刻[21]；而對于緩慢失效的故障信號，通常對應(yīng)軸承退化的多個階段，很難確定故障開始時刻[22]。因此針對緩慢失效故障信號，對檢測到的狀態(tài)變化點所在滑動窗數(shù)據(jù)時域波形和頻譜分析，并與正常數(shù)據(jù)進行比較，以確定方法檢測的準確性。本文以Bearing1_3數(shù)據(jù)為例進行說明，該軸承在試驗時發(fā)生了軸承外圈故障，正常數(shù)據(jù)和檢測到的故障數(shù)據(jù)的時域波形和頻譜圖，如圖7、圖8所示。從時域來看，檢測到的故障數(shù)據(jù)與正常數(shù)據(jù)相比波形上有較大差距，出現(xiàn)了較多故障發(fā)生時典型的調(diào)幅和沖擊成分。該數(shù)據(jù)集軸承一倍轉(zhuǎn)頻理論值為34.38 Hz，從頻譜圖來看，正常數(shù)據(jù)中主要包含35.85 Hz，102.40 Hz，291.90 Hz，327.70 Hz頻率成分，分別對應(yīng)軸承轉(zhuǎn)動的基頻和倍頻；檢測到的故障數(shù)據(jù)頻譜中，除了一倍轉(zhuǎn)頻成分35.85 Hz外，還出現(xiàn)了107.50 Hz及其二倍頻215.10 Hz、三倍頻322.60 Hz等倍頻成分，該頻率成分接近外圈故障特征頻率理論值107.91 Hz，由此可以獲知本文方法成功做出狀態(tài)變化預(yù)警。

圖5 Bearing1_3檢測結(jié)果

圖6 Bearing3_3檢測結(jié)果

圖7 正常數(shù)據(jù)波形及頻譜

圖8 故障數(shù)據(jù)波形及頻譜

為了進一步驗證方法的優(yōu)勢，本文將所提方法與廣泛采用的均方根值(root mean square, RMS)、峭度(Kurtosis)、偏度(Skewness)、均方誤差(mean square error, MSE)、標準差(standard deviation, SD)、方差(variation, Var)、特征能量(signature energy, SE)、特征能量比(signature energy ratio, SER)共8種指標進行了比較。將每種方法首次檢測到的異常警報作為軸承早期故障，利用查準率、查全率和綜合測度(F1_score)3個指標對每個方法進行綜合評價，計算公式如下

(9)

(10)

(11)

所有方法的評價結(jié)果，如表3所示。其中，本文方法同SE、SER、MSE、Var指標均達到了查全率100%，查準率100%，綜合測度1.0，優(yōu)于RMS、Kurtosis、SD、Skewness指標。進一步地，本文在表4給出了本文方法同SE、SER、MSE、Var 4種方法預(yù)警時刻的對比。以表4中Bearing1_3為例，本文方法在第386個滑動窗處發(fā)出預(yù)警，早于其他方法預(yù)警時刻。從表4可知，除Bearing1_2和Bearing3_2，本文方法發(fā)出故障預(yù)警的時刻均早于其他方法，說明本文方法能夠更有效提取滾動軸承振動信號中的早期微弱故障信息。

表3 試驗對比結(jié)果

表4 不同方法的預(yù)警時刻

在實際應(yīng)用中，計算時間同樣對早期故障檢測有著重要影響[23]。本文以“Bearing1_3”為測試數(shù)據(jù)，在MATLAB R2017b，i5-5200U CPU 2.2 GHz，RAM 4.00 GB操作環(huán)境下，對不同方法計算時間進行測試。測試數(shù)據(jù)前386個滑動窗內(nèi)數(shù)據(jù)每個滑動窗內(nèi)的平均計算時間，如表5所示。幾種方法前386個滑動窗數(shù)據(jù)總的計算時間，如圖9所示。相比4種對比方法，本文方法計算時間較高，并且隨著數(shù)據(jù)量增加升高明顯。但是，值得提出的是，本文方法在第386個滑動窗處已成功發(fā)出預(yù)警，而其他方法在較晚時刻，甚至是完全失效階段才發(fā)出預(yù)警。相比其他方法，本文方法對早期故障階段的信息更為敏感，能在故障發(fā)生后的較早階段發(fā)出故障預(yù)警，并且對于緩慢失效過程尤為明顯。因此，本文通過LMD和圖論方法提取信號多尺度信息，相比其他方法能夠更有效提取早期微弱故障特征。在實際應(yīng)用中，在設(shè)備故障的早期階段發(fā)出預(yù)警可幫助操作人員及時排查設(shè)備故障原因、避免損失進一步擴大，為停機操作和制定故障維修計劃爭取更多的時間。作為未來工作方向之一，我們將繼續(xù)改進所提方法，例如采用并行計算等方式提高方法計算效率，使其更加適用于實際應(yīng)用。

表5 平均計算時間

圖9 測試數(shù)據(jù)計算時間

5.3 故障識別試驗

首先利用t-SNE方法將故障數(shù)據(jù)映射到低維空間進行三維可視化。不同形狀的散點表示不同類型的軸承數(shù)據(jù)，如圖10所示。如IF7為軸承內(nèi)圈故障，故障尺寸為0.177 8 mm，其他具體圖例含義同表2。由圖10可知，相同類型的數(shù)據(jù)具有較好的聚集性，并且不同類型數(shù)據(jù)之間具有較好的分離性。這說明所建圖模型可以有效的對不同類數(shù)據(jù)進行描述，便于分類器進行分類實現(xiàn)故障診斷。

圖10 故障數(shù)據(jù)在特征空間的三維可視化

本文將構(gòu)建的圖模型的50%用于訓(xùn)練，50%用于測試。在數(shù)據(jù)集給出的4種負載條件下，對這10種類型數(shù)據(jù)利用KNN算法進行分類識別。本文方法在4種負載條件下，分別能達到90%、100%、100%、100%的準確率，為了進一步說明方法的優(yōu)越性，本文同現(xiàn)有的方法進行了對比，本文方法和對比方法的試驗結(jié)果，如表6所示。從表6可知，對于在負載745.69 W、1 491.39 W、2 237.07 W條件下，本文方法的表現(xiàn)為最優(yōu)。在0條件下，楊云博等[24-26]準確率優(yōu)于本文方法。楊云博等利用一種改進的人工魚群算法對SVM進行優(yōu)化，通過提取到的振動信號時域特征進行故障診斷；王太勇等提出了一種基于DenseNet的卷積核dropout(KD)智能故障診斷模型KD-DenseNet。其中，楊云博等只考慮了正常、外圈磨損、內(nèi)圈磨損、滾動體磨損4種故障類型，并且并未將故障破壞尺寸考慮在內(nèi)；王太勇等試驗中的故障類型只有7種，沒有考慮故障破壞尺寸為0.533 4 mm時的故障類型。本文認為故障類型數(shù)目的不同是導(dǎo)致這兩種方法準確率均高于所提方法的原因。并且，需要指出，在實際應(yīng)用中，軸承的故障類型多種多樣，不同的工況會導(dǎo)致不同的故障破壞尺寸，本文方法將CWRU數(shù)據(jù)集中的10種故障類型全部考慮在內(nèi)，分類更加精確，具有更好的現(xiàn)實意義。Zeng等提出了一種基于柔性凸殼的最大余量分類方法，并將其應(yīng)用于滾動軸承的故障診斷，該方法將83%的數(shù)據(jù)用于訓(xùn)練，17%的數(shù)據(jù)用于測試，在負載0時，其準確率達到99.60%，并且故障種類數(shù)為10。相比該方法，本文方法使用了更少的訓(xùn)練數(shù)據(jù)，具有更高的實際應(yīng)用價值。在實際應(yīng)用過程中，軸承故障類型多種多樣，軸承故障類數(shù)據(jù)又通常難以大量獲取，因此，使用較少的故障數(shù)據(jù)完成故障診斷任務(wù)更加符合實際應(yīng)用需求。因此，相比其他方法，本文所提方法具有較大優(yōu)勢。

表6 故障識別結(jié)果

6 結(jié) 論

針對滾動軸承健康狀態(tài)監(jiān)測問題，本文提出了一種新的圖譜域滾動軸承早期故障檢測與識別方法。通過基于PF分量構(gòu)建圖模型進行圖譜分析，能夠有效獲取滾動軸承早期故障特征，通過計算相鄰圖模型間的相似性指標對軸承健康狀態(tài)進行低維度動態(tài)描述，實現(xiàn)了滾動軸承早期故障檢測。利用LMD對滾動軸承振動信號進行多尺度分析，通過圖建模的方式對得到的多尺度數(shù)據(jù)增強，能夠很好地用于故障類型描述，達到滾動軸承故障識別的目標。在后續(xù)工作中，將研究提高方法的運算效率，使其更好的應(yīng)用于實際場合中。