Numerical analysis of soil vibrations due to trains moving at critical speed

High-speed train induced vibrations of track structure and underlying soils differ from that induced by low-speed train. Determining the critical speed of train operation remains difficult due to the complex properties of the track, embankment and ground. A dynamic analysis model comprising track, embankment and layered ground was presented based on the two-and-half-dimensional (2.5D) finite elements combining with thin-layer elements to predict vibrations generated by train moving loads. The track structure is modeled as an Euler–Bernoulli beam resting on embankment. The train is treated as a series of moving axle loads; the embankment and ground are modeled by the 2.5D finite elements. The dynamic responses of the track structure and the ground under constant and vibrating moving loads at various speeds are presented. The results show that the critical speed of a train moving on an embankment is higher than the Rayleigh wave velocity of the underlying soil, attributed to the presence of the track structure and the embankment. It is found that the dynamic response of ground induced by moving constant loads is mostly dominated by train speed. While for the moving load with vibration frequency, the ground response is mostly affected by the vibration frequency instead of train speed. Mach effect appears when the train speed exceeds the critical speed of the track–embankment–ground system.     

列车运行时由轨道不平顺引起的地基振动研究

采用半解析法研究了列车荷载作用下列车-轨道-饱和地基系统的耦合振动问题。研究模型共分为3部分：车体简化为一个多刚体系统,在车轮与钢轨之间引入线性Hertizian弹簧接触模型模拟轮轨动力相互作用、采用离散轨枕支撑的弹性Euler梁来模拟轨道系统、下卧土体采用多孔饱和半空间模型。列车荷载分为轴重和由轨道不平顺引起的轮轨动力相互作用力。采用Fourier变换分别求解各子系统的控制方程,并通过动力子结构法对各子系统进行耦合。土体在时域内的动力响应通过快速Fourier变换求得。在分析了轮轨动力相互作用力的基础上,研究了轮轨动力作用力和列车轴重作用下饱和地基的动力响应,并分析了轨枕间距和土体渗透系数对饱和地基振动响应的影响。研究表明,轮轨动力作用力对地基远场振动有重要贡献,同时枕木间距对轨道与地基振动响应有较大影响。     

Three‐dimensional steady‐state response of a railway system on layered half‐space soil medium subjected to a moving train

The paper presents a three‐dimensional study on the steady‐state response of a track system and layered half‐space soil medium subjected to the load induced by the passages of a moving train with the substructure method. Practically, due to the ground water table being several meters beneath the ground surface, the soil profile can be divided into two layers: the upper layer modeled by an elastic medium and the lower layer by a fully saturated poroelastic medium governed by Biot's theory. In the meanwhile, the rails are regarded as an infinitely long Euler beam, and the sleepers are represented by a continuous mass. The effect of the ballast is accounted for by introducing the Cosserat model for granular medium, and the train is described by a series of moving axle point loads, depending on the geometry of the train. The influences of the thickness, the mass and the rigidity of the elastic layer and the mass of the ballast on rail's displacement responses are carefully investigated. Numerical results show that the influences of these parameters are significant for high train velocity, while vanishes for low velocity. Copyright © 2008 John Wiley & Sons, Ltd.     

饱和土体中排桩对移动荷载的被动隔振效果分析

利用Muki和Sternberg的虚拟桩法,研究了饱和土体地基中排桩对移动荷载引起振动的被动隔振效果。隔振桩作为一维杆,饱和土体满足Biot理论。利用已有的移动荷载作用下的饱和土体的自由波场解及饱和土体内部受竖向圆形分布荷载作用下的基本解,建立了频域内土-桩的第2类Fredholm积分方程。通过Fourier逆变换得到时间域内评价隔振效果的振幅比。与已知文献结果相比较,验证了方法的正确性。数值结果表明,荷载速度对排桩的隔振效果有一定影响,即在相同隔振系统情况下,单排桩对低速荷载引起振动的隔振效果比高速移动荷载效果好。同时,较高速时的最佳桩长比低速时要短。     

Dynamic responses of unsaturated half‐space soil to a moving harmonic rectangular load

The problem of the dynamic responses of a semi‐infinite unsaturated poroelastic medium subjected to a moving rectangular load is investigated analytical/numerically. The dynamic governing equations are obtained with consideration of the compressibility of solid grain and pore fluid, inertial coupling, and viscous drag as well as capillary pressure in the unsaturated soil, and they can be easily degraded to the complete Biot's theory. Using the Fourier transform, the general solution for the equations is derived in the transformed domain, and then a corresponding boundary value problem is formulated. By introducing fast Fourier transform algorithm, the unsaturated soil vertical displacements, effective stresses, and pore pressures induced by moving load are computed, and some of the calculated results are compared with those for the degenerated solution of saturated soils and confirmed. The influences of the saturation, the load speed, and excitation frequency on the response of the unsaturated half‐space soil are investigated. The numerical results reveal that the effects of these parameters on the dynamic response of the unsaturated soil are significant.     

Vertical amplitude reduction of Rayleigh waves by a row of piles in a poroelastic half‐space

A row of rigid piles is addressed as the countermeasures for isolating Rayleigh waves in a poroelastic half‐space. The complex characteristic equations for Rayleigh waves are derived via Biot's theory and their existence conditions are given. The piles are modeled as Euler–Bernoulli beams with longitudinal displacements and the diffracted field by each pile is constructed only with Rayleigh waves. Six infinite linear systems of algebraic equations are obtained in terms of the equilibrium of forces and continuity of displacements at the pile–soil interfaces. The systems are subsequently solved in the complex least‐squares sense. The influence of certain pile and soil characteristics such as the permeability of poroelastic soil, spacing between the piles and length of the piles on the isolating performance of a pile barrier is investigated. Computed results show that the permeability of poroelastic soil displays a significant effect on the vertical amplitude reduction of Rayleigh waves. Copyright © 2009 John Wiley & Sons, Ltd.     

Three‐dimensional dynamic response of ground with a poroviscoelastic interlayer to a harmonic moving rectangular load

Stratification is a basic characteristic of ground. Due to the influence of ground water, saturated weak interlayers widely exist, particularly in soft soil area. An interlayer of high compressibility and low strength has a substantial effect on dynamic response of the ground, especially under high speed moving load. Thus, a comprehensive investigation in the influence of interlayer is essential and useful in geotechnical and transportation‐related engineering. This paper presents a three‐dimensional semi‐analytical approach to study the dynamic response of a layered ground with a soft saturated interlayer. The ground is modelled as a half‐space consisting of three parts: a viscoelastic upper layer, a saturated poroviscoelastic interlayer governed by Biot's theory and a viscoelastic half‐space. An ‘adapted stiffness matrix’ is proposed to obtain the semi‐analytical solution to the system. Comprehensive parametric study is conducted to investigate the influences of existence, geometrical and physical properties of the interlayer. Depth, thickness, hydraulic permeability of the interlayer, load speed and frequency significantly influence the dynamic response of the ground, among which the interlayer depth plays a dominant role. Resonant frequency exists, which is highly affected by the interlayer thickness, especially in low speed regime. Both hydraulic permeability and boundary conditions of the interlayer influence the characteristics of pore pressure distribution. Copyright © 2017 John Wiley & Sons, Ltd.     

空沟对列车运行引起的地基振动隔振效果研究

为研究空沟对高速列车引起地基振动的隔振效果,在已有研究基础上改进出一种新的解析研究模型。模型中首次运用了饱和半空间模型来研究地基上隔振沟对高速列车的隔振效果;隔振沟则通过在饱和半空间土体上设置3个合适宽度、截面为矩形的弹性层来模拟;中间矩形弹性层为路堤,路堤上方放置了枕木与轨道。枕木与轨道分别通过纵向异性Kirchhoff薄板与Euler梁来模拟,饱和土地基采用Biot多孔饱和介质理论来描述。控制方程通过傅立叶变换与傅立叶级数展开,在变换域中进行求解。研究表明,随着列车运行速度的提高,空沟的隔振效果明显提高;饱和土体固-液相的耦合作用对隔振沟的隔振效果的影响明显,尤其当列车运行速度超过土体表面Rayleigh波速时,随着土体渗透系数的增加,空沟隔振效果显著降低。此外,列车运行速度超过土体表面波速时,饱和土地基上空沟的隔振效果明显优于相应单相弹性地基上空沟的隔振效果。     

Rocking vibrations of a rigid embedded foundation in a poroelastic half‐space

This paper presents an analysis of the rocking vibrations of a rigid cylindrical foundation embedded in poroelastic soil. The foundation is subjected to time‐harmonic rocking excitation and is perfectly bonded to the surrounding soil. The soil underlying the foundation base is represented by a homogeneous poroelastic half‐space, whereas the soil along the side of the foundation is modeled as an independent poroelastic stratum composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory. The contact surface between the foundation base and the poroelastic soil is assumed to be smooth and either fully permeable or impermeable. The dynamic interaction problem is solved by employing a simplified analytical method. Some numerical results for the nondimensional rocking dynamic impedance and nondimensional angular displacement amplitude of the foundation are presented to show the effect of nondimensional frequency of excitation, poroelastic material parameters, hydraulic boundary condition, depth ratio and mass ratio of the foundation. Copyright © 2009 John Wiley & Sons, Ltd.     

The influence of the depth of the ground water table on free field road traffic-induced vibrations

This paper describes the influence of seasonal variations of the ground water table on free field traffic-induced vibrations. The passage of a truck on two types of road unevenness is considered: a joint in a road pavement consisting of concrete plates and a speed bump with a sinusoidal profile. Free field vibrations are computed with a two-step solution procedure, where the computation of the vehicle axle loads is decoupled from the solution of the road–soil interaction problem. The impedance of the soil is calculated using a boundary element method, based on the Green's functions for a dry layer on top of a saturated half-space. It is demonstrated that, in the low-frequency range of interest, wave propagation in the saturated half-space can be modelled with an equivalent single phase medium as an alternative to Biot's poroelastic theory for saturated porous media. The relation between the free field velocity and the depth of the ground water table is dominated by three phenomena: (1) the compressibility of the soil decreases due to the presence of the pore water, (2) the ground water table introduces a layering of the soil which may cause resonance of the dry layer and (3) refracted P-waves in the dry layer interfere with surface waves. If the depth of the ground water table is large with respect to the wavelength of the vibrations in the soil, the response tends to the response of a dry half-space. The average free field velocity is equal to the velocity in the absence of ground water. If the depth of the ground water table is small with respect to the wavelength of the vibrations in the soil, the response tends to the response of a saturated half-space and resonance of the dry layer does not occur. The average free field velocity is smaller than the velocity in the absence of ground water. The interference of refracted P-waves and surface waves causes an additional oscillation of the response as a function of the excitation frequency and the distance between the road and the receiver. Copyright © 2004 John Wiley & Sons, Ltd.     

Poroelastic model for pile–soil interaction in a half‐space porous medium due to seismic waves

In this paper, frequency domain dynamic response of a pile embedded in a half‐space porous medium and subjected to P, SV seismic waves is investigated. According to the fictitious pile methodology, the problem is decomposed into an extended poroelastic half‐space and a fictitious pile. The extended porous half‐space is described by Biot's theory, while the fictitious pile is treated as a bar and a beam and described by the conventional 1‐D structure vibration theory. Using the Hankel transformation method, the fundamental solutions for a half‐space porous medium subjected to a vertical or a horizontal circular patch load are established. Based on the obtained fundamental solutions and free wave fields, the second kind of Fredholm integral equations describing the vertical and the horizontal interaction between the pile and the poroelastic half‐space are established. Solution of the integral equations yields the dynamic response of the pile to plane P, SV waves. Numerical results show the parameters of the porous medium, the pile and incident waves have direct influences on the dynamic response of the pile–half‐space system. Significant differences between conventional single‐phase elastic model and the poroelastic model for the surrounding medium of the pile are found. Copyright © 2007 John Wiley & Sons, Ltd.     

Impedances of rigid cylindrical foundations embedded in transversely isotropic soils

A complete formulation and implementation for assessment of the response to dynamic loads of cylindrical rigid structures embedded in transversely isotropic elastic half‐spaces is presented. The analysis is performed in the frequency domain and the steady‐state structure response is obtained. The method is based on a non‐singular version of the indirect boundary element method which uses influence functions, instead of Green's functions, as fundamental solutions. These influence functions are the response of an elastic half‐space to distributed, internally applied loads. The proposed method imposes full bonding contact between the foundation and the surrounding soil. Numerical results for displacement (vertical and horizontal) and rotation (twisting and rocking) impedances, showing the influence of the soil anisotropy, are presented. Results for the soil–structure interface tractions and for the displacement field throughout the half‐space are also shown. Copyright © 2006 John Wiley & Sons, Ltd.     

移动荷载作用下饱和土中单桩的动力响应

根据Biot动力理论,采用Fourier和Hankel变换方法得到了半空间饱和土受移动载荷及土体内受垂直简谐载荷作用下频域内基本解。根据虚拟桩法,得到了移动载荷作用下桩基的第2类Fredholm积分方程,并应用IFFT方法得到时间、空间域内单桩的动力响应。数值结果表明,移动荷载会引起桩身的负摩擦力;桩身最大轴力、孔压随移动荷载速度增加而增大;此外,在桩上端部会出现孔压集中现象。     

层状饱和地基上刚性圆板的扭转振动及参数分析

用解析的方法首次研究了层状饱和地基上刚性圆板的扭转振动特性.运用Hankel变换求解饱和介质动力问题的控制方程,按混合边值条件建立了层状饱和地基上刚性圆板扭转振动的第二类Fredholm积分方程.数值算例给出了动力柔度系数和扭转角幅值随无量纲频率的变化曲线,与单相弹性及匀质饱和介质情况进行了对比分析,并进行了参数研究.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface

By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic response of pavements on poroelastic half-space soil medium to a moving traffic load

This paper presents an investigation of the steady-state response of pavement systems subjected to a moving traffic load. The traffic loads are simulated by four rectangular load pressures, and the rigid and flexible pavement systems are regarded as an infinite plate resting on a poroelastic half-space soil medium. The contact surface between the plate and the poroelastic half-space is assumed to be smooth and fully permeable. Kirchhoff small-deflection thin-plate theory is employed to analyze the plate, while Biot’s fully dynamic poroelastic theory is used to characterize the poroelastic half-space. The frequency wave-number domain solution of the pavement system is obtained by the compatibility condition between the plate and the poroelastic half-space. By applying the inverse fast Fourier transform, the time domain solution is obtained. Also, the influences of the load speed, the permeability of the soil, and the flexural rigidity of the plate on the response of the pavement system are investigated. The numerical results show that the influences of these parameters on the dynamic response of the pavement system are significant.     

A stepwise damping‐solvent extraction method for large‐scale dynamic soil–structure interaction analysis in time domain

Time‐domain analysis of dynamic soil–structure interaction based on the substructure method plays an increasing role in practical applications as compared with the frequency‐domain analysis. Efficient and accurate modelling of the unbounded soil or rock medium has been a key issue in such an analysis. This paper presents a subregional stepwise damping‐solvent extraction formulation for solving large‐scale dynamic soil–structure problems in the time domain. Accuracy and efficiency of the formulation are evaluated in detail for a classical problem involving a rigid strip foundation embedded in a half‐space. A practical large‐scale soil–structure interaction problem, which represents a high concrete gravity dam subjected to seismic load, is then analysed using the proposed method. Various responses of the dam, including time histories of the crest displacement and acceleration and contours of the peak principal stresses within the dam body, are presented. Comparisons are also made between these results with those obtained using other models for the unbounded medium. Copyright © 2007 John Wiley & Sons, Ltd.     

列车荷载下轨道系统-层状横观各向同性饱和地基动力响应

基于Biot波动理论,构建列车荷载-轨道系统-双层状横观各向同性饱和地基模型,将模型分为上覆路轨系统和地层系统。对上覆路轨系统和地层系统处理,并利用双重Fourier变换技术,在变换域中将横观各向同性饱和地基动力响应的求解简化为求解一个6阶控制方程的特征值问题,进而得到了列车荷载作用下双层横观各向同性饱和地基力响应的解析结果。利用离散Fourier逆变换得到数值计算结果,重点分析了上下土层的刚度和泊松比对位移和孔隙水压力和剪切应力响应的影响,结果表明,上、下土层刚度差异对地基动力响应有较大影响,土层各向异性参数中模量的影响较泊松比大。计算结果可为软土路基加固深度的确定提供理论依据。     

横观各向同性饱和地基竖向振动的衰减特性

基于提出的横观各向同性饱和多孔介质Biot波动方程的一般解,研究了饱和半空间地基在竖向点源简谐激振荷载作用下地表振动的衰减特征,分析了激振频率以及横观各向同性饱和土介质的各向异性参数和孔隙渗透系数对地表振动特征的影响。计算结果表明,低频和高频激振时,地表位移衰减特性存在明显差异;在饱和土的各向异性参数中,纵向和水平方向动态渗透系数比值和刚度系数比值对地表位移衰减影响最大,这也说明采用各向同性饱和介质的动力学模型不能准确地描述具有明显各向异性特性的饱和土地基的动力特性。