# Effect of the Longitudinal Contact Location on Vehicle Dynamics Simulation

Authors: N. Burgelman; Z. Li; R. Dollevoet
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

This paper investigates the effect of the calculation of the longitudinal location of a wheel rail contact point on the wheelset’s motion in a vehicle dynamic simulation. All current vehicle dynamic software programs assume that the contact between wheel and rail takes place in the vertical plane through the wheelset’s rolling axis. However, when the yaw angle of the wheelset is nonzero, the contact point is situated up to 10 mm from that plane. This difference causes a difference in the yaw moment on the wheelset which is used in the vehicle dynamic simulation. To such an end, an existing analytical method to determine the longitudinal method was validated using a numerical approach. Then vehicle dynamic simulations with both the classic and the new contact location were performed, concluding that using a more accurate contact point location results in a smaller wheelset yaw angle in a vehicle dynamic simulation, although the effect is small.


---

## Body

## 1. Introduction

To assess the wheel rail contact conditions (e.g., at a turnout) a vehicle simulation and a contact method are needed. There are two possible approaches. Either one can use contact models online in the vehicle simulation, which are fast and robust, yet reasonably accurate. Or one can first use a simple method contact method to evaluate the contact forces during the vehicle simulation and then a sophisticated method offline to resolve the contact problems and obtain a more accurate calculation of the local stress and slip distribution in the contact area. In the latter case, the balance of forces is not necessarily fulfilled as the contact forces calculated with the sophisticated method are not coupled back to the vehicle dynamic simulation, therefore producing an error. The difference between the simple and the sophisticated contact method can be in the calculated contact forces or in the calculated contact location. This paper will focus on the effect of the contact location, whereas the effect of the calculation of the tangential contact force is covered by [1].A common assumption in vehicle dynamics simulations is that the contact between the wheel and the rail takes place in the vertical plane through the rolling axis of the wheelset. This is correct when there is no angle of attack; however, when the angle of attack is not zero, the actual contact point will shift in longitudinal direction, opposite to the direction of the yaw rotation of the wheelset; see Figure1. Therefore the yaw moment caused by the lateral force of the rail on the wheel will be smaller than expected; that is, not accounting for this longitudinal shift in the contact point would lead to an overestimation of this yaw moment. As this moment acts to increase the yaw angle, an overestimation of this moment will lead to an overestimation of the angle of attack of the wheelset, which in turn would give an overestimation of the derailment risk, as some derailment criteria are based on the angle of attack (e.g., [2, 3]).Figure 1
A top view of a half wheelset with an angle of attack,γ. The lateral force of the wheel on the rail occurs at the longitudinal contact location defined (xcp), whereas current vehicle dynamic simulations assume that the contact occurs at the principal profile (longitudinal location defined by xpp).In this paper, the effect of the longitudinal location of the contact point on the wheelset’s yaw angle is quantified through vehicle dynamic simulations. Therefore a method is needed to calculate the longitudinal location of the contact point. We propose to use the principle of contact locus, a line of potential contact points on the wheel profile (see Section2.1). In this way, we avoid the need for a true 3D method such as that found in [4, 5]. When the angle of attack is zero the contact locus is the principal wheel profile. However, when the angle of attack is not zero, a method is needed to calculate the contact locus. Here, we use an analytical approximation: Wang’s method [6] (translated as Appendix D of Li’s thesis [7]), based on work of de Pater [8, 9], explained in Section 2.2. This method is fast and so suitable for online evaluation during the vehicle dynamic simulation. Wang’s method is compared to results obtained from the projection of a 3D mesh, explained in Section 2.3. Due to the interpolations and the numerical minimum search required, the 3D meshing method is considered too slow to be used in vehicle dynamics simulation. However, the accuracy of the 3D meshing method only depends on the mesh size; therefore it allows validating Wang’s method.


## 2. Determination of the Contact Location

### 2.1. Contact Locus

To determine the contact point we assume that the rail is straight and smooth in the longitudinal (rolling) direction, whereas in the lateral direction (i.e., along the rail cross section) we still allow rail surface irregularities. This assumption is widely used in vehicle dynamic simulation, although violated in the case of wheel rail contact at the point of a frog or at the gap of an insulated joint. If these are the areas of interest, it is better to use a finite element simulation instead, for example, Oregui et al. [10] for FEM of contact at rail joints or Zhao and Li [11] and Telliskivi and Olofsson [12] for general wheel/rail contact. However, due to the large computational effort, finite element methods cannot replace simpler contact methods for the evaluation of the contact forces in a vehicle dynamics simulation.The assumption of straight and smooth rail allows simplifying the problem by introducing the concept of “contact locus.” A contact locus is a line of points on the wheel surface where contact is possible, given a certain angle of attack and roll angle of the wheelset. A contact locus can be used instead of a wheel profile, so that instead of searching for contact between the wheel and the rail surfaces one can search for contact between the contact locus and the rail profile, so that the contact point search problem is reduced to being in the 2D domain.


### 2.2. Wang’s Method

Wang’s method is an analytical approximation; it calculates a correction of the wheel rotated principal profile, that is, the wheel profile in a vertical plane through the axle of the wheelset. Principal wheel profiles are commonly defined as a list of points (yO(i), zO(i)), where y is the lateral and z the vertical direction, and with the origin in O, which is located at half the tape circle distance (T) away from the wheelset’s center (O′) in lateral direction and at vertical distance R0 from the wheel axis; see Figure 2. This way the local rolling radius in each point on the wheel profile is R(i)=R0-z(i). The rotated wheel principal profiles (xpp(i), ypp(i), zpp(i)) are the principal profile (yO(i), zO(i)) rotated around the centre of the wheelset (O′) with the angle of attack and the roll angle of the wheelset. The wheel principal rotated profile can be calculated by translating the reference frame in which the wheel profile is defined to the center of the wheelset and then rotate it around the yaw angle (= angle of attack) (γ) and the roll angle (θ):(1)xppiyppizppi=cos⁡γsin⁡γ0-sin⁡γcos⁡γcos⁡θsin⁡θ0-sin⁡θcos⁡θ0y0i+Tz0i-R0.Figure 2
Side view of a half wheelset with the reference frame of the wheel profile (O) and the reference frame at the center of the wheelset (O′) with T, half the tape circle distance, and R0, the rolling radius.Wang’s method assumes that at each point on the wheel profile the wheel surface can be approximated as a cone [6]. These cones pass through a point on the wheel profile, with the axes coinciding with the wheelset’s rolling axis and a cone angle equal to the local contact angle. According to Wang [6], the contact locus (xcp(i), ycp(i), zcp(i)) is then found as a correction (Δx(i), Δy(i), Δz(i)) of the rotated principal profiles as follows [6]:(2)Δxi=xcpi-xppi=Risin⁡θsin⁡γ+lxRitan⁡δi,Δyi=ycpi-yppi=Risin⁡θcos⁡γ-Ri1-lx2lx2lytan⁡δi+lzm,Δzi=zcpi-zppi=Ricos⁡θ-Ri1-lx2lx2lytan⁡δi-lym,where δ(i) is the local angle of the wheel surface with respect to the wheelset rolling axis: (3)lx=cos⁡θsin⁡γ,ly=cos⁡θcos⁡γ,lz=sin⁡θ,m=1-lx21+tan2⁡δi.For the derivation of (2)-(3), including the physical meaning of the parameters in (3), we refer to Wang [6] or the English translation in the appendix of Li’s thesis [7].


### 2.3. Spline Interpolation on a 3D Mesh

To validate Wang’s method, a 3D mesh of the wheel surface has been constructed. A matrix of points on the wheel surface is obtained by revolving the wheel profile around the wheelset’s rolling axis, using the parametric equation for a circle through each point on the wheel profile. All the points are then rotated using (1); then a 3D cubic spline surface is fitted to the mesh points so that the resulting surface is the rotated wheel surface. On the intersection line of this spline surface and a plane vertical and parallel to the rails (semitransparent purple plane in Figure 3), the minimum of z is searched, for a given lateral location, using a golden section search. The resulting z and the corresponding x are the coordinates of the contact locus. The surface as well as the contact locus (in red) and the rotated principal profile (blue) are shown in Figure 3 for an angle of attack of 5°. This is done for a better visualization as common angles of attack are 0° to 0.5° for trains and 0° to 2° for trams. The advantage of the meshing method is that it does not make geometric simplifications; thus the accuracy only depends on the mesh size, which can always be refined. The disadvantage is that the minimum search requires a substantial computational effort; therefore the method is not suitable for evaluating the contact point location online in a vehicle dynamic simulation.Figure 3
The 3D mesh with the contact locus as a red line and the principal profile as a blue line, and the plane in which the minimum is to be found as a semitransparent purple surface. For visualization it is calculated here with a very large angle of attack of 5° and a coarse mesh.


## 2.1. Contact Locus

To determine the contact point we assume that the rail is straight and smooth in the longitudinal (rolling) direction, whereas in the lateral direction (i.e., along the rail cross section) we still allow rail surface irregularities. This assumption is widely used in vehicle dynamic simulation, although violated in the case of wheel rail contact at the point of a frog or at the gap of an insulated joint. If these are the areas of interest, it is better to use a finite element simulation instead, for example, Oregui et al. [10] for FEM of contact at rail joints or Zhao and Li [11] and Telliskivi and Olofsson [12] for general wheel/rail contact. However, due to the large computational effort, finite element methods cannot replace simpler contact methods for the evaluation of the contact forces in a vehicle dynamics simulation.The assumption of straight and smooth rail allows simplifying the problem by introducing the concept of “contact locus.” A contact locus is a line of points on the wheel surface where contact is possible, given a certain angle of attack and roll angle of the wheelset. A contact locus can be used instead of a wheel profile, so that instead of searching for contact between the wheel and the rail surfaces one can search for contact between the contact locus and the rail profile, so that the contact point search problem is reduced to being in the 2D domain.


## 2.2. Wang’s Method

Wang’s method is an analytical approximation; it calculates a correction of the wheel rotated principal profile, that is, the wheel profile in a vertical plane through the axle of the wheelset. Principal wheel profiles are commonly defined as a list of points (yO(i), zO(i)), where y is the lateral and z the vertical direction, and with the origin in O, which is located at half the tape circle distance (T) away from the wheelset’s center (O′) in lateral direction and at vertical distance R0 from the wheel axis; see Figure 2. This way the local rolling radius in each point on the wheel profile is R(i)=R0-z(i). The rotated wheel principal profiles (xpp(i), ypp(i), zpp(i)) are the principal profile (yO(i), zO(i)) rotated around the centre of the wheelset (O′) with the angle of attack and the roll angle of the wheelset. The wheel principal rotated profile can be calculated by translating the reference frame in which the wheel profile is defined to the center of the wheelset and then rotate it around the yaw angle (= angle of attack) (γ) and the roll angle (θ):(1)xppiyppizppi=cos⁡γsin⁡γ0-sin⁡γcos⁡γcos⁡θsin⁡θ0-sin⁡θcos⁡θ0y0i+Tz0i-R0.Figure 2
Side view of a half wheelset with the reference frame of the wheel profile (O) and the reference frame at the center of the wheelset (O′) with T, half the tape circle distance, and R0, the rolling radius.Wang’s method assumes that at each point on the wheel profile the wheel surface can be approximated as a cone [6]. These cones pass through a point on the wheel profile, with the axes coinciding with the wheelset’s rolling axis and a cone angle equal to the local contact angle. According to Wang [6], the contact locus (xcp(i), ycp(i), zcp(i)) is then found as a correction (Δx(i), Δy(i), Δz(i)) of the rotated principal profiles as follows [6]:(2)Δxi=xcpi-xppi=Risin⁡θsin⁡γ+lxRitan⁡δi,Δyi=ycpi-yppi=Risin⁡θcos⁡γ-Ri1-lx2lx2lytan⁡δi+lzm,Δzi=zcpi-zppi=Ricos⁡θ-Ri1-lx2lx2lytan⁡δi-lym,where δ(i) is the local angle of the wheel surface with respect to the wheelset rolling axis: (3)lx=cos⁡θsin⁡γ,ly=cos⁡θcos⁡γ,lz=sin⁡θ,m=1-lx21+tan2⁡δi.For the derivation of (2)-(3), including the physical meaning of the parameters in (3), we refer to Wang [6] or the English translation in the appendix of Li’s thesis [7].


## 2.3. Spline Interpolation on a 3D Mesh

To validate Wang’s method, a 3D mesh of the wheel surface has been constructed. A matrix of points on the wheel surface is obtained by revolving the wheel profile around the wheelset’s rolling axis, using the parametric equation for a circle through each point on the wheel profile. All the points are then rotated using (1); then a 3D cubic spline surface is fitted to the mesh points so that the resulting surface is the rotated wheel surface. On the intersection line of this spline surface and a plane vertical and parallel to the rails (semitransparent purple plane in Figure 3), the minimum of z is searched, for a given lateral location, using a golden section search. The resulting z and the corresponding x are the coordinates of the contact locus. The surface as well as the contact locus (in red) and the rotated principal profile (blue) are shown in Figure 3 for an angle of attack of 5°. This is done for a better visualization as common angles of attack are 0° to 0.5° for trains and 0° to 2° for trams. The advantage of the meshing method is that it does not make geometric simplifications; thus the accuracy only depends on the mesh size, which can always be refined. The disadvantage is that the minimum search requires a substantial computational effort; therefore the method is not suitable for evaluating the contact point location online in a vehicle dynamic simulation.Figure 3
The 3D mesh with the contact locus as a red line and the principal profile as a blue line, and the plane in which the minimum is to be found as a semitransparent purple surface. For visualization it is calculated here with a very large angle of attack of 5° and a coarse mesh.


## 3. Results

### 3.1. Comparison of Methods for Contact Locus

Comparing the results obtained from Wang’s method with the results from the 3D mesh for an angle of attack of 1° shows a very good fit; see Figures4 and 5. It is therefore concluded that Wang’s method is sufficiently accurate for use in vehicle dynamic simulations.Figure 4
The projection of the contact locus on theyz-plane, vertical and perpendicular to the rail, for a s1002 profile and an angle of attack of 1°. The whole profile (a) and a zoom-in on the wheel flange (b).
(a)
(b)Figure 5
The projection of the contact locus on theyx-plane (horizontal plane), for a s1002 profile and an angle of attack of 1°.Looking at Figure5, one can see that the largest longitudinal difference between the rotated principal profile and the real contact locus is around 23 mm. However, this largest difference occurs at the highest contact angle (70°), which will normally not be reached in a common vehicle dynamic simulation.


### 3.2. Effect on Vehicle Dynamics

Current vehicle dynamic simulations assume the contact takes place on the principal profile at (xpp, ypp, zpp), whereas it takes place at (xcp, ycp, zcp), calculated in the previous section. The location of the contact point has an effect on the moment on the wheelset transferred from the rail to the wheel (see Figure 1). The moment from the lateral wheel/rail contact force on the wheelset is normally calculated as Mz=Fyxpp, whereas it should be Mz=Fyxcp, where Fy is the lateral force from the wheel on the rail consisting of the lateral components of the normal contact force and the lateral creep force.A cosimulation was set up between the MultiBody Software (MBS) VIRail and Matlab in a similar way as presented in [1]. At each timestep of the simulation the angle of attack, the roll angle, the lateral creep force, the normal contact force, and the contact angle are passed from the MBS to Matlab. Matlab then calculates the lateral force and the difference in moment (ΔMz=ΔxFx) for each wheel. The resulting moment is then fed back to the MBS where it is applied through an actuator.The differences in longitudinal location and the resulting moment for a case of a vehicle entering a curve with radius 100 m (realistic for tram networks) are plotted in Figures6(a) and 6(b). This moment results in a change in the angle of attack in the simulation, seen in Figure 7. It can be seen that the influence of this moment on the angle of attack is very small at around 0.2%. This can be explained by the observation that the moment applied through the actuator is small compared to the other yaw moments that act on the wheelset, more specifically, the yaw moment caused by the longitudinal creep force. The longitudinal creep force causes a yaw moment equaling FxT, where T is half the tape circle distance, and even though Fx is generally smaller than Fy, Δx is much smaller than T, so that the moment originating from the lateral contact force becomes negligible.Figure 6
The longitudinal distance between the contact locus obtained with Wang’s method and the principal profile at the point of contact (a), and the difference in yaw moment due to the assumed location of the contact patch (b). This difference in moment is applied through the actuator in the MBS.
(a)
(b)Figure 7
The angle of attack of the first wheelset of a vehicle entering a curve with radius 100 m from a simulation with the contact on the rotated principal profile and for a simulation with the contact location obtained through Wang’s method (a) and a zoom-in on the moment steady curving is reached (b).
(a)
(b)We conclude that, for the case of a vehicle entering a 100 m curve, resulting in an angle of attack of 1° and a contact angle of 40°, the influence of the longitudinal contact location is negligible. In tramway the angle of attack can be as high as 2°, and in case of a derailment the contact angle can be up to 70°. However, such angles of attack will only occur in very tight curves so that the difference in rolling radius between the left and right wheel will not compensate for the difference in traveled distance on the left and right rail; therefore in such cases there will always be a high longitudinal creep force so that the change in moment due to the longitudinal contact location can safely be neglected.


## 3.1. Comparison of Methods for Contact Locus

Comparing the results obtained from Wang’s method with the results from the 3D mesh for an angle of attack of 1° shows a very good fit; see Figures4 and 5. It is therefore concluded that Wang’s method is sufficiently accurate for use in vehicle dynamic simulations.Figure 4
The projection of the contact locus on theyz-plane, vertical and perpendicular to the rail, for a s1002 profile and an angle of attack of 1°. The whole profile (a) and a zoom-in on the wheel flange (b).
(a)
(b)Figure 5
The projection of the contact locus on theyx-plane (horizontal plane), for a s1002 profile and an angle of attack of 1°.Looking at Figure5, one can see that the largest longitudinal difference between the rotated principal profile and the real contact locus is around 23 mm. However, this largest difference occurs at the highest contact angle (70°), which will normally not be reached in a common vehicle dynamic simulation.


## 3.2. Effect on Vehicle Dynamics

Current vehicle dynamic simulations assume the contact takes place on the principal profile at (xpp, ypp, zpp), whereas it takes place at (xcp, ycp, zcp), calculated in the previous section. The location of the contact point has an effect on the moment on the wheelset transferred from the rail to the wheel (see Figure 1). The moment from the lateral wheel/rail contact force on the wheelset is normally calculated as Mz=Fyxpp, whereas it should be Mz=Fyxcp, where Fy is the lateral force from the wheel on the rail consisting of the lateral components of the normal contact force and the lateral creep force.A cosimulation was set up between the MultiBody Software (MBS) VIRail and Matlab in a similar way as presented in [1]. At each timestep of the simulation the angle of attack, the roll angle, the lateral creep force, the normal contact force, and the contact angle are passed from the MBS to Matlab. Matlab then calculates the lateral force and the difference in moment (ΔMz=ΔxFx) for each wheel. The resulting moment is then fed back to the MBS where it is applied through an actuator.The differences in longitudinal location and the resulting moment for a case of a vehicle entering a curve with radius 100 m (realistic for tram networks) are plotted in Figures6(a) and 6(b). This moment results in a change in the angle of attack in the simulation, seen in Figure 7. It can be seen that the influence of this moment on the angle of attack is very small at around 0.2%. This can be explained by the observation that the moment applied through the actuator is small compared to the other yaw moments that act on the wheelset, more specifically, the yaw moment caused by the longitudinal creep force. The longitudinal creep force causes a yaw moment equaling FxT, where T is half the tape circle distance, and even though Fx is generally smaller than Fy, Δx is much smaller than T, so that the moment originating from the lateral contact force becomes negligible.Figure 6
The longitudinal distance between the contact locus obtained with Wang’s method and the principal profile at the point of contact (a), and the difference in yaw moment due to the assumed location of the contact patch (b). This difference in moment is applied through the actuator in the MBS.
(a)
(b)Figure 7
The angle of attack of the first wheelset of a vehicle entering a curve with radius 100 m from a simulation with the contact on the rotated principal profile and for a simulation with the contact location obtained through Wang’s method (a) and a zoom-in on the moment steady curving is reached (b).
(a)
(b)We conclude that, for the case of a vehicle entering a 100 m curve, resulting in an angle of attack of 1° and a contact angle of 40°, the influence of the longitudinal contact location is negligible. In tramway the angle of attack can be as high as 2°, and in case of a derailment the contact angle can be up to 70°. However, such angles of attack will only occur in very tight curves so that the difference in rolling radius between the left and right wheel will not compensate for the difference in traveled distance on the left and right rail; therefore in such cases there will always be a high longitudinal creep force so that the change in moment due to the longitudinal contact location can safely be neglected.


## 4. Conclusions

The exact longitudinal location of the contact point can be calculated for normal wheel/rail contact fast and accurately using Wang’s method, which was validated using a 3D mesh. The contact location calculated with Wang’s method was used in a vehicle dynamic simulation. The results from this simulation were compared with results where the contact location was calculated in the classical way. It was concluded that more accurate contact location from Wang’s method results in a smaller angle of attack.However, even though the extra computational time for obtaining such accurate longitudinal contact location is small, also the benefits are small, as simulations have shown that the effect of the longitudinal location of the contact point on the vehicle simulation is negligible for contact angles up to 40° and angles of attack up to 1°; for derailment scenarios with contact angles up to 70° further investigation is necessary.



---

*Source: 1901089-2016-01-27.xml*