Abstract
In order to accurately analyze the influence of electromagnetic transient signals on the jointless track circuit when the electromagnetic transient signal propagates in the rail, it is necessary to consider the frequency-variable load terminated in the ZPW-2000A jointless track circuit and the frequency-variable loss inside the rail. A method is proposed for calculating the transient response of transmission lines system with frequency-variable end load of jointless track circuit. Firstly, the transmission lines model of jointless track circuit is established, based on multiconductor transmission lines theory, the model equation is deduced and discretized by finite difference time domain (FDTD). The vector fitting method is used to express the admittance of the tuning region in the track circuit, and the rational approximation function of the tuning region is derived from the poles, residues, and constants. The voltage and current at nodes in the tuning region are calculated by piecewise linear recursive convolution algorithm. Combined with the discrete transmission line equation, the current and voltage expression of the transient electromagnetic signal at the receiving end of the track circuit in time domain is obtained. Compared with state variable method, the error is less than 6%, which verifies the correctness of the proposed method. Finally, this paper studies the influence laws of different factors on the overvoltage at the receiving end of jointless track circuit and the weak links of jointless track circuit under the influence of transient electromagnetic signal. It provides theoretical reference for fault research and anti-interference analysis of ZPW-2000A jointless track circuit.
1. Introduction
Jointless track circuit is the key equipment used to detect the occupancy of train section, the integrity of rail, and the transmission of running information in railway line, which is an important basic equipment to ensure the safe operation of high-speed railway [1]. With the increasing complexity of railway construction environment, the rail and signal transmission equipment laid outdoors are easily disturbed by strong electromagnetic signals. Electromagnetic interference signals use rail as transmission medium to produce conduction interference or radiation coupled interference on jointless track circuits [2]. The investigations show that lightning strikes can cause failures of trackside tuning matching units, signal transmission cables and lightning protection analog network disks, transmitters and receivers, etc. Jointless track circuit is easy to make wrong responses by this disturbance, which will threaten the safety of train running. Therefore, it is important to analyze the interference of transient electromagnetic signal to track circuit.
The theory of multiconductor transmission lines is the theoretical basis of time-domain transient analysis of track circuits. Kunz and Luebbers [3] first solved the transmission lines equation by FDTD method, which provided a theoretical basis for time-domain solution of transmission lines. With the complex changes in transmission lines termination equipment, Wang et al. [4–6] respectively proposed a numerical algorithm based on Thevenin equivalent circuit method, state variable method, etc., to solve the terminal frequency-varying load and combined with the transmission line equation to solve the electromagnetic impulse response of the system. In the study of transient response of jointless track circuits. Mazloom et al. [7] proposed a method based on transmission line equation and ATP-EMTP simulation model to analyze the voltage response of track circuit at the moment of train entering and clearing. Zhao et al. [8] obtained the general solution of current and voltage of rail line by Laplace transformation of the transmission line equation of track circuit, which simplified the calculation model and ignored the electrical insulation section; Yong-jian [9] used the state variable method to analyze the receiving end of the DC track circuit and solved equation of the track circuit by combining the transmission line equation, but it was difficult to write the state equation when the circuit structure was complex. Zhichao et al. [10] analyzed the overvoltage of each signal equipment during lightning strike catenary by establishing lightning strike simulation model of signal system.
ZPW-2000A jointless rack circuit includes parallel nonlinear equipment and rail with frequency variation characteristics. In order to accurately calculate the electromagnetic transient response of track circuit transmission line system, the influence of terminating frequency variation load and rail frequency variation loss must be considered. Firstly, according to the transmission characteristics of ZPW-2000A jointless track circuit, the transmission line model of track circuit is established. The model cannot be calculated directly due to the presence of parallel or terminated nonlinear lumped elements in the track circuit. Based on the multiconductor transmission lines theory, the model equation is solved by FDTD. The vector fitting method is a powerful tool for the modeling analysis of frequency-dependent systems, which is used to fit the admittance of the tuning region. The rational approximation function of the track circuit is derived from the poles, residues, and constants. The relationship between voltage and current at both ends of the admittance is derived by piecewise linear recursive convolution, which is submitted into the discrete equation of transmission lines. Finally, the time-domain expression of the current and voltage at the receiving end of the track circuit is obtained. Compared with the state variable method, the error is less than 6%, which verifies the correctness of the proposed method. This method is used to analyze the influence of lightning stroke distance, electromagnetic signal frequency, and track bed resistance on the overvoltage at the receiving end of track circuit; it provides a theoretical basis for the anti-interference analysis of track circuit.
2. Modeling of Transmission Line of Track Circuit with Frequency-Variable Load Terminated
The outdoor equipment of jointless track circuit is composed of steel rail, tuning unit (BA), air-core coil (SVA), signal cable, matching unit and compensation capacitor, etc. The tuning unit and the air-core coil part form an electrical insulation section to prevent cross-zone signal transmission from adjacent sections. The compensation capacitor is used to improve the transmission performance of the track circuit and increase the transmission distance of the track circuit information. According to the transmission characteristics of track circuit, the rails are equivalent to uniform transmission lines [11]; it is characterized by two asymmetrical currents, one of which leaks into the ground, and the other flows from one rail to the other through ballast and sleepers. Taking a certain interval of track circuit as an example, the MTL equivalent model of jointless track circuit is established [12], as shown in Figure 1.

Here, r1 and r2 are the self-resistance of the rail; l11 and l22 are the self-inductance of the rail; l12 is the mutual inductance between rails; and are the leakage conductance between rails and Earth; and are the ballast resistance parameter of track characterized by conductance; c11 and c22 are the self-capacitance of the rail; c12 is the mutual capacity of the rails; and z1 and z2 represent the tuning region, respectively.
Transmission line losses are reflected in the conductor and the surrounding mediums; in general, these losses are frequency dependent, and the increase of frequency will lead to additional high-frequency losses. The most important is caused by the skin effect in the conductor and the polarization in the dielectric. Dielectric loss has little influence on transmission lines and can be ignored [12]. Therefore, the frequency domain expression of transmission line equation based on multiconductor transmission lines theory iswhere impedance and admittance are, respectively, , . .
Converting the frequency domain equations of (1) and (2) is transformed by Laplace transform into the time domain given the convolution equation as
In the case of high frequency, unit length resistance, internal inductance, and unit length conductance of conductor loss parameters are functions of frequency. The approximate formula for conductor losses is based on the transformation pair conversion; the inverse Laplace transform of the impedance in the conductor iswhere , , rdc is DC resistance.
The convolution of the calculation in (3) is expressed as
The convolution term (6) is discrete with time step and divided by Δt segment. Discrete convolution is approximated in the following way:where , .
According to the FDTD algorithm, the position length of the transmission line is discretized by Δz, and the time length is discretized by Δt. The whole transmission lines are divided into NDZ segments, and the total solution times are divided into NDT segments. In order to ensure the stability of discrete calculation, the voltage of NDZ + 1 at discrete point and the discrete point current I1,I2, …, INDT does interleaving calculation. Then, (3), (4), and (7) can be discretized as
The recursive relation of current at discrete points inside the transmission line is
3. Tuning Region Transient Modeling
In the ZPW-2000A jointless track circuit, the electrical insulation section realizes the electrical isolation of adjacent track and the stable output of local signal by forming series and parallel resonance for different carrier frequency signals. In the matching part, the rail impedance and cable impedance are matched to realize the output of high power signal to the rail [14]. The structure of electrical insulation section in jointless track circuit is shown in Figure 2.

Firstly, the transient model of tuning matching unit is established, and the drive point admittance of the tuning unit port is obtained by the short-circuit test method. The frequency domain admittance is fitted to a rational function expression by the vector matching method, as shown in (10). Secondly, the frequency domain rational function is transformed into time-domain expression by inverse Fourier transform. Finally, based on the piecewise linear recursive convolution method, the recursive iterative relation of admittance current is derived in (14). The admittance can be fitted to a rational function in the complex frequency domain by the vector matching method and the objective function Y(s) expresses its s-domain admittance in pole and residue form [13]:where n is the order, i is an integer (i = 1,2, …, n); ai is the residue, and ci is the pole, both are real numbers or conjugate complex number pairs, sh is the coefficient of the first term, and is the constant term.
The fitting order of the poles in the tuning area is 12 to achieve appropriate fitting effect. The amplitude characteristic and phase angle characteristic function are fitted according to the pole and residue values. Compared with the admittance model of the tuning unit, the results are shown in Figures 3 and 4.


The results of admittance fitting of the tuning element are consistent with those of the exact value and meet the requirements of calculation. According to the admittance rational function, the current equation of the connection node is derived. The voltage current relationship of the tuning unit in time domain is expressed as . Defining two parameter variables [5],
Using piecewise recursive convolution integral, we know that if the two variables are satisfied,
Then, the time-domain current can be written in the recursive iterative form:
There are number of real numbers Nr and conjugate pairs Ng in the fitting residues of loads in residue regions and the poles, then the total current of the connected nodes iswhere , ;
4. Transmission Line Voltage Equation of Track Circuit
When the transmission line load is pure resistance load, the current is . The central difference method is used to discretize it.
According to the average value of power supply IS current, (3) and (4) can be discretized:
The initial voltage expression of the transmission line equation of the track circuit is
The internal point voltage iteration formula of the equation is
According to the average value of load current IL at the load at the end of the transmission line, the transmission line equation can be discretized as follows:
When the nonlinear lumped element of a track circuit is terminated with a tuning unit as the terminal condition, the terminal condition gives the relation IL between voltage and current at the same time and position. According to (14), the current and voltage equation IL at its endpoints is
The terminal voltage equation of the transmission line equation of the track circuit is
5. Method Validation
To verify the validity of the constructed track circuit model and the transient response formula, taking P60 rail carrier frequency signal 2300 Hz and 1.4 km track section as an example, the lightning wave adopted the double exponential lightning electromagnetic signal with wave head of and wavelength of , recorded as 1.2/50 µs standard lightning wave [15], and the amplitude was 30 kA for simulation. According to the transmission line model of track circuit, the overvoltage of track surface at the receiving end of track circuit is solved. The state variable method is used to compare the response results with the proposed method.
It can be seen from Figures 5 and 6 that the numerical algorithm for transient response of track circuit proposed in this paper is consistent with the calculation results of state variable method. The calculation error is less than 6%.


5.1. Impact Tolerance Level of Signal Equipment
When the electromagnetic field induced by lightning reaches a certain level, the signal equipment will work improperly or be damaged. The lightning tolerance level of signal equipment terminals is obtained through the destructive lightning shock tolerance experiment, which is the key factor to determine the lightning resistance of signal system equipment [10], as shown in Table 1.
6. Analysis of Influencing Factors
6.1. Influence of Lightning Transmission Distance on Rail Surface Voltage
When lightning strikes the rail, the change of track surface overvoltage with distance between the lightning point and the tuning area of the track circuit is shown in Figure 7.

It can be seen from Figure 7 that with the increase of the distance between the lightning point and the receiving end of the track circuit, the rail surface overvoltage decreases rapidly with the increase of the distance. Especially when the transmission distance is less than 700 m, the overvoltage generated by lightning will exceed the withstand voltage range of the tuning unit end, and there is a risk of breakdown of trackside electronic equipment.
6.2. Considering the Influence of Different Frequencies on Rail Surface Voltage
In ZPW-2000A jointless track circuit, the rail is used as the transmission channel of frequency shift signal, and the rail parameters have the characteristics of frequency change. In the high-frequency electromagnetic environment, the load at the end of the transmission line usually shows the frequency change effect. In order to analyze the variation rules of signal overvoltage at the receiving end of track circuit when interference signals of different frequencies act on track circuit, the calculation results are shown in Figure 8.

In Figure 8, the simulation results show that the frequency of electromagnetic interference signal has a great influence on the voltage amplitude at the receiving end of the rail surface. At the same time, the higher the frequency of the interference signal is, the smaller the amplitude of the rail surface voltage is at the receiving end. With the increase of signal frequency, the rail impedance increases, the loss of interference signal on the track circuit increases, and the rail voltage amplitude decreases at the receiving end. The internal frequency loss of the transmission line is very little affected by the change of frequency and can be ignored in the future research.
6.3. Influence of Track Bed Resistance on Rail Surface Voltage
The primary side parameters of the jointless track circuit will affect the transmission of interference signals on the rail. The P60 rail is selected, its impedance value is a certain value, while the resistance of the track bed will change with the difference of ambient temperature and humidity [16]. In order to study the influence rule of track bed resistance on rail surface voltage at the receiving end of track circuit, when the resistance value of track bed is 2 Ω·km, 2.5 Ω·km, 4 Ω·km, and 5 Ω·km, the influence rule of different track bed resistance on rail surface overvoltage is analyzed, as shown in Figure 9.

As shown in Figure 8, the greater the ballast resistance is, the greater the overvoltage is. Meanwhile, the faster the rate of rise and attenuation of overvoltage are. With the decrease of track ballast resistance, the amplitude of rail surface voltage decreases at the receiving end.
7. Conclusion
(1)The transient response model of track circuit is established. A numerical calculation method of transient electromagnetic response of track circuit transmission line system considering terminated frequency-varying load and frequency-varying loss is proposed, which is suitable for the analysis of transmission line model with terminated frequency-varying load. It provides theoretical reference for the transmission performance of jointless track circuits and the anti-interference design of track circuits.(2)The frequency of electromagnetic interference signal has great influence on the amplitude of rail voltage at the receiving end. However, the frequency has little effect on internal loss of rail. The overvoltage and the time of signal attenuation to stability increase with increasing of the ballast resistance.(3)According to withstand voltage level of the signal equipment, when the distance between the lightning strike point and the receiving end of the track circuit is small, there is a risk of insulation breakdown on the rail side of the tuning unit, so measures need to be taken to strengthen the protection of the signal equipment.Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work has been financially supported by the Natural Science Foundation of China (No. 5196070019), the Natural Science Foundation of Gansu Province (No. 21JR7RA292), and the Youth Scientific Research Fund of Lanzhou Jiaotong University (No. 1200060833).