Abstract
In allusion to the issue of rolling bearing degradation feature extraction and degradation condition clustering, a logistic chaotic map is introduced to analyze the advantages of C0 complexity and a technique based on a multidimensional degradation feature and Gath–Geva fuzzy clustering algorithmic is proposed. The multidimensional degradation feature includes C0 complexity, root mean square, and curved time parameter which is more in line with the performance degradation process. Gath–Geva fuzzy clustering is introduced to divide different conditions during the degradation process. A rolling bearing lifetime vibration signal from intelligent maintenance system bearing test center was introduced for instance analysis. The results show that C0 complexity is able to describe the degradation process and has advantages in sensitivity and calculation speed. The introduced degradation indicator curved time parameter can reflect the agglomeration character of the degradation condition at time dimension, which is more in line with the performance degradation pattern of mechanical equipment. The Gath–Geva fuzzy clustering algorithmic is able to cluster degradation condition of mechanical equipment such as bearings accurately.
Highlights
Performance degradation feature extraction based on C0 is proposed. A Curved time parameter is proposed for describing the performance degradation process. Sequence dispersion is proposed to comprehensively evaluate the clustering effect. The Gath–Geva method is introduced in performance degradation condition clustering. The effectiveness and superiority is proved by analysis on whole life data from intelligent maintenance system (IMS).
1. Introduction
Port crane is important material handling equipment in ports whose power mainly comes from its hoisting mechanism which is used for containers moving. Rolling bearing is the key rotating supporting component in hoisting mechanism. In case of sudden failure, economic losses even casualties will be caused. It is meaningful for reducing sudden failure probability and improving reliability by analyzing rolling bearing monitoring signals and assessing performance degradation condition accurately. Port cranes usually work in heavy-duty, high-speed environment, and the signals usually show nonlinear, nonstationary, and nonperiodic character. It is a hot topic that how to get an effective degradation condition assessment technique.
Degradation condition assessment mainly comprises two key procedures, including degradation feature extraction and assessment model constructing (Han et al., 2019). The purpose of degradation feature extraction is calculating health indicators reflecting performance degradation pattern. The merits of the health indicators will determine the accuracy of degradation condition assessment. At present, some typical degradation features are studied based on time domain analysis, frequency domain analysis, and time–frequency domain analysis methods; some health indicators are proposed as degradation features such as root mean square (RMS), variance, average frequency, and so on (Araneo et al., 2017; Xiao et al., 2015). Considering that most engineering vibration signals are nonlinear and unstable, a complexity analysis method based on information entropy and fractal theory are usually introduced to provide an effective way of degradation feature extraction. Some commonly used methods include fuzzy entropy (Wang et al., 2018a, 2018b), sample entropy (Yang et al., 2018), approximate entropy (Saif et al., 2018), and fractal dimension (Wang et al., 2015). Along the train of thought, the C0 complexity analysis method which is proposed based on the Fourier transform provides a new idea for degradation feature analysis. The C0 complexity method has achieved good effects on electroencephalogram and speech signals; the advantage of this method lies in convenient calculation and low signal requirement (Cai et al., 2018; Ye et al., 2018a, 2018b). Ye et al. (2018a, 2018b) proved that spectrum entropy and C0 complexity curves are able to describe dynamic character of continuous chaotic systems accurately and effectively. Xu and Zhang (2015) proposed an endpoint detection algorithm based on the combination of improved C0 complexity and Mel-scale frequency cepstral coefficients similarity. At present, there are few studies on degradation analysis of mechanical equipment using the C0 algorithm which contains potential value in nonlinear signal analysis.
In the health condition assessment procedure, degradation condition division is a difficult step; the condition numbers and boundaries are difficult to determine, and it is apparent that subjective division methods are not scientific. In the physical degradation process, signals within the same condition are usually continuous at time dimension, so it is pivotal to divide degradation conditions as its physical feature; in the meantime, it will lay an effective foundation for the on-site online health condition assessment.
An unsupervised clustering method is an effective technique in degradation condition division and some typical clustering methods include K-means clustering (Li et al., 2019), fuzzy C-means clustering (FCM) (Wang et al., 2017), Gath–Geva (GG) fuzzy clustering (Wang et al., 2018a, 2018b), and Gustafson–Kessel (GK) clustering (Chen et al., 2018) are usually applied. The GG clustering algorithm has some improvement comparing with FCM and GK clustering methods (Abonyi et al., 2005; Yu et al., 2017). The hot spots of the GG clustering algorithm mainly focus on fault diagnosis, and the clustering number and time continuity constraints in degradation process are less considered.
In summary, in allusion to degradation feature extraction and degradation condition division, a technique based on the multidimensional degradation feature and GG fuzzy clustering algorithm is proposed. The multidimensional degradation feature includes C0 complexity, RMS, and curved time (CT) parameter which is more in line with the performance degradation process. The GG algorithm is introduced to cluster different conditions in the degradation process. The bearing lifetime vibration signal from the IMS bearing test center is used to verify the method.
The article is organized as follows: Section 2 introduces the basic theory in this article. Section 3 expounds C0 complexity and its analysis on the logistic chaotic map. The procedure of degradation condition clustering technique is expounded in Section 4. In Section 5, the technique is verified and the results are discussed. Finally, the conclusion is given in Section 6.
2. Basic theory
2.1. C0 complexity
The C0 complexity algorithm is a nonlinear analysis method which has the advantage of fast calculation speed (Guo et al., 2017). The main idea lies in describing signal complexity with the proportion of nonregular components quantitatively in the sequence.
Assuming one-dimensional time series is x(t), t = 0, 1, 2, …, M−1, the basic procedures are as follows: Making the discrete Fourier transform on x(t) Defining Calculating mean square value of Retaining the spectrum, which is larger than the mean square value Processing the Fourier inverse transform on Calculating C0 complexity as
Generally speaking, C0 complexity is introduced for describing complexity of time sequence in quantity. The value of C0 is proportional to the degree of fluctuation pattern. In this article, C0 complexity is introduced to describe the performance degradation feature of rolling bearing.
2.2. GG fuzzy clustering
Fuzzy maximum likelihood distance is introduced in GG fuzzy clustering, and the basic procedures are as follows: Supposing that the clustering sample set is X = {X1, X2, …, Xn} and the number of feature vector is n. First, a clustering number is initialized as c (c ≥ 2), and the membership matrix is set as U = [u
ip
]c×n, u
ip
∈ [0, 1], i = 1, 2, …, c, p = 1, 2, …, n, where u
ip
represents the degree of the pth sample belonging to the ith clustering. Setting terminal parameter ε(ε > 0) and initializing membership matrix U. Calculating clustering centers Calculating fuzzy maximum likelihood distance as the following formula Setting the objective function J and achieving the optimal classification by minimizing the target Updating the membership matrix
where i = 1, 2, …, c and k = 1, 2, …, n; the updating will stop when ||U(l) − U(l−1)|| < ε.
2.3. Clustering evaluation indexes
It is a key issue that how to evaluate the degradation condition clustering result. Considering maximum membership degree principle is usually adopted in this type of clustering methods; therefore, three evaluation indexes are introduced toward membership degree matrix U are proposed in this section.
Among the indexes, the classification coefficient α and average fuzzy entropy β are introduced, and the definitions are as follows
The better the clustering effect, the closer the classification coefficient to 1 and the closer the average fuzzy entropy to 0.
Furthermore, sequence dispersion index (SD) is proposed in this section to measure the continuity at time dimension during the same degradation condition. The definition is as follows.
Assuming I is the label sequence in clustering result, and n is the number of samples of this cluster. m − 1 is defined as the D-value between the maximum and minimum values of I. The sequence dispersion of this cluster is defined as the following formula
Obviously, the value of b will be zero when I is a continuous sequence. The more discontinuous the sequence I, the greater the value of b.
Assume that the entire sequence is divided into number of c classes. The SD index of the clustering result is calculated as follows
From the definition above, the clustering method will get a good effect if SD is closer to zeros. The larger the index, the lower the time aggregation and the worse the clustering effect.
3. Degradation feature analysis based on C0 complexity
A logistic chaotic map sequence is used as the simulation signal for analysis of C0 complexity. The simulation signal expression is as follows (Jabbari and Mohasefi, 2019) Logistic chaotic map bifurcation diagram.
Algorithm parameters and calculating results.

Trends for different complexity methods: (a) C0 complexity, (b) fuzzy entropy, (c) approximate entropy, and (d) sample entropy.
From the figure above, it is evident that C0 complexity and fuzzy entropy are able to reflect the increscent complexity during the interval (3.45, 3.57) as well as the downtrend of the complexity caused by the mixed periodic cycle; however, the approximate entropy and sample entropy algorithms cannot reflect this phenomenon effectively. Considering the calculation speed, C0 complexity is much faster than the other three methods because it involves only the Fourier transform and inverse transform operations. Therefore, the C0 complexity parameter has a good correlation with the sequence complexity. In addition, this method has few parameters and fast calculation speed.
4. Flowchart of degradation condition clustering
A degradation condition clustering technique based on multidimensional degradation features and GG fuzzy clustering is proposed. The flowchart is shown in Figure 3. Flowchart for degradation condition clustering.
From the flowchart above, a lifetime vibration signal of mechanical equipment is obtained firstly, Multidimensional degradation features are extracted, respectively, and an RMS indicator is calculated to characterize signal energy accumulation. C0 complexity is introduced to reflect complexity changing information. To consider the agglomeration degree of the same degradation condition on time dimension, the CT parameter is introduced to reflect time character during the performance degradation process. The value of the CT indicator is calculated by normalizing lifetime time T and mapping it as the following formula Comparison of T and CT. Note: CT: curved time.

After feature extraction, GG fuzzy clustering is performed on the degradation division so that different degradation condition is able to be clustered with the fuzzy matrix calculated by the GG fuzzy clustering method. Finally, the clustering effect is evaluated by the three indicators expounded above.
5. Instance analysis
5.1. Bearing lifetime dataset
The whole lifetime dataset used in this section is from the IMS Center at the University of Cincinnati (Antoni and Borghesani, 2019). The test rig is shown in Figure 5. The testing bearing is Rexnord ZA-2115 double-row roller bearing, and the number of rollers is 16. The roller group pitch diameter is 75.501 mm, the roller diameter is 8.4074 mm, and the contact angle is 15.17°. Map of accelerating test rig.
One lifetime dataset is introduced for instance analysis which contains 984 groups of samples, and the final fault occurs at the outer ring. A time domain waveform is shown in Figure 6 in which the sampling interval has been filtered. It seems that the amplitude begins to increase at about the 700th group sample, and a quantitative feature needs to be extracted to describe the degradation process. Time domain wave of the lifetime dataset.
5.2. Multidimensional degradation feature
Degradation feature analysis is performed on each group to calculate a multidimensional feature vector [RMS, C0, CT]. The parameter r is set as 10 for the C0 complexity method. The normalized feature vector is shown in Figure 7. It is clear that C0 complexity gradually decreases as the performance degradation deepens, whereas the trend of the RMS curve is opposite. The reason lies that the RMS indicator represents energy accumulation, whereas C0 complexity describes the complexity degree. It shows that complexity degree of the signal will gradually decrease as the performance degradation degree increases, and signal energy will increase with the degradation degree. The two curves show a phased character which reflects different conditions of bearing performance degradation. In addition, the “CT” parameter curve marked in green is more effectively related to the overall trend of performance degradation through the mapping operation. Multidimensional degradation features.
5.3. Degradation condition clustering
GG fuzzy clustering after the normalization process on a feature vector [C0, RMS, CT] is performed. The degradation condition is divided into four categories in this article, and the parameters are set as c = 4, m = 2, and ε = 0.00001. The contour diagram is shown in Figure 8. It is clear that the contour line is arbitrary which means the clustering technique has a low requirement for data source distribution and has a good adaptability. Contour map for Gath–Geva clustering.
The clustering result is shown in Figure 9. The entire degradation process is divided into four conditions marked as normal, slight, severe, and failure; each condition has preferable continuity on the time scale. It is apparent in Figure 9(a) that the bearing is kept in normal condition for a long time before the 520th group and C0 complexity is maintained at around 0.7. When degradation occurs slightly at the 520th group, C0 is very sensitive and drops rapidly. The bearing performance was seriously degraded and the value of C0 was maintained at around 0.4 with little value rebounded in the third condition basically. In the failure condition, C0 is generally low and some discrete points of numerical anomalies appear. Clustering effect for Gath–Geva clustering: (a) two-dimensional clustering effect and (b) three-dimensional clustering effect.
5.4. Influence analysis of CT feature
To analyze the influence of the CT feature clearly, The GG fuzzy clustering model and parameter are kept unchanged; two different degradation feature vector schemes are introduced for comparative analysis including two-dimensional features [C0, RMS] and three-dimensional features [C0, RMS, T], where T is the time parameter without mapping. The clustering effect diagram is shown in Figure 10, and the quantitative evaluation is shown in Table 2. Compared with the clustering effect in Figure 9, the continuity of the same condition is intermittent which is not in line with the real physical process so that the same condition will continue a consequent stage whenever the indicators fluctuate at a certain interval has preferable continuity on the time scale. Gath–Geva clustering effect for different features: (a) two-dimensional feature [C0, RMS] and (b) three-dimensional feature [C0, RMS, T]. Note: RMS: root mean square. Quantitative evaluation for different features. Note: RMS: root mean square; CT: curved time.
It is clear that three schemes have similar values on the classification coefficient (α); the proposed method have the smallest value in fuzzy entropy parameter (β) and sequence dispersion (γ), indicating that the clustering time aggregation and clustering effect are good because of the CT parameter. The reason lies that the CT parameter is able to enhance the condition relevance in time dimension and imply a gentle tendency in early stage and a dramatic trend in the later stage.
5.5. Influence analysis of clustering methods
Keeping three-dimensional feature vectors [C0, RMS, CT] unchanged, GK clustering and FCM clustering are used for comparative analysis. The clustering effect of two clustering methods is shown in Figure 11. The quantitative comparison is shown in Table 3. Clustering effect for different clustering methods: (a) GK clustering effect and (b) FCM clustering effect. Note: GK: Gustafson–Kessel; FCM: fuzzy C-means. Quantitative evaluation results for different methods. Note: GK: Gustafson–Kessel; FCM: fuzzy C-means; GG: Gath–Geva.
In the table above, the proposed GG clustering technique has the largest classification coefficient (α) and the lowest average fuzzy entropy (β), showing a good clustering effect. In addition, it is clear that sequence dispersion of the three methods is relatively low because of the “CT” parameter, indicating a good time aggregation effect. However, as to the index of average fuzzy entropy (β), the chosen two contrastive methods are higher than GG clustering which means that the membership matrix is easy to cause condition misjudgment. For example, in Figure 11(a), some sampling groups of a slight degradation condition from 500th to 600th are misjudged as a normal condition by the GK clustering algorithm. Above all, the introduced GG clustering algorithm has a good clustering effect because of fuzzy maximum likelihood distance introduced in this clustering technique.
In summary, the C0 complexity parameter has a good performance degradation indication capability. Combining with the proposed multidimensional degradation feature, the GG clustering algorithm has a good clustering effect for degradation condition division.
6. Conclusions
GG fuzzy clustering for rolling bearing degradation condition division using multidimensional degradation features is used in this article. By means of instance data analysis, the following conclusions are obtained: C0 complexity is able to reflect the irregular components proportion and describe the regularity in the process of performance degradation effectively; in addition, it is sensitive to complexity changing and the calculation speed is fast. The analysis on the logistic chaotic map sequence and the instance signal show a good effect. The proposed “CT” parameter is able to reflect the temporal character of the dataset, and it is more in line with the performance degradation regular of the mechanical device by the mapping operation of the exponential function. The GG clustering method is able to cluster data samples with arbitrary shapes; by using the multidimensional features, the clustering accuracy will be improved as well as the time aggregation degree within the same condition. The proposed SD parameter is able to reflect time aggregation effect of clustering. The research results show that the degradation process of rolling bearing performance shows obvious stage, and different degradation condition shows certain agglomeration inside. More scientific degradation feature extraction method is apt to be studied based on the multidimensional features. At the same time, an effective real-time degradation condition assessment technique is also necessary in the following studies.
Footnotes
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by National High Technology Research Development Plan (2013AA041106), National Natural Science Foundation (31300783), and China Postdoctoral Science Foundation (2014M561458).
