Next Article in Journal
Identification of Pulmonary Hypertension Using Entropy Measure Analysis of Heart Sound Signal
Previous Article in Journal
Research on Weak Fault Extraction Method for Alleviating the Mode Mixing of LMD
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Teager Energy Entropy Ratio of Wavelet Packet Transform and Its Application in Bearing Fault Diagnosis

Department of Mechanical Engineering, North China Electric Power University, Baoding 071003, China
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(5), 388; https://doi.org/10.3390/e20050388
Submission received: 27 April 2018 / Revised: 15 May 2018 / Accepted: 18 May 2018 / Published: 21 May 2018

Abstract

:
Kurtogram can adaptively select the resonant frequency band, and then the characteristic fault frequency can be obtained by analyzing the selected band. However, the kurtogram is easily affected by random impulses and noise. In recent years, improvements to kurtogram have been concentrated on two aspects: (a) the decomposition method of the frequency band; and (b) the selection index of the optimal frequency band. In this article, a new method called Teager Energy Entropy Ratio Gram (TEERgram) is proposed. The TEER algorithm takes the wavelet packet transform (WPT) as the signal frequency band decomposition method, which can adaptively segment the frequency band and control the noise. At the same time, Teager Energy Entropy Ratio (TEER) is proposed as a computing index for wavelet packet subbands. WPT has better decomposition properties than traditional finite impulse response (FIR) filtering and Fourier decomposition in the kurtogram algorithm. At the same time, TEER has better performance than the envelope spectrum or even the square envelope spectrum. Therefore, the TEERgram method can accurately identify the resonant frequency band under strong background noise. The effectiveness of the proposed method is verified by simulation and experimental analysis.

1. Introduction

The rolling bearing is one of the most common components in rotating machinery, and it is also one of the components with the highest failure rate. The running condition of bearings is directly related to the safe and stable operation of the whole rotating machinery system. Extracting fault feature information from bearing vibration signals under strong background noise is a hot issue in the field of mechanical fault diagnosis. The process of bearing fault diagnosis is divided into two parts: vibration parameter (displacement, velocity and acceleration) acquisition and vibration data processing. Adamczak et al. analyzed the existing bearing fault testing system. They strived for more accurate bearing fault measurement equipment and fault measurement parameters by comparing and analyzing existing measuring devices [1]. Meanwhile, a large number of vibration signal processing methods have been put forward, among which the kurtogram method is able to adaptively recognize the resonance frequency band of the vibration signal. By filtering the resonance frequency band and integrating the spectrum analysis of filtered signals, the fault characteristics of vibration signals can be effectively extracted.
Kurtosis is the 4th-order cumulant of random variables. This index cannot reflect the change of specific signals, and is not suitable for condition monitoring under strong noise environments. In order to overcome the shortage of kurtosis in engineering applications, Antoni made a formal definition of kurtosis [2]. Subsequently, Antoni further proposed the kurtogram. The core idea of the algorithm is that the spectral kurtosis index of the unstable signal is maximized by reasonable selection of frequency resolution [3]. In order to reduce the operation time of the kurtogram and enable it to be widely used in engineering practice, Antoni further proposed a fast spectral kurtosis (FSK) calculation method based on a fast algorithm for computing the kurtogram over a grid. FSK can simplify the process of selecting the optimal band ( f , Δ f ) [4]. In the FSK algorithm, the horizontal axis represents the frequency, the vertical axis represents the number of decomposition layers, and the depth of the color indicates the spectral kurtosis value of each sub band.
Since then, many algorithms based on the kurtogram have been proposed. Lei et al. proposed an improved kurtogram method [5], which uses wavelet packet transform (WPT) instead of the short time Fourier transform (STFT) and FIR filters used in the traditional kurtogram. A filter based on WPT can accurately divide the frequency bands and control the noise effectively [5]. Then, on the basis of Lei et al., Wang et al. made further improvements. The kurtosis calculation method in the traditional kurtogram was changed to calculate the power spectrum of the envelope of different nodes of WPT. The power spectrum of the envelope signal reflects the sparsity of the signal and helps to capture the resonance band in the signal [6]. Based on the same idea, Tse and Wang proposed a new method called sparsogram. The power spectrum of the envelope was applied to improve the kurtosis index in the kurtogram, so that it reflected the sparsity of signals [6,7]. Gu proposed a method for analyzing kurtograms in the frequency domain [8]. The correlated kurtosis of the envelope was used to replace the traditional calculating method in the kurtogram, which makes the calculation index sensitive to periodic impulses and has better robustness. To obtain the correct relationship between the node and frequency band in WPT, a vital process called frequency ordering is conducted to solve the frequency folding problem due to down-sampling. Tian et al. proposed the modulation signal bispectrum (MSB)-based robust detector to reduce the noise and accurately identify the resonance frequency band [9]. The high-magnitude features that result from the use of the MSB also enhance the modulation effects of a bearing fault. In gearboxes, the kurtogram method would fail when gears and bearings partially malfunctioned at the same time. In view of this, Wang proposed an improved method by combining kurtogram and menshing resonance (MRgram) [10]. Random impulses often occur in the frequency bands of vibration signals of bearings, making the selection of kurtogram resonance frequency bands inaccurate. Xu et al. proposed a periodicity-based kurtogram to deal with random impulse resistance [11]. The periodic component to aperiodic component ratio (PAR) was utilized in this method to differentiate the types of impulses. Li et al. improved the kurtogram based on an impulse step dictionary and a reweighted minimizing nonconvex penalty Lq regular for rolling bearing fault diagnosis [12].
Entropy was originally a thermodynamic concept describing the regularity of information. Because entropy has some statistical properties superior to kurtosis in some respects, the application scope of entropy was gradually extended to the field of mechanical fault diagnosis. Antoni introduced Shannon entropy into the kurtogram algorithm and proposed infogram [13]. Infogram contains three parts, the envelope (SE) infogram, squared envelope spectrum (SES) infogram, and SE1/2/SES1/2 infogram. The calculation index of the SE infogram is the negative entropy of the square envelope, and the calculation index of the SES infogram is the negative entropy of the square envelope spectrum. On the basis of infogram, the multiscale clustering grey infogram (MCGI) was proposed by Li et al. [14]. Hemmati et al. proposed an index that combines kurtosis and Shannon entropy, and the index was used to obtain the optimal band pass filter utilizing wavelet packet transform (WPT) and envelope detection [15]. The MCGI combined the negentropies of the time and frequency domains in a grey fashion using multiscale clustering. Chen et al. proposed a fault classification method combining time-spectral kurtosis, entropy and support vector machines (SVM). Combining the advantages of time-spectral kurtosis (T-SK) and entropy, this method makes fault extraction more obvious [16].
To sum up, the improvement of kurtogram has mainly focused on two aspects. On the one hand, the decomposition methods (STFT, FIR) of the traditional kurtogram algorithm have been improved, so that the frequency band segmentation will have better results. On the other hand, the calculation index of each sub-band has been improved. The aim is to make the calculation index sensitive to the periodic impact caused by the fault, and to have strong robustness to the interference caused by noise or accidental impact. These improvements have made a significant contribution in some respects. In order to further improve the accuracy and efficiency of the algorithm in determining the center frequency and bandwidth, a new adaptive method called TEERgram is proposed in this article.

2. The Index of Teager Energy Entropy

In 2016, Antoni proposed the infogram method. Infogram theories include the squared envelope (SE) infogram, the squared envelope spectrum (SES) infogram, and the SE1/2/SES1/2 infogram. The SE and SES were defined as Equations (1) and (2).
Δ I e = SE ( n ; f , Δ f ) 2 SE ( n ; f , Δ f ) 2 In ( SE ( n ; f , Δ f ) 2 SE ( n ; f , Δ f ) 2 )
Δ I E = SES ( n ; f , Δ f ) 2 SES ( n ; f , Δ f ) 2 In ( SES ( n ; f , Δ f ) 2 SES ( n ; f , Δ f ) 2 )
where n Z , f is the frequency, and Δ f is the frequency resolution. Δ I e and Δ I E are entropy values in the time domain and entropy in the frequency domain respectively. That is, they were used to represent the impact characteristics and cyclostationarity characteristics of fault signals, respectively. Antoni further proposed the SE1/2/SES1/2 infogram to characterize the impact characteristics and cyclostationarity. The SE1/2/SES1/2 infogram was defined by Equation (3):
Δ I 1 / 2 = Δ I e / 2 + Δ I E / 2

2.1. Teager Energy Operator

Teager energy operator (TEO) can enhance transient impact components, which is suitable for detecting impact characteristics in signal, and has good effect in fault impact feature extraction. For signal x ( t ) , the Teager energy operator ψ can be defined as Equation (4):
ψ ( x ( t ) ) = ( x ( t ) ) 2 x ( t ) x ( t )
where x ( t ) and x ( t ) are the first- and second-order derivatives of signal x relative to time t. The general expression of the amplitude modulation frequency modulation signal is shown in Equation (5):
s ( t ) = a ( t ) cos ( ψ ( t ) )
where a(t) is the instantaneous amplitude of s ( t ) and ψ ( t ) is the instantaneous phase of s ( t ) . The instantaneous frequency of the signal s ( t ) is ω = d ψ / d t . The effect of the energy operator on s ( t ) is shown in Equation (6):
ψ ( s ( t ) ) = ( a ( t ) ψ ( t ) ) / 2 + a 2 ( t ) ψ ( t ) × sin ( 2 ψ ( t ) ) / 2 + cos 2 ( ψ ( t ) ) ψ ( a ( t ) )
Because the change of the carrier signal is much faster than the modulated signal, the instantaneous amplitude and the instantaneous frequency of the modulated signal are relatively slow relative to the high-frequency carrier, and can be approximately equal to the constant. Therefore, if ψ ( t ) 0 , ψ ( a ( t ) ) 0 , then the Formula (6) can be approximately expressed as Equation (7).
ψ ( s ( t ) ) ( a ( t ) ψ ( t ) ) 2 = a 2 ( t ) ω 2 ( t )
Meanwhile,
ψ ( s ( t ) ) a 2 ( t ) ω 4 ( t )
Therefore, the envelope signal (instantaneous amplitude) | a ( t ) | obtained by demodulation and separation of signal s can be expressed as Equation (9):
| a ( t ) | ψ ( s ( t ) ) ψ ( s ( t ) )
The signal energy in transmission is defined as the square of the amplitude of the signal. If the impact amplitude is small, the impact component may be drowned by other components. The Teager energy operator is the product of squared instantaneous amplitude and instantaneous frequency. Due to the high vibration frequency of transient impact, the Teager energy operator can effectively enhance transient impact components and significantly suppress noise.
In order to verify the superiority of the Teager energy operator spectrum (TEOS) compared with the envelope spectrum (ES) and the square envelope spectrum (SES), the following simulation signals are defined as Equation (10).
y 1 = i = 1 n A e ξ ( t q i ( t ) / f outer ) 2 × sin ( 2 π f o t ) + n ( t )
where ξ is the damping ratio, A is the amplitude, q i ( t ) = t × f i ( i = 1 , 2 , , n ), f o is the natural frequency, f outer = 25 Hz is the characteristic frequency and n ( t ) is the Gauss white noise. The parameters of the simulation signal are shown in Table 1.
The waveforms of the synthetic signal with fault impulses and noise are shown in Figure 1, Figure 1a shows a periodic signal, which is used to simulate the fault impulses. Figure 1b shows a synthetic signal with fault impulses and noise (–12 dB). The envelope spectrum (ES), the square envelope spectrum (SES) and the Teager energy operator spectrum (TEOS) of the simulation signal are analyzed, respectively. The results are shown in Figure 2. It can be seen from Figure 2 that the amplitude of TEOS is higher than the amplitude of ES and the amplitude of SES at the fault characteristic frequency spectrum line. However, in the rest of the spectrum, the result is just the opposite. That is to say, the signal-to-noise ratio (SNR) of the Teager energy operator spectrum is much higher than that of the square envelope spectrum and the envelope spectrum.

2.2. Shannon Entropy

In 1948, Shannon put forward the concept of “Shannon entropy” and solved the problem of quantitative measurement of information. There is a direct relationship between the amount of information in a signal and its uncertainty, so from this point of view, we can think that the measure of the amount of information is equal to how much uncertainty there is. Shannon entropy reflects the degree of disordering (ordering) of a system, the more orderly a system is, the lower the entropy, and vice versa. In view of this, Shannon entropy can reflect the order of vibration signals. Given a random sequence { x 1 , x 2 , , x n } , the calculation method of Shannon entropy is shown in Equation (11):
S h = c i = 1 n p ( x i ) × log ( p ( x i ) )
where S h is the value of Shannon entropy, p ( x i ) is the probability mass associated with the value x i and c is an arbitrary positive constant that dictates the units.
The Shannon entropy index can reflect the regularity of the vibration signal, and its value has a certain relation with the SNR, so it can be used as an indicator to reflect the intensity of the impact of the fault. To verify this conclusion, the simulation signal model in Equation (10) is added with different degrees of noise to simulate the numerical changes of Shannon under different SNR. It can be seen from Figure 3 that the numerical regularity of Shannon entropy decreases with the increase of SNR.

2.3. Teager Energy Entropy

In view of the superiority of the Teager energy operator relative to the square envelope spectrum, we propose a Teager Energy Entropy (TEE) index. Inspired by infogram, we define the Teager energy entropy in the time domain (TEEt) and the Teager energy entropy in the frequency domain (TEEf) respectively as Equations (12) and (13).
Δ I TEEt = TEO t ( n ; f , Δ f ) 2 TEO t ( n ; f , Δ f ) 2 In ( TEO t ( n ; f , Δ f ) 2 TEO t ( n ; f , Δ f ) 2 )
Δ I TEEf = TEO f ( n ; f , Δ f ) 2 TEO f ( n ; f , Δ f ) 2 In ( TEO f ( n ; f , Δ f ) 2 TEO f ( n ; f , Δ f ) 2 )
In order to optimize the proportion of time domain and frequency domain components, an adaptive weighting method is proposed. Firstly, the regularization of TEEt and TEEf is carried out using Equations (14) and (15).
T i t = Δ I TEEt i min i = 1 N Δ I TEEt i max i = 1 N Δ I TEEt i min i = 1 N Δ I TEEt i
T i f = Δ I TEEf i min i = 1 N Δ I TEEf i max i = 1 N Δ I TEEf i min i = 1 N Δ I TEEf i
In these, m = 1, 2, …, M (n = 1, 2, …, N) represents the M(N)-th signal component. Then the mean value ( μ ) and standard deviation ( σ ) of the two regularization sequences are calculated separately using Equations (16)–(19).
μ t = 1 M i = 1 M T i t
μ f = 1 N i = 1 N T i f
σ t = 1 M i = 1 M ( T i t μ t ) 2
σ f = 1 N i = 1 N ( T i f μ f ) 2
Then the weight coefficients of the regularized sequence are calculated as Equations (20) and (21).
C t = 1 σ t 1 σ t + 1 σ f
C f = 1 σ f 1 σ t + 1 σ f
The Teager energy entropy, combined with the time domain and the frequency domain, is finally expressed as Equation (22):
TEE = C t Δ I TEEt + C f Δ I TEEf

2.4. Teager Energy Entropy Ratio of Wavelet Packet

Lei et al. introduced wavelet packet transform (WPT) into the kurtogram algorithm to replace the traditional decomposition method (FIR, STFT) in kurtogram. Meanwhile, Stępień et al. confirmed that the entropy-based WPT is the most suitable approach that is not limited to one level of decomposition, and it allows us to find an optimal decomposition tree [17]. Inspired by this idea, the algorithm is further improved by the TEE index in this paper. The paving of the kurtogram and the TEEgram are shown in Figure 4.
In many cases, the components that reflect the characteristics of fault impact are likely to be drowned out by noise. This can easily cause misdiagnosis during the confirmation of the resonance frequency band. To illustrate these problems, different degrees of noise are added to the simulated signals in Equation (10). Figure 5a shows the simulated fault signal with SNR = −12 dB, Figure 5b shows the simulated fault signal with SNR = −16 dB, Figure 5c shows the TEEgram of (a), while Figure 5d shows the TEEgram of (b). When the signal-to-noise ratio (SNR) is −12 dB, the TEEgram can accurately locate the resonance frequency band (Table 1; the natural frequency is about 1500 Hz), but when the SNR is reduced to −16 dB, the resonance frequency band is inaccurate.
In view of the defect that kurtogram is too sensitive to strong background noise and singular points, Wang et al. proposed an improved method named SKRgram [18]. Subsequently, Miao et al. used the Gini coefficient as an index to further improve the SKRgram [19]. The essence of the SKRgram is to compare the change of spectral kurtosis before and after failure. Based on this, the potential bands of fault impact hidden in noise and other constant interference components can be found. Inspired by this idea, we further modify the TEE index as Equation (23):
TEER = TEE f a i l u r e / TEE h e a l t h

2.5. Process of the Proposed Method

  • The bearing vibration signals are collected under healthy state and fault state conditions (Figure 6a,b). According to the bearing parameters, the characteristic frequency of bearing fault is calculated based on the theoretical equations (Equation (24)).
    f i = Z 2 ( 1 + d D cos β ) N 60 ,   f o = Z 2 ( 1 d D cos β ) N 60 f ball = D 2 d ( 1 ( d D ) 2 cos 2 β ) N 60 ,   f cage = 1 2 ( 1 d D cos β ) N 60
    where f i , f o , f ball and f cage respectively represent inner ring fault, outer ring fault, rolling element fault and cage fault. d is the ball diameter, D is the pitch diameter, β is the contact angle, N is the speed of rotation, Z is the number of the balls.
  • Wavelet packet transform (WPT) is applied, in turn, to vibration signals in healthy state and fault state. The time domain Teager energy entropy ( Δ I TEEt ) and the frequency domain Teager energy entropy ( Δ I TEEf ) are calculated for each sub-band of the signal under the healthy condition. Δ I TEEt and Δ I TEEf are weighted to obtain Teager energy entropy ( TEE h e a l t h = C t Δ I TEEt + C f Δ I TEEt ) in the healthy state. Similarly, the energy entropy of Teager ( TEE f a i l u r e ) under the fault condition is obtained. Teager energy entropy diagram (TEEgram) of wavelet sub-band of health status and fault state signals is constructed respectively (Figure 6c,d).
  • According to the Teager energy entropy of the signal in the healthy state ( TEE h e a l t h ) and the Teager energy entropy of the signal in the fault state ( TEE f a i l u r e ), the Teager energy entropy rate ( TEER = TEE f a i l u r e / TEE h e a l t h ) is calculated. Teager energy entropy ratio diagram (TEERgram) of wavelet sub-band is constructed as shown in Figure 6e. It can be seen from Figure 6e that the location of the relative change of TEE value is the fault resonance frequency band. The wavelet subband with the largest TEER value is selected as the optimal band. The center frequency and bandwidth are extracted to set the parameters of the bandpass filter, and the filter is used to filter the selected wavelet sub-bands. Then the envelope spectrum of the filtered signal is analyzed (Figure 6f). The fault type is analyzed by comparing the characteristic frequency of the spectrum line with the theoretical characteristic frequency calculated in step 1.

3. Simulation Analysis

The effectiveness of the proposed algorithm is verified by the mathematical model of the inner ring fault simulation signal in Equation (25) [20].
{ x ( t ) = s ( t ) + n ( t ) = i A i h ( t i T τ i ) + n ( t ) A i = A 0 cos ( 2 π f r t + ϕ A ) h ( t ) = e B t cos ( 2 π f n t + ϕ ω )
where τ i is the ith shock; T is the period of shock; f r is the rotating frequency; A i is the amplitude modulation with 1 / f r as a cycle; h ( t ) is an exponential decay pulse; B is the attenuation coefficient of vibration caused by damping; and ϕ A , ϕ ω represent the phase. The other parameters of the simulation signal are shown in Table 2.
Figure 7a shows the waveform of the inner fault impulses; Figure 7b shows the noise; Figure 7c shows the synthetic signal with fault impulses and noise (−16 dB). Figure 7d shows the spectrum of Figure 7c. It can be seen from Figure 7d that the fault feature information cannot be found, and there is a suspected resonance frequency band near 2500 Hz.
The simulation signal is respectively analyzed by kurtogram, kurtosis of the wavelet packet transform (KWPT) and the proposed method. The results of kurtogram analysis are shown in Figure 8. The results of KWPT analysis are shown in Figure 9. The two methods are inaccurate in locating the resonant frequency band, and are significantly affected by noise. In the envelope spectrum of their sub-bands, the characteristic frequencies of faults can be found. Figure 10 is the method proposed in this paper. By calculating the Teager energy operator ratio (TEER) of the signal before and after the fault, the resonant frequency band that changes before and after the fault is obtained, that is, the resonance frequency band caused by the fault.

4. Experiment Analysis

4.1. Experiment 1

Data from the bearing data center, Case Western Reserve University, are widely used in the field of fault diagnosis as common open data, and have a high degree of recognition (Figure 11) [21]. The experimental data with the weakest degree (0.007 inch) of failure are used to verify the effectiveness of the proposed method in experimental applications. The sampling frequency is 12 kHz, and the rotation speed is 1772 r/min. The bearing parameter dimensions (JEMSKF6023-2RS) are shown in Table 3, and the characteristic frequency of the bearing fault calculated by the theoretical formula is shown in Table 4.
8192 points of vibration signal of rolling element fault (Figure 12) are analyzed. The vibration signal of rolling element fault is respectively analyzed by kurtogram, KWPT and the proposed method. The results of kurtogram analysis are shown in Figure 13. The results of KWPT analysis are shown in Figure 14. The two methods are inaccurate for locating the resonant frequency band, and are significantly affected by noise. In the envelope spectrum of their sub-bands, the characteristic frequencies of faults can be found. Figure 15 shows the method proposed in this paper. The center of the resonance band determined by the TEEgram is about 3500 Hz. The partial band is filtered, and then the envelope spectrum of the filtered signal is analyzed. It can be seen from the envelope spectrum (Figure 15c) that the fault characteristic frequency and its first harmonic are effectively extracted.

4.2. Experiment 2

Experiment 2 is designed on the experiment platform to further verify the robustness of the proposed method. Experiment 2 contains two parts, which are the outer ring fault experiment and the inner and outer ring compound fault experiment. The bearing failure experiments are carried out on the QPZZ-II rotating machinery failure experiment platform. The structure of the experiment platform is shown in Figure 16. The type of bearing used in the experiment is LYC6205E (Figure 17). The parameters of the rolling bearing are shown in Table 5. The spark erosion technique is used to process faults in inner rings and outer rings of bearings. The degrees of failure are 1.5 mm deep and 0.2 mm wide. The acceleration sensor arranged on the bearing seat is used to collect vibration signals generated during the rotation of the bearing. The sampling frequency is 12,800 Hz, and the motor speed is 1466 r/min.
The characteristic frequency of inner ring and outer ring fault can be calculated from the following equations:
f o = Z 2 ( 1 d D cos α ) N 60 = 87.7 Hz ,   f i = Z 2 ( 1 + d D cos β ) N 60 = 132 . 2 Hz

4.2.1. Case 1: Analysis of Outer Ring Fault

8192 points of vibration signal of outer ring fault (Figure 17a) are analyzed, and the time domain waveform of the outer ring fault is shown in Figure 18.
The outer ring fault signal is analyzed via the kurtogram (Figure 19). As shown in Figure 19a, the central frequency of the resonant band determined by the kurtogram method is 860 Hz. The band pass filtering is applied to the resonant frequency band and envelope spectrum analysis of filtered signals is carried out. As shown in Figure 14b, fault feature information is not found in the envelope spectrum. This means that the resonant frequency identified by kurtogram does not contain fault feature information; that is, kurtogram is not accurate at locating the resonant frequency band. Similarly, the resonant frequency band in Figure 20a is filtered and envelope spectrum analysis is performed, and the fault characteristic spectrum is not found in Figure 20b. This indicates that the KWPT method is also invalid. Finally, the vibration signal is analyzed by the proposed method. The center of the resonant frequency band is about 2600 Hz from the TEERgram shown in Figure 21a. The envelope spectrum analysis of the band is shown in Figure 21c. The characteristic frequency and its harmonics of the bearing outer ring fault can be found.

4.2.2. Case 2: Analysis of Compound Fault in Inner and Outer Ring

8192 points of vibration signal of inner and outer ring fault are analyzed (Figure 17b), and the time domain waveform is shown in Figure 22. Figure 23 shows the processing result of KWPT. It can be seen that due to noise interference, the sub-band determined by KWPT is a sub-band containing more noise, rather than a sub-band containing more fault feature information. The sub-band containing the characteristic information of the outer ring fault is not obvious, and the sub-band containing the fault information of the inner ring is completely submerged. The kurtogram analysis of compound fault vibration signal is shown in Figure 24. It can be seen from Figure 24a that the resonant frequency band of the outer ring fault is found, but the resonant frequency band of the relatively weak inner ring fault is submerged. The resonance frequency band is filtered, and the characteristic frequency and frequency doubling of the outer ring fault can be found in the envelope spectrum shown in Figure 24b, but there is no characteristic information of the inner ring fault. Then, the composite fault signal is analyzed by the method proposed in this paper. With the benefits of the relative change rate of the resonance band before and after failure, the main frequency band and the secondary frequency band can be found in the TEERgram in Figure 25a. The two frequency bands are filtered respectively, and the envelope spectrum analysis of filtered signals are shown in Figure 25b. The fault characteristic frequency and its harmonics of inner ring and outer ring fault are separated and extracted. The advantage of TEERgram is that the ratio of Teager energy entropy before and after fault is calculated as an index. This ratio reflects the relative change value caused by faults, not the size of the absolute quantity; that is, the change caused by the failure. This will avoid the influence of the absolute amount caused by random impulses and noise to the analysis process, and to a certain extent, it avoids the influence of stronger faults on the weaker faults. Therefore, we think that this method is effective for single faults submerged in random impulses and noise, and also has certain effect for some compound faults.

5. Conclusions

In this paper, the TEERgram method is proposed. On one hand, WPT is used to replace the STFT and FIR filters in traditional kurtogram. WPT can divide the frequency band meticulously and suppress noise at the same time. On the other hand, the TEER is used as the computing index of the wavelet sub-band to find the optimal resonance frequency band. The time domain entropy ( Δ I TEEt ) and frequency domain entropy ( Δ I TEEf ) of wavelet sub-bands are processed statistically so that the TEE index can reflect the impact characteristics and periodic characteristics adaptively. By calculating the change of the TEER value before and after the fault, it can highlight the effect of the fault. Finally, simulations and experiments are carried out to demonstrate the effectiveness of the proposed method and its superiority over kurtogram.

Author Contributions

S.W. directed the experimental analysis and paper writing; X.Z. performed the experiments; X.Z. analyzed the data; X.Z. contributed reagents/materials/analysis tools and wrote the paper.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities (No.2017XS133) and the National Natural Science Foundation of China (No. 51777075).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adamczak, S.; Stępień, K.; Wrzochal, M. Comparative Study of Measurement Systems Used to Evaluate Vibrations of Rolling Bearings. Procedia Eng. 2017, 192, 971–975. [Google Scholar] [CrossRef]
  2. Antoni, J. The spectral kurtosis: A useful tool for characterising non-stationary signals. Mech. Syst. Signal Process. 2006, 20, 282–307. [Google Scholar] [CrossRef]
  3. Antoni, J.; Randall, R.B. The spectral kurtosis: Application to the vibratory surveillance and diagnostics of rotating machines. Mech. Syst. Signal Process. 2006, 20, 308–331. [Google Scholar] [CrossRef]
  4. Antoni, J. Fast computation of the kurtogram for the detection of transient faults. Mech. Syst. Signal Process. 2007, 21, 108–124. [Google Scholar] [CrossRef]
  5. Lei, Y.; Lin, J.; He, Z.; Zi, Y. Application of an improved kurtogram method for fault diagnosis of rolling element bearings. Mech. Syst. Signal Process. 2011, 25, 1738–1749. [Google Scholar] [CrossRef]
  6. Wang, D.; Tse, P.W.; Tsui, K.L. An enhanced Kurtogram method for fault diagnosis of rolling element bearings. Mech. Syst. Signal Process. 2013, 35, 176–199. [Google Scholar] [CrossRef]
  7. Tse, P.W.; Wang, D. The automatic selection of an optimal wavelet filter and its enhancement by the new sparsogram for bearing fault detection. Mech. Syst. Signal Process. 2013, 40, 520–544. [Google Scholar] [CrossRef]
  8. Gu, X.; Yang, S.; Liu, Y.; Hao, R. Rolling element bearing faults diagnosis based on kurtogram and frequency domain correlated kurtosis. Meas. Sci. Technol. 2016, 27, 125019. [Google Scholar] [CrossRef]
  9. Tian, X.; Gu, J.X.; Rehab, I.; Abdalla, G.M.; Gu, F.; Ball, A.D. A robust detector for rolling element bearing condition monitoring based on the modulation signal bispectrum and its performance evaluation against the Kurtogram. Mech. Syst. Signal Process. 2018, 100, 167–187. [Google Scholar] [CrossRef]
  10. Wang, T.; Chu, F.; Han, Q.; Kong, Y. Compound faults detection in gearbox via meshing resonance and spectral kurtosis methods. J. Sound Vib. 2017, 392, 367–381. [Google Scholar] [CrossRef]
  11. Xu, X.; Zhao, M.; Lin, J.; Lei, Y. Periodicity-based kurtogram for random impulse resistance. Meas. Sci. Technol. 2015, 26, 085011. [Google Scholar] [CrossRef]
  12. Li, Q.; Liang, S.Y. Incipient Fault Diagnosis of Rolling Bearings Based on Impulse-Step Impact Dictionary and Re-Weighted Minimizing Nonconvex Penalty Lq Regular Technique. Entropy 2017, 19, 421. [Google Scholar] [CrossRef]
  13. Antoni, J. The infogram: Entropic evidence of the signature of repetitive transients. Mech. Syst. Signal Proc. 2016, 74, 73–94. [Google Scholar] [CrossRef]
  14. Li, C.; Cabrera, D.; de Oliveira, J.V.; Sanchez, R.-V.; Cerrada, M.; Zurita, G. Extracting repetitive transients for rotating machinery diagnosis using multiscale clustered grey infogram. Mech. Syst. Signal Process. 2016, 76–77, 157–173. [Google Scholar] [CrossRef]
  15. 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]
  16. Chen, B.; Yan, Z.; Chen, W. Defect Detection for Wheel-Bearings with Time-Spectral Kurtosis and Entropy. Entropy 2014, 16, 607–626. [Google Scholar] [CrossRef]
  17. Stępień, K. Defining the criteria to select the wavelet type for the assessment of surface quality. Tehn. Vjesn. Tech. Gazet. 2015, 22, 781–784. [Google Scholar] [CrossRef]
  18. Wang, T.; Han, Q.; Chu, F.; Feng, Z. A new SKRgram based demodulation technique for planet bearing fault detection. J. Sound Vib. 2016, 385, 330–349. [Google Scholar] [CrossRef]
  19. Miao, Y.; Zhao, M.; Lin, J. Improvement of kurtosis-guided-grams via Gini index for bearing fault feature identification. Meas. Sci. Technol. 2017, 28, 125001. [Google Scholar] [CrossRef]
  20. Wan, S.; Zhang, X.; Dou, L. Shannon Entropy of Binary Wavelet Packet Subbands and Its Application in Bearing Fault Extraction. Entropy 2018, 20, 260. [Google Scholar]
  21. Bearing Data Center-Case Western Reserve University, Cleveland, OH, USA. Available online: https//csegroups.case.edu/bearingdatacenter/home (accessed on 15 July 2013).
Figure 1. The simulated signals of (a) fault impulses; (b) fault impulses and noise (SNR = −12 dB).
Figure 1. The simulated signals of (a) fault impulses; (b) fault impulses and noise (SNR = −12 dB).
Entropy 20 00388 g001
Figure 2. The analysis results of envelope spectrum, square envelope spectrum and Teager energy spectrum.
Figure 2. The analysis results of envelope spectrum, square envelope spectrum and Teager energy spectrum.
Entropy 20 00388 g002
Figure 3. The curve of the Shannon entropy with different SNR.
Figure 3. The curve of the Shannon entropy with different SNR.
Entropy 20 00388 g003
Figure 4. The paving of the (a) kurtogram and (b) TEEgram.
Figure 4. The paving of the (a) kurtogram and (b) TEEgram.
Entropy 20 00388 g004
Figure 5. The processing results of TEEgram: (a) the synthetic signal of case 1 (SNR = −12 dB); (b) the synthetic signal of case 2 (SNR = −16 dB); (c) the TEEgram of (a); and (d) the TEEgram of (b).
Figure 5. The processing results of TEEgram: (a) the synthetic signal of case 1 (SNR = −12 dB); (b) the synthetic signal of case 2 (SNR = −16 dB); (c) the TEEgram of (a); and (d) the TEEgram of (b).
Entropy 20 00388 g005
Figure 6. The process of the proposed method. (a) Healthy signal; (b) Fault signal; (c) TEEgram of the healthy signal; (d) TEEgram of the fault signal; (e) TEERgram; (f) Envelope spectrum.
Figure 6. The process of the proposed method. (a) Healthy signal; (b) Fault signal; (c) TEEgram of the healthy signal; (d) TEEgram of the fault signal; (e) TEERgram; (f) Envelope spectrum.
Entropy 20 00388 g006
Figure 7. The simulated signals of (a) the inner fault impulses; (b) the noise signal; (c) the synthetic signal (SNR = −16 dB) and (d) the spectrum of (c).
Figure 7. The simulated signals of (a) the inner fault impulses; (b) the noise signal; (c) the synthetic signal (SNR = −16 dB) and (d) the spectrum of (c).
Entropy 20 00388 g007aEntropy 20 00388 g007b
Figure 8. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 8. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g008
Figure 9. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 9. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g009
Figure 10. The processing results of proposed method: (a) the TEERgram; (b) the filtered signal and (c) the envelope spectrum of (b).
Figure 10. The processing results of proposed method: (a) the TEERgram; (b) the filtered signal and (c) the envelope spectrum of (b).
Entropy 20 00388 g010
Figure 11. The bearing test stand at Case University.
Figure 11. The bearing test stand at Case University.
Entropy 20 00388 g011
Figure 12. The experimental signal of Experiment 1.
Figure 12. The experimental signal of Experiment 1.
Entropy 20 00388 g012
Figure 13. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 13. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g013
Figure 14. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 14. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g014
Figure 15. The processing results of proposed method: (a) the TEERgram; (b) the filtered signal, and (c) the envelope spectrum of (b).
Figure 15. The processing results of proposed method: (a) the TEERgram; (b) the filtered signal, and (c) the envelope spectrum of (b).
Entropy 20 00388 g015
Figure 16. The experiment platform of QPZZ-II. 1. The drive motor; 2. The coupling; 3. The bearing housing with the normal bearing seat; 4. The loading device; and 5. The bearing housing with the fail bearing seat.
Figure 16. The experiment platform of QPZZ-II. 1. The drive motor; 2. The coupling; 3. The bearing housing with the normal bearing seat; 4. The loading device; and 5. The bearing housing with the fail bearing seat.
Entropy 20 00388 g016
Figure 17. Bearing with (a) outer ring fault; and (b) inner and outer compound fault.
Figure 17. Bearing with (a) outer ring fault; and (b) inner and outer compound fault.
Entropy 20 00388 g017
Figure 18. The experimental signal of outer ring fault in experiment 2.
Figure 18. The experimental signal of outer ring fault in experiment 2.
Entropy 20 00388 g018
Figure 19. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 19. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g019
Figure 20. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 20. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g020
Figure 21. The processing results of the proposed method: (a) the TEERgram; and (b) the filtered signal; and (c) the envelope spectrum of (b).
Figure 21. The processing results of the proposed method: (a) the TEERgram; and (b) the filtered signal; and (c) the envelope spectrum of (b).
Entropy 20 00388 g021
Figure 22. The experimental signal of inner and outer ring fault in Experiment 2.
Figure 22. The experimental signal of inner and outer ring fault in Experiment 2.
Entropy 20 00388 g022
Figure 23. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 23. The processing results of KWPT: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g023
Figure 24. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Figure 24. The processing results of kurtogram: (a) the kurtogram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g024
Figure 25. The processing results of proposed method: (a) the TEERgram; and (b) the envelope spectrum of the selected band.
Figure 25. The processing results of proposed method: (a) the TEERgram; and (b) the envelope spectrum of the selected band.
Entropy 20 00388 g025
Table 1. Parameters of the simulated signal.
Table 1. Parameters of the simulated signal.
AmplitudeRotating Frequency f r (Hz)Natural Frequency (Hz)Sampling Frequency (Hz)Fault Frequency (Hz)
2201500819225
Table 2. Parameters of the inner fault simulation.
Table 2. Parameters of the inner fault simulation.
AmplitudeRotating Frequency f r (Hz)Natural Frequency (Hz)Sampling Frequency (Hz)Fault Frequency (Hz)
22025008192130
Table 3. Structure parameters of the inner fault simulation.
Table 3. Structure parameters of the inner fault simulation.
Bearing TypePitch Diameter (mm)Ball Diameter (mm)Number of BallsContact Angle (°)
SKF6023-2RS28.56.780
Table 4. Fault feature frequency of the rolling bearing.
Table 4. Fault feature frequency of the rolling bearing.
Types of FailuresOuter RingInner RingRolling Element
Defect frequencies/Hz90145118
Table 5. Structure parameters of the inner fault simulation.
Table 5. Structure parameters of the inner fault simulation.
Bearing TypePitch Diameter (mm)Ball Diameter (mm)Number of BallsContact Angle (°)
LYC6205E397.990

Share and Cite

MDPI and ACS Style

Wan, S.; Zhang, X. Teager Energy Entropy Ratio of Wavelet Packet Transform and Its Application in Bearing Fault Diagnosis. Entropy 2018, 20, 388. https://doi.org/10.3390/e20050388

AMA Style

Wan S, Zhang X. Teager Energy Entropy Ratio of Wavelet Packet Transform and Its Application in Bearing Fault Diagnosis. Entropy. 2018; 20(5):388. https://doi.org/10.3390/e20050388

Chicago/Turabian Style

Wan, Shuting, and Xiong Zhang. 2018. "Teager Energy Entropy Ratio of Wavelet Packet Transform and Its Application in Bearing Fault Diagnosis" Entropy 20, no. 5: 388. https://doi.org/10.3390/e20050388

APA Style

Wan, S., & Zhang, X. (2018). Teager Energy Entropy Ratio of Wavelet Packet Transform and Its Application in Bearing Fault Diagnosis. Entropy, 20(5), 388. https://doi.org/10.3390/e20050388

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