# Mitigation of Railway Traffic Induced Vibrations: The Influence of Barriers in Elastic Half-Space

Authors: Michele  Buonsanti; Francis Cirianni; Giovanni Leonardi; Adolfo Santini; Francesco Scopelliti
Journal: Advances in Acoustics and Vibration
 (2009)

---

## Abstract

In this paper, the problem of vibrations induced by trains and their propagation through the soil is studied. Particular attention is focused on the vibration induced by trains in motion and on the effects of such vibrations on the foundations of buildings in proximity of the tracks. The interaction between propagating waves induced by trains in motion and buildings foundations is a problem which does not admit a straightforward analytical solution; thus a solution is given by the use of a model based on the finite elements method. Firstly, we analyze the theoretical aspects of the problem by considering constant or harmonic loads moving along a straight railway track; then, we define a transfer function soil-railway and the response function of the entire system. The study aims to address the wave propagation in an elastic semi-space and the presence in the ground of a discontinuity element, such as a barrier of a given depth is considered. The efficiency variation of barriers is analyzed in function of the different materials used, and different numerical simulations are analyzed in order to study how the wave propagation and the track-soil interaction are influenced by the membrane, seen as damping barrier.


---

## Body

## 1. Introduction

The  foundations  of structures are subject to vibrations due to moving masses, as vehicles in transit, vibrations which can cause damage locally to the foundations, as they can also travel up to the structures in elevation (Figure1).Figure 1
Transmission of induced vibrations.In this paper we consider the effects of vibrations due to the transit of locomotives on a railway, analyzing how the presence of material discontinuity, in an elastic semispace has an influence on the transmission of waves generated by a train passing by. The presence of barriers, made in polymeric material, rubber, or in concrete, seems to cause reflection and refraction of the vibrations generated in the ground. There are some references in scientific literature. Francois et al. [1], approached interaction between ground and infrastructures, in presence of vibrations generated by road traffic, while Sheng et al. [2], modelled the propagation of vibrations due to constant or harmonic loads, produced by the movement of trains.Numerical aspects are developed in detail in Lombaert et al. [3] and Fiala et al. [4], while in Hubert et al. [5], a detailed analysis using the Boundary Element Method, in the dominion of frequency and of time, is developed on a model of elastic semispace. In the work of Andersen and Nielsen [6], a model for the reduction of vibrations of the ground through the use of barriers dipped in the ground and laid between rail track and structure is proposed. In the paper, after some general concepts of theory of the vibrations, the aim is set on the problem in homogenous elastic and isotropic semispace, where  a two-dimensional constitutional discontinuity is localized, representing the artificial barrier in the hypothesis set for different materials (rubber, polyurethane, and concrete).Once the case is defined in its theoretical formulation, a numerical simulation with the support of computer application is carried out. What we want to analyze is the variation of the mechanism of propagation of the vibration waves which happens in the semispace, regarding the conditions of quiet, and anyway without the influence of a traveling load. In particular, we observe as the nature of the material used for the realization of the barrier vary,  how the conditions of propagation of the waves vary, and their influence, in an equilibrium state, on the contour conditions of an eventual foundation laid under the barrier.


## 2. Vibrations Theory

In  regard to the fundamental aspects in the theory of the vibrations, in all this paper we follow the lines expressed in Hartog [7] and Clough and Penzien [8]. Our problem can have a first general organization in the context of the response for harmonic type loads. It is possible to characterize the following fundamental equations in the case of systems without damping: (1)mü(t)+ku(t)=p0sinωt,
and with damping,(2)mü(t)+cu(t)+ku(t)=p0sinωt.
The solution of (1) is admitted, from the combination of a complementary solution with a particular solution, in the shape(3)u(t)=Asinωt+Bcosωt+(p0βκ)(1-β2)sinωt.
In the case of (2) the general solution is of the type (4)u(t)=e-ξωt(AsinωDt+BcosωDt)+p0k×1(1-β2)2+(2ξβ)2[(1-β2)sinωt-2ξβcosωt].
The first term of (4) represents the transitory result, while the second term represents the stationary result, which has the pulsation of the solicitation but is out of sync with the last wave. The factor of dynamics amplification D is defined as the relationship between the amplitude of the stationary response and the static movement produced from the external load. In a generalized manner D is function of the damping ratio and of the frequency relationship and usually it tends to infinite when damping is absent.


## 3. Theoretical Aspects in the Model  Characterization

In an elastic semispace we consider that the source of disturbance (train) generates a set of waves which are distributed in the plane system and therefore we can classify the waves plane and classify our problem as a motion of waves within a semi-infinite body, which is assumed, for simplicity, elastic, homogenous, and isotropic. This type of phenomenon is shown in its generality and specificity in the case of the seismic waves as treated in Boschi and Dragoni [9]. We then consider a longitudinal elastic plane wave that propagates in x-direction, if we assume the movement u associated to it as a periodic function of the time, u can be placed in the form(5)u(x,t)=xacos(kx-ωt+φ),
where a is the amplitude, k the wave number, φ a constant.Using the relation between the exponential and trigonometrical functions, and the complex notation, each component of the movement, for a monochromatic wave, can be defined as(6)u(x,t)=Aei(kx-ωt),
that represents an oscillation of amplitude A and wavelength λ=2π/κ, which propagates in the direction of the wave vector k at constant speed (phase speed):(7)vf=ωκ,
any system of waves can be considered as a superimposition of monochromatic plane waves that move with the group velocity: (8)vg=dωdκ.
Equation (8) represents the speed with which the energy is transported and, if vf  is different from vg, the constant phase surfaces move from an extremity to another in the package of waves while this last one is the same one in motion.Given such assumption, we consider the elastic semispace as composed of different homogenous elements. Let’s assume to form a put in elastic discontinuity, in the cross-sectional sense to the semispace, so to be placed in cross-sectional sense to the wave movement. This eliminates the direct flow of the waves from the source to the site of the structure placing itself as element of filter to the vibrations generated from the movement of the train. Regarding the interaction with the wave, or the set of waves, the discontinuity assumes the role of reflection element and refraction for two-dimensional waves. If we consider two elastic semispaces joined with an element of material discontinuity, we may say thatρ,α,β are, respectively, the density and the speed of the j elastic waves in the first semispace while ρ*, α*, β* are the same amounts in the second semispace (Figure 2).Figure 2
Incident wave and reflection.Refraction and reflection of the elastic waves will happen on the discontinuity surface, longitudinal ones, and of cross-sectional ones which can transform one in the other when inciting on the surface that separates the two semispaces. In accordance with Graff [10], the coefficient of reflection or refraction cij  assumes the form: (9)cij=2ρiβicosθiρiβicosθi+ρj*βj*cosθj.
The  contour conditions, on plan x3=0 are the continuity of stress and the movement.We assume that given a monochromatic wave incident with frequencyω, number of wave k1  and travelling unitary amplitude in the semispace in direction-x3, it produces one cross-sectional movement, in direction x2 expressed in the form: (10)u(x3,t)=exp[ωt+k1(x3-h)],x3≥h,
after the phenomenon of refraction and reflection the field of movements takes the form:(11)u(x3,t)=exp[ωt+k1(x3-h)]+R*exp[ωt-k2(x3-h)],x3≥h,
where R* is a coefficient that holds account of all the waves that travel in positive direction x3. The contour conditions require that traction on free surface x3 is equal to zero, x3=0, and that the movement and the traction are continuous on the surface of discontinuity x3=h. In agreement with Graff [10], given Z1,Zthe2 as stiffness of the two semispace and A(ω) the coefficient of the exponential, the module of the latter supplies the amplification of the waves: (12)|A(ω)|=2sqr[cos2(ωh/βj)+(Z2/Z1)2sin2(ωh/βj)].
Equation (12) is a periodic function of the frequency ω  with maximums in correspondence of the values: (13)ωn=(2n-1)(π2)(βjh).
Figure 3 shows the diagram of amplification of the waves.Figure 3
Amplification waves forZ1/Z2=4.


## 4. Modeling and Simulation: The Track-Subsoil Model

In this paragraph we show a numerical analysis, using a software for defined elements, on a two-dimensional modeling of the studied system.


### 4.1. Model Characteristics

In the proposed model the rail track is the conventional one, and the elements that compose it are (Figure4)(i)
a protective compacted layer 30 cm thick of sand/gravel (called, in the Italian railways,supercompact);(ii)
a 12 cm subballast layer made of bituminous (hot-mix asphalt) concrete;(iii)
a traditional ballast layer 35 cm thick;(iv)
an embankment 60 cm high.Figure 4
Modelling of the railway embankment.The dimensional characteristics of the elements above indicated are those requested by the Italian standard (RFI) for high speed lines.Moreover, to analyze the propagation process of the vibrations in the soil, it has been considered a thickness of the support soil of the embankment equal to a total 5 m.The materials properties used in the model were derived from tests and available experimentations in literature. Table1 lists the physical-mechanical characteristics for the soil and track structure.Table 1
Properties of structural materials.
Mechanical characteristicsRail (UIC60)Sleeper (concrete)Traditional ballastHMA subballastProtective layerEmbankmentDensityρ (kg/m3)785024001250220020001000ModulusE (MPa)21000030000130600016080Poisson’s Ratioν0.300.150.300.400.450.50Finally a discontinuity (10 cm thick and 2 m deep) in the semispace at a given distance from the rail track (2 m from the embankment) is considered, considering the passage of a high speed locomotive (ETR 500) with 4 wheels for car.The effects of discontinuity is estimated considering three different materials whose properties are shown in Table2.Table 2
Properties of barriers.
Mechanical characteristicsDensityρ (kg/m3)ModulusE (MPa)Poisson’s ratioνDumping ratioηConcrete2500250000.150.04Polyurethane (elastomer)1170250.500.08Rubber chip500200.250.10


### 4.2. Finite Element Analyses of Train-Induced Ground Vibrations

The finite element simulations of train-induced ground vibrations was developed, using the ADINA software, considering a two-dimensional finite element model perpendicular to the track (two-dimensional perpendicular model).The two-dimensional model is composed of 7517 elements and 8533 nodes (Figure5). The elements consisted of four-node solid element for ballast, HMA subballast, protective layer and embankment and nine-node solid element for rail and sleeper. Furthermore, the rail pad was modeled with parallel discrete springs and dampers (spring constant K=11.2·108N/m and damper constant c=12·104Ns/m).Figure 5
Finite element mesh of the two-dimensional model.For these elements, under the accepted hypothesis of viscose-elastic behavior, the input parameters for the materials characterization, are the modulus of elasticity (E), Poisson’s ratio (v), density (ρ), α and β Rayleigh’s coefficients for the definition of the damping matrix.The acting forces are a sequence of axial loads, moving like the train (the loads are similar to those of an ETR 500 locomotive) (Figure6).Figure 6
Deflection of the track structure from one wheel.The forces are applied using a time function that represents the time history of the force in the considered node. During the simulations the loads were applied four times, reproducing the passing of four axes. The loads can be thought of as triangular pulses distributed over the wheel-rail surface of contact as shown in Figure7.Figure 7
Load distribution of a moving wheel travelling at speed of 200 km/h.All the finite element analyses in this study were performed in the time domain. The time step of the analyses was fixed and set to 0.001 second.Figure8 shows the variation of velocity (y-direction and z-direction) of the middle node under truck and of a node at the distance of 1.5 m from the embankment. The orbits of vertical to horizontal particle motions were mostly counter clockwise; this could indicate the presence of Rayleigh waves [11].Nodes velocity (m/s) responses for a train speed of 200 km/h.
(a)(b)


## 4.1. Model Characteristics

In the proposed model the rail track is the conventional one, and the elements that compose it are (Figure4)(i)
a protective compacted layer 30 cm thick of sand/gravel (called, in the Italian railways,supercompact);(ii)
a 12 cm subballast layer made of bituminous (hot-mix asphalt) concrete;(iii)
a traditional ballast layer 35 cm thick;(iv)
an embankment 60 cm high.Figure 4
Modelling of the railway embankment.The dimensional characteristics of the elements above indicated are those requested by the Italian standard (RFI) for high speed lines.Moreover, to analyze the propagation process of the vibrations in the soil, it has been considered a thickness of the support soil of the embankment equal to a total 5 m.The materials properties used in the model were derived from tests and available experimentations in literature. Table1 lists the physical-mechanical characteristics for the soil and track structure.Table 1
Properties of structural materials.
Mechanical characteristicsRail (UIC60)Sleeper (concrete)Traditional ballastHMA subballastProtective layerEmbankmentDensityρ (kg/m3)785024001250220020001000ModulusE (MPa)21000030000130600016080Poisson’s Ratioν0.300.150.300.400.450.50Finally a discontinuity (10 cm thick and 2 m deep) in the semispace at a given distance from the rail track (2 m from the embankment) is considered, considering the passage of a high speed locomotive (ETR 500) with 4 wheels for car.The effects of discontinuity is estimated considering three different materials whose properties are shown in Table2.Table 2
Properties of barriers.
Mechanical characteristicsDensityρ (kg/m3)ModulusE (MPa)Poisson’s ratioνDumping ratioηConcrete2500250000.150.04Polyurethane (elastomer)1170250.500.08Rubber chip500200.250.10


## 4.2. Finite Element Analyses of Train-Induced Ground Vibrations

The finite element simulations of train-induced ground vibrations was developed, using the ADINA software, considering a two-dimensional finite element model perpendicular to the track (two-dimensional perpendicular model).The two-dimensional model is composed of 7517 elements and 8533 nodes (Figure5). The elements consisted of four-node solid element for ballast, HMA subballast, protective layer and embankment and nine-node solid element for rail and sleeper. Furthermore, the rail pad was modeled with parallel discrete springs and dampers (spring constant K=11.2·108N/m and damper constant c=12·104Ns/m).Figure 5
Finite element mesh of the two-dimensional model.For these elements, under the accepted hypothesis of viscose-elastic behavior, the input parameters for the materials characterization, are the modulus of elasticity (E), Poisson’s ratio (v), density (ρ), α and β Rayleigh’s coefficients for the definition of the damping matrix.The acting forces are a sequence of axial loads, moving like the train (the loads are similar to those of an ETR 500 locomotive) (Figure6).Figure 6
Deflection of the track structure from one wheel.The forces are applied using a time function that represents the time history of the force in the considered node. During the simulations the loads were applied four times, reproducing the passing of four axes. The loads can be thought of as triangular pulses distributed over the wheel-rail surface of contact as shown in Figure7.Figure 7
Load distribution of a moving wheel travelling at speed of 200 km/h.All the finite element analyses in this study were performed in the time domain. The time step of the analyses was fixed and set to 0.001 second.Figure8 shows the variation of velocity (y-direction and z-direction) of the middle node under truck and of a node at the distance of 1.5 m from the embankment. The orbits of vertical to horizontal particle motions were mostly counter clockwise; this could indicate the presence of Rayleigh waves [11].Nodes velocity (m/s) responses for a train speed of 200 km/h.
(a)(b)


## 5. Results

The graphs of the outcomes of the various analyses are shown in Figures9 and 10.Figure 9
Trend of the velocity with the increase of the distance from embankment at the depth of 1 m.Figure 10
Trend of the acceleration with the increase of the distance from embankment at the depth of 1 m.To  estimate the effectiveness of the barrier we have estimated the dynamic characteristics (speed and acceleration) before and after the discontinuity, in correspondence of an eventual foundation plan of a building placed at a distance of 1.5 meters from the barrier.As a general result of the analyses, as illustrated in Figures9 and 11, for the considered loads, it is concluded that, when the load is travelling with the speed of 200 Km/h, concrete barriers provide a more efficient vibration shielding than the rubber and polyurethane barriers.Trend of displacement, velocity, and acceleration in two nodes located before and after the barrier.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Figure11 shows, for the three kinds of barrier, the course of the displacement, velocity, and acceleration in time.The results given in Figures9 and 11 supply useful information for the making of barriers, especially when it is required to respect regulation’s limit for vibration, as expressed in terms of displacements and velocities (i.e., Italian regulation UNI 9916/2220).


## 6. Conclusions

The results of the proposed elaborations at the finite elements  have been finalized to the assessment of the incidence of the barrier on the vibration state induced from the passage of a high speed locomotive. The response of three different materials for the barrier  was confronted and the following conclusions can be made:(i)
concrete seems to provide a better reduction of the vibration. In spite of the greater density of the material it involves an increase of the reflection phenomena and a consequent increase of the phenomenon at the top side of the barrier (Figures9 and 10);(ii)
polyurethane and rubber chip materials seem to respond in a similar way to the solicitations; however their damping contribution, in the analyzed geometric configuration, do not contribute meaningfully.The results obtained were given with the use of different materials for the barrier, in identical geometrical conditions, in size and position of the proposed damping barriers. A follow up to the presented study, requiring further analyses, is in course, to estimate the influence of the position and of the thickness of the barriers on their dampening ability.



---

*Source: 956263-2009-07-15.xml*