Five-Phase Modular External Rotor PMMachines with Different Rotor Poles : A Comparative Simulation Study

The performance of fault-tolerant modular permanent magnet machines depends on the proper selection of the pole and slot numbers which result in negligible coupling between phases. The preferred slot and pole number combinations eliminate the effect of low-order harmonics in the stator magnetomotive force and thereby the vibration and stray loss are reduced. In this paper, three external rotor machines with identical machine dimensions are designed with different slots per phase per pole ratios. A simulation study is carried out using finite element analysis to compare the performance of the three machines in terms of machine torque density, ripple torque, core loss, and machine efficiency. A mathematical model based on the conventional-phasemodel approach is also used for the comparative study. The simulation study is extended to depict machine performance under fault conditions.


Introduction
Thanks to their high power density, fault tolerant permanent magnet brushless drives have shown promise in critical applications, such as in the electric vehicles (EV), aerospace, and automotive sectors [1].Multiphase modular permanent magnet machines comprise m magnetically and physically isolated phases, each having one-per-unit inductance which limits the short circuit current to its rated value under fault conditions.Hence, a fault in any phase will not affect others.Many papers address the design and the operating principle of fault-tolerant permanent magnet (PM) machines [2,3].
These features have made the multiphase modular PM machine a strong candidate in the field of fault tolerant applications [4][5][6].Optimum control under healthy and faulted scenarios is proposed to minimize the machine torque ripple [7,8].Favorable slot and pole number combinations for fault-tolerant PM machines have been proposed in [3], which ensure inherently negligible coupling between phases.For a five-phase modular machine, there are three favorable slot and pole number combinations [3], namely, 20 slots/14 poles, 20 slots/18 poles, 20 slots/22 poles and their multiples.In [3], the first two combinations are compared in terms of flux distribution and core loss.In [9], a machine with the second combination is designed by means of the star of slots [10] and using a doublelayer fractional-slot winding.Current control strategies, in case of one-phase or two-phase faults are proposed.An analytical model has been adopted to individuate the most suitable current references, without increasing the motor or inverter cost.In [11]; a comparison between two 5phase PM machine winding configurations is carried out.A new fractional-slot winding, which drastically decreases the cogging torque, is proposed and compared with a classical full-pitched concentrated winding.This method was applied to improve the matching between a naval propulsion fivephase PM machine and its converter.The proposed winding configuration significantly improved torque quality (higher density and lower ripple) and simplified current control without oversizing the electronic components.
In this paper, a complete comparison between the three favorable combinations for five-phase modular machines is carried out in terms of machine torque density, torque ripple, core loss, and overall efficiency.The comparison is based on a simulation study using FEA.Additionally, a mathematical model based on the conventional-phase-model approach [12] is used to corroborate the results.Three machines with an external rotor [13] are designed with the same stator and

Machine Design
Three modular PM machines have been designed and modeled using finite element analysis (FEA).The magnets are NdFeB with B r = 1.1 Tesla.The stator has five phases, each comprising two coils with 87 turns per coil.The winding configuration from [3] has been used and is shown in Figure 1(a).It is a 20 stator slot single-layer winding.The stator outer diameter, shaft diameter, rotor inner diameter, rotor outer diameter, and stack length are 260 mm, 144 mm, 272 mm, 312 mm, 100 mm, respectively.The magnet thickness is 5 mm and the air gap length is 1 mm.The proposed power rating is 8 kW at a speed of 240 rpm. Figure 1(b) shows the layout of the three machines.

Favorable Slot/Pole Combinations and Corresponding Stator Phase Sequence
It has been shown in [3] that certain slot and pole number combinations give negligible coupling between phases.For any modular winding with m phases and S slots, the required number of poles is given by (1), and the resulting phase angle between adjacent coils is given by (2); where n = 1, or n = any nonzero odd integer less than m, such that n and m do not share any common factors.The variable n is related to the applied phase sequence number as will be shown.For m = 5, the available values for n will be 1 and 3. Substituting in (1) and assuming that the number of slots is 20, the corresponding number of poles will be as follows: (i) for n = 1 → 2p = 18 or 22 and the corresponding angle between phases will be 36 • and −36 • , respectively; (ii) for n = 3 → 2p = 14 or 26 and the corresponding angle between phases will be 108 • and −108 • , respectively.
To prove the following relations, the stator MMF should be analyzed to investigate the harmonic content of the resulting MMF for a given winding current.For the winding layout shown in Figure 1 [3], the corresponding winding function for any phase is shown in Figure 2(a), and the corresponding harmonic spectrum is shown in Figure 2(b).
For a multiphase system with m phases, the corresponding available current sequences also correspond to m.This includes one or two zero-sequence vectors for odd and even numbers of phases, respectively.For a five-phase system, there will be five (4 + 1 zero) available sequences.The fourth sequence is the negative sequence of the first, while the second sequence is the negative sequence of the third.The angle between phase currents in the fundamental sequence is α = 2π/m.The stator MMF corresponding to currents with fundamental sequence is shown in Figure 3(a).The corresponding harmonic spectrum is shown in Figure 3(c).The fourth sequence, which corresponds to an angle 4α between phase currents, will give rise to the same flux distribution as the fundamental sequence; however, the MMF wave travels in the reverse direction.For the third sequence, corresponding to an angle 3α between phase currents, the MMF distribution is shown in Figure 3(b), and the corresponding harmonic spectrum is shown in Figure 3(d).The second sequence gives rise to the same flux distribution as the third, with a reverse travelling MMF.
Comparing the two cases, it is evident that the fundamental sequence produces the following space harmonics (3, 7, 13, 17, 23, 27, etc.), while the third harmonic sequence produces the following harmonics (1, 9, 11, 19, 21, 29, 31, etc.).Combining the two harmonic spectra will result in all harmonics except the fifth and its multiples (5, 10, 15, etc.).Table 1 summarizes the results obtained up to the 21st harmonic.The air gap flux densities due to the armature current, without the effect of the rotor magnets, are obtained using FEA for the fundamental and third sequences as shown

FEA Simulation Results
To obtain a fair comparison, the three machines are simulated at the same rotor speed of 240 rpm.For maximum torque per ampere operation, the direct current component is set to zero, while the quadrature component is 30 A. A transient analysis of the three machines was carried out.The resulting torque, core loss, and emf voltage are plotted in Figure 4. Finite element analysis was also used to estimate the machine efficiency by calculating the input power from the phase voltages and phase currents.Table 2 summarizes the main simulation results.It is clear that the machine torque increases as the number of rotor poles increases for the same rms stator current, which, in turn, improves the machine torque density.Moreover, the machine efficiency increases with the increase in rotor poles.However, the required phase voltage, and hence, the corresponding DC link voltage, increases with the number of poles.This is a critical design factor in electric vehicle applications.Although in [3] it was concluded that the 14 pole design has lower core loss, better flux distribution, and eliminates low order armature MMF harmonics, it does not consider the superior advantages of the other combinations, namely, the 20 slot/18 pole and 20 slot/22 pole, which offer higher efficiency and torque density for the same machine volume.

Mathematical Model of PM Modular Machines
A conventional phase model is developed here to represent the five-phase PM modular machine [12] considering fault conditions.The main advantage of mathematical model over the FE model is that the computational time is largely optimized.Moreover, this model can be embedded into a larger control system which will be very difficult with FE simulation.In the phase model, the voltage equations of the machine windings can be written in vector-matrix form as follows: the flux linking stator phases due to the rotor permanent magnets (PM).The machine torque is calculated using A coil carrying a current of 1 A represents the rotor permanent magnets.The equivalent number of turns is calculated such that the same flux is produced.Inductances are calculated using the winding function approach [15].The calculated inductances are the stator self-inductances, and the mutual inductances between the stator windings and the virtual field winding representing the PM effect.

Winding Function Method.
According to winding function theory [15], the mutual inductance between any two arbitrary windings i and j in any electric machine can be computed using where, μ o is permeability of free space, r is the rotor radius, l is the rotor length, θ m is the rotor angular position with respect to the stator reference phase a axis, φ is the angular position along the stator inner surface, g −1 (θ m , φ) is the inverse gap function, N i (θ m , φ) is the winding function of winding i, and N j (θ m , φ) is the winding function of winding j.
The inverse gap function is the reciprocal of the air gap length assuming radial flux and neglecting the effect of iron.Since in modular machines the air gap is approximately uniform, the inverse air gap function is assumed constant.

Winding Function of the Virtual Field
Winding.The winding function of the virtual field winding is shown in Figure 6.The effective number of turns is found by assuming that this virtual field winding carries 1 A, and equating the flux produced by the PM with the coil flux calculated from R pole is the reluctance of magnetic circuit seen by each pole calculated from (7), and A pole is the pole area as follows:

Inductance Calculation.
The inductance of different windings can be calculated by substituting by the winding function of different windings into (5).The stator inductances can be evaluated as follows: where, β is the coil span (for this example it equals π/10).Since the mutual inductance between phases is zero, the direct axis inductance is equal to the phase inductance given by (9).
The mutual inductance between the stator winding and the virtual field winding will vary sinusoidally with rotor angle θ if the space harmonics are neglected.However, the actual variation has an approximate trapezoidal wave shape.A plot of the corresponding flux linkage for a virtual field current of 1 A is shown in Figure 7(a), while its derivative with respect to position is shown in Figure 7(b).
The machine's time domain simulation model is shown in Figures 8(a) and 8(b), while the controller's block diagram is shown in Figure 8(c).As previously mentioned, a fivephase winding has 4 different possible sequences, two of them are forward, and two of them are backward sequences.Based on the MMF distributions, the first sequence is applied to the 20 slot/14 pole machine, while the second sequence is applied to the 20 slot/18 pole machine.The third sequence is applied to the 20 slot/22 pole machine, and the forth sequence can be used to drive the 20 slot/14 pole machine backward.Also, swapping the 2nd and 3rd sequences for 20 slot/18 pole and 20 slot/22 pole machine cases can be used to reverse the machine direction.The direct axis current components for both fundamental and third planes are set to zero.The quadrature current component for the fundamental is derived from the speed error using a PI controller, while that of the third is set to zero.The machine is simulated at a constant speed of 240 RPM.The transformation matrix that converts the dq command values to the corresponding phase values is given by where α = 2π/5, and s 1 and s 2 are the sequence order.
The constants s 1 and s 2 depend on the number of rotor poles.This is because each value of rotor pole corresponds to a specific stator MMF harmonic with a given phase sequence.The values of s 1 and s 2 are given in Table 3. Harmonic injection is possible with modular machine designs if the magnets give a quasi-rectangular spatial flux distribution [4].In this design, FE results show that the third harmonic component in the induced voltage is small.This makes the torque enhancement by stator current higher-order harmonic injection of negligible effect.Hence, the reference values for the third harmonic sequence are set to zero, as shown in the controller block diagram shown in Figure 8(c).The controller is based on conventional vector control of PM machines.Based on the machine speed error, a PI controller is used to decide the required quadrature component for the fundamental sequence plane.Unless field weakening is employed, the direct axis component is set to zero.The machine phase currents are measured and transformed to their sequence values.Four PI controllers are then used to obtain the required sequence voltage components, from which different machine phase voltages are calculated.
Figure 9(a) illustrates the machine speed, torque, phase currents, and phase voltages, for different numbers of rotor poles.In this case, the machine is accelerated at no load in 0.05 s to a speed of 240 rpm, then a 300 Nm load is applied at 0.1 s.It is evident that the 20 slot/18 pole machine has the lowest torque ripple.Moreover, the rms phase current decreases as the number of rotor poles increases.Finally, the rms stator phase voltage increases as the number of rotor poles increases.A comparison between phase currents and phase voltages for different numbers of rotor poles are shown in Figures 9(b) and 9(c), respectively.The machine performance indices are summarized in Table 4.It is evident from Table 4 that the machine torque per RMS current and efficiency improve with the increase in number of poles.However, the required DC link voltage increases with the increase in rotor poles, due to the increase in required phase voltage.
As general conclusions derived from simulation results, while the 20 slot/22 pole selection gives better efficacy, Lookup table

Rotation transformation
Rotation transformation  the 20 slot/18 pole selection corresponds to minimum torque ripple and the 20 slot/14 pole selection requires less DC-link voltage for the same machine speed.Hence, the optimum selection is mainly application dependent.

Fault Tolerant Applications
Multiphase machines have promising fault tolerant capabilities because they present additional degrees of freedom and can steadily operate with one or two phases open [11].When one or two phases are open, due to a device failure or a fault in the phase windings, a forward rotating field can be still obtained by modifying the currents in the healthy phases to ensure some constraints, namely, the same output torque and minimum copper loss [12,14].In this paper, the condition of one phase open is considered.
6.1.Optimum Phase Currents for One Phase Open.When one phase is open due to a device failure or a fault in the phase windings, a forward rotating field can be still obtained by setting the currents in the faulted phase to zero, and keeping the MMF and torque unchanged.During the fault conditions, a new set of currents for the healthy phases is applied [16][17][18].The current control strategies are chosen so as to obtain no zero-sequence current ( i = 0), and a reasonable average torque.If phase a is open, the current i a is zero and the currents in the remaining phases should satisfy the condition in (12), also illustrated in Figure 10   i To maintain an undisturbed rotating MMF in a fivephase machine with one phase open, the fundamental current amplitude of the healthy phases should increase by 1.38 times the value when all five phases are functional [18].
Consider the i d and i q axes currents corresponding to the torque and flux producing components.The current commands for the remaining four phases are the following: The phase currents' commands are labelled with subscripts x, y, z, and w.The relation between the actual reference current and these four currents will depend on the number of rotor poles, because the current sequence should be adapted to the corresponding number of rotor poles, as given in Table 5.  Phase currents and command currents relationship Sequence Fundamental Second sequence Third sequence 6.2.Simulation Results.In this section, the simulation results for the modular machine with one phase open and different numbers of rotor poles are carried out using both FEA and mathematical model in Matlab/Simulink.Figure 11 shows the simulation results for the three selections.Comparing the simulation results for the three motors, similar conclusions, as in the healthy case, will be obtained.It is evident that the 18 pole selection gives lower torque ripple than the other selections in the case of one phase open and it will be more suitable for applications that require smooth developed torque.On the other hand, as the number of rotor poles increases, the magnitude of the corresponding phase voltages increases.This makes the 14 pole selection a better selection to limit the required DC-link voltage.Finally, the 22 pole will correspond to higher average torque and, consequently, better machine torque density but with a relatively higher torque ripples.

Conclusion
This paper presents a comparison study between three different slot/pole combinations for a five phase external rotor modular PM machine.The study shows that machine torque density and efficiency are enhanced as the number of rotor poles increases.However, this is at the expense of an increase in the DC-link voltage.The 20 slot/22 pole combination is found to be the best selection.However, better flux distribution and less core loss are reported for the 20 slot/14 pole combination.From a fault tolerant perspective, the performance of the three slot/pole selections is similar; however, the 20 slot/18 pole combination gives minimum torque ripple for both healthy and faulted cases.Analysis of the healthy and fault case is provided using a phase model approach developed to simulate the machine dynamics.

Figure 1 :
Figure 1: (a) Stator winding layout and (b) three different designs for modular PM machines with external rotors.

Figure 2 :
Figure 2: (a) Phase a turn function and (b) Harmonic Spectrum.
) and 3(f), respectively.The machine flux distribution and the corresponding flux density distribution for both sequences are also shown in Figures3(g), 3(h), 3(i), and 3(j).

5. 1 . 1 .
Winding Function of the Stator Winding.The winding function of a single-layer stator winding with fractional SPP can be derived by plotting the flux distribution produced by phase a, calculated using FEA, assuming other phases carry zero current.The flux distribution is shown in Figure 5(a), and the corresponding air gap flux density is shown in Figure 5(b).

Figure 5 :Figure 6 :
Figure 5: Flux produced by one stator phase using FEA, (a) flux distribution and (b) air gap flux density.

Figure 7 :
Figure 7: (a) Phase a flux linkage due to PM and (b) flux linkage derivative with respect to rotor angle.

Table 1 :
Relation between different harmonics and the corresponding sequence for a five-phase system.

Table 2 :
Comparison between the three machines based on FE results.

Table 4 :
Machine performance indices for different number of rotor poles.

Table 5 :
Relation between phase currents and command currents.