露天台阶爆破爆堆形态的PFC模拟.pdf

第2 9 卷第3 期 2 0 1 2 年9 月 爆破 B L A S T D 呵G V 0 1 ．2 9 N o ．3 S e p ．2 0 1 2 d o i 1 0 ．3 9 6 3 ／j ．i s s n ．1 0 0 1 - 4 8 7 X ．2 0 1 2 ．0 3 , 0 0 9 露天台阶爆破爆堆形态的P F C 模拟术 苏都都a , b ，严鹏a , b9 卢文波““，陈明“6 武汉大学a ．水资源与水电工程科学国家重点实验室；b ．水工岩石力学教育部重点实验室，武汉4 3 0 0 7 2 摘要尝试采用P F C2 D 数值方法预测台阶爆破的爆堆形态，并采用该方法对爆堆形态与爆破参数的关 系进行了研究。通过与相关观测、模拟结果的比较，论证了该研究方法的可靠性。模拟结果表明抛距随着 单耗q 和台阶高度日的增加而增大，随着抵抗线形的增大而减小；爆堆高度与各爆破参数的关系与抛距相 反。 关键词举堆形态’P F C2 D 方法； 台阶爆破 ，‘ 中图分类号T D 2 3 5 ．4 1文献标识码A文章编号1 0 0 1 4 8 7 X 2 0 1 2 } 0 3 0 0 3 5 0 7 P r e d i c t i o no fM u c k p i l eP r o f i l ef o rO p e nB e n c hB l a s t i n gw i t hP F C S UD U - d u 8 ⋯，Y A NP e n 9 8 一，L UW e n b 0 8 ⋯，C H E NM i n 9 8 ，6 a ．S t a t eK e yL a b o r a t o r yo fW a t e rR e s o u r c e sa n dH y d r o p o w e rE n g i n e e r i n gS c i e n c e ； b ．K e yL a b o r a t o r yo fR o c kM e c h a n i c si nH y d r a u l i cS t r u c t u r a lE n g i n e e r i n gM i n i s t r y o fE d u c a t i o n ，W u h a nU n i v e r s i t y ，W u h a n4 3 0 0 7 2 ，C h i n a A b s t r a c t T h eP F C2 Dn u m e r i c a lm e t h o dw a su s e dt op r e d i c tt h em u c k p i l ep r o f i l eo fb e n c hb l a s t i n ga n dt os t u d y t h er e l a t i o n s h i pb e t w e e nm u e k p i l ep r o f i l ea n db l a s t i n gp a r a m e t e r s ．B yc o m p a r i s o nt h er e l a t e do b s e r v a t i o nw i t I lt h e s i m u l a t i o nr e s u l t s ．t h i sm e t h o di sr e l i a b l e ．T h es i m u l a t i o nr e s u l t ss h o wt h a tt h et h r o w i n gd i s t a n c ee n l a r g e sw i t h 山e s p e c i f i cc h a r g eqa n dt h eb e n c hh e i g h tHa n dd e c r e a s e sw i t ht h eb u r d e n 彬i n c r e a s i n g ．T h eh e i g h to fm u e k p i l ep r o f i l e a n dt h eb l a s t i n gp a r a m e t e r si sc o n t r a r yt ot h et h r o w i n gd i s t a n c e ． K e yw o r d s m u c k p i l ep r o f i l e ；P F C2 D ；b e n c hb l a s t i n g 露天台阶爆破在矿山、水利、公路等工程中用途 极为广泛 2J ，随着各种土石方工程难度的增大和 爆破成本的增加，爆堆形态预测越来越成为爆破工 作者所关心的问题。爆堆形态是衡量爆破效果的重 要指标之一，它不仅反映了爆破参数、装药结构的合 理性，而且直接影响铲装、运输效率和经济效益。 针对这一问题，前人从物理数学方法、现场观 察、人工智能等方面做了许多研究，例如于亚伦提 出了台阶爆破爆堆形态预测的弹道理论模型和 收稿日期2 0 1 2 0 4 2 0 作者简介苏都都 1 9 8 8 一 ，男，硕士，主要从事水利水电工程施工 技术方面的研究工作， E m a i l s u d u d u l 5 0 4 0 1 q q ．t o m 。 基金项目国家重点基础发展规划 9 7 3 计划项目 2 0 1 0 C B 7 3 2 0 0 3 ； 国家自然科学基金资助项目 5 1 0 0 9 0 1 3 、5 1 1 7 9 1 3 8 W e i b u l l 模型旧o ；栾龙发等通过对现场台阶爆破过程 的高速摄影监测，探讨了深孔台阶爆破表面岩石的 移动规律H 1 ；赵春艳等应用遗传神经网络模型对台 阶压碴爆破效果进行了预测，并应用专家知识对台 阶爆破效果进行预测J 。然而，由于爆破作用过程 的复杂性和现有物理数学模型的不完善和其他模型 的局限性，目前绝大多数爆破工程实践中，对爆破效 果的预测还是依靠爆破设计人员个人的经验判断， 准确性及可靠性往往都难以满足要求。 随着计算机仿真技术的发展，数值模拟方法逐 渐成为研究爆破的主要工具。杨军、朱传云等采用 D D A 方法模拟分析了台阶爆破抛掷的过程，实现爆 生气体作用下岩石抛掷过程的计算机动态模 拟∞’刊；程晓君等采用3 d sm a x 软件实现了爆破过程 万方数据 3 6爆破2 0 1 2 年9 月 和爆堆形态可视化再现1 。 由于3 d sm a x 软件并不能真实反映爆破作用的 力学机理，D D A 方法需预先将岩体离散为规整的块 体，与实际情况存在一定差异，而颗粒流方法 P F C 采用细观离散元理论开发而成，将岩体离散为在最 大最小半径之间随机分布的二维圆盘或三维球形， 可以有效地反映岩体的微观本质，并且允许离散的颗 粒单元发生平移和旋转，颗粒也可以彼此全部分离， 所以它是研究爆破作用的理想工具。目前，该系列软 件已广泛应用于岩石力学行为模拟研究一’1 1 1 ，并被尝 试应用予爆破破岩机理方面的研究u 2 J 引，但应用于 台阶爆破爆堆形态模拟方面的研究尚不多见。 试图通过P F C 方法来预测台阶爆破的爆堆形 态，并研究台阶爆破参数对爆堆形态的影响。 1 P F C 方法中阻尼和粘结设置 P F C 系列软件以细观离散元理论为基础开发而 来，主要用于进行数值试验，将其用于爆堆形态的研 究，则需要结合爆破作用的特点，对影响岩体爆破过 程和爆堆形态的敏感参数，即颗粒系统阻尼和粘结 特性进行相关设置。 1 ．1 阻尼的设置 在颗粒流模型中，局部阻尼和粘性阻尼被引入 运动方程中来实现动能的消散，以使颗粒系统的动 能在一个合理的循环次数下达到一个稳定状态解。 根据局部阻尼的计算机理 具体见参考文献 [ 1 4 ] ，局部阻尼对只受重力作用的抛掷运动和在 抛掷过程中相互碰撞的颗粒是不适合的。因此，为 了使得能量消散过程更好的反应实际情况，比较好 的方法是在爆炸荷载施加结束，岩体破碎并开始运 动时，将待爆区颗粒的局部阻尼系数设置为零，并激 活粘性阻尼，即将局部阻尼置于不激活状态，使粘性 阻尼发挥阻尼作用。 如此，系统中的颗粒在抛掷过程中不受到阻尼 作用，而相互碰撞颗粒之间接触处的相对运动受到 阻尼作用，并产生能量损失。 1 ．2 粘结的设置 P F C 系列软件采用离散元理论，克服了有限元 软件单元不能脱离的问题，但根据文献[ 1 4 ] 可知， 在P F C 程序中．，颗粒分离后，在计算过程中如果颗 粒之间距离足够小时，将重新构成新的接触，这明显 与爆破作用过程中，岩块在抛掷过程中相互碰撞，使 岩块进一步破碎特性不符。 为更真实地反映爆破作用，使模拟过程更接近 实际，模拟过程中在岩体完成破碎并开始移动时，将 待爆区岩体的粘结半径和粘结力设置为零。 2 爆炸荷载的简化及施加 目前对爆炸荷载的模拟方法主要有爆炸荷载 曲线和流固耦合2 种方式，其中爆炸荷载曲线由于 简单方便被广泛使用，采用给炮孔壁粒子施加爆炸 荷载曲线的方式来模拟爆破作用。 根据凝聚炸药爆轰波的c J 理论，耦合装药条 件下炮孔壁上的初始平均爆轰压力为 p 一型⋯ ‘o 一2 7 1 、17 式中P 0 为炸药的爆轰压力；p 。、D 分别为炸药的密 度和爆轰速度；y 为等熵指数。 对于不耦合装药，若装药时的偶合系数b ／a 值 较小，则爆生气体的膨胀只经过P P 。一个状态， 此时炮孔初始平均压力P 0 为 P 0 2 考％1 护 2 1 0 一 1 ， ＼i 』 卜7 式中口为装药直径；b 为炮孔直径。 根据爆破工程广泛采用的多方气体状态方程可 以得到初步膨胀以后的炮孔压力 P ￡ 2 高 P o 3 式中K 为炮孔初始体积；7 1 ．4 。 裂缝贯穿瞬间，相邻炮孔连线上的裂缝呈锥形 扩展，如图1 所示。结合文献[ 1 5 ] 的研究成果，取t ， 为2m s ，t 。为7m s 。爆炸荷载作用过程的曲线形式 如图2 所示，P 的为爆炸荷载峰值，t ，为爆炸荷载上升 时间，t a 为爆炸荷载正压作用时间。 声 j H - 墨盅 j _ _ E 蛊 图1 炮孔连线上裂缝扩展图 F i g ．1 C r a c kp r o p a g a t i o nm a p t ，b f 图2 爆炸荷载曲线 F i g ．2 B l a s tl o a dc u r v e 3 爆堆形状的P F C 模拟 3 ．1 材料P F C 微参数的选择 与有限元软件采用宏观力学参数不同，P F C 软 万方数据 第2 9 卷第3 期苏都都，严鹏，卢文波，等露天台阶爆破爆堆形态的P F C 模拟3 7 件采用细观参数以表征颗粒间的力学行为。本节模花岗岩的P F C 微参数如表1 所示。 型所用岩体材料参数取花岗岩参数，根据文献[ 3 ] 表1 花岗岩宏、微观参数 T a b l e1M a c r oa n dm i c r o - p a r a m e t e r so fl a cd ub o n n e tg r a n i t e 宏观参数取值微观参数取值 微观参数取值 p ／ g c m 。3 2 6 3 0 p ／ g c m ‘3 2 6 3 0 A l E ／G P a6 9 R 。．，／R 。i 。 1 ．6 6 E ／G P a 6 2 o r ／M P a1 5 0 E 。／G P a 6 2 七4 ／∥ 2 ．5 o r “／M P a 9 0 l j } n ／％，2 ．5 盯 f ． 1 5 7 口0．26 p 0 ．5 o r r 3 6 3 ．2 、岩体破碎时间确定 根据王文龙等的实测资料【l6 | ，巷道爆破中，单 自由面条件下，岩体开始运动的时间为6 ～5 8m s ； 双自由面条件下为3 2 7m s 。F e l i c eJJ 与P r e e c eD s 对中等台阶高度的岩石爆破高速摄影观测结果表 明Ⅲ’埔1 待爆区岩体脱离母岩的时间约为几十毫 秒。张奇对爆破过程中岩石破碎速度的分行理研究 结果也表明【19 1 ，当抵抗线为0 ．5 5m 时，岩体的破 碎时间为5 ～5 0m s 。以上研究表明，岩块开始运动 的时间与自由面、台阶高度和抵抗线等参数有关。 且其长短主要受台阶高度控制，台阶高度越大，滞后 时间越长。 由于模拟的工况为深孔爆破，最大台阶高度为 8 ．0m ，因此取岩体破碎运动时间为1 5m s 。 3 ．3 排间延时时间确定 确定合适的微差间隔时间和准确地控制此间隔 时间是做好微差爆破的关键㈣J 。 段海峰提出了爆破机理的推墙假说和回弹假 说瞄1 1 ，并认为合理的微差时间应大于前排孔产生足够 补偿空间的时间，小于后排孔回弹结束时间，并介绍了 一种简单的确定微差时间的办法。王戈等通过工程实 际和数值试验认为旧J ，合理的排间微差1 0m 。 一般矿山爆破工作中实际采用的微差间隔时间 为1 5 ～7 5m s ，通常用1 5 ～3 0m s 。国外一般认为微 差间隔时间只要大于1 0m s 即可旧J ，排间微差间隔 时间可达到2 0 0m s 。在模拟过程中取排间微差时 间为6 0m s 。 3 ．4 单孔台阶爆破的P F C 模拟 选取典型的单排台阶爆破模型，如图3 所示。 在本例中，取炮孔深L 8 ．0m ，装药长度L 。 5 ．5m ，堵塞长度厶 2 ．5m ，炮孔直径为b 1 1 0 ／I l i a ，采 用密度P 1 0 0 0k g ／m 3 ，爆轰波速D 4 0 0 0m ／s 的乳化 炸药，荷载峰值P o 2 0 0 0M P a 。为简化计算，假设裂缝 的平均宽度为2c m 女g 图1 所示 ，则裂缝贯穿的瞬间 边界上的气体压力为P I 9 5 ．2M P a 。 f “i 【二生一 图3 台阶爆破模型图 单位m F i g ．3 B e n c hg r a t i n gm o d e ld i a g r a m u n i t m 图4 给出了爆破过程的模拟结果，整个爆破过 程共历时1 ．8 7s 。由图4 可以看出，在爆炸荷载作 用下，裂缝开始出现。由于台阶上部堵塞段不加载 而且台阶下部抵抗线较大，从而台阶中部出现鼓起 突出。最后台阶上下两部分由于速度差的存在出现 了脱离，台阶上部分由于下部速度大上部速度小而 向上旋转，台阶下部岩体受重力作用下落。．这一过 程与Y o u j u nN i n g 等采用D D A 对台阶爆破的模拟过 程和栾龙发等采用高速摄影手段，对台阶爆破岩石 移动规律的研究成果相似【4 ．6 】。 图5 为杨军采用D D A 方法模拟的单孔台阶爆破 爆堆形态∞J ，通过对比可以发现，其结果与所用模型 的模拟结果基本相同，爆堆接近炮孔处的不同是由于 后冲向岩体破坏程度不同造成的。采用文献[ 3 ] 所提 出的弹道理论模型 具体公式参考文献[ 3 ] 计算单 排孔工况的爆堆形态，其相关参数也与本节模拟结果 基本符合，具体结果如表2 所示。弹道法抛矩大并且 高度高的原因是由于其假设爆堆形态的基本形状为 三角形造成的。这些都证明本文采用P F C 方法来模 拟台阶爆破爆堆形态思路的可行性。 万方数据 爆破2 0 1 2 年9 月 正匝瞳 瞌』瞌一 图4 单排台阶爆破过程模拟图 F ；g ，4 , q i m u l a t e db e n c hb l a s t i n gp r o c e s s 图5D D A 计算结果图 F i g ．5 T h es i m u l a t e dr e s u l to fD D A 裹2 理论结果与实际结果对照表 T a b I e2S i m u l a t e da n db 叫陆I d cr e s u l t s 3 ．5 双排孔台阶爆破 双排孔模型亦采用表1 的参数，模型相关尺寸 如图6 所示。采用上述的计算模型和相应的参数进 行计算，图7 依次给出了0 ．0 0 0s , 0 ．0 1 4s , 0 ．0 9 78 、 0 ．2 1 5s , 0 ．8 2 1s 和1 ．9 9 7S 时刻的爆破过程图。由 图7 b 可知，第l 排孔首先起爆，在爆炸荷载作用 下，前冲向岩体向自由面移动，岩体被拉裂，裂纹开 始出现，后冲向岩体向相反方向运动，第2 排孔被压 缩。由图7 c 可知，第2 排孔起爆，在爆炸荷载作 用下，爆腔增大，第2 排孔抵抗线内岩体开始向自由 面方向运动，由于台阶上部堵塞段不加载而且台阶 下部抵抗线较大，因此、台阶中部速度最大，从而鼓 起突出。 由图7 d 可知，岩体继续向自由面方向移动，2 次爆破岩体在运动中碰撞。由图7 e 可知，受重力 作用岩块继续下落，并产生堆积。台阶爆破爆堆形 态的最终结果如图7 f 所示。 图6 双孔台阶爆破模型图 单位m F i g ．6 B e n c hb l a s t i n gm o d e ld i a g r a m u n i t m 4 单排炮孔爆堆形态与爆破参数关系 针冤 改变爆破条件，合理选择岩体参数，则可利用 P F C 软件来预测不同爆堆形态，研究其规律，以便进 行合理的台阶爆破设计。本节采用3 ．4 节台阶爆破 模型作为基本模型，通过分别改变相应的爆破参数， 来研究该参数对爆堆形态的影响。 共讨论了9 种工况，分别比较了不同抵抗线 形 2 ．5m 、3 ．0m 、3 ．5m 、炸药单耗 q 0 ．6k g ／m 3 、 0 ．6 5k g ／m 3 、0 ．7k g ／m 3 和台阶高度 H 8 ．0m 、 1 0 ．0m 、1 2 ．0m 条件下爆堆的形态变化情况。 为便于比较不同参数条件下的爆堆形态，对各 种工况的计算结果进行了简化处理。首先以模型中 同一点为原点，描出各工况最终爆破爆堆形态轮廓 曲线，如图8 所示；再将其换算为统一的比例，置于 同一张图中进行对比。 万方数据 第2 9 卷第3 期苏都都，严鹏，卢文波，等露天台阶爆破爆堆形态的P F C 模拟 3 9 h 卜 1 ．H ls{ i ’l ，一lu u 、 图7 双孔台阶爆破过程图 F i g ．7S i m u l a t e db e n c hb l 越t mp ／D C e 8 8 图8 爆破形态的简化处理 F i g ．8S i m p l i f i c a t i o no fM u c k p i l eP r o f i l e 4 ．1 爆堆形态与抵抗线关系 抵抗线是台阶爆破的重要参数之一，利用3 ．4 节中建立的模型，模拟了抵抗线形为2 ．5m 、3 ．0m 、 3 ．5m3 种工况 分别记为W - 2 ．5 、弘3 ．O 、肛3 ．5 下 的爆堆形态。其它参数与第3 节中相同，模拟结果 如图9 所示，爆堆的相关特征参数如表3 。 g 5 ．0 藿2 ．5 鲻 瞪 02 ．55 ．07 ．51 0 ．01 2 ．51 5 ．0 抛掷距离L ／m 图9 不同抵抗线工况下爆堆形态曲线 F i g ．9 T h em u c k p i l ep r o f i l eo fd i f f e r e n tr e s i s tl i n e 从图9 和表3 可以看到，随着抵抗线形的增 加，岩石抛掷距离减小，爆堆高度增大，整个爆堆形 状由扁平向高尖过渡。这是因为随着抵抗线的增 加，台阶爆破需要破碎和抛掷的岩体体积增加，导致 抛距减小，但同时由于被爆岩体的体积增加，岩体破 碎和抛掷需要更长时间，导致更高比例的爆破能量 传给岩体，使爆堆的高度有所增加。 表3 爆堆特征参数 T a b l e3T h ec h a r a c t e r i s t i cp a r a m e t e r so fm 删 吣p r o f i l e 4 ．2 爆堆形态与炸药单耗口的关系 ．炸药单耗q 是台阶爆破中的重要指标。它是指 爆破1m 3 岩石所需要的炸药量。本小节以3 ．4 节模 型为基础，分别计算q 等于0 ．6k g ／m 3 、0 ．6 5k g ／m 3 、 0 ．7k g ／m 33 种工况 分别记为q ．0 ．6 、g - 0 ．6 5 、 g 一0 ．7 下的爆堆形态，计算结果和相关特征参数分 别如图1 0 和表4 所示。 口5 ．0 蓑2 ．5 鳞 躜 7 ．51 0 ．o1 2 ．51 5 ．01 7 ．5 抛掷距离L ／m 图1 0 不同炸药单耗工况下爆堆形态曲线 F i g ．1 0 T h em u c k p i l ep r o f i l eo fd i f f e r e n t e x p l o s i v ec o n s u m p t i o n 表4 爆堆特征参数 T a b l e4T h ec h a r a c t e r i s t i cp a r a m e t e r so fm u c k p f l ep r o f i l e 万方数据 爆破 2 0 1 2 年9 月 通过图1 0 和表4 可以得出以下结论随着炸药 单耗q 的增加，抛掷距离增大，爆堆高度减小，爆堆 的最高点逐步远移。这是由于随着炸药单耗的增 加，爆炸荷载峰值增大，使得岩石抛掷距离增大，爆 堆的高度相应减小。 4 ．3 爆堆形态与台阶高度日的关系 台阶高度日是台阶爆破设计的重要参数，不同 的台阶高度不仅影响爆破生产效率而且对爆堆形态 也有比较太的影响。本小节以3 ．4 节模型为基础， 分别计算H 等于8 ．0m 、1 0 ．0m 、1 2 ．0I l l3 种工况 分别记为皿8 ．0 、H - 1 0 。0 、H - 1 2 ．0 下的爆堆形态， 计算结果及特征参数分别如图1 1 和表5 所示。 { 划 惶 辫 蹬 图1 1 不同台阶高度工况下爆堆形态曲线 F i g ．11 T h em u c k p i l ep r o f i l eo fd i f f e r e n tb e n c hh e i g h t 表5 爆堆特征参数 T a b l e5T h ec h a r a c t e r i s t i cp a r a m e t e r so fm u c k p i l ep r o f i l e 通过图1 1 和表5 的分析得出以下结论随着台 阶高度胃的增加，岩石抛掷距离逐渐增大。与抵抗 线、炸药单耗相比，台阶高度对爆堆形态的影响不 大，爆堆高度和抛掷距离分别有所增加。这是因为 由于台阶高度的增加，底盘抵抗线随之增大，下部岩 体抛掷作用减弱，导致爆堆最高点向台阶移动，但随 着台阶高度的增加，岩石抛掷运动时间有所增大，上 部岩体抛掷距离有所增加。 5 结论与展望 1 P F C 数值方法预测结果与相关模拟结果和 理论模型计算结果有较好的一致性，是预测露天台 阶爆破爆堆形态行之有效的方法。 2 岩石抛掷距离随着抵抗线形、炸药单耗q 和台阶高度日的增加而增大；爆堆高度随抵抗线形 的增加而增大，随q 的增大而减小；并且爆堆形态对 单耗q 与抵抗线形比对台阶高度日更为敏感。 3 实际岩体比较复杂，如何更好地反映结构 面对爆堆形态的影响仍需进一步研究。 参考文献 R e f e r e n c e s [ 1 ] 胡贵如．露天台阶爆破C A D 开发模型的研究[ J ] ．有 色金属 矿山部分 ，2 0 1 0 ，6 2 2 6 0 - 6 3 ． [ 1 ] H UG u i - r u ．R e s e a r c ho fo p e n p i tb e n c hb l a s t i n go fC A D d e v e l o p m e n tm o d e l [ J ] ，N o n f e r r o u sM e t a l s M i n i n gS e c - t i o n ，2 0 1 0 ，6 2 2 6 0 6 3 ． i nC h i n e s e [ 2 ] 张正宇．现代水利水电工程爆破[ M ] ．北京中国水利 水电出版社，2 0 0 3 ． [ 3 ]于亚伦，高焕新，张云鹏，等．用弹道理论模型和W e i b u l l 模型预测台阶爆破的爆堆形态[ J ] ．工程爆破，1 9 9 8 ， 4 2 1 9 - 2 2 ． [ 3 ] Y uY a - l u n ，G A OH u a n x i n ，Z H A N GY u n p e n g ，e ta 1 ．P r e - d i c t i o no fm u c k p i l ep r o f i l ef o rb e n c hb l a s t i n gw i t hb a l l i s t i c sm o d e la n dW e i b u l l ’Sd i s t r i b u t i o nm o d e l [ J ] ．E n g i ． n e e f i n gB l a s t i n g ，1 9 9 8 ，4 2 1 9 - 2 2 ． i nC h i n e s e [ 4 3 栾龙发，庙延钢．深孔台阶爆破岩石移动规律的研究 [ J ] ．中国工程科学，2 0 0 5 ，7 s 1 2 4 8 - 2 5 1 ． [ 4 ] L U A NL o n g f a ，M I A OY a n g a n g ．T h er o c km o v i n gr e g u l a r i t yi nd e e ph o l eb e n c he x l j l o s i o n [ J ] ．E n g i n e e r i n gS c i e n c e ，2 0 0 5 ，7 S 1 2 4 8 - 2 5 1 ． i nC h i n e s e [ 5 ] 赵春艳，常春，张继春，等．台阶压碴爆破效果遗传 神经网络预测[ J ] ．岩土力学．2 0 0 3 ，2 4 S 1 8 8 - 9 0 ． [ 5 ] Z H A OC h u n y a n ，C H A N GC h u n ，Z H A N GJ i e h u n ，e ta 1 ． P r e d i c t i o no fb e n c ht i g h tb l a s t i n ge f f e c t st h r o u g hg e n e t i c n e u r Mn e t w o r k [ J ] ．R o c ka n dS o i lM e c h a n i c s ，2 0 0 3 ， 2 4 S 1 8 8 - 9 0 ． i nC h i n e s e [ 6 ] N I N GY o - j u n ，Y A N GJ u n ，M AG u o - w e i ，e ta 1 ．M o d e l l i n g r o c kb l a s t i n gc o n s i d e r i n ge x p l o s i o ng a sp e n e t r a t i o nu s i n g d i s c o n t i n u o u sd e f o r m a t i o na n a l y s i sfJ ] ．R o c kM e c hR o c k E n g ，2 0 1 1 4 4 4 8 3 - 4 9 0 ． [ 7 ] 朱传云，戴晨，姜清辉．D D A 方法在台阶爆破仿真模 拟中的应用[ J ] ．岩石力学与工程学报，2 0 0 2 ，2 1 s 1 2 4 6 1 ．2 4 6 4 ． [ 7 ] Z H UC h u a n - y u n ，D A IC h e n ，J I A N GQ i n g - h u i ．N u m e r i c a l s i m u l a t i o no fb e n c hb l a s t i n gb yd i s c o n t i n u o u sd e f o r m a t i o n a n a l y s i sm e t h o d [ J ] ．C h i n e s eJ o u r n a lo fR o c kM e c h a n i c s a n dE n 6 n e e f i n g ，2 0 0 2 ，2 1 S 1 2 4 6 1 - 2 4 6 4 ． i nC h i - n e s e [ 8 ] 程晓君，朱传云．3 D SM A X 在台阶爆破计算机模拟中 的应用[ J ] ．中国农村水利水电，2 0 0 6 1 0 8 7 8 9 ． [ 8 ] C H E N GX i a o - j u n ，Z H UC h u a n y u n ．A p p l i c a t i o no f3 D S M A Xi nb e n c hb l a s t i n gfJ ] ．C h i n aR u r a lW a t e ra n dH y ， d r o p o w e r ，2 0 0 6 1 0 8 7 - 8 9 ． i nC h i n e s e [ 9 ] 。P O T Y O N D YDO ，C U N D A L LPA ．Ab o n d e d p a r t i c l em o d e l f o rr o c k 『J1 ．I n t e r n a t i o n a lJ o u r n a lo fR o c kM e c h a n i c sa n d M i n i n gS c i e n c e s ，2 0 0 4 ，4 1 8 1 3 2 9 1 3 6 4 ． 万方数据 第2 9 卷第3 期苏都都，严鹏，卢文波，等露天台阶爆破爆堆形态的P F C 模拟 4 l 上接第3 0 页 卢文波．拆除爆破中裸露钢筋骨架的失稳模型[ J ] ．爆 破，1 9 9 2 ，9 2 3 1 3 4 ． L UW e n b o ．Am e c h a n i c sm o d e lo nt h el o s ss t a b i l i t yo f e o l u F n n so fr e i n f o r c e dc o n c r e t ef r a m ei nb l a s t i n gd e m d i t i o n [ J ] ．B l a s t i n g ，1 9 9 2 ，9 2 3 1 - 3 4 ． i nC h i n e s e 关志中．用控制爆破拆除2 l 米高钢筋混凝土框架大楼 [ C ] ／／爆破与安全．武汉湖北科技出版社，1 9 8 4 1 2 1 6 ． 冯叔瑜．城市控制爆破[ M ] ．北京中国铁道出版社， 1 9 8 5 5 2 - 5 8 ． 李肇胤，袁澍萱．拆除爆破中钢筋混凝土框架柱的失 稳判据[ J ] ．爆破，1 9 9 3 ，1 0 2 1 0 9 1 1 2 ． L IZ h a o - y i n ，Y U A NP e n g x u a a ．C r i t e r i o no fl o s ss t a b i l i t y f o rc o l u m n so fRCf r a m e si nb l a s t i n gd e m o l i t i o n 『J ] ． B l a s t i n g ，1 9 9 2 ，l O 2 3 l 一3 4 ． i nC h i g e s e 吴枫．框架结构爆破拆除力学模型与数值分析[ D ] ． 天津天津大学，2 0 0 3 1 2 1 6 ． W UF e n g ．N u m e r i c a la n a l y s i s ＆m e c h a n i c a lm o d e la n a l y s i so fb l a s t i n gd e m o l i t i o no nf r a m e w o r ks t r u c t u r e 『D ] ． T i a n j i n g M a s t e rd e g r e et h e s i so f T i a n j i nU n i v e r s i t y ，2 0 0 3 1 2 ．1 6 ．fi nC h i n e s e 龚相超，胡百鸣．液压缸稳定性设计的模糊可靠性计 算[ J ] ．武汉科技大学学报，2 0 1 1 ，3 4 5 3 4 2 3 4 4 ． G O N GX i a n g - e h a o ，H UB a i - m i n g ．F u z z yc a l c u l a t i o no fh y - d r a u l i ec y l i n d e rr e l i a b i l i t ya b o u ts t a t i cs t a b i l i t y [ J ] ．J o u r - h a lo fW u h a nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y ， 2 0 1 1 。3 4 5 3 4 2 - 3 4 4 ． i nC h i n e s e 陈大勇，杨军．大型薄壁钢筋混凝土连体筒仓群爆 破拆除[ J ] ．工程爆破，2 0 1 2 ，1 8 1 5 2 - 5 5 ． C H E ND a y o n g ，Y A N GJ a n ．B l a s t i n gd e mo l i t i o no fl a r g ei n - t e m o a n e c t e dm i n f o m e dc o n c r e t es i l o sw i t ht h i nw a U I j ] ．E n g i n e e r i n gB l 碰I l g ，2 0 1 2 ，1 8 1 5 2 - 5 5 。 i nC h i n e s e 言志信，刘培林，叶振辉．框架一剪力墙结构建筑物定 向倒塌的研究[ J ] ．爆破，2 0 1 1 ，2 8 2 5 - 8 ． Y A NZ h i - x i n ，L I UP e i - l i n ，Y EZ h e n - h u i ．R e s e a r c ho no r i ‘ e n t a t i o nc o l l a p s eo ff r a m e - - s h e a rs t r u c t u r eb u i l d i n g [ J ] ． B l a s t i n g ，2 0 1 1 ，2 8 2 5 - 8 ． i nC h i n e s e 谢春明，杨军，薛里．高耸筒形结构爆破拆除的数