Next Article in Journal
Modeling and Performance Analysis of Opportunistic Link Selection for UAV Communication
Previous Article in Journal
Walking Secure: Safe Routing Planning Algorithm and Pedestrian’s Crossing Intention Detector Based on Fuzzy Logic App
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimal Parameter Selection Method for MOMEDA Based on EHNR and Its Spectral Entropy

1
Fauclty of Information Engineering and Automation, Kunming University of Science and Technology, Kunming 650500, China
2
Yunnan Key Laboratory of Artificial Intelligence, Kunming University of Science and Technology, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(2), 533; https://doi.org/10.3390/s21020533
Submission received: 13 December 2020 / Revised: 6 January 2021 / Accepted: 10 January 2021 / Published: 13 January 2021
(This article belongs to the Section Fault Diagnosis & Sensors)

Abstract

:
As a vital component widely used in the industrial production field, rolling bearings work under complicated working conditions and are prone to failure, which will affect the normal operation of the whole mechanical system. Therefore, it is essential to conduct a health assessment of the rolling bearing. In recent years, Multipoint Optimal Minimum Entropy Deconvolution Adjusted (MOMEDA) is applied to the fault feature extraction for rolling bearings. However, the algorithm still has the following problems: (1) The selection of fault period T depends on prior knowledge. (2) The accuracy of signal denoising is affected by filter length L. To solve the limitations, an improved MOMEDA (IMOMEDA) method is proposed in this paper. Firstly, the envelope harmonic-to-noise ratio (EHNR) spectrum is adopted to estimate the fault period of MOMEDA. Then, the improved grid search method with EHNR spectral entropy as the objective function is constructed to calculate the optimal filter length used in the MOMEDA. Finally, a feature extraction method based on the improved MOMEDA (IMOMEDA) and Teager-Kaiser energy operator (TKEO) is applied in the field of rolling bearing fault diagnosis. The effectiveness and generalization performance of the proposed method is verified through comparison experiment with three data sets.

1. Introduction

Rolling bearings are vital components of mechanical systems and have been widely used in rotating machinery. However, rolling bearings often work in more demanding environments for a long time, such as corrosive, high temperature, large impact and high load conditions, which increases the uncertainty of the mechanical system, resulting in various fault defects of the bearings [1]. In order to avoid the environmental pollution, economic loss and even safety problems caused by rolling bearing shutdown, the fault diagnosis method based on signal analysis is widely used in industrial production, in which the bearing condition is recognized according to the effective information collected by the vibration signal [2,3], pressure signal [4], temperature signal [5] and acoustic signal [6,7].
The fault diagnosis technique based on vibration signal has been widely used in practice because of its intuitive operation and reliable diagnosis results, where the vibration signal of the bearing is monitored by the acceleration sensors installed in the appropriate direction of the bearing seat or the box, and the signal is analyzed and processed to determine the bearing condition. Furthermore, some advanced measurement systems are applied in rolling bearings. Wang [8] designed a data acquisition system of high speed railway bearing based on ZigBee technology. Gao et al. [9] designed a wireless sensor vibrating signal acquisition and processing system, and the researchers could collect and analyze the signal timely and conveniently. Li et al. [10] proposed a smart metasurface shaft (SMST) to achieve single-sensor identification of vibration sources, which showed potential applications in rotating machinery condition monitoring. Verstraete et al. [11] presented multi-sensor fusion systems for integrated remaining useful life prognostic capabilities. Meng et al. [12] developed a fully enclosed bearing-structured self-powered rotation sensor (SPRS) for multi-tasking motion measurement. Furthermore, the device displayed high precision, broad range, quick response, and excellent durability.
Fault feature extraction based on signal analysis is critical after the reliable vibration data source is available. In industrial production, when the bearing fails to work, the collected vibration signal will generate periodic transient pulses at a certain characteristic frequency [13]. However, in fact, in addition to the pulses caused by fault defects, there are also some modulation harmonics caused by external excitations, such as axial unbalance, misalignment, critical speed, structural resonance and uncertain factors on vibration signal transmission paths, which make the spectrum more complex and increases the difficulty of bearing fault diagnosis [14]. Therefore, recovering and enhancing effective fault period impulse components from vibration signals is of great significance and importance for fault feature extraction. Some signal processing methods have been successfully employed in capturing fault-related transient components in the non-stationary signal for rolling bearings [15], such as spectral kurtosis [16,17], wavelet transform theory [18] and so on. Unfortunately, most of the methods mainly concentrate on emphasize random impulse with large amplitude and thus ignores the periodicity of fault signal.
To tackle the limitation, the deconvolution principle has been introduced for bearing fault periodic composition extraction, which achieved satisfactory results. In this line, Wiggins [19] proposed the minimum entropy deconvolution (MED) in 1978, which took kurtosis as the objective function and the optimal inverse filter was iteratively found to enhance the transient ingredients in the original signals. In 2007, Randall successfully applied the MED to fault detection [20]. Since then, the MED has been widely used in the fault feature extraction of rotating machinery [21,22,23].
However, the MED method aims at maximizing the kurtosis value, and thus it is easily affected by strong noise. In the iterative process, only the pulse component in the vibration signal can be enhanced and this defect limits the development of MED method [24]. To address this issue, McDonald et al. [25] proposed the maximum correlation kurtosis deconvolution (MCKD), which took the correlation kurtosis as the objective function and iteratively searched for the optimal reverse filter to enhance the periodic transient components in the original signal. The MCKD algorithm has achieved remarkable results in the fault feature extraction of rotating machinery such as rolling bearings [26,27], gearboxes [25], and motor rotors [28]. Nonetheless, there are still some issues in the practical application of MCKD. In fact, its accuracy is affected by four parameters: filter length, fault period, the number of displacement and the number of iterations [29]. Moreover, the MCKD can only address the single period pulse signal.
To extract the continuous periodic pulse signal, McDonald et al. [30] proposed the multipoint optimal minimum entropy deconvolution adjusted (MOMEDA) method. This method is a non-iterative deconvolution process and utilizes infinite pulse sequence as the optimization object function to search for the optimal filter directly. With these advantages, the MOMEDA has been successfully applied to the rotating machinery fault characteristics extraction [31,32]. Ma et al. [33] introduced MOMEDA to detect planet bearings faults and the experiments show that MOMEDA presented specific applicability compared with the traditional kurtogram filtering method. Li et al. [34] suggested a method based on the MOMEDA and long and short time memory (LSTM) for fault diagnosis. Zhu et al. [35] introduced a method based on the MOMEDA and TKEO for the rolling bearing fault diagnosis with good experimental qualities. However, the following weaknesses of the MOMEDA method are also well noted in the community: (1) the performance of the MOMEDA mainly depends on the prior knowledge of fault period T. It is necessary to identify the periodic impulses signal through the multipoint kurtosis spectrum [30] in the default fault period T search interval. However, when the original signal is affected by complex background noise, the multipoint kurtosis spectrum may not provide correct fault period information according to the location of maximum multipoint kurtosis. (2) The effect of noise reduction for MOMEDA is also affected by the filter length L. Hence, unreasonable filter length L will weaken or enhance the energy of the original signal, and then cause misdiagnosis. To address the above-mentioned issues, scholars have introduced some optimization algorithms [36,37] to determine the fault period T and filter length L of MOMEDA, such as particle swarm optimization (PSO) algorithm [27,38], grasshopper optimization algorithm (GOA) [39] and so on. These parameter optimization algorithms mostly take kurtosis, envelope spectrum kurtosis (ESK) and other similar indexes as objective functions to improve MOMEDA. However, the methods also have some limitations. These algorithms are usually applied to multi-objective optimization, i.e., training parameters such as fault period T, filter length L and window function simultaneously, which make the implementation process complex and time-consuming. Therefore, it is necessary to explore other effective methods to quickly calculate the fault period T and filter length L of MOMEDA for improving the deconvolution performance.
An improved MOMEDA (IMOMEDA) method is investigated for fault feature extraction in this paper. The proposed method incorporates the envelope harmonic-to-noise ratio (EHNR) [40] into the MOMEDA scheme. By calculating the EHNR spectrum of the original signal, the fault period T can be quickly determined without using the prior information. Furthermore, kurtosis spectral entropy is proposed as a grid search objective function to search for the optimal filter length L of MOMEDA, for which the single objective optimization can greatly simplify the procedures [24]. However, kurtosis spectral entropy tends to emphasize random impulse with large amplitude (it inherits the essential characteristic of kurtosis) and thus ignores the periodicity of fault signals. Therefore, EHNR spectral entropy is presented as the target function to overcome this shortcoming and obtain the improved grid search method to select the optimal filter length L of MOMEDA. Compared with the existing methods, the overall approach has the following advantages:
(1) A new index-based EHNR spectrum is proposed to calculate fault period T of MOMEDA and require less prior knowledge and calculation time.
(2) EHNR spectral entropy is introduced to construct the objective function of grid search method and select the optimal filter length L of MOMEDA.
(3) A novel feature extraction method-based IMOMEDA-TKEO is presented and applied in the field of rolling bearing fault diagnosis.
(4) The effectiveness and generalization performance of the presented method is verified through comparison experiment with three data sets.
The remaining part of the papers is organized as follows: Section 2 introduces the basic principles of MOMEDA and EHNR. Section 3 elaborates on the implementation process of IMOMEDA-TKEO. Section 4 presents the experimental results and analysis. Section 5 gives the discussion and summary.

2. Theoretical Background

2.1. MOMEDA Method

Suppose x is the vibration signal, y is the impulse sequence, h is the transfer function, e is noise, and “*” is the convolution product. The discrete-time response signal is as follows [41]:
x = h * y + e
The MOMEDA algorithm aims at finding an optimal filter f in a non-iterative way to realize the reconstruction of the original impulse sequence y and try to reduce the impact of noise. The major process of MOMEDA algorithm is given as follows:
y = f * x = k = 1 N - L f k x k + L - 1
Or in matrix form:
y = X 0 T f
X 0 = [ x 1 x 2 x 3 ... x N 0 x 1 x 2 ... x N 1 0 0 x 1 ... x N 2 ... ... ... ... ... 0 0 0 ... x N L + 1 ]
Among them, k = 1, 2,…, NL, according to the periodicity of the fault signal for rotating machinery, MOMEDA takes the multi D-norm of the filtered signal as the objective function, the concept of MOMEDA is presented in Equations (5) and (6):
MDN ( y , t ) = 1 | | t | | t T y | | y | |
MOMEDA = m a x f MDN ( y , t ) = m a x f t T y | | y | |
The target vector t is a constant vector that defines the location and weightings of the goal impulses to be deconvolved, which is suitable for the features extraction [19]. When t is completely consistent with original impulse signal y, the multi D-norm reaches the maximum value and the deconvolution effect works the best. Therefore, the non-iterative process of Equation (6) is formulated as follows:
d d f ( t T y | | y | | ) = 0
d d f ( t T y | | y | | ) = y - 1 ( t 1 M 1 + t 2 M 2 + ... + t k M k ) - y - 3 t T y X 0 y = 0
where f = f1, f2, f3,…fL, t = t1, t2, t3,… tN−L. Mk is the kth sequence pulse to be deconvolution, it can be expressed as M k = [ x k + L - 1 x k + L - 2 x k ] , Suppose X0 = [M1,M2,…MK], the Equation (8) can be simplified as:
y - 1 X 0 t - y - 3 t T y X 0 y = 0
t y y 2 X 0 y = X 0 t
We then substitute y = X 0 T f into Equation (10), and then calculate:
t y y 2 f = ( X 0 X 0 T ) - 1 X 0 t
The f = ( X 0 X 0 T ) 1 X 0 t can be denoted as a set of optimal filters, and then reconstruct the original impulse signal y:
Y = [ y 1 y 2 y M ] y = X 0 T f
Finally, the multipoint kurtosis (MKurt) is defined as follows:
MKurt = ( n = 1 N - L t n 2 ) 2 n = 1 N - L ( t n y n ) 4 / ( n = 1 N - L t n 8 ( n = 1 N - L y n 2 ) 2 )
However, the construction of target vector t needs to take T as a prior knowledge. In practical engineering, it is usually impossible to meet this condition due to the complex working environments, which in turn limits the application of the MOMEDA method. In addition, the periodic pulse energy is determined by the filter length L. If the parameter setting is unreasonable, it will weaken or enhance the energy of the periodic impulse signal which will lead to misdiagnosis. Therefore, how to calculate the fault period T and filter length L for the MOMEDA method quickly is one of the major issues, aiming at improving its performance.

2.2. EHNR Principle

In the field of condition-based maintenance, a large number of indicators are introduced to determine whether the measured signal is informative or not. These indicators are widely used, including kurtosis, entropy to L2/L1 norm, Gini index and other norm indexes [42]. Most of the indexes mainly concentrate on the peak value, non­­-gaussianity distribution and other general statistical distribution, thus ignoring the particularity of mechanical signals. Therefore, this paper introduces the EHNR method. According to the specific characteristics of rolling bearing fault signal, we can determine the parameter T of MOMEDA quickly.
The calculation process of EHNR is as follows [40]:
(1) The envelope signal Envx(t) of vibration signal x(t) is obtained by Hilbert transform and the DC component is removed.
x ^ ( t ) = H { x ( t ) } = 1 π - x ( τ ) t - τ d τ
E n v x ( t ) = x 2 ( t ) + x ^ 2 ( t )
E n v x ( t ) = E n v x ( t ) - mean ( E n v x ( t ) )
(2) The autocorrelation function of Envx(t) is calculated as:
r E n v x ( τ ) = E n v x ( t ) E n v x ( t + τ ) dt
(3) The maximum position of the autocorrelation function of the signal is found in the lag domain, where the point is shown in Figure 1.
Then, the EHNR can be defined as Equation (18):
E H N R = r E n v x ( τ m a x ) r E n v x ( 0 ) - r E n v x ( τ m a x )
where τmax is the time lag that makes the autocorrelation spectrum of Envx(t) reach its local maximum, r E n v x ( τ m a x ) is the energy of harmonic, r E n v x ( 0 ) is the total energy of envelope signal.

2.3. Setting Filter Length L Based on Grid Search Method

The grid search [43] is an exhaustive search method for specifying the accurate parameter values. According to the definite step size, the grids are divided in a certain range and the constraint function values on each grid point are calculated. Then, the independent variables that make the objective function value optimal (maximum or minimum) are selected as the solution of the problem.
In this paper, the EHNR spectral entropy is proposed as the objective function and the expression is shown in Equation (19):
E H N R _ E ( x ) = E H N R ( x ) E s
Among them, i = 1 N P i ( S R ( ω ) ) = 1 and
E s = - i = 1 N P i ( S R ( ω ) ) ln P i ( S R ( ω ) )
where S R ( ω ) is the envelope spectrum of the original signal, N is the data length, EHNR(x) represents the EHNR value of the input signal, and EHNR_E(x) represents the EHNR spectral entropy of the signal. For the periodic pulse signal, the envelope spectrum is concentrated on the low-frequency region which results in a smaller Es value. EHNR(x) is a function to describe the periodic change of faults. The larger the EHNR(x) is, the stronger the periodicity. Therefore, EHNR_E(x) can reflect the uniformity of the signal period.
In this paper, EHNR_E(x) is used as the objective function of grid search method, and then the optimal filter length L is determined.

2.4. TKEO Algorithm

The Teager-Kaiser energy operator is a nonlinear operator that can track the instantaneous energy of signal. The extraction of envelope by TKEO is to calculate the three adjacent sampling points of the measured waveform. TKEO can track the change of the measured signal waveform in real time. The algorithm has the advantages of excellent time resolution, simple and fast calculation. The specific mathematical expression and implementation process of the TKEO algorithm was proposed in [44].

3. The Proposed Method

Inspired by the above-introduced methods, in this section, an IMOMEDA-TKEO rolling bearing fault feature extraction method is proposed. The framework of the proposed method is shown in Figure 2. The general implementation procedure is summarized as follows:
Step 1: The autocorrelation function r E n v x ( τ ) = E n v x ( t ) E n v x ( t + τ ) dt of rolling bearing vibration signal x(t) is obtained by the autocorrelation analysis.
Step 2: The envelope harmonic-to-noise spectrum EHNR(x) of signal x(t) is obtained and the fault period T is calculated.
Step 3: According to [45], the filter length L meets the principle of L > 2fs /fc, where fs is the sampling frequency and fc is the resonance frequency. We initialize the input parameters of the grid search, set the search range of the filter length L to (250, 4000), the search step length as 50, and then the fault period T is obtained by Step 2.
Step 4: The parameters L are updated by EHNR spectral entropy EHNR_E(L). When EHNR_E(L) reaches maximum, the optimal filter length L is obtained.
Step 5: The parameters of MOMEDA algorithm are set according to step 2 and step 4, and then the original signal x(t) is deconvolved to obtain the analysis signal xcov(t).
Step 6: The signal xcov(t) is demodulated by TKEO to get its Teager-Kaiser energy spectrum, and then complete fault feature extraction.

4. Experimental and Comparative Analysis

In order to verify the effectiveness of the IMOMEDA-TKEO method, the experimental database of rolling bearing from Case Western Reserve University (CWRU) [46], National Aeronautics and Space Administration (NASA) [47] which are authoritative in the field of bearing fault diagnosis and self-made experiment platforms (or KUST-SY for short) are used for experimental analysis. Comparative experiments among IMOMEDA-TKEO, MOMEDA-TKEO and fast kurtogram (FK) [48] are carried out to show the superiority of the proposed method. In each case, we utilize the same input data to complete each comparative experiment. The brief overview of the comparative experiment is shown in Table 1. In each case, we utilize the same input data to complete each comparative experiment.
(a) Compared with the other experiments, the Test 1 demonstrates the effectiveness of the IMOMEDA-TKEO method.
(b) The comparison between Test 2 and Test 3 proves that it is necessary and effective to optimize the filter length of MOMEDA.
(c) Compared with fast kurtogram, the other tests verify that the MOMEDA method can enhance effective fault period impulse component from vibration signals.
(d) The demonstration of pseudo period in Test 3 further proves the importance to calculate fault period T exactly.

4.1. Case 1: CWRU Data Analysis

Figure 3 shows the experiment platform and bearings. Some specific measured parameters are presented in Table 2. The test-bed consists of a two horsepower motor, a torque sensor/decoder, a power tester and an electronic controller. The vibration signals were collected by accelerometers, where the accelerometers are installed at 12 o’clock position at drive end of the motor and fan end of motor housing, respectively. The collected signals were stored in a 16-channel DAT recorder [46] and post-processed in MATLAB 2018a on an IntelI CITM) i5-7300HQ CPU @ 2.50 GHz Laptop. The motor speed is 1797rpm. The sampling frequency fs is 12 kHz, and the number of data points N is 4096.
Outer race defect frequency:
BPFO = Z 2 ( 1 - d D cos θ ) × f r
Inner race defect frequency:
BPFI = Z 2 ( 1 + d D cos θ ) × f r
According to the above parameters and formulas, BPFO = 107.36 Hz and BPFI = 162.19 Hz are calculated, respectively. Next, we use the datasets in CWRU to complete the comparative experiments.

4.1.1. Outer Race Fault Feature Extraction

Figure 4 shows the time-domain and frequency-domain analysis results of the original outer race signal, respectively. It is seen from Figure 4 that these waveforms cannot directly obtain the information related to the outer race fault. Further analysis is required to determine the fault characteristic frequency and other details. Therefore, IMOMEDA-TKEO, MOMEDA-TKEO and the FK methods are proposed to complete subsequent analysis.
(A) IMOMEDA-TKEO method;
(a) Using the EHNR spectrum to calculate fault period T.
By carrying out the autocorrelation analysis of outer ring, the autocorrelation function (AC) spectrum of the original signal is obtained as shown in Figure 5. It can be seen from Figure 5 that the AC spectrum can reflect the periodic change of fault signal, which provides the basis for the estimation of fault period T. In order to estimate T, we further calculate the EHNR spectrum of signal as shown in Figure 6. It can be seen from Figure 6 that the fault signal shows obvious periodic change and gradually decays with time, which meets the trend of fault signal. Then, the fault period T = 112 is obtained.
(b) Using the EHNR spectral entropy to get filter length L
Given the fault period T = 112, the filter length L of the MOMEDA is quickly determined by the grid search method. Firstly, the filter length L satisfies L > 2fs /fc, where fs is the sampling frequency and fc is the resonance frequency [45]. To accurately cover the whole fault frequency band, the initialization L search range is (250, 4000) and the search step is 50. Secondly, the EHNR spectral entropy is proposed as the target function. Finally, the trend of EHNR spectral entropy is obtained by constantly updating the parameter L as shown in Figure 7. When obtaining the maximum output of EHNR_E(L), the optimal filter length is 2050.
(c) Outer race fault feature extraction
According to (a) and (b), the fault period T and filter length L of MOMEDA are calculated as 112 and 2050, respectively. Then, the periodic impulse signal xcov(t) of the original signal x(t) is extracted by using the MOMEDA method, as shown in Figure 8. Then, we demodulate xcov(t) with TKEO and get the corresponding spectrum as shown in Figure 9. It can be seen from Figure 9 that the characteristic frequency (105.5 Hz), with its frequency doubling, approaches BPFO (2~9BPFO). Then, it can be inferred that the outer race fault has occurred.
(B) MOMEDA-TKEO method
According to [38] and [30], we analyze the outer race fault signal with the fault period T = (90,134) and filter length L = 500. The multi-point kurtosis diagram is shown in Figure 10. When the multi-point kurtosis reaches the maximum value, the fault period T is 111.5. The optimal deconvolution signal xcov(t) is presented in Figure 11.
The xcov(t) is demodulated by the TKEO and the corresponding spectrum is shown in Figure 12 where the result shows the characteristic frequency (108.4 Hz), with its frequency doubling, approaches BPFO (2~9BPFO). Further, it is inferred from the original signal that there exists an outer ring fault.
(C) FK method
Fast kurtogram (FK) [48] is a time-frequency analysis method for spectral kurtosis that shows sufficient efficiency in transient fault detection and has been widely used in fault diagnosis field. Therefore, FK is used as a reference for comparative analysis in this paper. The FK and fault characteristic spectrum are shown in Figure 13 and Figure 14, respectively. The result shows the characteristic frequency (107.2 Hz) with its frequency doubling (2~6BPFO) can be extracted.
Compared with the MOMEDA-TKEO method and FK method, the amplitude obtained by the IMOMEDA-TKEO is clearer. Test 3 is carried out through the filter length optimized by the proposed method, without changing the fault period T search interval of traditional MOMEDA. Then, xcov(t) and its Teager-Kaiser energy spectrum are shown in Figure 15 and Figure 16, respectively.
It can be seen from the Mkurt spectrum that the search interval of fault period T directly affects the denoising effect of MOMEDA. However, the interval of T depends on prior knowledge and has trial and error. When the original signal is affected by complex background noise, the Mkurt spectrum may not provide correct fault period information according to the location of maximum multipoint kurtosis. For example, we get the pseudo fault period as T = 97.2. Furthermore, the optimal periodic pulse signal xcov(t) cannot be obtained, as shown in Figure 17. The deconvolution effect is not ideal. The wrong fault characteristic frequency is 123 Hz and its spectrum is shown in Figure 18.
In order to further confirm the effectiveness of the proposed method quantitatively, a new signal-to-noise ratio (SNR) [49] is used as an index for Teager-Kaiser energy spectrum to complete the quantitative comparative analysis among the IMOMEDA-TKEO, MOMEDA-TKEO and the FK. The new SNR is defined as follows.
S N R f d = 10 log 10 i = 1 M p ( i f d ) f = 0 M f d p ( f ) - i = 1 M p ( i f d )
where fd is the fault characteristic frequency, M is the multiple of the fault characteristic frequency; p(f) is the Teager-Kaiser energy spectrum [50]. It can be seen from Table 3 that the SNRfd index of the presented method is the highest, which indicates that the proposed method has better ability in weak feature representation.

4.1.2. Inner Race Fault Feature Extraction

Figure 19 shows the time-domain and frequency-domain analysis results of the original inner race signal, respectively. It is seen from Figure 19 that these waveforms cannot directly obtain the information related to the inner race fault. Further analysis is required to determine the fault characteristic frequency and other details. Therefore, IMOMEDA-TKEO and MOMEDA-TKEO are proposed to complete subsequent analysis.
The specific experimental process of IMOMEDA-TKEO, MOMEDA-TKEO and the FK method is the same as Section 4.1.1. The experimental comparison results are shown in Figure 20 and Figure 21.
The proposed method can achieve comparable effects as the MOMEDA-TKEO method as shown in Figure 20a,b. They both can extract the inner race fault features from the original vibration signals. It further proves the effectiveness of two methods.
Under the condition that the fault period T is invariant, filter length L is obtained by the proposed method and presented in Figure 20c. It can be seen from Figure 20c that the amplitude of TKEO spectrum apparently increases, which indicates filter length L directly affects the deconvolution performance of the MOMEDA method. In addition, it can be seen from Figure 20d that the selection of fault period T directly affects the performance of the method.
FK analysis results are shown in Figure 21. Compared with Figure 20, similar conclusions can be obtained, but the fault characteristic frequency with its frequency doubling characteristic obtained by IMOMEDA-TEO method are more prominent and clearer than those obtained by FK method.
From Table 4, the SNRfd for the IMOMEDA method is higher, which further implies the IMOMEDA has the qualitative and quantitative advantage over the original MOMEDA algorithm and FK method.

4.2. Case 2: NASA Data Analysis

In order to verify the good generalization ability of the IMOMEDA-TKEO method, the further experiments were carried out with NASA bearing datasets.
The test-bed shown in Figure 22 is composed of four bearings, rotating shaft, two acceleration sensors and electronic control unit. The test bearings are located in the motor shaft. Some accelerometers sensors are placed to the motor housing. The motor speed is 2000 rpm and the digital data is collected by sampling with the frequency of 20 kHz. The detailed measured parameters of the bearing are shown in Table 5.
From the bearing parameters shown in Table 5, according to Equation (19) and Equation (20), the fault characteristic frequency BPFO = 236.4 Hz for outer race and BPFI = 296.93 Hz for inner race can be calculated, respectively, and the data length N is 4096 points.
Experimental comparative analysis is complete by utilizing the vibration signals of outer race and inner race. Furthermore, the vibration signal time and frequency analysis is presented in Figure 23. It is seen from Figure 23 that these waveforms cannot directly determine the fault characteristic frequency and other details. In order to further analyze the fault information of the bearing, the same implementation process as Section 4.1 is taken to achieve the deconvolution results of the outer race and inner race, respectively, as shown in Figure 24 and Figure 25.
The experimental results of Figure 24 and Figure 25 show that the proposed IMOMEDA-TKEO method effectively extracts fault features of outer race and inner race and achieves significantly better processing performance than the MOMEDA-TKEO method.
FK based method can extract fault features of outer race and inner race from the original vibration signals. Unfortunately, it can be seen from Figure 26 that the obtained fault characteristic frequency is not obvious (the amplitude is far less than the other two methods) and contains more noise components. Compared with Figure 24 and Figure 25, the filtering performance is poorer than the IMOMEDA-TKEO method.
In order to intuitively compare the effectiveness of the method, the SNRfd index is calculated and shown in Table 6. By comparing and analyzing the data in the Table 6, the conclusions similar to Section 4.1 can be drawn.

4.3. Application Testing

In order to further prove the validity of the proposed method, the practical test data from the self-made KUST-SY experiment platform is selected for subsequent analysis. The platform is able to simulate different failure conditions of rolling bearings. The structure of the test rig is shown in Figure 27, which consists of driving motor, rotating shaft, support bearing, hydraulic loading system, accelerometers and tested bearing. The device can realize radial hydraulic loading and the rotation speed is controlled by computer. Accelerometers are installed in the vertical and horizontal directions of the test bearing base to collect the bearing vibration signal. Because the load is applied in the horizontal direction, the accelerometer placed in this direction is capable of capturing more fault information of the tested bearings. Therefore, the horizontal vibration signals are selected to analyze bearing fault characteristics.
The failure behaviors of bearing system are diverse, including stripping, wear, fracture, crack, electrical erosion, indentation for the inner race and outer race. As displayed in Figure 28, in this experiment, inner ring crack and out ring crack of the bearing have been mainly simulated considering the existing experimental conditions. The raceway surface cracks (the fault diameter is 0.2 mm) of the inner ring and outer ring have an artificial linear cutting with fiber laser cutter. The sampling frequency is 51.2 kHz and the motor speed is 1797 rpm. The sampling interval is 20 s and the sampling time is 10 s. The tested bearing model is 6205-2RSH, and its relevant parameters are shown in Table 7. According to Table 7, the fault characteristic frequency BPFO = 107.22 Hz for the outer race and BPFI = 162.33 Hz for the inner race can be calculated, respectively.
The time and frequency analysis is shown in Figure 29 for the vibration signal of the inner ring and outer ring (The length of data-set N is 32768). We cannot get valuable information related to fault characteristic frequency from Figure 29, which is the same as previous cases. Therefore, the fault information is determined by the same implementation process with Section 4.1 and plotted in Figure 30, Figure 31 and Figure 32, respectively.
It can be seen from Figure 30 and Figure 31 that the fault features extracted by the proposed IMOMEDA-TKEO method are more obvious in the Teager energy spectrum. At the same time, it can be seen from Figure 30d and Figure 31d that the improper period T will lead to erroneous results in the MOMEDA method.
The fault characteristic frequency for the inner race and outer race extracted by FK is shown in Figure 32. Comparing with the fault characteristic spectrum extracted by MOMEDA based method as shown in Figure 30 and Figure 31, it can be seen that the fault characteristic frequencies can be accurately obtained by FK. However, there are many interference components in the spectrum, which result in the faint characteristic amplitude. In particular, the outer race seems to be more seriously disturbed by noise, and only three times of double frequency is extracted, while the other two methods can extract nine times of double frequency. In addition, through the comparison of performance indicator SNRfd presented in Table 8, it is also found that the IMOMEDA-TKEO method can achieve better performance than the MOMEDA-TKEO method and FK method in noise reduction.

4.4. Comparative Analysis with Three Cases

Through the comparative experimental analysis of IMOMEDA-TKEO, MOMEDA-TKEO and FK, the following can be summarized:
(a) Combining with three case studies, the experimental results show that the IMOMEDA-TKEO method is better than the other method described in Table 1, which enhances the application prospect of the proposed method.
(b) In this experiment, we select execution time and iterations to evaluate the complexity of the optimal parameter selection method and all the experiments are programmed in MATLAB 2018a on an Intel(R) Core(TM) i5-7300HQ CPU @ 2.50 GHz Laptop with 16.00 GB RAM. The parameter optimization time and iterations of IMOMEDA algorithm and the main program execution time of the comparative methods are shown in Table 9 and Table 10, respectively. It is seen from Table 9 that we can calculate fault period T of MOMEDA quickly with little time. Compared with the traditional MOMEDA method in which the parameters are set to the default value, the grid search method based on EHNR spectral entropy for selecting optimal filter length L does consume some time. However, optimization algorithm itself is not the main reason of time-consuming. The main reason is that as the iterations increases, the MOMEDA filter time keeps accumulating, leading to more training time of filter length. The purpose of this paper is to accurately calculate the fault period T of MOMEDA and select the optimal filtering length L of MOMEDA. The rationality and necessity of MOMEDA parameter selection are demonstrated on three different datasets. Although the calculation needs a certain time, it is within the acceptable range.
(c) It can be seen from Table 10 that the FK method and the traditional MOMEDA method (default filter length and fault period) take less time. However, when dealing with abnormal or unknown vibration signals, the default parameters of MOMEDA are obviously not feasible. It is verified that the time-consuming by the trial and error to determine the MOMEDA parameters is almost the same as that of the proposed method. However, the result of unreasonable parameter setting is much worse than IMOMEDA, such as Figure 24d and Figure 25d, which also proves the significance of optimizing MOMEDA parameters. The proposed IMOMEDA-TKEO method can effectively estimate the period T of vibration signal, and set the filter length L reasonably by using the improved grid search method. The execution time is within the acceptable range.

5. Discussion

(a) For the traditional MOMEDA method, the range of T has to be set by a trial and error method. Furthermore, the range of T selected in the above experiments is the best result obtained by multiple attempts. If the selected range changes, the pseudo period will cause misdiagnosis. The proposed method introduces the EHNR spectrum to calculate the signal period T quickly and optimizes the filter length L by the improved grid search method. It transforms the iterative optimization problem into a single dimension optimization problem. Hence, the proposed method increases the operation efficiency and improves the deconvolution performance of MOMEDA method.
(b) In the proposed method, the search range of the filter length L is set as [34]. Only the lower limit of the filter length L is given, while there is no upper limit, which makes the search range larger and the algorithm run for a long time. According to the specific characteristics of the signal, optimizing the search range can improve the efficiency of the algorithm. In addition, if we can obtain the corresponding relationship and change rule between the MOMEDA filter length L and period T on the basis of (a), it will further improve the operation efficiency and the practical application value of the proposed method.
(c) Compared with other types of signal characteristics, the practical test vibration signals for inner race from the rolling bearing fault simulation experiment platform are not only disturbed by noise, but also by harmonic signal of the outer race and rolling element. Under strong background noise, the IMOMEDA-TKEO method has been demonstrated to significantly outperform traditional signal processing methods in terms of noise reduction performance and fault feature representation robustness, which boosts the generalization ability.

6. Conclusions

In this paper, a rolling bearing fault feature extraction method based on the IMOMEDA-TKEO is proposed. Some conclusions are obtained as follows:
(1) The EHNR method is introduced to obtain the EHNR spectrum of the original signal and solves the problem that the fault period T of MOMEDA depends on prior knowledge.
(2) The EHNR spectral entropy is proposed as the objective function, and the improved grid search method is introduced to solve the filter length L selection in MOMEDA.
(3) Comparative experimental analysis of the proposed method with the MOMEDA-TKEO method and the FK method are completed. The proposed method can achieve the best filtering performance, which verifies the effectiveness and feasibility of the IMOMEDA-TKEO.
In the proposed method, the upper limit of the filter length L is not given, which makes the execution time and iterations of the algorithm take a certain time (it is within the acceptable range). In future studies, we will focus on improving the calculation accuracy of fault period T for complex signals and optimizing the search range of filter length L. We will also explore the change rule between the filter length L and fault period T to improve the performance of MOMEDA.

Author Contributions

Conceptualization, J.M.; funding acquisition, X.W. methodology, J.M.; project administration, X.W.; resources, X.W. and J.M.; software, Z.L. and X.L.; supervision, J.M. and X.W.; validation, Z.L.; visualization, Z.L.; writing original draft, Z.L.; writing—review and editing, Z.L. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

The National Natural Science Foundation of China (Grant No. 51765022 and 61663017) and Science & Research Program of Yunnan Province (Grant No.2019FD042).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available on Reference [46] and Reference [47].

Acknowledgments

The contributions of all my coauthors are gratefully acknowledged. Moreover, this manuscript is supported by the National Natural Science Foundation of China (No.51765022 and 61663017) and Science & Research Program of Yunnan Province (No.2019FD042).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MOMEDA-multipoint optimal minimum entropy deconvolution adjustedTFMSR-time-delayed feedback monostable stochastic resonance
IMOMEDA-improved MOMEDALSTM-long and short time memory
EHNR-envelope harmonic-to-noise ratioPSO-particle swarm optimization
GOA-grasshopper optimization algorithmCVS-continuous vibration separation
MED-minimum entropy deconvolutionESK-envelope spectrum kurtosis
FK-fast kurtogramSTFT-Short Time Fourier Transform
MCKD-maximum correlation kurtosis deconvolutionNASA-National Aeronautics and Space Administration
AC-autocorrelation functionCWRU-Case Western Reserve University
TKEO-Teager-Kaiser energy operatorBPFO-Outer ring defect frequency
SNR-signal-to-noise ratioBPFI-Inner ring defect frequency

References

  1. Lei, Y.; Lin, J.; He, Z.; Zuo, M.J. A review on empirical mode decomposition in fault diagnosis of rotating machinery. Mech. Syst. Signal Process. 2013, 35, 108–126. [Google Scholar] [CrossRef]
  2. Cerrada, M.; Sánchez, R.-V.; Li, C.; Pacheco, F.; Cabrera, D.; De Oliveira, J.L.V.; Vásquez, R.E. A review on data-driven fault severity assessment in rolling bearings. Mech. Syst. Signal Process. 2018, 99, 169–196. [Google Scholar] [CrossRef]
  3. Bastami, A.R.; Vahid, S. A comprehensive evaluation of the effect of defect size in rolling element bearings on the statistical features of the vibration signal. Mech. Syst. Signal Process. 2021, 151, 1–20. [Google Scholar]
  4. Li, J.; Meng, Z.; Yin, N.; Pan, Z.; Cao, L.; Fan, F. Multi-source feature extraction of rolling bearing compression measurement signal based on independent component analysis. Measurement 2021, 172, 1–10. [Google Scholar] [CrossRef]
  5. Wang, H.; Ni, G.; Chen, J.; Qu, J. Research on rolling bearing state health monitoring and life prediction based on PCA and Internet of things with multi-sensor. Measurement 2020, 157, 1–9. [Google Scholar] [CrossRef]
  6. Rzeszucinski, P.; Orman, M.; Pinto, C.T.; Tkaczyk, A.; Sulowicz, M. Bearing Health Diagnosed with a Mobile Phone: Acoustic Signal Measurements Can be used to Test for Structural Faults in Motors. IEEE Ind. Appl. Mag. 2018, 24, 17–23. [Google Scholar] [CrossRef]
  7. Orman, M.; Rzeszucinski, P.; Tkaczyk, A.; Krishnamoorthi, K.; Pinto, C.T.; Sulowicz, M. Bearing fault detection with the use of acoustic signals recorded by a hand-held mobile phone. In Proceedings of the 2015 International Conference on Condition Assessment Techniques in Electrical Systems (CATCON), Bangalore, India, 10–12 December 2015; pp. 252–256. [Google Scholar]
  8. Wang, Q. Design and application of high-speed railway bearing vibration signal acquisition system based on ZigBee technology. Electron. Des. Eng. 2018, 26, 97–100. [Google Scholar]
  9. Gao, J.; Hao, L.; Li, Y.; Guo, H.; Hou, X. Development of Wireless Sensor Vibrating Signal Acquisition and Processing System Based on VC ++. Mach. Tool Hydraul. 2012, 40, 69–72. [Google Scholar]
  10. Chong, L.; Tianxi, J.; Qingbo, H.; Zhike, P. Smart metasurface shaft for vibration source identification with a single sensor. J. Sound Vib. 2021, 493, 1–18. [Google Scholar]
  11. Verstraete, D.; Droguett, E.L.; Modarres, M. A Deep Adversarial Approach Based on Multi-Sensor Fusion for Semi-Supervised Remaining Useful Life Prognostics. Sensors 2019, 20, 176. [Google Scholar] [CrossRef] [Green Version]
  12. Meng, X.S.; Li, H.Y.; Zhu, G.; Wang, Z.L. Fully enclosed bearing-structured self-powered rotation sensor based on electrification at rolling interfaces for multi-tasking motion measurement. Nano Energy 2015, 12, 606–611. [Google Scholar] [CrossRef]
  13. Huang, J.; Chen, G.; Shu, L.; Wang, S.; Zhang, Y. An Experimental Study of Clogging Fault Diagnosis in Heat Exchangers Based on Vibration Signals. IEEE Access 2016, 4, 1800–1809. [Google Scholar] [CrossRef] [Green Version]
  14. Howard, I.M. A Review of Rolling Element Bearing Vibration "Detection, Diagnosis and Prognosis"; NASA STI/Recon Technical Report N: Sydney, Australia, 1994; pp. 1–10.
  15. Radoslaw, Z.; Walter, B. Application of Adaptive Filtering for Weak Impulsive Signal Recovery for Bearings Local Damage Detection in Complex Mining Mechanical Systems Working under Condition of Varying Load. Solid State Phenom. 2011, 180, 250–257. [Google Scholar]
  16. Jérôme, A. The spectral kurtosis: A useful tool for characterising non-stationary signals. Mech. Syst. Signal Process. 2006, 20, 282–307. [Google Scholar]
  17. Tomasz, B.; Robert, B.R. 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]
  18. Jimeng, L.; Hui, W.; Xiangdong, W.; Yungang, Z. Rolling bearing fault diagnosis based on improved adaptive parameterless empirical wavelet transform and sparse denoising. Measurement 2020, 152, 1–14. [Google Scholar]
  19. Ralph, A.W. Minimum entropy deconvolution. Geoexploration 1978, 16, 21–35. [Google Scholar]
  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. Chen, H.; Wang, J.; Tang, B.; Li, J. Helicopter rolling bearing hybrid faults diagnosis using minimum entropy deconvolution and Teager energy operator. J. Vib. Shock. 2017, 36, 45–50, 73. [Google Scholar]
  22. Cheng, Y.; Zhou, N.; Zhang, W.; Wang, Z. Application of an improved minimum entropy deconvolution method for railway rolling element bearing fault diagnosis. J. Sound Vib. 2018, 425, 53–69. [Google Scholar] [CrossRef]
  23. Zhang, L.; Hu, N. Fault diagnosis of sun gear based on continuous vibration separation and minimum entropy deconvolution. Measurement 2019, 141, 332–344. [Google Scholar] [CrossRef]
  24. Wang, Z.; Du, W.; Wang, J.; Zhou, J.; Han, X.; Zhang, Z.; Huang, L. Research and application of improved adaptive MOMEDA fault diagnosis method. Measurement 2019, 140, 63–75. [Google Scholar] [CrossRef]
  25. 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]
  26. Ding, C.; Zhao, M.; Lin, J.; Jiao, J. Multi-objective iterative optimization algorithm based optimal wavelet filter selection for multi-fault diagnosis of rolling element bearings. ISA Trans. 2018, 82, 199–215. [Google Scholar] [CrossRef]
  27. Cheng, Y.; Wang, Z.; Zhang, W.; Huang, G. Particle swarm optimization algorithm to solve the deconvolution problem for rolling element bearing fault diagnosis. ISA Trans. 2019, 90, 244–267. [Google Scholar] [CrossRef]
  28. Hong, L.; Liu, X.; Zuo, H. Compound fault diagnosis of rotating machinery based on adaptive maximum correlated kurtosis deconvolution and customized multiwavelet transform. Mech. Syst. Signal Process. 2019, 146, 87–100. [Google Scholar] [CrossRef]
  29. Zhou, H.; Chen, J.; Dong, G. Application of maximum correlated kurtosis deconvolution on rolling element bearing fault diagnosis. Eng. Asset Manag. Syst. 2015, 159–172. [Google Scholar] [CrossRef]
  30. Mcdonald, G.L.; Zhao, Q. Multipoint Optimal Minimum Entropy Deconvolution and Convolution Fix: Application to vibration fault detection. Mech. Syst. Signal Process. 2016, 82, 461–477. [Google Scholar] [CrossRef]
  31. Cai, W.; Wang, Z. Application of an improved multipoint optimal minimum entropy deconvolution adjusted for gearbox composite fault diagnosis. Sensors 2018, 18, 2861. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, Z.; Wang, J.; Zhao, Z.; Wang, R. A novel method for multi-fault feature extraction of a gearbox under strong background noise. Entropy 2018, 20, 10. [Google Scholar] [CrossRef] [Green Version]
  33. Ma, H.; Feng, Z. Planet bearing fault diagnosis using multipoint Optimal Minimum Entropy Deconvolution Adjusted. J. Sound Vib. 2019, 449, 235–273. [Google Scholar] [CrossRef]
  34. Li, Y.; Cheng, G.; Chen, X.; Pang, Y. Research on Bearing Fault Diagnosis Method Based on Filter Features of MOMLMEDA and LSTM. Entropy 2019, 21, 25. [Google Scholar] [CrossRef] [Green Version]
  35. Zhu, X.; Wang, Y. Fault diagnosis of rolling bearings based on the MOMEDA and Teager energy operator. J. Vib. Shock 2018, 37, 1–8. [Google Scholar]
  36. Deng, W. An effective improved co-evolution ant colony optimization algorithm with multi-strategies and its application. Int. J. Bio-Inspired Comput. 2020, 2020, 1–10. [Google Scholar]
  37. Li, J.; Jiang, J.; Fan, X.; Wang, H.; Song, L.; Liu, W.; Yang, J.; Chen, L. A New Method for Weak Fault Feature Extraction Based on Improved MED. Shock Vib. 2018, 2018, 1–11. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, C.; Li, H.; Ou, J.; Hu, R.; Hu, S.; Liu, A. Identification of planetary gearbox weak compound fault based on parallel dual-parameter optimized resonance sparse decomposition and improved MOMEDA. Measurement 2020, 165, 1–18. [Google Scholar] [CrossRef]
  39. Zhou, C.; Ma, J.; Wu, J.; Feng, Z. A Parameter Adaptive MOMEDA Method Based on Grasshopper Optimization Algorithm to Extract Fault Features. Math. Probl. Eng. 2019, 2019, 1–12. [Google Scholar] [CrossRef] [Green Version]
  40. Xu, X.; Zhao, M.; Lin, J.; Lei, Y. Envelope harmonic-to-noise ratio for periodic impulses detection and its application to bearing diagnosis. Measurement 2016, 91, 385–397. [Google Scholar] [CrossRef]
  41. Li, Z.; Ma, J.; Wang, X.; Wu, J. MVMD-MOMEDA-TKEO model and its application in feature extraction for rolling bearings. Entropy 2019, 21, 1–25. [Google Scholar]
  42. Borghesani, P.; Pennacchi, P.; Chatterton, S. The relationship between kurtosis- and envelope-based indexes for the diagnostic of rolling element bearings. Mech. Syst. Signal Process. 2014, 43, 25–43. [Google Scholar] [CrossRef]
  43. Jonathan, A.W.; Paul, K.R. Optimizations of the energy grid search algorithm in continuous-energy Monte Carlo particle transport codes. Comput. Phys. Commun. 2015, 196, 134–142. [Google Scholar]
  44. Maragos, P.; Kaiser, J.; Quatieri, T. On amplitude and frequency demodulation using energy operators. IEEE Trans. Signal Process. 1993, 41, 1532–1550. [Google Scholar] [CrossRef]
  45. Miao, Y.; Zhao, M.; Lin, J.; Lei, Y. Application of an improved maximum correlated kurtosis deconvolution method for fault diagnosis of rolling element bearings. Mech. Syst. Signal Process. 2017, 92, 173–195. [Google Scholar] [CrossRef]
  46. Loparo, K. Bearings Vibration Data Set, Case Western Reserve University. Available online: http://csegroups.case.edu/bearingdatacenter/pages/download-data-file (accessed on 20 August 2020).
  47. Lee, J.; Qiu, H.; Yu, G.; Lin, J. Rexnord Technical Services. Ims. University of Cincinnati. “Bearing Data Set”, NASA Ames Prognostics Data Repository. Moffett Field, CA: NASA Ames Research Center; 2007. Available online: https://ti.arc.nasa.gov/tech/dash/groups/pcoe/prognostic-data-repository/ (accessed on 14 September 2020).
  48. Wang, L.; Liu, Z.; Miao, Q.; Zhang, X. Time frequency analysis based on ensemble local mean decomposition and fast kurtogram for rotating machinery fault diagnosis. Mech. Syst. Signal Process. 2018, 103, 60–75. [Google Scholar] [CrossRef]
  49. Wang, J.; He, Q.; Kong, F. Multiscale envelope manifold for enhanced fault diagnosis of rotating machines. Mech. Syst. Signal Process. 2015, 52–53, 376–392. [Google Scholar] [CrossRef]
  50. Jiang, X.; Cheng, X.; Shi, J.; Huang, W.; Shen, C.; Zhu, Z. A new l0-norm embedded MED method for roller element bearing fault diagnosis at early stage of damage. Measurement 2018, 127, 414–424. [Google Scholar] [CrossRef]
Figure 1. Autocorrelation function spectrum.
Figure 1. Autocorrelation function spectrum.
Sensors 21 00533 g001
Figure 2. Implementation flow chart of IMOMEDA-TKEO method.
Figure 2. Implementation flow chart of IMOMEDA-TKEO method.
Sensors 21 00533 g002
Figure 3. Bearing test rig of CWRU data center (a) The simulation experiment platform (b) The rolling bearing.
Figure 3. Bearing test rig of CWRU data center (a) The simulation experiment platform (b) The rolling bearing.
Sensors 21 00533 g003
Figure 4. Vibration signal time and frequency analysis.
Figure 4. Vibration signal time and frequency analysis.
Sensors 21 00533 g004
Figure 5. Autocorrelation function spectrum.
Figure 5. Autocorrelation function spectrum.
Sensors 21 00533 g005
Figure 6. EHNR spectrum.
Figure 6. EHNR spectrum.
Sensors 21 00533 g006
Figure 7. Trend chart of EHNR spectral entropy.
Figure 7. Trend chart of EHNR spectral entropy.
Sensors 21 00533 g007
Figure 8. Signal deconvolution (IMOMEDA).
Figure 8. Signal deconvolution (IMOMEDA).
Sensors 21 00533 g008
Figure 9. Teager-Kaiser energy spectrum (x_cov(t)).
Figure 9. Teager-Kaiser energy spectrum (x_cov(t)).
Sensors 21 00533 g009
Figure 10. Multipoint kurtosis spectrum.
Figure 10. Multipoint kurtosis spectrum.
Sensors 21 00533 g010
Figure 11. Signal deconvolution (MOMEDA).
Figure 11. Signal deconvolution (MOMEDA).
Sensors 21 00533 g011
Figure 12. Teager-Kaiser energy spectrum (x_cov(t)).
Figure 12. Teager-Kaiser energy spectrum (x_cov(t)).
Sensors 21 00533 g012
Figure 13. Fast kurtogram (outer race).
Figure 13. Fast kurtogram (outer race).
Sensors 21 00533 g013
Figure 14. Fault characteristic spectrum (FK).
Figure 14. Fault characteristic spectrum (FK).
Sensors 21 00533 g014
Figure 15. Signal deconvolution (L = 2050).
Figure 15. Signal deconvolution (L = 2050).
Sensors 21 00533 g015
Figure 16. Teager-Kaiser energy spectrum (L = 2050).
Figure 16. Teager-Kaiser energy spectrum (L = 2050).
Sensors 21 00533 g016
Figure 17. Signal deconvolution (pseudo period T = 97.2).
Figure 17. Signal deconvolution (pseudo period T = 97.2).
Sensors 21 00533 g017
Figure 18. Teager-Kaiser energy spectrum (T = 97.2).
Figure 18. Teager-Kaiser energy spectrum (T = 97.2).
Sensors 21 00533 g018
Figure 19. Vibration signal time and frequency analysis.
Figure 19. Vibration signal time and frequency analysis.
Sensors 21 00533 g019
Figure 20. Teager-Kaiser energy spectrum (inner race) (a) IMOMEDA-TKEO, T = 75, L = 2100 (b) MOMEDA-TKEO, T = (59,88), L = 500 (c) MOMEDA-TKEO, T = (59,88), L = 2100 (d) MOMEDA-TKEO, pseudo fault period T = 86, L = 2100.
Figure 20. Teager-Kaiser energy spectrum (inner race) (a) IMOMEDA-TKEO, T = 75, L = 2100 (b) MOMEDA-TKEO, T = (59,88), L = 500 (c) MOMEDA-TKEO, T = (59,88), L = 2100 (d) MOMEDA-TKEO, pseudo fault period T = 86, L = 2100.
Sensors 21 00533 g020
Figure 21. Vibration signal analysis results based on FK (a) Fast kurtogram (inner race) (b) Fault characteristic spectrum (FK).
Figure 21. Vibration signal analysis results based on FK (a) Fast kurtogram (inner race) (b) Fault characteristic spectrum (FK).
Sensors 21 00533 g021
Figure 22. Bearing test rig of NASA data center (a) Test platform (b) Sensor layout.
Figure 22. Bearing test rig of NASA data center (a) Test platform (b) Sensor layout.
Sensors 21 00533 g022
Figure 23. Vibration signal time and frequency analysis (a) Time domain waveform (outer race) (b) Frequency domain waveform (outer race) (c) Time domain waveform (inner race) (d) Frequency domain waveform (inner race).
Figure 23. Vibration signal time and frequency analysis (a) Time domain waveform (outer race) (b) Frequency domain waveform (outer race) (c) Time domain waveform (inner race) (d) Frequency domain waveform (inner race).
Sensors 21 00533 g023
Figure 24. Teager-Kaiser energy spectrum (outer race) (a) IMOMEDA-TKEO, T = 86, L = 2100 (b) MOMEDA-TKEO, T = (68,102), L = 500. (c) MOMEDA-TKEO, T = (68,102), L = 2100 (d) MOMEDA-TKEO, pseudo fault period T = 99, L = 2100.
Figure 24. Teager-Kaiser energy spectrum (outer race) (a) IMOMEDA-TKEO, T = 86, L = 2100 (b) MOMEDA-TKEO, T = (68,102), L = 500. (c) MOMEDA-TKEO, T = (68,102), L = 2100 (d) MOMEDA-TKEO, pseudo fault period T = 99, L = 2100.
Sensors 21 00533 g024
Figure 25. Teager-Kaiser energy spectrum (inner race) (a) IMOMEDA-TKEO, T = 69, L = 2050 (b) MOMEDA-TKEO, T = (54,81), L = 500 (c) MOMEDA-TKEO, T = (54,81), L = 2050 (d) MOMEDA-TKEO, pseudo fault period T = 73, L = 2050.
Figure 25. Teager-Kaiser energy spectrum (inner race) (a) IMOMEDA-TKEO, T = 69, L = 2050 (b) MOMEDA-TKEO, T = (54,81), L = 500 (c) MOMEDA-TKEO, T = (54,81), L = 2050 (d) MOMEDA-TKEO, pseudo fault period T = 73, L = 2050.
Sensors 21 00533 g025
Figure 26. Vibration signal analysis results based on FK (a) Fault characteristic spectrum(outer race) (b) Fault characteristic spectrum(inner race).
Figure 26. Vibration signal analysis results based on FK (a) Fault characteristic spectrum(outer race) (b) Fault characteristic spectrum(inner race).
Sensors 21 00533 g026
Figure 27. Bearing testbed.
Figure 27. Bearing testbed.
Sensors 21 00533 g027
Figure 28. Fault bearings.
Figure 28. Fault bearings.
Sensors 21 00533 g028
Figure 29. Vibration signal time and frequency analysis (a) Time domain waveform (outer race) (b) Frequency domain waveform (outer race) (c) Time domain waveform (inner race) (d) Frequency domain waveform (inner race).
Figure 29. Vibration signal time and frequency analysis (a) Time domain waveform (outer race) (b) Frequency domain waveform (outer race) (c) Time domain waveform (inner race) (d) Frequency domain waveform (inner race).
Sensors 21 00533 g029
Figure 30. Teager-Kaiser energy spectrum (outer race) (a) IMOMEDA-TKEO, T = 472, L = 2200 (b) MOMEDA-TKEO, T = (382,573), L = 500 (c) MOMEDA-TKEO, T = (382,573), L = 2200 (d) MOMEDA-TKEO, pseudo fault period T = 512, L = 2200.
Figure 30. Teager-Kaiser energy spectrum (outer race) (a) IMOMEDA-TKEO, T = 472, L = 2200 (b) MOMEDA-TKEO, T = (382,573), L = 500 (c) MOMEDA-TKEO, T = (382,573), L = 2200 (d) MOMEDA-TKEO, pseudo fault period T = 512, L = 2200.
Sensors 21 00533 g030
Figure 31. Teager-Kaiser energy spectrum (inner race) (a) IMOMEDA-TKEO, T = 319, L = 2400 (b) MOMEDA-TKEO, T = (252,378), L = 500 (c) MOMEDA-TKEO, T = (252,378), L = 2400 (d) MOMEDA-TKEO, pseudo fault period T = 340, L = 2400.
Figure 31. Teager-Kaiser energy spectrum (inner race) (a) IMOMEDA-TKEO, T = 319, L = 2400 (b) MOMEDA-TKEO, T = (252,378), L = 500 (c) MOMEDA-TKEO, T = (252,378), L = 2400 (d) MOMEDA-TKEO, pseudo fault period T = 340, L = 2400.
Sensors 21 00533 g031
Figure 32. Vibration signal analysis results based on FK(a) Fault characteristic spectrum (outer race) (b) Fault characteristic spectrum (inner race).
Figure 32. Vibration signal analysis results based on FK(a) Fault characteristic spectrum (outer race) (b) Fault characteristic spectrum (inner race).
Sensors 21 00533 g032
Table 1. The brief overview of the comparative experiment.
Table 1. The brief overview of the comparative experiment.
TestAlgorithmParameters OptimizationData Sets
Test 1IMOMEDA-TKEOThe EHNR spectrum is utilized to calculate fault period T and the improved grid search method base on EHNR spectral entropy is proposed to determine optimal filter length L.(a) CWRU data sets
(b) NASA data sets
(c) KUST-SY data sets
Test 2MOMEDA-TKEO-class 1According to [38], set the default range of T is (0.8Tα,1.2Tα), Tα = fs/f,Where fs is sampling frequency and f is fault frequency. The default filter length L [30] is put to 500.
Test 3MOMEDA-TKEO-class 2The value range of T is consistent with Test 2 and the method to determine L follows Test 1.
Test 4Fast kurtogramSpectral kurtosis is taken as the function of Short Time Fourier Transform (STFT) window width to search the optimal filter parameters.
Table 2. Bearing basic parameters.
Table 2. Bearing basic parameters.
Rolling Bearing-6205-2RSJEMSKFRolling Element Number (Z)Inner Diameter
(Inches)
Outer Diameter
(Inches)
Rolling Element Diameter d (Inches)Limiting Speed (rpm)
90.98432.04720.312618,000
Contact angle (θ)Pitch circle diameter D (inches)Dynamic load rating (N)Static load rating (N)/
1.53714,8007800/
Table 3. The value of SNRfd with different methods.
Table 3. The value of SNRfd with different methods.
MethodIMOMEDA-TKEO
(T = 112; L = 2050)
MOMEDA-TKEO-Class 1
(T = (90,134); L = 500)
MOMEDA-TKEO-Class 2
(T = (90, 134); L = 2050)
Fast Kurtogram
Value of SNRfd2.66000.23522.50140.1784
Table 4. The value of SNRfd with different methods.
Table 4. The value of SNRfd with different methods.
MethodIMOMEDA-TKEO
(T = 75; L = 2100)
MOMEDA-TKEO-Class 1
(T = (59,88); L = 500)
MOMEDA-TKEO-Class 2
(T = (59, 88); L = 2100)
Fast Kurtogram
Value of SNRfd1.45380.76641.42000.3632
Table 5. Bearing foundation parameters.
Table 5. Bearing foundation parameters.
Rolling Element Number (Z)Contact Angle (θ)Rolling Element Diameter d (mm)Pitch Diameter D (mm)Rotational Speed (rpm)Dynamic Load Rating (N)Static Load Rating (N)Limiting Speed (rpm)
1615.170.3312.8152000650074702500
Table 6. The value of SNRfd with different methods.
Table 6. The value of SNRfd with different methods.
Outer raceMethodIMOMEDA-TKEO
(T = 86; L = 2100)
MOMEDA-TKEO-class 1
(T = (68,102); L = 500)
MOMEDA-TKEO-class 2 (T = (68,102); L = 2100)Fast kurtogram
Value of SNRfd1.93440.16751.81670.1139
Inner raceMethodIMOMEDA-TKEO
(T = 69;L = 2050)
MOMEDA-TKEO-class 1
(T = (54,81); L = 500)
MOMEDA-TKEO-class 2 (T = (54,81); L = 2050)Fast kurtogram
Value of SNRfd1.42210.11271.25270.1105
Table 7. Bearing foundation parameters.
Table 7. Bearing foundation parameters.
Rolling Element Number (Z)Contact Angle (θ)Bearing Pitch Diameter (mm)Roller Diameter (mm)Dynamic Load Rating (N)Static Load Rating (N)Limiting Speed (rpm)
9039814,800780018,000
Table 8. The value of SNRfd with different methods.
Table 8. The value of SNRfd with different methods.
Outer RaceMethodIMOMEDA-TKEO
(T = 472; L = 2200)
MOMEDA-TKEO-class 1 (T = (382,573); L = 500)MOMEDA-TKEO-class 2 (T = (382,573); L = 2200)Fast kurtogram
Value of SNRfd2.88430.52281.19360.2830
Inner RaceMethodIMOMEDA-TKEO
(T = 319; L = 2400)
MOMEDA-TKEO-class 1 (T = (252,378); L = 500)MOMEDA-TKEO-class 2 (T = (252,378); L = 2400)Fast kurtogram
Value of SNRfd2.94230.51791.66560.2961
Table 9. The parameter optimization time and iterations of IMOMEDA algorithm.
Table 9. The parameter optimization time and iterations of IMOMEDA algorithm.
Fault TypeAlgorithmFault Period T/TimeFault Period T/IterationsFilter Length L/TimeFilter Length L/Iterations
CWRU(outer)IMOMEDA0.465s857.114s76
CWRU(inner)0.751s957.950s77
NASA(outer)0.594s964.015s76
NASA(inner)0.690s957.617s77
KUST-SY(outer)0.744s681.343s61
KUST-SY(inner)0.756s692.087s67
Table 10. Main program execution time of the comparative methods.
Table 10. Main program execution time of the comparative methods.
DatasetsCWRU
(Outer)
CWRU
(Inner)
NASA
(Outer)
NASA
(Inner)
KUST-SY
(Outer)
KUST-SY
(Inner)
Algorithm
IMOMEDA65.385s63.792s68.859s61.770s83.616s 102.741s
(T = 112)(T = 75)(T = 86)(T = 68)(T = 472)(T = 319)
(L = 2050)(L = 2100)(L = 2100)(L = 2050)(L = 2200)(L = 2400)
MOMEDA-TKEO-class 11.127s 1.334s 1.181s1.171s10.179s 7.389s
(T = (90,134))(T = (59,88))(T = (68,102))(T = (54,81))(T = (382,573))(T = (252,378))
(L = 500)(L = 500)(L = 500)(L = 500)(L = 500)(L = 500)
MOMEDA-TKEO-class 24.221s4.357s 4.685s4.562s21.490s 22.958s
(T = (90,134))(T = (59,88))(T = (68,102))(T = (54,81))(T = (382,573))(T = (252,378))
(L = 2050)(L = 2100)(L = 2100)(L = 2050)(L = 2400)(L = 2000)
FK0.809s0.763s0.825s0.837s9.543s8.906s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Ma, J.; Wang, X.; Li, X. An Optimal Parameter Selection Method for MOMEDA Based on EHNR and Its Spectral Entropy. Sensors 2021, 21, 533. https://doi.org/10.3390/s21020533

AMA Style

Li Z, Ma J, Wang X, Li X. An Optimal Parameter Selection Method for MOMEDA Based on EHNR and Its Spectral Entropy. Sensors. 2021; 21(2):533. https://doi.org/10.3390/s21020533

Chicago/Turabian Style

Li, Zhuorui, Jun Ma, Xiaodong Wang, and Xiang Li. 2021. "An Optimal Parameter Selection Method for MOMEDA Based on EHNR and Its Spectral Entropy" Sensors 21, no. 2: 533. https://doi.org/10.3390/s21020533

APA Style

Li, Z., Ma, J., Wang, X., & Li, X. (2021). An Optimal Parameter Selection Method for MOMEDA Based on EHNR and Its Spectral Entropy. Sensors, 21(2), 533. https://doi.org/10.3390/s21020533

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