基于特征頻譜對齒輪箱的振動監(jiān)測和故障診斷

肖湘平，陳立立，高 波

（荊州市巨鯨傳動機械有限公司，湖北荊州 434000）

0 引言

齒輪箱是現(xiàn)代工業(yè)生產(chǎn)制造中的一個重要部件，齒輪箱大量應(yīng)用在礦山、冶金、煤炭、電力、建材及軍工等行業(yè)中。但是由于齒輪箱結(jié)構(gòu)緊湊和運行工況的復(fù)雜性，并且齒輪箱通常是使用在高速重載的環(huán)境中，因此齒輪箱中的零部件極易因為疲勞失效而導(dǎo)致局部失效，使得設(shè)備無法正常運行。齒輪箱中齒輪故障占總故障的60%，因此對于齒輪的振動監(jiān)測和故障診斷很有必要[1-2]。

當(dāng)前，齒輪箱齒輪的振動監(jiān)測依賴于專業(yè)技術(shù)人員進行，這種方式的人工成本、時間成本高，并且隨著振動監(jiān)測技術(shù)的發(fā)展，監(jiān)測的齒輪箱越來越多，采集的數(shù)據(jù)量也越來越大，人工成本、時間成本也因此越來越高。在傳統(tǒng)的齒輪箱齒輪的診斷中，主要是通過收集齒輪箱關(guān)鍵部位的振動信號，然后對振動信號進行時域分析和頻域分析，實現(xiàn)對齒輪箱齒輪狀態(tài)的監(jiān)測，但是在監(jiān)測過程中對于信號的處理、識別、分析以及診斷結(jié)果的準(zhǔn)確性十分依靠專業(yè)人員的個人能力。因此提出了一種基于特征頻段的診斷深度神經(jīng)網(wǎng)絡(luò)，在傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程中往往需要大量的樣本去訓(xùn)練，但是通常在診斷早期是十分缺少樣本的，因此結(jié)合現(xiàn)有的齒輪診斷理論以及專業(yè)人員的經(jīng)驗，具體的是利用現(xiàn)有數(shù)據(jù)進行頻譜與包絡(luò)譜分析，然后提取相應(yīng)的故障頻段與無故障頻段。這樣，少量信號數(shù)據(jù)便可以提取足夠數(shù)量的頻段數(shù)據(jù)用以訓(xùn)練神經(jīng)網(wǎng)絡(luò)。在對故障診斷時，首先將信號提取為若干個頻段，然后利用訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)分別診斷每一個頻段的故障概率，最后綜合所有頻段的故障概率，輸出最終診斷結(jié)果。

1 基于振動信號的齒輪診斷

在通常情況下一對嚙合齒輪的重合度都不是一個整數(shù)，這就導(dǎo)致兩個齒輪在嚙合過程中參與嚙合的齒數(shù)會在一個值附近波動，因此嚙合點的結(jié)構(gòu)也會出現(xiàn)周期變化，導(dǎo)致齒輪嚙合點齒輪接觸力的周期變化，并且在嚙入和嚙出時會產(chǎn)生一個沖擊，對齒輪系統(tǒng)產(chǎn)生一個周期性激勵，從而形成了齒輪的嚙合點產(chǎn)生一個規(guī)律性振動。在齒輪系統(tǒng)無故障的情況下，嚙合點的振動信號可以表示為[3-4]：

式中 fZ——齒輪嚙合頻率，Hz

n——自然數(shù)，取1，2，3，…，M 為最大諧波次數(shù)

N——齒輪的旋轉(zhuǎn)速度，r/min

Z——齒輪齒數(shù)

定軸輪系齒輪故障的振動測試模型可以表示為：

其中fg為故障齒輪的故障頻率。根據(jù)恒等式：

其中，Jm（z）是自變量為z 的m 階第一類貝塞爾函數(shù)。以及歐拉公式，上式可以展開為：

對式（5）只考慮正頻率部分進行傅里葉變換為：

由式（6）可知，當(dāng)定軸輪系的齒輪發(fā)生故障時，振動信號以嚙合點的嚙合頻率fm為中心，在fm周圍分布以齒輪的故障頻率fg為間隔邊帶，因此頻譜峰值在fm±m(xù)fg處。在此處計算的均是只考慮了振動信號的基頻部分，在實際測試中還需要考慮齒輪嚙合頻率的倍頻，即考慮nfg，則振動信號可以表示為：

對式（7）只考慮正頻率部分，進行傅里葉變換為：

可以看出，在定軸輪系的嚙合頻率fm以及它的倍頻nfm周圍都分布以故障齒輪故障特征頻率fg為間隔的邊帶，頻譜峰值在nfm±m(xù)fg處。

2 基于頻段的樣本信號

2.1 故障信號樣本

根據(jù)前面基于振動信號的齒輪診斷，可以知道，在齒輪的故障診斷過程中主要的關(guān)心點為齒輪的嚙合頻率fm，以其邊帶與其之間的間隔頻率fg，通過對這2 個關(guān)鍵頻率值處識別和判斷是否存在明顯峰值，即可實現(xiàn)對齒輪的診斷。

一組振動信號的時域信號為：

對式（9）通過傅里葉變換為頻域信號：

將頻域信號按幅值從大到小排列，認(rèn)為前10%的幅值是明顯的峰值，并且將峰值頻率設(shè)定為“1”，非峰值頻率成分設(shè)定為“0”。

一對齒輪的參數(shù)見表1，輸入轉(zhuǎn)速為1488 r/min，根據(jù)故障頻率計算公式，計算可得到這對齒輪的故障頻率。

表1 齒輪參數(shù)

式中 d——滾動軸承滾動體直徑，mm

D——滾動軸承節(jié)徑，mm

α——滾動軸承壓力角，°

n——滾動體滾子數(shù)

高速級：foc1=159.1 Hz，fic1=237.7 Hz，fbc1=117.2 Hz，fcg1=159.1 Hz。

低速級：foc2=55.2 Hz，fic1=81.2 Hz，fbc1=41.9 Hz，fcg1=3.4 Hz。

只考慮嚙合頻率的1 倍頻及故障頻率3 倍頻的邊帶時，當(dāng)齒輪出現(xiàn)故障時，會在f={fm-3fg，fm-2fg，fm-fg，fm+fg，fm+2fg，fm+3fg}在對f 處的幅值均進行“0”和“1”處理，因此對于圖1 高速級齒輪故障頻譜的可以表示為A（f）={1，1，1，1，1，1，1，}，A（f）只能取“0”或“1”。對于齒輪故障A{f={fm-3fg，fm-2fg，fm-fg，fm，fm+fg，fm+2fg，fm+3fg}}的所有組合為：

上述的組合為高速級齒輪故障的一種頻率分布，但是在實際中齒輪出現(xiàn)故障時，不一定在fm周圍fm-3fg，fm-2fg，fm-fg，fm+fg，fm+2fg，fm+3fg都會出現(xiàn)。

取f={fm-3fg，fm-2fg，fm-fg，fm+fg，fm+2fg，fm+3fg}，即當(dāng)A（f′）≠{0，0，0，0，0，0}時都可以認(rèn)為齒輪出現(xiàn)了故障。高速級齒輪故障頻譜如圖1 所示。

圖1 高速級齒輪故障

在實際的頻譜中，由于傅里葉變換頻譜并不是連續(xù)的，而是一組離散數(shù)值，數(shù)值的間隔與頻譜的分辨率有關(guān)[5-6]，導(dǎo)致理論特征頻率可能會出現(xiàn)在2 個頻率值之間，并且轉(zhuǎn)速測量誤差也會使得計算得到的特征頻率值與實際頻率值存在誤差。這些影響都會導(dǎo)致實際的峰值并不是剛好在f 處的，因此需要對f 進行擴展。如嚙合頻率fm處，應(yīng)當(dāng)以最靠近fm離散頻率點frm 為中心，以頻譜的分辨率△Rf 為間隔，向正負(fù)量個方向各擴展3 個點，并對+3△Rf}七點的幅值進行“0”“1”處理，若A（）≠{0，0，0，0，0，0}，即存在5039 種組合則認(rèn)為A（fm）=1。對于f={fm-3fg，fm-2fg，fm，fm+fg，fm+2fg，fm+3fg}，均進行該擴展和“0”“1”處理。因此對于矩陣中所有的“1”都有5039 種組合，所以從A{f={fm-3fg，fm-2fg，fm-fg，fm，fm+fg，fm+2fg，fm+3fg}} 的所有組合中采用拉丁超立方抽樣抽取20 000 個樣本作為齒輪故障樣本。

2.2 無故障信號樣本

當(dāng)齒輪無故障時，在齒輪的嚙合頻率的周圍是不會出現(xiàn)邊帶的，因此當(dāng)出現(xiàn)A{f={fm-3fg，fm-2fg，fm-fg，fm，fm+fg，fm+2fg，fm+3fg}}所有的組合為：

因此也從無故障樣本中采用超立方拉丁抽樣，抽取20 000 個樣本用以去訓(xùn)練神經(jīng)網(wǎng)絡(luò)。所以將齒輪故障的20 000 個A {f} 樣本和齒輪無故障的20000 個A{f}作為輸入，齒輪出現(xiàn)故障時的輸出為“1”，無故障是的輸出為“0”，共40 000 個樣本訓(xùn)練神經(jīng)網(wǎng)絡(luò)。

3 基于特征頻帶的神經(jīng)網(wǎng)絡(luò)

在獲得充足的樣本后，需要通過樣本對神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練，本文采用的是深度神經(jīng)網(wǎng)絡(luò)，對40 000 個樣本，取60%作為訓(xùn)練數(shù)據(jù)，20%作為驗證數(shù)據(jù)，20%作為測試數(shù)據(jù)，訓(xùn)練結(jié)果如圖2所示。從圖中可以看出，基本上訓(xùn)練數(shù)據(jù)、驗證數(shù)據(jù)、測試數(shù)據(jù)均在零誤差線附近，說明該神經(jīng)網(wǎng)絡(luò)可用。

圖2 誤差直方圖

對于一組齒輪箱的振動信號模型的診斷流程如圖3 所示。首先對采集到的振動信號進行傅里葉變換，傅里葉頻譜的幅值進行“0”和“1”處理?？傻玫綄τ诿總€齒輪對應(yīng)的A{f={fm-3fg，fm-2fg，fm-fg，fm，fm+fg，fm+2fg，fm+3fg}}并對其中每個頻率進行擴展處理，通過訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)進行處理，根據(jù)輸出結(jié)果判斷各齒輪是否存在故障，如果輸出結(jié)果為“1”則齒輪已發(fā)生故障，反之輸出為“0”則齒輪無故障。

圖3 智能診斷流程

在計算各個軸承的故障頻率后，本系統(tǒng)采用包絡(luò)譜診斷神經(jīng)網(wǎng)絡(luò)來估計每一種可能的軸承故障發(fā)生的概率。

用此神經(jīng)網(wǎng)絡(luò)計算軸承概率時，分別將以其前3 倍故障概率為中心頻率的頻帶輸入神經(jīng)網(wǎng)絡(luò)求得前3 倍頻的故障概率，然后綜合前3 倍頻的故障概率得到該軸承的故障概率。

需要注意的是，由于計算到的軸承故障頻率是在理想情況下求得的。而齒輪箱實際工作情況很復(fù)雜，并且由于采樣原理，傅里葉變換得到的頻譜尖峰頻率往往與理論求得的故障頻率或故障頻率存在偏差。因此，在診斷中，會判斷故障頻率的窄邊帶內(nèi)所有采樣頻率點對應(yīng)的頻帶的故障概率。

4 結(jié)束語

通過結(jié)合齒輪和軸承故障診斷的理論方法，并將理論進行模糊化，確定了齒輪和軸承出現(xiàn)故障時的數(shù)值模型，并且考慮到在傅里葉變換過程中由于頻譜分辨率以及測量誤差，導(dǎo)致理論計算的故障特征頻率在傅里葉中的離散值中并不是完全對應(yīng)的，因此需要對這些特征頻率進行擴展，由此可得到充足的樣本，解決了在對神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練時需要大量樣本的問題，提高了在齒輪箱故障診斷中智能性，并且由于對樣本的模糊化，該診斷模型可對任意齒數(shù)、轉(zhuǎn)速的齒輪箱進行實時監(jiān)測，在齒輪或軸承出現(xiàn)故障時及時得出結(jié)果。

由于對于特征頻段僅僅是進行了一維的模糊處理，后面的科研人員，可對特征頻段進行二維的模糊化處理，使得該診斷系統(tǒng)能夠?qū)X輪的故障程度，甚至拓展到三維的模糊化，在時間軸上對特征進行模糊化，以便能預(yù)估故障的發(fā)展趨勢，以及齒輪的剩余壽命。