爆破震动对采空区稳定性影响的FLAC^3D分析.pdf

第 2 4卷第 1 6期 2 0 0 5 年 8月 岩石 力 学与 工程 学报 C h i n e s e J o u r n a l o f R o c k Me c h a n i c s a n d E n g i n e e r i n g Vb 1 ． 2 4 N0． 1 6 Au g ． ， 2 0 0 5 爆破震动对采空区稳定性影响的 F L A C 3 D 分析 闫长斌 ，徐 国元，李夕兵 中南大学 资源 与安全工程学院，湖南长沙4 1 0 0 8 3 摘要针对厂坝铅锌矿受到乱采乱挖破坏的情况，运用 F I A c ∞ 对爆破震动作用下采空区的稳定性进行了数值分 析。根据现场 实测数据 ，选取 了 2条不 同的实测速度时程 曲线作为爆破震动输入 ，分别进行了静力和动力计算 。 分析 结果表明 1 采 空区开挖后 ，最 大位移量 为 4 ． 7 5 m m 2 施加爆破动 载荷后，最大位移量增至 5 ． 4 7 mm， 增幅与爆破震动速度成正比； 3 围岩位移分布的变化受震动波形的影响。 爆破震动作用下围岩应力场发生二次重 分布 ，塑性区分布面积 明显扩大 ，采 空区右上 方露天边 坡脚亦受到影响 。分析结果与实际观测结果相吻 合，并给 出了避免爆破震动诱发采空 失稳塌陷的几条建议 。 关键词爆破工程；采矿工程；爆破震动；采空区；稳定性F L AC ；数值模拟 中圈分类号T D 2 3 5 ；T D 8 5 文献标识码A 文章编号1 0 0 0 6 9 1 5 2 0 0 5 1 6 2 8 9 4 0 6 S TABI LI TY ANALYS I S OF M I NED． OUT AREAS I NFLUENCED BY BLAS TI NG VI BRATI ON W I TH FLAC3 D YAN Cha n g- b i n， XU Gu o - y u a n， LI Xi - bi ng S c h o o l o f R e s o u r c e s a n d S a f e t y E n g i n e e r i n g ，C e n t r a l S o u t h U n i v e r s i t y ，C h a n g s h a 4 1 0 0 8 3 ，C h i n a Abs t r a c t Fo r t he r e a s o n t h a t Cha ng b a l e a d - z i n c mi n e i s s e v e r e l y d e s t r o ye d b y di s or de r e d mi n i n g， t h e n ume ric a l a na l ys i s o f mi n e d ． o u t a r e a s un d e r bl a s t i n g v i br a t i o n i S c a r r i e d o ut wi t h FLAC。 ． Two di fie r e nt me a s u r e d v e l oc i t y t i me ． d e p e n d e n t c u r ve s a r e s e l e c t e d a s b l a s t i n g v i b r a t i o n i n p ut a c c o r d i n g t o t h e i n- s i t u t e s t da t a； a nd s t a t i c c a l c ul a t i o n a n d d y n a mi c c a l c u l a t i o n a r e c a r r i e d o ut ， r e s p e c t i v e l y ．The a na l y t i c a l r e s ul t s s h ow t ha t t h e ma x i mum d i s p l a c e m i ． ． ， ．e n t o f s ur ro u n d i n g r oc k S 4 7 5 m m a f t e r mi n e d o ut a r e a s are e x c a v a t e d wh i l e t h e ma x i mum d i s pl a c e m e n t o f s u rro u n di n g r o c k i n c r e a s e s t o 5 ． 47 mm a f t e r t h e b l a s t i n g d y na mi c l o a d i n g i s a p p l i e d， a n d a mp l i t u d e of i n c r e a s i n g d i s p l a c e me n t i s pr op o r t i on a l t o t he v e l o c i t y o f bl a s t i n g v i br a t i on． Th e s ha p e o f v i b r a t i o n wa ve i n flue nc e s t h e v a ria t i on o f di s p l a c e me n t d i s t rib ut i o n f o r s urrou n di n g r oc k．W h e n t h e b l a s t i n g v i b r a t i on d yn a mi c l o a di n g i S a pp l i e d， t h e s e c on d r e d i s t ribu t i o n o f s t r e s s fie l d o f s u rro un d i n g r o c k oc c u r s a n d t h e di s t rib ut i on are a o f p l a s t i c z o ne i S e xp a n d e d o bv i o us l y．Th e fo o t o f t h e fig h t op e n s l o pe a b ov e t h e mi n e d - ou t are a i s a l s o i nflu e n c e d a s we l 1 ．Th e a na l y t i c a l r e s u l t s a r e a c c o r da nt wi t h t he ob s e rve d r e s u l t s ．S e v e r a 1 s ug g e s t i o n s t o a v oi d i ns t a bi l i t y a n d s u bs i d e o f mi n e d． ou t a r e a s i n d u c e d b y b l a s t i n g v i br a t i on are br o u gh t for ward ． Ke y wo r ds bl a s t i ng e n gi ne e rin g；mi n i ng e n g i n e e ring；b l a s t i n g v i b r a t i o n；mi n e d- o ut a r e a s ；s t a bi l i t y；F LAC ； n u me r i c a 1 s i mul a t i on 1 引 言 爆破震动是爆破 公害之一 ，不但对地面建 构 筑物 的安全产生影响，而且对采空区等地下洞室的 围岩产生损伤和扰动，甚至引起 失稳破坏 卜 。对 爆破震动危害的研究，主要是进行现场测试 ，根据 实测数据 ，结合爆破安全规程，确定最大震速和安 收藕日期l 2 0 0 5 0 42 1 修回日期2 0 0 5 0 6 1 4 基金项 目l国家自然科学基金重人项 目 5 049 0 2 7 2 ，5 0 4 9 0 2 7 4 国家 “ 十五”科技攻关项 目 2 0 0 3 B A 6 1 2 A一1 0 2 ；中南大学博 士学位论文创新工程 项 目 04o 1 0 9 作t■介t闰长斌 1 9 8 O一 ， 男， 2 0 0 1 年毕业于中南1 业大学建筑工程专业， 现为博士研究生， 主要从事岩土工程与地下工程方面的研究工作。 E m ml y a n c h a n g b i n 一2 0 01 1 6 3 ．c o r n。 维普资讯 第 2 4卷第 1 6期 闰长斌等． 爆破震动对采 空区稳定性影响的 F L AC 。 分析 2 8 9 5 全距离 。该方法能直观地反映爆破震动危害，但无 法对爆破震动引起的围岩动应力场、位移和塑性区 的变化以及动力失稳机制进行研究。而数值模拟方 法可以系统地研究爆破作用下围岩的动应力场分布 规律及动静应力场的叠加作用机理，因此成为近年 来研究爆破震动作用的有效方法L 5 】 。 鉴于 F L A C 在岩土力学分析中的优势和求解动力问题 的特点， 使得其可很好地解决非线性动力分析问题。F L A C Ⅲ 已成功应用于岩土开挖、边坡稳定分析及地震动力 响应分析等许多领域 。 厂坝铅锌矿区位于甘肃省成县的黄渚镇，由厂 坝、小厂坝和李家沟组成，是我国特大型铅锌矿床 之一 。2 0世纪 8 0年代后期以来，当地乱采乱挖 留 下的数百万平方米未处理群采空区，给矿山安全生 产构成了严重威胁。群采空区诱发的地质灾害时有 发生， 致使矿山曾连续发生重大事故I 叫引 。 特别是 近年来，爆破地震波作用下大范围采空区发生动力 失稳 的趋势 日益 明显。为保障露天转地下过渡层的 安全生产，进行 了爆破地震效应现场测试研究，获 得 了较准确的基础数据。为进一步研究复杂采空区 在爆破震动作用下产生动力失稳的过程和机理 ，本 文在现场测试的基础上，同时考虑开挖过程和爆破 震动对采空区稳定性的影响，利用 F L AC 进行了 三维数值模拟研究。 2 F L A C 3 D 动力计算步骤 F L A C D是二维的有限差分程序 F L AC D的扩 展，能够进行岩石 、土质和其他材料在达到屈服极 限后经历塑性变形的三维空间行为分析，为岩土工程 领域求解三维 问题提供了一种理想的分析工 。 运用 F L AC 叩进行动力计算，必须首先进行静 力分析 。在完成静力分析的基础上 ，才能施加动力 载荷进行动力分析。F L AC 叩 动力计算大致可以分为 以下几个步骤 1 确定计算区域，并进行网格划分 2 选择动力计算模式，定义本构模 型和材料 的物理力学参数； 3 定义计算所需的边界条件和初始条件； 4 进行计算 ，获得初始平衡状态，即开挖前 的原岩应力状态 ； 5 进行工程 采空区 开挖计算分析， 得到开挖 后的静力计算结果 6 检查静力计算结果 ，认为满意后 ，设置动 力计算边界条件和所需的阻尼； 7 施加动力载荷，进行动力计算分析，得到 动力计算结果。一般的计算流程如图 1 所示。 图 1 F L A C 。 动力分析流程图 F i g ． 1 Ho w c h a r t o f F LAC’ 。d y n a mi c a n a l y s i s 3 计算分析 3 ． 1 计算模型 3 ． 1 ． 1模型范围 根据前期采空区探测成果 ，2 矿体存在范 围较 大、互相贯通的采空区。因此，本次计算区域选择 2 矿体及周边围岩。计算区域尺寸为长 1 2 0 m，宽 1 8 0 m，高 1 2 5 I n 。采空区位于该区域的中部，即 1 2 0 2 1 2 5 0分段。贯通 的采空区尺寸为高 2 4 I n ， 横断面为下宽上窄，底宽 3 6 I n ，顶宽 1 1 I n ，沿矿体 延伸 3 2 I n 。上边界取至地表 ，下边界距离采空区底 板 7 0 I I 1 ， 左右边界踞离采空区分别为 9 0 I n和 6 0 I n ， 沿矿体方向 y方 向 取 8 8 I n 。根据圣维南原理，模 型大小满足计算精度要求 。 3 ． 1 ． 2模 型介绍 采空区开挖后，模型的 F L AC 3 D 剖分网格如图 2 所示。该模型共有 7 0 2 0 0个单元和 7 5 7 3 3个节点 维普资讯 岩石 力学 与j 程 学报 2 0 0 5缸 图 2 开挖后的 F L AC 。 模型 Fi g ． 2 F LAC 。mo d e l a f t e r e x c a v a t i o n 按岩性参数的不同，可将模型分为 3个区域 ，即采 空区、矿体和围岩 。模型中的采空区的尺寸 、形态 和位置以前期探测结果为准 。根据现场调查 ，模型 范围内岩性 比较简单，可分为围岩和矿石 2大类 ， 矿体上下盘围岩均为黑云母片岩，矿石为铅锌矿石。 为便于计算，模型简化为准三维情况，即不考虑断 面形状和尺寸在长度方 向 1 ， 方 向 上的变化 。 3 ． 1 ． 3初始应力场和边界条件 初始地应力场 由于采空区埋深较浅，且矿区 构造应力不明显 ，因此未进行现场地应力测试 ，初 始地应力场仅按 自重应力场考虑。根据弹性力学原 理可知，竖向应力和水平应力分别为 在 F L AC D计算中，岩体变形参数采用的是体 积模量 和剪切模量 G 。因此， 必须将弹性模量 或 变形模量 D和泊松 比 转化成体积模量 和剪 切模量 G ，转化公式如下 E 2 3 1 2 G 3 2 1 3 ．3 计算方案 为研究采空区开挖形成的围岩应力重分布和产 生的位移变化，找出爆破震动作用下采空区围岩应 力发生二次重分布的规律和位移变化 ，确定爆破震 动对采空区稳定的危害程度，运用 F L AC ∞对开挖 和爆破震动作用下采空区的稳定性进行了较为系统 的模拟分析 。分析结果可为预防采空区失稳塌陷和 减震、加固、治理提供依据 。根据震动载荷形式， 共有以下 3种计算方案 1 白重应力场中，采空 区开挖过程计算 ； 2 采 空区在中深孔爆破载荷作 用下的动力计算 ； 3 采空区在斜坡道掘进爆破载 荷作用下的动力计算。本次计算采用 Mo h r - C o u l o mb 准则，即弹塑性模型，并选取大应变计算模式。 4 爆破震动输入 v } 1 4 l 爆 破 震 动 速 度 时 程 t7 h k t r J 为获得第～手资料，对厂坝铅锌矿厂坝矿区进 式中 为埋置深度 ；k 为侧 压力系数 ，且 k ／ 1 一 ， 为泊松 比。 边界条件采用位移边界条件，即模型的左右 方向 边界、前后 1 ， 方 向 边界和底边界均施加位 移约束条件，上边界为 自由边界。 3 ． 2 岩体物理力学参数 本次计算模型中岩性较简单 ，仅分为矿石和围 岩 2种 。计算所需的岩体物理力学参数是根据室内 试验确定的岩石物理力学参数I l 6 J 折减得来的，具体 参数结果见表 1 。 裹 1 厂坝铅锌矿岩体物理力学参数 T a b l e 1 Ph y s i c o - me c h a n i c a l p a r a me t e r s o f t h e r o c k s i n Ch a n g b a l e a d - z i n c mi n e 行 了爆破震动现场测试研究 。测试分两个部分，一 是中深孔爆破强制崩矿地震效应测试 ，二是斜坡道 掘进爆破地震效应测试 。有关研究和实践证 明，地 震波的震动强度可用质点振动位移、速度、加速度 3个物理量来表达 。而震动速度既能反映震动 能量 大小，又与爆破产生的动应 力成正比，且测试方便 可靠，因此观测物理量 以质点振动速度为主。 测试仪器为成都中科动态仪器有限公司研制 的 I D T S 3 8 5 0二通道和三通道爆破振动记录仪，传感器 采用 S C D Z型垂直和 S C D P型水平速度传感器。 本次计算所选用的速度时程 曲线共两条，一条 来 自中深孔爆破强制崩矿地震效应测试 ，作为动载 荷一，如图 3所示；另一条来 自斜坡道掘进爆破斜 坡道掘进爆破地震效应测试，作为动载荷二，如图 4 所示。爆破震动速度时程曲线的主要参数为 1 动 载荷一的主频为 5 ． 6 4 Hz ，震动时问为 1 ． 0 0 S ，最大 幅值为 0 ． 1 8 1 2 m／ s ； 2 动载荷二的主频为 1 2 ． 2 8 H z ， 维普资讯 第 2 4卷托 I 6 K置 lI 等爆破 震曲 对乘 稳定性影响的 F L A C 。 分析 2 8 9 7 l 兰 l 一一二一 3中 孔爆破盎讪吐l 董 I I lih 线 Fi g 3 v e l o c i t y t i me d e p e n d e n l C L t F , e R me d i u m de p t h b o r e b l a s t i n g c o n c u s s i on 『 4 荆坡逋批逊爆破 动造埋时 [1} 线 Fi g 4 V e l o c i t y t i me e p e n d e m c u r v e s o t d e c l i n e d r i v i n g h l t L , t i n g v i b r a t i on 振动 m J - I’L lj 为 5 ． 1 5 s ．址人l i J 4 6 0 9 ri g s 4 ．2 ■破震动载荷施加 在 F LA C“ 动 力计赞 1 l ，动栽 倚输 入【 『J 以采 】 加 速度 时桴、速度 时 、f 移 n ,l 耵l 应力时 4千 1 t 方式 符果 粘滞 边 条件 必顽输入速度叫 进 行分析 。 f 1 F删试数掷 ● 质 振 动速J 篷，因此 动载 荷输 入选 川 _] 爆破震动i 虫崆叫 『 _ 卜测 试数 均 来 门采场或 甚道 的底扳．』 、 L 将墚破震动载 施 加 于采卒 区的底板 ，沿 z轴 负扪 0 播 5 模拟结果分析 黎 个汁符过程 分为静 J 汁算f 动 』 汁算 动 J j 分析 以静 计算 为前握 ，{l lj 在完成静力汁算的基础 } ，冉施 加动 力时程上 生 i’ 讣 尊 静力分析时 ，整个 采空区分 8个步骤 7 1 挖形成， Y 进行开挖 ， 每一步“挖 4m 每 ‘ 步升挖后剐 并争平衡状态， 卣到整个空区 开挖 完成 为止 将 婀蕞速度 时程通 过 袁的形式施圳 了 _． 采窄I 底板 ．舒圳进 行动 力汁算， 得到 力舒析结l粜。 5 ． 1 应力分析 阿 5 6分别为 采审区开挖后和施加爆破动载后 最大主麻 J 最小丰 力分 布图。可以霜 出，采空 I 毫歼挖后廊 』 发生重分 J ，边 界利捐 角处Ⅲ现 较大 』 J 鬃c l I ，秉卒 卜 匝 】 羁lJ l蓑扳 出现较大 的拉 应力． Ⅱ延伸至露 坑“ J 边坡脚I} 忖 近 罔岩 q - 最 人拉 』、 力 为 l 0 9 6 MP a ， J ‘ 岩体的抗批强度 0 ． 2 4 5 9 MP a ，右上 方 岩将 发 拉 裂破坏施加动 载后 ． 采 底板 处的拉 力集- 1 已经 小I蚶显，但 是右上 力⋯岩 I 仍然存在较大扎腑 力粲- 目_最 大拉 戍力 f J I 大 r岩体I’ 抗拉强度． 敞应采墩有敛的支护措施， 坝防采空 _ l ‘ 方旧岩失稳塌陷 动救荷一和动戴 苘二刈最人 E应 、最小 J ． J 分布 陶的影 响基本 类似，所以没有 t 目r I 具体进f T 分， 【 a J 挖 I b t咕茕 E 嗣后 } 5 最 凡j应 力静 电 “一 MP a F i 5 Di s t r i b u t i on o f ma x i mt l m p r i n c i p a l s t r e s s e s f u n i t MPa 强 蕾 I a I干托后 1 h J卸薮作用后 6 最小 1 山丹 n 】 c t f MP a F i c 6 Di s t fih u t i on o f rai n i n mm p r i n c i p a l s t r e s s e s { u n i t MPa ， 5 ．2 位移分析 闰 7为采空 区开挖和施 加动 载后围 位 移分布 冈 町见，最 大位 移出现 在采空【 上边 界，该边界 形状毫 墁 ／ f 圳划 ，儿已基本 彤成 “冠岩”形势。采空 开挖后阳 最大位移量大约 为4．7 5 mm，最 大位 移发 1 在打 f ‘ 与边界 l储昔处 ．同时底板 和包 括右上 方婊天边坡 内的 域 吼 了较 大的值移 量采 空 受动裁茼作用后 ，最大化移仍然发 生在右上 方 边 岩处．最 大位移量较静 茸 歼挖计算 结 疑，仃所增人 。采窀区受动拔 茼作用屙 ，『 岩最 、位移量为 5 ．4 0 mln．而受动救简 I 二 作用后，围岩 螭大位 移鼙选到 5 ．4 7 m m，增幅 与最大震动速度 成 一 ■■ ■ J ● ] r b _ _ ■ 维普资讯 2 8 9 8 _ _ 一 Jl 比 ．曲找 川．1l ．f 移 I pl； l 门l x域范 } ． 耍蜒 } 也 求 ⋯ ． I一 也 铺 “ 亦K 边坡 删线 6 结论与建议 荷． ．位移 l il 4 ‘门 域 范⋯ 1耍架 来 空【 t 1 ． ， J l if 板，f I 底扳处m 较静 计 t 尊小 这 说l ．施Ii f1 J j 城 ．小 Z J 人 岩 II{ I 最 犬位 移量，I L 变 } 的 ‘ l ℃ 范 l l 5 ．3 塑性区分析 『 8 震 卜 挖J _ I I 施JJ 艘 城 r,i塑陀 竹市 。从 [ 分 情况 求 ． 挖后型性 分 布较J 以 剪切 为土． 采 J} 和底扳l}『于 近 ．I_ I 露凡坑 fI J 边坡脚处 比H} 现型性l 似叫 ，采j 区受爆破动载 倚作用 ．型陆 分布 Y也m 比特 订 凡， 动j 戴茼 一作用I 比动哦 f f | H 』 j ．型恤 Ix范m 凡，说I ．IJ r 【 分 布 祝 【 lJ J ．／ -,J j 拔 , r Ib iii 值戚 I I 比 l I采- 【 仆挖 ．引起 麻 力重舒m和应 J 集t } I ；砸,IJ I l 爆破动载苘 ， 7 , i 蘸．s3 J _身发生_ 二 改 甄分m；⋯肯壕大似移 显情 夫，n增幅 鳅 大震动 建度 成 I- 爆艘震 ；叮 J fJ 小 仪增 太 ，} 协抖最大位 移量，而 l _ 改变 了n穆的 7- 布 式 1 7 h ． 2 3[ 连通 、范 Ⅲ人、 J 1 肜 j № 则 采空 lx J I 挖 性 范山较人． _ I 延伸 右 1 ． ， J 需灭坑边 垃喇l 施加爆破动钱简后． I l 会进步扩丈 ． 扩人峨度 功救甜幅位成 i l 比 I 3 复杂采空 Ix ． 爆破震动 倘我怍 F发生失 稳塌 的卮阶- 增加 4 避免爆破动 找荷引发 采 顶 板塌 陷的儿 基建 泌 ① 毓减少城 凡段药 ．阡 采 用敬羞爆破等臧 震措施。 ② 刘 ⋯观的塑性 『 ， 采取IJ贵m f 混凝 土、锚杆 、 锚索等有效 支 - 措施 ⑧ 对磷灭坑 i 边坡 进行雌测，必要时采取加 同擗施 参考文献I R e f e r e n c e s ．， I ． 1小 I ． ． 懈眭J t 锩删 扯址蔬 趸牟 } f t N l 。 1 I ,r j P I 报 ． 2 1 I I } 3 ． 2 2 I I J I 9 【 “I 9】 1 I Y a h Zhi xi nⅦ - ⋯ “ ； k ng h e- J i a nP i n o． e t a l,L L I d y c e a I uf h l a l i n d u c e d ’ m a n d b u i I d i n g l n r i l c n u [ ．I l C h i n e s e J o u n i a l t , f R, c k Me ；h a ni cs l ind En c a ec r i n g-2 X 3t 22 l I I 1 I 9 q7 I 9I I l i n Chi ne s e” Si n l l P K Bl a st v i b r a uu n d a mag e t o u nd e r t , z ou nd c oa l ndn e s f r om a d j a c e n t o p e n p i t b l a s t i n g ⋯l i l t c i l a t i l i n a l J o u rna l o r R k Mec ha n i c M i ni 。 Sc i e n c c 一 2 1 J l 1 2 门 9 50 9 7 3 Ya ng R I -R qu e P-Ka t s a b i L n i -P- I al M e a ui men l 0 『 训 { I hh 】 ， x ． ■_ _ 二 工 H ]曰■ ； ～ 帅 ㈨ ● ． 一 一一 一- 廿 _ m 扎 斯 8 g F 维普资讯 第 2 4卷第 l 6期 闩长斌等．爆破震动埘采 L父稳 定性影 响的 F L A C 。分析 2 8 9 9 【 4 】 【 5 】 【 6 】 【 7 】 【 8 】 【 9 】 【 l 0】 o f n e a r fi e l d b l a s t v i b r a t io n a n d d a ma g e [ J ] ．Ge o t e c h n ic al a n d En v ir o n me n t a l E n g i n e er i n g， 1 9 9 4， l 2 f 3 1 6 91 8 2 ． 言志信， 吴德伦，l ， 漪 ， 等地震效应及安全研究⋯．岩土力学， 2 0 0 2， 23 2 2 0 12 0 3 ． Y a n Zh i x i n，Wu De l u n，Wa n g Y i ，e t a 1 ． Th e s t u d y o f b l a s t v i b r a t i o n e ff e c t a n d s a f e t y [ J I ．Ro c k a n d S o i l Me c h a n i c s ，2 0 0 2．2 3 2 2 0 22 0 4 ． i n Ch i n e s e 谭忠盛，杨小林， 梦恕复线隧道施工爆破对既有隧道的影响 分析【 J 1 ．岩 力学 j 1 程 学报 ．2 0 0 3 ，2 2 2 2 8 1 2 8 5 T a n Zh o n g s h e n g，Y a n g Xi a o l i n ， W a n g Me n g s h u ． Ef f e c t o f b l a s t i n d o u b l e l i n e t u n n e l o n e x i s t i n g t u n n e l [ J ]Ch i n e s e J o u r n a l o f Ro c k Me c h a n i c s a n d E n g i n e e rin g，2 0 0 3 ，2 2 2 2 8 12 8 5 ． i n Ch i n e s e 刘亚群，李海波，赵 ．确定岩质边坡安全闽值的新方法⋯爆 炸与冲击， 2 0 0 4 ． 2 4 5 4 4 8 4 5 2 ． L i u Y a q u n ， L i H a i b o ， Z h a o J i a n A n e w me t h o d t o d e te r mi n e s a f e t y t h r e s h o l d v e l o c i t y o f s l o p e [ J ] Ex p l o s i o n a n d S h o c k Wa v e s ，2 0 0 4，2 4 5 44 84 5 2 ． i n Ch i n e s e 李宏男， 王炳乾， 林皋 爆破地震效应苦 r 题的探讨【 J J _ 爆炸 与冲击，1 9 9 6 ．1 6 1 6 67 2 ． f u H o n g n a n ．Wa n g B i n g q i a n ．L i n Ga o ． S t u d i e s o n s o me p r o b l e ms o n e x p l o s i o n e a r t h q u a k e r e s p o n s e s [ J ] E x p l o s i o n and S h o c k Wa v e s ， 1 9 9 6， 1 6 1 616 7 ． i n Ch i n e s e Ku h l e me y e r R L， L y s me r J ． F i n i t e e l e me n t me t h o d a c c u r a c y f o r wa v e p r o p a g a t i o n p r o b le ms [ J ] ． J ． S o i l Me c h ． and F o u n d a t i o n s Di v，AS CE， 1 9 7 3 ，9 9 S M5 4 2 14 2 7 ． I t a s c a Co n s u l t i n g Gr o u p，I n c ． F L AC 。 Fa s t L a g r a n g i a n An a ly s i s o f C o n t i n u a i n 3 D Ve r s i o n 2 ． 1 0． Us e r s Ma n u a l [ M] ．US A I t a s c a Co n s u lt i n g Gr o u p ， I n c ． ， 2 0 0 2 李廷春，李术／ j ，邱祥波 维快速拉格朗日法在安全顶板厚度 研究中fi cJ ／g 2 [ J ] ．岩 力学， 2 0 0 4 ， 2 5 6 9 3 5 9 3 9 r L i T i n g c h u n ． L i S h u c a i ，Qi u Xia n g b o Ap p l i c a t i o n o f 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 t o r e s e a r c h i n g o n s a f e r o c k c o v e r s o f Xi a me n s u b s e a t u n n e l [ J ] ．Ro c k a n d S o i l Me c h a n i c s ，2 0 0 4，2 5 6 9 3 59 3 9． i n Ch i n e s e 朱维申，李晓静，郭彦双，等．地下大型洞室群稳定性的系统性 [ 1 2 [ 1 3 J [ 1 4 】 【 1 5 【 1 6 】 研究【 J l l 岩 力学 j I 程学报，2 0 0 4 ，2 3 1 0 1 6 8 91 6 9 3 ． Z h u We i s h e n ，L i Xi a o j i n g ，Gu o Y a n s h u a n g ． e t a l S y s t e m a t i c a l s t u d y o n s t a b i l i t y o f l a r g e u n d e r g r o u n d c a v e r n s [ J I Ch i n e s e J o u rna l o f Ro c k Me c h a n i c s a n d E n g i n e e rin g， 2 X 4， 2 3 1 0 1 6 8 91 6 9 3 ． i n Ch i n e s e 刘春玲，祁牛文，童 强，等．利用 F L A C 分析某边坡的稳定 性【 J J _岩石力学 j 工程学报．2 0 04 ．2 3 1 6 2 7 3 02 7 3 3 ． L i u C h u n l i n g ，Qi S h e n g w e n T o n g L i q i a n g ．e t a l S t a b i l i t y a n a l y s i s o f s l o p e u n d e r e a nh q u a k e wi t h F LAC 。 【 J J l Ch i n e s e J o u rna l o f Ro c k Me c h a n i c s a n d En g i n e e r i n g．2 0 0 4， 2 3 1 6 2 7 3 02 7 3 3 ． i n Ch i n e s e 董玉书，徐友宁．』 坝铅锌矿 【 采J } 发引发的问题思考【 J 1 l华 南地 质 与矿产 ， 2 0 0 2 ， 4 2 3 2 7 ． D o n g Y u s h u ． Xu Y o u n i n g S o m e p r o b l e ms r e s u l t e d f r o m d is o r d e r e d mi n i n g o f l e a d z i n c o r e s i n Ch a n g b a o r e fi e l d [ J ] ． Ge o l o g y a n d Mi n e r a l Re s o u r c e o f S o u t h C