# A Spectrum Detection Approach for Bearing Fault Signal Based on Spectral Kurtosis

Authors: Yunfeng Li; Liqin Wang; Jian Guan
Journal: Shock and Vibration
 (2017)

---

## Abstract

According to the similarity between Morlet wavelet and fault signal and the sensitive characteristics of spectral kurtosis for the impact signal, a new wavelet spectrum detection approach based on spectral kurtosis for bearing fault signal is proposed. This method decreased the band-pass filter range and reduced the wavelet window width significantly. As a consequence, the bearing fault signal was detected adaptively, and time-frequency characteristics of the fault signal can be extracted accurately. The validity of this method was verified by the identifications of simulated shock signal and test bearing fault signal. The method provides a new understanding of wavelet spectrum detection based on spectral kurtosis for rolling element bearing fault signal.


---

## Body

## 1. Introduction

The rolling element bearing is one of the most important parts in aircraft engines. Many bearing incidents in aircraft engines have indicated that the initial damage of bearing occurred on the surface and subsurface. This initial damage appears as the peeling and develops into spalls eventually [1]. During the performance and reliability testing, the vibration signal can reveal the failure state in time and frequency domain promptly. In order to obtain the fault characteristic parameters of the rolling bearing, a variety of signal processing methods are applied to process the rolling bearing vibration signal. Since wavelet transform has the property of time-frequency localization, it is utilized to detect the signal transients from the raw signal extensively [2–6]. The wavelet transform provides an analysis method both in time domain and in frequency domain, which becomes an important method in the nonstationary signal analysis. It is well known that the similarity exists between Morlet wavelet and the impact response of the attenuation components, and the Morlet wavelet is Gaussian in the frequency domain. In addition to the above factors, the time-frequency structure of the Morlet wavelet is optimal match with the typical transients. Therefore, the Morlet wavelet is more suitable for the extraction of rotating machinery fault signal and the rolling bearing fault diagnosis.The purpose of wavelet transform is to select the proper frequency band of the correct faulty information and the best filter. Therefore, it is necessary and urgent to select the reasonable resonance frequency band for the demodulation technology. At the same time, an index unaffected by experimental conditions is needed to measure the transform result. Kurtosis is a normalized time domain statistics parameter, which is sensitive to the instantaneous signal. Thus, the spectral kurtosis (SK) is an effective tool to get the fault signals of rolling bearing and the reasonable bandwidth [7, 8]. SK is used to analyze the Gaussian components of the signal and locate these components in the frequency domain. It can be seen as a measurement of the energy distribution at each frequency. The peeling of the faulty bearing produces a non-gauss distribution signal, which can be described by the time averaging statistic. The final purpose of SK is to detect the fault signal characteristics in the filtering results [9–11].The concept of SK was proposed by Dwyer firstly [12]. Then, it was successfully applied to the rolling bearing resonance peak and resonance band selection by Antoni and Randall [13]. Fast Kurtogram based on band-pass filter and short-time Fourier transform (STFT) has made SK a practical tool for the fault diagnosis [14]. Nevertheless, in order to obtain a maximum value of the spectral kurtosis, the STFT window must be shorter than the spacing between the pulses and longer than the individual pulses. The maximum spectral kurtosis value corresponding to the band is not always the optimal band by the above method. In order to solve these problems, various kinds of improved algorithms based on spectral kurtosis and other signal processing technologies were presented, such as searching all the possible resonance frequency bands adaptively by shifting and expanding predetermined Morlet wavelet [15].In this paper, the Morlet wavelet and SK are combined. The spline curve function is utilized to obtain the optimal wavelet filter range, and different wavelet clusters are designed to detect the impact frequency. Firstly, a group of Morlet wavelet cluster with different quality factorQ is designed and the spectral kurtosis of the filtered signal is calculated, respectively. Then, the trend of spectral kurtosis is drawn by the spline curve based on a limited number of points. The signal detection range and filter bandwidth are gradually reduced. As a consequence, the accuracy of fault diagnosis is improved by the optimized filter.


## 2. Adaptive Spectral Kurtosis Filtering Based on Morlet Wavelet

### 2.1. Morlet Wavelet

The complex Morlet wavelet is the product of a complex exponential function multiplied by a Gaussian function. The expression of Morlet wavelet is(1)Morlt=σπe-σ2t2ei2πf0t.Meanwhile, the shape of Gaussian window in the frequency domain can be derived as [16](2)Morlf=Morl∗f=e-π2/σ2f-f02,where Morl(f) is the fast Fourier transformation of Morl(t). Since Morl(f) is real, Morl(f)=Morl∗(f), and the superscript ∗ denotes the complex conjugate. f0 is the center frequency of the wave and σ is waveform parameter that determines its width. For Morlet wavelet, the center frequency f0, and shape factor σ affect the shape and location simultaneously. By (1), the corresponding wave shape of Morlet wavelet is shown in Figure 1.Figure 1
Structure of Morlet wavelet.In Figure1, Morl1 is the wave of the Gaussian function, and Morl2 is the complex exponential wave function. They are expressed as(3)Morl1=e-t2σ2,Morl2=ei2πf0t.According to the Euler formula,Morl2 can be expressed as Morl2=cos⁡(2πf0t)+isin⁡(2πf0t), where Re(Morl2)=cos⁡(2πf0t). If t1 is half a cycle of Morl1 and t2 is a quarter period of Re(Morl2), t1 and t2 can be deduced as(4)t1=ln⁡m-1σ,t2=arccos⁡m2πf0,where m is the value of Morlet wave endpoint and it is infinitesimal (as shown in Figure 1). It is necessary to ensure that the Morlet wavelet is fluctuant enough and compactly supported. For Morl1 and Morl2, the period ratio T1/T2 is(5)T1T2=2t14t2=Qπln⁡m-1arccos⁡m=k,where Q is the quality factor, Q=f0/σ. When k1≤k≤k2, Morlet wavelet at most includes k2Morl2 and at least k1Morl2. If (5) is introduced into k1≤k≤k2, the expression of the quality coefficient Q can be obtained as follows:(6)k1arccos⁡mπln⁡m-1≤Q≤k2arccos⁡mπln⁡m-1.The quality coefficientQ is proportional to the fluctuation of the Morlet wavelet. The larger the quality coefficient Q is, the smaller the volatility of the Morlet wavelet is and vice versa. Therefore, when 0.5≤k≤2 and m=0.00001, the value range of Q is 0.0217≤Q≤0.0869.The corresponding family of wavelets consists of a series of son wavelets, which are generated by dilation and translation from the mother wavelet shown as follows:(7)Morla,bt=1aMorlt-ba,where a is the scale factor and b is wavelet displacement factor (the constant 1/a is used for energy normalization). The wavelet transform coefficient is defined by the following equation:(8)WTa,b=∫-∞+∞xt·Morla,b∗t·dt=xt,Morla,bt,where Morla,b∗(t) is the conjugate complex numbers of Morla,b(t) and · represents solving convolution.


### 2.2. Spline Curve Interpolation

Since the SK values are not continuous, the spline curve is used to obtain the variation trend of SK on the basis of discrete points. It can provide the basis for optimizing the analysis range.The function of spline curve interpolation is an approximation of piecewise three polynomials, but it is continuous. Its first- and second-order derivatives are both continuous. The value ofxj is yj=f(xj)(j=0,1,…,n), spline curve interpolation function is S(x), and its first-order derivative S′(xj)=mj. S(x) can be expressed as(9)Sx=∑j=0nyjαjx+mjβjx,where αj(x) and βj(x) are the interpolation function, expressed by(10)αjx=x-xj-1xj-xj-121+2x-xjxj-1-xjxj-1≤x≤xjj≠0x-xj+1xj-xj+121+2x-xjxj+1-xjxj≤x≤xj+1j≠n0Others,βjx=x-xj-1xj-xj-12x-xjxj-1≤x≤xjj≠0x-xj+1xj-xj+12x-xjxj≤x≤xj+1j≠n0Others.Under the condition of corresponding boundary, the functionS(x) can be calculated when the equations of mj(j=0,1,…,n) were obtained [17].


### 2.3. Transient Signal Detection Based on SK

The nonstationary signal is normally expressed by Wold-Cramér decomposition; that is,(11)Yt=∫-∞+∞ej2πftHt,fdXf,where the time-varying transfer function H(t,f) can be interpreted as the complex envelope or complex demodulate of process Y(t) at frequency f.The signal characteristics information inH(t,f) can be described by the spectral moments, the expression of which is(12)S2nYf≜ES2nYt,f=EHt,fdXf2ndf=EHt,f2n·S2nx.Under the conditions of ergodicity and stationarity ofH(t,f), it is easy to prove that(13)S2nYf=S2nYt,ft.The significant property of nonstationary processes is non-Gaussian. For a non-Gaussian process, spectrum accumulation that is greater than or equal to the fourth order is a nonzero value. Therefore, the fourth-order spectral accumulation is adopted, which is defined as(14)C4Yf=S4Yf-2S2Y2f.The greater the degree of signal deviation from Gaussianity is, the greater its fourth-order spectral accumulation is. Therefore, the energy normalized fourth-order spectral accumulation can be used to measure the peak of the signal process probability density at the frequencyf, that is, spectral kurtosis.The wavelet coefficients can be obtained by calculating the signal convolution by each row of the filter. The higher value of the wavelet coefficient is, the more similar the wavelet function with the signal waveform will be. The values of the wavelet coefficients are equal to the square of energy. SK is used to measure the energy spectrum. Therefore, SK achieves the maximum value when the waveforms are most similar [16]. Equation (15) gives the value of spectral kurtosis for each filter [18]:(15)kurtosisy≜C4YfS2Y2f=S4YfS2Y2f-2=meanH4t,fmeanH2t,f2-2.


### 2.4. Construction of Wavelet Clusters

Adaptive wavelet method uses spectral kurtosis as the evaluation index of Morlet wavelet filtered results. This method can optimize the latter filter parameters according to the calculation results in the previous step. Each step can improve the calculation accuracy and reduce the searching scope effectively. Compared with other bearing fault diagnosis method based on kurtosis, the accuracy of fault frequency identification and the detection bandwidth are both optimized.During the construction of octave filter with complex translation Morlet wavelet, the number of filters that cover the entire analysis band isM=K×N, where K is the number of the filters in each octave and N is the number of the octaves [19].The central frequency of the filter is expressed as(16)fmk=fmax-fmin21/Knn=1,2,…,M,where [fmin,fmax] is the frequency band range of the bearing characteristic frequencies in the operating mode. Bandwidth of each filter is expressed by(17)σnK=fnK·2K.The value ofK and N and the frequency analysis range [fmin,fmax] containing all possible fault frequencies are initially set. The center frequencies and bandwidths can be obtained by (16) and (17). The results are introduced into (1); then, N Morlet wavelet clusters for signal processing are obtained.When the frequency analysis range is optimized by the spline curve, the same number of filters is allocated within a smaller frequency range. As a result, the fault center frequency is more accurate and the bandwidth of the filter is smaller.


### 2.5. Adaptive Wavelet Method Based on Spectrum Diagnosis

Obviously, the proportion of the impact component in each filter window is different and it directly influences the spectral kurtosis of the filtering results. The former one is large means that the latter one is also large. The value of the quality factorQ affects the shape of the wavelet clusters. Therefore, each quality factor corresponds to an optimization process. For a rolling element bearing, the signal components related to the rotational speed usually cannot be extracted from the time-frequency representation directly. A new proposed method is presented to solve the above problems, as shown in Figure 2:(1)
In the frequency band range containing the bearing characteristic frequencies[fmin1,fmax1], it can calculate the spectral kurtosis of the filtered result and describe the spectral kurtosis by the spline curve interpolation along the frequency axis.(2)
By solving the derivation of the spline curve, the peak containing maximum spectral kurtosis and the valley on both sides of the maximum spectral kurtosis could be obtained, respectively, which is the frequency range of the next optimizing calculation[fmin2,fmax2].(3)
The quality factorQ is changed and the range of Morlet wavelet is reduced to [fmin2,fmax2]; the spectral kurtosis of the filter bandwidths can be calculated continually. The spline curve of the spectral kurtosis can be obtained by the spline curve interpolation. However, the bandwidth of the filter is smaller than the last step.(4)
Steps (2) and (3) are repeated, until the characteristic frequency and its bandwidth are identified. The target can be achieved usually after several optimization calculations.Figure 2
Program flow chart.The calculation procedure of quadratic optimized filter is depicted iconically in Figure3. The optimized result of Step 1 becomes the process object in Step 2. At the same time, the parameters of the Morlet wave are also adjusted. According to the result of Step 2, the fault signal was extracted in Step 3.Figure 3
Quadratic optimization flow chart.


## 2.1. Morlet Wavelet

The complex Morlet wavelet is the product of a complex exponential function multiplied by a Gaussian function. The expression of Morlet wavelet is(1)Morlt=σπe-σ2t2ei2πf0t.Meanwhile, the shape of Gaussian window in the frequency domain can be derived as [16](2)Morlf=Morl∗f=e-π2/σ2f-f02,where Morl(f) is the fast Fourier transformation of Morl(t). Since Morl(f) is real, Morl(f)=Morl∗(f), and the superscript ∗ denotes the complex conjugate. f0 is the center frequency of the wave and σ is waveform parameter that determines its width. For Morlet wavelet, the center frequency f0, and shape factor σ affect the shape and location simultaneously. By (1), the corresponding wave shape of Morlet wavelet is shown in Figure 1.Figure 1
Structure of Morlet wavelet.In Figure1, Morl1 is the wave of the Gaussian function, and Morl2 is the complex exponential wave function. They are expressed as(3)Morl1=e-t2σ2,Morl2=ei2πf0t.According to the Euler formula,Morl2 can be expressed as Morl2=cos⁡(2πf0t)+isin⁡(2πf0t), where Re(Morl2)=cos⁡(2πf0t). If t1 is half a cycle of Morl1 and t2 is a quarter period of Re(Morl2), t1 and t2 can be deduced as(4)t1=ln⁡m-1σ,t2=arccos⁡m2πf0,where m is the value of Morlet wave endpoint and it is infinitesimal (as shown in Figure 1). It is necessary to ensure that the Morlet wavelet is fluctuant enough and compactly supported. For Morl1 and Morl2, the period ratio T1/T2 is(5)T1T2=2t14t2=Qπln⁡m-1arccos⁡m=k,where Q is the quality factor, Q=f0/σ. When k1≤k≤k2, Morlet wavelet at most includes k2Morl2 and at least k1Morl2. If (5) is introduced into k1≤k≤k2, the expression of the quality coefficient Q can be obtained as follows:(6)k1arccos⁡mπln⁡m-1≤Q≤k2arccos⁡mπln⁡m-1.The quality coefficientQ is proportional to the fluctuation of the Morlet wavelet. The larger the quality coefficient Q is, the smaller the volatility of the Morlet wavelet is and vice versa. Therefore, when 0.5≤k≤2 and m=0.00001, the value range of Q is 0.0217≤Q≤0.0869.The corresponding family of wavelets consists of a series of son wavelets, which are generated by dilation and translation from the mother wavelet shown as follows:(7)Morla,bt=1aMorlt-ba,where a is the scale factor and b is wavelet displacement factor (the constant 1/a is used for energy normalization). The wavelet transform coefficient is defined by the following equation:(8)WTa,b=∫-∞+∞xt·Morla,b∗t·dt=xt,Morla,bt,where Morla,b∗(t) is the conjugate complex numbers of Morla,b(t) and · represents solving convolution.


## 2.2. Spline Curve Interpolation

Since the SK values are not continuous, the spline curve is used to obtain the variation trend of SK on the basis of discrete points. It can provide the basis for optimizing the analysis range.The function of spline curve interpolation is an approximation of piecewise three polynomials, but it is continuous. Its first- and second-order derivatives are both continuous. The value ofxj is yj=f(xj)(j=0,1,…,n), spline curve interpolation function is S(x), and its first-order derivative S′(xj)=mj. S(x) can be expressed as(9)Sx=∑j=0nyjαjx+mjβjx,where αj(x) and βj(x) are the interpolation function, expressed by(10)αjx=x-xj-1xj-xj-121+2x-xjxj-1-xjxj-1≤x≤xjj≠0x-xj+1xj-xj+121+2x-xjxj+1-xjxj≤x≤xj+1j≠n0Others,βjx=x-xj-1xj-xj-12x-xjxj-1≤x≤xjj≠0x-xj+1xj-xj+12x-xjxj≤x≤xj+1j≠n0Others.Under the condition of corresponding boundary, the functionS(x) can be calculated when the equations of mj(j=0,1,…,n) were obtained [17].


## 2.3. Transient Signal Detection Based on SK

The nonstationary signal is normally expressed by Wold-Cramér decomposition; that is,(11)Yt=∫-∞+∞ej2πftHt,fdXf,where the time-varying transfer function H(t,f) can be interpreted as the complex envelope or complex demodulate of process Y(t) at frequency f.The signal characteristics information inH(t,f) can be described by the spectral moments, the expression of which is(12)S2nYf≜ES2nYt,f=EHt,fdXf2ndf=EHt,f2n·S2nx.Under the conditions of ergodicity and stationarity ofH(t,f), it is easy to prove that(13)S2nYf=S2nYt,ft.The significant property of nonstationary processes is non-Gaussian. For a non-Gaussian process, spectrum accumulation that is greater than or equal to the fourth order is a nonzero value. Therefore, the fourth-order spectral accumulation is adopted, which is defined as(14)C4Yf=S4Yf-2S2Y2f.The greater the degree of signal deviation from Gaussianity is, the greater its fourth-order spectral accumulation is. Therefore, the energy normalized fourth-order spectral accumulation can be used to measure the peak of the signal process probability density at the frequencyf, that is, spectral kurtosis.The wavelet coefficients can be obtained by calculating the signal convolution by each row of the filter. The higher value of the wavelet coefficient is, the more similar the wavelet function with the signal waveform will be. The values of the wavelet coefficients are equal to the square of energy. SK is used to measure the energy spectrum. Therefore, SK achieves the maximum value when the waveforms are most similar [16]. Equation (15) gives the value of spectral kurtosis for each filter [18]:(15)kurtosisy≜C4YfS2Y2f=S4YfS2Y2f-2=meanH4t,fmeanH2t,f2-2.


## 2.4. Construction of Wavelet Clusters

Adaptive wavelet method uses spectral kurtosis as the evaluation index of Morlet wavelet filtered results. This method can optimize the latter filter parameters according to the calculation results in the previous step. Each step can improve the calculation accuracy and reduce the searching scope effectively. Compared with other bearing fault diagnosis method based on kurtosis, the accuracy of fault frequency identification and the detection bandwidth are both optimized.During the construction of octave filter with complex translation Morlet wavelet, the number of filters that cover the entire analysis band isM=K×N, where K is the number of the filters in each octave and N is the number of the octaves [19].The central frequency of the filter is expressed as(16)fmk=fmax-fmin21/Knn=1,2,…,M,where [fmin,fmax] is the frequency band range of the bearing characteristic frequencies in the operating mode. Bandwidth of each filter is expressed by(17)σnK=fnK·2K.The value ofK and N and the frequency analysis range [fmin,fmax] containing all possible fault frequencies are initially set. The center frequencies and bandwidths can be obtained by (16) and (17). The results are introduced into (1); then, N Morlet wavelet clusters for signal processing are obtained.When the frequency analysis range is optimized by the spline curve, the same number of filters is allocated within a smaller frequency range. As a result, the fault center frequency is more accurate and the bandwidth of the filter is smaller.


## 2.5. Adaptive Wavelet Method Based on Spectrum Diagnosis

Obviously, the proportion of the impact component in each filter window is different and it directly influences the spectral kurtosis of the filtering results. The former one is large means that the latter one is also large. The value of the quality factorQ affects the shape of the wavelet clusters. Therefore, each quality factor corresponds to an optimization process. For a rolling element bearing, the signal components related to the rotational speed usually cannot be extracted from the time-frequency representation directly. A new proposed method is presented to solve the above problems, as shown in Figure 2:(1)
In the frequency band range containing the bearing characteristic frequencies[fmin1,fmax1], it can calculate the spectral kurtosis of the filtered result and describe the spectral kurtosis by the spline curve interpolation along the frequency axis.(2)
By solving the derivation of the spline curve, the peak containing maximum spectral kurtosis and the valley on both sides of the maximum spectral kurtosis could be obtained, respectively, which is the frequency range of the next optimizing calculation[fmin2,fmax2].(3)
The quality factorQ is changed and the range of Morlet wavelet is reduced to [fmin2,fmax2]; the spectral kurtosis of the filter bandwidths can be calculated continually. The spline curve of the spectral kurtosis can be obtained by the spline curve interpolation. However, the bandwidth of the filter is smaller than the last step.(4)
Steps (2) and (3) are repeated, until the characteristic frequency and its bandwidth are identified. The target can be achieved usually after several optimization calculations.Figure 2
Program flow chart.The calculation procedure of quadratic optimized filter is depicted iconically in Figure3. The optimized result of Step 1 becomes the process object in Step 2. At the same time, the parameters of the Morlet wave are also adjusted. According to the result of Step 2, the fault signal was extracted in Step 3.Figure 3
Quadratic optimization flow chart.


## 3. Methods Validation

### 3.1. Simulation Analysis

The validity of proposed method is verified by the bearing fault simulation signal under strong noise background. The fault simulation signal including a natural frequency of 7000 Hz and a fault signal frequency of 307.6 Hz is generated. A white Gaussian noise of 25 dB is added to the signal and the noise ratio is 5.49, as shown in Figure4.Figure 4
(a) Simulation signal in time domain and (b) simulation signal in frequency domain.
(a) (b)The bearing fault simulation signal is detected by the above method. First of all, the spectral kurtosis is calculated along the initial pass band. The impact frequency is 268.8 Hz, andSKmax=310.8, as shown in Figure 5(a). Then, the impact frequency is 297.5 Hz, and SKmax=682.1 in the frequency range narrowed in the previous step, as shown in Figure 5(b). Thirdly, the frequency range is further reduced. The impact frequency is 307 Hz, and SKmax=5652, as shown in Figure 5(c). The optimization is close to the impact frequency gradually and the bandwidth is gradually reduced at each time.Figure 5
(a) The first optimization result of simulation signal. (b) The second optimization result of simulation signal. (c) The third optimization result of simulation signal.
(a) (b) (c)Finally, the bearing fault frequency is detected in 307.6 Hz and the bandwidth is about 200 Hz as shown in Figure6(a). Figure 6(b) shows the result of STFT-based SK [14]. The maximum spectral kurtosis (Kmax=0.3) is at Leave 4.5, the center frequency fc=208.33Hz, and the bandwidth BW=416.67Hz. The center frequency fc does not match better with impact frequency and the bandwidth BW is too wide to cover the feature frequency band precisely compared with the proposed method.Figure 6
(a) Impact signal in frequency domain. (b) Fast Kurtogram of simulation signal verification.
(a) (b)


### 3.2. Experimental Tests

This method is verified by detecting the test data of the Case Western Reserve University Bearing Data Center [20]. In the experiment process, the bearing (SKF6205) was seeded with faults by electrodischarge machining (EDM) previously. A 0.007-inch fault was introduced at the outer raceway of faulted bearings. Then, this bearing was reinstalled into the test motor. Vibration data was recorded for motor loads of 1 horsepower and motor speeds of 1772 RPM. The specific test data is shown in Table 1.Table 1
The specific test data.
Parameter Dimension of fault Shaft rotation speed Motor loads Sampling frequency Value 0.007 in 29.53 Hz 1 HP 12 kHzThe bearing characteristic frequency under this condition is shown in Table2. Therefore, the initial analysis band range is designed as [5,300] Hz. The recorded vibration signal is shown in Figure 7.Table 2
The bearing characteristic frequency.
Inner ring Outer ring Rolling element Cage 159.93 Hz 105.73 Hz 139.21 Hz 11.76 HzFigure 7
(a) The vibration signal in time domain. (b) The vibration signal in frequency domain.
(a) (b)The initial maximal spectral kurtosis and the impact frequency are 3104 and 131.3 Hz, respectively. Then, the spectral kurtosis is calculated in the optimized frequency band, and it can be extracted that the impact frequency is 99.72 Hz,SKmax=5959 in Figure 8.Figure 8
(a) The first optimization result of vibration signal. (b) The second optimization result of vibration signal.
(a) (b)The fault frequency is 106.2 Hz and the fault signal can be identified clearly in the frequency domain and time domain as shown in Figure9.Figure 9
(a) Bearing fault signal in frequency domain. (b) Bearing fault signal in time domain.
(a) (b)


## 3.1. Simulation Analysis

The validity of proposed method is verified by the bearing fault simulation signal under strong noise background. The fault simulation signal including a natural frequency of 7000 Hz and a fault signal frequency of 307.6 Hz is generated. A white Gaussian noise of 25 dB is added to the signal and the noise ratio is 5.49, as shown in Figure4.Figure 4
(a) Simulation signal in time domain and (b) simulation signal in frequency domain.
(a) (b)The bearing fault simulation signal is detected by the above method. First of all, the spectral kurtosis is calculated along the initial pass band. The impact frequency is 268.8 Hz, andSKmax=310.8, as shown in Figure 5(a). Then, the impact frequency is 297.5 Hz, and SKmax=682.1 in the frequency range narrowed in the previous step, as shown in Figure 5(b). Thirdly, the frequency range is further reduced. The impact frequency is 307 Hz, and SKmax=5652, as shown in Figure 5(c). The optimization is close to the impact frequency gradually and the bandwidth is gradually reduced at each time.Figure 5
(a) The first optimization result of simulation signal. (b) The second optimization result of simulation signal. (c) The third optimization result of simulation signal.
(a) (b) (c)Finally, the bearing fault frequency is detected in 307.6 Hz and the bandwidth is about 200 Hz as shown in Figure6(a). Figure 6(b) shows the result of STFT-based SK [14]. The maximum spectral kurtosis (Kmax=0.3) is at Leave 4.5, the center frequency fc=208.33Hz, and the bandwidth BW=416.67Hz. The center frequency fc does not match better with impact frequency and the bandwidth BW is too wide to cover the feature frequency band precisely compared with the proposed method.Figure 6
(a) Impact signal in frequency domain. (b) Fast Kurtogram of simulation signal verification.
(a) (b)


## 3.2. Experimental Tests

This method is verified by detecting the test data of the Case Western Reserve University Bearing Data Center [20]. In the experiment process, the bearing (SKF6205) was seeded with faults by electrodischarge machining (EDM) previously. A 0.007-inch fault was introduced at the outer raceway of faulted bearings. Then, this bearing was reinstalled into the test motor. Vibration data was recorded for motor loads of 1 horsepower and motor speeds of 1772 RPM. The specific test data is shown in Table 1.Table 1
The specific test data.
Parameter Dimension of fault Shaft rotation speed Motor loads Sampling frequency Value 0.007 in 29.53 Hz 1 HP 12 kHzThe bearing characteristic frequency under this condition is shown in Table2. Therefore, the initial analysis band range is designed as [5,300] Hz. The recorded vibration signal is shown in Figure 7.Table 2
The bearing characteristic frequency.
Inner ring Outer ring Rolling element Cage 159.93 Hz 105.73 Hz 139.21 Hz 11.76 HzFigure 7
(a) The vibration signal in time domain. (b) The vibration signal in frequency domain.
(a) (b)The initial maximal spectral kurtosis and the impact frequency are 3104 and 131.3 Hz, respectively. Then, the spectral kurtosis is calculated in the optimized frequency band, and it can be extracted that the impact frequency is 99.72 Hz,SKmax=5959 in Figure 8.Figure 8
(a) The first optimization result of vibration signal. (b) The second optimization result of vibration signal.
(a) (b)The fault frequency is 106.2 Hz and the fault signal can be identified clearly in the frequency domain and time domain as shown in Figure9.Figure 9
(a) Bearing fault signal in frequency domain. (b) Bearing fault signal in time domain.
(a) (b)


## 4. Conclusion

The method was based on the similarity between Morlet wavelet and faulty impulse signal. The spectral kurtosis maximum principle was employed into wavelet spectrum detection method for rolling element bearing fault signal as a criterion. At the initial filtering of the raw signal, the frequency range of the fault signal was found by the spline curve interpolation and its derivation. Different wavelet clusters that produced by the quality factorQ could optimize the searching range gradually and reduce the filter window width. The fault diagnosis ability of rolling bearing was improved. In conclusion, according to the treatment of the simulation signal and the bearing experimental test signal by the proposed method, it showed that the proposed method recognized the fault signal accurately and had better bearing fault time-frequency recognition feature.



---

*Source: 6106103-2017-02-19.xml*