Bearing fault diagnosis based on TEO and SVM

Qingzhu Liu1 , Yujie Cheng2

1, 2School of Reliability and Systems Engineering, Beihang University, Beijing 100191, China

1, 2Science and Technology on Reliability and Environmental Engineering Laboratory, Beijing 100191, China

1Corresponding author

Vibroengineering PROCEDIA, Vol. 4, 2014, p. 206-210.
Accepted 13 September 2014; published 3 November 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

A fault method for bearing based on Teager energy operator (TEO) and support vector machine (SVM) is proposed in this paper. First, the total energy of the vibration signal of the bearing is estimated by the TEO technique, which has good time resolution for the instantaneous signal. Then, the Teager spectrums are obtained by applying fast Fourier transform (FFT) to the Teager energy signal. The feature frequencies of different fault modes, as well as the ratio of resonance frequency band energy to total energy in the Teager spectrum are extracted to form the feature vectors. Finally, these vectors are introduced into SVM to realize fault classification for the bearing. Experiments are conducted to verify the feasibility of the proposed method, the results show that the proposed method performs effectively to identify the failure mode of the bearing under variable conditions.

Keywords: fault diagnosis, bearing, TEO, SVM, FFT.

1. Introduction

Mechanics is becoming more and more important to the country with the development of the industry. There is no doubt that bearing, supporting shafts, is a critical component of rotating machinery and is easily damaged in the form of pitting corrosion, abrasion, erosion, rolling body injury, damage, and so on. If bearing is breakdown, it will cause enormous economic loss, influence industrial coherence and may even cause personal injury. So it is necessary to perform fault diagnosis for the bearing. Fortunately, many algorithms have been put forward to conduct analysis to the bearing fault.

Wavelet analysis and Empirical Mode Decomposition (EMD) have been used in signal processing and the feature extraction [1, 2]. The basic theory of wavelet analysis is to use a function to represent or approximate a signal. It has strong adaptive ability to the whole frequency domain. However, it is difficult to choose wavelet basis function. EMD can overcome this shortcoming and is adaptive to the change of signal, but it also has the disadvantage of end effect.

TEO can overcome all of the above shortcomings and be is available to the nonlinear signal. It has a high time resolution for the transient signal and output the total energy of the signal [3]. Kaiser improved the TEO to process voice signals and the result is very good [4]. Wang Tian-Jin used TEO to identify the characteristic frequency [5]. Hou Guo-Ping proposed a new bearing fault detection method based on SVM [6]. However, the training of neural networks needs large amount of data, which is hardly acquired in practice. Thus, SVM is proposed to realize fault classification for the bearing.

To verify the effectiveness of the proposed method, experiments have been conducted. The results show that the proposed method performs well in fault diagnosis for the bearing.

2. Methodology

2.1. Principle of TEO

2.1.1. TEO for continuous random signal

The TEO Ψ for continuous random signal is defined as [7]:

(1)
Ψ ( x ( t ) ) = [ x ˙ t ] 2 - x t x ¨ t ,

where x˙(t) and x¨(t) represent the first-order and the second-order derivative of x(t) respectively.

Actually, the total energy of the vibration signal of the bearing is estimated by the TEO technique. Now the formula in Eq. (1) is proved with a realistic model. This model is an undamped system composed of a mass block whose quality is m and a spring whose stiffness is k. According to Newton theorem, the motion equation of the mass is:

(2)
m x ¨ ( t ) + k x ( t ) = 0 ,

where x(t) is displacement of the mass relative to the equilibrium position. The solution of the equation is:

(3)
x t = A c o s ω t + φ ,
(4)
x ˙ t = - A ω s i n ω t + φ ,
(5)
x ¨ t = - A ω 2 c o s ω t + φ ,

where A is the vibration amplitude, ω is the inherent frequency, φ is the initial phase.

Then, mechanical energy is calculated, which is the sum of potential and kinetic:

(6)
E = 1 2 k [ x ( t ) ] 2 + 1 2 m [ x ˙ t ] 2 .

Drive Eqs. (3), (4) into Eq. (6), get:

(7)
E = 1 2 m A 2 ω 2 .

Now, use this model to solve the instantaneous energy output of the TEO. Substitute Eqs. (3), (4), (5) into Eq. (1), get:

(8)
Ψ x t = Ψ A c o s ω t + φ = A 2 ω 2 .

In contrast to Eq. (7), (8) missed m/2. So we can conclude that TEO track the total energy of simple harmonic motion signal (m/2 is a constant).

2.1.2. TEO for discrete signal

TEO Ψ for discrete signal is defined as:

(9)
Ψ ( x ( n ) ) = [ x n ] 2 - x n - 1 x n + 1 .

The energy of the signal only needs three data which can be calculated at any moment by TEO. So TEO have good time resolution for the instantaneous signal, and can detect transient component of signal.

2.2. SVM

The main idea of SVM is as shown in Fig. 1 [8]. The triangles represent one kind of samples, while circles represent another. H represents the classifying line. H1, H2 are lines parallel to H and they are through the samples which are the closest to H at the same time. The distance between H1 and H2 is defined as class interval. The role of H is to correctly separate two kinds of samples and reaches zero of training error rate. The optimal classifying line is reaching classification margin maximization in Fig. 1.

The equation for classification line is:

(10)
ω x + b = 0 .

Eq. (10) uses normalized treatment and gets linearly separable data setsxi,yi,(i=1, 2,..., n), xRd, y{1, –1} meet the following definitions

(11)
y i ω i x + b 1 ,       i = 1,   2,...,   n .

2/ω is class interval. The classification margin maximization is the minimum of ω2. The optimal classification line is to minimize the ω2/2 and meats the type Eq. (11). Support vector is the training samples on H1 and H2.

SVM is developed for a value of two proposed classification problems. However, the bearing fault classification is a multi-class classification problem. So the multi-class classification problem should be converted to two classification problems [9]. This experiment adopts a multi-classifier made of three SVMs.

Fig. 1. Optimal classifying line

Optimal classifying line

Fig. 2. The flow chart of the proposed method

 The flow chart of the proposed method

2.3. The procedure of the proposed method

This paper proposes a fault diagnosis method for bearing based on TEO and SVM. First, the total energy of the vibration signal of the bearing is estimated by the TEO technique. Second, the teager spectrums are obtained by applying fast Fourier transform (FFT) to the Teager energy signal [10, 11]. Third, feature vectors are extracted. Finally, these vectors are introduced into SVM to realize fault classification for the bearing.

3. Case study

3.1. Feature extraction

All the rolling bearing data used in this study are acquired by the Bearing Data Center of Case Western Reserve University. In the experiment, natural frequency of bearing used in this paper is FIn= 157.94 Hz, FOut= 104.56 Hz.

The total energy of the vibration signal of the bearing is estimated by the TEO technique. The time domain – amplitude waveforms are obtained as shown in Fig. 3 by inputting normal data, inner fault data and outer fault data to Teager program, respectively.

In order to obtain Teager spectrum, this paper uses FFT to transform TEO, the results are shown in Fig. 4. Resonance frequency band energy E1 is the energy from 2000 to 4000 Hz. Because energy leakage is inevitable. According to the practice that Δf= 2 Hz, this experiment chooses E(FIN±2), E(FOUT±2), and E1/ETOTAL as the feature vectors.

Fig. 3. Time domain waveform at outer ring fault condition, inner fault condition, normal condition

Time domain waveform at outer ring fault condition, inner fault condition, normal condition

a)

Time domain waveform at outer ring fault condition, inner fault condition, normal condition

b)

Time domain waveform at outer ring fault condition, inner fault condition, normal condition

c)

Fig. 4. Frequency spectrum at outer ring fault condition, inner fault condition, normal condition

Frequency spectrum at outer ring fault condition, inner fault condition, normal condition

a)

Frequency spectrum at outer ring fault condition, inner fault condition, normal condition

b)

Frequency spectrum at outer ring fault condition, inner fault condition, normal condition

c)

Table 1. Feature vectors

Inner ring (Hz)
Outer ring (Hz)
Ratio of energy (Hz)
Feature vectors
E ( F I N ± 2 )
E ( F O U T ± 2 )
E 1 / E T O T A L

Table 2. Bearing fault identification

E ( F I N ± 2 )
E ( F O U T ± 2 )
E 1 / E T O T A L
S V M 1
S V M 2
S V M 3
Outer
4.6848
0.8221
0.0567
1
Correct
Ring
4.6823
0.8210
0.0564
1
Correct
Fault
4.6817
0.8201
0.0565
1
Correct
Inner
0.3234
0.8723
0.1538
0
1
Correct
Ring
0.3239
0.8721
0.1538
0
1
Correct
Fault
0.3236
0.8724
0.1538
0
1
Correct
Normal
0.0005
0.0005
0.1759
0
0
1
Correct
0.0006
0.0005
01756
0
0
1
Correct
0.0007
0.0005
0.1758
0
0
1
Correct

3.2. Fault diagnosis based on SVM

This experiment adopts a multi-classifier composed of three SVMs. The experiment uses a total of 20×3 groups data (each mode is 20 groups) to classify the faults. Set the targets output as Eq. (1), (0, 1) and (0, 0, 1), corresponding to the outer ring fault, the inner ring fault, and the normal condition, respectively. Because the place is limited, only parts of data (3×3 groups) is showed in Table 2.

4. Conclusions

This paper proposes a fault diagnosis method for bearing based on TEO and SVM. Experiments have been conducted to verify the effectiveness of the proposed method. From the results it can be seen that: 1) TEO can be used to calculate the total energy of signal, and have a good time resolution. 2) SVM is suitable to diagnose and is extremely accurate as well. 3) The proposed method has a good adaptability, and high accuracy.

Acknowledgements

This research was supported by the National Natural Science Foundation of China (Grant No. 61074083, 50705005, and 51105019), as well as the Technology Foundation Program of National Defense (Grant No. Z132013B002).

References

  1. Zhao Xiao-Ling The vibrating diagnosis method for rolling bearing fault. Journal of Chongqing University of Science and Technology, Vol. 9, Issue 1, 2007, p. 41-44. [Search CrossRef]
  2. Li H. Order bi-spectrum for bearing fault monitoring and diagnosis under run-up condition. Journal of Computers, Vol. 6, Issue 9, 2011, p. 1994-2000. [Search CrossRef]
  3. Teager H. Some observations on oral air flow during phonation. IEEE Transactions on Acoustics Speech and Signal Processing, Vol. 28, Issue 5, 1980. [Search CrossRef]
  4. Kaiser J. F. On a simple algorithm to calculate the “energy” of a signal. Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing, Vol. 1, 1990, p. 381-384. [Search CrossRef]
  5. Wang Tian-Jin, Feng Zhi-Peng, Hao Ru-Jiang, Chu Fu-Lei Fault diagnosis of bearings using TEO. Journal of Vibration and Shock, Vol. 31, Issue 2, 2012. [Search CrossRef]
  6. Hou Guo-Ping, Ma Xuan, Tang Ming Application research of improved SVM model for intrusion detection. Machinery Design and Manufacture, Vol. 9, 2012. [Search CrossRef]
  7. Potamianos A., Maragos P. A comparison of the energy operator and the Hilbert transform approach to signal and speech demodulation. Signal Processing, Vol. 37, Issue 1, 1994, p. 95-120. [Search CrossRef]
  8. Vapnik V. N. Statistical Learning Theory. New York, Wiley, 1998. [Search CrossRef]
  9. Shen Zhi-Xi, Huang Xi-Yue, Ma Xiao-Xiao Fault diagnosis of diesel engine on EMD and SVM. Journal of Vibration, Measurement and Diagnosis, Vol. 30, Issue 1, 2010, p. 19-22, (in Chinese). [Search CrossRef]
  10. Khan A. F. Condition monitoring of rolling element bearing. A Comparative Study of Vibration Based Techniques, Canada, 1990. [Search CrossRef]
  11. Su Y. T., Sheen Y. T. On the detectability of roller bearing damage by frequency analysis. Proceedings of the Institution of Mechanical Engineers C, Vol. 207, Issue 1, 1993, p. 23-32. [Search CrossRef]