How to Perform a Quasi-Static Curving Analysis with SIMPACK Wheel/Rail

Transcription

1 How to Perform a Quasi-Static Curving Analysis with SIMPACK Wheel/Rail Version Table of Contents 1 Introduction Methods Solving a nonlinear equation system Time integration Single calculation vs. parameter variation 2 2 Single calculation Setting up the model Track settings Vehicle speed Calculation parameters Result outputs Common results Wear index Roll coefficient Unbalanced lateral acceleration 7 3 Parameter variation Setting up the model Track settings Vehicle velocity Setting up the parameter variation Defining the parameters in the GUI Using parameter variation steering files Ensuring equal travelled distances with different speeds Calculation parameters Result outputs Plotting the results 13 4 Check list for set up 13 5 Example model 15 INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 1

2 1 Introduction Quasi-static curving analysis is a common task in the design process of railway vehicles. The aim is to know how the vehicle behaves in curves, to know the wheel-rail forces and an index of the wheel wear to be expected. This documentation describes how to perform this analysis with SIMPACK. It is strongly recommended that you have attended a SIMPACK Basics training as well as a SIMPACK Wheel/Rail training before using this documentation. The parameter variation requires a SIMPACK Virtual Testing Lab license. 1.1 Methods The basic method to assess the curving behaviour is to find the vehicle s equilibrium state in the curve. For an equilibrium to be found, the track must not have disturbances or irregularities, and the curve radius or the superelevation must not change within the curve. Since the vehicle is not standing but travelling with constant speed, we are not talking about a real equilibrium but a moving equilibrium where the speed is constant but not zero. Thus it is called a quasi-static rather than a static equilibrium. In general, there are two basic methods for finding an equilibrium Solving a nonlinear equation system This method describes the multi-body system with a non-linear algebraic equation system and finds the equilibrium by searching iteratively for a root of this equation system: & z ( z&, z, p, λ, t) 0 (see SIMPACK documentation to static equilibrium module). This method is very fast, but the equation system may have several different equilibrium states, and it depends on the initial configuration which equilibrium state is found; thus it is not recommended for curving analysis Time integration This method uses a time integration, letting the vehicle travel from a straight track into the curve. It is therefore very close to the reality and ensures that the correct equilibrium state is found in the curve, even considering bumpers and other strongly nonlinear kinematics, e.g. in tilting mechanisms, as well as active elements. In this document we concentrate therefore on the time integration method. 1.2 Single calculation vs. parameter variation In SIMPACK there are two possibilities to perform a quasi-static curving analysis. The first method is the manual start of a single calculation. The second method uses the parameter variation and allows several calculations with different parameters to be executed automatically, along with an automatic filtering of the output values in up to three stages. This document describes both methods. INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 2

3 2 Single calculation 2.1 Setting up the model The calculation can be performed with any vehicle model and any wheel-rail contact model. The one-point contact, however, has the best performance and is suitable for main-line profiles and not too narrow curve radii. The only necessary settings are curve radius R, superelevation u, vehicle velocity v. From these parameters results the unbalanced lateral acceleration a q as follows: 2 v u a q = cosϕ g sin ϕ with ϕ = arcsin R e 0 where e 0 is the rail base (in SIMPACK: reference length for superelevation) Track settings Use the standard track type curve entry. It consists of a straight track, a transition curve/superelevation ramp and a curve with constant radius and superelevation. Choose preferredly s-shaped ramp to avoid the setting of the smoothing length. (The main reason for using a standard track rather than a cartographic track is that the standard tracks can be modified in a parameter variation.) Choose a curve radius and a superelevation. Set the length of the straight section so that it is a little longer than the vehicle itself. Usually a length suffices which is five meters longer than the vehicle. Set an appropriate length of the transition curve. For high speeds, high superelevations and/or small curve radii, use longer transition curves. The following table gives some (generously estimated) examples, assuming a velocity yielding a q 1 m/s². If in doubt, tend to longer transition curves. Curve radius Transition length for u=0 m for u = 0.1 m for u = 0.15 m 50 m 10 m 75 m 120 m 100 m 15 m 75 m 120 m 200 m 25 m 75 m 120 m 300 m 50 m 75 m 120 m 500 m 75 m 75 m 120 m 1000 m 100 m 100 m 120 m > 1000 m m m m Set the total track length. Consider that the vehicle will need some time, usually about 10 to 15 seconds, to reach the equilibrium in the curve Vehicle speed Set a vehicle speed that yields the intended unbalanced lateral acceleration. Since the vehicle will lose speed due to frictional energy loss in the curve, it is recommended to change the joint type of the car body from 07 to 09 (v=const.). This resembles that the INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 3

4 vehicle is pulled through the curve by a locomotive. Keep in mind that then the first joint state (s) gets lost and the others (y, z, ϕ, ψ, γ) are shifted one position up. 2.2 Calculation parameters The curving analysis should be executed with the SIMPACK module Static Equilibrium rather than with the normal time integration. The Static Equilibrium module allows a global (inertia) damping in roll and lateral direction (ϕ and y) to be applied. These directions follow the track even in the curve. The damping prevents the vehicle from going into a limit cycle in the curve, which is a small periodic lateral and yaw movement, comparable to an instability in straight track. In the curve it is not generally harmful but corrupts the analysis results. Usually a rather small inertia damping is sufficient. Recommended values are about 100 to 1000 Ns/m lateral and 10 to 100 Nms/rad in roll direction, respectively. 2.3 Result outputs Common results The following table explains how to get the most important results. Since the equilibrium is the intended end state of the calculation, only the end values of the outputs are important as results. Description Q Vertical wheel-rail force Y Lateral wheel-rail force T x Longitudinal creep force Y/Q Derailment coefficient From where Output value 26 of wheel-rail force elements ($F_RW_Friction_of_ ) Output value 25 of wheel-rail force elements ($F_RW_Friction_of_ ) Output value 5 of wheel-rail force elements ($F_RW_Friction_of_ ) Output value 28 of wheel-rail force elements ($F_RW_Friction_of_ ) Remarks Sum of all contact points of the wheel; in track system, acting towards the rail Sum of all contact points of the wheel; in track system, acting towards the rail Attention: only for the current contact point. In contact patch system, acting towards the wheel Sum of all contact points of the wheel W s Wear index in [Nm/m] not directly available May be calculated by an expression, see below ΣY Lateral wheelset force ψ Yaw angle ( angle of attack) Output value 29 of wheel-rail force elements ($F_RW_Friction_of_ ) Joint state 5 of wheelset joint In track system, acting towards the rails y Lateral wheelset position Joint state 2 of wheelset joint in track s R Roll coefficient not directly available May be calculated by an expression, see below INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 4

5 Description a q Unbalanced lateral acceleration From where Control element 167 ( AQ- Acceleration sensor ) Remarks At wheelset marker, see below Wear index The wear index, which relates the frictional energy dissipated in the contact patch to the travelled distance, can be calculated by dividing the frictional power P by the vehicle speed v: W s P 1, Wheel = Wheel, v either in the General Plot, or in an expression within the model: FORCEOV($F_RW_Friction_, 27) / JOINTST($J_Wheelset, 1, 1) Here the frictional power is obtained from the appropriate output value of the wheelrail tangential force element. The vehicle speed comes directly from the wheelset joint. The best way to visualise an expression result in the General Plot is to define a dummy force element of type 50 (force by expression) and plot the output value of this force element: INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 5

6 2.3.3 Roll coefficient The roll coefficient is in kinematic or dynamic terms, respectively: ϕ Carbody, Track s R = or ϕ Track, Ground s R && y = & y Carbody It should be calculated by an expression as well. The easiest form is the kinematic one: Track IF( ABS(AX($M_Isys_Track_Frame_of_,$M_Isys_Track_Camera_of_ ))-1.e-2: 0, 0, AX($M_Body,$M_Isys_Track_Frame_of_ ) / AX($M_Isys_Track_Frame_of_,$M_Isys_Track_Camera_of_ ) ) It uses the fact that the Track Frame marker is canted with the track whilst the Track Camera marker is not. (The IF avoids a division by zero when the track itself is not canted.) Since the Track Frame and Track Camera marker are not automatically created for a car body or a bogie frame, you should either use the markers of the nearest wheelset or define new markers of this kind, see the following pictures. Here $J_Body is the carbody joint. INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 6

7 2.3.4 Unbalanced lateral acceleration If the unbalanced lateral acceleration a q itself is desired as output value, it has to be measured using the control element 167 ( AQ-Acceleration sensor ), which takes the gravity into account. This sensor should be defined in a way that the measuring marker is located on a wheelset, see the picture below. (Although the unbalanced lateral acceleration is originally defined on track level, the Track Frame marker must not be used here because it does not experience an acceleration.) Note that the sensor outputs the acceleration seen by the wheelset, the sign of which is contrary to that of the acceleration experienced by the passengers. INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 7

8 3 Parameter variation 3.1 Setting up the model In a parameter variation the settings for curve radius R, superelevation u and vehicle velocity v can be varied automatically, leading to different configurations in the calculation runs. Whilst in a manual calculation the parameters had to be calculated manually for the desired residual lateral acceleration, the parameter variation allows these parameters to be determined automatically. Generally there are four possible scenarios: Given (fix or varied) parameters Resulting parameter 1 v, R, u Unbalanced lateral acceleration a q 2 a q, R, u Vehicle speed v 3 a q, v, u Track radius R 4 a q, R, v Superelevation u The first scenario uses only given parameters that are explicitly defined in the model. The other scenarios have the implicit parameter a q as given parameter. In the parameter variation module it is possible to handle the scenarios 2 and 3 automatically, i.e. to have vehicle speed or track radius set by SIMPACK, yielding the desired a q. For the scenario 4 this function is not yet available. The necessary settings will be described in the the following Track settings The track settings are the same as for the single calculation, see section If values will be varied in the parameter variation, set arbitrary dummy values for them. Since the length of the transition curve cannot be varied in the parameter variation, set it to the longest curve resulting from the parameter combinations that will be used. Ensure that the track is long enough even for the highest vehicle speed used in the simulation. Otherwise the calculation will fail with an error message Vehicle velocity Set a vehicle velocity that yields the intended unbalanced lateral acceleration. If the velocity is to be varied in the parameter variation, use an arbitrary dummy value (but not zero). 3.2 Setting up the parameter variation There are two different methods to set up the parameter variation. The first method is to define the varied parameters and their variation range in the parameter variation GUI, the second method makes use of parameter variation steering files. The first method is easier to set up, the second one is more flexible. INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 8

9 Although the SIMPACK parameter variation allows nested variation loops, we concentrate on the case of one loop. Feel free to use the middle loop for variation of other parameters. (The outer loop is not useable here.) Defining the parameters in the GUI SIMPACK is able to calculate the vehicle velocity or the track radius automatically according to a given unbalanced lateral acceleration. The following table shows some common scenarios. Since the definition in the GUI allows only equidistant variation steps, it is not reasonable to vary more than one parameter, e.g. R and u, at once. Given, constant Given, varied Automatically calculated lated 1 R, u v 2 u, a q R v = f(a q, R, u, e 0 ) 3 R, a q u v = f(a q, R, u, e 0 ) 4 v, u a q R = f(a q, v, u, e 0 ) All these parameters are of the main type wheel-rail global. They have the following sub-types: Vehicle velocity v: type 01 in [m/s] or 04 in [km/h]. Track radius R: type 70. Superelevation u: type 71. R Unbalanced lateral acceleration a q, leading to v: type 72, v = ( aq g sinϕ) cosϕ. 2 v cosϕ Unbalanced lateral acceleration a q, leading to R: type 73, R =. a g sinϕ q ϕ is here the track roll angle, and g is the gravitational acceleration. The resulting values of type 72, 73 and 74 are not rounded and may therefore be quite unaccustomed. If R, u or v are constant, they are set in the model setup rather than in the parameter variation. If a q shall be constant, however, it has to be defined as a varied parameter, but with equal start and end value. The automatically calculated parameters have to be defined as varied parameters as well. The following pictures show the scenario 2 from the table above, with seven track radii from 300 to 1000 m and a constant a q of 1 m/s² (the superelevation is taken from the model setup): INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 9

10 3.2.2 Using parameter variation steering files Parameter variation steering files allow arbitrary variation steps and parameter combinations. It is possible to define all parameters to be varied at once in the steering file or to combine the steering file with additional parameters defined in the GUI. Parameter variation steering files have the extension.pvs and are stored in the SIM- PACK data base directory parvar_steering_files. To keep the data structure clean it is recommended to use the model or project specific data base rather than the user specific data base when working with steering files. The steering files used here must be of type 10 ( input parameters ). This type allows arbitrary model parameters in.sys file format to be set. The wheel-rail global parameters, which have no according.sys file entries, have the following names: Vehicle speed v: #Vehicle.v_m/s or #Vehicle.v_km/h Track radius R: #Vehicle.track.R Superelevation u: #Vehicle.track.u Automatically calculated v: #Vehicle.aq:v(aq,R,u) INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 10

11 Automatically calculated R: #Vehicle.aq:R(aq,v,u) This is an example for a parameter variation steering file: 1! Example steering file for static equilibrium in curves 10! file format = 10: input parameters 4! number of variations 3! number of input parameters '#Vehicle.v_km/h' '#Vehicle.track.R' '#Vehicle.track.u'! varied parameters ! R=500 m, u=0.1 m ! R=1000 m ! u=0.15 m Once the steering file is finished, it has to be selected in the parameter variation configuration GUI. The number of variations is automatically read from the file and cannot be changed: INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 11

12 3.2.3 Ensuring equal travelled distances with different speeds Usually, when the vehicle speed is varied, the vehicle travels different distances in each parameter variation calculation. Often this is undesirable. In a parameter variation there are two different methods to ensure that the travelled distance is the same, regardless of the vehicle speed. 1. Set the.num6 file parameters tend (simulation end time) and ntout (number of output time steps) in the parameter variation steering file according to the vehicle speed: 1! Number of comment lines! Comment line 10! Steering file format type 4! Number of variations 4! Number of varied parameters #Vehicle.track.R #Vehicle.v_m/s tend ntout Use the force element 234 Stop Integration (available with SIMPACK 8.616) to stop the integrator when the vehicle has reached a given position in the track. See the SIMPACK documentation to force element 234 and the following example, which stops the integration when the first wheelset has reached a position in the track of 500 m. The configuration is shown in the following picture: INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 12

13 3.3 Calculation parameters The curving analysis should be executed with the SIMPACK parameter variation module Static Equilibrium rather than with the normal time integration. This parameter variation uses the settings from the manual static equilibrium module, see section Result outputs Unlike the manual time integration, the parameter variation yields only outputs that have been previously defined in the parameter variation configuration. Otherwise even short calculations could produce huge amounts of output data. The type of outputs is generally the same as it is in a manual calculation see section 2.3 for their definition. They are configured on the notepage Results : (The output definitions may be re-used in different models by means of the Templates and the Global MBS Elements feature of SIMPACK.) 3.5 Plotting the results The results can be plotted in the parameter variation plots Static Equilibrium. The results may there also be exported to a.csv file (for Excel). Another possibility is to plot the results in the General Plot. Ensure to select the appropriate parameter variation case for the plot. 4 Check list for set up Here is a short check list for setting up a quasi-static curving analysis with parameter variation and steering files. Please refer to the above sections for further explanations. Step Description Remarks Section INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 13

14 Step Description Remarks 1 Define a track Standard track Curve entry. Appropriate straight track and transition curve/superelevation ramp. S- shaped ramp. Remove all excitations. 2 Set constant parameters Parameters (R, u, v) that will not be varied in the parameter variation. 3 Define additional measurements 4 Configure parameter variation Track Frame/Camera markers, expressions and dummy force elements, additional sensors etc. Close Model Setup before because parameter variation configuration is written to.sys file. Create parameter variation steering file in model specific data base and select it in configuration GUI. Configure further parameters, if necessary. Section , Configure results Use filters if necessary 2.3, Configure static equilibrium calculation Method time integration. Appropriate simulation time and inertia damping. 7 Start calculation ParVariation Perform Static Equilibrium Calculation + Measurements. 8 Plot results Parameter variation plots or General Plot. Define additional measurements, if necessary Sensors, expressions with dummy force elements, markers etc. Start measurements ParVariation Perform Static Equilibrium Measurements only. If there are new active elements (providing forces or torques etc.), the calculation must be restarted instead. Plot results 2.2 INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 14

15 5 Example model The example model Example_StaticEquilibriumInCurves is a standard railway vehicle with pre-configured static equilibrium calculation using the parameter variation. It uses a steering file. Expressions have been defined for the wear indices and the roll coefficient. The model uses the Stop Integration element to ensure equal travelled distances. INTEC GmbH, Argelsrieder Feld 13, D Wessling, Tel , Fax , intec@simpack.de, page 15