丁樹業(yè)， 葛云中， 陳衛(wèi)杰， 徐殿國， 苗立杰

(1.哈爾濱工業(yè)大學電氣工程及自動化學院，黑龍江哈爾濱150001;2.哈爾濱電站設備集團公司，黑龍江哈爾濱150040;3.哈爾濱理工大學電氣與電子工程學院，黑龍江哈爾濱150080;4.中國人民解放軍96215部隊，廣西柳州545616)

0 引言

我國風能資源十分豐富，地面的風能總量約有16億千瓦，其中近地面可供開發(fā)利用的約占16%，為2.53億千瓦［1］。在風能利用方面，我國至2010年國內風電裝機容量達到3 000萬千瓦，占國內總裝機容量的2.8%左右。目前，國內風力發(fā)電機的開發(fā)生產品種繁多［2］，單機容量不斷增大，2～3.0 MW以上的風電機組已經商業(yè)化，5 MW的風電機組也已經投入試運行，開始由陸地向海上和高海拔地區(qū)發(fā)展［3－4］。

隨著風力發(fā)電機單機容量的日益增大，使得電機運行時產生的單位體積損耗顯著增加，引起電機各部分溫度升高，直接影響電機的壽命和運行的可靠性。因此，在發(fā)電機的運行過程中，必須采用相應的冷卻介質有效地帶走電機內的熱量，使得發(fā)電機長期、可靠地運行。然而，發(fā)電機的通風結構以及冷卻介質在電機內的分布特性均比較復雜，因此對發(fā)電機內溫度場、溫度分布特性進行有效地分析，對發(fā)電機穩(wěn)定運行以及更大容量的發(fā)電機的設計均具有重要意義。

為了準確的對電機內溫度場進行研究，國內外專家學者提出了許多計算方法［5］，其中數值計算方法精度較高。目前運用數值計算的方法對不同種類電機內溫度場［6－10］、不同工況下鐵心與繞組間的熱量傳遞［11］、轉子通風方式和電機風量分配對電機溫升的影響進行了諸多研究［12－13］，為電機內的物理場的計算奠定了一定的基礎。上述文獻多是對電機的局部區(qū)域或相關因素對電機溫度場的影響情況進行研究，本文在此基礎上，對電機周向1/8區(qū)域內定子、轉子以及流體區(qū)域的三維溫度場進行研究，避免了對氣隙和轉子通風溝流量分配。

本文以一臺1.5 MW風力發(fā)電機為例，采用有限體積元法對電機內流體場和溫度場進行計算，得出內冷卻介質在電機內流量分布規(guī)律和特性，并對電機三維溫度場進行研究，分析電機溫升分布規(guī)律和傳熱特性;最后將耦合計算結果與實測結果進行比較分析。

1 數學模型

電機冷卻系統(tǒng)內流體流動與傳熱滿足質量守恒定律、動量守恒定律、能量守恒定律。如果流動處于紊流狀態(tài)，還要遵守附加的湍流運動方程，其守恒定律可以采用相應的控制方程進行描述［14－16］。

1.1 控制方程

根據流體力學理論，采用標準k－ε方程模型對流體進行描述。當流體為不可壓縮且處于穩(wěn)定流動狀態(tài)時，采用通用控制方程［17］為

式中:φ、V為通用變量;ρ為流體密度，kg/m3;Γ 為擴展系數;S為源項。

1.2 三維熱傳導方程

針對各向異性材料，由傳熱學基本原理可以寫出求解域內穩(wěn)態(tài)溫度場基本方程及其邊界條件為式中:T 為固體待求溫度，℃;λx、λy、λz為求解域內各種材料沿x、y以及z方向的導熱系數，W/(m·K);qV為求解域內各體熱源密度之和，W/m3;α為散熱表面的散熱系數，W/(m2·K);Tf為散熱面周圍流體的溫度，℃;Sj、Ss分別為絕熱和散熱面。

1.3 能量守恒方程

能量守恒定律是含熱交換過程的流動系統(tǒng)須滿足的基本定律，流體能量守恒方程為

式中:U為流體的速度矢量;u、v、w分別為流體的速度分量;h為流體的比焓;T為流體溫度;λ為流體的導熱系數;p為流體壓力;Sh為流體的內熱源;Φ為由于粘性作用機械能轉化為熱能的部分，稱為耗散函數。

2 求解模型確定

2.1 電機基本數據及冷卻系統(tǒng)結構

本文以1.5 MW雙饋風力發(fā)電機為例進行研究，電機基本參數:額定電壓為620 V;額定電流為1192 A;額定功率為1.5 MW;額定轉速為1 800 r/min;功率因數為0.85;頻率為50 Hz。

該發(fā)電機采用的是雙介質冷卻方式，分為內外兩路冷卻。

1)內冷卻系統(tǒng)，即風力發(fā)電機內部設置風扇，空氣由設在電機端部的風扇強迫通風實現(xiàn)封閉式冷卻，電機內通風系統(tǒng)結構如圖1所示。

2)外冷卻系統(tǒng)，發(fā)電機機殼內設有水循環(huán)系統(tǒng)，冷卻水經過與一次冷卻的熱空氣進行能量傳遞，將發(fā)電機運行因損耗而產生的熱量帶到外部。

圖1 發(fā)電機通風結構Fig.1 Ventilation structure of generator

2.2 假定條件

根據電機結構及電機內流體特性，作如下假設［9－13］:

1)電機內流場中，流體流速遠小于聲速，即馬赫數(Ma數)很小，故把流體作為不可壓縮流體處理。

2)由于電機中流體的雷諾數很大(Re＞2 300)，屬于紊流，因此采用紊流模型對電機內流場進行求解。

3)由于電機結構周向的對稱性，認為定、轉子槽中心面和定、轉子齒中心面是絕熱的。

4)本電機轉子采用隱極結構，轉子表面光滑，不考慮電機旋轉。

2.3 求解模型的確定

圖2 求解域物理模型Fig.2 Physical model of solved region

2.3.1 求解域模型

考慮電機在圓周上幾何結構周期性分布以及流體流動和傳熱特點，并依據基本假設，取電機周向1/8區(qū)域為研究對象，求解域物理模型如圖2所示，發(fā)電機各部位的名稱在圖中均已標出，圖中x軸表示定子周向方向，y軸表示定子徑向方向，z軸表示定子軸向方向。

2.3.2 網格生成

對模型進行網格劃分，采用四面體和六面體混合網格，檢查網格質量，輸出模型數據文件，求解域剖分如圖3所示。

圖3 求解域剖分圖Fig.3 Subdivision map of solved region

2.3.3 邊界條件的設置

流體場及溫度場耦合求解的邊界條件為:

1)環(huán)境溫度為27℃(300 K)，入口1為空氣入口，入口2為水入口，兩入口均給定為速度入口邊界，根據給定的質量流量數值，可求得入口速度。

2)出口1為空氣出口，出口2為水出口，均給定為壓力出口邊界條件，壓力均為1個標準大氣壓。

3)發(fā)電機機殼外圓及轉子內圓為散熱面，給定散熱系數和周圍的空氣溫度:散熱系數參考一般空冷電機設計經驗給出［18］，求解域兩側徑向截面，設為周期性邊界，其余與固體的接觸的面均視為無滑移邊界。

4)電機熱源密度分布(W/m3):定子股線為226 416;定子鐵心為38 451;轉子股線為313 747;轉子鐵心80 726。

3 內冷卻系統(tǒng)流變特性分析

由于發(fā)電機內熱性能的分布特性完全由冷卻介質的分布性能決定，首先對發(fā)電機內流體場進行分析。圖4給出了該發(fā)電機三維求解域內空氣速度分布。

通過對圖4流體速度分布以及計算數值結果分析，可知:

空氣入口截面較大，入流速度較低，冷空氣繞過非風扇端的端部繞組進入轉子通風孔和氣隙中;在電機中空氣沿軸向流過轉子通風孔和氣隙，由于氣隙入口狹長表面積較小，空氣流動入口風阻大，能量損失嚴重，空氣在氣隙中最大速度為21.22 m/s，風量占空氣總流量的10.43%;空氣在流入轉子通風孔時空氣的流動阻力系數相對較小，空氣入口面積大，風阻較小，空氣在轉子通風孔的最大速度為49.85 m/s，風量占空氣總流量的89.57%;當空氣到達風扇端區(qū)域，空間變大，空氣的速度降低，風速沿軸向截面的位置不同而不同，最高速度為43.26 m/s;在電機機殼風筋之間的風道中平均速度為32.52 m/s，空氣流經一個循環(huán)，從求解域出口流出。

圖4 空氣速度分布圖Fig.4 Velocity distribution of air

4 電機溫度場特性分析

4.1 電機溫度場分布

圖5為電機額定運行時溫升分布圖，圖6為定、轉子槽中心截面和齒中心截面溫升分布云圖。

圖5 電機溫升分布圖Fig.5 Temperature distribution of generator

圖6 電機徑向不同截面溫升分布云圖Fig.6 Temperature distribution of generator at different section

結合圖5、圖6可以得出，電機內最大溫升分布在轉子鐵心段內的轉子股線上，最大值為74.3 K;定子股線在風扇端溫升較高，鐵心段內溫升較低，并沿軸向產生一定的變化規(guī)律，定子股線的最高溫升為64.6 K;轉子股線的熱密度遠大于轉子鐵心的熱密度，兩者之間緊密接觸，轉子股線以傳導的方式將熱量傳導給轉子鐵心，轉子齒部溫升較大，最大值為62.0 K;定子鐵心傳熱條件好，能夠及時的把定子股線和鐵心產生的能量傳導給機殼，因而溫升較低，最高溫升為26.5 K。

4.2 定、轉子股線沿軸向長度溫升分布

為詳細分析股線溫升的分布規(guī)律，分別給出定、轉子股線沿軸向長度的分布，如圖7和圖8所示。

圖7 定子股線沿軸向長度溫升分布Fig.7 Temperature distribution of stator along with the axial length

圖8 轉子股線沿軸向長度溫升分布Fig.8 Temperature distribution of generator rotor along with the axial length

定子股線沿軸向長度溫升變化規(guī)律為:定子上層股線溫升高于下層股線溫升;最高溫升出現(xiàn)在風扇端，距離鐵心一定距離處，最大值為64.6 K;最低溫升分布在鐵心段的進風端。在定子鐵心段中股線溫升沿軸向長度變化趨勢基本一致，呈線性分布，由進風端至出風端溫升變化逐漸增大，相差約2～5 K，在鐵心段出風端溫升略有上升。在端部氣腔中股線溫升急劇上升并達到最大值，下層股線最大值低于上層股線約2.5 K，最高溫升處至股線端部，股線沿軸向長度溫升急劇下降，且上層股線的下降速度比下層股線快。

從圖8中可以看出轉子股線最高溫升位于0.55 m處的轉子鐵心段內，最低溫升在電機轉子鐵心段的進風端，在鐵心段內轉子股線溫升變化較快，相差約25 K。

4.3 電機齒部溫升分布

在研究股線沿軸向長度溫升分布的同時有必要對鐵心沿軸向長度溫升分布進行研究，圖9為定子鐵心徑向不同高度溫升分布，圖10為轉子齒徑向不同高度溫升分布。

圖9 定子鐵心徑向不同高度溫升分布Fig.9 Radial Temperature distribution of generator stator core at different height

圖10 轉子齒徑向不同高度溫升分布Fig.10 Radial temperature distribution of generator rotor at different height

圖9給出了定子鐵心內軛背處、距軛背1/2軛高處、齒根處、距齒頂1/2高處及齒頂處的溫升分布情況。由圖9可知，在定子鐵心徑向不同位置溫升相差很大，這是由于定子中產生的熱量大部分向軛背部傳導，硅鋼片是熱的良導體，越靠近軛背越有利于導熱;不同位置沿軸向長度溫升呈上升趨勢，只是徑向不同位置溫升變化趨勢不同，軛背處、距軛背部1/2軛高處、齒根處、距齒頂1/2高處變化趨勢一致，沿軸向長度逐漸增大;氣隙中空氣在流通過程中要帶走定子產生的一小部分熱量，空氣在氣隙中流動的過程溫升變化較快，對齒頂的散熱有影響，齒頂處溫升沿軸向長度溫升變化較大。

圖10中分別為轉子齒頂處、齒1/4處、齒1/2處、齒3/4處及齒根處沿軸向長度溫升分布，從圖中可以看出轉子齒不同高度沿軸向長度變化一致，最高溫升所在同一軸向截面上，最高約為62K;轉子齒部沿軸向不同截面最高溫升約在距齒頂1/4處;雖然齒頂距離轉子通風孔最遠，向轉子通風孔傳熱最為困難，齒頂處溫升卻不是太高，最大值為58.3K，是由于齒頂處和槽楔緊密接觸，槽楔雖然導熱能力差，但其不產生熱量，齒頂的一部分熱量向槽楔傳導，又由于氣隙中流動的空氣也能帶走一部分熱量;齒根處距離轉子通風孔最近，散熱條件好，能夠及時將自身產生的熱量和齒部傳導過來的熱量傳導給轉子通風孔中的空氣，在轉子齒軸向截面中溫升最低。

4.4 轉子股線最高溫升所在截面溫升分布

因電機轉子散熱條件較惡劣，電機轉子溫升較高，有必要對電機最高溫升處溫升分布進行分析。由于電機轉子各個相同位置溫升分布規(guī)律一致，取電機轉子任一股線最熱面為研究對象，其溫升分布如圖11所示。

圖11 轉子股線軸向最熱面溫升分布圖Fig.11 Temperature distribution of axial hottest surface of generator rotor

由圖11可以看出，轉子槽內溫升分布隨槽中心線對稱分布;最高溫升位于股線中心線上，沿銅塊周圍絕緣層溫升下降明顯，槽底和槽楔頂端散熱條件較好，溫升略低;上層股線溫升高于下層股線，上下層股線之間相差3.2 K。

4.5 計算結果與實驗值比較

發(fā)電機測溫元件埋置在風扇端，距離發(fā)電機端部一定距離的定子繞組層間絕緣處，6個測溫元件沿圓周方向均勻分布。表1給出了不同位置處的溫度測量結果。

表1 溫度測量結果Table 1 Measured data of temperature

根據表1的測量結果以及電機的實驗環(huán)境溫度，可以求出相應測量點的溫升數值，進而獲得電機的平均溫升，如表1所示。

結合溫度場求解域設定，可知求解域內有10處繞組層間絕緣位置處的溫升可與實驗值進行比較，溫升數值計算結果及其平均值與實驗值的對比如表2所示。

表2 計算值與實驗值比較Table 2 Comparison of calculated value with measured data

通過對表2中的數值進行對比分析以及對溫升計算結果的研究，可知:不同繞組層間絕緣在測量位置處的溫升計算結果差值較小;不同位置處的計算值均比實驗平均值低，最大相差5.34 K，最小誤差為0.17 K;溫升計算平均值比實驗平均值低2.63 K，誤差為4.24%，計算結果與實驗結果相吻合。

以上結果可能由于電機端部流體流動紊亂以及將電機靜止處理而轉子通風孔并非完全對稱的緣故造成。

5 結論

本文通過對1.5 MW雙饋發(fā)電機三維溫度場的計算以及對流體流動特性的分析，可以得出以下結論:

1)計算結果與實驗數值相比較誤差為4.24%，能夠滿足工程需要，驗證了結果的準確性，說明文中所采用有限體積法計算溫度場和流體場的方法可行，為發(fā)電機的運行以及大容量發(fā)電機設計有指導義。

2)從本文的分析和計算，電機轉子熱量要經過兩次傳導才能將熱量帶到電機外部，發(fā)電機轉子在電機中的溫升較高。

3)得出電機各部分溫升分布及電機內流體分布規(guī)律，為測溫元件的布置做理論指導。

4)在實際運行過程中應加強對端部股線溫度監(jiān)測，此部位由于絕緣材料導熱系數小，空氣流體繞流復雜，散熱條件最為惡劣，溫升較高，是電機溫度監(jiān)測的重點區(qū)域。

［1］ 侯俠，王靜.21世紀的綠色能源［J］.內蒙古石油化工，2006(12):52－54.

HOU Xia，WANG Jing.Green energy of 21 centuries［J］.Inner Mongolia Petrochemical Industry，2006(12):52 －54.

［2］ 程明，張運乾，張建忠.風力發(fā)電機發(fā)展現(xiàn)狀及研究進展［J］.電力科學與技術學報，2009，24(3):2－9.

CHENG Ming，ZHANG Yunqian，ZHANG Jianzhong.Development and research progress of wind power generator［J］.Journal of Electrical Power Science and Technology，2009，24(3):2 －9.

［3］ 尹超.淺談我國風力發(fā)電的現(xiàn)狀和前景［J］.山東電力技術，2009(6):14－15.

YIN Chao.Brief discussion on the present situation and prospects of wind power generation in china［J］.Shandong Electric Power，2009(6):14－15.

［4］ 施躍文，高輝，陳鐘.特大型風力發(fā)電機組技術現(xiàn)狀與發(fā)展趨勢［J］.神華科技，2009，7(2):34 －38.

SHI Yuewen，GAO Hui，CHEN Zhong.Technological development trends of super-large wind generating set［J］.Shenhua Science and Technology，2009，7(2):34 －38.

［5］ 李和明，李俊卿.電機中溫度計算方法及其應用綜述［J］.華北電力大學學報，2005，32(1):1 －5.

LI Heming，LI Junqing.Summarization on the temperature computation and application in electric machines［J］.Journal of North China Electric Power University，2005，32(1):1 －5.

［6］ 魯滌強，黃學良，胡敏強.汽輪發(fā)電機端部三維溫度場的有限元計算［J］.中國電機工程學報，2001，21(3):82 －85.

LU Diqiang，HUANG Xueliang，HU Minqiang.Using finite element method to calculate 3d thermal distribution in the end region of turbo generator［J］.Proceedings of the CSEE，2001，21(3):82 －85.

［7］ BOGLIETTI A，CAVAGNINO A，STATON D，et al.Evolution and modern approaches for thermal analysis of electrical machines［J］.IEEE Transactions on Industrial Electronics，2009，56(3):871－882.

［8］ 曹君慈，李偉力，程樹康.復合籠條轉子感應電動機溫度場計算及相關性分析［J］.中國電機工程學報，2008，28(30):96－102.CAO Junci，LI Weili，CHENG Shukang.Temperature field calculation and associated factor analysis of induction motor with compound cage rotor［J］.Proceedings of the CSEE，2008，28(30):96－102.

［9］ 溫志偉，顧國彪，王海峰.浸潤式與強迫內冷結合的蒸發(fā)冷卻汽輪發(fā)電機定子三維溫度場計算［J］.中國電機工程學報，2006，26(23):133 －138.

WEN Zhiwei，GU Guobiao，WANG haifeng.Calculation of 3D thermal field in the stator of turbo-generator with immersion evaporative-cooling system and forced inner-cooling［J］.Proceedings of the CSEE，2006，26(23):133 －138.

［10］ 霍菲，張曉晨，程鵬.屏蔽電機定子溫度場的數值計算與分析［J］.電機與控制應用，2006，33(4):7 －12.

HUO Fei，ZHANG Xiaochen，CHENG Peng.Numerical calculation and analysis of temperature field of canned motor stator［J］.Electric Machines ＆ Control Application，2006，33(4):7 －12.

［11］ 劉長紅，姚若萍.自循環(huán)蒸發(fā)冷卻電機定子鐵心與繞組間的熱量傳遞［J］.中國電機工程學報，2008，28(11):107 －112.

LIU Changhong，YAO Ruoping.Heat transfer between stator core and windings in hydro-generators with self pump circulated evaporation cooling system ［J］.Proceedings of the CSEE，2008，28(11):107－112.

［12］ 路義萍，陳朋飛，李俊亭，等.某新型空冷汽輪發(fā)電機轉子通風方式的流場分析［J］.中國電機工程學報，2010，30(6):63－68.

LU Yiping，CHEN Pengfei，LI Junting，et al.Flow field analysis of new type ventilation method in one air-cooled turbo-generator rotor［J］.Proceedings of the CSEE，2010，30(6):65 －67.

［13］ 李偉力，楊雪峰，顧德寶.空冷汽輪發(fā)電機冷卻氣流風量對定子內流體的影響［J］.中國電機工程學報，2009，30(21):53－61.

LI Weili，YANG Xuefeng，GU Debao.Influence of air current flow change on fluid flow and heat transfer of air-cooled turbo-generator with multipath ventilation［J］.Proceedings of the CSEE，2009，30(21):53 －61.

［14］ 章梓雄，董曾南.粘性流體力學［M］.北京:清華大學出版社，1998:255－279.

［15］ 鮑里先科 A N.電機中的空氣動力學與熱傳遞［M］.北京:機械工業(yè)出版社，1985:121－183.

［16］ 陶文銓.數值傳熱學［M］.西安:西安交通大學出版社，2001:347－353.

［17］ 周俊杰，徐國權，張華俊.FLUENT工程技術與實例分析［M］.北京:中國水利水電出版社，2010:14－17.

［18］ 丁舜年.大型電機的發(fā)熱與冷卻［M］.北京:科學出版社，1992.