。510· 北京科技大学学报 2006年第6期 13-2 143T 141 17 28 17-2* F断只600 W987 T075 -2 地下开挖 10 50 400 E100 050 图1三维边坡计算简图（单位：m) Fig I Schematic plot of the 3D model of the slope Unit m) 的质黄灰岩 模型（如图4）.考虑到岩石的尺寸效应，对岩石实 验参数进行了适当的折减，计算采用的岩石力学 度秋岩 参数和断层力学参数分别见表1和表2. 量色泥岩 同填材料 简母贝 断层 回采工作面 蓝从岩 图2地质剖面图 Fig.2 Gedlogic sectional profile 2模拟计算 图3地下开采前边坡图 2.1模型建立 Fig.3 Slope before underground mining 计算采用三维有限差分程序FLAC3D.计算 0-0 模型x方向长1400m,y方向宽1050m,z方向 高600m.共划分为156118个单元，168188个节 莫尔一库仑准 点.模型侧面限制水平移动，底面限制垂直移动， 则与应变软化 上部为自由面.图3为地下开挖前三维数值 模型. 依据现场取样和岩石力学实验结果，岩石在 不同围压条件下具有明显的弹塑性变形特征.当 岩石发生破坏后，表现为应变软化特征，即随着变 e 形的继续发展，残余强度逐渐降低，最终下降到峰 图4岩石的力学特性 值强度的20%左右.因此，采用非线性应变软化 Fig.4 Mechanical character of rocks图 1 三维边坡计算简图( 单位:m) Fig.1 Schematic plot of the 3D model of the slope ( Unit:m) 图 2 地质剖面图 Fig.2 Geologic sectional profile 2 模拟计算 2.1 模型建立 计算采用三维有限差分程序 FLAC 3D .计算 模型 x 方向长 1 400 m, y 方向宽 1 050 m, z 方向 高600 m .共划分为 156 118 个单元, 168188 个节 点.模型侧面限制水平移动, 底面限制垂直移动, 上部为自由面 .图 3 为地下开挖前三维数值 模型 . 依据现场取样和岩石力学实验结果, 岩石在 不同围压条件下具有明显的弹塑性变形特征.当 岩石发生破坏后, 表现为应变软化特征, 即随着变 形的继续发展, 残余强度逐渐降低, 最终下降到峰 值强度的 20 %左右 .因此, 采用非线性应变软化 模型( 如图 4) .考虑到岩石的尺寸效应, 对岩石实 验参数进行了适当的折减, 计算采用的岩石力学 参数和断层力学参数分别见表 1 和表 2 . 图 3 地下开采前边坡图 Fig.3 Slope before underground mining 图 4 岩石的力学特性 Fig.4 Mechanical character of rocks · 510 · 北 京 科 技 大 学 学 报 2006 年第 6 期