Abstract

Because of the characteristics of weak signal and strong noise, the low-speed vibration signal fault feature extraction has been a hot spot and difficult problem in the field of equipment fault diagnosis. Moreover, the traditional minimum entropy deconvolution (MED) method has been proved to be used to detect such fault signals. The MED uses objective function method to design the filter coefficient, and the appropriate threshold value should be set in the calculation process to achieve the optimal iteration effect. It should be pointed out that the improper setting of the threshold will cause the target function to be recalculated, and the resulting error will eventually affect the distortion of the target function in the background of strong noise. This paper presents an improved MED based method of fault feature extraction from rolling bearing vibration signals that originate in high noise environments. The method uses the shuffled frog leaping algorithm (SFLA), finds the set of optimal filter coefficients, and eventually avoids the artificial error influence of selecting threshold parameter. Therefore, the fault bearing under the two rotating speeds of 60 rpm and 70 rpm is selected for verification with typical low-speed fault bearing as the research object; the results show that SFLA-MED extracts more obvious bearings and has a higher signal-to-noise ratio than the prior MED method.

1. Introduction

Rolling bearings are one of the most widely used elements in rotary machines. The mechanical failure of bearings may result in financial loss, and even death. Therefore, condition monitoring and fault diagnosis of rolling bearings are critically important for production efficiency and plant safety in modern enterprises [1, 2]. When a rolling bearing fault is weak at an early stage or at a low shaft speed, weak fault features are often embedded in background noise. So it is not an easy task to extract the representative features from the original signal [3]. For instance, the fast Fourier transform cannot obtain ideal transient extraction, even as an increased signal-to-noise ratio assists weak fault feature extraction. Many scholars have reviewed the problem of fault diagnosis taken from weak signals and have made progress. In [4], wavelet transformation is applied to the problem of weak, abnormal vibration signal extraction, but reliance upon a single wavelet transform makes the extraction of complex vibration signals difficult [5, 6]. Many approaches to improved wavelet transform signal extraction have been proposed, by combination with other diagnostic tools: neutral networks, hidden Markov models, or singular value decomposition. In some cases, good results were obtained, but these approaches require the setting of a great many parameters, and yet when the signal is quite weak or the noise strong, these approaches give poor results [710]. Huang et al. in 1998 proposed empirical mode decomposition (EMD), an adaptive signal processing method that yields intrinsic mode functions (IMFs) [11]. EMD algorithms are widely used for rolling bearing fault type identification. However, there is a serious pattern mixing and endpoint effect in the EMD method. Wu and Huang proposed ensemble empirical mode decomposition (EEMD) to overcome the shortcoming of mode mixing, but end effects remain a problem. So the EMD and EEMD methods cannot be qualified for weak fault diagnosis [12, 13]. Stochastic resonance is usually used to enhance weak fault signals and to extract rolling bearing characteristic frequency components [14, 15]. Given the difficulty of choosing the parameters of a band-pass filter, spectral kurtosis (SK) is another useful tool; it provides a moderate way of detecting weak fault in rotation even under the background of strong noise [16, 17]. Liu et al. put forward a new method to extract new features using kernel joint approximate diagonalization of eigenmatrices (KJADE) [18]. Sparse decomposition derived from application of a matching pursuit algorithm to fault signal data can detect the amplitude, the frequency, and the phase of weak signals under the background of strong noise [19]. These results are not obvious because of the loss of signal details.

Wiggins introduced minimum entropy deconvolution (MED) in 1978. MED has been applied to seismic signal processing, especially for the separation of convolutional components of reflection signals [20]. In 2007, Sawalhi et al. first applied MED to detection of faults in rolling bearings and other gear. MED is usually combined with spectral kurtosis or AR filtering to enhance rolling bearing fault detection and diagnosis [21, 22]. MED is also a valid approach to the period estimation problem [23]. However, given strong noise, the failure impact component will be lost by MED, and this is because of the following reasons [2426].

MED is used as a linear operator and so is ill-adapted to process a limited frequency bandwidth. For noisy data, the limitation of the linear operator is difficult to overcome.

The filter coefficients obtained by the objective function method (OFM) are local optima, not global optima.

To address the latter limitation, an improved MED method based on the shuffled frog leaping algorithm (SFLA-MED) is proposed. In the new method, a shuffled frog leaping algorithm (SFLA) is applied to select an optimal set of filter coefficients that provide maximum kurtosis. So, the filtered signal will be acquired by taking advantage of a global optimal inverse filter. The envelope spectrum of the filtered signal is calculated with Hilbert transform (HHT) and fast Fourier transform (FFT). In the envelope spectrum, the characteristic frequency of a rolling bearing is quite obvious, and the bearing faults can be diagnosed.

2. Overview of the Improved Algorithm

2.1. Minimum Entropy Deconvolution Algorithm

Deconvolution is an inverse filter that is typically applied to seismic data, for the purpose of recovering reflection coefficients. The main benefit is compression of the seismic reflection pulse and to improve estimation of the subsurface reflection interface reflection coefficient. Deconvolution is particularly important to the bright spot technique of oil-gas exploration and to the analysis of seismic data associated with lithological study on the formation of seismology. Deconvolution can remove the interference from short cycle reverberation and multiple waves [27, 28]. Deconvolution can be divided into deterministic deconvolution and predictive deconvolution. There are many commonly known deconvolutions, including east square deconvolution, predictive deconvolution, homomorphic deconvolution, surface-consistent deconvolution, maximum entropy deconvolution, minimum entropy deconvolution, change mould deconvolution, deconvolution, Noah deconvolution, and minimum information deconvolution.

The MED technique is a type of system identification method that was originally developed to aid extraction of reflectivity information in seismic data. The reflectivity information can be used to identify and locate layers of subterranean minerals [23]. The essence of the method is to acquire a linear operator that maximizes the characteristics of the spike pulse [29, 30]. Maximum kurtosis defines a stop condition [30].

Rolling bearing fault vibration signal can be expressed as where is the impact sequence of the rolling bearing vibration signal, is the transfer function, and is the noise component. To facilitate analysis, can be ignored. The purpose of MED is to find a deconvolution optimal inverse transfer function to calculate an input signal from the output signal through .

From (3) one may find the best inverse filter that maximizes the kurtosis index of original signal.

Therefore, MED makes the recovery of signal by construction of an “optimal” filter possible.

MED can be summarized in five steps [26]:(1)Initialization:(2)Iteration computing:(3)Calculate the cross-correlation:(4)Calculate the coefficient of filtering:(5)Terminate condition.

Assume that the given threshold is 0.01. If , the iteration ends; else , and iterate.

2.2. Shuffled Frog Leaping Algorithm

SFLA combines the global breadth search of the entire frog populations and the local depth search of frogs individual information, leading the algorithm towards the global optimum direction [25, 31]. The frog representing the solution is divided into several subgroups, and each population has its own culture. Then the frogs evolve according to fitness in the population. After searching in each memeplex, the frogs from all memeplexes are shuffled and then frogs are redistributed forming new memeplexes, which makes searching process less possible to be trapped in local optimum [31].

Its steps include the following:(1)Initialize population size.(2)Initialize population generation ; calculate the performance of : (3)Rank the frogs , and arrange the array according to :(4)The frogs are placed in different populations in turn.(5)Within populations, frogs begin to evolve (update the worst frogs in the population):(6)After the evolution of the population, the subgroups are combined and arranged in descending order. By constantly iterating, then find out the best frogs to update.

The output signal of a linear time invariant system (LTI) system can be expressed as follows: is the input signal, and represent filtering coefficients, and is the filter order. In this paper, ,The best filter coefficient obtained by the leapfrog algorithm is substituted into the formula ; then the denoising signal becomes

Then the denoising signal is demodulated by HILBERT-FFT envelope. Finally, the validity of SFLA-MED is verified by comparing the envelope spectrum with the fault characteristic value.

2.3. Diagnosis Method by the SFLA-MED

Figure 1 shows the differences between MED and the SFLA-MED. The MED filter is implemented by the objective function method (OFM) given in reference [26]. The OFM is an optimization process designed to maximize the kurtosis of the MED output. It is pointed out that the filter coefficients obtained by OFM are merely local optima, not global optima. So the SFLA-MED uses SFLA to find these optimal coefficients.

The process of fault diagnosis is shown in Figure 2. Vibration signals are measured in a low-speed bearing experiment. The SFLA-MED method does not need to set threshold parameters during the iteration and directly filters the original signal so that a new signal formed. Envelope demodulation is then applied to the new signal. Envelope demodulation consists of a Hilbert transform (HHT) and a fast Fourier transform (FFT). To further verify the reliability and feasibility, the proposed method is compared to MED and theoretical characteristic frequency.

3. Experimental Validation

3.1. Simulation Signal Analysis

Assume the expression of the bearing simulation signal:

is the amplitude of the shock component; the value was set to 1.5. The characteristic frequency  Hz,  Hz, and carrier frequency  Hz. The noise component is , amplitude was 0.2. The time-domain waveform is shown in Figure 3.

Hilbert-FFT was used to demodulate the signals obtained by the MED method and the SFLA-MED method, and then the amplitude of the frequency domain was normalized for a clear contrast. As can be seen from Figure 4(a), after the signal had been reduced and processed by the MED method, the characteristic frequency was drowned in the noise signal, which cannot accurately extract the fault characteristics of the bearing. After the SFLA-MED method was adopted, the signal can be accurately extracted with 12 Hz and 30 Hz, and the harmonic characteristics were obvious, as shown in Figure 4(b). In addition, as can be seen from Figure 3, the modulation characteristics of the simulation signal were obvious, so the envelope analysis was carried out here, as shown in Figure 4(c). After the signal was enveloped, the submerged feature of the mixed white noise was well detected. It should be pointed out that the noise contained in the actual bearing fault signal was usually rather complicated. So the paper uses the method of SFLA-MED and envelope analysis to detect the weak fault signal of bearing.

3.2. Experimental Verification of Low-Speed Bearings
3.2.1. Experimental System

The experimental signals to be analyzed were collected from low-speed bearing experiments as shown in Figure 5. In order to better verify the effectiveness of the method put forward in this paper, the outer-race fault data and inner-race fault data of a rolling bearing were taken, and those flaws were artificially made with the use of a wire-cutting machine as shown in Figure 6. The sizes of the flaws are as follows:Outer-race flaw: 5.0 mm 0.5 mm (width depth)Inner-race flaw: 5.0 mm 0.5 mm (width depth).

3.2.2. Bearing Data Set Description

In order to fully analyze signal features and acquire more comprehensive information for the research of the fault diagnosis, the sampling frequency is set at 100 kHz, the sampling time is 10 seconds, and the rotation speed of spindle is 60 rpm and 70 rpm. So four kinds of vibration signal are obtained and they are outer-race fault vibration signal at 60 rpm, inner-race fault vibration signal at 60 rpm, outer-race fault vibration signal at 70 rpm, and the inner-race fault vibration signal at 70 rpm, respectively. By varying rotation speed, the time-domain wave form of vibration signals is shown in two figures. Figure 7 indicates the time-domain wave form of outer-race and inner-race fault vibration signals when the rotation speed is 60 rpm and Figure 8 shows the time-domain waveform of outer-race and inner-race vibration signals when the rotation speed is 70 rpm.

Bearing defects will generate a series of impact vibrations that emit at bearing characteristic frequencies every time a running roller passes over a defect [32]. Therefore the fault type can be determined by identifying the characteristic frequency [33].

Computed pass-frequencies are shown in Table 1. Pure rolling motion is a condition of the above equations. As there may be some sliding motion in practice, the results in Table 1 are approximate values.

4. Results and Discussion

4.1. Feature Extraction of MED and SFLA-MED

Envelop spectrum analysis based on the Hilbert transform has been widely used in bearing fault diagnosis [34]. It separates modulating signal from the collected signal, and this contains information that indicates the fault type. Envelope spectrum analysis is obtained from the FFT of the envelop signal. Figure 9 shows the normalized envelop spectrum of outer-race and inner-race fault vibration signals at the speed of 60 rpm filtered by MED method. As shown in Figure 9(a), the theoretical characteristic frequency is not clearly shown by the spectrum peaks. Although there is a clear spectrum peak appearing near the theoretical characteristic frequency as shown in Figure 9(b), it is not the highest spectral peak and cannot be identified for the first time. Hence, the MED method extracts the fault characteristic frequency with contingency.

Shuffled Frog Leaping Algorithm optimization improves the MED method. Based on theoretical study, the frequency equation of bearing fault characteristics is based on the assumption of a pure rolling motion. These vibrations occur at bearing characteristic frequencies, which are estimated based on the geometry of the bearing, its rotational speed, and the location of the defect. However, in practice, some unexpected status may occur, which causes slight deviation of the characteristic frequency locations. Figure 10 shows the normalized envelop spectrum of outer-race and inner-race fault vibration signals at the speed of 60 rpm filtered by SFLA-MED. From Figure 10, it can be seen that there is an obvious spectrum peak corresponding to the theoretical characteristic frequency, so the results obtained from SFLA-MED meet the diagnostic accuracy requirement. In summary, the SFLA-MED method can better extract the fault characteristic frequency.

According to the experimental results in Figure 11, the traditional MED method can recognize the bearing failure, which shows the effectiveness of the MED method. The SFLA-MED are used to analyze the experimental data, and it can also identify bearing failure, as shown in Figure 12. The two methods are good to identify fault bearing fault characteristics under 70 rpm. Compared with the 60 rpm experimental data, the SFLA-MED algorithm recognizes the characteristic frequency of the bearing fault under different speeds, which shows that it has a better universality and is more suitable for application.

4.2. Effect Comparison of MED and SFLA-MED

As is known, signal-to-noise ratio (SNR) is one of the most basic indicators to measure an algorithm. Combining the experimental data of the four groups, it is found that the characteristic frequency signals of the SFLA-MED method are not drowned by noise, and they are more easily identified with the fault characteristics (highest spectrum peaks). In addition, compared with the experimental result of signals at the speed of 60 rpm, the SFLA-MED method has higher signal-to-noise ratio than the original MED method as shown in Table 2. Therefore, the SFLA-MED can extract the failure frequency more accurately and reduce the possibility of miscalculation.

The calculation formula of frequency domain signal-to-noise ratio is shown in. is the sampling length, is the corresponding amplitude of the first three-order characteristic frequency, and is the frequency domain signal amplitude.

5. Conclusions

In this paper, an SFLA-MED combined with envelope demodulation is applied to the weak fault diagnosis of rolling bearings. Compared with the previous MED method, the proposed method extracts more obvious fault features of rolling bearings and is better adapted to engineering application. It is pointed out that the improved MED obtains a global optimal solution but that MED does not. SFLA makes the selection of filtering coefficients more flexible. As can be seen from the algorithm flow diagram, there is no error accumulation in the SFLA-MED, because the process of finding filter coefficients is not affected by the last filtered signal. There are no extra frequency components beyond the characteristic frequency that appears in the envelope spectrum.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

The presented work was carried out by collaboration of all authors. Huaqing Wang conceived and designed the experiments; Junlin Li and Jingshen Jiang performed the experiments and the simulations; Xiaohong Fan and Liuyang Song analyzed the data; Wenbin Liu, Jianfeng Yang, and Liangchao Chen performed the simulation and revised the paper; Junlin Li and Huaqing Wang wrote the paper. All authors read and approved the final manuscript.

Acknowledgments

This paper is supported by the National Natural Science Foundation of China (Grant no. 51675035) and PetroChina Science and Technology Innovation Fund (Grant no. 2015D50060606).