Abstract

Methods based on vibration analysis are currently regarded as the most conclusive means for fault diagnosis and health prognostics in rotary machinery. However, changing working conditions mean that the vibration signals originating from rotary machinery exhibit different levels of complexity. This complexity leads to increased difficulty in constructing health indicators (HIs). In this paper, we propose a multiscale Tsallis permutation entropy (MTPE) to construct the HIs of rotary machinery under different working conditions. MTPE values are a function of an entropy index and scale, which have the universality for handling the complexity of a permutated time series. The health condition of the rotary machinery was effectively represented by the MTPEs in conditional monitoring; the initial point of the unhealthy stage was found using the interval. This was set as the alarm threshold according to the varying HI trend. Once this was established, dividing the stages into two-stage health stages (HS) was straightforward. Using a rolling bearing, a run-to-failure experiment was conducted and results suggested that the proposed method effectively assessed the status of the rotary machinery. Taken together, this study provided a novel complexity measure based on a methodology for constructing the HIs of rotary machinery and enriches conditional monitoring theory.

1. Introduction

Condition-based maintenance (CBM) is effective in reducing unnecessary maintenance operations and improving the reliability of rotary machinery [1]. One of the major tasks in CBM is health prognostics, which aims to predict the remaining useful life (RUL) of rotary machinery based on its historical and ongoing degradation trends. Extensive research has been done into data-driven methods for CBM; from this work, vibration analysis has emerged as the most conclusive method for fault diagnosis and health prognostics in rotary machinery [26]. The prognostic program generally consists of four technical processes: data acquisition, health indicator (HI) construction, health stage (HS) division, and RUL prediction [1]. As a representative of the health condition of the machinery, HI performance has a great influence on the effect of prognostics. In recent years, synthesized HIs have received much attention, owing to their impressive results for RUL prediction [7, 8]. These are usually constructed by a combination of signal processing techniques and artificial intelligent (AI) techniques [9]. Along with the measured vibration signals, both time- [10] and frequency-domain features [11] are widely used to analyze the status of the rotary machinery.

Conventional indicators, such as kurtosis [12] and skewness, fluctuate drastically and affect the assessment of the machinery running status [9]. Root mean square (RMS) [1316] values are stable during the healthy stage and will usually experience an “increase-decrease-increase” trend during the unhealthy stage. When only using the RMS values, the fault threshold is difficult to determine. With the development of the FFT algorithm, frequency-domain features are also utilized as criteria to identify the health condition of the rotary machinery. By identifying the aroused fault frequency by the fault across a wide frequency band, the defect is better visualized. However, in some cases (e.g., inner ring defect of bearing and complex defect), the complex spectrum makes the frequency-domain features less effective [9]. Moreover, an additional disadvantage is that the construction of frequency-domain features requires significant expertise.

AI techniques have also been used to establish HIs and assess the overall health conditions. These indicators are generally constructed by fusing multiple features, presenting a virtual description of the degradation trends of the machinery itself. For instance, Chen [17] constructed fused feature vectors using multiple, two-layer sparse autoencoder neural networks for bearing fault classification. Guo et al. [8] constructed a health index for bearing prognostics by fusing multiple features using a recurrent neural network. Using the input of an autoregression model to predict bearing states, Qian [18] extracted features by implementing a recurrence quantification analysis. The vibration features processed by the isometric feature mapping technique (ISOMAP) were proposed by Benkedjouh [19] to be the inputs of a support vector regression (SVR) model to predict the bearing RUL. Ocak et al. [20] fused multiple features using the hidden Markov model (HMM) in the healthy stage. Keshtegar et al. [21] developed the multi-extremum-modified response basis model for dynamic, nonlinear responses of failure capacities for a turbine blisk. Fei et al. [22] introduced an enhanced network learning method to improve the modeling precision and simulation efficiency in flexible mechanism reliability evaluation. Han et al. [23] proposed a probability-based prediction method to predict the service safety life of aero engine turbine blades. As indicated here, AI techniques have been widely used; however, past approaches have revealed that their shallow architectures restrict their capability of learning complicate nonlinear relationships [24, 25]. Moreover, given the highly nonlinear characteristics in the degradation process of rolling bearings [26], more powerful methodologies for the construction of HIs are needed.

During the degradation process of rotary machinery, increasingly complicated behaviors have been shown in the measured vibration signals used for health prognosis. The advent of better computational efficiency has allowed the complexity measure to calculate the complexity values in the vibration signals. This then allows for the assessment of the working status of the rotary machines. The complexity of a signal time series can be described by the Lempel-Ziv Complexity (LZC) [27, 28], Approximation Entropy (ApEn) [29], permutation entropy (PE) [3032], and multifeature entropy [33, 34]. These LZC values describe the degree of regularity of the vibration signals, and the quantitative relationship between the LZC values and the defect size, and the severity level. Moreover, this approach has already been established in rolling bearings. Despite this, using the definite calculation of the vibration signals based on the LZC index will lead in incorrect judgements for health conditions of rotary machinery in practical engineering. This is because the fault feature information is often submerged in noise [28]. For instance, Yan et al. [29] presented the ApEn to give the measurement of the regularity for mechanical structure signals. However, the ApEn greatly relies on data length; as a result, it will produce inaccurate estimation by using the ApEn to deal with a short time series. Compared to either LZC or ApEn, PE (multiscale PE) is a novel method used to quantify the degree of complexity and was originally proposed by Bandt and Pompe [30]. Fei et al. [33] introduced the multifeature entropy distance method for the process character analysis and subsequent diagnosis of rolling bearing faults. This was achieved by integrating four information entropies and two signal types. Moreover, the WCFSE-FSVM method [34] has also been proposed; it consists of fusing the advantages of the wavelet correlation feature scale entropy with those of the fuzzy support vector machine. However, both the PE and multiscale PE [35] have undeniable shortcoming [36]. More specifically, it has difficulty in handling the complexity of a permutated time series for rotary machinery under different working conditions.

In this paper, we proposed a MTPE analysis for conditional monitoring in the degradation process of rotary machinery. The MTPE values are a function of both the entropy index and scale, which have the universality for dealing with the complexity of a permutated time series. MTPE analysis has the capacity to construct HIs of the rotary machinery under different working conditions. The MTPE was used to construct these HIs; the effects of entropy index, scale, embedding dimension, and time delay for the time series from the standard logistic map [37] were also assessed. Using the interval as the alarm threshold, the initial point of the unhealthy stage was determined; the process was then divided into two health stages according to the varying trends of the HIs. For the unhealthy stage, the HI trend regarding the degradation process was established using a double exponential function together with an ARMAX model. A run-to-failure experiment for the rolling bearing was then designed to evaluate the capability of the proposed method. Results suggested that that the proposed method was effective for the assessment of the rotary machinery.

2. Performance Degradation Assessment Based on the MTPE Method

2.1. Multiscale Tsallis Permutation Entropy

Tsallis permutation entropy (TPE) is a generalization of the Boltzmann–Gibbs–Shannon entropy and has been proved to be efficient and suitable for multifractal systems [36]. The main difference between TPE and PE based on Shannon entropy is that the entropic index embedded in TPE is used to measure the nature of complex time series under different conditions [36]. The TPE is a special form of PE and derives from ordinal patterns. Given the original time series , the subtime series is constructed through the scale . Each subtime serieswhere denotes the round-off number of and is the scale. For , the time series is the original time series . The length of each subtime series is .

For any segment , the value of TPE is computed according to previously published work [31]. According to the Takens-Maine theorem, the phase space of time series is reconstructed aswhere is the embedded dimension and is the time delay. Then, the number of real values contained in each are arranged in increasing order as

If two or more elements in each have the same value, for example,, their original position is sorted for , and is recorded. Then, the vector is mapped onto a group of symbols aswhere and . , is one of the symbol permutations, if are used to denote the probability distribution of each symbol sequence, respectively. The PE of is defined based on the symbol permutations as

Then, the multiscale permutation entropy (MPE) for the original time series is given based on the Shannon entropy. However, the TPE is not simply the generalization of the Shannon permutation entropy. Given the entropy index , the MTPE is defined as follows:where is calculated based on the time series . Obviously, the value of MTPE in equation (6) is the function of the entropy index and the scale . The purpose of MTPE is to study the universality for dealing with the complexity of the permutated time series. Since its values are the function of the entropy index and scale, MTPE analysis has the capacity to construct the HIs of the rotary machinery under different working conditions.

2.2. Parameter Effects on MTPEs

To study the effects of different parameters on MTPE values, the logistic map shown in [37] was used as an exemplary object. The logistic map is often used to describe the change of quantity with increasing time, which is written as the simple mathematical model as follows:where , and is the control parameter. The bifurcate evolutional process of the logistic map with and chaotic behavior with are shown in [31, 37], respectively. In this paper, the time series using equation (7) is generated with 4096 points by starting from , where chaotic behavior occurred with . The effects of entropy index , embedding dimension , and time delay on the MTPEs are discussed with increasing scale in Figures 15. First, the MPE and MTPE were compared in Figure 1. Given these results, it was concluded that MTPE was more sensitive than MPE with changing scales.

The effects of the positive entropy index on MTPEs are shown in Figure 2; the effects of the negative entropy index are shown in Figure 3. All MPTEs and MPEs first increased, after which they slowly decreased with increasing scales. It is obvious that the values of the entropy index had a great impact on MTPE amplitude. All MTPEs had maximum amplitudes when using a scale of two; amplitudes then decreased with increasing scales larger than two. Other parameters that affected the MTPE calculations were the embedding dimension and the time delay . The relationship between the MTPE values and embedding dimension is illustrated in Figure 4, where the embedding dimensions were selected to be . Similarly, the relationship between the MTPE values and time delay is illustrated in Figure 5, where the time delays were selected to be . According to the analysis presented in previous work [31], we selected and to compute the MTPEs of the vibration signals in the degradation process of the rolling bearing.

2.3. Process of Assessing Performance Degradation

The bearing performance degradation assessment based on our MTPE analysis included several steps: (1) acquiring vibration signals, (2) calculating MTPEs, and (3) calculating HIs by normalization. The procedures for assessing the degradation in bearing performance are shown in Figure 6.

2.3.1. Vibration Signals Acquisition

During the degradation process of the rotary machinery, the time series of the vibration signals were collected at time points , where denotes the maximum number of vibration signals sampled in the process. Then, a set of vibration signals is collected to reflect the dynamic changes of the rotary machinery during the degradation process. The length of each is assumed to be and .

2.3.2. Calculation of MTPEs

The nonlinear dynamics of the rotary machinery are discriminated by computing the MTPEs of the time series of measured vibration signals. Assuming that the time series at the time is , the value of MTPE with entropy index and scale is calculated as according to equation (6). Then, the MTPE values for are given as the set .

2.3.3. Calculation of HIs by Normalization

The HIs are obtained by the normalized process for the set of MTPEs:where denotes the HI, ,. The HIs are used to represent the health condition during the degradation process of rolling bearing.

2.3.4. Healthy Stage Division

According to the varying degradation trend of HIs, the whole lifetime of the rolling bearing was divided into two or more different health stages. The HIs in the unhealthy stage were fitted using the analysis of degradation trends. A strategy for the two-stage division was used by identifying whether the HIs exceeded a constant alarm threshold.

3. Experimental Verification

3.1. Experimental System

A run-to-failure experiment for the rolling bearing was next designed and conducted to investigate the proposed method. The ABLT-7 bearing test system is shown in Figure 7, which consisted of a supporting pillow block, a test bearing with supporting housing, control electronics, and other features. The input-shaft of the test system is connected with an electromotor through a V-type belt, where the speed of the motor was adjusted using the integrated speed controller.

An accelerometer (B&K4507B) and an acoustic emission sensor (B&K4190L) were installed on the testing bearing housing to measure the vibrations. The resulting data were acquired through a multichannel analyzer (B&K3560C) and Pulse software. A virtual instrumentation system was established by VC to obtain the time-domain feature of the vibrations, thereby ensuring the safety of the testing system.

In the run to failure experiment, a UCP207 insert bearing was tested. The radial load went from 5040 N for the first two hours to 8020 N for the remainder of its life; rotational speed was 3000 rpm. The sample frequency was 20 kHz and the length of each vibration time series for sampling was set to 4096. Vibration signals were obtained during the experiment at an interval of every 5 min. In total, 530 data segments were collected to conduct the MTPE analysis and determine the status of the rolling bearing.

3.2. Vibration Signals Denoising and Analysis

As a widely used data denoising technique, we used wavelet transform to denoise the set of vibration signals [27]. The transform was a linear transform that analyzed nonstationary signals at varying resolutions by decomposing the signals into their frequency bands. In this study, the Symmlet family of wavelet was chosen as the mother wavelet to implement the denoising. During the process of wavelet decomposition, the soft threshold was used and a denoised version of the signal was obtained using a wavelet coefficients threshold with a global positive threshold. The effects of denoising by Symmlet 4 wavelet and db 1 wavelet with level 5 for a set of vibration signals are compared in Figure 8.

To verify the degradation trend of rolling bearing, several collected vibration signals (Figure 9) at different sampling times were transformed to frequency domains using Fast Fourier Transform (FFT). The frequency spectrums of the vibration signals are shown in Figure 10. During the degradation process of the rolling bearing, more frequency components were aroused and the amplitudes of these frequencies were notably increased. Moreover, the dynamics of the rolling bearing during the degradation process became more complicated and nonlinear. In particular, the last vibration signal had a great deal of harmonics in its frequency spectrum.

3.3. Healthy Indicators Construction

In this section, the health indicators of the rolling bearing were constructed based on the MTPEs of the vibration signals obtained during the process of degradation. For the measured vibration signals , a suitable entropy index and scale were determined to calculate the MTPEs.

To investigate the influence of these parameters, the relationship between health indicators and the entropy index is given in Figure 11. The relationship between health indicators and scale is shown in Figure 12. In the following discussion, entropy index and scale were chosen to assess the status characterization for the rolling bearing during the degradation process. The embedding dimension and time delay also had significant influences on the MTPE values during the degradation process of the rolling bearing. The results of the health indicators on these different parameters are shown in Figures 13 and 14, respectively.

Health indicators’ trends were greatly changed with an embedding dimension of (Figure 12) and a time delay of (Figure 14). Based on the above discussion, the embedding dimension and time delay were determined to be and . These results were similar to a previous analysis performed by permutation entropy [31]. The time-domain features (e.g., kurtosis, RMS, and skewness) were also computed and compared with the proposed health indicators in Figure 15. The kurtosis and skewness of the measured vibration signals fluctuated dramatically during the process of degradation for the rolling bearing; the RMS values were stable during this process.

However, the values of RMS decreased to a certain extent at the unhealthy stage of the rolling bearing. This gave an unreliable estimate of the health condition of the test bearing. The performance of prognostic HIs had a great influence on the complexity of prognostic modeling and the prediction accuracy. It should be noted that while the proposed health indicators seriously fluctuated at certain stages, the overall trend of the proposed HIs was easily divided into healthy and unhealthy stages.

3.4. Healthy Stage Division

With the varying degradation trends exhibited by the HIs, the degradation processes of the machinery should be divided into different HSs before RUL prediction. The aim of HS division is to divide the degradation processes of rolling processes into different HSs according to the varying trends of the HIs. Generally, the degradation trend of a rolling element bearing has two distinct stages—healthy and unhealthy [1, 8, 38].

During the degradation process of the bearing, the proposed health indicators exhibited trends that were divided into two stages. A strategy for this two-stage division was used, whereby the HIs that exceeded a constant alarm threshold were identified [1]. Past work has used both constant [3941] and adaptive [42, 43] alarm thresholds to divide into healthy and unhealthy stages. This paper used the interval as the alarm threshold to determine the initial starting point of the unhealthy stage; results are compared with the RMS values in Figure 16. Since no fault occurred in the rolling element bearing, the values of the proposed healthy indicators presented random fluctuation during the healthy stage. Since there was no information regarding the degradation trend in the healthy stage, it was difficult and unnecessary to predict the RUL during this stage. However, HIs increased quickly with the deterioration of the tested rolling bearing at the unhealthy stage. However, the RMS values shown in Figure 16 indicate they were stable during the healthy stage. After that, they experienced an ‘‘increase-decrease-increase” trend during the unhealthy stage. Using the interval as the alarm threshold to determine the initial point of the unhealthy stage, two HIs increased rapidly and reached the summit in a short time. However, the RMS values decreased during the last part of the unhealthy stage and this trend was difficult to model for RUI prediction.

4. Conclusions

This paper proposed the MTPE method for conditional monitoring during the degradation process of rotary machinery. The HIs were calculated by employing the MTPE to assess the complexity of the short and noising vibration signals. Using the interval as the alarm threshold, two stages were identified—healthy and unhealthy. Since the values of MTPE were the function of the entropy index and scale, this method avoided the difficulties of constructing health indicators under changing working conditions. With the universality for assessing the complexity of permutated time series, this study provided a methodology using a novel complexity measure for conditional monitoring of rotary machinery.

Data Availability

The data and the MATLAB programs used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the National Natural Science Foundation of China (Grant no. 51775182) and Guangdong Provincial Natural Science Foundation of China (2020B1515120006).