趙大勇

（中航工業(yè)沈陽(yáng)發(fā)動(dòng)機(jī)設(shè)計(jì)研究所，沈陽(yáng) 110015）

渦輪彎扭葉片氣動(dòng)設(shè)計(jì)工程方法研究

趙大勇

（中航工業(yè)沈陽(yáng)發(fā)動(dòng)機(jī)設(shè)計(jì)研究所，沈陽(yáng) 110015）

采用彎扭葉片設(shè)計(jì)是提高燃?xì)鉁u輪效率的主要手段之一，而構(gòu)建彎扭葉片積疊線的參數(shù)化模式是設(shè)計(jì)渦輪彎扭葉片的重要基礎(chǔ)。為了滿足各種渦輪葉片彎扭設(shè)計(jì)的工程需要，提出了采用圓弧型、拋物線型、分段樣條曲線積疊線規(guī)律的渦輪彎扭葉片設(shè)計(jì)方法，編制了積疊線生成代碼，并推薦了各方法的應(yīng)用條件和方式；提出的回歸正交設(shè)計(jì)方法可以對(duì)分段樣條曲線積疊線的渦輪彎扭葉片進(jìn)行尋優(yōu)設(shè)計(jì)。所提出的彎扭葉片工程設(shè)計(jì)方法已經(jīng)在多個(gè)渦輪氣動(dòng)設(shè)計(jì)中得到應(yīng)用。

彎扭葉片；優(yōu)化；渦輪；氣動(dòng)設(shè)計(jì)；積疊線；航空發(fā)動(dòng)機(jī)

0 引言

隨著航空發(fā)動(dòng)機(jī)性能不斷提高，對(duì)渦輪葉片氣動(dòng)設(shè)計(jì)要求也越來(lái)越高。為了提高渦輪效率，滿足總體性能指標(biāo)，渦輪葉片設(shè)計(jì)不斷采用新技術(shù)，包括需要冷卻的低壓渦輪導(dǎo)葉和動(dòng)葉，大展弦比的渦輪葉片，具有大擴(kuò)張角端壁的渦輪導(dǎo)葉等。設(shè)計(jì)方法經(jīng)歷了初步（1維）設(shè)計(jì)方法、軸對(duì)稱通流（2維）設(shè)計(jì)

方法、求解3維非穩(wěn)態(tài)N-S方程的CFD方法[1]。在渦輪3維氣動(dòng)設(shè)計(jì)中，根據(jù)研究對(duì)象不同可細(xì)分為單排葉片、單級(jí)渦輪和多級(jí)渦輪的3維計(jì)算分析；根據(jù)冷卻不同可分為無(wú)冷氣渦輪流場(chǎng)、加源項(xiàng)方式的有冷氣流場(chǎng)和冷卻孔分網(wǎng)格考慮冷氣的流場(chǎng)分析；進(jìn)一步引入時(shí)間項(xiàng)，可分為穩(wěn)態(tài)流場(chǎng)和非穩(wěn)態(tài)流場(chǎng)分析。

渦輪葉片彎扭設(shè)計(jì)使葉片3維成型設(shè)計(jì)邁上了1個(gè)新臺(tái)階，采用彎扭葉片是提高渦輪效率的主要手段之一。近20年來(lái),彎扭葉片設(shè)計(jì)已經(jīng)被廣泛研究[2-3]，在國(guó)內(nèi)外多型航空發(fā)動(dòng)機(jī)上得到應(yīng)用，如湍達(dá)系列發(fā)動(dòng)機(jī)等[4]。構(gòu)建彎扭葉片積疊線是渦輪彎扭葉片設(shè)計(jì)的重要基礎(chǔ)。目前國(guó)內(nèi)彎扭葉片設(shè)計(jì)中，設(shè)計(jì)師一般根據(jù)經(jīng)驗(yàn)設(shè)計(jì)積疊線形式。由于積疊線空間解范圍很大，所以設(shè)計(jì)結(jié)果隨機(jī)性很大，所得設(shè)計(jì)解是否為最優(yōu)解或較優(yōu)解很難評(píng)判；而且需要花費(fèi)大量時(shí)間。渦輪彎扭葉片研究公開發(fā)表的論文多采用標(biāo)準(zhǔn)葉型（如NACA65葉型）進(jìn)行機(jī)理性研究，沒有構(gòu)建葉片積疊線的參數(shù)化程序，設(shè)計(jì)結(jié)果受設(shè)計(jì)工程師經(jīng)驗(yàn)和人工因素的影響較大，而且積疊線的構(gòu)建形式也相對(duì)單一，限定了葉片的可能形式，對(duì)其后葉片優(yōu)化設(shè)計(jì)工作有較大影響。在彎扭葉片設(shè)計(jì)中對(duì)采用的葉片積疊線設(shè)計(jì)方法進(jìn)行總結(jié)，編制程序代碼，定義積疊方式，并界定各種方法的適用情況和方法，對(duì)消除設(shè)計(jì)師人工隨機(jī)性影響尤為重要，而且可以提高效率，縮短設(shè)計(jì)周期，對(duì)設(shè)計(jì)解的覆蓋面有保障，大幅提高獲得較優(yōu)解的幾率，達(dá)到提高渦輪工作效率的目的。

1 渦輪彎扭葉片積疊線參數(shù)化設(shè)計(jì)方法

某一渦輪葉柵出口能量損失典型分布如圖1所示。根據(jù)低壓渦輪導(dǎo)葉出口能量損失分布規(guī)律，葉柵損失區(qū)域可以分為3部分：中部低損失區(qū)和根、尖逐漸增大的高損失區(qū)。葉柵中部主要是葉型的摩擦損失，靠近端壁為摩擦損失和2次流損失的疊加。

葉片彎扭設(shè)計(jì)的作用主要是降低端壁附近區(qū)域的2次流損失區(qū)[2]。為了降低端壁附近區(qū)域2次流損失，同時(shí)考慮結(jié)構(gòu)裝配的工程問(wèn)題，提出了分段樣條曲線、圓弧型、拋物線型積疊線規(guī)律的渦輪彎扭葉片設(shè)計(jì)方法。

1.1 分段樣條曲線積疊線設(shè)計(jì)方法

分段樣條曲線積疊線如圖2所示。從圖中可見，積疊線共分為5段，包括3個(gè)直線段和2個(gè)3次樣條曲線段。圖中，D、XH、XT、BH、BM、BT為用葉片高度進(jìn)行無(wú)量綱處理后的值：D為相對(duì)葉根基準(zhǔn)葉尖沿切線方向的偏移量，或中部為徑向的開關(guān)；XH、XT分別為根部和尖部直線段所占高度；BH、BM、BT分別為各部分直線段所占高度；AH、AT分別為葉根和葉尖的傾斜角度。

有時(shí)希望生成的積疊線的中部直線段方向保持為徑向，為了節(jié)省自變量個(gè)數(shù)，用D來(lái)實(shí)現(xiàn)。當(dāng)D取值為1.5時(shí)（葉尖偏移量為1.5倍的葉高，選擇1個(gè)不具備工程設(shè)計(jì)意義的數(shù)值），D值只代表計(jì)算開關(guān)，生成的積疊線將保持中部直線段為徑向。

采用該方法可以生成如圖3所示的各種積疊線，包括不同程度的正彎、反彎，以及J、S型的彎曲規(guī)律。傾斜直線、3次樣條等方式也可以作為本方法的1個(gè)特例。

1.2 圓弧型積疊線設(shè)計(jì)方法

葉片彎扭設(shè)計(jì)增加了冷卻葉片結(jié)構(gòu)的設(shè)計(jì)難度，因?yàn)槠浔仨毧紤]結(jié)構(gòu)通氣導(dǎo)管裝配問(wèn)題，而影響導(dǎo)管安裝的不僅是葉片根部、尖部?jī)A角的具體數(shù)值，還包括葉片沿徑向的彎曲規(guī)律分布。為此，本文提出了圓弧型積疊線設(shè)計(jì)思想。

渦輪葉片圓弧型積疊線可以通過(guò)求解2次曲線方程得到。求解2次曲線方程需要3個(gè)約束條件，首先是定義葉片根部為坐標(biāo)圓點(diǎn)，再附加2個(gè)約束條件，沿用設(shè)計(jì)人員的使用習(xí)慣，選擇根部和尖部?jī)A角。進(jìn)行導(dǎo)葉積疊設(shè)計(jì)時(shí)按照最大厚度積疊，再附加圓弧型積疊規(guī)律，可以降低冷卻葉片結(jié)構(gòu)設(shè)計(jì)難度。采用圓弧型規(guī)律生成的1族積疊線如圖4所示。該圓弧型積疊線設(shè)計(jì)，成功解決了在某型發(fā)動(dòng)機(jī)低壓渦輪冷卻葉片設(shè)計(jì)中遇到的通氣導(dǎo)管裝配問(wèn)題，還成功應(yīng)用于某小型燃?xì)廨啓C(jī)高壓渦輪導(dǎo)葉設(shè)計(jì)中。

1.3 改進(jìn)的拋物線型積疊線設(shè)計(jì)方法

拋物線型積疊線控制參數(shù)少，易于實(shí)現(xiàn)優(yōu)化，但不能用于設(shè)計(jì)S型葉片，對(duì)葉片內(nèi)導(dǎo)管安裝也不是很有利。采用拋物線型規(guī)律生成的1族積疊線如圖5所示。為了彌補(bǔ)該缺陷，用分段直線和樣條方法相結(jié)合的積疊方式設(shè)計(jì)與拋物線型類似的葉片。采用該方法的優(yōu)點(diǎn)是自變量數(shù)量少，但可能會(huì)丟失最優(yōu)解，可在實(shí)心葉片初步設(shè)計(jì)時(shí)使用，然后再用分段方法尋優(yōu)。

2 葉片積疊線規(guī)律尋優(yōu)方法

目前沒有通用的經(jīng)驗(yàn)公式和方法可以直接得到與葉型相匹配的積疊規(guī)律。對(duì)于工程設(shè)計(jì)任務(wù)，時(shí)間節(jié)點(diǎn)的限制是重要的約束條件，如何在有限時(shí)間內(nèi)快速得到最優(yōu)積疊規(guī)律是彎扭葉片設(shè)計(jì)的重點(diǎn)和難點(diǎn)。

2.1 圓弧型和拋物線型積疊線尋優(yōu)方法

采用圓弧型和拋物線方法設(shè)計(jì)彎扭葉片，因?yàn)樽宰兞恐挥?個(gè)，可以采用黃金分割法[5]或精化網(wǎng)格法實(shí)施優(yōu)化。經(jīng)過(guò)多個(gè)算例的計(jì)算，相鄰算例的目標(biāo)函數(shù)差小于程序的計(jì)算精度時(shí)，尋優(yōu)過(guò)程結(jié)束，可以得到1個(gè)等值線圖。優(yōu)化過(guò)程示例見表1。表中各行、列分別表示積疊線根部、尖部?jī)A斜角度。依據(jù)表1數(shù)據(jù)繪制的效率等值線如圖6所示。圖中橫、縱坐標(biāo)分別表示積疊線根部、尖部?jī)A斜角度（αH和αT）。各計(jì)算點(diǎn)積疊線分布規(guī)律如圖7所示，圖中黑線對(duì)應(yīng)表1和圖6中的效率取得最高值的積疊規(guī)律。圓弧型和拋物線型積疊線的葉片3維示意圖如圖8所示。

2.2 分段樣條曲線積疊線尋優(yōu)方法

分段樣條曲線方法的自變量數(shù)目較多，包含8個(gè)獨(dú)立變化的自變量，計(jì)算域空間較大，在有限時(shí)間內(nèi)快速尋優(yōu)是1個(gè)難點(diǎn)，可采用回歸正交設(shè)計(jì)方法尋優(yōu)[6-7]。

正交表是根據(jù)組合理論按照一定規(guī)律構(gòu)造的表格，在試驗(yàn)設(shè)計(jì)中應(yīng)用廣泛。以正交表為工具安排試驗(yàn)方案和進(jìn)行結(jié)果分析的試驗(yàn)稱為正交試驗(yàn)，適用于多因素、多指標(biāo)、多因素間存在交互作用、具有隨機(jī)誤差的試驗(yàn)。通過(guò)正交試驗(yàn)，可以分析各因素及其交互作用對(duì)試驗(yàn)指標(biāo)的影響，按其重要程度找出主次關(guān)系，并確定對(duì)試驗(yàn)指標(biāo)的最優(yōu)工藝條件[7]。

在彎扭葉片設(shè)計(jì)中，積疊規(guī)律受葉片葉型、子午流道形狀和工作條件等多方面的影響，目前尚不能用解析方法求解。葉片積疊線按照上述參數(shù)化表示以后，其解的空間仍然較大，且各參數(shù)之間存在較強(qiáng)的相互影響?？紤]到受加工水平和流場(chǎng)計(jì)算水平的影響，設(shè)定葉片優(yōu)化過(guò)程結(jié)束的衡量標(biāo)準(zhǔn)為相近計(jì)算結(jié)果的誤差小于精度目標(biāo)。上述特征決定了宜采用正交設(shè)計(jì)方法。根據(jù)葉片積疊線參數(shù)化的方法，選擇L18（2×37）正交表（見表2）設(shè)計(jì)比較適宜。

表1 精化網(wǎng)格法（拋物線積疊線）初步優(yōu)化示例

表2中8列數(shù)字表示積疊線參數(shù)化后的自變量，分別表示 D、AH、AT、XH、XT、BH、BM、BT參數(shù)；18行數(shù)字代表1輪試驗(yàn)（計(jì)算）次數(shù)。根據(jù)數(shù)理統(tǒng)計(jì)方法確定數(shù)值試驗(yàn)指標(biāo)之分項(xiàng)和、指標(biāo)均值和極差，從而評(píng)定各項(xiàng)參數(shù)對(duì)級(jí)效率影響的重要程度，依照各項(xiàng)參數(shù)指標(biāo)值的極差大小排出重要性順序，推論出最優(yōu)參數(shù)組合方案作為下一輪尋優(yōu)的基準(zhǔn)。

表2 L18（2×37）正交表

在前述拋物線優(yōu)化的基礎(chǔ)上，合理選擇分段樣條曲線積疊線參數(shù)，見表3。按照表3取值，得到18個(gè)需要計(jì)算（數(shù)值試驗(yàn)）的方案，如圖9所示。對(duì)各方案分別進(jìn)行3維流場(chǎng)計(jì)算，得到渦輪級(jí)效率計(jì)算結(jié)果，見表4。帶有分段樣條曲線積疊線渦輪導(dǎo)葉的渦輪級(jí)（3維）如圖10所示。

表3 積疊線參數(shù)優(yōu)化設(shè)計(jì)取值

表4 流場(chǎng)分析計(jì)算結(jié)果

對(duì)表4中數(shù)據(jù)進(jìn)行分析，結(jié)果見表5。從表5中可見，對(duì)效率影響最大的幾何參數(shù)為XT（即尖部直線段部分所占高度的比例），其它影響較大的參 數(shù) 還 有 XH、AH、BT，而參數(shù)D的影響最小。

根據(jù)正交設(shè)計(jì)方法和計(jì)算結(jié)果可以確定下一輪優(yōu)化推薦基準(zhǔn)取值，見表6。

表5 積疊線參數(shù)重要性分析

表6 經(jīng)過(guò)1輪優(yōu)化計(jì)算后的推薦基準(zhǔn)取值

經(jīng)過(guò)多輪次尋優(yōu)計(jì)算，可以改善葉片積疊線與葉片基元型線的匹配關(guān)系，達(dá)到提高渦輪級(jí)效率的目的。該方法還可以應(yīng)用于掠形葉片設(shè)計(jì)。

3 積疊線曲線生成程序

為便于工程設(shè)計(jì)應(yīng)用，采用上述方法編寫了積疊線曲線生成代碼。為保證程序使用范圍，積疊線生成程序保留了傳統(tǒng)的拋物線積疊方法，作為2次曲線，自變量選擇葉片根部和尖部的傾角。代碼中也包含了積疊線生成程序與葉片造型程序、渦輪流場(chǎng)3維計(jì)算軟件的接口程序，同時(shí)也包含了正交回歸優(yōu)化過(guò)程的數(shù)據(jù)處理程序。結(jié)合積疊線生成程序、葉片設(shè)計(jì)軟件和渦輪流場(chǎng)3維計(jì)算軟件，目前已用于多型渦輪導(dǎo)葉設(shè)計(jì)中。

4 從尋優(yōu)解中選擇最終方案

通過(guò)對(duì)空間內(nèi)多個(gè)方案進(jìn)行尋優(yōu)分析，可以在眾多方案中鎖定可供選擇的幾個(gè)方案。而如何選擇最終方案，需要綜合考慮目標(biāo)函數(shù)的計(jì)算結(jié)果和設(shè)計(jì)對(duì)象的具體設(shè)計(jì)要求。比如，在如圖11所示的設(shè)計(jì)結(jié)果曲線中，按照數(shù)學(xué)優(yōu)化方法，設(shè)計(jì)結(jié)果將選擇A點(diǎn)，如果設(shè)計(jì)要求為A1，由于B點(diǎn)沒有達(dá)到，則A點(diǎn)是1個(gè)合適的設(shè)計(jì)結(jié)果；但是如果設(shè)計(jì)要求為B1，那么選擇A點(diǎn)將存在不足，因?yàn)锳點(diǎn)右側(cè)曲線段目標(biāo)值衰減過(guò)快，假設(shè)累計(jì)誤差符合正態(tài)分布，且在選擇A或B點(diǎn)情況下誤差分布相同，與A點(diǎn)相比，選擇B點(diǎn)將提高產(chǎn)品合格率，所以最終方案應(yīng)該選擇B點(diǎn)。也就是說(shuō)，受多種因素影響的設(shè)計(jì)方案，需要參考計(jì)算結(jié)果分布曲線，并綜合考慮各方面影響做出選擇。

在渦輪彎扭葉片設(shè)計(jì)中，無(wú)論采用何種積疊線形式和優(yōu)化方式，都需要根據(jù)參數(shù)化和優(yōu)化的具體情況，綜合考慮氣動(dòng)、結(jié)構(gòu)、強(qiáng)度等方面的具體要求，做出合理選擇。

5 結(jié)論

（1）采用圓弧型、拋物線型、分段樣條曲線積疊線規(guī)律的渦輪彎扭葉片設(shè)計(jì)方法，可以用于設(shè)計(jì)有特定要求的彎扭葉片，不僅可以提高渦輪效率，還可以解決工程設(shè)計(jì)中存在的一些實(shí)際問(wèn)題；

（2）采用回歸正交設(shè)計(jì)方法對(duì)分段樣條曲線積疊線的渦輪彎扭葉片進(jìn)行優(yōu)化設(shè)計(jì)，可以有效提高彎扭葉片設(shè)計(jì)水平，貼近最優(yōu)解；

（3）根據(jù)本文所述方法編寫的積疊線曲線生成代碼，可以自動(dòng)生成圓弧型、拋物線型、分段樣條曲線積疊線，減少了設(shè)計(jì)工具的人工因素影響，縮短了設(shè)計(jì)周期，有效提高了設(shè)計(jì)效率，可滿足工程設(shè)計(jì)需要。

[1]Shahrokh Shahpar.Three-dimensional design and optimisation of turbomachinery blades using the Navier-Stokes equations[R].ISABE-2001-1053.

[2]王仲奇.采用彎扭靜葉片的燃?xì)馔钙皆O(shè)計(jì)計(jì)算方法[R].CSAA-90-PT-12.

[3]Chen N X.An exhaustive studyon skewingeffect ofsabre-shape blades in a turbine stator[R].ISABE-95-7055.

[4]陳光.“湍達(dá)”發(fā)動(dòng)機(jī)設(shè)計(jì)特點(diǎn)[J].國(guó)際航空，1991（3）：40-41.

[5]陳立周，周培德，高云章.機(jī)械優(yōu)化設(shè)計(jì)方法[M].北京：冶金工業(yè)出版社，1985：50-52.

[6]朱偉勇，關(guān)穎男.最優(yōu)設(shè)計(jì)理論與應(yīng)用[M].沈陽(yáng)：遼寧人民出版社，1981：146-249.

[7]數(shù)學(xué)手冊(cè)編寫組.數(shù)學(xué)手冊(cè)[M].北京：人民教育出版社,1979：809-859.

Aerodynamic Design Engineering Method of Bending-torsion Turbine Blade

ZHAO Da-yong
（AVIC Shengyang Aeroengine Research Institute,Shenyang 110015,China）

Bending-torsion blade is one of the main methods to improve turbine efficiency,construction of bending-torsion stacking line parametric mode is the important base for design of the turbine bending-torsion blade.To meet the engineering requirement of various turbine bending-torsion blade,the design method about arc,parabola,subsection spline curve stacking line rules were introduced,stacking line generating code was complied,and the application conditions and methods were suggested.The regression orthogonality design method can optimize the design of turbine bending-torsion blade with subsection spline curve stacking line.The bending-torsion balde engineering design method has been used in several turbine aerodynamic design.

bending-torsion blade;optimization;turbine;aerodynamic design;stacking line；aeroengine

趙大勇（1970），男，高級(jí)工程師，從事渦輪性能設(shè)計(jì)工作。