Next Article in Journal
Synchronization in Finite-Time of Delayed Fractional-Order Fully Complex-Valued Dynamical Networks via Non-Separation Method
Previous Article in Journal
Bath Engineering Enhanced Quantum Critical Engines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fault Feature-Extraction Method of Aviation Bearing Based on Maximum Correlation Re’nyi Entropy and Phase-Space Reconstruction Technology

1
School of Mechanical and Electrical Engineering, Henan University of Technology, Zhengzhou 450001, China
2
Henan Key Laboratory of Superabrasive Grinding Equipment, Henan University of Technology, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(10), 1459; https://doi.org/10.3390/e24101459
Submission received: 3 September 2022 / Revised: 26 September 2022 / Accepted: 27 September 2022 / Published: 13 October 2022

Abstract

:
To address the difficulty of extracting the features of composite-fault signals under a low signal-to-noise ratio and complex noise conditions, a feature-extraction method based on phase-space reconstruction and maximum correlation Re’nyi entropy deconvolution is proposed. Using the Re’nyi entropy as the performance index, which allows for a favorable trade-off between sporadic noise stability and fault sensitivity, the noise-suppression and decomposition characteristics of singular-value decomposition are fully utilized and integrated into the feature extraction of composite-fault signals by the maximum correlation Re’nyi entropy deconvolution. Verification based on simulation, experimental data, and a bench test proves that the proposed method is superior to the existing methods regarding the extraction of composite-fault signal features.

1. Introduction

As a core component ensuring the safety of aircraft flight, aeroengines have a complex structure and need to operate continuously for a prolonged time under extremely severe conditions such as high temperature, high pressure, high speed, high strength, and variable load. The failure of rolling bearings, as a common component of aeroengines, seriously affects the aircraft flight safety. Hence, the early fault diagnosis of key aeroengine components is of great significance to ensuring flight safety and lowering economic and life losses [1].
Given the tight arrangement and complex structure of aeroengines, it is difficult to arrange vibration-monitoring sensors at the proximal end of the core components. As a result, the dynamic response of the core component failure will be subjected to the modulation of complex transfer paths, as well as the impact of other excitations [2,3]. Moreover, the early bearing fault may appear as a composite fault with multiple coexisting faults. Due to the coupling and interference among different faults and between faults and other excitations, the identification and separation of fault-signal features are hardly achievable, which produces huge challenges regarding the fault diagnosis of aircraft bearings [4].
To separate and extract the composite bearing fault features under complex excitations, it is necessary to deal with three aspects, namely, the selection of the sensitive feature norm, the suppression method of noise and other excitations, and the separation and decoupling method of composite noise, in order to attain a preferable performance.
Kurtosis (Sk), as a sensitive feature capable of detecting instantaneous impacts, has received extensive attention in the fault-diagnosis field in recent years. Its successful application has been seen in the diagnosis of wind-turbine gear faults [5]; the vibration source identification for offshore wind turbines [6]; and the detection of gearbox composite faults [7], planetary gearbox faults [8], and rolling bearing faults [9], achieving a good effect. Antoni [10] elaborated on the relevant theory of Sk and formally gave its mathematical definition, i.e., the energy-normalized fourth-order spectral cumulant. Later, in 2007, a fast kurtogram (Fk) method based on short-time Fourier transform or finite impulse-response bandpass filter was proposed capable of adaptively obtaining the appropriate filter center frequency and bandwidth and effectively extracting the fault features. This method has greatly promoted the Sk application in the fault-diagnosis field. Subsequently, Lee [11] introduced a weighted kurtogram approach, whereas Wang [12] adaptively determined the filter bandwidth and center frequency by maximizing the filtered signal Sk through the incorporation of a right-expanding window. Although the above methods have all improved the Fk method to varying degrees and enhanced the fault-signal monitoring ability, all of them have encountered an unavoidable problem in practical use. That is, the primary function of Sk is to find instantaneous impacts. Thus, they are excessively sensitive to random impulses. Moreover, the inevitable sporadic strong impulse signals in actual production will greatly interfere with the ability of Sk to distinguish the real fault signals, which will even lead to the filtering failure [13]. In recent years, as an extended form of Shannon entropy, Re’nyi entropy (Re), has been widely used in the estimation of signal information content [14]. Boškosk and Juričić [15] proposed a new method to diagnose gearbox faults, using Re’nyi entropy as the characteristic index. Subsequently, it is also widely used in the fault diagnosis of bearings [16,17].
Given the structural peculiarity of the rotating machinery and bearings, the bearing signals are generated essentially by some periodic or cyclic mechanisms, so the dynamic response of bearings has cyclostationary features, and their fault signals exhibit typical repetitive transient characteristics [18]. Deconvolution is effective at eliminating the influence of complex transfer paths and enhancing the impact vibration characteristics of periodic faults. At its core, this method designs a finite impulse-response (FIR) filter with the aim of maximizing the value of the post-filtering signal sensitive characteristic index. For example, the minimum entropy deconvolution (MED) proposed by Wiggins [19], which finds the optimal filter based on the filtering result Sk, enables the iterative extraction of the impact signals from the test signals with unknown transfer paths. With advantages such as requiring less parameter adjustment and possessing fast convergence, MED has been widely applied in the field of mechanical fault-diagnosis [20,21]. However, since MED aims at Sk maximization in the iterative process, the iterations are more inclined to extract unit-impulse signals with larger values, which often fall into the local optimal solution. In 1984, Cabrelli [22] proposed the optimal minimum entropy deconvolution (OMED) and proved that it was an accurate global optimal solution. To improve the sensitivity of the MED technique to periodic impulse signals, McDonald [23] proposed the maximum correlated kurtosis deconvolution (MCKD) in 2012, with which the problem of periodic impulse deconvolution is solved. On the downside, MCKD requires the presence of priori knowledge about the fault cycle, which easily falls into the local optimal solution and is susceptible to strong noise interference. Later, McDonald put forward the multipoint optimal minimum entropy deconvolution adjusut (MOMEDA) in 2017, which eliminated the discontinuity problem between input signals and addressed some problems with OMED and MCKD, achieving preferable effects. However, all of the aforementioned deconvolution methods are considered for single-fault conditions, which are ineffective in the case of impulse signals with different periods caused by composite faults.
Phase-space reconstruction (PSR) is suitable for typical nonlinear and nonstationary mechanical impulse faults, which reflect the dynamic characteristics of the system under various fault states in a high-dimensional space. Despite receiving extensive attention from scholars in the industry [2], the use of sequential deconvolution methods of phase-space reconstruction for fault diagnosis is problematic since the sequential use of multiple methods inevitably produces signal distortion. Moreover, the use of more sequential methods indicates higher unpredictability of inter-method coupling. Singular-value decomposition (SVD), as a phase-space reconstruction algorithm with zero-phase and -time shifts, has been widely applied in extracting features of early weak-fault signals [24]. Zhao [25] discussed the similarity between SVD and the wavelet decomposition in signal processing, believing that SVD has the noise-suppression characteristics, as well as decomposition and extraction characteristics.
On the basis of the above literature analysis, because the kurtosis is too sensitive to the occasional pulse signal, it affects the accuracy of fault identification. A new method is proposed that integrates the phase-space reconstruction technique into the maximum correlation Re’nyi entropy deconvolution; in this method, the Re’nyi entropy, which can better balance the sensitivity of fault and the stability of accidental noise, is chosen as the performance index, the maximum correlation Re’nyi entropy is taken as the optimization objective, and deconvolution is chosen as the optimization method, in order to solve the problem of the low signal-to-noise ratio (SNR) and the composite-fault recognition rate, while the phase-space reconstruction technique is combined with Re’nyi entropy.
The rest of the paper is organized as follows. In Section 2, the concept and process of Re’nyi entropy, the maximum correlation Re’nyi entropy, the maximum correlation Re’nyi entropy deconvolution, and the phase-space reconstruction are introduced. In the third part, the method of phase-space reconstruction combined with maximum correlation Re’nyi entropy deconvolution is introduced. In Section 3 and Section 4, the algorithm is compared with the existing algorithms by using simulation, experimental data, and the bench test. The final conclusions are provided in Section 5.

2. Materials and Methods

2.1. Re’nyi Entropy (Re)

It is assumed that the incomplete probability set of a random event X is p = { p 1 , p 2 , , p n } , and its overall probability sum is ω ( p ) : = i p i 1 . Then, the generalized Re parameterized by order a can be defined as
H α R ( p ) = 1 1 α l o g 2 i p i α i p i   α > 0 , α 1
The parameter α in Rényi entropy can be used to make it more or less sensitive to particular segments of the probability distributions. The exponent α helps to provide flexibility by highlighting the values closer to the edges of the probability distribution [26]. The definition of generalized Re is introduced into the actual vibration detection, and the number set x = { x 1 , x 2 , , x N } is assumed as the discrete observation of the actual vibration process. Let there be a non-negative number set ξ = { ξ 1 , ξ 2 , , ξ n } ; x corresponds to { ξ i }   in a one-to-one manner and satisfies 0 ξ i 1 and i = 1 N ξ i 1 . According to the definition of generalized Re, the set { ξ i } can be regarded as the probability distribution function of random variable x ( t ) , and ξ i is the probability of the instantaneous amplitude x i for x ( t ) . In actual production, normal bearings often produce Gaussian vibration characteristics. In the presence of a fault, the vibrations collected from defective bearings exhibit a non-Gaussian distribution due to the fault-excited relative increase in the number of large-amplitude components. Moreover, the energy of the periodic impulse signal and the level of defect-induced excitation will increase with the defect development, ultimately resulting in the change of H α R value. Hence, Re is able to monitor the bearing-health status.
To ensure the non-negativity of ξ i , its value can be derived from the following formula:
( ξ i ) r = f ( x i ) = y i i y i ,   y i = { | x i μ | , r = a v ( x i μ ) 2 , r = s v
where μ denotes the mean value of { x i } and the subscript r represents the conversion method. Obviously, Formula (2) satisfies two basic conditions of generalized Re: 0 < ξ i < 1 and i ( ξ i ) r 1 . By assuming μ = 0 without a loss of generality, the generalized Re of order α can be obtained as
( H α ) r = 1 1 α l o g 2 ( 1 / N ) i = 1 N y i α ( ( 1 / N ) i = 1 N y i ) α + log 2 N ,   α > 0 , α 1  
We can let
  G M r α = ( 1 / N ) i = 1 N y i α ( ( 1 / N ) i = 1 N y i ) α = { ( 1 / N ) i = 1 N | x i | α ( ( 1 / N ) i = 1 N | x i | ) α , r = a v , α > 0 , ( 1 / N ) i = 1 N ( x i 2 ) α ( ( 1 / N ) i = 1 N x i 2 ) α ,         r = s v , α > 0
It is clear from Formulas (3) and (4) that when measuring the amplitude distribution of { x i } , G M r α is equivalent to ( H α ) r and is the kernel function of ( H α ) r . The two cases of Formula (4) can be written as a unified expression as
M n = ( 1 / N ) i = 1 N x i n ( ( 1 / N ) i = 1 N x i 2 ) n = ( 1 / N ) i = 1 N x i n σ n
As is clear from Formula (5),   G M r α can be regarded as the generalization of Re. Obviously, by choosing different r and α values, various statistical parameters can be derived from G M r α . When r = a v ,     a = 3 , r = s v , a = 3 / 2 , and a = 2 , the following three statistical indicators can be derived separately:
G M a v 3 = ( 1 / N ) i = 1 N | x i | 3 ( ( 1 / N ) i = 1 N | x i | ) 3 = R e ,                                                                                         r = a v ,   a = 3   G M s v 3 / 2 = ( 1 / N ) i = 1 N ( ( x i 2 ) ) 3 / 2 ( ( 1 / N ) i = 1 N ( x i 2 ) ) 3 / 2 = ( 1 / N ) i = 1 N | x i | 3 σ 3 = S r                 r = s v ,   a = 3 / 2 G M s v 2 = ( 1 / N ) i = 1 N ( ( x i 2 ) ) ( ( 1 / N ) i = 1 N ( x i 2 ) ) 2 = ( 1 / N ) i = 1 N x i 4 σ 4 = S k                 r = s v ,   a = 2  
According to Formula (6), G M s v 2 is equivalent to Sk and G M s v 3 / 2 is equivalent to the third-order moment Skewness (Sr). Thus, G M a v 3 has become a new statistical indicator R e , i.e., the narrow-sense Re proposed herein. This means that R e has a similar mathematical expression to S k and Sr. Therefore, R e , S r , and S k can be considered as different expressions derived from the generalized Re, all of which can be explained by generalized Re theory.
Suppose the fault at the inner ring of rolling bearing is
x ( t ) = j M a j e g γ ( t ) c o s [ ω 0 1 g 2 ( γ ( t ) τ j )
where   a j   is the magnitude of the Jth fault shock, g   denotes the attenuation coefficient of the bearing, M is the number of excitation of the bearing shock,   γ ( t )   is the pseudo-cycle time, T denotes the period of the bearing fault shock, τ j   stands for the time delay due to relative slip, and   f e   stands for the fault characteristic frequency. The values of the simulation parameters are shown in Table 1.
Figure 1 depicts the variations of S k , S r , and R e with the fault defect at the inner ring of the rolling bearing. To facilitate understanding and display, the defect evolution is converted into the SNR change. According to Figure 1, S r is insensitive to the defect size, while R e and S k share similar variation trends. Thus, clearly, R e and S k are highly sensitive to the alterations of fault defects, which can rather accurately indicate the bearing faults. Since this phenomenon may present similar variation trends for the outer ring and the rolling element faults, it is not described in detail here.
Figure 2 depicts the variations of three performance indicators S k , S r , and R e —with the sporadic impulse-response. Data on the figure demonstrates that S r has excellent robustness to the sporadic impulse-response, whereas S k is highly sensitive to the sporadic impulse-response. The S k value increases by over three-fold upon sporadic interference, indicating that the aero-engine produces a greater impact on Sk when it has sporadic impulses. R e is somewhat sensitive to sporadic impulses, although the overall variation is not large. As suggested by Figure 1 and Figure 2, R e has a preferable ability to trade off between the bearing defects and the sporadic impulse sensitivity.

2.2. Correlation Re’nyi Entropy (CRe)

Statistical properties of rotating components such as bearings, gears, shafts, and propellers change periodically over time. Thus, the signals generated by such components are called cyclostationary signals. Conventional signal processing can extend and exploit this feature. A signal can be nth-order cyclostationary when its nth-order statistic is periodic [27]. To exploit the cyclostationarity of rotating components, this study defines the correlation Re’nyi entropy (CRe) deconvolution on the basis of Re, with a view to extracting periodic impulse signals. The first-order CRe and Mth-order CRe are defined separately.
C R e 1 ( T ) = n = 1 N ( y n y n T ) 3 2 ( n = 1 N y n ) 3
C R e M ( T ) = n = l N ( m = 0 M y n m T ) 3 2 ( n = 1 N y n ) M + 2
Noise signal, sporadic impulse + noise signal, cosine signal + noise signal, and periodic impulse + noise signal (Figure 3) are defined to analyze S k , R e ,   C S k 1 , and C R e 1 . Since Gaussian noise signals are stable for the same higher-order statistics, the normalized sensitivity ρ i j is defined to explain the sensitivity of the aforementioned indicators to typical signals.
ρ i j = p i j p 0 j i = [ 0 , 1 , 2 , 3 ] = [ a , b , c , d ] ;             j = [ 0 , 1 , 2 , 3 ] = [ S k , R e , C S k 1 , C R e 1 ]
It is clear from Figure 4 that Sk has high sensitivity to both sporadic and periodic impulse signals and is easily affected by sporadic noise. C S k [20] and C R e have good robustness to both harmonic and sporadic impulse signals, and they are highly sensitive to periodic impulse signals, which can thus distinguish the periodic impulse signals. A horizontal comparison reveals that C R e has better sensitivity to C S k under the periodic impulse condition.

2.3. Maximum Correlation Re’nyi Entropy Deconvolution (MCReD)

Given the complex and compact interior of aeroengines, the vibration sensors cannot be arranged near the faulty bearings since the fault signals are easily affected by strong noise and intricate transmission paths. Thus, the dynamic response of faulty bearings can be regarded as the linear convolution of vibration signal and channel. The deconvolution method can effectively eliminate the influence of intricate transmission paths and enhance the fault-impact vibration characteristics. In this study, a novel MCReD-based method is proposed to overcome the limitations of heavily Gaussian and non-Gaussian background noises.
An actual fault signal of the aeroengine bearing will contain multiple components, which can be expressed as
x ( t ) = h d d + h u u + h e e x = [ x 1 x 2 x N ] ,   d = [ d 1 d 2 d N ] ,   u = [ u 1 u 2 u N ] ,   e = [ e 1 e 2 e N ]
where x stands for the vibration signal acquired by remote vibration sensor; d stands for the impulse signal generated by a faulty aero-engine bearing; u represents the interference signal generated by the other aero-engine components; e refers to the background noise; and h d ,     h u , and h e , respectively, stand for the transfer functions corresponding to different inputs.
The core idea of the MCReD algorithm is to find the global optimal FIR filter f by the deconvolution operation on the basis of eliminating the background noise and other interfering elements to the maximum extent, in order to highlight the pulse sequence in the fault signal, it can be expressed mathematically as:
{ y = f ( h u * u ) + f ( h d d ) + f ( h e e ) St.   f ( h u u ) 0                                                                                                                                         where   f = [ f 1   f 2   f L ] T f ( h d d ) + f ( h e e ) 0
From Formula (7) and the discrete signal convolution formula, the following can be deduced:
y n = k = 1 L f k x m ,   m = n k + 1
Without a loss of generality, the relevant conclusion can be illustrated by the Formula (7) for first-order Cre (M = 1). The optimal filter can be obtained by the following formula:
C f       m a x R e 1 ( T ) =       f       m a x n = 1 N ( y n y n T ) 3 2 ( n = 1 N y n ) 3 ,   f = [ f 1   f 2   f L ] T
To solve the filter coefficient of maximum CRe, Formula (13) can be solved and expressed as
d d f k C R e 1 ( T ) = 0 ,   k = 1 , 2 , 3 , L
Calculating the derivative C R e 1 N of Formula (14)’s denominator and substituting the d d f k y n = x m , m = n k + 1 yield:
d d f k C R e 1 N = 3 2 m = 0 1 X m T a     where   X i = [ x 1 + i r ,   x 2 + i r , , x N + i r ] T ,       X r = [ X 1 , X 2 , , X L ]                                                                               a = [ y 1 m T 1 ( y 1 y 1 T   y 1 M T ) 3 2 y 2 m T 1 ( y 2 y 2 T   y 2 M T ) 3 2 y N m T 1 ( y n y n T   y n M T ) 3 2 ]
Similarly, calculating the derivative d d f k C S 1 D of numerator yields
d d f k C R e 1 D = 3 ( n = 1 N y n ) 2 x m
By synthesizing Formulas (14) and (15), the following can be obtained according to Formula (13)
d d f k C R e 1 ( T ) = 3 ( n = 1 N y n ) 2 C R e 1 N 6 b 3 n = 1 N y n ( x m )                                                                                                                 where : b = [ y 1 y 1 T     y 1 M T y 2 y 2 T     y 1 M T           y N y N T   y N M T   ] T
Additionally, since
y = X 0 T f
the arrangement yields
( n = 1 N y n ) 2 m = 0 1 X m T a = 2 b 3 X 0 y
If ( X T X 0 T ) 1 exists, the solution formula for first-order C R e 1 can be obtained as
f = ( n = 1 N y n ) 2 2 b 3 ( X 0 X 0 T ) 1 m = 0 1 X m T a
This formula can be generalized to the Mth-order as:
f = ( n = 1 N y n ) 2 2 b 3 ( X 0 X 0 T ) 1 m = 0 M X m T a

2.4. Phase-Space Reconstruction(PSR)

MCReD can enhance the periodic impulse signals and eliminate the negative effects of transfer function and sporadic impulse-response. However, similar to the case of the MCKD method, its ability to extract composite faults with different periods is unsatisfactory at low SNR.
PSR, as a time-series analysis technique, recovers important system-component information from the high-dimensional space that is extended from a one-dimensional time series. The phase-space trajectory matrix X r composed of original signals is a Hankel matrix, which can be reconstructed by selecting the components of fault information based on the noise-suppression and decomposition characteristics of SVD.
For a phase-space trajectory matrix X r R L × N , regardless of whether its rows and columns are correlated or not, there must be orthogonal matrices U = u 1 , u 2 u n R L × N and V   =   v 1 , v 2 v n R L × N , so that
X r = U S V T = i = 1 i = L δ i u i v i = i = 1 i = L A i
is established. In the formula, the diagonal matrix is S = diag[ δ 1 , δ 2 δ m ], where δ 1   δ 2 ≥ … ≥   δ m > 0. Formula (20) is called the SVD of x(t). The main diagonal element of S is the singular value of the matrix X r . A i is defined as the subspace of matrix X, and the signal x i ( l ) reconstructed by the anti-diagonal method is defined as the sub-signal discretization.
In PSR, different subspaces can be selected for reconstruction to represent different signal characteristics. After the selection of sub-spaces, signals can be reconstructed to obtain a reconstructed phase-space that can represent the signal characteristics.
X r c = i P σ i u i v i
where P stands for the sub-signal screening
Substituting the reconstruction matrix into Formula (19) yields
f = ( n = 1 N y n ) 2 2 b 3 [ X 0 c ( X 0 c ) T ] 1 m = 0 M X m T c a
After the integration of the PSR, substantial priori knowledge required by MCReD will no longer be important. The rather important filter length L and period T in MCKD, for instance, only require suitable value-taking in PSR-MCReD due to the noise-suppression, decomposition, and extraction characteristics of PSR, where no accurate fault eigenperiod is needed. Thus, PSR-MCReD also has a certain blind-solution property.

2.5. Determination of Screening Space

With the PSR algorithm, screening of the subspace containing more fault information for reconstruction is critical. Using the sub-signal Re as the performance indicator, this study carries out reconstruction by selecting the subspace that contains subsignal mapping with a maximum Re value.
  • Phase-space trajectory matrix Xr is subjected to SVD.
  • Ai values of various orders are derived separately, and the anti-diagonal method is used to reconstruct the sub-signals x i ( l ) of various orders for Ai.
  • R e i = R e ( x i ( l ) ) is calculated and arranged in descending order as [Rei1, Rei2, …, Reij], Rei1Rei2 ≥ … ≥ Reij > 0.
  • By adopting the difference method Dj−1 = ReijReij−1 and solving max(Dj−1), the number of elements in the P space is determined to be j1.
  • Screening space P is obtained based on the i mapped by Seij.

3. PSR–MCReD-Based Fault-Diagnosis Method

In the phase-space reconstruction, the key of the phase-space reconstruction algorithm is to filter the subspace that contains more fault information. In this paper, the Re’nyi entropy of the subsignal is taken as the performance index, and the subspace containing the maximum Re’nyi entropy subsignal mapping is selected for reconstruction. The specific flowchart is shown in Figure 5.
  • Fault signal is input.
  • X r c is calculated according to the screening space selection method for PSR.
  • Filter length is assumed as L, and to prevent the filter from falling into local optimal solution, the initial filter is assumed to be f = [ 0   0   1   1     0   0 ] T .
  • y = X 0 a T f , the filtered signal y is calculated.
  • a and b are calculated according to Formula (14).
  • Filter f is updated according to Formula (19).
  • If d d f k C S M ( T ) > ε is determined to be true, step 4 is returned to; if false, the process is ended. To avoid the iteration falling into an infinite loop, the ε in the formula is chosen as a tiny positive number.
  • Final filter order result is calculated by formula   y = X 0 a T f .
  • Envelope spectrum analysis is performed, and fault information is obtained.

4. Results

4.1. Simulation Analysis

Suppose that the faulty bearing has both inner and outer ring faults and is disturbed by signals such as random impacts, discrete harmonics, and Gaussian white noise.
{ x ( t ) = x 1 ( t ) + μ x 2 ( t ) + x 3 ( t ) + x 4 ( t ) + n ( t ) x 1 ( t ) = - + A ( K T 1 ) e ε ( t K T 1 ) ) s i n [ 2 π f s ( t K T 1 + φ k 1 ) ] x 2 ( t ) = - + e ε ( t K T 2 T 2 2 ) ) s i n [ 2 π f s ( t K T 2 T 2 2 + φ k 2 ) ] x 3 ( t ) = i = 1 M 1 D i e ε ( t i T r ) ) s i n [ 2 π f s ( t i T r + φ k 3 ) ] x 4 ( t ) = j = 1 M 2 p j s i n ( 2 π f 4 j t φ k 4 )
x 1 ( t ) denotes the dynamic response of bearing with the inner ring fault, x 2 ( t ) denotes the dynamic response of bearing with the outer ring fault, x 3 ( t ) represents the external random impact interference during the measurement, x 4 ( t ) stands for the discrete harmonics from shafts or other components received by the remote sensors, and n ( t ) represents the Gaussian white noise with SNR= −8. Table 2 details the specific parameter selection.
Figure 6 displays the time domain waveforms of various component signals.
To evaluate the performance, this section compares the proposed FSR-MCReD method with the MCKD and FK methods [26] using simulation-signals, where the MCKD and FSR-MCReD filter lengths are L= 100, and the number of iterations is set at 30.
It is clear from Figure 7 that under complex operating conditions, the Sk value of random impulses is often larger than the periodic impulse sequence of signals, so that the filter fails to perform filtering in the frequency bands selected by the FK method, leading to the failure of fault identification. According to Figure 8, the filtering effect with the MCKD method is greatly compromised under the composite-fault condition, where the effective fault identification is hardly achievable. Figure 9a displays the differential Re result of PSR signals by the FSR-MCReD method when M = 1. As is clear, the differential Re values for two cycles are similar, with the maximum values both being 2. According to Figure 9b, since FSR-MCReD adopts the Re as the optimization condition, which is more robust to periodic simple harmonics and sporadic impulses, the decomposition and noise-suppression characteristics of PSR are fully exploited and integrated into the MCReD calculation, thereby enabling the preferable identification of composite faults.

4.2. Aeroengine Fault-Diagnosis

Bearings are the core components of aeroengines. Since the aeroengines are often tightly arranged, the vibration sensors can hardly be installed at the proximal end, which causes great difficulties in the feature-extraction of fault signals. Aeroengine vibration data were acquired from an accessory gearbox (Safran, France) [28]. Figure 10 displays the structural schematic of the aeroengine and the locations of sensors. At the L5 shaft location, outer ring spalling damage and retainer failure were present for the rolling bearing structure. Vibration sensors were arranged near the L1 and L5 shafts of the gearbox, while the rotational speed sensors were arranged on the L4 shaft. Measured data of vibration sensor 2 was selected for analysis, and the rotational speed and fault eigenfrequency were calculated according to the L4 shaft speed, the numbers of teeth on L4 and L5 shafts, and the bearing parameters of L5 rolling bearing (Table 3). Through comparison with the FK and FSC [29] algorithms, the ability of FSR-MCReD to extract features of aeroengine fault signals was verified.
Figure 11 displays the envelope spectrum analysis results with the FK algorithm. As is clear, the FK algorithm has a good filtering effect. However, filtering errors are caused due to excessive and complex external noise, making it difficult to distinguish the fault features of the outer ring signals. Figure 12 presents the analysis results with the FSC algorithm. Clearly, the FSC algorithm has a good enhancement effect on the periodic signals. However, since it has no filtering effect, the fault features of outer ring signals are covered in substantial irrelevant signal features, making it difficult to diagnose faults accurately. Figure 13 displays the envelope analysis results with the FSR-MCReD. As is clear, after integration of PSR into the MCReD algorithm, both the signal filtering and decomposition effects are improved, thus enabling accurate fault-diagnosis.

4.3. Validation of Composite-Fault Experimental Data

The accelerated life cycle test data of bearings provided by Professor Lei’s team from Xi’an Jiaotong University [30] enriches the experimental data on the fault-diagnosis and performance-degradation research. Under the experimental conditions of composite outer and inner ring faults, this study used the experimental data collected at t = 20 min, which were not clearly characterized by the CMS index, to theoretically calculate the fault frequency f i = 171   Hz ,   f o = 109   Hz . Further, a comparative analysis was made against the deconvolution algorithms OMEDA and MOMEDA [2], in order to demonstrate the advantages of our proposed algorithm in solving composite faults.
Figure 14 presents the envelope spectrum analysis results with the OMED, the MOMEDA, and our FSR-MCRD methods. As demonstrated by the results, both OMED and MOMEDA can preferably analyze the outer ring faults having obvious failures, with the MOMEDA exhibiting a better filtering effect since it addresses the local optimal solution problem in the OMED method. The common defect of the two, however, is that the inner ring faults are not displayed preferably. The method proposed herein, on the other hand, allows for better extraction of the fault signals in composite faults due to the incorporation of FSR’s noise-suppression and decomposition characteristics into the MCRD method. Thus, it has better diagnostic ability for composite faults.

4.4. Experimental Verification

Regarding the experimental setup, the rotating machinery fault test bench from Spectra Quest, the acceleration sensor from Yangzhou Kedong, the displacement sensor from Bently, and the LMS data acquisition system were adopted. The rotational speed was set at 2700 rpm, and the faulty bearing model was MB ER-12 K. Composite-fault bearings with a ball fault, an inner ring fault, and an outer ring fault were set up separately. The rotational speed was kept constant at 2700 rpm, and the vibration signals of acceleration sensors were collected by the LMS SCADAS mobile data acquisition system at a sampling frequency of 25.6 kHz. Table 4 details the eigenfrequency calculations of various faults, whereas Figure 15 depicts the experimental panorama. MCKD, OMED, MOMOEDA, and FSR-MCRD were separately used for comparison.
Figure 16 displays the envelope spectrum analysis results with the aforementioned four methods. As is clear, the MCKD method has a good effect on signal feature-extraction under single-fault conditions, which, however, easily falls into a local optimum. Thus, it hardly achieves feature extraction or enhancement for other coexisting periodic impulse signals. The non-iterative method OMED exhibits the worst effect of fault feature-extraction, which can hardly diagnose faults. In the case of MOEDA, the results of fault feature-extraction are incomplete since it is an improved version of OMED, which does not consider the condition of multiple fault coexistence. With the FSR-MCReD method, the simultaneous existence of multiple faults is considered, and Re, which has good stability for both sporadic impulses and harmonics, is used as the performance index. Moreover, PSR is integrated for performing deconvolution, so that the subspace signals containing more fault information are extracted, thus preferably achieving the feature extraction under composite-fault condition. However, it is noteworthy that the second- and third-order inner ring fault features are lost during feature extraction since part of the subspace is truncated in PSR.

5. Conclusions

This paper indicatively introduces the Re’nyi entropy, which can better balance the fault sensitivity and the stability of the accidental pulse signal, into the fault diagnosis of an aero-engine. The signal transmission noise and accidental noise are solved by MCReD. In order to better decouple the complex fault signals and improve the extraction ability of MCReD under low SNR conditions, the PSR technique is integrated into the MCReD, not a simple combination. The proposed method has achieved good results in simulation and experimental verification.
  • In this study, a novel method for extracting the fault-signal features of aeroengine core components based on maximum CRe deconvolution integrating PSR is proposed. Through simulation analysis, experimental data, and a bench test, the advantages of the proposed method over the existing signal feature-extraction methods are verified regarding the extraction of complex fault-signal features under low SNR and complex noise interference.
  • The mathematical formulas for Sk, Sr, and the narrow-sense Re proposed in this study are deduced based on the generalized Re. The three performance indicators can all be regarded as different expressions for generalized Re. As revealed by a simulation experiment on the correlations of the three with fault sensitivity and sporadic noise stability, both Re and Sk are highly sensitive to the faults, and Re is more stable to sporadic noise.
  • Inspired by CK, a definition of CRe is given, which has a better suppression effect on the sporadic and harmonic noises. By deriving the solution formula of MCReD, a non-iterative method for solving the MCReD of the optimal filter is proposed.
  • For the impulse signals with different periods under composite-fault conditions, a maximum Re-based subspace-filtering method is proposed by integrating the PSR technique into the MCReD calculation, which utilizes the noise-suppression and decomposition characteristics of SVD.
  • Through simulation, experimental data validation, and bench test verification, the method proposed herein is proven to be fairly effective at extracting composite-fault signal features under low SNR and complex noise conditions.
In this paper, firstly, the maximum correlation Re’nyi entropy deconvolution algorithm is similar to the maximum correlation kurtosis deconvolution algorithm, which is dependent on prior knowledge. Secondly, the theoretical derivation of the decoupling characteristics of the phase-space technology needs to be improved, and the selection method of the screening space of the singular-value decomposition needs to be further optimized. Finally, some experiments on the application of PSR-MCReD are carried out, and more experiments need to be completed in the future work.

Author Contributions

Resources, H.Z.; Software, Y.L.; Writing—original draft, Z.Z.; Writing—review & editing, B.L. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the supports from the National Natural Science Foundation of China under Grant Nos. 12072106 and U1604254.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AbbreviationImplication
SkKurtosis
FkFast kurtogram
FIRFinite impulse response
MEDMinimum entropy deconvolution
OMEDOptimal minimum entropy deconvolution
MCKDMaximum correlated kurtosis deconvolution
MOMEDAMultipoint optimal minimum entropy deconvolution adjustment
PSRPhase-space reconstruction
SVDSingular-value decomposition
ReRe’nyi entropy,
SrSkewness
CReCorrelation Re’nyi entropy
CSkCorrelation Kurtosis
MCReDMaximum correlation Re’nyi entropy deconvolution
SNR Signal-to-noise ratio

References

  1. Ma, J.; Han, S.; Li, C.; Zhan, L.; Zhang, G.-Z. A New Method Based on Time-Varying Filtering Intrinsic Time-Scale Decomposition and General Refined Composite Multiscale Sample Entropy for Rolling-Bearing Feature Extraction. Entropy 2021, 23, 451. [Google Scholar] [CrossRef] [PubMed]
  2. Zhao, Q.; McDonald, G.L. Multipoint Optimal Minimum Entropy Deconvolution and Convolution Fix: Application to vibration fault detection. Mech. Syst. Signal Process. 2017, 82, 461–477. [Google Scholar]
  3. Zhang, X.; Miao, Q.; Zhang, H.; Wang, L. A parameter-adaptive VMD method based on grasshopper optimization algorithm to analyze vibration signals from rotating machinery. Mech. Syst. Signal Process. 2018, 108, 58–72. [Google Scholar] [CrossRef]
  4. Qin, Y.; Jin, L.; Zhang, A.; He, B. Rolling Bearing Fault Diagnosis with Adaptive Harmonic Kurtosis and Improved Bat Algorithm. IEEE Trans. Instrum. Meas. 2021, 70, 3508112. [Google Scholar] [CrossRef]
  5. Barszcz, T.; Randall, R.B. Application of spectral kurtosis for detection of a tooth crack in the planetary gear of a wind turbine. Mech. Syst. Signal Process. 2009, 23, 1352–1365. [Google Scholar] [CrossRef]
  6. Dong, X.; Lian, J.; Wang, H. Vibration source identification of offshore wind turbine structure based on optimized spectral kurtosis and ensemble empirical mode decomposition. Ocean Eng. 2019, 172, 199–212. [Google Scholar] [CrossRef]
  7. Wang, T.; Chu, F.; Han, Q.; Kong, Y. Compound faults detection in gearbox via meshing resonance and spectral kurtosis methods. J. Sound Vib. 2017, 392, 367–381. [Google Scholar] [CrossRef]
  8. Smith, W.A.; Fan, Z.; Peng, Z.; Li, H.; Randall, R.B. Optimised Spectral Kurtosis for bearing diagnostics under electromagnetic interference. Mech. Syst. Signal Process. 2016, 371–394. [Google Scholar] [CrossRef] [Green Version]
  9. Xiang, J.; Zhong, Y.; Gao, H. Rolling element bearing fault detection using PPCA and spectral kurtosis. Measurement 2015, 75, 180–191. [Google Scholar] [CrossRef]
  10. Antoni, J. Fast computation of the kurtogram for the detection of transient faults. Mech. Syst. Signal Process. 2007, 21, 108–124. [Google Scholar] [CrossRef]
  11. Lee, J.H.; Kim, D.H.; Shin, Y.H. Hyperbolic localization of incipient tip vortex cavitation in marine propeller using spectral kurtosis. Mech. Syst. Signal Process. 2018, 110, 442–457. [Google Scholar] [CrossRef]
  12. Wang, Y.; Ming, L. An adaptive SK technique and its application for fault detection of rolling element bearings. Mech. Syst. Signal Process. 2011, 25, 1750–1764. [Google Scholar] [CrossRef]
  13. Zhang, K.; Chen, P.; Yang, M.; Song, L.; Xu, Y. The Harmogram: A periodic impulses detection method and its application in bearing fault diagnosis. Mech. Syst. Signal Process. 2022, 165, 108374. [Google Scholar] [CrossRef]
  14. Baraniuk, R.G.; Flandrin, P.; Janssen, A.J.; Michel, O.J. Measuring time-frequency information content using the Renyi entropies. IEEE Trans. Inf. Theory 2001, 47, 1391–1409. [Google Scholar] [CrossRef] [Green Version]
  15. Boškoski, P.; Juričić, Đ. Fault detection of mechanical drives under variable operating conditions based on wavelet packet Rényi entropy signatures. Mech. Syst. Signal Process. 2012, 31, 369–381. [Google Scholar] [CrossRef]
  16. Boškoski, P.; Gašperin, M.; Petelin, D.; Juričić, Đ. Bearing fault prognostics using Rényi entropy based features and Gaussian process models. Mech. Syst. Signal Process. 2015, 52–53, 327–337. [Google Scholar] [CrossRef]
  17. Tao, B.; Zhu, L.; Han, D. An alternative time-domain index for condition monitoring of rolling element bearings—A comparison study. Reliab. Eng. Syst. Saf. 2007, 92, 660–670. [Google Scholar] [CrossRef]
  18. Antoni, J. The infogram: Entropic evidence of the signature of repetitive transients. Mech. Syst. Signal Process. 2016, 74, 73–94. [Google Scholar] [CrossRef]
  19. Wiggins, R.A. Minimum entropy deconvolution. Geophys. Prospect. Pet. 1980, 16, 21–35. [Google Scholar] [CrossRef]
  20. Endo, H.; Randall, R.B. Enhancement of autoregressive model based gear tooth fault detection technique by the use of minimum entropy deconvolution filter. Mech. Syst. Signal Process. 2007, 21, 906–919. [Google Scholar] [CrossRef]
  21. Yang, F.; Kou, Z.; Wu, J.; Li, T. Application of Mutual Information-Sample Entropy Based Med-IceemdanDe-Noising Scheme for Weak Fault Diagnosis of Hoist Bearing. Entropy 2018, 20. [Google Scholar]
  22. Cabrelli, C.A. Minimum entropy deconvolution and simplicity: A noniterative algorithm. Geophysics 2012, 50, 394–413. [Google Scholar] [CrossRef] [Green Version]
  23. Mcdonald, G.L.; Zhao, Q.; Zuo, M.J. Maximum correlated Kurtosis deconvolution and application on gear tooth chip fault detection. Mech. Syst. Signal Process. 2012, 33, 237–255. [Google Scholar] [CrossRef]
  24. Miao, Y.; Zhao, M.; Lin, J. Application of an im-proved maximum correlated kurtosis deconvolution method for fault diagnosis of rolling element bearings. Mech. Syst. Signal Process. 2017, 2, 173–195. [Google Scholar] [CrossRef]
  25. Guo, M.; Li, W.; Yang, Q.; Zhao, X.; Tang, Y. Amplitude filtering characteristics of singular value decomposition and its application to fault diagnosis of rotating machinery. Measurement 2020, 154, 107444. [Google Scholar] [CrossRef]
  26. Singh, J.; Darpe, A.K.; Singh, S.P. Bearing damage assessment using Jensen-Rényi Divergence based on EEMD. Mech. Syst. Signal Process. 2017, 87, 307–339. [Google Scholar] [CrossRef]
  27. Zhao, X.; Ye, B. Singular value decomposition packet and its application to extraction of weak fault feature. Mech. Syst. Signal Process. 2016, 70, 73–86. [Google Scholar] [CrossRef]
  28. Antoni, J.; Xin, G.; Hamzaoui, N. Fast computation of the spectral correlation. Mech. Syst. Signal Process. 2017, 92, 248–277. [Google Scholar] [CrossRef]
  29. Antoni, J.; Griffaton, J.; André, H.; Avendaño-Valencia, L.D.; Bonnardot, F.; Cardona-Morales, O.; Castellanos-Dominguez, G.; Daga, A.P.; Leclere, Q.; Vicuña, C.M.; et al. Feedback on the Surveillance 8 challenge: Vibration-based diagnosis of a Safran aircraft engine. Mech. Syst. Signal Process. 2017, 97, 112–144. [Google Scholar] [CrossRef]
  30. Wang, B.; Lei, Y.; Li, N.; Li, N. A Hybrid Prognostics Approach for Estimating Remaining Useful Life of Rolling Element Bearings. IEEE Trans. Reliab. 2018, 69, 401–412. [Google Scholar] [CrossRef]
Figure 1. S k ,   S r , R e schematic diagram of the variation with defects.
Figure 1. S k ,   S r , R e schematic diagram of the variation with defects.
Entropy 24 01459 g001
Figure 2. S k , S r , R e sensitivity to accidental noise; (a)   S k ,   S r ,   R e under noise condition; (b) S k ,   S r ,    R e under noise and occasional pulse conditions.
Figure 2. S k , S r , R e sensitivity to accidental noise; (a)   S k ,   S r ,   R e under noise condition; (b) S k ,   S r ,    R e under noise and occasional pulse conditions.
Entropy 24 01459 g002
Figure 3. Typical signal of (a) noise condition, (b) occasional pulse with noise, (c) harmonic signal with noise, and (d) periodic pulse with noise.
Figure 3. Typical signal of (a) noise condition, (b) occasional pulse with noise, (c) harmonic signal with noise, and (d) periodic pulse with noise.
Entropy 24 01459 g003
Figure 4. Normalized sensitivity.
Figure 4. Normalized sensitivity.
Entropy 24 01459 g004
Figure 5. Flowchart of PSR-MCReD.
Figure 5. Flowchart of PSR-MCReD.
Entropy 24 01459 g005
Figure 6. Component signal time domain waveform of (a) inner fault x 1 ( t ) time domain waveform, (b) outer fault x 2 ( t ) time domain waveform, (c) outer fault x 2 ( t ) time domain waveform, and (d) random pulse   x 3 ( t ) time domain waveform.
Figure 6. Component signal time domain waveform of (a) inner fault x 1 ( t ) time domain waveform, (b) outer fault x 2 ( t ) time domain waveform, (c) outer fault x 2 ( t ) time domain waveform, and (d) random pulse   x 3 ( t ) time domain waveform.
Entropy 24 01459 g006
Figure 7. Fast spectral kurtosis results of (a) fast kurtosis, and (b) fast spectral kurtosis envelope spectrum.
Figure 7. Fast spectral kurtosis results of (a) fast kurtosis, and (b) fast spectral kurtosis envelope spectrum.
Entropy 24 01459 g007
Figure 8. MCKD results of (a) time-domain results of MCKD, and (b) FSR-MCReD envelope spectrum.
Figure 8. MCKD results of (a) time-domain results of MCKD, and (b) FSR-MCReD envelope spectrum.
Entropy 24 01459 g008
Figure 9. FSR-MCReD results of (a) determination of filter space, and (b) FSR-MCReD envelope spectrum.
Figure 9. FSR-MCReD results of (a) determination of filter space, and (b) FSR-MCReD envelope spectrum.
Entropy 24 01459 g009
Figure 10. A schematic diagram of an aero-engine structure.
Figure 10. A schematic diagram of an aero-engine structure.
Entropy 24 01459 g010
Figure 11. Fast spectral kurtosis results of (a) fast kurtosis, and (b) fast spectral kurtosis envelope spectrum.
Figure 11. Fast spectral kurtosis results of (a) fast kurtosis, and (b) fast spectral kurtosis envelope spectrum.
Entropy 24 01459 g011
Figure 12. FSC results of (a) fast spectral correlation, and (b) FSC enhanced envelope spectra.
Figure 12. FSC results of (a) fast spectral correlation, and (b) FSC enhanced envelope spectra.
Entropy 24 01459 g012
Figure 13. FSR-MCReD results of (a) determination of filter space, and (b) FSR-MCReD envelope.
Figure 13. FSR-MCReD results of (a) determination of filter space, and (b) FSR-MCReD envelope.
Entropy 24 01459 g013
Figure 14. Three methods result of (a) OMED envelope spectra (b) MOMED envelope spectra, and (c) MOMED envelope spectra.
Figure 14. Three methods result of (a) OMED envelope spectra (b) MOMED envelope spectra, and (c) MOMED envelope spectra.
Entropy 24 01459 g014aEntropy 24 01459 g014b
Figure 15. Experimental layout. 1. Motor; 2. coupling; 3. bearing housing I; 4. principal axis; 5. rotors; 6. bearing housing II; 7. bearing I; and 8. bearing II.
Figure 15. Experimental layout. 1. Motor; 2. coupling; 3. bearing housing I; 4. principal axis; 5. rotors; 6. bearing housing II; 7. bearing I; and 8. bearing II.
Entropy 24 01459 g015
Figure 16. Four methods of envelope results of (a) MCKD envelope spectra, (b) OMED envelope spectra, (c) MOMEDA envelope spectra, and (d) FSR-MCReD envelope spectra.
Figure 16. Four methods of envelope results of (a) MCKD envelope spectra, (b) OMED envelope spectra, (c) MOMEDA envelope spectra, and (d) FSR-MCReD envelope spectra.
Entropy 24 01459 g016
Table 1. Simulation-signal parameters.
Table 1. Simulation-signal parameters.
f 0 / Hz a j Mg f e / Hz τ j
500.710.7102%T
Table 2. Simulation-signal parameters.
Table 2. Simulation-signal parameters.
ParameterImplicationValueParameterImplicationValue
f s Sampling   frequency 12,000 Hz D i Amplitude of shockRandom variable
f i Inner   ring 213   Hz T 3 Random impact timeRandom variable
f o Outer   ring 143   Hz p j Harmonic amplitude p 1 = 0.03 ;       p 2 = 0.04
f r Rotating   frequency 2000   Hz f 4 j Harmonic frequency f 41 = 45.2 ;   f 42 = 56.7
μ Ratio   coefficient 0.3 φ k i Initial phase Random   [ π , π ]
Table 3. Characteristic frequency of engine bearings.
Table 3. Characteristic frequency of engine bearings.
Rotating Frequency: frRolling Element: fbOuter Ring: foInner Ring: fiCage: ft
182.9650.51419.51843.578.8
Table 4. Characteristic frequency of bearings.
Table 4. Characteristic frequency of bearings.
Rotating Frequency: frRolling Element: fbOuter Ring: foInner Ring: fiCage: ft
4589.6137222.717.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Z.; Liu, B.; Liu, Y.; Zhang, H. Fault Feature-Extraction Method of Aviation Bearing Based on Maximum Correlation Re’nyi Entropy and Phase-Space Reconstruction Technology. Entropy 2022, 24, 1459. https://doi.org/10.3390/e24101459

AMA Style

Zhang Z, Liu B, Liu Y, Zhang H. Fault Feature-Extraction Method of Aviation Bearing Based on Maximum Correlation Re’nyi Entropy and Phase-Space Reconstruction Technology. Entropy. 2022; 24(10):1459. https://doi.org/10.3390/e24101459

Chicago/Turabian Style

Zhang, Zhen, Baoguo Liu, Yanxu Liu, and Huiguang Zhang. 2022. "Fault Feature-Extraction Method of Aviation Bearing Based on Maximum Correlation Re’nyi Entropy and Phase-Space Reconstruction Technology" Entropy 24, no. 10: 1459. https://doi.org/10.3390/e24101459

APA Style

Zhang, Z., Liu, B., Liu, Y., & Zhang, H. (2022). Fault Feature-Extraction Method of Aviation Bearing Based on Maximum Correlation Re’nyi Entropy and Phase-Space Reconstruction Technology. Entropy, 24(10), 1459. https://doi.org/10.3390/e24101459

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