# Rolling Bearing Fault Feature Extraction Based on Adaptive Tunable Q-Factor Wavelet Transform and Spectral Kurtosis

Authors: Jianlong Zhao; Yongchao Zhang; Qingguang Chen
Journal: Shock and Vibration
 (2020)

---

## Abstract

The fault feature of the rolling bearing is difficult to extract when weak fault occurs and interference exists. The tunable Q-factor wavelet transform (TQWT) can effectively extract the weak fault characteristic of the rolling bearing, but the manual selection of the Q-factor affects the decomposition result and only using TQWT presents interference. Aiming at the above problems, an adaptive tunable Q-factor wavelet transform (ATQWT) and spectral kurtosis (SK) method is proposed in this paper. Firstly, the method applies particle swarm optimization (PSO) to seek the optimized Q-factor for avoiding manual selection, which uses the kurtosis value of the transient impact component as the particle fitness function. The rolling bearing fault signal is decomposed into continuous oscillation component and transient impact component containing fault feature by the optimized Q-factor. Then, due to the presence of interference in the decomposition result of ATQWT, the SK analysis of the transient impact component is used to determine the frequency band of periodic impact component characterizing fault feature by fast kurtogram. Finally, the band-pass filter established through the above frequency band is employed to filter the interference in the transient impact component. Simulation and experimental results indicate that the ATQWT can highlight the periodic impact component characterizing rolling bearing fault feature, and the SK can filter interference in the transient impact component, which improves feature extraction effect and has great significance to enhance fault diagnosis accuracy of the rolling bearing. Compared with EEMD-TQWT and TQWT-SK, the fault feature extracted by the proposed method is prominent and effective.


---

## Body

## 1. Introduction

Rolling bearing, a harsh working environment and high incidence of failure, is a critical part of rotating machinery. When fault occurs, the measured vibration signal is the combination effect from its structure, machining and assembly errors, movement and force of other parts, and operating failure, etc. [1–5]. It contains abundant feature information and the excitation with different fault leads to distinct response. Therefore, it is the most important task in the rolling bearing fault diagnosis, extracting the fault feature from the measured vibration signal [6].At the moment of rolling bearing failure, the continuous contact between the failure surface and other surfaces generates a series of impact pulse force with broadband [7]. The resonance is excited on account of impact pulse force covered high-frequency natural vibration of the rolling bearing, and its attenuation time is much smaller than the interval between the impulse forces, thus generating a series of periodic impact components [8]. In the early stage of rolling bearing fault, the periodic impact component that characterizes the rolling bearing fault is weak, has interference, and is difficult to extract, so a lot of research has been conducted on the fault feature extraction method of rolling bearing in recent years. Time-frequency analysis methods are used to extract fault feature of rolling bearing such as short-time Fourier transform (STFT), Winger-Ville distribution (WVD), wavelet transform (WT), and empirical mode decomposition (EMD) [9]. The size of the time-frequency window in STFT does not change and lacks adaptability. When the WVD is used for multi-component signals, it is easy to generate cross term and difficult to extract feature frequency [10]. The WT needs to select the basis function and decomposition layer in advance, and the selection has a greater impact on the final decomposition result. EMD has a certain mode confusion [11].On account of the problems existing in the above feature extraction methods, Selesnick proposed the tunable Q-factor wavelet transform (TQWT), which can decompose the vibration signal into the continuous oscillation component and the transient impact component according to the difference of pre-set Q-factor [12, 13]. Wang used the TQWT to decompose the IMF of the maximum kurtosis value obtained by EEMD, and the low resonance component was carried out by envelope spectrum analysis to realize the rolling bearing fault diagnosis [14]. Li decomposed the vibration signal of rolling bearing through ICD and extracted the fault feature of rolling bearing using TQWT [15]. Bharath utilized TQWT to extract rolling bearing fault feature and realized the classification of rolling bearing fault via classification algorithm [16]. Consequently, the TQWT can extract the rolling bearing fault feature. Spectral kurtosis (SK) can detect the non-Gaussian component in the signal and provide its statistical indicator in the frequency position [17–19]. Abbas founded the resonance band of the defective rolling bearing through SK and achieved the rolling bearing fault identification [20]. Wan applied SK to select the resonance frequency band and extract the fault feature using the envelope demodulation method [21]. After the above analysis, we can find that the periodic impact component characterizing the rolling bearing fault feature is not obviously directly obtained by the TQWT, has a lot of interferences, and is influenced by the selected Q-factor. The SK can locate the frequency band of the periodic impact component owing to its non-Gaussian property, but the filtered rolling bearing fault signal has many interferences only using the above frequency band. Therefore, the fault feature extraction method needs to be further improved.In view of the advantages of TQWT and SK in extracting rolling bearing fault feature, this paper proposes a fault feature extraction method based on adaptive tunable Q-factor wavelet transform (ATQWT) and SK. The high and low factors are optimized by particle swarm optimization (PSO), and the optimized Q-factor is used to perform tunable Q-factor wavelet transform to highlight the periodic impact component in the transient impact component; then the SK analysis of the transient impact component determines the frequency range of the periodic impact component. The filter established by the above frequency band can filter out the interference component in the transient impact component, highlight the fault feature, and improve the accuracy of fault diagnosis of rolling bearing.


## 2. Fault Feature Extraction Method of Rolling Bearing

### 2.1. Tunable Q-Factor Wavelet Transform

To adapt to the feature of the analyzed signal, the Q-factor of the wavelet transform should be selected according to its property. It should have a high Q-factor in analyzing continuous oscillation signal; otherwise it should own low Q-factor in the transient impact signal, but its Q-factor is difficult to adjust. The TQWT proposed by Selesnick can select the Q-factor according to the characteristic of the analyzed signal and possess perfect reconstruction and moderate overcompleteness and is completely discrete [13].The decomposition and reconstruction filter bank of the TQWT is shown in Figure1 and a total of 4 sub-band signals are obtained, where fs is the sampling frequency, α is the low-pass scaling factor, and β is the high-pass scaling factor. Jj≤J is the decomposition layer, and the high-frequency and low-frequency sub-band signal obtained by the j-level filter bank are cj and wj, respectively, and the sampling rates are αfs and βαj−1fs, respectively. cj is the input of the j + 1 level filter bank. H0ω and H1ω are low-pass and high-pass filter response function. H0∗ and H1∗ are the complex conjugate of H0ω and H1ω.Figure 1
Decomposition and reconstructed filter bank of TQWT. (a) Decomposition filter bank. (b) Reconstructed filter bank.
(a)(b)In order to ensure that the wavelet transform is not excessively redundant and has perfect reconstruction property, the scale factor should satisfy(1)0<α<1;0<β≤1;α+β>1,where α is the low-pass scaling factor and β is the high-pass scaling factor.For perfect reconstruction, the low-pass filter and the high-pass filter should be constructed by(2)H0ω=1,θω+β−1πα+β−1,0,(3)H1ω=1,θαπ−ωα+β−1,0,(4)θω=0.5×1+cosω×2−cosω,ω≤π,where H0ω and H1ω are low-pass and high-pass filter response function, respectively, and θω is the 2π-periodic power complementary function.The center frequency and bandwidth of the decomposition layerj can be calculated by equation (5), and the scale factor αandβ can be calculated by equation (6):(5)fz=αj2−β4αfs,BW=12βαj−1fs,(6)β=2Q+1,α=1−βr,where fz is the center frequency, BW is the bandwidth, Q is the quality factor, r is the oversampling rate, and fs is the sampling frequency of the input signal.The specified Q-factor should be chosen subject toQ≥1. Setting Q=1 leads to a wavelet transform for which the wavelet resembles the second derivative of a Gaussian. Higher values of Q lead to more oscillatory wavelets. The specified oversampling rate r must be strictly greater than 1. If r is close to unity, then the transition bands of H0ω and H1ω will be relatively narrow and the time-domain response will not be well localized. For r≈1, the wavelet will resemble the sinc wavelet. In order to avoid this issue, it is sufficient to select r≥3.DifferentQ,r correspond to different decomposition layer J. It is determined by the maximum decomposition layer. The maximum decomposition layer can be obtained by(7)Jmax=logβN/8log1/α=logN/4/Q+1logQ+1/Q+1−2/r,where N is the length of the input signal and Jmax is the maximum decomposition layer J≤Jmax.The rolling bearing fault signal is composed of continuous oscillation component and transient impact component, and it is difficult to separate them owing to their overlap frequency band or existing interference. Morphological component analysis can effectively separate the continuous oscillation component and the transient impact component according to their oscillation property; the signalx is supposed to be a linear combination of high resonance component x1 and low resonance component x2 with morphological difference; then, x is expressed as(8)x=x1+x2x,…x1,x2∈RN.In order to achieve the separation ofx1 and x2 employing morphological component analysis, S1 and S2 are assumed as overcomplete dictionaries for signals x1 and x2, respectively, and x1 can only be represented sparsely by S1 and x2 can only be represented sparsely by S2, and w1 and w2 are the coefficient vectors of sparse representation; then, the solution to the sparse decomposition coefficient of the signal is transformed into an optimization problem as in(9)w1opt,w2opt=argminw1,w2w10+w20s.t.x=S1w1+S2w2.The solution of the resonance sparse decomposition coefficient is transformed into optimization problem of equation (10) through morphological component analysis:(10)Jw1,w2=x−S1w1−S2w222+λ1w11+λ2w21.By iteratively solving the decomposition coefficientw1∗ and w2∗ corresponding to the minimum value of equation (10), the high resonance component x1 and low resonance component x2 obtained by the resonance sparse decomposition are, respectively, expressed as(11)x1=S1w1∗,x2=S2w2∗.


### 2.2. Adaptive Tunable Q-Factor Wavelet Transform

The Q-factor of the TQWT should be selected according to the oscillation property of the decomposed signal. According to the oscillation property of the signal, the high Q-factor of TQWT is greater than or equal to 3, and the low Q-factor is between 1 and 3. The greater high Q-quality does not mean the better decomposition effect. The search space is relatively small, and long search time will cause casualty for fault diagnosis. Moreover, the PSO is simple to implement and has fast convergence speed and strong practicability. Considering the above factors, the PSO was selected among several algorithms such as genetic algorithm and ant colony algorithm, etc. The PSO is used to optimize the Q-factor of the tunable Q-factor wavelet transform. The kurtosis value of the low resonance component is used as a particle fitness function. After continuous iteration, the Q-factor adapted to the decomposed signal is finally obtained. The PSO is an algorithm for solving practical optimization problem obtained by studying the migration and foraging behavior of birds. To randomly initialize a group of particles without mass and volume, each particle is regarded as a potential solution of the optimization problem and continuously flying in the search space at a certain speed, and the particle can find the optimal solution through iteratively updating fitness [22].Suppose there arem particles in the n-dimensional target search space, the position and velocity of the i-th particle are yi=yi1,yi2,…,yin and vi=vi1,vi2,…,vin, and the i-th particle individual extremum is pi=pi1,pi2,…,pi3, and the global extremum of the entire population is pg=pg1,pg2,…,pg3. Due to the iterative update of the current particle velocity and position according to its velocity and position, individual extremum, and global extremum, the k + 1 iteration of the i-th particle is carried out according to(12)vidk+1=wvidk+c1×b1×pid−yidk+c2×b2×pgd−yidk.(13)yidk+1=yidk+vidk,where w is the inertial weight, i=1,2,…,m,m is the total number of particles, n is the dimension of the solution space, c1 and c2 are the learning factors, and b1 and b2 are the random numbers between (0, 1).


### 2.3. Spectral Kurtosis

Spectral kurtosis can detect Gaussian components contained in the signal [17] and can provide statistical indicators of its position in the frequency domain. It calculates the kurtosis of time domain data corresponding to each spectral line to determine the position of the transient impact components [23]. The SK is defined as follows.Spectral moment is a very valuable statistical indicator for measuring nonstationary processes, which is defined as(14)S2nYf≜ES2nYt,f=EHt,fdXf2n/df=EHt,f2nS2nX,where S2nYf is the spectral moment and Ht,f is the time-varying function.The calculation of the spectral moment needs to be repeated and the alternative method of spectral moment averages the instantaneous moment along the time axis, as shown in(15)S2nYft≜limT⟶∞1T∫−T/2T/2S2nYt,fdt,where …t represents the time average operator and T is time.Under the conditions of stationarity and ergodicity of the complex envelopeHt,f, equation (16) can be proved:(16)S2nYf=S2nYt,ft.The spectral cumulant is very sensitive to nonstationary processes. The non-Gaussian process’s2nn≥2 order spectral cumulant is nonzero, and the fourth-order spectral cumulant is defined as(17)C4Yf=S4Yf−2S2Y2f,f≠0,where C4Yf is the fourth-order spectral cumulant.The fourth-order spectral accumulation is larger when the signal much more deviates the Gaussian. The energy-normalized fourth-order spectral accumulation will give the peak measure of the probability density function of the signal at frequencyf, also the so-called SK. It is shown in(18)KYf≜C4YfS2Y2f=S4YfS2Y2f−2,f≠0,where KYf is the spectral kurtosis.


### 2.4. Fault Feature Extraction Based on Adapted TQWT and SK

Using kurtosis as the objective function to optimize the Q-factor through PSO, the tunable Q-factor wavelet transform can highlight the periodic impact component that characterizes the fault of the rolling bearing in the low resonance component. The SK can locate the frequency band of the periodic impact component. The filter established by the frequency band can filter out interference component and highlight fault feature. The rolling bearing fault type can be determined by the envelope spectrum analysis. The algorithm flow chart is shown in Figure2. The specific algorithm is as follows:(1)
InitializingQ1i,Q2i, m, Loopcount, xiandvi, the fitness function is established through the kurtosis of the low resonance component of the tunable quality wavelet transform(2)
Setting the redundancy factorr1 and r2 and the decomposition layer J1 and J2, the Q-factor of the tunable factor wavelet transform is optimized by PSO(3)
The rolling bearing fault signal is decomposed by tunable quality wavelet transform using optimized Q-factor, and the high and low resonance components are acquired(4)
The fast kurtogram gained by SK of the low resonance component finds out central frequencyfc and bandwidth BW of the periodic impact component characterizing rolling bearing fault feature(5)
The finite impulse response filter established by the central frequencyfc and bandwidth BW filter the low resonance component(6)
The envelope spectrum of the filtered signal is gained by Hilbert envelope spectrum analysis, and the rolling bearing fault type is identified by the amplitude distribution in envelope spectrumFigure 2
The flow chart of ATQWT and SK.


## 2.1. Tunable Q-Factor Wavelet Transform

To adapt to the feature of the analyzed signal, the Q-factor of the wavelet transform should be selected according to its property. It should have a high Q-factor in analyzing continuous oscillation signal; otherwise it should own low Q-factor in the transient impact signal, but its Q-factor is difficult to adjust. The TQWT proposed by Selesnick can select the Q-factor according to the characteristic of the analyzed signal and possess perfect reconstruction and moderate overcompleteness and is completely discrete [13].The decomposition and reconstruction filter bank of the TQWT is shown in Figure1 and a total of 4 sub-band signals are obtained, where fs is the sampling frequency, α is the low-pass scaling factor, and β is the high-pass scaling factor. Jj≤J is the decomposition layer, and the high-frequency and low-frequency sub-band signal obtained by the j-level filter bank are cj and wj, respectively, and the sampling rates are αfs and βαj−1fs, respectively. cj is the input of the j + 1 level filter bank. H0ω and H1ω are low-pass and high-pass filter response function. H0∗ and H1∗ are the complex conjugate of H0ω and H1ω.Figure 1
Decomposition and reconstructed filter bank of TQWT. (a) Decomposition filter bank. (b) Reconstructed filter bank.
(a)(b)In order to ensure that the wavelet transform is not excessively redundant and has perfect reconstruction property, the scale factor should satisfy(1)0<α<1;0<β≤1;α+β>1,where α is the low-pass scaling factor and β is the high-pass scaling factor.For perfect reconstruction, the low-pass filter and the high-pass filter should be constructed by(2)H0ω=1,θω+β−1πα+β−1,0,(3)H1ω=1,θαπ−ωα+β−1,0,(4)θω=0.5×1+cosω×2−cosω,ω≤π,where H0ω and H1ω are low-pass and high-pass filter response function, respectively, and θω is the 2π-periodic power complementary function.The center frequency and bandwidth of the decomposition layerj can be calculated by equation (5), and the scale factor αandβ can be calculated by equation (6):(5)fz=αj2−β4αfs,BW=12βαj−1fs,(6)β=2Q+1,α=1−βr,where fz is the center frequency, BW is the bandwidth, Q is the quality factor, r is the oversampling rate, and fs is the sampling frequency of the input signal.The specified Q-factor should be chosen subject toQ≥1. Setting Q=1 leads to a wavelet transform for which the wavelet resembles the second derivative of a Gaussian. Higher values of Q lead to more oscillatory wavelets. The specified oversampling rate r must be strictly greater than 1. If r is close to unity, then the transition bands of H0ω and H1ω will be relatively narrow and the time-domain response will not be well localized. For r≈1, the wavelet will resemble the sinc wavelet. In order to avoid this issue, it is sufficient to select r≥3.DifferentQ,r correspond to different decomposition layer J. It is determined by the maximum decomposition layer. The maximum decomposition layer can be obtained by(7)Jmax=logβN/8log1/α=logN/4/Q+1logQ+1/Q+1−2/r,where N is the length of the input signal and Jmax is the maximum decomposition layer J≤Jmax.The rolling bearing fault signal is composed of continuous oscillation component and transient impact component, and it is difficult to separate them owing to their overlap frequency band or existing interference. Morphological component analysis can effectively separate the continuous oscillation component and the transient impact component according to their oscillation property; the signalx is supposed to be a linear combination of high resonance component x1 and low resonance component x2 with morphological difference; then, x is expressed as(8)x=x1+x2x,…x1,x2∈RN.In order to achieve the separation ofx1 and x2 employing morphological component analysis, S1 and S2 are assumed as overcomplete dictionaries for signals x1 and x2, respectively, and x1 can only be represented sparsely by S1 and x2 can only be represented sparsely by S2, and w1 and w2 are the coefficient vectors of sparse representation; then, the solution to the sparse decomposition coefficient of the signal is transformed into an optimization problem as in(9)w1opt,w2opt=argminw1,w2w10+w20s.t.x=S1w1+S2w2.The solution of the resonance sparse decomposition coefficient is transformed into optimization problem of equation (10) through morphological component analysis:(10)Jw1,w2=x−S1w1−S2w222+λ1w11+λ2w21.By iteratively solving the decomposition coefficientw1∗ and w2∗ corresponding to the minimum value of equation (10), the high resonance component x1 and low resonance component x2 obtained by the resonance sparse decomposition are, respectively, expressed as(11)x1=S1w1∗,x2=S2w2∗.


## 2.2. Adaptive Tunable Q-Factor Wavelet Transform

The Q-factor of the TQWT should be selected according to the oscillation property of the decomposed signal. According to the oscillation property of the signal, the high Q-factor of TQWT is greater than or equal to 3, and the low Q-factor is between 1 and 3. The greater high Q-quality does not mean the better decomposition effect. The search space is relatively small, and long search time will cause casualty for fault diagnosis. Moreover, the PSO is simple to implement and has fast convergence speed and strong practicability. Considering the above factors, the PSO was selected among several algorithms such as genetic algorithm and ant colony algorithm, etc. The PSO is used to optimize the Q-factor of the tunable Q-factor wavelet transform. The kurtosis value of the low resonance component is used as a particle fitness function. After continuous iteration, the Q-factor adapted to the decomposed signal is finally obtained. The PSO is an algorithm for solving practical optimization problem obtained by studying the migration and foraging behavior of birds. To randomly initialize a group of particles without mass and volume, each particle is regarded as a potential solution of the optimization problem and continuously flying in the search space at a certain speed, and the particle can find the optimal solution through iteratively updating fitness [22].Suppose there arem particles in the n-dimensional target search space, the position and velocity of the i-th particle are yi=yi1,yi2,…,yin and vi=vi1,vi2,…,vin, and the i-th particle individual extremum is pi=pi1,pi2,…,pi3, and the global extremum of the entire population is pg=pg1,pg2,…,pg3. Due to the iterative update of the current particle velocity and position according to its velocity and position, individual extremum, and global extremum, the k + 1 iteration of the i-th particle is carried out according to(12)vidk+1=wvidk+c1×b1×pid−yidk+c2×b2×pgd−yidk.(13)yidk+1=yidk+vidk,where w is the inertial weight, i=1,2,…,m,m is the total number of particles, n is the dimension of the solution space, c1 and c2 are the learning factors, and b1 and b2 are the random numbers between (0, 1).


## 2.3. Spectral Kurtosis

Spectral kurtosis can detect Gaussian components contained in the signal [17] and can provide statistical indicators of its position in the frequency domain. It calculates the kurtosis of time domain data corresponding to each spectral line to determine the position of the transient impact components [23]. The SK is defined as follows.Spectral moment is a very valuable statistical indicator for measuring nonstationary processes, which is defined as(14)S2nYf≜ES2nYt,f=EHt,fdXf2n/df=EHt,f2nS2nX,where S2nYf is the spectral moment and Ht,f is the time-varying function.The calculation of the spectral moment needs to be repeated and the alternative method of spectral moment averages the instantaneous moment along the time axis, as shown in(15)S2nYft≜limT⟶∞1T∫−T/2T/2S2nYt,fdt,where …t represents the time average operator and T is time.Under the conditions of stationarity and ergodicity of the complex envelopeHt,f, equation (16) can be proved:(16)S2nYf=S2nYt,ft.The spectral cumulant is very sensitive to nonstationary processes. The non-Gaussian process’s2nn≥2 order spectral cumulant is nonzero, and the fourth-order spectral cumulant is defined as(17)C4Yf=S4Yf−2S2Y2f,f≠0,where C4Yf is the fourth-order spectral cumulant.The fourth-order spectral accumulation is larger when the signal much more deviates the Gaussian. The energy-normalized fourth-order spectral accumulation will give the peak measure of the probability density function of the signal at frequencyf, also the so-called SK. It is shown in(18)KYf≜C4YfS2Y2f=S4YfS2Y2f−2,f≠0,where KYf is the spectral kurtosis.


## 2.4. Fault Feature Extraction Based on Adapted TQWT and SK

Using kurtosis as the objective function to optimize the Q-factor through PSO, the tunable Q-factor wavelet transform can highlight the periodic impact component that characterizes the fault of the rolling bearing in the low resonance component. The SK can locate the frequency band of the periodic impact component. The filter established by the frequency band can filter out interference component and highlight fault feature. The rolling bearing fault type can be determined by the envelope spectrum analysis. The algorithm flow chart is shown in Figure2. The specific algorithm is as follows:(1)
InitializingQ1i,Q2i, m, Loopcount, xiandvi, the fitness function is established through the kurtosis of the low resonance component of the tunable quality wavelet transform(2)
Setting the redundancy factorr1 and r2 and the decomposition layer J1 and J2, the Q-factor of the tunable factor wavelet transform is optimized by PSO(3)
The rolling bearing fault signal is decomposed by tunable quality wavelet transform using optimized Q-factor, and the high and low resonance components are acquired(4)
The fast kurtogram gained by SK of the low resonance component finds out central frequencyfc and bandwidth BW of the periodic impact component characterizing rolling bearing fault feature(5)
The finite impulse response filter established by the central frequencyfc and bandwidth BW filter the low resonance component(6)
The envelope spectrum of the filtered signal is gained by Hilbert envelope spectrum analysis, and the rolling bearing fault type is identified by the amplitude distribution in envelope spectrumFigure 2
The flow chart of ATQWT and SK.


## 3. Simulation Analysis

In order to verify the effectiveness of the ATQWT and SK in rolling bearing fault feature extraction, a simulation signalx is simulated to the rolling bearing fault vibration signal, where x1 with a frequency of 100 Hz simulates the periodic impact component characterizing rolling bearing fault feature, and the sinusoidal signals with a frequency of 30 Hz and 2358 Hz express the rotational frequency and harmonic, and the Gaussian white noise w is used as the background noise with the intensity of −9db:(19)x1t=∑i=1MAie−1200tsin4800πt+π2,Ai=0.1cos20πt+π2+0.2,xt=x1+0.05sin60πt+0.2sin4716πt+w.Setting the sampling frequency to 10 kHz and the sampling points to 2048,x1 and x are shown in Figures 3(a) and 3(b), respectively. The periodic impact component representing the fault feature of rolling bearing was submerged, and it was impossible to judge whether the rolling bearings failed in Figure 3(b). As shown in Figure 4, the envelope spectrum of x1 is most prominent in 100 Hz and its frequency doubling in Figure 4(a), which is consistent with the rolling bearing fault characteristic frequency, but the envelope spectrum of x has no noteworthy spectral lines in Figure 4(b), so it is not possible to directly apply Hilbert envelope to extract rolling bearing fault feature frequency.Figure 3
Periodic impact signalx1 and simulation signal x.
(a)(b)Figure 4
Hilbert envelope spectrum ofx1 and x.
(a)(b)The Q-factor optimized by PSO isQ1=11.4andQ2=1.8, and the rolling bearing fault signal is carried out ATQWT at parameters Q1=11.4,Q2=1.8,r1=r2=3,J1=32,J2=13, and the result is shown in Figure 5. In Figure 5(b), the periodic impact component in the low resonance component is not prominent and has interference, and its Hilbert envelope spectrum is presented in Figure 6. The 100 Hz and its frequency doubling emerge in Figure 6, but its amplitude is small and possesses interference in the frequency domain. Therefore, the periodic impact component in the low resonance component acquired by ATQWT is not prominent, has interference, and affects the accuracy of the fault diagnosis of the rolling bearing.Figure 5
The high and low resonance component ofx.
(a)(b)Figure 6
The Hilbert envelope spectrum of low resonance component ofx.The low resonance component is performed by the SK analysis, and the fast kurtogram is shown in Figure7. Its peak value presents in fc=2500Hz,Δf=1666.7Hz, and the above frequency band includes the 2000 Hz high-frequency natural vibration characterizing faults rolling bearing fault, so the SK can accurately locate the frequency band of the periodic impact component characterizing rolling bearing fault. Building a finite impulse response filter Hf,Δf with f=2500Hz,Δf=1666.7Hz, order = 30, the low resonance component is filtered by the above Hf,Δf, and the time domain and envelope spectrum of filter signal is shown in Figure 8. The periodic impact component in Figure 8(a) is prominent and the fault characteristic frequency of 100 Hz and its multiple frequency in Figure 8(b) have same outcome. Hence, the ATQWT and SK can extract the fault feature of rolling bearing, filter out the interference, and highlight the fault feature.Figure 7
Fast kurtogram of low resonance component ofx.Figure 8
Low resonance component and its envelope spectrum ofx after filtering.
(a)(b)In order to highlight the advantages of the proposed method in rolling bearing fault feature extraction, the decomposition results of the TQWT-SK and EEMD-TQWT methods are shown in Figures9 and 10; it can be found that the 100 Hz and its multiple frequency are not prominent and have very small amplitude, so the effect of the rolling bearing fault feature extraction proposed in this paper is significantly better than TWQT-SK and EEMD-TQWT.Figure 9
TQWT-SK decomposition result.
(a)(b)Figure 10
EEMD-TQWT decomposition result.
(a)(b)


## 4. Experimental Verification

To prove the effectiveness and feasibility of ATQWT and SK method in rolling bearing fault feature extraction, the bearing-rotor testbed is established and is presented in Figure11. The fault rolling bearing is made by electrical discharge machining and then using it to drill 0.2 mm diameter holes simulating the inner and outer race fault, respectively. The bearing designation is 6308, and the specific parameters are shown in Table 1. The inner and outer race fault signals of the rolling bearing are, respectively, collected in the rotational frequency fr=24.2Hz, the sampling frequency of 10 kHz, and the sampling points of 2048. According to the calculation formula of the fault characteristic frequency and the data in the Table 1, the theoretical fault characteristic frequency of the inner and outer race of the rolling bearing is fi=115.2Hz and fo=78.1Hz.Figure 11
Bearing-rotor testbed: (a) electromotor; (b) signal conditioner; (c) bearing pedestal; (d) piezoelectric accelerometer.Table 1
Rolling bearing parameters.
Outer race diameter (mm)Inner race diameter (mm)Pitch diameter (mm)Ball diameter (mm)Ball numberContact angle90406512.580


### 4.1. Rolling Bearing Inner Race Fault

Figure12 shows the fault signal and the envelope spectrum of the inner race fault; there is no periodic impact component in Figure 12(a) and no prominent spectral line of fi in Figure 12(b); it does not judge the rolling bearing fault type. The Q-factors Q1=8.5 and Q2=1 optimized by PSO are used to conduct ATQWT of the inner race fault signal and other parameters are set to r1=r2=3 and J1=32,J2=13. The decomposition result is shown in Figure 13, and there exist impact component and a lot of interference in the low resonance component. Figure 14 is not prominent in fi and its multiple frequency; therefore, directly performing ATQWT on the signal cannot accurately extract the periodic impact component of the inner race fault.Figure 12
Fault signal and envelope spectrum of inner race fault.
(a)(b)Figure 13
High and low resonance component of inner race fault.
(a)(b)Figure 14
Envelope spectrum of low resonance component of inner race fault.Figure15 shows the fast kurtogram of low resonance component, its peak present in the central frequency of 3750 Hz and bandwidth of 2500 Hz; it indicates that the periodic impact component is contained in the frequency band. The low resonance component is filtered by a finite impulse response filter Hf,Δf with f=3750Hz,Δf=2500Hz, order = 30, and its result is presented in Figure 16, the periodic impact component in Figure 16(a) and the spectral lines of fi,2fi,3fi in Figure 16(b) are prominent, and the fault feature of inner race was extracted successfully.Figure 15
Fast kurtogram of low resonance component of inner race fault.Figure 16
Fault signal and envelope spectrum of inner race after filter.
(a)(b)The decomposition results of the TQWT-SK and EEMD-TQWT in inner race fault are shown in Figures17 and 18; it can be found that the time domain waveform has many periodic impact components, but fi and its multiple frequency are not prominent and have small amplitude and present a series of fr modulation frequency. So the feature extraction effect of the inner race fault proposed in this paper is significantly better than TWQT-SK and EEMD-TQWT.Figure 17
TQWT-SK decomposition result of inner race fault.
(a)(b)Figure 18
EEMD-TQWT decomposition result of inner race fault.
(a)(b)


### 4.2. Rolling Bearing Outer Race Fault

The time domain waveform and envelope spectrum of the outer race fault are shown in Figure19, the periodic impact component in the time domain waveform is obliterated, and no prominent spectral lines are found in the envelope spectrum, and the failure state cannot be determined. Setting the parameters of Q1=4.1,Q2=1.1,r1=r2=3,J1=32,J2=13 performs the ATQWT to obtain high and low resonance components as shown in Figure 20, where Q-factor is obtained by the PSO; there presents periodic impact component and exists a series of interference. The envelope spectrum of the low resonance component in Figure 21 is prominent in fo−fr; it affects the accuracy of the rolling bearing fault diagnosis.Figure 19
Fault signal and envelope spectrum of outer race fault.
(a)(b)Figure 20
High and low resonance component of outer race fault.
(a)(b)Figure 21
Envelope spectrum of low resonance component of outer race fault.The peak value in Figure22 presents in the central frequency of 3125 Hz and bandwidth of 1250 Hz; it indicates that the periodic impact component is contained in the frequency band. It is presented in Figure 23 that time domain waveform and envelope spectrum of the low resonance component filtered by Hf,Δf with f=3125Hz,Δf=1250Hz, order = 30 and the periodic impact component in Figure 23(a) and the spectral lines of fo,2fo,3fo in Figure 16(b) are prominent; the fault feature of outer race was extracted successfully.Figure 22
Fast kurtogram of low resonance component of outer race fault.Figure 23
Fault signal and envelope spectrum of outer race after filter.
(a)(b)As present in Figures24 and 25, the outer race fault decomposition result of TQWT-SK and EEMD-TQWT presents a series of harmonic interference, and the periodic impact periodic component is not prominent and submerged. The frequency domain has frequency modulation and has a small amplitude in outer fault feature frequency fo and its multiple frequency. Therefore, the feature extraction effect of the outer race fault proposed in this paper is significantly better than TWQT-SK and EEMD-TQWT.Figure 24
TQWT-SK decomposition result of outer race fault.
(a)(b)Figure 25
EEMD-TQWT decomposition result of outer race fault.
(a)(b)


## 4.1. Rolling Bearing Inner Race Fault

Figure12 shows the fault signal and the envelope spectrum of the inner race fault; there is no periodic impact component in Figure 12(a) and no prominent spectral line of fi in Figure 12(b); it does not judge the rolling bearing fault type. The Q-factors Q1=8.5 and Q2=1 optimized by PSO are used to conduct ATQWT of the inner race fault signal and other parameters are set to r1=r2=3 and J1=32,J2=13. The decomposition result is shown in Figure 13, and there exist impact component and a lot of interference in the low resonance component. Figure 14 is not prominent in fi and its multiple frequency; therefore, directly performing ATQWT on the signal cannot accurately extract the periodic impact component of the inner race fault.Figure 12
Fault signal and envelope spectrum of inner race fault.
(a)(b)Figure 13
High and low resonance component of inner race fault.
(a)(b)Figure 14
Envelope spectrum of low resonance component of inner race fault.Figure15 shows the fast kurtogram of low resonance component, its peak present in the central frequency of 3750 Hz and bandwidth of 2500 Hz; it indicates that the periodic impact component is contained in the frequency band. The low resonance component is filtered by a finite impulse response filter Hf,Δf with f=3750Hz,Δf=2500Hz, order = 30, and its result is presented in Figure 16, the periodic impact component in Figure 16(a) and the spectral lines of fi,2fi,3fi in Figure 16(b) are prominent, and the fault feature of inner race was extracted successfully.Figure 15
Fast kurtogram of low resonance component of inner race fault.Figure 16
Fault signal and envelope spectrum of inner race after filter.
(a)(b)The decomposition results of the TQWT-SK and EEMD-TQWT in inner race fault are shown in Figures17 and 18; it can be found that the time domain waveform has many periodic impact components, but fi and its multiple frequency are not prominent and have small amplitude and present a series of fr modulation frequency. So the feature extraction effect of the inner race fault proposed in this paper is significantly better than TWQT-SK and EEMD-TQWT.Figure 17
TQWT-SK decomposition result of inner race fault.
(a)(b)Figure 18
EEMD-TQWT decomposition result of inner race fault.
(a)(b)


## 4.2. Rolling Bearing Outer Race Fault

The time domain waveform and envelope spectrum of the outer race fault are shown in Figure19, the periodic impact component in the time domain waveform is obliterated, and no prominent spectral lines are found in the envelope spectrum, and the failure state cannot be determined. Setting the parameters of Q1=4.1,Q2=1.1,r1=r2=3,J1=32,J2=13 performs the ATQWT to obtain high and low resonance components as shown in Figure 20, where Q-factor is obtained by the PSO; there presents periodic impact component and exists a series of interference. The envelope spectrum of the low resonance component in Figure 21 is prominent in fo−fr; it affects the accuracy of the rolling bearing fault diagnosis.Figure 19
Fault signal and envelope spectrum of outer race fault.
(a)(b)Figure 20
High and low resonance component of outer race fault.
(a)(b)Figure 21
Envelope spectrum of low resonance component of outer race fault.The peak value in Figure22 presents in the central frequency of 3125 Hz and bandwidth of 1250 Hz; it indicates that the periodic impact component is contained in the frequency band. It is presented in Figure 23 that time domain waveform and envelope spectrum of the low resonance component filtered by Hf,Δf with f=3125Hz,Δf=1250Hz, order = 30 and the periodic impact component in Figure 23(a) and the spectral lines of fo,2fo,3fo in Figure 16(b) are prominent; the fault feature of outer race was extracted successfully.Figure 22
Fast kurtogram of low resonance component of outer race fault.Figure 23
Fault signal and envelope spectrum of outer race after filter.
(a)(b)As present in Figures24 and 25, the outer race fault decomposition result of TQWT-SK and EEMD-TQWT presents a series of harmonic interference, and the periodic impact periodic component is not prominent and submerged. The frequency domain has frequency modulation and has a small amplitude in outer fault feature frequency fo and its multiple frequency. Therefore, the feature extraction effect of the outer race fault proposed in this paper is significantly better than TWQT-SK and EEMD-TQWT.Figure 24
TQWT-SK decomposition result of outer race fault.
(a)(b)Figure 25
EEMD-TQWT decomposition result of outer race fault.
(a)(b)


## 5. Conclusion

In order to overcome the difficulty of fault feature extraction of the rolling bearing, the selection of Q-quality, and existing interference in TQWT, the ATQWT and SK extraction method is proposed. The method applies PSO to seek the optimized Q-factor of TQWT and the SK analysis is used to determine the frequency band of periodic impact component characterizing fault feature by fast kurtogram. Its effectiveness is verified through simulation and experiment, and the following conclusions are obtained:(1)
The Q-factor optimization through PSO can be adapted to the decomposed signal, highlighting the periodic impact component characterizing the rolling bearing fault feature in the low resonance component(2)
The SK analysis of the low resonance component can accurately locate the frequency band of the periodic impact component, and the filter established by this band can filter out the interference in the rolling bearing fault signal(3)
The method based on ATQWT and SK can extract weak fault feature, improve the feature extraction effect, and enhance the fault diagnosis accuracy of rolling bearing



---

*Source: 8875179-2020-07-30.xml*