A Fault Diagnosis Model for Coaxial-Rotor Unit Using Bidirectional Gate Recurrent Unit and Highway Network

: A turbojet engine is the most signiﬁcant part of an Internal Combustion Engine (ICE) for Hybrid Electric Vehicles. Speciﬁcally, the coaxial-rotor unit is the key component, whose performance largely affects the working efﬁciency. Thereby, the fault diagnosis methods for coaxial-rotor units is a main focus. In line with our test results, the bearing circlip is the most vulnerable element while rotating. Moreover, the low-speed rotating fault diagnosis is even challenging for current methods. Since the fault diagnosis on the bearing circlip of coaxial-rotor units is absent, this paper establishes a test rig on a running coaxial-rotor unit under different working conditions. The three-directional vibration signals are collected and analyzed to demonstrate the working states. On the task of bearing circlip failure classiﬁcation, a deep-learning-based model using the Bidirectional Gate Recurrent Unit and the Highway Network is developed, which is capable of capturing hidden features and removing unrelated information. For working performance evaluation, experiments on the data of different rotating speeds are carried out. Among all the fault diagnosis methods, our model is the best approach and achieves an average accuracy of 99.4%. The encouraging results reveal that the proposed model is effective in both the high-speed and low-speed fault diagnosis of bearing circlip malfunction. The three-direction vibration signals are measured and dedicatedly analyzed to characterize the working status. The HHT analysis is performed to study the failure mechanism of the bearing circlips. In line with the property of bearing circlip failure, a fault classiﬁcation method, based on the integration of Bi-GRU and the Highway Network, is proposed. Basically, the GRU cells are capable of capturing the hidden information from inputs, while the Highway Network removes the unrelated features. The fault diagnosis performance of the proposed model is evaluated on nine different working conditions. An average accuracy of 99.4% of the proposed model is obtained, which outperforms the state-of-the-art methods. Speciﬁcally, our model also shows its superiority in low-speed fault classiﬁcation. The experimental results indicate that our model obtains sufﬁcient classiﬁcation accuracy. That is, our model shows its distinctive superiority in dealing with the detected sensing signals precisely and thus can be applied to coaxial-rotor unit fault diagnosis tasks.


Introduction
The automobile industry faces the challenge of saving energy and eliminating environmental damage, whilst delivering more goods and passengers to achieve a longer driving range [1]. In line with the ongoing advancements, electric vehicles, whose primitive form is proposed in 1900 s, have caused a revival of interest in vehicle engineering [2]. Based on the developing potential of electric vehicles and batteries, pure electric vehicles will ultimately meet the demands of high energy efficiency and low-pollutant emission. Despite its high energy efficiency and low-pollutant emission, the pure electric vehicle is currently a secondary choice for long-range driving restricted to battery storage capacity [3]. That is, these vehicle types are not likely to be widely used within the short term. Encouragingly, Hybrid Electric Vehicles (HEVs) bridge the gap between conventional vehicles and pure electric vehicles, which are viable means of transportation at present [4]. An HEV has an Internal Combustion Engine (ICE) as the main power source and an electric drive (i.e., a battery and an electric motor) to improve fuel economy [5]. Specifically, the battery required to power the electric motor is charged by the ICE without any external source. As such, almost all auto makers provide a hybrid version of their popular models such as Audi A8 Hybrid, Ford Fusion, Honda Civic Hybrid, and Hyundai Sonata Hybrid [6].
Typically, a turbojet engine, integrating with an electric generator, is employed in an ICE to supply stable and reliable power [7]. As long as the ICE aims to supply the (1) A test rig for coaxial-rotor unit fault generation is established, based on which the working properties under different conditions are analyzed. Specifically, the bearing circlip failure is pronounced. (2) A novel fault diagnosis method is established, integrating the Bidirectional Gated Recurrent Unit and Highway Network. Our model aims to not just capture the internal features of the working state but also identify the faults of the coaxial-rotor unit with various rotating speeds. (3) In line with the practical working conditions, experiments are carried out on coaxialrotor unit vibration signals to verify the effectiveness of the proposed model. Our model achieves decent working performance in the test. Figure 1 presents the schematic of a turbojet engine. The major components of the engine are illustrated. The air from upstream comes through the centrifugal impeller of the compressor and reaches the diffuser. Via both the axial and radial diffusing, the air enters into the combustion chamber by passing the airhole. The combustion process is triggered by the ignition of the igniter, with the air mixing with the fuel. Hereafter, the high-pressure and high-temperature gas expands towards the turbine, which provides power to drive the spindle, the turbine, and the compressor rotating at high speed. Subsequently, the gas goes into the jet nozzle at lower pressure and a higher speed. In this way, the thrust is produced whilst the exhaust is discharged to the atmosphere [21]. spindle, the high-speed rotors, and the bearings. Apparently, in this system, the rotors are mainly the turbine and the compressor. Meanwhile, there is a ball bearing on each end of the rotor. Generally, the bearing is made up of high-speed ceramic with the inner and outer ring of bearing steel. An example of a coaxial-rotor unit is exhibited in Figure 3. The rotating speed of the coaxial-rotor unit is over 10,000 r/min, while the highest speed reaching 100,000 r/min. Due to the working principle of the turbojet engine, the bearings are set to protect the rotors and, in turn, are most likely to break down during working. Based on the previous investigations, the bearing failure is the most probable failure mode in practical use [22][23][24].  We shall thus detach the structure of turbojet engine to resolve the coaxial-rotor unit. According to Figure 2, regardless of the outer case, the coaxial-rotor unit consists of a spindle, the high-speed rotors, and the bearings. Apparently, in this system, the rotors are mainly the turbine and the compressor. Meanwhile, there is a ball bearing on each end of the rotor. Generally, the bearing is made up of high-speed ceramic with the inner and outer ring of bearing steel. An example of a coaxial-rotor unit is exhibited in Figure 3. The rotating speed of the coaxial-rotor unit is over 10,000 r/min, while the highest speed reaching 100,000 r/min. Due to the working principle of the turbojet engine, the bearings are set to protect the rotors and, in turn, are most likely to break down during working. Based on the previous investigations, the bearing failure is the most probable failure mode in practical use [22][23][24].

Structure of Coaxial-Rotor Unit
Machines 2022, 10, x FOR PEER REVIEW 3 of 16 triggered by the ignition of the igniter, with the air mixing with the fuel. Hereafter, the high-pressure and high-temperature gas expands towards the turbine, which provides power to drive the spindle, the turbine, and the compressor rotating at high speed. Subsequently, the gas goes into the jet nozzle at lower pressure and a higher speed. In this way, the thrust is produced whilst the exhaust is discharged to the atmosphere [21]. We shall thus detach the structure of turbojet engine to resolve the coaxial-rotor unit. According to Figure 2, regardless of the outer case, the coaxial-rotor unit consists of a spindle, the high-speed rotors, and the bearings. Apparently, in this system, the rotors are mainly the turbine and the compressor. Meanwhile, there is a ball bearing on each end of the rotor. Generally, the bearing is made up of high-speed ceramic with the inner and outer ring of bearing steel. An example of a coaxial-rotor unit is exhibited in Figure 3. The rotating speed of the coaxial-rotor unit is over 10,000 r/min, while the highest speed reaching 100,000 r/min. Due to the working principle of the turbojet engine, the bearings are set to protect the rotors and, in turn, are most likely to break down during working. Based on the previous investigations, the bearing failure is the most probable failure mode in practical use [22][23][24].

System Establishment
As pointed out in the Introduction, the testing on the coaxial-rotor unit of the turbojet engine is a complex and interrelated process. For the task of fault diagnosis, the actual coaxial-rotor unit is built up and operated in various working states; see Figure 4. Note that we aim to concentrate on the condition of the coaxial-rotor unit, other components of the turbojet engine, including the combustion chamber and the ignition devices, are removed. The outlook of the coaxial-rotor unit in our experiment is presented in Figure 5a.
Considering the impact from the combustion, high-temperature gas is generated and insufflated to simulate the practical working environment. The gas generator is attached to the coaxial-rotor unit via the pipelines seen in Figure 5b. Notably, the rotating speed varies due to the control of the air-gas ratio inside the gas generator. That is, by adjusting the proportional valve of the flow meter, the inputs of gas and air are regulated via a controlling module (Figure 5d), and thus their mixing ratio is revised. In this way, the rotation of the coaxial-rotor can be determined. For instance, the rotating speed reaches 20,000 r/min at an air-gas ratio of 26.8:4.5. A lubricating oil condition unit and a cooling water condition unit are employed to facilitate the system operation. Besides, the temperature of the exhaust is set as 50 • C at the initial stage.

System Establishment
As pointed out in the Introduction, the testing on the coaxial-rotor unit of the turbojet engine is a complex and interrelated process. For the task of fault diagnosis, the actual coaxial-rotor unit is built up and operated in various working states; see Figure 4. Note that we aim to concentrate on the condition of the coaxial-rotor unit, other components of the turbojet engine, including the combustion chamber and the ignition devices, are removed. The outlook of the coaxial-rotor unit in our experiment is presented in Figure 5a. Considering the impact from the combustion, high-temperature gas is generated and insufflated to simulate the practical working environment. The gas generator is attached to the coaxial-rotor unit via the pipelines seen in Figure 5b. Notably, the rotating speed varies due to the control of the air-gas ratio inside the gas generator. That is, by adjusting the proportional valve of the flow meter, the inputs of gas and air are regulated via a controlling module (Figure 5d), and thus their mixing ratio is revised. In this way, the rotation of the coaxial-rotor can be determined. For instance, the rotating speed reaches 20,000 r/min at an air-gas ratio of 26.8:4.5. A lubricating oil condition unit and a cooling water condition unit are employed to facilitate the system operation. Besides, the temperature of the exhaust is set as 50 °C at the initial stage.
Sensing elements for detecting the working parameters are installed alongside the working components within the system. In the domain of measuring the rotating machinery, the vibration signal is most pronounced for characterizing the operational status directly. According to statistics, more than 70% of faults are manifested in the form of vibration [25]. In this way, to diagnose the potential faults, a three-axis vibrating acceleration sensor is attached to the rotor case. The specific installation position and sensor structure are shown in Figure 5c. The vibration of three directions-axial, radial, and tangentialis detected during working phases. With the power supplied, the vibrating signals are gathered by Integrated Circuits Piezoelectric (ICP) sensors 356A33 designed by PCB Piezotronics Inc. [26]. Sensing signals are transmitted through the 10-m low-noise cables. The data collecting and recording are implemented on an LMS830A analyzer. The measurement frequency is as high as 50 k Hz for dynamic signal detection. The signal acquisition time is 120 s and involves a whole work phase of the coaxial-rotor unit. The three-directional vibration signals are displayed in real time and kept for further analysis. In addi-    Sensing elements for detecting the working parameters are installed alongside the working components within the system. In the domain of measuring the rotating machinery, the vibration signal is most pronounced for characterizing the operational status directly. According to statistics, more than 70% of faults are manifested in the form of vibration [25]. In this way, to diagnose the potential faults, a three-axis vibrating acceleration sensor is attached to the rotor case. The specific installation position and sensor structure are shown in Figure 5c. The vibration of three directions-axial, radial, and tangential-is detected during working phases. With the power supplied, the vibrating signals are gathered by Integrated Circuits Piezoelectric (ICP) sensors 356A33 designed by PCB Piezotronics Inc. [26]. Sensing signals are transmitted through the 10-m low-noise cables. The data collecting and recording are implemented on an LMS830A analyzer. The measurement frequency is as high as 50 k Hz for dynamic signal detection. The signal acquisition time is 120 s and involves a whole work phase of the coaxial-rotor unit. The three-directional vibration signals are displayed in real time and kept for further analysis. In addition, a tungsten-rhenium thermocouple temperature sensor measures the temperature of the hot gas whilst a Hall speed sensor monitors the rotating. Figure 5 shows the photograph of the major components of the system. The specifications of the corresponding measurement devices are given in Table 1.  We find that the circlip of the bearing is the most vulnerable part during working; see Figure 6. Notably, according to recent publications, a great deal of attention is paid to the fault diagnosis of bearings of rotating machinery [27]. Research on bearing circlip failures is currently absent. Hence, we mainly focus on the identification of the bearing circlip malfunction of the coaxial-rotor unit in this work.
We find that the circlip of the bearing is the most vulnerable part during working; see Figure 6. Notably, according to recent publications, a great deal of attention is paid to the fault diagnosis of bearings of rotating machinery [27]. Research on bearing circlip failures is currently absent. Hence, we mainly focus on the identification of the bearing circlip malfunction of the coaxial-rotor unit in this work.

Preliminary Analysis
Generally, the coaxial-rotor unit rotating speed of turbojet engine ranges from 20,000 r/min to 60,000 r/min, with an interval of 5000 r/min. The working-states-based varyingspeed signals can be measured. As an example, the time-domain vibration of the three dimensions in the normal state and bearing circlip failure, with a speed of 40,000 r/min, is illustrated in Figure 7a and 7b, respectively. Firstly, the vibrating data properties are investigated to describe the damage mechanism. For rotating machinery, the output power varies within the natural frequency band

Preliminary Analysis
Generally, the coaxial-rotor unit rotating speed of turbojet engine ranges from 20,000 r/min to 60,000 r/min, with an interval of 5000 r/min. The working-states-based varying-speed signals can be measured. As an example, the time-domain vibration of the three dimensions in the normal state and bearing circlip failure, with a speed of 40,000 r/min, is illustrated in Figures 7a and 7b, respectively.
We find that the circlip of the bearing is the most vulnerable part during working; see Figure 6. Notably, according to recent publications, a great deal of attention is paid to the fault diagnosis of bearings of rotating machinery [27]. Research on bearing circlip failures is currently absent. Hence, we mainly focus on the identification of the bearing circlip malfunction of the coaxial-rotor unit in this work.

Preliminary Analysis
Generally, the coaxial-rotor unit rotating speed of turbojet engine ranges from 20,000 r/min to 60,000 r/min, with an interval of 5000 r/min. The working-states-based varyingspeed signals can be measured. As an example, the time-domain vibration of the three dimensions in the normal state and bearing circlip failure, with a speed of 40,000 r/min, is illustrated in Figure 7a and 7b, respectively. Firstly, the vibrating data properties are investigated to describe the damage mechanism. For rotating machinery, the output power varies within the natural frequency band Firstly, the vibrating data properties are investigated to describe the damage mechanism. For rotating machinery, the output power varies within the natural frequency band during the damage, which indicates the fault condition. At this stage, the Hilbert marginal spectrum is established to analyze the vibrating signals of various rotating speeds. To the best of our knowledge, the HHT is capable of tackling structure damage data in a variety of tasks [28,29].
Basically, the HHT consists of two main procedures, namely, Empirical Mode Decomposition (EMD) and Hilbert Spectrum Analysis (HSA). The EMD technique is taken to decompose a given signal into a number of Intrinsic Mode Functions (IMF) that provide much time-frequency information to the narrowband [30]. Each IMF component contains local features of different time scales and meets certain conditions, based on which fault features can be extracted [31]. Subsequently, the Hilbert transform is carried out on each IMF to obtain the corresponding Hilbert spectrum. That is, these IMFs are conveyed into the time-frequency domain. Then, the IMFs' entire Hilbert spectrum is aggregated to order the marginal spectrum derived in feature extraction [32]. In such a manner, the Hilbert marginal spectrum can be further used in fault diagnosis.
Let x(t) be the vibration signals of the coaxial-rotor unit. The EMD is employed to decompose x(t) to obtain a set of IMFs, which can be written as: where c i (t) stands for the ith IMF; n is the number of IMF components; and r n (t) is the residual function. Then, the HHT is applied to each c i (t), so we have From Equation (2), H[c i (t)] and c i (t) can be integrated to obtain where and The parameters A i (t), ϕ i (t), and w i (t) represent the instantaneous amplitude, instantaneous phase, and instantaneous frequency of the ith IMF, respectively.
Considering that r n (t) is either a monotonic function or a constant, and no IMF component can be extracted from r n (t), the original signal can be rewritten as: The time-frequency-amplitude distribution is described using Hilbert spectrum, which is: Further, the Hilbert marginal spectrum is obtained by the time integral of the Hilbert spectrum (Equation (9)), which delivers the amplitude distribution of each instantaneous frequency.
Machines 2022, 10, 313 8 of 16 Based on the Hilbert marginal spectrum, we define the Hilbert marginal spectrum energy as: According to Equation (10), we see that the signal frequency characteristics is magnified in E(w).
Since the vibration data of three directions are collected, the Hilbert marginal spectrum energy of the three-axis signal is the integration all these energies: where ∂ i is the weighting coefficient of the ith Hilbert marginal spectrum energy and satisfies: The energy can be normalized to: where E k represents the marginal spectrum energy value of the kth IMF component; p k is the energy normalization value; and m is the order of the IMF. Specifically, the Hilbert marginal spectrum energy of different rotating speeds is presented in Figure 8.
Based on the Hilbert marginal spectrum, we define the Hilbert marginal spectrum energy as: According to Equation (10), we see that the signal frequency characteristics is magnified in ( ).
Since the vibration data of three directions are collected, the Hilbert marginal spectrum energy of the three-axis signal is the integration all these energies: where is the weighting coefficient of the ℎ Hilbert marginal spectrum energy and satisfies: The energy can be normalized to: where represents the marginal spectrum energy value of the ℎ IMF component; is the energy normalization value; and is the order of the IMF. Specifically, the Hilbert marginal spectrum energy of different rotating speeds is presented in Figure 8. In Figure 8, one can easily see that the marginal spectral energy drops dramatically with the increase of IMF order. Notably, the marginal spectral energies of IMF2, IMF3, In Figure 8, one can easily see that the marginal spectral energy drops dramatically with the increase of IMF order. Notably, the marginal spectral energies of IMF2, IMF3, and IMF4 are negligible since their values are smaller than 10-5. For IMF1, the Hilbert marginal spectrum energies of the normal state and bearing circlip failure are distinguished when the rotating speed reaches 35,000 r/min. We thus infer that the fault identification at a low speed is challenging due to the similarity of marginal spectrum energy distribution.

Model Architecture
Aiming to identify the working state of the coaxial-rotor unit under different conditions, a fault diagnosis approach based on a Bidirectional Gated Recurrent Unit (Bi-GRU) and the Highway Network is proposed. The architecture of our model is presented in Figure 9. and the Highway Network is proposed. The architecture of our model is presented in Figure 9.
Let , , and represent the set of axial, radial, and tangential vibrating signals, respectively. [ , , ] ∈ can be the inputs of the proposed fault diagnosis model. Each input vector is sent to a fully connected layer for feature extraction. Specifically, the active function Rectified Linear Unit (ReLU) is exploited to transform the input into feature-based sequence, which is expressed as: = ′ (15) where and are weight matrix and bias vector derived from model training. Similarly, and can be obtained in the same manner. The features of three-direction vibration are thoroughly distilled.  Let x i , y i , and z i represent the ith set of axial, radial, and tangential vibrating signals, respectively. [x i , y i , z i ] ∈ R n can be the inputs of the proposed fault diagnosis model. Each input vector is sent to a fully connected layer for feature extraction. Specifically, the active function Rectified Linear Unit (ReLU) is exploited to transform the input into feature-based sequence, which is expressed as: where w and b are weight matrix and bias vector derived from model training. Similarly, Y i and can be Z i obtained in the same manner. The features of three-direction vibration are thoroughly distilled.
Then, the Bi-GRUs are taken to extract the internal properties of the input data. We will now briefly describe the theory of GRU. For details, see, for example, [17,19]. GRU is a variation of the Recurrent Neural Network (RNN), which is widely applied to learn sequential characteristics from a sequential data stream [33]. Concretely, a GRU is established on the gating mechanism to modulate the existing memory and new memory: a reset gate determines the combination of new input and previous information, and an update gate selects previous information to preserve [34]. The sigmoid activation function computes both gate outputs, while tanh is used as the candidate hidden activation function. Distinctively, the GRU controls the information flow and exposes the full hidden information without a memory unit [35]. In this way, the GRU is generally employed to obtain the hidden states of inputs.
In our model, a Bi-GRU related to sequence data is established. The input is fed into both a forward GRU and a backward GRU. These two GRUs of reverse timing are connected to one output. For the input X i , we have that: where h X i is the Bi-GRU output, consisting of forward hidden state → h X i and backward hidden state ← h X i . Likewise, for inputs Y i and Z i , we can also compute h Y i and h Z i . Furthermore, the Highway Network is introduced to deal with the outputs of Bi-GRU (i.e., h X i , h Y i , and h Z i ). In a Highway Network, a part of the inputs can be pass through the layer straightly, while others are revised using gating units [36].
Supposing that T is the transform gate and C is the carry gate, for an input vector h X , the output of Highway Network is: The computation of X H i is facilitated by setting C = 1 − T. Thereby, Equation (17) is modified to: where H is an affine transform followed by a non-linear activation function, and W H , W T , and W C are the parametric matrices optimized during training. Correspondingly, the outcomes Y H i and Z H i are computed as well. With the application of Highway Network, more significant properties from vibrating signals are retained whilst irrelevant information properties are removed. The Highway Network outputs are concatenated and sent to the two-layer activation. The information propagation during this process is delivered as: then Lastly, fault identification is performed via a softmax classifier to obtain the probability distribution over the different working states, which is: where w o is a trainable matrix, b o is the offset vector, and C is the number of working states.

Model Training
The training of our model is conducted by using the cross entropy and L 2 regularization as the loss function: where g i is the real working state distribution, and p i is the predicted one. The parameter λ r is the weight of L 2 regularization. The gradients and other parameters are updated through back propagation with the learning rate λ l :

Experimental Setting
Fault classification experiments are conducted on data measured from coaxial-rotor units of real working conditions. The dataset contains three-directional vibration signals of both normal-state and bearing-circlip failure, with nine distinctive rotating speeds. The data are randomly split into training, validating, and testing according to the ratio 8:1:1.
Our model is finetuned using an ADAM optimizer with β 1 = 0.9, β 2 = 0.999, and ε = 1 × 10 −8 . All the parameter matrices are generated within the distribution of U(−0.1, 0.1) randomly, while the bias is initialized as 0. The initialized parameters take 1000 times iteration for parameter optimization. Moreover, the hidden states dimension of the GRU cell is set as 256 with a learning rate of 0.001. The batch size and the training epoch are set to 10. The dropout rate is 0.5 to prevent overfitting. With respect to model training, the L 2 -regularization weight is set as 0.0001.

Experimental Protocol
In this experiment, the confusion matrix is taken as the evaluating metric to verify the working performance of our model. In a confusion matrix (Table 2), the basic parameters, accuracy, precision, and recall can be computed [37,38]. On the one hand, the index accuracy denotes the correctlyclassified samples of all the samples, which is: On the other hand, precision and recall also provide a richer measure of classification performance. According to Equation (27), precision stands for the proportion of predicted positive cases that are correctly real positives: In contrast, recall refers to the proportion of real positives that are correctly predicted as positive: In this experiment, 'precision' can be taken as a key metric for working performance demonstration. Normal working condition indicates 'positive', and failure indicates 'negative'; it is more severe to classify a fault as a normal state than the opposite. Therefore, the accuracy is adopted as the primary evaluation index and precision as the auxiliary.

Main Results
The main purpose of the experiment is to verify the effectiveness of the proposed model for the aspect-level sentiment analysis. The following baseline methods are taken for comparison. KNN (K-Nearest Neighbor): A simple machine learning algorithm for both linear and non-linear classification, which is insensitive to outliers. SVM (Support Vector Machine): A supervised learning algorithm for generalized linear classification, which can prevent the issue of over fitting.
CNN (Convolutional Neural Network): A classical deep learning method, consisting of a convolutional layer, a pooling layer, and a fully connected layer for feature extraction.
LSTM (Long Short-Term Memory): A variant of RNN, which can detect the hidden states of the inputs.
GRU: As mentioned above, GRU is a variant of RNN that captures the hidden states using a gating mechanism and an activation function.
Bi-GRU: Bi-GRU is built based on a forward GRU and a backward GRU, to encode and decode the inputs, respectively.
Notably, five parameters are picked up as features and applied to KNN and SVM, which are the average value, standard deviation, maximum value, kurtosis, and IMF1 marginal spectral energy.
The fault diagnosis results of different rotation speeds are obtained. The test accuracy of all the working conditions is given in Table 3. The confusion matrix that presents test outcomes is shown in Figure 10. Notably, the values on the diagonal line refer to all the correctly recognized numbers of samples, corresponding to the accuracy of the test. The darker the color is, the higher the accuracy is.
In nine evaluation settings, our model obtains the best and most consistent performance in classification accuracy. As presented in Table 3, an average accuracy as high as 99.4% is obtained due to its distinctiveness in both high-speed and low-speed rotating fault diagnosis. Among all the approaches, SVM and KNN have relatively lower accuracy than the deep learning models because current machine learning models still have challenges in feature extracting and exploiting. A deep neural network is taken as a better approach than machine learning algorithms in capturing the working properties under different working conditions. Specifically, one can easily see that the accuracy of low-speed diagnosis is even lower, which validates the analysis results based on Hilbert marginal spectrum energy.  Figure 10. The confusion matrix for the fault classification.

Concluding Remarks
In this work, a deep-learning based fault diagnosis model is devised and deplo on the task of coaxial-rotor unit working state identification in a turbojet engine. To c ture the essential characteristics of different states, the test rig is developed to reprod the working condition and regenerate the fault data of the coaxial-rotor unit. Based on test results, attention is paid to the faults on the bearing circlip, which are typically d aged at first. The three-direction vibration signals are measured and dedicatedly analy to characterize the working status. The HHT analysis is performed to study the fai mechanism of the bearing circlips. In line with the property of bearing circlip failur fault classification method, based on the integration of Bi-GRU and the Highway N work, is proposed. Basically, the GRU cells are capable of capturing the hidden in mation from inputs, while the Highway Network removes the unrelated features. fault diagnosis performance of the proposed model is evaluated on nine different work conditions. An average accuracy of 99.4% of the proposed model is obtained, which performs the state-of-the-art methods. Specifically, our model also shows its superio in low-speed fault classification. The experimental results indicate that our model obt sufficient classification accuracy. That is, our model shows its distinctive superiorit

Concluding Remarks
In this work, a deep-learning based fault diagnosis model is devised and deployed on the task of coaxial-rotor unit working state identification in a turbojet engine. To capture the essential characteristics of different states, the test rig is developed to reproduce the working condition and regenerate the fault data of the coaxial-rotor unit. Based on the test results, attention is paid to the faults on the bearing circlip, which are typically damaged at first. The three-direction vibration signals are measured and dedicatedly analyzed to characterize the working status. The HHT analysis is performed to study the failure mechanism of the bearing circlips. In line with the property of bearing circlip failure, a fault classification method, based on the integration of Bi-GRU and the Highway Network, is proposed. Basically, the GRU cells are capable of capturing the hidden information from inputs, while the Highway Network removes the unrelated features. The fault diagnosis performance of the proposed model is evaluated on nine different working conditions. An average accuracy of 99.4% of the proposed model is obtained, which outperforms the state-of-the-art methods. Specifically, our model also shows its superiority in low-speed fault classification. The experimental results indicate that our model obtains sufficient classification accuracy. That is, our model shows its distinctive superiority in dealing with the detected sensing signals precisely and thus can be applied to coaxial-rotor unit fault diagnosis tasks.