Abstract

Aiming at reducing the production downtime and maintenance cost, prognostics and health management (PHM) of rotating machinery often includes the remaining useful life (RUL) prediction of bearings. In this paper, a method combining the generalized Weibull failure rate function (WFRF) and radial basis function (RBF) neural network is developed to deal with the RUL prediction of bearings. A novel indicator, namely, the power value on the sensitive frequency band (SFB), is proposed to track bearing degradation process. Generalized WFRF is used to fit the degradation indicator series to reduce the effect of noise and avoid areas of fluctuation in the time domain. RBF neural network is employed to predict the RUL of bearings with times and fitted power values at present and previous inspections as input. Meanwhile, the life percentage is selected as output. The performance of the proposed method is validated by an accelerated bearing run-to-failure experiment, and the results demonstrate the advantage of this method in achieving more accurate RUL prediction.

1. Introduction

Prognostics and health management (PHM) has received extensive attention in recent years for its effectiveness in improving reliability and reducing costs [1, 2]. PHM generally combines condition monitoring, fault diagnostics, fault prognostics, and decision support [3]. Remaining useful life (RUL) prediction aims at assessing the performance degradation of equipment and detecting the impending failure [4]. Therefore, accurate RUL prediction is regarded as one of the most central components in performing PHM, since necessary maintenance actions are implemented based on the RUL prediction result [5].

Bearing is the basic component of rotating machinery and the health condition of it plays an important role in the machinery. In order to predict the RUL of bearings, many kinds of measurements are considered, including acoustic, temperature, oil analysis, and pressure. Among them, vibration is most commonly used because it is easy to measure and analyze. Besides, it can reflect the status of bearings roundly and timely [6]. It mainly takes two steps to realize accurate RUL prediction of bearings in vibration-based methods. Firstly, an appropriate degradation indicator should be selected from various features. This step is the premise of RUL prediction and the selection result is crucial to the prediction accuracy. Secondly, an optimal prediction model should be developed. For the former step, there are three conventional feature extraction methods including time domain, frequency domain, and time-frequency domain. Time domain methods are directly based on the time waveform to obtain the factors such as root mean square (RMS), kurtosis, and peak. Most frequency domain methods utilize fast Fourier transformation (FFT) to convert time domain signal into frequency domain. Thus it is easy to identify the frequency components and recognize the fault characteristic frequency. Time-frequency domain methods investigate signal in both time and frequency domain, and the change of frequencies with time is described clearly. Typical time-frequency domain methods are short time Fourier transformation (STFT), Wavelet analysis, Wigner-Ville distribution, and Hilbert Huang transformation. After feature extraction, the degradation indicator is selected from them by certain standard. Various indicators have been utilized to depict the degradation process in previous research. Mahamad et al. [7] selected RMS and kurtosis obtained from time waveform, while Huang et al. [8] utilized the minimum quantisation error (MQE) index produced by three time features and three frequency features. Ocak et al. [9] and Pan et al. [10] decomposed the signals into wavelet packets and used the node energies of decomposition tree to indicate bearing degradation. Zhang et al. [11] firstly decomposed the preprocessed signals and then selected the peak values of FFT for decomposed signals to judge the degradation of the monitored machine. In order to extract the most relevant features to the degradation process, Liao and Lee [12] and Lu et al. [13] used principal component analysis (PCA) to eliminate the superfluous and less-relevant features.

After the selection of degradation indicator, a prediction model should be established to perform the RUL prediction. There are mainly two categories of prediction methods: model-based methods and data-driven methods. An accurate physical model and a correct fault propagation model are essential to a model-based method. Liu and Mahadevan [14] proposed a unified multiaxial fatigue damage model for rolling contact fatigue, while Xu et al. [15] developed an improved Paris model to predict residual fatigue life of bearings online. Jin et al. [16] built a degradation model for lubricant loss in bearing to predict the RUL of an individual momentum wheel bearing. If precise models are obtained, the model-based methods can provide satisfying prediction results. However, due to the complex structure, it is hard to build the physical-based models. In contrast, the data-driven methods are practical and easy to be operated, since they predict the RUL only based on the condition monitoring data. Moreover, the prediction result can be updated with new inspection data available. Yan and Lee [17] utilized logistic regression to achieve performance degradation assessment. Pham and Yang [18] developed a liner autoregressive moving average (ARMA) model and a nonlinear generalized autoregressive conditional heteroscedasticity (GARCH) model to explain the fault condition of machine. Fei and He [19] applied multiple-kernel relevance vector machine (MkRVM) as an intelligent system to predict the state of bearings. Dong et al. [20] employed the support vector machine (SVM) to track the degradation process of bearings and utilized Markov model to improve the prediction accuracy. Artificial neural network (ANN) has also been widely used to RUL prediction due to the well performance of function approximation. Tian [21] proposed an ANN to predict the RUL of pump bearings, with ages and velocity measurements as inputs and the life percentage as output. Wu et al. [22] and Huang et al. [8] developed a back-propagation (BP) neural network for estimating the failure time of bearings. Lee et al. [23] utilized an Elman neural network for health condition prediction. In order to achieve more accurate prediction, Shao and Nezu [24] applied different neural networks to different running stages of bearings. Gebraeel et al. [25] proposed two classes of neural network models to perform the RUL prediction, and the result testified the advantages.

This paper proposes a hybrid method for RUL prediction utilizing the generalized Weibull failure rate function (WFRF) and radial basis function (RBF) neural network. The proposed method differs from other methods which have been used in past literature in adopting a novel indicator to capture the performance degradation of bearings. A new concept named sensitive frequency band (SFB) is defined and the power value on the SFB is selected as the degradation indicator. The SFB power value series along the lifetime is calculated and fitted by the generalized WFRF, then the fitted measurements are fed into RBF neural network. The remainder of this paper is organized as follows. Section 2 presents the different steps of the proposed method. In this section, the reason of degradation indicator selection is given and the basic principles of generalized WFRF and RBF neural network are detailed. Experimental setup is shown in Section 3 and the datasets used for validation are described. Section 4 is dedicated to evaluating the proposed method with the experimental data. Finally, some conclusions are drawn in Section 5.

2. Description of the Proposed Approach

As mentioned above, it takes several steps to implement the proposed method. In this section, the procedure of the proposed method is discussed and the relevant theoretical bases are introduced in detail.

2.1. Indicator Selection

An appropriate indicator that can capture the evolution of the bearing’s degradation is vital to perform RUL prediction [8]. A good indicator should be sensitive to initial degradation and be consistent with degradation process [30]. As a bearing fatigue occurs, a spall develops on the surface of the race and the defective frequency is excited. During the process of the bearing’s degradation, the amplitude of this frequency is proportional to the spall size [31]. In other words, the amplitude of the defective frequency is the most intuitive and accurate reaction of the bearing’s degradation process. In view of this, the average amplitude of the defective frequency and its first six harmonics is selected as indicator in [25] to predict the RUL of bearings. However, there exist several problems. Firstly, the diagnostics will be performed many times from beginning to initial failure occurrence, in order to identify the fault pattern and determine the defective frequency. It is meaningless and time-consuming when the bearing is normal for a long time. Secondly, it is not easy to calculate the defective frequency at the stage of incipient failure because of the heavy noise and the weak signature. Finally, there is always a frequency modulation phenomenon due to the bearing slackness [8], resulting in dramatic changes of amplitudes on modulated frequency band rather than defective frequency. To overcome these problems, the power value on the SFB is developed as a novel degradation indicator.

SFB is defined as a frequency band on which the FFT spectrums along the bearing degradation process are significantly different. It is an isolation of the frequencies associated with bearing defects from overall frequencies. As an example, Figure 1 presents a plot of FFT spectrums for a bearing at different stages of its life cycle. It is obvious that the amplitudes of frequency lines on [3000 Hz, 5000 Hz] increase more dramatically than other frequency bands with time until the bearing failure.

From the physical interpretation of the SFB power value we can see this indicator is robust to failure categories. Comparing with several statistical features usually adopted in previous studies, such as RMS and kurtosis, the SFB power value is superior. The disadvantage of the statistical features lies in that they depict vibration signals from overall range, which introduces insensitive factors to the degradation process.

In conclusion, the power value on the SFB is typically a good indicator which is used to predict the RUL of bearing in this work.

2.2. The Generalized Weibull Failure Rate Function

The generalized WFRF is derived from Weibull distribution and failure rate function. Weibull distribution is one of the most commonly used models in reliability assessment for a mechanical component. It is flexible to represent different life distributions with different parameters. Failure rate function is the measurement which indicates the failure probability of a certain unit at a given time [21]. Defined by three parameters, the WFRF is mathematically represented as [32]where is the failure rate of equipment at time , is the scale parameter denoting characteristic life, is the shape parameter representing Weibull slope, and is the position parameter. It is obvious in a degradation process that increases from 0 to 1 as time goes. Therefore, is usually greater than 1 [33]. 2-parameter WFRF, which is widely used in bearing life prognosis, is derived by setting .

The generalized WFRF is the updated version of WFRF. The equation of the generalized WFRF is as follows [21, 32].where is the fitted measurement value and is the age of the equipment. and are introduced to adapt the equation to any situation.

2.3. Radial Basis Function Neural Network

RBF neural network as a feedforward network [34] is shown in Figure 2. It consists of three layers: an input layer with neurons, a hidden layer with neurons, and an output layer with neurons. The input values are transferred from the input layer to the hidden layer. The basis functions in the hidden layer will give out a response to the inputs and generate the outputs in the neurons of the output layer. The response is nonzero only when the input falls within a limited region of the input space. Hence RBF neural network is also regarded as a local range network [35, 36].

The transfer function of the hidden layer is a Gaussian function. The output value of the th neuron in the hidden layer is expressed as [34]where and are the central and width parameters of the th Gaussian function, respectively.

In RBF neural network, the hidden layer realizes the nonlinear mapping of while the output layer realizes the linear mapping of . The projection between input and output can be expressed as [37]where is the connection weight from hidden layer to output layer.

Although both RBF and BP neural networks can approximate any nonlinear function, RBF neural network has the advantage of a fast converging time without local minimum, since it only performs linear adjustment of the weights.

2.4. RUL Prediction Method

The framework of the proposed method is presented in Figure 3.

Firstly, SFB is determined by comparing FFT spectrums along the bearing degradation process. Secondly, the power value on the SFB is calculated for each FFT spectrum. As an example shown in Figure 4, although the trend of the power value is roughly increasing, there are large fluctuations appearing especially after the occurrence of an initial defect. This may be because of the degradation mechanism that the spall is rounded after grinding for many times [10]. However, it is obvious that the health condition of a bearing deteriorating with time is a monotonic process [38]. For the sake of consistency with degradation process, the generalized WFRF is applied to fit the obtained power value series and four parameters of the generalized WFRF, namely, , , , and , are derived.

Thereafter, RBF neural network is used as an intelligent system to predict the RUL of a bearing with the fitted measurements. In this case, four parameters including , , , and are selected as the inputs of the RBF neural network. and are the ages at current and previous inspections, respectively, while and are the fitted SFB power values at and . These four measurements are selected as inputs because they are crucial to represent the health condition of bearings. On the one hand, the fitted SFB power values at current and previous inspections can not only import the degradation status into RBF neural network but also introduce the change information of the measurements. On the other hand, we choose the fitted measurements at two inspections instead of more because it can give the ANN better generalization capability. As it comes to the output, the life percentage is preferred. The failure occurs when the life percentage reaches 100%. We denote the life percentage as at the time inspection of . Then the RUL of the bearing is derived by

3. Experimental Setup

The vibration signals used in this paper are provided by Center for Intelligent Maintenance Systems (IMS), University of Cincinnati [26]. Figure 5 presents the experimental setup of IMS.

As it is shown in Figure 5, four Rexnord ZA-2115 double row bearings were installed on the shaft which was driven by an AC motor via rub belts. The rotation speed of the shaft was kept constant at 2000 rpm. Through a spring mechanism, a radial load of 6000 lbs was exerted on the shaft and bearings. One PCB 353B33 High Sensitivity Quartz ICP accelerometer for each bearing was installed on the bearing housing to collect the vibration signal. With the sampling rate of 20 KHz, a set of data consisting of 20480 points was collected for one second every 10 minutes. Lubricant was driven by an oil circulation system to lubricate the test rig. The bearing degradation process was indicated by the debris adhered to a magnetic plug installed in the oil feedback pipe. The test would stop when the accumulated debris arrived to a preset amount. More detailed information about this experiment can be found in [39].

In this research, two vibration datasets are applied for the validation of the proposed RUL prediction method. Dataset 2 Bearing 1 of 984 vibration signals with an outer race failure is selected as an example to illustrate the proposed method in detail, while Dataset 1 Bearing 3 of 2156 vibration signals with an inner race defect is adopted to perform a comparative analysis. Figure 6 shows the failure bearing components after run-to-failure experiment [4].

4. Results and Discussion

The development of FFT spectrum for Dataset 2 Bearing 1 is shown in Figure 1. Intuitively, we can observe that the spectrum evolves gradually at the beginning and dramatically at the end. In order to identify the SFB mathematically, we compute the accumulated amplitude difference of FFT spectrums on each frequency line as follows.where is the amplitude of the th frequency line on the th FFT spectrum and is the number of FFT spectrums which is equal to 984. It is clear that can represent the variation of each frequency line throughout the lifetime. Based on we define an indicator , which is

It is the summation of , which stands for the total difference from the 1st frequency line to the th frequency line. If the differences of the corresponding frequency line amplitudes along with lifetime are small, the change of will be smooth. If such differences are large, the value of will increase rapidly. We would like to choose the frequency band as the SFB, on which the slope of is relatively larger.

In Figure 7, and are depicted. It is easy to identify that the slope of is relatively larger from 3000 Hz to 5000 Hz. In other words, the amplitudes of frequency lines evolve dramatically on this frequency band. Hence the SFB in this case is selected as [3000 Hz, 5000 Hz].

It should be noted that the amplitudes on the SFB hardly appear to be different until the degradation process of bearing reaches to certain degree. For Dataset 2 Bearing 1, the amplitudes on [3000 Hz, 5000 Hz] do not increase abruptly until the life percentage comes to nearly 60%. Nevertheless, this matter does not become a barrier to the proposed prediction method. In comparison to the early stage of the bearing’s lifetime, what counts more in the RUL prediction is the accurate prediction late in its life, because this directly affects whether and when to replace the bearing [21].

To depict the advantage of the proposed degradation indicator, we draw a comparison among the SFB power value and several most used features for diagnosis and prognosis of bearings defined, respectively, in Table 1.

In order to describe the signal from multiangles, the average amplitude of the defective frequency and its first six harmonics (AADFSH) is picked as the frequency domain feature [8, 25]. Furthermore, as it was proposed in [40], the vibration signal was decomposed by empirical mode decomposition (EMD) and the first three intrinsic mode function (IMF) energy datasets were selected. Together with the five time domain features presented in Table 1, the original feature set was constructed. Then PCA was conducted and the first principal component was chosen as an ultimate indicator.

To evaluate the foregoing features and select the most sensitive indicator, three goodness metrics, namely, correlation (Corr), monotonicity (Mon), and robustness (Rob) are utilized [41]. Before then, each of the feature series is decomposed into two parts.where is the degradation feature value at time , is the trend value, and is the residual value.

The feature goodness metrics are computed by the following equations [32, 4143]:where denotes the number of observations during the lifetime, and are the corresponding means of and , and is the difference between two successive points.

Corr is the measurement of linearity between feature and time. Mon evaluates the consistency of feature variation tendency and Rob reflects the tolerance of the feature to outliers. From the equations, we can see that all three goodness metrics are fallen into the range and are positively correlated with performance of features. In order to consider three metrics comprehensively, a weighted linear combination is proposed in [41] as the degradation indicator selection criteria:where , , and are the weights of the metrics, respectively, and determined by the contribution to the degradation process [41]. However, this fusion method by directly formulating on the metrics brings great distortion and thus leads to false selection of degradation indicator because the scales of these three metrics are not at the same level. See Figure 8. To eliminate this inequality, each of the metric dataset is normalized through before the criteria are computed.

The feature datasets obtained above are processed by (9)-(10), and the results of three metrics and criteria are tabulated in Table 2.

The characteristic frequency of outer race defect of Dataset 2 Bearing 1 used in AADFSH is determined in [44].

From Table 2, we can see the SFB power value has the largest Mon, which indicates that the tendency of this feature is closest to the degradation process. Moreover, the computing time of each feature dataset is also obtained. In this paper, a computer equipped with an Intel® E5-1620 processor CPU @ 3.60 GHz, 16 Gb RAM, and the 64-bit Windows 7 operating system (OS) was utilized. The feature PC1 took considerably more time than other features because EMD was conducted 984 times which required high computing burdens. This makes PC1 unsuitable for online prognostics. With the largest Cri and relatively less running time, the SFB power value is the most sensitive feature and selected as degradation indicator in this paper.

Figure 4 presents the SFB power value series of Dataset 2 Bearing 1, in which 1 data point denotes 10 minutes and this stipulation is valid throughout this paper. The generalized WFRF is adopted to fit the measurement line. Table 3 shows the results of the estimated parameters computed by nonlinear least squares method. To evaluate the fitting performance, Figure 9 shows the actual SFB power value series and the fitted series. It is obvious that the fitted series gives a better indication of the bearing degradation.

According to the vibration signals of Dataset 2 Bearing 1, there are 984 inspections throughout the lifetime. Due to two fitted measurements at current and previous inspections combining as one input of RBF neural network, 983 pairs of data are obtained. In order to verify the effectiveness of the proposed method, we segment these 983 datasets into two categories: a training set and a test set. 500 pairs of data are chosen randomly from 983 pairs as training set and 100 pairs are selected arbitrarily from the remaining 483 pairs as test set. Then RBF neural network is trained utilizing the training set. During the training process, the inherent configure of the RBF neural network is adjusted to minimize the mean square error (MSE) between the network outputs and the actual life percentages.where is the number of data, is the network output, and is the actual life percentage.

As the training process completed, the test set is fed into the trained network to testify its performance. The training and test performance is illustrated in Figures 10(a) and 10(b). The actual life percentage values are represented by “+”, and the predicted life percentage values are denoted by “blue circle”. The MSE for the test set in this paper is , which is smaller than both training and validation MSE in [7]. From Figure 10, we can see that the actual percentage values and the predicted values coincide almost at each selected data point, which indicates the good performance of the RBF neural network.

The proposed method has been compared with other methods in previous research. The performance of prediction can be simply calculated aswhere is the actual RUL and is the predicted RUL.

Adopting the same run-to-failure datasets collected from IMS, the results are shown in Table 4. In order to compare the performance fairly, the same failure threshold is utilized in the corresponding reference.

It is clear that the accuracy of the proposed method is higher, which indicates a good performance of this method in predicting the RUL of bearings.

5. Conclusion

Accurate RUL prediction of bearing is crucial to perform PHM of rotating machinery for preventing fatal failure and reducing maintenance cost. Most literature concentrating on RUL prediction mainly adopts two steps: indicator selection and life prediction model construction. This paper proposes a novel degradation indicator, namely, the SFB power value; it is the direct reflection of bearing health condition. The advantage of this indicator is verified by comparing with several commonly used features. Subsequently, the SFB power value series along the lifetime is fitted with the generalized WFRF to keep consistent with the degradation process. RBF neural network is utilized to perform bearing RUL prediction with the ages and the fitted power values at present and previous inspections as inputs and the life percentage as output.

In the experimental verification, three prognostic methods were compared with the proposed method. The results showed that the proposed method provided a more accurate prediction, which demonstrated the effectiveness and superiority of this method in predicting the RUL of bearings. Moreover, the proposed method is suitable for practical application and industrial field, because the novel indicator in this paper is easy to obtain and robust to fault categories.

Conflicts of Interest

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

Acknowledgments

This work is supported by National Key R&D Program of China (2016YFC0802908), National Natural Science Foundation of China (no. 51475455), and the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).