摘要
针对高铁大型永磁同步牵引电机转子区域散热困难、内部温升分布不均匀的问题,在机壳水冷电机结构基础上的定子铁芯外表面处增设轴向矩形风道,并结合气隙、转子减重孔形成内外双循环散热结构,探究降低定子、绕组部位温升和提升电机内部散热均匀性的影响规律.首先通过Ansoft Maxwell平台仿真得到双循环散热结构在额定工况下各部件损耗值,同时为了更好模拟转子旋转带动气隙中的空气流动,对气隙进行分层处理,采用流-固耦合的有限元分析法模拟研究单、双循环散热结构下电机内部空气流动特性以及温升规律.结果表明,内循环风冷结构使电机内部空气流速显著提高,表面平均换热系数也显著提升,转子区域的热量会随着空气的流动更多地传给温度相对较低的定子区域及机壳,同时减少热量向转子部位传递,从而使转子和永磁体的温升降低.在此基础上,采用正交分析法对矩形通风孔的截面积、数量、高宽比进行结构参数优化,并采用温升分布均匀性系数对电机进行温升评价,得到最优方案下电机最高温升相比单循环散热结构降低12.1K,电机整体温升分布均匀性提升16.54%.
近年来,永磁同步电机凭借其高效性、高功率密度等优势在轨道交通、风电和汽车等领域得到大力发展和应
目前国内外学者针对永磁同步电机散热问题已开展大量研究,Jaege
针对风冷结构电机,Huang
本文以一台1 000 kW全封闭大型永磁同步牵引电机为例,针对永磁体温升过高问题提出内/外双循环散热方案,在机壳水冷方式基础上增加内部循环风散热,通过优化定子风道结构参数提升电机散热性能,为全封闭大型永磁同步电机温升分析方法及散热结构优化设计研究提供参考.
1 电机求解模型确立
1.1 电机的基本参数
本文所用电机为1 000 kW、1 500 r/min大型永磁同步牵引电机,基本结构参数如
参数 | 数值 | 参数 | 数值 |
---|---|---|---|
额定电压/V | 1 140 | 定子槽数 | 72 |
额定功率/kW | 1 000 | 转子内外径/mm | 180/444 |
额定转速/(r·mi | 1 500 | 气隙长度/mm | 3 |
电机极数 | 6 | 铁芯长度/mm | 500 |
定子内外径/mm | 450/670 | 绝缘等级 | F级 |
永磁同步电机物理模型如

(a) 样机内部结构(1/2模型)

(b) 双循环散热结构(3/4模型)
图1 电机结构示意图
Fig.1 Schematic diagram of motor structure
针对机壳水冷单循环样机转子区域散热困难、内部温升分布不均匀现象,提出内/外双循环(水冷和风冷)散热方案,即在电机一端主轴上增设离心风扇,在定子铁芯的外表面开设轴向矩形通风孔,并将转子铁芯原有的减重孔利用为转子轴向风道.在离心风扇作用下,转子铁芯区域高温空气由转子铁芯轴向风道和气隙向电机另一端移动,从定子铁芯表面开设的矩形通风孔返回,形成内部循环风冷,实现热量传递给定子铁芯和机壳(
1.2 基本假设和边界条件
根据永磁电机有限元分析计算流程要求,结合其散热结构特点和文献[
边界条件设置为:1)水道入口定义为速度入口,水速选用0.8 m/s,水温338 K;水道出口定义为标准压力出口.2)环境温度给定300 K.3)转子与旋 转流体域交界面设置为旋转壁面,转速设置为 1 500 r/min.4)利用UDF程序将风路出口温度附加给风路入口温度,以此来模拟循环风冷.
1.3 热源计算方法
在Ansoft Maxwell中使用二维瞬态磁场模型计算得到电机额定工况下各部件损耗
名称 | 损耗值/kW | 热生成率/(W· |
---|---|---|
定子铁芯 | 7.469 | 173 568.2 |
定子绕组 | 5.245 | 184 333.3 |
转子铁芯 | 0.526 | 114 166.7 |
永磁体 | 0.02 | 1 908.54 |
杂散损耗 | 5.128 | — |
1.4 气隙导热系数的确定
为了更准确地求解气隙导热系数,将气隙处空气进行分层处理,靠近定子铁芯的设为静止区域,导热系数按正常空气导热系数确定;靠近转子铁芯设为旋转区域,导热系数则引入气隙等效导热系数.气隙流场截面如
(1) |
式中:,为转子铁芯外径,为定子铁芯内径;为气隙厚度,mm;为电机额定转速;为空气运动黏度系数.

图2 气隙剖面示意图
Fig.2 Schematic diagram of air gap section
结合本文研究电机实际结构参数及额定转速值,由
1.5 网格无关性和模型验证
保证计算边界条件相同前提下,分别使用534万、622万、706万、798万四种数量网格开展温升模拟计算.通过比较同一条外层绕组沿轴向各点温度值变化来验证网格无关性.温度计算结果如

图3 网格无关性验证
Fig.3 Grid independence verification
企业所搭建的永磁同步电机温升测试平台如

图4 试验测试平台
Fig.4 Test and test platform

图5 实验温升测点示意图
Fig.5 Experimental temperature rise measurement point diagram
测温位置 | 温升值/K | 误差/% | |
---|---|---|---|
试验结果 | 模拟结果 | ||
绕组端部 | 77.5 | 80.78 | 4.2 |
转子端部 | 122.6 | 128.59 | 4.9 |
永磁体端部 | 121.7 | 127.54 | 4.8 |
2 散热结构仿真分析及优化
2.1 单循环结构流场和温度场仿真分析
在额定工况下电机内部空气速度矢量图如

图6 电机内部空气速度矢量图
Fig.6 Vector diagram of air velocity inside the motor
电机中心截面温度分布情况如

图7 中心截面温度云图
Fig.7 Central section cloud image
2.2 双循环散热结构流场与温度场分析计算
通过对内外双循环散热结构电机的温升特性模拟,得到电机内部空气速度分布,如

图8 电机内部空气速度矢量图
Fig.8 Vector diagram of air velocity inside the motor

图9 电机中心截面温度分布
Fig.9 Temperature distribution of motor center section

(a) 转子铁芯

(b) 永磁体

(c) 定子铁芯

(d) 绕组
图10 双循环方案电机各部件温度分布
Fig.10 Temperature distribution of each part of the double cycle motor
2.3 两种冷却方案对比
由

图11 不同方案各区域空气流速对比
Fig.11 Comparison of air flow rates in different regions
对2种冷却方案中电机各部件的最高温升计算结果进行对比分析,由

图12 不同方案各部件最高温升对比
Fig.12 Comparison of the maximum temperature rise of each component in different schemes
2.4 定子风道结构参数优化
基于双循环散热结构具有良好性能,对定子风道结构参数做进一步优化.以电磁转矩的减小量低于3%为约束条件,结合电机实际结构尺寸和散热要求,选取3种方案.其中定子风道结构参数的改变会引起定子铁耗及电磁性能的变化,因此需要利用Ansoft Maxwell软件重新进行电磁分析计算,并将所得损耗值转化为温度场求解热源,其余温度场求解边界条件均相同.具体参数见
风道总截面积S/m | 风道数量 | 高宽比 | 定子铁耗/W | 电磁转矩/(N·m) | 电磁转矩减小百分比/% |
---|---|---|---|---|---|
2 500 | 6 | 3∶1 | 7 539.15 | 6 094.7 | 1.24 |
3 500 | 12 | 1∶1 | 7 751.65 | 6 056.4 | 1.7 |
4 500 | 18 | 1∶3 | 8 325.37 | 5 981.9 | 2.59 |
采用正交因素法确定的27种组合方案的温度场计算结果如

(a) 电机整体最高温升

(b) 定子铁芯平均温升

(c) 定子绕组平均温升

(d) 永磁体平均温升
图13 不同结构参数电机温升分布
Fig.13 Temperature rise distribution of motor with different structural parameters
分析
2.5 电机温升特性评价
最佳定子风道结构方案下电机各区域的平均对流换热系数如
位置 | 平均对流换热系数/W/( | 提升百分比/% | |
---|---|---|---|
单循环结构 | 双循环最佳 冷却结构 | ||
左端部 | 1.49 | 6.53 | 338.3 |
气隙 | 26.74 | 30.41 | 13.7 |
右端部 | 1.89 | 6.97 | 268.8 |
最佳定子风道结构方案电机和单循环结构电机内部径向温升分布对比如

图14 不同散热结构电机中心截面中间位置径向温升分布
Fig.14 Radial temperature rise distribution in the middle position of the center section of motors with different cooling structures
为了量化评价不同方案下电机整体温升分布均匀性,引入温升分布均匀性系数,具体计算方法如
(2) |
式中:为温升分布均匀性系数;为电机整体平均温升值;为采样点个数;为采样点温升值.
通过仿真数据计算得到,单循环散热方案下电机温升分布均匀性系数为32.46,而基于内外双循环散热结构,同时优化定子风道结构参数,所得到最佳定子风道结构方案均匀性系数为27.09,相比单循环方案极大改善了电机整体温升分布状况,使得电机温升分布更加均匀,电机整体温升分布均匀性提升了16.54%.
3 结 论
本文采用有限元分析法对一台1 000 kW大型永磁同步牵引电机进行研究,设计了一种内外双循环散热结构,并进一步对定子风道结构参数进行优化分析,得出如下结论:
1)单循环机壳水冷电机存在温升分布不均匀、转子区域散热困难以及永磁体温升过高问题,其中定/转子最高温升差值为40.6 K,永磁体最高温升129 K.
2)双循环散热方案极大程度改善电机内部空气的流动性,有利于转子区域散热、降低永磁体温升,定/转子最高温升差值减小23.6 K,电机整体温升分布均匀性得到改善.
3)通过定子风道结构参数多因素综合分析,最佳定子风道结构参数方案为:风道总截面积3 500 m
参考文献
付兴贺,江政龙,吕鸿飞,等.电励磁同步电机无刷励磁与转矩密度提升技术发展综述[J].电工技术学报,2022,37(7):1689-1702. [百度学术]
FU X H,JIANG Z L,LÜ H F,et al.Review of the blushless excitation and torque density improvement in wound field synchronous motors[J].Transactions of China Electrotechnical Society,2022,37(7):1689-1702.(in Chinese) [百度学术]
周志刚,杨文豪,孟祥明.不同冷源轮毂电机多模式切换温度场研究[J].湖南大学学报(自然科学版),2021,48(8):51-58. [百度学术]
ZHOU Z G,YANG W H,MENG X M.Research on temperature field of in-wheel motors with different cold sources in multi-mode switching[J].Journal of Hunan University (Natural Sciences),2021,48(8):51-58.(in Chinese) [百度学术]
赵纪龙, 逯卓林, 韩青峰, 等. 轴向磁通永磁电机系统及关键技术前沿发展综述[J]. 中国电机工程学报, 2022, 42(7):2744-2765. [百度学术]
ZHAO J L, LU Z L, HAN Q F, et al. An overview on development of axial flux permanent magnet motor system and the key technology[J]. Proceedings of the CSEE, 2022, 42(7):2744-2765.(in Chinese) [百度学术]
龙吟,任晓辉,张珂,等.基于模态实验的轨道牵引电机整机有限元模型的建立[J].铁道科学与工程学报,2019,16(6):1553-1559. [百度学术]
LONG Y,REN X H,ZHANG K,et al.Finite element modeling of rail traction motor based on modal experiments[J].Journal of Railway Science and Engineering,2019,16(6):1553-1559.(in Chinese) [百度学术]
DU G H, XU W, ZHU J G, et al. Power loss and thermal analysis for high-power high-speed permanent magnet machines[J]. IEEE Transactions on Industrial Electronics, 2020, 67(4): 2722-2733. [百度学术]
JAEGE R M,RUF A,HAMEYER K,et al.Thermal analysis of an electrical traction motor with an air cooled rotor[C]//2018 IEEE Transportation Electrification Conference and Expo (ITEC).Long Beach,CA,USA.IEEE,2018:467-470. [百度学术]
DONG B T,WANG K,HAN B C,et al.Thermal analysis and experimental validation of a 30 kW 60 000 r/min high-speed permanent magnet motor with magnetic bearings[J].IEEE Access,2019,7:92184-92192. [百度学术]
赵纪龙, 全小伟, 林明耀.双转子混合励磁轴向磁通切换永磁电机设计与分析[J].中国电机工程学报,2020,40(24):7860-7868. [百度学术]
ZHAO J L,QUAN X W,LIN M Y.Design and analysis of a double-rotor hybrid excited axial switched-flux permanent magnet machine[J].Proceedings of the CSEE,2020,40(24):7860-7868.(in Chinese) [百度学术]
VANSOMPEL H, SERGEANT P.Extended end-winding cooling insert for high power density electric machines with concentrated windings[J].IEEE Transactions on Energy Conversion,2020, 35(2): 948-955. [百度学术]
刘威.非晶合金高速永磁电机热管理与转子动力学研究[D].太原:中北大学,2018. [百度学术]
LIU W.Thermal management and rotor dynamics of high speed permanent magnet motor with amorphous alloy[D]. Taiyuan:North University of China,2018.(in Chinese) [百度学术]
TÜYSÜZ A, MEYER F, STEICHEN M, et al.Advanced cooling methods for high-speed electrical machines[J]. IEEE Transactions on Industry Applications,2017,53(3):2077-2087. [百度学术]
HUANG Z Y,FANG J C,LIU X Q,et al.Loss calculation and thermal analysis of rotors supported by active magnetic bearings for high-speed permanent-magnet electrical machines[J].IEEE Transactions on Industrial Electronics,2016,63(4):2027-2035. [百度学术]
ZHANG X T,ZHANG C M,FU P R,et al.A novel cooling technique for the windings of high-torque-density permanent magnet machines[C]//2018 21st International Conference on Electrical Machines and Systems (ICEMS).Jeju,Korea (South).IEEE,2018:332-337. [百度学术]
FAWZAL A S,CIRSTEA R M,GYFTAKIS K N,et al.Fan performance analysis for rotor cooling of axial flux permanent magnet machines[J].IEEE Transactions on Industry Applications,2017,53(4):3295-3304. [百度学术]
佟文明, 程雪斌, 孙静阳, 等.转子风刺对高速永磁电机永磁体温升的抑制作用[J]. 中国电机工程学报, 2017, 37(5): 1526-1535. [百度学术]
TONG W M, CHENG X B, SUN J Y,et al.Suppression effect of rotor wafters on permanent magnet temperature rise for high-speed permanent magnet motor[J]. Proceedings of the CSEE,2017, 37(5): 1526-1535.(in Chinese) [百度学术]
鞠宇宁.全封闭永磁同步电机温度场分析及冷却结构设计[D].天津: 天津大学,2018. [百度学术]
JU Y N.Temperature field analysis and cooling structure design of fully enclosed permanent magnet synchronous motor[D].Tianjin:Tianjin University,2018.(in Chinese) [百度学术]
朱高嘉, 朱英浩,佟文明,等.基于流固耦合的高速永磁电机冷却结构分析与改进[J].电工电能新技术,2017,36(12):1-7. [百度学术]
ZHU G J,ZHU Y H,TONG W M,et al.Analysis and enhancement of cooling system of high-speed permanent magnet motor based on computational fluid dynamics[J].Advanced Technology of Electrical Engineering and Energy,2017,36(12):1-7.(in Chinese) [百度学术]
CHEN W,JU Y N,YAN D,et al.Design and optimization of dual-cycled cooling structure for fully-enclosed permanent magnet motor[J].Applied Thermal Engineering,2019,152:338-349. [百度学术]
VEG L,LAKSAR J.Comparison of two types of cooling of axial flux permanent magnet machines by CFD simulation[C]//2019 International Conference on Electrical Drives & Power Electronics (EDPE).The High Tatras,Slovakia.IEEE,2019:303-306. [百度学术]
徐显贺.基于多物理场的高铁永磁同步电机散热设计与性能优化[D].湘潭:湘潭大学.2023 [百度学术]
XU X H .Heat dissipation design and performance optimization of high-speed rail permanent magnet synchronous motor based on multi-physics[D]. Xiangtan:Xiangtan University, 2023. [百度学术]
ACQUAVIVA A,SKOOG S,THIRINGER T.Design and verification of In-slot oil-cooled tooth coil winding PM machine for traction application[J].IEEE Transactions on Industrial Electronics,2021,68(5):3719-3727. [百度学术]
HATZIATHANASSIOU V,XYPTERAS J,ARCHONTOULAKIS G. Electrical-thermal coupled calculation of an asynchronous machine[J].Archiv Für Elektrotechnik,1994,77(2):117-122. [百度学术]
彭德其, 周晓辉,晏才松,等.中型高压电机轴径向混合通风结构的强化散热特性模拟[J].机械科学与技术,2023,42(1):75-82. [百度学术]
PENG D Q,ZHOU X H,YAN C S,et al.Simulation study on heat dissipation enhancement characteristics of axial and radial mixed ventilation structure for medium-sized high-voltage motor[J].Mechanical Science and Technology for Aerospace Engineering,2023,42(1):75-82.(in Chinese) [百度学术]