Secciones
Referencias
Resumen
Servicios
Descargas
HTML
ePub
PDF
Buscar
Fuente


Use of Multiscale Permutation Entropy Feature Selection and Supervised Classifiers for Bearing Failures Diagnosis
Uso de la Entropía con Permutación Multiescalar, técnicas de Selección de Características y Clasificadores Supervisados para el Diagnóstico de Fallas en Rodamientos
Scientia Et Technica, vol. 26, núm. 4, pp. 443-448, 2021
Universidad Tecnológica de Pereira

Mecánica



Recepción: 30 Octubre 2020

Aprobación: 20 Noviembre 2021

DOI: https://doi.org/10.22517/23447214.24579

Abstract: Entropy measurements are an accessible tool to perform irregularity and uncertainty measurements present in time series. In signal processing, the Multiscale Permutation Entropy is recently presented as a methodology of characterization capable of measuring randomness and non-linear dynamics existing in non-stationary time series, such as mechanical vibration signals. In this article, the Multiscale Permutation Entropy is combined with diverse feature selection techniques and multiple classifiers based on machine learning aiming to detect different operative states in an internal combustion engine. The best combination is selected from the evaluation of parameters like precision and computational time. Finally, the proposed methodology is established as an effective tool to diagnose failures in bearing systems with a high precision rate and a reduced calculation time.

Keywords: Dynamics, Entropy, Machine, Multiescale, Permutation, Vibration.

Resumen: Las mediciones de entropía son una herramienta accesible para realizar mediciones de irregularidades e incertidumbres presentes en series de tiempo. En el procesamiento de señales, la Entropía de Permutación Multiescalar se presenta recientemente como una metodología de caracterización capaz de medir la aleatoriedad y la dinámica no lineal existente en series de tiempo no estacionarias, como las señales de vibración mecánica. En este artículo, la entropía de permutación multiescalar se combina con diversas técnicas de selección de características y múltiples clasificadores basados ​​en el aprendizaje automático con el objetivo de detectar diferentes estados operativos en un motor de combustión interna. La mejor combinación se selecciona a partir de la evaluación de

parámetros como precisión y tiempo computacional. Finalmente, la metodología propuesta se establece como una herramienta eficaz para diagnosticar fallas en sistemas de rodamientos con una alta tasa de precisión y un tiempo de cálculo reducido.

Palabras clave: Dinámica, Entropía, Máquina, Multiescalar, Permutación, Vibración..

I. INTRODUCTION

Bearing systems play an important role in rotary machines and in the modern manufacturing industry. Different methodologies have been developed for the detection and diagnosis of faults in its main components [1]. Generally, these diagnoses are made from the capture and processing of vibration signals, since they contain relevant information about the state of the machine [2]. However, these signals have a large number of non-stationary and non-linear characteristics, since their capture inevitably takes place with friction and impacts. To overcome this problem, a series of techniques have been developed for processing and classifying these signals. A widely used approach is based on the analysis of temporal and spectral characteristics of vibration signals [3][4]. However, analysis in time, frequency and time-frequency domains are seriously affected by the signal length and sampling frequency of the capture [5]. To solve this problem another approach has been presented in recent years, which is based on entropy of different natures, such as Simple Entropy (ApEn) [6], Approximate Entropy (SampEn) [7], Multiscale Entropy (MSE) [8], Permutation Entropy (PE) [9][31] and Multiscale Permutation Entropy (MPE) [10]. All of the above mentioned have been used successfully for the characterization of signals of different nature. For instance, in [6] the ApEn is used for diagnosis and clinical monitoring in the area of physiology. In [11] PE is used for the classification of patients from EEG signals. Finally, in [1][8][32] is used the MPE and MSE for the identification and diagnosis of faults in bearing systems based on vibration signals. It should be noted that the MPE is an evolution of the MSE and the PE since it gives a much more complete measure of the nonlinear dynamic parameters of a system. MPE includes a combination of different scales and time delays, which identifies particularities that are not perceptible within other entropies [10][33]. All the previous methods of characterization allow to obtain a great quantity of information of the system that is being analyzed. When classifying, many of these characteristics can be redundant or irrelevant, which can create redundancy or over-training of the classifier, for which characteristics selection techniques are implemented, improving the quality and efficiency of the model. In the bearing fault diagnosis application, Variance based on Relevance Analysis (VRA) [1], Laplacian Score (LS) [12] and Relief (REL) [5] are normally used. After the selection of features, a classification process is performed with machine learning algorithms such as Multiple Vector Support Machines (SVMM) [1],[13], [14] and [15], Hidden Markov Models (HMM) [16] and Artificial Neural Networks (ANN) [17], [18] and [19]. However, these classifiers have a high degree of complexity, computation time and initial parameters that must be optimized. Few works have attempted to exploit the potential of less complex conventional classifiers, such as Nearest K-Neighbors (KNN) [20], [21], Decision Trees (TREE) [22] or Naive Bayes (BAYES) [21]. This paper proposes a methodology for the diagnosis of bearing failures based on characteristics of the MPE. To make it more efficient and effective, it is combined with some feature selection techniques and multiple supervised classifiers. The verification of the advantage of the chosen parameters is done by comparison with different ones used in the state of the art.

The rest of the article is organized as follows: In section II.A is presented the mathematical formulation of the MPE and in section II.B is detailed the different feature selection techniques used. In section II.C the classifier used in the state of the art are exposed. Then, a novel methodology for the diagnosis of bearing failures based on MPE, feature selection techniques and supervised classifier is located in the section III. Finally, results are presented in section IV and the conclusions of this approach in section V.

II. Methodology

The proposed methodology combines a characterization method, a feature selection technique and a supervised classifier as show in Fig. 1. Each part of the methodology is described in the following sections. The proposed methodology begins with the characterization of vibration signals, then using an automatic classifier combined with a characteristic selection technique, the state of the system is estimated.


Fig 1
Methodology for vibration classification

A. Multiscale Permutation Entropy and Acronyms

Multiscale Permutation Entropy (MPE) is used in this paper as a signal characterization method, since it is a measure that allows to detect dynamic changes in the time series. It is based on the comparison in neighboring values without taking into account the size of the values and, therefore, has a calculation simple and fast [23]. The above, allows to position the MPE as a particularly useful and robust tool in the presence of dynamic noise [1]. In order to describe the MPE measure, it is important to review the entropy proposed by Shannon [14], described as follows: Considering a time series in a space of representation of characteristics, where T is the length of the time series, the entropy is represented as in (1):

(1)

Where xi ∈ R and p(xi) is the marginal probability. The time series can be represented with a delay of time and dimension given by (2):

(2)

Where i = 1, 2,...,T − (m−1)τ, m is the dimension and τ the delay. To perform the computation of MPE, the signal must be truncated in N = T − (m − 1)τ subvectors. Then, for each subvector is calculated the entropy mapped in a space of m! different symbols denoted as by (3):

(3)

The probability is calculated by (4):

(4)

Where the judgment type denotes the map from pattern space to symbol space. Also, 1A(u) = 1 if u ∈ and 1A(u) = 0 if u A. The MPE can take values between the ranges [0, ln(m!)] and it is invariant under nonlinear monotonic transformations. The values of m and τ vary from 1 to 8 and the values, which are values used for calculating MSE [8] and PE [9]. In Fig 2 is plotted the behavior of the MPE, when are varying m and τ for a vibration signal.


Fig 2.
MPE varying m and τ

When choosing the delay and dimension parameters, the nature of the signals must be taken into account. If the parameters are too small, the nonlinear dynamic of the characteristics from signals will not be analyzed effectively. If the parameters are too large, useful information will be deleted in consequence, which will result in a wicked analysis.

B. Featured Selection Technique

For the comparison of the characteristic selection technique, the ones most used in the literature were implemented.

· Relief (REL). Unsupervised method to generate a ranking based on the predictors give neighbors of the same class or different class.

· Laplacian Score (LS). Unsupervised method to generate a ranking based on the input characteristics based on a variability criterion [1].

· Variance-based Relevance Analysis (VRA). Unsupervised method, it generates a ranking with the input features based on a variability criterion [12].

· Non-negative matrix factorization (NF). Unsupervised method to analyze the relevance of the characteristics from a reduction in dimensionality by non-negative factorization techniques [20].

· Self-Weight Ranking (SW). Unsupervised method, technique, it codes the feature relevance in terms of a self-similarity measures [5].

· Distance-Weight Ranking (DW). Supervised method, this technique quantifies the relevance of the distance between samples from different clusters by using supervised information.

C. Supervised Classifiers

Characteristic selection techniques are combined with different supervised classifiers to find the best combination that suits the nature of the signal.

· K-Nearest Neighbors (KNN). Supervised and non-parametric classification method that estimates the posterior probability that an element belongs to the class from a set of information provided [24].

· Decision Trees (TREE). Supervised classifier based on prediction systems based on rules, which serve to represent and categorize a series of conditions that occur successively, for the resolution of a problem [22].

· Naive bayes (BAYES). Supervised classifier for multiclass learning. This classifier is based on estimation of prediction and re-substitution [21].

· Multiple Support Vector Machines (SVMM). Supervised classifier based on a hyperplane or set of hyperplanes in a space of very high (or even infinite) dimensionality that can be used in problems of classification or regression [1], [13], [15] and [14].

III. experimental data

The validation of the proposed methodology is carried out by evaluation performed when classifying bearing fault signals obtained from the Case Western Reserve database [25]. In this database, signals were collected for the normal bearings (Nor), faults in the internal train (IR1), external train (IR2) and ball (BE). Faults are also found in order of severity, 0.007 inches in diameter to 0.040 inches in diameter and at variable engine speeds of 1720 to 1797 RPM. Each experiment was repeated three times and the data was collected at 12 kHz for 5 seconds. Each signal was divided into 10 sub-signals in order to have more samples per class and imitating the experimental framework established in the literature [26]. A sub-signal of each of the faults can be seen in Fig.3.


Fig 3.
Vibration signals for different states

The signals are characterized with the MPE and then dividing it into training and testing. A feature selection technique is applied to the training group to obtain a ranking of relevance and it is applied to the test group. Then, with the reorganized training group, a classifier is trained and the test group is evaluated. The evaluation process is explained in Fig. 4.


Fig 4.
. Methodology proposed for the detection of faults in bearings systems

IV. results and discussion

This paper exposes a methodology for classify bearing fault in vibration signal. The classification carried out through a cross validation of (k = 5) and repeated by changing the number of training characteristics. After the validation, the best results were chosen through the shortest distance to the ideal point [100% (Acc); 0 (Std); 0 (Ca)], where Acc is the accuracy, Std is the standard deviation of the accuracy of the cross validation and Ca is the number of characteristics. In tables I and II is exposed the accuracy with a 95% confidence interval for the combinations of classifiers and feature selection techniques for 4 (Nor, IR1, IR2 and BE) and 10 (Nor, IR1, IR2 and BE, combined with different speeds) classes are appreciated.

Table I.
Accuracy Of The Classification Of 4 Classes

TABLE II.
ACCURACY OF THE CLASSIFICATION OF 10 CLASSES

The best results are obtained with the KNN classifier, regardless of the feature selection technique. Specifically, the best accuracy was obtained with the KNN classifier combined with the Relief feature selection technique. The quantity of characteristics used for the classification were 9 and 10 for 4 and 10 classes respectively. It should be noted that the results are obtained thanks to the characterization made with the MPE, which achieves a high level of separability of the classes that allows the classifiers to adapt and solve the proposed application. Finally, a summary of the best classifications can be seen in the table III.

TABLE III
COMPARISON OF THE BEST CLASSIFICATION RESULTS.

The proposed experimental methodology achieved a mean accuracy of 99.72% for classification of 4 different bearing system failures. Table IV shows the confusion matrix for this classifier.

TABLE IV.
CONFUCTION MATRIX OF THE BEST RANKING.

The proposed methodology is capable of classifying the failures with a high success rate, to the point that only one sample avoided obtaining 100% accuracy. The effectiveness of the methodology can be seen with also achieved a mean accuracy of 99.55% for the classification of 10 different bearing system failures.

V. conclusions

This article presents a methodology for the diagnosis of bearing failures based on the Multiscalar Permutation Entropy (MPE) technique. The MPE proves to be a highly effective characterization methodology to find information that allows to separate classes. Specifically, in the mechanical vibration signals that have a high non-stationary, the MPE manages to find characteristics that would not be detected by other methodologies. The MPE measures the non-linear dynamics existing in non-stationary time series and when combined with Relieff as a feature selection technique, a robust tool for classification applications is obtained. For the classification a method of K-Neighbors Nearest (KNN) was used, which manages to adapt to the nature of the characteristics. The results confirm a classification accuracy of more than 99:9% with a computation time of 16:37 seconds, which exceeds the results established in the literature.

Acknowledgments

The authors thank the Master in Electrical and Mechanical Engineering of the Technological University of Pereira for their support throughout the investigation. In addition, we would like to thank COLCIENCIAS for supporting the project entitled: “Desarrollo de un sistema de monitoreo para el análisis energético y de condición de emisiones en motores de combustión interna diésel con base en técnicas no destructivas” with code 1110-776-57801, through which the research described in this article was developed.

References

[1]. P Zheng, and H Yang, “Generalized composite multiscale permutation entropy and laplacian score based rolling bearing fault diagnosis,” Mech. Syst. Signal Process., vol. 99, pp. 229–243, 2018. DOI: https://doi.org/10.1016/j.ymssp.2017.06.011

[2]. R Tiwari, V Gupta, and P. K Kankar, “Bearing fault diagnosis based on multi-scale permutation entropy and adaptive neuro fuzzy classifier,” J. of Vib. Control, vol. 21, no 3, pp. 461–467, 2015. DOI: https://doi.org/10.1177/1077546313490778

[3]. H. Martin, and F. Honarvar, “Application of statistical moments to bearing failure detection,” Appl. Acoust, vol. 44, pp. 67–77, 1995. DOI: https://doi.org/10.1016/0003-682X(94)P4420-B

[4]. R. B. W. Heng, and M. J. M. Nor, “Statistical analysis of sound and vibration signals for monitoring rolling element bearing condition,” Appl. Acoust., vol. 53, 211–226, 1998. DOI: https://doi.org/10.1016/S0003-682X(97)00018-2

[5]. Z. Wei, Y. Wang, and H. Shuilong, “A novel intelligent method for bearing fault diagnosis based on affinity propagation clustering and adaptive feature selection,” Knowledge-Based Systems, vol. 116, pp. 1–12, 2017. DOI: https://doi.org/10.1016/j.knosys.2016.10.022

[6]. J. S. Richman, and J. R. Moorman, “Physiological time-series analysis using approximate entropy and sample entropy,” American J. of Physiology Heart and Circulatory Physiology, vol. 278: pp. 39–49, 2000. DOI: https://doi.org/10.1152/ajpheart.2000.278.6.H2039

[7]. M. Pincus, “Approximate entropy as a measure of system complexity,” Proceedings of the National Academy of Sciences, vol. 88, no 6, pp. 2297–2301, 1991. DOI: 10.1073/pnas.88.6.2297

[8]. M. Costa, A. Goldberger, and P. Chung, “Multiscale entropy analysis of physiologic time series,” Physical Review Letters, vol. 89, pp. 220–232, 2001. DOI: 10.1103/PhysRevLett.89.068102

[9]. C. Bandt, and B. Pompe, “Permutation entropy: A natural complexity measure for time series,” Physical Review Letters, vol. 88, pp. 74–102, 2002. DOI: 10.1103/PhysRevLett.88.174102

[10]. W. Aziz, and M. Arif, “Multiscale permutation entropy of physiological time series,” 2005 Pakistan Section Multitopic Conference, Karachi 2005, vol. 53, pp.1–6, 2005. DOI: https://doi.org/10.1109/ inmic.2005.334494

[11]. K. Zeng, and G. Ouyang, “Characterizing dynamics of absence seizure EEG with spatial-temporal permutation entropy,” Neurocomputing, vol. 275, pp. 577 – 585, 2018. DOI: https://doi.org/10.1016/j.neucom.2017.09.007

[12]. S. Daza, J. Arias, and J. Llorente, “Dynamic feature extraction: an application to voice pathology detection,” Intelligent Automation and Soft Computing, vol. 15, pp. 667–682, 2009. DOI: 10.1080/10798587.2009.10643056

[13]. X. Zhang, Y. Liang, J. Zhou, and Y. Zang, “A novel bearing fault diagnosis model integrated permutation entropy, ensemble empirical mode decomposition and optimized SVM,” Measurement, vol. 69, pp. 164–179, 2015. DOI: https://doi.org/10.1016/j.measurement.2015.03.017

[14]. J. Zheng, J. Cheng, and Y. Yang, “A rolling bearing fault diagnosis approach based on LCD and fuzzy entrop,” Mech. and Mach. Theory, vol. 70, pp. 441–453, 2013. DOIO: https://doi.org/10.1016/j.mechmachtheory.2013.08.014

[15]. H. Minghong, and P. Jiali, “A fault diagnosis method combined with LMD, sample entropy and energy ratio for roller bearings,” Measurement, vol. 76, pp. 7–19, 2015. DOI: https://doi.org/10.1016/j.measurement.2015.08.019

[16]. M. Yuwono, Y. Qin, and J. Zhou, “Automatic bearing fault diagnosis using particle swarm clustering and hidden Markov model,” Engin. Appl. of Artif. Intell., vol. 47: pp. 88–100, 2016. DOI: https://doi.org/10.1016/j.engappai.2015.03.007

[17]. B: Muruganatham, S. Krishnakumar, and S. Murty, “Roller element bearing fault diagnosis using singular spectrum analysis. Mech. Syst. Signal Process., vol. 35, no 1, pp. 150–166, 2013. DOI: https://doi.org/10.1016/j.ymssp.2012.08.019

[18]. B. Jaouherl, N. Fnaiech, and L. Saidi, “Application of empirical mode decomposition and artificial neural network for automatic bearing fault diagnosis based on vibration signals,” Applied Acoustics, 2015, vol. 89, pp. 16–27, 2015. DOI: https://doi.org/10.1016/j.apacoust.2014.08.016

[19]. P. William, and M. Hoffman, “Identification of bearing faults using time domain zero-crossings,” Mech. Syst. Signal Process., 2011, vol. 25, no 8, pp. 3078–3088, 2011. DOI: https://doi.org/10.1016/j.ymssp.2011.06.001

[20]. L. Liang, F. Liu, and M. Li, “Feature selection for machine fault diagnosis using clustering of non-negation matrix factorization,” Measurement, vol. 94, pp. 295–305, 2016. DOI: https://doi.org/10.1016/j.measurement.2016.08.003

[21]. O. Vencalek, O. Pokotylo, “Depth-weighted bayes classification,” Computational Statistics Data Analysis, vol. 123, pp. 1–12, 2018. DOI: https://doi.org/10.1016/j.csda.2018.01.011

[22]. B. Penny, and J. Geoff, “Classification trees for poverty mapping, “ Computational Statistics Data Analysis, vol. 115, pp. 53 – 66, 2017. DOI: https://doi.org/10.1016/j.csda.2017.05.009

[23]. C. Yinhe, T. Wen, and G. Jue, “Detecting dynamical changes in time series using the permutation entropy,” Phys. Rev. E, vol. 70, pp. 174–195, 2004. DOI: https://doi.org/10.1103/PhysRevE.70.046217

[24]. M. Philipp, S. Katri, and N. Ville, “Scent classification by k nearest neighbors using ion-mobility spectrometry measurements,” Expert Systems with Appl., vol. 15, pp. 199–210, 2018. DOI: https://doi.org/10.1016/j.eswa.2018.08.042

[25]. Bearing Data Center. https://csegroups.case.edu/bearingdatacenter/home

[26]. R. Nayana, and P. Geethanjali, “Analysis of statistical time-domain features effectiveness in identification of bearing faults from vibration signal,” IEEE Sensors J., vol., 17, pp. 5618–5625, 2017. DOI: 10.1109/JSEN.2017.2727638

[27]. Z. Keheng, S. Xigeng, and X. Dongxin, “A roller bearing fault diagnosis method based on hierarchical entropy and support vector machine with particle swarm optimization algorithm,” Measurement, vol. 47, pp. 669–675, 2014. DOI: https://doi.org/10.1016/j.measurement.2013.09.019

[28]. L. Zhiwen, C. Hongrui, and C. Xuefeng, “Multi-fault classification based on wavelet SVM with PSO algorithm to analyze vibration signals from rolling element bearings,” Neurocomputing, vol. 99, pp. 399–410, 2013. DOI: https://doi.org/10.1016/j.neucom.2012.07.019

[29]. H. Ocak, and K. Loparo, “HMM-based fault detection and diagnosis scheme for rolling element bearings,” J. of Vib. Acoust., vol. 127, no 4, pp. 299–306, 2005. DOI: https://doi.org/10.1115/1.1924636

[30]. H. Shao, H. Jiang, F. Wang, and H. Zhao, “An enhancement deep feature fusion method for rotating machinery fault diagnosis,” Knowledge-Based Systems, vol. 119, pp. 200–220, 2017. DOI: https://doi.org/10.1016/j.knosys.2016.12.012

[31]. M. Cristina and H. Anne, “Multiscale permutation entropy for two-dimensional patterns”, Pattern Recognition Letters, vol 150, pp. 139-146. 2021. DOI https://doi.org/10.1016/j.patrec.2021.06.028.

[32]. Z. Mei, Z. Zhang and L. Xie, “Permutation entropy-based detection scheme of replay attacks in industrial cyber-physical systems”, Journal of the Franklin Institute, vol. 358, pp. 4058-4076, 2021. DOI https://doi.org/10.1016/j.jfranklin.2021.02.024.

[33]. G. Vashishtha, S. Chauhan, M. Singh and R. Kumar, “Bearing defect identification by swarm decomposition considering permutation entropy measure and opposition-based slime mould algorithm”, Measurement, vol. 178, pp. 109-389, 2021, DOI https://doi.org/10.1016/j.measurement.2021.109389.

Notas de autor

J. C. Mejía

J. C. Mejía
Juan C. Mejia

Electronic Engineer of the University of Quindío. Professor at the Technological University of Pereira. He has extensive experience as a researcher in fields of instrumentation, digital systems programming, signal processing and machine learning algorithms. He has multiple bibliographic productions in indexed journals and a wide variety of participation in scientific events. He is part of the research group Manufacturing Processes and Machine Design attached to the Technological University of Pereira.

H. F. Quintero

H. F. Quintero
Héctor Fabio Quintero Riaza

received Ph.D. degree in Mechanical Engineering from Cataluña Polytechnic University, Barcelona, Spain, in 2006. Now he works at Technological University of Pereira. His current research interests include analysis and synthesis of mechanisms, mechanical vibrations, and fault diagnosis. He is part of the research group Manufacturing Processes and Machine Design attached to the Technological University of Pereira.

J. D. Echeverry-Correa

J. D. Echeverry-Correa,
Julian D. Echeverry-Correa

received his PhD in Electronic Systems from Universidad Politécnica de Madrid, Spain, in 2015. Since 2007 he is professor at the Program of Electrical Engineering at Universidad Tecnológica de Pereira, currently as an Associate Professor. His main research areas include signal processing, data analysis and pattern recognition.



Buscar:
Ir a la Página
IR
Visor de artículos científicos generados a partir de XML-JATS4R por