Next Article in Journal
A DFT Study of Phosphate Ion Adsorption on Graphene Nanodots: Implications for Sensing
Next Article in Special Issue
Intelligent Diagnosis of Rolling Bearings Fault Based on Multisignal Fusion and MTF-ResNet
Previous Article in Journal
Synthesis of ZnO Nanoflower Arrays on a Protrusion Sapphire Substrate and Application of Al-Decorated ZnO Nanoflower Matrix in Gas Sensors
Previous Article in Special Issue
Deconvolved Fractional Fourier Domain Beamforming for Linear Frequency Modulation Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dual-Optimization Fault Diagnosis Method for Rolling Bearings Based on Hierarchical Slope Entropy and SVM Synergized with Shark Optimization Algorithm

1
School of Automation and Information Engineering, Xi’an University of Technology, Xi’an 710048, China
2
Shaanxi Key Laboratory of Complex System Control and Intelligent Information Processing, Xi’an University of Technology, Xi’an 710048, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(12), 5630; https://doi.org/10.3390/s23125630
Submission received: 2 June 2023 / Revised: 11 June 2023 / Accepted: 13 June 2023 / Published: 16 June 2023
(This article belongs to the Special Issue Detection and Feature Extraction in Acoustic Sensor Signals)

Abstract

:
Slope entropy (SlopEn) has been widely applied in fault diagnosis and has exhibited excellent performance, while SlopEn suffers from the problem of threshold selection. Aiming to further enhance the identifying capability of SlopEn in fault diagnosis, on the basis of SlopEn, the concept of hierarchy is introduced, and a new complexity feature, namely hierarchical slope entropy (HSlopEn), is proposed. Meanwhile, to address the problems of the threshold selection of HSlopEn and a support vector machine (SVM), the white shark optimizer (WSO) is applied to optimize both HSlopEn and an SVM, and WSO-HSlopEn and WSO-SVM are proposed, respectively. Then, a dual-optimization fault diagnosis method for rolling bearings based on WSO-HSlopEn and WSO-SVM is put forward. We conducted measured experiments on single- and multi-feature scenarios, and the experimental results demonstrated that whether single-feature or multi-feature, the WSO-HSlopEn and WSO-SVM fault diagnosis method has the highest recognition rate compared to other hierarchical entropies; moreover, under multi-features, the recognition rates are all higher than 97.5%, and the more features we select, the better the recognition effect. When five nodes are selected, the highest recognition rate reaches 100%.

1. Introduction

Rolling bearings, as a key component in rotating machinery, serve a very significant role in modern industry. However, because of the increasingly sophisticated and complex structure of bearings and their common use in harsh working environments, rolling bearings are very prone to failures, which can lead to economic losses and even endanger personal safety [1,2,3]. Therefore, aiming to ensure the normal work of rotating machinery and reduce maintenance costs, it is of great importance to carry out fault diagnoses of rolling bearings [4,5,6].
Since bearing vibration signals contain rich state information about the bearing during operation, a vibration analysis method is broadly applied to rolling bearing faults [7,8]. In general, the method mainly consists of two steps: feature extraction and fault classification, in which valid feature extraction is crucial for accurate fault diagnosis. As the bearing vibration signal has nonlinear dynamic characteristics, traditional feature extraction methods based on Fourier transform and statistical analysis only characterize features from the time domain or frequency domain, and they cannot detect potential faults through changes in the complexity of the system to achieve effective and accurate extractions of fault features [9,10].
In recent years, nonlinear dynamic methods, such as sample entropy (SE) [11], permutation entropy (PE) [12], and dispersion entropy (DE) [13], have been widely used in the feature extraction of bearing signals and have presented superior performance. Han et al., used SE to extract bearing fault feature information effectively [14], but SE calculation is complicated and not suitable for real-time monitoring. As PE has the strengths of fast calculation and good stability, Xue et al., proposed a bearing fault diagnosis method based on PE and further improved the effectiveness of fault diagnosis [15]. While PE does not consider amplitude information, Dhandapani et al., applied DE to the feature extraction osf rolling bearing faults and considered the amplitude information of bearing signals [16]. Unlike the above entropies, slope entropy (SlopEn) is a new entropy estimator proposed based on symbolic patterns and magnitude information [17] and has been applied in the underwater acoustic field and medical field many times [18,19,20,21]. In 2022, SlopEn was introduced into the field of bearing fault diagnosis for the first time, and experimental results showed that, compared with PE and DE, SlopEn could better extract fault information [22]. However, all the above-mentioned entropy-based bearing fault diagnosis methods suffer from two defects: (i) the methods extract only the fault information of the low-frequency component for the bearing signal, and (ii) there is the problem of threshold selection for SlopEn, and the thresholds usually need to be optimized using an optimization algorithm.
Aiming at extracting the bearing fault information more comprehensively, some scholars have proposed fault diagnosis methods based on hierarchical entropy [23,24,25]. The authors of [26] proposed the concept of hierarchical permutation entropy (HPE) and employed it successfully for fault diagnosis. Moreover, Ref. [27] used hierarchical permutation entropy (HDE) to extract the fault information in both high- and low-frequency components. Diagnosis methods based on hierarchical entropy have proved that they can obtain diagnosis-related information of the whole frequency band and have strong noise resistance and stability; in addition, no scholars have introduced the concept of hierarchy to the SlopEn and used optimization algorithms to optimize the thresholds.
After feature extraction, the next step is fault classification. Commonly used fault classification methods mainly include k-nearest neighbor (KNN) [28], random forest (RF) [29], and the support vector machine (SVM) [30]. The SVM has been widely used in fault diagnosis because of its suitability for small sample classification and its simple structure [31]. Yet, since the parameter penalty factors and kernel functions of SVM models have an impact on fault diagnosis, some existing optimization algorithms have optimized the parameters of SVMs and improved the performance of fault classification, such as the genetic algorithm (GA) [32], particle swarm algorithm (PSO) [33], and whale optimization algorithm (WOA) [34]. Compared to common optimization algorithms, the white shark optimizer (WSO) is a new meta-heuristic optimization algorithm based on deep-sea foraging by great white sharks, proposed in 2022 for solving optimization problems on continuous search spaces [35]; in addition, the results of the basis function tests show that WSO is better than the common optimization algorithm in terms of optimization and has not yet been applied to optimize SVMs.
Based on the analysis above, a dual-optimization fault diagnosis method for rolling bearings is put forward, and the main novelties and contributions of this paper are presented as follows:
(1)
To extract the fault information of bearing signals more comprehensively, on the basis of SlopEn, this paper adds the concept of hierarchy and is the first to propose hierarchical slope entropy (HSlopEn).
(2)
Since the thresholds of HSlopEn have a relatively large impact on the entropy value and the selection of suitable parameters of an SVM is particularly important for the classification, this paper applies the WSO to optimize the parameters of HSlopEn and an SVM and proposes WSO-HSlopEn and WSO-SVM, respectively.
(3)
Targeting the application of bearing fault diagnosis under different operating conditions, this paper proposes a dual-optimization fault diagnosis method for rolling bearings based on HSlopEn and an SVM synergized with the WSO.
The remaining parts of this paper are structured as follows. Section 2 presents the basic concepts of algorithms. Section 3 introduces the steps of the proposed fault diagnosis method. Section 4 carries out the single-feature and multi-feature extraction experiments for bearing signals, and Section 5 summarizes the conclusions of this study.

2. Methodology

2.1. Slope Entropy

Slope entropy (SlopEn) is an algorithm proposed in 2019 to calculate the complexity of time series. It is based on symbolic patterns and magnitude information. The main calculation steps are listed below:
(1)
For a given time series X = x 1 , x 2 , , x N , according to the embedding dimension m , extract the subsequences: X 1 = x 1 , x 2 , , x m , X 2 = x 2 , x 3 , , x m 1 , , X k = x k , x k + 1 , , x N , of which k = N m + 1 .
(2)
Two threshold parameters, γ and δ , are used to delimit the symbolic patterns including (+2, +1, 0, −1, −2) between two consecutive samples x i + 1 x i of the subsequence X 1 , X 2 , , X k . When | x i + 1 x i | δ , the symbol is 0; when δ < x i + 1 x i < γ , the symbol is +1; when γ < x i + 1 x i < δ , the symbol is −1; when γ < x i + 1 x i , when symbol is +2; when x i + 1 x i < γ , the symbol is −2, where γ > δ > 0 . Figure 1 displays the symbol allocation of SlopEn.
(3)
According to the symbol pattern of step (2), divide all the subsequences X 1 , X 2 , , X k , and obtain the symbol pattern sequences Y 1 = y 1 , y 2 , , y m 1 , Y 2 = y 2 , y 3 , , y m , , Y k = y k , y k + 1 , , y N 1 , where y 1 is obtained from x 2 x 1 through step (2), y 2 is obtained from x 3 x 2 through step (2), … , and y N 1 is obtained from x N x N 1 through step (2).
(4)
The total number of types of symbol pattern sequences is recorded as r = 5 m 1 , and the number of occurrences of each symbol pattern sequence is s 1 , s 2 , , s r , respectively. Hence, the probability of each symbol pattern sequence appearing is P 1 = s 1 r , P 2 = s 2 r , , P r = s r r , respectively.
(5)
On the basis of Shannon information entropy, the definition formula of SlopEn is
S l o p E n ( X , m , γ , δ ) = i = 1 r P r l n P r

2.2. Hierarchical Slope Entropy

Since SlopEn only considers the low-frequency components of the time series, aiming to describe the time series more comprehensively, on the basis of SlopEn and combined with the concept of hierarchy, this paper proposes a new complexity feature, namely hierarchical slope entropy (HSlopEn). The specific process of HSlopEn is as follows:
(1)
First, given a time series X = x i , i = 0 , 1 , , N , N = 2 n of length N , define an average operator Q 0 x and a difference operator Q 1 x , which can be expressed as
Q 0 = x 2 j + x 2 j + 1 2 , j = 0 , 1 , 2 , , 2 n 1
Q 1 = x 2 j x 2 j + 1 2 , j = 0 , 1 , 2 , , 2 n 1
where the Q 0 x and Q 1 x operators are the low-frequency part and high-frequency part, respectively, of the original given time series after hierarchical decomposition and n is a positive integer.
(2)
The operators Q j ( j = 0 o r 1 ) in matrix form is defined as
Q j = 1 2 1 2 j 0 0 0 0 0 0 1 2 1 2 j 0 0 0 0 0 0 1 2 1 2 j 2 n 1 × 2 n
(3)
The l -dimension vectors u 1 , u 2 , , u l { 0,1 } ( l ϵ N ) are constructed, and an integer e can be expressed:
e = j = 1 l u j 2 l j
where, for a positive e , there is a unique set of l -dimension vectors u 1 , u 2 , , u l { 0,1 } and the positive integer e represents the sequence number of the node at each layer, where 0 e 2 n 1 .
(4)
The hierarchical decomposition of a given time series X yields a hierarchical component corresponding to the node e at the K th level, defined as
X K , e = Q u l Q u l 1 Q u 1 x X
(5)
By calculating the SlopEn of nodes on different layers, the HSlopEn can be expressed as
H S l o p E n X , m , K , γ , δ = S l o p E n X n , e , m , γ , δ
where K is the number of layers of decomposition, γ and δ are the two thresholds of SlopEn, and m is the embedding dimension.
As displayed in Figure 2, the hierarchical decomposition structure diagram when K = 3 is shown. SlopEn is calculated on each node after the hierarchical decomposition.
In Figure 2, X indicates the original time series, x 1,1 is the first node of the first layer, x 2,1 is the first node of the second layer, and so on.

2.3. Analysis of the Parameters for HSlopEn

The main parameters of HSlopEn include the number of decomposition layers K , embedding dimension m , two threshold parameters γ a n d δ , and time delay d . First, the number of decomposition layers K determines the number of nodes in the hierarchical decomposition. When the number of decomposition layers is too large, the number of nodes decomposed is too large, resulting in a large number of calculations for SlopEn values of all nodes; when the value is too small, resulting in a small number of decomposed nodes, there are insufficient frequency bands for the given time series. Referring to other references, the default number of decomposition layers K is 3 in this paper. Then, the embedding dimension m is used to extract the subsequence of a given time series. If it is too small, it is difficult to determine the dynamic changes of the time series; if it is too large, it is difficult to capture the subtle changes in the time series. After that, the two threshold parameters γ a n d δ are used to divide the symbol pattern of a given subsequence, which affects the change in entropy value. Lastly, the default time delay d is 1, as important information about frequency may be lost at that time if d > 1 . The effect of embedding dimension and thresholds on the performance of the HSlopEn is investigated below by analyzing the noisy signals.
To investigate the effect of embedding dimension on the entropy value of hierarchical slope entropy, 50 sets of white Gaussian noise (WGN) of signal length 2048 are used, with the embedding dimension m varying from 2 to 5 and the two threshold parameters γ a n d δ defaulting to 0.1 and 0.001, respectively. Figure 3 shows the mean and standard deviation (SD) of the HSlopEn values for different embedding dimensions in every node.
As shown in Figure 3, as the embedding dimension m becomes larger, the entropy value of the HSlopEn also becomes larger, but the entropy value of each node for HSlopEn is close to others at different embedding dimensions, and the difference between the mean and SD is small, which indicates that the change in the embedding dimension affects the size of the entropy value, but the stability of the HSlopEn hardly changes. The embedding dimension m is set to 3 in this paper.
In addition, to further study the effect of thresholds γ a n d δ on the entropy of the HSlopEn, 50 independent pink noise (PN) and WGN signals are selected, where each noise is sampled at 2048 Hz and the embedding dimension m is 3. The three sets of thresholds ( γ , δ ) for HSlopEn are manually set, which are (0.1, 001), (0.3, 0.1), and (0.8, 0.3), and the mean and standard deviation (SD) of the HSlopEn values for the three sets of thresholds in every node are displayed in Figure 4.
It can be seen in Figure 4 that, as the threshold changes, the entropy values of the two types of noise signals change; at the same time, the ability to discriminate between the noise signals is constantly changing, so the threshold has a significant effect on the entropy of the HSlopEn. The WSO is used in the paper to optimize the thresholds to avoid taking values based on artificial experience and further improve the fault diagnosis.

2.4. WSO-HSlopEn and WSO-SVM

Following the principle of the HSlopEn algorithm, the two threshold parameters γ and δ of the HSlopEn are used to divide the sign pattern of a given time sequence subsequence. Thus, the two threshold parameters have a great influence on the HSlopEn value. At the same time, the classification effect of the support vector machine (SVM) mainly depends on the selection of the penalty factor ( C ) and kernel function parameters ( g ), and it is generally difficult to take the values based on manual experience. Hence, the selection of an appropriate penalty factor and kernel function parameters is also particularly important for the classification and recognition accuracy of the SVM.
To enhance the performance of the fault diagnosis effect, in this paper, taking the recognition rate as the fitness function, the white shark optimizer (WSO) is used to optimize the parameters of HSlopEn and the SVM, and WSO-HSlopEn and WSO-SVM are proposed, respectively, where the WSO is a new meta-heuristic optimization algorithm based on deep-sea foraging by great white sharks, proposed in 2022 for solving optimization problems on continuous search spaces. The main process of optimizing the parameters of HSlopEn and the SVM is shown in Figure 5, and the specific process is as follows:
(1)
Set the initial parameter ranges of HSlopEn ( γ , δ ) and the SVM ( C , g ) .
(2)
Initialize the WSO parameters, such as population size, number of iterations I , position, and speed of white sharks.
(3)
Calculate the fitness function, and update the white sharks’ position and speed.
(4)
Evaluate the fitness function, and update the optimal white shark position
(5)
Update the position and speed of the white shark.
(6)
Judge whether the current iteration number reaches the maximum iteration number. If so, return to update the speed and position of the white shark and repeat the above steps; otherwise, output the best-optimized parameters ( γ , δ ) and ( C , g ) .

3. The Proposed Method for Fault Diagnosis of Rolling Bearing

Combining the concept of hierarchical structure, the new complexity feature HSlopEn is proposed, and the parameters of both HSlopEn and the SVM are optimized using the WSO algorithm, and WSO-HSlopEn and WSO-SVM are proposed, respectively. Then, a dual-optimization fault diagnosis method for rolling bearings based on WSO-HSlopEn and WSO-SVM is proposed. Figure 6 presents the flowchart of the proposed fault diagnosis method, and the method mainly includes the following steps:
(1)
The different bearing signals are input. In this paper, each type of bearing signal has 100 samples with 1024 data points.
(2)
The WSO algorithm is applied to optimize the parameters of HSlopEn ( γ , δ ) and the SVM ( C , g ) by taking the final recognition rate as the fitness function, and the optimized parameters are obtained. At the same time, other optimization algorithms, including SO, marine predator algorithm (MPA), and sparrow search algorithm (SSA), are used for comparison.
(3)
Different types of bearing signals are decomposed into several layers, and the nodes are obtained. In this paper, bearing signals are decomposed into three layers.
(4)
The nodes of WSO-HSlopEn are calculated, and then single-feature and multi-feature extraction experiments for bearing signals are carried out. Meanwhile, comparisons with some classical entropies, such as HFE, HPE, HSE, and HRDE, are conducted.
(5)
WSO-SVM is applied to classify bearing signals, and the recognition results are output. In this paper, for each type, select 25 sample signals as training samples and 75 sample signals as test samples.

4. Experiments and Results

In this chapter, two comparative experiments are implemented to examine the effectiveness of the proposed method in fault diagnosis: (1) In optimizing both HSlopEn and SVM parameters using the WSO, we compare different optimization algorithms, including SSA, MPA, and SO. (2) In extracting the WSO-HSlopEn of nodes, we compare classical hierarchical entropy metrics, including HPE, HSE, HFE, and HRDE.

4.1. Fault Diagnosis of Rolling Bearing Signal

The dataset used in this section was derived from the Bearing Data Center of Case Western Reserve University [36], which is an internationally recognized standard dataset for fault diagnosis of rolling bearings. The schematic of the test rig (Cleveland, USA) is shown in Figure 7.
As shown in Figure 7, the test rig consisted of an induction motor, drive-end bearing, self-aligning coupling, and accelerometer dynamometer. An accelerometer was installed on the base of the motor, which was used to detect the vibration acceleration of the faulty bearing at a sampling frequency of 12 kHz. The dataset divided the fault data into four categories: normal data (NOR), ball faults (BFs), outer race faults (ORFs), and inner race faults (IRFs). Among them, BFs, ORFs, and IRFs were simulated faults with single-point damage as an electric spark. The damage diameters were divided into 0.007, 0.014, and 0.021 inches. At the same time, the processed faulty bearing was reloaded into the test motor, and the vibration acceleration signal data were recorded under the load working conditions of 0, 1, 2, and 3 horsepower.
In this section, bearing signals with ten conditions were collected from the drive-end bearings, including rolling bearings in normal condition and those with damage to the inner race, the outer race, and the ball element. Bearings with various damage diameters were considered under a speed of 1730 rpm with a load of 3 horsepower. Table 1 illustrates the fault diagnosis sample collection of bearing signals. Each fault signal was divided into three types according to the fault diameter. We sampled from point 1001, and each condition had 100 samples with 1024 sampling points. Time-domain waveforms for each state bearing signals are displayed in Figure 8.

4.2. Comparison of Different Optimization Algorithms

Designed to verify the performance advantages of the WSO in optimizing HSlopEn and the SVM, this section introduces different optimization algorithms to optimize the parameters of HSlopEn and the SVM, and compares recognition rates of single-feature and multi-feature extractions with those of other optimization algorithms [37,38,39]. In this experiment, 10 different bearing signal conditions were sampled from the 1001 point as the starting point, and 100 samples were selected. Each sample had 1024 data points. First, the parameters of HSlopEn were set as follows: hierarchical layer K = 3 , embedding dimension m = 3 , and threshold parameters γ and δ were adaptively determined using different optimization algorithms. HSlopEn with optimized parameters of bearing signals was extracted. Then, the sample set was divided into the training set and test set, and the select single feature or multi-features were input to optimize the SVM. The penalty factor and kernel function parameters of the SVM were also adaptively determined using the WSO algorithm. Figure 9 presents the fitness iteration curves of different optimization algorithms to optimize HSlopEn and the SVM. These are the fitness iteration curves of different optimization algorithms in the case of extracting five nodes.
It can be found in Figure 9 that, in the condition of extracting five nodes, the highest recognition rate of these ten types of bearing signals reached 100% using the WSO to optimize HSlopEn. At the same time, the convergence speed of the WSO was relatively faster than other optimization algorithms. In addition, the early convergence of the WSO is quick. Its fitness curve eventually converged to a bigger value. To further demonstrate the significant advantages of WSO, we calculated the recognition rate of bearing signals based on using different optimization algorithms to optimize HSlopEn and the SVM under the situation of extracting the single feature and multi-features. The recognition rates of HSlopEn for all single-feature nodes are shown in Table 2 and Table 3, and the highest recognition rates of HSlopEn for multi-features for the four types of optimization algorithms are shown in Table 4.
According to the recognition rate of different types of bearing signals, we can find that no matter how many features are extracted, the advantages of the WSO algorithm are obvious. In the case of extracting a single feature, the recognition rate of the WSO reaches 79.33% on node 6, which is much higher than that of other optimization algorithms. Under the circumstances of extracting multi-features, as the number of selected nodes increases, the recognition rate also improves. When we select five features, it realizes the correct identification of all samples. The recognition rate of other optimization algorithms, including SO, MPA, and SSA, is, respectively, 3.8%, 10.53%, and 17.73% lower than that of WSO. Above all, we prove that using the WSO to optimize HSlopEn and the SVM is feasible. Therefore, in this paper, the WSO is used to optimize HSlopEn and SVM parameters.

4.3. Comparison of Different Hierarchical Entropies

Aiming to demonstrate the superiority of WSO-HSlopEn in fault diagnosis, we compared it to other classical hierarchical entropies, including HSE, HFE, HPE, and HRDE. The single-feature approach was first used to extract the fault feature and compare it with HFE, HSE, HPE, and HRDE. The parameters of HSlopEn were as follows: hierarchical layer K = 3 , embedding dimension m = 3 , time delay d = 1 , and threshold parameter γ and δ were adaptively determined using the WSO algorithm. For a fair comparison, the parameter settings of other hierarchical entropies were the same as those in the HSlopEn method. Among them, the similarity tolerances of HSE and HFE were set as r = 0.2 , and the category number of HRDE was set as c = 3 . The entropy distributions of an optimal node for the single-feature extraction of bearing signals are shown in Figure 10.
Figure 10a presents that there was no aliasing phenomenon between the features of the three fault types of NOR, ORF3, and IRF1 and other fault types in the entropy distribution of WSO-HSlopEn. ORF2 only had a few samples of entropy close to BF2, and the entropy values of the samples of the other fault types show a severe overlap. Based on the single-feature extraction shown in Figure 10b–e, compared with other hierarchical entropies, WSO-HSlopEn is not as serious as the aliasing of the distributions of other hierarchical entropies. The entropy values of several types of samples in the ten conditions of bearing signals are quite different from the other types. Relatively speaking, the distance between the various types of samples in the distribution of the entropy of WSO-HSlopEn is relatively large.
After using the WSO-HSlopEn as the fault feature of the bearing signal, the bearing fault diagnosis sample set was divided into a training set and a test set, and the training set was input into the WSO-SVM to train the model, and then the test set was input into the model to finish the fault diagnosis of bearings. The Gaussian kernel function was selected as the kernel function of the SVM. The penalty factor and kernel function parameters of the SVM were also adaptively determined by the WSO algorithm. Recognition rates of single features for the five types of hierarchical entropies are displayed in Table 5 and Table 6.
Table 5 and Table 6 illustrate that, when using WSO-HSlopEn, the recognition rate of node 6 was the highest, which was 79.33%. Compared with other hierarchical entropies, under each node, the recognition rate based on WSO-HSlopEn was always the highest, which shows the effectiveness of WSO-HSlopEn as a fault diagnosis feature of bearing signals.
Through observation, when single-feature extraction is used to extract the fault feature, there is still overlap between the features of different conditions of the bearing signals. Furthermore, the recognition rate of the best node was low, and there were many misclassified samples based on single-feature extraction. Aiming to further improve the recognition rate of different conditions of the bearing signals, double features were used to extract the bearing signals. All parameters used in the experiments were the same as those listed in the single-feature extraction. The entropy distribution on the optimal node for double-feature extractions of bearing signals is shown in Figure 11, where the abscissa and ordinate are the entropy values of the two nodes, respectively. For example, in Figure 11a, the abscissa is the SlopEn of node 1, and the ordinate is the SlopEn of node 5.
As can be observed from Figure 11, in the case of double-feature extraction, the WSO-HSlopEn distribution of sample signals belonging to the same type is relatively concentrated compared to other hierarchical entropies; for the other four types of hierarchical entropies, the bearing signals between different types are more divergent, and the entropy values of different types of bearing signals are very close.
To further improve the recognition performance, triple features were used to extract bearing fault features on various hierarchical entropies. The parameters for calculating various hierarchical entropies were the same as those of double features. Figure 12 presents the triple-feature distributions of ten types of bearing signals for different hierarchical entropies.
It can be seen from Figure 12 that there is almost no overlap based on the WSO-HSlopEn, but the feature distributions of the BF2 and IRF2 samples are relatively low in clustering; for the other hierarchical entropies, the clustering of the feature distributions of the samples are very poor because of their approximate entropy distributions. Nevertheless, the entropy distribution of WSO-HSlopEn is more dispersed, and WSO-HSlopEn of different fault types are quite different, which effectively verifies the validity of WSO-HSlopEn as a feature extraction method for ten types of bearing signals.
Next, WSO-SVM is used to construct a fault diagnosis model. The highest recognition rate is calculated for the five types of hierarchical entropies under multi-feature extraction, as shown in Table 7, where (1,5) indicates the combination of nodes with the highest recognition rate for two features are node 1 and node 5, (1,5,6) indicates the combination of nodes with the highest recognition rate for three features are node 1, node 5 and node 6, and so on.
Table 7 shows that no matter how many features are extracted, the recognition rate of these ten types of bearing signals using WSO-HSlopEn is higher than that of other hierarchical entropies; additionally, the more features we select, the better the recognition effect we obtain; in the circumstances of multi-features, the recognition rates of WSO-HSlopEn are all higher than 97.5%, yet the highest recognition rates of other hierarchical entropies are all significantly below 97.5%; for WSO-HSlopEn, when five nodes are selected, that is, choosing nodes (1,5,6,7,11), the highest recognition rate of these ten types of bearing signals reaches 100%; however, the highest recognition rate of other entropies is, respectively, 3.80%, 10.53%, 16.73%, and 4.13% lower than that of WSO-HSlopEn. Through the above comparison, we can clearly find the significant advantages of the proposed method based on WSO-HSlopEn, and the recognition results applied to diagnose faults of rolling bearings are higher than those of classic methods.

5. Conclusions

This paper puts forward a dual-optimization fault diagnosis method for rolling bearings based on WSO-HSlopEn and WSO-SVM. The effectiveness of the proposed methods is verified by comparing them with the classical methods. The main innovations and conclusions are as follows:
(1)
On the basis of SlopEn, combined with the idea of hierarchical decomposition, HSlopEn is proposed and introduced into the feature extraction of bearing signals for the first time; at the same time, WSO is used to optimize both HSlopEn and the SVM, and WSO-HSlopEn and WSO-SVM are proposed.
(2)
In the case of single-feature extraction, the proposed method based on WSO-HSlopEn has the highest recognition rate of 79.33% on node 6, which is, respectively, 19.47%, 26.67%, 31.60%, and 13.73% higher than those of HFE, HPE, HSE, and HRDE.
(3)
In the case of extracting multi-features, the recognition rates are higher than 97.5%, which is a significant improvement compared with the single-feature extraction method; moreover, with the different number of features, the recognition rate based on WSO-HSlopEn is always high than the other hierarchical entropies.
(4)
For the proposed dual-optimization fault diagnosis method for rolling bearings, based on WSO-HSlopEn and WSO-SVM, the more features we select, the better the recognition effect we obtain. When five nodes are selected, the highest recognition rate reaches 100%.
The proposed WSO-HSlopEn and WSO-SVM solve the problem of dependent parameter settings for SlopEn and the SVM, respectively, and their superiority has been confirmed in fault diagnosis. Therefore, WSO-HSlopEn and WSO-SVM are expected to be applied to other fields in future work, such as underwater acoustic signal processing and medical signal classification.

Author Contributions

Conceptualization, Y.L.; Methodology, Y.L.; Validation, B.T.; Data curation, B.H.; Writing—original draft, B.T.; Visualization, B.T., B.H. and X.X.; Supervision, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Natural Science Foundation of Shaanxi Province] grant number [2022JM-337], [Xi’an University of Technology Excellent Seed Fund] grant number [252082220].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. Li, X.; Xu, Y.; Li, N.; Yang, B.; Lei, Y. Remaining Useful Life Prediction with Partial Sensor Malfunctions Using Deep Adversarial Networks. IEEE CAA J. Autom. Sin. 2023, 10, 121–134. [Google Scholar] [CrossRef]
  3. Saucedo-Dorantes, J.J.; Arellano-Espitia, F.; Delgado-Prieto, M.; Osornio-Rios, R.A. Diagnosis Methodology Based on Deep Feature Learning for Fault Identification in Metallic, Hybrid and Ceramic Bearings. Sensors 2021, 21, 5832. [Google Scholar] [CrossRef]
  4. Li, Y.; Tang, B.; Geng, B.; Jiao, S. Fractional Order Fuzzy Dispersion Entropy and Its Application in Bearing Fault Diagnosis. Fractal Fract. 2022, 6, 544. [Google Scholar] [CrossRef]
  5. Li, X.; Yu, S.; Lei, Y.; Li, N.; Yang, B. Intelligent Machinery Fault Diagnosis with Event-Based Camera. IEEE Trans. Ind. Inform. 2023, 1–10. [Google Scholar] [CrossRef]
  6. Yadav, E.; Chawla, V.K. An explicit literature review on bearing materials and their defect detection techniques. Mater. Today Proc. 2022, 50, 1637–1643. [Google Scholar] [CrossRef]
  7. Hu, Q.; Li, Y.; Sun, X.; Chen, M.; Bu, Q.; Gong, B. Integrating test device and method for creep failure and ultrasonic response of methane hydrate-bearing sediments. Rev. Sci. Instrum. 2023, 94, 025105. [Google Scholar] [CrossRef]
  8. Li, Y.; Tang, B.; Jiang, X.; Yi, Y. Bearing Fault Feature Extraction Method Based on GA-VMD and Center Frequency. Math. Probl. Eng. 2022, 2022, 2058258. [Google Scholar] [CrossRef]
  9. Li, Y.; Jiao, S.; Geng, B. Refined composite multiscale fluctuation-based dispersion Lempel–Ziv complexity for signal analysis. ISA Trans. 2023, 133, 273–284. [Google Scholar] [CrossRef]
  10. Li, Y.; Geng, B.; Tang, B. Simplified coded dispersion entropy: A nonlinear metric for signal analysis. Nonlinear Dyn. 2023, 111, 9327–9344. [Google Scholar] [CrossRef]
  11. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 6, 2039–2049. [Google Scholar] [CrossRef] [Green Version]
  12. Bandt, C.; Pompe, B. Permutation entropy: A natural complexity measure for time series. Phys. Rev. Lett. 2002, 88, 174102. [Google Scholar] [CrossRef]
  13. Rostaghi, M.; Azami, H. Dispersion Entropy: A Measure for Time Series Analysis. IEEE Signal Process. Lett. 2016, 23, 610–614. [Google Scholar] [CrossRef]
  14. 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]
  15. Xue, X.; Li, C.; Cao, S.; Sun, J.; Liu, L. Fault Diagnosis of Rolling Element Bearings with a Two-Step Scheme Based on Permutation Entropy and Random Forests. Entropy 2019, 21, 96. [Google Scholar] [CrossRef] [Green Version]
  16. Dhandapani, R.; Mitiche, I.; McMeekin, S.; Mallela, V.S.; Morison, G. Enhanced Partial Discharge Signal Denoising Using Dispersion Entropy Optimized Variational Mode Decomposition. Entropy 2021, 23, 1567. [Google Scholar] [CrossRef] [PubMed]
  17. Cuesta-Frau, D. Slope Entropy: A New Time Series Complexity Estimator Based on Both Symbolic Patterns and Amplitude Information. Entropy 2019, 21, 1167. [Google Scholar] [CrossRef] [Green Version]
  18. Li, Y.; Tang, B.; Jiao, S. SO-slope entropy coupled with SVMD: A novel adaptive feature extraction method for ship-radiated noise. Ocean Eng. 2023, 280, 114677. [Google Scholar] [CrossRef]
  19. Li, Y.; Tang, B.; Yi, Y. A novel complexity-based mode feature representation for feature extraction of ship-radiated noise using VMD and slope entropy. Appl. Acoust. 2022, 196, 108899. [Google Scholar] [CrossRef]
  20. Cuesta-Frau, D.; Dakappa, P.H.; Mahabala, C.; Gupta, A.R. Fever Time Series Analysis Using Slope Entropy. Application to Early Unobtrusive Differential Diagnosis. Entropy 2020, 22, 1034. [Google Scholar] [CrossRef]
  21. Cuesta-Frau, D.; Schneider, J.; Bakštein, E.; Vostatek, P.; Spaniel, F.; Novák, D. Classification of Actigraphy Records from Bipolar Disorder Patients Using Slope Entropy: A Feasibility Study. Entropy 2020, 22, 1243. [Google Scholar] [CrossRef]
  22. Shi, E. Single Feature Extraction Method of Bearing Fault Signals Based on Slope Entropy. Shock. Vib. 2022, 2022, 6808641. [Google Scholar] [CrossRef]
  23. Jiang, Y.; Peng, C.K.; Xu, Y. Hierarchical entropy analysis for biological signals. J. Comput. Appl. Math. 2011, 236, 728–742. [Google Scholar] [CrossRef] [Green Version]
  24. Xing, J.; Xu, J. An Improved Incipient Fault Diagnosis Method of Bearing Damage Based on Hierarchical Multi-Scale Reverse Dispersion Entropy. Entropy 2022, 24, 770. [Google Scholar] [CrossRef]
  25. Peng, C.; Zhao, X.; Jiang, H. A New Method of Fault Feature Extraction Based on Hierarchical Dispersion Entropy. Shock. Vib. 2021, 2021, 8824901. [Google Scholar]
  26. Wang, X.; Si, S.; Li, Y. An integrated method based on refined composite multivariate hierarchical permutation entropy and random forest and its application in rotating machinery. J. Vib. Control 2020, 26, 3–4. [Google Scholar] [CrossRef]
  27. Xue, Q.; Xu, B.; He, C.; Liu, F. Feature Extraction Using Hierarchical Dispersion Entropy for Rolling Bearing Fault Diagnosis. IEEE Trans. Instrum. Meas. 2021, 70, 3521311. [Google Scholar] [CrossRef]
  28. Keller, J.M.; Gray, M.R.; Givens, J.A. A fuzzy K-nearest neighbor algorithm. IEEE Trans. Syst. Man Cybern. 1985, SMC-15, 580–585. [Google Scholar] [CrossRef]
  29. Svetnik, V.; Liaw, A.; Tong, C.; Culberson, J.C.; Sheridan, R.P.; Feuston, B.P. Random Forest:  A Classification and Regression Tool for Compound Classification and QSAR Modeling. J. Chem. Inf. Comput. Sci. 2003, 43, 1947–1958. [Google Scholar] [CrossRef]
  30. Samanta, B.; Al-Balushi, K.R.; Al-Araimi, S.A. Artificial neural networks and support vector machines with genetic algorithm for bearing fault detection. Eng. Appl. Artif. Intell. 2003, 16, 657–665. [Google Scholar] [CrossRef]
  31. Wang, Z.; Yao, L.; Cai, Y.; Zhang, J. Mahalanobis semi-supervised mapping and beetle antennae search based support vector machine for wind turbine rolling bearings fault diagnosis. Renew. Energy 2020, 155, 1312–1327. [Google Scholar] [CrossRef]
  32. Tang, H.; Yuan, Z.; Dai, H.; Du, Y. Fault Diagnosis of Rolling Bearing Based on Probability box Theory and GA-SVM. IEEE Access 2020, 8, 170872–170882. [Google Scholar] [CrossRef]
  33. 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]
  34. Jin, Z.; Chen, G.; Yang, Z. Rolling Bearing Fault Diagnosis Based on WOA-VMD-MPE and MPSO-LSSVM. Entropy 2022, 24, 927. [Google Scholar] [CrossRef]
  35. Braik, M.; Hammouri, A.; Atwan, J. White Shark Optimizer: A novel bioinspired meta-heuristic algorithm for global optimization problems. Knowl. Based Syst. 2022, 243, 108457. [Google Scholar] [CrossRef]
  36. Case Western Reserve University Bearing Data Center. Available online: https://engineering.case.edu/bearingdatacenter (accessed on 26 February 2023).
  37. Fatma, A.H.; Abdelazim, G.H. Snake Optimizer: A novel meta-heuristic optimization algorithm. Knowl. Based Syst. 2022, 242, 108320. [Google Scholar]
  38. Afshin, F.; Mohammad, H.; Seyedali, M.; Amir, H.G. Marine Predators Algorithm: A nature-inspired metaheuristic. Expert Syst. Appl. 2020, 152, 113377. [Google Scholar]
  39. Xue, J.; Shen, B. A novel swarm intelligence optimization approach: Sparrow search algorithm. Syst. Sci. Control Eng. 2020, 8, 22–34. [Google Scholar] [CrossRef]
Figure 1. The symbol allocation of SlopEn.
Figure 1. The symbol allocation of SlopEn.
Sensors 23 05630 g001
Figure 2. Hierarchical decomposition structure diagram when K = 3.
Figure 2. Hierarchical decomposition structure diagram when K = 3.
Sensors 23 05630 g002
Figure 3. The mean and standard deviation (SD) of HSlopEn values for different embedding dimensions in every node.
Figure 3. The mean and standard deviation (SD) of HSlopEn values for different embedding dimensions in every node.
Sensors 23 05630 g003
Figure 4. The mean and SD of HSlopEn values for three sets of thresholds in every node.
Figure 4. The mean and SD of HSlopEn values for three sets of thresholds in every node.
Sensors 23 05630 g004
Figure 5. The main process of optimizing the parameters of the HSlopEn and SVM.
Figure 5. The main process of optimizing the parameters of the HSlopEn and SVM.
Sensors 23 05630 g005
Figure 6. The flowchart of proposed dual-optimization fault diagnosis method for rolling bearings, based on WSO-HSlopEn and WSO-SVM.
Figure 6. The flowchart of proposed dual-optimization fault diagnosis method for rolling bearings, based on WSO-HSlopEn and WSO-SVM.
Sensors 23 05630 g006
Figure 7. The schematic of test rig.
Figure 7. The schematic of test rig.
Sensors 23 05630 g007
Figure 8. Time-domain waveforms for each state bearing signals.
Figure 8. Time-domain waveforms for each state bearing signals.
Sensors 23 05630 g008
Figure 9. The fitness iteration curve of different optimization algorithms.
Figure 9. The fitness iteration curve of different optimization algorithms.
Sensors 23 05630 g009
Figure 10. Entropy distribution of optimal node for single-feature extraction of bearing signals.
Figure 10. Entropy distribution of optimal node for single-feature extraction of bearing signals.
Sensors 23 05630 g010aSensors 23 05630 g010b
Figure 11. Double-features distribution of ten types of bearing signals.
Figure 11. Double-features distribution of ten types of bearing signals.
Sensors 23 05630 g011aSensors 23 05630 g011b
Figure 12. Triple-features distribution of ten types of bearing signals for different entropies.
Figure 12. Triple-features distribution of ten types of bearing signals for different entropies.
Sensors 23 05630 g012aSensors 23 05630 g012b
Table 1. Fault diagnosis sample collection of bearing signals.
Table 1. Fault diagnosis sample collection of bearing signals.
Bearing ConditionDefect Size (Inches)Sample Label
NOR-1
IRF10.0072
IRF20.0143
IRF30.0214
BF10.0075
BF20.0146
BF30.0217
ORF10.0078
ORF20.0149
ORF30.02110
Table 2. The recognition rate of HSlopEn for each single-feature node (nodes 1–7).
Table 2. The recognition rate of HSlopEn for each single-feature node (nodes 1–7).
Optimization AlgorithmRecognition Rate for Each Node (%)
1234567
WSO78.4079.2070.1378.0073.4779.3345.73
SO77.4774.8067.8776.3272.2775.4343.28
MPA74.0075.8762.9371.2066.9376.5340.40
SSA74.6773.0768.5371.8770.6775.6045.07
Table 3. The recognition rate of HSlopEn for each single-feature node (nodes 8–14).
Table 3. The recognition rate of HSlopEn for each single-feature node (nodes 8–14).
Optimization AlgorithmsRecognition Rate for Each Node (%)
891011121314
WSO51.6072.2759.6066.9362.8069.0770.93
SO48.3770.2658.2365.2760.8068.2769.73
MPA39.8770.8056.5360.1356.1367.8765.47
SSA50.4069.7357.7363.8760.9368.4269.57
Table 4. The highest recognition rate of HSlopEn for multi-features for four types of optimization algorithms.
Table 4. The highest recognition rate of HSlopEn for multi-features for four types of optimization algorithms.
Optimization AlgorithmsParameterNumber of Extracted Features
2345
WSOHighest recognition rate (%)97.8799.6099.87100
SOHighest recognition rate (%)87.6090.8093.8796.20
MPAHighest recognition rate (%)67.8784.3786.6789.47
SSAHighest recognition rate (%)64.1370.4074.6783.27
Table 5. Recognition rates of single features for the five types of hierarchical entropies (nodes 1–7).
Table 5. Recognition rates of single features for the five types of hierarchical entropies (nodes 1–7).
EntropiesRecognition Rates of Each Node (%)
1234567
WSo-HSlopEn78.4079.2070.1378.0073.4779.3345.73
HFE57.0755.0637.7355.3342.6759.7332.53
HPE43.7352.5324.8040.9338.1346.5318.80
HSE39.0747.6023.8736.2721.0746.6718.40
HRDE55.7365.4733.8753.8742.4054.5330.67
Table 6. Recognition rates of single features for the five types of hierarchical entropies (nodes 8–14).
Table 6. Recognition rates of single features for the five types of hierarchical entropies (nodes 8–14).
EntropiesRecognition Rates of Each Node (%)
891011121314
WSo-HSlopEn51.6072.2759.6066.9362.8069.0770.93
HFE27.0743.7336.4038.9336.4045.7340.27
HPE14.8029.6022.2729.0720.8032.2726.40
HSE20.4028.6720.4021.6020.5337.3325.87
HRDE25.0745.6037.8742.9337.0750.1348.93
Table 7. Highest recognition rate for the five types of hierarchical entropies for multi-features.
Table 7. Highest recognition rate for the five types of hierarchical entropies for multi-features.
EntropyParameterNumber of Extracted Nodes
2345
WSo-HSlopEnHighest recognition rate (%)97.8799.6099.87100
Choose the node(1,5)(1,5,6)(1,5,6,7)(1,5,6,7,11)
HFEHighest recognition rate (%)87.6090.8093.8796.20
Choose the node(2,6)(2,5,13)(1,3,4,6)(2,5,6,12,14)
HPEHighest recognition rate (%)67.8784.3786.6789.47
Choose the node(2,4)(1,2,14)(2,3,5,10)(1,5,7,12,14)
HSEHighest recognition rate (%)64.1370.4074.6783.27
Choose the node(2,5)(1,2,6)(1,2,3,6)(2,3,7,12,14)
HRDEHighest recognition rate (%)80.6089.8793.8795.87
Choose the node(1,2)(1,2,3)(1,2,4,6)(2,5,6,11,12)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Y.; Tang, B.; Huang, B.; Xue, X. A Dual-Optimization Fault Diagnosis Method for Rolling Bearings Based on Hierarchical Slope Entropy and SVM Synergized with Shark Optimization Algorithm. Sensors 2023, 23, 5630. https://doi.org/10.3390/s23125630

AMA Style

Li Y, Tang B, Huang B, Xue X. A Dual-Optimization Fault Diagnosis Method for Rolling Bearings Based on Hierarchical Slope Entropy and SVM Synergized with Shark Optimization Algorithm. Sensors. 2023; 23(12):5630. https://doi.org/10.3390/s23125630

Chicago/Turabian Style

Li, Yuxing, Bingzhao Tang, Bo Huang, and Xiaohui Xue. 2023. "A Dual-Optimization Fault Diagnosis Method for Rolling Bearings Based on Hierarchical Slope Entropy and SVM Synergized with Shark Optimization Algorithm" Sensors 23, no. 12: 5630. https://doi.org/10.3390/s23125630

APA Style

Li, Y., Tang, B., Huang, B., & Xue, X. (2023). A Dual-Optimization Fault Diagnosis Method for Rolling Bearings Based on Hierarchical Slope Entropy and SVM Synergized with Shark Optimization Algorithm. Sensors, 23(12), 5630. https://doi.org/10.3390/s23125630

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