# Vibration Analysis of Partially Damaged Rotor Bar in Induction Motor under Different Load Condition Using DWT

Authors: Jose Rangel-Magdaleno; Hayde Peregrina-Barreto; Juan Ramirez-Cortes; Roberto Morales-Caporal; Israel Cruz-Vega
Journal: Shock and Vibration
 (2016)

---

## Abstract

The relevance of the development of monitoring systems for rotating machines is not only the ability to detect failures but also how early these failures can be detected. The purpose of this paper is to present an experimental study of partially damaged rotor bar in induction motor under different load conditions based on discrete wavelet transform analysis. The approach is based on the extraction of features from vibration signals at different level of damage and three mechanical load conditions. The proposed analysis is reliable for tracking the damage in rotor bar. The paper presents an analysis and extraction of vibration features for partially damaged rotor bar in induction motors. The experimental analysis shows the change in behavior of vibration due to load condition and progressive damage.


---

## Body

## 1. Introduction

The increasing need for production at low cost in a safe manner and high quality has accelerated a change in rotating machine maintenance from corrective to preventive [1]. The relevance of the development of monitoring systems for rotating machines is not only about the ability to detect failures but also about how early these failures can be detected [2]. As one of the major equipment items in modern industry, induction motor constitutes around 85% of the motors [3]. A common fault in induction motors is a broken rotor bar; this fault corresponds approximately to 10% of the failures [4]. In addition, this failure significantly increases power consumption and it is responsible for further damage to the machinery because of the undesired vibrations and shattering that are produced.The most popular techniques for fault detection in induction motors are the motor current signature analysis (MCSA) and the motor vibration signature analysis (MVSA). There are several current analysis techniques in literature about the identification of faults in induction motors such as broken bars [2, 5–11] and also vibration analysis techniques [12–25]. When it comes to vibration analysis techniques, there are several methodologies for the detection of broken bars. Su et al. analyzed the vibration signal using an Artificial Neural Network (ANN) for motor fault detection [12]. Gritli et al. investigated the impact of the control system related to fault diagnosis; they concluded that current signal is affected and vibration analysis showed better result for fault diagnosis [13]. Betta et al. proposed a vibration analyzer based on a digital signal processor (DSP) for fault diagnosis in rotating machine [16]. Sadoughi et al. proposed an intelligent diagnosis based on ANN and new features in vibration spectrum [19]. Climente-Alarcon et al. studied the effect of interbar currents attenuation and how they are greatly reduced if the bar is severed at both ends. They used vibration spectrum analysis during stationary operation [25]. Several authors have used wavelet decomposition for broken bar detection. Antonio-Daviu et al. [20] developed an evaluation of different starting methods for motor fault detection; they showed how the advanced analysis based on startup technique relies on identifying the time-frequency evolutions of the rotor-fault-related components contained in the startup current signal. Miceli et al. [24] presented a time-frequency analysis from axial and radial core vibration signal for broken rotor bar detection. Ordaz-Moreno et al. [6] developed a methodology based on discrete wavelet transform (DWT) for broken bar detection during transient state using current signature analysis. Delgado-Arredondo et al. [15] presented a comparative study of time-frequency decomposition techniques for fault detection in induction motor using vibration analysis. A review of wavelet analysis applied to induction motor fault detection is presented in [26]. The works presented above are calibrated to one or two broken bar detections, in the better case a half broken bar can be detected, and the reason is due to the fact that amplitude of incipient broken bar is small and can be absorbed by the fundamental frequency in traditional MCSA methodology. Recently DWT has been used remarking the importance of wavelets analysis in induction motor fault detection, but the majority detects only one or half broken bar. Garcia-Perez et al. [27] proposed the use of MUSIC (Multiple Signal Classification) for analysis of incipient broken rotor bar in induction motor using current analysis. Their work represents a first study in partially broken bars and their detection using motor current signature analysis during the transient state. Nevertheless, their methodology was not tested with load variation, and then more investigation related to incipient broken bars is mandatory.In this paper, the use of DWT and MVSA for the study of partially damaged rotor bar is presented. The contribution of the paper is the analysis of some features extracted from wavelets decomposition that allows detecting the change in vibration behavior due to progressive damage. The novelty of the paper is the detection of incipient broken bars (3 mm, 5 mm, 7 mm, and 10 mm) using vibration analysis during startup transient. The experimental analysis under different mechanical load levels shows the change in behavior and vibration due to progressive damage.The remainder of this paper is organized as follows: Section2 presents the theoretical background of motor vibration signature analysis, discrete wavelet transform, autocorrelation, and some estimation parameters; the experimental setup is presented in Section 3; Section 4 presents results from several cases of study. The conclusions are presented in Section 5.


## 2. Theoretical Background

The identification of incipient broken bars implies an appropriate feature extraction. The background of broken bar detection and the techniques used for it are presented in this section.


### 2.1. Motor Vibration Signature Analysis (MVSA)

Several asymmetries on the rotor can be determined by the spectrum analysis of the motor vibration. A broken bar in an induction motor can be detected through MVSA. A broken bar induces a spurious component in the vibration spectrum given by [2](1)fBB=1±2ksfss,where k is an integer, s is the per-unit slip, and fss is the supply frequency. Nevertheless, due to the small amplitude of frequencies of partially broken bar its detection is compromised. For this reason different methodologies have been proposed [2, 27].


### 2.2. Discrete Wavelet Transform (DWT)

Discrete wavelet transform has widely been used for signal processing of nonstationary or transient signals. Mallat’s algorithm is used to compute the DWT based on a set of discrete-time low- and high-pass filters. Figure1 shows Mallat’s algorithm for DWT decomposition [28]. The decomposition produced the approximation (low-frequency component) and the detail (high-frequency component) coefficients. According to the DWT filter bank properties, the frequency bands for approximation ACL and detail DCL at level L are given by(2)ACL=0,fs2L+1,DCL=fs2L+1,fs2L,where fs is the sampling frequency. Table 1 shows an example of the frequency bands for detail and approximation levels.Table 1
DWT decomposition level frequency bands for a sampling frequencyfs = 3200 Hz.
Level (L)
Approximationbandwidth (Hz)
Detailbandwidth (Hz)
1
AC1  [0, 800]
DC1  [800, 1600]
2
AC2  [0, 400]
DC2  [400, 800]
3
AC3  [0, 200]
DC3  [200, 400]
4
AC4  [0, 100]
DC4  [100, 200]
5
AC5  [0, 50]
DC5  [50, 100]
6
AC6  [0, 25]
DC6  [25, 50]
7
AC7  [0, 12.5]
DC7  [12.5, 25]Figure 1
Mallat algorithm for DWT decomposition.Figure2 shows an example of DWT decomposition of an induction motor vibration signal at healthy condition and full mechanical load. Seven levels were used for the decomposition. Figure 2(a) shows the approximation levels and Figure 2(b) shows the detail levels. Ordaz-Moreno et al. [6] and Antonino-Daviu et al. [20] proposed the analysis of current and vibrations signals during transient state, respectively. They analyzed the variability of the signal at different levels; nevertheless, for incipient broken bars the direct analysis of the signal is not enough. In this paper, an additional preprocessing step by means of autocorrelation of the signals is proposed; this helps to reduce the noise in acquired vibration signals.Figure 2
DWT decomposition, (a) approximation, and (b) detail.
(a)
(b)


### 2.3. Autocorrelation

The autocorrelation is defined as the cross-correlation of the signal with itself. The autocorrelation allows finding repeating patterns and measures the linear relationship between time series observations separated by a lag ofn time units. The autocorrelation has shown good results in fault detection [29, 30]. The autocorrelation is given by(3)Cxx=EXn+mXn∗,where X is the vibration signal and E is the spectral estimation. Figure 3 shows an example of DWT decomposition of a healthy induction motor vibration signal at full mechanical load previously preprocessed by autocorrelation. Figure 3(a) shows the approximation levels and Figure 3(b) shows the detail levels.Figure 3
Autocorrelation and DWT decomposition for a healthy motor, (a) approximation, and (b) detail.
(a)
(b)An example of a DWT decomposition of a damage induction motor vibration signal is shown in Figure4; it was obtained at full mechanical load and previously preprocessed by means of autocorrelation. The motor has one broken rotor bar. The approximation and the detail levels are shown in Figures 4(a) and 4(b), respectively. Figure 3(a) at level 7 shows the expected triangle due to the autocorrelation process, whereas Figure 4(a) at level 7 shows a noisy triangle. The difference between a healthy and damaged motor signals is visually clear. The only difference in condition between the two motors is the broken rotor bar. Then, the signal contamination is due to the broken bar. In order to obtain a better interpretation the use of some basic statistic parameters like mean, zero-crossing, and standard deviation at level 3 of decomposition is proposed. In this level the range of frequencies is from 0 to 200 Hz; then the fundamental and some harmonics are included.Figure 4
Autocorrelation and DWT decomposition for a motor with one broken bar, (a) approximation, and (b) detail.
(a)
(b)


### 2.4. Estimation Parameters

Three basic statistic parameters were used, the mean (x^) defined in (4), standard deviation (SD) defined by (5), and zero-crossing (ZC) that counts the number of times the signal crosses the zero value or presents a change in sign:(4)x^=1N∑i=0Nxi,(5)SD=1N∑i=0Nxi-x^21/2.Also the variability, norm, kurtosis, skewness, and RMS value were computed; nevertheless, parameters like the mean, standard deviation, and zero-crossing provided better results of the information we expected to obtain.


### 2.5. Classification:k-Nearest-Neighbor Algorithm (kNN)

A common reference among the classifiers due to the simplicity of the learning rule is thek-nearest-neighbor classifier. This classifier is memory-based and requires no model to be fit. The basic idea is that, given a query point x0, it is found that k training points xr, r=1,…,k, are closest in distance to x0, and they are classified using the majority vote among the k neighbors [31]. Assuming that the features of the objects to be classified are real-value, a common distance measure in the feature space is the Euclidean, and all objects in the training set participate to determine k-nearest-neighbors to the query point. A common classifier is the 1-nearest-neighbor, which is used in this paper.


## 2.1. Motor Vibration Signature Analysis (MVSA)

Several asymmetries on the rotor can be determined by the spectrum analysis of the motor vibration. A broken bar in an induction motor can be detected through MVSA. A broken bar induces a spurious component in the vibration spectrum given by [2](1)fBB=1±2ksfss,where k is an integer, s is the per-unit slip, and fss is the supply frequency. Nevertheless, due to the small amplitude of frequencies of partially broken bar its detection is compromised. For this reason different methodologies have been proposed [2, 27].


## 2.2. Discrete Wavelet Transform (DWT)

Discrete wavelet transform has widely been used for signal processing of nonstationary or transient signals. Mallat’s algorithm is used to compute the DWT based on a set of discrete-time low- and high-pass filters. Figure1 shows Mallat’s algorithm for DWT decomposition [28]. The decomposition produced the approximation (low-frequency component) and the detail (high-frequency component) coefficients. According to the DWT filter bank properties, the frequency bands for approximation ACL and detail DCL at level L are given by(2)ACL=0,fs2L+1,DCL=fs2L+1,fs2L,where fs is the sampling frequency. Table 1 shows an example of the frequency bands for detail and approximation levels.Table 1
DWT decomposition level frequency bands for a sampling frequencyfs = 3200 Hz.
Level (L)
Approximationbandwidth (Hz)
Detailbandwidth (Hz)
1
AC1  [0, 800]
DC1  [800, 1600]
2
AC2  [0, 400]
DC2  [400, 800]
3
AC3  [0, 200]
DC3  [200, 400]
4
AC4  [0, 100]
DC4  [100, 200]
5
AC5  [0, 50]
DC5  [50, 100]
6
AC6  [0, 25]
DC6  [25, 50]
7
AC7  [0, 12.5]
DC7  [12.5, 25]Figure 1
Mallat algorithm for DWT decomposition.Figure2 shows an example of DWT decomposition of an induction motor vibration signal at healthy condition and full mechanical load. Seven levels were used for the decomposition. Figure 2(a) shows the approximation levels and Figure 2(b) shows the detail levels. Ordaz-Moreno et al. [6] and Antonino-Daviu et al. [20] proposed the analysis of current and vibrations signals during transient state, respectively. They analyzed the variability of the signal at different levels; nevertheless, for incipient broken bars the direct analysis of the signal is not enough. In this paper, an additional preprocessing step by means of autocorrelation of the signals is proposed; this helps to reduce the noise in acquired vibration signals.Figure 2
DWT decomposition, (a) approximation, and (b) detail.
(a)
(b)


## 2.3. Autocorrelation

The autocorrelation is defined as the cross-correlation of the signal with itself. The autocorrelation allows finding repeating patterns and measures the linear relationship between time series observations separated by a lag ofn time units. The autocorrelation has shown good results in fault detection [29, 30]. The autocorrelation is given by(3)Cxx=EXn+mXn∗,where X is the vibration signal and E is the spectral estimation. Figure 3 shows an example of DWT decomposition of a healthy induction motor vibration signal at full mechanical load previously preprocessed by autocorrelation. Figure 3(a) shows the approximation levels and Figure 3(b) shows the detail levels.Figure 3
Autocorrelation and DWT decomposition for a healthy motor, (a) approximation, and (b) detail.
(a)
(b)An example of a DWT decomposition of a damage induction motor vibration signal is shown in Figure4; it was obtained at full mechanical load and previously preprocessed by means of autocorrelation. The motor has one broken rotor bar. The approximation and the detail levels are shown in Figures 4(a) and 4(b), respectively. Figure 3(a) at level 7 shows the expected triangle due to the autocorrelation process, whereas Figure 4(a) at level 7 shows a noisy triangle. The difference between a healthy and damaged motor signals is visually clear. The only difference in condition between the two motors is the broken rotor bar. Then, the signal contamination is due to the broken bar. In order to obtain a better interpretation the use of some basic statistic parameters like mean, zero-crossing, and standard deviation at level 3 of decomposition is proposed. In this level the range of frequencies is from 0 to 200 Hz; then the fundamental and some harmonics are included.Figure 4
Autocorrelation and DWT decomposition for a motor with one broken bar, (a) approximation, and (b) detail.
(a)
(b)


## 2.4. Estimation Parameters

Three basic statistic parameters were used, the mean (x^) defined in (4), standard deviation (SD) defined by (5), and zero-crossing (ZC) that counts the number of times the signal crosses the zero value or presents a change in sign:(4)x^=1N∑i=0Nxi,(5)SD=1N∑i=0Nxi-x^21/2.Also the variability, norm, kurtosis, skewness, and RMS value were computed; nevertheless, parameters like the mean, standard deviation, and zero-crossing provided better results of the information we expected to obtain.


## 2.5. Classification:k-Nearest-Neighbor Algorithm (kNN)

A common reference among the classifiers due to the simplicity of the learning rule is thek-nearest-neighbor classifier. This classifier is memory-based and requires no model to be fit. The basic idea is that, given a query point x0, it is found that k training points xr, r=1,…,k, are closest in distance to x0, and they are classified using the majority vote among the k neighbors [31]. Assuming that the features of the objects to be classified are real-value, a common distance measure in the feature space is the Euclidean, and all objects in the training set participate to determine k-nearest-neighbors to the query point. A common classifier is the 1-nearest-neighbor, which is used in this paper.


## 3. Methodology

The methodology consists mainly in the following steps:(1)
Vibration signals of the three axes are sensed online from the motor and saved in a PC.(2)
Feature extraction stage is as follows:
(a)
An autocorrelation function is applied to each axis.(b)
The decomposition at level three of each signal is obtained using DWT algorithm.(c)
Statistical parameters like mean, standard deviation, and zero-crossings are calculated for the decomposition coefficients.(3)
Statistical analysis of the parameter is computed.(4)
Classification stage is as follows: the motor condition is classified based onkNN algorithm.


### 3.1. Experimental Setup

The experimental test was performed on a 1-HP motor, the vibration signal was acquired using MEMS-based triaxial accelerometer with a user-selectable full scale (±2 g/±6 g); a resolution of 12 bits and an oversampling frequency of 3.2 KHz was used. The signals were acquired at the startup transient state of the motor during 5 s. Figure5 shows the location and orientation of the accelerometer over the motor. The axis “x” is perpendicular to the rotor, where the best results to the analysis were obtained, axis “y” provided acceptable results, and axis “z” provided less satisfactory results.Figure 5
Accelerometer placed on induction motor.The data were recorded in the following states: healthy, 3 mm, 5 mm, 7 mm, and 1 broken bar. For each state, 50 runs for unload, 50 runs for half load, and 50 runs for full load were recorded. In order to obtain the incipient conditions (3 mm, 5 mm, 7 mm, and 1 br) the rotor was externally drilled at different depths. A total of 750 runs were analyzed in this study.


## 3.1. Experimental Setup

The experimental test was performed on a 1-HP motor, the vibration signal was acquired using MEMS-based triaxial accelerometer with a user-selectable full scale (±2 g/±6 g); a resolution of 12 bits and an oversampling frequency of 3.2 KHz was used. The signals were acquired at the startup transient state of the motor during 5 s. Figure5 shows the location and orientation of the accelerometer over the motor. The axis “x” is perpendicular to the rotor, where the best results to the analysis were obtained, axis “y” provided acceptable results, and axis “z” provided less satisfactory results.Figure 5
Accelerometer placed on induction motor.The data were recorded in the following states: healthy, 3 mm, 5 mm, 7 mm, and 1 broken bar. For each state, 50 runs for unload, 50 runs for half load, and 50 runs for full load were recorded. In order to obtain the incipient conditions (3 mm, 5 mm, 7 mm, and 1 br) the rotor was externally drilled at different depths. A total of 750 runs were analyzed in this study.


## 4. Analysis Results

As it was described above, the three axes were previously preprocessed by autocorrelation and then the DWT was applied and level 3 of decomposition was obtained. At this level the mean, ZC, and SD were applied. Each database (without load, half load, and full load) are constituted by 250 runs. The results obtained for each condition are shown in Figures6–17.Figure 6
Approximation coefficients at no mechanical load,Ax-axis.Figure 7
Unload condition analysis,Az-axis.Figure 8
Unload condition analysis,Ay-axis.Figure 9
Unload condition analysis,Ax-axis.Figure 10
Approximation coefficients at half-mechanical load,Ax-axis.Figure 11
Half-load condition analysis,Az-axis.Figure 12
Half-load condition analysis,Ay-axis.Figure 13
Half-load condition analysis,Ax-axis.Figure 14
Approximation coefficients at full mechanical load,Ax-axis.Figure 15
Full-load condition analysis,Az-axis.Figure 16
Full-load condition analysis,Ay-axis.Figure 17
Full-load condition analysis,Ax-axis.


### 4.1. Without Mechanical Load

Figure6 shows a section from the approximation coefficients at level 3 of decomposition. As can be seen, the coefficients show differences that can be evaluated by the proposed estimation parameters. Figure 6 shows the five conditions evaluated in this test. The three axes were analyzed to obtain better results.In Figure7, the analysis of the five conditions under test is shown for the Az-axis. The analysis is performed by means of the SD and mean value parameters. One can see that the five regions present some degree of overlap; nevertheless, for the healthy, 7 mm, and one-broken bar conditions the separation between classes is more evident, allowing a better classification.TheAy-axis analysis is depicted in Figure 8; the results are now obtained with the SD and ZC estimation parameters. As it can be seen, the conditions under test present a high dispersion degree and overlapping between them, which complicates the classification because none of the conditions are clearly separated.Finally, Figure9 shows the analysis of the Ax-axis, where the mean and ZC values were selected as estimation parameters. It is considerably clear that the Ax-axis shows the best results for classification; all the conditions are considerably grouped and separated. It is possible to detect not only the damage condition but also the level of such damage. The level of overlapping is too low showing that the analysis in the Ax-axis is the best option to obtain a clear classification of the damage level under a nonmechanical load condition.


### 4.2. Half-Mechanical Load

Figure10 presents a section of the approximation coefficients at level 3 of decomposition for Ax-axis; all the conditions are shown. Differences among amplitudes of the signals can be seen. For better results the three axes were analyzed; mean and SD were selected for the analysis. Figure 11 shows the analysis for Az-axis; in this case all the conditions are totally overlapped. This axis cannot be used for condition classification; even healthy and one-broken bar conditions are not identifiable.Figure12 shows the analysis for Ay-axis; ZC and SD were selected for the analysis. Two main groups are presented; one for healthy and 3 mm conditions and the other one for 5 mm, 7 mm, and 1 br conditions. Under this condition of half-mechanical load and for this axis, the 3 mm depth drilled is considered as healthy whereas the damage condition is detected at a depth of 5 mm. In the second group, using an appropriate classification algorithm, the one-broken bar condition can be detected. The 5 mm and 7 mm have more overlapping degree.Figure13 presents the analysis for Ax-axis; ZC and SD were selected for the analysis. In this case, healthy condition has less overlap with 3 mm. For the rest conditions it is difficult to obtain a clear result due to the highly overlapping degree. Then, Ay and Ax axes can be used together for a better classification; using both, the healthy and damage conditions (3 mm, 5 mm, 7 mm, and 1 br) can be better classified. Depending on the classification technique other conditions can be estimated to obtain healthy, medium damage (3 mm, 5 mm and 7 mm), and one-broken bar conditions.


### 4.3. Full-Mechanical Load

Figure14 shows a section of the approximation coefficients of Ax-axis for full-mechanical load. In the figure the differences among the five conditions are presented. In order to measure these differences the three axes were analyzed as is shown in Figures 15, 16, and 17.Figure15 shows the analysis of Az-axis; mean and SD were selected for the analysis. The conditions are totally overlapped. It is clear that this axis does not help in the analysis of partially broken bars.Figure16 presents the analysis for Ay-axis; mean and SD were selected for the analysis. Four groups in the figure can be seen, for the 1-bar condition and 5 mm of a depth drilled condition, a clearly division between classes is obtained, which does not happen with the 3 mm and healthy conditions and in certain degree of about 25% of the 7 mm condition; the Ax-axis is developed for this last group.Figure17 shows the analysis in Ax-axis; ZC and mean were the estimation parameters selected for the analysis. Three groups are showed in the figure, the first for healthy condition; the second for 3 mm; and the third is an overlap of 5 mm, 7 mm, and 1 br. Then the five conditions can be detected using a combination of Ay and Ax, Ay detects 1 br and 5 mm and 7 mm conditions with some overlap, and Ax detects healthy and 3 mm conditions. In the case of the overlap in 7 mm in Ay-axis (group four) it is necessary to analyze the third group in Ax (5 mm, 7 mm, and 1 br), if the sample belongs to this group then it can be classified as 7 mm, and we know that it is not 5 mm or 1 br because these conditions are perfectly separable in Ay-axis.


### 4.4.kNN

A simple classification technique was used to prove the usefulness of the analysis.k-nearest-neighbors (kNN) algorithm was used. Also tests adding white noise were made; the level of noise was based on the signal-to-noise ratio (SNR). The values used were SNR = 10, 5, and 1. Table 2 shows the results of the classification for the Ax-axis with k=1 and for the three cases of added noise. These results are obtained from a cross-validation process with 10-folds; the average loss for each cross-validation model is used to obtain the accuracy level in classification. We have to remark that the simple value of k=1 was sufficient to obtain a desirable accuracy level in classification.Table 2
Classification percentage for each condition usingkNN.
Condition/axis
A
x
SNR10
SNR5
SNR1
Unload
96.55%
96.70%
92.85%
89.8%
Half load
80.5%
81.5%
80.85%
80.55%
Full load
87.6%
85.85%
84.4%
85.15%


## 4.1. Without Mechanical Load

Figure6 shows a section from the approximation coefficients at level 3 of decomposition. As can be seen, the coefficients show differences that can be evaluated by the proposed estimation parameters. Figure 6 shows the five conditions evaluated in this test. The three axes were analyzed to obtain better results.In Figure7, the analysis of the five conditions under test is shown for the Az-axis. The analysis is performed by means of the SD and mean value parameters. One can see that the five regions present some degree of overlap; nevertheless, for the healthy, 7 mm, and one-broken bar conditions the separation between classes is more evident, allowing a better classification.TheAy-axis analysis is depicted in Figure 8; the results are now obtained with the SD and ZC estimation parameters. As it can be seen, the conditions under test present a high dispersion degree and overlapping between them, which complicates the classification because none of the conditions are clearly separated.Finally, Figure9 shows the analysis of the Ax-axis, where the mean and ZC values were selected as estimation parameters. It is considerably clear that the Ax-axis shows the best results for classification; all the conditions are considerably grouped and separated. It is possible to detect not only the damage condition but also the level of such damage. The level of overlapping is too low showing that the analysis in the Ax-axis is the best option to obtain a clear classification of the damage level under a nonmechanical load condition.


## 4.2. Half-Mechanical Load

Figure10 presents a section of the approximation coefficients at level 3 of decomposition for Ax-axis; all the conditions are shown. Differences among amplitudes of the signals can be seen. For better results the three axes were analyzed; mean and SD were selected for the analysis. Figure 11 shows the analysis for Az-axis; in this case all the conditions are totally overlapped. This axis cannot be used for condition classification; even healthy and one-broken bar conditions are not identifiable.Figure12 shows the analysis for Ay-axis; ZC and SD were selected for the analysis. Two main groups are presented; one for healthy and 3 mm conditions and the other one for 5 mm, 7 mm, and 1 br conditions. Under this condition of half-mechanical load and for this axis, the 3 mm depth drilled is considered as healthy whereas the damage condition is detected at a depth of 5 mm. In the second group, using an appropriate classification algorithm, the one-broken bar condition can be detected. The 5 mm and 7 mm have more overlapping degree.Figure13 presents the analysis for Ax-axis; ZC and SD were selected for the analysis. In this case, healthy condition has less overlap with 3 mm. For the rest conditions it is difficult to obtain a clear result due to the highly overlapping degree. Then, Ay and Ax axes can be used together for a better classification; using both, the healthy and damage conditions (3 mm, 5 mm, 7 mm, and 1 br) can be better classified. Depending on the classification technique other conditions can be estimated to obtain healthy, medium damage (3 mm, 5 mm and 7 mm), and one-broken bar conditions.


## 4.3. Full-Mechanical Load

Figure14 shows a section of the approximation coefficients of Ax-axis for full-mechanical load. In the figure the differences among the five conditions are presented. In order to measure these differences the three axes were analyzed as is shown in Figures 15, 16, and 17.Figure15 shows the analysis of Az-axis; mean and SD were selected for the analysis. The conditions are totally overlapped. It is clear that this axis does not help in the analysis of partially broken bars.Figure16 presents the analysis for Ay-axis; mean and SD were selected for the analysis. Four groups in the figure can be seen, for the 1-bar condition and 5 mm of a depth drilled condition, a clearly division between classes is obtained, which does not happen with the 3 mm and healthy conditions and in certain degree of about 25% of the 7 mm condition; the Ax-axis is developed for this last group.Figure17 shows the analysis in Ax-axis; ZC and mean were the estimation parameters selected for the analysis. Three groups are showed in the figure, the first for healthy condition; the second for 3 mm; and the third is an overlap of 5 mm, 7 mm, and 1 br. Then the five conditions can be detected using a combination of Ay and Ax, Ay detects 1 br and 5 mm and 7 mm conditions with some overlap, and Ax detects healthy and 3 mm conditions. In the case of the overlap in 7 mm in Ay-axis (group four) it is necessary to analyze the third group in Ax (5 mm, 7 mm, and 1 br), if the sample belongs to this group then it can be classified as 7 mm, and we know that it is not 5 mm or 1 br because these conditions are perfectly separable in Ay-axis.


## 4.4.kNN

A simple classification technique was used to prove the usefulness of the analysis.k-nearest-neighbors (kNN) algorithm was used. Also tests adding white noise were made; the level of noise was based on the signal-to-noise ratio (SNR). The values used were SNR = 10, 5, and 1. Table 2 shows the results of the classification for the Ax-axis with k=1 and for the three cases of added noise. These results are obtained from a cross-validation process with 10-folds; the average loss for each cross-validation model is used to obtain the accuracy level in classification. We have to remark that the simple value of k=1 was sufficient to obtain a desirable accuracy level in classification.Table 2
Classification percentage for each condition usingkNN.
Condition/axis
A
x
SNR10
SNR5
SNR1
Unload
96.55%
96.70%
92.85%
89.8%
Half load
80.5%
81.5%
80.85%
80.55%
Full load
87.6%
85.85%
84.4%
85.15%


## 5. Conclusions and Future Work

In this paper a study of incipient broken bar detection based on vibration analysis is presented. The obtained results show that information regarding partial damage in broken rotor bar is reflected on vibration signals. An initial damage of 3 mm was detected using the autocorrelation and DWT. A simplekNN classification considering the healthy motor and the four levels of damage (3 mm, 5 mm, 7 mm, and 1 br) showed results of 96.55% for unload conditions and 80.5% for half-load and 87.6% for full-load condition.Future work aims to employ different types of classifiers like the Multilayer Perceptron (MLP), Support Vector Machines (SVM), and Artificial Neural Networks (ANN), in order to improve the accuracy level in classification especially in the half- and full-load motor conditions.



---

*Source: 3530464-2016-02-03.xml*