童旭松，宋德強(qiáng)，劉平超

?

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

童旭松1，宋德強(qiáng)2，劉平超1

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

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

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

0 前言

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

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

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

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

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

圖1 線棒橫截面尺寸圖

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

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

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

表1 邊界條件及線圈損耗

圖3 邊界條件示意圖

2 計(jì)算結(jié)果

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

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

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

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

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

3 結(jié)論

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[14] 李子健. 電動(dòng)車用水冷感應(yīng)電機(jī)溫度場(chǎ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é)，工程碩士，主要從事大型電機(jī)技術(shù)研究，高級(jí)工程師。