彭茂林，劉 波，孫吉宏

（海軍裝備技術(shù)研究所，青島266012）

在將高溫高速燃氣的熱力學(xué)能轉(zhuǎn)換為渦輪機械能的過程中，渦輪葉片起到至關(guān)重要的作用。能量轉(zhuǎn)換的效率與葉片的形狀密切相關(guān)［1］，氣動效率最大化一直是研究者努力追求的目標［2－4］，但是氣動效率提高的同時往往會犧牲其它優(yōu)良性能，比如強度。因此，在對渦輪葉片進行氣動設(shè)計優(yōu)化的同時，如何保證葉片的強度，滿足可靠性要求，已成為難點重點問題。作者在文獻［5］中詳細介紹了渦輪葉片的參數(shù)化建模方法，在文獻［6］和文獻［7］中進一步研究了渦輪葉片氣動效率的計算方法和結(jié)構(gòu)可靠度計算方法，在文獻［8－9］中研究了基于響應(yīng)面法的可靠性穩(wěn)健設(shè)計優(yōu)化方法，在這些工作的基礎(chǔ)上，本文將葉片的結(jié)構(gòu)系統(tǒng)可靠度作為基本約束條件，對某船用渦輪葉片進行了三維氣動設(shè)計優(yōu)化，在提高效率的同時，保證了葉片的結(jié)構(gòu)可靠性，取得了較好的效果。

1 基于結(jié)構(gòu)可靠性的渦輪葉片氣動設(shè)計優(yōu)化模型

渦輪葉片氣動設(shè)計優(yōu)化模型見下式：

式中：f（μH，σH｜X）為渦輪葉片設(shè)計優(yōu)化的代價函數(shù)，μH，σH分別為結(jié)構(gòu)多性能目標函數(shù)的均值、標準方差向量，X為結(jié)構(gòu)的nxv維隨機變量向量；d為設(shè)計變量向量；dL，dU，σd分別為設(shè)計變量d的許可上、下限以及標準方差向量；gi（X；d）為結(jié)構(gòu)的第i個極限狀態(tài)函數(shù)；nc，nxv，ndv和βti分別為結(jié)構(gòu)可靠度約束、隨機變量、設(shè)計變量的個數(shù)以及第i個失效模式對應(yīng)的結(jié)構(gòu)可靠度指標目標值；P（gi（X；d）＜0）為gi（X；d）的概率；Φ（·）為標準正態(tài)累積概率分布函數(shù)；為連乘符號；βt0為多失效模式下的結(jié)構(gòu)系統(tǒng)可靠度指標；f為設(shè)計變量的穩(wěn)健性水平，以保證設(shè)計變量的最終優(yōu)化結(jié)果超出取值范圍的邊界的概率小于等于Φ（－f），通常取3≤f≤6。多目標可靠性穩(wěn)健設(shè)計代價函數(shù)f（μH，σH｜X）為各目標函數(shù)的加權(quán)和的形式，根據(jù)目標的不同特性而分為以下幾類：

第1類：對于望小特性：

第2類：對于望大特性：

第3類：對于望目特性：

式中：hit和分別為第i個目標函數(shù)的望目特性目標函數(shù)的目標均值和初始均值。

第4類：混合型

從表達式（2）－（5）可以看出，為了解決由于各個目標函數(shù)的量綱不同的問題，代價函數(shù)的各項都通過除以其初始設(shè)計值來進行無量綱化處理。權(quán)值的大小反映了設(shè)計者對于各項性能指標的偏好程度。多目標可靠性穩(wěn)健設(shè)計優(yōu)化代價函數(shù)不但尋求結(jié)構(gòu)性能目標均值的最優(yōu)，也期望結(jié)構(gòu)性能的標準方差能達到最小。而傳統(tǒng)設(shè)計優(yōu)化只追求性能目標的均值最優(yōu)解，但不考慮變量的隨機性。

2 渦輪葉片氣動設(shè)計優(yōu)化代價函數(shù)

2.1 性能目標函數(shù)的選取

1）對于靜葉片的設(shè)計優(yōu)化，選取靜葉氣動效率和靜葉片質(zhì)量為主要優(yōu)化目標，其中，氣動效率為望大特性目標函數(shù)，葉片質(zhì)量為望小特性目標函數(shù)；

2）對于動葉片，由于其在高溫高速氣流推動下高速旋轉(zhuǎn)，承受的應(yīng)力應(yīng)變及振動情況十分復(fù)雜，必須對其進行可靠性設(shè)計優(yōu)化。因此，選取動葉片的氣動效率作為望大特性目標函數(shù)，葉片質(zhì)量為望小特性目標函數(shù)，結(jié)構(gòu)可靠性為基本約束條件。

2.2 動葉片結(jié)構(gòu)可靠性約束

靜強度可靠性約束：

疲勞可靠性約束：

剛度可靠性約束：

式（6）－式（8）中：Smax、εmax、ξmax分別為動葉片運行中承受的最大應(yīng)力、最大應(yīng)變、最大變形；σs為葉片材料工作溫度下的屈服強度；ξs為葉片變形允許值；εP為渦輪設(shè)計安全壽命下概率－應(yīng)力應(yīng)變－疲勞壽命曲線對應(yīng)的循環(huán)應(yīng)變值；βts、βtP、βtξ分別為靜強度、疲勞、剛度可靠度指標的目標（約束）值。

2.3 渦輪葉片氣動設(shè)計優(yōu)化代價函數(shù)

式中：μE、μV分別為葉片氣動效率、葉片體積的均值；σE、σV分別為葉片氣動效率、葉片體的標準方差；μE0、μV0分別為μE、μV在葉片初始設(shè)計參數(shù)下的取值；σE0、σV0分別為σE、σV在葉片初始設(shè)計參數(shù)下的取值。本文取w11＝w12＝0.9，w21＝w22＝0.1。

3 渦輪葉片氣動設(shè)計優(yōu)化方法及優(yōu)化結(jié)果

3.1 優(yōu)化分析方法及步驟

1）首先基于正交試驗設(shè)計選取試驗點并建立相應(yīng)的實體及有限元模型，詳細建模過程請見參考文獻［5］，進行熱－流－固耦合有限元分析，分析流程如圖1所示；

2）根據(jù)有限元分析結(jié)果得到各試驗點下的葉片動能效率、最大應(yīng)力應(yīng)變、最大變形、體積等響應(yīng)參數(shù)，采用二次多項式函數(shù)擬合各性能函數(shù)和極限狀態(tài)函數(shù)對隨機變量的表達式。

3）基于響應(yīng)面函數(shù)，計算渦輪葉片的結(jié)構(gòu)可靠度，詳細分析方法參考文獻［8－9］。采用多目標序列二次規(guī)劃優(yōu)化算法對設(shè)計參數(shù)進行尋優(yōu)，得到最優(yōu)設(shè)計點；

4）在最優(yōu)點處再次進行熱－流－固耦合有限元分析，并將分析結(jié)果加入到試驗數(shù)據(jù)中，并更新響應(yīng)面模型，再次計算最優(yōu)設(shè)計點，如此循環(huán)迭代計算，當計算收斂到滿足精度要求時，優(yōu)化結(jié)果即為最優(yōu)設(shè)計。

圖1 渦輪葉片熱－流－固耦合有限元分析流程圖

3.2 優(yōu)化分析結(jié)果

優(yōu)化計算過程經(jīng)過6次迭代，結(jié)果達到收斂要求，歷時約120h。圖2為氣動設(shè)計優(yōu)化前后流場各截面燃氣壓力分布對比圖，從圖中可以看出，優(yōu)化后各截面燃氣流速分布比優(yōu)化前更為均勻，說明流動更平穩(wěn)，流動損失就會相對減少。圖3為設(shè)計優(yōu)化前后50%葉高截面葉片表面湍流運動能變化對比圖，湍流運動能的大小反映了流動損失的大小，從圖中可以看出，經(jīng)過設(shè)計優(yōu)化，靜動葉片的能量損失都得到明顯減少，以動葉片為例，50%葉高截面表面湍流運動能最高值比優(yōu)化前下降大約11.1%。圖4 為設(shè)計優(yōu)化前后流場及葉片模型對比圖。6.61%，靜、動葉片體積均值分別增加2.4%和0.64%。葉片體積增大，在制造葉片時原材料成本就會有所提高。

圖2 設(shè)計優(yōu)化前后流場各截面燃氣壓力分布對比圖

圖3 設(shè)計優(yōu)化前后50%葉高截面葉片表面湍流運動能變化對比圖

圖4 設(shè)計優(yōu)化前后流場及葉片模型對比圖

表1 渦輪靜葉片性能目標設(shè)計優(yōu)化結(jié)果

表1和表2分別為靜、動葉片性能目標設(shè)計優(yōu)化結(jié)果。從表中結(jié)果可知，經(jīng)過優(yōu)化后，靜、動葉的氣動效率均值分別提高了2.45%和6.17%，且氣動效率的標準方差也有一些減小。經(jīng)過設(shè)計優(yōu)化，動葉片的各項結(jié)構(gòu)可靠度指標都得到提高，且達到設(shè)計優(yōu)化對結(jié)構(gòu)可靠度的要求。優(yōu)化前，動葉片的結(jié)構(gòu)系統(tǒng)可靠度約為91.3%，而優(yōu)化后動葉片的結(jié)構(gòu)可靠度達到了99.87%。然而，優(yōu)化后靜、動葉片體積的均值和標準方差都有所增大，其中，靜、動葉片體積均值分別增加7.14%和

表2 渦輪動葉片性能目標設(shè)計優(yōu)化結(jié)果

4 結(jié)論

建立了渦輪葉片結(jié)構(gòu)可靠性三維氣動設(shè)計優(yōu)化模型和代價函數(shù)，采用多學(xué)科耦合方法對葉片進行有限元分析，并給出了設(shè)計優(yōu)化流程。最后，對其進行了多學(xué)科可靠性氣動設(shè)計優(yōu)化。結(jié)果表明，基于結(jié)構(gòu)可靠性的渦輪葉片三維氣動設(shè)計優(yōu)化在提高葉片氣動效率的同時，又可以保證葉片的結(jié)構(gòu)可靠性。因此，該方法具有較高的理論及工程應(yīng)用價值。但是，在保證葉片結(jié)構(gòu)可靠性、提高葉片氣動效率的同時，卻可能要在某些性能（如葉片體積和制造成本）上付出更高的代價，因此基于結(jié)構(gòu)可靠性的渦輪葉片三維氣動設(shè)計優(yōu)化過程也是一個權(quán)衡折中的過程。

［1］岳珠峰，李立州，王婧超.航空發(fā)動機渦輪葉片多學(xué)科設(shè)計優(yōu)化［M］.北京：科學(xué)出版社，2007.

［2］劉勛，周遜，溫風(fēng)波.F級中低熱值燃料氣輪機透平部分氣動方案優(yōu)化設(shè)計與分析［J］.中國電機工程學(xué)報，2011，31（26）：1－8.

［3］馬洪波，朱劍，席平.基于參數(shù)化的渦輪葉片三維氣動優(yōu)化仿真［J］.計算機仿真，2008，25（10）：27－30.

［4］虞跨海，岳珠峰，楊茜.渦輪冷卻葉片氣動與傳熱設(shè)計優(yōu)化［J］.計算力學(xué)學(xué)報，2010，27（4）：310－314.

［5］彭茂林，楊自春，曹躍云，等.基于貝賽爾曲線和粒子群算法的渦輪葉片型線參數(shù)化建模［J］.中國電機工程學(xué)報，2012，32（35）：101－108.

［6］彭茂林，楊自春，曹躍云.渦輪葉片三維氣動分析方法研究［J］.熱力透平，2012，41（4）：245－250.

［7］彭茂林，楊自春，曹躍云.渦輪盤片的結(jié)構(gòu)可靠度計算及靈敏度分析［J］.熱力透平，2011，40（3）：192－198.

［8］彭茂林，楊自春，曹躍云，等.渦輪葉片低周疲勞可靠性穩(wěn)健設(shè)計優(yōu)化研究［J］.中國電機工程學(xué)報，2013，33（11）：104－111.

［9］彭茂林，楊自春，曹躍云.基于響應(yīng)面法的可靠性穩(wěn)健設(shè)計優(yōu)化方法研究［J］.航空動力學(xué)報，2013，28（8）：1784－1790.