Next Article in Journal
Robustness of Cloud Manufacturing System Based on Complex Network and Multi-Agent Simulation
Next Article in Special Issue
An Entropy-Based Condition Monitoring Strategy for the Detection and Classification of Wear Levels in Gearboxes
Previous Article in Journal
Event-Triggered Tracking Control for Adaptive Anti-Disturbance Problem in Systems with Multiple Constraints and Unknown Disturbances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Broken Bar Fault Detection Using Taylor–Fourier Filters and Statistical Analysis

by
Sarahi Aguayo-Tapia
1,
Gerardo Avalos-Almazan
1,
Jose de Jesus Rangel-Magdaleno
1,* and
Mario R. A. Paternina
2
1
Digital Systems Group, National Institute for Astrophysics, Optics and Electronics, Puebla 72840, Mexico
2
Department of Electrical Engineering, National Autonomous University of Mexico, Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(1), 44; https://doi.org/10.3390/e25010044
Submission received: 14 November 2022 / Revised: 14 December 2022 / Accepted: 24 December 2022 / Published: 27 December 2022

Abstract

:
Broken rotor bars in induction motors make up one of the typical fault types that are challenging to detect. This type of damage can provoke adverse effects on the motors, such as mechanical and electrical stresses, together with an increase in electricity consumption, causing higher operative costs and losses related to the maintenance times or even the motor replacement if the damage has led to a complete failure. To prevent such situations, diverse signal processing algorithms have been applied to incipient fault detection, using different variables to analyze, such as vibrations, current, or flux. To counteract the broken rotor bar damage, this paper focuses on a motor current signal analysis for early broken bar detection and classification by using the digital Taylor–Fourier transform (DTFT), whose implementation allows fine filtering and amplitude estimation with the final purpose of achieving an incipient fault detection. The detection is based on an analysis of variance followed by a Tukey test of the estimated amplitude. The proposed methodology is implemented in Matlab using the O-splines of the DTFT to reduce the computational load compared with other methods. The analysis is focused on groups of 50-test of current signals corresponding to different damage levels for a motor operating at 50% and 75% of its full load.

1. Introduction

High efficiency, low cost, and low maintenance make up some of the main features of induction motors (IM) that make them the preferred machine type in the industry. For this reason, early fault detection becomes critical to schedule maintenance and prevent further damage to the motor [1].
Various factors can produce a fault in an IM, such as manufacturing defects, corrosion, overloading, and overheating [2]. The most common damages in IMs are related to bearing faults, stator faults, and broken rotor bar faults [3,4]. Among these, the broken bar damages are the most difficult to detect because this condition does not show any effect on the motor at the early stages of the fault.
Suppose the broken bar damage is not conveniently detected and treated. In that case, the unbalanced currents flowing through the rest of the rotor bars lead to thermal and mechanical stresses that may conduct the IM to complete failure. Moreover, this condition affects the operative cost of the machine, considering that this type of fault can significantly increase power consumption.
To achieve a convenient detection of broken bars, diverse techniques have been studied in recent years; some of them are based on vibration signals, some others on current analysis, and magnetic flux, among others [5,6,7,8,9,10,11]. For example, Morales-Perez et al. [5] explored the use of trained dictionaries to decompose vibration signals using the orthogonal matching pursuit algorithm. Panagiotou et al. [6] proposed the detection of broken rotor bars at steady state by analyzing stray flux signals. Garcia-Calva et al. [8] analyzed the speed signature of a faulted motor in the time-frequency domain. Overall, motor current signature analysis (MCSA) is widely used because current signals acquisition does not require additional sensors in the IM.
Previous works have studied the broken bar problem using MCSA, relying on traditional methods such as the fast Fourier transform (FFT) or the wavelet transform. For instance, da Costa et al. [12] used two diagnostic methods based on the FFT and Wavelet Transform, to compare them. García-Perez et al. [13] employed the multiple signal classification (MUSIC) method applied to current and vibration signals to identify multiple faults. Trujillo-Guardado et al. [14] proposed an algorithm for broken bar detection hinged on the use of two Taylor-Kalman filters in cascade.
For classification purposes, most of the recent works use artificial intelligence (AI) methods: Support vector machines (SVM), artificial neural networks (ANN), and convolutional neural networks, among others. For instance, Reljic et al. [15] introduced ANN and SVM to perform broken bar detection. Xiao et al. [16] utilized a recurrence quantification analysis and long short-term memory neural network to detect different faults in the IM. Valtierra-Rodriguez et al. [17] proposed using convolutional neural networks to detect broken rotor bars automatically. Kumar and Hati et al. [18] employed a dilated convolutional neural network model to detect bearing and broken rotor faults based on vibration signals. These techniques have proven high accuracy, but their main drawback is the high computational resources required for their use, which may be difficult to implement for monitoring purposes. Moreover, even with the use of AI techniques, most recent works are able to detect no less than half-broken bars.
Aiming to improve the detection rate without consuming large computational resources or requiring long computational times, the Taylor–Fourier Transform (TFT) emerges as a good option for signal processing based on its excellent filtering characteristics. In general terms, the TFT is an expansion of the Fourier subspace given by the modulation of each Fourier harmonic by a Taylor polynomial of order K greater than zero. The TFT has been mainly applied to power systems monitoring, and oscillatory signals analysis [19,20]. This transform delivers dynamic phasor estimation, which opens the possibility of working with transitory signals. Applying the TFT to a current signal offers the possibility to estimate mainly but not only the amplitude and phase of the signal. On the other hand, the O-splines implementation of the TFT considerably reduces the computation resources required for signal processing.
In this paper, MCSA for broken bar fault classification is implemented using Taylor–Fourier filters and taking advantage of the processing facilities given by the use of the O-splines. The feature extraction is performed based on the amplitude estimation given by the TFT. Finally, an analysis of variance is performed on each distribution to detect the damage using the Tukey procedure.
This paper is organized as follows. In Section 2, fundamentals of the TFT and its O-splines are presented, as well as the basis for MCSA and statistical analysis. Then, the proposed methodology for broken bar fault classification is given in Section 3. Results are presented in Section 4, and finally, conclusions are given in Section 5.

2. Fundamentals

2.1. Digital Taylor–Fourier Transform (DTFT)

The TFT is an expansion of the Fourier Transform. It can be seen as a Finite Impulse Response (FIR) filter bank formed by a set of maximally flat filters with less inter-harmonic interference than a Fourier Transform [21].
The signal model of the DTFT is given by
x ( n ) = h = H H p h ( n ) e j 2 π h f 1 n , C N 2 n C N 2
where n = 0 , 1 , , N 1 are the samples of the signal being N the number of samples per cycle, p h ( n ) = a ( n ) e j ϕ ( n ) represents the dynamic phasor, a ( n ) and ϕ ( n ) are the amplitude and phase respectively, and f 1 stands for the central frequency. C is the number of cycles defined as C = K + 1 , with K being the number of Taylor polynomials used to approach p h ( n ) as:
p h K ( n ) = p h ( n 0 ) + p h ˙ ( n 0 ) ( n n 0 ) + + p h K ( n 0 ) ( n n 0 ) K K !
A least square problem appears with the inclusion of (2) into (1) and leads to a factorized form of a matrix B that includes the Taylor Y and Fourier Ω terms, respectively [19]:
B = Y Ω = I T 1 1 K ! T 1 K I T 2 1 K ! T 2 K I T C 1 K ! T C K W N 0 0 0 W N 0 0 0 W N
where T C K are diagonal submatrices containing the sequence t = [ L h L h + 1 L h ] T / ( N F 1 ) where L h = ( L 1 ) / 2 and L = C N ; W N is a Fourier matrix given by:
W N = 1 1 1 1 1 ω N ω N 2 ω N ( N 1 ) 1 ω N 2 ω N 4 ω N 2 ( N 1 ) 1 ω N ( N 1 ) ω N 2 ( N 1 ) ω N ( N 1 ) 2
where ω N = e j 2 π / N is the harmonic phase factor for the N t h harmonic.
The solution of the least square problem leads to the expression for the analysis equation, which is the core of the filtering algorithm:
ξ ^ = B ˜ H p ,
where ξ ^ encloses the Taylor–Fourier coefficients; the dual matrix, B ˜ , is defined by:
B ˜ = B ( B H B ) 1
Then, considering that B = Y Ω , we have the following expression [19]:
B ˜ = Y ( Y H Y ) 1 Ω N = Y ˜ Ω N
Given that Y is an invertible square matrix, its dual is the transposed inverse matrix [22], this is:
Y ˜ = ( Y 1 ) T = ( C o f ( Y ) T ) T | Y | = C o f ( Y ) | Y |
where C o f ( Y ) and | Y | are the matrix of cofactors and the determinant of the matrix Y , respectively.
Finally, the amplitude estimation α ^ is obtained from the Taylor–Fourier coefficients ξ ^ through:
α ^ = 2 | ξ ^ |

2.2. O-Splines

In this section, the Kth O-spline and its first derivatives are calculated in closed-form based on the expression stated in (8), where the inverse matrix of Y is calculated. Its rows are concatenated in one vector B ˜ 0 K . For amplitude and phase estimation, only the first vector of the dual matrix B ˜ 0 ( K ) is required.
Hence, for K = 0 we need one observation cycle with t 1 = t [ T 1 2 , T 1 2 ) ; in this case B 0 = 1 and consequently B 0 0 ˜ = 1 . The O-spline is a rectangular pulse given by [23]:
B ˜ 0 ( 0 ) ( t ) = 1 f o r T 1 2 t T 1 2 , 0 o t h e r w i s e .
For K = 1 two cycles are required with t 1 = t [ T 1 , 0 ) and t 2 = t [ 0 , T 1 ) = t 1 + T 1 ; then, we have
B o ( 1 ) = 1 t 1 1 t 2
with | B o ( 1 ) | = t 2 t 1 = T 1 . Then:
B ˜ o ( 1 ) = t 2 1 t 1 1 T 1 = u 1 + 1 F 1 ( u 2 1 ) + F 1
where u n represents the normalized time u = t n / T 1 . The columns of the matrix B o ( 1 ) describe a triangular pulse:
B ˜ o ( 1 ) ( u ) = u + 1 f o r 1 u < 0 , 1 u f o r 0 u < 1 , 0 o t h e r w i s e .
This process is repeated with the consecutive values of K; however, it is noteworthy to point out that only the first two polynomials are needed to obtain the rest of the elements of the matrix B ˜ 0 ( K ) by recursive differentiation.

2.3. Motor Current Signature Analysis (MCSA)

The MCSA method has been applied in several works to detect different types of faults in IMs. Morever, fault detection systems are typically based on the online monitoring of stator current signals [3]. In the particular case of broken bars ( f b b ), this type of fault induces a spurious component with a frequency pattern given by
f b b = f 1 ( 1 ± 2 k s )
where f 1 is the supply frequency, k is the number of the harmonic, and s represents the slip of the motor [24]. Considering this feature, the localization of the fault component is known; however, it is important to note that this component is very small compared to the fundamental frequency. Moreover, it lies very close to the supply frequency, making its detection challenging.

2.4. Analysis of Variance and Tukey Procedure

The analysis of variance (ANOVA) allows testing equality among different data distributions by comparing the variance between their distributions relative and the variance within distributions. This test is used in several fields, such as agricultural experiments and medical research. Depending on the type of relationship among groups or distributions, a certain model of ANOVA is applied. For instance, to determine if any differences exist in the means of 3 or more groups, the one-way ANOVA test is the most common method [25].
The first step for applying this method consists in formulating the null hypothesis, which normally is “the population means of the groups are equal” or H 0 : μ 1 = μ 2 = μ 3 , whether the alternative hypothesis would be “at least one of the group means is different” or H 1 : μ 1 μ 2 or μ 1 μ 3 or μ 2 μ 3 . If at least one group differs, the null hypothesis can be rejected [25].
The core of the ANOVA test is the F calculation, which is comprised of the variance ratios [26]:
F = Between groups variance Within groups variance = S S R n ( t ) / ( k 1 ) S S E n ( t ) / ( n k )
where S S R n ( t ) and S S E n ( t ) given by
S S R n ( t ) = i = 1 k n i ( X ¯ i ( t ) X ¯ ( t ) ) 2
and
S S E n ( t ) = i = 1 k J = 1 n i ( X i j ( t ) X ¯ i ( t ) ) 2
represents the sum of squares between groups and within groups, respectively, and X ¯ and X ¯ i ( t ) represent the mean of all groups and the mean of each group, which are expressed as follows:
X ¯ ( t ) = 1 n i = 1 k j = 1 n i X i j ( t )
X ¯ i ( t ) = 1 n i j = 1 n i X i j ( t ) ,
with i = 1 , , k , being k the number of groups and n the number of observations of the group.
A large value of F provides evidence to reject H 0 . Thus, a subsequent method for multiple comparisons can be applied [27]. A typical table of results for an ANOVA test includes the sum of squares, the mean sum of squares, F, and P calculated between groups and within groups.
Usually, ANOVA is followed by another technique to identify the pattern of differences in the distributions. Often, this analysis consists of pairwise comparisons between the means to determine if there is a significant difference. The Tukey test is one popular method to perform this comparison.
This test is based on the computation of Q, which is a criterion that evaluates the difference between the means of two groups [28].
Q = M i M j S S E n ( 1 S )
where M i and M j are the mean of the groups being compared, S S E n is the sum of squares within groups, and S is the number of observations per group. This value is then compared against a Q c r i t i c a l value from a table, which depends on the degrees of freedom, the number of means being tested, and the significance level α [28].

3. Proposed Methodology

The diagram of Figure 1 illustrates the steps that structure the proposed methodology for broken bar detection using the DTFT. First, current signals for different health conditions of the motor need to be acquired. Then, the characteristics of the Taylor–Fourier filters can be selected to perform an amplitude estimation in a frequency of interest. Afterward, a statistical analysis is performed over each group of amplitude estimations. Then the variance analysis is carried out. Hence, discrimination between healthy and damaged rotor bars can be conducted, resulting in the detection of the fault. Next, each process is explained.

3.1. Current Signals

These signals are associated with different damage levels in the rotor bar. Typically, the analysis is focused on one broken bar, half broken bar, and so on; the healthy motor signal also needs to be acquired. These signals do not require any pre-processing and are analyzed according to the sampling frequency at which they were acquired.
In the presence of a level of damage in the bars, spurious components typically appear around the fundamental frequency and its harmonics, as stated by the theory in Section 2.3. According to this expected location of the fault component, the central frequency f c of the filter is selected.

3.2. Filter Design

The design’s parameters need to be selected according to the characteristics of the current signal. The first important parameter is the central frequency f c , which is selected in the estimated location of the fault to prevent any attenuation of the signal.
Another important feature is the bandwidth of the filter f b w , which should be selected in a way that the undesirable energy peaks in the analysis are not considered. For example, a bandwidth narrow enough for analyzing a spurious component related to a broken bar fault will let out the energy of the supply frequency. In this way, the filtered signal will preserve only the features of the faulty component.
Finally, the number of Taylor polynomials (K) is another important design parameter, considering that as this number increases, it is necessary to analyze more windows in the time domain. Still, at the same time, the filter response around the frequency of interest tends to be flatter and, therefore, more likely to be an ideal filter [29]. As can be observed from the normalized power spectral density (PSD) in Figure 2, the filter corresponding to K = 5 presents a flatter response and better harmonic rejection than, for example, the one with K = 1 , Nonetheless, the first one requires more cycles of the signal to implement it.
Once the filter is designed, a sliding window process of size L = C × N is performed on the signal, and a single value is taken at the center of each window, forming the filtered signal. It is important to consider that the filtering results will have a delay corresponding to half of the number of cycles C required for the analysis; remember that this is given by selecting the value of K.

3.3. Amplitude Estimation

After the filtering process, the amplitude estimation can be obtained using (9). A single amplitude value is taken at each processed signal, forming a vector of size M, where M is the number of measurements for a certain case of damage in the bar.

3.4. Statistical Analysis

The variance analysis is performed, and if the null hypothesis is successfully rejected, the Tukey procedure is applied to detect the damage in the rotor’s bar. This can be achieved if the damaged groups differ significantly from the healthy group. An exaggeration of the expected behavior of the groups is depicted in Figure 3, where the broken bar induces components that grow in amplitude according to the severity of the damage; the depicted distributions fulfill the 3 standard deviations ( σ ) criteria of the Chebyshev theorem.

4. Experimental Results

This work analyzed 10 different damage levels for an IM operating under load conditions of 75% and 50%. Each damage level was artificially created by drilling a hole in the bar, as can be observed in Figure 4, starting from 1 mm hole perforation and increasing by 1 mm until reaching a completely broken bar, as illustrated in Figure 5. The experimental test was carried out through a three-phase IM WEG 00118ET3EM143TW, 1HP, 1800RPM, 220 VAC/60 Hz, and 2.98 A; the data was acquired at a sampling frequency of 3200 Hz. For each case of damage, 50 steady-state current measurements were analyzed in this work. For the ANOVA, a significance level of α = 0.01 was selected.
Current signals of the healthy motor and broken bar IM illustrated in Figure 6 show that the effect of this fault is hardly noticeable in the time domain. Thus, the analysis is focused on its frequency spectrum.
Figure 7 depicts the frequency spectra of 50 signals corresponding to a half-broken bar condition. Fault-related components were identified around 54 Hz. Thereby, the central frequency f c of the filter was selected accordingly. Continuing with the parameters for the filter design, a value of K = 5 was selected for all study cases, given that the filter response is flatter enough without requiring larger segments of the signal for the analysis. This selection implies a delay in the results of 3 cycles of the signal sampled at 3200 Hz. Then, the associated O-spline is depicted in Figure 8.
It is worth noticing from Figure 7 that the spurious component represents a very small part of the signal’s energy. Besides, it lies near the center frequency of 60 Hz. For this reason, a bandwidth of 2 Hz was selected for the filters to avoid including other components that are irrelevant to the analysis. To illustrate this, Figure 9 shows the same 2 Hz bandwidth filter with a f c = 54 Hz and the frequency spectrum of a synthetic 60 Hz signal to remark that the reduced bandwidth excludes the energy of the fundamental component.
ANOVA results based on the amplitude estimations of Figure 10 are given in Table 1. The value of p < 0.01 rejects the null hypothesis, allowing us to continue with the Tukey procedure.
On the other hand, the Tukey test results are summarized in Table 2, suggesting that the detection is achieved when the broken bar is greater or equal to 1 mm, with a 99% of confidence level. Remember that the groups significantly differ when Q > Q c r i t i c a l .
For a 50% load condition, the spurious component appears around 56 Hz; then, the filter is centered at this frequency, similar to the one presented in Figure 7. The ANOVA is applied to the amplitude distributions in Figure 11, resulting in Table 3, which allows the rejection of H 0 . Hence, the Tukey test results in an incipient detection starting at 1 mm, as summarized in Table 4.

Comparison Using Fourier filters

To showcase the improvement using Taylor–Fourier filters, the proposed methodology is repeated but using a value of K = 0 , i.e., Fourier filters. As illustrated in Figure 12, the filter exhibits a poor harmonic rejection, giving as a result, the inclusion of undesired energy in the analysis.
Results for the ANOVA based on the distributions of Figure 13 and Figure 14 are summarized in Table 5 and Table 6, respectively. Both cases pass the ANOVA (P < α ). Nevertheless, the Tukey test results in Table 7 and Table 8 reveal that the detection is achieved at higher damage levels than Taylor–Fourier filters. In the case of 75% of load, the start of the detection goes from 1 mm to 4 mm; meanwhile, for the 50% load condition, the detection decreases from 1 mm to 8 mm.
These results demonstrate that Taylor–Fourier filters effectively detect damages in the bars of the rotor at earlier stages.

5. Conclusions

A methodology for early IMs broken bar detection based on the DTFT was presented and demonstrated. A total of 1100 current signals were analyzed, considering scenarios going from a healthy bar to a broken bar. For both load conditions, the Tukey test allowed an incipient detection of 1 mm damage in the bar with a confidence level of 99%; for the 75% loaded motor, only the 3 mm case could not be detected. Although the median value for the 3 mm distribution is larger than that for the healthy bar, the data dispersion tends to have smaller amplitude values. This drawback can be attributed not to the 3 mm case but to the healthy bar distribution, whose current signals presented a slight damping that tends to increase its estimated amplitude. However, the signals were purposefully not preprocessed to assess the reach of the proposed method for early fault detection. Overall, it is considered that the DTFT implementation for current signal analysis was successfully carried out to perform early fault detection.
The use of the O-splines of the DTFT enabled smooth filtering without an increase in the processing time. This is an important advantage regarding future implementation in an FPGA. In addition, the estimated fault amplitude was the only feature used for detection, which is a straightforward calculation due to the use of the O-splines.
In future works, we are interested in extracting more features with the TFT, such as the phase of the signals, in performing a classification of the damage levels in the rotor bar. Besides, this system coded in Matlab will be traduced to its equivalent in a hardware description language (VHDL) to implement it in an FPGA, aiming to assess the system’s performance for an IM operating under real circumstances.

Author Contributions

Conceptualization, S.A.-T., G.A.-A., J.d.J.R.-M. and M.R.A.P.; methodology, S.A.-T. and G.A.-A.; software, S.A.-T.; validation, S.A.-T.; formal analysis, S.A.-T. and J.d.J.R.-M.; investigation, S.A.-T.; resources, J.d.J.R.-M. and M.R.A.P.; writing—original draft preparation, S.A.-T. and G.A.-A.; writing—review and editing, S.A.-T., G.A.-A., J.d.J.R.-M. and M.R.A.P.; visualization, S.A.-T.; supervision, J.d.J.R.-M. and M.R.A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

Authors acknowledge financial support from CONACYT in the program of Master of Science with grant numbers 776219 and 776238.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
IMInduction motors
MCSAMotor current signature analysis
FFTFast Fourier transform
AIArtificial intelligence
SVMSupport vector machines
ANNArtificial neural networks
TFTTaylor–Fourier transform
DTFTDigital Taylor–Fourier transform
FIRFinite impulse response
MCSAMotor current signature analysis
ANOVAAnalysis of variance
PSDPower spectral density

References

  1. de Jesús Rangel-Magdaleno, J. Induction Machines Fault Detection: An Overview. IEEE Instrum. Meas. Mag. 2021, 24, 63–71. [Google Scholar] [CrossRef]
  2. Kumar, P.; Hati, A. Review on Machine Learning Algorithm Based Fault Detection in Induction Motors. Arch. Comput. Methods Eng. 2020, 28, 1929–1940. [Google Scholar] [CrossRef]
  3. Liang, X.; Ali, M.Z.; Zhang, H. Induction Motors Fault Diagnosis Using Finite Element Method: A Review. IEEE Trans. Ind. Appl. 2020, 56, 1205–1217. [Google Scholar] [CrossRef]
  4. Li, Z.; Li, Y.; Sun, Q.; Qi, B. Bearing Fault Diagnosis Method Based on Convolutional Neural Network and Knowledge Graph. Entropy 2022, 24, 1589. [Google Scholar] [CrossRef] [PubMed]
  5. Morales-Perez, C.; Rangel-Magdaleno, J.; Peregrina-Barreto, H.; Amezquita-Sanchez, J.P.; Valtierra-Rodriguez, M. Incipient Broken Rotor Bar Detection in Induction Motors Using Vibration Signals and the Orthogonal Matching Pursuit Algorithm. IEEE Trans. Instrum. Meas. 2018, 67, 2058–2068. [Google Scholar] [CrossRef]
  6. Panagiotou, P.A.; Arvanitakis, I.; Lophitis, N.; Antonino-Daviu, J.A.; Gyftakis, K.N. A New Approach for Broken Rotor Bar Detection in Induction Motors Using Frequency Extraction in Stray Flux Signals. IEEE Trans. Ind. Appl. 2019, 55, 3501–3511. [Google Scholar] [CrossRef] [Green Version]
  7. Moríñigo Sotelo, D.; Garcia-Calva, T.; Romero-Troncoso, R.; Fernandez-Cavero, V.; Garcia-Perez, A. Early Detection of Broken Rotor Bars in Inverter-Fed Induction Motors Using Speed Analysis of Startup Transients. Energies 2021, 14, 1469. [Google Scholar] [CrossRef]
  8. Ayala, J.; Meira, M.; Verucchi, C.; Bossio, G.; Acosta, G. Combined Vibration and Stator Current Techniques for Induction Motors Fault Detection—A Review. In Proceedings of the 2021 XIX Workshop on Information Processing and Control (RPIC), San Juan, Argentina, 3–5 November 2021; pp. 1–6. [Google Scholar] [CrossRef]
  9. Kabul, A.; Ünsal, A. An Alternative Approach for the Detection of Broken Rotor Bars and Bearing Faults of Induction Motor Based on Vibration Signals. In Proceedings of the 2021 8th International Conference on Electrical and Electronics Engineering (ICEEE), Antalya, Turkey, 9–11 April 2021; pp. 126–131. [Google Scholar] [CrossRef]
  10. Rangel-Magdaleno, J.d.J.; Romero-Troncoso, R.d.J.; Osornio-Rios, R.A.; Cabal-Yepez, E.; Contreras-Medina, L.M. Novel Methodology for Online Half-Broken-Bar Detection on Induction Motors. IEEE Trans. Instrum. Meas. 2009, 58, 1690–1698. [Google Scholar] [CrossRef]
  11. Valles-Novo, R.; de Jesus Rangel-Magdaleno, J.; Ramirez-Cortes, J.M.; Peregrina-Barreto, H.; Morales-Caporal, R. Empirical Mode Decomposition Analysis for Broken-Bar Detection on Squirrel Cage Induction Motors. IEEE Trans. Instrum. Meas. 2015, 64, 1118–1128. [Google Scholar] [CrossRef]
  12. da Costa, C.; Kashiwagi, M.; Mathias, M.H. Rotor failure detection of induction motors by wavelet transform and Fourier transform in non-stationary condition. Case Stud. Mech. Syst. Signal Process. 2015, 1, 15–26. [Google Scholar] [CrossRef] [Green Version]
  13. Garcia-Perez, A.; Romero-Troncoso, R.d.J.; Cabal-Yepez, E.; Osornio-Rios, R.A. The Application of High-Resolution Spectral Analysis for Identifying Multiple Combined Faults in Induction Motors. IEEE Trans. Ind. Electron. 2011, 58, 2002–2010. [Google Scholar] [CrossRef]
  14. Trujillo-Guajardo, L.A.; Rodriguez-Maldonado, J.; Moonem, M.A.; Platas-Garza, M.A. A Multiresolution Taylor–Kalman Approach for Broken Rotor Bar Detection in Cage Induction Motors. IEEE Trans. Instrum. Meas. 2018, 67, 1317–1328. [Google Scholar] [CrossRef]
  15. Reljič, D.; Jerkan, D.; Kanovič, Ž. Broken Rotor Bar Fault Detection Using Advanced IM Model and Artificial Intelligence Approach. In Proceedings of the IEEE EUROCON 2019—18th International Conference on Smart Technologies, Novi Sad, Serbia, 1–4 July 2019; pp. 1–6. [Google Scholar] [CrossRef]
  16. Xiao, D.; Huang, Y.; Qin, C.; Shi, H.; Li, Y. Fault Diagnosis of Induction Motors Using Recurrence Quantification Analysis and LSTM with Weighted BN. Shock Vib. 2019, 2019, 8325218. [Google Scholar] [CrossRef] [Green Version]
  17. Valtierra-Rodriguez, M.; Rivera-Guillen, J.R.; Basurto-Hurtado, J.A.; De-Santiago-Perez, J.J.; Granados-Lieberman, D.; Amezquita-Sanchez, J.P. Convolutional Neural Network and Motor Current Signature Analysis during the Transient State for Detection of Broken Rotor Bars in Induction Motors. Sensors 2020, 20, 3721. [Google Scholar] [CrossRef] [PubMed]
  18. Kumar, P.; Hati, A.S. Dilated convolutional neural network based model for bearing faults and broken rotor bar detection in squirrel cage induction motors. Expert Syst. Appl. 2022, 191, 116290. [Google Scholar] [CrossRef]
  19. de la O Serna, J.A. Analyzing Power Oscillating Signals With the O-Splines of the Discrete Taylor–Fourier Transform. IEEE Trans. Power Syst. 2018, 33, 7087–7095. [Google Scholar] [CrossRef]
  20. Avalos-Almazan, G.; Castillo-Garcia, G.; Paternina, M.; Zamora, A.; Fuentes, J.; Rodriguez-Rodriguez, J.; Guillen, D. Real-Time Phasor Estimation via the Taylor-Fourier’s Subspace. In Proceedings of the 2018 North American Power Symposium (NAPS), Fargo, ND, USA, 9–11 September 2018; pp. 1–6. [Google Scholar] [CrossRef]
  21. Platas-Garza, M.A.; de la O Serna, J.A. Dynamic Harmonic Analysis Through Taylor–Fourier Transform. IEEE Trans. Instrum. Meas. 2011, 60, 804–813. [Google Scholar] [CrossRef]
  22. Serna, J.A.d.l.O. Dynamic Harmonic Analysis With FIR Filters Designed With O-Splines. IEEE Trans. Circuits Syst. I Regul. Pap. 2020, 67, 5092–5100. [Google Scholar] [CrossRef]
  23. de la O Serna, J.A.; Paternina, M.R.A.; Zamora-Méndez, A.; Tripathy, R.K.; Pachori, R.B. EEG-Rhythm Specific Taylor–Fourier Filter Bank Implemented With O-Splines for the Detection of Epilepsy Using EEG Signals. IEEE Sens. J. 2020, 20, 6542–6551. [Google Scholar] [CrossRef]
  24. Toma, R.N.; Prosvirin, A.E.; Kim, J.M. Bearing Fault Diagnosis of Induction Motors Using a Genetic Algorithm and Machine Learning Classifiers. Sensors 2020, 20, 1884. [Google Scholar] [CrossRef] [Green Version]
  25. Kyun, K.T. Understanding one-way ANOVA using conceptual figures. kja 2017, 70, 22–26. [Google Scholar] [CrossRef] [Green Version]
  26. Górecki, T.; Smaga. A comparison of tests for the one-way ANOVA problem for functional data. Comput. Stat. 2015, 30, 987–1010. [Google Scholar] [CrossRef] [Green Version]
  27. Larson, M.G. Analysis of Variance. Circulation 2008, 117, 115–121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Abdi, H.; Williams, L.J. Newman-Keuls test and Tukey test. Encycl. Res. Des. 2010, 2, 897–902. [Google Scholar]
  29. de la O Serna, J.A.; Paternina, M.A.; Zamora-Mendez, A. Assessing Synchrophasor Estimates of an Event Captured by a Phasor Measurement Unit. IEEE Trans. Power Deliv. 2021, 36, 3109–3117. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed methodology for broken bar detection in IM.
Figure 1. Flowchart of the proposed methodology for broken bar detection in IM.
Entropy 25 00044 g001
Figure 2. Filter response according to the value of K with a central frequency of 60 Hz.
Figure 2. Filter response according to the value of K with a central frequency of 60 Hz.
Entropy 25 00044 g002
Figure 3. Example of classification using statistical distribution.
Figure 3. Example of classification using statistical distribution.
Entropy 25 00044 g003
Figure 4. Rotor bar drilling to create the damage levels.
Figure 4. Rotor bar drilling to create the damage levels.
Entropy 25 00044 g004
Figure 5. Bar damage levels under analysis.
Figure 5. Bar damage levels under analysis.
Entropy 25 00044 g005
Figure 6. Current signals of healthy and broken bar IM.
Figure 6. Current signals of healthy and broken bar IM.
Entropy 25 00044 g006
Figure 7. Frequency spectra of the signals for a half-broken bar IM.
Figure 7. Frequency spectra of the signals for a half-broken bar IM.
Entropy 25 00044 g007
Figure 8. O-spline for K = 5 .
Figure 8. O-spline for K = 5 .
Entropy 25 00044 g008
Figure 9. Filter with f c = 54 Hz and a 60 Hz synthetic signal.
Figure 9. Filter with f c = 54 Hz and a 60 Hz synthetic signal.
Entropy 25 00044 g009
Figure 10. Boxplot of the fault amplitudes’ distribution data in a 75% loaded IM.
Figure 10. Boxplot of the fault amplitudes’ distribution data in a 75% loaded IM.
Entropy 25 00044 g010
Figure 11. Statistical distribution for a 50% loaded motor centered at 56 Hz.
Figure 11. Statistical distribution for a 50% loaded motor centered at 56 Hz.
Entropy 25 00044 g011
Figure 12. Fourier filter response.
Figure 12. Fourier filter response.
Entropy 25 00044 g012
Figure 13. Statistical distribution for a 75% loaded motor using a Fourier filter centered at 54 Hz.
Figure 13. Statistical distribution for a 75% loaded motor using a Fourier filter centered at 54 Hz.
Entropy 25 00044 g013
Figure 14. Statistical distribution for a 50% loaded motor using a Fourier filter centered at 56 Hz.
Figure 14. Statistical distribution for a 50% loaded motor using a Fourier filter centered at 56 Hz.
Entropy 25 00044 g014
Table 1. ANOVA results for 75% loaded motor.
Table 1. ANOVA results for 75% loaded motor.
Between Groups Within GroupsOverall
Sum of squares0.28280.00550.2883
Degrees of freedom10539549
Mean sum of squares0.02831.03 × 10 5
F2758.8433
P0
Table 2. Tukey test results for 75% loaded motor, Q c r i t i c a l = 5.227 .
Table 2. Tukey test results for 75% loaded motor, Q c r i t i c a l = 5.227 .
ComparisonQResults
Healthy vs. 1 mm 17.2663Significantly different
Healthy vs. 2 mm42.4258Significantly different
Healthy vs. 3 mm22.8229Significantly different
Healthy vs. 4 mm12.6416Significantly different
Healthy vs. 5 mm17.7610Significantly different
Healthy vs. 6 mm29.2188Significantly different
Healthy vs. 7 mm52.2714Significantly different
Healthy vs. 8 mm113.1296Significantly different
Healthy vs. 9 mm113.2316Significantly different
Healthy vs. broken147.8344Significantly different
Table 3. ANOVA results for 50% loaded motor.
Table 3. ANOVA results for 50% loaded motor.
Between GroupsWithin GroupsOverall
Sum of squares0.10320.00140.1046
Degrees of freedom10539549
Mean sum of squares0.01032.60 × 10 6
F3968.8528
P0
Table 4. Tukey test results for 50% loaded motor, Q c r i t i c a l = 5.227 .
Table 4. Tukey test results for 50% loaded motor, Q c r i t i c a l = 5.227 .
ComparisonQResults
Healthy vs. 1 mm5.9463Significantly different
Healthy vs. 2 mm14.95354Significantly different
Healthy vs. 3 mm1.3155Not significantly different
Healthy vs. 4 mm11.6383Significantly different
Healthy vs. 5 mm10.7525Significantly different
Healthy vs. 6 mm41.5717Significantly different
Healthy vs. 7 mm44.2730Significantly different
Healthy vs. 8 mm101.2017Significantly different
Healthy vs. 9 mm166.3308Significantly different
Healthy vs. broken161.2966Significantly different
Table 5. ANOVA results for 75% loaded motor using a Fourier filter.
Table 5. ANOVA results for 75% loaded motor using a Fourier filter.
Between GroupsWithin GroupsOverall
Sum of squares 0.25150.10380.3553
Degrees of freedom10539549
Mean sum of squares0.02510.0002
F130.5811
P5.6 × 10 137
Table 6. ANOVA results for 50% loaded motor using a Fourier filter.
Table 6. ANOVA results for 50% loaded motor using a Fourier filter.
Between GroupsWithin GroupsOverall
Sum of squares 0.14820.18100.3292
Degrees of freedom10539549
Mean sum of squares0.01480.0003
F44.1252
P 9.48 × 10 64
Table 7. Tukey test results for 75% loaded motor using a Fourier filter, Q c r i t i c a l = 5.227 .
Table 7. Tukey test results for 75% loaded motor using a Fourier filter, Q c r i t i c a l = 5.227 .
ComparisonQResults
Healthy vs. 1 mm 1.3594Not significantly different
Healthy vs. 2 mm1.3737Not significantly different
Healthy vs. 3 mm1.5214Not significantly different
Healthy vs. 4 mm5.6934Significantly different
Healthy vs. 5 mm4.7304Not significantly different
Healthy vs. 6 mm6.2172Significantly different
Healthy vs. 7 mm11.3767Significantly different
Healthy vs. 8 mm10.4901Significantly different
Healthy vs. 9 mm29.3997Significantly different
Healthy vs. broken32.3060Significantly different
Table 8. Tukey test results for 50% loaded motor using a Fourier filter, Q c r i t i c a l = 5.227 .
Table 8. Tukey test results for 50% loaded motor using a Fourier filter, Q c r i t i c a l = 5.227 .
ComparisonQResults
Healthy vs. 1 mm 0.9307Not significantly different
Healthy vs. 2 mm2.2843Not significantly different
Healthy vs. 3 mm0.1505Not significantly different
Healthy vs. 4 mm0.0857Not significantly different
Healthy vs. 5 mm2.6871Not significantly different
Healthy vs. 6 mm2.0041Not significantly different
Healthy vs. 7 mm4.8410Not significantly different
Healthy vs. 8 mm10.3565Significantly different
Healthy vs. 9 mm18.6478Significantly different
Healthy vs. broken14.5947Significantly different
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aguayo-Tapia, S.; Avalos-Almazan, G.; Rangel-Magdaleno, J.d.J.; Paternina, M.R.A. Broken Bar Fault Detection Using Taylor–Fourier Filters and Statistical Analysis. Entropy 2023, 25, 44. https://doi.org/10.3390/e25010044

AMA Style

Aguayo-Tapia S, Avalos-Almazan G, Rangel-Magdaleno JdJ, Paternina MRA. Broken Bar Fault Detection Using Taylor–Fourier Filters and Statistical Analysis. Entropy. 2023; 25(1):44. https://doi.org/10.3390/e25010044

Chicago/Turabian Style

Aguayo-Tapia, Sarahi, Gerardo Avalos-Almazan, Jose de Jesus Rangel-Magdaleno, and Mario R. A. Paternina. 2023. "Broken Bar Fault Detection Using Taylor–Fourier Filters and Statistical Analysis" Entropy 25, no. 1: 44. https://doi.org/10.3390/e25010044

APA Style

Aguayo-Tapia, S., Avalos-Almazan, G., Rangel-Magdaleno, J. d. J., & Paternina, M. R. A. (2023). Broken Bar Fault Detection Using Taylor–Fourier Filters and Statistical Analysis. Entropy, 25(1), 44. https://doi.org/10.3390/e25010044

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop