Next Article in Journal
Ageing, Shocks and Wear Mechanisms in ZTA and the Long-Term Performance of Hip Joint Materials
Next Article in Special Issue
Fault Diagnosis for Rolling Bearings under Variable Conditions Based on Visual Cognition
Previous Article in Journal
Improvement in Char Strength with an Open Cage Silsesquioxane Flame Retardant
Previous Article in Special Issue
The Feasibility of Structural Health Monitoring Using the Fundamental Shear Horizontal Guided Wave in a Thin Aluminum Plate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bearing Fault Detection Based on Empirical Wavelet Transform and Correlated Kurtosis by Acoustic Emission

1
Shanxi Key Laboratory of Mechanical Product Quality Assurance and Diagnostics, Xi’an Jiaotong University, Xi’an 710049, China
2
State Key Laboratory of Manufacturing System Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Materials 2017, 10(6), 571; https://doi.org/10.3390/ma10060571
Submission received: 15 April 2017 / Revised: 12 May 2017 / Accepted: 17 May 2017 / Published: 24 May 2017
(This article belongs to the Special Issue Structural Health Monitoring for Aerospace Applications 2017)

Abstract

:
Rolling bearings are widely used in rotating equipment. Detection of bearing faults is of great importance to guarantee safe operation of mechanical systems. Acoustic emission (AE), as one of the bearing monitoring technologies, is sensitive to weak signals and performs well in detecting incipient faults. Therefore, AE is widely used in monitoring the operating status of rolling bearing. This paper utilizes Empirical Wavelet Transform (EWT) to decompose AE signals into mono-components adaptively followed by calculation of the correlated kurtosis (CK) at certain time intervals of these components. By comparing these CK values, the resonant frequency of the rolling bearing can be determined. Then the fault characteristic frequencies are found by spectrum envelope. Both simulation signal and rolling bearing AE signals are used to verify the effectiveness of the proposed method. The results show that the new method performs well in identifying bearing fault frequency under strong background noise.

1. Introduction

Rolling bearings are widely used in industry equipment and are key components of rotating machinery. Bearings are inevitably damaged because of overload and long-time operation. Fault bearings make vibration of equipment exceeding operation standard. Badly damaged bearings even lock the whole system and destroy the rotating axle, which initiates downtime and financial loss [1]. Therefore monitoring bearing fault is of great importance to ensure a safe operation of a machine.
Various kinds of monitoring technologies are employed to detect bearing faults, such as encoder, oil, ferro-graphic, vibration and acoustic emission (AE) [2]. Among them, the vibration signal is the most widely used. However, since AE signal can detect early faults of machinery components, the bearing fault monitoring based on AE has been of interest to researchers in recent years [3]. AE is transient stress wave rapidly released by localized energy [4]. During the operation of bearing, main sources of AE include impacting and friction of bearing rollers passing over outer and inner race and fatigue damage of bearing components [5]. Both traditional AE characteristic parameters and newer signal processing methods are employed to analyze AE signals.
Traditional AE characteristic parameters are studied thoroughly to detect rolling bearing faults [6,7,8]. The research results indicate that AE technology has better performance over vibration analysis for monitoring early deterioration of rolling bearing [3]. Morhain and Mba studied the relationship of RMS, energy and count values with rotational speed, load and the size of the faults [8]. Abdullah first used duration time of AE pulses to indicate the defect size of the outer race [9]. Since AE signal is always disturbed by noise, Khamis first proposed the Energy Index (EI) technique to enhance the AE burst buried in random noise. The EI index is utilized to detect incipient faults of bearings [5].
A variety of signal processing methods based on AE technology are used to monitoring bearing condition. Kurtogram and spectral kurtosis are applied to enhance AE events [10]. Kilundu first investigated cyclo-stationarity characteristic of AE bearing signals and used it to find outer race defects [11]. Pandya utilized Hilbert-Huang Transform (HHT) to obtain the time-frequency distribution of AE bearing signals. Then the K-nearest neighbor algorithm based on asymmetric proximity function was employed to classify the signal features, by which the reliability of bearing fault diagnosis increased [12]. Zvokelj proposed a new method, Ensemble Empirical Mode Decomposition-based multiscale Independent Component Analysis, to monitoring bearing fault. The algorithm can not only detect bearing faults but also de-noise AE signals [13]. Farzad Hemmati further developed the bearing fault detection technology by optimizing the ratio of Kurtosis and calculating Shannon entropy [14]. Mourad Kedadouche first used EWT to AE bearing signals to find frequency band which is excited more than other bands. The proposed method successfully detected defects as low as 40 microns on bearing outer race [15].
Although a lot of investigations are devoted to detect bearing faults, there are still two problems to be solved. Firstly, most of the methods are only effective in outer race defect, but ineffective in inner race fault due to signal attenuation. Secondly, AE signals are easily disturbed by background noise, so method with good anti-noise performance is needed. Therefore, in this paper a new method based on Empirical Wavelet Transform (EWT) and correlated kurtosis (CK) is proposed. Analysis on simulated signal proves that the resonant frequency band and fault characteristic frequency can be effectively obtained with this method. The rolling bearing experiment is carried out and AE signals of defect bearing are used to prove the efficiency of the proposed method.
The rest of the paper is organized as follows. The EWT algorithm and CK are introduced in Section 2. In Section 3, the proposed method is applied to simulation signal to validate its effectiveness. Then, in Section 4, rolling bearing experiment is carried out and the method is employed to detect the inner race fault of bearing. Section 5 summarizes the conclusion.

2. Introduction to the Theory

2.1. Empirical Wavelet Transform

EWT divides signals into frequency sub-bands adaptively according to the information located in the spectrum of the original signal [16,17]. The signal can be decomposed into sub-bands by a set of wavelet construction. The detail of the EWT algorithm is as follows.
Suppose that the signal has a Fourier support within [ 0 , π ] . The spectrum can be divided into N th contiguous partitions. Each component has the frequency bandwidth as: Λ n = [ ω n 1 , ω n ] , where ω 0 = 0 , ω N = π , n = 1 N Λ n = [ 0 , π ] . The empirical scaling function and the empirical wavelets are shown in Equations (1) and (2) respectively.
Φ n ( ω ) = { 1 if   | ω | ω n τ n cos [ π 2 β ( 1 2 τ n ( | ω | ω n + τ n ) ) ] if   ω n τ n | ω | ω n + τ n 0 otherwise } ,
Ψ n ( ω ) = { 1 if   ω n + τ n | ω | ω n 1 τ n + 1 cos [ π 2 β ( 1 2 τ n 1 ( | ω | ω n + 1 + τ n + 1 ) ) ] if   ω n + 1 τ n + 1 | ω | ω n + 1 + τ n + 1 sin [ π 2 β ( 1 2 τ n ( | ω | ω n + τ n ) ) ] if   ω n τ n | ω | ω n + τ n 0 otherwise } ,
Here, T n is the transition phase corresponding to the ω n , which has the width of 2 τ n . The β ( χ ) is an arbitrary function. In this paper, it is defined as follows:
β ( χ ) = χ 4 ( 35 84 χ + 70 χ 2 20 χ 3 ) ,
Similar to the continuous wavelet transformation (CWT), the detail coefficients and approximation coefficients can be obtained by the empirical wavelets and scaling function, the functions are shown as follows:
W χ ( n , t ) = χ ( τ ) ψ n ( τ t ) d τ = F 1 ( X ( ω ) ψ n ( ω ) ¯ ) ,
W χ ( 0 , t ) = χ ( τ ) ϕ 1 ( τ t ) d τ = F 1 ( X ( ω ) ϕ 1 ( ω ) ¯ ) ,
The extracted modes can be obtained as below:
f 0 ( t ) = W χ ( 0 , t ) × ϕ 1 ( t ) ,
f k ( t ) = W χ ( k , t ) × ψ k ( t ) ,
In order to ensure the decomposition has a tight frame, Equation 8 must be satisfied.
γ < min ( ω n + 1 ω n ω n + 1 + ω n ) ,
where γ = τ n ω n . More detailed information can be found in reference [16].

2.2. Correlated Kurtosis

According to reference [18], the definition of CK is shown as follows:
C K M ( T ) = n = 1 N ( m = 0 M y n m T ) 2 ( n = 1 N y n 2 ) M + 1 ,
Here, M behalves the order of shift. N is the length of the input signal. T is the sampling point which can be calculated as following:
T = f s f c ,
where f s is the sampling frequency; f c is the fault character frequency. In practice, the T varies in a certain range due to the fluctuation of rotational speed.
From the definition, it is known that CK reveals the periodicity of signal. If T is the period of fault, the CK can indicate the occurrence of fault by its magnitude. The frequency band with maximum value of CK is the frequency band which is excited most by the fault characteristic frequency corresponding to T.
For rolling bearing, the fault frequencies are always modulated to the resonant frequency and the higher order resonant frequencies. Since the fault frequencies of rolling bearing have obvious periodicity, the CK at time interval T that equals to the fault period will achieve higher value at resonant frequency bands. Therefore, the frequency bands with the highest CK value can be considered as the resonant frequency or the higher order resonant frequencies.
According to the analysis above, a new method is proposed to find resonant frequency and fault frequency by EWT and CK. The whole procedure of the proposed algorithm is shown as a flowchart (Figure 1).

3. Simulation Verification

To verify the proposed method, a simulation signal is employed. When a local fault in a rolling bearing emerges, harmonic components and periodical impulses will be generated. These frequency components are modulated by resonance frequency of bearing and decay exponentially [19].
The simulation signal of defect bearing is expressed as follows:
x ( t ) = k = 1 M A 0 s ( t k T ) + n ( t ) ,
s ( t ) = e ζ t cos ( 2 π f n t + φ 0 ) ,
where A 0 is the impulse response amplitude; ζ is the damping ratio of the impulse; M is the number of impulses; T is the period of the impulses which is related to the fault characteristic frequency f c ; f n is the resonance frequency; φ 0 is the primary phase; n ( t ) is white noise. The value of these parameters are displayed as following:
A 0 = 1 ;   ζ = 160 ;   M = 3 ;   T = 0.05 s ;   f c = 20 Hz ;   f n = 1300 Hz ;   φ 0 = 0 ;   S N R = 20 dB .
First, the EWT is utilized to decompose the simulated signal adaptively. Then the CK of these modes are calculated at the given value of T. The waveforms of the simulated signal with and without noise are shown in Figure 2. Figure 3 displays the spectrum and envelope spectrum of the simulation signal with noise. It shows that the frequency 20 Hz is buried in strong background noise. The spectrum segmentation by EWT is shown in Figure 4. It is hard to find resonant frequency only by the amplitude of the spectrum. Figure 5 is the CK value of each component at time interval T. The peak frequency appears at 1275 Hz. It approximates to the true value of resonant frequency, which is 1300 Hz. The error between them can be ascribed to the strong noise. Figure 6 is the envelope spectrum of frequency band centered around 1275 Hz. In Figure 6, the frequency 20 Hz and its harmonics are obvious. It is clearly that the CK performs well in identifying resonant frequency and periodic frequency under low SNR.

4. Experimental Verification

Our rolling bearing experiment is explained in this section. The experimental setup is displayed in Figure 7. It consisted of motor diver, shaft, bearings, bearing housing, inertia disc and control system. The shaft was supported by two bearings. The left one was normal and the right one defective. An inner race fault rolling bearings was used in this experiment. The defect was caused artificially. The rotation speed of the experiment was 3000 rpm.
The AE signals were picked up in the experiment. By comparing performance of three kinds of AE sensors, Nano 30 was selected. The AE sensor was mounted on the top surface of bearing housing (Figure 7). The parameters of Nano 30 is shown in Table 1. The original AE signals were amplified by 40 dB. The AE signal collection system was PCI-II (Physical Acoustics Corporation, Boston, MA, USA). The sample rate was 500 kHz and the AE threshold level was set to 40 dB.
When defects occur in the rolling bearings, some certain frequency components will be generated according to the rotational speed, geometric parameters of bearings and the type of faults [20]. The formulae of the defective frequencies are as follows:
Ball-pass frequency of outer race:
BPFO   =   n f r 2 { 1 d D cos ϕ } ,
Ball-pass frequency of inner race:
BPFI   =   n f r 2 { 1 + d D cos ϕ } ,
Fundamental train frequency (cage speed):
FTF   =   f r 2 { 1 d D cos ϕ } ,
Ball (roller) spin frequency:
BSF = D f r 2 d { 1 ( d D cos ϕ ) 2 }
Here, fr is the shaft rotation frequency, n is the number of rolling elements, and a is the contact angle of the load. D and d represent the diameter and pitch diameter of rolling element respectively.
The fault feature frequencies of the rolling bearing are listed in Table 2.
Figure 8 shows waveform and spectrum of rolling bearing with inner race faults. Figure 9 is the envelope spectrum of AE bearing signal. It can be seen that the rotating frequency and its harmonics are obvious, but the BPFI which is defined in Equation (14) and provided in Table 2 cannot be found in the figure. Then the EWT is carried out to decompose the bearing AE signal. The segmentation of the spectrum is shown in Figure 10. Figure 11 displays CK value of components at time interval T which equals to inner fault period. The maximum value appears at 66k Hz. It means that this frequency band is excited most by the inner race fault frequency. Figure 12 is the spectrum envelope of this frequency band. It is clear that the rotating frequency is also modulated to this band. However, in Figure 12, the amplitude of the inner race fault frequency increases and the higher order of inner race fault frequencies are clear. The result shows that the proposed method can improve the inner race fault detection of rolling bearing.
It must be pointed out that the time interval T is the reciprocal of fault character frequency and the fault character frequency of rolling bearing is proportional to rotating speed, so the speed variation will affect the time interval T. If the speed variation is small, the accurate time interval T changes within certain range centered around the reciprocal of the fault character frequency. If the speed variation is large, the proposed method will not work.

5. Conclusions

In this paper, a new method is proposed to detect defects in rolling bearing. EWT is utilized to decompose AE signals into different components. CK of these components at a certain time interval T are calculated. Components with the largest CK value is considered to be the resonant frequency. Then envelope demodulation is carried out to obtain the fault feature characteristics of corresponding frequency bands. Simulated and experimental signals were employed to verify the effectiveness of the method. Two conclusions can be drawn as follows:
  • With the right time interval T, the resonant frequency of the rolling bearing system can be obtained by calculating the CK value.
  • Weak fault features can be identified by the proposed method under low SNR conditions.
Future research work will be devoted to enhance the fault frequencies modulated by different resonant frequencies and estimation of the accurate time interval T under slight speed fluctuation will also be carried out.

Acknowledgments

This work was funded by the National Natural Science Foundation of China [Grant No. 51421004; 51505365]; and the Fundamental Research Funds of the Central Universities of China [Grant No. CXTD2014001].

Author Contributions

Zheyu Gao developed the methodology and wrote the manuscript. Jing Lin provided the professional guidance to complete the research work and manuscript revision. Xiufeng Wang provided great assistance in doing experiments and provided professional guidance. Xiaoqiang Xu provided great advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hecke, B.E.V. Development of Novel Acoustic Emission Based Methodology and Tools for Bearing Fault Diagnostics. Ph.D. Thesis, University of Illinois, Chicago, IL, USA, 2015. [Google Scholar]
  2. Cao, H.R.; Fan, F.; Zhou, K.; He, Z.J. Wheel-bearing fault diagnosis of trains using empirical wavelet transform. Measurement 2016, 82, 439–449. [Google Scholar] [CrossRef]
  3. Mba, D.; Rao, R. Development of Acoustic Emission technology for the condition monitoring and diagnosis of rotating machines: bearings, pumps, gearboxes, engines and rotating structures. Shock Vib. Dig. 2006, 38, 3–16. [Google Scholar] [CrossRef]
  4. Tandon, N.; Choudhury, A. A review of vibration and acoustic measurement methods for the detection of defects in rolling element bearings. Tribol. Int. 1999, 32, 469–480. [Google Scholar] [CrossRef]
  5. Al-Balushi, K.R.; Addali, A.; Charnley, B.; Mba, D. Energy Index technique for detection of Acoustic Emissions associated with incipient bearing failures. Appl. Acoust. 2010, 71, 812–821. [Google Scholar] [CrossRef]
  6. Balerston, H.L. The detection of incipient failure in bearings. Mater. Eval. 1969, 27, 121–128. [Google Scholar]
  7. Yoshioka, T.; Korenaga, A.; Mano, H.; Yamamoto, T. Diagnosis of rolling bearing by measuring time interval of Acoustic Emission generation. J. Tribol. 1999, 121, 468–472. [Google Scholar] [CrossRef]
  8. Morhain, A.; Mba, D. Bearing Defect Diagnosis and Acoustic Emission. Proc. Inst. Mech. Eng. Part J: J. Eng. Tribol. 2003, 217, 272–275. [Google Scholar] [CrossRef]
  9. Al-Ghamd, A.M.; Mba, D. A comparative experimental study on the use of acoustic emission and vibration analysis for bearing defect identification and estimation of defect size. Mech. Syst. Signal Process. 2006, 20, 1537–1571. [Google Scholar] [CrossRef]
  10. Eftekharnejad, B.; Carrasco, M.R.; Charnley, B.; Mba, D. The application of spectral kurtosis on Acoustic Emission and vibrations from a defective bearing. Mech. Syst. Signal Process. 2011, 25, 266–284. [Google Scholar] [CrossRef]
  11. Kilundu, B.; Chiementin, X.; Duez, J.; Mba, D. Cyclostationarity of Acoustic Emissions (AE) for monitoring bearing defects. Mech. Syst. Signal Process. 2011, 25, 2061–2072. [Google Scholar] [CrossRef]
  12. Pandya, D.H.; Upadhyay, S.H.; Harsha, S.P. Fault diagnosis of rolling element bearing with intrinsic mode function of acoustic emission data using APF-KNN. Expert Syst. Appl. 2013, 40, 4137–4145. [Google Scholar] [CrossRef]
  13. Žvokelj, M.; Zupan, S.; Prebil, I. EEMD-based multiscale ICA method for slewing bearing fault detection and diagnosis. J. Sound Vib. 2016, 370, 394–423. [Google Scholar] [CrossRef]
  14. Hemmati, F.; Orfali, W.; Gadala, M.S. Roller bearing acoustic signature extraction by wavelet packet transform, applications in fault detection and size estimation. Appl. Acoust. 2016, 104, 101–118. [Google Scholar] [CrossRef]
  15. Kedadouche, M.; Thomas, M.; Tahan, A. Empirical Mode Decomposition combined with Empirical Wavelets for extracting bearing frequencies in a noisy environment and early detection of defects. In Advances in Condition Monitoring of Machinery in Non-Stationary Operations; Springer: Lyon, France, 2014; pp. 151–165. [Google Scholar]
  16. Gilles, J. Empirical Wavelet Transform. IEEE Trans. Signal Process. 2013, 61, 3999–4002. [Google Scholar]
  17. Pan, J.; Chen, J.L.; Zi, Y.Y.; Li, Y.M.; He, Z.J. Mono-component feature extraction for mechanical fault diagnosis using modified empirical wavelet transform via data-driven adaptive Fourier spectrum segment. Mech. Syst. Signal Process. 2016, 72–73, 160–183. [Google Scholar] [CrossRef]
  18. McDonald, G.L.; Zhao, Q.; Zuo, M.J. Maximum correlated Kurtosis deconvolution and application on gear tooth chip fault detection. Mech. Syst. Signal Process. 2012, 33, 237–255. [Google Scholar] [CrossRef]
  19. McFadden, P.D.; Smith, J.D. Model for the vibration produced by a single point defect in a rolling element bearing. J. Sound Vib. 1984, 96, 69–82. [Google Scholar] [CrossRef]
  20. Randall, R.B.; Antoni, J. Rolling element bearing diagnostics—A tutorial. Mech. Syst. Signal Process. 2011, 25, 485–520. [Google Scholar] [CrossRef]
Figure 1. Flowchart of algorithm.
Figure 1. Flowchart of algorithm.
Materials 10 00571 g001
Figure 2. The simulated impulse signal and simulated signal with noise.
Figure 2. The simulated impulse signal and simulated signal with noise.
Materials 10 00571 g002
Figure 3. The spectrum and envelope spectrum of simulation signal with noise.
Figure 3. The spectrum and envelope spectrum of simulation signal with noise.
Materials 10 00571 g003
Figure 4. The segmentation of spectrum of simulation signal.
Figure 4. The segmentation of spectrum of simulation signal.
Materials 10 00571 g004
Figure 5. Correlated kurtosis (CK) value of components with time interval equals to T.
Figure 5. Correlated kurtosis (CK) value of components with time interval equals to T.
Materials 10 00571 g005
Figure 6. Envelope spectrum of component with maximum CK value.
Figure 6. Envelope spectrum of component with maximum CK value.
Materials 10 00571 g006
Figure 7. Experimental set up.
Figure 7. Experimental set up.
Materials 10 00571 g007
Figure 8. The waveform and spectrum of AE bearing signal.
Figure 8. The waveform and spectrum of AE bearing signal.
Materials 10 00571 g008
Figure 9. The envelope spectrum of AE bearing signal.
Figure 9. The envelope spectrum of AE bearing signal.
Materials 10 00571 g009
Figure 10. The segmentation of spectrum of bearing signal.
Figure 10. The segmentation of spectrum of bearing signal.
Materials 10 00571 g010
Figure 11. CK value of components at time interval T.
Figure 11. CK value of components at time interval T.
Materials 10 00571 g011
Figure 12. Envelope spectrum of component with maximum CK value.
Figure 12. Envelope spectrum of component with maximum CK value.
Materials 10 00571 g012
Table 1. Parameters of Nano30.
Table 1. Parameters of Nano30.
ParametersRange
Operating Frequency Range125–750 KHz
Resonant Frequency300 KHz
Temperature Range−65–177 °C
Table 2. Fault character frequency of rolling bearing.
Table 2. Fault character frequency of rolling bearing.
Rotation Speed/rpmBPFO/HzBPFI/HzBSF/Hz
3000152.4247.599.6

Share and Cite

MDPI and ACS Style

Gao, Z.; Lin, J.; Wang, X.; Xu, X. Bearing Fault Detection Based on Empirical Wavelet Transform and Correlated Kurtosis by Acoustic Emission. Materials 2017, 10, 571. https://doi.org/10.3390/ma10060571

AMA Style

Gao Z, Lin J, Wang X, Xu X. Bearing Fault Detection Based on Empirical Wavelet Transform and Correlated Kurtosis by Acoustic Emission. Materials. 2017; 10(6):571. https://doi.org/10.3390/ma10060571

Chicago/Turabian Style

Gao, Zheyu, Jing Lin, Xiufeng Wang, and Xiaoqiang Xu. 2017. "Bearing Fault Detection Based on Empirical Wavelet Transform and Correlated Kurtosis by Acoustic Emission" Materials 10, no. 6: 571. https://doi.org/10.3390/ma10060571

APA Style

Gao, Z., Lin, J., Wang, X., & Xu, X. (2017). Bearing Fault Detection Based on Empirical Wavelet Transform and Correlated Kurtosis by Acoustic Emission. Materials, 10(6), 571. https://doi.org/10.3390/ma10060571

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