钢筋混凝土双烟囱交叉折叠倒塌爆破数值模拟.pdf

第2 8 卷第4 期 爆破 V o l 2 8 N o 4 2 0 1 1 年1 2 月B L A S T I N G D e c ．2 0 1 l 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 1 ．0 4 ．0 0 7 钢筋混凝土双烟囱交叉折叠倒塌爆破数值模拟 罗倦／, 尘I ，池恩安2 ”，詹振锵3 ，王丹丹3 ，张修玉3 1 ．贵州久联民爆器材发展股份有限公司，贵阳5 5 0 0 1 ；2 ．贵州新联爆破工程有限公司，贵阳5 5 0 0 0 2 ； 3 ．贵州大学矿业学院，贵阳5 5 0 0 0 3 摘要针对2 座钢筋混凝土烟囱爆破拆除倒塌场地长度小于烟囱高度的特殊情况，拟采用双向折叠交叉 倒塌的爆破技术方案。为了对所采取方案进行验证，在爆破前首先采用A N S Y S 建立分离式共节点实体模 型，并对有限元网格进行划分，然后在D Y N A 环境下完成对烟囱倒塌过程的数值模拟计算，模拟结果表明所 采取的方案可行、有效。实际烟囱倒塌过程、倒塌长度与数值模拟结果较为吻合。因此可认为拆除爆破前，采 用数值模拟方法对爆破方案进行验证，并优化爆破设计方案具有重要意义。 关键词钢筋混凝土双烟囱；爆破拆除；数值模拟 中图分类号T U 7 4 6 ．5文献标识码A文章编号1 0 0 1 4 8 7 X 2 0 1 1 0 4 0 0 2 7 0 4 N u m e r i c a lS i m u l a t i o no fB l a s t i n go fT w oR e i n f o r c e dC o n c r e t e C h i m n e y w i t hC r o s sF o l d i n gC o l l a p s e L U OD e 厕1 ，C H IE n a n 2 ⋯，Z H A NZ h e n q i a n 9 3 ，W A N GD a n d a n 3 ，Z H A N GX i u y u 3 1 ．G u i z h o uJ i u l i a nI n d u s t r i a lE x p l o s i v eM a t e r i a l sD e v e l o p m e n tC oL t d ，G u i y a n g5 5 0 0 0 1 ，C h i n a ； 2 ．G u i z h o uX i n l i a nB l o tE n g i n e e r i n gL i m i t e dC o r p ，G u i y a n g5 5 0 0 0 2 ，C h i n a ； 3 ．S c h o o lo fM i n i n g ，G u i z h o uU n i v e r s i t y ，G u i y a n g5 5 0 0 0 3 ，C h i n a A b s t r a c t A i ma tt h es p e c i a ls i t u a t i o no ft h el e n g t hf o rc o l l a p s i n gl e s st h a nt h eh e i g h to fc h i m n e y 。t l l eb l a s t i n g s c h e m eo fc “瑚f o l d i n gc o l l a p s eW a su s e dt od e m d i s ht h et w or e i n 董o r e ee o n c l q e t ec h i m n e y ．B e f o r et h eb l a s t i n g ．f o rt e s - r i n ga n dv e r i f y i n gt h es c h e m e ，t h eA N S Y Sw a su s e dt oe s t a b l i s ht h es e p a r a t ee n t i t ym o d e lo fn o d e sa n dm e s h i n g ，t h e n s i m u l a t i n gt h ec o l l a p s eC O U E ％o ft h ec h i m n e yi nD Y N A ．T h es i m u l a t i o nr e s u l ts h o w st h a tt h i ss c h e m ei sf e a s i b l ea n d v a l i d ．A f t e rb l a s t i n g 。t h ec o l l a p s ec 0 岫ea n dt h ec o l l a p s el e n g t ho f c h i m n e y 删c o i n c i d ew i t ht h es i m u l a t i o nr e s u l t ，8 0 u s i n gt h es i m u l a t i o nt e s t sa n dv e r i f y sa n do p t i m i z e st h eb l a s t i n gs c h e m eb e f o r eb l a s t i n gi sv e r yi m p o r t a n t ． K e yw o r d s t w or e i 坷o r e e dc o n c r e t ec h i m n e y ；e x p l o s i v ed e m o l i t i o n ；n u m e r i c a ls i m u l a t i o n 爆破拆除在国内外高耸建筑物拆除工程中得到 广泛应用，技术上日益完善，但其理论仍严重落后于 实践。如果烟囱爆破拆除前能采用数值模拟的方法 模拟烟囱爆破拆除全过程，将会为爆破拆除提供有 力的理论技术支持，从而保证烟囱拆除爆破的安全 可靠引。针对2 座钢筋混凝土烟囱倒塌空间有限 收稿日期2 0 1 l 0 9 1 6 作者简介罗德丕 t 9 6 5 一 ，男，贵州贵阳人，高级工程师，华中科技 大学博士研究生，主要从事爆破技术研究工作， E m a i l n l i “g s h e n g z I I ∞ 1 6 3 ．c o m 。 的情况，采用交叉折叠倒塌的技术方案，并利用有限 元程序软件L S D Y N A 对爆破拆除方案进行模拟 预演，以期找到更为科学、安全的爆破拆除设计。 1 工程概况 某厂有2 座8 5m 高的钢筋混凝土烟囱需要爆 破拆除。烟囱周围的环境极其复杂，其中东面为厂 区道路，南面1 0 ．7m 为电除尘室，西面6 5m 处有南 北走向的架空管道，北面5 0m 处有东西走向的架空 万方数据 爆破2 0 1 1 年1 2 月 管道，2 座烟囱相距1 2m 。l 。烟囱到西北角的最远 距离为8 6m ；2 。烟囱到西南角的最远距离为8 3m ， 到西北角的最远距离为7 9m ，可供烟囱倒塌的场地 极其有限。2 座烟囱的结构完全一样，底部外周长 2 1 ．5m 、烟囱筒身底部外径6 ．8 4m 、壁厚0 ．3 5m 、空 气隔热层0 ．0 5m 、内衬0 ．1 2m ；顶部外径3 ．1 2m 、 顶部周长9 ．8m ，东面各有一高2 ．0m 、宽1 ．1m 的 出灰口。见图l 。 公路 图l 双烟囱四邻环境示意图 单位m l r , g ．1T h e 鲺m Ⅻ蛾e n v i r o r a n e n to ft h et w oe h i r r 瑚y a u n i t m 2 爆破方案 由于该厂是连续化生产，业主要求在烟囱爆破 的施工准备和实施阶段其毗邻的各厂房的生产不能 停工，生产工人不能撤离，这就对工程提出了很高的 技术要求和安全要求，在确保2 座烟囱在按设计方 向倒塌的同时，还必须完全杜绝爆破飞石、震动、粉 尘危害。因此针对特殊的周围环境、特殊的工程要 求，通过数值模拟制定科学的技术方案是很重要的， 并在技术方案基础上对各爆破参数合理选取。 1 。烟囱选择向西北定向倒塌的爆破技术方案， 为了确保烟囱倒塌触地后其爆碴不危及到架空管道 的安全，采用抬高爆破切口的施工技术。爆破切口 选在 l O ．5m 处，即 1 0 ．5m 以上部分采用定向爆 破拆除，0 ．0 0 ～ 1 0 ．5m 部分待上部爆破后采用机 械拆除。2 ’烟囱选择双向两折叠定向倒塌爆破技术 方案，其中 2 1 ．0 0m 以上部分向西南方向倒塌， 0 ．0 0 一 2 1 ．0 01 1 1 部分向东北方向倒塌，开2 个爆 破切口，第1 个爆破切口在 2 1 ．0 0m 处，切口高度 2 ．4m ，第2 个切口在 0 ．5 0m 处，切口高度3 ．01 1 3 。 2 。烟囱第2 个爆破切口滞后于第1 个爆破切口1 ．2s 起爆，1 ’烟囱的爆破切口滞后于2 。烟囱第2 个爆破 切口68 起爆，2 座烟囱的倒塌方向垂直交叉。 工程的重大技术风险主要是2 。烟囱上部分能 否按设计方向倒塌，爆破后瞬间失稳很关键，因此在 技术设计中对2 ’烟囱的切口圆心角的取值比较大 2 1 6 。 ，切口高度取了2 ．4m ，定向窗也开的比较 大 底长1 ．51 1 1 ，高2 ．4m ，这为2 。烟囱爆破后尽快 形成失稳倒塌趋势创造了条件。但是由于其结构特 点，如爆破不能达到预期的效果，对周围环境的安全 影响较大。 3 烟囱整体倒落数值模拟过程 3 ．1 模型简化 在数值模拟前首先对烟囱进行模型简化 1 烟 囱为均质混凝土材料； 2 除爆破切口外，倾倒部分烟 囱为刚体； 3 烟囱在倾倒过程中满足动量守恒； 4 计算中只考虑平面运动形式，不考虑碰撞后烟囱 转动以外的其他运动情况。根据爆破方案设计、烟囱 的结构特征、混凝土的材料参数，合理确定炸药爆炸 产生的载荷对模拟结果的准确性有着非常重要的作 用，包含确定爆破激振力的大小、作用位置和方向、作 用时刻和持续时间等对模拟结果的影响h6 。。 3 ．2 模拟倒塌过程 根据2 ’烟囱在双向折叠采用的是“先上后下” 的起爆顺序，在计算模型中，在计算机处理过程中， 0 ．3s 时对上切口处块体单元首先作用1 个冲击载 荷，1 ．5s 时，再在下缺口块体单元上作用另1 个冲 击载荷。6s 时，对1 。烟囱缺口块体单元上作用另一 个冲击载荷。下列图2 为2 座钢筋混凝土烟囱倒塌 整体过程。模拟中体现了烟囱在08 、3 ．7s 、6 ．0 5s 、 8 ．5 58 、1 08 、1 2 ．4 8 、1 3 ．5s 、1 5s 时间内过程中的倒 塌情况。 3 ．3 烟囱倒塌长度计算 2 烟囱由于周边环境限制，对倒塌长度有很高 的要求。利用A N S Y S 软件对本方案的模拟研究中， 可计算2 烟囱的倒塌长度。A N S Y S ／L s D Y N A 后 处理软件，空间坐标轴为 髫，，，， 。Y 轴为纵轴，分 析烟囱倒塌距离时，取髫，轴，后处理软件默认节点 单元起始坐标为 0 ，0 。建立模型时以茹轴的正方 向为北，z 轴正方向为西。1 。烟囱向西北方向倒塌， 2 。烟囱向东北方向倒塌。 图3 为2 4 烟囱顶部节点单元2 2 0 4 6 0 ，从图中可 知髫、彳的坐标为 2 9 ，2 9 。因此2 ’烟囱的倒塌距离 工 ／石2 z 2 4 1m 。图4 为1 ’烟囱顶部节单元 26 6 2 ，从图中可知髫、z 的坐标为 4 6 ，5 0 。因此l ’ 烟囱的倒塌距离L 以2 ， 6 7 ．9 4m 。因此认为 此控制爆破方案，烟囱倒塌长度满足要求，不会对周 边建筑物造成损坏。 万方数据 第2 8 卷第4 期 罗德丕，池恩安，詹振锵，等钢筋混凝土双烟囱交叉折叠倒塌爆破数值模拟 2 9 夕 - 翔- I 图22 烟囱模拟倒塌过程图 F i g ．2 S i m u l a t i o no ft w oc h i m n e y so fc o l l a p s ep r o c e s s 4 爆破效果与数值模拟过程比较 2 。烟囱上切口起爆后，上部6 2m 站立，1 ．2s 后 下部切口爆破，上部已形成向西南倾倒的趋势； 2 ．1s 时，上部倾斜约1o ，下部2 3m 已形成向东北 方向倾倒的趋势；3 ．3s 时，上下部倾斜均约2 0 ；3 ．7 s 时，上部6 2m 倾斜约1 5 。，下部倾斜约1 7 。；4 ．5s 时上部倾斜约2 5 。，下部基本倒地。6 ．0 58 时，上部 倾斜约6 5o ，此时1 。烟囱切口爆破，7 ．4s 时2 ’烟囱 基本全部倾倒着地。2 4 烟囱整个折叠倒塌过程历时 7s 。模拟计算的爆破后的效果表明烟囱上部6 2m 与下段2 3m 完全重叠在一起，2 ’烟囱实际倒塌距离 为4 3 ．5m 。7s 时l4 烟囱上部形成向西北倾倒的趋 势；8 ．5 5s 时，l4 烟囱倾斜约1 。；1 0s 时，爆破切口 上部分跟下部分重合，烟囱上部7 2m 开始对烟囱下 部1 0m 施加荷载；1 1 ．3s 时，烟囱上部倾斜约1 0o ； 1 2 ．4s 时，烟囱上部倾斜约2 0o ；1 3 ．5s 时，烟囱上 部倾斜约4 5 。；1 4 ．2s 时，烟囱上部倾斜约7 0 。；1 5s 时，烟囱上部基本全部倾倒着地。l 。烟囱倒塌过程 历时9s ，1 ’烟囱实际倒塌距离为7 1m 。烟囱的实际 倒塌过程和倒塌后状态见图5 。模拟计算爆破后的 效果表明，烟囱下部1 0m 会被烟囱上部坐塌。爆破 没有对周围建筑物、人员、设备造成危害，爆体破碎 彻底，因此可认为利用数值模拟的方法对方案进行 预演、校核有着重要意义。 _‘|J{J_l -．_J，11，，●●J --_li ．i 萤蓦 萤墓 翟蓦 韵弱 L 万方数据 3 0爆破2 0 1 1 年1 2 月 3 0 。2 5 l2 0 3 1 5 菖1 0 毛5 “0 - 5 5 、0 ；一5 8 1 0 喜一1 5 啐一2 0 “一2 5 3 0 图32 。烟囱2 2 0 4 6 0 单元位移时程曲线 F i g 3T h ed i s p l a c e m e n t t i m eCHIVeo f2 2 0 4 6 0u n i to f2 4c h i m n e y 5 结语 图41 。烟囱2 6 6 2 单元位移时程曲线 F i g ．4 T h ed i s p l a c e m e n t - t i m ecurveo f2 6 6 2u n i to f1 ’c h i m n e y a ／ - - 59s 时烟囱倒塌状态 b T 1 5s 时烟囱倒塌后状态 a 1 T h ec h i m n e yo f c o l l a p s es t a t e w i t h t h e t i m eo f 59s b T h ec h i m n e yo f c o l l a p s es t a t e w i t h t h e t i m ea f t e r l 5 日 图5 烟囱的倒塌过程和倒塌后状态 F i g 5T h ec o l l a p s ep r o c e s sa n dc o l l a p s es t a t eo ft h ec h i m n e y 1 根据钢筋混凝土烟囱折叠控制爆破的力学 原理，采用A N s Y S ／L s D Y N A 显式动力分析软件， 建立烟囱爆破拆除有限元模型，钢筋混凝土单元采 用分离式共节点模型，通过数值模拟验证了爆破技 术方案的可行性。从爆破效果看，2 座烟囱定向准 确，烟囱坍塌完全，未对周围的厂房、机械和电器设 备、生产工人造成危害。 2 通过数值模拟确定的2 。烟囱上、下爆破切 口延期时间为1 ．2s ，烟囱折叠效果明显，说明选取 1 ．2S 是合理的。 3 随着计算机性能的提高和大型计算软件的 发展，使得对结构倒塌的数值模拟成为可能。如能 采用合理的模型与材料参数，确定合理的边界条件， 并对单元细化，必将可以指导此类工程实践，使爆破 效果与预计设计目标更加吻合。 参考文献 R e f e r e n c e s [ 1 ] 杨勇，郭子雄基于A N S Y S 程序的钢筋混凝土梁非 线性数值模拟[ J ] 福建工程学院学报，2 0 0 5 ，3 6 6 2 8 - 6 3 2 [ 1 ] Y A N GY o n g ，G U OZ i x i o n gN u m e r i c a ls i m u l a t i o no fr e i n f o r c e dc o n c r e t eb e a m sb a s e dO i lA N S Y Sp r o g r a m [ J ] ． J o u m a lo fF u j i a nU n i v e r s i t yo fT e c h n o l o g y ．2 0 0 5 ，3 6 6 2 8 - 6 3 2 i nC h i n e s e 下转第3 6 页 h￡～S～u≈、nqH≈．N 万方数据 爆破2 0 1 1 年1 2 月 可知m 0 ，在此范围之内的某个值为最佳匹配值。 钢板与铜板的声阻抗比为1 ．1 4 ，钢板与铝板的 声阻抗比为3 ．1 6 ，钢板与橡胶板的声阻抗比为 3 9 ．2 7 ，因此，三者的匹配系数值分别为一0 ．0 7 、 一1 ．0 8 和一1 9 ．1 3 5 。虽然铝板的衰减系数小于橡胶 板，但试验中选取铝板为“软材料”时破坏程度最 小。因此设想，并不是衰减越大防护效果越好，对复 合结构中的“软材料”还有一定硬度要求，组合材料 应达到一个最佳匹配值。由此设想m 最佳取值范 围可能在一1 附近，但仍需进一步研究予以确定。 4结语 通过毁伤机理分析，运用有限元仿真技术及试 验研究相结合的手段对铝、铜、橡胶等材料作为“软 介质”的防护效果进行了研究，设想复合结构中应 具有一个最佳匹配值，证明了复合壳体结构具有良 好的抗爆性能。同时，为复合壳体结构的水下应用 提供了理论和技术上的支持。但是，由于试验相对 简单，仍需后续试验予以充分验证。下一步应通过 大量的数值模拟和试验研究，寻求具有最佳防护效 果的材料以及复合壳体的组合，还可以展开3 层复 合壳体结构水下爆炸损伤特性的研究，为水下船舶 结构设计提供理论和试验支持。 参考文献 R e f e r e n c e s [ 1 ] T E D E S C OJW ，H A Y E SJR ，L A N D I SDW ．D y n a m i c 陪 8 p o 舾eo fl a y e r e ds t r u c t u r e ss u b j e c tt ob l a s te f f e c t so fn o n n u c l e a rw e a p o n r y [ J ] ．C o m p u t e r ＆S t r u c t u r e ，1 9 8 7 ， 2 6 1 ／2 7 9 - 8 6 ． [ 2 ]周方毅，张可玉，张玮．复合壳体水下船舶抗爆性能 研究[ C ] ∥第二届全国舰艇抗冲击技术交流会，2 0 0 7 ． [ 2 ] Z H O UF a n g ．y i ，Z H A N GK e y u ，Z H A N GW e i ．R e s e a r c h f o rc a p a b i l i t yo fm u l t i p l eh u l l ss t r u c t u r eO i lr e s i s t i n ge x - p l o s i v ed a m a g e [ c ] ／／T h e2 n dP r o s e m i n a ro fS h i p s a ． g a i n s tS h o c kW a v e si nC h i n a ．2 0 0 7 ． i nC h i n e s e 姚熊亮，倪宝玉．水下爆炸对舰船结构毁伤研究进展 [ C ] ∥第三届全国舰艇抗冲击技术交流会，2 0 0 9 ． Y A OX i o n g h a n g ，N IB a o y u ．P e r s p e c t i v e so fd a m a g e e f f e c tt os h i pc o n s t r u c t i o ns u b j e c t e dt oU N D E X [ C ] ／／T h e 3 r de r o g e m i n a ro fS h i p sa g a i n s tS h o c kW a v e si nC h i n a ， 2 0 0 9 ． i nC h i n e s e 杨树涛，朱永凯，张阿漫．舰船舷侧防护结构抗冲击性能 研究[ c ] ∥第三届全国舰艇抗冲击技术交流会，2 0 0 9 ． Y A N GS h u t a o ，Z H UY o n g - k a i ，Z H A N GA - m a n ．R e s e a r c h f o rc a p a b i l i t yo fs h i pb r o a d s i d ed e f e n s i v es t r u c t u r es u h j e c t e dt o8 h o c kw a y e s [ C ] ∥T h e3 r dP r o s e m i n a ro fS h i p sa - g a i n s tS h o c kW a v e s i nC h i n a ，2 0 0 9 ． i nC h i n e s e 陈卫东，于诗源，王飞，等．多层板壳结构在水下接 触爆炸载荷作用下的试验研究[ J ] ．哈尔滨工程大学 学报，2 0 0 9 ，3 0 1 1 9 - 2 2 ． C H E NW e i - d o n g ，Y US h i - y u a n ，W A N GF e i ，e la 1 ．E x p e r i - m e n t a ls t u d yo fm u h i l a y e rs h e l ls t r u c t u r e si nu n d e r w a t e r c o n t a c te x p l o s i o n s [ J ] ．J o u r n a lo fH a r b i nE n g i n e e r i n gU - n i v e m i t y ，2 0 0 9 ，3 0 1 1 9 - 2 2 ． i nC h i n e s e 徐定海，盖京波，王善，等．防护模型在接触爆炸作 用下的破坏[ J ] ．爆炸与冲击，2 0 0 8 ，2 8 5 4 7 6 ．4 8 0 ． X UD i n g - h a i ，G A IJ i n g - b o ，W A N GS h a h ，e ta 1 ．D e f o r m a - t i o na n df a i l u r eo fl a y e r e dd e f e n s em o d e l ss u b j e c t e dt o c o n t a c te x p l o s i v el o a d [ J ] ．E x p l o s i o na n dS h o c kW a v e s ， 2 0 0 8 ，2 8 5 4 7 6 - 4 8 0 ． i nC h i n e s e 周方毅，张可玉，陈晓强，等．复合壳体抗水下爆炸作 用机理研究[ J ] ．爆破，2 0 0 6 ，2 3 2 9 1 - 9 2 ． Z H O UF a n g y l ，Z H A N GK e y u ，C H E NX i a o - q i a n ge ta 1 ． M e c h a n i s m0 fm u l t i p l eh u l l sr e s i s t i n gu n d e r w a t e re x p l o s i v e d a n l a g e [ J ] ．B l a s t i n g ，2 0 0 6 ，2 3 2 9 1 观． i nC h i n e s e 叶序双．爆炸作用基础[ M ] ．南京解放军理工大学，2 0 0 1 ． 库尔P ．水下爆炸[ M ] ．北京国防工业出版社，1 9 6 0 ． 上接第3 0 页 [ 2 ]万震雄．烟囱水塔类钢筋混凝土结构爆破拆除数值模 拟研究[ D ] ．北京北京理工大学，2 0 0 6 ． [ 2 ] W A N GZ h e n g x i o n g ．N u m e r i c a ls i m u l a t i o no fr e i n f o r c e d c o n c r e t ec h i m e ya n dw a t e r t o w e re x p l o s i v ed e m o l i t i o n [ D ] ．B e i j i n g B e i j i n gI n s t i t u t e o fT e c h n o l o g y ，2 0 0 6 ． i I l C h i n e s e [ 3 ] 池恩安．公路桥梁爆破拆除技术及数值模拟[ D ] ．武 汉武汉理工大学，2 0 1 1 ． [ 3 ] C H IE n - a n ．N u m e r i c a ls i m u l a t i o no fh i g h w a yb r i d g ee x - p l o s i v ed e m o l i t i o nt e c h n o l o g y [ D ] ．W u h a n W u b a nU n i v e r s i t yo fT e c h n o l o g y ，2 0 11 ． i nC l l i n e s e [ 4 ]魏兴，池恩安，邹锐．复杂环境下钢筋混凝土双烟 囱爆破拆除技术[ J ] ．爆破，2 0 0 8 ，2 5 2 5 7 - 6 0 ． [ 4 ] W E IX i n g ，C H IE n ．a n ，Z O UR u i ．B l a s t i n gd e m o l i t i o no f2 r e i n f o r c e dc o n c r e t ec h i m n e y si n c o m p l e xe n v i r o n m e n t [ J ] ．B L A S T I N G ，2 0 0 8 ，2 5 2 5 7 - 6 0 ． i nC h i n e s e 潘强，张继春，林凯．基于接触力学计算钢筋砼烟囱 定向爆破中土垫层厚度[ J ] ．爆破，2 0 1 1 ，2 8 2 2 8 - 3 0 ． P A NQ i 锄g ，Z H A N GJ i c h u n ，L I NK a i ．C a l c u l a t i o no fs o i l c u s h i o nt h i c k n e s si nr e i n f o r c e dc o n c r e t ec h i m n e y ’Sd i r e c ． t i o n a lb l a s t i n gb a s e do nc o n t a c tm e c h a n i c s [ J ] ．B L A S T - I N G ，2 0 1 l ，2 8 2 2 8 - 3 0 ． i nC h i n e 赞 夏卫国，曾政．海口华能电厂1 5 0 m 高钢筋混凝土烟 囱控制爆破拆除[ J ] ．爆破，2 0 1 l ，2 8 1 7 1 - 7 3 ． X I AW e i g u o ，Z E N GZ h e n g ．C o n - o H Me x p l o s i v ed e m o l i ． t i o no f1 5 0mr e i n f o r c e dc o n c r e t ec h i m n e yi nH a i k o uH u a ． n e n gp o w e rp l a n t [ J ] ．B L A S T I N G ，2 0 1 1 ，2 8 1 7 1 - 7 3 ． i nC h i n e s e 1 J 1 J 1 J 1 J 1 J 1 J 1 J 1 J 1 J 1 J 1 J 1 J 口 B H H 随 随 № № 口 口 随p 1 J 1 J 1 J 1 J 卵 朝 胡 乱 ． r L r L r L r L 万方数据