基于Ansys的11kW无刷直流电动机温度场分析.pdf

基于Ansys的11kW无刷直流电动机温度场分析 高彦骋,等中图分类号:T M36+1 文献标志码:A 文章编号:100126848(2008)0920013203基于Ansys的11k W无刷直流电动机温度场分析高彦骋,刘卫国(西北工业大学 自动化学院,西安 710072)摘 要:利用Ansys软件研究无刷直流电动机在运行时的发热情况和温度场分布。仿真结果和实验,结果基本一致,说明本文在电机热场分析中采取的方法对于电机设计和分析具有一定的工程应用价值。关键词:Ansys;瞬态热分析;温度场;无刷直流电动机Analysis of Temperature Field for BLDCM Based on AnsysGAO Yan2cheng,L I U Wei2guo(Northwestern PolytechnicalUniversity,Xian 710072,China)Abstract:Made a transient state thermal analysis of the temperature field of BLDCM with Ansys soft2ware,and studied the heat2emitting conditions of running electricalmotor and simulation of its distribu2tion of temperature field.By contrasting the si mulation result and the experimental data,demonstratedthe method which used the ther mal2analysis module of Ansys software had a certain project applicationvalue in the designing and analysis ofmotor.KeyW ords:Ansys;Transient state ther mal analysis;Temperature field;BLDCM收稿日期:20072102300 引 言稀土永磁无刷直流电机温升过高会影响它的运行状态,会影响耐热能力薄弱的绝缘材料,使其寿命缩短,甚至烧毁,所以在电机设计和分析中了解其发热情况,比较准确地计算温度场具有实际意义。但是由于电机结构、散热条件等因素,给温度的计算带来一定的困难,而且目前传统的温度场分析大多数都从整体考虑,计算的结果多为平均值。利用Ansys有限元计算软件并通过其后处理器post26和post1的应用就可以比较好地解决这一问题。本文通过Ansys软件中热分析模块的应用对实际电机温度场进行二维分析;通过温度分布云图和温升曲线来了解电机在运行时的温度场分布情况。1 电机材料和热参数选取111 材料选取电机各部分材料分别如下所示。轴:钢;绕组导线牌号:QY-2;磁轭材料:DW465-50;磁钢材料:NSC27G;其物理特性为:硅钢片材料密度:7165千克/立方厘米;磁钢材料密度:813千克/立方厘米;导线密度:819千克/立方厘米;空气的密度:11165千克/立方厘米;硅钢片导热系数:40;磁钢导热系数:119;导线导热系数:376;空气导热系数:010267112 对流率由于对流和导热是电机热量传递的主要途径,应用Ansys软件时把对流率作为主要的边界条件。对流传热现象不仅取决于流动物体的运动性质、运动速度和固体表面形状,还取决于流体的物理性质。运用的基本公式为牛顿换热定律:q=a(T1-T2)=aT(1)式中,q为热流密度,T1为发热体表面温度,T2为流体介质温度,a为对流换热系数。(1)定子内圆与线圈表面的传热系数在此过程中可以近似无轴向气流,所以可按照气流掠过大平壁的计算方法来处理1:ReL=LV(2)Nu=0168(ReL)12(Pr)13(3)ah=NugL(4)31微电机2008年第41卷第9期式中,为风速;L为风道长度;Pr为空气的普朗特常数;ReL为雷诺数,包含着流速的影响;Nu为雷诺数,包含着流体流速的影响;Nu为努塞特数;ah为传热系数;g为重力加速度;为空气的热传导系数。(2)电机气隙表面散热系数电机气隙的冷却一面受转子切向运动的影响,一面则受定子内圆表面的阻滞作用。因此,气隙中的气流切向速度呈双曲线分布。根据经验公式,其表面的传热系数可表示为:a=28(1+015)(5)式中,a为气隙表面传热系数W/(m2k);为气隙平均风速,一般12u2;u2为转子圆周速度。(3)当电机机座壁向周围空间自然传热时,假设机座壁外表面温度定义和室外温度相同,均为初始温度,则表面传热系数应按下式计算:a=a0(1+kv)325(6)式中,a0为发热表面在平静空气中的换热系数,近似等于14;v为空气吹拂外表面的速度;k为气流吹拂效率的系数,近似等于015;为机座壁外表面的温度。代入以上数据,公式(6)写为:a=14(1+015i)325如果吹拂机座壁的风速为0 m/s,则传热系数公式(6)变为:a=14325113 生热率生热率的物理定义是单位时间内单位体积中内热源的生成热:Q=Wq/V(7)式中,Wq为热损耗;V为材料体积。在电机发热过程中主要的内热源就是铜耗和铁耗。(1)铜耗计算PCu=(I2xRx)(8)式中,Ix为相绕组电流的有效值;Rx为相绕组电阻值。(2)铁耗的计算:PFe=KapFeGFe(9)式中,Ka为由于硅钢片加工,磁通密度分布不均,以及其不随时间正弦变化等原因而引起损耗增加的系数;pFe为单位质量的损耗,也称比损耗;GFe为铁心净用铁量。2 电机发热过程计算本文应用Ansys软件中的热模块对实际电机进行瞬态热分析,其中利用APDL语言进行编程。图1为热分析流程图。图1 热分析流程图文中所用电机为11 kW稀土永磁无刷直流电动机。额定电压:270V-DC;额定转速:8 000 r/min。为了简化分析,假设:(1)本文是二维温度场分析,无法考虑轴承损耗对温度变化起的作用,所以机械损耗忽略不计;(2)铁心损耗主要为定子铁心损耗,忽略不计转子铁心中产生的铁心损耗。对该样机进行的瞬态热分析,具体步骤为:(1)建模确定工作文件名,热分析单元、单位,进入前处理器。(2)构建几何模型主要尺寸:电枢内径:918 cm;电枢外径:1515 cm;极对数:2;气隙宽度:0112 cm;转子铁心长:9 cm;磁钢磁化方向厚:1 cm结构参数:定子齿数:36;磁钢类型:瓦片径向充磁;定子齿槽形:开口底半梨形槽。图2为依据以上参数所建立的几何模型。(3)输入材料属性:定义单元类型并设置单元参数,按照上节中的材料参数设定。本文中温度场分析用的是plane55热分析单元,其中维数是2维,自由度是温度,节点数目是4个。(4)剖分:41基于神经网络的永磁同步电动机新型矢量控制器设计 王 剑,等 划分网络,生成有限元模性。电机的网格划分比较复杂,不同的区域的精度不同,温升变化大的地方网格需要细致划分,其他地方可以粗略一些。经过反复尝试,最终的网格划分如图3所示。(5)施加边界条件和载荷,加载求解直接在实体模型上施加的热荷载包括初始温度、热流率、对流、热流密度和生热率。其中温度作为自由度约束施加于温度已知的边界上;对流边界条件作为面荷载施加于实体的外表面,计算流体的热交换;生热率作为体荷载施加于单元上,可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。边界条件主要是先赋各个节点的初始温度值。此时取常温24。完成温度赋值之后,添加对流条件,还有生热率这些热载荷。通过上节公式,可分别求出各部分的传热系数和生热率,在运算时,根据实际电机设计的额定铜耗铁耗算出生热率。其中,Pcu电枢铜损:136瓦,PFe电枢铁损:535瓦。将结果代入上节中的公式,施加到边界条件中,进行下一步计算。(6)设置瞬态热分析初始条件设定荷载步选项。荷载步选项包括每一个荷载步结束时的时间、每个荷载步的荷载子步数、迭代次数、自动时间步长、时间积分效果等。每一个子步默认的迭代次数为25次。这对大多数非线性热分析足够用,所以自动时间步长选项打开。Ansys进行热分析的过程实际上是先将待处理的对象划分成有限个单元,其中包括了若干的节点,然后基于传热学经典理论能量守恒原理的热平衡方程,求解一定边界条件和初始条件下每一节点处的热平衡方程,由此计算出各节点的温度,进而求解出其它相关量。(7)后处理进入通用后处理器,对瞬态热分析的结果进行后处理。Ansys提供两种后处理方式进行瞬态热分析的结果后处理。通过通用后处理器(post1),可以对整个模型在某一时间点的结果进行处理。通过时间历程后处理器(post26)可以看出变量随时间变化的曲线。见图4。3 仿真结果与实际情况比较分析从仿真云图(图5)中可以看出温度场仿真的最高温度出现在线圈绕组附近,大约94。这是因为绕组绝缘及绝缘介质的导热系数很小,而且线圈绕组是一个大功率损耗的发热源。在中间的气隙温度也很高,因为电机在高速运转时,产生很大的风磨损耗,而且靠绕组比较近,不能及时将热量散出。其温度约为95。电机运行时发热均来自其损耗。其中定子和绕组既是发热部件又是传热部件,别的部件是传热部件。随着电机温度的升高,它与周围介质之间的温度差逐渐增大,散失到周围介质中去的热量相应的也逐渐增多,本身温升升高的速度则逐渐变慢,最后电机产生的全部热量都传给周围介质,达到了热稳定状态。通常每小时温度变化小于1 时,即认为已达到热稳定状态。从温升曲线(图6)中可以看出样机在运行一段时间后的温升最大值达到了94。实验中测出的数据为运行10分钟以后也达到9611,所以计算结果与实验结果基本相一致,因此文中应用Ansys有限元软件热分析模块处理电机温升的方法可为电机的设计提供参考。图6 电机温升仿真曲线(下转第48页)51微电机2008年第41卷第9期3 结 论目前,国内在对超声电机的研究集中在研制、理论分析、实验研究、性能提高等方面。控制技术方面与国外相比,还有很大的差距,因此我们认为今后应着重于以下几个方面的研究。(1)从控制角度出发,根据超声电动机的速度特性,进行系统辨识,直接建立数学模型,并根据这种数学模型开发超声电动机的控制系统。(2)基于现代控制理论,寻求更好、更简单的控制方法,如最优控制、基于参考模型的模糊神经网络控制方法、遗传算法、混合控制等,使多种控制技术有机集成。(3)为实现良好的控制,对超声电机模态特性和输出特性的测试,是实现各种控制策略的先决条件,因此今后应加强对超声电机测试技术方面的研究。参考文献1 赵向东.旋转型行波超声电动机动力学模型及性能仿真的研究D.南京:南京航空航天大学,2000.2 胡敏强,石斌,钱俞寿.超声波电动机控制技术研究的发展及现状J.微电机,1998,31(2).3 雷伏容,程昱,钟宜生,等.超声马达数学建模研究现状与展望J.压电与声光,2003,25(1).4 张铁民,赵淳生.超声电机伺服控制系统的研究评述 J.振动、测试与诊断,2001,20(3).作者简介:袁义坤(1978-),女,山东潍坊人,硕士,讲师,从事超声电动机动力学建模及其控制技术研究。(上接第15页)4 结 语利用Ansys软件,根据本文介绍的方法对电机进行热场分析和计算,全面了解其内部具体温升情况。从样机温度场的仿真分析结果看出,电机模型各部分的温度分布趋势与理论分析基本一致。由此可见,用文中所提出的方法对电机进行热场分析是行之有效的。今后还可以继续应用Ansys自带的APDL语言编写代码,为同类型电机的温度场分析提供便利。对于其他电机也能用类似的分析方法来处理,这对于深入认识电机内部的温升情况有一定的指导意义。参考文献1 魏永田,孟大伟.电机内热交换M.北京:机械工业出版社,1998.2 唐兴伦,范群波.ANSYS工程应用教程M.北京:中国铁道出版社,2003.3 Xyptras,J.,Hatziathanassiou,V.Thermal Analysis of an E2lectricalMachine Taking into Account the Iron Losses and theDeep2bar effect J.IEEE Transactions on Energy Conversion,1999,14(4):99621003.4 SaeedMoaveni.有限元分析 2ANSYS理论与应用M.北京:电子工业出版社,2003.5 薛蕾,刘卫国.ANSYS软件在高功率密度异步电机定子热分析中的应用J.微电机,2006,39(1).6 丁文,周会军,鱼振民.基于Ansys的开关磁阻电机温度场分析J.微电机,2005,38(5).7 靳廷船,李伟力,陈文彪.感应电机运行时的定子温度场分析与计算J.防爆电机,2006,(2).8 孙建宏,丁文,鱼振民.扁平形直线感应电动机温度常的计算与分析J.电机与控制应用,2006,(1).作者简介:高彦骋,女,硕士研究生,研究方向为电机与电器。刘卫国,教授,博导,研究方向为运动控制,电气伺服控制。84