Research Article  Open Access
Stator Fault Detection in Induction Motors by Autoregressive Modeling
Abstract
This study introduces a novel methodology for early detection of stator short circuit faults in induction motors by using autoregressive (AR) model. The proposed algorithm is based on instantaneous space phasor (ISP) module of stator currents, which are mapped to statorfixed reference frame; then, the module is obtained, and the coefficients of the AR model for such module are estimated and evaluated by order selection criterion, which is used as fault signature. For comparative purposes, a spectral analysis of the ISP module by Discrete Fourier Transform (DFT) is performed; a comparison of both methodologies is obtained. To demonstrate the suitability of the proposed methodology for detecting and quantifying incipient short circuit stator faults, an induction motor was altered to induce differentdegree fault scenarios during experimentation.
1. Introduction
The fault diagnosis for electrical machines has become far more critical as the population of electric machines has greatly increased in recent years. The total number of operating electrical machines in the world was around 16.1 billion in 2011, with a growth rate of about 50% in the last five years [1].
Threephase induction motors represent the principal source of movement in the electrical and processing industry [2]. This type of motors has an exclusive position in the energy conversion and they are responsible of almost 90% of the electric energy consumed by electric machines. Threephase induction motors correspond to nearly 60% of all electric machines [3].
The success of threephase induction is due to the low cost, robustness, and high performance in variable speed tasks; this has been achieved through the development of new control laws and more versatile semiconductor devices [4]. However, most control algorithms may not be reliable when a fault condition occurs [1, 2].
Stator short circuit fault diagnosis in threephase induction motors represents a significant percentage of the electrical machine defects [5]. In this category there are faults in the statorwinding such as short circuits among turns and the magnetic circuit. In the first case, the internal asymmetry will cause the circulation of high currents in the portion of the statorwinding affected by the fault, this contributes to the degradation of other portions of the winding, and the remaining time between the onset of fault and the failure depends on various factors such as the initial number of shorted turns, winding configuration, the power and voltage achieved, and environmental conditions among others [6]. Therefore early detection is critical.
Unfortunately the induction motor may fail due to other mechanisms, such as bearings [7], the rotor due cracks in shortcircuiting rings or broken bars [8], and eccentricity in the air gap [9]. In [10, 11], some of the techniques commonly used for this purpose are shown.
A comparative fault distribution informed by [12] is shown in Table 1; this comparison was carried out in several surveys of faults of large induction motors, conducted by the Institute of Electrical and Electronics Engineers (IEEE) [13–15], Electric Power Research Institute (EPRI) [16], and Allianz [17]. Differences in distributions are because survey conducted by IEEE and EPRI focuses on mediumsized induction machines, while Allianz survey focuses more on medium to highvoltage large induction machines.

Spectral analysis has been widely used for diagnosis of electric and mechanical faults [18, 19]. Autoregressive (AR) models characterize spectral and energy properties of a system by representing them as time series, which consider each sample in the time domain as a linear combination of previous sampled values in the presence of white noise with a Gaussian distribution [20]. Therefore, the coefficients of a univariate AR model are the autocorrelations of the time series, which represent a dependence measure of the current value with respect to previous values, and provide information of the system dynamics in time. When a fault occurs in a system, its dynamics changes; therefore, the AR model coefficients change too. This suggests that a parameter obtained from these coefficients can be used as a fault diagnosis tool for induction machines.
The use of AR model coefficients has been reported in some works; for instance, in [21, 22] bearing faults are detected, and in [23, 24] faults in gearbox coupling systems are detected. In these works, coefficients are used to estimate the power spectral density (PSD) from the transfer function of the AR model. The harmonic analysis allows detecting those faults that cause atypical frequencies and parasitic oscillations [25].
This paper presents a novel stator fault detection technique based on AR model of the ISP module for the stator reference frame from the Concordia transformation, which is evaluated via the Akaike Final Prediction Error (FPE) criterion to obtain an indicator of the induction motor operational condition.
2. Fault Detection Scheme
2.1. AR Model
The proposed electric stator fault diagnosis scheme consists in acquiring samples of the induction motor three phase stator currents , , and :The three currents are mapped into the statorfixed reference frame using Concordia transformation: Then, ISP module [26] is obtained for each sample of and by
From the ISP module the coefficients of a order AR model are estimated. The AR model is a finite impulse response filter (FIR) that can be defined as follows:where contains the coefficients of the AR model and is the value of the ISP module delayed samples. Finally, represents uncorrelated random noise.
AR model coefficients can be estimated by diverse methods like the solution of YuleWalker equations [27], the Burg method [28], or the covariance method [29]. In this work the coefficients are calculated by least squares approximation algorithm [30], where the order AR model (4) is estimated in the form of a regression model expressed aswith a coefficients vectorand with predictors In terms of momentsThe least squared estimates of coefficients are estimated byThe residual varianceis an estimate of the noise covariance and can be expressed as The residual variance is proportional to Schur complement of the matrix which is the moment matrix belonging to the data matrixwhere the data matrix .
The least squares estimation (9) is obtained from a QR factorization of through the Householder algorithm:where and are orthogonal upper triangular matrices, and is given bywhere , , , and .
The QR factorization of the data matrix (14) leads to the Cholesky factorization of the moment matrix [31]The least squares estimation of the residual variance and the parameter vector are computed from the Cholesky factorization by
Regardless of the used method, it is necessary to define the model order that best fits the signal of interest (3). The chosen order selection criterion was the FPE which was evaluated for a order AR model by
Finally, the fluctuations in FPE values for different cases with and without fault are compared.
2.2. Frequency Scheme
The proposed method was compared with a frequency domain analysis using the Discrete Fourier Transform (DFT) algorithm. The stator fault spectral signature of the ISP module occurs in the integer multiples of twice the voltage frequency applied to the motor :where [32]. Since components of the fault spectral signature decrease while increases, just the main component () is considered; hence
3. Experimental Setup
The proposed methodology was validated on an experimental platform. The implemented test bed for performing the experimental analysis is depicted in Figure 1. It consists in a 3HP, threephase induction motor altered to introduce different electric stator fault scenarios; a regeneration system based on a direct current (DC) motor was used for adding load torque; a full bridge threephase inverter was used for varying the motor rotating speed; a Halleffect based data acquisition system was used for acquiring the stator currents; the induction motor was instrumented with an optical encoder to acquire the rotor angular position and speed shaft; and a central processing unit (CPU) based on TI TMS320F28335 32bit floating point digital signal processor (DSP) was used for performing the necessary control and measurement functions.
The CPU received a group of commands from a personal computer (PC) to acquire and control the instrumentation of a power electronic system on the test bed.
The stator currents were acquired and stored in the memory of an IRAM DSP, taking 4096sample frames from each phase. The acquired signals were sent to a PC via high speed USB. A computerprogram based on Delphi language was used to perform the fault diagnosis analysis.
In order to insert different types of stator faults, like interturn short circuit or open circuit in one phase, the induction machine was rewinded. The schematic diagram of the induction motor is shown in Figure 2.
The electric structure for the statorwinding comprises 6 groups of coils, 2 groups per phase (A, A′, B, B′, and C, C′) distributed in 36 slots, and three concentric coils per group of 45 turns each, with three 22AGW wires for a total of 7290 turns.
The lowvoltage operational condition for the induction motor was achieved by connecting in parallel the double star windings and supplying the phase voltages to terminals 1 and 7, to 2 and 8, and to 3 and 9 and joining terminals 4, 5, and 6 to conform the neutral star , as specified by the manufacturer.
The fault insertion in the statorwinding is accomplished by connecting or disconnecting some segments of phase C winding to phase voltage or neutral , using tap derivations . This is achieved by an insulated electric circuit containing a group of relays that are managed by the DSPbased control system to select the fault severity in the motor. Possible connections are depicted with dotted lines in Figure 3.
The characterization of the induction motor parameters was performed at line frequency (60 Hz) via the threevoltmeter method, which consists in inserting a pattern resistance in each phase of the motor and measuring stator voltages , , and and currents as shown in Figure 4; then, the impedance , inductive reactance , resistance , and inductance for each phase are calculated. The results are shown in Table 2.

4. Experimental Results
A comparative analysis was performed between the FPE of the AR model coefficients and fault spectral signature for ISP module using DFT. Twenty trials were performed on each one of the treated conditions. Three different states were considered: healthy case and two distinctdegree stator fault conditions as shown in Table 3. The healthy condition is obtained by to terminal 3 and to terminal 6 (see Figure 3). In the first faulty condition (three turns), is connected to terminal. In the second faulty condition (nine turns) is connected to terminal. Both faulty cases correspond to incipient short circuit stator faults.

The defect level for each faulty condition is indicated through a change in the percentage ratio (Δ%) of their corresponding coefficient magnitudes with respect to the healthy case.
Experimental data were obtained from the ISP module. For each acquired stator current, 4096 samples were obtained at the sampling frequency Hz. A load torque of 0.4 Nm was applied for all experiments at different speeds: 1200 rpm ( Hz), 1500 rpm ( Hz), and 1800 rpm ( Hz).
4.1. AR Model Approach
The mean FPE values of the ISP module from the three treated cases at rotating speed of 1800 rpm ( Hz), evaluated for different orders (), are shown in Figure 5. It can be observed that the mean FPE value decreases for the three cases when the AR model order increases. The gap among the values for orders () is minimum; however, the gap increases from order and it remains as the order goes up.
The mean of FPE value and their percentage rate of change Δ% among different fault degrees with respect to the healthy case at three different speeds for a 35order AR model are summarized in Table 4.

4.2. Spectra Analysis
Figure 6 shows the stator fault spectral signature corresponding to the three considered conditions reported in Table 5, at the rotating speed of 1800 rpm ( Hz). From this figure, it can be observed that the frequency component at 120 Hz increases its magnitude with respect to the fault level. Table 5 summarizes Δ% average percentage rate of the fault magnitude of characteristic frequency components for the different fault scenarios in comparison with the healthy case at three different speeds.

4.3. Discussion of Results
From the obtained results, the proposed methodology does not require a powertwo or adjusting the sampling frequency depending on the motor rotating speed. On the other hand, the common spectral analysis technique for detecting fault signatures does require to be a multiple of [28]; furthermore, in the case of not being a power of two, zero padding must have to be done, and the technique will fail.
For a fixed order AR model with a defined sampling period, for instance, and , the FPE computation and estimation of the AR model coefficients will require ( = 5,061,700 operations [30] for solving the linear equation system with real numbers. The proposed approach was executed on an Intel 2.5 GHz processor i74710HQ, and it took around 0.00401 seconds, whereas the DFT technique requires = 16,777,216 operations and takes 0.0129 seconds, under exactly the same conditions, involving the product of exponential functions with complex numbers.
The downward trend in the mean FPE value as the model order increases indicates that the higher the order of the AR model, the more accurate the fault estimation; however, the computational burden will also increase. Nevertheless, for a model order (see Figure 5) a tendency to keep each operational condition separated is clear. Hence, for diagnosis purposes, this indicates that it is not necessary to select a very highorder AR model.
On the other hand, it was possible to detect the presence of fault conditions even with low order models , but it was not possible to quantify the magnitude of the fault, whereas, for a higherorder model and given the scope of this research that only characterized the electrical stator parameters for Hz, it was shown that there is approximated proportionality between the fault level and the inductance value.
For instance, in the first faulty case the increment for the three the FPE value 8.06% with respect to the healthy condition the inductance value and the 7.6% obtained by three voltmeter method. While for the second case, it was of 11.5% against 8.6%, so that in a future work it will be necessary to characterize the induction motor at different frequencies and in order to verify whether it is possible to quantify the level of damage, without evaluating the transfer function of the AR model in the frequency domain as it is commonly used in literature.
To demonstrate the effectiveness of the proposed methodology, a 3sigma rule was applied on 20 experiments for each of the different treated conditions. The obtained results are depicted in Figure 7; from this figure and due to the fact that the Gaussian distributions do not overlap in their corresponding 3sigma rule interval for (healthy case and case 1), a confidence level of more than 99.7% is endured on the detection and 95% in the identification of stator faults.
5. Conclusions
Electric motors are responsible for around 90% of the electricmachine energy consumption. Stator faults are very important because of the effects on the performance of the induction machine; hence, several studies are available for induction motor fault detection. In this regard, this work proposed a novel methodology based on the use of FPE for stator fault detection in induction motors. The introduced approach was experimentally validated considering three different operational conditions: healthy case and two differentdegree stator faults. From the obtained results the usefulness and efficiency of the proposed method for fault detection and quantification of incipient short circuit stator faults, considering low AR model orders for different faulty conditions, were demonstrated considering a 3sigma rule.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was partially supported by the Doctorado en Ciencias de la Ingeniería of the Instituto Tecnologico de Aguascalientes and the PROMEP under Grants PROMEP/ITA/103.5/13/8890 and PROMEP/UASLP/CARepPrel15 and Universidad Autónoma de San Luis Potosí under Grant C14FAI0442.42. Mario A. GonzalezGarcia is a Research Fellow of the National Council of Science and Technology (CONACyT), Mexico.
References
 H. A. Toliyat, S. Nandi, S. Choi, and H. MeshgingKelk, Electric Machines, CRC Press, Boca Raton, Fla, USA, 2013.
 H. Henao, G.A. Capolino, M. FernandezCabanas et al., “Trends in fault diagnosis for electrical machines: a review of diagnostic techniques,” IEEE Industrial Electronics Magazine, vol. 8, no. 2, pp. 31–42, 2014. View at: Publisher Site  Google Scholar
 S. Nandi, H. A. Toliyat, and X. Li, “Condition monitoring and fault diagnosis of electrical motors—a review,” IEEE Transactions on Energy Conversion, vol. 20, no. 4, pp. 719–729, 2005. View at: Publisher Site  Google Scholar
 F. Betin, G.A. Capolino, D. Casadei et al., “Trends in electrical machines control: samples for classical, sensorless, and faulttolerant techniques,” IEEE Industrial Electronics Magazine, vol. 8, no. 2, pp. 43–55, 2014. View at: Publisher Site  Google Scholar
 C. H. De Angelo, G. R. Bossio, S. J. Giaccone, M. I. Valla, J. A. Solsona, and G. O. García, “Online modelbased statorfault detection and identification in induction motors,” IEEE Transactions on Industrial Electronics, vol. 56, no. 11, pp. 4671–4680, 2009. View at: Publisher Site  Google Scholar
 S. M. A. Cruz and A. J. M. Cardoso, “Multiple reference frames theory: a new method for the diagnosis of stator faults in threephase induction motors,” IEEE Transactions on Energy Conversion, vol. 20, no. 3, pp. 611–619, 2005. View at: Publisher Site  Google Scholar
 Z. Yang, “Automatic condition monitoring of industrial rollingelement bearings using motor's vibration and current analysis,” Shock and Vibration, vol. 2015, Article ID 486159, 12 pages, 2015. View at: Publisher Site  Google Scholar
 H.Y. Zhu, J.T. Hu, L. Gao, and H. Huang, “Practical aspects of broken rotor bars detection in PWM voltagesourceinverterfed squirrelcage induction motors,” Journal of Applied Mathematics, vol. 2013, Article ID 128368, 11 pages, 2013. View at: Publisher Site  Google Scholar
 R. Shnibha, A. Albarbar, A. Abouhnik, and G. Ibrahim, “A more reliable method for monitoring the condition of threephase induction motors based on their vibrations,” ISRN Mechanical Engineering, vol. 2012, Article ID 230314, 9 pages, 2012. View at: Publisher Site  Google Scholar
 A. H. Bonnett and C. Yung, “Increased efficiency versus increased reliability,” IEEE Industry Applications Magazine, vol. 14, no. 1, pp. 29–36, 2008. View at: Publisher Site  Google Scholar
 M. HernandezVargas, E. CabalYepez, and A. GarciaPerez, “Realtime SVDbased detection of multiple combined faults in induction motors,” Computers & Electrical Engineering, vol. 40, no. 7, pp. 2193–2203, 2014. View at: Publisher Site  Google Scholar
 P. Zhang, Y. Du, T. G. Habetler, and B. Lu, “A survey of condition monitoring and protection methods for mediumvoltage induction motors,” IEEE Transactions on Industry Applications, vol. 47, no. 1, pp. 34–46, 2011. View at: Publisher Site  Google Scholar
 R. N. Bell, D. W. McWilliams, P. O'Donnell, C. Singh, and S. J. Wells, “Report of large motor reliability survey of industrial and commercial installations, part I,” IEEE Transactions on Industry Applications, vol. IA21, no. 4, pp. 853–864, 1985. View at: Publisher Site  Google Scholar
 “Report of large motor reliability survey of industrial and commercial installations, Part II,” IEEE Transactions on Industry Applications, vol. 21, no. 4, pp. 865–872, 1985. View at: Publisher Site  Google Scholar
 “Report of large motor reliability survey of industrial and commercial installations: part 3,” IEEE Transactions on Industry Applications, vol. 21, no. 4, pp. 865–872, 1987. View at: Publisher Site  Google Scholar
 P. F. Albrecht, J. C. Appiarius, R. M. McCoy, E. L. Owen, and D. K. Sharma, “Assessment of the reliability of motors in utility applications—updated,” IEEE Transactions on Energy Conversion, vol. 1, no. 1, pp. 39–46, 1986. View at: Publisher Site  Google Scholar
 H. O. Seinsch, “Monitoring und diagnose elektrischer maschinen und antriebe,” in Proceedings VDE Workshop, Allianz Schadensstatistik an HS Motoren, 1996–1999, VDE, 2001. View at: Google Scholar
 L.L. Jiang, H.K. Yin, X.J. Li, and S.W. Tang, “Fault diagnosis of rotating machinery based on multisensor information fusion using SVM and timedomain features,” Shock and Vibration, vol. 2014, Article ID 418178, 8 pages, 2014. View at: Publisher Site  Google Scholar
 P. A. DelgadoArredondo, A. GarciaPerez, D. MorinigoSotelo et al., “Comparative study of timefrequency decomposition techniques for fault detection in induction motors using vibration analysis during startup transient,” Shock and Vibration, vol. 2015, Article ID 708034, 14 pages, 2015. View at: Publisher Site  Google Scholar
 J. Morton, P. A. Torrione, and L. M. Collins, “Variational Bayesian learning for mixture autoregressive models with uncertainorder,” IEEE Transactions on Signal Processing, vol. 59, no. 6, pp. 2614–2627, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 G. Wang, Z. Luo, X. Qin, Y. Leng, and T. Wang, “Fault identification and classification of rolling element bearing based on timevarying autoregressive spectrum,” Mechanical Systems and Signal Processing, vol. 22, no. 4, pp. 934–947, 2008. View at: Publisher Site  Google Scholar
 C. Liu, X. Zhou, S. Yang, W. Liang, and Q. Miao, “Cooling fan bearing diagnosis based on AR& MED,” in Proceedings of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering (ICQR2MSE '12), pp. 622–626, Chengdu, China, June 2012. View at: Publisher Site  Google Scholar
 G.H. Zhou, C.C. Zuo, J.Z. Wang, and S.X. Liu, “Gearbox fault diagnosis based on waveletAR model,” in Proceedings of the 6th International Conference on Machine Learning and Cybernetics (ICMLC '07), vol. 2, pp. 1061–1065, Hong Kong, August 2007. View at: Publisher Site  Google Scholar
 N. K. Nikhar, S. S. Patankar, and J. V. Kulkarni, “Gear tooth fault detection by autoregressive modelling,” in Proceedings of the International Conference on Computing, Communications and Networking Technologies (ICCCNT '13), pp. 1–6, Tiruchengode, India, July 2013. View at: Google Scholar
 F. J. VillalobosPiña and R. AlvarezSalas, “Robust algorithm for electric fault diagnosis in the three phase induction machine based on spectral and wavelet tools,” Revista Iberoamericana de Automatica e Informatica Industrial, vol. 12, no. 3, pp. 292–303, 2015. View at: Publisher Site  Google Scholar
 D. L. Milanez and A. E. Emanuel, “The instantaneousspacephasor: a powerful diagnosis tool,” IEEE Transactions on Instrumentation and Measurement, vol. 52, no. 1, pp. 143–148, 2003. View at: Publisher Site  Google Scholar
 X. Lin, Y. Xie, D. Zhao, and S. Xu, “Estimation of observer parameters for dynamic positioning ships,” Mathematical Problems in Engineering, vol. 2013, Article ID 173603, 7 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 J. G. Proakis and D. G. Manolakis, Digital Signal Processing, Pearson Prentice Hall, New Jersey, NJ, USA, 2007.
 M. H. Hayes, Statistical Digital Signal Processing and Modeling, John Wiley & Sons, New York, NY, USA, 1996.
 A. Neumaier and T. Schneider, “Estimation of parameters and eigenmodes of multivariate autoregressive models,” ACM Transactions on Mathematical Software, vol. 27, no. 1, pp. 27–57, 2001. View at: Publisher Site  Google Scholar
 H.C. Wu, S. Y. Chang, T. LeNgoc, and Y. Wu, “Efficient rankadaptive leastsquare estimation and multipleparameter linear regression using novel dyadically recursive hermitian matrix inversion,” International Journal of Antennas and Propagation, vol. 2012, Article ID 891932, 10 pages, 2012. View at: Publisher Site  Google Scholar
 M. Drif and A. J. M. Cardoso, “Airgapeccentricity fault diagnosis, in threephase induction motors, by the complex apparent power signature analysis,” IEEE Transactions on Industrial Electronics, vol. 55, no. 3, pp. 1404–1410, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Francisco M. GarciaGuevara et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.