刚度和阻尼频率相关的等效线性化方法

在实际工程中发现,在强地震动作用下等效线性化方法低估了地表地震动的峰值加速度反应。通过与场地效应竖向台阵观测记录的对比分析证实了这一现象,其主要原因是等效剪应变水平的取值高于高频段实际剪应变水平,从而抑制了土层地震反应中的高频成份。若干分析表明,动力反应过程中土体的剪应变是随频率而变化的,且试验表明土体的刚度和阻尼是与剪应变相关的,据此假设土体刚度和阻尼是与频率相关的。基于这一假设,对传统等效剪应变的选取如何使得土层对基底地震动输入中的高频成分产生滤波效应的原理进行探讨,且通过频率相关的等效线性化方法对该原理进行证明。将该方法应用于Port Island场地效应竖向台阵的分析,并与实际地震动记录进行了对比,结果表明,较之传统等效线性化方法,能更加合理地考虑地震动输入中的高频成分。     

基于子结构法的成层场地中沉积河谷的散射分析

土层场地对地震动的影响明显,研究成层场地局部复杂地形的波动场特性,对于抗震设防具有重要意义。基于土-结构相互作用理论,含有局部复杂地形的成层场地,可分解成广义结构（即近场复杂地形及其周围土体）和具有规则开挖边界的成层无限地基（即远场）。远场的格林函数可通过精细积分求解频域-波数域内的对偶波动方程获得,利用傅里叶逆变换得到频域-空间域内的柔度矩阵。近场可采用比例边界有限元进行模拟,通过高性能连分式的传递边界求解动刚度。自由场波动可通过传递矩阵法求得,将动刚度和自由场波动位移代入控制方程即可得到散射场的动力响应。数值算例验证了方法的准确性,并利用提出的方法讨论含有软夹层场地局部复杂场地的波动场特性。     

固结比对砂土动剪切模量及地表反应谱的影响

采用共振柱试验方法,研究了固结比对砂土最大动剪切模量以及对动剪切模量比与剪应变非线性关系的影响。通过土层地震反应计算,初步给出了这种影响对地表反应谱作用的估计。对福建标准砂和哈尔滨砂进行了系统试验,提出了计算不同固结比砂土最大动剪切模量的建议公式,公式显示固结比kc>1、最大动剪切模量较kc=1时的增长程度明显比Hardin公式要大。同时,给出了固结比对动剪切模量比与剪应变非线性关系及其对地表反应谱的影响。结果表明,固结比对砂土动剪切模量比与剪应变非线性关系有一定影响。随着固结比增加,模量比非线性关系曲线有一定提高,在强地震动作用下的土层反应分析中应考虑实际固结比的影响。     

纵向非一致激励下自由场土体的非线性地震反应研究

土体地震反应分析对结构的地震反应和抗震安全性评价具有重要的意义,以大型振动台模型试验为手段,研究一致激励和非一致激励下自由场土体的非线性地震反应规律及其影响因素,能够为结构的地震破坏机制分析提供支撑。通过对场地土体的动力特性、加速度响应、剪应力-剪应变曲线以及土体沉降的分析,研究了不同地震动、不同地震强度作用下自由场土体的地震动反应特性及其变化规律。研究结果表明:土体的非线性发展程度不仅与地震动记录有关,还与地震动输入方式、加载等级有关;在大震或纵向非一致激励作用下,土体运动的不一致性导致了更强的土体结构性变化,频率降低,阻尼比增大,土体刚度弱化,土体的非线性发展相对较快;土体剪应力-剪应变曲线和沉降曲线的变化规律一定程度上也反映了土体的塑性发展情况。所得结论和自由场地震反应宏观现象一致,彼此佐证了结论的合理性。     

考虑土体动力特征参数相关性的工程场地随机地震反应分析

研究等效线性化方法中土体动力参数不确定性对场地地震反应的影响,提出了可考虑土体动力特征参数间相关性的土体动剪切模量和动阻尼曲线随机样本生成方法,基于Matlab开发一维等效线性化场地地震反应分析计算程序,用于开展场地地震反应随机动力分析。以Ⅱ类成层场地为例建立一维自由场分析模型,采用不同地震设防水准下露头基岩场地加速度反应谱为目标反应谱,合成人工地震记录折半后作为场地下卧基岩处的输入。计算结果表明,考虑土体动力特征参数不确定性对场地地震反应具有较大的影响,且影响程度与地震动强度、频谱成分及场地基本周期均密切相关;场地最大峰值加速度和地表加速度反应谱的波动范围随地震动强度等级的增加而变大,且最大峰值应变和最大峰值加速度的波动范围可达10%和14%;目标反应谱平台段及场地基本周期处对应的场地地表加速度反应谱的波动范围超过20%。     

土体动非线性黏弹性模型及其ABAQUS软件的实现

基于土的动应力-应变关系--Davidenkov骨架曲线,采用破坏剪应变幅上限值作为分界点,对Davidenkov骨架曲线进行了修正,即当剪应力值大于破坏剪应力值时,土体产生破坏,土体的动剪切模量采用破坏后的动剪切模量;参考Mashing法则,构造了修正后的Davidenkov骨架曲线土体加卸载对应的应力-应变滞回圈曲线。基于ABAQUS软件的操作平台,开发了土体动弹塑性本构模型的子程序。选择南京某个典型软弱场地为研究对象,首先输入余弦加速度时程曲线,然后选取一条强震加速度记录作为基岩输入地震动,对不同峰值加速度水平下该场地的地震反应进行了二维有限元非线性分析,计算结果验证了所建模型的准确性和可用性。     

基于自由场大型振动台试验的饱和砂土固-液相变特征研究

基于饱和砂土自由场地大型振动台试验,利用加速度记录的线性方法反演得到模型地基土动剪应力-动剪应变,并根据流体动力学相关理论,引入非牛顿流体中的表观黏度、零剪切黏度等概念,对场地液化后的剪切稀化特征进行探讨,研究了饱和砂土在地震荷载下的固-液相变特征。研究结果表明:大震激励时,饱和砂土层上部土体孔压比达到1,饱和砂土液化,从动剪应力-动剪应变曲线变化中可知液化后的砂土动剪切模量下降明显,这表明土体逐渐软化;从动剪应力、动剪应变推导得出的饱和砂土剪应力-剪应变率变化趋势与非牛顿流体的流变曲线相近,且饱和砂土上部液化土层的表观黏度大幅下降,液化后砂土表现出"剪切稀化"的假塑性流体特性。     

基于场地特征周期的瑞利阻尼确定方法

从地震影响系数曲线的特征周期 出发,提出了分别取 和 作为时域计算瑞利(Rayleigh)阻尼的两个目标频率。通过3种不同厚度的单一土层和1种多层土层模型在3条不同地震动作用下的计算,比较了基于Annie方法和本文方法确定Rayleigh阻尼的时域解,以及基于等效线性化的频域解,结果表明：两种Rayleigh阻尼形式的地表加速度峰值误差都在10%以内,地表加速度反应谱基本一致,但本文方法对土层厚度和输入的地震动不太敏感,即误差随着土层厚度及输入地震动的不同离散性很小,且应用更为简便,无需求土层的自振频率,只需要根据场地类别获得场地特征周期,就可确定Rayleigh阻尼值。     

固结比对黏性土动剪切模量影响的研究

采用共振柱试验方法,对固结比对黏性土最大动剪切模量和对动剪切模量比与剪应变非线性关系的影响进行了研究,并通过土层地震反应计算,初步给出了这种影响对地表反应谱作用的估计。通过对粉质黏土和粉土的系统试验,提出了计算不同固结比粉质黏土和粉土最大动剪切模量的建议公式,并与Hardin公式进行了对比。结果显示,固结比kc >1.0时的最大动剪切模量较kc=1.0时的增长程度明显比Hardin公式要大。与砂土试验结果相比,固结比对黏性土最大动剪切模量影响程度更大。同时,分析了固结比对动剪切模量比与剪应变非线性关系及对地表反应谱的影响。结果表明,固结比对黏性土动剪切模量比与剪应变非线性关系有一定影响,随固结比增加,动剪切模量比非线性关系曲线有一定提高,动剪切模量比值增大,在强地震动作用下的土层反应分析中应考虑实际固结比的影响。     

隧道开挖引起土层沉降槽曲线形态的分析与计算

假定由于隧道开挖在地表以下土层所形成的沉降槽的体积等于地层损失以及各土层沉降槽曲线仍可用正态分布函数表示。在此假定下,首次提出可以通过对隧道中心线以上土体的沉降值的分析,研究不同深度土层沉降槽曲线宽度系数,并且通过回归分析,提出用于计算不同深度土层沉降槽曲线宽度系数的公式,为计算隧道开挖引起的地表以下土体任一点的竖向沉降提供了一种有效方法,通过工程实例的计算证明了该方法的适用性。     

采矿爆破振动波在岩溶区的传播影响因素分析

为研究不同情况下爆破振动强度衰变规律和振动对岩溶塌陷的影响,为当地矿山合理开采及减少地面塌陷灾害的发生提供可靠依据,对湖南水口山铅锌矿区进行现场爆破振动测试,并利用古丹铅锌矿实测数据作对比分析;试验矿区共布设4条测线,接收8组爆破振动数据。采用萨道夫斯基修正公式对采集的数据进行计算,以爆破产生的振动波频率及振动速度作为测试指标,对实测数据进行提取、处理,确定爆破振动波的频率及其在介质中的传播速度及地震波引发的质点振动峰值振速。试验结果表明:采矿活动是岩溶地面塌陷的主要影响因素;爆破振动波的频率衰减强度与其在岩土体中的传播距离和断层有关,振动波的传播速度受到岩土体性质、岩层结构特征、岩层走向等因素的影响。      

缓倾地层微震定位算法探讨及其数值验证

倾斜层状地层中速度模型对于震源精确定位具有重要意义。基于水平层状地层两点间的快速射线追踪法,分析了应力波在倾斜地层中的传播路径及规律。结合网格搜索法,提出了针对倾斜地层的微震定位算法,并与传统的基于单一速度模型的定位算法进行对比分析;探索性地采用了颗粒流理论构建地层数值模型,通过颗粒间的相互作用模拟应力波传播进行定位算法有效性验证。研究结果表明：（1）对于倾斜地层,基于简化弹性波等速传播规律的定位算法精度较差,单一速度模型不能满足复杂介质的定位要求;（2）采用倾斜分层速度模型并改进应力波传播路径的计算方法能有效提高倾斜地层的震源定位精度,同时分析了地层倾角以及地层数对算法精度的影响,表明定位误差与倾角、地层数成正比;（3）通过减小网格搜索法的网格尺寸可在一定程度上提高定位精度,合理的尺寸大小和搜索策略有助于网格搜索法在实际工程中的应用。建立的倾斜地层震源定位及其验证方法可为进一步研究复杂地层定位监测技术提供重要的理论和技术支撑。     

A study of ground‐structure interaction in dynamic plate load testing

A mathematical treatment is presented for the forced vertical vibration of a padded annular footing on a layered viscoelastic half‐space. On assuming a depth‐independent stress distribution for the interfacial buffer, the set of triple integral equations stemming from the problem is reduced to a Fredholm integral equation of the second kind. The solution method, which is tailored to capture the stress concentrations beneath footing edges, is highlighted. To cater to small‐scale geophysical applications, the model is used to investigate the near‐field effects of ground‐loading system interaction in dynamic geotechnical and pavement testing. Numerical results indicate that the uniform‐pressure assumption for the contact load between the composite disc and the ground which is customary in dynamic plate load testing may lead to significant errors in the diagnosis of subsurface soil and pavement conditions. Beyond its direct application to non‐intrusive site characterization, the proposed solution can be used in the seismic analysis of a variety of structures involving annular foundation geometries. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical analysis of soil vibrations due to trains moving at critical speed

High-speed train induced vibrations of track structure and underlying soils differ from that induced by low-speed train. Determining the critical speed of train operation remains difficult due to the complex properties of the track, embankment and ground. A dynamic analysis model comprising track, embankment and layered ground was presented based on the two-and-half-dimensional (2.5D) finite elements combining with thin-layer elements to predict vibrations generated by train moving loads. The track structure is modeled as an Euler–Bernoulli beam resting on embankment. The train is treated as a series of moving axle loads; the embankment and ground are modeled by the 2.5D finite elements. The dynamic responses of the track structure and the ground under constant and vibrating moving loads at various speeds are presented. The results show that the critical speed of a train moving on an embankment is higher than the Rayleigh wave velocity of the underlying soil, attributed to the presence of the track structure and the embankment. It is found that the dynamic response of ground induced by moving constant loads is mostly dominated by train speed. While for the moving load with vibration frequency, the ground response is mostly affected by the vibration frequency instead of train speed. Mach effect appears when the train speed exceeds the critical speed of the track–embankment–ground system.     

Dynamic soil properties for microzonation of Delhi,India

Delhi, the capital of India, has experienced mild seismic shaking during several earthquakes in the past. The large variations of depth to bedrock and ground water table coupled with different soil types at different locations of Delhi necessitate a seismic microzonation study. Dynamic soil properties such as shear wave velocity, modulus reduction and damping characteristics of local soils are the basic and essential input parameters for conducting even a preliminary ground response analysis which is an essential input in microzonation studies. Shear wave velocity is not measured routinely due to its high cost and lack of the required expertise. Several researchers in the past developed correlations between shear wave velocity (V _s ) and routinely measured N values. In the present study, shear wave velocity profiles measured in the field at more than 80 borehole locations to a depth of about 20 to 32m using Spectral Analysis of Surface Waves (SASW) are presented and correlations between shear wave velocity and N values are also presented for use by engineers and designers. Results of strain and stress controlled cyclic triaxial tests on remoulded samples of sand-silt mixtures in the high strain range are used for generating the modulus reduction and damping curves and are compared with the well-known curves in the literature. The results presented in this article can be used for microzonation studies as well as site specific ground response analyses at Delhi.     

倾斜层位偏移规律研究

层状介质模型是对实际地质情况较为准确的近似。而在实际的数据处理过程中,以等效均匀介质代替层状介质,波传播的几何路径、介质速度的差别会使对倾斜界面位置的判断偏离实际情况。针对设计的理论模型和从测井数据提取的实际介质模型,笔者运用几何射线理论和波场数值模拟技术,定量计算并提取了层状介质模型和等效均匀介质模型两种情况下反射点在时间和空间上的偏移量。结果证明,等效均匀介质模型近似会造成对倾斜界面位置的判断不准确; 准确的速度模型是偏移归位研究的关键,应尽量使用由地震测井资料求得的瞬时速度。     

饱和软土大应变自重固结非线性分析

饱和软土的自重固结涉及许多复杂的非线性问题,而现有分析方法一般基于线性假设,将固结控制方程简化,在沉降变形较大的情况下并不合理。基于Gibson大变形固结方程,推导在单面排水和双面排水条件下的边界条件,综合考虑饱和软土的自重特性与非线性压缩、渗透关系,采用稳定性较好的Crank-Nicolson型差分格式离散并求解非线性控制方程,并从孔隙比、沉降、孔隙水流速3个方面将数值解与CS2固结数值模型进行对比验证。结合工程实例,研究土层不同初始厚度和初始孔隙比对自重固结的影响。结果表明：初始厚度与初始孔隙比较大土层最终沉降也较大。与传统方法相比,该方法所得出的数值解精度较高,实现了高度非线性固结方程的快速求解,其成果更能广泛运用于实际工程的分析之中。     

Estimation of freeboard requirements against overtopping of surface impoundments under seismic action

The excitation of structural components and liquid contents of surface impoundments by seismic waves can generate turbulence that is large enough to overtop the bounding berms. In cases in which the liquids are wastes from industrial/municipal operations, their release from impoundments can pose significant risks to the environment. In this analysis, the freeboard magnitudes that can accommodate liquid head levels in impoundments are determined through linkage of configuration of waves in the liquid surface to incident seismic wave characteristics, liquid characteristics and impoundment design. For an impoundment site in a region of ground acceleration levels ranging from 0.2 to 1.0 g and impacted by seismic shear wave velocity of 180 m/s, freeboard requirements are in the range of 0.004–2.0 m on soft soil; 0.008–0.7 m on medium-dense soil; and 0.002–0.1 m for dense soil. For the same impoundment design, ground acceleration and incident wave characteristics, freeboard requirements are directly proportional to the depth of the soil mantle over bedrock. The impoundment slope, which is a key parameter with regards to liquid holding volumetric capacity of the impoundment, is a less significant parameter than depth to bedrock with regard to the size of the required freeboard. This implies that siting of an impoundment should be considered to be critical to impoundment performance in seismic zones.     

运动荷载附近有限层厚软土地基的振动研究

基于Biot多孔弹性介质的波动理论,研究了运动荷载附近软土地基的振动问题。假设一条形均布荷载作用在地基表面,则该模型可视为平面应变问题进行分析。通过引入4个势函数和Helmholtz原理,并利用Fourier变换及逆变换技术,获得了运动荷载作用下软土地基的应力、位移和孔隙水压力的解答。利用离散Fourier逆变换得到数值计算结果,分析了荷载速度,频率以及软土的渗透系数及多孔弹性参数对运动荷载作用下地表竖向位移及土体中任一点孔隙水压力分布的影响。     

Evaluating the structural model from ambient vibration recordings and shallow seismic refraction techniques. A case study of Zagazig City,Egypt

This study uses the seismic refraction and noise measurements to investigate the velocity structure of the subsurface and emphasize the advantage of ambient vibration over the conventional seismic refraction technique. Field measurements were carried out at nine sites in and around Zagazig city. Shallow seismic refraction data were interpreted using the delay time method to obtain the two-dimension ground model at each site. Ambient vibration arrays are used to infer the one-dimensional compressional and shear wave velocity profiles through two main steps. The first step is to derive the dispersion curve from the recorded signals using the frequency-wavenumber method. The second is to invert the dispersion curve to obtain the site velocity profiles. The results of the compressional wave velocities obtained from seismic refraction technique showed that the subsurface consists of a number of layers ranging from two to four layers and give a good agreement with the results of the seismic wave velocities obtained from the ambient vibration arrays. The ambient vibration arrays gave a deeper depth of penetration than the other method, providing more information on the subsurface structure without any disturbance to the environment. This work provides reliable estimates of the seismic velocity structures of both shallow and deep sedimentary layers within the area of interest.