Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Stereo Video Quality Metric Based on Multi-Dimensional Analysis
Next Article in Special Issue
A Machine Learning Approach for Gearbox System Fault Diagnosis
Previous Article in Journal
A Bit Shift Image Encryption Algorithm Based on Double Chaotic Systems
Previous Article in Special Issue
Investigation of Cyber-Security and Cyber-Crimes in Oil and Gas Sectors Using the Innovative Structures of Complex Intuitionistic Fuzzy Relations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intelligent Fault Diagnosis of Rolling-Element Bearings Using a Self-Adaptive Hierarchical Multiscale Fuzzy Entropy

1
School of Mechatronics Engineering, Nanjing Forestry University, Nanjing 210037, China
2
School of Mechanical Engineering, Southeast University, Nanjing 211189, China
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(9), 1128; https://doi.org/10.3390/e23091128
Submission received: 6 August 2021 / Revised: 22 August 2021 / Accepted: 23 August 2021 / Published: 30 August 2021

Abstract

:
The fuzzy-entropy-based complexity metric approach has achieved fruitful results in bearing fault diagnosis. However, traditional hierarchical fuzzy entropy (HFE) and multiscale fuzzy entropy (MFE) only excavate bearing fault information on different levels or scales, but do not consider bearing fault information on both multiple layers and multiple scales at the same time, thus easily resulting in incomplete fault information extraction and low-rise identification accuracy. Besides, the key parameters of most existing entropy-based complexity metric methods are selected based on specialist experience, which indicates that they lack self-adaptation. To address these problems, this paper proposes a new intelligent bearing fault diagnosis method based on self-adaptive hierarchical multiscale fuzzy entropy. On the one hand, by integrating the merits of HFE and MFE, a novel complexity metric method, named hierarchical multiscale fuzzy entropy (HMFE), is presented to extract a multidimensional feature matrix of the original bearing vibration signal, where the important parameters of HMFE are automatically determined by using the bird swarm algorithm (BSA). On the other hand, a nonlinear feature matrix classifier with strong robustness, known as support matrix machine (SMM), is introduced for learning the discriminant fault information directly from the extracted multidimensional feature matrix and automatically identifying different bearing health conditions. Two experimental results on bearing fault diagnosis show that the proposed method can obtain average identification accuracies of 99.92% and 99.83%, respectively, which are higher those of several representative entropies reported by this paper. Moreover, in the two experiments, the standard deviations of identification accuracy of the proposed method were, respectively, 0.1687 and 0.2705, which are also greater than those of the comparison methods mentioned in this paper. The effectiveness and superiority of the proposed method are verified by the experimental results.

1. Introduction

The rolling-element bearing is one of the essential parts of rotating machinery; its health condition directly affects the safe and steady operation of mechanical equipment [1]. When bearing failure occurs, it can easily cause economic losses for enterprises and even catastrophic accidents. Besides, the properties of the fault-bearing vibration signal are usually nonlinear and nonstationary, which indicates that it is difficult to obtain helpful bearing fault information using traditional methods [2,3,4]. Therefore, exploring a new and effective fault feature extraction method to ensure the normal running of rolling-element bearings is of great significance.
Currently, many signal complexity metric methods have been proposed for describing the complexity and irregularity of bearing vibration signals and extracting the corresponding bearing fault feature information, including correlation dimension [5], Lempel–Ziv complexity [6], pattern spectrum entropy [7], approximate entropy [8], sample entropy (SE) [9], permutation entropy (PE) [10], symbol dynamic entropy (SDE) [11], dispersion entropy (DE) [12], fuzzy entropy (FE) [13], and so on. These methods have received great attention in the field of bearing fault diagnosis and have obtained many achievements. However, they only excavate bearing fault information on a single scale. Concerning this problem, based on single-scale signal complexity metric methods, many multiscale analysis methods (e.g., multiscale approximate entropy [14], multiscale sample entropy (MSE) [15], multiscale permutation entropy (MPE) [16,17], multiscale symbol dynamic entropy (MSDE) [18], multiscale dispersion entropy (MDE) [19], and multiscale fuzzy entropy (MFE) [20]) are presented for bearing fault diagnosis. For instance, Han et al. [21] adopted hierarchical Lempel–Ziv complexity to extract the fault information hidden in both low-frequency and high-frequency components and achieve intelligent fault diagnosis of rotating machinery. Yan et al. [22] used multiscale pattern gradient spectrum entropy (MPGSE) to extract bearing fault information, then extreme learning machine (ELM) to classify different bearing fault patterns. Gao et al. [23] combined the improved local mean decomposition (LMD), multiscale permutation entropy (MPE), and hidden Markov model (HMM) for identifying different fault types of rolling element bearings. Wu et al. [24] utilized MPE to obtain bearing fault features, then used support vector machine (SVM) to finish the fault identification of a rolling bearing. Gan et al. [25] presented a composite multiscale fluctuation dispersion entropy (CMFDE) to extract multiscale fault information of rolling bearings and realized a classification process for bearing fault types. Zheng et al. [26] proposed a generalized composite multiscale permutation entropy (GCMPE) to extract bearing fault information, and then used SVM to identify bearing fault types. Wang et al. [27] adopted the generalized refined composite multiscale sample entropy (GRCMSE) to extract bearing fault feature information, and then, the optimized support vector machine was used to finish bearing fault diagnosis. Tang et al. [28] proposed a hierarchical instantaneous energy density dispersion entropy (HIEDDE) to extract fault features, then the extracted features were regarded as the input of dynamic time warping (DTW) to automatically identify the fault types for a gearbox. Li et al. [29] proposed hierarchical fuzzy entropy (HFE) to extract bearing fault features and adopted an improved support-vector-machine-based binary tree to identify bearing fault patterns. Meanwhile, Li et al. [30] combined the local mean decomposition (LMD) and improved multiscale fuzzy entropy (IMFE) to complete intelligent fault diagnosis of a rolling bearing. However, there are two shortcomings in the above-reported methods. On the one hand, the above-mentioned methods can only excavate bearing fault feature information on the different scales or different levels of the original vibration signal, which indicates that multiscale fault features of different levels or different frequency bands of the original vibration signal are not taken into account simultaneously. In other words, the fault feature information obtained by the above-mentioned methods is not very comprehensive and rich. On the other hand, the above-mentioned methods rely on human experience to select the important parameters of entropy, so they are not adaptive, and this can easily affect the performance of fault feature extraction. Hence, to solve these issues, this paper proposes a novel signal complexity metric method named hierarchical multiscale fuzzy entropy (HMFE). Specifically, the proposed HMFE method consists of three modules (i.e., hierarchical decomposition, improved coarse-grained process, and fuzzy entropy calculation). Firstly, the hierarchical decomposition process of the collected original signal is conducted to obtain a series of hierarchical components. Subsequently, the improved coarse-grained process of each hierarchical component is carried out to effectively obtain the improved coarse-grained time series of each hierarchical component at different scales. Finally, the fuzzy entropy of each improved coarse-grained series is calculated to obtain the final hierarchical multiscale fuzzy entropy (it is equivalent to a feature matrix at different levels and scales). Besides, an effective intelligent optimizer named bird swarm algorithm (BSA) is employed to self-adaptively determine the parameters of HMFE. In summary, the contributions and novelties of this paper are summarized as follows:
(1)
A novel signal complexity metric method named hierarchical multiscale fuzzy entropy (HMFE) is proposed via integrating hierarchic analysis into multiscale fuzzy entropy, which can effectively obtain more comprehensive and richer bearing fault feature information.
(2)
The important parameters of HMFE are selected automatically by using the bird swarm algorithm (BSA) method, which can effectively avoid the disadvantages of manual selecting of the important parameters of the existing entropy.
(3)
An intelligent bearing fault diagnosis method based on HMFE is presented, which can improve the identification accuracy of different health conditions of rolling-element bearings.
(4)
Two experimental cases and comparative analysis are conducted to verify the effectiveness and superiority of the proposed method.
The organization of the rest of this paper is as follows. Section 2 introduces the self-adaptive hierarchical multiscale fuzzy entropy, including the basic theory, adaptive parameter selection, and simulation analysis of HMFE. In Section 3, the concept of the support matrix machine is summarized. Section 4 shows the flowchart of the proposed method for bearing fault diagnosis. Section 5 validates the effectiveness and superiority of the proposed method through two experimental examples and contrastive analysis. The conclusions are drawn in Section 6.

2. Self-Adaptive Hierarchical Multiscale Fuzzy Entropy

In this section, self-adaptive hierarchical multiscale fuzzy entropy, which is mainly composed of HMFE and its adaptive parameter selection, is introduced.

2.1. HMFE

In this section, based on the advantages of HFE and MFE, a novel signal complexity metric method named hierarchical multi-scale fuzzy entropy (HMFE) is proposed. Figure 1 shows the flowchart of the proposed HMFE method. For a given time series  { x ( i ) , i = 1 , 2 , , N } , the calculation process of HMFE is described in detail as follows:
(1)
Firstly, the average and difference operator are respectively defined as follows:
Q 0 ( x ) = x ( 2 i ) + x ( 2 i + 1 ) 2 , i = 0 , 1 , 2 , , 2 n 1
Q 1 ( x ) = x ( 2 i ) x ( 2 i + 1 ) 2 , i = 0 , 1 , 2 , , 2 n 1
where n denotes the positive integer, 2 n 1 represents the length of the two operators, and Q 0 ( x ) and Q 1 ( x ) represent the low-frequency and high-frequency components of the original time series in the first layer decomposition, respectively.
(2)
Then, to implement the hierarchical analysis of a time series, when j = 0 or 1, the matrix form of the k-th layer operator Q j k is expressed as:
Q j k = 1 2 ( 1 ) j 2 0 0 0 0 0 0 1 2 ( 1 ) j 2 0 0 0 0 0 0 1 2 ( 1 ) j 2 2 n 1 × 2 n
(3)
To obtain the hierarchical components X k , e  of each layer in the process of hierarchical decomposition, here we define a one-dimensional vector as γ 1 , γ 2 , , γ k and an integral value as  e = p = 1 k 2 k p γ p , where { γ p , p = 1 , 2 , , k } { 0 , 1 } represents the average or difference operator at the p-th layer. Accordingly, the hierarchical component of the e-th node in the k-th layer can be expressed as:
X k , e = Q γ k k · Q γ k 1 k 1 · · Q γ 1 1 · x
where x represents the given time series.
(4)
Next, the improved coarse-grained time series y k , e ( τ ) = y k , e , 1 ( τ ) , y k , e , 2 ( τ ) , of each hierarchical component at the τ scale factor can be calculated by
y k , e , j ( τ ) = 1 τ i = j j + τ 1 X k , e i , 1 j N τ + 1
where N represents the length of the given time series x.
(5)
According to the definition of fuzzy entropy, calculate the fuzzy entropy of each improved coarse-grained series y k , e ( τ ) , so the final hierarchical multiscale fuzzy entropy can be obtained by the following:
H M F E ( x , m , k , e , τ , r ) = F E ( y k , e ( τ ) , m , τ , r )
where F E ( · ) represents the fuzzy entropy operation, m is the embedded dimension, k is the decomposition level, e represents the hierarchical node, τ denotes the scale factor, r = 0.15 × S D is the similarity tolerance controlling the width of membership function, and SD is the standard deviation of the original time series. In Figure 1, τ m  represents the predefined largest scale factor.

2.2. Adaptive Parameter Selection of HMFE

The previous studies have shown that fault feature extraction performance of the existing entropies (e.g., MSE, MPE, MFE, and HFE) is greatly affected by their parameter settings. As with traditional entropies (e.g., MSE, MPE, MFE, and HFE), the parameters of HMFE have a great impact on its performance, which indicates that it is necessary to adopt an effective method to achieve adaptive selection of the parameters of HMFE. At present, the intelligent swarm optimization algorithm has been proven to be effective in solving parameter selection problems, including particle swarm optimization (PSO), ant colony algorithm (ACO), bat algorithm (BA), and so on. Bird swarm algorithm (BSA) is a new bionic optimization algorithm proposed by Meng et al. [31] in 2016. Compared with other optimization algorithms, BSA has the advantages of high accuracy, strong stability, and fast convergence in the parameter optimization. Therefore, in this paper, BSA is adopted to automatically select several important parameters (i.e., the embedding dimension m, the decomposition level k and the scale factor τ ) of HMFE. Figure 2 plots the flowchart of parameter optimization of HMFE. The specific process of parameter optimization of HMFE based on BSA method is summarized as follows:
(1)
Initialize the population and set the BSA parameters. When the number of iterations t = 0, set the bird swarm size to N = 30 and the maximum iteration number to M = 50, initialize the flight frequency FQ, foraging frequency P, and several constants (i.e., C, S, FL, a1 and a2).
(2)
Calculate and compare the fitness value. According to the fitness function shown in Equation (7), the fitness value of bird swarm is calculated and compared to determine the optimal position of the individual and whole bird swarm.
fitness ( i ) = 1 x i x c + x i
where xi is the number of misclassified samples, xc is the number of samples correctly classified, and fitness(i) is the current fitness value of the i-th bird. When fitness(i) achieves the maximum value, pi,j is the corresponding optimal position of the individual bird swarm, and gj is the corresponding optimal position of the whole bird swarm.
(3)
The iterative operation is performed repeatedly, and the position update formula is determined by judging whether the operation t% × FQ has a remainder. Specific rules are summarized as follows:
If there is a remainder for t% × FQ, a uniformly distributed number is randomly generated. When the random number is less than the foraging frequency P, the foraging behavior is performed, and Equation (8) is used to update the position. Otherwise, the sentinel behavior is performed, and Equation (9) is used to update the position.
x i , j t + 1 = x i , j t + ( p i , j x i , j t ) × C × rand ( 0 , 1 ) + ( g j x i , j t ) × S × rand ( 0 , 1 )
where rand(0,1) is a uniformly distributed random number between 0 and 1, pi,j is the current best position of the i-th bird, gj is the current best position of bird swarm, and C and S are two positive numbers, which are called the cognitive acceleration coefficient and social acceleration coefficient, respectively.
x i , j t + 1 = x i , j t + A 1 ( mean j x i , j t ) × rand ( 0 , 1 ) + A 2 ( p k , j x i , j t ) × rand ( 1 , 1 ) A 1 = a 1 × exp ( pfit i sumfit + ε × N ) A 2 = a 2 × exp ( pfit i pfit k pfit k pfit i + ε N × pfit k sumfit + ε )
where x i , j t  is the position of individual birds in the t-th iteration, x i , j t + 1  is the position of individual birds in the t + 1 iteration, a 1 and a 2 are the positive number between 0 and 2, pfit i is the optimal fitness value of the i-th bird, sumfit is the sum of the optimal fitness values of bird swarm, ε is the smallest constant that avoids dividing the denominator by zero, k ( k i ) is a positive integer between 1 and N, and mean j  is the j-th element of the average position of the entire bird swarm.
If there is no remainder for t% × FQ when individual birds are the producers, Equation (10) is used to update the position. When individual birds are beggars, Equation (11) is used to update the position.
x i , j t + 1 = x i , j t + randn ( 0 , 1 ) × x i , j t
x i , j t + 1 = x i , j t + ( x k , j t x i , j t ) × F L × randn ( 0 , 1 )
where x i , j t  is the position of individual birds in the t-th iteration; x i , j t + 1  is the position of individual birds in the t + 1 iteration; randn(0,1) represents the random number of Gaussian distribution with mean value 0 and standard deviation 1, k [ 1 , 2 , 3 , , N ] , k i ; and FL denotes the integer between 0 and 2.
(4)
Update the position of each bird swarm according to the rules in step (3). If the individual of the current bird swarm is better than the individual of the previous bird swarm, the current individual bird swarm is regarded as the optimal position. Otherwise, the previous individual bird swarm is retained as the optimal position to continue the update of bird swarm.
(5)
Judge whether the stop condition is met. If the maximum number of iterations or the minimum error rate is reached, the whole optimization process will be stopped, and the optimal position of bird swarm (i.e., the optimal combination parameters of HMFE) will be outputted. Otherwise, the iteration process will continue to be conducted until the stop condition is satisfied.

2.3. Comparison Analysis Using Simulation Signal

To investigate the performance of the proposed HMFE method, according to the literature [32], one bearing vibration signal x O R ( t ) (i.e., the simulation signal 1) containing outer race fault and one bearing vibration signal x I R ( t ) (i.e., the simulation signal 2) containing inner race fault are established by using Equations (12) and (13), respectively.
x O R ( t ) = f 1 ( t ) * h 1 ( t ) + n ( t ) f 1 ( t ) = k = 0 N 1 δ ( t k T 1 ) h 1 ( t ) = e 600 π t cos ( 2 π f n t )
x I R ( t ) = f 2 ( t ) * h 2 ( t ) + n ( t ) f 2 ( t ) = k = 0 N 1 cos ( 2 π f r t ) δ ( t k T 2 ) h 2 ( t ) = e 800 π t cos ( 2 π f ^ n t )
where the asterisk * indicates the convolution operation; f 1 ( t ) and f 2 ( t ) are, respectively, the cyclic impact signal caused by bearing outer race fault and inner race fault; n(t) is the white Gaussian noise with a mean of 0 and standard deviation of 1; δ ( ) indicates the Dirac delta function, k = 0 , 1 , , N 1 ; N indicates the total number of impact impulses; and T 1 and T 2 are the time intervals between two adjacent impulses under bearing outer race fault and inner race fault, respectively. Bearing outer race fault frequency and inner race fault frequency are f O R = 1 / T 1 = 50 Hz and f I R = 1 / T 2 = 90 Hz, respectively. f r = 25 Hz is the rotation frequency of the driving shaft, h 1 ( t ) and h 2 ( t ) are, respectively, the impulse response signal caused by the bearing outer race fault and that caused by the inner race fault; and f n = 3000 Hz and f ^ n = 5000 Hz are the natural frequencies of the excitation system under bearing outer race fault and inner race fault, respectively. The sampling frequency and data length of simulation signal are set as 16,384 Hz and 4096 points, respectively.
Figure 3 shows the time domain waveforms and their corresponding amplitude spectra for two simulation signals. As can be seen from Figure 3, the periodic impact impulses of two simulation signals are drowned in the random noise. That is, it is difficult to identify bearing fault types by the direct observation of time domain waveform, which means that an effective method is needed to extract the fault features of the two simulation signals for distinguishing and recognizing them. Hence, the Euclidean distance (ED) of three entropies (i.e., HMFE, HFE, and MFE) of two simulation signals are calculated to compare their feature extraction ability. In the three methods (i.e., HMFE, HFE, and MFE), the embedded dimension is m = 3 and the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the simulation signal. Besides, for HMFE and HFE, the decomposition level is k = 3. For HMFE and MFE, the scale factor is τ = 8. Figure 4a–d shows the FE value obtained by the three methods (i.e., HMFE, HFE, and MFE) for two simulation signals. Apparently, the FE value at different levels and scales can be extracted by using HMFE, whereas HFE and MFE can only obtain the FE value of different levels or scales. This indicates that HMFE can obtain more comprehensive feature information compared to HFE and MFE. Table 1 lists the Euclidean distance (ED) and calculation time of different methods. As shown in Table 1, the ED of the proposed HMFE is largest, which shows that the feature information obtained by HMFE is more differentiated compared with HFE and MFE. Besides, the calculation time of HMFE is significantly less than that of MFE, but it is greater than that of HFE. This is mainly due to the hierarchical decomposition and multiscale analysis process being integrated in the HMFE method, which increases the computational burden of HMFE to some extent, but it is generally acceptable for practical bearing vibration signal analysis.

3. Support Matrix Machine

The intelligent fault identification step is required after the entropy-based fault feature extraction. At present, many linear or nonlinear classification models have been proposed for intelligent bearing fault identification, including linear discriminant analysis (LDA), BP neural network (BPNN), K-nearest neighbor (KNN), extreme learning machine (ELM) [33], partial least squares (PLS) [34], and support vector machine (SVM). However, these methods are only applicable to the classification of multidimensional feature vectors. When the above-mentioned classification model is used for the processing of the multidimensional feature matrix, their performance will be reduced. Hence, to solve this problem, a new nonlinear classification model named support matrix machine (SMM) was recently proposed by Luo et al. [35] in 2015, which can automatically learn the useful discriminant information from the multidimensional feature matrix. In view of this, in this paper, our plan is to adopt SMM to process the multidimensional feature matrix constructed by HMFE and achieve intelligent fault identification of the rolling bearing. The basic theory of SMM is described as follows:
Suppose that { X i , y i } i = 1 n is one given training set, where X i R d 1 × d 2  is the i-th input matrix, y i { 1 , 1 } is the training label, and d1 and d2 represent the number of rows and columns of the input matrix, respectively. Simply speaking, SMM implements the model training and classification process through the hinge loss function and spectral elastic network penalty function, as shown below.
arg min W , b 1 2 tr ( W T W ) + λ W * + C i = 1 n ξ i s . t . y i [ tr ( W T X i ) + b ] 1 ξ i , i = 1 , 2 , , n
where tr ( W T W ) = W F 2 is the squared F-norm [36], W R d 1 × d 2 is the regression coefficient matrix, W * is the kernel norm, λ and C are respectively kernel parameter and penalty parameter, and b is the bias term. Spectral elastic network penalty function consists of the squared F-norm and kernel norm. The squared F-norm can be used to control the complexity of model and prevent the over-fitting phenomenon. The kernel norm can approximately replace the rank of the regression coefficient matrix. Therefore, based on the classification property of spectral elastic networks, SMM can effectively capture intrinsic feature information and correlation in the input matrix. The alternating direction multiplier method (ADMM) can be used to optimize SMM. In particular, after introducing auxiliary variables Z R d 1 × d 2 , Equation (14) can be rewritten as:
arg min W , b , Z P ( W , b ) + Q ( Z ) s . t . Z W = 0
where P ( W , b ) and Q ( Z )  can be expressed as follows:
P ( W , b ) = 1 2 t r ( W T W ) + C i = 1 n { 1 y i [ t r ( W T X i ) + b ] } + Q ( Z ) = τ Z *
Specifically, Equation (14) can be solved by the augmented Lagrange multiplier method [37], that is:
L ( Z , W , b , M ) = P ( W , b ) + Q ( Z ) + t r ( M T ( Z W ) ) + β 2 ( Z W ) F 2
where M R d 1 × d 2  is the Lagrange multiplier and β  is a positive hyperparameter. The updated equations of Z, W, B and the Lagrange multiplier M can be expressed as:
Z t + 1 = arg min Z L ( Z , W t , b t , M t ) ( W t + 1 , b t + 1 ) = arg min ( W , b ) L ( Z t + 1 , W , b , M t ) M t + 1 = M t β ( W t + 1 Z t + 1 )
where t is the number of iterations. According to [35], the optimal solution of Equation (18) can be written as:
W ^ = 1 β + 1 ( M + β Z + i = 1 n α ^ i y i X i ) b ^ = 1 n i = 1 n { y i tr ( W ^ T X i ) } Z ^ = 1 β D τ ( β W M )
Finally, the decision function of SMM can be expressed as:
f ( X ) = sgn ( tr ( W T X ) + b )
where sgn ( · ) indicates the sign function and tr ( · ) indicates the trace of a matrix.

4. Proposed Method

To obtain abundant fault information and improve bearing fault identification accuracy, a new approach based on self-adaptive hierarchical multiscale fuzzy entropy is proposed for intelligent fault diagnosis of rolling-element bearings. The proposed method consists of three steps (i.e., bearing vibration data collection, hierarchical multiscale feature extraction, and bearing health condition identification). Figure 5 shows the flowchart of the proposed method, and its main steps are described below:
Step1:Bearing vibration data collection. Bearing vibration data under different health conditions are obtained by installing the accelerometer on the experimental equipment.
Step2:Hierarchical multiscale feature extraction. Using the HMFE method, the FE value of bearing vibration signals at different levels and scales are calculated, where the important parameters of the HMFE method are automatically determined by using the bird swarm algorithm method. Meanwhile, the calculated HMFE of each bearing health condition is adopted to construct a multidimensional feature matrix.
Step3:Bearing health condition identification. The constructed multidimensional feature matrix is randomly divided into the training sample matrix and testing sample matrix, where the training sample matrix is adopted to train the SMM classification model, and the testing sample matrix is fed into the well-trained SMM classification model to automatically identify different bearing health conditions.

5. Experimental Verification

In this section, two experiments about bearing fault diagnosis are conducted to show the effectiveness of the proposed method. Furthermore, comparisons among the proposed method and several representative entropies are adopted to highlight the superiority of the proposed method. Finally, we discuss the future research direction.

5.1. Case 1: Bearing Benchmark Data from CWRU

Bearing benchmark data from Case Western Reserve University (CWRU) is applied to validate the effectiveness of the proposed method. Figure 6 shows the experimental platform and its schematic diagram, which are mainly composed of induction motor, testing bearing, torque transducer, and load motor. Table 2 lists the specifications of the test bearing. In the experiment process, the sampling frequency is set as 12 kHz and the motor speed is set as 1797 rpm. Besides, nine single-point faults (i.e., inner race slight fault (IRSF), inner race medium fault (IRMF), inner race heavy fault (IRHF), outer race slight fault (ORSF), outer race medium fault (ORMF), outer race heavy fault (ORHF), ball slight fault (BSF), ball medium fault (BMF), and ball heavy fault (BHF)) are manufactured on normal bearings by using electric discharging machining technique [38]. One accelerometer is mounted on the bearing block of the drive-end of the induction motor to collect, respectively, bearing vibration data containing the ten health conditions (i.e., normal, IRSF, IRMF, IRHF, ORSF, ORMF, ORHF, BSF, BMF, and BHF). For each bearing health conditions, 50 data samples with a length of 2048 points are intercepted by using a non-overlapping sliding window. That is, a total of 500 data samples are obtained for the whole health conditions, where 25 data samples of each bearing health condition are randomly selected as the training set and the remaining 25 data samples are treated as the testing set. Table 3 presents a detailed description of bearing vibration data under ten health conditions. Figure 7 shows the time domain waveform of bearing vibration data under different health conditions in case 1. Seen from Figure 7, due to the measured bearing vibration data containing some background noises and having nonlinear non-stationary characteristics, it is difficult to accurately distinguish different health states of bearings by directly observing the time domain waveforms.
In order to effectively extract bearing fault features and improve the identification accuracy, the proposed method is adopted to analyze the experimental data listed in Table 3. Firstly, the BSA method is used to automatically determine the optimal parameters of HMFE at the embedded dimension m = 3, the decomposition level k = 3, and the scale factor τ = 8. Subsequently, the HMFE method with the optimal parameters is conducted to calculate the FE value of different bearing health conditions at different levels and scales. Meanwhile, the calculated HMFE of each bearing health condition is adopted to build the multidimensional feature matrix with a size of 8 × 8 × 500. Figure 8 shows the calculation results of HMFE of bearing vibration data under different health conditions. It can be seen from Figure 8 that the FE value of different bearing health conditions is different at some levels or scales due to the proper integration of hierarchical decomposition and multiscale coarse-grained analysis, which helps for the subsequent identification of bearing health conditions. Finally, the constructed multidimensional feature matrix is randomly divided into the training sample matrix with a size of 8 × 8 × 250 and the testing sample matrix with a size of 8 × 8 × 250. Besides, the SMM classification model is trained by the training sample matrix, and the testing sample matrix is fed into the well-trained SMM model to identify different bearing health conditions and automatically report the fault diagnosis results. Figure 9 shows the identification results of the proposed method in the first trial. Seen from Figure 9, the proposed method can achieve identification accuracy of 100%, which means that the proposed method is effective in identifying different fault categories and severities of rolling-element bearings.
To show the effectiveness of the parameter optimization process used in the proposed HMFE, the fault identification accuracy of the proposed method containing different parameters of HMFE is calculated, and the results are shown in Table 4. Seen from Table 4, when the BSA method is used to select the optimal parameters (i.e., the embedded dimension m = 3, the decomposition level k = 3, and the scale factor τ = 8) of HMFE, the proposed method with the parameter optimization process can obtain a high identification accuracy. This means that the BSA-parameter optimization process in the proposed HMFE is very useful for bearing fault identification. Besides, as shown in Table 4, the identification accuracy is greater than 95% only when the scale factor τ is higher than or equal to 7. For the phenomenon of the lower accuracy at low scale factor τ, here an explanation is given. Concretely, in the HMFE method, when the embedded dimension m and decomposition level k stay constant, if the scale factor τ is set as smaller, the feature matrix with a smaller dimension will be extracted, which indicates that the feature information obtained by HMFE is relatively less and the identification accuracy of the proposed method is reduced. Theoretically, the larger scale factor τ has more feature information and higher identification accuracy. However, the bigger the scale factor τ is, the longer the calculation time of the proposed method is, and the extracted features will contain some redundant information. Hence, regarding the scale factor τ, bigger is not necessarily better, and it needs to be selected appropriately.
Due to the number of training and testing samples having a great influence on the identification performance of the proposed method, the identification results of the proposed method are further investigated at different percentages of training samples. Specifically, identification results of the proposed method are calculated with the number of training samples set at, respectively, 50, 100, 150, 200, 250, 300, 350, 400, and 450. Identification accuracy and training time of the proposed method are shown in Figure 10 for when the percentage of training samples in the whole dataset is, respectively, 10% (50/500), 20% (100/500), 30% (150/500), 40% (200/500), 50% (250/500), 60% (300/500), 70% (350/500), 80% (400/500), and 90% (450/500). As can be seen from Figure 10, when the percentage of training samples is equal to or greater than 0.5 (i.e., 50%), the proposed method achieves an identification accuracy of 100%. However, as the percentage of training samples increases, the training time of the proposed method has an upward trend. In other words, the higher the number of training samples, the better trained the SMM model is, but the training time of the proposed method will be longer. Therefore, to strike a balance between identification accuracy and training time, the number of training and testing samples is set to the same percentages used in this paper.
To further show the effectiveness of the SMM classification model used in the proposed method, the proposed HMFE is combined with four classification models (i.e., SMM, SVM, ELM, and BPNN) to identify bearing fault patterns. The parameters of these classification models are set based on the previous work [39]. Specifically, in the SMM and SVM classifier, the kernel parameter λ =1/n and the penalty parameter C = 1, where n represents the dimensions of the extracted feature matrix. In the ELM classifier, the activation function uses a Sigmoid function, and the number of hidden nodes is N = 20. In the BPNN classifier, the number of hidden nodes is N = 20, the maximum training number is I = 500, and the learning rate is σ = 0.1. Each method is conducted for 5 trials to compare the results objectively. Table 5 lists the identification results of combining HMFE and different classification models (i.e., SMM, SVM, ELM, and BPNN). It can be clearly seen from Table 5 that the average recognition accuracy of the proposed method (i.e., the combination of HMFE and SMM) is significantly higher than the other combination methods, which indicates that the validity of using SMM classification model in the proposed method is verified.
To show the superiority of the proposed method, comparisons are made between the proposed method and six existing representative entropies (i.e., GCMPE [26], GRCMSE [27], HFE [29], MFE [30], refined composite multiscale dispersion entropy (RCMDE) [40], and hierarchical sample entropy (HSE) [41]). To avoid the randomness and occasionality of the recognition results of different methods, 10 trials of each method are conducted to observe and compare the identification results. Besides, to ensure a fair comparison, the important parameters of all methods are selected by the BSA method, and the classification process is completed by the SMM model. Table 6 lists the setting of the optimal parameters of different methods. Figure 11 shows the identification results of different methods in 10 trials. Besides, Table 7 gives the final fault diagnosis results of the different methods for 10 trials, including the maximum accuracy, minimum accuracy, average accuracy, and standard deviation. As shown in Figure 11 and Table 7, the average identification accuracy of the proposed method is 99.92%, which is higher than that of the six comparative methods (i.e., HFE, MFE, RCMDE, GCMPE, GRCMSE, and HSE), which are 97.08%, 95.44%, 97.72%, 96.20%, 94.28%, and 91.92%, respectively. This indicates that the identification ability of the proposed method for bearing health condition is better than the other methods. Moreover, the standard deviation of the proposed method is 0.1687, which is smaller than that of the six comparative methods (i.e., HFE, MFE, RCMDE, GCMPE, GRCMSE and HSE), which are 0.2700, 0.3373, 0.2700, 0.2828, 0.3293, and 0.2530, respectively. That is, compared with the six comparative methods (i.e., HFE, MFE, RCMDE, GCMPE, GRCMSE, and HSE), the proposed method has a better stability in identifying bearing fault category and severity. Hence, the effectiveness and advantages of the proposed method have been verified in the comparative analysis of case 1.

5.2. Case 2: Bearing Vibration Data from Laboratory

In this section, bearing vibration data collected from an ABLT-1A experimental device is used to verify the effectiveness of the proposed method in identifying bearing fault patterns. Figure 12a,b show, respectively, the bearing accelerated life test bench and its corresponding structure diagram, which is mainly composed of motor, drive belt, bearing test module, and coupling. The bearing test module is installed with four bearings. Type and size of testing bearings are listed in Table 8. Five bearing health conditions (i.e., outer race fault (ORF), inner race fault (IRF), ball fault (BF), outer and inner race compound fault (OIRCF), outer race and ball compound fault (ORBCF)) on the normal bearing 1 are generated by using wire electrical discharge machining. Figure 13 shows the photos of different faulty bearings. In this experiment, the sampling frequency is set as 10,240 Hz, and motor speed is stable at 1050 rpm. In order to simulate the weak bearing fault signal brought by the long transmission path, one accelerometer is mounted on a position away from faulty bearing 1 to collect bearing vibration data of different health conditions. For each bearing health condition, 60 data samples with a length of 2048 points are obtained via the sliding window, where 30 data samples of each bearing health condition are randomly selected as the training set, and the other 30 data samples are regarded as the testing set. That is to say, the training set and the testing set each have 180 samples. Table 9 provides the detailed description of the experimental dataset. Figure 14 plots time domain waveforms of bearing vibration data under different health conditions in case 2. As shown in Figure 14, due to the environmental noise and transmission path interference, it is difficult to identify different bearing fault patterns by directly observing the time domain waveforms of bearing vibration signals, which means that it is urgent to adopt an effective method to identify different health conditions of rolling bearings.
Firstly, to verify the effectiveness of the proposed method, the proposed method is utilized to analyze the experimental data listed in Table 9. According to the flowchart of the proposed method, the important parameters of HMFE are selected adaptively by using the BSA method. Concretely, in the HMFE method, the embedded dimension m = 4, the decomposition level k = 3, and the scale factor τ = 8 are used. Then, for each bearing health condition, the proposed HMFE method with the optimized parameters is conducted to obtain a multidimensional feature matrix with a size of 8 × 8 × 360. Figure 15 plots the calculation results of HMFE of one data sample of different bearing health conditions. Obviously, as seen from Figure 15, a comprehensive and recognizable feature matrix can be obtained by calculating HMFE. Finally, according to the sample percentage of 1:1, the obtained multidimensional feature matrix is randomly and averagely divided into the training sample matrix and testing sample matrix. The training sample matrix is adopted for the training of the SMM model, and the testing sample matrix is entered into the well-trained SMM model for identifying different bearing health conditions. Figure 16 gives the identification results of the proposed method in the first trial. It can be seen from Figure 16 that six bearing health conditions (i.e., Normal, ORF, IRF, BF, OIRCF, ORBCF) are all correctly identified by using the proposed method, which implies that the efficacy of the proposed method in bearing fault identification is preliminarily proven.
As with in case 1, to illustrate the validity and necessity of the parameter optimization process of the proposed method, we analyzed the identification results of the proposed method under different parameter settings, as listed in Table 10. Seen from Table 10, when the raw bearing vibration data are analyzed by the HMFE method with the optimal parameters (i.e., the embedded dimension m = 4, the decomposition level k = 3, and the scale factor τ = 8), the highest identification accuracy (100%) can be obtained, which means that the parameter selection of HMFE has a great influence on its diagnosis performance. Meanwhile, this also indicates that it is very useful for bearing fault identification by using BAS to select the important parameters of HMFE.
Similarly, to show the influence of the number of training samples on the proposed method, the identification results of the proposed method at different percentage of training samples are considered. The identification accuracy and training time of the proposed method are calculated, and the calculation results are plotted in Figure 17, for when the percentage of the training samples compared to all samples is, respectively, 10% (36/360), 20% (72/360), 30% (108/360), 40% (144/360), 50% (180/360), 60% (216/360), 70% (252/360), 80% (288/360), and 90% (324/360). Seen from Figure 17, the identification accuracy and training time of the proposed method are increased with the increase of the percentage of the training samples. Therefore, in this paper, the percentage of training samples is set as 50% (i.e., the number of training and testing samples is the same) to achieve a compromise between the identification accuracy and training time, which is relatively reasonable.
To further verify the effectiveness of the SMM model used in the proposed method, the identification results of combining HMFE and different classifiers (i.e., SMM, SVM, ELM, and BPNN) are calculated, and the results are listed in Table 11. Note that the model parameters of different classifiers are set the same as case 1. It can be seen from Table 11 that the proposed method (i.e., the combination of HMFE and SMM) can achieve a higher identification accuracy compared with other combination methods. This further validates the effectiveness of using the SMM model in the proposed method.
As in case 1, to further prove superiority of the proposed method, the proposed method and six entropies (i.e., HFE, MFE, RCMDE, GCMPE, GRCMSE, and HSE) are used to analyze the same experimental data. Again, 10 trials are performed to avoid the contingency of the diagnosis results, and the BSA method is used to determine the important parameters of all methods to ensure the fairness of the comparison. Table 12 gives the parameter settings of different methods. Figure 18 plots the identification accuracy of different methods over 10 trials. Moreover, Table 13 gives the specific comparison results of the identification accuracy of the different methods. Seen from Figure 18 and Table 13, compared with the other six representative methods (i.e., HFE, MFE, RCMDE, GCMPE, GRCMSE, and HSE), the proposed method can obtain the highest average identification accuracy (i.e., 99.83%), which further demonstrates the superiority of the proposed method for bearing health condition identification. Besides, it can be also found from Table 13 that the proposed method can obtain a smaller standard deviation (0.2705) than the other six representative methods (i.e., HFE, MFE, RCMDE, GCMPE, GRCMSE, and HSE), which implies that the proposed method has a more stable fault identification performance than the other six representative methods. Therefore, by the experiment and the comparative analysis, the superiority and efficacy of the proposed method in identifying bearing fault patterns are highlighted.

5.3. Further Discussion

Although the effectiveness and superiority of the proposed method in bearing health condition identification is verified by the above two experiments, some limitations and future works related to the proposed method are still left to discuss. These future works and limitations are summarized as follows:
(1)
In the proposed method, the key parameters of HMFE are determined by the bird swarm algorithm (BSA) method, which is helpful for bearing fault feature extraction. We all know that some other advanced optimizer (e.g., Grey wolf optimization (GWO), Whale optimization algorithm (WOA), Grasshopper optimization algorithm (GOA)) can also be introduced to automatically determine the key parameters of HMFE. Therefore, in our future work, the parameter selection problem of HMFE will continue to be studied by adopting other advanced optimizers instead of the BSA method.
(2)
In the final step of the proposed method, although the support matrix machine (SMM) is employed to achieve the automatic identification of fault patterns of rolling-element bearing and obtain a good diagnosis result, there are many other advanced classification models in the previously reported literature, including the improved versions (e.g., the non-parallel least squares support matrix machine [42], nonlinear kernel support matrix machine) of SMM and deep learning models (e.g., convolutional neural network [43], deep regularized variational autoencoder [44], deep belief network [45], and other deep learning methods [46]). Hence, in our future work, SMM of the proposed method will be replaced by these advanced classification models to automatically obtain the identification results of bearing fault patterns.
(3)
The proposed method is proven to be effective for identifying the bearing health condition at constant speed, but it is unknown for bearing fault identification under variable speed. Hence, the proposed method will be extended to solve the problem of bearing fault identification under variable speed, which is regarded as our future focus. Besides this, in our future work, other faults (e.g., gear, rotor, and blade) of rotating machinery will also be diagnosed by applying the proposed method.
(4)
In this paper, the proposed HMFE is only applied for single-channel sensor data analysis of rolling-element bearings. For bearing vibration data analysis of multi-channel sensors, with the help of the idea of multichannel data processing used in the existing multivariate multiscale entropy, a new multichannel data processing method named multivariate hierarchical multiscale fuzzy entropy (MHMFE) will be designed to solve the problem of multivariate fault diagnosis in our future work.
(5)
Due to the addition of the parameter optimization algorithm, the biggest limitation of the proposed method lies in the large calculation time. Therefore, to solve this issue and improve the computational efficiency of the proposed method, in our future work, some sensitive indicators (e.g., Chebychev distance and Mahalanobis distance) can be used to instead of the complex optimizer to automatically select the key parameters of HMFE. Besides, in the running of the algorithm, graphic processing units (GPU) can be adopted instead of the central processing unit (CPU) to accelerate the calculation process of the proposed method.

6. Conclusions

This paper proposes a new intelligent bearing fault diagnosis method based on self-adaptive hierarchical multiscale fuzzy entropy, which can not only solve the disadvantages of artificial selection of important parameters of most existing entropies, but also obtain richer and more comprehensive bearing fault feature information. Two experimental cases validate the effectiveness of the proposed method in bearing fault identification. Furthermore, compared with some reported entropy methods, the superiority of the proposed method in bearing fault identification is verified. Some specific conclusions and contributions are summarized as follows:
(1)
A new signal complexity metric method named hierarchical multiscale fuzzy entropy is developed by integrating the hierarchical decomposition into multiscale fuzzy entropy, which is aimed at improving fault feature extraction performance.
(2)
An effective parameter optimizer called bird swarm algorithm is introduced to automatically choose several important parameters of hierarchical multiscale fuzzy entropy, which can avoid the dependence of parameter selection of the existing entropy on specialist experience.
(3)
The effectiveness of the proposed method in the identification of bearing fault types and fault severity is verified by experimental and contrastive analysis. The experimental results show that compared with some existing representative multiscale entropies or hierarchical entropies, the proposed method can achieve broader and richer fault feature information and its identification accuracy has been greatly improved, which indicates that the proposed method has a certain competitiveness in bearing health condition identification. This study provides a new perspective for intelligent fault diagnosis for rolling-element bearings.

Author Contributions

Data curation, conceptualization, methodology, software, visualization, writing—original draft, writing—review and editing, funding acquisition: X.Y.; Software, formal analysis: Y.X.; Project administration; resources: M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 52005265), in part by the Natural Science Fund for Colleges and Universities in Jiangsu Province (No. 20KJB460002), in part by Scientific Research Foundation of Nanjing Forestry University (No. 163040095 and 163040117), in part by Jiangsu Provincial Key Research and Development Program (No. BE2019030637), and by the Jiangsu Provincial High-level Innovation and Entrepreneurship Talent Introduction Plan.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study are all owned by the research group and will not be transmitted.

Acknowledgments

The authors would like to thank the anonymous reviewers and the editor for their valuable and insightful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yan, X.; Liu, Y.; Jia, M. Multiscale cascading deep belief network for fault identification of rotating machinery under various working conditions. Knowl. Based Syst. 2020, 193, 105484. [Google Scholar] [CrossRef]
  2. Wang, X.; Si, S.; Wei, Y.; Li, Y. The optimized multi-scale permutation entropy and its application in compound fault diagnosis of rotating machinery. Entropy 2019, 21, 170. [Google Scholar] [CrossRef] [Green Version]
  3. Ye, M.; Yan, X.; Jia, M. Rolling bearing fault diagnosis based on VMD-MPE and PSO-SVM. Entropy 2021, 23, 762. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, F.; Sun, W.; Wang, H.; Xu, T. Fault diagnosis of a wind turbine gearbox based on improved variational mode algorithm and information entropy. Entropy 2021, 23, 794. [Google Scholar] [CrossRef] [PubMed]
  5. Chen, C.; Tian, H.; Wu, D.; Pan, Q.; Wang, H.; Liu, X. A fault diagnosis method for satellite flywheel bearings based on 3D correlation dimension clustering technology. IEEE Access 2018, 6, 78483–78492. [Google Scholar] [CrossRef]
  6. Hong, H.; Liang, M. Fault severity assessment for rolling element bearings using the Lempel–Ziv complexity and continuous wavelet transform. J. Sound Vib. 2009, 320, 452–468. [Google Scholar] [CrossRef]
  7. Hao, R.; Peng, Z.; Feng, Z.; Chu, F. Application of support vector machine based on pattern spectrum entropy in fault diagnostics of rolling element bearings. Meas. Sci. Technol. 2011, 22, 45708. [Google Scholar] [CrossRef]
  8. Yan, R.; Gao, R.X. Approximate entropy as a diagnostic tool for machine health monitoring. Mech. Syst. Signal Process. 2007, 21, 824–839. [Google Scholar] [CrossRef]
  9. Han, M.; Pan, J. A fault diagnosis method combined with LMD, sample entropy and energy ratio for roller bearings. Measurement 2015, 76, 7–19. [Google Scholar] [CrossRef]
  10. Bandt, C.; Pompe, B. Permutation entropy: A natural complexity measure for time series. Phys. Rev. Lett. 2002, 88, 174102. [Google Scholar] [CrossRef] [PubMed]
  11. Yan, X.; Liu, Y.; Jia, M. A fault diagnosis approach for rolling bearing integrated SGMD, IMSDE and multiclass relevance vector machine. Sensors 2020, 20, 4352. [Google Scholar] [CrossRef]
  12. Azami, H.; Escudero, J. Amplitude-and fluctuation-based dispersion entropy. Entropy 2018, 20, 210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Deng, W.; Zhang, S.; Zhao, H.; Yang, X. A novel fault diagnosis method based on integrating empirical wavelet transform and fuzzy entropy for motor bearing. IEEE Access 2018, 6, 35042–35056. [Google Scholar] [CrossRef]
  14. Ma, J.; Li, Z.; Li, C.; Zhan, L.; Zhang, G.Z. Rolling bearing fault diagnosis based on refined composite multi-scale approximate entropy and optimized probabilistic neural network. Entropy 2021, 23, 259. [Google Scholar] [CrossRef] [PubMed]
  15. Wu, S.D.; Wu, C.W.; Lin, S.G.; Wang, C.C.; Lee, K.Y. Time series analysis using composite multiscale entropy. Entropy 2013, 15, 1069–1084. [Google Scholar] [CrossRef] [Green Version]
  16. Tang, G.; Wang, X.; He, Y. A novel method of fault diagnosis for rolling bearing based on dual tree complex wavelet packet transform and improved multiscale permutation entropy. Math. Probl. Eng. 2016, 2016, 5432648. [Google Scholar] [CrossRef] [Green Version]
  17. Zheng, J.; Cheng, J.; Yang, Y. Multiscale permutation entropy based rolling bearing fault diagnosis. Shock Vib. 2014, 2014, 1–8. [Google Scholar] [CrossRef]
  18. Han, M.; Wu, Y.; Wang, Y.; Liu, W. Roller bearing fault diagnosis based on LMD and multi-scale symbolic dynamic information entropy. J. Mech. Sci. Technol. 2021, 35, 1993–2005. [Google Scholar] [CrossRef]
  19. Yan, X.; Jia, M. Intelligent fault diagnosis of rotating machinery using improved multiscale dispersion entropy and mRMR feature selection. Knowl. Based Syst. 2019, 163, 450–471. [Google Scholar] [CrossRef]
  20. Li, Y.; Miao, B.; Zhang, W.; Chen, P.; Liu, J.; Jiang, X. Refined composite multiscale fuzzy entropy: Localized defect detection of rolling element bearing. J. Mech. Sci. Technol. 2019, 33, 109–120. [Google Scholar] [CrossRef]
  21. Han, B.; Wang, S.; Zhu, Q.; Yang, X.; Li, Y. Intelligent fault diagnosis of rotating machinery using hierarchical Lempel-Ziv complexity. Appl. Sci. 2020, 10, 4221. [Google Scholar] [CrossRef]
  22. Yan, X.; Liu, Y.; Ding, P.; Jia, M. Fault diagnosis of rolling-element bearing using multiscale pattern gradient spectrum entropy coupled with Laplacian score. Complexity 2020, 2020, 1–29. [Google Scholar] [CrossRef] [Green Version]
  23. Gao, Y.; Villecco, F.; Li, M.; Song, W. Multi-scale permutation entropy based on improved LMD and HMM for rolling bearing diagnosis. Entropy 2017, 19, 176. [Google Scholar] [CrossRef] [Green Version]
  24. Wu, S.D.; Wu, P.H.; Wu, C.W.; Ding, J.J.; Wang, C.C. Bearing fault diagnosis based on multiscale permutation entropy and support vector machine. Entropy 2012, 14, 1343–1356. [Google Scholar] [CrossRef] [Green Version]
  25. Gan, X.; Lu, H.; Yang, G. Fault diagnosis method for rolling bearings based on composite multiscale fluctuation dispersion entropy. Entropy 2019, 21, 290. [Google Scholar] [CrossRef] [Green Version]
  26. Zheng, J.; Pan, H.; Yang, S.; Cheng, J. Generalized composite multiscale permutation entropy and Laplacian score based rolling bearing fault diagnosis. Mech. Syst. Signal Process. 2018, 99, 229–243. [Google Scholar] [CrossRef]
  27. Wang, Z.; Yao, L.; Cai, Y. Rolling bearing fault diagnosis using generalized refined composite multiscale sample entropy and optimized support vector machine. Measurement 2020, 156, 107574. [Google Scholar] [CrossRef]
  28. Tang, G.; Pang, B.; He, Y.; Tian, T. Gearbox fault diagnosis based on hierarchical instantaneous energy density dispersion entropy and dynamic time warping. Entropy 2019, 21, 593. [Google Scholar] [CrossRef] [Green Version]
  29. Li, Y.; Xu, M.; Zhao, H.; Huang, W. Hierarchical fuzzy entropy and improved support vector machine based binary tree approach for rolling bearing fault diagnosis. Mech. Mach. Theory 2016, 98, 114–132. [Google Scholar] [CrossRef]
  30. Li, Y.; Xu, M.; Wang, R.; Huang, W. A fault diagnosis scheme for rolling bearing based on local mean decomposition and improved multiscale fuzzy entropy. J. Sound Vib. 2016, 360, 277–299. [Google Scholar] [CrossRef]
  31. Meng, X.; Gao, X.; Lu, L.; Liu, Y.; Zhang, H. A new bio-inspired optimisation algorithm: Bird swarm algorithm. J. Exp. Theor. Artif. In. 2016, 28, 1–15. [Google Scholar] [CrossRef]
  32. Yan, X.; Jia, M.; Xiang, L. Compound fault diagnosis of rotating machinery based on OVMD and a 1.5-dimension envelope spectrum. Meas. Sci. Technol. 2016, 27, 075002. [Google Scholar] [CrossRef]
  33. Haidong, S.; Hongkai, J.; Xingqiu, L.; Shuaipeng, W. Intelligent fault diagnosis of rolling bearing using deep wavelet auto-encoder with extreme learning machine. Knowl. Based Syst. 2018, 140, 1–14. [Google Scholar] [CrossRef]
  34. Yan, X.; Liu, Y.; Xu, Y.; Jia, M. Multichannel fault diagnosis of wind turbine driving system using multivariate singular spectrum decomposition and improved Kolmogorov complexity. Renew. Energ. 2021, 170, 724–748. [Google Scholar] [CrossRef]
  35. Luo, L.; Xie, Y.; Zhang, Z.; Li, W.J. Support matrix machines. In Proceedings of the 32nd International Conference on Machine Learning, Lille, France, 6–11 July 2015. [Google Scholar]
  36. Lancaster, P.; Tismenetsky, M. The Theory of Matrices with Application; Academic Press: Cambridge, MA, USA, 1985. [Google Scholar]
  37. Tan, H.; Cheng, B.; Feng, J.; Feng, G.; Wang, W.; Zhang, Y.J. Low-n-rank tensor recovery based on multi-linear augmented Lagrange multiplier method. Neurocomputing 2013, 119, 144–152. [Google Scholar] [CrossRef]
  38. Ho, K.; Newman, S. State of the art electrical discharge machining (EDM). Int. J. Mach. Tool. Manuf. 2003, 43, 1287–1300. [Google Scholar] [CrossRef]
  39. Yan, X.; Jia, M. A novel optimized SVM classification algorithm with multi-domain feature and its application to fault diagnosis of rolling bearing. Neurocomputing 2018, 313, 47–64. [Google Scholar] [CrossRef]
  40. Luo, S.; Yang, W.; Luo, Y. Fault Diagnosis of a rolling bearing based on adaptive sparest narrow-band decomposition and refined composite multiscale dispersion entropy. Entropy 2020, 22, 375. [Google Scholar]
  41. Zhu, K.; Song, X.; Xue, D. A roller bearing fault diagnosis method based on hierarchical entropy and support vector machine with particle swarm optimization algorithm. Measurement 2014, 47, 669–675. [Google Scholar] [CrossRef]
  42. Li, X.; Yang, Y.; Pan, H.; Cheng, J. Non-parallel least squares support matrix machine for rolling bearing fault diagnosis. Mech. Mach. Theory 2020, 145, 103676. [Google Scholar] [CrossRef]
  43. Xie, C.; Liu, Y.; Zeng, W.; Lu, X. An improved method for single image super-resolution based on deep learning. Signal Image Video Process. 2019, 13, 557–565. [Google Scholar] [CrossRef]
  44. Yan, X.; She, D.; Xu, Y.; Jia, M. Deep regularized variational autoencoder for intelligent fault diagnosis of rotor-bearing system within entire life-cycle process. Knowl. Based Syst. 2021, 226, 107142. [Google Scholar] [CrossRef]
  45. Liu, F.; Liu, B.; Sun, C.; Liu, M.; Wang, X. Deep belief network-based approaches for link prediction in signed social networks. Entropy 2015, 17, 2140–2169. [Google Scholar] [CrossRef] [Green Version]
  46. Yu, Y.; Liu, Y.; Chen, J.; Jiang, D.; Zhuang, Z.; Wu, X. Detection method for bolted connection looseness at small angles of timber structures based on deep learning. Sensors 2021, 21, 3106. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the HMFE method.
Figure 1. Flowchart of the HMFE method.
Entropy 23 01128 g001
Figure 2. Flowchart of parameter optimization of HMFE.
Figure 2. Flowchart of parameter optimization of HMFE.
Entropy 23 01128 g002
Figure 3. Time domain waveforms and their corresponding amplitude spectra for two simulation signals.
Figure 3. Time domain waveforms and their corresponding amplitude spectra for two simulation signals.
Entropy 23 01128 g003
Figure 4. FE value obtained by different methods for two simulation signals: (a) HMFE of simulation signal 1, (b) HMFE of simulation signal 2, (c) HFE of simulation signal 1, (d) HFE of simulation signal 2, (e) MFE of simulation signal 1, and (f) MFE of simulation signal 2.
Figure 4. FE value obtained by different methods for two simulation signals: (a) HMFE of simulation signal 1, (b) HMFE of simulation signal 2, (c) HFE of simulation signal 1, (d) HFE of simulation signal 2, (e) MFE of simulation signal 1, and (f) MFE of simulation signal 2.
Entropy 23 01128 g004
Figure 5. Flowchart of the proposed method for bearing fault identification.
Figure 5. Flowchart of the proposed method for bearing fault identification.
Entropy 23 01128 g005
Figure 6. (a) The experimental platform and (b) its schematic diagram.
Figure 6. (a) The experimental platform and (b) its schematic diagram.
Entropy 23 01128 g006
Figure 7. Time domain waveforms of bearing vibration data under different health conditions in case 1.
Figure 7. Time domain waveforms of bearing vibration data under different health conditions in case 1.
Entropy 23 01128 g007
Figure 8. HMFE of bearing vibration data under different health conditions in case 1.
Figure 8. HMFE of bearing vibration data under different health conditions in case 1.
Entropy 23 01128 g008
Figure 9. Classification results of the proposed method in case 1.
Figure 9. Classification results of the proposed method in case 1.
Entropy 23 01128 g009
Figure 10. The identification results of the proposed method under different percentages of training samples in case 1.
Figure 10. The identification results of the proposed method under different percentages of training samples in case 1.
Entropy 23 01128 g010
Figure 11. The identification results of 10 trials of different methods in case 1.
Figure 11. The identification results of 10 trials of different methods in case 1.
Entropy 23 01128 g011
Figure 12. (a) Bearing accelerated life test bench and (b) its corresponding schematic structure drawing.
Figure 12. (a) Bearing accelerated life test bench and (b) its corresponding schematic structure drawing.
Entropy 23 01128 g012
Figure 13. Photo of the faulty bearing: (a) ORF, (b) IRF, (c) OIRCF, and (d) BF.
Figure 13. Photo of the faulty bearing: (a) ORF, (b) IRF, (c) OIRCF, and (d) BF.
Entropy 23 01128 g013
Figure 14. Time domain waveforms of bearing vibration data under different health conditions in case 2.
Figure 14. Time domain waveforms of bearing vibration data under different health conditions in case 2.
Entropy 23 01128 g014
Figure 15. HMFE of bearing vibration data under different health conditions in case 2.
Figure 15. HMFE of bearing vibration data under different health conditions in case 2.
Entropy 23 01128 g015
Figure 16. Classification results of the proposed method in case 2.
Figure 16. Classification results of the proposed method in case 2.
Entropy 23 01128 g016
Figure 17. The identification results of the proposed method under different percentages of training samples in case 2.
Figure 17. The identification results of the proposed method under different percentages of training samples in case 2.
Entropy 23 01128 g017
Figure 18. The identification results of 10 trials of different methods in case 2.
Figure 18. The identification results of 10 trials of different methods in case 2.
Entropy 23 01128 g018
Table 1. Performance comparison among different methods.
Table 1. Performance comparison among different methods.
MethodsEuclidean DistanceCalculation Time (s)
HMFE0.34215.167
HFE0.20482.163
MFE0.06139.768
Table 2. The specification of test bearing.
Table 2. The specification of test bearing.
Bearing TypeRoller Diameter (mm)Pitch Diameter (mm)Number of the RollerContact Angle (°)
SKF6205-2RS7.9439.0490
Table 3. The detailed description of experimental dataset in case 1.
Table 3. The detailed description of experimental dataset in case 1.
Bearing Health ConditionsAbbreviationFault Size (inches)Number of
Training Samples
Number of
Testing Samples
Class Label
NormalNORM025251
Inner race slight faultIRSF0.00725252
Inner race medium faultIRMF0.01425253
Inner race heavy faultIRHF0.02125254
Outer race slight faultORSF0.00725255
Outer race medium faultORMF0.01425256
Outer race heavy faultORHF0.02125257
Ball slight faultBSF0.00725258
Ball medium faultBMF0.01425259
Ball heavy faultBHF0.021252510
Table 4. The identification results of the proposed method with different parameters in case 1.
Table 4. The identification results of the proposed method with different parameters in case 1.
Parameter Setting of HMFEIdentification
Accuracy (%)
The Embedded
Dimension m
The Decomposition Level kThe Scale
Factor τ
12487.60
22693.20
23796.80
33798.80
338100
43899.20
441298.40
541698.80
Table 5. The identification results of combining HMFE and different classifiers in case 1.
Table 5. The identification results of combining HMFE and different classifiers in case 1.
ClassifierIdentification Accuracy Obtained by Combining HMFE and Different Classifier Methods in 5 TrialsAverage
Accuracy (%)
12345
SMM10099.6010010010099.92
SVM97.6098.0097.6097.6097.2097.60
ELM98.4098.8098.4098.0098.4098.40
BPNN96.4096.8096.4097.2096.8096.72
Table 6. Parameter settings for different methods in case 1.
Table 6. Parameter settings for different methods in case 1.
MethodsThe Optimal Parameter Setting
HMFEThe embedded dimension m = 3, the decomposition level k = 3, and the scale factor τ = 8, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
HFEThe embedded dimension m = 3, the decomposition level k = 3, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
MFEThe embedded dimension m = 3, the scale factor τ = 10, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
RCMDEThe embedded dimension m = 3, time delay d = 1, the number of classes c = 5, the scale factor τ = 10.
GCMPEThe embedded dimension m = 3, time delay d = 1, the scale factor τ = 12.
GRCMSEThe embedded dimension m = 3, the scale factor τ = 10, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
HSEThe embedded dimension m = 3, the decomposition level k = 3, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
Table 7. Comparison results of different methods in case 1.
Table 7. Comparison results of different methods in case 1.
MethodsMaximum Accuracy (%) Minimum Accuracy (%)Average Accuracy (%) Standard Deviation
HMFE10099.6099.920.1687
HFE97.6096.8097.080.2700
MFE96.0095.2095.440.3373
RCMDE98.0097.2097.720.2700
GCMPE96.8096.0096.200.2828
GRCMSE94.8094.0094.280.3293
HSE92.4091.6091.920.2530
Table 8. Parameters of rolling-element bearing.
Table 8. Parameters of rolling-element bearing.
Bearing TypeBall Diameter (mm)Pitch Diameter (mm)Number of BallsContact Angle (°)
HRB62057.9439.0490
Table 9. Detailed description of the experimental dataset in case 2.
Table 9. Detailed description of the experimental dataset in case 2.
Bearing Health ConditionsAbbreviationNumber of
Training Samples
Number of
Testing Samples
Class
Labels
NormalNORM30301
Outer race faultORF30302
Inner race faultIRF30303
Ball faultBF30304
Outer and inner race compound faultOIRCF30305
Outer race and ball compound faultORBCF30306
Table 10. The identification results of the proposed method with different parameters in case 2.
Table 10. The identification results of the proposed method with different parameters in case 2.
Parameter Setting of HMFEIdentification
Accuracy (%)
The Embedded
Dimension m
The Decomposition
Level k
The Scale Factor τ
12489.44
22693.33
23797.22
33798.33
33899.44
438100
441298.33
541698.89
Table 11. The identification results of combining HMFE and different classifiers in case 2.
Table 11. The identification results of combining HMFE and different classifiers in case 2.
ClassifierIdentification Accuracy Obtained by Combining HMFE and Different Classifier Methods in 5 TrialsAverage
Accuracy (%)
12345
SMM10010099.4410099.4499.77
SVM97.2296.6797.2296.6796.1196.78
ELM98.3397.7898.3397.2298.3397.99
BPNN95.5596.1195.5596.1195.0095.66
Table 12. Parameter settings for different methods in case 2.
Table 12. Parameter settings for different methods in case 2.
MethodsThe Optimal Parameter Setting
HMFEThe embedded dimension m = 4, the decomposition level k = 3, and the scale factor τ = 8, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
HFEThe embedded dimension m = 4, the decomposition level k = 3, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
MFEThe embedded dimension m = 3, the scale factor τ = 12, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
RCMDEThe embedded dimension m = 3, time delay d = 1, the number of classes c = 6, the scale factor τ = 12.
GCMPEThe embedded dimension m = 3, time delay d = 1, the scale factor τ = 15.
GRCMSEThe embedded dimension m = 4, the scale factor τ = 10, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
HSEThe embedded dimension m = 3, the decomposition level k = 3, the similarity tolerance r = 0.15 × S D , where SD is the standard deviation of the original signal.
Table 13. Comparison results of different methods in case 2.
Table 13. Comparison results of different methods in case 2.
MethodsMaximum
Accuracy (%)
Minimum
Accuracy (%)
Average
Accuracy (%)
Standard
Deviation
HMFE10099.4499.830.2705
HFE96.6795.5596.270.3780
MFE94.4493.3394.160.3913
RCMDE96.6795.5596.160.4132
GCMPE95.5594.4494.830.3758
GRCMSE92.7791.6792.220.3667
HSE90.5589.4090.040.4226
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yan, X.; Xu, Y.; Jia, M. Intelligent Fault Diagnosis of Rolling-Element Bearings Using a Self-Adaptive Hierarchical Multiscale Fuzzy Entropy. Entropy 2021, 23, 1128. https://doi.org/10.3390/e23091128

AMA Style

Yan X, Xu Y, Jia M. Intelligent Fault Diagnosis of Rolling-Element Bearings Using a Self-Adaptive Hierarchical Multiscale Fuzzy Entropy. Entropy. 2021; 23(9):1128. https://doi.org/10.3390/e23091128

Chicago/Turabian Style

Yan, Xiaoan, Yadong Xu, and Minping Jia. 2021. "Intelligent Fault Diagnosis of Rolling-Element Bearings Using a Self-Adaptive Hierarchical Multiscale Fuzzy Entropy" Entropy 23, no. 9: 1128. https://doi.org/10.3390/e23091128

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