#### Abstract

The paper tries to clarify the problem of solution and interpretation of railway track dynamics equations for linear models. Set of theorems is introduced in the paper describing two types of equivalence: between static and dynamic track response under moving load and between the dynamic response of track described by both the Euler-Bernoulli and Timoshenko beams. The equivalence is clarified in terms of mathematical method of solution. It is shown that inertia element of rail equation for the Euler-Bernoulli beam and constant distributed load can be considered as a substitute axial force multiplied by second derivative of displacement. Damping properties can be treated as additional substitute load in the static case taking into account this substitute axial force. When one considers the Timoshenko beam, the substitute axial force depends additionally on shear properties of rail section, rail bending stiffness, and subgrade stiffness. It is also proved that Timoshenko beam, described by a single equation, from the point of view of solution, is an analogy of the Euler-Bernoulli beam for both constant and variable load. Certain numerical examples are presented and practical interpretation of proved theorems is shown.

#### 1. Introduction

The problem of track dynamic response under moving load is the subject of many theoretical and experimental investigations. Under some assumptions, the beam on elastic foundation can be considered as a typical track model. It is worth mentioning the initial study of beams on the Winkler foundation subjected to a concentrated force moving with constant speed that was initiated by Timoshenko . The first solution to a simple stationary case of the Euler-Bernoulli beam on elastic foundation was properly obtained by Ludwig . Mathews formulated and partly solved the case of moving and oscillating force . The case of varying moving force was studied, for example, by Fryba  and by Bogacz and Krzyzynski .

Many papers are devoted to study various effects of generalized models:(1)Analysis of Timoshenko beam under moving constant and varying loads (presented, e.g., in )(2)Analysis of a beam on elastic half-space [10, 11](3)Response of beam on nonlinear foundation (e.g., [8, 1214])(4)Dynamic response of beam on random foundation; see (5)Dynamic response of track as multilayered structure (see [18, 19], analytical approach; , numerical approach);(6)Analysis of set of distributed moving forces, described by Heaviside functions (e.g., ), rectangular function [9, 19], cosine square formula [8, 12, 13], or Gauss function (7)Effect of axial force on dynamic response [19, 23](8)Analysis of set of forces varying harmonically and associated with track imperfections including the phase of sine function for particular axles (numerically ) and analytical approach [13, 19]

In all described generalizations of classical approach, the track response model was composed of rail (as the Euler-Bernoulli or Timoshenko beam) and viscoelastic or elastic foundation. The sleepers and the ballast were modelled as additional layers.

Analytical closed form solutions were obtained for infinite Euler-Bernoulli beam on elastic foundation under single concentrated load moving along the beam with constant speed using so-called matching conditions for homogeneous solution in the point related to a load position [2, 3]. Using the same approach, the damping properties of foundation or axial force in beam (rail) were also included in the closed form solutions (comp., e.g., [4, 5, 23]). For oscillating force or set of oscillating forces, the solution was obtained by using the Fourier transform and the inverse Fourier transform [5, 6, 8]. Difficulties arising from the integration of solution in the frequency domain lead to the various method of approximation of the transformed function. An interesting approach based on wavelets application is developed by Koziol and coauthors (comp. [8, 10]). Applications of Fourier series to obtain the solution of railway track response are practically limited to the case of bridges vibrations so far [4, 22]. An example of use of the Fourier series in bounded interval to study the railway track response is presented in . It is shown that solutions obtained in this manner can approximate the track response with very high accuracy depending on a number of Fourier series coefficients and the assumed length of interval.

Analysis of Timoshenko beam is mainly carried out for a set of equations describing coupled beam displacements and rotations [6, 8]. The characteristics of the dynamic response of Timoshenko beam described by a single equation are rarely analysed. Hunt  presents numerical approach to solve the inverse Fourier transform in the case of certain simplified equation for a beam on viscoelastic foundation loaded by a single oscillating force. In , the problem of critical speeds for both constant and varying distributed loads is studied analytically without detailed interpretation of final solution. Single equation of dynamics of Timoshenko beam without foundation was analysed by Majkut .

The above-mentioned papers present many interesting results. However, one can observe a lack of simple interpretation of the dynamic railway track response under moving distributed load in terms of differences between static and dynamic solutions and also between the Euler-Bernoulli and Timoshenko beams. This paper tries to clarify the formulated problem described by linear equation of rail motion. The equivalence between static and dynamic railway track response for foundation without damping is proved. Using the Fourier series in finite interval for both distributed load and solution it is shown that damping properties of rail foundation can be interpreted as additional, substitute load. The equivalence of dynamic track solution, in terms of mathematical physics equations, between the Euler-Bernoulli and Timoshenko beams is also proved. Numerical examples are presented leading to practical interpretation of the formulated theorems in wide range of train speed and foundation damping properties. The introduced set of theorems can be treated as an important contribution to a proper arrangement and classification of knowledge that can be recognized as an essential part of fundamentals of the theory of railway track dynamics. One can observe a lack of such papers trying to formulate basics of railway track analysis in a way similar to subjects recognized more systematized research fields.

#### 2. Track Response to a Set of Constant Forces Moving in Longitudinal Direction

##### 2.1. Rail Modelled by the Euler-Bernoulli Beam Equation

Basic equation of motion of track modelled by the Euler-Bernoulli beam on viscoelastic foundation has the following form (comp., e.g., ):where is rail Young modulus [], is moment of inertia of beam (rail) in vertical plane [], is axial force in rail (positive sign means compressive force) [], is unit mass of beam [], is damping coefficient of foundation [], is foundation stiffness [], and is distributed load [].

In the moving coordinate system (), if load is constant in time, (1) may be transformed to an ordinary differential equation:

Theorem 1. For linear model of the track, described by (1), if load does not change in time, the steady-state solution of (2) is an equivalent static case for damping by the substitution:where is substitute axial force [N]. Effect of foundation damping, with an accuracy determined by approximation of rail displacements by Fourier series in finite interval , can be considered as a static case with substitution (3) and additional (substitute) load described as follows:where is substitute load [N/m]; is function depending on real load q and track model parameters P (i.e., EI; N; m; etc.).

Proof. The proof of the first part of the theorem (see (3)) is relatively simple. The equation of static deflection of beam on elastic foundation (parameter ) with compressive axial force and load has the following form (see, e.g., ):As can be seen, for assumed train speed and unit mass of beam , in the moving coordinate system, the sum of realistic axial force and in (2) for expresses the substitute axial force for the beam on elastic foundation in (5), that is, .
For proving second part of the theorem one assumes that realistic load and the solution can be described by the Fourier series in the assumed interval [], that is,where are the coefficients of Fourier series.
Then (2) leads to the following formulas (for cosine and sine coefficients of the Fourier series):where (comp. (3)). One can derive the form of parameters and by solving the system of (7):whereand the constant value y0 is described by formula (8).
Let other coefficients and , together with the same value (formula (8)), describe the cosine and sine coefficients of solution of (2) in the static case with substituted axial force and the damping coefficient . In this case, the solution of (7) has a simple form:If the solution of dynamic case, for , and static solution are denoted by and , respectively, then the following formula can be written in terms of Fourier series in the assumed interval [] (comp. (6)):and the difference between dynamic and static solutions with isThis difference represents function which can be considered as linked to an additional load in the static case with . Then one can obtain the cosine and sine parts of substitute load in the static case for any index : with constant value equal to zero due to the fact that the difference between dynamic and static solution includes only cosine and sine Fourier series coefficients (the constant is the same for both dynamic and static solutions (comp. (2))). Using formulas (9), (11), and (14), one can obtain the following:The above formulas (15) describe the cosine and sine coefficients of substitute load:It means that substitute load which comes from damping properties of the rail foundation is determined as a function of realistic load (coefficients , and speed ) and model parameters and which describe beam and foundation properties (see notation (10)). The second part of the theorem is then proved.

For practical interpretation of Theorem 1, certain remarks should be formulated.

() The equivalence between static and dynamic response of railway track described by simple model, that is, the Euler-Bernoulli beam on viscoelastic foundation under moving load invariant in time, concerns only mathematical solutions.

() In the case when damping properties of rail foundation can be neglected, due to the load being distributed on very small spans and the distance between wheels being relatively large (practically more than about 5-6 m), dynamic factor for maximum rail displacements can be described by simple formula (based on [19, 23]):where formula (5), is critical value of axial force.

The sign of the substitute axial force depends on relation between realistic compression/tension force in rail and the term which has positive value. Knowing that real axial force in rail depends mainly on the rail temperature changes in relation to the laying rail temperature ,where is coefficient of thermal expansion of rail steel [1/°C] and the part of substitute axial force (depending on speed) can be calculated in relation to the equivalent rail temperature increase . For typical rail, being in use in Europe, that is, 60E1, and unit mass of rail 60 kg/m (with the unit mass of half sleeper, usually equal to around  kg/m), one can obtain the following equivalent rail temperature increase: or ; see . It means that(i)°C or °C for  km/h (27.77 m/s);(ii)°C or °C for  km/h (55.55 m/s);(iii)°C or °C for  km/h (83.33 m/s).

For typical track structure, the critical axial force (formula (17)) is around 25–35 MN. The equivalent rail temperature increase for these values is on the level of 1500°C. Using expression (17), it follows that, for realistic rail temperature increase (maximum 50–60°C) and speed up to 300 km/h, the dynamic factor for track without imperfections and damping is lower than 1%, in the case of typical track structure.

() The dynamic factor associated with damping properties of rail foundation in relation to the static response with substitute axial force cannot be described by simple formula. Also, the substitute load linked to the foundation damping (formulas (15)) cannot be simply explained. Currently, these problems can be numerically solved. The following track and load parameters are assumed in further study :(i)Rail: 60E1, Young modulus,  N/m2; cross-sectional area,  m2; thermal expansion coefficient, 1/°C; unit mass,  kg/m (or  kg/m if rail and a half of concrete sleeper PS-94 type are taken into account); rail temperature increase, °C(ii)Rail foundation: substitute foundation stiffness due to fasteners and sleeper foundation elastic properties  N/m2, damping properties – fasteners  Ns/m2 (reference value obtained during laboratory tests and also and ), sleepers foundation  Ns/m2 (reference value and )(iii)Load: 4 wheels of EMU-250 Pendolino train with configuration: 2.7 m, 4.5 m, and 2.7 m; see , wheel load 80 kN; load distribution represented by the Gauss function with the parameter  m, number of Fourier series coefficients equal to 3000

Figure 1 presents the dynamic factor due to rail foundation damping for two options: higher values describe dynamic factor for viscous properties of fasteners and sleeper foundation properties ; lower values describe dynamic factor for the reference values and . As can be observed, these relations have various curvatures and show sensitivity related to the effect of damping on dynamic factor. In the considered range of operating train speeds and viscous coefficients, the maximum value of dynamic factor can reach 4%.

Figure 2 presents rail deflection lines for two rail foundation damping properties (as shown in Figure 1) and speed of 300 km/h. Figures 3 and 4 show the substitute load associated with similar options of rail foundation damping properties. It can be seen that the substitute load due to damping properties of rail foundation strongly depends on train speed and damping coefficients.

##### 2.2. Rail Modelled by the Timoshenko Beam and Comparative Study of the Obtained Results

It is known that various authors use the term Timoshenko beam for different equations or systems of equations. This problem is more discussed in Appendix, at the end of this paper. On the basis of [7, 25] the following single differential equation of rail vertical motion described by the Timoshenko beam theory can be written:where the used notation is similar to ones used in (1) and, additionally, is rail mass density [kg/m3], is rail cross-sectional shear coefficient, is steel shear modulus [N/m2], is rail foundation viscoelastic reaction [N/m], and is substitute load [N/m].

The foundation reaction and the substitute load are characterized by the following expressions:where [m2] is a rail cross-sectional area.

In the moving coordinate system (; ), for the load constant in time, (19) and expression (20) lead to the ordinary differential equation:whereThe load and rail displacements function are described by Fourier series in the assumed interval [] (see formulas (6)). Following the methodology used in previous subsection, one can obtain the closed form solution in the case of the infinite Fourier series:whereand , , denote Fourier coefficients of distributed load function (comp. formulas (6)). One can observe that formulas (23), with a use of symbols (24), have the same mathematical form as formulas (9) with notations (10). Therefore, using procedure like the one used in the case of the Euler-Bernoulli beam, one can prove relatively easy the following theorem.

Theorem 2. For linear model of the track, modelled by the Timoshenko beam on viscoelastic foundation (see (19)), if load does not change in time, the steady-state solution of (21) is an equivalent static case with substitution (22). Effect of foundation damping , with accuracy determined by approximation of rail displacements using Fourier series in finite interval , can be considered as a static case with substitution (22) and additional (substitute) load described as follows:where

Proof. For a sake of reading simplicity the proof of Theorem 2 is left to readers as simple symbols manipulation using remarks formulated above.

It is also worth formulating the following theorem being a consequence of the studies carried out.

Theorem 3. For the linear model of the track, if load does not change in time, the steady-state solution of the Timoshenko beam equation (expression (21)) is equivalent to the steady-state solution of the Euler-Bernoulli beam (2), for any speed v and foundation damping , with the following substitution:whereEffect of foundation damping for the Timoshenko beam, with an accuracy determined by approximation of rail displacements using the Fourier series in finite interval , is equivalent to the Euler-Bernoulli beam in terms of algebraic equations for cosine and sine part of solution (comp. expressions (7)) with substitution

Proof. The theorem can be relatively easy proved based on the previous studies (see Theorem 1) and, here, similar to Theorem 2, its proof is omitted in order to avoid repetitions. Instead, the practical interpretation is presented as more important for understanding of the theorems formulation in terms of scientific arrangement of the subject.

For practical interpretation of Theorems 2 and 3, certain remarks and conclusions are presented on the basis of numerical examples using the track and load parameters assumed in Section 2.2 for the Euler-Bernoulli beam. Values of additional parameters are taken from : mass density of rail steel  kg/m; shear coefficient of rail cross section ; shear modulus of rail steel  N/m2:(1)The dynamic factor for maximum rail displacements for the Timoshenko beam without damping, in the speed range [0; 300] km/h, is less than 1%, which is similar to the Euler-Bernoulli beam case.(2)Rail foundation damping for the Timoshenko beam has significantly smaller effect in relation to the Euler-Bernoulli beam case. This is shown in Figure 5 in the case of typical damping and high values of damping coefficient (the same values are used in Figure 1). The dynamic factor for maximum rail displacements due to damping reaches a level of a few percent (3-4%) in the case of the Euler-Bernoulli beam and, at the same time, it reaches only 1% in the case of the Timoshenko beam.The peak corresponding to this one observed in Figure 5 for the Timoshenko beam case and higher damping value appears also in the case of the Euler-Bernoulli beam. However, this peak is reached for higher speed, that is, around 400 km/h (150 km/h in the case of the Timoshenko beam). In both cases, the dynamic factor decreases above these speed values but it remains positive. It reaches the level of 0.6 for both the Timoshenko and the Euler-Bernoulli beams around 2000 km/h. It also becomes slightly higher in the case of the Timoshenko beam compared to the Euler-Bernoulli model. For higher speeds, reaching 3000 km/h, the dynamic factor can be estimated around 0.5 in both cases.One should note that, in the case of high damping value, unsymmetrical response of the track in relation to the centre of the load system can be observed. This effect is minimal when one deals with low speeds. For these parameters, the maximal displacements increase along with the speed increase up to 150 km/h (for the Timoshenko beam) and 400 km/h (for the Euler-Bernoulli beam). Above these values, maximal deflections become smaller.In the case of the reference damping value, the dynamic factor increases without any extremal values. One can show that it reaches the level of 1.2 near 700 km/h.(3)The coefficient describing the ratio between maximum rail displacements for the Timoshenko and Euler-Bernoulli beams is on the level of 5.4% for the analysed speed range and it is practically independent of speed. Hence one can conclude that the static effect associated with the shear force can be treated as the main factor characteristic for the analysed model with moving load invariant in time.

This observation can be additionally justified by the fact that the dynamic factor for the Timoshenko beam is relatively small in both cases, that is, with and without damping (comp. Figures 6 and 7).

#### 3. Track Response for a Set of Oscillating Forces Moving along the Track with Constant Speed

##### 3.1. Rail Modelled by the Euler-Bernoulli Beam Equation

In the case of forces oscillating with circular frequency , (2) obtains the following form in the moving coordinate system ():The assumed form of solution is as follows (comp. [5, 19]):One should mention that, in general, a set of forces can be also described as follows:By differentiating (31) and substituting (32), (30) becomes a set of ordinary equations associated with cosine () and sine parts ():The functions , , and are expressed in terms of Fourier series in the assumed interval []:After differentiation and rearranging one obtains the following set of algebraic equations determining unknown parameters , , , and , :where

##### 3.2. Rail Modelled by the Timoshenko Beam Equation

In the case of forces oscillating with circular frequency , (19) and (20) obtain the following form in the moving coordinate system (, ):where the rail foundation reaction and the load are described by the following expression:Under assumption that the load and the unknown function of rail displacements are represented by cosine and sine parts,differentiating functions (39) and applying the following notation:lead to a possibility of transformation of formulas (38a) into the system of two ordinary differential equations related to cosine and sine parts:The right-hand sides of these two equations, representing the load, obtain the following form:The quantities and represent cosine and sine part of realistic load. is defined above (see formulas (40)). The functions , , , and are expressed in terms of Fourier series in the assumed interval [ (see (34)). By substituting the functions , , , and their derivatives to (41) and (42), one can obtain, after rearranging, the following system of algebraic equations for unknown values , , , :whereThe quantities for the right-hand side of (45) are defined by (40), whereas load coefficients (right-hand side of (43) and coefficients in expression (44)) are described as follows:

###### 3.2.1. Equivalence between the E-B and the Timoshenko Beams in the Case of Varying Load

It is easy to observe that solution for rail displacements determined in the case of the Timoshenko beam ((43) and (44) with notation (45) and (46)) has the same mathematical form as solution for rail displacements determined in the case of the Euler-Bernoulli beam (formulas (35) and (36) with notation (37)). Therefore, the following theorem can be formulated.

Theorem 4. For the linear model of the track, if load is varying in time, the steady-state solution for the Timoshenko beam (see (38a)) is equivalent to the steady-state solution for the Euler-Bernoulli beam (see (30)), with accuracy determined by approximation of rail displacements by using Fourier series in finite interval , in terms of the following substitution in matrix (35) and (36):where () defines formulas (45) and aT0, cT0, and () are described by (46).

Proof. The proof is omitted as a direct observation arising from results and analyses presented in Sections 3.1 and 3.2.

Certain remarks can be formulated for practical interpretation of the above theorem:(1)The Timoshenko and the Euler-Bernoulli beams analogy concern only mathematical form of solution of the steady-state response.(2)For comparative analysis of the steady-state response of Timoshenko and Euler-Bernoulli beams under moving load oscillating with circular frequency and having an amplitude , the same track parameters as used in Sections 2.1 and 2.2 are considered. Firstly, the substitute bending stiffness of Timoshenko beam is analysed:This quantity expresses the constant associated with the fourth derivative (comp. (41)). The comparison concerns the bending stiffness for the Euler-Bernoulli beam and formula (48). Figure 8 shows the effect of train speed on the relation . As can be seen, in the case of 60E1 rail typically used in Europe for main railway lines decrease of rail bending stiffness for Timoshenko beam in relation to the bending stiffness for Euler beam is on the level less than 0.3% for speed up to 350 km/h. It means that, in practical calculations, simple rail bending stiffness can be used, also if Timoshenko beam is taken into account.(3)Figure 9 presents the steady-state response of Timoshenko and Euler-Bernoulli beams to a set of distributed forces with amplitudes  kN circular frequencies  rad/s (833.33 Hz) and train speed  km/h. One can see that the difference between displacements of rail modelled by Timoshenko and Euler-Bernoulli beams is quite significant: the maximum value exceeds 10%.

#### 4. Conclusions

In this paper, analytical solutions of railway track response to moving forces are studied in the case of linear models. The obtained results lead to the following conclusions:(1)Dynamic equation of vertical track motion in the case of the Euler-Bernoulli beam resting on elastic foundation and subjected to a constant moving load is equivalent to a static beam equation with substitute axial force for fixed unit mass of beam and train speed. In the case of the Timoshenko beam, the substitute axial load depends additionally on shear properties of rail section, rail bending stiffness, and subgrade stiffness. This conclusion is valid for each method applicable to solution in the moving coordinate system in the case of constant speed, due to its formal justification based on mathematical formulation of the problem taking into account equations of mathematical physics only. Depending on analysed features of the modelled railway system, an appropriate approach based on the Euler-Bernoulli or the Timoshenko beam formulation can be applied. The obtained results allow suggesting that the Euler-Bernoulli beam should be preferred in the scope discussed in this paper, due to its simpler formulation leading to relatively easier computations.(2)Damping properties of track (beam) foundation can be treated as an additional substitute load taking into account a substitute axial force in the static case. The equivalence between dynamic and static responses of beam resting on elastic foundation with damping is clarified in terms of mathematical method of solution for fixed unit mass and train speed. This equivalence is shown by using the Fourier series representation for both the distributed vertical load and the vertical rail displacement. The equivalence remains correct for both track models: the Euler-Bernoulli and the Timoshenko beams.(3)The Timoshenko beam described by a single equation is an analogy of the Euler-Bernoulli beam for both constant and variable moving distributed loads, in terms of form of solution in the moving coordinate system.(4)The set of theorems formulated in the paper can be treated as an important part of fundamentals of the theory of railway track dynamics, being at the same time a trial approach to the systematization of the subject principles.

#### Appendix

Original approach presented by Timoshenko  provides the beam equation considering rotary inertia and shear effect in the following form:where is transverse displacement [m], is Young’s modulus of beam [N/m2], is moment of inertia of cross-sectional area [m4], beam mass density [kg/m3], is cross-sectional area [m2], is shear modulus [N/m2],and is shear coefficient of beam cross section.

Equation (A.1) does not include an external load. Forced vibrations of Timoshenko beam are presented in, for example, . The beam equation is formulated in this case as followsThe left-hand side of (A.2) is similar to this one in (A.1), whereas [N/m] is an external load distribution.

If viscoelastic continuous foundation of beam is taken into account (railway track case), then one usually uses two equations describing the Timoshenko beam (e.g., [26, 29]):where is beam slope due to bending, is elastic coefficient of beam foundation [N/m2], and is damping coefficient of foundation [Ns/m2].

After rearrangement of the two above equations, one can obtain the following equation for Timoshenko beam deflection:One should note that various authors use the term Timoshenko beam for different equations. For example,  describes the Timoshenko beam by using the following expression:It can be seen that four terms related to damping and elastic properties of beam foundation, as well as derivatives of the load, are neglected in the above (A.5) compared to (A.4). Further doubts arise from consideration of the axial force included in Timoshenko beam equation. In , equation of Timoshenko beam without foundation is considered:Here in (A.6), in relation to (A.1)-(A.2) and (A.4)-(A.5), the fourth derivative with respect to time is neglected and a simplified expression of fourth-order mixed derivative is proposed. It can be shown that these changes lead to the error of order around 13% in the case of typical track parameters. In the contrast, in , two equations are considered:After rearrangement of these two equations, one can obtain the part of equation for the beam deflection depending on axial force:Assuming that the axial force comes from the rail temperature changes (in relation to neutral temperature), one can write the following :This shows that the above term is independent of a cross-sectional area of the beam (rail). For typical rail steel parameters:  MN/m2;  MPa  MN/m2; ; [1/°C]; maximal rail temperature increase/decrease = ±50°C; and the thermal expansion coefficient  1/°C, one obtainsThus, the considered terms can be neglected in terms of railway tracks applications. Knowing that steel materials () and geometrical parameters (especially ) may vary for various types of rails, expression (A.8) can be written as follows:One should underline that, in the case of the Euler-Bernoulli beam, the problem of the axial force introduction is clearly interpretable, whereas in the case of the Timoshenko beam the axial force can be introduced in various ways, by inclusion of appropriate term in either equation for vertical deflection or torsional vibrations. Hence it is worth investigating how important is the difference between these two approaches, depending on the application analysed. In present paper, the presented analysis deals with typical parameters of railway track only. One can see that also other effects, for example, those associated with fourth derivative with respect to the time variable or derivatives of external load, are taken into account or omitted depending on the actual Timoshenko beam formulation.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors would like to thank Mr. Dariusz Kudla for his help in the paper preparation.