A Monotonic Degradation Assessment Index of Rolling Bearings Using Fuzzy Support Vector Data Description and Running Time

Performance degradation assessment based on condition monitoring plays an important role in ensuring reliable operation of equipment, reducing production downtime and saving maintenance costs, yet performance degradation has strong fuzziness, and the dynamic information is random and fuzzy, making it a challenge how to assess the fuzzy bearing performance degradation. This study proposes a monotonic degradation assessment index of rolling bearings using fuzzy support vector data description (FSVDD) and running time. FSVDD constructs the fuzzy-monitoring coefficient ε̄ which is sensitive to the initial defect and stably increases as faults develop. Moreover, the parameter ε̄ describes the accelerating relationships between the damage development and running time. However, the index ε̄ with an oscillating trend disagrees with the irreversible damage development. The running time is introduced to form a monotonic index, namely damage severity index (DSI). DSI inherits all advantages of ε̄ and overcomes its disadvantage. A run-to-failure test is carried out to validate the performance of the proposed method. The results show that DSI reflects the growth of the damages with running time perfectly.


Introduction
Rolling bearings, as important components of rotating machinery, not only support the load but also allow relative motion [1][2][3][4]. They are also a common failure unit due to their complex running conditions. The performance of bearings directly affects the operation reliability of the whole equipment [5][6][7][8][9], therefore, fault diagnosis and degradation assessment based on condition monitoring have been a key means to ensure the reliable operation of equipment, reduce the production downtime and save maintenance costs, etc. [10].
Vibration analysis is a powerful tool for fault diagnosis and degradation assessment [11][12][13][14]. The performance of rolling bearings is influenced by many factors, such as rotation speed, temperature, and lubrication conditions [15,16]. The running state is the final results of those factors. The vibrations reflect synthetically the present running state of the bearing. Many researches have been developed about the feature extraction from the vibration signals to assess the performance degradation. Time-domain features and frequency-domain features are the first choice because of the easy calculation and the definite physical meaning. Root mean square (RMS), Kurtosis Factor, the average amplitude of the defective frequency and its first six harmonics were successively used in [17][18][19], but the original features have different advantages and limitations as well. For example, RMS is a good stability feature which steadily grows with the fault development of the bearing. However, it is difficult to discover the incipient defects by RMS. On the contrary, Kurtosis Factor is sensitive to impulse faults. It distinctly appears at the initial fault stage while it decreases with the fault development. Kurtosis Factor has high sensitivity for incipient faults, but the poor stability for the serious damage. Few original features satisfy the conditions of sensitivity and stability simultaneously. A composite index is necessary which is both sensitive to the initial defect and rises stably as the damage grows.
One challenge is how to structure an intelligent assessment model based on the original features. Several scholars have proposed some comprehensive indexes and obtained impressive results. Qiu et al. proposed the minimum quantization error (MQE) for the assessment of the performance degradation of rolling bearings based on Self Organizing Map (SOM) and optimal wavelet filters [20]. Huang et al. applied MQE to the residual life for ball bearings using a back propagation neural network [21]. Pan et al. developed three indicators spectral entropy, health index and degradation indicator using information entropy, wavelet packet-support vector data description and lifting wavelet packet decomposition-fussy c-means, respectively [22][23][24]. Yu utilized a locality preserving projection for feature extraction, and provided the two complementary indexes, exponential weighted moving average (EWMA) and negative log likelihood probability-based EWMA statistic (NLLPEWMA) to assess the bearing performance degradation [25,26]. Caesarendra et al. proposed the combination of relevance vector machine (RVM) and logistic regression (LR) for the performance degradation assessment and prediction from incipient defects to final failure [27,28]. These indicators reflect the degradation trends of bearing performance to different degrees. However, sometimes a monotonic index is necessary to reflect the damage development more accurately. Little research about monotonic indexes is reported. The monotonic index should be constructed for two reasons. Firstly, the damage is irreversibly growing with the operation time. If the assessment index is oscillating, it is unable to reflect the damage development trends perfectly. Secondly, the damage shows accelerated growth with the running time. Few indicators reflect these accelerated relationships.
This study proposes a monotonic degradation assessment index for rolling bearings using fuzzy support vector data description (FSVDD) and running time. FSVDD is the combination of fuzzy mathematics theory and SVDD, and deals well with the fuzzy matter in small samples. The performance degradation of bearings belongs to this situation. The performance degradation with the strong fuzziness is the intermediate process between the normal running and the final failure. The initial defect, the final failure and the degrees of damage severity with different moments are hard to identify. Meanwhile, the dynamic information which reflects the change of bearing states is fuzzy and random. Some deviations appear when the fuzzy matter is dealt with by the deterministic mathematic method. Besides, the bearing monitored during its whole life is rare due to the difficulties in engineering. Support vector data description (SVDD) is an excellent method of one-class classification for small samples, with the advantages of robustness and high computation [29,30]. However, SVDD without the fuzzy identification function is unable to discriminate the damage severity degrees of samples. FSVDD is the fuzzy computing method based on SVDD. It introduces the fuzzy membership degree into the kernel function [31,32]. The fuzzy membership degree indicates the importance according to the damage severity degree. FSVDD constructs the fuzzy-monitoring membership which describes the accelerating relationships between the damages and running time. In addition, the running time is used to construct the monotonous growth of the bearing damage development.
This paper is organized as follows: In Section 2, the basic theory of FSVDD is introduced, and Section 3 presents the performance assessment method based on FSVDD and running time. The bearing run-to-failure tests and the related analysis are provided in Section 4. Section 5 provides the conclusions from the above studies.

Support Vector Data Description
Support vector data description proposed by Tax and Duin [29] is an excellent method of one-class classification. It is inspired by the idea of the support vector machine machines put forth by Vapnik [33]. The basis idea of SVDD is to find a spherically shaped decision boundary with minimum volume containing all (or most) targets, as described in Figure 1.  Assuming a training set contains n vectors of objects {x i , i = 1,2,…,n}. We try to find the minimum-volume hypersphere which contains all or most normal objects. This hypersphere is described by center d and radius R, and satisfies the following function: (1) where ξ i is slack variable to enhance the robustness and C is a penalty constant which gives the trade-off between volume of the hypersphere and errors. Then constraint condition is incorporated into the objective function by Lagrange multipliers: (2) where α i ≥ 0, γ i ≥ 0 are the Lagrange multipliers. L should be minimized with respecting to R, d, ξ i . Taking the partial derivatives of L and R, d, ξ i to zero, new constraints are given as follows: All α i are obtained by solving Equation (5) and only a small set of them are not zero. The objects with α i > 0 are the support vectors. The radius R is acquired by any support vector x k : For a new object z, its distance to the center d is described as follows: (7) when the data in the input space are not linearly predicted, a non-linear vector function (x) is needed to map them onto a high-dimensional feature space. Consequently, the kernel function, K(x i ,x j )  x i x  , is brought in SVM which maps the original data points to a high-dimensional feature space and transforms the nonlinear problem to a linear model. When the kernel function is introduced, Equations (5)-(7) are transformed to the following form: Then the monitoring coefficient of new object z to the hypersphere boundary is given by: (11) If ε ≤ 0, z is target, or else it is an outlier.

Fuzzy Support Vector Data Description
SVDD only simply identifies the normal samples and the fault samples, but is unable to accurately distinguish the samples of different degrees of damage severity. Then FSVDD is generated by introducing the fuzzy mathematics theory into SVDD to describe the development process from the initial defect to final failure [30][31][32].
A fuzzy membership degree s i , 0 ≤ s i ≤ 1, is introduced with each training sample x i , which is regarded as the importance according to the damage severity degree. The bigger the fuzzy membership degree s i is, the more important the training sample x i becomes. Suppose we are given a training set S of labeled training samples with the fuzzy membership degrees: The fuzzy non-linear vector function is redefined as , and kernel function is also recounted as Equation (12): Equations (8)- (11) are transformed to the following form: where is the fuzzy-monitoring coefficient.
From the comparison of Equations (8)- (11) and Equations (13)-(15), we find that they are same when s i = 0. Hence, SVDD could be seen as a special type of FSVDD.

Degradation Assessment Method Based on FSVDD and Running Time
SVM，SVDD and FSVDD are excellent classifiers for small samples, and have been widely applied to fault diagnosis [30][31][32]34,35]. However, few studies have reported the performance degradation assessment based on them, because the degradation assessment is different from fault diagnosis. The purpose of fault diagnosis is to identify the bearing faults, while the performance degradation assessment is to find an index which truly reflects the damage growth process. The performance degradation assessment is the development and the extension of fault diagnosis. In this study, FSVDD is applied to construct a degradation assessment indicator with the running time.
The bearing damage undergoes sustained development with the running time after the initial defect. The damage development trend should be an increasing function of time. Hence we consider that the running time is used to construct a degradation assessment indicator which could reflect the damage development of the bearing more accurately. On the other hand, the vibration signal of bearings usually contains lots of random information which may cause some deviations between the signal and the true damage. It is able to reduce the influence of the vibration randomness that the running time is introduced into the construction of the assessment indicator.
The performance degradation of the bearing is a continuous process of change. A new bearing is installed into the rotating shaft. After a short running-in, it enters a long-term stable working period. Then the minor fault appears, and the defect gradually increases as the bearing fault develops. Finally the bearing fails with a serious defect. According to the running process, the bearing life is divided into three parts: normal stage, degeneration stage and failure stage. In the beginning of the normal stage, there may be a short run-in period. The degradation assessment model is provided based FSVDD and running time as follows: Step 1: Feature extraction and selection: The time-domain features of vibration signals are extracted. The stability features, such as RMS, square-root amplitude (SRA) and absolute average values (AAV), and the sensitive feature as Kurtosis factor are selected as the inputs of FSVDD in the meantime.
Step 2: Training and testing based on FSVDD The function of fuzzy membership degree s i associated with each sample x i is given by a rising ridge distribution such as Equation (17) due to the rising statistic features, the fuzzy bounds of degradation beginning and final failure: (17) where E(x i ) is the energy of x i , a 1 is the maximum energy of normal state and a 2 is the minimum energy of failure state.
The normal samples are trained to construct the hypersphere. Then the fuzzy-monitoring coefficient for each testing sample z is obtained by testing all samples. The fuzzy matter, degradation assessment, is made clear by the fuzzy-monitoring coefficient .
Step 3: Degradation assessment A monotonic degradation assessment index of rolling bearing, damage severity index (DSI), is described based on FSVDD and running time as follows: (18) where i is the sample number, t i is its corresponding running time, t id is the initial defect time and DSI i is the damage severity index of the present moment.
Two assumptions are necessary for Equation (18). Firstly, the bearing has already used up the run-in period, and moved into the latter stage. Secondly, the bearing runs continuously and has a steady working condition, that is, the fuzzy-monitoring coefficient does not have wild points. The performance degradation begins at the initial defect time, t id , which is found by the constraint > 0. The non-equality > 0 means that its corresponding sample is an outlier. If the bearing is still in the run-in period, the fuzzy-monitoring coefficient may be more than zero, but the corresponding sample is not the initial defect. After the run-in period, one point corresponds to > 0, but its neighbor points are all < 0. The point is taken for the wild point. Only when the points with > 0 continuously appear three times, the first point among them is considered as the initial defect time. After the initial defect, the wild points of the parameter should also be rejected. The deletion of the wild point assures that DSI gets rid of the influence of the vibration randomness. If the amplification of one point is more than 150% compared with the amplitude of former point except the initial defect, the point is considered as the wild point and its amplitude should be replaced by the mean amplitude of the neighbor points.
In Equation (18), the fuzzy-monitoring coefficient and the running time are used to construct the novel indicator DSI which has at least two meanings. On the one hand, DSI i = when is bigger, which means the emergence of new fault. On the other hand, the defect of the current moment is not less than that of the former moment, even if no new defect appears, because the damage development is irreversible. The running time makes DSI increase slowly when the running state of the bearing is balanced. DSI inherits the defuzzification ability as good as the fuzzy-monitoring coefficient , and means the quantitative measurement of the fault severity degree. A higher DSI value indicates more serious faults.

Bearing Run-to-Failure Test
The bearing run-to-failure test is carried out under constant load conditions on the bearing tester to reflect the defect propagation processes. The test bench shown in Figure 2 is made up of power and drive system, main body, hydraulic loading system, lubrication system, control system and an It is designed as a simply supported beam structure, as described in Figure 3.  The two test bearings, 30311 tapered roller bearings, are installed on both ends of the shaft, while the two steady bearings, N312 cylindrical roller bearings, are fixed at the middle of the shaft. The axial load F a is directly inflicted on bearing 4, and transmitted to bearing 1 by the shaft. The radial load F r is exerted on the steady bearings, transferred to the shaft, and finally inflicted on the test bearings. Table 1 lists the parameters of test bearings and steady bearings. The vibration signals are acquired with a sampling frequency of 10 kHz per channel by YE6267 dynamic data collection and analysis Data acquisition computer system system by Sinocera Piezotronics Inc. The data recorder is equipped with low-pass filters at the input stage for anti-aliasing. Each sample with 32,768 points is collected every five min. The acceleration sensor, as shown in Figure 4, is the Lance LC0401 High Sensitivity ICP accelerometer.  The vibration signals are transmitted by the screw which touches the outer-race of bearing. The spring exerts a pre-tightening force to ensure the contact of the screw and the outer-race. The acceleration sensor is fixed on the screw with the insulation spacer that insulates electromagnetic interference. In addition, the temperatures of four bearings are monitored by the thermocouple sensors. The axial load F a and the radial load F r of single bearing are 15 kN and 27 kN, respectively. The rotation speed is kept constant at 1,500 r/min. The characteristic frequencies are calculated by Equations (19) where f r , f c , f o , f i , f e are the shaft frequency, the cage defect frequency, the outer-race defect frequency, the inner-race defect frequency and the rolling element defect frequency, respectively. d, D p , are the rolling element diameter and pitch diameter. is the number of rolling elements and is the contact angle. Besides, the ball spin frequency is f b = f e . Five characteristic frequencies are shown in Table 2.

Results and Discussion
There are three successful experiments in the present study. In each test, there is only one failure bearing which has huge vibration signals overwhelming those of the other three normal bearings. It may result from the individual factors of each bearing. The test 1 failure bearing has inner-race defects while the tests 2 and 3 failure bearings exhibit rolling element defects.
Firstly, the original time-domain features are researched. Three stability features, such as RMS, square-root amplitude (SRA), absolute average values (AAV), and one sensitive feature as Kurtosis factor of three bearings are listed in Figures 5-7. In each Figure, the right subgraph is the local enlargement of the left subgraph to describe the variance of each feature in degradation period more distinctly. These figures at least tell us the following: (1) The normal periods usually are obviously longer than the degradation period which is verified in reference [36]. (2) In the normal period, the three stability features have placid trends while they grow continuously with the development of faults in the degradation period, but there is not an obvious impulse for each stability feature when the initial defect occurs. Three features of test 2 failure bearing have large changes at 7,600 min. It could be resulted from a dismounting and reinstallation. (3) Kurtosis factor is bumping up when the incipient fault appears, but its subsequent behaviors are bad. Kurtosis factor could be used to roughly discriminate the initial defect times, that is, the beginning moments of degradation, which are 10,000 min, 10,710 min and 3,075 min, respectively. However, the beginning moments of degradation may be unfaithful because of the randomness of Kurtosis factor, especially when the impulse in Kurtosis factor is very weak, such as in the test 3 failure bearing. (4) The same feature for different bearings varies greatly because of the individual differences, even though at the same period. For examples, RMS of test 1 failure bearing is clearly less than that of tests 2 and 3 failure bearings at the normal period. (5) The overall trend of the original features is fuzzy, while their every point has the strong randomness. Therefore, none of the original features are suitable to assess the performance degradation over the whole lifetime.   Secondly, the vector x i is constructed by the four original features and imported into SVDD to compute the monitoring coefficient ε. Figure 8 provides the parameter ε and its local enlargement of three bearings. The monitoring coefficient ε is a comprehensive index which is sensitive to the initial defect and steadily increases with the damage development. The monitoring coefficient ε is utilized to find the initial defect time. The degradation beginning threshold is defined as ε > 0. The samples with ε > 0 are outside of the sphere acquired by training the normal samples. The two assumptions of Equation (18) are necessary in the processes of the initial defect time determined by the threshold ε > 0. They assure  The left subgraphs are the time domain waveforms at different moments, and the right subgraghs are their corresponding Hilbert spectra. The moment of 4,000 min is the normal stage for the test 1 failure bearing, and the fault characteristic frequencies are not seen in Figure 9(b). A characteristic frequency of 234.7 Hz close to the inner-race defect frequency appears in Figure 9(d), and gradually increases along with the damage development in Figure 9(f,h). Therefore, the moment of 9,860 min could be more suitable as the initial defect time for the test 1 failure bearing. The two initial defect times of test 2 failure bearing are equivalent. The rolling element defect frequency and its second harmonic appear at 10,710 min in Figure 10(d). The contrast experiment of test 3 failure bearing provides the same results as that of test 1 failure bearing. The time determined by ε i > 0, 3,040 min, is more suitable as the degradation beginning time. However, the monitoring coefficient ε is not the ideal indicator. First of all, the parameter is oscillating, while the actual damages are irreversible. Then the monitoring coefficient ε increases slowly with time, which is not able to reflect the accelerated relationship between the damage development and the running time perfectly. That could be as a result of the fuzziness of the damage quantitative. SVDD might need an improvement to deal with the fuzzy damage development.
Thirdly, the fuzzy membership degree s i is computed by Equation (17). The fuzzy-monitoring coefficients of three bearings are given by FSVDD. The parameter and its local enlargement are described in Figure 12. The fuzzy-monitoring coefficient is an improvement of the monitoring coefficient ε, and adds the function of defuzzification. The parameters and ε have a certain degree as well as lots of differences. The comparisons between ε and are carried out as shown in Figure 13. The parameter ε is blue while the parameter is green. They are same at the beginning of degradation, and much different with the damage development. The overall trend of − ε is accelerated, but the acceleration is not strict and not able to be proved by the second derivative. The fuzzy-monitoring coefficient has a similar accelerated trend, which could agree with the relationship of the damage development and running time. On the other hand, the D-value of the neighbor gets bigger as the damage increases. The increasing D-value means that the damage severity of the neighbor moments is differentiated easily. However, the fuzzy-monitoring coefficient is still oscillating and not consistent with the irreversible damage development.  Finally, a new index, DSI, is calculated as Equation (18), and Figure 14 shows the DSI and its local enlargement. The wild points of fuzzy-monitoring coefficient should be rejected before the computation of DSI. The parameter in Figure 12 contains lots of exceptional points. The elimination of wild points assures that DSI does not suffer from vibration randomness. DSI inherits all advantages of the fuzzy-monitoring coefficient , and overcomes its shortcomings. The monotonic index, DSI, reflects the increase of the bearing damages with running time perfectly. Sometimes DSI is jumping which means the occurrence of spalling or the emergence of new damages. On the other hand, the placid DSI implies that the running state of bearing is balanceable and no new damage occurs. The degradation beginning threshold is defined as DSI > 0 according to the monitoring coefficient ε. The initial defect occurs at the degradation beginning moment. And the confirmation of the failure threshold should consider the running state of the bearing and the synthetic performance of the whole rotary machine. In this study, the failure threshold is delimited as DSI ≤ 4 The failure moments of three bearings are 12,260 min, 11,385 min and 4,265 min, as shown in Figure 14.

Outer-Race Defect and Inner-Race Defect Simulation
In the run-to-failure tests, there are only two faults types, the inner-race defect and the rolling element defect. Then the outer-race defect and inner-race defect are simulated by increasing the impulse amplitude of the mathematical model to verify again the effectiveness of DSI.
The rolling bearing transfers the main load through elements in rolling contract than in sliding contract [37]. When the defect appears, the characteristic frequencies are seen easily. McFadden developed a model to simulate the vibrations of single point defect [38]. The model contains the effects of rolling element bearing geometry, shaft speed, load distribution, transfer function, and the exponential decay vibration. Then it is extended to describe the vibrations of multiple point defects [39]. The above models provide the cyclical shock, and have been successfully utilized to describe faults in rolling element bearings [40][41][42][43]: (23) where T is the period of impulse, s(t) is the vibration waveform and n(t) is external noise, A i is the amplitude modulator to simulate the possible modulated situation, τ i is the random fluctuation around average period T.
Simplifying the vibration waveform s(t) to an exponential damping cosine signal: (24) where B is the an appropriate value to simulate the attenuation of oscillation waveforms, f n is the natural frequency related to bearing or system. The amplitude modulator, A i , is simplified as a cosine signal: where A 0 is the resonance intensity, f m = 1/Q is the shaft speed for inner-race fault and the cage speed for rolling element fault, and Q is the modulated period. Let f m = 0 Hz, T = 1/f i = 0.061 s to simulate the outer-race defect and let f m = f r = 25 Hz, T = 1/f i = 0.042 s to simulate the inner-race defect. The natural frequency f n = 2,000 Hz. Fifty samples are simulated for each defect. Each simulation sample contains 32,768 points. The initial defects of the outer-race defect and inner-race defect are at about the 40th sample. The original features and comprehensive indexes of outer-race defect are shown in Figure 15. The initial defect is considered as the 41st sample by Kurtosis factor. The randomness of original features is big due to the strong noise in the time-domain waveform. The new initial defect location is determined at the 37th sample by the threshold ε > 0 and earlier than that determined by Kurtosis factor. The failure of outer-race defect appears at the 48th sample by the threshold DSI ≤ 4. The similar rule is seen for the inner-race defect. Kurtosis factor determines the initial defect at 41st sample, while the fuzzy-monitoring coefficient confirms the time at 39th sample. The final failure of inner-race defect is confirmed at 48th sample by DSI. Besides, the parameters ε and both have many wild points as result of the randomness of original features. DSI is an excellent index which effectively reflects the irreversible development of the bearing defect. In a word, DSI is an excellent degradation assessment index for rolling bearings and has at least the following advantages: (1) DSI is sensitive to the initial defect and grows stably with the development of faults. The stability features, such as RMS, SRA, AAV, reflect the damage development, but it is hard to find the initial defect. The sensitivity feature as Kurtosis Factor is the opposite. DSI is an excellent indicator which is increasing with the damage development and sensitive to the initial defect. Moreover, the run-to-failure experiment and simulation both verify that DSI determines the initial defect earlier.
(2) DSI with the defuzzification ability reflects the accelerating relation between the damage development and running time. This advantage of DSI comes from the fuzzy-monitoring coefficient . The parameter describes the accelerated development of bearing damage as a result of the fuzzy membership degree s i . However the acceleration is not strict and unable to be verified by the second derivative. (3) The monotonic DSI reflects the irreversible development of the bearing defects. The damage development is irreversibly increasing after the initial defect occurs. The monotonic indicator is more effective to reflect the degradation process.
The type and width parameter of the kernel function affect the parameters ε and , even DSI, as described in reference [24]. In this study, Gaussian Raial Basis Function (GRBF) function was adopted due to its good property and universality. The width parameter was calculated carefully by the ex-ante fault diagnosis experiments. At beginning, we selected some normal samples and several of fault samples from the vibration of each failure bearing. The range of width parameters was estimated by the rules of thumb [44]. In the range, the width parameter was sampled with the same interval. The classification error with each sample of the width parameter was calculated by the classifier. The best value range was decided by the minimum error. The width parameter was selected as width = 2 which is in the overlay region of the best ranges for the three failure bearings. The selection of type and width parameter is worth further study for the better degradation assessment.

Conclusions
This paper presents a study of degradation assessment based on FSVDD and the running time. SVDD constructs the monitoring coefficient ε with the advantages of sensitivity to the initial defect and stable growth as the damages develop, but SVDD could not deal well with fuzzy damage severity degree. FSVDD introduces the fuzzy membership degree s i to SVDD and provides a better assessment index . It could describe the accelerated relation between the development of the damages and running time. However, the index with a oscillating trend disagrees with the actual damage development. Finally a monotonic degradation assessment index, DSI, is constructed. DSI inherits all the advantages of and reflects the irreversible growth of the bearing damages with running time well.
Bearing run-to-failure tests and simulation experiments were carried out to validate the proposed method. The results show that DSI reflects the running state of rolling bearing in real time and effectively guarantees the operation reliability of bearings. Although the analysis results in this study are acceptable, more tests are necessary for further analysis. Meanwhile, other mathematical methods, such as Principal Component Analysis (PCA), can be used to construct some new indexes. PCA uses an orthogonal transformation to convert a set of possibly correlated variables into linearly uncorrelated variables. The transformation is defined in such a way that the first principal component has the largest possible variance, and each succeeding component in turn has the highest possible variance under the constraint which is orthogonal to the preceding components. If the statistical features are input into PCA, many separate composite characteristics are obtained. Some of the composite characteristics may be excellent assessment indexes.
In addition, DSI could be applied to the performance degradation assessment of other key machine components, such as gears, shafts and ball screws, because their dynamic performance is similar to that of bearings. The remaining life prediction could be carried out by the combination of DSI and the regression methods such as SVM, ANN etc. DSI is used to determine the time of the initial defect and the finial failure, and the regression methods predict the remaining life by inputting one or many fine indexes like DSI. It could take all advantages of DSI and the regression methods for accurate life prediction.