胡燕華,周七二,蔡顯新
(1.南京航空航天大學(xué)能源與動力學(xué)院,江蘇南京210016;2.中航工業(yè)航空動力機械研究所,湖南株洲 412002)
渦軸發(fā)動機整機S2參數(shù)化流道造型及數(shù)值仿真
胡燕華1,2,周七二2,蔡顯新2
(1.南京航空航天大學(xué)能源與動力學(xué)院,江蘇南京210016;2.中航工業(yè)航空動力機械研究所,湖南株洲 412002)
建立了一種渦軸發(fā)動機整機S2參數(shù)化流道造型方法,采用某S2數(shù)值仿真軟件完成了某渦軸發(fā)動機的整機計算。計算表明S2參數(shù)化流道造型方法簡單、有效,仿真軟件對渦軸發(fā)動機進行整機穩(wěn)態(tài)計算的精度基本滿足計算要求,為渦軸發(fā)動機整機多學(xué)科設(shè)計優(yōu)化奠定了基礎(chǔ)。
渦軸發(fā)動機;整機;參數(shù)化流道造型;數(shù)值仿真
航空發(fā)動機數(shù)值仿真技術(shù)是將現(xiàn)代相關(guān)學(xué)科最新研究成果、多年積累的實踐經(jīng)驗及現(xiàn)代計算機技術(shù)、信息技術(shù)等融合貫通,綜合應(yīng)用于發(fā)動機研制的一種手段,是現(xiàn)代先進發(fā)動機研制技術(shù)的一個重要發(fā)展方向,是計算流體力學(xué)(CFD)、計算結(jié)構(gòu)力學(xué)(CSD)、計算機仿真技術(shù)和現(xiàn)代信息技術(shù)等最新成果在航空發(fā)動機上的綜合應(yīng)用。航空發(fā)動機的數(shù)值仿真模型根據(jù)精度不同,分為零維、一維、準(zhǔn)三維和三維。隨著計算機通信技術(shù)的迅猛發(fā)展和計算機計算能力的突飛猛進,大規(guī)模的發(fā)動機部件復(fù)雜流場、溫度場、強度分析和優(yōu)化設(shè)計成為可能,實現(xiàn)了發(fā)動機部件的三維數(shù)值仿真。但是,單獨的發(fā)動機部件數(shù)值仿真與優(yōu)化設(shè)計存在著邊界條件不真實、各部件之間匹配聯(lián)系不緊密的弱點,而這將極大的限制數(shù)值仿真技術(shù)的應(yīng)用效果。[1]因此,將發(fā)動機部件數(shù)值仿真技術(shù)的成果進一步發(fā)展,結(jié)合數(shù)值傳遞技術(shù)、多學(xué)科優(yōu)化技術(shù)與必要的實驗修正,實現(xiàn)統(tǒng)一平臺上發(fā)動機整機的數(shù)值仿真就成為技術(shù)上可行的航空發(fā)動機仿真與優(yōu)化設(shè)計的發(fā)展方向。本文應(yīng)用的某數(shù)值仿真軟件可以對部件、整機進行數(shù)值模擬,考核部件、整機的性能指標(biāo),以及各部件在整機環(huán)境中的匹配關(guān)系,將發(fā)動機總體分析能力從零維、一維提高到準(zhǔn)三維的水平。
在建立了一種渦軸發(fā)動機整機S2參數(shù)化流道造型的基礎(chǔ)上,利用S2軟件進行了某渦軸發(fā)動機總體計算工作,包括發(fā)動機整機的網(wǎng)格劃分、計算、損失模型初步給定和結(jié)果分析等。
渦軸發(fā)動機整機S2參數(shù)化流道造型由直線、圓弧、樣條曲線首尾相連組成。[2]其中直線由前后端點的X(軸向坐標(biāo))、R(徑向坐標(biāo))控制;圓弧由前后端點加中心點的X、R坐標(biāo)控制;樣條曲線由分布點的X、R坐標(biāo)控制。再在控制點中插值生成新的流道,滿足優(yōu)化要求。[3]
壓氣機部分的流道劃分和流道參數(shù)控制如圖1所示。
圖1 壓氣機流道劃分示意圖
本文利用該軟件進行了某渦軸發(fā)動機總體性能計算。發(fā)動機整機模型的計算網(wǎng)格如圖2所示,圖3至圖6分別為計算得到的馬赫數(shù)、壓力、溫度等值線圖及速度矢量圖。
圖2 計算網(wǎng)格圖
圖3 馬赫數(shù)等值線圖
圖4 壓力等值線圖
圖5 溫度等值線圖
圖6 速度矢量圖
發(fā)動機各個部件、總體的計算結(jié)果及計算結(jié)果與實驗值的比較見表1。
表1 發(fā)動機各個部件、總體的計算結(jié)果
從計算結(jié)果看,壓氣機S2計算的結(jié)果與試驗值相比較,流量偏高,效率和壓比偏低;這主要是由于該發(fā)動機軸流部分的轉(zhuǎn)子后靜子葉片帶端彎,并為串列葉柵結(jié)構(gòu),為了能夠計算對計算模型進行了處理,導(dǎo)致計算軸流部分的計算結(jié)果比試驗值偏低。
在本次計算中,燃燒室由于結(jié)構(gòu)比較復(fù)雜,在建模過程中對模型進行了簡化處理。從計算的圖形中看溫度、壓力分布參數(shù)分布比較合理。
渦輪S2計算的結(jié)果與試驗值相比較,效率和流量偏高,膨脹比偏高。
由于以上原因使得總體S2計算的功率比試驗值偏低。
通過該渦軸發(fā)動機的整機S2仿真計算驗證,各個部件的計算值與試驗值比較接近(3%以內(nèi)),能夠初步滿足設(shè)計中的應(yīng)用。程序本身自帶的損失模型用在整機的計算中還存在較大的誤差,需要通過經(jīng)驗數(shù)據(jù)修正損失模型,得到相對準(zhǔn)確的各部件的性能參數(shù)。
航空發(fā)動機整機S2數(shù)值仿真模型計算結(jié)果與試驗數(shù)據(jù)對比分析結(jié)果表明該模型的計算結(jié)果的變化規(guī)律正確,能夠初步用于航空發(fā)動機的總體性能機部件匹配性能計算,但計算精度有待于進一步提高。
發(fā)動機整機的S2數(shù)值仿真和多學(xué)科優(yōu)化是從“傳統(tǒng)設(shè)計”向“預(yù)測設(shè)計”轉(zhuǎn)換的重要途徑,對提升現(xiàn)階段航空發(fā)動機設(shè)計技術(shù)手段意義重大。[4]本文是關(guān)于發(fā)動機整機S2數(shù)值仿真的起步研究工作,為深入研究航空發(fā)動機的數(shù)值仿真技術(shù)和整機多學(xué)科設(shè)計優(yōu)化方法奠定了初步基礎(chǔ),但仍存在一些問題,如計算結(jié)果收斂性不好,需要大量的經(jīng)驗數(shù)據(jù)修正損失模型等。
針對現(xiàn)階段存在的主要問題,需要在下一個階段開展以下方面的工作:S2計算軟件的修改和完善;研究損失系數(shù)的給定,進一步提高計算精度;將S2仿真軟件與發(fā)動機整機多學(xué)科優(yōu)化平臺融合。[5]
[1]CIAM.S2 Euler code for gas flow calculation through gas turbine engine flow passage User’s guide[R].2006.
[2]王劍.S2仿真程序壓氣機部件初步計算報告[R].中國航空動力機械研究所科研報告,2006.
[3]彭劍勇.S2仿真程序燃燒室部件初步計算報告[R].中國航空動力機械研究所科研報告,2006.
[4]許如琦.S2仿真程序渦輪部件初步計算報告[R].中國航空動力機械研究所科研報告,2006.
[5]衛(wèi)剛.發(fā)動機整機仿真在工程設(shè)計中的應(yīng)用[R].中國航空學(xué)會第一屆航空發(fā)動機數(shù)值仿真與數(shù)字優(yōu)化設(shè)計學(xué)術(shù)交流會,2008.
[編校:張芙蓉]
S2 Parametric Channel Modeling and Numerical Simulation for a W hole Turbine-shaft Engine
HU Yanhua1,2ZHOU Qi'er2,CAIXianxin2
(1.Nanjing University of Aeronautics and Astronautics,Nanjing Jiangsu 210016; 2.AVIC Aero-engine Powerplant Research Institute,Zhuzhou Hunan 412002)
The paper establishes a S2 parametric channelmodelingmethod for a whole turbineshaft engine by using a certain S2 numerical simulation software to complete the calculation of a whole turbine-shaft engine.Calculations show that the S2 parametric channelmodelingmethod is simple and effective,and the accuracy after simulation software's calculating turboshaftengine'steady state can basicallymeet the calculating requirement.Thus thismay lay a good foundation formultidisciplinary design optimization of a whole turbine-shaft engine.
turboshaft engine;machine;parametric modeling flow;simulation
V235.12
A
1671-9654(2011)01-033-04
2010-12-06
胡燕華(1981-),女,湖南株洲人,工程師,工學(xué)碩士,研究方向為發(fā)動機結(jié)構(gòu)強度。