李斓堃 朱万成 代 风 刘溪鸽 邓文学

（东北大学资源与土木工程学院，辽宁 沈阳 110819）

近年来，数值模拟被广泛应用于采矿工程中的岩体稳定性分析。大型岩土工程计算区域常达千米级，并且地质构造复杂、岩石类型多样以及涉及到多场耦合等问题。Park D等［1］建立了长达1 000 m的隧道模型，并得到隧道沿轴线方向的变形响应特征；尚振华等［2］采用FLAC3D软件对尺寸为1 000 m×1 600 m×660 m的某钨矿模型进行模拟，将采空区稳定性问题量化；刘春等［3］建立了尺寸为2 800 m×1 500 m×1 600 m的茂县新磨村滑坡模型，模拟了新磨村滑坡启动、高速下滑和堆积的全过程；张成良等［4］采用三维有限元模拟了3种开采方案对地表公路的影响，得到当开采高度超过1 300 m水平时，公路上部围岩将处于不稳定状态；漆祖芳等［5］建立630 m×1 850 m×900 m的有限元模型，模拟了大岗山水电站坝肩边坡开挖支护过程；马克等［6］用RFPA3D对大岗山水电站坝址区右岸实际微震监测范围（690 m×400 m×1 600 m）等比例的三维模型进行模拟，对抗剪洞加固前、后的边坡稳定性进行分析。

考虑到计算成本及技术方面的问题，前人对计算模型进行了大量简化，且受到计算模型网格数量受限，致使模型网格尺寸较大，导致了数值模型难以反映实际工程的现场真实情况［7］。因此，对于大规模、复杂岩石工程的数值模拟，大多数情况下会将关键位置网格和非关键位置网格差异化对待，即关键位置网格剖分需要细化，而非关键位置可采用过渡单元，将单元尺寸从核心研究区域向边界区域逐渐粗化［8-9］。……