On the Computer Formulations of the Wheel/Rail Contact Problem

In this investigation, four nonlinear dynamic formulations that can be used in the analysis of the wheel/rail contact are presented, compared and their performance is evaluated. Two of these formulations employ nonlinear algebraic kinematic constraint equations to describe the contact between the wheel and the rail (constraint approach), while in the other two formulations the contact force is modeled using a compliant force element (elastic approach). The goal of the four formulations is to provide accurate nonlinear modeling of the contact between the wheel and the rail, which is crucial to the success of any computational algorithm used in the dynamic analysis of railroad vehicle systems. In the formulations based on the elastic approach, the wheel has six degrees of freedom with respect to the rail, and the normal contact forces are defined as function of the penetration using Hertz’s contact theory or using assumed stiffness and damping coefficients. The first elastic method is based on a search for the contact locations using discrete nodal points. As previously presented in the literature, this method can lead to impulsive forces due to the abrupt change in the location of the contact point from one time step to the next. This difficulty is avoided in the second elastic approach in which the contact points are determined by solving a set of algebraic equations. In the formulations based on the constraint approach, on the other hand, the case of a non-conformal contact is assumed, and nonlinear kinematic contact constraint equations are used to impose the contact conditions at the position, velocity and acceleration levels. This approach leads to a model, in which the wheel has five degrees of freedom with respect to the rail. In the constraint approach, the wheel penetration and lift are not permitted, and the normal contact forces are calculated using the technique of Lagrange multipliers and the augmented form of the system dynamic equations. Two equivalent constraint formulations that employ two different solution procedures are discussed in this investigation. The first method leads to a larger system of equations by augmenting all the contact constraint equations to the dynamic equations of motion, while in the second method an embedding procedure is used to obtain a reduced system of equations from which the surface parameter accelerations are systematically eliminated. Numerical results are presented in order to examine the performance of var