Lifting load monitoring of mine hoist through vibration signal analysis with variational mode decomposition
Fan Jiang^{1} , Zhencai Zhu^{2} , Wei Li^{3} , Shixiong Xia^{4} , Gongbo Zhou^{5}
^{1, 2, 3, 5}School of Mechanical and Electrical Engineering, China University of Mining and Technology, Xuzhou, 221116, P. R. China
^{1, 2, 3}Key Laboratory of Mechanical and Electrical Equipment of Jiangsu Province, Xuzhou, 221116, P. R. China
^{1, 4}Postdoctoral research station of computer science and technology, China University of Mining and Technology, Xuzhou, 221116, P. R. China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 8, 2017, p. 60216035.
https://doi.org/10.21595/jve.2017.18859
Received 15 July 2017; received in revised form 29 September 2017; accepted 6 December 2017; published 31 December 2017
JVE Conferences
Mine hoists play a crucial role in verticalshaft transportation, and one of the main causes of their faults is abnormal lifting load. However, direct measurement of the load value is difficult. Further, the original structure must be destroyed for sensor installation. To facilitate efficient and accurate monitoring of the lifting load of mine hoist, this paper presents a novel conditionmonitoring method based on variational mode decomposition (VMD) and support vector machine (SVM) through vibration signal analysis. First, traditional empirical mode decomposition (EMD) is used to analyze the vibration signal collected by an acceleration sensor, and the number of obtained intrinsic mode functions (IMFs) is employed to set the VMD mode number. Second, the obtained vibration signal is processed by the parameterized VMD, and the useful IMFs of VMD are selected through correlation analysis for feature extraction. Third, the obtained features are used to train an SVM model, and the trained SVM is used to monitor the minehoist lifting load. In this study, experiments on an operated mine hoist are also conducted to verify the reliability and validity of the proposed method. The experimental results show that the proposed method can accurately identify the considered lifting load conditions.
Keywords: VMD, EMD, lifting load, mine hoist, SVM.
1. Introduction
In verticalshaft transportation, mine hoist is a very important item of equipment, which is used to transport personnel, materials, equipment, and coal between the underground and surface areas. Under operation involving frequent startstopping, high speed, and heavy loads, the accidental occurrence of a malignant fault, such as skidding, overloading, or potsticking, may cause a fatal crash and affect the safety of the mine and its production efficiency. The common faults of mine hoist are related to its lifting load. However, the mine hoist system is employed in a semiclosed deep and narrow space, and the ropes and hoisting conveyance are largespan moving parts. Therefore, it is difficult to directly measure the minehoist lifting load. Vibration signals can effectively provide health information on a large rotary machine, and many studies on vibration usage have been performed by researchers throughout the world [14]. That is, vibration signal analysis constitutes a new method of monitoring the lifting load of mine hoist.
Vibration signal is a kind of highfrequency transient signal, and it is difficult to directly use such a signal to accurately identify the lifting load conditions of mine hoist. Thus, signal processing methods should be used to extract the conditionrelated features for lifting load identification. To date, many techniques have been introduced for the analysis of vibration signals in the time, frequency, and timefrequency domains [5]. Timedomain analysis is the simplest method for feature extraction, and statistical indicators are usually obtained to monitor the mechanical health [68]. For frequencydomain analysis, the timedomain signals are mainly transformed into frequencydomain signals through Fouriertransform (FT)based signal processing methods [913], and the frequency related information is then used for health monitoring. Regardless of whether time or frequencydomain analysis is employed, however, only partial information can be obtained. To achieve more detailed analysis, timefrequency analysis methods have been introduced for signal processing, such as the shorttime FT (STFT) [14, 15], WignerVille distribution (WVD) [16], wavelet transform (WT)based methods [1719], empirical mode decomposition (EMD) [20], ensemble EMD (EEMD) [21], and variational mode decomposition (VMD) [22].
Among these methods, VMD can adaptively find the frequency center and bandwidth of each component by iteratively searching for the optimal solutions of variational models. This approach has frequently been used to analyze complicated signals in recent years. For example, in [23], an underdetermined blind source separation is proposed, in which VMD is employed to analyze the compound rollerbearing fault signals and the faultrelative features are accurately extracted. Yao et al. [24] have employed VMD and robust independent component analysis to achieve noise source identification for a diesel engine; their results show that the proposed method can accurately separate and identify the combustion noise and piston slap noise. A fault diagnosis method for multistage centrifugal pump based on VMD is designed in [25], and a comparison of the VMD and EMD fault feature extraction performance shows that the former can accurately extract the principal mode of the bearing faults more effectively. Liu et al. [26] have proposed a novel signal denoising method that combines VMD and detrended fluctuation analysis (DFA), labeled “DFAVMD”; the order of time complexity of this method is equivalent to that of EMD. In [27], a novel fault diagnosis method for wind turbine based on VMD and Teager energy operator is proposed. The above methods indicate that VMD is an effective signal processing method for complicated signals.
Through application of a signal processing method, features can be extracted. However, development of a method for employing these techniques for lifting load condition identification is also important. In this work, a novel lifting load monitoring method for a mine hoist, which is based on VMD and a support vector machine (SVM) through vibration signal analysis, is proposed. The remainder of this paper is structured as follows. Section 2 introduces the theoretical basis, which includes overviews of EMD and VMD. The recognition model of the minehoist lifting load is described in Section 3. Section 4 presents a verification of the proposed signal processing method, which is achieved through simulation analysis. Finally, the experimental study is presented and discussed in Section 5, and the conclusion is given in Section 6.
2. Theoretical basis
2.1. EMD
EMD, developed by Huang et al. [20], can selfadaptively decompose a complicated signal into a series of intrinsic mode functions (IMFs), based on the local characteristic time scale of signal. In EMD application, there are three assumptions: (1) the target signal has at least two extrema (one maximum and one minimum); (2) the characteristic time scale is defined by the time lapse between the extrema; and (3) if the data are totally devoid of extrema but contain inflection points only, they can then be differentiated one or more times to reveal the extrema. The EMD algorithm process is described as follows [20, 2830].
1) The upper and lower envelopes are generated using a cubic spline line to fit the local extrema points of the signal $x\left(t\right)$, and the mean of the upper and lower envelopes is defined as ${m}_{1}$. Then, we set:
For a strict IMF, the two requirements given in [20] should be met. If ${h}_{1}$ meets the two standards of the strict IMF, the result of Eq. (1) is then defined as the first IMF; otherwise, ${h}_{1}$ is regarded as the original signal $x\left(t\right)$, and the above steps are repeated until the result ${h}_{1k}$ is a strict IMF. Subsequently, the first IMF component is obtained and labeled:
2) The first IMF ${c}_{1}$ is separated according to:
where ${r}_{1}$ is the residue. The above steps are repeated $N$ times until the algorithm satisfies the termination condition. Then, we can obtain:
3) Finally, the signal $x\left(t\right)$ can be decomposed to:
2.2. VMD
VMD, a new selfadaptive and quasiorthogonal signal processing method, was proposed by Dragomiretskiy and Zoss [22] based on Wiener filtering, onedimensional Hilbert transform, and heterodyne demodulation. VMD can decompose a complicated multicomponent signal into a serious of subsignals, named IMF components, each of which is mostly compact around a center pulsation and limitedfrequency bandwidth. To evaluate the bandwidth of a mode, the following constrained variational problem should be solved:
where ${u}_{k}$ and ${\omega}_{k}$ are shorthand notations for the set of all modes and their center frequencies, respectively; $t$ is the time script, and $\delta $ is the Dirac distribution.
Theoretically, various methods to solve the reconstruction constraint can be introduced. In [22], a quadratic penalty term and Lagrangian multipliers are employed to render the problem unconstrained, and the following new expression of the solution can be obtained:
where $\alpha $ is the datafidelity constraint parameter and $\lambda $ is the Lagrangian multiplier.
Here, a sequence of iterative suboptimizations called the alternate direction method of multipliers (ADMM) is introduced to address the above formula. The solution to the original minimization problem of Eq. (6) is now found as the saddle point of Eq. (7). Then, the modes ${u}_{k}$ and their corresponding center frequency ${\omega}_{k}$ can be updated as:
Then, the following two formulas can be obtained:
Through continuous iteration updates, all subsignals, called IMF models, can be decomposed from the solution and are described as follows:
The complete VMD algorithm is summarized as follows:
Step 1: Initialize $\left\{{\widehat{u}}_{k}^{1}\left(\omega \right)\right\}$, $\left\{{\widehat{\omega}}_{k}^{1}\left(\omega \right)\right\}$, and ${\lambda}^{1}$, and set $n=$ 0;
Step 2: Update $n$ with $n+$ 1 and perform a loop;
Step 3: Update ${\widehat{u}}_{k}$ for all $\omega >$ 0 with the formula shown in Eq. (12);
Step 4: Renew ${\omega}_{k}$ according to ${\omega}_{k}^{n+1}$ based on the relation:
Step 5: Implement dual ascent for all $\omega >$ 0, where:
Step 6: Repeat steps 2 to 5 until the termination criterion of the following form is satisfied:
3. Proposed method
Vibration signal analysis is an effective means of accurately diagnosing mechanical faults and monitoring mechanical conditions, because these signals usually carry abundant staterelative information. Mine hoist is a typical item of rotating machinery used for verticalshaft transportation, and its common fault is abnormal lifting load. Here, vibration signals are analyzed to monitor the lifting load of mine hoist. However, these vibration signals are both nonlinear and unsteady. Further, it is difficult to extract the loadrelative information from these signals; thus, monitoring the lifting load of mine hoist using tradition signal processing methods is challenging. For EMD, its main purpose is to adaptively decompose a mixed signal into a series of IMFs. Theoretically, if the mixed signal is main consist of $k$ components, it can get $k$ IMFs by EMD. That is to say, EMD can get the number of components of a mixed signal. For VMD, it is also an algorithm to decompose a complex signal into its compositions with IMFs. The mode number is an important parameter for VMD and which is mean the number of the components in a complex signal. Therefore, a new signal processing approach, combined the advantages of EMD and VMD, is tentatively proposed in this paper.
SVM has been developed from statistical learning theory to solve pattern recognition and regression analysis problems using structural risk minimization principles [31, 32]; these principles can realize the optimal combination between learning precision and identification accuracy with only a small number of samples. Further, SVMbased classification can be regarded as a quadratic optimization problem. This technique first maps input samples into a particular highdimensional space with nonlinear transformations and then finds an optimal plane for the classes in highdimensional space to classify the input samples. In addition, SVM has stronger generalization ability than the nonlinear function approximation method, and can also overcome difficulties in determining the neural network structure and the local minimum.
Based on the above analysis, this study presents a new condition monitoring method, based on EMD, VMD and SVM with vibration signal analysis, to precisely monitor the lifting load of mine hoist. The flow chart of the proposed method is shown in Fig. 1. First, traditional EMD is used to process the vibration signals collected by an acceleration sensor, and the number of obtained IMFs is used to set the VMD mode number. Second, the vibration signal is processed by the parameterized VMD and the useful IMFs of the VMD are selected through correlation analysis for feature extraction. Third, the obtained features are evaluated to train an SVM model, and the trained SVM is used to monitor the minehoist lifting load.
Fig. 1. Flow chart of proposed method
4. Simulation analysis
To test the superiority of the proposed method, a simulation experiment was conducted, which is reported in this section. The test signals used in the simulation are as follows:
where ${T}_{1}=$ 1000, ${T}_{2}=\text{0.25}*i$, $i=$ 1, 2, 3, and $e\left(t\right)$ is noise that is equivalent to 0.2 times the standard deviation of the synthetic signal of ${y}_{1}$, ${y}_{2}$, and ${y}_{3}$.
Fig. 2 shows the timedomain waves of the simulation signals, i.e., the highfrequency oscillation signal (${y}_{3}$), lowfrequency waves (${y}_{1}$ and ${y}_{2}$), and composite signal (${y}_{com}$). From this figure, the difficulty in identifying the components of the simulated ${y}_{com}$ is apparent. The decomposition results provided by the traditional EMD and EEMD are plotted in Fig. 3. Directly observing the plots in Fig. 3(a), it is apparent that the main components are decomposed into ${c}_{3}$ to ${c}_{6}$. From the plots in Fig. 3(b), it is apparent that most of the information on ${y}_{1}$, ${y}_{2}$, and ${y}_{3}$ is decomposed into ${c}_{6}$, ${c}_{4}$, and ${c}_{3}$. However, part of ${y}_{1}$ is falsely decomposed into IMF ${c}_{5}$, because the parameter settings are not optimal. By comparing Figs. 3(a) and (b), it is apparent that EEMD can yield a superior decomposition result for a complex signal than traditional EMD.
Fig. 2. Timedomain waves of simulation signals
a)
b)
c)
d)
The above analysis reveals that traditional EMD and EEMD have certain limitations in the context of complicated signal processing, such as the mode mixing problem and the selection of optimal parameters. It has been verified that VMD, a new technique, can effectively manage complicated signals; here, it was used to decompose the composite signal shown in Fig. 2. According to the principles of the proposed method, the number of IMFs in the decomposition process of traditional EMD is 10; thus, the mode number of VMD was set to 10. Based on the results reported in [25], the balancing parameter of the datafidelity constraint was set to 1500. Fig. 4 shows the decomposition results of the EMDbased VMD method, where ${y}_{1}$, ${y}_{2}$, and ${y}_{3}$ are successfully decomposed into the front three IMFs ${c}_{1}$, ${c}_{2}$, and ${c}_{3}$.
To further analyze the decomposition effect of the EMDbased VMD method, the decomposition errors of the above three simulation signals are shown in Figs. 5(ac). Hence, we conclude that the EMDbased VMD method can effectively decompose the complicated simulated ${y}_{com}$ and yield superior results to the traditional EMD and EEMD.
Fig. 3. Decomposition results obtained using: a) traditional EMD, b) traditional EEMD with default parameters defined in [21]
a)
b)
Fig. 4. Decomposition results of the EMD based VMD
Fig. 5. Decomposition errors of simulation signals: a) ${y}_{1}$, b) ${y}_{2}$, c) ${y}_{3}$
a)
b)
c)
5. Experimental study
5.1. Experimental settings and data collection
To verify the efficacy of the proposed method in the context of lifting load monitoring with vibration signals, an experiment on a running mine hoist was conducted. Fig. 6 shows a schematic diagram of a multirope friction hoist and the data acquisition system used in this experiment. As shown in this figure, the mine hoist contained several ropes, two containers, and a main shaft apparatus. Usually, a multirope friction hoist is operated based on friction between ropes and pads installed in the main shaft apparatus.
Fig. 6. Schematic diagram of mine hoist and data acquisition system
Thus, the main shaft apparatus is an important component of mine hoist, because it is not only a key transmission system of the mine hoist, but also a main loadcarrying part. Therefore, the vibration signals collected from the main shaft apparatus can reveal characteristic information on the lifting load. In this work, the vibration signals were collected by the data acquisition equipment and a sensor placed on the top of the bearing housing of the main shaft apparatus with a strong magnetic foundation, as shown in Fig. 6. The data acquisition equipment included a data acquisition board and data processor, and the acquisition software system was selfdeveloped using Visual Basic Language.
Here, five different lifting loads were simulated for the mine hoist, as listed in Table 1. The sampling frequency used to collect the vibration signals was set to 1000 Hz. Fig. 7 shows the vibration signals and their FFT spectrums for these five considered lifting load conditions. From this figure, it is apparent that some nuances exist between the amplitudes of the vibration signals obtained under different lifting load conditions. However, it is difficult to use these vibration signals directly to identify the lifting load conditions without any other obvious characteristics. Thus, features must be extracted from the collected vibration signals to reveal the information on the minehoist lifting load. Then, the correlations between the features and lifting load conditions can be constructed to determine if the mine hoist is operating under a normal lifting load. Recall that the main purpose of the proposed method is to facilitate monitoring and determination of the minehoist lifting load through vibration signal analysis alone. Further it is meaningful to determine whether the mine hoist is operating under normal or critical lifting load conditions.
Table 1. Conditions of five lifting experiments
Experiment no.

Label

Lifting load (t)

Conditions

1

S1

0

Noload

2

S1

0

Noload

3

S2

8.5

Normal

4

S2

10.5

Normal

5

S3

11.6

Over load

Fig. 7. Vibration signals and their FFT spectrums for considered conditions
a)
b)
5.2. Feature extraction
Feature extraction based on signal processing tools has a significant influence on condition monitoring or classification efficacy. Figs. 8, 9 show real decomposed results for the vibration signals, which were obtained using EMD, EEMD, and the proposed method. Hence, it is apparent that the conditionrelative information is not decomposed from the complicated vibration signals when traditional EMD and EEMD are used. However, the first several IMFs for the EMDbased VMD method exhibit superior consistency and periodicity. That is, our method can effectively analyze the collected vibration signals.
Fig. 8. Decomposition results of real vibration signals obtained using: a) traditional EMD, b) traditional EEMD with default parameters defined in [21]
a)
b)
Through use of the EMDbased VMD method, a set of IMFs can be obtained and their sensitivity for condition monitoring can be assessed. Theoretically, not all IMFs contain sufficient useful information for lifting load identification. To reduce the feature dimensions and improve the processing speed of the SVM model, a correlation analysis was conducted to quantify the usefulness of each IMF and the five most sensitive IMFs were selected to calculate permutation entropy (PE), defined in [33], as foundation features. That is, each selected vibration signal segment can product five considered IMFs and each is able to obtain one PE values, so five PE values can be calculated from each selected vibration signal segment.
However, noise commonly exits in vibration signals. To reduce the effect of noise interference on the monitoring of lifting load, we used statistical PE features, each is the mean of 15 foundation features of the same IMF, to construct samples. Therefore, in this paper, the used sample are consistent of five statistical PE features. In this work, 60 samples were extracted for each considered lifting load, 20 of which were used to train the SVM model and the remainder of which were set for testing. The used samples and target outputs of the SVM are listed in Table 2. Note that the other parameters were set to the default values given in [34].
Fig. 9. Decomposition results of real vibration signals obtained using EMDbased VMD method
Table 2. Monitoring results obtained using proposed method
Conditions

Training samples

Test samples

Target outputs

S1

20

40

1

S1

20

40

1

S2

20

40

2

S2

20

40

2

S3

20

40

3

5.3. Results and discussion
In this study, the vibration signals collected from a mine hoist under the stage of uniform motion were used to establish a correlation model between the features and lifting load conditions. Fig. 10 shows the monitoring results of the proposed method, in which the abscissa indicates the number of test samples and the ordinate expresses the outputs of the trained SVM. According to the designed SVM model, the theoretical outputs of the test samples numbered 180, 81160, and 161200 should be “1”, “2,” and “3”, respectively. From the plots in Fig. 10, only one of these outputs of the test samples numbered 180 is not “1”; that is to say, these test samples, except one, are identified as having a noload condition, which is consistent with the actual situation.
Fig. 10. Test results for proposed method
However, for the test samples numbered 81160, some outputs are not “2”, but “1”and “3”. That is, some test samples originally corresponding to a normal load condition are mistakenly identified as corresponding to noload conditions and overload conditions. Similar situations also appear for the test samples numbered 161200, and some test samples are mistakenly identified as corresponding to the normal load condition. Further analysis revealed that the phenomenon in which more than two continuous adjacent test samples are mistakenly identified did not occur, and an incorrect test result for two continuous adjacent test samples occurred only once. Therefore, we can use the newest result of two or more continuous adjacent test samples as the final classification result.
To confirm the superiority of the proposed method, comparison experiments were conducted, which were labeled “SVM+EMD”, “SVM+EEMD”, and “SVM+VMD”. In these comparisons, the classifiers were all SVM models and the number of samples used to train and test the SVM model were also the same as those for the proposed method. The permutation entropy values were also all extracted from five sensitive IMFs selected via correlation analysis, and the statistical features were all calculated using 15 permutation entropy values. For the SVM+EEMD method, the EEMD was employed using the default parameters defined in [21]. For the SVM+VMD method, the mode number and balancing parameter of the datafidelity constraint were discretionarily set to 5 and 1500, respectively. Table 3 presents the test results for these comparison methods. From this table, it is apparent that the SVM+VMD method has the worst monitoring result. That is to say, the feature extraction of SVM+VMD is not suitable as the classifier for these comparison methods are also SVM. Further analysis the monitoring results of SVM+VMD reveals that the VMD with random configuration parameters cannot decompose lifting load related information into IMFs successfully.
Table 3. Monitoring results obtained using proposed method
Experiments

Training samples

Test samples

Total accuracy (%)

SVM+EMD

100

200

87

SVM+EEMD

100

200

89

SVM+VMD

100

200

81

Proposed method

100

200

92

Therefore, setting suitable mode number and balancing parameter of the datafidelity constraint has great influence on feature extraction and condition monitoring. Further, the analysis of the monitoring results of SVM+EMD and SVM+EEMD shown that their monitoring accuracies are almost the same but are all lower than that of the proposed method. From this comparative analysis, we can conclude that the proposed method can successfully recognize the lifting load conditions of the mine hoist.
6. Conclusions
Lifting load monitoring is an effective approach to achieve a more secure operation of a mine hoist. Focusing on this issue, in this paper, a novel classification method based on variational mode decomposition (VMD) and support vector machine (SVM) through vibration signal analysis was proposed to facilitate accurate monitoring of the lifting load of mine hoist. In this work, an empirical mode decomposition (EMD)based VMD was developed and the features of the selected sensitive intrinsic mode functions (IMFs) were extracted to construct an SVMbased classifier for monitoring the minehoist lifting load. Experiments on minehoist operation were conducted to verify the reliability and validity of the proposed method. The experimental results shown that the proposed method can accurately identify the considered lifting load conditions. Some advantages of the proposed method are the following:
1) It is difficult to directly measure the lifting load of a running mine hoist with a force transducer, because the original structure must be destroyed to install the sensor. However, vibration signals can be collected easily by installing sensors on the main shaft apparatus of the mine hoist. In this study, the minehoist lifting load was successfully identified by analyzing the vibration signals; thus, a new approach for lifting load monitoring was realized.
2) Complex signal processing is challenging in condition monitoring. It has been shown that VMD exhibits superior performance with respect to the processing of nonlinear and nonstationary signals compared to EMD and EEMD. However, it is difficult to preset the VMD number mode. This work presented an EMDbased VMD method to solve this problem, and the experimental results confirmed its improvement over traditional VMD.
Acknowledgements
This work was supported by the Fundamental Research Funds for the Central Universities (2017QNA17).
References
 Kim S. M., Brennan M. J., Abreu G. L. C. M. Narrowband feedback for narrowband control of resonant and nonresonant vibration. Mechanical Systems and Signal Processing, Vol. 76, Issue 77, 2016, p. 4757. [CrossRef]
 Xu X. G., Liu H. X., Zhu H., Wang S. L. Fan fault diagnosis based on symmetrized dot pattern analysis and image matching. Journal of Sound and Vibration, Vol. 374, 2016, p. 297311. [CrossRef]
 Sawalhi N., Randall R. B. Vibration response of spalled rolling element bearings: Observations, simulations and signal processing techniques to track the spall size. Mechanical Systems and Signal Processing, Vol. 25, Issue 3, 2011, p. 846870. [CrossRef]
 Barszcz T., Randall R. B. Application of spectral kurtosis for detection of a tooth crack in the planetary gear of a wind turbine. Mechanical Systems and Signal Processing, Vol. 23, Issue 4, 2009, p. 13521365. [CrossRef]
 Jardine A., Lin D., Banjevic D. A review on machinery diagnostics and prognostics implementing conditionbased maintenance. Mechanical Systems and Signal Processing, Vol. 20, Issue 7, 2007, p. 14831510. [CrossRef]
 Jiang F., Zhu Z. Z., Li W., Chen G. A., Zhou G. B. Identification and diagnosis of concurrent faults in rotorbearing system with WPT and zero space classifiers. Journal of Vibroengineering, Vol. 16, Issue 2, 2014, p. 901912. [CrossRef]
 Lei Y. G., He Z. J., Zi Hu Y. Y. Q. Fault diagnosis of rotating machinery based on multiple ANFIS combination with Gas. Mechanical Systems and Signal Processing, Vol. 21, Issue 5, 2007, p. 22802294. [CrossRef]
 Abboud D., Antoni J., Eltabach, Sieg Zieba S. Angle\time cyclostationarity for the analysis of rolling element bearing vibrations. Measurement, Vol. 75, 2015, p. 2939. [CrossRef]
 Miao Q., Cong L., Pecht M. Identification of multiple characteristic components with high accuracy and resolution using the zoom interpolated discrete Fourier transform. Measurement Science and Technology, Vol. 22, Issue 5, 2011, p. 112. [CrossRef]
 Borghesani P., Pennacchi P., Randall R. B., Sawalhi N., Ricci R. Application of cepstrum prewhitening for the diagnosis of bearing faults under variable speed conditions. Mechanical Systems and Signal Processing, Vol. 36, Issue 2, 2013, p. 370384. [CrossRef]
 Figlus T., Stanczyk M. A method for detecting damage to rolling bearings in toothed gears of processing lines. Metalurgija, Vol. 55, Issue 1, 2016, p. 7578. [CrossRef]
 Antoni J. Cyclic spectral analysis of rollingelement bearing signals: facts and fictions. Mechanical Systems and Signal Processing, Vol. 304, Issues 35, 2007, p. 497529. [CrossRef]
 Antoni J., Randall R. B. Differential diagnosis of gear and bearing faults. Journal of Vibration and AcousticsTransactions of the ASME, Vol. 124, Issue 2, 2002, p. 165171. [CrossRef]
 Antoni J., Randall R. B. The spectral kurtosis: application to the vibratory surveillance and diagnostics of rotating machines. Mechanical Systems and Signal Processing, Vol. 20, Issue 2, 2006, p. 308331. [CrossRef]
 Sawalhi N., Randall R. B. Gear parameter identification in a wind turbine gearbox using vibration signals. Mechanical Systems and Signal Processing, Vol. 42, Issues 12, 2015, p. 368376. [CrossRef]
 Zhou Y., Chen J., Dong G. M., Xiao W. B., Wang Z. Y. WignerVille distribution based on cyclic spectral density and the application in rolling element bearings diagnosis. Proceedings of the Institution of Mechanical Engineers Part C: Journal of Mechanical Engineering Science, Vol. 225, Issue 12, 2011, p. 28312847. [Publisher]
 Kankar P. K., Sharma S. C., Harsha S. P. Rolling element bearing fault diagnosis using autocor relation and continuous wavelet transform. Journal of Vibration and Control, Vol. 17, Issue 3, 2011, p. 20812094. [CrossRef]
 Zhong G. S., Ao L. P., Zhao K. Influence of explosion parameters on wavelet packet frequency band energy distribution of blast vibration. Journal Central South University of Technology, Vol. 19, Issue 9, 2012, p. 26742680. [CrossRef]
 Figlus T., Stanczyk M. Diagnosis of the wear of gears in the gearbox using the wavelet packet transform. Metalurgija, Vol. 53, Issue 4, 2014, p. 673676. [CrossRef]
 Huang N., Shen Z., Long S., Wu M., Shih H., Zheng Q., Yen N., Tung C., Liu H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Proceedings of the Royal Society of London, Vol. 454, Issue 1971, 1998, p. 903995. [Publisher]
 Wu Z., Huang N. E. Ensemble empirical mode decomposition: a noiseassisted data analysis method. Advances in Adaptive Data Analysis, Vol. 1, Issue 1, 2009, p. 141. [CrossRef]
 Dragomiretskiy K., Zoss O. D. Variational mode decomposition. IEEE Transactions on Signal Processing, Vol. 62, Issue 3, 2014, p. 531544. [CrossRef]
 Tang G., Luo G., Zhang W., Yang C., Wang H. Underdetermined blind source separation with variational mode decomposition for compound roller bearing fault signals. Sensors, Vol. 16, Issue 6, 2016, p. 117. [CrossRef]
 Yao J., Xiang Y., Qian S., Wang S., Wu S. Noise source identification of diesel engine based on variational mode decomposition and robust independent component analysis. Applied Acoustics, Vol. 116, 2017, p. 184194. [CrossRef]
 Zhang M., Jiang Z., Feng K. Research on variational mode decomposition in rolling bearings fault diagnosis of the multistage centrifugal pump. Mechanical Systems and Signal Processing, Vol. 93, Issue 15, 2017, p. 460493. [CrossRef]
 Liu Y., Yang G., Li M., Yin H. Variational mode decomposition denoising combined the detrended fluctuation analysis. Signal Processing, Vol. 125, 2016, p. 349364. [CrossRef]
 Zhao H., Li L. Fault diagnosis of wind turbine bearing based on variational mode decomposition and Teager energy operator. IET Renewable Power Generation, Vol. 11, 2016, p. 453460. [CrossRef]
 Ricci R., Pennacchi P. Diagnostics of gear faults based on EMD and automatic selection of intrinsic mode functions. Mechanical Systems and Signal Processing, Vol. 25, Issue 3, 2011, p. 821838. [CrossRef]
 Xiong Q., Xu Y. H., Peng Y. Q., Zhang W. H., Li Y. J., Tang L. Lowspeed rolling bearing fault diagnosis based on EMD denoising and parameter estimate with alpha stable distribution. Journal of Mechanical Science and Technology, Vol. 31, Issue 4, 2017, p. 15871601. [CrossRef]
 Ahn J. H., Kwak D. H., Koh B. H. Fault detection of a rollerbearing system through the EMD of a wavelet denoised signal. Sensors, Vol. 14, Issue 8, 2014, p. 15022–15038. [CrossRef]
 Vapnik V. N. The Nature of Statistical Learning Theory. Springer, New York, 1999. [CrossRef]
 Bordoloi Tiwari D. J. R. Optimum multifault classification of gears with integration of evolutionary and SVM algorithms. Mechanism and Machine Theory, Vol. 73, 2014, p. 4963. [CrossRef]
 Bandt C., Pompe B. Permutation entropy: a natural complexity measure for time series. Physical Review Letters, Vol. 88, Issue 17, 2002, p. 174102. [CrossRef]
 Chang C. C., Lin C. J. LIBSVM – a library for support vector machines. 2016, http://www.csie.ntu.edu.tw/~cjlin/libsvm/. [CrossRef]