Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 52, 2, pp. 395-404, Warsaw 2014 NUMERICAL SIMULATIONS OF THE WHEEL-RAIL TRACTION FORCES USING THE ELECTROMECHANICAL MODEL OF AN ELECTRIC LOCOMOTIVE Sławomir Duda Silesian University of Technology, Department of Theoretical and Applied Mechanics, Gliwice, Poland e-mail: slawomir.duda@polsl.pl This paper presents amethodology formodelling themovements of an electrical rail vehicle on any track. Vehicle movement is the result of the drive system of a vehicle transferring electromechanical torque generated by the electric motor onto the wheels of the vehicle, while the appropriate driving force is a result of forces resulting from contact between the wheel surface and the rail.The followingapplicablemodelswereprepared: electromechanical (for electricmotor connections that constitute rail vehicle drive) of the vehiclewith the drive system and of the contact wheel-rail. These models were implemented inMatlab/Simulink. The presented simulation tool for rail vehicles has been used to test a variety of cases, demonstrating capabilities of the proposed method in a wide range of situations. These cases involved analysis of dynamics of the start-up of theEU07 rail vehicle in a straight line. The results obtained for the analysed rail are valuable information that makes it possible to determine both the loading of the locomotive drive system and to estimate the vehicle travelling speed. Keywords: rail vehicle, electromechanical model 1. Introduction Theanalysis of the existentmodeling and simulating applications for complex electromechanical systems in the context of durability and reliability leads to a conclusion that searching for a new algorithm supported by numerical tools has a purpose. Electric railway vehicles constitu- te a separate group. Modern vehicles of this type make up a large group of electromechanical machines which require a specific and interdisciplinary approach to modeling. Dynamic pheno- mena occurring in the above may be described by division of these systems into two groups: electromechanical andmechanical systems. The laws of movement governing such systemsmay be presented in a form of differential algebraic, that is electromechanical equations. The simula- tion research on railway vehicles movement requires preparing both, models of the vehicles and models of the coupled electric motors which constitute the drive, while taking the phenomena which occur at the contact of the rail and wheel into consideration. Currently, complex testing of dynamics of the wheel-rail mechanical system (Kisilowski, 1991) may be conducted using computer numerical simulations (Duda, 2011; Grzesikiewicz and Seńko, 2009). Thenumerical simulations ofmotion of railroad vehicles require implementation of amathematical model of the vehicle, which would describe the vehicle correctly in the software. It is very convenient to apply a methodology based on multibody dynamics in creation of such models (Ambrósio et al., 2012; Kortum and Sharp, 1993). The solution to the problem of the wheel-rail contact requires an accurate description of contact kinematics by calculatingmicroslips, that is the normalized relative speeds at the point of contact. Scientific literature provides broad descriptions of the wheel-rail contact models of various simplification levels, based on the rolling contact theory (Kalker, 1982, 1991). 396 S. Duda Plots of time-histories for dynamic values in these systems depend, to a large extent, on the characteristics and power of the motor driving the machine. The analyses of the impact of the electric motor on themechanical system can be performed bymodelling the drive system as an electro-mechanical system (Mężyk, 2001; Świtoński et al., 2001). This is the way to determine the function of electric motor torque Me(t) that depends on the construction features of the mechanical system and the drive motor. Investigating the dynamic phenomena of rail vehicle drive systems while disregarding the vehicle itself leads to oversimplification. The load of the drive system is very complex. On one hand, there are forces coming from the electric motor constituting the drive of the vehicle, and being also, in some locomotive models, electrically coupled with the drives of remaining axles; on the other hand, there are forces coming from the interaction between the wheel and the rail. These forces are determined by the travel of the vehicle itself, while its motion results from the application of just these forces. It is undoubtedly the coupling between dynamics of the drive system with dynamics of the vehicle, which makes it necessary to take into account the real seating of the rail vehicle drive system in the body constituting the vehicle chassis. Whereas models of railway vehicles, wheel-rail contact or motors – both, direct current commutator and asynchronous that constitute drives of electric locomotives, are long present in the scientific literature, complexmodels which consider mutual coupling of the electrical and mechanical part are an insufficiently described scientific area. This paper presents the methods that can be used to analyze the dynamic phenomena in rail vehicle drive systems including their real seating. The tests are performed by providing a detailed description of electromagnetic phenomena in drive motors and contact phenomena at the interface between the wheel and the rail. The forces determined by using these models represent the load of a rail vehicle drive system. 2. Modelling of the electromechanical system of an electric locomotive To analyze dynamic phenomena of complex electromechanical systems, it is necessary to adopt a physical model of a real object under analysis, representing the most relevant features and phenomena involved, required for dynamical analysis to be performed. Such amodel frequently constitutes a compromise between the accuracy of the object representation and the complexity of descriptionof the involved phenomenaaffecting credibility of the obtained solution, simulation duration, or – in extreme cases – the possibility to obtain a solution. In the case of electrome- chanical systems, the powertrain components – both mechanical and electrical – are mutually coupled dynamic systems. By analyzing dynamic phenomena, particularly in transient states, it is necessary to use a model that makes it possible to provide the electromechanical feedback (Mężyk, 2001). The electromagnetic and mechanical systems of the electric rail vehicle powertrain couple mutually through electromagnetic torque Me and angular rotor velocity ω. For the rail vehicle system presented in Fig. 1, analysis of dynamics of the vehicle and its drive system must take into account the electromagnetic andmechanical elements. Taking into account the unique character of the solution of the electric locomotive drive system (type of its electricmotor, number ofmotors and the configuration of their connections), the paper presents an electromagnetic model of the locomotive chosen for further analysis. For this purpose, the electric locomotive EU07 has been chosen (Duda, 2011). In spite of being aware that this is a bit obsolete model, this one was chosen because of easy access to technical documentation making it possible to model this vehicle in a virtual space to determine basic features of thevehicle, including itsweight andmomentof inertianecessary tobuildamechanical system model. Numerical simulations of the wheel-rail traction forces... 397 Fig. 1. The structure of the drive system 2.1. Developing a model of the non-linear motor connection In their operation mode, the traction motors of the EU07 locomotive work in two configu- rations. During start-up four motors are connected in series; then, to increase the voltage, the motors are switched over to a parallel circuit, twomotors per branch. A serial locomotive motor connection was analyzed. The figure below (Fig. 2) presents a schematic diagram for the motor circuits connected in series. Fig. 2. Assumed substitute diagram of the electric motors connected in series For the above assumed schematic diagram, the following voltage equations and physical relations have been formulated d dt It = 1 4 ∑ i=1 Lti [ Ut− 4 ∑ i=1 eri− It ( 4 ∑ i=1 Rti+RD ) −Rb12(It− Isz12)−Rb34(It− Isz34) ] (2.1) d dt Isz12 = 1 2 ∑ i=1 Lszi ( Rb12(It−Isz12)− Isz12 2 ∑ i=1 Rszi ) (2.2) er1 = cE1f(Isz12)φn1ω1 er2 = cE2f(Isz12)φn2ω2 (2.3) d dt Isz34 = 1 4 ∑ i=3 Lszi ( Rb34(It−Isz34)− Isz34 4 ∑ i=3 Rszi ) (2.4) er3 = cE3f(Isz34)φn3ω3 er4 = cE4f(Isz34)φn4ω4 (2.5) Me1 = It er1 ω1 Me2 = It er2 ω2 Me3 = It er3 ω3 Me4 = It er4 ω4 (2.6) 398 S. Duda d dt ω1 = 1 J1 (Me1−Mob1) d dt ω2 = 1 J2 (Me2−Mob2) d dt ω3 = 1 J3 (Me3−Mob3) d dt ω4 = 1 J4 (Me4−Mob4) (2.7) where It is the tractionmotor armature current, Isz12, Isz34 – currents conducted by excitation windings in motors S1, S2 and S3, S4, kEi – machine constants, φni – streams at normal excitation, f(Iszi) – relative non-linear magnetizing characteristics, Rti – armature resistances, Rb12, Rb34 – shunt resistances for excitation winding, Rszi – serial circuit resistances, RD – total of additional resistances, Lszi,Lt – inductance for excitation windings, ωi – rotor angular velocities, Tei –motor electromagnetic torques, eri – voltages induced in armature circuits. The mathematical model of the motor given by equations (2.1)-(2.7) includes the additio- nal resistance RD and shunt resistance Rb of the excitation winding, whereas the non-linear characteristics of magnetization (Fig. 3) are expressed in equation (2.5). Fig. 3. Relative magnetization characteristic The characteristic feature of a seriesmotor is the fact that the armature current equal to the motor load current is at the same time its magnetizing current. To study dynamics of the drive system,anon-linearmodelof theanalyzedmotorhasbeendevelopedresulting fromitsnon-linear magnetization curve (Duda et al., 2004) that represents the relationbetween themagnetic stream and excitation. This characteristics in relative units, is essentially equal to the characteristic idle operation of themachine at n= const. Being in possession ofmagnetization characteristics and voltage characteristics induced in armature windings by the current E = f(It) at n = const and knownmachine resistances, it is possible to determine themechanical characteristics of the series motor. Manufacturers usually provide the mechanical characteristics of their motors. Themagnetization characteristics are determined from themechanical characteristics at the non-shunted excitation winding, i.e. from the curvewith 100%.By knowing the resistance in the armature circuit, it is possible to determine the voltage drops in the armature circuit at various current values and a constant supply voltage. Taking into account the determined characteristics for the calculated rated stream, it is possible to calculate relative magnetization characteristics presented in Fig. 3. Assuming that the series excitation winding current will change during machine operation in the range of 0.4-1.7Itn, regarding the characteristics presented in Fig. 3, it is possible to approximate a line, calculating the value of the excitation winding inductance coefficient. 2.2. Dynamical model of the rail vehicle At present, a complex study of the mechanical system – the rail-vehicle can be provided by numerical simulations by using a computer as a tool. Numerical simulations of rail vehicles require implementation ofamathematicalmodelof thevehicle inacomputerprogramthatwould describe it with a reasonable precision. Methods based on the multibody dynamics constitute Numerical simulations of the wheel-rail traction forces... 399 a very convenient approach for developing such models. Over the last few years we have been observing an immense evolution in the modelling technique based on the multibody systems, ranging from a simple program to high specialty codesmaking it possible to analyze kinematics and dynamics of any mechanical system, sometimes coupled with models of electric, hydraulic or pneumatic motors. With the development of advanced numerical methods, more and more effective and reliable numerical programs allowing for formulating and solving highly complex problems are created. As a result of analyzing the technical documentation of the rail vehicle, a physicalmodel has been created in form of a multi-rigid-body system mutually coupled with proper kinematic pa- irs and elastic and damping components. Subsequently, its interpretation in the SimMechanics application was obtained. The necessary parameters describing the model, including masses, moments of inertia and dimensions determining the position of kinematic pairs were obtained from a 3Dmodel created by using the Autodesk Inventor software. The remaining parameters, i.e. stiffness and damping suspension components for the first and second unsprung mass re- duction, were found in the documentation submitted by the Zakłady Naprawcze Lokomotyw Elektrycznych (Electric Locomotive Repair Plant) in Gliwice, Poland. Fig. 4. Calculation algorithm for analysis of dynamics of the rail vehicle travelling on any railway track The study of dynamics of electric rail vehicles requires creation of three intercoupledmodels: a vehicle including the drive system, rail and the wheel-rail interface. At the first stage of the rail vehicle modelling process during its run on the railway track, the subsystem models are built separately. Then, themodels are interconnected tomake a complete system. Thismethod was implemented in a proprietary software created in theMatlab environment. The calculation algorithm for analysis of the rail vehicle travelling on any railway track is presented in form of a schematic diagram in Fig. 4. Thepresented algorithmused to develop the computer program for analyzing the rail vehicle travel dynamics on any railway track can be expressed in a few steps: • Assuming the initial conditions for generalized coordinates q(t0) and generalized velocities q̇(t0) as well as determining the initial surface parameter values sr(t0), ur(t0), sw(t0) and uw(t0) related to a specific wheel-rail pair; • Solving the systemof non-linear equations to obtain the surfaceparameters that determine the contact point coordinates related to the specific wheel-rail pair; 400 S. Duda • Calculating normal forces during the contact that are generated as a result of thewheel-rail interaction and depend on the contact surface size; • Calculating creepages, tangent micro slides and spinmoments generated as a result of the wheel-rail interaction; • Determining drive torques for each axle of the wheelset separately; • Adding the forces and torques occurring in the contact related to each wheel as well as adding drive torques to the vector of external forces acting on the system. Using the multibody system formalism to obtain a solution, a new system of generalized positions and velocities for the subsequent time step t+∆t; • Updating the system for another moment by adopting the initial data from the previous step to determine non-generalized surfaces related to each wheel-rail pair; • Continuing the whole process for a new time step until final the time of the analysis is obtained. Detailed relationshipsused for thecalculation carried out inpoints of thepresentedalgorithm can be found in Lankarani and Nikravesh (1990), Polach (1999), Pombo and Ambrósio (2003), Shabana et al. (2008). The examinations have been performed using adequate models of: • the vehicle based on themultibody system formalism in theMatlab/SimMechanics softwa- re, • discrete dynamic models of the electromechanical drive system developed by using the Matlab/SimMechanics software and theMatlab application script, • electric motor connections constituting the rail vehicle drive system developed by using theMatlab/Simulink software, • the wheel-rail contact used to determine the support and guide forces developed by using the proprietary script package from theMatlab software. 3. Verification of the assumptions made while modelling the systems Theverification of the elaboratedmodel of thedrive systemof the electric locomotivewas carried out by comparison of the results obtained from measurements taken during the travel of the real locomotive EU07 over the Katowice-Gliwice route with the passenger rolling stock and the results obtained from numerical simulation. During the travel of the train, values of the current in themain circuit with the resistance control (Fig. 6) of the locomotive were recorded bymeans of a digital camera (Fig. 5). The indications of the instruments were read during the start-up of the rolling stock. In the simulations of the start-up of the locomotive, the mass of 224t of the whole draft of carriagespulledbythe locomotivewas taken intoconsideration.Thismasswas reduceduniformly to particular axles of the locomotive. Figure 6 shows comparison of the curves obtained from the measurements (dashed line)with the results of computer simulation (solid line) for the resistance control carried into effect in accordance with Fig. 7. 4. Numerical simulations for of motion of a rail vehicle As a result of the carried out numerical calculations, the curves of drive system parameters, displacements, vehicle velocity as well as forces and contact torques for two resistance controls were obtained. To compare the curves, the results were listed in two columns (Fig. 7). Numerical simulations of the wheel-rail traction forces... 401 Fig. 5. A view of the operator’s control desk Fig. 6. Course of variations in the main circuit current depending on the start-up time Fig. 7. Resistance control Fig. 8. Course of changes in the electric current of the stator in the auxiliary circuit of motor connections in a serial system as a function of time 402 S. Duda Fig. 9. Course of variations in the electromagnetic torque on the first motor as a function of time Fig. 10. Displacement of the center of mass of the first wheelset as a function of time Fig. 11. Velocity of the center of mass of the first wheelset as a function of time Fig. 12. Longitudinal force at the point of contact of the rail with the left wheel of the first wheelset as a function of time Numerical simulations of the wheel-rail traction forces... 403 Fig. 13. Lateral force at the point of contact of the rail with the left wheel of the first wheelset as a function of time Fig. 14. Traction torque at the point of contact of the rail with the left wheel of the first wheelset as a function of time Because of the similarity of the results obtained for the left and right wheel, the boundary conditions are almost symmetrical – themovement occurs centrally along the railway track axis at a velocity lower than critical, only those for the left wheel are presented. The simulation performed shows that reduction of the operation time of additional resistors results in an increment in the electromechanical torquegeneratedbymotors in theirmain circuits (Fig. 7), which, in turn, yields higher traction forces (Fig. 14). 5. Final conclusions The tool for simulating motion of rail vehicles has been used to test various cases, showing the potential of the presentedmethod in a variety of situations. These cases included analysis of the EU07 rail vehicle dynamics during its travel on a straight railway track at various resistance controls. The created vehicle model, in form of an electromechanical system including dedicated pro- grams for determining the vehicle support and guide forces, provides the possibility of obtaining the following results: • Curves of kinematic parameter changes at any selected point of model components (di- splacements, velocities, accelerations), interaction forces in particular kinematic pairs and at the wheel-rail interface. • The point of contact between the wheel and rail (separately for the rolling base and the flange), sizes of the contact ellipse at the wheel-rail surface as well as the place of their occurrence. 404 S. Duda The results of presented computer simulations allow for coming to the conclusion that the developed algorithms and computer programs are general in their character and can be used to determine structural features of electromechanical systems with similar design forms. Thepresentedmodels enable testing of dynamic phenomena occurring in power transmission systems, especially in transient states such as start-up or change of loading conditions. The models serve for determination of optimal traction parameters for locomotives by, e.g., selection of the gear ratio of the power transmission system. Acknowledgement This paper was realized within the framework of research project No. 6700/B/T07/2011/40 funded by National Science Centre in Poland. References 1. Ambrósio J., Pombo J., PereiraM., AntunesM.,MóscaA., 2012,A computational proce- dure for the dynamic analysis of the catenary-pantograph interaction in high-speed trains, Journal of Theoretical and Applied Mechanics, 50, 3, 681-699 2. Choromański W., 1999, Simulation and optimization in dynamics of railway vehicles (in Po- lish),Prace Naukowe Politechniki Warszawskiej. Transport, 42, OficynaWydawnicza Politechniki Warszawskiej,Warszawa 3. DudaS., 2011,Electromechanicalmodel of rail vehicleEU07(inPolish),Modelowanie Inżynierskie, 41, 47-54 4. DudaS.,KciukS.,MężykA.,TrawińskiT., 2004,Dynamicmodel of the drive systemelectric locomotive (in Polish), ZN KMS, XLIII Sympozjon Modelowanie w Mechanice, 23, 117-123 5. Grzesikiewicz W., Seńko J., 2009, Simulation of traffic safety of wagons during braking long train (in Polish),Logistyka, 4 6. Kalker J.J., 1982, A fast algorithm for the simplified theory of rolling-contact, Vehicle System Dynamics, 11, 1, 1-13 7. Kalker J.J., 1991,Wheel-rail rolling contact theory, Journal of Wear, 144, 243-261 8. Kisilowski J., 1991, The Dynamics of the Mechanical System Rail Vehicle – Track (in Polish), PWN,Warszawa 9. Kortum W., Sharp R.S. (edit.), 1993,Multibody computer codes in vehicle system dynamics, Supplement to Vehicle System Dynamics, 22, 73-84 10. Lankarani H.M., Nikravesh P.E.X., 1990, A contact force model with hysteresis damping for impact analysis of multibody systems,AMSE Journal of Mechanical Design, 112, 369-376 11. Mężyk A., 2001,Analysing and Shaping Dynamic Characteristics of Electromechanical Drives (in Polish), Monografia, Gliwice 12. Polach O., 1999, A fast wheel-rail forces calculation computer code, Vehicle System Dynamics, Supplement, 33, 728-739 13. Pombo J., Ambrósio J., 2003,General spatial curve joint for rail guided vehicles: kinematics and dynamics,Multibody Systems Dynamics, 9, 237-264 14. ShabanaA.A., ZaazaaK.E., SugiyamaH., 2008,Railroad Vehicle. Dynamics. A Computatio- nal Approach, Taylor & Francis Group 15. ŚwitońskiE.,DudaS.,MężykA., 2001,Modelling andoptimizationof the selecteddynamic fe- atures of complex electromechanical systems (in Polish),Materiały Konferencyjne XL Sympozjonu Modelowanie w Mechanice, 255-256 Manuscript received October 1, 2012; accepted for print October 24, 2013