结构面对岩体边坡地震动影响的数值模拟研究Word文件下载.docx

1、顺倾结构面会使坡肩部形成更强烈的地震动响应；坡角越陡，坡顶和坡肩部的地震动响应越强烈。上述规律对于岩质边坡地震稳定性分析具有指导意义。关键词：岩体边坡，结构面，地震动加速度放大系数，数值模拟Numerical Analysis of the Influence of Structure Planes on the Earthquake Motion in Rock SlopesSUN Jin-zhong TIAN Xiao-fuChina University of Geosciences （Beijing）, Beijing 100083, ChinaAbstract: The numeri

2、cal simulation with software UDEC for the influence of structure planes on the earthquake motion of rock slopes are carried out in this paper. The study deals with the influence of the variations in the physical and mechanical parameters and attitudes as well as positions of structure planes on the

3、intensification of earthquake motion in rock slopes. According to general situations of structure planes developing in rock slopes, two attitude types of structure planes are generalized in this paper: 1） penetrating structure planes intersecting with the slope surface; 2） structure planes parallel

4、to the slope surface. The first type of structure planes is relative to strata layers and long penetrating weak structure planes, and the second type mainly stands for unloaded joints of the rock slopes. It is revealed by numerical analyses that the variation in the damp ratio of rock mass can affec

5、t the intensification of earthquake motion in the slope to some extent, but it has little influence on the distribution of the earthquake motion in the slope; a structure plane will have significant influence on the earthquake response of the slope only if the stiffness of the structure plane is at

6、a proper value; the shorter the distance between a structure plane and the top surface of the slope, the stronger the earthquake response of the slope; the structure plane across and dip to the slope surface will cause stranger earthquake motion in the upper slope body than that caused by structure

7、planes horizontal or anti-dip to the slope surface; the stepper the angle of the slope surface, the stronger the earthquake motion in the top or shoulder of the slope. The laws above-mentioned will be valuable references to the of earthquake stability analysis of rock slopes.Keywords: rock slope, st

8、ructure plane, amplification coefficient of seismic acceleration, numerical analysis.1 引言我国是一个地震多发的国家，而地震是边坡失稳的主要诱因之一1，例如，汶川5.12大地震触发了大量的边坡崩滑，带来了不可估量的损失。所以，边坡地震稳定性问题已成为岩土工程研究的热点问题28。目前，对于边坡的地震稳定性评价主要采用拟静力法和有限滑动位移法，这两种方法都需要确定地震系数。地震系数 定义为一点的地震动加速度 a 与重力加速度 g的比值（=a/g），可见，边坡地震系数的确定需要明确边坡对于地震作用的响应规律。祁生文

9、、何蕴龙、徐光兴等针对边坡的地震响应问题开展了数值模拟研究，得到了一些有价值的研究成果911。这些成果主要针对的是均质边坡情况，对于存在着大量结构面的岩体边坡尚有待进一步研究，所以，本文展开了结构面对岩体边坡地震动影响的二维数值模拟研究。鉴于离散元软件UDEC可以很好地解决离散岩体的非线性动力分析问题，已经广泛应用于采矿工程、地下结构工程和边坡工程等领域12，本文的数值模拟采用了UDEC软件。为便于边坡内不同位置地震动强度的比较，研究中定义边坡岩体中一点地震动峰值加速度 aP与坡脚处地震动峰值加速度aP0之比为边坡的地震动加速度放大系数（=aP/aP0）。本文研究内容涉及结构面物理力学性状变化

10、、产状和位置变化对岩体边坡地震动加速度放大系数空间展布的影响。2 含结构面边坡动力响应模拟方案2.1 基本假定本文在结构面对于岩质边坡动力响应影响的研究中遵循如下假定：1）应力应变关系：平面应变状态；2）岩石材料：均匀的各向同性材料；3）结构面本构关系：莫尔-库伦滑动模型；4）入射地震波：从模型底面垂直入射的剪切波。2.2 模型工况选择岩体中结构面的分布具有不确定性，而且其发育的规模和特征也有很大的不同，考虑现实中所有情况下结构面对于岩质边坡动力响应影响存在着很大的困难。为了突出主要矛盾，考虑岩体边坡中结构面发育的一般情况，本文概括出两种类型的结构面：前者可以用以模拟边坡岩体中的层面、贯穿性长

11、大软弱结构面等；后者主要模拟边坡岩体卸荷节理。本文主要研究以下两种工况：工况1在边坡模型中布置一条与坡面相交的贯穿性结构面，考虑坡角、结构面与坡面交点的位置、结构面刚度以及结构面倾角这些参数的变化对边坡动力响应的影响。工况2在垂直于坡面向内一定范围内布置几条平行于坡面的结构面，上部切穿坡顶面，下部截止于一条水平贯穿性结构面。结构面刚度值设为较低值，结构面抗拉强度设为零。此种工况主要模拟卸荷带对于岩质边坡动力响应的影响。模拟中主要考虑卸荷裂隙发育密度和卸荷带深度的变化。两种工况的模型示意图如图1、图2所示。图中：D为边坡水平长度，H为边坡高度，h为结构面与坡面交点位置高度， 为坡角， 为结构面倾

12、度，d为卸荷裂隙间距，l为卸荷裂隙带的高度，s为卸荷带深度， f为入射地震波频率。模拟中主要关注结构面附近加速度放大系数等值线的分布与坡肩放大系数的量值。图1 工况1边坡动力响应分析模型示意图图2 工况2边坡动力响应分析模型示意图2.3 模型边界条件与参数选取为了消除截断边界对于应力波的反射效应，模型底面边界采用粘滞边界，模型两侧采用自由场边界13。为便于边坡内不同位置地震动强度的比较，引入本文开始定义的无量纲参数，边坡地震动加速度放大系数。为了全面了解边坡各个部位的动力响应，利用程序自带的FISH语言编制了自动布置监测点的函数，记录计算过程中监测点上地震动加速度的变化过程，并自动计算每个监测

13、点上地震动加速度的最大值。计算持续时间应该保证入射地震波传播到坡顶并发生反射叠加作用，待波场稳定后再进行加速度最大值的提取。各种工况模型参数按表1选取。表1 模型参数表材料参数符号取值块体密度2500kg/m3体积模量K27.7GPa剪切模量G20.8GPa泊松比0.2纵波速度Vp4714m/s横波速度Vs2886m/s结构面粘聚力c2MPa结构面内摩擦角402.4 地震荷载地面地震动过程一般以水平方向振动为主，频率成分复杂，地震动加速度主频一般在210Hz。为了进行一般规律性研究，本次数值模拟地震动力荷载输入采用水平向简谐振动剪切波，从模型底部边界垂直入射。首先确定输入地震动加速度时程，然后

14、对其进行积分转化成速度时程，再将其转化成应力时程从边界输入。选择加速度振幅为1m/s2、频率为5Hz的简谐水平地震动作为输入地震荷载，进行边坡地震动响应的数值模拟研究。同一测点上地震动的位移、速度或加速度均随输入地震荷载的强度的增加而增加，三者的强弱变化一致，在边坡中的强弱分布形式相同9,14。所以，在边坡地震响应数值模拟中可选用地震动加速度放大系数来表示边坡地震动响应强度的分布。3 边坡介质物性参数的影响3.1 阻尼比应力波在阻尼比不同的材料中传播时，随着传播距离的增大，其幅值和频率特性的差异也会随之增大15。因此，在进行含结构面边坡动力响应分析前，需要了解阻尼比对模拟结果的影响。建立坡高60m、坡角30的均质边坡模型，阻尼比 取值在1%5%的范围内，模拟结果如图3所示。从图3可以看出，阻尼比在1%5%范围内变化时对边坡的加速度放大系数的分布基本上没有影响，只是边坡加速度放大系数值随着阻尼比的降低而升高，但量值上变化并不显著。这说明在模拟阻尼范围内，边坡的动力响应加速度分布对于阻尼比的变化并不敏感，所以在以下的数值模拟过程中边坡介质的阻尼比 固定为3%。（a） =5%（b） =3%