尾水管渦作用下水輪機(jī)轉(zhuǎn)輪葉片動(dòng)水壓力特性

李兆軍，王玉江，劉福秀，黃鈺鈺，毛息軍

(1.廣西大學(xué) 機(jī)械工程學(xué)院，廣西 南寧 530004；2.廣西科技大學(xué) 廣西汽車(chē)零部件與整車(chē)技術(shù)重點(diǎn)實(shí)驗(yàn)室，廣西 柳州 545006)

0 引言

混流式水輪機(jī)彎肘型尾水管在部分負(fù)荷工況下會(huì)產(chǎn)生尾水渦流，造成轉(zhuǎn)輪葉片動(dòng)水壓力中存在著異常壓力脈動(dòng)[1]。這種異常壓力脈動(dòng)可能使轉(zhuǎn)輪葉片產(chǎn)生強(qiáng)烈振動(dòng)，進(jìn)而造成轉(zhuǎn)輪葉片出現(xiàn)疲勞裂紋，嚴(yán)重危害機(jī)組的安全運(yùn)行[2-3]。因而在研究水輪機(jī)轉(zhuǎn)輪葉片動(dòng)力學(xué)特性的過(guò)程中，有必要研究尾水管渦對(duì)轉(zhuǎn)輪葉片動(dòng)水壓力的影響。

目前有關(guān)尾水管渦對(duì)轉(zhuǎn)輪葉片動(dòng)水壓力的影響研究，主要通過(guò)計(jì)算流體力學(xué)(computational fluid dynamics，CFD)數(shù)值模擬和試驗(yàn)測(cè)試分析的方法。通過(guò)數(shù)值模擬方法，鐘林濤等[4]研究了轉(zhuǎn)輪出口速度、出口旋流數(shù)和渦帶形狀等因素與尾水管渦帶的關(guān)系。FOROUTAN等[5]研究了泄水錐射流與尾水管中心的軸向動(dòng)量以及尾水管渦的死水區(qū)的關(guān)系。通過(guò)試驗(yàn)測(cè)試方法，王小龍等[6-7]研究了轉(zhuǎn)輪葉片入口處的壓力脈動(dòng)幅值和頻率。以上兩種方法研究了不同工況時(shí)的動(dòng)水壓力變化規(guī)律，沒(méi)有得到反映動(dòng)水壓力與水力參數(shù)和結(jié)構(gòu)參數(shù)的數(shù)學(xué)模型，難以進(jìn)行轉(zhuǎn)輪葉片動(dòng)水壓力的暫態(tài)過(guò)程分析[8]。而數(shù)學(xué)模型能夠反映動(dòng)水壓力與水力參數(shù)和結(jié)構(gòu)參數(shù)的內(nèi)在關(guān)系，可對(duì)轉(zhuǎn)輪葉片動(dòng)水壓力的暫態(tài)過(guò)程進(jìn)行分析，且暫態(tài)過(guò)程分析又是研究轉(zhuǎn)輪葉片動(dòng)力學(xué)特性的基礎(chǔ)工作，而目前有關(guān)此方面的研究鮮見(jiàn)文獻(xiàn)報(bào)道。

本文以水輪機(jī)轉(zhuǎn)輪葉片為研究對(duì)象，根據(jù)水輪機(jī)轉(zhuǎn)輪葉片準(zhǔn)三元反設(shè)計(jì)的特點(diǎn)，應(yīng)用有限單元法建立水輪機(jī)轉(zhuǎn)輪葉片的準(zhǔn)三元有限元模型，構(gòu)建尾水管渦作用下的轉(zhuǎn)輪葉片上動(dòng)水壓力的數(shù)學(xué)模型，并通過(guò)實(shí)例分析進(jìn)行驗(yàn)證。

1 水輪機(jī)轉(zhuǎn)輪葉片的準(zhǔn)三元有限元模型

圖1 水輪機(jī)轉(zhuǎn)輪葉片坐標(biāo)系Fig.1 Global coordinate system of a turbine runner blade

PJij(r,θ,z,t)=N1Pij+N2P(i+1)j+N3P(i+1)(j+1)+N4Pi(j+1)，

(1)

式中：Pij、P(i+1)j、P(i+1)(j+1)、Pi(j+1)為節(jié)點(diǎn)Cij、C(i+1)j、C(i+1)(j+1)、Ci(j+1)的尾水管作用下的壓力脈動(dòng)，N1、N2、N3、N4為節(jié)點(diǎn)Cij、C(i+1)j、C(i+1)(j+1)、Ci(j+1)動(dòng)水壓力的形函數(shù)，根據(jù)水輪機(jī)轉(zhuǎn)輪葉片準(zhǔn)三元理論的動(dòng)水壓力變化規(guī)律[9]，尾水管作用下的壓力脈動(dòng)形函數(shù)為

(2)

式中：a、b分別為單元長(zhǎng)度和寬度的一半。

2 尾水管渦作用下的水輪機(jī)轉(zhuǎn)輪葉片動(dòng)水壓力

水輪機(jī)機(jī)組在部分負(fù)荷時(shí)，由于水輪機(jī)流量的減小，水流在轉(zhuǎn)輪出口產(chǎn)生具有與轉(zhuǎn)輪旋轉(zhuǎn)方向相同的圓周分量，導(dǎo)致在轉(zhuǎn)輪出口的正環(huán)量增加。當(dāng)環(huán)量達(dá)到一定程度，水流在尾水管產(chǎn)生一個(gè)與轉(zhuǎn)輪旋轉(zhuǎn)方向相一致的偏心渦帶[4]。尾水管渦引起了轉(zhuǎn)輪葉片的出口回流，在回流作用的影響下，轉(zhuǎn)輪葉片節(jié)點(diǎn)Cij的壓力脈動(dòng)pij可表示為[10]

pij=kijAijsin(2πft+φ)，

(3)

式中：Aij為節(jié)點(diǎn)Cij所在截面尾水管渦的壓力脈動(dòng)幅值，φ為初相位，與初始時(shí)刻的葉片相對(duì)位置有關(guān)；kij為尾水管渦引起的回流作用下葉片節(jié)點(diǎn)Cij的動(dòng)水壓力系數(shù)，與節(jié)點(diǎn)位置參數(shù)有關(guān)，為

(4)

式中：ε為葉片翼型常數(shù)，且-1≤ε≤1，可根據(jù)轉(zhuǎn)輪葉片的具體翼型確定。

式(3)中f為尾水管頻率，在部分負(fù)荷工況下，尾水管頻率是轉(zhuǎn)頻和設(shè)計(jì)水頭的函數(shù)，可表示為[11]

(5)

式中：n為水輪機(jī)轉(zhuǎn)輪的轉(zhuǎn)速，Q為水輪機(jī)實(shí)際流量，Q0為水輪機(jī)額定流量，ra為中央流線出口處半徑，β為葉片出口水流角，S為轉(zhuǎn)輪出口過(guò)水?dāng)嗝婷娣e，且S=2πdra，d為出水邊葉片弧長(zhǎng)長(zhǎng)度。

根據(jù)布拉修斯定理[12]，式(3)中節(jié)點(diǎn)Cij所在截面尾水管渦的壓力脈動(dòng)幅值A(chǔ)ij可用尾水管渦帶圓周速度分量的動(dòng)量表示為

(6)

(7)

式中：rw為尾水管入口截面半徑;z0為泄水錐的深度;τ1、τ2為待定正常數(shù)，可根據(jù)試驗(yàn)結(jié)果的經(jīng)驗(yàn)數(shù)據(jù)確定。

根據(jù)文獻(xiàn)[14]可知，當(dāng)Q=0.9Q0，H=H0時(shí)，尾水管渦帶消失；而在Q=0.25Q0，H=H0部分負(fù)荷運(yùn)行區(qū)時(shí)，下環(huán)處尾水管渦帶半徑接近于尾水管入口截面半徑，依據(jù)具體的流量和水頭，可計(jì)算出τ1、τ2。

式(6)中渦帶偏心距e可根據(jù)仿真結(jié)果和試驗(yàn)數(shù)據(jù)獲得[2]

(8)

(9)

式中:Cv為尾水管渦引起的節(jié)點(diǎn)Cij圓周速度與節(jié)點(diǎn)Cij絕對(duì)速度圓周分量的相關(guān)系數(shù)，與尾水管渦帶運(yùn)行的工況有關(guān)。uij為節(jié)點(diǎn)Cij的絕對(duì)速度的圓周分量，可表示為

(10)

將式(4)～(10)代入到式(3)中，并整理得節(jié)點(diǎn)Cij的尾水管渦作用下轉(zhuǎn)輪葉片的動(dòng)水壓力為

(11)

由式(11)可知，尾水管渦作用下轉(zhuǎn)輪葉片的動(dòng)水壓力幅值隨著流量的增加而逐漸減小，隨著水頭的增加而增加；而尾水管渦作用下轉(zhuǎn)輪葉片的動(dòng)水壓力頻率隨著轉(zhuǎn)速的增加而逐漸增加，隨著水頭的增加而減少。

根據(jù)式(11)尾水管渦作用下轉(zhuǎn)輪葉片節(jié)點(diǎn)Cij的動(dòng)水壓力，可得到尾水管渦作用下的Jij單元另外3個(gè)節(jié)點(diǎn)C(i+1)j、C(i+1)(j+1)、Ci(j+1)的動(dòng)水壓力，分別表示為

(12)

(13)

(14)

將式(2)、(11)至(14)分別代入式(1)中，即可得到t時(shí)刻尾水管渦作用下轉(zhuǎn)輪葉片Jij單元中任意位置(r,θ,z)動(dòng)水壓力，并依據(jù)式(3)，可以表示為

(15)

式中：Bij、B(i+1)j、B(i+1)(j+1)和Bi(j+1)分別為尾水管渦作用下水輪機(jī)轉(zhuǎn)輪葉片節(jié)點(diǎn)Cij、C(i+1)j、C(i+1)(j+1)和Ci(j+1)的動(dòng)水壓力幅值，且

(16)

(17)

(18)

(19)

由式(15)可知，尾水管作用下轉(zhuǎn)輪葉片的瞬時(shí)動(dòng)水壓力是一個(gè)時(shí)間和空間的函數(shù)，不僅自身的位置參數(shù)有關(guān)，如：Z軸坐標(biāo)、半徑等；也與水力參數(shù)有關(guān)，如：流量、水頭、轉(zhuǎn)速等；而且與水輪機(jī)結(jié)構(gòu)參數(shù)有關(guān)，如：泄水錐長(zhǎng)度、尾水管入口截面半徑、轉(zhuǎn)輪出口過(guò)水?dāng)嗝婷娣e、葉片出口水流角等，反映了動(dòng)水壓力與結(jié)構(gòu)參數(shù)和水力參數(shù)之間的關(guān)系，可以通過(guò)改變水力參數(shù)進(jìn)行暫態(tài)過(guò)程分析，這也為進(jìn)一步研究尾水管渦作用在水輪機(jī)轉(zhuǎn)輪葉片上的動(dòng)力學(xué)特性提供理論依據(jù)。

3 實(shí)例分析

3.1 實(shí)例

根據(jù)目前的研究成果，尾水管渦作用下的轉(zhuǎn)輪葉片動(dòng)水壓力與運(yùn)行工況相關(guān)。動(dòng)水壓力頻率是轉(zhuǎn)頻和設(shè)計(jì)水頭的函數(shù)；部分負(fù)荷運(yùn)行區(qū)，渦帶偏心嚴(yán)重，形成圓柱螺旋渦，壓力脈動(dòng)幅值大；而在額定負(fù)荷附近時(shí)，渦帶不明顯甚至消失，壓力脈動(dòng)很小[16]。本文選取Francis-99混流式水輪機(jī)轉(zhuǎn)輪葉片為研究對(duì)象，驗(yàn)證數(shù)學(xué)模型的可行性，并對(duì)尾水管渦作用下的轉(zhuǎn)輪葉片動(dòng)水壓力暫態(tài)過(guò)程變化規(guī)律進(jìn)行研究。

根據(jù)準(zhǔn)三元反設(shè)計(jì)的特點(diǎn)，采用5條徑向線與9條軸向線將葉片分成45個(gè)單元，其中第1條徑向線為上冠處，第5條徑向線為下環(huán)處，則長(zhǎng)葉片上的各個(gè)節(jié)點(diǎn)如圖2所示。進(jìn)水邊的節(jié)點(diǎn)為C00,C01,…,C04；出水邊的節(jié)點(diǎn)為C80,C81,…,C84。

圖2 試驗(yàn)研究對(duì)象結(jié)構(gòu)示意圖Fig.2 A structural sketch of the runner blade

3.2 數(shù)學(xué)模型驗(yàn)證

根據(jù)式(7)可知，尾水管渦初始于泄水錐的尾部，故以超過(guò)泄水錐深度的節(jié)點(diǎn)C71為研究對(duì)象。通過(guò)式(11)計(jì)算出節(jié)點(diǎn)C71的動(dòng)水壓力，其動(dòng)水壓力值時(shí)域仿真曲線與文獻(xiàn)[17]的試驗(yàn)數(shù)據(jù)，如圖3所示。

由于測(cè)試的試驗(yàn)數(shù)據(jù)中壓力不僅包含了尾水管渦作用下的壓力脈動(dòng)，也包含了其他壓力脈動(dòng)，所以部分?jǐn)?shù)據(jù)誤差比較大，由圖3可知，節(jié)點(diǎn)C71處的平均誤差為2.84 %。最大誤差值為8.32 %，從數(shù)據(jù)對(duì)比中，可以證明公式(11)在計(jì)算節(jié)點(diǎn)動(dòng)水壓力的有效性，進(jìn)而說(shuō)明了數(shù)學(xué)模型的可行性。

通過(guò)式(15)即可得到任意位置處的動(dòng)水壓力，其中t=0.08 s時(shí)葉片任意位置的動(dòng)水壓力如圖4所示。從圖4中可以看出，瞬時(shí)的動(dòng)水壓力沿軸線方向逐漸減少，瞬時(shí)的動(dòng)水壓力沿流線方向也逐漸減少。與文獻(xiàn)[17]的尾水管渦作用下轉(zhuǎn)輪葉片動(dòng)水壓力CFD仿真圖相比較，瞬時(shí)時(shí)刻的最大誤差出現(xiàn)在葉片的出水邊近下環(huán)處，誤差為7.92 %。從數(shù)據(jù)對(duì)比中，可以證明公式(15)在計(jì)算任何位置動(dòng)水壓力的有效性，這也證明數(shù)學(xué)模型的有效性。

圖3 節(jié)點(diǎn)C71的試驗(yàn)值與仿真值Fig.3 Time domain simulation curves and the experimental data of node C71

圖4 尾水管渦作用下t=0.08 s時(shí)轉(zhuǎn)輪葉片動(dòng)水壓力分布圖Fig.4 Hydrodynamic pressure diagram of the runner blade at t=0.08 sunder the vortex rope in draft tube

3.3 動(dòng)水壓力暫態(tài)過(guò)程變化規(guī)律

為了研究尾水管渦作用下轉(zhuǎn)輪葉片的動(dòng)力學(xué)特性，對(duì)暫態(tài)過(guò)程中動(dòng)水壓力變化規(guī)律進(jìn)行研究。通過(guò)改變水輪機(jī)運(yùn)行工況，研究水力參數(shù)對(duì)尾水管渦作用下轉(zhuǎn)輪葉片節(jié)點(diǎn)C71的動(dòng)水壓力幅值的影響。

圖5為額定轉(zhuǎn)速和額定水頭時(shí)，不同流量下尾水管渦作用下轉(zhuǎn)輪葉片節(jié)點(diǎn)C71的動(dòng)水壓力幅值仿真曲線，從圖5中可以看出，當(dāng)轉(zhuǎn)速和水頭不變時(shí)，隨著流量的增加，尾水管渦作用下的轉(zhuǎn)輪葉片動(dòng)水壓力幅值逐漸減小；當(dāng)運(yùn)行工況高于額定功率的90.1 %時(shí)，尾水管渦帶的半徑接近于零，轉(zhuǎn)輪葉片動(dòng)水壓力幅值也接近于零，此時(shí)可忽略動(dòng)水壓力對(duì)轉(zhuǎn)輪葉片的影響。而在高負(fù)荷下，又重新出現(xiàn)圓柱形死水區(qū)，此時(shí)尾水管渦作用下轉(zhuǎn)輪葉片動(dòng)水壓力幅值小于6 Pa，也可以忽略不計(jì)。

圖5 不同流量下尾水管渦作用下的節(jié)點(diǎn)C71的動(dòng)水壓力幅值仿真曲線Fig.5 Simulation curve of the hydrodynamic pressure amplitude of node C71 for different discharge

圖6為額定轉(zhuǎn)速和額定流量時(shí)，不同水頭下尾水管渦作用下的轉(zhuǎn)輪葉片節(jié)點(diǎn)C71動(dòng)水壓力幅值仿真曲線。從圖6中可以看出，當(dāng)轉(zhuǎn)速和流量不變時(shí)，隨著水頭的增加，尾水管渦作用下的轉(zhuǎn)輪葉片的動(dòng)水壓力逐漸增加。當(dāng)運(yùn)行工況低于額定功率的26.4 %時(shí)，死水區(qū)幾乎充滿整個(gè)尾水管，動(dòng)水壓力幅值接近零，此時(shí)可忽略動(dòng)水壓力對(duì)轉(zhuǎn)輪葉片的影響。結(jié)合圖5和圖6可知，尾水管渦作用下轉(zhuǎn)輪葉片的動(dòng)水壓力主要產(chǎn)生在額定功率的26.4 %～90.1 %的運(yùn)行工況。

圖6 不同水頭下尾水管渦作用下的節(jié)點(diǎn)C71的動(dòng)水壓力幅值仿真曲線Fig.6 Simulation curve of the hydrodynamic pressure amplitude of node C71 for different water heads

以上研究表明，應(yīng)用本文所建立的尾水管渦作用下的轉(zhuǎn)輪葉片動(dòng)水壓力數(shù)學(xué)模型，對(duì)Francis-99混流式水輪機(jī)轉(zhuǎn)輪葉片進(jìn)行實(shí)例分析，驗(yàn)證了數(shù)學(xué)模型的可行性，并研究了尾水管渦作用下的水輪機(jī)轉(zhuǎn)輪葉片的動(dòng)水壓力暫態(tài)過(guò)程變化規(guī)律。

4 結(jié)論

根據(jù)水輪機(jī)轉(zhuǎn)輪葉片準(zhǔn)三元反設(shè)計(jì)的特點(diǎn)，應(yīng)用有限單元法建立了尾水管渦作用下的水輪機(jī)轉(zhuǎn)輪葉片的動(dòng)水壓力數(shù)學(xué)模型，揭示了尾水管渦作用下水輪機(jī)轉(zhuǎn)輪葉片的動(dòng)水壓力暫態(tài)過(guò)程變化規(guī)律，研究表明：

① 在暫態(tài)過(guò)程變化中，尾水管渦作用下的轉(zhuǎn)輪葉片動(dòng)水壓力主要產(chǎn)生在額定功率26.4 %～90.1 %的運(yùn)行工況；當(dāng)運(yùn)行工況低于額定功率的26.4 %時(shí)，死水區(qū)幾乎充滿整個(gè)尾水管，動(dòng)水壓力幅值接近于零；當(dāng)運(yùn)行工況高于額定功率的90.1 %時(shí)，尾水管渦帶的半徑接近于零，動(dòng)水壓力幅值小于6 Pa。

② 尾水管渦初始于泄水錐的尾部，尾水管渦作用下轉(zhuǎn)輪葉片動(dòng)水壓力自上冠處到下環(huán)處逐漸增大，且動(dòng)水壓力最大值出現(xiàn)在轉(zhuǎn)輪葉片下環(huán)處靠近出水邊的位置。

③ 通過(guò)改變數(shù)學(xué)模型中的轉(zhuǎn)輪葉片位置參數(shù)以及水頭、流量等水力參數(shù)，揭示了轉(zhuǎn)輪葉片復(fù)雜的動(dòng)水壓力進(jìn)行暫態(tài)過(guò)程變化機(jī)理，為進(jìn)一步研究尾水管渦作用下水輪機(jī)轉(zhuǎn)輪葉片的動(dòng)力學(xué)特性提供理論依據(jù)。