Open AccessArticle

An Effective Bearing Fault Diagnosis Technique via Local Robust Principal Component Analysis and Multi-Scale Permutation Entropy

Mao Ge

^1,2

Yong Lv

^1,2,*,

Yi Zhang

^1,2

Cancan Yi

^1,2 and

Yubo Ma

^1,2

Key Laboratory of Metallurgical Equipment and Control Technology, Wuhan University of Science and Technology, Ministry of Education, Wuhan 430081, China

Hubei Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, Wuhan 430081, China

Author to whom correspondence should be addressed.

Entropy 2019, 21(10), 959; https://doi.org/10.3390/e21100959

Submission received: 11 September 2019 / Revised: 27 September 2019 / Accepted: 28 September 2019 / Published: 30 September 2019

(This article belongs to the Special Issue Multiscale Entropy Approaches and Their Applications)

Download

Browse Figures

Versions Notes

Abstract

The acquired bearing fault signal usually reveals nonlinear and non-stationary nature. Moreover, in the actual environment, some other interference components and strong background noise are unavoidable, which lead to the fault feature signal being weak. Considering the above issues, an effective bearing fault diagnosis technique via local robust principal component analysis (LRPCA) and multi-scale permutation entropy (MSPE) was introduced in this paper. Robust principal component analysis (RPCA) has proven to be a powerful de-noising method, which can extract a low-dimensional submanifold structure representing signal feature from the signal trajectory matrix. However, RPCA can only handle single-component signal. Therefore, in order to suppress background noise, an improved RPCA method named LRPCA is proposed to decompose the signal into several single-components. Since MSPE can efficiently evaluate the dynamic complexity and randomness of the signals under different scales, the fault-related single-components can be identified according the MPSE characteristic of the signals. Thereafter, these identified components are combined into a one-dimensional signal to represent the fault feature component for further diagnosis. The numerical simulation experimentation and the analysis of bearing outer race fault data both verified the effectiveness of the proposed technique.

Keywords:

bearing fault diagnosis; weak fault; multi-component signal; local robust principal component analysis; multi-scale permutation entropy

1. Introduction

The bearing as an essential element has been widely used in rotating machinery [1,2]. Due to the severe working conditions, such as long and uninterrupted operation, alternating loads, and corrosion, the probability of bearing failure increases greatly, which may cause heavy economic losses or even serious personal injury [3]. Hence, an available diagnosis technique for bearing faults is highly valuable [4]. The bearing faults can be classified into three main types: inner race fault, outer race fault, and rolling element fault [5]. The vibration signals of the bearings are often used for fault diagnosis for their containing abundant equipment operation information [6]. When the bearing faults occur, the corresponding vibration signals will produce periodic impulses, and the feature of this signal behaves in a typical nonlinear and non-stationary nature, which increases in spectral complexity [5,7,8]. In the actual industrial production environment, the signals usually contain some interference vibrations caused by other mechanical components and strong background noise besides useful fault feature. Especially in the early stages of the bearing fault, the fault feature is weak and completely drowned by the strong noise and interferences. Therefore, in order to realize accurate diagnosis of bearing faults, suppressing the background noise and extracting weak fault features from multi-component signals are becoming an urgent work to be solved.

For the diagnosis of the bearing fault signals, some researchers have proposed many methods. Ciabattoni et al. [9] proposed a novel bearing fault classification method by adopting the empirical cumulative distribution functions (ECDFs) of the signal statistical spectral images as the fault feature vectors. The wavelet transform (WT) is the inner product operation between a translated and dialed wavelet basis function and the raw time domain signal. The different feature components and noise in the signal can be separated by the obtained wavelet coefficients [10,11]. Wang et al. [12] extracted the weak fault feature of the rolling element via wavelet packet transform method. Deng et al. [13] presented a novel fault diagnosis method for a motor bearing based on integrating empirical wavelet transform (IEWT) and fuzzy entropy. Xiao et al. [14] applied the wavelet threshold denoising method to effectively de-noise a rolling bearing signal. However, the diagnostic performance of these methods depends on the selection of the wavelet basis functions and the threshold. The Wigner–Ville distribution (WVD) [15] can extract ridges representing feature information from two-dimensional time-frequency plane of the non-stationary signals. Ming et al. [16] applied the cyclic Wiener filter to detect the rolling bearing fault, which uses the spectral coherence theory induced by the second-order cyclostationary signal to extract the weak fault feature. Nevertheless, WVD will produce cross-terms when analyzing multi-component signals. As far as the adaptive signal processing techniques, the empirical mode decomposition method (EMD) and the local mean decomposition (LMD) can be used to deal with nonlinear and nonstationary signals. LMD can decompose any signal into product functions (PFs) representing different feature components [17]. Li et al. [18] introduced a fault diagnosis scheme based on local mean decomposition and an improved multi-scale fuzzy entropy to realize the automatic identification of the bearing fault patterns. LMD is inadequate in processing the signals containing narrow bandwidth components. The multi-component signals can be decomposed into a series of intrinsic mode functions (IMFs) with physical meaning by EMD [19]. Imaouchen et al. [20] employed some demodulation analysis methods based on frequency-weighted energy operator and complementary ensemble empirical mode decompositions to identify the early weak faults of the bearing. Bustos et al. [21] successfully identified the operating state of the gears in high-speed trains through the EMD-based methodology. But the EMD and its improved version always suffer from modal aliasing and boundary effect. Moreover, they are also sensitive to noise. The above research provides rich reference information for bearing fault diagnosis.

It has become quite familiar to view the dynamic characteristics of different features of the raw system by reconstructing the observed time series from nonlinear non-stationary systems into a high dimensional phase space [22,23]. In that way, extracting fault features from high a dimensional phase space is a feasible scheme. The singular value decomposition (SVD) method [24] can decompose the signal trajectory matrix into series interpretable components. The singular values obtained can effectively display the intrinsic properties of different feature components and noise in raw signal. The fault feature can be extracted by setting the singular values representing interference components and background noise to zeros. Currently, the selection of singular values representing fault feature components still depends on experience, which may lead to considerable error. Especially for the early weak faults of the bearings, the singular values representing different feature components are almost impossible to be identified [25]. The classical manifold learning theory holds that the feature component of the signal matrix has a lower intrinsic dimension, which is distributed in a low-dimensional submanifold of a high dimensional phase space [26,27]. As a widely used dimension reduction method, the RPCA can extract this submanifold structure through a rank function constraint based on low-rank matrix approximation (LRMA) and simultaneously suppress background noise through a

l_{0}

-norm regularization strategy [28]. RPCA has proven to be a powerful de-noising tool in image processing, computer vision, and so on [29,30]. However, RPCA is inoperative for the separation of submanifold structures composed of multiple feature components; that is, it cannot process a multi-component signal. Recently, a novel, convex, locally sensitive, low rank matrix approximation (CLSLRMA) method [31] was introduced into the data completion problem, which significantly relaxes the assumption in LRMA that the feature component in matrix has a low-rank submanifold structure. CLSLRMA can decompose a matrix drawn from linear mixture of multiple low-rank manifold subspaces into their respective single subspaces. Hence, it is a feasible way to decompose the trajectory matrix composed of multi-component signals by CLSLRMA.

Permutation entropy (PE) [32] can efficiently evaluate the dynamic complexity and randomness of the signal time series through measuring similarity among the ordinal patterns extracted from the series, which has been widely used for the fault diagnosis of mechanical equipment [33,34]. The dynamic complexity of the bearings will change with the occurrence of faults, resulting in the changing PE values of the vibration signals [33]. However, because of the strong nonlinear and non-stationary characteristic of the acquired mechanical fault feature signals, their complex dynamic characteristics can usually hardly be fully displayed on the original scale, while some important information may also exit over multiple spatial-temporal levels (scales) [32]. Fortunately, based on PE, the multi-scale permutation entropy (MSPE) [32,35] has proven to be one of the most effective methods for which one can explicitly explain the characteristic information from the multiple time scales present in complex time series. Therefore, the MSPE of the signal was adopted in this paper to identify the feature component signal representing the bearing faults.

In this paper, an effective bearing fault diagnosis technique via local robust principal component analysis (LRPCA) and MSPE is introduced. Firstly, on the basis of noise suppression, we proposed an improved RPCA method to decompose the signal into several single-components, which was termed LRPCA. According to CLSLRMA, in the phase space of the weighted matrix associated with different anchor point, we assume that the signal trajectory matrix behaves as a combination of a noise component and a low-dimensional submanifold component, and those submanifold components represent different feature components in the raw signal. LRPCA shows that those submanifold components can be approximated by low-rank matrices through solving a convex program about a weighted combination of the matrix rank constraint function and the

l_{0}

-norm regularization [36,37]. After that, the MPSE was adopted to identify the low-rank matrices corresponding to the fault feature component. Finally, the identified low-rank matrices were transformed into a one-dimensional signal to represent the global approximation of the fault feature component for further diagnosis via weighted Nadaraya–Watson regression model [38]. The processing of the numerical simulation data and the experimental bearing fault data both verified that the proposed technique can provide a great diagnostic performance for bearing faults.

The rest of the paper is organized as follows: Section 2 introduces the theory description, wherein Section 2.1 defines some notations and abbreviations used in this paper; Section 2.2 illustrates the proposed LPRCA method; Section 2.3 describes the MSPE; the detailed step of the proposed bearing fault diagnosis technique is presented in Section 2.4. The analysis of the simulated signal and the experimental signal are performed in Section 3. Section 4 draws the conclusions.

2. Theory Description

2.1. Notations and Abbreviations

Throughout this paper, we use lowercase letters for scalars, e.g.,

x

; bold and lowercase letters for vectors, e.g.,

x \in ℝ^{n_{1}}

; boldface and uppercase letters for matrices, e.g.,

X \in ℝ^{n_{1} \times n_{2}}

. The anchor points in

X

are marked as

e_{i} = (a_{i}, b_{i}), a_{i} = 1, \dots, n_{1}; b_{i} = 1, \dots, n_{2}; i = 1, \dots, m

; the i-th element of

x

can be expressed as

x_{i}

, the (i,j)-th element of

X

is denoted as

X (i, j)

and the i-th row of

X

is expressed as

X (i, :)

. The Hadamard product of two matrix

A \in ℝ^{n_{1} \times n_{2}}

B \in ℝ^{n_{1} \times n_{2}}

is defined as

C = A ⊙ B \in ℝ^{n_{1} \times n_{2}}

with its element

C (i, j) = A (i, j) B (i, j)

. The

l_{0}

-norm

{‖ X ‖}_{0}

represents the sum of nonzero elements of

X

and the

l_{1}

-norm

{‖ X ‖}_{1}

the sum of absolute value of all elements in

X

. The SVD of

X

is defined as

X = U \sum V

, where

U

and

V

are the left and right singular value matrices;

Σ = d i a g ({σ_{i}}_{1 \leq i \leq r})

represents the singular value matrix. The nuclear norm of

X

is denoted as

{‖ X ‖}_{*} = \sum_{i} σ_{i}

In the rest of this paper, the following abbreviations are used: RPCA—robust principal component analysis; LRPCA—local robust principal component analysis; MSPE—multi-scale permutation entropy; CLSLRMA—convex local sensitive low rank matrix approximation; SVD—singular value decomposition; SSA—singular spectrum analysis; EMD—empirical mode decomposition; ADMM—alternating direction method of multipliers; and SNR—signal to noise ratio.

2.2. Decomposing a Signal into Single-Components via LRPCA

2.2.1. RPCA

The acquired one-dimensional bearing fault signal

x \in ℝ^{n}

can be converted into a high dimensional signal trajectory matrix

X \in ℝ^{n_{1} \times n_{2}}

by phase space reconstruction, which is based on a embedding process with the parameter of the embedding dimension

n_{1}

and the delay time

τ

(where

(n_{1} - 1) τ + n_{2} = n

) [22]:

X = [\begin{matrix} x_{1} & x_{2} & \dots & x_{n_{2}} \\ x_{1 + τ} & x_{2 + τ} & \dots & x_{n_{2} + τ} \\ ⋮ & ⋮ & \dots & ⋮ \\ x_{1 + (n_{1} - 1) τ} & x_{2 + (n_{1} - 1) τ} & ⋱ & x_{n} \end{matrix}]

(1)

Except for the strong background noise component,

X

is composed of multiple feature components, which include the fault feature component and the unwanted interfering components. Hence, on the basis of noise suppression, how to separate the useful feature component from these mixed multi-components and back it to the one-dimensional signal to represent the extracted fault feature component was an inevitable task for extracting weak faults in this paper.

By referring to manifold learning theory, it can be found that the feature component in

X

has a low-dimensional submanifold structure [26,27]. RPCA can extract this structure by solving the following low-rank matrix and sparse matrix decomposition model (shown as Figure 1):

\underset{L, S}{\min \arg} r a n k (L) + α {‖ S ‖}_{0}, s . t . X = L + S

(2)

where

α = 1 / \sqrt{\max (n_{1}, n_{2})}

is the optimal weighted parameter. The first term (

L \in ℝ^{n_{1} \times n_{2}}

r (L) ≪ r (X)

) of this model is a rank constraint function, which uses a low-rank matrix to estimate the low-dimensional submanifold structure, and the second term (

S \in ℝ^{n_{1} \times n_{2}}

) is a regularization strategy, which is mainly used to correct the deviation of matrix data caused by noise. In fault feature extraction,

L

represents the signal feature component and

S

captures the noise [28]. Thus, RPCA can effectively separate the feature component of the signal from the noise. However, it can be seen from the theoretical basis of Equation (2) that RPCA can only deal with the matrix data composed of a single feature component and noise, which means this method cannot separate the data containing the multiple feature components.

2.2.2. LRPCA

To tackle the drawback of RPCA in processing the multi-component signals, we proposed a novel LRPCA method based on CLSLRMA to decompose

X

into several single-components. In LRPCA, the following fundamental assumption is introduced.

Fundamental Assumption: In addition to the background noise signal matrix

S

X

contains a feature signal matrix

L

composed of a linear mixture of several feature components. Furthermore, each feature component corresponds to a low-dimensional submanifold structure hidden in the high-dimensional phase space and has a characteristically of low-rank.

Figure 2 depicts the main ideal of the LRPCA. Specifically, each submanifold structure is generally hidden in the different high-dimensional phase space

T (e_{i})

associated with the local selected anchor point (

e_{i} = (a_{i}, b_{i}), a_{i} = 1, \dots, n_{1}; b_{i} = 1, \dots, n_{2}; i = 1, \dots, m

), and

T (e_{i})

is derived from the weighting of

X

T (e_{i}) = W_{e_{i}} ⊙ X

(3)

where

W_{e_{i}} \in ℝ^{n_{1} \times n_{2}}

is a local weighted coefficient matrix.

Thereafter,

T (e_{i})

is decomposed into a low-rank component

L_{e_{i}}

and a sparse component

S_{e_{i}}

. The resulting low-rank matrices

L_{e_{1}}, \dots, L_{e_{m}}

represent the low-dimensional submanifolds, which correspond to different feature components in raw signal, that of the fault feature component and the unwanted interfering components. Conclusively,

L

is expressed as the weighted combination of these local low-rank matrices:

L = R_{1} ⊙ L_{e_{1}} + \dots + R_{m} ⊙ L_{e_{m}}

(4)

where

R_{i} \in ℝ^{n_{1} \times n_{2}}

is a weighted regression matrix. Thus, the fault-related feature components can be extracted from these low-rank matrices obtained by decomposing

L

. For this task, a novel local low-rank matrix and sparse matrix decomposition model was proposed to obtain these low-rank matrices from

X

L_{e_{i}} = \underset{L_{e_{i}}, S_{e_{i}}}{\min \arg} r a n k (L_{e_{i}}) + α {‖ S_{e_{i}} ‖}_{0}, s . t . T (e_{i}) = L_{e_{i}} + S_{e_{i}}

(5)

The RPCA is actually a special case of this model when

W_{e_{i}}

is a unit matrix. Therefore, the de-noising performance of this model can be guaranteed.

• Model Construction and Algorithm Solving

The derivation of the phase space

T (e_{i})

is the key step in LRPCA. Firstly, a distance function

d ((a, b), (a^{'}, b^{'}))

is defined to describe the similarity between any two elements

X (a, b)

and

X (a^{'}, b^{'})

X

. A smaller value of

d

means the probability of the two elements being in the same phase space is higher. According to the theory of CLSLRMA, The standard incomplete SVD

X = U V^{T}

[39] is employed to divide

d

into two independent terms and the arc-cosine function is utilized to calculate these two terms:

d (a, a^{'}) = \arccos (\frac{〈 U (a, :), U (a^{'}, :) 〉}{‖ U (a, :) ‖ ‖ U (a^{'}, :) ‖}), d (b, b^{'}) = \arccos (\frac{〈 V (b, :), V (b^{'}, :) 〉}{‖ V (b, :) ‖ ‖ V (b^{'}, :) ‖})

(6)

Then, the following non-parametric smoothing Epanechnikov kernel function [40] is adopted to define the local weighted coefficient matrix

W_{e_{i}}

W ((a, b), (a^{'}, b^{'})) = W (a, a^{'}) W (b, b^{'}) = (1 - d ({(a, a^{'})}^{2}) (1 - d {(b, b^{'})}^{2})

(7)

The (i, j)-element of

W_{e_{i}}

is expressed as

W (e_{i}, (i, j))

. The Epanechnikov kernel function is a typical unilateral quadratic decrement function related to

d

. Consequently, a lager value of the weighted coefficient means that the weight for

X (i, j)

belonging to

T (e_{i})

is bigger. Besides, for different local anchor point

e_{i}

, there theoretically should exist a unique corresponding phase space hidden the single submanifold. However, as the fundamental assumption states, there should be finite number of single submanifold structures in

X

. Gratifyingly, this function has the mathematical property of changing slowly, which means as long as the similarity of two anchor points

e_{i}

and

e_{j}

is large enough, their corresponding phase spaces

T (e_{i})

and

T (e_{i})

will have quite high similarity. Moreover, both of those two phase spaces may hide the same single submanifold structure.

Through the above theoretical analysis process, the

T (e_{i})

satisfying the fundamental assumption was successful created. Finally,

T (e_{i})

is decomposed to extract the hidden low-rank matrix

L_{e_{i}}

by solving the Equation (5). Since the discrete combination nature of the

l_{0}

-norm and rank function, the solution of this model is an non-deterministic polynomial (NP)-hard problem [41]. Take this into consideration, some recent studies [36,37,42,43,44,45] pointed out that an equivalent convex optimization program of this model can be obtained from the convex hull of the two constraints; that is, the

l_{1}

-norm and the nuclear norm are employed to replace those two constraints respectively:

L_{e_{i}} = \underset{L_{e_{i}}, S_{e_{i}}}{\min \arg} {‖ L_{s_{i}} ‖}_{*} + α {‖ S_{e_{i}} ‖}_{1}, s . t . T (s_{i}) = L_{e_{i}} + S_{e_{i}}

(8)

The solution of Equation (8) is a typical convex optimization problem, whose minimizer is globally unique [43]. Algorithm 1 provides a precise and convergent solution to this equation via ADMM [44,46]. Note that steps 3 and 4 are both convex problems, they both have closed-form solutions via singular value thresholding operator [47].

Algorithm 1 solve (8) by ADMM

Input: signal trajectory matrix

X \in ℝ^{n_{1} \times n_{2}}

;

Parameter: number of anchor points:

q

; regularization parameter:

α = 1 / \sqrt{\max (n_{1}, n_{2})}

;

for all i = 1:

m

, parallel do;

1. select

e_{i} (a_{i}, b_{i})

uniformly in

X

, and calculate

W_{e_{i}}

by Equation (7);

T (e_{i}) = W_{e_{i}} ⊙ X

;

Initialize:

L_{e_{i}}^{0} = S_{e_{i}}^{0} = Y^{0} = 0

γ_{0} = e^{- 3}

γ_{\max} = e^{10}

μ = 1.1

ε = 1 e - 8

;

while not converged do;

3. fix the others and update

L_{e_{i}}^{k + 1}

by:

L_{e_{i}}^{k + 1} = \underset{L_{e_{i}}}{\arg \min} : {‖ L_{e_{i}} ‖}_{*} + \frac{γ_{k}}{2} {‖ L_{e_{i}} + S_{e_{i}}^{k} - T (e_{i}) + Y^{k} ‖}_{F}^{2}

4. fix the others and update

S_{e_{i}}^{k + 1}

by:

S_{e_{i}}^{k + 1} = \underset{S_{e_{i}}}{\arg \min} : λ {‖ S_{e_{i}} ‖}_{1} + \frac{γ_{k}}{2} {‖ L_{e_{i}}^{k + 1} + S_{e_{i}} - T (e_{i}) + Y^{k} ‖}_{F}^{2}

5. update Lagrange multiplier

Y

Y^{k + 1} = Y^{k} + T (e_{i}) - L_{e_{i}}^{k + 1} - S_{e_{i}}^{k + 1}

;

6. update

τ

τ_{k + 1} = \min (μ τ_{k}, τ_{\max})

;

7. check the convergence conditions:

{‖ L_{e_{i}}^{k + 1} - L_{e_{i}}^{k} ‖}_{\infty} \leq ε, {‖ S_{e_{i}}^{k + 1} - S_{e_{i}}^{k} ‖}_{\infty} \leq ε, {‖ L_{e_{i}}^{k + 1} + S_{e_{i}}^{k + 1} - T (e_{i}) ‖}_{\infty} \leq ε

end;

output:

L_{e_{1}}, \dots, L_{e_{m}}

Through Algorithm 1, we can obtain

m

low-rank matrices corresponding to different feature components in the raw signal. Besides, it needs to be emphasized that the performance of the final fault diagnosis is highly susceptible to the location selection and the number of the anchor points. Once the chosen anchor points are inappropriate, there may be multiple low-rank matrices corresponding to the same feature component. More seriously, the fault-related low-rank matrices may be completely missed. In view of the above problems, for the selection of the anchor points, we do not have a good solution for now. But, in the process of analyzing the experimental data, the following two principles are feasible. The first one is that

m

should be large enough to ensure that the low-rank matrices corresponding to all feature components can be extracted specifically. This can be explained from the perspective of the probability. When the number of the anchor points is more, the probability of the extracted fault-related low-rank matrices is higher. However, this inevitably requires a lot of computing time. The other one is that these anchor points should be uniformly chosen from the elements set and the distance between any two anchor points should be made large enough. This principle is to make extracted low-rank matrices corresponding to the different feature components differ as much as possible. In the numerical simulation experiment, since the raw simulated signal contained three feature components, and we found that when

m

was set as six, the final decomposition performance was quite good. Additionally, when analyzing the experimental signal, it was found that

m = 6

is also appropriate. Therefore, according to the analysis results of the experimental data, we assumed that the number of the main feature components in the acquired bearing fault signal generally would not exceed three, and set

m

to be twice that number; that is

m = 6

• Global Approximation of Fault Feature Component

These low-rank matrices can be backed to

m

one-dimensional component signals by inverse transform. In Section 2.2, it will show that the one-dimensional components related to fault feature can be identified from these signal through the MPSE characteristic of the signal. Thus, there are

o

(

o < m

) identified one-dimensional signals and their corresponding low-rank matrices

L_{e_{1}}, \dots, L_{e_{o}}

. These low-rank matrices are actually local sensitive, which can only be used to describe the fault feature information contained in the corresponding low-dimensional submanifolds. Hence, the Nadaraya-Watson regression model [38] is adopted to combine these low-rank matrices into a global approximation

{\hat{L}}_{f} \in ℝ^{n_{1} \times n_{2}}

of the fault feature component, which can be expressed as:

{\hat{L}}_{f} = \sum_{i = 1}^{o} \frac{W_{e_{i}}}{\sum_{j = 1}^{m} W_{e_{j}}} ⊙ L_{e_{i}}

(9)

Note that if

o = m

, this equation is equivalent to Equation (4). Thus, the estimator of

L

can be obtained, which is actually a de-noising process. Thereby, we returned

{\hat{L}}_{f}

to the one-dimensional time series

{\hat{l}}_{f} \in ℝ^{n}

, which is the expected fault feature signal. Hence, the task to suppress the strong background noise and extract the fault feature signal was completed.

2.3. Identification of Bearing Fault-Related Signal through MSPE

2.3.1. Basic Theory of MSPE

The MPSE can efficiently evaluate the dynamic complexity and randomness of the time series under different scales. The calculation of the MSPE depends on three parameters: scale factor

ε

, embedding dimension

d

and time-lag

δ

[32,36]. For a signal time series

x = [x_{1}, x_{2}, \dots, x_{n}]

, the main calculation steps can be divided as follows:

(1): Transform $x$ into a successive coarse-grained time series $y^{ε} \in ℝ^{n_{ε}}$ ( $n_{ε} = [n / ε]$ ) by averaging the time data points in $x$ with the given non-overlapping time slice of the increasing length, $ε$ . Then, each element of $y^{ε}$ is defined as:

$y_{i}^{ε} = \frac{1}{ε} \sum_{j = (i - 1) ε + 1}^{i ε} x_{j}, i = 1, 2, \dots, n_{ε}$

(10)
(2): For each coarse-grained time series $y^{ε}$ , the PE value needs to be calculated. Firstly, $y^{ε}$ is cut into a series of data segments through $d$ and $δ$ :

$g^{i} = [y_{i}^{ε}, y_{i + δ}^{ε}, \dots, y_{i + (d - 1) δ}^{ε}], i = 1, 2, \dots, n_{ε} - (d - 1) δ$

(11)

There will be $n_{ε} - (d - 1) δ$ data segments in total. Then, there are $d!$ different types of ordinal patterns ( $ψ_{i}, i = 1, \dots, d!$ ) in the data segments. Then, count the frequency of each pattern and denote them as $f {(ψ}_{i}), i = 1, \dots, d!$ . Thus, the relative frequency of each pattern can be written as:

$p (ψ_{i}) = f (ψ_{i}) / (n_{ε} - (d - 1) δ)$

(12)

Finally, the PE of $y^{ε}$ is expressed as:

$P (ε) = - \sum_{i = 1}^{d!} p (ψ_{i}) \log_{2} p (ψ_{i})$

(13)

For convenience, we normalize $P (ε)$ by dividing its maximum value $\log_{2} d!$ :

$0 \leq P (ε {) / \log}_{2} d! \leq 1$

(14)
(3): The PE values of different coarse-grained time series can be obtained and plotted as a function of the scale factors. The vector $P = [P (1), P (2), \dots]$ , formed by the set of PE values, is the MSPE of the original time series. Figure 3 shows the process of coarse granulation and the data segmentation of a time series.

2.3.2. Mathematical Model of Bearing Fault Feature Signal

The bearing fault feature signal can be modeled as the combination of finite pulsed excitations [8,48,49,50]:

x (t) = a (t) \sum_{i = 0}^{I} b_{i} (t) \cos (2 π f_{e} t - c_{i} (t) + θ_{i})

(15)

where

b_{i} (t)

c (t)

, and

a_{i} (t)

indicate the amplitude modulation component, the frequency modulation component, and the modulation effect caused by vibration transfer path, respectively.

f_{e}

and

θ_{i}

indicate the system resonance frequency and the initial phase, respectively.

b_{i} (t)

and

c_{i} (t)

can be expressed as:

b_{i} (t) = B_{i} e^{- ξ (t - i T_{d} - υ_{i})} u (t - i T_{d} - υ_{i})

(16)

c_{i} (t) = \sum_{l = 1}^{L} C_{i l} \sin (2 π l f_{c} t + θ_{i l})

(17)

where

ξ

f_{c}

, and

T_{d}

indicate the resonance attenuation coefficient, the fault feature frequency, and the time period of fault, respectively.

u (t)

represents a unit step function.

υ_{i}

and

θ_{i l}

represent the random slip of the i-th pulse and the initial phase, respectively.

B_{i}

and

C_{i l}

are amplitudes.

Generally, the vibration sensor is mounted at the bearing seat, which is fixed with the outer race. In the case of the rolling element fault or the inner race fault, the vibration propagation generated by the signal transfer path from the fault location to the sensor is varies with time, resulting in an amplitude modulation effect in the signal:

a (t) = A [1 + \cos (2 π f_{r} t)]

(18)

where

f_{r}

indicates the rotation frequency of the shaft where the fault bearing is located, and

C

is a constant.

In the case of the outer race fault, the vibration propagation only has a scaling effect on the signal amplitude due to the transfer path being fixed:

a (t) = A

(19)

One accepted approach to fault identification is to identify the fault-related frequency contents from the signal spectrum. For example, Figure 4 shows the waveform and the spectrum of a simulated bearing’s inner race fault signal, and the related frequency content includes: the fault feature frequency

f_{c}

, its harmonic frequencies

n f_{c}

, the rotational frequency

f_{r}

, and the modulated side band formed by their combination

n f_{c} \pm f_{r}

; Figure 5 shows the waveform and the spectrum of a bearing outer race fault signal. The frequency content includes: the fault feature frequency

f_{c}

and its harmonic frequencies

n f_{c}

. It can be observed that nonlinear and non-stationary nature of these fault signals increases the complexity of the spectrum. Especially in the fault of inner race, where the complex modulation side band appears. As a result, the dynamic characteristic of the inner race fault should be more complex than that of outer race fault.

2.3.3. Identification of Fault-Related Components through MSPE

When faults occur in a mechanical system, such as gear or bearing parts, the MSPE value of the vibration signal will change, and the value varies with different type of the faults. Thus, MSPE can be adopted to identify the fault-related components.

The common interference components in bearing signals include the shock signals generated by other parts and the harmonic signals. Therefore, without loss of generality, we discuss the MSPE characteristics of different types of components bearing fault signals, including the fault feature signals of the bearing’s inner race (shown as Figure 4) and the bearing outer race (shown as Figure 5), a harmonic signal, a shock signal, and a Gaussian white noise signal. The results of these five type signals are shown as Figure 6a. The results demonstrate that the MSPE value of a regular time series, such as harmonic signal or shock signal, is smaller, the examples being basically below 0.3; in contrast, since advent of the complex dynamic characteristic, the MPSE values of the bearing fault feature signals are larger, ranging from 0.4 to 0.7. In addition, it can be observed that the MSPE value of the inner race fault is larger than that of the outer race. In particular, the randomness of the noise signal is the strongest with the MSPE value more than 0.9, which proves that the MSPE is sensitive to noise.

According to the above analysis results, we can set a threshold range (0.5–0.85) of the MSPE value or select the larger value (but no more than 0.9) to identify the components related to bearing fault feature from the

q

one-dimensional component signals obtained by LRPCA.

It needs to be emphasized that although MSPE has a good anti-noise ratio, it may not be able to effectively identify the early weak fault feature signals under strong background noise. Figure 6b shows the MPSE of the four feature signals in Figure 6a, after Gaussian white noise with a SNR of −5 is added. Due to the existence of the strong background noise, the randomness of signals is greatly enhanced, resulting in the MSPE values of the four feature signals increasing to more than 0.9 and mixing together. As a result, the fault feature components are impossible to be identified. Therefore, it is an urgent problem to reduce the noise before identifying the weak fault feature components.

2.4. The Process of the Effective Fault Diagnosis Technique

Figure 7 depicts the flowchart of the proposed effective fault diagnosis technique based on the above theoretical description. And the main steps are summarized as follows:

(1): Using the proposed LRPCA method to decompose the trajectory matrix consisting of the acquired bearing fault signal into multiple low-rank matrices and to suppress the noise synchronously;
(2): Convert the low-rank matrices obtained into one-dimensional component signals by inverse transformations and identify the fault-related components from these signals through the MPSE characteristic of the signal;
(3): Using the weighted Nadaraya–Watson regression model and inverse transform to combine the low-rank matrices corresponding to identified components into a one-dimensional signal to represent the extracted fault feature component;
(4): Confirm the bearing fault by identifying the fault-related frequency contents from the signal spectrum.

3. Experiments

3.1. Numerical Simulation Experiment

In order to verify the diagnostic performance of the proposed technique, without loss of generality, a simulated signal composed of multiple components was generated:

x = x_{1} (t) + x_{2} (t) + x_{3} (t) + n (t)

(20)

where

x_{1} (t)

is the simulated fault feature signal of a bearing’s inner race, as shown in Figure 4, and its detailed parameters are listed in Table 1.

x_{2} (t)

and

x_{3} (t)

are the interferences of a shock signal and a harmonic signal, respectively. Figure 8a,b shows the waveforms of these two signals and both of their feature frequencies are 150 Hz.

n (t)

represents the strong background white Gaussian noise. The signal sampling frequency and the sampling point are

f_{s} =

20,000 and

n =

20,000, respectively.

Firstly, the de-noising performance of the proposed LRPCA method was tested. The strong background noise with SNR from −5 to 5 db was added to the simulated signals to imitate the early weak fault. For visually displaying the de-noising performance of LRPCA, the methods of RPCA, SSA, and EMD were employed to make a comparative analysis. The phase space reconstruction parameters used in LRPCA, RPCA, and SSA were all set as

n_{1} =

200 and

τ =

100. EMD employed the energy difference tracking method [51] to select desired IMF components. The hard threshold method [52] was adopted in SSA to select the best combination of singular values to reconstruct the signal. Figure 9 illustrates the de-noising result of the above four methods. In the vertical axis of the graphics, the higher SNR value indicates a better de-noising performance. It is clear that the proposed LRPCA method can provide the best noise suppression performance.

Then, the fault feature extraction performance of the proposed technique was investigated. In this experiment, a noise with a SNR of −5 db was added into the simulation signal. Figure 10 shows the waveform and the spectrum of the noise-signal mixture. It can be seen that the fault feature is completely submerged by noise and the other interferences, which inevitably increases the difficulty of recognizing the fault feature. The methods of wavelet shrinkage denoising [14], basis pursuit denoising [53], EMD, and SSA were selected for comparative analysis.

In the wavelet shrinkage denoising method, the decomposition layer was selected as 4 and the wavelet basis function was set as “db15.” The waveform and the spectrum of the final extracted fault feature signal are shown in Figure 11. In the resulting spectrum, the fault feature frequency (

f_{c}

) and some of its harmonic frequencies (

2 f_{c}

−

4 f_{c}

) can be found. However, due to a large number of interference components and strong background noise, the identification of these fault-related frequencies was seriously affected, and the modulation sidebands were completely submerged. hence, this method was insufficient in feature extraction of a simulated signal.

In the basis pursuit denoising method, we adopted the compressed sensing reconstruction algorithm [54] to extract the fault feature signal and Figure 12 shows the final analysis result. In the resulting spectrum, the Fault harmonic frequencies (

3 f_{c}

−

5 f_{c}

) and their modulated sidebands (

3 f_{c} \pm f_{r}

4 f_{c} \pm f_{r}

, and

5 f_{c} \pm f_{r}

) could be identified. But, there were still a lot of interference frequency components and strong background noise, leading to the fault of the signal not being directly determined.

In the EMD, fourteen IMFs can be obtained by decomposing the signal and the waveforms of their top twelve are shown in Figure 13a. By applying Fourier transform to these IMFs, in the resulting spectrums of IMF3 and IMF4, partial fault-related frequency contents can be found; i.e., fault harmonic frequencies (

2 f_{c}

−

5 f_{c}

) and the modulated sidebands (

3 f_{c} - f_{r}

4 f_{c} - f_{r}

, and

5 f_{c} + f_{r}

). The waveform and the spectrum of these two IMFs are shown in Figure 13b–e. However, there were still plenty of interference frequency components and noise, which are adverse to the identification of fault feature.

In the SSA, by setting the threshold value of the energy of singular value to reach 95% of the total, four singular subspaces representing different feature components can be obtained. Figure 14a illustrates the waveforms of the corresponding one-dimensional feature component signals of these four subspaces obtained through inverse transform. Then, through similarity analysis, it was concluded that components 2 and 4 had higher similarity with the raw fault feature signal. Thus, the signal obtained by adding them represents the extracted fault feature component. The waveform and the spectrum of this extracted signal are shown in Figure 14b,c. In the spectrum, some fault-related frequency contents could be found. But, similar to the analysis result of EMD, there were also some interference frequency components and noises, which may have had an adverse effect on the final diagnosed result.

The diagnosed result of the proposed technique is shown in Figure 15. Figure 15a shows the waveforms of the one-dimensional component signals obtained by LRPCA, which correspond to the six low-rank matrices. Additionally, the MSPEs of those components are displayed in Figure 13b. It can be observed that the MSPE values of component 3 and 6 are relatively higher and evenly distributed between 0.5 and 0.8. Thus, we could combine the low-rank matrices corresponding to those two components into a one-dimensional signal to represent the extracted fault feature signal. Figure 15c,d show the waveform and the spectrum of this extracted signal. All fault related frequencies were clearly discernable in the spectrum. Furthermore, it can be observed that the interference frequency components and the background noise were basically eliminated. According to the above analysis information, we can undoubtedly confirm that the fault occurred at the inner race.

The above simulation results indicate that EMD and SSA perform when decomposing a complicated multi-component signal into finite single-components. However, their anti-noise ability against strong noise is insufficient and there are still some interference components in the final extracted feature signals. On the contrary, the proposed technique is more effective in eliminating the strong background noise and extracting the weak fault feature component, and its diagnostic performance for the simulated signal is obviously better than the other two methods.

3.2. Experimental Signal Analysis

The vibration signal acquired on the spot is more complex than the simulation signal. In order to further verify the practicability of the proposed technique, a pitting fault signal of the bearing’s outer race sampling from a bearing-gear fault’s experimental table was analyzed. Figure 16 shows the experimental table, which consists of an AC motor, couplings, a gearbox with two pairs of meshing gears, and a magnetic powder brake. The test bearing is a single-row tapered roller bearing of the type of 32206 and its fault was handled by electrical discharge machining (EDM) method. The red arrow in Figure 16a displays the mounting position of the bearing. The experiment parameter was listed in Table 2. The vibration data of this experiment are measured by PCB acceleration sensor in the vertical direction of the bearing.

Figure 17 is the diagram of the waveform and the spectrum of the acquired signal. It can be seen that the weak fault feature was impossible to be identified. Then, the proposed technique was utilized to diagnose the signal, and the methods of wavelet shrinkage denoising, basis pursuit denoising, EMD, and SSA were chosen for a comparative analysis.

The analysis results of the wavelet shrinkage denoising method and basis pursuit denoising method are shown as Figure 18 and Figure 19, respectively. In their result spectrums, although the fault characteristic frequency (

f_{o}

) and its harmonic frequencies (

2 f_{o}

−

3 f_{o}

) can be found, the identification of these fault-related frequencies is seriously affected by a large number of interference components and strong background noise. Therefore, the fault type of bearing cannot be directly determined from these analysis results.

Figure 20a displays the waveforms of top 12 IMFs obtained through EMD. By applying the Fourier transform to them, the fault-related frequency contents (

f_{o}

and

2 f_{o}

) can be found in the spectrums of IMF 1 and IMF 8, which are shown in Figure 20b,d,e. However, that fault feature information is hard to be recognized due to the presence of many interference frequency components and strong background noise. Therefore, it is difficult to determine whether the outer race of the bearing is faulty.

Figure 21 shows the fault feature extraction result of SSA. The peaks of fault feature frequency (

f_{o}

) and its triple frequency (

3 f_{o}

) were obvious in the spectrum. But there are still many interference peaks and noise, which affects the identification of fault feature. These above analysis results indicate that neither EMD nor SSA can provide good fault diagnosis performances for the experimental fault signal.

The proposed technique was utilized to diagnose the signal. Figure 22a shows the waveforms of the one-dimensional component signal corresponding to the six low-rank matrices obtained by LRPCA. Furthermore, their MSPE values were displayed in Figure 22b. It can be observed that the MSPE values of component 1, 3, 4, and 6 are relatively higher. Meanwhile, the MSPE value of component 4 was basically above 0.9, which may be the feature signal of other component with more complex dynamics characteristic. The MSPE values of components 1, 3, and 6 ranged from 0.7 to 0.9, and their trends were similar, so we combined the low-rank matrices corresponding to these three components into a one-dimensional feature signal representing the extracted the fault feature signal. Figure 22c,d show the waveform and the spectrum of this extracted signal. In the spectrum, the rotation frequency (

f_{r}

), the fault feature frequency (

f_{o}

), and its double frequency (

2 f_{o}

) and triple frequency (

3 f_{o}

) can be easily found. Moreover, it can be seen that the energy of the noise was suppressed at a low level and there were a few interference frequency components, which had little effect on the final diagnosed result. Consequently, we could confirm that the outer race of bearing had fault. These analysis results demonstrate that the proposed technique can effectively suppress the strong background noise and extract weak fault feature component from the multi-component signal, which means the proposed technique can provide a great diagnostic performance when dealing with the experimental signal.

4. Conclusions

In general, the bearing’s weak fault feature exhibits the nature of nonlinear and non-stationary, which is hard to be extracted under the situation of existing strong background noise and interference components. In considering this problem, an effective bearing fault diagnosis technique via LRPCA and MSPE was introduced in this paper. The LRPCA can decompose the signal trajectory matrix into multiple low-rank matrices, and meanwhile, suppress the noise. The MSPE was used to identify the low-rank matrices corresponding to bearing’s fault feature. Thereafter, those identified low-rank matrices were combined into a one-dimensional signal to represent the extracted fault feature component for further fault diagnosis. The principle and the effectiveness of this technique was verified by the analysis of both the simulation signal and the acquired bearing fault signal. The analysis results indicate that the proposed technique can effectively detect and locate the bearing faults accurately.

The threshold range representing the bearing fault feature component was determined based on the results of simulation experiments. However, the feature components of the acquired bearing fault signals should be much more complex than simulated signals. Therefore, we will focus on determining a more appropriate threshold range for acquired signal in future work. In addition, the de-nosing performance of proposed LRPCA method was verified by the signal simulation Equation (20), and Figure 9 shows that it performs better than other methods. In fact, the robustness of the method can be improved by statistical evidence. By extending the evaluation to greater number of rounds, the results of each round are collected, and through meaningful statistical tests, the parameters of the method can be optimized, leading to a significant improvement in noise reduction performance. Therefore, we will also carry out more deep research to that end in future work.

Author Contributions

Conceptualization: Y.L.; Methodology: Y.L. and M.G.; Investigation: C.Y. and Y.Z.; Data curation: Y.Z.; Funding acquisition: Y.L.; Project administration: C.Y.; Resources: C.Y. and Y.Z.; Software: M.G. and Y.M.; Supervision: Y.L.; Validation: M.G. and Y.M.; Writing-original draft: M.G.; Writing-review and editing: Y.L. and Y.M.

Funding

This research was funded by the National Natural Science Foundation of China under Grants No.51875416 and No.51475339, the Wuhan Science and Technology Project under Grants No.2017010201010115.

Conflicts of Interest

The authors declare that there is no conflict of interests.

References

Zhang, W.; Li, C.; Peng, G.; Chen, Y.; Zhang, Z. A deep convolutional neural network with new training methods for bearing fault diagnosis under noisy environment and different working load. Mech. Syst. Signal Process. 2018, 100, 439–453. [Google Scholar] [CrossRef]
Wang, S.; Selesnick, I.; Cai, G.; Feng, Y.; Sui, X.; Chen, X. Nonconvex sparse regularization and convex optimization for bearing fault diagnosis. IEEE Trans. Ind. Electron. 2018, 65, 7332–7342. [Google Scholar] [CrossRef]
Lv, Y.; Yuan, R.; Song, G. Multivariate empirical mode decomposition and its application to fault diagnosis of rolling bearing. Mech. Syst. Signal Process. 2016, 81, 219–234. [Google Scholar] [CrossRef]
Guo, X.; Chen, L.; Shen, C. Hierarchical adaptive deep convolution neural network and its application to bearing fault diagnosis. Measurement 2016, 93, 490–502. [Google Scholar] [CrossRef]
Wang, Z.; Wang, J.; Kou, Y.; Zhang, J.; Ning, S.; Zhao, Z. Weak fault diagnosis of wind turbine gearboxes based on MED-LMD. Entropy 2017, 19, 277. [Google Scholar] [CrossRef]
Chen, J.; Pan, J.; Li, Z.; Zi, Y.; Chen, X. Generator bearing fault diagnosis for wind turbine via empirical wavelet transform using measured vibration signals. Renew. Energy 2016, 89, 80–92. [Google Scholar] [CrossRef]
Bab, S.; Najafi, M.; Sola, J.F.; Abbasi, A. Annihilation of non-stationary vibration of a gas turbine rotor system under rub-impact effect using a nonlinear absorber. Mech. Mach. Theory 2019, 139, 379–406. [Google Scholar] [CrossRef]
Liang, M.; Bozchalooi, I.S. An energy operator approach to joint application of amplitude and frequency-demodulations for bearing fault detection. Mech. Syst. Signal Process. 2010, 24, 1473–1494. [Google Scholar] [CrossRef]
Ciabattoni, L.; Ferracuti, F.; Freddi, A.; Monteriu, A. Statistical Spectral Analysis for Fault Diagnosis of Rotating Machines. IEEE Trans. Ind. Electron. 2017, 65, 4301–4310. [Google Scholar] [CrossRef]
Tang, B.; Liu, W.; Song, T. Wind turbine fault diagnosis based on Morlet wavelet transformation and Wigner-Ville distribution. Renew. Energy 2010, 35, 2862–2866. [Google Scholar] [CrossRef]
Chen, J.; Li, Z.; Pan, J.; Chen, G.; Zi, Y.; Yuan, J.; Chen, B.; He, Z. Wavelet transform based on inner product in fault diagnosis of rotating machinery: A review. Mech. Syst. Signal Process. 2016, 70, 1–35. [Google Scholar] [CrossRef]
Wang, Y.; Xu, G.; Liang, L.; Jiang, K. Detection of weak transient signals based on wavelet packet transform and manifold learning for rolling element bearing fault diagnosis. Mech. Syst. Signal Process. 2015, 54, 259–276. [Google Scholar] [CrossRef]
Deng, W.; Zhang, S.; Zhao, H.; Yang, X. A novel fault diagnosis method based on integrating empirical wavelet transform and fuzzy entropy for motor bearing. IEEE Access 2018, 6, 35042–35056. [Google Scholar] [CrossRef]
Xiao, M.; Wen, K.; Zhang, C.; Zhao, X.; Wei, W.; Wu, D. Research on fault feature extraction method of rolling bearing based on NMD and wavelet threshold denoising. Shock Vib. 2018. [Google Scholar] [CrossRef]
Pachori, R.B.; Nishad, A. Cross-terms reduction in the Wigner–Ville distribution using tunable-Q wavelet transform. Signal Process. 2016, 120, 288–304. [Google Scholar] [CrossRef]
Ming, Y.; Chen, J.; Dong, G. Weak fault feature extraction of rolling bearing based on cyclic Wiener filter and envelope spectrum. Mech. Syst. Signal Process. 2011, 25, 1773–1785. [Google Scholar] [CrossRef]
Liu, W.Y.; Gao, Q.W.; Ye, G.; Ma, R.; Lu, X.N.; Han, J.G. A novel wind turbine bearing fault diagnosis method based on Integral Extension LMD. Measurement 2015, 74, 70–77. [Google Scholar] [CrossRef]
Li, Y.; Xu, M.; Wang, R.; Huang, W. A fault diagnosis scheme for rolling bearing based on local mean decomposition and improved multiscale fuzzy entropy. J. Sound Vib. 2016, 360, 277–299. [Google Scholar] [CrossRef]
Lv, Y.; Yuan, R.; Wang, T.; Li, H.; Song, G. Health degradation monitoring and early fault diagnosis of a rolling bearing based on CEEMDAN and improved MMSE. Materials 2018, 11, 1009. [Google Scholar] [CrossRef] [PubMed]
Imaouchen, Y.; Kedadouche, M.; Alkama, R.; Thomas, M. A frequency-weighted energy operator and complementary ensemble empirical mode decomposition for bearing fault detection. Mech. Syst. Signal Process. 2017, 82, 103–116. [Google Scholar] [CrossRef]
Bustos, A.; Rubio, H.; Castejón, C.; García-Prada, J. EMD-based methodology for the identification of a high-speed train running in a gear operating state. Sensors 2018, 18, 793. [Google Scholar] [CrossRef] [PubMed]
Kennel, M.B.; Brown, R.; Abarbanel, H.D.I. Determining embedding dimension for phase-space reconstruction using a geometrical construction. Phys. Rev. A 1992, 45, 3403. [Google Scholar] [CrossRef] [PubMed]
Yi, C.; Lv, Y.; Dang, Z.; Xiao, H.; Yu, X. Quaternion singular spectrum analysis using convex optimization and its application to fault diagnosis of rolling bearing. Measurement 2017, 103, 321–332. [Google Scholar] [CrossRef]
Kouchaki, S.; Sanei, S.; Arbon, E.L.; Dijk, D.J. Tensor based singular spectrum analysis for automatic scoring of sleep EEG. IEEE Trans. Neural Syst. Rehabil. Eng. 2015, 23, 1–9. [Google Scholar] [CrossRef] [PubMed]
Liao, Z.; Song, L.; Chen, P.; Guan, Z.; Fang, Z.; Li, K. An effective singular value selection and bearing fault signal filtering diagnosis method based on false nearest neighbors and statistical information criteria. Sensors 2018, 18, 2235. [Google Scholar] [CrossRef] [PubMed]
Wang, H.; Nie, F.; Huang, H.; Makedon, F. Fast nonnegative matrix tri-factorization for large-scale data co-clustering. In Proceedings of the Twenty-Second International Joint Conference on Artificial Intelligence, Barcelona, Spain, 16–22 July 2011. [Google Scholar]
Lin, T.; Zha, H. Riemannian manifold learning. IEEE Trans. Pattern Anal. Mach. Intell. 2008, 30, 796–809. [Google Scholar]
Erichson, N.B.; Voronin, S.; Brunton, S.L.; Kutz, J.N. Randomized Matrix Decompositions Using R. J. Stat. Softw. 2019, 89, 1–48. [Google Scholar] [CrossRef]
Sun, W.; Du, Q. Graph-regularized fast and robust principal component analysis for hyperspectral band selection. IEEE Trans. Geosci. Remote Sens. 2018, 56, 3185–3195. [Google Scholar] [CrossRef]
Dang, C.; Radha, H. RPCA-KFE: Key frame extraction for video using robust principal component analysis. IEEE Trans. Image Process. 2015, 24, 3742–3753. [Google Scholar] [CrossRef]
Li, C.Y.; Zhu, L.; Bao, W.Z.; Jiang, Y.L.; Yuan, C.A.; Huang, D.S. Convex local sensitive low rank matrix approximation. In Proceedings of the 2017 International Joint Conference on Neural Networks (IJCNN), Anchorage, AK, USA, 14–19 May 2017. [Google Scholar]
Morabito, F.C.; Labate, D.; La Foresta, F.; Bramanti, A.; Morabito, G.; Palamara, I. Multivariate multi-scale permutation entropy for complexity analysis of Alzheimer’s disease EEG. Entropy 2012, 14, 1186–1202. [Google Scholar] [CrossRef]
Yi, C.; Lv, Y.; Ge, M.; Xiao, H.; Yu, X. Tensor singular spectrum decomposition algorithm based on permutation entropy for rolling bearing fault diagnosis. Entropy 2017, 19, 139. [Google Scholar] [CrossRef]
Zhang, X.; Liang, Y.; Zhou, J. A novel bearing fault diagnosis model integrated permutation entropy, ensemble empirical mode decomposition and optimized SVM. Measurement 2015, 69, 164–179. [Google Scholar] [CrossRef]
Zheng, J.; Pan, H.; Yang, S.; Cheng, J. Generalized composite multiscale permutation entropy and Laplacian score based rolling bearing fault diagnosis. Mech. Syst. Signal Process. 2018, 99, 229–243. [Google Scholar] [CrossRef]
Ge, M.; Lv, Y.; Yi, C.; Zhang, Y.; Chen, X. A Joint Fault Diagnosis Scheme Based on Tensor Nuclear Norm Canonical Polyadic Decomposition and Multi-Scale Permutation Entropy for Gears. Entropy 2018, 20, 161. [Google Scholar] [CrossRef]
Candès, E.J.; Recht, B. Exact matrix completion via convex optimization. Found. Comput. Math. 2009, 9, 717. [Google Scholar] [CrossRef]
Cai, Z. Weighted nadaraya–watson regression estimation. Stat. Probab. Lett. 2001, 51, 307–318. [Google Scholar] [CrossRef]
Lee, J.; Kim, S.; Lebanon, G.; Singer, Y.; Bengio, S. LLORMA: Local low-rank matrix approximation. J. Mach. Learn. Res. 2016, 17, 442–465. [Google Scholar]
Jones, M.C. Simple boundary correction for kernel density estimation. Stat. Comput. 1993, 3, 135–146. [Google Scholar] [CrossRef]
Liu, G.; Lin, Z.; Yan, S.; Sun, J.; Yu, Y.; Ma, Y. Robust recovery of subspace structures by low-rank representation. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 35, 171–184. [Google Scholar] [CrossRef] [PubMed]
Zhang, H.; Lin, Z.; Zhang, C.; Gao, J. Robust latent low rank representation for subspace clustering. Neurocomputing 2014, 145, 369–373. [Google Scholar] [CrossRef]
Wright, J.; Ganesh, A.; Rao, S.; Peng, Y.; Ma, Y. Robust principal component analysis: Exact recovery of corrupted low-rank matrices via convex optimization. In Proceedings of the Advances in Neural Information Processing Systems 2009, Vancouver, BC, Canada, 7–10 December 2009; pp. 2080–2088. [Google Scholar]
Candès, E.J.; Li, X.; Ma, Y.; Wright, J. Robust principal component analysis? J. ACM (JACM) 2011, 58, 11. [Google Scholar] [CrossRef]
Lv, Y.; Ge, M.; Zhang, Y.; Yi, C.; Ma, Y. A Novel Demodulation Analysis Technique for Bearing Fault Diagnosis via Energy Separation and Local Low-Rank Matrix Approximation. Sensors 2019, 19, 3755. [Google Scholar] [CrossRef] [PubMed]
Hong, M.; Luo, Z.Q. On the linear convergence of the alternating direction method of multipliers. Math. Program. 2017, 162, 165–199. [Google Scholar] [CrossRef]
Lin, Z.; Chen, M.; Wu, L.; Ma, Y. The Augmented Lagrange Multiplier Method for Exact Recovery of Corrupted Low-Rank Matrices; Coordinated Science Laboratory Report no. UILU-ENG-09-2215, DC-247; Coordinated Science Laboratory: Urbana, IL, USA, 2009. [Google Scholar]
Feng, Z.; Ma, H.; Zuo, M.J. Amplitude and frequency demodulation analysis for fault diagnosis of planet bearings. J. Sound Vib. 2016, 382, 395–412. [Google Scholar] [CrossRef]
Randall, R.B.; Antoni, J.; Chobsaard, S. The relationship between spectral correlation and envelope analysis in the diagnostics of bearing faults and other cyclostationary machine signals. Mech. Syst. Signal Process. 2001, 15, 945–962. [Google Scholar] [CrossRef]
Feng, Z.; Zuo, M.J.; Qu, J.; Tian, T.; Liu, Z. Joint amplitude and frequency demodulation analysis based on local mean decomposition for fault diagnosis of planetary gearboxes. Mech. Syst. Signal Process. 2013, 40, 56–75. [Google Scholar] [CrossRef]
Cheng, J.S.; Yu, D.J.; Yu, Y. Research on the intrinsic mode function (IMF) criterion in EMD method. Mech. Syst. Signal Process. 2006, 20, 817–824. [Google Scholar]
Tanner, J.; Wei, K. Normalized iterative hard thresholding for matrix completion. SIAM J. Sci. Comput. 2013, 35, S104–S125. [Google Scholar] [CrossRef]
Chen, S.S.; Donoho, D.L.; Saunders, M.A. Atomic decomposition by basis pursuit. SIAM Rev. 2001, 43, 129–159. [Google Scholar] [CrossRef]
Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]

Figure 1. Illustration of the robust principal component analysis (RPCA); A signal trajectory matrix

X \in ℝ^{n_{1} \times n_{2}}

can be decomposed into a low rank feature component

L \in ℝ^{n_{1} \times n_{2}}

and a sparse component

S \in ℝ^{n_{1} \times n_{2}}

Figure 1. Illustration of the robust principal component analysis (RPCA); A signal trajectory matrix

X \in ℝ^{n_{1} \times n_{2}}

can be decomposed into a low rank feature component

L \in ℝ^{n_{1} \times n_{2}}

and a sparse component

S \in ℝ^{n_{1} \times n_{2}}

Figure 2. Illustration of the proposed LRPCA method; in the different high-dimensional phase space

T (e_{i})

associated with the local selected anchor point

e_{i} = (a_{i}, b_{i})

, the signal trajectory matrix

X \in ℝ^{n_{1} \times n_{2}}

can be decomposed into a low rank component

L_{e_{i}} \in ℝ^{n_{1} \times n_{2}}

and a sparse component

S_{e_{i}} \in ℝ^{n_{1} \times n_{2}}

Figure 2. Illustration of the proposed LRPCA method; in the different high-dimensional phase space

T (e_{i})

associated with the local selected anchor point

e_{i} = (a_{i}, b_{i})

, the signal trajectory matrix

X \in ℝ^{n_{1} \times n_{2}}

can be decomposed into a low rank component

L_{e_{i}} \in ℝ^{n_{1} \times n_{2}}

and a sparse component

S_{e_{i}} \in ℝ^{n_{1} \times n_{2}}

Figure 3. Illustration of the coarse-grained and data segments of the time series with

ε = 2

d = 3

, and

δ = 2

, as well as the all

d! = 6

type of ordinal patterns.

Figure 3. Illustration of the coarse-grained and data segments of the time series with

ε = 2

d = 3

, and

δ = 2

, as well as the all

d! = 6

type of ordinal patterns.

Figure 4. A simulated fault feature signal of a bearing’s inner race; (a) signal waveform; (b) signal spectrum.

Figure 5. The simulated fault feature signal of bearing outer race; (a) signal waveform; (b) signal spectrum.

Figure 6. The MSPE of different types of simulated signals; (a) five types of signal: the two bearing fault feature signals shown as Figure 4 and Figure 5, a harmonic signal with the main frequency of 150 Hz, a shock signal with the main frequency of 150 Hz, and a Gaussian white noise signal; (b) the MSPEs of four feature signals in Figure 6a mixed with Gaussian white noise (signal to noise ratio (SNR) = −5).

Figure 7. The flowchart of the proposed effective fault diagnosis technique via LRPCA and MSPE.

Figure 8. The simulated shock signal and harmonic signal; (a) shock signal waveform; (b) harmonic signal waveform.

Figure 9. Comparison of the de-nosing performance of four methods when white Gaussian noise of varying SNR is added to the multi-component signal.

Figure 10. The simulated multi-component signal contains a strong white Gaussian noise with the SNR of −5 db; (a) signal waveform; (b) signal spectrum.

Figure 11. The analysis result of the wavelet shrinkage denoising method; (a) waveform of the extracted fault feature signal; (b) spectrum of the extracted fault feature signal.

Figure 12. The analysis result of the basis pursuit denoising method; (a) waveform of the extracted fault feature signal; (b) spectrum of the extracted fault feature signal.

Figure 13. The analysis result of the EMD method; (a) waveforms of top 12 IMFs; (b,d) waveform and spectrum of IMF 3; (c,e) waveform and spectrum of IMF 4.

Figure 14. The analysis result of the SSA method; (a) waveform of the one-dimensional component signals; (b) waveform of the extracted fault feature signal; (c) spectrum of the extracted fault feature signal.

Figure 15. Analysis result of the proposed technique; (a) waveform of the one-dimensional component signals; (b) the MSPE of the components; (c) waveform of the extracted fault feature signal; (d) spectrum of the extracted fault feature signal.

Figure 16. Bearing-gear fault experiment table; (a) physical photograph; (b) structural drawing: 1—AC motor; 2—the mounting position of fault bearing; 3—magnetic powder brake; 4—gearbox.

Figure 17. The acquired bearing fault signal; (a) signal waveform; (b) signal spectrum.

Figure 18. Analysis result of the wavelet shrinkage denoising method; (a) waveform of the extracted fault feature signal; (b) spectrum of the extracted fault feature signal.

Figure 19. Analysis result of the basis pursuit denoising method; (a) waveform of the extracted fault feature signal; (b) spectrum of the extracted fault feature signal.

Figure 20. Analysis result of the EMD method; (a) waveforms of top 12 IMFs; (b,d) waveform and spectrum of IMF 1; (c,e) waveform and spectrum of IMF 8.

Figure 21. Analysis result of the SSA method; (a) waveform of the extracted fault feature signal; (b) spectrum of the extracted fault feature signal. The fault feature extraction results of SSA. The peaks of the fault feature frequency (

f_{o}

) and its triple frequency (

3 f_{o}

) were obvious in the spectrum. But there are still many interference peaks and noise, which affect the identification of the fault feature. These above analysis results indicate that neither EMD nor SSA can provide a good fault diagnosis performance for the experimental fault signal.

f_{o}

) and its triple frequency (

3 f_{o}

Figure 22. Analysis result of the proposed technique; (a) waveform of the one-dimensional component signals; (b) the MSPE of the components; (c) waveform of the extracted fault feature signal; (d) spectrum of the extracted fault feature signal.

Table 1. Parameters in simulated fault feature signal of a bearing’s inner race.

$ξ$	$υ_{i}$	I	L	A	$B_{i}$	$C_{i l}$	$φ_{i}$	$φ_{i l}$	$f_{n}$	$f_{r}$	$f_{c}$
800	0.02/ $f_{c}$	250	100	1	0.0004	2/ $l^{2}$	$0^{\circ}$	$0^{\circ}$	2000 Hz	30 Hz	125 Hz

Table 2. Experiment parameters of fault bearing.

Number of Roller Elements	Roller Diameter (mm)	Medium Diameter (mm)	Contact Angle	Rotation Frequency (Hz)	Fault Frequency (Hz)	Sampling Points	Sampling Frequency (Hz)
$z = 17$	$d = 8$	$D = 46$	$α = {14.04}^{\circ}$	$f_{r}$ = 7.225	$f_{o}$ = 51.05	N = 20,000	$f_{s}$ = 10,000

© 2019 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 (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Ge, M.; Lv, Y.; Zhang, Y.; Yi, C.; Ma, Y. An Effective Bearing Fault Diagnosis Technique via Local Robust Principal Component Analysis and Multi-Scale Permutation Entropy. Entropy 2019, 21, 959. https://doi.org/10.3390/e21100959

AMA Style

Ge M, Lv Y, Zhang Y, Yi C, Ma Y. An Effective Bearing Fault Diagnosis Technique via Local Robust Principal Component Analysis and Multi-Scale Permutation Entropy. Entropy. 2019; 21(10):959. https://doi.org/10.3390/e21100959

Chicago/Turabian Style

Ge, Mao, Yong Lv, Yi Zhang, Cancan Yi, and Yubo Ma. 2019. "An Effective Bearing Fault Diagnosis Technique via Local Robust Principal Component Analysis and Multi-Scale Permutation Entropy" Entropy 21, no. 10: 959. https://doi.org/10.3390/e21100959

APA Style

Ge, M., Lv, Y., Zhang, Y., Yi, C., & Ma, Y. (2019). An Effective Bearing Fault Diagnosis Technique via Local Robust Principal Component Analysis and Multi-Scale Permutation Entropy. Entropy, 21(10), 959. https://doi.org/10.3390/e21100959

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

Article Menu

An Effective Bearing Fault Diagnosis Technique via Local Robust Principal Component Analysis and Multi-Scale Permutation Entropy

Abstract

1. Introduction

2. Theory Description

2.1. Notations and Abbreviations

2.2. Decomposing a Signal into Single-Components via LRPCA

2.2.1. RPCA

2.2.2. LRPCA

2.3. Identification of Bearing Fault-Related Signal through MSPE

2.3.1. Basic Theory of MSPE

2.3.2. Mathematical Model of Bearing Fault Feature Signal

2.3.3. Identification of Fault-Related Components through MSPE

2.4. The Process of the Effective Fault Diagnosis Technique

3. Experiments

3.1. Numerical Simulation Experiment

3.2. Experimental Signal Analysis

4. Conclusions

Author Contributions

Funding

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI