Next Article in Journal
Numerical Prediction of Turbulent Drag Reduction with Different Solid Fractions and Distribution Shapes over Superhydrophobic Surfaces
Next Article in Special Issue
Quantitative Comparisons of Outer-Rotor Permanent Magnet Machines of Different Structures/Phases for In-Wheel Electrical Vehicle Application
Previous Article in Journal
Wind Power Generation Forecast Based on Multi-Step Informer Network
Previous Article in Special Issue
An Improved PLL-Based Speed Estimation Method for Induction Motors through Harmonic Separation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Online Extraction Method of Slot Harmonics for Multiphase Induction Motor

College of Electrical Engineering, Qingdao University, Qingdao 266071, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(18), 6643; https://doi.org/10.3390/en15186643
Submission received: 20 July 2022 / Revised: 29 August 2022 / Accepted: 2 September 2022 / Published: 11 September 2022
(This article belongs to the Special Issue Advanced Design and Control of Multiphase Machines)

Abstract

:
The accurate extraction and analysis of slot harmonics caused by slotting in an induction motor are important for the motor’s performance evaluation and state monitoring. However, the frequency distribution of rotor slot harmonics (RSHs) varies along with the operating states of the motor, such as motor speed and slip ratio, and the voltage and current signals of the motor only contain small-amplitude RSHs compared with other harmonics; both make it difficult to extract and analyze the RSHs accurately online. While offline extraction and filters with constant parameters are mainly utilized in available works, a novel adaptive extraction method for RSHs in a multiphase induction motor is proposed here to realize online RSH extraction under different speed and load conditions. In this paper, the RSHs in the multiphase induction motors are firstly modeled by using the magnetic potential permeability method, and the influence of a skewed rotor on RSHs is analyzed through a multisection method. Then, an adaptive extraction method of RSHs is proposed, which can effectively realize the online processing of RSHs of stator current. Finally, the experimental platform of a nine-phase induction motor has been used to verify the effectiveness of the proposed method under different speeds and load conditions, with a relative error of less than 1% in identifying the RSH frequency distribution.

1. Introduction

Multiphase induction motors (IMs) have been widely used in electric transportation applications, such as electric vehicle drive, rail transport and ship propulsion [1]. Although the multiphase IM features a reliable structure, lower phase voltage level [2], convenience in maintenance [3] and good tolerant control capability [4,5], the working conditioning measurement is still vital for the healthy operation of IMs, in which the real-time extraction and analysis of voltage or current harmonics are important for providing essential or inherent information for motor control or monitoring.
On one hand, slot harmonics are generally regarded as disturbances for motor operation; they not only cause motor cost but also lead to torque ripples [6] and current fluctuation [7]. When the rotor is eccentric, the slot harmonic will also affect the magnetic pull force [8]. It is also necessary to optimize the stator and rotor slot number and slot dimensions of the motor within its design stage [9]. Fractional slot [10], closed-slot [11], skewed slot [12] and assisting magnetic pole [13] are all feasible methods for reducing slot harmonics, but the above methods cannot completely eliminate tooth harmonics. Therefore, slot harmonics need to be modeled and analyzed, which is usually combined with motor design and iron loss analysis [14] and efficiency estimation [15].
On the other hand, the time and frequency properties of slot harmonics are effective indicators that directly reveal the working condition of motor and can be used for fault detection, including broken rotor bar fault diagnosis [16] due to the irregular frequency distribution of rotor slot harmonics (RSHs), sensor-less rotor temperature estimation [17] by evaluating the rotor loss through RSHs, and the detection of stator winding faults [18], as well as the measurement of motor speed in real time because the slip ratio is inherently embedded in the frequencies of RSHs [19]. However, the modeling of slot harmonics cannot comprehensively consider the actual operating conditions such as temperature [20] and eccentric non-ideal conditions [21], as well as the other manufacturing errors [22], so the online extraction and analysis of motor slot harmonics are actually of priority before they can be further analyzed and utilized for the motor state monitoring and motor control methods. However, the motor tooth harmonics, especially the frequency distribution of RSHs, change with the motor running state, which makes it difficult to extract them online. In fact, the offline extraction and filters with constant parameters are mainly utilized in available works, which restricts the performance and accuracy of the online evaluation of RSH properties of induction motors.
In this paper, we introduce an adaptive rotor slot harmonic extraction algorithm for multiphase induction motors. The time and frequency properties of RSHs are firstly modeled and analyzed, especially under the transients of load or speed change. Then an adaptive online extraction and analysis method for RSHs is proposed to cope with the current harmonics of the motor in operation, especially for different speed and load conditions. Finally, a nine-phase induction motor platform was utilized to verify the proposed adaptive online extraction and analysis method for primary RSHs (PSHs).

2. The Model and Analysis of Rotor Slot Harmonics for Multiphase Induction Motors

2.1. The Analytical Model of RSHs for Multiphase IMs

The rotor slot harmonics are generated by the cogging effect and are difficult to completely eliminate by skewed slot and other methods. The open slot on the stator and rotor has a modulation effect on the air gap magnetic field, as shown in Figure 1, and is furtherly reflected in the motor voltage or current.
To obtain the theoretical model of RSHs for multiphase IMs, the magnetic potential and permeability method is used here. For a multiphase induction motor in normal operation, its synthetic magnetic potential of stator and rotor windings is first listed as follows:
f ( θ , t ) = f 0 ( θ , t ) + ν f ν ( θ , t ) + μ f μ ( θ , t )
where θ is mechanical rotor position, t is the time, f0 is the fundamental component, fv is the vth magnetic potential of stator windings and fμ is the μth magnetic potential of rotor windings.
f 0 ( θ , t ) = F 0 cos ( p θ ω 1 t φ 0 ) f ν ( θ , t ) = F ν cos ( ν θ ω v t φ v ) f μ ( θ , t ) = F μ cos ( μ θ ω μ t φ μ )
in which p is the pole pair. ω1, ωv and ωμ and φ1, φv and φμ are respectively the angular frequencies and phases of each component. F0, Fv and Fμ are their magnetic potential amplitudes.
As an example, if the number of slots per pole per phase is an integer for a nine-phase induction motor, the harmonic number of stator windings is
ν = ( 18 k 1 + 1 ) p ,       k 1 ± 1 ,   ± 2 ,  
If k1 is an integral multiple of the number of slots per pole per phase q1, such as k1 = ±q1, ±2q1, ±3q1 …, the order of stator slot harmonics will be obtained as
ν z = ( 18 q 1 k 1 + 1 ) p = k 1 Z 1 + p     k 1 = ± 1 , ± 2 ,
where Z1 is the stator slot number.
Then, the magnetic potential harmonic order for a squirrel-cage rotor winding can be defined as
μ = k 2 Z 2 + ν     k 2 = ± 1 , ± 2 ,
where Z2 is the rotor slot number.
Meanwhile, the fundamental magnetic field in the air gap will induce the electromotive force of rotor slot harmonics with the frequencies as
μ z = k 2 Z 2 + p     k 2 = ± 1 , ± 2 ,
As the magnetic potential has been modeled, the characteristics of air gap permeance due to stator and rotor slots should also be analyzed. With slots on the stator and rotor, the air gap permeability can be approximately expressed as
λ ( θ , t ) = Λ 0 ( 1 + k 1 λ k 1 Λ 0 ) ( 1 + k 2 λ k 2 Λ 0 ) Λ 0 + k 1 λ k 1 + k 2 λ k 2 + k 1 k 2 λ k 1 λ k 2
where Λ0 is the constant or DC part, Λk1 is the k1th magnetic permeance harmonics only with stator slots, Λk2 is the k2th magnetic permeance harmonics only with rotor slots and Λk1Λk2 is due to the relative effect between stator and rotor slots.
Λ 0 = μ 0 δ K c
Λ k 1 = Λ 0 ( 1 ) k 1 ( K c 1 1 ) | sin k 1 K c 1 1 K c 1 π k 1 K c 1 1 K c 1 π |
Λ k 2 = Λ 0 ( 1 ) k 2 ( K c 2 1 ) | sin k 2 K c 2 1 K c 2 π k 2 K c 2 1 K c 2 π |
λ k 1 k 2 = Λ k 1 Λ k 2 Λ 0 cos [ ( ± k 1 Z 1 k 2 Z 2 ) θ k 2 Z 2 p ( 1 s ) ω t ]
where δ is the nominal air gap length; Kc, Kc1 and Kc2 are the Karter coefficients, which are related to the dimensions of stator and rotor slots.
According to the properties of magnetic potential and permeance discussed above, the air gap flux distribution can be calculated as follows:
b ( θ , t ) = f ( θ , t ) λ ( θ , t ) = { f 0 ( θ , t ) + ν f ν ( θ , t ) + μ f μ ( θ , t ) } { Λ 0 + k 1 λ k 1 + k 2 λ k 2 } { F 0 Λ 0 cos ( p θ ω 1 t φ 0 ) + ν F ν Λ 0 cos ( ν θ ω 1 t φ 1 ) + μ F μ Λ 0 cos ( μ θ ω μ t φ 2 ) } + { k 1 F 0 Λ k 1 2 cos [ ( ± k 1 Z 1 + p ) θ ω 1 t φ 0 ) ] } + { k 2 F 0 Λ k 2 2 cos [ ( ± k 2 Z 2 + p ) θ ( ± k 2 Z 2 p ( 1 s ) + 1 ) ω 1 t φ 0 ) ] }
in which the Λk1Λk2 component is omitted for simplicity; s is the slip ratio. The rotor slot harmonics for windings is
b f μ z ( θ , t ) = μ z B f μ z cos ( μ z θ ω μ t φ 2 ) B f μ z = F μ z Λ 0 = p μ K δ I 2 I 0 B δ
and the rotor slot harmonics for magnetic permeance is
b Λ μ z ( θ , t ) = μ z B Λ μ z cos ( μ z θ ω μ t φ 0 ) B Λ μ z = 1 2 F 0 K δ Λ k 2 = 1 2 Λ k 2 Λ 0 B δ
where μ z = k 2 Z 2 + p k 2 = ± 1 , ± 2 ,
Taking the primary rotor slot harmonic (PSH) as an example, its time domain expression can be divided into two parts: the first one Br1 is generated by the fundamental magnetic potential of the motor, along with the magnetic conductivity due to slot cogging effect, and the other one Br2 is caused by the interaction between the RSH induced magnetic potential and the DC component of air gap magnetic conductivity. According to Equations (12)–(15), and by combining the components with the same frequency, the primary rotor slot harmonics (PSHs) can be modeled as follows:
B r ( θ , t ) = B r 1 ( θ , t ) + B r 2 ( θ , t ) = ( F r ( Z 2 / p 1 ) Λ 0 + 1 2 F 0 Λ r ) cos [ ( Z 2 p ) θ ( Z 2 p ( 1 s ) 1 ) ω 1 t ] + ( F r ( Z 2 / p + 1 ) Λ 0 + 1 2 F 0 Λ r ) cos [ ( Z 2 + p ) θ ( Z 2 p ( 1 s ) + 1 ) ω 1 t ] = B r s h 1 ( θ , t ) + B r s h 2 ( θ , t )
where Bt1−(θ,t) and Bt1+(θ,t) are two side frequency components of PSH, F r ( Z 2 / p 1 ) and F r ( 2 , p + 1 ) are amplitudes of rotor cogging MMFs and Λ r is the amplitude of sinusoidal air gap permeance due to rotor cogging effect.
The frequencies of PSH ft1− and ft1+ are respectively
f t 1 = ( Z 2 p ( 1 s ) 1 ) f 1 = ( Z 2 n 60 f 1 1 ) f 1
f t 1 + = ( Z 2 p ( 1 s ) + 1 ) f 1 = ( Z 2 n 60 f 1 + 1 ) f
in which f1 is the frequency of the fundamental magnetic field.
It can be found from Equations (17) and (18) that the frequencies of primary rotor slot harmonics change along with the driving frequency of the inverter and inherently contain the information of slip ratio s, which is adjusted during motor operation dynamics, thus making it difficult to extract the RSHs online.

2.2. The Effect of Skewed Slot Rotor on PSHs

The skewed slot of the rotor can effectively reduce the PSHs, thus improving the cogging torque as well as the electromechanical vibration and noise for an operating motor. In the magnetic field analysis, the effect of skewed slot can be equivalent to the decrease in winding coefficient and influence the distribution of PSHs, thus then changing the phase of electric potential harmonic for PSHs, so that their amplitudes will decrease. In addition, the skewed slot rotor also weakens F1.
Here we use a multisection method [23] to analyze the properties of PSHs with skewed rotor slots. In this method, the motor with a skewed slot rotor is segmented into several parts, and each motor part is regarded with no slot skewness. Then, the electromagnetic field of each motor part is calculated through the 2D finite element method [24], and the calculated results of each part are summed up and then averaged, so as to analyze the effect of rotor slot skewness.
The equivalent model of the multisection method mentioned above can be illustrated in Figure 2, where the motor with a skewed slot rotor is divided into K segments, with γ as the rotor slot slope. Thus, the rotor slot skewed angle between adjacent motor parts can be obtained as θ = γ/k.
By coping with the elements related to the vector magnetic potential in the finite element model of each motor segment, the finite element calculation model can be obtained. As an example, a nine-phase IM with its main parameters listed in Table 1 and the rotor slot skewed at a single stator pitch (5°) is divided into seven segments here. Each segment is regarded as a motor part without slot skewness and analyzed through FEA simulation. The 2D FEA simulation results of magnetic field distribution of a single segment are plotted in Figure 3.
After the 2D FEA simulation of each motor segment, the calculated results are further summed and averaged; the simulated amplitudes of PSH without and with the rotor slot skewness in stator induced phase voltage are shown in Figure 4.
It can be seen from Figure 4 that when the rotor slot is skewed by a stator tooth pitch, the amplitude of PSHs decreases significantly in stator induced phase voltage, compared with the non-skewed case. At the same time, the other harmonics are also reduced.

3. Adaptive Online Extraction of PSHs from Stator Currents

Because the rotor slot harmonic changes dynamically with the running state of the motor, and its amplitude is small and the signal−to−noise ratio is very low, it is actually difficult to use a filter with constant coefficients to realize the online extraction and analysis for these harmonics. Therefore, an adaptive filter is needed to process the stator current signal, so as to eliminate the unrequired harmonics but retain the PSH components.

3.1. The Online Parameter Updated Algorithm for the Adaptive Online Extraction Filter for PSHs

According to Equations (16) and (17), the center frequency of adaptive passband filter fpre0 and passband width fpw are designed as
f p r e 0 = Z 2 f 1 p ,   f p w = 4 f 1
Therefore, the lower and upper angular cut-off frequencies are respectively
ω 1 = 4 π ( f p r e 0 2 f 1 ) F s ,   ω 2 = 4 π ( f p r e 0 + 2 f 1 ) F s
A four-order Butterworth lowpass filter is then designed, with the transfer function of the Butterworth lowpass filter being
H L P ( s ) = 1 s 4 + 2.6131259 s 3 + 3.4142136 s 2 + 2.6131259 s + 1
Then the digital IIR bandpass filter can be obtained through bilinear transform [25],
H B P ( z ) = H L P ( s ) | s = D z 2 E z 1 + 1 1 z 2 = b 0 + b 1 z 1 + b 2 z 1 + b 3 z 3 a 0 a 1 z 1 a 2 z 1 a 3 z 3 a 4 z 4
where the coefficients of the filter can be obtained as
b 0 = 1 D 2 + 2 D + 1 ,   b 2 = 2 D 2 + 2 D + 1 ,   b 1 = b 3 = 0
a 0 = 1 ,   a 1 = ( 2 D 2 E + 2 D E ) D 2 + 2 D + 1 ,   a 2 = D 2 E 2 + 2 D 2 2 D 2 + 2 D + 1 , a 3 = 2 D E 2 D 2 E 2 D 2 + 2 D + 1 ,   a 4 = D 2 2 D + 1 D 2 + 2 D + 1
D = cot ( ω 2 ω 1 2 ) ,   E = 2 cos ( ω 2 + ω 1 2 ) cos ( ω 2 ω 1 2 )
Taking the direct type II transposed IIR filter as an example, the realization structure of the adaptive filter is shown in Figure 5, where K is the gain of the adaptive filter. The fundamental frequency track can be realized through phase-locked loop (PLL) or directly obtained from the inverter [26].

3.2. The Online Frequency Identification of PSHs and Verification

As the disturbances have been eliminated from the output of the adaptive filter mentioned above, the signal-to-noise ratio of the PSH signal will be significantly enhanced and can be further utilized for further analysis.
Additional processing is still needed to accurately identify the frequency of PSH. FFT and phase-locked loop (PLL) are the commonly used methods for calculating signal frequency. Since the output signal of the adaptive extraction filter still contains multiple frequency components and their frequencies change in real time, the PLL method is not applicable here for achieving stable output. Therefore, the FFT algorithm is used in this paper to analyze the spectrum of the output signal of the adaptive filter. The complete processing algorithm is diagramed in Figure 6, in which the adaptive extraction filter has already been illustrated in Figure 4 in detail. The anti-aliasing processing, decimation and signal reconstruction are the standard pre-processing of the signal before the FFT. The signal reconstruction reuses the sampling points discarded in the process of decimation, so as to increase the total signal length and improve the frequency identification accuracy. The frequency verification assures that the identified frequencies of PSHs are reasonable; for example, the difference between ft1− and ft1+ should be 2f1 under any working conditions.
The convergence rate is generally a critical concern of most adaptive filters due to the inherent signal feedback. However, the proposed online parameter update algorithm is a sequentially computed one, thus resulting in no convergence problem. Furthermore, the signal output transients due to filter parameter update will be averaged and diminished during the following FFT calculation and thus will not affect the RSH extraction results.
Another verification of the proposed method can be carried out by referring to actual motor speed. By combining Equations (17) and (18), the relationship between the frequency of PSHs and motor speed can be obtained as
f t 1 + f t 1 + = Z 2 n 15 p

4. Experimental Results

4.1. Experiment Platform Setup

A nine-phase induction motor platform, as shown in Figure 7, was established to verify the proposed adaptive online extraction and analysis method for PSHs. In this platform, a nine-phase squirrel-cage induction motor is driven by a multiphase inverter. The induction motor is coaxially connected to a DC generator, which is excited by a Siemens SINAMICS DCM regulator and serves for load adjustment. Speed and torque sensors are also installed with accuracies of 0.1 N·m and 1 r/min.
The stator current is measured through an LEM LA-25P current sensor and then A/D converted by ADS1203 and anti-aliasing processed by Xilinx xc4vlx25 FPGA at 10MHz rate before it is obtained by TI TMS320F28335 DSP. Differential signal measurement and Sigma-Delta ADC sample are utilized to improve the SNR. All the main parameters of the platform are listed in Table 1.

4.2. The PSH Extraction Results

Based on the above experimental device, the PSH characteristics of the motor were extracted under different rotational speeds and loads by utilizing the proposed adaptive online filter. In order to compare and analyze the time domain and frequency domain characteristics of the stator current and PSH of the nine-phase induction motor before and after adaptive filtering, The comparisons of the input and output of the adaptive PSH filter in time and frequency domain are given under different working conditions of speed 239–1396 rpm and torque 4.9–16.5 N·m, as shown in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15. It can be found that the proposed filter can effectively extract the PSHs from the stator current signal, which is rich in harmonics. All the above extraction and analysis results are summarized in Table 2.

4.3. Verification and Errors

The key metric of the proposed method is its universal applicability under different motor conditions. Compared with the available RSH extraction method that is designed for a certain motor speed and load condition, the proposed method is adjustable for different speed and load conditions with only one adaptive filter, which contributes to the scene working convenience. To better evaluate the accuracy of the proposed extraction method, the analytical model illustrated as Equation (25) is utilized with measured motor speeds at different working conditions. As the motor speeds can be easily measured online with enough accuracy through the resolvers or encoders, the proposed method can be more directly and effectively evaluated. The verification results are listed in Table 3, which shows that the errors of PSH extraction and frequency analysis are less than 1% for all tested working conditions. The error tends to increase with the decrease in speed. This is mainly because the PSHs are concentrated in the low-frequency domain and are close to each other at a low speed, which makes the extraction and analysis much more difficult. In fact, when the rotational speed continues to decrease, the bandwidth of the adaptive filter proposed is too narrow, which easily makes the output of the IIR filter oscillate or become unstable, making the extraction process invalid.

5. Conclusions

In this paper, an analytical model of PSHs was utilized to update the coefficients of an adaptive filter so that the PSHs of multiphase IMs could be extracted from a harmonic-rich motor current and analyzed online. It has been validated through experiments that the proposed method can be effectively utilized to cope with the changing frequency distribution of PSHs under different speed and load conditions.
The accuracy of RSH frequency identification has also been verified by referring to the measured rotor speed, with the relative error smaller than 1% in a wide speed range of 239–1396 rpm. In this case, the online-extracted properties of RSHs can be further utilized for motor cost and torque ripple evaluation and monitoring. However, further investigations are still needed for online extraction and analysis of RSHs at the ultra-low speed condition, since the passband of the IIR filter will be too narrow to be implemented.
Based on the proposed adaptive online extraction method, more RSH-related studies for multiphase induction motors, such as studies involving speed estimation, fault detection and state monitoring.

Author Contributions

Conceptualization, J.Y. and X.Z.; methodology, J.Y. and Y.Z.; software, H.S.; validation, J.Y. and H.S.; writing—original draft preparation, J.Y. and Y.Z.; writing—review and editing, J.Y. and Y.Z.; funding acquisition, J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

Natural Science Foundation of Shandong Province (Grant No. ZR2021QE118), Key Project of Innovative Teaching Laboratory of Qingdao University in 2020 (Grant No. CXSYZD202001) and National College Students Innovation and Entrepreneurship Training Program of China (Grant No. S202011065023).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study are available from the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Levi, E. Multiphase Electric Machines for Variable-Speed Applications. IEEE Trans. Ind. Electron. 2008, 55, 1893–1909. [Google Scholar] [CrossRef]
  2. Wang, H.; Wu, X.; Zheng, X.; Yuan, X. Virtual Voltage Vector based Model Predictive Control for a Nine-phase Open-end Winding PMSM with a Common DC bus. IEEE Trans. Ind. Electron. 2022, 69, 5386–5397. [Google Scholar] [CrossRef]
  3. You, R.; Lu, X. Voltage Unbalance Compensation in Distribution Feeders using Soft Open Points. J. Mod. Power Syst. Clean Energy 2022, 10, 1000–1008. [Google Scholar] [CrossRef]
  4. Gonzalez-Prieto, A.; Aciego, J.J.; Gonzalez-Prieto, I.; Duran, M.J. Automatic Fault-Tolerant Control of Multiphase Induction Machines: A Game Changer. Electronics 2020, 9, 938. [Google Scholar] [CrossRef]
  5. Kellner, J.; Kaščák, S.; Ferková, Ž. Investigation of the Properties of a Five-Phase Induction Motor in the Introduction of New Fault-Tolerant Control. Appl. Sci. 2022, 12, 2249. [Google Scholar] [CrossRef]
  6. Li, X.; Xue, Z.; Zhang, L.; Hua, W. A low-complexity three-vector-based model predictive torque control for SPMSM. IEEE Trans. Power Electron. 2021, 36, 13002–13012. [Google Scholar] [CrossRef]
  7. Joksimović, G.; Levi, E.; Kajević, A.; Mezzarobba, M.; Tessarolo, A. Optimal Selection of Rotor Bar Number for Minimizing Torque and Current Pulsations Due to Rotor Slot Harmonics in Three-Phase Cage Induction Motors. IEEE Access 2020, 8, 228572–228585. [Google Scholar] [CrossRef]
  8. Frauman, P.; Burakov, A.; Arkkio, A. Effects of the Slot Harmonics on the Unbalanced Magnetic Pull in an Induction Motor with an Eccentric Rotor. IEEE Trans. Magn. 2007, 43, 3441–3444. [Google Scholar] [CrossRef]
  9. Bień, A.; Drabek, T.; Kara, D.; Kołacz, T. Reduction in the Cogging Torques in the DCEFSM Motor by Changing the Geometry of the Rotor Teeth. Energies 2022, 15, 2455. [Google Scholar] [CrossRef]
  10. Guo, L.; Xu, J.; Wu, S.; Xie, X.; Wang, H. Analysis and Design of Dual Three-Phase Fractional-Slot Permanent Magnet Motor with Low Space Harmonic. IEEE Trans. Magn. 2022, 58, 8100112. [Google Scholar] [CrossRef]
  11. Zhou, G.; Shen, J. Current Harmonics in Induction Machine with Closed-Slot Rotor. IEEE Trans. Ind. Appl. 2017, 53, 134–142. [Google Scholar] [CrossRef]
  12. Asad, B.; Vaimann, T.; Belahcen, A.; Kallaste, A.; Rassolkin, A.; Van Khang, H.; Ghahfarokhi, P.S.; Naseer, M.U.; Iqbal, M.N. The Modeling and Investigation of Slot Skews and Supply Imbalance on the Development of Principal Slotting Harmonics in Squirrel Cage Induction Machines. IEEE Access 2021, 9, 165932–165946. [Google Scholar] [CrossRef]
  13. Chen, Q.; Xu, G.; Liu, G.; Zhai, F.; Eduku, S. Reduction of Torque Ripple Caused by Slot Harmonics in FSCW Spoke-Type FPM Motors by Assisted Poles. IEEE Trans. Ind. Electron. 2020, 67, 9613–9622. [Google Scholar] [CrossRef]
  14. Zhang, D.; Liu, T.; Zhao, H.; Wu, T. An Analytical Iron Loss Calculation Model of Inverter-Fed Induction Motors Considering Supply and Slot Harmonics. IEEE Trans. Ind. Electron. 2019, 66, 9194–9204. [Google Scholar] [CrossRef]
  15. Diarra, M.N.; Yao, Y.; Li, Z.; Niasse, M.; Li, Y.; Zhao, H. In-Situ Efficiency Estimation of Induction Motors Based on Quantum Particle Swarm Optimization-Trust Region Algorithm (QPSO-TRA). Energies 2022, 15, 4905. [Google Scholar] [CrossRef]
  16. Benouzza, N.; Boudinar, A.H.; Bendiabdellah, A.; Khodja, M. Slot harmonic frequency detection as a technique to improve stator current spectrum approach for broken rotor bars fault diagnosis. In Proceedings of the 2015 Intl Aegean Conference on Electrical Machines & Power Electronics (ACEMP), 2015 Intl Conference on Optimization of Electrical & Electronic Equipment (OPTIM) & 2015 Intl Symposium on Advanced Electromechanical Motion Systems (ELECTROMOTION), Side, Turkey, 2–4 September 2015; pp. 118–122. [Google Scholar] [CrossRef]
  17. Gao, Z.; Habetler, T.G.; Harley, R.G.; Colby, R.S. A Sensorless Rotor Temperature Estimator for Induction Machines Based on a Current Harmonic Spectral Estimation Scheme. IEEE Trans. Ind. Electron. 2008, 55, 407–416. [Google Scholar] [CrossRef]
  18. Ray, S.; Dey, D. Development of a Comprehensive Analytical Model of Induction Motor under Stator Intern Turn Faults Incorporating Rotor Slot Harmonics. IEEE Trans. Ind. Electron. 2022. [Google Scholar] [CrossRef]
  19. Degner, M.W.; Lorenz, R.D. Position estimation in induction machines utilizing rotor bar slot harmonics and carrier-frequency signal injection. IEEE Trans. Ind. Appl. 2000, 36, 736–742. [Google Scholar] [CrossRef]
  20. Duan, X.; Zhang, X.; Tang, Y.; Hao, M. Cogging Torque Reduction in PMSM in Wide Temperature Range by Response Surface Methodology. Symmetry 2021, 13, 1877. [Google Scholar] [CrossRef]
  21. Park, J.-C.; Park, S.-H.; Kim, J.-H.; Lee, S.-G.; Lee, G.-H.; Lim, M.-S. Diagnosis and Robust Design Optimization of SPMSM Considering Back EMF and Cogging Torque due to Static Eccentricity. Energies 2021, 14, 2900. [Google Scholar] [CrossRef]
  22. Lee, C.-S.; Kim, H.-J. Harmonic Order Analysis of Cogging Torque for Interior Permanent Magnet Synchronous Motor Considering Manufacturing Disturbances. Energies 2022, 15, 2428. [Google Scholar] [CrossRef]
  23. Brescia, E.; Costantino, D.; Massenio, P.R.; Monopoli, V.G.; Cupertino, F.; Cascella, G.L. A Design Method for the Cogging Torque Minimization of Permanent Magnet Machines with a Segmented Stator Core Based on ANN Surrogate Models. Energies 2021, 14, 1880. [Google Scholar] [CrossRef]
  24. Palangar, M.F.; Mahmoudi, A.; Kahourzade, S.; Soong, W.L. Electromagnetic and Thermal Analysis of a Line-Start Permanent-Magnet Synchronous Motor. In Proceedings of the 2020 IEEE Energy Conversion Congress and Exposition (ECCE), Detroit, MI, USA, 11–15 October 2020; pp. 502–508. [Google Scholar] [CrossRef]
  25. Yu, J.; Zhu, C.; Wu, X. The Analysis of Displacement Estimation Errors of Demodulation Type Self-sensing Active Magnetic Bearings. Proc. Chin. Soc. Electron. Eng. 2019, 10, 5876–5884. [Google Scholar] [CrossRef]
  26. Wang, H.; Zheng, X.; Yuan, X.; Wu, X. Low-Complexity Model-Predictive Control for a Nine-Phase Open-End Winding PMSM with Dead-Time Compensation. IEEE Trans. Ind. Electron. 2022, 37, 8895–8908. [Google Scholar] [CrossRef]
Figure 1. The cogging effect of motor slotting on air gap magnetic field.
Figure 1. The cogging effect of motor slotting on air gap magnetic field.
Energies 15 06643 g001
Figure 2. The equivalent multisegmented model for a skewed rotor slot.
Figure 2. The equivalent multisegmented model for a skewed rotor slot.
Energies 15 06643 g002
Figure 3. The 2D FEA simulation results of a segmented nine-phase IM: (a) the distribution of A (Wb/m); (b) the distribution of B (T).
Figure 3. The 2D FEA simulation results of a segmented nine-phase IM: (a) the distribution of A (Wb/m); (b) the distribution of B (T).
Energies 15 06643 g003
Figure 4. The effect of skewed rotor slot on the induced phase voltage harmonics of a nine-phase IM (p.u.).
Figure 4. The effect of skewed rotor slot on the induced phase voltage harmonics of a nine-phase IM (p.u.).
Energies 15 06643 g004
Figure 5. The diagram of the adaptive filter for PSH extraction.
Figure 5. The diagram of the adaptive filter for PSH extraction.
Energies 15 06643 g005
Figure 6. The online extraction and frequency identification of PSH.
Figure 6. The online extraction and frequency identification of PSH.
Energies 15 06643 g006
Figure 7. The nine-phase induction motor platform.
Figure 7. The nine-phase induction motor platform.
Energies 15 06643 g007
Figure 8. The comparison of the input and output of adaptive PSH filter in time and frequency domain (239 r/min, TL = 5.2 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 8. The comparison of the input and output of adaptive PSH filter in time and frequency domain (239 r/min, TL = 5.2 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g008
Figure 9. The comparison of the input and output of adaptive PSH filter in time and frequency domain (342 r/min, TL = 4.9 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 9. The comparison of the input and output of adaptive PSH filter in time and frequency domain (342 r/min, TL = 4.9 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g009
Figure 10. The comparison of the input and output of adaptive PSH filter in time and frequency domain (440 r/min, TL = 6.0 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 10. The comparison of the input and output of adaptive PSH filter in time and frequency domain (440 r/min, TL = 6.0 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g010
Figure 11. The comparison of the input and output of adaptive PSH filter in time and frequency domain (538 r/min, TL = 7.0 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 11. The comparison of the input and output of adaptive PSH filter in time and frequency domain (538 r/min, TL = 7.0 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g011
Figure 12. The comparison of the input and output of adaptive PSH filter in time and frequency domain (636 r/min, TL = 8.3 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 12. The comparison of the input and output of adaptive PSH filter in time and frequency domain (636 r/min, TL = 8.3 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g012
Figure 13. The comparison of the input and output of adaptive PSH filter in time and frequency domain (832 r/min, TL = 10.5 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 13. The comparison of the input and output of adaptive PSH filter in time and frequency domain (832 r/min, TL = 10.5 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g013
Figure 14. The comparison of the input and output of adaptive PSH filter in time and frequency domain (1028 r/min, TL = 12.7 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 14. The comparison of the input and output of adaptive PSH filter in time and frequency domain (1028 r/min, TL = 12.7 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g014
Figure 15. The comparison of the input and output of adaptive PSH filter in time and frequency domain (1396 r/min, TL = 16.5 N·m): (a) time domain waveform; (b) frequency distribution.
Figure 15. The comparison of the input and output of adaptive PSH filter in time and frequency domain (1396 r/min, TL = 16.5 N·m): (a) time domain waveform; (b) frequency distribution.
Energies 15 06643 g015
Table 1. Main parameters of nine-phase induction motor platform.
Table 1. Main parameters of nine-phase induction motor platform.
SymbolParameterValue
PNRated motor power8 kW
mPhase number9
fNRated power supply frequency50 Hz
PPole pairs2
nNRated speed1477 r/min
UNRated phase voltage85 V
INRated phase current12.87 A
Z1Stator slot number72
Z2Rotor slot number54
βPitch ratio5/6
Table 2. The frequency distribution of PSHs under different speed and load conditions.
Table 2. The frequency distribution of PSHs under different speed and load conditions.
Speed n (rpm)Torque (N·m)f1 (Hz)ft1− (Hz)ft1+ (Hz)
2395.28.3208.4225.1
3424.911.7296.8320.2
4406.015.0381.7411.7
5387.018.3466.6503.3
6368.321.6551.2594.6
83210.528.3721.1777.8
102812.735.0891.0961.0
139616.548.31208.01304.6
Table 3. Verification and errors.
Table 3. Verification and errors.
Speed n (rpm)Torque (N·m)Z2n/15p (Hz)ft1− + ft1− (Hz)Relative Error(%)
2395.2430.2433.50.77
3424.9615.6617.00.23
4406.0792.0793.40.18
5387.0968.4969.90.15
6368.31144.81145.80.09
83210.51497.61498.90.09
102812.71850.41852.00.09
139616.52512.82512.60.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, J.; Zhang, Y.; Shen, H.; Zheng, X. Adaptive Online Extraction Method of Slot Harmonics for Multiphase Induction Motor. Energies 2022, 15, 6643. https://doi.org/10.3390/en15186643

AMA Style

Yu J, Zhang Y, Shen H, Zheng X. Adaptive Online Extraction Method of Slot Harmonics for Multiphase Induction Motor. Energies. 2022; 15(18):6643. https://doi.org/10.3390/en15186643

Chicago/Turabian Style

Yu, Jie, Youjun Zhang, Hongyuan Shen, and Xiaoqin Zheng. 2022. "Adaptive Online Extraction Method of Slot Harmonics for Multiphase Induction Motor" Energies 15, no. 18: 6643. https://doi.org/10.3390/en15186643

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