郭振偉，于瀾濤，陳卓，黃首清，路彤

（1.中國空間技術(shù)研究院，北京 100094； 2.航天機(jī)電產(chǎn)品環(huán)境可靠性試驗(yàn)技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室，北京 100094；3.北京衛(wèi)星環(huán)境工程研究所，北京 100094）

引言

滾動(dòng)軸承是旋轉(zhuǎn)機(jī)械的關(guān)鍵零部件，軸承的可靠性高低往往會(huì)影響整個(gè)機(jī)械設(shè)備的運(yùn)行狀態(tài)[1]。導(dǎo)致軸承失效的因素[2]有很多，例如安裝不當(dāng)、潤(rùn)滑失效、污染、磕碰、惡劣的使用環(huán)境等，因此軸承實(shí)際壽命與設(shè)計(jì)壽命之間常常有較大偏差，且具有分散性。因此，監(jiān)測(cè)軸承的工作狀態(tài)并據(jù)此進(jìn)行故障診斷，可以及時(shí)發(fā)現(xiàn)軸承失效的征兆，預(yù)防和減少軸承失效造成的損失。

目前關(guān)于軸承故障診斷的文獻(xiàn)較多[3]，根據(jù)故障信號(hào)類型[4]進(jìn)行分類包括溫度檢測(cè)法、油液檢測(cè)法、振動(dòng)分析法、聲發(fā)射法、電流檢測(cè)法，其中振動(dòng)分析法是目前使用最多、最有效的方法之一。根據(jù)故障特征提取方式分類，包括時(shí)域分析和頻域分析兩大類，常見的時(shí)域特征參數(shù)包括均值、均方根、峰值、峭度等，常見的頻域特征參數(shù)包括平均能量、頻譜集中程度、主頻帶位置等。根據(jù)故障診斷模型進(jìn)行分類，包括數(shù)據(jù)驅(qū)動(dòng)方法和失效物理方法，目前數(shù)據(jù)驅(qū)動(dòng)方法[5-7]發(fā)展迅猛，主要包括神經(jīng)網(wǎng)絡(luò)方法、隨機(jī)森林法、支持向量機(jī)等，特點(diǎn)是準(zhǔn)確性較高，但常常需要?dú)v史數(shù)據(jù)和經(jīng)驗(yàn)支撐；失效物理方法包括特征頻率法[8]、譜峭度法[9]，具有不依賴歷史數(shù)據(jù)的特點(diǎn)，更適用于工程實(shí)際中的小樣本故障診斷場(chǎng)景。綜上，大部分文獻(xiàn)側(cè)重于對(duì)故障信號(hào)處理、故障診斷模型的研究，所研究的軸承常常只是早期輕微故障且具有多故障耦合的特點(diǎn)，針對(duì)典型單一故障模式開展故障注入軸承試驗(yàn)件的研究并不多見。

航天領(lǐng)域針對(duì)軸承產(chǎn)品的可靠性有極高要求，意外故障可能會(huì)導(dǎo)致重大損失，例如美軍HESSI衛(wèi)星在振動(dòng)試驗(yàn)中就因?yàn)樵囼?yàn)臺(tái)軸承故障造成整星損毀[10]。有很多文獻(xiàn)重點(diǎn)關(guān)注軸承加速壽命試驗(yàn)的試驗(yàn)方法研究。徐東等[11]深入分析滾動(dòng)軸承加速壽命試驗(yàn)過程中遇到的各種問題, 提出一套完整的滾動(dòng)軸承加速壽命試驗(yàn)方法。李新立等[12]針對(duì)某航天器天線指向機(jī)構(gòu)固體潤(rùn)滑球軸承，開展了恒定應(yīng)力加速壽命試驗(yàn)方法研究。黃小凱等[13]建立了基于響應(yīng)面方法辨識(shí)不同應(yīng)力與軸承失效響應(yīng)值之間的加速模型，并對(duì)其進(jìn)行了試驗(yàn)驗(yàn)證。基于加速壽命試驗(yàn)的數(shù)據(jù)，還可利用各種模型對(duì)滾動(dòng)軸承的壽命與可靠性開展研究[14-16]。

本文針對(duì)多個(gè)故障注入軸承試驗(yàn)件，通過搭建的軸承試驗(yàn)平臺(tái)，對(duì)各故障注入軸承試驗(yàn)件進(jìn)行對(duì)比試驗(yàn)，利用故障特征頻率分析方法對(duì)比分析了不同轉(zhuǎn)速下各故障軸承的頻率特征圖像，實(shí)現(xiàn)對(duì)滾動(dòng)軸承故障位置和故障類型的診斷，積累了故障位置為內(nèi)圈、外圈、滾珠以及故障模式為擦傷、裂紋的軸承對(duì)應(yīng)的振動(dòng)信號(hào)特征數(shù)據(jù)。

1 試驗(yàn)系統(tǒng)

參考圖1，本文采用的試驗(yàn)系統(tǒng)包括硬件和軟件兩部分，硬件部分主要包括兩個(gè)軸系，每個(gè)軸系有兩個(gè)軸承腔，故障軸承安裝在軸系一端的軸承腔（1-2和2-2）中，便于拆裝。兩個(gè)軸系由PC端控制轉(zhuǎn)速的電機(jī)驅(qū)動(dòng)，電機(jī)驅(qū)動(dòng)端與軸承腔間安裝了摩擦力矩傳感器，軸承腔外表面安裝了熱電偶和振動(dòng)傳感器，可實(shí)時(shí)監(jiān)測(cè)兩個(gè)軸系的摩擦力矩和各軸承腔的溫度和振動(dòng)信號(hào)，其中振動(dòng)信號(hào)是本案例研究中重點(diǎn)分析的對(duì)象。硬件部分還包括變頻器、溫度變送器、信號(hào)調(diào)理器、數(shù)據(jù)采集板卡等，用于試驗(yàn)數(shù)據(jù)測(cè)控。

圖1 軸承試驗(yàn)系統(tǒng)的軸系和傳感器

試驗(yàn)系統(tǒng)的監(jiān)控軟件界面如圖2所示，由Labview程序編寫，可操作電機(jī)按照設(shè)定的轉(zhuǎn)速運(yùn)轉(zhuǎn)并實(shí)現(xiàn)無級(jí)變速，記錄兩個(gè)軸系振動(dòng)、摩擦力矩以及溫度信號(hào)。

圖2 軸承試驗(yàn)系統(tǒng)的監(jiān)控軟件界面

2 故障注入軸承試驗(yàn)件

軸承試驗(yàn)件為HRB7004 CETA/P5型角接觸球軸承，軸承節(jié)圓直徑取D = 31 mm，滾珠個(gè)數(shù)Z = 13，滾珠直徑d = 6.35 mm，接觸角α = 15 °。按照故障類型和位置的不同，將各故障軸承試驗(yàn)件如表1進(jìn)行編號(hào)，故障位置照片見圖3。

圖3 各故障軸承試驗(yàn)件的故障照片

表1 故障注入軸承試驗(yàn)件列表

試驗(yàn)時(shí)，將故障軸承和用于對(duì)照的正常軸承同時(shí)安裝在同一個(gè)軸系上，分別在1 000 rpm、1 500 rpm、3 000 rpm三種轉(zhuǎn)速下對(duì)軸承的振動(dòng)信號(hào)進(jìn)行采集，并對(duì)采集到的振動(dòng)信號(hào)進(jìn)行FFT（Fast Fourier Transform，快速傅里葉變換）處理得到對(duì)應(yīng)的頻率特征圖。1 000 rpm時(shí)無故障軸承L00的頻率特征圖如圖4所示。頻率特征圖用于下一步的頻率特征分析和故障定位。

圖4 1 000 rpm時(shí)L00軸承的頻率特征圖

從圖中可以看出無故障時(shí)，軸承的振動(dòng)信號(hào)頻域分布相對(duì)均勻，幅值隨頻率有一定的起伏波動(dòng)，但是并沒有表現(xiàn)出在某個(gè)特定頻率下的異常極大值，最高點(diǎn)仍小于400。后續(xù)的案例診斷主要針對(duì)特征頻率附近軸承振動(dòng)信號(hào)的幅值變化。

3 典型故障位置的故障特征頻率計(jì)算方法

軸承部件如外圈、內(nèi)圈或滾珠，當(dāng)其接觸的工作面產(chǎn)生了擦傷、裂紋等局部的損傷故障時(shí)，會(huì)以特定的通過頻率引發(fā)一連串的沖擊振動(dòng)，這一通過頻率就為故障特征頻率。不同故障位置的故障特征頻率如下：

內(nèi)圈故障的特征頻率：

滾珠故障的特征頻率：

式中：

d—滾珠直徑；

D—軸承節(jié)圓直徑；

α—接觸角；

Z—滾珠個(gè)數(shù)；

n—轉(zhuǎn)速。

代入本次試驗(yàn)中的軸承數(shù)據(jù)，得到各個(gè)轉(zhuǎn)速下軸承的故障特征頻率，如表2所示。

表2 各個(gè)轉(zhuǎn)速下的故障特征頻率

4 故障診斷結(jié)果和討論

基于L06、L08和L10共3個(gè)故障軸承和1個(gè)無故障軸承L00的頻率特征圖，分別繪制內(nèi)圈故障特征頻率if，外圈故障特征頻率of和滾珠故障特征頻率rf處的幅值對(duì)比柱狀圖分別如表4~6所示。

無故障軸承L00的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)如表3所示，作為其他故障軸承數(shù)據(jù)的對(duì)比基準(zhǔn)。

表3 無故障軸承L00在各故障特征頻率處幅值

故障軸承L06的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)如表4所示。結(jié)合非故障軸承和其它故障軸承的振動(dòng)信號(hào)數(shù)據(jù)可以看出，故障軸承L06內(nèi)圈故障特征頻率（fi）幅值在3 000 rpm下明顯偏大，在1 000 rpm、1 500 rpm下相較于外圈故障特征頻率（fo）和滾動(dòng)體故障特征頻率（fr）也偏大，故可以推斷，L06的主要故障類型為外圈故障。

表4 故障軸承L06的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)

故障軸承L08的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)如表5所示?？梢钥闯?，故障軸承L08在1 000 rpm，1 500 rpm和3 000 rpm三個(gè)轉(zhuǎn)速下各個(gè)特征頻率對(duì)應(yīng)的振幅都比無故障軸承高，但是滾動(dòng)體故障特征頻率（fr）幅值明顯偏高。推斷其故障類型為滾動(dòng)體故障。

表5 故障軸承L08的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)

故障軸承L10的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)如表6所示?？梢钥闯觯收陷S承L10在三個(gè)轉(zhuǎn)速下外圈故障特征頻率（fo）幅值異常偏高，滾動(dòng)體故障特征頻率（fr）幅值在1 000 rpm轉(zhuǎn)速下偏大，在1 500 rpm及3 000 rpm下逐漸回歸正常范圍內(nèi)。故推斷L10主要故障類型為外圈損壞。

表6 故障軸承L10的振動(dòng)信號(hào)在各故障特征頻率處的幅值數(shù)據(jù)

針對(duì)三個(gè)軸承故障診斷結(jié)果進(jìn)行匯總，如表7所示。綜上所述，采用本文試驗(yàn)平臺(tái)和故障特征頻率分析方法能較為準(zhǔn)確地對(duì)滾動(dòng)軸承的故障進(jìn)行診斷。

表7 故障注入軸承案例診斷結(jié)果

5 結(jié)論

本文針對(duì)多個(gè)故障注入滾動(dòng)軸承試驗(yàn)件，設(shè)計(jì)了軸承試驗(yàn)平臺(tái)，利用故障特征頻率分析方法對(duì)故障位置和故障類型的診斷，主要結(jié)論如下：

1）所設(shè)計(jì)的軸承試驗(yàn)系統(tǒng)可實(shí)時(shí)監(jiān)測(cè)兩個(gè)軸系的摩擦力矩和各軸承腔的溫度和振動(dòng)信號(hào)。

2）針對(duì)HRB7004 CETA/P5型角接觸球軸承設(shè)計(jì)了L06（擦傷故障，內(nèi)圈溝道）、L08（擦傷故障，滾珠）、L10（裂紋故障，外圈溝道）3個(gè)故障注入試驗(yàn)件，測(cè)試了1 000 rpm、1 500 rpm和3 000 rpm三種轉(zhuǎn)速下的振動(dòng)數(shù)據(jù)。

3）計(jì)算HRB7004 CETA/P5型角接觸球軸承在三種轉(zhuǎn)速下的故障特征頻率，以1 000 rpm為例，內(nèi)圈、外圈、滾動(dòng)體故障特征頻率分別是129.77 Hz、86.9 Hz和39.09 Hz。

4）基于故障特征頻率頻率法推斷出故障軸承L06、L08和L10的故障位置和嚴(yán)重程度，與實(shí)際情況吻合。

后續(xù)將針對(duì)更多故障位置和故障類型進(jìn)行研究，細(xì)化對(duì)故障類型的診斷精度。