**Modeling, Simulation, and Results of Their Use in Railway Vehicle Dynamics Studies**

Krzysztof Zboinski

[6] Tsunashima, H., Kojima T., Matsummoto, A. and Mizuma T. Condition monitoring of railway tracks using in-service vehicles. Japanese Railway Engineering. 2008;48(2):

[7] Tsunashima, H., Naganuma, Y., Matsumoto, A., Mizuma T. and Mori, H. Condition monitoring of railway track using in-service vehicle. In: Xavier Perpinya, editor. Reli‐

[8] Tsunashima, H., Mori, H., Mizuma, T. and Saito, Y. Condition monitoring of railway tracks using portable on-board sensing system. In: The 6th International Symposium on Speed-up, Safety and Service Technology for Railway and Maglev Systems

[9] Mori, H., Ohno, H., Tsunashima, H. and Saito, Y. Development of compact size on‐ board device for condition monitoring of railway tracks. Journal of Mechanical Sys‐ tems for Transportation and Logistics. 2013;6(2):142–149. DOI: 10.1299/jmtl.6.142. [10] Ogino, M., Tsunashima, H., Yanagisawa, K., Mori, H, Asano, A. and Wakai, S. Devel‐ opment of track condition monitoring system based on compact onboard sensing de‐ vices (in Japanese). IEEE Transactions on Industry Applications. 2015;135(4):395–402.

[11] Waston, P. F., Roberts, C., Goodman, C. J. and Ling, C. S. Condition monitoring of railway track using in-service trains. In: Railway Condition Monitoring 2006; 29–30

[12] Waston, P. F., Ling, C. S., Roberts, C., Goodman, C. J., Li, P. and Goodall, R. M. Moni‐ toring vertical track irregularity from in-service railway vehicles. Proceedings of the Institution of Mechanical Engineers, Part F, Journal of Rail and Rapid Transit.

[13] Waston, P. F., Ling, C. S., Goodman, C. J., Roberts, C., Li, P. and Goodall, R. M. Moni‐ toring lateral track irregularity from in-service railway vehicles. Proceedings of the Institution of Mechanical Engineers, Part F, Journal of Rail and Rapid Transit.

ability and Safety in Railway. InTech; 2012. 333–356. DOI: 10.5772/35205.

6–10

(STECH2012); 2012.

DOI: 10.1541/ieejias.135.395.

2007;221(F1):75–88.

2007;221(F1):89–100.

Nov. 2006; Birmingham: IET; 2006. 26–31.

164 Railway Research - Selected Topics on Development, Safety and Technology

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/62105

## **Abstract**

This chapter focuses on problems related to building mathematical and numerical models of railway vehicle dynamics and then using these models in the process of vehicle dy‐ namics simulation. Finally, the results of such simulations devoted to selected dynamical problems are presented, highlighting the importance of powerful tools such as both the modeling and the simulation. The dynamical problems selected for the presentation con‐ cern railway vehicle stability and importance of kinematics accuracy for the description of the dynamics. These selected problems focus on the vehicle dynamics in a curved track, both in the circular and transition sections. Type of the chapter should be defined as the review paper, however, based on the authors' own results in the main.

**Keywords:** Railway vehicle, vehicle dynamics, curved track, transition curve, numerical simulation

## **1. Introduction**

The present review chapter is based on the authors' results gathered and published in subsequent parts for many years of his work in the field of railway vehicle dynamics with focus on a curved track motion. The idea of the chapter is to combine all the results on the one hand and to select them suitably on the other hand. Both these demands are fulfilled in order to give a picture of comprehensiveness of the combined issues and not overload the reader with excessive details that could spoil the presentation of the chapter's main aim. The aim is to present the outstanding role of modeling, simulation, and results of their use, in a shortened way, in the contemporary research questions of the railway vehicle dynamics. However, this is going to be done through exploitation of the author's own results. The reference could be done here to [1], the comprehensive monograph in Polish (371 pages). This chapter is profiled

© 2015 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

differently, however, and results are presented in the concise way and internationally acces‐ sible form, i.e. in English.

## **1.1. The fundamentals in rail vehicle dynamics**

It is rather a well-known fact (e.g. Refs. [1–3]), however sometimes being forgotten, that in rail vehicle dynamics, perturbations of vehicle motion relative to vehicle general motion are of the primary interest. The general motion of vehicle is of lesser interest as it is predefined by the track alignment (shape). Hence, in the general view, the vehicle reproduces the motion imposed by the track centerline shape, as the vehicle is guided by the track (rails). Such a guided motion is already known. Instead, the perturbed motion relative to the track centerline appears to be of real importance in the rail vehicle dynamics problems.

Therefore, the description of such a perturbed motion makes the bases in the rail vehicle dynamics. As recently shown by the author in Ref. [1, 2], there are three options to perform such description.

**Option 1.** In this option, one adopts coordinates relative to inertial reference system of type *A* in Figure 1 and any formalism of equations building valid for motion relative to *A*:

$$fd\mathbf{B} = fd\mathbf{Z} \quad \Rightarrow \quad \{\mathbf{x} = \mathbf{x}(t), \dot{\mathbf{x}} = \dot{\mathbf{x}}(t)\} \quad \Rightarrow \quad \{\mathbf{x}' = \mathbf{x} - \mathbf{u}(t), \dot{\mathbf{x}}' = \dot{\mathbf{x}} - \mathbf{v}(t)\} \tag{1}$$

where *fd* is the operator representing any chosen formalism; B is the inertia forces in relation to *A*; Z is the external (active) forces; *x*, *x*˙ are the coordinates and velocities, respectively; and *u*(*t*), *v*(*t*), *w*(*t*) are the known functions of time representing displacements, velocities, and accelerations of transportation, that is, of system *A'* in relation to *A*.

**Option 2.** In this option, one adopts coordinates relative to noninertial reference system of type *A'* in Figure 1 and any formalism adapted to motion description relative to *A'* (the direct methods of relative motion dynamics):

$$f d \mathbf{B}' = f d \mathbf{Z} + f p \mathbf{P}(\dot{u}, \ddot{u}) \quad \Rightarrow \quad \{\mathbf{x}' = \mathbf{x}'(t), \dot{\mathbf{x}}' = \dot{\mathbf{x}}'(t)\} \tag{2}$$

where B' is the inertia forces in relation to *A'*; P is the imaginary forces (inertia forces arising from and dependent on the transportation); and *fp* is the operator of the imaginary forces, appropriate for the formalism adopted.

**Option 3.** In this option, one takes coordinates relative to *A'* and any formalism valid in *A*. In practice, this is any of the variational principles of mechanics or the formalism arising directly from it:

$$\begin{aligned} f \mathbf{z} \mathbf{B} &= f \mathbf{z} \mathbf{Z} \\ \mathbf{x} = \mathbf{x}' + \boldsymbol{\mu}(t) \\ \dot{\mathbf{x}} = \dot{\mathbf{x}}' + \boldsymbol{\nu}(t) \\ \ddot{\mathbf{x}} = \dot{\mathbf{x}}' + \boldsymbol{\nu}(t) \end{aligned} \implies \begin{aligned} \{\mathbf{x} = \mathbf{x}(t), \dot{\mathbf{x}} = \dot{\mathbf{x}}(t)\} \\ \{\mathbf{x}' = \mathbf{x}'(t), \dot{\mathbf{x}}' = \dot{\mathbf{x}}'(t)\} \\ \end{aligned} \tag{3}$$

where *fz* is the operator representing the variational principle.

differently, however, and results are presented in the concise way and internationally acces‐

It is rather a well-known fact (e.g. Refs. [1–3]), however sometimes being forgotten, that in rail vehicle dynamics, perturbations of vehicle motion relative to vehicle general motion are of the primary interest. The general motion of vehicle is of lesser interest as it is predefined by the track alignment (shape). Hence, in the general view, the vehicle reproduces the motion imposed by the track centerline shape, as the vehicle is guided by the track (rails). Such a guided motion is already known. Instead, the perturbed motion relative to the track centerline appears

Therefore, the description of such a perturbed motion makes the bases in the rail vehicle dynamics. As recently shown by the author in Ref. [1, 2], there are three options to perform

**Option 1.** In this option, one adopts coordinates relative to inertial reference system of type *A*

where *fd* is the operator representing any chosen formalism; B is the inertia forces in relation to *A*; Z is the external (active) forces; *x*, *x*˙ are the coordinates and velocities, respectively; and *u*(*t*), *v*(*t*), *w*(*t*) are the known functions of time representing displacements, velocities, and

**Option 2.** In this option, one adopts coordinates relative to noninertial reference system of type *A'* in Figure 1 and any formalism adapted to motion description relative to *A'* (the direct

where B' is the inertia forces in relation to *A'*; P is the imaginary forces (inertia forces arising from and dependent on the transportation); and *fp* is the operator of the imaginary forces,

**Option 3.** In this option, one takes coordinates relative to *A'* and any formalism valid in *A*. In practice, this is any of the variational principles of mechanics or the formalism arising directly

> ( ) { ( ), ( )} ( ) { ( ), ( )}

*x x ut x xt x xt x x vt x x t x x t*

<sup>ï</sup> = +¢ <sup>ï</sup> = = <sup>ý</sup> <sup>Þ</sup> = +¢ ¢¢ ¢¢ = = ï

& & & &

*fd fd x x t x x t x x u t x x v t* B Z { ( ), ( )} { ( ), ( )} = Þ = = Þ =- =- & & ¢ ¢& & (1)

*fd fd fp u u x x t x x t* Β Z P( , ) { ( ), ( )} ¢ =+ Þ = = & && ¢¢ ¢¢ & & (2)

& &

(3)

in Figure 1 and any formalism of equations building valid for motion relative to *A*:

sible form, i.e. in English.

such description.

**1.1. The fundamentals in rail vehicle dynamics**

166 Railway Research - Selected Topics on Development, Safety and Technology

to be of real importance in the rail vehicle dynamics problems.

accelerations of transportation, that is, of system *A'* in relation to *A*.

Β Z

= ü

*fz fz*

*x x wt*

&& &&

= +¢ ï

( )

þ

methods of relative motion dynamics):

appropriate for the formalism adopted.

from it:

**Figure 1.** Multibody system *S* and useful coordinate systems

Commenting on Eqns. (1–3), one should note that the left-hand side equalities represent dynamical equations of motion. In case of Eq. (3), the equations of motion are supplemented with kinematical relations. They express absolute variables with the relative ones and need to be introduced into Eq. (31) when equations of motion are being solved. The general form of Eqns. (1–3) serves any form of the vectorial and scalar equations as well as the matrix form in certain cases. Meaning of the dynamical equations of motion is typical. It becomes obvious as B represents inertia forces relative to absolute (inertial) system, Z represents external forces, B' represents inertia forces relative to moving (noninertial) system, and P represents imaginary forces (inertia forces depending on the transportation, i.e., motion of noninertial system relative to the inertial one). Moreover, the operators *fd*, *fp*, and *fz* do not change the meaning of the forces. The operators are introduced as a reminder that particular formalisms may need application of some specific coordinates, velocities, recording, and operations to make use of the formalism.

The curly brackets represent solutions of the equations of motion. As the equations are secondorder ordinary differential equations, the solutions are displacements (coordinates) and velocities, both the linear and angular ones in general.

It is worth emphasizing that Eq. (2) is the only one that needs explicit form of the imaginary forces to be recorded. In contrast, in case of Eqns. (1) and (3), the inertia forces are also taken into account, however, in an inexplicit way. This is because B=B' + P. In addition, if one uses the same formalism of building the equations to build Eq. (2) or to build either Eqns. (1) or (3), then form of the forces B and B' is identical. The only difference is with their meanings that refer to absolute and relative motions. In order to highlight this, the superscript " **'** " is used. It is used, in fact, to distinguish all the relative variables in Eqns. (1–3) and also in the subse‐ quent text.

Based on the denotation explained above, one can note in the curly brackets that Eq. (2) is the only one that leads directly to the solution representing the relative coordinates and velocities. Eq. (2) is advantageous as the relative variables are those of interest in most of rail-vehicle dynamics problems. In case of Eq. (1), the relative solution can be obtained indirectly, whereas in the case of Eq. (3) a part of the solution represents absolute variables and another the relative ones.

Discussing briefly the practical use of particular options, it has to be stated at first that numbers and types of application differ. The approach defined with Eq. (1) is definitely used most rarely. The approach defined with Eq. (3) is the one among all three that is used most often in the commercial software packages for automatic generation of equations of motion (AGEM). This includes the software suitable for the systems in a rail vehicle type. Currently, such software (e.g., MEDYNA, VAMPIRE, SIMPACK, or VI-RAIL codes) are used quite often mainly in solving many contemporary engineering problems. However, their use in scientific research is questionable. It is contested because of the so-called black box problem. Thus, following the formal scientific methodology, the researcher himself and the others should precisely know how the problem was resolved. Therefore, the assumptions and the way of modeling have to be elucidated, which allows others to repeat the study or to compare it with results for other similar but different approaches (for different assumptions and methods of modeling). Examples of use of this approach to build the models and the simulation software for the scientific purposes can also be found. However, for such purposes, the option defined with Eq. (2) seems to be used more frequently, where different levels of accuracy are used in practice. Thus, some authors neglect selected terms of the imaginary forces. The importance of these terms in the rail vehicle dynamics is comprehensively discussed by the author of the present chapter in Ref. [4].

## **2. Equations of motion based on dynamics of relative motion approach**

In his studies, the author of the present chapter himself practices and is a supporter of option 2 described earlier. This option defined with Eq. (2) exploits direct results of the dynamics of relative motion. Both, the variables used and the equations of motion, are defined directly in relation to moving coordinate systems in type *A'* in Figure 1. In his studies, the author applied a few formalisms of equation building. He applied them in both the traditional and numerical approaches. In the first approach, equations of motion were derived on paper first and then implemented in the simulation software in order to solve them numerically. In the second approach, the AGEM was applied by the author to build and solve equations numerically. Lagrange type II equations (e.g., see Ref. [5]) were in use in the traditional approach. In the AGEM approach, the author exploited Kane's equations as presented in Refs. [1, 2, 6]. Here, five examples of the equations valid in noninertial systems will be presented. These are matrix equations for a single rigid body based on Newton and Euler equations; state–space form of Newton–Euler equations for a single rigid body; state–space form of Newton–Euler equations for the multibody system (MBS) with constraints; Kane's equations in its general form; and Huston's form of Kane's equation suitable in AGEM. Besides, the method valid for Kane's equations will be presented, which enables building the equations for the holonomic and nonholonomic systems based on the equations for the free (unconstrained) system. Derivations for all these equations can be found in Refs. [1, 2].

#### **2.1. The Newton–Euler equations**

then form of the forces B and B' is identical. The only difference is with their meanings that refer to absolute and relative motions. In order to highlight this, the superscript " **'** " is used. It is used, in fact, to distinguish all the relative variables in Eqns. (1–3) and also in the subse‐

168 Railway Research - Selected Topics on Development, Safety and Technology

Based on the denotation explained above, one can note in the curly brackets that Eq. (2) is the only one that leads directly to the solution representing the relative coordinates and velocities. Eq. (2) is advantageous as the relative variables are those of interest in most of rail-vehicle dynamics problems. In case of Eq. (1), the relative solution can be obtained indirectly, whereas in the case of Eq. (3) a part of the solution represents absolute variables and another the relative

Discussing briefly the practical use of particular options, it has to be stated at first that numbers and types of application differ. The approach defined with Eq. (1) is definitely used most rarely. The approach defined with Eq. (3) is the one among all three that is used most often in the commercial software packages for automatic generation of equations of motion (AGEM). This includes the software suitable for the systems in a rail vehicle type. Currently, such software (e.g., MEDYNA, VAMPIRE, SIMPACK, or VI-RAIL codes) are used quite often mainly in solving many contemporary engineering problems. However, their use in scientific research is questionable. It is contested because of the so-called black box problem. Thus, following the formal scientific methodology, the researcher himself and the others should precisely know how the problem was resolved. Therefore, the assumptions and the way of modeling have to be elucidated, which allows others to repeat the study or to compare it with results for other similar but different approaches (for different assumptions and methods of modeling). Examples of use of this approach to build the models and the simulation software for the scientific purposes can also be found. However, for such purposes, the option defined with Eq. (2) seems to be used more frequently, where different levels of accuracy are used in practice. Thus, some authors neglect selected terms of the imaginary forces. The importance of these terms in the rail vehicle dynamics is comprehensively discussed by the author of the present

**2. Equations of motion based on dynamics of relative motion approach**

In his studies, the author of the present chapter himself practices and is a supporter of option 2 described earlier. This option defined with Eq. (2) exploits direct results of the dynamics of relative motion. Both, the variables used and the equations of motion, are defined directly in relation to moving coordinate systems in type *A'* in Figure 1. In his studies, the author applied a few formalisms of equation building. He applied them in both the traditional and numerical approaches. In the first approach, equations of motion were derived on paper first and then implemented in the simulation software in order to solve them numerically. In the second approach, the AGEM was applied by the author to build and solve equations numerically. Lagrange type II equations (e.g., see Ref. [5]) were in use in the traditional approach. In the AGEM approach, the author exploited Kane's equations as presented in Refs. [1, 2, 6]. Here,

quent text.

ones.

chapter in Ref. [4].

Based on the fundamental kinematical relations for relative motion, relating absolute and relative velocities and accelerations,

$$\mathbf{p}\_{\circ} = \mathbf{a}\_{\circ}^{\prime} + \mathbf{a}\_{\circ 1} + \mathfrak{s} \times \mathbf{r}\_{\circ}^{\prime} + \mathfrak{sp} \times (\mathfrak{sp} \times \mathbf{r}\_{\circ}^{\prime}) + \mathfrak{Loss} \times \mathbf{v}\_{\circ}^{\prime} \tag{4}$$

$$
\mathfrak{G} = \mathfrak{op} + \mathfrak{G}'\tag{5}
$$

$$d\mathbf{o} = \frac{d\mathbf{o}}{dt} = \frac{d\mathbf{o}}{dt} + \frac{d\mathbf{o}'}{dt} = \boldsymbol{\varepsilon} + \frac{d'\mathbf{o}'}{dt} + \boldsymbol{\mathfrak{op}} \times \boldsymbol{\mathfrak{op}'} = \boldsymbol{\mathfrak{e}} + \boldsymbol{\mathfrak{e}}' + \boldsymbol{\mathfrak{op}} \times \boldsymbol{\mathfrak{op}'} \tag{6}$$

and vectorial forms of Newton and Euler equations, the vectorial equations of relative motion for translation and rotation of a single free rigid body can be obtained as those presented in Refs. [7, 1, 2]. Based on these vectorial equations, their matrix forms can be recorded in several ways as shown, e.g., in Refs. [1, 2, 8–10]. The forms as in Refs. [1, 2] are as follows:

$$
\hat{\mathbf{a}} \cdot \mathbf{m} \mathbf{a}'\_{\mathcal{C}} = -m(\mathbf{a}\_{\boldsymbol{\phi}1} + \hat{\mathbf{r}}'\_{\mathcal{C}} \mathbf{a} + \hat{\mathbf{o}} \hat{\mathbf{o}} \mathbf{\hat{u}} r'\_{\mathcal{C}} + 2 \hat{\mathbf{o}} \mathbf{\hat{v}} r'\_{\mathcal{C}}) + \mathbf{R}\_{\mathcal{C}} = \mathbf{Q}\_1 + \mathbf{R}\_{\mathcal{C}} \tag{7}$$

$$\mathbf{J}\mathbf{g}' = -\hat{\mathbf{\Theta}}\mathbf{J}\mathbf{\Theta} - \mathbf{J}\mathbf{\mathcal{e}} - \mathbf{J}(\hat{\mathbf{\Theta}}\mathbf{\bullet}\mathbf{\mathcal{O}}') + \mathbf{T}\_{\mathbb{C}} = \mathbf{Q}\_2 + \mathbf{T}\_{\mathbb{C}} \tag{8}$$

The meanings of the denotations present in the above Eqns. (4–6) are as follows:*pC*, *θ*, *α* are the absolute acceleration of body mass centre *C*, absolute angular velocity of body *B*, and absolute angular acceleration of body *B*; *r* **′** *<sup>C</sup>* is the radius vector of *C* in *O*1*xyz* system; *v* **′** *<sup>C</sup>*, *a* **′** *C* are the relative velocity and relative acceleration of *C* in relation to *O*1*xyz*; *ω* **′** , *ε* **′** are the relative angular velocity and relative angular acceleration of *B*; *ao*1, *ω*, *ε* are the absolute acceleration of *O*1 origin (i.e., in relation to *OXYZ*), transportation angular velocity, and transportation angular acceleration (i.e., of *O*1*xyz* in *OXYZ*). In Eqs. (7, 8), *m*, **J** are the mass and inertia tensor (its matrix representation) of body *B* in relation to *C*; **R***C*, **T***C* are the matrix of resultant external forces acting on *B* and matrix of resultant external torques acting on *B* in relation to *C*; and **Q**1, **Q**2 are the matrices representing sums of the inertia terms. The remaining denotations are matrices that have their counterparts in the vectors defined in Eqns. (4–6).

One should realize here that in order to express all matrices in Eqns. (7, 8) through their elements selection of the vector bases *im* (*m*=1, 2, 3) for Eqns. (7) and (8) is necessary.

Equations (7, 8) can easily be recorded in one matrix equation:

$$\mathbf{I}\mathbf{i}\mathbf{i}'\_{\mathrm{ll}} = \mathbf{Q} + \mathbf{A} \tag{9}$$

where

$$\dot{\mathbf{x}}'\_{\rm II} = \begin{bmatrix} \dot{\mathbf{x}}'\_{\rm III} \end{bmatrix}^{\rm T} = \begin{bmatrix} \mathbf{a}'\_{\rm C}, \ \mathbf{g}'\_{\rm j} \end{bmatrix}^{\rm T} = \begin{bmatrix} \dot{\mathbf{v}}'\_{\rm C}, \ \dot{\mathbf{d}} \end{bmatrix}^{\rm T} \text{ (\$i = 1, \ldots, 6; j = 1, \ \mathbf{2}, \ \mathbf{3});}\tag{10}$$

$$\mathbf{I} = \begin{bmatrix} m\mathbf{E} & \mathbf{0} \\ \mathbf{0} & \mathbf{J} \end{bmatrix}; \qquad \mathbf{Q} = \begin{bmatrix} \mathbf{Q}\_1 \\ \mathbf{Q}\_2 \end{bmatrix}; \qquad \mathbf{A} = \begin{bmatrix} \mathbf{R}\_c \\ \mathbf{T}\_c \end{bmatrix} \tag{11}$$

and **E** is the unit matrix. Besides, for any vector *c* the denotation **c ⌢** represents the vector's skewsymmetric matrix, while generally the matrix **I** is not the symmetrical one.

On analyzing Eq. (10), the variables **x′** II and **x′** I (velocities and coordinates) are the same as those in Newton and Euler equations. They define translation of the centre *C* and rotation of *B* around *C*. In their case, the linear kinematical relation holds **x˙**I= **K**(**x**<sup>I</sup> )**x**II.

It can be shown that the equation similar to Eq. (9) can be obtained when arbitrarily chosen set of variables **w′** II and **w′ <sup>I</sup>** are adopted. If matrix **K** is not singular then for variables

$$\mathbf{x}'\_1 = \mathbf{X}\_{\mathsf{WI}}(\mathbf{w}'\_1, \mathbf{t}); \dot{\mathbf{x}}'\_1 = \mathbf{H}(\mathbf{w}'\_{1'}\mathbf{t})\dot{\mathbf{w}}'\_1 + \mathbf{h}(\mathbf{w}'\_{1'}\mathbf{t})\mathbf{t}$$

and linear kinematical relations

$$\dot{\mathbf{w}}'\_1 = \mathbf{G}(\mathbf{w}'\_1)\mathbf{w}'\_{\mathbb{II}}$$

the equation can be written as follows:

$$\mathbf{x}'\_{\rm II} = \mathbf{K}^{-1} \dot{\mathbf{x}}'\_{\rm I} = \mathbf{K}^{-1} (\mathbf{H} \dot{\mathbf{w}}'\_{\rm I} + \mathbf{h}) = (\mathbf{K}^{-1} \mathbf{H} \mathbf{G}) \mathbf{w}'\_{\rm II} + \mathbf{K}^{-1} \mathbf{h} = \mathbf{Q} \mathbf{h} \mathbf{w}'\_{\rm II} + \mathbf{g} \tag{12}$$

while the derivative equals

Modeling, Simulation, and Results of Their Use in Railway Vehicle Dynamics Studies http://dx.doi.org/10.5772/62105 171

$$\dot{\mathbf{x}}'\_{\rm II} = \mathbf{DA}\dot{\mathbf{w}}'\_{\rm II} + \mathbf{DA}\mathbf{w}'\_{\rm II} + \dot{\mathbf{g}} \tag{13}$$

After the introduction of Eq. (13) in (9), one gets

$$
\hat{\mathbf{I}}\dot{\mathbf{w}}'\_{\mathrm{II}} = \hat{\mathbf{Q}} + \hat{\mathbf{A}} \tag{14}
$$

where

**Q**2 are the matrices representing sums of the inertia terms. The remaining denotations are

One should realize here that in order to express all matrices in Eqns. (7, 8) through their

( ) <sup>T</sup> <sup>T</sup> ' T

1 2 ; ; *<sup>C</sup>*

in Newton and Euler equations. They define translation of the centre *C* and rotation of *B* around

It can be shown that the equation similar to Eq. (9) can be obtained when arbitrarily chosen set

I WI I I II I **x X w t x Hw tw hw t** ¢ ¢¢ ¢¢ ¢ = =+ ( , ); ( , ) ( , ) & &

<sup>I</sup> I II **w Gw w** & ¢ ¢¢ = ( )

1 1 1 1

**<sup>I</sup>** are adopted. If matrix **K** is not singular then for variables

II <sup>I</sup> <sup>I</sup> II II ( )( ) - - - - **x K x K Hw h K HG w K h w** ¢¢ ¢ ¢ ¢ = = += + = + & & W V (12)

é ù *<sup>m</sup>* é ù é ù = == ê ú ê ú ê ú ë û ë û ë û **E 0 Q R I Q 0 J <sup>Q</sup> <sup>T</sup>**

II IIi [] [ ] [ ] 1,...,6; 1, 2, 3 ; *Cj j Cj j* **xx a v** & & ¢ ¢ ¢¢ ¢ = e = & *,* w& *i,* = *j* == (10)

*C*

**⌢**

)**x**II.

I (velocities and coordinates) are the same as those

L (11)

represents the vector's skew-

II **Ix Q** &¢ = + L (9)

elements selection of the vector bases *im* (*m*=1, 2, 3) for Eqns. (7) and (8) is necessary.

matrices that have their counterparts in the vectors defined in Eqns. (4–6).

Equations (7, 8) can easily be recorded in one matrix equation:

170 Railway Research - Selected Topics on Development, Safety and Technology

and **E** is the unit matrix. Besides, for any vector *c* the denotation **c**

*C*. In their case, the linear kinematical relation holds **x˙**I= **K**(**x**<sup>I</sup>

On analyzing Eq. (10), the variables **x′**

II and **w′**

and linear kinematical relations

while the derivative equals

the equation can be written as follows:

of variables **w′**

symmetric matrix, while generally the matrix **I** is not the symmetrical one.

II and **x′**

where

$$
\hat{\mathbf{I}} = \mathbf{I}\mathbf{Q}; \quad \hat{\mathbf{Q}}\_{\parallel} = [\mathbf{Q} - \mathbf{I}(\dot{\mathbf{Q}}\mathbf{w}\_{\parallel}^{\prime} + \dot{\mathbf{q}})]; \quad \hat{\mathbf{A}} = \mathbf{A} \tag{15}
$$

As shown explicitly in Refs. [1, 2], the result represented by Eq. (14) can be generalized to the case of constraint system with holonomic and nonholonomic constraints (e.g., see Ref. [7]). Then it can be extended so that the inertia matrix becomes symmetrical one what enables to get the corresponding inverse matrix and finally to solve the equations. This can be done through the left-hand side multiplication of the equation by the transpose matrix **Ω**T. As a result of the described operations, one arrives at

$$
\mathbf{\bar{I}}\mathbf{\bar{w}}'\_{\mathrm{II}} = \mathbf{\bar{Q}} + \mathbf{\bar{A}} + \mathbf{\bar{A}}\_{\mathrm{z}} \tag{16}
$$

where

$$\bar{\mathbf{I}} = \mathbf{D}^{\mathrm{tr}} \hat{\mathbf{I}} = \mathbf{D}^{\mathrm{tr}} \hat{\mathbf{I}} \mathbf{D} \qquad ; \quad \bar{\mathbf{Q}} = \mathbf{D}^{\mathrm{tr}} \hat{\mathbf{Q}} ; \tag{17}$$

$$\bar{\mathbf{A}} = \mathbf{\bar{D}}^T \hat{\mathbf{A}}; \qquad \bar{\mathbf{A}}\_x = \mathbf{\bar{D}}^T \hat{\mathbf{A}}\_x = \mathbf{\bar{D}}^T \boldsymbol{\Phi}^T(\mathbf{w}\_{1'}', \mathbf{t}) \tag{18}$$

In the equations above, **t** is the indicator of the dependence on time *t*; **Φ, λ** are the so-called constraint matrix and column matrix of Lagrange's multipliers (e.g., see Refs. [8, 1, 2]). If one is going to solve Eq. (16) then all equations of constraints have to be added to the system of equations to make it possible.

When one is not interested in values of the constraint (internal) forces **Λ ⌣ <sup>z</sup>** then it is reasonable to express equations of motion for the reduced set of independent variables **y'I** and y'II. This reduces the number of equations by the number of constraints. The relation at the velocity level between maximum and reduced (minimum) set of the variables is represented through matrix representation of explicit constraint equations at the velocity level:

$$\mathbf{w}'\_{\mathrm{II}} = \mathbf{q} \mathbf{y}'\_{\mathrm{II}} + \mathbf{f} \tag{19}$$

After Eq. (19) is exploited in Eq. (14), the state–space equation in the independent generalized velocities takes the following form:

$$
\bar{\mathbf{I}} \dot{\mathbf{y}}'\_{\text{II}} = \bar{\mathbf{Q}} + \bar{\mathbf{A}} \tag{20}
$$

where

$$\tilde{\mathbf{I}} = \mathbf{q}^{\mathrm{T}} \hat{\mathbf{I}} \mathbf{q} \quad ; \qquad \tilde{\mathbf{Q}} = \mathbf{q}^{\mathrm{T}} [\hat{\mathbf{Q}} - \hat{\mathbf{I}}(\hat{\mathbf{q}} \mathbf{y}\_{\mathrm{II}}^{\prime} + \dot{\mathbf{q}})] \quad ; \qquad \tilde{\mathbf{A}} = \mathbf{q}^{\mathrm{T}} \hat{\mathbf{A}} \tag{21}$$

Note that inertia matrix **I ˜** is symmetrical due to the left-hand side multiplication by **φ**T. The kinematical relations corresponding to Eq. (20) in the linear matrix form can be expressed as **y˙′** I=**Y ^** I (**y′** I)**y′** II. In case nonholonomic constraints exist in the system, they have to be provided to make the solution of Eq. (20) possible.

The final Eqns. (16, 20) are valid for a single rigid body. It was shown in Ref. [8] for inertial systems that forms valid for a single body can directly be generalized to any number of rigid bodies. As in principle the structure of equations for noninertial systems differs in additional inertia terms of correction character only (see, e.g., Eqns. (2, 3)), this result can be extended to noninertial systems. In terms of the notation, it is trivial and means that Eqns. (16) and (20) valid in noninertial systems remain unchanged for any number of rigid bodies.

#### **2.2. Kane's equations**

Any of the formalisms of analytical mechanics valid in inertial systems can be adapted to describe the relative motion in the noninertial system. The author of the chapter performed such an adaptation for Kane's equation. It was done in a formal manner, that is, corresponding equations of relative motion were derived as shown in Refs. [6, 2, 1].

The partial velocities are fundamental to original Kane's approach [11]. The corresponding relative linear and angular partial velocities were introduced by the author as shown in Refs. [6, 1, 2]. Let us introduce them for the simple nonholonomic system [11] of *l* degrees of freedom composed of *j* (*j=*1*,...,n*) rigid bodies by defining linear *v* **′** *<sup>j</sup>* and angular *ω* **′** *<sup>j</sup>* relative velocities of the mechanical system:

$$\mathbf{v}'\_{j} = \sum\_{\rho=1}^{l} \frac{\partial \mathbf{v}'\_{j}}{\partial \mathbf{u}\_{j\rho}} \mathbf{u}\_{\rho} + \frac{\partial \mathbf{r}'\_{j}}{\partial \mathbf{t}} = \sum\_{\rho=1}^{l} \mathbf{v}'\_{|\rho} \boldsymbol{\mu}\_{\rho} + \mathbf{v}'\_{|t} = \sum\_{m=1}^{3} \left( \sum\_{\rho=1}^{l} \mathbf{v}'\_{|\rho m} \boldsymbol{\mu}\_{\rho} + \mathbf{v}'\_{|tm} \right) \mathbf{i}\_{m}^{(j)} \tag{22}$$

$$\mathbf{u}\mathbf{o}'\_{\rangle} = \sum\_{\rho=1}^{l} \frac{\partial \mathbf{o}'\_{\rangle}}{\partial \mathbf{u}\_{/\rho}} \mathbf{u}\_{\rho} + \frac{\partial \mathbf{o}'\_{\rangle}}{\partial \mathbf{t}} = \sum\_{\rho=1}^{l} \mathbf{o}\mathbf{o}'\_{/\rho} \mathbf{u}\_{\rho} + \mathbf{o}\mathbf{o}'\_{\boldsymbol{\beta}} = \sum\_{m=1}^{3} \left( \sum\_{\rho=1}^{l} \alpha'\_{/\rho m} \mathbf{u}\_{\rho} + \alpha'\_{/m} \right) \mathbf{i}^{(\boldsymbol{\beta})}\_{\boldsymbol{u}} \tag{23}$$

where *v* **′** *<sup>j</sup>ρ*, *ω* **′** *<sup>j</sup>ρ*, *v* **′** *jt*, *ω* **′** *jt* are the functions of generalized coordinates *q*1,...,*qk* and time *t*. At the same time, *v* **′** *<sup>j</sup>ρ*, *ω* **′** *<sup>j</sup>ρ* are the *ρ*th relative partial velocity of body *Bj* mass centre *Cj* (or any point *Pj* ) and *ρ*th relative angular partial velocity of *Bj* in relation to *A'* , respectively; *v* ′ *<sup>j</sup>ρ<sup>m</sup>* , *ω* ′ *jρm* are the scalar components of *v* **′** *<sup>j</sup>ρ*, *ω* **′** *<sup>j</sup>ρ*, called coefficients of partial velocities; *uρ* is the quasivelocity (*ρ* =1,...,*l*) ; *m* is the indicator of unit vectors (base vectors) *i<sup>m</sup>* ( *j*) (*m=*1, 2, 3) defining axes directions of reference systems, where *v* **′** *<sup>j</sup>* and *ω* **′** *<sup>j</sup>* are expressed (there is no need for *i<sup>m</sup>* ( *j*) to have the same directions for translations and rotations and for all bodies *Bj* ).

Form of the adapted Kane's equations is as follows:

After Eq. (19) is exploited in Eq. (14), the state–space equation in the independent generalized

T T T

kinematical relations corresponding to Eq. (20) in the linear matrix form can be expressed as

The final Eqns. (16, 20) are valid for a single rigid body. It was shown in Ref. [8] for inertial systems that forms valid for a single body can directly be generalized to any number of rigid bodies. As in principle the structure of equations for noninertial systems differs in additional inertia terms of correction character only (see, e.g., Eqns. (2, 3)), this result can be extended to noninertial systems. In terms of the notation, it is trivial and means that Eqns. (16) and (20)

Any of the formalisms of analytical mechanics valid in inertial systems can be adapted to describe the relative motion in the noninertial system. The author of the chapter performed such an adaptation for Kane's equation. It was done in a formal manner, that is, corresponding

The partial velocities are fundamental to original Kane's approach [11]. The corresponding relative linear and angular partial velocities were introduced by the author as shown in Refs. [6, 1, 2]. Let us introduce them for the simple nonholonomic system [11] of *l* degrees of freedom

> *ll l j j <sup>j</sup> j j jt j m jtm m*

 = = = ¢ ¢ æ ö ¢ = + = += + ¢¢ ¢ ¢ ç ÷ è ø å å åå *v r*

*u u vu v*

*ll l j j <sup>j</sup> j j jt j m jtm m*

w

*uu u*

 = = = ¢ ¢ æ ö ¢ = + = += + ¢¢ ¢ ¢ ç ÷ è ø å å åå *<sup>i</sup>*

valid in noninertial systems remain unchanged for any number of rigid bodies.

equations of relative motion were derived as shown in Refs. [6, 2, 1].

1 1 1 1

1 1 1 1

w

r r

*j m*

r r

*j m*

composed of *j* (*j=*1*,...,n*) rigid bodies by defining linear *v* **′**

*u t* r

 ¶

> ¶

 ¶

j

 ¶

rr

*u t* r

rr

r

r

¶

¶

w

= ¶

= ¶

II ˆ ˆ <sup>ˆ</sup> <sup>ˆ</sup> **I I Q QIy** = =- + = ; [ ( )] ; ¢ % & % % <sup>j</sup> <sup>j</sup> <sup>j</sup> <sup>j</sup>& <sup>x</sup> <sup>L</sup> <sup>j</sup> <sup>L</sup> (21)

II. In case nonholonomic constraints exist in the system, they have to be provided

**˜** is symmetrical due to the left-hand side multiplication by **φ**T. The

II **Iy Q** ¢ = + % % % & <sup>L</sup> (20)

*<sup>j</sup>* and angular *ω* **′**

<sup>3</sup> ( )

<sup>3</sup> ( )

 w

(23)

r r

w

r r

 r

> r

*v v v i* (22)

*<sup>j</sup>* relative velocities

velocities takes the following form:

172 Railway Research - Selected Topics on Development, Safety and Technology

Note that inertia matrix **I**

**2.2. Kane's equations**

of the mechanical system:

w

to make the solution of Eq. (20) possible.

where

**y˙′** I=**Y ^** I (**y′** I)**y′**

( ) ( ) ( ) ( ) 1 1 1 1 1 1 1 \* 1 1 2 2 0,5 *n n j jj j Cj j j n j j oj j j j j j j j j j j j n n j jj j j j j Cj j j n j jj j j j j j j j j n jj jj j j m mm m m* r r r r r r r r J = = = = = = = = ¢ ¢¢ ×- + × + ×- - ´ - ´ ´ - ´ ¢ ¢ ¢¢ <sup>é</sup> <sup>ù</sup> <sup>ë</sup> <sup>û</sup> + ×- × - ´ × + × × ¢ ¢¢ ¢ ¢ + ×- × - ´ × - ´ - × ¢ ¢ é ù ë û = ×+ × ¢ ¢ å å å å å å å **J J J J JE** *v a vR va r r v T vR vR* e w w w w e w w w w e w w w w \*\* 1 \* \*\* 1 11 0 ( 1,..., ) *n n j j j n nn j j jj j j j jj l* r rrr r = = == + × ¢ + × + ×+ × = = ¢ ¢¢ å å å åå *v R* w *TTT* w w (24)

Most of the denotations in Eq. (24) are as already defined. Here, the meanings in Eqns. (4–8) are helpful. Note, however, that in Eq. (24) the *j* rigid bodies are described, while in Eqns. (4– 8) just one. Besides, subscript *C* was omitted in *r* **′** *<sup>C</sup>* resulting in **r ′** *Cj* ≡**r ′** *<sup>j</sup>*. The undefined is *ϑ* – trace of **J**, that is, *ϑ* = *J*<sup>11</sup> + *J*<sup>22</sup> + *J*33. In addition, *R <sup>j</sup>* \* , *T <sup>j</sup>* \* are the inertia forces and torques, respec‐ tively, of motion in relation to noninertial systems of *A'* type; *R <sup>j</sup>* \*\*, *<sup>T</sup> <sup>j</sup>* \*\* are the resultant imaginary forces and torques, respectively; while *RCj* ≡ *R <sup>j</sup>* and *TCj* ≡*T <sup>j</sup>* .

Equation (24) can be expressed in the shortened matrix form as follows [6, 2, 1]:

$$\left(F\_{\rho}\right)^{\*} + F\_{\rho} + F\_{\rho}\stackrel{\*\*}{\rightharpoonup} = 0 \qquad\qquad\left(\rho = 1, \ldots, l\right)\tag{25}$$

where scalar components are defined with

$$\mathbf{F}\_{\rho}^{\prime} = \sum\_{j=1}^{u} (\mathbf{R}\_{\rangle}^{\prime} \cdot \mathbf{v}\_{\rangle\rho}^{\prime} + \mathbf{T}\_{\rangle}^{\prime} \cdot \mathbf{o} \mathbf{o}\_{\rangle\rho}^{\prime}) \ \ = \mathbf{R}\_{\langle u}^{\prime} \mathbf{o}\_{\rangle\rho u}^{\prime} + \mathbf{T}\_{\langle m}^{\prime} \mathbf{o}\_{\rangle\rho u}^{\prime} \tag{26}$$

$$F\_{\rho} = \sum\_{j=1}^{n} (\mathbf{R}\_{j} \cdot \boldsymbol{\sigma}\_{\boldsymbol{\rho}\rho}^{\prime} + \mathbf{T}\_{j} \cdot \boldsymbol{\mathfrak{op}}\_{\boldsymbol{\rho}\rho}^{\prime}) \quad = \mathbf{R}\_{jm} \boldsymbol{\sigma}\_{\boldsymbol{\rho}\rho u}^{\prime} + \mathbf{T}\_{jm} \boldsymbol{\alpha}\_{\boldsymbol{\rho}\rho u}^{\prime} \tag{27}$$

$$\mathbf{F}\_{\rho}^{\prime \ast} = \sum\_{j=1}^{n} (\mathbf{R}\_{j}^{\ast \ast} \cdot \mathbf{v}\_{\cdot \vert \rho}^{\prime} + \mathbf{T}\_{j}^{\ast \ast} \cdot \mathbf{a}\_{\cdot \vert \rho}^{\prime}) = \mathbf{R}\_{\vert m}^{\ast \ast} \mathbf{v}\_{\cdot \vert \rho m}^{\prime} + \mathbf{T}\_{\vert m}^{\ast \ast} \mathbf{o}\_{\cdot \vert \rho m}^{\prime} \tag{28}$$

In the above Eqns. (26–28), *F<sup>ρ</sup>* \* , *Fρ* i *F<sup>ρ</sup>* \*\* are the generalized inertia forces (relative to non-inertial system(s)), generalized external forces (identical with those in inertial system(s)), and gener‐ alized imaginary forces; *R jm* \* , *<sup>R</sup> jm*, *<sup>R</sup> jm* \*\* , *<sup>T</sup> jm* \* , *<sup>T</sup> jm* and *<sup>T</sup> jm* \*\* are the scalar components of resultant forces and torques *R <sup>j</sup>* \* , *R <sup>j</sup>* , *R <sup>j</sup>* \*\*, *<sup>T</sup> <sup>j</sup>* \* , *T <sup>j</sup>* and *T <sup>j</sup>* \*\*. The final forms of Eqns. (26–28) exploit results for inertial systems by Huston (e.g., see Refs. [12–14]). Also note that signs of sums over *j* and *m* are omitted in these final forms. The original Huston's convention is used here, in which summations over repeated indices have to be performed instead.

Equations (24–28) are not useful in AGEM approach as they are too close to the vectorial origin of Kane's equations. Nevertheless, final forms of Eqns. (26–28) can be treated as the initial ones for the process of deriving the useful form. The process is based on the Huston results [12– 14] for the inertial scleronomic systems and present author's extensions [6, 2, 1] to non-inertial and rheonomic systems. Consequently, invoking results from Refs. [6, 2, 1] one can write

$$
\Theta\_{\rho\mathfrak{a}} \dot{u}\_{\mathfrak{a}} = f\_{\rho} \qquad \qquad \text{ ( $\rho = 1, ..., l$ )}\tag{29}
$$

where

$$
\Theta\_{\rho\pi} = m\_{\rangle} \upsilon\_{\rangle\_{\rho m}}' \upsilon\_{\rangle\_{\pi m}}' + \mathcal{J}\_{\rangle \ast m} \alpha\_{\rangle \rho s}' \alpha\_{\rangle \pi m}' \tag{30}
$$

\*\* ( ) [ ( )] *j j m j m jsm j s j m rsm jsw j w jqr j m q j j m jtm jsm j s jtm rsm jsw j m j w jtr jqr jtw q jtr jtw f F F mv v u J u e J u u mv v J e J u u* r r r r pp r pp p rp r r rp p w w w ww w w w w w ww ww =+ - + + ¢¢ ¢¢ ¢ ¢¢ - + ¢¢ ¢¢ + ¢ ¢ ¢ ¢¢ ¢¢ + + & & & & (31)

Besides, in Refs. [6, 2, 1], the form of imaginary forces corresponding to Eqns. (30) and (31) was presented as follows:

$$\begin{split} & F^{\prime \prime}\_{\rho} = -m\_{\slash} a\_{\u1m} \upsilon^{\prime}\_{\slash \rho m} - e\_{\text{sum}} m\_{\slash} \alpha\_{\not\rho} \langle e\_{\text{uv}\alpha} \alpha\_{\not\rho} r^{\prime}\_{\not\rho} \rangle \upsilon^{\prime}\_{\not\rho m} \\ & - e\_{\text{snu}} J\_{\mid\text{sw}} \alpha\_{\not\rho} \alpha\_{\not\rho} \alpha^{\prime}\_{\not\rho m} - e\_{\text{snu}} m\_{\not\rho} \alpha\_{\not\rho} r^{\prime}\_{\not\rho m} \upsilon^{\prime}\_{\!\rho m} - J\_{\mid\text{sw}} \alpha\_{\not\rho} \alpha^{\prime}\_{\not\rho m} \\ & - 2e\_{\text{snu}} m\_{\slash} \alpha\_{\not\rho} \upsilon^{\prime}\_{\not\rho m} \upsilon^{\prime}\_{\!\rho m} \upsilon^{\prime}\_{\!\pi} \mu\_{\!\pi} - 2e\_{\text{snu}} (J\_{\mid\text{sw}} - 0.5 \mathcal{S}\_{\!\!b} E\_{\!\!w \!\!u}) \alpha\_{\not\rho} \alpha^{\prime}\_{\!\rho m} \alpha^{\prime}\_{\!\rho m} \upsilon\_{\!\!u} \\ & - 2e\_{\text{snu}} m\_{\not\rho} \alpha\_{\not\rho} \upsilon^{\prime}\_{\not\rho m} \upsilon^{\prime}\_{\!\rho m} - 2e\_{\text{snu}} (J\_{\mid\text{sw}} - 0.5 \mathcal{S}\_{\!\!b} E\_{\!\!u}) \alpha\_{\not\rho} \alpha^{\prime}\_{\!\rho m} \alpha^{\prime}\_{\!\nu} \end{split} \tag{32}$$

The supplements to Eq. (29) are kinematical relations:

1

1

\* , *Fρ* i *F<sup>ρ</sup>*

\* , *<sup>R</sup> jm*, *<sup>R</sup> jm*

summations over repeated indices have to be performed instead.

rp p r& = =

rp

*mv v J*


r

<sup>1</sup> ( )

*e J e m rv J*

=- - ¢ ¢ ¢

*F ma v e m e r v*

r

r pp

w

r

 r

w

w ww

*j o jm j m swm j js wrq jr jq j m rsm jsw jw jr j m swm j js jw j m jms js j m*


w

2 2 ( 0,5 ) 2 2 ( 0,5 )

*em vv e J E*

*swm j js j m jtw rsm jsw j sw*

 a

*swm j js j m j w rsm jsw j sw jw j m j r*

 w


*em vvu e J E u*

[

 r pp

*j*

=

*n*

*j*

=

r

r

\* , *R <sup>j</sup>* , *R <sup>j</sup>* \*\*, *<sup>T</sup> <sup>j</sup>* \* , *T <sup>j</sup>*

In the above Eqns. (26–28), *F<sup>ρ</sup>*

alized imaginary forces; *R jm*

forces and torques *R <sup>j</sup>*

where

r r r

presented as follows:

\*\*

r

*n*

174 Railway Research - Selected Topics on Development, Safety and Technology

( )

( )

\*\* , *<sup>T</sup> jm*

and *T <sup>j</sup>*

 r

 rr

*F Rv T*

\*\* \*\* \*\* \*\* \*\*

*F Rv T*

= × +× = + å *Rv T* ¢¢ ¢ ¢ w

*jj j j jm j m jm j m*

*j j j j jm j m jm j m*

system(s)), generalized external forces (identical with those in inertial system(s)), and gener‐

for inertial systems by Huston (e.g., see Refs. [12–14]). Also note that signs of sums over *j* and *m* are omitted in these final forms. The original Huston's convention is used here, in which

Equations (24–28) are not useful in AGEM approach as they are too close to the vectorial origin of Kane's equations. Nevertheless, final forms of Eqns. (26–28) can be treated as the initial ones for the process of deriving the useful form. The process is based on the Huston results [12– 14] for the inertial scleronomic systems and present author's extensions [6, 2, 1] to non-inertial and rheonomic systems. Consequently, invoking results from Refs. [6, 2, 1] one can write

*Θu f* ( 1,..., )*l*

*j j m j m jsm j s j m*

*f F F mv v u J u e J u u*

=+ - + + ¢¢ ¢¢ ¢ ¢¢

w w

w w

 r pp

 r

rp

*Θ mv v J*

r p

*j j m jtm jsm j s jtm*

& &

w w w

r

w wr p

*j j m j m jsm j s j m rsm jsw j w jqr j m q*

*rsm jsw j m j w jtr jqr jtw q jtr jtw*

 ww

 p

+ ¢ ¢ ¢ ¢¢ ¢¢ + +

\*\* ( )

*e J u u*

Besides, in Refs. [6, 2, 1], the form of imaginary forces corresponding to Eqns. (30) and (31) was

r


rr

J

J

a w

> ww w

 w

\* , *<sup>T</sup> jm* and *<sup>T</sup> jm*

= ×+× = + å *Rv T* ¢ ¢¢¢ w

 r

rrr

 r w

 w

\*\* are the generalized inertia forces (relative to non-inertial

(27)

(28)

\*\*. The final forms of Eqns. (26–28) exploit results

(29)

 rp

= + ¢¢ ¢¢ (30)

p

& & (31)

 ww

w ww

r pp

ww

r

 r (32)

¢ ¢ =

( )]

\*\* are the scalar components of resultant

$$\dot{q}\_{\sigma} = P\_{\epsilon p}(q\_{\uparrow}t)u\_{\rho} + p\_{\sigma}(q\_{\uparrow}t) \quad ; \qquad (\sigma = 1, \ldots, k), \ (\rho = 1, \ldots, l) \tag{33}$$

where in the above equations *ersm*, *eswm* are the permutation symbols; *k* is the number of generalized coordinates of the nonholonomic system; *Pσρ*, *pσ* are the functional matrix and vector (column matrix). Besides, the square bracket in Eq. (31) (being the extension to rheo‐ nomic systems) and *pσ* vanish for scleronomic systems.

Equations (24, 25, 29, 32, 33) serve nonholonomic, holonomic, and free systems. In Eqns. (25, 29, 32, 33), the difference between particular systems is taken into account by suitable values of partial velocities' coefficients *v* ′ *<sup>j</sup>ρm*, *ω* ′ *<sup>j</sup>ρm*. For holonomic systems *l*=*k*, while for free ones *l*=*k*=6*n*. In addition, equations for free and holonomic systems and kinematical relations are sufficient to be solved. For nonholonomic systems, equations of nonholonomic constraints have to be provided additionally to get the solution.

It was shown by the author in Refs. [6, 2, 1] for the rail vehicle systems, taking account of their moderate dimension, nonoccurrence of nonholonomic constraints, limited number of holo‐ nomic constraints, identical coordinates and transportation for each rigid body within vehicle model, that it can be reasonable to make use of equations (or particular type of forces) for free systems. Assuming holonomic and nonholonomic constraint equations as

$$
\mu\_{\sigma} = \sum\_{\rho=1}^{l} B\_{\sigma\rho} \mu\_{\rho} + D\_{\sigma} \qquad \qquad \{\sigma = l+1, \ldots, k\} \tag{34}
$$

$$
\mu\_s = \sum\_{\rho=1}^k A\_{s\rho} \mu\_\rho + \mathbb{C}\_s \qquad \qquad \text{(s = k+1,...,6m)}\tag{35}
$$

the following formula enables to build equations (or selected forces type) for constraint nonholonomic system based on the equations (or selected forces type) for free system [6, 2, 1]:

$$F\_{\rho} = \tilde{F}\_{\rho} + \sum\_{s=k+1}^{6n} \tilde{F}\_{s} A\_{s\rho} + \sum\_{\sigma=l+1}^{k} \tilde{F}\_{\rho} B\_{\sigma\rho} + \sum\_{\sigma=l+1}^{k} \sum\_{s=k+1}^{6n} \tilde{F}\_{s} A\_{s\rho} B\_{\sigma\rho} \tag{36} = 1,...,l} \tag{36}$$

If the system is holonomic then Eq. (34) vanishes, while Eq. (36) is reduced to just two first addends. In Eqns. (34–36), the following so far undefined notations appear: *uρ*, *uσ*, *us* are the independent, dependent through nonholonomic constraints, and dependent through holo‐ nomic constraints quasi-velocities, respectively; *Asρ*, *Bσρ*, *Cs*, and *Dσ* are the functional coeffi‐ cients depending on generalized coordinates *q*1,...,*qk* and time *t*; and *<sup>F</sup>* , *<sup>F</sup>*˜ are the representations of any type of generalized forces *F* \* , *F* , *F* \*\* as defined in Eqns. (25–28) or of their sum for constraint and free systems, respectively.

## **2.3. Discussion of the equations**

Taking account of Refs. [15, 16], both the Newton–Euler and Kane's equations are those most often used in AGEM. This is an interesting fact because major differences between both approaches exist. Therefore, the equations are worthy of discussion in this context.

It can be seen in Section 2.1 that building Eq. (16) and (20) that represent receipt useful in AGEM are operations on matrices. To perform these operations, the matrices have to be defined with their components. To do this, one needs to adopt base vectors *i<sup>m</sup>* ( *<sup>j</sup>*) for the chosen reference systems. Unfortunately, operation must be performed at the beginning of the equations building, that is, jet in Eqns. (7, 8). Despite simple forms of Eqns. (16, 20), the method of obtaining equations for a given mechanical system is not short. Every time one builds equa‐ tions, he has to choose vector bases, define matrices, and perform matrix calculations as described in Eqns. (12, 13, 15, 17, 18, 19, 21). These calculations are multiple, often in each of the mentioned equations. In the case of nonholonomic systems and interest in constraint forces, Eq. (16) is used, which together with constraint equations form differential-algebraic equation (DAEs) set. Such a set is difficult to solve in general case. On the other hand, when constraint forces are of no interest then a reduced number of ordinary differential equation (ODE) set is used (Eq. (20)).

Equations (29–33) form the receipt useful in AGEM of Kane's equations. Their forms might seem discouraging as they are complex, especially while comparing with Eqns. (16, 20). This impression is misleading as equations possess serious advantages as well. The advantage is the same form of either scalar (Eq. (24)) or matrix equations (Eqns. (29–33)) for free, holonomic, and nonholonomic systems. The most distinguishing advantage of Kane's equations, appear‐ ing just for this formalism, is the moment the components of vectors and tensors are defined through the vector bases selection. This is done at the very end of the equations building. This is possible since Eqns. (24, 29–33) originate directly from Eq. (24) based on vectors. Therefore, Eqns. (24, 29–33) are valid for any vector basis. An additional advantage is the type of Kane's equations that are always ODEs.

## **3. Example objects, nominal models, and numerical models**

## **3.1. The objects and their nominal models**

Equations of motion of mechanical system are referred to as its mathematical model. In order to make use of the general methods of building equations of motion as shown in Section 2, the nominal model of particular vehicle (object) has to be first determined. In case of railway vehicle, its nominal model projects its structure and selected physical features of the structure elements. The nominal models of two example objects will be presented below. Both are of British origin and possess relatively simple structure that is suitable for the basic character of the author's research. The first model [1, 5, 17] corresponds to 2-axle hsfv1 freight car and is shown in Figure 2. The second one [1, 18] corresponds to 4-axle MKIII passenger car and is shown in Figure 3.

Modeling, Simulation, and Results of Their Use in Railway Vehicle Dynamics Studies http://dx.doi.org/10.5772/62105 177

**Figure 2.** Nominal model of 2-axle freight car hsfv1 [1, 5, 17]

**2.3. Discussion of the equations**

used (Eq. (20)).

equations that are always ODEs.

shown in Figure 3.

**3.1. The objects and their nominal models**

Taking account of Refs. [15, 16], both the Newton–Euler and Kane's equations are those most often used in AGEM. This is an interesting fact because major differences between both

It can be seen in Section 2.1 that building Eq. (16) and (20) that represent receipt useful in AGEM are operations on matrices. To perform these operations, the matrices have to be defined with

systems. Unfortunately, operation must be performed at the beginning of the equations building, that is, jet in Eqns. (7, 8). Despite simple forms of Eqns. (16, 20), the method of obtaining equations for a given mechanical system is not short. Every time one builds equa‐ tions, he has to choose vector bases, define matrices, and perform matrix calculations as described in Eqns. (12, 13, 15, 17, 18, 19, 21). These calculations are multiple, often in each of the mentioned equations. In the case of nonholonomic systems and interest in constraint forces, Eq. (16) is used, which together with constraint equations form differential-algebraic equation (DAEs) set. Such a set is difficult to solve in general case. On the other hand, when constraint forces are of no interest then a reduced number of ordinary differential equation (ODE) set is

Equations (29–33) form the receipt useful in AGEM of Kane's equations. Their forms might seem discouraging as they are complex, especially while comparing with Eqns. (16, 20). This impression is misleading as equations possess serious advantages as well. The advantage is the same form of either scalar (Eq. (24)) or matrix equations (Eqns. (29–33)) for free, holonomic, and nonholonomic systems. The most distinguishing advantage of Kane's equations, appear‐ ing just for this formalism, is the moment the components of vectors and tensors are defined through the vector bases selection. This is done at the very end of the equations building. This is possible since Eqns. (24, 29–33) originate directly from Eq. (24) based on vectors. Therefore, Eqns. (24, 29–33) are valid for any vector basis. An additional advantage is the type of Kane's

Equations of motion of mechanical system are referred to as its mathematical model. In order to make use of the general methods of building equations of motion as shown in Section 2, the nominal model of particular vehicle (object) has to be first determined. In case of railway vehicle, its nominal model projects its structure and selected physical features of the structure elements. The nominal models of two example objects will be presented below. Both are of British origin and possess relatively simple structure that is suitable for the basic character of the author's research. The first model [1, 5, 17] corresponds to 2-axle hsfv1 freight car and is shown in Figure 2. The second one [1, 18] corresponds to 4-axle MKIII passenger car and is

**3. Example objects, nominal models, and numerical models**

( *<sup>j</sup>*) for the chosen reference

approaches exist. Therefore, the equations are worthy of discussion in this context.

their components. To do this, one needs to adopt base vectors *i<sup>m</sup>*

176 Railway Research - Selected Topics on Development, Safety and Technology

**Figure 3.** Nominal model of 4-axle passenger car MKIII [1, 18]

If the flexibility of track is going to be taken into consideration in the study, then some track nominal models should be adopted. In case one considers low-frequency dynamics of the vehicle motion (not more than 50 Hz) but not the higher-frequency phenomena in the track itself, then it is reasonable to build the track model that is composed of rigid bodies. Then the whole vehicle-track system model is a multibody model. Separate models for lateral (Figure 4) and vertical (Figure 5) directions as used by the author in his studies are presented below [1, 5, 17, 18].

**Figure 4.** Nominal model of track flexible laterally [1, 17, 18]

**Figure 5.** Nominal model of track flexible vertically [1, 17, 18]

The denotations used in Figures 2–5 are as follows: *m* is the mass; *Jξ*, *Jη*, *Jζ* are the inertia moments in longitudinal, lateral, and vertical directions, respectively; *k, c* are the stiffness and damping in the flexible elements in general; *kz , cz , kp , cp* are the stiffness and damping in the flexible elements of primary and secondary suspension, respectively; *x, y, z* are the linear displacements (coordinates) in longitudinal, lateral, and vertical directions, respective‐ ly (if used as indices, they indicate the corresponding directions); *ϕ, χ, ψ* are the rotations (angular coordinates) around axes of longitudinal, lateral, and vertical directions, respective‐ ly (if used as indices, they indicate the corresponding rotations); 2*a*, 2*b* are the vehicle or bogie base; *h, r*t are the vertical dimensions and wheel rolling radius*; N, T, W* are the normal, tangential, and load forces in wheel/rail contact, respectively; *γ* is the contact angle. The indices used indicate *p* as the passenger car body or right-hand side, *l* the left-hand side; *b* the bogie or freight car body; and *t* the track. Values of the abovementioned parameters of models can be found in Refs. [1, 5, 17].

The whole hsfv1 car-track system model has 18 degrees of freedom (DOFs). The whole MKIII car-track model has as many as 38 degrees of freedom.

#### **3.2. Numerical models**

When mathematical models derived traditionally are available, based on general methods of the equations building as described in Section 2 and adopted nominal models as, for example, those shown in Section 3.1, then they have to be converted to numerical models. When AGEM approach is used to build the equations, based again on the results of Section 2 and Section 3.1, then equations automatically form the numerical model. The numerical models are indispensable, since rail vehicle systems are multidimensional ones. In addition, their nonlin‐ ear versions are mainly used at present. The only way to solve dynamical differential equations of motion of such systems is numerical integration. In order to do this, the numerical model representing the mathematical model in a form understandable to computers must be built, which means that the equations have to be coded in some of the programming languages. In the case of this author, it is Fortran. Then, the equations are solved with use of one of the numerical methods of equation integration. The author uses Gear's method [19]. The software containing traditionally derived equations of motion or generating the equations automatically combined with the integration procedure is called the simulation software in railway vehicle dynamics.

**Figure 4.** Nominal model of track flexible laterally [1, 17, 18]

178 Railway Research - Selected Topics on Development, Safety and Technology

**Figure 5.** Nominal model of track flexible vertically [1, 17, 18]

models can be found in Refs. [1, 5, 17].

**3.2. Numerical models**

car-track model has as many as 38 degrees of freedom.

The denotations used in Figures 2–5 are as follows: *m* is the mass; *Jξ*, *Jη*, *Jζ* are the inertia moments in longitudinal, lateral, and vertical directions, respectively; *k, c* are the stiffness and damping in the flexible elements in general; *kz , cz , kp , cp* are the stiffness and damping in the flexible elements of primary and secondary suspension, respectively; *x, y, z* are the linear displacements (coordinates) in longitudinal, lateral, and vertical directions, respective‐ ly (if used as indices, they indicate the corresponding directions); *ϕ, χ, ψ* are the rotations (angular coordinates) around axes of longitudinal, lateral, and vertical directions, respective‐ ly (if used as indices, they indicate the corresponding rotations); 2*a*, 2*b* are the vehicle or bogie base; *h, r*t are the vertical dimensions and wheel rolling radius*; N, T, W* are the normal, tangential, and load forces in wheel/rail contact, respectively; *γ* is the contact angle. The indices used indicate *p* as the passenger car body or right-hand side, *l* the left-hand side; *b* the bogie or freight car body; and *t* the track. Values of the abovementioned parameters of

The whole hsfv1 car-track system model has 18 degrees of freedom (DOFs). The whole MKIII

When mathematical models derived traditionally are available, based on general methods of the equations building as described in Section 2 and adopted nominal models as, for example, In order to make use of the procedure, shortly described above, three additional major elements have to be provided. The first one is to adopt the model and then to build its numerical module to take account of relative kinematics arising from description in moving coordinate systems. Strictly speaking, the linear and angular velocities and accelerations of transportation, represented with functions *u*(*t*), *v*(*t*), *w*(*t*) in Section 1.1 and with vectors *ao*1, *ω*, *ε* in Section 2, have to be determined. The author of this paper worked out a general method of their determination, which is a numerical calculation-oriented method [2, 1, 20]. Its generality arises from any three-dimensional track shape acceptable in the method. The most important requirement in the method is the description of the three-dimensional curve by parametric equations, with its length as the parameter. In fact, the three-dimensional case corresponds to transition curves (TCs) with the superelevation ramps, while circular (regular) curves (CC) and straight track (ST) are two- and one-dimensional special cases, respectively. The compo‐ nents, expressed analytically, of the velocities and accelerations of transportation for several types of TCs, CC, and ST can be found in Refs. [1, 20]. The sources for the algorithm to numerically determine the components for polynomial TCs are Refs. [21, 1].

The next element is tangential contact forces calculation. These forces arise from wheel/rail relative slip (creepage). Consequently, the longitudinal, lateral, and spin creepages are the inputs in the existing methods of the tangential forces calculation. Generally, one can talk about the exact and simplified numerical methods of nonlinear tangential contact forces calculation. The problem with the exact methods, as, for example, Kalker's CONTACT program [22], lies in the very slow calculations for the simulation purposes of rail vehicle dynamics. Thus, simplified methods are in use in the simulation programs (numerical models). One of them is Kalker's FASTSIM program [23]. This is just the software used by the present author in his models and so in the studies. Usually, the adopted value of the friction coefficient, necessary for the program to run, in the author's studies is *μ*=0.3. Other possibilities of tangential contact forces calculation are discussed in Refs. [24–27].

The last among three elements to be discussed is the problem of nonlinear contact geometry. This author uses ArgeCare RSGEO program [28] to address this. The program resolves normal contact problems and purely geometrical ones. As a result, contact areas become known and geometrical variables in the contact. The most important geometrical variables are instanta‐ neous contact angles, rolling radii, and contact point positions on the wheel and rail for any lateral relative shift between wheel and rail. In addition, the program is capable of taking account of the influence of wheelset yaw angle on the parameters. Both the left- and right-hand side wheel and rail parameters are calculated in the same step. The results of ArgeCare RSGEO program (the author uses) are tabulated before simulations for a given wheel and rail pair of profiles as a function of the discrete lateral wheel/rail displacements and yaw angles. To get the contact parameters for any relative displacement and yaw rotation between wheel and rail, the tabulated data are linearly interpolated.

At the end, let us return to the issue already discussed at the beginning of Section 2. It is the way in which the equations of motion in the simulation software (numerical model) are introduced. Here, this refers to the traditionally derived and automatically generated equations of motion (AGEM approach). Now, we will provide some references correspond‐ ing to both approaches. If one is interested in samples of the equations traditionally de‐ rived by this author, then he could find them in Ref. [29]. If one is interested in details of the package ULYSSES and its core program TITAN co-built by this author, then he can find them in Refs. [1, 2, 6]. The software packages based on AGEM built by other authors are dis‐ cussed in Refs. [15, 16].

## **4. Example results of the selected simulation studies**

The author of the present chapter has exploited simulation in his studies for more than 20 years. There is no possibility to refer to all of them. Readers interested in earlier applications of simulation by the author can find reviews of the corresponding references in Ref. [20]. Similar review, however, also including current applications, is done in Ref. [1]. Total number of all author's applications runs to tens. Due to the limit of space, just two applications of simulation used in currently studied problems will briefly be discussed in the following sections. Their content will be focused more on the achievements, contributions, and final results than on the details of the methods or procedures used in those studies. Interested readers will be informed about the essential publications where the details can be found. These publications can also be treated as a main base source for the results discussed further. The secondary source is Ref. [1]. The applications in view concern the stability and kinematics issues.

Two other applications extensively studied at present are connected with TC shape optimi‐ zation with the use of simulation and optimization methods as well as with dynamics of railway vehicles in TCs at velocities around the critical *vn*. The example publications on these problems are Refs. [21, 30] and [31, 32], respectively.

#### **4.1. Results of simulation use in rail vehicle nonlinear lateral stability studies**

The serious interest of the author on the stability problems started with reference [17]. The most important results were presented many years later in Refs. [33, 34]. The fundamental reference to the methods used by this author, based on bifurcation approach and on usage of the simulation, is [35]. Recently, the remarkable contributions have been by H. True and his coauthors (e.g., see Refs. [36, 37]). Works by O. Polach (e.g., see Ref. [38]) are also interesting. Much more comprehensive and detailed literature review can be found in Refs. [1, 33, 34].

geometrical variables in the contact. The most important geometrical variables are instanta‐ neous contact angles, rolling radii, and contact point positions on the wheel and rail for any lateral relative shift between wheel and rail. In addition, the program is capable of taking account of the influence of wheelset yaw angle on the parameters. Both the left- and right-hand side wheel and rail parameters are calculated in the same step. The results of ArgeCare RSGEO program (the author uses) are tabulated before simulations for a given wheel and rail pair of profiles as a function of the discrete lateral wheel/rail displacements and yaw angles. To get the contact parameters for any relative displacement and yaw rotation between wheel and rail,

At the end, let us return to the issue already discussed at the beginning of Section 2. It is the way in which the equations of motion in the simulation software (numerical model) are introduced. Here, this refers to the traditionally derived and automatically generated equations of motion (AGEM approach). Now, we will provide some references correspond‐ ing to both approaches. If one is interested in samples of the equations traditionally de‐ rived by this author, then he could find them in Ref. [29]. If one is interested in details of the package ULYSSES and its core program TITAN co-built by this author, then he can find them in Refs. [1, 2, 6]. The software packages based on AGEM built by other authors are dis‐

The author of the present chapter has exploited simulation in his studies for more than 20 years. There is no possibility to refer to all of them. Readers interested in earlier applications of simulation by the author can find reviews of the corresponding references in Ref. [20]. Similar review, however, also including current applications, is done in Ref. [1]. Total number of all author's applications runs to tens. Due to the limit of space, just two applications of simulation used in currently studied problems will briefly be discussed in the following sections. Their content will be focused more on the achievements, contributions, and final results than on the details of the methods or procedures used in those studies. Interested readers will be informed about the essential publications where the details can be found. These publications can also be treated as a main base source for the results discussed further. The secondary source is Ref. [1]. The applications in view concern the stability and kinematics

Two other applications extensively studied at present are connected with TC shape optimi‐ zation with the use of simulation and optimization methods as well as with dynamics of railway vehicles in TCs at velocities around the critical *vn*. The example publications on these

The serious interest of the author on the stability problems started with reference [17]. The most important results were presented many years later in Refs. [33, 34]. The fundamental

**4.1. Results of simulation use in rail vehicle nonlinear lateral stability studies**

the tabulated data are linearly interpolated.

180 Railway Research - Selected Topics on Development, Safety and Technology

**4. Example results of the selected simulation studies**

problems are Refs. [21, 30] and [31, 32], respectively.

cussed in Refs. [15, 16].

issues.

The meaning of stability is fundamental in the discussed problem. This author could describe its meaning in most of his works as follows. The used method of nonlinear lateral stability analysis regards formal theory of stability. By contrast, it is not so formal as compared with the methods based strictly on the theory. Instead, one finds the tool to be more practical. The simplification is based on adoption of certain assumptions as well as on behavioral expecta‐ tions for the studied system. They arise from commonly existing knowledge about the systems in the type of rail vehicle. Therefore, one can skip formal adoption of some solution as the reference in the simplified approach. Moreover, he can skip introducing some perturbations into the system to examine whether the new solution stays close enough (in narrow vicinity) to the reference solution. Such a stay is required by the formal definition of the stability theory. Taking account of that, bifurcation plot building for the system is the main task in the method. However, formal verification whether the solutions enabling to build this plot are stable in completely formal sense is not such a task. The assumption exists in the method that any typical solution of rail vehicle system (either stationary or periodic) is stable. It might be accepted with care basing on understanding that periodic solutions in rail vehicle dynamics are self-exciting vibrations being governed through the wheel–rail tangential contact forces. Following that, the theory of self-exciting vibrations can be useful in order to predict/expect typical periodic behavior of the system. The adopted assumption enables to limit the number of simulations for different initial conditions but the same given velocity. Such simplified approach is described in Ref. [29]. In case of any doubts if obtained solution is stable or about the possibility of multiple solutions existence, more formal check for the stability is a must. Then, reasonably a denser sweeping over the initial conditions, to introduce the perturbation, is performed. Such a more accurate approach is called the extended analysis and is described in Ref. [30].

The author has been especially interested in the stability of rail vehicles in a curved track from the very beginning of his interest in stability problems. The studies were initiated as a result of the observed periodic solutions in a circular curve of track. They appear for velocities *v* greater than the nonlinear critical velocity *vn*, exactly as in a straight track. The only difference relative to ST case is asymmetry of the results in relation to track centerline. Results of such kind are shown in Figure 6. Based on that fundamental observation, the studies were carried out focused on the influence of chosen factors on the stability in CC.

They involved curve radii *R* from the small to large ones and ST as well, where *R*=∞. The principles of these nonlinear stability analyses were adopted from the analysis in a ST. However, the important supplements to ST analysis were necessary. These supplements are shown in Figure 7. They finally resulted in the formally presented method of the nonlinear stability analysis in a curved track (CC).

It was successively improved for years and recently published in Refs. [33, 34]. The essence of this method is included in Figure 8.

**Figure 6.** Stable stationary and periodic solutions for *R*=600 m and velocities *v* below and above *vn* [1]

The meaning of the symbols and acronyms present in Figure 8 is as follows: *v*c, *v*n, *v*s is the linear critical velocity, nonlinear critical velocity, velocity of calculation stop due to unbounded growth of the solution, respectively; sss, sps are the stable stationary and periodic solutions, respectively; uss, ups are the unstable stationary and periodic solutions, respectively; SNB, HSB are the saddle-node and Hopf's bifurcation point, respectively. All these notions are well known in the nonlinear stability studies in ST. The method consists in building a pair of bifurcation plots as in Figures 8a and 8b based on the results of simulation as shown in Figures 8c and 8d and obtained subsequently for whole range of velocity of motion *v* and given curve radius *R*. If one decides to include in both bifurcation plots results for whole range of *R*, then one obtains complex bifurcation plots. Such a pair of the complex bifurcation plots was named stability map in Ref. [39]. The example stability maps, as in Refs. [33, 1], are shown in Figures 9 and 10. Note that the unstable solutions are omitted in these figures as those are less important because they are not observable in the real systems.

**Figure 7.** Supplements to the nonlinear lateral stability analysis in ST, necessary in CC analysis [33, 1]

The meaning of the symbols and acronyms present in Figure 8 is as follows: *v*c, *v*n, *v*s is the linear critical velocity, nonlinear critical velocity, velocity of calculation stop due to unbounded growth of the solution, respectively; sss, sps are the stable stationary and periodic solutions, respectively; uss, ups are the unstable stationary and periodic solutions, respectively; SNB, HSB are the saddle-node and Hopf's bifurcation point, respectively. All these notions are well known in the nonlinear stability studies in ST. The method consists in building a pair of bifurcation plots as in Figures 8a and 8b based on the results of simulation as shown in Figures 8c and 8d and obtained subsequently for whole range of velocity of motion *v* and given curve radius *R*. If one decides to include in both bifurcation plots results for whole range of *R*, then one obtains complex bifurcation plots. Such a pair of the complex bifurcation plots was named stability map in Ref. [39]. The example stability maps, as in Refs. [33, 1], are shown in Figures 9 and 10. Note that the unstable solutions are omitted in these figures as those are less important

**Figure 6.** Stable stationary and periodic solutions for *R*=600 m and velocities *v* below and above *vn* [1]

because they are not observable in the real systems.

182 Railway Research - Selected Topics on Development, Safety and Technology

**Figure 8.** The essence of method of stability analysis in a CC, i.e., building the bifurcation plot [34, 1]

**Figure 9.** Stability map for hsfv1 car for S1002/UIC60 wheel and rail profiles: (a) leading wheelset lateral displacements *ylw* and (b) peak to peak values of displacements *ylw*

Modeling, Simulation, and Results of Their Use in Railway Vehicle Dynamics Studies http://dx.doi.org/10.5772/62105 185

**Figure 10.** Stability map for hsfv1 car for SZDwheel/R65 wheel and rail profiles: (a) leading wheelset lateral displace‐ ment absolute values *ylw* and (b) peak to peak values of displacements *ylw*

The influence of several factors on the stability was studied based on the stability map technique. Among the factors studied are accuracy of wheelset's angle of attack determination, track superelevation, types of nominal wheel and rail profiles, type and value of the wheel and rail wear, stiffness and damping values in suspension, type of vehicle (car and bogie), rail inclination, way of mean rolling radius modeling, track gauge, and value of coefficient of friction. Figures 9 and 10 can represent part of the results for one of the factors, namely the wheel profile shape. Both maps were obtained for the same object, however, with different nominal wheel profiles. One can observe significant differences between both maps. One among the most important differences are nonlinear critical velocities *v*n of 43 and 24.1 m/s, respectively. Another is the existence of quasi-static solutions for *R*<2150 m in Figure 10. This can be easily recognized as no vibrations appear, resulting in p-t-p *y*lw=0, for such values of *R* in Figure 10b. In Figure 9, vibrations appear for all *R* values represented there. More differences are discussed in Refs. [33, 1].

The results in Figures 9 and 10 refer to 2-axle freight car hsfv1 of the model shown in Figure 2. Example results of stability studies for the model of MKIII passenger car shown in Figure 3 can be found in Ref. [18].

#### **4.2. Use of the simulation in studying the influence of kinematics accuracy on vehicle dynamics in a curved track at variable velocity**

Use of Option 2 from Section 1.1 to build mathematical and numerical models of vehicle-track systems by the present author, made it natural that he was always interested in the importance of accurate modeling relative kinematics connected with description in moving reference systems. This interest was amplified since an additional work is necessary as compared with description in absolute reference systems. Besides, some of the authors neglect additional terms (imaginary forces) in their equations of motion, without proper justification. Good examples of such works might be provided in Refs. [3, 40–42]. Author of this book chapter undertook two stage attempt to finally resolve the problem of particular inertia terms impor‐

**Figure 9.** Stability map for hsfv1 car for S1002/UIC60 wheel and rail profiles: (a) leading wheelset lateral displacements

**Figure 8.** The essence of method of stability analysis in a CC, i.e., building the bifurcation plot [34, 1]

184 Railway Research - Selected Topics on Development, Safety and Technology

*ylw* and (b) peak to peak values of displacements *ylw*

tance. The first one concerned vehicle motion with constant velocity. It finished with publica‐ tion [5]. The second attempt concerned the motion with variable velocity and finished with publication [4]. Publication [4] is also verification of the results from [5], as motion with constant velocity was just a special case in Ref. [4]. Therefore, Ref. [4] is the main source for results forthcoming in the current section. These results were also published in Ref. [1]. Both Refs. [4] and [1] present the issue in a comprehensive form. Here, only samples of the results and most general conclusions will be represented.

The idea of the study was to compare the results for the vehicle model with all imaginary forces included with those for the model in which the imaginary forces were omitted. In order to precisely determine importance of particular terms types, the forces (and torques) were selectively, rather than totally, omitted. Let us now recognize the generalized imaginary forces according to Ref. [7]. Here, equation (24) will be useful. The forces' terms are present in the second line of this equation and the torques' terms are present in the fourth line of it. Note that all the terms are multiplied by the corresponding linear and angular partial velocities. Therefore, we can refer directly to content of the square brackets in lines 2 and 4. Making use of that idea, the 1st term in the square bracket in line 2 represents inertia forces of translation. The 2nd term in line 2 and 1st one in line 4 form inertia forces of rotation. The 3rd term in line 2 and 2nd one in line 4 correspond to centrifugal forces of inertia. The 4th term in line 2 and the 3rd in line 4 make gyroscopic forces. The abovementioned four categories of imaginary forces' terms appear in the equations as the terms' components corresponding to the direction a particular equation of motion describes. Finally, the omissions of all components for given directions were performed. In fact, these were longitudinal, lateral, and vertical directions. At the end, for the given direction, the analyses were performed to find the components of the terms of particular (practical) importance among all.

The variable velocity was realized in the studies with use of the uniform variable motion. The change in velocity is indirectly represented by the corresponding acceleration value *a*. Consequently, negative values of acceleration *a* mean deceleration, corresponding to vehicle braking, the positive values of *a* mean vehicle acceleration, corresponding to its speeding up, while zero value of *a* means motion with constant velocity, that is, *v*=const. Different intensity of braking or speeding up were realized with different acceleration values.

## *4.2.1. Scope of the studies and example results of simulation*

Generally, the tests described in Ref. [4] represent seven different routes and two different vehicles. The routes included radii *R*=300, 600, 1200, and 2000 m and were composed in continuation either as ST, TC, and CC or as ST, TC, CC, TC, and ST. Thus, the entrance TCen and exit TCex transition curves were considered. The TCs were of the 3rd degree parabola type. Characterizing all the routes briefly, it has to be stated that their two general types existed. The first type corresponded to real railway conditions. The second one represented the case with unnaturally shortened TC in order to get stronger effects for research purposes. As concerns the two vehicles used in Ref. [4], the first one was hsfv1 car described in Section 3.1. Second was vehicle referred to as the 2-axle freight car in Ref. [5]. Both models have the same structure as shown in Figure 2, and they are supplemented with the same track models from Figures 3 and 4. This also results in the same number of DOFs for both of them. The major difference is that first model represents empty car, while the second one the laden car. The other minor differences in the parameters (see [5]) are negligible, here.

tance. The first one concerned vehicle motion with constant velocity. It finished with publica‐ tion [5]. The second attempt concerned the motion with variable velocity and finished with publication [4]. Publication [4] is also verification of the results from [5], as motion with constant velocity was just a special case in Ref. [4]. Therefore, Ref. [4] is the main source for results forthcoming in the current section. These results were also published in Ref. [1]. Both Refs. [4] and [1] present the issue in a comprehensive form. Here, only samples of the results

The idea of the study was to compare the results for the vehicle model with all imaginary forces included with those for the model in which the imaginary forces were omitted. In order to precisely determine importance of particular terms types, the forces (and torques) were selectively, rather than totally, omitted. Let us now recognize the generalized imaginary forces according to Ref. [7]. Here, equation (24) will be useful. The forces' terms are present in the second line of this equation and the torques' terms are present in the fourth line of it. Note that all the terms are multiplied by the corresponding linear and angular partial velocities. Therefore, we can refer directly to content of the square brackets in lines 2 and 4. Making use of that idea, the 1st term in the square bracket in line 2 represents inertia forces of translation. The 2nd term in line 2 and 1st one in line 4 form inertia forces of rotation. The 3rd term in line 2 and 2nd one in line 4 correspond to centrifugal forces of inertia. The 4th term in line 2 and the 3rd in line 4 make gyroscopic forces. The abovementioned four categories of imaginary forces' terms appear in the equations as the terms' components corresponding to the direction a particular equation of motion describes. Finally, the omissions of all components for given directions were performed. In fact, these were longitudinal, lateral, and vertical directions. At the end, for the given direction, the analyses were performed to find the components of the

The variable velocity was realized in the studies with use of the uniform variable motion. The change in velocity is indirectly represented by the corresponding acceleration value *a*. Consequently, negative values of acceleration *a* mean deceleration, corresponding to vehicle braking, the positive values of *a* mean vehicle acceleration, corresponding to its speeding up, while zero value of *a* means motion with constant velocity, that is, *v*=const. Different intensity

Generally, the tests described in Ref. [4] represent seven different routes and two different vehicles. The routes included radii *R*=300, 600, 1200, and 2000 m and were composed in continuation either as ST, TC, and CC or as ST, TC, CC, TC, and ST. Thus, the entrance TCen and exit TCex transition curves were considered. The TCs were of the 3rd degree parabola type. Characterizing all the routes briefly, it has to be stated that their two general types existed. The first type corresponded to real railway conditions. The second one represented the case with unnaturally shortened TC in order to get stronger effects for research purposes. As concerns the two vehicles used in Ref. [4], the first one was hsfv1 car described in Section 3.1. Second was vehicle referred to as the 2-axle freight car in Ref. [5]. Both models have the same structure as shown in Figure 2, and they are supplemented with the same track models from Figures 3 and 4. This also results in the same number of DOFs for both of them. The major

of braking or speeding up were realized with different acceleration values.

and most general conclusions will be represented.

186 Railway Research - Selected Topics on Development, Safety and Technology

terms of particular (practical) importance among all.

*4.2.1. Scope of the studies and example results of simulation*

Two types of simulation results are presented below. The first type is imaginary forces being omitted in the study. They are a cause for the solution differences being of interest. The imaginary forces are shown in Figures 11, 14, and 16. The second type of simulation results is those representing the differences between solutions for the whole model and the model with particular imaginary force omitted. They are represented with Figures 12, 13, 15, and 17. Results for complete model are drawn with the solid line, while for the model with imaginary torques omitted with the dashed line.

Just two routes appear in the mentioned figures. The route I is as follows: ST (*l*=10 m, *H*=0 m); TCen (*l*=7 m, *R*min=300 m, *H*max=0.15 m); CC (*l*=30 m, *R*=300 m, *H*=0.15 m). Initial velocities for route I were between *v*0=10.25 and 15.0 m/s, where the smallest value corresponded to *a*=0 m/s2 and the biggest one to *a*=6 m/s2 . Route I includes TC with length unnaturally shortened for research purposes. Route II is as follows: ST (*l*=10 m, *H*=0 m); TCen (*l*=102, 4 m, *R*min=600 m, *H*max=0.128 m); CC (*l*=60 m, *R*=600 m, *H*=0.128 m); TCex (*l*=102.4 m, *R*min=600 m, *H*max=0.128 m); ST (*l*=40 m, *H*=0 m). Initial velocities for the simulations were *v*0=19.91, 22.0, and 22.09 m/s. They corresponded to *a*=2, 0, and −2 m/s<sup>2</sup> , respectively. Route II entirely corresponds to real railway conditions. In the above given brackets, the denotations are *l* is denoted as the length of the track section, *R* as the circular curve radius, and *H* as the circular curve superelevation. Values of the corresponding accelerations are given directly in the figures in round brackets.

It is worth noting in Figure 11 that it includes the biggest values of the imaginary torque obtained in the author's studies. As is seen, it is the roll imaginary torque *Pϕb* of hsfv1 car body related to longitudinal direction. Character of the torque changes is rectangular, that is, nonzero values of the torque exist for TC only. Values of this torque depend strongly on value of the acceleration *a*. The value is biggest for the biggest acceleration *a=6 m/s*<sup>2</sup> and smallest as well as equal to zero for *a*=0 *m/s*<sup>2</sup> . The last fact means no influence of this torque on motion with constant velocity *v*=const. In the case of Figures 12 and 13, strong influence of omission of torque *Pϕb* from Figure 11 can be observed. The two mostly influenced coordinates are shown. They are roll *ϕ*b and yaw *ψ*b angles of hsfv1 car body, respectively. Note that the influence on angle *ψ*b (Figure 13) is approximately one order of magnitude smaller than on angle *ϕ*b (Figure 12).

**Figure 11.** Roll imaginary torque of hsfv1 car body on route I for different accelerations

**Figure 12.** Roll angle of hsfv1car body on route I for different accelerations and roll imaginary torque of hsfv1car body omitted

**Figure 13.** Yaw angle of hsfv1car body on route I for different accelerations and roll imaginary torque of hsfv1car body omitted

**Figure 14.** Yaw imaginary torque of hsfv1 car body on route I for different accelerations

Modeling, Simulation, and Results of Their Use in Railway Vehicle Dynamics Studies http://dx.doi.org/10.5772/62105 189

**Figure 15.** Yaw angle of hsfv1car body on route I for different accelerations and yaw imaginary torque of hsfv1car body omitted

**Figure 16.** Roll imaginary torque of 2-axle freight car body on route II for different accelerations

**Figure 12.** Roll angle of hsfv1car body on route I for different accelerations and roll imaginary torque of hsfv1car body

188 Railway Research - Selected Topics on Development, Safety and Technology

**Figure 13.** Yaw angle of hsfv1car body on route I for different accelerations and roll imaginary torque of hsfv1car body

**Figure 14.** Yaw imaginary torque of hsfv1 car body on route I for different accelerations

omitted

omitted

**Figure 17.** Lateral displacement of 2-axle freight car body on route II for different accelerations and roll imaginary tor‐ que of 2-axle freight car body omitted

In Figure 14, the value of the imaginary torque *Pψb* of hsfv1 car body, related to vertical direction, is represented. Character of the torque changes corresponds to rectangle supple‐ mented with triangle. The rectangle part is independent of the acceleration *a*. Therefore, it exists also for *a*=0 *m/s*<sup>2</sup> *.* Increase of the supplementary triangle part grows with increase of *a.* Therefore, the torque value is biggest for the biggest acceleration *a=6 m/s*<sup>2</sup> and smallest for *a*=0 *m/s*<sup>2</sup> . Overall, this means that in case of *v*=const. and *a*=0 *m/s*<sup>2</sup> influence of *Pψb* omission exists; however, for *a*>0, it is bigger. This is also the case in Figure 15 for *ψb*, where very strong influence of the omission can be seen. The biggest influence for *a*=6 m/s2 is of the same order of magnitude as in Figure 13 related to *Pϕb* omission. Therefore, it can be concluded that influences of *Pϕb* and *Pψb* omissions on *ψb* values accumulate.

In Figure 16, the character of torque *Pϕb* for 2-axle freight car is in accordance with that in Figure 11. The difference exists, however, that results from route II configuration. One can observe opposite signs of *Pϕb* for the entrance TC (TCen) and exit TC (TCex). The sign becomes also opposite when opposite acceleration *a* sign is adopted. In Figure 16, this case corresponds to *a*=2 m/s2 and *a*=−2 m/s2 . In Figure 17, for 2-axle freight car very strong influence of *Pϕb* omission on vehicle body lateral displacements *yb* can be observed. The changes in signs for *y<sup>b</sup>* corre‐ spond to those for the *Pϕb* described above provided one remembers that if *Pϕb*>0 then *yb*<0 and vice versa.

#### *4.2.2. General conclusions from the study*

Except the important influences described in Section 4.2.1, some other important influences were found and indicated in Ref. [4]. Therefore, except influences of *Pϕb* on *ϕb*, *yb*, and *ψb*, the important influence on the wheelsets' lateral displacements *ylw* and *ytw* exists. Except the influence of *Pψb* on *ψb*, very serious influence on the wheelsets yaw angles *ψlw* and *ψtw* exists. Assuming that lateral component (parallel to track plane in curved track) is known for its importance, the vertical component (perpendicular to the track plane) was studied. It appeared that the influence of this component omission on the results, especially vehicle body vertical displacements *zb*, is also significant.

Conclusions from the studies, including track section, character of velocity, vehicle elements, direction defining equation, and purpose of calculations, were summarized in table form in Ref. [4]. Here, we present them in Table 1. This table is analogous to that in Ref. [4] in terms of merits. Nevertheless, Table 1 is differently arranged.

Generally, the most important are longitudinal and vertical direction terms related to car body. The influence of eventual term omission is particularly important in TCs. The centrifugal force is important in both the TCs and the CC sections. The influence in TCs increases as compared with CC when the vehicle moves with variable velocity. The stronger the change in velocity (both accelerating and braking), the greater the influence. The last column represents a direct recommendation of terms that should not be neglected.


**Table 1.** Importance and recommendations for imaginary forces omission in rail vehicle dynamics

## **5. Concluding remarks**

In Figure 14, the value of the imaginary torque *Pψb* of hsfv1 car body, related to vertical direction, is represented. Character of the torque changes corresponds to rectangle supple‐ mented with triangle. The rectangle part is independent of the acceleration *a*. Therefore, it

Therefore, the torque value is biggest for the biggest acceleration *a=6 m/s*<sup>2</sup> and smallest for *a*=0

however, for *a*>0, it is bigger. This is also the case in Figure 15 for *ψb*, where very strong influence

as in Figure 13 related to *Pϕb* omission. Therefore, it can be concluded that influences of *Pϕb* and

In Figure 16, the character of torque *Pϕb* for 2-axle freight car is in accordance with that in Figure 11. The difference exists, however, that results from route II configuration. One can observe opposite signs of *Pϕb* for the entrance TC (TCen) and exit TC (TCex). The sign becomes also opposite when opposite acceleration *a* sign is adopted. In Figure 16, this case corresponds to

on vehicle body lateral displacements *yb* can be observed. The changes in signs for *y<sup>b</sup>* corre‐ spond to those for the *Pϕb* described above provided one remembers that if *Pϕb*>0 then *yb*<0 and

Except the important influences described in Section 4.2.1, some other important influences were found and indicated in Ref. [4]. Therefore, except influences of *Pϕb* on *ϕb*, *yb*, and *ψb*, the important influence on the wheelsets' lateral displacements *ylw* and *ytw* exists. Except the influence of *Pψb* on *ψb*, very serious influence on the wheelsets yaw angles *ψlw* and *ψtw* exists. Assuming that lateral component (parallel to track plane in curved track) is known for its importance, the vertical component (perpendicular to the track plane) was studied. It appeared that the influence of this component omission on the results, especially vehicle body vertical

Conclusions from the studies, including track section, character of velocity, vehicle elements, direction defining equation, and purpose of calculations, were summarized in table form in Ref. [4]. Here, we present them in Table 1. This table is analogous to that in Ref. [4] in terms of

Generally, the most important are longitudinal and vertical direction terms related to car body. The influence of eventual term omission is particularly important in TCs. The centrifugal force is important in both the TCs and the CC sections. The influence in TCs increases as compared with CC when the vehicle moves with variable velocity. The stronger the change in velocity (both accelerating and braking), the greater the influence. The last column represents a direct

. Overall, this means that in case of *v*=const. and *a*=0 *m/s*<sup>2</sup>

of the omission can be seen. The biggest influence for *a*=6 m/s2

190 Railway Research - Selected Topics on Development, Safety and Technology

*Pψb* omissions on *ψb* values accumulate.

and *a*=−2 m/s2

*4.2.2. General conclusions from the study*

displacements *zb*, is also significant.

merits. Nevertheless, Table 1 is differently arranged.

recommendation of terms that should not be neglected.

*.* Increase of the supplementary triangle part grows with increase of *a.*

. In Figure 17, for 2-axle freight car very strong influence of *Pϕb* omission

influence of *Pψb* omission exists;

is of the same order of magnitude

exists also for *a*=0 *m/s*<sup>2</sup>

*m/s*<sup>2</sup>

*a*=2 m/s2

vice versa.

The author of this chapter has been involved in technological problems of railways for about 30 years. More precisely, his involvement concerns part of the railway system that are railway vehicles. Within the railway vehicles, many areas of studies can be mentioned as, for example, material issues, production technical issues, braking systems issues, traction systems issues, construction issues, exploitation issues, control issues, safety systems issues, etc. The author is interested in dynamics of vehicle motion with a special regard to the motion in curved sections of the track. It is obvious that the chapter was devoted to these last aspects. The author focused his efforts on giving the idea to the readers how important and powerful tools are methods of modeling and simulation when used for the purposes of the rail vehicle dynamics. In Sections 1 and 2, the general methods of modeling dynamics in moving coordinate systems, useful for multibody systems of rail vehicle type, were represented. The Newton–Euler and Kane's equations in the form suitable in AGEM approach were discussed. At the end of Section 2.2, the method profiled for rail vehicle systems was also presented. These methods help build mathematical models of rail vehicle dynamics. In Section 3, example nominal models of rail vehicles were introduced, as used by the author. In that section, the issue of building numerical models corresponding to the mathematical ones was also discussed.

In Section 4, two examples of the author's simulation studies, with use of the numerical models, were discussed. The first example represented in Section 4.1 concerned the problem of rail vehicle stability in a curved track. Due to the author's systematic and consequent simulation studies, the method was formulated suitable for rail vehicle nonlinear lateral stability studies in a curved section of the track. The examples of the so-called stability map, being the most important result in the studies, were shown and briefly commented on. The stability maps approach was used by this author to study the influence of different factors on the results of stability analysis in a curved track.

The second example is represented in Section 4.2. It was devoted to the use of the simulation in studying the influence of kinematics accuracy on vehicle dynamics in a curved track at variable velocities. In fact, both the constant and variable velocity cases were of interest. Due to the consequent author's interest in this issue, it appeared possible, based on the direct results of simulation and equation terms analysis, to resolve the problem in full. The author presents tabular conclusion of the results in Section 4.2.2. The table states precisely which imaginary forces terms, in what conditions, and for which vehicle major elements (bodies) may have an important influence on the simulation results. The last column states that imaginary force of translation, imaginary torque of rotation, and centrifugal force should not be neglected in both theoretical and practical analyses (calculation). The gyroscopic forces should not be neglected in the theoretical issues as well.

Both described results of the simulation studies are original and important contributions to the knowledge on rail vehicle dynamics. These contributions would not be possible without the mathematical, nominal, and numerical models of vehicles build based on the modeling methods discussed in the chapter. Therefore, it is hoped that the readers have obtained an idea of the efficiency and importance of both the modeling and the simulation for the development of contemporary rail vehicle dynamics, and thus the railways in general.

## **Author details**

Krzysztof Zboinski\*

Address all correspondence to: kzb@wt.pw.edu.pl

Faculty of Transport, Warsaw University of Technology, Poland

## **References**

In Sections 1 and 2, the general methods of modeling dynamics in moving coordinate systems, useful for multibody systems of rail vehicle type, were represented. The Newton–Euler and Kane's equations in the form suitable in AGEM approach were discussed. At the end of Section 2.2, the method profiled for rail vehicle systems was also presented. These methods help build mathematical models of rail vehicle dynamics. In Section 3, example nominal models of rail vehicles were introduced, as used by the author. In that section, the issue of building numerical

In Section 4, two examples of the author's simulation studies, with use of the numerical models, were discussed. The first example represented in Section 4.1 concerned the problem of rail vehicle stability in a curved track. Due to the author's systematic and consequent simulation studies, the method was formulated suitable for rail vehicle nonlinear lateral stability studies in a curved section of the track. The examples of the so-called stability map, being the most important result in the studies, were shown and briefly commented on. The stability maps approach was used by this author to study the influence of different factors on the results of

The second example is represented in Section 4.2. It was devoted to the use of the simulation in studying the influence of kinematics accuracy on vehicle dynamics in a curved track at variable velocities. In fact, both the constant and variable velocity cases were of interest. Due to the consequent author's interest in this issue, it appeared possible, based on the direct results of simulation and equation terms analysis, to resolve the problem in full. The author presents tabular conclusion of the results in Section 4.2.2. The table states precisely which imaginary forces terms, in what conditions, and for which vehicle major elements (bodies) may have an important influence on the simulation results. The last column states that imaginary force of translation, imaginary torque of rotation, and centrifugal force should not be neglected in both theoretical and practical analyses (calculation). The gyroscopic forces should not be neglected

Both described results of the simulation studies are original and important contributions to the knowledge on rail vehicle dynamics. These contributions would not be possible without the mathematical, nominal, and numerical models of vehicles build based on the modeling methods discussed in the chapter. Therefore, it is hoped that the readers have obtained an idea of the efficiency and importance of both the modeling and the simulation for the development

of contemporary rail vehicle dynamics, and thus the railways in general.

Address all correspondence to: kzb@wt.pw.edu.pl

Faculty of Transport, Warsaw University of Technology, Poland

models corresponding to the mathematical ones was also discussed.

192 Railway Research - Selected Topics on Development, Safety and Technology

stability analysis in a curved track.

in the theoretical issues as well.

**Author details**

Krzysztof Zboinski\*


and resulting activities. Oberpfaffenhofen: Society for Engineering and Scientific Ed‐ ucation; 1991.


ence on Railway Technology: Research, Development and Maintenance; 2014; Corsi‐ ca. Stirlingshire, UK: Civil-Comp Press; 2014. Paper 206. DOI: 10.4203/ccp.104.206

[31] Zboinski K, Golofit-Stawinska M. Simulation study into non-linear properties of rail‐ way vehicles in transition curve at velocities around the critical one. In: Materials of 23rd IAVSD Symposium on Vehicles on Roads and Tracks (electronic form); August 2013; Qingdao, China. 2013. Paper 2013-20.3-ID452.

and resulting activities. Oberpfaffenhofen: Society for Engineering and Scientific Ed‐

[16] Schiehlen W, editor. Multibody systems handbook. Berlin, Heidelberg, New York,

[17] Zboinski K. Dynamical investigation of railway vehicles on a curved track. European

[18] Dusza M, Zboinski K. Simulation study of 4-axle railway vehicle stability analysis in a curved track. In: Materials of 23rd IAVSD Symposium on Vehicles on Roads and Tracks (electronic form); August 2013; Qingdao, China. 2013. Paper 2013-20.3-ID429.

[19] Gear CW. Numerical initial value problems in ordinary differential equations. Engle‐

[20] Zboinski K. Numerical and traditional modelling of dynamics of multi-body system

[21] Zboinski K, Woznica P. Optimisation of the railway transition curves' shape with use of vehicle-track dynamical model. Archives of Transport. 2010;22(3):387–407.

[22] Kalker JJ. User's manual for program CONTACT vers. 1 and 3. Delft: Delft Universi‐

[23] Kalker JJ. A fast algorithm for the simplified theory of rolling contact. Vehicle System

[24] Iwnicki S, editor. Handbook of rail vehicle dynamics. Boca Raton, London, New

[25] Shabana AA, Zaazaa KE, Sugiyama H. Railroad vehicle dynamics: a computational approach. Boca Raton, London, New York: CRC Press, Taylor and Francis Group;

[26] Lee SY, Cheng YC. Hunting stability analysis of higher-speed railway vehicle trucks

[27] Kik W, Piotrowski J. A fast, approximate method to calculate normal load between wheel and rail and creep forces during rolling. In: Zobory I, editor. Proceedings of 2nd Mini-Conference Contact Mechanics and Wear of Rail/Wheel Systems ; Buda‐

[28] Kik W. Comparison of behaviour of different wheelset-track models. Vehicle System

[29] Kisilowski J, Knothe K, editors. Advanced railway vehicle system dynamics. War‐

[30] Zboinski K, Woznica P. Optimisation of railway transition curves with regard to wheel/rail wear. In: Pombo J, editor. Proceedings of the Second International Confer‐

on tangent track. Journal of Sound and Vibration. 2005;282(3–5):881–898.

in type of a railway vehicle. Archives of Transport. 2004;16(3):81–106.

London, Paris, Tokyo, Hong Kong: Springer-Verlag; 1990.

194 Railway Research - Selected Topics on Development, Safety and Technology

Journal of Mechanics, part A Solids. 1998;213(3):1001–1020.

wood Cliffs, New Jersey: Prentice -Hall Inc.; 1971. p. 253.

York: CRC Press, Taylor and Francis Group; 2006. 526 p.

ucation; 1991.

ty of Technology; 1986.

Dynamics. 1982;11:1–13.

pest. TU Budapest; 1996.

saw: WNT; 1991. 442 p.

Dynamics. 1992;20(suppl.):325–339.

2008. 331 p.


## *Edited by Krzysztof Zboinski*

This book focuses on selected research problems of contemporary railways. The first chapter is devoted to the prediction of railways development in the nearest future. The second chapter discusses safety and security problems in general, precisely from the system point of view. In the third chapter, both the general approach and a particular case study of a critical incident with regard to railway safety are presented. In the fourth chapter, the question of railway infrastructure studies is presented, which is devoted to track superstructure. In the fifth chapter, the modern system for the technical condition monitoring of railway tracks is discussed. The compact on-board sensing device is presented. The last chapter focuses on modeling railway vehicle dynamics using numerical simulation, where the dynamical models are exploited.

Railway Research - Selected Topics on Development, Safety and Technology

Railway Research

Selected Topics on Development,

Safety and Technology

*Edited by Krzysztof Zboinski*

Photo by zukanowa13 / iStock