渦輪葉片二維冷卻結(jié)構(gòu)參數(shù)化設(shè)計(jì)技術(shù)研究

虞跨海，楊茜，羅昌金，方坤

(1.河南科技大學(xué)規(guī)劃與建筑工程學(xué)院，河南洛陽471003；2.洛陽光電技術(shù)發(fā)展中心，河南洛陽471009；3.中航商用航空發(fā)動(dòng)機(jī)有限責(zé)任公司，上海200241)

渦輪葉片二維冷卻結(jié)構(gòu)參數(shù)化設(shè)計(jì)技術(shù)研究

虞跨海1，2，楊茜1，羅昌金3，方坤1

(1.河南科技大學(xué)規(guī)劃與建筑工程學(xué)院，河南洛陽471003；2.洛陽光電技術(shù)發(fā)展中心，河南洛陽471009；3.中航商用航空發(fā)動(dòng)機(jī)有限責(zé)任公司，上海200241)

研究了渦輪葉片二維冷卻結(jié)構(gòu)的參數(shù)化設(shè)計(jì)技術(shù)。采用參數(shù)控制點(diǎn)方法實(shí)現(xiàn)冷卻葉片壁面的變厚度設(shè)計(jì)，采用隔肋數(shù)量、位置參數(shù)、偏轉(zhuǎn)角度實(shí)現(xiàn)任意數(shù)量和形式的冷卻腔造型，根據(jù)前緣縮進(jìn)參數(shù)確定冷卻通道前緣切線弧位置，通過尾緣切割參數(shù)實(shí)現(xiàn)半劈縫和全劈縫尾緣結(jié)構(gòu)設(shè)計(jì)。結(jié)合葉片外形造型技術(shù)開發(fā)了造型設(shè)計(jì)程序，該程序可建立包含任意形式冷卻通道和常用尾緣結(jié)構(gòu)的變壁厚二維冷卻葉片模型。

渦輪冷卻葉片；參數(shù)化設(shè)計(jì)；隔肋；尾緣劈縫

1 引言

渦輪前燃?xì)馊肟跍囟入S著燃?xì)廨啓C(jī)性能要求的提高而不斷提高，已遠(yuǎn)遠(yuǎn)超過葉片材料的耐受溫度，給葉片設(shè)計(jì)帶來巨大挑戰(zhàn)。目前，冷卻葉片技術(shù)作為主要解決手段，已逐漸應(yīng)用于先進(jìn)燃?xì)廨啓C(jī)，冷卻葉片設(shè)計(jì)也越來越復(fù)雜，先進(jìn)冷卻技術(shù)已成為燃?xì)廨啓C(jī)渦輪研制的關(guān)鍵[1，2]。

隨著各學(xué)科分析、數(shù)值仿真技術(shù)和優(yōu)化理論的發(fā)展，特別是多學(xué)科設(shè)計(jì)優(yōu)化理論的提出和興起，基于數(shù)值仿真模型的多學(xué)科設(shè)計(jì)優(yōu)化成為渦輪葉片設(shè)計(jì)的重要手段[3～5]。筆者也開展了直通道冷卻葉片、葉片外形的多學(xué)科設(shè)計(jì)優(yōu)化研究[6，7]。由于多學(xué)科設(shè)計(jì)優(yōu)化必須以參數(shù)化造型為基礎(chǔ)，且鑒于渦輪冷卻葉片設(shè)計(jì)的實(shí)際需求，葉片三維參數(shù)化建模技術(shù)的研究逐漸興起[6，8～11]。但渦輪冷卻葉片結(jié)構(gòu)復(fù)雜、形式多樣，很難直接建立其三維參數(shù)化模型。而基于三維模型仿真的渦輪冷卻葉片多學(xué)科設(shè)計(jì)優(yōu)化，往往由于計(jì)算成本和優(yōu)化系統(tǒng)運(yùn)行可靠性等問題導(dǎo)致優(yōu)化失敗。我國科研院所和高校，積累了大量關(guān)于二維葉片氣動(dòng)、傳熱、強(qiáng)度等學(xué)科的設(shè)計(jì)理論和成熟算法程序[12]，若能充分利用其技術(shù)和經(jīng)驗(yàn)，先開展二維冷卻結(jié)構(gòu)的設(shè)計(jì)優(yōu)化，再以此為基礎(chǔ)進(jìn)行三維設(shè)計(jì)，可極大地縮短冷卻葉片的設(shè)計(jì)周期，降低設(shè)計(jì)成本。

目前，國內(nèi)開展的葉片冷卻結(jié)構(gòu)參數(shù)化設(shè)計(jì)技術(shù)研究，基本上是基于成熟CAD軟件二次開發(fā)來實(shí)現(xiàn)[9～11]，其三維造型結(jié)果很難直接用于二維冷卻結(jié)構(gòu)設(shè)計(jì)，造成三維設(shè)計(jì)與二維設(shè)計(jì)脫節(jié)。為此，本文開展了渦輪葉片二維復(fù)雜冷卻結(jié)構(gòu)的參數(shù)化設(shè)計(jì)技術(shù)研究，并基于數(shù)值方法和自編程序，實(shí)現(xiàn)了包含任意形式冷卻通道和常用三種尾緣結(jié)構(gòu)的變壁厚二維冷卻葉片造型設(shè)計(jì)，可為開展冷卻葉片二維結(jié)構(gòu)設(shè)計(jì)優(yōu)化提供自動(dòng)造型支撐。

2 葉片冷卻通道型線設(shè)計(jì)

變壁厚渦輪冷卻葉片可根據(jù)壁面的應(yīng)力、溫度分布情況和冷卻需求，合理設(shè)計(jì)壁面厚度沿弦向的分布，有利于葉片既注重氣動(dòng)、傳熱性能，又能兼顧強(qiáng)度、結(jié)構(gòu)重量的合理設(shè)計(jì)。本文基于函數(shù)解析方法實(shí)現(xiàn)變壁厚渦輪冷卻葉片的參數(shù)化設(shè)計(jì)，通過在葉片外型線上定義若干參數(shù)點(diǎn)，用三次樣條插值方法構(gòu)造葉片壁面厚度函數(shù)，建立變壁厚渦輪冷卻葉片。

根據(jù)葉型主要參數(shù)建立方程矩陣，可得到葉身壓力面和吸力面型線[10]。圖1所示為渦輪冷卻葉片截面參數(shù)化設(shè)計(jì)示意圖，隔肋將葉片冷卻腔分割成若干回流腔，定義若干壁厚參數(shù)點(diǎn)(a0，b0，c0，d0，e0，f0，a1，b1，c1，d1，e1，f1)，其中a0，a1，f0，f1分別為外型線壓力面、吸力面型線與前緣和尾緣的切點(diǎn)。根據(jù)壁厚參數(shù)構(gòu)造葉片壁厚函數(shù)，再結(jié)合葉片外型線函數(shù)，可計(jì)算得到葉片冷卻通道壁面型線[13]。

圖1 冷卻通道參數(shù)化壁面控制點(diǎn)Fig.1 Parametric controlling points of cooling passage

由于參數(shù)三次樣條具有二階連續(xù)性，曲率變化小，構(gòu)造較為簡單，因此本文采用三次樣條對(duì)造壁面厚度進(jìn)行插值構(gòu)造函數(shù)。其構(gòu)造形式為：

3 隔肋參數(shù)化設(shè)計(jì)

得到冷卻通道外形型線后，以冷卻通道壓力面型線前緣、尾緣端點(diǎn)為起點(diǎn)和終點(diǎn)，沿x軸方向投影，并歸一化處理，即隔肋位置坐標(biāo)區(qū)間為(0，1)。本文研究了三類隔肋的參數(shù)化設(shè)計(jì)：①隔肋中心線垂直于葉片壓力面外型線，②隔肋中心線垂直于冷卻通道壓力面型線，③隔肋中心線與葉片壓力面型線法線成任意角度。其中①可認(rèn)為是③的特例，因此本文主要對(duì)后兩種造型方法展開討論。

圖2、圖3分別為垂直于冷卻通道壓力面型線和任意角度隔肋造型設(shè)計(jì)示意圖。根據(jù)輸入的各隔肋寬度(b1，b2，……，bn)、隔肋數(shù)量、偏轉(zhuǎn)角度φbi和中心位置坐標(biāo)(x1，x2，……，xn)，從冷卻通道壓力面出發(fā)，結(jié)合冷卻通道壓力面型線函數(shù)，建立垂直于壓力面的各隔肋中心線方程?(x)；以該中心線為對(duì)稱軸，根據(jù)隔肋寬度，建立隔肋的左右兩條端線方程?1(x)和?2(x)；采用數(shù)值優(yōu)化算法搜索得到端線與冷卻通道壓力面和吸力面的交點(diǎn)坐標(biāo)，最終可輸出隔肋造型離散數(shù)據(jù)。任意角度隔肋時(shí)，偏轉(zhuǎn)角度方向以通道型線內(nèi)法向順時(shí)針轉(zhuǎn)動(dòng)為正，反之為負(fù)。

圖2 垂直于冷卻通道壓力面型線隔肋造型設(shè)計(jì)示意圖Fig.2 Parametric design of the ribs vertical to pressure side of cooling passage profile

圖3 任意角度隔肋造型設(shè)計(jì)示意圖Fig.3 Parametric design of the ribs with variable angles

隔肋厚度、偏轉(zhuǎn)角度等數(shù)據(jù)根據(jù)葉片具體形狀確定。造型過程中需關(guān)注各隔肋的相對(duì)位置以保證各冷卻腔的幾何尺寸，避免相鄰隔肋相交。特別應(yīng)關(guān)注第一個(gè)和最后一個(gè)隔肋的位置，確保前緣第一個(gè)腔和尾緣最后一個(gè)腔的空間。

4 通道前緣與尾緣參數(shù)化設(shè)計(jì)

4.1 前緣切線弧連接設(shè)計(jì)

冷卻通道前緣結(jié)構(gòu)由圓弧構(gòu)造。為保證幾何形狀的光滑性，在連接點(diǎn)位置要求具有一階連續(xù)性，因此文中冷卻通道前緣由切線弧構(gòu)造。圖4所示為冷卻通道前緣切線弧結(jié)構(gòu)示意圖。根據(jù)前緣縮進(jìn)參數(shù)確定冷卻通道壓力面一點(diǎn)為圓弧連接點(diǎn)A，作壓力面在A點(diǎn)位置的法線AB′。由于前緣圓弧與冷卻通道型線在A點(diǎn)相切，故前緣圓弧圓心必然在過A點(diǎn)的法線上。建立AB′直線方程φ(x)，設(shè)冷卻通道吸力面型線方程為f(x)，則圓心坐標(biāo)和切線弧半徑的計(jì)算可描述為以下優(yōu)化問題：在f(x)上尋找一點(diǎn)B，使得從B點(diǎn)引出的法線與AB′相交，若交點(diǎn)為O，使得OB=OA，則該點(diǎn)為冷卻通道前緣切線弧的圓心，OA為切線弧半徑。

圖4 冷卻通道前緣切線弧結(jié)構(gòu)示意圖Fig.4 Tangent arc parametric design at cooling passage leading edge

切線弧連接冷卻通道前緣，能保證前緣圓弧與冷卻通道壓力面、吸力面型線的一階連續(xù)性，使得整段冷卻通道型線光滑，具有較好的氣動(dòng)和加工性能。但切線弧前緣構(gòu)造方式，可能會(huì)由于A點(diǎn)位置選取過于靠近葉片外形前緣，而使得求解失效或葉片前緣位置厚度過小。因此，在造型過程中，需關(guān)注對(duì)葉片前緣縮進(jìn)參數(shù)的控制。

4.2 尾緣劈縫結(jié)構(gòu)設(shè)計(jì)

本文開展了冷卻葉片尾緣圓弧、半劈縫和全劈縫三種常用結(jié)構(gòu)的參數(shù)化造型技術(shù)研究。只需定義尾緣切割參數(shù)和尾緣形式，即可獲得所需幾何造型。以下主要討論半劈縫和全劈縫的參數(shù)化設(shè)計(jì)技術(shù)。

(1)全劈縫結(jié)構(gòu)設(shè)計(jì)

圖5為尾緣全劈縫結(jié)構(gòu)造型示意圖。根據(jù)尾緣切割參數(shù)δtx和葉片壓力面型線方程yp=fp(x)，可得尾緣壓力面切割點(diǎn)A的坐標(biāo)。以過A點(diǎn)的壓力外型線法線為尾緣切割線，根據(jù)葉片外型線吸力面方程、葉片冷卻通道型線吸力面方程和冷卻通道壓力面型線方程，采用優(yōu)化算法分別搜索相關(guān)交點(diǎn)B、D、C的坐標(biāo)，去除尾緣其余部分，建立尾緣全劈縫結(jié)構(gòu)。圖5虛線顯示部分即為造型去除幾何部分，只需δtx，即可實(shí)現(xiàn)尾緣全劈縫造型。

圖5 尾緣全劈縫結(jié)構(gòu)造型Fig.5 Trailing edge with whole slot

(2)半劈縫結(jié)構(gòu)設(shè)計(jì)

圖6為尾緣半劈縫結(jié)構(gòu)造型示意圖。根據(jù)δtx和yp=fp(x)，計(jì)算得到尾緣壓力面切割點(diǎn)A。以過A點(diǎn)的壓力面外型線法線方向?yàn)槲簿壐罹€方向，根據(jù)葉片冷卻通道壓力面型線方程，采用優(yōu)化算法搜索冷卻通道壓力面型線與割線交點(diǎn)D的坐標(biāo)，對(duì)吸力面壁厚部分做圓弧處理，建立尾緣半劈縫結(jié)構(gòu)。

圖6 尾緣半劈縫結(jié)構(gòu)造型Fig.6 Trailing edge with half slot

另外，尾緣割線方向也可根據(jù)過A點(diǎn)法線方向偏轉(zhuǎn)一定角度建立；尾緣劈縫結(jié)構(gòu)造型過程中，要避免尾緣切割參數(shù)過小，導(dǎo)致割線不能與相關(guān)型線相交的問題；由于加工工藝的因素，一般要求尾緣氣流出口寬度不小于0.5 mm，若出口寬度過小，可通過調(diào)整尾緣附近的壁厚或增大尾緣切割參數(shù)來增大出口厚度。

5 軟件開發(fā)

本文在渦輪冷卻葉片外型線、冷卻通道型線、隔肋造型、前緣造型和尾緣造型參數(shù)化設(shè)計(jì)的基礎(chǔ)上，編制了渦輪葉片二維冷卻結(jié)構(gòu)造型程序，其界面如圖7所示，只需輸入相關(guān)設(shè)計(jì)參數(shù)即可根據(jù)用戶需求輸出二維幾何模型。圖8為渦輪冷卻葉片三種不同造型結(jié)果。其中，圖8(a)為尾緣封閉七腔冷卻葉片造型，隔肋中心線垂直于葉片壓力面外型線；圖8 (b)為尾緣半劈縫八腔冷卻葉片，隔肋中心線垂直于冷卻通道壓力面型線；圖8(c)為尾緣全劈縫七腔冷卻葉片，隔肋中心線為任意角度造型。此外，為配合Stan5、Flownet等程序，軟件可計(jì)算輸出葉片面積、各冷卻腔面積、各冷卻腔壁面厚度均值等參數(shù)。

圖7 渦輪葉片二維冷卻結(jié)構(gòu)造型設(shè)計(jì)軟件界面Fig.7 Software interface of turbine blade 2D cooling structure

圖8 渦輪冷卻葉片造型結(jié)果Fig.8 Design examples of cooling turbine blade

6 結(jié)論

本文對(duì)渦輪葉片二維冷卻結(jié)構(gòu)參數(shù)化設(shè)計(jì)方法開展了研究，實(shí)現(xiàn)了任意數(shù)量和形式冷卻腔變壁厚二維渦輪冷卻葉片造型設(shè)計(jì)，并且包含常用的三種尾緣造型形式，可為開展渦輪冷卻葉片二維結(jié)構(gòu)設(shè)計(jì)優(yōu)化提供自動(dòng)造型支撐。本文的主要工作為：

(1)提出了冷卻葉片的變壁厚設(shè)計(jì)方法，可根據(jù)葉片需求合理設(shè)計(jì)葉片壁面尺寸；

(2)實(shí)現(xiàn)了二維冷卻葉片任意數(shù)量和角度隔肋的參數(shù)化設(shè)計(jì)；

(3)實(shí)現(xiàn)了圓弧、半劈縫和全劈縫三種形式的尾緣冷卻結(jié)構(gòu)參數(shù)化設(shè)計(jì)；

(4)最終形成了渦輪葉片二維冷卻結(jié)構(gòu)的參數(shù)化設(shè)計(jì)方法，并編制了造型程序，開發(fā)了造型設(shè)計(jì)軟件。

[1]Iacovides H，Raisee M.Recent Progress in the Computa?tion of Flow and Heat Transfer in Internal Cooling Passag?es of Turbine Blades[J].International Journal of Heat and Fluid Flow，1999，20：320—328.

[2]Koff B L.Gas Turbine Technology Evolution-A Design?er's Perspective[R].AIAA 2003-2722，2003.

[3]Talya S S，Rajadas J N.Multidisciplinary Analysis and De?sign Optimization Procedure for Cooled Gas Turbine Blades[R].AIAA 2000-4877，2000.

[4]Hasenjager M，Sendhoff B，Sonoda T，et al.Three Dimen?sional Aerodynamic Optimization for an Ultra-Low Aspect RatioTransonicTurbineStatorBlade[R].ASME GT2005-68680，2005.

[5]楊俊杰，王榮橋，樊江，等.渦輪葉片的氣動(dòng)-熱-結(jié)構(gòu)多學(xué)科設(shè)計(jì)優(yōu)化研究[J].航空動(dòng)力學(xué)報(bào)，2010，25(3)：617—622.

[6]虞跨海，岳珠峰.渦輪冷卻葉片參數(shù)化建模及多學(xué)科設(shè)計(jì)優(yōu)化[J].航空動(dòng)力學(xué)報(bào)，2007，22(8)：1346—1351.

[7]虞跨海，王金生，楊茜，等.基于近似的渦輪冷卻葉片外形多學(xué)科設(shè)計(jì)優(yōu)化[J].機(jī)械工程學(xué)報(bào)，2011，47(10)：106—112.

[8]KIM Yeoung，KIM Li-ra，JUN Cha-soo.Parametric De?sign of a Part with Free-Form Surfaces[J].Journal of Zheji--ang University SCIENCEA，2006，7(9)：1530—1534.

[9]宋玉旺，席平.基于特征造型技術(shù)的渦輪葉片參數(shù)化設(shè)計(jì)[J].北京航空航天大學(xué)學(xué)報(bào)，2004，30(4)：321—324.

[10]虞跨海，李立州，岳珠峰.基于解析成型及特征造型技術(shù)的渦輪冷卻葉片參數(shù)化設(shè)計(jì)[J].推進(jìn)技術(shù)，2007，28 (6)：637—640.

[11]曹率，莫蓉，宮中偉.航空渦輪葉片氣膜孔參數(shù)化建模研究[J].航空計(jì)算技術(shù)，2011，41(3)：77—81.

[12]卿雄杰，潘炳華，曾軍.渦輪葉片外換熱計(jì)算程序?qū)Ρ妊芯縖J].航空動(dòng)力學(xué)報(bào)，2010，25(4)：741—747.

[13]虞跨海，楊茜，岳珠峰.變壁厚渦輪冷卻葉片參數(shù)化造型設(shè)計(jì)[J].機(jī)械設(shè)計(jì)，2012，29(2)：5—7.

Parametric Design Method of 2D Turbine Blade Cooling Structure

YU Kua-hai1，2，YANG Xi1，LUO Chang-jin3，F(xiàn)ANG Kun1
(1.School of Architecture Engineering，Henan University of Science and Technology，Luoyang 471003；2.Luoyang Opt-Electro Development Center，Luoyang 471009；3.AVIC Commercial Aircraft Engine CO.LTD，Shanghai 200241)

The parametric design method for 2D turbine blade cooling structure is studied.Parametric con?trolling point method is used to realize variable thickness design of blade wall.With the parameters of rib number,rib position and deflection angles,a variable number cooling chamber structure can be built.The retraction parameter decides the position of tangent circular arc at cooling passage leading edge.A division parameter is defined to get the trailing edge structure with whole slot or half slot.Then,combined with the blade profiles design technique,a 2D cooling turbine blade parametric design procedure is developed.With the procedure,a 2D turbine blade with variable thickness,including cooling passage in any forms and com?mon-used trailing edge can be easily built.

turbine cooling blade；parametric design；rib；trailing edge slot

V231.1

A

1672-2620(2013)01-0012-04

2012-06-11；

2012-12-12

國家自然科學(xué)基金(51105132)；河南省基礎(chǔ)與前沿技術(shù)研究計(jì)劃項(xiàng)目(112300410166)；河南省教育廳自然科學(xué)基金(2011B590001)

虞跨海(1982-)，男，浙江義烏人，博士，副教授，研究方向?yàn)楹娇胀七M(jìn)系統(tǒng)、多學(xué)科設(shè)計(jì)優(yōu)化。