- 1、本文档共84页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
UDEC模拟实例及解析
UDEC 实例翻译与命令解析
翻译:珠穆朗玛
1 地震诱发地层坍塌 Seismic-Induced Groundfall
1.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 数据文件列表
title
Example 1.1 SEISMIC.DAT
SEISMIC 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,-25 jset 45,0 200,0 0,0 5.0,0 (0,0) range jreg 1 jset -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.51 crack -5.01, 2.51 5.01, 2.51 crack -5,-2.5 -5,2.5
crack 5,-2.5 5,2.5
crack 2.25,2.5 1.93,5.0
; generate fdef zones and assign joint properties (mat=1
您可能关注的文档
- TDD-LTE 干扰分析及排查.docx
- TEACH-急性冠脉综合征抗血小板治疗.pptx
- TDD-LTE上行干扰定位方法及排查指导手册.doc
- TD-SCDMAMR定位试验总结报告.doc
- Topic-15_Image_Coding-I 及图像处理有关的代码matlab.pptx
- TPM工具11-设备完好及机台状态管理140218.ppt
- TMT红外热扫简介及中医体质辨识.ppt
- TSI探头工作原理及安装.ppt
- Tuxedo及Weblogic互连指南.docx
- ug 曲线编辑操作.ppt
- 国开景区管理作业2试题及答案.pdf
- 国开景区管理作业1-4试题及答案.pdf
- 河南开放大学本科《地域文化(本)》作业练习1-3试题及答案.pdf
- 2024年大型游乐设施操作证考试题库及答案很全.pdf
- 2024年门座式起重机司机考试题库及答案.pdf
- 2022-2023学年河北省衡水市武强中学高二(下)期末数学试卷【答案版】.docx
- 2022-2023学年河北省保定市崇德实验中学高二(下)期末数学试卷【答案版】.docx
- 江西省2017年中小学教师招聘考试高中化学试卷及答案.docx
- 2024年河北省八年级中考生物真题(解析版).docx
- 2024年南阳市社会保险中心(唐河县企业养老保险分中心)(参公)一级科员招录1人《行政职业能力测验》高频考点、难点(答案详解版).docx
文档评论(0)