基于近似模型的變循環(huán)發(fā)動(dòng)機(jī)穩(wěn)態(tài)性能分析及優(yōu)化

韓 佳，蘇桂英，張躍學(xué)

(中國(guó)航發(fā)沈陽(yáng)發(fā)動(dòng)機(jī)研究所，沈陽(yáng)110015)

基于近似模型的變循環(huán)發(fā)動(dòng)機(jī)穩(wěn)態(tài)性能分析及優(yōu)化

韓 佳，蘇桂英，張躍學(xué)

(中國(guó)航發(fā)沈陽(yáng)發(fā)動(dòng)機(jī)研究所，沈陽(yáng)110015)

針對(duì)變循環(huán)發(fā)動(dòng)機(jī)可調(diào)節(jié)變量多的特點(diǎn)，基于變循環(huán)發(fā)動(dòng)機(jī)穩(wěn)態(tài)性能計(jì)算結(jié)果，提出了通過(guò)構(gòu)建和利用近似模型對(duì)穩(wěn)態(tài)性能進(jìn)行分析、優(yōu)化的研究方法，同時(shí)還建立了優(yōu)化求解程序。以超聲速巡航工況穩(wěn)態(tài)性能的分析及優(yōu)化為例，解析了主要變幾何部件調(diào)節(jié)對(duì)推力和耗油率的影響，獲得了滿足優(yōu)化條件的性能方案。與傳統(tǒng)大范圍變幾何部件參數(shù)研究方法對(duì)比驗(yàn)證表明，該方法在提高設(shè)計(jì)效率方面作用顯著。為變循環(huán)發(fā)動(dòng)機(jī)性能分析及控制規(guī)律設(shè)計(jì)提供了一種研究方法。

航空發(fā)動(dòng)機(jī)；變循環(huán)；變幾何部件；雙外涵；性能優(yōu)化；近似模型；響應(yīng)曲面

1 引言

變循環(huán)發(fā)動(dòng)機(jī)(VCE)能通過(guò)改變可變幾何部件/機(jī)構(gòu)的形狀、尺寸或位置來(lái)改變發(fā)動(dòng)機(jī)的熱力循環(huán)參數(shù)，從而使發(fā)動(dòng)機(jī)在各種工作條件下都具有優(yōu)良的熱力循環(huán)狀態(tài)[1-3]。但由于變循環(huán)發(fā)動(dòng)機(jī)可調(diào)部件的增加和參變量的增多，也導(dǎo)致其性能優(yōu)化較常規(guī)渦扇發(fā)動(dòng)機(jī)更為復(fù)雜[4]。因此，VCE的特性分析及性能優(yōu)化對(duì)其研究具有重要的意義。

從上世紀(jì)50年代末60年代初VCE概念提出至今，美國(guó)從未間斷對(duì)變循環(huán)技術(shù)的研究。同時(shí)，法國(guó)和日本也緊隨美國(guó)的研究步伐開(kāi)展了相關(guān)研究，積累了一定經(jīng)驗(yàn)。國(guó)內(nèi)對(duì)VCE的研究起步較晚、基礎(chǔ)較薄，對(duì)其特性的認(rèn)識(shí)也較淺。工程上，主要基于大范圍的可調(diào)部件參數(shù)研究結(jié)果，根據(jù)約束條件進(jìn)行篩選，觀察分析變化規(guī)律，細(xì)化調(diào)節(jié)范圍及步長(zhǎng)、反復(fù)計(jì)算，最終得到工程上的最優(yōu)方案?；A(chǔ)研究方面，朱之麗、王占學(xué)等開(kāi)展了變幾何部件調(diào)節(jié)對(duì)VCE性能影響的研究[5-7]，黃紅超、劉振德等開(kāi)展了基于Isight平臺(tái)的發(fā)動(dòng)機(jī)性能優(yōu)化分析[8-9]，隨陽(yáng)等提出了基于遺傳算法的VCE穩(wěn)態(tài)性能優(yōu)化方法[10]。但研究也發(fā)現(xiàn)：基于可調(diào)部件參數(shù)研究的優(yōu)化方法，存在計(jì)算量大、設(shè)計(jì)效率低且無(wú)法考慮到空間內(nèi)所有解的問(wèn)題；基于遺傳算法等的優(yōu)化方法[11-12]，雖然能快速得到VCE穩(wěn)態(tài)性能的優(yōu)化結(jié)果，但卻很難獲取變幾何部件調(diào)節(jié)對(duì)穩(wěn)態(tài)性能的影響規(guī)律，不利于對(duì)VCE特性的研究。

為此，本文提出了基于近似模型的發(fā)動(dòng)機(jī)性能分析及優(yōu)化方法。以VCE超聲速巡航狀態(tài)為例，通過(guò)構(gòu)建和求解近似模型，分析了變幾何部件調(diào)節(jié)對(duì)發(fā)動(dòng)機(jī)穩(wěn)態(tài)性能的影響規(guī)律；同時(shí)建立優(yōu)化程序，求出較優(yōu)的性能方案，并與參數(shù)研究的尋優(yōu)結(jié)果進(jìn)行對(duì)比，以期為VCE的設(shè)計(jì)和性能優(yōu)化提供幫助。

2 研究方案

以一雙外涵VCE的設(shè)計(jì)方案為基礎(chǔ)，對(duì)其超聲速巡航性能進(jìn)行分析、優(yōu)化。首先根據(jù)典型工作狀態(tài)，確定變幾何部件尋優(yōu)分析的調(diào)節(jié)范圍；應(yīng)用適當(dāng)?shù)奶摂M試驗(yàn)設(shè)計(jì)方法確定樣本方案點(diǎn)，并利用VCE穩(wěn)態(tài)性能計(jì)算程序?qū)颖痉桨更c(diǎn)進(jìn)行計(jì)算。隨后根據(jù)樣本方案點(diǎn)的計(jì)算結(jié)果，建立發(fā)動(dòng)機(jī)性能參數(shù)及約束參數(shù)關(guān)于變幾何部件調(diào)節(jié)參數(shù)的近似模型，利用近似模型對(duì)部件調(diào)節(jié)的耦合作用進(jìn)行分析。最后根據(jù)方案優(yōu)化設(shè)計(jì)中的約束條件，對(duì)近似模型直接求解，得到優(yōu)選方案，指導(dǎo)VCE設(shè)計(jì)和研究。分析及優(yōu)化流程如圖1所示。

圖1 基于近似模型的穩(wěn)態(tài)性能分析及優(yōu)化流程Fig.1 Analysisand optim ization processof VCE steady-state performance based on approximatemodels

3 近似模型的建立

超聲速巡航時(shí)，要求發(fā)動(dòng)機(jī)在不加力狀態(tài)下工作。為獲取足夠的推力，在該狀態(tài)點(diǎn)選用單外涵小涵道比模式。此時(shí)核心機(jī)驅(qū)動(dòng)風(fēng)扇(CDFS)導(dǎo)葉角度開(kāi)大，同時(shí)低壓渦輪導(dǎo)向器面積開(kāi)大、后涵道引射器(RVABI)外涵面積及噴管喉部面積縮小。本文在方案分析時(shí)，CDFS導(dǎo)葉角度開(kāi)至最大，在高/低壓物理轉(zhuǎn)速、渦輪前溫度不超限的情況下，改變低壓渦輪導(dǎo)向器面積、RVABI外涵面積及噴管喉部面積進(jìn)行性能匹配，控制規(guī)律表達(dá)式見(jiàn)式(1)。

式中：n1為低壓物理轉(zhuǎn)速，n2為高壓物理轉(zhuǎn)速，T4為渦輪前溫度，下標(biāo)max表示最大限制值。

變幾何部件參數(shù)調(diào)節(jié)范圍見(jiàn)表1，表中 ALPT代表低壓渦輪導(dǎo)向器面積，ARVABI代表后涵道引射器外涵面積，A8代表噴管喉部面積。為方便數(shù)據(jù)處理及分析，文中推力Fn、耗油率SFC均以相對(duì)變幾何部件不調(diào)節(jié)方案性能變化百分比表示。

表1 變幾何部件參數(shù)調(diào)節(jié)范圍Table 1 Parameter adjustmentdomain of variablegeometry components

在變幾何部件參數(shù)調(diào)節(jié)范圍內(nèi)，選取85個(gè)樣本方案點(diǎn)，利用VCE穩(wěn)態(tài)性能程序?qū)颖痉桨更c(diǎn)進(jìn)行計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行二次響應(yīng)曲面(RSM)[13]建模分析，建立尋優(yōu)目標(biāo)及限制參數(shù)的近似模型。

近似模型的復(fù)相關(guān)系數(shù)R2及修正的復(fù)相關(guān)系數(shù)AdjR2較高，表明模型具有高的響應(yīng)度。圖2給出了表1參數(shù)變化范圍內(nèi)，步長(zhǎng)為單位1的所有點(diǎn)組合方案的近似模型計(jì)算值相對(duì)穩(wěn)態(tài)性能程序計(jì)算值的誤差分布。從圖中可以看出：所建立的近似模型的誤差在3%以內(nèi)，且誤差在-1%～1%之間的方案占總方案的百分比高于50%。這表明該近似模型誤差較低，分布趨勢(shì)近似正態(tài)分布，可以用于參數(shù)分析及性能優(yōu)化。

4 基于近似模型的性能分析及優(yōu)化結(jié)果

4.1 變幾何部件聯(lián)合調(diào)節(jié)對(duì)性能的影響

近似模型能反映變幾何部件調(diào)節(jié)對(duì)總體性能影響的趨勢(shì)。圖3給出了任意兩個(gè)變幾何部件調(diào)節(jié)對(duì)超聲速巡航狀態(tài)推力及耗油率的影響，其中第三個(gè)參數(shù)設(shè)為取值范圍內(nèi)的中間值。

圖2 近似模型的求解誤差分布Fig.2 Solution error distribution ofapproximatemodels

圖3 任意兩個(gè)變量對(duì)推力、耗油率的影響Fig.3 Effectofany two variableson Fn and SFC

分析圖3可知，變幾何部件調(diào)節(jié)對(duì)超聲速巡航狀態(tài)推力和耗油率的影響趨勢(shì)基本一致；噴管喉部面積和低壓渦輪導(dǎo)向器面積對(duì)性能的影響較大，而后涵道引射器外涵面積對(duì)性能的影響較??；各變幾何部件調(diào)節(jié)對(duì)發(fā)動(dòng)機(jī)性能存在耦合作用。噴管喉部面積和低壓渦輪導(dǎo)向器面積匹配調(diào)節(jié)的耦合作用明顯，下面就兩者耦合作用對(duì)推力的影響進(jìn)行詳細(xì)分析。

隨著噴管喉部面積的增加，低壓渦輪膨脹比增大，渦輪做功能力增強(qiáng)，高、低壓物理轉(zhuǎn)速呈現(xiàn)增加的趨勢(shì)，而低壓物理轉(zhuǎn)速先達(dá)到限制值。隨著噴管喉部面積的進(jìn)一步增加，低壓渦輪膨脹比進(jìn)一步增大，此時(shí)必須通過(guò)減小燃燒室供油量降低渦輪前溫度，以保持低壓物理轉(zhuǎn)速不超出限制值，而此時(shí)高壓物理轉(zhuǎn)速開(kāi)始減小，推力呈下降趨勢(shì)。

隨著低壓渦輪導(dǎo)向器面積增加，高壓渦輪膨脹比增大，高壓渦輪功增加，高壓物理轉(zhuǎn)速增加，發(fā)動(dòng)機(jī)涵道比降低，推力呈增加趨勢(shì)。然而，高壓物理轉(zhuǎn)速達(dá)到限制值后，需要通過(guò)減小燃燒室供油量降低渦輪前溫度，以保持高壓物理轉(zhuǎn)速不超出限制值，此時(shí)發(fā)動(dòng)機(jī)推力呈降低趨勢(shì)。

噴管喉部面積較大時(shí)，低壓物理轉(zhuǎn)速較高，甚至達(dá)到限制值。此時(shí)，高壓物理轉(zhuǎn)速較低，高壓物理轉(zhuǎn)速可增加的裕度較大。因此，隨著低壓渦輪導(dǎo)向器面積的增加，高壓物理轉(zhuǎn)速增大，推力增加。隨著噴管喉部面積的減小，低壓物理轉(zhuǎn)速不再受限，渦輪前溫度提高，高壓物理轉(zhuǎn)速增大。此時(shí)，高壓物理轉(zhuǎn)速可增加的裕度較小，隨著低壓渦輪導(dǎo)向器面積的增加，高壓物理轉(zhuǎn)速增加直到達(dá)到限制值。因此，推力呈先增加后下降的趨勢(shì)。噴管喉部面積進(jìn)一步減小，高壓物理轉(zhuǎn)速此時(shí)已達(dá)到限制值，隨著低壓渦輪導(dǎo)向器面積的增加，渦輪前溫度降低，推力呈下降趨勢(shì)。

4.2 穩(wěn)態(tài)性能優(yōu)化結(jié)果

不加力超聲速巡航狀態(tài)對(duì)發(fā)動(dòng)機(jī)推力提出了較高的要求，同時(shí)還要求發(fā)動(dòng)機(jī)具有較低的巡航耗油率。因此，在VCE方案設(shè)計(jì)階段，選取在一定耗油率限制內(nèi)使得推力最優(yōu)的模式，對(duì)超聲速巡航狀態(tài)性能進(jìn)行優(yōu)化。

優(yōu)化目標(biāo)：Fn約束條件：

式中：x為變幾何部件調(diào)節(jié)變量。

采用傳統(tǒng)大范圍參數(shù)研究的方法，選取的參數(shù)見(jiàn)表2。需計(jì)算1 573個(gè)方案點(diǎn)。經(jīng)計(jì)算及篩選，得出超聲速巡航狀態(tài)工作點(diǎn)參數(shù)，見(jiàn)表3中傳統(tǒng)方法的優(yōu)化方案。根據(jù)所建立的近似模型及優(yōu)化限制條件，對(duì)超聲速巡航狀態(tài)方案進(jìn)行優(yōu)化，得到變幾何部件的最佳參數(shù)組合，并以VCE穩(wěn)態(tài)性能程序?qū)υ搮?shù)組合進(jìn)行驗(yàn)算，得出最終優(yōu)化方案的推力及耗油率。優(yōu)化方案1的優(yōu)化目標(biāo)是推力最大；優(yōu)化方案2的優(yōu)化目標(biāo)是推力與傳統(tǒng)方法優(yōu)化結(jié)果相等，減小變幾何部件調(diào)節(jié)量。優(yōu)化結(jié)果見(jiàn)表3中基于近似模型的優(yōu)化方案1和優(yōu)化方案2。

表2 變幾何部件參數(shù)研究范圍Table 2 Parameter study range of variable geometry components

表3 傳統(tǒng)方法和基于近似模型優(yōu)化方法的優(yōu)化結(jié)果Table 3 Optimization resultswith traditionalmethod and themethod based on approximatemodels

從表3中可以看出，在耗油率、壓縮部件穩(wěn)定裕度滿足要求的條件下，優(yōu)化方案1的無(wú)量綱推力較傳統(tǒng)優(yōu)化方法的推力高0.57個(gè)百分點(diǎn)，低壓渦輪導(dǎo)向器面積調(diào)節(jié)量降低8.39個(gè)百分點(diǎn)，RVABI外涵面積調(diào)節(jié)量降低8.24個(gè)百分點(diǎn)；推力相當(dāng)時(shí)，優(yōu)化方案2的低壓渦輪導(dǎo)向器面積調(diào)節(jié)量較傳統(tǒng)優(yōu)化方法降低8.73個(gè)百分點(diǎn)，RVABI外涵面積調(diào)節(jié)量降低10.46個(gè)百分點(diǎn)。由此可見(jiàn)，采用基于近似模型的優(yōu)化方法，可以以更少的計(jì)算量得到更優(yōu)的方案計(jì)算結(jié)果。

5 結(jié)論

以超聲速巡航狀態(tài)優(yōu)化為例，驗(yàn)證了基于近似模型的變循環(huán)發(fā)動(dòng)機(jī)優(yōu)化方法的可行性及準(zhǔn)確性。通過(guò)研究得出以下結(jié)論：

(1)超聲速巡航狀態(tài)，噴管喉部面積減小、低壓渦輪導(dǎo)向器面積增大會(huì)造成推力和耗油率有增加的趨勢(shì)，RVABI外涵面積對(duì)發(fā)動(dòng)機(jī)性能的影響較小。

(2)相比常規(guī)渦扇發(fā)動(dòng)機(jī)，變循環(huán)發(fā)動(dòng)機(jī)性能模型的非線性、強(qiáng)耦合特點(diǎn)更加明顯，變幾何部件交互作用對(duì)性能的影響較大；隨著低壓渦輪導(dǎo)向器面積的增大，推力、耗油率增加，但這種趨勢(shì)隨著噴管喉部面積的減小而減弱，在噴管喉部面積小至一定程度后，推力甚至有降低趨勢(shì)。

(3)基于近似模型的優(yōu)化算法，一方面能夠在保證發(fā)動(dòng)機(jī)穩(wěn)定裕度要求的前提下，匹配出較優(yōu)的變幾何部件調(diào)節(jié)方案，將變循環(huán)發(fā)動(dòng)機(jī)的性能發(fā)揮得更優(yōu)；另一方面能夠在滿足一定性能要求的前提下，匹配出變幾何部件調(diào)節(jié)量更小的方案，降低變幾何部件和控制系統(tǒng)的設(shè)計(jì)難度。

(4)本文的優(yōu)化流程、方法同樣適用于變循環(huán)發(fā)動(dòng)機(jī)的其他飛行狀態(tài)，有利于充分挖掘變循環(huán)發(fā)動(dòng)機(jī)的性能優(yōu)勢(shì)，提高設(shè)計(jì)研究效率，為變循環(huán)發(fā)動(dòng)機(jī)的特性研究及控制規(guī)律設(shè)計(jì)提供新思路、新方法。

[1]Vyvey P，Bosschaerts W，Villace V F，et al.Study of an airbreathing variable cycle engine[R].AIAA 2011-5758，2011.

[2]Johnson JE.Variable cycle engines-the next step in pro?pulsion evolution[R].AIAA-P-76-758，1976.

[3]方昌德.變循環(huán)發(fā)動(dòng)機(jī)及其關(guān)鍵技術(shù)[J].國(guó)際航空，2004，(7)：49—51.

[4]朱之麗，王曉波.高推重比渦扇發(fā)動(dòng)機(jī)性能尋優(yōu)分析研究[J].航空動(dòng)力學(xué)報(bào)，1999，14(7)：260—264.

[5]朱之麗，李 東.變幾何渦扇發(fā)動(dòng)機(jī)幾何調(diào)節(jié)對(duì)性能的影響[J].航空動(dòng)力學(xué)報(bào)，1999，14(1)：35—38.

[6]周 紅，王占學(xué)，劉增文，等.雙外涵變循環(huán)發(fā)動(dòng)機(jī)可變幾何特性研究[J].航空學(xué)報(bào)，2014，35(8)：2126—2135.

[7]劉增文，王占學(xué)，黃紅超，等.變循環(huán)發(fā)動(dòng)機(jī)性能數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào)，2010，25(6)：1310—1315.

[8]黃紅超，王占學(xué)，劉增文，等.基于iSIGHT的變循環(huán)發(fā)動(dòng)機(jī)性能優(yōu)化[J].機(jī)械設(shè)計(jì)與制造，2012，(2)：217—219．

[9]劉振德，薛 亮，陳玉春，等.渦扇發(fā)動(dòng)機(jī)性能優(yōu)化調(diào)整的理論與試驗(yàn)研究[J].推進(jìn)技術(shù)，2006，27(3)：211—215.

[10]隨 陽(yáng)，葉志鋒，薛益春.基于遺傳算法的變循環(huán)發(fā)動(dòng)機(jī)穩(wěn)態(tài)性能優(yōu)化[J].航天控制，2013，31(6)：17—21.

[11]劉旭東，郭迎清.基于遺傳算法的渦扇發(fā)動(dòng)機(jī)最大狀態(tài)性能尋優(yōu)[J].航空計(jì)算技術(shù)，2006，36(6)：54—58.

[12]時(shí)瑞軍，樊思齊.基于遺傳算法的渦扇發(fā)動(dòng)機(jī)多變量加速尋優(yōu)控制[J].推進(jìn)技術(shù)，2003，24(4)：357—360.

[13]薛 亮，韓萬(wàn)金.基于遺傳算法與近似模型的全局氣動(dòng)優(yōu)化方法[J].推進(jìn)技術(shù)，2008，29(3)：360—366.

Analysisand optim ization of variab le cycle engine steady-state perform ance based on app roxim atem odels

HAN Jia，SUGui-ying，ZHANG Yue-xue
(AECCShenyang Engine Research Institute，Shenyang 110015，China)

Focus on features of variable cycle enginewhich hasmany variable geometry components,based on results of steady-state engine performance calculation as well,approximatemodelwas constructed and itsutilization on performance analysisand optimization was proposed as researchmethod.Optimization solv?er program was also established.Taking performance analysis and optimization of supersonic cruise condi?tion asexample,influence of themain variable geometry components on thrustand fuel consumption wasan?alyzed and performance scheme satisfying the optimization condition wasobtained.Compared with tradition?alwide-range variable geometry components parameters studymethod,this one shows significantmerits on improving design efficiency.Itprovides a researchmethod for the performance analysis of variable cycle en?gine and control law design.

aero-engine；variable cycle；variable geometry component；double bypass；performance optimization；approximatemodel；response surface

V231.3

A

1672-2620(2017)03-0016-05

2016-11-06；

2017-02-22

韓 佳(1987-)，男，黑龍江齊齊哈爾人，工程師，碩士，主要從事航空發(fā)動(dòng)機(jī)總體性能設(shè)計(jì)工作。