Numerical prediction of ground vibrations induced by high-speed trains including wheel–rail–soil coupled effects

A simplified analytical model including the coupled effects of the wheel–rail–soil system and geometric irregularities of the track is proposed for evaluation of the moving train load. The wheel–rail–soil system is simulated as a series of moving point loads on an Euler–Bernoulli beam resting on a visco-elastic half-space, and the wave-number transform is adopted to derive the 2.5D finite element formulation. The numerical model is validated by published data in the literature. Numerical predictions of ground vibrations by using the proposed method are conducted at a site on the Qin-Shen Line in China.     

Dynamic analysis of a layered half-space subjected to moving line loads

Using a thin-layer method enhanced by continued-fraction absorbing boundary conditions, dynamic responses of a layered half-space subjected to a series of constant and time-harmonic line loads moving at a constant speed are studied. The thin-layer method for moving line loads is formulated for plane-strain as well as antiplane-shear conditions and is verified by comparison of computed responses of a homogeneous half-space subjected to a single constant load on its surface against available analytical solutions. Next, time-harmonic loads on a homogeneous half-space are examined. The study continues with both constant and time-harmonic loads on a layered half-space. Finally, multiple constant and time-harmonic loads are considered. The formulation and results demonstrate the effectiveness and versatility of the method in problems of dynamic response of layered media to moving loads.     

Dynamic responses of a saturated poroelastic half-space generated by a moving truck on the uneven pavement

A truck–pavement–ground coupling model was established to study the dynamic responses of a saturated poroelastic half-space generated by a moving heavy truck on the uneven pavement. The ground was simulated as a fully saturated poroelastic half-space governed by Biot’s theory. The overlying pavement was simplified as a Kirchhoff thin plate. With the assumption of a sinusoidal pavement surface, the dynamic wheel–pavement force was obtained through a linear Hertizian contact model. The numerical results showed that this dynamic load could make considerable contributions to the stress and excess pore water pressure responses in the ground. Furthermore, the effective stress path of the soil unit beneath the pavement caused by the moving truck was firstly calculated and presented. It was found that the differences between the total stress path and the effective stress path became significant as the truck speed increased, thus the effective stress path was more suitable than total stress path to reflect the stress history of soil elements in the saturated ground during the passage of high-speed traffics.     

Three-dimensional simulation of track on poroelastic half-space vibrations due to a moving point load

An analytical approach is used to investigate dynamic responses of a track system and the poroelastic half-space soil medium subjected to a moving point load under three-dimensional condition. The whole system is divided into two separately formulated substructures, the track sub-system and the ground. The ballast supporting rails and sleepers is placed on the surface of the ground. The rail is modeled by introducing the Green function for an infinitely long Euler beam subjected to the action of the moving point load and the reaction of sleepers represented by a continuous mass. Using the double Fourier transform, the governing equations of motion are then solved analytically in the frequency–wave-number domain. The time domain responses are evaluated by the inverse Fourier transform computation for a certain load velocities. Computed results show that dynamic responses of the soil medium are considerably affected by the fluid phase as well as the load velocity.     

移动荷载对桥梁动力响应的影响研究——以开裂预应力混凝土简支梁为例

混凝土桥梁在工作过程中会产生裂缝,为分析移动荷载对开裂混凝土桥梁结构刚度的影响,对开裂梁动力响应进行分析。建立简支T梁桥有限元模型,并将移动荷载施加至有限元模型中。根据简支T梁桥破坏横向分布位置和强度的不同,研究不同工况下各梁荷载横向分布及不同移动速度对裂缝扩展宽度的影响。结果表明,数值模拟结果能较好地验证计算模型的准确性;在较大的移动荷载作用下,混凝土开裂,导致结构刚度减小、位移增大;随着移动荷载和速度的增加,开裂时间增加,结构刚度降低,持续时间增加,位移增大,使结构响应呈现明显非线性。     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.     

Dynamic effects of moving loads on road pavements: A review

This review paper deals with the dynamic response of road pavements to moving loads on their surface. The road pavement can be modeled as a beam, a plate, or the top layer of a layered soil medium. The foundation soil can be modeled as a system of elastic springs and dashpots or a homogeneous or layered half-space. The material behavior of the pavement can be elastic or viscoelastic, while that of the foundation layers elastic, viscoelastic, water-saturated poroelastic or even inelastic. The loads are concentrated or distributed of finite extent, may vary with time and move with constant or variable speed. The analysis is done by analytical, analytical/numerical and purely numerical methods, such as finite element and boundary element methods, under conditions of plane strain or full three-dimensionality. A number of representative examples is presented in order to illustrate the problem and the methods of analysis, demonstrate the dynamic effects of moving loads on the layered soil medium and indicate the implications of the results on road and airport pavement design.     

三维层状黏弹性半空间中球面SH、P和SV波源自由场

采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称(柱面SH波)和轴对称(柱面P-SV波)情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数.     

列车引起的地基土应力状态变化的三维有限元分析

将列车移动荷载简化为多个移动轮轴荷载,基于列车-轨道-路基解析模型推求的列车运行时不同时刻、不同位置时作用于路基的振动荷载时程,采用多点输入方式实现列车荷载的移动施加方式,建立路(地)基-场地体系三维有限元动力分析模型,基于Abaqus软件的并行计算集群平台,对轨道交通振动荷载下路(地)基-场地体系的动力反应进行数值模拟,研究了列车荷载作用线正下方地基中的动应力特征及土单元应力状态变化,分析了列车轮轴荷载移动过程中不同深度处土单元的应力路径和主应力轴的旋转。     

A semi-analytical approach for analyzing ground vibrations caused by trains moving over elevated bridges

A semi-analytic approach is presented for the three-dimensional analysis of ground vibrations induced by trains moving over elevated bridges. The train is modeled as two sets of moving loads, with one for the front wheels and the other for the rear ones, the elevated bridge as a series of elastically supported beams, and the ground as a viscoelastic half space. Three key elements are considered in the solution: (1) the analytic solution for the vibration of an elastically supported beam under a series of moving loads, (2) the impedance of the foundation–soil system, and (3) Green's function for an elastic half space under a harmonic point load. Such an approach allows us to consider the structural dynamics of the elevated bridge, the foundation–soil interaction, and the wave propagation characteristics in the half space. From the numerical examples studied, the proposed approach was demonstrated to be accurate and efficient. The framework of analysis described herein can be generalized to solve problems with complex foundations and layered soils.     

Dynamic analysis of a transversely isotropic multilayered half-plane subjected to a moving load

In this paper, the analytical layer-element method is utilized to analyze the plane strain dynamic response of a transversely isotropic multilayered half-plane due to a moving load. We assume that the studied system moves synchronously with the moving load, then the moving load relative to the moving system is considered to be motionless. Therefore, the vertical stress and the vertical displacement under the moving load need not update for the variation of the load position. Based on the governing equations of motion in the moving system, the analytical layer-element solutions for a finite layer and a half-plane in the Fourier transform domain are derived by using the algebraic operations in Ref. [7]. The global matrix of the problem can be obtained by assembling the analytical layer-elements of all layers. The corresponding solution in the frequency domain is further recovered by the inverse Fourier transform. Several examples are given to confirm the accuracy of the proposed method and to illustrate the influence of material properties.     

土-结构动力相互作用问题分析中地震动输入的一种新方法

为实现地震作用下土-结构动力相互作用问题的有限元模拟,需要在人工边界上完成地震动的有效输入,目前工程和科研中常用的地震动输入方法有两种:波动输入方法和振动输入方法。波动输入方法的模拟精度高,但实施上相对复杂且耗时,而振动输入方法处理简单,但模拟精度较低。针对应力型人工边界提出一种在人工边界上实现地震动输入的新方法,该方法通过对土-结构有限元模型中由人工边界节点及相邻节点组成的局部子结构施加自由波场位移时程并进行动力分析,从而直接获得可实现地震波动有效输入的等效地震荷载,然后在土-结构有限元模型的人工边界节点上施加等效输入地震荷载并完成动力计算,由此完成土-结构动力相互作用问题的地震动输入和地震反应计算。与原有波动输入方法相比,新方法避免了原方法需分别计算人工边界上自由场应力和由引入人工边界条件引起的附加力,以及需要根据不同人工边界面的外法线方向确定荷载作用方向等较为复杂的处理过程,具有等效地震荷载计算简便、地震动输入过程更易于实施的特点。采用均匀弹性半空间和成层弹性半空间一维地震反应算例初步验证新方法的正确性和可靠性。     

Dynamic soil-tunnel interaction in layered half-space for incident P-and SV-waves

The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space,which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves.The indirect boundary element method is used,combined with the Green's function of distributed loads acting on inclined lines.It is shown that the dynamic characteristics of soil-tunnel interaction in layered half-space are different much from that in homogeneous half-space,and that the mechanism of soil-tunnel interaction is also different much from that of soil-foundation-superstructure interaction.For oblique incidence,the tunnel response for in-plane incident SV-waves is completely different from that for incident SH-waves,while the tunnel response for vertically incident SV-wave is very similar to that of vertically incident SH-wave.     

三维层状场地的精确动力刚度矩阵及格林函数

本文对Wolf二维层状场地精确动力刚度矩阵进行推广,给出了三维层状场地的精确动力刚度矩阵。刚度矩阵具有对称的特点,且因刚度矩阵是精确的,计算结果不受土层单元厚度的影响,可以大大提高计算效率。文中对刚度矩阵进行了数值验证。利用三维层状场地动力刚度矩阵,计算分析了基岩上单一土层场地的动力响应。最后作为动力刚度矩阵的另一重要应用,给出了表面或埋置矩形均布荷载或集中荷载的动力格林函数计算方法。     

Vibrations of the elastic half-space due to vertical surface loads

The displacements of the surface of the homogeneous elastic half-space subjected to a vertical concentrated load acting at the surface are calculated. The solution is shown to approach Lamb's simple formula for the far-field. The solution for the concentrated load can be used to calculate displacements due to distributed loads with high accuracy and economically. The results are applied to define the regions ‘near-field’ and ‘far-field’. Some near-field phenomena are described.     

The vibration frequencies of an elastically supported cantilever column with variable cross-section and distributed axial load

The solution for transverse vibration frequencies of an elastically supported cantilevered column with continuously varying cross-section under distributed axial load is developed. Boundary conditions are formulated for fixed and elastic support. Values of the fundamental frequency are tabulated for various cross-sections and axial load distributions due to self-weight as functions of the ratios between the rigidities and also the axial loads at the fixed and free ends, the ratio between the actual and critical buckling axial load, the degree of polynomial approximation and the rotational stiffness of the elastic support.     

Train-induced wave propagation in layered soils using finite/infinite element simulation

In this paper, the transmissibility of soils for vibrations induced by trains moving at different speeds is studied. The 2.5 D finite/infinite element approach adopted herein allows us to consider the load-moving effect of the train in the direction normal to the two-dimensional profile of the soils considered, and, therefore, to obtain three-dimensional responses for the soils using only plane elements. The moving train is simulated by a sequence of moving wheel loads that may vibrate with certain frequency. Two train speeds are considered, one is smaller and the other is greater than the Rayleigh wave speed of the layered soils, to represent the effects of speed in the sub-critical and super-critical ranges. In order to evaluate the effect of each parameter on the ground response induced by moving trains, parametric studies are conducted for the following parameters: the shear wave speed, damping ratio and stratum depth of the supporting soils, and the moving speed and vibration frequency of the traveling trains. Conclusions concerning the mechanism of wave propagation in layered soils are drawn from the parametric studies, which should prove useful to practicing engineers.     

瑞雷波斜入射下层状半空间中沉积谷地周围的三维散射研究

针对层状半空间中沉积谷地对斜入射瑞雷波的三维散射问题, 采用直接刚度法计算自由场波场, 以层状半空间中移动斜线均布荷载动力格林影响函数求解三维散射波场, 建立了求解该问题的间接边界元方法. 通过与已有结果的比较, 验证了该方法的正确性, 并以均匀半空间以及弹性基岩上单一土层场地中沉积谷地为例进行了计算分析. 研究结果表明: 层状半空间与均匀半空间中沉积谷地对瑞雷波的散射存在显著差别; 层状半空间中瑞雷波的振动模态对沉积附近位移幅值有着重要影响; 土层刚度和厚度等参数也对沉积附近位移幅值大小及空间分布有着显著的影响.      

Soil–structure interaction in the time domain using halfspace Green's functions

A time-domain formulation is proposed for the transient response analysis of general, three-dimensional structures resting on a homogeneous, elastic halfspace subjected to either external loads or seismic motions. The formulation consists of two parts: (a) the time domain formulation of the soil behaviour and (b) the coupling of the corresponding soil algorithms to the Finite Element Code ANSYS. As far as the structure is concerned, this coupling opens the way for the analysis of non-linear soil–structure interaction. The approach is based on halfspace Green's functions for displacements elicited by Heaviside time-dependent surface point loads. Hence, the spatial discretisation can be confined to the contact area between the foundation and the soil, i.e. no auxiliary grid beyond the foundation as for conventional boundary element formulations is required. The method is applied to analyse the dynamic response of a railway track due to a moving wheel set by demonstrating the influence of ‘through-the-soil coupling’.     

Surface vibration due to a sequence of high speed moving harmonic rectangular loads

The transmission of vibrations over the surface of the ground, due to high-speed moving, vertical harmonic rectangular loads, is investigated theoretically. The problem is three-dimensional and the interior of the ground is modelled as an elastic half-space or a multilayered ground. The transformed solutions are obtained using the Fourier transform on the space variable. A new damping model in the spatial wavenumber domain, presented in Lefeuve-Mesgouez et al. [J. Sound. Vibr. 231 (2000) 1289] is used. Numerical results for the displacements on the surface are presented for loads moving with speeds up to and beyond the Rayleigh wave speed of the half-space.