The Wheel-Rail Contact Problem in Vehicle Dynamic Simulation [thesis]

N. Burgelman
One of the major costs incurred by railway companies is the maintenance of turnouts. This situation occurs because the large dynamic forces between the wheels of a train and the rails of a turnout cause excessive wear, rolling contact fatigue and rapid degradation of other components. A thorough understanding of the dynamic interaction between a train and a turnout could lead to a better design of the vehicle and track structures, deeper insight into the damage mechanisms and subsequently to
more » ... rter maintenance planning. To properly model the interaction between a train and a turnout, three important issues need to be considered: 1. Track flexibility: The rails, the fastening, the sleepers and the ballast are not rigid but can move and absorb or transmit vibrations. This factor is especially relevant when impact between the wheel and the rail occurs, which is always the case in a turnout. 2. Wheel/rail contact: In a turnout the rail geometry is complex at the switch blade and at the frog. Therefore, some of the assumptions commonly made in contact models used in vehicle dynamic software are no longer valid when considering those contacts. 3. Effect of train coupling: In a turnout the forces acting between the vehicles through the couplers have a significant effect on the vehicle dynamics. This is especially true when traction or braking is considered. The first issue has been extensively covered in the literature, and flexible track models are available in modern commercial software for vehicle dynamic simulation. Therefore, the issue is not further investigated in this dissertation. The second issue can be approached in two ways. One can use contact models online for the vehicle simulation, which are advanced while still being fast and robust. Alternatively, one can use a simple contact method to evaluate the contact forces online during the vehicle simulation. Then, a more sophisticated method can be employed offline resolving the contact problems and obtaining a more accurate estimate 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 force calculated with the more advanced method are not coupled back to the vehicle dynamic simulation; therefore generating an error. This error from the offline approach is investigated in this thesis for the determination of the contact point location and for the calculation of the tangential contact forces. A fast analytical method for calculating the longitudinal contact location was validated, but the effect of the contact point location on the wheelset's yaw angle was found to be negligible. vii viii SUMMARY To quantify the effect of different solution methods of the tangential contact forces, we compared four models: FastSim (the reference model, currently used most often in multibody simulations) and three non-elliptical models based on interpenetration: Kik-Piotrowski, Linder and Stripes. These models were applied in a co-simulation between the vehicle dynamic software and Matlab, simulating hunting motion and steady curving. It was concluded that the interpenetration methods predict a better curving behaviour (lower creepages and lower creep forces). The simulations performed using the Kik-Piotrowski method resulted in a slightly shorter hunting wavelength. Although the above-mentioned interpenetration methods are more advanced than FastSim, they still make a number of assumptions that may not be satisfied under some certain contact conditions. Therefore, a new contact method is introduced based on Kalker's full theory with the addition of varying spin creepage. This new method was applied offline in a multibody simulation of a wheel-rail contact in a turnout. The resulting stress and slip distributions were realistic. Moreover, the method converged even in the most difficult contact cases. Thus the method is suitable for online application. The question regarding how the contact models can be validated through measurements is a difficult one, as the local stress and slip in the contact area cannot be measured directly. Many methods measuring the consequences of the stress and slip in the contact area have been proposed in the literature. One possibility is to measure the energy dissipated in the wheel/rail contact. This dissertation presents and compares several methods calculating the dissipated energy, it then compares the results obtained from vehicle dynamic simulations with measured energy values. It is concluded that although the approach is promising, there is room for improvement for both the simulations and the measurements. The third issue is especially relevant for long trains with a braking locomotive at the front or a tractive locomotive at the rear. In these configurations the couplers will be compressed. Because of the curvature of the track in a curve or a turnout, there will be a misalignment between the couple and the adjacent carbodies. This misalignment produces an outwards lateral force, which combines with the lateral centripetal force to create a large force on the wheels. The magnitude of this force needs to be below a certain limit for derailment safety. In this dissertation, an existing method for calculating the coupler angles in curves with a constant radius is extended to calculate the coupler angle in curve transients and turnouts. With the knowledge of the coupler angle, quasi-statics can be used to estimate the lateral wheel rail force and thus the risk of derailment. The results obtained from the quasi-static approach are compared to results obtained from vehicle dynamic simulations. This comparison allows for the definition of a dynamic multiplication factor to be applied to the quasi-static results to obtain a first estimate of the dynamic lateral forces. Such a fast method could be useful when a quick and dirty estimate of the derailment risk is required for a non-daily train configuration or in the early stages of a new vehicle/track design. The analysis of the contact models in this dissertation will help researchers choose between the different models and between implementing the model offline or online. This should enable researchers to accurately model the contact in turnouts so that the deterioration mechanisms underlying the excessive wear and RCF in turnouts can be understood. In turn, this understanding should lead to better maintenance and associated SUMMARY ix cost savings. SAMENVATTING Het onderhoud van spoorwissels is een belangrijke kostenpost voor spoorweguitbaters. Dit komt doordat de hoge krachten tussen wielen van de trein en de spoorstaaf in het wissel hoge slijtage en vermoeiing veroorzaken. Een diepgaand inzicht in de dynamische interactie tussen de trein en het wissel kan een beter ontwerp van het wissel mogelijk maken. Ook zorgt een beter inzicht in de dynamische interactie voor een beter begrip van de slijtagemechanismen, wat op zijn beurt weer leidt tot een slimmere onderhoudsplanning. Om de dynamische interactie tussen trein en wissel accuraat te modelleren moeten de volgende drie factoren in acht genomen worden: 1. de elasticiteit van het spoor: de spoorstaven, de dwarsliggers en het ballastbed kunnen bewegen en zo trillingen absorberen of doorgeven. Dit is zeker van belang als er sprake is van impact tussen het wiel en de spoorstaaf, wat in een wissel altijd het geval is. 2. het wiel/spoorstaafcontact: in een wissel is de spoorstaafgeometrie complex, zowel bij de tong als bij het puntstuk. Hierdoor zijn sommige veronderstellingen die vaak gemaakt worden in de contactmodellen voor voertuigsimulaties niet langer geldig. 3. de invloed van de treinkoppeling: in een wissel hebben de krachten in de koppeling tussen de wagenbakken een significante invloed op de voertuigdynamica. Dit wordt nog belangrijker als ook rekening gehouden wordt met de tractie of het remmen van het voertuig. Het eerste punt, de elasticiteit van het spoor, is uitvoerig beschreven in de vakliteratuur. Verscheidene modellen voor het modelleren van deze elasticiteit zijn beschikbaar in commerciële software voor voertuigsimulaties. Daarom wordt dit punt niet verder behandeld in dit proefschrift. Het tweede punt, het wiel/spoorstaafcontact in het wissel, kan op twee manieren benaderd worden. Ofwel wordt een tamelijk geavanceerd contactmodel gebruikt, online in de voertuigsimulatie. Dit contactmodel moet dan snel en robust genoeg zijn. Een tweede benadering bestaat eruit om een voertuigsimulatie te doen met een relatief eenvoudig contactmodel, waarna de contactproblemen opnieuw worden opgelost met een geavanceerd contact model, om zo een betere benadring van de contactspanningen en de locale slip te verkrijgen. In dit geval worden de contactkrachten berekend met het geavanceerde model niet teruggekoppeld naar de voertuigsimulatie, waardoor het krachtenevenwicht niet gegarandeerd wordt. De fout die hierdoor ontstaat werd in dit proefwerk onderzocht voor wat betreft de locatie van het contact punt en voor de berekening van de tangentiële contactkrachten. Er werd een snelle analytische methode voor xi xii SAMENVATTING de berekening van de locatie van het contact punt gevalideerd door vergelijking met een numerieke methode. Hieruit werd geconcludeerd dat alhoewel de anlytische methode inderdaad een meer accurate locatie van het contact punt vindt, de invloed van de contact locatie op de voertuigsimulatie miniem is. Om de invloed van het tangentiële contactmodel op de vuirtuigsimulatie te quantificeren werden vier contactmodellen toegepast in een co-simulatie. De modellen waren: FastSim (gebruikt als referentiemodel, tegenwoordig het meest gebruikte model voor voetuigsimulaties) en drie niet-elliptische modellen gebaseerd op interpenetratie: Kik-Piotrowski, Linder and Stripes. In de co-simulatie werd de sinusloop en de passage door een boog gesimuleerd, waarbij VIRail gebruikt voor de voertuigsimulatie en Matlab voor de contactmodellen. De conclusie was dat de interpenetratiemodellen een beter booggedrag simuleren (minder kruip en kleinere kruipkrachten). Het Kik-Piotrowski model voorspelde een iets kortere golflengte van de sinusloop. Hoewel de bovenvernoemde modelen geavanceerder zijn dan FastSim gaan zij nog steeds uit van een aantal veronderstellingen die geschonden worden in het geval van contact in een wissel. Daarom werd een nieuw contactmodel geïntroduceerd, dat gelijkvormig contact en een variërende rotatiekruip toelaat. Eerst werd een simpel model (FastSIM) gebruikt tijdens de voertuigsimulatie, vervolgens weren alle contactsituaties herberekend met het nieuwe model. De resulterende spanning-en slipverdeling was realistisch, bovendien convergeerde de methode in alle gevallen. Dit laatste is het bewijs dat het model robuust genoeg is om tijdens de voertuigsimulaties gebruikt te worden. Een empirische validatie van de contactmodellen is moeilijk, omdat de spanning en slipverdeling in het contactoppervlakte niet rechtstreeks gemeten kunnen worden. Daarom zijn er veel methodes ontwikkeld on de gevolgen van de spanning en slip in het contact oppervlak te meten. Een van deze gevolgen is het energieverlies door wrijving in het wiel/spoorstaafcontact. Dit proefschrift presenteert en vergelijkt een aantal methoden om dit energieverlies te berekenen, waarna het berekende energieverlies vergeleken werd met gemeten waarden. De conclusie was dat, alhoewel de methode veelbelovend is, er verbetering nodig is van zowel de inputwaarden van de simulatie, als van de meetmethoden. Het derde punt, de invloed van de koppeling tussen de wagenbakken, is in het bijzonder relevant voor lange treinen met een duwende locomotief achteraan of een remmende locomotief vooraan. In deze configuraties worden de koppelingen tussen de wagons ingedrukt. In een boog of een wissel ontstaat dan een hoek tussen de koppeling en de aangrenzende wagons. Deze hoek zorgt ervoor dat de drukkrachten in de koppeling bijdragen tot de zijdelingse kracht van de wielen op het spoor. Om het ontsporingsrisico te beperken moet deze zijdelingse kracht onder een bepaald niveau blijven. In dit proefschrift werd een bestaande methode om de koppelingshoek te berekenen in bogen met een constante boogstraal, uitgebreid voor transitiebogen en wissels. Als de koppelingshoek bekend is, kan quasi-statica aangewend worden om de zijdelingse krachten en zo het ontsporingsrisico te berekenen. De resultaten verkregen met de quasi-statica werden vergeleken met resultaten van voertuigsimulaties. Deze vergelijking laat ons toe om een vermenigvuldigingscoëfficiënt te bepalen die kan gebruikt worden met de resultaten van de quasi-statica om een eerste schatting te bekomen van de dynamische zijdelingse krachten. Deze werkwijze kan gebruikt worden om een snelle schatting te krijgen van SAMENVATTING xiii het ontsporingsrisico, bijvoorbeeld voor een treinsamenstelling die van de dagelijkse afwijkt of om een eerste schatting te hebben in de eerste fase van een nieuw ontwerp van het train/spoor systeem. De hier gepresenteerde studie van de contactmodellen zal onderzoekers helpen met die keuze tussen de verschillende contactmodellen. Dit zal onderzoekers toelaten om het contact tussen wiel en spoorstaaf in een wissel correct te modelleren en zo de schademechanismen beter te begrijpen. Dit beter begrip kan bijdragen aan het verminderen van de onderhoudskosten.
doi:10.4233/uuid:3cee1830-4179-41f7-b032-ff48aaa9cad3 fatcat:mxvltumsljbjtfgxz6yeoadmxq