某轴向无刷励磁机多方案通风三维热流场分析

1、某轴向无刷励磁机多方案通风三维热流场分析韩家德张羽楠刘延浩李会兰路义萍付秀兰摘要：以7. 8MW无刷励磁机现运行方案为基础，对通风冷却系统进行优化，提 出了以改变现方案下通风结构、空气入出口位置的多种方案。建立其整机流固 耦合物理模型，基于计算流体动力原理，采用有限体积法，给定边界条件，对 多方案下励磁机热流场进行数值计算并进行比照分析。结果说明：改变现方案 通风结构，励磁机内空气压力、速度、温度分布规律与现运行方案基本相同， 但定转子部件处平均温度不同，最高温度仍位于远离整流盘侧的转子线棒端 部，均低于原方案的最高温度；现方案下互换空气入出口位置，对励磁机的热流 场有较大影响，最高温度在靠近

2、整流盘侧的定子线棒端部。研究结果为结构类 似的励磁机的优化设计提供理论指导。Key：励磁机;全轴向通风；多方案;数值模拟;热流场DOI： 10. 15938/j. jhust. 2019. 04. 004 原方案基本相同，仅数值大小不同，转子部件的温度降低；另外，对励磁机增加 转子风沟数量和在转子部件前方绝缘板上额外增加通风孔，所到达的冷却效果 要比减小转子风沟数量同时增大风沟直径的方法略好。2)互换空气入出口位置，即冷却空气从原方案的出风筒侧流入定转子风道内， 经整流盘侧风筒流出，定子端部两侧空气区的压差变小，进入定子风沟的动力 和克服空气在定子风沟内流动阻力的动力变小，最终使得总风量和定子

3、风量变 少，相应地转子风量增大;定子部件温度高于转子部件温度，最高温度在靠近整 流盘侧的定子线棒端部。参考文献：1陈云华，解锦辉，彭兵.无刷励磁同步电机用旋转整流器设计船电技 术，2002 (2) ： 12.2朱荣，孙晖.阿尔斯通TKJ型无刷励磁机磁极线圈烧损原因及处理对策 J.大电机技术，2008 (5) ： 56.3戴彦德，戴琳.中国电机系统能源效率与市场潜力分析M.北京：机械 工业出版社，2001.4杨丽伟，张奕黄.机车用异步牵引电动机温度场分析J.微电机，2008,41 (8) ： 35.5邓日江，王冬梅.定子轴向通风孔对三相异步牵引电机温度场的影响J. 防爆电机，2012, 47 (

4、2) ： 30.6廖毅刚，侯小全.全空冷汽轮发电机通风冷却研究J.东方电气评论，2008(1) ： 1.7王芳，董惠，安志华等,大型汽轮发电机定子变结构对温度场的影响J.哈 尔滨理工大学学报，2013, 18 (6) ： 95.8周黎民，黄鹏程，郑国丽，等.防爆电机换热分析与结构改进J.电机与控 制应用，2012 (29) ： 7.9路义萍，孙雪梅，杜鹏，王佐民.CFD设置对凸极同步电动机转子温度场影 响J,哈尔滨理工大学学报，2015, 20 (3) ： 91.10孙维炎.应用在空冷汽轮发电机上的新型无刷励磁机组J.科技创新导报， 2009 (23) ： 133.11周凤武，龙欣苗，吴小波.

5、水轮发电机组无刷励磁机的研发J.湖北工业大 学学报，2013, 28 (1) ： 94.12史德利，王庭山，胡建波.配核电百万千瓦级发电5800kW无刷励磁机研制J.大电机技术，2014 (1) ： 69.13侯君，刘兆江.交流励磁发电机及旋转整流器的设计要点J.上海大中型电 机,2012, 3： 44.14张超，夏立.发电机旋转整流器故障的分形和动态测度诊断J.电机与控制 学报,2009, 13 (1) ： 6.15MCARDLE M G, MORROW D J. Noninvasive Detection of Brushless Exciter Rotating Diode Fai lu

6、reJ. IEEE Institute of Electrical and Electronic Engineers, 2004, 19 (2) ： 378.161CHHAVI N, PRATEEK K J, VIKRAM K R V.Rotor Vent Hole Shape Studies on Totally Enclosed Fan Cooled MotorC/Drives and EnergySystems. 2012 IEEE International Conference on Power Electronics, Bengaluru, 2012： 4672.17韩雪岩，祁坤.

7、特殊结构起重机用PMSM冷却系统的优化J.起重运输机 械，2015, (1) ： 1.18李伟力，王耀玉，黄东洙，等.转子通风结构对永磁电机转子流体场和温 度场的影响J北京交通大学学报，2015, 39 (2) ： 1.19孙铁源，蔡茂林.压缩空气系统的运行现状和节能改造J.机床与液压, 2010, 38 ( 13) ： 108.20路义萍，付余，张万全，等.一种新型无刷励磁机的三维热流场分析J.中国电机工程学报，2014, 34 (27) ： 4699.21路义萍，张东学，等.某无刷励磁机通风冷却数值模拟研究J .电机与控制 学报，2016, 20 (6) ： 26.22路义萍，戴景禄，张东

8、学，等.多旋转部件无刷励磁机改进方案三维流场 研究J.23路义萍，李会兰，韩家德，等.某轴向通风无刷励磁机三维流场数值模拟 J.哈豳滨理工大学学报,2016, 21 (4) ： 1.24陶文铃.数值传热学M.西安：西安交通大学出版社，2001： 1-9, 347.(编辑：关毅)一全文完一TM343；TK121文献标志码：A：1007-2683 (2019) 04-0022-07Abstract： According to the existing scheme of 7. 8MW brushless exciter , the ventilation and cooling system is

9、 optimized and a variety of schemes are proposed to change the ventilation structure, rotation direction and inlet position of the present scheme in this paper. The fluid-solid coupling model of the whole machine is established. Numerical calculation of heat flow field of exciter under multi-scheme

10、is made and compared by adopting the finite volume method and giving boundary conditions based on the principle of computational fluid dynamics. Changing the ventilation structure of the present scheme, the results show that the distribution law of air pressure, the velocity and temperature in excit

11、er are basically the same as that of current operation, the average temperature of stator and rotor parts is different, and the maximum temperature remains at the end of the rotor bar far from the rectifier panel and is lower than the maximum temperature of the original scheme. Exchanging inlet and

12、outlet position under present scheme, it has a great influence on the heat flux field of the exciter, and the maximum temperature is at the end of the stator bar near the rectifier plate. The research results provide theoretical guidance for the optimal design of the exciter with similar structure.K

13、eywords： exciter； full axial ventilation； multi-scheme； numerical simulation； thermal flow field0引言辗刷励磁机是利用整流盘进行整流为同步电机提供直流励磁电源的一种电机结 构1。在其运行过程中，整流盘装置产生热量、定转子铁心和线棒的损耗等使 得部件温度升高；通风冷却系统的不合理设计也会导致励磁机散热效果较差，以 上两方面均易造成磁极线圈过热烧损和整流元件高温失效2,影响运行性能和 经济指标3。因此，像研究异步电机、凸极电机等其他类型电机一样4-9, 在开发新型无刷励磁机时，为了保证励磁机内绕组、绝缘

14、和整流元件的温升在 其允许的温升极限范围内，需要研究冷却系统各因素变化对励磁机温升的影 响，进而有针对性地对励磁机进行方案优化。目前，关于励磁机设计原那么、结构特点口0-12以及励磁机中整流电路设计 13、故障诊断14-15方面的研究涉及较多，公开发表的文献针对励磁机通风 冷却系统多方案优化设计的热流场分析，即励磁机热流场影响因素分析较少， 大多是关于空冷汽轮发电机、异步电机、永磁电机等其他类型电机的热流场研 究。文16对采用轴向通风结构的感应电机在圆形和梯形两种转子通风孔形状 下的电磁场和温度场进行数值模拟，并分析了通风孔形状对其影响。文17对 端盖引风、定子铁心布置15的孔和轴内引风、定子

15、铁心布置中20的孔以及 空心轴强迫风冷三个通风方案下的永磁同步电机的温度场进行比照分析，得到 了较好的通风方案。文18为了降低转子温度，提出了在电机转子上开轴径向 通风沟的多种结构方案，研究了转子径向通风沟数量与位置对转子内温度分布 及空气运动的影响。电机通风冷却系统的合理设计可以有效带走电机内部产生 的热量，使电机的温度均匀分布19。本文研究的无刷励磁机，主轴上布置转子和整流盘两个旋转部件，采用轴向通 风系统即定转子铁心仅开设轴向通风沟（简称定、转子风沟），其通风结构 和采用的通风系统均不同于文20-22,属于国外引进二次开发的新产品。另 外，文23仅研究了 一种励磁机的流场。本文针对一种全

16、轴向通风励磁机提出 改变通风结构、互换空气入出口位置的多方案冷却系统，研究在不同冷却系统 下励磁机热流场分布规律，为结构相似的无刷励磁机通风冷却系统的设计提供 理论上的指导。1物理模型现运行方案下（以下简称为“原方案”）研究的7.8MW无刷励磁机与12MW凸极 同步电机同轴布置，该励磁机转速与主机相同。为显示该励磁机内部结构，图 1给出了省略进风筒、出风筒和方筒后热流场计算域物理模型的剖面图，图中 Z轴与转轴重合，计算域的原点与靠近整流盘侧的转轴端面的中心点重合，其 中，固体部件以实体显示，空气区以线框显示。2数学模型及求解条件数学模型对该励磁机内冷却空气的流动特征进行数值模拟时，仅研究稳态时

17、励磁机内的 空气热流场，其内部马赫数小于0.7,视空气为不可压缩流体，且其物性参数 恒定，励磁机中冷却空气受重力作用远远小于旋转科氏力作用，在数值模拟 时，忽略重力对空气流动的影响，热源在发热体中均匀布置。经屡次试算，励磁机入口处雷诺数Re=66334,定、转子风沟入口雷诺数Re分 别为8217、11142,气隙内主要是旋转剪切流，采用泰勒数描述流态，由式(1)算得气隙中旋转泰勒数Ta=54189,因此，励磁机内空气流动均为湍流。边界条件及求解设置根据工厂提供数据，现运行方案中，将励磁机的空气入口设置为压力入口，入 口表压为742Pa；空气出口设置为压力出口，出口表压为OPa；主轴转速为 10

18、00r/min；励磁机是从主机风扇后的流道中取风，按照电机规范，励磁机入口 的空气温度取50。绝缘、定转子直段铜绕组、转子支架、二极管座的热导率人均为各向同性，分 别为0.22、387.6、31.8, 202.4,定转子铁心的热导率人为各向异性，沿叠 片径向、切向和轴向的数值分别为42. 5, 42. 5和0.57,热导率单位均为W/(m , K) o3多方案数值模拟结果分析多方案介绍及风量分布计算结果方案A和方案B分别将原方案非均匀布置的22个20圆形轴向转子风沟变为 转子支架周向均匀布置24个020和12个32的转子风沟；方案C在原方案的 基础上将转子部件前方的绝缘板上布置6个中15的通风

19、孔；互换原方案的空气 入出口位置，冷却空气从原方案的出风筒侧流入励磁机内，沿轴向风道流动， 最后从原方案的进风筒侧流出，即为方案D。 在进行数值模拟时，为保证网格的统一，使各方案结果具有可比性，方案A、B 采用与原方案相同的网格类型和节点间距值对转子风沟、转子支架和相邻的区 域进行网格划分;方案C中绝缘板、整流盘空气区和相邻的区域使用和原方案相 同的网格类型和节点间距值进行网格划分。由于转子支架和绝缘板均不是发热 部件，因此，方案AC中发热体的热源密度与原方案相同。方案D只需将原方 案的入口设置为压力出口边界条件，原方案的出口设置为压力入口边界条件， 其他和原方案保持一致。在原方案基础上，对该

20、励磁机风冷系统分别采取方案AD措施后，通过数值模 拟分析，得出了各方案风量分布结果，图4给出了原方案和方案AD中励磁机 总风量、定转子风量和气隙风量分配及各风量所占总风量比的比拟。与原方案比照，方案AC总风量增大，方案D总风量降低，其中方案B相对各 方案总风量增加幅度最大；方案AD转子风量均增加，气隙风量均减小，其中 方案D相对各方案转子风量增加幅度和气息风量减少幅度最大，方案B次之； 方案AC定子风量与原方案相比变化不大，而方案D明显下降。关于各局部流量所占比例，通过数值分析，方案A定子和气隙风量占总风量的 比例减小，整体而言，各局部风量占总风量的比例变化不大；方案B进入定子风 沟的风量最多

21、，约占总风量的59.1%,转子风量和气隙风量基本相同，分别占 总风量的20. 4%和20. 5%；方案C转子风量占总风量的比例有所增大，但仍是三 局部风量中占总风量比例最少的，各局部风量所占风量的比例变化不大；方案 D中定子风沟的风量最多，占总风量的53.0%,其次是转子风沟的风量，占总风 量的31. 6%,气隙风量最少，仅占总风量的15. 4虬通风结构改变方案热流场分析 通过数值模拟，得出方案AC的压力分布规律均与原方案基本相同，仅数值大 小不同。其中，图5给出了原方案与方案A极角60和极角-90截面的压力 分布云图。由图5可知，在方案A和原方案中，旋转空气区的压力低于静止空 气区的压力，整

22、流盘侧定子端部空气区呈现随半径增大，压力升高的规律等。另外，结果还显示方案AC中进风筒中心处的压力分别约为298. 4Pa、301. 2Pa. 305. 9Pa,均略低于原方案进风筒中心处的压力315. 4Pa；结果说明， 方案A的励磁机空气入口和进风筒中心处的压差最大，使得更多的冷却空气流 入励磁机。与压力分布规律类似，通过数值模拟得出方案AC的速度分布规律均与原方案 基本相同，仅数值大小不同。通过以上风量分布分析，其中方案B总风量与转 子风量相对各方案均明显增大，但空气速度却明显降低，如图6所示。方案B 中定子风沟、转子风沟和气隙内的空气平均速度分别约为6.38m/s, 11.46m/s

23、和7.66m/s,均小于原方案。尽管方案B的转子风量为原方案的1.3倍，但转 子风沟的总通风面积由原方案的6908mm2增大为9646. 08mm2,约为原方案的 L4倍，因而，转子风沟内的空气平均流速减小；方案B的定子风量和气隙风量 比原方案小，两方案的定子风沟和气隙的结构相同，流量与流速成正比，因而 方案B的定子风沟和气隙内的空气平均速度均比原方案小。通过对温度场的分析，发现沿轴向和沿径向，方案AC与原方案的定子部件和 转子部件的温度分布规律基本相同，靠近整流盘侧的固体部件和空气的温度 高，另一侧温度低，仅数值大小不同等。如图7给出了原方案与方案B极角 90的温度分布，方案B的定子各局部的

24、温度升高，转子部件的温度略有降低，定转子部件的温度差异较小。为分析励磁机重要部件温度情况，图8给出了原方案和方案AC定子线棒直段(stator-zd)和定子棒端部(stator-db)及转子线棒直段(rotor-zd)和 转子线棒端部(rotor-db)的最低温度tmin、最高温度tmax和体平均温度 tavg大小比照图。方案A、B均是通过增大转子风沟的总通风截面，提高转子风量和转子风沟对流 换热面积，进而降低转子部件温度而提出的方案；方案C是通过绝缘板布置通风 孔，使整机风阻减小，进风量增加。由图8可知，方案A、B、C的转子部件最 高温度均较原方案有所降低，定子部件温度均有所升高，定转子部件

25、的温度差 异变小，总体来说，冷却效果均比原方案有所提高。模拟发现，通过减少转子 风沟数量而增大风沟孔的直径后，方案B的转子部件的最高温度仅比原方案降 低0. 4,定子最高温度升高5. 6。0而在原有结构基础上增加转子风沟数量， 最终方案A转子线棒最高温度比原方案降低0.9T,同时定子最高温度升高1.4,方案B定转子部件温度高于方案A；方案C转子最高温度降低1.4七，定 子最高温度升高2. 29。总体而言，方案A、C冷却效果相近，比方案B的冷却 效果略好。3.3入出口位置改变方案热流场分析通过数值分析，图9给出了方案D极角60。和极角-90。截面压力和速度分布 云图。从压力云图可知，方案D中的两侧定子端部空气区的压力分布较均匀， 不再呈现类似旋转空