Periodic non-sinusoidal time-delay stochastic resonance weak fault diagnosis method and its application
Zhixing Li1 , Huiru Bao2 , Songjiu Han3
1Beijing University of Civil Engineering and Architecture, Beijing, China
1, 3School of Mechanical Engineering, University of Science and Technology Inner Mongolia, Baotou, China
2Baotou Vocational and Technical College, Baotou, China
Journal of Vibroengineering, Vol. 23, Issue 4, 2021, p. 847-860.
Received 2 June 2020; received in revised form 5 November 2020; accepted 24 November 2020; published 22 June 2021
Stochastic resonance is the use of nonlinear systems to synchronize an original signal with noise. This method is commonly used to extract useful signals by reducing noise and has been widely used for mechanical weak fault diagnosis. This paper analyzes the characteristics of a periodic non-sinusoidal potential function, considers the shape of the model, and introduces a time-delay. The steady-state probability density function, effective potential function, and signal-to-noise ratio are then analyzed. As a result, a signal detection method for periodic non-sinusoidal time-delay stochastic resonance (PNTSR) is proposed. Experimental and engineering data are used to explain the PNTSR through the simulation. It is found that the PNTSR method has better fault detection results when compared to the classic bi-stable stochastic resonance method.
- A periodic non-sinusoidal potential function
- A signal detection method for time-delay stochastic resonance
- The model considers the shape of the model and introduces a time-delay
Keywords: stochastic resonance, periodic non-sinusoidal potential, time-delay feedback, weak signal, fault detection.
Fault detection technology in industrial production for mechanical equipment has received recent attention. This is because equipment failure can cause equipment damage and also affect personnel safety [1, 2]. Because large mechanical equipment is often operated at low-speed in heavy-duty and noisy working environments, it can be difficult to identify early warnings of failure . Early fault signals in key parts (e.g., bearings) are often difficult to detect in noisy settings and can lead to inaccurate fault diagnosis . Therefore, for early fault detection, the ability to extract weak fault characteristics within a noisy background is important .
A variety of signal processing methods have been proposed and studied. These include empirical mode decomposition [6, 7], wavelet analysis [8, 9], and singular value decomposition [10, 11]. These methods are primarily based on noise reduction, having the potential to weaken the weak fault characteristic signal while processing noise . Some scholars have developed fault detection systems based on traditional signal processing algorithms. These include supervised learning systems that predict performance and remaining life of bearings and other components [13-15]. Stochastic resonance (SR), on the other hand, uses a different approach for filtering noise. It adds moderate noise that can excite the particles in the nonlinear system and enhance the amplitude and power of the weak signal.
SR was first proposed by Benzi  as a method to detect weak faults in paleo-meteorological glaciers. Since then, SR has been applied to meteorology, biology, and physics [17-21]. Recently, SR has been used for mechanical fault detection [22, 23]. Early fault detection SR studies aimed to overcome small parameter limits (e.g., fault frequency < 1 Hz) [24-26]; however, the actual fault frequency is typically greater than 1 Hz. Lu et al.  developed an embedded system based on multi-scale noise-tuned SR and applied it to train signal detection. Lei et al.  proposed an adaptive signal processing mechanism and applied it to fault detection of missing or broken teeth. Qiao et al.  researched an unsaturated SR to overcome the output saturation problem of classic bi-stable SR (CBSR). This application effectively improved the signal-to-noise ratio (SNR). Lei et al.  investigated under-damped SR on the synergistic effect of the damping factor and vibration signal. It was shown to obtain both a higher amplitude and larger SNR output. Shi et al.  proposed multiple fault frequency detection to the background of color noise. Lu et al.  investigated an SR mechanism for the Woods-Saxon potential function. This could detect weak signals at different noise levels. Lu et al.  proposed a tri-stable SR model that could further amplify the SNR in low SNR signal detection settings.
All of the aforementioned SR studies are considered to be short-term memory systems as they do not consider historical information. A time-delay feedback SR uses historical information of the original system and applies it to the negative feedback to form a new system. This may improve the detection of weak periodic signals; however, few studies have been conducted using this method. Two examples are Lu et al.  and Shi et al. . Lu et al.  used a time-delay CBSR weak fault detection method for the purpose of analyzing bearings; Shi et al.  investigated the SR of a three-stable system with a time-delay feedback. Although these two methods considered time-delay, both were directed to bi-stable or tri-stable nonlinear systems with a single potential structure. Therefore, it is impossible to form a perfect nonlinear system structure and match complex vibration signals. The synergy between time-delay feedback, system parameters, and vibration signals is also not considered which can affect the enhancement ability of SR. Furthermore, the global optimization ability of genetic algorithms, time-delay feedback, system parameters, and vibration signals is synergistic. Therefore, these should also be included to achieve the best SNR output using the most optimal parameter settings.
This paper presents a periodic non-sinusoidal time-delay SR (PNTSR) model and applies it to bearing fault detection. The PNTSR model has a richer potential model structure with coefficients that can be adjusted to better match with the complex vibration signals. The time-delay is introduced on the basis of a traditional SR system, allowing the model to incorporate the influence of historical factors. The second part of this paper analyzes the characteristics of the PNTSR model and discusses the effective potential function (EPF), steady-state probability density function (SPDF), and SNR. The third part describes the weak fault detection strategy and simulation for the PNTSR model. The paper then verifies the PNTSR method using experimental and engineering data and concludes with a summary of the findings.
2. PNTSR introduction
2.1. Model of PNTSR
SR is the synergy between noise and weak signals in a nonlinear system. The appropriate amount of noise and weak signals are used to drive the particles to move periodically into the potential wells on both sides to detect weak signals. CBSR methods are based on the bi-stable potential model:
where and are non-zero values. By adjusting and , different potential model shapes can be obtained. However, this article introduces a new potential model :
Fig. 1. The shape of depends on the parameter. Assuming , is plotted for a) –0.6, b) 0, c) 0.4, and d) 0.6
Fig. 2. The peak value of depends on : a) 0.1, b) 0.25, c) 0.5, d) 0.75, where 0.4
The potential has periodicity but belongs to a non-sinusoidal function, and the period is 2 (Fig. 1, Fig. 2). As you can see from Fig. 1 that the parameter determines the roundness of the barrier and the width of the potential well. The smaller parameter is, the smoother the barrier will be, and the potential well becomes wider. Fig. 2 illustrates that the parameter controls the peak of the barrier, the larger , the higher the barrier. Therefore, adjusting the size of the parameters and determines the shape of the potential structure. Comparing with the CBSR system, the periodic non-sinusoidal potential model has multiple potential wells, when the noise is too large, multiple potential wells can absorb part of the noise energy and reduce noise interference; the periodic non-sinusoidal potential model can achieve independent adjustment of the well width and the barrier height, and a richer potential model structure can be formed by adjusting the parameters.
2.2. Impact of the time-delay feedback on EPF, SPDF, and SNR
Time-delay feedback SR is based on CBSR, which introduces historical information from the original short-term memory system to the long-term memory system, and the classical bi-stable delay system is:
is a extracted signal, is signal strength, and represents the phase. represents the feedback strength, represents the time-delay. stands for noise intensity, stands for external noise. We consider a delay period non-sinusoidal potential system driven by a weak periodic signal. Substituting Eq. (2) into Eq. (3) gives:
Eq. (4) is a non-Markov that needs to be simplified to a Markov process. The Fucker Plane is as follows :
where represents the conditional average drift rate to satisfy:
is the zero-order approximate Markov transition probability density :
Eq. (10) brings in Eq. (6) and we get:
Further derivation can get an Equivalent Langevin equation:
Without considering the periodic signal, the EPF can be derived as:
The effect of different and on the EPF is shown in Fig. 3.
Fig. 3. Different and for EPF
Fig. 3(a) shows that as increases, the potential wells on both sides of the zero point become lower, the barrier is relatively lower, and the round trip of the particles between the potential wells becomes easier. In Fig. 3(b), with the increase of , the potential wells on both sides of the zero point rise, and it is easier for the particles to transition from one side well to the other. Comparing Fig. 3(a) and 3(b), it can be noted that the is more sensitive to the adjustment EPF than the . The time delay varies from 1 to 10, and the EPF changes significantly, while the feedback strength requires only a small parameter variation range. With further calculation, we obtain the SPDF expression:
where is the normalized constant. The SPDF is shown in Fig. 4.
Fig. 4. Different and for SPDF
Different values of have no effect on the peak of the SPDF (Fig. 4(a)). As increases, the SPDF is more concentrated. Figs. 4(b) and 4(a) have similar properties. As increases, the peak of SPDF does not change. Comparing the two figures, the effect of the feedback strength is more sensitive than the effect of the time-delay time on the SPDF. Subsequently the power spectral density can be derived as :
where and is the power spectrum of the useful signal and external noise, respectively. and are as follows:
Subsequently, the output SNR can be expressed as . Substituting and gives:
As shown in Fig. 5, as the feedback strength increases, the SNR tends to increase first and then decrease. Since noise intensity at the detection of a weak fault signal is fixed, the optimum SNR should be obtained by changing the delay term.
3. Detection strategy and simulation for the PNTSR
3.1. Detection strategy for the PNTSR method
Because the realization condition of SR is under the adiabatic approximation theory, it satisfies the small parameter requirement. Therefore, the frequency shift variable scale to change signal scale, and then use the fourth-order Runge-Kutta method and get the output. SNR is the evaluation index. The definition of SNR is:
where represent the data length, and and correspond to the amplitude of the extracted signal and the external noise. The genetic algorithm is used to adjust the parameters , , and to match the complex vibration signals to obtain the optimal output. Therefore, we get the PNTSR fault detection method flow:
(1) Data preprocessing. The signal demodulation methods such as band pass filtering and envelope extraction are used to obtain the driving signal. The small parameter limitation is solved by the frequency shift variable scale.
(2) Parameter initialization. Initialize the potential system parameters, time extension and feedback strength, and use SNR as the output evaluation index.
(3) Parameter optimization. Using the global optimization ability of the genetic algorithm, the system parameters, time-delay feedback and complex vibration signals work together to obtain the optimal values of system parameters, time-delay, and feedback strength.
(4) Output calculation. The result output is calculated by Eq. (14), and then the SNR is calculated by Eq. (15) to ensure the best combination of parameters.
(5) Signal post processing. The optimal parameters are substituted into the potential model to obtain the optimal filtering, which is the spectrum, to realize the identification and diagnosis of fault characteristics.
Fig. 5. Effect of time-delay and feedback strength on SNR
The specific flow chart is as follows (Fig. 6).
3.2. Simulation verification of PNTSR method
To verify the PNTSR method, the PNTSR method is analyzed by simulating the rolling bearing fault signal. The simulation signal is generated by the following formula:
where represents the signal strength, represents the reaction signal attenuation, which control pulse period. 79 Hz is the fault frequency. Moreover, the PNTSR method is contrast CBSR to emphasize that PNTSR is more effective.
The simulated clean signal and the noisy signal are shown in Fig. 7(a) and Fig. 7(b), respectively, and the corresponding spectrogram is seen in Fig. 7(c). Fig. 7(d) shows the resulting envelope spectrum. The spectrum and envelope extraction reveal that the useful signal is completely hidden, and fault signal cannot be detected.
Fig. 6. Weak signal detection strategy for PNTSR
Fig. 7. Bearing simulation data: a) clean data; b) noisy signal; c) spectrogram; d) envelope spectrogram
Fig. 8. a) PNTSR processing signal, b) spectrogram, c) CBSR processing signal, d) spectrogram
In order to extract the simulated rolling bearing fault signal, the simulated fault signal is applied to the CBSR and the PNTSR. Here, the distance between the amplitude of surrounding noise and the amplitude of simulated data is represented by the recognition degree. Comparing Figs. 8(b) and 8(d), the recognition degree of PNTSR method is 0.153, and the recognition degree of CBSR method is 0.01271. In Fig. 8(b), it is apparent that the primary frequency is the fault characteristic frequency and the other noises are smaller. In Fig. 8(d), the signal is more seriously affected by noise. Thus, we can conclude that the PNTSR has better extraction results.
4. Experimental verification of PNTSR method
In the next step, we use the proposed method to detect bearing inner ring with slight wear faults in the laboratory. The experimental platform is a comprehensive bench for mechanical equipment (Fig. 9). The rotational frequency is 40 Hz, the experimental sampling frequency is 5120 Hz, and the bearing inner ring frequency is 217.28 Hz. The processed bearing fault signal is displayed in Fig. 10. In Fig. 10(b), the fault frequency is difficult to recognize due to interference from multiple high frequency and low frequency noise. However, in Fig. 10(c), although the useful signal frequency of 216 Hz, and seriously interfered by low frequency noise.
Fig. 9. Bearing fault simulation platform
We apply the proposed method to the extraction fault feature of the bearing. Fig. 11(b) shows the spectrum output result. The fault frequency is the dominant frequency and the noise component is less affected. The recognition degree of the PNTSR method is 0.21. In Fig. 11(d), the CBSR output shows that the recognition degree of the faulty constant frequency is 0.04807. Therefore, we can conclude that the PNTSR is further better effective than CBSR in the bearing signal diagnosis.
Fig. 10. a) Bearing signal, b) spectrogram, c) envelope spectrogram
Fig. 11. a) PNTSR processing signal, b) spectrogram, c) CBSR processing signal, d) spectrogram
5. Engineering verification of PNTSR method
In view of the good results, we apply the proposed method to engineering practice. The experimental object used is the rolling bearing of a steel mill and its inner ring fault characteristics are identified (Fig. 12). The bearing inner ring useful signal frequency is 24 Hz. The processed data is shown in Fig. 13. Because noise is large, the signal contour in the time domain graph is very inconspicuous. In Fig. 13(b), the extracted signal frequency of 24.3 Hz is hidden in the noise and cannot be discerned. Therefore, applying the PNTSR method processing signal can be seen in Fig. 14(a) and (b), which indicates that the signal profile in the time domain graph is relatively clear. In Fig. 14(b), the extracted signal frequency of 24.32 HZ, other noise interference is small, and the recognition degree of PNTSR method is 0.325. The collected data are applied to the CBSR method, the processed signal can be obtained, as shown in Figs. 14(c) and (d). The CBSR method has a large noise interference and the recognition degree is only 0.023 (Fig. 14d). Therefore, we can conclude that the proposed method not only has good signal extraction effects in simulation and experiments, but is also effective in practical engineering.
Fig. 12. Fault detection of the bearing inner ring of the steel mill.
Fig. 13. a) Engineering verification bearing signal, b) spectrogram, c) envelope spectrogram
Fig. 14. a) PNTSR processing signal, b) spectrogram, c) CBSR processing signa, d) spectrogram
In this paper, the PNTSR system was investigated and a method for diagnosing weak faults was proposed. Following are the conclusions from this paper:
1) The introduced periodic non-sinusoidal potential function consists of two structural parameters and . These can be adjusted to obtain a variety of potential models. Notably, transitions are more abundant.
2) The PNTSR system is a multiple potential well negative feedback system. Because of the presence of multiple potential wells, this system can absorb excess noise energy and reduce noise interference. An additional long memory feedback term was introduced to this system. Under suitable conditions, the feedback item can improve the effect of the periodic weak signal detection by adding the historical information to the current output.
3) The system parameters, time-delay, feedback strength, and complex vibration signals work synergistically while using the global optimization ability of the genetic algorithm. As a result, the SR effect is optimized.
4) The feasibility of the PNTSR method is verified by simulation, experiments, and rolling mill gearbox. In contrast to the CBSR output results, PNTSR is shown to be more effective on fault feature extraction with better recognition and less noise frequency interference.
This work was supported by the National Natural Science Foundation of China (NSFC) under Grant 51805275, 51865045, 51565046. And Young teachers’ scientific research ability improvement plan of Beijing University of Civil Engineering and Architecture, the project number is X21053.
- Sawalhi N., Randall R. B. Vibration response of spalled rolling element bearings: observations, simulations and signal processing techniques to track the spall size. Mechanical Systems and Signal Processing, Vol. 25, Issue 3, 2011, p. 846-870. [Publisher]
- Chen Y. Y., Zhen Z. M., Yu H. H., et al. Application of fault tree analysis and fuzzy neural networks to fault diagnosis in the internet of things (loT) for aquaculture. Sensors, Vol. 17, Issue 12, 2017, p. 153. [Publisher]
- Zhang G., Yi T., Zhang T. Q., et al. A multiscale noise tuning stochastic resonance for fault diagnosis in rolling element bearings. Chinese Journal of Physics, Vol. 56, Issue 1, 2018, p. 145-157. [Publisher]
- Zhu D. C., Zhang Y. X., Liu S. Y., et al. Adaptive combined HOEO based fault feature extraction method for rolling element bearing under variable speed condition. Journal of Mechanical Science and Technology, Vol. 32, Issue 10, 2018, p. 4589-4599. [Publisher]
- Wang Z. J., Wang J. Y., Kou Y. F., et al. Weak fault diagnosis of wind turbine gearboxes based on MED-LMD. Entropy, Vol. 19, Issue 6, 2017, p. 277. [Publisher]
- Mohanty S., Gupta K. K., Raju K. S. Hurst based vibro-acoustic feature extraction of bearing using EMD and VMD. Measurement, Vol. 117, 2018, p. 200-220. [Publisher]
- Vernekar K., Kumar H., Gangadharan K. V. Engine gearbox fault diagnosis using empirical mode decomposition method and Naïve Bayes algorithm. Sādhanā, Vol. 42, Issue 7, 2017, p. 1-11. [Publisher]
- He S. L., Liu Y. K., Chen J. L., et al. Wavelet transform based on inner product for fault diagnosis of rotating machinery. Mechanical Systems and Signal Processing, Vol. 70, Issue 2016, 2016, p. 1-35. [Publisher]
- Bessous N., Zou S. E., Bentrah W., et al. Diagnosis of bearing defects in induction motors using discrete wavelet transform. International Journal of System Assurance Engineering and Management, Vol. 9, Issue 4, 2018, p. 1-9. [Search CrossRef]
- Yang H. G., Lin H. B., Ding K. Sliding window denoising K-Singular value decomposition and its application on rolling bearing impact fault diagnosis. Journal of Sound Vibration, Vol. 421, 2018, p. 205-219. [Publisher]
- Yang B. Y., Liu R. N., Chen X. F. Fault diagnosis for a wind turbine generator bearing via sparse representation and Shift-Invariant K-SVD. IEEE Transactions on Industrial Informatics, Vol. 13, Issue 3, 2017, p. 1321-1331. [Publisher]
- Qiao Z. J., Lei Y. G., Lin J., et al. An adaptive unsaturated bi-stable stochastic resonance method and its application in mechanical fault diagnosis. Mechanical Systems and Signal Processing, Vol. 84, 2017, p. 731-746. [Publisher]
- Deng W., Liu H. L., Xu J. J. An improved quantum-inspired differential evolution algorithm for deep belief network. IEEE Transactions on Instrumentation and Measurement, Vol. 69, Issue 10, 2020, p. 7319-7327. [Publisher]
- Zhao H. M., Zheng J. J., Deng W. Semi-supervised broad learning system based on manifold regularization and broad network. IEEE Transactions on Circuits and Systems I-Regular Papers, Vol. 67, Issue 3, 2020, p. 983-994. [Publisher]
- Zhao H. M., Liu H. L., Xu J. J. Performance prediction using high-order differential mathematical morphology gradient spectrum entropy and extreme learning machine. IEEE Transactions on Instrumentation and Measurement, Vol. 69, Issue 7, 2020, p. 4165-4172. [Publisher]
- Benzi R., Parisi G., Sutera A., et al. A theory of stochastic resonance in climatic change. SIAM Journal on Applied Mathematics, Vol. 43, Issue 3, 1983, p. 565-578. [Publisher]
- Lorito S., Schmitt D., Consolini G., et al. Stochastic resonance in a bi-stable geodynamo model. Astronomische Nachrichten, Vol. 326, Issues 3-4, 2005, p. 227-230. [Publisher]
- Alibegov M. Stochastic resonance, the rayleigh test, and identification of the 25-day periodicity in the solar activity. Astronomy Letters, Vol. 22, Issue 4, 1996, p. 564-572. [Search CrossRef]
- Kaur D., Filonenko I., Mourokh L., et al. Stochastic resonance in a proton pumping Complex I of mitochondria membranes. Scientific Reports, Vol. 7, Issue 1, 2017, p. 12405. [Publisher]
- Qiao Z. J., Lei Y. G., Lin J., et al. Stochastic resonance subject to multiplicative and additive noise: The influence of potential asymmetries. Physical Review E, Vol. 94, Issues 5-1, 2016, p. 52214. [Publisher]
- Yu H. T., Galán R. F., Wang, et al. Stochastic resonance, coherence resonance, and spike timing reliability of Hodgkin–Huxley neurons with ion-channel noise. Physica A Statistical Mechanics and Its Applications, Vol. 471, 2017, p. 263-275. [Publisher]
- Lu S. L., He Q. B., Zhang H. B., et al. Rotating machine fault diagnosis through enhanced stochastic resonance by full-wave signal construction. Mechanical Systems and Signal Processing, Vol. 85, 2017, p. 82-97. [Publisher]
- Liu X. L., Liu H., Yang J., et al. Improving the bearing fault diagnosis efficiency by the adaptive stochastic resonance in a new nonlinear system. Mechanical Systems and Signal Processing, Vol. 96, 2017, p. 58-76. [Publisher]
- Leng Y. G., Leng Y. S., Wang T. Y., et al. Numerical analysis and engineering application of large parameter stochastic resonance. Journal of Sound and Vibration, Vol. 292, Issue 3, 2006, p. 788-801. [Publisher]
- Tan J. Y., Chen X. F., Wang J. Y., et al. Study of frequency-shifted and re-scaling stochastic resonance and its application to fault diagnosis. Mechanical Systems and Signal Processing, Vol. 23, Issue 3, 2009, p. 811-822. [Publisher]
- Cui Y., Zhao J., Guo T. T., et al. Bearing fault diagnosis based on scale-transformation stochastic resonance. Journal of Beijing Information Science and Technology University, Vol. 8916, Issue 1, 2015, p. 36. [Search CrossRef]
- Lu S. L., He Q. B., Hu F., et al. Sequential multiscale noise tuning stochastic resonance for train bearing fault diagnosis in an embedded system. IEEE Transactions on Instrumentation and Measurement, Vol. 63, Issue 1, 2013, p. 106-116. [Publisher]
- Lei Y. G., Han D., Lin J., et al. Planetary gearbox fault diagnosis using an adaptive stochastic resonance method. Mechanical Systems and Signal Processing, Vol. 38, Issue 1, 2013, p. 113-124. [Publisher]
- Lei Y. G., Qiao Z. J., Xu X. F., et al. An underdamped stochastic resonance method with stable-state matching for incipient fault diagnosis of rolling element bearings. Mechanical Systems and Signal Processing, Vol. 94, 2017, p. 148-164. [Publisher]
- Shi P. M., Ding X. J., Han D. Y. Study on multi-frequency weak signal detection method based on stochastic resonance tuning by multi-scale noise. Measurement, Vol. 47, Issue 1, 2014, p. 540-546. [Publisher]
- Lu S. L., He Q. B., Hu F. Stochastic resonance with Woods-Saxon potential for rolling element bearing fault diagnosis. Mechanical Systems and Signal Processing, Vol. 45, Issue 2, 2014, p. 488-503. [Publisher]
- Lu S. L., He Q. B., Zhang H. B., et al. Note: Signal amplification and filtering with a tristable stochastic resonance cantilever. Review of Scientific Instruments, Vol. 84, Issue 2, 2013, p. 026110. [Publisher]
- Dong C. M., Lu C. D., Wang C. J. The effects of time delay on stochastic resonance in a bistable system with correlated noises. Journal of Statistical Physics, Vol. 137, Issue 4, 2009, p. 625-638. [Publisher]
- Shi P. M., Yuan D. Z., Han D. Y., et al. Stochastic resonance in a time-delayed feedback tristable system and its application in fault diagnosis. Journal of Sound and Vibration, Vol. 424, 2018, p. 1-14. [Publisher]
- Djuidjé, Kenmoé G., Ngouongo Y. J. W., Kofané T. C. Effect of the potential shape on the stochastic resonance processes. Journal of Statistical Physics, Vol. 161, Issue 2, 2015, p. 475-485. [Publisher]
- Frank T. D. Delay Fokker-Planck equations, Novikov’s theorem, and Boltzmann distributions as small delay approximations. Physical Review E, Vol. 72, Issue 1, 2005, p. 011112. [Publisher]
- Frank T. D. Delay Fokker-Planck equations, perturbation theory, and data analysis for nonlinear stochastic systems with time delays. Physical Review E, Vol. 71, Issue 3, 2005, p. 031106. [Publisher]
- Mcnamara B., Wiesenfeld K. Theory of stochastic resonance. Physical Review A, Vol. 39, Issue 9, 1989, p. 4854-4869. [Publisher]