To extract weak faults under strong noise, a method for feature extraction of weak faults with time-delayed feedback mixed potential stochastic resonance (TFMSR) is proposed. This method not only overcomes the saturation characteristics of classical bistable stochastic resonance (CBSR), but also verifies a new potential function model. Based on this model, considering the short memory characteristics of the CBSR method, a method is proposed that can add historical information to the negative feedback process of the stochastic resonance (SR). Through the combination of the above two methods, the weak fault extraction under strong background noise is realized. The article analyzes the effects of the delay term, feedback term, and system parameter on the effect of SR and uses the ant colony algorithm (ACA) to optimize the above parameters. Finally, through simulated and engineering experimental results, it is proved that the proposed method has more advantages than the CBSR method in weak fault feature extraction.
National Natural Science Foundation of China51805275Beijing University of Science and Technology Innovation Talents Fund1. Introduction
With the development of technology, equipment becomes large scale and intelligent with high speed [1, 2]. Failure of any component in the system can jeopardize the operational safety of the entire device. Therefore, early prediction of weak faults in mechanical systems is particularly important for the extraction of the weak fault feature. It is an emerging science to study how to extract weak signals from strong noise. It requires multiple disciplines to work together, so it is also an interdisciplinary [3, 4].
Benzi et al. [5] first proposed SR when explaining the problem of ancient glaciers. It transfers some parts of the noise energy to the weak signal. While enhancing the faint fault characteristics, some noise becomes weakened, thereby improving the signal-to-noise ratio (SNR). Existing noise reduction methods, such as wavelet analysis [6–8], empirical mode decomposition (EMD) [9], variational mode decomposition [10], and dictionary learning [11, 12], can weaken the useful weak signal when filtering the noise. As a signal processing tool, SR has an interesting feature; that is, the noise added into the nonlinear system can excite and amplify the early weak fault features. Utilizing this principle, Yang et al. [13] obtained an amazing result. They recover unknown signals completely submerged in strong noise through a new SR method. Due to this unique characteristic of SR, it has attracted wide attention in the fields of physics [14], medicine [15], and machinery.
Initially, SR theory suffers from the small frequency limitation under adiabatic approximate conditions [16]. Therefore, it cannot be applied in actual engineering. In order to make SR universal [17, 18], many researchers have made their own contributions. Leng et al. [19] studied a recalibration of the frequency SR method to achieve a large parameter signal SR. Tan et al. [20] used frequency shift variable scale processing SR to alleviate the contradiction between the number of sample points and sampling frequency. Later, the CBSR model became the most widely used SR model. Li et al. [21] proposed a Morlet wavelet transform CBSR method that can obtain a finer frequency, which is used to extract wind turbine fault characteristics. Lei et al. [22] proposed an adaptive SR method for mechanical fault diagnosis. He et al. [23] proposed a new multiscale noise adjustment method to improve SR performance. Lu et al. [24] proposed a multiscale noise adjustment SR method, which was proved in the detection of actual sound signals from train bearings. Shi et al. [25] found that feedback parameters, time-delay parameters, and asymmetric parameters all have an effect on SR, and the time-delay parameters are more sensitive. Li and Shi [26] proposed an SR with a bistable constrained potential and applied to weak fault feature extraction of bearings. Nishiguchi and Fujiwara [27] suggested boosting an effect of CBSR and silhouette noise from noise via a field-effect transistor. The above method mainly uses the CBSR method to lay a foundation for the rapid development of SR. Later, Zhao et al. [28] proved the existence of the output saturation characteristics of the CBSR system. Gosak et al. [29] studied the negative effects of output saturation on the system output in an ultradamped oscillator. Zhang et al. [30] studied SR output saturation based on the additive and multiplicative noises. These studies have recognized the effect of output saturation characteristics on SR, but have not proposed an effective solution. When saturation occurs, the amplitude of the output signal will not increase as the amplitude of the input signal increases. This makes it difficult to extract weak fault characteristics. Therefore, this paper proposes a mixed potential SR model to solve the saturation problem of the CBSR model. Further, the CBSR system is also a short-memory system, without considering the effects of time delay and feedback. This makes that the output of the system depend on the output of the previous one. Lv et al. [31, 32] studied the special role of delay in nonlinear systems. Zhuo et al. [33, 34] studied the application of memory and feedback models in fault diagnosis. In summary, adding historical information to the system’s feedback process in a certain form will greatly improve the output signal-to-noise ratio (SNR).
This paper proposes a method for extracting weak fault features based on TFMSR. This method not only overcomes the saturation characteristics of the CBSR, but also adds historical information to the system’s feedback process in a certain form, effectively solving the two major shortcomings of the CBSR method. Through simulated and experimental results, it proves that the proposed method has more advantageous than the CBSR method in weak fault feature extraction.
2. Theoretical Foundation
The CBSR model is expressed as(1)dxdt=−dU1xdx+st+ξt,where st=AcosΩt is a weak periodic signal, A is the amplitude, and ξt is a zero mean Gaussian white noise. The following condition is satisfied:(2)Eξtξt+τ=2Dδt−τ,where D is the noise intensity. The potential function of the CBSR system is expressed as(3)U1x=−12ax2+14bx4,where a and b represent the system parameters. In Figure 1, U1x is the potential function curve of the CBSR model. Without loss of generality, we set a = 1 and b = 1. Random boost is related to the specific nonlinear form of the potential function. By selecting different nonlinear potential functions and adjusting related parameters, the random boost can be effectively enhanced. Therefore, in order to avoid the output saturation characteristic of the CBSR potential function, we propose a mixed potential stochastic resonance model. It is expressed as(4)Ux=exp−x2−lne+14mx4,where m is the system parameter. The mixed potential function model curve U(x) and the CBSR potential function curve U1x are shown in Figure 1. The potential function of the mixed potential model is wider than that of CBSR, which reflects that the proposed method has more advantages in antisaturation characteristics than CBSR [30].
The mixed potential function model curve U(x) and the CBSR potential function curve U1x (a = 1, b = 1).
Considering the effects of the delay term and the feedback term, the Langevin equation is obtained as follows:(5)dxtdt=−dUxdx+Kxt−σ+st+ξtdxtdt=−2xe−x2−mx3+Kxt−σ+AcosΩt+ξt.
The probability density of Fokker–Planck with approximate delay feedback can be expressed as [35](6)∂Px,t∂t=−∂hnPx,t∂x+D∂2Px,t∂x2,where hn is the conditional average drift. hn can be obtained [36] as(7)hn=−1+Kσ2xe−x2+mx3−Kx+1+KσAcosΩt.
Without considering the effects of periodic signals, equation (7) can be deduced as(8)Uex=1+Kσe−x2−lne+m4x4−K2x2.
Compared with the CBSR, the potential well of TFMSR is affected by the system parameter m, the feedback term K, and the delay term σ. As shown in Figure 2, by adjusting the values of m, K, and σ, we can observe that the shape of the potential well is constantly changing.
Curve of potential function with different parameters: (a) when m = 1 and σ = 1.2, the relationship between U, x, and K; (b) detailed relationship diagram between U and K; (c) when m = 1 and K = 0.7, the relationship between U, x, and σ; (d) detailed relationship diagram between U and σ; (e) when K = 0.7 and σ = 1.2, the relationship between U, x, and m; (f) detailed relationship diagram between U and m.
In Figures 2(a) and 2(b), as the feedback intensity K decreases from 1.1 to 0.7, the height of the potential barrier remains unchanged, but the bottom of the potential well moves slowly upward. In Figures 2(c) and 2(d), by adjusting the delay term σ, the depth of the potential well continuously changes. As σ increases from 0.4 to 1.2, the steepness of the potential wall increases. In Figures 2(e) and 2(f), as m decreases from 1 to 0.6, the steepness of the potential wall does not change, but the depth of the potential well increases.
According to the theory of the two-state model [37], the probability density of particles moving between two potential wells can be expressed as(9)R±t=K1+Kσ2πexp−K21+Kσ4Dmexp∓K/mA1+Kσ4DcosΩt.
Under small delay condition, equation (9) can be deduced as(10)λ=K1+Kσ2πexp−K21+Kσ4Dm.
Under adiabatic approximation [38], the power spectral density expression is(11)sω=s1ω+s2ω,where the output power spectral density of the signal can be expressed as(12)s1ω=Kπ2m1+KσAK/mD24λ24λ2+Ω2δω−Ω+δω+Ω.
The power spectral density of the noise can be expressed as(13)s2ω=1−A2K1+Kσ24λ22mD24λ2+Ω24Kλm4λ2+ω2.
The output SNR of the system can be expressed as [39](14)SNR=2A2K24mD21+Kσ3exp−K2+K3σ4mD.
In Figure 3(a), the system parameter m and delay term σ remain unchanged. As the noise intensity D increases, a single peak appears in the output SNR. As analyzed in Section 2, as K decreases, the depth of the potential well is reduced. The energy required for a particle to transition from one potential well to another decreases, the output signal becomes larger, and the effect of SR is obvious. In Figure 3(b), the system parameter m and the feedback intensity K remain unchanged. When σ = 0.4, because σ is small, the particles cannot break through the barrier of the potential and realize the transition from one potential well to another. The effect of SR is not obvious. Increasing σ to 1.2, the particles have enough energy to cross the potential barrier. At this time, the particles oscillate back and forth between the two potential wells. In Figure 3(b), we can see that the effect of SR is obvious. In Figure 3(c), the feedback intensity K and the delay term σ remain unchanged. When the system parameter m = 0.6, the potential wall is too steep, and the particles cannot obtain sufficient acceleration to break through the barrier constraint. Increasing m to 1, the potential well depth decreases, and the potential function structure can form the best match with the periodic force, so as to achieve the best SR effect.
The change trends of SNR versus D: (a) m = 1 and σ = 1.2; (b) m = 1 and K = 0.7; (c) K = 0.7 and σ = 1.2.
3. Proposed Model
Through the above analysis, we can conclude that the effect of SR is most obvious when the delay term, feedback intensity, and system parameters are optimally matched. This article uses ACA to match and optimize the above parameters. The detection strategy of the TFMSR method is as follows:
Signal Preprocessing. First, the acquired signal is demodulated by the Hilbert transform, and then the small parameters (with signal amplitude, noise intensity, and signal frequency all being far less than 1) required for SR are obtained by frequency shift variable scale processing [20].
Parameter Initialization and Optimization. First, initialize the parameters and optimize delay term, feedback intensity, and system parameter. The range of parameter optimization is set to [0, 5].
Calculating the SNR. The optimal match of parameters are inputted into the SR system. The maximum output SNR of the TFMSR method is calculated by the Runge–Kutta equation [40].
Postprocessing. Through the postprocessing of the signal by the TFMSR system, the faint fault features are extracted. Figure 4 shows the detection strategy of the TFMSR method.
TFMSR weak signal detection strategy.
4. Simulation
To verify the effort of the TFMSR method, we decided to conduct a simulation experiment. First, a periodic signal with frequency f = 55 Hz and amplitude A = 1 is simulated. The sampling frequency is 10000 Hz. Figures 5(a)–5(c) show the time domain, the spectrum, and the envelope spectrum of the simulated signal. Since the background noise is too strong and the characteristic frequency is submerged, we cannot get the fault characteristic frequency. Here, we use the TFMSR method to process the simulated signal. The ACA was used to adaptively optimize the system parameters. The obtained parameters are σ = 1.4623, K = 0.3321, and m = 0.9928. The time domain and spectrum obtained by the simulated signal are shown in Figures 6(a) and 6(b). Figure 6(b) shows that the spectral peak value Amax obtained by using the TFMSR method is 0.07659. The value of the second highest point is 0.03546. The difference value ΔA is 0.04113. We propose a new evaluation index to evaluate the effectiveness of the two methods. We use ΔA/Amax to get a dimensionless constant. The value of the dimensionless constant is 0.537.
(a) Simulation signal without noise; (b) simulation signal with noise; (c) spectrum; (d) envelope spectrum.
The results of the TFMSR method: (a) time domain waveform of the simulated signal; (b) spectrum. The results of the CBSR method: (c) time domain waveform of the simulated signal; (d) spectrum.
For comparison, we will use the CBSR method to process the simulated signal. The initial parameters are settled consistent with the TFMSR method. The obtained parameters are a = 0.0369 and b = 0.117. The parameters are substituted into the CBSR system. Figures 6(c) and 6(d) show the time domain and spectrum of the simulated signal. In Figure 6(d), the spectral peak value A1max is 0.04329. The value of the second highest point is 0.03364. The difference value ΔA1 is 0.00965. ΔA1/A1max=0.2229. In summary, the proposed TFMSR method can enhance weak signals effectively.
5. Engineering Experiment Verification
In the field of metallurgy, rolling mills are widely used. Due to the harsh working environment, the fault signals of the components in the equipment are easily overwhelmed by strong noise. Some minor failures, if cannot be found in time, may create serious consequences over time. The hot rolling production site of a thin plate factory is shown in Figure 7. The working environment of the equipment is very bad. Through routine inspections, we found that the gearbox of the rolling mill failed, but it was not possible to determine which component was faulty. The on-site signal acquisition device is the ZonicBook/618E, and the sampling position is shown in Figure 8.
Thin-plate mill hot rolling production site.
Sensor placement location.
Table 1 lists the basic parameters of the gearbox. The acquired signal is processed to obtain the original signal, spectrum, and envelope spectrum as shown in Figures 9(a), 9(b), and 9(c), respectively. Due to the strong background noise, we were unable to clearly distinguish the fault characteristics and the specific location of the damaged components. In order to find the specific location of the fault from the strong background noise, the proposed TFMSR method is used to process the signal. The obtained parameters are σ = 1.0691, K = 0.1048, and m = 0.8175. Figures 10(a) and 10(b) show the time domain waveform and the spectrum obtained by the acquired signal. Figure 10(b) shows that the fault characteristic frequency is 25 Hz. By comparing the characteristic frequency calculated in Table 1, we can see the characteristic frequency of the bearing 1. The characteristic frequency of the inner ring is 24.3 Hz. Therefore, we boldly conclude that the inner ring of the bearing 1 has failed.
Parameters of the gearbox.
Name
Speed (r/min)
Parameters
Characteristic frequency (Hz)
Z1
63
32
33.6
Z2
20.5714
98
33.6
I axis
63
—
1.05
II axis
20.5714
—
0.343
Rolling
Outer
Inner
Bearing 1
63
DP = 610, d = 62, n = 42
4.9
19.7
24.3
Bearing 2
63
DP = 692, d = 118, n = 53
2.9
22.7
32.6
Bearing 3
20.5714
DP = 610, d = 62, n = 42
1.6
6.4
7.9
Bearing 4
20.5714
DP = 692, d = 118, n = 53
0.9
7.4
10.6
Rolling bearing signal: (a) original signal; (b) spectrum; (c) envelope spectrum.
The results of the TFMSR method: (a) time domain waveform of the original signal; (b) spectrum. The results of the CBSR method: (c) time domain waveform of the original signal; (d) spectrum.
Further, we use CBSR method to process the signal. The obtained parameters are a = 0.0708 and b = 0.195. Figures 10(c) and 10(d) show the time domain waveform and the spectrum of the output signal. Figure 10(d) shows that the fault characteristic frequency is 27.6 Hz. By comparing with the data in Table 1, we cannot find the characteristic frequency which is close to 27.6 Hz. This further proves the advancement of the proposed TFMSR method. To finally confirm the fault location, we use the endoscope to check the inside of the gearbox. It was finally confirmed that cracks occurred in the inner ring of the bearing 1, and the crack is shown in Figure 11. This is in complete agreement with our speculation, further validating the usability of the TFMSR method.
Gearbox bearing 1 with crack failure.
From the above analysis, we can conclude that the TFMSR method has a significant detection effect at l-speed single faults. However, the availability of this method has yet to be verified in the face of weak fault detection under high-speed and multifault conditions. This will also be our next research work.
6. Conclusion
Considering the saturation characteristics of the CBSR model, this paper explores a new potential function model. Because of the short memory characteristics of the CBSR method, we add historical information to the feedback process of the SR. Through the combination of the above two methods, the weak fault extraction under strong background noise is realized. Based on the theoretical analysis and experimental results, the main conclusions are as follows:
A potential function model with antisaturation characteristics is proposed, and the characteristics of the potential function model are analyzed.
Considering the short memory characteristics of the CBSR method, a delay term and a feedback term are introduced on the basis of the proposed potential function model, and the output SNR is theoretically derived. The analysis of the effect of parameters on the output SNR is completed.
The proposed method is used to simulate bearing fault extraction. Simulated results show that the proposed method effectively extracts fault features and is superior to the CBSR method.
The proposed TFMSR method was used to verify the failure of the rolling mill gearbox. The experimental results show that the TFMSR method has better antisaturation characteristics and feedback characteristics than the CBSR method in extracting weak fault characteristics.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (51805275) and in part by the Beijing University of Science and Technology Innovation Talents Fund.
RandallR. B.AntoniJ.Rolling element bearing diagnostics-a tutorial201125248552010.1016/j.ymssp.2010.07.0172-s2.0-78649600770BaydarN.ChenQ.BallA.KrugerU.Detection of incipient tooth defect in helical gears using multivariate statistics200115230332110.1006/mssp.2000.13152-s2.0-0035271127HuN.ChenM.QinG.XiaL.PanZ.FengZ.Extended stochastic resonance (SR) and its applications in weak mechanical signal processing20094445046110.1007/s11465-009-0072-32-s2.0-70350347121YanS.MaB.ZhengC.ChenJ.An optimal lubrication oil replacement method based on selected oil field data20197921109211810.1109/access.2019.29274262-s2.0-85073888967BenziR.SuteraA.VulpianaA.The mech-anism of stochastic resonance1981141145345710.1088/0305-4470/14/11/0062-s2.0-20644470598YanR.GaoR. X.ChenX.Wavelets for fault diagnosis of rotary machines: a review with applications20149611510.1016/j.sigpro.2013.04.0152-s2.0-84887433963LinJ.QuL.Feature extraction based on morlet wavelet and its application for mechanical fault diagnosis2000234113514810.1006/jsvi.2000.28642-s2.0-0034206997ChenJ.LiZ.PanJ.Wavelet transform based on inner product in fault diagnosis of rotating machinery: a review201670-7113510.1016/j.ymssp.2015.08.0232-s2.0-84961055987LeiY.LinJ.HeZ.ZuoM. J.A review on empirical mode decomposition in fault diagnosis of rotating machinery2013351-210812610.1016/j.ymssp.2012.09.0152-s2.0-84870404381ZhangS.WangY.HeS.JiangZ.Bearing fault diagnosis based on variational mode decomposition and total variation denoising201627710.1088/0957-0233/27/7/0751012-s2.0-84976322329CempelC.Generalized singular value decomposition in multidimensional condition monitoring of machines-A proposal of comparative diagnostics200923370171110.1016/j.ymssp.2008.07.0042-s2.0-58049145306ThanhN. T.NiamsupP.PhatV. N.Finite-time stability of singular nonlinear switched time-delay systems: a singular value decomposition approach201735483502351810.1016/j.jfranklin.2017.02.0362-s2.0-85015371650HuangD.YangJ.ZhouD.SanjuánM. A. F.LiuH.Recovering an unknown signal completely submerged in strong noise by a new stochastic resonance method20196615616610.1016/j.cnsns.2018.06.0112-s2.0-85048792634GammaitoniL.The long run of the stochastic resonance idea2016494510.1088/1751-8113/49/45/4510052-s2.0-84992751900FloresA.ManillaS.HuidobroN.Stochastic resonance in the synaptic transmission between hair cells and vestibular primary afferents in development201632241642910.1016/j.neuroscience.2016.02.0512-s2.0-84960406296TangJ.ShiB.LiZ.Asymmetric delay feedback stochastic resonance detection method based on prior knowledge particle swarm optimization20185652104211810.1016/j.cjph.2018.08.0192-s2.0-85054176114KlameckiB. E.Use of stochastic resonance for enhancement of low-level vibration signal components200519222323710.1016/j.ymssp.2004.03.0062-s2.0-4344561562LiJ.PanX.A nonlinear monostable filter for bipolar pulse signal detection20072131223123210.1016/j.ymssp.2006.03.0122-s2.0-33845442457LengY. G.LengY. S.WangT. Y.GuoY.Numerical analysis and engineering application of large parameter stochastic resonance20062923-578880110.1016/j.jsv.2005.09.0402-s2.0-33644596185TanJ.ChenX.WangJ.Study of frequency-shifted and re-scaling stochastic resonance and its application to fault diagnosis200923381182210.1016/j.ymssp.2008.07.0112-s2.0-58049183063LiJ.ChenX.DuZ.FangZ.HeZ.A new noise-controlled second-order enhanced stochastic resonance method with its application in wind turbine drivetrain fault diagnosis20136071910.1016/j.renene.2013.04.0052-s2.0-84877340574LeiY.HanD.LinJ.HeZ.Planetary gearbox fault diagnosis using an adaptive stochastic resonance method201338111312410.1016/j.ymssp.2012.06.0212-s2.0-84876936237HeQ.WangJ.LiuY.DaiD.KongF.Multiscale noise tuning of stochastic resonance for enhanced fault diagnosis in rotating machines20122844345710.1016/j.ymssp.2011.11.0212-s2.0-84857362290LuS.HeQ.HuF.KongF.Sequential multiscale noise tuning stochastic resonance for train bearing fault diagnosis in an embedded system201463110611610.1109/tim.2013.22752412-s2.0-84890550776ShiP.XiaH.HanD.FuR.Dynamical complexity and stochastic resonance in an asymmetry bistable system with time delay201755113314110.1016/j.cjph.2016.10.0132-s2.0-85010638232LiZ.ShiB.Fault diagnosis of rotating machinery based on stochastic resonance with a bistable confining potential2018201812845250910.1155/2018/84525092-s2.0-85061650189NishiguchiK.FujiwaraA.Detecting signals buried in noise via nanowire transistors using stochastic resonance20121011919310810.1063/1.47669462-s2.0-84869027336ZhaoW.WangJ.WangL.The unsaturated bistable stochastic resonance system20132333311710.1063/1.48163972-s2.0-84884925758GosakM.PercM.KraljS.Stochastic resonance in a locally excited system of bistable oscillators201180451952810.1140/epjb/e2011-10573-82-s2.0-79955828469ZhangG.HuD.ZhangT.Stochastic resonance in unsaturated piecewise nonlinear bistable system under multiplicative and additive noise for bearing fault diagnosis20197584355844810.1109/access.2019.29141382-s2.0-85065871364LvM. L.ShenG.WangH. L.Is the high-frequency signal necessary for the resonance in the delayed system2015321151710.1088/0256-307x/32/1/0105012-s2.0-84920836244YangJ. H.SanjuánM. A. F.LiuH. G.Signal generation and enhancement in a delayed system2015221-31158116810.1016/j.cnsns.2014.08.0052-s2.0-85027942994ZhuoP.ZhuY.WuW.ShuJ.XiaT.Real-time Fault diagnosis for gas turbine blade based on output-hidden feedback elman neural network201823S19510210.1007/s12204-018-2028-42-s2.0-85059169931XiaT.SongY.ZhengY.PanE.XiL.An ensemble framework based on convolutional bi-directional LSTM with multiple time windows for remaining useful life estimation202011510318210.1016/j.compind.2019.103182GuillouzicS.L’HeureuxI.LongtinA.Small delay approximation of stochastic delay differential equations19995943970398210.1103/physreve.59.39702-s2.0-0000680355FrankT. D.Delay Fokker-Planck equations, perturbation theory, and data analysis for nonlinear stochastic systems with time delays20057133110610.1103/physreve.71.0311062-s2.0-41349118922McnamaraB.WiesenfeldK.Theory of stochastic resonance19893994854486910.1103/physreva.39.48542-s2.0-0000183328GuillouzicS.L’HeureuxI.LongtinA.Small delay approximation of stochastic delay differential equations199959439703982NguyenS. D.NguyenQ. H.ChoiS.-B.A hybrid clustering based fuzzy structure for vibration control-part 2: an application to semi-active vehicle seat-suspension system201556-5728830110.1016/j.ymssp.2014.10.0192-s2.0-84919721347ØksendalB.2003Berlin, GermanySpringer