Assessment of the Energy Consumption and Drivability Performance of an IPMSM-Driven Electric Vehicle Using Different Buried Magnet Arrangements

: This study investigates the inﬂuence of the buried magnet arrangement on the efﬁciency and drivability performance provided by an on-board interior permanent magnet synchronous machine for a four-wheel-drive electric car with two single-speed on-board powertrains. The relevant motor characteristics, including ﬂux-linkage, inductance, electromagnetic torque, iron loss, total loss, and efﬁciency, are analyzed for a set of six permanent magnet conﬁgurations suitable for the speciﬁc machine, which is controlled through maximum-torque-per-ampere and maximum-torque-per-voltage strategies. Moreover, the impact of each magnet arrangement is analyzed in connection with the energy consumption along four driving cycles, as well as the longitudinal acceleration and gradeability performance of the considered vehicle. The simulation results identify the most promising rotor solutions, and show that: (i) the appropriate selection of the rotor conﬁguration is especially important for the driving cycles with substantial high-speed sections; (ii) the magnet arrangement has a major impact on the maximum motor torque below the base speed, and thus on the longitudinal acceleration and gradeability performance; and (iii) the conﬁgurations that excel in energy efﬁciency are among the worst in terms of drivability, and vice versa, i.e., at the vehicle level, the rotor arrangement selection is a trade-off between energy efﬁciency and longitudinal vehicle dynamics.


Introduction
In recent years, interior permanent magnet synchronous machines (IPMSM) have been widely used for electric vehicle (EV) powertrains [1,2], due to their superior electromagnetic performance, e.g., in terms of high airgap-flux and power densities, low rotor losses, and high efficiency for a wide speed range. Because of the variety of operating conditions of electric powertrains, the electromagnetic maps of the motors are essential tools for designing and evaluating the drives, and have been widely used in the literature [3][4][5].
In IPMSMs, the rotor is the most complex part, because of the permanent magnet (PM) arrangement and its effect on the overall electromagnetic performance. The buried design of the magnets prevents their separation from the rotor, despite the significant centrifugal force at high speed, and produces a hybrid torque, combination of magnet and reluctance contributions, which results in high efficiency [6]. Many complex magnet arrangements have been proposed for achieving specific characteristics [7][8][9]. References [7][8][9][10][11][12][13] discuss the performance improvement of IPMSMs, associated with various magnet arrangements and shaping techniques. For example, in [10], Kim et al. investigate the demagnetization performance of permanent magnets within three kinds of rotor types for IPMSMs, as well as the resulting motor performance in terms of rated torque, output current, and efficiency. In [11], Cirani et al. consider a novel rotor design for IPMSMs. The advantage of their solution is a major reduction in magnetic flux leakage in the rotor iron bridges, which represents the main limitation of this machine topology. Kano [12] studies a novel flux barrier design for torque ripple reduction in saliency-based sensorless drive concentrated-winding IPMSMs. Zhu et al. [13] adopt magnet shaping techniques to improve the torque capability in brushless AC (alternating current) machines, by optimizing the third order harmonic in an arrangement with inverse cosine shape airgap and sinusoidal shape magnet, which results in 11% average torque increase.
The literature extensively analyzes the influence of the electric motor topology on the resulting efficiency maps [14]. Several methodologies have been proposed for generating efficiency maps, with the purpose of predicting the system performance through vehicle simulations, and select the drive/s accordingly [3,4,14]. For example, in [15], the motor efficiency map is generated from a time-domain two-dimensional (2-D) finite element analysis (FEA) model in ANSYS Maxwell, which calculates the dand q-axis currents, and is coupled with a multi-objective optimization. During high-speed EV operation, IPMSMs are typically subject to flux weakening; in this condition, the motor iron loss, more precisely the eddy-current loss, tends to significantly increase, thus reducing efficiency. In [16], Mohammadi and Lowther consider the 2010 Toyota Prius IPMSM and a PMassisted synchronous reluctance machine, and use nonlinear motor control formulations, i.e., maximum-torque-per-ampere (MTPA), flux weakening, and maximum-torque-pervoltage (MTPV), in the computation of the efficiency maps, while accounting for both saturation and cross-coupling effects. MTPA and MTPV algorithms [17][18][19][20][21][22][23][24] are widely used to enhance the efficiency and control performance of IPMSMs, below and above the base speed, respectively. These strategies can be implemented through offline-generated maps, reducing the online computational burden but increasing the memory requirements, or through online optimization. The conventional optimality criteria for MTPA and MTPV do not consider magnetic core saturation, and therefore result in deviations from the optimal trajectory, especially at higher load currents. Recent research [25] proposes mathematical formulations to account for magnetic core saturation in the MTPA and MTPV criteria, and obtains the corresponding optimal solutions.
Although the available literature discusses various PM arrangements for EV traction motors [14], to the best of the authors' knowledge there is a lack of comprehensive comparisons of the resulting motor efficiency characteristics, and evaluations of the implications at the EV level, e.g., in terms of energy consumption during driving cycles and drivability performance. For a case study IPMSM for a four-wheel-drive EV with two single-speed on-board powertrains, this research targets the identified gap with the following main novel contributions: • Simulation analysis, mapping, and comparison of the electromagnetic characteristics of the IPMSM with six of the most common buried magnet arrangements, for the same target magnet volume, demagnetization rating, rotor yoke volume per magnet volume, and stator configuration; • Comparison of the resulting EV energy consumption along a comprehensive set of driving cycles, as well as drivability performance, and identification of the most suitable rotor configurations.
The remainder of the manuscript is organized as follows: Section 2 presents the considered IPMSM rotor topologies and their parametric design process; Section 3 describes the adopted analytical models; Section 4 compares the motor maps; Section 5 analyzes the EV energy efficiency and drivability performance, and is followed by the conclusions.

Considered Rotor Configurations
The simulated three-phase IPMSM with inner rotor topology has 48 slots and 8 poles, with a total mass of 10.5 kg. Each phase has 110 conductors, with a net fill factor equal to 0.51. The magnet and core materials are NdFeB and M330-35A, respectively. Figure 1 presents the three-phase winding distribution of the stator, and the slot-related parameters (given in Figure 1b), while Table 1 reports the geometric stator specifications, which are the same for all considered arrangements. Figure 1c shows the rotor parameters. The rotor arrangements depend on the PM layouts, which affect the electromagnetic behavior of the machines. For a fair comparison, the total PM volume, V m , of the considered configurations is identical, and set to 24 cm 3 . All magnets are radially magnetized. Figure 2 illustrates the considered PM arrangements, commonly investigated by researchers, e.g., see [10,11,13,[26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41], in a cross-sectional view: (i) a rectangular magnet with a rectangular and open window (referred to as R-ROW in the remainder) [10,26,27], see Figure 2a; (ii) an arrangement similar to the previous one, but with a semicircular open window (R-SOW) [11], see Figure 2b; (iii) a rectangular magnet with an isosceles-trapezoidshaped open window (R-IOW) [28][29][30][31][32][33], see Figure 2c; (iv) two separate rectangular magnets, divided by a central leg, with a semicircular open window (RCL-SOW), similarly to [35][36][37], see Figure 2d. The material used in the central leg is the same as for the rotor yoke; (v) a rectangular magnet with a central leg, and an isosceles-trapezoid-shaped open window (RCL-IOW) [34,38,40], see analyzes the EV energy efficiency and drivability performance, and is followed by th conclusions.

Considered Rotor Configurations
The simulated three-phase IPMSM with inner rotor topology has 48 slots and 8 poles with a total mass of 10.5 kg. Each phase has 110 conductors, with a net fill factor equal t 0.51. The magnet and core materials are NdFeB and M330-35A, respectively. Figure presents the three-phase winding distribution of the stator, and the slot-related parameters (given in Figure 1b), while Table 1 reports the geometric stator specifications which are the same for all considered arrangements. Figure 1c shows the rotor parameters The rotor arrangements depend on the PM layouts, which affect the electromagneti behavior of the machines. For a fair comparison, the total PM volume, , of th considered configurations is identical, and set to 24 cm 3 . All magnets are radially magnetized. Figure 2 illustrates the considered PM arrangements, commonly investigated by researchers, e.g., see [10,11,13,[26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41], in a cross-sectional view: (i) a rectangular magne with a rectangular and open window (referred to as R-ROW in the remainder) [10,26,27] see Figure 2a; (ii) an arrangement similar to the previous one, but with a semicircular open window (R-SOW) [11], see Figure 2b; (iii) a rectangular magnet with an isosceles trapezoid-shaped open window (R-IOW) [28][29][30][31][32][33], see Figure 2c; (iv) two separat rectangular magnets, divided by a central leg, with a semicircular open window (RCL SOW), similarly to [35][36][37], see    The first three magnet arrangements, i.e., R-ROW, R-SOW, and R-IOW, use differ shapes of the open windows to reduce the leakage flux at both ends of the magnet. T second and third arrangements are more practical than the first one from the viewpo of the manufacturing process to fix the magnets. The central legs in RCL-SOW and RC IOW improve the ratio of reluctance torque to magnet torque and saliency, in addition reducing the cogging torque. Although they increase the leakage flux, the higher d-a inductance can improve the flux weakening capability. The last magnet arrangeme H-NOW, is designed to exploit the harmonic components of the magnetomotive force a permeance.  The first three magnet arrangements, i.e., R-ROW, R-SOW, and R-IOW, use different shapes of the open windows to reduce the leakage flux at both ends of the magnet. The second and third arrangements are more practical than the first one from the viewpoint of the manufacturing process to fix the magnets. The central legs in RCL-SOW and RCL-IOW improve the ratio of reluctance torque to magnet torque and saliency, in addition to reducing the cogging torque. Although they increase the leakage flux, the higher d-axis inductance can improve the flux weakening capability. The last magnet arrangement, H-NOW, is designed to exploit the harmonic components of the magnetomotive force and permeance.  Figure 3 presents the effect of the main geometric magnet parameters, i.e., magnet width and pole arc. The sensitivity analysis, using the design of experiment/regression (DOE/R) method [42], shows that the increase in the pole arc amplitude brings a reduction in the electromagnetic torque ( Figure 3a) and torque ripples level (Figure 3b). The output power increases with increasing magnet width ( Figure 3c); however, the magnet width also provokes an increase in the peak-to-peak cogging torque, see Figure 3d. The demagnetization effects are covered in Figure 3e,f. The increase in the magnet width results in higher magnetic flux density, and consequently higher torque and power. However, this increase can critically cause demagnetization risks, which vary among the rotors, see Figure 3f, i.e., the R-IOW, R-SOW, and R-ROW arrangements have the highest magnetic field capability. This advantage becomes vital for high-speed applications because the risk of demagnetization significantly increases during full-load and high-speed operation. Figure 3 presents the effect of the main geometric magnet parameters, i.e. width and pole arc. The sensitivity analysis, using the design of experiment/re (DOE/R) method [42], shows that the increase in the pole arc amplitude brings a r in the electromagnetic torque ( Figure 3a) and torque ripples level (Figure 3b). Th power increases with increasing magnet width ( Figure 3c); however, the magnet w provokes an increase in the peak-to-peak cogging torque, see Figure  demagnetization effects are covered in Figure 3e,f. The increase in the magnet widt in higher magnetic flux density, and consequently higher torque and power. How increase can critically cause demagnetization risks, which vary among the rotors, s 3f, i.e., the R-IOW, R-SOW, and R-ROW arrangements have the highest magn capability. This advantage becomes vital for high-speed applications because th demagnetization significantly increases during full-load and high-speed operation For fairness of comparison, the main magnet parameters for the co arrangements have been selected through sensitivity analysis, as a trade-off in electromagnetic performance, while accounting for saturation rating. The parameters are reported in Table 2, which includes ranges for the parameters t among the configurations. The axial length of the magnets is the same as the leng rotor for all cases. For fairness of comparison, the main magnet parameters for the considered arrangements have been selected through sensitivity analysis, as a trade-off in terms of electromagnetic performance, while accounting for saturation rating. The selected parameters are reported in Table 2, which includes ranges for the parameters that vary among the configurations. The axial length of the magnets is the same as the length of the rotor for all cases.

Finite Element Analysis Model
The objective of the simulation study is the electromagnetic evaluation of the magnet arrangements and their influence on dq-axis flux linkage, dq-axis synchronous inductance, back-electro-motive force (EMF), electromagnetic torque, as well as Joule and iron losses. The electromagnetic parameters are calculated with a 2-D FEA approach using the software package Flux. In the FEA simulations, constant settings are defined for: (i) the maximum root mean square (RMS) line current, provided by the inverter to the motor, set to 103 A, i.e., it is assumed that the operating limits of the powertrain are caused by the inverter; (ii) the maximum motor speed, i.e., ω max = 10,000 rpm; and (iii) the current coefficient, set to 0.8. Moreover, the speed range is discretized with steps of 10 rpm; the adopted airgap mesh density is 1.5; 100 operating points are considered per electrical torque period. Finally, steady-state thermal conditions are assumed. By supplying sinusoidal currents into the steady-state FEA model, the influence of the various magnet arrangements on the IPMSM performance is observed from the motor maps, which are generated over the whole motor speed (ω mot ) range, by applying the MTPA and MTPV strategies. The FEA simulation outputs several variables, the main ones being the magnetic flux density, B, and PM flux linkage, λ m , which-in the post-processing phase-are used to calculate the dand q-axis components of the inductance, L d and L q , as well as the different motor loss components, through analytical formulations.

Main Analytical Formulations
In the steady-state dq-framework rotating with electrical frequency ω e [22], the voltage equations of the IPMSM are: where v d , v q , i d , and i q are the dq-axis voltages and currents; and R s is the stator winding resistance. i d and i q are related to the stator current, i s , through the current angle β: The electromagnetic and reluctance torque, T e , is computed as: In the generation of the results, MTPA is adopted below the base speed, ω b , while MTPV is used above ω b . The MTPA algorithm is based on the solution of an optimization problem [22][23][24][25], which minimizes the continuous stator current, i s , that generates the reference torque, T re f , while meeting current and voltage constraints: where I max is the stator current limit, which depends on the capability of the inverter; V max is the stator voltage limit; V dc is the DC (direct current) voltage of the inverter; and v s is the magnitude of the stator voltage, given by: The solution of the MTPA optimization problem results in the following current angle, β MTPA , which was imposed in the generation of the motor maps through the FEA model: MTPV is an effective strategy to obtain maximum torque when operating the motor in the flux weakening region, i.e., above the base speed [22][23][24][25]. The MTPV optimization problem can be formulated as: which, without considering core saturation, leads to: Condition (8) was imposed in the map generation process for ω mot ≥ ω b . The total motor power loss, W mot,loss , is computed as the sum of the Joule loss W J , iron loss W f e , and mechanical loss W mech : where W J is: In (10), m is the number of phases; R ph is the phase resistance; and I ph is the RMS value of the phase current. W f e depends on the magnetic flux density, B. For each node of the iron loss mapping, B is output by a semi-numerical method based on the integration of the flux density in the airgap. B is computed as a function of the angular position θ, by considering all relevant parts of the machine, i.e., the stator foot teeth, stator teeth, and yokes. A mathematical transformation is applied to obtain the time derivative of B for each speed interval: W f e is given by: where K h , K c , K e , and K f are hysteresis, eddy-current, excessive, and stacking coefficients; α h , α c , and α e are fitting constants for the hysteresis, eddy-current, and excessive terms; V f e is the iron lamination volume; B max is the maximum magnetic flux density; and t e is the electrical period. W mech is defined as: where W 0 is the mechanical loss at a reference speed; ω mot,0 is the nominal rotor speed; and K is the speed exponent. The magnitude of the mechanical loss is estimated for all motor configurations, and is the lowest for the R-SOW arrangement. The motor efficiency, η mot , is computed as: where T mot is the mechanical torque at the motor shaft during steady-state operation.

Electromagnetic Performance Mapping of the IPMSMs
Based on the models in Section 3, the electromagnetic characteristics are generated as functions of flux linkage, across the whole speed-torque range. Table 3 reports the main electromagnetic properties for a current density of 6 A/mm 2 . The produced magnetic field strength H (mean value) due to the permeances (with relative permeability of 1.1) and the flux per pole have their highest values in R-SOW, R-IOW, and H-NOW. The PM flux λ m is particularly high for the R-IOW, R-SOW, and H-NOW arrangements, which are characterized by the lowest flux leakage. High values of the saliency ratio, ρ max , such as those of the R-IOW and R-SOW arrangements, imply higher power capability and lower demagnetization, i.e., better high-speed performance. L q has the highest values for RCL-SOW and RCL-IOW; moreover, the d-axis inductance L d of RCL-SOW and RCL-IOW is rather high because the magnet is cut into two parts, which ultimately results in lower ρ max . Figure 4 reports the flux linkage as a function of rotor position, with the highest magnitude shown by R-IOW and RCL-SOW, followed by R-ROW. Results of fast Fourier transform analysis show total harmonic distortions for the studied magnet arrangements ranging from 9.28 to 11.26%, with the lowest values belonging to R-SOW and R-IOW. Figure 5 includes the iron loss maps at 6000 rpm, as functions of the dq-axis currents. Each map is divided into four quadrants, defined by dashed lines in the plots. In particular, for all layouts, the maximum iron losses are in the quadrant with dq-axis currents in the ranges −100 A ≤ I d ≤ −50 A and 0 A ≤ I q ≤ 50 A. The variations of the shape of the iron loss iso-lines highlight that the current amplitude and its angle are highly dependent on the magnet arrangement, which affects the power-related computations in the entire operating range. In general, the lowest maximum iron losses occur for R-ROW, RCL-SOW, and RCL-IOW.    Figure 6 presents the total loss maps as functions of torque, while Table 4 reports th RMS value of the phase current , as well as the Joule loss , iron loss , mechanica loss , and total loss , , resulting from (9)-(13), at the maximum speed and torque of each configuration. For facilitating the visual analysis, the maps in Figure 6 sub divide the operating speed range of the motor into the following three regions, namely (i) the starting-torque low-speed (ST-LS) region, below the base speed of the machine; (ii the torque attenuation phase (AP), where the maximum torque is subject to a rather stee reduction with speed. This region covers the intermediate range of available speeds; and (iii) the highway driving high-speed region (HD-HS), covering the top range of achievabl speeds of the machine, where the maximum torque reduction with speed is more moder ate than in (ii). Depending on the vehicle type and application, the performance can b prioritized for different regions, e.g., the AP region can be the most significant one for racing car, while the ST-LS region can be the most important one for an urban EV.
For example, the map for the R-ROW configuration in Figure 6a shows a constan torque range capability of 135 Nm, with a maximum loss of 1084 W at the base speed. I the AP speed range, the torque is subject to a significant reduction, with a maximum los of 1807 W at 65 Nm and 6000 rpm. Under HD-HS operation, the total loss reaches 3750 W at 10,000 rpm, speed at which the maximum torque is ~30 Nm. In the R-SOW configura tion in Figure 6b, the total losses are higher than for the R-ROW layout across the entir speed range; however, the torque production has considerably increased to 158 Nm, wit a ~1300 W maximum loss at the base speed. During HD-HS operation, the torque capabil ity is slightly higher than for R-ROW, but this corresponds to an increase in total loss reaching a value in excess of 4000 W. Similar analyses can be made for all configurations Figure 7 reports the resulting power efficiency maps. All configurations have a regio with efficiency values in excess of 0.97 at medium-to-high torque, for a speed range acros the ST-LS and AP sectors. Such premium efficiency region covers a relatively small torqu range in R-ROW, RCL-SOW, and RCL-IOW, and a wider torque range in R-SOW, R-IOW   Table 4 reports the RMS value of the phase current I ph , as well as the Joule loss W J , iron loss W f e , mechanical loss W mech , and total loss W mot,loss , resulting from (9)-(13), at the maximum speed and torque of each configuration. For facilitating the visual analysis, the maps in Figure 6 subdivide the operating speed range of the motor into the following three regions, namely: (i) the starting-torque low-speed (ST-LS) region, below the base speed of the machine; (ii) the torque attenuation phase (AP), where the maximum torque is subject to a rather steep reduction with speed. This region covers the intermediate range of available speeds; and (iii) the highway driving high-speed region (HD-HS), covering the top range of achievable speeds of the machine, where the maximum torque reduction with speed is more moderate than in (ii). Depending on the vehicle type and application, the performance can be prioritized for different regions, e.g., the AP region can be the most significant one for a racing car, while the ST-LS region can be the most important one for an urban EV.
For example, the map for the R-ROW configuration in Figure 6a shows a constant torque range capability of 135 Nm, with a maximum loss of 1084 W at the base speed. In the AP speed range, the torque is subject to a significant reduction, with a maximum loss of 1807 W at 65 Nm and 6000 rpm. Under HD-HS operation, the total loss reaches 3750 W at 10,000 rpm, speed at which the maximum torque is~30 Nm. In the R-SOW configuration in Figure 6b, the total losses are higher than for the R-ROW layout across the entire speed range; however, the torque production has considerably increased to 158 Nm, with a~1300 W maximum loss at the base speed. During HD-HS operation, the torque capability is slightly higher than for R-ROW, but this corresponds to an increase in total loss, reaching a value in excess of 4000 W. Similar analyses can be made for all configurations. Figure 7 reports the resulting power efficiency maps. All configurations have a region with efficiency values in excess of 0.97 at medium-to-high torque, for a speed range across the ST-LS and AP sectors. Such premium efficiency region covers a relatively small torque range in R-ROW, RCL-SOW, and RCL-IOW, and a wider torque range in R-SOW, R-IOW, and H-NOW, while the speed range of the top efficiency region is rather similar for all maps. During high-speed operation, i.e., in the HD-HS range, the power efficiency of all configurations significantly decreases to values well below 0.90, and in some cases below 0.80. Table 5 summarizes the electromagnetic performance of the considered configurations, with the following highlights: (i) the peak phase voltage is the highest in the R-SOW configuration, followed by R-IOW and H-NOW; (ii) the energy production of the permanent magnets is the highest for R-IOW, followed by R-SOW and H-NOW; (iii) at ω b , the highest power is produced by R-SOW, R-IOW, and H-NOW. The maximum output power at ω b is subject to a major variation depending on the rotor configuration, and ranges from~42.3 tõ 51.7 kW, i.e., the EV with the considered set of motors will have very different longitudinal acceleration performance characteristics; (iv) at ω max , the top power capability is provided by R-SOW, followed by R-IOW; (v) the usable torque per magnet volume, i.e., T e /V m , is the highest for R-SOW at ω b and ω max ; and (vi) R-SOW has the highest maximum efficiency at ω b , while H-NOW has the highest maximum motor efficiency at ω max . In general, based on the table and Figures 6 and 7, the R-SOW and R-IOW configurations, followed by H-NOW, appear to have overall desirable characteristics in terms of torque density production and premium efficiency area, which would make them good candidates for EV implementation.
To complete the analysis, Figure 8 shows the steady-state temperature distributions in the rotors, based on the FEA model results at the peak torque and base speed, and includes iso-temperature lines. The maximum temperatures are rather uniform among the magnet arrangements, with approximately 4 • C variations among the configurations. The lowest maximum temperature, equal to 64.5 • C, occurs for the H-NOW arrangement. and H-NOW, while the speed range of the top efficiency region is rather similar for all maps. During high-speed operation, i.e., in the HD-HS range, the power efficiency of all configurations significantly decreases to values well below 0.90, and in some cases below 0.80. Table 5 summarizes the electromagnetic performance of the considered configurations, with the following highlights: (i) the peak phase voltage is the highest in the R-SOW configuration, followed by R-IOW and H-NOW; (ii) the energy production of the permanent magnets is the highest for R-IOW, followed by R-SOW and H-NOW; (iii) at , the highest power is produced by R-SOW, R-IOW, and H-NOW. The maximum output power at is subject to a major variation depending on the rotor configuration, and ranges from ~42.3 to ~51.7 kW, i.e., the EV with the considered set of motors will have very different longitudinal acceleration performance characteristics; (iv) at , the top power capability is provided by R-SOW, followed by R-IOW; (v) the usable torque per magnet volume, i.e., / , is the highest for R-SOW at and ; and (vii) R-SOW has the highest maximum efficiency at , while H-NOW has the highest maximum motor efficiency at . In general, based on the table and Figures 6 and 7, the R-SOW and R-IOW configurations, followed by H-NOW, appear to have overall desirable characteristics in terms of torque density production and premium efficiency area, which would make them good candidates for EV implementation.
To complete the analysis, Figure 8 shows the steady-state temperature distributions in the rotors, based on the FEA model results at the peak torque and base speed, and includes iso-temperature lines. The maximum temperatures are rather uniform among the magnet arrangements, with approximately 4 °C variations among the configurations. The lowest maximum temperature, equal to 64.5 °C, occurs for the H-NOW arrangement.

Case Study Vehicle
The case study EV is a light four-wheel-drive electric car, with two identical onboard IPMSMs, one per axle. The four-wheel-drive powertrain layout with central motors is becoming rather common in recent electric vehicles, see [14,43,44], as it conjugates a simple configuration of the individual powertrains (see [45] for an example of powertrain optimization study) with the enhanced performance of a four-wheel-drive system. Each IPMSM is connected to the two wheels of the axle through a single-speed mechanical transmission with open differential, half-shafts, and constant velocity joints, according to the schematic in Figure 9. The main vehicle parameters are in Table 6.

Case Study Vehicle
The case study EV is a light four-wheel-drive electric car, with two identical on-board IPMSMs, one per axle. The four-wheel-drive powertrain layout with central motors is becoming rather common in recent electric vehicles, see [14,43,44], as it conjugates a simple configuration of the individual powertrains (see [45] for an example of powertrain optimization study) with the enhanced performance of a four-wheel-drive system. Each IP-MSM is connected to the two wheels of the axle through a single-speed mechanical transmission with open differential, half-shafts, and constant velocity joints, according to the schematic in Figure 9. The main vehicle parameters are in Table 6.  Rolling resistance coefficient (factor to be multiplied by vehicle speed squared) s 2 /m 2 6.5 × 10 −6

Driving Cycle Simulation
By using the power loss maps in Figure 6, the motor performance for each magnet arrangement is analyzed in terms of energy efficiency during traction along multiple driving cycles, representative of a wide range of operating conditions, namely: (i) the New European Driving Cycle (NEDC); (ii) the Worldwide harmonized Light vehicle Test Procedure (WLTP) Class 2; (iii) the ARTEMIS urban cycle; and (iv) the ARTEMIS motorway cycle [46].
To emulate the operation along a driving cycle, the torque of each machine is calculated under the following assumptions: (i) the road has zero longitudinal gradient; (ii) the regenerative braking capability of the electric powertrains is neglected, i.e., only the friction brakes are used for slowing down the vehicle, and thus the efficiency performance of Figure 9. Schematic of the considered EV layout.

Driving Cycle Simulation
By using the power loss maps in Figure 6, the motor performance for each magnet arrangement is analyzed in terms of energy efficiency during traction along multiple driving cycles, representative of a wide range of operating conditions, namely: (i) the New European Driving Cycle (NEDC); (ii) the Worldwide harmonized Light vehicle Test Procedure (WLTP) Class 2; (iii) the ARTEMIS urban cycle; and (iv) the ARTEMIS motorway cycle [46].
To emulate the operation along a driving cycle, the torque of each machine is calculated under the following assumptions: (i) the road has zero longitudinal gradient; (ii) the regenerative braking capability of the electric powertrains is neglected, i.e., only the friction brakes are used for slowing down the vehicle, and thus the efficiency performance of the motors is analyzed solely for traction conditions; (iii) the transmission efficiency is constant; and (iv) the torque is evenly distributed between the two motors.
Each motor torque, T mot , is calculated through a backward facing approach: where R W is the wheel radius; η trans is the transmission efficiency; G r is the transmission gear ratio; and F t is the longitudinal force, given by the sum of the aerodynamic drag force, F d , the rolling resistance force, F r , and the inertial force, F a : where the individual force contributions are: In (17)- (19), ρ air is the air density; v is the vehicle speed; C r is the rolling resistance coefficient, expressed as C r = C r0 + Kv 2 , see [47]; g is the gravitational acceleration; m app is the apparent mass of the rotating components of the electric powertrains; and a x is the longitudinal vehicle acceleration. In the phases of the driving cycles in which v = 0, the condition F r = 0 is imposed.
Through (15)- (19), the operating motor profiles along the considered driving cycles, i.e., the time histories of motor torque and speed, are calculated and represented on the power loss map. The interpolation of the speed and traction torque values on the motor maps brings the total motor power loss profile in traction, W mot,loss (t), from which, through time integration, the motor energy loss along the driving cycle, E mot,loss , is obtained: where t 0 = 0 s and t N are the initial and final times of the driving cycle; the subscript i defines the considered time step; and the index N represents the last step of the driving cycle. The total motor energy consumption along the cycle, E mot,cons , is given by: where E mot,traction is the net traction energy of one motor, and ω mot = vG r /R W . The specific traction energy, E mot,traction , the specific energy loss, E mot,loss , and the specific energy consumption of the motor, E mot,cons , are obtained by dividing E mot,traction , E mot,loss , and E mot,cons by the driving cycle distance, d cycle : E mot,loss = E mot,loss d cycle (23) The average energy efficiency of the traction motors along the driving cycle, η mot , is calculated as:

Longitudinal Acceleration and Road Gradient Performance Simulation
Under the assumption of high tire-road friction condition, i.e., with friction coefficient µ max = 1.0, the vehicle performance associated with the considered motor arrangements is evaluated in terms of:

•
Maximum longitudinal acceleration from standstill, a x,max , on a road with zero longitudinal gradient: where the subscripts 'F' and 'R' refer to the front and rear axles; the notation F x,F/R indicates the front or rear longitudinal tire force in traction, calculated from the respective motor torque; and F z,F/R is the vertical load on the front or rear axle, which is a function of a x , see [47] for the details. According to (26), the optimization results are calculated under the assumption that the front-to-total motor torque distribution can vary, according to the intervention of a traction controller that prevents wheel spinning on the critical axle. The results showed that for the specific vehicle the intervention of the traction controller is not required for maximum longitudinal acceleration in high tire-road friction conditions. • Maximum longitudinal road gradient, ϑ max , on which the vehicle can travel at very low speed (approaching zero) and zero longitudinal acceleration: In this case, it was verified that the traction controller has to intervene on the front axle for most motor configurations, i.e., the motor torque on the front axle has to be reduced to the level compatible with the tire-road friction condition. Table 7 illustrates the motor and resulting vehicle performance characteristics. For the driving cycles, the table reports: (i) E mot,traction ; (ii) E mot,loss ; (iii) E mot,cons ; (iv) η mot ; (v) η mot,avg , i.e., the average efficiency of each motor configuration along the four cycles; (vi) σ η mot , i.e., the standard deviation of η mot of each configuration along the cycles. Low values of σ η mot highlight consistent energy efficiency for the considered variety of driving conditions; (vii) η DC , i.e., the average value of η mot of all motor configurations, computed for each driving schedule; and (viii) σ η DC , i.e., the standard deviation of η mot for all motor configurations, computed for each cycle. (i)-(vi) assess the efficiency performance of each motor configuration, while (vii) and (viii) provide an indication of the impact of the driving schedule on the results. Moreover, for the four cycles, the torque-speed diagrams in Figures 10-13 show the motor energy loss distributions for all cases. The white markers represent the operating points of the individual motors along the cycles. The magnitude of the energy losses in the corresponding torque-speed regions is indicated by the values in red as well as by the color maps, obtained by summing the energy losses of all operating points located in the region. Notes: Bold text indicates the best performance among the considered arrangements, whereas underlined text indicates the worst performance.

Vehicle Simulation Results
During the cycles, the motors mainly operate at low torque, i.e., at less than half of their maximum torque. The variety of driving conditions is confirmed by the spread of the E mot,cons and η mot values. For example, E mot,cons along the WLTP is less than half than along the ARTEMIS motorway. The energy efficiency along the ARTEMIS urban is the highest among the four driving cycles for each magnet arrangement. In general, motor efficiency is characterized by a decrease in high speed conditions. The effect is confirmed by η DC , which ranges from 0.9248 for the ARTEMIS motorway to 0.9520 for the ARTEMIS urban, indicating that the low speed operation typical of urban driving is favorable to efficiency. The highest significance of the magnet arrangement on energy consumption is obtained for the ARTEMIS motorway, with~3% energy efficiency increase from the R-IOW to the RCL-IOW configurations.
Among all arrangements, RCL-IOW provides the highest energy efficiency along each selected driving cycle, and has the lowest energy loss not only in the peak-loss region, but also in most operating areas, see Figures 10-13, followed by RCL-SOW and R-ROW. The RCL-IOW set-up also shows the lowest energy efficiency standard deviation among the different cycles, and thus provides consistently low consumption in various driving conditions.
At the bottom end of the ranking, R-SOW has the lowest energy efficiency along the NEDC and WLTP, while H-NOW and R-IOW rank last along the ARTEMIS urban and ARTEMIS motorway. The H-NOW arrangement, which was very promising according to the power efficiency map in Figure 7, is not competitive when realistic driving scenarios are investigated, although it exhibits high maximum efficiency for both the base and top speeds, see Table 5. On the contrary, the RCL-IOW arrangement, which in Table 5 has significantly lower maximum values of power efficiency than H-NOW, provides higher energy efficiency during driving cycle operation.
Interestingly, the results show that the observation of the power efficiency maps may be misleading for predicting the energy consumption performance along different driving conditions. This can be explained based on Figures 7 and 10, Figures 11-13, as most of the operating points of the considered driving cycles do not fall in the premium efficiency regions, but are rather at lower torque and higher speeds within the AP region, where the gradient of the motor efficiency characteristic, with respect to the variation of both speed and torque, becomes the most important aspect. The RCL-SOW and RCL-IOW configurations present indeed lower gradients of the power efficiency characteristic in comparison with the solutions with the widest premium efficiency regions, i.e., R-SOW and R-IOW.
The bottom half of Table 7 deals with the longitudinal acceleration and gradeability performance associated with each rotor configuration, and reports the values of a x,max and ϑ max . The most energy-efficient configurations, e.g., RCL-IOW, are characterized by low values of torque for the considered maximum inverter current, thus providing sub-optimal drivability at low speed. Vice versa, the configurations with the highest maximum torque score worse in the energy consumption simulations. For example, below the base speed, a layout like R-IOW provides a maximum torque that is~30 Nm higher than that for RCL-IOW, which is the most energy-efficient configuration. The correlation of the driving cycle results and electromagnetic analysis suggests that the low gradient of the efficiency map of the best performing motors along the cycles is associated with low flux densities, and thus low maximum torque capability. The important conclusion is that the rotor arrangement selection is a compromise between: (i) desirable acceleration and gradeability performance, requiring motor behavior optimized for high torque conditions, with high flux density; and (ii) energy efficiency along driving cycles, requiring low gradients in the relevant parts of the efficiency map, even at the price of narrower premium efficiency regions. the highest among the four driving cycles for each magnet arrangement. In general, motor efficiency is characterized by a decrease in high speed conditions. The effect is confirmed by , which ranges from 0.9248 for the ARTEMIS motorway to 0.9520 for the ARTEMIS urban, indicating that the low speed operation typical of urban driving is favorable to efficiency. The highest significance of the magnet arrangement on energy consumption is obtained for the ARTEMIS motorway, with ~3% energy efficiency increase from the R-IOW to the RCL-IOW configurations.
Among all arrangements, RCL-IOW provides the highest energy efficiency along each selected driving cycle, and has the lowest energy loss not only in the peak-loss region, but also in most operating areas, see Figures 10-13, followed by RCL-SOW and R-ROW. The RCL-IOW set-up also shows the lowest energy efficiency standard deviation among the different cycles, and thus provides consistently low consumption in various driving conditions.
At the bottom end of the ranking, R-SOW has the lowest energy efficiency along the NEDC and WLTP, while H-NOW and R-IOW rank last along the ARTEMIS urban and ARTEMIS motorway. The H-NOW arrangement, which was very promising according to the power efficiency map in Figure 7, is not competitive when realistic driving scenarios are investigated, although it exhibits high maximum efficiency for both the base and top speeds, see Table 5. On the contrary, the RCL-IOW arrangement, which in Table 5 has significantly lower maximum values of power efficiency than H-NOW, provides higher energy efficiency during driving cycle operation.
Interestingly, the results show that the observation of the power efficiency maps may be misleading for predicting the energy consumption performance along different driving conditions. This can be explained based on Figures 7 and 10-13, as most of the operating points of the considered driving cycles do not fall in the premium efficiency regions, but are rather at lower torque and higher speeds within the AP region, where the gradient of the motor efficiency characteristic, with respect to the variation of both speed and torque, becomes the most important aspect. The RCL-SOW and RCL-IOW configurations present indeed lower gradients of the power efficiency characteristic in comparison with the solutions with the widest premium efficiency regions, i.e., R-SOW and R-IOW.      The bottom half of Table 7 deals with the longitudinal acceleration and gradeability performance associated with each rotor configuration, and reports the values of , and . The most energy-efficient configurations, e.g., RCL-IOW, are characterized by low values of torque for the considered maximum inverter current, thus providing suboptimal drivability at low speed. Vice versa, the configurations with the highest maximum torque score worse in the energy consumption simulations. For example, below the base speed, a layout like R-IOW provides a maximum torque that is ~30 Nm higher than that for RCL-IOW, which is the most energy-efficient configuration. The correlation of the driving cycle results and electromagnetic analysis suggests that the low gradient of the efficiency map of the best performing motors along the cycles is associated with low flux densities, and thus low maximum torque capability. The important conclusion is that the rotor arrangement selection is a compromise between: (i) desirable acceleration and gradeability performance, requiring motor behavior optimized for high torque conditions, with high flux density; and (ii) energy efficiency along driving cycles, requiring low gradients in the relevant parts of the efficiency map, even at the price of narrower premium efficiency regions.

Conclusions
This study evaluated the influence of six buried magnet arrangements on the performance of a case study 48-slot 8-pole IPMSM for automotive traction applications, controlled through maximum-torque-per-ampere (MTPA) and maximum-torque-per-voltage (MTPV) strategies. The magnet sizing for each arrangement was based on the DOE/R sensitivity method, for a constant total magnet volume. The electromagnetic parameters for the rotor configurations were calculated through a two-dimensional finite element analysis approach, under the assumption of steady-state thermal conditions in the rotors. dqaxis modeling was adopted to analyze the performance of the IPMSMs with different rotor configurations, by using circuit equations that consider stator and PM flux linkages, as well as d-and q-axis inductances that are nonlinear functions of the flux level. The respective power loss and efficiency maps were generated over the entire torque-speed range. The performance of the considered IPMSMs was evaluated along four driving cycles, representative of urban and extra-urban driving conditions, as well as in terms of acceleration

Conclusions
This study evaluated the influence of six buried magnet arrangements on the performance of a case study 48-slot 8-pole IPMSM for automotive traction applications, controlled through maximum-torque-per-ampere (MTPA) and maximum-torque-per-voltage (MTPV) strategies. The magnet sizing for each arrangement was based on the DOE/R sensitivity method, for a constant total magnet volume. The electromagnetic parameters for the rotor configurations were calculated through a two-dimensional finite element analysis approach, under the assumption of steady-state thermal conditions in the rotors. dq-axis modeling was adopted to analyze the performance of the IPMSMs with different rotor configurations, by using circuit equations that consider stator and PM flux linkages, as well as dand q-axis inductances that are nonlinear functions of the flux level. The respective power loss and efficiency maps were generated over the entire torque-speed range. The performance of the considered IPMSMs was evaluated along four driving cycles, representative of urban and extra-urban driving conditions, as well as in terms of acceleration and gradeability, by simulating a light four-wheel-drive electric vehicle with two on-board motors and even wheel torque distribution among the two axles.
The main conclusions of the electric motor and vehicle simulation analyses are: • The power efficiency values show major variations for the different motor configurations, especially at medium-to-high speed values. The premium efficiency region, i.e., with efficiency values in excess of 0.97, covers a small range, both in terms of torque and speed, for the R-ROW, RCL-SOW, and RCL-IOW configurations, and a slightly wider speed range with a significantly wider torque range in the R-SOW, R-IOW, and H-NOW layouts.

•
For the considered maximum inverter current, the torque capability of the motor is significantly affected by the magnet arrangement; in fact, the maximum torque ranges from~130 Nm for the RCL-IOW configuration to~160 Nm for the R-SOW layout.

•
As the power efficiency of the case study motors tends to significantly decrease with speed, the cycles with substantial urban driving sections provide higher energy efficiency, e.g., the average energy efficiency of the considered motors is 0.9248 along the ARTEMIS motorway cycle, and 0.9520 for the ARTEMIS urban.

•
The energy efficiency along driving cycles is affected by the magnet arrangement, and can vary up to~3% depending on the rotor configuration. The WLTP class 2 cycle originates the maximum standard deviation of the motor energy efficiency with the considered set of rotor layouts. The magnet configuration selection is crucial for traction motors of electric vehicles operating in a wide range of speeds, while it is less significant in urban driving, e.g., the standard deviation of the energy efficiency along the ARTEMIS urban is less than one third of the respective value along the ARTEMIS motorway.

•
In typical driving schedules, the motors operate at relatively low torque, i.e., normally at less than half of their peak torque in the considered electric vehicle. Hence, the magnet arrangements providing high values of maximum efficiency for high torque at the base and top speeds, such as H-NOW, may not yield practical benefits along normal driving, while configurations with lower maximum efficiency at high torque, such as RCL-IOW, can be rather effective in realistic scenarios.

•
With an average efficiency of 0.9466, the RCL-IOW layout is the most efficient configuration for the specific vehicle, followed by RCL-SOW and R-ROW, having average energy efficiencies of 0.9461 and 0.9375 along the selected schedules.

•
For the given maximum inverter current, the longitudinal acceleration and gradeability performance is very strongly affected by the rotor arrangement, i.e., the maximum longitudinal acceleration at zero road gradient ranges from 5.61 to 6.87 m/s 2 , and the maximum road gradient ranges from 34.7 to 37.5 deg, in high tire-road friction conditions.
Interestingly, the rotor configurations providing the best drivability performance are those with the lowest energy efficiency, and vice versa. Hence, the selection of the motor rotor configuration must be the result of a trade-off between drivability and energy efficiency, depending on the specific vehicle requirements, and cannot be based only on motor performance analysis results.
Future investigations will include: (i) more accurate modeling techniques for the machines, e.g., three-dimensional FEA approaches accounting for the axial length effects; (ii) optimal sizing of the magnets for each rotor arrangement; (iii) more advanced vehicle simulations, e.g., to include tire slip power losses; (iv) energy-efficient front-to-rear motor torque distribution algorithms, including regenerative braking functionalities; and (v) different electric powertrain topologies.
Author Contributions: For this research article, P.A., R.B., A.L. D.T., L.S. and A.S. have contributed to the investigation, methodology, resources, writing and editing, and supervision. The software simulations and visualizations were mainly carried out by P.A. and L.S. All authors have read and agreed to the published version of the manuscript.