劉振峰，楊偉新

（中國(guó)航發(fā)湖南動(dòng)力機(jī)械研究所，湖南株洲 412002）

0 引言

航空發(fā)動(dòng)機(jī)故障診斷要求做好故障定位準(zhǔn)確快速，以方便后續(xù)的維修；另外，故障的診斷還可以對(duì)發(fā)動(dòng)機(jī)的實(shí)際運(yùn)行狀態(tài)進(jìn)行記錄與判斷，方便后續(xù)的繼續(xù)使用，確保發(fā)動(dòng)機(jī)始終處于安全狀態(tài)[1]。

EMD分解可將信號(hào)分解為多個(gè)IMF（Intrinsic Mode Function，本征模態(tài)函數(shù)）分量與一個(gè)趨勢(shì)的代數(shù)和，適用于非線性、非平穩(wěn)信號(hào)分析。基于EMD理論，以LabVIEW為開(kāi)發(fā)平臺(tái)，設(shè)計(jì)并開(kāi)發(fā)出一套功能完善的航空發(fā)動(dòng)機(jī)故障診斷系統(tǒng)。

1 彈支動(dòng)應(yīng)力測(cè)試技術(shù)

在彈支水平與垂直方向的彈條上粘貼應(yīng)變片。彈支自由端徑向位移與彈支彈條的應(yīng)變量的關(guān)系見(jiàn)式（1）[2]。式中Y為彈支自由端的徑向位移；l為彈條長(zhǎng)度；ε為彈支彈條的應(yīng)變量；h為彈條厚度。根據(jù)式（1），彈支自由端徑向位移與彈支彈條的應(yīng)變量成正比。根據(jù)彈條的因變量（動(dòng)應(yīng)力信號(hào)）可監(jiān)測(cè)發(fā)動(dòng)機(jī)轉(zhuǎn)子的故障。

發(fā)動(dòng)機(jī)轉(zhuǎn)子系統(tǒng)發(fā)生故障時(shí)，因受到剛度的非線性、油膜間隙及載荷變化等因素的影響，彈支的動(dòng)應(yīng)力信號(hào)必然呈現(xiàn)出非線性、非平穩(wěn)特性[3]，另外，由于測(cè)試過(guò)程中的噪聲污染，使得通過(guò)彈支動(dòng)應(yīng)力信號(hào)獲取的軸心軌跡形狀復(fù)雜，難以辨別出轉(zhuǎn)子系統(tǒng)的運(yùn)行狀態(tài)。因此，如何從含有背景噪聲且具有非線性、非平穩(wěn)特性的彈支動(dòng)應(yīng)力信號(hào)中提取出轉(zhuǎn)子系統(tǒng)故障特征信息，是航空發(fā)動(dòng)機(jī)轉(zhuǎn)子故障診斷系統(tǒng)的關(guān)鍵。

2 EMD基本原理[4]

EMD（Empirical Mode Decomposition，經(jīng)驗(yàn)?zāi)B(tài)分解）分解是將信號(hào)分解成多個(gè)IMF分量，可表示為式（2）。EMD分解由信號(hào)自身的變化決定，是一種自適應(yīng)信號(hào)分解方法，適用于分析航空發(fā)動(dòng)機(jī)轉(zhuǎn)子系統(tǒng)彈支非平穩(wěn)動(dòng)應(yīng)力信號(hào)。

3 總體設(shè)計(jì)

（1）硬件系統(tǒng)設(shè)計(jì)。硬件系統(tǒng)采用基于數(shù)據(jù)采集卡的虛擬儀器系統(tǒng)，采用日本共和電業(yè)應(yīng)變計(jì)，型號(hào)KFH-1-120-C1-11 H1M3，數(shù)據(jù)采集卡采用美國(guó)NI公司的PXI-4471B動(dòng)態(tài)信號(hào)采集卡。

（2）軟件構(gòu)成。軟件的設(shè)計(jì)是核心部分，LabVIEW語(yǔ)言擁有高級(jí)信號(hào)分析子函數(shù)以及各種儀器驅(qū)動(dòng)程序，用戶只需要將輸入控制、輸出顯示的圖標(biāo)調(diào)出來(lái)，就可以完成參數(shù)的配置，然后再與程序的框圖連接起來(lái)，就可以完成對(duì)應(yīng)的工作[5]。該系統(tǒng)主要由數(shù)據(jù)采集、數(shù)據(jù)分析、數(shù)據(jù)處理及故障診斷組成。

4 系統(tǒng)軟件設(shè)計(jì)

（1）數(shù)據(jù)采集模塊。支持工具中，不能缺少數(shù)據(jù)采集的驅(qū)動(dòng)程序及編程接口，其能夠提供硬件設(shè)備與應(yīng)用層開(kāi)發(fā)工具之間的接口，確保應(yīng)用層的開(kāi)發(fā)工具能被不同的硬件設(shè)備使用[6]，NI公司的NI-DAQmx驅(qū)動(dòng)程序，可以與LabVIEW和NI公司的硬件產(chǎn)品無(wú)縫結(jié)合使用，因此成為選用LabVIEW構(gòu)建數(shù)據(jù)采集系統(tǒng)時(shí)驅(qū)動(dòng)程序的首選。

（2）數(shù)據(jù)分析模塊。數(shù)據(jù)分析模塊可以對(duì)在線信號(hào)波形進(jìn)行動(dòng)態(tài)顯示，也能將數(shù)據(jù)作為文件存儲(chǔ)，方便后續(xù)使用。數(shù)據(jù)分析給出時(shí)域與頻域波形，界面主要有數(shù)據(jù)采集時(shí)的日期、時(shí)間及經(jīng)過(guò)的時(shí)間，另外給出一些參數(shù)設(shè)置，如分辨率、窗函數(shù)、坐標(biāo)軸顯示方式等，此外還顯示一些主要特征參數(shù)，如峰值、均值、有效值等。

（3）數(shù)據(jù)EMD處理。在Matlad的環(huán)境下實(shí)現(xiàn)EMD的使用，可將Matlab與LabVIEW結(jié)合進(jìn)行描述[2]。文本語(yǔ)言的可維護(hù)性與易讀性比較差，影響軟件后續(xù)的調(diào)試。針對(duì)EMD工具箱缺少EMD算法問(wèn)題，對(duì)LabVIEW進(jìn)行二次開(kāi)發(fā)[7]，EMD分解方法部分程序框圖如圖1所示。

5 應(yīng)用實(shí)例

某型航空發(fā)動(dòng)機(jī)轉(zhuǎn)子振動(dòng)信號(hào)通過(guò)安裝在發(fā)動(dòng)子轉(zhuǎn)子柔性支承－彈支彈條上的應(yīng)變片采集[3]，采樣頻率10 kHz。轉(zhuǎn)子轉(zhuǎn)速42 770 r/min。從時(shí)域圖（圖2）和頻譜圖（圖3）很難判斷轉(zhuǎn)子出現(xiàn)何種故障。

直接進(jìn)行故障信號(hào)EMD分解，基于其各個(gè)分量信號(hào)特性，通過(guò)4個(gè)分量的分解，采樣點(diǎn)數(shù)為1024.其中，信號(hào)最高頻成分

圖1 EMD分解程序框圖

圖2 時(shí)域波形圖

圖3 頻域波形圖

6 結(jié)語(yǔ)

基于LabVIEW的航空發(fā)動(dòng)機(jī)轉(zhuǎn)子故障診斷，目前主要是將計(jì)算機(jī)軟硬件技術(shù)、非平穩(wěn)信息處理技術(shù)以及應(yīng)變片測(cè)試技術(shù)進(jìn)行綜合運(yùn)用，然后在齒輪箱故障診斷中直接的使用虛擬儀器技術(shù)，可以滿足振動(dòng)信號(hào)的采集與處理要求，進(jìn)而將故障診斷的精度提高，縮短系統(tǒng)開(kāi)發(fā)周期，降低硬件維護(hù)費(fèi)用。通過(guò)實(shí)驗(yàn)證明，基于這種方式的故障診斷系統(tǒng)，工作穩(wěn)定可靠，能滿足在線與離線的故障分析，可以快速準(zhǔn)確地診斷發(fā)動(dòng)機(jī)轉(zhuǎn)子出現(xiàn)的故障。為INF1呈現(xiàn)出來(lái)的信息。EMD分解如圖4所示。

圖4 故障信號(hào)EMD分解圖

在主頻f=712.89 Hz（折合轉(zhuǎn)速為42 770 r/min）左右兩邊分別出現(xiàn)f1=628.66 Hz和f2=797.12 Hz的邊頻帶。又因f-f1=f2-f=84.23 Hz，根據(jù)文獻(xiàn)[8]，可以判斷轉(zhuǎn)子的故障為動(dòng)靜碰摩故障。該實(shí)例驗(yàn)證了轉(zhuǎn)子故障診斷系統(tǒng)的有效性。