滑坡对不同特性地震波的动力响应规律
发布时间：2019-08-08
摘要: 地震边坡是边坡工程和地震工程的重要对象之一，而边坡在地震作用下的动力响应规律则是研究的核心内容。本文以玉树机场路滑坡群 0#滑坡第七条块为研究对象，以大型有限元分析软件 FLAC3D为研究工具，通过对正弦波、玉树水平波、玉树竖向波不同加速峰值作用下滑坡动力响应规律的研究，得到如下结论: ( 1) 加速度放大系数在滑床内基本呈线性增大，但在滑面处有明显间断。( 2) 在正弦波、玉树水平波不同加速度峰值条件下滑带对加速度的影响规律一致。当加速度峰值较小时，滑带对加速度的影响较弱，随着加速度峰值的增大，滑带的滤波作用越来越显著。( 3) 相较于输入水平地震动，竖向地震动在坡体内引起的水平地震峰值加速度放大系数规律不太明显，仍需要进一步研究。
关键词: 机场滑坡; 动力响应; 放大系数; 加速度
0 引言
我国是大陆地震最集中、活动性最高的地区之一，而多山地、多地震的地理地质条件不可避免带来了大量和地震作用相关的边坡问题。同时地震边坡的研究是边坡工程和地震工程研究的重要内容，而边坡在地震作用下的动力响应规律则是研究的核心。滑坡的地震动力响应分析主要探究地震加速度峰值在坡体内的放大系数发展规律，以及震后滑坡残余变形特征，前者为滑坡动力稳定性分析提供基础，后者为滑坡动力失稳机理分析提供依据[1 - 3]。
对于边坡在地震作用下的动力响应规律，国内外学者做了大量研究。CELEBI 等[4]提出地震波具有明显的地形放大效应。刘洪兵等[5]进行国内外相关文献总结时，提出了相同的认识，认为地形放大效应与方向效应关系密切。张江伟等[6]建立了二维均质边坡有限元模型，模拟计算了其在地震作用下的动力响应得出在坡顶和坡面处，与输入地震动加速度时程比较，输出加速度峰值出现的时刻有滞后现象。邹烨等[7]从三维均质边坡出发，给出了一种将矢量和法与动力有限元法相结合的三维边坡动力稳定性分析方法，并给出其具体计算过程，最后将该方法运用到大型堆积体边坡的动力稳定性分析中。 HUTCHINSON[8]基于现场滑坡模型试验以及室内试验，提出了导致土质液化形成滑坡的认识。陈臻林等[9]通过数值模拟研究了地震波作用下含反倾软弱夹层岩质边坡动力响应规律。
上述学者在滑坡动力失稳机理研究方面较多，对滑坡动力稳定性研究较少，尤其是地震波特性对滑坡动力响应规律的影响研究更少。因此，本文以玉树机场路滑坡群 0#滑坡第七条块为原型，运用有限元分析软件 FLAC3D建立数值分析模型，研究正弦波、玉树水平波、玉树竖向波的加速度放大系数在坡体内的发展规律，为地震作用下滑坡稳定性分析提供借鉴。
1 工程概况
青海玉树机场路滑坡群 0#滑坡坡面坡度 20°，距离线路中心约 212 m，高程在 3 776 m 附近为滑坡体后缘，其地形较缓。滑坡右侧发育有流向为 NE14° 的“U”型沟，沟宽 26 m。第七块滑坡滑动方向 NE46°，主滑段倾角约 18° ～ 19°，沿线路宽 133 m，垂直线路长度 166 m，滑体体积 3. 82 × 105 m3 ，滑坡分为三级。堆积层主要由角砾土、碎石土以及黏土构成。由于 0#第七条堆积层滑坡坡体结构及构成相近，简化为二维模型 ( 图 1) 。
0#滑坡第七条块的地球物理特性如图 2 所示，其滑坡介质电阻率分布较均匀且物性层位清楚，能够反映边坡内部情况及滑坡发育情况。两剖面有非常好的对应关系，表层局部阻值较高，而中层较低且形成连续的闭合低阻闭合圈，为低阻软弱含水层或过湿带，后部高阻区为基岩。
2 计算模型建立
以青海玉树机场路滑坡群 0#滑坡断面为对象，利用 FLAC3D软件建立模型，工程地质断面如图 3 所示。该滑坡为多滑面堆积层典型坡体结构，根据计算需要，数值模型共离散为 2 407 个单元和 5 106 个节点( 图 4) 。计算时，模型顶部自由，左、右边界施加吸收边界条件，以减小地震波的反射，底部施加水平地震加速度时程，以模拟地震荷载影响。
3 模型设计参数
3. 1 材料参数
滑带抗剪强度参数黏聚力取 30 kPa，内摩擦角取 20°，动力作用结束后的残余强度参数中黏聚力为 15 kPa，内摩擦角 18°。其余各计算参数如表 1 所示。
3. 2 地震荷载参数
为了说明不同地震波特性对边坡的地震反应规律，选用正弦波、玉树水平波、玉树竖向波来进行模拟分析。其中正弦波固定频率为 2 Hz，峰值取不同值，分别为 0. 05g、0. 1g、0. 2g、0. 4g。
玉树地震中记录到大量的强震动加速度时程曲线。其中，飞机场布设有强震观测台站，距离结古镇约 19 km。机场路滑坡的临空面指向 E，故本文选用 EW 向地震波作为玉树水平波的动力输入，EW 向地震波如图 5( a) 所示。玉树飞机场强震记录历时较长，而对滑坡稳定性影响最为明显的是 30 ～ 40 s 间的地震荷载，从中截取加速度时程曲线如图 5 所示。
4 计算结果及分析
为便于比较峰值加速度放大系数在竖向上的变化，在计算域内指定了三条参考线( 图 6) ，分别标示滑体前部( F) 、滑体中后部( M) 和滑体后部( R) 。为了消除模型边界( 尤其是底部边界) 对计算结果的影响，计算结果( 峰值加速度放大系数) 分析中，计算域各量均是相对于参考点( RP) 的成果，其中，放大系数是计算域各点峰值与参考点峰值的比值。
4. 1 输入正弦波的影响
输入的正弦波峰值加速度分别为 0. 05g，0. 1g， 0. 2g，0. 4g，通过比较计算结果( 图 7) ，可以得出加速度放大系数随高度增加而增大，最大值出现在坡顶后缘。峰值加速度为 0. 05g 时，滑块加速度放大系数分布与滑床相同。峰值加速度为 0. 1g、0. 2g 和 0. 4g 时，滑床上加速度放大系数发生无规律的变化，变化速率随正弦波峰值加速度增加而增大。
从图 7 中数据所示的参考线峰值加速度放大系数来看，滑坡地震响应的特点: ( 1) 加速度放大系数在滑床内基本呈线性增大，但在滑面处有明显间断。显示出滑带对加速度具有强烈的影响，其强度较低，滑带出现塑性屈服后，加速度峰值不能有效的传播到滑体内。 ( 2) 当加速度峰值较小时，滑带的影响不明显。放大系数随高度基本呈线性增大。随加速度峰值逐渐增大，滑带的滤波作用越来越明显。当峰值加速度达到 0. 4g 时，甚至出现了滑体内放大系数小于滑带的现象。( 3) 从空间上来看，滑带的滤波作用不是同时发挥出来的。就本滑坡而言，滑体前部首先出现塑性屈服，滤波作用显现。如图 7( a) 可见。
4. 2 输入玉树水平波的影响
玉树水平波峰值加速度输入值分别为 0. 05g， 0. 1g，0. 2g，0. 3g，0. 4g。通过图 8 计算结果的比较，得出加速度放大系数最小值为 0. 6 左右，出现在坡顶下方距底面 1 /3 滑床高处，向四周辐射增加。对于滑块，当峰值加速度为 0. 05g，滑块并未完全破坏，放大系数随高度增加而增大。随着峰值加速度的增加，滑块放大系数变化剧烈，分布无规律。而整个坡体在此过程中，放大系数在 0. 8 ～ 1. 2 的区域不断扩大，峰值加速度达到 0. 4g 时，已有 4 /5 的土体的放大系数处于这一区域。
从图 8 中数据来看，玉树水平地震波对加速度放大系数的影响规律与正弦波基本类似。加速度放大系数在滑床内呈线性增大，但在滑面处间断，说明滑带对加速度具有强烈的影响。当加速度峰值较小时，滑带对加速度的影响较弱，随着加速度峰值的增大，滑带的滤波作用越来越显著。当峰值加速度达到 0. 3g 及以上时，甚至出现了滑体内放大系数小于滑带的现象。
4. 3 输入玉树竖向波的影响
输入的地震动竖向波峰值加速度分别为 0. 1g， 0. 2g，0. 3g，0. 4g。在边坡滑块没有发生较大塑性变形以致滑移破坏时，即加速度峰值为 0. 1g 工况，在垂直方向上，加速度放大系数表现出随着高度的增加而增大的变化趋势; 在水平方向上，当峰值加速度为0. 1g 时，加速度等值线向内倾斜，即相同高程距离坡面越近，放大 系 数 越 小; 当峰值加速度为 0. 2g，0. 3g 和 0. 4g 时，滑块发生较大塑性变形，有整体滑移的倾向，加速度放大系数由 0. 6 逐渐减小到 0. 4( 图 9) 。
从图 9 中数据所示的参考线峰值加速度放大系数来看，相较于输入水平地震动，竖向地震动在坡体内引起的水平地震峰值加速度放大系数规律不太明显: ( 1) 当加速度峰值较小时( 0. 1g) ，计算域内引起的水平加速 度 较 为 稳 定，随高度放大系数有一定增大。 ( 2) 当竖向加速度峰值较大时，水平加速度峰值变化规律性不强，放大系数较小( 一般小于 0. 5) 。( 3) 竖向加速度所引起的水平加速度变化，仍需要进一步研究。
5 结 论
通过对玉树机场路滑坡群 0 #滑坡第七条块坡体在正弦波、玉树水平波、玉树竖向波的不同加速度峰值条件下边坡的地震反应特性研究，得出如下结论:
( 1) 加速度放大系数在滑床内基本呈线性增大，但在滑面处有明显间断。
( 2) 在正弦波、玉树水平波不同加速度峰值条件下滑带对加速度的影响规律一致。当加速度峰值较小时，滑带对加速度的影响较弱，随着加速度峰值的增大，滑带的滤波作用越来越显著。
( 3) 相较于输入水平地震动，竖向地震动在坡体内引起的水平地震峰值加速度放大系数规律不太明显，仍需要进一步研究。
滑坡对不同特性地震波的动力响应规律相关论文推荐：震陷滑坡建筑震害与设计
