汽车工程师

参数	数值
工作电压/V	350
峰值功率/kW	70
峰值转矩/N·m	160
最高工作转速/r·min^-1	15 500

参数	数值
工作电压/V	350
峰值功率/kW	70
峰值转矩/N·m	160
最高工作转速/r·min^-1	15 500

转速 /r·min^-1	转矩 /N·m	定子铁损	转子铁损	铜损	磁钢损耗
3 196	35.1	113a	7a	187a	a
5 000	49.4	170a	10.5a	328a	1.5a
10 000	30.7	266.5a	22.5a	284.5a	4a

转速 /r·min^-1	转矩 /N·m	定子铁损	转子铁损	铜损	磁钢损耗
3 196	35.1	113a	7a	187a	a
5 000	49.4	170a	10.5a	328a	1.5a
10 000	30.7	266.5a	22.5a	284.5a	4a

工况	绕组温度/℃
3 196r/min@35 N·m	105.5	100.9
5 000r/min@50 N·m	115.0	119.3
10 000r/min@30 N·m	116.5	119.8

工况	绕组温度/℃
3 196r/min@35 N·m	105.5	100.9
5 000r/min@50 N·m	115.0	119.3
10 000r/min@30 N·m	116.5	119.8

车用永磁同步电机电磁-热耦合分析

PDF下载

魏颖颖 , 王金昊 , 郭守仑 , 王斯博 , 王宇

汽车工程师 | 红旗新能源电驱系统专刊 2023,(7): 32-36

收起

汽车工程师 | 红旗新能源电驱系统专刊 2023, (7): 32-36

车用永磁同步电机电磁-热耦合分析

全屏

魏颖颖, 王金昊, 郭守仑, 王斯博, 王宇

作者信息

中国第一汽车股份有限公司研发总院，长春 130013

Electromagnetic-Thermal Coupled Analysis of Permanent Magnet Synchronous Motor for Vehicles

Yingying Wei, Jinhao Wang, Shoulun Guo, Sibo Wang, Yu Wang

Affiliations

Global R&D Center, China FAW Corporation Limited, Changchun 130013

出版时间: 2023-07-15 doi: 10.20104/j.cnki.1674-6546.20230126

文章导航

摘要

收起

以车用永磁同步电机为研究对象，根据电机定转子冲片的结构尺寸参数和性能参数，在电磁场有限元分析软件中建立二维仿真模型，求解出电机在特殊工况条件下的电磁场及损耗分布情况，将电磁场仿真得到的损耗密度导入温度场仿真软件，计算该工况条件下达到稳态的温度分布情况，由耦合仿真和单物理场仿真对比结果可知，耦合仿真的热点分布更集中，更有利于电机热风险的及时发现和改进。工况温升试验结果表明，耦合仿真温升数据与试验温升数据基本一致，证明了电机耦合仿真计算结果的准确性。

关键词

永磁同步电机 / 电磁-热耦合 / 有限元分析

Abstract

收起

In this paper, a permanent magnet synchronous motor for vehicles was used as the research object. According to the specific structural size parameters and performance parameters of the stator and rotor blades of the motor, a two-dimensional simulation model was established in the electromagnetic field finite element analysis software to solve the electromagnetic field and loss distribution of the motor under special operating conditions. The loss density obtained by electromagnetic field simulation was imported into the temperature field simulation software to calculate the temperature distribution reaching the steady state under the operating conditions. Through the comparison of coupled simulation and single-physics simulation, it is concluded that the hot spots distribution of the coupled simulation is more concentrated, and the thermal risk of the motor can be discovered and improved in time. The working condition temperature rise test of this motor indicates that the coupled simulation temperature rise data is basically consistent with the test data, proving correctness of the motor coupled simulation results.

Key words

Permanent Magnet Synchronous Motor (PMSM) / Electromagnetic and thermal coupled / Finite element analysis

引用本文

魏颖颖, 王金昊, 郭守仑, 王斯博, 王宇. 车用永磁同步电机电磁-热耦合分析. 汽车工程师, 2023 , (7) : 32 -36 . DOI: 10.20104/j.cnki.1674-6546.20230126

Yingying Wei, Jinhao Wang, Shoulun Guo, Sibo Wang, Yu Wang. Electromagnetic-Thermal Coupled Analysis of Permanent Magnet Synchronous Motor for Vehicles[J]. Automotive Engineer, 2023 , (7) : 32 -36 . DOI: 10.20104/j.cnki.1674-6546.20230126

正文

收起

1 前言

收起

永磁同步电机具有功率密度高、效率高、体积小、噪声低等特点，在混合动力汽车领域受到广泛关注。电机运转过程中产生的高温会直接降低电机的工作性能^[1]。对于永磁电机，载荷较大导致的高温可能导致绕组和定子铁芯间的绝缘材料失效或者转子铁芯包裹的永磁体不可逆退磁，使电机功能性损毁^[2-3]，因此，精准分析电机温度场分布非常重要。

电机多物理场耦合仿真计算已成为计算机辅助设计中的关键技术之一^[4-8]。国内外学者采用磁热耦合方法在电机的损耗和温升方面开展了大量研究^[9-12]。本文利用电磁仿真软件建立二维有限元的电机几何模型，计算电机在特殊工况下的损耗密度分布情况，并在温度场仿真软件中建立三维仿真模型，将电磁仿真得到的损耗密度导入三维温度场仿真模型，计算特殊工况下达到稳定运行状态的温度场分布情况，最后，将电磁-热耦合仿真温升数据与试验实测结果进行对比，验证仿真结果的准确性。

2 电磁模型的建立及仿真

收起

本文以永磁同步电机为研究对象，根据电机主要性能参数在电磁仿真软件中建立电磁模型。选取整车常用工况中不同转速下负荷较大的3个工况，根据电机损耗MAP表差值获得不同整车工况的具体参数，计算整车工况的铁损损耗密度分布情况。电机主要性能参数如表1所示，整车工况具体损耗参数如表2所示，其中，将3 196 r/min@35.1 N·m工况下的磁钢损耗设为a，其他损耗由a的倍数表征。

采用JMAG电磁仿真软件计算得到3 196 r/min@35.1 N·m、5 000 r/min@49.4 N·m和10 000 r/min@30.7 N·m 3种工况下的电机铁损损耗密度，其中10 000 r/min@30.7 N·m工况下的电机铁损损耗密度分布如图1所示。

由图1可知，靠近气隙处的定子齿部和转子的损耗密度较大，远离气隙的定子轭部和转子的损耗密度非常小，定子和转子的整体损耗密度由靠近气隙侧向远离气隙侧逐渐减小。

3 损耗密度映射

收起

耦合仿真分为单向耦合和双向耦合。单向耦合是将第1个求解器的求解结果作为条件，输入到第2个求解器，求解数据单向传递，具有计算速度快、节约资源和仿真效率高的优势。双向耦合是将各求解器的求解结果在求解器之间双向传递，反复迭代直到运行结果达到收敛后停止计算。相较于单向耦合，双向耦合计算结果更精确，但计算量显著增大，对计算机要求较高，计算速度低，不利于提高仿真效率，因此本文采用单向耦合的方法对电机的电磁场和温度场进行耦合仿真分析。通过JMAG电磁仿真软件仿真得到损耗结果，以表格形式输出定子铁芯和转子铁芯上的各坐标位置的损耗密度。

将JMAG电磁仿真软件导出的定子铁芯和转子铁芯的损耗密度导入STAR CCM+温度场仿真软件中，利用电磁仿真软件和温度场仿真软件中一致的模型坐标位置，将损耗密度映射到温度场仿真模型上，实现在温度场模型上加载损耗密度的效果。与温度场模型上均匀分布的损耗对比，温度场模型上映射损耗密度更能反映详细的热点位置。在不受其他热源影响的条件下，损耗密度越大的位置温度越高。因此，在其他热源不变的条件下，推测映射损耗密度的耦合仿真得到的温度将会比损耗均匀分布的单物理场仿真得到的温度高，本文对其进行仿真对比验证。

4 电机温度场仿真及分析

收起

根据传热学基本理论，直角坐标系下电机的稳态温度场求解可以归结为如下边值问题^[13]：

（1）$\left\{\begin{array}{l}\frac{\partial }{\partial x}\left({\lambda }_{x}\frac{\partial T}{\partial x}\right)+\frac{\partial }{\partial y}\left({\lambda }_{y}\frac{\partial T}{\partial y}\right)+\frac{\partial }{\partial z}\left({\lambda }_{z}\frac{\partial T}{\partial z}\right)=-{q}_{v}\\ \frac{\partial T}{\partial n}=0\\ -{\lambda }_{z}\frac{\partial T}{\partial n}=h(T-{T}_{e})\end{array}\right.$

式中，λ_x、λ_y、λ_z分别为x、y、z向导热系数；T为温度；q_v为热流密度，即电机各部件单位体积的损耗；n为长度；λ为导热系数；h为对流换热系数；T_e为环境温度。

输入给定损耗与边界条件，可以计算出电机内部温度场分布情况。

本文将一个定子槽中的多层方导线视为一个导线整体，导线与定子槽之间的绝缘材料包括导线漆皮、绝缘纸和灌封的绝缘漆。将绝缘材料等效处理为导线与定子铁芯热传导过程中的热阻。这些绝缘材料的等效导热系数为^[13-14]：

（2）${\lambda }_{eq}=\sum _{i=1}^{n}{\delta }_{i}/\left(\sum _{i=1}^{n}{\delta }_{i}/{\lambda }_{i}\right)$

式中，λ_eq为等效导热系数；δ_i为各绝缘材料的等效厚度；λ_i为各绝缘材料的导热系数。

定转子铁芯是由多层硅钢片叠加形成的，所以径向与轴向导热系数不同。在径向上，定转子铁芯是完整的硅钢材料，导热系数为硅钢材料的导热系数。在轴向上，由于定转子是由多层硅钢片叠加而成的，叠加的硅钢片之间存在空气、灰尘等，导热系数远小于径向导热系数。轴向等效导热系数为^[15]：

（3）${\lambda }_{r}=\frac{{\delta }_{Fe}+{\delta }_{0}}{\frac{{\delta }_{Fe}}{{\lambda }_{1}}+\frac{{\delta }_{0}}{{\lambda }_{0}}}=\frac{1}{\frac{{K}_{Fe}}{{\lambda }_{1}}+\frac{1-{K}_{Fe}}{{\lambda }_{0}}}$

式中，δ_Fe为铁心叠厚；λ₁为硅钢片的导热系数；δ₀为绝缘介质的净长度；λ₀为绝缘介质的导热系数；K_Fe为铁心的叠片系数。

电机机壳和端盖表面与外部环境进行自然对流换热，其换热系数为^[2,14]：

（4）${\alpha }_{c}=14\times {(1+0.5\sqrt{\omega })}^{3}\sqrt{\frac{{T}_{c}}{25}}$

式中，ω为外部环境风速，本文设定为0；T_c为外部环境温度。

转子转动时会引起电机内部空气扰动，从而增强端部绕组与内部空气之间的换热，其对流换热系数可由经验公式^[2,14]获得：

（5）${\alpha }_{s}=\frac{1+0.04{v}_{r}}{0.045}$

式中，v_r为转子外径上的线速度。

转子端部与内部空气之间的对流换热系数的经验公式^[2,14]为：

（6）${\alpha }_{r}=28\times (1+\sqrt{0.45{v}_{r}})$

将电磁仿真得到的电机铁损损耗密度和其他损耗导入三维温度场仿真软件，施加到各零部件上。考虑电机中存在的热对流和热传导，在有接触的零部件之间建立接触热阻和对流换热面，将边界条件设置到建立的接触热阻和对流换热面上。通过温度场仿真计算，得到如图2、图3所示结果。

由图2可知，低转速下，最高温度位于绕组处，高转速下，最高温度位于转子处。根据电磁仿真得到的损耗可知，随着转速升高，铜损逐渐降低，铁损逐渐升高，导致低转速时绕组温度高，高转速时转子温度高。同时，转子无有效的直接散热途径，导致热量集聚在转子内，转子温度升高。

由图2可知，耦合场仿真得到的零部件温度高于单物理场仿真得到的零部件温度，证明了在温度场中映射损耗密度的耦合仿真得到的温度更高的推论。这是因为单物理场中的铁损均匀分布在定转子铁心上，其各坐标位置的损耗密度一致，使其发热均匀。而耦合场中的铁损按损耗密度分布在定转子铁心上，其温度分布受到损耗密度分布的影响，损耗密度越大的位置温度也越高，使耦合仿真得到的热点温度更高。

由图3单物理场与耦合场仿真获得的温度分布情况对比可知，定子和转子的热点位置向气隙侧移动。从图1的电磁仿真结果来看，铁损不是均匀分布的，主要集中分布在气隙附近，温度分布与铁损损耗密度分布对应。与单物理场仿真相比，耦合仿真得到的温度分布热点更为集中，更能反映出电机的散热缺陷，以规避电机运行中的散热风险。

5 试验验证

收起

本文开展试验对仿真结果进行验证。在电机台架上进行工况点稳态温升试验，运行至基本热平衡，在电机绕组端部布置K型温度热电偶进行测温。电机持续运行30 min，每隔1 s记录一次热电偶测温数据，热电偶测得的最高温度如表3所示。

耦合仿真计算结果与试验实测温度非常接近，误差小于4.5%，在可接受范围内。出现误差的原因是试验与仿真的损耗和边界条件存在差异，还需进一步优化。

6 结束语

收起

本文针对某永磁同步电机在特殊工况下进行了电磁-热耦合仿真计算及稳态温升试验，通过对比分析，得到如下结论：

a. 电磁-热耦合仿真计算得到的温度场分布是与电磁仿真得到的损耗密度分布相映射的，热点主要集中分布在气隙附近。相比单物理场仿真中设置的均匀分布的损耗，映射损耗密度的耦合仿真得到了更为集中的热点分布，从而能够精确捕捉热点位置，准确识别电机运行中的散热风险。

b. 电磁-热耦合仿真计算得到的温升数据与试验实测温升数据基本一致，电机耦合仿真计算结果能够作为电机热失效风险预判的重要依据。

下一步计划通过电磁-热-流体耦合仿真方法优化计算结果，进一步缩小仿真结果与试验结果的差距。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

马媛, 李建贵, 刘璟轩. 基于混合动力电动汽车的永磁无刷电机热分析及研究[J]. 工程热物理学报, 2018, 39(11): 2525-2531.

, LI

J G

, LIU

J X

. Thermal Analysis and Research of Permanent Magnet Brushless Motor Based on Hybrid Electric Vehicle Application[J]. Journal of Engineering Thermophysics, 2018, 39(11): 2525-2531.

[2]

兰志勇, 王琳, 徐琛, 等. 高速永磁同步电机温度场及水冷分析[J]. 电气传动, 2018, 48(12): 81-85.

LAN

Z Y

, WANG

, XU

, et al. Analysis of Temperature Field and Water Cooling of High Speed Permanent Magnet Synchronous Motor[J]. Electric Drive, 2018, 48(12): 81-85.

[3]

刘威, 刘慧娟, 杜晋文, 等. 定子水冷永磁同步电机多工况热特性的研究[J]. 微特电机, 2019, 47(4): 1-4+12.

LIU

, LIU

H J

, DU

J W

, et al. Thermal Performance Analysis of Water-Cooled PMSM Running in Multi-Conditions[J]. Small & Special Electrical Machines, 2019, 47(4): 1-4+12.

[4]

路义萍, 李伟力, 马贤好, 等. 大型空冷汽轮发电机转子温度场数值模拟[J]. 中国电机工程学报, 2007, 27(12): 7-13.

Y P

, LI

W L

, MA

X H

, et al. Numerical Simulation of Temperature Field in Rotor of Large Turbo Generator with Air-Coolant[J]. Proceedings of the Chinese Society for Electrical Engineering, 2007, 27(12): 7-13.

[5]

GLIERE

, MASSE

R H

, FAUTRELLE

. A Finite Element-Finite Difference Computation of Magenetic and Turbulent Flow Coupled Problem[J]. IEEE Transactions on Magnetic, 2002, 24(1): 252-255.

[6]

ALBERTI

, BIANCHI

. A Coupled Thermal-Electromagnetic Analysis for a Rapid and Accurate Prediction of IM Performance[J]. IEEE Transactions on Industrial Electronics, 2008, 55(10): 3575-3582.

[7]

COULOMB

J L

. A Methodology for the Determination of Global Electromechanical Quantities from a Finite Element Analysis and Its Application to the Evaluation of Magnetic Forces, Torque and Stiffness[J]. IEEE Transactions on Magnetic, 1983, 19(6): 2514-2519.

[8]

李伟力, 赵志海, 侯云鹏. 大型同步发电机定子同相槽和异相槽的温度场计算[J]. 电工技术学报, 2002, 17(3): 1-6+92.

W L

, ZHAO

Z H

, HOU

Y P

. Calculation of Stator Temperature Field of Generator[J]. Transactions of China Electrotechnical Society, 2002, 17(3): 1-6+92.

[9]

刘奕新, 张志军. 油冷电机磁热耦合仿真研究[J]. 汽车技术, 2022(11): 41-46.

LIU

Y X

, ZHANG

Z J

. Simulative Research on Magneto Thermal Coupling of Oil-Cooled E-Motor[J]. Automobile Technology, 2022(11): 41-46.

[10]

J K

, ZHAO

S Z

, FENG

H C

, et al. Analysis of Temperature Field for a Surface- Mounted and Interior Permanent Magnet Synchronous Motor Adopting Magnetic-Thermal Coupling Method[J]. CES Transactions on Electrical Machines and Systems, 2018, 2(1): 166-174.

[11]

TIKADAR

, KUMAR

, JOSHI

, et al. Coupled Electro-Thermal Analysis of Permanent Magnet Synchronous Motor for Electric Vehicles[C]// 2020 19th IEEE Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems (ITherm). Orlando, FL, USA: IEEE, 2020.

[12]

陈裕明, 陈鹏满, 王柱, 等. KX-ZM01涵道风扇无刷直流电机磁热耦合仿真分析[J]. 机电工程技术, 2022, 51(9): 165-169.

CHEN

Y M

, CHEN

P M

, WANG

, et al. Magnetic Thermal Coupling Simulation Analysis of KX-ZM01 Ducted Fan BLDC Motor[J]. Mechanical & Electrical Engineering Technology, 2022, 51(9): 165-169.

[13]

吴尧辉, 张建强, 刘思. 参数修正的场路结合法仿真异步机起动性能[J]. 计算机仿真, 2017, 34(7): 188-194.

X H

, ZHANG

J Q

, LIU

. Simulation of Starting Performance of Asynchronous Machine Based on Parameter Modified Field-Circuit Combined Method[J]. Computer Simulation, 2017, 34(7): 188-194.

[14]

刘慧军, 陈芬放, 黄瑞, 等. 车用驱动电机冷却系统仿真研究[J]. 中南大学学报(自然科学版), 2020, 51(7): 2002-2012.

LIU

H J

, CHEN

F F

, HUANG

, et al. Simulation Study on Cooling System of Automotive Driving Motor[J]. Journal of Central South University (Science and Technology), 2020, 51(7): 2002-2012.

[15]

陈杰, 王明盛. 高速永磁同步电机冷却系统优化分析[J]. 微特电机, 2020, 48(5): 25-27+31.

CHEN

, WANG

M S

. Optimization Analysis of Cooling System for High-Speed Permanent Magnet Synchronous Motor[J]. Small & Special Electrical Machines, 2020, 48(5): 25-27+31.

2023年第卷第7期

PDF下载

214

引用本文

BibTeX

文章信息

doi: 10.20104/j.cnki.1674-6546.20230126

首发时间：2025-11-25
出版时间：2023-07-15

补充材料

相关文章

文章信息

作者

出版历史

修回日期：2023-06-08

基金

作者信息

中国第一汽车股份有限公司研发总院，长春 130013

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/qcgcs/CN/10.20104/j.cnki.1674-6546.20230126

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

工况	绕组温度/℃
工况	耦合仿真结果	试验实测结果
3 196r/min@35 N·m	105.5	100.9
5 000r/min@50 N·m	115.0	119.3
10 000r/min@30 N·m	116.5	119.8