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

将列车移动荷载简化为多个移动轮轴荷载,基于列车-轨道-路基解析模型推求的列车运行时不同时刻、不同位置时作用于路基的振动荷载时程,采用多点输入方式实现列车荷载的移动施加方式,建立路(地)基-场地体系三维有限元动力分析模型,基于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.     

A transmitting boundary for the dynamic finite element analysis of cross anisotropic soils

The problem of the transmitting boundary used in the dynamic finite element analysis of layered axisymmetric soil models with non-axisymmetric displacements is considered. Each layer is modelled as a homogeneous, viscoelastic cross anisotropic medium with a vertical axis of material symmetry.     

Dynamic soil–structure interaction of coupled shear walls by boundary element method

For a class of civil engineering structures, that can be accurately represented by ‘coupled shear walls’ (CSWs), a discrete model for the analysis of the dynamic interaction with the underlying soil is proposed. The CSWs, with one or more rows of openings, rest on a rigid foundation embedded in the elastic or viscoelastic half-space. A hierarchical finite element model based on an equivalent continuum approach is adopted for the structure. A frequency-domain boundary element method is used to represent the half-space. Finally, the set of equations governing the response of the coupled soil-structure system to harmonic lateral loads acting on the structure is also given. The frequency deviation effect with respect to the fixed-base structure and the effects of radiation and material damping in the soil are presented for different characteristics of the structure and different soil properties.     

Numerical prediction of low-frequency ground vibrations induced by high-speed trains at Ledsgaard, Sweden

The ground vibrations induced by a passenger train at the test site of Ledsgaard, Sweden, have been analysed and numerically simulated through a spectral element discretization of the soil. To calculate the spatial distribution of loading due to train passage, the train is decoupled from the track, and a suitable series of static forces is applied. The track and the embankment are modeled as a beam on elastic foundation, using analytical solutions for loads moving at constant velocity. The results of both 2D and 3D modelling assumptions are thoroughly discussed, in terms of prediction of track motion and of attenuation of peak ground velocity with distance.     

A numerical model for the dynamic analysis of inclined pile groups

The paper presents a numerical model for the dynamic analysis of pile groups with inclined piles in horizontally layered soil deposits. Piles are modelled with Euler–Bernoulli beams, while the soil is supposed to be constituted by independent infinite viscoelastic horizontal layers. The pile–soil–pile interaction as well as the hysteretic and geometric damping is taken into account by means of two‐dimensional elastodynamic Green's functions. Piles cap is considered by introducing a rigid constraint; the condensation of the problem permits a consistent derivation of both the dynamic impedance matrix of the soil–foundation system and the foundation input motion. These quantities are those used to perform inertial soil–structure interaction analyses in the framework of the substructure approach. Furthermore, the model allows evaluating the kinematic stress resultants in piles resulting from waves propagating in the soil deposit, taking into account the pile–soil–pile interactions. The model validation is carried out by performing accuracy analyses and comparing results in terms of dynamic impedance functions, kinematic response parameters and pile stress resultants, with those furnished by 3D refined finite element models. To this purpose, classical elastodynamic solutions are adopted to define the soil–pile interaction problem. The model results in low computational demands without significant loss of precision, compared with more rigorous approaches or refined finite element models. Copyright © 2015 John Wiley & Sons, Ltd.     

Impedance of flexible suction caissons

The dynamic response of offshore wind turbines is affected by the properties of the foundation and the subsoil. The aim of this paper is to evaluate the dynamic soil–structure interaction of suction caissons for offshore wind turbines. The investigations include evaluation of the vertical and coupled sliding–rocking vibrations, influence of the foundation geometry and examination on the properties of the surrounding soil. The soil is simplified as a homogenous linear viscoelastic material and the dynamic stiffness of the suction caisson is expressed in terms of dimensionless frequency‐dependent coefficients corresponding to different degrees of freedom. The dynamic stiffness coefficients for the skirted foundation are evaluated using a three‐dimensional coupled boundary element/finite element model. Comparisons with known analytical and numerical solutions indicate that the static and dynamic behaviours of the foundation are predicted accurately using the applied model. The analysis has been carried out for different combinations of the skirt length, Poisson's ratio of the subsoil and the ratio of the soil stiffness to the skirt stiffness. Copyright © 2007 John Wiley & Sons, Ltd.     

高速铁路桥梁-桥墩-桩基础-地基耦合系统的地震反应

高速铁路中的桥梁常采用灌注桩基础以控制沉降,地震作用是桩基础的设计工况之一。建立桥梁-桥墩-桩基础-地基为一体的耦合系统非线性三维数值分析模型,以典型地震波为输入,考虑上部结构和基础的共同工作、土-结构动力相互作用、材料非线性和土层对桩的侧阻及端阻作用,开展三向地震作用下的动力有限元计算,并对地基主要土层压缩模量、桩体材料弹性模量、桩径和桩长进行参数敏感性分析。计算结果表明:现行的桩基础设计方案能有效控制地震荷载作用下桥梁的变形;地震过程中的不同时刻,桩侧阻发挥程度不同且不可忽略,以单纯的梁单元模拟桩的动力学行为的适用性值得商榷;桩长和地基主要土层压缩模量对桥梁地震反应影响最大,桩体材料弹性模量的影响次之,桩径的影响最小。     

Dynamic analysis of large 3-D underground structures by the bem

An up to date literature survey on the dynamics of underground structures is presented briefly. The dynamic response of large three-dimensional underground structures to external or internal dynamic forces or to seismic waves is numerically determined by the frequency domain boundary element method. This method is used to model both the structure and the soil medium, which are assumed to behave as linear elastic or viscoelastic bodies. The full-space dynamic fundamental solution is employed in the formulation and this requires a free soil surface discretization, confined to a finite portion around the area of interest, in addition to soil—structure interface and free structural surface discretizations. The dynamic disturbances can have a harmonic or a transient time variation. The transient case is treated with the aid of numerical Laplace transforms with respect to time. Various numerical examples involving lined cavities and long lined tunnels buried in the full- or the half-space subjected to harmonic or transient external forces or seismic waves are presented to illustrate the method and demonstrate its advantages.     

基于高阶双渐近透射边界的重力坝-层状地基动力相互作用分析

时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析.     

水平向不均匀土对箱基动力阻抗的影响

采用薄层元素法和有限单元法,建立了地基-箱形基础动力相互作用的三维分析模型,该模型可考虑箱形基础周围土的不均匀性影响.利用该模型,分析了水平方向不均匀土对箱形基础的水平和摇摆阻抗的影响.分析了水平方向不均匀土的厚度、剪切弹性模量以及材料阻尼比的影响.分析结果表明:水平方向不均匀土使得箱形基础的各种阻抗均降低,特别是阻抗的虚部.水平方向不均匀土对摇摆阻抗的影响较大.水平方向不均匀土的剪切弹性模量的影响较大而其材料阻尼的影响较小.随着水平方向不均匀土的厚度的增加和基础埋深的增加,其影响增大.     

Numerical modelling of free field traffic-induced vibrations

This paper deals with the numerical modelling of free field traffic-induced vibrations during the passage of a vehicle on an uneven road. The road unevenness subjects the vehicle to vertical oscillations that cause dynamic axle loads. The latter are calculated from the vehicle transfer functions and the frequency content of the road profile as experienced by the vehicle axles. A transfer function between the source and the receiver that accounts for the dynamic interaction between the road and the soil is used to calculate the free field response. Its calculation is based on a dynamic substructure method, using a boundary element method for the soil and an analytical beam model for the road. The methodology is validated with analytical results and is finally illustrated by a numerical example where the free field vibrations during the passage of a vehicle on a traffic plateau are considered.     

结构-基础-土-河流相互作用地震反应分析

为了研究不同基础形式下河水对成层软土地基上建筑物地震反应的影响,构建了结构-基础-土-河水体系动力相互作用的非线性完全有限元计算模型,在考虑土体重力的前提下,对筏片基础、桩筏基础、箱形基础、桩箱基础相互作用体系进行了时域数值分析。计算结果表明,建筑物靠近河水时,筏基和箱基体系顶部位移发生了向河水一侧的偏移,而采用桩基和桩箱基础时其偏移度明显减小。同种基础形式下,建筑物距河水越远,其框架剪力峰值越大,而不同基础形式的同种工况,桩箱基础时柱剪力最大,筏片基础时最小。单纯筏基和箱基时,基础周围土体均大量进入塑性,上部结构仍保持弹性;而采用下部有桩基础的筏基或箱基时,地基土体只有很少量进入塑性,上部框架结构进入了塑性。     

地铁列车荷载作用下黏弹性地基动应力特征分析

基于Kelvin空间半无限体黏弹性解,以单个轮轴荷载为例,研究地铁列车荷载下地基土动应力状态、应力路径及主应力轴旋转等应力变化规律,并与相应弹性解下的各动应力特征进行对比分析,探讨更适用于工程实际的地基附加动应力特征解。结果表明:弹性地基与黏弹性地基中的各应力状态变化总体趋势、主应力轴旋转变化趋势相同,但极限应力状态各应力分量大小及应力路径形状有明显差别,考虑黏弹性解更符合软土区地铁长期变形实际工况;并对单轮荷载下动应力特征的黏弹性解进行参数分析,发现特征系数对应力分量影响具有不同特性。     

The effect of road unevenness on the dynamic vehicle response and ground-borne vibrations due to road traffic

In this paper, the relation between road unevenness, the dynamic vehicle response, and ground-borne vibrations is studied. In situ measurements of road unevenness and the dynamic vehicle response for six roads with different types of pavement are supplemented by numerical predictions of ground vibrations. The predictions are performed in two stages. In the first stage, the dynamic vehicle response is computed based on the measured road unevenness. The vehicle model is validated by comparing the predicted and measured vehicle response and subsequently used to predict the dynamic vehicle loads. In the second stage, the dynamic road–soil interaction problem is considered and the transfer functions between the road and the soil are computed. The effect of the pavement type (continuous, jointed, or composed of individual pavers) on the road–soil transfer functions is investigated and the free field vibrations are calculated using the dynamic vehicle loads computed in the first stage. The predicted free field vibrations are validated by measurements at one of the measurement sites before and after rehabilitation of a deteriorated concrete pavement. Finally, the results are used to investigate the relation between indicators of road unevenness such as the ISO 8608 road class, the International Roughness Index, and the coefficient of evenness on one hand, and the dynamic vehicle response and level of ground-borne vibration on the other hand.     

Displacement solutions for dynamic loads in transversely-isotropic stratified media

Solutions for the displacements caused by dynamic loads in a viscoelastic transversely-isotropic medium are derived. The medium extends horizontally to infinity, but is bounded below by a rigid base. Stratification of the medium presents no difficulties. The medium is discretized in the vertical direction only; discretization in the horizontal direction is obviated by use of analytical solutions to the equations of motion. Application of the displacement solutions to soil-structure interaction is illustrated. A soil flexibility matrix (and hence, a stiffness matrix) for a surface foundation follows directly from the displacement solutions. A simple modification to obtain the soil stiffness for an embedded foundation of arbitrary geometry is described. Stiffnesses of rigid surface and embedded foundations are computed and compared with previously published results. In addition, the dynamic stiffness of a rigid surface foundation on a soil layer with linearly increasing shear modulus is compared to that for a homogeneous soil layer. A reduction in radiation damping is found to result from the inhomogeneity.     

Hamilton体系及弹性波在层状介质中的传播问题

利用结构力学与最优控制的模拟理论,研究弹性波在层状介质中传播的数值计算方法. 将弹性波传播问题导向哈密顿(Hamilton)体系,在哈密顿体系中,推导出一种新的半解析单元,称之为动力-部分杂交元,由此导出一套哈密顿体系下的半解析数值计算方法. 本文给出了该方法在层状正交各向异性材料介质的弹性波传播问题的数值算例,分析了一定频率的弹性波在层状介质中传播时的位移、应力的模式. 计算结果展现了Hamilton体系和辛几何在弹性波传播问题研究的应用前景.     

Dynamic soil–structure interaction analysis via coupled finite-element–boundary-element method

In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.     

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.