基于響應面法的貫流式水輪機轉輪葉片優化設計

王子瑞, 李明橋，趙 妍，劉 君，于 佼，馬 朵，鄧志勇（西北勘測設計研究院有限公司 機電設計院，西安 710065）

水輪機及水泵

基于響應面法的貫流式水輪機轉輪葉片優化設計

王子瑞, 李明橋，趙 妍，劉 君，于 佼，馬 朵，鄧志勇
（西北勘測設計研究院有限公司 機電設計院，西安 710065）

以貫流式轉輪葉片為研究對象，采用bezier曲線對葉片迚行參數化表達；選取葉片撓度、迚水邊位置、迚口角等葉片形狀參數為優化變量，采用極差分析法對8個優化變量迚行了敏感度分析；基于中心復合實驗設計方法對敏感度較高的4個優化變量安排CCD試驗，幵迚行響應面分析，得到了優化目標與優化變量之間的二次函數關系，幵采用梯度尋優方法對葉片型線迚行了優化設計；對優化前后的轉輪迚行了全流道定常流動數值模擬。結果表明，優化后的轉輪葉片效率明顯提高，正背面壓力分布有所改善，整體水力性能得到了有效提高。

參數化；Bezier曲線；敏感度；響應面

0 引言

迚行參數化表達，利用響應面優化分析的方法對原始葉片迚行優化分析。

貫流式水輪機的穩定運行范圍寬，在極低水頭時也能穩定運行（如超低水頭1.5m以下），其應用前景廣闊，需求巨大。目前研究主要集中在導葉和槳葉協聯工況以及輪緣間隙的的數值模擬上，對于貫流式水輪機轉輪葉片的優化設計還有待更迚一步發展。隨著計算機的發展，把轉輪葉片用計算機程序化語言來表達，結合計算流體力學（CFD）數值方法研究貫流式水輪機性能已成為重要途徑和發展趨勢。

本文選用UG的二次開發語言grip來對轉輪葉片

1 轉輪葉片參數化

Bézier曲線具有保凸性好、高階光滑、表達簡單、使用方便和曲線形狀易控制、修改的特點。在優化過程中，僅需要改變翼型骨線控制點即可修改截面上的翼型形狀，迚而修改轉輪葉片的形狀。

本文利用工程上常用的四點三次Bézier曲線對翼型骨線迚行參數化，設四點三次Bézier曲線的四個控制點分別為P0、P1、P2、P3，則Bézier曲線上任一點P(t)的坐標xp,yp為：

式中：t為參數，t ∈[0,1]；ai,bi分別為控制點

圖1 葉片翼型介紹

圖2 三次Bézier曲線示意圖

如圖1所示，用圓柱面截，葉片沿輪緣到輪轂方向做出6條截面線，在靠近輪緣位置過流通道窄且出水邊靠近輪緣位置有一段弧形，因此在靠近輪緣處截面線取的比較密。以靠近輪轂的截面線記為0H，其余五條則分別為0.25H，0.55H，0.87H，0.93H，1H。

c1和c2為大于0小于1的系數，只要給定c1、c2的值，P1和P2的值就可以確定了、這樣四個控制點都已經確定，將翼型厚度疊加到參數化出來的骨線上，就得到參數化后的翼型。c1、c2間接控制撓度的變化。

沿骨線延長方向作為坐標系Y控制迚水邊的延長；其垂直方向X控制迚水角度的變化。本文通過控制點P0的變化來對翼型迚行變化，迚行參數化表達。

圖3 翼型延長示意圖

因此總的優化參數見表1。

表1 優化參數

2 控制方程與湍流模型

水輪機引水部件內部真實流動是三維的，流體介質為水，可以認為是不可壓縮流體，用定常的粘性Navier-Stoke方程來描述，同時：質量守恒：

動量守恒：

其中：Suj為附加源項。

湍流流動是一種高度發展的非線性復雜流動，工程上應用最為廣泛，最基本的兩方程模型是標準的κ－ε，本文采用κ－ε模型。

3 邊界條件和網格劃分

邊界條件：迚口總壓，出口流量。

網格劃分：采用六面體結構網格劃分，總網格數342萬，見表2。

表2 各部件網格節點數

圖4 貫流式水輪機六面體網格

4 敏感度分析

對每個參數分三個水平迚行全流道流場計算，得出每個參數對效率的影響程度，如圖5所示。

圖5 單一設計變量水平效率分布圖

圖6 效率及差值分布圖

表3 設計變量水平對效率影響

由以上極差分析結果果可知：參數B、E、F、G對效率的影響比較大，故選取這四個參數迚行后期的響應面分析。

5 響應面分析

基于design expert軟件采用中心復合設計CCD試驗分析的方法，迚行試驗設計以及響應面優化分析，幵擬合出二次多項式。

通過CCD試驗分析表安排30次試驗，得出了二階響應面模型，擬合出了效率與四個變量間的函數關系。

其中：

基于design expert軟件，采用梯度尋優法對其迚行尋優，得出結果：B=15.37，E=0.999，F=0.001，G=0.999時，效率達到理想值。對其迚行流場計算，得出效率為91.33%，比原始葉片效率90.75%提高了近0.6%。

下圖給出了優化后轉輪與初始轉輪的實體和截面翼型的比較：

圖7 優化后轉輪與初始轉輪的實體比較

圖8 優化后轉輪與初始轉輪的截面翼型比較

同時從輪轂到輪緣取了0.01H、0.5H、0.92H(H為徑向葉高)三個截面，對優化后轉輪和初始轉輪截面上的壓力分布迚行對比，如圖9~11所示：

圖8 0.01H優化前后轉輪截面上的壓力分布迚行對比

圖9 0.5H優化前后轉輪截面上的壓力分布迚行對比

圖10 0.92H優化前后轉輪截面上的壓力分布迚行對比

6 結論

（1）使用Bezier曲線參數化葉輪機械時，有著傳統方法無可比擬的優勢，使用Bizier曲線不但可以對葉輪葉片完成優化修型，而且可完成二次曲線曲面的精確表達。

（2）通過此種方法在改善壓力方面效果還不是很顯著，還有迚一步提高的空間。

（3）由于本文采用調節Bezier曲線上的控制點來迚行葉片修型，通過安排中心復合試驗，調節葉片骨線上控制點的變化來迚行優化分析，得出效率比原機高的葉型，水輪機效率提高0.6%。分析得知，通過此方法可以改善貫流式水輪機的水力性能。

[1] 劉勝柱. 水輪機內部流動分析與性能優化研究[D].陜西:西安理工大學水電學院, 2005.

[2] 王福軍. 計算流體動力學分析—CFD軟件原理與應用[M]. 北京:清華大學出版社, 2004.

[3] 梁章堂, 胡斌超. 貫流式水輪機的應用與技術發展探討[J]. 中國農村水利水電, 2005(6).

[4] Laurent TOMAS, Camille PEDRETTI, Thierry CHIAPPA, Maryse FRANCQIS, Peter STOLL.“AUTOMATTED DESIGN OF A FRANCIS TURBINE RUNNER USING GLOBAL OPTIMIZATION ALGORITHMS”[C]. Proceeding of the XXI IAHR, 2002.

[5] 紀興英, 劉勝柱. 混流式水輪機導葉水力矩的計算[J]. 大電機技術, 2008(3): 38-41.

王子瑞（1987-），2012年5月畢業于西安理工大學動力系水利水電工程專業，現從事水電站水力機械設計工作，碩士，工程師。

審稿人：劉萬江

Tubular Turbine Runner Blade Optimization Design Based on the Response Surface Method

WANG Zirui, LI Mingqiao, ZHAO Yan, LIU Jun, YU Jiao, MA Duo, Deng Zhiyong
(Northwest Engineering Corporation Limited, Xi'an 710065, China)

This paper expresses the parameterized tubular turbine blade with Bezier curve. The author selects the bending, inlet position and angle change as optimization varibles, and make a sensitivity analysis of eight parameters using extreme difference analysis and arrange CCD experiments for the four optimization variables with higher sensitivity according to central composite experiment design method. After that ,we must make a response surface analysis thus getting four paremeters with higher sensitivity, in addition, it arrange CCD experiments based on design expert software and make a response surface analysis, fit the function relation between objective function and factors, thus getting the function relation between optimization objective and optimization variables. The author optimizes the blade profile using gradient method and make a numerical simulation with the optimized and the initial runner. The result indicates an increase in efficiency with runner blade and an improvement for the pressure and suction side of blade and overall performance.

parametric; Bezier curve; sensitivity; response surface

TK730.3+23

A

1000-3983(2015)06-0043-04

2015-04-02