Vibration Analysis as a Diagnosis Tool for Health Monitoring of Industrial Machines
View this Special IssueResearch Article  Open Access
Peiming Shi, Cuijiao Su, Dongying Han, "Fault Diagnosis of Rotating Machinery Based on Adaptive Stochastic Resonance and AMDEEMD", Shock and Vibration, vol. 2016, Article ID 9278581, 11 pages, 2016. https://doi.org/10.1155/2016/9278581
Fault Diagnosis of Rotating Machinery Based on Adaptive Stochastic Resonance and AMDEEMD
Abstract
An adaptive stochastic resonance and analytical mode decompositionensemble empirical mode decomposition (AMDEEMD) 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 signaltonoise 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.
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 defectinduced 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 signaltonoise 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 noiseassisted 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 AMDEEMD 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 SRAMDEEMD method is proposed in Section 2. Section 3 describes the advantage of the proposed method by simulation experiment. The adaptive stochastic resonance and AMDEEMD 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 AMDEEMD 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:where and are the structural parameters of the bistable system, is the input signal, and is the noise signal. Let , , and , where is white noise with zero mean and unit variance and is the noise intensity. The potential function for the above bistable system can be denoted as
The above equation has two minimum values in and has a maximum value in , 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 .
2.2. Particle Swarm Optimization Algorithm
Particle swarm optimization (PSO) algorithm is an evolutionary computing technology proposed by Eberhart and Kennedy. Assuming particle position and velocity in dimensional search space are and , respectively, in each iteration, particles update their own by tracking two optimal solutions, one is the optimal solution found by particle itself, and the other is the optimal solution of current population . When finding the two best values, the particle updates its velocity and position according to the following formula:where is inertia weight, and are the learning factors, and and 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:where and are maximum value and minimum value of , is the current number of iterations, and is the maximum number of iterations. Usually, and .
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 series is divided into two signals by bisecting frequency We can get thatwhere 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 AMDEEMD method are as follows:(1)The signal is extracted by AMD method firstly. Assuming there is a time series , the frequencies are , respectively, and . If one wants to extract frequency component of the signal, take between and as the bisecting frequency. If one wants to extract frequency component of the signal, decompose the first part of the signal which is taking frequency value between and as the bisecting frequency, and then subtract the first half part of the signal which is taking frequency value between and as the bisecting frequency.(2)Assuming that the signal after the extraction is , one uses EEMD method to decompose it.(3)Initialize the times of EMD decomposition and white noise amplitude coefficient , and let .(4)Perform the time of EMD decomposition. Consider the following:(a)Add a random sequence of Gauss white noise on the input signal , and then get the signal with noise .(b)Use EMD to decompose and get , being the decomposition of the first IMF in the time of EMD test.(c)If , return to step (4), and let .(5)Calculate the mean value of each IMF in times’ test.(6)Take the output of as the first IMF decomposed by EEMD.Usually, is 100, takes 0.01~0.5 times the standard deviation of the original signal which is more appropriate.
2.4. Adaptive Stochastic Resonance and AMDEEMD Algorithm
Stochastic resonance output SNR formula is defined as follows:where is the signal power, is the noise power, is data length, and 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:Equation (8) shows that the output SNR is relevant to signal amplitude , noise intensity , and system parameters and . 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 and .
A method of rotating machinery fault diagnosis based on adaptive stochastic resonance and AMDEEMD 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 AMDEEMD method. The main steps are as follows:(1)Population initialization: Set up the PSO maximum generation , population quantity , learning factors and , upper limit of inertia weight , lower limit of inertia weight , upper limit of flight speed , lower limit of flight speed , and the search range of the system parameters and . Initialize a group of particles’ position and velocity randomly.(2)The fitness evaluation of each particle: Calculate every particle’s fitness value according to the SNR. Store the current position and fitness value of each particle in pbest of each particle. Then, store the position and fitness value of the best individual of all pbest in the gbest.(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 value pbest or global optimal fitness value gbest of the individual particle is better than the pbest and gbest 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 number stop searching. Then output the optimal system parameters and and put the original signal into stochastic resonance system with the optimal parameters.(5)Using the AMDEEMD method to intercept and decompose the signal which is output from the stochastic resonance system.
The specific process is shown as Figure 1.
3. Simulation Experiment
Assuming the input signal is where , Hz, the strength of Gauss white noise is 0.31, the sampling frequency 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.
First of all, use particle swarm optimization algorithm to optimize parameters and of stochastic resonance system. Set the population number as 40, the search range of and as , 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 and . Use stochastic resonance system with the optimal parameters to process the original signal, and then get timedomain 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.
(a) The original signal waveform
(b) The original signal spectrum
(c) The SR signal waveform
(d) The SR signal spectrum
The 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 Figures 5 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.
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 Tables 1 and 2, respectively.

 
Note. and 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 timedomain waveform of the bearing signal and its spectrum are shown as Figure 7. The original signal is decomposed by EEMD method and the first 8 components are taken; the result is shown as Figure 8.
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 frequency Hz; that is to say the transformation scale 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 and . 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 Figures 9 and 10. Comparing with Figure 8, the component is significantly reduced in Figure 9.
Figure 9 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 AMDEEMD 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 AMDEEMD 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 Tables 3 and 4, respectively.

 
Note. and is bearing speed. 
Fortyfour groups of bearing inner fault data were analyzed by EMD, EEMD, SREMD, SREEMD, AMDEMD, AMDEEMD, and SRAMDEEMD method, respectively. Among them, fortytwo groups of data can effectively increase the amplitude and reduce the number of components after SRAMDEEMD 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 Tables 5 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 SREMD and SREEMD methods increase the amplitude, but the number of the components is not changed, the AMDEMD and AMDEEMD methods reduce the number of components, but the amplitude is almost unchanged, and the SRAMDEEMD method has the largest amplitude, while the component is the least. This shows the effectiveness of the proposed method in this paper.

 
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 Figure 11. And the average number of components is shown in Figure 12. Figure 11 shows that the amplitude increased by about 10% after EEMD, AMDEMD, and AMDEEMD method compared with EMD method. And the amplitude increased by 96.74%, 116.10%, and 122.58% after SREMD, SREEMD, and SRAMDEEMD method, respectively. From Figure 12 we can see that the minimum average number of components is 3 processed by the SRAMDEEMD method. To sum up we know that the SRAMDEEMD method has the largest amplitude and the minimum number of components.
5. Conclusions
Consider the following:(1)This paper presents an adaptive stochastic resonance and AMDEEMD 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 SRAMDEEMD method and EMD, EEMD, SREMD, SREEMD, AMDEMD, and AMDEEMD method, respectively, the result shows that the SRAMDEEMD 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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant no. 51475407), Hebei Provincial Natural Science Foundation of China (no. E2015203190), and Key Project of Natural Science Research in Colleges and Universities of Hebei Province (Grant no. ZD2015050).
References
 Y. Tian, J. Ma, C. Lu, and Z. Wang, “Rolling bearing fault diagnosis under variable conditions using LMDSVD and extreme learning machine,” Mechanism and Machine Theory, vol. 90, pp. 175–186, 2015. View at: Publisher Site  Google Scholar
 J. Zheng, J. Cheng, and Y. Yang, “A rolling bearing fault diagnosis approach based on LCD and fuzzy entropy,” Mechanism & Machine Theory, vol. 70, no. 6, pp. 441–453, 2013. View at: Publisher Site  Google Scholar
 H. Liu and M. Han, “A fault diagnosis method based on local mean decomposition and multiscale entropy for roller bearings,” Mechanism & Machine Theory, vol. 75, pp. 67–78, 2014. View at: Publisher Site  Google Scholar
 J. Zheng, J. Cheng, Y. Yang, and S. Luo, “A rolling bearing fault diagnosis method based on multiscale fuzzy entropy and variable predictive modelbased class discrimination,” Mechanism and Machine Theory, vol. 78, pp. 187–200, 2014. View at: Publisher Site  Google Scholar
 R. Benzi, A. Sutera, and A. Vulpiani, “The mechanism of stochastic resonance,” Journal of Physics A: Mathematical and General, vol. 14, no. 11, pp. L453–L457, 1981. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Hang and H. Yan, “Research on weak signal detection based on cascaded stochastic bistable resonance,” Computer Measurement & Control, vol. 20, no. 7, pp. 1821–1836, 2012. View at: Google Scholar
 P. Shi, X. Ding, and D. Han, “Study on multifrequency weak signal detection method based on stochastic resonance tuning by multiscale noise,” Measurement, vol. 47, no. 1, pp. 540–546, 2014. View at: Publisher Site  Google Scholar
 Y. G. Leng and T. Y. Wang, “Numerical research of twice sampling stochastic resonance for the detection of a weak signal submerged in a heavy noise,” Acta Physica Sinica, vol. 52, no. 10, pp. 2432–2437, 2003. View at: Google Scholar
 S. Lu, Q. He, and F. Kong, “Effects of underdamped stepvarying secondorder stochastic resonance for weak signal detection,” Digital Signal Processing, vol. 36, pp. 93–103, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 Q. He and J. Wang, “Effects of multiscale noise tuning on stochastic resonance for weak signal detection,” Digital Signal Processing, vol. 22, no. 4, pp. 614–621, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 L. Zhang, F. Wang, and Z. Hui, “Detecting high frequency weak signal using parameteradjusted frequency shifting based on stochastic resonance (SR),” Physics Procedia, vol. 24, pp. 382–389, 2012. View at: Publisher Site  Google Scholar
 N. E. Huang, Z. Shen, S. R. Long et al., “The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis,” Proceedings of the Royal Society of London, Series A: Mathematical, Physical and Engineering Sciences, vol. 454, no. 1971, pp. 903–995, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Wu and N. E. Huang, “Ensemble empirical mode decomposition: a noiseassisted data analysis method,” Advances in Adaptive Data Analysis, vol. 1, no. 1, pp. 1–49, 2011. View at: Publisher Site  Google Scholar
 Y.H. Wang, C.H. Yeh, H.W. V. Young, K. Hu, and M.T. Lo, “On the computational complexity of the empirical mode decomposition algorithm,” Physica A: Statistical Mechanics & Its Applications, vol. 400, no. 2, pp. 159–167, 2014. View at: Publisher Site  Google Scholar
 W. JinDun, C. Lue, and Q. Wei, Research on Fault Weak Signal Feature Extraction Based on EEMD Method, Electronic Design Engineering, Norco, Calif, USA, 2012.
 Y. Lei, N. Li, J. Lin, and S. Wang, “Fault diagnosis of rotating machinery based on an adaptive ensemble empirical mode decomposition,” Sensors, vol. 13, no. 12, pp. 16950–16964, 2013. View at: Publisher Site  Google Scholar
 F. Zhao, X. Mei, T. Tao, G. Jiang, and Y. Zhou, “Fault diagnosis of a machine tool rotary axis based on a motor current test and the ensemble empirical mode decomposition method,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 225, no. 5, pp. 1121–1129, 2011. View at: Publisher Site  Google Scholar
 G. Chen and Z. Wang, “A signal decomposition theorem with Hilbert transform and its application to narrowband time series with closely spaced frequency components,” Mechanical Systems and Signal Processing, vol. 28, pp. 258–279, 2012. View at: Publisher Site  Google Scholar
 S. Fauve and F. Heslot, “Stochastic resonance in a bistable system,” Physics Letters A, vol. 97, no. 12, pp. 5–7, 1983. View at: Publisher Site  Google Scholar
 Y.G. Leng, T.Y. Wang, Y. Guo, Y.G. Xu, and S.B. Fan, “Engineering signal processing based on bistable stochastic resonance,” Mechanical Systems & Signal Processing, vol. 21, no. 1, pp. 138–150, 2007. View at: Publisher Site  Google Scholar
 Y. Qin, Y. Tao, Y. He, and B. Tang, “Adaptive bistable stochastic resonance and its application in mechanical fault feature extraction,” Journal of Sound and Vibration, vol. 333, no. 26, pp. 7386–7400, 2014. View at: Publisher Site  Google Scholar
 J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proceedings of the IEEE International Conference on Neural Networks, vol. 4, pp. 129–132, Perth, Australia, December 1995. View at: Publisher Site  Google Scholar
 Y. Shi and R. C. Eberhart, “Parameter selection in particle swarm optimization,” in Evolutionary Programming VII, vol. 1447 of Lecture Notes in Computer Science, pp. 591–600, Springer, Berlin, Germany, 1998. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Peiming Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.