姜榮俊，高建華，劉永葆

艦船燃?xì)廨啓C氣路測量參數(shù)的理論選擇方法

姜榮俊，高建華，劉永葆

（海軍工程大學(xué)動力工程學(xué)院，武漢430033）

針對艦船燃?xì)廨啓C氣路性能監(jiān)測診斷的有效性要求，歸納總結(jié)了氣路測量參數(shù)選擇要求，提出了1套系統(tǒng)完整的燃?xì)廨啓C氣路測量參數(shù)的理論選擇方法。在可測量和易測量要求及測試精度要求的初步選擇基礎(chǔ)上，依次采用影響系數(shù)、相關(guān)系數(shù)和條件數(shù)對氣路測量參數(shù)的敏感性、相關(guān)性和診斷誤差等進行了分析，最終從理論上選擇了某艦船3軸燃?xì)廨啓C氣路性能監(jiān)測診斷的測量參數(shù)。實例表明該理論方法具有一定的通用性和工程應(yīng)用價值，對燃?xì)廨啓C氣路測量參數(shù)的選擇有一定的指導(dǎo)作用。

氣路；測量；監(jiān)測；診斷；艦船；燃?xì)廨啓C

0 引言

發(fā)動機監(jiān)測診斷是保證發(fā)動機安全運行、減少壽命費用、防止惡性事故的有效手段。進行燃?xì)廨啓C監(jiān)測診斷必須根據(jù)實際目標(biāo)要求首先對其測量參數(shù)進行優(yōu)化選擇[1]。針對燃?xì)廨啓C的氣路性能監(jiān)測診斷，Urban早在20世紀(jì)70年代就提出了氣路分析（gas path analysis-GPA）技術(shù)[2]，其原理是基于燃?xì)廨啓C各部件間的氣動熱力學(xué)關(guān)系式建立小偏差故障模型，構(gòu)建燃?xì)廨啓C性能變化量與測量參數(shù)變化量的函數(shù)關(guān)系，從而實現(xiàn)利用帶有噪聲或偏置的氣路參數(shù)對燃?xì)廨啓C氣路性能進行監(jiān)測診斷。

為了較好地對燃?xì)廨啓C的氣路性能實施監(jiān)測和診斷，基于各種要求和考慮，國內(nèi)外學(xué)者對燃?xì)廨啓C氣路測量參數(shù)的選擇進行了研究。Urban基于線性氣路分析技術(shù)（LGPA—linear gas path analysis），提出采用診斷誤差作為篩選測量參數(shù)的1個重要準(zhǔn)則，并明確提出測量參數(shù)的數(shù)目不能少于性能參數(shù)數(shù)目[3]；Consumi M和Agostino L d基于數(shù)據(jù)統(tǒng)計推理，根據(jù)不同的氣路測量參數(shù)組合的診斷精度對測量參數(shù)進行了選擇[4]；Ogaji SOT等基于敏感度和診斷誤差，分別利用線性和非線性氣路分析方法（NLGPA—non-linear gas path analysis），對燃?xì)廨啓C測量參數(shù)的選擇問題做了進一步研究，結(jié)果表明在模擬單一氣路故障情況下，根據(jù)NLGPA方法選擇出的測量參數(shù)方案更為合理；而在多種氣路故障情況下，NLGPA方法的優(yōu)越性并不十分明顯[5]；唐耿林基于敏感性、相關(guān)性和診斷誤差對航空發(fā)動機的性能測量參數(shù)進行了選擇[6]；孫祥逢等針對氣路故障診斷主因子模型，基于敏感性、相關(guān)性和診斷誤差對測量參數(shù)進行了選擇，證明測量參數(shù)太多對采用故障方程的診斷不利，并嘗試給出了解釋[7]。

本文主要針對監(jiān)測診斷結(jié)果的有效性，從理論上對某艦船3軸燃?xì)廨啓C的氣路監(jiān)測診斷參數(shù)進行了優(yōu)化選擇，對于工程中的安全性、可靠性等因素則不予太多考慮。

1 氣路測量參數(shù)的初步選擇

某艦船3軸燃?xì)廨啓C為帶雙轉(zhuǎn)子燃?xì)獍l(fā)生器的簡單開式循環(huán)燃?xì)廨啓C，其具體結(jié)構(gòu)如圖1所示，圖中標(biāo)識出其氣路通道的各特征截面?？紤]到穩(wěn)態(tài)過程的氣路容積效基本可以忽略，因此對于同一特征截面處的上1個部件出口和下1個部件的進口，可以認(rèn)為其氣路參數(shù)相同。

圖1 3軸燃?xì)廨啓C結(jié)構(gòu)

燃?xì)廨啓C的工作狀態(tài)可用一定的特征量來描述，由這些特征量可以確定燃?xì)廨啓C的完好性、工作能力和其動作的正確性。在燃?xì)廨啓C氣路監(jiān)測診斷中，特征量通常指需要分析的部件特性參數(shù)或幾何參數(shù)，如壓氣機效率、渦輪導(dǎo)向器面積等。在實際運行的不解體情況下，特征量通常是不可測量的。從特征量可以引出描述故障狀態(tài)的故障因子[8]。而測量參數(shù)主要指可觀測的且能區(qū)分發(fā)動機技術(shù)狀態(tài)的狀態(tài)量，這些參數(shù)主要是指氣路各截面上能直接測量的參數(shù)（圖1），如轉(zhuǎn)子轉(zhuǎn)速、溫度、壓力等[9]。燃?xì)廨啓C的氣路性能監(jiān)測診斷主要是通過對其氣路通道中壓力、溫度和轉(zhuǎn)速等狀態(tài)量的監(jiān)測，實現(xiàn)對壓氣機、渦輪和燃燒室等部件的效率、流量等特性參數(shù)或幾何參數(shù)特征量的診斷。

考慮到艦船燃?xì)廨啓C中動力渦輪的氣路性能故障通常較少，且其氣路參數(shù)受負(fù)載（如螺旋槳）的影響較大，變化較復(fù)雜，因此一般不作為氣路性能監(jiān)測參數(shù)，但為了防止超速，仍需對其轉(zhuǎn)速進行監(jiān)控。由于燃燒室掉塊、燒裂等故障不會引起發(fā)動機性能參數(shù)的明顯變化[10]，因此只分析與燃油噴嘴堵塞故障有關(guān)的燃燒室效率特征量。選取9個需要監(jiān)測診斷的特征量：高、低壓壓氣機的效率η和流量W；高、低壓渦輪及動力渦輪的效率和流量；燃燒室效率。

艦船3軸燃?xì)廨啓C，在一定的環(huán)境參數(shù)（大氣壓力、溫度和濕度）和油門位置下，可以被用作氣路監(jiān)測診斷的參數(shù)主要有以下幾種：（1）轉(zhuǎn)速：包括高、低壓壓氣機轉(zhuǎn)速；（2）壓力：包括高、低壓壓氣機，燃燒室，高、低壓渦輪的出口壓力等；（3）溫度：包括高、低壓壓氣機，燃燒室，高、低壓渦輪的出口溫度等；（4）其它：包括燃油流量、發(fā)動機扭矩或功率等。

首先從可測量和易測量要求出發(fā)來分析以上氣路測量參數(shù)：由于高速、高溫、高壓和氣流影響，燃燒室出口、高壓渦輪出口等處溫度和壓力往往比較難測出且成本較高，即使在臺架試車非全流程參數(shù)測試情況下[11]，這些參數(shù)也往往不測量，通常不作為氣路監(jiān)測診斷參數(shù)。為了檢驗不測量此處參數(shù)對氣路監(jiān)測診斷效果的影響，這里暫時考慮測量燃燒室出口、高壓渦輪出口處溫度和壓力。

對這些參數(shù)的測試精度進行分析，相對臺架試驗條件，實際使用中通常會降低測試精度[11]，但隨著發(fā)動機測控技術(shù)的進步和測試標(biāo)準(zhǔn)的提高，在排除人為誤差的情況下，認(rèn)為其測試精度滿足測量誤差小于故障變化引起的測量參數(shù)偏差的要求。

基于以上步驟，除環(huán)境溫度和壓力、運行工況參數(shù)（油門位置等）的運行環(huán)境參數(shù)外，初步確定了13個氣路測量參數(shù)：低壓壓氣機的出口溫度T1.5、出口壓力P1.5、轉(zhuǎn)速n1，高壓壓氣機的出口溫度T2、出口壓力P2、轉(zhuǎn)速 n2，燃燒室的出口溫度 T3、出口壓力 P3，高壓渦輪的出口溫度T3.5、出口壓力P3.5，低壓渦輪的出口溫度T4、出口壓力P4，發(fā)動機功率N。這13個參數(shù)是否能夠滿足敏感性、相關(guān)性和診斷誤差等測量參數(shù)要求，還需進一步分析。在繼續(xù)分析之前，首先對所采用的GPA方法進行介紹。

2 氣路分析法介紹

GPA的原理是映射燃?xì)廨啓C測量參數(shù)和特征量間的氣動熱力學(xué)關(guān)系[2,13]，其一般形式的數(shù)學(xué)模型為

式中：輸入矢量P∈Rp，為運行環(huán)境參數(shù)；P∈Rn，為特征量參數(shù)；輸出矢量Y∈Rm，為氣路性能測量參數(shù)。

某基準(zhǔn)運行工況（用下標(biāo)“0”表示），式（1）利用泰勒級數(shù)展開

式中：HOTs為偏差的高階項，在小偏差的情況下忽略其影響，成為線性化方程。若運行環(huán)境參數(shù)與基準(zhǔn)環(huán)境參數(shù)相同，則（P-P0）=0。式（2）簡化為

以小偏差量來表示

A∈Rm×n，是影響系數(shù)矩陣（ICM-influence coefficient matrix），是 F（X，P） 關(guān)于特征量 X 的Jacobian矩陣

式（4）即經(jīng)典的GPA小偏差方程，它反映了測量參數(shù)和特征量參數(shù)的關(guān)系。若矩陣A可逆，則得

式中：A-1為故障系數(shù)矩陣（FCM-fault coefficient matrix）?；谠摼€性化的關(guān)系式，通過測量參數(shù)變化量可得到燃?xì)廨啓C部件特征量參數(shù)的變化量，從而實現(xiàn)對燃?xì)廨啓C氣路性能的監(jiān)測診斷。

實踐中往往很難得到實際發(fā)動機模型F（X，P）的精確解析形式，因此現(xiàn)在通常采用精確仿真模型的數(shù)值線性化的方法來得到影響系數(shù)矩陣A。且A通常為非正定矩陣，因此實際中很少直接采用式（6）來求解式（4），更多采用最小二乘等方法。

如果考慮測量誤差因素ν，則式（4）為

以上是Urban早在1960年提出的基于小偏差方程的LGPA[14]。

3 氣路測量參數(shù)的敏感性選擇

對特征量敏感的測量參數(shù)才能暴露故障，因此可根據(jù)敏感性原則直接選擇測量參數(shù)。表征系統(tǒng)的特征量的各狀態(tài)量，都可以作為對系統(tǒng)進行技術(shù)監(jiān)測診斷的主要依據(jù)。每個狀態(tài)量對系統(tǒng)的故障部位和故障程度的敏感性不同，可以從這些參數(shù)中挑選出對故障特征量最敏感的狀態(tài)量，也可以用某幾個狀態(tài)量組合成1個敏感性較大的綜合參數(shù)。氣路性能測量參數(shù)可通過測量與特定氣路故障有關(guān)的狀態(tài)量來直覺猜測，但基于GPA的方法往往更可取。根據(jù)GPA，利用小偏差原理，可以計算出特征量偏差和測量參數(shù)偏差之間的影響系數(shù)矩陣A?；谌?xì)廨啓C仿真模型，通過數(shù)值計算，可以得到某型艦船燃?xì)廨啓C在1.0設(shè)計工況下的影響系數(shù)矩陣，具體數(shù)值見表1，該系數(shù)矩陣反應(yīng)了測量參數(shù)對特征量的敏感性。測量參數(shù)對不同特征量變化的2種敏感性如圖2所示。在描述敏感程度時，依據(jù)影響系數(shù)的絕對值大小，可分為：較敏感（≥0.50）和非常敏感（≥1.00）。

表1 某型艦船燃?xì)廨啓C1.0設(shè)計工況下的影響系數(shù)矩陣

圖2 不同特征量變化下測量參數(shù)的敏感性

（1）低壓壓氣機對 δWCL，n1、P1.5較敏感；對 δηCL，P1.5、P2、P3.5、P3、N、P4、T4、T3.5較敏感。T1.5對包括低壓壓氣機在內(nèi)的性能特征量變化都不敏感，但對低壓渦輪的特征量變化有比較突出的表現(xiàn)（絕對值仍較小），因此從低壓渦輪氣路性能的監(jiān)測診斷需要出發(fā)，是否測量T1.5值得進一步分析。

（2）在高壓壓氣機中，對于δWCH，幾乎所有的參數(shù)都不很敏感。壓氣機葉片結(jié)垢、葉頂間隙增大會導(dǎo)致流量特征量變化，考慮故障監(jiān)測診斷的需要，暫時保留相對較敏感的幾個參數(shù)：n2、P2、P3、N（按敏感度從大到小排序，取影響系數(shù)絕對值≥0.25）；對于δηCH，N、P2、P3、P3.5、P4、T4、T3.5、T3都較敏感。T2對包括高壓壓氣機自身在內(nèi)的特征量變化也不敏感，所以不測量該參數(shù)。

（3）在高壓渦輪中，對于 δWTH，P3、P2、N 較敏感；對于 δηTH，P2、P3、N、P3.5、P4、T4、T3.5都較敏感。結(jié)合前述初步選擇時對高溫截面處測量參數(shù)的選擇，確定采用溫度相對較低且對高壓渦輪特征量變化也非常敏感的 N、P2、P4、T4，而不測量 T3.5和 P3.5。

（4）在低壓渦輪中，對于 δWTL，P3.5、P1.5、T1.5較敏感；對于 δηTL，P1.5、P2、P3.5、P3、N、T4都較敏感。結(jié)合前述分析，選擇不測量 P3、P3.5。

（5）在燃燒室中，對于 δηB，N 非常敏感，P3.5、P3、P2、P4、P1.5也較敏感，根據(jù)前述分析，選擇舍棄 P3.5、P3。通過敏感性分析，選擇了 7 個測量參數(shù)：P1.5、n1、P2、n2、T4、P4、N。

4 氣路測量參數(shù)的相關(guān)性選擇

相關(guān)是指變量之間的線性關(guān)系。參數(shù)的相關(guān)性是指在不同故障作用下測量參數(shù)組變化之間的線性關(guān)系。對測量參數(shù)的相關(guān)性要求是：對不同的故障特征量，測量參數(shù)有不同的變化特征。即測量參數(shù)之間的相關(guān)性要盡量小，使所選擇的測量參數(shù)不僅能夠表現(xiàn)各部件特征量的變化，且對于不同的故障，測量參數(shù)能表現(xiàn)出不同的變化特征，達(dá)到隔離故障的目的[7]。相關(guān)性的強弱常用相關(guān)系數(shù)γxy表示，其值越接近于0，表示相關(guān)性越弱。另文獻[6]提出：若某2個不同故障對同1組測量參數(shù)的影響相當(dāng)，則該組參數(shù)對這2個故障是相關(guān)的，而某2個故障對同1組測量參數(shù)的影響方向相反，則不相關(guān)。因此相關(guān)性分析可以通過ICM中的正負(fù)符號進行判斷。

綜合這2種測量參數(shù)相關(guān)性判斷方法，本文提出：

（1）若γxy接近0，則測量參數(shù)相關(guān)性?。?/p>

（2）若 γxy＜0（負(fù)相關(guān)），則測量參數(shù)不相關(guān)；

（3）若γxy接近1，且能找到變化相反并足夠敏感的測量參數(shù)，也認(rèn)為測量參數(shù)不相關(guān)。

根據(jù)上面提出的3種判斷方法，可以對氣路測量參數(shù)進行相關(guān)性選擇。

基于式（8），不同特征量下7個測量參數(shù)變化間的相關(guān)系數(shù)值見表2。

表2 不同特征量下測量參數(shù)變化間的相關(guān)系數(shù)值

從表2中可見：

（1）δWCL與 δWTL、δηCL與 δηTL、δηCH與 δηTH這 3 組特征量對應(yīng)的測量參數(shù)變化間的的正相關(guān)性較強，因此基于現(xiàn)有的測量參數(shù)較難區(qū)分。

（2）觀察 δWCL與 δWTL組合，P1.5對 δWTL變化非常敏感，其影響系數(shù)值是P1.5對δWCL的影響系數(shù)值的2.744倍，有可能據(jù)此來區(qū)分這2個特征量。另外依據(jù)判斷方法（3）分析，從表1中可以發(fā)現(xiàn)T2的變化方向相反，但注意到它對δWCL和δWTL的變化都不敏感，影響系數(shù)分別為+0.029和-0.146，因此只有在保證足夠測量精度的情況下才有可能區(qū)分這2個特征量的變化，所以是否選擇T2值得考慮。

（3）觀察 δηCL與 δηTL組合，從表 1 可以發(fā)現(xiàn) T1.5變化方向相反，且對δηTL足夠敏感，影響系數(shù)為+0.407，因此選擇T1.5。

（4）觀察 δηCH與 δηTH組合，從表 1 可以發(fā)現(xiàn) T2的變化方向相反，且二者的影響系數(shù)值也足夠敏感，綜合考慮（2），因此確定選擇T2。

通過以上相關(guān)性分析，在原來7個測量參數(shù)基礎(chǔ)上增加了2個測量參數(shù)T1.5和T2。根據(jù)計算，測量參數(shù)增加后，以上3組特征量對應(yīng)的測量參數(shù)變化間的正相關(guān)性都有所減弱，氣路性能監(jiān)測診斷的有效性有所提高。

5 氣路測量參數(shù)的最終選擇

通過以上步驟得到比較滿意的測量參數(shù)組合。測量參數(shù)數(shù)目達(dá)到了氣路監(jiān)測診斷的最低要求，但監(jiān)測診斷的效果還有待繼續(xù)分析。

氣路診斷的過程實際就是基于小偏差方程（4），根據(jù)已知的測量參數(shù)偏差求解未知性能特征量偏差的數(shù)值計算問題。用數(shù)值方法求解線性方程組（4），計算結(jié)果存在誤差，原因主要有2個：一是選用的計算方法，為了保證算法的穩(wěn)定性，采用選主元的高斯消去法，選主元的三角分解法等；二是方程組本身。一般來說，由于存在各種誤差，線性方程組（4）的系數(shù)矩陣A與測量項δY的數(shù)值并不絕對正確[15]。對于這種偏離正確值的擾動導(dǎo)致的計算精度問題，根據(jù)數(shù)值計算理論，可以用矩陣A的條件數(shù)cond(A)來衡量[16]

式中：“||||”是矩陣的某種范數(shù)。

條件數(shù)反映了矩陣A和方程本身的性態(tài)，條件數(shù)小，擾動引起的解的相對誤差?。粭l件數(shù)大，擾動引起的解的相對誤差就可能大。理論上，條件數(shù)確定了解的相對誤差上限，即擾動可引起誤差的最大放大倍數(shù)。

經(jīng)以上步驟確定的小偏差方程（4），影響系數(shù)矩陣A的性態(tài)將對計算結(jié)果有直接影響。對此，Urban指出[5]，如果測量參數(shù)選取較差，將導(dǎo)致故障系數(shù)矩陣（FCM）中的系數(shù)較大，從而導(dǎo)致方程迭代過程不收斂或者導(dǎo)致較大的計算誤差。由此，可以直接通過影響系數(shù)矩陣的條件數(shù)cond(A)來判斷測量參數(shù)誤差對特征量計算精度的影響。

表1、表3列出了主要幾種不同測量參數(shù)組合下影響系數(shù)矩陣的譜范數(shù)||||2條件數(shù)。表中符號“×”表示選擇該參數(shù)。

對表3進行分析：

（1）前面選擇的9個測量參數(shù)組合的條件數(shù)相對7個基本測量參數(shù)的條件數(shù)較大，即增加測量參數(shù)后，區(qū)分了3組相關(guān)的特征量，使故障定位更有效，同時使方程組（4）的性態(tài)變差了，增加了誤診斷的可能性。

表3 不同測量參數(shù)組合的條件數(shù)

（2）測量參數(shù)的多少反映了對發(fā)動機運行狀況的了解程度。從性能監(jiān)測診斷的角度出發(fā)，傾向于合理選擇更多的測量參數(shù)。但增加測量參數(shù)后可能導(dǎo)致方程組（4）的性態(tài)變差，如（1）所述，因此具體效果需深入分析。

（3）適當(dāng)犧牲故障發(fā)現(xiàn)性能，即減少預(yù)診斷的特征量的數(shù)目，可以改善方程組的性態(tài)。對于選擇的9個測量參數(shù)，如果不考慮δηB特征量變化，則計算的條件數(shù)只有90.8，方程組的性態(tài)得到了很大改善。實際δηB特征量變化有關(guān)的燃油噴嘴堵塞故障也可以通過噴嘴前燃油壓力或燃油流量等反映出來。因此氣路監(jiān)測診斷時還應(yīng)綜合采用其它監(jiān)測診斷手段。

（4）適當(dāng)犧牲故障定位性能，即恰當(dāng)減少測量參數(shù)，可提高方程組（4）的性態(tài)。因此在監(jiān)測診斷時，可減少測量參數(shù)的方程組先行判斷，提高故障發(fā)現(xiàn)性能。在出現(xiàn)易混淆特征量變化，即故障定位性能不夠時，再引入更多測量參數(shù)的方程組進行判斷。

（5）方程組的病態(tài)程度是相對的。病態(tài)的影響隨著計算機字長的增長減輕。病態(tài)方程的有效解法可參看相關(guān)文獻。在消去法的基礎(chǔ)上，對已取得的近似解不斷改進，就可以達(dá)到十分準(zhǔn)確的結(jié)果。因此，由于測量參數(shù)組合導(dǎo)致的診斷誤差問題，可以通過改進診斷算法予以彌補。

（6）加入P3.5測量參數(shù)，可以使參數(shù)組合的條件數(shù)減少。因此，如果高溫下壓力測試可行，可以用參數(shù)P3.5代替測試相對困難的N參數(shù)，該組合不僅滿足診斷的最低數(shù)量測量參數(shù)要求，而且組合的條件數(shù)也相對前面所選9個參數(shù)組合的條件數(shù)要小得多，數(shù)值參見表3。

通過以上對不同測量參數(shù)組合下影響系數(shù)矩陣的條件數(shù)分析，可以認(rèn)為前述步驟選擇的9個測量參數(shù)比較符合某3軸艦船燃?xì)廨啓C的氣路監(jiān)測診斷需求。對于條件數(shù)較大的情況，可以采用（3）、（4）或（5）的方法予以改進。該測量參數(shù)組合包含了實際發(fā)動機控制系統(tǒng)現(xiàn)有的全部5個氣路測量參數(shù)，實現(xiàn)了最大程度的共享，所以該選擇更為優(yōu)化合理。

6 總結(jié)

（1）針對燃?xì)廨啓C氣路性能監(jiān)測診斷需求，歸納總結(jié)了選擇氣路測量參數(shù)的6個要求。

（2）重點利用影響系數(shù)、相關(guān)系數(shù)和條件數(shù)，分析了某艦船3軸燃?xì)廨啓C的氣路測量參數(shù)的敏感性、相關(guān)性和診斷誤差，最終選擇了適用于氣路性能監(jiān)測診斷的9個測量參數(shù)。

（3）該理論方法具有一定的通用性和工程應(yīng)用價值，對燃?xì)廨啓C氣路測量參數(shù)的選擇有一定的指導(dǎo)作用。

[1]ИВ.凱巴著.航空燃?xì)鉁u輪發(fā)動機技術(shù)診斷[M].王起鵬,劉德磬，譯.北京:航空工業(yè)出版社，1990：2-5.ИВ.Keba.Aero gas turbine engine technology diagnosis[M].WANG Qipeng,LIU Dexin,translated.Beijing:Aviation Industry Press,1990:2-5.(in Chinese)

[2]Urban L A.Gas path analysis:applied to turbine engine condition monitoring[J].Journal of Aircraft,1973,10（7）：400-406.

[3]Urban L A.Parameter selection for multiple fault diagnosis of gas turbine engines [J].Journal of Engineering for Power，1975，97（2）：225-230.

[4]Consumi M,Agostino L.A statistical inference approach to gas path analysis of a turbofan[R]．AIAA-1998-3551.

[5]Ogaji S O T，Sampath S，Singh R，et al.Parameter selection for diagnosing a gas turbine’s performance deterioration [J].Applied Ener gy，2002，73（1）：25-46.

[6]唐耿林.航空發(fā)動機性能監(jiān)視參數(shù)選擇的研究 [J].推進技術(shù)，1998，19（2）：38-42.TANGGenglin.Investigationson selectingperformance monitoring pa rameters of aeroengine[J].Journal of Propulsion Technology，1998，19（2）：92-96.(in Chinese)

[7]孫祥逢，陳玉春，胡福.發(fā)動機故障診斷主因子模型的測量參數(shù)選擇[J].航空動力學(xué)報，2010，25（1）：129-135.SUN Xiangfeng，CHEN Yuchun，HU Fu.Research on selection of measurement parameters of engine fault diagnosis based on primary factor mode[J].Journal of Aerospace Power，2010，25（1）：129-135.(in Chinese)

[8]范作民，孫春林，白杰.航空發(fā)動機故障診斷導(dǎo)論[M].北京：科學(xué)出版社，2004：134-138.FANZuomin,SUNChunlin,BAIJie.Introduction toaeroengine fault diagnosis[M].Beijing：Science Press，2004：134-138.(in Chinese)

[9]鄧明，金業(yè)壯.航空發(fā)動機故障診斷[M].北京:北京航空航天大學(xué)出版社，2012：95-96.DENG Ming，JIN Yezhuang.Aeroengine fault diagnosis[M].Beijing：Behang University Press，2012：95-96.(in Chinese)

[10]Ogaji S O T.Advanced gas path fault diagnostics for stationary gas turbines[D].Bedford：Cranfield University，2003.

[11]陳益林.航空發(fā)動機試車工藝[M].北京：北京航空航天大學(xué)出版社，2010：212-215.Chen Yilin.Aeroengine test technology[M].Beijing：Beihang University Press,2010：212-215.(in Chinese)

[12]陳衛(wèi)，程禮，李全通，等.航空發(fā)動機監(jiān)控技術(shù)[M].北京:國防工業(yè)出版社，2011：156-159.Chen Wei，Cheng Li，LI Quantong,et al.Aircraft engine monitoring technology[M].Beijing:National Defense Industry Press,2011：156-159.(in Chinese)

[13]LI Yiguang.A GPA diagnostic system for aeroengine applications[J].Journal of Aerospace Power,2007,22（6）：898-908.

[14]LI Yiguang.Performance analysis based gas turbine diagnostics:a review[J].Journal of Power and Energy，2002，216（5）：363-377.

[15]何漢林，魏汝祥，李衛(wèi)軍.數(shù)值分析[M].武漢:湖北科學(xué)技術(shù)出版社，1999：245-246.HE Hanlin，WEI Ruxiang，LI Weijun.Numerical analysis[M].Wuhan：Hubei Scienceand Technology Press，1999：245-246.(in Chinese)

[16]鄧建中，葛仁杰，程正興.計算方法[M].西安:西安交通大學(xué)出版社，1985：50-52.DENG Jianzhong，GE Renjie，CHENG Zhengxing.Calculation method[M].Xi'an：Xi'an Jiaotong University Press，1985：50-52.(in Chinese)

Selection of Gas Path Measurement Parameters for Marine Gas Turbine in Theory

JIANG Rong-jun,GAO Jian-hua,LIU Yong-bao
（College of Power Engineering,Naval University of Engineering,Wuhan 430033,China）

According to the validity demand of gas turbine gas path performance monitoring and diagnosis,requirements for gas path measurement parameters selection were generalized and summarized,and a systematic and complete theoretical method of gas turbine gas path measurement parameter selection was put forward.Based on the preliminary selection with the requirements of measurability,testability and measuring accuracy,the sensitivity,correlation,diagnosis error of the gas path measurement parameters were analyzed in combine with influence coefficient,correlation coefficient and condition number respectively.Finally,the measurement parameters were chosen for the gas path performance monitoring and diagnosis of a marine gas turbine with three shafts.Example shows that the theoretical method has certain universality and engineering application value,and has certain instruction function to the selection of gas path measurement parameters of gas turbine.

gas path;monitoring;measurement;diagnosis；marine;gas turbine

V231.1

A

10.13477/j.cnki.aeroengine.2015.01.015

2013-10-09

姜榮?。?973），男，博士，研究方向為燃?xì)廨啓C的狀態(tài)監(jiān)測和故障診斷；E-mail：jiang_rj@163.com。

姜榮俊,高建華，劉永葆.艦船燃?xì)廨啓C氣路測量參數(shù)的理論選擇方法[J].航空發(fā)動機，2015，41（1）：74-79.JIANGRongjun，GAOJianhua，LIUYongbao.Selection of gaspath measurement parametersfor marinegasturbinein theory[J].Aeroengine，2015,41（1）：74-79.

（編輯：趙明菁）