半水冷水輪發(fā)電機定子線圈穩(wěn)態(tài)溫度場分析

童旭松，宋德強，劉平超

?

半水冷水輪發(fā)電機定子線圈穩(wěn)態(tài)溫度場分析

童旭松1，宋德強2，劉平超1

（1. 哈爾濱電機廠有限責(zé)任公司，哈爾濱 150040；2. 黑龍江牡丹江抽水蓄能有限公司，黑龍江 牡丹江 150070）

本文利用有限體積法對一臺460MW半水冷水輪發(fā)電機定子線圈穩(wěn)態(tài)溫度場進行了共軛傳熱計算研究，減少了經(jīng)驗性參數(shù)對計算結(jié)果的影響。計算結(jié)果表明：水溫、銅溫沿水流方向逐漸增加；線圈中各截面水、銅平均溫度及銅最高溫度與截面和水路入口的間距成近似線性關(guān)系，且三者沿水流方向變化程度相近。

半水冷水輪發(fā)電機；穩(wěn)態(tài)溫度場；有限體積法；共軛傳熱

0 前言

作為發(fā)電機設(shè)計的關(guān)鍵指標(biāo)之一，溫升的大小直接關(guān)系到發(fā)電設(shè)備的壽命及運行可靠性，所以良好的冷卻系統(tǒng)起著至關(guān)重要的作用。目前主要的冷卻方式有空冷[1-2]、氫冷[3-4]、水冷[5-6]及蒸發(fā)冷卻[7-8]等方式。本文中的水輪發(fā)電機定子線圈電流密度高，宜采用半水冷冷卻方式，即定子線圈為水冷方式，其他各部件由空氣冷卻。當(dāng)前對于水冷、半水冷冷卻技術(shù)及其冷卻作用下定轉(zhuǎn)子溫度場計算的研究較多：魯長彬等對水內(nèi)冷定子溫度場進行了有限元分析，并以劉家峽機組為例將計算結(jié)果與實驗數(shù)據(jù)進行了比較，二者較為吻合，驗證了有限元方法的可靠性[9]。董韶峰等通過計算的方式研究了水內(nèi)冷汽輪發(fā)電機的定子溫度場分布，并模擬了水路堵塞對定子溫升的影響[10]。王紅宇等分別采用了熱-流體耦合有限元方法和網(wǎng)絡(luò)模型方法對三峽定子溫度場進行了分析研究[11-12]。除了在發(fā)電設(shè)備領(lǐng)域的研究，和偉超、李子健及Ping Zheng等對于水冷技術(shù)在電動車領(lǐng)域的應(yīng)用也進行了一系列關(guān)于冷卻系統(tǒng)設(shè)計及溫升計算的研究[13-15]。

就目前的研究情況可以發(fā)現(xiàn)，對于采用半水冷方式的定子線圈溫度場研究主要還是通過有限元及網(wǎng)絡(luò)模型的方法進行的，即使是采用熱-流體耦合的方法也由于模型尺寸及網(wǎng)格劃分等限制只能進行局部分析。本文則是假設(shè)定子線圈產(chǎn)生的損耗全部由水帶走，利用有限體積法對一條水路長度上的銅與水之間的共軛傳熱進行了直接的整場分析計算，減少了人為經(jīng)驗參數(shù)對計算結(jié)果的影響。

1 模型、網(wǎng)格劃分及計算方法

由于定子線棒采用水冷方式，計算時假設(shè)線棒產(chǎn)生的損耗全部由水帶走。線棒水路由三根上層線棒和三根下層線棒組成。為方便建模及網(wǎng)格劃分，建模時將線棒展開拉直且不考慮匝間的沿程換位，其中上層線棒的拉直長度為4.88m，下層線棒的拉直長度為4.92m，總長度為29.39m。線棒橫截面尺寸如圖1，由于其截面為典型的對稱結(jié)構(gòu)，所以為了節(jié)省計算資源，取其1/4作為計算域，即圖中虛線標(biāo)示處。

將圖1中計算域沿水路方向拉伸29.39m即得到對應(yīng)的三維模型，如圖2（a）所示，建模時匝間絕緣利用無厚度面代替。為提高網(wǎng)格質(zhì)量以保證計算準(zhǔn)確性，對實體模型進行網(wǎng)格劃分時均采用六面體結(jié)構(gòu)網(wǎng)格。計算中考慮到對流體附面層的求解以及銅與流體的換熱作用，所以在近壁面對網(wǎng)格進行了加密，第一層網(wǎng)格節(jié)點與壁面的距離為0.1mm，網(wǎng)格總數(shù)約為450萬，網(wǎng)格剖分如圖2（b）所示。

圖1 線棒橫截面尺寸圖

圖2 模型及網(wǎng)格剖分示意圖

共軛傳熱問題的求解在Fluent中完成，且認(rèn)為此傳熱過程為定常的。因為涉及傳熱問題，所以需要對能量方程進行求解。湍流模型選用Realizable k-ε，對于粘性底層的處理，這里采用Standard Wall Functions壁面函數(shù)。

邊界條件如圖3所示，水路入口為質(zhì)量流量入口，對應(yīng)側(cè)水路出口為壓力出口，并將Target Mass Flow Rate設(shè)置成與入口處相同的數(shù)值，以保證水路流量為定值。水路中間面設(shè)置為對稱邊界條件。不同材料間的每對接觸面設(shè)置為Coupled耦合邊界。建模時，線圈中的匝間絕緣利用無厚度面代替，在求解器中將其設(shè)置為0.5mm厚的絕緣材料即可，采用這種設(shè)置及求解方法可以帶來三個好處：第一，可以有效降低建模及網(wǎng)格劃分的復(fù)雜性；第二，能夠減少網(wǎng)格數(shù)量；第三，不用考慮銅材料傳熱的各向異性問題。圖中其余外表面均為絕熱面。邊界條件及線圈損耗設(shè)定值見表1。

表1 邊界條件及線圈損耗

圖3 邊界條件示意圖

2 計算結(jié)果

圖4給出了沿水流方向不同截面的溫度分布云圖，其中=0m和=29.39m兩個截面分別為水路入口端及出口端。從云圖中可以看出，水流和線圈溫度在入口端溫度均最低，隨著水流沿程逐漸升溫?fù)Q熱性能變差，線圈溫度沿程也在升高，且在出口端達(dá)到最大值。此外，如圖4所示，圈出位置即各截面的高溫區(qū)，結(jié)合圖1及圖2（a）可發(fā)現(xiàn)此高溫區(qū)發(fā)生在頂匝空心線圈相鄰的實心線圈處，且靠近主絕緣拐角。高溫區(qū)域均發(fā)生在此位置的原因有兩個：第一，頂匝空心線圈為單匝，其他3處中部空心線圈均為雙匝，這導(dǎo)致了頂匝空心線圈中水量僅為其他中部空心線圈水量的一半而使其帶走的熱量少；第二，頂匝空心線圈位于中部，而其他空心線圈均布于兩側(cè)，這使二者之間的熱量傳遞不均勻，靠近主絕緣位置處的熱量不易被水帶走。

從云圖中定性直觀地分析了沿水路方向（即Y方向）各截面溫度的變化情況。下面將各截面中的溫度變化進行定量分析，如圖5所示。圖5中給出了沿水流方向不同截面水的平均溫度（加權(quán)平均值）、銅的平均溫度（加權(quán)平均值）以及銅的最高溫度值，各溫度值的變化趨勢與云圖中的變化情況一致，均有所升高。其中出口的水溫平均值、銅溫平均值及銅最高溫度值相對入口的分別增長了35.65%、32.54%、29.51%。此外，從圖中可以看出各截面溫度同該截面與水路入口的距離成近似的線性關(guān)系，水平均溫度變化、銅平均溫度變化、銅最高溫度變化的斜率分別為0.61、0.61、0.59，三者斜率相近，即變化程度相近，截距即為入口處各自溫度值。

圖4 沿水流方向不同截面溫度云圖分布

圖5 各截面水平均溫度、銅平均溫度及銅最高溫度

綜上所述，線圈中水、銅平均溫度及銅最高溫度沿水路流向均會增加，各值與截面和水路入口間距成近似線性關(guān)系，且三者沿程變化程度相近。

3 結(jié)論

利用Fluent進行定子線圈水冷共軛傳熱計算可以直接對水與銅之間的換熱進行求解，減小了人為輸入經(jīng)驗性參數(shù)對結(jié)果產(chǎn)生影響的問題，并且通過建立無厚度面來模擬匝間絕緣，從而不用考慮銅傳熱的各向異性問題。從計算結(jié)果的分析還得到了以下結(jié)論：

（1）水溫、銅溫沿水流方向會逐漸增加；

（2）各截面線圈中水、銅平均溫度及銅最高溫度與此截面和水路入口的間距成近似線性關(guān)系；

（3）同一線圈截面中水、銅平均溫度及銅最高溫度相對入口截面處變化程度相近。

[1] 鮑里先科, 丹科. 電機中的空氣動力學(xué)與熱傳遞[M]. 北京: 機械工業(yè)出版社, 1985.

[2] 溫嘉斌, 孟大偉, 周美蘭, 等. 大型水輪發(fā)電機通風(fēng)發(fā)熱場模擬研究及通風(fēng)結(jié)構(gòu)優(yōu)化計算[J]. 電工技術(shù)學(xué)報, 2000, 15(6): 1-4.

[3] 李德基, 曹國宣. 汽輪發(fā)電機氫氣直接冷卻轉(zhuǎn)子的三維溫度場計算[J]. 電工技術(shù)學(xué)報, 1986, (3): 3-10.

[4] 胡晟, 金如麟, 譚茀娃. 1000MW級軸向氫冷汽輪發(fā)電機定子鐵心通風(fēng)系統(tǒng)[J]. 上海交通大學(xué)學(xué)報, 2002, 36(7): 1004-1007.

[5] 吳琳, 王宏光. 水冷電機冷卻系統(tǒng)設(shè)計與計算[J]. 機械設(shè)計與制造, 2008, (8): 40-41.

[6] 吳桂珍, 孟大偉, 許明宇. 高能量密度水冷電機冷卻系統(tǒng)設(shè)計與熱力計算[J]. 防爆電機, 2008, 43(142): 1-4.

[7] 顧國彪. 蒸發(fā)冷卻技術(shù)的發(fā)展與應(yīng)用[J]. 電器工業(yè), 2003, (1): 5-6.

[8] 董海虹, 顧國彪. 汽輪發(fā)電機蒸發(fā)冷卻系統(tǒng)的熱力過程及評價[J]. 中國電機工程學(xué)報, 2008, 28(20): 137-141.

[9] 魯長彬, 苗立杰. 水輪發(fā)電機水內(nèi)冷定子溫度場的分析方法[J]. 大電機技術(shù), 2000, (6): 1-4.

[10] 董韶峰, 張少鋒, 李蔭堂, 等. 發(fā)電機定子水冷穩(wěn)態(tài)的方針分析[J]. 電力與能源, 2011, 32(2): 131-133.

[11] 王紅宇, 蘇鵬聲, 王祥珩. 采用熱-流體耦合單元計算三峽半水冷水輪發(fā)電機定子三維穩(wěn)態(tài)溫度場[J]. 大電機技術(shù), 2006, (6): 1-12.

[12] 王紅宇, 李和明, 羅應(yīng)立, 等. 三峽水輪發(fā)電機全定子非線性熱-流體耦合網(wǎng)絡(luò)模型研究[J]. 中國電機工程學(xué)報, 2008, 28(8): 136-142.

[13] 和偉超. 電動汽車用永磁同步電機水冷系統(tǒng)設(shè)計及溫升分析[D]. 杭州：浙江大學(xué), 2011.

[14] 李子健. 電動車用水冷感應(yīng)電機溫度場及其水冷系統(tǒng)的研究[D]. 哈爾濱：哈爾濱工業(yè)大學(xué), 2008.

[15 ] Ping Zheng, Ranran Liu, P. Thelin, et al. Research on the Cooling System of a 4QT Prototype Machine Used for HEV[J]. IEEE Transactions on Energy Conversion, 2008, 23(1): 61-67.

Analysis on the Steady-state Temperature Field of a Half-water-cooled Hydro-generator Stator Coil

TONG Xusong1, SONG Deqiang2, LIU Pingchao1

(1.Harbin Electric Machinery Company Limited, Harbin 150040, China;2. Heilongjiang Mudanjiang Pumped Storage Co., Ltd., Mudanjiang 157000, China)

In this paper, a conjugate heat transfer on the steady-state temperature field of a 460MW half-water-cooled hydro-generator stator coil is studied by using finite volume method, and the impact of empirical parameters is effectively decreased. The results show that the temperatures of water and copper gradually increase along the direction of water flow. The average temperature of cooling water and copper, and the maximum temperature of copper in each coil cross section have an approximate linear relation with the distance between the inlet and the cross section, and they change similarly along the direction of the water flow.

half-water-cooled hydro-generator; steady state temperature field; finite volume method; conjugate heat transfer

TM312

A

1000-3983(2018)03-0024-03

2017-08-10

童旭松（1972-），2005年畢業(yè)于哈爾濱理工大學(xué)，工程碩士，主要從事大型電機技術(shù)研究，高級工程師。