Comparative Study of Multiphase IPMSMs For EV Applications
Comparative Study of Multiphase IPMSMs For EV Applications
EV Applications
Furkan Tokgoz1, Z.Q. Zhu1*, Xiao Chen1, Chaohui Liu2, Lianghui Yang3, Yiqing Yuan3, Han Yang3, and Yiming Huang3
1
Department of Electronic and Electrical Engineering, University of Sheffield, United Kingdom
2
National New Energy Vehicle Technology Innovation Center (NEVC), China
2023 26th International Conference on Electrical Machines and Systems (ICEMS) | 979-8-3503-1758-9/23/$31.00 ©2023 IEEE | DOI: 10.1109/ICEMS59686.2023.10345120
3
Xiamen NEVC Advanced Electric Powertrain Technology Innovation Centre, China
Email: [email protected]
1
Authorized licensed use limited to: University of Ulsan. Downloaded on December 18,2023 at 03:34:31 UTC from IEEE Xplore. Restrictions apply.
ripple and PM volume. The Pearson correlation coefficients
of the design variables are calculated by
𝜌𝑥𝑖,𝑦𝑖
𝑁 ∑ 𝑥𝑖𝑦𝑖−∑ 𝑥𝑖 ∑ 𝑦𝑖
= �𝑁 ∑ 𝑥2−(∑ 𝑥 )2�𝑁 ∑ 𝑦2−(∑ 𝑦 )2
wr Lc 𝑖 𝑖 𝑖 𝑖
where 𝛾𝛾 is the split ratio, 𝑟𝑟 and 𝑟𝑠 are the rotor and stator
outer radii, respectively. The maximum value of the back-
iron width is limited by the other design parameters such as
split ratio. As the split ratio becomes lower, the maximum
value of the back-iron width also decreases to eliminate
infeasible solutions from the optimization search space. The
rotor structure is defined by the split ratio, 𝛾𝛾, the width of
the ribs, 𝑤𝑟 , the centre post width, 𝑤𝑐 , the angle between
two magnets, 𝛼𝑚, the magnet thickness, 𝑡𝑚 , the magnet
Fig. 2. Pearson correlation coefficients of design variables.
width,
𝑤𝑚, and the distance of the magnets from the centre post, Another important part of the optimization is the design
𝑑𝑚. The maximum width of the magnet is limited by the constants which are shown in Table III. The outer diameter
other of the stator, air-gap length, slot opening width, and inner
parameters as the width of the duct changes with these diameter of the rotor are kept constant since they are limited
parameters. The iron bridge width is kept constant at 1.3 by mechanical constraints. The number of stator slots is kept
mm since it is limited by the mechanical stress created by constant for the three-phase IPMSM while it is changed for
the centrifugal force. In addition to these parameters, phase multiphase machines. The current density in the coils and
advance angle of the stator current, 𝛼𝑎𝑑𝑣 , is also added to the fill factor in the slot are constants to have a fair
the optimization algorithm since the maximum torque angle comparison between designs. Phase current is limited by the
will change as the rotor structure is changed. The upper and following equation:
lower limits of the design variables are presented in Table II.
Lower Upper
𝐼𝐿𝑚 3𝐼𝐿
Variables Limit Limit = 𝑚
Tooth width, tm, (mm) 3 4.75 maximizing the average torque and minimizing the torque
Back-iron width, 𝐿𝑐, (mm) 10 35
Magnet angle, 𝛼𝑚, (deg) 130 170
Split ratio, 𝛾𝛾 0.55 0.75
Rib width, 𝑤𝑟, (mm) 1 10
Magnet distance, 𝑑𝑚, (mm) 0 18
Magnet width, 𝑤𝑚, (mm) 10 40
Magnet thickness, 𝑡𝑚, (mm) 2 8
Centre post width, 𝑤𝑐, (mm) 3 10
Phase advance, 𝛼𝑎𝑑𝑣, (deg) 40 60
3
Authorized licensed use limited to: University of Ulsan. Downloaded on December 18,2023 at 03:34:31 UTC from IEEE Xplore. Restrictions apply.
𝐿𝑒𝑛𝑑 = 𝜋 𝑝𝑠𝑙𝑜𝑡 Table IV which shows that the 60-slot/6-pole machine has
the highest winding factor compared to other combinations.
where 𝑀𝐿𝑇 is the mean length per turn, 𝐿𝑚𝑜𝑡𝑜𝑟 and 𝐿𝑒𝑛𝑑 are It also has a high winding factor for the third harmonic that
the machine axial length and the end winding length. 𝑝𝑠𝑙𝑜𝑡 is can be used to increase the average output torque of five-
the radius of the slot pitch. phase machines with the third harmonic injection method.
TABLE III. DESIGN CONSTANTS USED IN OPTIMIZATION
However, the effectiveness of the third harmonic injection is
determined by the winding configuration and the air-gap flux
Parameters Value density distribution. Although the winding used in the five-
Stator outer diameter (mm) 225 phase machine is suitable for having a high third harmonic
Rotor inner diameter (mm) 70 injection thanks to the high winding factor, the third
Axial length (mm) 134 harmonic at the air-gap flux density cannot be increased in
Number of stator slots 54
this IPMSM, since the pole arc is limited. The limited pole
arc results in cancelling out the third harmonic at the air gap
Slot opening width (mm) 2.5
flux density. Thus, the third harmonic injection is not used in
Airgap (mm) 0.7
the optimization of the five-phase IPMSM.
2
Current density (A/mm ) 37 The feasible slot/pole number combinations for the six-
Fill factor 0.4 phase machine are presented in Table IV. Although the 36-
Phase current (A) <1100 slot/6-pole machine has a higher winding factor it is not
Copper loss (kW) <20 chosen in this paper as it has a unity winding factor for all
harmonics and cogging torque frequency of the 36-slot/6-
The end winding shape is assumed to be semi-circular. pole machine is 12. This becomes an issue with a six-phase
𝑝𝑠𝑙𝑜𝑡 is calculated by machine since the torque ripple caused by the
electromagnetic torque is at the same frequency which results
in significantly high torque ripples in the output torque.
𝑝𝑠𝑙𝑜𝑡 = 2𝜋 𝑟𝑠𝑙𝑜𝑡
𝑡𝑐𝑜𝑖𝑙
𝑄𝑄𝑠 TABLE IV. WINDING FACTORS OF DIFFERENT WINDING CONFIGURATIONS
where 𝑟𝑠𝑙𝑜𝑡 is the radius of the middle of the slot, 𝑡𝑐𝑜𝑖𝑙 is the
Number of Stator
throw of the coil, and 𝑄𝑄𝑠 is the number of stator slots. The Slots
Fundamental Third Fifth Seventh
resistivity of the copper is dependent upon the temperature of Three-phase –
0.960 0.667 0.218 0.177
the coils with (7). 54-slot/6-pole
Five-phase –
0.973 0.772 0.462 0.176
40-slot/6-pole
𝜌𝑐𝑢,𝑓𝑖𝑛𝑎𝑙 = 𝜌𝑐𝑢(1 + 0.0039(𝑇𝑚𝑜𝑡𝑜𝑟 − 20))
Five-phase –
0.980 0.830 0.577 0.300
45-slot/6-pole
where 𝜌𝑐𝑢,𝑓𝑖𝑛𝑎𝑙 is the resistivity of the copper at the selected Five-phase –
0.982 0.848 0.616 0.344
50-slot/6-pole
temperature 𝜌𝑐𝑢 is the resistivity of the copper at 20 °C, and Five-phase –
0.988 0.891 0.707 0.454
𝑇𝑚𝑜𝑡𝑜𝑟 is the temperature of the windings. The operating 60-slot/6-pole
temperature of the machine is chosen as 80 °C. The Five-phase –
0.980 0.830 0.578 0.303
remanence flux density and magnetic coercivity of the 70-slot/6-pole
Six-phase –
NdFeB magnets are changed accordingly. Finally, phase 36-slot/6-pole
1 1 1 1
resistance, 𝑅𝑝ℎ𝑎𝑠𝑒, of the machine can be calculated by Six-phase –
0.991 0.924 0.793 0.609
72-slot/6-pole
Nine-phase –
𝑅𝑝ℎ𝑎𝑠𝑒 𝜌𝑐𝑢,𝑤𝑤𝑖𝑛𝑑𝑖𝑛𝑔 𝑀𝐿𝑇 𝑁𝑡𝑢𝑟𝑟𝑛 54-slot/6-pole 1 1 1 1
= 𝐴𝑐𝑜𝑖𝑙
4
Authorized licensed use limited to: University of Ulsan. Downloaded on December 18,2023 at 03:34:31 UTC from IEEE Xplore. Restrictions apply.
(a) (b)
(c) (d)
Fig. 3. Optimization results of machines. (a) Three-phase machine. (b) Five-phase machine. (c) Six-phase machine. (d) Nine-phase machine.
Fig. 4. Designed multiphase machines. (a) Tesla Model 3 IPMSM. (b) Three-phase. (c) Five-phase. (d) Six-phase. (e) Nine-phase.
5
Authorized licensed use limited to: University of Ulsan. Downloaded on December 18,2023 at 03:34:31 UTC from IEEE Xplore. Restrictions apply.
TABLE VI. COMPARISON OF MULTIPHASE MACHINES AT 5000
R/MIN AND MAXIMUM PHASE CURRENTS IN MAXWELL torque ripple when compared to other multiphase
machines. The reason for this discrepancy can be
explained by the slot/pole number combinations. The
cogging torque frequency of a PM machine can be
Tesla Three- Five- Six- Nine-
Parameter
Model 3 phase phase phase phase
calculated by
Output power (kW) 230 248 255 261 254
Output torque (Nm) 439 474 488 499 485 𝑓𝑐𝑜𝑐𝑐𝑐𝑐𝑖𝑛𝑐𝑐
Torque ripple (Nm/%) 16.2/3.7 55/11.6 18/3.7 24.4/4.9 43/8.9 2𝑁𝑠𝑓
= 𝐺𝐶𝐷(𝑁𝑠,2𝑝)
PM volume (cm3) 18.8 16.8 16.8 16.73 16.83
Torque per PM volume
23.4 28.2 29.1 29.8 28.8 where 𝑓𝑐𝑜𝑐𝑐𝑐𝑐𝑖𝑛𝑐𝑐 is the cogging torque frequency, 𝑁𝑠 is the
(Nm/cm3)
RMS phase current (A) 1100 1100 635 547 355 number of stator slots, 𝑓 is the electrical frequency,
Iron losses (kW) 0.87 1.04 1.11 1.05 1.11 𝐺𝐶𝐷(𝑁𝑠, 2𝑝) is the greatest common divisor of the
Efficiency (%) 91.7 92.2 92.37 92.56 92.33
number of stator slots, 𝑁𝑠 , and the number of poles, 2𝑝 .
TABLE VII. CHANGE IN IMPORTANT MACHINE PARAMETERS AFTER
The cogging frequency of a 54-slot/6-pole machine is the
OPTIMIZATION 18th harmonic of the electrical frequency. Similarly, the
torque ripple created by the electromagnetic torque of a
Tesla Model Optimum nine- phase machine is at the same frequency which
Variables 3 Design
creates an even higher torque ripple when combined with
Electrical loading (A/mm) 166 187
Magnetic loading (T) 0.8 0.96 the cogging torque.
300 300
torque improvement under third harmonic current injection
200 200
condition," IEEE Transactions on Magnetics, vol. 58, no. 8, pp. 1-
100 100
0
6, 2022.
0
0 3000 6000 900
Speed (r/ min)
12000 15000 18000 [7] S. Sadeghi, L. Guo, H. A. Toliyat, and L. Parsa, "Wide operational
Speed (r/min)
(d) speed range of five phase permanent magnet machines by using
(c)
500 Efficiency different stator winding configurations," IEEE Transactions on
400 0.5 0.75 0.98
Industrial Electronics, vol. 59, no. 6, pp. 2621-2631, 2012.
Torque (Nm)
7
Authorized licensed use limited to: University of Ulsan. Downloaded on December 18,2023 at 03:34:31 UTC from IEEE Xplore. Restrictions apply.