UDEC模拟
大众矿围岩运移的UDEC数值模拟
ｇ 为 重 力加 速 度（ 力 ） 体 。 将 上 式 进 行 积 分 ， 得块 体 新 的状 态 ： 可
（＋ ｔ２ ｔ Ａ／ ）
（） ２ 节理 和 角 点 圆弧 化 。 岩石 的 节 理 可 以 认 为 是 块 体 之 间 的 接
触面 。 由于 粗 糙 度 的 关 系 ， 触 面 系 由个 别 接
值 模 拟 的 发 展 ， 用 这 种 方 法 来 取 代 一 部 采
分的相似模拟 ， 某些方面 ， 在 能得 到相 似 模
拟 所 达 不 到 的 效 果 。 了 更 全 面 地 从 各 个 为 方 面 研 究 大 众 煤 业 １ ０８ 开 采 工 作 面 的 岩 ｌ ２
∑ 为块体上的转动惯量和；
当考 虑体 力时 ， 于 二维 块 体 ， 据 牛 移 约 束 ， 边 界 采 用 固 定 ｘ方 向位 移 约 束 ， 对 根 右
并 方 面 的规 律 ， 于这 些 观 测 研 究 ， 能 对 顶 顿 第 二 定 律 ， 由中心 差 分 式 得 速 度 方 程 ： 基 不
“ ： 啦ｆ一 △ ｆ —＋１ ＋ ｍ ｇｆ ｌ Ｉ
层移 动 及 应 力分 布 规律 ， 用ＵＤ Ｃ 采 场 采 Ｅ 对
围岩进行模 拟。
的 接 触 点 所 构 成 ， 般 可 用 两 个 端 点 的 接 一 触 来 代 表 ， 运 动 时 端 点 有 相 对 的 法 向 位 在
移 和切向位移 。
（ ＝ ’ 『 ） ＋ｕ ＋
‘ ） ＝ （ ＋０ ｆ ’
＝ｚ ｆ 十二 一
ｍ
６ 。 文 计 算 所 选 用 计 算 机 的 内 存 为 ｍ 本
１２ Ｍ Ｂ 。 ８
Ａｔ
地震作用对采空区塌陷的UDEC模拟
收稿日期:2007-05-20作者简介:伍永田(1969-),男(汉族),广西桂林人,中南大学地学与环境工程学院在读博士研究生,主要从事企业管理与有色金属相关产业研究。
地震作用对采空区塌陷的UDEC 模拟伍永田1,2 张旭生2 李晓芸3(11中南大学,湖南长沙410083;21广西堂汉锌铟有限公司,广西南丹547200;31长沙有色研究设计院,湖南长沙410001) 摘要:采空区的稳定性问题除了岩体强度、节理、裂隙、层面等地质缺陷自身的因素外,地震、爆破、地下水等外界的因素也不容忽视。
本文利用离散元软件UDEC 模拟了地震波对天德石膏矿采空区塌陷的影响,分析了节理围岩在地震波作用下的运动情况,指出地震波的作用是采空区产生塌陷的诱发因素。
关键词:采空区;地震;离散元中图分类号:TD 32513 文献标识码:A 文章编号:1671-8550(2007)06-0019-030 引言当采矿工程中的采空区、大型水利水电工程的地下硐室等处于中等及以上强度的地震活动带时,必须考虑地震作用的影响。
与爆炸荷载相似,地震荷载也是以波的形式在岩土介质中传递。
在传递过程中,地震波产生的震动效应造成岩体的损伤,影响岩体的完整性,进而影响地下工程的安全[1,2]。
因此,分析地震荷载作用下地震波在采空区围岩岩体中的传递规律,进一步评价采空区的稳定是迫切需要解决的问题。
由于岩体中含有层面、断层、节理等不连续面,在整体上表现为一种非均质、不连续体,一般情况下很难准确地确定地震波在岩体介质中的传播规律,常用的方法就是利用地震监测仪根据过去的实测数据进行分析处理,进而确定地震波的衰减规律。
但岩体的性质具有很大的离散性,不同地方具有不同的岩体属性,从而地震波的衰减规律也不同[3～6]。
由美国Itasca Consulting Group ,Inc 1推出的UD EC (Universal Distinct Element Code )软件[7]是一种基于非连续体模拟离散元的二维数值计算程序,不仅能方便地模拟节理岩体的运动(包括移动和转动),还可模拟地震波等在岩体中的传播规律,本文用UDEC 程序对石膏矿开挖后围岩在地震波作用下的情况进行模拟分析,旨在验证采空区在地震作用下的稳定性。
用UDEC模拟爆炸波在节理岩体中的传播
万方数据 需 指 出 的 是% 用 3)4 $ 5* 678 *计算爆炸波
第 (期
赵 坚 等 ~ 用 @C t ,模 拟 爆 炸 波 在 节 理 岩 体 中 的 传 播
> 界和右边界 #$% 均设为不反射边界 & 而左边界和底边界则被设为轴 对称边界 ’
" " R 量的消耗 Q ’ Q J R 0 + X . & & 提 出的 离散元 数值计 算法 是 专 门 用 来 解决 节 理 岩 体 不 连 续 性 问 题 的 方 法 ’ 然 而I 现有
数量和 位 移 量 均 很 小 的 情 况 ’ 离散元法D G H 9E
Q J R
是专门用来解决不连续介质问题的数值模拟方法 I 该法 把 节理岩体 看 成 是 由 离 散 的 岩 块 和 岩 块 间 的 节理面所组成 ’ 岩块能移动 P 转动和变形 I 而节理面 可被压缩 P 分离或滑动 ’ 因此 I 岩体被看作一种不连 续的离散介质 ’ 离散介质内可存在大位移 P 旋转 P 滑 动乃至块体的分离 I 从而可以较真实地模拟岩体中 的不连续面 ’ 在地下岩石硐室中爆炸是一个快速而复杂的 过程’ 当大 量 高 性 能 爆 炸 材 料 在 岩 硐 中 被 引 爆 时 I 爆炸作用产生的高温 P 高压使岩硐内周围气体体积 迅 速膨胀 I 形 成 一 种 爆 炸 冲 击 波I 冲击波随后在岩 体中传播
图 # 3)4 $ 5* 678 * 与 )* + ,结合进行 计算的程序框图 ( # , = > ? @ AB C D E D > @ F E G B H @ I J K C H JC L > F ? $ 3)4 5* 678 *E F K)* + ,A@ K J G L
udec模拟实例
6.4喷射混凝土UDEC模拟6.4.1 UDEC简介刚体离散单元法一般认为Cundall于1971年提出来的。
该法适用于研究在准静力或动力条件下的节理系统或块体集合的力学问题，最初用来分析岩石边坡的运动。
该法是在牛顿第二定律的基础上建立起来的, 假设块体为准刚度体，块体运动主要受节理或弱面控制。
刚性块体的假设对于应力水平较低的问题，如边坡稳定是合理的。
将节理岩体视为由裂隙切割的非连续介质，相互切割的裂隙将岩体分成相互独立的块体单元，单元之间可以看成是角-角接触、角-边接触或边-边接触。
块体间的边-边接触可分解为由两个角-边接触而成，并且随着单元的平移和转动，允许调整各个单元之间的接触关系，最终块体单元可能达到平衡状态，也可能一直运动下去。
这些块体在平衡条件发生变化时，块体之间就产生相互作用力，从而导致块体产生一定的加速度和位移，使块体的空间位置和状态发生变化。
运动的块体之间，由于差异位移矢量的存在，从而使块体之间又发生新的作用力，根据新的力系，又可以计算出来各个块体在新的力系下的加速度、位移及新的运动位置。
如此反复迭代直到整个体系在新的力系作用下达到平衡状态为止，这样整个岩体的破坏运动过程就被真实的模拟出来。
离散单元法可以对由不同块体构成的整体进行应力、应变的分析计算，各不同块体之间通过接触点的耦合而互相连接在一起。
就大多数岩体来说，其构造弱面的刚度和强度均比岩石本身要小得多，从这点出发，为了减少研究对象的不确定性（自由度）的数量，通常假定各不同岩石块体为刚性，结构产生的总位移仅仅是由各接触点（面）的变形所引起。
这里的研究对象被认为是各种离散块体的堆砌，块体之间的相互作用力可根据位移和力的关系式来求解，单个块体的运动遵循牛顿运动定律，即力和力矩的平衡。
数值分析模型的建立必须满足平衡方程、变形协调方程和本构方程，此外，还需要满足一定的边界条件。
但离散元块体之间不存在变形协调的约束，因为块体之间是彼此互不约束的，因而仅需满足物理方程和运动方程。
UDEC模拟
5.3.1 采高 1m 时，采空区下坚硬顶板覆岩运动特征 在采深为 400m 的条件下，采高为 1m 时分别对采空区下部煤层采动引起的塑性区、
垂直应力和采动位移等覆岩运动参数的变化情况进行了详细的数值模拟计算，随着工作 面的不断推进，工作面在模拟开挖第 12 步（L=120m）时出现初次来压，之后在 160m、
工作面推进 100m 时，出现了大面积的塑性和拉伸破坏单元，工作面初次来压。当工作
面继续推进到 200m 时，顶板岩层在工作面随后的推进过程中呈现出缓慢的整体弯曲下
沉，并根据三次周期来压时工作面塑性区的分布统计可以发现，采空区上部的采动影响
形成的塑性区和拉伸破坏单元在空间上呈现层次性。
在工作面开采初期，下煤层工作面的应力集中依旧不明显，当工作面初次来压时，
2730
2.34
1.76
2.4
2620
1.93
1.45
2.2
2560
1.74
1.25
1.8
2500
1.56
1.12
1.48
1350
1.05
0.65
1.2
2730
2.34
1.76
2.4
2620
1.93
1.45
2.2
2560
1.74
1.25
1.8
1350
1.05
0.65
1.2
2500
1.56
1.12
1.48
(a) L=20m 时塑性区图
(b) L=60m 时塑性区图
(c) L=100m 时塑性区图
(d) L=200m 时塑性区图
基于UDEC数值模拟的大断面回采巷道数值模拟分析
基于 UDEC 数值模拟的大断面回采巷道数值模拟分析摘要：针对大采高工作面大断面回采巷道围岩稳定性的问题，结合岩石力学、结构力学、弹性力学、矿山压力等理论，对大断面回采巷道围岩变形机理进行分析，结合该煤矿15号主采煤层S8310大采高综采工作面的地质特征和生产技术条件，系统分析了大采高综采大断面回采巷道变形破坏的基本形式，并分析了影响大采高综采大断面回采巷道变形破坏的各种因素，运用UDEC数值模拟软件模拟分析锚杆不同参数对大采高综采大断面回采巷道支护的影响，提出现场回采巷道支护参数及措施。
关键词:大断面；回采巷道；支护参数；厚煤层；数值模拟中图分类号：TD 823文献标识码：A1.工程概况某矿地理位置优越，交通便利。
井田走向长9.6km，倾斜长5km，面积约48km2，15号煤层位于太原组下部，是组内及区内最主要的煤层，煤厚最大值为9.03m，最小值为4.77m，平均6.91m，属稳定煤层，一般含夹石1~4层，夹石在0.1~1.0m之间，属简单结构煤层， 15号煤层S8310工作面走向长1023m，工作面倾斜长度为220m。
工作面地质储量210万吨，可采储量182万吨。
该工作面内15号煤层总厚度最大6.60m，最小6.35m，平均厚度6.51m。
2.数值模拟结果分析（1）模拟方案1中，固定锚杆间排距800mm×900mm，锚杆直径20mm，通过对不同的锚杆长度进行分析，确定合理的锚杆长度。
随着锚杆长度的增加，巷道围岩塑性区范围减小，塑性区单元区域减少，塑性区单元由2.0m时的328减少到2.6m时的304。
表1方案1中位移值比较方案左帮移近量/mm右帮移近量/mm顶板下沉量/mm底板鼓起量/mm1-119319681.42811-217517876.72601-316817073.72511-416516670.9246由表1可知，随着锚杆长度增加，巷道围岩变形量逐渐减小，当锚杆长度为2.4m时，再增加锚杆长度，围岩变化量基本不变，因此，锚杆长度选为2.40m。
UDEC 数值模拟(入门学习)
UDEC 入门;new 是刷新udec窗口，从新调用一个程序;title 与heading代表标题，后面紧跟标题的名称。
如：titlehang dao mo ni;round 指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一。
如：round 0.05set ovtol=0.5;此命令是指层与层之间的嵌入厚度block x1,y1 x2,y2 x3,y3 x4,y4;建立模型框架，crack x1,y1 x2,y2;两点划一线jregion id n x1,y1 x2,y2 x3,y3 x4,y4 deletejset 90,0 4,0 4,0 6,0 0,-50 range jreg 3;jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n;其中xm,ym为起始点坐标,n为设置的区域标号gen quad 10 range xl xu yl yu;在指定的区域生成一定宽度的单元（xu为x方向的取值）zone model mo range xl xu yl yu;使指定的区域材料采用摩尔--库仑本构关系计算（即弹塑性）change jcons=2 range xl xu yl yu;使指定的区域节理遵循摩尔--库仑准则计算（即弹塑性）change mat=1 range xl xu yl yuchange mat=2 range xl xu yl yuchange mat=3 range xl xu yl yu;指定各岩层的材料标号change jmat=1 range xl xu yl yuchange jmat=2 range xl xu yl yuchange jmat=3 range xl xu yl yu;指定各岩层的节理标号prop mat=1 dens=2000prop mat=2 dens=2650prop mat=3 dens=2700;指定各材料的密度，比如1号材料dens=2000，即1立方米重2吨zone k=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6 range mat=1zone k=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6 range mat=2zone k=6.9e9,g=6.6e9,fric=38.62,coh=5.63e6,ten=3.20e6 range mat=3;k为材料的法向刚度，g为材料的切向刚度，friction为材料的内摩擦角，;cohesion为材料的内聚力，tension为材料的抗拉强度prop jmat=1 jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0prop jmat=2 jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0prop jmat=3 jkn=20e8,jks=16e7,jcoh=0.4e6,jfric=15,jten=0;jkn为节理的法向刚度，jks为节理的切向刚度，jfriction为节理的内摩擦角，;jcohesion为节理的内聚力，jtension为节理的抗拉强度set gravity 0,-9.81;设置重力加速度，x方向为0，y方向为-9.8bound xvel=0 range -0.1 3.00 -60.1 20.1bound xvel=0 range 97 100.1 -60.1 20.1bound yvel=0 range 0.1 100.1 -60.1 -58;采用位移法固定边界solve\step 5000\cycle 5000;执行计算save pingheng.sav;保存文件，文件的后缀为.sav，文件名可以自己命名。
地震作用对采空区塌陷的UDEC模拟
地震作用对采空区塌陷的UDEC模拟500kV输电线路、59 条 220kV输电线运 4 条前言随着我国西电东送战略的实施，输础的破坏十分严重，而由于陡峭的地形、区域地路、122条 110kV输电线路。
严重损毁变电站 16 电线路质构造、岩体风化破碎等引发的滑坡、崩塌及泥的电压等级越来越高,超高压、特高压交、直流输座，其中 500kV1座、220kV3座、110kV5座。
四石流等地震地质次生灾害对电网工程的破坏更川全网损失负荷近 400 万 kW，负荷损失率为电已成为主电网调配区域能源的主要手段。
由为严重。
于我国西部地区受到印度板块向北推移挤压， 31.8,，6 市州负荷损失率达75.7,。
在复杂的区域地质构造及地质地震背景青藏高原强烈变形，高原内部及其边缘的活断重灾区主要受损的线路有:220kV 平回线、下，线路路径及塔位选择时，如何合理利用地形 220kV福回一线、220kV 福回二线、220kV 耿山层上经常发生强烈地震，我国西部地区已经是地质条件并考虑地震地质作用，将是线路勘测世界大陆内部活跃的强烈地震区，因此很多西南(北)线、220kV 映山线、220kV 渔山东(西)线、设计中地质工作者的重要任务。
220kV平山线、220kV源山南(北)线、220kV 丹 1,汶川大地震对震区输电线路的破坏现状部山区输电线路将从强震区穿过，跨越可能的 2008年 5 月 12 日 14 点 28 分，汶川发生山线、110kV 映灌线、220kV 福银线、220kV 竹发震断层不可避免。
确保西部山区输电线路的了举世震惊的里氏 8.0 级特大地震。
地震不但给茂线、220kV 红薛线、220kV 薛州线、220kV 州地震安全十分重要，不仅可以避免造成巨额的经济损失,还可以确保其他基础设施如通讯、交茂线共 5 条线路、500kV 茂谭一二线及220kV 上百万家庭带来巨大灾难，也摧毁了大量的电通、供水等的正常运行,避免次生灾害的发生,并茂永线等。
UDEC 数值模拟(入门学习)
UDEC 入门;new 是刷新udec窗口，从新调用一个程序;title 与heading代表标题，后面紧跟标题的名称。
如：titlehang dao mo ni;round 指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一。
如：round 0.05set ovtol=0.5;此命令是指层与层之间的嵌入厚度block x1,y1 x2,y2 x3,y3 x4,y4;建立模型框架，crack x1,y1 x2,y2;两点划一线jregion id n x1,y1 x2,y2 x3,y3 x4,y4 deletejset 90,0 4,0 4,0 6,0 0,-50 range jreg 3;jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n;其中xm,ym为起始点坐标,n为设置的区域标号gen quad 10 range xl xu yl yu;在指定的区域生成一定宽度的单元（xu为x方向的取值）zone model mo range xl xu yl yu;使指定的区域材料采用摩尔--库仑本构关系计算（即弹塑性）change jcons=2 range xl xu yl yu;使指定的区域节理遵循摩尔--库仑准则计算（即弹塑性）change mat=1 range xl xu yl yuchange mat=2 range xl xu yl yuchange mat=3 range xl xu yl yu;指定各岩层的材料标号change jmat=1 range xl xu yl yuchange jmat=2 range xl xu yl yuchange jmat=3 range xl xu yl yu;指定各岩层的节理标号prop mat=1 dens=2000prop mat=2 dens=2650prop mat=3 dens=2700;指定各材料的密度，比如1号材料dens=2000，即1立方米重2吨zone k=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6 range mat=1zone k=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6 range mat=2zone k=6.9e9,g=6.6e9,fric=38.62,coh=5.63e6,ten=3.20e6 range mat=3;k为材料的法向刚度，g为材料的切向刚度，friction为材料的内摩擦角，;cohesion为材料的内聚力，tension为材料的抗拉强度prop jmat=1 jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0prop jmat=2 jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0prop jmat=3 jkn=20e8,jks=16e7,jcoh=0.4e6,jfric=15,jten=0;jkn为节理的法向刚度，jks为节理的切向刚度，jfriction为节理的内摩擦角，;jcohesion为节理的内聚力，jtension为节理的抗拉强度set gravity 0,-9.81;设置重力加速度，x方向为0，y方向为-9.8bound xvel=0 range -0.1 3.00 -60.1 20.1bound xvel=0 range 97 100.1 -60.1 20.1bound yvel=0 range 0.1 100.1 -60.1 -58;采用位移法固定边界solve\step 5000\cycle 5000;执行计算save pingheng.sav;保存文件，文件的后缀为.sav，文件名可以自己命名。
UDEC模拟实例与解析
^`UDEC 实例翻译与命令解析翻译：珠穆朗玛1 地震诱发地层坍塌 Seismic-Induced Groundfall1.1 问题描述本例展示使用 UDEC 模拟分析地震诱发地层坍塌的一类的问题,模型见图 1.1,该模型基 于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿 34-1-554 切割断面的一个剖面图的结构和 尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高 5m,宽 10m.假定两个连续节理交叉平面分析:一个角度为 45 度,另一个为-9 度,两者节理间距均为 5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳 定性。
围岩参数来自试验室平均测试数值，假定岩石块体参数如下：假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节 理参数，如下：初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m 的岩层产生）。
1.2 UDEC 分析UDEC 模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1.2.超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型顶部y 方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加1.25Mpa 的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是1.25 MPa/2, or 0.625 MPa.0.02 秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1 点位于开挖面的左角,点2 位于拱顶块体的右角, 图1.4 的位移时间曲线显示两点本质上是弹性反应.本例关心的问题是在模型顶部施加的速度和计算速度的对比,下面的公式可以用以估计施加的波速.使用这个方程,施加的最大波速大概是0.04m/sec,图1.5 显示的峰值波速小于0.06m/sec. 估计的波速和监测波速的不同在于使用的围岩模量.而是没有考虑节理变形的相等变形模量.在第二个案例中,施加应力波峰值12.5 Mpa(有效应力6.25Mpa).0.02 秒后的开挖拱顶应力分布见图1.6.该图显示出拱顶岩体不受力,表面该块体已经松散并正在下落.对于关心的问题,后来三个时间的几何体和应力分布见图1.8 至图1.10.在问题的顶部预测的波速(从上面的方程)是0.4m/sec.从模型中计算的波速见图1.11,再次,由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.1.3 节包含了该模型的数据列表,该列表包含了一个FISH 函数(show)被用来创建坍塌的动画文件,每隔0.02 秒俘获一个显示的图片.通过改变FISH 参数time_int 可以改变动画帧的间隔.视图的总数也可以通过改变snap_shot 的数值进行改变.为了显示80 帧的显示图片而创建的该电影文件需要大概13MB 的硬盘空间.1.3 数据文件列表Example 1.1 SEISMIC.DATtitleSEISMIC INDUCED ROOF COLLAPSE 地震诱发拱顶坍塌;round 0.01; define original boundary of modeled region 定义模型区域的原始边界block -25,-20 -25,20 25,20 25,-20; generate joint pattern over entire original region 在整个原始区域生成节理形态jregion id 1 -25,-25 -25,25 25,25 25,-25jset 45,0 200,0 0,0 5.0,0 (0,0) range jreg 1jset -9,0 200,0 0,0 5.0,0 (0,0) range jreg 1; put in joints needed for the later excavation 为了后面开挖而设置的节理crack -5.01,-2.51 5.01,-2.51crack -5.01, 2.51 5.01, 2.51crack -5,-2.5 -5,2.5crack 5,-2.5 5,2.5crack 2.25,2.5 1.93,5.0; generate fdef zones and assign joint properties (mat=1 & jmat=1;default) 生成单元和设置节理参数generate edge 9.0 range -30,30 -30,30prop mat=1 d=0.00300 k=39060 g=31780prop jmat=1 jkn=20000 jks=20000prop jmat=1 jf=30.0; apply boundary conditions and initial conditions to 在地应力下施加边界条件和初始条件; consolidate model under field stressesbound stress=-24.0, 0.0, -24.0 ygrad=-0.3 0 -0.3insitu stress=-24.0, 0.0, -24.0 ygrad=-0.3 0 -0.3bound yvel 0.0 range -26,26 -21,-19grav 0.0 -10.0; track the x-displacement, and y-displacement over time 追踪位移hist solvehist xdis=0,7 ydis=0,7 type 1solve rat 1e-5; save consolidated statesave seismic1.sav; make excavationdelete range -5,5 -2.5,2.5solve rat 1e-5; save excavated statesave seismic2.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.04 m/sec);; set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -1.25 yhist=cos(100.0,0.0195) range -26 26 19 21 ;reset time hist disp rothist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 mass; 0.02 sec.cyc time 0.02save seismic3.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.4 m/sec); set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -12.5 yhist=cos(100.0,0.0195) range -26 26 19 21 reset time hist disphist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 masssave seismov.sav;; 0.02 sec.cyc time 0.02save seismic4.sav; 0.25 sec.cyc time 0.23save seismic5.sav; 0.50 sec.cyc time 0.25save seismic6.sav; 0.75 seccyc time 0.25save seismic7.sav;rest seismov.sav; make a movie of the groundfall;wind -12 12 -12 12set ovtol 0.05plot block vel max 2.0 blue stress max 50movie onmovie file = seismic.dcxmovie step 1000step 400003 隧道支护荷载Tunnel Support Loading3.1 问题陈述本例模拟展示了UDEC 在检查衬砌隧道方面的应用,着重强调了荷载在混凝土衬砌中的发展,本例也解释了模拟连续建造操作中独立阶段的模拟程序.隧道系统的理想几何体见图3.1.系统包含在海床下大约70m(中线)深度,中线间距12m 的两个隧道, 初始水位在隧道中线上方110m 处.服务隧道直径5.24m,衬砌厚度37cm.主隧道直径8.22m,衬砌厚度46cm.服务隧道先于主隧道开挖和衬砌.随后设置主隧道衬砌,水位上升增加到100m.施工顺序是:(1)开挖服务隧道excavation of the service tunnel;(2)衬砌服务隧道lining of the service tunnel; (3)开挖主隧道excavation of the main tunnel; (4)衬砌主隧道lining of the main tunnel; and (5)升高水位raising of the water level.分析的目的是评价每个施工阶段服务隧道和主隧道支护状况.本例的材料参数见下:岩体——开挖隧道的围岩参数为：弹性模量elastic modulus 0.89 GPa泊松比Poisson’s ratio 0.35单轴抗压强度uniaxial compressive strength 3.5 MPa粘聚力cohesion 1 MPa密度density 1340 kg/m3混凝土衬砌——弹性模量为24 GPa ，泊松比为0.19. 假定衬砌为线弹性材料。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2560
1.74
1.25
1.8
2730
2.34
1.76
2.4
内摩擦角（˚）
41 37 35 33 27 41 37 35 27 33 35 41
5.2.3 对比方案及模拟步骤 为了研究采空区下坚硬顶板顶分层开采时，工作面的应力分布及覆岩运动规律，采
用 UDEC4.0 对神华宁煤汝箕沟煤矿 3213 工作面进行数值模拟，研究下煤层顶分层开采 活动对于采动应力分布的影响。主要研究下煤层采高、开采深度、煤层间距以及上部采 空区保护煤柱等因素对工作面附近煤体以及工作面上覆岩层的塑性区、垂直应力和采动 位移的分布情况进行分析。模拟方案如下：
2730
2.34
1.76
2.4
2620
1.93
1.45
Hale Waihona Puke 2.225601.74
1.25
1.8
2500
1.56
1.12
1.48
1350
1.05
0.65
1.2
2730
2.34
1.76
2.4
2620
1.93
1.45
2.2
2560
1.74
1.25
1.8
1350
1.05
0.65
1.2
2500
1.56
1.12
1.48
5.3 采高对采空区下坚硬顶板顶分层开采覆岩运动的影响
5.3.1 采高 1m 时，采空区下坚硬顶板覆岩运动特征 在采深为 400m 的条件下，采高为 1m 时分别对采空区下部煤层采动引起的塑性区、
垂直应力和采动位移等覆岩运动参数的变化情况进行了详细的数值模拟计算，随着工作 面的不断推进，工作面在模拟开挖第 12 步（L=120m）时出现初次来压，之后在 160m、
5.1 UDEC2D 简介
UDEC2D 是一个利用离散元来模拟不连续介质（存在节理裂隙的岩体等）在承载静 压载荷或者动压载荷的作用下，岩体变形、破坏或者运动的程序。岩体被划分成多个块 体，变形块体又被划分成有限个单元网格，且每一单元根据给定的“应力－应变”准则， 力学形式表现为线性或非线性特性。在 UDEC 中，非连续介质是通过离散的块体集合体 加以表示，不连续面发生法向和切向的相对运动也由线性或非线性“力－位移”的关系控 制。
综合柱状图，数值模型结构如图 5.1 示。
图 5.1 数值模拟模型 Fig.5.1 Numerical simulation model
5.2.2 岩层属性及边界条件 根据神华宁煤汝箕沟煤矿煤层的地质柱状图，取埋深为 400m，按平均密度为
2.7g/cm3 计算出来的重力 γH，模型上部边界施加 10.8Mpa 的均布载荷，垂直加速度为 9.8m/s2。模拟采用位移控制边界条件的方法，左右边界的 x 方向速度为零，固定底部边 界，上边界作为自由承载面。岩石力学参数见表 5.1 示。
5.4.2 采高 3m 时，采空区下坚硬顶板覆岩运动特征
当采高为 3m 时，通过数值模拟表明，工作面在模拟开挖第 10 步（L=100m）时出 现初次来压，初次来压步距为 100m，并分别在 L=140m、170m、200m 出现周期来压， 平均周期来压步距 33m。与采高为 1m 时相比单次周期来压在时间上提前近一个开挖步 距，平均周期来压步距缩短了 7m。
位移相比采高 1m 时向上扩展了 12m，大位移区域面积增大。当 L=200m 时，采动影响
引起的上覆岩层垂直位移沿走向方向进一步的扩展，但在垂直方向上不再有明显的变
化，不再向上扩展，“三带”高度较采高 1m 时扩展了 16m，最终“三带”扩展高度为
50m。
5.4.3 采高 5m 时，采空区下坚硬顶板覆岩运动特征
200m 分别出现周期来压，初次来压步距为 120m，周期来压步距为 40m。下面就工作面 推进过程中各项应力场参数分布特征作详细分析。
(a) L=40m 时塑性区图
(b) L=80m 时塑性区图
(c) L=120m 时塑性区图
(d) L=200m 时塑性区图
(e L=40m 时垂直应力图
(f) L=80m 时垂直应力图
5 采空区下厚煤层顶分层开采坚硬顶板运动规律数值模拟分析
在厚煤层分层开采过程中，尤其是上部存在采空区的条件下，随着工作面不断地推 进上覆岩层的运移在形态和结构上呈现出非线性的特征，其运动演化规律与上部不存在 采空区的坚硬顶板条件下表现出强烈的不同，上部采空区以及内部存在的保护煤柱影响 着下煤层回采工作面的应力分布及顶板运移规律。所以本章采用二维离散元法进一步模 拟采空区下坚硬顶板顶分层在正常推进过程中的岩层运动规律及工作面应力分布特征， 分析上部采空区对下煤层顶分层开采的垮落方式、覆岩运动及应力分布的影响。
工作面前方形成了“弱应力集中区”，但与采高 1m 时相比起应力集中的位置更靠近煤壁，
应力集中的范围更广。L=200m 时，工作面前方相当于“强应力集中区”，由于采动影响
引起的上覆岩层的应力分布更加集中，叠加效果更加明显。
从垂直位移云图中可以看出，当 L=100m 时，由于采动影响引起的上覆岩层的垂直
当采高为 5m 时，通过数值模拟表明，工作面在模拟开挖第 9 步（L=90m）时出现 初次来压，初次来压步距为 90m，并分别在 L=120m、150m、170m、200m 出现周期来 压，平均周期来压步距 28m。与采高为 3m 时相比单次周期来压在时间上提前了半个开 挖步距，平均周期来压步距缩短了 5m。
4、上部采空区保护煤柱对采空区下坚硬顶板顶分层开采的影响。上煤层在开采过 后会留下一定的保护煤柱，根据上煤层保护煤柱与下煤层工作面开切眼处的距离，在初 始模型不变的情况下，变换上煤层残留煤柱与下煤层开切眼的距离为 20m、60m、100m， 进行数值模拟研究。
模型开挖步骤：模拟过程中，先开采上煤层，模拟上煤层回采工作面每次的开挖步 距 ΔL=20m，每次要进行单独的平衡计算，模型在应力平衡后再进行下一步的开挖，上 煤层工作面连续推进∑L=240m，当上煤层开采完之后，再模拟下煤层的开挖过程，模拟 下煤层回采工作面每次的开挖步距 ΔL=10m，每次要进行单独的平衡计算，模型在应力 平衡后再进行下一步的开挖，下煤层工作面连续推进∑L=200m，由于文章研究的重点是 下煤层的上覆岩层运动规律及应力分布特征，因此数值模型在开挖上煤层时左右边界煤 柱留设 80m，开采下煤层研究工作面前方的采动应力及运动步距，因此下煤层在开挖时 左右边界煤柱留设 100m。
5.2 数值模拟方案及参数选取
5.2.1 数值模型 本次模拟是以神华宁煤汝箕沟煤矿 3213 工作面的工程地质条件为基础。模型尺寸
400*300m，上煤层顶板距上部边界 150m，上煤层厚度 6m，下煤层顶板距上部边界 216m， 煤层厚度 9m，下煤层分层开采顶分层开采厚度为 3m，底板岩层厚度取 75m。根据地层
工作面推进 100m 时，出现了大面积的塑性和拉伸破坏单元，工作面初次来压。当工作
面继续推进到 200m 时，顶板岩层在工作面随后的推进过程中呈现出缓慢的整体弯曲下
沉，并根据三次周期来压时工作面塑性区的分布统计可以发现，采空区上部的采动影响
形成的塑性区和拉伸破坏单元在空间上呈现层次性。
在工作面开采初期，下煤层工作面的应力集中依旧不明显，当工作面初次来压时，
1、下煤层开采高度对于采空区下坚硬顶板顶分层开采的影响。在采深为 400m 的条 件下，变换采高为 1m、3m、5m，进行数值模拟研究；
2、开采深度对于采空区下坚硬顶板顶分层开采的影响。在采高为 3m 条件下，变换 初始模型中的采深分别为 400m、600m、800m，进行值模型研究；
3、煤层间距对采空区下坚硬顶板顶分层开采的影响。在采高为 3m 和采深为 400m 的条件下，变换初始模型中的煤层间距为 40m、60m、100m，进行数值模拟研究；
(g) L=120m 时垂直应力图
(h) L=200m 时垂直应力图
(i) L=40m 垂直位移云图
(j) L=80m 垂直位移云图
(k) L=120m 垂直位移云图
(l) L=200m 垂直位移云图
图 5.2 采高 M=1m 时数值模拟结果
Fig.5.2 Mining height M=1m, numerical simulation results
而导致顶板的断裂，出现了工作面初次来压。工作面推进到 200m 时，顶板周期来压，
上部覆岩出现了较大面积的拉伸破坏，采动影响范围基本稳定。
对垂直应力的分析主要在于应力集中区的形成范围及发展变化。L=40m 时，开切眼 和煤壁处的应力分布较均匀，并未出现应力集中的现象，这是由于上部采空区的卸压作 用使得下煤层工作面的应力较小。L=120m 时，应力集中区并未随着工作面悬顶距的不 断增大而激增，反而在超前工作面的上方岩层中应力较低，在工作项目煤壁前方 40m 处 出现了较不明显的应力集中区，本文称之为“弱应力集中区”。其特点是：位于采空区 下，应力相对于周围存在集中的现象，应力低于原岩应力，这是由于上部采空区的卸压 作用。L=200m 时，下煤层工作面靠近到上部采空区边界，在上煤层应力向下传播过程 中与下煤层边界应力叠加，导致了在下煤层工作面前方出现了较强的应力集中区，本文 称之为“强应力集中区”。 其特点是：应力集中系数很大（可达到 6.5），应力集中范围 较小（仅仅有 6m 左右），属于边界煤柱应力相互的叠加的结果。
岩性
粗砂岩 中砂岩 细砂岩 粉砂岩 上煤层 粗砂岩 中砂岩 细砂岩 下煤层 粉砂岩 细砂岩 粗砂岩
厚度 (m) 60 50 24 16 6 16 20 16 9 8 15 60
表 5.1 模型岩石属性
Table.5.1 Model of rock properties
密度
体积模量 剪切模量
内聚力