# Rolling Bearing Fault Signal Extraction Based on Stochastic Resonance-Based Denoising and VMD

Authors: Xiaojiao Gu; Changzheng Chen
Journal: International Journal of Rotating Machinery
 (2017)

---

## Abstract

Aiming at the difficulty of early fault vibration signal extraction of rolling bearing, a method of fault weak signal extraction based on variational mode decomposition (VMD) and quantum particle swarm optimization adaptive stochastic resonance (QPSO-SR) for denoising is proposed. Firstly, stochastic resonance parameters are optimized adaptively by using quantum particle swarm optimization algorithm according to the characteristics of the original fault vibration signal. The best stochastic resonance system parameters are output when the signal to noise ratio reaches the maximum value. Secondly, the original signal is processed by optimal stochastic resonance system for denoising. The influence of the noise interference and the impact component on the results is weakened. The amplitude of the fault signal is enhanced. Then the VMD method is used to decompose the denoised signal to realize the extraction of fault weak signals. The proposed method was applied in simulated fault signals and actual fault signals. The results show that the proposed method can reduce the effect of noise and improve the computational accuracy of VMD in noise background. It makes VMD more effective in the field of fault diagnosis. The proposed method is helpful to realize the accurate diagnosis of rolling bearing early fault.


---

## Body

## 1. Introduction

Rolling bearing is one of the most critical and easily damaged components in rotating machinery. Its running state is directly related to the performance of the whole mechanical system [1–3]. Therefore, it is significant to implement fault diagnosis for rolling bearings so as to prevent fatal malfunction of rotating machinery [4]. Vibration analysis is one of the most frequently used methods for health monitoring and fault diagnosis of rotating machinery. The vibration signal of mechanical equipment is rich in status information [5, 6]. When a rotating machine has incipient faults or works in a wicked environment, the useful fault signals are often submerged in strong background noise. It will seriously affect the accuracy of fault diagnosis. Therefore, the extraction of weak fault signal is a research hotspot in the field of signal processing [7].In order to extract useful signals from complex original signals, many scholars have proposed a lot of effective weak signal extraction methods. One category is to extract weak signal from the perspective of denoising, such as wavelet transform, chaos theory, empirical mode decomposition (EMD), local mean decomposition (LMD), ensemble empirical mode decomposition (EEMD), and singular value decomposition (SVD). Under the influence of sampling frequency, the EEMD decomposition error is large. Meng and Xiang [8] proposed an improved EEMD which extracts the real intrinsic mode functions (IMF) by using the correlation coefficients between the original signal and IMFs. With this method, the pseudo low-frequency IMFs can be eliminated. Yi et al. [9] proposed the augmented quaternion singular spectrum analysis multichannel denoising method. This method has a better ability than multivariate EMD method in multisignal processing. Xie et al. [10] proposed an improved LMD based on extension characteristic wave method to eliminate the end effect. Yi et al. [11] proposed a convex optimization algorithm using nonconvex penalty functions based on SVD for extracting weak fault characteristics. Variational mode decomposition (VMD) is a new adaptive signal processing method proposed by Dragomiretskiy and Zosso in 2014 [12]. In this method, the frequency center and bandwidth of each component are determined by iteration searching for the optimal solution of the variational mode during the process of obtaining decomposition components. Therefore, the effective separation of the frequency domain can be realized adaptively. The existing study shows that in complicated signal decomposition the VMD has more advantage than EMD [13]. Compared with the recursive filtering mode of EEMD and LMD, VMD is a nonrecursive and variational mode decomposition. It has a solid theoretical basis. Its essence is a number of adaptive Wiener filter groups, which show better noise robustness. Through the reasonable control for the convergence condition, the influence of sampling on VMD is much smaller than EEMD and LMD. In terms of mode separation, the VMD can successfully separate two pure harmonic signals with similar frequencies. At present, the research results of the variational mode decomposition in the field of mechanical fault diagnosis show that VMD has better performance than LMD and EEMD in many aspects [14]. The fault feature extraction based on variational mode decomposition has important application value.The other category is to use the noise to enhance the fault signal by means of stochastic resonance (SR). Because of its obvious advantages in weak signal detection, SR has become a hotspot in the field of signal processing in recent years [15, 16]. Han et al. [17] combined wavelet transform with stochastic resonance theory, which can detect several high frequency weak signals in strong noise background by adjusting the amplitude of wavelet multiscale decomposition. Rother et al. [18] combined stochastic resonance with EMD, which is applied to detect deviations in strip travel of hot strip mill. The results show that this method can reduce downtimes and extend the lifetime of components. Barbini et al. [19] focused on the use of SR in a two-dimensional system of gradient type for detection of weak signals submerged in Gaussian noise. The two-dimensional system is more accurate and more reliable than the traditional one-dimensional system. Castiglione et al. [20] used bistable stochastic resonance and Woods-Saxon stochastic resonance in bearing faults detection, respectively. Experimental data shows that each SR method has its own characteristic features. Gao et al. [21] constructed corresponding Duffing oscillator based on chaos theory to monitor weak signal, which reduces the SNR threshold compared with the traditional monitoring method. Lu et al. [22] enhanced stochastic resonance by full-wave signal construction, which surpass the traditional stochastic resonance in the field of fault signal processing and can be used in areas related to weak signal detection. Lin et al. [23] presented a new additive noise driving and aperiodic chirp signal stochastic dynamical system, which can improve the system performance of localization and moving target detection. Breen et al. [24] used the mechanical stochastic resonance to detect weak signals, which harvest weak periodic signals via the movement of an inverted pendulum. Shi et al. [25] combined analytical mode decomposition (AMD) and EEMD with adaptive stochastic resonance, which can improve the decomposition results of EEMD. This method has a good effect in the rolling bearing fault diagnosis. Xiang and Zhong [26] presented a two-step method of Duffing oscillator and stochastic resonance. And this method combined stochastic resonance with the artificial bee colony algorithm. It is favorable to detect weak signal from a strong noise environment. Liu et al. [27] combined wavelet reconstruction and stochastic resonance. The reconstruction scale was selected by comparing the wavelet entropy of each decomposition scale. This fault detection method is more accurate than the conventional wavelet reconstruction. The above studies which combine the stochastic resonance with the signal processing method or the optimization method show that the appropriate application of the stochastic resonance method can make the signal processing method have a better effect in fault diagnosis. Stochastic resonance has important application value in fault diagnosis under strong noise background.Combining the respective advantages of VMD and SR, this study presents a weak signal extraction method of rolling bearing fault based on VMD and quantum particle swarm optimization (QPSO) adaptive stochastic resonance. First, the parameters of stochastic resonance system are optimized according to the original signal feature and quantum particle swarm optimization algorithm. Thereafter, the original signal is processed by the stochastic resonance so that the noise frequency component is weakened and the fault frequency component is enhanced. Lastly, the signal which is processed by stochastic resonance is decomposed by VMD method. The useful signal features are extracted. The rolling bearing measured signal analysis shows that the proposed method can discharge most of the interfering signal generated by the background noise so that the components obtained by the VMD can better reflect the fault signal information and enhance the useful signal amplitude. The detection results are more accurate and reliable.The rest of this paper is organized as follows. Section2 provides a brief description of the principal theory of VMD, SR, and QPSO. Section 3 introduces the proposed fault diagnosis process and discusses the simulation study to validate the performance of the proposed method. Section 4 discusses the case study to validate the practical application value of this method. Lastly, Section 5 presents the conclusion.


## 2. Basic Theory

### 2.1. Variational Mode Decomposition

VMD algorithm is a new nonrecursive variational mode signal decomposition method, which can decompose the complex input signal into a set of discrete mode components [28]. The implementation steps of the VMD algorithm are as follows.Step 1.
The Hilbert transform is performed for each mode functionuk to obtain its analytical signal, as follows:(1)δt+jπt∗ukt.Step 2.
The estimated center frequencyejωkt of uk is mixed. The spectrum of each mode is modulated to the corresponding base band, as follows:(2)δt+jπt∗ukt∗e-jωkt.Step 3.
Calculate the gradient squareL2 norm of the above demodulated signal and estimate the bandwidth of each mode signal. The variational problem is constructed as follows: (3)minuk,ωk∑k∂tδt+jπt∗ukt∗e-jωkt2,s.t.∑kuk=f,where {uk}≔{u1,…,uk} represents K variational mode components. {ωk}≔{ω1,…,ωk} represents the center frequency of each component. K is the number of variational mode components.t represents the time. ∂t represents the partial derivative of t. δt represents the impulse function. j is the square root of −1. ω represents the circular frequency. f represents the input signal.Step 4.
In order to transform the constraint variational problem into the unconstrained problem, the quadratic penalty factorα is introduced to ensure the accuracy of signal reconstruction in noise condition. The Lagrangian multiplier λ is introduced to ensure the stringency of the constraint condition. Through introducing the two parameters, the extended Lagrangian expression is obtained as follows:(4)Luk,ωk,λ≔α∑k∂tδt+jπt∗ukt∗e-jωkt22+ft-∑kukt22+λt,ft-∑kukt.Step 5.
Solve the extended Lagrangian expression by using the alternate direction method of multipliers (ADMM). Specific steps are as follows:(1) Initialize {u^k1}, {ωk1}, λ^1, n.(2) Repeat n=n+1.(3) Update uk, ωk, and λ, and do(5)u^kn+1ω⟵f^ω-∑i≠ku^iω+λ^ω/21+2αω-ωk2,ωkn+1⟵∫0∞ωu^kω2dω∫0∞u^kω2dω,λ^n+1ω⟵λ^nω+τf^ω-∑ku^kn+1ω.(4) Judge whether or not the following conditions are met. If met, the iteration is stopped; otherwise, return to execute Step (2).(6)∑ku^kn+1-u^kn22u^kn22<e,where e is the convergence stopped condition. k is the iteration number.


### 2.2. Stochastic Resonance

Benzi et al. [29] first introduced the mechanism of stochastic resonance in 1981 and used SR in climatic change. In the field of signal processing, stochastic resonance can transfer the energy of some noise signals to useful signals, which is beneficial to the identification of weak signals submerged in noise [30–32]. The noise signal is input to the bistable or multistable nonlinear system. The system performance can be optimized by adjusting the parameters so that the effect of stochastic resonance achieves the best. The equation of the bistable Langevin system is as follows:(7)x˙=-dUxdx+st+nt,where U(x) is the potential function.(8)Ux=-12ax2+14bx4,where a and b are system parameters. The potential function (8) has one unstable equilibrium point and two stable equilibrium points. The height of the potential barrier is a2/(4b). st=∑i=1nAisin⁡2πfit represents an input signal with amplitude Ai and frequency fi. nt=2Dεt represents the Gaussian white noise signal with intensity D, mean 0, and variance 1. The bistable Langevin system structure is shown in Figure 1.Figure 1
SR structure for the bistable Langevin system.The stochastic resonance problem is equivalent to a multidimensional-continuous optimization problem [33]. The same set of input signals using different parameters of the system for stochastic resonance processing will produce different effect [34]. Setting st=Asin⁡2πf0t, Figure 2 shows the response amplitude x- changes with SR system parameters when D=0.2, f0=0.1 Hz, and A=0.1. It is shown that there exist optimal system parameters a, b to make the stochastic resonance performance optimal.Figure 2
Influence of SR system parameters on amplitude.


### 2.3. Quantum Particle Swarm Optimization

Quantum Particle Swarm Optimization was proposed on trajectory analysis of particle swarm optimization (PSO) and the quantum mechanism [35]. In the quantum space, the position and velocity of the particle cannot be determined simultaneously, so the wave function is used to describe the particle state [36]. The square of the wave function represents the probability of finding a particle at a particular position in the quantum space. By employing the Monte Carlo method, the position of a particle can be expressed as follows:(9)xt=P±L2ln⁡1υ,where υ is a random number between 0 and 1. L is the characteristic length of potential well which changes with the time as follows:(10)Lt+1=2β∑i=1NPiN-xt,where parameter β is the contraction expansion coefficient and is used to control the convergence speed of the particle. N is the number of particle swarms. Pi is the local optimal value of particle i. Finally, the position of the particle can be expressed as follows:(11)xt+1=P±β∑i=1NPiN-xtln⁡1υ.


## 2.1. Variational Mode Decomposition

VMD algorithm is a new nonrecursive variational mode signal decomposition method, which can decompose the complex input signal into a set of discrete mode components [28]. The implementation steps of the VMD algorithm are as follows.Step 1.
The Hilbert transform is performed for each mode functionuk to obtain its analytical signal, as follows:(1)δt+jπt∗ukt.Step 2.
The estimated center frequencyejωkt of uk is mixed. The spectrum of each mode is modulated to the corresponding base band, as follows:(2)δt+jπt∗ukt∗e-jωkt.Step 3.
Calculate the gradient squareL2 norm of the above demodulated signal and estimate the bandwidth of each mode signal. The variational problem is constructed as follows: (3)minuk,ωk∑k∂tδt+jπt∗ukt∗e-jωkt2,s.t.∑kuk=f,where {uk}≔{u1,…,uk} represents K variational mode components. {ωk}≔{ω1,…,ωk} represents the center frequency of each component. K is the number of variational mode components.t represents the time. ∂t represents the partial derivative of t. δt represents the impulse function. j is the square root of −1. ω represents the circular frequency. f represents the input signal.Step 4.
In order to transform the constraint variational problem into the unconstrained problem, the quadratic penalty factorα is introduced to ensure the accuracy of signal reconstruction in noise condition. The Lagrangian multiplier λ is introduced to ensure the stringency of the constraint condition. Through introducing the two parameters, the extended Lagrangian expression is obtained as follows:(4)Luk,ωk,λ≔α∑k∂tδt+jπt∗ukt∗e-jωkt22+ft-∑kukt22+λt,ft-∑kukt.Step 5.
Solve the extended Lagrangian expression by using the alternate direction method of multipliers (ADMM). Specific steps are as follows:(1) Initialize {u^k1}, {ωk1}, λ^1, n.(2) Repeat n=n+1.(3) Update uk, ωk, and λ, and do(5)u^kn+1ω⟵f^ω-∑i≠ku^iω+λ^ω/21+2αω-ωk2,ωkn+1⟵∫0∞ωu^kω2dω∫0∞u^kω2dω,λ^n+1ω⟵λ^nω+τf^ω-∑ku^kn+1ω.(4) Judge whether or not the following conditions are met. If met, the iteration is stopped; otherwise, return to execute Step (2).(6)∑ku^kn+1-u^kn22u^kn22<e,where e is the convergence stopped condition. k is the iteration number.


## 2.2. Stochastic Resonance

Benzi et al. [29] first introduced the mechanism of stochastic resonance in 1981 and used SR in climatic change. In the field of signal processing, stochastic resonance can transfer the energy of some noise signals to useful signals, which is beneficial to the identification of weak signals submerged in noise [30–32]. The noise signal is input to the bistable or multistable nonlinear system. The system performance can be optimized by adjusting the parameters so that the effect of stochastic resonance achieves the best. The equation of the bistable Langevin system is as follows:(7)x˙=-dUxdx+st+nt,where U(x) is the potential function.(8)Ux=-12ax2+14bx4,where a and b are system parameters. The potential function (8) has one unstable equilibrium point and two stable equilibrium points. The height of the potential barrier is a2/(4b). st=∑i=1nAisin⁡2πfit represents an input signal with amplitude Ai and frequency fi. nt=2Dεt represents the Gaussian white noise signal with intensity D, mean 0, and variance 1. The bistable Langevin system structure is shown in Figure 1.Figure 1
SR structure for the bistable Langevin system.The stochastic resonance problem is equivalent to a multidimensional-continuous optimization problem [33]. The same set of input signals using different parameters of the system for stochastic resonance processing will produce different effect [34]. Setting st=Asin⁡2πf0t, Figure 2 shows the response amplitude x- changes with SR system parameters when D=0.2, f0=0.1 Hz, and A=0.1. It is shown that there exist optimal system parameters a, b to make the stochastic resonance performance optimal.Figure 2
Influence of SR system parameters on amplitude.


## 2.3. Quantum Particle Swarm Optimization

Quantum Particle Swarm Optimization was proposed on trajectory analysis of particle swarm optimization (PSO) and the quantum mechanism [35]. In the quantum space, the position and velocity of the particle cannot be determined simultaneously, so the wave function is used to describe the particle state [36]. The square of the wave function represents the probability of finding a particle at a particular position in the quantum space. By employing the Monte Carlo method, the position of a particle can be expressed as follows:(9)xt=P±L2ln⁡1υ,where υ is a random number between 0 and 1. L is the characteristic length of potential well which changes with the time as follows:(10)Lt+1=2β∑i=1NPiN-xt,where parameter β is the contraction expansion coefficient and is used to control the convergence speed of the particle. N is the number of particle swarms. Pi is the local optimal value of particle i. Finally, the position of the particle can be expressed as follows:(11)xt+1=P±β∑i=1NPiN-xtln⁡1υ.


## 3. Proposed Fault Identification Scheme

### 3.1. Fault Diagnosis Process

Different from the traditional fault identification method, this fault diagnosis process considers the noise feature of different input signal. In various moments, the fault features are varying degrees submerged in noise. Hence, fixed SR system parameters may be unable to suppress the noise and properly enhance the fault feature. This study proposes an adaptable and varying-scale stochastic resonance to overcome the shortcomings of the traditional method. The values of SR system parameters are optimized by using the quantum particle swarm optimization algorithm. The fitness of each particle is determined by the objective function. The performance of stochastic resonance system is mainly reflected in the signal to noise ratio (SNR) of the output signal [37]. Hence, the fitness function of QPSO algorithm is shown as follows:(12)fitnessx=10log⁡PSPN,where PS≔2X[k0]2 is the active power of the signal. PN≔∑k=0L-1Xk2-PS is the active power of the noise. Xk is the discrete Fourier transform (DFT) of the sample sequence, which is defined according to the following formula:(13)Xk=∑t=0L-1xtexp⁡-j2πLkt,k=0,1,…,L-1.The procedures for this method of fault feature extraction are shown in Figure 3.Figure 3
Flowchart of the proposed scheme.


### 3.2. Simulation Signals Analysis

To demonstrate the differences between VMD and the improved VMD based on QPSO-SR, the synthesized signal which includes 5 Hz and 50 Hz sinusoidal signal and the Gauss white noise is used to simulate the fault signal. The analog signal expression is established as follows:(14)xt=0.15sin⁡2π×5t+0.15sin⁡2π×50t+nt,where t ≔ 0.001 : 0.001 : 1 is the sample time points. n(t) is the Gauss white noise. The intensity of the Gauss white noise is 0.8. For comparison, the simulated signal is first subjected to VMD decomposition. The quadratic penalty factor α selected the default value 2000. The bandwidth τ selected the default value 0. The number of variational mode components is determined by the method of observation center frequency [38]. Set the initial value of K to 2. The VMD is carried out with the number of variational mode components K=K+1 and judged whether the center frequency of each mode component is close. When K=6 began to show the phenomenon of overdecomposition, the optimal mode number K is 5. The variational mode components can be obtained through VMD algorithm. The VMD results and the Fast Fourier Transform (FFT) are shown in Figure 4.Figure 4
VMD analysis result of the simulated signal.From Figure4, it can be seen that the fault information cannot be identified from time spectrum and frequency spectrum of the simulated signal. After the VMD algorithm, the amplitude of 5 Hz (represented by point 1) in u1 FFT is 0.1416 and the amplitude of 50 Hz (represented by point 2) in u2 FFT is 0.1439, showing the simulated fault information. And u3, u4, and u5 show the noise signal information. However, u1, u2 mixed with the interference noise cannot accurately identify the simulated fault features.In order to eliminate the effect of interference noise, the QPSO-SR method is used to denoise the simulated signal before VMD. The stochastic resonance system parameters a and b are optimized by quantum particle swarm optimization algorithm. The optimization range of system parameters is from 0.01 to 2. The optimization results show that SNR achieves the maximum value whena=0.05 and b=0.17. These system parameters are used to denoise the simulated signal. And then the signal processed by stochastic resonance is decomposed by VMD method. After QPSO-SR noise reduction, the result of VMD decomposition is shown in Figure 5.Figure 5
VMD analysis result of simulated signals based on QPSO-SR denoising.From Figure5, it can be seen that the interference noise in u1 and u2 is weakened. The amplitude of 5 Hz (represented by point 1) in u1 FFT is increased from 0.1416 to 0.1820. The amplitude of 50 Hz (represented by point 2) in u2 FFT is increased from 0.1439 to 0.1922. The amplitude of the fault information frequency is obviously enhanced and the useful information is accurately identified. The amplitude of the mode components u3, u4, and u5 which represent the noise signal is weakened. After the noise reduction by QPSO-SR, a part of the noise energy can be transmitted to the useful signal, which is helpful for VMD to obtain the useful information more accurately.


## 3.1. Fault Diagnosis Process

Different from the traditional fault identification method, this fault diagnosis process considers the noise feature of different input signal. In various moments, the fault features are varying degrees submerged in noise. Hence, fixed SR system parameters may be unable to suppress the noise and properly enhance the fault feature. This study proposes an adaptable and varying-scale stochastic resonance to overcome the shortcomings of the traditional method. The values of SR system parameters are optimized by using the quantum particle swarm optimization algorithm. The fitness of each particle is determined by the objective function. The performance of stochastic resonance system is mainly reflected in the signal to noise ratio (SNR) of the output signal [37]. Hence, the fitness function of QPSO algorithm is shown as follows:(12)fitnessx=10log⁡PSPN,where PS≔2X[k0]2 is the active power of the signal. PN≔∑k=0L-1Xk2-PS is the active power of the noise. Xk is the discrete Fourier transform (DFT) of the sample sequence, which is defined according to the following formula:(13)Xk=∑t=0L-1xtexp⁡-j2πLkt,k=0,1,…,L-1.The procedures for this method of fault feature extraction are shown in Figure 3.Figure 3
Flowchart of the proposed scheme.


## 3.2. Simulation Signals Analysis

To demonstrate the differences between VMD and the improved VMD based on QPSO-SR, the synthesized signal which includes 5 Hz and 50 Hz sinusoidal signal and the Gauss white noise is used to simulate the fault signal. The analog signal expression is established as follows:(14)xt=0.15sin⁡2π×5t+0.15sin⁡2π×50t+nt,where t ≔ 0.001 : 0.001 : 1 is the sample time points. n(t) is the Gauss white noise. The intensity of the Gauss white noise is 0.8. For comparison, the simulated signal is first subjected to VMD decomposition. The quadratic penalty factor α selected the default value 2000. The bandwidth τ selected the default value 0. The number of variational mode components is determined by the method of observation center frequency [38]. Set the initial value of K to 2. The VMD is carried out with the number of variational mode components K=K+1 and judged whether the center frequency of each mode component is close. When K=6 began to show the phenomenon of overdecomposition, the optimal mode number K is 5. The variational mode components can be obtained through VMD algorithm. The VMD results and the Fast Fourier Transform (FFT) are shown in Figure 4.Figure 4
VMD analysis result of the simulated signal.From Figure4, it can be seen that the fault information cannot be identified from time spectrum and frequency spectrum of the simulated signal. After the VMD algorithm, the amplitude of 5 Hz (represented by point 1) in u1 FFT is 0.1416 and the amplitude of 50 Hz (represented by point 2) in u2 FFT is 0.1439, showing the simulated fault information. And u3, u4, and u5 show the noise signal information. However, u1, u2 mixed with the interference noise cannot accurately identify the simulated fault features.In order to eliminate the effect of interference noise, the QPSO-SR method is used to denoise the simulated signal before VMD. The stochastic resonance system parameters a and b are optimized by quantum particle swarm optimization algorithm. The optimization range of system parameters is from 0.01 to 2. The optimization results show that SNR achieves the maximum value whena=0.05 and b=0.17. These system parameters are used to denoise the simulated signal. And then the signal processed by stochastic resonance is decomposed by VMD method. After QPSO-SR noise reduction, the result of VMD decomposition is shown in Figure 5.Figure 5
VMD analysis result of simulated signals based on QPSO-SR denoising.From Figure5, it can be seen that the interference noise in u1 and u2 is weakened. The amplitude of 5 Hz (represented by point 1) in u1 FFT is increased from 0.1416 to 0.1820. The amplitude of 50 Hz (represented by point 2) in u2 FFT is increased from 0.1439 to 0.1922. The amplitude of the fault information frequency is obviously enhanced and the useful information is accurately identified. The amplitude of the mode components u3, u4, and u5 which represent the noise signal is weakened. After the noise reduction by QPSO-SR, a part of the noise energy can be transmitted to the useful signal, which is helpful for VMD to obtain the useful information more accurately.


## 4. Experimental Validations

The bearing vibration signals are supplied by the Shenyang University of Technology. The rolling bearing vibration signals were measured by QPZZ-II rotating machinery vibration fault experimental platform. The experimental platform is shown in Figure6. The motor rotational speed is at 1720 r/min to 1797 r/min. The digital data used in this analysis were collected at 16,000 samples per second and the sampling number is 5468. A 0.2 mm wide and 0.1 mm deep groove was wire-cut on the surface of the bearing inner ring to simulate early fault. The detailed information of the tested rolling bearing is presented in Table 1. The fault feature frequency of the tested rolling bearing is presented in Table 2.Table 1
Geometry information of the tested rolling bearing.
Inside diameter(mm)
Outside diameter(mm)
Bearing thickness(mm)
Ball diameter(mm)
Pitch diameter(mm)
25.00
52.00
15.00
7.94
39.04Table 2
Fault feature frequency of the tested rolling bearing.
Inner ring(Hz)
Outer ring(Hz)
Rolling element(Hz)
Cage element(Hz)
5.4152fr
3.5848fr
4.7135fr
0.39828fr
f
r
=
29.95 Hz.Figure 6
Experimental facility.The time domain spectrum and frequency domain spectrum of rolling bearing vibration signal with inner ring fault are shown in Figure7. The SNR of raw vibration signal is −15.3097 dB. It can be seen from the waveform that there are obvious interference noise and shock in the rolling bearing vibration signal and the frequency component is rich. Therefore, the fault information cannot be obtained from the spectrum directly.Figure 7
Raw vibration signal.
(a)
(b)The raw vibration signal is decomposed by VMD. The decomposition which is obtained by VMD is processed by using Envelope analysis and fast Fourier transform algorithm. The VMD results and the FFT are shown in Figure8. In contrast, replace VMD with EEMD and perform the same processing on the raw signal. The EEMD results and the FFT are shown in Figure 9. The EEMD decomposition result reveals that the phenomena of mode mixing make the early fault indistinguishable. The VMD decomposition result reveals that the VMD method solves the phenomena of mode mixing but the component still contains noise disturbance and shock. VMD algorithm has not realized the potential that was expected of it. Characteristic frequency 162 Hz (represented by point 1), double frequency 324 Hz (represented by point 2), and third frequency 486 Hz (represented by point 3) of bearing inner ring fault are not obvious in the decomposition result (see Figure 8). This shows that it is difficult to ensure the accuracy of the fault diagnosis in the practical engineering if we use the VMD algorithm independently.Figure 8
VMD analysis result of the rolling bearing vibration signal.Figure 9
EEMD analysis result of the rolling bearing vibration signal.The traditional stochastic resonance method is employed to process the raw vibration signal. SR system parameters are selected according to experience. The parametera takes on 0.01, 0.03, 0.05, 0.10, and 1.19, respectively. The parameter b takes on 0.07, 0.06, 0.16, 0.12, and 1.97, respectively. The raw vibration signal is subjected to SR processing by the above five sets of parameters, respectively. The SNR of the output signal is shown in Table 3.Table 3
The SNR of the output signal.
a
b
S
N
R
o
u
t (dB)
0.01
0.07
3.8942
0.03
0.06
1.8319
0.05
0.16
3.0896
0.10
0.12
2.5896
1.19
1.97
5.0046The proposed fault weak signal extraction method based on VMD and QPSO-SR denoising is employed to process the raw vibration signal. Firstly, the raw signal is preprocessed with QPSO-SR method. According to the raw signal characteristics, the quantum particle swarm optimization algorithm is used to optimize the SR system parameters. The signal is subjected to SR processing based on the optimal parameters. The SNR of the output signal is 7.6465 dB. Evidently, the processed signals have been well denoised with the aid of QPSO-SR method. The frequency domain information exchange method is used to overcome the limitation of SR to large parameter signal. Liu et al. [39] described the specific calculation process of the frequency domain information exchange method in detail. The signal is subjected to VMD after denoising. Center frequency corresponding to different K is shown in Table 4. When K is 6, overdecomposition began to happen, so K was chosen as 5. The time domain and frequency domain waveforms of the results are shown in Figure 10.Table 4
Center frequency corresponding to differentK.
K
Center frequency (Hz)
2
1767
4784
—
—
—
—
3
1767
2789
4784
—
—
—
4
824
1767
2789
4784
—
—
5
824
1767
2789
3641
4784
—
6
824
1767
2789
3641
4668
4784Figure 10
VMD analysis result of the rolling bearing vibration signal based on QPSO-SR denoising.From Figure10, it can be seen that the VMD results are less affected by interference noise. Characteristic frequency (represented by point 1), double frequency (represented by point 2), and third frequency (represented by point 3) of bearing inner ring fault are obvious in the result. The fault feature is highlighted and the noise component is weakened. The bearing fault type can be analyzed intuitively and effectively. Compared with the original signal which has not been processed by QPSO-SR method, the shock wave of the processed signal is obviously reduced in time domain waveform. The energy distribution of the processed signal spectrum is also more concentrated in the bearing fault frequency band. This proofs that a portion of the noise energy can be passed to the useful signal by the appropriate SR method, which is helpful for VMD to obtain the useful information more accurately.


## 5. Conclusion

This study presents a novel weak fault extraction method. This method includes QPSO-SR method and variational mode decomposition. QPSO-SR method is an adaptive denoising method that can transmit a part of the noise energy to the weak fault signal. First, QPSO-SR method is proposed to enhance the fault characteristic signal and reduce the influence from the interference noise and the shock. According to the characteristics of the original fault vibration signal, the parameters of SR system are optimized by QPSO algorithm dynamically. Thereafter, the optimal SR system is used for denosing and the processed signal is decomposed by VMD method. The failure feature information can be extracted through the FFT of the VMD component. The results obtained in the investigation of the simulated fault signal and the rolling bearing data set show that the proposed scheme avoids the influence of the interference noise and the shock. The proposed scheme reaches a higher degree of accuracy in machinery fault diagnosis compared with previous VMD method.



---

*Source: 3595871-2017-11-01.xml*