感应电机定子温度场的数值计算.docx
《感应电机定子温度场的数值计算.docx》由会员分享,可在线阅读,更多相关《感应电机定子温度场的数值计算.docx(23页珍藏版)》请在冰豆网上搜索。
感应电机定子温度场的数值计算
第10卷第5期2006年9月
电机与控制
学报
ELECTRICMACHINES
AND
CONTROL
V01.10No.5
Sep.2006
感应电机定子温度场的数值计算
靳廷船1,
李伟力2,
李守法2
(1.天津阿尔斯通水电设备有限公司,天津300400;
2.哈尔滨理工大学电气与电子工程学院,黑龙江哈尔滨150040)
摘要:
针对中小型感应电机结构复杂热模型难以建立的问题,以一台小型感应电机为例,采用有
限元方法,建立了电机的二维温度场数值计算模型。
给出了定子二维温度场计算的泛函以及相关热性能参数的计算方法;计算了不同负载运行时样机定子的稳态温度场。
计算结果与实测值的比较,验证了所采用计算模型及方法的合理性。
该电机温度场计算模型可以应用到其他同类电机定
子温度场的计算与分析。
在该温度场计算模型的基础上,分析比较了机壳散热翅高度的变化对电
机定子温度场的影响,研究了气隙温度对定子温度场的影响。
关键词:
感应电机;数值计算;有限元;温度场
中图分类号:
TM301
文献标识码:
A
文章编号:
1007—449X(2006)05一0492一06
Numericalcalculationandanalysisof
stator
thermal
fieldin
an
inductionmachine
JINTing—chuanl,LI
Wei~li,LI
Shou—fa3
(1.Tianjin
ALSTOMHYDRO
Co.Ltd,TiaIljin
300400,China;
2.Collegeof
Electrical&ElectronicsEngineering,HarbinUniv,ofSci.Tech.,Harbin150040,China)
Abstract:
Aimingattheproblemofdevelopingthethermalfieldmodeldifficuldyofsmallandmedium—sizedinductionmotorsfortheircomplexstructures,the"two—dimensional(2-D)statorthermal
fieldmodelof
a
small—sizedinductionmotoriSestablishedbyusingthefiniteelement
method(FEM).The
functionalequationofcalculatingthe2-Dstator
thermalfieldispresentedandthemethodsofdetermining
somether—malparameters
are
given.Thesteadystator
thermal
fieldsofthemotorunderdifferentloads
arc
calculat—
ed:
andthevalidityofthe
thermal
fieldmodelandtherationalityofthecalculatedmethod
are
verifiedby
experimentalresultsderivedfrom
a
prototype
motor.Thethermalfieldmodel
can
beappliedtothestator
thermal
fieldcalculationandanalysisofotherinductionmotorswithsimilarstructures.Based
on
thepro—
posed
thermal
fieldmodel,someworkisgiven
to
theanalysisofthermalsensitivitytotheheightofthe
coolingribs.Theinfluenceoftheair—gaptemperature
on
the
stator
fieldisstudied.
Keywords:
inductionmachine;numericalcalculation;FEM;thermalfield
1引言
拿蓑惹裂毳銎笛鬣嚣蒜蒺慧奏篙
中小型感应电机的各项技术指标中,电机的温
题。
传统的温度计算方法主要是热路法…,它是用
收稿日期:
2006—03—23;修订日期:
2006一07一03
作者简介:
靳廷船(1980一),男,硕士,主要从事电机电磁场理论和电机温度场理论的研究;
李伟力(1962一),男,博士、教授、博士生导师,从事特种电机的研制、理论分析以及大型电机综合物理场的研究工作;
李守法(1982一),男,硕士,主要从事电机故障诊断和电机温度场计算与分析。
万方数据
第5期
感应电机定子温度场的数值计算
493
平均温升来衡量电机的发热情况,因此不能确定各部分的温度分布情况,仅是一个非常粗糙的解。
A卜mor和Chari在1976年首次将有限元法用在计算大型汽轮发电机定子铁心的稳态温度场的分析中心J,此后电机温度场的数值计算方法主要集中于大中型感应电机口。
41、特种电机”1和大型同步发电机上∞一71,由于中小型电机热模型结构复杂、绕组难以处理和对流换热系数不易确定等因素,数值计算方法很少采用。
本文利用二维有限元方法,在相应的假设条件及等效条件下,计算了小型感应电机不同负载情况下定子的稳态温度场,分析了气隙温度和机壳表面散热翅的变化对定子温度分布的影响,得出一些有益的结论。
2定子二维温度场的数值计算
2.1基本方程
普通Y系列中小型感应电机普遍采用全封闭自扇冷结构,电机内部无轴向通风冷却系统,各部分产生的热量主要通过定子铁心由机壳传递给周围的
空气,且转子旋转时电机内部空气随转子运动,热传
递过程复杂。
为简化分析,假设:
1)电机铁心沿轴向的温度梯度为零;2)机座的散热效果用散热翅散热近似体现;3)接线盒部位定子铁心外圆与静止空气接触部分为绝热面。
考虑实际安装尺寸,机壳和定子铁心之问加一层装配间隙,可用表示"1为
f=(0.5+3.Dn)×10一’
(1)
式中Do为定子铁心外径。
本文所采用的样机型号为Y80:
一2,额定功率为1.1kW,在上述条件下,可以取轴向的中心截面作为温度场的求解域,如图1所示。
图中Z:
、Z,为求解域的边界。
由于接线盒的存在,使接线盒区域风速基本为零,定子硅钢片可视为裸露在静止的空气中。
根据能量守恒原理和热传递的基本定律,求解域内二维温度场方程为
3’r
V・(A
V
r)+qy—pc导三=0
(2)
d7-
式中:
A为材料的导热系数(W/(m・℃));T为各点的温度(℃);g,为热源密度(W/m3);P为材料的密度(kg/m3);c为材料的比热(J/(kg・℃));丁为时间(S)。
本文计算电机的稳态温度场,温度不随时间变
3叮1
化,所以}=0。
在直角坐标系中式(2)可以表示为
图1求解域
Fig.1
Solution
model
A。
粤+A,磐:
Ⅶ
(3)
a石
oy
式中,A;、A,分别为材料在戈、Y方向的导热系数,本
文A。
=A,。
式(3)为热传导方程,它是用微分方程的形式来表达各物理量在相邻空间的数值之间关系。
此方程必须在给出空间的边界条件时,才有定解,常见的边界条件表达式可以分为三类:
第一类边界条件是已知任何时刻物体边界的温度,即:
T
1。
,=兀
(4)
式中2,为边界线,%为给定的温度值;
第二类边界是已知边界上的热流密度,即
A嚣卜_go
(5)
式中f:
为边界,q。
为通过边界的热流密度(W/m2),对于绝热边界条件,可令g。
=0,凡为边界上的法向;
第三类边界条件为边界周围流体的温度乃和散热系数d已知,即:
一A兰I=嫂(T一0)
(6)
Oft,I
f3
式中f,为边界,理为求解域边界的散热系数(W/(m2・℃)),0为周围介质的温度。
在电机求解域用到第二类和第三类边界条件,综合式(3)、(5)、(6),可得求解域内二维稳态温度场的边值问题为
击(A。
面OT)+喜(A,等)=一gr
A乳…。
A而OTb=一a(r-re)
(7)
由变分原理可知,方程(7)可写成如下的等价
变分方程
万方数据
494
电机与控制学报第10卷
们)=瓤【A,(警)2u(等)2】as一
上功y棚+扣上(r一2rs)rdL=min
丁=∑Nk瓦
(9)
[募aJ'二/;aT≥。
]=[k蔓k蔓曼][i]一[曼]=
a厂/a巧f:
I略
磅
‰ff乃I—I马I=
L
j
L耐
呵后一.|LL_j
Lp。
j
k“=咖(b;+c:
)
b=空(碍+弓)+警尼。
=①(6:
+c:
)+警
ko.=k=中(bibi+c。
q)
k抽=k。
i=西(bf
b。
+ci
c。
)
‰=‰=痧(bj
b。
+%)警
qvA。
伽i弓.qvA。
pi。
丁,功。
pm2彳+丁
其中:
多2去;△。
为三角形单元的面积;
Ⅱi
2了,一了m’呵2Ym—Yi,om2yf一乃’
b。
=省m—xj,q=戈;一菇m,bm=勺一戈i,
ci=巧ym一戈m”,cj=菇mYi—xiYm,
cm2瓮i,i—xi,t。
(筏,Yi)、(菇。
,Y。
)、(石f,乃)为三角形单元各顶点的坐标。
当泛函.,达到极值,即a∥aT=0时,由上述推导和单元分析知:
坼T=F
(12)
式中:
r为求解域内全部节点温度所形成的温度列
阵;坼为总体系数矩阵;F为总体右端列矢量。
再
将边界条件代入上式修改,最后求解修改后的方程组,即可得到各个节点的温度值。
2。
2热性能参数分析
一般中小型感应电机的定子绕组都采用圆形散线,线径细且各根导线在槽内排列不规则,实体建模极其困难。
为简化分析,可将槽内的铜线(不包括漆膜)等效地看作整个铜条;浸渍漆、槽绝缘和铜线
的漆膜近似看作绝缘层,等效之后的铜条位于槽的中心,四周与槽壁平行,浸渍漆和槽绝缘均匀的分布在铜条周围。
根据文[9],各种绝缘材料的等效导热系数为
8
,
“
譬
Aq=∑6,/∑}
(13)
i=1
’
i=1’‘i
式中:
A。
为等效导热系数;盈(i=1,2,3,…,n)为各导
热体厚度;Ai为各导热体的平均导热系数。
气隙中和机壳上的对流换热系数的确定可参阅文[9—10]。
电机运行时,各部分产生的电磁损耗和机械损耗直接导致了电机温度的升高,这些损耗主要有铁
心损耗、绕组铜耗、机械损耗和杂散损耗等,参数可由相关分析和试验方法得到。
温度场计算时分别将电机各部分热密施加到电机的有限元求解模型中,作为温度场计算时的热源。
3计算结果及分析
以Y80:
一2电机为例,求解时在不同的部位施加相应的热密,机壳表面和定子内圆赋第三类边界条件,接线盒对应的定子区域施加第二类(q。
=0)边界条件。
其中定子内圆气隙温度可用热敏电阻测
得:
在经过特殊处理的槽楔上放置体积很小的热敏
电阻,不与槽楔接触,使其仅与定子内侧的空气接触,可测量出定子内圆气隙的温度,作为添加第三类边界条件时的环境温度。
测量电机内部温度的热敏电阻放置点如图2所示。
A在定靠近槽口部位,B在槽楔和绝缘层中间,紧靠槽绝缘,A、B均位于铁心的轴向中部。
图2测温点分布Fig.2
Measured
pointsof
temperature
3.1不同负载情况下定子温度场的计算
额定负载下,电机的稳态温度分布如图3所示(环境温度17℃)。
47%额定负载下,电机定子的温度分布如图4所示。
其中Max和Min表示温度分布的最高点和最低点。
万方数据
第5期感应电机定子温度场的数值计算
495
图3额定负载时定子二维温度分布
Fig.3
Stator
temperaturedistributionatratedload
图4
47%额定负载时定子二维温度分布
Fig.4
Stator
temperaturedistributionat47%ofratedload
从上两图可以看出:
接线盒区域的绕组均比其他的部分偏高,与实际测量的结果一致。
其他几种负载情况下电机定子的温度分布趋势与图3和图4类似。
定子槽内温度的分布趋势如图5所示,其中温度最高的部分为等效铜条,整个区域几乎具有相同的温度,这是因为铜的导热系数较大,而且其本身也是热源;沿两侧绝缘层温度有明显的下降。
不同负载情况下,A点与B点的计算值与测量值之间的的关系如表I所示。
表1不同负载情况下定子的温度场
Table1
Stator
temperatureatdifferentloading
之间的误差较小,满足工程计算精确度。
采用有限元法计算中小型电机的定子温度场具有较高的准确性。
3.2气隙温度对电机定子温度的影响
电机运行时,定转子的各种损耗共同作用使气隙中的温度升高,因此定子温度场的准确分析与定子内圆空气温度的确定有很大的关系,不同负载下定子最高温度与定子内圆气隙温度的关系如图6~9所示。
p、
越赠m删
/最高值
/’
\
/
.
\/
//
√
//\最低值/
/定子内圆气隙温度,℃
图6空载时定子最高温度与定子内圆气隙温度的关系
Fig.6
Relationshipbetweenstatortemperatureand
airgap
temperature
near
thestatoratno-load
p、
毯赠m删
蚤离宿
./1
\,7
\/一.‘
一j,r
-一一一一
/
./
一/・7\目fm&
-/“‘”“
r.
.....
。
,图5.槽窒温竺2布。
.
警茎訾翥釜磊翥墨熹嚣窆鋈篝麓嚣茎妻
Fig.5
nmpe豫m弛豳倒加舶n村嘲‘盯310‘
花:
苫税爵磊爱羞圣孚二蓑磊基泛苁。
图苫‘。
。
9可以
从表1中可看出,A点和B点的计算值和测量值
看出:
随着定子内表面气隙温度的增大,电机定子的
万方数据
496
电机与控制学报第10卷
p、
型赠
H'
删
最I锃I倡
./’
~●一一一、
厂
.、一r
—,,
/1
-,,一
I—_r一,7\一...
-—一r
’最1或但
I.
定子内圆气隙温度/℃
图8
70%负载时定子最高温度与定子内圆气隙温度的关系
Fig.8
Relationshipbetweenstator
temperatureandairgap
temperature
near
thestatorat
70%ofratedload
p、
魁赠
怫
趔
最高值
.——一◆
\
.∥—一
.y。
最低值
/..一—一-
-I一—I一一
●—Jr
定子内圆气隙温度/0C
图9额定负载定子最高温度与定子内圆气隙温度的关系
Fig.9
Relationshipbetweenstator
temperatureandairgap
temperature
near
thestatoratratedload
最高温度和最低温度在数值上都有升高,最高温的升幅略大于最低温的升幅;当不同负载率下气隙温度都升高14℃时,空载、47%额定负载、70%额定负载和额定负载下定子区域最高温度分别上升了11.50℃、11.49℃、11.50℃和11.40℃,定子温度随气隙温度呈线性变化。
气隙的温度对定子温度场有显著的影响,初始条件下定子内圆气隙温度的准确
分析是精确求解定子温度的关键因素之一。
3.3散热翅高度对定子温度的影响
在全封闭式电机中,机座上设有纵向的冷却翅,借助于端部风扇的旋转,使具有冷却作用的空气沿散热翅流动,逐步减速扩散,冷却机座表面,降低电机内部的温度。
热流的大小与机座散热表面的面积成比例,散热翅的高度对定子温度的影响如图10~13所示。
当散热翅高度变化时,机壳上的对流换热系数也随之变化。
从求解的结果知,改变散热翅的高度并没有改变电机内温度的分布趋势,即最高温度点和最低温度点的位置不变,只是随着散热翅高度的变化,电机的温度发生了一定的变化。
从图10~13可以看出:
各种负载率下,在相应的定子内圆气隙内温度不变的情况下,当散热翅的高度变低时,定子各点温度变大,各部分温度呈上升的趋势,与散
p、
越赠
H’
删
p、
越赠№删
~~.
/最高值
~、、<
一’~—、.
\、、。
~.最低值
、\
/
~、≮\
、\.、
无散热翅l,3翅高2/3翅高全翅高1.5倍翅高2倍翅高
散热翅高度
图10空载时机座散热翅高度对定子温度的影响
Fig.10
EffectsofcooHngribsheighton
stator
temperature
atno-load
~~‰、1
r~\.
最高值
~’\
、、、。
,最低值
、’~/
\、、,、
.\、!
无散热翅l/3翅高2/3翅高全翅高1.5倍翅高2倍翅高
散热翅高度
图1l
47%额定负载时机座散热翅高度对定子温度的影响
Fig.11
EffectsofcooHngribsheight
on
stator
temperature
at47%ofn040ad
58
54
芝50
嚣46渔42
38
34
k、÷,j
\最高值
一—“~———~.
.
_\\、
~“、~~-、.
/最低值
、’。
“、<。
,一
~”\,
无散热翅l,3翅高2/3翅高全翅高1.5倍翅高2倍翅高
散热翅高度
图12
70%额定负载时机座散热翅高度对定子温度的影响
Fig.12
Effectsofcoolingribs
height
onstator
temperatureat70%ofno-load
80
7672
芝68
赵64
鲣
48
-~..
。
~、卜~
、最高值—~~’ ̄~
.
^“*、、。
\~、
、、。
、
,最低值、≮,、。
、一\-
.、’\-、、、..
弘
娩
∞勰拍抖
稻钳舵∞弛拍M亚如
万方数据
第5期感应电机定子温度场的数值计算参考文献:
[1]497热翅高度呈线性关系;空载时,在气隙温度不变的情况下,在无散热翅时定子最高温度比全翅高时定子最高温度升高了1.77℃,升幅为5.74%;定子最高杨国强.封闭式小型感应电机的等值热流图和温升计算方法[J].中小型电机,1978,(4):
12—19.温度在两倍散热翅高时比全翅高时降低了0.77℃,降幅为2.50%;47%额定负载时,对应的升幅和降幅分别为5.26%和3.92%;70%额定负载时,对应的升幅和降幅分别为6.46%和4.46%;额定负载时,对应的升幅和降幅分别为8.25%和5.46%。
从以上分析可知:
散热翅的高度对定子温度场有较大的影响,与标准散热翅高度相比,当散热翅高度降低时,定子温度上升明显,散热翅高度升高时,定子温度下降,下降的幅度有限。
合理的设计电机散热翅的高度,并从外观方面综合考虑,选择合适的散热翅高度,对电机的散热性能有一定的改善作用,从而提高电机的运行性能。
[2]ARMORAF,CHADMVK.Heatflowinthestatorcoreoflargeel--turbine—-generatorsbythe"methodements:
Partofthree—-dimensionalfinitel:
Analysisinbyscalarpotentialformulmion:
PartII:
Temperatllmdistributionthestatoriron[J].IEEETrans,PAS一95,1976,(5):
1648—1668.[3]李伟力,宋凯,付敏,等.用三维有限元法计算大中型感应电动机定子三维温度场[J].黑龙江大学自然科学学报,1999,16(2):
33—36.[4]李伟力,付敏,周封,等.基于流体相似理论和三维有限元法计算大中型感应电动机的定子三维温度场[J].中国电机工程学报,2000,20(5):
14—17,21.[5]黄东洙,沈稼丰,侯云鹏.高压防暴电机定子三维温度场的计算与分析[J].防爆电机,2002,(4):
11—13.[6]栾茹,傅德平,唐龙尧.新型浸润式蒸发冷却电机定子三维温4结论通过对一台小型感应电机不同负载时定子稳态温度场的计算,以及气隙温度和散热翅高度对定子温度的影响的分析,得出如下结论:
1)在求解过程中做了适当的假设,可以把电机定子复杂的三维温度场问题简化为二维问题求解,减小了工作量,同时采用二维有限元数值计算方法计算电机绕组的温升,具有较高的计算精确度;2)随着定子内圆气隙温度的升高,定子各部分温度呈线性上升趋势;3)定子温度与散热翅高度是线性变化关系,散热翅的高度对定子温度有较大的影响。
[9]度场的研究[J].中国电机工程学报,2004,24(8):
205-209.[7]李伟力,侯云鹏,周封,等.大型水轮发电机定子股线导热的数值分析[J].中国电机工程学报,2001,21(7):
115—119.[8]黄国治,傅丰礼.Y2系列三相感应电动机技术手册[M].北京:
机械工业出版社,2004.AM鲍里先科,Br丹科.电机中的空气动力学与热传递[M].魏书慈,邱建甫译.北京:
机械工业出版社,1985.[10]XYPTERASJ,HATZIATHANASSIOUV.Thermalanalysisofanelectricalmachinetakingintoaccounttheironlossesandthedeep—bareffect[J].IEEETransactionsOltEnergyConversion,1999,14(4):
996—1003.[11]龚晓峰,刘长红,饶方权,等.特种感应电机转子温度场的计算[J].大电机技术,2004,(5):
13一16.万方数据