岩石高边坡爆破振动局部放大效应的有限元分析.docxVIP

岩石高边坡爆破振动局部放大效应的有限元分析 中国有很多大型和大型水电工程，如金沙江的西洛渡、白鹤滩、乌东、江家水库和雅鸡江的锦川一级和二级。这些项目分布在中国西南部的崇山溪和高山峡谷。这些电工工程中存在着复杂的设计斜坡问题。这类高边坡具有如下特征:天然边坡陡峭, 坡角一般超过45°, 大多在70°到90°之间;工程边坡高, 可达300~600 m以上。边坡岩体爆破开挖过程中, 爆破振动的高程放大效应是边坡上振动速度传播规律的重要研究内容之一。 大量现场爆破振动监测研究成果均表明高程放大效应是高边坡爆破振动动力响应的典型特征, 它受到边坡的高差、坡比、岩体地质条件及爆破规模等因素的影响。许多学者采用理论分析、试验分析、数值模拟对边坡的爆破振动高程放大效应进行了研究。石崇等利用平面射线理论和方法, 推导出单面边坡地震的高程放大效应理论解。钟冬望等通过对实验室边坡模型试验研究表明, 临近边坡爆破时, 边坡具有振动速度放大作用, 且振动速度的衰减速度变慢。杨风威等结合振动监测和数值分析方法研究了爆破开挖时边坡的振动速度衰减规律, 发现边坡爆破振动速度随高程增加存在放大效应, 并且这种放大效应在水平径向比垂直向更为显著。唐海等的研究表明, 凸形地形、地貌对爆破振动波垂直向的放大效应大于水平向, 放大效应具有方向性, 并且通过量纲分析法得到一个能反映正高程差放大效应的爆破振动公式。陈明等通过理论及数值分析得到在边坡坡形骤变、坡度增大的坡段会产生显著的振动速度高程放大效应。国外方面, Marrara、Havenith及Vladimir等也对边坡的振动放大效应进行了研究。 在爆破振动测试方面, 测点的布置占有极其重要的地位, 它直接影响爆破振动测试的效果以及观测数据的应用价值, 测点数目过少, 观测数据不足以说明问题;测点数目过多, 则所需一起数量及测试工作量较大。如果测点布置不当, 对不合理的测点进行分析则不能反映工程实际情况, 甚至导致错误的结果。 基于LS-DYNA的动力有限元计算, 结合溪洛渡高边坡开挖背景建立数值模型, 研究岩石高边坡振速高程放大效应, 进而考虑实际工程中的岩体开挖损伤与夹层等条件下爆破振动的分布规律与局部放大特征, 分析了不同典型位置作为实际工程爆破监测点的利弊, 为爆破振动测点的布置与优化提供依据和参考。 1 透射边界模拟 以溪洛渡高边坡的开挖为工程背景, 建立如图1所示的边坡数值计算模型, 根据工程资料显示, 一次开挖的台阶高度为15 m, 马道宽度为3 m, 模拟其中某一台阶光面爆破的振动响应。模型整体示意图以及炮孔附近网格放大图如图1所示, 光爆孔孔径90 mm, 药径32 mm, 堵塞2 m, 除自然坡面外其余边界均施加透射边界以模拟无限岩体。计算中本构模型采用弹性模型, 基于相关文献, 岩体密度为2600 kg/m3, 弹性模量为20 GPa, 泊松比为0.228。 计算中采用LS-DYNA的流-固耦合 (ALE) 算法模拟炸药的动力冲击作用, 其主要原理是实现材料在网格中流动, 从而可以处理流体流动、爆炸冲击等大变形问题。炸药用mat high explosive burn材料实现, 结合JWL状态方程来模拟炸药爆炸过程中压力与体积的关系, 表达式如下 空气采用mat null材料实现, 同时结合式 (2) , 采用多线性方程描述空气的作用 式中:C0=C1=C2=C3=C6=0;C4=C5=0.4;μ=ρ/ρ0, 其中, ρ, ρ0分别为初始与当前的材料密度。 2 爆心高差对爆破振动的放大作用影响 众所周知, 当建筑物受地震作用时, 若其顶部存在小突出部分, 则小突出部分由于质量和刚度比较小, 在每一个来回转折的瞬间, 形成较大的振动速度, 该现象在工程抗震中被称为鞭梢效应。对于岩石高边坡, 其本身是一个复杂的岩石结构体, 但如果将边坡不同高程台阶坡脚一一连接起来, 其边坡台阶部位的岩体相当于边坡主体结构的突出物, 而实际工程中这些部位是爆破振动监测点的主要布置区域, 因此有必要研究台阶突出物的局部爆破振动特征。 统计如图2所示的坡面控制点的水平向与竖直向的爆破振动峰值, 得到各控制点爆破峰值振速与爆心高差的关系如图3所示。 由图可知, 随着爆心高差的增加, 单孔爆破振动的峰值速度大体上呈现出衰减趋势, 马道外缘点和马道中点出现明显的振动放大效应, 而马道内侧点和坡面中点则表明在该处爆破振动放大并不明显, 由此可见, 岩石高边坡的爆破振动放大效应针对单层马道而言只是在局部部位出现放大, 主要体现在马道外缘附近;同时从图3中可以看出, 随着爆心距的增加, 局部的这种爆破振动放大效应逐渐减小。 为进一步对这种局部放大效应的分布特征进行分析, 沿坡面和马道大量取点, 图4给出了这些密集取点的水平径向爆破