随着计算机技术的高速发展, 数值模拟技术成为分析岩土工程问题的重要手段[1].但在很多岩土工程数值模拟中, 建模等前处理工作复杂耗时, 特别是面对复杂不规则的地质体, 常规的前处理方法往往无能为力, 这影响了数值模拟技术在岩土工程领域的应用.
RFPA3D是在RFPA2D的基础上扩展开发的一种基于有限元应力分析和统计损伤理论的材料破裂过程分析软件, 能够模拟岩土工程中岩质材料渐进破裂直至失稳的过程[2-4].文献[5-7]分别用RFPA3D软件对云阳县云阳大道边坡和锦屏水电站左岸岩质边坡稳定性进行模拟, 成功得到两种边坡失稳破坏过程.但RFPA3D软件前处理功能薄弱, 自身无法建立复杂三维模型, 同时其单元必须为六面体单元, 这又限制了很多前处理软件在RFPA3D模型建立中的应用.
SURPAC是一款先进的矿业软件, 具有强大的复杂地质体精确建模能力, 可以准确表现出不同地质结构体及其在空间中的分布, 还可建立六面体单元的块体模型[8].将SURPAC和RFPA3D联合应用, 可以各取所长, 实现对复杂地质体的RFPA3D建模.类似的SURPAC与岩土分析软件的联合应用已经得到实现, 林杭等[9]采用Fortran语言编制SURPAC-FLAC3D接口程序, 实现模型转换; Li等[10]通过MATLAB软件编制接口程序, 将SURPAC模型成功导入FLAC3D中.
本文在SURPAC中建立块体模型, 输出块体模型的单元数据, 利用MATLAB软件的m文件编制SURPAC-RFPA3D的模型转换接口STOR.m, 将SURPAC块体模型转换为RFPA3D计算模型.利用该转换方式, 进行某露天铁矿东边坡挂帮矿开采的模拟, 模拟结果与现场观测结果基本吻合, 证明了SURPAC与RFPA3D耦合计算的有效性.
1 SURPAC模型构建SURPAC对于地质体模型有两种表现形式:一种是实体模型, 另一种是块体单元.实体单元可以形象地展示出地质结构体的形态和分布, 例如岩层、空区和断层等, 但是这种实体模型本质仅为一层蒙皮, 其内部没有划分单元, 即内部无有效的单元数据.基于实体模型可以建立一种更准确的三维块体模型, 它包括大量离散的规则六面体单元, 这些块体模型不仅可以准确地建立地质体特性与空间位置的关系, 而且其本质上是一个数据库, 包含顶点坐标、岩性种类等数据, 恰好可以对应RFPA3D建模所要求的六面体单元.
SURPAC构建块体模型的步骤如图 1所示, 利用地质平面图、剖面图和地表等高线图等二维图纸, 建立三维实体模型, 根据实体模型和不同岩性的约束情况建立块体模型, 并为块体赋属性.
将SURPAC块体模型的单元信息导出为EXCEL格式文件, 该文件包含的信息为块体模型每一单元的质心三维坐标、单元三边长度和该单元的属性.
2 SURPAC与RFPA3D模型转换SURPAC中, 块体单元模型的形状为长方体, 这满足RFPA3D所要求的单元必须为六面体的原则.RFPA3D通过4个文本文件接受外部导入的模型, 分别为coor0.txt, elem0.txt, id0.txt, disp0.txt.
要完成模型的转换, 需要对所提取的SURPAC块体模型数据进行处理, 变为RFPA3D可接受的4个文本文件.
coor0.txt文件为RFPA3D模型所有节点的三维空间坐标信息文件.RFPA3D的单元为六面体, 其8个顶点(P1~P8)的编号顺序见图 2.SURPAC块体单元输出的信息为单元质心的坐标(X0, Y0, Z0)、三维尺寸L1, L2, L3和块体单元属性的编号[9].图 2展现了2种模型数据的关系, 根据几何关系得到两者坐标的转换方程:
式中,XP1, YP1, ZP1分别表示点P1的X, Y, Z坐标.
同理可得P2~P8各点的三维坐标.
将所得到的单元顶点数据进行筛选, 删掉重复的顶点数据, 将筛选后的顶点数据按照Gn, Xn, Yn, Zn的格式写入文本coor.txt, n为筛选后的顶点数目, Gn为顶点编号, Xn, Yn, Zn为相应顶点的坐标.
elem0.txt为记录模型所有单元数据的文件, 包括RFPA3D模型中每个单元的顶点编号和该单元的属性编号.在筛选后的顶点中, 查询每个单元包含的顶点编号, 按照Zm, P1, P2, P3, P4, P5, P6, P7, P8, V的格式将顶点编号写入elem0.txt文件, m为单元数目, Zm为单元编号, P1~P8为单元顶点编号(如图 2所示), V为单元属性编号.
id0.txt为模型边界条件的信息文件, 包含每个顶点X, Y, Z三个方向的约束情况.RFPA3D规定, 节点位移约束为-1, 应力约束为1.搜索边界的顶点、区分边界处顶点和模型内部顶点, 按照Gn, Xn, Yn, Zn的格式, 将顶点信息写入id0.txt文件, Xn, Yn, Zn为相应顶点的三个坐标方向的约束条件, 为1或-1.
disp0.txt为模型节点的初始位移信息文件, 包含每个节点X, Y, Z三个方向的位移情况, 一般默认为0.将筛选后的顶点数据按照Gn, Xn, Yn, Zn(Xn, Yn, Zn一般默认为0)的格式写入disp0.txt.
根据上述的SURPAC到RFPA3D模型转换的方法, 在MATLAB中利用m语言编制模型转换程序STOR.m.该程序将SURPAC导出的块体模型单元数据转换为RFPA3D可接受的4个数据文件, 再导入RFPA3D中,分别给不同组岩体赋予力学参数, 即可进行计算.具体转换流程见图 3.
某大型露天铁矿露天生产进入末期, 为延续产能, 利用无底柱分段崩落法开采东部边坡下部的挂帮矿, 由此引起东部边坡失稳破坏, 边坡及地表产生大规模岩移.本文应用上述方法, 在SURPAC中对该露天矿东部边坡进行建模, 将模型导入RFPA3D中模拟计算边坡破坏失稳的情况.
选取该矿东边坡作为模拟区域, 所选区域南北方向长700 m, 东西方向长600 m, 由地表至-200 m水平; 开采区域位于边坡下部矿体中, 垂直范围为-51~-123 m, 水平范围为南北走向100~140 m, 东西走向约130 m.利用该矿的平面图和剖面图, 在SURPAC中建立东边坡的实体模型, 结合约束条件将实体模型转化为包含3组属性的块体模型.如图 4所示, 模型北部为千枚岩, 南部为混合岩, 中部为矿体; 边坡下挂帮矿开采区域见图 5.所建模型单元为尺寸统一的立方体单元, 单元南北走向边长7 m, 东西走向边长6 m, 垂直走向边长5 m, 模型单元共计542 158个.
将块体模型的质心数据文件导出, 采用本文中提出的转换方法, 在MATLAB中通过SURPAC- RFPA3D接口程序STOR.m, 导入块体质心文件, 导出coor0.txt, elem0.txt, id0.txt和disp0.txt 4个文件, 将这4个文件导入RFPA3D中, 赋予各岩性单元相应的岩石力学参数, 如表 1所示.生成的RFPA3D计算模型如图 6所示.
在RFPA3D中, 对东边坡下的挂帮矿进行开挖(开挖区域见图 5).由图 7所示的计算结果分析可知, 空区距离边坡最近的一侧(面向露天坑一侧)会首先发生破坏, 破坏沿空区边缘向南北两边发展, 产生一道连续的南北走向的破坏带后, 空区上方边坡岩体失去下部的拱脚支撑后, 发生整体滑移, 向采空区塌陷, 地表形成崩落区.
为验证模型转换方法的可靠性, 建立与RFPA3D模型尺寸和单元划分形式相同的Flac3D模型, 采用表 1的岩石力学参数, 进行挂帮矿开挖计算, 得到如图 8所示的计算结果, 选取位移≥2 m的区域为崩落区域[11], 得到挂帮矿开采引起的崩落范围.
将RFPA3D和Flac3D模拟计算得到的边坡崩落范围与现场实际破坏情况进行对比分析, 如图 9所示, 模拟结果与现场实际观察的破坏情况基本吻合, 很好地解释了破坏发展的过程, 证明了SURPAC-RFPA3D接口的有效性.
1) SURPAC块体单元和RFPA3D计算单元的形状均为六面体, 以此为基础, 将SURPAC地质模型转化为RFPA3D计算模型, 利用SURPAC强大的建模能力来弥补RFPA3D在复杂模型建模中的不足.
2) 根据所建立的SURPAC和RFPA3D模型单元数据关系, 在MATLAB中使用m语言编制了转换程序STOR.m.
3) 使用SURPAC和RFPA3D模型转换程序STOR.m进行某大型露天矿东边坡挂帮矿开采的模拟计算, 结果表明:所编制的转换程序能够快速准确地进行模型转换, 并与Flac3D计算结果吻合, 转换程序有效、可靠.
[1] |
侯恩科, 吴立新, 李建民, 等.
三维地学模拟与数值模拟的耦合方法研究[J]. 煤炭学报, 2002, 27(4): 388–392.
( Hou En-ke, Wu Li-xin, Li Jian-min, et al. Study on the coupling of 3D geoscience modeling with numerical simulation[J]. Journal of China Coal Society, 2002, 27(4): 388–392. ) |
[2] |
Tang C A, Tang A B.
Applications of rock failure process analysis (RFPA) method[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2011, 3(4): 352–372.
DOI:10.3724/SP.J.1235.2011.00352 |
[3] |
Liang Z Z, Tang C A, Li H X, et al.
Numerical simulation of 3D failure process in heterogeneous rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(1): 323–328.
|
[4] |
Dai F, Wei M D, Xu N W, et al.
Numerical assessment of the progressive rock fracture mechanism of cracked chevron notched Brazilian disc specimens[J]. Rock Mechanics and Rock Engineering, 2015, 48(2): 463–479.
DOI:10.1007/s00603-014-0587-8 |
[5] |
Li L C, Tang C A, Zhu W C, et al.
Numerical analysis of slope stability based on the gravity increase method[J]. Computers and Geotechniques, 2009, 36(7): 1246–1258.
DOI:10.1016/j.compgeo.2009.06.004 |
[6] |
徐奴文, 唐春安, 周钟, 等.
基于三维数值模拟和微震监测的水工岩质边坡稳定性分析[J]. 岩石力学与工程学报, 2013, 32(7): 1371–1381.
( Xu Nu-wen, Tang Chun-an, Zhou Zhong, et al. Stability analysis of hydraulic rock slope based on three-dimensional numerical simulation and microseismic monitoring[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(7): 1371–1381. ) |
[7] |
Xu N W, Tang C A, Li L C, et al.
Microseismic monitoring and stability analysis of the left bank slope in Jinping first stage hydropower station in southwestern China[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(6): 950–963.
DOI:10.1016/j.ijrmms.2011.06.009 |
[8] |
Luo Z Q, Xie C Y, Zhou J M, et al.
Numerical analysis of stability for mined-out area in multi-field coupling[J]. Journal of Central South University, 2015, 22(2): 669–675.
DOI:10.1007/s11771-015-2569-8 |
[9] |
林杭, 曹平, 李江腾, 等.
基于SURPAC的FLAC3D三维模型自动构建[J]. 中国矿业大学学报, 2008, 37(3): 339–342.
( Lin Hang, Cao Ping, Li Jiang-teng, et al. Automatic generation of FLAC3D model based on SURPAC[J]. Journal of China University of Mining & Technology, 2008, 37(3): 339–342. ) |
[10] |
Li X B, Li D Y, Liu Z X, et al.
Determination of the minimum thickness of crown pillar for safe exploitation of a subsea gold mine based on numerical modeling[J]. International Journal of Rock Mechanics & Mining Sciences, 2013, 57(1): 42–56.
|
[11] |
Sainsbury D P, Sainsbury B L, Lorig L J.
Investigation of caving induced subsidence at the abandoned Grace Mine[J]. Mining Technology, 2013, 119(3): 151–161.
|