Analytical Calculation of Armature Reaction Field of the Interior Permanent Magnet Motor

The energy crisis and environmental concerns worldwide have helped usher in the age of electric vehicles (EVs) and hybrid EVs (HEVs). The interior permanent magnet motors (IPMMs) are widely used in these vehicles. The analysis of the armature reaction field is the most critical issue in the study of IPMMs since it determines the characters of torque, efficiency, vibration, and the radiated acoustic noise. This paper provides a calculation method of the armature reaction magnetic field (ARMF) of an IPMM. First, the formulas of ARMF without magnetic barrier are derived. Second, the relative permeance function of an IPMM is calculated. Third, the analytical solution of the ARMF of an IPMM is derived by applying the armature reaction magnetic field with unsaturated rotor multiplied by relative permeance function. Finally, several results of comparisons between the calculation method proposed in this paper and the finite element method are presented. Based on the calculation method proposed in this paper, the magnetic barrier’s influence on the ARMF is studied. The spatial harmonic orders and time harmonic orders of the ARMF of IPMM are revealed respectively.


Introduction
The energy crisis and environmental concerns worldwide have helped usher in the age of electric vehicles and hybrid electric vehicles.Since the interior permanent magnet motors (IPMMs) can utilize the reluctance torque, they have the advantages of high torque density, high efficiency, extended speed range.Therefore, they are widely used in these vehicles [1][2][3].The analysis of the armature reaction field is the most critical issue in the study of IPMMs since it determines the characters of torque, efficiency, vibration, and the radiated acoustic noise [4,5].
A 2D analytical model of a slotless surface permanent magnet motor (SPMM) has been developed [6][7][8][9].In Reference [10], a slotless model was developed to predict the radial component of flux density distribution; a relative permeance function was developed to consider the influences of stator slot; the calculation for the distribution of flux density in the slotted motor was carried out through multiplying the slotless model with the permeance function.The method proposed in Reference [10] has received extensive applications in the SPMM for predicting radial force [11][12][13].
However, the above theoretical analysis models of airgap magnetic field are all aimed at SPMMs.Because of the local saturation of the rotor magnetic circuit in IPMMs, radial attractive forces are mostly analyzed based on finite element method.
In References [14][15][16][17][18], the airgap flux density and electromagnetic pressure are obtained by finite element model (FEM) to investigate the effect of the combination of pole and slot on vibration and noise.While the eccentricity exists, the magnetic unbalanced force and radial force density of IPMM and SPMM are investigated in Reference [19].The magnetic field considering the saturation of the magnetic field and rotor eccentricity is analyzed by the finite element method.The radial force density and magnetic unbalanced force are calculated by the maxwell stress tensor.The calculation results are analyzed by using discrete fourier transform.In Reference [20], the FEM is used to analyze the magnetic flux density of a 120-slot-116-pole machine.The radial electromagnetic forces and tangential electromagnetic forces in airgap are then calculated by the maxwell stress tensor method.
Although the airgap flux density can be directly calculated by the finite element (FE) software, the spatial harmonic orders and time harmonic orders of the airgap flux density cannot be investigated straightforwardly.This paper provides a calculation method of the ARMF of an IPMM.First, the formulas of ARMF with unsaturated rotor are derived.Second, the relative permeance of IPMM is calculated.Third, the analytical solution of the armature reaction magnetic field (ARMF) of IPMM is derived by applying the ARMF with unsaturated rotor by relative permeance.Finally, several results of comparisons between the calculation method proposed in this paper and the finite element method are presented.Based on the calculation method proposed in this paper, the magnetic barrier's influence on the ARMF is studied.The spatial harmonic orders and time harmonic orders of the ARMF of IPMM are revealed respectively.

Analytical Model
Figure 1 shows an interior permanent magnet motor.The rotor iron is partially saturated because the permanent magnet was inserted into the rotor.The assumptions of the infinite permeable iron material in the stator and the rotor except saturation iron area are used on the analytical model in this paper.noise.While the eccentricity exists, the magnetic unbalanced force and radial force density of IPMM and SPMM are investigated in Reference [19].The magnetic field considering the saturation of the magnetic field and rotor eccentricity is analyzed by the finite element method.The radial force density and magnetic unbalanced force are calculated by the maxwell stress tensor.The calculation results are analyzed by using discrete fourier transform.In Reference [20], the FEM is used to analyze the magnetic flux density of a 120-slot-116-pole machine.The radial electromagnetic forces and tangential electromagnetic forces in airgap are then calculated by the maxwell stress tensor method.
Although the airgap flux density can be directly calculated by the finite element (FE) software, the spatial harmonic orders and time harmonic orders of the airgap flux density cannot be investigated straightforwardly.This paper provides a calculation method of the ARMF of an IPMM.First, the formulas of ARMF with unsaturated rotor are derived.Second, the relative permeance of IPMM is calculated.Third, the analytical solution of the armature reaction magnetic field (ARMF) of IPMM is derived by applying the ARMF with unsaturated rotor by relative permeance.Finally, several results of comparisons between the calculation method proposed in this paper and the finite element method are presented.Based on the calculation method proposed in this paper, the magnetic barrier's influence on the ARMF is studied.The spatial harmonic orders and time harmonic orders of the ARMF of IPMM are revealed respectively.

Analytical Model
Figure 1 shows an interior permanent magnet motor.The rotor iron is partially saturated because the permanent magnet was inserted into the rotor.The assumptions of the infinite permeable iron material in the stator and the rotor except saturation iron area are used on the analytical model in this paper.Analysis of the magnetic field in the IPMM is difficult due to the magnetic saturation in the rotor.The analytic formulas of the armature reaction magnetic field of IPMM can be derived as follows: (1) As shown in Figure 2, the armature reaction field distribution is obtained assuming a smooth airgap and without considering the interior permanent magnet.The solution of the simplified model has been discussed in the References [6,7].To apply and understand this solution more conveniently and accurately, some coefficients in this solution have been revised, the solution details and results of the model have been extended in this paper.(2) A relative permeance function is introduced to account for the effect of the magnetic barrier t in the rotor.Analysis of the magnetic field in the IPMM is difficult due to the magnetic saturation in the rotor.The analytic formulas of the armature reaction magnetic field of IPMM can be derived as follows: (1) As shown in Figure 2, the armature reaction field distribution is obtained assuming a smooth airgap and without considering the interior permanent magnet.The solution of the simplified model has been discussed in the References [6,7].To apply and understand this solution more conveniently and accurately, some coefficients in this solution have been revised, the solution details and results of the model have been extended in this paper.(2) A relative permeance function is introduced to account for the effect of the magnetic barrier t in the rotor.
As shown in Figure 2, the armature reaction field distribution is obtained assuming a smooth airgap and without considering the interior permanent magnet.The solution of the simplified model has been discussed in the References [6,7].To apply and understand this solution more conveniently and accurately, some coefficients in this solution have been revised, the solution details and results of the model have been extended in this paper.(2) A relative permeance function is introduced to account for the effect of the magnetic barrier t in the rotor.The simplified magnetic field model of a coil is shown in Figure 3.In polar coordinates, the governing Laplacian equation can be expressed as: where r is the radius of the calculation air gap, α is the circumferential angle, φ is the vector magnetic potential.The simplified magnetic field model of a coil is shown in Figure 3.In polar coordinates, the governing Laplacian equation can be expressed as: where r is the radius of the calculation air gap, α is the circumferential angle, ϕ is the vector magnetic potential.α α= 0 Assuming that the permeability of the stator core and rotor hub is infinite, the boundary conditions to be satisfied are: where J is the current density, i is the current of the coil, Rs is the stator inner radius, Rr is the rotor outer radius, b0 is the width of slot opening.Since the general solution for ϕ(α,r) is: And Assuming that the permeability of the stator core and rotor hub is infinite, the boundary conditions to be satisfied are: where J is the current density, i is the current of the coil, R s is the stator inner radius, R r is the rotor outer radius, b 0 is the width of slot opening.Since the general solution for φ(α,r) is: And where µ 0 is the permeability of the vacuum.
Based on the current sheet model, the radial component of the magnetic flux density is derived as: where B rv is the amplitude of the magnetic flux density, Q 1 is the number of stator slots that a coil is covering, Q s is the number of stator slots.

Armature Reaction Field of 1-Phase Winding
According to the superposition method, the flux density of the one phase winding is derived.
where K dυ is the winding distribution factor.For integer slot motors, K dυ can be expressed as: where q = Q s /2pm; p is the number of poles, m is the number of phases.
Taking an 8-pole 48-slot single-layer distributed winding as an example: For cos So, the Formula (12) can be simplified as: Similarly, Equation (10) can be simplified as: where B rv = NK dv B rv N is the number of series turns per phase.

ARMF of 3-Phase Winding
There are significant harmonics contained in the phase winding current waveform of a brushless dc motor.These can be expressed as a Fourier series, viz: where ω r is the electrical angular velocity, t is time, θ u is the initial phase angle.
In a brushless dc motor, the flux density produced by 3-phase winding can be derived from (15) as: where Therefore, where (u Energies 2018, 11, 2375 6 of 12

Complex Relative Permeance in Rotor Motion Coordinate
To calculate the airgap permeance per area unit, Reference [21] developed a numerical-analytical coupling model.To determine the fluctuations in the average airgap magnetic field caused by the rotor saturation of IPMM, two series of two magnetostatic FEs simulations are necessary as shown in Figure 4.

Complex Relative Permeance in Rotor Motion Coordinate
To calculate the airgap permeance per area unit, Reference [21] developed a numericalanalytical coupling model.To determine the fluctuations in the average airgap magnetic field caused by the rotor saturation of IPMM, two series of two magnetostatic FEs simulations are necessary as shown in Figure 4.While unconsidering the influence of the magnetic barrier, the radial airgap magnetic flux density and tangential airgap magnetic flux density are calculated using the model shown in Figure 4a.While considering the influence of the magnetic barrier, the radial airgap magnetic flux density and tangential airgap magnetic flux density are calculated using the model shown in Figure 4b.The real (radial) and imaginary (tangential) parts of the global airgap permeance per unit can be expressed as: where λreal and λimag are the radial and tangential components of relative permeance, Br_mb is the radial flux density with magnetic barrier, and Bθ_mb is the tangential flux density with a magnetic barrier, Br_unmb is the radial flux density without a magnetic barrier, and Bθ_unmb is the tangential flux density without a magnetic barrier.The relative permeances are calculated by (20) and (21).By analyzing the curves of Figures 5 and 6, the radial relative permeance and tangential relative permeance introduced by magnetic barrier can also be expressed as follows: where p is the number of pole pairs.While unconsidering the influence of the magnetic barrier, the radial airgap magnetic flux density and tangential airgap magnetic flux density are calculated using the model shown in Figure 4a.While considering the influence of the magnetic barrier, the radial airgap magnetic flux density and tangential airgap magnetic flux density are calculated using the model shown in Figure 4b.The real (radial) and imaginary (tangential) parts of the global airgap permeance per unit can be expressed as: where λ real and λ imag are the radial and tangential components of relative permeance, B r_mb is the radial flux density with magnetic barrier, and B θ_mb is the tangential flux density with a magnetic barrier, B r_unmb is the radial flux density without a magnetic barrier, and B θ_unmb is the tangential flux density without a magnetic barrier.The relative permeances are calculated by (20) and (21).By analyzing the curves of Figures 5 and 6, the radial relative permeance and tangential relative permeance introduced by magnetic barrier can also be expressed as follows: where p is the number of pole pairs.
Energies 2018, 11, 2375 where p is the number of pole pairs.

Complex Relative Permeance in Stator Static Coordinate
The stator of the motor is stationary, but the rotor is rotated opposite the stator in the stator static coordinate.Therefore, the spatial angle and the time are both variables in the function of the complex relative permeance.Fourier series expansions of the λreal and λimag, in the stator static coordinate are: where, an and bn is the amplitude of relative permeance.

Armature Reaction Field of IPMM
The radial airgap magnetic fields of IPMM are calculated in the stator static coordinate.( ,r) The spatial harmonic orders and time harmonic orders of the ARMF can be obtained from (26), as shown in Table 1.
Table 1.The spatial harmonic orders and time harmonic orders of the armature reaction magnetic field (ARMF).

Items Amplitude Spatial Order Time Order
Radial airgap magnetic fields unconsidering magnetic saturation

Complex Relative Permeance in Stator Static Coordinate
The stator of the motor is stationary, but the rotor is rotated opposite the stator in the stator static coordinate.Therefore, the spatial angle and the time are both variables in the function of the complex relative permeance.Fourier series expansions of the λ real and λ imag , in the stator static coordinate are: where, a n and b n is the amplitude of relative permeance.

Armature Reaction Field of IPMM
The radial airgap magnetic fields of IPMM are calculated in the stator static coordinate.
The spatial harmonic orders and time harmonic orders of the ARMF can be obtained from (26), as shown in Table 1.
Table 1.The spatial harmonic orders and time harmonic orders of the armature reaction magnetic field (ARMF).

Amplitude Spatial Order Time Order
Radial airgap magnetic fields unconsidering magnetic saturation  2 shows the IPMM's parameters.The material of the Iron core is M19-29G.To eliminate the effect of the permanent magnets on the armature-reaction field, the permanent magnets in the FEM are replaced by air.The calculation radius of airgap magnetic field is chosen close to the rotor surface to reduce the effect of the stator slots.To realize the effect of rotor saturation on airgap magnetic field, three models are adopted, which are the model unconsidering rotor saturation, the model considering rotor saturation in static rotor condition, and the model considering rotor saturation in rotating rotor condition.Finite-element simulation results are shown in Figure 7.The rotational speed is 3000 r/min, the current is 1000 A, the number of turns per slot in the stator is 1, the radius of calculation airgap is 71 mm.effect of rotor saturation on airgap magnetic field, three models are adopted, which are the model unconsidering rotor saturation, the model considering rotor saturation in static rotor condition, and the model considering rotor saturation in rotating rotor condition.Finite-element simulation results are shown in Figure 7.The rotational speed is 3000 r/min, the current is 1000 A, the number of turns per slot in the stator is 1, the radius of calculation airgap is 71 mm.

Verification of Analytical Method
As shown in Figure 8, the analytical model results are the same as the results of the finite element method without considering the magnetic barrier.As can be seen from Figure 9, while considering the magnetic barrier, there is a small error between the analytical model analysis results and the finite element calculation results, but most of them are consistent.The analytical results can reflect the spatial harmonic orders and time harmonic orders of the ARMF.

Verification of Analytical Method
As shown in Figure 8, the analytical model results are the same as the results of the finite element method without considering the magnetic barrier.As can be seen from Figure 9, while considering the magnetic barrier, there is a small error between the analytical model analysis results and the finite element calculation results, but most of them are consistent.The analytical results can reflect the spatial harmonic orders and time harmonic orders of the ARMF.

Verification of Analytical Method
As shown in Figure 8, the analytical model results are the same as the results of the finite element method without considering the magnetic barrier.As can be seen from Figure 9, while considering the magnetic barrier, there is a small error between the analytical model analysis results and the finite element calculation results, but most of them are consistent.The analytical results can reflect the spatial harmonic orders and time harmonic orders of the ARMF.

Spatial Order Distribution Characteristics
The spatial harmonic orders of the radial component of the flux density of ARMF are shown in Figure 10, which is obtained by FEM.Without considering the magnetic barrier, the spatial orders of armature reaction magnetic field are 1, 5, 7, 11… |6c ± 1|; while the magnetic barrier is considered, the spatial orders of armature reaction magnetic field are 1, 3, 5, 7, 11….|6c ± 1 ± 2n|.These are consistent with the results of Table 1.And as can be seen from Figure 10, while considering interior permanent magnet, the fundamental component of the airgap magnetic field is significantly reduced, and there is a large third harmonic component.That is because the 0th harmonic and the 1st harmonic of the relative permeance have a great influence on the ARMF due to the saturation of the rotor.

Spatial Order Distribution Characteristics
The spatial harmonic orders of the radial component of the flux density of ARMF are shown in Figure 10, which is obtained by FEM.Without considering the magnetic barrier, the spatial orders of armature reaction magnetic field are 1, 5, 7, 11 . . .|6c ± 1|; while the magnetic barrier is considered, the spatial orders of armature reaction magnetic field are 1, 3, 5, 7, 11 . . . .|6c ± 1 ± 2n|.These are consistent with the results of Table 1.And as can be seen from Figure 10, while considering interior permanent magnet, the fundamental component of the airgap magnetic field is significantly reduced, and there is a large third harmonic component.That is because the 0th harmonic and the 1st harmonic of the relative permeance have a great influence on the ARMF due to the saturation of the rotor.
armature reaction magnetic field are 1, 5, 7, 11… |6c ± 1|; while the magnetic barrier is considered, the spatial orders of armature reaction magnetic field are 1, 3, 5, 7, 11….|6c ± 1 ± 2n|.These are consistent with the results of Table 1.And as can be seen from Figure 10, while considering interior permanent magnet, the fundamental component of the airgap magnetic field is significantly reduced, and there is a large third harmonic component.That is because the 0th harmonic and the 1st harmonic of the relative permeance have a great influence on the ARMF due to the saturation of the rotor.

Frequency Distribution Characteristics
The frequency distribution characteristics of the radial magnetic flux density of the ARMF are shown in Figure 11, which is obtained by FEM.Without considering the magnetic barrier, the frequency order of the ARMF is 1 (u = 1); while the magnetic barrier is considered, the frequency orders of the armature reaction magnetic field are 1, 3, 5, 7….|u ± 2n)|.These are consistent with the results of Table 1.And as shown in Figure 11, while considering the magnetic barrier, the fundamental component of the airgap magnetic field is significantly reduced, and there are large 3rd, 5th, and 7th harmonics components.That is because the 0th harmonic, 1st, and 3rd harmonic of the relative permeance have a great influence on the ARMF due to the saturation of the rotor.

Frequency Distribution Characteristics
The frequency distribution characteristics of the radial magnetic flux density of the ARMF are shown in Figure 11, which is obtained by FEM.Without considering the magnetic barrier, the frequency order of the ARMF is 1 (u = 1); while the magnetic barrier is considered, the frequency orders of the armature reaction magnetic field are 1, 3, 5, 7 . . . .|u ± 2n)|.These are consistent with the results of Table 1.And as shown in Figure 11, while considering the magnetic barrier, the fundamental component of the airgap magnetic field is significantly reduced, and there are large 3rd, 5th, and 7th harmonics components.That is because the 0th harmonic, 1st, and 3rd harmonic of the relative permeance have a great influence on the ARMF due to the saturation of the rotor.

Conclusions
This paper provides a calculation method of the armature reaction magnetic field (ARMF) of an IPMM.The analytical model results are consistent with the finite element results.The influence of the magnetic barrier on the armature reaction field is analyzed by the proposed analytical model.The spatial harmonic orders and time harmonic orders of the ARMF of IPMM are revealed respectively.
(1) Without considering the magnetic barrier, the spatial orders of ARMF are |6c ± 1|; while the magnetic barrier is considered, the spatial orders of ARMF are |6c ± 1 ± 2n|.For spatial order distribution characteristics of an IPMM, the fundamental component of the airgap magnetic field is significantly reduced, and there is a large 3rd harmonic component.That is because the 0th harmonic and the 1st harmonic of the relative permeance have a great influence on the ARMF due to the saturation of the rotor.
(2) Without considering the magnetic barrier, the frequency orders of the ARMF is |u|; while the magnetic barrier is considered, the frequency orders of the armature reaction magnetic field are

Conclusions
This paper provides a calculation method of the armature reaction magnetic field (ARMF) of an IPMM.The analytical model results are consistent with the finite element results.The influence of the magnetic barrier on the armature reaction field is analyzed by the proposed analytical model.The spatial harmonic orders and time harmonic orders of the ARMF of IPMM are revealed respectively.
(1) Without considering the magnetic barrier, the spatial orders of ARMF are |6c ± 1|; while the magnetic barrier is considered, the spatial orders of ARMF are |6c ± 1 ± 2n|.For spatial order distribution characteristics of an IPMM, the fundamental component of the airgap magnetic field is

2. 1 .
The Armature Reaction Field Distribution of the Simplified Model 2.1.1.Magnetic Field Produced by a Coil

Figure 3 .
Figure 3. Distributed current sheet model of coil.

Figure 3 .
Figure 3. Distributed current sheet model of coil.

Figure 4 .
Figure 4. Finite element (FE) simulations of airgap permeance (a) motor without magnetic barrier (b) motor with a magnetic barrier.

3 .
Comparison with the Results of Finite-Element Calculation Take an 8-pole-48-slot IPMM for example, a finite element model is established by the maxwell.Which is a kind of finite element software.The finite element model is used to verify the accuracy of the analytical model proposed in this paper.Table

Figure 7 .
Figure 7. Finite-element simulation model.(a) motor without magnetic barrier; (b) motor with a magnetic barrier.

Figure 7 .
Figure 7. Finite-element simulation model.(a) motor without magnetic barrier; (b) motor with a magnetic barrier.

Figure 9 .
Figure 9. Radial magnetic flux density of armature reaction field considering the magnetic barrier.

Figure 10 .
Figure 10.Spatial order distribution characteristics of the radial magnetic flux density of armature reaction field (Calculation with finite element model (FEM)).(a) radial flux density; (b) Fourier transform of radial flux density.

Figure 10 .
Figure 10.Spatial order distribution characteristics of the radial magnetic flux density of armature reaction field (Calculation with finite element model (FEM)).(a) radial flux density; (b) Fourier transform of radial flux density.

Figure 11 .
Figure 11.Frequency distribution characteristics of radial magnetic flux density of armature reaction field (Calculation with FEM).(a) radial flux density; (b) Fourier transform of radial flux density.

Figure 11 .
Figure 11.Frequency distribution characteristics of radial magnetic flux density of armature reaction field (Calculation with FEM).(a) radial flux density; (b) Fourier transform of radial flux density.