金山店铁矿地下开采引起地表变形规律的离散元模拟研究.pdfVIP

第21卷 第8期 岩石力学与工程学报 21(8)：1130～1135 2002年 8月 Chin Aug．，2002 ． ese Journal ofRock Mechanics and Engineering 金山店铁矿地下开采引起地表变形 规律的离散元模拟研究 郑榕明 陈文胜 葛修润 冯夏庭 (’香港理工大学土木与结构工程系 香港) ( 中国科学院武汉岩土力学研究所 武汉 430071) 摘要 用离散单元法对金山店矿山进行了二维模拟研究，得到了其围岩变形规律，且计算结果为圈定地表的移动 和陷落范围提供了可靠的依据：同时，对离散元方法的应用，特别是参数选取进行了有益的探讨。 关键词 离散单元法，地表变形，地下采矿 分类号 TD 31，O 242 文献标识码 A 文章编号 1000．6915(2002)08．1 130．06 发出的最新的离散元分析程序，2D．Block共包括前 1 前 言 处理器、求解器和后处理器3部分。用户只需将设 计图、材料常数、载荷、边界条件等参数以及填土、 因地下开挖而引起地表变形甚至塌陷是一个复 挖掘等与施工工艺相关的数据直接输入给前处理 杂的岩体大变形和大位移问题Il 】。武钢金山店铁 器，就可从后处理器获得块体接触关系图、块体情 矿是武钢的主要矿石基地之一，随着地下采矿深度 况图等易于理解的直观结果。用户不必关心与离散 的不断增加和采空区的不断扩大，由此引起的对地 元理论相关的中间处理过程Io】。 表的影响问题己成为矿山能否继续生存及生产经济 效益的关键问题。 2 计算模型和计算参数 对于该矿山出现的此类地表问题，用基于连续 性假设的有限元等方法研究己很难满足矿山生产和 2．1计算模型 设计的要求。由于离散单元法特别适合于因岩体塌 本次研究工作选择了3个计算剖面，计算剖面 陷而出现的大位移这种岩石力学问题I3~51，故本文 所取的范围考虑到为了满足计算要求，计算范围应 采用该方法对金山店矿山进行模拟计算。 超出预估的开采影响范围。一般可按采矿设计的错 离散单元法的计算原理虽然很简单，但在计算 动角、陷落角预估影响范围I7】，即从开采到的最低 机上实施起来却非常复杂，涉及到很多问题。离散 高程按错动角延伸至地表，再作适当扩延即可作为 单元法可采用动态松弛法求解，其基本运动方程为 横向的计算范围；纵向范围根据开采深度，取从地 mii(t)+ca(t)+ku(t)= (，) (1) 表至最低开采高程并略作延伸。最终水平方向长度 式中：m，c，k，厂分别为质量、阻尼系数、弹性系 为 l 200 m，垂直方向为从地表至V一600 m高程。 数及载荷函数；U为块体位移。 由于计算工作量大，根据经验采用了下窄上宽的边 动态松弛法对上式用中心差分法求解，是一种 界。被减掉的块体是基本不发生位移的块体，对最 显式求解，缺点在于计算时步△f要求很小，且需要 终计算结果影响小。本文取其中的 34 剖面为例， 合理地确定计算系数等。在算法中，假定力在一个 其岩性分区及离散元块体模型见图l。 时步内只能传递到一