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 α-β stator-fixed 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 different-degree 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].
Three-phase 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. Three-phase induction motors correspond to nearly 60% of all electric machines [3].
The success of three-phase 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 three-phase induction motors represents a significant percentage of the electrical machine defects [5]. In this category there are faults in the stator-winding 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 stator-winding 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 short-circuiting 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 medium-sized induction machines, while Allianz survey focuses more on medium- to high-voltage large induction machines.
Comparison of IEEE, EPRI, and Allianz surveys.
Component
Fault percent distribution
IEEE
EPRI
Allianz
Bearings
44
41
13
Stator
26
36
63
Rotor
8
9
13
Others
22
14
8
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 Scheme2.1. AR Model
The proposed electric stator fault diagnosis scheme consists in acquiring N samples of the induction motor three phase stator currents ia, ib, and ic:(1)iaibic=ia1ia2⋯iaNib1ib2⋯ibNic1ic2⋯icN.The three currents are mapped into the α-β stator-fixed reference frame using Concordia transformation: (2)iαiβio=121cos-2π3cos2π30-sin-2π3-sin2π3121212iaibic.Then, ISP module [26] is obtained for each sample of iα and iβ by(3)ISP=iαβ=iα2+iβ2.
From the ISP module the coefficients of a p-order AR model are estimated. The AR model is a finite impulse response filter (FIR) that can be defined as follows:(4)ISPn=∑k=1pAkISPn-k+εn,n=1,…,N,where Ak contains the coefficients of the AR model and ISPn-k is the value of the ISP module delayed k samples. Finally, εn represents uncorrelated random noise.
AR model coefficients can be estimated by diverse methods like the solution of Yule-Walker 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 p-order AR model (4) is estimated in the form of a regression model expressed as(5)ISPn=Apun+εnwith εn=noiseCwith a coefficients vector(6)A1,…,Apand with predictors (7)un=1ISPn-1⋯ISPn-pT.In terms of moments(8)U=∑n=1NununT,ISP=∑n=1NISPnISPnT,W=∑n=1NISPnunT.The least squared estimates of coefficients A^ are estimated by(9)A^=WU-1.The residual variance(10)C^=1N-p∑n=1Nε^nε^nTwith ε^n=ISPn-A^punis an estimate of the noise covariance C and can be expressed as (11)C^=1N-p+1ISP-WU-1WT.The residual variance is proportional to Schur complement of the matrix (12)Γ=UWTWISP=∑n=1NunISPnunTISPnTwhich is the moment matrix Γ=KTK belonging to the data matrix(13)K=u1TISP1T⋮⋮uNTISPNT,where the data matrix K∈RN-p×p+1.
The least squares estimation (9) is obtained from a QR factorization of K through the Householder algorithm:(14)K=QR,where Q and R are orthogonal upper triangular matrices, and R is given by(15)R=R11R120R22,where R∈RN+1×p, R11∈Rp×p, R12∈R1×p, and R22∈R.
The QR factorization of the data matrix (14) leads to the Cholesky factorization (16)Γ=KTK=RTRof the moment matrix [31](17)UWTWISP=RTR=R11TR11R11TR12R12TR11R12TR12+R222.The least squares estimation of the residual variance C^ and the parameter vector A^ are computed from the Cholesky factorization by(18)C^=R222N-p+1.A^=R11-1R12T.
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 p-order AR model by(19)FPE=logR222-logNN-p+1N+p+1.
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 fest of the ISP module occurs in the integer multiples of twice the voltage frequency applied to the motor fl:(20)fest=2mfl,where m=1,2,…,∞ [32]. Since components of the fault spectral signature decrease while m increases, just the main component (m=1) is considered; hence(21)fest=2fl.
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, three-phase 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 three-phase inverter was used for varying the motor rotating speed; a Hall-effect 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 32-bit floating point digital signal processor (DSP) was used for performing the necessary control and measurement functions.
Block diagram of experimental setup.
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 4096-sample frames from each phase. The acquired signals were sent to a PC via high speed USB. A computer-program 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.
Three-phase induction motor stator-winding structure.
The electric structure for the stator-winding 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 22-AGW wires for a total of 7290 turns.
The low-voltage operational condition for the induction motor was achieved by connecting in parallel the double star windings and supplying the phase voltages VAN to terminals 1 and 7, VBN to 2 and 8, and VCN to 3 and 9 and joining terminals 4, 5, and 6 to conform the neutral star YN, as specified by the manufacturer.
The fault insertion in the stator-winding is accomplished by connecting or disconnecting some segments of phase C winding to phase voltage VCN or neutral YN, using tap derivations C1,…,C8. This is achieved by an insulated electric circuit containing a group of relays that are managed by the DSP-based control system to select the fault severity in the motor. Possible connections are depicted with dotted lines in Figure 3.
Modified stator-winding with tap derivations insertion.
The characterization of the induction motor parameters was performed at line frequency fl (60 Hz) via the three-voltmeter method, which consists in inserting a pattern resistance RP in each phase of the motor and measuring stator voltages V1, V2, and V3 and currents as shown in Figure 4; then, the impedance ZL, inductive reactance XL, resistance RL, and inductance L for each phase are calculated. The results are shown in Table 2.
Stator-winding parameters.
Parameter
Phase A
Phase B
Phase C
L (mH)
5.055
5.025
5.035
ZL (Ω)
1.909
1.899
1.902
XL (Ω)
1.905
1.894
1.898
RL (Ω)
0.127
0.128
0.126
Characterization of the induction stator-winding parameters.
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 distinct-degree stator fault conditions as shown in Table 3. The healthy condition is obtained by VCN to terminal 3 and YN to terminal 6 (see Figure 3). In the first faulty condition (three turns), VCN is connected to C1 terminal. In the second faulty condition (nine turns) YN is connected to C8 terminal. Both faulty cases correspond to incipient short circuit stator faults.
Stator-winding parameters for healthy and fault cases.
Phase Cparameter
Healthy
Case 1Δ%
Case 2Δ%
L (mH)
5.035
6.494
8.639
ZL (Ω)
1.902
6.519
8.622
XL (Ω)
1.898
6.480
8.640
RL (Ω)
0.126
7.936
12.698
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 fs=4000 Hz. A load torque of 0.4 Nm was applied for all experiments at different speeds: 1200 rpm (fl=40 Hz), 1500 rpm (fl=50 Hz), and 1800 rpm (fl=60 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 (fl=60 Hz), evaluated for different orders (p=2,…,70), 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 (p=3,…,8) is minimum; however, the gap increases from order p=11 and it remains as the order goes up.
Mean FPE value for three fault scenarios at a speed of 1800 rpm.
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 35-order AR model are summarized in Table 4.
Mean FPE value of order p=35 for three fault scenarios at different speeds.
fl (Hz)
Healthy
Case 1Δ%
Case 2Δ%
60
−5.1151
8.06
11.59
50
−5.2309
5.70
7.63
40
−5.3140
5.51
7.26
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 (fl=60 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.
Mean value of the spectrum signature for healthy and faulty cases at different rotating speeds.
fl (Hz)
Healthy
Case 1Δ%
Case 2Δ%
60
0.045
583.59
858.21
50
0.033
568.77
856.36
40
0.047
511.59
756.01
Spectrum for healthy and faulty cases at 1800 rpm.
4.3. Discussion of Results
From the obtained results, the proposed methodology does not require a power-two N or adjusting the sampling frequency fs depending on the motor rotating speed. On the other hand, the common spectral analysis technique for detecting fault signatures does require fs to be a multiple of fl [28]; furthermore, in the case of N not being a power of two, zero padding must have to be done, and the technique will fail.
For a fixed p-order AR model with a defined sampling period, for instance, p=35 and N=4096, the FPE computation and estimation of the AR model coefficients will require (N+p+2)(p)2 = 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 i7-4710HQ, and it took around 0.00401 seconds, whereas the DFT technique requires N2 = 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 p≥12 (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 high-order AR model.
On the other hand, it was possible to detect the presence of fault conditions even with low order models (p=1), but it was not possible to quantify the magnitude of the fault, whereas, for a higher-order model and given the scope of this research that only characterized the electrical stator parameters for fl=60 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 3-sigma 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 3-sigma 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.
Experimental FPE normal distributions.
5. Conclusions
Electric motors are responsible for around 90% of the electric-machine 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 different-degree 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 3-sigma 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/CARepPrel-15 and Universidad Autónoma de San Luis Potosí under Grant C14-FAI-04-42.42. Mario A. Gonzalez-Garcia is a Research Fellow of the National Council of Science and Technology (CONACyT), Mexico.
ToliyatH. A.NandiS.ChoiS.Meshging-KelkH.2013Boca Raton, Fla, USACRC PressHenaoH.CapolinoG.-A.Fernandez-CabanasM.FilippettiF.BruzzeseC.StrangasE.PuscaR.EstimaJ.Riera-GuaspM.Hedayati-KiaS.Trends in fault diagnosis for electrical machines: a review of diagnostic techniques201482314210.1109/mie.2013.22876512-s2.0-84903176130NandiS.ToliyatH. A.LiX.Condition monitoring and fault diagnosis of electrical motors—a review200520471972910.1109/tec.2005.8479552-s2.0-29244467991BetinF.CapolinoG.-A.CasadeiD.KawkabaniB.BojoiR. I.HarneforsL.LeviE.ParsaL.FahimiB.Trends in electrical machines control: samples for classical, sensorless, and fault-tolerant techniques201482435510.1109/mie.2014.23137522-s2.0-84903151489De AngeloC. H.BossioG. R.GiacconeS. J.VallaM. I.SolsonaJ. A.GarcíaG. O.Online model-based stator-fault detection and identification in induction motors200956114671468010.1109/TIE.2009.20124682-s2.0-70350413799CruzS. M. A.CardosoA. J. M.Multiple reference frames theory: a new method for the diagnosis of stator faults in three-phase induction motors200520361161910.1109/tec.2005.8479752-s2.0-24644433498YangZ.Automatic condition monitoring of industrial rolling-element bearings using motor's vibration and current analysis201520151248615910.1155/2015/4861592-s2.0-84944199211ZhuH.-Y.HuJ.-T.GaoL.HuangH.Practical aspects of broken rotor bars detection in PWM voltage-source-inverter-fed squirrel-cage induction motors201320131112836810.1155/2013/1283682-s2.0-84893861066ShnibhaR.AlbarbarA.AbouhnikA.IbrahimG.A more reliable method for monitoring the condition of three-phase induction motors based on their vibrations20122012923031410.5402/2012/230314BonnettA. H.YungC.Increased efficiency versus increased reliability2008141293610.1109/MIA.2007.9098022-s2.0-39849088301Hernandez-VargasM.Cabal-YepezE.Garcia-PerezA.Real-time SVD-based detection of multiple combined faults in induction motors20144072193220310.1016/j.compeleceng.2013.12.0202-s2.0-84922794071ZhangP.DuY.HabetlerT. G.LuB.A survey of condition monitoring and protection methods for medium-voltage induction motors2011471344610.1109/tia.2010.20908392-s2.0-78751666101BellR. N.McWilliamsD. W.O'DonnellP.SinghC.WellsS. J.Report of large motor reliability survey of industrial and commercial installations, part I1985IA-21485386410.1109/tia.1985.3495322-s2.0-0022206607Report of large motor reliability survey of industrial and commercial installations, Part II198521486587210.1109/TIA.1985.349533Report of large motor reliability survey of industrial and commercial installations: part 3198721486587210.1109/TIA.1985.349533AlbrechtP. F.AppiariusJ. C.McCoyR. M.OwenE. L.SharmaD. K.Assessment of the reliability of motors in utility applications—updated198611394610.1109/TEC.1986.4765668SeinschH. O.Monitoring und diagnose elektrischer maschinen und antriebe2001VDEJiangL.-L.YinH.-K.LiX.-J.TangS.-W.Fault diagnosis of rotating machinery based on multisensor information fusion using SVM and time-domain features20142014841817810.1155/2014/4181782-s2.0-84898026077Delgado-ArredondoP. A.Garcia-PerezA.Morinigo-SoteloD.Osornio-RiosR. A.Avina-CervantesJ. G.Rostro-GonzalezH.Romero-TroncosoR. D. J.Comparative study of time-frequency decomposition techniques for fault detection in induction motors using vibration analysis during startup transient201520151470803410.1155/2015/7080342-s2.0-84937798836MortonJ.TorrioneP. A.CollinsL. M.Variational Bayesian learning for mixture autoregressive models with uncertain-order20115962614262710.1109/TSP.2011.2128310MR28572752-s2.0-79957464625WangG.LuoZ.QinX.LengY.WangT.Fault identification and classification of rolling element bearing based on time-varying autoregressive spectrum200822493494710.1016/j.ymssp.2007.10.0082-s2.0-39649085377LiuC.ZhouX.YangS.LiangW.MiaoQ.Cooling fan bearing diagnosis based on AR& MEDProceedings of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering (ICQR2MSE '12)June 2012Chengdu, China62262610.1109/icqr2mse.2012.62463102-s2.0-84866606189ZhouG.-H.ZuoC.-C.WangJ.-Z.LiuS.-X.Gearbox fault diagnosis based on wavelet-AR model2Proceedings of the 6th International Conference on Machine Learning and Cybernetics (ICMLC '07)August 2007Hong Kong1061106510.1109/icmlc.2007.43703002-s2.0-37849026776NikharN. K.PatankarS. S.KulkarniJ. V.Gear tooth fault detection by autoregressive modellingProceedings of the International Conference on Computing, Communications and Networking Technologies (ICCCNT '13)July 2013Tiruchengode, India16Villalobos-PiñaF. J.Alvarez-SalasR.Robust algorithm for electric fault diagnosis in the three phase induction machine based on spectral and wavelet tools201512329230310.1016/j.riai.2015.04.0032-s2.0-84941899838MilanezD. L.EmanuelA. E.The instantaneous-space-phasor: a powerful diagnosis tool200352114314810.1109/tim.2003.8090692-s2.0-0037310218LinX.XieY.ZhaoD.XuS.Estimation of observer parameters for dynamic positioning ships20132013717360310.1155/2013/173603MR3035629ProakisJ. G.ManolakisD. G.2007New Jersey, NJ, USAPearson Prentice HallHayesM. H.1996New York, NY, USAJohn Wiley & SonsNeumaierA.SchneiderT.Estimation of parameters and eigenmodes of multivariate autoregressive models2001271275710.1145/382043.3823042-s2.0-0002537923WuH.-C.ChangS. Y.Le-NgocT.WuY.Efficient rank-adaptive least-square estimation and multiple-parameter linear regression using novel dyadically recursive hermitian matrix inversion201220121089193210.1155/2012/8919322-s2.0-84870213164DrifM.CardosoA. J. M.Airgap-eccentricity fault diagnosis, in three-phase induction motors, by the complex apparent power signature analysis20085531404141010.1109/TIE.2007.9090762-s2.0-41349083807