Fully fluid-solid coupling dynamic model for seismic response of underground structures in saturated soils

The seismic response characteristics of underground structures in saturated soils are investigated. A fully fluid-solid coupling dynamic model is developed and implemented into ABAQUS with a user-defined element to simulate the dynamic behavior of saturated soils. The accuracy of the model is validated using a classic example in literature. The performance of the model is verified by its application on simulating the seismic response characteristics of a subway station built in saturated soils. The merits of the model are demonstrated by comparing the difference of the seismic response of an underground structure in saturated soils between using the fully coupling model and a single-phase medium model. The study finds that the fully coupling model developed herein can simulate the dynamic response characteristics of the underground structures in saturated soils with high accuracy. The seismic response of the underground structure tends to be underestimated by using the single-phase medium model compared with using the fully coupling model, which provides a weaker confining action to the underground structure.     

饱和砂土地铁车站的振动台试验研究

地震作用下,饱和砂土地层地铁车站的动力反应特征是城市轨道工程抗震的关键问题。以太原地铁新近沉积粉细砂地层地铁工程为对象,通过模拟地震运动输入的饱和砂土地基地下结构的振动台模型试验,分析了不同峰值加速度地震作用下饱和砂土与地下结构相互作用的动力反应性状。研究了地震波作用的放大效应与频率特征,动孔压比增长发展过程和液化区域分布,以及动土压力的变化规律。表明加速度放大系数为1.5~2.0;0.1~0.25g峰值加速度地震作用下饱和砂土均产生动孔隙水压力累计发展;0.3g峰值加速度地震作用下饱和砂土产生液化,抑制了土与地下结构的振动放大效应,地表面大量冒水,结构模型出现了明显上浮,地下结构两侧产生震陷。     

Evaluation of the effect of earthquake frequency content on seismic behavior of cantilever retaining wall including soil–structure interaction

A three-dimensional backfill–structure–soil/foundation interaction phenomenon is simulated using the finite element method in order to analyze the dynamic behavior of cantilever retaining wall subjected to different ground motions. Effects of both earthquake frequency content and soil–structure interaction are evaluated by using five different seismic motions and six different soil types. The study mainly consists of three parts. In the first part, following a brief review of the problem, the finite element model with viscous boundary is proposed under fixed-base condition. In the second part, analytical formulations are presented by using modal analysis technique to provide the finite element model verification, and reasonable agreement is found between numerical and analytical results. Finally, the method is extended to further investigate parametrically the effects of not only earthquake frequency content but also soil/foundation interaction, and nonlinear time history analyzes are carried out. By means of changing the soil properties, some comparisons are made on lateral displacements and stress responses under different ground motions. It is concluded that the dynamic response of the cantilever wall is highly sensitive to frequency characteristics of the earthquake record and soil–structure interaction.     

平面P-SV波入射时非均匀饱和土 自由场地的响应

基于Biot多孔介质波动模型,研究了非均匀饱和土层对平面P-SV波入射时的动力响应.考虑饱和土地基的物理力学特性沿厚度方向连续变化,利用亥姆霍兹矢量分解原理和动力刚度法,分析了平面入射P-SV波在非均匀饱和土层中的反射和透射,并给出了基岩表面和自由表面处反射系数和透射系数的计算表达式.基于理论推导结果,数值分析了平面SV波入射下非均匀饱和土自由场地的动力响应,其中假设饱和土地基的物理力学性质沿土层深度按幂律梯度变化.数值结果表明,平面SV波入射所引起的地面位移与基岩位移之比均随土层厚度和土体的非均匀程度、波的入射角和入射频率的增加而减小,且其竖向位移比的减小更为显著,厚土层对地震波的耗散作用尤为明显.      

Three-dimensional finite element analysis of the seismic behavior of inclined micropiles

This paper presents a thorough study of the behavior of inclined micropiles under seismic loading. Analysis is carried out using a full three-dimensional finite element modeling. The soil media is assumed to be elastic with Rayleigh damping, while micropiles are modeled as 3D elastic beam elements. The structure is described by a single degree of freedom system composed of a concentrated mass and a column. The paper is composed of four parts. The first part includes a literature survey on the behavior of inclined micropiles. The second part presents the numerical model used in this study. The third part concerns analysis related to the influence of micropiles inclination on the seismic behavior of a group of micropiles embedded in a homogeneous soil with a uniform stiffness. The last part deals with the seismic behavior of inclined micropiles embedded in a soil layer with a depth-based increasing stiffness. The results of this study provide valuable information about the influence of micropiles inclination on dynamic amplification and on the seismic-induced internal forces in micropiles.     

Adaptive element free Galerkin method applied to analysis of earthquake induced liquefaction

An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils.The method is based on the element free Galerkin method (EFGM) and the fission procedure that is part of h-refinement,indicated by error estimation. In the proposed method, a posteriori error estimate procedure that depends on the energy normof stress and the T-Belytschko (TB) stress recovery scheme is incorporated. The effective cyclic elasto-plastic constitutivemodel is used to describe the nonlinear behavior of the saturated soil. The governing equations are established by u-pformulation. The proposed method can effectively avoid the volumetric locking due to large deformation that usually occursin numerical computations using the finite element method (FEM). The efficiency of the proposed method is demonstratedby evaluating the seismic response of an embankment and comparing it to results obtained through FEM. It is shown that theproposed method provides an accurate seismic analysis of saturated soil that includes the effects of liquefaction     

土、结构特性对饱和土-地下综合管廊动力相互作用影响分析

将土体视为固-液两相介质,基于饱和土体有效应力原理,建立饱和土体-地下综合管廊结构体系相互作用动力模型：在地应力平衡的静力状态下,采用Duncan-Chang非线性弹性本构模型,在地震波作用的动力状态下,采用Davidenkov非线性黏弹性本构模型;考虑饱和土体黏弹性动力人工边界条件,并将地震动作用转化为作用在人工边界节点上的动力荷载。模型考察不同土体材料、结构特性以及土-结构接触摩擦对结构地震响应的影响,得出如下结论：（1）地震波的卓越周期与场地卓越周期相近时,引起结构上的变形最大;（2）综合管廊结构管廊壁厚越薄,埋深越深,结构尺寸越大,结构刚度越小,结构变形越大;（3）不考虑土-结构接触面的状态非线性将会增大结构变形。     

饱和土体-地下综合管廊结构地震响应分析

将土体视为固-液两相介质,基于饱和土体有效应力原理,建立饱和土体-地下综合管廊结构体系相互作用动力模型:在地应力平衡的静力状态下采用Duncan-Chang非线性弹性本构模型,在地震波作用的动力状态下采用Davidenkov非线性黏弹性本构模型;考虑饱和土体黏弹性动力人工边界条件,将地震动作用转化为作用在人工边界节点上的动力荷载。模型考察不同地震波时程、地震波加速度峰值、入射角度、孔隙率以及地应力场的影响,得出如下结论:(1)地震波的卓越周期与场地卓越周期相近时引起结构上的变形最大;随着地震波加速度峰值的增大结构变形增大;随着地震波入射角度的增加结构变形增大,地震波斜入射情况下产生的行波效应使得结构变形最大。(2)土体材料的孔隙水压力是影响地震中结构变形的主要因素之一。(3)将土体材料考虑为单相介质时结构上的变形要比考虑为固-液两相介质时大得多,直接将饱和土体场地中得到的地震波等效荷载施加到单相土介质-结构动力相互作用模型上,能够得到与完全基于有效应力法一致的结果。     

基于流固耦合动力模型的饱和土体-隧道体系地震反应研究

基于ABAQUS软件平台,应用自行开发的流固耦合动力模型孔压单元模拟场地土体,并通过黏弹性人工边界方法实现地震动的输入,对饱和土体场地中的双孔隧道结构在地震荷载作用下的动力反应进行研究。计算结果表明:在地震反应结束时刻,场地土体位移幅值在两隧道之间以及两隧道的附近区域较大,而远离隧道的区域则较小;场地底部区域土体的孔压幅值较大,而场地顶部区域土体则较小;隧道左右两侧拱腰部位的衬砌的应力较大,而拱顶部位则较小。计算结果同时表明了流固耦合动力模型孔压单元在饱和土体-隧道体系地震反应研究中的适用性。     

饱和砂土中箱型隧道的地震反应分析

饱和砂土地基在地震作用下存在液化的潜在危险,液化引起的地基失效可能会导致地下结构的严重震害。以Opensees作为计算平台,对饱和砂土中带中柱箱型隧道的地震反应进行输入不同幅值地震动时的动力数值计算,研究场地和结构的加速度反应及其频谱特性、场地的永久变形、隧道的震后位移以及隧道的内力分布。计算结果表明,饱和砂土中箱型隧道的地震附加内力仍受周围土体的相对位移控制,此外在震后隧道可能会产生侧移和上浮的永久位移,并且可能存在残余内力。     

饱和土体瞬态响应有限元分析

给出基于Biot多孔介质理论分析饱和土体在动载荷作用下瞬态响应的有限元公式,数值计算部分采用本文有限元法分别计算一维饱和土柱在两种不同类型动载荷作用下的瞬态响应,并将数值计算结果与文献中的解析解进行比较,二者结果十分吻合,从而验证本文方法的可行性。     

A simplified 3 D.O.F. model of A FEM model for seismic analysis of a silo containing elastic material accounting for soil–structure interaction

The purpose of this study is the evaluation of dynamic behavior induced by seismic activity on a silo system, containing bulk material, with a soil foundation. The interaction effects between the silo and bulk material, as well as the effects produced between the foundation of the silo and the soil, were taken into account. Proposed simplified approximation, as well as the finite model, were used for analysis. The results, from the presented approximation, were compared with a more rigorous obtainment method. Initially, the produced simplified approximation, with elastic material assumption for the grain, could determine the pressures on the dynamic material along with displacements along the height of the silo wall and base shear force, etc., with remarkable precision. Some comparisons, via a change of soil and/or foundation conditions, were also made regarding the seismic pressure of the dynamic material pressure, displacement and base shear forces for both squat and slender silos. Comparing the analytical predictions to results from the numerical simulations produced good results. It can be concluded that the model can be used effectively to perform a broad suite of parametric studies, not only at the design stage but also as a reliable tool for predicting system behavior under the limit state of the system. The results and comprehensive analysis show that displacement effects and base shear forces generally decreased when soil was softer; however, soil structure interaction (SSI) did not have any considerable effects on squat silos and therefore need not be taken into practice.     

遮帘式板桩码头地基地震液化破坏机理

遮帘式板桩码头作为一种新型的板桩结构型式,其抗震性能研究是设计建造过程中的重要环节。在FEM-FDM水土耦合计算的平台上引入循环弹塑性本构模型,借助FORTRAN编程软件形成饱和砂土动力液化分析的数值方法,可有效模拟饱和砂土在地震动力作用下的非线性及大变形特性,同时也可模拟砂土液化流动对遮帘桩和前墙的动土压力。研究表明:地震作用下可液化土层超孔隙水压力比增长并发生较大的水平流动变形,对前墙的水平破坏大于竖向破坏;前墙剪力最大值位于海床与前墙交界处;遮帘桩剪力最大值位移与前墙底平行的位置;后拉杆拉力逐渐变大,前拉杆拉力逐渐变小。通过对板桩码头地震液化灾害的分析,可为抗震和抗液化设计提供参考依据。     

Nonlinear 3D finite element analysis of soil–pile–structure interaction system subjected to horizontal earthquake excitation

The dynamic response of a seismic soil–pile–structure interaction (SSPSI) system is investigated in this paper by conducting nonlinear 3D finite element numerical simulations. Nonlinear behaviors such as non-reflecting boundary condition and soil–pile–structure interaction modeled by the penalty method have been taken into account. An equivalent linear model developed from the ground response analysis and the modified Drucker–Prager model are separately used for soil ground. A comparison of the two models shows that the equivalent linear soil model results in an underestimated acceleration response of the structure under this ground shaking and the soil behavior should be considered as a fully-nonlinear constitutive model in the design process of the SSPSI system. It was also observed that the dynamic response of the system is greatly affected by the nonlinearity of soil–pile interface and is not sensitive to the dilation angle of the soil. Furthermore, the effect of the presence of pile foundations on SSPSI response is also analyzed and discussed.     

基于Comsol的饱和多孔介质动力方程的数值模拟及应用

两相饱和多孔介质的动力响应问题在地震工程领域具有重要的研究意义,由于涉及到固相和液相的动力耦合,使得该问题的求解尤为复杂。本文利用Comsol在求解多场耦合问题上的优点,针对Biot饱和多孔介质u-U耦合形式下的波动方程特征,经过一系列微分算子运算和矩阵变换得到导数形式下的波动方程,基于Comsol Multiphysics提供的广义偏微分方程模式对变形后的波动方程进行求解,并把改进后的无限元边界应用到无限域动力问题的模拟中。通过与饱和多孔介质动力响应的解析解进行对比,验证模型求解技术的可行性和正确性,并在此基础上讨论饱和土地基中空沟隔振效果与饱和土体参数孔隙率、泊松比的关系。通过研究分析,可以为饱和土地基中空沟隔振设计提供一些有价值的参考。     

A frequency‐dependent and intensity‐dependent macroelement for reduced order seismic analysis of soil‐structure interacting systems

The computational demand of the soil‐structure interaction analysis for the design and assessment of structures, as well as for the evaluation of their life‐cycle cost and risk exposure, has led the civil engineering community to the development of a variety of methods toward the model order reduction of the coupled soil‐structure dynamic system in earthquake regions. Different approaches have been proposed in the past as computationally efficient alternatives to the conventional finite element model simulation of the complete soil‐structure domain, such as the nonlinear lumped spring, the macroelement method, and the substructure partition method. Yet no approach was capable of capturing simultaneously the frequency‐dependent dynamic properties along with the nonlinear behavior of the condensed segment of the overall soil‐structure system under strong earthquake ground motion, thus generating an imbalance between the modeling refinement achieved for the soil and the structure. To this end, a dual frequency‐dependent and intensity‐dependent expansion of the lumped parameter modeling method is proposed in the current paper, materialized through a multiobjective algorithm, capable of closely approximating the behavior of the nonlinear dynamic system of the condensed segment. This is essentially the extension of an established methodology, also developed by the authors, in the inelastic domain. The efficiency of the proposed methodology is validated for the case of a bridge foundation system, wherein the seismic response is comparatively assessed for both the proposed method and the detailed finite element model. The above expansion is deemed a computationally efficient and reliable method for simultaneously considering the frequency and amplitude dependence of soil‐foundation systems in the framework of nonlinear seismic analysis of soil‐structure interaction systems.     

A practical numerical method for large strain liquefaction analysis of saturated soils

Accurate prediction of the liquefaction of saturated soils is based on strong coupling between the pore fluid phase and soil skeleton. A practical numerical method for large strain dynamic analysis of saturated soils is presented. The u–p formulation is used for the governing equations that describe the coupled problem in terms of soil skeleton displacement and excess pore pressure. A mixed finite element and finite difference scheme related to large strain analysis of saturated soils based on the updated Lagrangian method is given. The equilibrium equation of fluid-saturated soils is spatially discretized by the finite element method, whereas terms associated with excess pore pressure in the continuity equation are spatially discretized by the finite difference method. An effective cyclic elasto-plastic constitutive model is adopted to simulate the non-linear behavior of saturated soils under dynamic loading. Several numerical examples that include a saturated soil column and caisson-type quay wall are presented to verify the accuracy of the method and its usefulness and applicability to solutions of large strain liquefaction analysis of saturated soils in practical problems.     

地下RC结构物地震响应特征土工离心试验的模拟

地下结构物抗震稳定性的研究具有重要实际意义。基于动态土工离心试验,应用有效应力法分析了不同水平地震作用下土层和地下RC结构物的最大水平位移、不同埋深的土层加速度和超孔隙水压力变化规律,以及结构物的破坏发展特征,并进一步与试验结果进行了相互验证和对比,取得了较好效果,为重要地下结构物的抗震设计提供了依据和参考。     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

可液化土层隧道围岩的动力响应及液化区发展分析

隧道可液化土层围岩对地震动作用非常敏感,可液化土层动孔压的产生和发展使得地下结构受到上浮作用,从而影响地下结构的稳定性.通过对可液化土层中隧道动力响应计算,研究了不同静应力场隧道围岩动孔压场分布、围岩液化区域分布以及衬砌结构仰拱底与拱顶的动孔压差变化.研究结果表明,不同静应力场对围岩可液化土的动孔压分布、液化区域分布及...