Next Article in Journal
Comprehensive Survey of Various Energy Storage Technology Used in Hybrid Energy
Next Article in Special Issue
Robust and Multi-Objective Pareto Design of a Solenoid
Previous Article in Journal
Wearable IoTs and Geo-Fencing Based Framework for COVID-19 Remote Patient Health Monitoring and Quarantine Management to Control the Pandemic
Previous Article in Special Issue
Application of Particle Swarm Optimization in the Design of an ICT High-Voltage Power Supply with Dummy Primary Winding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Study and Comparison of Electromagnetic Characteristics of 8-Pole 9-Slot and 8-Pole 12-Slot Permanent Magnet Synchronous Machines Considering Rotor Eccentricity

Department of Electrical Engineering, ChungNam National University, Daejeon 34134, Korea
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(16), 2036; https://doi.org/10.3390/electronics10162036
Submission received: 27 July 2021 / Revised: 17 August 2021 / Accepted: 19 August 2021 / Published: 23 August 2021
(This article belongs to the Special Issue Robust Design Optimization of Electrical Machines and Devices)

Abstract

:
In this study, a magnetic field is analyzed using an analytical method and compared with the electromagnetic characteristics of 8-pole 9-slot and 8-pole 12-slot permanent magnet synchronous machines considering rotor eccentricity. The magnetic flux density and back electromotive force (EMF) are derived using perturbation theory and electromagnetic theory. First, the Fourier modeling of a permanent magnet is performed through magnetization modeling, and two analysis regions are set based on several assumptions for applying the analytical method. Accordingly, the governing equations of the analysis regions are derived in the form of Poisson and Laplace equations. In addition, the undefined coefficients of the general solutions are calculated through general solutions and appropriate boundary conditions, and the magnetic flux density and back EMF of the air gap region are derived based on the definition of the magnetic vector potential. The results obtained using the analytical method are compared with the finite element method and experimental results. In addition, we perform a torque analysis considering rotor eccentricity and analyze the torque ripple based on rotor eccentricity for two cases involving the pole/slot combination.

1. Introduction

Permanent magnet synchronous machines (PMSMs) are becoming popular as a key technology for applications such as home appliances, industrial tools, and electric vehicles because of their high efficiency, high power density, and low maintenance cost; in addition, interest in resolving motor malfunction has increased. Data provided in [1,2,3] indicate that 41% of motor faults are bearing faults, 37% are stator faults, 10% are rotor faults, and 12% are other faults. One of the main causes of malfunction is rotor eccentricity, wherein the center of the rotor axis deviates from the center of the stator, resulting in a non-uniform air gap. Static eccentricity, which is a type of rotor eccentricity, is a condition where the position of the minimum radial air gap is fixed [2,3,4,5]. This can be caused by stator core ovality, incorrect positioning of the stator core, or bearing at commissioning or following a repair, and its level does not change over time. Furthermore, the magnetic flux density in the air gap is an important characteristic of the machine performance. Therefore, the impact of rotor eccentricity on the magnetic field distribution must be analyzed for predicting the characteristics of PMSMs [6,7]. The finite element method (FEM) and analytical methods are employed to design PMSMs or perform characteristic analysis. The FEM is a numerical process that can be performed using a commercial tool [8]. However, it is necessary to learn to use each commercial tool, which is disadvantageous as the experience of the designer is necessitated for an accurate analysis. In the analytical method, electromagnetic field characteristics are analyzed using Maxwell’s equation. Magnetization modeling is applied to the design of the machine using the Fourier series and partial differential equations derived using the magnetic vector potential. It must be preceded by deriving a solution that considers various boundary conditions. Several studies are being conducted because the analytical method can rapidly predict the characteristic variation based on the design parameters [9].
The magnetic field is analyzed using an analytical method and compared with the electromagnetic characteristics of 8-pole 9-slot and 8-pole 12-slot PMSMs considering rotor eccentricity. The magnetic flux density and back-electromotive force (EMF) are derived using perturbation theory and electromagnetic theory. The results obtained using the analytical method are compared with the FEM and experimental results. In addition, we performed a torque analysis based on rotor eccentricity and analyzed the torque ripple based on rotor eccentricity for two cases involving the pole/slot combination.

2. Electromagnetic Analysis of a PMSM Using Analytical Method

2.1. Analytical Model

Figure 1 shows the analysis models of 8-pole 9-slot and 8-pole 12-slot PMSMs. First, the analysis regions are defined for applying the analytical method, and the governing equations are derived using electromagnetic theory regarding each region. The magnetic flux density can be derived through the general solutions of the defined governing equations and the boundary conditions of each region. To apply the analytical method, the analysis model is simplified to the slotless model, the relative permeability of the iron core is infinite, the relative permeability of the permanent magnets (PMs) is equal to that of air, and the eddy effect is disregarded to apply the analytical method. Figure 2 shows the simplified model for applying the analytical method. The analysis regions are classified as an air gap region (Region I) and PM region (Region II), where rotor eccentricity is determined based on the stator coordinate (r,θ), and rotor coordinate (ξ,ψ). Os and Or are the centers of the stator and rotor, respectively; Rs, Rm, and Rr are the radii of the stator inner, PM surface, and rotor inner, respectively. The relationship between the parameters of the r-θ coordinate and ξ-ψ coordinate is organized based on Taylor’ expansion, as follows [9].
ξ = r ε cos ( θ ϕ ) + O ( ε 2 ) ψ = ( θ w t ) + ε r sin ( θ ϕ ) + O ( ε 2 ) ,    ε = e c × g
where r is the length from the stator center to the observation point, ε is the length of rotor eccentricity, ec is the eccentricity ratio, and g is the nominal air gap length. O(ε2) is the second order of the perturbation term. Only the first-order term of the perturbation term is considered, and the second-order term is disregarded in this study, since rotor eccentricity is usually small in comparison to the air gap length.

2.2. Magnetic Field Analysis Using Analytical Method

To derive the magnetic flux density for PMs, Fourier modeling must be performed. Fourier modeling for PMs is expressed as follows:
M = n = 1 , 3 , 5 M r n cos ( n p θ ) i r + M θ n sin ( n p θ ) i θ
where M is the magnetization vector, Mrn and Mθn are Fourier coefficients of the radial component and tangential component of the r- and θ-direction, respectively. n and p are the space harmonic order and number of pole pairs, respectively. Rotor eccentricity can be treated as one of kind of perturbation phenomenon. Perturbation theory is necessitated to consider rotor eccentricity. Perturbation theory is a theory that expresses the solution of a problem that cannot be solved analytically as a Taylor series of parameters that can be considered extremely small. When this is applied to the magnetic vector potential, it is expressed as follows [9]:
A z 1 ( r , θ , ε ) = A z 1 ( 0 ) ( r , θ ) + ε A z 1 ( 1 ) ( r , θ ) + ε 2 A z 1 ( 2 ) ( r , θ ) + A z 2 ( r , θ , ε ) = A z 2 ( 0 ) ( r , θ ) + ε A z 2 ( 1 ) ( r , θ ) + ε 2 A z 2 ( 2 ) ( r , θ )
where Az1 and Az2 are the magnetic vector potential of Region I and II, respectively. Azn(0) and Azn(1) are the zeroth and first orders of perturbation, respectively. Azn(0) occurs when the rotor is not eccentric. Azn(1) is generated by the rotor eccentricity effect without PMs or current excitation. The reason why even the first term is considered is that more accurate results can be derived when the second term is considered, but the calculation process becomes complicated, and reliable results can be obtained even when only the first term is considered. Therefore, the governing equation for each region is expressed as Poisson and the Laplace equations.
2 A z 1 ( 0 ) = 2 A z 1 ( 0 ) r 2 + 1 r A z 1 ( 0 ) r q 2 r 2 A z 1 ( 0 ) = 0
2 A z 2 ( 0 ) r 2 + 1 r A z 2 ( 0 ) r q 2 r 2 A z 2 ( 0 ) = μ 0 × M
2 A z 1 ( 1 ) r 2 + 1 r A z 1 ( 1 ) r q 2 r 2 A z 1 ( 1 ) = 0
2 A z 2 ( 1 ) r 2 + 1 r A z 2 ( 1 ) r q 2 r 2 A z 2 ( 1 ) = 0
The general solution of each region is derived from the Cauchy–Euler equation. The general solutions of each analysis region can be obtained as follows:
A z 1 ( 0 ) = n = 1 , 3 , 5 ( A 1 r n p + B 1 r n p ) sin ( n p θ ) i z
A z 2 ( 0 ) = n = 1 , 3 , 5 ( A 2 r n p + B 2 r n p + μ 0 r n p M n ( n p ) 2 1 ) sin ( n p θ ) i z
A z 1 ( 1 ) = ε n = 1 , 3 , 5 ( W 1 r n p 1 + X 1 r n p + 1 ) sin [ ( n p 1 ) θ + ϕ ] + ( Y 1 r n p + 1 + Z 1 r n p 1 ) sin [ ( n p + 1 ) θ ϕ ] i z
A z 2 ( 1 ) = ε n = 1 , 3 , 5 ( W 2 r n p 1 + X 2 r n p + 1 ) sin [ ( n p 1 ) θ + ϕ ] + ( Y 2 r n p + 1 + Z 2 r n p 1 ) sin [ ( n p + 1 ) θ ϕ ] i z
The effect of the rotor eccentricity can be observed at the rotor yoke surface and the interface between the PMs and air gap in the r-θ coordinate. In the ξ-ψ coordinate, the radius corresponding to the surface and interface are Rr and Rm, respectively. By applying ρ = Rr and Rm in Equation (1), the following equations can be obtained:
f R r ( r , θ , ε ) = r ε cos ( θ ϕ ) R r
f R m ( r , θ , ε ) = r ε cos ( θ ϕ ) R m
The normal direction vector of the rotor core surface and PM surface can be derived by applying the gradient (∇) to Equations (12) and (13).
n R m = f R m = i r + ε r sin ( θ ϕ ) i θ
n R s = f R s = i r + ε r sin ( θ ϕ ) i θ
The boundary conditions can be derived using Equations (14) and (15) and electromagnetic theory. The boundary conditions at the rotor core surface and PM–air gap interface can be expressed as follows:
n R r × H 2 = 0
n R m × ( H 1 H 2 ) = 0
n R m ( B 1 B 2 ) = 0
The boundary conditions of the zeroth and first order can be derived by applying × A = B to Equation (3) and substituting it into Equations. (16)–(18). A1, B1, A2, B2, W1, X1, Y1, Z1, W2, X2, Y2, and Z2 of Equations (8)–(11) can be calculated using the derived boundary conditions. The magnetic flux density of the r- and θ-directions can be express as follows [6]:
B r = 1 r A θ , B θ = A r
The zeroth order of the perturbation term implies no rotor eccentricity, and the first-order term refers to the amount of change in magnetic flux density based on rotor eccentricity; therefore, the final magnetic flux density is the sum of the zeroth- and first-order magnetic flux density. The magnetic flux density considering rotor eccentricity can be expressed as:
B = B ( 0 ) + ε B ( 1 )

2.3. Back-EMF

To obtain the back-EMF, the magnetic flux linking the coil must be determined as follows [10]:
Φ c o i l , ( 1 , , Q ) = θ c / 2 θ c / 2 B d S ,      E ( t ) = N c d ( Φ 1 + + Φ Q ) d t
where B is the magnetic flux density vector in the stator surface region, and dS is an element of the coil surface area vector. Q is the slot number/phase. θ c denotes the mechanical angle per slot. The back-EMF can be obtained using Faraday’s law. Nc is the number of turns in each coil. Figure 3 is an example explaining the calculation of the back-EMF of the 8-pole 9-slot. The back-EMFs in other coils of the same phase are not necessarily similar because of the eccentricity. They must be calculated individually by performing an appropriate shift.

3. Results and Discussion

Figure 4 shows the experimental set. Table 1 shows the design specifications. The stator and rotor were made of 50PN470, whereas the permanent magnet was made of N42SH; the experiment was performed at room temperature. The servo motor was installed on the opposite side of the test motor to measure the back-EMF. Figure 5a shows the magnetic flux density in the air gap when no eccentricity occurred, and Figure 5b is the variation in the magnetic flux density according to rotor eccentricity. Figure 6 shows the magnetic flux density in the air gap when 25% rotor eccentricity occurred. A 25% eccentricity is based on nominal air gap length. Since the air gap of the analysis model is 4 mm, it means that 1 mm eccentricity occurs. The analytical results were verified by comparing them with the FEM results. FEM analysis was performed using ANSYS Electronics 2020R2. It was confirmed that the results agreed well. Figure 7 is the mesh plot of non-simplified analysis models. Non-simplified model in Figure 8 and Figure 9 are the slotted models. FEM analysis and experiments are performed using slotted models. Figure 8a,b show the back-EMFs of the 8-pole 9-slot PMSM when no rotor eccentricity and 25% rotor eccentricity occurred, respectively. The analytical results of the back-EMFs of the 8-pole 9-slot PMSM were compared with the FEM and experimental results. It was observed that the results exhibited good agreement with each other. In the 8-pole 9-slot model, the winding arrangement distributed asymmetrically to one side, as shown in Figure 8b; therefore, if eccentricity occurs, the result of the unbalanced back-EMF can be confirmed. Figure 9a,b show the back-EMF of the 8-pole 12-slot PMSM based on rotor eccentricity. As shown in Figure 9b, the winding arrangement of the 12-slot stator was distributed symmetrically. Consequently, the effect of rotor eccentricity was nullified, and the results of the cases with and without rotor eccentricity were compared. Table 2 and Table 3 summarize the analysis results. The error of the analytical results is confirmed to be within 5%. Figure 10a, b show the results of the output torque and torque ripple analysis based on the eccentricity of the 8-pole 9-slot and 8-pole 12-slot PMSMs. When comparing those PMSMs, it was observed that the cogging torque of the 8-pole 12-slot PMSM was larger than that of the 8-pole 9-slot PMSM; therefore, the torque ripple ratio of the 8-pole 12-slot PMSM was larger than that of the 8-pole 9-slot under a low input current. However, in the 8-pole 9-slot model, the increasing rate of the torque ripple based on eccentricity was larger than that of the 8-pole 12-slot owing to the asymmetric winding arrangement. Under the low output power condition, the 8-pole 9-slot PMSM indicated a sinusoidal back-EMF, resulting in a small torque ripple ratio in the output torque. However, as the output power increased, the 8-pole 12-slot PMSM was less affected by the rotor eccentricity.

4. Conclusions

In this study, we performed a magnetic field analysis based on rotor eccentricity using an analytical method, and the torque characteristics of a two pole/slot combination were compared. Based on perturbation theory and electromagnetic theory, the magnetic flux density and back-EMF in the air gap region were derived, and the analytical result was verified by comparing it with the FEM and the experimental results. Both models with slotted and simplified models without slots were analyzed using the FEM. We conformed that even if the magnetic flux density is predicted as slotless and the back-EMF is derived using it, a slotted model can obtain similar results. Subsequently, torque analysis was performed based on rotor eccentricity to analyze the variation in the torque ripple for the two pole/slot combinations. The analysis results from this study will benefit the design of motors.

Author Contributions

J.-Y.C.: conceptualization, review and editing; H.-K.L.: analysis, original draft preparation; T.-K.B.: magnetic field calculation; J.-I.L.: derivation of governing equation; J.-H.W.: experiment; H.-S.S.: torque analysis; I.-J.Y.: review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Basic Research Laboratory (BRL) of the National Research Foundation (NRF-2020R1A4A2002021) funded by the Korean government.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goktas, T.; Zafarani, M.; Lee, K.W.; Akin, B.; Sculley, T. Comprehensive Analysis of Magnet Defect Fault Monitoring Through Leakage Flux. IEEE Trans. Magn. 2017, 53, 1–10. [Google Scholar] [CrossRef]
  2. Prieto, M.D.; Espinosa, A.G.; Ruiz, J.R.; Urresty, J.C.; Ortega, J.A. Feature extraction of demagnetization faults in PMSMs based on box-counting fractal dimension. IEEE Trans. Ind. Electron. 2011, 58, 1594–1605. [Google Scholar] [CrossRef]
  3. Hyun, D.; Hong, J.; Lee, S.B.; Kim, K.; Wiedenbrug, E.J.; Teska, M.; Nandi, S.; Chelvan, I.T. Automated Monitoring of Airgap Eccentricity for Inverter-Fed Induction Motors Under Standstill Conditions. IEEE Trans. 2011, 47, 1257–1266. [Google Scholar]
  4. Park, Y.; Fernandez, D.; Lee, S.B.; Hyun, D.; Jeong, M.; Kommuri, S.K.; Cho, C.; Reigosa, D.; Briz, F. On-line detection of rotor eccentricity for PMSMs based on Hall-effect field sensor measurements. In Proceedings of the 2017 IEEE Energy Conversion Congress and Exposition (ECCE), Cincinnati, OH, USA, 1–5 October 2017; pp. 4678–4685. [Google Scholar]
  5. Thomson, W.; Barbour, A. On-line current monitoring and application of a finite element method to predict the level of static airgap eccentricity in three-phase induction motors. IEEE Trans. Energy Convers. 1998, 13, 347–357. [Google Scholar] [CrossRef]
  6. Amirat, Y.; Benbouzid, M.E.H.; Bensaker, B.; Wamkeue, R. Condition monitoring and fault diagnosis in wind energy conversion systems: A review. In Proceedings of the IEEE PEMDC Conference, Antalya, Turkey, 3–5 May 2007. [Google Scholar]
  7. Lu, B.; Li, Y.; Wu, X.; Yang, Z. A review of recent advances in wind turbine condition monitoring and fault diagnosis. No. 24026. In Proceedings of the IEEE Conference on Power Electronics and Machines in Wind Applications (PEMWA), Lincoln, NE, USA, 24–26 June 2009; pp. 1–7. [Google Scholar]
  8. Yang, W.; Tavner, P.J.; Wilkinson, M. Wind Turbine Condition Monitoring and Fault Diagnosis Using Both Mechanical and Electrical Signatures. In Proceedings of the 2008 IEEE/ASME International Conference on Advanced Intelligent Mechatronics, Xi’an, China, 2–5 July 2008; pp. 1296–1301. [Google Scholar] [CrossRef]
  9. Kim, U.; Lieu, D.K. Magnetic field calculation in permanent magnet motors with rotor eccentricity: Without slotting effect. IEEE Trans. Magn. 1998, 34, 2242–2252. [Google Scholar] [CrossRef]
  10. Rahideh, A.; Korakianitis, T. Analytical open-circuit magnetic field distribution of slotless brushless permanent magnet machines with rotor eccentricity. IEEE Trans. Magn. 2011, 47, 4791–4808. [Google Scholar] [CrossRef]
Figure 1. Analysis model: (a) 8 pole, 12 slot; (b) 8 pole, 9 slot.
Figure 1. Analysis model: (a) 8 pole, 12 slot; (b) 8 pole, 9 slot.
Electronics 10 02036 g001
Figure 2. Simplified model for using analytical method.
Figure 2. Simplified model for using analytical method.
Electronics 10 02036 g002
Figure 3. Winding distribution of 8-pole 9-slot.
Figure 3. Winding distribution of 8-pole 9-slot.
Electronics 10 02036 g003
Figure 4. Experimental set: (a) experiment configuration; (b) manufactured model.
Figure 4. Experimental set: (a) experiment configuration; (b) manufactured model.
Electronics 10 02036 g004
Figure 5. Flux density: (a) without rotor eccentricity; (b) variation in radial flux density according to rotor eccentricity.
Figure 5. Flux density: (a) without rotor eccentricity; (b) variation in radial flux density according to rotor eccentricity.
Electronics 10 02036 g005
Figure 6. Flux density with 25% rotor eccentricity at the air gap: (a) radial flux density; (b) tangential flux density.
Figure 6. Flux density with 25% rotor eccentricity at the air gap: (a) radial flux density; (b) tangential flux density.
Electronics 10 02036 g006
Figure 7. Mesh plot of analysis models: (a) 8−pole 9−slot; (b) 8−pole 12−slot.
Figure 7. Mesh plot of analysis models: (a) 8−pole 9−slot; (b) 8−pole 12−slot.
Electronics 10 02036 g007
Figure 8. Back−EMF of 8−pole 9−slot PMSM: (a) without rotor eccentricity; (b) with 25% eccentricity.
Figure 8. Back−EMF of 8−pole 9−slot PMSM: (a) without rotor eccentricity; (b) with 25% eccentricity.
Electronics 10 02036 g008
Figure 9. Back-EMF of 8−pole 12-slot PMSM: (a) without rotor eccentricity; (b) with 25% eccentricity.
Figure 9. Back-EMF of 8−pole 12-slot PMSM: (a) without rotor eccentricity; (b) with 25% eccentricity.
Electronics 10 02036 g009
Figure 10. Torque characteristics: (a) 8−pole 12−slot and (b) 8−pole 9−slot.
Figure 10. Torque characteristics: (a) 8−pole 12−slot and (b) 8−pole 9−slot.
Electronics 10 02036 g010
Table 1. Design specification of analysis model.
Table 1. Design specification of analysis model.
ParameterValueUnit
Number of slots/poles9/8, 12/8mm
Outer radius of stator75mm
Inner radius of stator47mm
Outer radius of rotor43mm
Thickness of PMs5mm
Axial length30mm
Magnet remanence1.27T
Rated speed1000Rpm
Pole arc ratio0.9
Table 2. Analysis results of back-EMF without rotor eccentricity.
Table 2. Analysis results of back-EMF without rotor eccentricity.
ParameterAnalyticalFEMMeasurement
8p9s8p12s8p9s8p12s8p9s8p12s
Number of meshes--31,49630,541--
Analysis time (s)0.110.121123--
Back-EMF (Vmax)17.6318.4517.5218.3517.6819.05
Error (%)0.283.140.93.67--
Table 3. Analysis results of back-EMF with 25% rotor eccentricity.
Table 3. Analysis results of back-EMF with 25% rotor eccentricity.
ParameterAnalyticalFEMMeasurement
8p9s8p12s8p9s8p12s8p9s8p12s
Number of meshes--31,83933,910--
Analysis time (s)0.110.121225
Back-EMF (Vmax)18.3618.4518.218.3518.619.04
Error (%)1.293.12.153.62--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, H.-K.; Bang, T.-K.; Lee, J.-I.; Woo, J.-H.; Shin, H.-S.; Yoon, I.-J.; Choi, J.-Y. Analytical Study and Comparison of Electromagnetic Characteristics of 8-Pole 9-Slot and 8-Pole 12-Slot Permanent Magnet Synchronous Machines Considering Rotor Eccentricity. Electronics 2021, 10, 2036. https://doi.org/10.3390/electronics10162036

AMA Style

Lee H-K, Bang T-K, Lee J-I, Woo J-H, Shin H-S, Yoon I-J, Choi J-Y. Analytical Study and Comparison of Electromagnetic Characteristics of 8-Pole 9-Slot and 8-Pole 12-Slot Permanent Magnet Synchronous Machines Considering Rotor Eccentricity. Electronics. 2021; 10(16):2036. https://doi.org/10.3390/electronics10162036

Chicago/Turabian Style

Lee, Hoon-Ki, Tae-Kyoung Bang, Jeong-In Lee, Jong-Hyeon Woo, Hyo-Seob Shin, Ick-Jae Yoon, and Jang-Young Choi. 2021. "Analytical Study and Comparison of Electromagnetic Characteristics of 8-Pole 9-Slot and 8-Pole 12-Slot Permanent Magnet Synchronous Machines Considering Rotor Eccentricity" Electronics 10, no. 16: 2036. https://doi.org/10.3390/electronics10162036

APA Style

Lee, H. -K., Bang, T. -K., Lee, J. -I., Woo, J. -H., Shin, H. -S., Yoon, I. -J., & Choi, J. -Y. (2021). Analytical Study and Comparison of Electromagnetic Characteristics of 8-Pole 9-Slot and 8-Pole 12-Slot Permanent Magnet Synchronous Machines Considering Rotor Eccentricity. Electronics, 10(16), 2036. https://doi.org/10.3390/electronics10162036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop