上海深厚饱和覆盖土层的动力耦合地震反应分析

上海的第四纪沉积土层呈水平层状分布，市区的土层厚度约为270－290m，年平均地下水位0.5-0.7m。上海地区上覆深厚、饱和的沉积土层，对土层的地震反应特性具有重要的影响。基于动力耦合理论，将上海深厚饱和覆盖土层视为由固相和液相组成的两相饱和多孔介质，建立了能够反映上海土层深厚、饱和特点的地震反应计算模型。然后应用该模型，以El Centro,Taft，苏南和唐山地震波作为基岩输入加速度，对深度280m的上海深厚覆盖土层进行了动力耦合地震反应计算，并对加速度、振动孔隙水压力和地基震陷结果进行了分析。     

巨厚沉积土夹火山岩场地非线性地震反应特性

地震基岩深度和土体动力本构模型的选取对核岛场地地震效应评价结果的合理性具有重要影响。以拟建某沿海核电厂深度470 m沉积土夹火山岩层场地的3个钻孔剖面为研究对象,采用一维等效线性波传分析（ELA）法、基于Matasovic本构模型和Davidenkov-Chen-Zhao（DCZ）本构模型的一维非线性分析（NLA）法,选取不同剪切波速的5个岩土层作为地震基岩,研究了输入地震动特性、地震基岩深度和土体动力本构模型的选取对巨厚沉积土夹火山岩层场地非线性地震反应特性的影响。结果表明：（1）以浅层硬岩夹层或深部土层作为地震基岩,NLA法计算的5%阻尼比的地表谱加速度SA的短周期部分较之ELA法的计算值大,但两者计算的地表SA谱的长周期部分几乎一致;（2）基于Matasovic模型和DCZ模型的NLA法计算的地表SA谱谱形和峰值加速度随深度的变化趋势基本一致;（3）从NLA法计算的地表峰值加速度和累积绝对速度而言,以剪切波速约2 500 m/s的浅层硬岩夹层作为地震基岩是适宜的。     

远场大地震作用下大尺度深软场地的非线性地震效应分析

基于ABAQUS软件自行研制的并行计算显式算法集群平台,针对苏州城区典型地层剖面,建立了大尺度深软场地的二维精细化非线性有限元分析模型,对人工地震波和大地震远场地震动作用下深软场地的非线性地震效应进行了比较研究。研究结果表明：（1）与人工地震波作用时深软场地的地表峰值加速度放大效应相比,大地震远场地震波作用时的放大效应尤为显著,由于土介质的横向不均匀性及其非线性,使不同地表的峰值加速度放大效应存在明显的变异性。（2）深软场地对周期小于0.3 s的高频地震波均具有显著的滤波效应;大地震远场地震波作用时,深软场地对周期0.85～1.65 s的长周期地震波的放大效应非常显著,但对2.5～7.0 s的长周期地震波呈现出明显的滤波效应。（3）地震动峰值加速度PGA值沿土层深度和横向的分布形态呈现出明显的高低起伏现象,在不同成因的土层更迭面附近及土介质横向不均匀性显著的区域,地震波的局部聚焦放大和过滤减小现象尤为明显,且大地震远场地震动作用时,20 m以浅土层的PGA值呈现出非常显著的放大效应。（4）地震波的频谱特性、土层的横向不均匀性对深软场地地表加速度反应谱? 谱的谱形有显著影响;给出了描述加速度反应谱沿土层深度变化特征的三维谱形曲线,可以直观地展示出深软场地中细长地下结构地震反应可能存在类共振现象的土层深度。     

How the response spectrum of non-liquefied loose-to-medium sand deposits is affected by the groundwater level

This study examines the water level influence on the seismic response of sand deposits. The seismic ground response of a sand deposit for two cases was computed: groundwater level near the surface and at great depth. The ratio of soil factors computed for both cases evaluates the influence of having the water table near the surface. This parameter is small when the full range of periods is considered, but rises to ∼1.15 when narrowing the range of periods to 0.2–0.6 s, corresponding to the plateau of Eurocode 8 ground type C response spectrum. So, the water level effect is of the same magnitude of the existing soil factor in seismic codes.     

First Order Seismic Microzonation of Delhi,India Using Geographic Information System (GIS)

A first order seismic microzonation map of Delhi is prepared using five thematic layers viz., Peak Ground Acceleration (PGA) contour, different soil types at 6 m depth, geology, groundwater fluctuation and bedrock depth, integrated on GIS platform. The integration is performed following a pair-wise comparison of Analytical Hierarchy Process (AHP), wherein each thematic map is assigned weight in the 5-1 scale: depending on its contribution towards the seismic hazard. Following the AHP, the weightage assigned to each theme are: PGA (0.333), soil (0.266), geology (0.20), groundwater (0.133) and bedrock depth (0.066). The thematic vector layers are overlaid and integrated using GIS. On the microzonation theme, the Delhi region has been classified into four broad zones of vulnerability to the seismic hazard. They are very high (> 52%), high (38–52%), moderate (23–38%) and less ( < 23%) zones of seismic hazard. The “very high” seismic hazard zone is observed where the maximum PGA varies from 140 to 210 gal for a finite source model of M_w 8.5 in the central seismic gap. A site amplification study from local and regional earthquakes for Delhi region using Delhi Telemetry Network data shows a steeper site response gradient in the eastern side of the Yamuna fluvial deposits at 1.5 Hz. The ‘high’ seismic hazard zone occupies most of the study area where the PGA value ranges from 90 to 140 gal. The ‘moderate’ seismic hazard zone occurs on either side of the Delhi ridge with PGA value varying from 60 to 90 gal. The ‘less’ seismic hazard zone occurs in small patches distributed along the study area with the PGA value less than 60 gal. Site response studies, PGA distribution and destruction pattern of the Chamoli earthquake greatly corroborate the seismic hazard zones estimated through microzonation on GIS platform and also establishes the methodology incorporated in this study.     

苏州第四纪深厚沉积层一维等效线性和非线性地震反应对比分析

合理描述土体动力本构关系对场地设计地震动参数取值的准确性有显著影响。以苏州城区200 m深的钻孔剖面为研究对象,对比分析了一维等效线性波传分析法（ELA法）和基于修正Matasovic本构模型的多自由度集中质量非线性分析法（NLA法）给出的深厚场地地震反应,研究了基岩输入地震动特性和地震基岩面的选取对深厚场地地震反应的影响。结果表明：（1）随基岩地震动强度（PGA）的增大,ELA法给出的地表PGA呈单调递增的特征,而NLA法给出的地表PGA呈先快速增大后缓慢减小或几乎不变的趋势;（2）ELA法和NLA法给出的地表加速度反应谱（Sa谱）在短周期范围内存在明显差异,ELA法对基岩高频地震动具有显著的滤波效应,而NLA法对基岩高频地震动的影响表现为随基岩PGA的增大先放大后减小的特征;（3）随地震基岩面深度的增大,地表Sa谱的谱值呈现出稍许增大的趋势,但对周期T<0.1 s部分,NLA法给出的地表Sa谱的谱值则呈现出稍许减小的现象;（4）中、大震作用下,地表地震动持时不仅与基岩地震动特性有关,还与地震反应分析方法和地震基岩面的选取密切相关,NLA法能更合理地反映基岩地震动强度和上覆土层厚度对地表地震动持时的影响。     

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

Nonlinear Site Response and Liquefaction Analysis in the New Madrid Seismic Zone

Many existing highway bridges in the New Madrid Seismic Zone are located in the Mississippi Embayment, consisting of deep soil deposits and liquefaction susceptible near surface soils. It is important to understand the comprehensive impact of deep soil deposits and liquefaction on the response of the bridge foundations under seismic loading. A nonlinear soil model is then presented to study the impacts of the deep soil deposit and liquefaction on response analysis. The soil model has the advantage of using input parameters that can be obtained from conventional field and laboratory testing methods, which makes it attractive to engineering practice. The model calibrations used field recorded motions and laboratory test data, which indicate that the model provides an acceptable outcome based on simple input parameters. The model is implemented into the site response analysis for a typical Missouri highway bridge site in this seismic zone. The effect of the deep soil deposit and liquefaction on the site response analyses is discussed.     

深厚土层流固多向地震全耦合动力响应分析

根据上海地区典型土层的力学特性,拟合得到简化亚塑性边界面模型的动力计算参数,同时考虑流、固耦合以及地震多向作用的耦合,以San Fernando1971年地震记录的加速度时程曲线为输入地震波,对上海280 m深典型地层的地震反应进行了分析,得到了有意义的结果。     

上海软土场地三维非线性地震反应分析

考虑土的非线性,采用弹塑性边界面模型,对上海软土场地进行三维地震反应分析。利用多维地震动作用下的水平地层弹塑性动力反应分析程序,以具有水平和垂直三向完整加速度记录的Taft波作为地震输入,对比分析单向和多向地震输入场地水平和竖向地震反应特征,并分析不同强度地震动输入下竖向和水平加速度峰值比特征及地下水位变化对场地土层地震反应特征的影响。计算分析表明,竖向与水平向地震反应特征有较大差异;与水平单向地震输入相比,三向地震输入场地土层放大效应明显增大;地下水位上升对水平向和竖向峰值加速度的放大效应影响差异显著,地下水位上升,地表水平峰值加速度放大效应增大,竖向峰值加速度放大效应减小,研究结果对上海地区的工程抗震设计具有参考作用。     

基于土体小应变本构模型的TRD工法成墙试验数值模拟

等厚度水泥土搅拌墙技术即TRD工法,近年来在深基坑工程中得到了广泛应用。以上海国际金融中心基坑工程开展的0.7 m厚、8 m宽、56.7 m深TRD成墙试验为背景,采用有限元方法,并基于土体小应变本构模型对其成墙过程进行了模拟,得到了土体侧向位移和地表沉降曲线,并与实测数据进行了对比。结果表明,距离墙体5 m处两者的土体侧向位移曲线基本一致,而距离墙体1.4 m处的土体侧向位移在深度大于20 m后的计算结果较实测值偏小;地表沉降在靠近墙体处最大,随着距墙体的距离增大而逐渐减小。最后分析了成墙深度对地表沉降和土体侧向变形的影响,结果表明,深度越深,引起的土体侧向变形和地表沉降也越大。通过不同成墙深度引起的地表沉降归一化曲线可看出,TRD成墙引起的最大地表沉降约为0.05%H（H为成墙深度）,沉降影响区域约为1.8H。     

Characterization of site effects in Montreal,Canada

Recent destructive earthquakes have clearly shown that near-surface geological conditions play a major role in the level of ground shaking in urban areas. In Canada, Montreal is ranked second for seismic risk after Vancouver considering its population and regional seismic hazard. The city is largely built on recent unconsolidated marine and river deposits and most of its infrastructure is old and deteriorated. A seismic risk project that includes a combined methodology for site effects zoning in large cities, using microtremor measurements (H/V method) coupled with 1D numerical modelling (SHAKE91), has been initiated. The experimental approach gives good estimates of the fundamental frequency of soft deposits, while the numerical approach provides good estimates of the soil response in terms of amplification factor related to frequency. Main mechanical properties of soft soils were compiled from various data available, and a sample of input rock motions from real and synthetic earthquakes was used to compute soil response. The influence of marine clays on soil response is significant and is well correlated with thickness of these deposits. PGA amplification factors range from 2 to 4 at frequencies from 2 to 7 Hz, with some occasional larger values. The results demonstrate that the methodology used for our study is both fast and efficient to determine the influence of soft soils in urban environments. Such studies are essential for the effective deployment of seismic instrumentation, land-use planning and seismic mitigation.     

Ramberg-Osgood土动力非线性模型在 ABAQUS软件上的开发及应用

基于ABAQUS有限元软件开发平台,对Ramberg-Osgood土动力非线性本构模型进行了修正。通过合理的假设得出了模型修正参数?、R的取值方法和算法。滞回曲线加卸载规律符合修正后的Masing准则。并进行了三维场地地震动非线性分析,分析了不同输入峰值的地震波对土体应力-应变关系的影响,并与一维场地动力分析Shaking程序的计算结果进行了比较。分析结果表明,非线性分析时的水平位移和相对位移要大于等效线性情况下的相应位移,这表明采用何种本构关系将影响到得到的场地的地震动反应结果,因此,在工程应用时需加以合理的选择,并对最终结果做出合理的判断     

饱和软土层地铁列车运行引起的环境振动研究

我国沿海地区大量的地铁隧道都修建在饱和软土地层,饱和软土地层中地铁列车运行引起的环境振动逐渐成为社会关注的问题。以上海2号线某地铁区间隧道为研究对象,采用循环运动模型本构模型和水土耦合动力分析方法,分析了饱和软土中地铁列车运行引起的地表振动加速度以及振动位移响应规律,采用加速度振级对地表的振动强度进行了评价,分析了单次列车通过时引起的隧道下卧土层超孔隙水压力（EPWP）响应规律。计算结果表明,饱和土的振动响应比干土要小得多,地表水平向与竖向的振动规律明显的不同,超孔隙水压力沿隧道下方沿向是不断减少的,在距离隧道中心的纵断面上是先增加后减少。分析饱和土中地铁运行引起的环境响应规律,可为地铁沿线建筑物减隔振设计及沉降控制提供一定依据。     

Liquefaction hazard mapping at the town of Edessa,Northern Greece

The town of Edessa is located on Northern Greece at a region that is characterized as low seismicity zone due to the fact that few moderate events of M < 6 occurred during the last century. According to the Greek Seismic Code, the expected acceleration having a 10% probability of being exceeded in 50 years is equal to 0.16g. However, an amplification of ground motion is likely to occur due the local geology that is consisted of Holocene fluvio-torrential deposits. The basic aim of this paper is to evaluate the site amplification due to geological conditions and to assess the liquefaction hazard. In order to achieve this, 1-D site response analyses were performed. The data that were employed for the construction of the numerical models have been collected from borings with standard penetrations tests (SPT) that were drilled for construction purposes. Afterward, the liquefaction potential of the subsoil layers was evaluated taking into consideration two seismic scenarios. The first scenario was based on the seismic parameters, earthquake magnitude and PGA, assigned by the Greek Seismic Code. On the second seismic model, we employed the values of acceleration, resulted from the 1-D analyses and the earthquake magnitude as it was defined by the Greek Seismic Code. In order to compile the liquefaction hazard maps, we initially estimated the liquefaction potential index (LPI) of the soil columns using the parameters provided by SPT, for both seismic loadings, and afterward we correlated these values with the proposed classification of the severity of liquefaction-induced deformations. In addition, having computed the value of probability based on the LPI, liquefaction manifestations probability maps were compiled for both scenarios. The result of this study was that liquefaction-induced ground disruptions are likely to occur at the center of the city, among the branches of Voda River, only when the amplified values of acceleration are taken into account to the computation of liquefaction potential.     

基于随机地震反应的桥墩-群桩-土相互作用研究

考虑地震动的随机性,利用复反应分析技术,采用随机地震反应计算方法对某一特大型桥梁群桩基础与土动力相互作用效应进行了数值试验研究。将土与群桩体系视为一个整体进行有限元离散,采用等效线性化方法考虑土体的动力非线性性能。将桥墩-群桩-土相互作用体系与自由场随机地震反应进行比较,结果表明：相互作用效应的影响与桩土模量差异以及土体与群桩基础距离有关;软土层剪应变水平及分布发生了显著变化,群桩基础两侧附近土域剪应变呈现明显的弧形分布;地表及浅层土体最大地震加速度反应有所减小,但覆盖层中下部土体加速度反应峰值明显增大,增幅达5 %～30 %左右;此外,地震地面运动的频谱成分存在显著差别。桥梁桩基础抗震设计中应充分考虑桥梁结构-群桩-土相互作用效应。     

Influence of seismic motions on behavior of piles in liquefied soils

In the present study an analytical procedure based on finite element technique is proposed to investigate the influence of vertical load on deflection and bending moment of a laterally loaded pile embedded in liquefiable soil, subjected to permanent ground displacement. The degradation of subgrade modulus due to soil liquefaction and effect of nonlinearity are also considered. A free headed vertical concrete elastic nonyielding pile with a floating tip subjected to vertical compressive loading, lateral load, and permanent ground displacement due to earthquake motions, in liquefiable soil underlain by nonliquefiable stratum, is considered. The input seismic motions, having varying range of ground motion parameters, considered here include 1989 Loma Gilroy, 1995 Kobe, 2001 Bhuj, and 2011 Sikkim motions. It is calculated that maximum bending moment occurred at the interface of liquefiable and nonliquefiable soil layers and when thickness of liquefiable soil layer is around 60% of total pile length. Maximum bending moment of 1210 kNm and pile head deflection of 110 cm is observed because of 1995 Kobe motion, while 2001 Bhuj and 2011 Sikkim motions amplify the pile head deflection by 14.2 and 14.4 times and bending moment approximately by 4 times, when compared to nonliquefiable soil. Further, the presence of inertial load at the pile head increases bending moment and deflection by approximately 52% when subjected to 1995 Kobe motion. Thus, it is necessary to have a proper assessment of both kinematic and inertial interactions due to free field seismic motions and vertical loads for evaluating pile response in liquefiable soil.     

Reverse fault slip through soft rock and sand strata by centrifuge modeling tests

The devastating damage after the 1999 Chi-Chi and 1999 Izmit earthquakes has greatly motivated soil–reverse fault interaction studies. However, most centrifuge modeling studies have employed a single homogeneous soil layer during testing, which does not represent in situ conditions. Indeed, while geological conditions vary spatially, engineering soils are often underlain by soft rocks. Therefore, four centrifuge models were developed to evaluate the effect of soft rock layers on the ground surface and subsurface deformation. Sand–cement mixtures of varying thicknesses with a uniaxial compressive strength of 0.975 MPa, simulating extremely soft rock, were overlain by pluviated sandy soil. The model thickness was 100 mm, corresponding to 8 m in the prototype scale when spun at 80 g. Every model was subjected to a vertical offset of 50 mm/4 m (0.5 H; H: total sedimentary deposit thickness) along a reverse fault with a 60° dip. The results indicate that the presence of a soft rock stratum results in the creation of a horst profile at the ground surface. Additionally, the thinner the soil layer on top of the soft rock stratum is, the longer and higher the horst created at the ground surface. Consequently, the fault deformation zone lengthens proportionally with the increasing thickness ratio of the soft rock. Furthermore, the presence of soft rock as an intermediary stratum between bedrock and soil causes the deformation zone boundary on the hanging wall side to move in the direction of fault movement.

     

Shear wave propagation along infinite slopes: A theoretically based numerical study

In this paper, the seismic response of ‘infinitely’ long slopes is numerically analysed via the formulation of a 1D analytical/numerical model, in which the soil mechanical behaviour is assumed to be elasto‐perfectly viscoplastic and simple shear (SS) kinematical constraints are imposed. In order to simplify the problem, a theoretically based procedure to set up a fully 1D shear constitutive model is defined, within which the mechanical response of a multiaxial relationship is condensed. The use of a 1D shear constitutive model is aimed at reducing the number of unknowns and, therefore, the computational costs. In particular, the case of the Mohr–Coulomb yield criterion is considered, while an enhanced Taylor–Galerkin finite element algorithm is employed to simulate the seismic wave propagation within the soil stratum. The proposed ‘condensation’/calibration procedure captures both the ‘pseudo’‐hardening pre‐failure behaviour and the influence of dilation on the occurrence of strain‐localization, which characterize, under SS conditions, the static response of virgin perfectly plastic soils. The effectiveness of the conceived method is shown with reference to freshly deposited deposits, while, in the case of highly overconsolidated strata, some difficulties arise because of the brittle behaviour induced both by unloading and non‐associativeness. Copyright © 2011 John Wiley & Sons, Ltd.