王鶯 葉菁

摘要：座環(huán)是水輪機的主要部件，由上環(huán)、下環(huán)和若干固定導(dǎo)葉組成。本文采用CFD數(shù)值模擬揭示了座環(huán)內(nèi)部流動規(guī)律，得到了最大設(shè)計水頭下座環(huán)受到的水壓力分布規(guī)律及數(shù)據(jù)。將壓力數(shù)據(jù)加載到UGNX7.5中的Nastran有限元分析模塊作為其約束條件，計算靜態(tài)剛強度，得出正常工況下座環(huán)以及其主要各部件的變形和應(yīng)力分布規(guī)律，其平均應(yīng)力和最大應(yīng)力都滿足強度要求，在材料許用范圍內(nèi)。并在此基礎(chǔ)上對影響座環(huán)的關(guān)鍵環(huán)板（上環(huán)板、下環(huán)板、筋板）尺寸進行正交試驗分析，找到環(huán)板厚度的優(yōu)化方法，從而為優(yōu)化座環(huán)結(jié)構(gòu)提供有效依據(jù)。

關(guān)鍵詞：水輪機座環(huán) ?有限元 ?正交分析 ?優(yōu)化

1 概述

水輪機的座環(huán)承受著水輪發(fā)電機組的重量、蝸殼上表面部分混泥土重量及過流面水壓力，故結(jié)構(gòu)上它要求有足夠的強度和剛度[1]。目前，中小機組常采用帶碟形邊的座環(huán)直接與金屬蝸殼連接，其特點是鋼板焊接，制造剛度好，與蝸殼的連接點遠離固定導(dǎo)葉中心，改善了受力情況，在上下環(huán)外圓焊接圓形導(dǎo)流板，改善流動條件。本文主要是針對這類座環(huán)進行有限元優(yōu)化設(shè)計，以往的有限元分析中，壓力場主要是根據(jù)機組水頭換算施加到座環(huán)，而此次先通過流場分析軟件解算得到壓力場作為邊界載荷，使數(shù)據(jù)分析更加真實、精確。

2 座環(huán)CFD流場分析

2.1 機組主要參數(shù)

型號：HLA531-LJ-145;額定功率：PN= 22338kW;額定轉(zhuǎn)速：n=600r/min;額定水頭：H=193.5m;最大水頭：Hmax=202.6m;升壓水頭：H'=265m;額定流量：Q=12.4m3/s;額定效率：η=95%;水輪機二維原理圖如圖1所示。

■

圖1 ?水輪機二維原理圖

2.2 物理模型

水輪機的整個流道由蝸殼、座環(huán)、活動導(dǎo)葉、轉(zhuǎn)輪及尾水管組成，其三維模型如圖2所示。定義結(jié)合面INTERFACE，對實際流體模型做了適當(dāng)簡化。離散方程采用SIMPLEC算法。同時計算域采用非結(jié)構(gòu)化網(wǎng)格。由于尾水流道比較大，為節(jié)省計算時間，只對轉(zhuǎn)輪進行了加密，整個計算區(qū)域的網(wǎng)格數(shù)如表1所示。

■

圖2 ?水輪機過流部件三維模型

表1 ?水輪機模型過流部件的網(wǎng)格數(shù)據(jù)

■

2.3 邊界條件

水輪機最大水頭工況下給定進口壓力，并對相應(yīng)的質(zhì)量流量給定進口速度，在出口定義OUTFLOW自由出流條件。在所有的壁面上，定義無滑移邊界條件。動靜之間的耦合采用多重旋轉(zhuǎn)坐標系隱式方法[2，3]。水輪機流道網(wǎng)格和邊界條件如圖3所示。

■

圖3 ?水輪機過流部件網(wǎng)格劃分和邊界條件

2.4 模擬結(jié)果

從水輪機流道內(nèi)的靜壓分布圖4可看出蝸殼進口壓力大。蝸殼及座環(huán)壓力分布均勻，靠近活動導(dǎo)葉處壓降不大，在座環(huán)內(nèi)表面得到的相應(yīng)壓力數(shù)據(jù)將用于后續(xù)計算。蝸殼內(nèi)部流速分布均勻（如圖5所示），進出口速度環(huán)量從小到大變化明顯。

■ ? ■

圖4 ?水輪機流道靜壓分布圖 ?圖5 ?水輪機流道速度云圖

3 座環(huán)有限元分析

3.1 模型及網(wǎng)格

座環(huán)的結(jié)構(gòu)是旋轉(zhuǎn)周期對稱型，且根據(jù)CFD流場規(guī)律看到座環(huán)內(nèi)表面壓力沿轉(zhuǎn)輪半徑方向線性一致，故可取1/6結(jié)構(gòu)進行有限元模分析，如圖6所示。按表2的數(shù)據(jù)對部件進行網(wǎng)格劃分，總單元大小為213352個，足以達到精確要求。

■

圖6 ?1/6座環(huán)模型圖

表2 ?網(wǎng)格基本數(shù)據(jù)

■

3.2 邊界條件

正常工況時，座環(huán)承受：導(dǎo)水機構(gòu)自重（除底環(huán)）、頂蓋、底環(huán)傳來的力和自重，蝸殼水壓力和上下環(huán)水壓力[4]。根據(jù)CFD計算結(jié)果及分布規(guī)律進行加載。加載后的邊界條件如圖7所示。

■

圖7 ?座環(huán)邊界條件

3.3 有限元計算結(jié)果及分析

由座環(huán)的應(yīng)力云圖（如圖8）可看出：上環(huán)碟形邊在蝸殼水壓力作用下出現(xiàn)較大應(yīng)力，而該碟形邊與蝸殼焊接處接觸良好;固定導(dǎo)葉進水邊頭部應(yīng)力較大，是無焊接圓角過渡造成的。從座環(huán)的位移云圖（如圖9）可看出：上環(huán)碟形邊與蝸殼焊接處出現(xiàn)較大變形，主要是蝸殼上部與混泥土層之間非剛性固結(jié)，而是存在彈性層。受其影響，使得蝸殼上部出現(xiàn)上浮，釋放湍流脈沖壓力。

■■

圖8 ?座環(huán)應(yīng)力云圖 ? ? ? ? ? 圖9 ?座環(huán)位移云圖

對上環(huán)、固定導(dǎo)葉和下環(huán)進行有限元分析，可得出上環(huán)最大應(yīng)力為289.3MPa，主要出現(xiàn)在筋板與上碟形邊焊接處。固定導(dǎo)葉的最大應(yīng)力237.7MPa，接近材料屈服極限275MPa，但由于固定導(dǎo)葉實際焊接過程中會存在過渡圓角，其最大應(yīng)力值將會降低。但是上下環(huán)的環(huán)板應(yīng)力只有50MPa左右，故通過上述有限元分析方法，可以對座環(huán)的關(guān)鍵環(huán)板尺寸進行優(yōu)化。

4 上下環(huán)及筋板壓力正交試驗分析

由于固定導(dǎo)葉的結(jié)構(gòu)和數(shù)量與轉(zhuǎn)輪型號、直徑、水利性能有關(guān)，其高度為定值，故可轉(zhuǎn)而對上環(huán)板尺寸、下環(huán)尺寸和筋板進行正交試驗分析，找到這三者中影響座環(huán)強度和剛度的主要因素，尋找最優(yōu)水平組合，從而對環(huán)板或筋板的厚度進行優(yōu)化。本文采用三因子三水平方法[5]測試，各環(huán)板對應(yīng)名稱見圖10。試驗尺寸如表3所示。

■

圖10 ?環(huán)板位置示意圖

表3 ?各環(huán)板試驗尺寸

■

為減少試驗次數(shù)，選擇9次試驗正交表，相應(yīng)試驗結(jié)果見表4，數(shù)據(jù)分析見表5。

表4 ?試驗結(jié)果數(shù)據(jù)

■表5 ?試驗極差數(shù)據(jù)

■

對試驗結(jié)果數(shù)據(jù)進行分析，整理得到極差表單，通過極差可以發(fā)現(xiàn)，上環(huán)板尺寸影響強度和剛度最顯著，主要是上環(huán)受到蝸殼水壓力和頂蓋浮力通過螺栓產(chǎn)生的拉力造成，故在優(yōu)化厚度時，盡可能不考慮，下環(huán)板影響座環(huán)的剛強度小，可以安全優(yōu)化;此次試驗最好的組合：滿足強度條件的是A3B3C3，即：上環(huán)板35mm，筋板25mm，下環(huán)板35mm;滿足剛度條件的最優(yōu)組合是A2B3C2，即：上環(huán)板30mm，筋板25mm，下環(huán)板30mm;筆者結(jié)合實際工作經(jīng)驗推薦采用：保證強度的組合A3B3C3，而后驗算其剛度是否在材料的許可范圍即可。

5 結(jié)論

①應(yīng)用CFD數(shù)值模擬，得到座環(huán)靠近蝸殼進口壓力大，中間壓力分布均勻，進出口速度環(huán)量從小到大變化明顯的規(guī)律。

②目前的座環(huán)有限元分析，壓力場主要是根據(jù)機組水頭換算施加到座環(huán)，而本文通過流場分析軟件得到壓力場作為邊界載荷，使數(shù)據(jù)分析更加真實、精確。

③和傳統(tǒng)算法結(jié)果相比，此分析更直觀的得到座環(huán)各部件的應(yīng)力分布情況，并采用強度理論進行校核，發(fā)現(xiàn)強度在材料許用范圍內(nèi)，其他部件有優(yōu)化空間。

④對座環(huán)的上環(huán)板、下環(huán)板、筋板進行壓力正交試驗分析，得出影響座環(huán)強度和剛度的主要因素是上環(huán)，并研究得到最優(yōu)組合，該方法為優(yōu)化座環(huán)結(jié)構(gòu)提供有效依據(jù)。

參考文獻：

[1]嚴錦麗，徐志明.水輪機座環(huán)與蝸殼結(jié)構(gòu)剛強度靜力分析[J].水電站機電技術(shù)，2001（1）：12-13.

[2]王春林，司艷雷，鄭海霞，彭娜，趙佰通，張浩.旋流自吸泵內(nèi)部流場的數(shù)值模擬[J].排灌機械，2008，26（2）：31-35.

[3]田愛民，許洪元.旋轉(zhuǎn)噴射泵集流管內(nèi)部流動計算[J].石油化工設(shè)備，2005.3.

[4]龐立軍，魏洪久.水輪機蝸殼座環(huán)的應(yīng)力分析與評定[J].大電機技術(shù)，2008.3：39-42.

[5]袁壽其，張金鳳，袁建平.正交試驗研究分流葉片主要參數(shù)對性能影響[J].排灌機械，2008，26（2）：1-5.

[6]JamesR.Farr，MaanH.Jawad著.ASME壓力容器設(shè)計指南[M].鄭津洋，方曉斌等譯.第二版.北京：化學(xué)工業(yè)出版社，2003.

基金項目：

2013年度浙江省水利廳科技項目——中小型水電站水輪機座環(huán)的優(yōu)化設(shè)計研究（RC1315）。

作者簡介：

王鶯（1978-），女，浙江杭州人，副教授，主要從事CAD、CAM的教學(xué)及研究工作。