UDEC模拟实例与解析Word版.pdf

传播优秀Word版文档，希望对您有帮助，可双击去除！

UDEC实例翻译与命令解析

翻译：珠穆朗玛

传播优秀Word版文档，希望对您有帮助，可双击去除！

1地震诱发地层坍塌Seismic-InducedGroundfall

1.1问题描述

本例展示使用UDEC模拟分析地震诱发地层坍塌的一类的问题,模型见图1.1,该模型基

于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿34-1-554切割断面的一个剖面图的结构和

尺寸.用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高5m,宽10m.

假定两个连续节理交叉平面分析:一个角度为45度,另一个为-9度,两者节理间距均为

5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳

定性。围岩参数来自试验室平均测试数值，假定岩石块体参数如下：

假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节

理参数，如下：

传播优秀Word版文档，希望对您有帮助，可双击去除！

初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m的岩层

产生）。

1.2UDEC分析

UDEC模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开

挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1.2.超采正上方和下方的块

体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题

域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型

顶部y方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加

1.25Mpa的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部,施加的有效影响

应力应该是1.25MPa/2,or0.625MPa.0.02秒后的开挖面拱顶的应力分布见图1.3,两点的

位移被监测,1点位于开挖面的左角,点2位于拱顶块体的右角,图1.4的位移时间曲线显示两

点本质上是弹性反应

.

传播优秀Word版文档，希望对您有帮助，可双击去除！

传播优秀Word版文档，希望对您有帮助，可双击去除！

本例关心的问题是在模型顶部施加的速度和计算速度的对比,下面的公式可以用以估计

施加的波速.

使用这个方程,施加的最大波速大概是0.04m/sec,图1.5显示的峰值波速小于0.06m/sec.

估计的波速和监测波速的不同在于使用的围岩模量.而是没有考虑节理变形的相等变形模量.

在第二个案例中,施加应力波峰值12.5Mpa(有效应力6.25Mpa).0.02秒后的开挖拱顶应

力分布见图1.6.该图显示出拱顶岩体不受力,表面该块体已经松散并正在下落.对于关心的问

题,后来三个时间的几何体和应力分布见图1.8至图1.10.

在问题的顶部预测的波速(从上面的方程)是0.4m/sec.从模型中计算的波速见图1.11,再

次,由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.

1.3节包含了该模型的数据列表,该列表包含了一个FISH函数(show)被用来创建坍塌的

动画文件,每隔0.02秒俘获一个显示的图片.通过改变FISH参数time_int可以改变动画帧的

传播优秀Word版文档，希望对您有帮助，可双击去除！

间隔.视图的总数也可以通过改变snap_shot的数值进行改变.为了显示80帧的显示图片而

创建的该电影文件需要大概13MB的硬盘空间.

传播优秀Word版文档，希望对您有帮助，可双击去除！

传播优秀Word版文档，希望对您有帮助，可双击去除！

传播优秀Word版文档，希望对您有帮助，可双击去除！

1.3数据文件列表