# A 1.5D Spectral Kurtosis-Guided TQWT Method and Its Application in Bearing Fault Detection

Authors: Xiong Zhang; Ming Zhang; Shuting Wan; Rujiang Hao; Yuling He; Xiaolong Wang
Journal: Shock and Vibration
 (2021)

---

## Abstract

Bearings are the key parts of rotating machinery, and their operation status is related to the operation safety of the whole equipment. Vibration signals often contain periodic impulse components which can reflect the fault state of bearings. However, due to the interference of signal transmission path and the influence of operating environment noise, the periodic impulse components in the signal are often submerged by the nonperiodic transient impulse components, modulation harmonic components, and noise components. Therefore, the core problem of bearing fault diagnosis theory is used to accurately extract the frequency band of bearing fault state information and suppress the frequency band of interference information. In this paper, the signal is processed by the tunable Q-factor wavelet transform (TQWT), the midfrequency band of the signal is tightly divided by selecting different Q-values, and the 1.5D spectral kurtosis defined in frequency domain is used to select the optimal subband. Simulated analysis shows that this method can avoid low-frequency harmonic interference, nonperiodic transient impulse components, and strong noise components in the time domain. Therefore, it can effectively realize the selection of the subbands of periodic impulse components and effectively extract fault feature information. Through experimental signal analysis, TQWT has good sparsity decomposition characteristics and can reasonably divide the signal frequency band, so as to separate the useful fault characteristic frequency band and interference frequency band. At the same time, compared with the kurtosis index, 1.5D spectral kurtosis has better robustness and resolution for low signal-to-noise ratio signals, which can achieve the purpose of fault characteristic band extraction.


---

## Body

## 1. Introduction

Rolling bearing is an important component of rotating machinery, and its operation status is directly related to whether the whole equipment can run safely and smoothly. Therefore, the fault analysis of rolling bearing has become a hot issue in the field of signal analysis [1–3]. Considering the problems of signal transmission path, signal acquisition convenience, and cost, vibration signal is a widely used analysis medium. So many analysis and processing methods based on vibration signal are derived [4–6]. Especially, in recent years, with the development of mathematical model construction and parameter optimization method of processing algorithm [7–9], signal processing algorithm has been greatly improved in terms of accuracy, efficiency, and adaptability.The key link of bearing vibration signal processing is filtering process. Through certain analysis methods, the fault information within the specific frequency band can be retained to the maximum extent, and the interference components outside the frequency band can be suppressed to the maximum extent, so as to realize the signal filtering and fault information extraction. There are two main problems: one is to use appropriate methods to segment the frequency band, which makes the subband after segmentation have certain sparsity, and the other is to use appropriate index to evaluate the subband after segmentation, so that the subband with fault impulse component has certain specificity. The kurtogram algorithm [10] proposed by Antoni is considered to be a typical signal processing method, which can divide the signal plane into subbands with different bandwidth and center frequency through filter banks and select the subbands by the kurtosis index. Subsequently, many scholars began to improve the frequency band segmentation process and index selection of this kind of algorithm. Lei et al. [11] used wavelet packet decomposition to replace the equivalent filter banks in the kurtogram algorithm. Wang and Tsui [12] used dynamic Bayes to optimize the wavelet parameters and adjust the center frequency and bandwidth of the equivalent filtering characteristics, so that it can more accurately locate the fault impact subband. Qin et al. [13] constructed a new M-band flexible wavelet transform to identify the potential fault features in the measured signal. This transform provides more intensive time-frequency plane sampling and retains the tunable filter parameters and dilation factors. The simulated and experimental results show that this transform has better denoising performance than analytic flexible wavelet transform, orthogonal wavelet transform, and biorthogonal wavelet transform. Meanwhile, Qin et al. [14] proposed a novel K-SVD algorithm to get the optimized transient dictionary for the one-dimensional signal. Simulated and experimental results show that this method has advantages over the traditional K-SVD method and infogram method in periodic transient impulse components detection. The kurtosis index is too sensitive to impulse components and is easily disturbed by nonperiodic transient impulses and noise. Based on these considerations, Antoni further defines the square envelope entropy index and the square envelope spectral entropy index in the time domain and frequency domain, respectively, and constructs the infogram algorithm [15]. Li et al. [16] proposed a multiscale clustering gray infogram, which combined two negative entropies with the gray method. Moshrefzadeh and Fasana [17] proposed an autogram method using square envelope unbiased autocorrelation instead of time-domain signal. Through this operation, uncorrelated interference components can be discarded.Selesnick designed the tunable Q-factor wavelet transform (TQWT) from the perspective of frequency domain filtering [18]. TQWT can make the sparsity of frequency band division more diversified by adjusting its Q-factor. The key problem of TQWT is how to determine its parameter value and make the selection of frequency band more accurate by adjusting the oscillation characteristics of wavelet. Zhao et al. [19] used the particle swarm optimization algorithm to optimize TQWT parameters, combined with the spectral kurtosis index to select the optimal subband. Ma et al. [20] also used particle swarm optimization to optimize the parameters of TQWT and further combined with frequency slice wavelet transform (FSWT) to select the frequency band. Ding et al. [21] constructed a filter bank with almost constant bandwidth by adjusting the Q-factor and decomposition level of TQWT.In this paper, a 1.5D spectral kurtosis-guided TQWT method is proposed. In this method, the 1.5D spectral kurtosis of the square envelope sequence of wavelet subband signals with different Q-values and scales is calculated to replace the traditional kurtosis index to select the optimal subband. This method can suppress the noise and the nonperiodic transient impulse interference without coupling relationship.The outline of this paper is as follows. Section2 introduces the 1.5D spectral kurtosis-guided TQWT method and its superiority in signal periodic feature extraction. Sections 3 and 4 verify the effectiveness of the proposed method with the aid of simulated and experimental data.


## 2. 1.5D Spectral Kurtosis-Guided TQWT

### 2.1. TQWT

Compared with the traditional discrete wavelet transform, TQWT can determine the wavelet basis function by designing the Q-factor (the Q-factor is the ratio of signal center frequency to bandwidth,Q=fc/Bω) and the redundancy factor according to the signal oscillation characteristics. This design makes the choice of wavelet function more flexible, makes the signal adaptability stronger, and makes it easier to control the sparsity of frequency band division. TQWT uses band-pass filter banks to realize signal decomposition and reconstruction in an iterative manner. The three-level decomposition and reconstruction is taken as an example. The corresponding dual channel filter banks are shown in Figure 1.Figure 1
The corresponding dual channel filter banks of (a) decomposition process and (b) reconstruction process.
(a)(b)In Figure1, H1ω and H0ω are the frequency response functions of the high-pass filter and the low-pass filter (H1ω and H0ω can be expressed as equations (1) and (2)), H1∗ω and H0∗ω are the conjugate complex numbers of H1ω and H0ω, respectively, β=2/Q+1 is the high-pass scaling factor, and α=1−β/r is the low-pass scaling factor, where r is the redundancy factor:(1)H0ω=1,ω≤1−β,θω+β−1πα+β−1,1−βπ≤ω<απ,0,απ≤ω<π,(2)H1ω=0,ω≤1−β,θαπ−ωα+β−1,1−βπ≤ω<απ,1,απ≤ω<π.The maximum number of decomposition layers of TQWT can be obtained by(3)Jmax=logβNs/8log1/α.The scale factorsα and β can be substituted into equation (3) to obtain(4)Jmax=logNs/4Q+1logQ+1/Q+1−2/r.In summary, the filtering characteristics of TQWT change with the Q-factor and the redundancy factorr. As the number of decomposition layers increases, the center frequency of the filter gradually decreases and the bandwidth also decreases. When the sampling frequency of the signal is fs, the center frequency fc and bandwidth Bω can be expressed as equations (5) and (6), respectively:(5)fc=αj2−β4αfs,(6)Bω=12βαj−1πfs.


### 2.2. 1.5D Spectral Kurtosis

For a discrete sequence of signalxt, the system response is represented by Xt, which can be expressed as(7)Xt=∫−∞+∞e2πfjtHt,fdxf,where Ht,f is the transfer function and dxf is the spectral increment.The second-order spectral distance ofXt can be expressed as(8)S2nyt,f=EHt,fdXf2ndf=Ht,f2n⋅S2nx.When 2n = 2, equation (8) can be abbreviated as(9)S2yt,f=EHt,f2⋅σx2.Furthermore, the fourth-order spectral cumulant can be expressed as(10)C4yf=S4yf−2S2y2f.The spectral kurtosis can be defined by fourth-order normalized cumulants, as shown in(11)Kx=C4yfS2y2f=S4yf−2S2y2fS2y2f=S4yfS2y2f−2=Ht,f4Ht,f22−2.For a signalxt, the diagonal slice of the third-order cumulant R3xτ1,τ2τ1=τ2=τ can be expressed by(12)R3xτ,τ=Extxt+τxt−τ,where E⋅ represents mathematical expectation.The 1.5D spectrum of the original signalxt can be obtained by performing Fourier transform on the three diagonal slices as(13)Bω=∫−∞∞R3xτ,τe−jωτdτ.The spectral distance of the spectrum discrete sequence is calculated and brought into equation (11), and the 1.5D spectral kurtosis value K1.5D is obtained as(14)K1.5D=B3t,f4B3t,f22−2.In order to compare the robustness of the kurtosis index defined in the time domain with the 1.5D spectral kurtosis index defined in the frequency domain against noise, different degrees of noise are added to the periodic impulse simulated signal. The mathematical model of the simulated signal can be expressed by(15)St=∑i=1nAe−ξt−qit/fouter2⋅sin2πfot+nt,where A is the amplitude, ξ is the damping ratio, fo is the natural frequency of the system, fouter is the characteristic frequency, and nt is the noise. The parameters of the simulated signal are listed in Table 1.Table 1
Parameters of the simulated signal.
Parameter typeParameter valuesAmplitude A2Natural frequencyfo (Hz)2100Sampling frequency (Hz)8192Sampling points8192Fault frequency (Hz)50Figures2(a)–2(p) are the time-domain waveforms and 1.5D spectrum of simulated signals with different SNR. Figure 2(q) shows the variation trend of the kurtosis index in the time domain and the 1.5D spectral kurtosis index in the frequency domain with the SNR. It can be seen from the figure that when the periodic impulse in the time-domain signal is gradually overwhelmed by noise, the spectral impact component composed of the fault characteristic frequency and its harmonics in the 1.5D spectrum is still obvious. This shows that the 1.5D spectrum has a relatively strong ability to suppress noise. By comparing the variation trends of the time-domain kurtosis index and the 1.5D spectral kurtosis index, it can also be found that when the SNR ratio is low, the 1.5D spectral kurtosis has a higher slope, indicating that it has a higher resolution.Figure 2
Curve of the kurtosis and 1.5D spectral kurtosis with different SNR. (a)SNR=−2dB. (b) SNR=−2dB. (c) SNR=−4dB. (d) SNR=−4dB. (e) SNR=−6dB. (f) SNR=−6dB. (g) SNR=−8dB. (h) SNR=−8dB. (i) SNR=−10dB. (g) SNR=−10dB. (k) SNR=−12dB. (l) SNR=−12dB. (m) SNR=−14dB. (n) SNR=−14dB. (o) SNR=−16dB. (p) SNR=−16dB. (q) SNRdB.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)(m)(n)(o)(p)(q)Then, the following simulated signals (Figures3(a) and 3(b)) are used to compare the robustness of the kurtosis index and the 1.5D spectral kurtosis index against nonperiodic transient impulse components. The time-domain kurtosis and 1.5D spectral kurtosis are calculated for Gaussian noise, and the results are 2.9587 and 12.1610, respectively. Then, a set of nonperiodic transient impulse are added at 0.3s of the simulated signal, and the kurtosis and 1.5D spectral kurtosis values are 9.3163 and 11.9932, respectively. The results show that the kurtosis index defined in the time domain is easily interfered by nonperiodic transient impulse, and the 1.5D spectral kurtosis value is more robust to nonperiodic transient impulse.Figure 3
Analysis of simulated signals. (a) Gaussian noise; (b) noise accompanied by nonperiodic transient impulse.
(a)(b)


### 2.3. Proposed Method

(1)
Firstly, the parameters of TQWT are selected. According to the characteristics of the frequency band where the bearing fault information is located, the number of signal decomposition layersJ is selected as 6, which is beneficial to abandon the low-frequency modulation harmonic characteristics and the high-frequency interference influence. The Q-factor is obtained through traversal optimization within a certain range, and this paper uses 0.5 as the step size to obtain the optimization between 2.5 and 5.(2)
The subbands under different Q-factors and decomposition scales are demodulated by Hilbert envelope, and the 1.5D spectrum of the demodulated signal is calculated, and then the 1.5D spectral kurtosis value is calculated for each discrete spectrum sequence.(3)
The optimal subband is selected, and its time-domain signal and the corresponding 1.5D spectrum are analyzed, which can extract fault feature information.


## 2.1. TQWT

Compared with the traditional discrete wavelet transform, TQWT can determine the wavelet basis function by designing the Q-factor (the Q-factor is the ratio of signal center frequency to bandwidth,Q=fc/Bω) and the redundancy factor according to the signal oscillation characteristics. This design makes the choice of wavelet function more flexible, makes the signal adaptability stronger, and makes it easier to control the sparsity of frequency band division. TQWT uses band-pass filter banks to realize signal decomposition and reconstruction in an iterative manner. The three-level decomposition and reconstruction is taken as an example. The corresponding dual channel filter banks are shown in Figure 1.Figure 1
The corresponding dual channel filter banks of (a) decomposition process and (b) reconstruction process.
(a)(b)In Figure1, H1ω and H0ω are the frequency response functions of the high-pass filter and the low-pass filter (H1ω and H0ω can be expressed as equations (1) and (2)), H1∗ω and H0∗ω are the conjugate complex numbers of H1ω and H0ω, respectively, β=2/Q+1 is the high-pass scaling factor, and α=1−β/r is the low-pass scaling factor, where r is the redundancy factor:(1)H0ω=1,ω≤1−β,θω+β−1πα+β−1,1−βπ≤ω<απ,0,απ≤ω<π,(2)H1ω=0,ω≤1−β,θαπ−ωα+β−1,1−βπ≤ω<απ,1,απ≤ω<π.The maximum number of decomposition layers of TQWT can be obtained by(3)Jmax=logβNs/8log1/α.The scale factorsα and β can be substituted into equation (3) to obtain(4)Jmax=logNs/4Q+1logQ+1/Q+1−2/r.In summary, the filtering characteristics of TQWT change with the Q-factor and the redundancy factorr. As the number of decomposition layers increases, the center frequency of the filter gradually decreases and the bandwidth also decreases. When the sampling frequency of the signal is fs, the center frequency fc and bandwidth Bω can be expressed as equations (5) and (6), respectively:(5)fc=αj2−β4αfs,(6)Bω=12βαj−1πfs.


## 2.2. 1.5D Spectral Kurtosis

For a discrete sequence of signalxt, the system response is represented by Xt, which can be expressed as(7)Xt=∫−∞+∞e2πfjtHt,fdxf,where Ht,f is the transfer function and dxf is the spectral increment.The second-order spectral distance ofXt can be expressed as(8)S2nyt,f=EHt,fdXf2ndf=Ht,f2n⋅S2nx.When 2n = 2, equation (8) can be abbreviated as(9)S2yt,f=EHt,f2⋅σx2.Furthermore, the fourth-order spectral cumulant can be expressed as(10)C4yf=S4yf−2S2y2f.The spectral kurtosis can be defined by fourth-order normalized cumulants, as shown in(11)Kx=C4yfS2y2f=S4yf−2S2y2fS2y2f=S4yfS2y2f−2=Ht,f4Ht,f22−2.For a signalxt, the diagonal slice of the third-order cumulant R3xτ1,τ2τ1=τ2=τ can be expressed by(12)R3xτ,τ=Extxt+τxt−τ,where E⋅ represents mathematical expectation.The 1.5D spectrum of the original signalxt can be obtained by performing Fourier transform on the three diagonal slices as(13)Bω=∫−∞∞R3xτ,τe−jωτdτ.The spectral distance of the spectrum discrete sequence is calculated and brought into equation (11), and the 1.5D spectral kurtosis value K1.5D is obtained as(14)K1.5D=B3t,f4B3t,f22−2.In order to compare the robustness of the kurtosis index defined in the time domain with the 1.5D spectral kurtosis index defined in the frequency domain against noise, different degrees of noise are added to the periodic impulse simulated signal. The mathematical model of the simulated signal can be expressed by(15)St=∑i=1nAe−ξt−qit/fouter2⋅sin2πfot+nt,where A is the amplitude, ξ is the damping ratio, fo is the natural frequency of the system, fouter is the characteristic frequency, and nt is the noise. The parameters of the simulated signal are listed in Table 1.Table 1
Parameters of the simulated signal.
Parameter typeParameter valuesAmplitude A2Natural frequencyfo (Hz)2100Sampling frequency (Hz)8192Sampling points8192Fault frequency (Hz)50Figures2(a)–2(p) are the time-domain waveforms and 1.5D spectrum of simulated signals with different SNR. Figure 2(q) shows the variation trend of the kurtosis index in the time domain and the 1.5D spectral kurtosis index in the frequency domain with the SNR. It can be seen from the figure that when the periodic impulse in the time-domain signal is gradually overwhelmed by noise, the spectral impact component composed of the fault characteristic frequency and its harmonics in the 1.5D spectrum is still obvious. This shows that the 1.5D spectrum has a relatively strong ability to suppress noise. By comparing the variation trends of the time-domain kurtosis index and the 1.5D spectral kurtosis index, it can also be found that when the SNR ratio is low, the 1.5D spectral kurtosis has a higher slope, indicating that it has a higher resolution.Figure 2
Curve of the kurtosis and 1.5D spectral kurtosis with different SNR. (a)SNR=−2dB. (b) SNR=−2dB. (c) SNR=−4dB. (d) SNR=−4dB. (e) SNR=−6dB. (f) SNR=−6dB. (g) SNR=−8dB. (h) SNR=−8dB. (i) SNR=−10dB. (g) SNR=−10dB. (k) SNR=−12dB. (l) SNR=−12dB. (m) SNR=−14dB. (n) SNR=−14dB. (o) SNR=−16dB. (p) SNR=−16dB. (q) SNRdB.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)(m)(n)(o)(p)(q)Then, the following simulated signals (Figures3(a) and 3(b)) are used to compare the robustness of the kurtosis index and the 1.5D spectral kurtosis index against nonperiodic transient impulse components. The time-domain kurtosis and 1.5D spectral kurtosis are calculated for Gaussian noise, and the results are 2.9587 and 12.1610, respectively. Then, a set of nonperiodic transient impulse are added at 0.3s of the simulated signal, and the kurtosis and 1.5D spectral kurtosis values are 9.3163 and 11.9932, respectively. The results show that the kurtosis index defined in the time domain is easily interfered by nonperiodic transient impulse, and the 1.5D spectral kurtosis value is more robust to nonperiodic transient impulse.Figure 3
Analysis of simulated signals. (a) Gaussian noise; (b) noise accompanied by nonperiodic transient impulse.
(a)(b)


## 2.3. Proposed Method

(1)
Firstly, the parameters of TQWT are selected. According to the characteristics of the frequency band where the bearing fault information is located, the number of signal decomposition layersJ is selected as 6, which is beneficial to abandon the low-frequency modulation harmonic characteristics and the high-frequency interference influence. The Q-factor is obtained through traversal optimization within a certain range, and this paper uses 0.5 as the step size to obtain the optimization between 2.5 and 5.(2)
The subbands under different Q-factors and decomposition scales are demodulated by Hilbert envelope, and the 1.5D spectrum of the demodulated signal is calculated, and then the 1.5D spectral kurtosis value is calculated for each discrete spectrum sequence.(3)
The optimal subband is selected, and its time-domain signal and the corresponding 1.5D spectrum are analyzed, which can extract fault feature information.


## 3. Simulated Analysis

In order to demonstrate the process of the proposed method and verify the effectiveness of the proposed method, a group of simulated signals are used for analysis. The mathematical model of simulated signal can be expressed by(16)y1=∑i=1nAe−ξt−qit/fouter2⋅sin2πfot,y2=A1e−b1×t−t1⋅sin2πf1⋅t−t1+A2e−b2×t−t2⋅sin2πf2⋅t−t2,y3=A3⋅sin2πfmesht+β11+A4sin2πfht+β2,Y=y1+y2+y3+nt.The simulated signal consists of four parts:y1 represents periodic impulse components, y2 represents nonperiodic transient impulse components, y3 represents modulation harmonic components, nt represents a certain degree of noise, and Y is the composite signal of four components. Where A=2 is the amplitude of the periodic impulse components, A1=3 and A2=6 are the amplitudes of the nonperiodic transient impulse components, A3=A4=0.6 are the amplitudes of the modulation harmonic components. ξ is the damping ratio, fo is the natural frequency of the system, and fouter is the characteristic frequency of the periodic impulse components. b1=260, b2=360, t1=0.2, t2=0.3, and f1=f2=3600Hz are damping parameters, time parameters, and natural frequencies of nonperiodic transient impulse components, respectively. nt is the noise (SNR = −12 dB). The parameters of the simulated signal are listed in Table 2.Table 2
Parameters of the simulated signal.
Parameter typeParameter values (Hz)Sampling frequency8192Natural frequencyfo2100Natural frequencyf1=f23600Carrier frequencyfmesh360Modulation frequencyfh15Fault frequencyfouter50Figures4(a)–4(e) show the time-domain waveforms of four components in the simulated signal, and they are periodic impulse components, nonperiodic transient impulse components, modulation harmonic components, and noise components, and Figure 4(e) shows the composite signal of four components.Figure 4
Time-domain waveform of simulated signals. (a) Periodic impulse componentsy1, (b) nonperiodic transient impulse components y2, (c) modulation harmonic components y3, (d) noise components nt, and (e) composite signal of four components Y.
(a)(b)(c)(d)(e)The simulated signal is analyzed using the method proposed in this paper. By adjusting the adjustable factorQ, the TQWT with different sparsity are constructed, and then the 1.5D spectral kurtosis values corresponding to different wavelet subbands are calculated. The result is shown in Figure 5. The horizontal axis represents different scales, the color of the line represents the Q-factor, and the vertical axis represents the 1.5D spectral kurtosis value. It can be seen that when the Q-factor is 5 and the scale is 5, the corresponding subband index value reaches the maximum, indicating that the subband is the optimal subband.Figure 5
Optimal subband selection process of the 1.5D spectral kurtosis-guided TQWT method.The process of finding the optimal Q-factor and scale is equivalent to finding the optimal center frequency and bandwidth of the filter bank, which can be represented by the following equivalent filtering process (Figure6).Figure 6
Equivalent filtering characteristics of subbands with different Q-values and scales.
(a)(b)(c)(d)(e)(f)The selected subband with Q-factor of 5 and scale of 5 is analyzed, and its time domain and 1.5D spectrum are shown in Figures7(a) and 7(b). And the fault characteristic frequency and its multiplication component can be extracted.Figure 7
(a) Time-domain waveform and (b) 1.5D spectrum of the optimal subband.
(a)(b)For comparison, the kurtosis-guided TQWT method is used to select subbands. The optimal subband is corresponding to the Q-factor of 5 and scale of 1 (Figure8(a)). The subband signal is analyzed (Figures 8(b) and 8(c)), and it can be found that the subband signal is the nonperiodic transient impulse components, and no fault characteristic information can be obtained.Figure 8
(a) Kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) envelope spectrum of (b).
(a)(b)(c)


## 4. Experimental Analysis

### 4.1. Experimental Data with Nonperiodic Transient Impulse Components

In order to further verify the effectiveness of the proposed method in processing bearing fault signals containing abnormal impulse components, the bearing vibration signal data of CWRU are analyzed. The experimental platform is shown in Figure9. The bearing model used is JEMSKF6203-2RS. The selected experimental data are the rolling element fault data with relatively small fault size and weak fault degree. The driving motor speed is 1750 rpm, and the sampling frequency is 12000 Hz. Through the above information, combined with the calculation equations of bearing fault characteristic frequency, the theoretical fault characteristic frequency under different types of bearing faults can be obtained, as shown in Table 3.Figure 9
Bearing fault experimental equipment of CWRU.Table 3
Theoretical fault characteristic frequency of bearings.
Type of faultParameter values (Hz)Outer ring fault88Inner ring fault143Rolling element fault115The time-domain waveform of the intercepted experimental data is shown in Figure10. It can be seen from the figure that this segment of data has obvious nonperiodic transient impulse components. The original signal is processed by the method proposed in this paper, and the results are shown in Figure 11. It can be seen from Figure 11(a) that the optimal subband is selected when the Q-factor is 5 and the scale is 2. Then, the subband is analyzed. The time-domain waveform is shown in Figure 11(b), which contains obvious periodic impulse components. The 1.5D spectrum is shown in Figure 11(c). The characteristic frequency of rolling element fault is extracted. This shows that the proposed method is robust to the nonperiodic transient impulse components.Figure 10
Time-domain waveform of the intercepted experimental data.Figure 11
(a) Optimal subband selection process of the 1.5D spectral kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) 1.5D spectrum of (b).
(a)(b)(c)As a comparison, the kurtosis index is used to select the subband signals under different Q-factors and scales. The optimal subband is shown in Figure12(a). The time-domain waveform of the subband is shown in Figure 12(b). It can be seen that its main component is nonperiodic transient impulse, and its 1.5D spectrum is shown in Figure 12(c), and no obvious fault feature information can be found.Figure 12
(a) Optimal subband selection process of the kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) envelope spectrum of (b).
(a)(b)(c)


### 4.2. Experimental Data with Low Signal-to-Noise Ratio

In order to verify the processing ability of the proposed method for low signal-to-noise ratio signal, the data of QPZZ-II bearing test-bed are used for analysis and verification. The structure of the experimental platform is shown in Figure13. The power source of the test-bed is the variable frequency drive motor. The motor is connected with the rotor through the coupling. The two ends of the rotor are bearing seats. The normal bearing is in the left bearing seat, and the fault bearing is in the right bearing seat. The structural parameters of the bearing are shown in Table 4. The bearing model is 6205E. The speed of driving motor is 1470 rpm. The theoretical fault characteristic frequency under different types of bearing faults can be obtained, as shown in Table 5.Figure 13
(a) Structure of the experimental platform; (b) bearing model of 6205E.
(a)(b)Table 4
Structural parameters of the bearing.
Parameter typeParameter valuesInside diameter (mm)25Outside diameter (mm)52Ball diameter (mm)7.9Pitch diameter (mm)39Number of balls9Contact angle (°)0Table 5
Theoretical fault characteristic frequency of bearings.
Type of faultParameter values (Hz)Outer ring fault87.8Inner ring fault132.6The time-domain waveform of the experimental signal is shown in Figure14. It can be seen from the time-domain waveform that the signal-to-noise ratio of the vibration signal is low, and the periodic impulse characteristics are submerged. The original signal is analyzed by the kurtogram algorithm, and the analysis results are shown in Figure 15. It can be seen from Figure 15(a) that the optimal frequency band is in the low-frequency band (the frequency band is actually noise band, and the frequency band where the periodic impulse characteristics of bearing fault is located is about 2000 Hz, which is submerged). The time-domain waveform and envelope spectrum of the optimal subband signal are shown in Figures 15(b) and 15(c), respectively. The method fails to extract the effective fault information.Figure 14
Time-domain waveform of the experimental data.Figure 15
(a) Optimal subband selection process of the kurtogram method, (b) time-domain waveform of the optimal subband, and (c) envelope spectrum of (b).
(a)(b)(c)The original signal is processed by the method proposed in this paper, and the results are shown in Figure16. It can be seen from Figure 16(a) that the optimal subband is selected when the Q-factor is 3 and the scale is 5. Then, the subband is analyzed. The time-domain waveform is shown in Figure 16(b). The 1.5D spectrum is shown in Figure 16(c). The characteristic frequency of rolling element fault is extracted. This shows that the proposed method is robust to low signal-to-noise ratio signal.Figure 16
(a) Optimal subband selection process of the 1.5D spectral kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) 1.5D spectrum of (b).
(a)(b)(c)


## 4.1. Experimental Data with Nonperiodic Transient Impulse Components

In order to further verify the effectiveness of the proposed method in processing bearing fault signals containing abnormal impulse components, the bearing vibration signal data of CWRU are analyzed. The experimental platform is shown in Figure9. The bearing model used is JEMSKF6203-2RS. The selected experimental data are the rolling element fault data with relatively small fault size and weak fault degree. The driving motor speed is 1750 rpm, and the sampling frequency is 12000 Hz. Through the above information, combined with the calculation equations of bearing fault characteristic frequency, the theoretical fault characteristic frequency under different types of bearing faults can be obtained, as shown in Table 3.Figure 9
Bearing fault experimental equipment of CWRU.Table 3
Theoretical fault characteristic frequency of bearings.
Type of faultParameter values (Hz)Outer ring fault88Inner ring fault143Rolling element fault115The time-domain waveform of the intercepted experimental data is shown in Figure10. It can be seen from the figure that this segment of data has obvious nonperiodic transient impulse components. The original signal is processed by the method proposed in this paper, and the results are shown in Figure 11. It can be seen from Figure 11(a) that the optimal subband is selected when the Q-factor is 5 and the scale is 2. Then, the subband is analyzed. The time-domain waveform is shown in Figure 11(b), which contains obvious periodic impulse components. The 1.5D spectrum is shown in Figure 11(c). The characteristic frequency of rolling element fault is extracted. This shows that the proposed method is robust to the nonperiodic transient impulse components.Figure 10
Time-domain waveform of the intercepted experimental data.Figure 11
(a) Optimal subband selection process of the 1.5D spectral kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) 1.5D spectrum of (b).
(a)(b)(c)As a comparison, the kurtosis index is used to select the subband signals under different Q-factors and scales. The optimal subband is shown in Figure12(a). The time-domain waveform of the subband is shown in Figure 12(b). It can be seen that its main component is nonperiodic transient impulse, and its 1.5D spectrum is shown in Figure 12(c), and no obvious fault feature information can be found.Figure 12
(a) Optimal subband selection process of the kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) envelope spectrum of (b).
(a)(b)(c)


## 4.2. Experimental Data with Low Signal-to-Noise Ratio

In order to verify the processing ability of the proposed method for low signal-to-noise ratio signal, the data of QPZZ-II bearing test-bed are used for analysis and verification. The structure of the experimental platform is shown in Figure13. The power source of the test-bed is the variable frequency drive motor. The motor is connected with the rotor through the coupling. The two ends of the rotor are bearing seats. The normal bearing is in the left bearing seat, and the fault bearing is in the right bearing seat. The structural parameters of the bearing are shown in Table 4. The bearing model is 6205E. The speed of driving motor is 1470 rpm. The theoretical fault characteristic frequency under different types of bearing faults can be obtained, as shown in Table 5.Figure 13
(a) Structure of the experimental platform; (b) bearing model of 6205E.
(a)(b)Table 4
Structural parameters of the bearing.
Parameter typeParameter valuesInside diameter (mm)25Outside diameter (mm)52Ball diameter (mm)7.9Pitch diameter (mm)39Number of balls9Contact angle (°)0Table 5
Theoretical fault characteristic frequency of bearings.
Type of faultParameter values (Hz)Outer ring fault87.8Inner ring fault132.6The time-domain waveform of the experimental signal is shown in Figure14. It can be seen from the time-domain waveform that the signal-to-noise ratio of the vibration signal is low, and the periodic impulse characteristics are submerged. The original signal is analyzed by the kurtogram algorithm, and the analysis results are shown in Figure 15. It can be seen from Figure 15(a) that the optimal frequency band is in the low-frequency band (the frequency band is actually noise band, and the frequency band where the periodic impulse characteristics of bearing fault is located is about 2000 Hz, which is submerged). The time-domain waveform and envelope spectrum of the optimal subband signal are shown in Figures 15(b) and 15(c), respectively. The method fails to extract the effective fault information.Figure 14
Time-domain waveform of the experimental data.Figure 15
(a) Optimal subband selection process of the kurtogram method, (b) time-domain waveform of the optimal subband, and (c) envelope spectrum of (b).
(a)(b)(c)The original signal is processed by the method proposed in this paper, and the results are shown in Figure16. It can be seen from Figure 16(a) that the optimal subband is selected when the Q-factor is 3 and the scale is 5. Then, the subband is analyzed. The time-domain waveform is shown in Figure 16(b). The 1.5D spectrum is shown in Figure 16(c). The characteristic frequency of rolling element fault is extracted. This shows that the proposed method is robust to low signal-to-noise ratio signal.Figure 16
(a) Optimal subband selection process of the 1.5D spectral kurtosis-guided TQWT method, (b) time-domain waveform of the optimal subband, and (c) 1.5D spectrum of (b).
(a)(b)(c)


## 5. Conclusions

The core problem in bearing signal processing is how to correctly select the frequency band of the fault information, so as to extract the information within the frequency band and filter out the interference information outside the frequency band. TQWT can provide a signal frequency band division method, and because its Q-factor is adjustable, it can adjust the wavelet oscillation characteristics to adapt to various signal characteristics. In order to improve the accuracy of subband selection and avoid the noise and nonperiodic transient impulse interference in the environment, this paper proposes a Q-factor and scale optimization method of TQWT with 1.5D spectral kurtosis as the evaluation index. Simulated and experimental results show that the proposed method has good analysis effect when SNR is low and when there is nonperiodic transient impulse interference.



---

*Source: 5554981-2021-05-17.xml*