Rolling Bearing Fault Diagnosis Based on Adaptive Multiparameter-Adjusting Bistable Stochastic Resonance

+e weak-signal detection technologies based on stochastic resonance (SR) play important roles in the vibration-based health monitoring and fault diagnosis of rolling bearings, especially at their early-fault stage. Aiming at the parameter-fixed vibration signals in practical engineering, it is feasible to diagnose the potential rolling bearing faults through adaptively adjusting the SR system parameters, as well as other generalized parameters such as the amplitude-transformation coefficient and scale-transformation coefficient. However, extant adaptive adjustment methods focus on the system parameters, while the adjustments of other adjustable parameters have not been fully studied, thus limiting the detection performance of the adaptive SR method. In order to further enhance the detection performance of adaptive SR methods and extend their application in rolling bearing fault diagnosis, an adaptive multiparameter-adjusting SR (AMPASR) method for bistable systems based on particle swarm optimization (PSO) algorithm is proposed in this paper. +is method can produce optimal SR output through adaptively adjusting multiparameters, thus realizing fault feature extraction and further fault diagnosis. Furthermore, the influence of algorithm parameters on the optimization results is discussed, and the optimization results of the Langevin system and the Duffing system are compared. Finally, we propose a weak-signal detection method based on the AMPASR of the Duffing system and employ three diagnosis examples involving inner ring fault, outer ring fault, and rolling element fault diagnoses to demonstrate its feasibility in rolling bearing fault diagnosis.


Introduction
Rolling bearings, as core components of rotating machinery, play a significant role in modern machines such as wind turbines, machine tools, centrifugal pumps, compressors, and motorized spindles. Due to their harsh and complex working environments, faults are easy to occur in the inner ring, outer ring, and rolling elements of rolling bearings, thus resulting in the breakdown of the entire machine, or even a disastrous accident. erefore, it is necessary to monitor the rolling bearings' operating conditions and diagnose their potential faults as soon as possible, thus avoiding potential accidents through timely maintenance.
In the past few decades, various diagnosis technologies based on vibration, acoustics, temperature, and liquid have been proposed for the health monitoring and fault diagnosis of rolling bearings. Among all these diagnosis technologies, the vibration-based diagnosis technologies, such as maximum correlated kurtosis deconvolution (MCKD) [1], spectral kurtosis (SK) [2], empirical mode decomposition (EMD) [3], wavelet transform (WT) [4], and chaos theory [5], have received a lot of attention because the vibration signals with rolling bearing faults contain rich condition information. However, due to tough working environments and remote measurement distance, the fault vibration signals are always weak, especially at the fault's early stage.
Under these conditions, the ability of weak-signal detection is a key index to evaluate the methods' diagnosis performance. Although the aforementioned diagnosis methods can achieve weak-signal detection to some extent through suppressing or canceling the noise embedded in fault vibration signals and highlighting the fault features, their weak-signal detection performances are limited because they inevitably damage the weak fault features submerged in heavy noise background in the denoising process. Compared with these noise cancellation-based fault diagnosis methods, a nonlinear phenomenon called stochastic resonance (SR) leads a type of noise utilization-based fault diagnosis methods, which have intrinsic superiority in weak-signal detection by taking advantage of the noise to enhance the weak fault features through some nonlinear systems [6]. SR was first proposed by Benzi to describe the periodicity associated with the Earth's ice ages in climatology in the 1980s [7]. is interesting nonlinear phenomenon makes it possible for the weak signals to strengthen their intensities by absorbing a certain fraction of the noise energy, thus highlighting the weak signals [8,9]. SR has attracted much attention from the physics and engineering community in the research areas such as dynamical nonlinearity [10][11][12], structure monitoring [13], and energy harvesting [14][15][16]. Moreover, in view of the good performance of SR using noise to enhance periodic signal features, SR-based weaksignal detection and fault diagnosis methods have been investigated and successfully applied in rolling bearing fault diagnosis [17,18]. e performance of an SR-based weak-signal detection method is mainly determined by the form of its nonlinear systems. Research has shown that SR can appear in both monostable systems [19,20] and bistable systems [21][22][23]. e latter ones, which involve the classical overdamped Langevin system and underdamped Duffing system, have much better SR performance than the former systems. erefore, the bistable systems are mostly used as weaksignal detection models in previous studies.
SR describes an optimal synergistic action of the input signal, noise, and nonlinear systems. erefore, aiming at the parameter-fixed noisy signals, such as the measured vibration fault signals of rolling bearings, it is possible, and the only way to adjust the system parameters, to match the signal parameters to produce SR and obtain further feature extractions.
is method to achieve SR is the so-called parameter-adjusting SR method [24]. e tunable system parameters in bistable SR systems contain the potential function parameters and damping ratio (only for the underdamped systems). Furthermore, researchers found that, for the signals with inappropriate amplitude and large frequency, it is necessary to introduce amplitude-transformation coefficient and scale-transformation coefficient, respectively, to transform the amplitude and frequency of the characteristic signal to an appropriate range [25]. Combining the amplitude-transformation coefficient and scale-transformation coefficient, multiparameter-adjusting SR methods are further proposed and the multiparameters adjustment rules are fully studied theoretically [26,27]. According to the qualitative adjustment rules, SR can be achieved by manually adjusting the multiparameters for different noisy signals.
However, the success of manual adjustments depends on the experience of engineers. e adjustment process may be time-consuming and the optimal SR output may not be obtained by engineers without sufficient experiences. is makes it an inefficient and even an unreliable way to realize SR and fault diagnosis in practical engineering. Moreover, the manual adjustment method can only be applied in offline detections. In order to solve these issues, by using some multiparameter optimization algorithms such as particle swarm optimization (PSO) [28] and genetic algorithm (GA) [22,29], some adaptive parameter-adjusting methods, which can realize SR in online conditions with high efficiency and high reliability, have been proposed and studied [30][31][32][33][34][35].
ese adaptive SR methods have been applied in the rolling bearing fault diagnosis successfully.
However, extant adaptive methods mainly study the adjustments of traditional system parameters, while the adjustments of other generalized parameters such as amplitude-transformation coefficient and scale-transformation coefficient have not been fully studied. erefore, the optimal output result may not be achieved especially for those signals with a low signal-to-noise ratio (SNR), thus limiting the detection performance of the adaptive SR methods. In order to further enhance the detection performance of adaptive SR methods and extend their application in fault diagnosis, it is necessary to take the multiparameters into account and propose corresponding adaptive multiparameter-adjusting SR method. Moreover, in most references [22,[32][33][34]36], the adaptive methods were adopted directly to diagnosis mechanical faults without the reliability of the optimization results being verified and the optimization parameters being discussed. Furthermore, the comparison between the adaptive fault diagnosis performances of Langevin system and Duffing system, which has not been studied yet, should be conducted that a better bistable system can be selected in further applications.
In this paper, we propose an adaptive multiparameteradjusting SR (AMPASR) method based on bistable systems. e application of this method on rolling bearing fault diagnosis is also presented. e remainder of this paper is organized as follows. In Section 2, we introduce the multiparameter-adjusting SR for both the Langevin system and the Duffing system. In Section 3, an adaptive multiparameter-adjusting SR method for bistable systems based on PSO algorithm is proposed and the adaptive optimization results are analyzed. In Section 4, we discuss the influence of algorithm parameters on the optimization results and compare the optimization results of the Langevin system and the Duffing system. Section 5 presents the procedure of the weak-signal detection based on AMPASR of the Duffing system, along with three rolling bearing fault diagnosis examples. Conclusion and summary are drawn in Section 6.

2
Shock and Vibration

Multiparameter-Adjusting SR for Bistable Systems
Bistable SR systems have two typical forms. One is described by equation (1) indicating an underdamped Duffing system [31] driven by a weak signal s(t) and a source of noise n(t), and the other is described by equation (2) indicating an overdamped Langevin system [37] whose inertia term is ignored. In both bistable systems, dU(x)/dx � − ax + bx 3 denotes the force of the potential field, whose potential function U(x) � − ax 2 /2 + bx 4 /4 with a and b being positive system parameters. In Duffing system equation (1) k denotes the damping ratio. sn(t) � s(t) + n(t) is defined as the system input signal; in the simulation analyses of this paper, s(t) � A cos(2πf 0 t) represents a harmonic weak signal with amplitude A and frequency f 0 , and ��� 2D √ ξ(t) denotes a source of Gaussian white noise with noise intensity D: In the absence of both s(t) and n(t), the potential function in equations (1) and (2) describes a bistable potential field with two stable equilibrium points at x � ± ��� a/b √ and one unstable equilibrium point at x � 0, as shown in Figure 1. e height of the potential barrier that separates two potential wells is ΔU � a 2 /(4b). is is why equations (1) and (2) are called bistable systems.
e output x(t) of equations (1) and (2) can be understood as the time-related position of a unit-mass Brownian particle moving in the bistable potential field U(x) under the coaction of damping force − k _ x(t) (only for the Duffing system equation (1)), potential field force − dU(x)/dx, driving force s(t) and random noise n(t).
rough theoretical and simulation analyses, previous research has demonstrated that SR can occur in both Duffing system equation (1) and Langevin system equation (2). at is, with the help of nonlinear systems, the noise can produce a positive effect on the weak signal s(t) and enhance its intensity in the output x(t). erefore, SR in bistable systems has potential advantages in weak-signal detection.
Actually, SR presents a synergy between the weak signal, noise, and nonlinear systems. SR can only occur when these three factors are matched. For the measured vibration signal sn(t) with mechanical faults, the signal parameters are fixed. erefore, in order to utilize SR of bistable systems to detect weak fault signals, the only way is to adjust the system parameters to match the vibration signals. e SR achieved through adjusting system parameters is the so-called parameter-adjusting SR. Traditionally, the system parameters to be adjusted involve k, a, b in Duffing system equation (1) and a, b in Langevin system equation (2). However, the adjustments of these system parameters are insufficient to match various noisy signals, especially for those signals with a largefrequency component, an inadequate amplitude, and a low signal-to-noise ratio (SNR). erefore, an amplitude-transformation coefficient ε and a scale-transformation coefficient R are further introduced in equations (1) and (2), which can be rewritten as Here, ε is used to linearly transform the amplitude of the input signal to an appropriate range, and R � Δt ′ /Δt is used to transform the time/frequency scale of the input signal. e scale transformation is simply realized by setting the time step as Δt ′ � R/f s instead of Δt � 1/f s in numerical calculations, where f s indicates the sampling frequency. us, the system output x(t ′ ) in terms of time scale t ′ � Rt can be obtained. Based on Duffing system equation (3) and Langevin system equation (4), a multiparameter-adjusting method [26] has been proposed to achieve SR by jointly adjusting the generalized parameters ε, R and system parameters k (only for Duffing system), a, b. e SR-oriented adjustment rules for each parameter have been studied and summarized. ese parameter adjustment rules have been successfully used in SR-based weak-signal detection and bearing fault diagnosis.
However, the parameter adjustment rules only show how to adjust the separate parameter to achieve SR in bistable systems. It does not reveal how to jointly adjust multiparameters, which is helpful to achieve SR for complicated and low-SNR input signals. Moreover, the existing multiparameter-adjusting method relies on the observation and manual adjustments of the engineer, and the optimal parameter sets and SR output may not be found out. erefore, it is necessary to study an adaptive multiparameter-adjusting SR (AMPASR) method for bistable systems, which can adjust multiparameters adaptively to achieve optimal SR output for different input signals, to fully utilize the advantages of multiparameter-adjusting method for weak-signal detection and rolling bearing fault diagnosis. Figure 1: Bistable potential function. Shock and Vibration 3

AMPASR for Bistable Systems Based on PSO Algorithm
In this section, the particle swarm optimization (PSO) algorithm will be introduced as an adaptive algorithm, and the output SNR for bistable systems will be defined as the objective function. us, an AMPASR method based on PSO algorithm will be proposed to achieve SR adaptively in bistable systems. e optimization results will be further analyzed in this section.

Introduction of PSO Algorithm.
e PSO algorithm is an efficient multiparameter optimization algorithm. In PSO, each particle represents a potential solution to the optimization issue, and it is determined by the fitness value of the objective function. Each particle has a velocity, which determines its flight direction and distance. e velocity is dynamically adjusted by the moving experiences of the particle itself and other particles. e optimal solution will be updated according to the particles' best fitness. e PSO algorithm can be described in mathematics. Assume that there is a population composed of N p particles in a d-dimensional solution space. e position and velocity of the ith particle in the population are e optimal solutions for each particle and the whole population are and gbest � P g � (p g,1 , p g,2 , . . . , p g,d ), respectively. Before the particles find out these two optimal solutions, the particle updates its velocity and position according to where q � 1, 2, 3, . . . , M is the iteration number with M the largest one; w is the inertia weight; c 1 and c 2 are the positive learning factors; r 1 and r 2 are evenly distributed random numbers among [0, 1]. In this paper, c 1 and c 2 are chosen as 2, and w is decided by an adaptive weighting method: where w max � 0.9 and w min � 0.6 are the given maximal and minimal values of w; f means the particle's current fitness value of the objective function, f avg and f min denote the average value and minimum value of the objective function of all particles. With regard to the multiparameter optimization process of adaptive SR based on PSO algorithm, d represents the number of multiparameters. at is, d � 5 for the AMPASR of Duffing system (adjusting parameters include k, a, b, ε and R), and d � 4 for the AMPASR of Langevin system (adjusting parameters include a, b, ε and R). As the PSO algorithm can only find out the minimal objective function among the solution space, the objective function of the optimization process is selected as the minus of output SNR, which is a commonly used SR performance index. In this paper, the input SNR and output SNR of the signals are defined as Here, SN(f) represents the single-side spectrum of system input sn(t), and X(f) the single-side spectrum of system output x(t), which can be solved by fourth-order Runge-Kutta algorithm in this paper. us, Am in � SN(f 0 ) and Am out � X(f 0 ) indicate the amplitude values corresponding to characteristic signal frequency at system input and output, respectively. Accordingly, the objective function in the multiparameter optimization of bistable SR systems is selected as f � − SNR out . us, the flowchart of the PSO algorithm used in this paper is presented in Figure 2.

Analysis of Adaptive Optimization Results.
In this subsection, the optimal SNR out of Duffing system equation (3) under given signal parameters will be obtained through solving the equation with different parameter sets and using the PSO algorithm, respectively. It is notable that, using the PSO algorithm, the Langevin system can produce similar simulated results, which are not presented in this paper to make it concise. us, the reliability and effectiveness of the PSO algorithm in AMPASR for bistable systems can be addressed.
First, the optimal SNR out is obtained by coadjusting the system parameters a and b. In Duffing system equation (3) Moreover, the largest system output SNR in Figure 3(a) appears at a � 0.24 and b � 0.1 with SNR out � − 4.8258 dB, and the largest output SNR in Figure 4(a) appears at a � 1.52 and b � 0.56 with SNR out � − 11.1684 dB. Again, the simulation results show the significance of parameter-adjusting method in improving the output SNR of bistable systems.
Next, the optimal SNR out is obtained through the PSO algorithm.
e same signal parameters of A � 0.1, f 0 � 0.01 Hz and the sampling frequency f s � 5 Hz are set in Duffing system equation (3). e proposed PSO algorithm, whose optimization parameters are set as N p � 20 and M � 20, is used to search the optimal SNR out and corresponding parameters set. e searching ranges for tunable parameters are a ∈ Update the position and velocity of each particle according to equations (5) and (6) Export the optimal SR results  It can be seen from Figure 5 that, in each optimization, the PSO algorithm can only find out a local optimal solution, resulting in a different optimal SNR out in each optimization. e optimal SNR out for D � 0.4 in 10 random optimizations range from − 8.035 dB to − 5.993 dB, which are within the 0.7% largest range of SNR out (SNR out ≥ − 8.83 dB) compared to Figure 3(b). Similarly, the optimal SNR out for D � 2.0 range from − 14.65 dB to − 11.68 dB, within the 0.4% largest range of SNR out (SNR out ≥ − 14.7 dB) compared to Figure 4(a). ese results indicate that although the PSO algorithm may not find the global optimal SNR out , the obtained local optimal SNR out is quite close to the global optimal one. us, the proposed PSO algorithm is feasible in the optimization of AMPASR. Moreover, the relationship between SNR out and the turntable parameters, as shown in Figures 3(a) and 4(a), is time-consuming to be obtained, while the PSO algorithm can find out an available optimal solution in a short time, thus providing a convenient method in practical applications such as online detections.
It should be pointed out that, in order to present the output results using clear three-dimensional images such as Figures 3(a) and 4(a), only two adjusting parameters (a and b) are considered in this subsection. Similar results are also obtained and show that PSO algorithm can find out an available optimal solution in a short time for all adjusting parameters.

Influence of PSO Algorithm Parameters on the Optimization Results.
In the PSO algorithm stated in Section 3.1, the particles' number N p and largest iteration number M are preset. In this subsection, through analyzing the optimal SNR out of a Duffing system, the influences of these two parameters on the output SR performance of bistable systems are discussed. e selection of the adjusting parameters' searching range should be stated first. e search range of k in Duffing system depends on the noise density of the input signal. e searching ranges of ε and R are determined by guaranteeing the amplitude and frequency of the characteristic signal component after transformation within the appropriate small-parameter range, while the searching ranges of systems a and b can be as large as possible. Without loss of generality, the upper limits of a and b are set as 5 in this section for the comparison of the optimization results.
Again, in Duffing system equation (3) First, the influence of N p on the system optimization results is studied. Table 1 presents the optimal SNR out and iteration time with M � 20 and different N p in 8 random optimizations. e optimization results are also plotted in Figure 6 with the data being in ascending order. It can be seen from Table 1 and Figure 6(a) that, in 8 random optimizations with different N p , the obtained optimal SNR out are within a narrow range around − 9.5 dB, which is significantly enhanced compared to SNR in , indicating that the value of N p barely affects the optimization results. However, Table 1 and Figure 6(b) show that, with the value of N p increasing, the iteration time significantly increases, thus reducing the optimization efficiency. erefore, a small value of N p , such as 10, is chosen in this paper in PSO optimizations with the acceptable optimization results and short iteration time being guaranteed.
Similarly, the influence of M on the system optimization results is studied as well. Table 2 and Figure 7 present the optimal SNR out and iteration time with N p � 10 and different M in 8 random optimizations. It can be seen from Figure 7(a) that the optimal SNR out with M � 20 and M � 30 are close and obviously better than those with M � 10, indicating that the largest iteration number M affects the optimization results when it is small, while when it is large enough, its influence can be ignored. erefore, in this paper M � 20 is chosen as an appropriate value of the largest iteration number that can produce acceptable optimization results in an appropriate iteration time.

Comparison between Optimization Results of Langevin
System and Duffing System. Duffing system (3) and Langevin system (4) are two common used bistable systems. Previous research shows that Duffing system has better adaptability to signals with large noise intensity due to the tunable damping ratio [25]. In this subsection, PSO algorithm is used for both systems to obtain their optimal output SNR under three groups of the same signals, which are given by A � 0.1, f 0 � 0.01 Hz, D � 0.4, 2, 5, the sampling frequency f s � 5 Hz, and the number of sampling points N � 20000 in equations (3) and (4). us, comparison of the optimization results is further presented.
In the optimization process of PSO algorithm, N p � 10 and M � 20 are selected based on the analysis in Section 4.1, and the search ranges for equations (3) and (4) Tables 3 and 4, respectively. It can be seen from Tables 3 and 4 that when the noise intensity is small (D � 0.4, SNR in � − 22.011 dB), the optimal SNR out and average SNR out in both systems are in the same level, indicating that the Duffing system and Langevin system provide similar optimization results under small noise intensity conditions. However, when the noise intensity increases (D � 2.0, SNR in � − 29.001 dB and D � 5.0, SNR in � − 32.980 dB), the optimal SNR out of Duffing system is larger than that of Langevin system, and more importantly, the average SNR out of Duffing system is much larger than that of Langevin system, indicating that the Duffing system can produce a desirable SNR out more easily by using PSO algorithm. ese can explain that the damping ratio k in   Shock and Vibration the Duffing system is helpful to deal with the large noise intensity signals. Moreover, the simulation processes show that divergence may easily occur in the optimization of Langevin system, and an optimal SNR out may not be obtained in some cases, which is another shortage of Langevin system compared to Duffing system.
Based on the previous analyses, the Duffing system, which will be selected in this paper for adaptive fault diagnosis, is suggested as a better bistable system in weaksignal detection.

Weak-Signal Detection Method Based on AMPASR of the
Duffing System. According to previous analyses, the AMPASR of Duffing system provides a solution for the adaptive detection of weak signals emerged in a strong noise background. One potential application of the proposed method is the diagnosis of rolling bearing incipient faults, whose vibration signals sn(t) are generally noisy signals with low SNR. Considering that a steady-state sn(t) in practical     engineering can be decomposed into some periodic signals and a source of noise, the AMPASR of Duffing system proposed in this paper can be used to produce optimal SR for sn(t) adaptively. us, the existence of the incipient fault can be judged. In practical vibration fault signals, the accurate characteristic frequency f 0 of the faults is always unknown. erefore, the definition of the output SNR should be modified by estimating the potential frequency domain [f l , f r ] of the fault signal. e modified SNR out , which is chosen as the objective function in the adaptive detection in practical applications, is given by SNR out � 10 log 10 max us, the procedure of bearing fault diagnosis by using the AMPASR of Duffing system, which corresponds to equation (3), is presented in Figure 8. In the procedure, the fault frequency domain [f l , f r ] should be estimated first, then an appropriate sampling frequency f s should be chosen to acquire the vibration signal sn(t). Before optimizations, the search range of multiparameters should be determined according to the vibration signal. In what is following in this section, three fault diagnosis examples of rolling bearings are presented to demonstrate the effectiveness and reliability of the proposed method.

Adaptive Diagnosis Examples of Rolling Bearing Faults.
e rolling bearing fault signals proposed in this subsection are cited from the Case Western Reserve University (CWRU) Bearing Data Center, which provides data presenting various degrees of difficulty for diagnosis [38]. e basic layout of the test rig is shown in Figure 9. It consists of a 2-horsepower reliance electric motor driving a shaft on which a torque transducer and encoder are mounted. Torque is applied to the shaft via a dynamometer and electronic control system.
In the tests, motor bearings were seeded with faults using electrodischarge machining (EDM). Faults ranging from 0.007 inches in diameter to 0.040 inches in diameter were introduced separately at the inner ring, outer ring, and rolling element. Faulted bearing was reinstalled into the test motor and vibration data was recorded for motor loads of 0 to 3 horsepower, and acceleration data was measured at locations near to and remote from the motor bearings.
In this subsection, the fault vibration signals from the drive end bearing with a motor load of 1 horsepower are used for verification. e rolling bearing is a deep groove ball bearing with the type of SKF 6205-2RS JEM. e geometric parameters and fault frequencies of the bearing are provided in 5. For the vibration signals with faults, characteristic frequencies f BPFI , f BPFO , and f BPFB are expected in the signal spectrums when the fault appears in inner ring, outer ring, and rolling element, respectively. During the test, an accelerometer was placed at the 12-o'clock position at the drive end of the motor housing; the approximate motor speed was 1774 rpm (29.57 Hz) and vibration signals were collected using a 16-channel DATrecorder at 48,000 samples per second (f s � 48000 Hz) for the drive end faults, and the data length is N � 10000.

Case 1: Inner Ring Fault Diagnosis.
e fault frequency of an inner ring fault can be calculated as f BPFI ≈ 160.1 Hz according to Table 5. e waveform and spectrum of the original vibration signal with inner ring fault, which is 0.021 inches in diameter indicating a severe fault condition, are presented in Figures 10(a) and 10(b). It can be seen that, due to the low input SNR (SNR in � − 34.99 dB), the signal component of inner ring fault frequency f BPFI ≈ 160.1 Hz is submerged in the waveform and is not evident compared to other signal components in the spectrum. us, the inner ring fault signal hidden in the original vibration signal cannot be identified.
In order to detect the potential inner ring fault of the rolling bearing, the proposed AMPASR method of Duffing system is adopted to analyze the vibration signal. In each optimization, the potential frequency domain is  Figures 10(c) and 10(d), where the optimal parameters are k � 9.3271, a � 0.1112, b � 0.5065, ε � 11.4685 and R � 6478, respectively. It can be seen from Figure 10(d) that the signal component of frequency f BPFI � 158.4 Hz is evident in the spectrum of output signal due to the greatly enhanced SNR out of − 12.78 dB. Moreover, the 3X and 4X higher-order harmonics can be observed as well. e output result indicates that an inner ring fault may exist in the drive end rolling bearing, which is accordant with the real physical truth of the rolling bearing. us, the adaptive diagnosis of an inner ring fault is realized by using the AMPASR method of Duffing system. Moreover, the same date has been analyzed in [27] using a manual adjusted tristable SR system. e output SNR in that reference (− 25.43 dB) is much smaller than that in Figure 10(d) (− 12.78 dB), thus proving the advantage of the proposed AMPASR method in adaptive rolling bearing inner ring fault diagnosis.
Moreover, the same vibration signal with inner ring fault can be also detected by the proposed AMPASR method of Langevin system. e optimal SNR out of the output signals from both Duffing system and Langevin system in 10 random adaptive detections are presented in Figure 11. It can be seen that although the Langevin system is helpful to enhance the SNR out of the detected signal, the optimal SNR out of the output signals from Langevin system is generally smaller than those from Duffing system. Again, the Duffing system is demonstrated as a better bistable system in    weak-signal detection compared with the Langevin system. us, in the following fault diagnosis cases, we only study the detection results from the Duffing system.

Case 2: Outer Ring Fault Diagnosis.
Similarly, the fault frequency of an outer ring fault can be calculated as f BPFO ≈ 106.0 Hz according to Table 5. e waveform and spectrum of the original vibration signal with outer ring fault, which is 0.021 inches in diameter, are presented in Figures 12(a) and 12(b). It can be seen that the signal component of outer ring fault frequency f BPFO ≈ 106.0 Hz cannot be identified in both the waveform and the spectrum due to the extremely low input SNR (SNR in � − 53.16 dB).
Again, the proposed AMPASR method of Duffing system is adopted to analyze the vibration signal. In each optimization, the potential frequency domain is    Once again, the proposed AMPASR method of Duffing system is adopted to analyze the vibration signal. In each optimization, the potential frequency domain is [135 Hz, 145 Hz], and the searching ranges of the tunable parameter are k ∈ [0.1, 15], a ∈ [0.1, 10], b ∈ [0.1, 10], ε ∈ [20, 50] and R ∈ [5000, 20000], respectively. e waveform and spectrum of the optimal output signal in 10 random optimizations are presented in Figures 13(c) and 13(d), where the optimal parameters are k � 8.0513, a � 6.1798, b � 9.6208, ε � 22.6837 and R � 12364, respectively. It can be seen that the signal component of frequency f BPFB � 139.2 Hz can be identified from the spectrum of output signal shown in Figure 13(d), where SNR out is enhanced to − 14.78 dB. is result indicates that a rolling element fault is in the drive end rolling bearing. us, the adaptive diagnosis of a rolling element fault is realized by using the AMPASR method of Duffing system.

Discussion of Practical Examples.
In this section, the feasibility of the proposed weak-signal detection method based on AMPASR of Duffing system has been demonstrated. It can be seen that the weak features of the inner ring, outer ring, and rolling elements faults can be greatly enhanced by adjusting the system parameters adaptively to produce an optimal system output SR. Moreover, the calculation time with the given optimization parameters is around 6 s, which is acceptable in practical engineering. erefore, the proposed approach provides a reliable solution for both online and offline adaptive diagnosis of rolling bearing faults.
Actually, the length of the detected vibration signal T has an influence on the detection results, especially for the online detections. It is well known that T � N/f s . e changes of T indicate two cases: (1) When N is fixed and f s is set according to the characteristic frequency f 0 of different vibration fault signals, we can simply change the searching range of the scale-transformation coefficient R to get a good optimization result; (2) when f s is fixed for a given vibration fault signal, the change of N indicates the change of T. On one hand, when N(T) is smaller, it costs less time for the adaptive optimization process and the detection is more efficient especially for online detection. On the other hand, when N(T) is larger, more information is involved in the signal and the spectrum resolution f s /N is larger, which is important for some complicated signals. erefore, we should set the value of N appropriately especially for online detections, or we can set different values of N to deal with the same vibration fault signal to find the best output result.
It should be pointed out that we always assumed that the fault types of rolling bearing are already known in the above examples. us, the corresponding fault frequencies, which are necessary for adaptive adjustments, can be calculated, and the feasibility of the proposed method can be verified. However, in most cases of practical engineering, the existence of fault and the potential fault type cannot be foreseen. erefore, one should calculate the fault frequencies of the inner ring, outer ring, and rolling elements based on the rotating frequency of the motor. Potential frequency domains are selected accordingly to obtain the optimal output result for each fault type. us, the existence of the three rolling fault types can be recognized and the fault detection can be realized.
In this section, we only consider the cases of single fault to demonstrate the feasibility of the proposed adaptive method. However, the weak vibration signals with a multifault signal are quite common in practical engineering. For these cases, one can select different frequency domains for different fault frequencies to realize the detections of multifaults. is is worthy of further study in the near future.

Conclusion
Aiming at the fault diagnosis of rolling bearings, an adaptive multiparameter-adjusting SR method for bistable systems based on PSO algorithm has been proposed and studied in this paper. e tunable parameters include the traditional nonlinear system parameters and other generalized parameters such as the amplitude-transformation coefficient and the scale-transformation coefficient. With the output SNR (SNR out ) as objective function, the optimal SR output of bistable systems can be obtained by adjusting the multiparameters adaptively for parameter-fixed vibration signals, and a larger SNR can be achieved in the output signal. us, potential weak rolling bearing fault features can be extracted using the proposed method and fault diagnosis can be realized. Simulation results show that although the adaptive method can only find out a local optimal solution in each optimization, the obtained SNR out is close to the optimal one and the weak signal feature can be extracted from the SR output. Based on the proposed AMPASR method, the influence of PSO algorithm parameters on the simulation results is studied, and a set of reasonable algorithm parameters including particles' number and largest iteration number is chosen for further simulations. Furthermore, the optimization results in the Langevin system and Duffing system are compared, and results show that the Duffing system, which is selected in this paper for adaptive fault diagnosis, is a superior bistable system in weak-signal detection. At last, a weak-signal detection method based on AMPASR of the Duffing system is proposed. is method is employed to adaptively diagnose the rolling bearing faults including inner ring fault, outer ring fault, and rolling element fault based on the data from CWRU, thus demonstrating the availability and reliability of the proposed method.
Data Availability e data analyzed in Section 5 can be accessed from the Case Western Reserve University Bearing Data Center Website (http://csegroups.case.edu/bearingdatacenter/pages/downloaddata-file).

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper. Shock and Vibration 13