溫嘉斌, 鄭軍, 于喜偉

(1.哈爾濱理工大學(xué)電氣與電子工程學(xué)院,黑龍江哈爾濱150080;2.南陽防爆電機廠,河南南陽473000)

YKK中型高壓異步電動機換熱分析與冷卻器優(yōu)化

溫嘉斌1, 鄭軍1, 于喜偉2

(1.哈爾濱理工大學(xué)電氣與電子工程學(xué)院,黑龍江哈爾濱150080;2.南陽防爆電機廠,河南南陽473000)

為了提高電機的散熱性能,改善電機的溫升條件,本文采用數(shù)值分析的方法對電機冷卻器的冷卻性能加以優(yōu)化。對于YKK結(jié)構(gòu)的封閉式交流異步電動機,因其結(jié)構(gòu)特點,溫升一直是設(shè)計時需要關(guān)注的重點問題之一。首先,以YKK450-4、800 kW中型高壓異步電機為例,建立電機冷卻器的三維物理模型和數(shù)學(xué)模型,給出假設(shè)條件和邊界條件,并進行數(shù)值分析,得到了冷卻器的流體場和溫度場是分布云圖;其次,通過改變檔風板的個數(shù)及位置、增加導(dǎo)風板并改變其位置,找到其最優(yōu)值;最后,對冷卻器進行重新建模,經(jīng)過數(shù)值分析發(fā)現(xiàn)冷卻器的性能有了明顯的改善,為以后冷卻器的優(yōu)化設(shè)計提供理論依據(jù)。

中型高壓異步電動機;冷卻器;流體場;溫度場;優(yōu)化設(shè)計

0 引 言

YKK系列中型高壓異步電動機廣泛應(yīng)用于冶金、電力、化工及水處理等場合。該系列電機具有高效、節(jié)能、起動力矩大、過載能力強、噪音低、振動小及運行可靠等特點。電機運行的可靠性是考察電機性能好壞的重要指標,而對電機運行可靠性的評判又主要體現(xiàn)在電機運行過程中的溫升問題。電機內(nèi)溫升過大,會導(dǎo)致絕緣出現(xiàn)分層、脫殼、老化等現(xiàn)象,從而使絕緣的介電性能下降引起絕緣的損壞,引發(fā)電機內(nèi)部的各種放電、短路故障,最終導(dǎo)致電機燒毀［1］。通過分析電機的實際運行狀態(tài),準確描述電機內(nèi)各部件溫度場的分布,為電機故障的實時預(yù)警和診斷提供可靠的依據(jù);并還可以與一些優(yōu)化設(shè)計理論相結(jié)合。在冷卻器的優(yōu)化設(shè)計中可以采用兩種方法。一種是做風洞實驗［2］,其可信度高但是成本較大;另一種方是利用計算機做仿真分析來對冷卻器進行優(yōu)化設(shè)計。本文采用第二種方法,這樣不但可以節(jié)約生產(chǎn)成本,還可以減少樣機的制造周期,具有重要的意義。

1 冷卻器物理模型與數(shù)學(xué)模型

YKK450-4、800 kW電動機通風結(jié)構(gòu)如圖1所示。

圖1 電機通風結(jié)構(gòu)Fig.1 Ventilation system of them otor

圖1 所示電機內(nèi)風路采用密閉循環(huán)結(jié)構(gòu),內(nèi)部安裝擋風板對流體流動方向進行約束,在電機左側(cè)安裝一個離心式風扇,通過離心式風扇產(chǎn)生足夠大的壓強迫使內(nèi)風路氣體進行循環(huán)流動。內(nèi)風路氣體依次經(jīng)過定子端部、軸向通風溝、轉(zhuǎn)子徑向通風溝、氣隙、定子徑向通風溝、內(nèi)風扇、冷卻器,最后回到定子端部形成一個封閉的循環(huán)系統(tǒng)。圖中虛線箭頭所指方向為內(nèi)風路氣體的流動方向。為了將內(nèi)風路氣體冷卻,在電機外部安裝離心式風扇,使溫度較低的大氣進人冷卻管從而達到冷卻內(nèi)風路氣體的目的。外風路冷卻氣體的流動方向為大氣、外風扇、冷卻器人風筒、冷卻管最后回到大氣。圖中實線箭頭所指方向為外風路氣體流動方向。

1.1 冷卻器區(qū)域流體場物理模型

本文電機的冷卻器有713根冷卻管,冷卻管排列示意圖如圖2所示。

圖2 冷卻管排列方式Fig.2 Cooling arrangement

冷卻器由冷卻管、人風筒、擋風板組成。本文中將對冷卻器的擋風板和人風筒進行優(yōu)化設(shè)計,將擋風板的個數(shù)由1個增加到3個,在人風筒處增加導(dǎo)風板并且改變斜板的角度,從而減小渦流的產(chǎn)生。冷卻器系統(tǒng)的整體模型如圖3所示。

圖3 冷卻器物理模型Fig.3 Physicalmodel of cooler

1.2 基本假設(shè)和邊界條件

1.2.1 基本假設(shè)

1)計算區(qū)域內(nèi)的雷諾數(shù)很大,流動屬于湍流,即采用湍流模型［7］;

2)常壓下空氣冷卻的電機,忽略流體域內(nèi)流體的浮力和重力的影響［7］;

3)流體流速遠小于聲速,即把電機內(nèi)流體作為不可壓縮流體處理［7］;

4)僅研究冷卻器內(nèi)流體場與溫度場的穩(wěn)定狀態(tài),即定常流動,因而方程中不含有時間項［8］。

1.2.2 邊界條件

1)外風路人口速度為4.7 m/s(人口處5個點速度的平均值),人口溫度為22.0℃(室溫),采用速度人口;外風路出口采用壓力出口;

2)內(nèi)風路人口速度為6.0 m/s(人口處測出3個點速度的平均值),人口溫度為80.0℃(人口處測量8個點溫度的平均值),采用速度人口;內(nèi)風路可以計算出電機外殼表面的對流換熱系數(shù)為30W/(m2·k)。

1.3 冷卻器區(qū)域流體場數(shù)學(xué)模型

流體流動所遵循的物理定律是建立流體運動基本方程組的依據(jù)。這些定律包括質(zhì)量守恒、能量守恒和動量守恒［3］三大定律。

質(zhì)量守恒方程為出口采用壓力出口;

3)所有流體和固體的交界面設(shè)置為藕合壁面,冷卻器外殼和外界空氣對流換熱,根據(jù)電機外殼內(nèi)壁的風速v(m/s)和電機外殼表面的溫度θ(K),即

動量守恒方程為

式中:p為流體微元體上的壓力;τxx、τxy和τxz為粘性應(yīng)力τ沿x、y和z方向的分量;Fx、Fy和Fz為微元體上的體積力。

因為本文設(shè)計到溫度場計算,還需要能量守恒方程和牛頓冷卻定律及表面?zhèn)鳠嵯禂?shù)。

能量守恒方程為

式中:keff=kt+k為有效導(dǎo)熱系數(shù);Jj'為組分j'的擴散系數(shù)動能耗散率;方程右邊為導(dǎo)熱項、組分擴散項和粘性耗散項;Sh是化學(xué)反應(yīng)熱和其它體積熱源。

牛頓冷卻定律表達式［4,10］為

式中:Φ為單位時間通過面積A對流熱流量;Tw為固體表面絕對溫度;T為流體平均溫度;A為傳熱表面積;q為對流換熱熱流面積;α為對流傳熱系數(shù)［W/(m2·K)］。

在Fluent中進行數(shù)值分析時,采用標準模型,其湍流方程［5］為式中:k為湍動能;ε為湍動能耗散率;μi為i方向速度;μ為動力粘度;Gk為由于平均速度梯度引起的湍動能k的產(chǎn)生項,Gb是由于浮力引起的湍動能k的產(chǎn)生項,YM代表可壓湍流中脈動擴張的貢獻,Sk和Sε為自定義的源項;μt為湍動粘度［6］,可以表示成k和ε的函數(shù),即μt=ρCμ在Fluent中取經(jīng)驗值C1ε=1.44,C2ε=1.92,C3ε=0.09;湍動能k與耗散率ε的湍流普朗特數(shù)分別為σk=1.0,σε=1.3。

2 冷卻器區(qū)域流體場溫度場耦合計算

2.1 原冷卻器區(qū)域流體場與溫度場計算

由于冷卻器有713根冷卻鋁管且整體尺寸較大,在進行網(wǎng)格剖分和數(shù)值計算時對計算機的內(nèi)存和運算有很大的要求,因為冷卻器的結(jié)構(gòu)是對稱的,所以本文選取1/9的冷卻器進行建模計算分析。整個計算區(qū)域剖分約290萬個單元,370萬個節(jié)點。因為本文要對冷卻管的表面散熱情況進行分析,所以對其采用結(jié)構(gòu)化網(wǎng)格剖分,剖分后的模型如圖4所示。

圖4 計算區(qū)域網(wǎng)格剖分Fig.4 Computational domain mesh dissection

剖分、定義邊界條件后,數(shù)值計算后得到冷卻器的溫度云圖如圖5所示、速度矢量圖如圖6所示、流體跡線圖如圖7所示和外管壁對流傳熱系數(shù)如圖8所示。

圖5 截面的溫度圖Fig.5 The temperature distribution of a section

圖6 截面的速度矢量圖Fig.6 The vector velocity of a section

圖7 截面的流體跡線圖Fig.7 Computational domain mesh dissection

圖8 鋁管內(nèi)表面對流傳熱系數(shù)圖Fig.8 A lum inum inner surface convective heat transfer coefficient

從圖5中可以看出內(nèi)風路人口溫度最高達到80.0℃,經(jīng)過與冷卻管接觸,在內(nèi)風路出口的溫度降低到52.5℃,外風路的低溫氣體經(jīng)過冷卻器后由22.0℃升高到51.0℃。但是從冷卻管溫度分布圖可以看出其溫度升高的不是很多,說明冷卻空氣沒有得到充分的利用。從圖6中看出人風筒區(qū)域的流體進人冷卻管后流速增大,這是由于流體流動的截面積突然變小造成的,流量不變,流速增大。從圖7中看出A處有明顯的渦流出現(xiàn),這將會導(dǎo)致能量損失加大并且產(chǎn)生很大的噪音［9］,但是擋風板處回流的存在有助于熱空氣與冷卻管的接觸,使內(nèi)風路溫度降低。從圖8中可以看出在內(nèi)、外風路人口處的對流傳熱系數(shù)比較大,這是因為空氣剛進人冷卻管,流速突然增大,流速增加又會使流體內(nèi)部相對運動加強,從而使對流作用加強,也就使對流換熱加強,對流換熱系數(shù)增大［10］。

為驗證數(shù)值分析合理性,對YKK450-4、800 kW電機進行試驗對比,使電機工作在額定功率下,達到穩(wěn)定運行時,測得多組數(shù)據(jù)來進行求取平均值作為實驗數(shù)據(jù)。其中外風路出口處的風速共測15個點。表1中給出的試驗數(shù)據(jù)與數(shù)值計算結(jié)果的對比圖。

表1 實驗與數(shù)值分析對比Table 1 The contrast of experience and numerical analysis

2.2 對冷卻器優(yōu)化設(shè)計計算

從原冷卻器模型中可以看到冷卻管的冷卻空氣沒有被充分的利用;經(jīng)過分析后,增加導(dǎo)風板來均衡分配每層管的流量;減少人風筒直角的使用可以明顯減少渦流的產(chǎn)生。基于此本文對冷卻器進行優(yōu)化設(shè)計。改變結(jié)構(gòu)后的冷卻器物理模型如圖9所示。將擋風板的個數(shù)由1塊變?yōu)?塊,即增加擋板1和3,將原擋板上移變?yōu)閾醢?,對各個部分的優(yōu)化設(shè)計進行分別做如下闡述。

圖9 優(yōu)化后的模型Fig.9 Themodel of optim ization

2.2.1 對導(dǎo)風板的優(yōu)化

增添一個導(dǎo)風板,并將導(dǎo)風板在原有的基礎(chǔ)上分別上移20、40、60和80 mm。分別建模進行數(shù)值分析計算,得出上移40mm時各層流量分配比較均勻。各層流量分配對比圖如圖10所示,內(nèi)風路出口溫度對比如圖11所示。從圖10中可以看出改進后冷卻管人口的流量分配更加均勻,這樣會使得管壁散熱效果更好。在圖11中能夠看出導(dǎo)風板上移40mm時內(nèi)風路出口有所降低,溫度為51.5℃。

圖10 各層流量分配對比Fig.10 Contrast flow distribution layers

圖11 內(nèi)風路出口溫度對比Fig.11 Contrast of the inner air duct outlet temperature

再將導(dǎo)風板在原有的基礎(chǔ)上分別外移10、20、30和40 mm。重新建模進行對比分析得出將導(dǎo)風板外移30mm時溫度下降教明顯,原因是改變導(dǎo)風板的位置后,從流體跡線圖進行對比分析得出將導(dǎo)風板外移30 mm時使渦流明顯減少,能量損失減少,進而溫度下降,其值為50.8℃。內(nèi)風路出對溫度對比圖如圖12所示。

圖12 導(dǎo)風板右移溫度對比圖Fig.12 Contrast of temperature of the wind deflector right

2.2.2 對擋板的優(yōu)化

優(yōu)化前的冷卻器模型只有一塊擋風板,內(nèi)風路出來的風沒有得到充分的冷卻,所以將擋風板改成3塊,如圖9所示。改變模型后進行數(shù)值計算分析,內(nèi)風路的出口溫度有明顯的降低,改進前后溫度對比圖如圖13所示。從圖中可以看到溫度降到50.2℃。但是通過數(shù)值分析后發(fā)現(xiàn)內(nèi)風路人口所需要的靜壓會提高,所以要想辦法將內(nèi)風路人口所需靜壓強減小。在圖9中將3塊擋板標號以便區(qū)分,首先擋板3的位置保持不變,改變擋板1的位置,并且使擋板2始終處于擋板1和3的中間位置。將擋板1從內(nèi)人口的右端點依次右移70、140和210mm。計算后發(fā)現(xiàn)改變擋板的位置對內(nèi)風路出口溫度幾乎沒有影響,只對內(nèi)風路人口所需要的靜壓有很大的影響。對比圖如圖14所示。從圖中看出將擋板1右移210mm時所需要的靜壓強最小為2 776.8 Pa。

圖13 增加擋板個數(shù)溫度對比圖Fig.13 Contrast of tem perature at increased the number of baffles

圖14 減靜壓對比圖Fig.14 Contrast of the less static pressure

2.2.3 對冷卻器直角形狀的優(yōu)化

流體在流動過程中,除了與管壁之間的摩擦阻力以外,由于管道形狀、尺寸變化或發(fā)生某些障礙時,流體質(zhì)點之間的相對速度、方向?qū)l(fā)生變化,碰撞將加劇,便產(chǎn)生一定的能量損失［10］。

以壓力降的形式表示為

式中:ζ表示局部阻力系數(shù);ζ'表示風阻系數(shù);Z表示風阻。

電機冷卻系統(tǒng)內(nèi),局部損失占很大比重,因此局部阻力系數(shù)ζ對電機通風是必不可少的。本文將圖中畫圈的M處的斜板用M'處的圓弧板來代替,這樣可以使ζ降低,從改進后的流體跡線圖中能夠看到M'處的渦流有明顯的減少。斜板改變前后的對比圖如圖15所示。

圖15 改變斜板形狀對比圖Fig.15 Contrast of changed the inclined p late shape

2.3 對冷卻器綜合優(yōu)化設(shè)計結(jié)果分析

經(jīng)過上文對冷卻器結(jié)構(gòu)的改變,按照各個部分的最優(yōu)值方法設(shè)計出冷卻器的物理模型,如表2所示。在Fluent數(shù)值迭代后,得到溫度云圖如圖16所示、速度矢量圖如圖17所示、流體跡線圖如圖18所示和內(nèi)管壁傳熱系數(shù)如圖19所示。

表2 整體優(yōu)化前后對比Table 2 The contrast result of holistic optim ization

圖16 改進后截面溫度圖Fig.16 Improved of tem perature distribution of a section

圖17 改進后速度矢量圖Fig.17 Improved of vector velocity of a section

圖18 改進后的流體跡線圖Fig.18 Improved of fluid traces of a section

圖19 改進后的內(nèi)管壁傳熱系數(shù)Fig.19 Improved of alum inum inner surface convective heat transfer coefficient

改進后冷卻器內(nèi)風路出口溫度為49.2℃,外風路出口溫度升高到54.9℃,內(nèi)風路人口所需要的靜壓為2 776.8 Pa。內(nèi)風路出口的溫度較改變前降低了3.3℃,外風路出口的溫度升高了3.9℃,較原結(jié)構(gòu)有明顯的改善。從速度矢量圖中可以看到改進后冷卻管中的流速大于改進前的,這是因為將人風筒中的渦流減少,從而使能量損失降低,速度增大,流速的增大會使得冷空氣快速流動帶走更多的熱量。從流體跡線圖中可以看到A'上下處的渦流明顯減少,這有助于減少噪音和能量的損失。通過管壁傳熱系數(shù)分布圖中我們可以看到增加擋板的個數(shù)使得熱空氣與冷卻管進行對流換熱的次數(shù)增加,熱空氣和冷卻管能夠更加充分接觸,平均的對流傳熱系數(shù)較改變模型前有很大的提高,從而將內(nèi)風路出口溫度降的更低。

3 結(jié) 論

1)通過數(shù)值分析計算得到冷卻器內(nèi)外風路的流體場和溫度場分布,從中可以看出人風筒渦流較大,能量損失大、噪音大。經(jīng)改進后渦流減少,能量損失減少,冷卻管中流速提高,使冷卻管能夠帶走更多的熱量,從而提高冷卻器的散熱性能;

2)改變擋風板的個數(shù)能夠明顯的改善內(nèi)風路出口的溫度,并且通過改變擋板的位置能夠減少內(nèi)風路人口所需要的靜壓;

3)通過本文分析得出,在以后設(shè)計模型時要避免直角和尖角的出現(xiàn)盡量用平滑的圓弧來代替,這樣可以減少渦流的產(chǎn)生。

［1］ 張新波,許承千,楊紅.電機內(nèi)冷卻器換熱效果的綜合計算［J］.電工技術(shù)學(xué)報,1998,8:8-9. ZHANG Xinbo,XU Chengqian,YANG Hong.Consolidation of themotor cooler of effectof heat transfer［J］.Transaction of China Electrotechnical Society,1998,8:8-9.

［2］ NAKAHAMA Takafumi,BISWAS Debasish,KAWANO Koichiro,et al.Improved cooling performance of largemoors using fans［J］. IEEE Transactions on Energe Conversion,2006,21(2):324 -331.

［3］ 王福軍.計算流體動力學(xué)分析CFD軟件原理與應(yīng)用［M］.北京:清華大學(xué)出版社,2004:7-10.

［4］ 孟大偉,溫嘉斌,魏永田,等.雙層筒結(jié)構(gòu)燈泡式水輪發(fā)電機熱交換與肋優(yōu)化計算［J］.電工技術(shù)學(xué)報.1998,13(1):18-22. MENG Dawei,WEN Jiabin,WEIYongtian,et al.Heat exchange and rib optimal calculation for bulb tubular turbine generator with layer structure［J］.Transaction of China Electrotechnical Society,1998,13(1):18-22.

［5］ 江帆,黃鵬.Fluent高級應(yīng)用與實力分析［M］.北京:清華大學(xué)出版社,2008:13-15.

［6］ 陶文銓.數(shù)值傳熱學(xué)［M］.第2版.西安:西安交通大學(xué)出版社,2001:348-350.

［7］ 路義萍,豐帆,孫明琦,等.同步電機定子與氣隙流場數(shù)值計算與分析［J］.電機與控制學(xué)報,2011,15(8):47-51. LU Yiping,FENG Fan,SUN Mingqi,et al.Numerical calculation and analysis of fluid flow field of stator and gap of a synchronousmachine［J］.Electric Machines and Control,2011,15(8): 47-51.

［8］ 李偉力,楊雪峰,顧德寶,等.多風路空冷汽輪發(fā)電機定子內(nèi)流體場與傳熱藕合計算與分析［J］.電工技術(shù)學(xué)報,2009,24 (12):24-30. LIWeili,YANG Xuefeng,GU Debao,et al.Stator windy road in the air-cooled heat transfer fluid field and coupling calculation and analysis［J］.Transaction of China Electrotechnical Society,2009,24(12):24-30.

［9］ 孟大偉,何金澤,夏云彥.箱式Y(jié)JKK緊湊型中型高壓電動機全域流體預(yù)測［J］.電機與控制學(xué)報,2012,16(12):66-71. MENG Dawei,HE Jinze,XIA Yunyan.All-around fluid forecast research of the compactmedium size high voltagemotor［J］.Electric Machines and Control,2012,16(12):66-71.

［10］ 魏永田,孟大偉,溫嘉斌.電機內(nèi)熱交換［M］.北京:機械工業(yè)出版社,1998:146-148,190-205.

(編輯:劉琳琳)

Heat transfer analysis and cooler optim ization of YKK medium size high voltage asynchronousmotor

WEN Jia-bin1, ZHENG Jun1, YU Xi-wei2
(1.College of Electrical&Electronic Engineering,Harbin University of Science and Technology,Harbin 150080,China; 2.Nanyang Explosion Protection Group Co.,LTD,Nanyang 473000,China)

In order to improve the thermal dissipating performance of themotor and to improve the conditions of themotor temperature rise,a numerical analysismethod was applied to optimize the cooling performance of the motor cooler.For AC induction motor YKK closed structure,because of its structural characteristics,the temperature rise is one of key points of the design that is concerned.Firstly,taking a YKK450-4,800kW medium size high voltage asynchronousmotor was taken as an example.By creating its3D physical andmathematicalmodel of cooler of themotor and giving its basic assumptions and boundary conditions according to CFD theories,the fluid field and temperature field were calculated and analyzed,and their contours of cooler were shown.Secondly,by changing the numbers and the positions of the windshields,increasing the deflector and change its position,its optimal valueswere got.Finally,the cooler was remodeled.Numerical analysis shows that the performance of cooler is significantly improved,which provides theoretical basis for optimizing future design of cooler.

medium size high voltage asynchronousmotors;cooler;flow field;temperature field;optimal design

10.15938/j.emc.2015.09.005

TM 301.4

A

1007-449X(2015)09-0033-07

2014-05-21

國家自然科學(xué)基金(51275137)

溫嘉斌(1961—),男,博士,教授,研究方向為電機冷卻技術(shù)、電機及電機控制;鄭 軍(1989—),男,碩士研究生,研究方向為電機多物理場藕合分析計算;于喜偉(1986—),男,博士,工程師,研究方向為高壓異步電機設(shè)計。

溫嘉斌