# Assessing the Open Trenches in Screening Railway Ground-Borne Vibrations by Means of Artificial Neural Network

Authors: Gaetano Di Mino; Marinella Giunta; Ciro Maurizio Di Liberto
Journal: Advances in Acoustics and Vibration
 (2009)

---

## Abstract

Reducing ground borne vibrations in urban areas is a very challenging task in railway transportation. Many mitigation measures can be considered and applied; among these open trenches are very effective.
This paper deals with the study of the effect, in terms of reduction of vertical and horizontal displacements and velocities, of the open trenches. 2D FEM simulations have been performed and several open trench configurations have been analysed varying the main geometric features such as   width and depth, distance from the rail, thickness of the soil layer over the rigid bedrock, type of the ground, ratio between the depth of the trench,  and the thickness of the soil layer. For quantifying the influence of the above specified parameters in reducing ground-borne vibrations an analysis using artificial neural networks (ANNs) has been carried out. Results show that among the geometric parameters the role of the depth of the trench is very significant; however the influence of the depth must be also evaluated in relation to the thickness of the soil layer.


---

## Body

## 1. Introduction

The problem of the ground-borne vibrations induced by railway traffic has received an increasing interest in the recent years becoming a relevant scientific and technical research area.Due to the frequent construction of high-speed railways and mass rapid transit systems worldwide, most highly developed cities or metropolitan areas have encountered the problem that the rail tracks inevitably come cross or close to vibration-sensitive sites, as discussed by Hung and Yang [1, 2].The economical and environmental aspects of the issue require a careful assessment of the problem before the construction of new rail tracks or the upgrading of the existing ones for heavier and faster traffic.The main source of excitation of the track is represented by the vertical force determined by the wheel-rail interaction during the passage of the train. Vibrations can be amplified by the passage of trains due to the surface irregularities of wheels and rails, by the rise and fall of axles over sleepers, and by the propagation of deformation patterns in the track and ground. Such vibrations are transmitted through the track structure, including rails, sleepers, ballast, and sublayers  and propagate as waves through the soil medium [3].The study of the ground-borne vibrations requires the consideration of four main components: the “source,”  that generate the vibration, that is, the excitation caused by the motion of the trains over rails with irregular surfaces, the “propagation path” through the soil medium, the “receiver,” that is, the nearby buildings, and finally the “interceptor,” that is, wave barriers, such as piles, in-filled and open trenches, and  isolation pads.At each component is related one specific phase of the process of the transmission of the vibrations and in particular  the generation, the transmission, the reception, and the interception.In regard to this last phase, three groups of mitigation measures can be considered:(i)
mitigations in the source, that include all types of interventions in the railway structures (active isolation);(ii)
mitigations in the path, such as barriers to the waves propagation from the source to the receiver (in-filled and open trenches, lime cements columns); these barriers provide an active isolation when they are close to the source, and  a passive one when they are far away;(iii)
mitigations to the receiver, that include all the measures aimed at reducing the effects on the buildings (passive isolation) or on the other vibration-sensitive sites.Among the mitigation measures in the path, the open trenches have exhibited a good performance on the screening of the vibrations but, to achieve this effect, it is very important to assign it the proper dimensions.


## 2. Previous Studies

The first experimental surveys on the effectiveness of open and in-filled  trenches were carried out by Barkam (1962) [4]. He showed that their effectiveness increases with increasing the depth and the distance for raising frequencies.In regard to the continuous barriers also Woods [3] showed the influence of the depth and the distance from the vibratory source. He also  demonstrated that the passive isolation is better than the active for screening the P-body and  S-body waves while the active isolation is more suitable for the screening of the Rayleigh waves; moreover he highlighted that the most relevant geometric parameter in the screening process is the ratio between depth barrier and Rayleigh wavelength.Several experimental surveys stated that the best screening performance takes place when the depth of the trench is equal to the wavelength while the width of the trench is small [3–10].Numerous researches on the effectiveness of the barriers have been carried out  using FEM and BEM modelling  [11–16].In particular, Beskos et al. [11] employed boundary element method (BEM) to study open and in-filled trenches as well as pile wave barriers. Ahmad and Al-Hussaini [13]proposed a simplified design methodology for vibration screening of machine foundations by trenches using 2D BEM.Adam and Estorff (2005) [6] inspected the effectiveness of open and in-filled trenches in reducing the six-storey building vibrations due to passing trains using a two-dimensional FEM analysis.Yang and Hung (1997) [1, 2] combined the finite and infinite elements to investigate the effect of trenches and elastic foundations in reducing train-induced vibration.Few works focus on the effectiveness of barriers in full scale.The most important application is based on Gas Cushion Method that consists in a vertical panels filled of gas and flexible cushion with very low impedance installed in a trench having a great depth [17]. In Gnarp, Sweden 1984, the first application of these barriers was carried out in order to screen a sensitive building from railway ground borne vibrations. The effectiveness of this application reached an attenuation value of about 70%.


## 3. Problem Formulation and Basic Assumptions

In this work the screening performance of open trenches, excavated in a soil layer over a bedrock, has been studied. By applying a 2D FEM model, an extensive analysis has been carried out to the aim of determining the contribution of the main geometric parameters of the trenches in the interception of the vibrations. The geometric characteristics assumed to be variable parameters are  width (w), depth (d), distance from the rail (l), distance from the vibration-sensitive site (s), thickness of the soil layer over the rigid base (h), and ratio between the depth of the trench and the thickness of the soil layer (d/h) (see Figure 1).Figure 1
Schematic of the position of the open trench.The study has been carried out considering different scenarios obtained by varying the above specified geometric parameters of the trenches and the mechanical and geometric characteristics of the soil medium. The values attributed to the geometric features are summarised in Table1, while the geometric conditions of the soil medium have been set up according to the Table 2.Table 1
Geometric parameters of the trenches.
Widthw (m)Distance from the raill (m)Distance from the sites (m)0.58191.012151.516112.5207—243Table 2
Depth of the trench and ratio depth/thickness of layer.
ParametersAssumed valuesd (m)2, 5, 8, 11,14d/h1/5, 1/2, 2/3Two types of unsaturated soil, labelled, according to the classification of Eurocode 8, respectively, ground type C (medium dense sand) and ground type D (firm clay), have been considered. These soils have been characterised by  the shear modulus (G), the Poisson’s ratio (ν), the mass density (ρ), the damping ratio (ξ), and the shear and the Rayleigh’s wave velocities (Vs, Vr). The values assigned to each of these parameters are shown in Table 3.Table 3
Mechanical characteristics of the soil medium.
Type   of  soilG  (MN/m2)ρ  (N/m4 sec2)νξVs (m/sec)Vr (m/sec)Clay20.71583.80.30.06114.3108.1Sand1381682.80.30.04286.3270.9The value of the damping ratio for clay and sand has been selected   according to the results of an Italian comprehensive review [18] that analyses a great number of experimental results (44 international researches carried out in the last 30 years) and concludes that  although a considerable scatter of the experimental points exists, for deformations lower than 0.001%, the damping ratio has a tendency to increase according to the plasticity index.By combining all the parameters, 540 scenarios have been obtained and investigated. The propagation and interception of ground-borne vibrations for each scenario have been modelled by means of the finite elements method. The scheme adopted is shown in Figure2. The 2D model has been assumed to be symmetric in respect to the railway line.Figure 2
Modelling scheme.It has been also established that the railway line is sited over an embankment 1.5 m high and with escarpment’s slope of 2/3; the geotechnical parameters of the embankment areG=125 MPa,  ρ=2000 N/m4 s2, and ν=0.33.For each scenario the following responses have been obtained: vertical (dir.z) and horizontal (dir. y) displacements and vertical (dir. z) and horizontal (dir. y) velocities. In particular, these parameters have been calculated in two areas:(a)
in a “vibration-sensitive site,” large 10 m and far away 27 m from the symmetry axis of the railway;(b)
in a zone, in the following called “site close to the trench,”  5 m wide and sited far away 0.5 m from the trench.The dimensions assigned to the two sites are well suited to be comparable with the dimensions of receptors such us buildings.The study of the response in this last site allowed to verify the effectiveness of the open trench as passive isolation.


## 4. Finite Element Modelling

The propagation of vibrations is a typical three-dimensional issue, specially in railway field where the train acts  as a series of incoherent point source rather than a fully coherent line source. So three-dimensional models are certainly more suitable to predict the absolute vibration but require computation times not always consistent with extensive analyses.Andersen and Jones [19] performed analyses 2D-3D FEM and BEM and stated that 2D models provide results qualitatively comparable with those obtained using 3D models in a wide range of frequencies.In the light of this, the 2D model has been considered appropriate to output results consistent with the aim of the present study.


### 4.1. Material Constitutive Model and Type of FEM Elements Applied

Linear viscoelastic constitutive model for the embankment and the soil medium materials has been adopted to carry out the investigation, according to some previous studies [20, 21]. Both the embankment and the soil medium have been meshed with triangular elements having 6 nodes. The mechanical parameters assigned to each element are  the Young modulus E,   Poisson’s ratio ν, and the mass density ρ, while the damping behaviour has been introduced by means of the Rayleigh coefficientsα and β. The damping matrix [C] has been defined on the basis of the quoted Rayleigh coefficients α and β.The presence in the analysed scheme of a rigid base has been taking into account by introducing in the model fixities at   bedrock.


### 4.2. Modelling of the Rail  Source

The railway source has been modelled applying in the surface of the embankment below the ballast a stress having a value equal to the dynamic stress produced in the contact surface ballast-embankment.In particular, according to [6], two concentrated linear loads acting over the rails have been considered. The time history of each load consists of four consecutive impulses; each impulse have a time duration of 0.02 seconds and an amplitude equal to 1000 KN. Since only the embankment has  been modelled, the dynamic pressure produced by the linear loads in the contact surface ballast/embankment, having dimensions of about 2,5m×1m, has been fixed in  400 kN/m2 (Figure 3) while its time history is the one plotted in Figure 4.Figure 3
Loading conditions.Figure 4
Time history of the rail force.The applied load function gives a significant contribution for the frequencies included in the range from 20 to 35 Hz, with a central frequency of about 27 Hz. This range, in the present study-case, is representative of the frequencies produced by the passage of a freight train at the velocity of 60 km/h, and having a boogie axle distance of 9 m. The load can appear overabundant but it takes into account the dynamic amplification factor by the  shortest wavelength [6, 22].


### 4.3. Finite Element Mesh

The definition of the model mesh is in general a compromise between the necessity to restrict the number of the elements and that to limit their dimensions. However, the element’s dimension must be accurately chosen in order to adequately represent the deformed shape associated to the wave lengthsλ at the different harmonic frequencies f imputed.Iffmax is the greatest frequency of the significant components of the Fourier spectrum and Vj is the propagation speed of the waves (P-waves, S-waves, or R-waves) in the material, the dimension of the mesh element for the material should satisfy the following relation:(1)d≤λmink=Vjk·fmax,
where k is a coefficient ranging from 4 to 10 according to the type of finite element and to its shape function [23].In the light of this, the maximum dimension of the element has been defined in the following way:(i)
on the basis of the material properties, there has been calculated the minimum wave propagation velocity for the considered domain; that is, for the superficial elements there has been considered the Rayleigh’s velocityVr, for the elements close to deep source the shear velocity Vs;(ii)
there has been defined the maximum value of the appropriate frequencyfmax; that is, in this study, according to the railway source model,  fmax has been set up equal to 35 Hz;(iii)
there has been calculated, using the expression (1), the minimal dimension of the element.By applying this procedure to the two types of soil medium, there has been obtainedClay⇒Vr=108.16 m/sec; fmax=35 Hz; k=6⇒dA≤0.51 m;Sand⇒Vr=270.94 m/sec; fmax=35 Hz; k=6⇒dS≤1.29 m.By considering these dimensions, the mesh has been made as follows:(i)
in the part of the model where the response is significant to the aim of the study, corresponding to a strip, large about 40 m and 40 m deep, closest to the rail track, a maximum dimension of 0.5 m has been assigned to the elements; this mesh has been designed to be denser in this focal area in order to achieve a good precision;(ii)
in the part of the model included between the quoted part and the lower edge a maximum dimension of 0.5 m has been assigned when the clayey soil has been analysed, while in the case of sandy soil the dimensions of the element gradually increase until 1.5 m near the lower edge;(iii)
externally to the above mentioned areas, the dimensions of the elements increase gradually reaching the value of 3 m in proximity of the vertical boundary far 300 m from the symmetry axis.The linear variation of the mesh dimensions, from the sidey=40 m to the boundary, has been chosen in accordance with other similar works [24, 25].The scheme of the finite element mesh is shown in Figure5.Figure 5
Scheme of FEM mesh.Figure 6
Rayleigh’s parameters determination.It can be observed that outside the zone of interest (40m×40m), the growth of the dimensions of the element  in connection with the time step does not produce effects in the part of the model where is sited the trench; in fact, the vibratory waves, also when undergo reflections in correspondence of the vertical boundary right hand, cannot reach this zone before the end of the analysis.


### 4.4. Courant’s Condition

The time step integration has been assigned taking into account the Courant condition, that defines the maximum time step as [26](2)Vj·Δth=C,C≤1
in which Vj is the velocity of the considered wave, h is the maximum spatial displacement admitted for the perturbation  that in FEM models is generally assumed equal to the dimension of the mesh, and C is the Courant number.It should be noted  that in the problems involving different type of waves, like the propagation of theP and S waves in the elastic space, the time step is imposed by the greatest propagation velocity, while the element’s dimensions depend on the lowest velocity. So, a great difference between the greatest and the lowest propagation velocity requires a small finite element and a short-time step.Having in mind this, in this task there has been assumed thatClay⇒Vp=191.55 m/sec, dA,min=0.5m⇒Δt≤0.0026 sec; there has been chosen Δt=0.002 sec.Sand⇒Vp=479.83 m/sec, dS,min=0.5m⇒Δt≤0.00104 sec; there has been  chosen Δt=0.001 sec.In both cases of clay and sand the FEM simulations have been performed in a time ofT=0.75 sec; this time is sufficient to permit the complete passage of the dynamic perturbation in the examined zone.


### 4.5. Boundary Condition

Boundary conditions have been imposed in the model in consideration of the geometric and mechanical conditions of symmetry and taking into account the need to limit the errors, namely, the difference existing between the theoretical value of the stress in a point of the boundary and the value obtained in presence of the restraints.FE calculations need a self-consistent simulation area to operate on. Especially in the case of wave propagation, special boundary conditions have to be incorporated into the calculation scheme [27]. Many authors have proposed special absorbing boundary conditions (ABCs) for the FE calculation of elastic waves. In these works the boundary conditionssuggested from Lysmer and Kuhlemeyer have been applied [28]:(3)σ=-ρcdVn,τ=-ρcsVt,
where σ and  τ are  the normal and the shear stress,  ρ  the material density, and  cd  and cs the longitudinal and the shear wave velocities of the transmitting media. These equations reveal that the magnitude of these stresses at the boundaries is proportional to the particle velocities in the normal (Vn) and in the tangential (Vt) directions.   Lysmer’s dampers placed on the artificial boundary are effective in reducing unwanted wave reflections if the boundary of the finite element mesh is sufficiently far outward.In the light of this, there has been applied(i)
in the symmetry edge, symmetry restraints like bi-pendulum iny direction;(ii)
in the lower edge, fixities at the bedrock;(iii)
in the right-side edge,   Lysmer’s dampers in the directionsy and z.The right-side edges have been applied at a distance of 300 m from the source; the errors obtained in the boundary have been(i)
ε(σn)=0.10%  and  ε(τ)=0.45%, in the case of the clayey soil medium;(ii)
ε(σn)=0.26% and  ε(τ)=1.12%, in the case of the sandy soil medium.


### 4.6. The Determination of the Rayleigh Damping Coefficientsα  and β

In an FEM model the damping behaviour can be modelled using   Rayleigh’s method. The damping matrix can be obtained from the following relation:(4)[C]=α[M]+β[K],
where [M] and [K] are, respectively, the mass matrix and the stiffness matrix.Theαand β damping coefficients [23, 29] have been determined following the procedure described below and shown in Figure 4.(i)
For each 2D model, a frequency analysis has been carried out and the first resonance frequencyf1 has been calculated.(ii)
There has been calculated the frequencyf2 by means of the relation f2=nf1, where n is the first odd number greater than the ratio fi/f1, being fi the dominant frequency. In this case fi has been fixed equal to 27 Hz, like the medium frequency of the load function.(iii)
Rayleigh’s parameters value has been calculated by establishing that, for frequencies equal tof1 and f2, the damping coefficient of the element must be equal to the value ξs, which is proper to the examined soil. In this way there has been introduced a subdamping for the frequencies ranging from f1 to f2 and a superdamping for the other frequencies not included in this range.For the clayey soil   Rayleigh’s parameters have been calculated assuming that a damping ratioξs   is equal to 0.06, while for the sandy one there has been assumed that  ξs  is  equal to 0.04.In   Figures7 and 8, we can see examples of damping diagrams for clayey and sandy soil medium.Figure 7
Damping ratioξs versus frequency for clay bank of 25 m depth.Figure 8
Damping ratioξs versus frequency for sand bank of 40 m depth.


### 4.7. Validation Model

In order to validate the model, results by Beskos et al. [11] have been considered and used. Results concern an open trench having a depth of 2.8 m (1λr) and a width of 0.28 m (0.1λr); it is made in an elastic half-space at a distance of 10.04 m (5λr) from the axis of a rigid foundation having a width of 0.7 m and subjected to a harmonic load with frequency of 50 Hz.In particular the comparison between the results given by the proposed 2D FEM model and the Beskos one has been performed considering the trend of the attenuation ratioAr versus the normalized distance  y/λr, where λr is the length of   Rayleigh’s waves generated from the source and y the horizontal distance of the generic point from the axis of the source.In each point of both vibration-sensitive site and site close to the trench the attenuation ratioAr  has been defined as the ratio between the maximum value of the kinematic parameters (y and z displacements and y and z velocities) in presence of trench and the same  value without trench:(5)Ar=(maxvalueofthekinematicparametersinapointwithtrench)/(maxvalueofthekinematicparametersinapointwithouttrench).
The attenuation ratio has been calculated not only for the vertical displacement (dir. z), as it is usual in   literature, but also for the horizontal displacement (dir y) and for the horizontal and vertical velocities.Since the sites are not   points but   large zones and considering the evolution of the attenuation ratio increasing distance (Figure7), the effectiveness of the trench in these sites   has been evaluated in terms of average attenuation ratio:(6)Am=1C∫ocAr(y)dy,
where C is the length of the investigated zone (equal to 10 m for vibration sensitive site and 5 m for site close to the trench).As it can see in Figure9, a good agreement between the numerical results from Beskos simulation and the ones obtained with the proposed model has been achieved. The model has given good results also in others experimental applications [29] and thus it has been applied in the present study.Figure 9
Model results versus Beskos results.


### 4.8. Results

From the numerical FEM simulations of the 540 scenarios the values of the displacements and velocities have been obtained; the minimum and the maximum values reached in the significant points are summarised in Table4.Table 4
Maximum and minimum values of displacements and velocities.
Site close to the trenchSoil  typeVelocity  dir.y min max  (m/sec)Velocity dir.z min max (m/sec)Displacement dir.y min max  (m)Displacement dir.z min max  (m)Sand0.130.090.000980.000560.350.250.002520.00619Clay0.160.130.005640.003690.290.480.010840.02509Site vibration-sensitiveSoil typeVelocity  dir.y min max (m/sec)Velocity dir.z min max (m/sec)Displacement dir.y min max  (m)Displacement dir.z  min  max  (m)Sand0.0520.0350.000780.000090.210.200.001980.00341Clay0.0350.0280.001770.000400.270.300.010840.01714In Figure10 two samples of time history of displacement and velocity in a point are plotted.Samples of results in a point in case of sandy soil.
(a)(b)When the horizontal propagation of in-plane waves is considered, the Rayleigh waves and theP-waves are the main causes of the ground motion.At quite distance (5λr) from a vibratory source the ground motion is due to the Rayleigh waves energy.If an open trench is excavated near the source, the incidence of the waves on this obstacle (discontinuity field) gives rise to reflected and transmitted body waves. Behind the open trench, as the distance increases, the transmittedP-waves get partially into Rayleigh waves (phenomenon known as mode conversion). Therefore, the Rayleigh waves are the main cause of  the ground motion [13].According to this phenomenon the results show a great attenuation of displacement and velocity passing from the site close to the trench to the vibration sensitive site (see Table4).To evaluate the effectiveness of the trench there has been adopted the very strict following criterion: a trench is effective in screening the ground borne vibration if the value of the attenuation ratio is lower than 0.4 for every considered displacement  and velocity.The analysis of the FEM results for the two types of ground has been performed having in mind this criterion.The results give some indications for the design of the open trenches in both cases of site close to the trench and vibration-sensitive site (see Tables5 and 6).Table 5
Design parameters for trench—site close to the trench.
d/hClayey soil layer[λr=4,00m]Sandy soil layer[λr=10,00m]dlwdlw1/5———0.2λr2λr≤0.1λr1/21.25λr>3λrall0.5λr≥1.2λrall2/3<2.75λr6λrall0.2λrallallTable 6
Design parameters for trench—vibration-sensitive site.
d/hClayey soil layer[λr=4,00m]Sandy soil layer[λr=10,00m]dlwdlw1/5———0.2λr0.8λr0.05λr1/21.25λr>3λrall0.5λr0.8λr0.05λr2/30.5λrallall0.2λr≥ 0.8λrallThe effect of frequency is fully taken into account by normalizing the open trench dimensions with respect to the Rayleigh wavelength.The obtained results are very strict in regard to the dimensions and the position of the trench because these features depend on the ratiod/h and so on the finite value of the layer thickness.For instance, in clayey soil results show that the trench is effective if it is 5 m deep and far away at least 12 m from the source, for thickness layer equal to 10 m or 2 m deep and far away at least 24 m from source, and  for layer thickness equal to 3 m.Some results of the parametric studies carried out are plotted in Figure11. It can be seen that the efficiency of the trench is in general independent from the width of the trench; therefore it can be established according to technical-economic criteria, while the depth affects significantly the attenuation process.Some results of parametric study (W=w/λr,D=d/λr).
(a)(b)


## 4.1. Material Constitutive Model and Type of FEM Elements Applied

Linear viscoelastic constitutive model for the embankment and the soil medium materials has been adopted to carry out the investigation, according to some previous studies [20, 21]. Both the embankment and the soil medium have been meshed with triangular elements having 6 nodes. The mechanical parameters assigned to each element are  the Young modulus E,   Poisson’s ratio ν, and the mass density ρ, while the damping behaviour has been introduced by means of the Rayleigh coefficientsα and β. The damping matrix [C] has been defined on the basis of the quoted Rayleigh coefficients α and β.The presence in the analysed scheme of a rigid base has been taking into account by introducing in the model fixities at   bedrock.


## 4.2. Modelling of the Rail  Source

The railway source has been modelled applying in the surface of the embankment below the ballast a stress having a value equal to the dynamic stress produced in the contact surface ballast-embankment.In particular, according to [6], two concentrated linear loads acting over the rails have been considered. The time history of each load consists of four consecutive impulses; each impulse have a time duration of 0.02 seconds and an amplitude equal to 1000 KN. Since only the embankment has  been modelled, the dynamic pressure produced by the linear loads in the contact surface ballast/embankment, having dimensions of about 2,5m×1m, has been fixed in  400 kN/m2 (Figure 3) while its time history is the one plotted in Figure 4.Figure 3
Loading conditions.Figure 4
Time history of the rail force.The applied load function gives a significant contribution for the frequencies included in the range from 20 to 35 Hz, with a central frequency of about 27 Hz. This range, in the present study-case, is representative of the frequencies produced by the passage of a freight train at the velocity of 60 km/h, and having a boogie axle distance of 9 m. The load can appear overabundant but it takes into account the dynamic amplification factor by the  shortest wavelength [6, 22].


## 4.3. Finite Element Mesh

The definition of the model mesh is in general a compromise between the necessity to restrict the number of the elements and that to limit their dimensions. However, the element’s dimension must be accurately chosen in order to adequately represent the deformed shape associated to the wave lengthsλ at the different harmonic frequencies f imputed.Iffmax is the greatest frequency of the significant components of the Fourier spectrum and Vj is the propagation speed of the waves (P-waves, S-waves, or R-waves) in the material, the dimension of the mesh element for the material should satisfy the following relation:(1)d≤λmink=Vjk·fmax,
where k is a coefficient ranging from 4 to 10 according to the type of finite element and to its shape function [23].In the light of this, the maximum dimension of the element has been defined in the following way:(i)
on the basis of the material properties, there has been calculated the minimum wave propagation velocity for the considered domain; that is, for the superficial elements there has been considered the Rayleigh’s velocityVr, for the elements close to deep source the shear velocity Vs;(ii)
there has been defined the maximum value of the appropriate frequencyfmax; that is, in this study, according to the railway source model,  fmax has been set up equal to 35 Hz;(iii)
there has been calculated, using the expression (1), the minimal dimension of the element.By applying this procedure to the two types of soil medium, there has been obtainedClay⇒Vr=108.16 m/sec; fmax=35 Hz; k=6⇒dA≤0.51 m;Sand⇒Vr=270.94 m/sec; fmax=35 Hz; k=6⇒dS≤1.29 m.By considering these dimensions, the mesh has been made as follows:(i)
in the part of the model where the response is significant to the aim of the study, corresponding to a strip, large about 40 m and 40 m deep, closest to the rail track, a maximum dimension of 0.5 m has been assigned to the elements; this mesh has been designed to be denser in this focal area in order to achieve a good precision;(ii)
in the part of the model included between the quoted part and the lower edge a maximum dimension of 0.5 m has been assigned when the clayey soil has been analysed, while in the case of sandy soil the dimensions of the element gradually increase until 1.5 m near the lower edge;(iii)
externally to the above mentioned areas, the dimensions of the elements increase gradually reaching the value of 3 m in proximity of the vertical boundary far 300 m from the symmetry axis.The linear variation of the mesh dimensions, from the sidey=40 m to the boundary, has been chosen in accordance with other similar works [24, 25].The scheme of the finite element mesh is shown in Figure5.Figure 5
Scheme of FEM mesh.Figure 6
Rayleigh’s parameters determination.It can be observed that outside the zone of interest (40m×40m), the growth of the dimensions of the element  in connection with the time step does not produce effects in the part of the model where is sited the trench; in fact, the vibratory waves, also when undergo reflections in correspondence of the vertical boundary right hand, cannot reach this zone before the end of the analysis.


## 4.4. Courant’s Condition

The time step integration has been assigned taking into account the Courant condition, that defines the maximum time step as [26](2)Vj·Δth=C,C≤1
in which Vj is the velocity of the considered wave, h is the maximum spatial displacement admitted for the perturbation  that in FEM models is generally assumed equal to the dimension of the mesh, and C is the Courant number.It should be noted  that in the problems involving different type of waves, like the propagation of theP and S waves in the elastic space, the time step is imposed by the greatest propagation velocity, while the element’s dimensions depend on the lowest velocity. So, a great difference between the greatest and the lowest propagation velocity requires a small finite element and a short-time step.Having in mind this, in this task there has been assumed thatClay⇒Vp=191.55 m/sec, dA,min=0.5m⇒Δt≤0.0026 sec; there has been chosen Δt=0.002 sec.Sand⇒Vp=479.83 m/sec, dS,min=0.5m⇒Δt≤0.00104 sec; there has been  chosen Δt=0.001 sec.In both cases of clay and sand the FEM simulations have been performed in a time ofT=0.75 sec; this time is sufficient to permit the complete passage of the dynamic perturbation in the examined zone.


## 4.5. Boundary Condition

Boundary conditions have been imposed in the model in consideration of the geometric and mechanical conditions of symmetry and taking into account the need to limit the errors, namely, the difference existing between the theoretical value of the stress in a point of the boundary and the value obtained in presence of the restraints.FE calculations need a self-consistent simulation area to operate on. Especially in the case of wave propagation, special boundary conditions have to be incorporated into the calculation scheme [27]. Many authors have proposed special absorbing boundary conditions (ABCs) for the FE calculation of elastic waves. In these works the boundary conditionssuggested from Lysmer and Kuhlemeyer have been applied [28]:(3)σ=-ρcdVn,τ=-ρcsVt,
where σ and  τ are  the normal and the shear stress,  ρ  the material density, and  cd  and cs the longitudinal and the shear wave velocities of the transmitting media. These equations reveal that the magnitude of these stresses at the boundaries is proportional to the particle velocities in the normal (Vn) and in the tangential (Vt) directions.   Lysmer’s dampers placed on the artificial boundary are effective in reducing unwanted wave reflections if the boundary of the finite element mesh is sufficiently far outward.In the light of this, there has been applied(i)
in the symmetry edge, symmetry restraints like bi-pendulum iny direction;(ii)
in the lower edge, fixities at the bedrock;(iii)
in the right-side edge,   Lysmer’s dampers in the directionsy and z.The right-side edges have been applied at a distance of 300 m from the source; the errors obtained in the boundary have been(i)
ε(σn)=0.10%  and  ε(τ)=0.45%, in the case of the clayey soil medium;(ii)
ε(σn)=0.26% and  ε(τ)=1.12%, in the case of the sandy soil medium.


## 4.6. The Determination of the Rayleigh Damping Coefficientsα  and β

In an FEM model the damping behaviour can be modelled using   Rayleigh’s method. The damping matrix can be obtained from the following relation:(4)[C]=α[M]+β[K],
where [M] and [K] are, respectively, the mass matrix and the stiffness matrix.Theαand β damping coefficients [23, 29] have been determined following the procedure described below and shown in Figure 4.(i)
For each 2D model, a frequency analysis has been carried out and the first resonance frequencyf1 has been calculated.(ii)
There has been calculated the frequencyf2 by means of the relation f2=nf1, where n is the first odd number greater than the ratio fi/f1, being fi the dominant frequency. In this case fi has been fixed equal to 27 Hz, like the medium frequency of the load function.(iii)
Rayleigh’s parameters value has been calculated by establishing that, for frequencies equal tof1 and f2, the damping coefficient of the element must be equal to the value ξs, which is proper to the examined soil. In this way there has been introduced a subdamping for the frequencies ranging from f1 to f2 and a superdamping for the other frequencies not included in this range.For the clayey soil   Rayleigh’s parameters have been calculated assuming that a damping ratioξs   is equal to 0.06, while for the sandy one there has been assumed that  ξs  is  equal to 0.04.In   Figures7 and 8, we can see examples of damping diagrams for clayey and sandy soil medium.Figure 7
Damping ratioξs versus frequency for clay bank of 25 m depth.Figure 8
Damping ratioξs versus frequency for sand bank of 40 m depth.


## 4.7. Validation Model

In order to validate the model, results by Beskos et al. [11] have been considered and used. Results concern an open trench having a depth of 2.8 m (1λr) and a width of 0.28 m (0.1λr); it is made in an elastic half-space at a distance of 10.04 m (5λr) from the axis of a rigid foundation having a width of 0.7 m and subjected to a harmonic load with frequency of 50 Hz.In particular the comparison between the results given by the proposed 2D FEM model and the Beskos one has been performed considering the trend of the attenuation ratioAr versus the normalized distance  y/λr, where λr is the length of   Rayleigh’s waves generated from the source and y the horizontal distance of the generic point from the axis of the source.In each point of both vibration-sensitive site and site close to the trench the attenuation ratioAr  has been defined as the ratio between the maximum value of the kinematic parameters (y and z displacements and y and z velocities) in presence of trench and the same  value without trench:(5)Ar=(maxvalueofthekinematicparametersinapointwithtrench)/(maxvalueofthekinematicparametersinapointwithouttrench).
The attenuation ratio has been calculated not only for the vertical displacement (dir. z), as it is usual in   literature, but also for the horizontal displacement (dir y) and for the horizontal and vertical velocities.Since the sites are not   points but   large zones and considering the evolution of the attenuation ratio increasing distance (Figure7), the effectiveness of the trench in these sites   has been evaluated in terms of average attenuation ratio:(6)Am=1C∫ocAr(y)dy,
where C is the length of the investigated zone (equal to 10 m for vibration sensitive site and 5 m for site close to the trench).As it can see in Figure9, a good agreement between the numerical results from Beskos simulation and the ones obtained with the proposed model has been achieved. The model has given good results also in others experimental applications [29] and thus it has been applied in the present study.Figure 9
Model results versus Beskos results.


## 4.8. Results

From the numerical FEM simulations of the 540 scenarios the values of the displacements and velocities have been obtained; the minimum and the maximum values reached in the significant points are summarised in Table4.Table 4
Maximum and minimum values of displacements and velocities.
Site close to the trenchSoil  typeVelocity  dir.y min max  (m/sec)Velocity dir.z min max (m/sec)Displacement dir.y min max  (m)Displacement dir.z min max  (m)Sand0.130.090.000980.000560.350.250.002520.00619Clay0.160.130.005640.003690.290.480.010840.02509Site vibration-sensitiveSoil typeVelocity  dir.y min max (m/sec)Velocity dir.z min max (m/sec)Displacement dir.y min max  (m)Displacement dir.z  min  max  (m)Sand0.0520.0350.000780.000090.210.200.001980.00341Clay0.0350.0280.001770.000400.270.300.010840.01714In Figure10 two samples of time history of displacement and velocity in a point are plotted.Samples of results in a point in case of sandy soil.
(a)(b)When the horizontal propagation of in-plane waves is considered, the Rayleigh waves and theP-waves are the main causes of the ground motion.At quite distance (5λr) from a vibratory source the ground motion is due to the Rayleigh waves energy.If an open trench is excavated near the source, the incidence of the waves on this obstacle (discontinuity field) gives rise to reflected and transmitted body waves. Behind the open trench, as the distance increases, the transmittedP-waves get partially into Rayleigh waves (phenomenon known as mode conversion). Therefore, the Rayleigh waves are the main cause of  the ground motion [13].According to this phenomenon the results show a great attenuation of displacement and velocity passing from the site close to the trench to the vibration sensitive site (see Table4).To evaluate the effectiveness of the trench there has been adopted the very strict following criterion: a trench is effective in screening the ground borne vibration if the value of the attenuation ratio is lower than 0.4 for every considered displacement  and velocity.The analysis of the FEM results for the two types of ground has been performed having in mind this criterion.The results give some indications for the design of the open trenches in both cases of site close to the trench and vibration-sensitive site (see Tables5 and 6).Table 5
Design parameters for trench—site close to the trench.
d/hClayey soil layer[λr=4,00m]Sandy soil layer[λr=10,00m]dlwdlw1/5———0.2λr2λr≤0.1λr1/21.25λr>3λrall0.5λr≥1.2λrall2/3<2.75λr6λrall0.2λrallallTable 6
Design parameters for trench—vibration-sensitive site.
d/hClayey soil layer[λr=4,00m]Sandy soil layer[λr=10,00m]dlwdlw1/5———0.2λr0.8λr0.05λr1/21.25λr>3λrall0.5λr0.8λr0.05λr2/30.5λrallall0.2λr≥ 0.8λrallThe effect of frequency is fully taken into account by normalizing the open trench dimensions with respect to the Rayleigh wavelength.The obtained results are very strict in regard to the dimensions and the position of the trench because these features depend on the ratiod/h and so on the finite value of the layer thickness.For instance, in clayey soil results show that the trench is effective if it is 5 m deep and far away at least 12 m from the source, for thickness layer equal to 10 m or 2 m deep and far away at least 24 m from source, and  for layer thickness equal to 3 m.Some results of the parametric studies carried out are plotted in Figure11. It can be seen that the efficiency of the trench is in general independent from the width of the trench; therefore it can be established according to technical-economic criteria, while the depth affects significantly the attenuation process.Some results of parametric study (W=w/λr,D=d/λr).
(a)(b)


## 5. ANNs Analysis

### 5.1. Introduction

Artificial neural networks (ANNs) are a tool that simulates the biological processes and, as demonstrated by the wide applications in engineering field, they are able to solve functional mapping problems. ANNs are an assemblage of mathematical simple computational elements called neurons.In particular the back propagation neural network (BPNN) is a collection of neurons distributed over an input layer, that contain the input variables of the problem, one or more hidden layers with a certain number of nodes, and an output layer with a number of nodes equal to the output variables. The nodes between layers are connected by the links having a weight that describes quantitatively the strength of the connection.A significant effort is required in the selection of the ANN architecture, particularly, as it is obvious, in the definition of the hidden layers and the corresponding nodes.The learning process utilized for this type of ANN is the “back-propagation learning” which consists in an error minimization technique [30]. Inputs from the mapping examples are propagated through each layer and give the outputs that are compared with the correct answers; the difference represents the error. This error is propagated backward through the network and the connection weights are individually adjusted so as to reduce and minimize the error.To ensure an efficient convergence and the desired performance of the trained network, several parameters are incorporated in the training phase. These parameters include the learning rate, the momentum term, and the number of training iterations. The learning rate is a factor that proportions the amount of the adjustment applied at each time the weight is updated. The use of a momentum term could carry the weight change process through one or more local minima and get it into global minima. The early stopping method and the number of training epochs determine the training stop criteria [31].


### 5.2. Neural Network Application

For the present application the neural networks have been trained and tested with the data given from the FEM modelling of the 540 scenarios.In details, the input data chosen for the neural analysis have been 6 and precisely: width of the trench (w), depth of the trench (d), ratio of depth of the trench to thickness of the soil layer (d/h), shear wave velocity (Vs), and damping ratio (ξ), while the output data have been the average attenuation ratio of the horizontal and vertical displacements, and the average attenuation ratio of the horizontal and the vertical velocity. In fact, the neural network application aimed at evaluating the contribution of each input factor in the attenuation of the displacements and velocities.To this purpose a great number of trial ANNs with one, two, and three hidden layer have been trained to evaluate the performance of different network architectures in the comprehension and generalisation of the problem.The training procedure has been tested using a test set of examples having a percentage of 10% of the training set; the prediction performance of the developed model has been exhibited in the test set.The error term, represented by the mean square error (RMS error), computed for the validation input-output pairs, has been monitored during the training process of the networks. The error normally decreases during the initial phase of training. However, when the network begins to overfit the data (a situation arising when ANN works well only with the training data), the error on the validation set will typically begin to increase. When the validation error increased for a specified number of iterations, training was stopped and the weights at the minimum of the validation error were saved. This point is the point of maximum generalization [30].Considering this, networks with errors comprising in the following ranges have been selected: RMS error training set: 0.02–0.1; RMS error test set: 0.02–0.12; Correlation training set: 0.5–1. It has been stated that the neural networks having the training parameters included in the quoted ranges provide acceptable results.The analysis of the training and the testing phases of these neural networks allowed to evaluate the percentage of contribution of each input factor in the output determination. In other terms the study allowed to assess qualitatively and quantitatively the influence of the features of the trench and the soil medium in the attenuation process of ground-borne vibrations.


### 5.3. Results

The results obtained for both vibration-sensitive site and site close to the trench are explained in the following.


#### 5.3.1. Site Close to the Trench

Regarding the attenuation ratio of the horizontal and vertical displacements, the neural network analysis has given the percentages of contribution summarised in Tables7 and 8.Table 7
Percentage of contribution to the attenuation of the displacements—site close to the trench.
Attenuation ratio for horizontal displacementsAttenuation ratio for vertical displacementsParametersPercentage of contributionParametersPercentage of contributiond28.00%d31.00%d/h26.00%d/h28.00%Vs17.00%l14.50%ξ15.50%Vs13.00%l10.50%ξ9.00%w3.00%w4.50%Table 8
Percentage of contribution to the attenuation of the velocities—site close to the trench.
Attenuation ratio for horizontal velocitiesAttenuation ratio for vertical velocitiesParametersPercentage of contributionParametersPercentage of contributiond31.00%d27.00%d/h23.00%d/h26.00%ξ17.00%l15.00%l14.00%Vs13.50%Vs11.00%ξ11.50%w4.00%w7.00%By analysing the results it can be deduced that for both vertical and horizontal displacements the percentage of contribution to the attenuation ratio of the depth of the trenchd and of the ratio d/h is close  to the 60%. The influence of the terrain, expressed by  Vs and ξ, ranges from 32% for horizontal displacements to about 22% for vertical displacements. The influence of the distance of the trench from rail l is evaluated around the 10%–14% while the width of trench whas a marginal weight of about 3–5%.Regarding the attenuation of the velocities the hierarchy of influence of the parameters seems to be the same of the displacements: the influence of the depth of the trenchd and of the ratio d/h is preponderant (about 50%), while the influence of the terrain is evaluated around the 25%–28% while, in contrast to the results found for the displacements, there appears significant the distance of the trench from rail l that contributes for around the 15% to the attenuation process. The width of the trench w, even if with the different weights  for horizontal (4%) and vertical (7%) velocities, gives a poor contribution to the attenuation phenomenon.


#### 5.3.2. Vibration-Sensitive Site

The percentages of contribution to the attenuation ratio of the horizontal and vertical displacements are summarised in Tables9 and 10.Table 9
Percentage of contribution to the attenuation of the displacements—vibration-sensitive site.
Attenuation ratio for horizontal displacementsAttenuation ratio for vertical displacementsParametersPercentage of contributionParametersPercentage of contributiond36.00%d42.00%d/h30.00%d/h26.00%Vs13.00%Vs12.00%ξ11.00%ξ11.00%l7.00%l7.00%w3.00%w2.00%Table 10
Percentage of contribution to the attenuation of the velocities—vibration-sensitive site.
Attenuation ratio for horizontal velocitiesAttenuation ratio for vertical velocitiesParametersPercentage of contributionParametersPercentage of contributiond34.00%d35.00%d/h27.00%d/h22.50%Vs17.00%Vs16.00%ξ14.00%ξ15.50%l5.00%l8.00%w3.00%w3.00%In this site the influence of the depth of the trenchd on the attenuation ratio of the vertical and horizontal displacements is very significant (36%–42%) and together with the ratio d/h seems to govern preponderantly the attenuation process. The other factors (Vs, ξ, l) have an influence comparable (each ranged from 7 to 13%) and lesser than d and d/h. Once again the width w of the trench is not much significant.Regarding the attenuation ratio of the horizontal and vertical velocities the trend is the same of the displacements but the contribution of the terrain’s parameters (Vs and ξ) is higher (about 30%). The distance of the vibration sensitive site from the source l and even more the width w of the trench have a low contribution.


## 5.1. Introduction

Artificial neural networks (ANNs) are a tool that simulates the biological processes and, as demonstrated by the wide applications in engineering field, they are able to solve functional mapping problems. ANNs are an assemblage of mathematical simple computational elements called neurons.In particular the back propagation neural network (BPNN) is a collection of neurons distributed over an input layer, that contain the input variables of the problem, one or more hidden layers with a certain number of nodes, and an output layer with a number of nodes equal to the output variables. The nodes between layers are connected by the links having a weight that describes quantitatively the strength of the connection.A significant effort is required in the selection of the ANN architecture, particularly, as it is obvious, in the definition of the hidden layers and the corresponding nodes.The learning process utilized for this type of ANN is the “back-propagation learning” which consists in an error minimization technique [30]. Inputs from the mapping examples are propagated through each layer and give the outputs that are compared with the correct answers; the difference represents the error. This error is propagated backward through the network and the connection weights are individually adjusted so as to reduce and minimize the error.To ensure an efficient convergence and the desired performance of the trained network, several parameters are incorporated in the training phase. These parameters include the learning rate, the momentum term, and the number of training iterations. The learning rate is a factor that proportions the amount of the adjustment applied at each time the weight is updated. The use of a momentum term could carry the weight change process through one or more local minima and get it into global minima. The early stopping method and the number of training epochs determine the training stop criteria [31].


## 5.2. Neural Network Application

For the present application the neural networks have been trained and tested with the data given from the FEM modelling of the 540 scenarios.In details, the input data chosen for the neural analysis have been 6 and precisely: width of the trench (w), depth of the trench (d), ratio of depth of the trench to thickness of the soil layer (d/h), shear wave velocity (Vs), and damping ratio (ξ), while the output data have been the average attenuation ratio of the horizontal and vertical displacements, and the average attenuation ratio of the horizontal and the vertical velocity. In fact, the neural network application aimed at evaluating the contribution of each input factor in the attenuation of the displacements and velocities.To this purpose a great number of trial ANNs with one, two, and three hidden layer have been trained to evaluate the performance of different network architectures in the comprehension and generalisation of the problem.The training procedure has been tested using a test set of examples having a percentage of 10% of the training set; the prediction performance of the developed model has been exhibited in the test set.The error term, represented by the mean square error (RMS error), computed for the validation input-output pairs, has been monitored during the training process of the networks. The error normally decreases during the initial phase of training. However, when the network begins to overfit the data (a situation arising when ANN works well only with the training data), the error on the validation set will typically begin to increase. When the validation error increased for a specified number of iterations, training was stopped and the weights at the minimum of the validation error were saved. This point is the point of maximum generalization [30].Considering this, networks with errors comprising in the following ranges have been selected: RMS error training set: 0.02–0.1; RMS error test set: 0.02–0.12; Correlation training set: 0.5–1. It has been stated that the neural networks having the training parameters included in the quoted ranges provide acceptable results.The analysis of the training and the testing phases of these neural networks allowed to evaluate the percentage of contribution of each input factor in the output determination. In other terms the study allowed to assess qualitatively and quantitatively the influence of the features of the trench and the soil medium in the attenuation process of ground-borne vibrations.


## 5.3. Results

The results obtained for both vibration-sensitive site and site close to the trench are explained in the following.


### 5.3.1. Site Close to the Trench

Regarding the attenuation ratio of the horizontal and vertical displacements, the neural network analysis has given the percentages of contribution summarised in Tables7 and 8.Table 7
Percentage of contribution to the attenuation of the displacements—site close to the trench.
Attenuation ratio for horizontal displacementsAttenuation ratio for vertical displacementsParametersPercentage of contributionParametersPercentage of contributiond28.00%d31.00%d/h26.00%d/h28.00%Vs17.00%l14.50%ξ15.50%Vs13.00%l10.50%ξ9.00%w3.00%w4.50%Table 8
Percentage of contribution to the attenuation of the velocities—site close to the trench.
Attenuation ratio for horizontal velocitiesAttenuation ratio for vertical velocitiesParametersPercentage of contributionParametersPercentage of contributiond31.00%d27.00%d/h23.00%d/h26.00%ξ17.00%l15.00%l14.00%Vs13.50%Vs11.00%ξ11.50%w4.00%w7.00%By analysing the results it can be deduced that for both vertical and horizontal displacements the percentage of contribution to the attenuation ratio of the depth of the trenchd and of the ratio d/h is close  to the 60%. The influence of the terrain, expressed by  Vs and ξ, ranges from 32% for horizontal displacements to about 22% for vertical displacements. The influence of the distance of the trench from rail l is evaluated around the 10%–14% while the width of trench whas a marginal weight of about 3–5%.Regarding the attenuation of the velocities the hierarchy of influence of the parameters seems to be the same of the displacements: the influence of the depth of the trenchd and of the ratio d/h is preponderant (about 50%), while the influence of the terrain is evaluated around the 25%–28% while, in contrast to the results found for the displacements, there appears significant the distance of the trench from rail l that contributes for around the 15% to the attenuation process. The width of the trench w, even if with the different weights  for horizontal (4%) and vertical (7%) velocities, gives a poor contribution to the attenuation phenomenon.


### 5.3.2. Vibration-Sensitive Site

The percentages of contribution to the attenuation ratio of the horizontal and vertical displacements are summarised in Tables9 and 10.Table 9
Percentage of contribution to the attenuation of the displacements—vibration-sensitive site.
Attenuation ratio for horizontal displacementsAttenuation ratio for vertical displacementsParametersPercentage of contributionParametersPercentage of contributiond36.00%d42.00%d/h30.00%d/h26.00%Vs13.00%Vs12.00%ξ11.00%ξ11.00%l7.00%l7.00%w3.00%w2.00%Table 10
Percentage of contribution to the attenuation of the velocities—vibration-sensitive site.
Attenuation ratio for horizontal velocitiesAttenuation ratio for vertical velocitiesParametersPercentage of contributionParametersPercentage of contributiond34.00%d35.00%d/h27.00%d/h22.50%Vs17.00%Vs16.00%ξ14.00%ξ15.50%l5.00%l8.00%w3.00%w3.00%In this site the influence of the depth of the trenchd on the attenuation ratio of the vertical and horizontal displacements is very significant (36%–42%) and together with the ratio d/h seems to govern preponderantly the attenuation process. The other factors (Vs, ξ, l) have an influence comparable (each ranged from 7 to 13%) and lesser than d and d/h. Once again the width w of the trench is not much significant.Regarding the attenuation ratio of the horizontal and vertical velocities the trend is the same of the displacements but the contribution of the terrain’s parameters (Vs and ξ) is higher (about 30%). The distance of the vibration sensitive site from the source l and even more the width w of the trench have a low contribution.


## 5.3.1. Site Close to the Trench

Regarding the attenuation ratio of the horizontal and vertical displacements, the neural network analysis has given the percentages of contribution summarised in Tables7 and 8.Table 7
Percentage of contribution to the attenuation of the displacements—site close to the trench.
Attenuation ratio for horizontal displacementsAttenuation ratio for vertical displacementsParametersPercentage of contributionParametersPercentage of contributiond28.00%d31.00%d/h26.00%d/h28.00%Vs17.00%l14.50%ξ15.50%Vs13.00%l10.50%ξ9.00%w3.00%w4.50%Table 8
Percentage of contribution to the attenuation of the velocities—site close to the trench.
Attenuation ratio for horizontal velocitiesAttenuation ratio for vertical velocitiesParametersPercentage of contributionParametersPercentage of contributiond31.00%d27.00%d/h23.00%d/h26.00%ξ17.00%l15.00%l14.00%Vs13.50%Vs11.00%ξ11.50%w4.00%w7.00%By analysing the results it can be deduced that for both vertical and horizontal displacements the percentage of contribution to the attenuation ratio of the depth of the trenchd and of the ratio d/h is close  to the 60%. The influence of the terrain, expressed by  Vs and ξ, ranges from 32% for horizontal displacements to about 22% for vertical displacements. The influence of the distance of the trench from rail l is evaluated around the 10%–14% while the width of trench whas a marginal weight of about 3–5%.Regarding the attenuation of the velocities the hierarchy of influence of the parameters seems to be the same of the displacements: the influence of the depth of the trenchd and of the ratio d/h is preponderant (about 50%), while the influence of the terrain is evaluated around the 25%–28% while, in contrast to the results found for the displacements, there appears significant the distance of the trench from rail l that contributes for around the 15% to the attenuation process. The width of the trench w, even if with the different weights  for horizontal (4%) and vertical (7%) velocities, gives a poor contribution to the attenuation phenomenon.


## 5.3.2. Vibration-Sensitive Site

The percentages of contribution to the attenuation ratio of the horizontal and vertical displacements are summarised in Tables9 and 10.Table 9
Percentage of contribution to the attenuation of the displacements—vibration-sensitive site.
Attenuation ratio for horizontal displacementsAttenuation ratio for vertical displacementsParametersPercentage of contributionParametersPercentage of contributiond36.00%d42.00%d/h30.00%d/h26.00%Vs13.00%Vs12.00%ξ11.00%ξ11.00%l7.00%l7.00%w3.00%w2.00%Table 10
Percentage of contribution to the attenuation of the velocities—vibration-sensitive site.
Attenuation ratio for horizontal velocitiesAttenuation ratio for vertical velocitiesParametersPercentage of contributionParametersPercentage of contributiond34.00%d35.00%d/h27.00%d/h22.50%Vs17.00%Vs16.00%ξ14.00%ξ15.50%l5.00%l8.00%w3.00%w3.00%In this site the influence of the depth of the trenchd on the attenuation ratio of the vertical and horizontal displacements is very significant (36%–42%) and together with the ratio d/h seems to govern preponderantly the attenuation process. The other factors (Vs, ξ, l) have an influence comparable (each ranged from 7 to 13%) and lesser than d and d/h. Once again the width w of the trench is not much significant.Regarding the attenuation ratio of the horizontal and vertical velocities the trend is the same of the displacements but the contribution of the terrain’s parameters (Vs and ξ) is higher (about 30%). The distance of the vibration sensitive site from the source l and even more the width w of the trench have a low contribution.


## 6. Conclusions

In order to intercept the elastic waves and in particular   Rayleigh’s ones generated by the moving source from reaching buildings between the railways and the sensitive site, vibration screening may be established by means of the open trenches. Several factors determine the effectiveness of the open trenches in the screening of the ground borne vibrations in soil medium.In the present work, as outcomes of an FEM modelling of the problem and of an analysis carried out by means of the neural network, the weights of the main geometric features of both the trench and the soil medium have been established.The main conclusions of this application are the following:(i)
first of all, the screening performance of the trench is mainly affected by the depth of the trench, in both cases of site close to the trench and vibration-sensitive site. From ANN analysis the percentage of contribution of this geometric parameter in reducing displacements and velocities of vibration reached the 25%–35%;(ii)
layer thickness increasing, the attenuation ratio increases according to the depth of the trench; to this reason, it appears very noteworthy in the topic not only the depthd of the trench but also the ratio d/h (depth of the trench/ thickness of the bank); the two parameters d and d/h govern the attenuation process with a percentage of contribution more than 50%;(iii)
to the aim of evaluating the effectiveness of the trench we have to analyse all the kinematic parameters and not only the displacement in one direction, generally the directionz; in fact, some geometric configurations able to reduce one kinematic parameters  can be inefficient with regard to the attenuation of the other kinematic parameters.



---

*Source: 942787-2010-01-13.xml*