Diagnosis of Defective Rotor Bars in Induction Motors

: This paper proposes a diagnosis method, combining signal analysis and classification models, to the rotor defect problems of motors. Two manufacture technologies, nonmagnetic high-temperature resistant ceramic adhesive and electrical discharge machining (EDM), are applied to make testing samples, including blowhole and perforation defects of rotor bars in this study. The typical multiresolution analysis (MRA) model is used to analyze acquired source current signals of motors. The features are extracted from the signals of each column of MRA-matrix, including maximum, mean, standard deviation, root-mean-square, and summation. The typical back-propagation neural network (BPNN) model is used to diagnose the rotor bar defects of motors, and then the various signal-to-noise ratio (SNR) of white Gaussian noise (WGN), 30, 25, and 20 dB, are added to the signals to verify the robustness of the proposed method. The results show the availability of the proposed method to diagnose the rotor bar defects of motors.


Introduction
Three-phase squirrel-cage induction motors (SCIM) are used widely for industrial applications. Additionally, they are the important equipment for the whole production line. Quality control (QC) must be built in the process of the motor manufacture. Faults in an induction motor are contributed by bearing, stator, rotor, and others, with percentages of 41%, 37%, 10%, and 12%, respectively [1]. In [2] the best characteristics of the infinite impulse response (IIR) filter are exploited to observe the broken rotor bar (BRB) frequencies with good legibility in current and voltage the spectrum of the grid and inverter-fed motor, respectively. Rotor bars are one of the critical components in motors, and defective rotor bars always seriously decay the operation performance, thus the detection of defective rotor bars has primarily been considered. Blowholes of rotor bars are always detected in the die casting process of aluminum alloy when parameters of die casting machines are improper, as shown Figure 1. When a motor has been assembled and quality tested completely, the motor manufacture cuts the rotor destructively to observe the inside of the cross section of defective rotor bars. Each kind of fault always leaves certain specific frequency components in the current or voltage spectrum. It can be used as a potential indicator of fault nature and its severity [3]. There are some methods for detecting defective rotor bars. In [4][5][6][7][8], the literature use a sensing coil, such as Hall transducer or Rogowski coil, to measure the variations of magnetic flux in induction motors. A few studies are based on the operating performances of broken rotor bars, such as voltage [9,10], current [11][12][13], vibration [14,15], noise [16,17], temperature [18], and electromagnetic torque [19]. The aforementioned methods are proposed for the rotor bars malfunction where the induced voltages and currents in the rotor bars are incorrect [20,21]. The voltages and currents cannot be induced to the situation of the broken rotors defect. It is difficult to diagnose the blowholes defect because the inductive voltages and currents of the motors can be easily detected to the situation of the broken rotors defect, but not to the situation of the blowholes rotor bars defect. A good time to diagnose the defective rotor bars is in the QC process, but motor manufacturers usually use a destruct approach or an infrared thermometer [22] to diagnose those. Those detecting methods are in addition to use of voltage, current, and magnetic field, others are not applied in the QC process or need a more expensive detecting sensor, and although the magnetic field sensor is cheap [4], but the sensing coil needs to adjust the position to coordinate with a motor, thus increasing the measurement inconvenience. Nowadays, the common methods of analyzing signals include the fast Fourier transform (FFT) applied as stable and linear, expressed in frequency and energy. The multiresolution analysis (MRA) is applied as nonstable and linear expressed in frequency, energy, and time. The Hilbert-Huang transform (HHT) is applied as nonstable and nonlinear. In the past, FFT [23,24] was used as the main method for the signal analysis of motors. FFT is commonly applied in industry, the spectrum analyzer transforms original signals from the time domain to the frequency domain, in order to obtain the real-time spectrum, but the drawback is inability to get the occurrence time of the signal spectrum. In recent years, signal analysis methods that can resist environmental noise interference have been developed, such as MRA [25][26][27] and HHT. This paper proposes two manufacture technologies, nonmagnetic high-temperature resistant ceramic adhesive and electrical discharge machining (EDM), that are applied to make testing samples to simulate the rotor defects problem, and then acquire the signals of defective rotor bars. Secondly, this study plans an experiment for measuring source current signals of motors and a servo motor as load, the acquired signals, respectively, use typical FFT, MRA, and HHT models to analyze. The features are respectively extracted from the signals of FFT, MRA, and HHT matrix. The normalized features, respectively, use back-propagation neural network (BPNN) and probability neural network (PNN) as classification algorithm models to recognize whether rotor bars are healthy or not and the accuracy of the proposed method is verified by cross validation. Finally, the various signal-to-noise ratio (SNR) of white Gaussian noise (WGN), 30, 25, and 20 dB, are, respectively, added to the original signals to verify the robustness of the proposed method. Finally, the results show the availability of the proposed method to diagnose the rotor bars defects of motors.

Common Signals Analysis Methods
The mother wavelet of MRA is added to the original signals to scale and translate for observing the local feature, but the signals analysis distortion can be caused by the different basis (wavelet family). HHT is a time-frequency analysis tool for analyzing the nonstationary and nonlinear data, it is unique and different from typical FFT and MRA, that can availably reduce the adverse effect of noise interference on signal analysis [28][29][30][31]. MRA uses the scaling function (low-pass filter) and the mother wavelet function (high-pass filter) of the same analytical layer to decompose the original signal into different frequencies in a hierarchical manner, that is, to analyze or reconstruct the signal in different analytical layers. This decomposes the original signal into a monotonous signal for analyzing the information contained in the original signal. HHT is used to analyze nonstable and nonlinear signals. The filtering process of empirical mode decomposition (EMD) is mainly used to eliminate the signal carrier and make the waveform more symmetrical. The signal to be analyzed is decomposed to obtain the intrinsic mode function (IMF) and trends of the signal. This section will introduce the above three signal analysis methods in order.

FFT
FFT is a fast algorithm of discrete Fourier transform (DFT) that can obtain the real-time signals spectrum. A given periodic function can be represented as Fourier series that combines cosine and sine functions [32], is defined as (1).
Fourier transform (FT) applies the Fourier series to a nonperiodicity function, ( ) is called the continuous-time FT of ( ), defined as (5), and ( ) is called the discrete-time FT, defined as (6), where is time index of discrete-time signal, k = 0, 1, …, k-1, is frequency index of DFT, n = 0, 1, …, n-1, is separation of sampling time, is frequency separation, is number of sampling.
As operational factor ⋅ of DFT has the periodicity, therefore, FFT uses Cooley-Tukey type algorithms to make a simple computation. If a number of sampling is and the complexity of computation is , but the complexity of time is that computes by DFT, the former computation complexity is / of the latter. If there is a large number of samples, then more and more computation time is saved [33]. Figure 2 shows the butterfly diagram of 4-element Cooley-Tukey type algorithms that output ( ), defined as (7).
where operational factor = , output (0) = element + ⋅ element , Many motors defective diagnoses are analyzed by FFT, that is, motor current signature analysis (MCSA) [24], the advantage can only measure a lead wire of the balanced three-phase system. When rotors produce asymmetry, and then cause additional side band harmonic components at the fundamental frequency, this harmonic frequency is defined as (8).

MRA
Wavelet transform (WT) is a multiresolution analysis tool which extends from FT, also called as the mathematical microscope, including the continuous wavelet transform (CWT) and the discrete wavelet transform (DWT) [34]. CWT of any function ( ) ∈ ( ) is defined as (9).
where ( ) is mother wavelet function, is scaling parameter, is translation parameter, √ is normalization factor, to ensure energy of and are constant.
DWT can change the scaling-translation parameter of CWT to reduce the complexity computation of DWT, and then uses the scaling function and the mother wavelet as a low-pass filter and a high-pass filter, respectively, to decompose signals in different frequency, that is MRA proposed by Mallat [35,36]. MRA can layer-by-layer decompose signals to observe the implicit message. Formula (10) is MRA of time-domain signal function ( ), ( ) is the first layer scaling parameter after WT, and then uses it as the input signal of the next layer to proceed to the second layer WT, and so on.
is the layer scaling parameter, ( ) is the layer mother wavelet parameter, is sequence data of signal, is MRA layer, = 1, 2, 3 ⋯ .
In Formula (10), ( ) and ( ) are defined as (11) and (12), respectively, ℎ is the coefficient of low-pass filtering, and is the coefficient of high-pass filtering; ( ) and ( ) are composed of linear basis function √2 (2 − ). Figure 3 shows MRA frame of signals, ↓ 2 represents sampling rate of bandwidth halved.

HHT
HHT was proposed by Dr. Norden E. Huang, it is based on the mathematical theory of mathematician Hilbert in the 19th century that was designed to analyze the nonstationary and nonlinear signals. HHT consists of empirical mode decomposition (EMD) and Hilbert spectral analysis (HSA) [37], EMD uses the sifting process to obtain the intrinsic mode function (IMF) and the trend of signal, the sifting process can eliminate a carrier wave and make a more symmetrical waveform, the IMF must be adequate for the following conditions: (1) The number of local extreme points and zero-crossing points of signal shall be equal or one difference at most. (2) Average value of the high and low envelope curves must be close to zero.
The sifting process steps of IMF are as below: Step 1 Find out the local maximum ℎ( ) and the local minimum ( ) input signal ( ).
Step 2 Connect the maximum and minimum envelope curve, respectively, and obtain the mean envelope curve ( ).

Step 4
Whether ( ) is IMF or not? If it is true, then save it into ( ), otherwise, repeat the previous steps 1 to 4.

Step 6
If trend function ( ) is monotone function or constant, then the decomposition is complete, otherwise come back to the beginning step to re-sifting.
The EMD process shown in Figure 4. The original signal decomposes as IMF components, and then each IMF component uses Formula (13) to do Hilbert transform (HT). ( ) is IMF and is the mainly value of the singular integral. The converted data use Formula (14) to represent Hilbert spectral.
( ) and ( ) are the layer intrinsic mode against time-amplitude-frequency functions. Hilbert spectral ( , ) integrates against time.

Neural Network Classifier
Neural network (NN) has been successfully applied in different domains, such as recognizing the faulty electrical machinery [38], monitoring the power quality [39], and classifying the fault signal [40]. This paper proposes a recognition approach for the defect diagnosis of motors, signal analysis methods decompose original signals, automatic recognition of defective rotor bars use acquired signal features by the classifier. The classifier uses the different classification algorithms, such as a multilayer perceptron [41,42], BPNN or PNN, to deal with different kinds of problems. BPNN needs setting the number of neurons and the learning rate of network, this supervised neural network has a high accuracy more than an unsupervised neural network. However, there are also limitations in BPNN. This network requires a lot of input and target for training the network. PNN does not need trained, but needs to find out a smoothing parameter. There are some advantages of using PNN for classification. For example, the run time of PNN is faster than BPNN, and it is robust to noise. Furthermore, the training manner of PNN is simple and instantaneous.

BPNN
BPNN is a multilayer feedforward neural network [41,42], also known as a multilayer perceptron which has a fast self-learning speed and high accuracy features for application to classification and predication. The network architecture shown in Figure 5, including input layer, hidden layer, and output layer. , , . . . , are the input of network, , , , are the predicted output values of network, neurons interconnect with each other by a weight, but they do not connect with the same layer neurons. The hidden layer can be a single layer or a multilayer architecture which depends on the complexity degree of functions, in most of the cases, as long as the number of neurons are enough, a single hidden layer can conduct any function.
Signals are feedforward propagation (feedforward phase) in BPNN, the predication of output error are back propagation (back propagation phase). Input vectors sequentially propagate to the output layer via the input layer, the hidden layer of network in the feedforward phase. The weight value of the network is corrected by the error of output and the set target value in the back propagation phase, the purpose is forcing the output value close the expected target value. The basic principle of BPNN uses the gradient descent method [43] which continuously corrects the weight and the bias to promote the output result close to the expected value. NN application process is as follows: firstly, the signal data are obtained from the experiment, and then divided into training and testing data, training data are used to train the network to set up network model, testing data are used to test the accuracy of the network. When the recognizable result is not ideal, in addition to training the network again, the training data can be increased, it would be beneficial to reach the expected value of network output. Besides, if we want to build a successful network model that shall confirm whether the input and output of network exist in relationship, the algorithm steps are as below: Step 1 Set the number of neurons of input layer, hidden layer, and output layer.
Step 2 Set the initial value of network, randomly produce the weight and the bias , and then set the trainable learning rate which usually use 0.1-0.9, determine which type of activation function, as represented is sigmoid function and is piecewise linear function [43]. Step 3 Set the input vector as , the output of hidden layer and the output layer in the feedforward phase, defined as (15) and (16), respectively.
Step 4 Calculate the error function ( ) = ( ) − ( ) which is the number of training times, is the target value, is the output value.
Step 5 Calculate the gradient of weight and bias, defined as (17) to (20), respectively.
Step 7 Repeat steps 1 to 6 to reach the number of training times or compliance with the convergent criteria.

PNN
PNN was proposed by Dr. D. F. Specht, it replaces the commonly used sigmoid activation function by exponential function, and then uses Bayes optimal decision rule to calculate the nonlinear decision boundary. After that, PNN uses linear learning algorithm to improve the Bayesian classifier's weakness that it is not easy to create a probability density function in each class [44].
PPN belongs to a forward, supervised, and learnable network architecture which has both fast training and learning speed, it uses the training data as the weight of neurons with each other, and then needs only to adjust the smoothing parameter of probability density function (PDF), to reduce the error of the network output and the target vector value with each other that let the output vector close to the target vector value. PDF ( ) defined as (25), shows that if the input vector of classification is different from the training vector data of category , and the smoothing parameter does not adjust appropriately, then prompts PDF ( ) close to zero, which caused a poor recognition, therefore, to adjust the smoothing parameter is the key for reaching the optimal NN. The basic architecture of PNN shown in Figure 6. The first layer is the input layer for classifying the eigen vector , , . . . , and then sends them to the second layer. The second layer is the feature layer or the pattern units, it undertakes the input vector of the above layer, and then proceeds to the nonlinear calculation by the exponential activation function. The third layer is the summation layer which is the sum of the different types of feature units to send to next layer. The fourth layer is the output layer or the competitive layer where the healthy rotor bar competes with three kinds of defective rotor bars each, and then obtains the recognizable result, that is called G.
The four-layer neural network architecture can handle classification problems of any input type and can modify the decision boundary range through new input signal data in real-time. Therefore, the speed of the parallel operation is relatively better than that of the BPNN that requires continuous weight adjustment. Additionally, it is more suitable for the conversion of the neuron system into hardware architecture for application [45].
where ( ) is PDF of category at point, is the number of modes, is dimension of measure space, is smoothing parameter, is amount of category training vector, is input classificatory vector, is training vector data of category , is transposed vector.

Experimental Setup for Defective Rotor Bars
This section discusses how to manufacture rotor bar defects, blowhole defects, perforation defects-EDM, experiment setup, and measurement data.

Manufacture of Rotor Bar Defects
This study mainly addresses defects of rotor bars as the average of rotor fault rate of induction motors in industrial applications was approximately 10% [1,[45][46][47]. The defective rotor bars are divided into two parts: 1. Manufacture of blowhole defect. 2. Manufacture of perforation defect-EDM. Figure 7 depicts the blowhole defect of rotor bars. When the silicon steel of the rotor is laminated completely and then the ceramic adhesive is fixed into the rotor slot so that it becomes the blowhole defect. Moreover, Figure 8 shows the rotor die-cast aluminum alloys by the die-casting machine which becomes an aluminum-rotor, and then the rotor cuts by a cutting machine. The result shows that the availability of the proposed manufacture technology to simulate the blowhole defect of rotor bars as shown in Figure 9. Therefore, besides making one healthy 2HP SCIM, two rotors in which one has a blowhole defect (in one rotor bar), and the other one has four blowhole defects (in four rotor bars) are made also. Moreover, a computer numerical control (CNC) is used to make an EDM process to simulate a perforation defect of the rotor bar (in one rotor bar). Thus, this experiment has four testing samples: one healthy SCIM, one rotor with a blowhole defect, one rotor with four blowhole defects, and one rotor with a perforation defect.

Manufacture of Blowhole Defect
A measuring cup measures 1535 °C resistant ceramic powder and clean water as shown in Figure  10a and 10b, using agitation mixing by the ratios of 100:15 to become ceramic adhesive (Figure 10c), and then measures 20 mm in length of plastic pipe as a mold (Figure 10d), ceramic adhesive injected into the plastic pipe to be cut to take out ceramic adhesive after it hardens, that is a simulated blowhole (Figure 10e). The simulated blowhole was put into the rotor slot (Figure 10f), the diameter of the rotor is 95 mm, one blowhole in a slot ( Figure 10g) and four blowholes in four slots ( Figure  10i), the dimension of the rotor slot shown in Figure 10h and 10j. Finally, rotor die-cast aluminum alloys by the die casting machine become the aluminum-rotor finished product with 72 mm core length ( Figure 11). 72mm Figure 11. Aluminum rotor.

Manufacture of Perforation Defect-EDM
A rotor bar of one healthy rotor proceeds EDM by CNC to simulate a perforation defect, the CNC control panel is shown in Figure 12a. Firstly, the CNC set 9 mm in vertical depth (Figure 12b), the EDM platform is shown in Figure 12c, EDM processing shown in Figure 12d, discharge electrode shown in Figure 12e, the finished product of EDM has one perforation defect that is 2 mm in diameter and 9 mm in vertical depth as shown in Figure 13.

Experiment Setup and Measurement Data
This paper proposes two defective rotor bars which are blowhole and perforation. Additionally, the signal analysis proceeds with FFT, MRA, and HHT, respectively, for understanding the abnormal conditions and the serious degree of the defective rotor bars. Experiment setup dynamometer platform combines one servo motor (3 phase AC 69 Hz, 11 kW, 2000 rpm, Volts: 190-330, Amps: 41.2-23.8) with the torque sensor as load of IM (3 phase SCIM 4 Pole 1.49 kW, 60 Hz, Volts: 220, Amps: 5.58, 1715 rpm). Any one lead wire of three phase IM connects to (NI)PIX-1033 data acquisition, and adjusts the motor speed and torque by the control panel of the dynamometer platform, and then a current probe measures the single-phase motor current that the current signal records by the data acquisition, and sampling frequency is 1000 Hz/s. The experimental system architecture is shown in Figure 14 and equipment of the experiment are shown in Figure 15.

Features Extraction and Defects Recognition
This study uses (NI) PXI-1033 current acquisition of National Instruments to set the sampling rate per second equal to 1000 Hz, and takes the current signals of four IM samples which run, respectively, for 100 s. The measured current signals are divided into 100 records by Matlab, in example four samples amount to 400 records, and then the signal analysis proceeds with FFT, MRA, and HHT, respectively, BPNN and PNN classification algorithm as the defective recognition classifier, but after transform analysis the data is too much to compute directly by classifier which caused the high computation time. Therefore, the six features are extracted from the signals in advance, including the maximum values, the minimum values, the root mean square values, the average values, the standard deviation values, and the summation values of the raw matrix. The features are normalized between 0 and 1, 80 of the 100 recorders are the training data, and the remaining 20 recorders are the testing data to input in classifier, finally, the accuracy of the proposed method is verified by 10 times cross validation. In addition, the various SNR of WGN, 30, 25, and 20 dB, are, respectively, added to the signals to verify the robustness of the proposed method.

FFT-Based Recognition
Firstly, the healthy rotor and three samples of defective rotors (with one blowhole rotor bar, four blowhole rotor bars, and one perforation defective rotor bar) take 100 recorders each sample, four samples amount to 400 records, and then the signal analysis proceeds with FFT to obtain the frequency-amplitude spectrum. Secondly, the features are extracted from the signals, including maximum, minimum, mean, standard deviation, root-mean-square, and summation values of the FFT matrix, and the features are normalized between 0 and 1. The features energy diagram of 100% loading rate and no noise is shown in Figure 16a, where the vertical axis is the number of features and the horizontal axis is the sample number. The enlarge diagram of features 4 and 5 are shown in Figure 16b and 16c, which can recognize clearly whether the rotor bars are healthy or defective (blowhole or perforation), but it is hard to distinguish between one blowhole and four blowholes. In addition, Figure 17a is shown the features energy diagram of 100% loading rate and 30 dB noise is added to the signals, where the vertical axis is the number of features and the horizontal axis is the sample number. The enlarge diagram of features 4 and 5 are shown in Figure 17b and 17c, which also can recognize clearly whether the rotor bars are healthy or defective (blowhole or perforation), but it is still hard to distinguish between one blowhole and four blowholes.

MRA-Based Recognition
The Haar wavelet decomposes the fault signals (stator currents) into one to five levels of highfrequency components, moreover, five signals features are extracted from the high-frequency signals, including maximum (Max), mean (Mean), standard deviation (Std), root-mean-square (Rms), and summation (Sum) values of the raw matrix, five layer high-frequency signals amount to 25 features, and the features extraction process is shown in Figure 18. The original signals are decomposed into the low-pass filtered (ℎ ) and the high-pass filtered ( ) coefficients, and extracts five features, including Sum (F1), Std (F2), Rms (F3), Mean (F4), and Max (F5), and the low-pass filtered (ℎ ) signals to do the second layer decomposition to obtain the second layer low-pass filtered (ℎ2) and high-pass filtered ( 2) coefficients. Additionally, five features are extracted, including Sum (F6), Std (F7), Rms (F8), Mean (F9), and Max (F10), and so on, after the fifth layer is decomposed to obtain F21 to F25 features. Finally the 25 features are input into the neural network classifier for training and cross validation.   The features energy diagram of 100% loading rate and no noise is shown in Figure 19a, where the vertical axis is the number of features and the horizontal axis is the sample number. The enlarged diagrams of features 8 and 9 are shown in Figure 19b and 19c, which can recognize clearly whether the rotor bars are healthy or defective (blowhole or perforation), but it is hard to distinguish between one blowhole and four blowholes. In addition, Figure 20a shows the features energy diagram of 100% loading rate and 30 dB noise is added to the signals, where the vertical axis is the number of features and the horizontal axis is the sample number. The enlarged diagrams of features 8 and 9 are shown in Figure 20b and 20c, which also can recognize clearly whether the rotor bars are healthy or defective (blowhole or perforation), but it is still hard to distinguish between one blowhole and four blowholes.

HHT-Based Recognition
EMD method decomposes the original signal (stator current) into five layers of IMF, a HHT analysis of each extracted IMF to obtain the Hilbert matrix with time-energy (Time scale) and timefrequency (frequency scale), respectively, to find the maximum characteristic curve, average characteristic curve, root-mean-square characteristic curve, standard deviation characteristic curve, energy sum characteristic curve. Additionally, five signals features are extracted, including Max (Tmax, F-max), Mean (T-mean, F-mean), Std (T-std, F-std), Rms (T-rms, F-rms), and Sum (T-e, F-e) values of the raw matrix which amount to 50 features, 10 characteristic curves multiplied by five features (F1, F2, …, F50). The features extraction process of HHT is shown in Figure 21. The features energy diagram of 100% loading rate and no noise is shown in Figure 22a, where the vertical axis is the number of features and the horizontal axis is the sample number. The enlarged diagrams of features 32 and 33 are shown in Figure 22b and 22c, which can recognize clearly whether the rotor bars are healthy or defective (blowhole or perforation), but it is hard to distinguish between one blowhole and four blowholes. In addition, Figure 23a shows the features energy diagram of 100% loading rate and 30 dB noise is added to the signals, where the vertical axis is the number of features and the horizontal axis is the sample number. The enlarged diagrams of features 32 and 33 are shown in Figure 23b and 23c, which also can recognize clearly whether the rotor bars are healthy or defective (blowhole or perforation), but it is still hard to distinguish between one blowhole and four blowholes.  The various SNR of WGN, 30, 25, and 20 dB are added to the signals, three kinds of signal analysis methods are used to analyze the signal, and obtain features, including six FFT-based features, 25 MRA-based features, and 50 HHT-based features, and then PNN and BPNN, respectively, used as classifiers for classification and recognition, the recognition rate of PNN classification is shown in Table 1. The recognition rate of FFT is lower than MRA and HHT which have more feature numbers. The recognition rate of HHT is 94.8% under ∞dB noise and 84.5% under 25 dB noise, but it decreases to 68.4% under 20 dB noise which shown lack of anti-noise capability. The recognition rate of MRA is similar to HHT under low noise, but 74.6% is higher than HHT under 20 dB noise. The trend of recognition rate toward PNN and noise interference is shown in Figure 24. From ∞ dB to 25 dB noise, the trend of anti-noise recognition rate is still appearing as a linear trend on FFT. It is similar on HHT and MRA, both of two method show a linear trend from ∞ dB to 25 dB, but HHT has quite an obvious descent that it is not sufficient on the anti-noise ability from 25 to 20 dB. However, MRA always shows a stationary linear trend. Figure 25 represents the recognition rate toward PNN and the load conditions, FFT has a lower recognition rate than HHT and MRA under each load condition, although the recognition rate of HHT is similar to MRA under light, medium, and heavy load conditions, but MRA has a higher recognition rate under no load condition. Q.C. (Quality Control) is routine test in the manufacturing process of motors.   Table 2 shows the recognition rate of using BPNN classification, it is quite obvious that FFT is lower than MRA and HHT at various degrees of noise, and the recognition rate of MRA approximates to HHT at 25 dB noise, but it is 85.3% at 20 dB noise.  Figure 26 represents the trend of the recognition rate toward BPNN and the noise interference, which the anti-noise tendency of MRA and HHT still show a linear trend from ∞ dB to 25 dB, but HHT has quite an obvious descent that it is not sufficient on the anti-noise ability from 25 to 20 dB. However, MRA always shows a stationary linear trend. Figure 27 represents the recognition rate toward BPNN and the load conditions, HHT has a higher recognition rate under medium load and heavy load conditions, but MRA has a higher recognition rate under light load and no load conditions, it has the advantage of the no load condition especially.

Discussion
This study proposed a diagnosis method, combining the MRA signal analysis and the BPNN classification model, to recognize rotor defective problem of motors. The experimental evidence demonstrated that this proposed method has three advantages below: 1. As mentioned above, a good time to diagnose defective rotor bars is in the QC process for motor manufacture, and this study proposes the method has the high recognition rate of 97.0% under no-load and 97.6% under light-load conditions, respectively, it more advantage compared with another two kinds of signal analysis methods, especially under no-load condition. Under noload condition, MRA recognition rate of 97.0% is higher than the HHT of 87.2%. 2. Suitable for high noise workshops: the recognition rate still has 85.3% under 20 dB high noise interference environment, for example it has a good anti-noise ability. 3. Only a few features can make the defective recognition; MRA recognizes only half of the features of HHT that can effectively reduce the recognition time.

Conclusions
This paper proposed a diagnosis method, combining signal analysis and NN classification algorithms under different loading rate conditions, to rotor defective problems of IM motors. FFT, MRA, and HHT were, respectively, used to analyze the motor current signals where six signals features were extracted from the signals, including maximum, minimum, mean, standard deviation, root-mean-square, and summation values of the raw matrix. The typical PNN and BPNN models were used to recognize the rotor bar defects of motors, and this paper studied the advantages and disadvantages, combining three kinds of signal analysis methods with two kinds of NN, to apply in the rotor bars defects of motors.
Furthermore, this paper considered the anti-noise ability for the environment where the various SNR of WGN, 30, 25, and 20 dB, were added to the signals to verify the robustness of the proposed method.
The results of experiment and Matlab simulation show that the recognition rate of BPNN is higher than PNN, and MRA is has a more linear stationary recognition rate than HHT for anti-noise interference. As mentioned above, this paper recommends that MRA combined with BPNN is the comparatively good method for diagnosing the rotor bar defects of motors.