Feature Analysis of Degradation Signal of Rolling Stock Lateral Damper

: The lateral damper is one of the key components of rolling stock. Establishing the relationship between the degraded signal and the health state of the lateral damper is important in order to perform timely performance detection and fault diagnosis. This paper proposes a wavelet packet cross-correlation method (WPCC) that is based on wavelet packet transform (WPT) and cross-correlation analysis (CCA). First, the vibration signals under di ﬀ erent running speeds, di ﬀ erent running conditions, and di ﬀ erent track excitations were collected and analyzed. Second, the wavelet packet transform was used to select larger energy band signals for reconstruction. Subsequently, the WPCC coe ﬃ cient was calculated between the reference signal and the signal to be measured. The proposed method was applied to analysis of vibration signals of the lateral damper performance degradation. The lateral damper health condition was divided into four intervals, and the average accuracy calculated under di ﬀ erent running speeds, di ﬀ erent running conditions, and di ﬀ erent track excitation was 95%.


Introduction
The research and development of rolling stock technology has made the running speed of trains continue to increase. At the same time, the requirements for the stability and comfort of the train are becoming more stringent. The lateral damper is one of the key components that affect the running performance of the train. The vibration can be effectively suppressed in the lateral direction of the train, and to ensure riding performance [1]. In the course of use, the lateral damper will undergo a process of failure from normal to degradation, due to various internal and external factors. Different running states of rolling stock will show different data characteristics. The current maintenance method of lateral dampers is regular maintenance. As the running speed increases, the future development trend will inevitably be to carry out state maintenance to identify the health status of components. The establishment of the relationship between different states and signal data is what is missing from regular to state repair. The key to research is feature extraction and state recognition, and the problem is the identification of early degradation states. There is an urgent need to explore new feature extraction and feature selection methods to more effectively realize the health state recognition of lateral dampers. The faults of high-speed train systems and key components are often expressed in the form of various signals. Using appropriate feature analysis methods, the working status of components can be directly obtained from the signals to achieve the purpose of status recognition and fault diagnosis. Early signal analysis research focused on classic signal analysis methods. Traditional signal processing methods, such as spectrum analysis, Fourier analysis, envelope analysis, correlation analysis, and maximum entropy spectrum analysis, played a huge role in system and component state recognition. Li Tan [2] adopted principal component analysis (PCA) and multi-fractal analysis method in order to analyze the vibration characteristics of trains during single failure and cross-mixing failure of key components of high-speed trains. Yu Y. [3] adopted empirical mode decomposition (EMD) and adaptive short-time Fourier transform (STFT) diagnosis method to diagnose the faults of railway freight rolling bearings. The results showed that the method could effectively obtain fault information in the time-frequency domain. Zeng Hua et al. [4] used time-domain statistical methods to calculate the root mean square of the bogie vibration signal and analyzed the frequency response characteristics of the vibration transfer function between each position. Guo-Liang S [5] proposed a state estimation method of wavelet packet energy moment that was effective for accurately estimating different faults of high-speed train. When compared with other time domain and frequency domain analysis, the wavelet packet analysis method can effectively eliminate the interference of noise in the fault signal, and it has better sensitivity and stability to features. Therefore, this paper used wavelet packet analysis to extract the features of the fault signal.
Presently, a variety of state recognition methods have been developed in the field of fault state recognition, mainly divided into methods that are based on qualitative empirical knowledge, methods based on analytical models and methods based on data-driven. Lin S. [6] used principal component analysis (PCA) and rough set-based fault identification method for train door fault identification. The comparison found that, although the PCA method can determine whether a fault has occurred, it cannot perform fault identification. Li Yong [7] combined Hilbert transform and Berg spectrum estimation method to extract the time-frequency and spectrum characteristics of the fault current signal, and used the binary classification support vector machines as a high-impedance ground fault classifier for the traction network, which can effectively identify the fault type of the traction network, and obtain a high classification accuracy rate. Since the degree of performance degradation is gradual, it cannot be divided into several fixed fault types, so it is unreasonable to use a classifier, such as a support vector machine, to evaluate the degree of performance degradation. References [8,9] applied cross correlation analysis to transformer internal fault identification and bearing fault diagnosis, and achieved certain results. When compared with classification methods, such as support vector machines, cross-correlation analysis can effectively extract minor differences between signals, so it is suitable for early degraded signal feature analysis.
Most of the existing researches on the fault diagnosis of train parts are to identify the failure status of the components. There is no in-depth analysis on the identification of the early degradation of the components. In this paper, a high-speed train lateral damper was used as an example, and a feature analysis method of wavelet packet cross-correlation analysis was proposed in order to identify the early degradation state. It provided a reference for the identification of the degradation status of other components, and introduced new ideas for the early warning calibration of key parts of high-speed trains.

Lateral Vibration Analysis
The body will produce severe lateral vibration during the operation of high-speed trains, due to the uneven track, the working state of train components, bends, turnouts, and other factors. Body vibrations are often complex and they contain multiple frequency components, multiple directions, changing over time. The frequency of vibration will have various effects on the car body and passengers, especially the occurrence of resonance, which will lead to physical discomfort in all aspects of passengers and even cause harm to health [10]. Therefore, it is necessary to evaluate the vibration acceleration of the train in order to process the vibration design of the high-speed train. Riding comfort is one of the important items that show the quality of running vehicles. Taking the lateral damper as the research object, the lateral sperling indices of the vehicle body are analyzed when the damping value changes from 0 to 100 kN.s.m −1 , as shown in Table 1. Figure 1 shows the relation curve of the sperling index versus the damping coefficient.   It can be seen from Figure 1 and Table 1: when the damping value is small, such as less than 20 kN.s.m −1 , the sperling index of cab decreases as the damping value increases; when the value is larger, the sperling index increases as the damping value increases. At the same speed, the sperling index of the front and rear cab is similar with the variation law of damping value, but the amplitudes are slightly different. It is shown that the lateral damper is a key component of the high-speed train bogie that affects the running performance of the train. It can effectively suppress the vibration of the train in the lateral direction and ensure the riding quality of the train and the comfort of the ride. When the lateral damper is out of order, it will produce abnormal vibration. From practice, it is known that the frequency of the vehicle body is mainly concentrated in the low frequency band. In order to extract the subtle characteristics of the train under different faults, the wavelet packet decomposition has great advantages, and the wavelet packet cross-correlation analysis method is proposed to reflect the health condition of the lateral damper that corresponds to different vibration signals.

Wavelet Packet Transformation
The wavelet packet transformation (WPT) method is developed on the basis of the multiresolution wavelet transform (WT), which is, the signal is decomposed in a series of different levels of space. It retains the characteristics of the wavelet transform, with the increase of resolution, the widened spectral window further fines the fine quality, and can perform comprehensive timefrequency decomposition of the signal, which has unique advantages for vibration signal analysis. As an extension of the WT, WPT can be achieved on any frequency of infinite refinement [11], and Figure 2 shows its frequency domain distribution. It can be seen from Figure 1 and Table 1: when the damping value is small, such as less than 20 kN.s.m −1 , the sperling index of cab decreases as the damping value increases; when the value is larger, the sperling index increases as the damping value increases. At the same speed, the sperling index of the front and rear cab is similar with the variation law of damping value, but the amplitudes are slightly different. It is shown that the lateral damper is a key component of the high-speed train bogie that affects the running performance of the train. It can effectively suppress the vibration of the train in the lateral direction and ensure the riding quality of the train and the comfort of the ride. When the lateral damper is out of order, it will produce abnormal vibration. From practice, it is known that the frequency of the vehicle body is mainly concentrated in the low frequency band. In order to extract the subtle characteristics of the train under different faults, the wavelet packet decomposition has great advantages, and the wavelet packet cross-correlation analysis method is proposed to reflect the health condition of the lateral damper that corresponds to different vibration signals.

Wavelet Packet Transformation
The wavelet packet transformation (WPT) method is developed on the basis of the multi-resolution wavelet transform (WT), which is, the signal is decomposed in a series of different levels of space. It retains the characteristics of the wavelet transform, with the increase of resolution, the widened spectral window further fines the fine quality, and can perform comprehensive time-frequency decomposition of the signal, which has unique advantages for vibration signal analysis. As an extension of the WT, WPT can be achieved on any frequency of infinite refinement [11], and Figure 2 shows its frequency domain distribution.
The signal can be decomposed into the interested frequency bands, which is advantageous for extracting the characteristic signal, as can be seen from Figure 2. The mathematical description of the wavelet packet transform is as follows: The bandwidth of analysis frequency in Figure 1 is written as [10]: where B is the bandwidth of analysis frequency, f h is the maximum analysis frequency of the signal, and i is the number of layers decomposed by the wavelet packet, i = 0, 1, 2, 3, . . . Additionally, the upper and lower limits of frequency band be written as [10]: where m is the lower limit of frequency band, n is the upper limit, and j is the index of the sub-band. The signal can be decomposed into the interested frequency bands, which is advantageous for extracting the characteristic signal, as can be seen from Figure 2. The mathematical description of the wavelet packet transform is as follows: The bandwidth of analysis frequency in Figure 1 is written as [10]: where B is the bandwidth of analysis frequency, is the maximum analysis frequency of the signal, and i is the number of layers decomposed by the wavelet packet, i = 0, 1, 2, 3, … Additionally, the upper and lower limits of frequency band be written as [10]: where m is the lower limit of frequency band, n is the upper limit, and j is the index of the sub-band.
The traditional wavelet transform method does not consider the energy distribution on the time axis of each decomposition frequency band, so that the extracted characteristic parameters cannot accurately reflect the characteristics of the fault features. Signals in different states will exhibit different frequency band energy distribution characteristics. The wavelet packet energy moment is introduced in order to calculate the energy of each frequency band. By comparing the sensitivity of each frequency band to component degradation, the frequency band that best reflects the signal characteristics is extracted. The energy moment parameter is the energy moment of the j-th frequency band in the i-th layer after the wavelet packet decomposition of the signal. The calculation formula of is [12]: The traditional wavelet transform method does not consider the energy distribution on the time axis of each decomposition frequency band, so that the extracted characteristic parameters cannot accurately reflect the characteristics of the fault features. Signals in different states will exhibit different frequency band energy distribution characteristics. The wavelet packet energy moment M ij is introduced in order to calculate the energy of each frequency band. By comparing the sensitivity of each frequency band to component degradation, the frequency band that best reflects the signal characteristics is extracted. The energy moment parameter M ij is the energy moment of the j-th frequency band in the i-th layer after the wavelet packet decomposition of the signal. The calculation formula of M ij is [12]: where ∆t is the sampling time interval, n is the total number of sampling points, and k is the sampling point.

Wavelet Packet Cross-Correlation Analysis Method
The cross-correlation analysis method is one of the common methods for time-domain analysis of mechanical vibration signals and an important means of fault diagnosis [13]. The correlation refers to the linear relationship or correlation between variables.
The cross-correlation function embodies the correlation between the two signals, so the cross-correlation analysis method can be used in order to determine the correlation of the two signals in frequency. Assuming that both signals x(t) and y(t) contain the same frequency component and the same period, then the cross-correlation function R xy also contains the frequency component of the same period. Therefore, for the signal feature changes that are caused by component degradation, this feature of the cross-correlation function can be used to effectively extract features and identify states.
The cross-correlation function of discrete vibration signals is [14]: where x i and y i+k are the data sequence of the response signal, K is the length of the data sequence, and k is the time delay.
Taking advantage of the WPT in the sensitive frequency band, combined with the correlation function to judge the signal correlation. This method is applicable to evaluate the health condition for rolling stock lateral damper. The specific steps are as follows: (1) Wavelet packet transform. Five-layer wavelet packet transform transformed the collected lateral vibration acceleration signals of the vehicle body. (2) Calculate cross-correlation coefficient. The cross-correlation functions of each frequency band of the reference signal and the signal to be tested are calculated, respectively, and the cross-correlation coefficients of eight frequency bands are calculated according to Equation (5).
The correlation coefficient ρ xy gives an indication on the strength of the linear relationship between the two time series x(t) and y(t) [15]. As the correlation coefficient approaches zero, there is less of a relationship between the two time series (closer to uncorrelated). The closer the coefficient to either −1 or 1, the stronger is the correlation between the two time series. Figure 3 shows the process. The cross-correlation function embodies the correlation between the two signals, so the crosscorrelation analysis method can be used in order to determine the correlation of the two signals in frequency. Assuming that both signals x(t) and y(t) contain the same frequency component and the same period, then the cross-correlation function also contains the frequency component of the same period. Therefore, for the signal feature changes that are caused by component degradation, this feature of the cross-correlation function can be used to effectively extract features and identify states.
The cross-correlation function of discrete vibration signals is [14]: where i x and k i y  are the data sequence of the response signal, K is the length of the data sequence, and k is the time delay.
Taking advantage of the WPT in the sensitive frequency band, combined with the correlation function to judge the signal correlation. This method is applicable to evaluate the health condition for rolling stock lateral damper. The specific steps are as follows: (1) Wavelet packet transform. Five-layer wavelet packet transform transformed the collected lateral vibration acceleration signals of the vehicle body. (2) Calculate cross-correlation coefficient. The cross-correlation functions of each frequency band of the reference signal and the signal to be tested are calculated, respectively, and the crosscorrelation coefficients of eight frequency bands are calculated according to Equation (5).
The correlation coefficient gives an indication on the strength of the linear relationship between the two time series x(t) and y(t) [15]. As the correlation coefficient approaches zero, there is less of a relationship between the two time series (closer to uncorrelated). The closer the coefficient to either −1 or 1, the stronger is the correlation between the two time series. Figure 3 shows the process.

Results and Discussion
Based on vehicle track coupling dynamics, with Multi-Body Simulation (MBS) software SIMPACK, the nonlinear dynamics model of the China standard EMU vehicle system is established, which fully considered the geometric nonlinearity of wheel-rail contact, the creep nonlinearity of wheel-rail, and the nonlinear suspension. LMA type treads and 60 kg/m steel rail are adopted, and the distance between backs of wheel flanges is in accordance with the Chinese standard of 1353 mm, as shown in Figure 4. Based on vehicle track coupling dynamics, with Multi-Body Simulation (MBS) software SIMPACK, the nonlinear dynamics model of the China standard EMU vehicle system is established, which fully considered the geometric nonlinearity of wheel-rail contact, the creep nonlinearity of wheel-rail, and the nonlinear suspension. LMA type treads and 60 kg/m steel rail are adopted, and the distance between backs of wheel flanges is in accordance with the Chinese standard of 1353 mm, as shown in Figure 4. The experimental study is carried out in order to identify the degradation state of the lateral damper under different running speeds, different running conditions and different track excitations. The rolling stock was run for 3 min. with the performance of the lateral damper from normal to gradually degraded (deterioration degree of the normal 90% to 10%), recording vehicle body lateral vibration acceleration sensor data at a sampling frequency of 256 Hz. Six seconds (1536 sample points) were taken from the experimental data as one sample, and 270 samples at each speed were obtained.

Selection of Sensitive Signal Bands
It is necessary to select a sensitive signal band that reflects the degradation trend of the lateral damper performance in order to accurately establish the relationship between the signal and the degradation state of the lateral damper. Considering the analysis of signals within 30 Hz, only the wavelet packet coefficients of the first eight bands in the fifth layer are reconstructed. The corresponding frequency range is: 0~3.75 Hz, 3.75~7.5 Hz, 7.5~11.25 Hz, 11.25~15 Hz, 15~18.75 Hz, 18.75~22.5 Hz, 22.5~26.25 Hz, and 26.25~30 Hz. Taking the running speed of 200 km/h under the linear running condition on the Beijing-Tianjin track excitation as an example, the signals of 90~10% degradation degree are decomposed by wavelet packet, and the wavelet packet energy moment Mij are shown in Table 2. Figure 5 shows the variation curves of the wavelet packet cross-correlation method (WPCC) coefficient of lateral acceleration signals of vehicle body in each degenerate degree.  The experimental study is carried out in order to identify the degradation state of the lateral damper under different running speeds, different running conditions and different track excitations. The rolling stock was run for 3 min. with the performance of the lateral damper from normal to gradually degraded (deterioration degree of the normal 90% to 10%), recording vehicle body lateral vibration acceleration sensor data at a sampling frequency of 256 Hz. Six seconds (1536 sample points) were taken from the experimental data as one sample, and 270 samples at each speed were obtained.

Selection of Sensitive Signal Bands
It is necessary to select a sensitive signal band that reflects the degradation trend of the lateral damper performance in order to accurately establish the relationship between the signal and the degradation state of the lateral damper. Considering the analysis of signals within 30 Hz, only the wavelet packet coefficients of the first eight bands in the fifth layer are reconstructed. The corresponding frequency range is: 0~3.75 Hz, 3.75~7.5 Hz, 7.5~11.25 Hz, 11.25~15 Hz, 15~18.75 Hz, 18.75~22.5 Hz, 22.5~26.25 Hz, and 26.25~30 Hz. Taking the running speed of 200 km/h under the linear running condition on the Beijing-Tianjin track excitation as an example, the signals of 90~10% degradation degree are decomposed by wavelet packet, and the wavelet packet energy moment M ij are shown in Table 2. Figure 5 shows the variation curves of the wavelet packet cross-correlation method (WPCC) coefficient of lateral acceleration signals of vehicle body in each degenerate degree.
The higher energy components are M 50 and M 53 , as can be seen from Table 2. It can be seen from Figure 3 that different frequency bands have different sensitivity to the degradation of lateral dampers, and the M 53 frequency band reflects the most obvious change trend of its performance degradation. Therefore, the M 53 frequency band (11.25~15 Hz) is selected as the sensitive frequency band for calculating WPCC coefficients.  The higher energy components are M50 and M53, as can be seen from Table 2. It can be seen from Figure 3 that different frequency bands have different sensitivity to the degradation of lateral dampers, and the M53 frequency band reflects the most obvious change trend of its performance degradation. Therefore, the M53 frequency band (11.25~15 Hz) is selected as the sensitive frequency band for calculating WPCC coefficients.

Solution of Wavelet Packet Cross-Correlation Coefficients
Normal performance signal is selected as the reference signal, and each degradation degree is taken as the signal to be tested. The wavelet packet cross-correlation coefficients of different degradation degrees is calculated according to the steps in 3.2. The results under different running speeds, different running conditions, and different track excitation are shown in Tables 3-5, respectively. It can be seen from the above table that the number of the coefficient under different degrees of degradation is different. As the degree of degradation deepens, the value of the coefficient also

Solution of Wavelet Packet Cross-Correlation Coefficients
Normal performance signal is selected as the reference signal, and each degradation degree is taken as the signal to be tested. The wavelet packet cross-correlation coefficients of different degradation degrees is calculated according to the steps in Section 3.2. The results under different running speeds, different running conditions, and different track excitation are shown in Tables 3-5, respectively.  It can be seen from the above table that the number of the coefficient under different degrees of degradation is different. As the degree of degradation deepens, the value of the coefficient also becomes smaller, which is mutually consistent with the physical meaning of the WPCC coefficients that are defined above. The physical meaning of the number of relationships is consistent. This indicates that the WPCC coefficients can well reflect different degrees of degradation.
It can be seen from Tables 4 and 5 that the WPCC coefficient of degeneration degree of the lateral damper is significantly reduced under the curve condition and the track excitation of German railway spectra of high irregularity (GRSHI), indicating that the similarity of the signals is weakened more. Wavelet packet cross-correlation coefficients can better reflect different degradation states in harsh environments. Therefore, the wavelet packet cross-correlation coefficient is a reliable reference base.

Curve Fitting of Wavelet Packet Cross-Correlation Coefficients
It can be obtained from the above simulation analysis that the wavelet packet cross-correlation analysis method can effectively identify the health status of the lateral damper. The wavelet packet cross-correlation coefficient can reflect the degradation state of the lateral damper and it is a reliable reference. The above simulation data are fitted with the least squares fitting polyfit function in MATLAB, and the following exponential curve can be obtained, as shown in Figure 6.
Information 2020, 11, x FOR PEER REVIEW 8 of 11 becomes smaller, which is mutually consistent with the physical meaning of the WPCC coefficients that are defined above. The physical meaning of the number of relationships is consistent. This indicates that the WPCC coefficients can well reflect different degrees of degradation. It can be seen from Tables 4 and 5 that the WPCC coefficient of degeneration degree of the lateral damper is significantly reduced under the curve condition and the track excitation of German railway spectra of high irregularity (GRSHI), indicating that the similarity of the signals is weakened more. Wavelet packet cross-correlation coefficients can better reflect different degradation states in harsh environments. Therefore, the wavelet packet cross-correlation coefficient is a reliable reference base.

Curve Fitting of Wavelet Packet Cross-Correlation Coefficients
It can be obtained from the above simulation analysis that the wavelet packet cross-correlation analysis method can effectively identify the health status of the lateral damper. The wavelet packet cross-correlation coefficient can reflect the degradation state of the lateral damper and it is a reliable reference. The above simulation data are fitted with the least squares fitting polyfit function in MATLAB, and the following exponential curve can be obtained, as shown in Figure 6.  Table 6 shows the results of goodness of fit. Goodness of fit is the degree to which the regression model fits the observations. It can be seen from the above table that the RMSE is 0.0154 and the Rsquare is 0.9916, indicating that the regression curve fits the observations very well.

Evaluation of the Health Condition for Lateral Damper
The assessment of the health conditions requires prior knowledge as a criterion. The gradient between 100% and 70% is denoted as stage 1 (Health), 70%~40% as stage 2 (Sub-health), 40% to 10% as stage 3 (Fault), and 10% as stage 4 (Scrap), according to the practical experience of engineering projects. Three specific stages of performance degradation degree, namely 70%, 40%, and 10%, are predetermined.
The WPCC coefficients of these three stages are extracted as the critical value, which could constitute four degradation degree intervals. When evaluating the signal to be tested, the WPCC coefficients is first calculated, and the nearest two critical values are found, and the degradation degree of the signal is divided into the degradation degree interval that is formed by the two critical  Table 6 shows the results of goodness of fit. Goodness of fit is the degree to which the regression model fits the observations. It can be seen from the above table that the RMSE is 0.0154 and the R-square is 0.9916, indicating that the regression curve fits the observations very well.

Evaluation of the Health Condition for Lateral Damper
The assessment of the health conditions requires prior knowledge as a criterion. The gradient between 100% and 70% is denoted as stage 1 (Health), 70%~40% as stage 2 (Sub-health), 40% to 10% as stage 3 (Fault), and 10% as stage 4 (Scrap), according to the practical experience of engineering projects. Three specific stages of performance degradation degree, namely 70%, 40%, and 10%, are predetermined.
The WPCC coefficients of these three stages are extracted as the critical value, which could constitute four degradation degree intervals. When evaluating the signal to be tested, the WPCC coefficients is first calculated, and the nearest two critical values are found, and the degradation degree of the signal is divided into the degradation degree interval that is formed by the two critical values. Six degradation degrees (90%, 80%, 60%, 50%, 30%, and 20%) were selected for the signals to be tested, with 10 samples per degradation degree, at the speed of 200 km/h, 250 km/h, 300 km/h, and 350 km/h, respectively. Figure 7 shows the test results of the signals that were collected by the vehicle body at the speed of 200 km/h under the linear running condition on the Beijing-Tianjin track excitation.
values. Six degradation degrees (90%, 80%, 60%, 50%, 30%, and 20%) were selected for the signals to be tested, with 10 samples per degradation degree, at the speed of 200 km/h, 250 km/h, 300 km/h, and 350 km/h, respectively. Figure 7 shows the test results of the signals that were collected by the vehicle body at the speed of 200 km/h under the linear running condition on the Beijing-Tianjin track excitation. It can be seen from Table 7 that the algorithm has a high recognition rate at all speeds. As the train running speed increases, the correct recognition rate for the health state of the lateral damper increases, that is, as the speed increases, the characteristics of wavelet packet cross-correlation coefficients are more prominent in the representation of health status. It can be concluded from Table 8 that the wavelet packet cross-correlation coefficient feature can obtain a higher recognition rate than the traditional wavelet packet energy and wavelet entropy feature, especially at high speed. The recognition rate of wavelet packet cross-correlation coefficient feature reaches a satisfactory result of more than 96% when the speed is 350 km/h.  Tables 7, 9, and 10 shows the accuracy of the lateral acceleration signal analysis collected by the sensor under different running speeds, different running conditions, and different track excitation. It can be seen from the table that there is a high recognition rate under each condition, and the overall It can be seen from Table 7 that the algorithm has a high recognition rate at all speeds. As the train running speed increases, the correct recognition rate for the health state of the lateral damper increases, that is, as the speed increases, the characteristics of wavelet packet cross-correlation coefficients are more prominent in the representation of health status. It can be concluded from Table 8 that the wavelet packet cross-correlation coefficient feature can obtain a higher recognition rate than the traditional wavelet packet energy and wavelet entropy feature, especially at high speed. The recognition rate of wavelet packet cross-correlation coefficient feature reaches a satisfactory result of more than 96% when the speed is 350 km/h.  Tables 7, 9, and 10 shows the accuracy of the lateral acceleration signal analysis collected by the sensor under different running speeds, different running conditions, and different track excitation. It can be seen from the table that there is a high recognition rate under each condition, and the overall average accuracy rate is 95%, which verifies the effectiveness and feasibility of the wavelet packet cross-correlation coefficient in evaluating the health condition of the lateral damper. Table 9. Accuracy rates at each speed.

Accuracy Rates Front Cab Rear Cab
Curve 96.67% 96.11% Table 10. Accuracy rates at each speed.

Conclusions
This paper realizes the quantitative description of the health state of the lateral damper by solving the wavelet packet cross-correlation coefficient in order to accurately and deeply evaluate of health condition for rolling stock lateral damper. Using the simulation data of the performance degradation of the lateral damper under different running speeds, different running conditions, and different track excitation, a satisfactory evaluation accuracy is obtained, and the applicability and effectiveness of the method are verified.

Conflicts of Interest:
The authors declare no conflict of interest.