A Monotonic Degradation Assessment Index of Rolling Bearings Using Fuzzy Support Vector Data Description and Running Time Academic research paper on "Mechanical engineering"

Sensors 2012, 12, 10109-10135; doi:10.3390/s120810109

OPEN ACCESS

sensors

ISSN 1424-8220

www.mdpi.com/j ournal/sensors

Article

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

Zhongjie Shen, Zhengjia He *, Xuefeng Chen, Chuang Sun and Zhiwen Liu

State Key Laboratory for Manufacturing System Engineering, School of Mechanical Engineering, Xi'an Jiaotong University, Xi'an 710049, China; E-Mails: zjshen.2007@stu.xjtu.edu.cn (Z.S.); chenxf@mail.xjtu.edu.cn (X.C.); sunchuang522@163.com (C.S.); Lzw1682007@126.com (Z.L.)

* Author to whom correspondence should be addressed; E-Mail: hzj@mail.xjtu.edu.cn; Tel./Fax: +86-29-8266-3689.

Received: 6 June 2012; in revisedform: 4 July 2012 /Accepted: 9 July 2012 / Published: 26 July 2012

Abstract: 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.

Keywords: performance degradation assessment; fuzzy support vector data description; running time; run-to-failure test; operation reliability

1. Introduction

Rolling bearings, as important components of rotating machinery, not only support the load but also allow relative motion [1-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-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-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-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-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.

2. Fuzzy Support Vector Data Description

2.1. 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.

Figure 1. Schematic diagram of support vector data description.

hyper sphere (boundary)

target

Assuming a training set contains n vectors of objects (x,, i = 1,2,...,«}. 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:

min L(R,d,£) = R2 + C,

st f(x -d)TX -d)< R2 + £, (1)

. [ £ > 0, i = 1,2, —, n

where £ 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:

L{R, d,a, ,6 ,rt ) = R2 + C Z 6-± at [r2 + 6 -(x2 - 2dx, + d2 )]-£ y¿t (2)

i=i i=i i=i

where ai > 0, yi > 0 are the Lagrange multipliers. L should be minimized with respecting to R, d, Taking the partial derivatives of L and R, d, £ to zero, new constraints are given as follows:

= 1, d = Z^iXi, y, = C-ai (3)

i=1 i=1

Equation (3) deduces:

0<a <C (4)

Reconstituting Equation (2) by Equation (3) results in:

max L(a) = 2 a(x- • X-) - 2 aaj(x- • xj)

i=1 i=1, j=1

s.t. <

Za = i,

0 <a< C

All ai are obtained by solving Equation (5) and only a small set of them are not zero. The objects with ai > 0 are the support vectors. The radius R is acquired by any support vector xk:

j2 n n

R =(xk ■ xk) - 2Z a(x- • xk) + Z aaj(x- • xj) (6)

i=1 i, j =1

For a new object z, its distance to the center dis described as follows:

2 II i. 2 n n

Rz HIz - 4 =(z •z) - 2z a(z • x, ) + Z aaj(X- • xj) (7)

i =1 i, j =1

when the data in the input space are not linearly predicted, a non-linear vector function (j(x) is needed to map them onto a high-dimensional feature space. Consequently, the kernel function, K(xi,xj) = (xi)(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:

max L(a) = S atK(x • x)- S aiajK(X ' X)

¿=1 i=1, j=1

s.t. <

S a¿ = 1,

0 <a< C

ry n n

R2 = K(xk • xk) -2ea¿k(x¿ • xk) + s a^jK(x¿ • x.) (9)

i=1 i, j=1

o II n2 n n

r2 =\\z - d|| = K(z • z) - 2 Sa,K(z • x¿) + S a1ajK(xi ■ xj )

i = 1 i, j =1

Then the monitoring coefficient of new object z to the hypersphere boundary is given by:

s = (Rz - R)/ R (11)

If e < 0, z is target, or else it is an outlier. 2.2. 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-32].

A fuzzy membership degree si, 0 < si <1, is introduced with each training sample xi, which is regarded as the importance according to the damage severity degree. The bigger the fuzzy membership degree si is, the more important the training sample xi becomes. Suppose we are given a training set S of labeled training samples with the fuzzy membership degrees:

S = {(x1, s1(x2, s2 " •, (xn, sn )}

2, s2)

The fuzzy non-linear vector function is redefined as ^(x, s ) = (s + 1)^(x), and kernel function is also recounted as Equation (12):

K(x¿ , xj )=(s, + 1)^(x, )'(sj + ifaj )=(s, + + 1)k(x¿ , xj ) (12)

Equations (8)-(11) are transformed to the following form:

max L(a) = S a, (s, +1)2 K(x¿ • x¿ ) - S a'aj (s' +1) (sj +1)K(x¿ ' xi)

i=1 i=1, j=1

S ai = 1,

0 < a < C

s.t. <

R2 =(sk +1)2 K(xk ■ xk ) - 2(sk + 1)S a (s¿ + 1)k(x¿ • xk ) + S a^ (s¿ + ^ + 1)k(x¿ • xj) (14)

i=1 i,j=1

R2 = (s +1)2K(z • z) - 2(s + 1)S a (s¿ + 1)K(z • x¿ ) + S aaj (s¿ + ^ + 1)k(xt • x.)

i=1 i, j=1

e = (Rz - R)/ R (16)

where is the fuzzy-monitoring coefficient.

From the comparison of Equations (8)—(11) and Equations (13)—(15), we find that they are same when Sj = 0. Hence, SVDD could be seen as a special type of FSVDD.

3. 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-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 st associated with each sample xt 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:

where E(xj) is the energy of xu a1 is the maximum energy of normal state and a2 is the minimum energy of failure state.

0< e (x)<a

E (X )> a1

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:

where i is the sample number, t, is its corresponding running time, tid is the initial defect time and DSI 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, tid, 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, DSIi = £ 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.

4. Experiment and Simulation Verification

4.1. Bearing Run-to-Failure Test

DSL =\

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

independent data acquisition system, etc. It is designed as a simply supported beam structure, as described in Figure 3.

Figure 2. Test bench for bearing run-to-failure.

Figure 3. Load diagram of four bearings.

Radial load Fr

Axial load Fa

Bearing 1 Bearing 2 Bearing 3 Bearing 4

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 Fa is directly inflicted on bearing 4, and transmitted to bearing 1 by the shaft. The radial load Fr 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

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.

Table 1. Parameters of test bearings and steady bearings.

Bearing Inner Outer Roller Roller Calculating Coefficient Basic

Diameter Diameter Diameter Dynamic

Type Number e Y

(mm) (mm) (mm) Load (kN)

30311 55 120 16.25 16 0.35 1.7 152

N312 60 130 19.1 16 — — 212

Figure 4. Location of accelerometer sensors and temperature sensors.

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 Fa and the radial load Fr 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)-(22):

f — fr Jc ~ 2

1--cos 0

f= -Jo 2

f= ~ f 2

t d a 1--cos 0

v D J i d

1 + — cos0

where fr, fc, fo, fi, fe 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, Dp, are the rolling element diameter and pitch diameter. n is the number of rolling elements and 0 is the contact angle. Besides, the ball spin frequency is fb = fe. Five characteristic frequencies are shown in Table 2.

4.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.

Table 2. Bearing characteristic frequencies.

fr (Hz) fc (Hz) fo (Hz) f (Hz) fe (Hz)

25 10.2393 163.8288 236.1712 65.1061

Figure 5. Four time-domain features of test 1 failure bearing: (a) RMS during the whole life; (b) local enlargement of RMS; (c) SRA during the whole life; (d) local enlargement of SRA; (e) AAV during the whole life; (f) local enlargement of AAV; (g) Kurtosis factor during the whole life; (h) local enlargement of Kurtosis factor.

Figure 6. Four time-domain features of test 2 failure bearing: (a) RMS during the whole life; (b) local enlargement of RMS; (c) SRA during the whole life; (d) local enlargement of SRA; (e) AAV during the whole life; (f) local enlargement of AAV; (g) Kurtosis factor during the whole life; (h) local enlargement of Kurtosis factor.

Figure 7. Four time-domain features of test 3 failure bearing: (a) RMS during the whole life; (b) local enlargement of RMS; (c) SRA during the whole life; (d) local enlargement of SRA; (e) AAV during the whole life; (f) local enlargement of AAV; (g) Kurtosis factor during the whole life; (h) local enlargement of Kurtosis factor.

Secondly, the vector xt is constructed by the four original features and imported into SVDD to compute the monitoring coefficient e. Figure 8 provides the parameter e and its local enlargement of three bearings.

Figure 8. Results of SVDD for tests 1-3 failure bearing: (a) e of test 1 failure bearing during the whole life; (b) local enlargement of e for test 1 failure bearing; (c) e of test 2 failure bearing during the whole life; (d) local enlargement of e for test 2 failure bearing; (e) e of test 3 failure bearing during the whole life; (f) local enlargement of e for test 3 failure bearing.

The monitoring coefficient e is a comprehensive index which is sensitive to the initial defect and steadily increases with the damage development. The monitoring coefficient e is utilized to find the initial defect time. The degradation beginning threshold is defined as e > 0. The samples with e > 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 e > 0. They assure

that the initial defect time is free from the influence of the vibration randomness, so the moments of 0 ~ 220 min of test 1 failure bearing are not the degradation beginning time though e, > 0, because of the bearing in the run-in period. Likewise, the moments of 5,585 min, 7,340 min, 7,830 min are not the initial defect time because they are the exceptional points. Their e, > 0 but their neighbor points are all e, < 0. The similar situations appear in the tests 2 and 3 failure bearings. The degradation beginning times of three bearings are determined as 9,860 min, 10,710 min and 3,040 min with et > 0. The new initial defect moments are 40 min, 0 min and 35 min earlier than that determined by Kurtosis factor. Which initial defect time is more accurate? The contrast experiments shown in Figures 9-11 are carried out to verify them.

Figure 9. The Time-domain waveforms and Hilbert spectrums of test 1 failure bearing: (a) Time-domain waveform on 4,000 min; (b) Hilbert spectrum on 4,000 min; (c) Time-domain waveform on 9,860 min; (d) Hilbert spectrum on 9,860 min; (e) Time-domain waveform on 10,000 min; (f) Hilbert spectrum on 10,000 min; (g) Time-domain waveform on 12,250 min; (h) Hilbert spectrum on 12,250 min.

-1 oj-

Time(minutes)

Time(minutes)

Time(minutes)

tu> 0.08

a 0.06

tp 0.04

(d) 0.12

i> 0 08

S 0.04

i> 0.1

< 0.05

1000 1500 2000

Frequency(Hz)

X:234.7 Y:0.0893

500 1000 1500

Frequency(Hz)

X:234.7 Y:0.120

500 1000 1500

Frequency(Hz)

90898948239089009190892391908901919190892348908948919089009190900000919089239190895391179089489117894891908953919089894891908900919089234890898923919089539190895391905353912390895348908900239089489191488989919089539190234891

Figure 9. Cont.

153 1°-

H S ■

< mL -i° -

1 1.5 2

Time(minutes)

5°° 1°°° 15°° 2°°° 25°°

Frequency(Hz)

Figure 10. The Time-domain waveforms and Hilbert spectrums of test 2 failure bearing: (a) Time-domain waveform on 4,000 min; (b) Hilbert spectrum on 4,000 min; (c) Time-domain waveform on 10,710 min; (d) Hilbert spectrum on 10,710 min; (e) Time-domain waveform on 11,400 min; (f) Hilbert spectrum on 11,400 min.

ü Ü.5 1 1.5 2 2.5 3

Time(minutes)

° °.5 1 1.5 2 2.5 3

Time(minutes)

15 F 1°

1 °1111

° °.5 1 1.5 2 2.5 3

Time(minutes)

°.15r

<1> °.1

B °.°5

3 °.15

fp S °.1

< °.°5

5°° 1°°° 15°° 2°°°

Frequency(Hz)

25°°

5°° 1°°° 15°° 2°°°

Frequency(Hz)

X:13°.6 Y:°.21°

25°°

X:65.31 Y:°.168

1°°° 15°°

Frequency(Hz)

2°°°

25°°

Figure 11. The Time-domain waveforms and Hilbert spectrums of test 3 failure bearing: (a) Time-domain waveform on 1,000 min; (b) Hilbert spectrum on 1,000 min; (c) Time-domain waveform on 3,040 min; (d) Hilbert spectrum on 3,040 min; (e) Time-domain waveform on 3,075 min; (f) Hilbert spectrum on 3,075 min; (g) Time-domain waveform on 4,685 min; (h) Hilbert spectrum on 4,685 min.

ig CD T3

Time(minutes)

ig <D T3

i— 500

1000 1500 2000

Frequency(Hz)

ig 1> T3

Time(minutes)

1 1.5 2 2.5 3

Time(minutes)

500 1000 1500 2000

Frequency(Hz)

Y:0.1781

~^X:65.31 Y:0.1276

1000 1500 2000

Frequency(Hz)

0.5 1 1.5 2 2.5 3

Time(minutes)

1000 1500

Frequency(Hz)

"7UPPP!"!7U UPPU""+"!UPU7U PU!77P PU7" PU7UU!7P77P P7"7!UP7UU7"77""UUP PUUUU PPPP

01910223482301238900029123905389480002028900239153895348910190488948539023892348904889230191319100895348914889530191489053230148904889002391489002022600900189235390020300029101890153912390005391239002895300530289004891238948539100902389004890488900

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 £j > 0, 3,040 min, is more suitable as the degradation beginning time. However, the monitoring coefficient e is not the ideal indicator. First of all, the parameter is oscillating, while the actual damages are irreversible. Then the monitoring coefficient e 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 st is computed by Equation (17). The fuzzy-monitoring coefficients e of three bearings are given by FSVDD. The parameter e and its local enlargement are described in Figure 12. The fuzzy-monitoring coefficient is an improvement of the monitoring coefficient e, and adds the function of defuzzification. The parameters e and e have a certain degree as well as lots of differences. The comparisons between e and e are carried out as shown in Figure 13. The parameter e 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 e - e 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 Ej 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.

Figure 12. Results of FSVDD for tests 1-3 failure bearing: (a) e of test 1 failure bearing during the whole life; (b) local enlargement of e for test 1 failure bearing; (c) e of test 2 failure bearing during the whole life; (d) local enlargement of e for test 2 failure bearing; (e) e of test 3 failure bearing during the whole life; (f) local enlargement of e for test 3 failure bearing.

(C) 5 4 3

I CO 2 1 0 -1

Figure 12. Cont. (d) e,

0 2000 4000 6000 8000 10000

Time(minutes)

1.05 1.1

Time(minutes)

0 1000 2000 3000 4000

Time(minutes)

X:3040 Y:0.222

3000 3500 4000

Time(minutes)

Figure 13. Comparisons between s and e: (a) s & e of test 1 failure bearing; (b) e - s of test 1 failure bearing; (c) s & e of test 2 failure bearing; (d) e - s of test 2 failure bearing; (e) s & e of test 3 failure bearing; (f) e - s of test 3 failure bearing.

1.1 1.15

Time(minutes)

1 -0 --1 -

1.1 1.15

Time(minutes)

1.08 1.09 1.1 1.11 1.12 1.13 1.14 1.15

. .. «4

Time(minutes)

1.08 1.09 1.1 1.11 1.12 1.13 1.14 1.15

Time(minutes)

I to 4

Figure 13. Cont. (f)

3500 4000

Time(minutes)

I to 1 0

3500 4000 4500

Time(minutes)

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 e should be rejected before the computation of DSI. The parameter e 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 e, 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 e. 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.

Figure 14. Results of degradation assessment for tests 1-3 failure bearing: (a) DSI of test 1 failure bearing during the whole life; (b) local enlargement of DSI for test 1 failure bearing; (c) DSI of test 2 failure bearing during the whole life; (d) local enlargement of DSI for test 2 failure bearing; (e) DSI of test 3 failure bearing during the whole life; (f) local enlargement of DSI for test 3 failure bearing.

0 2000 4000 6000 8000 10000 1

Time(minutes)

X:12255

X:9860 J

- Y:0.132 -

1.05 1.1 1.15 1.2

Time(minutes) x 104

Figrue 14. Cont.

4.3. 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-43]:

x(t) = Z A ■ s(t-iT-Tt) + n(t) (23)

where T is the period of impulse, s(t) is the vibration waveform and n(t) is external noise, At is the amplitude modulator to simulate the possible modulated situation, t is the random fluctuation around average period T.

Simplifying the vibration waveform s(t) to an exponential damping cosine signal:

s (t) = ^Bt •cos (2f + ) (24)

where B is the an appropriate value to simulate the attenuation of oscillation waveforms, fn is the natural frequency related to bearing or system. The amplitude modulator, At, is simplified as a cosine signal:

4 = 4 • cos(2fffmt + yA ) (25)

where A0 is the resonance intensity, fm = 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 fm = 0 Hz, T = 1f = 0.061 s to simulate the outer-race defect and let fm = fr = 25 Hz, T = Hft = 0.042 s to simulate the inner-race defect. The natural frequency fn = 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 e > 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 e and e 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.

Figure 15. Original features and comprehensive indexes of outer-race defect: (a) RMS; (b) SRA; (c) AAV; (d) Kurtosis factor; (e) e; (f) e; (g) DSI.

Y:0.02

Figure 15. Cont. (f)

Number

Y:0.02 - 2

Number

X:48 _ Y:5.618

X:37 Y:0.02 '

Number

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,. 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 e and e, 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.

Figure 16. Original features and comprehensive indexes of inner-race defect: (a) RMS; (b) SRA; (c) AAV; (d) Kurtosis factor; (e) e; (f) e; (g) DSI.

0.9 0.85 0.8

0.8 0.75 0.7 0.65

Number

Number

Number

Number

0.65 -

< 0.6-

0.5 L 50 0

— 2.8 50 0

X:39 -

Y:0.01

I to 2 0

X:48 '

Y:4.99 ^

Y:0.01 -

Number

Number

Number

X:41 '

Y:3.09

VVvA-n^V -

Y:0.01

5. Conclusions

This paper presents a study of degradation assessment based on FSVDD and the running time. SVDD constructs the monitoring coefficient e 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 Sj 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 e 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 e 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.

Acknowledgements

This work is supported by National Natural Science Foundation of China (No. 51035007 and No. 51175401), and Doctoral Fund of Ministry of Education of China (No. 20110201130001), and Fork YING TUNG Education Foundation (No. 121052).

References

1. Tse, P.W.; Peng, Y.H.; Yam, R. Wavelet analysis and envelope detection for rolling element bearing fault diagnosis-their effectiveness and flexibilities. J. Vib. Acoust. 2001, 123, 303-310.

2. Yan, R.Q.; Gao, R.X. A nonlinear noise reduction approach to vibration analysis for bearing health diagnosis. J. Comput. Nonlinear Dynam. 2011, 7, 1555-1423.

3. Wang, Y.X.; Liang, M. Identification of multiple transient faults based on the adaptive spectral kurtosis method. J. Sound Vib. 2012, 331, 470-486.

4. Gao, L.X.; Yang, Z.J.; Cai, L.G.; Wang, H.Q.; Chen, P. Roller bearing fault diagnosis based on nonlinear redundant lifting wavelet packet analysis. Sensors 2011, 11, 260-277.

5. Pusca, R.; Romary, R.; Ceban, A.; Bridny, J.F. An online universal diagnosis procedure using two external flux sensors applied to the AC electrical rotating machines. Sensors 2010, 10, 10448-10466.

6. Guo, L.X.; Zai, F.L.; Su, S.B.; Wang, H.Q.; Chen, P.; Liu, L.M. Study and application of acoustic emission testing in fault diagnosis of low-speed heavy-duty gears. Sensors 2011, 11, 599-611.

7. Gui, L.L.; Kang, C.H.; Wang, H.Q.; Chen, P. Application of composite dictionary multi-atom matching in gear fault diagnosis. Sensors 2011, 11, 5981-6002.

8. Lei, Y.G.; Lin, J.; He, Z.J.; Kong, D. A method based on multi-sensor data fusion for fault detection of planetary gearboxes. Sensors 2012, 12, 2005-2017.

9. Wang, S.B.; Huang, W.G.; Zhu, Z.K. Transient modeling and parameter identification based on wavelet and correlation filtering for rotating machine fault diagnosis. Mech. Syst. Signal Process. 2011, 25, 1299-1320.

10. Du, S.C.; Lv, J.; Xi, L.F. Degradation process prediction for rotational machinery based on hybrid intelligent model. Robot. Comput.-Integr. Manuf. 2012, 28, 190-207.

11. Malhi, A.; Yan, R.Q.; Gao, R.X. Prognosis of defect propagation based on recurrent neural networks. IEEE Trans. Instrum. Meas. 2011, 60, 703-711.

12. Lei, Y.G.; He, Z.J.; Zi, Y.Y.; Hu Q. Fault diagnosis of rotating machinery based on multiple ANFIS combination with Gas. Mech. Syst. Signal Process. 2007, 21, 2280-2294.

13. Wang, H.Q.; Chen, P. A feature extraction method based on information theory for fault diagnosis of reciprocating machinery. Sensors 2009, 9, 2415-2436.

14. Gao, L.X.; Ren, Z.Q.; Tang, W.L.; Wang, H.Q.; Chen, P. Intelligent Gearbox diagnosis methods based on SVM, Wavelet Lifting and RBR. Sensors 2010, 10, 4602-4621.

15. Qiu, H.; Lee, J.; Lin, J.; Yu, G. Wavelet filter-based weak signature detection method and its application on rolling bearing prognostics. J. Sound Vib. 2006, 289, 1066-1090.

16. Feng, Z.P.; Wang, T.J.; Zuo, M.J.; Chu, F.L.; Yan, S.Z. Teager energy spectrum for fault diagnosis of rolling element bearings. J. Phys.: Confer. Ser. 2011, 305, 012129:1-012129:7.

17. Gebraeel, N.; Lawley, M.; Liu, R.; Parmeshwaran, V. Residual life predictions from vibration-based degradation signals: A neral network approach. IEEE Trans. Ind. Electron. 2004, 51, 694-700.

18. Gebraeel, N.; Lawley, M. A neural network degradation model for computing and updating residual life distributions. IEEE Trans. Autom. Sci. Eng. 2008, 5, 154-163.

19. Shao, Y.; Nezu, K. Prognosis of remaining bearing life using neural networks. Proc. Inst. Mech. Eng. I-J. Syst. C. 2000, 214, 217-230.

20. Qiu, H.; Lee, J.; Lin, J.; Yu, G. Robust performance degradation assessment methods for enhanced rolling element bearing prognostics. Adv. Eng. Inform. 2003, 17, 127-140.

21. Huang, R.Q.; Xi, L.F.; Li, X.L.; Liu, C.R.; Qiu, H.; Lee, J. Residual life predictions for ball bearings based on self-organizing map and back propagation neural network methods. Mech. Syst. Signal Process. 2007, 21, 193-207.

22. Pan, Y.; Chen, J.; Li, X.L. Spectral entropy: A complementary index for rolling element bearing performance degradation assessment. J. Mech. Eng. Syst. 2009, 223, 1223-1231.

23. Pan, Y.; Chen, J.; Li, X.L. Bearing performance degradation assessment based on lifting wavelet packet decomposition and fuzzy c-means. Mech. Syst. Signal Process. 2010, 24, 559-566.

24. Pan, Y.; Chen, J.; Guo, L. Robust bearing performance degradation assessment method based on improved wavelet packet-support vector data description. Mech. Syst. Signal Process. 2009, 23, 669-681

25. Yu, J.B. Bearing performance degradation assessment using locality preserving projections, Expert Syst. Appl. 2011, 38, 7440-7450.

26. Yu, J.B. Bearing performance degradation assessment using locality preserving projections and Gaussian mixture models. Mech. Syst. Signal Process. 2011, 25, 2573-2588.

27. Caesarendra, W.; Widodo, A.; Yang, B.S. Application of relevance vector machine and logistic regression for machine degradation assessment. Mech. Syst. Signal Process. 2010, 24, 1161-1171.

28. Caesarendra, W.; Widodo, A.; Thom, P.H.; Yang, B.S.; Setiawan, J.D. Combined probability approach and indirect data-driven method for bearing degradation prognostics. IEEE Trans. Rel. 2011, 60, 14-20.

29. Tax, D.M.J.; Duin, R.P.W. Support vector domain description. Pattern Recognit. Lett. 1999, 20, 1191-1199.

30. Wang, D.; Tse, P.W.; Guo, W.; Miao, Q. Support vector data description for fusion of multiple health indicators for enhancing gearbox fault diagnosis and prognosis. Meas. Sci. Technol. 2011, 22, 1 -14.

31. Forghani, Y.; Sadoghi Yazdi, H.; Effati, S. An extension to fuzzy support vector data description (FSVDD*). Pattern Anal. Appl. 2011, 21, 1-11.

32. Luo, H.; Wang, Y.R.; Cui, J. A SVDD approach of fuzzy classification for analog circuit fault diagnosis with FWT as preprocessor. Expert Syst. Appl. 2011, 38, 10554-10561.

33. Vapnik, V.N. Statistical Learning Theory; Wiley: New York, NY, USA, 1998.

34. Hao, R.J.; Peng, Z.K.; Feng, Z.P.; Chu, F.L. Application of support vector machine based on pattern spectrum entropy in fault diagnostics of rolling element bearings. Meas. Sci. Technol. 2011, 22, 1-14.

35. Zhang, Y.; Liu, X.D.; Xie, F.D.; Li, K.Q. Fault classifier of rotating machinery based on weighted support vector data description. Expert Syst. Appl. 2009, 36, 7928-7932.

36. Williams, T.; Ribadeneira, X.; Billington, S.; Kurfess, T. Rolling element bearing diagnostics in run-to-failure lifetime testing. Mech. Syst. Signal Process. 2001, 15, 979-993.

37. Shigley, J.E.; Mischke, C.R.; Budynas, R.G. Mechanical Engineering Design; McGraw-Hill: New York, NY, USA, 2004.

38. McFadden, P.D.; Smith, J.D. Model for the vibration produced by a single point defect in a rolling element bearing. J. Sound Vib. 1984, 96, 69-82.

39. McFadden, P.D.; Smith, J.D. The vibration produced by multiple point defects in a rolling element bearing. J. Sound Vib. 1985, 98, 263-273.

40. Antoni, J.; Bonnardot, F.; Raad, A.; EI Badaoui, M. Cycloststionary modeling of rotating machine vibration signals. Mech. Syst. Signal Process. 2004, 18, 1285-1314.

41. Antoni, J.; Randall, R.B. Differential diagnosis of gear and bearing faults. J. Vib. Acoust. 2002, 124, 165-171.

42. Cong, F.; Chen, J.; Pan Y. Kolmogorov-Smirnov test for rolling bearing performance degradation assessment and prognosis. J. Vib. Control 2010, 17, 1337-1347.

43. Fan, B.; Hu, N.; Hu, L.; Gu, F. Application of phase space warping on damage tracking for bearing fault. J. Phys.: Confer. Ser. 2012, 364, 1-8.

44. Cherkassky, V.; Ma, Y. Practical selection of SVM parameters and noise estimation for SVM regression. Neural Netw. 2004, 17, 113-126.

© 2012 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Copyright of Sensors (14248220) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.