Next Article in Journal
Entropy as a Measure of Attractiveness and Socioeconomic Complexity in Rio de Janeiro Metropolitan Area
Next Article in Special Issue
A Cross Entropy Based Deep Neural Network Model for Road Extraction from Satellite Images
Previous Article in Journal
Optimal Digital Implementation of Fractional-Order Models in a Microcontroller
Previous Article in Special Issue
A Multiple-Input Multiple-Output Reservoir Computing System Subject to Optoelectronic Feedbacks and Mutual Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compound Fault Diagnosis of Rolling Bearing Based on Singular Negentropy Difference Spectrum and Integrated Fast Spectral Correlation

Department of Mechanical Engineering, North China Electric Power University, Baoding 071000, China
*
Author to whom correspondence should be addressed.
Submission received: 21 February 2020 / Revised: 18 March 2020 / Accepted: 21 March 2020 / Published: 23 March 2020
(This article belongs to the Special Issue Entropy-Based Algorithms for Signal Processing)

Abstract

:
Compound fault diagnosis is challenging due to the complexity, diversity and non-stationary characteristics of mechanical complex faults. In this paper, a novel compound fault separation method based on singular negentropy difference spectrum (SNDS) and integrated fast spectral correlation (IFSC) is proposed. Firstly, the original signal was de-noised by SNDS which improved the noise reduction effect of singular difference spectrum by introducing negative entropy. Secondly, the de-noised signal was analyzed by fast spectral correlation. Finally, IFSC took the fourth-order energy as the index to determine the resonance band and separate the fault features of different single fault. The proposed method is applied to analyze the simulated compound signals and the experimental vibration signals, the results show that the proposed method has excellent performance in the separation of rolling bearing composite faults.

1. Introduction

Rotating machinery is widely used in modern industries, such as helicopters, airplanes, machining centers, track loaders, mining tracks and wind turbines [1], as shown in Figure 1. Rolling bearing is one of the most common components in rotating machines. Whether the rolling bearing runs normally is directly related to the running state of the whole rotating machinery [2,3]. Due to the bad working environment and long-term working conditions, several key parts of the rolling bearing are easy to be damaged at the same time, resulting in composite faults [4]. Composite faults are more harmful to machinery than the single faults, so it is important to diagnose early complex faults. However, all kinds of fault features are closely coupled and interfere with each other, which makes the separation of rolling bearing composite faults more challenging [5].
A variety of technical system of rotating machinery fault diagnosis has been formed, such as vibration analysis [6], acoustic emission analysis [7,8], electrical signal analysis [9] and oil analysis [10]. Among them, fault diagnosis based on vibration analysis is the most widely used. The separation of rolling bearing composite faults based on vibration analysis is a hot issue. At present, researchers have proposed some methods to extract the composite fault features of rolling bearings, including resonance demodulation, spectral kurtosis, wavelet analysis, maximum correlated kurtosis deconvolution and variational mode decomposition. For example, Wang et al. [11] solved the compound fault problem by combining the meshing resonance and spectral kurtosis. Wang et al. [12] applied the adaptive spectral kurtosis to multi-fault detection. Dhamande et al. [13] proposed a fault diagnosis method combining continuous and discrete wavelet transform. He et al. [14] proposed an adaptive redundant multiwavelet packet method to diagnose compound fault. Teng et al. [15] established a novel vibration model and used empirical wavelet transform to find multiple fault feature. Lyu et al. [16] proposed an improved maximum correlated kurtosis deconvolution method based on quantum genetic algorithm to diagnose compound fault. Miao et al. [17] developed an improved parameter-adaptive variational mode decomposition for identification of compound fault. Pan et al. [18] utilized symplectic geometry mode decomposition to decompose complex signal and got a good result.
Singular value decomposition (SVD) is a nonlinear filtering method widely used in signal denoising and fault diagnosis [19,20,21]. Conventional noise reduction methods based on SVD require the feature signal to be the main component of the signal. In this way, the prominent singular values obtained by singular value decomposition correspond to the signal space, while the smaller singular values correspond to the noise space. As long as the dimensionality reduction matrix corresponding to the more prominent singular values of the first few orders is retained, the signal after noise reduction can be obtained. However, when the background noise is very strong, the characteristic signal is completely submerged by the noise, and the singular value cannot be obtained after singular value decomposition, the number of useful component signals cannot be determined by singular value difference spectrum method. In order to overcome this problem, singular negentropy difference spectrum (SNDS) was proposed in this paper.
Cyclostationarity is a widespread physical phenomenon in rotating machinery, which is manifested by the correlation between the spectrum lines of vibration signals [22]. The vibration signal of rolling bearing is a typical cycle stationary signal and its characteristic cyclic frequency with harmonic frequency components. According to the different statistical characteristic parameters, cyclostationary stochastic processes can be divided into first, second and higher order cyclostationary processes. Spectral correlation (SC) is the main tool of second-order cyclostationary analysis. SC shows the whole structure of modulation and carrier in the signal in the form of bispectrum [23]. Fast spectral correlation (FSC) is a compromised algorithm between performance and computational efficiency which is an improved version of spectral correlation based on short time Fourier transform [24]. FSC can effectively identify the periodic impact components in the bearing, and has a high calculation efficiency. Therefore, FSC is chosen as the fault feature extraction method in this paper. However, FSC can only extract the composite fault features together, it can not separate the composite fault features. Based on FSC, integrated fast spectral correlation (IFSC) is proposed. It can identify the resonance band through the fourth order energy, and separate the composite fault. In the research of bearing composite fault diagnosis, it is found that the existence of noise affects the separation of composite fault. Therefore, in this paper, SNDS was used to de-noise the composite fault signal, and then IFSC separated the de-noised signal.
The rest of the paper is organized as follows. Section 2 describes the principles of SVD, SNDS and IFSC. Section 3 shows the specific process of the proposed method. In Section 4 and Section 5, the effectiveness of the proposed method is verified by simulation signal and experimental signal, respectively. The conclusions are presented in Section 6.

2. Basic Theory

2.1. Singular Negentropy Difference Spectrum

2.1.1. SVD

Assuming that the original discrete signal X = [x(1), x(2), …, x(N)], based on the theory of phase space reconstruction, the Hankel matrix A is constructed as follows [19]:
A = x ( 1 ) x ( 2 ) x ( n ) x ( 2 ) x ( 3 ) x ( n + 1 ) x ( N n + 1 ) x ( N n + 2 ) x ( N ) ,
where N is the length of the signal X, 1 < n < N. Let m = N-n+1, ARm×n which is the reconstructed attractor orbital matrix, and then perform singular value decomposition on it.
Perform singular value decomposition (SVD) on the matrix A, then the following equation can be obtained:
A = U S V T = i = 1 k u i λ i v i T = i = 1 k A i ,
where U = [u1, u2, …, um]∈Rm×m, V = [v1, v2, …, vn]∈Rn×n and i = 1, 2, …, k, k = min(m, n). S = diag(λ1, λ2,..., λk) is a diagonal matrix arranged in descending order, and its diagonal element is the singular value of matrix A. Ai is the submatrix corresponding to singular value λi obtained via SVD decomposition.
The submatrix Ai is inversely transformed to obtain the component signal Pi. The result of linear superposition of k component signals obtained in this way is the original discrete signal X, that is:
X = P 1 + P 2 + + P k ,

2.1.2. Singular Negentropy Difference Spectrum

Shannon Entropy is an effective index to quantitatively evaluate the uncertainties of signal or system state [25]. When the impact component appears, the Shannon entropy becomes smaller. The negentropy was defined, which is the negative value of the Shannon entropy, in order to keep the two changing regularities.
The definition of negentropy is as follows:
E y = i = 1 n ( p i log 2 p i ) ,
where y = (y1, y2, …, yn) denotes a random variable, pi is the probability of yi.
In order to effectively preserve useful fault features and minimize the impact of noise, the singular negentropy difference spectrum (SNDS) based on SVD was introduced in this paper. The definition of SNDS is as follows:
S N D S i = E y ( i ) E y ( i + 1 ) ,
where Ey(i) means the negentropy of the first i-order reconstructed signal, i denotes the number of singular values.
The schematic diagram of SNDS is shown in Figure 2. When the negentropy of two reconstructed signals are quite different, a significant peak will appear in the difference spectrum. There will inevitably be a maximum peak in the whole difference spectrum, that is, the position SNDSq of the maximum mutation of the peak degree. The maximum break point SNDSq not only shows that there are abundant fault impact characteristics in the reconstructed signal, but also shows the boundary between the useful signal and the noise component. Thus, the first q-order reconstructed signal is selected.

2.2. Integrated Fast Spectral Correlation

2.2.1. Fast Spectral Correlation

For a non-stationary random signal x(tn), its spectral correlation is given by [24]:
S C x ( α , f ) = 1 F s 2 n = τ = R x ( t n , τ ) e j 2 π α n 1 F s e j 2 π f τ 1 F s ,
where Fs denotes the sampling frequency, tn means the sampling time, tn = n /Fs, R x ( t n , τ ) denotes the cyclic autocorrelation function of x(tn), τ denotes time delay, α means the cyclic frequency, f means frequency.
Fast spectral correlation (FSC) is an improved version of spectral correlation, which shorten the calculation time and improve the efficiency by short-time Fourier transform (STFT) [24]. STFT of signal x(tn) is as follows:
X STFT ( i , f k ) = n = 0 N w 1 x [ i R + n ] w [ n ] e j 2 π n f k F s ,
The phase-corrected STFT might be presented as follows:
X w ( i , f k ) = n = 0 L 1 x [ n ] w [ n i R ] e j 2 π n f k F s = X STFT ( i , f k ) e j 2 π i R f k F s
where Nw is the window length of STFT, R denotes the block shift, w[n] means the function of time index n, fk denotes the k-th discrete frequency and fk = kΔff is the frequency resolution and it is expressed as Δf = Fs/Nw), L is the length of signal x(tn).
Assume f = fk = kΔf and α = pΔf + δ, then f-α = fk-αfk-p, hence αpΔf. Substitute these results into Equation (8)
X w ( i , f k α ) X w ( i , f k p ) e j 2 π ( α F s p Δ f F s ) ( i R + N 0 ) ,
Substituting Equation (9) into Equation (6), the following equation is described as:
S x ( α , f k ; p ) = 1 K w 2 F s i = 0 K 1 X w ( i , f k ) X w ( i , f k p ) * e j 2 π ( α F s p N w ) ( i R + N 0 ) = 1 K w 2 F s D F T i α X S T F T ( i , f k ) X S T F T ( i , f k p ) * e j 2 π N 0 ( α F s p N w )
If p = 0, the band [fkf/2, fkf/2] is flowed by the energy. Otherwise, the energy will flow between bands [fkf/2, fkf/2] and [fk-pf/2, fk-pf/2].
Therefore, the expression of fast spectral correlation is given by [24]:
S x F a s t ( α , f ) = p = 0 p S x ( α , f ; p ) p = 0 p R w ( α p Δ f ) R w ( 0 ) ,
where R w ( α ) = n = 0 N w 1 w n 2 e j 2 π ( n N 0 ) α F s represents the kernel function, and when α=0 the expression of kernel function is as R w ( 0 ) = w 2 .

2.2.2. Integrated Fast Spectral Correlation

In order to effectively separate the fault features, the integrated fast spectral correlation (IFSC) theory was proposed in this paper. The spectral correlation matrix S obtained by fast spectral correlation is an I×J matrix, whose dimension is expressed as frequencies × cyclic frequencies. After the spectral correlation matrix is obtained, the fourth-order energy at all the cyclic frequencies is added together, and the resonance band can be identified by observing the distribution of the accumulated energy along the frequency axis. In general, different resonance bands in complex faults represent different characteristics of faults. By observing the resonance band, the frequency range [f1, f2] of the resonance band is determined, and the spectral correlation matrix is integrated to obtain the integrated fast spectral correlation results. The expression of fourth-order energy is given by:
E i = j = 1 J S i , j 4 ,

3. The Proposed Method

Since the actual collected rolling bearing signal is usually covered by many noise components which have an interference effect on the extraction of fault features. In order to effectively improve signal-to-noise ratio and separate composite faults, the following steps are taken:
Step 1: Collect the composite fault signal.
Step 2: Singular negentropy difference spectrum is applied to separate the trend component which contains the fault feature information from the interference component.
Step 3: Fast spectral correlation is utilized to obtain the fast spectral correlation spectrum, the fourth-order energy of whose frequency axis is added together. As the result, the resonance band can be identified.
Step 4: Fast spectral correlation is used to separate different faults and obtain the enhanced envelope spectrum of each fault frequency band.
The flow chart of the method presented in this paper is shown in Figure 3.

4. Simulation Analysis

In order to verify the effectiveness of the proposed method, a simulated composite fault signal of inner and outer ring is performed.
x 1 ( t ) = 3 e 350 t 1 sin ( 2 π f 1 t ) , t 1 = mod ( t , 1 / f i ) x 2 ( t ) = 2.5 e 400 t 2 sin ( 2 π f 2 t ) , t 2 = m o d ( t , 1 / f o ) x 3 ( t ) = x 1 ( t ) + x 2 ( t ) x ( t ) = x 3 ( t ) + n ( t ) .
where x1(t) is the simulation signal of rolling bearing with inner ring fault whose inner ring fault characteristic frequency is expressed as fi = 130 Hz and natural frequency f1 is expressed as f1 = 3000 Hz. x2(t) is the simulation signal of rolling bearing with outer ring fault whose outer ring fault characteristic frequency is expressed as fo = 90 Hz and natural frequency f2 is expressed as f2 = 1000 Hz. n(t) is white noise. In the simulation, the sampling frequency of the example signal is fs = 8192 Hz and the sampling number is N = 4096.
Figure 4a shows the time domain waveform of the simulation signal. It can be seen from Figure 4a that the simulation signals x1(t), x2(t) and x3(t) all contain periodic impulse components, while the periodic impulse components in x(t) have been submerged by noise. Figure 4b shows that the envelope spectrum of x(t) has no prominent frequency components. The Hankel matrix is constructed from the original signal and processed by SVD. The negentropy of first 50 points and singular negentropy difference spectrum are shown in Figure 5. It can be found that the 42th point is the maximum mutation points of the difference spectrum, retaining the first 42 singular values obtained by SVD processing, and the other singular values are all set to 0. The reconstruction signal shown in Figure 6 is obtained by singular value reconstruction. Compared with Figure 6 and Figure 4a, the periodic impulse components are more obvious in the time domain waveform after SVD reconstruction.
Then, the fast spectral correlation is performed on the reconstructed signal, and the fast spectral correlation spectrum shown in Figure 7 is obtained. Two resonance bands are found in Figure 7, among which the red box area means the resonance band of the inner ring fault and the green box area is the resonance band of the outer ring fault. Next, the fourth-order energy of the ordinate in Figure 7 is carried out to obtain the resonance band screening diagram shown in Figure 8. According to Figure 8, the two resonance bands are 512–1408 Hz and 2304–3712 Hz, respectively. Fast spectral correlation is conducted for the frequency ranges of the two resonance bands respectively, the outer and inner ring faults after separation are finally obtained as shown in Figure 9 and Figure 10. Figure 9 describes the fast spectral correlation and enhanced envelope spectrum of outer ring fault, from which fo, 2fo, 3fo, 4fo, 5fo, 6fo and 7fo can be recognized obviously. Figure 10 describes the fast spectral correlation and enhanced envelope spectrum of inner ring fault. fi and its 2X, 3X, 4X, 5X frequency doubling can be clearly identified from Figure 10.
In order to show the superiority of the proposed method, the integrated fast spectral correlation based on singular difference spectrum and spectral kurtosis based on wavelet [26] are used to analyze the simulation signal. Figure 11 shows the singular difference spectrum of the simulation signal. As shown in Figure 11, i = 44 is the mutation point. The singular values after i= 44 are set to 0, and the time domain diagram of the reconstructed signal is shown in Figure 12. While the periodic impulse in the reconstructed signal is more obvious than that in the simulation signal, it still contains some noise components. The fast spectral correlation spectrum is shown in Figure 13. Compared with Figure 7, the resonance bands of outer and inner ring faults cannot be clearly identified.
Figure 14 shows the separation of resonance bands, from which, the inner fault resonance band can be accurately identified, while the outer fault resonance band cannot be clearly identified. Figure 15 shows the fast spectral correlation and enhanced envelope spectrum of outer ring fault. From Figure 15a,b, only fo, 2fo, 3fo, 6fo and 7fo can be recognized. It can be also found from Figure 15b that there are many interference components. The fast spectral correlation and enhanced envelope spectrum of inner ring fault is shown in Figure 16. fi and its 2X, 3X, 4X, 5X frequency doubling can be clearly identified from Figure 16. Compared with Figure 9 and Figure 10, the separation effect of the proposed method is better than that of the integrated fast spectral correlation based on singular difference spectrum.
As shown in Figure 17, periodic impact components cannot be noticed in the time-domain waveform after wavelet denoising. Kurtogram of the denoised signal is shown in Figure 18, and two frequency bands was selected according to it. It can be seen from Figure 19a,b, the envelope spectrums contain the characteristic frequency of inner ring and outer ring, but the composite fault is not separated into single ones. The results show that the proposed method has a better separation performance than the method combing wavelet transform and spectral kurtosis.

5. Experimental Analysis

5.1. Introduction of Experiment

In order to prove the effectiveness of the proposed method in practical application, QPZZ model test-bed is applied to simulate the composite fault of inner and outer rings of rolling bearings. Figure 20a, b and c show the experimental devices, sensor distribution and composite faults of inner and outer rings, respectively. The test bench is mainly composed of a motor, a shaft support, coupling, bearings, block and accelerometers. The geometric parameters of rolling bearings are shown in Table 1. The speed of motor is 1470 r/min during the experiment. The sampling frequency fs is 12,800 Hz, the sampling point N is 6400 and the rotating frequency fr is 24.5 Hz. The characteristic frequency of inner ring fault fi is 132 Hz and the outer ring fault characteristic frequency fo is 88 Hz.

5.2. Compound Fault Diagnosis

Figure 21a,b are the time domain waveform and envelope spectrogram of the original signal, respectively. From the time-domain waveform of the original signal, we can see that it contains many noise components, which cannot be identified from the envelope spectrum. What is more, the fault type cannot be determined from Figure 21b. Figure 22 shows the negentropy and singular negentropy difference spectrum of the signal reconstructed from the first i-th singular values after SVD decomposition. From Figure 22, we can see that i = 47 is the mutation position of negentropy value. Hence, the negentropy values after i= 47 are set to 0. The signal is reconstructed, then the time domain waveform of Figure 23 is obtained.
The fast spectral correlation analysis of the reconstructed signal is carried out, and the fast spectral correlation spectrum of Figure 24 is obtained. From Figure 24, two resonance bands are found, among which the green rectangular frame denotes the resonance band of the outer ring fault and the red rectangular frame means the resonance band of the inner ring fault. By calculating the fourth-order energy of the frequency of Figure 24, the resonance bands of the outer and inner ring faults can be separated. It can be seen from Figure 25 that the resonance band of outer ring fault is 1600 Hz–2600 Hz and the inner ring fault is 3200 Hz–6400 Hz. Fast spectral correlation analysis of the two resonance bands is carried out, and then the fast spectral correlation results shown in Figure 26 and Figure 27 are obtained.
The result of integrated fast spectral correlation method based on singular difference spectrum is compared to show the superiority of SNDS. Figure 28 shows the singular difference spectrum of the reconstructed signals which is composed of the first i-th singular values. As shown in Figure 28, i=44 is the mutation point. The singular values after i = 44 are set to 0, and the time domain diagram of reconstructed signal is shown in Figure 29. Although the periodic pulse in reconstructed signal is more obvious than that in simulated signal, it still contains some noise components. The fast spectral correlation spectrum is shown in Figure 30. Compared with Figure 24, the resonance bands of outer and inner ring faults cannot be clearly identified. Figure 31 shows the separation of the resonance band. It can be seen from Figure 31 that the internal fault resonance band can be accurately identified, while the external fault resonance band cannot be clearly identified. Figure 32 shows the fast spectral correlation and enhanced envelope spectrum of the outer ring fault. From Figure 32a,b, we can only identify fo, 2fo, 3fo and 4fo. Figure 33 describes the fast spectral correlation and enhanced envelope spectra of inner ring faults. fi and its 2X, 3X, 4X and 5X multiples can be clearly identified in Figure 33. Compared with Figure 26 and Figure 27, the separation effect of this method is better than that of the integrated fast spectral correlation method based on singular difference spectrum.
Similar to simulation analysis, the method based on wavelet transform and spectral kurtosis was performed to compare with the proposed method. It can be seen from Figure 34 that the period impact components are more obvious after wavelet denoising. The two frequency bands with the strongest impact were selected from Kurtogram shown in Figure 35. The fault with weak impact is easy to be submerged by the fault with strong impact. Hence, as shown in Figure 36, only the outer ring fault is separated. The frequency of inner ring fault cannot be identified from both Figure 36a,b. Comparing Figure 26 and Figure 27 with Figure 36, the proposed method has superior capability to separate the composite faults than the method based on wavelet transform and spectral kurtosis.

6. Conclusions

In this paper, a compound fault feature separation method based on singular negentropy difference spectrum and integrated fast spectral correlation was proposed. Singular negentropy difference spectrum was applied to de-noise and then integrated fast spectral correlation was performed. The results show that the composite fault signal can be separated effectively and accurately. The following conclusions can be drawn:
(1)
Singular negentropy difference spectrum (SNDS) can adaptively determine the effective singular value, effectively remove the noise components and retain useful fault information. What is more, the comparison between SNDS and the singular difference spectrum shows that SNDS has better denoising performance.
(2)
The fourth-order energy was used as the index by integrated fast spectrum correlation (IFSC) to select different resonance bands, so as to realize the separation of different faults. The method combining wavelet transform and spectral kurtosis was used to compare with the proposed method in this paper, the results show that the proposed method can separate the composite faults better.
(3)
Limited to the experimental conditions, the composite fault diagnosis of rolling bearing is only discussed in this paper. There are many other composite faults of rotating machines, such as compound fault of gear, compound fault of gear and rolling bearing. In the future, we will continue to study the difficult problem of fault diagnosis of other fault modes.

Author Contributions

T.T. conceived the method, performed the experiment and completed the paper writing; G.T. provided guidance. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Y.; Wang, X.; Liu, Z.; Liang, X.; Si, S. The Entropy algorithm and its variants in the fault diagnosis of rotating machinery: A Review. IEEE Access 2018, 6, 66723–66741. [Google Scholar] [CrossRef]
  2. Yan, X.; Liu, Y.; Jia, M.; Zhu, Y. A Multi-stage Hybrid fault diagnosis approach for rolling element bearing under various working conditions. IEEE Access 2019, 7, 138426–138441. [Google Scholar] [CrossRef]
  3. Wang, Z.; Wang, J.; Cai, W.; Zhou, J.; Du, W.; Wang, J.; He, G.; He, H. Application of an improved ensemble local mean decomposition method for gearbox composite fault diagnosis. Complexity 2019. [Google Scholar] [CrossRef] [Green Version]
  4. Wang, Z.J.; Zheng, L.K.; Du, W.H. A novel method for intelligent fault diagnosis of bearing based on capsule neural network. Complexity 2019. [Google Scholar] [CrossRef] [Green Version]
  5. Wang, Z.J.; Du, W.H.; Wang, J.Y. Research and application of improved adaptive MOMEDA fault diagnosis method. Measurement 2019, 140, 63–67. [Google Scholar] [CrossRef]
  6. Zhao, D.; Li, J.; Cheng, W.; Wen, W. Compound faults detection of rolling element bearing based on the generalized demodulation algorithm under time-varying rotational speed. J. Sound Vib. 2016, 378, 109–123. [Google Scholar] [CrossRef]
  7. Glowacz, A. Acoustic fault analysis of three commutator motors. Mech. Syst. Signal Process. 2019, 133, 106226. [Google Scholar] [CrossRef]
  8. Glowacz, A. Recognition of acoustic signals of loaded synchronous motor using FFT, MSAF-5 and LSVM. Arch. Acoust. 2015, 40, 197–203. [Google Scholar] [CrossRef] [Green Version]
  9. Singh, S.; Kumar, N. Detection of bearing faults in mechanical systems using stator current monitoring. IEEE Trans. Ind. Inform. 2016, 13, 1341–1349. [Google Scholar] [CrossRef]
  10. Newell, G.E. Oil analysis cost-effective machine condition monitoring technique. Ind. Lubr. Tribol. 1999, 51, 119–124. [Google Scholar] [CrossRef]
  11. Wang, T.Y.; Chu, F.L.; Han, Q.K.; Kong, Y. Compound faults detection in gearbox via meshing resonance and spectral kurtosis methods. J. Sound Vib. 2017, 392, 367–381. [Google Scholar] [CrossRef]
  12. Wang, Y.X.; Liang, m. Identification of multiple transient faults based on the adaptive spectral kurtosis method. J. Sound Vib. 2012, 331, 470–486. [Google Scholar] [CrossRef]
  13. Laxmikant, S.D.; Mangesh, B.C. Compound gear-bearing fault feature extraction using statistical features based on time-frequency method. Measurement 2018, 125, 63–77. [Google Scholar]
  14. He, S.; Chen, J.; Zhou, Z.; Zi, Y.; Wang, Y.; Wang, X. Multifractal entropy based adaptive multiwavelet construction and its application for mechanical compound-fault diagnosis. Mech. Syst. Signal Process. 2016, 76–77, 742–758. [Google Scholar] [CrossRef]
  15. Teng, W.; Ding, X.; Cheng, H.; Han, C.; Liu, Y.; Mu, H. Compound faults diagnosis and analysis for a wind turbine gearbox via a novel vibration model and empirical wavelet transform. Renew. Energy 2019, 136, 393–402. [Google Scholar] [CrossRef]
  16. Lyu, X.; Hu, Z.Q.; Wang, Q. Application of improved MCKD method based on QGA in planetary gear compound fault diagnosis. Measurement 2019, 139, 236–248. [Google Scholar] [CrossRef]
  17. Miao, Y.H.; Zhao, M.; Lin, J. Identification of mechanical compound-fault based on the improved parameter-adaptive variational mode decomposition. ISA Trans. 2019, 74, 82–95. [Google Scholar] [CrossRef]
  18. Pan, H.; Yang, Y.; Li, X.; Zheng, J.; Cheng, J. Symplectic geometry mode decomposition and its application to rotating machinery compound fault diagnosis. Mech. Syst. Signal Process. 2019, 114, 189–211. [Google Scholar] [CrossRef]
  19. Luo, B.; Hancock, E. Structural graph matching using the EM algorithm and singular value decomposition. IEEE Trans. Pattern Anal. 2001, 23, 1120–1136. [Google Scholar]
  20. Pang, B.; Tang, G.J.; Tian, T. Rolling bearing fault diagnosis based on SVDP-based kurtogram and Iterative autocorrelation of Teager energy operator. IEEE Access 2019, 7, 77222–77237. [Google Scholar] [CrossRef]
  21. Vivekananda, B.K.; Sengupta, I.; Das, A. An adaptive audio watermarking based on the singular value decomposition in the wavelet domain. Digit. Signal Process. 2010, 20, 1547–1558. [Google Scholar]
  22. Antoni, J. Cyclic spectral analysis in practice. Mech. Syst. Signal Process. 2007, 21, 597–630. [Google Scholar] [CrossRef]
  23. Gardner, W.A. Spectral correlation theory of cyclostationary time-series. Signal Process. 1986, 11, 13–36. [Google Scholar] [CrossRef]
  24. Antoni, J.; Ge, X.; Hamzaoui, N. Fast computation of the spectral correlation. Mech. Syst. Signal Process. 2017, 92, 248–277. [Google Scholar] [CrossRef]
  25. Wang, Y.; Alajaji, F.; Linder, T. Hybrid digital-analog coding with bandwidth compression for gaussian source-channel pairs. IEEE Trans. Commun. 2009, 57, 997–1012. [Google Scholar] [CrossRef] [Green Version]
  26. Liu, H.Y.; Huang, W.G.; Wang, S.B.; Zhu, Z.K. Adaptive spectral kurtosis filtering based on Morlet wavelet and its application for signal transients detection. Signal Process. 2014, 96, 118–124. [Google Scholar] [CrossRef]
Figure 1. Some application examples of rotating machinery [1].
Figure 1. Some application examples of rotating machinery [1].
Entropy 22 00367 g001
Figure 2. The schematic diagram of singular negentropy difference spectrum (SNDS).
Figure 2. The schematic diagram of singular negentropy difference spectrum (SNDS).
Entropy 22 00367 g002
Figure 3. The flow chart of the proposed method.
Figure 3. The flow chart of the proposed method.
Entropy 22 00367 g003
Figure 4. (a) Time domain waveform of the simulation signals; (b) the envelope spectrum of x(t).
Figure 4. (a) Time domain waveform of the simulation signals; (b) the envelope spectrum of x(t).
Entropy 22 00367 g004
Figure 5. Negentropy of first 50 points and the singular negentropy difference spectrum (SNDS).
Figure 5. Negentropy of first 50 points and the singular negentropy difference spectrum (SNDS).
Entropy 22 00367 g005
Figure 6. Time domain waveform of the reconstructed signal.
Figure 6. Time domain waveform of the reconstructed signal.
Entropy 22 00367 g006
Figure 7. The fast spectral correlation spectrum of the reconstructed signal.
Figure 7. The fast spectral correlation spectrum of the reconstructed signal.
Entropy 22 00367 g007
Figure 8. Screening of resonance bands
Figure 8. Screening of resonance bands
Entropy 22 00367 g008
Figure 9. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 9. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g009
Figure 10. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 10. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g010
Figure 11. Singular value of first 50 points and singular difference spectrum (SDS).
Figure 11. Singular value of first 50 points and singular difference spectrum (SDS).
Entropy 22 00367 g011
Figure 12. Time domain waveform of the reconstructed signal.
Figure 12. Time domain waveform of the reconstructed signal.
Entropy 22 00367 g012
Figure 13. The fast spectral correlation spectrum of the reconstructed signal.
Figure 13. The fast spectral correlation spectrum of the reconstructed signal.
Entropy 22 00367 g013
Figure 14. Screening of resonance bands
Figure 14. Screening of resonance bands
Entropy 22 00367 g014
Figure 15. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 15. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g015
Figure 16. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 16. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g016
Figure 17. Time-domain waveform after wavelet denoising.
Figure 17. Time-domain waveform after wavelet denoising.
Entropy 22 00367 g017
Figure 18. Kurtogram of the denoised signal.
Figure 18. Kurtogram of the denoised signal.
Entropy 22 00367 g018
Figure 19. The results of separation: (a) the envelope spectrum of band 1; (b) the envelope spectrum of band 2.
Figure 19. The results of separation: (a) the envelope spectrum of band 1; (b) the envelope spectrum of band 2.
Entropy 22 00367 g019
Figure 20. (a) Experiment stand; (b) sensors distribution and (c) rolling bearing with inner ring fault and outer ring fault.
Figure 20. (a) Experiment stand; (b) sensors distribution and (c) rolling bearing with inner ring fault and outer ring fault.
Entropy 22 00367 g020
Figure 21. (a) Time domain waveform of the original signal; (b) the envelope spectrum of the original signal.
Figure 21. (a) Time domain waveform of the original signal; (b) the envelope spectrum of the original signal.
Entropy 22 00367 g021
Figure 22. Negentropy of the first 50 points and singular negentropy difference spectrum (SNDS).
Figure 22. Negentropy of the first 50 points and singular negentropy difference spectrum (SNDS).
Entropy 22 00367 g022
Figure 23. Time domain waveform of the reconstructed signal.
Figure 23. Time domain waveform of the reconstructed signal.
Entropy 22 00367 g023
Figure 24. The fast spectral correlation spectrum of the reconstructed signal.
Figure 24. The fast spectral correlation spectrum of the reconstructed signal.
Entropy 22 00367 g024
Figure 25. Screening of resonance bands.
Figure 25. Screening of resonance bands.
Entropy 22 00367 g025
Figure 26. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 26. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g026
Figure 27. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 27. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g027
Figure 28. Singular value of the first 50 points and singular difference spectrum (SDS).
Figure 28. Singular value of the first 50 points and singular difference spectrum (SDS).
Entropy 22 00367 g028
Figure 29. Time domain waveform of the reconstructed signal.
Figure 29. Time domain waveform of the reconstructed signal.
Entropy 22 00367 g029
Figure 30. The fast spectral correlation spectrum of the reconstructed signal.
Figure 30. The fast spectral correlation spectrum of the reconstructed signal.
Entropy 22 00367 g030
Figure 31. Screening of resonance bands.
Figure 31. Screening of resonance bands.
Entropy 22 00367 g031
Figure 32. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 32. Outer ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g032
Figure 33. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Figure 33. Inner ring fault after separation: (a) the fast spectral correlation spectrum; (b) the enhanced envelope spectrum.
Entropy 22 00367 g033
Figure 34. Time-domain waveform after wavelet denoising.
Figure 34. Time-domain waveform after wavelet denoising.
Entropy 22 00367 g034
Figure 35. Kurtogram of denoised signal.
Figure 35. Kurtogram of denoised signal.
Entropy 22 00367 g035
Figure 36. The results of separation: (a) the envelope spectrum of band 1; (b) the envelope spectrum of band 2.
Figure 36. The results of separation: (a) the envelope spectrum of band 1; (b) the envelope spectrum of band 2.
Entropy 22 00367 g036
Table 1. Specific information about the experiment bearings.
Table 1. Specific information about the experiment bearings.
TypeDiameter of Balls, d (mm)Pith Diameter, D (mm)Number of Balls, zContact Angle, α (°)Damage Size of Inner Ring, (mm)Damage Size of Outer Ring, (mm)
SKF62057.538.5900.0080.059

Share and Cite

MDPI and ACS Style

Tang, G.; Tian, T. Compound Fault Diagnosis of Rolling Bearing Based on Singular Negentropy Difference Spectrum and Integrated Fast Spectral Correlation. Entropy 2020, 22, 367. https://0-doi-org.brum.beds.ac.uk/10.3390/e22030367

AMA Style

Tang G, Tian T. Compound Fault Diagnosis of Rolling Bearing Based on Singular Negentropy Difference Spectrum and Integrated Fast Spectral Correlation. Entropy. 2020; 22(3):367. https://0-doi-org.brum.beds.ac.uk/10.3390/e22030367

Chicago/Turabian Style

Tang, Guiji, and Tian Tian. 2020. "Compound Fault Diagnosis of Rolling Bearing Based on Singular Negentropy Difference Spectrum and Integrated Fast Spectral Correlation" Entropy 22, no. 3: 367. https://0-doi-org.brum.beds.ac.uk/10.3390/e22030367

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