钢筋混凝土简支板梁桥爆破拆除数值模拟.pdf

第 3 1 卷 第 1 期 2 0 1 4年 3月 建筑科 学与工程 学报 J o u r n a l O f Ar c h i t e c t u r e a n d C i v i l E n g i n e e r i n g Vol I 31 M a r ． NO． 1 2 01 4 文章 编号 ： 1 6 7 3 — 2 0 4 9 ( 2 0 1 4 ) O l 一 0 0 5 6 - 0 7 0 钢筋混凝 土简支板梁桥爆破拆 除数值模 拟 黄平 明 ， 冯剑 平 ， 王 蒂 ， 朱 郑 ( 1 ．长安大学 公路学 院， 陕西 西安 7 1 0 0 6 4 ；2 ．山东高速股份有限公司京台分公司 枣庄养 护所 ， 山东 济南2 5 0 0 1 4 ) 摘要 ： 为 了实现对钢筋混凝土 简支板 梁桥爆 破拆 除过 程的有 限元数值 模拟， 采用动 力分析软 件 ANS YS ／ L S — D YNA 建立 了分 离式 共 节点模 型 ； 对 比分 析 了在板 梁 两侧 切 口与跨 中切 口采 用微 差 爆破 与 不采 用微 差爆破 2种 方案 下 的爆 破 效果 ， 并 对 2块 板 梁上 、 下缘 混凝 土的 受力过 程进 行 了分 析 。结 果表 明 ： 分 离式共 节 点模 型 可以体现 出钢 筋和 混凝 土 的力 学性 能差异 ； 采 用微 差爆破 可 以减 小爆堆 范 围 ， 使 得 更 多的钢 筋和 混凝 土单 元 失效被删 除 。 关 键词 ： 桥 梁工 程 ； 爆破 拆 除 ； 数 值模 拟 ； 分 离式共 节点模 型 ； 钢 筋混凝 土 简支板 梁桥 ； 微 差爆破 中图分 类号 ： T U7 4 6 ． 5 文献 标志 码 ： A Nu me r i c a l Si m u l a t i o n o f Bl a s t i ng De mo l i t i o n f o r RC Br i d g e s o f S i m pl y — s u pp o r t e d Pl a t e Gi r d e r HUANG P i n g — mi n g 。FE NG J i a n — p i n g ，W ANG Di ，ZHU Zh e n g ( 1 ．Sc ho ol of Hi g hwa y，Cha ng ’ a n Uni ve r s i t y，Xi ’ a n 71 0 06 4，Sha a nx i ，Ch i n a；2．Za o z hu a ng M a i nt e n a nc e De p a r t me n t ，B e i j i n g — Ta i b e i B r a n c h o f S h a n d o n g Hi — s p e e d Co mp a n y Li mi t e d ， J i n a n 2 5 0 0 1 4，S h a n d o n g，Ch i n a ) Ab s t r a c t ： I n o r de r t o r e a l i z e t he nu m e r i c a l s i m u l a t i on o f bl a s t i n g d e mol i t i o n f o r r e i nf or c e d c on c r e t e ( RC)br i d ge s o f s i mpl y — s up po r t e d pl a t e gi r d e r ，t he s e pa r a t e mod e l o f c o m mo n no d e s wa s e s t a b l i s h e d u s i n g t h e d y n a mi c a n a l y s i s s o f t wa r e ANS YS ／ I S — DYNA．Th e b l a s t i n g e f f e c t i v e n e s s be t we e n u s i n g mi l l i s e c o nd bl a s t i ng a n d wi t h out us i ng mi l l i s e c o n d bl a s t i n g a t t h e c ut s whi c h wa s on b o t h s i de s a nd t he pl a t e g i r d e r a t mi d — s p a n we r e c o m p a r a t i ve l y a na l y z e d， t he n t he l oa d i n g p r oc e s s e s o f t he c o nc r e t e b ot h o n t he u pp e r a n d l owe r e d ge we r e s t u d i e d ． The r e s ul t s s ho w t ha t t he s e D a r a t e mod e l o f c o mm o n no de s c a n r e f l e c t t h e di f f e r e nc e s o f t he me c h a ni c a 1 p r op e r t y b e t we e n t he s t e e l a n d c o nc r e t e；t h e s c o p e o f t he bl a s t i n g h e a p c a n be de c r e a s e d a nd m o r e s t e e l a nd c o nc r e t e e l e me nt s c a n b e i nv a l i d b y us i ng t he s c he me o f m i l l i s e c on d bl a s t i n g． Ke y wo r d s： br i d g e e ng i ne e r i ng； b l a s t i ng de mo l i t i o n； nu m e r i c a l s i m u l a t i o n； s e p a r a t e m o de l o f c o mmo n n o d e ；RC b r i d g e o f s i mp l y — s u p p o r t e d p l a t e g i r d e r ；mi l l i s e c o n d b l a s t i n g 引 日 随着 经 济 的 高速 发 展 和桥 梁 保有 量 的增 加 ， 越 来越多的桥梁 的通行能力已经不能满足车流量的增 加 ， 许 多桥梁 由于各种 原 因需要 拆 除或重 建 ， 爆破拆 除是拆除钢筋混凝土结构 的主要方法之一 ， 因其快 收稿 日期 ： 2 0 1 3 1 0 2 2 基金项 目： 国家 自然科学基金项 目( 5 1 2 7 8 0 6 4 ) ； 云南省交通运输厅科技 教育项 目( 云交科教[ 2 O 1 1 1 3 5 5 g - ) ； 云南省交通运输厅科技计划项 目( 云交科[ 2 0 1 1 ] 1 7 一 a ) 作者简介 ： 黄平 明( 1 9 6 5 一 ) ， 男 ， 湖北当阳人 ， 教授， 博士研究生导师 ， 工学博士 ， E — ma i l ： h p mi n g @v i p ． s i n a ． c o i n 。 学兔兔 w w w .x u e t u t u .c o m 第 1期 黄平 明 ， 等 ： 钢 筋 混凝 土 简支板 梁桥爆 破拆 除数 值模 拟 5 7 捷 、 高效的特点得到 了广泛 的应用。在 以往的爆破 拆除中， 人们多是通过经验来进行建 ( 构) 筑物 的爆 破拆除设计 ， 这样不能全面地反 映爆破拆除的力学 性能 。利用计算机模拟技术则可以对建( 构) 筑物的 爆破拆除过程进行数值模拟 ， 优化爆破设计方案 ， 预 演爆 破拆 除 过 程 。本 文 中笔 者 尝 试 采 用 分 离 式 共节点钢筋混凝土模型对钢筋混凝土简支板梁桥上 部结构的爆破拆除过程进行数值模拟 ， 为数值模拟 技术 在钢 筋混 凝土 桥梁 爆破 拆 除 中的应 用提供 理 论 依据和技术支持 。 1 数值计 算方法 在 AN S YS ／ L S — DYNA 中 ， 采 用 的是 一 种 显 式 方 法—— 改进 的 中心 差 分 时 间积 分 法 ] 。在 时 刻 的加 速度 定义 为 a = = = M F e x I — F【 n I ( 1 ) 式 中 ： a 为 加 速 度 矢 量 ； “为 施 加 的 外 部 力 矢 量 ； F ； 为内部力矢量 ； M 为总体质量矩阵 。 定义F == = ∑( IB d n+F “ ) =F ， 其 一 中 ， B为 应变 位移 矩 阵 ， F 为 沙 漏 阻力 ， F 为 接 触 力 ， 为单 元应 力 ， n 为单元 体积 。 速 度 y 和位移 U计 算式 为 V件△ ， ／ 2 一 Vl -A t／ 2 4 - a z At f ] l l 一 H f + Vf +△ 2 ／ 2 At +△ f／ 2 ( 2) At f +△ = =5 ( At ￡ — At + )J 式 中 ： △ 为时 间增 量 。 式( 1 ) ， ( 2 ) 不需要对刚度矩阵进行分解 ， 内部矢 量包 含所 有 的非 线 性 问题 ， 不需 要对 其 进 行 收 敛 检 查 。该方 法 在处理 接 触碰撞 、 爆 炸等 大 变形 、 大位 移 问题 中具 有 优势 。 2有限元模 拟 2 ． 1 材 料 的选择 本模型采用的是 ANS YS ／ L S — DYNA 中提供 的 * M AT —PLASTI C—KI NE MAT I C 材 料 模 型 7 ] 。 该材料模型是与应变率有关 的塑性随动强化模型， 具有 失效 应变 ， 是 一种 各 向同性 、 随 动强 化或各 向同 性 与随动 强 化 的混 合模 型 。在 冲击 或爆 炸荷 载作 用 下 ， 结构材料会 产生快速 变形 ， 应 变率也会 明显提 高。应变率采用 C o wp e r — S y mo n d s模型来考 虑， 采 用与应变率有关的因数来表示屈服应力 ， 即 一 [ 1 +( ) 寺 ] ( + ￡ ) ( 3 ) 式 中： ￡ 为应变率 ； C， P均为 C o wp e r — S y mo n d s应变 率参数 ； [ 1 +( ) 古 ] 为材料 放大因子； 为初始屈 U 服应力 ； 为可调强化参数 ， 一0为塑性随动强化模 型， 一1为等 向强化 模 型； E 为 塑 性 强化 模量 ， F F E 一 1 2 ~ t a n l & ，E， E 分 别 为 弹性 模 量 和 切 线 模 量 ； rt 0 e 为等效塑性应变， e pff — I( ／ , 。 p ￡ 3 ) d t ， ￡ 8 为塑性 J 0 0 应变率 ， ￡ 8 一。 _j —e ； ， e 。 j 为总 应 变 率， e ； 为 弹 性应 变 率 。 模型中钢筋和混凝土材料均为塑性随动硬化材 料， 钢筋和混凝土材料的力学性能参数如表 1 所示 。 表 1 材料的力学性能参数 Ta b ．1 M e c ha ni c a l Be h a vi o r Pa r ame t e r s o f M a t e r i al s 密度／ 弹性模 抗拉强 抗压 强 材料类型 泊松 比 ( k gm ) 量／ G P a 度／ MP a 度／ MP a 钢筋 7 8 0 0 2 0 0 0 ． 2 7 0 0 2 8 0 2 8 O 混凝土 3 2 0 0 3 0 0 ． 1 6 6 7 6 2 O 2 ． 2 模 型 的建立 2 ． 2 ． 1 基 本假 设 ( 1 ) 采用分离式共节点模型模拟结构 的倒塌过 程， 并且不考 虑钢筋 和混凝 土 材料 之间 的粘结 和 滑移口 ] 。 ( 2 ) 不 考虑爆 破 切 口破 坏 形 成 的过 程 和 炸 药 爆 炸对整体结构的作用 ， 采用直接删除爆破切 口内混 凝 土单 元 的方法 来模 拟爆 破切 口的形成 。 ( 3 ) 爆 破切 口是 平滑 的 。 ( 4 ) 不 考虑 钢筋 混 凝 土 结 构 倒 塌过 程 中与 地 面 的相互 作用 。 2 ． 2 ． 2 模 型 概 况 该模型为单跨钢筋混凝土简支板梁桥 ， 跨径为 1 0 m， 桥面宽度为 9 m。对于该桥梁结构， 采用有限 元 软件 ANS YS ／ L S — DYAN建 立分 离式 共 节点 钢 筋 混凝土有 限元模 型 ， 在 L S — DY AN环境下 求解 ， 在 L S — P R E P O S T 中 进 行 后 处 理 分 析 。 钢 筋 采 用 B e a m1 6 1单元 ， 混凝 土 采用 S o l i d 1 6 4单元 。 B e a ml 6 1单元 用 3个 节点 定 义 ， 即节 点 ， J均 确定梁的轴向， 节点 K确定横截面的主轴方位。 S o l i d 1 6 4单元是三维实体单元 ， 具有 8个节点 ， 每个节点有 9个 自由度， 只有位 移是具有实际意义 上的物理 自由度 。 由于板梁为空心板结构 ， 采用扫掠生成体 网格 的方式 对 板梁结 构 进 行 网 格划 分 ， 网格 划 分 尺 寸 为 0 ． 1 m， 有限元模型如图 1 所示。 学兔兔 w w w .x u e t u t u .c o m 5 8 建 筑科 学与 工程 学报 2 0 1 4年 图 1 有 限 兀模 型 Fi g ． 1 Fi n i t e El e m e n t M o de l 由于 板 梁 在 下 落 的过 程 中 既包 含 与 桥 台 的接 触 ， 又包 含与地 面 的 接 触 ， 接 触关 系 比较 复 杂 ， 采 用 L S — D YNA 提供 的 AUT O— S I NGL E — S URF AC E 的 接触 方式 来模 拟 ， 该 接 触方式 可 以 自动 搜索 接触 面 、 判断接触面并可以处理侵蚀 、 断裂等复杂的边界变 化情 况 ， 材料 的摩 擦 因数 设 定 为 0 ． 6 。该 接 触 方 式 自动确定壳单元接触表面方 向， 自动检查壳单元 的 每一面 ， 通常限制搜索深度 。 2 ． 3爆破 模拟 方案 钢 筋 和混 凝 土 材料 的破 坏 准则 非 常 复 杂 ， 涉及 到材料的本构关系和各种失效准则等因素， 本 文中 钢 筋 和混凝 土材 料 的 失效 均 由应 变 控 制 ， 其 失 效应 变 分别 为 5 ． 0 ， 0 ． 4 8 。在 有 限元 数 值 模 拟 中 ， 材 料 的失 效是 靠删 除单 元来 实现 的 。根据 微差爆 破 和 重 力作 用 的原理 ， 利 用 关 键 字 *MAT— AD D— E R O— S I ON来 删除 混凝 土单 元 网格 并形 成 爆 破切 口。板 梁在重力作用下发生失稳、 破坏、 倒塌落地， 通过关 键 字 *L OA D — B ( ) DY — Y( Y方 向) 来 实 现 。 本 文 中拟 对钢筋 混凝 土简 支板 梁进 行微 差爆 破 拆除 ， 保留两侧桥台， 采用 2种方案进行数值模 拟。 方案 1 ： 爆 破 拆 除 板 梁 两侧 各 1 ． 0 m 及跨 中 1 ． 0 m 的板 梁 结 构 ， 采 取 从 板 梁 一 侧 到 板 梁 另 一 侧 微 差 0 ． 5 s 的爆破 方案 ， 切 口形 成 的时 间分别 为 0 ． 5 ， 1 ． 0 ， 1 ． 5 S 。方 案 2 ： 爆破 拆 除 板梁 两 侧 各 1 ． 0 m 及 跨 中 1 ． 0 m 的板梁 结 构 ， 采 取 板 梁两 侧 切 口同 时形 成 并 与 中 间切 口微 差 0 ． 5 S的爆 破 方 案 ， 切 口形成 时 问 分别 为 0 ． 5 ， 1 ． 0 S 。 3数值 结果分 析 3 ． 1倒塌 过程 模拟 结果 比较 3 ． 1 ． 1 方案 1数值 模拟 效 果 方 案 1的爆 破拆 除倒 塌效 果如 图 2所示 。跨 中 切 口两侧 第 1个 切 口、 第 2个切 口附近 板 梁 单元 的 方向位移一 时间曲线分别如图 3 ( a ) ， ( b ) 所示 。 从 图 2 ， 3可 以看 出 ： ( 1 ) 当 t = = = 0 ． 6 S时， 第 1个 爆 破切 口刚形 成 之 后 。 ( 2 ) 当 一1 ． 1 S时 ， 跨 中爆破 切 口刚形 成之 后 。 ( 3 ) 当 ￡ 一 1 ． 6 S时 ， 第 2个 爆 破 切 口形 成 ， 第 1 个 切 口处 的钢筋 由于 与桥 台强 烈 的摩擦 ， 出 现 了混 凝 土 和钢筋 单元 失 效 破碎 被 删 除 的情 况 ， 跨 中切 口 处 的钢 筋也 开始 出 现变 形 ， 并 拉 扯 着 第 2切 口处 的 板 梁下 落 。 ( 4 ) 当 t 一3 ． 9 S 时 ， 两侧 切 口处 的钢 筋都 在与 桥 台 的摩 擦作 用 中失 效被 删 除 ， 第 1个 切 口处 板 梁 的 下落速 度 明显快 于 第 2个 切 口处 的板 梁 ， 跨 中切 口 处 的钢 筋 在 两 侧 混 凝 土 板 梁 的拉 扯 作 用 下 缓 慢 变 形 ， 并 拉着 2块 板梁 互相靠 拢 。 ( 5 ) 当 t 一5 ． 2 S时 ， 第 1 个 切 口处 的钢筋 开 始着 地， 第 2 个切 口处的钢筋仍在与桥台接触摩擦 ， 大量 钢筋 失效 被删 除 。 ( 6 ) 当 t 一6 ． 0 S时， 第 1个切 口处的板梁落地 ， 先落地的板梁通过跨中裸露出来 的钢筋拉扯着第 2 个切 口处的板梁下落 ， 同时第 1块板梁向跨 中靠拢 的速度加快[ 图 3 ( a ) ] 。 ( 7 ) 当 t 一7 ． 0 S时 ， 第 2个 切 口处 的板梁 落 地 ， 在跨 中裸 露钢 筋 的拉 力 作 用 下 ， 第 1块 板 梁靠 近桥 台的一边 翘起 ， 并逐 渐靠 近第 2 块 板梁 ， 第 2 块 板 梁 同样 在裸 露钢 筋 的作 用 下 ， 出现 了远 离第 1块 板 梁 的趋势[ 图 3 ( b ) ] ， 但是速度不及第 1块板梁大。 ( 8 ) 当 ￡ 一9 ． 0 S时， 2块板梁基本靠拢 ， 跨 中裸 露 钢筋 严重 失效 。 ( 9 ) 当 ￡ = = = 1 0 ． 0 S 时， 2块板梁落地靠拢完毕， 两 者 移动 的位 移之 和基 本上 等于跨 中切 口长 度 1 ． 0 m ( 图 3 ) 。 3 ． 1 ． 2 方 案 2数 值 模 拟 效 果 方 案 2的爆 破拆 除倒 塌效 果如 图 4所示 。跨 中 切 口两侧第 1 个切 口、 第 2个切 口附近板梁单元 的 方 向位移一 时 间 曲线 分别 如 图 5 ( a ) ， ( b ) 所 示 。 从 图 4 ， 5可 以看 出 ： ( 1 ) 当 一0 ． 6 S 时 ， 板梁两 侧 的爆破 切 口刚形成 之后 。 ( 2 ) 当 t 一1 ． 1 s 时 ， 跨 中爆破切 口刚形成之后。 ( 3 ) 当 t 一 1 ． 6 S 时 ， 两侧 切 口处 的钢筋 开始 与桥 台产 生摩 擦作 用 ， 出 现 了钢 筋 和 混 凝 土 单元 失 效 破 碎 被删 除的情 况 ， 跨 中切 口处 钢 筋也 开始 出现 变形 。 ( 4 ) 当t 一 3 ． 9 S 时 ， 伴 随 着 板 梁 的 不 断 下 落 ， 两 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 第 1期 黄 平明 ， 等 ： 钢 筋 混凝 土 简支板 梁桥 爆破 拆 除数值 模 拟 6 1 时 间／ s ( a ) 第1 个切 口处 板梁 上 缘混 凝 土 ， ＼ 。 时 间／ s ( b ) 第 1 个切 13处 板梁 下 缘混 凝 土 时间／ s ( c ) 第2 个 切 口处 板 梁上 缘 混 凝土 时 间／ s ( d ) 第2 个 切 口处 板 梁 下缘 混凝 土 图 7 方案 1混凝土 的应 力一 时间 曲线 Fi g ． 7 S t r e s s - t i me Cur v e s o f Co nc r e t e i n Sc he me One 缘大 ， 直 到板 梁 落地 ( ￡ 一6 ． 0 s ) 之 前 ， 混 凝 土几 乎 都 是受到拉应力作用， 至板梁落地 时达到 了一个拉应 力 极 值 ； 板 梁 落地 时混 凝土 达到 了一 个压 应力 极值 ， 之后混凝土不再受任何作用力 。从 图 8 ( b ) 可以看 出 ： 在 t 一2 ． O s之 前 ， 下 缘 混 凝 土 几 乎 不 受 作 用 力 ； 当 2 ． 0 s <￡ 1 ． 0 s ) ， 上缘混凝土受拉 ； 直到 t 一5 ． 0 S时 ， 第 1块 山 皇 -叵 一 对 山 羔 匠 时 间／ s ( a ) 第 1 个 切 口处 板梁 上 缘混 凝 土 时 间／ s ( b ) 第 1 个 切 口处 板梁 下 缘混 凝 土 R -叵 时 间／ s ( c ) 第2 个切 口处板 梁 上 缘混 凝 土 时 『司／ s ( d ) 第2 个切 口处板 粱 下 缘混 凝 土 图 8 方案 2混凝土 的应 力- 时间 曲线 Fi g ． 8 St r e s s - t i me Cu r v e s o f Con c r e t e i n S c h e me Two 板梁切 口处的钢筋落地， 通过跨 中裸露钢筋的作用 ， 混凝土受到一个压应力极值 ， 而后 又受拉 ； 直到第 2 块板 梁 落地 ( t 一 7 ． O s ) ， 混 凝 土 受 到 一 个 拉 应 力 极 值 ， 后 又缓 慢 减小 ， 直 到 整个 板 梁 完 全 落 地。从 图 8 ( c ) 可知 ： 在 t 一5 ． 5 S之前 ， 上缘 混 凝 土一 直 受 压 ； 当 t 一5 ． 5 S时 ， 混 凝 土受 到轻微 的拉 应力 ， 随后 处于不断变化 的拉应力与压应力交替状态 ； 直到t 一 7 ． 2 S时 ， 板梁 开 始落 地 ， 混凝 土受 到 一 个 拉应 力 极 值 ； 之 后 拉应 力 减小 ， 混凝 土 受压 直 到板 梁完 全 落地 。 从图 7 ( d ) 可以看出： 下缘混凝土受到微小 的压 学兔兔 w w w .x u e t u t u .c o m 6 2 建 筑科 学与 工程 学报 2 0 1 4生 应 力作 用 ； 第 1 块 板 梁落 地 时 ( t = ： = 6 ． 3 s ) ， 压应 力 有 稍 许增 加 ； 第 2块 板 梁落 地 时 ( t 一7 ． 0 s ) ， 混凝 土 受 到一个拉应力极值 ( 6 ． 5 MP a ) ， 此时混凝土达到 了 抗拉失效强度 ， 之后不再受力 ， 直到板梁完全落地。 从 图 8 ( d ) 可 以看 出 ： 在板梁 落地 之前 ( ￡ 7 ． 2 s ) ， 混凝 土 又 处 于幅 度变化比较大的拉应力与压应力交替状态 ； 到板梁 完全 落地 之 前 的近 2 S ， 混 凝 土 受 到 较 小 的 拉 应 力 作 用 。 4 结语 ( 1 ) 基 于动 力 分析 软 件 AN S YS ／ L S - DYNA， 采 用分离式共节点的钢筋混凝土模型可以模拟板梁结 构从下落到完全落地的过程 ， 从而实现 了桥梁结构 爆破 拆 除塌 落过 程 的连续仿 真模 拟 。虽然 没有 结合 工程 实 际或试 验 ， 但 可 以为今 后 类 似 桥 梁 的爆破 拆 除数值模拟提供借鉴 。 ( 2 ) 通过 对 比板 梁 两侧 微差 起 爆 和 同时 起爆 可 以得 出： 微差爆破使得跨中切 口两侧 的板梁在落地 后 完全 靠拢 ， 节 省 了场地 ； 微差 爆破 使得 2块 板梁 的 下 缘混 凝 土在其 落 地 后 不再 参 与受 力 ； 微 差 爆 破 使 得 更 多 的钢筋 和混凝 土单 元 由于下 落速 度 的反 复变 化 失效 被 删除 。 ( 3 ) 通 过对 钢筋 和混 凝 土 的应 力 随 时 间不 断 变 化 的分 析 可 以得 出 ， 利 用 分 离 式 共 节点 模 型能 够 较 好 地体 现 出这 2种 不 同材 料 在力学 性 能上 的差异 。 参考 文献 ： Re f e r e nc e s ： [1 ] 胡 文军 ， 卢艳华 ． 钢结 构建筑 爆破拆 除计 算机模 拟初 步研究[ J ] ． 工程爆破 ， 2 0 0 2 ， 8 ( 4 ) ： 1 1 — 1 5 ， 1 O ． HU W e n — j u n ，LU Ya n - h u a ．P r e l i mi n a r y S t u d y o n Compu t e r Si mul a t i o n o f Bl a s t i ng De mol i t i on of St e e l C o n s t r u c t i o n [ J ] ． E n g i n e e r i n g B l a s t i n g ， 2 0 0 2 ， 8( 4 ) ： 1卜 1 5， 1 0． [2 ] 余 业清 ， 钟 冬望 ， 汪洋． 数值模拟在爆破拆除 中的应 用E J ] ． 爆破 ， 2 0 0 6 ， 2 3 ( 2 ) ： 2 2 — 2 5 ． YU Ye - qi n g， ZHONG Don g— wa ng， W ANG Ya ng ．Ap — pl i c a t i on of Nu m e r i c a l Si mu l a t i on i n Expl o s i ve De mo — l i t i o n [ J ] ． B l a s t i n g ， 2 0 0 6 ， 2 3 ( 2 ) ： 2 2 — 2 5 ． [3] 贾金河 ， 于亚伦． 建 筑物拆 除爆破 数值 模拟 方法 的研 究E J ] ． 工程爆破 ， 1 9 9 9 ， 5 ( 1 ) ： 4 7 — 5 l l J I A J i n — h e 。 YU Ya — l u n ． S t u d y o n Nu me r i c a l S i mu l a — t i o n f o r D e mo l i t i o n B l a s t i n g o f B u i l d i n g [ J ] ． E n g i n e e r — i ng Bl a s t i n g， 1 99 9， 5( 1)： 47 — 51_ [4] [ 5] [6 ] [7 ] [8] [ 9] [ 1 O ] [ 1 1 ] [ 1 2 ] 谢春 明， 杨 军 ， 薛 里 ． 高耸筒形结构爆破拆 除的数 值模 拟[ J ] ． 爆 炸与冲击 ， 2 0 1 2 ， 3 2 ( 1 ) ： 7 3 — 7 8 ． XI E Ch u n - mi n g， YANG J u n， XUE Li ． Nu me r i c a l S i m- To we r i n g Tu b b y S t r u c t u r e B l a s t i n g D e mo l i t i o n [ J ] ． E x p l o s i o n a n d S h o c k W a v e s ， 2 01 2 ， 3 2 ( 1 ) ： 7 3 — 7 8 ． 宋天培 ， 谢春明 ， 杨 军 ． 公路双 曲拱桥结构爆 破拆除 数值模拟 [ J ] ． 爆破 ， 2 0 1 2 ， 2 9 ( 4 ) ： 9 0 — 9 4 ． S 0NG Ti a n - p e i ， XI E C h u n - mi n g， YANG J u n ． Nu me r — i c a l Si mu l a t i on of Exp l o s i ve De mo l i t i o n of H i gh wa y Do u b l e C u r v e d Ar c h B r i d g e[ J ] ．B l a s t i n g ， 2 0 1 2 ， 2 9 ( 4 ) ： 9 0 — 9 4 ． Li v e r mor e Sof t wa r e Te c h no l o gy Cor p or at i o n．LS— DY— N A Us e r ’ S Ma n u a l ： Ve r s i o n 9 7 0 [ M] ． L i v e r mo r e ： L i v — e r mo r e S of t war e Te c hn ol o gy， 2 00 3． 石少卿 ， 康建功 ， 汪 敏 ， 等． A NS YS ／ I S — DY NA 在爆 炸与冲击领域 内 的工 程应 用[ M] ． 北 京 ： 中国建筑 工 业 出 版 社 ， 2 0 1 1 ． S HI S h a o — q i n g ， KANG J i a n — g o n g， W ANG M i n， e t a 1 ． Th e E n g i n e e r i n g Ap p l i c a t i o n o f ANS YS ／ L S — DYNA i n t h e F i l e d o f E x p l o s i o n a n d I mp a c t [ M] ． B e i j i n g ： Ch i n a Ar c h i t e c t u r e & B u i l d i n g Pr e s s ， 2 0 1 1 ． 杨军 ， 杨 国梁 ， 张光雄． 建筑 结构 爆破 拆除数值模 拟 [ M] ． 北京 ： 科学 出版社 ， 2 0 1 2 ． YANG J u n，YANG Gu o — l i a n g。ZHANG Gu a n g — xi on g．Nu me r i c a l Si mu l a t i on of Bl a s t i ng Demol i t i o n f o r A r c h i t e c t u r a l S t r u c t u r e[ M ] ．B e i j i n g ：S c i e n c e Pr e s s。 20 12 ． 言志信 ， 于换小 ， 朱 辉． 框架结构爆破拆除 过程研究 [ J ] ． 爆破 ， 2 0 1 2 ， 2 9 ( 4 ) ： 1 9 — 2 2 ， 2 7 ． YAN Zh i — x i n， YU Hua n — x i a o， ZHU Hui ．Re s e a r c h o n E x p l o s i v e D e mo l i t i o n P r o c e s s o f F r a me S t r u c t u r e [ J ] ． Bl a s t i n g， 20 12， 2 9( 4)： 1 9 — 22， 2 7． 池 恩安 ． 公路 桥梁 组合拆 除 爆破及 数 值模 拟[ D] ． 武 汉 ： 武汉理工大学 ， 2 0 1 1 ． CHI En — a n．Nu me r i c a l Si mul a t i o n an d Appl i c a t i on of Hi g h wa y B r i d g e De mo l i t i o n B l a s t i n g [ J ] ． Wu h a n ： Wu — ha n Uni v e r s i t y of Te c hn ol o gy， 2 0 1 1 ． 杨 国梁． 钢筋混凝 土建筑 结构爆破 拆除数值模拟研究 [ D ] ． 北京 ： 北京理工大学 ， 2 0 0 9 ． YANG Gu o — l i a ng． Nume r i c al Si m u l a t i o n on Re i n— f o r c e d Co n c r e t e S t r u c t u r e B l a s t i n g De mo l i t i o n r D] ． B e i j i n g： B e i j i n g Un i v e r s i t y o f Te c h n o l o g y ， 2 0 0 9 ． 杨忠华． 钢 筋混 凝土桥梁结构爆破拆除数值模拟研 究 [ D] ． 北 京 ： 北京理工大学 ， 2 0 1 0 ． YANG Zho ng — hua ．Nu me r i c a l Si mul at i o n o n Bl a s t i n g De mol i t i o n o f Br i d ge Con s t r u c t i on of Re i nf o r c e d Con— c r e t e S t r u c t u r e[ D] ．B e ij i n g ： B e i j i n g Un i v e r s i t y o f Te c hn ol ogy， 201 2 ． 学兔兔 w w w .x u e t u t u .c o m