Abstract
The fastener system is an important component of the ballastless track, and the fastener looseness threatens the safety of the railway track. This paper proposes a Bayesian model updating method integrated with a Python–Abaqus interface framework to identify the fastener partial looseness of the ballastless track via vibration data. By following the technical standards of Chinese Railway Industry, a laboratory scaled ballastless track model was constructed for the demonstration and verification of the proposed methodology. Not only the damage location can be accurately identified, but also the damage severity can be evaluated. In addition, the associated posterior uncertainties of the identified results can be quantified by calculating the posterior probability density functions of all uncertain model parameters, which can provide valuable information to engineers for ballastless track system damage detection.
Introduction
The fastener system in ballastless track is an important component to guarantee the track geometries and safety. With the rapid development of the high-speed railways, the rail structures are suffering from the reciprocating heavy train loads, which may cause the elastic fasteners to become loose/failed. And the loose/failed fasteners will lead to a large variation in the lateral wheel/rail forces and increase the vibration of rails (Zhou and Shen, 2013; Weng et al., 2008) and reduce the stability of the rail track. What’s more, the fastener looseness will lead to the failure of the adjacent ones (Yuan et al., 2021). Therefore, if the loose fasteners are not detected and repaired in time, the track structure will be further damaged and the safety of the track structure will be threatened (Xu et al., 2017).
At present, regular visual inspection is a common method for the detection of the track fastener damage, especially in a country such as China with the world’s largest high speed railway network. However, visual inspection can only detect the apparent condition of the fasteners, it is difficult to discriminate the slight damage such as fastener partial looseness and clip slight shifting, which is no longer adapted to the development of railway tracks (Aytekin et al., 2015). A variety of loosening detection methods have arisen with the evolution of sensor, signal processing and computer vision technology (Huang et al., 2022), among which the machine vision and artificial neural networks based inspection methods were the major methods utilized for railway fastener detection. Gibert et al. (2017) put forward a new algorithm integrating a multitask learning framework with multiple detectors to detect the defects of railway ties and fasteners, the promising detection performance showed that the proposed approach can improve the detecting accuracy. While this approach is limited to identify invisible defects of the structure. Mao et al. (2018) proposed a promising high-speed railway fastener inspection approach utilizing the structured light sensors to identify the damaged fasteners with a high accuracy over 99.80%. While this method can only effectively identify the location of loose fasteners, but cannot effectively identify the severity of fastener looseness, what’s more, the 3-D point clouds from the structured light sensors are easily affected by some unrelated parts next to the rail track and the cost involved is high. Therefore, the visual-based identification methods can hardly detect the fastener looseness or the elastic clip damage in the early stage, a practical railway fastener looseness detection is significantly demanding nowadays.
Recently, the research on structural damage identification by solving engineering inverse problems has made great progress (Xia et al., 2008; Esfandiari et al., 2010; Zhu et al., 2013). However due to the influence of measurement noise and the limited number of sensors as well as the complexities of boundary conditions, the errors related to constructing a theoretical model of a structure are inevitable, which may lead to the non-convergence in model updating process and high uncertainties of damage identification results (Bakhary et al., 2007). Compared with the deterministic damage identification method, the probabilistic method can better deal with the uncertainties of the damage identification process (Saito et al., 2005). Beck and Katafygiotis (1998) have done pioneering work in the field of Bayesian model updating. Some researchers (Lam et al., 2014, 2018; Hu et al., 2018) proposed a ballast damage detection method based on Bayesian model updating for the first time in railway track. Noted that the mentioned works employ the modal-domain data to identify the uncertain model parameters. For the raw data is in time-domain (e.g., from accelerometers), an extra modal identification process is necessary to transform the time-domain data to modal-domain data, which may be time-consuming and easily cause probable loss of information (Lam et al., 2017). In order to avoid this additional process, a fastener partial looseness identification method based on the time-domain Bayesian model updating integrated with a Python–Abaqus interface framework is proposed for the first time in the realm of structural damage detection in the railway track.
The structure of this paper is outlined below. The modeling of the ballastless track structure is firstly introduced, afterwards the Bayesian approach is proposed and developed for fastener stiffness identification. Then the proposed approach is demonstrated and validated by using the experimental case studies. It turns out that the proposed approach is feasible to properly detect the fastener partial looseness and quantify the uncertainties involved in the identification process simultaneously.
Proposed ballastless track fastener looseness detection method
Modeling of a representative ballastless track system
The ballastless track structure is composed of the rails, fasteners, and the substructure as exhibited in Figure 1(a). The heavy train load was directly supported by the rails, and transformed to the substructure through the fasteners. Fasteners are significant to fix the rails, maintain the gauge distance and reduce the dynamics between the wheel and rail (Xu et al., 2013). A three-dimensional finite element (FE) model (as exhibited in Figure 1(b)) of a scaled ballastless track model in the laboratory (as shown in Figure 4) was established by using the commercial software Abaqus. The substructure (including the prefabricated slab, the cement-emulsified asphalt (CA) mortar layer and the concrete base) was established by using linear hexahedral elements (C3D8R), and tie contact was utilized to simulate the interface contact between the substructure of the ballastless track. The elastic foundation was considered in this study. As the purpose of this study is to investigate the feasibility of the proposed Bayesian model updating method integrated with a Python–Abaqus interface framework to identify the fastener partial looseness of the ballastless track, only half of the track was considered in this study. In the FE model, the length of the track system was set to 1.6 m (along the rails) and the width was set to 1.4 m (perpendicular to the rails), and the dynamic system is assumed to be classically damped and the damping ratio of the track system was assumed to be 8%, which was estimated by modal analysis of the measured impact hammer test. The fasteners were simulated by using the linear elastic spring elements, and the masses of the fasteners are ignorable compared to that of the whole vibration system. WJ-8 type fastener with the gauge apron plate (seen in Figure 4) was considered and the transverse and vertical stiffness of the fastener system were simulated by using spring sets. The spring stiffness reduction was considered to describe the looseness of the fastener in this study. In order to detect the fastener looseness, the scaled model (Hu et al., 2021) with the fastener system was utilized in this study, and more details about the values of the stiffness of the prefabricated slab, CA mortar layer and concrete base can be referred to Hu et al. (2021). The ballastless track structure. (a) Schematic of a representative ballastless track system. (b) A three-dimensional finite element model of the scaled ballastless track system.
The parameters which are likely to change by local conditions, load patterns, or damage should be considered as uncertain parameters (Lam et al., 2020; Adeagbo et al., 2021, 2022). The loose fastener only can influence the fastener stiffness a lot, and the values of the model parameters of the substructure (including the prefabricated slab, CA mortar layer and concrete base as well as the foundation) would not be changed too much, which can be referred to Hu et al. (2021), therefore, only the scaling factors of the fastener stiffness will be treated as the uncertain model parameters in this study. Impact force (as shown in Figure 2) was utilized to calculate the dynamic acceleration responses and was obtained from the impact hammer test on the scaled track model in the laboratory with the value of approximately 20,000N. The linear perturbation analysis was employed to calculate the dynamic responses of the FE model under impact excitations. Fastener looseness leads to the fastening system less stiff, which changes the dynamic characteristics of the track system (Zhan et al., 2020), therefore, fastener looseness identification can be treated as an inverse problem, in which the measured vibrational data is utilized in the process of back-calculating the fastener stiffness. In this study, only the vertical vibration of the track model and the uncertain model parameters of the fastening system were considered in the following analysis. A typical impact excitation.
Bayesian model updating
The main objective of Bayesian probabilistic approach is to calculate the posterior probability density function (PDF) of the uncertain model parameters
As the model class
Fastener looseness identification methodology
An overview of the proposed fastener partial looseness identification method of the ballastless track system based on a time-domain Bayesian model updating method integrated with a Python–Abaqus interface framework is schematically illustrated in Figure 3. The FEM of the ballastless track system was established in Abaqus, which is commonly used by engineers for structural analysis and calculation. And the measured impact force from the laboratory impact hammer test was put into the FE model. The linear perturbation analysis was adopted for the dynamic time domain analysis in this study to obtain the vibrational responses of the structure under impact excitation. Impact hammer tests were carried out on the target scaled ballastless track model in laboratory (as shown in Figure 4) to obtain the time domain responses (e.g., acceleration responses at different DOFs). A newly developed Python–Abaqus interface framework was firstly introduced for railway fastener partial looseness identification. Secondary development of the Abaqus model based on Python scripts was utilized to ensure the seamless interaction between the Bayesian model updating algorithm and FEM such that the read and write process of model-predicted time domain responses with different groups of uncertain model parameters can be easily implemented. After the Bayesian model updating was completed, the most probable model parameter vector Flowchart of the proposed partial fastener looseness identification method. The scaled model of the ballastless track structure.

Experimental case studies
Experimental setup
Main parameters of the ballastless track structure (Hu et al., 2021; Zhang et al., 2019).
Figure 6 shows the measurement equipment used in the impact hammer tests. A torque wrench (shown in Figure 6(a)) was used to set the preload torque of the fastener system. The impact hammer (Type 8210, seen in Figure 6(b)) with a sensitivity of 0.1985 mV/N and maximum force of 22,200 N was used to excite the track structure with different tips, softer tips were more suitable for exciting the lower modes of the system and stiffer tips were more appropriate for exciting the higher modes. And the medium hammer tip with the green color was employed, which was suitable for the current scaled model of the slab track, and the excitation location was shown in Figure 5. Six unidirectional IEPE sensors (with sensitivities around 100 mV/g, seen in Figure 6(c)) were placed on the rail seat near to the fasteners (as shown in Figure 5) to measure the vertical vibration responses of the track structure. The measured vibration signals from the accelerometers were transferred to the data acquisition equipment (Dong-Hua DH5922, seen in Figure 6(d)), and the equipment was connected to a computer with DHDAS dynamic signal acquisition and analysis system. The sampling frequency was set as 1000 Hz in the impact hammer tests and approximately five impulses were recorded for each set of measurements. A typical measured time-domain vibration response from one accelerometer is shown in Figure 7(a), and the fourth impulse is zoomed and shown in Figure 7(b), which was damped to zero in 0.2 sec implying that the damping of the system is relatively high. Schematic of the sensor locations and impact location in the impact hammer test. Experimental equipment. A typical measured time-domain vibration response. (a) A typical set of measured acceleration response. (b) zoomed response of the fourth impulse.


Health conditions of the fastener system.
Undamaged case
The elastic stiffness of the three fasteners is defined by the multiplication of its nondimensional scaling factor,
Bayesian model updating results in undamaged and damaged cases.
Figure 8 shows the matching between the model-predicted and measured time-domain vibrational responses from the six accelerometers in the undamaged case. It is obvious that the matching is relatively good. It is believed that the accuracy of the identified FE model is good enough to help in response prediction of the track structure, which is extremely significant in the process of model updating. Matching between the measured and model-predicted responses in the undamaged case.
Damaged cases
Then three damaged cases are utilized to test the proposed Bayesian fastener partial looseness identification methodology. The Bayesian model updating results of the three damaged cases are listed in Table 3. In damaged Case 1, the most probable model parameter vector is Matching between the model-predicted and measured responses in damaged cases. (a) Damaged Case 1. (b) Damaged Case 2.
In damaged Case 2, the identified most probable model parameter vector is
Uncertainty analysis
Generally, the less the amount of information available from the measurement, the larger the posterior uncertainties associated with the model parameters. Figure 10 shows the normalized marginal PDFs of the model parameters in the damaged cases 1 and 2. It can be seen that the posterior distributions of the scaling factor The normalized marginal PDFs. (a) Damaged Case 1. (b) Damaged Case 2.
Conclusions
This paper investigates the feasibility study on the identification of the fastener partial looseness of the ballastless track based on time-domain Bayesian model updating integrated with a Python–Abaqus interface framework using measured vibrational responses. A scaled ballastless track model was built in the laboratory for the verification and demonstration of the proposed methodology. The experimental case studies reveal that: (1) The proposed Bayesian probabilistic approach integrated with the Python–Abaqus interface framework is feasible to identify the fastener partial looseness of the ballastless track by calculating the most probable values (MPVs) of fastener stiffness scaling factors. (2) Not only the looseness location can be successfully detected, but also the looseness severity can be identified synchronously. (3) The uncertainties associated with the fastener looseness identification results can be properly quantified by calculating the posterior PDFs (and COVs) of all uncertain model parameters. (4) The closer the impact location is away from the sensor, the lower the uncertainties and the higher the accuracies of the corresponding identified model parameters are.
It can be concluded that the proposed method has a high potential of application in fastener partial looseness identification of the ballastless track, which can provide help for the detection of similar structural damage and make the health monitoring and early warning of the ballast track fastener system possible.
In this study, only three pairs of fasteners were considered in the experimental case studies, further work should be conducted on the field tests to ensure that the proposed methodology can be put into industrial applications, and not only the accelerations of the slab track, but the receptances of the rail can be utilized in the model updating in future work. It’s believed that the optimal locations and accelerometer(s) and impact are important to successfully detect the fastener partial looseness, which can be studied under the Bayesian framework by calculating the information entropy (Chow et al., 2011; Lam and Adeagbo, 2022). In addition, impact hammer tests were employed to collect vibrational measurement data of the ballastless track structure, for continuous monitoring of a long span of ballastless track, it is proposed to extend the proposed approach to utilize train-induced vibration for the identification of the fastener stiffness. As thousands of kilometers of the ballastless track are constructed of the same specification, a proper fastener looseness identification methodology can be repeatedly utilized along the entire railway line. Therefore, the long-term potential of the proposed approach cannot be underrated.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the This work was fully supported by the grant from the National Natural Science Foundation of China (52178287, 51708242); and Natural Science Foundation of Hubei Province (2020CFA047).
