李國(guó)鴻

(中國(guó)飛行試驗(yàn)研究院，陜西西安710089)

航空發(fā)動(dòng)機(jī)振動(dòng)數(shù)據(jù)階比分析應(yīng)用研究

李國(guó)鴻

(中國(guó)飛行試驗(yàn)研究院，陜西西安710089)

試飛中發(fā)動(dòng)機(jī)振動(dòng)信號(hào)處理時(shí)，需獲得各轉(zhuǎn)子階次諧波分量的準(zhǔn)確幅值，以便對(duì)發(fā)動(dòng)機(jī)振動(dòng)進(jìn)行評(píng)估和監(jiān)視。針對(duì)該問題，對(duì)階比分析技術(shù)在該領(lǐng)域的具體應(yīng)用進(jìn)行了研究，并給出了實(shí)現(xiàn)方法。該方法采用sinc函數(shù)插值的方法，將等時(shí)間間隔采樣的離散數(shù)據(jù)序列轉(zhuǎn)換為等角度采樣的離散數(shù)據(jù)序列，隨后對(duì)等角度采樣后的數(shù)據(jù)序列進(jìn)行快速傅里葉變換(FFT)和諧波分解，來獲得信號(hào)階次譜。分析結(jié)果表明，給出的實(shí)現(xiàn)方法消除了FFT分析的截?cái)嘈?yīng)，提高了階次譜分量幅值的分析精度。

航空發(fā)動(dòng)機(jī)；截?cái)嘈?yīng)；階比分析；等角度采樣；階次譜

1 引言

目前，航空發(fā)動(dòng)機(jī)振動(dòng)信號(hào)處理，一般采用通用的基于FFT的頻譜分析方法。但這種分析方法在執(zhí)行過程中會(huì)引起頻譜分析誤差，且這些誤差可用混疊現(xiàn)象、欄柵效應(yīng)和截?cái)嘈?yīng)三種現(xiàn)象來描述[1～5]。雖然混疊現(xiàn)象可通過采樣前的抗混疊濾波消除，欄柵效應(yīng)可通過提高信號(hào)采樣率的方法保證，截?cái)嘈?yīng)可通過加窗處理來盡可能消除，但由于目前一般的振動(dòng)數(shù)據(jù)采集記錄設(shè)備采用的都是等時(shí)間間隔采樣方式，F(xiàn)FT分析數(shù)據(jù)塊時(shí)長(zhǎng)固定，分析頻率分辨率不變，在航空發(fā)動(dòng)機(jī)狀態(tài)變化時(shí)，其主要特征分量由于截?cái)嘈?yīng)帶來的誤差變化，不能同時(shí)兼顧各特征分量的截?cái)嗾`差。

階比分析技術(shù)以旋轉(zhuǎn)機(jī)械的轉(zhuǎn)速而不是時(shí)間、頻譜上的階比而不是頻率作為關(guān)注重點(diǎn)，因此一般階比譜圖上，信號(hào)的幅值或頻率是旋轉(zhuǎn)機(jī)械轉(zhuǎn)軸頻率階次的函數(shù)。在旋轉(zhuǎn)機(jī)械狀態(tài)變化時(shí)，某一階次的諧波成分就固定在一獨(dú)立的譜線上，不隨轉(zhuǎn)速變化，避免了一般通用頻譜分析中出現(xiàn)的諧波混疊效應(yīng)[6～9]。此外，由于一般階比分析技術(shù)都采用等角度采樣[8～12]，理論上可完全消除各階次諧波分量的截?cái)嘈?yīng)，保證各階次分量的截?cái)嗾`差。

航空發(fā)動(dòng)機(jī)主要特征分量的頻率為轉(zhuǎn)子轉(zhuǎn)頻的倍數(shù)，這與旋轉(zhuǎn)機(jī)械階比分析中階比的含義一致。同時(shí)考慮到階比分析的優(yōu)點(diǎn)，本文將結(jié)合目前試飛中航空發(fā)動(dòng)機(jī)振動(dòng)信號(hào)處理現(xiàn)狀，對(duì)如何將階比分析技術(shù)應(yīng)用于航空發(fā)動(dòng)機(jī)振動(dòng)數(shù)據(jù)處理領(lǐng)域進(jìn)行研究，以期進(jìn)一步提高振動(dòng)數(shù)據(jù)分析的能力和質(zhì)量。

2 階比分析理論

2.1 階比定義

旋轉(zhuǎn)機(jī)械振動(dòng)時(shí)，通常關(guān)注的振動(dòng)特征分量的振動(dòng)頻率為該旋轉(zhuǎn)機(jī)械轉(zhuǎn)速的倍數(shù)或分?jǐn)?shù)，該倍數(shù)和分?jǐn)?shù)就是階比。振動(dòng)分析中使用階比，主要是因?yàn)檗D(zhuǎn)速與階比的對(duì)應(yīng)關(guān)系恒定。階比與旋轉(zhuǎn)機(jī)械轉(zhuǎn)速的關(guān)系通?？捎孟率奖硎荆?/p>

式中：f為頻率，E為階比，n為轉(zhuǎn)速。

2.2 等角度采樣

與一般振動(dòng)信號(hào)采樣不同，采用階比分析時(shí)，被測(cè)振動(dòng)信號(hào)要經(jīng)過等角度增量的階比采樣變?yōu)榻怯螂x散序列，再經(jīng)離散傅立葉變換到階比域進(jìn)行分析。同步采樣后所得數(shù)據(jù)的采樣單位，是每轉(zhuǎn)的分?jǐn)?shù)而不是每秒的分?jǐn)?shù)。若再用這些數(shù)據(jù)進(jìn)行傅里葉變換，得到的測(cè)試結(jié)果將是一種頻率尺度類型的函數(shù)，以階比而不是赫茲為增量。

對(duì)于平穩(wěn)信號(hào)，若信號(hào)為整周期等角度采樣，則每次采樣截取的信號(hào)長(zhǎng)度都等于振動(dòng)信號(hào)周期的整數(shù)倍，且信號(hào)為等間隔采樣。頻譜分析時(shí)，分析譜線將落在振動(dòng)信號(hào)的頻率上，從而獲得比較準(zhǔn)確的頻譜，以避免頻譜泄漏和柵欄效應(yīng)對(duì)頻譜分析精度的影響。對(duì)于非平穩(wěn)信號(hào)，整周期等角度采樣可使信號(hào)在一定程度上變得平穩(wěn)，可減弱頻譜泄漏和柵欄效應(yīng)對(duì)分析精度的影響。

3 階比分析應(yīng)用

目前，發(fā)動(dòng)機(jī)試飛中所獲得的發(fā)動(dòng)機(jī)轉(zhuǎn)速數(shù)據(jù)和振動(dòng)數(shù)據(jù)，來自兩個(gè)不同的數(shù)據(jù)采集記錄系統(tǒng)。振動(dòng)數(shù)據(jù)以等時(shí)間間隔采樣方式獲得，對(duì)其的分析主要依賴事后的二次處理。要開展針對(duì)性的階比分析算法研究工作，必須考慮上述現(xiàn)狀。

圖1所示為充分考慮目前試飛中發(fā)動(dòng)機(jī)振動(dòng)數(shù)據(jù)分析現(xiàn)狀，并根據(jù)實(shí)際應(yīng)用需要，重點(diǎn)實(shí)現(xiàn)了整周期等角度采樣的發(fā)動(dòng)機(jī)振動(dòng)數(shù)據(jù)階比譜實(shí)現(xiàn)方法示意圖。圖中，n(n)為轉(zhuǎn)速數(shù)據(jù)序列，a(n)為振動(dòng)數(shù)據(jù)序列，ns(n)為時(shí)間同步后的轉(zhuǎn)速數(shù)據(jù)序列，as(n)為時(shí)間同步后的振動(dòng)數(shù)據(jù)序列，asb(n)為緩沖區(qū)振動(dòng)數(shù)據(jù)序列，asb,rs(l)為完成了等角度重采樣的振動(dòng)數(shù)據(jù)序列，asb,rs,w(l)為完成了等角度重采樣和數(shù)據(jù)截取后的振動(dòng)數(shù)據(jù)序列，Asb,rs(m)為序列asb,rs,w(l)的復(fù)頻譜，Asb,rs(m)為序列asb,rs,w(l)的諧波分解結(jié)果。

圖1 發(fā)動(dòng)機(jī)振動(dòng)數(shù)據(jù)階比譜實(shí)現(xiàn)方法示意圖Fig.1 Order spectrum implementing method of aero-engine vibration

圖中首先對(duì)n(n)和a(n)進(jìn)行時(shí)間同步處理的主要原因，是n(n)和a(n)來自兩個(gè)不同的數(shù)采系統(tǒng)，采樣率不同且差別很大；對(duì)振動(dòng)數(shù)據(jù)設(shè)置一個(gè)緩沖區(qū)，主要是為后續(xù)的FFT做準(zhǔn)備，緩沖區(qū)大小除必須考慮FFT分析的塊大小要求外，同時(shí)還要考慮重采樣和截?cái)嗟囊恍┨厥庖?。?shù)據(jù)截取過程可通過對(duì)序列asb,rs(l)加矩形窗來實(shí)現(xiàn)。

3.1 重采樣

選用sinc函數(shù)插值法實(shí)現(xiàn)等角度采樣。設(shè)a(n)、as(n)的采樣率為fs，需跟蹤的轉(zhuǎn)速信號(hào)頻率為fn，則等角度重采樣過程如式(2)所示：

asb(n)長(zhǎng)度N與序列asb,rs(l)長(zhǎng)度L之間的關(guān)系如式(3)所示：

設(shè)重采樣過程中對(duì)頻率為fn的正弦信號(hào)每周期采樣點(diǎn)數(shù)為2M1，序列asb,rs,w(l)的長(zhǎng)度為2M，則asb,rs(l)的采樣率為：

M、M1、M2三者之間的關(guān)系如式(5)所示：

3.2 緩沖區(qū)長(zhǎng)度

由于sinc(t) 函數(shù)是一個(gè)主瓣寬度為2、旁瓣寬度為1且迅速衰減的函數(shù)，而緩沖區(qū)振動(dòng)數(shù)據(jù)長(zhǎng)度固定，因此利用式(2)插值時(shí)，會(huì)在asb,rs(l)首、尾兩端產(chǎn)生振蕩效應(yīng)。鑒于此，為保證重采樣結(jié)果的準(zhǔn)確性，采取了使緩沖區(qū)振動(dòng)數(shù)據(jù)的時(shí)間長(zhǎng)度大于asb,rs,w(l)的時(shí)間長(zhǎng)度，然后在后續(xù)數(shù)據(jù)截取中所截取數(shù)據(jù)序列避開as,rs(l)首、尾兩端振蕩區(qū)的方法。實(shí)際應(yīng)用中，N的取值可參照下式進(jìn)行：

3.3 階次譜

對(duì)asb,rs,w(l)實(shí)施常規(guī)發(fā)動(dòng)機(jī)振動(dòng)信號(hào)處理中的FFT和諧波分解操作，就可得到需要的階次譜序列Asb,rs(m)。此時(shí)Asb,rs(m)的頻率分辨率為fn/2M2，階次分辨率為1/2M2，可用階次范圍為1/2M2～2M1-1。

4 仿真分析

4.1 sinc函數(shù)插值重采樣

圖2示出了某仿真信號(hào)重采樣測(cè)試結(jié)果。可見，正如前文分析結(jié)果，利用sinc函數(shù)插值完成重采樣，數(shù)據(jù)的首、尾兩端有明顯的振蕩效應(yīng)，且振幅迅速衰減(大約5個(gè)sinc(fs,at)函數(shù)的旁瓣寬度)，振蕩波形與sinc函數(shù)類似；其余部分波形與重采樣前的數(shù)據(jù)波形一致。

圖3示出了不同采樣率對(duì)仿真結(jié)果的影響?？梢?，重采樣前信號(hào)采樣率fs,a，對(duì)重采樣后數(shù)據(jù)首、尾兩端振蕩效應(yīng)的影響不同，fs,a越大影響越明顯。

4.2 穩(wěn)態(tài)階次譜

仿真高低壓轉(zhuǎn)子頻率分別為：fn,h=181.25 Hz，fn,l=60.00 Hz。仿真振動(dòng)信號(hào)如式(7)所示：

式中：A1=1，f1=fn,l，φ1=0，A2=2，f2=fn,h，φ2= 0.3π。fs=5000，振動(dòng)仿真信號(hào)時(shí)長(zhǎng)T=5 s。

仿真獲得的低壓、高壓階次譜如圖4所示，常規(guī)FFT頻譜如圖5所示。對(duì)比兩圖可發(fā)現(xiàn)，由階次譜獲得的低壓轉(zhuǎn)子轉(zhuǎn)速nl和高壓轉(zhuǎn)子轉(zhuǎn)速nh基頻分量的幅值誤差，總體上明顯小于由常規(guī)FFT頻譜獲得的基頻分量的誤差。此外，常規(guī)FFT頻譜分析中FFT數(shù)據(jù)塊大小對(duì)轉(zhuǎn)子基頻分量幅值誤差的影響較大，且具有一定隨機(jī)性；但M對(duì)由階次譜獲得的nl和nh基頻分量幅值的影響有規(guī)律，即M越大誤差越小。

圖2 某仿真信號(hào)重采樣測(cè)試結(jié)果Fig.2 Test results of a simulating signal

圖3 不同采樣率對(duì)仿真結(jié)果的影響Fig.3 Effects of different sample rate on simulation results

5 實(shí)測(cè)數(shù)據(jù)分析

利用某渦扇發(fā)動(dòng)機(jī)試飛中加裝的加速度振動(dòng)傳感器記錄的實(shí)測(cè)數(shù)據(jù)，對(duì)前文所述階比分析算法進(jìn)行進(jìn)一步分析和驗(yàn)證。

5.1 穩(wěn)態(tài)階次譜

圖6所示為某渦扇發(fā)動(dòng)機(jī)某穩(wěn)態(tài)高、低壓階次譜分析結(jié)果?？梢?，利用階次譜可更清晰地分辨出發(fā)動(dòng)機(jī)轉(zhuǎn)子基頻和倍頻分量(低壓轉(zhuǎn)子階次分別為2、3、6、11、17的振動(dòng)分量，高壓轉(zhuǎn)子階次分別為1、2、11的振動(dòng)分量)。與常規(guī)的FFT頻譜相比，由于階次譜的橫坐標(biāo)為階比，因此利用階次譜識(shí)別高、低壓轉(zhuǎn)子基頻和倍頻分量時(shí)更直觀。此外，由前文的分析可知，由階次譜獲得轉(zhuǎn)子基頻、倍頻分量幅值的誤差更小。

圖4 穩(wěn)態(tài)仿真數(shù)據(jù)階次譜Fig.4 Order spectrum simulating results in stable state

圖5 穩(wěn)態(tài)仿真數(shù)據(jù)常規(guī)FFT頻譜Fig.5 FFT spectrum simulating results in stable state

5.2 過渡穩(wěn)態(tài)階次譜

圖7所示為利用階比分析方法，對(duì)某渦扇發(fā)動(dòng)機(jī)某過渡態(tài)振動(dòng)數(shù)據(jù)進(jìn)行分析的高、低壓轉(zhuǎn)子階次譜云圖?？梢?，過渡態(tài)階次譜可以很清楚地反映轉(zhuǎn)子基頻、倍頻分量幅值隨時(shí)間的變化。圖7(a)中，對(duì)應(yīng)時(shí)間歷程曲線nl=41%、nh=74%的時(shí)刻(此時(shí)刻發(fā)動(dòng)機(jī)轉(zhuǎn)速比約為3)，nl的3倍頻或nh的1倍頻、nl的6倍頻或nh的2倍頻、nl的14倍頻分量比較明顯(在云圖上顯示為一條跟背景噪聲明顯區(qū)分的亮線)；其余時(shí)刻，nl的5倍頻分量比較明顯。與穩(wěn)態(tài)階次譜類似，利用過渡態(tài)階次譜可更直觀地識(shí)別過渡態(tài)高、低壓轉(zhuǎn)子基頻和倍頻分量的變化，且精度更高。

5.3 分量跟蹤

圖8所示為在5.2節(jié)基礎(chǔ)上，對(duì)高、低壓轉(zhuǎn)子過渡態(tài)基頻分量進(jìn)行跟蹤的結(jié)果?？梢?，通過階次譜獲得的高、低壓轉(zhuǎn)子基頻分量在整個(gè)時(shí)間歷程中的變化趨勢(shì)，分別與高、低壓轉(zhuǎn)子轉(zhuǎn)速的變化趨勢(shì)類似。在相對(duì)時(shí)間35～50 s期間，有一個(gè)比較明顯的增大，這與轉(zhuǎn)子基頻分量主要反映轉(zhuǎn)子系統(tǒng)的不平衡量、一般情況下其幅值隨發(fā)動(dòng)機(jī)轉(zhuǎn)速增大而增大的原理一致。

圖6 某渦扇發(fā)動(dòng)機(jī)穩(wěn)態(tài)階次譜Fig.6 Order spectrum of an aero-engine in stable state

圖7 某渦扇發(fā)動(dòng)機(jī)過渡態(tài)階次譜Fig.7 Rotor order spectrum of an aero-engine in transition state

圖8 某渦扇發(fā)動(dòng)機(jī)過渡態(tài)轉(zhuǎn)速與基頻分量幅值的時(shí)間歷程Fig.8 The relation between rotating speed and basic order amplitude of an aero-engine in transition state

與穩(wěn)態(tài)階次譜和過渡態(tài)階次譜分析結(jié)果類似，利用階比分析技術(shù)跟蹤轉(zhuǎn)子分量，所獲得的分量幅值由于徹底消除了截?cái)嘈?yīng)帶來的誤差，其分析誤差更小，精度更高。

6 結(jié)論

(1)階比分析方法由于采用了整周期等角度采樣技術(shù)，理論上可徹底消除常規(guī)FFT頻譜分析中的截?cái)嘈?yīng)，提高所跟蹤轉(zhuǎn)子基頻、倍頻分量的精度。

(2)利用sinc函數(shù)插值的方法，可實(shí)現(xiàn)重采樣前、后數(shù)據(jù)序列采樣率之比為任意值的整周期等角度采樣，但處理時(shí)必須注意重采樣后數(shù)據(jù)序列首、尾兩端的振蕩效應(yīng)。

(3)階比分析中，M的取值對(duì)處理結(jié)果有一定影響，M取值越大，階比譜分量的幅值誤差就越小，但效率會(huì)下降。

(4)與常規(guī)FFT頻譜相比，利用階比譜診斷發(fā)動(dòng)機(jī)振動(dòng)故障時(shí)，高、低壓轉(zhuǎn)子的基頻和倍頻分量更容易識(shí)別，這主要是因?yàn)殡A比與發(fā)動(dòng)機(jī)高、低壓轉(zhuǎn)子轉(zhuǎn)速之間存在著恒定的對(duì)應(yīng)關(guān)系。

(5)利用階比譜分析結(jié)果，可跟蹤過渡態(tài)發(fā)動(dòng)機(jī)轉(zhuǎn)子基頻、倍頻振動(dòng)分量，且所獲得的振動(dòng)分量的幅值誤差，總體上要小于基于常規(guī)FFT的分量跟蹤結(jié)果。

[1]胡廣書.數(shù)字信號(hào)處理理論、算法與實(shí)現(xiàn)[M].2版.北京：清華大學(xué)出版社，2003.

[2]丁玉美，高西全.數(shù)字信號(hào)處理[M].2版.西安：西安電子科技大學(xué)出版社，2000.

[3]屈玉池，魏海濤，趙述元.動(dòng)態(tài)信號(hào)處理中的問題及改進(jìn)[J].飛行試驗(yàn)，2001，17(2)：5—10.

[4]李杭生，陳丹.頻譜分析中窗函數(shù)的研究[J].微計(jì)算機(jī)信息，2008，24(1-4)：272—273.

[5]劉廣臣，張惠安，賈愛賓.數(shù)字信號(hào)處理中的加窗問題研究[J].長(zhǎng)沙大學(xué)學(xué)報(bào)，2003，17(4)：59—62.

[6]宋文健.階比分析方法的研究[D].河北：燕山大學(xué)，2006.

[7]楊志堅(jiān)，丁康，梁茜.基于頻譜校正理論的階比跟蹤分析[J].機(jī)械工程學(xué)報(bào)，2009，45(12)：41—45.

[8]楊通強(qiáng)，鄭海起，龔烈航，等.計(jì)算階次分析中的采樣率設(shè)置準(zhǔn)則[J].中國(guó)工程機(jī)械學(xué)報(bào)，2011，9(1)：14—18.

[9]高艷，鄭華文，陳永星.階比跟蹤中的惠特克插值[J].昆明理工大學(xué)學(xué)報(bào)(理工版)，2007，32(3)：30—34.

[10]汪偉，楊通強(qiáng)，王紅，等.非穩(wěn)態(tài)信號(hào)計(jì)算階次分析中的重采樣率研究[J].振動(dòng)、測(cè)試與診斷，2009，29(3)：349—351.

[11]溫政鋼，欒軍英.基于定點(diǎn)DSP的計(jì)算階次跟蹤研究[J].裝甲兵工程學(xué)院學(xué)報(bào)，2010，24(4)：62—64.

[12]郭瑜，秦樹人，湯寶平.基于分段重疊零相位濾波的階比跟蹤濾波法[J].振動(dòng)工程學(xué)報(bào)，2003，16(4)：399—403.

Order Tracing Application Study on Aero-Engine Vibration Data Processing

LI Guo-hong
(Chinese Flight Test Establishment，Xi’an 710089，China)

While vibration data is processed in aero-engine flight tests,the accurate amplitude of the rotors is required for vibration analysis and monitor.In this paper,how to use order tracing technology to solve aero-engine vibration data processing problem in flight tests was studied carefully;moreover the realizable or?der tracing numerical algorithm was presented too.Sinc function interpolation method was used to get the sampling angles discrete data sequence from sampling time span.Then FFT and harmonic decomposition was complied on the new data sequence to obtain the signal order spectrum.The analysis results show that the presented method can avoid the truncation effect,and improve the accuracy of the order amplitude analysis.

aero-engine；truncation effect；order tracing；uniform angle sampling；order spectrum

V231.92

：A

：1672-2620(2014)04-0042-06

2013-09-19；

：2014-01-14

李國(guó)鴻(1977-)，男，陜西鳳翔人，高級(jí)工程師，碩士，主要從事航空動(dòng)力裝置飛行試驗(yàn)相關(guān)的結(jié)構(gòu)強(qiáng)度試飛，以及相關(guān)的發(fā)動(dòng)機(jī)振動(dòng)測(cè)試、振動(dòng)數(shù)據(jù)分析等研究工作。