Next Article in Journal
Effective Construction of a High-Capacity Boronic Acid Layer on a Quartz Crystal Microbalance Chip for High-Density Antibody Immobilization
Next Article in Special Issue
A Switched-Element System Based Direction of Arrival (DOA) Estimation Method for Un-Cooperative Wideband Orthogonal Frequency Division Multi Linear Frequency Modulation (OFDM-LFM) Radar Signals
Previous Article in Journal
GPS Multipath Analysis Using Fresnel Zones
Previous Article in Special Issue
Integration of Terrestrial Laser Scanning and NURBS Modeling for the Deformation Monitoring of an Earth-Rock Dam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Multiclass Mahalanobis Taguchi System for Bearing Fault Diagnosis under Variable Conditions

1
State Key Lab of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing 100044, China
2
National Engineering Laboratory for System Safety and Operation Assurance of Urban Rail Transit, Guangzhou 510000, China
3
Beijing Research Center of Urban Traffic Information Sensing and Service Technologies, Beijing Jiaotong University, Beijing 100044, China
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(1), 26; https://doi.org/10.3390/s19010026
Submission received: 16 November 2018 / Revised: 18 December 2018 / Accepted: 19 December 2018 / Published: 21 December 2018
(This article belongs to the Special Issue Sensor Signal and Information Processing II)

Abstract

:
Bearings are vital components in industrial machines. Diagnosing the fault of rolling element bearings and ensuring normal operation is essential. However, the faults of rolling element bearings under variable conditions and the adaptive feature selection has rarely been discussed until now. Thus, it is essential to develop a practicable method to put forward the disposal of the fault under variable conditions. Considering these issues, this paper uses the method based on the Mahalanobis Taguchi System (MTS), and overcomes two shortcomings of MTS: (1) MTS is an effective tool to classify faults and has strong robustness to operating conditions, but it can only handle binary classification problems, and this paper constructs the multiclass measurement scale to deal with multi-classification problems. (2) MTS can determine important features, but uses the hard threshold to select the features, and this paper selects the proper feature sequence instead of the threshold to overcome the lesser adaptivity of the threshold configuration for signal-to-noise gain. Hence, this method proposes a novel method named adaptive Multiclass Mahalanobis Taguchi system (aMMTS), in conjunction with variational mode decomposition (VMD) and singular value decomposition (SVD), and is employed to diagnose the faults under the variable conditions. Finally, this method is verified by using the signal data collected from Case Western Reserve University Bearing Data Center. The result shows that it is accurate for bearings fault diagnosis under variable conditions.

1. Introduction

Rolling element bearings have wide applications in industrial machines and are one of the most critical components. If faults occur in bearings, equipment could be damaged and disasters might happen consequently. Therefore, it is essential to monitor the health conditions of bearings. The analysis of vibration signals has been a hot research pot and used to detect faults of bearings. It is crucial to recognize faults occurring in bearings and avoid fatal breakdowns as early as possible. For decades, many researchers have conducted extensive research on fault diagnosis. At present, the fault diagnosis methods are divided into model-based methods and data-driven methods. The model-based methods generally build on the physics of the process, generating the residuals between the measure process variables and estimates [1], such as Hidden Markov Modeling (HMM) which is successfully applied to bearing fault detection and diagnosis [2]. Autoregressive modelling [3] also has had excellent performance in bearing fault diagnosis. The accelerated degradation testing (ADT) [4] method is useful in fault diagnosis and lifespan prediction, and reference [5] presents a new approach using observer-based residual generation with no complicated design constraints to establish the relationship between the state estimation error and the fault signal. For the data-driven method, which is based on the historical data and does not need accurate mathematic and priori knowledge, it has wide applications in fault diagnosis. For example, Bayesian network is an excellent data-driven diagnosis method [6,7,8]. Machine learning algorithms, such as K-nearest Neighbor (KNN), Deep Convolutional Neural Networks (DCNN), and auto-encoders are effective in the fault diagnosis of various bearings [9,10,11,12]. There is also decision tree [13], Support Vector Machine (SVM) [14], and wavelet transform [15]. Artificial neural network (ANN) and Trace Ratio Linear Discriminant Analysis are other methods of diagnosing the bearings fault [16,17]. In addition, the vibration signal under the various operating conditions (especially in low rotational speed) is non-stationary and non-linear, the characteristic defect frequencies move continuously with the change of rotating speed, and is the same as the bearing that is going to break down. If the bearing progresses toward failure, the nonlinear features also start to be dominated by stochastic signal. With the vibration signal measured for diagnosing the bearing fault under variable condition, it will be difficult to diagnose the fault of bearing by using the traditional methods. Reference [18] proposed a signal selection scheme based upon two order tracking techniques from complicated non-stationary operational measured vibrations. Reference [19] reviewed features extraction methods and its application on bearing vibration signal and presents an empirical study of feature extraction methods in low rotational speed. Reference [20] proposes the estimation of instantaneous speed relative fluctuation (ISRF) in a vibration speed. Reference [21] proposes Stacked Convolutional Autoencoders (SCAE) together with DCNN in stationary and non-stationary speed operation. Graph-based rebalance semi-supervised learning (GRSSL) [22], weighted self-adaptive evolutionary extreme learning machine (WSaE-ELM) [23] and Singular Spectrum Analysis [24] are effective in diagnosing the fault under variable conditions.
However, although the aforementioned methods are effective for bearing faults diagnosis, the feature selection part is often non-adaptive or unexplainable. The conventional methods mainly extract and select features manually, which relies heavily on the experts’ knowledge and experience. Since the signals acquired in the real world might be various in many different aspects, the features selected manually might be sensitive in the variations of operation conditions and import inevitable errors for fault diagnosis. The deep learning algorithms can extract features automatically and overcome this drawback. However, deep learnings are data-hungry and require plenty of training data which are hardly acquired in practice, especially the faulty data under different conditions. Besides, the features acquired by deep learnings are unexplainable. Therefore, the aforementioned methods can hardly diagnose the faults under various operation conditions in practice. In contrast, the Mahalanobis Taguchi System (MTS) is more robust than the other methods in various operating conditions [25].
MTS offers a tool to determine important features and optimize the system. It is a different form compared with the other classification methods, because this classification model of measurement scale is constructed by using the class samples. It is useful to diagnose bearing faults under various conditions because the different pattern could be identified by using the Mahalanobis distance (MD) and Taguchi method. In this paper, MD is used to calculate the distance of the correlations between the benchmark and others, and the distance could be measured without the volatility of data. The advantage of MD is that it takes into consideration the correlations between the features and this consideration is very important in pattern analysis, which is why MTS is suitable for bearing fault diagnosis under various conditions [26]. On the other hand, the Taguchi method is used to select features without manual intervention, which could improve the robustness of the algorithm. MTS also offers an effective tool for multivariate analysis [27], considering that the bearing faults can be classified according to locations, such as inner race, outer race and rolling element [28]. However, when the conventional MTS is used for bearing fault diagnosis, misclassifications might occur due to less adaptivity of the threshold configuration for signal-to-noise gain. During the feature selection, the threshold is normally set as a constant, which might result in the overfitting problem. If the threshold value is too large, some critical features might be eliminated. On the contrary, if the threshold value is too low, some useless or harmful feature might be selected. Therefore, if the threshold value does not march the training data sufficiently, misclassifications emerge.
To overcome this drawback, this paper presents a novel method named adaptive Multiclass Mahalanobis Taguchi system (aMMTS) for bearing fault diagnosis. This method employs the MTS for multi-classification by considering different conditions as different benchmarks respectively. The results are based on the minimum MDs between the data and each benchmark data, and the label of the data is determined to be consistent with the label of the benchmark data whose MD is minimum. The method can be described briefly as follows: Firstly, after the Mahalanobis space (MS) is constructed by the two-level orthogonal array of the Taguchi method, aMMTS calculates the MDs from the data to the benchmark data and obtains features’ signal-to-noise ratios (SNRs) and gain values by using two-level orthogonal array. Secondly, the features are selected adaptively by recalculating the MDs via rearranging the order of features’ gain values by ascending and descending. Therefore, the proposed method is able to select the best classification result according to the adaptive chosen sequence of features’ SNRs instead of a hard threshold. Here, the sequence of SNR is determined by a function, which selects several maximum or minimum features to calculated MDs. Finally, a set of features with the best results is selected as the final feature vector. In this method, two different sets of training samples are employed to calculate the SNRs respectively and obtain the final feature vectors respectively. By the aforementioned improvement, the proposed aMMTS is capable to overcome the drawback of the conventional MTS and prevent the over-fitting problem. Therefore, the aMMTS is insensitive to the operation conditions and can be employed for bearing fault diagnosis.
Moreover, this method is combined with variational mode decomposition (VMD) [29] and singular value decomposition (SVD) to diagnose the faults. VMD is an entirely non-recursive algorithm, and is used to decompose the signal. It has been proven that due to the characteristics of nonlinear vibration in the bearings, VMD is more efficient than empirical mode decomposition (EMD) and Fourier transform (FT) under variable condition. SVD is used to extract the features. Therefore, VMD and SVD are employed in this paper.
The rest of this article is organized as follows: Section 2 introduces the algorithms involved in this paper. Section 3 illustrates the experiments to validate the proposed method. Section 4 is the conclusions.

2. Methodology

In this paper, the main steps of fault diagnosis are signal decomposition, feature extraction and fault detection. The detailed process and method are as follows:
Step1: The VMD in conjunction with wavelet denoising is employed to eliminate the noises and decompose the raw signals;
Step2: Extracting features from the decomposed signals by SVD;
Step3: The proposed aMMTS is employed for the fault diagnosis.
The steps of this method are shown in Figure 1.

2.1. VMD

After the wavelet denoising is used to estimate the noise of the raw signal, this paper employs VMD to decompose non-stationary signals. VMD can decompose a signal into different simple intrinsic mode functions, whose frequency center and bandwidth are band-limited and determined by iterative searching for the optimal solution of the variational model. The constrained formula is given as:
min μ k , ω k { k | | t [ ( δ ( t ) + j π t ) * μ k ( ) t ] | e j ω k t | | 2 2 } s . t . μ k = f
where μ k is the sub-signals, ω k represents the center frequency of sub-modes. The optimal solution can be solved as the minimization problem, which could be addressed by introducing a quadratic penalty and Lagrangian multipliers [30]:
L ( { μ k } { ω k } , λ ) = α k | | t [ ( δ ( t ) + j π t ) * μ k ( t ) ] | e j ω k t | | 2 2 + | | f ( t ) k μ k ( t ) | | 2 2 + λ ( t ) , f ( t ) k μ k ( t )
α denotes the balancing parameter of the constraint.
All sub-signals μ k are updated for all ω 0 as follow:
μ ^ k n + 1 f ^ i < k μ ^ i n + 1 i > k μ ^ i n + λ ^ n 2 1 + 2 α ( ω ω k n ) 2
μ ^ k 1 , ω ^ k 1 and λ ^ 1 are initialized to all zeroes.
All center frequency of sub-modes ω k are updated as follow:
ω k 0 ω | μ ^ k n + 1 ( ω ) | 2 d ω 0 | μ ^ k n + 1 ( ω ) | 2 d ω
End for
λ ^ n + 1 λ ^ n + τ ( f ^ k μ ^ k n + 1 )
Until:
k μ ^ k n + 1 μ ^ k n 2 2 μ ^ k n 2 2 < ε

2.2. SVD

After the signals are decomposed into several modes by VMD, the features are extracted from the modes by SVD, and can be constructed as a feature matrix. SVD is a powerful tool for feature extraction in linear algebra. According to SVD, the matrix could be decomposed as follow:
X = U ω V T
where X represents a m × n matrix. There are two orthogonal matrices: matrix U ( m × m ) and V ( n × n ) , and a singular diagonal matrix ω ( ω i j 0 , i = j and ω 11 ω 22 0 ), the diagonal element ω 11 , ω 22 , , ω m m is the singular value of X. U is called left singular vector, and the columns of U . V is called right singular vector. The columns of U and V especially are orthogonal to each other, and are base vector. To obtain more intrinsic information in the matrix, the singular vectors are selected. As a consequence, SVD is employed to decompose the eigenmatrix, and obtained the singular value vectors ( ω 11 , ω 22 , , ω m m ) .

2.3. Mahalanobis–Taguchi System

Mahalanobis–Taguchi System is a pattern recognition method integrated by the MD, orthogonal table and other tools such as SNR that are proposed by Taguchi [31], who introduces the experimental design of the field SNR to pattern recognition, which can reduce the dimensions of data, and use the orthogonal table to construct the MS. The MSs are used to calculate the MDs of the experimental data, and the valid features are distinguished by the SNR. Then, the MSs are recalculated by using the valid features. Finally, the results are obtained. The calculation of the MD is described as follows:

2.3.1. Mahalanobis Distance

The MD is a method of using normal data to normalize the fault data to compute the average distance between points and groups using normal data, the calculation formula of MD is as follows:
M D j = 1 k Z i j T C 1 Z i j
Z i j = { z 1 j , z 2 j , z 3 j , , z k j } , z i j = x i j x ¯ i s i
M D j represents the MD of the j th sample, k represents the number of the feature, x i j represents the i th feature’s value of the of the j th sample x ¯ i represents the mean of i th feature, s i represents the Standard deviation of i th feature, C 1 represents the inverse matrix of the correlation coefficient matrix.

2.3.2. Taguchi Method

In the Mahalanobis–Taguchi System, the MD measures the deviation of the test value from the normal value. The Taguchi method is able to select the features which have a larger contribution to identifying bearing faults, and then use the selected valid features to calculate the MD; the method is as below:

Orthogonal Array

Selecting the appropriate two-level orthogonal array, and then the k -original features obtained by VMD and SVD, are arranged into each column of the orthogonal array. In the orthogonal array, “1” indicates that the feature is selected, “2” indicates that the feature is not selected, and a MS is generated according to each row of the orthogonal array.

SNR and Its Gain

The main function of signal noise ratio is to select a valid feature, the calculation formula of generating SNR η i in the i line based on the orthogonal array.
η i = 10   lg 1 N j = 1 N M D i j
j [ 1 , N ] represents the number of training samples.
η i represents the recognition effects of the characteristic feature, the valid feature is selected by comparing the mean of SNR of each characteristic feature at two levels. The formula is as follow:
η j ¯ = η i m j = 1 , 2
Δ η j ¯ = ( η j ¯ m )
j = 1 , 2 represents two levels, i represents the number of rows in the MS, η 1 ¯ represents that in the level of ‘1’, the recognition effects to identify abnormal conditions of using this feature. η 2 ¯ represents that in the level of ‘2’, the recognition effects to identify abnormal conditions of not using this feature. If Δ η j ¯ represents the SNR gain, R = { Δ η j ¯ | η i ¯ > 0 } that indicates that this feature is a valid feature, if Δ η j ¯ < 0 that indicates that this feature is not a valid feature.

2.3.3. Adaptive Multiclass MTS

After SNR gain is calculated, to overcome the drawback of the conventional MTS during the threshold selection, this paper presents the adaptive multiclass MTS. There are several following improvements:
(1) Solving the multiple-classification problem. Selecting samples from each kind as the benchmark data. Then, the distances between other data and each benchmark data are calculated by MTS. Therefore, the label of the benchmark data with the minimum distance is selected as the label of the training data.
(2) Selecting the features adaptively. The feature sequence is selected several times, and the best fault recognition is the one which is minus MDs between it and benchmark. It solves the error problem caused by hard threshold selection.
(3) Avoiding the overfitting problem. Since the SNR gains are calculated by training samples, different training samples are employed to calculate the MDs and identify bearing faults, and the difference validation samples are set to validate the identified result and recalculate MDs.
This method is shown in Figure 2.
This adaptive multiclass MTS can be described as follow:
First, the data are labeled, the multi-class MSs are constructed, MDs between MSs are calculated and SNR gains are obtained. This step is shown as Figure 3; m is the number of samples, n is the number of features, and j is the number of the label. The training data are divided into three parts, and one of them is named as Benchmark. The data A j and C j represents one kind of data respectively (such as normal, fault of inner race, outer race, rolling element), and data B includes all kinds of data.
Second, new sequences of feature parameters are generated by the sequence of features’ SNR gains in ascending and descending order, then two collections are obtained from the above sequences based on the ascending or descending order, with the ascending and descending collection as follows:
r k = ( η k i ¯ , , η k j ¯ | η k i ¯ > η k j ¯ , i < j , i [ 1 , N ] )
q k = ( η k i ¯ , , η k j ¯ | η k i ¯ < η k j ¯ , i > j , i [ 1 , N ] )
R = { r k | k [ 1 , N ] }
Q = { q k | k [ 1 , N ] }
This step is shown in Figure 4:
Third, and the positions of feature’s SNR gain are the same as the positions of corresponding features in the sequences, the features that are used to recalculate the MDs are selected by the corresponding feature’s SNR gain. This step is shown in Figure 5. The MDs that are between each kind of A and C i are calculated.
Forth, the proper sequence of SNR gain is chosen by a function, which is according to the minimum MD. If two labels of data corresponding to minimum MD are the same, the recognition result is right, and accumulate the number of right recognition results. S is the recognition result. This step is shown in Figure 6.
Finally, the recognition result is verified if there is a unique optimal recognition result, and the SNR gain’s position in the sequence is the feature’s order. If there are several optimal recognition results, repeat step 3 to recalculate the result. Afterwards, a set of sequences with the best recognition effect is determined as the feature sequence, which solves the self-adaptation problem of thresholds.

3. Results

In this paper, the experimental data are from Case Western Reserve University Bearing Data Center. This experiment involved three different faults that occurred on three components: inner race, outer race and rolling element. The vibration signals were acquired under four different speeds: 1797 r/min, 1772 r/min, 1750 r/min, and 1730 r/min, and the sampling frequency was set to 12 kHz. To demonstrate the aMMTS, this study randomly selected the data in the dataset under the defect of 0.07 inches. The number of samples are shown in Table 1.
There were 2192 samples; 548 for inner race, 548 for outer race, 548 for rolling element and 548 for normal. The data are divided into three parts: training data, validation data and test data. In order to avoid the overfitting caused by the training data, training samples were used to construct MS, generate the SNR gain and calculate MDs by using the sequences of SNR gains, and were divided into three parts, with one of the parts set as benchmark group. To avoid the over-fitting problem, group A was used to construct MS and generate the SNR gain, and group B were used to calculate the MDs with the sequences of SNR gain and identify faults.
Validation samples were used to verify the recognition result if there exists the same minimum MDs, and the sequence was selected according to the best result.
Test samples were used to validate the proposed method.

3.1. Signal Decomposition by Using VMD and Wavelet Denosing

Above all, this study employed wavelet denoising to remove the noise from the raw signals. First, the Daubechies 5 (db5) was used to decompose the signal, and obtained the wavelet decomposition vector and the bookkeeping vector. Second, thresholds wavelet coefficient was calculated by setting the detail vector which would be compressed as [1,2,3] and the vector which is the corresponding percentages of lower coefficients as [100,90,80], and using the wavelet decomposition vector and the bookkeeping vector. Lastly, the thresholds, Daubechies 5 (db5) and decomposed signals were used to reconstruct the denoising signal. Then the VMD was used to decompose the signal, and was needed to give the preset IMF component number K and penalty parameter α which constrained the moderate bandwidth. The value of α toke the default value 1024, the value of K was 8. An example is shown in Figure 7.
IMFs are as shown in Figure 8.

3.2. Feature Extraction by Using SVD

SVD was used to analyze the IMFs. After the signal decomposition, the IMF matrix was decomposed by SVD, and obtained singular value vectors. The singular value vectors were considered as features and formed the feature matrix. Then, the feature matrix was used to diagnose the fault by aMMTS. To avoid the over-fitting problem, the features were divided into training samples, validation samples and test samples. The features of the above IMFs of those were shown in Table 2.
The features obtained by SVD are shown in Table 3.

3.3. Fault Diagnosis Using aMMTS

After the feature extraction, the aMMTS was used to identify and diagnose fault modes. The steps of aMMTS are as follow:
Firstly, the MS of training and benchmark were constructed, the eight-factor and two-level orthogonal array is shown in Table 4, and the MS based on Table 2 is shown in Table 5;
Secondly, the MD was calculated, and SNR gain was also obtained by benchmark samples and training samples. The SNR gain is shown in Table 6;
Thirdly, the MDs between the benchmark samples and validation samples were calculated by using the ascending and descending order of SNR;
Fourthly, the validation samples were used to verify the correctness of feature selection which existed more than one smallest MD;
Fifthly, the best sequence was chosen and set as the sequence of features.
Lastly, the best sequence was used to identify the test samples. Took the benchmark is outer race as the example shown in Figure 9.
Finally, the test sample was used to test the result of the method, and the benchmarks were inner race, rolling element, outer race and normal. The results are shown in Table 7 and the MDs between benchmark and test sample are shown in Figure 10.
As shown in Table 7 and Figure 10, this method accurately classified and diagnosed the fault of the bearing by using the different benchmarks. The recognition results of normal and outer race reached 100%. However, it is not accurate enough to diagnose the fault of inner race and rolling element. However, in the normal and the fault of inner race, it is effective in industrial application.

4. Discussion

Rolling element bearings are one of the most frequently used components in rotating machineries. This paper presents the method based on the wavelet denoising VMD-SVD-aMMTS to diagnose the fault of bearings under the variable conditions. Firstly, VMD is used to decompose the signal. Secondly, SVD is used to extract the feature. The adaptive aMMTS uses the feature sequences and multi-benchmarks to overcome the drawback of MTS for adaptive feature selection, multi-classification and over-fitting. The experimental result shows that the method could accurately diagnose faults effectively.
However, in the actual situation, there is an imbalance between fault data and normal data. In this method, aMMTS lacks research on the imbalanced study. The absence of faulty data may create a new problem, such over-fitting. Therefore, additional experiments under imbalanced data should be done to improve the method.

Author Contributions

N.W. collected and analyzed the data, made charts and diagrams, conceived and performed the experiments and wrote the paper; Z.W. and L.J. conceived the structure, provided guidance and modified the manuscript; X.C. analyzed the data and contributed analysis tools. Y.Q. provided guidance. Y.Z. revised the reviews.

Funding

This research was funded by National Key R&D Program of China grant number (No.2016YFB1200402).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tidriri, K.; Chatti, N.; Verron, S.; Tiplica, T. Bridging data-driven and model-based approaches for process fault diagnosis and health monitoring: A review of researches and future challenges. Annu. Rev. Control 2016, 42, 63–81. [Google Scholar] [CrossRef] [Green Version]
  2. Xin, G.; Hamzaoui, N.; Antoni, J. Semi-automated diagnosis of bearing faults based on a hidden Markov model of the vibration signals. Measurement 2018, 127, 141–166. [Google Scholar] [CrossRef]
  3. Al-Bugharbee, H.; Trendafilova, I. A fault diagnosis methodology for rolling element bearings based on advanced signal pretreatment and autoregressive modelling. J. Sound Vib. 2016, 369, 246–265. [Google Scholar] [CrossRef]
  4. Park, B.; Jeong, H.; Huh, H.; Kim, M.; Lee, S. Experimental study on the life prediction of servo motors through model-based system degradation assessment and accelerated degradation testing. J. Mech. Sci. Technol. 2018, 32, 5105–5110. [Google Scholar] [CrossRef]
  5. Jeong, H.; Park, B.; Park, S.; Min, H.; Lee, S. Fault detection and identification method using observer-based residuals. Reliab. Eng. Syst. Saf. 2018. [Google Scholar] [CrossRef]
  6. Cai, B.; Huang, L.; Xie, M. A data-driven fault diagnosis methodology in three-phase inverters for PMSM drive systems. IEEE Trans. Power Electron. 2017, 32, 5590–5600. [Google Scholar] [CrossRef]
  7. Cai, B.; Huang, L.; Xie, M. Bayesian Networks in Fault Diagnosis. IEEE Trans. Ind. Inform. 2017, 13, 2227–2240. [Google Scholar] [CrossRef]
  8. Cai, B.; Kong, X.; Liu, Y.; Lin, J.; Yuan, X.; Xu, H.; Ji, R. Application of Bayesian networks in reliability evaluation. IEEE Trans. Ind. Inform. 2018. [Google Scholar] [CrossRef]
  9. Chen, Z.; Deng, S.; Chen, X.; Li, C.; Sanchez, R.-V.; Qin, H. Deep neural networks-based rolling bearing fault diagnosis. Microelectron. Reliab. 2017, 75, 327–333. [Google Scholar] [CrossRef]
  10. Pandya, D.H.; Upadhyay, S.H.; Harsha, S.P. Fault diagnosis of rolling element bearing with intrinsic mode function of acoustic emission data using APF-KNN. Expert Syst. Appl. 2013, 40, 4137–4145. [Google Scholar] [CrossRef]
  11. Shao, H.; Jiang, H.; Lin, Y.; Li, X. A novel method for intelligent fault diagnosis of rolling bearings using ensemble deep auto-encoders. Mech. Syst. Signal Process. 2018, 102, 278–297. [Google Scholar] [CrossRef]
  12. Jia, F.; Lei, Y.; Guo, L.; Lin, J.; Xing, S. A neural network constructed by deep learning technique and its application to intelligent fault diagnosis of machines. Neurocomputing 2018, 272, 619–628. [Google Scholar] [CrossRef]
  13. Jagath Sri Lal Senanayaka, H.V.K.; Kjell, G. Robbersmyr Towards online bearing fault detection using envelope analysis of vibration signal and decision tree classification algorithm. In Proceedings of the 2017 20th International Conference on Electrical Machines and Systems (ICEMS), Sydney, Australia, 11–14 August 2017. [Google Scholar] [CrossRef]
  14. Rama Krishna, K.; Ramachandran, K.I. Machinery Bearing Fault Diagnosis Using Variational Mode Decomposition and Support Vector Machine as a Classifier. IOP Conf. Ser. Mater. Sci. Eng. 2018, 310, 012076. [Google Scholar] [CrossRef]
  15. Yan, R.; Gao, R.X.; Chen, X. Wavelets for fault diagnosis of rotary machines: A review with applications. Signal Process. 2014, 96, 1–15. [Google Scholar] [CrossRef]
  16. Ben Ali, J.; Fnaiech, N.; Saidi, L.; Chebel-Morello, B.; Fnaiech, F. Application of empirical mode decomposition and artificial neural network for automatic bearing fault diagnosis based on vibration signals. Appl. Acoust. 2015, 89, 16–27. [Google Scholar] [CrossRef]
  17. Muruganatham, B.; Sanjith, M.A.; Krishnakumar, B.; Satya Murty, S.A.V. Roller element bearing fault diagnosis using singular spectrum analysis. Mech. Syst. Signal Process. 2013, 35, 150–166. [Google Scholar] [CrossRef]
  18. Ke, F.; KeSheng, W.; Mian, Z.; Qing, N.; Ming, J.Z. A diagnostic signal selection scheme for planetary gearbox vibration monitoring under non-stationary operational conditions. Meas. Sci. Technol. 2017, 28, 035003. [Google Scholar]
  19. Caesarendra, W.; Tjahjowidodo, T. A Review of Feature Extraction Methods in Vibration-Based Condition Monitoring and Its Application for Degradation Trend Estimation of Low-Speed Slew Bearing. Machines 2017, 5, 21. [Google Scholar] [CrossRef]
  20. Francois Combet, R.Z. A new method for the estimation of the instantaneous speed relative fluctuation in a vibration signal based on the short time scale transform. Mech. Syst. Signal Process. 2009, 23, 16. [Google Scholar] [CrossRef]
  21. Cabrera, D.; Sancho, F.; Li, C.; Cerrada, M.; Sánchez, R.V.; Pacheco, F.; de Oliveira, J.V. Automatic feature extraction of time-series applied to fault severity assessment of helical gearbox in stationary and non-stationary speed operation. Appl. Soft Comput. 2017, 58, 53–64. [Google Scholar] [CrossRef]
  22. Chen, X.; Wang, Z.; Zhang, Z.; Jia, L.; Qin, Y. A Semi-Supervised Approach to Bearing Fault Diagnosis under Variable Conditions towards Imbalanced Unlabeled Data. Sensors 2018, 18, 2097. [Google Scholar] [CrossRef]
  23. Wang, Z.; Jia, L.; Kou, L.; Qin, Y. Spectral Kurtosis Entropy and Weighted SaE-ELM for Bogie Fault Diagnosis under Variable Conditions. Sensors 2018, 18, 1705. [Google Scholar] [CrossRef]
  24. Ciabattoni, L.; Ferracuti, F.; Freddi, A.; Monteriù, A. Statistical Spectral Analysis for Fault Diagnosis of Rotating Machines. IEEE Trans. Ind. Electron. 2018, 65, 4301–4310. [Google Scholar] [CrossRef]
  25. Chao-Ton Su, Y.-H.H. An Evaluation of the Robustness of MTS for Imbalanced Data. IEEE Trans. Knowl. Data Eng. 2007, 19, 12. [Google Scholar] [CrossRef]
  26. Cudney, E.A.; Hong, J.; Jugulum, R.; Paryani, K.; Ragsdell, K.M.; Taguchi, G. An Evaluation of Mahalanobis-Taguchi System and Neural Network for Multivariate Pattern Recognition. J. Ind. Syst. Eng. 2007, 1, 12. [Google Scholar]
  27. Soylemezoglu, A.; Jagannathan, S.; Saygin, C. Mahalanobis Taguchi System (MTS) as a Prognostics Tool for Rolling Element Bearing Failures. J. Manuf. Sci. Eng. 2010, 132. [Google Scholar] [CrossRef]
  28. Idriss El-Thalji, E.J. A summary of fault modelling and predictive health monitoring of rolling element bearings. Mech. Syst. Signal Process. 2015, 60, 252–272. [Google Scholar] [CrossRef]
  29. Dragomiretskiy, K.; Zosso, D. Variational Mode Decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  30. Wang, Y.; Markert, R.; Xiang, J.; Zheng, W. Research on variational mode decomposition and its application in detecting rub-impact fault of the rotor system. Mech. Syst. Signal Process. 2015, 60–61, 243–251. [Google Scholar] [CrossRef]
  31. Taguchi, G.; Konishi, S.; Konishi, S. Taguchi Methods: Orthogonal Arrays and Linear Graphs. Tools for Quality Engineering; American Supplier Institute: Dearborn, MI, USA, 1987. [Google Scholar]
Figure 1. The scheme of the proposed fault diagnosis method.
Figure 1. The scheme of the proposed fault diagnosis method.
Sensors 19 00026 g001
Figure 2. The step of Multiclass–Mahalanobis–Taguchi system (aMMTS).
Figure 2. The step of Multiclass–Mahalanobis–Taguchi system (aMMTS).
Sensors 19 00026 g002
Figure 3. The first step of aMMTS.
Figure 3. The first step of aMMTS.
Sensors 19 00026 g003
Figure 4. The second step of aMMTS.
Figure 4. The second step of aMMTS.
Sensors 19 00026 g004
Figure 5. The third step of aMMTS.
Figure 5. The third step of aMMTS.
Sensors 19 00026 g005
Figure 6. The forth step of aMMTS.
Figure 6. The forth step of aMMTS.
Sensors 19 00026 g006
Figure 7. The intercepted signal.
Figure 7. The intercepted signal.
Sensors 19 00026 g007
Figure 8. The IMFs.
Figure 8. The IMFs.
Sensors 19 00026 g008
Figure 9. The classification result of outer race.
Figure 9. The classification result of outer race.
Sensors 19 00026 g009
Figure 10. The Mahalanobis distances (MDs) between benchmark and testing data: (a) The MDs between benchmark (Inner race) and testing data; (b) The MDs between benchmark (Outer race) and testing data; (c) The MDs between benchmark (Rolling element) and testing data; (d) The MDs between benchmark (Normal) and testing data.
Figure 10. The Mahalanobis distances (MDs) between benchmark and testing data: (a) The MDs between benchmark (Inner race) and testing data; (b) The MDs between benchmark (Outer race) and testing data; (c) The MDs between benchmark (Rolling element) and testing data; (d) The MDs between benchmark (Normal) and testing data.
Sensors 19 00026 g010aSensors 19 00026 g010b
Table 1. The number of samples.
Table 1. The number of samples.
LabelMotor Load (hp)Speed (r/min)Training Samples Validation SamplesTest Samples
BenchMarkGroup AGroup B
Inner Race10 (0W)17972727292730
11 (735W)17722727292730
12 (1470W)17502727292730
13 (2205W)17302727292730
Outer Race20 (0W)17972727292730
21 (735W)17722727292730
22 (1470W)17502727292730
23 (2205W)17302727292730
Rolling Element30 (0W)17972727292730
31 (735W)17722727292730
32 (1470W)17502727292730
33 (2205W)17302727292730
Normal00 (0W)17972727292730
01 (735W)17722727292730
02 (1470W)17502727292730
03 (2205W)17302727292730
Table 2. The feature of IMFs.
Table 2. The feature of IMFs.
IMF1IMF2IMF3IMF4IMF5IMF6IMF7IMF8
Features2.0191.9681.8361.5821.4911.4591.2641.168
Table 3. The features of decomposed signals.
Table 3. The features of decomposed signals.
Features
12345678
Normal1.5541.1160.8480.4790.4120.3240.2140.095
1.4721.2300.4730.4250.2920.2580.1970.083
1.0941.0340.8970.4280.3210.2700.1850.082
1.1730.9390.9310.6830.3870.2870.2290.080
Inner Race2.0201.9681.8361.5821.4911.4591.2641.168
2.0341.9931.9161.6031.5831.3061.1950.913
2.1152.0231.9181.7631.5701.4101.3761.143
1.9411.8411.7981.5921.4821.4301.3931.151
Rolling Element0.8110.7930.6670.5770.5620.5420.4510.223
0.9320.7020.5850.5830.5200.5050.4640.408
0.7400.6570.5560.5150.5010.4870.4350.356
0.9680.7730.6860.6230.5910.5280.4760.467
Outer Race5.3494.2093.8333.7723.1461.9431.4691.163
5.3124.4684.0033.5052.9882.7511.4161.145
4.1413.2743.0242.9452.4401.7431.6311.019
3.3342.9012.5602.3822.0361.6971.2360.941
Table 4. The eight-factor and two-level orthogonal array.
Table 4. The eight-factor and two-level orthogonal array.
ABCDEFGH
111111111
211111222
311222111
412122122
512212212
612221221
721221122
821212221
921122212
Table 5. The Mahalanobis space (MS) based on the inner race.
Table 5. The Mahalanobis space (MS) based on the inner race.
ABCDEFGH
12.0201.9681.8361.5821.4911.4591.2641.168
22.0201.9681.8361.5821.491
32.0201.968 1.4591.2641.168
42.020 1.836 1.459
52.020 1.582 1.264
62.020 1.491 1.168
7 1.968 1.4911.459
8 1.968 1.582 1.168
9 1.9681.836 1.264
Table 6. The signal-to-noise ratio (SNR) gain of features.
Table 6. The signal-to-noise ratio (SNR) gain of features.
Features
12345678
Normal2.512−0.4710.3900.7260.789−0.0491.5191.469
Inner Race2.512−0.4710.3900.7260.789−0.0491.5191.469
Rolling element2.3970.1371.0070.4321.0780.4821.2830.770
Outer Race2.4390.9703.801−0.1040.6762.2161.501−1.809
Table 7. The recognition result.
Table 7. The recognition result.
Inner RaceOuter RaceRolling ElementNormalTotal
Result100%99.16%95%100%98.54%

Share and Cite

MDPI and ACS Style

Wang, N.; Wang, Z.; Jia, L.; Qin, Y.; Chen, X.; Zuo, Y. Adaptive Multiclass Mahalanobis Taguchi System for Bearing Fault Diagnosis under Variable Conditions. Sensors 2019, 19, 26. https://doi.org/10.3390/s19010026

AMA Style

Wang N, Wang Z, Jia L, Qin Y, Chen X, Zuo Y. Adaptive Multiclass Mahalanobis Taguchi System for Bearing Fault Diagnosis under Variable Conditions. Sensors. 2019; 19(1):26. https://doi.org/10.3390/s19010026

Chicago/Turabian Style

Wang, Ning, Zhipeng Wang, Limin Jia, Yong Qin, Xinan Chen, and Yakun Zuo. 2019. "Adaptive Multiclass Mahalanobis Taguchi System for Bearing Fault Diagnosis under Variable Conditions" Sensors 19, no. 1: 26. https://doi.org/10.3390/s19010026

APA Style

Wang, N., Wang, Z., Jia, L., Qin, Y., Chen, X., & Zuo, Y. (2019). Adaptive Multiclass Mahalanobis Taguchi System for Bearing Fault Diagnosis under Variable Conditions. Sensors, 19(1), 26. https://doi.org/10.3390/s19010026

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