爆破荷载作用下岩石边坡动态响应的FLAC^3D模拟研究.pdf

第 2 2卷第4期 2 0 0 5年 1 2月 爆破 BLAS T I NG V0 1 ． 2 2 N o ． 4 D e c ． 2 0 o 5 文章编号 1 0 0 1 4 8 7 X 2 0 0 5 0 4 0 0 0 8 0 6 爆破荷载作用下岩石边坡动态 响应的 F L AC 3 D 模拟研究 陈占 军 ， 朱 传 云， 周 小 恒 武汉大学水资源与水电工程科学国家重点实验室 ， 武汉 4 3 0 0 7 2 摘要 已有的岩石边坡爆破动力响应分析多用有限单元法。作者尝试运用F L A C 的动力分析模块进行 岩石边坡爆破动力响应分析。建立了能够反映主要地质构造的三维岩石边坡数值模型。介绍了三维有限差 分程序 F L A C 动力分析的理论基础， 包括边界条件的设定、 动力荷载的输入、 阻尼的选取以及岩体本构模型 的选择。采用 F L A C 模拟了该岩石边坡在爆破荷载作用下的动态响应。分析了爆破结束后边坡体 内位移 场、 应力场、 速度场以及塑性 区的分布情况， 并与爆破前进行对比。将计算结果与实测数据进行了比较 ， 结果 表明F L A C 用于爆破荷载作用下岩石边坡动态响应的数值模拟是可行的。 关键词 爆破荷载； 岩石边坡； 动态响应； F L AC 模拟 中图分类号 T D 2 3 5 ． 1 文献标识码 A FLACa D S i mu l a t i o n f o r Dy n a mi c Re s p o n s e o f Ro c k S l o p e u nd e r Ex p l o s i o n C H E N Z h a n - j u n 。 Z HU C h u a n y u n 。 Z HO U X i a o kn g S t a t e K e y L a b o r a t o r y o f Wa t e r R e s o u r c e s a n d H y d r o p o w e r E n g i n e e r i n g S c i e n c e 。Wu h a n U n i v e r s i t y 。 Wu h a n 4 3 0 0 7 2 - C h i n a Abs t r a c t T h e d yn a c r e s p o n s e o f b l a s t i n g f o r r o c k s l o p e i s s i mu l a t e d b y d y n a mi c a l l y a n a l y z i n g t h e mo d ule o f F L AC p r o g r am ．A 3 - D r o c k s l o p e mode l t h a t c a n refl e c t i t s ma i n g e o l o g i c a l s t r u c t u re i s e s t a b l i s h e d ，Th e the o ret i c a l b a s i s o f d y n am i c a n aly s i s b y the F L AC p r o g r am i s i n t r o d u c e d．i n c l u d i n g s e t t i n g b o u n d a r y c o n d i t i o n s -i n p u tt i n g a y n a mi c l o a d i n g s 。c h o o s i n g d am p i n g and the p h y s i c al mode l o f t h e r o c k s l o p e ．Th e F L AC p r o g r a m i s u ti li z e d t o s i re ula t e t h e d y n am i c r e s p o n s e o f t h e r o c k s l o p e u n d e r e x p l o s i o n ．T h e d i s t ri b u t i o n s o f s t res s fi e l d，dis p l a c e me n t fi e l d， v e l o c i t y fi e l d and p l ast i c z o n e s i n the r o c k s l o pe a re an aly zed a f t e r b l ast i n g ．Th e s i mula t i o n r e s u l t s a r e c o mp a r e d w i th t h e s i t e mo n i t o ri n g res ult s ．R e s ult s s h o w s t h a t t h e F L AC D p r o g r am C an b e u s e d t o s i mula t e the d y n am i c res p o nse o f r o c k s l o p e u n d e r e x p l o s i o n e ff e c t i v e l y ． Ke y wo r d s e x p l o s i o n l o a din g；r ock s l o p e ；d y n am i c r e s p o n s e； F L AC 。s i mula t i o n 1 引 言 爆破荷载对岩石边坡稳定性 的影响是极其复杂 的岩石工程地质力学问题 ， 涉及到岩石动力学 、 工程 收稿 日期 2 0 o 50 9一l 1 ． 作者简介 陈占军 1 9 8 0 一 ， 男； 武汉 武汉大学水利水电学院研究 生． 地质学等相关学科 。其 影响主要有两个方面⋯ 一 是“ 弱化作用” ， 即由于爆破荷载的反复作用， 引起 岩体结构面的张开、 扩展和拉裂。 以及岩体的松动变 形， 造成岩体强度指标降低； 二是“ 附加荷载” 作用， 即爆破振动引起的惯性力导致边坡整体的下滑力加 大， 降低了边坡的安全系数。因此， 研究爆破荷载下 岩石边坡的动力响应特征， 并据此提出相应的控制 维普资讯 第 2 2卷第 4期 陈占军等爆破荷载作用下岩石边坡动态响应的 F L A C 如模拟研究 9 措施 ， 对于实际工程有重要 的指导意义。 F L A C F a s t L a g r a n g i a n A n a l y s i s o f C o n t i n u a ， 连续介质快速拉格朗日分析 是由 C u n d a l l 和美国 I T A S C A公司开发 出的有限差分数值计算程序 ， 主要 适用于地质和岩土工程的力学分析。由于其在材料 的弹塑性分析 、 大变形分析和施工过程模拟等领域 有独到的优点 ， 近年来 ， 广泛应 用于边坡稳 定性评 价、 支护设计及评价、 地下洞室、 施工设计 开挖、 填 筑等 、 河谷演化进程再现、 拱坝稳定分析、 隧道工 程 、 矿 山工程等 多个领域。 目前 ， 利用 F L A C 如进行 静力稳定分析及开挖过程模拟的文献[ 3 ～ 6 ] 较多， 而用其进行岩石边坡动力响应研究 的文献较少 ， 为 此 ， 作者尝试利用 F L A C 如对岩石边坡在爆破荷载下 的动态响应进行数值模拟。讨论了边坡动态响应分 析时边界条件 的设定、 阻尼 的选取以及爆 破荷载的 输入等关键 问题 ， 计算结果表 明 F L A C 。 用于爆破 荷载下岩石边坡动态响应的数值模拟是可行 的。 2 工程地质概 况 图 l为某矿 山采场岩石边坡的地质剖面图。该 边坡走向与爆破开采山体岩层 的走向一致。根据相 关地质资料， 靠近坡角的岩体为变粒、 浅粒构造变质 岩， 其中含有倾角约为 一 2 9 。 的层理面， 并含有锰土 层面和炭质片岩两处软弱结构面； 远离坡角的岩体 为花岗片麻岩， 强度较构造变质岩高， 其中含有倾角 约为 一 2 6 。 的层理面。完整岩石材料和节理的动态 力学参数分别 见表 1 ， 表 2 。 ①眼球状花 岗片麻岩；②炭锰质软弱层；③云英片岩； ④构造片麻岩；⑥含磷浅粒岩；⑦含磷变粒岩；P 为主含矿层 图 1 边坡地质剖面图 F i g 1 G e o l o g i c a l s e c t i o n o f t h e r o c k s l o p e 表 1 完整岩石的力学参数 Ta b l e 1 Pr o p e r t i e s o fi n t a c t r o c k ma t e r i a l s 表 2岩体结构面的力学参数 T a b l e 2 P r o p e r ti e s o f r o c k e cti O l l kq 3 计算模型的建立 根据该岩石边坡的地质剖面图， 截取一定范围 如图2 所示， 建立三维数值模型如图3 所示。根据 文献[ 7 ] 的研究结果， 网格剖分的尺寸受输入波动 的最短波长控制。设网格的最大尺寸为△ f ， 输入波 维普资讯 1 0 爆破 2 0 0 5年 1 2月 动的最短波长为 A， 则△ 2 必须小于 1 ／ 1 0 1 ／ 8 A。 考虑到此因素， 网格剖分时最小尺寸取为 2 m， 经验 算 ， 满足条件 △ 2 1 ／ 1 0 1 ／ 8 A 1 模型共剖分节点 5 8 2 0个 ， 单元 4 8 0 6个 。建模 时， 先采用前处理功能强大的有限元软件 A N S Y S 建 立三维模型 ， 然后通过读取 自编 的“ a n s y s t o fl a c ” 程 序， 将模型导人 F L A C 如。模型中采用不同的颜色表 示不同岩层， 从左至右依次为 花 岗片麻岩 、 层理 面 1 、 层理面2 、 锰 土层 面、 炭质 片岩和构造变质岩 ， 分 别命名为 g r o u p l g r o u p 6 。 层 理 图2 计算模型的边界条件 F i g 2 B o u n d a r y c o n d i t i o n s o f t h e c a l c u l a t i n g mo d e l 图3 边坡岩体动力分析模型 F i g 3 T h e d y n a mi c a n aly t i c al mo d e l o f t h e s l o p e 3 ． 1 边界条件的设定 F L A C 中求解动力问题的边界条件设置有远 置人工边界和粘滞边界2 种_8 ] 。远置人工边界条件 的办法在处理弹性波的能量逸散时， 沿用静力条件 下边界的设置， 这样的边界对人射波起着完全反射 的作用， 即既不传播， 也不吸收任何能量。而粘滞边 界其功能和 L S D Y N A中无反射边界相同， 即将模型 的边界取得足够远， 把模型的范围取得足够大， 使边 界对人射波的反射尽可能减少， 如此一来便可减小 边界条件对计算结果产生的影响。 在 F L A C 动力分析中， 粘滞边界是通过在边界 的法线方向和水平方向施加独立的粘滞力来实现 的 。其法向粘滞力 和剪切粘滞力 的计算式 为 一p c V 2 一p c V 3 式中， ， 分别为边界上的法向和切向速度分量； p 为岩体密度； c 。 ， C s 分别为 P波和 s波的波速。当波 动人射角小于 6 0 。 时， 粘滞边界对于体波的吸收是 有效的。 3 ． 2 动力荷载的输入 爆破动力荷载的确定是进行爆破荷载下岩石边 坡动力响应研究的基础， 也是一大难点。目 前， 在实 际的工程应用 中， 炮孔压力变化历程的选取有半经 验、 半理论的炮孑 L 压力指数衰减公式， 如文献[ 1 0 ] 将作用在结构上的爆破荷载简化为具有线性上升段 和下降段的三角形荷载。也有的将炮孑 L 压力变化历 程简化为级数形式。 然而， 简化的三角形或级数形式爆破荷载与实 际情况有较大的误差。为避免此误差， 数值模拟中， 爆破动力荷载直接采用现场实测的爆破振动速度时 程曲线 输人， 如图4 所示。 0- 3 0． 2 0 1 ‘ 0 目0 一0_ 一 0． 2 -0． 3 0．4 水平径 向速度曲线 ⋯⋯一 垂直向速度曲线 ‘ ．A 。 。 ～一 3 0 { D ‰ ． 0 ． 5 0 ． _ f ’ 时间， s _ - 图4 实测爆破振动速度时程曲线 F i g4 T i me h i s t o ry of v e l o c i t y o b t a i n e d b y s i t e mo n i t o r i n g 值得注意的是， 对于粘性边界， 在输人动荷载 时， 须将速度时程曲线转化成应力时程曲线， 可通过 下式进行转化 o r 2 p c v 4 or 2 p c 5 式中， or 、 or 分别为正应力和剪应力；p为岩体密度； c 、 c 分别为 P 波和s波的波速； 、 分别为质点在 法向和切向的振动速度， 即分别为现场实测的水平 和竖直向振动速度。 维普资讯 第2 2卷第 4期 陈占军等爆破荷载作用下岩石边坡动态响应的 F L A C 如模拟研究 l 1 3 ． 3 阻尼的选取 进行动力响应分析时， F L A C 。 。 提供了2种可供 选择的阻尼型式， 即瑞利 R a y l e i g h 阻尼和局部 1 o - c a 阻尼。其中， 瑞利阻尼应用较多， 其通过下式确 定 Ca M ／ 3 K 6 式中， O ／ 、 J8 分别为质量阻尼系数和刚度阻尼系数。 设置瑞利阻尼时， 须确定一个中心频率 c e n t e r f r e q u e n c y j 。 ， 其值由式 7 给出 厂 m i i ／ 2 a t 7 本例中， 输入动荷载的主频为4 6 ． 7 H z ， 其 i 为 7 ． 8 ， 在加载时， 用命令“ s e t d y n d a m p r a y l e i g h 1 7 ． 8 ” 来实现。 3 ． 4 岩体本构模型的选取 F L A C 。 中提供了由空单元模型、 弹性模型和塑 性模型组成的共 1 0种基本的本构模型。模拟时选 取 M o h r C o u lo m b 莫尔一 库伦 弹塑性模型。该模型 通常用 F 描述土体和岩石的剪切破坏。模型的破坏 包络线和 M o h r C o u l o m b 强度准则 剪切屈服函数 以及拉破坏准则 拉屈服函数 相对应。图5所示 即为 M o h r C o u l o m b 强度准则示意图。 图5 M o h r ． C o u l o m b 强度准则示意图 F i g 5 I n t e n s i t y p r i n c i p l e o f Mo h r - C o ulo mb mo d e l 4 模拟结果及分析 应用 F L A C 如 进行动力响应分析前， 须进行静力 条件下的边坡稳定分析。为此， 该数值模拟分二步 进行 第一步， 施加自重， 分析该岩石边坡在初始应 力条件下的应力和变形 ； 然后， 将第一步计算的位移 初始化为零， 再输入动力荷载， 分析其在爆破荷载作 用下的动态响应特征。以下为部分计算结果。 4 ． 1 初始状态 图6 为边坡岩体在初始状态下的位移矢量场分 布示意图。从图6中可以看出， 初始状态下边坡体 的位移很小， 最大处仅有 8 ． 5 n l／ n ， 其在初始地应力 作用下是稳定的。相比之下， 边坡体在坡脚和坡顶 处的位移略大于其他地方， 因此这 2个区域可看作 边坡体的危险区域。为研究边坡体在爆破荷载作用 下的动态响应， 下面的结果分析特选取这 2个区域 的质点作为研究对象。 图6 初始状态下边坡的位移矢量图 F i g 6 T h e d i s p l a c e me n t v e c t o r o f r o c k s l o p e i n i n i t i a l s t a t e 4 ． 2 动态响应 爆破荷载作用下岩石边坡动态响应的分析包括 应力、 位移、 速度矢量的分布， 以及塑性区的分布特 征。下面分别就这几项对模拟结果进行动力分析。 4 ． 2 ． 1 应力场分布 在爆破荷载作用下， 边坡体内竖向应力分布如 图7 所示。从图中可以看出 在爆破荷载作用下， 边 坡体内各点应力均较小， 最大处 软弱层理面 仅为 0 ． 3 0 M P a 。这说明， 爆破后， 边坡整体仍处于稳定 状态， 这与爆破后边坡稳定状态的宏观观测结果吻 合 。 图7 爆后坡体竖向应力分布示意图 F i g 7 Ve r t i c a l s t r e s s o f r o c k s l o p e a f t e r b l a s t i n g 4 ． 2 ． 2 位移场分布 计算终态时， 边坡体内的竖向位移分布如图8 所示。从图8 中可以看出 在爆破荷载作用下， 边坡 体内坡脚处和坡顶处的竖向位移较其他部位明显大 很多。坡脚处最大竖向位移为 2 ． 5 c m， 坡顶处最大 竖向位移也达到了2 ． 2 e m 。且随着时间的推移， 坡 体的位移呈现处稳定趋势， 即产生了永久位移。 a a a a a a a P a P 胁 胁 胁胁 胁 胁 胁 胁 如 加 2 m ∞ 维普资讯 1 2 爆破 2 0 0 5年 1 2月 0 ．3 c m O ．6 c m 0 ． 9 c m 1 ．2 c m l 5 c m 1 ．8 c m 2 ． 1 c m 2 ．2c m 2 ．4c m 2 ．5 c m 3 ．0G i n 图8 爆后坡体竖向位移分布示意图 F i g 8 V e r t i c a l d i s p l a c e me n t o f r o c k s l o p e a f t e r b l a s t i n g 4 ． 2 ． 3 塑性 区分布 计算终态时， 边坡体内塑性区分布如图9 所示。 从图9中可以看出 边坡体内塑性区主要集中在岩 体软弱带内， 且多为受剪应力屈服 图中的红色单 元 和受拉应力屈服 图中的绿色单元 ， 这说明在 爆破荷载作用下， 冲击应力波对边坡体的破坏作用 是从降低其抗拉强度和抗剪强度开始的， 这与岩体 的力学试验 所得结论吻合。 图 9 爆后边坡体内塑性 区分布示意 图 F i g 9 P l a s t i c z o n e s o f r o c k s l o p e a f t e r b l a s t i n g 4 ． 2 ． 4 速度场分布 计算终态时， 边坡体内Y向速度场分布如图 l O 所示。从图 1 0中可以看出 在爆破荷载作用下 ， 边 坡体内坡脚和坡顶处速度较其他地方明显大很多， 在坡脚处其最大值达0 ． 4 1 m ／ s 。计算中， 选取坡脚 处一点， 跟踪记录其振动速度峰值 ， 其结果与实测 值 】 的比较见表 3 。从表 3中可以看出， 计算值与 实测值间的误差在 2 0 ％以内， 属较小范畴， 因此可 以认为此次数值模拟计算是有效的。图 1 1 为该点 垂直向振动速度计算值与实测值的对比曲线图。 表3 坡脚处质点振动速度计算结果与实测结果比较 T a b l e 3 Co mp a r i s o n o f p e a k p a r t i c l e v e l o c i t y b e t we e n c a l c u l a t i n g r e s u l t s a n d s i t e m o n i t o r i n g r e s u l t s 图 1 0 爆后边坡体内竖向速度分布示意图 F i g 1 0 Ve rt i c al v e l o c i t y o f r o c k s l o p e aft e r b l a s t i n g 图 l 1 坡脚点垂直向计算与实测速度时程曲线 F i g 1 1 Vv e r t i c al v e l o c i t y h i s t o r y o f s l o p e t o e o b t a i n e d b y c o mp u t a t i o n a n d s i t e s mo n i t o r i n g 5 结 论 作者尝试运用 F L A C 如进行岩石边坡在爆破作 用下动态响应的数值模拟。详细讨论了数值模拟时 边界条件的设定 、 阻尼的选取 、 爆破动力荷载的输入 等关键问题。经计算分析 ， 爆后边坡体内位移、 应 力、 速度及塑性区分布等动态响应特征， 并与实测结 果进行比较， 误差较小， 结果说明 F L A C 如 用于进行 爆破作用下边坡动态响应的数值模拟是可行的。 参考文献 [ 1 ] 朱传云， 卢文 波． 岩质边坡爆破振动安全判据综述 [ J ] ． 爆破， 1 9 9 7 ， 1 4 4 1 31 7 ． Z h u C Y，L u W B． Ge n e r a l S t u d y o n C o n t r o l l i n g S t and a r d o f B l a s t i n g S a f e t y f o r R r o c k S l o p e [ J ] ． B l a s t i n g ，1 9 9 7 ， 1 4 4 1 3～1 7 ． [ 2 ] 黄润秋， 许强． 显式拉格朗 日差分分析在岩石边坡 工程中的应用[ J ] ． 岩石力学与工程学报， 1 9 9 5 ， 1 4 3 4 6 3 5 4． H u a n g R Q， X u Q． A p p l i c a t i o n o f E x p h c a t e d F ast L a - g r a n g i a n An aly s i s o f C o n t i n u a i n Ro c k S l o p e E n g i n e e r i n g [ J ] ． C h i n e s e J o u rnal o f R o c k M e c h a n i c s a n d E n g i n e e r - i n g ， 1 9 9 5 ， 1 4 3 4 6 3 5 4 ． 维普资讯 第2 2 卷第 4 期 陈占 军等爆破荷载作用下岩石边坡动态响应的 F L A C 如 模拟研究 1 3 [ 3 ] 丁秀美， 黄润秋． F L A C ． 3 D前处理程序开发及其工程 应用[ J ] ． 地质灾害与环境保护， 2 0 0 4 ， 2 3 3 3 5 ． D i n g X M， H u a n g R Q． D e v e l o p m e n t o f P r e - p r o c e s s i n g S o f twa r e f o r F L A C a n d I t s Ap p l i c a t i o n t o E n g i n e e ri n g [ J ] ． J o u r n a l o f G e o l o g i c al H a z a r d s a n d E n v i r o n m e n t P r e s ． e r v a t i o n ， 2 0 0 4， 1 5 2 3 3 3 5 ． [ 4 ] 杨典森， 陈卫忠． 龙滩地下洞 室群 围岩稳定性分析 [ J ] ． 岩土力学 ， 2 0 0 4， 2 5 3 5～ 8 ． Ya n g D S， C h e n W Z． S t a b i l i t y A n aly s i s o f S u r r o u n d i n g R ock Ma s s o f L o n g t an U n d e r g r o u n d C a v e rns [ J ] ． R o c k a n d S o i l Me c h a n i c s ， 2 0 0 4 ， 2 5 3 5 8 ． [ 5 ] 寇晓东， 周维垣， 杨若琼． 三维快速拉格朗日 法及其在 拱坝稳定分析中的应用 [ J ] ． 水利水电技术 ， 2 O 0 0 ， 3 1 ．7 4 7 ． K o u X D， Z h o u W Y， Y a n g R Q ． T h r e e d i m e n s i o n al F a s t L a gra n g i an Me t h o d a n d I t s Ap p l i c a t i o n o n S t a b i l i t y An aly - s i s If t h e A r c h D a m[ J ] ． Wa t e r R e s o u r c e s a n d H y d rop o w ． e r E n gin e e ri n g ， 2 0 0 0 ， 3 1 7 4～ 7 ． [ 6 ] 谢和平， 周宏伟 ， 王金安， 等． F L A C在煤矿开采沉陷 预测中的应用及对比分析[ J ] ．岩石力学与工程学报， 1 9 9 9 ， 1 8 4 3 9 7～ 4 0 1 ． X i e H P， Z h o u H W ， Wa n g J A， e t a1． C o n t r ast An a l y s i s and A p p l i c a t i o n o f F L A C i n I mme r s i o n P r e d i c t i o n f o r C o al Mi n e D e v e l o p m e n t [ J ] ． C h i n e s e J o u rnal o f R o c k Me c h a n ． i c s and E n gi n e e ri n g ， 1 9 9 9， 1 8 4 3 9 7 4 0 1 ． [ 7 ] K u t fl e m e y e r R L， L y s m e r J ． F i n i t e E l e m e n t M e t h o d A c c u ． r a c y for Wa v e P r o p a g a t i o n P r o b l e m s[ J ] ． J S o i l Me c h＆ F o u n d a t i o n s D i v ， A S C E， 1 9 7 3 ， 9 9 5 4 1 7～ 4 2 1 ． [ 8 ] 刘春玲， 祁生文． 利用 F L A C 如分析某边坡地震稳定性 [ J ] ．岩 石 力 学 与 工 程 学 报， 2 0 0 4 ，2 31 6 2 7 3 0～2 7 3 3 ． L i u C L ， Q i S W． S t a b i l i t y A n al y s i s o f S l o p e u n d e r E a r t h - q u a k e w i t h F L A C 如[ J ] ． C h i n e s e J o u r n al o f R o c k M e c h a n ． i c s F n gi n e e r i n g ， 2 0 0 4 ， 2 3 1 6 2 7 3 0～2 7 3 3 ． [ 9 ] I t a s c a C o n s u l t i n g G r o u p I n c ．F L A C - 3 D F ast L a gr ang i a n A n al y s i s o f C o n t i n u a i n 3 D i m e n s i o n s ， V e rs i o n 2 ． 1 0 ． U s e rs Manu al V o l u me V[ R ] ．U S A I t a s c a C o n s u l t i n g Gr o u p I n c ， 1 9 9 7 ． [ 1 O ] 刘国华， 王振宇． 爆破荷载作用下隧道的动态响应和 抗爆分析[ J ] ． 浙江大学学报 工学版 ， 2 0 0 4 ， 3 8 2 1 5～1 8． L i u G H， Wan g Z Y． Dy n a mi c Re s p o n s e an d B l a s t i n g ． R e s i s t a n c e A n al y s i s o f a T u n n e l S u b j e c t e d t o B l a s t i n g L o a d i n g s [ J ] ． J o u rnal o f Z h e j i ang U n i v e rs i t y ，2 0 0 4 ， 3 8 2 1 5～1 8 ． [ 1 1 ] 刘亚群， 李海波． 爆破荷载作用下黄麦岭磷矿岩质边 坡动态响应的 U D E C模拟研究[ J ] ． 岩石力学与工程 学报， 2 0 0 4， 2 3 1 1 4 2 4 5 ． L i u Y Q ． L i H B ． U D E C S i m u l a t i o n o n D y n a mi c R e ． s p o n s e o f Ro c k S l o p e o f Hu an g ma i h n g P h o s p h o r i t e Mi n e U n d e r E x p l o s i o n [ J ] ． C h i n e s e J o u r n al o f R o c k M e c h ani c s and E n g i n e e ri n g ， 2 0 04 ， 2 3 1 1 4 2～ 4 5 ． [ 1 2 ] 徐国元， 古德生． 爆破破岩机理的试验研究 [ J ] ． 中南 工业大学学报， 1 9 9 7 ， 1 2 6 1 ～ 6 3 ． Xu G Y， Gu D S ． E x p e ri me n t al I n v e s t i g a t i o n o n Me c h a ． n i s m o f R o c k F r a g m e n t a t i o n b y B l ast i n g [ J ] ． J o u rnal o f C e n t r a l S o u t h U n i v e r s i t y ，1 9 9 7 ， 1 2 6 1 ～ 6 3 ． [ 1 3 ] 中国科学院武汉岩土力学研究所岩土工程检测中 心．黄麦岭磷矿采场爆破振动监测、 声波检测及爆破 设计优化报告 [ R] ． 武汉 中国科学院武汉岩土力学 研究所， 2 0 0 0 ． I n s t i t u t e o f Ro c k a n d S o i l Me c h a n i c s ， T h e C h i n e s e A c a d ． e m y o f S c i e n c e ， T h e E n gi n e e ri n g T e s t i n g C e n t e r ． R e p o a o n Op t i mi z i n g B l a s t i n g De s i g n i n g a n d t h e Mo n i t o ri n g o f B l ast i n g Vi b r a t i o n o f Hu a n gm a i l i n g P h o s p h o r i t e Mi n e [ R] ． Wu h an I n s t i t u t e o f R o c k a n d S o i l M e c h a n i c s ． The C h i n e s e A c a d e my o f S c i e n c e ， 2 0 0 0． 上接第3页 4 介质参与后对爆生气体做功能力和冻土爆 破中爆破参数的影响及如何确定这些参数也有待研 究。 参考文献 [ 1 ] 周幼吾， 郭东信， 邱国庆． 中国冻土[ M] ． 北京 科学出 版社， 2 0 0 0 ． [ 2 ] 王省身． 矿井灾害防治理论与技术[ M] ． 徐州 中国矿 业大学出版社， 2 0 0 0 ． [ 3 ] 高尔新 ， 杨仁树． 爆破工程[ M] ． 徐州 中国矿业大学 出版社 ， 1 9 9 9 ． [ 4 ] 杨振声． 工程爆破文集[ c ] ． 深圳 海天出版社， 1 9 9 7 ． [ 5 ] 韩修栋 ， 白宝库． 刚性被筒炸药的研究[ J ] ． 爆破器材， 1 9 9 1 ， 1 1 7～1 9 ． 维普资讯