节理岩体变形参数数值模拟试验研究.pdfVIP

维普资讯 第7期 广东土木与建筑 No．7 20o6年7月 GUANGD0NG ARCHITECIURECIVILENGINEERING JUL20o6 节理岩体变形参数数值模拟试验研究 王 良奎 (广东省地质调查院 广州510080) 摘 要：通过获得岩石和结构面的物理力学参数，建立岩块和结构面的概化模型．运用三维有限元程序对不同岩块 和结构面组合成的岩体模型进行数值模拟试验 ，从而得到不连续面的不 同组合岩体变形参数的影响和变化规律 。 关键词 ：数值模拟试验 ；变形模量；岩体 ；岩块 ；节理 随着科技 的发展 ，计算数值法在岩土工程 中得 到了广泛应用，但 由于工程岩体的复杂性，使岩体计 算模型选择和岩体力学参数选取成为阻碍数值法发 展的两个瓶颈问题。其中，岩体力学参数的选择是否 合理，是各种数值方法在岩体稳定分析中能否成功 运用的关键 ，目前主要采用的研究方法包括传统的 (a)0。一3 (b)30o一9O。一5—5 图 1 三维计算模型 室内岩块试验法 、工程类比法 、经验法 、野外现场测 试法、声波测试法，以及位移反分析法、遗传算法和 神经网络法等 2 计算方法和参数的选取 针对室内岩石试块不能代表岩体属性 ，原位测 试周期长、费用大且数据少等缺点，人们在不断寻求 此次分析采用 日本某公司开发的三维有限元程 新的可行方法来解决这一问题 ．如遗传算法、神经网 序 3D—o-作为试验平台，选用Mohr—Culomb屈服准 络等。本文综合室内岩块和节理的资料，运用成熟 则 ，需要的计算参数主要有材料的变形模量、泊松 的有限元程序对岩体的变形参数进行数值模拟试 比、密度、内聚力、内摩擦角及抗拉强度等。 验 ，将是一种有益和探索性的研究。 以西南某坝址区的岩体作为研究对象，综合大 量的岩块室内试验获得岩石的物理力学参数 。计算 1 概化模型 模型中岩石的物理力学参数见表 l，节理作为软弱 层处理，其物理力学参数亦通过试验求取，见表 l。 在数值模型的建立过程 中，把岩块和节理看成 表 1 计算参数表 两种材料单元，前者采用空间六面体弹塑性单元，后 者假设为有一定厚度 的贯穿的不抗拉弹塑性单元 ， 节理厚度 2cm，间距设定为 20cm，同组节理相互平 行。按这样 的原则分别建立一组节理中包含 l一5条 及节理倾角 0~-90。变化的数值模型共 5O个 ，其中角 度变化间距为 lO。，两组节理中分别包含 l5条及节 理角度从 0~-90。变化的数值模型共 225个 。 3 计算边界的确定 试验数值模型为边长 2．0m的正方体，计算模型 编号分别为：①节理角度一节理条数 ，如 0。一3表示 采用无侧限单轴压缩试验 ，可获得极限强度下 模型中含有 l组 3条倾角 0。的节理，如图 Ia所示； 岩体破坏时的变形值 。为 了实现该试验条件 ，可从 ②第 l、2组节理角度一第 l、2组节理条数，如 3O。一 模型的边界条件上加 以设定来得到满足 ，这样在数 9Oo一5—5表示模型中3Oo角节理 5条和9O。角节理 5 值模拟过程 中模型计算边界条件被确定为 ：模型顶 条 ．如图 lb所示 部受压