# Fault Diagnosis of Rotating Machinery Based on Adaptive Stochastic Resonance and AMD-EEMD

Authors: Peiming Shi; Cuijiao Su; Dongying Han
Journal: Shock and Vibration
 (2016)

---

## Abstract

An adaptive stochastic resonance and analytical mode decomposition-ensemble empirical mode decomposition (AMD-EEMD) method is proposed for fault diagnosis of rotating machinery in this paper. Firstly, the stochastic resonance system is optimized by particle swarm optimization (PSO), and the best structure parameters are obtained. Then, the signal with noise is put into the stochastic resonance system and denoising and enhancing the signal. Secondly, the signal output from the stochastic resonance system is extracted by analytical mode decomposition (AMD) method. Finally, the signal is decomposed by ensemble empirical mode decomposition (EEMD) method. The simulation results show that the optimal stochastic resonance system can effectively improve the signal-to-noise ratio, and the number of effective components of EEMD decomposition is significantly reduced after using AMD, thus improving the decomposition results of EEMD and enhancing the amplitude of components frequency. Through the extraction of the rolling bearing fault signal feature proved that the method has a good effect.


---

## Body

## 1. Introduction

Rotating machinery plays a significant role in a wide range of industrial applications, such as aerospace, transportation vehicles, and power generators. The rotating machine faults can cause violent vibration on the machine and even endanger normal machine operation. Therefore, accurate health monitoring and diagnosis system is needed to identify incipient fault that may occur in a rotating machine [1–4]. However, the defect-induced fault signal of the rotating machine is often corrupted by the noise coming from other coupled machine components and working environment, which makes some incipient faults not easy to be recognized. The challenge of fault recognition requires enhancing the weak fault information from heavy background noise.Stochastic resonance (SR) theory was proposed by Benzi et al. in 1981 to explain the periodicity of the Earth’s ice ages [5]. It can make part of the noise energy transfer to the low frequency signal through a nonlinear system, so weak signal drowned in the noise has been greatly strengthened at the same time to reduce the noise, and greatly improve the output signal-to-noise ratio (SNR) [6–11]. Therefore the stochastic resonance method has strong immunity to noise and can effectively extract weak signal under strong noise background. Stochastic resonance is actually the optimal matching relationship between signal, noise, and nonlinear system. However, the amplitude of the signal and the noise intensity may change over time in the actual engineering, so it is necessary to automatically adjust the nonlinear system parameters according to the different signal, in order to achieve stochastic resonance and improve the output SNR.Ensemble empirical mode decomposition (EEMD) is a noise-assisted data analysis method which was proposed based on empirical mode decomposition (EMD) [12] by Huang et al. This method uses the Gauss white noise that has the statistical characteristics of frequency distribution; the signal after adding the Gauss white noise is continuous in different scales. It can effectively solve the mode mixing problem of EMD and successfully applied in some engineering practice [13–17].Analytical modal decomposition (AMD) method was proposed by Chen and Wang [18]. This method can separate different frequency components from the signal, but the premise is to know the frequency components in the signal and then determine the bisecting frequency. Therefore, if one wants to extract the known part of the signal, the signal can be processed by AMD method.A method of rotating machinery fault diagnosis based on adaptive stochastic resonance and AMD-EEMD is proposed in this paper. Firstly, the bistable stochastic resonance system is optimized by particle swarm optimization (PSO) and putting the signal which containing noise into the stochastic resonance system and then denoising and enhancing the signal. Secondly, the signal output from the stochastic resonance system is extracted by AMD. Finally, the signal is decomposed by EEMD method. The simulation results show that the optimal stochastic resonance system can effectively improve the SNR, and the number of effective components of EEMD decomposition is significantly reduced after using AMD, thus improving the decomposition results of EEMD and enhancing the amplitude of components frequency.The rest of the paper is organized as follows. The classical bistable SR theory, the PSO algorithm, and the theories of EEMD and AMD are given and the SR-AMD-EEMD method is proposed in Section2. Section 3 describes the advantage of the proposed method by simulation experiment. The adaptive stochastic resonance and AMD-EEMD method is applied to the fault diagnosis of a rolling bearing in Section 4, and the application result proves the effectiveness of the proposed method. Finally, Section 5 provides the conclusion.


## 2. Adaptive Stochastic Resonance and AMD-EEMD Algorithm

### 2.1. Bistable Stochastic Resonance Model

Stochastic resonance is synergistic effect of the input signal and the noise signal on the nonlinear system with certain parameters, and then part of the energy of noise signal transfers to the input signal, which can improve the SNR to achieve the purpose of weak signal recognition.The Langevin equation is to describe the typical model of nonlinear bistable system [19–21]. The Langevin equation of bistable system is considered as follows:(1)dxdt=-U′x+St+Nt=ax-bx3+St+Nt,where a and b are the structural parameters of the bistable system, S(t) is the input signal, and N(t) is the noise signal. Let N(t)=0, N(t)N(t-τ)=2Dδ(t), and N(t)=2Dξ(t), where ξ(t) is white noise with zero mean and unit variance and D is the noise intensity. The potential function for the above bistable system can be denoted as(2)Ux=-12ax2+14bx4.The above equation has two minimum values inx=±a/b and has a maximum value in x=0, corresponding to the system of two lowest points and a barrier point, corresponding to the two potential hydrazine points and a barrier point of the system. And its barrier height is ΔU=a2/(4b).


### 2.2. Particle Swarm Optimization Algorithm

Particle swarm optimization (PSO) algorithm is an evolutionary computing technology proposed by Eberhart and Kennedy. Assumingi particle position and velocity in d dimensional search space are Xi=(xi1,xi2,…,xid) and Vi=(vi1,vi2,…,vid), respectively, in each iteration, particles update their own by tracking two optimal solutions, one is the optimal solution Pi=(pi1,pi2,…,pid) found by particle itself, and the other is the optimal solution of current population Pg. When finding the two best values, the particle updates its velocity and position according to the following formula:(3)vi,jt+1=ωvi,jt+c1r1Pi,j-xi,jt+c2r2Pg,j-xi,jt,xi,jt+1=xi,jt+vi,jt+1,j=1,2,…,d,where ω is inertia weight, c1 and c2 are the learning factors, and r1 and r2 are random numbers between 0 and 1 uniform distribution [22, 23].Because the larger weights are conducive to jumping out of the local minimum value, they are suitable for global search. And smaller weights are conducive for accurate local search in the current search area, which is helpful in the convergence of the algorithm.Therefore, in order to prevent the phenomenon of premature convergence and oscillation occurred in the vicinity of the global optimal solution of PSO algorithm, you can use the linear decreasing weight as follows:(4)ω=ωmax-t∗ωmax-ωmintmax,where ωmax and ωmin are maximum value and minimum value of ω, t is the current number of iterations, and tmax is the maximum number of iterations. Usually, ωmax=0.9 and ωmin=0.4.


### 2.3. EEMD Based on AMD

A new signal decomposition method called AMD was proposed by Chen and Wang in 2012 [18]. This method can decompose signals with closely spaced frequency components, which means being able to extract the signal.If a time seriesx(t) is divided into two signals by bisecting frequency ωb(5)xt=s1t+s-1t.We can get that(6)s1t=sin⁡ωbtHxtcos⁡ωbt-cos⁡ωbtHxtsin⁡ωbt,(7)s-1t=xt-s1t,where H· is Hilbert transformation of the function in brackets.So the original signal will get the two signals by using the AMD method.The essence of EEMD method is repeated EMD decomposition added Gauss white noise. Finally, the original signal is decomposed into several intrinsic mode functions (IMFs) and a residue.The specific steps of AMD-EEMD method are as follows:(1)
The signal is extracted by AMD method firstly. Assuming there is a time seriesx(t)=x1(t)+x2(t)+x3(t), the frequencies are f1,f2,f3, respectively, and f1<f2<f3. If one wants to extract f1 frequency component of the signal, take between f1 and f2 as the bisecting frequency. If one wants to extract f2 frequency component of the signal, decompose the first part of the signal which is taking frequency value between f2 and f3 as the bisecting frequency, and then subtract the first half part of the signal which is taking frequency value between f1 and f2 as the bisecting frequency.(2)
Assuming that the signal after the extraction isx2(t), one uses EEMD method to decompose it.(3)
Initialize the times of EMD decompositionM and white noise amplitude coefficient K, and let m=1.(4)
Perform them time of EMD decomposition. Consider the following:
(a)
Add a random sequence of Gauss white noiseni(t) on the input signal x2(t), and then get the signal with noise x2(t)=x2(t)+k∗ni(t).(b)
Use EMD to decomposex2(t) and get I  IMFjm  (j=1,2,…,I), jm being the decomposition of the first j IMF in the m time of EMD test.(c)
Ifm<M, return to step (4), and let m=m+1.(5)
Calculate the mean value of each IMF inM times’ test.(6)
Take the output ofIMFj as the first j IMF decomposed by EEMD.Usually, M is 100, K takes 0.01~0.5 times the standard deviation of the original signal which is more appropriate.


### 2.4. Adaptive Stochastic Resonance and AMD-EEMD Algorithm

Stochastic resonance output SNR formula is defined as follows:(8)SNR=10lg⁡SN=10lg⁡2Xk02∑k=0L-1Xk2-2Xk02,where S is the signal power, N is the noise power, L is data length, and X(k0) is the unilateral amplitude of output signal frequency component.When the system input is a small parameter signal, the output SNR approximately is as follows:(9)SNR=2a2A2exp⁡-a2/4bD4bD2.Equation (8) shows that the output SNR is relevant to signal amplitude A, noise intensity D, and system parameters a and b. A signal amplitude as well as the noise intensity is known, so if one wants to get the maximum output SNR one needs to adjust the parameters of a and b.A method of rotating machinery fault diagnosis based on adaptive stochastic resonance and AMD-EEMD is proposed in this paper. Firstly, the bistable stochastic resonance system is optimized by particle swarm optimization (PSO). Secondly, the signal output from the stochastic resonance system is intercepted and decomposed by AMD-EEMD method. The main steps are as follows:(1)
Population initialization: Set up the PSO maximum generationtmax, population quantity M, learning factors c1 and c2, upper limit of inertia weight ωmax, lower limit of inertia weight ωmin, upper limit of flight speed Vmax, lower limit of flight speed Vmin, and the search range of the system parameters a and b. Initialize a group of particles’ position and velocity randomly.(2)
The fitness evaluation of each particle: Calculate every particle’s fitness valuefitness(x) according to the SNR. Store the current position and fitness value of each particle inpbestof each particle. Then, store the position and fitness value of the best individual of allpbest in thegbest.(3)
The velocity and position of particles being updated: The position and velocity of all particles are updated according to (3), and the weight value is updated according to (4). Recalculate every particle’s fitness value. If local optimal fitness valuepbest or global optimal fitness valuegbest of the individual particle is better than thepbest andgbest of the generation particles, then update the corresponding individual particle’s local optimal fitness value or global optimal fitness value.(4)
To determine whether the end condition is met: When the evolution times reached the maximum numbertmax stop searching. Then output the optimal system parameters a and b and put the original signal into stochastic resonance system with the optimal parameters.(5)
Using the AMD-EEMD method to intercept and decompose the signal which is output from the stochastic resonance system.The specific process is shown as Figure1.Figure 1
The flow chart of adaptive stochastic resonance and AMD-EEMD.


## 2.1. Bistable Stochastic Resonance Model

Stochastic resonance is synergistic effect of the input signal and the noise signal on the nonlinear system with certain parameters, and then part of the energy of noise signal transfers to the input signal, which can improve the SNR to achieve the purpose of weak signal recognition.The Langevin equation is to describe the typical model of nonlinear bistable system [19–21]. The Langevin equation of bistable system is considered as follows:(1)dxdt=-U′x+St+Nt=ax-bx3+St+Nt,where a and b are the structural parameters of the bistable system, S(t) is the input signal, and N(t) is the noise signal. Let N(t)=0, N(t)N(t-τ)=2Dδ(t), and N(t)=2Dξ(t), where ξ(t) is white noise with zero mean and unit variance and D is the noise intensity. The potential function for the above bistable system can be denoted as(2)Ux=-12ax2+14bx4.The above equation has two minimum values inx=±a/b and has a maximum value in x=0, corresponding to the system of two lowest points and a barrier point, corresponding to the two potential hydrazine points and a barrier point of the system. And its barrier height is ΔU=a2/(4b).


## 2.2. Particle Swarm Optimization Algorithm

Particle swarm optimization (PSO) algorithm is an evolutionary computing technology proposed by Eberhart and Kennedy. Assumingi particle position and velocity in d dimensional search space are Xi=(xi1,xi2,…,xid) and Vi=(vi1,vi2,…,vid), respectively, in each iteration, particles update their own by tracking two optimal solutions, one is the optimal solution Pi=(pi1,pi2,…,pid) found by particle itself, and the other is the optimal solution of current population Pg. When finding the two best values, the particle updates its velocity and position according to the following formula:(3)vi,jt+1=ωvi,jt+c1r1Pi,j-xi,jt+c2r2Pg,j-xi,jt,xi,jt+1=xi,jt+vi,jt+1,j=1,2,…,d,where ω is inertia weight, c1 and c2 are the learning factors, and r1 and r2 are random numbers between 0 and 1 uniform distribution [22, 23].Because the larger weights are conducive to jumping out of the local minimum value, they are suitable for global search. And smaller weights are conducive for accurate local search in the current search area, which is helpful in the convergence of the algorithm.Therefore, in order to prevent the phenomenon of premature convergence and oscillation occurred in the vicinity of the global optimal solution of PSO algorithm, you can use the linear decreasing weight as follows:(4)ω=ωmax-t∗ωmax-ωmintmax,where ωmax and ωmin are maximum value and minimum value of ω, t is the current number of iterations, and tmax is the maximum number of iterations. Usually, ωmax=0.9 and ωmin=0.4.


## 2.3. EEMD Based on AMD

A new signal decomposition method called AMD was proposed by Chen and Wang in 2012 [18]. This method can decompose signals with closely spaced frequency components, which means being able to extract the signal.If a time seriesx(t) is divided into two signals by bisecting frequency ωb(5)xt=s1t+s-1t.We can get that(6)s1t=sin⁡ωbtHxtcos⁡ωbt-cos⁡ωbtHxtsin⁡ωbt,(7)s-1t=xt-s1t,where H· is Hilbert transformation of the function in brackets.So the original signal will get the two signals by using the AMD method.The essence of EEMD method is repeated EMD decomposition added Gauss white noise. Finally, the original signal is decomposed into several intrinsic mode functions (IMFs) and a residue.The specific steps of AMD-EEMD method are as follows:(1)
The signal is extracted by AMD method firstly. Assuming there is a time seriesx(t)=x1(t)+x2(t)+x3(t), the frequencies are f1,f2,f3, respectively, and f1<f2<f3. If one wants to extract f1 frequency component of the signal, take between f1 and f2 as the bisecting frequency. If one wants to extract f2 frequency component of the signal, decompose the first part of the signal which is taking frequency value between f2 and f3 as the bisecting frequency, and then subtract the first half part of the signal which is taking frequency value between f1 and f2 as the bisecting frequency.(2)
Assuming that the signal after the extraction isx2(t), one uses EEMD method to decompose it.(3)
Initialize the times of EMD decompositionM and white noise amplitude coefficient K, and let m=1.(4)
Perform them time of EMD decomposition. Consider the following:
(a)
Add a random sequence of Gauss white noiseni(t) on the input signal x2(t), and then get the signal with noise x2(t)=x2(t)+k∗ni(t).(b)
Use EMD to decomposex2(t) and get I  IMFjm  (j=1,2,…,I), jm being the decomposition of the first j IMF in the m time of EMD test.(c)
Ifm<M, return to step (4), and let m=m+1.(5)
Calculate the mean value of each IMF inM times’ test.(6)
Take the output ofIMFj as the first j IMF decomposed by EEMD.Usually, M is 100, K takes 0.01~0.5 times the standard deviation of the original signal which is more appropriate.


## 2.4. Adaptive Stochastic Resonance and AMD-EEMD Algorithm

Stochastic resonance output SNR formula is defined as follows:(8)SNR=10lg⁡SN=10lg⁡2Xk02∑k=0L-1Xk2-2Xk02,where S is the signal power, N is the noise power, L is data length, and X(k0) is the unilateral amplitude of output signal frequency component.When the system input is a small parameter signal, the output SNR approximately is as follows:(9)SNR=2a2A2exp⁡-a2/4bD4bD2.Equation (8) shows that the output SNR is relevant to signal amplitude A, noise intensity D, and system parameters a and b. A signal amplitude as well as the noise intensity is known, so if one wants to get the maximum output SNR one needs to adjust the parameters of a and b.A method of rotating machinery fault diagnosis based on adaptive stochastic resonance and AMD-EEMD is proposed in this paper. Firstly, the bistable stochastic resonance system is optimized by particle swarm optimization (PSO). Secondly, the signal output from the stochastic resonance system is intercepted and decomposed by AMD-EEMD method. The main steps are as follows:(1)
Population initialization: Set up the PSO maximum generationtmax, population quantity M, learning factors c1 and c2, upper limit of inertia weight ωmax, lower limit of inertia weight ωmin, upper limit of flight speed Vmax, lower limit of flight speed Vmin, and the search range of the system parameters a and b. Initialize a group of particles’ position and velocity randomly.(2)
The fitness evaluation of each particle: Calculate every particle’s fitness valuefitness(x) according to the SNR. Store the current position and fitness value of each particle inpbestof each particle. Then, store the position and fitness value of the best individual of allpbest in thegbest.(3)
The velocity and position of particles being updated: The position and velocity of all particles are updated according to (3), and the weight value is updated according to (4). Recalculate every particle’s fitness value. If local optimal fitness valuepbest or global optimal fitness valuegbest of the individual particle is better than thepbest andgbest of the generation particles, then update the corresponding individual particle’s local optimal fitness value or global optimal fitness value.(4)
To determine whether the end condition is met: When the evolution times reached the maximum numbertmax stop searching. Then output the optimal system parameters a and b and put the original signal into stochastic resonance system with the optimal parameters.(5)
Using the AMD-EEMD method to intercept and decompose the signal which is output from the stochastic resonance system.The specific process is shown as Figure1.Figure 1
The flow chart of adaptive stochastic resonance and AMD-EEMD.


## 3. Simulation Experiment

Assuming the input signal is(10)xt=Asin⁡2πf0t+Nt,where A=0.3, f0=0.01 Hz, the strength D of Gauss white noise N(t) is 0.31, the sampling frequency fs is 5 Hz, the sampling points are 5000, the signal is decomposed by EEMD method directly and the first 8 components are taken, and the result is shown as Figure 2.Figure 2
The EEMD components of original signal and its spectrum.First of all, use particle swarm optimization algorithm to optimize parametersa and b of stochastic resonance system. Set the population number as 40, the search range of a and b as [0,2], the maximum value of flight speed as 0.01, and the maximum number of evolutions as 100. Because the signal amplitude as well as intensity of noise is known, (6) is regarded as the fitness function. From Figure 3 one can get that the maximum output SNR is 0.151 and the optimization result is a=0.7381 and b=0.4393. Use stochastic resonance system with the optimal parameters to process the original signal, and then get time-domain waveform and frequency spectrum which are shown as Figures 4(c) and 4(d). Figure 4 shows that the noise has been weakened and the 0.01 Hz frequency component has been greatly improved.Figure 3
The optimal convergence curve of PSO algorithm.Figure 4
The time-domain waveform and frequency spectrum of original signal and stochastic resonance signal.
(a)
The original signal waveform
(b)
The original signal spectrum
(c)
The SR signal waveform
(d)
The SR signal spectrumThe 0~0.02 Hz band part of the original signal and the signal output from stochastic resonance system are extracted by the AMD method, respectively. Then they are decomposed by EEMD method and the frequency spectrum of each component is calculated; the decomposition results are shown as Figures5 and 6. Comparing with Figure 2, the component is significantly reduced in Figure 5. Figure 5 shows that component C1 contains the signal whose frequency is 0.01 Hz and its frequency amplitude is 0.2716. Figure 6 shows that also component C1 contains the signal whose frequency is 0.01 Hz, but its frequency amplitude is 1.461. This shows that the original signal is processed by stochastic resonance and then EEMD decomposition, which can make the signal whose frequency is 0.01 Hz outstanding.Figure 5
The AMD-EEMD components of original signal and its spectrum.Figure 6
The AMD-EEMD components of stochastic resonance signal and its spectrum.


## 4. Application Example

For the Case Western Reserve University rolling bearing fault simulation test bench, the type of drive end bearing is SKF6205, the motor speed is 1797 r/min (29.95 Hz), and the sampling frequency is 12 kHz. The main parameters of the rolling bearing and the rolling bearing fault feature frequency of the different parts are shown as Tables1 and 2, respectively.Table 1
The main parameters of the rolling bearing of SKF6205.
Out diameterR (mm)
Inner diameterr (mm)
Pitch diameterD (mm)
Ball diameterd (mm)
Ball numberZ
Contact anglea
52.00
25.00
39.04
7.94
9
0°Table 2
Rolling bearing fault feature frequency of SKF6205.
Bearing element
Inner ring
Outer ring
The retainer
Rolling body
Feature frequency
5.4152fr
3.5848fr
0.39828fr
4.7135fr
Note.  fr=N/60 and N is bearing speed.Take the inner race fault of rolling bearing as an example; select a set of data for analysis and processing. By calculating the feature frequency is 162.18 Hz. The time-domain waveform of the bearing signal and its spectrum are shown as Figure7. The original signal is decomposed by EEMD method and the first 8 components are taken; the result is shown as Figure 8.Figure 7
The bearing fault signal and its spectrum.Figure 8
The EEMD components of original fault signal and its spectrum.Because the bearing fault signal does not meet the requirements of the small parameter, the signal needs sampling second time. Set up the second sampling frequencyfsr=fs/R=4 Hz; that is to say the transformation scale R is 3000. Use particle swarm algorithm to optimize the parameters of stochastic resonance system. The noise strength of the bearing fault signal is unknown, so (5) is regarded as the fitness function. After optimizing by PSO algorithm the optimal parameters are a=0.15 and b=2.14. Then put the bearing signal into stochastic resonance system with the optimal parameters.The 150~180 Hz band part of the original fault signal and the signal output from stochastic resonance system are extracted by the AMD method, respectively. Then they are decomposed by EEMD method and the frequency spectrum of each component is calculated. The decomposition results are shown as Figures9 and 10. Comparing with Figure 8, the component is significantly reduced in Figure 9.Figure 9
The AMD-EEMD components of original fault signal and its spectrum.Figure 10
The AMD-EEMD components of stochastic resonance signal and its spectrum.Figure9 shows that the main frequency component of C2 is 162 Hz, according to the inner race fault feature frequency, and its frequency amplitude is 0.01481. Figure 10 shows that the main frequency component of C2 is also 162 Hz, but its frequency amplitude is 0.02887, improving significantly compared with Figure 9.The comparative analysis outlined above indicates that the adaptive stochastic resonance and AMD-EEMD method may not only increase the output SNR but also reduce the effective components of EEMD and improve the amplitude of components frequency. Therefore, the adaptive stochastic resonance and AMD-EEMD method and its application in fault diagnosis of rotating machinery have practical significance.For the Case Western Reserve University rolling bearing fault simulation test bench, the other type of fan end bearing is SKF6203 and the sampling frequency is 12 kHz. The main parameters of the rolling bearing and the rolling bearing fault feature frequency of the different parts are shown as Tables3 and 4, respectively.Table 3
The main parameters of the rolling bearing of SKF6203.
Out diameterR (mm)
Inner diameterr (mm)
Pitch diameterD (mm)
Ball diameterd (mm)
Ball numberZ
Contact anglea
40.00
17.00
28.50
6.75
8
0°Table 4
Rolling bearing fault feature frequency of SKF6203.
Bearing element
Inner ring
Outer ring
The retainer
Rolling body
Feature frequency
4.9469fr
3.0530fr
0.3817fr
3.9874fr
Note.  fr=N/60 and N is bearing speed.Forty-four groups of bearing inner fault data were analyzed by EMD, EEMD, SR-EMD, SR-EEMD, AMD-EMD, AMD-EEMD, and SR-AMD-EEMD method, respectively. Among them, forty-two groups of data can effectively increase the amplitude and reduce the number of components after SR-AMD-EEMD method; that is, the accuracy of this method is 95.45%.Here take eight groups of bearing inner race fault of SKF6205 and SKF6203 that are analyzed under 1797 r/min, 1772 r/min, 1750 r/min, and 1730 r/min different speed, respectively. The frequency amplitude and component count are shown in Tables5 and 6. From Tables 5 and 6, we can clearly see that the EMD and EEMD methods have similar amplitude and number of components, the SR-EMD and SR-EEMD methods increase the amplitude, but the number of the components is not changed, the AMD-EMD and AMD-EEMD methods reduce the number of components, but the amplitude is almost unchanged, and the SR-AMD-EEMD method has the largest amplitude, while the component is the least. This shows the effectiveness of the proposed method in this paper.Table 5
The amplitude and component count of fault feature frequency of SKF6205 under different rotational speed.
Bearing type
SKF6205
Motor speed
1797 r/min
1772 r/min
1750 r/min
1730 r/min
Frequency value
162 Hz
159.6 Hz
157.2 Hz
154.8 Hz
Amplitude
Component count
Amplitude
Component count
Amplitude
Component count
Amplitude
Component count
SR-AMD-EEMD
0.02887
3
0.02164
2
0.02567
3
0.02608
2
AMD-EEMD
0.01481
3
0.01402
2
0.01518
3
0.00881
3
AMD-EMD
0.01536
3
0.01336
2
0.01528
4
0.00863
2
SR-EEMD
0.02880
11
0.02034
10
0.02555
10
0.02480
9
SR-EMD
0.02579
9
0.01738
9
0.02001
8
0.02729
10
EEMD
0.01324
10
0.01672
10
0.01314
10
0.00834
11
EMD
0.01163
10
0.01371
9
0.01224
9
0.00763
10Table 6
The amplitude and component count of fault feature frequency of SKF6203 under different rotational speed.
Bearing type
SKF6203
Motor speed
1797 r/min
1772 r/min
1750 r/min
1730 r/min
Frequency value
147.6 Hz
146.4 Hz
144 Hz
142.8 Hz
Amplitude
Component count
Amplitude
Component count
Amplitude
Component count
Amplitude
Component count
SR-AMD-EEMD
0.01287
4
0.00858
3
0.00754
3
0.00748
4
AMD-EEMD
0.00524
4
0.00342
3
0.00323
4
0.00375
4
AMD-EMD
0.00529
4
0.00345
3
0.00335
4
0.00382
4
SR-EEMD
0.01208
10
0.00854
10
—
—
—
—
SR-EMD
0.01228
9
0.00660
9
—
—
—
—
EEMD
0.00574
10
0.00364
10
—
—
—
—
EMD
0.00568
9
0.00469
9
—
—
—
—
Note. “—” expresses no result.Calculate the average value of the amplitude increment of eight groups which is compared with the EMD method. The percentage increase is shown in Figure11. And the average number of components is shown in Figure 12. Figure 11 shows that the amplitude increased by about 10% after EEMD, AMD-EMD, and AMD-EEMD method compared with EMD method. And the amplitude increased by 96.74%, 116.10%, and 122.58% after SR-EMD, SR-EEMD, and SR-AMD-EEMD method, respectively. From Figure 12 we can see that the minimum average number of components is 3 processed by the SR-AMD-EEMD method. To sum up we know that the SR-AMD-EEMD method has the largest amplitude and the minimum number of components.Figure 11
The column chart of percentage increase.Figure 12
The average number of components.


## 5. Conclusions

Consider the following:(1)
This paper presents an adaptive stochastic resonance and AMD-EEMD method for fault diagnosis of rotating machinery. The process is as follows: firstly, the bistable stochastic resonance system is optimized by PSO, and the best structure parameters are obtained; then, the signal with noise is put into the stochastic resonance system in order to improve the SNR of the signal; secondly, the signal output from the stochastic resonance system is extracted by AMD method; finally, the signal is decomposed by EEMD method.(2)
By processing the simulation signal and application examples and comparing the SR-AMD-EEMD method and EMD, EEMD, SR-EMD, SR-EEMD, AMD-EMD, and AMD-EEMD method, respectively, the result shows that the SR-AMD-EEMD method can effectively improve the SNR and also reduce the number of EEMD effective components, thus improving the quality and accuracy of EEMD decomposition and enhancing the amplitude of components frequency. Therefore, it is valuable in engineering applications.



---

*Source: 9278581-2016-02-09.xml*