朱玉鵬, 晉會杰，謝鵬飛, 牛青波，楊茹萍

(1.河南科技大學(xué) 土木工程學(xué)院,河南 洛陽 471003;2.商丘工學(xué)院，河南 商丘 476000；3.洛陽軸承研究所有限公司,河南 洛陽 471039)

在航空航天、汽車、工程機(jī)械領(lǐng)域，軸承常處于變速、變載等典型工況條件或極端環(huán)境，運(yùn)行過程中滾動體與保持架的碰撞不可避免[1-5]，系統(tǒng)表現(xiàn)出復(fù)雜的非線性動力學(xué)特性[6-7]，保持架成為軸承故障高發(fā)部位。目前，針對保持架的研究大多集中在保持架數(shù)值仿真模型的動力學(xué)特性研究領(lǐng)域[8-9]，關(guān)于保持架狀態(tài)識別及故障診斷研究的相關(guān)工作很少。

陀螺電機(jī)軸承主要元件由高強(qiáng)度合金鋼制成，保持架則采用非金屬材料，其剛度和質(zhì)量均遠(yuǎn)小于內(nèi)、外圈，導(dǎo)致軸承系統(tǒng)中保持架故障特征頻率振動量級小，常常堙沒在復(fù)雜噪聲背景下，傳統(tǒng)信號分析方法很難識別。且保持架一旦發(fā)生故障，振動傳遞路徑復(fù)雜，測試采集到的信號大多是多分量調(diào)幅調(diào)頻信號，如何在強(qiáng)噪聲、復(fù)雜調(diào)制背景下識別保持架頻率是保持架故障診斷的關(guān)鍵。

HHT(Hilbert-Huang Transform)具有良好的時(shí)頻特性，且其自適應(yīng)的分解模式非常適合非線性和非平穩(wěn)信號的分析[10-11]，具有很高的信噪比，在處理調(diào)制信號方面效果良好[12-14]。因此，采用基于自適應(yīng)噪聲完備集合經(jīng)驗(yàn)?zāi)B(tài)分解(Complete Ensemble Empirical Mode Decomposition With Adaptive Noise，CEEMDAN)算法對HHT進(jìn)行改進(jìn)，并與自動提取敏感IMF分量算法相結(jié)合，通過對實(shí)測保持架故障軸承振動信號的特征提取，驗(yàn)證該方法的有效性。

1 改進(jìn)的HHT方法

1.1 CEEMDAN

CEEMDAN是在經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)[15-16]的每個(gè)階段有限次添加自適應(yīng)白噪聲，通過計(jì)算唯一余量信號獲取固有模態(tài)函數(shù)(IMF)，在消除模態(tài)混疊[17-18]的同時(shí)降低了噪聲殘留。

假設(shè)Ek(·)為EMD處理后得到的第k個(gè)IMF，且令w(i)為一個(gè)零均值單位方差的白噪聲，則CEEMDAN算法的具體步驟為：

1)采用EMD算法對x(i)=x+β0w(i)進(jìn)行分解，得到第1個(gè)CEEMDAN的IMF，即

2)計(jì)算第1個(gè)殘差：r1=x-f1。

3)采用EMD計(jì)算r1+β1E1(w(i))的第1個(gè)模態(tài)，并定義為CEEMDAN的第2個(gè)IMF，即

4)當(dāng)k=1,…,K，計(jì)算第k個(gè)殘差rk=xk-1-fk。

5)采用EMD計(jì)算rk+βkEk(w(i))的第1個(gè)模態(tài)，并定義為CEEMDAN的第k+1個(gè)IMF，即

6) 當(dāng)k=k+1時(shí)，返回至第4步，計(jì)算CEEMDAN的下一個(gè)IMF。

為驗(yàn)證CEEMDAN算法的有效性，對仿真信號進(jìn)行CEEMDAN處理，添加白噪聲50次，幅值為仿真信號標(biāo)準(zhǔn)差的0.01倍，最大篩選迭代次數(shù)為100次。分解結(jié)果如圖1所示，從圖中可以看出,脈沖信號和正弦信號清晰地分解到了不同的IMF中。

圖1 CEEMDAN分解IMF圖Fig.1 Diagram of IMFs by CEEMDAN

1.2 敏感IMF選擇算法

對信號進(jìn)行CEEMDAN處理后故障特征成分會分解到某一個(gè)或幾個(gè)敏感IMF中，信號特征提取的關(guān)鍵在于敏感IMF的選取。因此，提出一種基于相關(guān)系數(shù)與峭度指標(biāo)結(jié)合的敏感IMF選擇算法。

以各IMF與原始信號的相關(guān)系數(shù)作為統(tǒng)計(jì)指標(biāo)，相關(guān)系數(shù)越大，則IMF與原始信號相關(guān)程度越高，即為目標(biāo)敏感IMF。為使結(jié)果清晰且減小幅值較小的真實(shí)IMF分量被剔除的可能性，將所有IMF與原始信號進(jìn)行歸一化處理，各IMF與原始信號的歸一化相關(guān)系數(shù)ci(i=1,2,…,K)為

閾值Q定義為歸一化相關(guān)系數(shù)ci的標(biāo)準(zhǔn)差。若ci>Q，則保留相對應(yīng)的IMF；若ci≤Q，則相應(yīng)的IMF不作分析。

峭度為量綱一的參數(shù)，與結(jié)構(gòu)尺寸、載荷分布無關(guān)，但受脈沖信號影響較大，軸承表面出現(xiàn)損傷類故障時(shí)峭度指標(biāo)變化明顯[19]。軸承無故障運(yùn)轉(zhuǎn)時(shí)，時(shí)域信號波形平穩(wěn)無明顯尖銳波峰，峭度指標(biāo)值接近3；一旦軸承存在局部故障，由周期旋轉(zhuǎn)引起的故障部位脈沖信號持續(xù)出現(xiàn)，信號波形尖峰程度加劇，峭度值會隨之增大，當(dāng)峭度指標(biāo)值大于3時(shí)，表明軸承已出現(xiàn)故障，大于8時(shí)故障已十分明顯。設(shè)定峭度閾值為3，進(jìn)行敏感IMF篩選。

綜上，采用歸一化相關(guān)系數(shù)和峭度值相結(jié)合的算法來進(jìn)行敏感IMF提取，可以凸顯潛在的故障成分，削弱噪聲信號影響，從而實(shí)現(xiàn)故障特征全面、有效地提取。

1.3 改進(jìn)的HHT方法

基于以上分析，改進(jìn)的HHT分析方法的主要步驟為：1)采用CEEMDAN算法對采集的保持架故障信號進(jìn)行本征模態(tài)分解，得到一組IMF；2)分別計(jì)算各階IMF與原始信號的歸一化相關(guān)系數(shù)并計(jì)算其標(biāo)準(zhǔn)差，計(jì)算各階IMF峭度指標(biāo)，根據(jù)歸一化相關(guān)系數(shù)和峭度指標(biāo)綜合篩選敏感IMF；3)對篩選得到的IMF進(jìn)行Hilbert變換包絡(luò)并進(jìn)行頻譜分析；4)根據(jù)包絡(luò)譜對振動信號進(jìn)行特征提取和故障診斷。

2 陀螺電機(jī)軸承保持架故障實(shí)例分析

研究對象為某撓性支承陀螺儀中存在保持架故障的軸承，其具有結(jié)構(gòu)精密、體積小、轉(zhuǎn)速高等特點(diǎn)，且不可拆卸，測試過程中振動傳遞路徑復(fù)雜，信號噪聲較大。其中的結(jié)構(gòu)固有振動與軸承工作狀況無關(guān)，僅取決于系統(tǒng)的質(zhì)量、材質(zhì)、形狀等，屬于噪聲信號；由Hertz彈性接觸理論產(chǎn)生的振動信號為特征信號，其大小取決于軸承結(jié)構(gòu)尺寸以及工作轉(zhuǎn)速，能反映軸承的損傷部位，是進(jìn)行軸承故障診斷的重要依據(jù)。

該陀螺電機(jī)軸承為角接觸球軸承，具體參數(shù)見表1，保持架材料為多孔聚酰亞胺。由于結(jié)構(gòu)復(fù)雜，采用在特制底座上粘貼PCB加速度傳感器的方法進(jìn)行振動信號采集，采集工作由LMS Signature Testing系統(tǒng)完成。分析工作則由MATLAB軟件完成。測試轉(zhuǎn)速為18 000 r/min，采樣頻率25 600 Hz，頻率分辨率0.781 25 Hz，采樣時(shí)間0.64 s，計(jì)算可得保持架故障特征頻率fc為112.1 Hz。

表1 軸承結(jié)構(gòu)參數(shù)Tab.1 Structural parameters of bearing

故障軸承時(shí)域即FFT頻譜圖如圖2所示。時(shí)域信號波動較大、沖擊振動明顯。頻譜圖中，低頻段僅能發(fā)現(xiàn)300 Hz的成分，其為轉(zhuǎn)子轉(zhuǎn)動頻率；在1 000～8 000 Hz內(nèi)頻率分布密集，其中3 374 Hz接近保持架故障特征頻率30倍頻，但峰值相對較低，整體難以發(fā)現(xiàn)保持架低頻故障成分，無法得出準(zhǔn)確結(jié)論。

圖2 軸承時(shí)域信號及頻譜圖Fig.2 Time domain signal and frequency spectrum of bearing

為實(shí)現(xiàn)保持架故障特征的準(zhǔn)確提取，利用改進(jìn)的HHT算法對圖中時(shí)域信號進(jìn)行分析處理。CEEMDAN算法參數(shù)設(shè)置：白噪聲幅值為振動信號標(biāo)準(zhǔn)差的0.1倍，白噪聲添加598次，最大迭代次數(shù)500。CEEMDAN處理后共得到19個(gè)IMF，包含主要信息的前5階IMF如圖3所示。各IMF與原始信號的歸一化相關(guān)系數(shù)以及峭度指標(biāo)如圖4所示。

圖3 CEEMDAN處理后部分IMF圖Fig.3 Diagram part of IMFs after CEEMDAN

圖4 各IMF相關(guān)系數(shù)及峭度指標(biāo)Fig.4 Correlation coefficient and kurtosis index of all IMFs

根據(jù)自相關(guān)系數(shù)法，篩選得到敏感分量為IMF1，IMF2，IMF3和IMF4；根據(jù)峭度指標(biāo)則得到敏感分量為IMF1和IMF19。由于IMF19已接近殘差分量且?guī)缀醪话l率信息，故最終選擇IMF1為敏感分量。

對IMF1進(jìn)行Hilbert包絡(luò)分析，結(jié)果如圖5所示。從圖中可以看出，Hilbert包絡(luò)譜中最大峰值出現(xiàn)在112.5 Hz處，與保持架故障特征頻率相對應(yīng)，另外還可發(fā)現(xiàn)保持架故障特征頻率的不同倍頻，由此可以得出保持架處于動不平衡狀態(tài)[20]，與實(shí)際情況相一致。

圖5 IMF1的包絡(luò)譜Fig.5 Envelope spectrum of IMF1

3 結(jié)論

1)CEEMDAN算法有效降低了模態(tài)混疊，結(jié)合自相關(guān)系數(shù)、峭度指標(biāo)篩選法則，可準(zhǔn)確篩選出敏感IMF，剔除了無故障IMF的干擾。Hilbert包絡(luò)譜分析能從高頻采樣的頻譜中識別出振動能量較小的低頻部分，使之不會湮沒在復(fù)雜噪聲背景中，適合強(qiáng)噪聲條件下振動信號分析。

2) 由于結(jié)構(gòu)原因，軸承保持架轉(zhuǎn)速低于轉(zhuǎn)子；且剛度遠(yuǎn)小于內(nèi)圈、外圈、球等元件，故障振動頻率小、量級低。如果振動測試信號中出現(xiàn)保持架故障特征頻率，則反映保持架出現(xiàn)缺陷，處于動不平衡狀態(tài)，屬于軸承保持架動不平衡故障現(xiàn)象。