UDEC 实例翻译与命令解析.doc

UDEC 实例翻译与命令解析
翻译:珠穆朗玛
1 地震诱发地层坍塌 Seismic-Induced Groundfall
问题描述
本例展示使用 UDEC 模拟分析地震诱发地层坍塌的一类的问题,模型见图 ,该模型基于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿 34-1-554 切割断面的一个剖面图的结构和尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高 5m,宽 10m.
假定两个连续节理交叉平面分析:一个角度为 45 度,另一个为-9 度,两者节理间距均为
5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳定性。围岩参数来自试验室平均测试数值,假定岩石块体参数如下:
假定块体仅具有弹性行为,节理假定符合库伦滑动准则,选择典型的教课书数值作为节理参数,如下:
初始应力状态按各向同性估计为 24Mpa(假定垂直荷载由覆盖深度大约 800m 的岩层
产生)。
UDEC 分析
UDEC 模拟顺序分三个阶段,首先,, ,,地震事件用施加到模型顶部 y ,施加
的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是 MPa/2, or 秒后的开挖面拱顶的应力分布见图 ,两点的位移被监测,1 点位于开挖面的左角,点 2 位于拱顶块体的右角, 图 的位移时间曲线显示两点本质上是弹性反应
.
本例关心的问题是在模型顶部施加的速度和计算速度的对比,下面的公式可以用以估计
施加的波速.
使用这个方程,施加的最大波速大概是 ,图 显示的峰值波速小于 .
.
在第二个案例中,施加应力波峰值 Mpa(有效应力 ). 秒后的开挖拱顶应力分布见图 ,,后来三个时间的几何体和应力分布见图 至图 .
在问题的顶部预测的波速(从上面的方程)是 ,再次,由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.
节包含了该模型的数据列表,该列表包含了一个 FISH 函数(show)被用来创建坍塌的动画文件,每隔 FISH 参数 time_int 可以改变动画帧的
snap_shot 80 帧的显示图片而
创建的该电影文件需要大概 13MB 的硬盘空间.
数据文件列表
title
Example
SEISMIC INDUCED ROOF COLLAPSE 地震诱发拱顶坍塌
;
round
; 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,-25 jset 45,0 200,0 0,0 ,0 (0,0) range jreg 1 jset -9,0 200,0 0,0 ,0 (0,0) range jreg 1
; put in joints needed for the later excavation 为了后面开挖而设置的节理
crack -,- ,- crack -, , crack -5,- -5,
crack 5,- 5,
crack , ,
; generate fdef zones and assign joint properties (m