Next Article in Journal
Betting, Selection, and Luck: A Long-Run Analysis of Repeated Betting Markets
Previous Article in Journal
Ordinal Patterns in Heartbeat Time Series: An Approach Using Multiscale Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Early Fault Diagnosis Method of Rolling Bearings on the Basis of Adaptive Frequency Window and Sparse Coding Shrinkage

Department of Mechanical Engineering, North China Electric Power University, Baoding 071000, China
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(6), 584; https://doi.org/10.3390/e21060584
Submission received: 6 May 2019 / Revised: 5 June 2019 / Accepted: 10 June 2019 / Published: 12 June 2019
(This article belongs to the Section Signal and Data Analysis)

Abstract

:
Early fault information of rolling bearings is weak and often submerged by background noise, easily leading to misdiagnosis or missed diagnosis. In order to solve this issue, the present paper puts forward a fault diagnosis method on the basis of adaptive frequency window (AFW) and sparse coding shrinkage (SCS). The proposed method is based on the idea of determining the resonance frequency band, extracting the narrowband signal, and envelope demodulating the extracted signal. Firstly, the paper introduces frequency window, which can slip on the frequency axis and extract the frequency band. Secondly, the double time domain feature entropy is proposed to evaluate the strength of periodic components in signal. The location of the optimal frequency window covering the resonance band caused by bearing fault is determined adaptively by this entropy index and the shifting/expanding frequency window. Thirdly, the signal corresponding to the optimal frequency window is reconstructed, and it is further filtered by the sparse coding shrinkage algorithm to highlight the impact feature and reduce the residue noise. Fourthly, the de-noised signal is demodulated by envelope operation, and the corresponding envelope spectrum is calculated. Finally, the bearing failure type can be judged by comparing the frequency corresponding to the spectral lines with larger amplitude in the envelope spectrum and the fault characteristic frequency. Two bearing vibration signals are applied to validate the proposed method. The analysis results illustrate that this method can extract more failure information and highlight the early failure feature. The data files of Case Western Reserve University for different operation conditions are used, and the proposed approach achieves a diagnostic success rate of 83.3%, superior to that of the AFW method, SCS method, and Fast Kurtogram method. The method presented in this paper can be used as a supplement to the early fault diagnosis method of rolling bearings.

1. Introduction

Rolling bearings, a mechanical component with compact structure and effective reduction of friction loss, are widely used in transportation, the energy and chemical industry, lifting machinery, and other fields. Rolling bearing failure will cause equipment shutdown, machine damage, and even serious personal injury. The fault monitoring and diagnosis of bearings is a useful way to decrease the loss of maintenance and ensure the safety of equipment operation [1,2].
In the initial stage of bearing failure, the defective size of the contact interface is small, the periodic impact force produced by the component interaction is weak, and the feature information is not obvious. Smith et al. [3] discussed the performance of envelope analysis, cepstrum prewhitening, and the benchmark method in extracting failure information from early bearing fault signals. In addition, vibration transmission attenuation, strong background noise, and the interference of multiple vibration sources increase the difficulty of feature extraction [4,5,6]. Aiming at this problem, Miao et al. [7] improved the variational mode decomposition method so that its parameters can be determined adaptively, Jiao et al. [8] proposed the hierarchical discriminating sparse coding method, while Hou et al. [9] put forward an integrated method on the basis of globally optimized sparse coding and approximate SVD. Zheng et al. [10] proposed an NLM-KSVD method for early fault diagnosis of rolling bearings based on shift-invariant K-singular value decomposition (K-SVD) with nonlocal means (NLM) sensitive atom enhancement. Although the above methods can extract weak faults of rolling bearings, complex prior knowledge is needed.
The periodic impact force produced by rolling bearing defects is very weak in the time domain, while the energy distribution is very wide in the frequency domain, arousing the high-frequency resonance of single or multiple mechanical structures such as bearing elements, bearing blocks, etc. The frequency band of resonance response contains abundant fault information, and the fault characteristics can be obtained by envelope demodulation analysis of its corresponding signal. Using resonance demodulation to extract impact response characteristics from the vibration signal is a fast and simple approach for rolling bearing fault diagnosis, the pivotal step of which is to precisely identify the resonance frequency band that contains plentiful fault information. The traditional resonance demodulation approach has three disadvantages. The first one is that the evaluation index cannot distinguish the random impact and periodic impact of the signal, which easily leads to identifying the erroneous resonance band. The second one is that the optimal frequency band may not include the whole resonant frequency band region. The third one is that the parameters of the bandpass filter need to be preset manually. Aiming at these problems, spectral kurtosis was put forward by Dyer [11] originally, then systematically defined by Antoni et al. [12,13,14], and successfully applied to the determination of the resonance frequency band of rolling bearings. In Reference [14], Fast Kurtogram (FK), a fast algorithm based on a finite impulse response filter and short-time Fourier transform, was used to adaptively identify the optimal resonance frequency band. However, the FK algorithm has two obvious shortcomings; one is that the frequency band is decomposed by dichotomy/trichotomy, causing the over-decomposition of the determined resonance band. The other is that the kurtosis index is susceptible to random noise and accidental shock; as a result, the determined frequency band is not the resonance frequency band aroused by the bearing defect. In order to improve the FK algorithm, Lei et al. [15], Gu et al. [16], and Wan et al. [17] proposed improved methods using wavelet packet to construct the filter that better matches the fault impact. Although such improvements can improve the identification accuracy, the resonance frequency band is easily segmented by the adjacent subset because of the pyramid frequency segmentation method. Su et al. [18], Chen et al. [19], and Wan et al. [20] utilized the genetic algorithm, particle swarm optimization algorithm, and water cycle algorithm to determine the function of the Morlet wavelet and asymmetric real Laplace wavelet, and constructed the bandpass filter with an arbitrary Q quality factor to realize adaptive matching of the resonance frequency band. However, such improvements require a lot of prior knowledge, as the optimization algorithms are relatively complex. Wang et al. [21] put forward an adaptive spectral kurtosis method using the maximum spectral kurtosis, continuously expanding the window function in the overall frequency range to determine the optimal resonance frequency band. This improvement can effectively preserve the whole resonance frequency band without the complex optimization process. In addition to the improvement of filter design, many scholars also discussed the evaluation criteria of bandpass filter signal. Yu et al. [22] and Obuchowski et al. [23] described the bearing fault response more accurately by introducing some new statistical models, such as alpha-stable, Jarque–Bera, Kolmogorov–Smirnov, etc. Considering the periodicity of fault impact, Zhang et al. [24] used correlation kurtosis instead of kurtosis to identify the resonance frequency band in diagnosis. Wang et al. [25,26] concluded that the envelope spectrum of the demodulated signal corresponding to the resonance frequency band is sparse, and proposed the Sparsogram method that calculates the kurtosis of the envelope spectrum. Antoni et al. [27] utilized the spectral negentropy of the time domain to represent the impact characteristic, and the spectral negentropy of the frequency domain to represent the cyclostationality characteristic. A method named Infogram based on average spectral negentropy was proposed.
Although most of the background noise in the signal can be filtered by extracting the resonance band caused by bearing fault through a bandpass filter, the noise in the filter passband cannot be filtered out. Especially in the analysis of strong noise signal, the effective removal of noise in the signal can highlight the relevant feature information and increase the failure identification accuracy. Sparse code shrinkage (SCS) was put forward on the basis of maximum-likelihood estimation theory [28]. This method was successfully applied not only to sparse analysis and image de-noising, but also to fault signal diagnosis of mechanical equipment. Lin et al. [29] utilized the maximum-likelihood estimation function of SCS as the soft threshold function of wavelet threshold de-noising for diagnosis pf gear faults, but the diagnosis effort was not ideal under the strong noise interference due to a lack of pre-filtering. Wang et al. [30] put forward an integrated method based on wavelet packet and SCS, which can filter the noise contained in the crack signal, but which needs manual selection of the wavelet packet coefficients to reconstruct the signal for subsequent SCS processing. Yu et al. [31] combined intrinsic time-scale decomposition (ITD) and SCS to diagnose the bearing failure type. However, the baseline of ITD is obtained by the linear transformation of the signal itself, making the waveform of the decomposed proper rotation components appear blurred and distorted, then leading to the inaccurate decomposition of a multi-component AM/FM signal.
Inspired by the above references, a fault diagnosis method on the basis of adaptive frequency window and sparse code shrinkage (i.e., AFW-SCS) is put forward to address the fact that the feature information of rolling bearing early faults is easily submerged by background noise. The AFW can slip along the whole frequency axis and extract the frequency band arbitrarily. Dual-time domain feature entropy is proposed to evaluate the strength of periodic components of the reconstructed signal corresponding to the AFW. According to this evaluation index, the optimal frequency window position is determined adaptively by shifting and expanding the window bandwidth, and covering the resonance band caused by the bearing fault. The reconstructed signal corresponding to the optimal frequency window is further de-noised by the SCS algorithm to make the impact characteristic more prominent and the diagnosis result more accurate. Two experimental signals and the data files of Case Western Reserve University for different operation conditions are applied to validate the effectiveness and the universality of AFW-SCS approach.
The paper consists of five sections. Section 2 illustrates the concept of the frequency window and its adaptive optimization process. Section 3 reviews the theoretical knowledge of SCS. Section 4 exhibits the implementation procedure of the AFW-SCS method. Section 5 validates the feasibility and superiority of the AFW-SCS method through two experimental signals and some existing analysis methods. Section 6 draws some conclusions.

2. Adaptive Frequency Window

According to the idea of Reference [32], the adaptive frequency window is proposed to overcome the shortcomings of spectral boundary division caused by frequency domain extreme points. In the next section, the concept of frequency window is firstly introduced. Then, the dual-time domain feature entropy is proposed and used as the evaluation index of the frequency window. Finally, the optimization process of the adaptive frequency window is introduced.

2.1. Concept of Frequency Window

Empirical wavelet transform (EWT) can extract the single AM/FM component with a compact spectrum by dividing the Fourier spectrum of the signal and performing orthogonal empirical wavelet transform. Assuming that the signal x(t) is composed of N single AM/FM components, the spectrum range is normalized to [0, π]. In order to extract all single components, it is necessary to divide [0, π] into N continuous intervals, and the divided spectrum is as shown in Figure 1. The N continuous intervals of the spectrum are expressed as Λn = [wn−1, wn] (n = 1, 2, ..., N). The shadow part with the center of wn and the width of 2τn is the transitional region of each segmentation interval, and satisfies the following requirement:
U n = 1 N Λ n = [ 0 , π ]
The distribution of extreme points in the signal frequency domain X(f) directly determines the process of signal spectrum segmentation. The extreme points caused by strong background noise will seriously affect the segmentation process, and then affect the EWT processing effect. For this reason, Deng et al. [33] introduced the concept of the frequency window and the signal spectrum segmentation process, which is described in Figure 2.
As seen, the frequency window Λ is expressed as [wa, wb], which represents the center frequency of the upper and lower cut-off bands of the window, where the shadow part is the transitional region of the segment with a width of 2τ.
After the signal spectrum is divided by the frequency window, the empirical wavelet function ψ ^ ( w ) is constructed according to the Meyer wavelet construction method, the equation of which is as follows:
ψ ^ ( w ) = { 1 , w a + τ | w | w b τ cos [ π 2 β ( 1 2 τ ( | w | w b + τ ) ) ] , w b τ | w | w b + τ sin [ π 2 β ( 1 2 τ ( | w | w a + τ ) ) ] , w a τ | w | w a + τ 0 , else
where the relevant parameters satisfy the following requirement:
{ β ( x ) = x 4 ( 35 84 x + 70 x 2 20 x 3 ) τ = γ w a γ < ( w b w a ) / ( w b + w a ) .
The wavelet coefficient can be calculated as
W f f ( t ) = x ( t ) , ψ ( t ) = x ( τ ) ψ ( τ t ) d τ = F 1 [ x ( w ) ψ ^ ( w ) ] ,
where represents the inner product, and F 1 ( ) represents the inverse Fourier transform.
The signal x*(t) corresponding to the frequency window can be reconstructed as follows:
x ( t ) = W f f ( t ) ψ ( w ) = F 1 [ W ^ f f ( w ) ψ ^ ( w ) ] ,
where represents the convolution operation, and W ^ f f ( w ) is the Fourier transform of W f f ( t ) .
The defect of rolling bearing parts can cause resonance, and the resonance band contains abundant fault information. If the frequency domain window covers the resonance band, the signal reconstructed based on this frequency window has obvious fault characteristics. The frequency window parameters wa and wb can slide freely, and the bandwidth of the frequency window is variable. How to determine the frequency window that covers the whole bearing fault resonance frequency band without manually doing so is a challenge.

2.2. Dual-Time Domain Feature Entropy

Dual-time domain transformation (DTDT) is an analysis method on the basis of generalized S-transform (GST) [34] and Fourier transform, representing the non-stationary feature in the dual-time domain. The GST is an extension of S-transform (ST). The ST of one-dimensional signal x(t) is illustrated as follows [35]:
S T x ( t , f ) = + x ( τ ) w ( τ t , σ ( f ) ) exp ( j 2 π f τ ) d τ ,
where τ represents the time-shift factor, f represents the carrier frequency, and w(τ − t) represents the Gauss window, the standard deviation of which is σ(f) = 1/∣f∣. The width of the window function (WF) decreases with the increase of carrier frequency. In order to enhance the adaptability of WF of ST, two adjusting factors k and p are introduced into thw standard deviation σ(f).
σ ( f ) = k / | f | p .
The window function is rewritten as
w ( τ t ) = | f | p k 2 π exp [ f 2 p ( τ t ) 2 2 k 2 ] .
The GST of one-dimensional signal x(t) can be illustrated as follows [12]:
GST x ( τ , f ) = | f | p k 2 π + x ( τ ) exp [ f 2 p ( τ t ) 2 2 k 2 2 π f τ ] d τ .
The essence of GST is to replace the fixed window function of short-time Fourier transform with the Gaussian window function, the standard deviation of which varies according to Equation (7). GST has a high time resolution for high-frequency signals and a high frequency resolution for low-frequency signals, showing a good time–frequency aggregation performance. When k = 1 and p = 1, the GST degenerates to ST.
GST can map one-dimensional time signals to the two-dimensional time–frequency plane, and Fourier inverse transform can transform the frequency domain to the time domain. Therefore, DTDT can be obtained by inverse Fourier transform of the results of GST at different times, the mathematical equation of which is described as
D T D T x ( τ 1 , τ 2 ) = + GST x ( τ 1 , f ) exp [ 2 π f τ 2 ] d f = + + x ( t ) w ( τ 1 t ) exp [ 2 π f ( τ 2 t ) ] d t d f .
Because the vibration signal x = {x(kT), k = 0, 1, 2, …, N − 1} (N represents the signal length and T represents the sample interval) collected by sensors is a discrete signal, the DTDT should be discretized. Using Fourier transform and convolution theory, Equation (4) can be rewritten as
GST x ( τ , f ) = + X ( α + f ) exp [ 2 π 2 α 2 k 2 f 2 p + i 2 π α τ ] d α ,
where X(·) is the Fourier transform of x(t), and α represents the shifting frequency. Setting f n / ( N T ) and τ k T , the discrete generalized S-transform can be expressed as
GST x ( k T , n N T ) = m = 0 N 1 X [ n + m N T ] exp [ 2 n 2 k 2 m 2 n 2 p + i 2 π m j N ] ,
where k, m, n = 0, 1, 2, …, N − 1.
Similarly, the discrete expression of DTDT can be deduced as follows:
D T D T x ( k T , i T ) = n = 0 N 1 GST x ( k T , n N T ) exp [ 2 π n k ] ,
where i, k, n = 0, 1, 2, …, N − 1.
After DTDT is performed on the signal, an N × N dimension complex matrix, named the dual-time domain transformation matrix (DTDTM), is obtained. The signal energy is mainly concentrated near the diagonal position of DTDTM; the low-frequency component of the signal will produce larger energy leakage, while the high-frequency component will not produce obvious energy leakage. In other words, at the diagonal position of DTDTM, the high-frequency component can obtain higher gain than the low-frequency component. According to this characteristic, the diagonal elements of DTDTM are extracted to construct a new one-dimensional signal, which can reflect the high-frequency characteristics with lower amplitude of the original signals.
With the application of entropy theory in various fields, various forms of entropy gradually evolved, such as information entropy, sample entropy, approximate entropy, etc. This present paper puts forward a new entropy method named dual-time domain feature entropy (DTDFE), which consists of information entropy and DTDT. According to the definition of information entropy, the equation for calculating DTDFE can be constructed as follows:
D T D F E = [ i = 1 N P i ln ( P i ) ] ln ( N ) P i = A ( i ) / i = 1 N A ( i ) i = 1 N P i = 1 } ,
where A(i) is the amplitude of the one-dimensional signal constructed by the diagonal elements of DTDTM. DTDFE can not only make full use of the information obtained by the dual-time domain analysis and better characterize the dynamic changes of the signal, but it can also suppress the low-frequency clutter interference.
Rolling bearing fault simulation signals are obtained using periodic impulse signals superimposed with Gaussian noise signals. The temporal waveforms of bearing fault simulation signals with signal-to-noise ratios (SNR) of 3 dB, 6 dB, and 10 dB are described in Figure 3. As shown in Figure 4, one, two, and three random impacts are added to the 6-dB fault simulation signal described in Figure 3b. The red and blue line waveforms represent the random impact and the 6-dB fault simulation signal, respectively. The value of DTDFE of the signals shown in Figure 3 are calculated, the variation tendency of which is shown in Figure 5a. As seen, DTDFE decreases with the increase of periodic pulse intensity in the signal. The value of DTDFE of the signals shown in Figure 4 are calculated, the variation tendency of which is shown in Figure 5b. As seen, DTDFE changes little with different random impacts. In summary, DTDFE can not only evaluate the strength of periodic components in signal, but also immunize them against random impact, with better robustness.

2.3. Optimization Process of Adaptive Frequency Window

The defect of rolling bearing parts can arouse resonance, and the resonance band contains abundant fault information. If the frequency domain window covers the resonance band, the analysis of the constructed signal corresponding to the resonance band will reduce the difficulty and improve the accuracy of diagnosis. The resonance band can be determined adaptively by shifting and expanding the frequency window, the specific steps of which are as follows:
Step 1: The fault frequency Ff of the diagnostic object is calculated, and the frequency band of the initial reference frequency window (RFW) is set as [wa, wb], wa = 0 and wb = wf = 3 × Ff. wf represents the initial bandwidth of the frequency window, the value of which is set according to the conclusion that the optimal fault bandwidth should not be less than three times the fault feature frequency [25,26].
Step 2: The shifting frequency window (SFW) is obtained by moving the RFW along the frequency axis, the frequency band of which is set as [wb, wb + wf]. The expanding frequency window (EFW) consists of RFW and SFW, the frequency band of which is set as [wa, wb + wf]. The DTDFE of the narrowband signals corresponding to RFW, EFW, and SFW is calculated, the values of which are recorded as FE1, FE2, and FE3, respectively.
Step 3: The values of FE1, FE2, and FE3 determine the shifting and expanding of the frequency domain window. The relationship of the above entropy values is judged by Equation (15).
FE 2 min ( FE     1 , FE     3 ) .
If the condition is satisfied, the narrowband signal corresponding to EFW contains more fault information and, thus, is conducive to fault diagnosis. The upper and lower cut-off frequency of EFW is used as that of the new RFW, and the new SFW and EFW are also constructed according to Step 2. The entropy values of the new RFW, EFW, and SFW are calculated, the relationships of which are judged by Equation (15). If the condition is not satisfied, the narrowband signal corresponding to RFW or SFW has the higher SNR and, thus, is conducive to fault diagnosis. The upper and lower cut-off frequency of SFW is used as that of the new RFW, and the new SFW and EFW are also constructed according to Step 2. The entropy values of the new RFW, SFW, and EFW are calculated, the relationships of which are judged by Equation (15).
Step 4: Step 3 is repeated until the upper limit of the RFW (i.e., wb + wf) is greater than fs/2. The RFW with the minimum value of DTDFE is the optimal frequency window, the extracted frequency band of which covers the whole resonance band aroused by the bearing fault. That is to say, the reconstructed signal of the optimal frequency window has the most fault information.
In order to clearly describe the influence of DTDFE criterion on RFW, SFW, and EFW, a schematic diagram is drawn in Figure 6. The frequency band of the initial RFW is [0, 100] Hz. The green, blue, and red lines correspond to RFE, EFW, and SFW, respectively.

3. Sparse Coding Shrinkage

The sparse coding shrinkage algorithm [28] is a de-noising method on the basis of sparse coding theory to deduce the hidden non-Gaussian components by using statistical characteristics of data. Hyvarinen proposed the following mathematical model to represent the probability density function of non-Gaussian signals:
p ( x ) = 1 2 d ( α + 2 ) [ α ( α + 1 ) / 2 ] ( α / 2 + 1 ) [ α ( α + 1 ) / 2 + | x / d | ] ( α + 3 ) ,
where x represents the signal with non-Gaussian statistical properties, d represents the standard deviation of the signal x, and α represents the parameter to control the sparsity of the probability density function. The greater the value of α is, the sparser the probability density function will be. For non-Gaussian signals with impact characteristics, the value of α can be set to 0.1.
The impulse component of the bearing fault vibration signal is a typical non-Gaussian component. Therefore, the algorithm is used to further process the results of AFW, making the impulse feature more obvious and the fault identification more accurate. The vibration signal of fault rolling bearing y, represented as Equation (17), consists of the defect impulse signal x and the Gauss white noise signal v.
y = x + v ,
where the statistical characteristics of x show non-Gaussian properties. The mean and variance of noise signal v are 0 and σ2, respectively.
According to the sparse probability density function model represented by Equation (13), the estimated signal x ^ of the impulse signal x can be deduced from the vibration signal y by using the maximum-likelihood estimation theory.
β = α ( α + 1 ) / 2 ,
σ = mad ( y ) ,
d = σ y 2 σ 2 ,
x ^ = sign ( y ) × max { 0 , | y | β d 2 + 1 2 ( | y | + β d ) 2 4 σ 2 ( α + 3 ) } ,
where σ represents the standard deviation of noise signal v, σ y 2 represents the standard deviation of vibration y, and mad represents the mean absolute deviation of the values in the signal.

4. Diagnosis Method Based on AFW and SCS

AFW can filter the most background noise in the signal by extracting the resonance band caused by bearing faults through a specific bandpass filter. However, the noise in the filter passband cannot be filtered out. The impulse components of bearing faults are usually non-Gaussian, while the noise components are Gaussian. SCS can effectively separate the Gaussian signal and non-Gaussian signal; thus, this algorithm is used to further process the narrowband signal obtained by AFW to highlight the impact component and then effectively eliminate the residual noise. Combining the advantages of AFW and SCS in signal processing, a fault diagnosis method named AFW-SCS is proposed in the present paper to solve the misdiagnosis or missed diagnosis caused by the difficulty of detecting the early fault feature information submerged in background noise. Comparing with the existed diagnosis methods based on resonance demodulation theory, this current work has three improvements. Firstly, the proposed AFW method can avoid the resonance band segmentation, keep the resonance band to the maximum extent, and extract more fault information. The adaptive process of AFW can be implemented by shifting and expanding the frequency window without the help of complex optimization algorithms. Secondly, the proposed DTDFE can effectively evaluate the strength of periodic impact in the signal and is not disturbed by random impact. The width and position of the resonance band can be determined more accurately by using DTDFE as the evaluation index of the frequency window. Thirdly, the SCS method is used to further eliminate non-Gaussian components of the extracted signal, so that the impact components can be highlighted, and the residue noise in the extracted signal can be reduced. The specific diagnostic steps and diagnostic flow chart of AFW-SCS are illustrated in Algorithm 1 and Figure 7, respectively.
Algorithm 1. The execution step of AFW-SCS
1: Load original vibration signal, set sampling frequency fs, and calculate fault feature frequency Ff.
2: Set the initial bandwidth of frequency window as wf = 3 × Ff and the analysis scope as [0, fs/2].
3: Construct RFW, EFW, and SFW. If RFW is [wa, wb], EFW is [wa, wb + wf] and SFW is [wb, wb + wf].
4: Run the judgment condition, wb + wf < fs/2. If condition is met, run Steps 4, 5, and 6. If condition is not satisfied, run Step 7.
5: Calculate DTDFE of RFW, EFW, and SFW, and record as FE1, FE2, and FE3.
6: Run the judgment condition, FE2 < min(FE1, FE3). If conditions are met, set the EFW as the new RFW, and return to Step 3. If condition is not satisfied, record FE1 value and RFW position (i.e., upper and lower cut-off frequency), set the SFW as the new RFW, and return to Step 3.
7: Select the frequency window corresponding to the minimum entropy and reconstruct the signal corresponding to the frequency window.
8: Perform SCS on the reconstructed signal to further reduce the noise.
9: Demodulate the de-noised signal by envelope operation, and calculate its envelope spectrum.
10: Diagnose the fault type by contrasting the frequency corresponding to the larger-amplitude spectral line in the envelope spectrum and the failure characteristic frequency.

5. Rolling Bearing Fault Signal Analysis Case

In this section, the vibration signals of the ball bearing under outer ring fault and the cylindrical bearing under inner ring fault are utilized to testify the feasibility, generality, and superiority of this diagnosis approach proposed in the present paper. The signal acquisition analyses are completed using the Case Western Reserve University (CWRU) [36] and North China Electric Power University (NCEPU) datasets.

5.1. Case 1: Ball Bearing

The vibration signals of the SKF6203 ball bearing under outer ring fault are collected from the experimental platform of CWRU. Figure 8a,b describe the photograph and structural sketch of the platform, respectively. Table 1 illustrates the structure parameters of the SKF6203 ball bearing. The acceleration sensors used to acquire the signal are fixed at the fan end and the drive end of the induction motor. In the present paper, the 2048 sample points of the fan-end accelerometer data from the 132nd data file (i.e., OR007@6_3) are selected for analysis. The corresponding defective position, defective diameter, shaft speed, sample frequency, and load are the outer ring, 0.007 inches, 1750 rpm, 12,000 Hz, and 2 horsepower (HP), respectively. The outer ring fault feature frequency fo, i.e., ballpass frequency of outer ring (BPFO) [37], is 105 Hz, calculated using Equation (22).
BPFO = Z 2 ( 1 d D cos α ) n 60 ,
where Z, d, D, α, and n represent the roller number, roller diameter, pitch diameter, contact angle, and shaft speed, respectively.
The temporal waveform and envelope spectrum of the ball bearing fault experimental signal are described in Figure 9. In the envelope spectrum of Figure 9, a spectral line with higher amplitude appears, and its frequency corresponds to 105 Hz (i.e., fo). The higher harmonic spectral lines corresponding to fo are not obvious and can easily be neglected.
The SCS method is utilized to analyze the ball bearing fault experimental signal, and the analysis result is illustrated in Figure 10. Although the signal temporal waveform becomes sparse in Figure 10, only a spectral line corresponding to fo appears in the envelope spectrum of Figure 10. The spectral lines corresponding to the higher harmonics of fo are absent. Using the AFW method to process this signal, the optimal frequency window is determined adaptively by the shifting and expanding operation, the upper and lower cut-off frequencies of which are 2835 Hz and 4725 Hz. The temporal waveform and fast fourier transform (FFT) spectrum of the reconstructed signal corresponding to the optimal frequency window are shown in Figure 11a,b. The reconstructed signal is demodulated by the envelope operation, and its envelope is obtained as described in Figure 11c. As seen, there are two spectral lines with high amplitude, the frequency of which correspond to 105 Hz (i.e., fo) and 25 Hz (i.e., fr, the rotation frequency). The higher harmonic spectral lines corresponding to fo are still not apparent. The diagnosis results of the SCS and AFW method may lead to misdiagnosis or missed diagnosis. The reconstructed signal obtained by the AFW method is further processed by the SCS method, and the diagnosis results of the AWF-SCS method are described in Figure 12. In the temporal waveform of Figure 12, the signal exhibits more sparsity than that of Figure 10a. In the envelope spectrum of Figure 12, the spectral lines corresponding to the fault frequency fo and its higher harmonic 2fo, 3fo, and 4fo can be clearly seen. In addition, the phenomenon of rotation frequency modulation is obvious, and the spectral lines corresponding to the sidebands of fo and 2fo (i.e., fo ± fr, 2fo ± fr) are also quite apparent. According to these results, the fault type can be definitively identified as an outer ring fault. For comparison, the Fast Kurtogram (FK) method [14] is used to process this ball bearing experimental signal. In the kurtogram of Figure 13, the resonance band determined by FK is marked by a red line, and the temporal waveform and envelope spectrum of its corresponding signal are illustrated in Figure 13b,c. As seen, the spectral lines corresponding to fr, fo, and 2fo can be seen.
In this case, comparing the diagnosis results of SCS, AFW, FK, and AFW-SCS, the proposed method can extract more fault information and, thus, is conducive to accurately diagnosing the fault type as an outer ring fault.

5.2. Case 2: Cylindrical Bearing

The vibration signals of the N205 cylindrical bearing under inner ring fault are collected on the QPZZ-II fault simulation platform for rotating machinery of the NCEPU. An eddy current sensor is used to acquire data. Figure 14a,b show the physical photo and the structural sketch of the QPZZ-II platform, respectively. The defective bearing is installed at the right bearing block. Table 2 illustrates the structure parameters of the N205 cylindrical bearing. To obtain the bearing vibration signal under an inner ring fault running station, a defect with a width of 0.2 mm and depth of 0.1 mm on the inner ring surface is machined using electric discharge machining technology, the physical photo of which is shown in Figure 15. Throughout the entire process of the experiment, the shaft speed, the sampling frequency, and load are set as 1440 rpm, 12,800 Hz, and 0 HP, respectively. The inner ring fault feature frequency fi, i.e., ballpass frequency of inner ring (BPFI) [37], is 172 Hz, calculated using Equation (23).
BPFI = Z 2 ( 1 + d D cos α ) n 60 ,
where the definitions of Z, d, D, α, and n are the same as those in Equation (23).
The temporal waveform and envelope spectrum of the cylindrical bearing fault experimental signal are described in Figure 16. As shown in the temporal waveform of Figure 16, six similar local waveforms appear in 0.25 s. This signal can be considered as a harmonic signal, the period of which is 0.041 s (0.25/6 ≈ 0.041). In the envelope spectrum of Figure 16, the spectral line corresponding to 24 Hz (i.e., fr, the rotation frequency) can be clearly seen, but the fault frequency cannot be detected. Fault characteristics are submerged by the rotation frequency carrier signal.
For the fact that the SCS method cannot analyze the harmonic signal, it is impossible to process this bearing signal to realize fault diagnosis. The AFW method is used to analyze the signal, and the upper and lower cut-off frequencies of the optimal frequency window are 516 Hz and 2064 Hz. The temporal waveform and FFT spectrum of the reconstructed signal corresponding to the optimal frequency window are described in Figure 17. As shown in the temporal waveform of Figure 17, the impact characteristics appear. Figure 17c is the envelope spectrum of the reconstructed signal, where the spectral lines corresponding to 24 Hz (i.e., fr), 172 Hz (i.e., fi), and 344 Hz (i.e., 2fi) can be detected, but the amplitude of those spectral lines are not prominent. The spectral line amplitude corresponding to the sideband fofr is greater than that of fo, and the spectral line corresponding to the sideband 2fo + fr is absent. These results may lead to misdiagnosis or missed diagnosis. The reconstructed signal obtained by the AFW method is further processed by the SCS method, and the diagnosis results of the AWF-SCS method are described in Figure 18. In the temporal waveform of Figure 18, the impact characteristics of the signal are more obvious and there is no noise interference. In the envelope spectrum of Figure 18, the spectral line amplitudes corresponding to 24 Hz, 172 Hz, and 344 Hz are more prominent than those of Figure 17c. In addition, the spectral lines corresponding to the sidebands of fo and 2fo (i.e., fo ± fr, 2fo ± fr) are more explicit than those of Figure 17c. According to these results, the fault type can be definitively identified as inner ring fault. For comparison, the Fast Kurtogram (FK) method [9] is used to process this cylindrical bearing experimental signal. In the kurtogram of Figure 19, the resonance band determined by FK is marked by the red line, and the temporal waveform and envelope spectrum of its corresponding signal are described in Figure 19b,c. As seen, the spectral lines corresponding to fr, fi, 2fi, and fo ± fr can be seen, but the amplitude of those spectral lines is less than that of the interference frequency.
In this case, comparing the diagnosis results of SCS, AFW, FSK, and AFW-SCS, the proposed method can extract more fault information and, thus, is conducive to accurately diagnosing the fault type as an inner ring fault.

5.3. Further Discussion

In order to evaluate the universality of the proposed method under different operation conditions, the 105th, 119th, 132nd, 170th, 186th, 199th, 209th, 222nd, 236th, 108th, 188th, and 237th data files of CWRU are selected for analysis. The data length, sample frequency, and acquisition location of these selected signals are as follows: 2048 sample points, 12,000 Hz, and drive end; the corresponding operation conditions of these signals, such as defect position (DP), defect diameter (DD), and motor load (ML), are listed in Table 3. The AFW-SCS, AFW, SCS, and FK methods are utilized to analyze these signals. If the fault characteristic frequency and its second- and third-order harmonics can be extracted clearly, the diagnosis is successful. Otherwise, the diagnosis is a failure. The diagnosis results of these methods are summarized in Table 4. The diagnostic success rates of the AFW-SCS, AFW, SCS, and FK methods are 83.3% (10/12), 41.7% (5/12), 16.7% (2/12), and 33.3% (4/12), respectively, indicating that the AFW-SCS method is effective in extracting more failure information, achieving early fault diagnosis of the rolling bearing, and having higher reliability.

6. Conclusions

This present paper put forward an adaptive approach, named AFW-SCS, for the diagnosis of rolling bearing weak faults. The proposed method is based on the idea of determining the resonance frequency band, extracting the narrowband signal, and envelope demodulating the extracted signal. The diagnosis results of the bearing vibration signals of CRWU and NCEPU illustrate that the proposed AFW-SCS method can extract more failure information, highlight the early failure feature, and accurately diagnose the bearing failure type. Compared with the published approach, the AFW-SCS method has three highlights. Firstly, the AFW method can avoid the resonance band segmentation, keep the resonance band to the maximum extent, and extract more fault information. The adaptive process of AFW can be implemented by shifting and expanding the frequency window without the help of complex optimization algorithms. Secondly, DTDFE can effectively evaluate the strength of periodic impact in the signal and is not disturbed by random impact. The width and position of the resonance band can be determined more accurately using DTDFE as the evaluation index of the frequency window. Thirdly, the SCS method can further eliminate non-Gaussian components of the signal extracted by the AFW method, so that the impact characteristics in the extracted signal can be highlighted, the residue noise in the extracted signal can be reduced, and the fault diagnosis accuracy can be improved.
Although the AFW-SCS method based on resonance demulation has advantages in the determination of resonance frequency band and the enhancement of fault feature information, its diagnostic success rate under strong noise interference is not ideal. Diagnostic methods based on dictionary learning, such as the NLM-KSVD method, perform better in anti-strong noise interference, but the process of acquiring dictionaries is complex and difficult to understand. Comparatively speaking, the AFW-SCS method is easier to understand and requires less prior knowledge. In any case, the AFW-SCS method and the NLM-KSVD method are complements to their respective algorithmic domains. In future research, we will systematically analyze the application scope of the AFW-SCS method, thoroughly study how to simplify the process of acquiring dictionaries in the NLM-KSVD method, try to combine the advantages of the two methods, and propose a bearing fault diagnosis method with an easy concept and a good diagnostic success rate.

Author Contributions

S.W. provided the experiment platform and sensors, checked the first draft, and revised some content; B.P. put forward the method, processed the experimental signals, and wrote the first draft.

Funding

This research was supported the National Natural Science Foundation of China (No. 51777075), the Fundamental Research Funds for the Central Universities (No. 2019QN131) and the Natural Science Foundation of Hebei Province, China (No. E2019502064).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, X.L.; Jia, M.P. Fault diagnosis of rolling bearing based on feature reduction with global-local margin Fisher analysis. Neurocomputing 2018, 315, 447–464. [Google Scholar] [CrossRef]
  2. Nibaldo, R.; Pablo, A.; Lida, B.; Guillermo, C. Combining multi-scale wavelet entropy and kernelized classification for bearing multi-fault diagnosis. Entropy 2019, 21, 152. [Google Scholar]
  3. Smith, W.A.; Randall, R.B. Rolling element bearing diagnostics using the Case Western Reserve University data: A benchmark study. Mech. Syst. Signal Process. 2015, 64–65, 100–131. [Google Scholar] [CrossRef]
  4. Wang, Z.J.; Zhou, J.; Wang, J.Y.; Du, W.H.; Wang, J.T.; Han, X.F.; He, G.F. A Novel Fault Diagnosis Method of Gearbox Based on Maximum Kurtosis Spectral Entropy Deconvolution. IEEE Access 2019, 7, 29520–29532. [Google Scholar] [CrossRef]
  5. Yan, X.A.; Jia, M.P. Application of CSA-VMD and optimal scale morphological slice bispectrum in enhancing outer race fault detection of rolling element bearings. Mech. Syst. Signal Process. 2019, 122, 56–86. [Google Scholar] [CrossRef]
  6. Zhao, X.J.; Qin, Y.; He, C.B.; Jia, L.M.; Kou, L.L. Rolling element bearing fault diagnosis under impulsive noise environment based on cyclic correntropy spectrum. Entropy 2019, 21, 50. [Google Scholar] [CrossRef]
  7. Miao, Y.H.; Zhao, M.; Lin, J. Identification of mechanical compound-fault based on the improved parameter-adaptive variational mode decomposition. ISA Trans. 2019, 84, 82–95. [Google Scholar] [CrossRef] [PubMed]
  8. Jiao, J.Y.; Zhao, M.; Lin, J.; Liang, K.X. Hierarchical discriminating sparse coding for weak fault feature extraction of rolling bearings. Reliab. Eng. Syst. Saf. 2019, 184, 41–54. [Google Scholar] [CrossRef]
  9. Hou, F.T.; Chen, J.; Dong, G.M. Weak fault feature extraction of rolling bearings based on globally optimized sparse coding and approximate SVD. Mech. Syst. Signal Process. 2018, 111, 234–250. [Google Scholar] [CrossRef]
  10. Zheng, K.; Li, T.L.; Su, Z.Q.; Wen, J.F.; Zhang, B.; Zhang, Y. Faults diagnosis of rolling bearings based on shift invariant K-singular value decomposition with sensitive atom nonlocal means enhancement. Measurement 2019, 135, 836–851. [Google Scholar] [CrossRef]
  11. Dyer, D.; Stewart, R.M. Detection of rolling element bearing damage by statistical vibration analysis. J. Mech. Des. 1978, 100, 229–235. [Google Scholar] [CrossRef]
  12. Antoni, J. The spectral kurtosis: A useful tool for characterising non-stationary signals. Mech. Syst. Signal Process. 2006, 20, 282–307. [Google Scholar] [CrossRef]
  13. 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]
  14. Antoni, J. Fast computation of the kurtogram for the detection of transient faults. Mech. Syst. Signal Process. 2007, 21, 108–124. [Google Scholar] [CrossRef]
  15. 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]
  16. Gu, X.H.; Yang, S.P.; Liu, Y.Q.; Hao, R.J. Rolling element bearing faults diagnosis based on kurtogram and frequency domain correlated kurtosis. Meas. Sci. Technol. 2016, 27, 125019. [Google Scholar] [CrossRef]
  17. Wan, S.T.; Peng, B. The FERgram: A rolling bearing compound fault diagnosis based on maximal overlap discrete wavelet packet transform and fault energy ratio. J. Mech. Sci. Technol. 2016, 33, 157–172. [Google Scholar] [CrossRef]
  18. Su, W.S.; Wang, F.T.; Zhu, H.; Zhang, Z.X.; Guo, Z.G. Rolling element bearing faults diagnosis based on optimal Morlet wavelet filter and autocorrelation enhancement. Mech. Syst. Signal Process. 2001, 24, 1458–1472. [Google Scholar] [CrossRef]
  19. Chen, X.L.; Zhang, B.Z.; Feng, F.Z.; Jiang, P.C. Optimal resonant band demodulation based on an improved correlated kurtosis and its application in bearing fault diagnosis. Sensors 2017, 17, 360. [Google Scholar] [CrossRef]
  20. Wan, S.T.; Peng, B. Adaptive asymmetric real Laplace wavelet filtering and its application on rolling bearing early fault diagnosis. Shock Vibr. 2019, 2019, 7475868. [Google Scholar] [CrossRef]
  21. Wang, Y.X.; Liang, M. An adaptive SK technique and its application for fault detection of rolling element bearings. Mech. Syst. Signal Process. 2011, 25, 1750–1764. [Google Scholar] [CrossRef]
  22. Yu, G.; Li, C.N.; Zhang, J.F. A new statistical modeling and detection method for rolling element bearing faults based on alpha-stable distribution. Mech. Syst. Signal Process. 2013, 41, 155–175. [Google Scholar] [CrossRef]
  23. Obuchowski, J.; Wylomanska, A.; Zimroz, R. Selection of informative frequency band in local damage detection in rotating machinery. Mech. Syst. Signal Process. 2014, 48, 138–152. [Google Scholar] [CrossRef]
  24. Zhang, X.H.; Kang, J.S.; Lei, X.; Zhan, J.M.; Teng, H.Z. A new improved Kurtogram and its application to bearing fault diagnosis. Shock Vibr. 2015, 2015, 385412. [Google Scholar] [CrossRef]
  25. Barszcz, T.; JabŁo’nski, A. A novel method for the optimal band selection for vibration signal demodulation and comparison with the Kurtogram. Mech. Syst. Signal Process. 2011, 25, 431–451. [Google Scholar] [CrossRef]
  26. 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]
  27. Antoni, J. The infogram: Entropic evidence of the signature of repetitive transients. Mech. Syst. Signal Process. 2016, 74, 73–94. [Google Scholar] [CrossRef]
  28. Hyvarinen, A. Sparse code shrinkage: De-noising of nonGaussian data by maximum likelihood estimation. Neural Comput. 1999, 11, 1739–1768. [Google Scholar] [CrossRef]
  29. Lin, J.; Zuo, M.J.; Fyfe, K.R. Mechanical fault detection based on the wavelet de-noising technique. J. Vib. Acoust. 2004, 126, 9–16. [Google Scholar] [CrossRef]
  30. Wang, X.H.; Luo, Z.M.; Hu, H.W.; Mao, H.L. Extraction of weak crack signals based on sparse code shrinkage combined with wavelet packet filtering. Appl. Acoust. 2016, 112, 53–58. [Google Scholar] [CrossRef]
  31. Yu, J.B.; Liu, H.Q. Sparse coding shrinkage in Intrinsic time-scale decomposition for weak fault feature extraction of bearings. IEEE Trans. Instrum. Meas. 2018, 67, 1579–1592. [Google Scholar] [CrossRef]
  32. Gilles, J. Empirical wavelet transform. IEEE Trans. Signal Process. 2013, 61, 3999–4010. [Google Scholar] [CrossRef]
  33. Deng, F.Y.; Qiang, Y.W.; Yang, S.P.; Han, R.J.; Liu, Y.Q. A fault diagnosis method of rolling element bearings with adaptive frequency window empirical wavelet transform. J. Xi’an Jiaotong Univ. 2018, 52, 22–29. [Google Scholar]
  34. Djurović, I.; Sejdić, E.; Jiang, J. Frequency-based window width optimization for S-transform. AEU Int. J. Electron. Commun. 2008, 62, 245–250. [Google Scholar] [CrossRef]
  35. Pinnegar, C.R. Time frequency and time-time filtering with the S-transform and TT-transform. Digit. Signal Process. 2005, 5, 604–620. [Google Scholar] [CrossRef]
  36. Case Western Reserve University Center Website. Available online: http:// csegroups.case.edu/bearingdatacenter/home (accessed on 8 July 2017).
  37. Randall, R.B.; Antoni, J. Rolling element bearing diagnostics—A tutorial. Mech. Syst. Signal Process. 2011, 25, 485–520. [Google Scholar] [CrossRef]
Figure 1. The partitioning of the frequency domain.
Figure 1. The partitioning of the frequency domain.
Entropy 21 00584 g001
Figure 2. The schematic diagram of the frequency window.
Figure 2. The schematic diagram of the frequency window.
Entropy 21 00584 g002
Figure 3. Temporal waveform of simulation fault signal under different signal-to-noise ratios (SNRs): (a) 3 dB; (b) 6 dB; (c) 10 dB.
Figure 3. Temporal waveform of simulation fault signal under different signal-to-noise ratios (SNRs): (a) 3 dB; (b) 6 dB; (c) 10 dB.
Entropy 21 00584 g003
Figure 4. Temporal waveform of 6-dB simulation fault signal under different random impacts: (a) one random impact (t); (b) two random impacts; (c) three random impacts.
Figure 4. Temporal waveform of 6-dB simulation fault signal under different random impacts: (a) one random impact (t); (b) two random impacts; (c) three random impacts.
Entropy 21 00584 g004
Figure 5. Variation tendency of dual-time domain feature entropy (DTDFE) with (a) SNR, and (b) random impact.
Figure 5. Variation tendency of dual-time domain feature entropy (DTDFE) with (a) SNR, and (b) random impact.
Entropy 21 00584 g005
Figure 6. The effect of DTDFE criteria on the frequency window (FW).
Figure 6. The effect of DTDFE criteria on the frequency window (FW).
Entropy 21 00584 g006
Figure 7. The diagnostic flow chart of the adaptive frequency window with sparse coding shrinkage (AFW-SCS).
Figure 7. The diagnostic flow chart of the adaptive frequency window with sparse coding shrinkage (AFW-SCS).
Entropy 21 00584 g007
Figure 8. Case Western Reserve University (CWRU) platform: (a) photograph; (b) structural sketch.
Figure 8. Case Western Reserve University (CWRU) platform: (a) photograph; (b) structural sketch.
Entropy 21 00584 g008
Figure 9. The ball bearing fault experimental signal: (a) temporal waveform; (b) envelope spectrum.
Figure 9. The ball bearing fault experimental signal: (a) temporal waveform; (b) envelope spectrum.
Entropy 21 00584 g009
Figure 10. The diagnosis results of SCS: (a) temporal waveform; (b) envelope spectrum.
Figure 10. The diagnosis results of SCS: (a) temporal waveform; (b) envelope spectrum.
Entropy 21 00584 g010
Figure 11. The diagnosis results of AFW: (a) extracted frequency band; (b) temporal waveform; (c) envelope spectrum.
Figure 11. The diagnosis results of AFW: (a) extracted frequency band; (b) temporal waveform; (c) envelope spectrum.
Entropy 21 00584 g011
Figure 12. The diagnosis results of AFW-SCS: (a) temporal waveform; (b) envelope spectrum.
Figure 12. The diagnosis results of AFW-SCS: (a) temporal waveform; (b) envelope spectrum.
Entropy 21 00584 g012
Figure 13. The diagnosis results of the Fast Kurtogram (FK) method: (a) kurtogram; (a) temporal waveform; (b) envelope spectrum.
Figure 13. The diagnosis results of the Fast Kurtogram (FK) method: (a) kurtogram; (a) temporal waveform; (b) envelope spectrum.
Entropy 21 00584 g013
Figure 14. North China Electric Power University (NCEPU) platform: (a) photograph; (b) sensor position; (c) structural sketch.
Figure 14. North China Electric Power University (NCEPU) platform: (a) photograph; (b) sensor position; (c) structural sketch.
Entropy 21 00584 g014
Figure 15. Physical photo of the N205 cylindrical bearing with outer ring defect.
Figure 15. Physical photo of the N205 cylindrical bearing with outer ring defect.
Entropy 21 00584 g015
Figure 16. The cylindrical bearing fault experimental signal: (a) temporal waveform; (b) envelope spectrum.
Figure 16. The cylindrical bearing fault experimental signal: (a) temporal waveform; (b) envelope spectrum.
Entropy 21 00584 g016
Figure 17. The diagnosis results of AFW: (a) extracted frequency band; (b) temporal waveform; (c) envelope spectrum.
Figure 17. The diagnosis results of AFW: (a) extracted frequency band; (b) temporal waveform; (c) envelope spectrum.
Entropy 21 00584 g017
Figure 18. The diagnosis results of AFW-SCS: (a) temporal waveform; (b) envelope spectrum.
Figure 18. The diagnosis results of AFW-SCS: (a) temporal waveform; (b) envelope spectrum.
Entropy 21 00584 g018
Figure 19. The diagnosis results of FK: (a) kurtogram; (b) temporal waveform; (c) envelope spectrum.
Figure 19. The diagnosis results of FK: (a) kurtogram; (b) temporal waveform; (c) envelope spectrum.
Entropy 21 00584 g019
Table 1. SKF6205 ball bearing structure parameters.
Table 1. SKF6205 ball bearing structure parameters.
Roller DiameterPitch DiameterRoller NumberContact Angle
7.9 mm39 mm9
Table 2. N205 cylindrical bearing structure parameters.
Table 2. N205 cylindrical bearing structure parameters.
Ball DiameterPitch DiameterBall NumberContact Angle
7.5 mm38.5 mm12
Table 3. Operation conditions of different data files. DP—defect position; DD—defect diameter; ML—motor load; HP—horsepower.
Table 3. Operation conditions of different data files. DP—defect position; DD—defect diameter; ML—motor load; HP—horsepower.
Data FileDPDD (inch)ML (HP)Data FileDPDD (inch)ML (HP)
No. 105IR0.0070No. 199OR0.0142
No. 119RE0.0071No. 188RE0.0143
No. 132OR0.0072No. 209OR0.0210
No. 108IR0.0073No. 223IR0.0211
No. 169IR0.0140No. 236RE0.0212
No. 186RE0.0141No. 237OR0.0213
IR represents inner ring; RE represents rolling element; OR represents outer ring.
Table 4. Diagnosis results of different methods. AFW—adaptive frequency window; SCS—sparse coding shrinkage; FK—Fast Kurtogram.
Table 4. Diagnosis results of different methods. AFW—adaptive frequency window; SCS—sparse coding shrinkage; FK—Fast Kurtogram.
Data FileAFW-SCSAFWSCSFKData FileAFW-SCSAFWSCSFK
No. 105×××No. 199××××
No. 119××××No. 188×××
No. 132×××No. 209
No. 108×No. 223×××
No. 169×No. 236××
No. 186×××No. 237
√ represents success, × represents failure.

Share and Cite

MDPI and ACS Style

Wan, S.; Peng, B. An Early Fault Diagnosis Method of Rolling Bearings on the Basis of Adaptive Frequency Window and Sparse Coding Shrinkage. Entropy 2019, 21, 584. https://doi.org/10.3390/e21060584

AMA Style

Wan S, Peng B. An Early Fault Diagnosis Method of Rolling Bearings on the Basis of Adaptive Frequency Window and Sparse Coding Shrinkage. Entropy. 2019; 21(6):584. https://doi.org/10.3390/e21060584

Chicago/Turabian Style

Wan, Shuting, and Bo Peng. 2019. "An Early Fault Diagnosis Method of Rolling Bearings on the Basis of Adaptive Frequency Window and Sparse Coding Shrinkage" Entropy 21, no. 6: 584. https://doi.org/10.3390/e21060584

APA Style

Wan, S., & Peng, B. (2019). An Early Fault Diagnosis Method of Rolling Bearings on the Basis of Adaptive Frequency Window and Sparse Coding Shrinkage. Entropy, 21(6), 584. https://doi.org/10.3390/e21060584

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