Next Article in Journal
Neuroevolution-Based Adaptive Antenna Array Beamforming Scheme to Improve the V2V Communication Performance at Intersections
Next Article in Special Issue
A Survey on Data-Driven Predictive Maintenance for the Railway Industry
Previous Article in Journal
A Novel Reduced Graphene Oxide Modified Carbon Paste Electrode for Potentiometric Determination of Trihexyphenidyl Hydrochloride in Pharmaceutical and Biological Matrices
Previous Article in Special Issue
Interpretable Detection of Partial Discharge in Power Lines with Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Automated Data Fusion-Based Gear Faults Classification Framework in Rotating Machines

Department of Mechanical, Aerospace and Civil Engineering, University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(9), 2957; https://doi.org/10.3390/s21092957
Submission received: 14 February 2021 / Revised: 12 April 2021 / Accepted: 20 April 2021 / Published: 23 April 2021

Abstract

:
The feasibility and usefulness of frequency domain fusion of data from multiple vibration sensors installed on typical industrial rotating machines, based on coherent composite spectrum (CCS) as well as poly-coherent composite spectrum (pCCS) techniques, have been well-iterated by earlier studies. However, all previous endeavours have been limited to rotor faults, thereby raising questions about the proficiency of the approach for classifying faults related to other critical rotating machine components such as gearboxes. Besides the restriction in scope of the founding CCS and pCCS studies on rotor-related faults, their diagnosis approach was manually implemented, which could be unrealistic when faced with routine condition monitoring of multi-component industrial rotating machines, which often entails high-frequency sampling at multiple locations. In order to alleviate these challenges, this paper introduced an automated framework that encompassed feature generation through CCS, data dimensionality reduction through principal component analysis (PCA), and faults classification using artificial neural network (ANN). The outcomes of the automated approach are a set of visualised decision maps representing individually simulated scenarios, which simplifies and illustrates the decision rules of the faults characterisation framework. Additionally, the proposed approach minimises diagnosis-related downtime by allowing asset operators to easily identify anomalies at their incipient stages without necessarily possessing vibration monitoring expertise. Building upon the encouraging results obtained from the preceding part of this approach that was limited to well-known rotor-related faults, the proposed framework was significantly extended to include experimental and open-source gear fault data. The results show that in addition to early established rotor-related faults classification, the approach described here can also effectively and automatically classify gearbox faults, thereby improving the robustness.

1. Introduction

Since the rise of machines and consequent industrial revolutions, rotating machines have become an integral and inevitable asset within virtually all industrial setups, irrespective of the sector. This is mainly due to their versatility and ability to adapt to the incredibly harsh operating environments that prevail in most industries. Components such as electric motors [1,2,3,4] and gearboxes are integral to the functioning of numerous industrial rotating systems and are often envisaged to be robust enough to withstand highly dynamic operations.
However, despite their often-resilient designs and configurations, gearboxes still experience failures, some of which have direct devastating consequences on income, environment, and human safety. Although the value of gearboxes has never been undermined, their criticality (especially due to their contributions to failure rates) to alternative energy systems, however, has further increased the intensity of the scrutiny they have received in recent times. For instance, the study by Spinato et al. [5] highlighted that wind turbine (WT) gearboxes account for the highest mean-time-to-repair among all onshore sub-assemblies. Similarly, Gray and Watson [6] also highlighted that as much as 33% of total operational downtime of energy generation systems can be attributable to their gearboxes. The findings from these studies [5,6] have also been corroborated by regional surveys and other technical reports. For instance, 3-year WT operational data from offshore wind farms in Netherlands stated that gearbox failures can account for more than 55% of total annual downtime, which sometimes corresponds to over 52% of energy not served [7]. Feng et al. [8] also conducted a similar survey for selected wind farms in the United Kingdom and reported nearly identical findings.
In general, incessant failure of the gearboxes of industrial systems are often attributed to inaccurate estimation of actual operating loads, unforeseen changes to loading conditions during operations, faulty component design, and/or inaccurate repair/replace decisions [9]. Traditionally, maintenance interventions (mainly repair and replace) have served as remedies to such failures since the amount of consideration allotted to downtime was insignificant in the past. In contrast to such times, modern-day operations are very lean and mainly customer-oriented, owing to fierce global competitions. This is perhaps the reason for the surge in the popularity of predictive and condition-based maintenance (CBM) strategies [10,11,12], whereby industrial assets dictate the frequency of maintenance interventions. Just as gearboxes have earned themselves the status of inevitability within most industrial operations, vibration monitoring (VM) [13,14,15] is arguably one of the most widely applied CBM techniques owing to the established fact that all structures (static or rotating) exhibit their own peculiar individual dynamic characteristics. The fundamental premise of VM is to adequately understand, track, and determine the trend of these characteristics for individual critical assets, so as to determine deviations at incipient stages before the occurrence of catastrophic failures. Despite the huge successes recorded with well-established VM techniques in time [16,17], frequency [18], and time–frequency [19] domains, the rigour often associated with individualised synthesis of large volumes of data acquired from each measurement location on a typical rotating machine can prolong decision-making, which may lead to fatal consequences when dealing with critical safety systems. To further compound this problem, most modern-day industrial rotating machines are multi-component (e.g., gears, bearings, drive belts, rotors, electric motors, couplings, etc.) and multi-state (e.g., speed and load variations), thereby requiring a holistic approach to VM.
In order to alleviate this limitation and better optimise conventional rotating machine VM approaches, research endeavours over the past few decades have been exploring mechanisms through which VM information can be harmonised into single but representative frameworks. Such approaches are generally referred to as data fusion, information fusion, or hybridisation. In general, data fusion involves the combination of data and information from different sources to obtain enhanced accuracy that may not be achievable from a single source. This approach has the potential to overcome some of the real-life problems that have plagued the use of VM for incipient fault detection and diagnosis. Data fusion can occur at different stages of a typical faults classification process, with sensor, feature, and decision levels being the most common stages. At the sensor or data acquisition level, data from various commensurate sensors are directly fused, after which the most representative features are used to take decisions [20,21,22,23,24]. At the feature level, individual sensors are used to acquire data separately, and the required feature vectors extracted. The feature vectors are eventually fused together and form the basis for decision-making of operational state classification [25,26,27]. Finally, the decision-level fusion approach harmonises the estimated decisions that have been separately drawn from the feature vectors of individual sensors [28,29].
Generally, machinery fault diagnosis approaches that embed machine learning techniques within them usually entail two key phases, namely feature extraction and fault identification. The feature vectors generated during the feature extraction phase are usually applied as inputs in the chosen machine learning technique(s) for the fault identification phase [30]. The fundamental essence of fault identification is to distinguish healthy from faulty machine conditions, based on the extracted features, which is also equivalent to the mapping of information from the feature space to the fault space [30]. Classifiers and statistical learning methods, such as k-nearest neighbour (k-NN) algorithms, Bayesian classifier, support vector machine (SVM), and artificial neural network (ANN) have been widely utilised in structural health monitoring systems of rotating machinery [30].
Kalman filtering [31], weighted average [32], algebraic functions [33], Bayesian estimators, nonlinear system fusion, and adaptive observers [34] are some of the earlier and well-established traditional sensor fusion approaches. Despite the usefulness of the academic research advancements achieved through these techniques, their mathematical intrusiveness could be a reason for their limited application for real-life rotating machine faults classification. For instance, Bayesian estimators are popular; the processing frameworks for dynamic models that are nonlinear often entail some multi-dimensional integrals that are often analytically intractable, thereby leading to estimate difficulties. In addition, the generated outputs are sometimes generic and associated with multimodality, asymmetries, and discontinuities [35]. With regards to Kalman filtering, linearised transformations are only proficient when error propagations can be estimated using a linear function; otherwise, the generated linearised outputs are negatively impacted and lead to complete divergence [36]. Adaptive observers can be very power and accurate when dealing with continuous time domain signals, but their effectiveness dwindles under low control and switching frequency, which is often an attribute of nonlinearities of several rotating machine components, especially gear systems [37]. Moreover, each of these classes of techniques [31,32,33,34] are often focussed on a single stage fusion as well as applied to a single machine component, thereby increasing computational complexity and, in turn, further complicating an already tedious faults classification process. On the contrary, the recently developed composite coherent spectrum (CCS) [38] and poly-coherent composite spectrum (pCCS) [38,39,40,41] significantly reduce computational stages associated with data fusion since it embeds both sensor-level and feature-level fusion into a single framework.
While previous studies on CCS [38] and pCCS [38,39,40,41] have yielded encouraging outcomes, they have only been applied to rotor-related machine faults (mainly misalignment, rub, crack, and bow), which does not adequately represent the multi-component configuration of most modern-day rotating machines. Besides the restriction of scope of the founding CCS and pCCS studies with respect to rotor-related faults, their diagnosis approach is manually implemented, which could be unrealistic when faced with routine condition monitoring of multi-component industrial rotating machines, which often entails high-frequency sampling at multiple locations. The automated framework applied here was recently presented by Yunusa-Kaltungo and Cao [42] to help address the laborious nature of manual faults classification of CCS. Although the framework encompassed initial feature generation through CCS data fusion, data dimensionality reduction via principal component analysis (PCA) and subsequent faults classification was achieved using several machine learning techniques, including ANN, SVM, k-NN, etc. The outcomes of the automated approach are usually a set of visualised decision maps representing individually simulated scenarios, which simplifies as well as minimises diagnosis-related downtime by allowing asset operators to easily identify anomalies at their incipient stages without necessarily possessing vibration monitoring expertise. Additionally, comparisons between other VM techniques indicated that CCS has significant advantages as a feature extraction method, owing to its ability to greatly reduce potential complexities that are sometimes associated with the machine learning input datasets. The results obtained from the study [42] were encouraging and showed that ANN was most compatible with CCS. However, the study was also limited by its application to only rotor-related faults.
This study adequately extends previous works on CCS and pCCS based on two main premises. Firstly, it establishes a framework by which the diagnosis of multiple classes of rotating machine faults can be automated through machine learning algorithms. The second major contribution of this study is that it significantly builds upon the encouraging results obtained from the preceding part of this approach that was limited to well-known rotor-related faults [42] by incorporating gearbox fault detection into a single framework. Hence, the extension provided here now considers an entirely different and unique class of rotating machine components—the gearbox, so as to complement earlier findings and ascertain robustness. Additionally, the proposed approach is primarily based on tools and features that are universally established across academia and industry (especially amplitude spectrum), thereby easing the transfer of theoretical knowledge into practice.
To accomplish this, the paper initially compares the proficiency of its approach to that of earlier related approaches in Section 2, after which a brief theoretical overview of the proposed framework is provided in Section 3 such that the current paper can be fully comprehended without the need to consult earlier articles. Section 4 provides full details of the experimental designs, with particular emphasis on the experimental rig configuration, types of machine operating conditions simulated, technical specifications of instruments, and signal processing parameters. Previous studies [42] have already recommended several rotor fault detection features, but it is uncertain that all of such features will adequately support the computational effectiveness of the current study. Therefore, in Section 5, the performance of relevant features is initially examined, after which the most influential features are then identified. In Section 6, the results of faults classification based on the proposed approach are presented as well as explanations of the implications of the findings for VM of rotating machines (in this case, gearbox faults). Section 7 provides the validation of the applied method with independent public datasets. Finally, Section 8 concludes the study and highlights possible future directions.

2. Comparison with Closely Related Works

The study of fault diagnosis in rotating machines is well-established and continues to generate spates of useful but sometimes closely related outputs in some cases, which makes it imperative to compare and contrast to identify niche areas. Therefore, this section is based on a comparative analysis between closely related approaches in fault diagnosis of rotating machines and the current study, so as to clearly highlight areas of potential interface, overlap, variation, limitation, and superiority. In order to better show the advantages of the proposed method and point out future research directions, Table 1 provides comparisons with other recent studies in a similar area. The main criteria used for comparison are the data types, classification algorithm(s), application of data fusion, and fault classes considered.
Yunusa-Kaltungo et al. developed CCS [38] and pCCS [38,39,40,41], which significantly rationalised computational stages associated with fault diagnosis through data fusion by embedding both sensor- and feature-level fusion into a single framework. However, their application has been limited to rotor-related machine faults and entail manual classification, which will increase the downtime related to fault diagnosis. This limitation led to the proposal of an automated framework [42] that still used features generated via CCS but further involved data dimensionality reduction by PCA and eventual machine learning-based faults classification. The outcomes were very encouraging especially that the study [42] exposed the compatibility of ANN with CCS but study was again confined to rotor-related faults thereby not all encompassing.
Cao et al. [43] developed a deep transfer learning approach based on a convolutional neural network (CNN) algorithm, and their study advocated the suitability of the approach for deep feature extraction and gear fault diagnosis. Similarly, Shao et al. [44] developed a CNN-based deep transfer learning framework for mechanical fault diagnosis and classification, while Soualhi et al. [45] proposed a health indicator fed into an adaptive neuro-fuzzy inference system (ANFIS) to detect the state of health of a typical system and then diagnose sources of anomalies. The data collected through this method are electrical signals, mainly current signals as opposed to mechanical signals such as vibrations for non-invasive benefits, since the current, voltage, or power sensors that are already integrated into the control systems of electrical machines can be used. Azamfar et al. [46] developed a novel multi-sensor data fusion methodology based on 2-D CNN for gearboxes fault diagnosis using motor current signature analysis. Zhang et al. [47] proposed a novel unsupervised learning algorithm named fast intrinsic component filtering (FICF) for the fault diagnosis of rotating machinery. These studies have no doubt enhanced the knowledge around the fault diagnosis of rotating machines. However, they are limited by either focus on the class of singular faults (e.g., rotor faults or gear faults or bearing faults alone) which implies that alternative approaches will need to be considered for other fault classes, thereby increasing rigour and downtime or computational intensiveness of CNN-based approaches. Table 1 provides more targeted merits, demerits, and coverage of individual study classes.

3. Theoretical Overview of the Approach

3.1. Mathematical Representation

It is vital to highlight that full details of the automated faults classification framework applied here have been provided in the preceding article that focussed on rotor-related faults [42]. However, the provision of high-level description, here and again, was adjudged useful, so as to allow this second part to be comprehensible as a standalone section. The process commences with the use of the CCS approach [38] to fuse amplitude spectra computed from the time domain datasets acquired by individual VM sensors installed on the rotating machine. The mathematical representation of the CCS process is iterated in Equations (1)–(3) [38,39,40,41]. The fundamental rationale behind the CCS is that it eliminates downtime associated with routine VM processes, especially when such monitoring returns a no-fault result from the studied machine. During such VM processes, technicians are required to analyse data from all measurement locations on the machine, but this process is minimised by the CCS as only one spectrum needs to be routinely observed. The only instance that would warrant analysis of individual spectra is when a deviation from the single CCS is observed. This stage of the fusion is referred to as Stage 1 in Figure 1 and further explained thus:
If the number of measurement points on a particular rotating machine is b , each of which is furnished with a VM sensor, then the vibration signals acquired from individual sensor can be divided into n s equal-length segments. The coherent cross-power spectral density of the signals from the p th and ( p + 1 ) th measurement points at a frequency f h can be defined as:
S x p γ p ( p + 1 ) 2 x p + 1 r ( f h ) = [ X p r ( f h ) γ p ( p + 1 ) 2 ( f h ) X p + 1 r * ( f h ) ]
where X p r ( f h ) is the discrete Fourier transform (FT) of the r th segment of the signal x p , and X p r * ( f h ) is its complex conjugate, for p = 1 ,   2 , , b 1 . γ p ( p + 1 ) 2 ( f h ) is the coherence of the signals x p and x p + 1 for background noise suppression.
Hence, each of the r t h segments from each signal can be fused into a single component, X C C S r ( f h ) , thus:
X C C S r ( f h ) = ( S x 1 γ 12 2 x 2 r ( f h ) S x 2 γ 23 2 x 3 r ( f h ) S x ( b 1 ) γ ( b 1 ) b 2 x b r ( f h ) ) 1 ( b 1 )
The CCS for the entire machine can then be calculated as:
S C C S ( f h ) = r = 1 n s X C C S r ( f h ) X C C S r * ( f h ) n s
The S C C S ( f h ) is a sequence of complex numbers that enables the estimation of the amplitude spectrum of the CCS according to Equation (4):
A C C S ( f h ) = | 2 N S C C S ( f h ) |                     ( h = 1 ,   2 , , N / 2 )
In addition to the earlier CCS harmonic amplitudes that only offer a single-point value of differentiation that could be similar for several harmonics, here, we also consider the spectral energy that can be estimated according to Equation (5).
For a typical A C C S ( f h ) computed as per Equation (4) at a frequency f h , where f h = ( h 1 ) d f , h = 1 , 2 , , N / 2 , N is the number of data points and d f is the frequency resolution, the SE between the selected harmonics at intervals of d f can be defined as:
A S E ( f h ) = i = h 10 h + 10 A C C S ( f i ) × d f   ( h = 11 , 12 , , N / 2 10 )
In order to analyse the gear fault, in addition to considering the harmonics of rotating speeds of gears and shafts, the harmonics of the gear mesh frequency (GMF) should also be included in CCS calculation because sidebands around the GMF and its harmonics contain information on gearbox faults [48]. The GMF can be calculated by:
f G M F = z t × n R
where z t is the number of teeth on pinion, and n R is the rotating speed of the pinion.
Once the CCS harmonic amplitudes and their corresponding SEs of interest (depending on the fault types considered, e.g., low frequency for rotor-related and higher frequency for gear faults) have been obtained, Stage 2 of fusion involves their standardisation, dimensionality reduction, and harmonisation based on PCA [49,50,51] and ANN [52,53,54]. The computational steps required for Stage 2 are described by Equations (7)–(11).
Owing to the variations in the amplitude ranges that may be associated with the diagnosis of a complex multi-component system such as that considered here (e.g., shaft and gear mesh frequency harmonic amplitudes), dimensionality reduction through PCA would require some prior standardisation of the input data A . In the matrix A R m × n , m is the number of samples and n is the number of features (dimensions), a i j represents typical elements of the matrix, while x i j is a corresponding element of the standardised matrix X .
The element x i j of the standardised matrix X is defined as:
x i j = ( a i j A ¯ j ) S j
where A ¯ j is the sample mean of the elements of the j th column of matrix A , and S j is the sample standard deviation of the j th column of A , which is mathematically represented as:
S j = i = 1 m ( a i j A ¯ j ) 2 m 1
The computation of PCs of X reduces to the solution of an eigenvalue–eigenvector problem:
C X V = V Λ
where C X is the covariance matrix of X , and V is the orthogonal matrix whose j th column is the j th eigenvector of C X , corresponding to the j th largest eigenvalue of C X which is the j th diagonal element of the diagonal matrix Λ .
The columns of the matrix V R n × n are orthogonal unit vectors and are referred to as the right singular vectors of X .
The calculation of the score matrix (result) T R m × n for a PCA can be mathematically represented as:
T = X V
After dimensionality reduction and selection of the PCs combination that offers the highest representation, ANN is then used to classify the different experimentally simulated machine conditions as per Equation (11) [42]:
y = f ( W T x ) = f ( i = 1 N W i x i + b )
where f is the activation function, W are the weights and b is the scalar bias term.

3.2. Operational Description of The Approach

In this section, we provide a step-by-step description of the individual operational stages of the applied approach, so as to foster better understanding.

3.2.1. Training Steps

  • Feature extraction
    Obtain the CCS harmonic amplitudes and/or GMF harmonic amplitudes as well as their corresponding SEs for numbers of segment averages from a dataset with known health conditions. This provides the input data matrix A R m × n , where m is the number of samples (averages) and n is the number of features (dimensions).
  • PCA application
    Compute A ¯ j , S j , T , V as stipulated by Equations (7)–(10).
  • Dimensionality reduction
    Here, we consider the L largest singular values to obtain the truncated score matrix T L R m × L , where PC1 is the first column of T L , PC2 is the second column of T L , and so on.
  • Training
    Use T L as the input of the ANN model. Train the model with labels (individual health conditions) so as to obtain the classifier as depicted by Equation (12):
    z = f ( x 1 ,   x 2 , ,   x n )
    In Equation (12), PC1 is equivalent to x 1 , PC2 is equivalent to x 2 , and PCn is equivalent to x n . Additionally, z denotes the resultant class of machine health condition.

3.2.2. Automatic Classification Steps

  • Feature extraction
    Obtain the CCS harmonic amplitudes and/or GMF harmonic amplitudes and their corresponding SEs for a new dataset that does not possess any labels so as to obtain the input data matrix B R 1 × n , where n is the number of features (dimensions) and b i j represents a typical element of the matrix B .
  • Linear transform
    Project the data matrix B into the same linear space as the PCA obtained from the training steps. The element y i j of the transformed matrix Y R 1 × n is defined as:
    y i j = ( b i j A ¯ j ) S j
    Similarly, the transformed score matrix S R 1 × n is defined as:
    S = Y V
  • Dimensionality reduction
    By considering only the L largest singular values, we obtain the transformed truncated score matrix S L R 1 × L , where PC1 is s 11 , PC2 is s 12 , and so on.
  • Classification
    With S L as the input to the trained classifier z = f ( x 1 ,   x 2 , ,   x 3 ) , we can obtain the health condition classification result z = f ( s 11 ,   s 12 , ,   s 1 n ) .
After training, by computing just 3 equations, we can determine the health condition from the CCS of a new vibration dataset. This operation will be automatically performed for all steps.

4. Experimental Design and Data Acquisition

Various operation conditions (mainly faults and speed variability) usually associated with typical industrial rotating machines were experimentally simulated on a laboratory scale rig, after which vibration datasets were acquired. This section offers full details of the experiments used to generate the VM datasets used in this study.

4.1. Rig Characteristics

The rig used for this experiment is a multi-component rotating machine with two main rotors coupled together by two helical-geared gearboxes. Rotational force to the entire rig is provided by a 2HP electric motor that runs at a maximum speed of 3600 RPM. The electric motor shaft is coupled to the drive end (DE) gearbox through a belt pulley system. The driven pulley is directly connected to the drive shaft of the intermediate gearbox through a stepped shaft, which then transmits motion to its driven shaft via a pair of helical gears. Finally, the driven shaft of the intermediate gearbox then serves as the drive shaft for the DE gearbox through another set of helical gears. The main structure of the rig is supported by five bearings (2 bearings for each shaft and an additional bearing for the driven pulley). The rig is fitted with a lubricating system that comprised of a pump, filter, radiator, and sump for oil circulation and cooling of DE and intermediate gearboxes. The rig rotation is regulated through a variable frequency drive. All of the rotating components of the experimental rig are covered by mesh to prevent injury due to entrapment. Figure 2 and Table 2 respectively provide an image of the experimental rig and the technical specifications of its main components.

4.2. Instrumentation

VM data were acquired through the aid of three accelerometers (one installed on bearings near each of the three gearboxes). The raw signals from the accelerometers pass through a signal condition that also powers the accelerometers, then to an analogue-to-digital converter (ADC). Table 3 provides a summary of the technical specifications of the main instrumentation.

4.3. Seeded Operating Scenarios

Considering that this was an existing multi-component rig without full knowledge of its state of health, the initial case that contained no seeded fault was termed the baseline case (BC). However, the amplitude spectra generated from the data obtained under BC displayed significantly high amplitudes at several harmonics of the machine speeds, which was adjudged to be due to inherent misalignment and unbalance faults. Therefore, BC can be classified as exhibiting rotor-related anomalies. The other two cases are the single fault (SF) and multiple fault (MF) cases. As shown in Figure 3, the SF case was simulated by introducing a slight notch on a gear tooth within the non-drive end (NDE) gearbox while the MF case was a combination of SF and additional wear on a gear tooth within the intermediate gearbox. Under all cases, VM data were collected at three distinct machine speeds (i.e., 7, 14, and 21 Hz), thereby yielding a total of nine experimentally simulated operating scenarios. In this study, a scenario represents one combination of case and speed (e.g., VM datasets for BC @ 7 Hz). The experimental flow is shown in Figure 4.

4.4. Data Acquisition and Signal Processing Parameters

For each of the scenarios described in Section 4.3, 3 VM datasets were acquired for approximately 120 s (totally 27 VM datasets in 9 scenarios). Here, two additional datasets are collected in each scenario to confirm that there are no anomalies in the experimental data. During spectrum and CCS calculation, the signal processing parameters used were 10,000 Hz sampling frequency ( f s ), 80% segment overlap, 0.5 Hz frequency resolution ( d f ), 448 number of segment averages, 20,000 as the number of FT data points ( N ), and Hanning window.

5. Feature Selection and Optimisation

A typical VM process of rotating machines is usually associated with the generation of various features, especially when dealing with those characterised by multiple components. The fundamental objective of CCS data fusion approach is to rationalise data such that the VM of rotating machines can be simplified. Since the faults considered here are rotor- and gear-related, three sets of features (1st–5th harmonics of shaft speed and 1st–5th harmonics of gear mesh frequency) were extracted after computing the CCS as per Equations (1)–(4). Owing to the high energy contents of typical GMFs, the SE of the resultant CCS was also computed as per Equation (5) so as to observe its performance as a feature. Figure 5 shows the amplitude distributions of shaft (α1–α5) and GMF (β1–β5) features across all scenarios, where it can be seen that β1 was the most consistently dominant feature at all speeds and the patterns of other features were inconsistent across different scenarios. It can also be observed that BC, at all speeds, contained a prominent β1 feature, which is unsurprising due to the existence of inherent gears/shaft misalignment and shaft unbalance. Although the selected shaft harmonic features (α1–α5) in Figure 5 were observable for all scenarios, their distribution is similarly inconsistent. With regards to harmonic distribution, the SE-based GMFs (denoted by γ1–γ5 in Figure 6) exhibited very similar trends, except for slightly higher amplitudes of higher harmonics at 21 Hz (Figure 6). This therefore implies that while all selected features are appropriate for identifying the presence of anomalies, reliance on such features alone for fault characterisation and separation is nearly impossible due to the identified inconsistencies for different scenarios. The preliminary performance comparisons performed here aided and formed the basis for the selection of the most influential features, which then formed the basis for the next stage of the faults classification framework.

6. Classification Results and Their Implications

Having established the most influential features for both shaft/rotor and gear faults in Section 5, the next stage of the analysis involves reducing the dimensionality of such features using PCA as well as examining the abilities of different combination of features to retain the highest variability. The implementation of PCA was based on the theories described in Equations (7)–(10). The four classes of features considered for this study are rotor/shaft only (α1–α5); GMFs only (β1–β5), combined rotor/shaft and GMFs (α1–α5) + (β1–β5) and SE-GMFs (γ1–γ5) features. The content distributions for 10 PCs were compared for all classes of features at all speeds as shown in Table 4. Since it is well established that the most significant information will usually reside within the first few PCs, the performance of combined PC1–2 and PC1–3 was compared, where it can be seen that PC1–3 held slightly superior information, which implies that it holds the potential to offer the most distinctive classification for all scenarios.
However, owing to the higher data requirements for the PC1–3 combination and correspondingly higher computational burden, the performance of the PC1–2 combination was additionally explored for comparative purposes as shown in Figure 7. As anticipated, the PC1–3 combination offered the best separation between the clusters that represent all machine conditions (Figure 7a–c), but the performances of several PC1–2 combinations were also encouraging, especially those that involved (α1–α5) + (β1–β5) and SE-GMFs (γ1–γ5) features in Figure 7e,f and Figure 7m–o, respectively. However, despite the good intercluster separations achieved with PCA, its manual approach makes it unsustainable for routine diagnosis of rotating machines, whereby huge amounts of data related to highly dynamic scenarios is involved. Based on this perceived limitation, there is a need for applying approaches that possess self-learning capabilities with minimal human intervention. One of such approaches is ANN, whose proficiency with the current framework has already been established with several rotor-related faults at various machine speeds.
The current study aims to consolidate as well as extend the robustness of the approach by investigating an entirely novel class of faults with regards to a CCS-based data fusion approach. The classification problem is defined as classifying the data into 3 classes (BC, SF, and MF) based on the selected features. To achieve this, 3 ANN architectures were examined for PC1–3 and PC1–2 combinations for all cases at all speeds. For the PC1–3 combination, the ANN architectures had 3–10–3, 3–20–3, and 3–30–3 configurations for ANN1, ANN2, and ANN3 respectively. For PC1-2 combinations, however, 2–10–3, 2–20–3, and 2–30–3 configurations were respectively applied for ANN1, ANN2, and ANN3. In order to ascertain the performance without PCA, ANN4 was computed without PCA and its outcome was also used for comparison (i.e., 10–30–3 for (α1–α5) + (β1–β5) and 5–30–3 for (γ1–γ5)). It is vital to note that 3–10–3, 3–20–3, 3–30–3, 2–10–3, 2–20–3, 2–30–3, 10–30–3, and 5–30–3 for individual ANN configurations, respectively representing the inputs, number of neurons for hidden layers, and outputs. The analysis was conducted based on a 70–15–15 random split of features extracted from measured VM data for training, validation, and testing, respectively. The PCA step described in Section 3.2.1 was then applied to 85% of the datasets (i.e., combined training and validation datasets), after which 15% of the datasets were then extracted from the 85% and used for validation. Subsequently, the classification steps described in Section 3.2.2 were then applied to the testing datasets. The transfer function adopted here is the sigmoid symmetric transfer function. Since the ANN type is backward propagation, scaled conjugate gradient (SCG) was used as a learning algorithm as well as for overfitting avoidance. Table 5, Table 6 and Table 7 provide full details of the configurations and performance at all speeds.
There are 2 aspects of evaluating the performance of ANNs: one is the accuracy of fitting and the other is whether overfitting occurs. As shown in Table 5, Table 6 and Table 7, the results of different ANN architectures are very similar for same scenarios (i.e., same speeds and same sets of features). For instance, at 21Hz, the accuracy of ANN with inputs of PC1–3 for (α1–α5) + (β1–β5) was significantly better than that of PC1–2 for (α1–α5) + (β1–β5) and PC1–2 for (γ1–γ5). However, PC1–2 for (γ1–γ5) has the best classification results at the other 2 speeds. The ANN computed based on inputs without PCA yielded similar results overall, except that it performed better than PC1–2 for both (α1–α5) + (β1–β5) and (γ1–γ5) at 21Hz. This was because the percentages of explained variance by PC1–2 at 21Hz were relatively small (i.e., 45.577 or 3% and 73.442 or 8%). In general, there was no significant difference in the accuracies of the ANNs trained based on these 3 features as inputs at the same speeds. Further evidence on the reason for not using (α1–α5) + (β1–β5) + (γ1–γ5) as a feature in this study are depicted in Table A1 within Appendix A. In order to demonstrate the rationale behind using ANN as the machine learning classifier in this study, the classification accuracy of ANN was compared to those obtained from three other machine learning classifiers, namely, k-NN (k = 10), naïve Bayes, and linear SVM as shown in Table 8. The comparisons were based on two input feature types. Figure 8 shows that k = 10 for k-NN had overall best results in a range of k from 1 to 15 for all considered scenarios. Therefore, k = 10 has been chosen for comparisons in Table 8. The results indicate that the ANN method outperformed all other classifiers for every scenario considered in this study.
In order to ensure good classification effects, overfitting must be avoided. Since the decision boundary of the classifier trained by the input sets with 3 dimensions or above is reasonably hyperplane in nature, it is difficult to visualise the decision rules in a 2-dimensional map. Thus, the difficulty of direct observation on whether there is an overfitting problem in an ANN with high-dimensional inputs could yield challenges in practice. On the contrary, the decision rules of ANNs trained by 2-dimensional input sets can be easily displayed. Based on this premise, it is fair to assume that 2-dimensional training input sets with PCA are advantageous when the variations in accuracy are minimal.
It is well known that overfitting is an immense threat to the abilities of machine learning algorithms to accurately detect and classify new data, owing to the incorporation of extrinsic details during the training process. In this study, it is envisaged that the application of SCG as a training algorithm will help mitigate potential problems. During individual trainings, the initial values of neurons will be reset randomly, with a corresponding random redivision of the data into 3 distinct groups for training, validation, and testing. This approach implies that training multiple times with a single input set will produce different results with slightly different decision boundaries. Figure 9 shows the decision rules of ANNs trained by PC1–2 for (α1–α5) + (β1–β5) and (γ1–γ5) at different speeds (i.e., typical results after a single round of training). The input datasets here correspond to (d–f) and (m–o) in Figure 7. The number of neurons of the hidden layer is considered as a variable for controlling potential overfitting problems. The emergence of complex boundary curves and narrow or slender envelope area within decision regions are likely indications of overfitting. For instance, the curvature of the decision boundary that exists between SF and MF regions is quite steep in Figure 9c as well as the visible elongated sharp strip area at the lower end of MF region in Figure 9i indicate that ANN3 could be associated with overfitting problems. With reference to Figure 9, decision maps generated from more neurons tend to be associated with overfitting problems. Therefore, since there are only 2 input values and 3 output classes in this classification, 10 neurons are adjudged sufficient for optimised, reasonably comprehensive, and complete classification of the cases considered in this study (i.e., increasing the number of neurons may not lead to better results). However, as more and more fault types emerge, it may be necessary to increase the number of neurons to boost the accuracy of classification. In general, the results presented here show that the initially proposed automated fault diagnosis framework is capable of identifying and classifying common gearbox faults using very simple and well-known features such as amplitude of rotor-related and GMF harmonics. This thereby provides good encouragement that the approach may be suitable for integrating rotor and gearbox fault diagnosis into a single framework in the near future.

7. Validation Dataset

In order to further examine the effectiveness of the applied method for classifying independent datasets, the study obtained publicly available gearbox fault datasets provided in an earlier study by Shao et al. [44] for validation. According to Shao et al. [44], the validation gearbox datasets were acquired from a drivetrain dynamic simulator, whereby two kinds of working conditions (i.e., rotating speed and load) were experimentally simulated. The rotating speed and load configurations were set to 20 Hz–0 V and 30 Hz–2 V. Vibration data were collected using 6 accelerometers mounted at 2 measuring positions. Position one (P1) datasets were acquired from the planetary gearbox measurement location in three directions (i.e., x, y, and z). Similarly, Position two (P2) datasets were acquired in three directions (i.e., x, y, and z) but from a parallel gearbox.
The different types of faults for both gearboxes are shown in Table 9. The datasets contain five different working conditions (i.e., four fault types and one healthy). Hence, the fault diagnosis here is based on a 5-class classification task. For each of the scenarios, 10 VM datasets were acquired for approximately 200 s. During spectrum and CCS calculation, the signal processing parameters used are 5120 Hz sampling frequency ( f s ), 80% segment overlap, 0.5 Hz frequency resolution ( d f ), 249 number of segment averages, 10,240 number of FT data points ( N ), and Hanning window. For CCS computation, two forms of data fusion approaches were considered. The former on the one hand was implemented to fuse the data from all six accelerometers mounted at the two measurement locations into a single spectrum (i.e., P1xyz+P2xyz). The latter on the other hand was implemented to fuse the data from the two accelerometers that had the same orientation (i.e., P1x+P2x, P1y+P2y, and P1z+P2z). ANN1 (2–10–5) was used as a classifier, and the PC1–2 of shaft harmonic features (α1–α5) was used as input. The analysis was also conducted based on a 70–15–15 random split of data for training, validation, and testing, respectively. Based on the linear space generated by the application of PCA to the training and validation datasets, linear transform was then implemented on the testing datasets.
The classification problem is defined as classifying the data into 5 classes (Health, Chipped, Miss, Root, and Surface) based on the selected features. Table 10 and Table 11 and Figure 10 show the results of the validation, where it can be observed that the applied approach effectively classifies all the considered validation datasets, thereby confirming the robustness. It was also observed that the outcomes obtained by integrating all six accelerometers are better than when only two accelerometers were used.

8. Concluding Remarks and Future Possibilities

Industrial rotating machines are multi-component assets which imply that a truly holistic faults classification framework should be capable of detecting anomalies associated with each component, since faults rarely occur in isolation. Previous studies on CCS data fusion have effectively rationalised vibration-based condition monitoring data of rotating machines as well as characterised most of the faults commonly encountered in practice. However, as valuable as the findings from those studies were, their applications have been restricted to rotor-related faults such as misalignment, rub, crack, looseness, and bend, which often raises questions about the efficacy of the technique. Additionally, the founding works on CCS were based on manual classifications, which may be unrealistic for routine VM that often involves the analysis of huge amounts of data on a continuous basis. In order to alleviate these challenges, this study enhances current knowledge through the following two main premises:
  • It establishes a framework by which the diagnosis of multiple classes of rotating machine faults can be automated through machine learning algorithms.
  • It incorporates gearbox fault detection into a single framework. Hence, the extension provided here now considers an entirely different and unique class of rotating machine components—the gearbox, so as to complement earlier findings and ascertain robustness.
The results observed further affirmed the proficiency of the framework for both rotor and gearbox faults. Additionally, the proposed approach is primarily based on tools and features that are universally established across academia and industry (especially amplitude spectrum), thereby easing the transfer of theoretical knowledge into practice. Considering that all studies related to the application of CCS and pCCS data fusion approaches to fault diagnosis feature generation have been purely experimental, future endeavours will focus on validating such experimental scenarios within a theoretical environment.

Author Contributions

Conceptualization, A.Y.-K. and R.C.; methodology, A.Y.-K. and R.C.; software, R.C.; validation, A.Y.-K. and R.C.; formal analysis, R.C.; investigation, R.C.; resources, A.Y.-K. and R.C.; data curation, A.Y.-K. and R.C.; writing—original draft preparation, R.C.; writing—review and editing, A.Y.-K.; visualization, A.Y.-K. and R.C.; supervision, A.Y.-K.; project administration, R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

A Matrix of the input data
A C C S Amplitude of composite coherent spectrum
A S E Spectral energy
a i j Element of the input data matrix A
A ¯ j Sample mean of elements of the j th column of A
B Matrix of the new input data without labels
b Total number of measurement points
b Scalar bias term of artificial neural network
b i j Element of the new input data matrix B
C X Covariance matrix of X
d f Frequency resolution
f G M F Gear mesh frequency
f k Frequency
f s Sampling frequency
h Ordinal number of data points
i Ordinal number of rows of matrix
j Ordinal number of columns of matrix
k Number of nearest neighbours
L Number of largest singular values
m Number of samples
N Number of data points
n Number of features
n R Rotating speed of the pinion
n s Number of vibration signal segments
p Ordinal number of measurement points
R Real matrix
r Ordinal number of vibration signal segments
S Transformed score matrix
S Coherent cross-power spectral density
S C C S Composite coherent spectrum
S j Sample standard deviation of the j th column of A
S L Transformed truncated score matrix
T Score matrix for principal components analysis
T L Truncated score matrix
V Orthogonal matrix
W Weights of artificial neural network
X Standardised matrix
X Discrete Fourier transform
X C C S Component of composite coherent spectrum
x i j Element of the standardised matrix X
Y Transformed matrix
y i j Element of the transformed matrix Y
z Resultant class of machine health condition
z t Number of teeth on pinion
α 1 α 5 Amplitudes of 1st–5th harmonics of shaft speed
β 1 β 5 Amplitudes of 1st–5th harmonics of gear mesh frequency
γ 1 γ 5 Spectrum energies of 1st–5th harmonics of gear mesh frequency
γ 2 Coherence
Λ Diagonal matrix

Appendix A

As shown in the Table A1, the performance of (α1–α5) + (β1–β5) + (γ1–γ5) is very similar to that obtained from using only (γ1–γ5) at 7 and 14 Hz, but performed even more poorly at 21 Hz. This was adjudged to be owing to the extraction of both (β1–β5) and (γ1–γ5) from GMF information. Therefore, combining (γ1–γ5) with (α1–α5) + (β1–β5) would lead to information overload and possible redundancy, which would be counterintuitive to the primary aims of the paper.
Table A1. ANN3 properties for different features using PC1–2.
Table A1. ANN3 properties for different features using PC1–2.
Feature1–α5) + (β1–β5)1–γ5)1–α5) + (β1–β5) + (γ1–γ5)
ParametersANN3ANN3ANN3
Network Structure2–30–32–30–32–30–3
Rotation Frequency7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz
Accuracy (%)Training98.999.195.499.510089.399.610087.3
Validation99.099.095.599.510089.599.510087.0
Testing99.599.093.599.510087.599.510084.0
Overall99.099.195.199.510089.099.610086.8

References

  1. Sun, X.; Shi, Z.; Lei, G.; Guo, Y.; Zhu, J. Multi-objective design optimization of an IPMSM based on multilevel strategy. IEEE Trans. Ind. Electron. 2020, 68, 139–148. [Google Scholar] [CrossRef]
  2. Sun, X.; Shi, Z.; Cai, Y.; Lei, G.; Guo, Y.; Zhu, J. Driving-Cycle-Oriented Design Optimization of a Permanent Magnet Hub Motor Drive System for a Four-Wheel-Drive Electric Vehicle. IEEE Trans. Transp. Electrif. 2020, 6, 1115–1125. [Google Scholar] [CrossRef]
  3. Sun, X.; Jin, Z.; Cai, Y.; Yang, Z.; Chen, L. Grey Wolf optimization algorithm based state feedback control for a bearingless permanent magnet synchronous machine. IEEE Trans. Power Electron. 2020, 35, 13631–13640. [Google Scholar] [CrossRef]
  4. Sun, X.; Wu, M.; Lei, G.; Guo, Y.; Zhu, J. An improved model predictive current control for PMSM drives based on current track circle. IEEE Trans. Ind. Electron. 2020, 68, 3782–3793. [Google Scholar] [CrossRef]
  5. Spinato, F.; Tavner, P.J.; van Bussel, G.J.; Koutoulakos, E. Reliability of wind turbine subassemblies. IET Renew. Power Gener. 2009, 3, 387–401. [Google Scholar] [CrossRef] [Green Version]
  6. Gray, C.S.; Watson, S.J. Physics of failure approach to wind turbine condition based maintenance. Wind Energy 2010, 13, 395–405. [Google Scholar] [CrossRef]
  7. Feng, Y.; Qiu, Y.; Crabtree, C.J.; Long, H.; Tavner, P.J. Monitoring wind turbine gearboxes. Wind Energy 2013, 16, 728–740. [Google Scholar] [CrossRef]
  8. Feng, Y.; Tavner, P.; Long, H. Early experiences with UK round 1 offshore wind farms. Proc. Inst. Civ. Eng. Energy 2010, 163, 167–181. [Google Scholar] [CrossRef] [Green Version]
  9. Musial, W.; Butterfield, S.; McNiff, B. Improving Wind Turbine Gearbox Reliability; National Renewable Energy Lab (NREL): Golden, CO, USA, 2007.
  10. Yunusa-Kaltungo, A.; Sinha, J.K.; Elbhbah, K. HOS analysis of measured vibration data on rotating machines with different simulated faults. In Advances in Condition Monitoring of Machinery in Non-Stationary Operations; Springer: Berlin/Heidelberg, Germany, 2014; pp. 81–89. [Google Scholar]
  11. Nembhard, A.D.; Sinha, J.K.; Yunusa-Kaltungo, A. Experimental observations in the shaft orbits of relatively flexible machines with different rotor related faults. Measurement 2015, 75, 320–337. [Google Scholar] [CrossRef]
  12. Yunusa-Kaltungo, A.; Kermani, M.M.; Labib, A. Investigation of critical failures using root cause analysis methods: Case study of ASH cement PLC. Eng. Fail. Anal. 2017, 73, 25–45. [Google Scholar] [CrossRef]
  13. Wang, J.; Xu, M.; Zhang, C.; Huang, B.; Gu, F. Online Bearing Clearance Monitoring Based on an Accurate Vibration Analysis. Energies 2020, 13, 389. [Google Scholar] [CrossRef] [Green Version]
  14. Li, X.; Elasha, F.; Shanbr, S.; Mba, D. Remaining useful life prediction of rolling element bearings using supervised machine learning. Energies 2019, 12, 2705. [Google Scholar] [CrossRef] [Green Version]
  15. Salomon, C.P.; Ferreira, C.; Sant’Ana, W.C.; Lambert-Torres, G.; da Silva, L.E.B.; Bonaldi, E.L.; Oliveira, L.E.D.L.D.; Torres, B.S. A Study of Fault Diagnosis Based on Electrical Signature Analysis for Synchronous Generators Predictive Maintenance in Bulk Electric Systems. Energies 2019, 12, 1506. [Google Scholar] [CrossRef] [Green Version]
  16. Ciabattoni, L.; Ferracuti, F.; Freddi, A.; Monteriu, A. Statistical spectral analysis for fault diagnosis of rotating machines. IEEE Trans. Ind. Electron. 2017, 65, 4301–4310. [Google Scholar] [CrossRef]
  17. Li, G.; McDonald, G.L.; Zhao, Q. Sinusoidal synthesis based adaptive tracking for rotating machinery fault detection. Mech. Syst. Signal Process. 2017, 83, 356–370. [Google Scholar] [CrossRef]
  18. Vishwakarma, M.; Purohit, R.; Harshlata, V.; Rajput, P. Vibration analysis & condition monitoring for rotating machines: A review. Mater. Today Proc. 2017, 4, 2659–2664. [Google Scholar]
  19. Wang, T.; Lu, G.; Yan, P. A Novel Statistical Time-Frequency Analysis for Rotating Machine Condition Monitoring. IEEE Trans. Ind. Electron. 2019, 67, 531–541. [Google Scholar] [CrossRef]
  20. Chen, Y.; Du, R.; Qu, L. Fault features of large rotating machinery and diagnosis using sensor fusion. J. Sound Vib. 1995, 188, 227–242. [Google Scholar] [CrossRef]
  21. Hansen, R.J.; Hall, D.L.; Kurtz, S.K. A new approach to the challenge of machinery prognostics. In Proceedings of the ASME 1994 International Gas Turbine and Aeroengine Congress and Exposition, The Hague, The Netherlands, 13–16 June 1994. [Google Scholar]
  22. Fang, X.; Yao, Y. In-process evaluation of the overall machining performance in finish-turning via single data source. J. Manuf. Sci. Eng. 1997, 119, 444–447. [Google Scholar] [CrossRef] [Green Version]
  23. Khan, A.; Ceglarek, D.; Shi, J.; Ni, J.; Woo, T. Sensor optimization for fault diagnosis in single fixture systems: A methodology. J. Manuf. Sci. Eng. 1999, 121, 109–117. [Google Scholar] [CrossRef]
  24. Dalpiaz, G.; Rivola, A. Condition monitoring and diagnostics in automatic machines: Comparison of vibration analysis techniques. Mech. Syst. Signal Process. 1997, 11, 53–73. [Google Scholar] [CrossRef] [Green Version]
  25. Azouzi, R.; Guillot, M. On-line prediction of surface finish and dimensional deviation in turning using neural network based sensor fusion. Int. J. Mach. Tools Manuf. 1997, 37, 1201–1217. [Google Scholar] [CrossRef]
  26. Leem, C.S.; Dornfeld, D.; Dreyfus, S. A customized neural network for sensor fusion in on-line monitoring of cutting tool wear. J. Eng. Ind. 1995, 117, 152–159. [Google Scholar] [CrossRef]
  27. Mayes, I. Use of neutral networks for on-line vibration monitoring. Proc. Inst. Mech. Eng. Part A J. Power Energy 1994, 208, 267–274. [Google Scholar] [CrossRef]
  28. Chen, Y.; Sha, J.; Wu, S. Diagnosis of the tapping process by information measure and probability voting approach. J. Eng. Ind. 1990, 112, 319–325. [Google Scholar] [CrossRef]
  29. Chen, Y.; Orady, E. An entropy-based index evaluation scheme for multiple sensor fusion in classification process. J. Manuf. Sci. Eng. 1999, 121, 727–732. [Google Scholar] [CrossRef]
  30. Liu, R.; Yang, B.; Zio, E.; Chen, X. Artificial intelligence for fault diagnosis of rotating machinery: A review. Mech. Syst. Signal Process. 2018, 108, 33–47. [Google Scholar] [CrossRef]
  31. Luo, R.C.; Kay, M.G. A tutorial on multisensor integration and fusion. In Proceedings of the IECON’ 90: 16th Annual Conference of IEEE Industrial Electronics Society, Pacific Grove, CA, USA, 27–30 November 1990; pp. 707–722. [Google Scholar]
  32. Staszewski, W. Structural and mechanical damage detection using wavelets. Shock Vib. Dig. 1998, 30, 457–472. [Google Scholar] [CrossRef]
  33. Mao, Z.; Jiang, B. Fault identification and fault-tolerant control for a class of networked control systems. Int. J. Innov. Comput. Inf. Control 2007, 3, 1121–1130. [Google Scholar]
  34. Banerjee, T.P.; Das, S. Multi-sensor data fusion using support vector machine for motor fault detection. Inf. Sci. 2012, 217, 96–107. [Google Scholar] [CrossRef]
  35. Chang, L.; Hu, B.; Li, A.; Qin, F. Unscented type Kalman filter: Limitation and combination. IET Signal Process. 2013, 7, 167–176. [Google Scholar] [CrossRef]
  36. Julier, S.J.; Uhlmann, J.K. Unscented filtering and nonlinear estimation. Proc. IEEE 2004, 92, 401–422. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, Q. Adaptive observer for multiple-input-multiple-output (MIMO) linear time-varying systems. IEEE Trans. Autom. Control 2002, 47, 525–529. [Google Scholar] [CrossRef]
  38. Yunusa-Kaltungo, A.; Sinha, J.K.; Elbhbah, K. An improved data fusion technique for faults diagnosis in rotating machines. Measurement 2014, 58, 27–32. [Google Scholar] [CrossRef]
  39. Yunusa-Kaltungo, A.; Sinha, J.K.; Nembhard, A.D. A novel fault diagnosis technique for enhancing maintenance and reliability of rotating machines. Struct. Health Monit. 2015, 14, 604–621. [Google Scholar] [CrossRef] [Green Version]
  40. Yunusa-Kaltungo, A.; Sinha, J.K. Sensitivity analysis of higher order coherent spectra in machine faults diagnosis. Struct. Health Monit. 2016, 15, 555–567. [Google Scholar] [CrossRef] [Green Version]
  41. Yunusa-Kaltungo, A.; Sinha, J.K.; Nembhard, A.D. Use of composite higher order spectra for faults diagnosis of rotating machines with different foundation flexibilities. Measurement 2015, 70, 47–61. [Google Scholar] [CrossRef]
  42. Yunusa-Kaltungo, A.; Cao, R. Towards Developing an Automated Faults Characterisation Framework for Rotating Machines. Part 1: Rotor-Related Faults. Energies 2020, 13, 1394. [Google Scholar] [CrossRef] [Green Version]
  43. Cao, P.; Zhang, S.; Tang, J. Preprocessing-free gear fault diagnosis using small datasets with deep convolutional neural network-based transfer learning. IEEE Access 2018, 6, 26241–26253. [Google Scholar] [CrossRef]
  44. Shao, S.; McAleer, S.; Yan, R.; Baldi, P. Highly accurate machine fault diagnosis using deep transfer learning. IEEE Trans. Ind. Inform. 2018, 15, 2446–2455. [Google Scholar] [CrossRef]
  45. Soualhi, M.; Nguyen, K.T.; Medjaher, K. Pattern recognition method of fault diagnostics based on a new health indicator for smart manufacturing. Mech. Syst. Signal Process. 2020, 142, 106680. [Google Scholar] [CrossRef]
  46. Azamfar, M.; Singh, J.; Bravo-Imaz, I.; Lee, J. Multisensor data fusion for gearbox fault diagnosis using 2-D convolutional neural network and motor current signature analysis. Mech. Syst. Signal Process. 2020, 144, 106861. [Google Scholar] [CrossRef]
  47. Zhang, Z.; Li, S.; Lu, J.; Wang, J.; Jiang, X. A novel intelligent fault diagnosis method based on fast intrinsic component filtering and pseudo-normalization. Mech. Syst. Signal Process. 2020, 145, 106923. [Google Scholar] [CrossRef]
  48. Scheffer, C.; Girdhar, P. Practical Machinery Vibration Analysis and Predictive Maintenance; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  49. Luwei, K.C.; Sinha, J.K.; Yunusa-Kaltungo, A. Comparison of Amplitude to Real and Imaginary Features of the poly-Coherent Composite Bispectrum (pCCB) Components in Machine Diagnosis. In Advances in Asset Management and Condition Monitoring; Springer: Cham, Switzerland, 2020; pp. 1–8. [Google Scholar]
  50. Luwei, K.C.; Yunusa-Kaltungo, A.; Sha’aban, Y.A. Integrated fault detection framework for classifying rotating machine faults using frequency domain data fusion and artificial neural networks. Machines 2018, 6, 59. [Google Scholar] [CrossRef] [Green Version]
  51. Luwei, K.C.; Sinha, J.K.; Yunusa-Kaltungo, A.; Elbhbah, K. Data fusion of acceleration and velocity features (dFAVF) approach for fault diagnosis in rotating machines. InMATEC Web Conf. 2018, 211, 21005. [Google Scholar] [CrossRef] [Green Version]
  52. Qiao, Q.; Yunusa-Kaltungo, A.; Edwards, R.E. Towards developing a systematic knowledge trend for building energy consumption prediction. J. Build. Eng. 2020, 5, 101967. [Google Scholar] [CrossRef]
  53. Qiao, Q.; Yunusa-Kaltungo, A.; Edwards, R. Predicting building energy consumption based on meteorological data. In Proceedings of the 2020 IEEE PES/IAS PowerAfrica 2020, Nairobi, Kenya, 25–28 August 2020; IEEE: Piscataway, NJ, USA, 2020; pp. 1–5. [Google Scholar]
  54. Qiao, Q.; Yunusa-Kaltungo, A.; Edwards, R. Hybrid method for building energy consumption prediction based on limited data. In Proceedings of the 2020 IEEE PES/IAS PowerAfrica 2020, Nairobi, Kenya, 25–28 August 2020; IEEE: Piscataway, NJ, USA, 2020; pp. 1–5. [Google Scholar]
Figure 1. Schematic representation of the faults classification process.
Figure 1. Schematic representation of the faults classification process.
Sensors 21 02957 g001
Figure 2. Experimental rig.
Figure 2. Experimental rig.
Sensors 21 02957 g002
Figure 3. Experimentally simulated cases: (a) notch on driver gear tooth; (b) worn-off surfaces on driver gear teeth.
Figure 3. Experimentally simulated cases: (a) notch on driver gear tooth; (b) worn-off surfaces on driver gear teeth.
Sensors 21 02957 g003
Figure 4. Experimental flowchart.
Figure 4. Experimental flowchart.
Sensors 21 02957 g004
Figure 5. Harmonic distributions for all cases at all speeds: (a) rotor-related harmonics at 7 Hz; (b) rotor-related harmonics at 14 Hz; (c) rotor-related harmonics at 21 Hz; (d) gear mesh harmonics at 7 Hz; (e) gear mesh harmonics at 14 Hz; (f) gear mesh harmonics at 21 Hz.
Figure 5. Harmonic distributions for all cases at all speeds: (a) rotor-related harmonics at 7 Hz; (b) rotor-related harmonics at 14 Hz; (c) rotor-related harmonics at 21 Hz; (d) gear mesh harmonics at 7 Hz; (e) gear mesh harmonics at 14 Hz; (f) gear mesh harmonics at 21 Hz.
Sensors 21 02957 g005
Figure 6. Spectral energy (SE) based gear mesh frequencies harmonics: (a) 7 Hz; (b) 14 Hz; (c) 21 Hz.
Figure 6. Spectral energy (SE) based gear mesh frequencies harmonics: (a) 7 Hz; (b) 14 Hz; (c) 21 Hz.
Sensors 21 02957 g006
Figure 7. PCs combinations for all scenarios: (ac) PC1–3 for (α1–α5) + (β1–β5) at 7, 14, and 21 Hz, respectively; (df) PC1–2 for (α1–α5) + (β1–β5) at 7, 14, and 21 Hz, respectively; (gi) PC1–2 for (α1–α5) at 7, 14, and 21 Hz, respectively; (jl) PC1–2 for (β1–β5) at 7, 14, and 21 Hz respectively; (mo) PC1–2 for (γ1–γ5) at 7, 14, and 21 Hz, respectively.
Figure 7. PCs combinations for all scenarios: (ac) PC1–3 for (α1–α5) + (β1–β5) at 7, 14, and 21 Hz, respectively; (df) PC1–2 for (α1–α5) + (β1–β5) at 7, 14, and 21 Hz, respectively; (gi) PC1–2 for (α1–α5) at 7, 14, and 21 Hz, respectively; (jl) PC1–2 for (β1–β5) at 7, 14, and 21 Hz respectively; (mo) PC1–2 for (γ1–γ5) at 7, 14, and 21 Hz, respectively.
Sensors 21 02957 g007
Figure 8. Classification accuracies for a range of k values of k-NN.
Figure 8. Classification accuracies for a range of k values of k-NN.
Sensors 21 02957 g008
Figure 9. Visualised decision maps for ANN1–ANN3 using PC1–2 as inputs for all cases: (ac) (α1–α5) + (β1–β5) at 7 Hz; (df) (α1–α5) + (β1–β5) at 14 Hz; (gi) (α1–α5) + (β1–β5) at 21 Hz; (jl) (γ1–γ5) at 7 Hz; (mo) (γ1–γ5) at 14 Hz; (pr) (γ1–γ5) at 21 Hz.
Figure 9. Visualised decision maps for ANN1–ANN3 using PC1–2 as inputs for all cases: (ac) (α1–α5) + (β1–β5) at 7 Hz; (df) (α1–α5) + (β1–β5) at 14 Hz; (gi) (α1–α5) + (β1–β5) at 21 Hz; (jl) (γ1–γ5) at 7 Hz; (mo) (γ1–γ5) at 14 Hz; (pr) (γ1–γ5) at 21 Hz.
Sensors 21 02957 g009
Figure 10. PC1–2 and visualised decision maps for ANN1 using PC1–2 as inputs for (α1–α5) of validation datasets at 20 and 30 Hz: (a) PC1–2 at 20 Hz; (b) PC1–2 at 30 Hz; (c) ANN1 at 20 Hz; (d) ANN1 at 30 Hz.
Figure 10. PC1–2 and visualised decision maps for ANN1 using PC1–2 as inputs for (α1–α5) of validation datasets at 20 and 30 Hz: (a) PC1–2 at 20 Hz; (b) PC1–2 at 30 Hz; (c) ANN1 at 20 Hz; (d) ANN1 at 30 Hz.
Sensors 21 02957 g010
Table 1. Comparison between the proposed framework and the related state-of-the-art works.
Table 1. Comparison between the proposed framework and the related state-of-the-art works.
ReferenceData TypeClassification AlgorithmUse of Data FusionFault Classes ConsideredLimitation
[38,39,40,41]Vibration dataPCA and spectrum-based liner classificationYesRotor faultsClassification approach used does not involve machine learning, thereby making the approach unable to self-learned from historical data. Additionally, all the simulated cases are rotor-related.
[42]Vibration dataANN, SVM, k-NN, naïve BayesYesRotor faultsAll the simulated cases are rotor-related.
[43]Vibration dataCNNNoGearsComputational intensiveness due to the application of deep learning approach. This study only considered a single fault class. Hence, there might be need to further investigate the applicability on other faults.
[44]Vibration dataCNNNoMotor, gears and bearing faultsAlthough this study considered multiple fault classes, the approach is also computationally intensive and would lead to costlier solutions.
[45]Current dataANFISNoGearsThis study will benefit from the implementation of multi-sensor data fusion (such as electrical and mechanical data), so as to enhance the reliability of fault diagnosis.
[46]Current dataCNNYesGearsComputational intensiveness due to the application of deep learning approach. This study only considered a single fault class. Hence, there might be a need to further investigate the applicability for other faults.
[47]Vibration dataFICFNoGears and bearing faultsFICF is suitable for multi-sample training but the convolution activation limits its performance during single sample operations. Additionally, FICF is often considered a high-efficiency technique, but the poor noise adaptability sometimes undermines its proficiency.
Present studyVibration dataANNYesRotor and gear faultsThe only limitation envisaged with the current approach is its initial requirement for training data for different fault types, which is also common to most of the aforementioned techniques.
Table 2. Technical specifications of main rig components.
Table 2. Technical specifications of main rig components.
Rig ComponentDescription (Abbreviation)Specification/Dimension
Electric motorHorsepower2.0
Speed3600 RPM
Maximum torque3.9 Nm
TypeShunt
Current6.8 Amps
Voltage200 V
ShaftsTypeMild steel
Length1 040 mm
Diameter35 mm
GearsTypeHelical (key-mounted)
Pitch circle diameter (PCD)107 mm
Face width37 mm
Circular thickness4.5 mm
Number of teeth35
Addendum2.35 mm
Dedendum1.95 mm
BearingsTypeAnti-friction ball bearings
MakeSKF
ModelSY20TF/RA SEY20/NP20
Number of rolling elements8
Diameter of rolling elements7.938 mm
Bearing width31 mm
External diameter47 mm
Internal diameter20 mm
Bearing pitch circle diameter33.5 mm
PulleysTypeToothed, taper lock-mounted
Thickness32 mm
Tooth thickness6.8 mm
Drive pulley diameter100 mm
Driven pulley diameter125 mm
BeltsTypeToothed, timing
ModelFenner 200H-100
Table 3. Technical specifications of main instruments.
Table 3. Technical specifications of main instruments.
VM InstrumentParameterSpecification
AccelerometerModel No.352C33
Sensitivity (±10%)100 mV/g
Frequency range (±5%)0.5 to 10,000 Hz
Resonant frequency≥50 kHz
Temperature range−65 to +200 °F
Settling time (within 10% of bias)<10 s
Signal conditionerInput sensor typeICP, voltage, charge
Voltage gain×0.1 to ×200
Voltage gain increment0.1
Charge conversion (selectable)0.1, 1.0, 10.0 mV/pC
Frequency range (gain <100)0.05 to 100 kHz
Frequency range (gain 100)0.05 to 75 kHz
ADCNumber of channels16 differential/32 single ended
ADC resolution16 bits
Sampling rate250 kS/s single channel; 250 kS/s multi-channel (aggregate)
Input range±10 V, ±5 V; ±1 V, ±0.2 V
Input FIFO size4095 samples
Table 4. Percentage of explained variance by each principal component (%).
Table 4. Percentage of explained variance by each principal component (%).
Feature1–α5) + (β1–β5)1–α5)1–β5)1–γ5)
Frequency7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz
PC162.657657.237429.939868.708755.750043.111167.691468.291238.735572.214975.004851.7278
PC27.806512.069215.637511.665214.006619.137513.859918.489322.104513.675717.849321.7150
PC36.44516.935410.96019.363112.043216.49819.70196.287216.78078.77333.127411.5762
PC45.57885.64639.51305.76209.899910.95784.92163.927613.06653.02542.42688.2653
PC54.56854.55338.17504.50108.300310.29543.82523.00479.31272.31071.59166.7156
PC64.09544.24806.3779------------------------------------------------------
PC72.70943.26985.5352------------------------------------------------------
PC82.29922.73405.1251------------------------------------------------------
PC92.06041.87714.7339------------------------------------------------------
PC101.77911.42954.0022------------------------------------------------------
PC1–270.464169.306645.577380.373969.756662.248681.551386.780560.840085.890692.854173.4428
PC1–376.909276.242056.537489.737081.799878.746791.253293.067777.620794.663995.981585.0190
Table 5. ANN properties for (α1–α5) + (β1–β5) using PC1–3 and without PCA.
Table 5. ANN properties for (α1–α5) + (β1–β5) using PC1–3 and without PCA.
ParametersANN1ANN2ANN3ANN4
Network STRUCTURE3–10–33–20–33–30–310–30–3 (without PCA)
Rotation Frequency7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz
Accuracy (%)Training99.198.898.299.199.098.299.399.398.598.910098.7
Validation99.598.598.599.599.598.599.599.598.599.510099.0
Testing99.599.098.599.599.599.099.510099.098.510097.0
Overall99.298.898.399.299.298.399.399.498.698.910098.5
Table 6. ANN properties for (α1–α5) + (β1–β5) using PC1–2 and without PCA.
Table 6. ANN properties for (α1–α5) + (β1–β5) using PC1–2 and without PCA.
ParametersANN1ANN2ANN3ANN4
Network Structure2–10–32–20–32–30–310–30–3 (without PCA)
Rotation Frequency7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz
Accuracy (%)Training98.999.195.198.999.195.398.999.195.498.910098.7
Validation99.099.095.099.099.095.599.099.095.599.510099.0
Testing99.099.095.099.099.094.599.599.093.598.510097.0
Overall98.999.195.098.999.195.299.099.195.198.910098.5
Table 7. ANN properties for (γ1–γ5) using PC1–2 and without PCA.
Table 7. ANN properties for (γ1–γ5) using PC1–2 and without PCA.
ParametersANN1ANN2ANN3ANN4
Network Structure2–10–32–20–32–30–35–30–3 (without PCA)
Rotation Frequency7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz7 Hz14 Hz21 Hz
Accuracy (%)Training99.410088.299.510089.399.510089.399.510097.5
Validation99.510089.099.510088.599.510089.599.010096.5
Testing99.510089.099.510090.599.510087.599.510097.0
Overall99.410088.499.510089.399.510089.099.410097.3
Table 8. Classifier properties for (α1–α5) + (β1–β5), and (γ1–γ5) using PC1–2.
Table 8. Classifier properties for (α1–α5) + (β1–β5), and (γ1–γ5) using PC1–2.
ParametersFeaturesRotation FrequencyANN1 (2–10–3)k-NN (k = 10)Naïve BayesLinear SVM
Accuracy (%)1–α5) + (β1–β5)7 Hz98.998.998.598.8
14 Hz99.198.898.598.7
21 Hz95.094.391.594.3
1–γ5)7 Hz99.499.498.099.3
14 Hz10010099.7100
21 Hz88.488.387.886.7
Table 9. Description of validation gearbox fault types [44].
Table 9. Description of validation gearbox fault types [44].
TypeDescription
ChippedCrack occurs in the gear teeth
MissMissing one gear tooth
RootCrack occurs in the root of the gear
SurfaceWear occurs in the surface of gear
Table 10. Percentage of explained variance by each principal component for (α1–α5) of public datasets (%).
Table 10. Percentage of explained variance by each principal component for (α1–α5) of public datasets (%).
Working Condition20 Hz–0 V30 Hz–2 V
Measuring PositionP1xyz + P2xyzP1x + P2xP1y + P2yP1z + P2zP1xyz + P2xyzP1x + P2xP1y + P2yP1z + P2z
PC154.905655.645941.089665.245742.423154.284239.287248.8888
PC219.559421.086324.545217.661327.536824.771623.009825.1579
PC313.845814.440719.534414.205113.120810.294417.976215.9588
PC410.42567.260713.34571.98479.87418.093811.52198.3389
PC51.26361.56641.48510.90327.04522.55608.20491.6557
PC1–274.465076.732265.634882.907069.959979.055862.297074.0467
Table 11. ANN1 (2–10–5) properties for (α1–α5) of public datasets using PC1–2.
Table 11. ANN1 (2–10–5) properties for (α1–α5) of public datasets using PC1–2.
Working Condition20 Hz–0 V30 Hz–2 V
Measuring PositionP1xyz + P2xyzP1x + P2xP1y + P2yP1z + P2zP1xyz + P2xyzP1x + P2xP1y + P2yP1z + P2z
Accuracy (%)Training10010091.510010099.587.596.7
Validation10010092.010010099.586.696.8
Testing10010090.499.510099.586.695.2
Overall10010091.499.910099.587.296.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cao, R.; Yunusa-Kaltungo, A. An Automated Data Fusion-Based Gear Faults Classification Framework in Rotating Machines. Sensors 2021, 21, 2957. https://doi.org/10.3390/s21092957

AMA Style

Cao R, Yunusa-Kaltungo A. An Automated Data Fusion-Based Gear Faults Classification Framework in Rotating Machines. Sensors. 2021; 21(9):2957. https://doi.org/10.3390/s21092957

Chicago/Turabian Style

Cao, Ruifeng, and Akilu Yunusa-Kaltungo. 2021. "An Automated Data Fusion-Based Gear Faults Classification Framework in Rotating Machines" Sensors 21, no. 9: 2957. https://doi.org/10.3390/s21092957

APA Style

Cao, R., & Yunusa-Kaltungo, A. (2021). An Automated Data Fusion-Based Gear Faults Classification Framework in Rotating Machines. Sensors, 21(9), 2957. https://doi.org/10.3390/s21092957

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