Comparative Study of Non-Rare-Earth and Rare-Earth PM Motors for EV Applications

: Recently, non-rare-earth motors are attracting more and more attention due to the booming of the electric vehicle (EV) market and, more importantly, the increasing price of the rare-earth magnet material. This paper focuses on the performance comparison among a commercial interior permanent magnet (IPM) motor and two non-rare-earth motors, including a synchronous reluctance motor (SynRM) and a permanent-magnet-assisted synchronous reluctance motor (PMaSynRM). The design procedure to develop a high-torque-density, low-torque-ripple and high-efﬁciency SynRM is presented. Combined with a developed automatic modeling and simulation procedure, the ﬁnite element analysis (FEA)-based differential evolution (DE) algorithm is introduced for the SynRM rotor optimization. In order to fully inspire the potential of the SynRM, a novel method to optimize the motor split ratio is proposed under the constraint of the copper loss. In addition, different slot–pole combinations are investigated to maximize the motor torque, and the rotor structure is also dealt with towards the centrifugal stress at the maximum operating speed. Finally, the motor performance comparison is carried out, and the results show that although the SynRM achieves almost 61% cost savings, its poor torque capability, power factor and ﬂux weakening (FW) capability are non-negligible defects. On the contrary, the PMaSynRM exhibits excellent features for the EV applications in terms of cost, torque density, efﬁciency and FW capability. This paper presents a novel split ratio optimization method for the optimal SynRM/PMaSynRM design and demonstrates the characteristics of the IPM motors, SynRMs and PMaSynRMs for EV applications.


Introduction
With the growth of global warming and air pollution, electrical traction systems as gaining more and more interest in both the academic and industrial fields. It is remarkable that the popularity of electric vehicles (EVs) and hybrid electric vehicles (HEVs) has expanding significantly in the last decade. The electrical traction motor represents the heart of the energy conversion process. Among various machine topologies, the interior permanent magnet (IPM) motor has been recognized as the most promising candidate [1][2][3]. On the other hand, some challenges exist in IPM motors, i.e., high production cost and high risk of irreversible demagnetization. Most importantly, a worldwide concern over the price volatility and availability of rare-earth magnets is growing. Therefore, many researchers are focused on alternative non-rare-earth traction drives [4,5].
The induction motor (IM) firstly comes to the researchers' minds, mainly due to its low cost, robust structure and well-established manufacturing techniques [6,7]. The wide flux-weakening range makes it very suitable for EV applications. Nevertheless, the low power density and inferior efficiency are potential challenges of IMs [4]. The switched reluctance machine (SRM) is also a good candidate in terms of simple and robust rotor The investigation of SynRM or PMaSynRM for EV applications has been carried out in literature. A comparative study of SynRM, PMaSynRM and the benchmark Prius 2010 IPM motor was conducted in [25], in terms of electromagnetic performance and material cost. The slot-pole combination and copper loss are kept constant for these designs. The SynRM design utilizing dual-phase magnetic materials can be competitive with an IPM motor when considering the power density and efficiency [26]. A carbon-fiber-wrapped SynRM was designed and tested according to the US DOE Freedom Car 2020 specifications, which aims to reduce the flux leakage and improve the saliency ratio [27]. A novel magnet-axis-shifted hybrid PM machine was proposed for EV application in [28]. It achieved a higher torque capability, wider high-efficiency operation range but higher cogging torque and iron losses than that of the benchmark BMW i3 IPM machine. In addition, the comparison of the ferrite and rare-earth-based PM machines was summarized in [22], and both the advantages and disadvantages were highlighted.
In [15], an automatic modeling and simulation procedure of SynRM was proposed, considering different flux barrier numbers and barrier end angles. This procedure is then extended to be coupled with the differential evolution (DE) algorithm to achieve an automatic SynRM rotor optimization. By using this tool, this paper aims to design a SynRM and a PMaSynRM according to the dimension of Lexus LS 600h IPM motor, so as to achieve a detailed performance comparison. The following sections will be organized as follows: The SynRM rotor optimization process using finite element analysis (FEA)based DE algorithm is introduced in Section 2. After that, the design procedure for the whole machine geometry, which includes split ratio optimization, slot-pole combination selection, rotor optimization and rotor structure improvement, is proposed in Section 3. Then, the performance comparison between the designed SynRM/PMaSynRM and the benchmark IPM motor is carried out in Section 4, and finally, the conclusions are drawn in Section 5.

SynRM Rotor Optimization with DE Algorithm
The Lexus LS 600h IPM motor is selected as the reference in this paper, with its performance deeply evaluated and tested in [29]. Main dimensions of the Lexus LS 600h IPM motor are listed in Table 1. These dimensions are the baselines of the SynRM design. It is noticed that the Lexus LS 600h IPM motor achieves 110 kW within the stator temperature limit of 150 • C for a transient time of 18 s when the current is 330 A peak at the speed of 4500 rpm [29]. Then, the current limit is set to be 330 A peak , and the calculated current density and stator fill factor are 25 A/mm 2 and 0.42, respectively. This operation point is chosen as the reference for the following design. In order to present a fair comparison with the Lexus LS 600h IPM motor, the SynRM is constrained to satisfy some dimension limits of the benchmark motor. Firstly, the rotor is replaced by a SynRM topology, while the stator remains the same. For the design of SynRM, the most challenging and creative aspect lies in the rotor design by dramatically increasing the rotor's anisotropy to improve the torque density and selecting proper barrier end angles for torque ripple suppression. Key efforts to achieve these objectives include parametric machine modeling and optimization strategy. In [15], an automatic modeling and simulation procedure has been developed. Hereafter, it is combined with differential evolution (DE) algorithm to achieve the optimal SynRM rotor design.

DE Algorithm
As is well-known, the DE algorithm is a population-based stochastic searching method, suitable for combinatorial optimization problems [30][31][32]. It was firstly proposed in the mid-1990s and has been applied for machine optimization problems in recent years [33,34].
The algorithm includes the processes of initialization, mutation, crossover and selection. At the beginning, the initial parameter vectors of the individuals are randomly generated, which covers the entire parameter space. After evaluating the fitness of each member, parameter vectors of a new generation are created through the operators of mutation and crossover. Then, an explicit one-to-one survivor selection is carried out. Once the new population is generated, the process of mutation, crossover and selection is repeated until the convergence of the algorithm or a specified number of generations is reached. This procedure is illustrated in Figure 2, and detailed explanations are referred to in [30,31]. Machine modelling, simulation and post-processing to obtain the objectives.
Find Pareto front of the individuals.
Creat the individuals of new generation.
Initialization of the individuals randomly.

No
Limit the parameters within the bounds. The process of mutation expands the search space by using a scale factor F (i.e., mutation ratio) to recombine a random vector differential with the base individual. A high value of F helps to avoid the local optimum, while it simultaneously sacrifices the speed of convergence. Therefore, the value of F is seldom greater than 1.0. The process of crossover, which is performed with a probability of CR (i.e., crossover ratio), mixes successful solutions from the previous generation with current donors. In the extreme circumstance, CR = 1 indicates a full consideration of the current donors, which undoubtedly increases the population diversity but makes the difference between the donor vector and the previous vector larger. The rule of thumb values of the mutation and crossover ratios is recommended as F ∈ [0.5, 1.0] and CR ∈ [0.5, 1.0] [31]. In this problem, they are set to be F = 0.9 and CR = 0.7, respectively.

SynRM Rotor with Three and Four Flux Barriers
The preliminary design of the SynRM is carried out by fixing the same stator as the Lexus LS 600h IPM motor. The rotor pole number and airgap thickness remain the same as the Lexus LS 600h IPM. The DE algorithm is combined with finite element analysis (FEA) for the rotor optimization, which aims to have a low-torque-ripple and high-torque design. The optimization procedure is illustrated in Figure 2.
As one of the key issues for geometry optimization, SynRM rotor parametrization usually presents many degrees of freedom. The automatic drawing procedure proposed in [15], however, requires only two main parameters for SynRM rotor modeling, i.e., insulation ratio k air and barrier end angle θ b . The streamlined flux barrier shape is adopted, which is set to be parallel to the d-axis flux lines and perpendicular to the q-axis flux lines when saturation is neglected. This allows one to theoretically maximize the d-axis inductance and minimize the q-axis inductance, thus achieving a high saliency ratio and high torque. Points on the q-axis, as shown in Figure 3, are selected as starting points of the flux barrier boundaries, and their coordinates are computed through the width of each flux barrier t bi and the width of each iron path w f ei . By introducing the insulation ratio and the barrier end angles, w f ei is calculated according to the d-axis flux density distribution, and t bi is obtained based on the q-axis magnetic voltage drop across the corresponding flux barrier. Related equations and some tuning tricks to overcome geometric issues can be referred to [15]. As a result, a robust and fast modeling framework for SynRM rotor design is created. Since this automatic modeling procedure is developed in a FEA tool, machine performance is then obtained by FEA simulation and post-processing. Therefore, the SynRM rotor optimization can be achieved by combining the automatic modeling and simulation procedure with DE algorithm. The evaluation of individuals is completed by calling the FEA tool and finding the Pareto front, which are described in Figure 2. As previously mentioned, the optimization parameters are the insulation ratio k air and the barrier end angles θ b . In fact, the number of barrier end angles depends on the number of flux barriers per pole. For example, the parameters definition of a three-barrier SynRM rotor is shown in Figure 3. The optimization objectives are focused on average torque and torque ripple, aiming at a low-torque-ripple and high-torque design. The torque ripple is defined as the ratio of peak-to-peak value to the average torque.
According to the previous investigation [15], three and four flux barriers per pole are studies in order to balance the output torque and manufacture difficulties. The optimization result for SynRM with three barriers per pole is plotted in Figure 4a, and the corresponding result of the four barriers per pole motor is shown in Figure 4b. Individuals are marked with dots on the average torque and torque ripple plane, where different colors indicate different generation numbers. For each individual, the current angle is tuned to achieve the maximum torque within the current limit. It is notable that the sign of the average torque is reversed because the aim of this optimization is to minimize each objective. The better solutions move towards the bottom left of the plane. The black diamonds form the Pareto front, which views clearly how each objective is penalized by the improvement in the other ones. The preferred compromise solution can be selected among the ones on the front.
The results show that the torque ripple of the SynRM can be significantly reduced by using the optimization algorithm, irrespective of the number of flux-barriers. The maximum torque of the four-barrier machine is slightly higher than that of the three-barrier machine, which are 167.5 Nm and 165.5 Nm, respectively. These values account for about 70% of the Lexus LS 600h IPM motor. For the next step, the machine torque will be further improved by redesigning the whole machine geometry. Only the three barriers per pole rotor is considered in the following analysis.

Reconsideration of the SynRM Design
To improve the performance of the SynRM, the stator is redesigned with respect that of the Lexus LS 600h IPM motor. The design has several degrees of freedom but still has some constrains. In particular, the stator's outer diameter, stack length, peak current and copper loss constraints are the same as the Lexus LS 600h IPM motor. Hereafter, the optimization procedure is discussed. It is notable that the optimization procedure starts with the 48s-8p configuration, referring to the Lexus LS 600h IPM benchmark.

Split Ratio Optimization
The split ratio (χ) is defined as the ratio between the stator's inner and outer diameters, i.e., χ = D/D e . The torque of an electric machine can be expressed as [35]: where L stk is the stack length,B g is the amplitude of the fundamental harmonic of airgap flux density, andK s is the linear current density, which is represented as: where k w is the winding factor, Q s is the number of stator slots, n cs is the number of series conductors per slot andÎ is the phase current peak value. Therefore, the torque equation can be expressed as: It is observed that the output torque is affected by the split ratio when the stator outer diameter D e is fixed. The focus of this optimization lies in finding the optimal split ratio that achieves the highest torque.
In this analysis, the rotor parameters (flux barrier number, barrier end angles and k air ) remain the same. The stator slot is simplified to a trapezoidal shape, as shown in Figure 5. The slot opening is calculated from the slot pitch and stator tooth width, while the slot wedge is neglected. The flux density ratio β is defined as the ratio of airgap flux densityB g and stator teeth flux density B t . Then, the stator teeth width w t is computed as [36]: This equation is also valid for non-linear conditions only if the leakage flux is neglected. The flux passing through the stator's back iron accounts for half of the main flux in the air gap, which is [36]: where B bi is the flux density of the stator's back iron, p is the number of pole pairs. Combing (4) and (5), the stator's back iron height h bi is then calculated as: Here, another parameter α is introduced, which is defined as α = B t /B bi . Considering the knee point of the iron lamination, this ratio is assumed to be constant and equal to 1.16. Using Equations (4) and (6), the slot area S slot is represented as: where f a , f b and f c are: As a result, S slot is simply related to β and χ, since α, D e , Q s and p are supposed to be constant for a given slot-pole combination. The following step is to determine the variation in the slot area S slot with split ratio χ, which is based on the assumption of a constant copper loss constraint.
The copper loss can be calculated as: where ρ cu is the copper resistivity, L cond is the copper conductor length, k f ill is the slot fill factor and m is the number of phases. L cond includes the stack length and the end-turn length. The end-turn length is technically simplified to 1.6 times of the pole pitch. Then, L cond is also related to the split ratio: From (9), all the parameters are constant, except n cs , χ and S slot . Therefore, the variation in the slot area S slot with χ can be directly obtained once n cs is given. The relationship of S slot , n cs and χ for a 48-slot 8-pole SynRM is plotted in Figure 6. The results show that the slot cross-area increases with the increase in the split ratio χ or n cs . After determining S slot , (7) and (8) are combined to calculate β, which results in: where: Then, w t and h bi are computed from (4) and (6), respectively. As a consequence, all the parameters to describe the stator's geometry, i.e., w t , h bi and D, are determined. By using the automatic modeling and simulation procedure, it is easy and fast to repeat these calculations and to investigate the torque behavior for different split ratios. It is noted that the torque is obtained through FEA simulations, and iron saturation is also included. In addition, the current angle is also varied to achieve the maximum average torque values.
The variation of torque with a split ratio is shown in the upper plot of Figure 7, while the corresponding current density is plotted below. Obviously, there exists an optimal split ratio to achieve the highest torque, which is simultaneously affected by n cs . On the other hand, the split ratio is constrained by the current density limit (shown by the dotted line in the lower plot of Figure 7), which is fixed to 22.8 A/mm 2 . This value is reduced from 25 A/mm 2 of the Lexus LS 600h IPM motor due to the increased slot area in the simplified trapezoidal slot. With the decrease in n cs , the current density tends to overcome the limit. Consequently, the calculation of the current density is required during the split ratio optimization process, and the candidates with current densities beyond the limit should be discarded. The maximum torque obtained from Figure 7 is 170.9 Nm, with n cs = 4 and χ = 0.65.

Slot-Pole Combinations
For the SynRM, it has been revealed that the saliency ratio and torque capability highly depend on the number of pole pairs [35]. Therefore, the slot-pole combination should be reconsidered for maximizing the machine torque. In the following, several combinations are selected and analyzed as listed in Table 2. For each slot-pole combination, the split ratio optimization is carried out with the aforementioned procedure. After plenty of calculations and FEA simulations, the maximum torque for each slot-pole combination is again reported in Table 2. The optimal n cs and χ are also shown. It is interesting to notice that the optimal split ratio is almost constant corresponding to the same number of pole-pairs. It is increased as the pole-pair number is increased, so as to reduce the tooth length and limit the slot leakage component. A similar conclusion was drawn for the PMaSynRM in [24]. The optimal n cs is found to be inversely proportional to the slot number due to the constant copper loss constraint. In addition, the highest torque is achieved when the pole-pair number is three instead of two, which does not match with the previous finding that the higher pole number the lower torque [37]. In fact, the machine saliency ratio is dramatically reduced as the magnetic saturation increases under overload condition. Therefore, the merit of a thinner stator back iron for higher pole number is becoming more significant.
The results show that the slot-pole combinations with q = 3 exhibit the highest torque capability. The torque behavior versus the split ratio with the optimal n cs of these combinations is plotted in Figure 8. As mentioned, the optimal split ratio increases with the number of pole pairs. The maximum torque of the 54s-6p combination, though not significantly improved from the 72s-8p combination, is the highest. In addition, a low number of pole pairs reduces the operating frequency, leading to low iron losses in the stator. Overall, the 54s-6p combination is chosen as the optimal competitor. As concluded in [35], the saliency ratio of the SynRM is increased as the pole pair number is reduced. Therefore, it is reasonable to have a lower number of pole pairs in a SynRM design compared with an IPM motor. Since this optimization process is based on a trapezoidal shape slot, the slot opening and wedge are then added to the stator's geometry (the same values as the Lexus LS 600h IPM motor). The slot area is no longer the same, and β also has to be changed to achieve the same copper loss. Figure 9 shows the torque versus split ratio comparison between the ideal "trapezoidal-shape" slot and the "actual-shape" slot (including wedge and slot opening). As expected, the trend of torque variation is almost the same. Finally, the optimal split ratio is selected to be 0.6. . Torque versus split ratio comparison between the SynRMs with simplified "trapezoidalshape" slot and "actual-shape" slot.

Rotor Optimization
The rotor optimization process is the same as that described in Section 2. Additionally, efficiency is considered as the third objective, which is obtained from the loss calculation. The copper loss is always the same, since it is kept constant. For simplicity, the iron losses are calculated as follows: P f e = k magg t p sp t G t + k magg bi p sp bi G bi (13) where k magg t and k magg bi are manufacturing loss increase factors in teeth and back iron, respectively. G t and G bi are the weight of teeth and back iron, respectively, and the specific losses (in teeth and back iron, respectively) are computed as: where the specific iron loss of iron lamination p sp f e is the value at reference flux density B re f and reference frequency f re f . k hy and k ec are the iron hysteresis loss coefficient and the eddy current loss coefficient. f is the operating frequency, which is 225 Hz at the speed of 4500 rpm. B refers to the maximum flux density of the teeth or back iron, which are obtained from the FEA simulations. The stray loss is considered as 10% of the total losses (P cu + P f e ). As before, the operation point for the optimization is under the maximum current limit and base speed. After optimization, a three-dimensional view of the result is shown in Figure 10. It is interesting that the variation in the torque ripple for the solutions on the Pareto front is more significant than the variations of average torque and efficiency. The optimal point is then selected to achieve the minimum torque ripple, with slight penalization on average torque and efficiency. The projections of this optimal point on each plane are drawn by the black points in Figure 10. The rotor parameter values of the optimal SynRM are shown in Table 3.

Rotor Structure Enhancement
Since the Lexus LS 600h IPM motor is operated as high as 10,230 rpm [29], the rotor structure analysis of the designed SynRM is conducted. Figure 11a shows the rotor mechanical stress distribution towards centrifugal stress. The maximum stress is noticed to be 1461 MPa, while the Yield stress of the iron lamination is only around 440 MPa. As a result, central bridges in the second and third flux barriers are introduced to reduce the stress concentration on the outer ribs. The thicknesses of these central bridges are optimized to 0.3 mm and 1.0 mm, respectively. As shown in Figure 11b, the maximum stress is then reduced to 354 MPa, which allows a safety factor of 1.24. As a sacrifice, the machine torque is slightly reduced due to the leakage flux path created by the central bridges. The cross-section of the designed SynRM is shown in Figure 1.

Performance Comparison among the Motors
According to [29], the maximum torque of the Lexus LS 600h IPM motor within the temperature limit of 150 • C for a transient time of 18 s is around 233 Nm. The maximum torque of the optimized SynRM under the same peak current is 178 Nm, which accounts for only 76%. Therefore, the stack length of the designed SynRM is increased in order to meet the same torque and power requirements. As described in the literature [23,24], the PMaSynRM can be obtained directly from the optimized SynRM considering the maximized reluctance torque. For simplicity, the ferrite magnets are fully inserted into the flux barriers of the SynRM rotor, which achieves the same torque as the benchmark IPM motor. The residual magnetic flux density of the ferrite magnet is 0.4 T, and the relative permeability is 1.05. In the following, a detailed performance comparison between the designed motors and the benchmark motor is presented. Table 4 summarizes the main design characteristics and performance of the designed SynRM/PMaSynRM and the benchmark IPM motor. The active stack length of the SynRM is increased from 135.4 mm to 182 mm, which leads to higher volume and mass. The copper consumed in the SynRM and PMaSynRM is considerably higher than that in the IPM motor, which is attributed to the higher slot area and the longer end winding due to the lower number of pole pairs. The magnet mass of the PMaSynRM is almost doubled from the benchmark IPM motor because of the low residual flux density of the ferrite magnet. It is also noticed that the SynRM requires a deeper rotor to allow for flux flow and a high saliency ratio, and thus, the shaft diameter is lower than the benchmark IPM motor. A shaft diameter of 44 mm is safe and refers to the 2012 Nissan Leaf motor [38].

Design Features
The unit prices for iron lamination, copper, ferrite and NdFeB magnet are assumed to be 1.4 $/kg [39], 9.8 $/kg [39], 3.2 $/kg [40] and 126.2 $/kg [40], respectively. Both of the SynRM and the PMaSynRM achieves significant cost reduction when compared to the benchmark IPM motor, which are 61.0% and 65.4%, respectively. The rare-earth magnets of the IPM motor contribute over 70% of the total cost, though they account for only a small fraction of the total mass (less than 5%). Given the low cost of the ferrite magnet, the PMaSynRM exhibits the best balance between machine volume and cost.

Torque and Power Capabilities
Generally speaking, the maximum torque and power capabilities of a traction motor are limited by both inverter maximum Volt-Amps rating and machine temperature rise. According to [29], the phase current peak value is set to be 330 A, and the base speed is 4500 rpm with a DC-bus voltage of 650 V. Due to the low power factor of the SynRM, either the peak phase current or DC-bus voltage should be increased, which can be achieved by adjusting the number of series conductors per phase. Since both motors show the same machine performance, the case with higher DC-bus voltage (874 V) but the same peak current is adopted.
As shown in Table 4, the performance of the three motors operated with the peak current at base speed is compared. For a fair comparison, the SynRM and PMaSynRM are designed to exhibit similar torque and power values as the benchmark IPM motor. As expected, the power density and specific power are lower for the designed SynRM. The other problem with the SynRM is the low power factor, which is only 0.58, as reported in Table 4. As a comparison, the power factor of the benchmark IPM motor is 0.78, since it is under the overload condition. The low power factor of the SynRM increases the inverter rating, which explains the increase in the DC-bus voltage. On the contrary, the PMaSynRM achieves almost the similar power density and specific power as the IPM benchmark.
The torque waveforms of the three motors are compared in Figure 12. The average torque of the benchmark IPM motor is very similar to the experimental results in [29], which validates the FEA simulations in this paper. Although the FEA-combined DE algorithm has been adopted for optimization, the torque ripple of the designed SynRM and PMaSynRM are both higher than that of the benchmark IPM motor. The torque ripple is slightly increased from 5.7% to 7.8% by introducing the magnets into the rotor flux barriers. The torque ripple of the PMaSynRM can be further reduced by skewing or introducing minor slits in the rotor surface.

High-Speed Flux-Weakening Capability
As far as the high-speed flux-weakening (FW) capability is considered, some control methodologies are adopted to achieve the best performance. The optimal current vector trajectory over base speed is achieved by means of FW control and maximum torqueper-voltage (MTPV) control [35]. When the inverter voltage limit is reached, the d-axis current is reduced to weaken the stator flux, i.e., the flux-weakening region. In some cases, for example, when the PM flux linkage is quite low, the current limit circle, the voltage limit ellipse and the MTPV trajectory intersect at a certain speed. Then, the motor moves to the MTPV region at higher speed, in order to maximize the torque under current and voltage limit. In this region, the current amplitude decreases with the increasing speed until reaching the maximum speed. The obtained current vector trajectory is plotted in Figure 13. Within the speed limit of 10,230 rpm, the IPM motor follows only FW operation, while the operations of the SynRM and PMaSynRM cover both FW and MTPV.
The comparison of torque and power capabilities at different speeds of the three motors are plotted in Figure 14. The dashed lines represent the characteristics of the benchmark IPM motor, while the solid lines show the performance of the designed SynRM and the dotted lines refer to the PMaSynRM. The results show that the power of the SynRM drops considerably and the power difference between the SynRM and the IPM motor increases with the speed. More specifically, the benchmark IPM motor delivers more than 120 kW at 10,000 rpm, whereas the SynRM outputs less than 60 kW. Therefore, the SynRMs are not suitable for applications that require a wide constant power speed range (CPSR). The high-speed FW capability of the SynRM can be improved by insetting PMs into flux barriers, as demonstrated by the dotted lines in Figure 14. As concluded in [15], the wider the PM width the better the FW capability.

Efficiency
The efficiency maps of the three motors are compared in Figure 15. They are obtained by considering copper loss, iron loss and stray loss. The effect of temperature rise on the stator resistance is neglected. From these plots, some useful conclusions can be drawn.
The maximum achievable efficiency of the PMaSynRM is higher than that of the benchmark IPM motor, which is as high as 96%. In addition, a large portion of the operation region achieves an efficiency above 90% in the PMaSynRM. It is also worth noticing that the efficiencies of both SynRM and PMaSynRM are higher than that of the IPM counterpart at high-speed low-torque conditions, which is mainly due to the lower flux density of the SynRM and PMaSynRM at this condition.

Conclusions
This paper presents a comparative study among a commercial IPM benchmark, a SynRM and a PMaSynRM. Simply starting from rotor replacement, the SynRM rotor is designed by an optimization procedure using FEA-based DE algorithm. It is shown that the achievable torque of the SynRM accounts for 70% of the IPM motor. To fully inspire the potential of the SynRM, the design procedure, which consists of split ratio optimization, slot-pole combination selection, rotor optimization and rotor structure improvement, is followed. A novel method to obtain the optimal split ratio and slot-pole combination with the copper loss constraint, which takes into account the saturation, is proposed. It is noticed that the optimal pole number for the SynRM is six instead of eight for the IPM motor, which implies that the slot-pole combination should be redesigned when machine topologies are changed.
The performance comparison between the optimized SynRM/PMaSynRM and the commercial benchmark is also carried out. The results show that the SynRM achieves the characteristics of significant cost saving (61% cheaper) and low torque ripple (5.7%). The low power factor reveals increased inverter rating and cost for the SynRM, and the poor FW capability limits its potential for wide CPSR applications. With added ferrite magnet in the rotor flux barriers, the torque density, power factor, efficiency and FW capability are improved. The PMaSynRM enables to achieve a good balance between the cost and the electromagnetic performance. Although the ferrite magnets may experience the risk of demagnetization, it can be well addressed by fine designs with the flux barrier. This paper clearly shows the gap between well-optimized SynRM/PMaSynRM and a commercial IPM motor. The SynRM can potentially be adopted in cost-effective specialty vehicles, while the competitive performance as the rare-earth PM machines makes the PMaSynRM very promising for passenger cars, buses and heavy-duty trucks.