Effect of periodic vibration-isolation piles on soil vibration in deep-buried tunnels of suburban railways

Abstract

This study aims to investigate the impact of isolation piles on soil vibrations in the environment surrounding suburban railways. Initially, a comprehensive numerical model of the train was established to simulate the wheel-rail interaction forces, which were then applied to a three-dimensional coupled track-soil model. The accuracy of the model was validated through comparison with measured data. The focus of the research is to analyze the vibration isolation effects of single-row, double-row, and triple-row piles, with particular attention to the influence of pile spacing in double-row piles on isolation effectiveness. The study found that the effectiveness of isolation piles in reducing soil vibration acceleration within a depth range of 0 m to 2 m decreases as the pile spacing increases. When the pile spacing reaches 2 times the pile radius, single-row piles exhibit the best vibration isolation effect. For reducing peak soil acceleration within a depth range of 4 m to 8 m, single-row piles are most effective when the pile spacing exceeds 1 times the pile radius. To attenuate peak soil acceleration at depths greater than 6 m, the pile spacing for single-row and double-row piles should be set at 0.5 times the pile radius, while triple-row piles should have a spacing of 1 times the pile radius. In terms of reducing peak soil velocity, single-row piles demonstrate better vibration isolation effects. To reduce peak velocity in soil at depths greater than 2 m, the pile spacing for single-row and double-row piles should be set at 0.5 times the pile radius, while triple-row piles should be spaced at 1 times the pile radius. For reducing peak soil displacement, single-row piles should be used, with optimal results achieved at a spacing of 1 times the pile radius.

Keywords

Vibration isolation piles peak acceleration peak velocity numerical model pile spacing

Introduction

The vibrations induced by rail transit significantly disturb human life and work and, in severe cases, can even lead to the destruction of buildings.^1–4 Typically, these issues are addressed through active and passive vibration isolation methods^5,6 Active vibration isolation systems are commonly employed to isolate vibration sources from the ground,^7,8 while passive vibration isolation systems are frequently used to reduce the energy transmitted to structures.⁹ This approach avoids altering the dynamic characteristics of the vibration source and the structure, thereby preventing potential safety issues at these locations. Passive vibration isolation methods can be broadly categorized into pile isolation and trench isolation techniques.^10–12 Extensive theoretical research and engineering practice have been conducted on pile isolation methods, revealing that isolation piles can effectively isolate environmental vibrations^13,14 and reduce the energy of vibrations transmitted to structures.^15,16 Consequently, isolation piles are widely applied in vibration mitigation measures.¹⁷

Various methods are currently employed to study the vibration isolation effects of row piles¹⁸ including field experiments, model tests, and numerical simulation methods. Field experiment results are authentic and reliable,^19,20 yet they are challenging to coordinate, are costly, and involve significant environmental disturbances.²¹ Model tests, based on similarity principles, employ scaled models and yield relatively authentic and reliable results. However, these systems require substantial floor space and entail high costs for construction, operation, and maintenance.^22,23 Numerical simulation, which uses computers to approximate model solutions, requires validation through field experiments or model tests. Nonetheless, numerical simulation presents no coordination difficulties, has no land requirements, and offers adjustable and controllable testing costs,²⁴ making it a widely adopted research method.^25–29 Zhang³⁰ employed a finite element model to analyze the vibration reduction performance of periodic piles under subway ground source impact. Using the measured vibration source intensity as the excitation input, the study investigated the influence of different burial depths and row numbers on the damping effect of the piles. Wang et al.³¹ established a coupled model of moving vehicle loads and subsoil construction structures, finding that the ground vibration in sensitive areas induced by moving vehicle loads is primarily vertical and low-frequency, with the main frequency range being 3 to 4 Hz, consistent with ground pulsation tests. Zhang et al.³² utilized the finite element analysis software ANSYS to analyze the vertical ground vibration amplitude behind a single row of cast-in-place concrete thin-walled pipe piles and investigated its influencing factors.

However, previous studies have mainly focused on the attenuation of surface vibrations, with less attention given to their impact on underground soil vibrations.^33–37 With the rapid urbanization in China and the extensive construction of underground structures (such as basements, underground parking lots, and underground shopping malls), studying soil vibrations has become increasingly important. Therefore, this paper combines numerical analysis and finite element methods to establish a coupled train-tunnel-soil model. The study analyzes the vibration isolation effects of isolation piles on peak acceleration, velocity, and displacement in the time domain under various influencing factors and explores the impact of pile spacing at different pile radii on soil vibrations (as shown in the process in Figure 1).

Figure 1.

Technical roadmap.

Model validation

The model establishment in this study involves three steps. First, a numerical model of the train was developed using Intel Visual Fortran to calculate the wheel-rail interaction forces.³⁸ Second, a finite element model of the rail-tunnel-soil coupling system was established. Third, the calculated wheel-rail interaction forces were applied to the finite element model.

Wheel-rail interaction force modeling

The train model adopts a two-degree-of-freedom suspension train model, where the vertical and pitch degrees of freedom are considered for both the car body and the bogie, while only the vertical degree of freedom is considered for the wheels. The primary suspension consists of a spring k_1s and a damper c₁ in series, which is then in parallel with another spring k_1p; the secondary suspension consists of a spring k_2p and a damper c₂ in parallel. Therefore, the complex stiffness of the suspension is expressed as:^39,40

k_{1} = \frac{k_{1 p} k_{1 s} + i Ω c_{1} (k_{1 p} + k_{1 s})}{k_{1 s} + i Ω c_{1}}

(1)

k_{2} = k_{2 p} + i Ω c_{2}

(2)

Under the action of the vertical wheel-rail force, the vertical motion equation of the vehicle body can be expressed as:

{\begin{matrix} M_{c} {\ddot{z}}_{c} - k_{2} z_{B 1} - k_{2} z_{B 2} = 0 \\ J_{c} {\ddot{θ}}_{c} + 2 k_{2} l_{B}^{2} θ_{c} - k_{2} l_{B} z_{B 1} + k_{2} l_{B} z_{B 2} = 0 \end{matrix}

(3)

where M_c is the mass of the car body, z_c is the vertical displacement of the car body, θ_c is the rotation angle of the car body, J_c is the moment of inertia of the car body, z_B1 and z_B2 are the vertical displacements of the front and rear bogies respectively, c is the yaw angle of the car body, and l_B is half the wheelbase. The motion equation for the front bogie is as follows:

{\begin{matrix} M_{B 1} {\ddot{z}}_{B 1} + (k_{2} + 2 k_{1}) z_{B 1} - k_{2} z_{c} - k_{2} l_{B} θ_{c} - k_{1} z_{w 1} - k_{1} z_{w 2} = 0 \\ J_{B 1} {\ddot{θ}}_{B 1} + 2 k_{1} l_{w}^{2} θ_{B 1} - k_{1} l_{w} z_{w 1} + k_{1} l_{w} z_{w 2} = 0 \end{matrix}

(4)

where z_w1 and z_w2 represent the vertical displacements of the first and second axles of the train, M_B1 is the mass of the front bogie, θ_B1 is the yaw angle of the front bogie, J_B1 is the moment of inertia of the front bogie, and l_w is half of the fixed axle distance of the bogie. The motion equation for the rear bogie of the train is given by:

{\begin{matrix} M_{B 2} {\ddot{z}}_{B 2} + (k_{2} + 2 k_{1}) z_{B 2} - k_{2} z_{c} + k_{2} l_{B} θ_{c} - k_{1} z_{w 3} - k_{1} z_{w 4} = 0 \\ J_{B 2} {\ddot{θ}}_{B 2} + 2 k_{1} l_{w}^{2} θ_{B 2} - k_{1} l_{w} z_{w 3} + k_{1} l_{w} z_{w 4} = 0 \end{matrix}

(5)

where z_w3 and z_w4 represent the vertical displacements of the third and fourth axles of the train, M_B2 is the mass of the rear bogie, θ_B2 is the yaw angle of the rear bogie, and J_B2 is the moment of inertia of the rear bogie. The vertical motion equations for the four axles are given by:

{\begin{matrix} M_{w 1} {\ddot{z}}_{w 1} + k_{1} z_{w 1} - k_{1} z_{B 1} - k_{1} l_{w} θ_{B 1} = - p_{1} \\ M_{w 2} {\ddot{z}}_{w 2} + k_{1} z_{w 2} - k_{1} z_{B 1} + k_{1} l_{w} θ_{B 1} = - p_{2} \\ M_{w 3} {\ddot{z}}_{w 3} + k_{1} z_{w 3} - k_{1} z_{B 2} - k_{1} l_{w} θ_{B 2} = - p_{3} \\ M_{w 4} {\ddot{z}}_{w 4} + k_{1} z_{w 4} - k_{1} z_{B 2} + k_{1} l_{w} θ_{B 2} = - p_{4} \end{matrix}

(6)

The motion equation for the steel rail is given by:

E I \frac{\partial^{4} w (x, t)}{\partial x^{4}} + m_{R} \frac{\partial^{2} w (x, t)}{\partial t^{2}} + k_{p} [w (x, t) - w_{2} (x, t)] = F_{0} (x, t)

(7)

where x represents the coordinate in the longitudinal direction of the rail; w represents the vertical displacement of the steel rail; and w₂ is the vertical displacement of the sleeper layer. m_R is the mass per unit length (linear density); k_p represents the stiffness per unit length of the rail pads, which simulate the trackbed layer; F₀ represents any load acting on the steel rail; EI represents the flexural stiffness of the cross-section.

The sleeper layer is still manufactured using the Euler beam model along the x-direction but without bending stiffness. Therefore, the motion equation is given by:

m_{s} \frac{\partial^{2} w_{2} (x, t)}{\partial t^{2}} + k_{p} [w_{2} (x, t) - w (x, t)] = - F_{1} (x, t)

(8)

where F1 represents the interaction force between the sleeper and the dam, and ms represents the sleeper's linear density when modeled as a flat beam with no bending stiffness extending along the track direction. The motion equation for the uniform mass approximation of the embankment layer is given by:

\frac{m_{B}}{6} [\begin{matrix} 2 & 1 \\ 1 & 2 \end{matrix}] {\begin{matrix} \frac{\partial^{2} w_{2} (x, t)}{\partial t^{2}} \\ \frac{\partial^{2} w_{3} (x, t)}{\partial t^{2}} \end{matrix}} + k_{B} [\begin{matrix} 1 & - 1 \\ - 1 & 1 \end{matrix}] {\begin{matrix} w_{2} (x, t) \\ w_{3} (x, t) \end{matrix}} = {\begin{matrix} F_{1} (x, t) \\ - F_{2} (x, t) \end{matrix}}

(9)

F₂ represents the interaction force between the dam and the foundation.

In equations (7)–(9), EI, k_B, and k_p are complex numbers, with their imaginary parts representing the respective damping coefficients.

The vertical vibration of the wheel generates a series of wheel-rail interaction forces that vary harmonically with time at the wheel-rail contact point. These forces can be represented as

P (t) = \tilde{P} (Ω) e^{i Ω t}

(10)

where Ω represents the vertical vibration frequency of the train,

i = \sqrt{- 1}

The representation of the track displacement amplitude under the action of resonant wheel-rail forces P(t) at the wheel-rail contact point is given by:

{\tilde{Z}}_{R} (Ω) = H_{R} (Ω) \tilde{P} (Ω)

(11)

where

H_{R} (Ω)

is the transfer function matrix of the track-ground system at the wheel-rail contact point.

During regular train operation, the wheel and rail remain in contact. By considering the continuity of vertical displacements between the wheel and rail, the wheel/rail interaction forces can be solved. The vertical displacement of the wheel axle is composed of the Hertz spring compression, wheel/rail irregularities, and vertical deformation of the rail. It can be expressed as follows:

{\tilde{Z}}_{w} (Ω) = {\tilde{Z}}_{R} (Ω) + \tilde{Z} (Ω) + K_{h}^{- 1} \tilde{P} (Ω)

(12)

The vector

\tilde{Z} (Ω) = {e^{i Ω a_{1} / c}, e^{i Ω a_{2} / c}, \dots, e^{i Ω a_{M} / c}}^{T}

represents the vertical harmonic roughness amplitudes at all M wheel-rail contact points of the train, and the diagonal matrix

K_{h} = diag (k_{h 1}, k_{h 2}, \dots, k_{h M})

represents the Hertz contact stiffness at all wheel-rail contact points.

Substituting equation (11) and z_W1/ z_W2/ z_W3/z_W4 into equation (12), the expression for the wheel-rail interaction force can be obtained as follows.³⁸

\tilde{P} (Ω) = - [H_{R} (Ω) + H_{v} (Ω) + K_{h}^{- 1}]^{- 1} \tilde{Z} (Ω)

(13)

Finite element model establishment

The measurement site was selected near Kaifa Avenue in Taizhou, Zhejiang Province, China (as shown in Figure 2), with data collected on the morning of April 28, 2023. The train was operating at a speed of 16.9 m/s.

Figure 2.

Schematic diagram of on-site measurements. (a) Map of the observation site. (b) Observation instruments.

The selected tunnel cross-section includes seven layers of soil. In the field data collection, six measurement points were arranged, with each point spaced 10 m apart. Additionally, the tunnel train is 93 m long, with a tunnel burial depth of 23 m, an outer diameter of 8.5 m, and an inner diameter of 7.7 m. The tunnel depth ratio is approximately 3d, as shown in Figure 3. Consequently, the finite element model size was chosen to be 200 m long, 100 m wide, and 30 m deep, as depicted in Figure 4.

Figure 3.

Schematic diagram of tunnel cross-section structure.

Figure 4.

Schematic diagram of finite element model.

Considering that the vibrations induced by the train result in small-strain deformations that do not produce plastic strain, an elastic model was adopted. To reduce computational costs, a linear elastic constitutive model was used for the soil. The rails were modeled using beam elements, and the track slabs were modeled using plate elements. The model was meshed using 10-node tetrahedral elements, comprising a total of 123,631 elements.

In reality, soil can extend infinitely. However, when applying finite element methods, a specific range of soil is artificially selected, introducing artificial boundaries on the front, back, sides, and bottom of the model. These boundaries can cause reflections of vibrations, which do not align with actual conditions. To minimize wave reflections near the boundaries, viscoelastic boundary conditions were applied at the model boundaries. All surfaces, except the ground surface, are fully constrained, and viscoelastic artificial boundaries are used to ensure the absorption of stress waves.

The contact relationships between the soil layers are set as rigid connections. The interfaces between the soil layers are continuous, meaning that stress and strain are continuously transmitted between the layers without significant separation or slippage. The calculated wheel-rail forces were applied as concentrated external loads to the finite element model of the track. This was achieved by gradually moving the position of the wheel-rail forces at different time steps. The Newmark-β method was used in solving the model. This method strikes a balance between stability and accuracy, making it suitable for most dynamic problems. The soil layer parameters were derived from geological survey data at the measurement site, as shown in Table 1.

Table 1.

Calculation parameters of the soil layer in the suburban railway underground section.

Soil layer	Soil thickness (m)	Shear wave (m·s⁻¹)	Densities (kN·m⁻³)	Poisson's ratio	Rayleigh α	Rayleigh β
Artificial fill	0.7	41.62	19	0.49	2.128	4.1e-5
Clay(1)	3.2	100.3	13.75	0.49	2.128	4.1e-5
Silt	8.1	77.97	11.1	0.49	2.128	4.1e-5
Clay(2)	1	116.3	15.09	0.49	2.128	4.1e-5
Clay(3)	1.7	163.6	13.65	0.49	2.128	4.1e-5
Silty clay	7.4	198.9	15.48	0.49	2.128	4.1e-5
Weathered tuff	7.9	309.9	24	0.49	2.128	4.1e-5

Model validation

To facilitate model validation, measurement points were set up both on the ground surface and inside the tunnel in the underground section of the suburban railway for data collection. The spatial relationships between the measurement points are shown in Figures 3 and 5, with a spacing of 30 m between the centerlines of the two tunnels. Data were collected using Donghua 2D001 accelerometers, which have a sensitivity of 0.3 V·s²/m and a maximum range of 20 m/s². The data acquisition equipment used was the Donghua DH5922D dynamic signal testing and analysis system.

Figure 5.

Schematic of measurement points on the site.

Since impulse loads have a broad frequency range, they can more comprehensively verify the model's reliability compared to simple harmonic loads. Therefore, a unit impulse load was applied to the center of the model's top surface, and the acceleration impulse response at the application point was obtained (as shown in Figure 6). Figure 6 shows that in the latter part of the response, the acceleration essentially returns to zero, indicating that the model has good reliability.

Figure 6.

Acceleration impulse response.

The wheel-rail interaction forces were calculated based on the vehicle parameters in Table 2 and the track structure parameters in Table 3, and these forces were applied to the track in the finite element model. The maximum vertical vibration levels obtained from the model calculations were then compared with the maximum vertical vibration levels from the measured data to validate the model's reliability, as shown in Figure 7 and Table 4.

Figure 7.

Comparison of maximum Z-vibration levels.

Table 2.

Parameter values for the vehicles.³⁹

Vehicle parameters	Unit	Numerical value
Mass of car body	kg	5.0 × 10⁴
Pitch inertia moment of car body	kg·m²	1.5 × 10⁶
Mass of bogie	kg	2.98 × 10³
Pitch inertia moment of bogie	kg·m²	3.605 × 10³
Wheelset mass	kg	1.35 × 10³
Primary vertical stiffness per axle	kN/m	2500
Secondary vertical stiffness per axle	kN/m	1000
Primary vertical viscous damping	kN·s/m	120
Secondary vertical viscous damping	kN·s/m	110
Primary viscous damping	kN·s/m	9.8
Secondary viscous damping	kN·s/m	6.0
Distance between bogies	M	13
Distance between axles under a bogie	M	2.6
Wheel radius	M	0.42
Hertzian stiffness	kN/m	2.4 × 10⁶

Table 3.

Parameter values for the track structure.³⁹

Rail components	Symbols	Parameter	Value
Rails	m_R	Mass per unit length	103 kg/m
	E₀I	Bending stiffness	1.32 × 10⁴ kN·m²
	η_R	Hysteric damping ratio	0.005
Rail pad	k_P0	Undamped stiffness	3.5 × 10⁵ kN/m²
Rail pad	η_P	Hysteric damping ratio	0.1
Sleeper layer	m_s	Mass per unit length	434 kg/m
Embankment layer	m_B	Mass per unit length	2000 kg/m
	k_B0	Undamped stiffness	5 × 10⁴ kN/m²
	η_B	Hysteric damping ratio	0.05
	2b_t	Contact width with the ground	5.0 m

Table 4.

Comparison of maximum Z-vibration levels.

Measuring point	Actual measurement	Simulation
P0	89.06	86.5
P1	57.93	58.55
P2	58.3	61.01
P3	57.94	63.44
P4	62.74	65.08
P5	62.18	66.01
P6	61.57	64.65

Analysis of vibration isolation effect of vibration isolation piles in time domain

To facilitate the analysis of the impact of isolation piles on soil vibrations under different factors, the number of pile rows and pile spacing were selected as the two variables for this study. The pile rows were set as single-row, double-row, and triple-row configurations. The pile radius was uniformly set to 1 m, and the pile spacing was set to 0.5 times, 1 time, 1.5 times, and 2 times the pile radius. Measurement point P6 was selected as the comparison point, with the edge of the isolation piles located 3 m from the tunnel. For clarity in analysis, we defined the legend in the figures such that the first digit represents the number of pile rows, and the second digit represents the pile spacing as a multiple of the pile radius. For example, “2-1.5” indicates a double-row pile configuration with a spacing of 1.5 times the pile radius. Additionally, the arrangement of the double-row and triple-row piles is shown in Figure 8. Additionally, in Figure 9, we present the time history of vibration acceleration at different depths of soil without isolation walls.

Figure 8.

Arrangement of piles. (a) Two row pile arrangement. (b) Three-row pile arrangement.

Figure 9.

Time history of vibration acceleration at different depths of soil without isolation piles. (a) 0 m. (b) −2 m. (c) −4 m. (d) −6 m. (e) −8 m. (f) −10 m.

Peak acceleration vibration isolation effect analysis of vibration isolation piles in time domain

As shown in Figure 10(a), it is observed that when using a single row of piles and the soil depth ranges from 0 to 6 m, the peak acceleration is minimized with a pile spacing of 2r. Conversely, for soil depths between 6 and 10 m, the peak acceleration is minimized with a pile spacing of 0.5r. Additionally, under the influence of a single row of isolation piles, the peak accelerations of soil at different depths are significantly lower than those without any isolation measures. Notably, at a soil depth of 10 m and a pile spacing of 0.5r, the peak acceleration is markedly lower than at other spacings, constituting only 43% of the peak acceleration without isolation.

Figure 10.

Shows the curves of peak acceleration varying with depth under the influence of single-row piles, double-row piles, and triple-row piles. (a) Single row of piles. (b) Two row of piles. (a) Three row of piles.

As shown in Figure 10(b), the peak acceleration is minimized when using double-row piles with soil depths ranging from 0 to 2 m and 6 to 10 m, and the pile spacing is 0.5r. Additionally, except for the case where the isolation pile spacing is 1r and the soil depth is 6 m, which closely approaches the peak acceleration without isolation measures, all other isolation measures significantly reduce the peak acceleration. According to Figure 10(c), the peak acceleration is minimized when using triple-row piles with soil depths between 4 and 8 m and a pile spacing of 1r.

According to Figure 10, it is observed that the peak ground acceleration is generally lower when the spacing between the vibration isolation piles is 2r compared to other spacings. In the absence of vibration isolation piles, the peak acceleration in the soil is highest at a depth of 8 m and lowest at a depth of 4 m. Under the influence of either single-row or double-row piles, the peak acceleration is highest at a depth of 10 m and lowest at a depth of 2 m. When using double-row piles, it is recommended to adopt a spacing of 0.5r, as this spacing results in peak accelerations that are lower than those of other pile spacings, except for a peak velocity that is larger at a soil depth of 4 m.

As shown in Figure 11(a), the vibration isolation effect of the triple-row pile is significantly better than that of the single-row and double-row piles when the soil depth ranges from 0 to 4 m. However, when the soil depth exceeds 4 m, the vibration isolation effect of the single-row pile is superior to that of the double-row and triple-row piles. According to Figure 11(b), the vibration isolation effect of the triple-row pile is significantly better than that of the single-row and double-row piles when the soil depth ranges from 4 to 8 m. However, at the ground surface and a soil depth of 10 m, the vibration isolation effect of the single-row pile is superior to that of the double-row and triple-row piles. As indicated in Figure 11(c), the vibration isolation effect of single-row piles is significantly better than that of double-row and triple-row piles when the soil depth ranges from 4 to 8 m. Similarly, Figure 11(d) shows that the vibration isolation effect of single-row piles is also superior to double-row and triple-row piles across a soil depth range of 2 to 10 m. It can be concluded that, contrary to expectations, the vibration isolation effect of multi-row piles decreases as the spacing between piles increases, falling below that of single-row piles.

Figure 11.

Shows the curves of peak acceleration varying with depth under the influence of different row numbers of vibration isolation piles. (a) The pile spacing is 0.5r. (b) The pile spacing is 1r. (c) The pile spacing is 1.5r. (d) The pile spacing is 2r.

According to Figure 11, when the spacing between piles is 0.5r, it is recommended to use single-row piles. This is because the peak accelerations at a soil depth of 10 m for double-row and triple-row piles are approximately 1.8 times greater than those of single-row piles. Moreover, the peak acceleration of triple-row piles at a soil depth of 6 m exceeds that of piles without vibration isolation. When the spacing between rows of piles is 1r, it is recommended to use a single row of piles, as the peak acceleration at any soil depth is smaller with a single row. When the spacing is 1.5r, either a single or two row of piles is suggested because the peak acceleration at a soil depth of 4 m with three rows of piles is greater than that without the effect of vibration isolation piles. When the spacing is 2r, a single row of piles is recommended, as the peak acceleration at any soil depth is smaller with a single row.

Peak velocity isolation effect analysis of vibration isolation piles in time domain

As indicated in Figure 12(a), the peak velocities of soil at various depths under the influence of single-row isolation piles are significantly lower than those without isolation measures. At a soil depth of 10 m, the peak velocity is notably lower when the pile spacing is 0.5r compared to other spacings, amounting to only 70% of the peak velocity without isolation. According to Figure 12(b), when double-row piles are employed and the soil depth ranges from 4 to 10 m, the peak velocity is minimized at a pile spacing of 1.5r. Additionally, except for the case where the isolation pile spacing is 1r, which shows an increase in peak velocities at the surface and at a soil depth of 6 m, all other isolation measures result in a significant reduction in peak velocities. As indicated in Figure 12(c), the peak velocity is minimized when using a three-row pile configuration with a soil depth of 4 to 8 m and a pile spacing of 1r. Additionally, an increase in peak velocity is observed at a soil depth of 6 m when the pile spacing is 0.5r and at the ground surface when the pile spacing is 1.5r.

Figure 12.

Illustrates the variation of peak velocity with depth under the influence of single-row piles, double-row piles, and triple-row piles. (a) Single row of piles. (b) Two row of piles. (c) Three row of piles.

According to Figure 12, both double-row and triple-row piles achieve the minimum peak velocity at a pile spacing of 0.5r when the soil depth is 2 m. When using single-row piles, it is recommended to set the spacing between the piles at 0.5 times r, as this spacing results in the minimum peak velocity at a soil depth of 10 m, with little variation in peak velocity at other soil depths. For double-row piles, it is advised to use a spacing of 0.5 and 1.5 times r, as this ensures the minimum peak velocity for any soil depth. In the case of triple-row piles, a spacing of 1 times r is recommended, as it yields the minimum peak velocity when the soil depth ranges from 4 to 8 m.

Based on Figure 13(a), it is recommended to use a single row of piles when the spacing between piles is 0.5r, as the vibration isolation effect of a single row of piles is significantly better than that of double or triple rows when the soil depth ranges from 4 to 10 m. However, at a soil depth of 8 m, the vibration isolation effect of two rows of piles is superior to that of a single row. According to Figure 13(b), when the spacing between piles is 1r, it is advisable to use a single row of piles, as the vibration isolation effect of a single row is significantly better than that of double or triple rows at any soil depth. However, under the influence of double-row piles at the ground surface and a soil depth of 6 m, a phenomenon of vibration enhancement was observed. As indicated in Figure 13(c), it is recommended to use single-row piles when the spacing between piles is 1.5r. This is because the vibration isolation effect of single-row piles is significantly better than that of double-row and triple-row piles within a soil depth range of 2 to 6 m. According to Figure 13(d), when the pile spacing is 2r, single-row piles are also recommended, as their vibration isolation effect is significantly superior to that of double-row and triple-row piles within a soil depth range of 2 to 8 m. However, the vibration isolation performance of double-row and triple-row piles is superior to that of single-row piles when the soil depth is 2 m and 10 m.

Figure 13.

Shows the curves of peak velocity varying with depth under the influence of different numbers of vibration isolation piles. (a) The pile spacing is 0.5r. (b) The pile spacing is 1r. (c) The pile spacing is 1.5r. (d) The pile spacing is 2r.

According to Figure 13, without considering the vibration impact at a soil depth of 2 m, single-row piles offer the best vibration isolation effect in terms of peak velocity in the time domain. When using double-row piles, the spacing between piles should avoid 1r. Double-row piles can also be employed for vibration isolation when the spacing is 1.5r, and their performance is better than single-row piles both at the ground surface and at the soil depth.

Peak displacement isolation effect analysis of vibration isolation piles in the time domain

As shown in Figure 14(a), the peak displacements of soil at different depths under the influence of single-row isolation piles are significantly smaller than those without isolation measures. Figure 14(b) indicates that when double-row piles are used with a spacing of 1r, an increase in vibration is observed within the soil depth range of 0 to 2 m, while peak displacements under other isolation measures are significantly reduced. Figure 14(c) reveals that when triple-row piles are employed, the peak displacements at the ground surface all exhibit an increase.

Figure 14.

Shows the curves of peak displacement varying with depth under the influence of single-row piles, double-row piles, and triple-row piles. (a) Single row of piles. (b) Two row of piles. (c) Three row of piles.

According to Figure 14, at a soil depth of 2 m, both double-row and triple-row piles achieve their minimum peak displacements when the spacing between the rows is 0.5r. When using single-row piles, it is recommended to set the spacing between the piles at 1r, as this spacing results in relatively small peak displacements at any soil depth. For double-row piles, the recommended spacing is 1.5r, which ensures minimal peak displacements at any soil depth and avoids displacement amplification. However, when considering displacement isolation, especially for isolation at a distance of 10 m from the isolation piles, it is not advisable to use triple-row piles.

As indicated in Figure 15(a), the vibration isolation effect of a single row of piles is significantly better than that of double or triple rows when the spacing between piles is 0.5r. However, when the soil depth is 8 m, the vibration isolation effect of two rows of piles is superior to that of single rows. According to Figure 15(b), when the pile spacing is 1r, the vibration isolation effect of a single row of piles is significantly better than that of double or triple rows. As shown in Figure 15(c), when the pile spacing is 1.5r, the vibration isolation effect of a single row of piles is significantly better than that of double or triple rows. However, when the soil depth is 8 m, the vibration isolation effect of double-row piles is superior to that of single-row piles. As shown in Figure 15(d), when the spacing between piles is 2r, the vibration isolation effect of single-row piles is significantly better than that of double-row and triple-row piles. However, when the soil depth is 10 m, the vibration isolation effect of double-row piles is superior to that of single-row piles.

Figure 15.

Shows the curves of peak displacement varying with depth under the influence of different row numbers of vibration isolation piles. (a) The pile spacing is 0.5r. (b) The pile spacing is 1r. (c) The pile spacing is 1.5r. (d) The pile spacing is 2r.

According to Figure 15, for peak displacement in the time domain, the vibration isolation effect of single-row piles is overall the best.

Discussion

Compared to recent studies, this research correlates pile spacing with the radius of the isolation piles. The primary focus of this study is on soil vibrations induced by train operations, rather than solely on surface vibrations. The study shows that the vibration isolation effectiveness of isolation piles varies with changes in pile spacing. Inappropriate pile spacing may even exacerbate soil vibrations. In terms of reducing vibration acceleration, single-row piles are more effective than double-row and triple-row piles in most cases, especially when the pile spacing is larger. For attenuating peak velocity, single-row piles consistently provide better vibration isolation at any pile spacing, while careful consideration of pile spacing is required for double-row and triple-row piles to avoid vibration amplification. Regarding peak displacement reduction, single-row piles are the most effective, particularly when the pile spacing is equal to the pile radius.

Due to the limitations of our study, we only investigated the impact of pile spacings of 0.5, 1, 1.5, and 2 times the pile radius on vibration. Future research could expand the range of pile spacings and explore the influence of varying front and rear pile spacings on vibration. Additionally, it is important to acknowledge that the current conclusions are applicable only to stratified soil with homogeneous composition within each layer. This limitation implies that the applicability of these findings may be affected in cases where the lateral geological structure is uneven.

Conclusions

This study first developed a numerical model that includes the train and track slab to address the wheel-rail interaction forces between the train and the rail. Subsequently, a finite element model with dynamic coupling of the track-tunnel-soil system was constructed, and the wheel-rail interaction forces were applied to it. The reliability of the model was validated by comparing the simulation results with measured data. The study then investigated the impact of isolation piles on soil vibrations under various influencing factors, leading to the following conclusions:

To reduce vibration acceleration in the soil within a depth range of 0 m to 2 m, the vibration isolation effectiveness decreases as pile spacing increases. When the pile spacing reaches 2 times the pile radius, single-row piles perform better. As the pile spacing increases beyond 1 time the pile radius, the vibration reduction effectiveness of double-row and triple-row piles in the 4 m to 8 m depth range becomes less effective than that of single-row piles. To attenuate peak acceleration in soils deeper than 6 m, the pile spacing should be set to 0.5 times the pile radius for single-row and double-row piles and to 1 time the pile radius for triple-row piles.

At any pile spacing, single-row piles generally provide better vibration isolation in terms of peak velocity. When using double-row piles, the spacing should not be set to 1 time the pile radius, as this may lead to amplification of peak velocity. To reduce peak velocity in soils deeper than 2 m, the pile spacing should be set to 0.5 times the pile radius for single-row and double-row piles and to 1 time the pile radius for triple-row piles.

To reduce peak soil displacement, single-row piles should be used, with the spacing set to 1 time the pile radius for optimal vibration isolation. When using double-row piles, the spacing should not be set to 1 time the pile radius, as this may result in displacement amplification.

Footnotes

Authors’ contribution

Shusen Cao and Dong Li contributed equally to this work and shared first authorship.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Scientific Research Fund of the Institute of Engineering Mechanics, China Earthquake Administration (Grant No. 2021D32), the Fund of State Key Laboratory of Bridge Engineering Structural Dynamics (Grant No. 202104), and the Heilongjiang Province Ecological Environment Protection Research Project (Grant No. HST2023JC010).

ORCID iD

Futong Wang

References

Liang

Luo

Zhou

, et al. Vibration filtering effect of a novel three-dimensional isolation bearing on metro vibration isolation. Eng Struct 2024; 301: 117304.

Connolly

Marecki

Kouroussis

, et al. The growth of railway ground vibration problems—A review. Sci Total Environ 2016; 568: 1276–1282.

Tao

Wang

Sanayei

, et al. Experimental study of train-induced vibration in over-track buildings in a metro depot. Eng Struct 2019; 198: 109473.

Masoumi

Hunt

. Assessment of building isolation performance against train induced vibrations by in situ measurement. J Acoust Soc Am 2022; 151: A226.

Wang

. Numerical analysis of isolation of the vibration due to moving loads using pile rows. J Sound Vib 2009; 319: 940–962.

Bouna

Nana Nbendjo

Woafo

. Isolation performance of a quasi-zero stiffness isolator in vibration isolation of a multi-span continuous beam bridge under pier base vibrating excitation. Nonlinear Dyn 2020; 100: 1125–1141.

Zhang

. Simultaneous topology optimization of supporting structure and loci of isolators in an active vibration isolation system. Comput Struct 2018; 194: 74–85.

Yang

Fang

, et al. Research on arrangement scheme of magnetic suspension isolator for multi-degree freedom vibration isolation system. J Ind Inf Integr 2017; 6: 47–55.

Xia

Sun

Chen

, et al. Analysis on multiple scattering by an arbitrary configuration of piles as barriers for vibration isolation. Soil Dyn Earthq Eng 2011; 31: 535–545.

10.

Mandal

Somala

. Periodic pile-soil system as a barrier for seismic surface waves. SN Applied Sciences 2020; 2: 1–8.

11.

Yao

Zhao

Zhang

, et al. Vibration isolation effect study of in-filled trench barriers to train-induced environmental vibrations. Soil Dyn Earthq Eng 2019; 125: 105741.

12.

Cai

Gao

, et al. The surface wave attenuation zone of periodic composite in-filled trenches and its isolation performance in train-induced ground vibration isolation. Comput Geotech 2021; 139: 104421.

13.

Wang

Gui

Meng

. Experiment study of vibration isolation by multi-rows of piles. IOP Conf Ser: Mater Sci Eng 2021; 634: 012047.

14.

Ding

Yang

, et al. Experimental study on the pile group-soil vibration induced by railway traffic under the inclined bedrock condition. Acta Geotech 2020; 15: 3613–3620.

15.

Kim

Yang

Hwang

, et al. Impact and blast mitigation using locally resonant woodpile metamaterials. Int J Impact Eng 2017; 101: 24–31.

16.

Amanat

Rafiee-Dehkharghani

Bitaraf

, et al. Analytical and numerical investigation of finite and infinite periodic lattices for mitigation of seismic waves in layered grounds. Int J Eng Sci 2022; 173: 103655.

17.

Cao

. Vibration isolation of row of piles embedded in transverse isotropic multi-layered soils. Comput Geotech 2018; 99: 115–129.

18.

Hui

Jue

, et al. Newly constructed subway on over-track bridge safety and vibration reduction measure. Adv Civ Eng 2024; 2024: 5851849.

19.

Liu

Y-Z

S-S

G-L

, et al. A field experiment on surface wave isolation by partially embedded periodic pile barriers. J Vib Control 2023; 29: 871–881.

20.

Luo

Deng

. A rubber-coated bored pile in isolating metro-induced vibrations: concept, analytical and in-situ investigation. Constr Build Mater 2024; 426: 136094.

21.

Jiang

Gao

Jian

, et al. Experimental study on attenuation zone of soil-periodic piles system. Soil Dyn Earthq Eng 2019; 126: 105738.

22.

Bin

Harvey Jr

. A dual-mode floor isolation system to achieve vibration isolation and absorption: experiments and theory. J Sound Vib 2022; 525: 116757.

23.

Zhou

Jia

Yuan

, et al. Vibration response prediction and vibration safety assessment method for building structures around viaducts under the action of high-speed railway trains. Constr Build Mater 2024; 422: 135877.

24.

Shi

, Surface-wave attenuation by periodic pile barriers in layered soils. Constr Build Mater 2018; 180: 177–187.

25.

Shi

. Feasibility of vibration mitigation in unsaturated soil by periodic pile barriers. Comput Geotech 2023; 164: 105798.

26.

Zhou

Wan

Wang

, et al. Novel periodic pile barriers in saturated soil and applications to propagation attenuation of shear plane waves. Comput Geotech 2022; 150: 104925.

27.

Saikia

. Numerical study on screening of surface waves using a pair of softer backfilled trenches. Soil Dyn Earthq Eng 2014; 65: 206–213.

28.

Sitharam

Sebastian

Fazil

. Vibration isolation of buildings housed with sensitive equipment using open trenches – case study and numerical simulations. Soil Dyn Earthq Eng 2018; 115: 344–351.

29.

Meng

Cheng

Shi

. Filtering property of periodic pile barriers under moving loads. Comput Geotech 2021; 136: 104244.

30.

Zhang

Wang

, et al. Study on vibration isolation effect of periodic pile arrangement on subway ground. IOP Conf Ser: Mater Sci Eng 2020; 455: 012093.

31.

Wang

Zhang

Jiang

, et al. Vibration characteristics and isolation in vibration-sensitive areas under moving vehicle load. Soil Dyn Earthq Eng 2022; 153: 107077.

32.

Zhang

Chen

Wei

, et al. Cast-in-place concrete thin-wall pipe pile as barrier for vibration isolation. J Cent South Univ Technol 2008; 15: 121–125.

33.

Zhou

Wan

Nian

, et al. Multi-objective topology optimization of periodic pile barriers with ambient body waves isolation performance. Mech Adv Mater Struct 2024; 31: 4254–4268.

34.

Hui

Hou

, et al. Newly constructed subway on over-track bridge safety and vibration reduction measure. Adv Civ Eng 2024; 2024: 5851849.

35.

Qiu

Zou

, et al. Prediction and mitigation of building vibrations caused by train operations on concrete floors. Appl Acoust 2024; 219: 109941.

36.

Tian

Zou

, et al. Train-induced vibration attenuation measurements and prediction from ground soil to building column. Environ Sci Pollut Res 2023; 30: 39076–39092.

37.

Qiu

Zou

, et al. Building vibration measurements induced by train operation on concrete floor. Constr Build Mater 2023; 394: 132283.

38.

Wang

Tao

Xie

, et al. Green’s function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic. Earthquake Eng Eng Vib 2017; 16: 311–328.

39.

Sheng

Jones

CJC

Thompson

D J

. A theoretical model for ground vibration from trains generated by vertical track irregularities. J Sound Vib 2004; 272: 937–965.

40.

Zhai

. Vehicle Track Coupled Dynamics. Beijing: Science Press, 2015.