汽车技术

参数	参数值	参数	参数值
定子槽数	48	额定功率/kW	60
极数	8	额定转速/r·min^-1	5 000
定子外径/mm	220	额定转矩/N·m	112
定子内径/mm	154	线径/mm	0.92
电机轴长/mm	140	铜径/mm	0.83
气隙/mm	0.8	并联支路数	8

参数	参数值	参数	参数值
定子槽数	48	额定功率/kW	60
极数	8	额定转速/r·min^-1	5 000
定子外径/mm	220	额定转矩/N·m	112
定子内径/mm	154	线径/mm	0.92
电机轴长/mm	140	铜径/mm	0.83
气隙/mm	0.8	并联支路数	8

因子水平	优化参数
水平1	8.4	0.65	2.5	23.0	2.5	13.0
水平2	8.7	0.70	3.0	23.5	3.0	13.5
水平3	9.0	0.75	3.5	24.0	3.5	14.0
水平4	9.3	0.80	4.0	24.5	4.0	14.5
水平5	9.6	0.85	4.5	25.0	4.5	15.0

因子水平	优化参数
水平1	8.4	0.65	2.5	23.0	2.5	13.0
水平2	8.7	0.70	3.0	23.5	3.0	13.5
水平3	9.0	0.75	3.5	24.0	3.5	14.0
水平4	9.3	0.80	4.0	24.5	4.0	14.5
水平5	9.6	0.85	4.5	25.0	4.5	15.0

试验次序	因子水平	P_Fe/W	T_avg/N·m	K_mb/%
1	1	1	1	1	1	1	840.5	294.9	1.88
2	1	2	2	2	2	2	847.1	277.2	1.15
3	1	3	3	3	3	3	815.6	284.8	1.94
4	1	4	4	4	4	4	856.7	290.2	2.03
5	1	5	5	5	5	5	826.6	296.9	2.14
6	2	1	2	3	4	5	929.2	290.4	1.71
7	2	2	3	4	5	1	886.2	277.3	3.55
8	2	3	4	5	1	2	821.1	286.5	2.03
9	2	4	5	1	2	3	835.5	281.9	2.09
10	2	5	1	2	3	4	854.5	270.7	1.91
11	3	1	3	5	2	4	849.2	292.6	1.93
12	3	2	4	1	3	5	832.4	291.2	1.88
13	3	3	5	2	4	1	975.3	281.1	2.24
14	3	4	1	3	5	2	863.1	272.1	1.87
15	3	5	2	4	1	3	785.8	271.6	1.93
16	4	1	4	2	5	3	943.2	290.9	2.03
17	4	2	5	3	1	4	903.6	290.9	2.03
18	4	3	1	4	2	5	836.2	277.9	1.64
19	4	4	2	5	3	1	863.2	272.6	1.95
20	4	5	3	1	4	2	854.6	273.9	2.03
21	5	1	5	4	3	2	842.4	294.7	2.01
22	5	2	1	5	4	3	901.4	279.6	1.69
23	5	3	2	1	5	4	901.9	281.7	1.86
24	5	4	3	2	1	5	820.2	278.9	1.82
25	5	5	4	3	2	1	835.5	272.9	2.18

试验次序	因子水平	P_Fe/W	T_avg/N·m	K_mb/%
1	1	1	1	1	1	1	840.5	294.9	1.88
2	1	2	2	2	2	2	847.1	277.2	1.15
3	1	3	3	3	3	3	815.6	284.8	1.94
4	1	4	4	4	4	4	856.7	290.2	2.03
5	1	5	5	5	5	5	826.6	296.9	2.14
6	2	1	2	3	4	5	929.2	290.4	1.71
7	2	2	3	4	5	1	886.2	277.3	3.55
8	2	3	4	5	1	2	821.1	286.5	2.03
9	2	4	5	1	2	3	835.5	281.9	2.09
10	2	5	1	2	3	4	854.5	270.7	1.91
11	3	1	3	5	2	4	849.2	292.6	1.93
12	3	2	4	1	3	5	832.4	291.2	1.88
13	3	3	5	2	4	1	975.3	281.1	2.24
14	3	4	1	3	5	2	863.1	272.1	1.87
15	3	5	2	4	1	3	785.8	271.6	1.93
16	4	1	4	2	5	3	943.2	290.9	2.03
17	4	2	5	3	1	4	903.6	290.9	2.03
18	4	3	1	4	2	5	836.2	277.9	1.64
19	4	4	2	5	3	1	863.2	272.6	1.95
20	4	5	3	1	4	2	854.6	273.9	2.03
21	5	1	5	4	3	2	842.4	294.7	2.01
22	5	2	1	5	4	3	901.4	279.6	1.69
23	5	3	2	1	5	4	901.9	281.7	1.86
24	5	4	3	2	1	5	820.2	278.9	1.82
25	5	5	4	3	2	1	835.5	272.9	2.18

参数因子	平均转矩方差	转矩脉动系数方差	铁损方差
r	240.80	80.27	0.021 5	15.6	189.46	10.2
δ	28.66	9.55	0.003 5	2.5	340.92	18.4
h_m1	16.95	5.65	0.022 4	16.3	106.27	5.8
ω_m1	4.35	1.45	0.013 5	9.8	265.05	14.3
h_m2	1.89	0.63	0.038 4	27.9	767.58	41.4
ω_m2	7.35	2.45	0.038 4	27.9	183.48	9.9
总	300	100	0.137 7	100	1 852.76	100

参数因子	平均转矩方差	转矩脉动系数方差	铁损方差
r	240.80	80.27	0.021 5	15.6	189.46	10.2
δ	28.66	9.55	0.003 5	2.5	340.92	18.4
h_m1	16.95	5.65	0.022 4	16.3	106.27	5.8
ω_m1	4.35	1.45	0.013 5	9.8	265.05	14.3
h_m2	1.89	0.63	0.038 4	27.9	767.58	41.4
ω_m2	7.35	2.45	0.038 4	27.9	183.48	9.9
总	300	100	0.137 7	100	1 852.76	100

转速 /r·min^-1	P_Fe/W	T_avg/N·m	K_mb/%
2 000	811.7	777.3	284.2	280.0	2.04	1.80
15 000	5 151.2	4 647.8	290.3	276.6	2.10	1.30

转速 /r·min^-1	P_Fe/W	T_avg/N·m	K_mb/%
2 000	811.7	777.3	284.2	280.0	2.04	1.80
15 000	5 151.2	4 647.8	290.3	276.6	2.10	1.30

转速 /r·min^-1	计算方法	磁滞损耗 /W	涡流损耗 /W	附加损耗 /W	总铁损 /W
3 000	方法1	321.9	716.1	103.2	1 141.2
方法2	357.7	795.8	114.0	1 267.5
方法3	335.8	747.1	107.1	1 189.6
5 000	方法1	519.2	662.8	206.6	1 388.4
方法2	642.6	834.4	283.2	1 661.2
方法3	610.8	809.2	290.4	1 610.4
9 000	方法1	669.6	874.9	240.2	1 784.7
方法2	846.1	1 107.5	342.9	2 296.5
方法3	769.5	1 056.1	280.0	2 105.6
12 000	方法1	1 010.5	1 333.3	390.4	2 734.2
方法2	1 339.5	1 767.6	475.1	3 582.2
方法3	1 244.6	1 584.7	487.1	3 316.4
15 000	方法1	1 558.8	2 022.4	575.6	4 156.8
方法2	1 996.2	2 482.8	570.8	5 049.7
方法3	1 870.2	2 403.1	576.2	4 849.5

转速 /r·min^-1	计算方法	磁滞损耗 /W	涡流损耗 /W	附加损耗 /W	总铁损 /W
3 000	方法1	321.9	716.1	103.2	1 141.2
方法2	357.7	795.8	114.0	1 267.5
方法3	335.8	747.1	107.1	1 189.6
5 000	方法1	519.2	662.8	206.6	1 388.4
方法2	642.6	834.4	283.2	1 661.2
方法3	610.8	809.2	290.4	1 610.4
9 000	方法1	669.6	874.9	240.2	1 784.7
方法2	846.1	1 107.5	342.9	2 296.5
方法3	769.5	1 056.1	280.0	2 105.6
12 000	方法1	1 010.5	1 333.3	390.4	2 734.2
方法2	1 339.5	1 767.6	475.1	3 582.2
方法3	1 244.6	1 584.7	487.1	3 316.4
15 000	方法1	1 558.8	2 022.4	575.6	4 156.8
方法2	1 996.2	2 482.8	570.8	5 049.7
方法3	1 870.2	2 403.1	576.2	4 849.5

转速/r·min^-1	∑p/W	P_Cu/W	p_FeR/W	P_fw/W	P_∆/W	p_Fe/W
3 003	3 667.2	1 650.2	48.3	660.1	66.3	1 242.3
5 000.2	4 538.7	1 852.4	53.9	824.3	121.4	1 686.7
9 000.5	5 940.6	2 299.7	47.2	892.9	224.6	2 476.2
12 030	7 995.6	2 864.1	39.8	1 186.2	305.3	3 600.2
14 998	11 803.7	4 395.6	33.4	1 818.7	389.4	5 166.6

转速/r·min^-1	∑p/W	P_Cu/W	p_FeR/W	P_fw/W	P_∆/W	p_Fe/W
3 003	3 667.2	1 650.2	48.3	660.1	66.3	1 242.3
5 000.2	4 538.7	1 852.4	53.9	824.3	121.4	1 686.7
9 000.5	5 940.6	2 299.7	47.2	892.9	224.6	2 476.2
12 030	7 995.6	2 864.1	39.8	1 186.2	305.3	3 600.2
14 998	11 803.7	4 395.6	33.4	1 818.7	389.4	5 166.6

转速 /r·min^-1	电机效率/%
3 000	95.52	95.32	94.37	94.82
5 000	96.78	96.26	95.67	95.84
9 000	96.85	96.54	95.86	96.12
12 000	95.68	95.62	93.94	95.43
15 000	95.24	95.12	93.12	94.60

转速 /r·min^-1	电机效率/%
3 000	95.52	95.32	94.37	94.82
5 000	96.78	96.26	95.67	95.84
9 000	96.85	96.54	95.86	96.12
12 000	95.68	95.62	93.94	95.43
15 000	95.24	95.12	93.12	94.60

电动汽车用高速永磁同步电机的优化设计与铁损分析^*

PDF下载

赵南南 ¹ , 许檬 ¹ , 杨延忠 ¹ , 韩常青 ¹ , 陈阁 ²

汽车技术 | 2023,(1): 21-28

收起

汽车技术 | 2023, (1): 21-28

电动汽车用高速永磁同步电机的优化设计与铁损分析^*

全屏

赵南南¹, 许檬¹, 杨延忠¹, 韩常青¹, 陈阁²

作者信息

¹ 西安建筑科技大学，西安 710055

² 三门峡速达交通节能科技股份有限公司，三门峡 472000

通讯作者:

许檬（1997—），硕士研究生，研究方向为新型电力传动系统、新能源汽车驱动系统与控制， xumeng@xauat.edu.cn。

Optimization Design and Iron Loss Analysis of High-Speed Permanent Synchronous Motor for Electric Vehicles

Nannan Zhao¹, Meng Xu¹, Yanzhong Yang¹, Changqing Han¹, Ge Chen²

Affiliations

¹ Xi’an University of Architecture and Technology, Xi’an 710055

² Sanmenxia Suda Transportation Energy Saving Technology Co., Ltd., Sanmenxia 472000

出版时间: 2023-01-24 doi: 10.19620/j.cnki.1000-3703.20220019

文章导航

摘要

收起

为了提升电动汽车用高速永磁同步电机电磁性能及效率，应用田口法对某额定功率60 kW、额定转速5 000 r/min的电动汽车用高速永磁同步电机的结构进行了优化，通过分析各优化变量对电机铁损、平均电磁转矩和转矩脉动的影响程度，得到了最佳优化方案；建立电机的铁损模型，分别用正弦等效法、谐波分量法和有限元法计算了电机的定子铁损；最后，通过对比定子铁损和效率的测量结果与计算结果，验证了电机铁损分析方法的合理性。

关键词

高速电机 / 永磁同步电机 / 田口法 / 铁损

Abstract

收起

The Taguchi method is adopted to optimize the structure of high-speed permanent magnet synchronous motor for an electric vehicle with rated speed of 5000 r/min and rated power of 60 kW, in an attempt to improve the electromagnetic performance and efficiency of high-speed permanent magnet synchronous motor for electric vehicles. By analyzing the effect of the optimization variables on the motor iron loss, the average electromagnetic torque and the torque ripple, this paper obtaines the optimal optimization solution. The iron loss model of motor is established, and the Sine-equivalence method, harmonic-component method and FEA method are used to calculate the iron loss of stator respectively. Finally, the rationality of iron loss analysis method of the motor is verified by comparing calculated results and measured results.

Key words

High-speed motor / Permanent magnet synchronous motor / Taguchi method / Iron loss

引用本文

赵南南, 许檬, 杨延忠, 韩常青, 陈阁. 电动汽车用高速永磁同步电机的优化设计与铁损分析^*. 汽车技术, 2023 , (1) : 21 -28 . DOI: 10.19620/j.cnki.1000-3703.20220019

Nannan Zhao, Meng Xu, Yanzhong Yang, Changqing Han, Ge Chen. Optimization Design and Iron Loss Analysis of High-Speed Permanent Synchronous Motor for Electric Vehicles[J]. Automobile Technology, 2023 , (1) : 21 -28 . DOI: 10.19620/j.cnki.1000-3703.20220019

正文

收起

1 前言

收起

目前，永磁同步电机是纯电动汽车使用的主要驱动电机，其转速的提高有助于提高电机的效率和功率密度。纯电动汽车用永磁同步电机的转速通常高于10 000 r/min，最高可达20 000 r/min。在对高速永磁同步电机进行优化时，采用粒子群算法等优化算法可以对全局参数进行自动寻优，但电机的优化变量较多，计算耗时较长。田口法因其能以较少的试验次数得到最优的设计变量组合，可快速实现电机结构多目标设计，在电机优化设计上具有优势。

田口法是一种通过建立正交表选择试验参数并安排试验的局部优化方法^[1]。文献[2]对永磁直线电机粒子群优化和田口优化的试验结果进行对比，得出田口法的优化结果更精确、更高效的结论；文献[3]采用田口法和响应面法对永磁同步电机进行优化，优化结果表明，电机在低速（2 700~4 000 r/min）状态下效率的误差小于1%；文献[4]采用田口法改进转子几何结构，减小了定子铁损。目前，田口法多用于低速状态下电机结构的优化，对电机全工况下的研究较少。

在较高频率下，高速永磁同步电机的磁密和损耗分布较普通电机更为复杂，损耗系数呈指数趋势增长^[5]，导致铜损、定子铁损大幅增加。因此，研究高速永磁同步电机损耗非常必要^[6]。

电机损耗主要包括铁心损耗和铜损，此外还有机械损耗、转子涡流损耗和其他附加损耗。文献[7]分析了异步电机空载运行时定、转子铁心中谐波分量对铁损分布的影响；文献[8]采用矩形磁铁和三角形屏障的转子，使转子谐波损耗减小，优化后的电机铁损下降了一半；文献[9]将磁密分量进行傅里叶分析以提高铁损分析的精确性。

本文针对某额定功率60 kW、额定转速为5 000 r/min的纯电动汽车用高速永磁同步电机，以电机平均铁损最小、平均转矩最大、转矩脉动最小、永磁体用量较少为优化目标，利用田口法对电机转子结构相关参数进行优化，对电机定子齿部、轭部进行细分，分析各区域典型位置的磁密变化规律，并对谐波进行分析，分别在电机低速、中速、高速状态下计算定子各区域的铁损和电机的总铁损，最后对样机开展试验，并对正弦法、谐波分析法和有限元法计算的铁损结果及效率进行对比。

2 电机分析与优化

收起

2.1 样机主要参数与有限元分析

图1所示为初步设计的纯电动汽车用高速永磁同步电机的1/8模型，定子采用48槽，永磁体采用内置“V一”型结构，极对数为4。表1所示为电机的主要规格参数。高速永磁同步电机高频状态下，在定子绕组中产生趋肤效应和临近效应，会造成较大的附加损耗，所以高速永磁同步电机定子绕组采用较细的导线并联绕制^[9]。

应用ANSYS Maxwell分析电机模型，在电机剖分时，设置定、转子网格大小均为4 mm，永磁体、绕组、气隙网格大小分别为2 mm、3 mm、2 mm，如图2所示。图3所示为电机额定负载下的磁密云图，图4所示为电机气隙磁密波形，气隙磁密幅值为0.84 T，所设计的电机磁力线分布合理，永磁体利用充分。

2.2 电机优化

电机体积受车内空间的限制，所以本文在电机体积、质量、功率密度一定的条件下对电机的参数进行优化，优化目标为永磁体的体积尽量小、电机铁损最小和转矩脉动最小。由于气隙长度、永磁体在转子中的位置和永磁体长度、宽度对优化目标影响较大，故本文选取第1层相邻永磁体的隔磁桥间距r，第1层永磁体的厚度h_m1、宽度ω_m1，第2层永磁体的厚度h_m2、宽度ω_m2和气隙宽度δ共6个参数进行优化，如图5所示。

本文针对每个优化参数选取5个因子水平，以电机各参数的初始设计值为中心，取4个相邻值建立正交表，如表2所示。建立试验矩阵L₂₅(5⁶)，选取的优化对象为铁损P_Fe、转矩脉动系数K_mb、平均转矩T_avg。表3所示为2 000 r/min转速下根据正交矩阵对电机进行有限元计算的结果。

K_mb可用于评价电机总体脉动程度^[10]：

（1）

K m b = T m a x - T m i n 2 T a v g × 100 %

式中，T_max、T_min分别为电机瞬态仿真下的最大、最小转矩。

2.2.1 计算平均值

计算可得2 000 r/min转速下各次试验电机铁损、平均转矩、转矩脉动系数的平均值分别为860.8 W、282.9 N·m、1.98%。然后计算出各水平下各性能指标的平均值。

2.2.2 方差分析及优化方案的确定

通过分析各优化因子在不同水平下某一性能的平均值对该性能的全体平均值的方差，可以评估各因子变化对该性能影响的程度。方差表达式为：

（2）

S A = 15 ∑ i = 1 5 m A (i) - m 2

式中，S_A为因子A各水平下某一优化对象的方差；m_A(_i₎为因子A在水平i下某一优化对象的值；m为某一优化目标全体试验平均值。

在2 000 r/min转速条件下计算各参数设置下定子铁损、平均电磁转矩和转矩脉动的方差，结果如表4所示。

2.2.3 电机最终优化方案的确定

由表4可以看出：r、δ对平均转矩的影响占比较大；h_m2对铁损的影响占比最大；h_m1、ω_m1、ω_m2对转矩脉动的影响占比较大。综合各优化参数对平均铁损、平均转矩、转矩脉动系数的影响比重，最终方案采用各参数的平均组合：r(1)、δ(5)、h_m1(3)、ω_m1(5)、h_m2(2)、ω_m2(2)。

表5所示为电机优化前、后性能指标对比结果，电机2 000 r/min转速下，铁损、平均转矩、转矩脉动分别降低了4.23%、1.5%、11.8%；15 000 r/min转速下，铁损、平均转矩、转矩脉动分别降低了9.8%、4.6%、38.1%。同时，永磁体体积减小了15.8%。可以看出，高速15 000 r/min时和2 000 r/min时的优化结果趋势相同。

铁损的形成机理较为复杂，与高频材料的损耗特性和电机定子不同位置磁密变化有关，需针对电机定子铁损进行精确计算。

3 电机定子铁损计算

收起

对于大功率、转速较高的永磁同步电机，定子铁损在损耗中占比最大，需要精确的计算模型。电机铁损主要由磁滞损耗P_h、涡流损耗P_c和附加损耗P_e构成。高速永磁同步电机铁心材料中包含旋转磁化和交变磁化两种磁化方式。旋转磁化分为圆形旋转磁化和椭圆形旋转磁化。圆形旋转磁化的磁场大小不变，磁化方向随时间而旋转；椭圆形旋转磁场的磁场大小和磁化方向均随时间变化。交变磁化的磁化方向不变，磁场大小按正弦规律变化^[11]。

3.1 定子铁损模型

为了使定子铁损计算结果更加精确，目前应用较多的铁损计算模型是文献[12]提出的Bertotti铁损模型：

（3）P_Fe=P_h+P_c+P_e=k_hf

B n α

+k_c(fB_n)²+k_e(fB_n)^3/2

式中，P_Fe为定子铁损；B_n为铁心磁通密度幅值；f为频率；k_h、α为磁滞损耗系数；k_c为涡流损耗系数；k_e为附加损耗系数。

电机铁心中磁密波形不一定为正弦，存在一系列谐波，考虑铁心中基波和各次谐波的影响，可对磁密波形进行傅里叶分析，带有谐波分量的磁密可等效为一系列椭圆形谐波磁密矢量，采用2个相互正交的交变磁密矢量B_r(t)、B_t(t)来等效椭圆形旋转磁化^[10]：

（4）B(t)=B_r(t)+B_t(t)

z次谐波的椭圆形旋转磁场可以等效为径向磁密矢量为B_r_z(t)和切向磁密矢量为B_t_z(t)的2个正交的磁密矢量，分别计算这2个磁密引起的铁损并相加即可得到电机的铁损。铁心损耗模型的表达式为：

（5）

P F e s = k h ∑ z N z f B r z x + B t z x + k c ∑ z N z f 2 B r z 2 + B t z 2 + k e ∑ z N z f 1.5 B r z 1.5 + B t z 1.5

式中，x为磁滞损耗系数；B_r_z、B_t_z分别为第z次谐波磁密基波的径向分量和切向分量；N为谐波次数。

3.2 高频损耗系数

电机定子需选取高频、高速时单位铁损较低的硅钢片材料或软磁复合材料，以降低定子铁损^[11]。本文铁心材料采用30DH，采用最小二乘法对该材料在多个频率下的磁密-铁损曲线进行拟合，表达式为^[13]：

（6）

f k h, k c, k e = ∑ i = 1 ∑ j = 1 P i j - k h f i B m i j 2 + k c f i 2 B m i j 2 + k e f i 1.5 B m i j 1.5 2 = k m i n

式中，P_ij为受i、j两个参数影响的铁损；f_i为第i次频率；B_m_ij为受i、j两个参数影响的磁密幅值；k_min损耗系数最小值。

在计算15 000 r/min转速下铁损的损耗系数时，对式（6）中的各铁损系数（k_h、k_c、k_e）分别进行偏导计算，得到各铁损系数的三元一次方程，求解方程组得到铁损的3个系数。

3.3 定子典型位置磁密波形与谐波分析

由式（3）~式（5）可知，精确计算铁损必须分析铁心中各区域的磁密分布情况，得出不同铁损模型对应的铁心磁密幅值。本文将定子齿部、轭部分为5个区域，分别选取点A、B、C、D、E，如图6所示，定子铁损为各区域铁损之和。

3.3.1 磁密波形分析

图7所示为正弦激励下定子不同区域中A~E点的磁密波形图，由图7可以看出，A、B、C点的磁密幅值变化较大，在1.1 T以上，D、E点的磁密幅值变化较小。

3.3.2 考虑谐波和旋转磁场的磁密波形分析

高速永磁同步电机运行时，由于定、转子之间匝链的磁链是不断旋转变化的，因此电机每个剖分单元的磁密都是不断变化的矢量。为了准确计算定子铁损，需求解各点的磁密分布，并将磁密等效为2个互相正交的磁密。在电机高速状态下（12 000 r/min），定子铁心5个不同区域（见图6）磁密径向分量B_r和磁密切向分量B_t随转子位置变化的波形曲线和磁密轨迹曲线如图8~图12所示。

由图8可知，定子齿顶部分（A点）的磁密分量主要是径向分布，切向磁密分量较少。A点磁通密度在径向和切向均有变化，主要受交变磁场的影响，旋转磁场对齿顶的磁化方式影响较小。

图9中定子齿身（B点）的磁密轨迹近似为直线，表明其磁密轨迹主要在径向变化，受交变磁场影响较大，且几乎不受旋转磁场的影响。所以，定子齿身的磁化方式主要是交变磁化。

图10中定子齿根（C点）的磁密轨迹为近似的椭圆形，切向的磁密分量变化较多，因此定子齿根主要受旋转磁场的影响。

图11、图12表明定子轭部（D、E点）磁密轨迹主要是切向分布，D点切向磁密分量远大于径向磁密，E点的磁密分量主要于切向变化，说明定子铁心轭部主要受交变磁场影响。

由图8~图12可知，A~E点的径向磁密与切向磁密在一个周期内呈非正弦分布。因此，各区域中的径向磁密、切向磁密含有大量的谐波分量，需要对其进行谐波分解。经过傅里叶分析，各点的径向磁密和切向磁密的基波及各次谐波的幅值如图13所示。

图13中：A点由于受开槽影响，磁密径向、切向分量变化明显，径向、切向谐波含量均较多；B点由于其磁力线主要在径向变化，傅里叶分析后径向谐波含量较多，切向谐波分量较少；C点和D点的磁力线在切向和径向均有变化，径向磁密含量较为明显；相比于C点，D点的切向磁密波形更接近正弦波，切向磁密谐波含量较少；E点由于其磁力线主要在切向变化且接近正弦波，切向磁密基波幅值较为明显，切向和径向磁密谐波均较少。

3.4 不同转速下的铁损计算

本文采用3种方法计算高速永磁同步电机在额定负载下的定子铁损^[14]。方法1为正弦等效法，将图7中各点磁密幅值B_m代入式（3）中计算定子铁损。方法2为谐波分析法，将磁密分解为相互正交的2个磁密（见图8~图12），分解得到各次谐波的B_r、B_t（见图13），将B_r、B_t分别带入式（5）计算其引起的铁损，相加即可求得电机的定子铁损。方法3为利用有限元仿真计算的定子铁损。如表6所示为不同转速下定子铁损的计算结果。

4 样机试验

收起

为了验证铁损模型计算的精确度，本文以一台样机为损耗试验测试对象，试验台架如图14所示，主要由测功机、样机、示波器、扭矩仪等组成。采用致远电子PA5000H测功机，测功机采集电机的输入、输出功率，精度为0.1%；测功机上的扭矩仪测量电机的转矩和转速，精度为0.05%；采用示波器查看电流波形，计算谐波，谐波精度为0.1%。

电机低速（3 000 r/min）、中速（5 000 r/min、9 000 r/min）、高速（12 000 r/min、15 000 r/min）条件下采集每个工况下对应的电机输入功率P_in、输出功率P_out，从而计算得到各工况下的电机总损耗。

4.1 电机铁损分离试验

定子铁损需要从电机总损耗中分离得到，分离通过下式完成^[15]：

（7）P_Fe=∑P-P_Cu-P_w-P_∆-P_FeR

式中，∑p为电机总损耗；p_Cu为铜损；p_fw为机械损耗；P_∆为杂散损耗；P_FeR为转子铁损。

铜损计算公式为：

（8）P_Cu=m∑I_i²R

式中，m为相数；I_i为电机相电流的第i次谐波含量有效值；R为绕组相电阻。

机械损耗包括风摩损耗P_mpn和轴承摩擦损耗P_mpb，与转速n_N有关，计算公式为^[15]：

（9）P_fw=P_mpn+P_mpb=K_mG_pn_N+2D₂

n N 3

式中，K_m为轴承摩擦损耗系数；G_p为电机转子所受重力；D₂为紧圈外径；L为电枢铁心长。

杂散损耗计算公式为^[16]：

（10）P_∆=P_out×(0.025-0.005×lgP_in)

转子铁损无法直接计算，因此，应用有限元软件对转子铁损进行仿真计算。表7所示为电机低速、中速和高速条件下的电机各种损耗的实测和计算结果。

4.2 铁损数据分析

将表6中用铁心模型计算出的定子铁损与表7中通过测试从总损耗中分离计算的定子铁损进行对比，如图15所示。从图15中可以看出，铁损实测值随着电机转速的升高而增大。正弦等效法基于磁密的变化计算定子铁损，仅考虑了磁密波峰到波谷的变化的峰值，没有考虑波形波动时的细微变化对铁损的影响，也没有考虑到谐波对铁损的影响，计算得到的定子铁损与试验值差别较大；有限元法在计算定子铁损过程中依赖不同频率的铁损曲线计算铁损，忽略了旋转磁化造成的谐波对铁损计算的影响；谐波分析法考虑了旋转磁场和谐波在铁心中的影响，将磁密分解为径向和切向分量，并利用傅里叶分析得出了各点的谐波分量，定子铁损的计算结果比较接近试验值。

利用正弦等效法、谐波分析法、有限元法计算出电机在低速、中速、高速条件下的效率，并与实测效率进行对比，如表8所示。由表8可以看出，谐波分析法计算出的电机效率与实测值较为接近，电机在各工况下都可高效运行。

5 结束语

收起

本文对一台电动汽车用高速永磁同步电机进行了局部优化。应用田口法优化电机的部分参数，优化后的电机定子铁损降低了8%、转矩脉动系数减小了13.3%，电机运行更加平稳。同时用3种方法计算出电机定子铁损，其中谐波分析法考虑了磁密和谐波，计算的定子铁损更加精确，电机在各工况下均可高效运行，证明了定子铁损计算模型的有效性。

电机运行时，各种损耗将引起电机的发热。因此需要设计相应的散热通道，这一部分内容将在后续的研究中开展。

基金

收起

* 陕西省自然科学基础研究计划项目(2019JM-443)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

韩爱国, 苑昭, 郑青星, 等. 田口法在永磁同步电机优化设计中的应用[J]. 汽车实用技术, 2019(6): 40-44.

HAN

A G

, YUAN

, ZHENG

Q X

, et al. Application of Taguchi Method in Optimum Design of Permanent Magnet Synchronous Motors[J]. Automobile Applied Technology, 2019(6): 40-44.

[2]

SONG

J C

, DONG

, ZHAO

J W

, et al. Optimal Design of Permanent Magnet Liner Synchronous Motors Based on Taguchi Method[J]. IET Electric Power Applications, 2017, 11(1): 41-48.

[3]

GOPE

, GOEL

S K

. Design Optimization of permanent Magnet Synchronous Motor Using Taguchi Method and Experimental Validation[J]. International Journal of Emerging Electric Power Systems, 2021, 22(1): 9-20.

[4]

YAMAZAKI

, KUMAGAI

, IKEMI

, et al. A Novel Rotor Design of Interior Permanent-Magnet Synchronous Motors to Cope with Both Maximum Torque and Iron-Loss Reduction[J]. IEEE Transactions Industry Applications, 2013, 49(6): 2478-2486.

[5]

董剑宁, 黄允凯, 金龙, 等. 高速永磁电机设计与分析技术综述[J]. 中国电机工程学报, 2014, 34(27): 4640-4653.

DONG

J N

, HUANG

Y K

, JIN

, et al. Review on High Speed Permanent Magnet Machines Including Design and Analysis Technologies[J]. Proceedings of the CSEE, 2014, 34(27): 4640-4653.

[6]

沈建新, 李鹏, 郝鹤, 等. 高速永磁无刷电机电磁损耗的研究概况[J]. 中国电机工程学报, 2013, 33(3): 62-74.

SHEN

J X

, LI

, HAO

, et al. Study on Electromagnetic Losses in High-Speed Permanent Brushless Motors the State of the Art[J]. Proceedings of the CSEE, 2013, 33(3): 62-74.

[7]

赵海森, 罗应立, 刘晓芳, 等. 异步电机空载铁耗分布的时步有限元分析[J]. 中国电机工程学报, 2010, 30(30): 99-106.

ZHAO

H S

, LUO

Y L

, LIU

X F

, et al. Analysis on No-Load Iron Losses Distribution of a Synchronous Motor with Time-Stepping Finite Element Method[J]. Proceedings of the CSEE, 2010, 30(30): 99-106.

[8]

YAMAZAKI

, ISHIGAMI

. Rotor-Shape Optimization of Interior-Permanent-Magnet Motors to Reduce Harmonic Iron Losses[J]. IEEE Transactions on Industry Electronics, 2010, 57(1): 61-69.

[9]

王凤翔. 高速电机的设计特点及相关技术研究[J]. 沈阳工业大学学报, 2006, 28(3): 258-264.

WANG

F X

. Study on Design Feature and Related Technology of High-Speed Electric Machines[J]. Journal of Shenyang University of Technology, 2006, 28(3): 258-264.

[10]

王艾萌, 温云. 田口法在内置式永磁同步电机优化设计中的应用[J]. 华北电力大学学报(自然科学版), 2016, 43(3): 39-44.

WANG

A M

, WEN

. Application for Optimal Designing of Interior Permanent Synchronous Motor by using Taguchi Method[J]. Journal of North China Electric Power University (Natural Science Edition), 2016, 43(3): 39-44.

[11]

江善林. 高速永磁同步电机的损耗分析与温度场计算[D]. 哈尔滨: 哈尔滨工业大学, 2010.

JIANG

S L

. High-Speed Permanent Synchronous Motor Loss Analysis and Temperature Field Calculation[D]. Harbin: Harbin Institute of Technology, 2010.

[12]

BERTOTTI

. General Properties of Power Losses in Soft Ferromagnetic Material[J]. IEEE Transactions on Magnetics, 1988, 24(1): 621-630.

[13]

于明湖, 张玉秋, 张波. 永磁同步电动机铁耗计算方法分析[J]. 微特电机, 2015, 43(7): 21-25.

M H

, ZHANG

Y Q

. Analysis on the Iron Loss Calculation Method of Permanent Magnet Synchronous Motor[J]. Small and Special Electrical Machines, 2015, 43(7): 21-25.

[14]

戈宝军, 罗前通, 王立坤, 等. 高速永磁同步电动机铁耗分析[J]. 电机与控制学报, 2020, 24(4): 32-39.

B J

, LUO

Q T

, WANG

L K

, et al. Analysis of Iron Losses of High-Speed Permanent Synchronous Motor[J]. Electric Motors and Control, 2020, 24(4): 32-39.

[15]

李钟明, 刘卫国, 刘景林, 等. 稀土永磁电机[M]. 北京: 国防工业出版社, 1997.

Z M

, LIU

W G

, LIU

J L

, et al. Rare Earth Permanent Magnet Motor[M]. Beijing: National Defense Industry Press, 1997.

[16]

宁榕, 李巧莲, 黄坚. 三相异步电动机的杂散损耗研究[J]. 电机与控制应用, 2014, 41(5): 14-18.

NING

, LI

Q L

, HUANG

. Research on Stray Loss of Three-Phase Asynchronous Motor[J]. 2014, 41(5): 14-18.

2023年第卷第1期

PDF下载

391

175

引用本文

BibTeX

文章信息

doi: 10.19620/j.cnki.1000-3703.20220019

首发时间：2025-12-07
出版时间：2023-01-24

补充材料

相关文章

文章信息

作者

出版历史

修回日期：2022-02-04

基金

* 陕西省自然科学基础研究计划项目(2019JM-443)

作者信息

¹ 西安建筑科技大学，西安 710055

² 三门峡速达交通节能科技股份有限公司，三门峡 472000

通讯作者:

许檬（1997—），硕士研究生，研究方向为新型电力传动系统、新能源汽车驱动系统与控制， xumeng@xauat.edu.cn。

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/qcjs/CN/10.19620/j.cnki.1000-3703.20220019

分享至

全文二维码

扫描看全文

引用本文

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

参数	参数值	参数	参数值
定子槽数	48	额定功率/kW	60
极数	8	额定转速/r·min^-1	5 000
定子外径/mm	220	额定转矩/N·m	112
定子内径/mm	154	线径/mm	0.92
电机轴长/mm	140	铜径/mm	0.83
气隙/mm	0.8	并联支路数	8

参数

参数值

参数

参数值

定子槽数

额定功率/kW

极数

额定转速/r·min^-1

5 000

定子外径/mm

220

额定转矩/N·m

112

定子内径/mm

154

线径/mm

0.92

电机轴长/mm

140

铜径/mm

0.83

气隙/mm

0.8

并联支路数

因子水平	优化参数
水平1	8.4	0.65	2.5	23.0	2.5	13.0
水平2	8.7	0.70	3.0	23.5	3.0	13.5
水平3	9.0	0.75	3.5	24.0	3.5	14.0
水平4	9.3	0.80	4.0	24.5	4.0	14.5
水平5	9.6	0.85	4.5	25.0	4.5	15.0

因子水平

优化参数

h_m1

ω_m1

h_m2

ω_m2

水平1

8.4

0.65

2.5

23.0

2.5

13.0

水平2

8.7

0.70

3.0

23.5

3.0

13.5

水平3

9.0

0.75

3.5

24.0

3.5

14.0

水平4

9.3

0.80

4.0

24.5

4.0

14.5

水平5

9.6

0.85

4.5

25.0

4.5

15.0

试验次序	因子水平	P_Fe/W	T_avg/N·m	K_mb/%
1	1	1	1	1	1	1	840.5	294.9	1.88
2	1	2	2	2	2	2	847.1	277.2	1.15
3	1	3	3	3	3	3	815.6	284.8	1.94
4	1	4	4	4	4	4	856.7	290.2	2.03
5	1	5	5	5	5	5	826.6	296.9	2.14
6	2	1	2	3	4	5	929.2	290.4	1.71
7	2	2	3	4	5	1	886.2	277.3	3.55
8	2	3	4	5	1	2	821.1	286.5	2.03
9	2	4	5	1	2	3	835.5	281.9	2.09
10	2	5	1	2	3	4	854.5	270.7	1.91
11	3	1	3	5	2	4	849.2	292.6	1.93
12	3	2	4	1	3	5	832.4	291.2	1.88
13	3	3	5	2	4	1	975.3	281.1	2.24
14	3	4	1	3	5	2	863.1	272.1	1.87
15	3	5	2	4	1	3	785.8	271.6	1.93
16	4	1	4	2	5	3	943.2	290.9	2.03
17	4	2	5	3	1	4	903.6	290.9	2.03
18	4	3	1	4	2	5	836.2	277.9	1.64
19	4	4	2	5	3	1	863.2	272.6	1.95
20	4	5	3	1	4	2	854.6	273.9	2.03
21	5	1	5	4	3	2	842.4	294.7	2.01
22	5	2	1	5	4	3	901.4	279.6	1.69
23	5	3	2	1	5	4	901.9	281.7	1.86
24	5	4	3	2	1	5	820.2	278.9	1.82
25	5	5	4	3	2	1	835.5	272.9	2.18

试验次序

因子水平

P_Fe/W

T_avg/N·m

K_mb/%

h_m1

ω_m1

h_m2

ω_m2

840.5

294.9

1.88

847.1

277.2

1.15

815.6

284.8

1.94

856.7

290.2

2.03

826.6

296.9

2.14

929.2

290.4

1.71

886.2

277.3

3.55

821.1

286.5

2.03

835.5

281.9

2.09

854.5

270.7

1.91

849.2

292.6

1.93

832.4

291.2

1.88

975.3

281.1

2.24

863.1

272.1

1.87

785.8

271.6

1.93

943.2

290.9

2.03

903.6

290.9

2.03

836.2

277.9

1.64

863.2

272.6

1.95

854.6

273.9

2.03

842.4

294.7

2.01

901.4

279.6

1.69

901.9

281.7

1.86

820.2

278.9

1.82

835.5

272.9

2.18

参数因子	平均转矩方差	转矩脉动系数方差	铁损方差
r	240.80	80.27	0.021 5	15.6	189.46	10.2
δ	28.66	9.55	0.003 5	2.5	340.92	18.4
h_m1	16.95	5.65	0.022 4	16.3	106.27	5.8
ω_m1	4.35	1.45	0.013 5	9.8	265.05	14.3
h_m2	1.89	0.63	0.038 4	27.9	767.58	41.4
ω_m2	7.35	2.45	0.038 4	27.9	183.48	9.9
总	300	100	0.137 7	100	1 852.76	100

参数因子

平均转矩方差

转矩脉动系数方差

铁损方差

方差值
/(N·m)²

比重
/%

方差值
/%²

比重
/%

方差值
/W²

比重
/%

240.80

80.27

0.021 5

15.6

189.46

10.2

28.66

9.55

0.003 5

2.5

340.92

18.4

h_m1

16.95

5.65

0.022 4

16.3

106.27

5.8

ω_m1

4.35

1.45

0.013 5

9.8

265.05

14.3

h_m2

1.89

0.63

0.038 4

27.9

767.58

41.4

ω_m2

7.35

2.45

0.038 4

27.9

183.48

9.9

总

300

100

0.137 7

100

1 852.76

100

转速 /r·min^-1	P_Fe/W	T_avg/N·m	K_mb/%
2 000	811.7	777.3	284.2	280.0	2.04	1.80
15 000	5 151.2	4 647.8	290.3	276.6	2.10	1.30

转速
/r·min^-1

P_Fe/W

T_avg/N·m

K_mb/%

优化前

优化后

优化前

优化后

优化前

优化后

2 000

811.7

777.3

284.2

280.0

2.04

1.80

15 000

5 151.2

4 647.8

290.3

276.6

2.10

1.30

转速 /r·min^-1	计算方法	磁滞损耗 /W	涡流损耗 /W	附加损耗 /W	总铁损 /W
3 000	方法1	321.9	716.1	103.2	1 141.2
方法2	357.7	795.8	114.0	1 267.5
方法3	335.8	747.1	107.1	1 189.6
5 000	方法1	519.2	662.8	206.6	1 388.4
方法2	642.6	834.4	283.2	1 661.2
方法3	610.8	809.2	290.4	1 610.4
9 000	方法1	669.6	874.9	240.2	1 784.7
方法2	846.1	1 107.5	342.9	2 296.5
方法3	769.5	1 056.1	280.0	2 105.6
12 000	方法1	1 010.5	1 333.3	390.4	2 734.2
方法2	1 339.5	1 767.6	475.1	3 582.2
方法3	1 244.6	1 584.7	487.1	3 316.4
15 000	方法1	1 558.8	2 022.4	575.6	4 156.8
方法2	1 996.2	2 482.8	570.8	5 049.7
方法3	1 870.2	2 403.1	576.2	4 849.5

转速
/r·min^-1

计算方法

磁滞损耗
/W

涡流损耗
/W

附加损耗
/W

总铁损
/W

3 000

方法1

321.9

716.1

103.2

1 141.2

方法2

357.7

795.8

114.0

1 267.5

方法3

335.8

747.1

107.1

1 189.6

5 000

方法1

519.2

662.8

206.6

1 388.4

方法2

642.6

834.4

283.2

1 661.2

方法3

610.8

809.2

290.4

1 610.4

9 000

方法1

669.6

874.9

240.2

1 784.7

方法2

846.1

1 107.5

342.9

2 296.5

方法3

769.5

1 056.1

280.0

2 105.6

12 000

方法1

1 010.5

1 333.3

390.4

2 734.2

方法2

1 339.5

1 767.6

475.1

3 582.2

方法3

1 244.6

1 584.7

487.1

3 316.4

15 000

方法1

1 558.8

2 022.4

575.6

4 156.8

方法2

1 996.2

2 482.8

570.8

5 049.7

方法3

1 870.2

2 403.1

576.2

4 849.5

转速/r·min^-1	∑p/W	P_Cu/W	p_FeR/W	P_fw/W	P_∆/W	p_Fe/W
3 003	3 667.2	1 650.2	48.3	660.1	66.3	1 242.3
5 000.2	4 538.7	1 852.4	53.9	824.3	121.4	1 686.7
9 000.5	5 940.6	2 299.7	47.2	892.9	224.6	2 476.2
12 030	7 995.6	2 864.1	39.8	1 186.2	305.3	3 600.2
14 998	11 803.7	4 395.6	33.4	1 818.7	389.4	5 166.6

转速/r·min^-1

∑p/W

P_Cu/W

p_FeR/W

P_fw/W

P_∆/W

p_Fe/W

3 003

3 667.2

1 650.2

48.3

660.1

66.3

1 242.3

5 000.2

4 538.7

1 852.4

53.9

824.3

121.4

1 686.7

9 000.5

5 940.6

2 299.7

47.2

892.9

224.6

2 476.2

12 030

7 995.6

2 864.1

39.8

1 186.2

305.3

3 600.2

14 998

11 803.7

4 395.6

33.4

1 818.7

389.4

5 166.6

转速 /r·min^-1	电机效率/%
3 000	95.52	95.32	94.37	94.82
5 000	96.78	96.26	95.67	95.84
9 000	96.85	96.54	95.86	96.12
12 000	95.68	95.62	93.94	95.43
15 000	95.24	95.12	93.12	94.60

转速
/r·min^-1

电机效率/%

实测值

谐波分析法

正弦等效法

有限元法

3 000

95.52

95.32

94.37

94.82

5 000

96.78

96.26

95.67

95.84

9 000

96.85

96.54

95.86

96.12

12 000

95.68

95.62

93.94

95.43

15 000

95.24

95.12

93.12

94.60