Two-Dimensional Exact Subdomain Technique of Switched Reluctance Machines with Sinusoidal Current Excitation

This paper presents a two-dimensional (2D) exact subdomain technique in polar coordinates considering the iron relative permeability in 6/4 switched reluctance machines (SRM) supplied by sinusoidal waveform of current (aka, variable flux reluctance machines). In non-periodic regions (e.g., rotor and/or stator slots/teeth), magnetostatic Maxwell’s equations are solved considering non-homogeneous Neumann boundary conditions (BCs). The general solutions of magnetic vector potential in all subdomains are obtained by applying the interface conditions (ICs) in both directions (i.e., rand θ-edges ICs). The global saturation effect is taken into account, with a constant magnetic permeability corresponding to the linear zone of the nonlinear B(H) curve. In this investigation, the magnetic flux density distribution inside the electrical machine, the static/dynamic electromagnetic torques, the magnetic flux linkage, the self-/mutual inductances, the magnetic pressures, and the unbalanced magnetic forces (UMFs) have been calculated for 6/4 SRM with two various non-overlapping (or concentrated) windings. One of the case studies is a M1 with a non-overlapping all teeth wound winding (double-layer winding with left and right layer) and the other is a M2 with a non-overlapping alternate teeth wound winding (single-layer winding). It is important to note that the developed semi-analytical model based on the 2D exact subdomain technique is also valid for any number of slot/pole combinations and for non-overlapping teeth wound windings with a single/double layer. Finally, the semi-analytical results have been performed for different values of iron core relative permeability (viz., 100 and 800), and compared with those obtained by the 2D finite-element method (FEM). The comparisons with FEM show good results for the proposed approach.


Introduction
Benefiting from the advantages of a simple mechanical structure-the rotor does not carry any windings, commutators, or permanent magnets (PMs)-and a robust, fault-tolerant nature, low-cost maintenance, high-thermal capability, and high-speed potential [1][2][3], SRM is receiving renewed attention as a viable candidate for various adjustable-speed and high-torque applications such as in the automotive and traction fields [4][5][6][7][8].
However, a major disadvantage of this machine is the undesirable electromagnetic vibration and acoustic noise, which are mainly excited by the radial UMF acting on the salient stator and rotor poles [9,10].Moreover, it poses a drawback for SRM in noise-sensitive applications and still creates bottlenecks in vehicle propulsion.It is important to consider noise and vibration problems during the process of electrical machine design.Electrical machine noise and vibration are mainly of electromagnetic, aerodynamic, and mechanical origin, the most important of which are generated by electromagnetic sources [11].Also, in SRM, the attraction magnetic force can be divided into tangential and radial components relative to the rotor.The tangential magnetic force is converted into rotational torque, and the radial magnetic force converts into magnetic pressure equal to the radial magnetic force per unit area of the stator tooth and UMFs, which contributes to the radial vibration behaviour and therefore the motor noise [9,12].A perfect machine with balanced stator windings should have net zero UMFs on the stator structure.However, UMFs can be present in machines having diametrically asymmetric disposition of slots and phase windings [13,14].This magnetic force acts on the stator of these machine configurations due to an asymmetric magnetic field distribution in the air gap.
In the interest for design and optimization of electrical machines, there are various modelling methods; the first step in these is the magnetic field calculation.Some comprehensive reviews of the models of electrical machines for magnetic field prediction along with their (dis)advantages can be found in [15][16][17][18][19][20][21][22][23][24] and their references.Currently, the Maxwell-Fourier method is one of the most used semi-analytical methods, and combines the very accurate electromagnetic performances calculation with a reduced computation time compared to numerical methods.In models from this method (viz., multi-layer models, eigenvalues model, and subdomain technique), the magnetic field solutions are based on the formal resolution of Maxwell's equations by using the separation of variables method and the Fourier's series.In electromagnetic devices, the major assumption is that an infinite permeability of iron parts has to be assumed [25].Therefore, the global and/or local saturation effect is neglected.It is interesting to note that an overview of the existing (semi-)analytical models in the Maxwell-Fourier method with a global and/or local saturation effect has been realized in [24], where some details and the (dis)advantages of these techniques can be found.To overcome that issue, Spranger et al. (2016) [21] and Dubas et al. (2017) [24,26] have recently developed new techniques to account for finite soft-magnetic material permeabilities: • multi-layer models using the convolution theorem (i.e., Cauchy's product theorem).The adjacent regions (e.g., rotor and/or stator slots/teeth) are assumed to be one homogeneous region with a relative permeability developed as a Fourier's series expansion; • the subdomain technique using a superposition that allows for any non-periodic subdomain.The subdomain connection is performed directly in both directions.The general solutions of Maxwell's equations are deduced by applying the principle of superposition by respecting the BCs on the various edges of subdomains.
For the same reason, another technique based on subdomain technique and Taylor polynomial has been developed and only applied in spoke-type PM synchronous machines (PMSM) [27,28].Spranger's approach has been extended and used in different machines with only the global saturation effect.It has been applied with the finite soft-magnetic material permeability in synchronous reluctance machine [29], surface-mounted PMSM [30], and many structures of PMSMs (i.e., for inset-/surface-/spoke-type PMSMs with different PM magnetization patterns and internal/external rotor) [31], with the nonlinear B(H) curve in switched reluctance machine [32,33].The Dubas superposition technique has been implemented in radial-flux electrical machines with(out) PMs supplied by a direct or alternate current (with any waveforms) [34].This technique has been extended to: (i) the thermal modelling for the steady-state temperature distribution in rotating electrical machines [35], and (ii) elementary subdomains in the rotor and stator regions for full prediction of magnetic field in rotating electrical machines with the local saturation effect solving by the Newton-Raphson iterative algorithm [36].The Dubas superposition technique is very interesting since, like Spranger's approach, it enables the magnetic field calculation in iron parts of slotted structures.Apart from its complexity, the main downfall of Spranger's approach is that it suffers from the Gibb's phenomenon at boundaries between slots and teeth.This introduces inaccuracies in the computation of the field and results in higher computational times [37].
In this paper, the authors propose applying the Dubas superposition technique in polar coordinates [26] to SRM with sinusoidal current excitation, which has not yet been realized in the literature.The soft magnetic material permeability is constant corresponding to the linear zone of the B(H) curve.Nevertheless, as in [33,36], it should be mentioned that the material properties could be updated iteratively to take the nonlinear B(H) curve of the material into account.However, this is beyond the scope of the paper.In this investigation, the magnetic flux density distribution inside the machine, electromagnetic performances and non-intrinsic UMFs have been calculated for 6/4 SRM supplied by sinusoidal waveform of current with two various non-overlapping (or concentrated) windings.One of the case studies is a M1 with a non-overlapping all teeth wound winding (double-layer winding with left and right layer) and the other is a M2 with a non-overlapping alternate teeth wound winding (single-layer winding).All results obtained with the proposed semi-analytical model are verified by 2D FEM [38] for different values of iron core relative permeability (viz., 100 and 800).The comparisons with FEM show good results.

Machine Geometry and Assumptions
Figure 1 represents the studied SRMs having two various non-overlapping windings: (i) M1 with a non-overlapping all teeth wound winding (double-layer winding with left and right layer) (see Figure 1a), and (ii) M2 with a non-overlapping alternate teeth wound winding (single-layer winding) (see Figure 1b).The three-phase SRMs have six stator slots and four rotor slots, and do not contain any stator tooth tips.The main geometrical parameters of two studied SMRs are shown in Figure 1 and are given in Table 1 for the semi-analytical and numerical comparisons.These machines have been partitioned into nine regions as shown on Figure 2, viz.,

•
Region I is the air gap;  Region XI is the non-periodic air gap (i.e., between the two-layer winding of the stator slots).
The semi-analytical model, based on the exact subdomain technique, is formulated in 2D, in polar coordinates, and in magnetic vector potential with the following assumptions:

•
The end-effects are neglected, i.e., A = {0; 0; A z }; • The eddy-current effects in the materials are neglected;

•
The current density in the stator slots has only one component along the z-axis, i.e., J = {0; 0; J z }; • The magnetic materials are considered as isotropic with constant magnetic permeability corresponding to linear zone of the B(H) curve;

•
The stator and rotor slots/teeth have radial sides (see Figure 2).
However, it accounts for:

•
The internal/external rotor topology;

General Solution with Non-Homogeneous Neumann BCs
Magnetic vector potential A is calculated analytically with solving the magnetostatic Maxwell's equations with the separation of variables method, viz., where µ 0 is the vacuum permeability.
According to [24,26], the solutions to A in all regions of conventional SRM are: • Air gap subdomain (Region I): The solution of (1) in Region I, r ∈ [R 3 ; R 4 ] & ∀θ, is defined by: where n is a positive integer, and {A 10 ; A 20 ; A 1n ; A 4n } are the integration constants of Region I.

•
Stator and rotor yoke subdomain (Region II and III): ∀θ, can be written as: where {A 50 ; A 5n ; A 6n } are the integration constants of Region II.

• i-th Stator slot subdomain (Region VI and VII):
The solution of (2 , is defined by: where m and k are positive integers, γ 1i = γ i − (e + f )/2 and f are respectively the position and opening width of first layer winding in the i-th stator slot, {C 1i0 ; C 2i0 ; C 1im ; C 2im ; C 3ik ; C 4ik } are the integration constants of Region VI, v m f = mπ/ f and λ ks = kπ/ln(R 5 /R 4 ) are respectively the periodicity of A zV Ii in θ-and r-edges.

•
i-th Non-periodic air gap and i-th stator tooth subdomain (Region XI and VIII): , can be obtained directly from (5) with J1(i) z = 0.

•
j-th Rotor slot and j-th rotor tooth subdomain (Region IV and V): The solution of (1) in Region IV, r ∈ [R 2 ; R 3 ] & θ ∈ α j − a/2; α j + a/2 , is defined by: where α j and a are respectively the position and opening width of j-th rotor slot, B 1j0 ; B 2j0 ; B 1jm ; B 2jm ; B 3jk ; B 4jk are the integration constants of Region IV, v ma = mπ/a and λ kr = kπ/ln(R 3 /R 2 ) are respectively the periodicity of A zIV j in θ-and r-edges.

Magnetic Flux Density
The field vectors B = {B r ; B θ ; 0} and H = {H r ; H θ ; 0} are coupled by: where µ rc is the relative recoil permeability of iron parts.Using B = ∇ × A, the components of B = ∇ × A can be deduced by (9)

Stator Current Density Source
The stator current densities in the stator slots for double-layer concentrated winding are defined as [33]: where i g = i a i b i c is the vector of phase currents whose currents' waveform is sinusoidal with is the surface of the stator slot coil, and are the transpose of the connection matrix between the three phases and the stator slots that represent the distribution of stator windings in the slots of the M1 with all teeth wound (double-layer winding with left and right layer) (see Figure 1a) is given by [33]: For the M2 with alternate teeth wound (single-layer winding) (see Figure 1b), the same model is used with few modifications:

•
The opening of the non-periodic air gap will be equal to zero (i.e., e = 0);

•
The stator current density will be equal to These connection matrices can be generated automatically by using the ANFRACTUS TOOL developed in [39].

Boundary Conditions
The ICs in this semi-analytical model can be divided into two types, viz., • θ-edges ICs: over angle interval for given radius value {R 2 ; R 3 ; R 4 ; R 5 }; • r-edges ICs: over radius interval for given angle α j ± a/2; Therefore, we obtain on the: • θ-edges ICs: -The ICs between Region II, IV and V at r = R 2 as: -The ICs between Region I, IV and V at r = R 3 are similar to ( 13)-( 16) by replacing II with I and R 2 with R 3 .

-
The ICs between Region I, VI, VII, VIII and XI at r = R 4 as: -The ICs between Region III, VI, VII, VIII and XI at r = R 5 are similar to ( 17)-( 24) by replacing I with III and R 4 with R 5 .

• r-edges ICs:
-The ICs between Region IV and V at A zIV(j+1) r, α j+1 − a/2 = A zV j r, -The ICs between Region VII and VIII at γ i + c/2 = δ i − d/2 and between Region VI and VIII at -The ICs between Region VI and XI at γ i − e/2 = γ i − c/2 + f and between Region VII and XI at The system of the 36 BCs matrix (Equations ( 13)-( 36)) is used to determine the coefficients of A in nine regions.
Figure 3 briefly represents a flowchart of the subdomain technique.-The ICs between Region IV and V at -The ICs between Region VII and VIII at 2 2 and between Region VI and VIII at -The ICs between Region VI and XI at The system of the 36 BCs matrix (Equations ( 13)-( 36)) is used to determine the coefficients of A in nine regions.
Figure 3 briefly represents a flowchart of the subdomain technique.To solve the Cramer's system, the number of integration constants is equal to where N , M and K are the finite numbers of spatial harmonic terms in the various regions.

Torque, Flux Linkage and Inductance Calculations
The semi-analytical model, based on the 2D exact subdomain technique and taking into account the iron core relative permeability, is used to determine the static/dynamic electromagnetic torque,  To solve the Cramer's system, the number of integration constants is equal to 2 where N, M and K are the finite numbers of spatial harmonic terms in the various regions.

Torque, Flux Linkage and Inductance Calculations
The semi-analytical model, based on the 2D exact subdomain technique and taking into account the iron core relative permeability, is used to determine the static/dynamic electromagnetic torque, the magnetic flux linkage, and the self-/mutual inductances whose various formulas have been clarified in [32,33].

Magnetic Pressure and UMF Calculations
The magnetic pressure P is the magnetic force per unit area of the stator tooth.It can be calculated both by Maxwell's stress tensor and by finite element analysis.The r-and θ-components of P are calculated from the spatial magnetic field in the air gap middle at R g = (R 4 + R 3 )/2 [14]: where θ r = Ω • t + θ rs0 is the temporal rotor angle with Ω the mechanical pulse and θ rs0 the initial mechanical angular position between the rotor and the stator at the instant t = 0 s.The xand y-components of non-intrinsic UMF F are calculated at R g over θ = [0; 2π] as [14]: where L u is the axial length of the machine.The acoustic noise and vibration is primarily due to the rotor eccentric position with respect to the stator bore, the UMF, if present in a motor even with perfectly aligned shaft, can create the rotor eccentricity.Moreover, P and F are transmitted through the teeth from the air gap to the yoke, which may cause deformation on the stator rings resulting from the rotor displacement and result in excessive acoustic noise and vibration.Different vibration modes are commonly called "mode shapes" having their own natural mode frequency.Any particular mode shape is excited when its natural mode frequency matches with any of the harmonics of P and F [40].

Results and Validations
The developed model (see Section 2) considering finite soft-magnetic material permeability is used to determine the magnetic flux density distribution inside the electrical machines as well as the electromagnetic performances for 6/4 SRM with two various non-overlapping (or concentrated) windings.The main dimensions and parameters of studied machines are given in Table 1.The results of the semi-analytical model are verified by 2D FEM.

Magnetic Flux Density Distribution
The waveforms of rand θ-components of the magnetic flux density in the various regions are computed with a finite number of harmonic terms, viz., N = 200 and M = K = 30.The analytic calculation of magnetic flux density distribution in all regions is done considering the same relative permeability in all iron parts (i.e., stator/rotor yoke and teeth).The soft magnetic material permeability is constant corresponding to the linear zone of the B(H) curve.However, it is possible to use a different relative permeability value for each region [28,33,36].
In Figures 4 and 5, a comparison between the numerical results and semi-analytical predictions is shown the rand θ-components of B in the air gap middle (i.e., Region I at R g = (R 4 + R 3 )/2) for two studied SRMs (i.e., M1 and M2).The simulations are done for two different values of iron core relative permeability (viz., 100 and 800).It can be seen an asymmetric distribution in the rand θ-component of the air gap magnetic flux density in M2 in contrary to M1, this is due to the diametrically asymmetric disposition of slots and phase windings as shown in Figure 1b.The simulations were done for both values of iron core relative permeability.One can see that a very good agreement is obtained for the various components of B in all regions.The simulations were done for both values of iron core relative permeability.One can see that a very good agreement is obtained for the various components of B in all regions.

Static/Dynamic Electromagnetic Torques
For two studied SRMs (i.e., M1 and M2) and for both values of iron core relative permeability, Figures 16-19 show the waveform as well as the harmonic spectrum of the static/dynamic electromagnetic torques for full-load condition (viz., 15 A @ 1500 rpm).The static electromagnetic torque represents the torque due to a single phase of the electrical machine (e.g., due to phase A).The dynamic electromagnetic torque represents the torque when the three phases are powered or due to the combination of three static electromagnetic torques.The good agreement between the results from 2D FEM and the proposed semi-analytical model can be seen.It is interesting to note that the ripple torques are more important for M2 (see Figure 19) with the same operating point.The simulations were done for both values of iron core relative permeability.One can see that a very good agreement is obtained for the various components of B in all regions.

Static/Dynamic Electromagnetic Torques
For two studied SRMs (i.e., M1 and M2) and for both values of iron core relative permeability, Figures 16-19 show the waveform as well as the harmonic spectrum of the static/dynamic electromagnetic torques for full-load condition (viz., 15 A @ 1500 rpm).The static electromagnetic torque represents the torque due to a single phase of the electrical machine (e.g., due to phase A).The dynamic electromagnetic torque represents the torque when the three phases are powered or due to the combination of three static electromagnetic torques.The good agreement between the results from 2D FEM and the proposed semi-analytical model can be seen.It is interesting to note that the ripple torques are more important for M2 (see Figure 19) with the same operating point.

Static/Dynamic Electromagnetic Torques
For two studied SRMs (i.e., M1 and M2) and for both values of iron core relative permeability, Figures 16-19 show the waveform as well as the harmonic spectrum of the static/dynamic electromagnetic torques for full-load condition (viz., 15 A @ 1500 rpm).The static electromagnetic torque represents the torque due to a single phase of the electrical machine (e.g., due to phase A).The dynamic electromagnetic torque represents the torque when the three phases are powered or due to the combination of three static electromagnetic torques.The good agreement between the results from 2D FEM and the proposed semi-analytical model can be seen.It is interesting to note that the ripple torques are more important for M2 (see Figure 19) with the same operating point.

Magnetic Flux Linkage and Self-/Mutual Inductances
For full-load condition (viz., 15 A @ 1500 rpm), the induced magnetic flux linkage per phase of two studied SRMs (i.e., M1 and M2) are given in Figure 20.The simulations were done for both values of iron core relative permeability.
Figures 21 and 22 show the self-and mutual inductance for M1 and M2, respectively; the simulation is done for nominal current.One can see that the self-inductance is slightly more important, while the mutual inductance is a much more important and negative value for M2.The obtained results confirm the accuracy of the proposed semi-analytical model, considering both amplitude and waveform.
of iron core relative permeability.
Figures 21 and 22 show the self-and mutual inductance for M1 and M2, respectively; the simulation is done for nominal current.One can see that the self-inductance is slightly more important, while the mutual inductance is a much more important and negative value for M2.The obtained results confirm the accuracy of the proposed semi-analytical model, considering both amplitude and waveform. of iron core relative permeability.Figures 21 and 22 show the self-and mutual inductance for M1 and M2, respectively; the simulation is done for nominal current.One can see that the self-inductance is slightly more important, while the mutual inductance is a much more important and negative value for M2.The obtained results confirm the accuracy of the proposed semi-analytical model, considering both amplitude and waveform.

Magnetic Pressure and Non-Intrinsic UMFs
Figures 23 and 24 show the rand θ-components of P in function of space angle for the two studied SRMs (i.e., M1 and M2) and both values of iron core relative permeability.The analytical radial magnetic pressure in function of temporal rotor angle and the spatial angle is represented in Figure 25 for a no-load condition.Figure 26 shows the fast Fourier transform in 2D (FFT2D) of analytical radial magnetic pressure P r for M1 and M2.The x-and y-components of F in M1 and M2, for a no-load condition, are shown in Figure 27.It is clear that the non-intrinsic UMFs can be significant in SRMs, having diametrically asymmetric disposition of non-overlapping winding and due to the asymmetric magnetic field distribution in the air gap (see Figure 5).
The UMFs in M1 are null due to the proper choice of the armature winding type in the stator with same phase windings in diametrically opposite slots (see Figure 1a).Figure 28 shows the locus of the non-intrinsic UMF in M2.It is interesting to note that the UMFs and the magnetic pressures increases with the increase of the iron core relative permeability.One can see that the proposed semi-analytical model taking into account the iron core relative permeability gives good results compared to FEM.The obtained results confirm the previous interpretation of the UMF.The component of the mode shapes one correspond to the wave number (i.e., r = 1) of the UMF does not appear in M1 due to the symmetric distribution of the air gap magnetic field with respect to the space angle p, contrary to M2.The component corresponds to the wave number (i.e., r = 2) is appeared in both SRMs, because of the asymmetric distribution magnetic field in the air gap with respect to the space angle p/2, as seen in Figures 4 and 5.Moreover, it can be seen that the other modes that appeared in the two SRMs are multiples of the least non-null wave number r min = 1 in M2 and r min = 2 in M1.

Conclusions
In this paper, we have developed a 2D exact subdomain technique in polar coordinates considering the iron relative permeability for SRM supplied by a sinusoidal waveform of the current (aka, variable flux reluctance machines).This semi-analytical model, based on the scientific works of [24,33], predicts the magnetic flux density distribution inside the electrical machine as well as the electromagnetic performances.It has been applied to 6/4 SRM with two various non-overlapping (or concentrated) windings.These two configurations of non-overlapping winding have been considered to show their effect on the UMFs.Moreover, the spectrum of UMF permit us to study the effect of each harmonic on the vibrations of these machines.However, this research proved that the SRM with an asymmetric disposition of winding is more prone to higher levels of vibration than the SRM with a symmetric disposition due to the UMFs' presence, which is the main source of vibration and acoustic noise.All results confirmed the accuracy of the proposed model.It can be considered a reliable alternative to FEM for analysis of SRMs.
Mesh machine & Definition of all regionsDefinitions of the solutions in all regions of conventional SRM Performed of ICs Resolution of linear system (Cramer's system) obtained from the ICs between the various regions to find the integration constants of A in all the regions.Calculation of B in the various regions

Figures 6 and 7
Figures6 and 7show the magnitude of B in all machine's regions for two studied SRMs (i.e., M1 and M2), we can notice that M2 is more saturated than M1.In Figures 8-15, for two studied SRMs (i.e., M1 and M2), a numerical and semi-analytical comparison is shown the r-and θ-components of B in the • stator yoke middle (see Figure 8 for M1 and Figure 9 for M2); • rotor yoke middle (see Figure 10 for M1 and Figure 11 for M2); • stator slots/non-periodic air gap/teeth middle (see Figure 12 for M1 and Figure 13 for M2); • rotor slots/teeth middle (see Figure 14 for M1 and Figure 15 for M2).

Figure 4 .
Figure 4. Waveform of the magnetic flux density in the air gap middle (i.e., Region I) for M1: (a) rand (b) θ-component.

Figures 6 and 7
Figures6 and 7show the magnitude of B in all machine's regions for two studied SRMs (i.e., M1 and M2), we can notice that M2 is more saturated than M1.In Figures 8-15, for two studied SRMs (i.e., M1 and M2), a numerical and semi-analytical comparison is shown the r-and θ-components of B in the • stator yoke middle (see Figure 8 for M1 and Figure 9 for M2); • rotor yoke middle (see Figure 10 for M1 and Figure 11 for M2); • stator slots/non-periodic air gap/teeth middle (see Figure 12 for M1 and Figure 13 for M2); • rotor slots/teeth middle (see Figure 14 for M1 and Figure 15 for M2).

Figure 5 .
Figure 5. Waveform of the magnetic flux density in the air gap middle (i.e., Region I) for M2: (a) rand (b) θ-component.

Figures 6 Figure 6 .
Figures 6 and 7 show the magnitude of B in all machine's regions for two studied SRMs (i.e., M1 and M2), we can notice that M2 is more saturated than M1.Math.Comput.Appl.2018, 23, x FOR PEER REVIEW 11 of 22

Figure 7 .
Figure 7. Flux density inside the machine in M2: (a) analytic and (b) FEM.In Figures 8-15, for two studied SRMs (i.e., M1 and M2), a numerical and semi-analytical comparison is shown the rand θ-components of B in the • stator yoke middle (see Figure 8 for M1 and Figure 9 for M2); • rotor yoke middle (see Figure 10 for M1 and Figure 11 for M2); • stator slots/non-periodic air gap/teeth middle (see Figure 12 for M1 and Figure 13 for M2); • rotor slots/teeth middle (see Figure 14 for M1 and Figure 15 for M2).

Figure 9 .
Figure 9. Waveform of the magnetic flux density in the stator yoke middle for M2: (a) rand (b) θ-component.

Figure 15 .
Figure 15.Waveform of the magnetic flux density in the rotor slots/teeth for M2: (a) rand (b) θ-component.

Figure 15 .
Figure 15.Waveform of the magnetic flux density in the rotor slots/teeth for M2: (a) r-and (b) θ-component.

Figure 16 .
Figure 16.The static electromagnetic torque due to phase-A for M1: (a) waveform; (b) harmonic spectrum.Figure 16.The static electromagnetic torque due to phase-A for M1: (a) waveform; (b) harmonic spectrum.

Figure 22 .
Figure 22.Waveform of the (a) self-and (b) mutual inductance for M2.

Figure 28 .
Figure 28.Locus of the non-intrinsic UMF in M2 for no-load condition.Figure 28.Locus of the non-intrinsic UMF in M2 for no-load condition.

Figure 28 .
Figure 28.Locus of the non-intrinsic UMF in M2 for no-load condition.Figure 28.Locus of the non-intrinsic UMF in M2 for no-load condition.

Table 1 .
Parameters of the studied SRMs.

Table 1 .
Parameters of the studied SRMs.

Table 1 .
Parameters of the studied SRMs.