考虑Gmax附加衰减的场地地震响应离心机试验数值模拟

在常规的场地地震响应有效应力分析法中,通常只考虑Gmax随有效应力和孔隙比的变化,而最新的试验研究表明,Gmax在大应变幅振动下存在附加衰减。收集了不同类型土体在大幅振动下的298组Gmax数据,拟合并给出了对应于两种Gmax附加衰减模式的计算模型,通过编制可考虑Gmax附加衰减影响的一维场地地震响应有效应力分析程序,分别对典型场地和离心机振动台模型场地的地震响应进行分析对比。结果表明：Gmax附加衰减对场地加速度和剪应力响应均有较大影响,将加速场地液化进程和扩大场地最终液化范围;考虑Gmax附加衰减可较好地模拟动力离心模型试验场地动力响应;但在强震情况下,场地迅速液化,Gmax附加衰减对场地动力响应的影响反而减弱。     

Finite element analyses of two-tier geosynthetic-reinforced soil walls: Comparison involving centrifuge tests and limit equilibrium results

This study presents the procedure and results of the finite element (FE) analyses of a series of centrifuge tests on geosynthetic-reinforced soil (GRS) two-tier wall models with various offset distances. The objectives of this study were to evaluate the applicability of FE for analyzing GRS two-tier walls with various offset distances and to investigate the performance and behavior of GRS two-tier walls in various stress states. The FE simulations were first verified according to the centrifuge test results by comparing the locations of failure surfaces. The FE results were then used to investigate the effective overburden pressure, mobilization and distribution of reinforcement tensile loads, and horizontal deformation at the wall faces. The interaction between two tiers was investigated based on the FE results, which were also used to examine the modeling assumption of reinforcement tensile loads in limit equilibrium (LE) analysis and to evaluate the design methods in current design guidelines. This study demonstrated favorable agreement between FE and the centrifuge model in locating the failure surface. The FE results indicated that as the offset distance increased, the reinforcement tensile load and wall deformation decreased in both the upper and lower tiers, suggesting the attenuation of interaction between the two tiers. The maximum tensile loads of all reinforcement layers at the wall failure predicted using FE analysis and LE method assuming uniform distribution of reinforced tensile loads were comparable. Compared with the FE results, the Federal Highway Administration (FHWA) design guidelines are conservative in determining the effect of overburden pressure, required tensile strength, location of maximum tension line (for designing the reinforcement length), and the critical offset distance. Furthermore, the FHWA design guidelines do not account for the influence of the lower tier on the upper tier that was observed in this study.     

Displacement based seismic assessment of base restrained retaining walls

The present work deals with assessment of earthquake-induced displacement of the base restrained retaining walls (RW’s). A detailed and rigorous finite element (FE) investigation has been carried out following the shaking table experiments on a scaled-down RW model. The FE simulations were performed by conducting several nonlinear time history analyses on a two-dimensional (2D) plane strain FE model of a prototype RW. The hardening and softening of backfill have been simulated by calibrating the Mohr Coulomb material model against the triaxial test results. Role of different backfill into the seismic performance of base restrained RW has also been investigated. It was observed that the cohesionless backfill has a slight influence on the earthquake induced displacement of base restrained RW’s. Amplification of horizontal acceleration in backfill has been observed with no direct correlation with the applied earthquake excitation. The understanding and findings based on shaking table experiment and FE simulations have been used for development of an analytical model for estimation of earthquake induced displacement of base restrained RW. The validity of proposed analytical model has also been examined against the shaking table experiment and FE simulation results.

     

A fully coupled non-linear dynamic analysis procedure and its verification using centrifuge test results

The equations governing the dynamic behavior of saturated porous media as well as a finite element spatial discretization of these equations are summarized. A three-parameter time integration scheme called the Hilber–Hughes–Taylor α-method is used together with a predictor/multi-corrector algorithm, instead of the widely used Newmark's method, to integrate the spatially discrete finite element equations. The new time integration scheme possess quadratic accuracy and desirable numerical damping characteristics. The proposed numerical solution and bounding surface plasticity theory to describe the constitutive behaviour of soil have been implemented as the computer code DYSAC2. Predictions made by DYSAC2 code are verified using dynamic centrifuge test results for a clay embankment. Importance of initial state of a soil on its dynamic behaviour is demonstrated.     

考虑水库岸坡特征的滑坡稳定性计算修正传递系数法

随着低碳能源的推行,大量水电站的兴建引发了大量库岸滑坡问题。库岸滑坡与普通滑坡有着显著区别,但在实际工程的稳定性计算过程中,工程人员常将二者归为一体,采用同样的思路和方法进行研究。为了更好地治理库岸滑坡,根据库岸自身的特点,在长期野外地质调查和室内文献调研的基础上,提出按结构面特征、坡体位置和滑体物质成分进行库岸斜坡分区的原则;针对库岸滑坡的特点,推导了被牵引段有推力+后缘拉裂面有剪力、被牵引段不存在推力+后缘拉裂面有剪力、被牵引段不存在推力+后缘拉裂面无剪力三种情况的水库岸坡滑坡稳定性计算公式;结合工程实例,验证基于水库岸坡特征滑坡稳定性计算修正传递系数法的正确性。     

A performance-based framework for assessing liquefaction potential based on CPT data

This article describes a new performance-based approach for evaluating the return period of seismic soil liquefaction based on standard penetration test (SPT) and cone penetration test (CPT) data. The conventional liquefaction evaluation methods consider a single acceleration level and magnitude and these approaches fail to take into account the uncertainty in earthquake loading. The seismic hazard analysis based on the probabilistic method clearly shows that a particular acceleration value is being contributed by different magnitudes with varying probability. In the new method presented in this article, the entire range of ground shaking and the entire range of earthquake magnitude are considered and the liquefaction return period is evaluated based on the SPT and CPT data. This article explains the performance-based methodology for the liquefaction analysis – starting from probabilistic seismic hazard analysis (PSHA) for the evaluation of seismic hazard and the performance-based method to evaluate the liquefaction return period. A case study has been done for Bangalore, India, based on SPT data and converted CPT values. The comparison of results obtained from both the methods have been presented. In an area of 220 km² in Bangalore city, the site class was assessed based on large number of borehole data and 58 Multi-channel analysis of surface wave survey. Using the site class and peak acceleration at rock depth from PSHA, the peak ground acceleration at the ground surface was estimated using probabilistic approach. The liquefaction analysis was done based on 450 borehole data obtained in the study area. The results of CPT match well with the results obtained from similar analysis with SPT data.     

饱和砂土地层中隧道结构动力离心模型试验

饱和砂土地层中隧道结构可能会因地震地基液化而发生破坏。通过对可液化地层中地铁隧道结构的地震反应进行动力离心模型试验,研究了饱和松砂地基在地震作用下的反应特性、可液化地层中地铁隧道结构的上浮及变形特性和设置截断墙对限制隧道结构上浮的效果等问题。研究结果表明,地基液化引起的隧道衬砌上的附加变形内力以及隧道上浮量主要受地基液化时土水压力的变化影响。截断墙的设置限制了隧道两侧土体向隧道下方流动的趋势,有效减小了隧道结构的上浮量     

基于饱水页岩地震频段岩石物理测试的地震响应频散特征分析

页岩储层由于骨架和孔隙结构的微观非均质性,能够在地震频段内引起波诱导的孔隙流体流动相关的弹性频散和衰减,但声波测井和实验室超声实验测量频率段远超地震频段,不可避免产生测量速度和衰减差异。这里基于澳大利亚惠灵顿地区的页岩岩心在干燥和水饱和条件下地震频段的弹性和衰减测量结果,考虑频变反透射系数和衰减对地震波反射的影响,针对不同层厚的页岩储层开展地震响应计算,并利用Wigner-Ville分布时频分析技术讨论了页岩频散、衰减性质及层厚对其地震响应规律的作用。对于无明显速度频散和衰减的页岩储层,顶、底反射振幅的频变属性受调谐效应控制,主要呈现出低频异常和陷频特征;而对于速度频散和衰减强烈的储层,顶、底反射振幅的频变属性同时受到调谐效应与储层本身频散衰减性质的作用。干燥和含水页岩储层地震响应的频变性质差异有助于加深对含流体页岩地震响应特征的理解,而基于地震频段实测数据正演和时频分析技术,对于提高页岩的地震勘探精度和流体识别具有重要的作用。     

基于能量桩现场试验的土体综合热导率系数研究

能量桩技术兼具支撑上部荷载和浅层地热能换热器双重功能;作为一种节能减排技术,近年来在国内外获得了一定的发展。然而,目前简单套用基于传统地埋管换热器获得的土体综合热导率系数,无法准确计算能量桩换热效率。依托河南理工大学某低承台3×3能量桩群桩工程应用,开展基于能量桩的土体综合热导率系数测试现场试验和数值模拟研究,分析加热时长、加热功率、流速及桩长等因素对土体综合热导率系数的影响规律,继而探讨能量桩在群桩中的布置形式对土体综合热导率系数的影响规律。研究结果表明,基于传统地源热泵测试所发展起来的土体综合热导率系数线热源分析方法,并不适用于分析基于能量桩现场实测所获得的相关数据;有必要推导一套考虑桩径影响的、适用于能量桩的土体综合热导率系数测试与计算分析方法。     

地震加速度动态分布及对高土石坝坝坡抗震稳定的影响

《水工建筑物抗震设计规范》（DL5073-1997）建议的土石坝地震加速度动态分布系数适合于150 m以下的土石坝,而前许多土石坝的设计高度已远大于150 m。与低坝相比,高坝受地基刚性的约束减弱,坝体的自振周期延长。在坝体地震反应中,高阶自振周期与地震卓越周期耦合的机率增大,高阶振型的振动易被激发放大,从而导致坝体地震加速度沿坝高的分布与低坝有所不同。据此,采用有限元法研究了高土石坝的加速度分布,提出了高于150 m土石坝的地震加速度动态分布系数图示。在此基础上,利用堆石材料的非线性强度准则,对高度超过150 m的土石坝坝坡抗震稳定性作了进一步的分析,得出随着坝体地震加速度动态分布系数的降低,坝坡临界破坏的安全系数有所提高。     

Implementation of a hydromechanical elastoplastic constitutive model for fully coupled dynamic analysis of unsaturated soils and its validation using centrifuge test results

Prediction of unsaturated soil behavior during earthquake loading has received increasing attention in geotechnical engineering research and practice in recent years. Development of a fully coupled analysis procedure incorporating a coupled hydromechanical elastoplastic constitutive model for dynamic analysis of unsaturated soils has, however, been limited. This paper presents the implementation of a coupled hydromechanical elastoplastic constitutive model into a fully coupled dynamic analysis procedure and its validation using a centrifuge test. First, the fully coupled finite element equations governing the dynamic behavior of unsaturated soils with the solid skeleton displacement, pore water pressure, and pore air pressure as nodal unknowns are briefly presented. The closest point projection method is then utilized to implement the coupled hydromechanical elastoplastic constitutive model into the finite element equations. The constitutive model includes hysteresis in soil–water characteristic curves, cyclic elastoplasticity of the solid skeleton, and the coupling mechanisms between the SWCCs and the solid skeleton. Finally, the analysis procedure is validated using the results from a dynamic centrifuge test on an embankment constructed of compacted unsaturated silt subjected to base shaking. Reasonable comparisons between the predicted and measured accelerations, settlements, and deformed shapes are obtained.

     

Insight into the seismic response of earth dams with an emphasis on seismic coefficient estimation

This paper presents a numerical investigation of the seismic response of earth dams by employing results from 110 nonlinear two-dimensional (2D) dynamic analyses of four different cross-sections with heights ranging from 20 to 120 m. The analyses were of a parametric nature, considering the effects of seismic excitation characteristics (intensity and frequency content), foundation soil stiffness, and the existence of typical stabilising berms and/or an impounded reservoir. The results of these studies indicated that the predominant period of a dam’s vibration was strongly affected by its height and the input motion characteristics. The results also indicated that the peak acceleration at the dam’s crest was governed by its height, the input motion characteristics, and the stiffness of the foundation soil, but not by the other parameters. These same analyses yielded results on pseudo-static seismic coefficients for a total of 1084 potential sliding masses within the analysed cross-sections, demonstrating that the seismic coefficients decreased as the sliding mass became deeper and bulkier, increased if the mass was located upstream rather than downstream, and were strongly affected by the seismic excitation characteristics and stiffness of the foundation soil. Moreover, these results allowed for a thorough evaluation of existing methodologies for seismic coefficient estimation, quantifying their accuracy and depicting their limitations. This evaluation process also illustrated the fact that there is currently no methodology accounting for all significant problem parameters.     

Exact seismic response of smooth rigid retaining walls resting on stiff soil

The assessment of forces exerted on walls by the backfill is a recurrent problem in geotechnical engineering, owing to its relevance for both retaining systems and underground structures. In particular, the work by Arias and colleagues, and later also the one by Veletsos and Younan, among others, becomes pertinent when considering pressure increments on underground structures triggered by seismic events. As a first step, they studied the response of a rigid retaining wall resting on rigid bedrock subjected to SV waves, introducing some simplifying assumptions. This paper presents the exact solution to this reference problem. The solution is given in horizontal wavenumber domain; hence, it comes in terms of inverse Fourier transforms, which can be approximated numerically in Mathematica, which in turn are verified against finite-element simulations. Specific features of this exact solution that were not captured by prior engineering approximations are highlighted and discussed.     

Finite element simulations of seismic effects on retaining walls with liquefiable backfills

Finite element simulations of two centrifuge tests on the same cantilever retaining wall model holding liquefiable backfill were conducted using the Biot formulation‐based program DIANA–SWANDYNE II. To demonstrate the effects due to different pore fluids in seismic centrifuge experiments, water was used as the pore fluid in one experiment whereas a substitute pore fluid was used in the second experiment. The cantilever wall model parameters were determined by comparing simulations with measurements from free‐vibration tests performed on the model wall without backfill. The initial stress conditions for dynamic analysis for the soil backfill were obtained by simulating static loads on the retaining wall from the soil backfill. Level‐ground centrifuge model results were used to select the parameters of the Pastor–Zienkiewicz mark III constitutive model used in the dynamic simulations of the soil. The effects due to different pore fluids were captured well by the simulations. The magnitudes of excess pore pressures in the soil, lateral thrust and its line of action on the wall, and wall bending strains, deflections, and accelerations were predicted well. Predictions of settlements and accelerations in the backfill were less satisfactory. Relatively high levels of Rayleigh damping were needed to be used in the retaining wall simulations in order to obtain numerically stable results, which is one of the shortcomings of the model. The procedure may be used for engineering purpose dealing with seismic analysis of flexible retaining walls where lateral pressures, bending strains and deflections in the wall are typically of importance. Copyright © 2008 John Wiley & Sons, Ltd.     

基于动力离心试验的软基尾矿库地震响应研究

基于离心机振动台,分别针对典型软基尾黏土尾矿库及加高扩容后的尾矿库开展了动力离心模型试验,重点探究了加速度分布规律、软土及尾矿内部孔压变化规律以及“软基?库内尾矿?尾矿坝”系统的变形规律等内容。试验结果表明：软基对地震动的放大效应较为微弱,而坝体加速度沿高程逐渐放大,高层子坝加速度响应最为强烈,加高扩容后的现坝顶加速度响应可达原坝顶的2.2倍;地震作用下软土与库内尾矿内部均会产生一定的孔压增量,但未达到液化状态。地震动强度、软基及库内尾矿的固结状态对尾矿库的变形模式影响较大。当固结不充分时,在强震作用下易发生尾矿的水平滑移,进而造成坝顶及下游软土隆起。在固结较为充分且地震动强度较弱的情况下,变形模式以震陷为主。该试验成果将为此类尾矿库的动力稳定分析及抗震加固设计提供一定参考。     

确定含水层渗透系数的冲击试验方法

运用Hvorslevr冲击试验方法测定CIRP试验场地下研究设施(URF)试验含水层的渗透系数。由于含水层厚度小、富水性差,正在进行核素迁移试验,以及操作空间有限,因此不能进行抽水或注水试验,同时考虑到试验范围较小,故采取冲击试验求含水层的渗透系数。含水层介质为亚粘土,隔水底板距URF底板面约6m。试验的冲击量小,可以忽略URF内南部含水层对试验结果的影响。试验钻孔的孔径和井管内径相近,减小了钻孔结构对试验结果的影响。试验结果(472×10-4、658×10-4cm/s)与采用自URF内试验含水层的未扰动土柱所作的室内渗透试验结果(41×10-4～82×10-4cm/s)在同一数量级。     

利用HCA研究黏土渗透系数的标准探索

黏性土渗透系数影响因素较多,甚至不同试验手段对试验结果都有较大影响。针对目前渗流试验中存在的问题,提出采用空心圆柱扭剪仪（HCA）测定软土渗透系数的方法,对其开展意义、可行性及试验步骤进行说明分析,并以重塑高岭土作为研究对象,开展一系列试验研究,探讨利用HCA开展渗流试验的标准,最终进行动力特性对软土渗透性影响的试验初步探索。结果表明：（1）以3 600 s作为单位时间,当连续24 h的单位时间内平均渗透系数差值小于2%,将该阶段中得到渗透系数平均值作为最终渗透系数较为准确;（2）综合试验效率、试验间的可比性以及试验结果准确性,建议利用HCA开展重塑高岭土渗流试验时可选水力梯度范围为10～15;（3）重塑高岭土经过交通荷载作用10万次后,其渗透系数比经静载压缩同等轴向应变后的渗透系数低12.79%,即研究动力特性对渗透系数的影响意义重大。     

不同运动模式的悬臂式挡墙地震永久位移算法EI北大核心CSCD

为了有效确定支挡边坡的悬臂式挡墙的地震永久位移,考虑悬臂墙踵板上方局部填土可能存在的不同滑裂特征,基于拟静力法、塑性极限分析上限定理与Newmark滑块法,针对地震条件下可能发生的墙-坡整体对数螺旋面转动、墙-局部填土体系水平滑动及绕墙趾转动3种运动模式,根据较为严格的力学定义分别推导了墙体地震永久位移计算公式。实例分析表明,前两种运动模式所得的地震永久位移相近,均远大于后一种模式,在工程设计中属于位移控制模式。与既有的经验公式法对比表明,文中方法与其误差在30%以内,且比概率置信水平取0.7的Ambraseys-Menu方法小7%。对于墙-坡系统整体旋转滑动模式,在水平地震影响系数一定的情况下,10个主要参数敏感性正交分析结果表明,其对水平屈服加速度影响的敏感性大小排序为:填土黏聚力、墙高、填土内摩擦角、墙体重度、立臂顶宽、填土重度、趾板长度、竖向地震影响系数、底板厚度和踵板长度。     

Regionalization of seismic hazard in Greece based on seismic sources

A semi-probabilistic approach to the seismic hazard assessment of Greece is presented. For this reason, a recent seismotectonic model for shallow and intermediate depth earthquake sources, based on historical as well as on instrumental data, was used. Different attenuation formulae were proposed for the macroseismic intensity and the strong ground motion parameters for the shallow and the intermediate focal depth shocks. The data were elaborated in terms of McGuire's computer program, which is based on the Cornell's method.A grid of equally spaced points at 20 km distance was made and the seismic hazard recurrence curves for various parameters of the seismic intensity was estimated for each point. Finally, seismic hazard maps for the area of Greece were compiled utilizing the entire range of recurrence curves. These maps depict areas of equal seismic hazard and for every area the analytical relations of the typeSI =f(Tm), whereSI is a seismic intensity parameter andTm is the mean return period, were determined.     

黄土路基震陷系数预测的模糊神经网络方法

以郑（州）-西（安）高速铁路为例,采集黄土样品,进行室内土工试验及动三轴震陷实验,以影响震陷系数主要因素为基础,建立了高速铁路黄土路基震陷系数的模糊神经网络模型,通过对样本的训练和预测,表明该模型预测的结果与动三轴试验震陷结果比较接近,是一种比较理想的预测方法。