Abstract
The 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. This 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.
1. 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. Therefore, 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]. This 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–12], structure monitoring [13], and energy harvesting [14–16]. Moreover, in view of the good performance of SR using noise to enhance periodic signal features, SR-based weak-signal detection and fault diagnosis methods have been investigated and successfully applied in rolling bearing fault diagnosis [17, 18].
The 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–23]. The latter ones, which involve the classical overdamped Langevin system and underdamped Duffing system, have much better SR performance than the former systems. Therefore, the bistable systems are mostly used as weak-signal detection models in previous studies.
SR describes an optimal synergistic action of the input signal, noise, and nonlinear systems. Therefore, 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. This method to achieve SR is the so-called parameter-adjusting SR method [24]. The 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. The adjustment process may be time-consuming and the optimal SR output may not be obtained by engineers without sufficient experiences. This 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–35]. These 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. Therefore, 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–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 multiparameter-adjusting SR (AMPASR) method based on bistable systems. The application of this method on rolling bearing fault diagnosis is also presented. The 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. 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 and a source of noise , and the other is described by equation (2) indicating an overdamped Langevin system [37] whose inertia term is ignored. In both bistable systems, denotes the force of the potential field, whose potential function with and being positive system parameters. In Duffing system equation (1) denotes the damping ratio. is defined as the system input signal; in the simulation analyses of this paper, represents a harmonic weak signal with amplitude and frequency , and denotes a source of Gaussian white noise with noise intensity :
In the absence of both and , the potential function in equations (1) and (2) describes a bistable potential field with two stable equilibrium points at and one unstable equilibrium point at , as shown in Figure 1. The height of the potential barrier that separates two potential wells is . This is why equations (1) and (2) are called bistable systems.
The output 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 under the coaction of damping force (only for the Duffing system equation (1)), potential field force , driving force and random noise . Through theoretical and simulation analyses, previous research has demonstrated that SR can occur in both Duffing system equation (1) and Langevin system equation (2). That is, with the help of nonlinear systems, the noise can produce a positive effect on the weak signal and enhance its intensity in the output . Therefore, 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 with mechanical faults, the signal parameters are fixed. Therefore, 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. The SR achieved through adjusting system parameters is the so-called parameter-adjusting SR.
Traditionally, the system parameters to be adjusted involve , , in Duffing system equation (1) and , 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 large-frequency component, an inadequate amplitude, and a low signal-to-noise ratio (SNR). Therefore, an amplitude-transformation coefficient and a scale-transformation coefficient 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 is used to transform the time/frequency scale of the input signal. The scale transformation is simply realized by setting the time step as instead of in numerical calculations, where indicates the sampling frequency. Thus, the system output in terms of time scale 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 , and system parameters (only for Duffing system), , . The SR-oriented adjustment rules for each parameter have been studied and summarized. These 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. Therefore, 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.
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. Thus, an AMPASR method based on PSO algorithm will be proposed to achieve SR adaptively in bistable systems. The optimization results will be further analyzed in this section.
3.1. Introduction of PSO Algorithm
The 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. The velocity is dynamically adjusted by the moving experiences of the particle itself and other particles. The optimal solution will be updated according to the particles’ best fitness.
The PSO algorithm can be described in mathematics. Assume that there is a population composed of particles in a -dimensional solution space. The position and velocity of the ith particle in the population are and , respectively. The optimal solutions for each particle and the whole population are and , respectively. Before the particles find out these two optimal solutions, the particle updates its velocity and position according towhere is the iteration number with the largest one; is the inertia weight; and are the positive learning factors; and are evenly distributed random numbers among . In this paper, and are chosen as 2, and is decided by an adaptive weighting method:where and are the given maximal and minimal values of ; means the particle’s current fitness value of the objective function, and 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, represents the number of multiparameters. That is, for the AMPASR of Duffing system (adjusting parameters include , , , and ), and for the AMPASR of Langevin system (adjusting parameters include , , and ). 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 asHere, represents the single-side spectrum of system input , and the single-side spectrum of system output , which can be solved by fourth-order Runge–Kutta algorithm in this paper. Thus, and 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 . Thus, the flowchart of the PSO algorithm used in this paper is presented in Figure 2.
3.2. Analysis of Adaptive Optimization Results
In this subsection, the optimal 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. Thus, the reliability and effectiveness of the PSO algorithm in AMPASR for bistable systems can be addressed.
First, the optimal is obtained by coadjusting the system parameters and . In Duffing system equation (3), the signal parameters are set as , , the sampling frequency , and the number of sampling points . By setting other tunable parameters as constant (, and ), Figures 3(a) and 4(a) present of Duffing system against system parameters and under ( dB) and ( dB), respectively. It can be seen from Figures 3(a) and 4(a) that the system output SNR generally presents a decreasing trend with and being increased, but the distribution rules of values are much more complicated. The amounts of within different ranges for Figures 3(a) and 4(a) are counted and presented in Figures 3(b) and 4(b), respectively. It can be seen that among values, 426 values appear in the range of dB in Figure 3(a) (probability of 0.7%) and 242 values appear in the range of dB in Figure 4(a) (probability of 0.4%). Moreover, the largest system output SNR in Figure 3(a) appears at and with dB, and the largest output SNR in Figure 4(a) appears at and with dB. Again, the simulation results show the significance of parameter-adjusting method in improving the output SNR of bistable systems.
(a)
(b)
(a)
(b)
Next, the optimal is obtained through the PSO algorithm. The same signal parameters of , Hz and the sampling frequency Hz are set in Duffing system equation (3). The proposed PSO algorithm, whose optimization parameters are set as and , is used to search the optimal and corresponding parameters set. The searching ranges for tunable parameters are , , while , and are set as those in Figures 3 and 4. The results of 10 optimizations for and are plotted in Figure 5 with presented in ascending order and the x-axis indicating the order number of each optimization result.
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 in each optimization. The optimal for in 10 random optimizations range from −8.035 dB to −5.993 dB, which are within the 0.7% largest range of ( dB) compared to Figure 3(b). Similarly, the optimal for range from −14.65 dB to −11.68 dB, within the 0.4% largest range of ( dB) compared to Figure 4(a). These results indicate that although the PSO algorithm may not find the global optimal , the obtained local optimal is quite close to the global optimal one. Thus, the proposed PSO algorithm is feasible in the optimization of AMPASR. Moreover, the relationship between 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 ( and ) 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.
4. Discussion
4.1. Influence of PSO Algorithm Parameters on the Optimization Results
In the PSO algorithm stated in Section 3.1, the particles’ number and largest iteration number are preset. In this subsection, through analyzing the optimal of a Duffing system, the influences of these two parameters on the output SR performance of bistable systems are discussed.
The selection of the adjusting parameters’ searching range should be stated first. The search range of in Duffing system depends on the noise density of the input signal. The searching ranges of and 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 and can be as large as possible. Without loss of generality, the upper limits of and are set as 5 in this section for the comparison of the optimization results.
Again, in Duffing system equation (3), the signal parameters are set as , Hz, ( dB), and the sampling frequency Hz, sampling points . The PSO algorithm is adopted to search the in the searching ranges of , , , and . It is necessary to first state that the computer used for numerical simulations has CPU of Intel(R) Core(TM) i7-6700 and RAM with 8.00 GB as its main configurations.
First, the influence of on the system optimization results is studied. Table 1 presents the optimal and iteration time with and different in 8 random optimizations. The 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 , the obtained optimal are within a narrow range around −9.5 dB, which is significantly enhanced compared to , indicating that the value of barely affects the optimization results. However, Table 1 and Figure 6(b) show that, with the value of increasing, the iteration time significantly increases, thus reducing the optimization efficiency. Therefore, a small value of , such as 10, is chosen in this paper in PSO optimizations with the acceptable optimization results and short iteration time being guaranteed.
(a)
(b)
Similarly, the influence of on the system optimization results is studied as well. Table 2 and Figure 7 present the optimal and iteration time with and different in 8 random optimizations. It can be seen from Figure 7(a) that the optimal with and are close and obviously better than those with , indicating that the largest iteration number affects the optimization results when it is small, while when it is large enough, its influence can be ignored. Therefore, in this paper is chosen as an appropriate value of the largest iteration number that can produce acceptable optimization results in an appropriate iteration time.
(a)
(b)
4.2. 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 , Hz, , the sampling frequency Hz, and the number of sampling points in equations (3) and (4). Thus, comparison of the optimization results is further presented.
In the optimization process of PSO algorithm, and are selected based on the analysis in Section 4.1, and the search ranges for equations (3) and (4) are (only for Duffing system (3)), , , and , respectively. The optimization results of Duffing system and Langevin system are recorded in Tables 3 and 4, respectively. It can be seen from Tables 3 and 4 that when the noise intensity is small (, dB), the optimal and average 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 (, dB and , dB), the optimal of Duffing system is larger than that of Langevin system, and more importantly, the average of Duffing system is much larger than that of Langevin system, indicating that the Duffing system can produce a desirable more easily by using PSO algorithm. These can explain that the damping ratio in 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 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 weak-signal detection.
5. Practical Examples
5.1. 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 are generally noisy signals with low SNR. Considering that a steady-state 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 adaptively. Thus, the existence of the incipient fault can be judged.
In practical vibration fault signals, the accurate characteristic frequency of the faults is always unknown. Therefore, the definition of the output SNR should be modified by estimating the potential frequency domain of the fault signal. The modified , which is chosen as the objective function in the adaptive detection in practical applications, is given by
Thus, 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 should be estimated first, then an appropriate sampling frequency should be chosen to acquire the vibration signal . 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.
5.2. Adaptive Diagnosis Examples of Rolling Bearing Faults
The 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]. The 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. The rolling bearing is a deep groove ball bearing with the type of SKF 6205-2RS JEM. The geometric parameters and fault frequencies of the bearing are provided in 5. For the vibration signals with faults, characteristic frequencies , , and 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 DAT recorder at 48,000 samples per second ( Hz) for the drive end faults, and the data length is
5.2.1. Case 1: Inner Ring Fault Diagnosis
The fault frequency of an inner ring fault can be calculated as Hz according to Table 5. The 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 ( dB), the signal component of inner ring fault frequency Hz is submerged in the waveform and is not evident compared to other signal components in the spectrum. Thus, the inner ring fault signal hidden in the original vibration signal cannot be identified.
(a)
(b)
(c)
(d)
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 , and the searching ranges of the tunable parameter are , , , and , respectively. It should be pointed out that the search ranges of and are determined by guaranteeing the amplitude and frequency of the fault signal component after transformation among the appropriate small-parameter range. The waveform and spectrum of the optimal output signal in 10 random optimizations are presented in Figures 10(c) and 10(d), where the optimal parameters are , , , and , respectively. It can be seen from Figure 10(d) that the signal component of frequency Hz is evident in the spectrum of output signal due to the greatly enhanced of −12.78 dB. Moreover, the 3X and 4X higher-order harmonics can be observed as well. The 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. Thus, 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. The 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. The optimal 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 of the detected signal, the optimal 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. Thus, in the following fault diagnosis cases, we only study the detection results from the Duffing system.
5.2.2. Case 2: Outer Ring Fault Diagnosis
Similarly, the fault frequency of an outer ring fault can be calculated as Hz according to Table 5. The 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 Hz cannot be identified in both the waveform and the spectrum due to the extremely low input SNR ( dB).
(a)
(b)
(c)
(d)
Again, the proposed AMPASR method of Duffing system is adopted to analyze the vibration signal. In each optimization, the potential frequency domain is , and the searching ranges of the tunable parameter are , , , and , respectively. The waveform and spectrum of the optimal output signal in 10 random optimizations are presented in Figures 12(c) and 12(d), where the optimal parameters are , , , and , respectively. One can see that the signal component of frequency Hz can be identified from the spectrum of output signal shown in Figure 12(d), where is enhanced to −14.64 dB. This result indicates the existence of an outer ring fault in the drive end rolling bearing. Thus, the adaptive diagnosis of an outer ring fault is realized by using the AMPASR method of Duffing system.
5.2.3. Case 3: Rolling Element Fault Diagnosis
At last, the vibration signal with a rolling element fault is considered. The fault frequency of a rolling element fault can be calculated as Hz according to Table 5. The waveform and spectrum of original vibration signal with rolling element fault, which is 0.021 inches in diameter, are presented in Figures 13(a) and 13(b), where the signal component of rolling element fault frequency Hz is completely submerged due to the extremely low input SNR ( dB).
(a)
(b)
(c)
(d)
Once again, the proposed AMPASR method of Duffing system is adopted to analyze the vibration signal. In each optimization, the potential frequency domain is , and the searching ranges of the tunable parameter are , , , and , respectively. The 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 , , , and , respectively. It can be seen that the signal component of frequency Hz can be identified from the spectrum of output signal shown in Figure 13(d), where is enhanced to −14.78 dB. This result indicates that a rolling element fault is in the drive end rolling bearing. Thus, the adaptive diagnosis of a rolling element fault is realized by using the AMPASR method of Duffing system.
5.3. 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. Therefore, 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 has an influence on the detection results, especially for the online detections. It is well known that . The changes of indicate two cases: (1) When is fixed and is set according to the characteristic frequency of different vibration fault signals, we can simply change the searching range of the scale-transformation coefficient to get a good optimization result; (2) when is fixed for a given vibration fault signal, the change of indicates the change of . On one hand, when () 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 () is larger, more information is involved in the signal and the spectrum resolution is larger, which is important for some complicated signals. Therefore, we should set the value of appropriately especially for online detections, or we can set different values of 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. Thus, 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. Therefore, 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. Thus, 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. This is worthy of further study in the near future.
6. 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. The 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 () 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. Thus, 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 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. This 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
The data analyzed in Section 5 can be accessed from the Case Western Reserve University Bearing Data Center Website (http://csegroups.case.edu/bearingdatacenter/pages/download-data-file).
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by National Natural Science Foundation of China (Grant no. 51905349), Natural Science Foundation of Jiangxi Province (CN) (Grant no. 20161BAB216111), Natural Science Foundation of Shenzhen University (Grant no. 2019036, and 860-000002110264), and Key Laboratory of Lightweight and High Strength Structural Materials of Jiangxi Province (Grant no. 20171BCD40003).