公路双曲拱桥结构爆破拆除数值模拟.pdf
第2 9 卷第4 期 爆破 V o l 2 9 N o 4 2 0 1 2 年1 2 月 B L A S T Ⅱ呵G D e c .2 0 1 2 d o i l O .3 9 6 3 /j .i s s n .1 0 0 1 4 8 7 X .2 0 1 2 .0 4 .0 2 2 公路双曲拱桥结构爆破拆除数值模拟 宋天培1 ,谢春明1 ,杨军2 . 1 .贵州新联爆破工程有限公司,贵阳5 5 0 0 0 2 ;2 .北京理工大学爆炸科学与技术国家重点实验室,北京l O 0 0 8 1 摘要为了确保双曲拱桥的爆破拆除能够按照设计方案顺利进行,采用共节点分离式钢筋混凝土模型,对 典型公路双曲拱桥结构爆破拆除过程进行了数值模拟,方案选择桥拱顶部预切缝0 .5m 以及副拱立柱爆破范 围1 .6m 。模拟结果表明对混凝土和钢筋单元的受力过程进行分析,共节点分离式模型可以体现混凝土和钢 筋材料的力学性能差异;预切缝采用o .5m 足以保证桥身相互叠加彻底倒塌,副拱立柱爆破范围采用1 .6m 可 使得桥体解体比较完全,并且可以减小整体倒塌时间和爆堆高度。研究认为结构倒塌过程的数值模拟。将成为 研究结构爆破拆除力学过程的重要手段并辅助指导结构爆破拆除设计,具有重要的工程实用价值。 关键词钢筋混凝土双曲拱桥;数值模拟;爆破拆除;共节点分离式模型 中圈分类号T u 7 4 6 .5文献标识码 A 文章编号1 0 0 1 4 8 7 x 2 0 1 2 0 4 0 0 9 0 0 5 N Ⅲn e r i c a lS i m l l l a t i o n H i g h w a yD o u b l e o fE x p l o s i V eD e m 础t i o no f C u r V e dA r d lB r i d g e S 叫崛丁缸凡- p e i l ,脚饥Ⅱn 一施n 9 1 ,y A 临血n 2 1 .G h i z h o uX i n l i a nB l a s tE n g i n e e r i n gU m i t e dG b r p ,G h i y a n g5 5 0 0 0 2 ,C h i n a ; 2 .S ta _ t eK e yL a b o r a t o r yo fE x p l o s i o nS c i e n e ea n dT e c h n o l o g y , B e i j i n gI n s t i t u t eo fT e c h n o l o g y ,B e i j i n g1 0 0 0 8 l ,C h i n a A b s t r a c t T oe n s u r et 壬l ed e m 0 1 i t i o nb l 觞t i l l go fd o u b l ec u r v e da r c hb r i d g es m o o t h l y ,an u m e r i c a ls i m u l a t i o no n t y p i c a lh i g h w a yd o u b l ec u n r e da r c hb r i d g es t n I c t u r eb l 鹊t i n gd e m o l i t i o np r o c e s sw 鹪c a Ⅱi e do u tb 鹤e do nc o m n l o n n o d es e p a r 砒e 阳i n f o r c e dc o n c r e t em d d e lw i t h0 .5mw i d ep r e - c u t t i n gs e 砌蚰dv i c e a r c hc o l u m nb l 鹪t i n gm n g e1 .6 m .R e s u h ss h o w 出a tt h ec o 玎l m o nn o d es e p a r a t em o d e lr 趣e c t m e c h 蛆i c a l 滞m e 8d 证I 玳n c eo fc o n c r e t e 柚dr e i n - f o r c e m e n tb 嬲e do nt l l e i rl o a d i n gp r o c e s s .T h e0 .5mw i d ep r e - c u t t i n gs e 啪i ss u f f i c i e n tt oe 璐u r et l l es u p e q 0 s i t i o n 蚰dc o u 叩8 ef u l l y .弧es c h e m eo fv i c e - a r c hb l a s t i n g 础g e0 f1 .6 mc 粕m a I c et I l e d g ed i 8 i n t e g m t e dm u c hm o r e , w h i o ha l s od e c r e 酗e dt l l ew h o l ec o l l 印s et i n l e 髓dm eh e i g h to fm u c kp i l e .T h e8 i m u l a t i o no f 出es t n l c t u 聆c o l l 印s e 访l lb e 锄i m p o r t a n tm e a n st os t u d yt l l em e c h a I l i c a lp r o c e s so fs t m c t u r eb l 髓t i n gd e m o l i t i ∞. K e yw o r d s r e i n f o r c e dc o n c r e t ed o u b l ec u r v e da r c hM d g e ;删m e r i c a ls i m u l a t i o n ;e x p h i v ed e m o l i t i o n ;c o m - m o nn o d es e Da r a _ t em o d e l 随着经济高速发展,越来越多的桥梁满足不了 车流量的要求,需要拆除在原址上建设新桥。钢筋 收稿日期2 0 1 2 1 1 一0 1 作者简介宋天培 1 9 5 8 一 ,男,贵州凯里人,主要从事工程爆破及 安全管理工作, E m a i l 5 1 1 7 9 8 5 6 0 q q .c o m 。 通讯作者谢春明 1 9 8 3 一 ,男,博士,从事爆破拆除数值模拟方向 研究, E m a i l x i e c h m i n g 殍I I a i l .c o m 。 混凝土结构是公路双曲拱桥的主要材料结构形式, 而爆破是对大型建构筑物进行拆除的主要办法。在 拆除爆破中,行之有效的预处理技术措施往往能达 到令人满意的爆破效果,极大地提高施工中的效费 比,因此很多学者从爆破方法、装药参数、爆破位置 等角度对桥梁结构的爆破拆除进行了研究,如文 献[ 1 引。尝试采用共节点分离式钢筋混凝土模型,对 万方数据 第2 9 卷第4 期 宋天培,谢春明,杨军公路双曲拱桥结构爆破拆除数值模拟 9 1 典型的公路双曲拱桥结构爆破拆除进行数值模拟研 究,具体分析设计方案对结构破坏倒塌的影响,为数 值模拟方法在爆破拆除的实际应用提供理论依据和 技术基础。 5 血。 缸。缸 。 方程不形成总体刚度矩阵,内部矢量包含所有 非线性问题,不需要进行收敛检查。在处理接触碰 撞、爆炸等大位移大变形问题具有优势。 1 数值计算方法2 模型建立及爆破方案选择 采用有限元软件D Y N A q D 进行数值计算。 程序中使用中心差分时间积分法∞J 。T 时刻的加速 度定义如下 { 口。} [ M ] - 1 [ F 芦I ] 一[ ∥“] 1 这里{ F } 为施加的外部体矢量;{ F I 埘} 为内部 力矢量。定义F 妇 ∑ I 艿7 盯。d 力 F k r 一‘, 白 式中,一为沙漏阻力,P ““为接触力。 速度和位移由下式计算 { K M } { K A 以} { 口;} △£。 2 这里{ M 。 缸} { M 。} { E 蚍} m 。 ∽;△f 。 们 2 .1 桥体概况 该公路桥在使用多年后主桥拱与复拱出现了多 处裂缝,复拱与盖梁连接处砼铰部位严重变形裂缝, 重载汽车已无法安全通过,被定为险桥,现决定爆破 拆除。 此桥为钢筋混凝土1 7 孔双曲拱桥,每拱跨度 3 1 .6m ,全桥共长5 5 4m 。桥梁宽度为1 5 .6m ,桥面 距离河床6 .3m 。主桥拱由1 0 排宽0 .5m 、厚1 .0m 的弧形拱肋支撑,拱肋间距为lm ,净跨3 0m ,净矢 高5m 。盖梁长1 5 .6m ,宽、高都为0 .4m ,由1 0 根 立柱支撑,立柱间距3 .9m ,截面尺寸为4 0c m 5 0c m 。桥体结构参见图1 。 1 6 1 7 西 图l 双曲拱桥结构示意图 单位m F i g .1 1 kd o 曲l ec u n ,e da r c h 硒d g e8 c h 咖觚cd i 8 班吼 u I l i t m 2 .2 爆破方案 采用半秒雷管分段延时起爆,模拟中具体起爆 位置和延迟时间如图2 所示,②一③一为起爆延 迟时间,图中①为桥面预切缝,可通过D Y N A 提供 的关键字M A T A D D E R O S I O N 在0 .0 1s 时从模型 中删除,从而实现实际工程中的预拆除处理。 ⋯l ’2 ‘3 ‘4 。5 6 ‘7 。8 ‘.9 ’ ②③⋯代表雷管分段 图2 起爆位置和延迟时间示意图 F i g .2 S k e t c ho ft h ed e t o n 砒i o np o s i t i o na n dd e l a yt i m e 模拟中采用设计方案如图3 所示。为了使桥身 爆破塌落后达到较好的破碎效果,在桥拱顶部位置 处,从桥面向下开凿一条宽5 0c m 的切缝,剔除已暴 露的钢筋,爆破后使桥体断开,以达到充分破碎的效 果。为了确保桥身彻底坍塌,副拱内的立柱需要爆破 处理,以增加爆破缺口高度,立柱爆破范围为1 .6m 。 2 .3 有限元模型的建立 模拟采用有限元软件A N S Y s /L s D Y N A 建立 有限元模型,运用大型显示动力分析软件L S D Y N A 进行求解。在本模拟中,采用分离式共节点模 型建立部分钢筋混凝土拱桥,充分体现钢筋在结构 倒塌过程中的拉应力作用。拱桥共为1 7 跨,取其2 跨进行建模分析,考虑到结构的对称性,在第9 跨中 间设置对称边界,沿桥宽度方向取1 /2 模型进行计 算。所建模型大小与实际结构完全相同,考虑到桥 梁部件的不规则性,采用四面体单元进行网格划分, 网格大小控制在o .4m 左右。具体有限元模型与钢 筋分布如图4 所示。 万方数据 爆破2 0 1 2 年1 2 月 图3 爆破拆除方案 单位m F i g .3 n eb l 鹪t i n gd e m o l i t i o np l 粕 u n i t m a 桥梁模型 2 .4 单元类型、材料及接触方式 模型中,混凝土与钢筋均采用随动硬化双线性 弹塑性材料,混凝土采用S O u D l 6 4 单元,钢筋采用 B E A M l 6 1 单元J 。钢筋混凝土材料的破坏非常复 杂,涉及材料本构关系、各种失效准则等因素,文中 钢筋和混凝土材料均由应变控制材料失效,失效应 变分别为5 %和0 .4 8 %。钢筋与混凝土的物理力学 参数见表1 。 图4 有限元模型 单位m F i g .4 ,n I ef i n i t ee l e m e n tm o d e l 岫i t m 表1 材料的物理力学参数 T l b k1 n 唧d c 岫彻旺m ∞蛔n Ⅻ舯涨t e 璐0 f m 曲蛐‘ 銎暴 0 髻, 弹繁影泊松比抗絮影抗譬影名称 k g m 。 P a ⋯⋯⋯ P aP a 钢筋 7 8 5 02 l O E 9O .2 73 5 0 E 63 1 0 E 6 混凝土 3 2 0 03 l E 9O .2 11 0 E 63 0 E 6 由于倒塌过程中接触非常复杂,因此本模拟中 选用了I s D Y N A 提供的A U T O S I N G L E S u R F A C E 接触方式∽j 。该接触方式可以自动搜索 接触面,判断接触面并可以处理侵蚀、断裂等复杂边 界变化情况。材料的摩擦系数设定为0 .6 。 3 数值模拟结果分析与比较 3 .1 倒塌过程模拟结果比较 图5 为不同时刻桥体第8 跨倒塌过程。 £ o .0 1s ,桥面预处理结束,桥体在重力作用下 应力重新分布,并集中在拱上立柱及桥墩处。第8 跨右侧桥墩及拱上立柱爆破时间为0 .5 5s ,由于桥 拱、桥面与拱肋形成的桥体结构整体刚性较大,此时 桥体整体处于稳定状态。t 1 .0s ,桥墩左侧拱肋进 行起爆,破坏了桥体的稳定性,由于跨中桥面的预处 理以及此前桥墩和立柱的爆破,此时第8 跨右半跨 在重力作用下开始整体下坠,副拱桥面应力集中且 受拉破坏。随着半秒延迟起爆,第8 跨左侧部分开 始向右侧倾斜。£ 1 .5s ,拱上立柱剩余部分与衬拱 发生接触,两者混凝土单元各有不同程度的侵蚀,桥 体以接触部位为支点形成了活动铰。此时已经能看 b 模型钢筋分布 出左侧半跨有叠加在右侧半跨之上的趋势,这也是 方案设计所需要的。扛1 .8s ,随着两个半跨逐渐以 各自的支点偏转,此时左侧半跨已经叠加在右侧半 跨上并开始接触。扛1 .9s ,第8 跨整体倒塌结束, 跨体分为两半并相互叠加。由于拱上立柱爆破部分 较多,结构触地冲击作用使得桥面破坏的较多。第 8 跨整个倒塌模拟过程从开始桥墩及拱上立柱起爆 到完全倒塌历时1 .3 5s 。 对模拟结果分析可知,预切缝采用0 .5m 可以 实现爆破方案设计的效果,保证跨体相互叠加彻底 倒塌,并且拱上立柱爆破范围采用1 .6m 可使得桥 体解体效果比较好。 3 .2 第8 跨右半跨衬拱混凝土、钢筋单元应力分析 如图6 所示,图6 为右半跨衬拱混凝土、钢筋单 元应力时程曲线。 图中单元2 5 8 9 8 0 所在部位为右半跨衬拱上表 面,提取其混凝土单元应力和与之连接的钢筋单元 轴向应力进行分析。t o .0 1s ,桥体切缝预处理完 成,混凝土单元初始受较小的压力,随着桥体在重力 作用下应力重新分布,单元由受压转为受拉。 t O .5 5s 时,第8 跨右侧桥墩及拱上立柱起爆,衬 拱上表面混凝土单元在重力作用下开始受压,在 0 .9s 达到其最大值1 .8M P 。随着1 .0s 桥墩左侧 拱肋起爆,右半跨整体下坠,拱肋单元触地侵蚀失 效,此时衬拱上表面混凝土单元持续受拉,最大值 8 .9M P ,发生在1 .3 5s ,达到其失效阀值,单元受拉 屈服从结构中删除;与混凝土连接部位钢筋单元最 万方数据 第2 9 卷第4 期宋天培,谢春明,杨军公路双曲拱桥结构爆破拆除数值模拟 9 3 大拉应力发生在1 .4s ,最大值为3 1M P a ,最大压应 力发生在1 .4 5s ,最大值o .2 2G P a 。由于钢筋材料 的可承受应力远大于混凝土,因此,桥体在倒塌过程 孽 ; 古一2 童 一4 蠢 一6 娶 一8 ~一1 0 f 1 .8 08 中,大部分钢筋单元并不会发生屈服,分析可知,分 离式共节点模型可以考虑钢筋与混凝土的不同受力 状态,进而对结构提供更精确的分析。 f 1 .9 08 图5 第8 跨倒塌过程 F i g .5n ec o l l 印s ep r o c e s so fd l ee i g l l t l l8 p 肌 E l e m e n tn o 彳2 5 8 9 8 0 j “s t | a 混凝土单元应力时间曲线 b 钢筋单元轴向应力曲线 图6 第8 跨右半跨混凝土、钢筋单元应力时程曲线 F i g .6 T h es t r e s s t i m ec u r v eo fc ∞c r e t ea n d 糟b 盯e l e m e n to nt l l er i g h th a l f _ s p a nl i m 站a 阳ho ft l l ee i g l l t l ls 耻mo np l 肌o n e 3 .3 结构顶点y 向位移分析 图7 为结构顶点倒塌过程y 向时间位移曲线。 t 0 .5 5s 爆破切口形成后,桥体在重力作用下开始 出现明显的下坠趋势,并以桥墩为支点前倾倒塌, y 向位移线性增加。t 1 .9s ,桥跨整体倒塌结束, 曲线斜率为O ,y 向位移不再增加,最大值5 .7 5m 。 位移曲线斜率在倒塌过程中保持稳定,表明结构倒 塌过程完全符合设计方案。 a 模型结构顶点 b 结构顶点】,向位移曲线 图7 结构顶点y 向位移分析 F i g .7 弧ey d i r e c t i o nd i s p l a c e m e n t 粕a l y s i s o f t I l et 叩o fs t m c t u r e 3 .4 结构堆散范围分析 起爆后桥体在起爆后迅速失稳逐次坍塌,完全 符合设计方案所需达到的效果。图8 所示为模拟中 双曲拱桥爆破拆除倒塌后的爆堆高度。桥体最高点 离地距离为3 .2 4m ,拱顶桥面距地2 .0 4 2m 。对拱 桥爆后剩余桥体现场测量得出桥体最高点离地距离 为2 .9 5m ,拱顶桥面距地1 .8 2m ,与模拟结果相比, 误差分别控制在8 .9 %与1 0 .8 %。因此可见通过计 算机数值模拟,可以预先估计结构的堆散高度,起到 辅助工程设计的作用。 4 爆破效果 采用此爆破方案,该桥梁爆破成功,桥梁按照设 计方案由东至西、从上向下顺利逐次逐拱原地坍塌。 从爆破后的现场观察,桥面、桥梁按设计要求解体、 坍塌、破碎。爆堆平坦,高度与方案设计误差较小, 爆渣块度较小,便于爆后机械或人工清渣。爆破后 周围建筑物和结构物完好,飞石控制很好,其下游方 向大部分爆渣抛掷在1 5m 左右,上游爆渣抛掷在 1 0m 左右,少量的不超过5 0m ,爆破完全达到了预 想效果。桥梁爆破坍塌的整个过程与数值模拟结果 极为相似,完全达到了预期的爆破效果,因此可证明 筋加”m2加晒 O O O 0 O .O萎6目『藿dl蓉鑫毫{毒 万方数据 9 4爆破 2 0 1 2 年1 2 月 拱桥爆破拆除数值模拟的预见性和可行性。 参考文献 R e f e 心n c 姻 r 嗡 一一_ ●■o 登磊拳_ 兰灞 n 1 篱芋瓮搿暑篇路桥组厶爆破拆除技术 寸 N m l C H IE n .舳.n e8 0mu n d 盯一s u p p o r t e da 咒h 曲m a d 图8 爆堆高度模拟结果 单位m t i 哩裔l o ,2 7 1 ;7 2 .7 } i nc h i n e s e “⋯ n g 葛 胁e8 1 m u 胁1 0 nr e 8 u 1 坞mm o d e l [ 2 ]谢先启,韩传伟,贾永胜,等.宜春大桥爆破拆除[ J ] . 8 眦kh 8 i g h 2 u r I i ‘m 爆破,2 0 1 2 ,2 9 1 6 2 .6 8 . [ 2 ] x I E 【i 蛐- q i ,H A NC h 啪w e i ,J L AY o n g - s h e n g ,e ta 1 .E c - 5 结论 p I o s i v ed e m o l i t i 。no fY i c h 吼晡d g e [ J ] .B l a s t i n g ,2 0 1 2 , 。一曼翌警值模跫之翌场壁餐除分堑妻望篓例u 姜盖嚣亲譬蒹;.大型石拱桥坍塌爆破拆 拱顶部预切缝O 5m 足以保证桥身相互叠加彻底倒 除[ J ] .爆破,2 0 0 5 ,2 2 2 6 1 .6 3 . 塌,模拟结果可用来指导爆破方案设计,减小施工强 [ 3 ] c A oY u e ,P A N GR u i .辨n ,Y A N Gc h 锄鼬。n 。e ta 1 .D e m . 度与作业时间。 o l i t i ∞0 fal a 噼8 t o 鹏眦h 嫡d g eb yc o u a p s eb l 踮t i n g 2 通过对副拱立柱爆破范围1 .6m 在倒塌时 [ J ] .B 妇咄,2 ∞5 ,2 2 2 6 l 6 3 . i nc h i n 啪 间和爆堆高度上的分析,可知加大副拱立柱爆破范 [ 4 ] 郭天天,杜耀志,马海鹏.钢筋混凝土公路桥爆破拆除 里譬蹩煅苎氅黛嬲墨罂⋯㈣裂嚣喜嚣’戮芸凇⋯ 果,从而减小结构堆散高度和整体倒塌时间。 。。藻苫,J ≤≤二≮2 焉3 一 ”“一⋯“⋯“’”“ 3 分析了支撑区混凝土和钢筋单元的受力过 [ 5 ] D E N G i ∞.q 嘲,ⅪA o 踟.w u ,L UY 抽.h o n g ,e ta 1 .‰ 程,共节点分离式模型可以体现混凝土和钢筋材料仃o u e db l 枷l l go fas t 0 ∞a r c hb r i d g em 州J ] .m 碰I l g , 的力学性能差异。数值模拟的结果与实际工程具有2 0 0 4 ,2 l 2 6 3 6 5 . i nC I l i 胝圯 较好的相似性,通过对结构倒塌过程的模拟,可以在 [ 6 ]L i V e m o 陀s 疵恤‰I I l l o l o g yc o r p 嘣t i o n L s D Y N A 爆破前对其倒塌过程、效果进行预测,从而指导爆破u ∞r ’sM 彻u a l [ M ] V e 鹞‘∞9 7 0 ,2 0 0 3 A l - A 5 翌躲黧黧些擎璧篓翌璧望竺曹望竺“ 1 茗蒜霎篇兰爹裟破拆除数值模拟研究 破拆除力学过程的重要手段并辅助指导结构爆破拆 [ 8 ] I 忑Y i 瑟;;磊i .A 磊s Y s /L S D Y N A 培训手册 除设计。 [ M ] .2 0 0 2 . 上接第4 l 页 [ 1 2 ] T A Y L O RLM ,c H E NEP ,K u s z M A u LJs .M i c r o c r a c k i n d 删d 锄a g ea c c u m I I l a t i o ni nb r i m er o c ku n d e rd y - n 锄i cl o a d i n g [ J ] .c 眦n p u t e rM e t l l o di nA P PM e c h & E n g n g ,1 9 9 6 ,5 5 3 0 l - 3 2 0 . [ 1 3 ]夏祥,李俊如,李海波,等.爆破荷载作用下岩体振动 特征的数值模拟[ [ J ] .岩土力学,2 0 0 5 ,拍 1 融5 6 . [ 1 3 ] ⅪA 【i 帅g ,uJ 蛐- m ,uH a i b o ,e ta I .u d e cm o d e l i 唱o f v i b r a t i ∞c h a r a c t e r i s t i c 8o fj o i n t e dr o c kI n a s su n d e re x . p 1 0 8 i o n [ J ] .R o c k 肌ds o i lM e c h 肌i c s ,2 0 0 5 ,2 6 1 5 0 - 5 6 . i nC h i n e ∞ [ 1 4 ] 邱流潮,金峰.地震分析中人工边界处理与地震动 输入方法研究[ J ] .岩土力学,2 0 0 6 ,2 7 9 1 5 0 2 一 1 5 0 4 . [ 1 4 】Q I uL i u c h ∞。J I NF e n g .s t u d yo fI 他t h o do fe a n h q u a k e [ 1 5 ] [ 1 5 ] [ 1 6 ] i n p u t 舳da r t m c i a lb 0 岫l l d a r yc o n d i t i o n sf o r ∞i s m i cs o i I s 砸l c t u mi n t e 腿c t i o n 卸衄I y s i s [ J ] .R o c k 锄d 鼬lM e c h 蛐. i c s ,2 0 0 6 ,2 7 9 1 5 0 2 - 1 5 0 4 . i nC h i ∞∞ 黄胜,陈卫忠,杨建平,等.地下工程地震动力响应 及抗震研究[ J ] .岩石力学与工程学报,2 0 0 9 ,2 8 3 4 8 3 .4 9 l o . H U A N GS h e n g ,C H E NW e i - z h o 唱,Y A N GJ i 舭- p i n g ,e t a 1 .R e s e a r c h 伽e a n I I q l l a k e - i n d 眦e dd y I I 锄i c 燃p 蚰螂 雠das e i s ⅡI i cm e 鹊u 他f o r 呻d e r g r 饥n de n g i n e 捌唱[ J ] . C h i n e s eJ 叫m a l0 fR o c kM e c l 啪i c 8 蛐dE n g i n e e r i n g , 2 0 0 9 ,2 8 3 4 8 3 4 9 0 . i nC h i n 鹅e 中华人民共和国行业标准编写组.J ’m I y 7 0 2 0 0 4 公 路隧道设计规范[ S ] .北京人民交通出版社,2 0 0 4 . 万方数据