Open AccessArticle

Vertical Oscillation of Railway Vehicle Chassis with Asymmetry Effect Consideration

Frantisek Klimenda

^1,*

Jan Skocilas

Blanka Skocilasova

¹,

Josef Soukup

¹ and

Roman Cizek

Faculty of Mechanical Engineering, University of Jan Evangelista Purkyne in Usti nad Labem, Pasteurova 1, 400 96 Usti nad Labem, Czech Republic

Faculty of Mechanical Engineering, Czech Technical University in Prague, 4 Technicka, 166 07 Prague, Czech Republic

Author to whom correspondence should be addressed.

Sensors 2022, 22(11), 4033; https://doi.org/10.3390/s22114033

Submission received: 20 April 2022 / Revised: 23 May 2022 / Accepted: 23 May 2022 / Published: 26 May 2022

(This article belongs to the Section Vehicular Sensing)

Download

Browse Figures

Review Reports Versions Notes

Abstract

This paper deals with the problem of vertical oscillation of rail and road vehicles under symmetrical and asymmetrical loading and symmetrical and asymmetrical kinematic excitation. The term asymmetry is understood as the asymmetric distribution of vehicle mass and elastic and dissipative elements with respect to the axes of geometric symmetry, including asymmetric kinematic excitation. The various models used (spatial, planar, quarter-plane) are discussed and their analytical solutions are outlined. The theory of the spatial model is applied to the chassis of a model railway vehicle. The basic relations for the calculation of the equations of motion of this vehicle are given. In the next section, the experimental solution of a four-axle platform rail car is described and the measurements of vertical displacement and accelerations when crossing wedges (representing unevenness) are given.

Keywords:

vehicle model; symmetry and asymmetry loading; kinematic excitation; motion equations; experimental results

1. Introduction

The problem of vibration in wheeled vehicles, both road and rail, is solved based on different assumptions and models, under different operating conditions. There is extensive literature on this issue, e.g., [1,2,3,4,5,6,7]. The solutions are represented by numerical methods [8] considering bogie track interaction with a system of nine degrees of freedom. The track unevenness excitation of vibration has been investigated by several authors [9] using a 68-degrees of freedom model, where results are compared with experiments. Various frequency ranges of dynamic system behaviour were investigated in [10]. One review summarises actual state of art of the results of railway vehicle vibration induced by dynamic impact loadings [11]. Furthermore, interaction between the wheel and the track has significant impact on the ride safety and control, resulting in the wearing of the wheel. A review of this problem is presented in another study [12]. The wheel flat defect caused by the above-mentioned phenomena and its impact on the wheel/rail dynamics can be described by a simple 3D model [13]. The excitation of the mechanical system (railway car) oscillation is usually modelled by passing the wheel over the uneven section of the road. Other types of irregularities such as stiffness irregularities, irregularities from different track positions, and irregularities in the wave propagation are investigated in [14] and also have an import role in drive comfort and safety. Other authors also paid attention to this issue from an experimental work perspective [15].

An overview of the solution to the vibrations of a spatially elastic body with consideration of various influences are given, for example, in [16]. The suppression of the influence of inadmissible vibrations and shocks in the investigation of the vibrations of suspended parts of rail and road vehicles is dealt with, for example, in [17,18,19,20]. The initial problem lies in investigating the vibrations of more complex systems of rigid bodies by elastic couplings. The most commonly used model for investigating the vibrations of a real structure is a system model with one or more degrees of freedom, solved as a one-dimensional chain of material points moving along a single line (or along straight parallel lines), possibly with the consideration of energy dissipation [21].

A major disadvantage of the current models, i.e., most often quarter-plane or planar, but also spatial, is the neglect of the effect of asymmetry on their driving characteristics. Therefore, only solutions to partial problems are given, mostly for half-plane (planar) models, which assume planar symmetry with one longitudinal axis of symmetry, respect the effect of the displacement of the centre of gravity from the centre of geometry, or consider different intensities of viscous damping and their combination. However, the most common solutions to the oscillation of vehicle models are based on quaternion models, usually vertical and elastic coupling with dissipative elements, and with different numbers of bodies—these are models with multiple degrees of freedom with coupled vertical displacements. In the application of the quaternary model, two axes of vehicle symmetry are assumed, i.e., full symmetry and its kinematic excitation.

It is known from practice that vehicle crashes, both road and rail, often occur for reasons that are difficult to explain. In the case of road vehicles, the cause is often cited as not adapting the speed to the road conditions, etc. In the case of railway vehicles, the causes are often more difficult to determine, apart from the obvious causes such as rail breakage, rail gauge collapse, etc. The cause is often due to incorrect load distribution and the resulting asymmetry of the load on the individual axles or wheels. In the case of vehicle oscillation, vertical vibrations can cause loss of wheel contact with the road (loss of vehicle controllability), and in the case of railway vehicles, loss of wheel contact with the rail due to unevenness of the rail, resulting in derailment.

Therefore, the investigation of the effect of asymmetry on the vertical oscillation of the model vehicle must be carried out by analysing the various causes and their consequences on spatial models.

When analysing the effect of asymmetry on the vertical oscillation of vehicles, we must distinguish three basic cases of the effect of asymmetry with respect to the axes of geometric symmetry, which are determined by two mutually perpendicular axes, the symmetry of the vehicle’s track and wheelbase, and the intersection at the geometric centre of the vehicle:

Asymmetry of the vehicle weight distribution with respect to the axes of geometric symmetry, the position of the centre of gravity, the directions of the main central axes of inertia, both of the vehicle structure itself and of the loaded vehicle
Asymmetry of the geometry of the distribution of elastic and dissipative elements of the couplings of the individual bodies of the vehicle system, their mechanical properties, spring stiffnesses, intensity of viscous damping assuming linear couplings of the individual quantities, and small displacements and rotations of the parts of the system
Asymmetry of kinematic excitation, i.e., the field of unevenness of the road or track surface, which define the kinematic excitation of the system at the point of contact between wheel-vehicle or wheel-track.

The above types of asymmetry can exist separately or together, and the third case is the most common.

The basic prerequisite for investigating the influence of asymmetry on the vibration of a system of bodies (a vehicle) is the choice of a suitable simple spatial model that allows for the investigation of different cases of asymmetry and the comparison or verification of the results obtained by different methods, i.e., theoretical (analytical, numerical, and simulation) and experimental. Most often, we compare the vertical displacement of the centre of gravity and the rotation about the axes passing through the centre of gravity and the vertical movements of any point of the vehicle.

In this study, a methodology for solving the vertical vibration of a vehicle was developed and this methodology was applied to several real vehicle cases (wheeled and tracked) to verify its applicability to a large set of vehicles. The methodology also allows for the determination of the magnitude and time history of wheel forces, the deformation magnitude or allowability of the individual elements of the spring links and parts of the structures, and the damping intensity criteria, including the fulfilment of the physiological requirements for the vehicle ride.

2. Theoretical Solution

We use the basic theory of motion equations of an elastically supported rigid body when solving the vertical vibration of the vehicles, both road and rail. A rigid body of general shape can be elastically supported or suspended in an inertial orthonormal Cartesian space, which is defined by a stationary component system. In our case, it is a more complex system where we solve the vertical oscillation of the vehicle (model) including the effect of asymmetry based on various assumptions. As already mentioned, a quarter, half, or full 3D model of the mechanical system with several degrees of freedom, with different mass distributions and with different elastic and dissipative elements (bindings) is used. The kinetic excitation is considered, and the solution will be demonstrated analytically, numerically, and experimentally.

For the basic analytical derivation, a simple but sufficiently general 3D model was designed to meet the requirements outlined in the previous section. The proposed model (Figure 1) is suitable for both the analytical derivation of relationships and experimental verification.

The model consists of a rigid plate, which is supported on four springs (k₁–k₄). The springs are supplemented by dampers (b₁–b₄). Weights z₁ and z₂ are placed on the plate to simulate the imbalance of the vehicle. The centre of gravity of the vehicle T is displaced in the x-axis direction by e_x and in the y-axis direction by e_y relative to the geometric centre C. The x and y axes passing through the geometric are the axle base and wheel spacing axes. The rotation of the model about the x-axis is φ_x and about the y-axis is φ_y. The value of h is the height of the road unevenness.

The time course of vertical displacement, i.e., the change in the position of the centre of gravity w_T(t) and the rotation of the model about the central axes of inertia φ_x(t), φ_y(t), was chosen as the criterion for comparing the individual cases.

In general, we can describe the motion of any point A (x_A, y_A) in terms of its vertical displacement, velocity, and acceleration versus time (e.g. at the suspension location) and determine the wheel pressures of the vehicle (Equations (1)–(3)).

w_{A} (t) {= w}_{T} (t) {+ x}_{A} φ_{y} (t) - y_{A} φ_{x} (t),

(1)

{\dot{w}}_{A} (t) = {\dot{w}}_{T} (t) + x_{A} {\dot{φ}}_{y} (t) - y_{A} {\dot{φ}}_{x} (t),

(2)

{\ddot{w}}_{A} (t) = {\ddot{w}}_{T} (t) {+ x}_{A} {\ddot{φ}}_{y} (t) - y_{A} {\ddot{φ}}_{x} (t) .

(3)

2.1. Individual Models Motion Equations

The solution was carried out for a spatial model (3D), a planar model, and a quartic model. The solution is found by determining the equations of motion of a mechanical system with three degrees of freedom.

2.1.1. Spatial Model

The equations of motion for the vibration of a spatial model of an elastically supported rigid body were determined from the Lagrange equations of the second kind

\frac{d}{d t} (\frac{\partial E_{k}}{\partial {\dot{q}}_{j}}) - \frac{\partial E_{k}}{\partial q_{j}} + \frac{\partial E_{p}}{\partial q_{j}} + \frac{\partial R_{d}}{\partial {\dot{q}}_{j}} = Q_{j} for j = 1, 2, \dots p,

(4)

where p—number of degrees of freedom, q_j—generalised coordinate, E_k—kinetic energy, E_p—potential energy, R_d—Rayleigh dissipated energy, Q_j—excitation forces.

The solution procedure is known from the literature. After substituting for the individual variables and performing the appropriate derivations, we obtain a system of linear inhomogeneous equations, which can be written in matrix form

M_{h} {\ddot{q}}_{j} + M_{b} {\dot{q}}_{j} + M_{k} q_{j} = Q_{j} (q_{j^{'}} {\dot{q}}_{j^{'}} t) for j = 1, 2, 3,

(5)

where M_h—mass matrix, M_b—damping matrix, M_k—stiffness matrix, Q_j—vector of generalised excitation function (in our case kinematic excitation of the type

Q = k_{j} q_{j} + b_{j} {\dot{q}}_{j}

To solve Equation (5), we use Laplace integral transforms, modifying the elements of the matrix M_h by multiplying it from the left by the diagonal matrix D = (d_ij) of the third order, whose elements are equal to the inverse of the diagonal elements of the matrix M_h.

In the case of a symmetrical mechanical system, the mass matrix is unitary. If we multiply Equation (5) by the diagonal matrix D from the left, we get

{DM}_{h} {\ddot{q}}_{j} + {DM}_{b} {\dot{q}}_{j} + {DM}_{k} q_{j} = {DQ}_{j} .

(6)

We rewrite this equation in the form

M {\ddot{q}}_{j} (t) + B {\dot{q}}_{j} (t) + {K q}_{j} ({t) = F}_{j} (t),

(7)

where M—mass matrix, B—damping matrix, K—stiffness matrix, F_j—vector of generalised kinematic excitation function.

The solution to Equation (7) is found by determining the elements of the damping matrix b_ij and the stiffness matrix a_ij. The calculation of these elements is given, e.g., in [22].

If the vertical component of the force at the m-th wheel location between the wheel and the ground is given by

Φ_{m} (t) {= k}_{m} h_{m} (t) {+ b}_{m} {\dot{h}}_{m} (t),

(8)

where k_m—spring stiffness of m-th wheel, h_m(t)—the height of the unevenness at the point of contact of m-th wheel, b_m—damping coefficient of m-th wheel, then the component of the generalised function F₁(t)—is physically the acceleration component of the vertical displacement [ms⁻²] and is given by the ratio of the product of the vertical components of the wheel forces and the mass of the vehicle model

F_{1} (t) {= m}^{- 1} \sum_{m = 1}^{4} Φ_{m} (t) .

(9)

The component of the generalised function F₂(t)—the physical component of the angular acceleration with respect to the x-axis [s⁻²]—is given by the ratio of the product of the moments of the components of the wheel forces to the x-axis and the moment of inertia of the vehicle model to this axis

F_{2} (t) {= J}_{x}^{- 1} [- \sum_{m = 1}^{2} Φ_{m} (t) l_{y m} + \sum_{m = 3}^{4} Φ_{m} (t) l_{y m} {+ e}_{y} \sum_{m = 1}^{4} Φ_{m} (t)] .

(10)

The component of the generalised function F₃(t)—physically the component of the angular acceleration with respect to the y-axis [s⁻²]—is given by the ratio of the product of the moments of the wheel force components to the y-axis and the moment of inertia of the vehicle model to this axis

F_{3} (t) {= J}_{y}^{- 1} [\sum_{m = 1, 4}^{} Φ_{m} (t) l_{x m} + \sum_{m = 2, 3}^{} Φ_{m} (t) l_{x m} {+ e}_{x} \sum_{m = 1}^{4} Φ_{m} (t)],

(11)

where l_xm a l_ym for m = 1, 2, 3, 4 are the distances of the support (elastic viscoelastic damping) from the axes

\bar{x}

and

\bar{y}

of the geometric symmetry with the origin at point C, J_x, J_y—moments of inertia of the vehicle mass to the respective axes l_xm, l_ym.

By defining the elements of the matrices M, B, K and the components of the vector F, a system of simultaneous linear inhomogeneous differential equations, the equations of motion of a spatial model of a vehicle with three degrees of freedom, with complete asymmetry: the mass distribution of the system (

e_{x} \neq {0, e}_{y} \neq {0, D}_{x y} \neq {0, D}_{y x} \neq 0

) of the geometry and stiffness of the elastic support, the geometry and intensity of the viscous damping, and with the asymmetry of the kinematic excitation defined by the road roughness h(x) → h(t) is determined.

It should be emphasised that the investigation of vehicle oscillations assuming a planar vehicle model with a longitudinal axis of symmetry, and especially assuming a quarter vehicle model with both longitudinal and transverse axes of symmetry, will substantially change the definition of the elements of the matrices of the equations of motion (7) and thus the required solution. This can be documented by comparing the matrix elements of the single-individual vehicle models.

2.1.2. Planar Model

For spatial mode of vehicle with a longitudinal axis of symmetry, which has only two degrees of freedom, w(t) and φ_y(t), the following can be assumed:

Mass matrix M is diagonal and unit, and for mass distribution is valid D_xy = D_yx = 0, => S₂₃ = S₃₂ = 0, e_y = 0.
The geometry of support is defined by dimensions l_y₁ = l_y₄, l_y₂ = l_y₃, l_x₁ = l_x₄, l_x₂ = l_x_3.
Damping intensity b₁ = b₄, b₂ = b₃.
Stiffness of elastic support is given by formulas k₁ = k₄, k₂ = k₃.

Then damping matrix B elements are

\begin{matrix} b_{11} {= 2 m}^{- 1} (b_{1} {+ b}_{2}), & b_{12} = 0, & b_{13} {= 2 m}^{- 1} (b_{1} l_{x 1} - b_{2} l_{x 2}) - e_{x} b_{11}, \\ b_{21} = 0, & b_{22} {= J}_{x}^{- 1} 2 \sum_{j = 1}^{3} b_{j} l_{y j}^{2} & b_{23} = 0, \\ b_{31} = m J_{x}^{- 1} b_{13}, & b_{32} = 0, & b_{33} {= J}_{y}^{- 1} [2 \sum_{j = 1}^{2} b_{j} l_{x j}^{2} {+ e}_{x}^{2} b_{11} m] . \end{matrix}

(12)

Stiffness matrix K elements are

\begin{matrix} a_{11} {= 2 m}^{- 1} (k_{1} {+ k}_{2}), & a_{12} = 0, & a_{13} {= 2 m}^{- 1} (k_{1} l_{x 1} - k_{2} l_{x 2}) - e_{x} a_{11}, \\ a_{21} = 0, & a_{22} {= 2 J}_{x}^{- 1} \sum_{j = 1}^{2} k_{j} l_{y j}^{2}, & a_{23} = 0, \\ a_{31} = m J_{x}^{- 1} a_{13}, & a_{32} = 0, & a_{33} {= J}_{y}^{- 1} [2 \sum_{j = 1}^{2} k_{j} l_{x j}^{2} {+ e}_{x}^{2} a_{11} m] . \end{matrix}

(13)

Based on the above relations, we can determine the components of the vertical forces acting at the location of the m-th wheel between the wheel and the road (rail) for crossing over unevenness of height h₁ = h₄, h₂ = h₃

\begin{matrix} F_{1} {= 2 m}^{- 1} [\sum_{j = 1}^{2} k_{j} h_{j} (t) + \sum_{j = 1}^{2} b_{j} {\dot{h}}_{j} (t)] {, F}_{2} = 0, \\ F_{3} (t) {= J}_{y}^{- 1} [\sum_{m = 1, 4}^{} Φ_{m} (t) l_{x m} - \sum_{m = 2, 3}^{} Φ_{m} (t) l_{xm} {+ e}_{x} \sum_{m = 1}^{4} Φ_{m} (t)] \end{matrix}

(14)

2.1.3. Quarter Model

For quarter model of vehicle with a longitudinal and transverse axis of symmetry, i.e., for a system of one degree of freedom w(t) the following is assumed:

Mass distribution D_xy = D_yx = 0 → S₂₃ = S₃₂ = 0, e_x = 0, e_y = 0.
Mass matric M is unit, the centre of gravity of the system is identical to the centre of geometric gravity C, the main central axes of inertia are identical to axes of geometrical symmetry.
The support geometry is determined by dimensions l_xj = l_x, l_yj = l_y for j = 1, 2, 3, 4.
Damping intensity _j = b, for j = 1, 2, 3, 4.
Stiffness of elastic support is given by k_j = k.

Then, damping matrix B elements are

\begin{matrix} b_{11} {= 4 m}^{- 1} b, & b_{12} = 0, & b_{13} = 0, \\ b_{21} = 0, & b_{22} {= 4 J}_{x}^{- 1} {k l}_{y}^{2}, & b_{23} = 0, \\ b_{31} = 0, & b_{32} = 0, & b_{33} {= 4 J}_{y}^{- 1} {k l}_{x}^{2} . \end{matrix}

(15)

Stiffness matrix K elements are

\begin{matrix} a_{11} {= 4 m}^{- 1} k, & a_{12} = 0, & a_{13} = 0, \\ a_{21} = 0, & a_{22} {= 4 J}_{x}^{- 1} {k l}_{y}^{2}, & a_{23} = 0, \\ a_{31} = 0, & a_{32} = 0, & a_{33} {= 4 J}_{y}^{- 1} {k l}_{x}^{2} . \end{matrix}

(16)

Components of the generalised kinematic excitation function vector for k_j = h

\begin{matrix} F_{1} {= 4 m}^{- 1} [k h (t) + b \dot{h} (t)], & F_{2} = 0, & F_{3} = 0 . \end{matrix}

(17)

2.1.4. Models Discussion

From the brief analysis of the equations of motion discussed above, and the forces on the individual wheels, it follows that both the plane and quarter models require complete symmetry along the longitudinal axis (plane model) or along both the longitudinal and transverse axes (quarter model).

Due to the distribution of the different structural groups in the vehicle, no vehicle is geometrically symmetrical, it is clear that neither the quarter nor the half model are able to determine the dynamic properties of the vehicle with sufficient accuracy. In addition, the asymmetry in the loading of the vehicle (both cargo and possibly passengers) must be taken into account. For this reason, it is advisable to use a full spatial model for the further solution, which may be generally unsymmetrical (inertia and geometric axes are not parallel) or only unsymmetrical but with parallel inertia and geometric symmetry axes.

2.2. Simple Model of Four Axles Wagon

We will show the solution of this system on the chassis of a model of a railway vehicle, considering the effect of asymmetry. We assume a forced oscillation of the model, which is represented by a spatially elastically supported rigid plate with single and multiple primary linear suspension by coil springs (two-axle railway vehicle model). It is a kinematically excited system of three rigid bodies elastically supported and bounded, considering the effect of asymmetry [1,23,24]. However, these authors did not address the issue of asymmetry.

A simple computational model of the vehicle chassis was chosen for the calculation. The model consists of two two-axle chassis with simple primary suspension and simple vehicle body suspension. The chassis are replaced by rigid plates of mass m₁ and m₂ with symmetric mass distribution (centre of gravity is identical to the kinematic centre, and the main central axes of inertia are identical to the axes of geometric symmetry). However, we consider the asymmetry of the spring stiffness parameter and their geometrical support. The model of the vehicle body is considered as a rigid plate of mass m for which an asymmetric mass distribution is considered (the position of the centre of gravity is deflected by e_x and e_y from the geometric centre, and the main central axes of inertia are rotated with respect to the axes of geometric symmetry). The computational model is shown in Figure 2.

The derivation of the following relationships is based on general assumptions, i.e., body stiffness, small displacements and rotations, and linear spring characteristics. We consider only the vertical displacements of arbitrary points of individual bodies. This vertical change in the position of a point of a body is determined by the displacement of the body’s centre of gravity w, w₁ and w₂, the rotations φ_x, φ_x1, φ_x2, φ_y, φ_y1, φ_y2 about the central axes of inertia of individual masses m, m₁, m₂ and its distance from these axes. Thus, we solve a system of bodies with nine degrees of freedom.

The equations of motion are derived from the Lagrange equations of the second kind (4). Kinetic energy of the system

\begin{matrix} E_{k} = \frac{1}{2} m {\dot{w}}^{2} + \frac{1}{2} (J_{x} {\dot{φ}}_{x}^{2} + J_{y} {\dot{φ}}_{y}^{2} - 2 D_{x y} {\dot{φ}}_{x} {\dot{φ}}_{y}) + \frac{1}{2} m_{1} {\dot{w}}_{1}^{2} + \frac{1}{2} J_{x 1} {\dot{φ}}_{x 1}^{2} + \frac{1}{2} J_{y 1} {\dot{φ}}_{y 1}^{2} + \\ + \frac{1}{2} m_{2} {\dot{w}}_{2}^{2} + \frac{1}{2} J_{x 2} {\dot{φ}}_{x 2}^{2} + \frac{1}{2} J_{y 2} {\dot{φ}}_{y 2}^{2}, \end{matrix}

(18)

where J_x, J_y—moments of inertia, D_xy—deviation moment to central axes of inertia of chassis with mass m, J_x₁, J_y₁, J_x₂ and J_y₂—inertia moments to the main central axes of bogies with mass m₁, m₂.

The potential energy of the system depends on the displacements of the individual masses in the places of their elastic support, i.e., in the places of their elastic bound. In our case, in the chassis model we have chosen the system of marking for points A_jki, their coordinates x_jki, y_jki, vertical displacements w_jki and stiffness constants k_jki. The individual indices correspond to masses j = 1, 2, quadrants k = 1, 2, 3, 4 and spring order i = 1, 2, ...., n (for the case of multiple support).

Similarly, the vertical change h_jki is marked, which is at the point of contact of the axle (spring) with the rigid base, to which we relate the position of the body. This vertical change represents crossing over an unevenness in the road (rail)—kinematic excitation.

We can therefore determine the displacements of the individual points of the chassis, which are given by the relations for body m₁, j = 1

Point	Vertical displacement	Stiffness constant
A₁₁₁	$w_{111} {= w}_{1} - y_{111} φ_{x 1} {+ x}_{111} φ_{y 1} - h_{111,}$	$k_{111}$	(19)
A₁₂₁	$w_{121} {= w}_{1} - y_{121} φ_{x 1} - x_{121} φ_{y 1} - h_{121},$	$k_{121}$
A₁₃₁	$w_{131} {= w}_{1} {+ y}_{131} φ_{x 1} - x_{131} φ_{y 1} - h_{131,}$	$k_{131}$
A₁₄₁	$w_{141} {= w}_{1} {+ y}_{141} φ_{x 1} {+ x}_{141} φ_{y 1} - h_{141 .}$	$k_{141}$

For body m₂, j = 2

Point	Vertical displacement	Stiffness constant
A₂₁₁	$w_{211} {= w}_{2} {+ y}_{211} φ_{x 2} {+ x}_{211} φ_{y 2} - h_{211,}$	$k_{211}$	(20)
A₂₂₁	$w_{221} {= w}_{2} - y_{221} φ_{x 2} - x_{221} φ_{y 2} - h_{221},$	$k_{221}$
A₂₃₁	$w_{231} {= w}_{2} {+ y}_{231} φ_{x 2} - x_{231} φ_{y 2} - h_{231,}$	$k_{231}$
A₂₄₁	$w_{241} {= w}_{2} {+ y}_{241} φ_{x 2} {+ x}_{241} φ_{y 2} - h_{241 .}$	$k_{241}$

In the case of a model of a body with a simple spring suspension in the axis of the chassis, the division into quadrants cannot be used to mark the individual points of action of the springs. Marking of the points, coordinates, and stiffness constants B_jki, x_jki, y_jki, k_jki, j = 0 body, mass m, k = 1, 2 marking of chassis 1 and 2, respectively, i = 1, 4 belonging to the half-point A₁₁₁, respectively, A₁₄₁ follows from Figure 2.

Displacement of individual point of body are given by

Point	Vertical displacement	Stiffness constant
B₀₁₁	$w_{011} {= w + y}_{011} φ_{x} {+ x}_{011} φ_{y} - w_{1} + (y_{011} {+ e}_{y}) φ_{x 1,}$	$k_{011}$	(21)
B₀₁₄	$w_{014} {= w + y}_{014} φ_{x} {+ x}_{014} φ_{y} - w_{1} - (y_{014} {+ e}_{y}) φ_{x 1,}$	$k_{014}$
B₀₂₁	$w_{021} {= w - y}_{021} φ_{x} {+ x}_{021} φ_{y} - w_{2} + (y_{021} - e_{y}) φ_{x 2,}$	$k_{021}$
B₀₂₄	$w_{024} {= w + y}_{024} φ_{x} {+ x}_{024} φ_{y} - w_{2} - (y_{024} - e_{y}) φ_{x 2 .}$	$k_{024}$

Potential energy of bogies m₁, m₂ and body m is

E_{p} = \frac{1}{2} \sum_{j = 1}^{2} \sum_{k = 1}^{4} \sum_{i = 1}^{k_{j}} k_{j k i} w_{j k i}^{2} + \frac{1}{2} \sum_{j = 0}^{} \sum_{k = 1, 2}^{} \sum_{i = 1, 4}^{} k_{j k i} w_{j k i}^{2} .

(22)

Kinetic energy from (18) and potential energy from (22) are substituted into the Equation (4) resulting in

\frac{d}{d t} \frac{\partial E_{k}}{\partial {\dot{q}}_{j}} - \frac{\partial E_{k}}{\partial q_{j}} + \frac{\partial E_{p}}{\partial q_{j}} = Q_{j} j = 1, \dots, 9,

(23)

where j—degree of system freedom, q_j = (φ_x, φ_y, w, w₁, w₂, φ_x₁, φ_y₁, φ_x₂, φ_y₂)^T and Q_j—generalised forces.

The solution (after derivation with respect to individual coordinates) is given by motion equations in matrix form

\begin{matrix} [\begin{matrix} J_{x} & - D_{x y} & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ - D_{x y} & J_{y} & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & m & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & m_{1} & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & m_{2} & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & J_{x 1} & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & J_{y 1} & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & J_{x 2} & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & J_{y 2} \end{matrix}] \cdot [\begin{matrix} {\ddot{φ}}_{x} \\ {\ddot{φ}}_{y} \\ \ddot{w} \\ {\ddot{w}}_{1} \\ {\ddot{w}}_{2} \\ {\ddot{φ}}_{x 2} \\ {\ddot{φ}}_{y 1} \\ {\ddot{φ}}_{x 2} \\ {\ddot{φ}}_{y 2} \end{matrix}] + \\ + [\begin{array}{ccccccccc} κ_{11} & κ_{12} & κ_{13} & κ_{14} & κ_{15} & κ_{16} & 0 & κ_{18} & 0 \\ κ_{21} & κ_{22} & κ_{23} & κ_{24} & κ_{25} & κ_{26} & 0 & κ_{28} & 0 \\ κ_{31} & κ_{32} & κ_{33} & κ_{34} & κ_{35} & κ_{36} & 0 & κ_{38} & 0 \\ κ_{41} & κ_{42} & κ_{43} & κ_{44} & 0 & κ_{46} & κ_{47} & 0 & 0 \\ κ_{51} & κ_{52} & κ_{53} & 0 & κ_{55} & 0 & 0 & κ_{58} & κ_{59} \\ κ_{61} & κ_{62} & κ_{63} & κ_{64} & 0 & κ_{66} & κ_{67} & 0 & 0 \\ 0 & 0 & 0 & κ_{74} & 0 & κ_{76} & κ_{77} & 0 & 0 \\ κ_{81} & κ_{82} & κ_{83} & 0 & κ_{85} & 0 & 0 & κ_{88} & κ_{89} \\ 0 & 0 & 0 & 0 & κ_{95} & 0 & 0 & κ_{98} & κ_{99} \end{array}] \cdot [\begin{matrix} φ_{x} \\ φ_{y} \\ w \\ w_{1} \\ w_{2} \\ φ_{x 2} \\ φ_{y 1} \\ φ_{x 2} \\ φ_{y 2} \end{matrix}] = [\begin{matrix} 0 \\ 0 \\ 0 \\ Q_{4} \\ Q_{5} \\ Q_{6} \\ Q_{7} \\ Q_{8} \\ Q_{9} \end{matrix}] \cdot \end{matrix}

(24)

The individual elements of the stiffness matrix κ_ij are functions of stiffness constants k_jki, the dimensions of support x_jki, y_jki and the eccentricities e_x and e_y. The detailed solution of the individual elements is given, e.g., in [22].

The generalised force functions Q_j ((23) and (24)) are again functions of the stiffness constants of the chassis springs k_jki and functions of the road unevenness (elevation) at time h(t) at the location of the wheel springs. Then the kinematic excitation functions in Equation (24) are given by

\begin{matrix} Q_{1} {= Q}_{2} {= Q}_{3} = 0, \\ Q_{4} {= k}_{111} h_{111} {+ k}_{121} h_{121} {+ k}_{131} h_{131} {+ k}_{141} h_{141}, \\ Q_{5} {= k}_{211} h_{211} {+ k}_{221} h_{221} {+ k}_{231} h_{231} {+ k}_{241} h_{241}, \\ Q_{6} = - k_{111} h_{111} y_{111} - k_{121} h_{121} y_{121} {+ k}_{131} h_{131} y_{131} {+ k}_{141} h_{141} y_{141}, \\ Q_{7} {= k}_{111} h_{111} x_{111} - k_{121} h_{121} x_{121} - k_{131} h_{131} x_{131} {+ k}_{141} h_{141} x_{141}, \\ Q_{8} = - k_{211} h_{211} y_{211} {+ k}_{221} h_{221} y_{221} - k_{231} h_{231} y_{231} {+ k}_{241} h_{241} y_{241}, \\ Q_{9} {= k}_{211} h_{211} x_{211} - k_{221} h_{221} x_{221} - k_{231} h_{231} x_{231} {+ k}_{241} h_{241} x_{241} \end{matrix}

(25)

where k_jkih_jki(t)—the shape of the time function of kinematic excitation at individual axle points and their sequence.

Equation (24) can be rearranged to

M_{h} {\ddot{q}}_{j} {+ M}_{k} q_{j} {= Q}_{j} (t) .

(26)

Then the equation is multiplied by matrix D (d_ij) from the left, where d_ij = 1/α_ij—diagonal elements of mass matrix M_h(α_ij). After arrangement, the result is [25,26]

M {\ddot{q}}_{j} {+ K q}_{j} {= F}_{j} (t) .

(27)

Mass matrix can be expressed in form M = E + S and is substituted into (26)

(E + S) {\ddot{q}}_{j} {+ K q}_{j} {= Q}_{j} (t),

(28)

where elements of matrix S are s_ij. These elements are given by equations s₁₁ = s₂₂ = 0; s_ij = 0 for i = 3, …, 9; j = 3, …, 9; s₁₂ =

-

D_xy/J_x; s₂₁ =

-

D_xy/J_y. These quantities represent the effect of asymmetric distribution of mass, i.e., rotation of main central inertia axes to central axes parallel with geometric axes of the plate (vehicle). Similarly, the elements a_ij of matrix K are determined by division of i-th row of stiffness matrix by i-th element at diagonal of inertia matrix. Excitation function of time F_j(t) is determined by the same procedure with function Q_j(t).

System of motion equations in case of symmetry is transformed into

E {\ddot{q}}_{j} {+ K q}_{j} {= F}_{j} for i = 1 \dots 9 .

(29)

where stiffness matrix is

K = [\begin{matrix} a_{11} & 0 & 0 & 0 & 0 & a_{16} & 0 & a_{18} & 0 \\ 0 & a_{22} & 0 & a_{24} & a_{25} & 0 & 0 & 0 & 0 \\ 0 & 0 & a_{33} & a_{34} & a_{35} & 0 & 0 & 0 & 0 \\ 0 & a_{42} & a_{43} & a_{44} & 0 & 0 & 0 & 0 & 0 \\ 0 & a_{52} & a_{53} & 0 & a_{55} & 0 & 0 & 0 & 0 \\ a_{61} & 0 & 0 & 0 & 0 & a_{66} & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & a_{77} & 0 & 0 \\ a_{81} & 0 & 0 & 0 & 0 & 0 & 0 & a_{88} & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & a_{99} \end{matrix}] .

(30)

Elements of stiffness matrix are function of spring stiffness k₁₂ and k₀ and function of displacement x₁₂, y₁₂, x₀, y₀. In the case of complete symmetry, the system of nine differentiated simultaneous equations is transformed to the system of seven simultaneous equations and two independent equations of harmonic motion for q₇ and q₉(y₉).

{\ddot{q}}_{7} {+ a}_{77} q_{7} {= F}_{7} (t) {\ddot{y}}_{9} {+ a}_{99} y_{9} {= F}_{9} (t) .

(31)

The solution of the system of Equation (26) can be performed by any of the numerical methods (MATLAB, MAPLE, etc.). Analytical solution is possible by applying matrix calculus, Lagrange’s method of variance of constants, or by Laplace transform. The Laplace transformation yields a system of linear algebraic equations for zero initial conditions [27].

As indicated, the solution can also be done by applying the Laplace transformation, where Equation (27) after modifications becomes

[(E + S) p^{2} + K] {\bar{y}}_{i} (p) = {\bar{F}}_{i} (p),

(32)

where

{\bar{y}}_{i} (p)

and

{\bar{F}}_{i} (p)

are images of function q_j(t) and F_j(t) for j = 4 … 9, p—parameter of Laplace transformation. By rewriting the Equation (32) into the matrix form, it is recommended to solve the system of linear algebraic equations by Cramer rule

{\bar{y}}_{i} (p) = \frac{D_{j} (p)}{D (p)} = \sum_{I = 4}^{n = 9} {(- 1)}^{j + i} {\bar{F}}_{i} (p) \frac{D_{j i} (p)}{D (p)},

(33)

where D(p)—determinant of matrix of above-mentioned system of equations, which is defined as

{D (p) = C}_{A} {I for n = 9 C}_{A} = 1 - s_{12} s_{21} .

(34)

Equations for evaluation of real coefficients A_{2(n − i)} can be found in [12]. These equations are generally valid for n ≥ 2 and 0 ≤ i < n.

Equations for evaluation of coefficients A_{2(n − i)} for 2 < i < n are very complex and should be determined using MAPLE or MATLAB.

To determine the search required function q_j(t) by reverse transformation of the images

{\bar{y}}_{i} (p)

, it is convenient to modify the ratio of the determinants D_ji(p) and D(p) so that the product of the expression and image

{\bar{F}}_{i} (p)

obtained by the modification allows the application of the image convolution theorem.

The polynomial in Equation (34) with real coefficients can be replaced by the product of quadratic binomials

\sum_{i = 0}^{n} A_{2 (n - 1)} p^{2 (n - 1)} = \prod_{i = 1}^{n} (p^{2} {+ ω}_{i}^{2}) .

(35)

After performing the product on the right side, a polynomial is obtained

\sum_{i = 0}^{n} A_{2 (n - 1)} p^{2 (n - 1)} = \sum_{i = 0}^{n} B_{2 (n - i)} p^{2 (n - i)},

(36)

where real coefficients

B_{2 (n - i)}

are given by

\begin{matrix} B_{2 n} = 1, & B_{2 (n - 1)} = \sum_{i = 1}^{n} ω_{i}^{2}, \\ B_{2 (n - 2)} = \sum_{i = 1}^{n - 1} \sum_{j = i + 1}^{n} ω_{i}^{2} ω_{j}^{2}, & B_{2 (n - 3)} = \sum_{i = 1}^{n - 2} \sum_{j = i + 1}^{n - 1} \sum_{l = j + 1}^{n} ω_{i}^{2} ω_{j}^{2} ω_{l}^{2}, \\ B_{2 (n - 4)} = \sum_{i = 1}^{n - 3} \sum_{j = l + 1}^{n - 2} \sum_{j = 1}^{n - 1} \sum_{k = l + 1}^{n} ω_{i}^{2} ω_{j}^{2} ω_{l}^{2} ω_{k}^{2}, & B_{2 (n - n)} {= B}_{0} \prod_{i = 1}^{n} ω_{i}^{2} . \end{matrix}

(37)

By comparing the coefficients of the polynomials for the same powers of the parameter p^{2(n − i)} on both sides of Equation (36) A_{2(n − i)} = B_{2(n − i)} a system of equations is obtained to determine ω_i, ω_j, ω_l, ….

By comparing Equation (37) with the relations determining the coupling of the coefficients of the algebraic equation and its root factors, the algebraic equation is obtained

\begin{matrix} f (ω^{2}) {= ω}^{2 n} - A_{2 (n - 1)} ω^{2 (n - 1)} {+ A}_{2 (n - 2)} ω^{2 (n - 2)} - A_{2 (n - 3)} ω^{2 (n - 3)} + \\ + {(- 1)}^{n - 1} A_{2} ω^{2} + {(- 1)}^{n} A_{0} = 0 . \end{matrix}

(38)

Determining the roots of the frequency Equation (38) is the most complicated part of the proposed solution, especially in terms of numerical accuracy, i.e., ω²—the circular frequencies of the functions y_j(t), for j = 1, 2, 3, …, n. The next solution, after modifications and back-transformation, would give the convolution integral

q_{j} (t) = \sum_{i = 4}^{n} {(- 1)}^{j + i} \sum_{k = 1}^{n} \frac{L_{j i, k}}{ω_{k}} \int_{0}^{t} F_{i} (τ) {\sin ω}_{k} (t - τ) d τ

(39)

where ω_k is the solution of Equation (38). From the known functions q_j(t), i.e., the solution of the system of Equation (26), the required quantities of the system of Equation (24) are determined.

This method makes it possible to determine the vertical displacements of arbitrary points of the chassis or body of the vehicle, for example, the displacements of the points described by Equations (19)–(21). Thus, a program can be developed to calculate the vertical displacements of arbitrary points of railway vehicle chassis models with a large range of design variations.

3. Experimental Methods

The above method was experimentally verified on a real railway vehicle. A four-axle railway freight car of the Smmps series was used for the experiment. The aim of the measurements was to determine the geometrical, mass, and stiffness parameters of the single-individual parts of the freight car. The main objective was to establish the oscillations of the individual parts of the wagon when crossing the wedges at different wedge locations and load placements on the wagon platform (wedges = kinematic excitation, load placement = symmetrical or asymmetrical load). The measurement requirements were specified, and the wagon was modified for these specific requirements.

3.1. Description of Vehicle Arrangement

Due to the requirements for suspension linearity, the original chassis with leaf springs (type 26-2.8) were replaced by chassis with coil springs (Y25). For the new chassis, the outer springs of the primary suspension were removed (Figure 3a), and a secondary suspension was created by placing three springs on the longitudinal beam of each chassis (Figure 3b).

The bogies were fixed to the body (car platform) by a network of steel cables with cable tensioners—Figure 4a.

For all springs (internal and external) of both chassis, their characteristics and stiffnesses were determined, and based on these characteristics the selection was made so that the resulting stiffnesses of the primary and secondary suspension on the wheel, wheelset, chassis, and car were distributed with maximum symmetry to the longitudinal and transverse axes of the car.

The flatcar was loaded with two loading units (Figure 4b), each weighing 4000 kg (the units were made up of concrete blocks with dimensions of 4900 × 1000 × 300 mm). For the chassis, the masses of the parts, their centre of gravity coordinates, and the moments of inertia to the axes passing through the centre of gravity were experimentally determined. Similar parameters were determined for the vehicle body. For these measurements, a proprietary methodology was developed.

3.2. Methodology

The modified vehicle was subjected to wedge crossing tests (Figure 5). The car was loaded with loading units (either symmetrically or unsymmetrically). The location of the wedges was varied during the tests depending on the kinematic excitation requirement (symmetrical, asymmetrical, single axle, both axles, single chassis, both chassis).

The relative vertical displacement of the chassis frame relative to the wheelset (9 sensors), the car body relative to the chassis frame (4 sensors), and the vertical acceleration of selected points on the body (5 sensors) and chassis (2 sensors) were measured. The location of the sensors can be seen in Figure 6. During the tests, time histories of the measured quantities were recorded and used to verify the theoretical model.

Strain gauge track sensors (chassis frame × wheelset) and cable track sensors (box frame × chassis frame) were used to measure vertical displacement. Acceleration was measured with B12/200 Hz inductive sensors (Figure 7).

The experiments were carried out according to an optimised test plan. Selected combinations of experiments arose from statistical analysis of combinations of all possible experimental settings. These combinations were sufficient to achieve the goals, i.e., to identify critical situations in which derailment may occur [28]. The tests were carried out for different cases of load distribution on the car platform and wedge placement on the rails, taking into account the wheelbase of the bogies (1800 mm) and the car wheelbase (9400 mm). The test plan can be seen in Figure 8, where the different variants of the car crossing the wedges are shown. A total of 30 combinations were measured. Each measurement was repeated 3 times and a total of 90 tests were performed. Prior to the actual car crossings over the wedges, the front and rear heights of the bogie frames and car body above the tops of both rails were measured. The static seating of the primary and secondary suspension springs under a given load (symmetrical, unsymmetrical) was checked. The car was pulled over the wedges by a motorised locomotion system using a steel cable. The wedges were attached to the rail sections that were part of the block for dynamic measurements.

3.3. Experimental Results

Due to the scope of the measurements, only selected results for symmetrical and asymmetrical loading and for symmetrical and asymmetrical excitation are presented in this paper. The labels of the locations where the relative deflections between the case and chassis, chassis, and wheelset were measured are shown in Figure 9.

(a): The time trend of the deflections when the car crosses the wedges under symmetrical loading of the car and symmetrical kinematic excitation is shown in Figure 10 and Figure 11.

In this case, the symmetrical load is represented by the symmetrical positioning of the load with respect to the geometric centre of the car. Figure 10 shows the time history of the vertical displacement of the first wheelset relative to the chassis frame. The waveform is symmetric about zero displacement. The only damping that affects the deflection is material damping. The variation in the vertical displacement is due to slight differences in the material properties of the springs of the first wheelset. The displacements show an excellent agreement.

Figure 11 shows the vertical displacement of the chassis frame at two locations, see Figure 9 (2L, 2R), relative to the car body. Both places are located off the geometric centre of the chassis, therefore, in agreement with Equation (1), the waveforms are displaced off the axis of symmetry. Again, the waveforms show a very good agreement, and it can be seen that the chassis frame is rotated relative to the longitudinal and transverse axes when crossing the wedges.

(b): Time trends of displacement when the car crosses the wedges under symmetrical loading, asymmetrical kinematic excitation—Figure 12 and Figure 13.

Figure 12 shows the vertical displacement between the first wheelset and the chassis frame under symmetrical loading and asymmetrical kinematic excitation (variant A–IV). In this kinematic excitation, there is a different excitation of the first wheelset on the left and right side, which can be seen in the deflection of waveforms. In this kinematic excitation, the wheelsets rotate about the longitudinal axis of the car. The whole chassis and body of the car are also rotated about the longitudinal and transverse axes. The figure also shows that at certain stages the chassis wheel loses contact with the rail (loss of adhesion).

From the above time trend (Figure 13), it is clear that due to the oscillation and the resulting rotation of both the wheelset and the chassis frame at the measuring points, there is also a rotation of the body. This kinematic excitation causes the wheel to bounce significantly off the rail and, under certain conditions (e.g., in a horizontal curve, on a turnout), could cause a derailment due to the large displacement.

(c): Time trends of displacement when the car crosses the wedges under asymmetrical loading, symmetrical kinematic excitation—Figure 14 and Figure 15.

The displacement of the front wheelset relative to the chassis frame (Figure 14) show a difference in position relative to the zero deflection axis even under symmetrical kinematic excitation and asymmetrical loading. This is due to the non-symmetrical load distribution and hence the rotation of the car body and the different loading on the left and right side of the chassis. The maximum displacement is less than that of variant A-IV.

The relative deflection between the chassis frame and the car body is again shifted relative to the axis through the centre of gravity in the case D-I, i.e., under asymmetrical loading and symmetrical kinematic excitation (Figure 15). The displacements are much larger at their maximum values than between the chassis frame and the wheelset. The difference in their trends is due to the fact that in one case (2R) the displacement is in the axis of transverse symmetry of the chassis frame. In 2L case the displacement is affected by rotation and displacement from the transverse axis—see Equation (1). The maximum deflections are again quite significant and reach a maximum of about 18 mm.

(d): Time trends of displacement when the car crosses the wedges under asymmetrical loading, asymmetrical kinematic excitation—Figure 16 and Figure 17.

The load distribution is the same as in the previous case, but the kinematic excitation is asymmetrical (Figure 16). Displacement trend between the chassis frame and the wheelset reaches more than 20 mm at the maximum (this is less than the wedge height). The decay to an acceptable value occurs relatively quickly, in about 20 s. It is also clearly visible in the figure that the first and second wheelsets are repeatedly dropped on each side (the wedges are directly connected to each other).

The difference in the position of the vertical displacement between the car body and the chassis frame (Figure 17) at points 2L and 2R is due to the different position of these points relative to the transverse axis of the chassis frame and the transverse axis of the car body. The magnitude of the displacement is influenced by the asymmetrical loading (variant D-IV) which caused the car body to rotate relative to its longitudinal and transverse axes. The consequence is also a different load on the left and right side of all wheelsets.

In addition to the measurements of the relative vertical displacement, measurements of the vertical acceleration were also carried out. Most of the acceleration measurement points were located on the vehicle body (five sensors in total). In this paper we present the vertical acceleration waveform of the vehicle body in its transverse axis on the left and right side (points 3L and 3R) for asymmetrical loading and asymmetrical kinematic excitation (variant B-III). To illustrate this situation, the time history of the relative displacements at points 1L, 1R (Figure 18), 2L, 2R (Figure 19) is also shown for this case.

Figure 18 shows the relative displacement between the chassis frame and the wheelset (B-III). This variant is characterised by an asymmetrical load (load placed on the right side of the car body symmetrical to the transverse axis). The kinematic excitation is also asymmetrical. This load has caused the car body to rotate about the longitudinal axis and thus also the load on the chassis to be unequal. From the displacement waveforms it is evident that the larger amplitude after dropping from the wedges are on the unloaded side. After the last wedge jump on the loaded side the displacement is maximum. The displacement on the unloaded side reached almost the same values immediately after the wedge crossings. A gradual damping of the deflections occurs after 27 s. On the unloaded side, there is a significant loss of wheel-rail contact after each wedge, while on the loaded side there are only two significant peaks, after the first and last wedges.

The relative displacement between the car body and the chassis frame (Figure 19) are relatively far apart (about 35 mm) due to the rotation of the car body and therefore the different loads on the left and right sides. The time histories of the deflections are almost identical, but it is clear that the left side was already lightened before the test. The maximum deflections of the loaded side are higher than those of the unloaded side when jumping off the wedges, which is consistent with the nature of the oscillation of the car body. Such a load distribution under asymmetrical kinematic excitation may also cause the car to derail.

The vertical acceleration in the middle of the car body on both the left and right sides (Figure 20) was higher on the unloaded side (3L) than on the loaded side (3R) at the peak of the wedge jump. The maximum value was reached after crossing the first wedge, and the maxima decreased when crossing the other wedges.

4. Conclusions

The paper deals with the theory of oscillations of mechanical systems of elastically supported and bound bodies. The application of this solution to vehicle models has been carried out. The methodology of the model solution is proposed and the application of this methodology to a real vehicle under different operating conditions has been carried out. The application has been made, among other things, to a railway chassis vehicle—its solution is the subject of this paper. Based on the developed methodology, an experimental solution of a four-axle railway vehicle without damping (shock absorbers) was implemented. The results of the experimental solution, i.e., the time course of vertical displacement, are presented in the previous section.

From the given results, it can be seen that in the case of asymmetrical kinematic excitation (e.g., rail fracture, unevenness on the rail), the vertical oscillation of the car will cause the vertical oscillation. The consequence of this oscillation is the loss of contact between the wheel and the rail (loss of adhesion). In any case where the critical value of vertical displacement of the wheel, which is equal or higher than height of axle, is achieved, the vehicle may be derailed. These vertical deflections are also strongly influenced by the asymmetry of the weight distribution of the vehicle and, in particular, by the asymmetry of the load. From the results obtained, it is necessary to pay particular attention to the distribution of the load (especially piece loads) on the railway carriage.

Author Contributions

Conceptualization, J.S. (Josef Soukup), J.S. (Jan Skocilas), B.S., F.K. and R.C.; methodology, J.S. (Josef Soukup) and B.S.; software, J.S. (Jan Skocilas); validation, F.K., R.C. and J.S. (Josef Soukup); formal analysis, F.K. and J.S. (Josef Soukup); investigation, J.S. (Josef Soukup), J.S. (Jan Skocilas), B.S., F.K. and R.C.; data curation, F.K. and J.S. (Josef Soukup); writing—original draft preparation, F.K. and J.S. (Josef Soukup); writing—review and editing, F.K.; visualization, F.K. and J.S. (Josef Soukup). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of the Ethics Committee of Faculty of Mechanical Engineering Jan Evangelista Purkyne University in Usti nad Labem (20 April 2022).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

Mitschke, M. Dynamik der Kraftfahrzeng, 4th ed.; Springer: Berlin/Heidelberg, Germany, 2004; p. 806. ISBN 978-3540420118. [Google Scholar]
Rus, L. Dynamic of Railway Vehicles. Ph.D. Thesis, VÚML ČKD, Praha, Czech Republic, 1988. (In Czech). [Google Scholar]
Rajalingham, C.; Rakheja, S. Influence of Suspension Damper Asymmetry on Vehicle Vibration Response to Ground Excitation. J. Sound Vib. 2003, 266, 1117–1129. [Google Scholar] [CrossRef]
Karnopp, D. Vehicle Stability; Marcel Dekker: New York, NY, USA, 2004; ISBN 0-8247-5711-4. [Google Scholar]
Biggs, J.M.; Testa, B. Introduction to Structural Dynamics; McGraw-Hill: New York, NY, USA, 1964. [Google Scholar]
Knothe, K.; Stichel, M. Rail Vehicle Dynamics; Springer International Publishing: Cham, Switzerland, 2017; ISBN 978-3-319-45374-3. [Google Scholar] [CrossRef]
Iwnicki, S.; Spiryagin, M.; Cole, C.; Mcsweeney, T. (Eds.) Handbook of Railway Vehicle Dynamics; CRC Press: Boca Raton, FL, USA, 2019; ISBN 9780429469398. [Google Scholar] [CrossRef]
Dumitriu, M.; Gheti, M.A.; Bondera, I.; Cofaru, N.F.; Inta, M. On the longitudinal vibration of the railway bogie. In Proceedings of the MATEC Web of Conferences, Online, 21 August 2019; Volume 290. [Google Scholar] [CrossRef] [Green Version]
Chelli, F.; Corradi, R. On rail vehicle vibrations induced by track unevenness: Analysis of the excitation mechanism. J. Sound Vib. 2011, 330, 3744–3765. [Google Scholar] [CrossRef]
Giorgio, D.; Cheli, F.; Collina, A.; Corradi, R.; Melzi, S. The Development of a Numerical Model for Railway Vehicles Comfort Assessment through Comparison with Experimental Measurements. Veh. Syst. Dyn. 2002, 38, 165–183. [Google Scholar] [CrossRef]
Lin, J.; Wang, K.; Zhai, W. Impact vibration behavior of railway vehicles: A state-of-the-art overview. Acta Mech. Sin. 2021, 37, 1193–1221. [Google Scholar] [CrossRef]
Gongquan, T.; Wen, Z.; Jin, X.; Yang, X. Polygonisation of railway wheels: A critical review. Railw. Eng. Sci. 2020, 28, 317–345. [Google Scholar] [CrossRef]
Ren, Z.S. An investigation on wheel/rail impact dynamics with a three-dimensional flat model. Veh. Syst. Dyn. 2018, 57, 269–285. [Google Scholar] [CrossRef]
Auersch, L. Different Types of Continuous Track Irregularities as Sources of Train-Induced Ground Vibration and the Importance of the Random Variation of the Track Support. Appl. Sci. 2022, 12, 1463. [Google Scholar] [CrossRef]
Aris, T.; Berg, M.; Stichel, S. Analysing the correlation between vehicle responses and track irregularities using dynamic simulations and measurements. Proc. Inst. Mech. Eng. Part F J. Rail Rapid Transit. 2020, 234, 170–182. [Google Scholar] [CrossRef]
Wen, R.K. Dynamic response of beams traversed by two-axle loads. J. Eng. Mech. Div. 1960, 86, 91–112. [Google Scholar] [CrossRef]
Švejnoch, V. Theory of Railway Vehicles; ČVUT Praha: Prague, Czech Republic, 1991. (In Czech) [Google Scholar]
Zhai, W. Vehicle–Track Coupled Dynamics; Springer: Singapore, 2020; ISBN 978-981-32-9282-6. [Google Scholar] [CrossRef]
Pacejka, H.B. The Dynamics of Vehicles on Roads and on Railway Tracks. In Proceedings of the IUTAM Sympozium; Pelft: Praha, Czech Republic, 1975. [Google Scholar]
Vlk, F. Motor Vehicle Dynamics; Nakladatelství Vlk: Brno, Czech Republic, 2001. (In Czech) [Google Scholar]
Timoshenko, S.P. On the forced vibrations of bridges. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1922, 43, 1018–1019. [Google Scholar] [CrossRef] [Green Version]
Soukup, J.; Volek, J. Vibration of Mechanical Systems—Vehicles; University of Jan Evangelista Purkyně: Ústí nad Labem, Czech Republic, 2008; Volume 227, ISBN 978-80-7414-020-4. [Google Scholar]
Bauer, R. Erzeugung von Räumlichen Fahrbahnmodellen für Komfortunter-Schungen im Labor. Diploma Thesis, Steyer-Daimler-Puch AG, Steyr, Austria, 1987. [Google Scholar]
Ellis, J.R. Road Vehicle Dynamics; USA, 1989. [Google Scholar]
Festa, M.; Stalter, F.; Tavornmas, A.; Gauterin, F. Human Response to Vehicle Vibrations and Acoustics during Transient Road Excitations. Vibration 2021, 4, 357–368. [Google Scholar] [CrossRef]
Abuabiah, M.; Dabbas, Y.; Herzallah, L.; Alsurakji, I.H.; Assad, M.; Plapper, P. Analytical Study on the Low-Frequency Vibrations Isolation System for Vehicle’s Seats Using Quasi-Zero-Stiffness Isolator. Appl. Sci. 2022, 12, 2418. [Google Scholar] [CrossRef]
Volek, J.; Soukup, J.; Skočilasová, B. Investigation of the Vibration of a Spatially Elastic Rigid Plate Considering the Effects of Asymmetry, Linear Damping and the General Excitation Function I: Internation Seminar Dynamic of Rigid and Elastic Bodies 2003, Ústí n. L. 12–13; 2003; Volume 11, pp. 83–90, UTRV UJEP Ústí n. L.; ISBN 80-7044-492-4. (In Czech) [Google Scholar]
Shi, H.; Wu, P.; Luo, R.; Zeng, J. Estimation of the Damping Effects of Suspension Systems on Railway Vehicles Using Wedge Tests. Proc. Inst. Mech. Eng. Part F J. Rail Rapid Transit. 2016, 230, 392–406. [Google Scholar] [CrossRef]

Figure 1. Simple general model of the vehicle.

Figure 2. The scheme of the model.

Figure 3. Chassis support: (a) Original springs; (b) substituted springs.

Figure 4. (a) Bogie fixed to the car platform; (b) loading unit.

Figure 5. Wedge locations at rail.

Figure 6. Sensor locations.

Figure 7. Acceleration sensor location: (a) at frame of backward car; (b) at platform of car—transversal axis.

Figure 8. Measured combination of loading and kinematic excitation.

Figure 9. Location of relative vertical deflection and vertical acceleration sensors at selected points. B—handbrake, 1L/R—chassis frame vertical deflection sensor, first wheelset, left/right side, 2L/R—chassis frame and body vertical deflection sensor, left/right side, 3L/R—car body vertical acceleration sensor, left/right side.

Figure 10. Relative vertical displacement between chassis frame and wheelset—variant A–I.

Figure 11. Relative displacement between car body and chassis frame—variant A–I.

Figure 12. Relative displacement between chassis frame and wheelset—variant A–IV.

Figure 13. Relative displacement between car body and chassis frame—variant A-IV.

Figure 14. Relative displacement between chassis frame and wheelset—variant D-I.

Figure 15. Relative displacement between car body and chassis frame—variant D-I.

Figure 16. Relative displacement between chassis frame and wheelset—variant D-IV.

Figure 17. Relative displacement between car body and chassis frame—variant D-IV.

Figure 18. Relative displacement between chassis frame and wheelset—variant B-III.

Figure 19. Relative deflections between car body and chassis frame—variant B-III.

Figure 20. Vertical acceleration in the middle of the car body on the left and right side—variant B-III.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Klimenda, F.; Skocilas, J.; Skocilasova, B.; Soukup, J.; Cizek, R. Vertical Oscillation of Railway Vehicle Chassis with Asymmetry Effect Consideration. Sensors 2022, 22, 4033. https://doi.org/10.3390/s22114033

AMA Style

Klimenda F, Skocilas J, Skocilasova B, Soukup J, Cizek R. Vertical Oscillation of Railway Vehicle Chassis with Asymmetry Effect Consideration. Sensors. 2022; 22(11):4033. https://doi.org/10.3390/s22114033

Chicago/Turabian Style

Klimenda, Frantisek, Jan Skocilas, Blanka Skocilasova, Josef Soukup, and Roman Cizek. 2022. "Vertical Oscillation of Railway Vehicle Chassis with Asymmetry Effect Consideration" Sensors 22, no. 11: 4033. https://doi.org/10.3390/s22114033

APA Style

Klimenda, F., Skocilas, J., Skocilasova, B., Soukup, J., & Cizek, R. (2022). Vertical Oscillation of Railway Vehicle Chassis with Asymmetry Effect Consideration. Sensors, 22(11), 4033. https://doi.org/10.3390/s22114033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu