基于數據的發動機模型匹配方法研究

薛文鵬，段小維，潘鵬飛

（中國飛行試驗研究院,西安710089）

基于數據的發動機模型匹配方法研究

薛文鵬，段小維，潘鵬飛

（中國飛行試驗研究院,西安710089）

針對在使用環境條件下發動機部件特性未知的問題，為獲得在實際裝機條件下的發動機部件特性，采用1種基于參考數據的發動機部件級模型匹配方法，在對測量數據和模型特征分析的基礎上，選取適當調整參數，以模型在非設計點的仿真輸出與參考數據的匹配精度為目標，通過迭代方法求解部件特性。仿真數據表明：采用基于數據的發動機匹配技術得到的發動機模型，其仿真輸出與參考數據的偏差均在允許范圍內。利用該方法可以得到在使用條件下的發動機部件特性，為裝機狀態發動機的仿真預測提供技術支撐。

模型匹配；部件特性；動態過程；航空發動機

0 引言

航空發動機性能仿真、故障診斷和健康管理高度依賴于發動機部件特性，在非設計狀態下，部件特性的形狀會嚴重影響發動機模型的輸出結果。通常，部件特性是發動機生產廠商通過部件試驗或者在已知部件幾何尺寸的情況下通過CFD仿真得到，在大部分情況下，通過部件試驗獲得的部件特性，由于試驗環境、測試方法、進出口流場等差異，使部件試驗的特性與在發動機上的裝機特性存在差異；由于發動機制造、裝配等原因，發動機部件特性與部件試驗獲得的特性之間存在差異；隨著發動機使用時間的增加，各部件特性發生退化，使得部件特性與試驗特性存在差異。綜上所述，對于發動機試驗和使用者而言，獲得準確的部件特性相對難度大。為了提高發動機仿真預測精度，需根據實際測量參數不斷調整通用部件參數，以反推、更新發動機部件特性，使得模型仿真結果與試驗測量結果相匹配。

國外Saravanamuttoo采用通用部件特性完成發動機非設計點仿真計算[1]；Kurzke采用輔助線法解決部件特性不匹配問題[2]；Kong et al.采用最優化算法使得部件特性與測量數據相匹配，但其要求初始的部件特性與測試發動機特性非常接近[3]；Elias Tsoutsanis采用二次曲線的形式表示部件特性，進而對發動機部件級模型進行匹配，獲得比較理想的結果[4]；李秋紅提出了1種基于變適應度函數的模型優化算法，使各截面的進、出口參數與部件特性直接關聯，避免了局部特性修正的誤差積累[5]。

為了提高發動機仿真和預測的準確性和可靠性，必須通過發動機整機試驗數據反推發動機的部件特性。本文通過最優化方法調整發動機部件特性，使得模型的仿真輸出與試驗測量結果相吻合，從而建立發動機穩態仿真模型。

1 模型匹配技術

1.1 穩態過程

發動機處于穩態時，壓氣機耗功與渦輪功平衡，發動機轉子加速率為零，即

式中：N為發動機轉子轉速；W為轉子功。

1.2 過渡態過程

當發動機工作在過渡狀態時，其轉子具有不平衡扭矩，從而實現發動機加、減速過程，如圖1所示。此時，發動機轉子加速率不再為零，表示為

其中ΔΝ采用向后差分方式。

1.3 模型匹配技術

發動機模型匹配過程是在已建立發動機部件模型基礎上，在環境條件和輸入相同的前提下，通過不斷調整模型部件特性，使得模型的仿真輸出與試驗結果相吻合的過程，其結構如圖2所示。模型匹配主要是將發動機部件特性匹配發動機測量數據，如壓氣機和渦輪壓比（πc，πt）、換算流量（mc,mt）和等熵效率（ηc,ηt）。在發動機模型匹配過程中，發動機的輸入輸出關系為

式中：Y為發動機各截面測量參數，包括截面溫度、壓力以及發動機推力等，即Y=[P,T,…FN]；X為發動機的部件特性（流量、壓比和效率），X=[m,π,η]；u 為發動機工作的環境條件和發動機狀態，可選取u=[Ta,Pa,N,Wf]。

當發動機處于過渡態時，壓氣機耗功與渦輪功不平衡，不平衡功率使得轉子加、減速，因而，轉子的加速率還可表示為

式中：J為轉子轉動慣量；N為轉子轉速，r/min；ΔW表示為

式中：Wc、Wt分別為壓氣機功和渦輪功。

在發動機模型匹配過程中，首先確保部件特性X對應的發動機部件模型的共同工作方程ψ=0成立，同時要求模型輸出Y與試驗數據Yr差值最小。因此取目標函數J為

當J取最小值時，X為u對應的部件特性。

1.4 壓氣機模型匹配

壓氣機部件特性包括[mc,ηc,πc]，壓氣機進、出口可直接測量參數包括壓氣機進出口溫度 [T2Ref,T3Ref]和壓力[P2Ref,P3Ref]。

根據測量參數，可直接計算的部件特性包括[ηc,πc]。而流經壓氣機的空氣流量需根據其它測量參數迭代計算確定，在進行發動機模型匹配時，可根據需要給定流量初始值mc。其中壓氣機壓比πc可通過為壓氣機出口壓力與進口壓力的比值計算，即

可通過迭代算法選取壓氣機效率ηc為局部調整變量，使得壓氣機出口仿真溫度與測量溫度的差值最小，即目標函數

另外，壓氣機效率ηc也可通過熱力學方法計算。計算過程可參考相關資料，在此不再敘述。

1.5 渦輪模型匹配

渦輪部件特性包括[mt,ηt,πt]，由于渦輪工作環境溫度高，參數測量難度大，因此，渦輪部件特性計算沒有可依賴的測量數據。但按照發動機流程，根據發動機工作原理，可間接得到渦輪部件特性[mt,πt]，而渦輪效率可通過迭代計算方法確定，在進行渦輪計算時，可給定效率初始值ηt。按照發動機流程可計算渦輪進口截面參數[T41,P41,m41]，則渦輪流量為

當發動機工作處于過渡態時，由式（4）可得渦輪功表達式

其中ΔW可由式（3）得到，即

氣流經過渦輪焓值變化為

當發動機處于穩態時，由式（1）、（4）可得渦輪功的表達式

根據渦輪焓值變化ΔHt、進口參數和給定效率迭代值，利用熱力學關系計算渦輪壓比，得到渦輪出口截面參數。

綜上分析，在進行部件模型匹配時，需取壓氣機流量和渦輪效率作為迭代變量，即X=[mc,ηt]，在壓氣機部件匹配時，取壓氣機效率ηc為壓氣機部件局部迭代變量。

2 仿真實例

選用單軸渦輪噴氣發動機作為仿真對象，進行發動機模型匹配研究。采用2個發動機模型進行發動機模型的匹配，運用Gasturb-10工業軟件仿真生成參考數據。采用上述方法，利用參考數據進行模型匹配。根據上述分析，在發動機模型的匹配過程中，選取調整參數 X=[mc,ηt]。模型輸入參數為 u=[Ta,Pa,N,Wf]，模型仿真輸出和參考數據Y=[P3,T3,…FN]。輸出參數和調整參數的選擇根據發動機參數敏感性分析，如圖3所示。輸出參數的變化應對調整參數的變化敏感，在調整參數進行1%變化時，輸出參數的變化應大于1%或至少接近1%。對圖3分析可知，對于上述發動機輸出參數為 Y=[P3，T3，P5，T5，FN]。

給定調整參數初始值[mc0,ηt0]，采用全局尋優算法進行部件特性調整，匹配結果為換算轉速及其對應的部件特性，即 y=[ncor,mc,πc,ηc]，模型匹配過程結構如圖4所示。

2.1 穩態過程匹配

采用上述方法，令ΔW=0,dN/dt=0，對發動機進行穩態模型匹配，匹配結果如圖5～8所示。從圖5中可見模型匹配過程中目標函數的變化，隨著迭代次數的增大，目標函數趨向于極小；從圖6、7中可見模型匹配結果與參考模型的對比，壓氣機部件特性（流量、壓比、效率）與參考模型特性相一致；從圖8中可見，發動機測量截面參數誤差，壓氣機出口截面參數仿真誤差較小，按照發動機流程，仿真誤差不斷累積，發動機推力預測誤差相對較大。

對于發動機不可直接測量參數（燃燒室出口溫度）的匹配偏差和匹配前、后對比結果如圖9所示。匹配的最大偏差不大于2%。

2.2 動態過程匹配

發動機某過渡態過程的供油規律如圖10所示。采用上述匹配技術進行模型匹配，匹配結果如圖11～14所示。

從圖11、12中可見模型匹配輸出y與參考模型部件特性的比較，壓氣機部件特性（流量、壓比、效率）的變化規律與參考模型一致，其中流量最大偏差為3.2×10-4；壓比的匹配偏差相對較小；從圖13中可見模型匹配得到的發動機加速率與參考模型的對比結果，加速率計算誤差均滿足精度要求；從圖14中可見發動機各測量截面參數Y的匹配結果，圖中發動機推力和各界面參數誤差均在10-3范圍內。對于非測量參數（如燃燒室出口溫度）如圖15所示，匹配偏差在10-4范圍內。

采用上述方法對發動機模型進行匹配，匹配前、后發動機部件特性的比較如圖16所示。

針對渦輪噴氣發動機的仿真實例表明：利用基于數據的發動機模型匹配技術可以還原發動機部件特性，同時實現發動機各截面參數（包括不可測量參數）的仿真、預測。

3 總結

（1）介紹了1種發動機模型匹配方法（包括穩態和動態過程）。是對發動機數據模型匹配技術的初步探索，按照發動機工作流程，根據可直接測量參數，對發動機模型匹配過程進行分析，選取調整參數和目標函數。采用迭代優化算法對匹配參數進行全局尋優，使得模型仿真輸出與試驗數據差值的目標函數最小。

（2）根據建立的單軸渦噴發動機部件級模型，利用數據對部件級模型進行模型匹配。結果顯示部件特性的匹配取得理想結果，匹配模型的輸出和參考數據相吻合。對于不可測量參數的仿真結果與參考數據的誤差在允許范圍內。

（3）文中僅對單軸發動機實現模型匹配過程，但提供的模型匹配方法具有一定的通用性，可根據實際發動機類型和測量參數選取適當的調整參數和目標函數，實現其它模型匹配。

[1]Zhu P,Saravanamuttoo H I H.Simulation of an advanced twin-spool industrial gas turbine[J].Journal of Engineering for Gas Turbines&Power,1992,114（2）:180-186.

[2]Kurzke J.How to get component maps for aircraft gas turbine performance calculations[C]//ASME 1996 International Gas Turbine and Aeroengine Congress and Exhibition.American Society of Mechanical Engineers,1996.

[3]Kong C,Ki J,Kang M.A new scaling method for component maps of gas turbine using system dentification[J].Journal of Engineering for Gas turbines and Power,2003,125（4）:979-985.

[4]Li Y G,Pilidis P,Newby M A.An adaptation approach for gas turbine design-point performance simulation[J].Journal of Engineering for Gas Turbines and Power,2006,128（4）:789-795.

[5]Tsoutsanis E,Meskin N,Benammar M,et al.An efficient component map generation method for prediction of gas turbine performance[C]//ASME Turbo Expo.2014.

[6]潘鵬飛,李秋紅,任冰濤,等.基于遺傳算法的航空發動機部件特性修正[J].北京航空航天大學學報,2014（5）:690-694.PAN Pengfei,LI Qiuhong,REN Bingtao，et al.Research on aeroengine design-point performance adaptation using a genetic algorithm[J].Journal of Beijing University of Aeronautics and Astronautics,2014（5）:690-694.（in Chinese）

[7]Tsoutsanis E,Meskin N,Benammar M,et al.Transient gas turbine performance diagnostics through nonlinear adaptation of compressor and turbine maps[J].Journal of Engineering for Gas Turbines&Power,2015,137（9）：1-2.

[8]Sethi V,Doulgeris G,Pilidis P,et al.The map fitting tool methodology:gas turbine compressor off-design performance modeling[J].Journal of Turbomachinery,2013,135（6）:522-535.

[9]Miste G A,Benini E.Turbojet engine performance tuning with a new map adaptation concept[J].Journal of Engineering for Gas Turbines&Power,2014,136（7）:V001T06A002.

[10]Li Y G.Aero gas turbine flight performance estimation using engine gas path measurements[J].Journal of Propulsion and Power,2015,31（3）:851-860.

[11]Ghorbanian K,Gholamrezaei M.An artificial neural network approach to compressor performance prediction[J].Applied Energy,2009,86（7-8）:1210-1221.

[12]Kong C,Ki J.Components map generation of gas turbine engine using genetic algorithms and engine performance deck data[J].Journal of Engineering for Gas Turbines&Power,2006,129（2）:377-383.

[13]Yoshida T,Morikawa T,Gotou J,et al.Gas turbine performance analysis method and gas turbine performance analysis system:US 20090055105[P].2010-9-14.

[14]Oyedepo S O,Fagbenle R O,Adefila S S,et al.Performance evaluation and economic analysis of a gas turbine power plant in Nigeria[J].Energy Conversion&Management,2014,79（3）:431-440.

[15]Ghisu T,Parks G T,Jarrett J P,et al.Adaptive polynomial chaos for gas turbine compression systems performance analysis[J].Aiaa Journal,2015,48（6）:1156-1170.

Research on Matching Method of Engine Model based on Test Data

XUE Wen-peng DUAN Xiao-wei PAN Peng-fei
（China Flight Test Establishment,Xi'an 710089,China）

In view of the problem of unknown components at service condition,an adapting method based on the rig data was proposed to obtain the characteristics of the engine components under the actual condition.Adjustment parameters were selected based on the analysis of model characteristics and test data.The error of model simulation and test data was selected as objective function and used in iterative algorithm for component characteristics.The simulation results show that the error between the model simulation which derived from the method and test data is within the maximum error.The components characteristics was obtained by using the engine matching technology based on the test data,which provide support for engine simulation and prediction under the actual condition.

model matching；component characteristics；transition state；aeroengine

V 231.1

A

10.13477/j.cnki.aeroengine.2017.05.017

2017-04-09

薛文鵬（1987），男，碩士，工程師，從事航空發動機整機試驗技術研究工作；E-mail:xwpfenyun@163.com。

薛文鵬，段小維，潘鵬飛.基于數據的航空發動機模型匹配方法研究[J].航空發動機，2017，43（5）：97-102.XUE Wenpeng，DUAN Xiaowei,PAN Pengfei.Research on matching method of engine model based on test data[J].Aeroengine，2017，43（5）：97-102.

（編輯：李華文）