Udec计算算例分析和总结.docxVIP

UDEC 算例 工程概况 某隧道位于一包含高角度连续节理岩体内，节理倾角为 50 度，平均间距为7m，隧道为一半径为 9m 的圆形隧道。贯穿于开挖面内的一垂直断层，在隧道拱顶形成了一个三角楔形体。 本算例使用 UDEC 的结构单元逻辑来模拟喷射混凝土和锚杆联合支护的圆形隧道开挖问题。 构建模型 隧道埋深 451m，为半径 9m 圆形隧道，本次计算模型左右边界取 41m，隧道至上下边界也取 41m。总的来说，模型长 100m，宽 100m。计算模型如图 1-1 所示。 计算参数 图 1-1 UDEC 计算模型 在包含高倾角节理和垂直断层的岩体内进行圆形开挖的 UDEC 模型岩体、节理和断层参数如下所示： 密度 Dens（kg/m3 密度 Dens （kg/m3） 体积模量K （Gpa） 剪切模量G （Gpa） 2500 1.5 0.6 法向刚度KN 表 1-2 切向刚度 KS 节理、断层物理力学参数 抗拉强度 jten 粘聚力 c 摩擦角φ （Gpa/m） （Gpa/m） (Pa) （Pa） （o） 2 2 100 100 10 密度 dens（Kg/m3） 密度 dens （Kg/m3） 2500 弹性模量 E/Gpa 21 泊松比ν 0.15 拉伸屈服强度 （Mpa） 2 残余屈服强度 （Mpa） 1 抗压屈服强度 （Mpa） 4 表 1-4 岩体和喷射混凝土接触面物理力学参数 法向刚度KN 切向刚度 KS 抗拉强度 jten （Gpa/m） （Gpa/m） (Mpa) 1 1 1 粘聚力 c （Mpa） 1 摩擦角φ （o） 45 锚杆模量 表 1-5 极限抗拉强度 锚杆物理力学参数 浆液粘结刚度 浆液黏聚强度 （Gpa） (MN) (GN/m/m) (MN/m) 100 10 1 1 作为演示的目的，隧道开挖和支护是瞬时发生的。本算例共两种支护分析被计算：第一，只施加喷射混凝土衬砌；第二，喷射混凝土和锚杆联合提供支护。 为了在第二种支护情况分析中可以更清晰的看到锚杆提供的支护，算例采用喷射混凝土的抗压强度被设置成一个很低的值，且厚度仅取为 10cm。 模拟步骤 建立模型 在 UDEC 中输入以下命令可建立隧道结构模型及边界。如图 1-2 所示。 new round 0.1 block -50,-50 -50,50 50,50 50,-50 jset -50,0 100,0 0,0 7,0 ；刷新 UDEC 窗口，重新调用一个新程序 ；块与块之间的圆角半径，必须小于块体最小边的1/2 ；建立模型框架 ；设置节理 crack -6 -50 -6 50 tunnel (0,0) 9,16 del range area 0.08 gen edge 10 ；设置断层 ；模拟开挖隧道边界 ；删除面积小于 0.08 的块体 ；自动划分单元，单元最大边长不超过 10 设置单元属性和材料特征 图 1-2 初步模型图 在 UDEC 中输入如下命令设置单元属性和材料特征。 prop mat 1 de 2500 k 1.5e9 g 0.6e9 ；定义围岩密度、体积模量、剪切模量（软件默认的是 线弹性模型，可以通过 change 命令改变本构关系） prop jmat=1 jkn 2e9 jks 1e9 jfric 10 jcoh 100 jtens 100 施加约束及荷载 施加约束条件 ；定义节理与断层的法向刚度、切向刚度、摩擦角、黏 聚力以及抗拉强度（默认节理面遵循摩尔库伦准则计算， 可以通过 change 命令设置） 输入如下命令对模型施加约束条件，对围岩左侧和右侧节点，设置约束条件为 UX 被约束，围岩底部节点，设置约束为 UY 被约束。如图 1-3、4 所示。bound stress 0 0 -10e6 range yran 49 51 ；上侧边界施加 400m 埋深等效荷载 bound yvel=0 range yran -51 -49 bound xvel=0 range xran -51 -49 ；约束下边界 ；约束左侧边界 bound xvel=0 range xran 49 51 ；约束右侧边界 图 1-3 左右侧边界施加约束 施加自重 图 1-4 下侧边界施加约束 设置 Y 方向的加速度 10m/s2，生成 Y 轴反方向的惯性力。 set grav 0 -10 ；设置重力加速度 计算结果分析 求解初始应力场 岩体在自重应力下会产生初始应力场，这时还未开挖，未施加支护结构，求解得到初始应力场。相应的命令如下： solve rat 1e-5 save supp1.sav ；求解 ；保存（默认为安装目录，可以设置保存路径） 得到的初始应力场（图 1-5）和初始位移场（图 1-6）如下所示，输入命令 如下：