Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

钢筋混凝土烟囱爆破拆除的下坐及早期断裂预测

孙金山 谢先启 贾永胜 姚颖康 刘昌邦 韩传伟 王洪刚 黄小武

孙金山, 谢先启, 贾永胜, 姚颖康, 刘昌邦, 韩传伟, 王洪刚, 黄小武. 钢筋混凝土烟囱爆破拆除的下坐及早期断裂预测[J]. 爆炸与冲击, 2022, 42(8): 085202. doi: 10.11883/bzycj-2021-0316
引用本文: 孙金山, 谢先启, 贾永胜, 姚颖康, 刘昌邦, 韩传伟, 王洪刚, 黄小武. 钢筋混凝土烟囱爆破拆除的下坐及早期断裂预测[J]. 爆炸与冲击, 2022, 42(8): 085202. doi: 10.11883/bzycj-2021-0316
SUN Jinshan, XIE Xianqi, JIA Yongsheng, YAO Yingkang, LIU Changbang, HAN Chuanwei, WANG Honggang, HUANG Xiaowu. Prediction of sinking down and early break in the air of reinforced concrete chimney during blasting demolition[J]. Explosion And Shock Waves, 2022, 42(8): 085202. doi: 10.11883/bzycj-2021-0316
Citation: SUN Jinshan, XIE Xianqi, JIA Yongsheng, YAO Yingkang, LIU Changbang, HAN Chuanwei, WANG Honggang, HUANG Xiaowu. Prediction of sinking down and early break in the air of reinforced concrete chimney during blasting demolition[J]. Explosion And Shock Waves, 2022, 42(8): 085202. doi: 10.11883/bzycj-2021-0316

钢筋混凝土烟囱爆破拆除的下坐及早期断裂预测

doi: 10.11883/bzycj-2021-0316
基金项目: 湖北省自然科学基金(2020CFA043);湖北省重点研发计划(2020BCA084);江汉大学科技创新专项
详细信息
    作者简介:

    孙金山(1980- ),男,博士,教授,sun99001@126.com

  • 中图分类号: O389; TU746.5

Prediction of sinking down and early break in the air of reinforced concrete chimney during blasting demolition

  • 摘要: 为分析钢筋混凝土烟囱在爆破拆除时发生下坐与空中断裂现象的机制并对其进行预测,对一高180 m烟囱的下坐和空中断裂过程进行了观测和分析。基于混凝土的压缩全应力-应变曲线特征,分析了烟囱支撑区的破坏过程,构建了烟囱失稳下坐的判别模型。通过建立烟囱下坐冲击作用下爆破切口以上烟囱的动力响应模型,分析了下坐冲击附加动应变波在烟囱中的传播特征。研究结果表明,考虑混凝土全应力-应变曲线特征和支撑区横截面应力和应变分布特征时,倾覆力矩与抵抗力矩的比值f可作为失稳下坐的判别条件之一;烟囱发生下坐的必要条件是支撑区最小残余承载力小于烟囱的重量。烟囱在下坐结束阶段,获得一定初速度的烟囱冲击基础时将产生冲击荷载,并在烟囱中部引起大于底端应变的应变,即产生动应变高程放大效应,该效应是导致烟囱发生早期断裂的主要原因。烟囱越高,下坐冲击历时越短,动应变高程放大效应越显著,发生断裂的风险也越大。随着烟囱高度的增加,烟囱最危险截面的位置也越高:由烟囱中下部移至烟囱中上部。
  • 目前,废弃的钢筋混凝土烟囱主要采用爆破技术拆除。近年来,在拆除约180 m的高烟囱时,部分烟囱会发生严重的下坐问题,并在起爆后不久发生中部断裂,且上、下段有时会发生分离,上段烟囱的倒塌方向可能失控,而由于下段烟囱的重心位置发生了改变,容易导致烟囱下半段向相反方向倒塌或“炸而不倒”(如图12所示),该现象可能导致严重的安全事故。

    图  1  萧山热电厂180 m高烟囱
    Figure  1.  The 180-m-high chimneyof Xiaoshan thermal power plant
    图  2  成都热电厂210 m高烟囱
    Figure  2.  The 210-m-high chimneyof Chengdu thermal power plant

    烟囱爆破拆除过程涉及爆破切口形成、支撑区破坏、应力重分布和结构动力响应等诸多方面,其力学机制较复杂。针对烟囱爆破切口形成后支撑部位的破坏机制问题,褚怀保等[1]观测了爆破切口形成后保留筒壁的受力状态,认为爆破切口形成后存在0.5~3.0 s的荷载重新分布和中性轴形成的过程。郑炳旭等[2]观测了6座钢筋混凝土烟囱切口的支撑区,认为支撑区在爆破后承受自重突加的荷载,容易引起烟囱下坐;切口端先受压破坏,而后中性轴受拉,呈现大偏心受压脆性断裂特征。郑炳旭等[3]还分析了切口自重突加载荷引起的支撑区受压范围,认为切口圆心角宜取210°~230°。徐鹏飞等[4]认为爆破切口形成后2~3 s的中性轴稳定时间是烟囱预防过早下坐和形成定向倾倒趋势的关键。言志信等[5]建立了钢筋混凝土烟囱爆破后支撑区的应力计算模型,提出引入冲压系数来考虑突加荷载的影响。

    针对高烟囱爆破过程中的“空中折断”问题,杨建华等[6]建立了烟囱任一截面的内力分布和极限承载力模型,认为高度超过150 m的钢筋混凝土烟囱在倾倒角度超过40°后,均可能在离顶部约1/3高度处发生断裂,烟囱越高,折断发生的时间越早。言志信等[7]认为烟囱折断发生的位置和时间与切口形状及材料的力学性能密切相关。唐海等[8]认为烟囱在倾倒过程中主要发生弯曲破坏,首次折断的部位约在距离顶部1/3高度处,强度不大的烟囱可能会有多次折断。侯吉旋等[9]认为质量均匀分布烟囱的断裂点距离顶部1/3处,而对于上细下粗的烟囱,断裂点将会下移。

    针对烟囱结构在纵波作用下的动力响应问题,王云剑[10]实验研究了烟囱断裂位置与冲击波作用周期和烟囱固有周期之间的关系。Pallarés等[11]采用三维有限元模型分析了砌体烟囱在地震作用下的破坏现象,得到了结构的破坏模式、最大应力和位移特征。Wolf等[12]研究了核电站典型烟囱在地震和冲击荷载作用下的响应特征。Wilson[13]根据10座烟囱在地震作用下的非线性特征,提出了钢筋混凝土烟囱的非线性动力分析方法。Huang等[14]根据高115 m钢筋混凝土烟囱的地震动力响应,提出了一种新的三维推覆分析方法。Minghini等 [15]对地震中砖砌烟囱破坏问题进行了分析,阐述了烟囱上部的剪切破坏机制。

    综上所述,对钢筋混凝土烟囱爆破拆除的失稳、倒塌、运动过程以及烟囱地震响应等问题已有大量研究,但对烟囱在爆破拆除过程中的下坐及其诱发的早期断裂问题研究较少。本文中将通过对一座高180 m的钢筋混凝土烟囱的运动和断裂过程的观测,对烟囱爆破拆除的下坐和早期空中断裂现象进行分析讨论,研究失稳、下坐和早期断裂的判别或预测方法。

    爆破拆除的钢筋混凝土烟囱高180 m,0~25 m高程混凝土标号为C40,25~180 m高程混凝土标号为C30,烟囱主要结构尺寸见表1。烟囱横截面轴向配筋为双层配筋,环向配水平箍筋。7.25 m高程处分布有2个烟道口并设有积灰平台,大烟道尺寸为5.9 m×10.4 m,小烟道尺寸为5.9 m×5.4 m;东、西方向各有2个检修门,尺寸均为2.4 m×2.4 m,如图3所示。

    表  1  烟囱主要结构尺寸
    Table  1.  Structure parameters of the chimney
    高程/m筒壁外半径/cm筒壁内半径/cm筒壁厚度/cm隔热层厚度/cm内衬厚度/cm
    08127575500
    7.25776721551023
    20.00712662501023
    30.00662617451023
    45.00617575421012
    60.00572533391012
    75.00527491361012
    90.00482449331012
    105.00460430301012
    120.00437410271012
    135.00415391241012
    150.00392370201012
    165.00392372201012
    180.00392372201012
    下载: 导出CSV 
    | 显示表格
    图  3  爆破方案
    Figure  3.  Blasting plan of the chimney

    烟囱采用正梯形爆破切口,如图3所示。切口高4.4 m,布设在烟囱底部0.5 m高程处,爆破切口的圆心角为220°,底边展开长30.95 m。切口两侧布设定向窗,底边长3.0 m,张开度为30°。爆破前用钢筋混凝土对西侧检修门进行封堵。爆破切口区共布置634个炮孔,乳化炸药总装药量为136 kg,炮孔内安装延期时间为3 400 ms的导爆管雷管。

    为分析烟囱的失稳破坏过程,在烟囱支撑区外侧布置一套动态摄影测量系统。该摄影测量系统由高精度工业相机、基准尺、测量标志、计算软件和电脑组成,先利用测量相机采集目标点的xyz坐标值,再利用目标点做点线面的标准拟合等。将2台工业相机对称安装在支撑区背侧地面的支座上(见图3),相机分辨率4872×3248,采集频率1000 s−1。经标定,动态摄影测量系统的测量误差标准值小于3 mm。

    同时,在支撑区外围地面上布置2个分辨率1280×960、采集频率25 s−1的普通监控摄像头。在爆破远区,则采用无人机在空中对爆破过程进行视频录像。

    由爆破过程的监测视频可知,起爆后切口两侧的混凝土受到强烈挤压而发生破坏。起爆后约0.5 s,支撑区产生了与水平方向夹角为45°的裂缝(见图4(a)),伴随着主裂缝的扩展,大量的混凝土不断从筒壁挤出、脱落;起爆后约1.2 s,裂缝沿着混凝土薄弱部位扩展(见图4(b));起爆后约1.8 s,支撑区的主裂缝贯通(见图4(c));起爆后约2.5 s,爆破切口完全闭合(见图4(d))。

    图  4  烟囱支撑区裂纹扩展过程
    Figure  4.  Crack propagation process in the support part

    监测结果显示,在支撑区后部裂缝扩展较快,混凝土破坏也较严重。在裂缝产生、扩展至贯通的过程中,支撑区后部未观测到明显的拉裂缝产生,表明该烟囱支撑区的破坏方式主要是压剪破坏。

    烟囱贯穿裂缝形成后,支撑区无法承受上部结构的荷载而发生整体下坐。下坐过程中,烟囱筒体不断与底部残余结构及地面发生剧烈碰撞,筒体底部混凝土被压碎、挤出,并堆积在烟囱周围。起爆约4.0 s后,烟囱下坐过程结束。整个下坐过程中,烟囱下坐速度先增大,然后迅速减速至0,共历时约2 s,下坐总高度约10 m。与此同时,烟囱在下坐过程中还产生了轻微转动,但总体转动角度不大,下座结束时烟囱转动角度大约为3°。

    通过摄影测量,获得了烟囱下坐过程的位移时程曲线,如图5所示。测量数据表明,烟囱爆破后下坐量先随时间缓慢增长,随后快速增加,最后则逐渐停止。将下坐位移时程曲线对时间t进行一次和二次微分后,可以近似获得下坐速度(见图6)和加速度时程曲线(见图7)。在下坐开始后的1.0~1.5 s,下坐速度可达10 m/s,而向下运动的加速度最大可接近10 m/s2。下坐运动减速时,向上的加速度最大可接近30 m/s2。因此,烟囱在下坐过程中可能经历多次“失重”和“超重”效应。而“失重”和“超重”的变换速度较快,将在烟囱中诱发压缩应力波,可能造成烟囱的持续下坐。当下坐量较大时,底部破碎的钢筋混凝土能起到显著的缓冲作用,“超重”效应逐渐减弱,直至下坐停止。

    图  5  烟囱下坐位移时程曲线
    Figure  5.  Displacement-time history curve of chimney sinking down
    图  6  烟囱下坐速度时程曲线
    Figure  6.  Velocity-time history curve of chimney sinking down
    图  7  烟囱下坐加速度时程曲线
    Figure  7.  Acceleration-time history curve of chimney sinking down

    烟囱在下坐结束时,在其中部发生了断裂,断裂位置约在90 m高程处。断裂后的上半段筒体继续向原倒塌方向发生转动。起爆后约12 s时,烟囱上、下两部分筒体分离,上部分筒体在其原初速度的基础上加速下落;下部分筒体由于重心位置未偏出筒体的投影区且倾倒动能不足,导致无法继续倾倒(见图8(d))。

    图  8  烟囱空中折断过程
    Figure  8.  Breaking in the air of the chimney

    根据烟囱支撑区的受力特征与现场观测结果,同时考虑混凝土受压破坏过程的全应力-应变曲线特征(见图9),支撑区的失稳下坐过程可分为以下4个阶段。

    图  9  混凝土典型应力-应变曲线
    Figure  9.  Typical full strain-stress curve of concrete

    (1)起爆后应力瞬态调整

    爆破切口形成瞬间,支撑区将发生应力瞬态调整,爆破切口部分承受的荷载将迅速向支撑区混凝土上转移。

    (2)支撑区正截面抗弯承载力达到峰值

    应力瞬态调整完成后,保留的支撑区变为大偏心受压构件(见图10(a))。钢筋与混凝土将发挥材料最大承载性能以抵抗倾覆力矩和竖向压力。根据大偏心受压构件的受力变形特点,在支撑区横截面上距离中性轴最远位置处的应力和应变最大,同时,由于混凝土的屈服应变远低于钢筋的屈服应变,靠近定向窗的支撑区混凝土最先发生压缩屈服,此时整个截面的抵抗力矩达到峰值或接近峰值(见图10(b))。若支撑区抵抗力矩大于倾覆力矩则混凝土将不再进一步破坏,烟囱则不能顺利倾倒;若支撑区抵抗力矩小于倾覆力矩则混凝土将进一步破坏,烟囱则可继续发生转动。

    图  10  支撑区不同区域混凝土的应力与应变状态示意图
    Figure  10.  Schematic of the stress and strain status of the concrete in the support part

    (3)支撑区正截面受压承载力达到峰值

    当烟囱失稳倾覆时,随着其转动角度的增大,支撑区截面上的应变将持续增大,且中性轴位置发生变化。当支撑区中性轴消失时,整个支撑区横截面均呈受压状态。由混凝土的全应力-应变曲线特征可知,混凝土达到屈服状态后,随着应变的增大其承载能力将不断降低直至达到一定的残余强度,变为塑性材料。因此,中性轴刚刚消失时,截面上各处混凝土的压应力分布与其应变相对应,分别处于混凝土全应力-应变曲线的不同应力值,此时支撑区断面达到其最大竖向承载力(见图10(c))。

    (4)支撑区塑性铰的形成或下坐

    受倾覆力矩的作用,混凝土的应变还将持续增大,因此一般情况下整个支撑区最终将变成“塑性铰”,此时,支撑区的混凝土仅剩余残余承载力(见图10(d)~(f))。当烟囱的重量小于或等于支撑区的残余承载力时,烟囱则匀速破坏、缓慢下坐、持续转动。当烟囱重量大于支撑区竖直向残余承载力时,烟囱将获得一定的竖向加速度进而向下加速下坐,并同步发生转动,但由于下坐速度快、历时短,下坐过程中烟囱偏转角度一般较小。

    爆破拆除烟囱时,为了提高爆破切口设计的可靠性,需对不同设计方案进行对比和失稳判别。为了简化计算,传统的失稳判别力学模型通常假定受拉区和受压区完全达到极限承载状态,这与支撑区的实际受力状态存在显著差异。支撑区中性轴附近的应力和应变均较小,而在前后边缘处应变最大(但应力不一定最大)。

    在起爆后的短时间内,支撑区筒壁处于低应变状态,因此,可认为中性轴处竖向应变为0,中性轴两侧混凝土的应变与到中性轴的距离呈线性关系。即混凝土受压区外边缘屈服应变为εc,将逐渐减小至中性轴处的0。

    图10(a)所示,设烟囱支撑区横截面上单侧受压区(受压区的1/2)对应的弧度为β,其在极坐标系下的起始弧度为α,单侧受拉区(受拉区的1/2)对应的弧度则为π/2−αβ,其起始弧度为α+β。设筒壁外半径为r1、内半径为r2,则在受压区中极角为α+ϕ处混凝土应变εϕ可表示为:

    εϕ = εcˉr[sin(α+β)sinα]ˉr[sin(α+β)sinϕ] (1)

    式中: ˉr为筒壁外半径r1和内半径r2的平均值,ϕ为支撑区横截面某点的极角弧度值。

    设支撑区筒壁微元体的弧度为dϕ,微元体面积则为δˉrdϕ,则受压区承受竖向压力的合力Ppc为:

    Ppc = α+βαεϕEcδˉrdϕ=α+βαsin(α + β)sinϕsin(α + β)sinαεcEcδˉrdϕ (2)

    式中:δ为支撑区筒壁的厚度,Ec为混凝土的弹性模量。

    同理,已知支撑区筒壁单根钢筋的极角时,也可得到单根钢筋的压力、拉力及其合力。但是在爆破设计过程中钢筋的定位难度较大,逐根钢筋的受力计算也较繁琐,因此可将钢筋网等效为薄壁圆筒进行计算。忽略钢筋保护层厚度时,将钢筋网等效为钢筒后,等效钢筒壁厚δs1δs2可表示为:

    {δs1=As1Δs1δs2=As2Δs2 (3)

    式中:As1As2为支撑区底部截面外侧和内侧竖向主筋截面面积,Δs1Δs2为外侧和内侧钢筋的间距。

    类似式(2),可得内、外侧钢筋筒受到的竖向合力Pps为:

    Pps = α+βαsin(α + β)sinϕsin(α + β)sinαεcEsδs1r1dϕ+α+βαsin(α + β)sinϕsin(α + β)sinαεcEsδs2r2dϕ (4)

    式中:Es为钢筋的弹性模量。

    同理,支撑区的受拉区混凝土受到的竖向拉力Ptc为:

    Ptc = π/2α+βsinϕsin(α + β)sin(α + β)sinαεcEcδˉrdϕ (5)

    受拉区钢筋受到的竖向拉力Pts为:

    Pts = π/2α+βsinϕsin(α + β)sin(α + β)sinαεcEsδs1r1dϕ+π/2α + βsinϕsin(α + β)sin(α + β)sinαεcEsδs2r2dϕ (6)

    需要指出的是,爆破后切口区裸露的钢筋骨架也具有一定的承载能力,但在工程实践中爆破切口高度往往较高,且爆炸荷载会使竖向钢筋发生弯曲。因此,爆破后切口区的钢筋骨架竖向承载能力较低时,其贡献的竖向承载力可忽略。

    假定应力重新调整后,烟囱处于临界失稳状态,且受压区边缘应变恰好达到屈服应变,此时烟囱处于静力平衡状态,且未发生加速转动,则烟囱的重量G与支撑区竖直向力的合力为0,得到平衡方程:

    2Ppc + 2Pps2Ptc2Pts = G (7)

    式中:G为烟囱的总重量。

    将式(2)、(4)~(6)代入式(7)并化简得:

    εc2cscβ2sec(α + β2)[(π 2α)sin(α + β)2cosα][δˉrEc + Es(r1δs1+r2δs2)]=G (8)

    式(8)为超越方程,无解析解。为了方便工程应用,利用泰勒级数对三角函数进行展开,化简后解得β的近似解为:

    βV+W+V210U(π 2α) (9)

    式中:

    U=δˉrEcεc + Esεc(r1δs1+r2δs2)
    V=48G5π U(8+3α)+4U(3+20α+6α2)
    W=40U2(π2α)(4+α)[4(3+α2)5πα]

    在工程设计过程中,先初步确定α后可通过试算法解得方程(8)的中受压区对应的弧度β,或通过近似解(9)计算β,进而可确定中性轴位置,并可进一步对烟囱能否失稳进行验算。

    同样,在与式(7)相同的基本假定下,设θ = α + β,1/2受压区混凝土所产生的抵抗力矩Mpc为:

    Mpc=2(αθ)(cos2θ2)+2cosα(sinα4sinθ)+3sin2θ4(sinθsinα)εcEcδˉr2 (10)

    将钢筋等效为钢筒后,1/2受压区钢筋产生的抵抗力矩Mps为:

    Mps = 4θ4α+2(αθ)cos2θ+2cosα(sinα4sinθ)+3sin2θ4(sinθsinα)(εcEsδs1r12 + εcEsδs2r22) (11)

    1/2受拉区钢筋产生的极限抵抗力矩Mts为:

    Mts = (π 2θ)(cos2θ2)+3sin2θ4(sinαsinθ)(εcEsδs1r12 + εcEsδs2r22) (12)

    烟囱要打破静力平衡状态而发生加速转动时,重力形成的倾覆力矩要大于抵抗力矩。忽略受拉区混凝土产生的抵抗力矩时,烟囱的失稳条件为:

    Gˉrsinθ2Mpc + 2Mps + 2Mts (13)

    定义倾覆失稳系数f为:

    f = Gˉrsinθ2Mpc + 2Mps + 2Mts (14)

    由式(10)~(12)、(14)可计算烟囱的失稳系数f

    根据混凝土应力应变特征及其在支撑区横截面上的分布规律,设支撑区边缘处混凝土刚进入塑性状态时,f=1;假定受压区全部进入塑性状态且混凝土材料为理想弹塑性材料时,f=2;假定受压区一半面积进入塑性状态且混凝土为理想弹塑性材料时,f=1.5。显然,f=1.5~2时可满足失稳的要求。而f>2时倾覆力矩过大,可能导致支撑区快速压溃而发生下坐。

    根据3.1节的分析,烟囱支撑区满足倾覆失稳条件后,将发生加速转动,而由于钢筋混凝土的不断屈服,其抵抗力矩将变得越来越小,因此截面的中性轴将不断发生移动以发挥材料的最大承载能力,当钢筋最终被拉断或者中性轴消失时,支撑区正截面将主要受压缩作用。

    混凝土在受压屈服后,承载能力将随着应变的增大而持续降低,最终将达到其残余承载力。由于烟囱爆破切口常常采用正梯形切口,会形成变截面的支撑区,因此设支撑区下侧最小截面面积为Ac,min,上侧最大截面面积为Ac,max,则支撑区残余承载力与烟囱重量存在以下3种关系。

    (1)烟囱重量大于支撑区最大残余承载力时,会发生显著加速下坐直至支撑区完全压溃消失。支撑区和爆破切口消失后,如新形成的断面较为平整,则烟囱重心在水平面的投影落在新截面的形心附近,倾覆力矩将急剧减小甚至消失,烟囱可能发生炸而不倒的现象,其判别条件表示为:

    GAc,maxσcr (15)

    式中:σcr为支撑区混凝土的残余强度,Ac,max为支撑区上侧最大截面面积(采用正梯形切口时)。

    (2)烟囱重量大于支撑区最小残余承载力,且小于最大残余承载力时,烟囱先发生加速下坐直至重量和支撑区残余承载力相等时再发生减速下坐,直至达到新的平衡。此时,支撑区可能部分保留,也可能因烟囱惯性力作用而全部压溃,因此,顺利倾倒和炸而不倒的现象均可能发生,其判别条件表示为:

    Ac,minσcrGAc,maxσcr (16)

    式中:Ac,min为支撑区下侧最小截面面积(采用正梯形切口时)。

    (3)烟囱重量小于支撑区最小残余承载力时,烟囱因支撑区混凝土的压溃挤出而发生向下的位移,可能会发生短暂的下坐,但支撑区残余截面的不断增大将阻止下坐的进一步发生,使烟囱得以顺利倾倒,其判别条件表示为:

    GAc,minσcr (17)

    烟囱空中断裂现象的宏观特征表明,发生下坐现象时支撑区底部截面将与基础碰撞,碰撞接触面上必然产生较高的动态压应力和应变并向烟囱顶部传播。假设烟囱筒体中有一压缩应力波传播,仅分析烟囱切口以上且远离切口区的部分时,该部分可简化为经典的变截面一维直杆力学模型。同时,仅考虑一个波动信号从底端到达顶端的应力波传播过程时,可忽略波在两端边界处的入射和反射等过程,进而可方便地对其动力学过程进行理论求解。

    因此,当烟囱下坐导致切口闭合、支撑区消失时,在烟囱烟道口以上部分的横截面上将作用一个轴向脉冲荷载,在满足平截面假定条件时,取烟囱的一段微元体进行分析,如图11所示。

    图  11  不受应力集中效应影响的烟道口以上烟囱的微元体模型
    Figure  11.  Microelement model of the chimney above the flue without stress concentration effect

    根据达伦贝尔原理,微元体的受力平衡方程可表示为:

    σA+(σ+σxdx)(A+dAdx)=12ρ[A+(A+dAdx)]dx2ut2 (18)

    式中:σ为烟囱横截面上的竖向应力,A为烟囱的横截面面积。为获得微分方程的解,设烟囱截面随烟囱高度L连续变化,可表示为:

    A=κeλxL (19)

    式中:κλ为系数,x为横截面距离地面的高度,L为烟囱总高度。

    将式(19)代入式(18)中,化简并略去高阶项后,式(18)简化为:

    Aσx+σdAdx=ρA2ut2 (20)

    化简得:

    1A(σA)x=ρ2ut2 (21)

    设钢筋混凝土材料为线弹性,则将σ=Ecux代入式(21),得到波动方程:

    2ux2+1AuxdAdx=1cp22ut2 (22)

    式中:cp为钢筋混凝土中声波波速,cp=Ec/ρ

    在任意时刻,烟囱顶部为自由边界,其应变始终为0,则烟囱上端边界条件为:

    ux|x=L=0 (23)

    t=0时,烟囱底端脉冲荷载为0,任意截面的位移为0;设t时刻烟囱底端的运动速度为vt,则烟囱下端边界条件为:

    {u|t=0=0ut|x=0=vt (24)

    在式(24)中,vt为变量,需确定其数值方可求解波动方程。设t=0时,烟囱的运动速度为v0,并逐渐降低至0。根据碰撞冲击作用过程的特征,且为了便于方程的求解,以三角函数拟合的vt为:

    vt=v02+v02cosπ tt0 (25)

    式中:t0为烟囱速度由v0降低至0时所经历的时间。

    根据3.3节的分析,烟囱在下坐过程中,需克服支撑区的残余承载力而发生加速运动,其加速度a为:

    aGAcσcrm (26)

    式中:Ac为支撑区截面面积(可取平均值)。

    设下坐高度为h,根据运动距离、速度和加速度间的关系,v0可表示为:

    v0=2ah (27)

    由式(26)~(27)可得:

    v0=2hGAcσcrm (28)

    烟囱速度由v0降低至0的过程中,设F为基础提供的减速阻力,根据动量守恒原理得:

    mv0=Ft0 (29)

    烟囱下坐后筒壁与基础的接触面为整个圆环形截面时,接触面上的力先增大,直至达到混凝土的动态屈服强度后,再减小至烟囱的重量。因此,取接触力的近似平均值时,F的取值范围为:

    F(ησcA0+G2ησcA0) (30)

    式中:σc为混凝土的抗压强度,η为混凝土的抗压强度动态提高系数,A0为烟囱底端筒壁的横截面面积。由式(28)~(30)得到t0的取值范围为:

    t0(mv0ησcA02mv0ησcA0+G) (31)

    以式(23)~(31)为边界条件,可解方程(22)得到u(x,t)的近似解,由ux可得到高度x处的烟囱横截面上的竖向应变。

    根据式(22)~(31),可对下坐一定高度的烟囱进行动力学分析,估算烟囱横截面上的应变值,当一定高度范围内筒壁的应变均大于混凝土的压缩屈服应变时,则可判定烟囱将发生早期空中断裂现象。

    以第1节中的高180 m的烟囱为例对预测模型进行验证。该烟囱质量约8600 t,重量约84.28 MN,截面面积函数中κ=24.5、λ=1.5。混凝土密度取2500 kg/m3,混凝土屈服强度标准值取σc=26.8 MPa,弹性模量取Ec=32.5 GPa,泊松比取0.17,钢筋弹性模量取200 GPa。切口烟囱底端外半径为8.12 m,内半径为7.57 m,壁厚0.55 m。爆破切口设计圆心角ω=220°,支撑区下边缘展开长度约33.56 m,截面面积18.46 m2;上边缘展开长度约18.31 m,截面面积10.07 m2。外侧钢筋截面面积为254.5 mm2,内侧钢筋截面面积为153.9 mm2,间距均为0.2 m。爆破切口参数如图3所示。

    根据式(9)~(10),可在初选切口圆心角后对支撑区中性轴位置进行计算,进一步计算支撑区受压区比例。计算结果如图12所示,圆心角ω取180°~230°时,受压区占整个支撑区的比例约32%~34%。由式(14)计算得到,切口圆心角超过190°时,失稳系数接近满足f>1.5的要求(见图13)。而该工程爆破方案确定的切口圆心角为220°,失稳系数f=2.53,说明切口圆心角偏大,倾覆力矩远大于抵抗力矩,一旦切口角部混凝土失效,受压区必然向受拉区快速扩展,很难避免下坐现象的出现。

    图  12  切口圆心角ω与支撑区受压占支撑区比例ζ的关系
    Figure  12.  Relationship between the blasting notch’s central angle and the compressive region ratio
    图  13  切口圆心角ω与烟囱倾覆失稳系数f的关系
    Figure  13.  Relationship between the blasting notch’s central angle and the instability coefficient

    已知烟囱自重G,根据爆破切口形状取支撑区下部的最小截面面积Ac,min=10.07 m2,上部最大截面面积Ac,max=33.56 m2,取混凝土残余强度σcr=0.2σc=5.36 MPa时,得到:

    G=84.28 MNAc,maxσcr=179.88 MNAc,minσcr=53.98 MN

    因此,烟囱支撑区静态残余承载力满足以下条件:

    Ac,minσcrGAc,maxσcr

    计算结果表明,该工程原设计方案将发生下坐现象,且支撑区横截面面积较小的下半部分的破坏将导致烟囱加速下坐,而横截面面积较大的上半部分则使下坐减速,但当减速产生的加速度超过−9.8 m/s2时,下段支撑区也将承受大于2G的荷载,整个支撑区可能完全被压溃,此时爆破切口将因下坐而完全闭合,倾覆力矩也将消失,烟囱存在“炸而不倒”的风险。因此,预测结果与烟囱实际下坐情况基本吻合。

    设烟囱自重与支撑区残余承载力相等时,对应的支撑区横截面面积为28.58 m2,截面距离切口底边1.44 m,切口总高度5 m,则可能发生加速下坐的高度为2.96 m,近似取3 m进行计算时,由式(28)计算得到加速下坐可能获得的最大初速度v0≈7.67 m/s。分别取η=1,1.5,由式(31)计算得到烟囱下坐由7.67 m/s减速至0时所需的时间区间为:

    t0(0.060 s,0.163 s)

    t0取值区间内分别取不同的时间值,由式(26)计算u(x,t)u(x,t)x

    计算结果表明(见图14),在烟囱下坐结束阶段,若其历时较短(如t0=0.06 s),随着烟囱横截面所处高度的增加,其竖向(烟囱轴向)峰值应变呈现出先增大后减小的趋势,呈现出“应变高程放大效应”。而若下坐结束历时较长(如t0=0.15 s),则“应变高程放大效应”消失,此时随着高度的增加,应变持续减小。

    图  14  减速过程对轴向应变放大系数ξ的影响
    Figure  14.  Distribution of amplification factor of longitudinal strain in deceleration process

    定义烟囱某处横截面与底部冲击断面上轴向峰值应变的比值为该处动应变放大系数ξ,下坐结束阶段历时较短时,峰值应变出现的相对位置较高,例如t0=0.06 s时ξ最大值出现在约90 m高度处, t0=0.1 s时ξ最大值出现在约60 m高度处。

    因此,当t0=0.06~0.1 s时,存在ξ>1的可能,发生空中压溃的概率较高(忽略动应变衰减情况下)。本试验中180 m烟囱空中折断位置大约位于距底端90 m处,与t0=0.06 s的情况较吻合。

    由上述180 m高烟囱下坐诱发空中断裂的案例可知,下坐产生的冲击荷载可造成烟囱的损伤甚至断裂。为了分析其他常见钢筋混凝土烟囱的断裂风险,对5种典型烟囱进行了分析。假定烟囱下坐碰撞基础的最大速度为3 m/s,经历0.06 s后速度降为0。忽略附加动应变沿烟囱轴线逐渐衰减时,计算结果(见表2)表明,随着烟囱高度的增加,烟囱的最大动应变放大系数ξ由1.093增加到1.728。烟囱越高,动应变放大系数越大,发生断裂的风险也越大。并且随着烟囱高度的增加,烟囱最大动应变放大系数所处的相对高度也随之增加,由烟囱中下部提升至烟囱中上部。实际上,当考虑附加动应变随传播距离发生衰减时,高烟囱发生空中断裂的可能性较高,而矮烟囱发生空中断裂的可能性较低。

    表  2  不同高度烟囱的最大动应变放大系数
    Table  2.  Maximum amplification factor of peak dynamic strain of chimneys with different heights
    烟囱高度/m烟囱截面面积随高度变化函数ξ最大动应变所处高度/m
    1 90A = 7.46e−0.017x1.09330
    2120A = 10.34e−0.014x1.24150
    3150A = 20.70e−0.012x1.38380
    4*180A = 22.35e−0.011x1.522110
    5210A = 32.05e−0.011x1.728140
    注: 表中180 m烟囱为一般的等截面烟囱,与第2节案例烟囱形状存在一定差异。
    下载: 导出CSV 
    | 显示表格

    由上述分析可知,避免烟囱下坐或者避免硬着陆是预防烟囱早期空中断裂的重要措施,具体可采用以下方法:

    (1) 优化爆破切口形状,改梯形切口为三角形切口或喇叭形切口,使支撑区截面面积逐渐过渡至整个圆环截面,增大支撑区的极限承载力,避免支撑区被瞬间整体压溃;

    (2) 优化爆破切口的圆心角,在满足倾覆条件时,需考虑支撑区的残余承载力,通过增大支撑区截面面积减小下坐的速度;

    (3) 必要时应对支撑区混凝土进行加固,提高其抗压强度。

    高烟囱爆破拆除过程中,常出现下坐和中部断裂的现象,严重影响爆破安全和爆破效果。结合典型钢筋混凝土烟囱爆破拆除工程,分析了烟囱下坐和空中断裂的宏观特征,并研究了其力学机制,构建了判别力学模型,得到以下结论。

    (1) 传统的烟囱爆破设计方法认为,支撑区轴心受压极限承载力大于烟囱的重量即可预防支撑区的压溃,然而支撑区实际处于大偏心受压状态,因此,正截面轴心受压稳定条件不能作为下坐判别条件。

    (2) 爆破切口形成后,基于混凝土的全应力-应变曲线特征以及应力和应变在支撑区横截面的分布特征、倾覆力矩与抵抗力矩的比值f应作为切口设计的控制条件之一;且f>2时支撑区受压截面面积过小可能使其发生强烈压缩破坏而导致烟囱下坐。

    (3) 大偏心受压作用下支撑区混凝土的压缩破坏不可避免,支撑区最小残余承载力小于烟囱的重量是下坐判别的必要条件。

    (4) 高度超过100 m的混凝土烟囱在下坐结束阶段,其底端与基础碰撞产生冲击荷载作用,在烟囱中部诱发的应变值可能大于其底端的失效应变,即诱发动应变高程放大效应,该现象是导致烟囱发生早期断裂的主要原因。

    (5) 烟囱越高,下坐冲击历时越短,动应变高程放大效应越显著,发生断裂的风险也越大。随着烟囱高度的增大,烟囱最大动应变放大系数所处的相对位置也随之变化,将由烟囱中下部移至烟囱中上部。

  • 图  1  萧山热电厂180 m高烟囱

    Figure  1.  The 180-m-high chimneyof Xiaoshan thermal power plant

    图  2  成都热电厂210 m高烟囱

    Figure  2.  The 210-m-high chimneyof Chengdu thermal power plant

    图  3  爆破方案

    Figure  3.  Blasting plan of the chimney

    图  4  烟囱支撑区裂纹扩展过程

    Figure  4.  Crack propagation process in the support part

    图  5  烟囱下坐位移时程曲线

    Figure  5.  Displacement-time history curve of chimney sinking down

    图  6  烟囱下坐速度时程曲线

    Figure  6.  Velocity-time history curve of chimney sinking down

    图  7  烟囱下坐加速度时程曲线

    Figure  7.  Acceleration-time history curve of chimney sinking down

    图  8  烟囱空中折断过程

    Figure  8.  Breaking in the air of the chimney

    图  9  混凝土典型应力-应变曲线

    Figure  9.  Typical full strain-stress curve of concrete

    图  10  支撑区不同区域混凝土的应力与应变状态示意图

    Figure  10.  Schematic of the stress and strain status of the concrete in the support part

    图  11  不受应力集中效应影响的烟道口以上烟囱的微元体模型

    Figure  11.  Microelement model of the chimney above the flue without stress concentration effect

    图  12  切口圆心角ω与支撑区受压占支撑区比例ζ的关系

    Figure  12.  Relationship between the blasting notch’s central angle and the compressive region ratio

    图  13  切口圆心角ω与烟囱倾覆失稳系数f的关系

    Figure  13.  Relationship between the blasting notch’s central angle and the instability coefficient

    图  14  减速过程对轴向应变放大系数ξ的影响

    Figure  14.  Distribution of amplification factor of longitudinal strain in deceleration process

    表  1  烟囱主要结构尺寸

    Table  1.   Structure parameters of the chimney

    高程/m筒壁外半径/cm筒壁内半径/cm筒壁厚度/cm隔热层厚度/cm内衬厚度/cm
    08127575500
    7.25776721551023
    20.00712662501023
    30.00662617451023
    45.00617575421012
    60.00572533391012
    75.00527491361012
    90.00482449331012
    105.00460430301012
    120.00437410271012
    135.00415391241012
    150.00392370201012
    165.00392372201012
    180.00392372201012
    下载: 导出CSV

    表  2  不同高度烟囱的最大动应变放大系数

    Table  2.   Maximum amplification factor of peak dynamic strain of chimneys with different heights

    烟囱高度/m烟囱截面面积随高度变化函数ξ最大动应变所处高度/m
    1 90A = 7.46e−0.017x1.09330
    2120A = 10.34e−0.014x1.24150
    3150A = 20.70e−0.012x1.38380
    4*180A = 22.35e−0.011x1.522110
    5210A = 32.05e−0.011x1.728140
    注: 表中180 m烟囱为一般的等截面烟囱,与第2节案例烟囱形状存在一定差异。
    下载: 导出CSV
  • [1] 褚怀保, 徐鹏飞, 叶红宇, 等. 钢筋混凝土烟囱爆破拆除倒塌与受力过程研究 [J]. 振动与冲击, 2015, 34(22): 183–186,198. DOI: 10.13465/j.cnki.jvs.2015.22.032.

    CHU H B, XU P F, YE H Y, et al. Collapse process and load-bearing process of reinforced concrete chimney during blasting demolition [J]. Journal of Vibration and Shock, 2015, 34(22): 183–186,198. DOI: 10.13465/j.cnki.jvs.2015.22.032.
    [2] 郑炳旭, 魏晓林, 陈庆寿. 钢筋混凝土高烟囱爆破切口支撑部破坏观测研究 [J]. 岩石力学与工程学报, 2006, 25(S2): 3513–3517. DOI: 10.3321/j.issn:1000-6915.2006.z2.026.

    ZHENG B X, WEI X L, CHEN Q S. Study on damage surveying of cutting-support of high reinforced concrete chimney demolished by blasting [J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(S2): 3513–3517. DOI: 10.3321/j.issn:1000-6915.2006.z2.026.
    [3] 郑炳旭, 魏晓林, 陈庆寿. 钢筋混凝土高烟囱切口支撑部失稳力学分析 [J]. 岩石力学与工程学报, 2007, 26(S1): 3348–3354. DOI: 10.3321/j.issn:1000-6915.2007.z1.114.

    ZHENG B X, WEI X L, CHEN Q S. Mechanical analysis of cutting-support destabilization of high reinforced concrete chimney [J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(S1): 3348–3354. DOI: 10.3321/j.issn:1000-6915.2007.z1.114.
    [4] 徐鹏飞, 刘殿书, 张英才. 烟囱高位组合切口定向爆破倒塌过程数值研究 [J]. 振动与冲击, 2017, 36(15): 265–270. DOI: 10.13465/j.cnki.jvs.2017.15.040.

    XU P F, LIU D S, ZHANG Y C. Numerical study on the directional blasting collapse process of chimney with high combined incision [J]. Journal of Vibration and Shock, 2017, 36(15): 265–270. DOI: 10.13465/j.cnki.jvs.2017.15.040.
    [5] 言志信, 叶振辉, 刘培林, 等. 钢筋混凝土高烟囱定向爆破拆除倒塌过程研究 [J]. 振动与冲击, 2011, 30(9): 197–201. DOI: 10.3969/j.issn.1000-3835.2011.09.041.

    YAN Z X, YE Z H, LIU P L, et al. Collapsing process of high reinforced concrete chimney in blasting demolition [J]. Journal of Vibration and Shock, 2011, 30(9): 197–201. DOI: 10.3969/j.issn.1000-3835.2011.09.041.
    [6] 杨建华, 马玉岩, 卢文波, 等. 高烟囱爆破拆除倾倒折断力学分析 [J]. 岩土力学, 2011, 32(2): 459–464. DOI: 10.3969/j.issn.1000-7598.2011.02.023.

    YANG J H, MA Y Y, LU W B, et al. Analysis of fracture mechanics for falling tall chimneys during demolition blasting [J]. Rock and Soil Mechanics, 2011, 32(2): 459–464. DOI: 10.3969/j.issn.1000-7598.2011.02.023.
    [7] 言志信, 叶振辉, 刘培林. 烟囱定向爆破拆除倒塌过程 [J]. 爆炸与冲击, 2010, 30(6): 607–613. DOI: 10.11883/1001-1455(2010)06-0607-07.

    YAN Z X, YE Z H, LIU P L. Collapsing process of chimney demolition by directional blasting [J]. Explosion and Shock Waves, 2010, 30(6): 607–613. DOI: 10.11883/1001-1455(2010)06-0607-07.
    [8] 唐海, 梁开水, 张成良. 烟囱爆破倾倒折断的力学浅析 [J]. 爆破, 2003, 20(1): 9–11. DOI: 10.3963/j.issn.1001-487X.2003.01.003.

    TANG H, LIANG K S, ZHANG C L. Mechanics analysis of fall-down process of chimney by blasting demolition [J]. Blasting, 2003, 20(1): 9–11. DOI: 10.3963/j.issn.1001-487X.2003.01.003.
    [9] 侯吉旋, 李志昂, 郭兴, 等. 质量非均匀分布的烟囱在倾倒过程中的力学分析 [J]. 大学物理, 2017, 36(6): 50–51,55. DOI: 10.16854/j.cnki.1000-0712.2017.06.013.

    HOU J X, LI Z A, GUO X, et al. Mechanical analysis of the non-uniform falling chimney [J]. College Physics, 2017, 36(6): 50–51,55. DOI: 10.16854/j.cnki.1000-0712.2017.06.013.
    [10] 王云剑. 烟囱纵向冲击断裂试验与分析 [J]. 力学与实践, 2000, 22(2): 41–43. DOI: 10.3969/j.issn.1000-0879.2000.02.012.

    WANG Y J. Longitudinal shock test and analysis on chimney models [J]. Mechanics in Engineering, 2000, 22(2): 41–43. DOI: 10.3969/j.issn.1000-0879.2000.02.012.
    [11] PALLARÉS F J, AGÜERO A, MARTÍN M. Seismic behaviour of industrial masonry chimneys [J]. International Journal of Solids and Structures, 2006, 43(7/8): 2076–2090. DOI: 10.1016/j.ijsolstr.2005.06.014.
    [12] WOLF J P, SKRIKERUD P E. Collapse of chimney caused by earthquake or by aircraft impingement with subsequent impact on reactor building [J]. Nuclear Engineering and Design, 1979, 51(3): 453–472. DOI: 10.1016/0029-5493(79)90133-X.
    [13] WILSON J L. Earthquake response of tall reinforced concrete chimneys [J]. Engineering Structures, 2003, 25(1): 11–24. DOI: 10.1016/S0141-0296(02)00098-6.
    [14] HUANG W, GOULD P L. 3-D pushover analysis of a collapsed reinforced concrete chimney [J]. Finite Elements in Analysis and Design, 2007, 43(11/12): 879–887. DOI: 10.1016/j.finel.2007.05.005.
    [15] MINGHINI F, MILANI G, TRALLI A. Seismic risk assessment of a 50m high masonry chimney using advanced analysis techniques [J]. Engineering Structures, 2014, 69: 255–270. DOI: 10.1016/j.engstruct.2014.03.028.
  • 期刊类型引用(9)

    1. 刘涛,李洪伟,刘伟,梁昊,吴延梦,黄昕旭,储亚坤. 钢混水塔爆破参数的安全校核及数值模拟验证分析. 工程爆破. 2025(01): 74-81 . 百度学术
    2. 袁方,罗鹏,王威,陈德志,黄小武,杨坤,吴赛格. 复杂环境下含人字形内衬的180 m烟囱爆破拆除. 爆破. 2024(01): 120-126 . 百度学术
    3. 孔庆亮,夏治园,王刚,刘明锋,钱明渊,杨帆,高朋飞. 双筒型旧硝铵造粒塔爆破拆除数值模拟与施工设计. 爆破. 2024(02): 151-159 . 百度学术
    4. 高文乐,闫桐庆,孙鸿,寇玉明,刘卫东,徐明,段耀奎. 爆破切口高度对钢筋混凝土烟囱爆破拆除效果的影响. 爆破器材. 2024(04): 45-51+57 . 百度学术
    5. 赵盟乔,陈磊,王猛,罗宁,权树恩,刘桐. 复杂环境下爆破拆除130 m多管式套筒烟囱的数值模拟研究. 爆破. 2024(04): 128-135 . 百度学术
    6. 黄小武,谢先启,贾永胜,刘昌邦,姚颖康,孙金山,伍岳. 钢筋混凝土水塔垂直原地爆破拆除试验研究. 爆破. 2023(02): 1-8+60 . 百度学术
    7. 叶海旺,吴家鹏,郑长青,许松,傅建秋,罗伟涛. 钢筋砼烟囱拆除爆破切口高度优选及效果分析. 工程爆破. 2023(03): 117-122 . 百度学术
    8. 徐顺香,吴赛格,罗鹏,王威,陈德志,谢广波. 复杂环境下110 m冷却塔控制爆破拆除. 爆破. 2023(03): 129-133 . 百度学术
    9. 柴亚博,罗宁,袁翊硕,杨振,韩涛,曹祺. 超高钢筋混凝土烟囱爆破切口角度选取的研究. 爆破. 2023(04): 124-131+160 . 百度学术

    其他类型引用(0)

  • 加载中
图(14) / 表(2)
计量
  • 文章访问数:  549
  • HTML全文浏览量:  218
  • PDF下载量:  100
  • 被引次数: 9
出版历程
  • 收稿日期:  2021-07-27
  • 修回日期:  2022-04-25
  • 网络出版日期:  2022-05-12
  • 刊出日期:  2022-09-09

目录

/

返回文章
返回