Abstract
Advanced research studies on industrial Internet of things require effective feature extraction and accurate machinery health state evaluation. For roller bearing, a well-known mechanical component most extensively used in the industry, its running status directly affects the operation of the entire machinery and equipment. For intelligent fault diagnosis of roller bearing, the selection of the intrinsic mode function (IMF) modes in approaches of ensemble empirical mode decomposition (EEMD)/variational mode decomposition (VMD) becomes a tricky problem. To solve this problem, this study proposed an efficient scheme on roller bearing fault diagnosis that combines the refined composite multivariate multiscale sample entropy (RCMMSE) with different classifiers. Firstly, the synthetic noise signals are introduced to compare the effectiveness of the multiscale sample entropy (MSE) and the RCMMSE models. Secondly, the random noise signals are used to compare the performance of EEMD and VMD methods, where the envelope spectrum characteristics of fault signals are well described. Moreover, EEMD/VMD methods are utilized to decompose the roller bearing vibration signals into various modes to get the entropy values. Finally, the obtained RCMMSE is adopted as a feature vector and subsequently employed as an input of support vector machine, random forest, and probabilistic neural network models to conduct roller bearing fault identification. The extensive experimental results prove that this proposed scheme performs well and the classification accuracy of VMD-RCMMSE is higher than EEMD-RCMMSE.
1. Introduction
With the rapid development of Internet of things (IoT) [1, 2] and Industry 4.0 [3], there are increasingly massive real-time data from various types of mechanical equipment [4]. The availability of these data that contain abundant information about machine health has attracted more and more enterprises’ attention. It has been proved that large volume, high velocity, and diversity mechanical big data are the major properties of mechanical big data [5, 6]. Effective feature extraction from these data and accurate machinery health state evaluation with ever-accelerated updating of schemes have become hot research issues in the prognostic and health management systems in the era of industrial IoT [7].
For roller bearing, a well-known mechanical component most extensively used in the industry, its operating status directly affects the operation of the entire machinery and equipment. Roller bearing failure is an important factor leading to mechanical equipment failure, so timely detection and fault diagnosis are of great significance, and analysing the vibration signal collected by the sensors to determine its failure is a commonly adopted scheme [8]. As the roller bearings are usually working in the vibration source environment, with the destruction of complex forms, the vibration signal is typically nonlinear and nonstationary. How to excavate the fault feature based on the vibration signals has been a research hotspot [9, 10].
The widely used signal decomposing schemes include wavelet analysis [11, 12] and EMD [13]. However, the wavelet transformation is essentially nonadaptive due to the configuration of wavelet basis and decomposing layers. In EMD, the complicated signal can be self-adaptively decomposed into some IMFs with a residual component. Nevertheless, mode mixing is a stumbling block in EMD. To deal with this challenge, EEMD is proposed [14], in which a complex signal can be disintegrated into IMFs in terms of the local time-scale characteristic of the signal. In recent, EEMD has been extensively employed in fault detection [15].
Using the EMD or EEMD schemes for fault feature extraction has been widely concerned. However, these two schemes have some disadvantages. The recursive mode decomposition in EMD/EEMD will propagate the envelope estimation error continuously. The signal contains no noise or intermittent signal, which leads to the decomposition of the mode mixing. Although the white noise scheme was used to suppress the mode mixing, the scheme needs to be carried out several hundreds of times of EMD/EEMD operation and will break out the signal composed of more than a real component. Moreover, EMD and EEMD cannot be separated correctly from the close frequencies, and it is quite challenging to select the suitable number of IMFs when applying EMD and EEMD.
To address the abovementioned problems, a neoteric scheme of signal decomposition estimation called VMD was proposed [16–18]. The whole framework of the alternating direction multiplier scheme was resorted to consecutively refresh the modes and their centre frequencies and gently demodulate the modes into the corresponding base frequency bands. Ultimately, each mode and corresponding centre frequencies were abstracted together. In comparison with the recursive filter pattern of EMD and EEMD, VMD converts the signal into variational and non-recursive decomposition patterns, and it consists of adaptive Wiener filter groups in nature. VMD can isolate two pure harmonic signals with similar frequencies [16–18]. The IMFs obtained by EEMD denote the natural oscillatory pattern inserted into the signals, where the entropy values of every IMF were often abstracted as a feature to discover the properties of the vibration signals [19, 20].
Entropy was one of the schemes to evaluate the time-series complexity. S. M. Pincus proposed the approximate entropy (AE) [21]. However, the length of the time series makes a critical influence on the performance of AE scheme. Hence, the value of AE is conformably lower than the intended one and fails relative coherence particularly when the data length is short. Based on AE, SE was proposed to deal with this challenge [22]. It has the advantages of short data, stable and low noise and interference capacity, and good consistency in the region of large parameter range. It can be noticed that the irregularity of time series can be only reflected by AE and SE individually. When the roller bearings fail, not only the frequencies but also the corresponding complexity of vibration signals has tremendous deviations. Accordingly, ME can be viewed as a property index for fault diagnosis [23]. Considering the MSE, the feature of the vibration signals can be abstracted under all kinds of conditions, and the eigenvectors were counted as the input of adaptive neuro-fuzzy inference system (ANFIS) for roller bearing fault recognition [22]. Moreover, the SE, to some extent, was undefined as no standard vectors, which were in accordance with one another. Undefined or imprecise SE leads to the degradation of authenticity of MSE algorithm. To overcome this drawback, the RCMMSE was revolved by Zhang et al. [24] to overcome these challenges. It is demonstrated that RCMMSE can not only ascend the precision of entropy assessment but also descend the probability of inducing undefined entropy.
However, after EEMD and VMD, a necessary step is how to select the number of modes. For example, the number of patterns was usually determined by the correlation coefficient and mutual information scheme between each component and the original signals [18, 25]. In addition, MSE and RCMMSE can only get the single-channel signal. Therefore, it cannot precisely display the overall signal information. The multivariate multiscale sample entropy (MMSE) was a scheme that gets different time series and conducts various embedding aspects, delay time, and amplitude ranges of data channels in a strict way. Hence, it can directly analyse multichannel data [26]. Therefore, a scheme that relied on VMD and RCMMSE was proposed in this study to get the vibration feature of roller bearings.
With the advent of computer techniques, a lot of fault recognition algorithms, such as SVM, RF, and PNN [27, 28], have been broadly applied to fault diagnosis. Therefore, the SVM, RF, and PNN were resorted to fulfill the fault recognition.
As mentioned above, a scheme based on RCMMSEVMD and VMD is designed. At first, the vibration signals were disintegrated into series IMF/BL-IMF modes by the EEMD/VMD. Secondly, the RCMMSE was utilized to figure out the entropy values. Finally, the SVM/RF/PNN models were utilized to achieve the fault recognition.
The rest of the study is organized as follows. It is shown that the review of VMD and RCMMSE schemes, respectively, is presented in Section 2. The comparison of MSE/RCMMSE and EEMD/VMD is presented in Section 3. Section 4 gives the procedures of the proposed scheme, experimental data sources, and parameter selection. Section 5 provides the experimental validation. Conclusions are given in Section 6.
2. Theoretical Framework of VMD and RCMMSE Models
2.1. Basic Principle of Variational Mode Decomposition
The VMD process is divided into establishment and solution of variational constraint problems. Assuming there is a limited bandwidth in each mode, the variational problem is formulated as mode functions , which could minimize the estimated bandwidth. The summation of each mode is constrained to be corresponded to the input signal . In particular, the structure of this problem could be divided into the following steps.
2.1.1. Variational Problem Formulation
①Hilbert transform: the analytic signal of every mode function is obtained, whose purpose is to get its spectrum:②The centre frequency of each modal signal is estimated. The spectrum of each mode is modified to be adaptive to its own baseband.③The is derived, and it denotes the signal gradient, and the signal’s bandwidth is obtained. Therefore, constrained variational questions can be obtained as follows:where is the decomposition of the mode components, , is the centre frequencies of each mode, and .
2.1.2. The Solution of the Problem
①The above constraint problem could be modified as a nonbinding problem by adding it to the penalty factor and the Lagrangian multiplication operator , and here, the second penalty factor is used to guarantee the accuracy of the signal with the Gaussian noise, and the Lagrangian operator is adopted to ensure the constraint condition satisfied strictly, and the extended Lagrangian description is as follows:②VMD used multiplication operator alternating direction method of multipliers (ADMMs) to solve the above variable problem. By alternately updating A, B, and C to seek the “saddle” of the Lagrangian expression:where is equal to , and is equal to . Equation (5) is transformed to another domain by frequency by the Parseval/Plancherel Fourier equidistant transformation:
The first item of with instead is as follows:
Equation (7) is converted into a nonnegative integral form with frequency interval:
In this case, the solution of the quadratic optimization problem is as follows:
Based on the above process, the centre frequency value is converted to the frequency domain:
The centre frequency updated scheme is as follows:where is equivalent to the current remaining amount of Wiener filtering. A is the central gravity of the spectrum of the current mode function. For inverse Fourier transmission of , then the real part is .
The procedure of VMD algorithm is as follows:(1)Initialize the , , , and .(2)Update and followed by (9) and (11).(3)Update , and(4)For a given accuracy , if , stop iteration; otherwise, return to Step 2.
2.2. Basic Principle of Refined Composite Multivariate Multiscale Sample Entropy
2.2.1. Sample Entropy
The SE uses the measurement of the exponential function for two sequences with a tolerance r from m points to remain r of each other at the next point. For a time series with N sample points, the procedure of SE calculation is shown as follows:(1)The m length vectors are formed: where m represents embedding dimension and has m consecutive values. Commencing with the ith point and generalized by removing a baseline:(2)For each , the similarity degree between the and its neighboring vector is calculated by : where denotes the maximum absolute difference in the corresponding scalar components of and .(3)For each and the fixed tolerance r, let be the number of vectors that satisfy, and then, is denoted as follows:(4)The average of the is designated as follows: where r represents the boundary width of the exponential function.(5)Increasing dimension to , steps (1)-(4) are repeated to calculate the corresponding of SE values to find , and SE is defined as follows:
When N is finite, the SE can be estimated as follows:
SE determines the time sequence irregularity on the single scale. The smaller the value of FE is, the higher time sequence self-similarity can be achieved. Conversely, the greater the FE value is, the more complicated time sequences without rules can be achieved.
2.2.2. Multiscale Entropy and Multiscale Sample Entropy
ME could be represented as a series of times with different scales, which is obtained through a coarse-grained process. The irregularity of time series is then presented by the ME, where the self-similarity of various scales is reflected. Considering the scale entropy, when there is a sequence with a higher entropy than another sequence, the corresponding complexity will also be higher than the other. That is, if there is a case that the increasing scale of a time series conversely results in the decreased entropy value, there are more chances that the time series holds a relatively simple sequence structure.
The MSE is obtained as follows:(1)Consider a discrete-time series . Firstly, the value of embedding dimension m and similar tolerance r of the SE are set. Another auxiliary time sequence is then constructed in a vector form, which is presented as and named coarse-grained vector. where is the scale factor. Actually, we can find that when , the coarse-grained time series is equal to the previous one , which means that by dividing the length of time series, we can decompose the original time series into a coarse-grained vector series . It is worth mentioning that the number of coarse-grained time series is and the length is .(2)Based on the different scales of the time sequence, the SE values are obtained. Generally, the in SE takes the standard deviation of the primary time series when calculating the value of SE. This procedure is called MFE analysis.
2.2.3. Refined Composite Multiscale Sample Entropy
In the process of applying coarse graining, the length of the original time series is declined by introducing a factor of , as shown in equation (13). For example, the authors developed the RCMMSE algorithm with the aim to boost the accuracy of the MSE. In particular, the SE of all the time series is obtained at a certain scale of factor . The RCMMSE value is defined as follows:where and . Thus, the RCMMSE value can be obtained as follows:
2.2.4. Refined Composite Multivariate Multiscale Sample Entropy
Given time series with variables, the RCMMSE is calculated according to the following procedures:(1)Each coarse-grained time series requires to compute the RCMMSE, the multiple embedding vectors , and therefore, the original is reconstructed using the multiple embedding vectors . here is the time delay vector.(2)The maximum distance of two vectors and is defined as follows:(3)For a composite time delay and threshold , the parameter is the number of the vector pairs, and therefore, , and then, the probability is computed as follows: here .(4)The mean of the is obtained as follows:
Increasing dimension to , steps (1)-(4) are repeated to get the corresponding SE values to find and .
The calculation of the RCMMSE is defined as follows:
3. Scheme Comparisons
3.1. Comparison of EEMD/VMD
To compare the EEMD and VMD models, a simulation signal is used to prove that the VMD model is better than EEMD. The simulation signal is superimposed by the multifrequency signal and the random noise signal with the standard deviation of 1, and the following formula is the signal definition.
The time-domain waveforms and its envelope spectrum of the original signal are given in Figure 1. As shown in Figure 1(b), the main frequency of the original signal focuses on 0–200 Hz, especially 30.27 Hz. Then, EEMD and VMD were used to decompose the aforementioned signal.

(a)

(b)
To achieve the VMD, firstly, the number of modes should be determined. It should be noted that if the value of in VMD is chosen too tiny, the original signal with time frequency could not be fully captured by the decomposition of the mode. Larger values produce the similar frequency between the BL-IMF components, and they may be over-decomposition. Therefore, we use the observing centre frequency of the signal scheme to select the applicable in this study [19]. The results of the centre frequency corresponding to different values are given in Table 1.
As shown in Table 1, there is a centre frequency difference of less than 0.008 kHz from BL-IMF3 to BL-IMF5 modes (0.0636 kHz in BL-IMF3, 0.0582 kHz in BL-IMF4, and 0.051 kHz in BL-IMF5) when . These three modes exist similar frequency. Therefore, the number of values is selected as 4 in this study to decompose the above simulation signal.
The EEMD/VMD models were used to decompose the original signal, and the results of IMFs (EEMD) and BL-IMFs (VMD) are given in Figures 1(a) and 2(a). Figures 1(b) and 2(b) are the corresponding spectrum analysis results. As shown in Figure 2, each BL-IMF component is mainly distributed around a single frequency (30.27 Hz), compared with VMD in Figure 1(b), and some IMF components in Figure 2(b) have a series frequency from IMF4 mode.

(a)

(b)
The mentioned result shows that the decomposition effect of EEMD algorithm is not ideal for the multicomponent synthesis simulation signal, and the mode mixing is serious. Because some slight signals drown in the signal, which is to be decomposed, EEMD in the selection process of the three-spline envelope fitting leads to decomposition bias. The weak signal is embedded in the vast majority of the strong signal where the EEMD can be filtered and extracted, but when the weak signal appears only in the maximum slope range of the strong signal, the weak signal will be in the form of wave frequency modulation and does not produce additional local extreme points. Therefore, EEMD is difficult to extract the useful components and easy to produce some components with mode mixing.
It can be seen from Figure 1(b) that VMD can not only effectively remove the dummy components, but also each BL-IMF exhibits a mode in the range of a certain scale, and there is no mode mixing problem between them. Therefore, the scale characterization and decomposition effect of VMD and better than EEMD, and this indicates that VMD has good robustness Figure 3.

(a)

(b)
3.2. Comparison of Validity of MSE/RCMMSE
In this section, the noise signal was employed to demonstrate the superior performance of the RCMMSE compared with MSE. Resorting to MSE and RCMMSE, the probabilities of inducing undefined entropy were presented to address the white noise and noises with 200 samples. For each sample, the length is chosen as 1000. The results of the corresponding probabilities are shown in Table 2. From Table 2, we can see that the probability of undefined entropy showed the same tendency as the time scale increased, while there was a decrease in the length of the time series. As given in Table 1, the probability of undefined entropy is about zero when they were used to analyse white noise in the MSE model. On the contrary, the probability of inducing undefined entropy is 0.01 when for all 200 noise samples. In the MSE algorithm, the entropy is undefined when or is zero. However, the RCMMSE can successfully get the values of entropy from 1 to 20 when the white and noises are considered. Hence, the RCMMSE model is superior to MSE.
3.3. Comparison of Accuracy in MSE/RCMMSE
In this section, the efficiency of the MSE/RCMMSE models was verified through case studies. In each simulation study, we employ noise samples with 1000 data points containing white and noises. The results of MSE/RCMMSEs are presented in Figure 4. Meanwhile, Figures 5 and 6 present the results of means and standard deviations of entropies.

(a)

(b)

(a)

(b)

(a)

(b)
As seen in Figure 4(a), the overall recently RCMMSE model shows higher performance compared with the MSE model in the condition of white noise, resulting from the tolerance in SE, which is utilized to evaluate the similarity between any two time series. Note that was often chosen as 0.15 × SD of the original time series.
It can be seen from Figures 5(a) and 5(b) that the means of the entropy values obtained using the MSE/RCMMSE are nearly equal in white noise. Nevertheless, the means of the entropy values of the MSE are higher than that of the RCMMSE (see Figure 5(b)). Figure 6 shows that the standard deviation of RCMMSE is all lower than that of MSE, and this result indicates that the entropies obtained using the RCMMSE algorithm were more consistent than those obtained using the MSE algorithm.
With the purpose of studying the relationship between the data length and the effectiveness of the MSE/RCMMSE models, several different data lengths (N = 600, 1200, 1800, and 2400) are used to get the entropies. In Table 3, the means and standard deviation of the entropies are presented. It can be found that the results of the means of entropies by MSE and RCMMSE are nearly equivalent with the long primary time series (N = 2400). However, with short original time series, the results of means show a significant difference. On the other hand, the results of the RCMMSE algorithm hold the lowest error when the discrimination between the entropy of noise with 500 data points and the entropy of noise with 2000 data points is considered. That is to say, the RCMMSE has higher performance compared with other algorithms at any length of the time series. In conclusion, the RCMMSE is more reliable than MSE according to the above discussions.
4. Experimental Setup
4.1. The Dataset Source
The performance of the proposed scheme is verified through experiments in this part. Case Western Reserve University Bearing Data are used, which was obtained from an induction motor. It should be noted that the data are generated by accelerometers at the drive end and fan end. Besides, EDM is employed to model the possible faults in the motor bearings. In particular, the fault tolerance is set as 0.1778 mm. After the fault detection, the data located in the outer race were aggregated in a state, which holds the time configuration at 6:00clock. To ensure the data acquisition system suitable for the vibration signals, the amplifier is particularly designed with a high bandwidth. On the other hand, to improve the accuracy, a sampling frequency of 12000 Hz of each channel was set for the data recorder.
Table 4 demonstrates different working conditions, which are considered in this experiment. Note that 0.1778 mm is the fault tolerance. “Normal” and “slight” denote the fault severity. In addition, considering the drive end of the motor 1785 rpm is chosen as the motor revolving speed. In total, 200 data samples are adopted with each fault condition consisting of 51 samples. In particular, in each data sample, there are 2048 data points.
4.2. Procedure of Our Proposed Scheme
The procedure of the proposed scheme can be shown as follows:(1)Resorting to EEMD and VMD models, the vibration signals under different cases were disintegrated into a sequence of IMF and BL-IMF modes.(2)Using the RCMMSE model to get the entropy values, the entropy values were regarded as the input of SVM/RF/PNN models for training and testing.(3)Testing samples were regarded as the input of the trained SVM/RF/PNN classifier, while the operating conditions can be recognized by the output of the SVM/RF/PNN classifiers. The procedure of our proposed fault diagnosis scheme is shown in Figure 7.

4.3. Parameter Selection
There are some parameters in different models that should be present before its calculation:(1)EEMD: considering EEMD, there are two parameters that need to be determined, namely the ensemble number m and related to the white noise. In particular, denotes the amplitude. It should be mentioned that there is a one-to-one correspondence between the result and the ensemble number. Note that the ensemble number should be chosen as a few hundred. In addition, consider the situation that the input signal holds the standard deviation, and a fraction of one percent of error will be caused by the remaining noise with the added noise. In particular, is set as 100.(2)VMD: the number of BL-IMF modes is according to the centre frequency. The second penalty factor in equation (4) is often set as 2000.(3)SE/RCMMSE:①SE: before the process of SE, there are three parameters that must be chosen suitably. Firstly, considering the different lengths of sequences in SE, embedding dimension should be designed. In particular, the dynamic process will be reconstructed in a more detailed manner with a larger . However, it should also be mentioned that if is chosen too large, the need of would be greatly limited, which will cause the loss of important information, and the general condition will be hard to satisfy. In the literature, is often set as 2. Another two parameters, namely, similarity tolerance and , are dependent on the gradient of the exponential function and its bound, respectively. Experimentally, , and is chosen as 0.15SD in SE.②RCMMSE: the multiple embedding vector in RCMMSE scheme is set as , and here, is the number of the BL-IMF components, and therefore, the time delay is set as 1.(4)SVM: the kernel function in SVM is selected as the radial basis function (RBF).(5)PNN: the distribution density of PNN is set as 1.5.(6)RF: there are two parameters that need to be confirmed before using RF model, such as the number of input variables is selected randomly based on the input variables. Regarding the scale factor as 20 and EEMD/VMD-RCMMSE as feature, the number of input variables , and the parameters often meet the condition [15]. Therefore, we let and the number of the DT is set as 4 and 1000 in this simulation.
5. Experimental Results and Analysis
5.1. EEMD/VMD
The data are selected from the experiments in which SKF bearings are used. The number of each sample is set as NR:1–50, IRF: 51–100, BF:101–150, ad ORF:151–200. As limited space, here with a sample of each state for an example, the time-domain waveforms of vibration signals under different working conditions are shown in Figure 8.

The vertical axis is the acceleration vibration amplitude. Because of the influence of noise, it is difficult to find significant differences in different states. As shown in Figure 8, it is hard to distinguish the four signals; in particular, there is no obvious regularity in two states of NR and BF signals.
The EEMD and VMD were used to disintegrate the vibration signals into series modes (IMF and BL-IMF). Here, we use the IRF vibration signals to observe the frequency change situation, and the waveforms of one IRF vibration signal and its spectrum envelope are shown in Figure 9. As shown in Figure 5(b), the IRF signal fault frequency is 164.1 Hz. The results of the centre frequency under different values with VMD are given in Table 5. There are some modes, such as BL-IMF4 (0.2787 kHz) and BL-IMF5 (0.3071 kHz), and with similar centre frequencies, when in Table 5, it is considered that there is over-decomposition. Therefore, the number of values is selected as 4 to decompose the vibration signals.

(a)

(b)
Employing EEMD/VMD, the roller bearing primary signals in Table 4 are decomposed into IMF1-10 and BL-IMF1-4. Figures 10 and 11 show the VMD and the corresponding spectrum of each mode for an IRF signal.

(a)

(b)

(a)

(b)
It can be shown from Figure 6(b) that the VMD can not only effectively remove the dummy components, but also each IMF in the range of a certain scale, and there is no mode mixing problem between them. Therefore, the scale characterization and decomposition effect of VMD are better than EEMD, and this result reveals that the VMD owns a good robustness.
The correlation coefficient scheme is utilized to prove the degree of correlation between each mode and the primary signal to further prove the effectiveness of VMD scheme. In Tables 6 and 7, the overall average correlation values of BL-IMF component and original IRF signals are higher than IMF. Therefore, the decomposition effect of VMD is superior to EEMD.
After the vibration signal feature extraction using VMD and EEMD, the RCMMSE was employed to get the entropy value. The value of VMD-RCMMSE values under different values and scale factors is shown in Figure 12.

Figure 12 shows that the RCMMSE values of the VMD effective mode decrease with the increasement in the scale factor and the number of modes , and the overall RCMMSE values are close to the steady state when the value of is more than 3, because the VMD algorithm decomposes the original signals into several BL-IMFs with information of different frequency bands. This can also be supported from Table 5. The correlation coefficient of most EEMD modes and the original signal is much smaller than that of VMD mode. This indicates that the VMD algorithm can overcome the mode mixing shortcoming in EEMD. The entropy values of EEMD-RCMMSE and VMD-RCMMSE (k = 4) are shown in Figure 13.

(a)

(b)
As shown in Figure 13(a), it is difficult to distinguish the four types of roller bearing signals with EEMD-RCMMSE, especially the NR and IRF, BF, and ORF. On the one hand, the complexity of signals under normal states and the conditions with ball fault is similar. In other words, the feature of the two kinds of signals is alike. On the other hand, much noise hidden in the signals may affect the identification of them. However, the RCMMSE curves of BL-IMF components can recognize the four states clearly in Figure 13(b).
5.2. Fault Identification
Regarding the extracted RCMMSE vectors as the inputs of the chosen fault classifier, such as SVM/RF/PNN, there are four kinds of data with 200 samples in total. 10, 20, and 30 samples are used under different statuses as the training samples, and the rest samples (20, 30, 40) as the testing samples are used to compare the classification accuracy. Part of the experiment classification results is shown in Figure 14 and Table 8.

(a)

(b)

(c)

(d)

(e)

(f)
As given in Table 8, the highest classification precision is up to 100% when in VMD-RCMMSE model. On the contrary, the lowest accuracy is 63.75% in EEMD-RCMMSE model. The overall classification accuracy in VMD-RCMMSE model is higher than EEMD-RCMMSE model. The overall classification accuracy in SVM model is higher than RF/PNN models. As mentioned above, the experimental results show that VMD-RCMMSE combination model can recognize different working conditions of the roller bearings effectively.
6. Conclusion
The VMD analysis model and RCMMSE are proposed in this study. The VMD model is utilized to analyse the complex vibration signals. Therefore, BL-IMF can be obtained from the decomposition of the roller bearing vibration signals, the RCMMSE can deal with multichannel data, and the EEMD/VMD-RCMMMSE values are regarded as eigenvectors. Because the VMD can achieve adaptive subdivision of each component in the frequency domain of signals, the roller bearing vibration signals are distinguished. Finally, the VMD-RCMMSE-SVM/RF/PNN combination models can distinguish roller bearing fault types effectively.
With the rapid development of Internet of things (IoT) [1, 2] and Industry 4.0 [3], there are increasingly massive real-time data from various types of mechanical equipment [4]. The availability of these data that contain abundant information about machine health has attracted more and more enterprises’ attention. It has been proved that large volume, high velocity, and diversity mechanical big data are the major properties of mechanical big data [5, 6]. Effective feature extraction from these data and accurate machinery health state evaluation with ever-accelerated updating of schemes have become hot research issues in the prognostic and health management systems in the era of industrial IoT [7].
7. Future Ideas
Most related roller bearing signal fault diagnoses depend on manual aids and expert knowledge. To diagnose more smartly, convolutional neural network (\textsc{CNN}), as well as long short-term memory (LSTM), will be considered for an end-to-end scenario. For ensuring high recognition accuracy, the fault classifier based on CNN may be used to extract the effective signal from the severely distorted signal; the LSTM could cope with time-varying bearing failure. Additionally, the aforementioned SVM fault classifier performs well when the data samples are insufficient. In the future, we may consider the generative adversarial networks (GANs) and transfer learning for lack of data. By simulating the distribution of fault samples, GAN could obtain more data. Besides, the transfer learning-related diagnosis model can reuse the previous information in the new task. Consequently, small-sized samples could achieve accurate fault identification.
Data Availability
The prior studies and data are cited at relevant places within the text as references [20].
Conflicts of Interest
The authors declare that they have no conflicts of interest.