Research Article  Open Access
Ya Gao, Jin Shi, Chenxu Lu, "A TwoStep Composite Time Integration Scheme for VehicleTrack Interaction Analysis considering Contact Separation", Shock and Vibration, vol. 2019, Article ID 1212069, 13 pages, 2019. https://doi.org/10.1155/2019/1212069
A TwoStep Composite Time Integration Scheme for VehicleTrack Interaction Analysis considering Contact Separation
Abstract
The twostep composite time integration scheme is combined with the linear contact relation between the wheel and rail (termed the twostep method) to solve coupled vehicletrack dynamic problems. First, two coupled vehicletrack models with different degreesoffreedom are constructed, in which a vehicletrack system is modeled as two subsystems and the interaction between the two subsystems is implemented via kinematic constraints. Then, the twostep method is applied in the models to simulate several cases, and the accuracy and efficiency of this method are discussed in comparison with additional methods of the Newmark family. Contact separation can also be simulated using the same scheme without causing spurious oscillations for large integration steps. The results indicate that track irregularities can cause wheels to momentarily lose contact with the track, causing large impacts between the wheels and rail. Additionally, with the increase of the running velocity, contact separation will occur more frequently. The adoption of the twostep integration method can ensure the accuracy of solutions and significantly increase the computational efficiency; moreover, the matrixes representing the model information will not change during the calculation process, so the substructure data exported from commercial finite element software can be directly adopted.
1. Introduction
The accuracy and efficiency of the integration method to solve dynamic responses of the vehicletrack coupled system are very important. The track is an elastic structure with damping. Vibrations of the vehicle are transmitted to the track through wheelrail contact, resulting in the excitation of track vibrations, which is in reverse act on the vehicle [1]. Meanwhile, the wheel flat [2] and track irregularities [3] are input to the coupled system via the wheelrail force which significantly amplifies the dynamic interaction and generally displays nonlinear properties. The computational cost is greatly increased.
The vehicletrack coupled system can be represented as two sets of equations of motion, for the vehicle and track subsystems. The first group of methods solves the two sets of equations of motion fully coupled by combining the vehicle and track subsystems as a single integrated system that can be defined as the direct coupled method. Yang et al. [4, 5] constrained the wheel to the bridge to eliminate the degreesoffreedom of the vehicles at the contact point between the wheel and rail. Nguyen et al. [6] proposed an approximate linear contact model, which defines that the wheelrail force can be multiplied by the elastic compression and the contact stiffness at the contact point. The equations of motion can be coupled by integrated wheelrail force into the stiffness matrix [7]. Zhai et al. [1] and Sun et al. [8] used the explicit algorithm [9, 10] to solve the equations of motion with Hertzian nonlinear contact theory for simulating the wheel and rail contact.
However, the explicit integration scheme is conditionally stable. To solve the nonlinear dynamic problems with an unconditionally stable integration scheme, the other kind of method treats the vehicletrack system as two subsystems, and the interaction between the two subsystems is implemented via additional constraints at the contact point between the wheel and rail. The equations of motion of the vehicle and the track are solved separately by using iterative procedures, so the kind of method can be defined as the iterative method. Baeza et al. [11] solved the nonlinear vehicletrack interaction problem with the Newmark method in the frequency domain. Zhang et al. [12] used a precise integration method for solving coupled vehicletrack dynamics to avoid calculation of the inversion of the stiffness matrix which can efficiently calculate the solution. Lei and Mao [13] solved two subsystems with the Newmark method by the iteration progress to meet the constraint condition at each time step. Neves et al. [14] and Montenegro et al. [15] used the Alpha method, an improvement of the Newmark method, for analyzing the nonlinear vehicle structure. Sadeghi et al. [16] used the Newmark Beta integration method together with full Newton–Raphson incremental iterative to solve the equations of motion of the vehicletrack system which can reduce computational cost. Uzzal et al. [17] used the Rayleigh–Ritz method to solve the coupled equations of the vehicletrack system by considering wheel flat. Zhang et al. [18] solved the nonlinear coupled system by the combination of the Newmark method and Newton method with consideration of contact separation and substructure surface roughness.
It can be seen that the selection of an integration method is closely related to the assumption of wheelrail contact relation. In general, the direct coupled method has high computational efficiency, but the coupled matrix is time varying, which is inconvenient for using the substructure data obtained from finite element software. Additionally, the solution of the responses is not accurate enough especially in the case of track irregularities or wheel flat [19]. The iterative method can better solve the equations of motion of the vehicle and the track with nonlinear contact relation between the wheel and rail, but for interaction equations of higher order, several iterations are required in the numerical integration process [20]. Hence, these integration schemes are computationally costly.
In this paper, a twostep integration scheme [21–23] associated with contact separation is applied to analyze the dynamics of a vehicletrack coupled model. Furthermore, the calculation efficiency is evaluated by comparison with the methods of the Newmark family. The effects of vehicle speed and track irregularities are investigated on a multiwheel vehicle and multilayer trackcoupled system with contact separation. It is expected that the algorithm can be more widely applied to different types of models.
2. Equations of the System of VehicleTrack Interaction
A conceptual twodimensional model of a vehicle running on a track is illustrated in Figure 1(a). The model consists of two parts: the upper vehicle model and the lower track model. The vehicle is modeled in three degreesoffreedom: , , and . The carbody and the bogie are connected by a springdamper of stiffness and damping , while the bogie and the wheelset are connected by a springdamper of stiffness and damping . The vehicle is moving along the track at a constant velocity .
(a)
(b)
The motion equation of the vehicle system can be written in matrix form as follows:where , , and are the vehicle mass, damping, and stiffness matrixes, respectively; is the vector describing the effects of wheelrail force on the vehicle; is the external force matrix acting on the vehicle; , , and are the acceleration, velocity, and displacement vectors, respectively; and is the wheelrail force vector.
Using a freebody diagram with vertical displacement of car body , vertical displacement of bogie , and vertical displacement of bogie , the motion equations of the vehicle model shown in Figure 1(a) are derived as follows:
The track is simplified as an Euler beam, and the motion equation can be written as follows:where , , and are the track mass, damping, and stiffness matrixes, respectively; is the external force matrix acting on the vehicle; , , and are the acceleration, velocity, and displacement vectors of the track nodes, respectively; and is the influence matrix of the wheelrail force acting on the track.
In the finite element model of the simplified track, the nonnodal wheelrail force is assigned to the nodes of track elements through the influence matrix . The influence matrix can be written as follows:where is the position transforming matrix consisting of 0 and 1 and is the element influence matrix at the nodes due to a unit point load acting at the location (a, b) of the track element, as shown in Figure 1(b), andwhere is the length of the track element. Then, the displacement of the track at the wheelrail contact point can be expressed as follows:
Since track irregularities amplify the dynamic impact of the wheel to the track, the multiwheel vehicle and multilayer track model, as shown in Figure 2(a), can reflect the stability of the solution more accurately. In this case, the equation of motion of the vehicle can be written as follows:where is the vertical displacement and rotation angle of the carbody and bogies (Figure 2(a)) and is the vertical displacement of the four wheelsets.where is the identity matrix; , , and are the masses of the carbody, the bogie, and the wheelset, respectively; and are the moments of inertia of the carbody and the bogie, respectively; and are the primary and secondary suspension dampings, respectively; and are the primary and secondary suspension stiffnesses, respectively; is the longitudinal distance of the primary suspension; is the longitudinal distance of the secondary suspension; and is the acceleration due to gravity.
(a)
(b)
Finite elements are also used to simplify the multilayer track model, and the track between two neighboring fasteners is divided into one element; each element contains 12 degreesoffreedom, as shown in Figure 2(b) [24]. In Figure 2(b), represents the nth (n = 1, 2, 3, and 4) wheelrail force. The stiffness matrix of the track element can be written as follows:where is the track element stiffness matrix; , , and are the stiffnesses of the rail, track slab, and base plate resisting elastic deformation, respectively; and , , and are the element stiffnesses of the fasteners, CA mortar, and subgrade, respectively.
The mass matrix of the track element can be written as follows:where is the mass matrix of track elements and , , and are the masses of the rail, track slab, and base plate elements, respectively.
The damping matrix of the track element can be written as follows:where is the damping matrix of the track elements; and are the proportional damping coefficients; and , , and are the dampings of the fasteners, CA mortar, and subgrade elements, respectively.
3. Nonlinear Contact Relation with Separation
The nonlinear contact relation between the track and the vehicle is represented by a Hertz nonlinear spring. According to the Hertz theory [25], the vertical force between the wheel and rail can be expressed aswhere is the wheelrail force; is the Hertz spring constant; and is the vertical relative compressive deformation at the wheelrail contact point.
The vertical elastic compression can be obtained by the following equation:where represents the verticaltrack irregularities. Track irregularities have significant dynamic impact on both vehicle and track and should be considered carefully.
When separation of the vehicle and track is considered, it is crucial to check whether separation occurs between the vehicle and track. Since only compressive stress possibly exists between the wheels and rail, the contact force is always nonnegative. Hence, when the relative compression deformation between a wheel and the rail is not less than zero, the wheelrail force is applied on the track, and when it is less than zero, the wheel is separated from the track. In this respect, equation (12) is changed as follows:
Then, the numerical integration method is used to solve equations (1) and (3), and the responses of the system “, ” at the n time steps are obtained. The vertical elastic compression at the n time steps is expressed as follows:where is the number of wheelsets. Once is derived, the wheelrail force is calculated with equation (14) and then used as the initial condition at the (n + 1) time steps:
At each time step, the calculation process is repeated several times until the difference (s is the times of iteration) between adjacent results is less than 10^{−6}.
4. Linear Contact Relation with Separation
The linear approximate contact model assumes that the wheel and rail are connected by a linear spring, representing elastic contact, and the vertical wheelrail force is equal to the product of vertical contact stiffness and vertical elastic compression at the wheelrail contact point [6, 25]. In this section, the relationship between the wheel and the rail is described by the linear contact relation with consideration of separation.
According to the linearized contact relation at the equilibrium point, the wheelrail force is obtained aswhere is the linear contact stiffness.
The contact force is always nonnegative. The wheelrail relative compression deformation is used to judge whether the wheel is separated from the track.
4.1. Solving the Kinematic Equations with the Direct Coupled Method
Combining equations (1), (3), and (17), the equations of motion can be obtained:where and are the additional stiffness matrix and additional wheelrail force matrix caused by the vertical wheelrail contact force, respectively.
Taking the moving vehicle model shown in Figure 1(a) as an example, the additional matrixes and can be expressed as follows:
Because of the linear equilibrium condition of equation (18), the equations of motion of the vehicle and track subsystems are coupled [14] into global equation (19). It can be seen that the wheelrail force is expressed by the additional matrixes and ; both matrixes are updated at each time step. At each time step, equation (18) is still valid. The value of the linear contact stiffness in equations (20) and (21) at the (n + 1) time steps is determined by the vertical elastic compression calculated at the n time steps.
4.2. Solving the Kinematic Equations with the TwoStep Integration Method
The twostep integration scheme [21, 22] is stabilized by the assumptions of the trapezoidal rule and backward difference approximation. This integration scheme is available in some nonlinear problems of structural and solid mechanics against the trapezoidal rule and the Wilson method [21]. In view of the flexibility of the track system, the inertia forces caused by the wheel mass play an important role in the interaction between the vehicle and track. To make the twostep integration scheme remain stable and to maintain the secondorder accuracy of the vehicletrack interaction analysis, the wheel is no longer constrained to the track [23], and the linear equilibrium equation is used as the constraint equation.
Assume that the response at the n time steps is known; the goal is to solve the response at the (n + 1) time steps. Different from the majority of time integration schemes, the twostep integration scheme calculates all unknown variables not only at the (n + 1) time steps but also at the (n + 1/2) time steps at the same time. Then, according to the known responses at the n and (n + 1/2) time steps, the response at the (n + 1) time steps can be solved by using a backward difference.
The kinematic equations consist of two parts: one part is the governing equations of the motion of the vehicle equation (1) and the track equation (3), and the other part is the constraint equation (18). Hence, the coupled equations of motion can be obtained as follows:where is the contact stiffness between the wheel and rail.
According to the twostep integration scheme, in the first substep, the response at the (n + 1/2) time steps is calculated using the trapezoidal rule:where is the time interval.
The coupled equations (equations (22a)–(22c)) at the (n + 1/2) time steps can then be rewritten as follows:where
In the second substep, the (n+1) time steps are calculated with the 3point backward difference method:
Substituting equations (26a) and (26b) into equations (22a)–(22c) yields the coupled equations at the (n + 1) time steps:where
The present algorithm can be implemented through the steps, as shown in Figure 3.
The timevarying coupled relation caused by the vehicletrack interaction is embodied in equations (24c) and (27c). The matrix is not changed in the stepbystep integration method, and iteration is not required in the solution process. Hence, the coupled equations can be solved without modifying the model formulation so that the substructure data exported from commercial finite element software can be directly used.
5. Applications regarding the TwoStep Integration Scheme
5.1. Responses of the VehicleTrack System
The vehicletrack model shown in Figure 1(a) is considered in this section. The track is simplified as an Euler beam continuously supported by springs and dampers, and the vehicle is represented by a multibody system comprising rigid bodies. The running speed of the vehicle is constant at 117 km/h, and only the vibration is considered in the vertical direction. Suppose that the vehicletrack model is laterally symmetric along the line; therefore, a half structure is studied. The parameters of the vehicletrack model are presented in Table 1.

The maximum error of the results [26] in this section is shown, which is defined as follows:where is the response and is the corresponding static response.
If is the displacement, is defined as follows:
If is the acceleration, is defined as follows:
If is the wheelrail force, is defined as follows:
5.1.1. Verification of the TwoStep Integration Method
The effectiveness of the twostep integration method (referred to as “Twostep” in the following text) in solving vehicletrack coupled problems has been verified by comparison with the results simulated with the model in Figure 1(a). The responses of the vehicletrack coupled model with linear contact relation (referred to as “NML” in the following text) and nonlinear contact relation (referred to as “NMN” in the following text) between the wheel and rail are solved by the Newmark Beta method. Moreover, the effect of track irregularities on the model investigated by Newton and Clark (referred to as “NCM” in the following text) has been calibrated by field data [27]. The result extracted from the NCM model is also used as a comparison [20].
The track irregularity selected in this case is representative of wheel flat. In an attempt to overcome the difficulties in locating the positions of the flats and their impacts relative to the track system, an equivalent track irregularity is used instead. The shape of the irregularity is determined based on a large number of test results, and the wheel diameter will be adjusted to fit various vehicle types. The equivalent irregularity can be defined as follows:where is the position along the track; is the maximum depth of the equivalent irregularity, 2.15 mm; is the length of the equivalent irregularity, 150 mm; and is the wheel radius.
Figure 4 presents the analysis results obtained for the model with the track irregularity shown in Figure 5(a). Comparisons are made among the effects of the wheel flat on the wheelrail force computed with the twostep method, NMN, NML, and NCM; all of the results, except for those for NML are in quite good agreement, as shown in Figure 4(a). The highfrequency modal responses of the NCM have been filtered out in Figure 4(a); therefore, compared to the wheelrail forces from the twostep method and NMN, oscillations do not appear in the NCM result after 0.929 s. To show the accuracy of the twostep method more clearly, the accelerations of the midpoint of the track and wheel calculated by the twostep method, NMN, and NML are compared in Figures 4(b) and 4(c). As seen, the results of the twostep method and NMN are in quite good agreement, whereas those for NML become unrealistically large.
(a)
(b)
(c)
(a)
(b)
However, the Newmark Beta method with linear contact relation between the wheel and rail becomes unstable when large deformations and longtime contact separations are considered. The twostep method can remain stable as long as an appropriate time step is chosen. More details about the twostep method regarding computational cost and application will be introduced in the following sections.
5.1.2. VehicleTrack System without Irregularity
In this subsection, the responses of the vehicletrack model without irregularity are solved with the twostep method, and the results are compared with those obtained by the methods of the Newmark family. Two other methods are used, that is, the Newmark Beta method (NML) and the Alpha method (referred to as “APL” in the following text) with linear contact relation between the wheel and rail. The Alpha method, a member of the Newmark family, has better stability and accuracy with introduction of alpha damping [28]. The value of the alpha damping is −0.1 for the case.
As shown in Figures 6(a)–6(c), the responses of the vehicletrack system calculated with the NML, APL, and twostep method are in great agreement. However, the calculation times required by the three methods differ, as shown in Figure 6(d). The calculation time required by the twostep method is much less than that for NML and significantly less than that for APL; it is 5.6% for NML and 61.8% for APL. The integration steps of the Newmark Beta method and Alpha method are 10^{−5} s and 10^{−4} s, respectively, whereas the integration step for the twostep method is 10^{−3} s. The smaller integration step required in the computational process greatly increases the computational cost.
(a)
(b)
(c)
(d)
Figure 7 shows the acceleration of the midpoint of the track for an integration step of 10^{−3} s. Abnormal highfrequency oscillations appear in the response of the track calculated with the Newmark family methods, probably because the solution at the time step of 10^{−3} s does not satisfy the dynamic equilibrium accurately. Calculation errors from each time step gradually accumulate to deteriorate the accuracy of the overall solution.
5.2. Multiwheel Vehicle and Multilayer Track System with Random Track Irregularities
In this section, a multiwheel vehicle and multilayer track model is considered subject to random track irregularities, as shown in Figure 2(a). The random rail irregularities selected are the measured dynamic track irregularity of track inspection car on the Beijing–Jiulong Railway with a design running at a speed of 200 km/h (Sample 1) and the Beijing–Shanghai highspeed railway with a design running speed of 350 km/h (Sample 2), as shown in Figure 5(b). When there is track irregularity, impact loads occur due to contact separation. The responses of the system are solved with the twostep method to show the algorithm stability with highfrequency vibrations induced by impact load. The parameters of the model are presented in Table 2. The length of the track is 100 m, and the integration step is 10^{−5} s.

5.2.1. Responses of the Multilayer Track with Contact Separation
The random rail irregularity of Sample 1 is used in this subsection, and the responses of the track with perfectly smooth rail (without irregularity) are used for comparison. The running speed of the vehicle is 300 km/h.
As shown in Figure 8, the calculated responses of the multilayer track with Beijing–Shanghai measured track irregularity and without irregularity have significant differences during the time history. When oscillations appear in the displacements of the rail, track slab, and base plate in the case with irregularity, the corresponding results are still smooth in the case without irregularity, as shown in Figures 8(a)–8(c). The contact force with irregularity continuously oscillates in the range of 0–250 kN and becomes zero from 0.49 s to 0.50 s, as indicated in Figure 8(d). It shows that, in the case with irregularity, the large amplitude of wheelrail force variation on the substructure leads to stronger dynamic interaction compared with the case without irregularity.
(a)
(b)
(c)
(d)
5.2.2. Influence of Random Track Irregularities on Contact Separation
To obtain more information about contact separation of a multiwheel and multilayer track system with irregularity, both the track irregularities of Sample 1 and Sample 2 are considered. The numbers of separations induced by random track irregularities are shown for all the wheelsets in Figure 9. Wheelrail contact separation is defined when the wheelrail force is zero, and it counts only once from wheelrail separation to recontact. As seen, the front wheelsets have more contact separations than the rear wheelsets of the same bogie during the simulation history. More importantly, the influence of velocity is very important to the number of contact separations. Greater numbers of separations, with increased impacts on the wheel and rail, will be found with increased running velocity.
(a)
(b)
(c)
(d)
6. Conclusions
In this paper, the twostep integration method with the linear contact relation is used to solve the response of a vehicletrack coupled system, and the results and the application conditions are discussed. The following conclusions are drawn.(i)Compared with responses from the Newmark family methods, the responses calculated with the twostep method, NML, and APL are in good accordance, but the calculation time required for the twostep method is much shorter. The adoption of a larger integration step and the avoidance of iterations significantly reduce the calculation time and increase the efficiency of the twostep method. Moreover, the situation of contact separation can be well simulated.(ii)The simulation results indicate that when contact separation occurs owing to random track irregularities, the impact between the wheel and rail will boost the interaction. With an increased vehicle speed, contact separation emerges more frequently.(iii)The high efficiency of the twostep method is also brought about by the integration of all the model information into matrix , which is the usual symmetric matrix. Hence, timevarying coefficients are absent on the left side of the substructure equations, and the substructure data exported from commercial finite element software can be directly adopted.
This method has potential for wide applications in dynamic problems of coupled vehicletrack systems.
Data Availability
The numerical data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest with respect to the research, authorship, and publication of this article.
Acknowledgments
This work was supported by the Fundamental Research Funds for the Central Universities (Grant no. 2018YJS116) and the National Natural Science Foundation of China (Grant no. 51578054).
References
 W. M. Zhai, K. Y. Wang, and C. B. Cai, “Fundamentals of vehicletrack coupled dynamics,” Vehicle System Dynamics, vol. 47, no. 11, pp. 1349–1376, 2009. View at: Publisher Site  Google Scholar
 R. U. A. Uzzal, S. Rakheja, and W. Ahmed, “Dynamic analysis of pitch plane railway vehicletrack interactions due to single and multiple wheel flats,” in Proceedings of the IEEE/ASME/ASCE 2008 Joint Rail Conference, pp. 89–97, American Society of Mechanical Engineers, Wilmington, DE, USA, April 2008. View at: Google Scholar
 J. Sadeghi, M. Khajehdezfuly, and D. Poorveis, “Investigation of rail irregularity effects on wheel/rail dynamic force in slab track: comparison of two and three dimensional models,” Journal of Sound and Vibration, vol. 374, pp. 228–244, 2016. View at: Publisher Site  Google Scholar
 Y.B. Yang and J.D. Yau, “Vehiclebridge interaction element for dynamic analysis,” Journal of Structural Engineering, vol. 123, no. 11, pp. 1512–1518, 1997. View at: Publisher Site  Google Scholar
 Y. B. Yang, M. C. Cheng, and K. C. Chang, “Frequency variation in vehiclebridges interaction systems,” International Journal of Structural Stability & Dynamics, vol. 13, no. 2, Article ID 1350019, 2013. View at: Publisher Site  Google Scholar
 D.V. Nguyen, K.D. Kim, and P. Warnitchai, “Simulation procedure for vehiclesubstructure dynamic interactions and wheel movements using linearized wheelrail interfaces,” Finite Elements in Analysis and Design, vol. 45, no. 5, pp. 341–356, 2009. View at: Publisher Site  Google Scholar
 Z. Li and F. T. K. Au, “Damage detection of bridges using response of vehicle considering road surface roughness,” International Journal of Structural Stability and Dynamics, vol. 15, no. 3, Article ID 1450057, 2015. View at: Publisher Site  Google Scholar
 Y. Sun, Y. Guo, Z. Chen, and W. Zhai, “Effect of differential ballast settlement on dynamic response of vehicle—track coupled systems,” International Journal of Structural Stability & Dynamics, vol. 18, no. 7, Article ID 1850091, 2018. View at: Publisher Site  Google Scholar
 W. M. Zhai, “The predictorcorrector scheme based on the Newmark iteration algorithm for nonlinear structure dynamic analysis,” Chinese Journal of Computational Mechanics, vol. 7, no. 2, pp. 51–58, 1990. View at: Google Scholar
 W. M. Zhai and Z. H. Huang, “Nonlinear numerical methods for train dynamic analysis,” Journal of Southwest Jiao Tong University, vol. 26, no. 1, pp. 82–90, 1991. View at: Google Scholar
 L. Baeza, A. Roda, and J. C. O. Nielsen, “Railway vehicle/track interaction analysis using a modal substructuring approach,” Journal of Sound and Vibration, vol. 293, no. 12, pp. 112–124, 2006. View at: Publisher Site  Google Scholar
 J. Zhang, Q. Gao, S. J. Tan, and W. X. Zhong, “A precise integration method for solving coupled vehicletrack dynamics with nonlinear wheelrail contact,” Journal of Sound and Vibration, vol. 331, no. 21, pp. 4763–4773, 2012. View at: Publisher Site  Google Scholar
 X. Y. Lei and L. Mao, “Analyses of dynamic response of vehicle and track coupling system with random irregularity of rail vertical profile,” China Railway Science, vol. 22, no. 6, pp. 38–43, 2001. View at: Google Scholar
 S. G. M. Neves, P. A. Montenegro, A. F. M. Azevedo, and R. Calçada, “A direct method for analyzing the nonlinear vehiclestructure interaction,” Engineering Structures, vol. 69, pp. 83–89, 2014. View at: Publisher Site  Google Scholar
 P. A. Montenegro, S. G. M. Neves, R. Calçada, M. Tanabe, and M. Sogabe, “Wheelrail contact formulation for analyzing the lateral trainstructure dynamic interaction,” Computers & Structures, vol. 152, pp. 200–214, 2015. View at: Publisher Site  Google Scholar
 J. Sadeghi, A. Khajehdezfuly, M. Esmaeili, and D. Poorveis, “An efficient algorithm for nonlinear analysis of vehicle/track interaction problems,” International Journal of Structural Stability and Dynamics, vol. 16, no. 8, Article ID 1550040, 2016. View at: Publisher Site  Google Scholar
 R. Uzzal, W. Ahmed, and R. Bhat, “Impact analysis due to multiple wheel flats in threedimensional railway vehicletrack system model and development of a smart wheelset,” Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, vol. 230, no. 2, pp. 450–471, 2016. View at: Publisher Site  Google Scholar
 Y. Zhang, H. Zhao, and S. T. Lie, “A nonlinear multispring tire model for dynamic analysis of vehiclebridge interaction system considering separation and road roughness,” Journal of Sound and Vibration, vol. 436, pp. 112–137, 2018. View at: Publisher Site  Google Scholar
 K. Liu, N. Zhang, H. Xia, and G. De Roeck, “A comparison of different solution algorithms for the numerical analysis of vehiclebridge interaction,” International Journal of Structural Stability and Dynamics, vol. 14, no. 2, Article ID 1350065, 2014. View at: Publisher Site  Google Scholar
 J. Sadeghi, A. Khajehdezfuly, M. Esmaeili, and D. Poorveis, “Dynamic interaction of vehicle and discontinuous slab track considering nonlinear Hertz contact model,” Journal of Transportation Engineering, vol. 142, no. 4, Article ID 0401601, 2016. View at: Publisher Site  Google Scholar
 K.J. Bathe and M. M. I. Baig, “On a composite implicit time integration procedure for nonlinear dynamics,” Computers & Structures, vol. 83, no. 3132, pp. 2513–2524, 2005. View at: Publisher Site  Google Scholar
 K.J. Bathe and G. Noh, “Insight into an implicit time integration scheme for structural dynamics,” Computers & Structures, vol. 9899, pp. 1–6, 2012. View at: Publisher Site  Google Scholar
 M. Fedorova and M. V. Sivaselvan, “An algorithm for dynamic vehicletrackstructure interaction analysis for highspeed trains,” Engineering Structures, vol. 148, pp. 857–877, 2017. View at: Publisher Site  Google Scholar
 X. Lei and L. Mao, “Dynamic response analyses of vehicle and track coupled system on track transition of conventional high speed railway,” Journal of Sound and Vibration, vol. 271, no. 3–5, pp. 1133–1146, 2004. View at: Publisher Site  Google Scholar
 J. J. Yang, N. Zhang, and H. Xia, “Comparative analysis on applicability of four wheelrail vertical contact models for coupling vibration of vehicletrack system,” China Railway Science, vol. 37, no. 6, pp. 11–20, 2016. View at: Google Scholar
 D. Y. Zhu, Y. H. Zhang, and H. Ouyang, “A linear complementarity method for dynamic analysis of bridges under moving vehicles considering separation and surface roughness,” Computers & Structures, vol. 154, pp. 135–144, 2015. View at: Publisher Site  Google Scholar
 S. G. Newton and R. A. Clark, “An investigation into the dynamic effects on the track of wheelflats on railway vehicles,” Journal of Mechanical Engineering Science, vol. 21, no. 4, pp. 287–297, 1979. View at: Publisher Site  Google Scholar
 S. G. M. Neves, A. F. M. Azevedo, and R. Calçada, “A direct method for analyzing the vertical vehiclestructure interaction,” Engineering Structures, vol. 34, pp. 414–420, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Ya Gao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.