Published: 22 June 2021

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
2Baotou Vocational and Technical College, Baotou, China
3, 1School of Mechanical Engineering, University of Science and Technology Inner Mongolia, Baotou, China
Corresponding Author:
Zhixing Li
Views 123
Reads 59
Downloads 863

Abstract

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.

Periodic non-sinusoidal time-delay stochastic resonance weak fault diagnosis method and its application

Highlights

  • 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

1. Introduction

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 [3]. Early fault signals in key parts (e.g., bearings) are often difficult to detect in noisy settings and can lead to inaccurate fault diagnosis [4]. Therefore, for early fault detection, the ability to extract weak fault characteristics within a noisy background is important [5].

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 [12]. 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 [16] 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. [27] developed an embedded system based on multi-scale noise-tuned SR and applied it to train signal detection. Lei et al. [28] proposed an adaptive signal processing mechanism and applied it to fault detection of missing or broken teeth. Qiao et al. [12] 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. [29] 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. [30] proposed multiple fault frequency detection to the background of color noise. Lu et al. [31] investigated an SR mechanism for the Woods-Saxon potential function. This could detect weak signals at different noise levels. Lu et al. [32] 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. [33] and Shi et al. [34]. Lu et al. [33] used a time-delay CBSR weak fault detection method for the purpose of analyzing bearings; Shi et al. [34] 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:

1
Ucx=12ax2+14bx4,

where a and b are non-zero values. By adjusting a and b, different potential model shapes can be obtained. However, this article introduces a new potential model [35]:

2
Ux=U01-r21-cosx1+r2+2rcosx.

Fig. 1The shape of Ux depends on the parameterr. Assuming U0=1, Ux is plotted for a) r= –0.6, b) r= 0, c) r= 0.4, and d) r= 0.6

The shape of Ux depends on the parameterr. Assuming U0=1, Ux is plotted for a) r= –0.6, b) r= 0, c) r= 0.4, and d) r= 0.6
The shape of Ux depends on the parameterr. Assuming U0=1, Ux is plotted for a) r= –0.6, b) r= 0, c) r= 0.4, and d) r= 0.6
The shape of Ux depends on the parameterr. Assuming U0=1, Ux is plotted for a) r= –0.6, b) r= 0, c) r= 0.4, and d) r= 0.6
The shape of Ux depends on the parameterr. Assuming U0=1, Ux is plotted for a) r= –0.6, b) r= 0, c) r= 0.4, and d) r= 0.6

Fig. 2The peak value of Ux depends on U0: a) U0= 0.1, b) U0= 0.25, c) U0= 0.5, d) U0= 0.75, where r= 0.4

The peak value of Ux depends on U0:  a) U0= 0.1, b) U0= 0.25, c) U0= 0.5, d) U0= 0.75, where r= 0.4
The peak value of Ux depends on U0:  a) U0= 0.1, b) U0= 0.25, c) U0= 0.5, d) U0= 0.75, where r= 0.4
The peak value of Ux depends on U0:  a) U0= 0.1, b) U0= 0.25, c) U0= 0.5, d) U0= 0.75, where r= 0.4
The peak value of Ux depends on U0:  a) U0= 0.1, b) U0= 0.25, c) U0= 0.5, d) U0= 0.75, where r= 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 r determines the roundness of the barrier and the width of the potential well. The smaller parameter r is, the smoother the barrier will be, and the potential well becomes wider. Fig. 2 illustrates that the parameter U0 controls the peak of the barrier, the larger U0, the higher the barrier. Therefore, adjusting the size of the parameters r and U0 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:

3
dxtdt=xt-x3t+λxt-τ+AcosΩt+2Dξt.

Acos(Ωt)is a extracted signal, A is signal strength, and Ω represents the phase. λ represents the feedback strength, τ represents the time-delay. D stands for noise intensity, ξ(t) 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:

4
dxtdt=U01-r2sinxt+2rsinxt1-cosxt1+r2+2rcosxt2
+λxt-τ+AcosΩt+2Dξt.

Eq. (4) is a non-Markov that needs to be simplified to a Markov process. The Fucker Plane is as follows [36]:

5
px,tt=heffpx,tx+D2px,tx2,

where heff represents the conditional average drift rate to satisfy:

6
heff=badxτhx,xτpxτ,t-τ|x,t,

where:

7
xτ=xt-τ,
8
hx,xτ=U01-r2sinxt+2rsinxt1-cosxt1+r2+2rcosxt2+λxτ+AcosΩt,
9
hx=U01-r2sinxt+2rsinxt1-cosxt1+r2+2rcosxt2+λx+AcosΩt.

pxτ,t-τ|x,t is the zero-order approximate Markov transition probability density [37]:

10
pxτ,t-τ|x,t=14πDτexp-xτ-x-hxτ24Dτ.

Eq. (10) brings in Eq. (6) and we get:

11
heff=1+λτU01-r2sinxt+2rsinxt1-cosxt1+r2+2rcosxt2
+λ1+λτ+1+λτAcosΩt.

Further derivation can get an Equivalent Langevin equation:

12
dxtdt=U01-r2sinxt+2rsinxt1-cosxt1+r2+2rcosxt2+λx+AcosΩt +λτU01-r2sinxt+2rsinxt1-cosxt1+r2+2rcosxt2+λx+AcosΩt+2Dξt.

Without considering the periodic signal, the EPF can be derived as:

13
Ueffx=1+λτU01-r1-cosx1+r2+2rcosx+12λ1+λτx2.

The effect of different λ and τ on the EPF is shown in Fig. 3.

Fig. 3Different λ and τ for EPF

Different λ and τ for EPF
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:

14
pst=Nexp-UeffD,

where N is the normalized constant. The SPDF is shown in Fig. 4.

Fig. 4Different λ and τ for SPDF

Different λ and τ for SPDF
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 [38]:

15
Sω=S1ω+S2ω=π1+λM22N2+Ω2δΩ-ω+δΩ+ω
+1-M22N2+Ω221+λNN2+ω2,

where S1(ω) and S2(ω) is the power spectrum of the useful signal and external noise, respectively. N and M are as follows:

16
N=21+λ1+λτπexp-1+λ21+λτ4D,
17
M=AD1+λτ1+λN.

Subsequently, the output SNR can be expressed as S1(ω)/S2(ω). Substituting S1(ω) and S2(ω) gives:

18
SNR=S1ωS2ω=2A24D21+λ21+λτ3exp-1+λ21+λτ4D.

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:

19
SNR=10log10Ad2i=1N/2Ai2-Ad2,

where N represent the data length, and Ad2 and i=1N/2Ai2-Ad2 correspond to the amplitude of the extracted signal and the external noise. The genetic algorithm is used to adjust the parameters U0, r, λ 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. 5Effect of time-delay τ and feedback strength λ on SNR

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:

20
St=exp-dt-ntTd2Asin2πft,

where A represents the signal strength, d represents the reaction signal attenuation, nt which control pulse period. f= 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. 6Weak signal detection strategy for PNTSR

Weak signal detection strategy for PNTSR

Fig. 7Bearing simulation data: a) clean data; b) noisy signal; c) spectrogram; d) envelope spectrogram

Bearing simulation data: a) clean data; b) noisy signal; c) spectrogram; d) envelope spectrogram

Fig. 8a) PNTSR processing signal, b) spectrogram, c) CBSR processing signal, d) spectrogram

a) PNTSR processing signal, b) spectrogram, c) CBSR processing signal, d) spectrogram
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. 9Bearing fault simulation platform

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. 10a) Bearing signal, b) spectrogram, c) envelope spectrogram

a) Bearing signal, b) spectrogram, c) envelope spectrogram

Fig. 11a) PNTSR processing signal, b) spectrogram, c) CBSR processing signal, d) spectrogram

a) PNTSR processing signal, b) spectrogram, c) CBSR processing signal, d) spectrogram
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. 12Fault detection of the bearing inner ring of the steel mill.

Fault detection of the bearing inner ring of the steel mill.

Fig. 13a) Engineering verification bearing signal, b) spectrogram, c) envelope spectrogram

a) Engineering verification bearing signal, b) spectrogram, c) envelope spectrogram

Fig. 14a) PNTSR processing signal, b) spectrogram, c) CBSR processing signa, d) spectrogram

a) PNTSR processing signal, b) spectrogram, c) CBSR processing signa, d) spectrogram
a) PNTSR processing signal, b) spectrogram, c) CBSR processing signa, d) spectrogram

6. Conclusions

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 U0 and r. 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.

References

  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Mcnamara B., Wiesenfeld K. Theory of stochastic resonance. Physical Review A, Vol. 39, Issue 9, 1989, p. 4854-4869.

About this article

Received
02 June 2020
Accepted
24 November 2020
Published
22 June 2021
SUBJECTS
Fault diagnosis based on vibration signal analysis
Keywords
stochastic resonance
periodic non-sinusoidal potential
time-delay feedback
weak signal
fault detection
Acknowledgements

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.