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.
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: Section 2 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. Figure 1 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.
DWT decomposition level frequency bands for a sampling frequency fs = 3200 Hz.
Level (L)
Approximation bandwidth (Hz)
Detail bandwidth (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]
Mallat algorithm for DWT decomposition.
Figure 2 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.
DWT decomposition, (a) approximation, and (b) detail.
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 of n 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.
Autocorrelation and DWT decomposition for a healthy motor, (a) approximation, and (b) detail.
An example of a DWT decomposition of a damage induction motor vibration signal is shown in Figure 4; 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.
Autocorrelation and DWT decomposition for a motor with one broken bar, (a) approximation, and (b) detail.
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.
A common reference among the classifiers due to the simplicity of the learning rule is the k-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:
Vibration signals of the three axes are sensed online from the motor and saved in a PC.
Feature extraction stage is as follows:
An autocorrelation function is applied to each axis.
The decomposition at level three of each signal is obtained using DWT algorithm.
Statistical parameters like mean, standard deviation, and zero-crossings are calculated for the decomposition coefficients.
Statistical analysis of the parameter is computed.
Classification stage is as follows: the motor condition is classified based on kNN 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. Figure 5 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.
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 Figures 6–17.
Approximation coefficients at no mechanical load, Ax-axis.
Unload condition analysis, Az-axis.
Unload condition analysis, Ay-axis.
Unload condition analysis, Ax-axis.
Approximation coefficients at half-mechanical load, Ax-axis.
Half-load condition analysis, Az-axis.
Half-load condition analysis, Ay-axis.
Half-load condition analysis, Ax-axis.
Approximation coefficients at full mechanical load, Ax-axis.
Full-load condition analysis, Az-axis.
Full-load condition analysis, Ay-axis.
Full-load condition analysis, Ax-axis.
4.1. Without Mechanical Load
Figure 6 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 Figure 7, 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.
The Ay-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, Figure 9 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
Figure 10 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.
Figure 12 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.
Figure 13 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
Figure 14 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.
Figure 15 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.
Figure 16 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.
Figure 17 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.
Classification percentage for each condition using kNN.
Condition/axis
Ax
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 simple kNN 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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
YanR.GaoR. X.ChenX.Wavelets for fault diagnosis of rotary machines: a review with applications201496part A11510.1016/j.sigpro.2013.04.0152-s2.0-84887433963Rangel-MagdalenoJ. D. J.Romero-TroncosoR. D. J.Osornio-RiosR. A.Cabal-YepezE.Contreras-MedinaL. M.Novel methodology for online half-broken-bar detection on induction motors20095851690169810.1109/TIM.2009.20129322-s2.0-67349131568GyftakisK. N.SpyropoulosD. V.KappatouJ. C.MitronikasE. D.A novel approach for broken bar fault diagnosis in induction motors through torque monitoring201328226727710.1109/tec.2013.22406832-s2.0-84877920650ZhangP.DuY.HabetlerT. G.LuB.A survey of condition monitoring and protection methods for medium-voltage induction motors2011471344610.1109/TIA.2010.20908392-s2.0-78751666101RazikH.de Rossiter CorrêaM. B.da SilvaE. R. C.A novel monitoring of load level and broken bar fault severity applied to squirrel-cage induction motors using a genetic algorithm200956114615462610.1109/TIE.2009.20295802-s2.0-70350426580Ordaz-MorenoA.Romero-TroncosoR. D. J.Vite-FriasJ. A.Rivera-GillenJ. R.Garcia-PerezA.Automatic online diagnosis algorithm for broken-bar detection on induction motors based on discrete wavelet transform for FPGA implementation20085552193220210.1109/TIE.2008.9186132-s2.0-44349192213NemecM.DrobničK.NedeljkovićD.FišerR.AmbrožičV.Detection of broken bars in induction motor through the analysis of supply voltage modulation20105782879288810.1109/tie.2009.20359912-s2.0-77954568356LiD. Z.WangW.IsmailF.A spectrum synch technique for induction motor health condition monitoring20153041348135510.1109/tec.2015.2454440FaizJ.GhorbanianV.EbrahimiB. M.EMD-based analysis of industrial induction motors with broken rotor bars for identification of operating point at different supply modes201410295796610.1109/tii.2013.22899412-s2.0-84900541494Sapena-BañoA.Pineda-SanchezM.Puche-PanaderoR.Martinez-RomanJ.MaticD.Fault diagnosis of rotating electrical machines in transient regime using a single stator current’s FFT201564113137314610.1109/tim.2015.2444240SadeghianA.YeZ.WuB.Online detection of broken rotor bars in induction motors by wavelet packet decomposition and artificial neural networks20095872253226310.1109/TIM.2009.20137432-s2.0-67649604414SuH.ChongK. T.Ravi KumarR.Vibration signal analysis for electrical fault detection of induction machine using neural networks201120218319410.1007/s00521-010-0512-32-s2.0-79951769436GritliY.Di TommasoA. O.MiceliR.FilippettiF.RossiC.Closed-loop bandwidth impact on MVSA for rotor broken bar diagnosis in IRFOC double squirrel cage induction motor drivesProceedings of the International Conference on Clean Electrical Power (ICCEP '13)June 2013Alghero, ItalyIEEE52953410.1109/iccep.2013.65869042-s2.0-84884634861DengL.ZhaoR.A vibration analysis method based on hybrid techniques and its application to rotating machinery20134693671368210.1016/j.measurement.2013.07.0142-s2.0-84880984970Delgado-ArredondoP. A.Garcia-PerezA.Morinigo-SoteloD.Osornio-RiosR. A.Avina-CervantesJ. G.Rostro-GonzalezH.Romero-TroncosoR. d.Comparative study of time-frequency decomposition techniques for fault detection in induction motors using vibration analysis during startup transient201520151470803410.1155/2015/708034BettaG.LiguoriC.PaolilloA.PietrosantoA.A DSP-based FFT-analyzer for the fault diagnosis of rotating machine based on vibration analysis20025161316132110.1109/TIM.2002.8079872-s2.0-4043048955SadoughiA.EbrahimiM.MoalemM.SadriS.Intelligent diagnosis of broken bars in induction motors based on new features in vibration spectrumProceedings of the IEEE International Symposium on Diagnostics for Electric Machines, Power Electronics and Drives (SDEMPED '07)September 2007Cracow, PolandIEEE10611110.1109/demped.2007.43930792-s2.0-51149099760BaccariniL. M. R.Rocha e SilvaV. V.de MenezesB. R.CaminhasW. M.SVM practical industrial application for mechanical faults diagnostic20113866980698410.1016/j.eswa.2010.12.0172-s2.0-79951580103SadoughiA.EbrahimiM.RezaeiE.A new approach for induction motor broken bar diagnosis by using vibration spectrumProceedings of the SICE-ICASE International Joint ConferenceOctober 2006Busan, Republic of Korea4715472010.1109/sice.2006.3147352-s2.0-34250782126Antonino-DaviuJ. A.Pons-LlinaresJ.Climente-AlarconV.RazikH.Evaluation of startup-based rotor fault severity indicators under different starting methodsProceedings of the 40th Annual Conference of the IEEE Industrial Electronics Society (IECON '14)October-November 2014Dallas, Tex, USA3361336610.1109/iecon.2014.7048995SuH.ChongK. T.Induction machine condition monitoring using neural network modeling200754124124910.1109/TIE.2006.8887862-s2.0-33947363928RodriguezP. J.BelahcenA.ArkkioA.Signatures of electrical faults in the force distribution and vibration pattern of induction motors2006153452352910.1049/ip-epa:200502532-s2.0-33745727602FilhoP. C. M. L.BritoJ. N.SilvaV. A. D.PederivaR.Detection of electrical faults in induction motors using vibration analysis201319436438010.1108/JQME-06-2013-00402-s2.0-84885143456MiceliR.GritliY.di TommasoA.FilippettiF.RossiC.Vibration signature analysis for monitoring rotor broken bar in double squirrel cage induction motors based on wavelet analysis201433516251641Climente-AlarconV.Antonino-DaviuJ. A.Vedreño-SantosF.Puche-PanaderoR.Vibration transient detection of broken rotor bars by PSH sidebands20134962576258210.1109/TIA.2013.22658722-s2.0-84890101487YanR.GaoR. X.ChenX.Wavelets for fault diagnosis of rotary machines: a review with applications20149611510.1016/j.sigpro.2013.04.0152-s2.0-84887433963Garcia-PerezA.Romero-TroncosoR. J.Cabal-YepezE.Osornio-RiosR. A.Rangel-MagdalenoJ. D. J.MirandaH.Startup current analysis of incipient broken rotor bar in induction motors using high-resolution spectral analysisProceedings of the IEEE International Symposium on Diagnostics for Electric Machines, Power Electronics & Drives (SDEMPED '11)September 2011Bologna, ItalyIEEE65766310.1109/demped.2011.60636942-s2.0-81255138786MallatS. G.Theory for multiresolution signal decomposition: the wavelet representation198911767469310.1109/34.1924632-s2.0-0024700097SametH.GhanbariT.AhmadiM.An auto-correlation function based technique for discrimination of internal fault and magnetizing inrush current in power transformers201543439941110.1080/15325008.2014.963265GhanbariT.Autocorrelation function-based technique for stator turn-fault detection of induction motor201510.1049/iet-smt.2015.0118HastieT.FriedmanJ.TibshiraniR.2001New York, NY, USASpringerSpringer Series in Statistics10.1007/978-0-387-21606-5MR1851606