# Rolling Bearing Fault Detection Based on the Teager Energy Operator and Elman Neural Network

Authors: Hongmei Liu; Jing Wang; Chen Lu
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This paper presents an approach to bearing fault diagnosis based on the Teager energy operator (TEO) and Elman neural network. The TEO can estimate the total mechanical energy required to generate signals, thereby resulting in good time resolution and self-adaptability to transient signals. These attributes reflect the advantage of detecting signal impact characteristics. To detect the impact characteristics of the vibration signals of bearing faults, we used the TEO to extract the cyclical impact caused by bearing failure and applied the wavelet packet to reduce the noise of the Teager energy signal. This approach also enabled the extraction of bearing fault feature frequencies, which were identified using the fast Fourier transform of Teager energy. The feature frequencies of the inner and outer faults, as well as the ratio of resonance frequency band energy to total energy in the Teager spectrum, were extracted as feature vectors. In order to avoid a frequency leak error, the weighted Teager spectrum around the fault frequency was extracted as feature vector. These vectors were then used to train the Elman neural network and improve the robustness of the diagnostic algorithm. Experimental results indicate that the proposed approach effectively detects bearing faults under variable conditions.


---

## Body

## 1. Introduction

Bearing failure is one of the main causes of breakdown in rotating machinery, often resulting in worker injury or equipment damage, mission-critical failure, and costly system downtime. For these reasons, the detection of bearing faults has been the subject of extensive research. Vibration signals are usually collected to diagnose bearing faults. In practical industrial fields, however, the vibration signals of bearings are typically complex because of the effects of the working environment; detecting bearing faults only by observing the fast Fourier transform (FFT) of bearing measurement signals is difficult. Researchers have proposed numerous techniques for analyzing the frequency domain of vibration signals. FFT is generally used to generate vibration signal spectra, and the application of reasonable window functions and smoothing techniques can improve the quality of the frequency analysis of vibration signals [1, 2]. The primary purposes of such analyses are to identify the feature frequencies of bearing faults and use the parameters related to feature frequencies to diagnose bearing faults. To solve the interference of weak bearing fault signals during the extraction of vibration signal features from strong broadband noise, scholars proposed envelopment analysis; this method is extensively used, particularly for early fault detection [3–5]. Envelopment analysis can extract periodic pulses or amplitude-modulated signals from vibration signals. At increasingly severe bearing damage, however, bearing vibration mode becomes irregular and exhibits more random features. This irregularity makes envelopment analysis invalid for identifying feature frequency. Logarithmic cepstrum analysis was originally proposed for processing seismic data in an echo problem and was then gradually extended to the analysis of bearing vibration signals. Similar to envelopment analysis, the function of the logarithmic cepstrum in vibration signal analysis is to detect the harmonics and subband corresponding to the feature frequencies of bearing faults [6]. The difference is that envelopment analysis is applied to a particular band of a vibration signal, whereas the logarithmic cepstrum is used to analyze the vibration signals in the entire frequency domain. Nonetheless, the logarithmic cepstrum is rarely used in bearing fault detection possibly because such analysis requires the use of a computer and involves relatively complex interpretation of results. A wavelet transform is used in time-frequency analysis that has the characteristics of multiresolution analysis. Given that its time-frequency localization characteristics more accurately represent nonstationary signals, wavelet transforms have been frequently applied to bearing fault detection in the last 20 years [7]. A discrete wavelet transform is based on two subfrequency splits, which causes man-made damage to the feature frequency bands of bearing faults; in practical bearing fault detection, a continuous wavelet transform is typically used to analyze vibration signals. To a certain extent, this limits the application of wavelet transforms in systems that demand high computing speed and large storage capacity. Furthermore, similar to envelopment analysis, identifying the feature frequency of the bearing fault of a vibration spectrum becomes impossible at severe bearing damage.The Teager energy operator (TEO) has recently been applied in processing nonlinear signals. It is a type of nonlinear differential operator that estimates the total energy (kinetic energy and potential energy) produced by a signal source through the instantaneous value of a signal and a differential nonlinear combination. Therefore, the TEO can enhance the transient features of a signal, making it suitable for detecting impact signals. The TEO has a high time resolution and a good adaptability to the instantaneous changes in signals. Kaiser put forward and improved the TEO, as well as processed voice signals by applying this criterion [8]. Hongxing et al. verified the validity of the TEO by using a compressor vibration signal, comparing the performance of TEO demodulation with that of Hilbert envelope transform demodulation [9]. Bassiuny and Li used empirical mode decomposition to separate the amplitude modulation signal and frequency modulation signal of a machine tool feed motor and then applied the TEO demodulation method to calculate the instantaneous amplitude and instantaneous frequency of intrinsic mode functions; time-frequency information was constructed to monitor the condition of the machine tool [10]. To extract the cyclical shock feature of a bearing fault, this paper presents a TEO-based method for extracting the features of bearing vibration signals. The TEO showed the advantages of detecting the transient components of a signal, an approach that enhances the impact characteristics of the vibration signals of bearing faults. The faults are detected through feature frequencies obtained from the FFT of Teager energy.As depicted in Figure1, considering the advantage of the TEO in detecting transient impact, this study presents a bearing fault detection method based on this criterion. First, the TEO output of the original vibration signal of the bearing fault was calculated to analyze the cyclical characteristic of the transient impact. The instantaneous Teager energy sequence was denoised by wavelet packet and Fourier transformed to obtain the Teager spectrum. The fault feature frequencies differ under various faults, and for various types of bearing faults, the energies in the resonance frequency band also differ. The amplitude of the Teager spectrum around the fault frequencies of the inner and outer rings (within ±ΔfHz) and the ratio of resonance band energy to total energy were extracted as the feature vectors used to train the Elman neural network.Figure 1
The flow chart of the proposed method.In Table1, Fin represents the feature frequency of the inner ring, Fout denotes the feature frequency of the outer ring, E1 is the sum of Teager energy from 2000 Hz to 4000 Hz, and E0 represents the sum of total Teager energy.Table 1
Feature vectors.
Inner ring
Outer ring
Ratio of energy
Feature vector
F
in
±
Δ
f
F
out
±
Δ
f
E
1
/
E
0This paper is organized as follows. Section2 introduces the principle of TEO-based feature extraction. The Teager energy signal denoising by wavelet packet and fault classification by the Elman neural network are presented in Sections 3 and 4, respectively. The experimental results are reported in Section 5. Conclusions are drawn in the final section.


## 2. Principle of TEO-Based Feature Extraction

For any signalx(t), TEO ψ is defined as [11]
(1)ψ[x(t)]=[x˙(t)]2-x(t)x¨(t),
where x˙(t) and x¨(t) represent the first- and second-order time differential of x(t). The output of the TEO monitors the total energy required to produce a signal. In a linear undamped vibration system that consists of a mass and a spring, the second-order differential equation according to Newton’s laws of motion is
(2)mx¨(t)+kx(t)=0,
where x(t) represents the mass displacement relative to the equilibrium position, x¨(t) denotes the acceleration, k represents a stiffness coefficient of the spring, and m is the weight of the mass. This equation describes the movement of the mass, and its solution is a simple harmonic motion:
(3)x(t)=Acos(ωt+φ).Correspondingly, its velocity and acceleration are(4)x˙(t)=-Aωsin(ωt+φ),(5)x¨(t)=-Aω2cos(ωt+φ),
where A represents the amplitude of vibration and φ represents the initial phase, and  (k/m)1/2 represents the natural frequency.At any time, the mechanical energy of the vibration system is the sum of the potential energy of the spring and the kinetic energy of the mass:(6)E=12k[x(t)]2+12m[x˙(t)]2.Substituting (3) and (4) into (6) yields the total energy thus
(7)E=12mA2ω2.Equation (7) shows the instantaneous total energy of the harmonic vibration, and the square of the amplitude is inversely proportional to the square of the frequency.From the aforementioned equations, we can obtain that(8)ψ[x(t)]=ψ[Acos(ωt+φ)]=A2ω2.Comparing (7) and (8), the only difference is the constant m/2 between the output of the TEO and the instantaneous total energy of the simple harmonic motion.For discrete-time signalx(t), the difference equation is used instead of the differential equation:
(9)ψ[x(n)]=[x(n)]2-x(n-1)x(n+1).For the discrete-time signal, the TEO needs only three samples to calculate the energy of the signal at any time. Therefore, the TEO has a good time resolution for the instantaneous change in the signal and can detect the transient components of the signal.If the inner rings, outer rings, or rolling elements of rolling bearings have abrasions, cracks, pointed corrosions, or other local damages, the impact phenomenon caused by the mutations in vibration energy will occur when the damage point goes through the other components. This impact will intensify the resonance of the bearing components, and the resonance will rapidly increase. Therefore, the transient shock vibration caused by the bearing fault not only dramatically changes the amplitude of the vibration signal but also generates a higher frequency. Traditional signal energy is defined as the square of the signal amplitude, and it can represent only kinetic or potential energy. Although this method can highlight transient impact features, impact signals may be masked by other components if the impact is of small amplitude. The TEO can monitor the total energy of signals, including kinetic and potential energy, whose output is the square of the product of the instantaneous amplitude and the instantaneous frequency of the vibration. Compared with traditional signal energy, the TEO considers the frequency of a signal. Given the high frequency of transient impact vibration, the TEO can more effectively highlight the transient characteristics of the impact.In order to show the advantage of the TEO method, the Hilbert-Huang was used to compare.These two figures in Figures2 and 3 show the Teager and Hilbert spectrums can both extract the feature frequency from the rolling bearing signal when the bearing have the ordinary inner ring fault.Figure 2
Teager energy spectrum of the ordinary inner ring fault bearing.Figure 3
Hilbert-Huang spectrum of the ordinary inner ring fault bearing.However, if the rolling bearing works in harsh conditions, the TEO will show its advantage. To demonstrate it, the white noise with an SNR of 0.6 is introduced into the original vibration signal.As shown in Figures4 and 5, the TEO can successfully identify the feature frequency when the rolling bearing works in harsh conditions, whereas the Hilbert-Huang exhibits low performance in this regard.Figure 4
Teager energy spectrum of the inner ring fault bearing working in the harsh condition.Figure 5
Hilbert-Huang spectrum of the inner ring fault bearing working in the harsh condition.


## 3. Teager Energy Signal Denoising Based on Wavelet Packet

The collected vibration signals from rolling bearings running in complex conditions usually contain strong noises. In order to extract the fault features with high SNR, a wavelet packet was used to reduce noise existing in the Teager energy signal in this study.Traditional denoising methods, such as low-pass filter high-pass filter, can eliminate only low-frequency or high-frequency noise, but wavelet packet transforms, which have the function of multiresolution analysis, can analyze the signal both in the time domain and the frequency domain. Wavelet packet transforms can decompose signals into different bands and effectively distinguish the noise from the useful signal in different frequency bands, eliminating the noise of the signal in various frequency bands.The TEO can calculate the total energy of the sample points in original vibration signal. Because the original signal of a rolling bearing usually contains a strong background noise, the time domain waveform of Teager energy is still influenced by the energy of the noise. Wavelet packet-based denoising can effectively eliminate the random noise of the transient signal; therefore, it was used to reduce the noise in the Teager energy signal and extract the transient impact component.As shown in Figure6, the transient impact components can be highlighted with the effect of the noise reduced.Time domain waveform of Teager energy without and with wavelet packet-based denoising.
(a)
(b)


## 4. Fault Pattern Classification Based on the Elman Neural Network

For the extracted feature vectors by the TEO, neural networks are always used to classify the fault mode. In this study, the Elman neural network was used for fault classification.The Elman neural network includes input, output, hidden, and context layers. When the output of the hidden layer neurons goes through the context layer, it associates with the input of the hidden layer [12]. This type of association is sensitive to the historical state of data; internal feedback increases the capacity of the network to deal with dynamic information, overcoming the drawback of the feed-forward network, which does not have dynamic characteristics. Therefore, compared with other neural networks, the Elman neural network is more applicable to classifying the fault mode. The structure of the Elman neural network is shown in Figure 7.Figure 7
Elman neural network.The transfer function of the hidden layer in the Elman neural network is a sigmoid function. The activation function of the output layer neuron is a linear function. The nonlinear state equations of the network are(10)y(k)=g[w3x(k)],x(k)=f[w1xc(k)+w2u(k-1)],xc(k)=x(k-1),
where y represents the output vector of the network, u is the input vector of the network, x denotes the output vector of the hidden layer, xc represents the feedback state vector, w3 is the connection matrix from the context layer to the output layer, w2 is the connection matrix from the hidden layer to the context layer, w1 represents the connection layer from the input layer to the hidden layer, g(·) denotes the activity function of the output neurons, and f(·) represents the transfer function of the hidden layer.


## 5. Experimental Results Analysis

A cyclical shock will be generated when a bearing works. The frequency of this cyclical shock reflects the cause of the bearing fault. The feature frequency caused by the inner ring failure is(11)FIn=12fn(1+dD·cos∅)Zb.The feature frequency caused by the outer ring failure is(12)FOut=12fn(1-dD·cos∅)Zb.The feature frequency caused by the rolling element failure is expressed as follows:(13)FBall=12fn[1-(dD·cos∅)2]Dd,
where d represents the inner diameter, D represents the outer diameter, Zb denotes the number of rolling elements, and ∅ is the contact angle.All the rolling bearing data used in this study were acquired by the Bearing Data Center of Case Western Reserve University. The bearing used was a 6205-2RS JEM SKF deep groove ball bearing. The motor speeds of the training and test data were 1750 and 1772 rpm, respectively. The vibration data were collected using accelerometers. The digital data were collected at a speed of 12000 samples per second. The detailed bearing information is shown in Table2.Table 2
Bearing specifications.
(a)
Size: (in)
Inner diameter
Outer diameter
Thickness
Ball diameter
Pitch diameter
0.9843
2.0472
0.5906
0.3126
1.537
(b)
Defect frequencies: (multiples of running speed in Hz)
Inner ring
Outer ring
Cage train
Rolling element
5.4152
3.5848
0.39828
4.7135On the basis of (11), (12), and (13), the feature frequency can be calculated as FIn=157.94 Hz, FOut=104.56 Hz.The figures discussed in the succeeding statements show the vibration signals in the time domain, the frequency spectrum, and the Teager spectrum of each failure mode. Figures8, 9, and 10 illustrate the time domain waveforms of each fault. In the time domain waveforms, the cyclical shock features are not observable because the given signals are affected by noise. The fault modes cannot be distinguished from the time domain waveforms. Figures 11, 12, and 13 show the frequency spectrum of each fault. The feature frequencies of the fault are not easily distinguishable because of the resonance frequencies in the frequency spectrum. Figures 14, 15, and 16 indicate the Teager energy time domain waveform of each fault after wavelet packet-based denoising. The TEO highlights the transient impact feature, so that such features caused by the fault become more visible. The time interval between two adjacent impacts corresponds to fault feature frequencies.Figure 8
Time domain waveform of the vibration signal at normal conditions.Figure 9
Time domain waveform of the vibration signal at inner ring fault conditions.Figure 10
Time domain waveform of the vibration signal at outer ring fault conditions.Figure 11
Frequency spectrum of the vibration signal at normal conditions.Figure 12
Frequency spectrum of the vibration signal at inner ring fault conditions.Figure 13
Frequency spectrum of the vibration signal at outer ring fault conditions.Figure 14
Time domain waveform of Teager energy in a normal bearing after wavelet packet-based denoising.Figure 15
Time domain waveform of Teager energy at inner ring fault conditions after wavelet packet-based denoising.Figure 16
Time domain waveform of Teager energy at outer ring fault conditions after wavelet packet-based denosing.Figures17, 18, and 19 show the Teager spectrum of each fault. The fundamental and multiplier frequencies are clearly observable from the Teager spectrum.Figure 17
Teager spectrum of the normal bearing.Figure 18
Teager spectrum of the inner ring fault bearing.Figure 19
Teager spectrum of the outer ring fault bearing.The frequency spectrum of each fault mode shows that the resonance frequencies of the different fault modes vary confirming that the ratio of the energy of resonance frequency to total energy can be used as a feature vector.The amplitudes of the Teager spectrum in the inner and outer ring fault frequencies, as well as the ratio of the energy of resonance frequency to total energy, were extracted as the feature vectors, as listed in Table3. The weighted Teager spectrum in the fault frequency and around Δf=2 Hz was extracted as a feature vector to avoid a frequency leak error.Table 3
Feature vectors.
Inner ring (Hz)
Outer ring (Hz)
Ratio of energy
Feature vector
E
(
F
in
±
2
)
E
(
F
out
±
2
)
E
1
/
E
0To confirm the effectiveness of the proposed method under variable conditions, we used two types of bearing vibration signals with different motor speeds to detect the fault modes (inner ring fault, outer ring fault, or normal) in this study.A total of 60 groups of training data (normal data, inner ring fault data, and outer ring fault data; each mode with 20 groups) were used. The motor speed was 1750 rpm. The feature vectors were used to train the Elman neural network. According to the dimensionality of the feature vectors, the number of input nodes was set as three. The inputs of the Elman neural network were the fault feature vectors. The target outputs were set as (1, 0, 0), (0, 1, 0), and (0, 0, 1), representing the normal bearing, inner ring fault, and outer ring fault, respectively. The training data components are shown in Table4.Table 4
Training data.
Signal type(1750 rpm)
Target output
Input of the Elman neural network
Normal
(
1,0
,
0
)
0.0005
0.0005
0.1769
(
1,0
,
0
)
0.0006
0.0005
0.1760
(
1,0
,
0
)
0.0007
0.0005
0.1757
(
1,0
,
0
)
0.0007
0.0005
0.1765
(
1,0
,
0
)
0.0007
0.0006
0.1760
Inner ring fault
(
0,1
,
0
)
0.3274
0.8753
0.1548
(
0,1
,
0
)
0.3269
0.8751
0.1548
(
0,1
,
0
)
0.3267
0.8749
0.1548
(
0,1
,
0
)
0.3269
0.8750
0.1548
(
0,1
,
0
)
0.3192
0.8833
0.1551
Outer ring fault
(
0,0
,
1
)
4.6878
0.8272
0.0577
(
0,0
,
1
)
4.6851
0.8230
0.0576
(
0,0
,
1
)
4.6846
0.8216
0.0576
(
0,0
,
1
)
4.6821
0.8199
0.0576
(
0,0
,
1
)
4.6840
0.8188
0.0576To further validate the effectiveness of the proposed method under variable conditions, 15 groups of test data (normal data, inner ring fault data, and outer ring fault data; each mode with 5 groups) were collected under a different motor speed. These test data were collected at a motor speed of 1772 rpm. The experimental results are listed in Table5.Table 5
Test results.
Signal type(1750 rpm)
Expected output
Actual output
Normal
(
1,0
,
0
)
1.00
23
0.00
12
0.00
08
0.998
7
0.00
08
0.000
6
0.999
3
0.000
5
-
0.0001
0.99
74
0.000
7
0.0004
0.99
86
0.00
09
0.00
06
Inner ring fault
(
0,1
,
0
)
0.00
2
3
1.00
13
0.00
12
0.00
1
7
1.000
1
-
0.00
17
0.002
1
0.99
92
0.00
08
0.00
1
1
0.999
5
0.00
1
5
0.00
1
3
1.00
12
0.00
04
Outer ring fault
(
0,0
,
1
)
-
0.00
21
-
0.00
12
1.00
09
-
0.001
5
-
0.0002
1.00
17
0.0032
0.00
27
0.99
56
-
0.00
34
-
0.001
8
1.00
3
2
0.00
85
0.00
36
0.98
78
Signal type(1772 rpm)
Expected output
Actual output
Normal
(
1,0
,
0
)
0.9989
-
0.0054
-
0.0214
0.9987
-
0.0058
-
0.0243
1.0002
-
0.0073
-
0.0317
0.9973
-
0.0066
-
0.0293
0.9997
-
0.0068
-
0.0312
Inner ring fault
(
0,1
,
0
)
0
.
2
521
0.99
67
-
0.0009
0
.
2
546
0.9949
-
0.0031
0
.
2
524
0.99
89
-
0.0015
0
.
2
521
0.999
5
-
0.0014
0
.
2
523
1.00
12
-
0.0014
Outer ring fault
(
0,0
,
1
)
0.0769
-
0.1424
1.00
02
0.0774
-
0.1429
0.9992
0.0777
-
0.1430
0.99
67
0.0779
-
0.1438
1.00
01
0.0781
-
0.1444
0.9
976
Signal type(1721 rpm)
Expected output
Actual output
Normal
(
1,0
,
0
)
0.9989
-
0.1411
0.2092
0.9962
-
0.1401
0.2066
0.9961
-
0.1330
0.1877
0.9961
-
0.1334
0.1886
0.9961
-
0.1347
0.1913
Inner ring fault
(
0,1
,
0
)
0.0082
0.99
43
0.1207
0.0050
0.9949
0.1206
0.0054
0.99
39
0.1199
0.0035
0.99
44
0.1200
-
0.0001
0.9951
0.1230
Outer ring fault
(
0,0
,
1
)
0.0608
0.0854
0.9963
0.0609
0.0855
0.9961
0.0614
0.0859
0.99
67
0.0617
0.0871
0.9972
0.0615
0.0862
0.9
970To verify the robustness and applicability of the proposed method, three Teager spectrums of the inner ring fault bearing working in different conditions are shown in Figure20.Teager spectrums of the inner ring fault bearing working in different conditions.
(a)
1750 r/min
(b)
1772 r/min
(c)
1721 r/minWhen the speed of the bearing increases from 1721 r/min to 1772 r/min, the fault frequency correspondingly changes from 156 Hz to 159 Hz. In this study, the weighted Teager spectrum in the fault frequency and around (±2 Hz) was an extracted feature vector. It cannot only avoid a frequency leak error but also cover the fault frequencies working in different conditions. Thus, the robustness of the proposed method is improved and suitable for different conditions.The TEO can extract the transient impact feature of the faulty bearing, and the Elman neural network can classify the bearing failure mode. In addition, because the motor speed of the test data was different from that of the training data, the method exhibited good robustness, indicating applicability to different operating conditions.


## 6. Conclusion

The TEO can calculate the instantaneous total energy of signals; thus, it is suitable for detecting transient signals and highlights the impact characteristics of transient peak signals. With consideration for the impact characteristics of the vibration signal of a faulty bearing, this paper presents a TEO and Elman neural network-based method for detecting bearing faults. The TEO is advantageous for detecting the impact characteristics of a signal and can therefore effectively detect the impact signal of a faulty bearing. Moreover, the TEO exhibits a good time resolution and self-adaptability to transient signals, showing better robustness in bearing fault detection.



---

*Source: 498385-2013-05-20.xml*