Nonlinear Varying-Network Magnetic Circuit Analysis of Consequent-Pole Permanent-Magnet Motor for Electric Vehicles

: To conserve rare earth resources, consequent-pole permanent-magnet (CPPM) machine has been studied, which employs iron-pole to replace half PM poles. Meanwhile, to increase ﬂux-weakening ability, hybrid excitation CPPM machine with three-dimensional (3-D) ﬂux ﬂow has been proposed. Considering ﬁnite element method (FEM) is time-consuming, for the analysis of the CPPM machine, this paper presents a nonlinear varying-network magnetic circuit (NVNMC), which can analytically calculate the corresponding electromagnetic performances. The key is to separate the model of CPPM machine into different elements reasonably; thus, the reluctances and magnetomotive force (MMF) sources in each element can be deduced. While taking into account magnetic saturation in the iron region, the proposed NVNMC method can accurately predict the 3-D magnetic ﬁeld distribution, hence determining the corresponding back-electromotive force and electromagnetic power. Apart from providing fast calculation, this analytical method can provide physical insight on how to optimize the design parameters of this CPPM machine. Finally, the accuracy of the proposed model is veriﬁed by comparing the analytical results with the results obtained by using FEM. As a result, with so many desired attributes, this method can be employed for machine initial optimization to achieve higher power density. ma- chines is presented, and magnetic saturation and flux leakage is considered. This paper considers, for the first time, the use of 2-D NVNMC model to predict the performance of 3-D machine. In order to verify the accuracy of the proposed model, the NVNMC results and FEM results are quantitatively compared. Consequently, the proposed NVNMC method can be employed as a tool to optimize the CPPM machine in future study. NVNMC CPPM FEM


Introduction
Hybrid excitation synchronous machines (HESM) incorporate both permanent magnets and field winding for field excitation [1].By employing these two excitation field sources, the flux weakening capability and hence the speed range can be significantly improved [2].According to the arrangement of PM and excitation coils, HESM can be represented as series hybrid excitation (SHE) and parallel hybrid excitation (PHE) machines [3].
For PHE machines, the excitation fluxes produced by PMs and field winding have different trajectories [4].Therefore, the risk of irreversible demagnetization of the PMs can be avoided.Among those PHE machines, the consequent-pole PM (CPPM) machine possesses inherent field weakening capability [5].The PM and field-winding are housed in the rotor and stator separately.Therefore, air-gap flux can be regulated, and the demagnetization risk can be avoided.Meanwhile, as the field control winding is housed in the stator, the brushes or slip rings can be eliminated.Besides, compared with traditional PM machine, less rare-earth PM is required, as CPPM machine employed iron-pole to replace half of PM poles [6], as shown in Figure 1.
To improve machine's power density and torque density, various design methods have been proposed in CPPM machines [7,8].For example, in [9], the fundamental of the air-gap flux density distribution can be improved by optimizing the width of PMs.To improve machine's power density and torque density, various design methods have been proposed in CPPM machines [7,8].For example, in [9], the fundamental of the air-gap flux density distribution can be improved by optimizing the width of PMs.
In these papers, the finite element method (FEM) is generally employed.Thanks to the enormous improvement in digital computing, FEM has developed into a mature modeling technique, capable of accurately predicting motor performance.Meanwhile, many available FEM software packages, such as JMAG [4], Ansys [10] and MotorSolve [11], have largely simplified the implementation.However, compared to the analytical models, FEM is very time-consuming, especially for three-dimensional (3-D) motors [12].Therefore, analytical methods, such as nonlinear varying-network magnetic circuit (NVNMC) method and Fourier modeling method, are proposed for the initial design.Based on analytical equations and often tuned empirically, analytical methods can provide comparatively accurate prediction of the motor performance, and take the definite advantages of fast computation and simple implementation [13].Based on previous research results, a comparison between FEM, Fourier modeling method and NVNMC is listed in Table 1.It can be seen that NVNMC model is more suitable for saturation condition than Fourier modeling method [13]; thus, the former one is employed in this paper.However, the available NVNMC model is ill-suited for advanced machines with 3-D flux flow.
The purpose of this paper is to develop the NVNMC model of CPPM machines, hence assessing their electromagnetic performances, including the magnetic field distribution and back-electromotive force (EMF).The two-dimensional (2-D) NVNMC for CPPM machines is presented, and magnetic saturation and flux leakage is considered.This paper considers, for the first time, the use of 2-D NVNMC model to predict the performance of 3-D machine.In order to verify the accuracy of the proposed model, the NVNMC results and FEM results are quantitatively compared.Consequently, the proposed NVNMC method can be employed as a tool to optimize the CPPM machine in future study.
In Section 2, the NVNMC model of CPPM machine in [14] will be described.Then, the model is analyzed in Section 3. Section 4 will be devoted to solving this model.Then, the results will be presented in Section 5, and these are compared with FEM results.Finally, a conclusion will be drawn in Section 6.In these papers, the finite element method (FEM) is generally employed.Thanks to the enormous improvement in digital computing, FEM has developed into a mature modeling technique, capable of accurately predicting motor performance.Meanwhile, many available FEM software packages, such as JMAG [4], Ansys [10] and MotorSolve [11], have largely simplified the implementation.However, compared to the analytical models, FEM is very time-consuming, especially for three-dimensional (3-D) motors [12].Therefore, analytical methods, such as nonlinear varying-network magnetic circuit (NVNMC) method and Fourier modeling method, are proposed for the initial design.Based on analytical equations and often tuned empirically, analytical methods can provide comparatively accurate prediction of the motor performance, and take the definite advantages of fast computation and simple implementation [13].Based on previous research results, a comparison between FEM, Fourier modeling method and NVNMC is listed in Table 1.It can be seen that NVNMC model is more suitable for saturation condition than Fourier modeling method [13]; thus, the former one is employed in this paper.However, the available NVNMC model is ill-suited for advanced machines with 3-D flux flow.
The purpose of this paper is to develop the NVNMC model of CPPM machines, hence assessing their electromagnetic performances, including the magnetic field distribution and back-electromotive force (EMF).The two-dimensional (2-D) NVNMC for CPPM machines is presented, and magnetic saturation and flux leakage is considered.This paper considers, for the first time, the use of 2-D NVNMC model to predict the performance of 3-D machine.In order to verify the accuracy of the proposed model, the NVNMC results and FEM results are quantitatively compared.Consequently, the proposed NVNMC method can be employed as a tool to optimize the CPPM machine in future study.
In Section 2, the NVNMC model of CPPM machine in [14] will be described.Then, the model is analyzed in Section 3. Section 4 will be devoted to solving this model.Then, the results will be presented in Section 5, and these are compared with FEM results.Finally, a conclusion will be drawn in Section 6.

NVNMC Model
Figure 1a shows the structure of a three-phase 12-slot 4-pole CPPM machine.It should be mentioned that the rotor has radially magnetized bipolar PMs distributed alternatively in the circumferential direction, and other parts are iron poles (also called rotor tooth in this paper) [6].Due to symmetry, this CPPM machine can be considered as two 6/2-pole machines.Hence, the modeling is needed only for the 6/2-pole machine.
Although the machine involves 3-D flux flow, the corresponding 2-D NVNMC model is depicted in Figure 2. The key is to translate the 3-D flux path into a 2-D NVNMC model properly, namely the upper half of rotor is drawn inside the rotor yoke, while the upper half of stator is drawn outside the stator yoke.

NVNMC Model
Figure 1a shows the structure of a three-phase 12-slot 4-pole CPPM machine.It should be mentioned that the rotor has radially magnetized bipolar PMs distributed alternatively in the circumferential direction, and other parts are iron poles (also called rotor tooth in this paper) [6].Due to symmetry, this CPPM machine can be considered as two 6/2-pole machines.Hence, the modeling is needed only for the 6/2-pole machine.
Although the machine involves 3-D flux flow, the corresponding 2-D NVNMC model is depicted in Figure 2. The key is to translate the 3-D flux path into a 2-D NVNMC model properly, namely the upper half of rotor is drawn inside the rotor yoke, while the upper half of stator is drawn outside the stator yoke.In this model, PTR, PYR, PCR, PTS1, PTS, PYS, and PCS are the permeances of rotor tooth, rotor yoke, rotor consequent-yoke, stator inner tooth, stator outer tooth, stator yoke, and stator consequent-yoke, respectively, which vary with the nonlinear saturation in the corresponding magnetic paths, while PPM, PLR, and PLS are the permeances of PM, rotor toothto-PM leakage flux, and stator tooth-to-tooth leakage flux, respectively, which are of constant permeability [18].In this model, P TR , P YR , P CR , P TS1 , P TS , P YS , and P CS are the permeances of rotor tooth, rotor yoke, rotor consequent-yoke, stator inner tooth, stator outer tooth, stator yoke, and stator consequent-yoke, respectively, which vary with the nonlinear saturation in the corresponding magnetic paths, while P PM , P LR , and P LS are the permeances of PM, rotor tooth-to-PM leakage flux, and stator tooth-to-tooth leakage flux, respectively, which are of constant permeability [18].
It should be emphasized that although there is only one P A module drawn in front of each stator tooth in the model, it can be considered as each stator tooth having two P A modules connecting the PM and rotor tooth.So, there is a total of 24 P A modules in this model, and apparently, the value of some P A modules could become zero when the rotor is rotating.Besides, the connection lines between the outer P A and inner PM, as well as rotor teeth, are omitted in this figure.
In addition, Φ PM , Φ A , Φ B , Φ C and Φ DC represent the magnetic flux source (MFS) of the PM, phase A winding, phase B winding, phase C winding and field winding, respectively.It should be noted that the armature MFS supplied by phase A, phase B, and phase C winding is represented by eight Φ A , Φ B , and Φ C modules separately, while the MFS provided by DC winding is represented by twelve Φ DC modules.

Model Analysis
In NVNMC model, the magnetic field distribution is highly related to the geometry of the PM and the reluctance of its path.In this case, the magnetic flux source, the air-gap, and magnet reluctances are the most important values.
The basic equation which governs the permeance of each element is given by [21] where µ 0 and µ iron are the permeability of air and iron region, while b, h, and l are the width, height, and length of the element.
Although most elements' permeances can be easily acquired from (1), the calculations of main air-gap, stator and rotor cores are difficult and subtle.Hence, the calculation procedure for these parts is detailed in the following.

The Permeance of Air-Gap
As the permeance of air-gap is varying with the rotor position, the whole mechanical period is separated into seven regions, which can be defined as Region 1 0 ≤ θ ≤ θ 1 , Region 2 θ 1 < θ ≤ θ 2 , . . ., and Region 7 θ 6 < θ, where θ is the angle between the central lines of rotor and stator teeth.As depicted in Figure 3, θ 1 to θ 6 can be calculated as , where τ r and τ s are the pole pitches of rotor and stator, and β r and β s are the tooth arc of rotor and stator.
In addition, ΦPM, ΦA, ΦB, ΦC and ΦDC represent the magnetic flux source (MFS) of the PM, phase A winding, phase B winding, phase C winding and field winding, respectively.It should be noted that the armature MFS supplied by phase A, phase B, and phase C winding is represented by eight ΦA, ΦB, and ΦC modules separately, while the MFS provided by DC winding is represented by twelve ΦDC modules.

Model Analysis
In NVNMC model, the magnetic field distribution is highly related to the geometry of the PM and the reluctance of its path.In this case, the magnetic flux source, the air-gap, and magnet reluctances are the most important values.
The basic equation which governs the permeance of each element is given by [21] 0 (air) = (silicon steel) where μ0 and μiron are the permeability of air and iron region, while b, h, and l are the width, height, and length of the element.
Although most elements' permeances can be easily acquired from (1), the calculations of main air-gap, stator and rotor cores are difficult and subtle.Hence, the calculation procedure for these parts is detailed in the following.

The Permeance of Air-Gap
As the permeance of air-gap is varying with the rotor position, the whole mechanical period is separated into seven regions, which can be defined as Region 1 0 ≤ θ ≤ θ1, Region 2 θ1 < θ ≤ θ2, …, and Region 7 θ6 < θ, where θ is the angle between the central lines of rotor and stator teeth.As depicted in Figure 3, θ1 to θ6 can be calculated as θ1 = 1/2 × (βr − τs), θ2 = 1/2 × (βr − βs), θ3 = 1/2 × (τr − βs), θ4 = 1/2 × (βr + βs), θ5 = 1/2 × (βr + τs), and θ6 = 1/2 × (τr + τs), where τr and τs are the pole pitches of rotor and stator, and βr and βs are the tooth arc of rotor and stator.By employing the calculation method in [18], the air-gap permeance in each region can be calculated.It should be noted that the permanence between rotor yoke and stator pole should not be considered when calculating the PA between PM and stator pole.In this paper, PA-θ curve, the relationship between air-gap permeance and rotor angle, is shown in Figure 4.By employing the calculation method in [18], the air-gap permeance in each region can be calculated.It should be noted that the permanence between rotor yoke and stator pole should not be considered when calculating the P A between PM and stator pole.In this paper, P A -θ curve, the relationship between air-gap permeance and rotor angle, is shown in Figure 4.By employing the calculation method in [18], the air-gap permeance in each region can be calculated.It should be noted that the permanence between rotor yoke and stator pole should not be considered when calculating the PA between PM and stator pole.In this paper, PA-θ curve, the relationship between air-gap permeance and rotor angle, is shown in Figure 4.

The Permeance of Tooth, Yoke, and Consequent-Yoke
Figure 5a,b show the field distribution of the stator and rotor core separately.In order to analyze these two iron regions, the stator is separated as stator tooth (including PTS1 and PTS), stator yoke, and stator consequent-yoke, so that the corresponding magnetic field distributions can be seen as radial, circumferential, and axial directions, respectively [22].Similarly, the magnetic field distributions in the rotor tooth, rotor yoke, and rotor consequent-yoke are separated in accordance with the relevant parts of the stator.
To separate the stator and rotor core preciously, the widths of PCR and PCS, as well as the lengths of PYR and PYS, are specified based on FEM [18], and change dynamically with field current.But it should be mentioned that the FEM simulation only needs to run a few angles to get the widths of PCR and PCS, thus the time can be saved.Then, as the magnetic field distribution in each element is simplified as one-directional [23], (1) can be employed to calculate the PTR, PYR, PCR, PTS1, PTS, PYS, and PCS.In order to analyze these two iron regions, the stator is separated as stator tooth (including P TS1 and P TS ), stator yoke, and stator consequent-yoke, so that the corresponding magnetic field distributions can be seen as radial, circumferential, and axial directions, respectively [22].Similarly, the magnetic field distributions in the rotor tooth, rotor yoke, and rotor consequent-yoke are separated in accordance with the relevant parts of the stator.

The Magnetic Flux Source
The magnetomotive force (MMF) obtained from PM can be simplified as [24] where hm and Br is the thickness and the remanence of PM, separately.Therefore, the MFS of PM, namely ΦPM, can be acquired by where PPM is the permeance of PM, which can be obtained through (1).
The armature windings are formed by symmetrical three-phase windings.Ta phase C winding as an example, the MMF can be acquired through [25] = where NC and IC are the number of turns and the current of phase C winding, respecti It should be noticed that as the armature field flows as Figure 6 shows, the w MMF is evenly distributed to every stator tooth near the armature windings.To separate the stator and rotor core preciously, the widths of P CR and P CS , as well as the lengths of P YR and P YS , are specified based on FEM [18], and change dynamically with field current.But it should be mentioned that the FEM simulation only needs to run a few angles to get the widths of P CR and P CS , thus the time can be saved.Then, as the magnetic field distribution in each element is simplified as one-directional [23], (1) can be employed to calculate the P TR , P YR , P CR , P TS1 , P TS , P YS , and P CS .

The Magnetic Flux Source
The magnetomotive force (MMF) obtained from PM can be simplified as [24] F PM = B r h m /µ 0 (2) where h m and B r is the thickness and the remanence of PM, separately.Therefore, the MFS of PM, namely Φ PM , can be acquired by where P PM is the permeance of PM, which can be obtained through (1).
The armature windings are formed by symmetrical three-phase windings.Taking phase C winding as an example, the MMF can be acquired through [25] where N C and I C are the number of turns and the current of phase C winding, respectively.It should be noticed that as the armature field flows as Figure 6 shows, the whole MMF is evenly distributed to every stator tooth near the armature windings.For the field-winding, it is formed by DC coils.The MMF of this winding can be acquired through [25] = where NDC and IDC are the numbers of turns and the current of DC windings, separately.
It should be emphasized that the flux excited by DC winding flows from one stator tooth to another axially adjacent tooth through the stator consequent-yoke.Therefore, the whole MMF is evenly distributed to every stator tooth.
The basic equation which governs the MFS of each stator tooth is given by [25] =

Establishment of the Magnetic Circuit Equations
The nodal analysis and Kirchhoff's Law are employed to establish the magnetic circuit, the basic equation which governs the whole NVNMC model is given by [24] ( ) where n number of nodes; P(i, j) for i, j = 1, 2, …, n − 1, is the permeance of branch, which connects node i and j; F(i) for i = 1, 2, …, n − 1, is the node magnetic potential; Φs(i) for i = 1, 2, …, n − 1, is the node magnetic flux source.
It should be mentioned that the number of independent nodes in this model is n − 1, and the magnetic potential for node n is zero.

Calculation of the Magnetic Circuit Equations
As the magnetic circuit equation is established in the former section, the magnetic potential of all nodes can be obtained by [26] 1 -S F=P Φ (9) where P, F, and Φs are the node permeance matrix, node magnetic vector potential, and node magnetic flux source, respectively.Then the magnetic flux in each branch can be deduced by [18] ( ) ( ) ( ) ( ) For the field-winding, it is formed by DC coils.The MMF of this winding can be acquired through [25] F DC = N DC I DC (6) where N DC and I DC are the numbers of turns and the current of DC windings, separately.
It should be emphasized that the flux excited by DC winding flows from one stator tooth to another axially adjacent tooth through the stator consequent-yoke.Therefore, the whole MMF is evenly distributed to every stator tooth.
The basic equation which governs the MFS of each stator tooth is given by [25] Φ DC = F DC P TS /2 (7)
It should be mentioned that the number of independent nodes in this model is n − 1, and the magnetic potential for node n is zero.

Calculation of the Magnetic Circuit Equations
As the magnetic circuit equation is established in the former section, the magnetic potential of all nodes can be obtained by [26] F= P −1 •Φ S (9) where P, F, and Φs are the node permeance matrix, node magnetic vector potential, and node magnetic flux source, respectively.Then the magnetic flux in each branch can be deduced by [18] Φ(i, j) = [F(i)−F(j)]/P(i, j) where Φ(i, j) is the magnetic flux flows branch, which connects nodes i and j.
However, it should be noticed that the permeance of each branch is varying; i.e., the permeance of the main air-gap varies as the rotor rotates, and the permeance of stator and rotor element varies as the flux density changes.Different from air-gap permeance, which can be acquired by employing an appropriate P A − θ curve (depicted in Figure 4), the precise permeance of the iron region can only be obtained by solving the NVNMC model iteratively [17].In detail, by employing a reliable initial value of iron permeability, such as 3500 µ 0 , the magnetic flux in each iron element can be acquired through (10), as well as the magnetic flux density can be obtained.Therefore, the permeability of this region can be updated through the material's B-H curve.
Nevertheless, directly employing the new permeability in ( 8) is ill-suitable, because the iteration procedure could be an endless loop in this case.For example, if the flux density acquired in (k − 1)th step is higher than 2T, then the permeability µ k−1 iron would be relatively small, such as 100 µ 0 , then with this value, the flux density acquired in kth step could be much lower, such as 1T, and the corresponding permeability, namely µ k iron , could be 4000 µ 0 , so the flux density acquired in (k + 1)th step could be 2T again.Hence, to avoid this problem, the permeability for kth iteration can be obtained by [24] where p is a damping factor that is between 0 to 1, and is varying with the iteration steps.For example, damping factor could be 0.3, 0.1 and 0.02 in the first, second and last ten iterations.Compared with fixed damping factor, the precise result can be acquired quickly with varying damping factor.Taking the stator tooth as an example, the permeance of stator tooth is shown in Figure 7.It can be observed that this model cannot acquire precise result when damping factor is 0.25, and it takes a long time to solve the NVNMC model when damping factor is 0.02.However, the precise result can be acquired quickly with varying damping factor.where Φ(i, j) is the magnetic flux flows branch, which connects nodes i and j.However, it should be noticed that the permeance of each branch is varying; i.e., the permeance of the main air-gap varies as the rotor rotates, and the permeance of stator and rotor element varies as the flux density changes.Different from air-gap permeance, which can be acquired by employing an appropriate PA-θ curve (depicted in Figure 4), the precise permeance of the iron region can only be obtained by solving the NVNMC model iteratively [17].In detail, by employing a reliable initial value of iron permeability, such as 3500 μ0, the magnetic flux in each iron element can be acquired through (10), as well as the magnetic flux density can be obtained.Therefore, the permeability of this region can be updated through the material's B-H curve.
Nevertheless, directly employing the new permeability in ( 8) is ill-suitable, because the iteration procedure could be an endless loop in this case.For example, if the flux density acquired in (k − 1)th step is higher than 2T, then the permeability μ k−1 iron would be relatively small, such as 100 μ0, then with this value, the flux density acquired in kth step could be much lower, such as 1T, and the corresponding permeability, namely μ k iron, could be 4000 μ0, so the flux density acquired in (k + 1)th step could be 2T again.Hence, to avoid this problem, the permeability for kth iteration can be obtained by [24] ( ) where p is a damping factor that is between 0 to 1, and is varying with the iteration steps.For example, damping factor could be 0.3, 0.1 and 0.02 in the first, second and last ten iterations.Compared with fixed damping factor, the precise result can be acquired quickly with varying damping factor.Taking the stator tooth as an example, the permeance of stator tooth is shown in Figure 7.It can be observed that this model cannot acquire precise result when damping factor is 0.25, and it takes a long time to solve the NVNMC model when damping factor is 0.02.However, the precise result can be acquired quickly with varying damping factor.When the discrepancy between μ k−1 iron and μ k iron is tolerable, the iteration can be finished, and the precise solution of this NVNMC model can be acquired [23].

Results and Discussion
Based on the aforementioned solving procedure, major parameters of the CPPM machine are determined as listed in Table 2, and marked in Figure 8, which will be used for predicting electromagnetic performance.It should be noted that the method of determining winding turns is given in [11].Besides, the results acquired by NVNMC models and 3-D FEM are compared in this section.When the discrepancy between µ k−1 iron and µ k iron is tolerable, the iteration can be finished, and the precise solution of this NVNMC model can be acquired [23].

Results and Discussion
Based on the aforementioned solving procedure, major parameters of the CPPM machine are determined as listed in Table 2, and marked in Figure 8, which will be used for predicting electromagnetic performance.It should be noted that the method of determining winding turns is given in [11].Besides, the results acquired by NVNMC models and 3-D FEM are compared in this section.

Coil Flux Linkage
Figure 9 shows the waveforms of φA, namely the coil flux linkage of phase A. Meanwhile, the strengthening and weakening actions of the armature current are tested under different field winding currents with the value of 0 A, 10 A and −10 A. As aforementioned, hybrid excitation machines can regulate the field flux by varying the current of field winding, as shown in Figure 10.It can be seen that more flux flows at circumferential direction with 10 A field winding current, thus the width of PCR and PCS should be smaller; however, when the system is weakening the field flux, more flux flows at axial direction; thus, the width of PCR and PCS should be larger.But it should be mentioned that the width of PCR and PCS are not linearly related to the field current because of the saturation.Therefore, based on FEM, the widths of PCR and PCS are set as 40%, 25% and 55% of the stator and rotor pitches, when the field current is 0 A, 10 A and −10 A, respectively.As can be seen in Figure 9, a good agreement can be obtained between NVNMC and 3-D FEM.

Coil Flux Linkage
Figure 9 shows the waveforms of ϕ A , namely the coil flux linkage of phase A. Meanwhile, the strengthening and weakening actions of the armature current are tested under different field winding currents with the value of 0 A, 10 A and −10 A. As aforementioned, hybrid excitation machines can regulate the field flux by varying the current of field winding, as shown in Figure 10.It can be seen that more flux flows at circumferential direction with 10 A field winding current, thus the width of P CR and P CS should be smaller; however, when the system is weakening the field flux, more flux flows at axial direction; thus, the width of P CR and P CS should be larger.But it should be mentioned that the width of P CR and P CS are not linearly related to the field current because of the saturation.Therefore, based on FEM, the widths of P CR and P CS are set as 40%, 25% and 55% of the stator and rotor pitches, when the field current is 0 A, 10 A and −10 A, respectively.As can be seen in Figure 9, a good agreement can be obtained between NVNMC and 3-D FEM.

Back-EMF
Figure 11 illustrates EA, the back-EMF of phase A, versus rotor angle under differen IDC and phase current, which can be acquired through [14]    where ω is the angular velocity of the rotor.
As can be seen, there are few discrepancies between back-EMF calculated by th NVNMC model and FEM.Because the NVNMC method employs the simplified 2model to replace the 3-D machine, and the aim of this method is to rapidly predict th dynamic performance of CPPM machines with different dimensions and excitation con ditions.Therefore, the degradation of accuracy is tolerable in view of the time-saving.
However, it should be mentioned that, as shown in Figure 11a, the results of NVNM are not so accurate in some regions.The main reason is that the air-gap magnetic field very complicated, especially in region 5 and region 6 (θ5 < θ), thus when the rotor turn around 30 and 110 degrees, the results of these two methods are different.

Back-EMF
where ω is the angular velocity of the rotor.
As can be seen, there are few discrepancies between back-EMF calculated by the NVNMC model and FEM.Because the NVNMC method employs the simplified 2-D model to replace the 3-D machine, and the aim of this method is to rapidly predict the dynamic performance of CPPM machines with different dimensions and excitation conditions.Therefore, the degradation of accuracy is tolerable in view of the time-saving.

Power
As the winding internal resistance is very small, the copper loss is not consider this paper.Thus, the input power of phase A and total input power be simplified as However, it should be mentioned that, as shown in Figure 11a, the results of NVNMC are not so accurate in some regions.The main reason is that the air-gap magnetic field is very complicated, especially in region 5 and region 6 (θ 5 < θ), thus when the rotor turns around 30 and 110 degrees, the results of these two methods are different.

Power
As the winding internal resistance is very small, the copper loss is not considered in this paper.Thus, the input power of phase A and total input power be simplified as [27] P A = I A E A (13) and Figure 12 shows the variation of the input power with rotor position, fed by threephase AC current with the peak of 3 A. As will be seen, the results acquired from NVNMC model compare well with those from 3-D FEM.

Discussion
Consistent with the theoretical analysis, the proposed NVNMC method can p comparatively reliable results in a shorter time.As a matter of fact, the NVNMC so Matlab only needs to run 2 min, while the FEM software JMAG needs to run 27 h cycle with same computer.Besides, it is easy to see which part is too saturated and the flow of machine's magnetic field; thus, we can optimize the motor paramete physical insight.
However, the result of NVNMC method is still not as accurate as FEM metho main reason is that the NVNMC method separates the machine into dozens of ele and assuming the magnetic field is evenly distributed in each element to calcul permeability.However, uneven distribution of magnetic field is unavoidable, th inappropriate to calculate the permeance of each element with single permeability.ory, this problem can be solved by separating the machine into smaller elements, a method separates the machine into numerous meshes, but the running time is ine increased.Therefore, it is better to use NVNMC method for initial design and ref design with FEM method, as shown in Figure 13.Considering there are too many p eters that can have an impact on the performance of the motor, how to optimiz parameters quickly and reasonably will become the future research direction of this

Discussion
Consistent with the theoretical analysis, the proposed NVNMC method can provide comparatively reliable results in a shorter time.As a matter of fact, the NVNMC software Matlab only needs to run 2 min, while the FEM software JMAG needs to run 27 h for one cycle with same computer.Besides, it is easy to see which part is too saturated and limits the flow of machine's magnetic field; thus, we can optimize the motor parameters with physical insight.
However, the result of NVNMC method is still not as accurate as FEM method.The main reason is that the NVNMC method separates the machine into dozens of elements, and assuming the magnetic field is evenly distributed in each element to calculate the permeability.However, uneven distribution of magnetic field is unavoidable, thus it is inappropriate to calculate the permeance of each element with single permeability.In theory, this problem can be solved by separating the machine into smaller elements, as FEM method separates the machine into numerous meshes, but the running time is inevitably increased.Therefore, it is better to use NVNMC method for initial design and refine the design with FEM method, as shown in Figure 13.Considering there are too many parameters that can have an impact on the performance of the motor, how to optimize these parameters quickly and reasonably will become the future research direction of this study.

3. 2 .
Figure5a,b show the field distribution of the stator and rotor core separately.In order to analyze these two iron regions, the stator is separated as stator tooth (including P TS1 and P TS ), stator yoke, and stator consequent-yoke, so that the corresponding magnetic field distributions can be seen as radial, circumferential, and axial directions, respectively[22].Similarly, the magnetic field distributions in the rotor tooth, rotor yoke, and rotor consequent-yoke are separated in accordance with the relevant parts of the stator.

Figure 5 .
Figure 5. Separation of stator and rotor core: (a) One stator segment; (b) One rotor segment.

Figure 6 .
Figure 6.Armature field excited by phase C winding.

Figure 6 .
Figure 6.Armature field excited by phase C winding.

Figure 7 .
Figure 7. Permeance of rotor tooth with different damping factors.

Figure 7 .
Figure 7. Permeance of rotor tooth with different damping factors.

Figure 9 .
Figure 9. Flux linkage of phase A with different DC currents.Figure 9. Flux linkage of phase A with different DC currents.

Figure 9 .Figure 10 .
Figure 9. Flux linkage of phase A with different DC currents.Figure 9. Flux linkage of phase A with different DC currents.

Figure 11
Figure 11  illustrates E A , the back-EMF of phase A, versus rotor angle under different I DC and phase current, which can be acquired through[14]

Figure 11 .
Figure 11.Back-EMF of phase A with different DC currents: (a) No phase current; (b) Feeding 3 phase AC current with the peak of 3 A.

Figure 11 .
Figure 11.Back-EMF of phase A with different DC currents: (a) No phase current; (b) Feeding 3 phase AC current with the peak of 3 A.

Figure 12 .
Figure 12.Input power with different DC current: (a) Power of phase A; (b) Total input power.

Figure 12 .
Figure 12.Input power with different DC current: (a) Power of phase A; (b) Total input power.

Table 1 .
Comparison of different analysis methods.