引用本文
刘慧娟, 刘威, 杜晋文, 宋腾飞. 物流电动车用80 kW永磁同步电机冷却系统分析[J]. 北京交通大学学报, 2019,43(2): 80-85.
LIU Huijuan, LIU Wei, DU Jinwen, SONG Tengfei. Investigating cooling system of an 80 kW PMSM for electric vehicle application[J]. JOURNAL OF BEIJING JIAOTONG UNIVERSITY, 2019,43(2): 80-85.
Doi:10.11860/j.issn.1673-0291.20180099  
Permissions
Copyright©2019, 北京交通大学学报编辑部
北京交通大学学报编辑部 所有
物流电动车用80 kW永磁同步电机冷却系统分析
刘慧娟, 刘威, 杜晋文, 宋腾飞
北京交通大学 电气工程学院,北京 100044

第一作者:刘慧娟(1967—),女,四川资阳人,教授,博士,博士生导师.研究方向为电机设计及其控制,温度场计算.email:hjliu@bjtu.edu.cn.

收稿日期: 2018-07-05
基金: 中央高校基本科研业务费专项资金(2018YJS161); 国家自然科学基金(51377008)
摘要

为了提高电动车用永磁同步电机冷却系统的效率,本文以一台80 kW永磁同步电机为例,利用3D有限元法计算该永磁同步电机在不同水道结构、不同注水方式、以及冷却液不同初始温度和不同流量下的温度场分布,分析研究电机冷却系统各种因素对电机温度场分布的影响,结果表明:相较于螺旋形水道,Z字型水道冷却效果更好,但后者对冷却泵的功率要求也更高;冷却液的初始温度每降低1 ℃,电机的温升平均下降0.9 ℃;单入单出型注水方式比单入双出冷却效果更佳,但是后者对冷却泵的功率要求相对更低.

关键词: 永磁同步电机; 有限元法; 温度场; 冷却系统; 有限体积法
中图分类号:TM351 文献标志码:A 文章编号:1673-0291(2019)02-0080-06
Investigating cooling system of an 80 kW PMSM for electric vehicle application
LIU Huijuan, LIU Wei, DU Jinwen, SONG Tengfei
School of Electrical Engineering, Beijing Jiaotong University, Beijing 100044,China
Fund:Fundemental Research Funds for the Central Universities (2018YJS161); National Natural Science Foundation of China(51377008)
Abstract

In order to improve the efficiency of the permanent magnet synchronous motor cooling system for electric vehicles, this paper takes a 80 kW permanent magnet synchronous motor as an example, and uses the 3D finite element method to calculate the temperature field distribution of the permanent magnet synchronous motor in different waterway structure, different water injection mode, and the different initial temperature and different flow rate of the coolant, and analyzes and studies the motor cooling. The results showed that compared with spiral channel, Z-shaped channel has better cooling effect, but the latter has higher power requirement for cooling pump; the average temperature rise of motor decreases by 0.9 ℃ when the initial temperature of coolant decreases by 1 ℃; and single-in and single-out mode of water injection has better cooling effect than single-in and double-out mode, but the latter has lower power requirement for cooling pump.

Keyword: permanent magnet synchronous motor; finite element method; temperature field; cooling system; finite volume method

随着电动车用永磁同步电机效率和功率密度的提高, 电机的温升问题日益突出[1], 过高的温升对电机的性能以及寿命都是很大的威胁, 为了解决这些问题, 需要从电机的冷却系统入手, 设计出冷却效果更加优异的冷却方案.对于电机冷却系统的研究, 国内外学者提出了很多的方案.目前商业用电机主要通过轴向风冷、机壳水冷、定子油冷以及水冷结合氢冷等方式冷却[1], 而对于电动车用永磁同步电机, 则主要通过水冷的方式来进行[2].

水冷永磁同步电机通常从两个方面优化设计其冷却系统:注水流速和水道结构.就注水方式而言, 单入单出、单入双出、双入双出的注水结构在实际工程中均有应用, 同时文献[2]得到电机温升的减小量与水流速度不成正比, 当水流速度达到一定值时, 电机温升的减小量基本趋于饱和, 过高的流速不仅会额定增加冷却水泵的功耗等, 而且对电机的温升性能改善较小, 从而可能导致整个电机的冷却系统效率降低.对于水道结构, 根据文献[3], 改变水道的长、宽以及内层机壳的厚度对电机散热以及机械强度都有不同程度的影响.单纯针对水道结构而言, 同样的长度、宽度、内层机壳厚度前提下, 轴向Z字形水道冷却效果较螺旋形好[4, 5, 6, 10], 而轴向Z字形水道比较于螺旋形水道对冷却泵的功耗要求更高.

本文作者以80 kW永磁同步电机额定工况下温度场分布为例, 基于Ansys Workbench多物理场仿真平台, 首先对实际的三维电机模型进行几何简化、网格剖分等前处理, 并对绕组等部分进行导热系数等效计算, 再应用Fluent计算电机不同冷却方式下的温度场.比较不同水道结构、不同注水方式、冷却液不同的初始温度和不同流量下电机的温度场分布, 确定出比较优越的冷却方案, 对电机的设计与生产有较大的指导意义.

1 电机模型、参数和网格剖分

考虑到网格前处理的可实现性及对温升的影响程度, 对电机模型的槽绝缘、匝间绝缘、永磁体粘胶等部分进行适当的简化, 建立1∶ 1温度场仿真模型, 电机的模型如图1所示.

图1 电机的三维模型Fig.1 3D model of motor

本文主要是确定电机各部分的最高温升, 核验其材料是否能够耐受极限温度, 转轴、端盖、机座等部件导热系数远大于空气, 实际中更利于散热.同时针对剖分前处理的可实现性, 主要考虑电机的电磁部件, 其基本参数如表1所示.

表1 永磁电机的主要部件基本参数 Tab.1 Basic parameters of permanent magnet motor

将电机模型导入到Workbench中进行网格划分, 模型的网格如图2所示.

图2 电机网格划分Fig.2 Grid division of motor

2 各向异性材料等效导热系数的确定

利用Fluent软件分析电机的温度场, 需要确定电机每部分的材料属性以及各部分的导热系数.各向同性的部位导热系数比较容易处理, 比如机壳; 各向异性的部位需要分别确定 x/y/z3个方向的导热系数, 比如定子绕组和定、转子.电机实体模型中定子槽内有铜导线、漆膜、空气、浸漆、槽绝缘这些结构, 这里将槽内所有的结构等效为一个实体, 轴向与径向的导热系数需要分别确定, 等效过程见图3.

图3 定子槽内部分简化示意图Fig.3 Simplified sketch map in the stator slot

假定槽内空间被导线、漆膜、绝缘层、浸漆与空气均匀填充, 其等效实体结构的径向导热系数恒定不变, 轴向导热系数也恒定不变.

电机绕组径向导热系数计算如下[7]

λj=SSCuλCu+Sqmλqm+Sjyλjy+Sjqλjq+Sairλair(1)

式中: λj、 λcu、 λqm、 λjy、 λjq、 λair分别为槽内等效实体径向、铜、漆膜、槽绝缘、浸漆、空气的导热系数, 单位均为W/ k·m; S、Scu、 Sqm、 Sjy、 Sjq、 Sair分别为等效实体、单个槽内的铜、漆膜、槽绝缘、浸漆、剩余空气的截面积, 单位均为m2.

绕组轴向导热系数的计算为[7]

λz=λCu×Sf(2)

硅钢片的轴向导热系数为[7]

λdg=100100Kcλg+1001-Kcλair(3)

端部铜的导热系数的确定为[7]

λdCu=φCu×λcu(4)

式中: λz、 λg、 λdg、 λdCu分别为轴向等效、硅钢片、轴向等效硅钢片、端部铜的等效导热系数, W/ k·m; Sf为实际槽满率; Kc为硅钢片的叠压系数; φu为端部实际用铜的体积比.

电机沿轴向为 z轴所在的方向, 计算得到电机各个部分的导热系数如表2所示.

表2 电机各部位材料的基本参数 Tab.2 Basic parameters of the material in each part of the motor
3 冷却系统的分析比较

车用永磁同步电机多运行于额定工况, 该工况下为了找到一种效果更好的冷却系统, 本文对不同冷却方式进行了对比分析.保持水道的总长度、截面形状、截面大小以及水道在机壳中的位置不变, 通过改变水道的结构、注水方式、冷却液的初始温度和注水流量, 比较分析额定工况下电机的温度场分布.水道的截面示意图见图4.

图4 水道截面形状示意图Fig.4 Schematic diagram of the channel shape

水道的基本参数见表3.

表3 螺旋水道的基本数据 Tab.3 Basic data of the spiral waterway
3.1 不同形状水道的比较

轴向螺旋形水道与轴向Z字形水道的形状如图5所示, 可以明显看出, 出水口长度明显大于入水口长度, 这样处理是为了防止出水口的回流现象[8, 9].仿真中采用50%乙二醇溶液为冷却液, 注水入口处采用流量入口边界条件, 保持注水的流量为20 L/min, 注水的初始温度为65 ℃; 注水的出口处采用压力出口边界条件.

图5 两种不同形状的水道结构Fig.5 Two different shapes of waterways

螺旋水道与Z字形水道的压强、流速分布和温度分布如图6所示.

图6 螺旋水道与Z字形水道的对比Fig.6 Comparisons of spiral waterways and Z shaped waterways

从图6(a)与6(b)可以看出, 水道的入水口压强较高, 出水口压强低, Z字形水道的最高压强为螺旋水道的4倍多, 所以Z字形水道的水阻大, 对于供水泵的功耗大.从图6(c)与6(d)可以看出, 螺旋形水道中液体的流速分布较为均匀, Z字形水道液体的流速分布不均匀, 流速较大的区域只占据水道截面的一半, 其他区域液体流速很小, 对应于图6(e)与6(f)水道的温度分布, 可以看出螺旋形水道随着水的流动, 温度逐渐升高, 温度分布均匀, Z字形水道在流速较大的位置液体的温度较低, 流速慢的位置温度高.具体的数据对比如表4所示.

表4 不同水道形状对应的温升与压差 Tab.4 Temperature rise and pressure differences corresponding to the different waterways shapes

从表4可以看出Z字形水道进出水口的压差约为螺旋形水道的4.5倍, 电机采用螺旋形水道, 温升较大, 所以Z字形水道的整体冷却效果稍好于螺旋形水道, 但是Z字形水道对于供水泵的功耗远大于螺旋形水道, 综合经济型等因素考虑, 文中的永磁同步电机采用螺旋形水道.

3.2 不同注水方式的比较

针对螺旋形水道, 保持总的注水流量为20 L/min不变, 注水的初始温度为65 ℃, 通过改变不同的注水方式, 比较电机的温升以及压强分布情况, 注水方式如图7所示.

图7 水道的注水示意图Fig.7 Water injection diagram of a waterway

表5为3种注水方式下, 水道进出水口的压差与电机的温升情况.可以看出单入单出型注水方式压差差不多是其他两种的3倍, 单入双出型压差最小, 对水泵的功耗最小.使用单入单出型注水方式电机的温升最小, 冷却效果最佳.双入单出型注水方式电机的温升最高, 冷却效果最差.本电机采用单入单出型注水方式.

表5 不同注水方式下压差与温升 Tab.5 Pressure differences and temperature rise under different water injection modes
3.3 不同注水流量的比较

根据经验可以确定单位时间内注入冷却液的流量越大, 电机的冷却效果越好.实际中考虑到如果注入冷却液的速度过快, 对于电机机壳的机械强度是一个考验, 更重要的是越高的注入速度消耗的水泵的功率就越高, 不利于节能; 再者, 根据文献[2], 当流速达到一定值时, 再提高冷却液的流量对于电机的温升影响不大, 综合考虑节能以及冷却效果等因素, 选取冷却效果曲线的拐点处所对应的流速最佳.为了找到相对经济的注水流量, 等距选取224 L/min之间10个注水流量值进行温度场仿真, 压强以及温升随注水流量的变化情况如表6所示, 变化曲线如图8所示.

表6 不同流量下电机温升以及压强的比较 Tab.6 Comparison of the temperature rise and pressure of the motor under different flow rates

图8 电机温升以及压差随注入速量的变化曲线Fig.8 Variation curves of temperature rise and pressure difference of motor with injection rate

从图8中的变化曲线可以看出, 随着注水流量的加大, 最开始电机温升下降的比较明显, 但是在流量大于12 L/min以后, 电机温升下降的相对缓慢, 所以12 L/min是比较经济的注水流量.压差与注水流量基本呈线性变化, 流量越大, 流体的压差越大.为了达到充分冷却的效果, 本电机使用20 L/min的注水流量.

3.4 不同注水初始温度的比较

为了比较不同初始温度的冷却液对于电机温升的影响, 设置冷却液的初始温度温度梯度为5 ℃, 通过仿真得到电机的温升随冷却液不同初始温度的变化情况如图9所示.

图9 电机温升随冷却液初始温度的变化Fig.9 Variations of motor temperature rise with initial temperature of coolant

从图9可以看出, 电机的温升与冷却液初始温度的变化曲线基本呈线性关系, 冷却液的初始温度越低, 电机的温升越小, 冷却效果越好.冷却液的初始温度每降低1 ℃, 电机的温升平均降低0.9 ℃, 冷却液的初始温度对于电机的冷却效果影响较大.为了模拟极端的工作环境, 本电机冷却液的初始温度采用65 ℃.

4 结论

本文以80 kW永磁同步电机为研究对象, 利用3D有限元法仿真计算了不同水道结构、不同注水方式、冷却液不同初始温度和不同流量下电机的温度场分布, 得出以下结论:

1)Z字形水道的冷却效果稍好于螺旋形水道, 但是Z字形水道的水阻约为螺旋形水道水阻的4.5倍, 对冷却泵等的功耗要求更高.

2)单入单出型注水方式的冷却效果最好, 单入双出型的水阻最小, 单入单出型的水阻约为单入双出型的3倍.

3)冷却液的初始温度对于电机的温升影响很大, 冷却液的温度每降低1 ℃, 电机的温升平均降低0.9 ℃, 12 L/min是最经济的注水流量.

4)对于本电机或者相似电机, 选用螺旋水道, 单入单出型注水方式, 注水流速选用12 L/min, 在条件许可的情况下尽量降低冷却液的初始温度时, 对于电机的综合冷却效果最佳.

参考文献
[1] 何思源. 全封闭永磁同步牵引电动机冷却系统设计[J]. 大功率变流技术, 2012(3): 48-53.
HE Siyuan. Cooling system design of fully enclosed permanent magnet synchronous traction motor[J]. High-Power Converter Technology, 2012(3): 48-53. (in Chinese) [本文引用:2]
[2] WANG W, YANG F, LUN S. The cooling system design and performance analysis of the disc type permanent magnet synchronous motor[C]//Proceedings of IEEE International Conference on Vehicular Electronics and Safety. Dongguan, 2013: 202-205. [本文引用:3]
[3] ZHANG B, QU R, FAN X, et al. Thermal and mechanical optimization of water jacket of permanent magnet synchronous machines for EV application[C]//IEEE International Electric Machines & Drives Conference . Coeur d’Alene, 2015: 1329-1335. [本文引用:1]
[4] 佟文明, 程雪斌. 高速水冷永磁电机冷却系统分析[J]. 电机与控制应用, 2016(3): 16-21.
TONG Wenmin, CHENG Xuebin. Analysis on cooling system of high speed water-cooled permanent magnet motor[J]. Motor and Control Application, 2016(3): 16-21. (in Chinese) [本文引用:1]
[5] 胡田, 唐任远, 李岩, 等. 永磁风力发电机三维温度场计算及分析[J]. 电工技术学报, 2013(3): 122-126.
HU Tian, TANG Renyuan, LI Yan, et al. Calculation and analysis of three-dimensional temperature field of permanent magnet wind turbines[J]. Journal of Electrical Technology , 2013(3): 122-126. (in Chinese) [本文引用:1]
[6] 刘蕾. 纯电动汽车水冷永磁同步电机多工况热特性及冷却系统研究[D]. 合肥: 合肥工业大学, 2015.
LIU Lei. Study on thermal characteristics and cooling system of water cooled permanent magnet synchronous motor in multiple operating conditions for pure electric vehicles[D]. Heifei: Hefei University of Technology, 2015. (in Chinese) [本文引用:1]
[7] A N 鲍里先科, B T丹科, A N 亚科夫列夫. 电机中的空气动力学与热传递[M]. 魏书慈, 邱建甫译. 北京: 机械工业出版社, 1985.
BORISHCHENKO A N, DANKO B, YAKOVLEV A N. Aerodynamics and heat transfer in electric machinery[M]. WEI Shuci, QIU Jianpu, trans. Beijing: Machinery Industry Press, 1985. (in Chinese) [本文引用:4]
[8] 王福军. 计算流体动力学分析——CFD软件原理与应用[M]. 北京: 清华大学出版社, 2004.
WANG Fujun. Computational fluid dynamics analysis - principles and applications of CFD software[M]. Beijing: Tsinghua University Press, 2004. (in Chinese) [本文引用:1]
[9] 贾月梅, 赵秋霞, 赵广慧. 流体力学[M]. 北京: 国防工业出版社, 2006.
JIA Yuemei, ZHAO Qiuxia, ZHAO Guanghui. Fluid mechanics[M]. Beijing: National Defense Industry Press, 2006. (in Chinese) [本文引用:1]
[10] 杨学威, 张小发. 电机壳体Z字型冷却水道设计[J]. 电机与控制应用, 2016, 43(9): 62-65.
YANG Xuewei, ZHANG Xiaofa. Design of Z-type cooling channel for motor housing[J]. Application of Motor and Control, 2016, 43(9): 62-65. (in Chinese) [本文引用:1]