Next Article in Journal
Quantification of a Low-Cost Stretchable Conductive Sensor Using an Expansion/Contraction Simulator Machine: A Step towards Validation of a Noninvasive Cardiac and Respiration Monitoring Prototype
Next Article in Special Issue
Rotating Electrical Machine Condition Monitoring Automation—A Review
Previous Article in Journal
Root Cause Identification of Machining Error Based on Statistical Process Control and Fault Diagnosis of Machine Tools
Previous Article in Special Issue
Drilling Rig Hoisting Platform Security Monitoring System Design and Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Review of Feature Extraction Methods in Vibration-Based Condition Monitoring and Its Application for Degradation Trend Estimation of Low-Speed Slew Bearing

by
Wahyu Caesarendra
1,2 and
Tegoeh Tjahjowidodo
3,*
1
Mechanical Engineering Department, Diponegoro University, Semarang 50275, Indonesia
2
School of Mechanical, Materials and Mechatronic Engineering, University ofWollongong, Wollongong, NSW 2522, Australia
3
School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798, Singapore
*
Author to whom correspondence should be addressed.
Machines 2017, 5(4), 21; https://doi.org/10.3390/machines5040021
Submission received: 9 August 2017 / Revised: 16 September 2017 / Accepted: 19 September 2017 / Published: 27 September 2017
(This article belongs to the Special Issue Machinery Condition Monitoring and Industrial Analytics)

Abstract

:
This paper presents an empirical study of feature extraction methods for the application of low-speed slew bearing condition monitoring. The aim of the study is to find the proper features that represent the degradation condition of slew bearing rotating at very low speed (≈1 r/min) with naturally defect. The literature study of existing research, related to feature extraction methods or algorithms in a wide range of applications such as vibration analysis, time series analysis and bio-medical signal processing, is discussed. Some features are applied in vibration slew bearing data acquired from laboratory tests. The selected features such as impulse factor, margin factor, approximate entropy and largest Lyapunov exponent (LLE) show obvious changes in bearing condition from normal condition to final failure.

1. Introduction

Slew bearings are large roller element bearings that are used in heavy industries such as mining and steel milling. They are particularly used in turntables, cranes, cranes, rotatable trolleys, excavators, reclaimers, stackers, swing shovels and ladle cars [1]. In this paper, the slew bearing being analyzed is the one typically used in BlueScope Steel Australia; specifically for the coal loading into the steel ladle turret, coal stacker and reclaimer. Slew bearing is different to typical rolling element bearing in terms of geometry, load and speed application. Slew bearings are usually operated at very low rotational speed, from 0.5 to 15 rpm in either continuous or reversible rotation and are designed with three axes to support high axial and radial loads. The bearings have large diameter (e.g., approximately from 200 mm to 5000 mm). Due to the high geometric and cost issue, industries purchase these bearings for replacement in a regular maintenance period. However, in some cases, these bearings may fail before the predicted time. To prevent catastrophic failure and to effectively schedule the bearing purchase for minimizing the replacement cost, a condition monitoring method for early failure detection is required.
Analyzing the vibration signatures of slew bearings that operate at a very low speed is more complex than that of typical rolling element bearings. Slew bearing signals have low energy. They are also non-stationary, non-linear, chaotic and practically hard to analyze. Due to their low energy, the signal is masked by a higher background noise level. Therefore, it is difficult to extract the pertinent bearing signal. It is necessary to have reliable condition monitoring parameters or features that could extract accurate information about the bearing condition from the signal.
The number of studies on slew bearing condition monitoring and early failure detection based on vibration analysis is lower than for rolling element bearing. Studies of slew bearing condition monitoring and life prediction utilizing vibration signal are presented in [2,3,4,5,6,7]. Žvokelj et al. [2] presents the condition monitoring method utilizing vibration signal based on ensemble empirical mode decomposition (EEMD) and principal component analysis (PCA). An extension work of [2] is presented in [3] which used combined EEMD with kernel PCA. Recent studies of the slew bearing fault detection and diagnosis is presented by Žvokelj et al. [4]. The paper presents combined EEMG and multiscale independent component analysis (ICA). Nevertheless, studies of slew bearing failure prediction and prognosis also increased recently [5,6,7]. Feng et al. [5] develops a reliability model for residual life prediction based on Weibull analysis. A method for degradation trend estimation of slew bearings was presented in [6]. The method employed PCA, particle swarm optimization (PSO) and least squares support vector machine (LSSVM). Another residual life prediction model study of slewing bearing can be found in [7]. The authors used PCA and support vector regression (SVR) to construct the model.
In order to monitor the condition of slew bearings, a feature that represents the trends related to degradation condition from normal to failure is necessary. Due to the unique vibration signal generated by the roller and raceway contact in slew bearing, a suitable and sensitive feature extraction method is mandatory in slew bearing fault diagnosis and prediction. There are numerous existing methods, which are divided into different categories. For instance, Henao et al. [8] presented diagnostic methods in a non-stationary condition. The methods are classified into four approaches: (1) frequency domain approach; (2) time domain analysis; (3) slip-frequency; (4) diagnosis in the time-frequency domain.
This paper particularly presents an empirical study of feature extraction methods for vibration signals acquired from naturally degraded slew bearings operating at very low rotational speed. In terms of feature extraction methods, six categories are reviewed, namely (i) time-domain features extraction; (ii) frequency-domain features extraction; (iii) time-frequency representation; (iv) phase-space dissimilarity measurement; (v) complexity measurement; and (vi) other features. Among the six categories, the time-domain features extraction method is found to be the most dominant method in typical rolling element bearings. In addition, the third and fourth methods are suitable for non-stationary, non-linear and chaotic signals.

2. Laboratory Slew Bearing Experiment

2.1. Slew Bearing Test-Rig and Sensor Location

The vibration data employed in this paper were acquired from a slew bearing test rig [1]. The test rig was designed to replicate real working conditions of a slew bearing in a steel company, such as a high-applied load and low rotational speed. The slew bearing used was a large thrust bearing with an inner and outer diameter of 260 mm and 385 mm, respectively. The rotational speed of the lab test rig can be operated from 1 to 12 r/min. A photo of the test rig and sensor location on the slew bearing is shown in Figure 1.

2.2. Data Acquisition Procedure

The vibration data for continuous rotation was collected from four accelerometers that were installed on the inner diameter surface at 90° to each other with 4880 Hz sampling rate. The accelerometers were IMI 608A11 ICP (IMI Sensors, PCB Piezotronics Division, New York, NY, USA) type sensors. The accelerometers were connected to a high-speed PicoScope DAQ (PS3424) (Pico Technology, St Neots, UK). The bearing was subjected to an axial load of 15 tons. The bearing began running in 2006 with rotational speed of 1 r/min. However, to provide continuous monitoring and produce run-to-failure bearing data, the bearing data was collected on a daily basis from February to August 2007 with total number of 139 days. To accelerate the bearing service life, coal dust was injected into the bearing in mid-April 2007 (58 days from the beginning). In practice, especially in steel-making companies, the slew bearing is located in open air where the bearing is exposed to a dusty environment and for this reason the coal dust was inserted in mid-April to simulate the real working conditions. The schematic of the vibration data acquisition from slew bearing test-rig is presented in Figure 2.

3. Features Extraction Methods and Its Application on Slew Bearing Vibration Signal

3.1. Category 1: Time-Domain Features Extraction

3.1.1. Statistical Time-Domain Features

Statistical time-domain features such as mean, root mean square (RMS), standard deviation and variance were usually used in past studies to identify the differences between one vibration signal and another. More advanced statistical-based features such as skewness and kurtosis can be applied to the signal which is not purely stationary. These features examine the probability density function (PDF) of the signal. It is a well-known fact that if the condition of the bearing changes, the PDF also changes, thus the skewness and kurtosis might also be affected. In particular, skewness is used to measure whether the signal is negatively or positively skewed, while kurtosis measures the peak value of the PDF and indicates if the signal is impulse in nature. For a signal with a normal distribution i.e., normal bearing signal has a skewness value of zero.
In addition, kurtosis is obtained from the peak of the PDF of the vibration signal, while skewness is obtained from the mean value of the PDF of the vibration signal. It is well known that the kurtosis value of vibration signal from normal bearing is approximately three [9] and the skewness value of approximately zero. Then, when the PDF of vibration signal changes due to faults, the kurtosis value will increase to greater than three and the skewness value will shift to either negative or positive. Other features which can be calculated from PDF of the vibration signal includes: entropy, which calculates the histogram of the PDF and measures the degree of randomness of the vibration signal; lower bound and upper bound histogram, which measure the lower and upper values of the PDF respectively.
Aside from the statistical time-domain features mentioned above, there are other non-dimensional features such as shape factor (defines as RMS divided by mean) and crest factor (defines as standard deviation divided by RMS). Both features will change as the mean, RMS and standard deviation change. The summary of the statistical time-domain features is presented in Table 1.

3.1.2. Upper and Lower Bound of Histogram

The application of the histograms or a discrete PDF for time-domain feature extraction is presented in [13]. Let d be the number of division that group the ranges, and h i with 0 i < d be the columns of the histogram for the time x i
h i = j = 0 n 1 n r i ( x i ) , i , 0 i < d
r i = { 1 , i f i ( max ( x i ) min ( x i ) ) d x < ( i + 1 ) ( max ( x i ) min ( x i ) ) d 0 , o t h e r w i s e
The upper bound h U and lower bound h L of histogram are defined as
h U = max ( x i ) + Δ / 2
h L = max ( x i ) Δ / 2
where Δ = ( max ( x i ) min ( x i ) ) / ( n 1 ) .
According to the features extraction results presented in Figure 3, some features, such as RMS, kurtosis, histogram upper bound and histogram lower bound, show the fluctuation in the last measurement day, approximately from the 90th day to the 139th day. These indicate that some features are sensitive to the slew bearing condition while others are less. The incipient damage occurred in day 90th that is clearly indicated by RMS, kurtosis, histogram upper bound and histogram lower bound feature, where the fluctuation is associated with roller and outer race defects. However, more appropriate and reliable features are necessary for the condition monitoring of slew bearing to detect an impending damage and to estimate the degradation index.

3.1.3. Autoregressive (AR) Coefficients

In other studies [13], autoregressive (AR) coefficients of the vibration signal are calculated as bearing vibration features. It is known that a faulty vibration signal of typical rolling element bearing produces different autoregressive coefficients compared to that of normal vibration signal. The following is a common univariate time series model to calculate the AR coefficients:
y t = a 1 y t 1 + a 2 y t 2 + + a n y t n + ε t = i = 1 n a i y t 1 + ε t
where a 1 to a n are the autoregressive coefficients, y t is the digitized vibration signal, n is the order of the AR model ( n = 8 ) and ε is the residual (Gaussian white noise).
The RMS, skewness, kurtosis, entropy and AR coefficients have been applied in vibration datasets from induction motor [14] and rolling element bearing [15]. These papers have made great contributions to the research areas on the bearing condition monitoring and fault diagnosis. However, the features were applied for artificially bearing damage. In practice, the damage is typically developed naturally rather than artificially.
The AR coefficients results are presented in Figure 4. It is shown that most of the features are less relevant to represent the slew bearing condition except for AR coefficient 4 which is shown fluctuating in the last measurement.
Other time domain features such as impulse factor (IF) and margin factor (MF) have been recently used [12,16,17]. Similar to CF, IF is used to measure how much impact is generated from the bearing defect. The CF measures the impact by dividing the maximum absolute value of vibration signal and RMS value, while IF divides the maximum absolute value by the mean of absolute value. The formula is presented as follows:
IF = max | x i | 1 N i = 1 N | x i |
Margin factor (MF) measures the level of impact between rolling element and raceway. The MF is calculated by dividing the maximum absolute value of vibration signal to the RMS of absolute value of vibration signal (7):
MF = max | x i | ( 1 N i = 1 N | x i | ) 2
The progression of slew bearing deterioration of IF and MF features are more apparent compared to RMS, kurtosis and histogram lower features as presented in Figure 5.

3.1.4. Hjorts’ Parameters

In addition to the time-domain features mentioned in previous subsections, Hjorts’ parameters also falls into this category [18]. These features are calculated based on the first and the second derivatives of the vibration signal. In the time series context, the numerical values for the derivatives are obtained as the differences between the current value and the prior value. The first Hjorth’s parameter is called ‘activity’. The activity of vibration signal can be computed by simply calculating the variance of the vibration signal amplitude σ x 2 , where σ x is the standard deviation of the sampled vibration signal x = ( x 1 , x 2 , x 3 , , x N ) .
The ‘mobility’ parameter of vibration signal is calculated as the square root of the ratio of the activity of the first derivative and the activity of the vibration signal.
mobility = σ x σ x
where σ x is the standard deviation of the first derivative of the vibration signal. The numerical values for the derivatives can be obtained as the differences between samples as follows.
x = x ( i + 1 ) x i        i = 1 , , N 1
The last parameter called ‘complexity’ which is calculated as the ratio of mobility of the first derivative and the mobility of the vibration signal.
complexity = σ x / σ x σ x / σ x
These parameters have been used in electroencephalography (EEG) signals to detect the epileptic seizures [19]. They have never been used in vibration bearing signal except for activity feature, which is similar to the variance feature in the statistical time-domain features extraction. It can be seen from Figure 6 that the progression of bearing condition is difficult to track using mobility and complexity feature. In contrast, the activity or variance feature shows one highest peak at approximately the 90th day and continue fluctuation on the last measurement days.

3.1.5. Mathematical Morphology (MM) Operators

In image processing, mathematical morphology (MM) was introduced as a non-linear method to analyze two dimensional (2D) image data including binary images and grey-level images based on set theory [20]. The fundamental principle of the MM method is to modify the shape of the original signal by transforming it through its intersection with another object called the ‘structuring element’ (SE). Four operators including erosion, dilation, closing and opening are used to achieve the transformation.
After having been applied successfully in 2D image processing, the MM method has been used in biomedical research field to analyze one dimensional (1D) EEG signals [21,22,23] and electrocardiography (ECG) signals [24]. The first known study of MM in 1D vibration signals was presented by Nikolaou and Antoniadis [25]. The authors used the MM method for vibration data with inner race and outer race faults. The effects of four MM operators on the enveloping impulse-type signals were examined. The result shows that the ‘closing’ operator is more suitable for extracting the impulse signal corresponding to the bearing fault type. The SE value was empirically investigated within the range of 0.2 T to 1.2 T. The optimum SE for a signal with additive noise was selected as 0.6 T (0.6 times the pulse repetition period). In addition, T is the period of each bearing component frequency such as BPFO, BPFI and BSF.
The original discrete one dimensional (1D) signal and structuring element (SE) must be first determined to use the four basic operators in the MM method. Suppose the discrete 1D vibration signal x = { x ( 1 ) , x ( 2 ) , , x ( N ) } , which functions over a domain x ( n ) = ( 1 , 2 , , N ) . Let S = { S ( 1 ) , S ( 2 ) , , S ( M ) } be the SE, which is the discrete function over a domain S ( m ) = ( 1 , 2 , , M ) , where N and M are integers and N M . Then, the four morphology operators for the discrete 1D vibration signal can be defined as follows:
  • Erosion: also refer to as min filter.
    ( x Θ S ) ( n ) = min ( x ( n + m ) S ( m ) )
  • Dilation: also refer to as max filter.
    ( x S ) ( n ) = max ( x ( n m ) + S ( m ) )
  • Closing: Dilates 1D signal and then erodes the dilated signal using the similar structuring element for both operations.
    ( x S ) ( n ) = ( ( x S ) Θ S ) ( n )
  • Opening: Erodes 1D signal and then dilates the eroded signal using the similar structuring element for both operations.
    ( x S ) ( n ) = ( ( x Θ S ) S ) ( n )
where n ( 1 , 2 , , N ) and m ( 1 , 2 , , M ) , the notations Θ , , , denote the erosion operator or Minkowski subtraction, dilation operator or Minkowski addition, closing operator and opening operator in the MM method, respectively.
The original simulated 1D signal is shown in the first row of Figure 7. The various effects of the basic MM operators application in the simulated 1D vibration signal are presented in the second to fifth row of Figure 7. The second row of Figure 7 to the fifth row of Figure 7 illustrate the effect of erosion, dilation, closing and opening, respectively. Please note that the solid blue line represents the original signal and the dashed red line is the signal after being processed or transformed by the MM operators. The simulated signal has zero mean value, which is typical to 1D vibration signals. The work principle of the basic MM operators can be summarized as follows: erosion of x ( n ) by the structuring element A ( m ) reduces the positive part and enlarges the negative part of x ( n ) , while dilation of x ( n ) by the structuring element A ( m ) reduces the negative part and enlarges the positive part of x ( n ) . In addition, closing of x ( n ) smoothens the signal x ( n ) by cutting down any impulses in the negative part and opening x ( n ) smoothens the signal x ( n ) by cutting down its impulses in the positive part.
In rolling element bearing, the MM method has been studied. The multi-scale morphological analysis was developed by Zhang et al. [26]. The method shows independency of empirical rules in terms of SE selection, but the technique is rather complex. Wang et al. [27] proposed an improved morphological filter as features of periodic impulse signals. An average weighted combination of open-closing and close-opening morphological operators was employed to eliminate the statistical deflection of amplitude and to extract the impulse component from the original signal. To optimize the SE, they used a new criterion called impulse attenuation (IMA). Dong et al. [28] presented the fault diagnosis of rolling element bearings using the modified morphological method. The morphological operator used in their method was the average weight combination of the closing and opening operator. They proposed a new criterion called SNR criterion to optimize the length of SE. Another recently published research on MM method was presented by Raj and Murali [29]. The authors combined the MM operators with fuzzy inference to classify the artificial bearing defects. The vibration data used in [29] were similar to the data used in Zhang et al. [26].
The application of the MM method in detecting the bearing fault frequencies from among impulsive-type signals have been demonstrated [26,27,28,29]. These signals are generated from known bearing conditions such as inner and/or outer race faults. It is acknowledged that the signals from artificial faults are easier to be identified than those from natural bearing faults. Besides, typically, the impulse-type signals are easily identified visually without any additional tool or processing step [28]. The application of MM method to the naturally damaged slew bearing is presented in Figure 8, Figure 9 and Figure 10.
The MM method has been used effectively for the damage detection of high speed rolling bearing. It is because when the damage occurred, the contact between the damage spot and rolling element generates detectable impulse due to the contact between the defect spot and race way or roller. In the case of naturally degraded low-speed slew bearing, the multiple damage e.g., outer race, inner race and roller are typically occurred at close sequence time and the magnitude of impulse is very low. Thus, it is difficult to identify the origin of impulse from outer race, inner race or roller damage. Therefore, the feature extraction result of MM for the three conditions: BPFO, BPFI and BSF present an identical figure.

3.2. Category 2: Frequency-Domain Features Extraction

To apply frequency-domain features, the time-domain vibration signals must initially be transformed into frequency-domain signals using fast Fourier transform (FFT). FFT is a common method in vibration analysis for bearing fault detection. FFT measures the dominant frequency of the repetitive impulse period of certain faults due to the contact between rolling elements and defective spot. Therefore, the type of bearing defects such as outer race, inner race or ball defect can be detected. FFT performed effectively in stationary periodic signals; however, it is less effective for non-stationary signals that arise from time-dependent events [30].

3.2.1. Statistical Frequency-Domain Features

Other frequency-domain features such as frequency centre (FC), root mean square frequency (RMSF) and root variance frequency (RVF) have been studied in [13]. When fault exists, the frequency element changes, and the values of the FC, RMSF and RVF also change. The FC and RMSF indicate the position changes of main frequencies, while the RVF shows the convergence of the power spectrum. Another feature called median frequency is sometimes used. The median frequency is measures by dividing the area of the amplitude spectrum into two equal part. In addition, the mean frequency (FC) can be defined as the first-order moment of the Fourier spectrum, RVF as the second-order moment, SS as the third-order moment and SK as the fourth-order moment of the Fourier spectrum [18]. These features can be calculated as follows:
FC = i = 2 N x i x i 2 π i = 1 N x i 2
MSF = i = 2 N ( x i ) 2 4 π 2 i = 1 N x i 2
RMSF = MSF
RVF = MSF FC 2
The result of FC, RMSF and RVF feature extraction is presented in Figure 11.

3.2.2. Spectral Skewness, Spectral Kurtosis, Spectral Entropy and Shannon Entropy Feature

Recently, new methods such as spectral skewness, spectral kurtosis, spectral entropy and Shannon entropy have been developed. Spectral skewness (SS) and spectral kurtosis (SK) are the advanced statistical measures applied to the magnitude spectrum. The SS measures the symmetry of the distribution of the spectral magnitude values around its mean [31]. It is defined by
SS ( n ) = 2 k = 0 Β L / 2 1 ( | X ( k , n ) | μ | X | ) 3 Β L σ | X | 3
The SK measures the distribution of the spectral magnitude values and compare to a Gaussian distribution. It is defined by
SK ( n ) = 2 k = 0 Β L / 2 1 ( | X ( k , n ) | μ | X | ) 4 Β L σ | X | 4 3
Antoni and Randall [32] presents the development of classical kurtosis analysis called spectral kurtosis (SK) to deal with the transient behaviour in a signal. SK can be used to detect incipient faults even in the presence of resilient masking noise. The concept of spectral statistic is adopted in this method to supplement the classical power spectral density (PSD). PSD ideally gives zeros values at those frequencies when the signal has stationary Gaussian noise, and it gives high positive values at those frequencies during the occurrence of the transients. When the noise-to-signal ratio is high, the transient signal will be buried in the background noise, thus the incipient fault cannot be clearly detected. SK can overcome this problem by analyzing the entire frequency band and to select the sensitive frequency band that correspond to the bearing condition.

3.3. Category 3: Time-Frequency Representation

Time-frequency domain representation methods such as short-time Fourier transform (STFT), wavelet transform and Wigner-Ville distribution (WVD) are commonly used for non-stationary or transient signal. These methods implement a mapping of one-dimensional time-domain signals to a two-dimensional function of time and frequency. The objective is to provide a true time-frequency representation of a signal. A review of time-frequency analysis methods for machinery fault diagnosis has been presented by Feng et al. [33].

3.3.1. Short-Time Fourier Transform (STFT)

Due to its low computational complexity and definite physical meaning, STFT is often used as an initial pre-processing tool to analyze non-stationary vibration signals. The STFT maps the vibration signal into a two-dimensional (2D) function of time and frequency [34]. STFT divides a non-stationary signal into small windows of equal time frame. The Fourier transformation is then applied to the time segment. Unlike FFT where the frequency transformation includes the whole original signal, in STFT original signal is decomposed using FFT at pre-defined time intervals. By doing this, the defect signal frequency can be identified in a particular window of time. The key to this method is the selection of the time interval. The STFT can be expressed as follows [35]:
STFT ( f , t ) = t T / 2 t + T / 2       w ( t τ ) x ( τ ) e j 2 π f τ d τ
where STFT ( f , t ) is short-term frequency spectrum, t is the time, f is frequency, w ( t τ ) is the shifted window, x ( τ ) denotes the input vibration signal, exp ( j 2 π f τ ) is the complex exponential, T is the window interval and w ( τ ) is the windowing function which satisfies w ( τ ) = 0 for | τ | = > T / 2 .
By taking the square magnitude of the STFT in (21), an energy density using the spectrogram S ( t , f ) method can be obtained
S ( t , f ) = | STFT ( t , f ) | 2
Equation (22) represents the energy density of the vibration signal x ( τ ) windowed by w ( τ ) .
The visualization of window length and hop size is presented in Figure 12. The vibration signal with the sample index i is split into overlapping blocks of length Β L and block index n . While, the hop size Η s is the distance from the start of one particular block and its subsequent block. The STFT result of one day slew bearing signal measurement is shown in Figure 13 for 10 s vibration signal with window length of 2048, hop size of 256 and number of FFT points of 16,384. It can be seen that during 10 s sampled signal, the amplitudes at lower frequencies were higher than the amplitude at high frequencies (approximately at 1300 Hz and 1750 Hz). This indicates that slew bearing frequencies were dominated by low frequency contents.

3.3.2. Wavelet Transform and Wavelet Decomposition

Wavelet transform is well-known as an appropriate method to analyze non-stationary and transient signal [14,36]. Wavelet transformation employs wavelet function instead of sinusoidal function as the basis function. It enhances a scale variable in addition to the time variable in the inner product transformation [33]. For any limited energy signal x ( t ) L 2 ( R ) , the wavelet function is defined as [26]
WT x ( t , a ) = 1 a + x ( τ ) ψ ( τ t a ) d τ
where wavelet ψ [ ( τ t ) / a ] is derived by dilating and translating the wavelet basis ψ ( t ) , a is the scale parameter ( a > 0 ) , t is the time shift and 1 / a is a normalization factor to maintain energy preservation.
A review of the wavelet transform applications in machine condition monitoring and fault diagnostics is presented by Peng et al. [37]. As a non-stationary signal analysis method, wavelet transform can be applied in machine fault diagnostics as a feature extraction method [37]. Since the wavelet transformations can provide an excellent energy concentration properties due to the uses of the wavelet basis function; the wavelet transform can represent the signal with a certain number of coefficients. The process is usually called wavelet decomposition.
Wavelet decomposition technique decomposes a non-stationary vibration signal into linear form of time-scale units. It decomposes and organizes original signal into several signal components according to the translation of the wavelet function. This changes the scale and shows the transition of each frequency component [38]. Through this decomposition, a selected coefficient can be used directly as the fault features. In wavelet decomposition, the input signal is basically decomposed into two coefficients. The input signal that is passed through the low-pass filter becomes the ‘approximation coefficient’. The input signal which is passed through the high-pass filter becomes the ‘detail coefficient.’ Since the low-frequency component is the most sensitive part to extract the useful information related to the bearing condition, the ‘approximate coefficient’ is passed to the next wavelet decomposition loop in multilevel wavelet decomposition process. Figure 14 shows three level wavelet decomposition structures. The structure illustrates the decomposition of slew bearing vibration signal into approximate coefficient (A) and detailed coefficient (D) at each level. For further processing, features such as mean, variance, skewness and kurtosis are calculated from the detailed coefficients of level 3 (D3) [39]. The features extraction results are presented in Figure 15, it can be seen that kurtosis can predict the damage compared to mean, variance and skewness calculated from signal D3 of the wavelet decomposition due to the significant changes after day 90 of the bearing life time.
The application of wavelet decomposition to the transient stator current signal acquired from induction motor is presented by Niu et al. [40]. Statistical features have also been used to reveal the fault information from the selected coefficients. Another study that uses the wavelet coefficients as the fault features for ball bearing is presented by Liu et al. [41].

3.3.3. Empirical Mode Decomposition-Based Hilbert Huang Transform

Empirical mode decomposition (EMD) [42] has been used in some applications that involve non-stationary signals (see, e.g., Braun and Feldman [43]). A review of EMD applications in fault diagnosis of rotating machinery is discussed in [44]. The review gives a detailed introduction to EMD application in various areas such as bearings, gears and rotors etc.
The EMD decomposes the vibration signal into several signals from high to low frequencies content, namely intrinsic mode functions (IMFs) based on the enveloping technique, i.e., Hilbert-Huang transform (HHT). From the specific frequencies content information, the bearing fault frequency can be identified. In some preliminary studies on condition monitoring [39,45], the original EMD [46] is used for non-stationary slew bearing data. A bearing fault such as outer race, inner race or rolling element fault has occurred when one of the IMF frequencies is identical to one of the bearing fault frequencies (as shown in Table 2) [45]. The EMD results for slew bearing data are presented in Table 3. A detail EMD result is presented in [39,45].
According to Table 3, the frequencies that correspond to the bearing fault condition did not appear on 24 February. In contrast, the fault frequency of rolling element (BSF axial) is identified on 3 May at 0.43 Hz. This frequency is associated to the roller defects as shown in Figure 16a. Interesting result was obtained from data on 30 August, where the BSF axial is hidden by BPFI axial of 1.37 Hz. The reason could be that the amplitude of signal 0.43 Hz is too weak and is covered by the background noise. The photo of the outer raceway damage is shown in Figure 16b. The roller and outer race damage photo was taken after final failure during the bearing dismantle on 2 September. The detailed vibration data used from the slew bearing test rig is presented in [39]. Statistical features from the decomposition results of EMD method is shown in Figure 17. The mean, variance, skewness and kurtosis features are calculated from the summation of the low IMF frequencies (from IMF 9 to the lowest IMF) of the EMD results. The selection of these IMFs is based on the characteristic of slew bearing signal being dominated by low frequency components as presented previously in STFT result.

3.3.4. Wigner-Ville Distribution (WVD)

The Wigner distribution (WD) is derived from the relationship between the power spectrum and the autocorrelation function for time-variant and non-stationary processes [47]. The correlation function is a product of the function at a posterior time with the function at a subsequent time [48]. The original WD for real signal x ( t ) is defined as [49]
W x ( t , f ) = + x ( t τ 2 ) x ( t + τ 2 ) e j 2 π f τ d τ
The symbol ( ) in Equation (24) denotes the complex conjugate.
When the WD is applied for the analytical signal x A ( t ) , it is called the Wigner-Ville distribution (WVD). The analytical signal is defined as
x A ( t ) = x ( t ) + j x h ( t )
where x h ( t ) is the Hilbert transform of x ( t ) :
H [ x ( t ) ] = x h ( t ) = 1 π + x ( τ ) 1 t τ d τ
The WVD has been used for gear fault detection [47,48] and it is recently used for rolling element bearing to represent the time-frequency features of vibration signals [50].

3.4. Category 4: Phase-Space Dissimilarity Measurement

Phase-space dissimilarity measurements, such as fractal dimension, correlation dimension approximate entropy and largest Lyapunov exponent (LLE), are developed based on the chaos theory. These methods have been applied in biomedical engineering studies [19,51] and vibration signals [52,53]. Present studies have conducted investigations into these methods in the slew bearing condition monitoring for features extraction. The motivation is that multiple defects occur more frequently than single defects in slew bearings, therefore the phase-space dissimilarity measurements can be used to analyze the irregular or chaotic vibration signals. The phase-space can be computed as follows:
X = [ y 1 y 1 + J y 1 + 2 J y 1 + ( m 1 ) J y 2 y 2 + J y 2 + 2 J y 2 + ( m 1 ) J y 3 y 3 + J y 3 + 2 J y 3 + ( m 1 ) J y M y M + J y M + 2 J y M + ( m 1 ) J ]
where J is reconstruction delay that can be calculated by time lag/ Δ t , m is the embedding dimension and M is the number of reconstructed vectors. The relation between N , J , m and M can be defined as: N = M + ( m 1 ) J or M = N ( m 1 ) J . Thus, the dimension of phase-space is a M × m matrix.

3.4.1. Fractal Dimension

Fractal dimension is a method to quantify the signal complexity by characterizing fractal patterns. The complexity can be quantified by examining the dynamical behaviours of phase-space matrix. The frequently used algorithm of fractal dimension was proposed by Higuchi [54]. The field is rapidly growing as estimated fractal dimensions for statistically self-similar phenomena. This method has been applied in many fields including image analysis [55], neuroscience [56,57], medicine [58], physics [59] and acoustics [60]. The application of fractal dimension in vibration signal of rolling element bearing has been presented in [10,61]. The algorithm of fractal dimension is shown in (28) [10]. Once the phase-space matrix, X is obtained (27), the mean absolute length between J t h and ( J t h 1 ) element of X is defined as follows:
f o r    a = 1 : J         f o r    b = 2 : m                 L 1 ( a , b 1 ) = | X ( a , b ) X ( a , ( b 1 ) ) |                 L m ( a , b 1 ) = m e a n ( L 1 )         e n d e n d
where J and m is reconstruction delay and embedding dimension, respectively. The fractal dimension D can be calculated using standard least-square fitting method of the logarithmic natural of ( 1 / m ) and L m .

3.4.2. Correlation Dimension

Correlation dimension can be used to quantify the self-similarity of vibration signal. A large value of correlation dimension corresponds to less-similarity of vibration signal. An example application of correlation dimension in vibration signal for bearing fault diagnosis was proposed by Logan and Mathew [52]. The correlation dimension was applied for vibration signals acquired from four different conditions: (1) normal (new bearing); (2) outer race fault; (3) inner race fault; and (4) roller fault. The damages were artificially introduced to the bearing parts. The method is able to identify four different bearing conditions.
The commonly used procedure to calculate the correlation dimension was introduced by Grassberger and Procaccia [62]. Similar to fractal dimension, the phase-space matrix (27) is fed to the correlation dimension algorithm as an input matrix. The formula for correlation dimensionis given by:
C ( l ) = lim x [ 2 M 2 i = 1 M k j = i + k M Θ ( l | X i X j | ) ]
where X i and X j are the position vectors on attractor, l is the distance under consideration, Θ ( x ) is the Heaviside step function, Θ ( x ) = 0 if X 0 , Θ ( x ) = 1 if X > 0 , k is the summation offset, M is the number of reconstructed vectors from the original vibration signal and C ( l ) is the correlation dimension.

3.4.3. Approximate Entropy

The approximate entropy quantifies the degree of regularity in the vibration signal. The smaller the value of the approximate entropy, the more the regularity of the behavior of the signal is. An application example of approximate entropy in vibration signal was proposed by Yan et al. [63] for high speed rolling element bearing. Moreover, Yan et al. [63] mentioned that the deterioration of bearing condition corresponds to the increase of number of frequency components. In addition, when the bearing starts to deteriorate the regularity will decrease and the corresponding approximate entropy value will increase. Similar to fractal dimension and correlation dimension, approximate entropy also uses phase-space matrix, X in (27) as an input. The initial step analysis of the approximate entropy algorithm is to measure the distance between two vectors X ( i ) and X ( j ) which can be defined as the maximum difference in their respective corresponding elements [63]:
d ( X ( i ) , X ( j ) ) = max k = 1 , 2 , , m ( | x ( i + k 1 ) x ( j + k 1 ) | )
where i = 1 , 2 , , N m + 1 , j = 1 , 2 , , N m + 1 and N is the number of data points in the vibration signal. A calculation that describes the similarity between the vector X ( i ) and all other vectors X ( j ) , where j i can be constructed as
C i m ( r ) = 1 N ( m 1 ) j 1 Θ { r d [ X ( i ) , X ( j ) ] }
the Heaviside step function, Θ ( x ) is similar to the symbol in correlation dimension where Θ ( x ) = 0 if X 0 and Θ ( x ) = 1 if X > 0 .
The symbol r in (31) denotes a predetermined tolerance value, defined as
r = k . s t d ( Y )
where s t d ( Y ) refers to the standard deviation of vibration signal, Y and k is a constant ( k > 0 ) . By defining
φ m ( r ) = 1 N m + 1 i ln [ C i m ( r ) ] ,       i = 1 , 2 , , N m + 1
the approximate entropy value of time series can be computed as
ApEn    ( m , r ) = lim N [ φ m ( r ) φ m + 1 ( r ) ]

3.4.4. Largest Lyapunov Exponent

The largest Lyapunov exponent (LLE) algorithm is an established method and has been used in some areas such as biomedical signal processing analysis especially for electroencephalography (EEG) signal [19]. The LLE algorithm calculates the exponential divergence (either positive or negative exponential) of two initial neighboring trajectories in a phase-space matrix. In other words, the LLE algorithm quantifies the degree of chaos of vibration signal in certain time. The degree of chaos corresponds to any local instability vibration signal due to the dynamic contact between rolling elements and defect spots during the data acquisition. The detailed discussion and the application of LLE in slew bearing vibration data is presented in [64].
The features extraction result for correlation dimension, fractal dimension and approximate entropy from slew bearing vibration data is presented in Figure 18. It can be seen that the approximate entropy shows better result in presenting the degradation condition of the slew bearing than the correlation dimension and fractal dimension. A fluctuation in the approximate entropy in the last measurement (after the 90th day) indicates that the slew bearing condition has changed.

3.5. Category 5: Complexity Measurement

3.5.1. Kolmogorov-Smirnov Test

The Kolmogorov-Smirnov (KS) test is a nonparametric test that can be used to compare or to measure the similarity of two cumulative distribution functions (CDFs), T ( x ) and R ( x ) [65]. T ( x ) is a target CDF and R ( x ) is a reference CDF. Because the signature of bearing condition can be identified by its CDF, therefore the KS test can be applied in bearing condition monitoring to distinguish bearing health signal from faulty signal. In order to distinguish the two CDFs, T ( x ) and R ( x ) , a statistical distance D is calculated. The distance D can be defined as the maximum absolute distance between T ( x ) and R ( x ) . This is defined as follows
D = max < x < | T ( x ) R ( x ) |
The KS feature of slew bearing vibration data is presented in Figure 19. The ‘0’ represents normal condition and ‘1’ represent abnormal condition. It can be seen that at the beginning of the measurement, the condition of bearing was still normal. Then it changed to abnormal condition when bearing had run for a couple of days. This result needs further investigation as the abnormal detection happened too early.

3.5.2. Sample Entropy

The sample entropy has been proposed by Richman and Moorman [66] to improve the approximate entropy. Chen et al. [67] presents an application of sample entropy for features extraction of gearbox vibration signals. However, the application of the sample entropy for bearing vibration signal is still limited [68]. The sample entropy feature for slew bearing vibration data is presented in Figure 20. It can be observed from Figure 20 that the sample entropy feature is not effective for slew bearing damage detection.

3.6. Category 6: Other Features

3.6.1. Singular Value Decomposition (SVD)

The singular value decomposition (SVD) is usually applied to quantify the periodicity of the time series data [69]. A comparison study [70] has been shown that SVD performs the Fourier decomposition in terms of information content and robustness. Formula of SVD can be found in [71].

3.6.2. Piecewise Aggregate Approximation (PAA) and Adaptive Piecewise Constant Approximation (APCA)

In online bearing condition monitoring cases, the changes of bearing condition can be effectively identified by measuring the similarity between two vibration signals (i.e., normal or reference and monitored signal) using Euclidean distance. However, the size of the measured data is a problem in Euclidean distance method. The similarity search methods are techniques that perform dimensionality data reduction methods such as piecewise aggregate approximation (PAA) [72,73] and adaptive piecewise constant approximation (APCA) [72]. The APCA [74], is another approach of PAA that allows arbitrary length segments. Either PAA or APCA have particular advantages; for example, PAA has twice as many approximating segments, and APCA is able to separate a single segment in an area of low activity and many segments in areas of high activity [74]. The PAA has been applied in previous study [39] for data pre-processing stage of circular domain features extraction. The detailed discussion and the application of PAA in slew bearing vibration data is presented in [39].

4. Conclusions

Features extraction methods in wide range application such as mechanical engineering (i.e., for bearing vibration signals), biomedical engineering area (i.e., for EEG and ECG signals) and finance area (i.e., for time series signals) have been presented. According to the certain characteristic of existing feature extraction methods, the methods in this paper were classified into six categories: (1) time-domain features extraction; (2) frequency-domain features extraction; (3) time frequency representation; (4) phase-space dissimilarity measure; (5) complexity measure; and (6) other features. Some of features extraction methods in each category have been empirically tested using actual normal to failure slew bearing vibration signal. It was shown that some features do not show bearing degradation trend clearly and the others are shown to be potential slew bearing condition parameters. The potential features are: impulse factor, margin factor, mathematical morphology operators, wavelet decomposition-kurtosis, approximate entropy and largest Lyapunov exponent.
Since slew bearing operating at very low speed with highly applied load, the vibration signals exhibit non-linear, non-stationary and chaotic behaviour. The methods such as impulse factor and margin factor measured the statistical properties of vibration signal and therefore could revealed the information related to the impulse level. Other potential methods based on the phase-space dissimilarity measurement such as approximate entropy and LLE algorithm can be used to extract pertinent bearing signal corresponding to the bearing damage progression. This is because the methods analyzed the vibration signal in certain time-frame to extract the useful information from the non-linear and non-stationary behaviour.

Acknowledgments

Wahyu Caesarendra thanks the University of Wollongong, Australia for the financial support through International Postgraduate Research Scholarship (IPRS) during this study.

Author Contributions

Wahyu Caesarendra conducted Conceptualization, resources, experimental setup, data collection, data processing, data analysis and original draft preparation. Tegoeh Tjahjowidodo provided Funding acquisition, review and editing.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Formula for Calculating Bearing Fault Frequencies [75]

  • Fault frequency of outer ring:
    F O R = | | I R r p m O R r p m | 2 [ 1 ( cos ( α ) ) d r d m ] | z
  • Fault frequency of inner ring:
    F I R = | | I R r p m O R r p m | 2 [ 1 + ( cos ( α ) ) d r d m ] | z
  • Fault frequency of rolling element:
    F R = | | I R r p m O R r p m | 2 [ d m d r ( cos ( α ) ) 2 d r d m ] |
where IRrpm and ORrpm are the rotational speeds of the inner ring and outer ring. For 1 r/min the value of IRrpm is 1 and the value of ORrpm is 0. dm denotes the mean bearing diameter, dr is diameter of the rolling element and z is number of rolling elements.

References

  1. Caesarendra, W. Vibration and Acoustic Emission-Based Condition Monitoring and Prognostic Methods for Very Low Speed Slew Bearing. Ph.D. Thesis, University of Wollongong, Wollongong, Australia, 2015. [Google Scholar]
  2. Žvokelj, M.; Zupan, S.; Prebil, I. Multivariate and multiscale monitoring of large-size low-speed bearings using ensemble empirical mode decomposition method combined with principal component analysis. Mech. Syst. Signal Process. 2010, 24, 1049–1067. [Google Scholar] [CrossRef]
  3. Žvokelj, M.; Zupan, S.; Prebil, I. Non-linear multivariate and multiscale monitoring and signal denoising strategy using kernel principal component analysis combined with ensemble empirical mode decomposition method. Mech. Syst. Signal Process. 2011, 24, 2631–2653. [Google Scholar] [CrossRef]
  4. Žvokelj, M.; Zupan, S.; Prebil, I. EEMD-based multiscale ICA method for slewing bearing fault detection and diagnosis. J. Sound Vib. 2016, 370, 394–423. [Google Scholar] [CrossRef]
  5. Feng, Y.; Huang, X.; Chen, J.; Wang, H.; Hong, R. Reliability-based residual life prediction of large-size low-speed slewing bearings. Mech. Mach. Theory 2014, 81, 94–106. [Google Scholar] [CrossRef]
  6. Lu, C.; Chen, J.; Hong, R.; Feng, Y.; Li, Y. Degradation trend estimation of slewing bearing based on LSSVM model. Mech. Syst. Signal Process. 2016, 76–77, 353–366. [Google Scholar] [CrossRef]
  7. Hua, W.; Yan, T.; Rongjing, H. Multiple physical signals based residual life prediction model of slewing bearing. J. Vibroeng. 2016, 18, 4340–4353. [Google Scholar]
  8. Henao, H.; Capolino, G.A.; Fernandez-Cabanas, M.; Filippetti, F.; Bruzzese, C.; Strangas, E.; Pusca, R.; Estima, J.; Riera-Guasp, M.; Hedayati-Kia, S. Trends in fault diagnosis for electrical machines: A review of diagnostic techniques. IEEE Ind. Electron. Mag. 2014, 8, 31–42. [Google Scholar] [CrossRef]
  9. Eftekharnejad, B.; Carrasco, M.R.; Charnley, B.; Mba, D. The application of spectral kurtosis on Acoustic Emission and vibrations from a defective bearing. Mech. Syst. Signal Process. 2011, 25, 266–284. [Google Scholar] [CrossRef] [Green Version]
  10. Caesarendra, W.; Kosasih, B.; Tieu, K.; Moodie, C.A.S. An Application of Nonlinear Feature Extraction—A Case Study for Low Speed Slewing Bearing Condition Monitoring and Prognosis. In Proceedings of the IEEE/ASME International Conference on Advanced Intelligent Mechatronics (AIM), Wollongong, Australia, 9–12 July 2013. [Google Scholar]
  11. Shen, K.; He, Z.; Chen, X.; Sun, C.; Liu, Z. A monotonic degradation assessment index of rolling bearings using fuzzy support vector data description and running time. Sensors 2012, 12, 10109–10135. [Google Scholar] [CrossRef] [PubMed]
  12. Yiakopoulos, C.T.; Gryllias, K.C.; Antoniadis, I.A. Rolling element bearing fault detection in industrial environments based on a K-means clustering approach. Expert Syst. Appl. 2011, 38, 2888–2911. [Google Scholar] [CrossRef]
  13. Yang, B.S.; Widodo, A. Introduction of Intelligent Machine Fault Diagnosis and Prognosis; Nova Science Publishers: New York, NY, USA, 2009. [Google Scholar]
  14. Widodo, A.; Yang, B.S. Wavelet support vector machine for induction machine fault diagnosis based on transient current signal. Expert Syst. Appl. 2008, 35, 307–316. [Google Scholar] [CrossRef]
  15. Widodo, A.; Kim, E.Y.; Son, J.D.; Yang, B.S.; Tan, A.C.C.; Gu, D.S.; Choi, B.K.; Mathew, J. Fault diagnosis of low speed bearing based on relevance vector machine and support vector machine. Expert Syst. Appl. 2009, 36, 7252–7261. [Google Scholar] [CrossRef]
  16. Yu, J. Bearing performance degradation assessment using locality preserving projections and Gaussian mixture models. Mech. Syst. Signal Process. 2011, 25, 2573–2588. [Google Scholar] [CrossRef]
  17. Xia, Z.; Xia, S.; Wan, L.; Cai, S. Spectral regression based fault feature extraction for bearing accelerometer sensor signals. Sensors 2012, 12, 13694–13719. [Google Scholar] [CrossRef] [PubMed]
  18. Rangayyan, R.M. Biomedical Signal Analysis: A Case-Study Approach; Wiley: New York, NY, USA, 2002. [Google Scholar]
  19. Päivinen, N.; Lammi, S.; Pitkanen, A.; Nissinen, J.; Penttonen, M.; Grӧnfors, T. Epileptic seizure detection: A nonlinear viewpoint. Comput. Methods Programs Biomed. 2005, 79, 151–159. [Google Scholar] [CrossRef] [PubMed]
  20. Serra, J. Image Analysis and Mathematical Morphology; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  21. Nishida, S.; Nakamura, M.; Miyazaki, M.; Suwazono, S.; Honda, M.; Nagamine, T.; Shibazaki, H. Construction of a morphological filter for detecting an event related potential P300 in single sweep EEG record in children. Med. Eng. Phys. 1995, 17, 425–430. [Google Scholar] [CrossRef]
  22. Nishida, S.; Nakamura, M.; Shindo, K.; Kanda, M.; Shibazaki, H. A morphological filter for extracting waveform characteristics of single sweep evoked potentials. Automatica 1997, 35, 937–943. [Google Scholar] [CrossRef]
  23. Nishida, S.; Nakamura, M.; Ikeda, A.; Shibazaki, H. Signal separation of background EEG and spike by using morphological filter. Med. Eng. Phys. 1999, 21, 601–608. [Google Scholar] [CrossRef]
  24. Sedaaghi, M.H. ECG wave detection using morphological filters. Appl. Signal Process. 1998, 5, 182–194. [Google Scholar] [CrossRef]
  25. Nikolaou, N.G.; Antoniadis, I.A. Application of morphological operators as envelope extractors for impulsive-type periodic signals. Mech. Syst. Signal Process. 2003, 17, 1147–1162. [Google Scholar] [CrossRef]
  26. Zhang, L.; Xu, J.; Yang, J.; Yang, D.; Wang, D. Multiscale morphology analysis and its application to fault diagnosis. Mech. Syst. Signal Process. 2008, 22, 597–610. [Google Scholar] [CrossRef]
  27. Wang, J.; Xu, G.; Zhang, Q.; Liang, L. Application of improved morphological filter to the extraction of impulsive attenuation signals. Mech. Syst. Signal Process. 2009, 23, 236–245. [Google Scholar] [CrossRef]
  28. Dong, Y.; Liao, M.; Zhang, X.; Wang, F. Faults diagnosis of rolling element bearings based on modified morphological method. Mech. Syst. Signal Process. 2011, 25, 1276–1286. [Google Scholar] [CrossRef]
  29. Santhana, R.; Murali, N. Early classification of bearing faults using morphological operators and fuzzy inference. IEEE Trans. Ind. Electron. 2013, 60, 567–574. [Google Scholar]
  30. Drummond, C.F.; Sutanto, D. Classification of power quality disturbances using the Iterative Hilbert Huang Transform. In Proceedings of the 14th International Conference on Harmonics and Quality of Power, Bergamo, Italy, 26–29 September 2010; p. 107. [Google Scholar]
  31. Lerch, A. An Introduction to Audio Content Analysis: Applications in Signal Processing and Music Informatics; Wiley-IEEE Press: Hoboken, NJ, USA, 2012. [Google Scholar]
  32. Antoni, J.; Randall, R.B. The spectral kurtosis: Application to the vibratory surveillance and diagnostics of rotating machines. Mech. Syst. Signal Process. 2006, 20, 308–331. [Google Scholar] [CrossRef]
  33. Feng, Z.; Liang, M.; Chu, F. Recent advances in time-frequency analysis methods for machinery fault diagnosis: A review with application examples. Mech. Syst. Signal Process. 2013, 38, 165–205. [Google Scholar] [CrossRef]
  34. Al-Badour, F.; Sunar, M.; Cheded, L. Vibration analysis of rotating machinery using time-frequency analysis and wavelet techniques. Mech. Syst. Signal Process. 2011, 25, 2083–2101. [Google Scholar] [CrossRef]
  35. Allen, J.B. Short term spectral analysis, synthesis, and modification by discrete Fourier transform. IEEE Trans. Acoust. Speech Signal Process. 1977, 25, 235–238. [Google Scholar] [CrossRef]
  36. Wang, S.; Huang, W.; 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. [Google Scholar] [CrossRef]
  37. Peng, Z.K.; Chu, F.L. Application of the wavelet transform in machine condition monitoring and fault diagnostics: A review with bibliography. Mech. Syst. Signal Process. 2004, 18, 199–221. [Google Scholar] [CrossRef]
  38. Daubechies, C.I. Ten Lectures on Wavelet; SIAM: Philadelphia, PA, USA, 1992. [Google Scholar]
  39. Caesarendra, W.; Kosasih, B.; Tieu, A.K.; Moodie, C.A.S. Circular domain features based condition monitoring for low speed slewing bearing. Mech. Syst. Signal Process. 2014, 45, 114–138. [Google Scholar] [CrossRef]
  40. Niu, G.; Widodo, A.; Son, J.D.; Yang, B.S.; Hwang, D.H.; Kang, D.S. Decision-level fusion based on wavelet decomposition for induction motor fault diagnosis using transient current signal. Expert Syst. Appl. 2008, 35, 918–928. [Google Scholar] [CrossRef]
  41. Liu, B.; Ling, S.F.; Meng, Q. Machinery diagnosis based on wavelets packets. J. Vib. Control 1997, 3, 5–17. [Google Scholar]
  42. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. A 1998, 454, 903–995. [Google Scholar] [CrossRef]
  43. Braun, S.; Feldman, M. Decomposition of non-stationary signals into varying time scales: Some aspects of the EMD and HVD methods. Mech. Syst. Signal Process. 2011, 25, 2608–2630. [Google Scholar] [CrossRef]
  44. Lei, Y.; Lin, J.; He, Z.; Zuo, M.J. A review on empirical mode decomposition in fault diagnosis of rotating machinery. Mech. Syst. Signal Process. 2013, 35, 108–126. [Google Scholar] [CrossRef]
  45. Caesarendra, W.; Kosasih, P.B.; Tieu, A.K.; Moodie, C.A.S.; Choi, B.K. Condition monitoring of naturally damaged slow speed slewing bearing based on ensemble empirical mode decomposition. J. Mech. Sci. Technol. 2013, 27, 1–10. [Google Scholar] [CrossRef]
  46. Rilling, G.; Flandrin, P.; Gonçalves, P. On empirical mode decomposition and its algorithms. In Proceedings of the IEEE-EURASiP Workshop on Nonlinear Signal and Image Processing, Grado, Italy, 8–11 June 2003. [Google Scholar]
  47. Staszewski, W.J.; Worden, K.; Tomlinson, G.R. Time-frequency analysis in gearbox fault detection using the Wigner-Ville distribution and pattern recognition. Mech. Syst. Signal Process. 1997, 11, 673–692. [Google Scholar] [CrossRef]
  48. Loutridis, S.J. Instantaneous energy density as a feature for gear fault detection. Mech. Syst. Signal Process. 2006, 20, 1239–1253. [Google Scholar] [CrossRef]
  49. Claasen, T.A.M.; Mecklenbrauker, W.F.G. The Wigner distribution—A tool for time-frequency analysis. Part 1: Continuous time signals. Philips J. Res. 1980, 35, 217–250. [Google Scholar]
  50. Dong, G.; Chen, J. Noise resistant time frequency analysis and application in fault diagnosis of rolling element bearings. Mech. Syst. Signal Process. 2012, 33, 212–236. [Google Scholar] [CrossRef]
  51. Ubeyli, E.D. Adaptive neuro-fuzzy inference system for classification of EEG signals using Lyapunov exponents. Comput. Methods Programs Biomed. 2009, 93, 313–321. [Google Scholar] [CrossRef] [PubMed]
  52. Logan, D.; Mathew, J. Using the correlation dimension for vibration fault diagnosis of rolling element bearings-I. Basics concepts. Mech. Syst. Signal Process. 1996, 10, 241–250. [Google Scholar] [CrossRef]
  53. Lu, C.; Sun, Q.; Tao, L.; Liu, H.; Lu, C. Bearing health assessment based on chaotic characteristics. Shock Vib. 2013, 20, 519–530. [Google Scholar] [CrossRef]
  54. Higuchi, T. Relationship between the fractal dimension and the power law index for a time series: A numerical investigation. Phys. D Nonlinear Phenom. 1990, 46, 254–264. [Google Scholar] [CrossRef]
  55. N’Diaye, M.; Degeratu, C.; Bouler, J.M.; Chappard, D. Biomaterial porosity determined by fractal dimensions, succolarity and lacunarity on microcomputed tomographic images. Mater. Sci. Eng. C 2013, 33, 2025–2030. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Accardo, A.; Affinito, M.; Carrozi, M.; Bouquet, F. Use of the fractal dimension for the analysis of electroencephalographic time series. Biol. Cybern. 1997, 77, 339–350. [Google Scholar] [PubMed]
  57. King, R.D.; George, A.T.; Jeon, T.; Hynan, L.S.; Youn, T.S.; Kennedy, D.N.; Dickerson, B.; Alzheimer’s Disease Neuroimaging Initiative. Characterization of atrophic changes in the cerebral cortex using fractal dimensional analysis. Brain Imaging Behav. 2009, 3, 154–166. [Google Scholar] [CrossRef] [PubMed]
  58. Liu, J.Z.; Zhang, L.D.; Yue, G.H. Fractal dimension in human Cerebellum measured by magnetic resonance imaging. Biophys. J. 2003, 85, 4041–4046. [Google Scholar] [CrossRef]
  59. Roberts, A.J.; Cronin, A. Unbiased estimation of multi-fractal dimensions of finite data sets. Phys. A Stat. Mech. Appl. 1996, 233, 867–878. [Google Scholar] [CrossRef]
  60. Maragos, P.; Potamianos, A. Fractal dimensions of speech sounds: Computation and application to automatic speech recognition. J. Acoust. Soc. Am. 1999, 105, 1925–1932. [Google Scholar] [CrossRef] [PubMed]
  61. Yang, J.; Zhang, Y.; Zhu, Y. Intelligent fault diagnosis of rolling element bearing based on SVMs and fractal dimension. Mech. Syst. Signal Process. 2007, 21, 2012–2024. [Google Scholar] [CrossRef]
  62. Grassberger, P.; Procaccia, I. Estimation of the Kolmogorov entropy from a chaotic signal. Phys. Rev. A 1983, 28, 2591–2593. [Google Scholar] [CrossRef]
  63. Yan, R.; Gao, R.X. Approximate entropy as a diagnostic tool for machine health monitoring. Mech. Syst. Signal Process. 2007, 21, 824–839. [Google Scholar] [CrossRef]
  64. Caesarendra, W.; Kosasih, B.; Tieu, A.K.; Moodie, C.A.S. Application of the largest Lyapunov exponent algorithm for feature extraction in low speed slew bearing condition monitorin. Mech. Syst. Signal Process. 2005, 50, 116–138. [Google Scholar]
  65. Cong, F.; Chen, J.; Pan, Y. Kolmogorov-Smirnov test for rolling bearing performance degradation assessment and prognosis. J. Vib. Control 2011, 17, 1337–1347. [Google Scholar] [CrossRef]
  66. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [PubMed]
  67. Chen, X.; Yin, C.; He, W. Feature extraction of gearbox vibration signals based on EEMD and sample entropy. In Proceedings of the 10th International Conference on Fuzzy Systems and Knowledge Discovery (FSKD), Shenyang, China, 23–25 July 2013; pp. 811–815. [Google Scholar]
  68. Wong, M.L.D.; Liu, C.; Nandi, A.K. Classification of ball bearing faults using entropic measures. Surveillance 2013, 7, 1–8. [Google Scholar]
  69. Kanjilal, P.P.; Palit, S. On multiple pattern extraction using singular value decomposition. IEEE Trans. Signal Process. 1995, 43, 1536–1540. [Google Scholar] [CrossRef]
  70. Kanjilal, P.P.; Bhattacharya, J.; Saha, G. Robust method for periodicity detection and characterization of irregular cyclical series in term of embedded periodic components. Phys. Rev. E 1999, 59, 4013–4025. [Google Scholar] [CrossRef]
  71. Qiu, H.; Lee, J.; Lin, J.; Yu, G. Wavelet filter-based weak signature detection method and its application on rolling element bearing prognostics. J. Sound Vib. 2006, 289, 1066–1090. [Google Scholar] [CrossRef]
  72. Keogh, E.; Chakrabarti, K.; Pazzani, M.; Mehrotra, S. Dimensionality reduction for fast similarity search in large time series databases. Knowl. Inf. Syst. 2001, 3, 263–286. [Google Scholar] [CrossRef]
  73. Yi, B.K.; Faloutsos, C. Fast time sequence indexing for arbitrary Lp norms. In Proceedings of the 26th International Conference on Very Large Data Bases Cairo (VLDB ’00), Cairo, Egypt, 10–14 September 2000; pp. 385–394. [Google Scholar]
  74. Chakrabarti, K.; Keogh, E.; Mehrotra, S.; Pazzani, M. Locally adaptive dimensionality reduction for indexing large time series databases. ACM Trans. Database Syst. 2002, 27, 188–228. [Google Scholar] [CrossRef]
  75. Eschmann, P.; Hasbargen, L.; Weigand, K. Die Wälzlagerpraxis: Handbuch für die Berechnung und Gestaltung von Lagerungen; R. Oldenburg: München, Germany, 1953. [Google Scholar]
Figure 1. Photo of slew bearing lab test rig and sensor location.
Figure 1. Photo of slew bearing lab test rig and sensor location.
Machines 05 00021 g001
Figure 2. Schematic diagram of slew bearing test-rig and data acquisition flow process.
Figure 2. Schematic diagram of slew bearing test-rig and data acquisition flow process.
Machines 05 00021 g002
Figure 3. Statistical time-domain features calculated from the vibration slew bearing signal.
Figure 3. Statistical time-domain features calculated from the vibration slew bearing signal.
Machines 05 00021 g003
Figure 4. AR features (up to 8 coefficients) calculated from the vibration slew bearing signal.
Figure 4. AR features (up to 8 coefficients) calculated from the vibration slew bearing signal.
Machines 05 00021 g004
Figure 5. Impulse factor (IF) and margin factor (MF) feature calculated from the vibration slew bearing signal.
Figure 5. Impulse factor (IF) and margin factor (MF) feature calculated from the vibration slew bearing signal.
Machines 05 00021 g005
Figure 6. Hjorts’ features calculated from the vibration slew bearing signal.
Figure 6. Hjorts’ features calculated from the vibration slew bearing signal.
Machines 05 00021 g006
Figure 7. Illustration of four MM operators (second to fifth row).
Figure 7. Illustration of four MM operators (second to fifth row).
Machines 05 00021 g007
Figure 8. RMS feature calculated from MM operators for 0.9 T (T = 1/BPFI): (a) erosion; (b) dilation; (c) closing and (d) opening.
Figure 8. RMS feature calculated from MM operators for 0.9 T (T = 1/BPFI): (a) erosion; (b) dilation; (c) closing and (d) opening.
Machines 05 00021 g008
Figure 9. RMS feature calculated from MM operators for 0.9 T (T = 1/BPFO): (a) erosion; (b) dilation; (c) closing and (d) opening.
Figure 9. RMS feature calculated from MM operators for 0.9 T (T = 1/BPFO): (a) erosion; (b) dilation; (c) closing and (d) opening.
Machines 05 00021 g009
Figure 10. RMS feature calculated from MM operators for 0.9 T (T = 1/BSF): (a) erosion; (b) dilation; (c) closing and (d) opening.
Figure 10. RMS feature calculated from MM operators for 0.9 T (T = 1/BSF): (a) erosion; (b) dilation; (c) closing and (d) opening.
Machines 05 00021 g010
Figure 11. Frequency-domain features: FC, RMSF and RVF calculated from the vibration slew bearing signal.
Figure 11. Frequency-domain features: FC, RMSF and RVF calculated from the vibration slew bearing signal.
Machines 05 00021 g011
Figure 12. Schematic visualization of hop size and overlapping block length in STFT algorithm.
Figure 12. Schematic visualization of hop size and overlapping block length in STFT algorithm.
Machines 05 00021 g012
Figure 13. STFT of the vibration slew bearing signal (data on 3 May).
Figure 13. STFT of the vibration slew bearing signal (data on 3 May).
Machines 05 00021 g013
Figure 14. Three level wavelet decompositions of the vibration slew bearing signal.
Figure 14. Three level wavelet decompositions of the vibration slew bearing signal.
Machines 05 00021 g014
Figure 15. Statistical-based features: mean, variance, skewness and kurtosis calculated from D3 signal of wavelet decomposition.
Figure 15. Statistical-based features: mean, variance, skewness and kurtosis calculated from D3 signal of wavelet decomposition.
Machines 05 00021 g015
Figure 16. (a) A photo of damaged axial rollers; (b) A photo of outer raceway damage.
Figure 16. (a) A photo of damaged axial rollers; (b) A photo of outer raceway damage.
Machines 05 00021 g016
Figure 17. Statistical-based features calculated from EMD decomposition.
Figure 17. Statistical-based features calculated from EMD decomposition.
Machines 05 00021 g017
Figure 18. Fractal dimension, correlation dimension and approximate entropy features calculated from the vibration slew bearing signal.
Figure 18. Fractal dimension, correlation dimension and approximate entropy features calculated from the vibration slew bearing signal.
Machines 05 00021 g018
Figure 19. KS test feature calculated extracted from the vibration slew bearing signal.
Figure 19. KS test feature calculated extracted from the vibration slew bearing signal.
Machines 05 00021 g019
Figure 20. Sample entropy feature extracted from the vibration slew bearing signal.
Figure 20. Sample entropy feature extracted from the vibration slew bearing signal.
Machines 05 00021 g020
Table 1. Brief review of statistical time-domain features extraction applied in condition monitoring of typical rolling element bearings [10].
Table 1. Brief review of statistical time-domain features extraction applied in condition monitoring of typical rolling element bearings [10].
Feature NameDescription
Brief DefinitionFormula
RMSThe RMS value increase gradually as fault developed. However, RMS is unable to provide the information of incipient fault stage while it increases with the fault development [11]. RMS = 1 N i = 1 N x i 2
VarianceVariance measures the dispersion of a signal around their reference mean value. Var = i = 1 N ( x i m ) 2 ( N 1 ) σ 2
SkewnessSkewness quantifies the asymmetry behavior of vibration signal through its probability density function (PDF). Sk = i = 1 N ( x i m ) 3 ( N 1 ) σ 3
KurtosisKurtosis quantifies the peak value of the PDF. The kurtosis value for normal rolling element bearing is well-recognized as 3. Ku = i = 1 N ( x i m ) 4 ( N 1 ) σ 4
Shape factorShape factor is a value that is affected by an object’s shape but is independent of its dimensions [12]. SF = 1 N i = 1 N x i 2 1 N i = 1 N | x i |
Crest factorCrest factor (CF) calculates how much impact occur during the rolling element and raceway contact. CF is appropriate for “spiky signals” [12]. CF = max | x i | 1 N i = 1 N x i 2
EntropyEntropy, e ( p ) , is a calculation of the uncertainty and randomness of a sampled vibration data. Given a set of probabilities, ( p 1 , p 2 , , p n ) , the entropy can be calculated using the formulas as shown in the right column. e ( p ) = i = 1 n p ( z i ) log 2 p ( z i )
Table 2. Fault frequencies of slew bearing test-rig (run at 1 r/min).
Table 2. Fault frequencies of slew bearing test-rig (run at 1 r/min).
Defect ModeFault Frequencies (Hz) (Calculation is Given in Appendix A)
AxialRadial
Outer raceway (BPFO)1.320.55
Inner raceway (BPFI)1.370.55
Rolling element (BSF)0.430.54
Table 3. Decomposition result of EMD method for the slew bearing data.
Table 3. Decomposition result of EMD method for the slew bearing data.
Vibration Data IMFs of EMD (Hz)
IMF2IMF3IMF4IMF11IMF12IMF13IMF14IMF15
24 February641.81694.52390.352.471.440.690.331.99
3 May702.88684.34346.232.041.230.640.430.19
30 August651.72679.86245.921.370.680.330.10

Share and Cite

MDPI and ACS Style

Caesarendra, W.; Tjahjowidodo, T. A Review of Feature Extraction Methods in Vibration-Based Condition Monitoring and Its Application for Degradation Trend Estimation of Low-Speed Slew Bearing. Machines 2017, 5, 21. https://doi.org/10.3390/machines5040021

AMA Style

Caesarendra W, Tjahjowidodo T. A Review of Feature Extraction Methods in Vibration-Based Condition Monitoring and Its Application for Degradation Trend Estimation of Low-Speed Slew Bearing. Machines. 2017; 5(4):21. https://doi.org/10.3390/machines5040021

Chicago/Turabian Style

Caesarendra, Wahyu, and Tegoeh Tjahjowidodo. 2017. "A Review of Feature Extraction Methods in Vibration-Based Condition Monitoring and Its Application for Degradation Trend Estimation of Low-Speed Slew Bearing" Machines 5, no. 4: 21. https://doi.org/10.3390/machines5040021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop