Train-Bridge Dynamic Behaviour of Long-Span Asymmetrical-Stiffness Cable-Stayed Bridge

Long-span cable-stayed bridge (LCB) with unequal-height towers is being designed and constructed inmetro lines due to its better adaptability to environment and terrain conditions compared to traditional cable-stayed bridge with equal-height towers. However, the asymmetrical arrangement of towers leads to obvious nonuniformity of the structural stiffness along the longitudinal direction, which intensifies the wheel-rail coupled vibration behaviour, and affects the running safety of operating trains and ride comfort. 0erefore, train-bridge dynamic behaviour of long-span asymmetrical-stiffness cable-stayed bridge is deeply investigated in this work. Primarily, considering the comprehensive index of frequency difference and modal assurance criterion (MAC), a nonlinear model updating technique (NMUT) based on penalty function theory is proposed, which can be used to optimize the bridge numerical model. Secondly, on the basis of the train-track-bridge dynamic interaction theory (TDIT), a traintrack-bridge coupled dynamic model (TCDM) is established. Finally, a LCB with unequal-height towers is applied as a case to illustrate the influence of asymmetrical stiffness on the train-track-bridge dynamic characteristics. Results show that the proposed NMUT is efficacious and practical. For the LCB with unequal-height towers, a significant difference between the bridge vibration at low tower location and that at high tower location appears.0e vertical displacement difference of the main beam on both sides of the bridge increases with the distance from the observation point to the bridge tower increasing. 0e variation of acceleration difference on both sides of the bridge is influenced by the speed of the train and the position of the observation point simultaneously. In general, vibrations of the main beam at low tower location are larger than those at high tower location.


Introduction
During the past few years, due to the dynamic development of China's urban long-span railway bridges, cable-stayed bridges are widely built because of their unique and beautiful shape, small weight of bridge structure, excellent spanning ability, and perfect construction technology. Compared with the traditional cable-stayed bridge with equal-height towers, the cable-stayed bridge with unequal-height towers (UHTCB) has better adaptability to the environment and terrain conditions. e UHTCB is an asymmetrical statically indeterminate structure. is asymmetrical arrangement of towers leads to great nonuniformity of the stiffness along the longitudinal direction, which results in additional stiffness irregularity [1]. When the train passes through the UHTCB, the rail geometric irregularities and additional stiffness irregularities worsen the wheel-rail nonlinear contact state and then intensify the wheel-rail coupled vibration behaviour, which affects the ride comfort and running safety of operating trains. Moreover, the increased wheel-rail dynamic interaction (WDI) further amplifies the vibration of the cablestayed bridge system, resulting in the great difference in vibration, deformation, and stress states between the bridge at low tower location and that at high tower location, which eventually reduces the overall stability of the structure and affects the normal and safe service of the bridge. So it is very essential to investigate the vibration of the UHTCB.
Train-bridge coupled vibration (TCV) is a complicated interaction process, which has been studied for a very long time since it is a typical topic of railway dynamics [2,3]. According to Zhai et al., the train, the track, and the bridge served as a whole dynamic system, with the purpose of offering an approach to analyse and assess the running safety of high-speed trains and the ride comfort. As a result, Zhai et al. established an integrated train-track-bridge dynamic interaction model and developed the software TTBSIM. Moreover, lots of field tests were carried out for validating this software [4,5]. Olmos et al. [6] developed an efficient nonlinear train-track-bridge dynamic interaction model to analyse the transverse dynamic response of high pier viaducts under the action of high-speed train. e influences of railway car-body flexibility on the dynamic performance of high-speed trains operating on bridges were investigated by Younesian et al. who found that the calculated ride comfort index can be influenced by the flexural mode shapes of the car-body structure dramatically particularly in the low frequency range [7]. Under the resonant excitations, an alternative analytical method was proposed by Yau et al. on the basis of an equivalent VBI model, providing insights of accounting for VBI impacts on short-span simply supported bridges in the railway by enabling a precise prediction of the additional amount of damping [8]. e models considering nonlinear wheel-rail contact forces were developed by Antolín et al. to research dynamic influences in the vehicles and in the bridge caused by the coupling and to analyse the dynamic interaction between bridges and high-speed trains [9]. In consideration of the vehicle-track-structure interaction, Melo et al. calibrated and validated a numerical model of a steel-concrete composite bridge [10]. Moreover, references [11][12][13][14][15][16][17][18] also investigated the train-track-bridge coupled interaction and proposed many beneficial conclusions. As known from existing studies, most works are focused on simple bridge structures such as continuous beam bridges or simply supported beam bridges, few researches on TCV involving complex bridges such as cable-stayed bridges, and even less on the UHTCB.
In the TCV simulation, the accuracy of the bridge FE model has an important influence on the accuracy of the simulation results. However, the FE model is different from the actual structure due to the influence of mesh generation, boundary conditions, and uncertainty of structural physical parameters [19]. erefore, the FE model updating should be carried out to make it as close as possible to the actual structure, so as to conduct the vibration characteristics analysis of train-bridge system more accurately. Scholars worldwide have done a lot of researches on the model optimization and updating. Ren et al. [20] put forward an FE model updating technique in virtue of a response surface method in structural dynamics for civil engineering structures. Liu et al. [21] presented a method in virtue of fuzzy FE to take into account the impact of the uncertainties of the measured modal parameters on the updated FE model. Chakraborty et al. [22] made an attempt to explore the effectiveness of moving least-squares technique based on response surface method to update the FE model. Asgarieh et al. [23] proposed updating parameters of hysteretic material models assigned to substructures of a nonlinear FE model by minimizing an OF. Astroza et al. [24] advanced a hybrid global optimization algorithm by combining the unscented Kalman filter with the simulated annealing. Wang et al. [25] proposed a model updating tactic for localized nonlinear structures. ese algorithms greatly improve the accuracy of the structural numerical model, but they are rarely used in long-span bridges, especially in the UHTCB.
Aiming at the complex train-bridge spatial coupled vibration problem caused by asymmetrical stiffness of the UHTCB, an NMUT is proposed, and a refined TCDM is established on the basis of the TDIT. A long-span UHTCB is applied as a case to illustrate the influence of asymmetrical stiffness on the train-bridge system dynamic behaviour.

Model Optimization Based on NMUT
When studying train-bridge dynamic behaviour of the LCB, remarkable discrepancies in respect of the mechanical and dynamic characteristics between the numerical model and the actual structure can be in existence because of the influence of mesh division, boundary conditions, and uncertainties related to physical parameters and assumptions. erefore, it is considered very necessary to improve the accuracy of the FE model. Before conducting the study on the train-bridge dynamic characteristics of the LCB, the natural vibration characteristics of the LCB are analysed first. By updating the parameters of the FE model, the natural vibration characteristics of the model are as close as possible to the actual structure, to ensure that the subsequent study of the dynamic characteristics of the train and bridge is more meaningful. To calculate and analyse the TCV more accurately, considering the comprehensive index of frequency difference and MAC, an NMUT based on penalty function theory is proposed to optimize the bridge numerical model. e flow chart of updating the bridge model is shown in Figure 1.

Establishment of Initial Bridge Model.
In virtue of the FEM, the bridge model is established, and dynamic equation can be expressed by where M b /C b /K b are the mass/damping/stiffness matrices of the bridge system; X b / _ X b / € X b are the displacement/velocity/ acceleration vectors of the bridge system; F b is the force vector for the bridge system.
According to (1), the free vibration motion equation of the undamped structural system is established: Under specific initial conditions, the simple harmonic vibration equation of the system at the same frequency is established: Taking the derivative of (3) and substitute it into (2) and noticing the arbitrariness of sin(ωt + θ), there is 2 Shock and Vibration en, the characteristic equation of the system is obtained as e root ω 2 i (i �1, 2, . . ., N) of (5) is called the ith eigenvalue. ω i (i �1, 2, . . ., N) is called the ith order natural circle frequency (rad/s). If the system has N degrees of freedom, the system has N order natural circle frequency ω i and N order mode shapes ϕ i (i � 1, 2, . . ., N). e relationship between natural circle frequency ω i and natural frequency f i (cycle/s or Hz) and natural period of vibration T i (s) is shown in . ., N) corresponding to ω i that satisfies (4) is called the ith order mode shapes.

Establishment of Optimization Model for Updating.
e nonlinear updating of model parameters is essentially the mathematical optimization design of the model. e core idea of optimization design is to get the minimum value of the objective function by establishing an optimization model and using an optimization algorithm under constraint conditions. is process involves design variables (DVs), state variables (SVs), and objective function (OF). DVs are independent variables, which are usually the geometric dimensions, structural physical parameters of the research object. ey are the optimization objects and they change constantly in the whole optimization process. SVs are not independent variables, which change with the change of DVs. SVs are functions of DVs and constraint conditions for optimization design. In general, the existence of SVs is to ensure that the optimization has a clear physical significance. e OF is the minimization variable in optimization design. It is the function of DVs and will change with the change of DVs too. It is the criterion to determine the optimization design scheme. DVs, SVs, and OF are described in the following sections in turn.

Design Variable.
e selection of appropriate DVs is crucial to the success of model updating. It is considered very necessary to conduct the sensitivity analysis of the parameters, which selects the appropriate DVs and improves the updating efficiency. In addition, the natural vibration frequency is only related to the stiffness matrix [K] and mass matrix [M] of the system according to (5). For nonlinear updating of bridge model, the sensitive parameters such as elasticity modulus, inertia moment, and density of each component of the structure can be selected as DVs. e number of DVs and their upper and lower limits are determined comprehensively by the requirements of model updating accuracy and calculation efficiency, and the DVs can be expressed as  x

State
Variable. e dynamic analysis should focus on the structural frequency and mode shape. e absolute value of the difference between f i (the analytical frequency) and f i ′ (the experimental frequency) is written as is selected as the constraint condition, which is required less than or equal to the upper limit value e f . e constraint condition can be expressed as As a quantitative correlation criterion in statistics, the MAC is adopted in mode pairing to correlate the results between the analytical mode shape and the experimental mode shape [26]. e MAC can be calculated as where MAC i is the ith order modal assurance criterion, ϕ i / ϕ i ′ are analytical/experimental mode shape vector, and i is modal order.
e MAC value is between 0 and 1. Generally, if the MAC value is greater than 0.7, it is considered that the two sets of mode shapes can match to a certain extent [26]. e closer the MAC value is to 1, the higher the correlation between the two sets of mode vectors is. Usually, the correlation degree is higher at lower modes, and the higher the order is, the more difficult it is to pair the modes. erefore, when the MAC is used as the SV, the lower limit e mac of the MAC can be set slightly higher at lower modes, while the lower limit e mac of the MAC should be lowered at higher modes. e constraint condition can be expressed as

Objective Function.
e OF can comprehensively consider the minimum sum of the two parts. Part one is the frequency difference between the analytical frequency and the experimental frequency. Part two is the MAC i . e OF can be expressed as where α i and β i are weight.

Optimization Model for Updating.
Once the DVs, SVs, and OF are determined, the optimization model for updating can be established. e essence of the optimization model is to solve the minimum value of a function under certain constraints. e optimization model for updating can be simplified as where x j and x j are the lower and upper limits on the DVs, respectively, and n is the number of the DV; s i and s j are the lower and upper limits on the SVs, respectively, and m is the number of the SV.

Model Optimization Algorithm Based on Penalty Function eory.
According to (12), the optimization model for updating established in this work is a constrained optimization problem. e penalty function method is used to convert the constrained optimization problem into an unconstrained optimization problem [27]. Furthermore, the partial derivatives of the dependent variables (OF and SV) to the DVs are used to determine the search direction, and the optimization results are obtained to minimize the true FE results. e steepest descent method and conjugate direction method are used in the iteration process. Each iteration consists of a series of subiterations, and an iteration is equivalent to multiple analysis cycles to achieve higher optimization accuracy. e constrained problem is written as where F 0 is the unconstrained OF, which is dimensionless; g 0 is the reference function weight, which is adopted to adjust the ratio between the various functions; q k is the response surface parameter of the optimization function, which can improve the optimization accuracy； P x is a penalty function for DV constraints； P s is a penalty function for state variable constraints； k is the number of iteration indicator. e function F 0 is expressed as the sum of two functions, which can be defined as en, F 0 can be replaced as where the functions F g and F q are related to the OF and penalty function, respectively. e basic idea of the unconstrained optimization method is to select a favorable search direction at an approximate point and search in one dimension along this direction to obtain a new approximation point. en the same process is performed for the new point, so that iteration is repeated until the predetermined accuracy requirements are met.
For each iteration (j), a search vector in one direction is defined, and the next iteration is shown as where the line search parameter α (j) is the minimum value of the function F taken in the d (j) direction.
In the initial iteration step, the search method is the steepest descent method, and then the conjugate direction method is used for calculation: where β (j− 1) is the conjugate gradient parameter [27]. When the following equations are satisfied, the iteration converges.
where g (j) and g (j− 1) are the results of the jth iteration and the (j-1)th iteration of the OF, respectively, g (b) is the current optimal OF value, and e g is the OF tolerance.

Establishment of Final Bridge
Model. Based on the NMUT, the structural geometry dimensions, physical parameters, mass matrix, damping matrix, and stiffness matrix of the initial bridge model are updated several times to form the final bridge model. Equation (1) can be written as

Train-Track-Bridge Coupled Dynamic Model
On the basis of the TDIT, the TCDM is established, which is composed of three submodels, including the train, the track, and the bridge, as shown in Figure 2. e train, the track, and the bridge form a whole dynamic system, where the track and bridge are linked through the track-bridge interaction (TBI) while the train and the track are coupled by the WDI. e train submodel is established on the basis of the multibody dynamics theory. e track submodel and the bridge submodel are established by adopting the FEM [5].

Train Submodel.
e train is regarded as a sequence of vehicles, and each vehicle is an MDOF system. It consists of the car body, two bogies, four wheelsets, and springdamping suspension devices. e vehicle is generally considered as a whole mass to analyse the dynamic impact of its vibration on the bridge in the dynamic interaction study of the train-track-bridge system [5]. erefore, to simplify the analysis process, the following assumptions of the vehicle model are used to analyse the train-track-bridge coupled vibration model: (i) e car body, bogie, and wheelset are all rigid bodies; i.e., the elastic deformation of the car body, bogie, and wheelset during vibration is neglected.
(ii) e connection of spring and damper between car body and bogie is called "secondary suspension device", while that between bogie and wheelset is called "primary suspension device". e primary spring and secondary spring can be either linear or nonlinear; while the primary damper and secondary damper can be either viscous and frictional.
(iii) Along the longitudinal direction of the vehicle, the vibration of car body, bogie, and wheelset is not considered.
(iv) Each car body and each bogie only consider two DOFs, including the vertical motion and the pitch motion, and each wheelset considers one DOF, i.e., the vertical motion. e primary and secondary suspensions are modeled as springs and dampers. erefore, for a two-series four-axle vehicle with two bogies, the calculated degree of freedom is ten. e vehicle equation can be obtained by applying D′ Alembert principle to each rigid body. Eventually, the dynamic model of train submodel can be expressed as where M v /C v /K v are the mass/damping/stiffness matrices of the vehicle system; X v / _ X v / € X v are the displacement/velocity/ acceleration vectors of the vehicle system; P v is the load vector of the vehicle system.

Track Submodel.
e ballastless track is generally used for traffic lines on urban rail transit bridges. e rail is simulated as a continuous beam, which is supported by fasteners discretely, and only one DOF of the rail is taken into consideration, i.e., the vertical motion [5]. e dynamic model of track submodel is expressed as where M t /C t /K t are the mass/damping/stiffness matrices of the track system；X t / _ X t / € X t are the displacement/velocity/ acceleration vectors of the track system; P t is the load vector of the track system.

Bridge Submodel.
Since the modeling process of the bridge submodel has been described in Chapter 2 (see (19)), it will not be repeated here.

Dynamic Interaction between Wheel and Rail.
e Hertz nonlinear elastic theory is used in order to describe the WDI [5], which is written as where P(t) is wheel-rail force; G is the contact constant of wheel and rail; R is the radius of wheel; Z w /Z r are the wheel/ rail displacement; Z 0 is rail random irregularity.
In the process of the train crossing the bridge, the train, the track, and the bridge form a whole vibration system. e three subsystems are spatially coupled by the WDI and the TBI. In other words, the load vectors at the right side of the equations above are coupled with each other, and the vibration of the large system satisfies the equilibrium condition of force and the deformation compatibility condition at any time. e Newmark-β numerical integration method can be chosen to solve (23) of the train-track-bridge model [5].

Investigation on Train-Bridge Dynamic
Behaviour of Long-Span UHTCB

Parameters Adopted in Calculation.
is section determines the parameters of the train-track-bridge system. e calculated train is Type-A metro vehicle; the key dynamic parameters of the train are listed in Table 1.
e double-block ballastless track, whose detailed dynamic parameters can be referred to in [28], is applied in this investigation.
A long-span UHTCB is selected as a case study to illustrate its vehicle-bridge dynamic behaviour. e bridge is the main span 480 m semifloating system cable-stayed bridge, and the span layout of the main bridge is 34.  Figure 4. e bridge tower adopts the form of gate type, the height of the low tower is 160 m, while the height of the high tower is 229 m, as shown in Figure 5, and the height difference between the high tower and the low tower is 69 m. e pier is in the form of a column pier, and the foundation is a pile group foundation. Galvanized high strength and low relaxation parallel steel wire whose diameter is 7.0 mm is used for stay cables. 54 pairs are set at the high tower side while 32 pairs at the low tower side.
With the purpose of magnifying the dynamic response difference between the bridge at the low tower location and that at the high tower location as far as possible, the track irregularity is not considered when the train runs the bridge.
When studying the dynamic behaviour of train-bridge system, the representative calculation conditions are selected, as listed in Table 2.

Nonlinear Updating of the Bridge Model.
e bridge model optimization design is carried out based on the NMUT mentioned above. Nine DVs are selected for the updating of the bridge model, respectively, represented by parameters x 1 to x 9 , as listed in Table 3. In this paper, the investigation on the train-bridge dynamic behaviour of the cable-stayed bridge mainly focuses on the main beam structure, so the first four vertical vibration modes and the first two lateral vibration modes of the main beam in the natural vibration characteristics of the bridge are selected as the constraint objects. e ratio (|f i − f i ′ |/f i ′ ) and the MAC i of the six modes are taken as the SVs, which are expressed by the parameters s 1 to s 12 , respectively, as listed in Table 4, and the OF is illustrated in Section 2.2.3. In addition, the upper and lower limits of DVs and SVs can be referred to in [26], as well as the real condition of engineering simultaneously.
After 35 iterations of updating, the modal analysis results of the final updated model are highly in agreement with the experimental modal results. e iterative process of the OF value is shown in Figure 6, and the variations of each SV before and after the updating are illustrated in Figures 7 and  8, respectively. It can be seen that the frequency difference is greatly reduced with the maximum reduction of 8.7%, and the frequency difference after the nonlinear updating is within 3.5%. e modal assurance criterion of each mode is also significantly improved. In addition, the OF value is reduced from 0.0393 to 0.0073 with a decrease of 81.3%. In conclusion, the optimization effect is remarkable, which can be considered as a successful model updating. Due to length limitations of the paper, only the main vibration modes are given here, as illustrated in Figure 9.
As seen from the results, the first-order natural vibration frequency of the bridge is 0.208 Hz, and the period is 4.81 s, which is less than 5 s, indicating that the bridge belongs to the short-period category and the overall stiffness of the bridge is higher. e first-order lateral vibration mode of the main beam appears before the first-order vertical vibration mode of the main beam, indicating that the lateral stiffness of the main beam is smaller than the vertical stiffness. e lateral stiffness is also the key to the bridge design, which conforms to the general law of railway cable-stayed bridge in metro lines. e first-order vertical/lateral vibration mode shape is in the form of half-wave, and the maximum displacement value of the bridge is located near the midspan of the main span, which should be worthy of being paid more attention to in bridge design. High tower first-order lateral vibration mode occurs before the low tower first-order lateral vibration mode, because the high tower is taller than the low tower and its lateral stiffness is lower.

Influence of Stiffness Asymmetry of the Cable-Stayed Bridge on Vibration of the Train System.
e asymmetrical arrangement of UHTCB results in greater stiffness inhomogeneity along the longitudinal direction, which leads to the additional irregularity of stiffness. When the train runs through the bridge, the additional irregularity of stiffness will aggravate the coupled vibration of wheel and rail, which will greatly influence the safety and stability of the train. erefore, in this section, the dynamic response difference of the vehicle between Section-A and Section-B during train  As can be seen from the figures above, with train speed changing, the dynamic responses of the train are different. As can be seen from Figure 10, when a train is running on From the low tower side to the high tower side Focused section e location on the main span, which is 40m-100 m away from the tower (Section-A and Section-B) Observation points L1-L4 and R1-R4     Shock and Vibration the low tower side, its wheel unloading rate is greater than when the train is running on the high tower side. In addition, with train speed increasing, the difference gradually increases. Furthermore, when the train speed is 120 km/h, the wheel unloading rate on the low tower side can reach nearly twice that on the high tower side. As can be seen from Figure 11, the vertical acceleration of the car body when the train is running on the low tower side is much greater than that of the train when it is running on the high tower side. With train speed increasing, the vertical acceleration difference of car body between the low tower side and the high tower side gradually increases; the vertical acceleration difference of the car body indicates the obvious difference in the stability of the train running on the low tower side and the high tower side, which will directly affect the riding comfort. According to Figures 12 and 13, when the train passes through the bridge, there are some differences in the vertical force of the wheelset and derailment coefficient between the low tower side and the high tower side, but the variation law of the difference is not apparent, and the difference is small on the whole.  Shock and Vibration 9

Influence of Stiffness Asymmetry of the Cable-Stayed
Bridge on Vibration of the Bridge System. In this section, the dynamic response difference of the main beam between Section-A and Section-B induced by the TCV and its variation laws are studied under different train speeds and at different positions from the bridge tower. e timehistory data of vertical displacement and vertical vibration acceleration at L1-L4 and R1-R4 are collected, respectively. When the train runs across the bridge at the speed of 120 km/h, the vertical vibration acceleration time-history curves of L1/L4 and R1/R4 are shown in Figures 14 and 15, respectively.
From the comparisons of the two sets of time-history curves, it is found that the trend of the vertical vibration acceleration time-history curves of L1 and L4 is similar under the same train speed and this conclusion can also apply to R1 and R4. In addition, the vertical vibration acceleration of the main beam increases obviously when the observation point is closer to the midpoint of main span, and the maximum value of the vertical vibration acceleration on the low tower side is larger than that on the high tower side under the same operating conditions. Furthermore, as the observation point is gradually away from the bridge tower, the difference of vertical vibration acceleration between the observation points near the high tower side and the observation points near the low tower side also increases obviously.
In order to further analyse the dynamic response of the main beam, the time-domain transformation of the vertical vibration acceleration of the main beam is carried out. e frequency-domain curves are shown in Figures 16 and 17.      From the comparisons of the two sets of frequencydomain curves, we can see that under the same train speed, the trend of the vertical vibration acceleration frequencydomain curves of L1/R1 and L4/R4 is similar. When the observation point is far away from the tower, the amplitude of the main beam is obviously larger. From the frequencydomain distribution of the vertical vibration acceleration of the main beam, when the frequency is low (0.02 Hz -0.6 Hz), the amplitude of the vertical vibration acceleration of the low tower side is generally greater than that of the high tower side. When the frequency is higher (more than 0.6 Hz), the amplitude of vertical vibration acceleration on the high tower side is larger than that on the low tower side. Figures 18 and 19 show the dynamic responses of L1-L4 and R1-R4 when the train runs through the bridge at different speeds.
When the train travels through the bridge at different speeds, the variation of the dynamic response difference between the observation points near the high tower side and the observation points near the low tower side is further shown in Figures 20 and 21. e analysis of the above-mentioned dynamic response data shows that under different running speeds, the dynamic response difference of the main beam at the observation points near the low tower and at the observation points near the high tower is significant. Under some operating conditions, the maximum vertical displacement of the observation points on the low tower side is 52% higher than that on the high tower side, and the maximum vertical vibration acceleration of the low tower side can be 139% higher than that on the high tower side. According to Figure 20, the vertical displacement difference of the observation points at the same distance from the bridge tower is not obviously fluctuated, which is not changed by the change of the speed. With the increasing distance from the observation points to the tower, the vertical displacement difference increases gradually, and the variation trend of the difference is very obvious. From the analysis of Figure 21, it is found that the vertical vibration acceleration difference of the observation points at the same distance from the bridge tower increases with increasing train speed, while the vertical vibration acceleration difference increases with the distance from the observation point to the bridge tower increasing when the speed of train is the same. Furthermore, the vertical vibration acceleration difference increases significantly with the speed increasing when the observation points are getting far from the tower. In general, the vibration of the main beam on the low tower side is larger than that on the high tower side.
From Figures 20 and 21, it can be seen that there is a mapping relationship between the vertical displacement difference and the distance from the observation point to the bridge tower, while there is also a mapping relationship between the vertical vibration acceleration difference and the train speed and the distance from the observation point to the bridge tower simultaneously. erefore, in order to describe the mapping relationship accurately, these data of the vertical displacement difference, the vertical vibration acceleration difference, the train speed, and the distance of observation point to bridge tower are nonlinear fitted. According to the law of data distribution, the vertical displacement difference and the vertical vibration acceleration difference are fitted by plane model and surface model, respectively (see Figure 22). In Figure 22, the differences of vertical displacement and vertical vibration acceleration are, respectively, represented by Δz and Δa, while v represents the train speed and d t indicates the distance from the observation point to the bridge tower.
From Figure 22, we can see that there is a mapping relationship between the dynamic response difference and the speed and the position of the observation point, and the nonlinear fitting coefficients have exceeded 0.95.

(24)
Based on the above equations, the difference of vertical displacement and vibration acceleration of the main beam on both sides of the bridge can be calculated at different train speeds and at different positions of the bridge, which provides the basis for reasonable distribution of the vertical stiffness along the longitudinal direction of the bridge.
For the UHTCB, the asymmetrical characteristics of its structure will lead to the asymmetry and inhomogeneity of stiffness distribution. Due to the difference of stiffness on both sides of the main span of the bridge, the stress on both sides of the structure is not uniform. As a railway cablestayed bridge in metro lines, under the impact of the train for many years, the problem of nonuniform force on both sides of the bridge is enlarged gradually, which causes further structural fatigue problems, resulting in the reduction of the ultimate bearing capacity of the bridge and accelerating the aging of the bridge. erefore, this paper can offer useful reference for similar bridge design, which can be used to control the bridge's nonuniform force, and provide reference for bridge's vibration mitigation and operation.

Conclusions
In this work, primarily, an NMUT based on penalty function theory is proposed, which can be used to optimize the bridge model. en the TCDM is established on the basis of the TDIT. Finally, a long-span UHTCB is applied as the case to investigate its train-track-bridge dynamic behaviour. From above research, the conclusions can be as follows.
(i) is proposed NMUT is efficacious and practical.
After 35 iterations of updating, the modal analysis results of the final updated model are highly consistent with the experimental modal results. (ii) For the long-span UHTCB, with train speed changing, the dynamic responses of the train are different. Generally speaking, when the train travels through the bridge, the difference of the vertical force of the wheelset and the derailment coefficient has no clear variation law, and the value is small. e difference of the wheel unloading rate and the vertical acceleration of the car body is large and intensifies with the increase of the train speed, which will significantly affect the safety and ride comfort of the train.
(iii) For the long-span UHTCB, there is a significant dynamic response difference of the main beam on both sides of the main span of the bridge. e vertical displacement difference of the main beam on both sides of the main span of the bridge increases with the distance from the observation point to the bridge tower increasing. e variation trend of vertical vibration acceleration difference on both sides of the main span of the bridge will be affected by the speed of the train and the distance from observation point to bridge tower simultaneously. In general, the vibration of the main beam on the low tower side is larger than that on the high tower side.

Data Availability
No data were used to support this study.