航空發(fā)動(dòng)機(jī)是航空器飛行的動(dòng)力，是航空器的“心臟”。目前應(yīng)用最為廣泛的航空發(fā)動(dòng)機(jī)為燃?xì)鉁u輪發(fā)動(dòng)機(jī)，它主要由進(jìn)氣口、壓氣機(jī)、燃燒室、渦輪和尾噴管組成。從進(jìn)氣口進(jìn)入的空氣在壓氣機(jī)中被壓縮后進(jìn)入燃燒室，在燃燒室中與噴入的燃油混合燃燒，生成高溫高壓的燃?xì)狻Ｈ細(xì)庠谂蛎涍^(guò)程中驅(qū)動(dòng)渦輪高速旋轉(zhuǎn)，將部分能量轉(zhuǎn)換為渦***。渦輪帶動(dòng)壓氣機(jī)不斷吸進(jìn)空氣并進(jìn)行壓縮，使發(fā)動(dòng)機(jī)能連續(xù)工作。由壓氣機(jī)、燃燒室和驅(qū)動(dòng)壓氣機(jī)的渦輪這3個(gè)部件組成燃?xì)獍l(fā)生器，它不斷輸出具有一定可用能量的燃?xì)狻?

在發(fā)動(dòng)機(jī)控制系統(tǒng)特別是多變量數(shù)字式電子控制系統(tǒng)的研制中，要解決大量的參數(shù)擇優(yōu)問(wèn)題。由于發(fā)動(dòng)機(jī)試車非常費(fèi)時(shí)和費(fèi)錢，而且具有一定的危險(xiǎn)性，所以發(fā)動(dòng)機(jī)控制系統(tǒng)的調(diào)試全部在真實(shí)的發(fā)動(dòng)機(jī)上進(jìn)行時(shí)不現(xiàn)實(shí)的。一個(gè)準(zhǔn)確的發(fā)動(dòng)機(jī)數(shù)字孿生體可以代替發(fā)動(dòng)機(jī)進(jìn)行控制系統(tǒng)的研究和試驗(yàn)。以F100-PW-229發(fā)動(dòng)機(jī)控制系統(tǒng)設(shè)計(jì)為例，數(shù)字孿生體適用于包括調(diào)節(jié)計(jì)劃、回路增益、提前量和修正調(diào)節(jié)計(jì)劃在內(nèi)的全部控制系統(tǒng)設(shè)計(jì)。將F100-PW-229發(fā)動(dòng)機(jī)的加速過(guò)渡態(tài)的發(fā)動(dòng)機(jī)/控制系統(tǒng)數(shù)學(xué)模型的預(yù)估結(jié)果與初步放飛（IFR）發(fā)動(dòng)機(jī)的試驗(yàn)數(shù)據(jù)進(jìn)行了比較。選擇的飛行狀態(tài)是高度6960m，飛行Ma為0.8，選擇的過(guò)渡狀態(tài)是空中慢車到最大加力。在此過(guò)程中，由于轉(zhuǎn)子轉(zhuǎn)速增加和同時(shí)接通加力燃燒室，從而使推力增加，使得推力的快速響應(yīng)。比較結(jié)果表明，無(wú)論是推力還是總增壓比，數(shù)字孿生體的結(jié)果與高空試驗(yàn)臺(tái)試驗(yàn)結(jié)果都十分接近。

構(gòu)建航空發(fā)動(dòng)機(jī)的數(shù)字孿生體，可以實(shí)時(shí)反映航空發(fā)動(dòng)機(jī)的運(yùn)行狀態(tài)，預(yù)測(cè)發(fā)動(dòng)機(jī)在未來(lái)工況變動(dòng)下的性能特征，是提高發(fā)動(dòng)機(jī)可靠性和設(shè)備健康管理的先進(jìn)手段。本文闡述了航空發(fā)動(dòng)機(jī)的數(shù)字孿生體解決方案，旨在為協(xié)助國(guó)內(nèi)航空發(fā)動(dòng)機(jī)及燃?xì)廨啓C(jī)企業(yè)落地?cái)?shù)字孿生業(yè)務(wù)。

參考架構(gòu)

航空發(fā)動(dòng)機(jī)的數(shù)字孿生體架構(gòu)包括發(fā)動(dòng)機(jī)的物理系統(tǒng)、數(shù)字孿生體、測(cè)量與控制設(shè)備、以及包含了人機(jī)界面的用戶域。如圖1所示。

圖1 航空發(fā)動(dòng)機(jī)的數(shù)字孿生體架構(gòu)

航空發(fā)動(dòng)機(jī)總體性能的數(shù)字孿生體

圖2航空發(fā)動(dòng)機(jī)總體系統(tǒng)的數(shù)字孿生體示意圖。數(shù)字孿生體中包括了進(jìn)氣口、風(fēng)扇、壓氣機(jī)、燃燒室、渦輪、尾噴管等。

圖2 航空發(fā)動(dòng)機(jī)總體系統(tǒng)的數(shù)字孿生體示意圖

航空發(fā)動(dòng)機(jī)的數(shù)字孿生體可以反映的主要性能參數(shù)包括發(fā)動(dòng)機(jī)單位推力和單位燃油消耗率，實(shí)時(shí)控制參數(shù)主要包括的壓氣機(jī)增壓比、風(fēng)扇增壓比、涵道比、渦輪進(jìn)口總溫等。數(shù)字孿生體可以實(shí)時(shí)預(yù)測(cè)發(fā)動(dòng)機(jī)總體性能，包括設(shè)計(jì)點(diǎn)熱力循環(huán)分析、非設(shè)計(jì)條件下的性能分析，如速度特性、高度特性以及節(jié)流特性，（在給定油門桿位置、飛行高度、大氣條件以及調(diào)節(jié)規(guī)律的情況下，推力和耗油率等參數(shù)隨飛行馬赫數(shù)的變化關(guān)系成為發(fā)動(dòng)機(jī)速度特性；給定油門桿位置、飛行馬赫數(shù)、大氣條件以及調(diào)節(jié)規(guī)律的情況下，推力和耗油率等參數(shù)隨飛行高度的變化關(guān)系成為發(fā)動(dòng)機(jī)高度特性；給定飛行馬赫數(shù)、飛行高度、大氣條件以及調(diào)節(jié)規(guī)律的情況下，推力和耗油率等參數(shù)隨油門桿位置的變化關(guān)系成為發(fā)動(dòng)機(jī)節(jié)流特性。通常用轉(zhuǎn)子轉(zhuǎn)速、渦輪前總溫或主燃燒室供油量等參數(shù)代表油門桿位置，用發(fā)動(dòng)機(jī)推力與耗油率之間的變化關(guān)系描述發(fā)動(dòng)機(jī)的節(jié)流特性）。發(fā)動(dòng)機(jī)性能分析還包括過(guò)渡工作狀態(tài)分析，即發(fā)動(dòng)機(jī)從一個(gè)穩(wěn)定的工作狀態(tài)迅速地過(guò)渡到另一個(gè)穩(wěn)定工作狀態(tài)的過(guò)程，包括起動(dòng)過(guò)程、加速過(guò)程、減速過(guò)程、接通加力與斷開加力等。

在數(shù)字孿生體的界面中，通過(guò)傳感器測(cè)量得到發(fā)動(dòng)機(jī)的飛行高度、飛行Ma數(shù)以及燃油供給量，這次實(shí)時(shí)參數(shù)傳輸給數(shù)字孿生體后，孿生體會(huì)計(jì)算得到發(fā)動(dòng)機(jī)進(jìn)口總溫、總壓以及出口背壓，以及相應(yīng)截面的參數(shù)，包括轉(zhuǎn)子轉(zhuǎn)速，發(fā)動(dòng)機(jī)進(jìn)口流量、壓氣機(jī)增總壓比，壓氣機(jī)喘振裕度，渦輪前溫度，發(fā)動(dòng)機(jī)推力等重要參數(shù)。

圖3 航空發(fā)動(dòng)機(jī)總體性能數(shù)字孿生的參數(shù)調(diào)試界面

此外，數(shù)字孿生體還可以實(shí)現(xiàn)對(duì)單個(gè)或多個(gè)參數(shù)的虛擬驗(yàn)證。對(duì)于發(fā)動(dòng)機(jī)復(fù)雜的調(diào)節(jié)規(guī)律，例如n=nmax=const發(fā)動(dòng)機(jī)物理轉(zhuǎn)速保持最大值不變；渦輪前總溫保持最大值不變；發(fā)動(dòng)機(jī)尾噴管臨界截面可調(diào)時(shí)的組合調(diào)節(jié)規(guī)律等，在孿生體運(yùn)行過(guò)程中，可以加入script腳本文件，控制燃油供給，也可以利用控制元件來(lái)調(diào)節(jié)，仿真發(fā)動(dòng)機(jī)在設(shè)計(jì)條件或非設(shè)計(jì)條件下的穩(wěn)態(tài)仿真過(guò)程，也可以進(jìn)行過(guò)渡工作狀態(tài)的非定常過(guò)程仿真。

?

圖4 ?11km、Ma=0.85時(shí)壓氣機(jī)和渦輪運(yùn)行狀態(tài)點(diǎn)

航空發(fā)動(dòng)機(jī)主要部件的數(shù)字孿生體

除了總體性能的數(shù)字孿生體外，還可以搭建各子部件的數(shù)字孿生體，從而反映各部件的詳細(xì)性能特征。

圖5為燃燒室的數(shù)字孿生體。根據(jù)火焰筒每段的具體參數(shù)，給定相應(yīng)的換熱面積，搭建燃燒室內(nèi)的流動(dòng)和換熱網(wǎng)絡(luò)：添加燃燒模型計(jì)算燃燒室燃?xì)獾臏囟龋惶砑訉?duì)流換熱元件計(jì)算火焰筒外壁面和冷空氣的換熱過(guò)程；分別用軸向?qū)嵩蛷较驅(qū)嵩?jì)算薄壁火焰筒的導(dǎo)熱過(guò)程；用氣膜冷卻元件模擬冷卻孔形成的氣膜冷卻換熱過(guò)程；采用輻射換熱元件模擬燃?xì)夂突鹧嫱仓g的輻射換熱過(guò)程。搭建完成后，燃燒室的數(shù)字孿生體可以根據(jù)實(shí)時(shí)的燃燒室入口參數(shù)，反映出不同面積的冷卻孔的冷卻效果及溫度分布。

圖5 T56發(fā)動(dòng)機(jī)燃燒室的數(shù)字孿生體

圖6是某葉片輪盤腔室的數(shù)字孿生過(guò)程示意圖，空氣通過(guò)圖中虛線所示的引流孔進(jìn)入葉片輪盤內(nèi)部腔室結(jié)構(gòu)后沿著設(shè)計(jì)流路流出：進(jìn)口總溫為650K，總壓為13bar，引流孔總流量為20kg/s；出口靜壓為8.1bar；轉(zhuǎn)速為3000rpm；30個(gè)傾斜引流孔沿輪盤周向均布，孔的結(jié)構(gòu)參數(shù)孔徑以及傾斜引流孔角度α和β未知。建模過(guò)程中，采用旋轉(zhuǎn)通道模擬傾斜引流孔；采用轉(zhuǎn)子-靜子盤腔元件模擬輪盤腔室結(jié)構(gòu)。采用數(shù)字孿生體可以進(jìn)行敏感性分析功能，對(duì)引流孔直徑d、進(jìn)氣角α、β進(jìn)行分析等。

圖6?渦輪盤腔二次氣流系統(tǒng)的數(shù)字孿生體