
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38622239
59328
10.1038/s41598-024-59328-5
Article
Study on monitoring broken rails of heavy haul railway based on ultrasonic guided wave
Xu Xining xuxining@bjtu.edu.cn

124
Wen Ziyu 2
Ni Yi 3
Shao Bohuai 2
Ma Xinyu 2
Pan Zecheng 2
1 https://ror.org/01yj56c84 grid.181531.f 0000 0004 1789 9622 State Key Laboratory of Advanced Rail Autonomous Operation, Beijing Jiaotong University, Beijing, 100044 China
2 https://ror.org/01yj56c84 grid.181531.f 0000 0004 1789 9622 School of Mechanical, Electronic, and Control Engineering, Beijing Jiaotong University, Beijing, 100044 China
3 grid.464214.1 0000 0001 1860 7263 Railway Engineering Research Institute, China Academy of Railway Sciences Corporation Limited, Beijing, 100081 China
4 grid.181531.f 0000 0004 1789 9622 Frontiers Science Center for Smart High-Speed Railway System, Beijing, 100044 China
15 4 2024
15 4 2024
2024
14 866722 1 2024
9 4 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Real-time monitoring of broken rails in heavy haul railways is crucial for ensuring the safe operation of railway lines. U78CrV steel is a common material used for heavy haul line rails in China. In this study, the semi-analytical finite element (SAFE) method is employed to calculate the dispersion curves and modal shapes of ultrasonic guided waves in U78CrV heavy steel rails. Guided wave modes that are suitable for detecting rail cracks across the entire cross-section are selected based on the total energy of each mode and the vibration energy in the rail head, web, and foot. The excitation method for the chosen mode is determined by analyzing the energy distribution of the mode shape on the rail surface. The ultrasonic guided wave (UGW) signal in the rail is analyzed using ANSYS three-dimensional finite element simulation. The group velocity of the primary mode in the guided wave signal is obtained through continuous wavelet transform to confirm the existence of the selected mode. It is validated that the selected mode can be excited by examining the similarity between the vibration shapes of a specific rail section and all modal vibration shapes obtained through SAFE. Through simulation and field verification, the guided wave mode selected and excited in this study demonstrates good sensitivity to cracks at the rail head, web, and foot, and it can propagate over distances exceeding 1 km in the rail. By detecting the reflected signal of the selected mode or the degree of attenuation of the transmitted wave, long-distance monitoring of broken rails in heavy-haul railway tracks can be achieved.

Keywords

Heavy haul railway
Ultrasonic guided wave
Monitoring broken rails
Mode selection
Mode excitation
Subject terms

Mechanical engineering
Characterization and analytical techniques
Electrical and electronic engineering
Fundamental Research Funds for the Central Universities of China2022JBQY002 Xu Xining Science and technology innovation project of Shuohuang Railway Development Co., LtdGJNY-21-65 Xu Xining issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Heavy haul railway transportation offers substantial advantages such as large transportation capacity, high efficiency, and low transportation costs, greatly fostering the development of China's national economy1. China's heavy haul railway mainly uses seamless rail; the phenomenon of thermal expansion and cold contraction of seamless rail is pronounced, coupled with the wheel-rail force of heavy haul railway, it is easy to lead to rail fracture. In addition, rail top abrasion, rail core damage, rail splint clamp damage, rail welding quality, and other problems can also lead to rail fracture2.

Currently, rail fracture detection methods are primarily categorized into mobile detection and fixed detection. Mobile detection involves assessing rail conditions using hand-propelled rail detection trolleys or large rail inspection vehicles. The detection methods mainly consist of ultrasonic, eddy currents, and magnetic flux leakage. Within ultrasonic detection, there are piezoelectric and electromagnetic methods. China's rail detection vehicle is mainly developed based on piezoelectric ultrasonic technology. For instance, the GT-2 rail detection trolley utilizes seven ultrasonic transducers with varying angles to simultaneously emit and receive ultrasonic waves on the rail for internal crack detection3. Based on electromagnetic ultrasonic detection technology, Canada's Tektrend company developed the RailPro rail detection vehicle4, which can detect defects in the rail head and web. The rail flaw detection system developed by VIGOR company5 in Russia realizes the detection of rail defects. Sperry company in the United States has developed a hand-pushed rail surface crack detector6 and a high-speed flaw detection vehicle7 based on eddy current detection technology and developed a rail defect detection system combined with ultrasonic detection technology and magnetic flux leakage detection technology, with an inspection speed of 35 km/h6. The mobile detection method has simple operation and high detection sensitivity. The advantage of mobile detection methods lies in their simple operation and high detection sensitivity, allowing for quick identification of rail defects. However, due to technological limitations, the detection distance of mobile detection methods is limited, and some areas or corners may not be fully covered and detected. Additionally, because operations must be carried out during idle periods of train transportation, mobile detection cannot monitor the condition of rails in real-time. This, to a certain extent, affects the ability to timely detect and repair rail issues.

Fixed detection refers to the installation of specific detection equipment for continuous online monitoring of rail conditions under all weather conditions. Currently, commonly used detection methods include strain gauge, optical fiber, ultrasonic guided wave, and track circuit8. The detection principle of strain gauge broken rail involves judging whether the rail is broken by detecting changes in rail stress. The detection distance is generally 30–60 m9, which is unsuitable for long-distance broken rail detection. Moreover, this method requires calibrate the position of the zero-pressure point before sensor installation10, which is not applicable to established lines. The optical fiber rail breaking detection involves installing optical fibers on the rail web; when the rail breaks, the optical fiber will also break so that the optical receiver cannot receive the effective signals. This method has high detection sensitivity, but due to the fragile nature of optical fibers, its installation and maintenance process is complex and delicate. It can only be used in short-range or difficult electrical isolation situations11,12, which is unsuitable for long-distance rail breaking monitoring. The track circuit detects the rail breaking in real-time by monitoring the changes in rail voltage, current, and other parameters, which have high accuracy and real-time performance. However, the track circuit cannot detect the partial fracture of the rail, and its detection sensitivity is low in the section with poor track bed condition. In summary, while fixed detection methods can meet the needs of online monitoring, they still have limitations. New full-section online monitoring technologies are needed for practical applications.

Ultrasonic guided wave is a unique ultrasonic wave generated through continuous reflection, refraction, and conversion of the longitudinal and transverse wave with the medium boundary when the ultrasonic wave propagates in the waveguide media13. Its principle for rail breaking detection is relatively simple: when the rail is completely broken, the transducer cannot receive the guided wave signal, and when the rail is partially damaged, the energy of the guided wave signal received will attenuate14. Compared to other detection methods, UGW detection has advantages such as long-distance propagation, covering the entire section of the rail, fast detection speed, and its detection process does not interfere with train transportation, making it suitable for long-distance rail breaking monitoring.

Lamb15 discovered the existence of Lamb waves and confirmed the multimodality in the early twentieth century, many scholars have conducted in-depth research on the propagation characteristics of UGW. In 1950, Thomson first proposed to use the matrix analytical method to solve the modal response results of UGW in layered solid media16. In the following decades, the solution method using matrix calculation was widely used to solve the propagation characteristics of UGW17–19. In the twenty-first century, as the waveguide medium’s cross-section becomes more complex, the solution method through matrix calculation cannot meet the needs gradually. Scholars have proposed various ways to study the propagation characteristics of the guided wave, including the finite element method, the semi-analytical finite element method, the boundary element method, etc.20. Among them, using the SAFE method to solve the dispersion curve of UGW in a waveguide medium has gradually become a mature method. Many scholars have done much research based on this and have made good progress. The following articles discussed numerical methods in SHM or wave propagation analysis. In 2002, Mukdadi et al.21 calculated the dispersion curve of a guided wave in a multi-layer rectangular section beam structure. In 2003, Hayashi et al.22 solved the dispersion curve of the guided waves in a waveguide medium with an arbitrary section. In 2017, Cong et al.23 solved the dispersion curve of UGW in spiral fin. In 2018, Duan and Kirby24 in the UK combined the SAFE method with the perfect matching layer to analyze the propagation characteristics of ultrasonic guided waves in fluid-filled pipes. They developed a SAFE method in 2019 to study the propagation characteristics of ultrasonic guided waves in composite laminates25.

In recent years, the application of UGW in railway nondestructive testing has been fully developed. In 2009, Chong et al.26 used the mixed design of analytical method and finite element method to design a guided wave sensor to detect the defects inside the rail so that the guided wave can propagate a longer distance and proposed the mixed analytical method of guided wave and finite element, which is used to analyze the propagation characteristics of guided wave in the rail. In 2011, Coccia et al.27 solved the free vibration solution and forced vibration solution of UGW propagation at the rail head under laser excitation. In 2014, Xu et al.28 solved the dispersion curve of UGW in CHN60 rail using the SAFE method and analyzed four modes at low frequency. In 2015, Loveday and Long29 used a laser vibrometer to detect the propagation ability of each UGW mode in the rail at a low frequency on the 5 m long rail in the laboratory and selected the guided wave mode that can propagate over a long distance. In 2018, Setshedi et al.30 estimated rail material and geometric characteristics based on SAFE and UGW and used Poisson's ratio and three geometric parameters to represent rail contour loss. In the same year, Ramatlo et al.31 designed the UGW transducer for detecting rail cracks by using the SAFE method and three-dimensional finite element method combined with a mathematical optimization algorithm. In 2020, they developed a modeling framework for the ultrasonic detection of waveguides with arbitrary discontinuities to simulate the detection of multiple tracks with discontinuities by ultrasonic guided waves32. With the continuous development of artificial intelligence technology, many researchers have begun to apply machine learning to railways33. In 2022, Harsh Mahajan et al.34 explored the application of UGW generated by pasting a piezoelectric sensor on the rail surface in high-frequency rail damage detection and proposed a machine learning algorithm framework for rail head damage detection. The experimental results show that the maximum error range of the algorithm framework in predicting damage is 2% to 16.16%. In 2022, Sresakoolchai et al.35 utilized track geometry data obtained by a track geometry car (TGC) to detect rail, switch and crossing, fastener, and rail joint defects. They developed accurate detection models by employing supervised and unsupervised machine learning techniques.

U78CrV rail is a wear-resistant heavy haul railway developed by Angang Steel Group Co. LTD. With good comprehensive performance, it has been widely laid and used in China since 200936. It is one of China's most commonly used rail materials for heavy haul railways. Currently, there is a lack of mature and reliable online broken rail monitoring system for the long-distance broken rail monitoring problem in heavy-haul railways. However, the advantages of ultrasonic guided wave testing, such as long detection distance and fast detection speed, make it have broad application prospects in the field of long-distance broken rail detection.

In recent years, the laboratory team has conducted extensive research in detecting rail fractures and rail stress using UGW14,20,28,35–38. Aiming at the application of UGW in the rail break detection of heavy haul railways in China, this paper combines the semi-analytical finite element method and three-dimensional finite element method to identify a guided wave mode suitable for the long-distance detection of broken rails made of U78CrV type heavy rail, as well as the excitation method for the selected mode. It is demonstrated that the excitation method can successfully induce the chosen mode by verifying the consistency of group velocity and the similarity of vibration shape. Simulation experiments show that the selected mode exhibits good sensitivity to the rail head, rail web, and rail foot. Finally, an attenuation test confirms that the chosen mode is suitable for long-distance rail break detection in heavy haul railways.

Dispersion curve of guided wave in U78CrV rail

According to the dimensional specifications of the 75 kg/m rail, the cross-section of the rail is drawn in SolidWorks. The rail section height is 192 mm, the rail head width is 75 mm, and the rail foot width is 150 mm. The coordinate data of the section is imported into MATLAB and plotted to represent the rail cross-section Ω. The cross-section is then discretized into triangular grid elements Ωe using the PDE discretization tool of MATLAB. The division effect is shown in Fig. 1. The section is divided into 550 units consisting of 340 nodes.Figure 1 Discrete drawing of 75kg/m rail section.

The direction of guided wave propagation is defined as the z-direction, and the rail section is defined as the Oxy plane, where the y-direction represents the vertical direction. According to the SAFE theory38, the harmonic displacement, stress, and strain field components of each point in the rail can be expressed as:1 u=uxuyuzTσ=σxσyσzσyzσxzσxyTε=εxεyεzγyzγxzγxyT

The relationship between stress and strain is given by σ = C ε, and C is the constant elastic matrix of rail.

The displacement field can express the strain at any point in the rail as:2 ε=Lx∂u∂x+Ly∂u∂y+Lz∂u∂z

where:Lx=100000000010000001,Ly=000010000100001000,Lz=000000001000010100

The amplitude variation of the section displacement field can be expressed by Eq. (3).3 ux,y,z,t=uxx,y,z,tuyx,y,z,tuzx,y,z,t=Uxx,yUyx,yUzx,yeiξz-ωt

where ξ is the wave number in the z-direction, and ω is the angular frequency.

According to Eq. (4), the displacement of any point within the discrete triangular element can be expressed using the shape function as:4 uex,y,z,t=∑k=13Nkx,yUxk∑k=13Nkx,yUyk∑k=13Nkx,yUzkeeiξz-ωt=Nx,yqeeiξz-ωt

where N(x, y) is the shape function matrix, q(e) is the node displacement vector. The shape function of each node can be calculated from the coordinates of the node in the section39. Based on the Hamiltonian principle, the general homogeneous wave equation of UGW can be deduced39:5 K1+iξK2+ξ2K-ω2MU=0

In the equations, K1, K2, K3 are three stiffness matrices of size 3n × 3n, where n is the total number of nodes. These matrices can be calculated directly according to the material parameters of U78CrV. M and U are the mass matrix and displacement matrix, respectively. The wave equation precisely describes the dispersion characteristics and multimodality of UGW. By introducing an auxiliary matrix of size 3n × 3n, the imaginary part in Eq. (5) can be eliminated.T=i11⋱i11

After introducing the auxiliary matrix, Eq. (5) becomes:6 K1+ξK^2+ξ2K3-ω2MU^=0

In the equations, K^2=TTK2T-i,U^=TU. Equation (7) can be reformulated as a first-order eigenvalue system:7 A-ξBU^ξU^=0

In the equations:A=0K1-ω2MK1-ω2MK^2,B=K1-ω2M00-K3

For each ω, we can obtain 6n eigenvalues of the wavenumber ξ from Eq. (7). Among these, the real eigenvalues correspond to propagating guided wave modes, while the complex eigenvalues correspond to evanescent waves. By iterating through the values of ω, we can obtain the ξ − ω relationship, from which we can solve for the phase velocity dispersion curve of ultrasonic guided waves. The phase velocity Cp is defined as follows:8 Cp=ωξ

Based on Eq. (8), we can obtain the phase velocity dispersion curve of ultrasonic guided waves in the cross-section of U78CrV steel rails, as shown in Fig. 2.Figure 2 Dispersion curve of phase velocity in U78CrV rail.

From the dispersion curve, it can be observed that there are multiple guided wave modes in the steel rail. As the frequency of the guided waves changes, the modes of the guided waves also change. Moreover, with an increase in frequency, the number of guided wave modes increases.

Group velocity is defined as the derivative of angular frequency to wave number :9 Cg=dωdξ

And finally, the calculation formula of group velocity is28:10 Cg=∂ω∂ξ=U^LTK^2+2ξK3U^R2ωU^LTMU^R

where K^2 is a symmetric matrix for undamped motion, U^ is a new nodal displacement vector, U^R represents the right eigenvector and U^L represents the left eigenvector ω,ξ,U^L,U^R can be obtained when solving the characteristic equation. The dispersion curve of guided wave group velocity in U78CrV rail is obtained by taking them into Eq. (10), as shown in Fig. 3.Figure 3 Group velocity dispersion curve in U78CrV rail.

Based on the node displacement information in the eigenvectors solved in the previous step, the vibration mode shapes of all guided wave modes can be plotted. Taking the guided wave frequency of 35kHz as an example, there are 23 propagation modes in the rail at this frequency, and the vibration shape diagram of each mode is shown in Fig. 4.Figure 4 Guided wave modes in rails at 35 kHz.

Mode selection

According to the previous research40, when using UGW to detect the long-distance rail breakage of U78CrV rail, selecting 35kHz as the guided wave frequency has a good detection effect. To ensure that the exciting mode can cover the whole section of the rail, the mode requires a large vibration amplitude at the rail’s head, web, and foot.

At the frequency of 35kHz, there are 23 guided wave modes in the rail. The rail section height of 138–192 mm is defined as the rail head area, 33–138 mm as the rail web area, and 0–33 mm as the rail foot area. Define the vibration displacement of the No. i node of mode m as (uxim,uyim,uzim), and the average vibration energy of mode m is defined as shown in Eq. (8).11 Em=1n∑i=1nuxim2+uyim2+uzim2

According to the vibration data of all nodes in each area, the average vibration energy of 23 modes at the rail head, rail web, and rail foot, as well as the vibration energy of each mode in x, y, and z-directions are calculated respectively. The corresponding results are presented in the form of histograms, as shown in Fig. 5.Figure 5 Vibration energy distribution of 23 guided wave modes.

We selected the five modes with the highest total energy for in-depth analysis. Theses modes are mode 7, mode 9, mode 12, mode 16 and mode 18. The mode shapes are shown in Fig. 6, where Cp is the phase velocity of the mode in m/s.Figure 6 Vibration shapes diagram of the selected five modes.

Observing the vibration shapes of the five modes, it can be seen that mode 7, mode 12, and mode 16 exhibit many distorted forms, making them difficult to excite in the field environment. Additionally, the energy of mode 18 is very limited at the rail head and rail web, which may result in the omission of defects. Therefore, mode 9 is finally selected. The modes in the dispersion curve are separated, and the phase velocity and group velocity dispersion curves of mode 9 are drawn separately, as shown in Fig. 7. The blue line and red line in the figure represent the phase velocity and group velocity curves of mode 9, respectively. When the frequency is 35kHz, the phase velocity and group velocity of mode 9 solved by SAFE are 3571m/s and 2669m/s, respectively.Figure 7 Phase velocity and group velocity dispersion curve of mode 9.

According to the field environment, the nodes that the UGW transducer can be installed on the outer surface of the rail section are selected, as shown in Fig. 8. There are 28 nodes symmetrically on the rail surface that meet the requirements. To choose the best excitation point, the vibration displacement of mode 9 in the x, y, and z-directions at 28 nodes is calculated, and the histogram is plotted, as shown in Fig. 9. The point with a large displacement is selected as the excitation point. It can be seen from the figure that the displacement of mode 9 in the z direction is the largest at nodes 118 and 84. Select node 84 (x = 10.4mm, y = 88mm) to excite in the z-direction (longitudinal direction of the rail). The following ANSYS simulation analysis is conducted to verify whether mode 9 can be excited here.Figure 8 The node position where the transducer can be installed.

Figure 9 Vibration displacement of mode 9 at each node.

Mode verification

ANSYS is a commonly used three-dimensional simulation and analysis software. The rail vibration data are obtained through the transient dynamic simulation calculation of ANSYS, the group velocity is solved by wavelet transform, and the vibration shape of the rail section is analyzed to determine the guided wave mode in the vibration signal.

Mode verification based on group velocity

A 5m long three-dimensional rail model is established by SolidWorks, as shown in Fig. 10.Figure 10 3D model of 5m long U78CrV rail.

The rail model is meshed using Hypermesh software, and each section is divided into 480 nodes. The meshing effect is shown in Fig. 11.Figure 11 3D mesh generation effect of rail.

The material properties of the U78CrV rail are assigned to all grid units, and the rail material parameters are shown in Table 1. Table 1 Material parameters of U78CrV rail.

Model	Density (kg/m3)	Young's modulus (Pa)	Poisson's ratio	
U78CrV	7850	2.08e11	0.33	

The established mesh model is imported into the ANSYS software for transient dynamic analysis. During the simulation, the element type used is the 3D solid element SOLID45 with 8 nodes. The material properties are as follows: material density is 7850 kg/m3, elastic modulus E is 208 GPa, Poisson's ratio ν is 0.33. The material is modeled as an isotropic material, MAT1, with properties that do not change with temperature. The excited ultrasonic guided wave modes are analyzed using ANSYS. According to the analysis results in section "Mode selection", at the middle of the rail web(x = 10.4mm, y = 88mm), a five-cycle sinusoidal signal modulated by the Hanning window is applied longitudinally along the rail, as shown in Fig. 12. The left figure shows the position of the excitation point and the right figure shows the excitation signal.Figure 12 Excitation position (the left) and excitation signal (the right).

The deformation process of rail with time obtained by simulation is shown in Fig. 13.Figure 13 Deformation process of rail.

The vibration data of 240 nodes at the same horizontal position as the excitation point and within the range of 1184–3104 mm from the excitation point are extracted, and continuous wavelet transform is performed on the data to obtain the wavelet time–frequency diagram, as shown in Fig. 14.Figure 14 Wavelet transform diagram.

The "distance" in the figure represents the distance from the receiving node to the excitation point. The point with the highest frequency component is marked in the figure. It can be seen from the figure that the frequency of the signal is concentrated near 35kHz, and the calculated speed of the guided wave energy concentration point reaching each node is shown in Table 2, which is very close to the group speed of 2669m/s of mode 9 calculated by the SAFE method. Therefore, the mode propagating in the rail is mainly mode 9. Table 2 Group speed of the guided wave.

Distance (mm)	Time (us)	Group velocity (m/s)	Bias (%)	
1312	491	2672.10	0.12	
1824	673	2710.25	1.55	
2336	879	2657.57	 − 0.43	
2592	964	2688.80	0.74	

Mode verification based on mode shapes

The Sect. 2080mm from the excitation point was selected as a specific cross-section for mode validation, named Ω. The theoretical group velocity of mode 9 solved by the SAFE method is 2669 m/s. According to the group velocity, it can be calculated that the time for mode 9 to reach the section Ω is about 0.000779s.t=2080mm2669m/s=0.000779s

The vibration shapes of section Ω near the time of T = 0.000779 s are drawn, as shown in Fig. 15. The vibration shape at t = 0.000773 is selected for analysis after comparison.Figure 15 Vibration shapes diagram of section Ω near 0.000779 s.

The vibration shapes of 23 modes existing in the rail at 35kHz solved by the SAFE method are drawn and compared with the vibration shapes of section Ω at t = 0.000773s. The comparison diagram is shown in Fig. 16, with the vibration shape of section Ω at 0.000773 s in the right part of each subplot and the vibration shapes of each mode in the left part.Figure 16 Comparison diagram of vibration shapes of section Ω and various modes at 0.000773 s.

The 340 nodes of the section are numbered, and the vibration data of the section Ω at that time are extracted according to the three directions of x, y, and z and formed into the vector qΩ:12 qΩ=ux1Ωuy1Ωuz1Ω⋯uxnΩuynΩuznΩ

Similarly, the vibration data of the 23 modes are extracted in the x, y, and z-directions to form one-dimensional data and normalized to form the vector qm:13 qm=ux1muy1muz1m⋯uxnmuynmuznm

The difference in vibration between the section Ω and the mode m is expressed in terms of the Manhattan distance, i.e.14 DΩ,m=∑i=1nuxim-uxiΩ+uyim-uyiΩ+uzim-uziΩ

A difference value is obtained for each mode, and the 23 difference values are sorted in ascending order to get the results shown in Table 3. Table 3 Difference value (D-value) between section Ω and the vibration shape of each mode at 0.000773 s.

Mode	9	2	1	5	
D-value	256.1907	263.3642	266.8521	270.5071	
Mode	3	8	14	15	
D-value	273.2616	288.3123	323.6595	338.7663	
Mode	4	18	11	19	
D-value	347.0474	368.7845	371.1023	377.0815	
Mode	23	13	22	6	
D-value	388.8818	392.8972	396.9228	399.5571	
Mode	12	7	10	21	
D-value	403.0213	406.0222	412.2435	414.2943	
Mode	17	20	16		
D-value	418.8620	425.0082	481.8620		

It can be seen from the table that when t = 0.000773, the difference between the vibration of section Ω and mode 9 is the smallest; that is, the vibration shape of section Ω is most similar to that of mode 9. Figure 17 shows the comparison between the vibration of mode 9 and the vibration of section Ω at the moment 0.000773s.Figure 17 Comparison diagram of vibration shape of section Ω and mode 9 at 0.000773 s.

By comparing the group velocity and vibration shape diagrams, it can be demonstrated that the guided wave mode shown in the mode 9 vibration shape can be excited in the 75 kg/m rail by applying a five-cycle sinusoidal signal modulated by the Hanning window along the longitudinal rail direction at the node position (x = 10.4 mm, y = 88 mm).

Simulation analysis of rail defect detection

To further verify whether the mode excited at the node (x = 10.4mm, y = 88mm) has good sensitivity to the defects at the rail’s head, web, and foot, finite element simulation experiments were conducted (Supplementary Information). The exciting, receiving, and defect positions for the finite element simulation are shown in Fig. 18.Figure 18 Exciting, receiving, and defect points for ANSYS simulation.

The guided wave excitation point T1 is located at the rail waist of the rail end, and the receiving point R1 is set 1m away from the excitation point to receive the echo signal generated by the guided wave at the defect. At a distance of 2m from the T1 point, the rail head, web, and foot defects are set, respectively. A receiving point R2 is set at 4m to receive the transmitted wave signal after the guided wave passes through the defect.

The five-cycle sinusoidal signal modulated by the Hanning window is applied at the rail web with a center frequency of 35 kHz. The simulation steps are shown in Table 4; four simulations were conducted, the first simulation was performed on the intact rail, and the remaining three were conducted with defects at the rail head, rail web, and rail foot, respectively. Table 4 Defect settings of simulation experiment.

Number	Defect position	Defect size (mm)	
1	No defect	0	
2	Rail head	50 × 18 × 8	
3	Rail web	10 × 50 × 8	
4	Rail foot	50 × 18 × 8	

The crack settings are shown in Fig. 19.Figure 19 Settings of rail cracks at different positions.

Four groups of signals are collected at the R1 point, the first group is the signal collected on the intact rail, and the second to fourth groups are echo signals of T1 excited guided wave signals after encountering defects at the rail head, rail web, and rail foot. By comparing the difference between the signals of the second to fourth groups and the signals of the first group, the sensitivity of the guided wave excited at the T1 position to the defects in different places was analyzed.

The left plot in Fig. 20 shows the crack reflection signals received by R1 when the defect is at the rail head, rail web and rail foot, respectively, indicated by the red line. It is compared with the reflected wave at the corresponding position when there is no defect. The right side shows the result of the difference made between the two sets of signals.Figure 20 Comparison and differential analysis of reflected waves from defect and no defect at different positions.

It can be seen that the guided wave signals excited at point T1 show clear echo signals at the locations of cracks on the rail head, web, and foot.

Four groups of signals are collected at the R2 point, the first group is the signal collected on the intact rail, and the 2nd–4th groups are the transmitted waves of the guided wave signal excited at T1 after passing through the defects at the rail head rail web and rail foot, respectively. Comparing the signal attenuation of groups 2–4 with group 1 of the T1-excited guided wave signal after passing through the defects at different locations is analyzed.

The left plot in Fig. 21 shows the crack transmission signals received by R2 when the defect is at the rail head, rail web and rail foot, respectively, indicated by the red line. It is also compared with the transmitted wave at the corresponding position when there is no defect. The right side shows the result of the difference made between the two sets of signals.Figure 21 Comparison and differential analysis of transmitted waves from defect and no defect at different positions.

Based on the signal analysis from Fig. 20, the selected modes exhibit reflection coefficients of 0.105, 0.21, and 0.15 for defects at the rail head, rail web, and rail foot, respectively. This indicates that the chosen modes are sensitive to defects at these locations. Additionally, it can be observed from the figures that the strength of the defect echoes follows the order of rail web, rail foot, and rail head, consistent with the energy analysis shown in Fig. 5, further confirming that the excited mode primarily corresponds to Mode 9.

From Fig. 21, it can be observed that the selected modes exhibit different degrees of attenuation when passing through defects at different positions. Considering both the reflection and attenuation signals, by applying the excitation signal at the node positions x = 10.4mm and y = 88mm ,it is possible to detect significant cracks in the rail by analyzing the echo of the defect or analyzing the attenuation of the transmitted wave..

Experimental results

Field experiments are conducted to further validate the effectiveness of the excitation mode applied to the actual rail line. The installation of the sandwich piezoelectric transducers in the field is illustrated in Fig. 22.Figure 22 Transducer installation diagram.

Firstly, the defect detection experiment was carried out. It is impossible to make holes or cut out other defects on the actual operating line. Therefore, the preset earthing hole on the line was taken as the defect point. The line earthing holes are shown in Fig. 23.Figure 23 Grounding holes for simulating defects.

The experimental setup is shown in Fig. 24. The exciting transducer is installed at point t1, as shown in Fig. 22. The receiving transducer r1 is 6m away from the exciting transducer, the first earthing hole h1 is 18.8m away from the exciting transducer, and the second grounding hole h2 is 119.13m away from the exciting transducer.Figure 24 Experimental line setup.

The field signal collected by the oscilloscope is shown in Fig. 25. The echo signals of defects h1 and h2 can be observed.Figure 25 Field signal waveform.

To verify the long-distance propagation performance of guided wave modes excited at the rail web, a 1km attenuation experiment was carried out on the Beijing Circular Railway (Supplementary Information). The transducer was installed at the rail web, and the excitation signal shown in Fig. 12 was applied along the longitudinal direction of the rail.

To verify the long-distance propagation performance of guided wave modes excited at the rail web, an attenuation experiment was conducted on the Beijing Circular Railway with a radius of 1.5 km. The transducer was installed at the rail web, and the excitation signal shown in Fig. 12 was applied along the longitudinal direction of the rail.

In the experiment, the amplitude of the excitation signal is ± 125V, and a receiving point is set at every 50m. We continued to use the excitation signal as shown in Fig. 12, exciting at a frequency of 35 kHz, and made efforts to ensure that the excitation point was the same as the selected excitation point in the simulation experiments. At each receiving point, three groups of data are measured, and the signal strength of the receiving point is determined by calculating the average peak energy value across the three groups of data. The data are shown in Table 5. Through experimental measurement, the amplitude of the guided wave signal can still reach more than 4.7mV after 1km propagation. Table 5 Data from the field experiments.

Distance (m)	Energy (V2)	
50	14.23606	
100	8.190338	
150	7.03298	
200	3.220529	
250	3.274978	
300	2.600834	
350	1.776809	
400	1.079978	
450	0.532596	
500	0.372291	
550	0.040284	
600	0.003655	
650	0.003576	
700	0.002464	
750	0.001573	
800	0.001172	
850	0.000266	
900	5.15e − 05	
950	6.11e − 05	
1000	2.24e − 05	

The data in Table 5 shows that the attenuation curve is fitted exponentially. The fitting effect is shown in Fig. 26, and the energy attenuation coefficient is calculated to be about 0.0078.Figure 26 Attenuation curve on the Beijing Circular Railway.

The attenuation experiment shows that the mode selected in this study can receive an effective signal at 1km when the excitation signal amplitude is ± 125V, indicating that this mode can be used for long-distance rail flaw detection.

To explore the influence of environmental factors on the attenuation of guided wave signals, we conducted additional experiments on two different railway tracks. One experiment was conducted on a straight track, while the other was conducted on a track with smaller radius curves. In these experiments, we used higher-precision three-axis accelerometer sensors and increased the excitation voltage to ± 300V. After excitation, signals were collected at different distances. The experimental data on a straight track were shown in Table 6 and the attenuation curves of the second experiment were plotted as shown in Fig. 27. Table 6 Experimental data on a straight track.

Distance (m)	Eheamplitude spectrum (dB)	
5	117.3	
100	99.99	
866	93.27	
2500	68.53	
3300	37.92	

Figure 27 Attenuation curve with curve section.

From Table 6, it can be observed that on the straight track segment with good operating conditions, the guided wave modes excited by this method exhibit smaller attenuation. Guided wave signals can still be detected at a distance of 3 km.

In Fig. 27, we can see that on curved track segments, the attenuation of guided wave signals increases. This indicates that the method proposed in this study is still influenced by environmental factors such as curves. To address this issue, we can increase sensor density in curved sections according to the actual condition of the track.

Conclusion

In recent years, ultrasonic guided wave (UGW) has become an essential tool for non-destructive testing of railways. Aiming at the application of UGW in long-distance rail breakage monitoring of U78CrV rail, the main contributions of this paper are as follows.The dispersion curve of UGW in the U78CrV rail is solved by the semi-analytical finite element method, and the guided wave mode suitable for monitoring the whole section monitoring of this rail is selected. The excitation point and mode are determined by calculating the vibration energy of this mode on the rail surface.

The three-dimensional finite element simulation of ANSYS is used to analyze the rail vibration, and the group velocity of the vibration energy concentration point is solved by wavelet transform. It is verified that the selected mode is the primary mode existing in the rail.

A specific rail section is selected, the arrival time of the selected mode is calculated, and the vibration mode shape data near this time is extracted and compared with the mode shape of the semi-analytical mode to verify the successful excitation of the selected mode.

Through simulation and field experiment, it is verified that the guided wave mode excited in this paper is sensitive to the cracks at the rail head, rail web, and rail foot and can travel more than 1km in heavy haul railway rails. Remote online monitoring of rail defects or broken rails can be realized by detecting the reflected echo of defects or attenuating the guided wave signals after defects.

In conclusion, through experimental validation, the ultrasonic guided wave detection method proposed in this study achieves a detection distance of over 3 km on straight sections, demonstrating promising prospects in long-distance rail break detection. An important development direction for future rail integrity health monitoring is ultrasonic guided wave-based long-distance rail break detection. It involves selecting and exciting modes with minimal influence from track conditions and environmental factors and low attenuation rates by deeply studying the propagation mechanism of guided waves.

Ethical approval

This article does not contain any studies with human participants or animals performed by any authors.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-59328-5.

Author contributions

X.X. supervised the experimental program, manuscript preparation and improved the manuscript. Z.W. and Y.N. wrote the main body of the manuscript, and B.S and X.M. carried out the experimental part, and Z.P. carried out the simulation work. All the authors reviewed the manuscript".

Funding

This study was supported by three grants: the Fundamental Research Funds for the Central Universities of China (Grant Number 2022JBQY002) and the Science, technology innovation project of Shuohuang Railway Development Co., Ltd. under China Energy (Grant Number GJNY-21-65) and State Key Laboratory of Advanced Rail Autonomous Operation (Contract No.RAO2023ZZ003).

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Song ZY Bian LP Research on Intelligent heavy haul railway system framework Railw. Constr. 2022 62 01 4 10
Song, Z. Y. & Bian, L. P. Research on Intelligent heavy haul railway system framework. Railw. Constr. 62(01), 4–10 (2022).
2. Li RJ Cause analysis and preventive measures of track breaking of Shenshuo heavy haul railway Shenhua Technol. 2011 9 04 81 84
Li, R. J. Cause analysis and preventive measures of track breaking of Shenshuo heavy haul railway. Shenhua Technol. 9(04), 81–84 (2011).
3. Shi HZ Ren Y Development status and research of track breaking detection technology at home and abroad Railw. Oper. Technol. 2010 16 4 1 7
Shi, H. Z. & Ren, Y. Development status and research of track breaking detection technology at home and abroad. Railw. Oper. Technol. 16(4), 1–7 (2010).
4. Chahbaz A Development of a mobile inspection system for rail Integrity assessment Transp. Dev. Centre Saf. Secur. Transp. Canada 2000 66 1 31
Chahbaz, A. Development of a mobile inspection system for rail Integrity assessment. Transp. Dev. Centre Saf. Secur. Transp. Canada 66, 1–31 (2000).
5. Vitaly, I., & Alexander, V. Development of instruments for NDT of rail-tracks with use of contactless electromagnetic acoustic emission transducers (EMATs). In WCNDT, Shanghai, China 255–262 (2008).
6. Papaelias MP Roberts C A review on non-destructive evaluation of rails: State-of-the-art and future development Proc. IMechE 2008 222 4 367 384 10.1243/09544097JRRT209
Papaelias, M. P. & Roberts, C. A review on non-destructive evaluation of rails: State-of-the-art and future development. Proc. IMechE 222(4), 367–384 (2008).10.1243/09544097JRRT209
7. Clark R Rail flaw detection: Overview and needs for future developments NDT&E Int. 2004 37 2 111 118 10.1016/j.ndteint.2003.06.002
Clark, R. Rail flaw detection: Overview and needs for future developments. NDT&E Int. 37(2), 111–118 (2004).10.1016/j.ndteint.2003.06.002
8. Du QL Zhao DH Research on real-time track breaking detection based on non-track circuit method Sci. Technol. Wind 2020 35 91 92
Du, Q. L. & Zhao, D. H. Research on real-time track breaking detection based on non-track circuit method. Sci. Technol. Wind 35, 91–92 (2020).
9. Zhao C Wang P Detection method for broken rail based on rate of change of strain mode J. Vib. Mea. Diagnos. 2012 32 5 723 729
Zhao, C. & Wang, P. Detection method for broken rail based on rate of change of strain mode. J. Vib. Mea. Diagnos. 32(5), 723–729 (2012).
10. Deng HL Lu Q Composite insulator detection based on ultrasonic guided wave High Volt. Tech. 2016 42 4 1236 1244
Deng, H. L. & Lu, Q. Composite insulator detection based on ultrasonic guided wave. High Volt. Tech. 42(4), 1236–1244 (2016).
11. Bi, Y. S. Research on Broken Track State of Track Circuit of High-Speed Railway (Ph.D. BJTU, Beijing, 2015).
12. Tian MX Chen YF Overview of real-time rail breaking detection methods J. Lanzhou Jiaotong Univ. 2011 30 01 122 126
Tian, M. X. & Chen, Y. F. Overview of real-time rail breaking detection methods. J. Lanzhou Jiaotong Univ. 30(01), 122–126 (2011).
13. He CF Zheng MF Development, applications and challenges in ultrasonic guided waves testing technology Chin. J. Sci. Instrum. 2016 37 8 66
He, C. F. & Zheng, M. F. Development, applications and challenges in ultrasonic guided waves testing technology. Chin. J. Sci. Instrum. 37(8), 66 (2016).
14. Zhu LQ Xu Xining XN Study on the railway integrity monitoring method based onultrasonic guided waves Chin. J. Sci. Instrum. 2016 37 7 1603 1609
Zhu, L. Q. & Xu Xining, X. N. Study on the railway integrity monitoring method based onultrasonic guided waves. Chin. J. Sci. Instrum. 37(7), 1603–1609 (2016).
15. Lamb H On waves in an elastic plate Proc. R. Soc. Lond. 1917 93 648 114 128
Lamb, H. On waves in an elastic plate. Proc. R. Soc. Lond. 93(648), 114–128 (1917).
16. Thomson WT Transmission of elastic waves through a stratified solid medium J. Appl. Phys. 1950 21 2 89 93 10.1063/1.1699629
Thomson, W. T. Transmission of elastic waves through a stratified solid medium. J. Appl. Phys. 21(2), 89–93 (1950).10.1063/1.1699629
17. Knopoff L A matrix method for elastic wave problems Bull. Seismol. Soc. Am. 1964 54 1 431 438 10.1785/BSSA0540010431
Knopoff, L. A matrix method for elastic wave problems. Bull. Seismol. Soc. Am. 54(1), 431–438 (1964).10.1785/BSSA0540010431
18. Dunkin JW Computation of modal solutions in layered, elastic media at high frequencies Bull. Seismol. Soc. Am. 1965 55 2 335 358 10.1785/BSSA0550020335
Dunkin, J. W. Computation of modal solutions in layered, elastic media at high frequencies. Bull. Seismol. Soc. Am. 55(2), 335–358 (1965).10.1785/BSSA0550020335
19. Pavlakovic B Lowe M Disperse: A general purpose program for creating dispersion curves Rev. Prog. Quant. Nondestruct. Eval. 1997 66 185 192 10.1007/978-1-4615-5947-4_24
Pavlakovic, B. & Lowe, M. Disperse: A general purpose program for creating dispersion curves. Rev. Prog. Quant. Nondestruct. Eval. 66, 185–192 (1997).10.1007/978-1-4615-5947-4_24
20. Xing, B. Research on Rail Crack Detection Method Based on Ultrasonic Guided Wave (M.S., BJTU, Beijing, 2020).
21. Mukdadi OM Desai YM Elastic guided waves in a layered plate with rectangular cross-section J. Acoust. Soc. Am. 2002 112 5 1766 1779 10.1121/1.1509426 12430791
Mukdadi, O. M. & Desai, Y. M. Elastic guided waves in a layered plate with rectangular cross-section. J. Acoust. Soc. Am. 112(5), 1766–1779 (2002).12430791 10.1121/1.1509426
22. Hayashi T Song W Rose JL Guided wave dispersion curves for a bar with an arbitrary cross-section, a rod and rail example Ultrasonics 2003 41 3 175 183 10.1016/S0041-624X(03)00097-0 12726938
Hayashi, T., Song, W. & Rose, J. L. Guided wave dispersion curves for a bar with an arbitrary cross-section, a rod and rail example. Ultrasonics 41(3), 175–183 (2003).12726938 10.1016/S0041-624X(03)00097-0
23. Ming C Wu XJ Dispersion analysis of guided waves in the finned tube using the semi-analytical finite element method J. Sound Vib. 2017 401 114 126 10.1016/j.jsv.2017.04.037
Ming, C. & Wu, X. J. Dispersion analysis of guided waves in the finned tube using the semi-analytical finite element method. J. Sound Vib. 401, 114–126 (2017).10.1016/j.jsv.2017.04.037
24. Duan WB Ray K Guided wave propagation in buried and immersed fluid-filled pipes: Application of the semi analytic finite element method Comput. Struct. 2019 212 236 247 10.1016/j.compstruc.2018.10.020
Duan, W. B. & Ray, K. Guided wave propagation in buried and immersed fluid-filled pipes: Application of the semi analytic finite element method. Comput. Struct. 212, 236–247 (2019).10.1016/j.compstruc.2018.10.020
25. Duan WB Gan TH Investigation of guided wave properties of anisotropic composite laminates using a semi-analytical finite element method Compos. Part B Eng. 2019 173 66 10.1016/j.compositesb.2019.106898
Duan, W. B. & Gan, T. H. Investigation of guided wave properties of anisotropic composite laminates using a semi-analytical finite element method. Compos. Part B Eng. 173, 66 (2019).10.1016/j.compositesb.2019.106898
26. Lee CM Rose JL A guided wave approach to defect detection under shelling in rail NDT&E Int. 2009 42 174 180 10.1016/j.ndteint.2008.09.013
Lee, C. M. & Rose, J. L. A guided wave approach to defect detection under shelling in rail. NDT&E Int. 42, 174–180 (2009).10.1016/j.ndteint.2008.09.013
27. Coccia S Bartoli I Alessandro M Numerical and experimental study of guided waves for detection of defects in the rail head NDT&E Int. 2011 44 1 93 100 10.1016/j.ndteint.2010.09.011
Coccia, S., Bartoli, I. & Alessandro, M. Numerical and experimental study of guided waves for detection of defects in the rail head. NDT&E Int. 44(1), 93–100 (2011).10.1016/j.ndteint.2010.09.011
28. Xu XN Guo BQ Solution of ultrasonic guided wave dispersion curve in rail by semi analytical finite element method Chin. J. Sci. Instrum. 2014 35 10 2392 2398
Xu, X. N. & Guo, B. Q. Solution of ultrasonic guided wave dispersion curve in rail by semi analytical finite element method. Chin. J. Sci. Instrum. 35(10), 2392–2398 (2014).
29. Loveday PW Long CS Laser vibrometer measurement of guided wave modes in rail track Ultrasonics 2015 57 209 217 10.1016/j.ultras.2014.11.010 25497003
Loveday, P. W. & Long, C. S. Laser vibrometer measurement of guided wave modes in rail track. Ultrasonics 57, 209–217 (2015).25497003 10.1016/j.ultras.2014.11.010
30. Setshedi II Loveday PW Estimation of rail properties using semi-analytical finite element models and guided wave ultrasound measurements Ultrasonics 2019 96 240 252 10.1016/j.ultras.2018.12.015 30685066
Setshedi, I. I. & Loveday, P. W. Estimation of rail properties using semi-analytical finite element models and guided wave ultrasound measurements. Ultrasonics 96, 240–252 (2019).30685066 10.1016/j.ultras.2018.12.015
31. Dineo AR Daniel NW Development of an optimal piezoelectric transducer to excite guided waves in a rail waist NDT&E Int. 2018 95 72 81 10.1016/j.ndteint.2018.02.002
Dineo, A. R. & Daniel, N. W. Development of an optimal piezoelectric transducer to excite guided waves in a rail waist. NDT&E Int. 95, 72–81 (2018).10.1016/j.ndteint.2018.02.002
32. Ramatlo DA Long CS A modelling framework for simulation of ultrasonic guided wave-based inspection of welded rail tracks Ultrasonics 2020 108 66 10.1016/j.ultras.2020.106215
Ramatlo, D. A. & Long, C. S. A modelling framework for simulation of ultrasonic guided wave-based inspection of welded rail tracks. Ultrasonics 108, 66 (2020).10.1016/j.ultras.2020.106215
33. Hua GF Zhu L Blockchain-based federal learning for intelligent control in heavy haul railway IEEE Access 2020 99 1
Hua, G. F. & Zhu, L. Blockchain-based federal learning for intelligent control in heavy haul railway. IEEE Access 99, 1 (2020).
34. Mahajan H Banerjee S A machine learning framework for guided wave-based damage detection of rail head using surface-bonded piezo-electric wafer transducers Mach. Learn. Appl. 2022 7 66
Mahajan, H. & Banerjee, S. A machine learning framework for guided wave-based damage detection of rail head using surface-bonded piezo-electric wafer transducers. Mach. Learn. Appl. 7, 66 (2022).
35. Sresakoolchai J Kaewunruen S Railway defect detection based on track geometry using supervised and unsupervised machine learning Struct. Health Monit. 2022 21 4 1757 1767 10.1177/14759217211044492
Sresakoolchai, J. & Kaewunruen, S. Railway defect detection based on track geometry using supervised and unsupervised machine learning. Struct. Health Monit. 21(4), 1757–1767 (2022).10.1177/14759217211044492
36. Li YF Liu FS Cause analysis of corrosion and fracture of u78crv rail on heavy vehicle line of Daqin Railway Railw. Constr. 2019 59 03 116 119
Li, Y. F. & Liu, F. S. Cause analysis of corrosion and fracture of u78crv rail on heavy vehicle line of Daqin Railway. Railw. Constr. 59(03), 116–119 (2019).
37. Xu XN Bo X A graphical analysis method of guided wave modes in rails Appl. Sci. 2019 9 8 66 10.3390/app9081529
Xu, X. N. & Bo, X. A graphical analysis method of guided wave modes in rails. Appl. Sci. 9(8), 66 (2019).10.3390/app9081529
38. Xu, X. N. Research on the Application of Ultrasonic Guided Wave on-Line Stress Monitoring Technology for Seamless Rail 43–48 (Ph.D, BJTU, Beijing, 2020).
39. Bartoli I Marzani A Modeling wave propagation in damped waveguides of arbitrary cross-section Proc. SPIE Int. Soc. Opt. Eng. 2006 295 3–5 685 707
Bartoli, I. & Marzani, A. Modeling wave propagation in damped waveguides of arbitrary cross-section. Proc. SPIE Int. Soc. Opt. Eng. 295(3–5), 685–707 (2006).
40. Rose JL Avioli MJ Guided wave inspection potential of defects in rail NDT&E Int. 2004 37 153 161 10.1016/j.ndteint.2003.04.001
Rose, J. L. & Avioli, M. J. Guided wave inspection potential of defects in rail. NDT&E Int. 37, 153–161 (2004).10.1016/j.ndteint.2003.04.001
