# Stator Fault Detection in Induction Motors by Autoregressive Modeling

Authors: Francisco M. Garcia-Guevara; Francisco J. Villalobos-Piña; Ricardo Alvarez-Salas; Eduardo Cabal-Yepez; Mario A. Gonzalez-Garcia
Journal: Mathematical Problems in Engineering
 (2016)

---

## 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α
-
β 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.


---

## Body

## 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.Table 1
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
8Spectral 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 acquiringN samples of the induction motor three phase stator currents i
a, i
b, and i
c:(1)
i
a
i
b
i
c
=
i
a
1
i
a
2
⋯
i
a
N
i
b
1
i
b
2
⋯
i
b
N
i
c
1
i
c
2
⋯
i
c
N
.The three currents are mapped into the α
-
β stator-fixed reference frame using Concordia transformation: (2)
i
α
i
β
i
o
=
1
2
1
cos
⁡
-
2
π
3
cos
⁡
2
π
3
0
-
sin
⁡
-
2
π
3
-
sin
⁡
2
π
3
1
2
1
2
1
2
i
a
i
b
i
c
.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 ap-order AR model are estimated. The AR model is a finite impulse response filter (FIR) that can be defined as follows:(4)
I
S
P
n
=
∑
k
=
1
p
A
k
I
S
P
n
-
k
+
ε
n
,
n
=
1
,
…
,
N
,where A
k contains the coefficients of the AR model and I
S
P
n
-
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)
IS
P
n
=
A
p
u
n
+
ε
n
with
ε
n
=
noise
Cwith a coefficients vector(6)
A
1
,
…
,
A
pand with predictors (7)
u
n
=
1
IS
P
n
-
1
⋯
IS
P
n
-
p
T
.In terms of moments(8)
U
=
∑
n
=
1
N
u
n
u
n
T
,
ISP
=
∑
n
=
1
N
IS
P
n
IS
P
n
T
,
W
=
∑
n
=
1
N
IS
P
n
u
n
T
.The least squared estimates of coefficients A
^ are estimated by(9)
A
^
=
W
U
-
1
.The residual variance(10)
C
^
=
1
N
-
p
∑
n
=
1
N
ε
^
n
ε
^
n
T
with
ε
^
n
=
IS
P
n
-
A
^
p
u
nis an estimate of the noise covariance C and can be expressed as (11)
C
^
=
1
N
-
p
+
1
ISP
-
W
U
-
1
W
T
.The residual variance is proportional to Schur complement of the matrix (12)
Γ
=
U
W
T
W
ISP
=
∑
n
=
1
N
u
n
IS
P
n
u
n
T
IS
P
n
Twhich is the moment matrix Γ
=
K
T
K belonging to the data matrix(13)
K
=
u
1
T
IS
P
1
T
⋮
⋮
u
N
T
IS
P
N
T
,where the data matrix K
∈
R
N
-
p
×
p
+
1.The least squares estimation (9) is obtained from a QR factorization of K through the Householder algorithm:(14)
K
=
Q
R
,where Q and R are orthogonal upper triangular matrices, and R is given by(15)
R
=
R
11
R
12
0
R
22
,where R
∈
R
N
+
1
×
p, R
11
∈
R
p
×
p, R
12
∈
R
1
×
p, and R
22
∈
R.The QR factorization of the data matrix (14) leads to the Cholesky factorization (16)
Γ
=
K
T
K
=
R
T
Rof the moment matrix [31](17)
U
W
T
W
ISP
=
R
T
R
=
R
11
T
R
11
R
11
T
R
12
R
12
T
R
11
R
12
T
R
12
+
R
22
2
.The least squares estimation of the residual variance C
^ and the parameter vector A
^ are computed from the Cholesky factorization by(18)
C
^
=
R
22
2
N
-
p
+
1
.
A
^
=
R
11
-
1
R
12
T
.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)
F
P
E
=
log
⁡
R
22
2
-
log
⁡
N
N
-
p
+
1
N
+
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 signaturef
e
s
t of the ISP module occurs in the integer multiples of twice the voltage frequency applied to the motor f
l:(20)
f
e
s
t
=
2
m
f
l
,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)
f
est
=
2
f
l
.


## 2.1. AR Model

The proposed electric stator fault diagnosis scheme consists in acquiringN samples of the induction motor three phase stator currents i
a, i
b, and i
c:(1)
i
a
i
b
i
c
=
i
a
1
i
a
2
⋯
i
a
N
i
b
1
i
b
2
⋯
i
b
N
i
c
1
i
c
2
⋯
i
c
N
.The three currents are mapped into the α
-
β stator-fixed reference frame using Concordia transformation: (2)
i
α
i
β
i
o
=
1
2
1
cos
⁡
-
2
π
3
cos
⁡
2
π
3
0
-
sin
⁡
-
2
π
3
-
sin
⁡
2
π
3
1
2
1
2
1
2
i
a
i
b
i
c
.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 ap-order AR model are estimated. The AR model is a finite impulse response filter (FIR) that can be defined as follows:(4)
I
S
P
n
=
∑
k
=
1
p
A
k
I
S
P
n
-
k
+
ε
n
,
n
=
1
,
…
,
N
,where A
k contains the coefficients of the AR model and I
S
P
n
-
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)
IS
P
n
=
A
p
u
n
+
ε
n
with
ε
n
=
noise
Cwith a coefficients vector(6)
A
1
,
…
,
A
pand with predictors (7)
u
n
=
1
IS
P
n
-
1
⋯
IS
P
n
-
p
T
.In terms of moments(8)
U
=
∑
n
=
1
N
u
n
u
n
T
,
ISP
=
∑
n
=
1
N
IS
P
n
IS
P
n
T
,
W
=
∑
n
=
1
N
IS
P
n
u
n
T
.The least squared estimates of coefficients A
^ are estimated by(9)
A
^
=
W
U
-
1
.The residual variance(10)
C
^
=
1
N
-
p
∑
n
=
1
N
ε
^
n
ε
^
n
T
with
ε
^
n
=
IS
P
n
-
A
^
p
u
nis an estimate of the noise covariance C and can be expressed as (11)
C
^
=
1
N
-
p
+
1
ISP
-
W
U
-
1
W
T
.The residual variance is proportional to Schur complement of the matrix (12)
Γ
=
U
W
T
W
ISP
=
∑
n
=
1
N
u
n
IS
P
n
u
n
T
IS
P
n
Twhich is the moment matrix Γ
=
K
T
K belonging to the data matrix(13)
K
=
u
1
T
IS
P
1
T
⋮
⋮
u
N
T
IS
P
N
T
,where the data matrix K
∈
R
N
-
p
×
p
+
1.The least squares estimation (9) is obtained from a QR factorization of K through the Householder algorithm:(14)
K
=
Q
R
,where Q and R are orthogonal upper triangular matrices, and R is given by(15)
R
=
R
11
R
12
0
R
22
,where R
∈
R
N
+
1
×
p, R
11
∈
R
p
×
p, R
12
∈
R
1
×
p, and R
22
∈
R.The QR factorization of the data matrix (14) leads to the Cholesky factorization (16)
Γ
=
K
T
K
=
R
T
Rof the moment matrix [31](17)
U
W
T
W
ISP
=
R
T
R
=
R
11
T
R
11
R
11
T
R
12
R
12
T
R
11
R
12
T
R
12
+
R
22
2
.The least squares estimation of the residual variance C
^ and the parameter vector A
^ are computed from the Cholesky factorization by(18)
C
^
=
R
22
2
N
-
p
+
1
.
A
^
=
R
11
-
1
R
12
T
.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)
F
P
E
=
log
⁡
R
22
2
-
log
⁡
N
N
-
p
+
1
N
+
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 signaturef
e
s
t of the ISP module occurs in the integer multiples of twice the voltage frequency applied to the motor f
l:(20)
f
e
s
t
=
2
m
f
l
,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)
f
est
=
2
f
l
.


## 3. Experimental Setup

The proposed methodology was validated on an experimental platform. The implemented test bed for performing the experimental analysis is depicted in Figure1. 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.Figure 1
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 Figure2.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 voltagesV
A
N to terminals 1 and 7, V
B
N to 2 and 8, and V
C
N to 3 and 9 and joining terminals 4, 5, and 6 to conform the neutral star Y
N, 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 voltageV
CN or neutral Y
N, using tap derivations C
1
,
…
,
C
8. 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.Figure 3
Modified stator-winding with tap derivations insertion.The characterization of the induction motor parameters was performed at line frequencyf
l (60 Hz) via the three-voltmeter method, which consists in inserting a pattern resistance R
P in each phase of the motor and measuring stator voltages V
1, V
2, and V
3 and currents as shown in Figure 4; then, the impedance Z
L, inductive reactance X
L, resistance R
L, and inductance L for each phase are calculated. The results are shown in Table 2.Table 2
Stator-winding parameters.
Parameter
Phase A
Phase B
Phase C
L (mH)
5.055
5.025
5.035
Z
L (Ω)
1.909
1.899
1.902
X
L (Ω)
1.905
1.894
1.898
R
L (Ω)
0.127
0.128
0.126Figure 4
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 Table3. The healthy condition is obtained by V
CN to terminal 3 and Y
N to terminal 6 (see Figure 3). In the first faulty condition (three turns), V
CN is connected to C
1 terminal. In the second faulty condition (nine turns) Y
N is connected to C
8 terminal. Both faulty cases correspond to incipient short circuit stator faults.Table 3
Stator-winding parameters for healthy and fault cases.
Phase Cparameter
Healthy
Case 1
Δ%
Case 2
Δ%
L (mH)
5.035
6.494
8.639
Z
L (Ω)
1.902
6.519
8.622
X
L (Ω)
1.898
6.480
8.640
R
L (Ω)
0.126
7.936
12.698The 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 frequencyf
s
=
4
000 Hz. A load torque of 0.4 Nm was applied for all experiments at different speeds: 1200 rpm (f
l
=
4
0 Hz), 1500 rpm (f
l
=
5
0 Hz), and 1800 rpm (f
l
=
6
0 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 (f
l
=
6
0 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.Figure 5
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 Table4.Table 4
Mean FPE value of orderp
=
35 for three fault scenarios at different speeds.
f
l (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

Figure6 shows the stator fault spectral signature corresponding to the three considered conditions reported in Table 5, at the rotating speed of 1800 rpm (f
l
=
6
0 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.Table 5
Mean value of the spectrum signature for healthy and faulty cases at different rotating speeds.
f
l (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.01Figure 6
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-twoN or adjusting the sampling frequency f
s depending on the motor rotating speed. On the other hand, the common spectral analysis technique for detecting fault signatures does require f
s to be a multiple of f
l [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 fixedp-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 N
2 = 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 orderp
≥
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 f
l
=
6
0 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 Figure7; 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.Figure 7
Experimental FPE normal distributions.


## 4.1. AR Model Approach

The mean FPE values of the ISP module from the three treated cases at rotating speed of 1800 rpm (f
l
=
6
0 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.Figure 5
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 Table4.Table 4
Mean FPE value of orderp
=
35 for three fault scenarios at different speeds.
f
l (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

Figure6 shows the stator fault spectral signature corresponding to the three considered conditions reported in Table 5, at the rotating speed of 1800 rpm (f
l
=
6
0 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.Table 5
Mean value of the spectrum signature for healthy and faulty cases at different rotating speeds.
f
l (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.01Figure 6
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-twoN or adjusting the sampling frequency f
s depending on the motor rotating speed. On the other hand, the common spectral analysis technique for detecting fault signatures does require f
s to be a multiple of f
l [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 fixedp-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 N
2 = 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 orderp
≥
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 f
l
=
6
0 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 Figure7; 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.Figure 7
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.



---

*Source: 3409756-2016-03-30.xml*