Next Article in Journal
Seismic Exploration Wireless Sensor System Based on Wi-Fi and LTE
Next Article in Special Issue
Probability-Based Algorithm for Bearing Diagnosis with Untrained Spall Sizes
Previous Article in Journal
A Threshold-Based Max-log-MPA Low Complexity Multiuser Detection Algorithm
Previous Article in Special Issue
A Hybrid Sensor Fault Diagnosis for Maintenance in Railway Traction Drives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

State-Degradation-Oriented Fault Diagnosis for High-Speed Train Running Gears System

1
School of Computer Science and Engineering, Changchun University of Technology, Changchun 130012, China
2
CRRC Changchun Railway Vehicles Co., Ltd., Changchun 130062, China
3
Department of Automation, Tsinghua University, Beijing 100084, China
4
Academy of Astronautics, Harbin Institute of Technology, Harbin 150001, China
5
School of Mechatronic Engineering, Changchun University of Technology, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(4), 1017; https://doi.org/10.3390/s20041017
Submission received: 5 January 2020 / Revised: 31 January 2020 / Accepted: 11 February 2020 / Published: 13 February 2020
(This article belongs to the Special Issue Sensors for Fault Diagnosis)

Abstract

:
As one of the critical components of high-speed trains, the running gears system directly affects the operation performance of the train. This paper proposes a state-degradation-oriented method for fault diagnosis of an actual running gears system based on the Wiener state degradation process and multi-sensor filtering. First of all, for the given measurements of the high-speed train, this paper considers the information acquisition and transfer characteristics of composite sensors, which establish a distributed topology for axle box bearing. Secondly, a distributed filtering is built based on the bilinear system model, and the gain parameters of the filter are designed to minimize the mean square error. For a better presentation of the degradation characteristics in actual operation, this paper constructs an improved nonlinear model. Finally, threshold is determined based on the Chebyshev’s inequality for a reliable fault diagnosis. Open datasets of rotating machinery bearings and the real measurements are utilized in the case studies to demonstrate the effectiveness of the proposed method. Results obtained in this paper are consistent with the actual situation, which validate the proposed methods.

Graphical Abstract

1. Introduction

Nowadays, railway construction and operation are high-speed, which considerably affects the sustainability and the rapid progress of the national economy. The safety and reliability of railway operation have drawn more and more attention. As a core system that directly affects the smooth running of high-speed trains, the running gears system is a key component of the train in providing dynamics and traction performance and has the functions of buffering, vibration isolation, generating power and supporting the vehicle body. As a matter of fact, the running gears system is a complex system that is coupled with many components. Among those components, the corrosion, shedding and degradation of one component may easily spread the local fault and propagate into a major fault in the system level, which may cause unexpected losses. It is very important to analyze system performance degradation [1,2] and prevent fault propagation. Therefore, the timely diagnosis of a running gears system plays a key role in ensuring the safe operation of trains.
In the past few decades, many work has been done on fault diagnosis technologies in different aspects. Traditionally, the fault diagnosis methods for high-speed train systems mainly consist of two factors: The model-based one and the signal-based one [3]. The signal-based fault diagnosis methods analyze the measurement data for fault diagnosis without knowing/identifying the system model. The common signal-based methods include multivariate statistical analysis [4,5], machine learning methods [6,7], signal processing methods [8], information fusion methods [9,10], etc. The above methods have high requirements on data volume and type. The measurements of the running gears system of high-speed trains are of high sampling frequency and long measurement cycle. Once the running gears system of the high-speed train is damaged, it will be repaired in a timely manner and regularly maintained. However, the adequate failure data for analysis are very rare. The most important step of the signal-based fault diagnosis method is the accurate extraction of fault characteristics [11]. If the measurements with inconspicuous fault characteristics are used for the learning, accurate feature extraction normally cannot be guaranteed, and the result of fault diagnosis will be directly affected. In contrast, model-based fault diagnosis methods, which are more sensitive to uncertainty/fault, have lower requirements for data quantity and data type. They are more suitable for a reliable fault diagnosis of high-speed train running gears systems.
Among the model-based fault diagnosis methods, a class of improved methods based on the Kalman filter is used to deal with the influence of system/measurement noise. Yan et al. [12] introduced the extended Kalman filter model and the improved support vector machine model into the online classification framework, where the nonlinearity is approximated by the first-order Taylor expansion and the online fault detection of the ventilation and air-conditioning system is later achieved. Ning et al. [13] proposed an enhanced robust Kalman filtering method based on the optimal radial basis function neural network, which reduces the impact of state errors and observation errors on state estimation. Wang et al. [14] reduced the influence of estimated noise of on-board fault diagnosis and fault-tolerant control through the Kalman filter and improved the robustness of fault diagnosis. Finally, the method is verified by an onboard prototype control system. Akai et al. [15] used the particle filter method to realize the state estimation of a target robot in order to realize system filtering under the influence of noise. Wu et al. [16] proposed a fault diagnosis algorithm for early failures based on the generalized estimation. Yin et al. [17] applied a state resonance algorithm-based stochastic resonance method to high-speed train traction systems to achieve optimal extraction and frequency recovery. Cadini et al. [18] used a particle filter to estimate the length of cracks to diagnose the early faults of aeronautical structure cracks. For state estimation with unknown inputs, some researchers use unknown input residual generators and estimation filters to achieve fault diagnosis [19]. Scholars also consider the decomposition of external disturbances and combine unknown input observers to complete composite fault diagnosis [20]. For nonlinear systems, He et al. [21] applied rank particle filtering to diagnosis the fault of underwater vehicles. In recent years, to improve the accuracy of state estimation and reduce the false alarm rate, scholars have carefully studied the measurement noise and unknown disturbance. Zhou et al. [22] applied an improved odorless Kalman filter and a radial basis function based on an adaptive noise factor to the pumping unit for fault diagnosis. Liu et al. [23] used an adaptive controller to process unknown system parameters and disturbances to achieve target state tracking. Gan et al. [24] used an adaptive factor to improve the accuracy of the unscented Kalman filter algorithm, thereby improving the accuracy of fault diagnosis of analog circuits. However, most of the existing methods are based on state estimation and improvement of a single node of the monitoring system. Such methods mostly rely on the reliability of information collection nodes. For an actual high-speed train information collection system, the measurements obtained in different locations of the same monitoring area contain different information, which will indirectly affect the consistency of information collection and reduce the accuracy of fault diagnosis. With the development of a multi-sensor system, distributed estimation has been studied and applied in different scenarios. Jing et al. [25] realized multi-sensor data fusion of the planetary gearbox by fusing vibration signal, sound signal, drive motor current, shaft speed, oil chip and other measurements, and then detected various system faults. Liu et al. [26] proposed a multi-sensor intelligent data fusion method based on the ant colony algorithm to realize fault detection of the gear box, aiming at the unreliability and instability of fault detection results of single sensor vibration signals. To avoid the disadvantage of combining a single information source with a decision method, a fault diagnosis method combining different types of sensor data sources and multiple classifier decisions is proposed in [27]. Banerjee et al. [28] proposed and investigated a hybrid method for fault signal classification based on sensor data fusion by using the support vector machine and short-term Fourier transform techniques. This method estimates the state of the target node by combining the information of surrounding nodes. This method has been successfully applied in individual formation [29], stability monitoring [30] and output coordination [31,32]. It can be concluded that distributed monitoring has better robustness and stability for single sensor monitoring in the sensor information fusion architecture. Considering multiple sensor information of the target monitoring system can improve the accuracy of state estimation more than considering single sensor information.
While the model-based method has been applied in fault diagnosis for high-speed train, there are still some problems. First, for the area monitored by multiple sensors, the estimation does not fully consider the influence of neighbor nodes which will further affect the later state estimation and fault diagnosis. Secondly, the statistical characteristics and the noise in the measurements collected in high-speed trains are normally unknown. The impacts of initial state and noise will be continuously amplified during the recursive estimation process, which may reduce the accuracy of fault diagnosis. Finally, the degradation of the system performance cannot be avoided. If the degradation is ignored, it will increase the false alarm rate of fault diagnosis and reduce the diagnostic accuracy. Aiming at solving the above problems, this paper proposes a distributed state estimation filter, which is combined with system state degradation characteristics to achieve fault diagnosis. This method is different from fault prognosis, the purpose of which is to analyze the possible faults in the future through the fault state feature [33], and then to provide guidance for future system health management. The purpose of this paper is to diagnose the existing faults of the degraded system. From the perspective of methods, the signal-based fault prognosis method [34] also has the constraint problem of data quantity and data type. The model-based fault prognosis method [35] mainly predicts the future fault generation time by analyzing the evolution trend of the fault. The method proposed in this paper is to diagnose the existing faults through the state evolution trend. They are similar, but still different.
Innovations and contributions of this article are as follows.
(1)
Comprehensively consider the measuring point position and information acquisition method of a composite sensor. A distributed topology structure is established by taking the axle box bearing of an actual running gears system as an example. Based on this structure, a bilinear distributed filter is proposed, and the gain parameters of the filter are designed to minimize the mean square error.
(2)
Unbiased constraint conditions are used to reduce the impact of the initial unknown information of nodes on state estimation. By constructing the difference to deal with the problem of colored noise in real measurements, estimation accuracy is improved.
(3)
A nonlinear degradation model of the Wiener process considering temperature change characteristics is built to describe the state degradation phenomenon during train operation. The solution of nonlinear degradation process parameters is given by maximum likelihood estimation and combined with distributed filters to increase the accuracy of fault diagnosis.
Finally, to demonstrate the effectiveness of the proposed method, this paper conducts simulation experiments on the open dataset of rotating machinery bearings. The fault diagnosis results given are the same as actual experimental results. The proposed method was applied to operational monitoring data of a certain type of high-speed train running gears system, and fault diagnosis results obtained were consistent with a real situation.
This paper is arranged as follows. The second part describes the problem and introduces the preparation work. The third part proposes a method based on the Wiener process and multi-sensor filtering. The fourth part uses two case studies to test and analyze the above method. The final part concludes the article.

2. Preliminaries and Problem Formulations

The running gears system is a complex mechanism, and its state information is an important indicator of traffic safety, as shown in Figure 1a,b. At present, the train fault diagnosis system is mainly composed of the train main engine, vehicle extension, preprocessing and composite sensors. The monitoring area of high-speed train operation consists of three or four sensors. Each sensor is integrated with two or more sensor units to detect different physical quantities and collect different types of data, such as temperature and vibration. Therefore, a single sensor is also called a composite sensor component. However, information acquisition in fault diagnosis based on single sensor information is not comprehensive or accurate. For the bearings of high-speed trains running a gears system, equipment status information from a single sensor is often significantly incomplete, which cannot provide an accurate and detailed data basis for subsequent fault diagnosis. Therefore, the fault diagnosis method based on single sensor information widely used in trains is limited in the accuracy of diagnosis and the completeness of information. To solve the problem, this paper abstracts a high-speed train running gears system into a multi-sensor system. The train control and management system processes the node information and gives full play to the advantages of a multi-node joint operation so that each node in a multi-sensor system can generate a consistent interpretation and description of the monitoring target. The reliability of sensor information is enhanced to achieve the purpose of collaborative monitoring and improve the accuracy of state estimation.
Eleven composite sensors are installed in the support area of the monitoring bearing, and the installation direction needs to be consistent with the direction of the impact signal. The specific position is shown as follows, A1–A4: Measuring points of the axle box bearing; B1–B3: Measuring points for traction motor bearings; C1–C4: Measuring points for gearbox bearings. Taking the axle box bearing and its composite sensor in the above system as an example, a networked sensor system as shown in Figure 2 is established. The topology directed graph is represented by G = ( V , E , D ) , where the vertex set is represented as V = { 1 , 2 , 3 , 4 } , the set of edges is represented as ε ν × ν .
This paper considers the effects of a sensor’s gain attenuation, colored measurement noise and accuracy reduction noise. The model of the target system is as follows.
x ( k + 1 ) = A ( k ) x ( k ) + B ( k ) u ( k ) + w ( k ) y i ( k ) = λ i C i ( k ) x ( k ) + D i ( k ) u ( k ) + v i ( k ) + s i ( k )
To better conform to the changes in the actual state, the running characteristics of axle box bearings of the running section need to be analyzed. The axle box bearings of the running section are made of metal, and the specific heat capacity of metal is a function of temperature. In general, its specific heat capacity increases with increasing temperature. Therefore, the change of temperature of the axle box bearing will no longer conform to the law of linear change, and a bilinear term will be generated in a system equation. Consider describing the above system as the following bilinear form.
x ( k + 1 ) = A ( k ) x ( k ) + B ( k ) u ( k ) + N ( k ) u ( k ) x ( k ) + w ( k ) y i ( k ) = λ i C i ( k ) x ( k ) + D i ( k ) u ( k ) + v i ( k ) + s i ( k )
where k is the discrete time index, u ( k ) is the control input signal, w ( k ) is the process noise and the mean is 0, the covariance is Q. v ( k ) represents measurement noise of the i-th sensor. s ( k ) is precision-degraded noise of the i-th sensor, the mean value is 0 and the covariance is O. λ is gain reduction of the i-th sensor. A ( k ) , B ( k ) , C ( k ) and D ( k ) are known coefficient matrices, and N ( k ) u ( k ) x ( k ) is a bilinear term.
In addition, as a key mechanical rotating part of the running gears system, the performance of the bearing plays an essential role in the safe running of the train. Generally speaking, mechanical equipment must go through a series of regular degradation stages from normal operation to complete failure. Therefore, taking the degradation process into account in the system model can more accurately describe the changing process of the actual system and increase the accuracy of state estimation. This paper considers the modeling of system state degradation based on the Wiener process. The more common nonlinear degradation model is as follows.
X ( t ) = X ( 0 ) + 0 t μ ( u ; θ ) d u + σ B ( t )
where μ ( u ; θ ) is a non-linear function used to characterize the non-linear characteristics of the system degradation. θ is vector of unknown parameters contained in the function. σ B ( t ) is the diffusion term used to describe the uncertainty of the degradation process. When the surface temperature of the friction pair rises to a certain extent for a short time, a secondary quenching layer and a high-temperature tempering layer will be produced on the surface of the bearing part. This burnt layer will cause obvious changes in the structure and performance of the surface of the bearing part, affecting the performance of the friction pair. This paper considers improving the above-mentioned Wiener process degradation model to achieve state degradation modeling.

3. Main Results and Discussion

The fault diagnosis method proposed in this paper is mainly divided into the following three stages: First, the distributed filtering design of multi-sensor systems. Second, model of the system state degradation based on the Wiener process. Third, fault diagnosis methods for the running gears system, shown in Figure 3.

3.1. Multi-Sensor Filter

Measurement noise during actual sensor operation is colored noise, the measurement noise in the differential form is constructed as follows.
v i ( k ) = ψ i ( k 1 ) v i ( k 1 ) + ζ i ( k 1 )
where ζ represents Gaussian white noise. To process v, an auxiliary signal z similar to the above is constructed.
z i ( k ) = y i ( k ) ψ i ( k 1 ) v i ( k 1 ) = λ i C i ( k ) x ( k ) + D i ( k ) u ( k ) + v i ( k ) + s i ( k ) ψ i ( k 1 ) v i ( k 1 ) = λ i C i ( k ) x ( k ) + D i ( k ) u ( k ) + ζ i ( k 1 ) + s i ( k )
The auxiliary signal z thus obtained will no longer contain colored noise terms, which is more convenient for processing. Continuing to use the equation of state of the original system, the measurement equation is described by auxiliary signal z, and the following equivalent equation of target system can be obtained:
x ( k + 1 ) = A ( k ) x ( k ) + B ( k ) u ( k ) + N ( k ) u ( k ) x ( k ) + w ( k ) z i ( k ) = λ i C i x ( k ) + D i ( k ) u ( k ) + ζ i ( k 1 ) + s i ( k )
Considering the mutual influence between multi-sensor nodes, a distributed state estimation structure for bilinear systems is proposed as follows.
x ^ i ( k + 1 ) = A ( k ) x ^ i ( k ) + L i ( k ) u ( k ) + N ( k ) u ( k ) x ^ i ( k ) + j N i H ( k ) [ z j ( k ) λ i C i ( k ) x ( k ) ]
where x ^ i ( k + 1 ) represent the value of the estimated state of the i-th sensor, N i is the set consisting of node i and its neighbors and H ( k ) is the filter gain of the multi-sensor.
In the process of state estimation, the uncertainty of the initial state will be transferred in the process of estimation iteration, leading to the deviation of state estimation. To satisfy the unbiasedness of the filter, the state estimation process is processed as follows, so that the mean of each state estimate is the same as the mean of the states, which is E [ x ^ ( k ) ] = E [ x ( k ) ] . The above formula is used as an unbiased constraint, so for the state estimation equation and the state equation at time k + 1 :
E [ x ^ i ( k + 1 ) ] = A ( k ) E [ x ^ i ( k ) ] + L i ( k ) u ( k ) + N ( k ) u ( k ) E [ x ^ i ( k ) ] + j N i H ( k ) D i ( k ) u ( k ) E [ x i ( k + 1 ) ] = A ( k ) E [ x i ( k ) ] + B i ( k ) u ( k ) + N ( k ) u ( k ) E [ x i ( k ) ] E [ x ^ i ( k + 1 ) ] E [ x i ( k + 1 ) ] = [ L i ( k ) B i ( k ) + j N i H ( k ) D i ( k ) ] u ( k ) = 0
Therefore, the unbiased constraint can be equivalent to:
L i ( k ) B i ( k ) + j N i H ( k ) D i ( k ) = 0
Next, the gain parameter H of the filter is obtained by minimizing the trace of estimated error covariance. The estimated error function is established as follows.
e i ( k + 1 ) = [ x ^ i ( k + 1 ) x i ( k + 1 ) ]
The gain parameter H of the filter is expressed as:
H B = arg min H E { e T ( k + 1 ) · e ( k + 1 ) } = arg min H E { t r a c e [ e ( k + 1 ) · e T ( k + 1 ) ] } = arg min H E { t r a c e [ ( j N i H ( k ) ζ j ( k 1 ) + H ( k ) S j ( k ) w j ( k ) ) · ( ) T ] } = arg min H t r a c e { j N i [ H ( k ) R j ( k 1 ) H T ( k ) + H ( k ) O j ( k ) H T ( k ) Q j ( k ) ] }
Using the unbiased constraint as the condition for minimizing the above formula, the Lagrange function is constructed.
J s ( H , λ ) = [ j N i H ( k ) R j ( k 1 ) H T ( k ) + H ( k ) O j ( k ) H T ( k ) Q i ( k ) ] + λ [ L i ( k ) B i ( k ) + j N i H ( k ) D j ( k ) ]
Then, H can be expressed as:
H ( k ) = j N i 1 2 λ [ R j ( k 1 ) + O j ( k ) ] 1 · D j ( k ) j N i H ( k ) D j ( k ) = B i ( k ) L i ( k )
In summary, the analytical solution of the gain parameter H of the distributed filter is finally obtained as:
H ( k ) = [ B i ( k ) L i ( k ) ] j N i { [ R j ( k 1 ) + O j ( k ) ] 1 D j ( k ) } 1 D j 1 ( k ) [ R j ( k 1 ) + O j ( k ) ] 1 D j ( k )
Remark 1.
This paper only considers the measurement noise of a system as colored noise, instead of colored process noise. The reason is that the state degradation process is the main reason that directly affects system change. The characteristics of process noise have little effect on the system. If process noise is processed in the same form, it will increase the computational complexity.

3.2. Parameter Estimation of State Degradation

Considering the effects of random shocks and environmental factors (temperature), this paper builds a nonlinear degradation process based on the Wiener degradation model and provides a solution for parameter estimation.
Construct a nonlinear state degradation model.
X ( t ) = X ( 0 ) + 0 t μ ( u ; θ ) d u + i = 1 N J i + u ( t ) x s t a t e ( t ) + σ B ( t )
where X ( t ) is the degradation process driven by the standard drift Brownian motion B ( t ) , μ ( u ; θ ) is a non-linear function, J is the amplitude of random shocks, N represent the random shock, and u ( t ) x ( t ) is the environmental factor on the system state degraded parameters. It can be seen that when μ ( u ; θ ) = μ , the degradation process becomes a linear degradation model. In order to describe the above models in detail, this paper gives a parameter estimation method for a class of nonlinear models.
μ ( t ; θ ) = a b t b 1 , X ( t ) = X ( 0 ) + a t b + i = 1 N J i + u ( t ) x s t a t e ( t ) + σ B ( t )
The unknown parameters of the drift and diffusion terms in the degradation process are estimated, assuming current time is t k , and the historical detection data of equipment degradation is x 1 , x 2 , , x k , where x is obtained by sampling the degradation process at equal intervals. X ˜ 1 : k represents the incremental degradation process.
Theorem 1.
The average value of increment X ˜ 1 : k of the above degradation process is:
a E [ t 1 b t 0 b t 2 b t 1 b t i + 1 b t i b ] T E [ u 1 x s t a t e 1 u 0 x s t a t e 0 u 2 x s t a t e 2 u 1 x s t a t e 1 u i + 1 x s t a t e i + 1 u i + 1 x s t a t e i + 1 ] T
The variance is σ B 2 · Q and follows a normal distribution.
Proof of Theorem 1. 
Tectonic degradation increment:
X ( i + 1 ) = x ( i + 1 ) x ( i ) = a [ t ( i + 1 ) ] b + u [ t ( i + 1 ) ] x s t a t e [ t ( i + 1 ) ] + σ B B [ t ( i + 1 ) ] a t ( i ) b u [ t ( i ) ] x s t a t e [ t ( i ) ] σ B B [ t ( i ) ] = a [ t ( i + 1 ) b t ( i ) b ] + Δ + σ B { B [ t ( i + 1 ) ] B [ t ( i ) ] }
where Δ can represent the influence of the environment (temperature) on state degradation, and the mean value of the above equation can be obtained:
E [ X ( i + 1 ) ] = a E [ t ( i + 1 ) b t ( i ) b ] E [ Δ ]
The covariance can be expressed as:
var [ X ( i + 1 ) · X ( j + 1 ) ] = var { x [ t ( i + 1 ) ] x [ t ( i ) ] } · { x [ t ( j + 1 ) ] x [ t ( j ) ] } = E { x [ t ( i + 1 ) ] x [ t ( i ) ] } { x [ t ( j + 1 ) ] x [ t ( j ) ] } E { x [ t ( i + 1 ) ] x [ t ( i ) ] } · E { x [ t ( j + 1 ) ] x [ t ( j ) ] } + σ B ( · ) = E { a 2 { [ t ( i + 1 ) ] b t ( i ) b } [ t ( j + 1 ) b t ( j ) b ] } + E { a [ t ( i + 1 ) b t ( i ) b ] · Δ ( j ) } + E { Δ ( i ) · a [ t ( j + 1 ) b t ( j ) b ] } + E [ Δ ( i ) Δ ( j ) ] E { a 2 [ t ( i + 1 ) b t ( i ) b ] [ t ( j + 1 ) b t ( j ) b ] } a E [ t ( i + 1 ) b t ( i ) b ] Δ ( i ) Δ ( j ) · a E [ t ( j + 1 ) b t ( j ) b ] + σ B ( · ) = E { σ B { B [ t ( i + 1 ) ] B [ t ( i ) ] } { a [ t ( j + 1 ) b t ( j ) b ] + Δ ( j ) } } + E { σ B { B [ t ( j + 1 ) ] B [ t ( j ) ] } { a [ t ( i + 1 ) b t ( i ) b ] + Δ ( i ) } } + E { σ B 2 { B [ t ( i + 1 ) ] B [ t ( i ) ] } · { B [ t ( j + 1 ) ] B [ t ( j ) ] } } = σ B 2 E { { B [ t ( i + 1 ) ] B [ t ( i ) ] } · { B [ t ( j + 1 ) ] B [ t ( j ) ] } }
where σ B ( · ) represents the term associated with σ B .
Set Q = E { { B [ t ( i + 1 ) ] B [ t ( i ) ] } { B [ t ( j + 1 ) ] B [ t ( j ) ] } } , the variance of the degenerate increment is σ B 2 · Q .
Construct the likelihood function:
l ( θ ; σ B | X ˜ 0 : k ) = 1 2 π · σ B Q · e ( X ˜ 1 : k a μ 1 : k + Δ 1 : k ) 2 2 σ B 2 Q
according to the above likelihood function, the desired a and σ B can be obtained:
d ln ( l ) d a = X ˜ 1 : k · μ 1 : k + Δ 1 : k · μ 1 : k 2 σ B 2 Q μ 1 : k 2 · a σ B 2 Q = 0 ; a ^ = X ˜ 1 : k + Δ 1 : k 2 μ 1 : k
d ln ( l ) d σ B = 1 σ B + ( X ˜ 1 : k a μ 1 : k + Δ 1 : k ) 2 σ B 3 Q = 0 ; σ ^ B = X ˜ 1 : k a μ 1 : k + Δ 1 : k Q
The estimated values of a and σ B are substituted into the log-likelihood function, and the log-likelihood function is maximized by the simplex method to get b. □
Remark 2.
The state degradation model needs to be combined with the state estimation model. The continuous function of the degradation process needs to be discretized.

3.3. Residual Generator Design

Discretizing the state variables in stage 2 and combining them with the state variables in stage 1:
x ( k ) = x ( k ) + X ( k ) x ^ ( k ) = x ^ ( k ) + X ( k )
Constructing a residual generator:
T i ( k ) = y i ( k ) y ^ i ( k )
where y ^ i ( k ) = λ i C i ( k ) x ^ ( k ) + D i ( k ) u ( k ) , and then:
T i ( k ) = λ i C i ( k ) [ x ( k ) x ^ ( k ) ] + v i ( k ) + s i ( k ) = λ i C i ( k ) · E + ( k 1 ) E ( k 1 ) · [ x ( k ) x ^ ( k ) ] + v i ( k ) + s i ( k ) = λ i C i ( k ) · E + ( k 1 ) { A ( k 1 ) x ( k 1 ) + N ( k 1 ) u ( k 1 ) x ( k 1 ) + B ( k 1 ) u ( k 1 ) + w ( k 1 ) [ A ( k 1 ) x ^ i ( k 1 ) + L i ( k 1 ) u ( k 1 ) + N ( k 1 ) u ( k 1 ) x ^ i ( k 1 ) + j N i H ( k 1 ) [ Z j ( k 1 ) λ j C j ( k 1 ) x ( k 1 ) ] ] } + v i ( k ) + s i ( k ) = λ i C i ( k ) · E + ( k 1 ) { A ( k 1 ) [ x i ( k 1 ) x ^ i ( k 1 ) ] + N ( k 1 ) u ( k 1 ) [ x i ( k 1 ) x ^ i ( k 1 ) ] + w ( k 1 ) j N i H ( k 1 ) [ ζ j ( k 2 ) + s j ( k 1 ) ] } + v i ( k ) + s i ( k )
Theorem 2.
The given residual generator satisfies the following distribution
T i ( k ) { E [ v i ( k ) ] , P ˜ ( k ) }
where the mean value and variance of v are derived by following equation:
v i ( k ) = y i ( k ) λ i C i ( k ) x ^ i ( k ) D i ( k ) u ( k ) s i ( k )
Proof of Theorem 2. 
E [ T i ( k ) ] = E [ v i ( k ) ]
P ˜ ( k ) = E [ T i ( k ) · T i T ( k ) ] = [ λ i C i ( k ) ] · E + ( k 1 ) { Q ( k 1 ) j N i H ( k 1 ) [ R ( k 2 ) + O ( k 1 ) ] H T ( k 1 ) } · { [ λ i C i ( k ) ] · E + ( k 1 ) } T + E [ v i ( k ) · v i T ( k ) ] + O ( k )
The fault diagnosis threshold is obtained by the Chebyshev inequality:
Pr ( μ ( X ) α σ ( X ) X μ ( X ) + α σ ( X ) ) 1 1 α 2

4. Practical Verifications and Discussion

4.1. State Estimation

Firstly, this paper presents the state estimation curve of four nodes in monitoring systems, as shown in Figure 4.
Figure 5a,b, respectively, shows the state estimation error value and state estimation error covariance of a single node. It can be seen that the state estimation error keeps within ± 0.5 , while the state estimation error covariance converges to a smaller value, proving the effectiveness of the filter.

4.2. Cincinnati Dataset

First, the 2003 Cincinnati test 1 bearing 2 channel 4 data was used for simulation experiments. The sampling time of this dataset is from 12:06:24 on 22 October 2003 to 23:39:56 on 25 November 2003. The test bench device contains four bearings, and the AC motor is coupled to the shaft through a friction band. The speed was kept at 2000 RPM. A 6000-pound radial load was applied to the shaft and bearings. Each dataset consists of a single file of 1 second. The sampling rate is 20 kHz, and 20,480 points are saved at a time.
The degree of bearing damage in the experiment is related to operating time of device, and the corresponding data format will also change. Generally speaking, it can be clearly seen from the amplitude of the data whether the monitoring system has failed. Therefore, to better illustrate the reliability of algorithm, this article selects some data from above data set, including non-faulty data points, early fault data points, and severe fault data points. In this paper, the initial wear state of the bearing is defined as early failure, but it does not affect the continued operation of the whole system. A serious fault may refer to the state of near failure of the bearing in the period before the end of the experiment. At this point, the serious fault has affected the operation of system and even caused some damage to bearing itself. Therefore, it is of greater practical significance to the diagnosis of serious fault, and diagnosis of early fault has certain guiding significance to prediction and maintenance of the future.To better demonstrate effectiveness of the proposed method, data with small amplitude differences were selected from above three types of data for verification, as shown in Figure 6.
Among them, time 0–200 is a non-fault data point; time 200–400 is an early fault data point; 400–600 is a severe fault data point. After passing the filter with state degradation, the innovation vector is obtained as shown in Figure 7a. Taking the adjustable parameter α = 7 in the Chebyshev inequality, the threshold value obtained are ± 0.53 . For the convenience of viewing, a decision function as shown in Figure 7b is given. The threshold exceeded is set to 1, and the threshold not exceeded is set to 0. It can be seen that the fault occurred from 200–400 moments, but because the fault occurred early, the number is relatively sparse, and the fault is a serious fault at 400–600 moments, and the number of faults is relatively large.

4.3. High-Speed Train Temperature Dataset

This section uses an actual running gears system to prove the effectiveness of proposed algorithm. Selecting data of the failed train and the data of non-faulted train. The data is derived from temperature value of temperature sensor on the right position of the 2 axis of one-run part of the train, in Figure 8.
Time 0–300 is a non-faulty data point, and time 300–600 is a faulty data point. After passing through the filter with state degradation, the innovation vector is obtained as shown in Figure 9a. Taking the adjustable parameter α = 7 of Chebyshev inequality, the obtained threshold is ± 0.56 . Similarly, for the convenience of viewing, the decision function shown in Figure 9b is given, and it can be seen that there is a fault in time 300–600, which is consistent with the error information given by actual high-speed train.

4.4. Performance Comparison

In this paper, the proposed algorithm is compared with algorithms commonly used in fault diagnosis. Partial data of the Cincinnati public dataset is selected for fault diagnosis. The diagnostic accuracy is shown in following table (see Table 1). Among them, FNR is the rate of false negatives and FPR represents the rate of false positives.

5. Conclusions

The stable operation of high-speed trains is closely related to the running gears system. This paper proposes a method based on the Wiener state degradation process and multi-sensor filtering for running gears system diagnosis. First, a distributed topology diagram of running gears system axle box bearings and its composite sensors is established by analyzing a running system of a high-speed train. Secondly, a distributed filter is proposed to the problem of unknown initial information of sensor nodes by the unbiased constraint condition, considering the problems of gain attenuation and precision decline, and the gain parameters of the filter are calculated under the premise of minimum error mean square error. Then, to improve the fault diagnosis accuracy of the running gears system, this paper built a nonlinear degradation model of the Wiener process considering the factor of temperature, and estimated parameters through maximum likelihood estimation. Finally, a fault diagnosis threshold is obtained by using the Chebyshev inequality. The performance of the method is verified by an open dataset of rotating mechanical bearings and temperature data of a certain type of high-speed rail running gears system. The result is as expected.
The method proposed in this paper is more meaningful for practical engineering, and improves the accuracy of fault diagnosis for systems that cannot obtain available data. This method aims to diagnose early faults. For other types of dynamic systems, thresholds can be set based on their actual characteristics and the definition of the corresponding fault level. Future research work starts from the correlation between states, decouples filter states, and considers a non-linear degradation problem of non-Markov processes, making the model closer to the actual process.

Author Contributions

C.C. proposed the idea and designed the experiment. W.W. conducted experiments and completed the manuscript. H.L. provides guidance for data analysis and thesis writing. B.Z. made recommendations and supervised the work. G.C. and W.T. provide the experimental environment/equipment. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grants 61903047, 61973046 and 61803044, Jilin Province Development and Reform Commission under Grant 2019C040-3.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ji, H.; He, X.; Shang, J.; Zhou, D. Incipient fault detection with smoothing techniques in statistical process monitoring. Control Eng. Pract. 2017, 62, 11–21. [Google Scholar] [CrossRef]
  2. Chen, H.; Jiang, B.; Lu, N. A multi-mode incipient sensor fault detection and diagnosis method for electrical traction systems. Int. J. Control 2018, 16, 1783–1793. [Google Scholar] [CrossRef]
  3. Su, L.; Ma, L.; Qin, N.; Huang, D.; Andrew, H.K. Fault diagnosis of high-speed rrain bogie by residual-squeeze net. IEEE Trans. Ind. Inform. 2019, 15, 3856–3863. [Google Scholar] [CrossRef]
  4. Huang, J.; Yan, X. Quality-driven principal component analysis combined with kernel least squares for multivariate statistical process monitoring. IEEE Trans. Control Syst. Technol. 2019, 27, 2688–2695. [Google Scholar] [CrossRef]
  5. Chen, H.; Jiang, B.; Ding, S.X.; Lu, N.; Chen, W. Probability-relevant incipient fault detection and diagnosis methodology with applications to electric drive systems. IEEE Trans. Control Syst. Technol. 2018, 27, 2766–2773. [Google Scholar] [CrossRef]
  6. Wang, Z.; Jia, L.; Kou, L.; Qin, Y. Spectral kurtosis entropy and weighted SaE-ELM for bogie fault diagnosis under variable conditions. Sensors 2018, 18, 1705. [Google Scholar] [CrossRef] [Green Version]
  7. Fu, J.; Chu, J.; Guo, P.; Chen, Z. Condition monitoring of wind turbine gearbox bearing based on deep learning model. IEEE Access 2019, 7, 57078–57097. [Google Scholar] [CrossRef]
  8. Zhou, R.; Bao, W.; Li, N.; Huang, X.; Yu, D. Mechanical equipment fault diagnosis based on redundant second generation wavelet packet transform. Digit. Signal Process. 2010, 20, 276–288. [Google Scholar] [CrossRef]
  9. Cai, B.; Liu, Y.; Fan, Q.; Zhang, Y.; Liu, Z.; Yu, S.; Ji, R. Multi-source information fusion based fault diagnosis of ground-source heat pump using Bayesian network. Appl. Energy 2014, 114, 1–9. [Google Scholar] [CrossRef]
  10. Irhoumah, M.; Pusca, R.; Lefevre, E.; Lefevre, E.; Mercier, D.; Romary, R.; Cristian, D. Information fusion with belief functions for detection of inter-turn short circuit faults in electrical machines using external flux sensors. IEEE Trans. Ind. Electron. 2017, 65, 2642–2652. [Google Scholar] [CrossRef]
  11. Gao, Z.; Cecati, C.; Ding, S. A survey of fault diagnosis and fault-tolerant techniques-part I: fault diagnosis with model-based and signal-based approaches. IEEE Trans. Ind. Electron. 2015, 62, 3757–3767. [Google Scholar] [CrossRef] [Green Version]
  12. Yan, K.; Ji, Z.; Shen, W. Online fault detection methods for chillers combining extended kalman filter and recursive one-class SVM. Neurocomputing 2016, 228, 205–212. [Google Scholar] [CrossRef] [Green Version]
  13. Ning, Y.; Wang, J.; Han, H.; Tan, X.; Liu, T. An optimal radial basis function neural network enhanced adaptive robust kalman filter for GNSS/INS integrated systems in complex urban areas. Sensors 2018, 18, 3091. [Google Scholar] [CrossRef] [Green Version]
  14. Wang, Y.; Sun, Y.; Chang, C.l.; Hu, Y. Model-based fault detection and fault tolerant control of SCR urea injection system. IEEE Trans. Veh. Technol. 2016, 65, 4645–4654. [Google Scholar] [CrossRef]
  15. Akai, N.; Morales, Y.; Murase, H. Simultaneous pose and reliability estimation using convolutional neural network and Rao–Blackwellized particle filter. Adv. Robot. 2018, 32, 930–944. [Google Scholar] [CrossRef]
  16. Wu, Y.; Jiang, B.; Lu, N. A descriptor system approach for estimation of incipient faults with application to high-speed railway traction devices. IEEE Trans. Syst. Man Cybern. Syst. 2017, 49, 2108–2118. [Google Scholar] [CrossRef]
  17. Yin, J.; Xie, Y.; Chen, Z.; Peng, T. Weak-fault diagnosis using state-transition-algorithm-based adaptive stochastic-resonance method. J. Cent. South Univ. 2019, 26, 1910–1920. [Google Scholar] [CrossRef]
  18. Cadini, F.; Sbarufatti, C.; Corbetta, M.; Giglio, M. A particle filter-based model selection algorithm for fatigue damage identification on aeronautical structures. Struct. Control Health Monit. 2017, 24, 2002. [Google Scholar] [CrossRef]
  19. Jain, T.; Yame, J. Fault-tolerant economic model predictive control for wind turbines. IEEE Trans. Sustain. Energy 2018, 10, 1696–1704. [Google Scholar] [CrossRef]
  20. Bai, W.; Dong, H.; Yao, X.; Ning, B. Robust fault detection for the dynamics of high-speed train with multi-source finite frequency interference. ISA Trans. 2018, 75, 76–87. [Google Scholar] [CrossRef] [PubMed]
  21. He, J.; Li, Y.; Jiang, Y.; Li, Y.; Li, A. Propeller fault diagnosis based on a rank particle filter for autonomous underwater vehicles. Brodogr. Brodogr. Theory Pract. 2018, 69, 147–164. [Google Scholar] [CrossRef]
  22. Zhou, W.; Li, X.; Yi, J.; He, H. A novel UKF-RBF method based on adaptive noise factor for fault diagnosis in pumping unit. IEEE Trans. Ind. Inform. 2018, 15, 1415–1424. [Google Scholar] [CrossRef]
  23. Liu, S.; Jiang, B.; Mao, Z.; Ding, S. Adaptive backstepping based fault-tolerant control for high-speed trains with actuator faults. Int. J. Control Autom. Syst. 2019, 17, 1408–1420. [Google Scholar] [CrossRef]
  24. Gan, X.; Gao, W.; Dai, Z.; Liu, W. Research on WNN soft fault diagnosis for analog circuit based on adaptive UKF algorithm. Appl. Soft Comput. 2016, 50, 252–259. [Google Scholar]
  25. Jing, L.; Wang, T.; Zhao, M.; Wang, P. An adaptive multi-sensor data fusion method based on deep convolutional neural networks for fault diagnosis of planetary gearbox. Sensors 2017, 17, 414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Liu, Z.; Guo, W.; Tang, Z.; Chen, Y. Multi-sensor data fusion using a relevance vector machine based on an ant colony for gearbox fault detection. Sensors 2015, 15, 21857–21875. [Google Scholar] [CrossRef] [Green Version]
  27. Niu, G.; Han, T.; Yang, B.; Tan, A. Multi-agent decision fusion for motor fault diagnosis. Mech. Syst. Signal Process. 2007, 21, 1285–1299. [Google Scholar] [CrossRef] [Green Version]
  28. Banerjee, T.; Das, S. Multi-sensor data fusion using support vector machine for motor fault detection. Inf. Sci. 2012, 217, 96–107. [Google Scholar] [CrossRef]
  29. Viegas, D.; Batista, P.; Oliveira, P.; Silvestre, C.; Chen, C.L.P. Distributed state estimation for linear multi-agent systems with time-varying measurement topology. Automatica 2015, 54, 72–79. [Google Scholar] [CrossRef]
  30. Chen, Z.; Cao, Y.; Ding, S.; Zhang, K.; Koenings, T.; Peng, T.; Yang, C.; Gui, W. A Distributed Canonical Correlation Analysis-based Fault Detection Method for Plant-wide Process Monitoring. IEEE Trans. Ind. Inform. 2019, 15, 2710–2720. [Google Scholar] [CrossRef]
  31. Cai, H.; Lewis, F.L.; Hu, G.; Huang, J. The adaptive distributed observer approach to the cooperative output regulation of linear multi-agent systems. Automatica 2017, 75, 299–305. [Google Scholar] [CrossRef]
  32. Mohammad, P.; Ehsan, H.; Amir, K.; Baris, F.; Bakhtiar, L.; Shih-ken, C.; Shreyas, S. Cooperative vehicle speed fault diagnosis and correction. IEEE Trans. Intell. Transp. Syst. 2018, 20, 783–789. [Google Scholar]
  33. Wu, C.; Zhu, C.; Ge, Y. A new fault diagnosis and prognosis technology for high-power lithium-ion battery. IEEE Trans. Plasma Sci. 2017, 45, 1533–1538. [Google Scholar] [CrossRef]
  34. Wang, X.; McArthur, S.; Strachan, S.; Kirkwood, J.; Paisley, B. A data analytic approach to automatic fault diagnosis and prognosis for distribution automation. IEEE Trans. Smart Grid 2018, 9, 6265–6273. [Google Scholar] [CrossRef] [Green Version]
  35. Yan, W.; Zhang, B.; Dou, W.; Liu, D.; Peng, Y. Low-cost adaptive lebesgue sampling particle filtering approach for real-time li-ion battery diagnosis and prognosis. IEEE Trans. Autom. Sci. Eng. 2017, 14, 1601–1611. [Google Scholar] [CrossRef]
Figure 1. (a) High-speed train test bench; (b) high-speed train running gears system.
Figure 1. (a) High-speed train test bench; (b) high-speed train running gears system.
Sensors 20 01017 g001
Figure 2. Distribution of measuring points of running gears system and topological structure of axle box bearing nodes.
Figure 2. Distribution of measuring points of running gears system and topological structure of axle box bearing nodes.
Sensors 20 01017 g002
Figure 3. Fault diagnosis method framework of high-speed train running gears system.
Figure 3. Fault diagnosis method framework of high-speed train running gears system.
Sensors 20 01017 g003
Figure 4. Distributed state estimation.
Figure 4. Distributed state estimation.
Sensors 20 01017 g004
Figure 5. (a) The estimated error of a node; (b) the estimated error covariance of a node.
Figure 5. (a) The estimated error of a node; (b) the estimated error covariance of a node.
Sensors 20 01017 g005
Figure 6. Cincinnati bearing vibration data values.
Figure 6. Cincinnati bearing vibration data values.
Sensors 20 01017 g006
Figure 7. (a) Innovation vector of Cincinnati vibration data; (b) fault decision function of Cincinnati bearing vibration data.
Figure 7. (a) Innovation vector of Cincinnati vibration data; (b) fault decision function of Cincinnati bearing vibration data.
Sensors 20 01017 g007
Figure 8. Temperature data of running gears system.
Figure 8. Temperature data of running gears system.
Sensors 20 01017 g008
Figure 9. (a) Innovation vector of temperature data; (b) fault decision function of temperature data.
Figure 9. (a) Innovation vector of temperature data; (b) fault decision function of temperature data.
Sensors 20 01017 g009
Table 1. Performance comparison.
Table 1. Performance comparison.
MethodsFNRFPR
State-degraded distributed filter0.26%0.52%
Kalman filter1.28%6.03%
Extended Kalman filter2.04%5.46%

Share and Cite

MDPI and ACS Style

Cheng, C.; Wang, W.; Luo, H.; Zhang, B.; Cheng, G.; Teng, W. State-Degradation-Oriented Fault Diagnosis for High-Speed Train Running Gears System. Sensors 2020, 20, 1017. https://doi.org/10.3390/s20041017

AMA Style

Cheng C, Wang W, Luo H, Zhang B, Cheng G, Teng W. State-Degradation-Oriented Fault Diagnosis for High-Speed Train Running Gears System. Sensors. 2020; 20(4):1017. https://doi.org/10.3390/s20041017

Chicago/Turabian Style

Cheng, Chao, Weijun Wang, Hao Luo, Bangcheng Zhang, Guoli Cheng, and Wanxiu Teng. 2020. "State-Degradation-Oriented Fault Diagnosis for High-Speed Train Running Gears System" Sensors 20, no. 4: 1017. https://doi.org/10.3390/s20041017

APA Style

Cheng, C., Wang, W., Luo, H., Zhang, B., Cheng, G., & Teng, W. (2020). State-Degradation-Oriented Fault Diagnosis for High-Speed Train Running Gears System. Sensors, 20(4), 1017. https://doi.org/10.3390/s20041017

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