# Analysis of Asymmetric Piecewise Linear Stochastic Resonance Signal Processing Model Based on Genetic Algorithm

Authors: Lina He; Chuan Jiang
Journal: Complexity
 (2020)

---

## Abstract

The stochastic resonance system has the advantage of making the noise energy transfer to the signal energy. Because the existing stochastic resonance system model has the problem of poor performance, an asymmetric piecewise linear stochastic resonance system model is proposed, and the parameters of the model are optimized by a genetic algorithm. The signal-to-noise ratio formula of the model is derived and analyzed, and the theoretical basis for better performance of the model is given. The influence of the asymmetric coefficient on system performance is studied, which provides guidance for the selection of initial optimization range when a genetic algorithm is used. At the same time, the formula is verified and analyzed by numerical simulation, and the correctness of the formula is proved. Finally, the model is applied to bearing fault detection, and an adaptive genetic algorithm is used to optimize the parameters of the system. The results show that the model has an excellent detection effect, which proves that the model has great potential in fault detection.


---

## Body

## 1. Introduction

When the nonlinear system interacts with the signal to be measured and the nonsignal to be measured (noise), resonance will occur. At this time, when the noise intensity increases, the noise will not inundate the signal but will increase the signal-to-noise ratio, and the noise energy will transfer to the signal. This shows that stochastic resonance is a powerful method to extract a weak signal from strong noise [1, 2]. However, a large number of studies have shown that the classical SR theory has obvious detection advantages only in the case of small parameters, i.e., adiabatic approximation. In engineering application, most target signals are not small parameter signals, which greatly limits the application in engineering. In recent years, a series of achievements have been made in the resonance phenomenon of large parameter signals, such as the use of secondary sampling, single sideband modulation, frequency-domain information exchange, and other methods, which can make large parameter signals produce resonance phenomenon [3]. In reference [4], it is known that the piecewise linear model has better performance than the classical model, and in reference [5], the asymmetric model can obtain better performance than the symmetric model by adjusting the asymmetric factor. Zhang et al. applied a stochastic resonance system for bearing fault detection [6, 7]. Zhang et al. found that the performance of the system is also affected by the system parameters [8], and the selection of parameters will directly affect the quality of the system. Zhang and He obtained better system parameters through an adaptive genetic algorithm and applied them to bearing fault detection [9].At present, some scholars are committed to exploring the effect of asymmetry on system performance [2, 4]. Some scholars have devoted themselves to the study of piecewise linear systems. Both of them can achieve better system performance [10]. However, there is no research on the combination of them to obtain better performance. Based on this, an asymmetric linear piecewise bistable model is proposed. The analytical expression and signal-to-noise ratio of the model are derived and compared with the symmetrical bistable piecewise linear stochastic resonance system and the continuous bistable system. At the same time, in order to obtain better performance in the application of bearing fault detection, a genetic algorithm is used to obtain better system parameters. In the second part of the paper, the concrete formula derivation is given. It provides a theoretical basis for the simulation and experiment. In the third part, numerical simulation is carried out to verify the correctness of the formula derivation. In the fourth part, the genetic algorithm is proposed to optimize the parameters, which can improve the system performance to the greatest extent. Finally, it is applied to fault detection. The results show that this method can effectively extract fault feature frequency and has better performance than the symmetric model.


## 2. Theoretical Analysis of Asymmetric Piecewise Linear Model

### 2.1. Asymmetric Piecewise Linear Model

The expression of potential functionUx is as follows:(1)Ux=−k3rk1−k2x+rk1,x<−rk2,k3rk2x,−rk2≤x<0,−k3k2x,0≤x<k2,−k3k1−k2x−k1k2≤x.Among them,r is asymmetric factor, k1, k2 and k3 are system parameters and are all greater than 0. Under the static condition, the system has two potential wells and one barrier. The bottom of the two wells is x+=k2 and x−=−rk2, respectively, and the barrier height is ΔU=c.Whenr=1 is a symmetric piecewise linear model, other cases are asymmetric linear systems. It can be seen from Figure 1 that the width of a potential well varies with the asymmetric factor.Figure 1
Potential functionUx.


### 2.2. System Response and Signal-to-Noise Ratio

The steady state of the system is+xm,−xn and rxm=xn. Let W±t be the probability of time t system in bistability +xm,−xn and define n±t as the probability of transition from steady state +xm,−xn at time t. Because of the asymmetry of the two potential wells, the transition probability of the two potential wells is not the same, that is, n+t≠n−t. When the asymmetry factor is r=1, the equation holds. According to the adiabatic approximation theory [2], the escape rate of the asymmetric linear bistable system is as follows:(2)ƛ−−1=1D∫−rk1−rk2e−1/Dk3/rk1−k2x+rk1dx∫−rk20e1/Dk3/rk2xdx,ƛ+−1=1D∫−k1−k2e−1/Dk3/k1−k2x+k1dx∫−k20e1/Dk3/k2xdx,where D is the noise intensity, r is asymmetric factor, k1, k2 and k3 are system parameters. The transition rate n±t is generally considered to have an exponential form. Under the action of periodic signal st=Acosw0t, it is expanded by Taylor series, as shown in the following equation:(3)n+=ƛ+1+AxmDcosω0t+12AxmD2cos2ω0t+⋯,n−=ƛ−1−AxnDcosω0t+12AxnD2cos2ω0t−⋯.Under the assumption of adiabatic approximation, the probability equation of the model can be established according to (2) and (3):(4)dw−tdt=−ƛ−W−t+ƛ+W+t,dw+tdt=−ƛ+W+t+ƛ−W−t.The equations of solution (4) can be obtained as follows:(5)W−t=121ƛ−+e−ƛ−+ƛ+t,W+t=121ƛ+−e−ƛ−+ƛ+t.The probability distribution function obtained from formula (5) is shown in the following equation:(6)w−t=W˙−t=−ƛ−+ƛ+2e−ƛ−+ƛ+,w+t=W˙+t=ƛ−+ƛ+2e−ƛ−+ƛ+.The frequency-domain probability distribution function can be obtained by Fourier transform of formula (6)(7)wω=w−ω+w+ω=ƛ−+ƛ+ƛ−+ƛ+2+ω2,φω=arctanωƛ−+ƛ+.SinceA is small, the response of the bistable system can be expressed as follows:(8)xt=x¯Dcosω0t+ϕ¯D,where(9)x¯D=1+r2Axm2Dƛ−+ƛ+ƛ−+ƛ+2+ω02,Φ¯D=arctanω0ƛ−+ƛ+.The power spectrum of noise can be written as follows:(10)SNω=1−Axm22Dƛ−+ƛ+2ƛ−+ƛ+2+ω02⋅ƛ−+ƛ+21+r2xm2ƛ−+ƛ+2+ω02.The output power spectrum is obtained by Fourier transform of the autocorrelation function, as follows:(11)Sω=π2x¯2Dδω−ω0+δω+ω0+SNω.Then SNR is as follows:(12)SNR=∫−∞+∞Sω−SNωdwSNω0=1+r2πA2xm2ƛ−+ƛ+4D21−1/2Axm/D2ƛ−+ƛ+2/ƛ−+ƛ+2+ω02.For the system,xm=b and obtained from formula (2):(13)ƛ−=k32r2k1−k2k2e−k3/D,ƛ+=k32k1−k2k2e−k3/D.If the higher order term of the denominator in formula (12) is not considered, the SNR can be approximated as follows:(14)SNR=πk21+r22k32A24D2k1−k2r2e−k3/D.Whenr=1, the SNR formula of the symmetrical piecewise linear model shown in [5] can be obtained. According to reference [2], the formula of signal-to-noise ratio of continuous bistable system is as follows:(15)SNR=2u2A24D2e−u2/4D,where u is the system parameter, the barrier height is ΔU=u2/4, and the well bottom is in x=±u, u=2. Select the system parameter k1=1, k2=1, k3=1. The height of the potential well is the same as the position of the potential well, and then the asymmetry factor is changed to make a comparative study. When the input signal and noise are identical, as shown in Figure 2, it is the change of SNR with the noise intensity D. It can be seen from Figure 2 that with the increase of the noise intensity D of the piecewise linear system and the classical continuous bistable system, the SNR also increases rapidly to the maximum value, but the noise continues to increase and decreases. The asymmetric piecewise linear system also has the same phenomenon, but in the process of noise intensity from small to large, the system is always larger than the other two systems. It can be seen that the higher SNR can be obtained for the asymmetric piecewise linear system with the same system parameters.Figure 2
Variation of SNR withD in the three systems.


## 2.1. Asymmetric Piecewise Linear Model

The expression of potential functionUx is as follows:(1)Ux=−k3rk1−k2x+rk1,x<−rk2,k3rk2x,−rk2≤x<0,−k3k2x,0≤x<k2,−k3k1−k2x−k1k2≤x.Among them,r is asymmetric factor, k1, k2 and k3 are system parameters and are all greater than 0. Under the static condition, the system has two potential wells and one barrier. The bottom of the two wells is x+=k2 and x−=−rk2, respectively, and the barrier height is ΔU=c.Whenr=1 is a symmetric piecewise linear model, other cases are asymmetric linear systems. It can be seen from Figure 1 that the width of a potential well varies with the asymmetric factor.Figure 1
Potential functionUx.


## 2.2. System Response and Signal-to-Noise Ratio

The steady state of the system is+xm,−xn and rxm=xn. Let W±t be the probability of time t system in bistability +xm,−xn and define n±t as the probability of transition from steady state +xm,−xn at time t. Because of the asymmetry of the two potential wells, the transition probability of the two potential wells is not the same, that is, n+t≠n−t. When the asymmetry factor is r=1, the equation holds. According to the adiabatic approximation theory [2], the escape rate of the asymmetric linear bistable system is as follows:(2)ƛ−−1=1D∫−rk1−rk2e−1/Dk3/rk1−k2x+rk1dx∫−rk20e1/Dk3/rk2xdx,ƛ+−1=1D∫−k1−k2e−1/Dk3/k1−k2x+k1dx∫−k20e1/Dk3/k2xdx,where D is the noise intensity, r is asymmetric factor, k1, k2 and k3 are system parameters. The transition rate n±t is generally considered to have an exponential form. Under the action of periodic signal st=Acosw0t, it is expanded by Taylor series, as shown in the following equation:(3)n+=ƛ+1+AxmDcosω0t+12AxmD2cos2ω0t+⋯,n−=ƛ−1−AxnDcosω0t+12AxnD2cos2ω0t−⋯.Under the assumption of adiabatic approximation, the probability equation of the model can be established according to (2) and (3):(4)dw−tdt=−ƛ−W−t+ƛ+W+t,dw+tdt=−ƛ+W+t+ƛ−W−t.The equations of solution (4) can be obtained as follows:(5)W−t=121ƛ−+e−ƛ−+ƛ+t,W+t=121ƛ+−e−ƛ−+ƛ+t.The probability distribution function obtained from formula (5) is shown in the following equation:(6)w−t=W˙−t=−ƛ−+ƛ+2e−ƛ−+ƛ+,w+t=W˙+t=ƛ−+ƛ+2e−ƛ−+ƛ+.The frequency-domain probability distribution function can be obtained by Fourier transform of formula (6)(7)wω=w−ω+w+ω=ƛ−+ƛ+ƛ−+ƛ+2+ω2,φω=arctanωƛ−+ƛ+.SinceA is small, the response of the bistable system can be expressed as follows:(8)xt=x¯Dcosω0t+ϕ¯D,where(9)x¯D=1+r2Axm2Dƛ−+ƛ+ƛ−+ƛ+2+ω02,Φ¯D=arctanω0ƛ−+ƛ+.The power spectrum of noise can be written as follows:(10)SNω=1−Axm22Dƛ−+ƛ+2ƛ−+ƛ+2+ω02⋅ƛ−+ƛ+21+r2xm2ƛ−+ƛ+2+ω02.The output power spectrum is obtained by Fourier transform of the autocorrelation function, as follows:(11)Sω=π2x¯2Dδω−ω0+δω+ω0+SNω.Then SNR is as follows:(12)SNR=∫−∞+∞Sω−SNωdwSNω0=1+r2πA2xm2ƛ−+ƛ+4D21−1/2Axm/D2ƛ−+ƛ+2/ƛ−+ƛ+2+ω02.For the system,xm=b and obtained from formula (2):(13)ƛ−=k32r2k1−k2k2e−k3/D,ƛ+=k32k1−k2k2e−k3/D.If the higher order term of the denominator in formula (12) is not considered, the SNR can be approximated as follows:(14)SNR=πk21+r22k32A24D2k1−k2r2e−k3/D.Whenr=1, the SNR formula of the symmetrical piecewise linear model shown in [5] can be obtained. According to reference [2], the formula of signal-to-noise ratio of continuous bistable system is as follows:(15)SNR=2u2A24D2e−u2/4D,where u is the system parameter, the barrier height is ΔU=u2/4, and the well bottom is in x=±u, u=2. Select the system parameter k1=1, k2=1, k3=1. The height of the potential well is the same as the position of the potential well, and then the asymmetry factor is changed to make a comparative study. When the input signal and noise are identical, as shown in Figure 2, it is the change of SNR with the noise intensity D. It can be seen from Figure 2 that with the increase of the noise intensity D of the piecewise linear system and the classical continuous bistable system, the SNR also increases rapidly to the maximum value, but the noise continues to increase and decreases. The asymmetric piecewise linear system also has the same phenomenon, but in the process of noise intensity from small to large, the system is always larger than the other two systems. It can be seen that the higher SNR can be obtained for the asymmetric piecewise linear system with the same system parameters.Figure 2
Variation of SNR withD in the three systems.


## 3. Numerical Simulation of Asymmetric Linear Model

In order to verify the above results, the symmetrical piecewise linear model has been compared with the continuous bistable system in reference [5], and the performance of the piecewise linear model is better. Therefore, only asymmetric and symmetric piecewise linear models are compared. According to formula (1), the system model can be written as follows:(16)x˙t=−dUxdx+Ht,where Ht=Acos2πft+εt, εt for noise. Its mean value is zero, and autocorrelation function can be expressed as εtεt+τ=2Dδt−τ. D is the noise intensity, τ is the delay time, f is the input signal frequency, and A is its amplitude. The fourth-order Runge Kutta of (1) is simulated [4], and the resonance behavior of r=1 and r≠1 is compared and explored.Figure3 is the input signal frequency f=0.01, sampling frequency, and amplitude noise intensity fs=5,A=0.1,D=4, system parameter k1=2,k2=1,k3=0.25. Figures 3(b) and 3(d) show the amplitude frequency characteristics of the input signal and the output signal, respectively. It can be seen that the spectrum energy concentrates on the low-frequency component. While Figure 4 is under the condition that other conditions remain unchanged, making the asymmetric factor r=1.5, it can be seen from Figure 4(a) that the reduction effect is better than Figure 3(c), and from the frequency domain, it can be seen that the energy concentrated on the input signal frequency is higher.Figure 3
Time-domain and frequency-domain graphics of input and output signals.
(a)(b)(c)(d)Figure 4
Time-domain and frequency-domain graphics of input and output signals.
(a)(b)Figure5 is a comparison of the average signal-to-noise ratio gain [10] of the asymmetric piecewise linear model and the symmetric piecewise linear model under the same other conditions as Figure 4. It can be seen from the figure that the average signal-to-noise ratio gain of the asymmetric model is always higher than that of the symmetric model with the increase of noise intensity, which is consistent with the conclusion of the formula derivation. Figure 5 shows that the performance of the asymmetric system is not superior to that of the symmetric model because the selection of asymmetric factors is not optimal. Figure 6 is a three-dimensional graph of noise intensity D, asymmetric factor r an average signal-to-noise ratio gain MSNRI under the same conditions as Figure 5. It can be seen that r=2 is not the optimal case. As you can see, the performance improvement is not very high. This is because the asymmetric coefficient is not very suitable. It is because of this that the adaptive algorithm mentioned below is needed to get the appropriate system parameters.Figure 5
Average SNR gain.Figure 6
Three-dimensional graph of MSNRI varying withD and R.


## 4. Adaptive Genetic Algorithm

The model proposed in this paper has four parameters, namelyk1, k2,k3,r and the dimension of the genetic algorithm is set as four dimensions [11–13]. The key of the algorithm is to transform the solution into the chromosome needed in a genetic algorithm. There are a variety of “chromosome” transformation methods, which can be divided into real number coding (parameter optimization problem) and integer coding (shortest path problem) according to requirements. The genetic algorithm has a wide range of applications. Usually, some changes will be made according to actual needs. In this study, the range of parameter optimization is determined by theoretical analysis, that is, the range of initial gene in genetic algorithm. By replacing the signal-to-noise ratio with fitness function, the purpose of improvement is achieved. The genetic algorithm can be well applied to the stochastic resonance system. The specific process is shown in Figure 7. Specific steps of the proposed algorithm:Figure 7
Flow chart of genetic algorithm.The fitness function of the paper is the output signal-to-noise ratio. Signal-to-noise ratio (SNR) is a common measure in signal processing research. The higher the SNR, the better the fitness. Compared with correlation and average signal-to-noise ratio, the genetic algorithm has lower time complexity and little difference in effect.Initialize population size and iteration: Random selection of individuals to build an initial population.The fitness function constructed in step (1) is used to calculate the fitness of all individuals. And keep the best individuals to the next generation.The individual genes in the population were crossed, and when the crossing process met the variation conditions, the cross was executed (5). The crossover operator of the following formula is used:(17)X1′=λ1X1+1−λ2X2,X2′=λ1X2+1−λ2X1.Variation, that is to say, the offspring produced genes that the parents did not have. The construction of the mutation operator is as follows:(18)X′=X+Δ.Generate offspring and replace any random individual in the offspring with the optimal solution individual in step (3). According to whether the termination condition reaches the maximum number of iterations, the algorithm's branch flow is determined.Adding optimization parameters to the genetic algorithm will not exponentially increase the time complexity like an ordinary iterative algorithm. The time complexity of the iterative algorithm isoNn, N is the number of iterations, and n is the number of parameters. The time complexity of GA is oN∗nn. N is the genetic algebra, n is the number of parameters.


## 5. Bearing Fault Detection

### 5.1. Fault Characteristics Extraction of the Bearing Type 6205-2RS JEM SKF

#### 5.1.1. Characteristic Frequency

The bearing model is 6205-2RS JEM SKF deep groove ball bearing. The main parameters are shown in Table1. Because the condition of a small parameter is not satisfied, the resonance phenomenon is produced by using the method of second sampling.Table 1
Main parameters of bearing.
Inner diameter (cm)Outer ring diameter (cm)Ball straight (cm)Thickness (cm)Stanza (cm)Number of ball bearings (one)2.50015.19990.79401.50013.9049Sampling frequencyfs=12000Hz, number of sampling points N=10000. Here, 5 Hz is selected as the second sampling frequency. The calculation formula of the characteristic frequency is as follows:(19)fBPFI=nrfr21+D1D2cosα,fBPFO=nrfr21−D1D2cosα,nr represents the number of rolling elements, D1 represents the diameter, D2 represents the bearing, the rotation frequency is fr, the contact angle is α, fBPFI and fBPFO are the characteristic frequencies of the inner and outer rings of the bearing, respectively. Substituting the data in the table into equation (19), it is known that the characteristic frequency of the outer ring fault is 107.28 Hz, and that of the inner ring fault is 162.11 Hz.


#### 5.1.2. Outer Ring Fault Detection

Figure8 shows the input signal of bearing fault to be detected and the output signal of the symmetrical piecewise linear system. It can be seen that the fault signal is completely covered by other high-frequency noises. After passing the SR system, the fault signal is detected. The parameters of this check are all found through the adaptive parameter optimization of the genetic algorithm. The population number of genetic algorithms is 200, the genetic algebra is 200, the crossover probability is set to 0.4, and the mutation probability is set to 0.2. (In the actual natural environment, the crossover probability is much greater than the mutation probability, but in order to jump out of the local optimum, the mutation probability is set to be larger here.)Figure 8
Time and frequency-domain diagrams of input and output of outer circle fault detection systema=0.422000206845763,b=0.01,c=2.645191139445148,r=1.
(a)(b)(c)(d)Figure9 is an output signal through an asymmetrical piecewise linear system. Compared with Figure 8, it can be seen that the accuracy of amplitude value and inspection at fault frequency is higher than that of the symmetrical segmented system. The time-domain waveforms of the stochastic resonance system are compared. It can be seen that the burr of time-domain waveform of an asymmetric system is obviously less than that of a symmetric system. It is consistent with higher amplitude at fault frequency of frequency-domain waveform.Figure 9
Time and frequency-domain diagrams of the output of outer circle fault detection systema=1.00619166102941,b=0.11,r=1.6,c=0.284256838355532.
(a)(b)


#### 5.1.3. Inner Ring Fault Detection

In order to prove that there is still a lot of room to improve the performance of the system, the number of population and the number of iterations are doubled when a genetic algorithm is used for parameter optimization. Figure10 shows the signal to be detected for the inner ring fault and the output signal after the piecewise linear system model. It can be found that the fault signal has been detected, but after a longer time of parameter optimization, the system performance has not been significantly improved.Figure 10
Time-domain and frequency-domain diagrams of the input and output of inner-loop fault detection systema=1.263827550905825,b=1.226069,r=1,c=0.253128458753468.
(a)(b)(c)(d)Figure11 shows the output waveform of the signal to be detected through the asymmetric piecewise linear model. It can be found that the performance is greatly improved after the time of parameter optimization is increased. It is the same as the previous conclusion and proves the superiority of the system.Figure 11
Time and frequency-domain diagrams of the output of outer circle fault detection systema=0.938548936195951,b=0.11,c=2.11048612862940,r=1.6.
(a)(b)


### 5.2. Fault Diagnosis of ID-25/30 Bearing Health Test Bench

In order to further verify the practicability of the system, it is applied to another group of experimental platform for verification. The physical figure of the ID-25/30 bearing health test bench is shown in Figure12. The calculation method of the characteristic frequency signal of bearing fault is the same as the previous one. The frequency of the inner ring and outer ring can be obtained from the formula. The inner raceway frequency of the bearing can be calculated theoretically to be fInner=117.14Hz.Figure 12
The entity diagram of the ID-25/30 test bench.Figure13 shows the fault signal to be tested collected by the platform. It can be seen that the fault model is completely submerged by noise and cannot be detected. Figure 14 is the output graph of the piecewise linear system. It can be seen that the fault signal is detected. This part of the parameters is also optimized by genetic algorithm, which proves the practical value of the system again. Figure 15 shows the output of the asymmetrical piecewise linear system. Figure 15 has a higher amplitude and the frequency doubling has been checked, indicating that the detection effect performance is very good.Figure 13
Signal to be detected.
(a)(b)Figure 14
Output signal of the piecewise linear system.
(a)(b)Figure 15
Output signal of the asymmetrical piecewise linear system.
(a)(b)


## 5.1. Fault Characteristics Extraction of the Bearing Type 6205-2RS JEM SKF

### 5.1.1. Characteristic Frequency

The bearing model is 6205-2RS JEM SKF deep groove ball bearing. The main parameters are shown in Table1. Because the condition of a small parameter is not satisfied, the resonance phenomenon is produced by using the method of second sampling.Table 1
Main parameters of bearing.
Inner diameter (cm)Outer ring diameter (cm)Ball straight (cm)Thickness (cm)Stanza (cm)Number of ball bearings (one)2.50015.19990.79401.50013.9049Sampling frequencyfs=12000Hz, number of sampling points N=10000. Here, 5 Hz is selected as the second sampling frequency. The calculation formula of the characteristic frequency is as follows:(19)fBPFI=nrfr21+D1D2cosα,fBPFO=nrfr21−D1D2cosα,nr represents the number of rolling elements, D1 represents the diameter, D2 represents the bearing, the rotation frequency is fr, the contact angle is α, fBPFI and fBPFO are the characteristic frequencies of the inner and outer rings of the bearing, respectively. Substituting the data in the table into equation (19), it is known that the characteristic frequency of the outer ring fault is 107.28 Hz, and that of the inner ring fault is 162.11 Hz.


### 5.1.2. Outer Ring Fault Detection

Figure8 shows the input signal of bearing fault to be detected and the output signal of the symmetrical piecewise linear system. It can be seen that the fault signal is completely covered by other high-frequency noises. After passing the SR system, the fault signal is detected. The parameters of this check are all found through the adaptive parameter optimization of the genetic algorithm. The population number of genetic algorithms is 200, the genetic algebra is 200, the crossover probability is set to 0.4, and the mutation probability is set to 0.2. (In the actual natural environment, the crossover probability is much greater than the mutation probability, but in order to jump out of the local optimum, the mutation probability is set to be larger here.)Figure 8
Time and frequency-domain diagrams of input and output of outer circle fault detection systema=0.422000206845763,b=0.01,c=2.645191139445148,r=1.
(a)(b)(c)(d)Figure9 is an output signal through an asymmetrical piecewise linear system. Compared with Figure 8, it can be seen that the accuracy of amplitude value and inspection at fault frequency is higher than that of the symmetrical segmented system. The time-domain waveforms of the stochastic resonance system are compared. It can be seen that the burr of time-domain waveform of an asymmetric system is obviously less than that of a symmetric system. It is consistent with higher amplitude at fault frequency of frequency-domain waveform.Figure 9
Time and frequency-domain diagrams of the output of outer circle fault detection systema=1.00619166102941,b=0.11,r=1.6,c=0.284256838355532.
(a)(b)


### 5.1.3. Inner Ring Fault Detection

In order to prove that there is still a lot of room to improve the performance of the system, the number of population and the number of iterations are doubled when a genetic algorithm is used for parameter optimization. Figure10 shows the signal to be detected for the inner ring fault and the output signal after the piecewise linear system model. It can be found that the fault signal has been detected, but after a longer time of parameter optimization, the system performance has not been significantly improved.Figure 10
Time-domain and frequency-domain diagrams of the input and output of inner-loop fault detection systema=1.263827550905825,b=1.226069,r=1,c=0.253128458753468.
(a)(b)(c)(d)Figure11 shows the output waveform of the signal to be detected through the asymmetric piecewise linear model. It can be found that the performance is greatly improved after the time of parameter optimization is increased. It is the same as the previous conclusion and proves the superiority of the system.Figure 11
Time and frequency-domain diagrams of the output of outer circle fault detection systema=0.938548936195951,b=0.11,c=2.11048612862940,r=1.6.
(a)(b)


## 5.1.1. Characteristic Frequency

The bearing model is 6205-2RS JEM SKF deep groove ball bearing. The main parameters are shown in Table1. Because the condition of a small parameter is not satisfied, the resonance phenomenon is produced by using the method of second sampling.Table 1
Main parameters of bearing.
Inner diameter (cm)Outer ring diameter (cm)Ball straight (cm)Thickness (cm)Stanza (cm)Number of ball bearings (one)2.50015.19990.79401.50013.9049Sampling frequencyfs=12000Hz, number of sampling points N=10000. Here, 5 Hz is selected as the second sampling frequency. The calculation formula of the characteristic frequency is as follows:(19)fBPFI=nrfr21+D1D2cosα,fBPFO=nrfr21−D1D2cosα,nr represents the number of rolling elements, D1 represents the diameter, D2 represents the bearing, the rotation frequency is fr, the contact angle is α, fBPFI and fBPFO are the characteristic frequencies of the inner and outer rings of the bearing, respectively. Substituting the data in the table into equation (19), it is known that the characteristic frequency of the outer ring fault is 107.28 Hz, and that of the inner ring fault is 162.11 Hz.


## 5.1.2. Outer Ring Fault Detection

Figure8 shows the input signal of bearing fault to be detected and the output signal of the symmetrical piecewise linear system. It can be seen that the fault signal is completely covered by other high-frequency noises. After passing the SR system, the fault signal is detected. The parameters of this check are all found through the adaptive parameter optimization of the genetic algorithm. The population number of genetic algorithms is 200, the genetic algebra is 200, the crossover probability is set to 0.4, and the mutation probability is set to 0.2. (In the actual natural environment, the crossover probability is much greater than the mutation probability, but in order to jump out of the local optimum, the mutation probability is set to be larger here.)Figure 8
Time and frequency-domain diagrams of input and output of outer circle fault detection systema=0.422000206845763,b=0.01,c=2.645191139445148,r=1.
(a)(b)(c)(d)Figure9 is an output signal through an asymmetrical piecewise linear system. Compared with Figure 8, it can be seen that the accuracy of amplitude value and inspection at fault frequency is higher than that of the symmetrical segmented system. The time-domain waveforms of the stochastic resonance system are compared. It can be seen that the burr of time-domain waveform of an asymmetric system is obviously less than that of a symmetric system. It is consistent with higher amplitude at fault frequency of frequency-domain waveform.Figure 9
Time and frequency-domain diagrams of the output of outer circle fault detection systema=1.00619166102941,b=0.11,r=1.6,c=0.284256838355532.
(a)(b)


## 5.1.3. Inner Ring Fault Detection

In order to prove that there is still a lot of room to improve the performance of the system, the number of population and the number of iterations are doubled when a genetic algorithm is used for parameter optimization. Figure10 shows the signal to be detected for the inner ring fault and the output signal after the piecewise linear system model. It can be found that the fault signal has been detected, but after a longer time of parameter optimization, the system performance has not been significantly improved.Figure 10
Time-domain and frequency-domain diagrams of the input and output of inner-loop fault detection systema=1.263827550905825,b=1.226069,r=1,c=0.253128458753468.
(a)(b)(c)(d)Figure11 shows the output waveform of the signal to be detected through the asymmetric piecewise linear model. It can be found that the performance is greatly improved after the time of parameter optimization is increased. It is the same as the previous conclusion and proves the superiority of the system.Figure 11
Time and frequency-domain diagrams of the output of outer circle fault detection systema=0.938548936195951,b=0.11,c=2.11048612862940,r=1.6.
(a)(b)


## 5.2. Fault Diagnosis of ID-25/30 Bearing Health Test Bench

In order to further verify the practicability of the system, it is applied to another group of experimental platform for verification. The physical figure of the ID-25/30 bearing health test bench is shown in Figure12. The calculation method of the characteristic frequency signal of bearing fault is the same as the previous one. The frequency of the inner ring and outer ring can be obtained from the formula. The inner raceway frequency of the bearing can be calculated theoretically to be fInner=117.14Hz.Figure 12
The entity diagram of the ID-25/30 test bench.Figure13 shows the fault signal to be tested collected by the platform. It can be seen that the fault model is completely submerged by noise and cannot be detected. Figure 14 is the output graph of the piecewise linear system. It can be seen that the fault signal is detected. This part of the parameters is also optimized by genetic algorithm, which proves the practical value of the system again. Figure 15 shows the output of the asymmetrical piecewise linear system. Figure 15 has a higher amplitude and the frequency doubling has been checked, indicating that the detection effect performance is very good.Figure 13
Signal to be detected.
(a)(b)Figure 14
Output signal of the piecewise linear system.
(a)(b)Figure 15
Output signal of the asymmetrical piecewise linear system.
(a)(b)


## 6. Conclusion

An unsymmetrical piecewise linear system has good efficiency. This paper discusses the performance of the system from three aspects: formula derivation, numerical simulation, and engineering application. The results show that it has a good performance both in theory and in practice, which is consistent with the formula derivation results. Of course, the performance of a system is closely related to its parameters, and it is difficult to find out the coordination between parameters. Therefore, a genetic algorithm is proposed to find better parameters.The asymmetry mentioned in this paper is the asymmetry of the width of the potential well and the steady state is assumed to be two steady states, which is not explored in the asymmetry of the depth of the potential well and the multisteady state, so we will continue to study in the following work.



---

*Source: 8817814-2020-10-30.xml*