Next Article in Journal
Nonlinear Neutral Delay Differential Equations of Fourth-Order: Oscillation of Solutions
Next Article in Special Issue
An Information-Theoretic Perspective on Proper Quaternion Variational Autoencoders
Previous Article in Journal
An Entropy Metric for Regular Grammar Classification and Learning with Recurrent Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy and Its Application to Multivariate Signal of Rotating Machinery

1
Institute of Agricultural Machinery, Hubei University of Technology, Wuhan 430068, China
2
Hubei Engineering Research Center for Intellectualization of Agricultural Equipment, Wuhan 430068, China
*
Author to whom correspondence should be addressed.
Submission received: 4 November 2020 / Revised: 6 January 2021 / Accepted: 7 January 2021 / Published: 19 January 2021

Abstract

:
In the fault monitoring of rotating machinery, the vibration signal of the bearing and gear in a complex operating environment has poor stationarity and high noise. How to accurately and efficiently identify various fault categories is a major challenge in rotary fault diagnosis. Most of the existing methods only analyze the single channel vibration signal and do not comprehensively consider the multi-channel vibration signal. Therefore, this paper presents Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy (RCMMFDE), a method which extracts the recognition information of multi-channel signals with different scale factors, and the refined composite analysis ensures the recognition stability. The simulation results show that this method has the characteristics of low sensitivity to signal length and strong anti-noise ability. At the same time, combined with Joint Mutual Information Maximisation (JMIM) and support vector machine (SVM), RCMMFDE-JMIM-SVM fault diagnosis method has been proposed. This method uses RCMMFDE to extract the state characteristics of the multiple vibration signals of the rotary machine, and then uses the JMIM method to extract the sensitive characteristics. Finally, different states of the rotary machine are classified by SVM. The validity of the method is verified by the composite gear fault data set and bearing fault data set. The diagnostic accuracy of the method is 99.25% and 100.00%. The experimental results show that RCMMFDE-JMIM-SVM can effectively recognize multiple signals.

1. Introduction

Rotating machinery plays an important role in the modern industry with a complex dynamic system. During its operation, its bearings and gears may fail due to fatigue, wear and corrosion. Due to the influence of friction, load and impact, its vibration signal often appears unsteady and nonlinear characteristics [1,2,3]. Therefore, a variety of nonlinear signal analysis methods have been widely used in bearing fault diagnosis. Due to its unique advantages in feature extraction, more and more attention has been paid to entropy by researchers in ever-increasing fields, and a series of research achievements have been made [4,5,6].
Yan et al. used Approximate Entropy (APE) [7] to extract bearing fault feature information, and the results showed that this method could effectively extract fault feature, but APE had poor stability and slow calculation speed. Han used sample entropy (SE) [8,9] for bearing fault diagnosis, but this method was easily affected by abrupt signal. Azami proposed dispersion entropy (DE) [10] to quantify the regularity of the time series. This algorithm is fast in calculation and less affected by mutation signals. Rostaghi identified the status of rolling bearing and gear by DE [11], and the results showed that the status representation of rotating machinery by DE was more stable. Subsequently, Azami et al. proposed the fluctuation dispersion entropy (FDE) [12], which considered the fluctuation of time series. FDE only describes the complexity of a nonlinear time series on a single scale, resulting in the loss of a lot of important information. Therefore, Azami et al. proposed multi-scale fluctuation dispersion entropy (MFDE) [13] by extending FDE at multiple scales. To conquer the shortcomings of the coarsening method in the multi-channel signals for MFDE, On the basis of MFDE, GAN et al. introduced the coarse-grained method and proposed the composite multi-scale wave dispersion entropy (CMFDE) [14], proving that the sliding coarse-grained method based on CMFDE has better entropy stability. Zhou et al. proposed that refined complex multi-scale fluctuation dispersion entropy (RCMFDE) [15] is stronger and more stable in extracting features, and RCMDE has a smaller dependence on the length of time series. Azami et al. proposed multi-variable multi-scale dispersion entropy (MMDE) [16] in order to quantify the complexity of multivariate time series. In order to analyze data consisting of more than one-time series, Ahmed et al. proposed multi-variable multi-scale sample entropy (mvMSE) [17]. mvMSE considered both time domain and space domain, reflecting the complexity of multi-channel signals.
In order to solve the shortage of RCMFDE in multivariable time series and the poor stability of mvMSE in feature extraction, in this paper, a refined composite multivariate multiscale fluctuation dispersion entropy is proposed. This method synthesizes the information of multiple coarse-grained sequences in each channel of the multi-variable time series, and uses a refine composite method to make it less dependent on the length of the time series and more stable and reliable in feature extraction. Meanwhile, a fault diagnosis method based on Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy, Joint Mutual Information Maximisation, and Support Vector Machine (RCMMFD-JMIM-SVM) is proposed. In this method, RCMMFDE extracts the multi-variable time-series information of rotating machinery faults, uses JMIM to extract sensitive features, reduce feature dimension, and reduce the time of fault diagnosis. Finally, the SVM classifier is used to identify the fault state of rotating machinery. Two experimental results show that this method can achieve good recognition accuracy.
The rest of this paper is organized as follows. Section 2 introduces the refined complex multi-variable multi-scale fluctuation dispersion entropy method, the correlation characteristic analysis of this method, and introduces the RCMMFDE-JMIM-SVM fault diagnosis method. Section 3 introduces the diagnosis results of RCMMFDE-JMIM-SVM in different rotating machinery fault data sets and discusses relevant results. Finally, the conclusion is drawn in Section 4.

2. Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy

2.1. Multivariate Fluctuating Dispersion Entropy (MFDE)

In order to quantify the complexity of multivariate time series, based on the multivariate embedding theory, the fluctuation dispersion entropy (FDE) was extended to the multivariate fluctuation dispersion entropy (MFDE) [18]. For a multiple time series, X = { x k , i } k = 1 , 2 , , n i = 1 , 2 , , N . The detailed steps of MFDE of X are as follows:
(1)
The multivariate time series X = { x k , i } k = 1 , 2 , , n i = 1 , 2 , , N is mapped to Y = { y k , i } k = 1 , 2 , , n i = 1 , 2 , , N using a normal distribution function
y k , i = 1 σ k 2 π x k , i e t μ k 2 σ k 2 d t
where μ is the expected value and σ 2 is the variance.
(2)
The linear transformation maps Y to Z = { z k , i } k = 1 , 2 , , n i = 1 , 2 , , N . 0 z k , i c
z k , i = i n t ( c · y k , i + 0.5 )
where c is an integer and int(·) represents the rounding function. According to the multi-variable embedding theory, the embedding vector G m ( j ) is calculated by:
G m ( j ) = [ z 1 , j , z 1 , j + τ 1 , , z 1 , j + ( m 1 1 ) τ 1 , z 2 , j , z 2 , j + τ 2 , , z 2 , j + ( m 2 1 ) τ 2 , , z n , j , z n , j + τ n , , z n , j + ( m n 1 ) τ n ]
where j = 1 , 2 , , N ( m 1 ) τ . The embedding dimension vector m = [ m 1 , m 2 , , m n ] , The delay coefficient vector τ = [ τ 1 , τ 2 , , τ n ] . Convert G m ( j ) to F m ( j ) ,
F m ( j ) = { z 1 , j + τ 1 z 1 , j + c , , z 1 , j + ( m 1 1 ) τ 1 z 1 , j + ( m 1 2 ) τ 1 + c , z 2 , j + τ 2 z 2 , j + c , , z 2 , j + ( m 2 1 ) τ 2 z 2 , j + ( m 2 2 ) τ 2 + c , , z n , j + τ n z n , j + c , , z n , j + ( m n 1 ) τ n z n , j + ( m n 2 ) τ n + c , }
For every G m ( j ) , all combinations of m−1 elements in G m ( j ) are called ϕ q , l ( j ) , where q [ 1 , C m 1 ( m 1 ) n ] , l [ 1 , m 1 ] . Map each ϕ q , l ( j ) to a scatter pattern π v 0 v 1 v m 2 ( v = 1 , 2 , , c ) . In this mode, ϕ q , 1 ( j ) = v 0 , ϕ q , 1 ( j ) = v 1 ,…,   ϕ q , l ( j ) = v m 2 . Since π v 0 v 1 v m 2 consists of m−1 digits and each m −1 has class c, there are c ( m 1 ) dispersion patterns. the total number of combinations per G m ( j ) is C m 1 ( m 1 ) n . Therefore, there are [ n ( m 1 ) d ] C m 1 ( m 1 ) n dispersion patterns for time series with n variables. Calculate the probability of each distribution mode π v 0 v 1 v m 2 by:
P ( π v 0 v 1 v m 2 ) = N u m ( π v 0 v 1 v m 2 ) ( N ( m 1 ) d ) C m 1 ( m 1 ) n
(3)
According to the definition of Shannon entropy, the multivariate fluctuation dispersion entropy of the original signal X is:
M F D E ( X , m , c , τ ) = π = 1 ( 2 c 1 ) m 1 P ( π v 0 v 1 v m 2 ) I n P ( π v 0 v 1 v m 2 )

2.2. Multivariate Multiscale Fluctuation Dispersion Entropy (MMFDE)

The steps of the MMFDE are first to establish the coarse-grained time series of the original multivariate signal and then calculate the MFDE for each coarse-grained multivariate time series. For a multiple time series, U = { μ k , b } k = 1 , 2 , , N b = 1 , 2 , , L . The detailed steps of MFDE of U are as follows:
(1)
For a multivariable time series, U = { μ k , b } k = 1 , 2 , , N b = 1 , 2 , , L with length L and number of signal channels N. Then, the coarse granulation time series of scale factor τ is:
x k , i τ = 1 τ b = ( i 1 ) τ + 1 i τ μ k , b ,   1 i L τ
(2)
MMFDE was obtained by calculating the MFDE of each coarse-grained multivariate time series { x k , i τ } under the same parameters. By extending the single-scale MFDE to the multi-scale, more information is obtained from the multi-scale coarse-grained time series of different scales to obtain the multi-scale MMFDE. However, in the coarse-grained multivariate time series whose MMFDE scale factor is τ , only the coarse-grained multivariate time series starting with μ k , 1 is considered, and the remaining τ 1 multivariable time series is missing. The relationship between the coarse-grained time series has not been taken into consideration, resulting in the lack of statistical information.

2.3. Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy (RCMMFDE)

(1)
For a multivariable time series U = { μ k , b } k = 1 , 2 , , N b = 1 , 2 , , L with length L and number of signal channels N, the a coarse-grained time series for the given scale factor τ is X a τ = { x k , i , 1 τ , x k , i , 2 τ , , } , where x k , i , a τ = 1 τ b = a + ( i 1 ) τ a + i τ 1 μ k , b , where 1 i L τ , 1 a τ .
(2)
The RCMMFDE of the original multivariable time series U = { μ k , b } k = 1 , 2 , , N b = 1 , 2 , , L is:
R C M M F D E ( U , m , c , d , τ ) = π ( 2 c 1 ) m 1 p ¯ ( π v 0 v 1 v m 2 ) I n p ¯ ( π v 0 v 1 v m 2 )
where m is the embedding dimension, c is the number of categories, d denotes the delay coefficient, τ represent the scale factor. p ¯ ( π v 0 v 1 v m 2 ) = 1 τ 1 τ p a τ is the average probability of the dispersion mode π v 0 v 1 v m 2 of coarse-grained sequence X a τ . p a τ is the frequency of scattering mode π v 0 v 1 v m 2 in the A multivariate coarse-grained time series X a τ .

2.4. RCMMFDE Feature Analysis

In this section, five sets of simulation signals are used to verify the relevant characteristics of RCMMFDE. Firstly, white noise and 1/f noise are used for analysis to study RCMMFDE’s ability to measure the complexity of multi-variable time series and its stability under different time series lengths. Finally, the synthetic signal is used to verify the anti-noise performance of RCMMFDE.

2.4.1. Mixed Analysis of White Noise and 1/f Noise

In order to illustrate the corresponding characteristics of multivariable time series, RCMMFDE is applied to the generated four-channel time series: (1) all four channels are 1/f noise; (2) three channels are 1/f noise and one channel is white noise; (3) two channels are 1/f noise and two channel is white noise; (4) one channels are 1/f noise and three channel is white noise; (5) all four channels are white noise. The parameter values used to calculate RCMMFDE in this section are set as: embedding dimension m = 3, category c = 5, delay coefficient d = 1, scale factor τ = 2. For the convenience of comparison, the calculated results of the 50 samples were averaged. Figure 1 shows the mean standard deviation curves of RCMMFDE, Refined Composite Multivariate Multiscale Fluctuation Entropy(RCMMFE) [19] and Refined Composite Multivariate Multiscale Sample Entropy(RCMMSE) [20]. As shown in Figure 1, the overall trend of changes in RCMMFDE, RCMMFE and RCMMSE is roughly similar, which all decrease with the increase of white noise channel and the increase of scale factor. This is because the time series consisting of 4-channel 1/f noise is rather complex. The complexity of time series decreases with the decrease of 1/f channel, and the complexity of time series composed of four channel white noise is the lowest. RCMMFE and RCMMSE overlap widely in the first category (all four channels are 1/f noise) and the second category (three channels are 1/f noise and one channel is white noise), and the separation is not obvious. but, RCMMFDE has a clear distinction between category 1 and category 2 with a lower standard deviation. Which proves that RCMMFDE has better resolution ability and stability and more suitable for the multi-channel signal feature extraction.

2.4.2. Analysis of Anti-Noise Performance

In order to study the anti-noise performance of different entropy algorithms, a two-channel composite signal with noise power variation is adopted. Each channel has a sampling frequency of 512Hz and a length of 32s. To make a fair comparison, the embedded dimensions in RCMMFDE, RCMMFE, RCMMSE and RCMMPE are set to equal m = 3, and r = 0.15 in RCMMFE and RCMMSE set the fuzzy power of RCMMFE to n = 2. The sensitivity of RCMMFDE to noise was studied by using the dual-channel periodic signal with noise power variation. This periodic signal had no noise in the first 4 s, and white noise was added after 4 s, and the noise power was increased every 1 s. In order to reduce the influence of entropy on the visualization of entropy change, the entropy values of different methods are scaled to the range of 0–1, as shown in Figure 2. As shown in Figure 2, the values of RCMMSE and RCMMPE increase rapidly with the increase of noise power, but RCMMPE increases faster than RCMMSE. The curve of RCMMFDE grows slowly, but RCMMFE performs better in variable noise synthesis signals. The results show that RCMMFE and RCMMFDE have better anti-noise capability.

2.4.3. RCMMFDE Data Length Sensitivity Analysis

In order to compare the stability of RCMMFDE under different data lengths, four-way white noise samples with data lengths of 512, 1024, 1536, 2048, 2560 and 3072 were used, and the number of samples of each data length was 100. RCMMFDE was calculated 20 times for each data length, and its mean value and standard deviation are shown in Figure 3. As shown in Figure 3a, as the data length increases, the value of RCMMFDE is relatively close. However, when the scale factor is greater than 5, the RCMMFDE with lower data length is slightly smaller than the RCMMFDE with higher data length, and this gap will increase with the increase of the scale factor. This indicates that RCMMFDE is less sensitive to data length and can obtain reliable measurements of samples of different data lengths. As shown in Figure 3b, the standard deviation of RCMMFDE increases with the increase of scale factor and the decrease of data length. In particular, when the data length is 500, the standard deviation increases significantly. When the data length is 3000, the standard deviation of RCMMFDE is the smallest, which indicates that the longer the data length is, the more stable RCMMFDE is. However, the data length is positively correlated with the computation time. Given the tradeoff between RCMMFDE stability and computation time, the data sample length is set to 2048.

2.5. Fault Diagnosis Based on RCMMFDE

2.5.1. JMIM Feature Selection

Joint Mutual Information Maximisation (JMIM) [21] is an effective feature selection algorithm. Considering the overall stability of Joint Mutual Information, JMIM is able to select a subset of features with the same or even better effect from the original feature set [22]. Compared with other algorithms, such as mRMR, JMI and IG, the feature subset selected by JMIM algorithm has better classification performance [23]. Therefore, the JMIM algorithm is used as the sensitive sub-feature selection method in this paper.

2.5.2. RCMMFDE-JMIM-SVM Fault Diagnosis Algorithm

In many fault diagnosis fields, Support Vector Machine (SVM) has been proved to be robust nonlinear multiple classifiers [24,25]. Combined with the advantages of RCMMFDE, JMIM, and SVM, this paper proposes a fault diagnosis method of rotating machinery based on RCMMFDE-JMIM-SVM. The flow chart is shown in Figure 4.
Step 1.
Signal acquisition: Through acceleration sensors installed at different positions of the rotating machinery to be monitored, multi-channel vibration signals of the rotating machinery under different speeds are collected. The data samples collected under different running speeds can form a fault sample data set.
Step 2.
Feature extraction: RCMMFDE algorithm is used to extract features from fault sample data to form feature data sets, which are divided into training data sets and test data sets.
Step 3.
Feature selection: JMIM algorithm is used to select the most sensitive sub-features of the feature.
Step 4.
SVM model training: SVM classifier is trained using sensitive features. The SVM classifier includes M-1 dichotomous SVM, where M is the type of fault sample.
Step 5.
Fault diagnosis: Input the samples of the test data set into the RCMMFDE-JMIM-SVM classifier to identify the fault types of different bearing samples.

3. Experimental Verification and Analysis

In this paper, a gearbox fault data set and a rolling bearing fault data set are respectively used to verify the effectiveness of the method. Collect the multivariable time series of gearbox and rolling bearing at constant speed. RCMMFE, RCMMSE, RCMMPE and RCMMFDE were compared, and the good performance of the method was verified by the identification accuracy and CPU time.

3.1. Validation of Gearbox Fault Data Set

3.1.1. Data Set Description

This data set is the 2009 PHM Challenge data set (in Supplementary Materials), which is the gearbox composite failure test data obtained on the gearbox test bench. The test platform is shown in Figure 5a, which mainly includes induction motor, gearbox, speed measuring device, solenoid brake, AC drive, experimental bench and data acquisition system. The gearbox includes input shaft, intermediate shaft and output shaft. The internal structure diagram is shown in Figure 5b. Two accelerometers are placed at the input end and output end of the gearbox respectively to collect vibration signals of the gear box, as shown in Figure 5c,d. This paper takes the meshing fault of helical gear as the research object, Based on Tooth Surface Crack (TSC), Tooth Root Fracture (TRF), Bearing Outer Fault (BOF), Bearing Inner Outer Fault (BIF), Bearing Balling Fault, As for fault types such as BBF), Gear Imbalance (GI) and Input Shaft Imbalance (ISI), a result of eight helical gears meshing composite fault states is simulated. Information about the composite fault state is illustrated in Table 1, and Gear of each fault is illustrated in Figure 6. The sampling frequency is 66.67 KHz, the rated speed of input shaft is 3000 r/min, and the sampling time is 4 s gear box vibration signal under each compound fault state. The signal is divided into several samples, each with a length of 2048. The labels and sample Numbers are shown in Table 2. The time-domain waveform of the gear in different states is shown in Figure 7. Because there are many fault states, it is difficult to distinguish the specific bearing fault from the time domain waveform.

3.1.2. RCMMFDE-JMIM-SVM Sensitive Feature Number Analysis

In each data sample of the above data set, 200 samples were randomly selected as test objects, of which 150 were training samples and 50 were test samples. The sensitive feature dimension p extracted by JMIM was set as 25. In order to verify the performance of the RCMMFDE-JMIM-SVM method, RCMMFDE, RCMMFE, RCMMSE, and RCMMPE were used to extract the features of vibration signals, and the optimal features were extracted through the JMIM algorithm. Finally, the optimal features were input into the SVM classifier to identify the fault type of the test sample. Each method was run 20 times respectively. The specific parameters of RCMMFDE can be known from Azami research that the number of categories c = 6, the delay coefficient d = 1, the embedding dimension m = 2, and the maximum scale factor τ = 25 [26]. The average diagnostic accuracy of different methods is shown in Figure 8. As can be seen from Figure 8, the recognition accuracy of RCMMFDE-JMIM-SVM method is better than other methods, because RCMMFDE takes into account the amplitude changes between two adjacent elements at different scales of each channel in the multi-channel time series. Secondly, too large or too small number of sensitive features will reduce the accuracy of fault recognition. This is because too small number of feature sensitivities will contain less fault feature information. On the contrary, too large number of sensitive features will lead to redundancy of fault feature information and reduce fault identification accuracy. Secondly, too large or too small several sensitive features will reduce the accuracy of fault recognition. When the number of sensitive features is less than 21, the accuracy of the RCMMFDE-JMIM-SVM method will increase with the increase of the number of sensitive features and eventually tend to be stable. When the number of sensitive features is greater than 21, the accuracy of the model will decrease with the increase in the number of sensitive features. Therefore, this paper sets the number of sensitive features at 20.

3.1.3. RCMMFDE-JMIM-SVM Performance Analysis at Constant Speed

To verify the effectiveness of the proposed method under the condition of constant rotation speed, RCMMFDE-JMIM-SVM, RCMMFE-JMIM-SVM, RCMMSE-JMIM-SVM, and RCMMPE-JMIM-SVM methods were used to analyze the gear fault data set. The main parameters of these methods include embedding dimension m, number of categories c, delay coefficient d, fuzzy power f, scale factor τ , threshold r and characteristic sensitive number p. For the sake of fair comparison, the scale factors of these methods are consistent, and the specific parameter settings are shown in Table 3.
Figure 9 and Table 4 show the identification accuracy of four different fault diagnosis methods within 30 running times. The computers used in this paper are configured as Windows 10 operating system, Intel Core i5-9400f CPU and 16 G RAM. As is shown in Figure 9, The average accuracy of RCMMFDE-JMIM-SVM is 98.757%, and the maximum classification accuracy is 99.252%. This shows that the method is effective for the diagnosis of gear fault data set. The classification accuracy of RCMMFE-JMIM-SVM is between 97.867% and 97.305%. Compared with the method presented in this paper, the diagnostic accuracy of RCMMFE-JMIM-SVM is slightly lower. The classification accuracy of RCMMPE-JMIM-SVM method was the lowest, between 93.217% and 94.045%. Furthermore, the processing time of RCMMFDE-JMIM-SVM for a single sample was 10.32s, 0.53s slower than that of RCMMSE-SVM, which was nearly 90 times slower than that of RCMMFDE-JMIM-SVM. The reason why RCMMFDE-JMIM-SVM achieves excellent performance is that the fluctuation between adjacent elements of each channel of the multivariate time series is considered, and the features can be extracted from the multivariate time series between different signals, and the influence of noise and fluctuation can be resisted.
To analyze the relationship between the recognition accuracy and the number of training samples, we randomly selected training samples from the total samples under different proportions (10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%) and calculated the average recognition accuracy of different diagnostic methods in the running time of 50 times respectively. Figure 10 shows the average classification accuracy under different proportions of training samples. As shown in Figure 10, the accuracy of different classification methods increases with the proportion of training samples. Although the proportions of training samples are different, RCMMFDE-JMIM-SVM can achieve the highest diagnostic accuracy. Even if 10% of the samples are used for training, the diagnostic accuracy can reach 97.896%, and when the proportion of training samples increases to 30%, the accuracy reaches 98.961%. It can be found that if sufficient training data samples are available, all the four methods can achieve high average diagnostic accuracy. However, too many training samples will lead to a significant increase in training time. Therefore, to balance the training time and accuracy, the proportion of training samples was set at 40%.

3.1.4. RCMMFDE-JMIM-SVM Performance Verification

To better present the performance of RCMMFDE-JMIM-SVM, this paper also uses the common feature of vibration signal and the SVM (CFVS-SVM) to verify performance of proposed method. Common features we take in this paper are time-domain features, frequency-domain features, and time-frequency features. We take 18 features including mean value, standard deviation, peak-to-peak value, median, interquartile range and etc. as the time domain features. The frequency-domain features are median mean, interquartile range, percentile deviation, and Fourier transform features. The time-frequency domain features are the six energy features and their normalized energy features obtained by the five-stage wavelet transform.
Table 5 shows the results of four entropy fault diagnosis methods and CFVS-SVM diagnosis methods trained with 40% training samples. The classification accuracy of Table 5 is the average value of 20 times’ result. The Table 5 shows that in G2, G4 and G5 and G8 condition, CFVS-SVM and RCMMPE-JMIM-SVM accuracy have plunged, the accuracy of the RCMMSE-JMIM-SVM and RCMMFE-JMIM-SVM have appeared in the different amplitude of falling, but in this paper, fault identification accuracy of the presented method still remain more than 99%, fully shows in different fault conditions, the proposed method of feature extraction is better than the other three methods.
In order to demonstrate the effectiveness of RCMMFDE-JMIM feature extraction, Figure 11 shows the feature distribution visualization diagram of all samples extracted by RCMMFDE-JMIM algorithm by T-SNE method [27]. It can be seen from Figure 11 that each fault sample has a high degree of intra-class clustering and inter-class separation, but the distance between classes is small in some states of G3 and G4, which leads to the relatively low accuracy of this method in G3 state, but it does not affect the performance of RCMMFDE-JMIM-SVM method. In order to demonstrate the effectiveness of RCMMFDE-JMIM feature extraction, Figure 11 shows the feature distribution visualization diagram of all samples extracted by the RCMMFDE-JMIM algorithm by the T-SNE method [28]. It can be seen from Figure 11 that each fault sample has a high degree of intra-class clustering and inter-class separation, but the distance between classes is small in some states of G3 and G4, which leads to the relatively low accuracy of this method in the G3 state, but it does not affect the performance of RCMMFDE-JMIM-SVM method.

3.2. Verification of Rolling Bearing Fault Data

3.2.1. Description of Rolling Bearing Fault Data Set

The test platform is the threshing cylinder part of the combine harvester, which is mainly composed of the motor driving part, threshing assembly, and data acquisition system. The test bearing is the 6307 deep-groove ball bearing at the back end of the threshing roller. The outer ring of the bearing is fixed on the bearing pedestal, and the inner ring rotates with the threshing roller shaft. Other bearings of the system are normal bearings. In the test, three common fault types of 6307 deep-groove ball bearings were considered, including Inner Ring Fault (IRF), Outer Ring Fault (ORF), and Composite Fault (CF) occurring simultaneously in the inner ring. The detailed parameters of each fault are shown in Table 6. Bearing faults of different types and sizes are simulated by EDM pitting technique. Figure 12 shown rolling bearings with combined faults of inner ring, outer ring and inner ring, respectively. The acceleration sensor is fixed at the measuring point through the magnetic base. The voltage sensitivity of the 1# to 4# acceleration sensors is 101.6 mV/g, 99.1 mV/g, 101.2 mV/g, and 101.3 mV/g, respectively. The installation locations of acceleration sensors and speed sensors are shown in Figure 13. The threshing cylinder rotation speed was set at 80 r/min, 160 r/min, 240 r/min, and 320 r/min to collect the bearing vibration signal through the data acquisition system, the sampling frequency is 5120 Hz, each group of sampling time 50 s, The sampling number of each sample is 2048, each bearing status get 500 samples, a total of 5500 samples. The time domain waveform of the original four-channel signals of different fault types is shown in Figure 14.

3.2.2. Rolling Bearing Fault Data Analysis

The 100 samples of each bearing state are taken as training data, and the other 25 samples are taken as test data, thus forming the rolling bearing fault data set. In the rolling bearing fault data analysis, the RCMMFDE-JMIM-SVM method parameters determined in 3.1 are used. The comparison results of RCMMFDE-JMIM-SVM, RCMMFE-JMIM-SVM, RCMMSE-JMIM-SVM, RCMMPE SVM and CFVS-SVM are shown in Table 7. As can be seen from Table 7, the accuracy of all methods on the rolling bearing data set is significantly higher than that on the gear data set. The accuracy of RCMMFDE-JMIM-SVM and RCMMFE-JMIM-SVM reaches 100%, and the accuracy of RCMMSE-JMIM-SVM RCMMPE-JMIM-SVM and CFVS-SVM is also greater than 96%. The reasons for this phenomenon may be the relatively ideal test environment and the good performance of the data acquisition system, which makes the characteristics of the signal data more easily to identify. It also shows that RCMMFDE-JMIM-SVM has good performance.
To demonstrate the effectiveness of RCMMFDE-JMIM-SVM in the bearing fault data set, Figure 15 shows the visualization of the feature distribution of all samples by t-SNE method. It can be seen from Figure 15 that the separation between classes of each fault feature is high and the separation is also high at different speeds of different faults.

4. Conclusions

In this paper, a new complexity assessment method of multivariable nonlinear time series based on RCMMFDE is proposed. Compared with RCMMFE, RCMMSE, and RCMMPE, RCMMFDE has a better stability. This method takes into account the information of different scale sequences of each channel in the original multicomponent signal and obtains more stable and reliable characteristics by using the fine composite analysis technique. Use the JMIM method to select the more sensitive feature. In the classification step, the SVM classifier is trained to classify sensitive feature samples. This method can effectively process multivariate time series data.
The simulation signals with different combinations of WGN and 1/f show that RCMMFDE has better resolution and stability. The research on the two-channel stimulation signal of noise power variation shows that RCMMFDE has excellent anti-noise performance. Meanwhile, the use of four-channel WGN signal indicates that RCMMFDE is not sensitive to the length of time series signal. Taking gear fault data set and bearing fault data set as examples, the effectiveness of RCMMFDE-JMIM-SVM method is verified. The results show that this method has the advantages of high classification accuracy and short calculation time.
The research in this paper shows that the RCMMFDE-JMIM-SVM method is suitable for the diagnosis of gear and bearing faults, and the effectiveness of this method on other data sets can be further studied.

Supplementary Materials

The following are available online at http://www.phmsociety.org/competition/09, Figure: Test diagram of gear failure, Data set: generic industrial gearbox data.

Author Contributions

Data curation, C.X. and L.L.; Formal analysis, C.X.; Investigation, G.Y.; Project administration, G.Y.; Resources, G.Y.; Software, C.X.; Supervision, G.Y.; Writing—original draft, C.X.; Writing—review & editing, H.J. and X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by National Key Research and Development Program of China (No. 2018YFB0105300 and No. 2017YFD07006).

Data Availability Statement

The datasets used or analysed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, T. Rolling-bearings fault diagnosis based-on empirical mode decomposition and least square support vector machine. J. Mech. Eng. 2007, 43, 88–92. [Google Scholar] [CrossRef]
  2. Duan, Z.; Wu, T.; Guo, S.; Shao, T.; Malekian, R.; Li, Z. Development and trend of condition monitoring and fault diagnosis of multi-sensors information fusion for rolling bearings: A review. Int. J. Adv. Manuf. Technol. 2018, 96, 803–819. [Google Scholar] [CrossRef] [Green Version]
  3. De Azevedo, H.D.M.; Araújo, A.M.; Bouchonneau, N. A review of wind turbine bearing condition monitoring: State of the art and challenges. Renew. Sustain. Energy Rev. 2016, 56, 368–379. [Google Scholar] [CrossRef]
  4. Zheng, J.; Huang, S.; Pan, H.; Jiang, K. An improved empirical wavelet transform and refined composite multiscale dispersion entropy based fault diagnosis method for rolling bearing. IEEE Access 2019, 8, 168732–168742. [Google Scholar] [CrossRef]
  5. Coifman, R.R.; Wickerhauser, M.V. Entropy-based algorithms for best basis selection. IEEE Trans. Inf. Theory 2002, 38, 713–718. [Google Scholar] [CrossRef] [Green Version]
  6. Porta, A.; Bari, V.; De Maria, B.; Cairo, B.; Vaini, E.; Malacarne, M.; Pagani, M.; Lucini, D. On the Relevance of Computing a Local Version of Sample Entropy in Cardiovascular Control Analysis. IEEE Trans. Biomed. Eng. 2019, 3, 623–631. [Google Scholar] [CrossRef]
  7. 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]
  8. 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]
  9. Richman, J.S.; Lake, D.E.; Moorman, J.R. Sample Entropy. Method Enzym. 2004, 384, 172–184. [Google Scholar]
  10. Rostaghi, M.; Azami, H. Dispersion Entropy: A Measure for Time-Series Analysis. IEEE Signal Process. Lett. 2016, 23, 610–614. [Google Scholar] [CrossRef]
  11. Rostaghi, M.; Ashory, M.R.; Azami, H. Application of dispersion entropy to status characterization of rotary machines. J. Sound Vib. 2019, 438, 291–308. [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. Azami, H.; Arnold, S.E.; Sanei, S.; Chang, Z.; Sapiro, G.; Escudero, J.; Gupta, A.S. Multiscale Fluctuation-based Dispersion Entropy and its Applications to Neurological Diseases. IEEE Access 2019, 7, 68718–68733. [Google Scholar] [CrossRef]
  14. 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] [PubMed] [Green Version]
  15. Zhou, F.; Yang, X.; Shen, J.; Liu, W. Fault Diagnosis of Hydraulic Pumps Using PSO-VMD and Refined Composite Multiscale Fluctuation Dispersion Entropy. Shock Vib. 2020, 2020, 1–13. [Google Scholar] [CrossRef]
  16. Azami, H.; Fernandez, A.; Escudero, J. Multivariate Multiscale Dispersion Entropy of Biomedical Times Series. Entropy 2017, 21, 913. [Google Scholar] [CrossRef] [Green Version]
  17. Yang, Y.; Zheng, H.; Yin, J.; Xu, M.; Chen, Y. Refined composite multivariate multiscale symbolic dynamic entropy and its application to fault diagnosis of rotating machine. Measurement 2020, 151, 107233. [Google Scholar] [CrossRef]
  18. Ahmed, M.U.; Mandic, D.P. Multivariate Multiscale Entropy Analysis. IEEE Signal Process. Lett. 2012, 19, 91–94. [Google Scholar] [CrossRef] [Green Version]
  19. Zheng, J.; Tu, D.; Pan, H.; Hu, X.; Liu, T.; Liu, Q. A Refined Composite Multivariate Multiscale Fuzzy Entropy and Laplacian Score-Based Fault Diagnosis Method for Rolling Bearings. Entropy 2017, 19, 585. [Google Scholar] [CrossRef] [Green Version]
  20. Azami, H.; Escudero, J.; Fernández, A. Refined composite multivariate multiscale entropy based on variance for analysis of resting-state magnetoencephalograms in Alzheimer’s disease. In Proceedings of the 2016 International Conference for Students on Applied Engineering (ICSAE), Newcastle upon Tyne, UK, 20–21 October 2016; pp. 413–418. [Google Scholar] [CrossRef] [Green Version]
  21. Bennasar, M.; Hicks, Y.; Setchi, R. Feature selection using Joint Mutual Information Maximisation. Expert Syst. Appl. 2015, 42, 8520–8532. [Google Scholar] [CrossRef] [Green Version]
  22. Mao, Y.; Cao, H.; Ping, P.; Li, X. Feature selection based on maximum conditional and joint mutual information. J. Comput. Appl. 2019, 39, 734–741. [Google Scholar]
  23. Rajlakshmi, S.; Sane, S. JMIM: A Feature Selection Technique using Joint Mutual Information Maximization Approach. Int. J. Comput. Appl. 2017, 975, 8887. [Google Scholar]
  24. Dao, T.K.; Nguyen, T.T.; Pan, J.S.; Qiao, Y.; Lai, Q.A. Identification Failure Data for Cluster Heads Aggregation in WSN based on Improving Classification of SVM. IEEE Access 2020, 8, 61070–61084. [Google Scholar] [CrossRef]
  25. Venkatesan, C.; Karthigaikumar, P.; Paul, A.; Satheeskumaran, S.; Kumar, R.J. ECG Signal Preprocessing and SVM Classifier-Based Abnormality Detection in Remote Healthcare Applications. IEEE Access 2018, 6, 9767–9773. [Google Scholar] [CrossRef]
  26. Mao, X.; Wang, L.; Li, C. SVM Classifier for Analog Fault Diagnosis Using Fractal Features. In Proceedings of the Second International Symposium on Intelligent Information Technology Application (IITA’08), Shanghai, China, 20–22 December 2008. [Google Scholar]
  27. Azami, H.; Rostaghi, M.; Escudero, J. Refined Composite Multiscale Dispersion Entropy: A Fast Measure of Complexity. arXiv 2016, arXiv:1606.01379. [Google Scholar]
  28. Laurens, V.D.M.; Hinton, G. Visualizing Data using t-SNE. J. Mach. Learn. Res. 2008, 9, 2579–2605. [Google Scholar]
Figure 1. RCMMFDE, RCMMFE and RCMMSE values of different methods under different noise combinations: (a) RCMMFDE; (b) RCMMFE; (c) RCMMSE.
Figure 1. RCMMFDE, RCMMFE and RCMMSE values of different methods under different noise combinations: (a) RCMMFDE; (b) RCMMFE; (c) RCMMSE.
Entropy 23 00128 g001
Figure 2. RCMMFDE of containing periodic signals of varying noise.
Figure 2. RCMMFDE of containing periodic signals of varying noise.
Entropy 23 00128 g002
Figure 3. Comparison of white noise with different data lengths: (a) Mean value of RCMMFDE; (b) Standard deviation of RCMMFDE.
Figure 3. Comparison of white noise with different data lengths: (a) Mean value of RCMMFDE; (b) Standard deviation of RCMMFDE.
Entropy 23 00128 g003
Figure 4. RCMMFDE-JMIM-SVM fault diagnosis flow chart.
Figure 4. RCMMFDE-JMIM-SVM fault diagnosis flow chart.
Entropy 23 00128 g004
Figure 5. Gearbox fault test platform: (a) gearbox test bench; (b) internal structure of gearbox; (c) installation position of accelerometer 1; (d) installation position of accelerometer 2.
Figure 5. Gearbox fault test platform: (a) gearbox test bench; (b) internal structure of gearbox; (c) installation position of accelerometer 1; (d) installation position of accelerometer 2.
Entropy 23 00128 g005
Figure 6. Gear fault type: (a) Normal; (b) TRF; (c) TSC.
Figure 6. Gear fault type: (a) Normal; (b) TRF; (c) TSC.
Entropy 23 00128 g006
Figure 7. Time domain waveform of vibration signal: (a) G1; (b) G2; (c) G3; (d) G4; (e) G5; (f) G6; (g) G7; (h) G8.
Figure 7. Time domain waveform of vibration signal: (a) G1; (b) G2; (c) G3; (d) G4; (e) G5; (f) G6; (g) G7; (h) G8.
Entropy 23 00128 g007aEntropy 23 00128 g007b
Figure 8. Average diagnostic accuracy under different number of sensitive features.
Figure 8. Average diagnostic accuracy under different number of sensitive features.
Entropy 23 00128 g008
Figure 9. Diagnostic accuracy of different diagnostic methods.
Figure 9. Diagnostic accuracy of different diagnostic methods.
Entropy 23 00128 g009
Figure 10. Recognition accuracy of different training sample proportions.
Figure 10. Recognition accuracy of different training sample proportions.
Entropy 23 00128 g010
Figure 11. Visualized by T-SNE features.
Figure 11. Visualized by T-SNE features.
Entropy 23 00128 g011
Figure 12. Different types of faulty bearings.
Figure 12. Different types of faulty bearings.
Entropy 23 00128 g012
Figure 13. Schematic diagram of experimental platform and sensor installation: (a) threshing cylinder front view; (b) threshing cylinder back view; (c) sensor mounting position 1; (d) sensor mounting position 2.
Figure 13. Schematic diagram of experimental platform and sensor installation: (a) threshing cylinder front view; (b) threshing cylinder back view; (c) sensor mounting position 1; (d) sensor mounting position 2.
Entropy 23 00128 g013
Figure 14. Time domain waveform of four different states: (a) normal; (b) IRF15; (c) ORF15; (d) CD10.
Figure 14. Time domain waveform of four different states: (a) normal; (b) IRF15; (c) ORF15; (d) CD10.
Entropy 23 00128 g014
Figure 15. Visualized features by T-SNE.
Figure 15. Visualized features by T-SNE.
Entropy 23 00128 g015
Table 1. Compound fault status information of helical gear meshing.
Table 1. Compound fault status information of helical gear meshing.
AbbreviationG1G2G3G4G5G6G7G8
ShaftInputNormalNormalNormalNormalNormalISINormalISI
OutputNormalNormalNormalNormalNormalNormalKeyway wearNormal
BearingInput shaft—Input endNormalNormalNormalBBFBIFBIFBIFNormal
Intermediate shaft—
Input shaft
NormalNormalNormalNormalBBFBBFNormalBBF
Output shaft—Input endNormalNormalNormalNormalBOFBOFNormalBOF
Input shaft—Output endNormalNormalNormalNormalNormalNormalNormalNormal
Intermediate shaft—
Output end
NormalNormalNormalNormalNormalNormalNormalNormal
Output shaft—Output endNormalNormalNormalNormalNormalNormalNormalNormal
Gear32TNormalTSCNormalNormalTSCNormalNormalNormal
96TNormalNormalNormalNormalNormalNormalNormalNormal
48TNormalGIGIGIGINormalNormalNormal
80TNormalNormalNormalTRFTRFTRFNormalNormal
Table 2. Gearbox fault data set label.
Table 2. Gearbox fault data set label.
NoStateTraining SamplesTesting SamplesClass Label
1G12001001,0,0,0,0,0,0,0
2G22001000,1,0,0,0,0,0,0
3G32001000,0,1,0,0,0,0,0
4G42001000,0,0,1,0,0,0,0
5G52001000,0,0,0,1,0,0,0
6G62001000,0,0,0,0,1,0,0
7G72001000,0,0,0,0,0,1,0
8G82001000,0,0,0,0,0,0,1
Table 3. Parameter setting table.
Table 3. Parameter setting table.
MethodParameter
RCMMFE-JMIM-SVM p = 17m = 2r = 0.15d = 1 τ = 25f = 2
RCMMSE-JMIM-SVM p = 12m = 2r = 0.15d = 1 τ = 25
RCMMPE-JMIM-SVM p = 17m = 2 τ = 17
RCMMFDE-JMIM-SVM p = 20m = 2c = 6d = 1 τ = 25
Table 4. Diagnostic accuracy of different diagnostic methods.
Table 4. Diagnostic accuracy of different diagnostic methods.
MethodMax AccuracyMin AccuracyAverage AccuracyRun Time
RCMMFDE-JMIM-SVM99.25298.51898.7571.67
RCMMSE-JMIM-SVM96.88296.21596.6481.14
RCMMPE-JMIM-SVM94.04593.21793.554148.79
RCMMFE-JMIM-SVM97.86797.30597.6948.28
Table 5. Different methods each fault state accuracy rate.
Table 5. Different methods each fault state accuracy rate.
StateTesting SamplesAccuracy/%
RCMMFDE-JMIM-SVMRCMMFE-JMIM-SVMRCMMSE-JMIM-SVMRCMMPE-JMIM-SVMCFVS-SVM
G118099.42598.62798.12996.95197.273
G218099.35497.05697.28191.48894.445
G318097.81295.62798.83597.53897.758
G418099.48997.56695.75692.15795.561
G518099.10597.63496.01892.26693.891
G618099.28198.46196.48594.13796.659
G718099.57797.34397.52894.51396.562
G818098.90796.93294.41493.06494.827
Total90099.17397.78296.61993.75194.658
Table 6. Table of parameters of different fault state of bearing.
Table 6. Table of parameters of different fault state of bearing.
Bearing StateAbbreviationFault Width (mm)Fault Depth (mm)Threshing Cylinder Speed (r/min)
NormalNormal0080/160/240/320
Outer Ring Fault 1IRF070.73.780/160/240/320
Outer Ring Fault 2IRF101.03.780/160/240/320
Outer Ring Fault 3IRF121.23.780/160/240/320
Outer Ring Fault 4IRF151.53.780/160/240/320
Outer Ring Fault 1ORF070.73.280/160/240/320
Outer Ring Fault 2ORF101.03.280/160/240/320
Outer Ring Fault 3ORF121.23.280/160/240/320
Outer Ring Fault4ORF151.53.280/160/240/320
Composite Fault 1CF050.5/0.53.2/3.780/160/240/320
Composite Fault 2CF101.0/1.03.2/3.780/160/240/320
Table 7. The results of different diagnostic methods were compared.
Table 7. The results of different diagnostic methods were compared.
Bearing StateTesting SamplesRCMMFDE-JMIM-SVMRCMMFE-JMIM-SVMRCMMSE-JMIM-SVMRCMMPE-JMIM-SVMCFVS-SVM
MCRAcc/%MCRAcc/%MCRAcc/%MCRAcc/%MCRAcc/%
Normal250100.000100.000100.000100.00CF0596.00
IRF07250100.000100.00IRF1096.00IRF1288.00IRF1092.00
IRF10250100.000100.000100.00IRF0796.000100.00
IRF12250100.000100.000100.000100.00IRF1096.00
IRF15250100.000100.00IRF1296.000100.000100.00
ORF07250100.000100.000100.000100.00ORF1088.00
ORF10250100.000100.000100.00ORF1596.000100.00
ORF12250100.000100.00ORF1292.00ORF1096.000100.00
ORF15250100.000100.000100.000100.00ORF1296.00
CF05250100.000100.000100.00CF1096.000100.00
CF10250100.000100.000100.000100.00CF0596.00
Total2750100.000100.00499.27797.46996.73
Note: MCR represent the misclassification result; Acc represent the Accuracy.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xi, C.; Yang, G.; Liu, L.; Jiang, H.; Chen, X. A Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy and Its Application to Multivariate Signal of Rotating Machinery. Entropy 2021, 23, 128. https://0-doi-org.brum.beds.ac.uk/10.3390/e23010128

AMA Style

Xi C, Yang G, Liu L, Jiang H, Chen X. A Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy and Its Application to Multivariate Signal of Rotating Machinery. Entropy. 2021; 23(1):128. https://0-doi-org.brum.beds.ac.uk/10.3390/e23010128

Chicago/Turabian Style

Xi, Chenbo, Guangyou Yang, Lang Liu, Hongyuan Jiang, and Xuehai Chen. 2021. "A Refined Composite Multivariate Multiscale Fluctuation Dispersion Entropy and Its Application to Multivariate Signal of Rotating Machinery" Entropy 23, no. 1: 128. https://0-doi-org.brum.beds.ac.uk/10.3390/e23010128

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