Research on the Simulation of Wheelset Response Characteristic Identification of Railway Fastener Loosening

. Rail fastener is a crucial component equipment to ensure the safe operation of the train, and it is very paramount to detect the loose state of the fastener. In this paper, the vertical vibration acceleration signal of wheelset is taken as the research object, and the loose state of fastener is identiﬁed by separating and calculating the key IMF energy entropy. Firstly, based on the ﬁnite element theory and the principle of multibody dynamics, the rigid-ﬂexible coupling simulation model of vehicle track is established. Then, the vertical vibration acceleration signals of the wheelset under the speed of 200km/h are obtained by setting the diﬀerent loosening degrees of the fastener. Finally, we use optimized HHT to process signals, and the orthogonal empirical mode decomposition method (OEMD) is proposed to optimize the orthogonality of the intrinsic mode function, to eliminate the IMF component having poor correlation with the original signal; the Hilbert time spectrum and information entropy theory are combined to calculate the energy entropy of the key IMF, and the HHT energy entropy evaluation algorithm of the vertical acceleration response signal of the train wheelset is proposed. The simulation results show that the HHTenergy entropy of 100% fastener looseness is less than 25%, 50%, and 75%, decreasing trend. The algorithm can recognize the looseness of track fastener through the experiment under diﬀerent working conditions.


Introduction
At present, the methods of detecting fasteners for railway lines in China mainly include manual inspection, track inspection car, and computer-aided visual inspection. When the vehicle is running on the rail, the interaction between wheel and rail will happen, and when the fastener is damaged or missing, the dynamic parameters of the structure will change to some extent [1]. ese changes will be responded to by the vibration signals of some vehicle components in some forms [2,3].
Liu et al. [4] used ANSYS/ LS-DYNA simulation software to establish the vertical coupling vibration model of vehicle-ballastless track-subgrade system, which verified that the sudden change of fastener stiffness had tiny effect on the vertical acceleration of bogie and had obvious effect on the wheelset vibration acceleration; Zhao and Tan [5] established a simulation model of metro vehicle-track flexible body and compared and analyzed the maximum dynamic responses of vehicle-track system under different conditions of fastener failure and different speeds. e research shows that fastener failure has a certain impact on the vibration response of track; Huang [6] analyzed the orthogonality of each intrinsic mode function by numerical simulation and applied the improved HHT method to the damage identification of large structure system; Zhang [7] combined HHT with various theories and applied it to identify rail vibration signals with different fastener looseness under moving load impact. e results show that the method can reflect the change of rail fastener looseness to a certain extent.
To sum up, there is no way to judge the loose state of fastener by separating the vibration acceleration of wheelset [8,9]. is paper presents a method to identify whether the fastener is loose from the response signal of the wheelset before and after the change of the fastener state. It is of great practical significance to identify the fastener.

Vehicle-Track Rigid-Flexible Coupling
Modeling and Dynamic Simulation 2.1. Flexible Track Modeling. In order to obtain the wheelset vibration signal, an accurate vehicle-track rigid-flexible coupling model needs to be established [10]. At first, the track system is considered as a flexible body on the basis of the existing rigid vehicle-track model [11,12]. e slab ballastless track system is shown in Figure 1. e finite element model of rail and rail slab is established in ANSYS. In the multibody dynamics software SIMPACK, force element is used to simulate CA mortar and fastener and assembled into the whole ballastless track system.

Track System Dynamics Model Based on ANSYS.
When the rail is modeled, the rail is regarded as a continuous elastic beam. China's standard 60 kg/section rail parameters are used for modeling and analysis. Because the number of marker points in SIMPACK is limited and if there are too many master DOF (degrees of freedom), it will lead to the incomputable situation in the SIMPACK; therefore, in order to ensure the calculation efficiency, and combined with the purpose of this paper, a 120-meter rail model in ANSYS is established. e distance between two fasteners is 0.6 m, and this length is divided into 20 small units. A total of 201 master DOF nodes are selected. e discrete model of the rail is shown in Figure 2(a). e solid model of the rail slab is shown in Figure 2(b).

Dynamics Analysis on Rail Subsystem.
In this paper, the substructure modal analysis of rail and rail slab needs to reduce the main degree of freedom of the model, and it is reduced by Guyan reduction method in ANSYS. e substructure model of rail and rail slab in this paper is shown in Figure 3. e first 100 modes of the rail and the top 20 modes of the rail slab are calculated, respectively. When selecting the model of rail and rail slab, it is indispensable to accurately reflect the vibration characteristics of rail and rail slab and consider the weight of each model [13]. When the train passes, the vertical vibration is the main vibration of rail, and Figure 4 shows the main mode shapes of the rail and rail slab.
e basic process flow of flexible data transmission between ANSYS and SIMPACK is shown in Figure 5. After modal analysis of rail and rail slab, ANSYS generates geometric model file (.CDB), mass stiffness matrix file (.Sub), and modal mode file (.RST). e interface of FEMBS (finite element multibody system) converts the flexible body data to the standard code (SID) in the format of ASCII readable by SIMPACK, that is, inputting the characteristics of flexible body of rail and rail slab into motion equation to generate SID file. e input information for FEMBS includes DOF and coordinates of nodes, mass attribute, translational and rotational vibration, modal mass matrix, stiffness matrix, damping matrix, and geometric stiffening matrix for initial loads.

Realization of Rigid-Flexible Coupling.
When assembling the track system in SIMPACK, the fastener and the CA mortar are modeled by linear force element. e flexible rails and the rigid wheels are connected with Hertz spring; however, the contact conditions of wheel and rail in SIM-PACK must be rigid wheelset and rigid track and longitudinal relative displacement cannot occur between wheel and rail, and the track model is flexible; the rail does not move with the wheel, so to define a virtual rail between each wheel and flexible rail and define a moving marker point allow the wheel rail force to be transmitted downward to the flexible rail system. e virtual rail body is an object with zero mass and zero-moment of inertia, which only acts as a connection in the dynamic system, and does not affect the other connection structures, but the mass and the moment of inertia of the object in SIMPACK cannot be zero, so the relevant parameters of the virtual rail body are fetched as small as possible, which is assumed to be 1.0e − 6 in this paper, whose impact on wheel rail force transmission is negligible.
In order to ensure that each virtual rail body moves longitudinally along the rail with its corresponding wheelset, the moving marker points on the virtual rail body and the flexible rail are required to define the restraint to hold them, and the wheelset is assembled to the virtual rail body through the hinges, so that the virtual rail body is equivalent to a part of the rail and the vehicle can run along the rail. e wheelrail rigid-flexible coupling model is shown in Figure 6(a). In simulation, the wheel rail force is calculated through Hertz contact between the virtual rail body and rigid wheel. Data exchange between the virtual rail body and the flexible rail is done through the deformation coordination condition and the force balance condition [14,15]. According to the above method, the vehicle and the flexible track are assembled       e established vehicle-track rigid-flexible coupling model is shown in Figure 6(b).

Model Verification When Fasteners Are Not Loose.
Before simulation, firstly we should verify whether the initial state of vehicle-track rigid-flexible coupling model is balanced, which is similar to the verification of rigid body vehicle model. In SIMPACK, we verify the model by calculating the nominal force maximum acceleration value of initial state. e maximum acceleration is calculated to be 6.811 * 10 − 6 m/s 2 , which is less than the evaluation order of magnitude 10 −4 m/s 2 , so it can be determined that the built model is balanced in the initial state and the coupling model is correct.
In [16], the actual measured values of the vehicle vibration acceleration when the high-speed vehicle runs through under different conditions of track irregularity are given, and the simulation calculation is carried out with the VICT simulation software under the same conditions. In this paper, the same test conditions are set in the model to carry out the simulation calculation of vehicle dynamics, and the simulation results are compared with the above results. See Table 1 for comparison. We can see that the simulation results of this paper are not much different from the test results and the VICT simulation results in [17]. e reason for the slight difference is because the vehicle simulation parameters of this paper are different from the other article, and the vertical correctness of the simulation model is verified, which lays a foundation for the subsequent simulation.

SIMPACK Simulation Excitation
Model. SIMPACK can generate excitation by means of inputting the coefficient of the track spectral density formula and the method of converting the measured data. e software can be set up for two types of excitations: one is track-related; the other is railrelated. Due to the fact that there is no unified track spectrum standard in our country, the low-interference track spectrum of German high-speed railway is adopted in this paper. e vertical track irregularities generated in SIMPACK are shown in Figure 7.

Simulation Conditions.
When two fasteners are loosened on the same section of the same track, and two fasteners are loosened continuously, the impact on the vehicle and track system is relatively large [18,19]. erefore, three typical simulation conditions are set up for the looseness of the fastener, as shown in the diagram. e cross line is the fastener loosening position in the picture, and there is one loose fastener on the first condition, there are two loose fasteners on the same section of the track in the second condition, and there are two loose fasteners on the same rail in the third condition.
Selecting the fastener classic stiffness of 50 kN/mm, the vehicle passes through the fastener in different degrees of looseness at 200 km/h speed (not loose, 25% loose, 50% loose, 75% loose, and completely loose), the looseness of the fastener is simulated by changing the stiffness value of the fastener force element in the SIMPACK, the corresponding stiffness value of the fastener loosening degree, and the equivalent stiffness value of the mortar unit length as shown in Table 2, and the looseness degree of fastener is the same in the three working conditions, shown in Figure 8. e vertical acceleration of wheelset and vehicle body is calculated by simulation.

Simulation Results of Wheelset Vertical Acceleration
(1) Wheelset vertical acceleration in condition 1 e vehicle response is the same before the vehicle arrives at the loosened part of the fastener [20,21], so only the vertical acceleration response of the 0.6 s∼1.7 s is given. e vertical acceleration comparison diagram of the wheelset in different degrees of fastener loosening is given, as shown in Figure 9; the acceleration at about 0.72 s∼0.76 s starts changing; with the increase of the degree of fastener looseness, the vibration signal also changes

Identification of Fasteners State in Varying Degrees.
HHT has good reliability in dealing with nonstationary and nonlinear signals, detecting structural faults. However, it is not based on a complete theory and has energy leakage problems. In order to solve this problem, this paper optimizes the algorithm according to the theory of orthogonality and correlation and combines the improved Hilbert Huang algorithm with the theory of information entropy [22,23], which is applied to the identification of rail fastener looseness.  e crux of HHT algorithm is EMD decomposition. e result of EMD decomposition directly affects the accuracy of subsequent signal processing, so it is necessary to ensure that the IMF component of the EMD decomposition should possess completeness and orthogonality, so that no energy leakage in the decomposition can be guaranteed [24,25]. From the practical point of view, Huang et al. consider that there exists orthogonality between all IMFs decomposed by EMD, but there is no rigorous theoretical derivation proving that they are rigorously orthogonal in the overall situation [29,30].
Huang proposed two orthogonality indicators to measure the orthogonality between IMF components, namely, the overall orthogonality index (IOT) and the orthogonality index between two components (IO jk ) [31]. When all the IMF components are rigorously orthogonal to each other, the IOT and IO jk should be zero. (1) In the formula, x i is the ith original signal, imf ji is the jth IMF component of x i , imf k is the kth IMF component of x i , and j is not equal to k. In addition, the energy can be used to measure the degree of orthogonality between IMFs. e energy of the original signal e energy of each IMF is In the formula, imf j (t) is the jth IMF component; if all the IMFs are strictly orthogonal to each other, E x � E j and the leaking energy between the IMFs is zero, that is, where E IOT is the energy of IOT, total energy after decomposition which is constant.
E jk is the leaking energy between imf j and imf k . From the processing steps of EMD decomposition, we can see that, in the actual decomposition process, the mean value curve is obtained by fitting the approximate value, and the decomposed IMF is a part of the original signal, so the different IMF components are approximately orthogonal [29,30]. Huang's theory proves that the EMD orthogonality is in the order of 10 −3 to 10 −2 , and the IMFs are not strictly orthogonal, resulting in the phenomenon of energy leakage and modal aliasing during signal analysis, which will bring errors in the later analysis, making the result inaccurate. In order to eliminate energy leakage, it is necessary to ensure that the orthonormal mode functions are strictly orthogonal to each other [31]. erefore, the traditional EMD algorithm is improved.
In equations (2)-(12), a j � n i�j β i,j , j � 1, 2, . . . , n, β i,j � 1(i � j). From the above calculation, we can see that there is a strict orthogonality in imf j (t), so the linear transformation on imf j (t) will not change its orthogonality, and imf * j (t) is also strictly orthogonal. e integration of the intrinsic mode components after the decomposition of the above EMD is called the orthogonal empirical mode decomposition (OEMD). According to the different orthogonal sequence of IMF components, the OEMD can be divided into three kinds. e way to orthogonalize the high frequency IMF component before orthogonalizing the low frequency IMF component is called OEMD1, the way to deal with IMF in reverse order is called OEMD2, and the way to perform orthogonalization from any IMF is called OEMD3. e essence of OEMD is still EMD; only the IMF components decomposed by EMD are orthogonalized and then reorganized to achieve strict orthogonality of each IMF component.

OEMD Decomposition of the Vertical Acceleration Signals of Wheelset.
e EMD is used to decompose the vertical acceleration signal (Figure 9) of the wheelset without loosening of the fastener under the working condition 1 and obtain 7 IMF components and 1 residual component, as shown in Figure 12.
e same signals are processed by OEMD1 and OEMD2, respectively; the results of 7 IMF components and 1 residual component can be obtained. e overall orthogonal index obtained by the decomposition of EMD is 0.7546, which can be calculated by formula (1). e overall orthogonal indices obtained by OEMD1 and OEMD2 are 0.006463 and 0.007347, respectively. It can be seen that the overall orthogonal index obtained by the OEMD is 2 orders of magnitude higher than that of the overall EMD. ere is little difference in accuracy of the orthogonality of the two orthogonal EMD algorithms. Since the 8th IMF is a residual component with a very small value, it has no orthogonality with the first 7 components. e orthogonality between the 7 IMF components after the decomposition of EMD, OEMD1, and OEMD2 is calculated, respectively, by formula (1); the results of the calculation are shown in Tables 3 and 4. e upper triangular data of Table 3 are the orthogonal indexes between IMF components obtained by EMD decomposition, and the lower triangular data are the orthogonal indexes between IMF components obtained by OEMD1 decomposition. e upper triangular data of Table 4 are the orthogonal indexes obtained by EMD decomposition, and the lower triangular data are the orthogonal indexes obtained by OEMD2 decomposition.
It can be seen from Tables 3 and 4 that the orthogonality index between IMFs by the traditional EMD decomposition is up to 10 −4 , and the IMF orthogonal index by orthogonalization can reach 10 −19 , and the accuracy is improved by 15 orders of magnitude, which is less than 10 −16 , the effective magnitude that the computer calculates in orthogonal index. e result shows that the IMF obtained by using the orthogonalized algorithm has strict orthogonality.
In addition, the energy index can also be used to analyze the orthogonality of IMF. Based on the total energy of EMD, OEMD1, and OEMD2 which can be drawn from formulas (2) and (3), and IMF-decomposed energy shown in Tables 5 and 6, it can be seen that the energy sum of IMF components obtained by EMD has a larger error than that of the original signal, which is 30.8%, and the energy sums of IMF components obtained by the decomposed algorithm of OEMD1 and OEMD2 have the errors of 0.76% and 1.14% against that of the original signal, so the application of OEMD1 and OEMD2 can effectively reduce the energy error, and the effect is obvious.
It is known from Table 5 that the energy of this signal is mainly concentrated in the first 5 orders, accounting for 96.05% of the total energy. erefore, OEMD3 is used to process the signal in the orthogonal order of c5-c4-c3-c2-c1-c6-c7-c8. e decomposed IMFs are shown in Figure 13. e overall orthogonal index is 0.001121 drawn from formula (1)-IOT. e orthogonal index of each component is calculated according to formula (1)-IOT jk , as shown in Table 7; the upper triangular data of this table are the orthogonal indexes between each IMF component obtained by EMD decomposition, and the lower triangular data are the orthogonal indexes by OEMD3 decomposition. From the table, it can be seen that the application of OEMD3 algorithm improves the IMF orthogonal order of magnitude as well as the first two algorithms. is algorithm ensures the strict orthogonality of each IMF component that is decomposed.
According to equations (2) and (3), the energy values of IMFs obtained by the decomposition of OEMD3 algorithm are shown in Table 8. e table shows that the energy error between IMFs decomposed by this algorithm is 0.38%, which is smaller than that of the first two algorithms. e energy errors between the IMF components reduce the energy leakage effectively.
To sum up, from the calculation of the wheelset vertical acceleration signals without fastener loosening, it can be seen that the orthogonal index of IMFs decomposed by OEMD1, OEMD2, and OEMD3 has a great improvement, and the strict orthogonality between the components is guaranteed compared with the orthogonal index of the IMFs decomposed by EMD. According to the results of energy index and orthogonality index calculated by the three algorithms, the total energy index error of each IMF component decomposed by OEMD3 and the error of each component energy index are less than the other two algorithms. erefore, this algorithm is used to decompose the vertical acceleration of the wheelset with different degrees of fasteners loosening in condition 1, condition 2, and condition 3. e wheelset vertical accelerations decomposed by EMD and OEMD3 with fasteners loosening 50% and 100% in condition 1 are shown in Figure 14. It can be seen that, in condition 1, there are some certain differences in the IMF obtained from OEMD3 with the vertical acceleration signals that have different degrees of loose fasteners. In particular, there is one more orthogonal IMF component when the fastener is completely loosened than in the other loosening degrees. is shows that the different degrees of fastener loosening have different effects on the frequency components of the vertical acceleration signal; the difference is most obvious when the fastener is completely loosened, so the time-frequency analysis is considered in the follow-up analysis to identify the loosening features of the fastener.

Orthogonal IMF Selection Based on Correlation.
e above analysis shows that the OEMD can well guarantee the strict orthogonality between the IMF components and ensure that the energy does not leak, but it can be seen from Figure 13 that some interfering signals are introduced by IMFs. ese interfering signals have an impact on the accuracy of the subsequent analysis and should be eliminated. Orthogonalized IMF should have a good correlation with the original signal, so this paper uses the correlation coefficient to select a large correlation IMF. Setting the threshold to 0.5, when the correlation coefficient of IMF and the original signal is greater than 0.5, it is believed that the reliability of the signal is relatively higher. e correlation coefficients between the IMF components and the original signals decomposed by the OEMD3 in the first condition when the fastener is not loosed are shown in Table 9. It can be seen    Method     Method  from the table that only IMF2 and IMF3 satisfy the conditions, and the rest of IMF components and residual components are discarded without any computing.

Fastener Looseness Signal Analysis Based on Hilbert-Huang Transform.
Hilbert transform is applied to the IMF signal processed by OEMD3 to obtain the corresponding Hilbert time spectrum, and the loose feature is extracted by quantitative analysis. Time-frequency spectrum can reflect the distribution of energy in time and frequency and can reflect the change of signal amplitude with time and frequency. However, due to the complexity of wheelset vertical acceleration signal components and Hilbert spectrum, it requires a lot of experience to directly identify the feature information from Hilbert spectrum. In the time and frequency plane, the difference of different degrees of loosening signals is reflected in the different energy distribution in the corresponding area, that is, the uniformity of energy distribution in the same area. erefore, the information entropy theory is applied to the quantitative recognition of Hilbert time spectrum, which simplifies the complexity of the algorithm in the calculation speed and arduousness.

Fastener Looseness Feature Recognition Based on Energy
Entropy. In the information theory, supposing the sample space of discrete random variable X is S � {x 1 , x 1 , . . ., x n } and the probability of random variable X � x j is P j , then the information quantity is e amount of information describes the relationship between the probability of occurrence and information content. e mean value of information quantity I (x j ) in S is the information entropy of X, written as In this paper, the information entropy and the energy distribution of Hilbert time-frequency spectrum are combined to calculate the entropy of fastener looseness in each condition. is method is called HHT energy entropy. e specific method is as follows.
e energy of the entire Hilbert spectrum plane is denoted as A, which is equally divided into N equal areas, the energy in each small area is denoted as W i (i � 1, 2, . . ., N), and the energy of each small area is normalized; thus, q i � W i /A. According to the formula of information entropy, the formula of HHT energy entropy of wheel acceleration signal is According to the basic properties of energy entropy of HHT, the value of E (q) can reflect the uniformity of the energy in the time-frequency plane. e more uniform the energy distribution in the plane, the larger the value of E (q); otherwise, E (q) becomes smaller.
According to the principle of correlation, the IMF components which are closely related to the original signal are selected for the Hilbert frequency spectrum analysis. Finally, the time-frequency spectrum is divided into 170 equal portions, and the entropy of HHT energy in timefrequency spectrum is calculated. e specific process is shown in Figure 15. e values of HHT energy entropy and their trends under various conditions are shown in Table 10 and Figure 16. It can be seen that the HHT energy entropy of the first condition (one loose fastener) is higher than that of condition 2 and condition 3 (two loose fasteners); in case of 100% looseness, the energy entropy of condition 1 is 4.2672, and those of condition 2 and condition 3 are 4.2043 and 4.2339, respectively. For the same working condition, the HHT energy entropy is less than that without looseness; for example, in working condition 1, the energy entropy decreases from 4.5678 to 4.2672 as looseness increases. erefore, the HTT energy entropy of wheelset vibration signals decreases with the increase of fastener loosening degree. When the fastener is completely loose, the HHT energy entropy is the smallest. is is because the stiffness of the fastener affects the vibration of various components of the vehicle and rail system. As the rigidity decreases, the rail vibration and the rail displacement will change, just like higher irregularity's impact on a vehicle in the running line, which will affect the vibration of the vehicle system, make the time-frequency distribution of the vibration signal inhomogeneous, and cause the value of energy entropy to drop. erefore, calculating HHT energy entropy of wheelset vertical acceleration can reflect the looseness of fastener to a certain extent, and it can provide reference for the identification of fastener looseness.

Conclusions
is paper starts with the idea of identifying the loosening of the fastener from the vibration response signal of the vehicle, and the traditional vehicle-track rigid-body model is improved. e vehicle-track rigid-flexible coupling model is established by using multibody dynamics software SIM-PACK and finite element analysis software ANSYS. e rigid-flexible coupling model simulates the vibration response signals of the wheelsets and the vehicle body in three operating conditions. Combined with improved HHT algorithm and the information entropy theory, the energy entropy is applied to the recognition of the vertical acceleration signal of wheelset with fastener loosening; a certain degree of recognition effects had been achieved. e main conclusions of this paper are listed as follows: (1) Firstly, the finite element model of rail and rail slab is established in ANSYS, and the model is discretized. e substructure modal analysis of rail and rail slab is carried out, and the modeling parameters are obtained. In order to solve the shortcoming that SIMPACK wheel rail module can only realize the physical contact between rigid body and rigid body, a virtual rail without mass and moment of inertia is added between the wheel and the rail, and the transfer of force is realized based on the balance of force and deformation compatibility condition, the fastener and CA mortar are simulated by force element, and the flexible track and rigid vehicle are assembled in the SIMPACK, and a vehicle-track rigid-flexible coupling model is established. Finally, the vertical accuracy of the model is verified by comparing the measured acceleration of quasi-high speed vehicles with different track irregularities.
(2) ree typical lines with loosened fasteners are set up and the five degrees of fastener loosening are set up. Based on the model built in this paper, the vehicle runs on the rails that the fasteners are loosened in varying degrees and under the effect of German lowinterference spectrum at the speed of 200 km/h, and the corresponding vertical acceleration of wheelset and the vibration response of vehicle vertical acceleration are obtained. Finally, the vertical acceleration of the vehicle wheel with obvious vibration response characteristics is selected for the subsequent analysis. (3) e Hilbert-Huang Transform is applied to the analysis of wheelset vibration signals. In order to solve the problems of low orthogonality and energy leakage in the first step EMD of HHT, the EMD is replaced by the orthogonal empirical mode decomposition in processing the vertical acceleration signal of wheelset. Based on this, the correlation  Figure 15: Calculation of the Hilbert energy entropy process. Firstly, varying vertical acceleration signals are gotten and then decomposed by OEMD3. Secondly, IMFS is selected based on correlation principle, which is calculated by Hilbert transform. Lastly, energy entropy theory is used to compute the energy in different looseness conditions.  Figure 16: Energy entropy trend of vertical vibration of wheelset. Condition 1 (black square), condition 2 (red dot), and condition 3 (blue triangle); the looseness is from 0% to 100%, increasing gradually at an interval of 25%.
principle is applied to eliminate the components whose correlation with the original signal in the decomposition results is weak. Owing to the fact that Hilbert time-frequency spectrum of vertical acceleration signals under different working conditions is complicated, the algorithm of Hilbert energy-spectrum entropy is proposed by combining the theory of HHT time-frequency spectrum with the theory of information entropy and is applied to the analysis and identification of wheelset vertical acceleration signals in different conditions. e results show that the method of energy entropy can be used to analyze the vertical vibration signals of the wheelset under different degrees of fastener looseness, and it can realize the identification of the looseness of rail fasteners. e simulation results show that the HHT energy entropy of the wheelset vertical vibration signals decreases with the increase of the degree of fastener loosening; the algorithm can recognize the loose state of rail fastener.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
W.Z. and X.G. conceptualized the study; W.Z., X.G., and S.Z. contributed to methodology; W.Z., Y.W., and L.P. provided software and validated the study; W.Z., L.P., Y.W., and S.Z. investigated the study; X.G. contributed to resources; W.Z., L.P., and S.Z. contributed to data curation; W.Z., L.P., S.Z., and X.G. wrote, reviewed, and edited the manuscript; W.Z. contributed to visualization, writing-original draft preparation, formal analysis, and project administration; L.P. supervised the study; S.Z. acquired funding. All authors have read and agreed to the published version of the manuscript.