李偉力， 袁世鵬， 霍菲陽， 張奕黃

(1.北京交通大學(xué)電氣工程學(xué)院，北京 100044;2.哈爾濱理工大學(xué)電氣與電子工程學(xué)院，黑龍江 哈爾濱 150080)

0 引言

風(fēng)能作為清潔的可再生能源為人們所關(guān)注，風(fēng)力發(fā)電技術(shù)也逐步成熟。隨著對風(fēng)力發(fā)電電機性能和成本的雙重要求提高，電機的功率密度越來越大，電磁負荷和熱負荷也越來越高，故溫升也成為電機設(shè)計中需要考慮的問題［1－7］。特別是對于大功率永磁同步發(fā)電機，溫升過高會使永磁體的磁通密度下降，影響電機性能的穩(wěn)定，甚至?xí)?dǎo)致永磁體的不可逆退磁，故電機溫度場計算尤為重要。

對于流體、傳熱等問題的分析，有限體積法有其獨特的優(yōu)勢。有限體積法得出的離散方程，要求因變量的積分守恒對任意一組控制體積都要滿足。

所以在電機溫升計算中，有限體積法也漸漸為研究者所重視［8－12］。

在傳統(tǒng)的研究方法中，大都把電機外環(huán)境的流場風(fēng)速視為一個常數(shù)，而相對應(yīng)的表面散熱系數(shù)也視為常數(shù)，但沒有將流動和傳熱過程進行耦合計算，給計算帶來很大的偏差［8－9］。

本文在流體力學(xué)理論、傳熱學(xué)理論和有限體積法的基礎(chǔ)上，首先，確定發(fā)電機熱源及等效導(dǎo)熱系數(shù)的計算方法，以及電機表面自然風(fēng)的紊流模型。其次，計算了一臺表面風(fēng)冷式小型風(fēng)力發(fā)電機的三維溫度場，并通過與實驗結(jié)果的對比，驗證了該電機溫度場計算模型的合理性。最后，利用相同方法對一臺1.5 MW永磁同步風(fēng)力發(fā)電機的三維溫度場進行計算與分析，在此基礎(chǔ)上，根據(jù)大功率永磁風(fēng)力發(fā)電機的實際特點，仿真計算了在不同環(huán)境溫度下以及不同風(fēng)速狀況下發(fā)電機的溫度分布。

1 電機熱性能參數(shù)的確定

對于全封閉表面風(fēng)冷式結(jié)構(gòu)電機，其內(nèi)部無通風(fēng)冷卻系統(tǒng)，定子繞組內(nèi)熱量散熱主要是經(jīng)定子鐵心由機殼傳遞給周圍的空氣和由鐵心向氣隙內(nèi)放熱。為了方便分析電機傳熱過程，可對電機內(nèi)部條件作相應(yīng)的假設(shè)和歸算，這對建模和有限元分析具有很大的幫助。

1.1 熱源的確定

電機運行過程中，必然要產(chǎn)生相應(yīng)的損耗，這些損耗轉(zhuǎn)換成熱，通過電機各部件與周圍環(huán)境進行熱交換，最終達到熱平衡。準確計算電機各部分損耗是準確計算電機溫度場的前提。

半直驅(qū)風(fēng)力發(fā)電機以同步轉(zhuǎn)速旋轉(zhuǎn)，主磁場在轉(zhuǎn)子鐵心的交變頻率很低，轉(zhuǎn)子鐵耗非常小;且由于電機轉(zhuǎn)速很低、頻率低，無集電環(huán)、無通風(fēng)系統(tǒng)、無阻尼繞組，故機械損耗與雜散損耗非常小，相對于銅耗與基本鐵耗近似可以忽略。

通過設(shè)計值或?qū)嶒炛荡_定定子繞組銅耗。利用場路耦合方法計算定子各部位鐵耗值:由電流幅值和內(nèi)功率角，對電機內(nèi)的電磁場進行數(shù)值計算;通過改變內(nèi)功率因數(shù)角，將計算結(jié)果代入電機負載向量關(guān)系，以電壓、功率、功率因數(shù)作為收斂目標進行迭代計算，即可確定電機的不同狀態(tài)下的單元內(nèi)磁密。定子鐵心分成齒部與軛部，提取每個單元磁密，利用公式1［13］確定定子每個單元的鐵耗，求和計算出定子總鐵耗。

式中:n為剖分的單元數(shù);p(W)為定子鐵耗值;pe(W)為每個單元的鐵耗值;Be(T)為定子鐵心內(nèi)單元的磁通密度;q為B=1 T、f=50 Hz時，硅鋼片的單位損耗;Ge(kg)為單元硅鋼片的重量;ka是鐵損增加系數(shù)，定子齒部取1.5，定子軛部取1.3。

1.2 等效導(dǎo)熱系數(shù)的確定

對于繞組內(nèi)的銅導(dǎo)體與絕緣材料的復(fù)雜分布，為簡化分析，對定子槽做如下假設(shè):

1)浸漬狀況良好，浸漬漆填充均勻;

2)銅線的絕緣漆分布均勻;

3)槽絕緣和鐵心緊密結(jié)合在一起;

4)槽內(nèi)各個導(dǎo)線的溫差忽略不計。

根據(jù)以上假設(shè)，可將定子槽內(nèi)的繞組、絕緣與氣體等各種材料效成一個與槽等體積的繞組，其導(dǎo)熱系數(shù)的計算可根據(jù)公式2確定［14－15］。

式中:λeq(W/m·K)為槽內(nèi)材料的等效導(dǎo)熱系數(shù);δi(i=1，2，3，…，n)(m)為槽內(nèi)各種材料的等效厚度;λi(W/m·K)為各種材料的導(dǎo)熱系數(shù)。

計算定轉(zhuǎn)子全域溫度場時，需要對氣隙進行特殊處理。轉(zhuǎn)子的旋轉(zhuǎn)帶動氣隙中空氣的流動，使得定子與氣隙之間以及轉(zhuǎn)子與氣隙之間主要以對流方式換熱。為簡化計算分析，引入有效導(dǎo)熱系數(shù)λeff(W/m·K)，它是用靜止流體的導(dǎo)熱系數(shù)來描述氣隙中流動空氣的熱交換能力。這樣可以用導(dǎo)熱方式換熱的效果等價氣隙中對流方式換熱的效果。

由于氣隙中的空氣流動為紊流，這種狀態(tài)下，氣隙的有效導(dǎo)熱系數(shù)可用下式計算［15］，即

1.3 耦合模型中流體的數(shù)學(xué)描述

由流體傳熱理論和采用有限體積法建立電機表面通風(fēng)冷卻過程的數(shù)學(xué)模型，流體流動過程通過以下方程組［11－12］進行描述和控制:

假設(shè)流體不可壓縮，故流體密度為常值，流體質(zhì)量守恒方程為

式中:ρ(kg/m3)為流體的密度;u(m/s)為流體的速度矢量;t(s)為時間。

流體動量守恒方程為

式中:F(N)為流體受到的質(zhì)量力;p(N)為流體壓力;μ為流體的粘性系數(shù)。

自然界的大氣流動以湍流流動為最常見。采用標準k－ε模型能夠很好的模擬一般的湍流［12］。標準模型使用湍流動能k和擴散因子ε的通用控制方程，即

式中:φ、V為通用變量;Γ為擴散系數(shù);S為源項。

2 樣機的三維溫度場耦合計算與實驗對比

2.1 實驗樣機的計算模型

以一臺小型離網(wǎng)式永磁同步風(fēng)力發(fā)電機為樣機，建立該電機的三維溫度場計算模型，描述了電機與大氣流動的耦合傳熱過程。為了簡化模型與計算難度，假設(shè)機殼表面通風(fēng)條件一致;電機定子繞組端部等效成與定子槽內(nèi)繞組等截面積的平行體;電機外圍的大氣為一層足夠大的氣層流域，邊界為絕熱面。其三維溫度場的計算模型如圖1所示，其中，面S1是自然風(fēng)入風(fēng)側(cè)，面S2是出風(fēng)側(cè)，面S3即機殼外表面是電機的散熱面。

實驗電機主要參數(shù)為:額定容量400 W;額定風(fēng)速8 m/s;定子外徑0.15 m;轉(zhuǎn)子外徑0.0922 m;軸向長度0.06 m;永磁體厚度4.5 mm;定子槽數(shù)33;三相，極數(shù)10。

流體模型在求解溫度中，除要滿足方程(4)～(6)外，還須滿足如下能量守恒方程

式中:U為流體的速度矢量;u、v、w為流體的速度分量;h為流體的比焓;T為流體溫度;λ為流體的導(dǎo)熱系數(shù);Sh為流體的內(nèi)熱源;Φ為由于粘性作用機械能轉(zhuǎn)化為熱能的部分，稱為耗散函數(shù)。

圖1 實驗電機三維溫度場的計算模型Fig.1 The 3-D temperature field calculated model of the prototype

根據(jù)流動和傳熱特點，給出模型求解邊界條件為:

1)電機自然風(fēng)冷卻，面S1定義為風(fēng)速入口邊界條件，根據(jù)外界條件設(shè)定入口風(fēng)速的大小與溫度;面S2定義為出流邊界條件。

2)模型中的定子齒部、定子軛部以及定子繞組均給定為熱源體。

在熱源、固體的導(dǎo)熱系數(shù)、邊界條件和流體狀態(tài)已知條件下，由流體傳熱理論來求解電機內(nèi)的溫度場。

2.2 計算結(jié)果與實驗值對比

為了印證計算方法的可行性，本文利用電機的實驗數(shù)據(jù)與計算的結(jié)果做對比驗證。實驗中機殼表面的風(fēng)速由EY3－2A電子微風(fēng)儀測定;采用熱敏電阻測量電機內(nèi)各部分的溫度;用Raytek紅外測溫儀測量機殼的溫度。

本文在不同負載情況和不同的風(fēng)速情況下，進行溫升實驗，表1顯示其中3種工況下實驗數(shù)據(jù)［16］。根據(jù)測量結(jié)果，進行銅耗和定子各部位鐵耗的計算。將各種損耗轉(zhuǎn)換成熱流密度施加于電機溫度場模型中，依照上述的計算方法，進而可計算電機內(nèi)溫度分布。

表1 3種不同工況實驗數(shù)據(jù)Table 1 Test data at three different states

圖2是第3種工況下，電機三維溫度場沿中軸線的橫切圖與縱切圖。表2是電機3個部位溫度的測量值與計算值對照。

圖2 樣機溫度場計算結(jié)果Fig.2 The calculated result of the prototype

表2 計算值與實驗值對比Table 2 Comparison in the calculated and measured results

3 1.5 MW永磁半直驅(qū)風(fēng)力發(fā)電機三維溫度場計算與分析

3.1 建立耦合計算模型

由表2可以看出，計算結(jié)果與實測值比較接近，說明計算方法可行。根據(jù)相同計算方法，本文對一臺采取表面自然風(fēng)冷卻的1.5 MW永磁半直驅(qū)風(fēng)力發(fā)電機進行三維溫度場計算與分析。該風(fēng)力發(fā)電機的主要參數(shù)為:額定容量1500 kW;額定轉(zhuǎn)速150 r/min;額定電壓690 V;定子外徑2.2 m;額定風(fēng)速15 m/s;軸向長度0.62 m;氣隙長度7 mm;定子槽型11.5 mm×64 mm;永磁體厚度64 mm。

根據(jù)電機對稱關(guān)系理論，取電機的一個單元進行三維溫度場計算，其計算模型如圖3所示，其中，面S1是自然風(fēng)入風(fēng)側(cè)，面S2是出風(fēng)側(cè)，面S3即機殼外表面是電機的散熱面。

圖3 1.5 MW PMSG三維溫度場的計算模型Fig.3 The 3-D temperature field calculated model of the 1.5 MW PMSG

3.2 計算結(jié)果與分析

在額定風(fēng)速下，風(fēng)力發(fā)電機功率達到額定值，也即最大值。圖4是在環(huán)境溫度為20℃，額定風(fēng)速15 m/s時，半直驅(qū)風(fēng)力發(fā)電機三維溫度場沿中軸線的橫切面圖與縱切面圖。

圖4 1.5 MW PMSG計算結(jié)果Fig.4 The calculated result of the 1.5 MW PMSG

由于定子繞組和定子鐵心是電機的主要熱源，轉(zhuǎn)子側(cè)損耗小，電機外表面是主要的散熱途徑，所以最高溫度出現(xiàn)在上層定子繞組。考慮到端蓋與軸具有一定的散熱作用，轉(zhuǎn)子鐵心由外向內(nèi)溫度遞減。

圖5是定子繞組沿軸向的最高溫度分布圖。由于封閉式風(fēng)冷的電機，內(nèi)部無散熱風(fēng)溝，繞組端部散熱條件差，溫度較繞組中部要高。由于入風(fēng)側(cè)的散熱條件好于出風(fēng)側(cè)的，所以上風(fēng)側(cè)的繞組端部溫度會略低于下風(fēng)側(cè)的端部繞組。

圖5 繞組軸向溫度分布圖Fig.5 The axial temperature distribution of windings

兆瓦級風(fēng)力發(fā)電機的工作條件往往在野外或海上，四季溫差和晝夜溫差都比較大。一般考慮環(huán)境溫度在－20～40℃，圖6是發(fā)電機工作在額定狀態(tài)下，不同環(huán)境溫度時，電機額定狀態(tài)下各部的溫度變化圖。

圖6 不同環(huán)境溫度下電機的溫度Fig.6 The temperature of generator at different ambient temperature

風(fēng)機的切入風(fēng)速為3 m/s，切出風(fēng)速為25 m/s，在低于額定轉(zhuǎn)速時，假定風(fēng)機運行在漿距角為0、最佳葉尖速比的狀態(tài)下;高于額定風(fēng)速時，變漿距機構(gòu)發(fā)揮作用，使電機轉(zhuǎn)速限制在額定值附近。故給定其漸變風(fēng)風(fēng)速與電機轉(zhuǎn)速的關(guān)系［17］如圖7所示。在額定轉(zhuǎn)速以下，電機輸出功率由式(8)［18］的比例關(guān)系確定。在額定風(fēng)速以上時，電機以額定功率輸出，表3是不同風(fēng)速下的損耗分布。由此計算在20℃環(huán)境溫度、不同風(fēng)速狀態(tài)下電機各部溫度，如圖8所示。

式中:Pm為風(fēng)力機所捕獲的能量;Cp為功率利用系數(shù);ρ為空氣密度;A為風(fēng)葉掃過的面積;υ0為風(fēng)速。

由于風(fēng)機吸收的功率與風(fēng)速的三次方成正比，所以在低于額定轉(zhuǎn)速時，電機溫度隨風(fēng)速的增大而快速增高;當(dāng)超過額定轉(zhuǎn)速后，電機恒功率輸出，機殼表面風(fēng)速增加，使其溫度隨著風(fēng)速的增加而緩慢降低。

圖7 給定風(fēng)速與電機轉(zhuǎn)速關(guān)系Fig.7 The given relation between wind speed and generator rotary speed

表3 不同風(fēng)速下電機損耗Table 3 The loss distribution at different wind speed

圖8 不同風(fēng)速下電機溫度Fig.8 The temperature of generator at different wind speed

4 結(jié)論

通過對實驗樣機溫度場的計算結(jié)果與實驗值對比，以及對1.5MW半直驅(qū)式永磁同步風(fēng)力發(fā)電機的三維溫度場計算，可以得出如下結(jié)論:

1)通過對小型實驗樣機三維溫度場的計算與實驗結(jié)果對比，驗證了計算的準確性，說明文中所采取熱源、等效散熱系數(shù)的計算方法以及流體傳熱計算模型的可行性。為表面自然風(fēng)冷卻電機的溫度場計算提供了參考。

2)通過三維溫度場的計算可以得出自然風(fēng)冷卻電機沿軸、徑向溫度的分布規(guī)律。在徑向上最高溫度出現(xiàn)在定子繞組上層。在軸向上，由于端部散熱能力較差，導(dǎo)致繞組端部溫度要高于繞組中部，并且電機入風(fēng)側(cè)的溫度要低略于出風(fēng)側(cè)的。

3)對于變速恒頻的大型永磁同步風(fēng)力發(fā)電機，在額定風(fēng)速附近電機溫度達到最高。低于額定風(fēng)速時，電機溫升隨風(fēng)速的增加而快速增大;高于額定風(fēng)速后，電機溫升隨風(fēng)速的增加而緩慢降低。

［1］張岳，王鳳翔.直驅(qū)式永磁同步風(fēng)力發(fā)電機性能研究［J］.電機與控制學(xué)報，2009，13(1):78 －82.ZHANG Yue，WANG Fengxiang.Performance research on permanent magnet synchronous generator directly driven by wind turbine［J］.Electric Machines and Control，2009，13(1):78 －82.

［2］姚若萍，饒方權(quán).蒸發(fā)冷卻水輪發(fā)電機定子溫度場研究［J］.中國電機工程學(xué)報，2003，23(6):87 －90.YAO Ruoping，RAO Fangquan.Analysis of 3D thermal field in the stator of large hydro-generator with evaporation-cooling system［J］.Proceedings of the CSEE，2003，23(6):87 －90.

［3］溫志偉，顧國彪.實心磁極同步電動機轉(zhuǎn)子溫度場計算［J］.大電機技術(shù)，2005(2):1 －5.WEN Zhiwei，GU Guobiao.Calculation of rotor temperature field for solid pole synchronous motor［J］.Large Electric Machine and Hydraulic Turbine，2005(2):1 －5.

［4］黃學(xué)良，胡敏強，周鶚.電機三維溫度場新的有限元模型［J］.中國電機工程學(xué)報，1998，18(2):78 －82.HUANG Xueliang，HU Minqiang，ZHOU E.Novel FEM method to compute 3D thermal distribution of electrical machine［J］.Proceedings of the CSEE，1998，18(2):78 －82.

［5］BOGLIETTI A，CAVAGNINO A，LAZZARI M，et al.A simplified thermal model for variable-speed self-cooled induction motor［J］.IEEE Transactions on Industry Applications，2003，39(4):945－952.

［6］魏永田，孟大偉.電機內(nèi)熱交換［M］.北京:機械工業(yè)出版社，1998:251－344.

［7］GUO Jianzhong，GUO Hui，HOU Zhe.Rotor temperature monitoring technology of direct-drive permanent magnet wind turbine［C］//2009 International Conference on Electrical Machines and Systems，November 15 －18，2009，Tokyo，Japan.2009:1 －4.

［8］路義萍，馬賢好，李偉力.空冷汽輪發(fā)電機轉(zhuǎn)子風(fēng)道結(jié)構(gòu)對傳熱的影響［J］.中國電機工程學(xué)報，2007，27(17):61 －66.LU Yiping，MA Xianhao，LI Weili.The effect of geometry factors on heat transfer of rotor’s cooling duct of large turbo generator with air-coolant［J］.Proceedings of the CSEE，2007，7(17):61 －66.

［9］霍菲陽，李勇，李偉力，等.大型空冷汽輪發(fā)電機定子通風(fēng)結(jié)構(gòu)優(yōu)化方案的計算與分析［J］.中國電機工程學(xué)報，2010，30(6):69－75.HUO Feiyang，LI Yong，LI Weili，et al.Calculation and analysis on stator ventilation structure of different optimum proposal in aircooled turbogenerator［J］.Proceedings of the CSEE，2010，30(6):69－75.

［10］孟大偉，劉宇蕾，張慶軍，等.潛油電機整體三維溫度場耦合計算與分析［J］.電機與控制學(xué)報，2010，14(1):52 －55.MENG Dawei，LIU Yulei，ZHANG Qingjun，et al.Calculation and analysis of 3D temperature of the submersible motor based on coupled field［J］.Electric Machines and Control，2010，14(1):52 －55.

［11］張梓雄，董曾南.粘性流體力學(xué)［M］.北京:清華大學(xué)出版社，1998:58－101.

［12］吳德銘，郜冶.實用計算流體力學(xué)基礎(chǔ)［M］.哈爾濱:哈爾濱工程大學(xué)出版社，2006:117－266.

［13］汪耕，李希明.大型汽輪發(fā)電機設(shè)計、制造與運行［M］.上海:上?？茖W(xué)技術(shù)出版社，2000:148－151.

［14］XYPTERAS J，HATZIATHANASSIOU V.Thermal analysis of an electrical machine taking into account the iron losses and the deep-bar effect［J］.IEEE Transactions on Energy Conversion，1999，14(4):996 －1003.

［15］HATZIATHANASSIOU V，XYPTERAS J，ARCHONTOULAKIS G.Electrical-thermal coupled calculation of an asynchronous machine［J］.Archiv Fur Electrotechnik，1994，77(2):117 －122.

［16］李偉力，仝世偉，程鵬.離網(wǎng)型永磁同步風(fēng)力發(fā)電機電磁場和溫度場數(shù)值計算與分析［J］.中國電機工程學(xué)報，2010，30(30):107－113.LI Weili，TONG Shiwei，CHENG Peng.Calculation and analysis of electromagnetic and temperature fields in off-grid type permanent magnet synchronous generator［J］.Proceedings of the CSEE，2010，30(30):107 －113.

［17］陳昆明，湯天浩，陳新紅，等.永磁半直驅(qū)風(fēng)力機控制策略仿真［J］.上海海事大學(xué)學(xué)報，2008，29(4):39 －44.CHEN Kunming，TANG Tianhao，CHEN Xinhong，et al.Control strategy simulation on half-direct permanent magnet wind energy conversion system［J］.Journal of Shanghai Maritime university，2008，29(4):39 －44.

［18］王承煦，張源.風(fēng)力發(fā)電［M］.北京:中國電力出版社，2003:56－57.