孫百紅,田 川

(1.西安航天動(dòng)力研究所，陜西 西安 710100;2.空軍裝備部駐西安地區(qū)軍事代表局，陜西 西安 710100)

0 引言

液體火箭發(fā)動(dòng)機(jī)研制過(guò)程中要進(jìn)行大量的地面熱試車(chē)考核，發(fā)動(dòng)機(jī)及其組件在試車(chē)過(guò)程中要承受其系統(tǒng)本身產(chǎn)生的振動(dòng)、沖擊、噪聲及壓力脈動(dòng)等復(fù)雜動(dòng)態(tài)載荷的作用，過(guò)大或是異常的動(dòng)力學(xué)環(huán)境有時(shí)會(huì)引起發(fā)動(dòng)機(jī)結(jié)構(gòu)和設(shè)備的失效或者破壞，嚴(yán)重時(shí)發(fā)動(dòng)機(jī)被燒毀、試車(chē)臺(tái)被嚴(yán)重破壞。準(zhǔn)確有效地測(cè)量、分析發(fā)動(dòng)機(jī)試車(chē)過(guò)程中的動(dòng)態(tài)響應(yīng)對(duì)發(fā)動(dòng)機(jī)研制至關(guān)重要,也是獲取發(fā)動(dòng)機(jī)動(dòng)力學(xué)環(huán)境最直接、最有效的手段，是發(fā)動(dòng)機(jī)結(jié)構(gòu)設(shè)計(jì)、改進(jìn)和定型最具有說(shuō)服力的依據(jù)。對(duì)出現(xiàn)故障的試車(chē)，不同故障以不同的特征形式表現(xiàn)出來(lái)，認(rèn)真處理分析包含著故障信息的試車(chē)信號(hào)，有助于確定故障類(lèi)型、性質(zhì)及其產(chǎn)生的原因，有助于掌握各種故障的特征，以便采取相應(yīng)的措施，減少和預(yù)防故障重現(xiàn)。對(duì)成功試車(chē)，處理分析這些動(dòng)態(tài)響應(yīng)信號(hào)有助于了解發(fā)動(dòng)機(jī)及其組件不同工況下的動(dòng)力學(xué)特征，從而為故障分析提供對(duì)比依據(jù),對(duì)發(fā)動(dòng)機(jī)可靠性增長(zhǎng)試車(chē)前風(fēng)險(xiǎn)分析、試車(chē)過(guò)程中狀態(tài)監(jiān)測(cè)、試車(chē)后性能評(píng)估、試車(chē)故障診斷分析及發(fā)動(dòng)機(jī)組件力學(xué)環(huán)境試驗(yàn)設(shè)計(jì)等方面都具有極其重要的作用[1]。

1 試車(chē)動(dòng)態(tài)信號(hào)處理與分析方法

1.1 試車(chē)信號(hào)的分類(lèi)及特征

液體火箭發(fā)動(dòng)機(jī)試車(chē)速變信號(hào)可以分為以下4類(lèi)：平穩(wěn)隨機(jī)信號(hào)、非平穩(wěn)隨機(jī)信號(hào)、瞬態(tài)信號(hào)及周期信號(hào)與隨機(jī)信號(hào)的疊加。發(fā)動(dòng)機(jī)工作過(guò)程中，傳感器測(cè)量到的信號(hào)主要以后兩類(lèi)最為常見(jiàn)。發(fā)動(dòng)機(jī)的起動(dòng)、關(guān)機(jī)過(guò)程就是比較明顯的時(shí)變、瞬態(tài)信號(hào)特征，發(fā)動(dòng)機(jī)轉(zhuǎn)級(jí)過(guò)程表現(xiàn)為非平穩(wěn)隨機(jī)信號(hào)，發(fā)動(dòng)機(jī)穩(wěn)定工作段通常近似為一個(gè)平穩(wěn)的各態(tài)歷經(jīng)過(guò)程，測(cè)量信號(hào)則以周期信號(hào)與隨機(jī)信號(hào)的疊加為主要特征。

液體火箭發(fā)動(dòng)機(jī)的心臟——泵系統(tǒng)，其試車(chē)振動(dòng)信號(hào)主要表現(xiàn)為旋轉(zhuǎn)機(jī)械轉(zhuǎn)速頻率及由其結(jié)構(gòu)葉片數(shù)相關(guān)的倍頻以及流體噪聲等，發(fā)動(dòng)機(jī)的燃燒組件如發(fā)生器、推力室等其試車(chē)振動(dòng)信號(hào)特征主要表現(xiàn)為燃燒主頻、燃燒振蕩、燃燒噪聲及由旋轉(zhuǎn)機(jī)械傳遞而來(lái)的轉(zhuǎn)速及其倍頻等。

1.2 試車(chē)信號(hào)的分析方法

液體火箭發(fā)動(dòng)機(jī)試車(chē)過(guò)程動(dòng)態(tài)信號(hào)一般需要進(jìn)行時(shí)域、頻域以及時(shí)頻域分析。時(shí)域主要關(guān)注發(fā)動(dòng)機(jī)起動(dòng)與關(guān)機(jī)動(dòng)作產(chǎn)生的沖擊時(shí)刻及幅值以及沖擊過(guò)后信號(hào)的變化趨勢(shì)走向、全程時(shí)域幅值及RMS(Root-Mean-Square,均方根)曲線的變化情況，頻域主要關(guān)注發(fā)動(dòng)機(jī)穩(wěn)態(tài)工作段的頻率成分及其幅值、指定頻段內(nèi)的RMS值大??；時(shí)頻分析主要關(guān)注發(fā)動(dòng)機(jī)工作過(guò)程頻率成分及幅值隨時(shí)間變化情況。

根據(jù)液體火箭發(fā)動(dòng)機(jī)熱試車(chē)不同工作段的信號(hào)特點(diǎn)及關(guān)注點(diǎn)，試車(chē)速變信號(hào)分析的內(nèi)容、步驟及方法如圖1所示[2-4]。

圖1 發(fā)動(dòng)機(jī)試車(chē)動(dòng)態(tài)信號(hào)分析方法Fig.1 The analysis method of dynamic data from rocket engine tests

2 某發(fā)動(dòng)機(jī)試車(chē)信號(hào)分析

某型號(hào)發(fā)動(dòng)機(jī)試車(chē)多次出現(xiàn)結(jié)構(gòu)焊縫開(kāi)裂、燃料泄漏并起火等故障，圖2是試后檢查發(fā)現(xiàn)焊縫開(kāi)裂情況，圖3是試車(chē)過(guò)程中的燃料泄漏情況。采用發(fā)動(dòng)機(jī)信號(hào)分析方法，對(duì)該發(fā)動(dòng)機(jī)16次試車(chē)某組件的振動(dòng)信號(hào)進(jìn)行對(duì)比分析，包括3次出現(xiàn)故障的試車(chē)和13次沒(méi)有出現(xiàn)故障的試車(chē)信號(hào)。已知該發(fā)動(dòng)機(jī)每次試車(chē)先試額定工況后試過(guò)載工況，起動(dòng)和關(guān)機(jī)信號(hào)分析沒(méi)有發(fā)現(xiàn)異常，以下是穩(wěn)態(tài)段信號(hào)的分析情況。

2.1 頻譜分析

某振動(dòng)參數(shù)穩(wěn)態(tài)段頻譜分析如圖4所示，圖4(a)是出現(xiàn)故障的試車(chē)的頻譜，圖4(b)是大部分沒(méi)有出現(xiàn)故障的試車(chē)的頻譜，橫軸為頻率，縱軸為幅值?？梢?jiàn)當(dāng)發(fā)動(dòng)機(jī)試車(chē)出現(xiàn)故障時(shí)，頻譜中位于970～1 000 Hz之間的分頻最為突出；沒(méi)有出現(xiàn)故障時(shí)，頻譜中發(fā)動(dòng)機(jī)1 010～1 030 Hz之間的分頻(發(fā)動(dòng)機(jī)主轉(zhuǎn)速某倍頻)占主導(dǎo)，而970～1 000 Hz之間的分頻沒(méi)有明顯表現(xiàn)。

圖2 某結(jié)構(gòu)焊縫開(kāi)裂Fig.2 The weld-line crack of engine component

圖3 試車(chē)過(guò)程燃料泄漏Fig.3 The leakage of engine fuel

圖4 穩(wěn)態(tài)段頻譜Fig.4 The spectrum of vibration stationary-signal

然而，并不是只要出現(xiàn)970～1 000 Hz之間的分頻發(fā)動(dòng)機(jī)就必然出現(xiàn)故障。圖5是16次試車(chē)970～1 000 Hz之間的分頻的最大幅值對(duì)比，可見(jiàn)：其中有5次試車(chē)分頻幅值明顯高于其他次試車(chē)，但實(shí)際試車(chē)只有其中3次試車(chē)出現(xiàn)了結(jié)構(gòu)局部損傷破壞等故障。

圖5 分頻幅值對(duì)比Fig.5 Frequency-amplitude comparison

2.2 振動(dòng)全頻段RMS值分析

該發(fā)動(dòng)機(jī)試車(chē)速變信號(hào)采樣頻率12.56 kHz，分析最大頻率5 kHz。則某頻段內(nèi)的RMS值為

(1)

式中Ai為fi對(duì)應(yīng)的幅值，i=1，2，…,k-1。

該參數(shù)16次試車(chē)全頻段RMS值對(duì)比分析如圖6所示，可見(jiàn)，從全頻段無(wú)法識(shí)別出故障試車(chē)。這是因?yàn)槿l段的振動(dòng)綜合值既包含了故障信號(hào)又含有正常信號(hào)的貢獻(xiàn)量。

圖6 全頻段RMS值對(duì)比Fig.6 RMS of the full-band comparison

3 信號(hào)故障特征提取

關(guān)注特征頻段970～1 000 Hz內(nèi)RMS值，16次試車(chē)對(duì)比分析如圖7所示，可見(jiàn)1，7，11次特征頻段內(nèi)的RMS值遠(yuǎn)大于其他次試車(chē)，而這3次試車(chē)恰好為出現(xiàn)局部破壞等故障的試車(chē)，因此，特征頻段內(nèi)的RMS值能夠很好地識(shí)別出故障試車(chē)。

圖7 特征頻段RMS值對(duì)比Fig.7 RMS of the characteristic frequency band for stationary-signal

對(duì)特征全頻段進(jìn)行全程瀑布圖細(xì)化分析，如圖8～圖10所示。由圖8可知，對(duì)于大部分沒(méi)有出現(xiàn)任何故障的試車(chē)，其信號(hào)頻域在起動(dòng)時(shí)刻也會(huì)出現(xiàn)970～1 000 Hz之間的明顯分頻，但起動(dòng)結(jié)束后，此分頻消失。由圖9可知，當(dāng)發(fā)動(dòng)機(jī)試車(chē)出現(xiàn)故障時(shí)，其頻域必然出現(xiàn)明顯的介于970～1 000 Hz之間的分頻，且其幅值隨著時(shí)間持續(xù)增大。由圖10可知：對(duì)于兩次出現(xiàn)明顯的幅值較大的特征頻率但沒(méi)有發(fā)生局部損壞的試車(chē)，其瀑布圖表現(xiàn)為幅值較大的970～1 000 Hz之間的分頻在150 s左右，幅值明顯減小。

圖8 正常試車(chē)特征頻段全程瀑布圖Fig.8 The full-time domain waterfall of the characteristic frequency band for normal tests

4 故障分析及預(yù)防

通過(guò)上述分析可見(jiàn)：發(fā)動(dòng)機(jī)試車(chē)速變數(shù)據(jù)包含著大量有用信息，是評(píng)估發(fā)動(dòng)機(jī)運(yùn)行狀態(tài)的不可缺少的組成部分。發(fā)動(dòng)機(jī)試車(chē)速變信號(hào)的分析不僅要關(guān)注信號(hào)的時(shí)域信息，更要關(guān)注其頻域信息，不僅要關(guān)注信號(hào)的振動(dòng)量級(jí)還要關(guān)注信號(hào)的頻率成份，不僅要關(guān)注頻率成份的幅值，更要關(guān)注幅值隨時(shí)間的變化情況，只有這樣，才能全方位、多角度對(duì)信號(hào)進(jìn)行深層次的挖掘和分析。

圖9 出現(xiàn)故障的試車(chē)特征頻段全程瀑布圖 Fig.9 The full-time domain waterfall of the characteristic frequency band for unnormal tests

圖10 某次試車(chē)特征頻段全程瀑布圖Fig.10 The full-time domain waterfall of the charac teristic frequency band for a certain test

該型號(hào)發(fā)動(dòng)機(jī)某次飛行數(shù)據(jù)也出現(xiàn)了量級(jí)較大的、處于970～1 000 Hz之間的分頻，且其飛行殘骸斷口處呈現(xiàn)疲勞條紋，另外通過(guò)計(jì)算發(fā)現(xiàn)該發(fā)動(dòng)機(jī)組件的一階縱向聲學(xué)頻率與該分頻相接近。結(jié)合試車(chē)速變數(shù)據(jù)分析結(jié)果，顯然，發(fā)動(dòng)機(jī)試車(chē)過(guò)程中存在的970～1 000 Hz之間的分頻是導(dǎo)致發(fā)動(dòng)機(jī)結(jié)構(gòu)破壞主要原因。當(dāng)試車(chē)信號(hào)出現(xiàn)此分頻且其幅值隨時(shí)間或者工況不斷增大時(shí)，發(fā)動(dòng)機(jī)結(jié)構(gòu)由于受到持續(xù)的、大量級(jí)的耦合共振而出現(xiàn)疲勞裂紋，從而會(huì)導(dǎo)致發(fā)動(dòng)機(jī)結(jié)構(gòu)焊縫開(kāi)裂或是燃料泄漏等故障。

隨著現(xiàn)代傳感器技術(shù)及計(jì)算機(jī)應(yīng)用技術(shù)的高度發(fā)展，在多個(gè)行業(yè)已經(jīng)實(shí)現(xiàn)對(duì)試驗(yàn)過(guò)程的測(cè)量信號(hào)進(jìn)行實(shí)時(shí)分析和監(jiān)測(cè)[5-10]。針對(duì)該發(fā)動(dòng)機(jī)試車(chē)信號(hào)故障特征，可以在發(fā)動(dòng)機(jī)試車(chē)過(guò)程中對(duì)信號(hào)特征頻段的RMS值進(jìn)行實(shí)時(shí)監(jiān)控，在發(fā)動(dòng)機(jī)試車(chē)至150 s左右時(shí)(工況調(diào)整)，一旦發(fā)現(xiàn)監(jiān)測(cè)對(duì)象沒(méi)有衰減反而有增大趨勢(shì)，可以馬上實(shí)施手動(dòng)緊急關(guān)機(jī)，故障進(jìn)一步發(fā)展惡化，對(duì)試車(chē)產(chǎn)品及試車(chē)臺(tái)造成更大的破壞和損失。

5 結(jié)語(yǔ)

給出了液體火箭發(fā)動(dòng)機(jī)信號(hào)分析方法，對(duì)某型號(hào)發(fā)動(dòng)機(jī)16次試車(chē)振動(dòng)數(shù)據(jù)進(jìn)行對(duì)比分析，找到了正常的試車(chē)與出現(xiàn)故障的試車(chē)的數(shù)據(jù)的異同，提取到了故障識(shí)別的特征量，提出了一種基于特征段RMS值的故障識(shí)別方法，此方法可以用于發(fā)動(dòng)機(jī)試車(chē)過(guò)程故障的實(shí)時(shí)監(jiān)測(cè)，對(duì)于確保試車(chē)產(chǎn)品、預(yù)防試車(chē)破壞性災(zāi)難具有一定的實(shí)際意義。