岩石中锚杆拔出试验数值模拟.pdf

水 　　利 　　学 　　报 　2002 年 12 月 SHUILI 　　XUEBAO 第 12 期 ( ) 文章编号 2002 012006806 岩石中锚杆拔出试验的数值模拟 1 1 1 杨 　强 , 任继承 , 张 　浩 ( 1 清华大学 水利系 , 北京　100084) 摘要 :采用二维格构模型分析了岩石中锚杆拔出试验。格构模型将材料在细观尺度上离散为由梁单元组成的格构系统。 单元采用简单的本构关系和破坏准则 ,假定单元的抗拉强度服从正态分布 , 以引入岩石的非均质性。在外载作用 下对整体网格进行线弹性分析 ,计算出格构中各单元的局部应力 ,超过断裂阈值的单元将破坏 ,材料的开裂过程通 过单元的依次破坏来模拟。研究中考察了材料非均质程度、约束条件对开裂机理的影响。数值结果表明 ,非均匀 性导致了开裂形态的非对称性 ,而增加水平约束后 ,试件的荷载 - 位移反应由脆性变为柔性。 关键词 :格构模型 ;锚杆 ;脆性材料 ;非均质 ;开裂 中图分类号:TU457 文献标识码 :A 岩体锚固技术在水利工程中被广泛运用 ,适用于地质条件较差的岩体边坡稳定加固、隧洞的支护和 混凝土结构的增强加固等。锚杆作为主要的加固手段 ,对岩体变形和强度起着重要的作用 ,但目前锚杆 和围岩相互作用的机理还不是很清楚 ,主要通过现场观测、实验测试或经验公式来分析锚杆对岩体力学 性质的影响。加锚岩体的数值模拟方法大都还是基于有限元法 , 以锚杆单元模拟锚杆,但一般都过低估 计锚固效果。作者使用格构模型来分析岩石中锚杆的拔出试验。格构模型[1～6 ] (Lattice Model) 是一种基 于细观结构考虑的数值方法 ,其特点在于引入简单的单元本构模型和简单的破坏准则 , 同时考虑介质的 非均匀性。虽然单元的力学行为比较简单 ,但宏观复杂的整体破坏现象可通过单元的相互作用反映出 来。 1 　模型和方法 ( ) 格构模型的基本原理是 : 1 连续介质可在细观尺度上离散成由弹性杆、梁等单元联结而成的格构 系统 ,每个单元代表系统的一小部分。网格一般为规则三角形或四边形 ,也可是随机的不规则网格。图 ( ) ( ) 1 为计算单轴拉伸的三角形格构网格。 2 单元采用简单的本构关系和破坏准则 如弹脆性本构关系 。 ( ) ( ) 3 引入统计分布规律以描述材料细观非均质性。单元之间参数 弹模或强度 可以互有差异 ,按某种规 ( ) 律分布 ,如正态分布 ,Weibull 分布等。 4 在外载作用下对整体系统进行线弹性分析 ,计算出格构中各 单元的局部应力 ,超过断裂阈值的单元将破坏 ,断裂后的单元从系统中除去以模拟材料的开裂 ,单元的 ( ) 断裂为不可逆过程。 5 单元破坏后 ,荷载将重新分配 ,再次计算以得出下个断裂单元。不断重复该计 算过程 ,直至整个系统完全破坏。单元的渐进破坏即可用于模拟材料的宏观开裂过程。 格构系统中的单元可采用杆、梁等不同单元。在二维情况下 ,每根梁能传递轴力 ,剪力和弯矩 ,每个 节点有三个自由度 ,即两个线位移和一个转角位移。从目前的研究来看 ,采用梁单元得到的结果较理 想 ,能够模拟出较为复杂的开裂形态。 收稿日期 ( ) 基金项目 : 国家自然科学基金