费鸿禄,张超逸,吕金齐,张龙飞

(1.辽宁工程技术大学爆破技术研究院,辽宁 阜新 123000;2.内蒙古宏大爆破工程有限责任公司,内蒙古 包头 014010)

在冷却塔爆破拆除倒塌中，由于塔体倾覆造成塔体发生下坐现象。对于在转动过程中的下坐问题可从两个角度进行分析：(1)保留区人字柱的最大承载能力，即认为当塔体在转动过程中，保留区人字柱所承受的最大垂向荷载大于人字柱可提供的最大承载能力，此时保留区人字柱受压破坏造成塔体失稳下坐；(2)保留区截面塑性铰转角达到转动极值，即当爆破切口形成瞬间，保留区人字柱与塔体接触部分形成塑性铰。当塔体转动角度大于塑性铰转角极值时，塑性铰部位发生压剪破坏而造成塔体下坐。通过观察以往冷却塔爆破拆除效果[1]，可以发现当塔体完全倒塌后保留区部分人字柱未发生整体受压破坏，而是柱端破碎。由此证明应用塑性铰理论分析塔体下坐前的偏转角度更符合工程实际。

鉴于此,本文中根据保留区人字柱的破坏形态对塑性铰转动机理进行分析，利用达朗贝尔原理和拉格朗日方程得到塑性铰转动极值，从而对冷却塔整个倾倒过程中的运动状态进行深入研究，并且通过建立塔体触地瞬间的力学模型分析塔壁的破碎情况及其影响因素；最后基于现场振动监测结果，对降低拆除过程中的振动危害提出合理化建议。

1 冷却塔倾倒过程及解体力学分析

冷却塔从爆破切口形成瞬间，到塔体触地前，塔体在倾覆力矩作用下绕中性轴做定轴转动[2]。在转动过程中，由于塔体底部圈梁作用，保证了冷却塔的整体性，因此认为在倾倒过程中塔体完好未发生解体破坏。在这个过程中，保留区人字柱与塔体圈梁接触部分形成塑性铰。为分析塔体运动机理，对塑性铰的转动极限进行分析。由以往冷却塔拆除现场观测情况[1]可知，由于圈梁截面尺寸较大且承载能力较高，因此在倾倒过程中破坏主要发生在人字柱与圈梁连接部位。因此判定塔体倾倒时保留区形成的塑性铰为弱柱型塑性铰[3]。为简化计算，假定在转动过程中，中性轴位置保持不变。

1.1 塔体倾倒塑性铰分析

设截面塑性铰转角极值为[θ]，则其计算公式为[4]：

图1 冷却塔受力分析简图Fig.1 Force analysis of the cooling tower

[θ]=θy+θu

(1)

式中：θy为屈服状态时塔体转动角度，θu为极限状态时塔体转动角度。起爆后塔体开始转动，倾覆角度增大，切口附近人字柱逐渐受压破坏失效，截面抗剪强度随人字柱的逐渐破坏而逐渐降低，当截面受到的剪应力达到人字柱截面提供的抗剪强度时，人字柱柱端开始发生压剪破坏[5]，塑性铰转动能力丧失[6]，塔体开始下坐，则认为此时倾覆角度达到混凝土塑性铰的转角极值[θ]。塑性铰转角极值越大，塔体开始下坐的时间越晚，则塔体在触地瞬间转动角度越大。转动过程中塔体受力分析简图如图1所示。

根据达朗贝尔原理，列出平衡方程：

(2)

(3)

根据动能定理列出方程：

式中：Mα为偏转到任意角度α时的塑性铰弯矩。则:

将上述计算结果代入式(2),得到当塔体转动时截面剪力关于转角的计算公式：

(4)

根据上述分析得到的塔体运动规律和几何关系，列出方程求解塔体从开始下坐到触地所经历的时间t和触地瞬间塔体转动角度θ：

(5)

式中：R为塔体底部外半径，e2为塔身重心至人字柱截面中性轴的偏心距。因t与θ的表达式过于繁琐，不便直接列出，因此代入MATLAB中进行求解分析。

1.2 塔体触地瞬间受力分析

在冷却塔与地面碰撞的过程中，塔体处于非剧烈冲击状态[8]，计算这类低速碰撞中塔壁钢筋混凝土的整体冲击作用时，可利用经典模型(弹性、刚塑性模型)以准静态进行分析[9-10]。在这种状态下塔体耗能远高于地面，因此忽略在碰撞过程中地面的能量损失。假设塔体触地瞬间，塔体静止且无任何运动趋势,塔体触地瞬间计算计算简图如图2所示。由于开设了卸荷槽，在塔体触地瞬间塔体的触地状态为塔壁柱触地[11]。因此，为简化分析，假设在塔体触地瞬间，塔体与地面间的冲击力均匀地作用在中间卸荷槽两侧的塔壁柱底部。

图2 触地瞬间塔体碰撞力分析简图Fig.2 Sketch for calculating the collision force of the tower body at the moment of touchdown

图2中C为塔体质心；vC为触地瞬间塔体质心速度；θ为塔体触地瞬间的偏转角度;P为塔体触地瞬间的碰撞力，其方向与塔体质心速度方向平行；Px和Py分别为碰撞力在水平和竖直方向的分量，根据几何关系列出方程：

(6)

1.3 冷却塔触地瞬间解体力学分析

由于开设了高卸荷槽，降低了切口下部塔壁的整体强度，而且卸荷槽顶部发生应力集中，因此在塔体触地碰撞瞬间，塔壁必将在卸荷槽顶部发生破碎[11]。但上部塔体完好，其是否因受冲击而发生破碎解体仍为未知。现任取上部塔壁一垂直于母线的截面为研究对象，且设该截面在母线轴上的坐标为L。由于在低峰值、长持时的准静态荷载作用下，钢筋混凝土倾向于发生整体性的弯曲破坏[12]，因此确定截面在倒塌方向的中心点为最易破坏点。则该点应力状态为：

(7)

塔壁混凝土受压而产生脆性断裂破坏，因此选用最大线应变理论[13]作为判断塔壁破坏的力学依据，为方便表述，设σ为等效应力，[σ]为钢筋混凝土许用拉应力，经整理得到：

(8)

式中：ν为混凝土的泊松比。将式(6)与(7)代入式(8)得到等效应力σ的表达式，因其表达式过于繁琐，不能直接列出，因此需要代入塔体具体参数并利用MATLAB软件进行求解，并对塔壁破碎条件及破碎位置进行分析。

2 工程实例

某冷却塔为钢筋混凝土结构，总高度为125 m，人字柱高度h1=8.328 m。塔体下人字柱共40对横截面直径为0.65 m的圆柱。筒壁设计为双曲线形，底部半径为42.746 m，最顶部半径为27.56 m，中间位于高程91.847 m处的圆周为双曲线的顶点，该处半径最小为25.365 m。通风筒的底部圈梁厚度为0.7 m，自底部至高程+26.815 m渐变至0.18 m，自此处以上的通风筒壁厚均为0.18 m。本次冷却塔拆除采用只对人字柱进行爆破并在塔体开设高卸荷槽的复合型切口的爆破方案。爆破区分为3个区，切口高度为26.8 m，圆心角为216°，对应需爆破的人字柱为24对，爆破切口和分区如图3所示。

图3 冷却塔爆破切口与爆破区示意图Fig.3 Diagram of blasting incision and area of the cooling tower

2.1 塔体倒塌过程分析

观察实际冷却塔爆破过程(图4)可知，冷却塔倒塌过程分为2个部分。首先，爆破切口形成瞬间，塔体在倾覆力矩的作用下前倾；在塔体前端触地前，塔体偏转至一定角度后，塔体开始下坐；在下坐过程中，伴随着塔体偏转。这与利用塑性铰理论分析塔体运动状态得到的结论相吻合。塔体触地后，塔壁上有2个位置的钢筋混凝土同时受到挤压而破碎，且破碎范围随着塔体偏转幅度的增大而增大。而观测冷却塔倒塌视频，测量得到上部破碎截面位置在母线轴上的坐标约为46.589 m，下部破碎截面在母线轴上的坐标为18.5 m。而后，破碎区下方塔体受纵向冲量作用导致其混凝土不断被压溃[14]，从而丧失承载能力，诱发塔体连续倒塌。

图4 冷却塔爆破倒塌过程图Fig.4 The blasting collapse process of the cooling tower

2.2 触地瞬间塔壁解体位置分析

根据塔体设计尺寸，将塔体外、内半径与母线长度拟合成函数关系：

(9)

(10)

图5 MATLAB数值计算结果Fig.5 Result of numerical calculation by MATLAB

将上述塔体参数代入式(8)中，利用MATLAB求解，求解结果如图5所示。

图5中，L为塔壁截面在母线轴上的坐标，[θ]为塔体塑性铰偏转角度。由图5可知，MATLAB数值计算结果为一曲面。这说明，在触地瞬间，塔壁截面等效应力值与母线轴上的坐标值及塔体塑性铰转角极值有关。在塔体触地瞬间，塔壁在应力值最大的截面处发生破坏，而后应力释放，其余位置保持完好。当18.488 m

为分析塔体触地瞬间上部塔体破碎位置与塑性铰转角极值[θ]之间关系，则令dσ/dL=0，则得到坐标值L与塑性铰转角极值[θ]之间的函数关系，利用MATLAB软件进行求解，求解结果如图6所示。根据图6可知，上部塔壁破碎截面在母线轴上的坐标随塑性铰转角极值的增大而减小。

图6 塔体破碎截面坐标与塑性铰转角极值的关系Fig.6 Relationship between the broken section coordinate and the deflection angle of the tower body

3 冷却塔爆破拆除振动监测及分析

为确定本次冷却塔爆破拆除对周围重点建筑物的影响，采用TC-4850测振仪对冷却塔爆破及触地塌落振动进行现场监测。监测点分别位于冷却塔南侧28 m的电厂医院处与东南侧47m的住宅楼处，振动监测结果如图7所示。

图7 两个检测点处质点振动速度的变化Fig.7 Particle vibration velocity-time curves at two vibration monitoring points

根据现场监测结果可知，采用高卸荷槽复式切口的冷却塔爆破拆除诱发的地表振动可分为3个部分。第1部分为0～1 s，该部分振动是由炸药爆炸引起的；第2部分为2.5～3.5 s，该部分振动是由爆破切口范围内塔壁柱陆续触地所诱发的；第3部分为4.5～10.0 s，该部分为塔体连续塌落诱发的塌落振动，这个过程中振动持时较长并且由以往研究成果[15]可知该过程的振动频率与周围建筑物频率相近，因此认为冷却塔连续倒塌时对周围建筑物的振动危害最严重[16]。通过比较现场监测结果与倒塌视频可知，在连续塌落过程地表的振速峰值是由塔壁解体后上部塔体塌落触地所诱发的。因此，为了有效降低整个拆除过程中产生的振动危害，应从控制连续塌落过程中的振速峰值入手。

根据建筑物爆破拆除时的塌落振动计算公式[17]：

(11)

式中：vt为塌落引起的地表振动速度；M为建筑物的质量；H为建筑物重心下落高度[18]；σd为建筑物材料的破坏强度；D为塌落中心距观测点的距离；Kt、β分别为衰减系数和指数，其中β为负值。

(12)

式中：l2为冷却塔总高度减去人字柱高度。由式(12)可知，当被爆物确定之后，解体后上部塔体的重力做功与塔壁解体位置L及塔体触地瞬间转角θ有关。经上述分析可知，塔体触地瞬间转动角度θ随塑性铰转角极值[θ]的增大而增大，即dθ/d[θ]>0；塔壁解体位置L随塑性铰转角极值[θ]增大而减小，即dL/d[θ]<0，因此dL/dθ=(dL/d[θ])(d[θ]/dθ)<0，证明塔壁解体位置L随塔体触地瞬间转角θ的增大而减少，则将式(12)对转动角度θ求导得：

(13)

对式(13)进行分析，由塔体参数可知，(c1-c2)πρg为正值，而由上述证明可知dL/dθ<0，则-(2aL2+2bL+c1+c2)(Rsinθ+Lcosθ)dL/dθ为正值。因此为判断d(MgH)/dθ的正负，只需对Rcosθ-lsinθ的数值进行分析，将塔体参数代入MATLAB中进行求解，求解结果如图8所示。

图8 MATLAB数值计算结果Fig.8 Result of numerical calculation by MATLAB

由图8可知，Rcosθ-lsinθ的数值恒为正，因此d(MgH)/dθ>0，即解体后的上部塔体在连续塌落过程中的重力做功随塔体触地瞬间的偏转角度增大而增大。由此证明，解体后上部塔体在连续塌落过程中诱发的地表振速峰值随触地瞬间塔体转动角度的增大而增大。由人字柱截面受力状态可知，当第一区起爆后，余下人字柱截面即承受偏心荷载。根据冷却塔倾覆条件[19]可知，在爆破切口形成的瞬间，塔体重力产生的倾覆力矩大于截面屈服弯矩。因此在第3区起爆前，塔体的转动角度小于屈服状态时塔体的转动角度，则冷却塔不会下坐且转动角度随整体延期时间的增长而增大。由此判断，当爆破切口圆心角确定之后，保留区塑性铰转动极值和冷却塔触地瞬间的转动角度随整体延期时间的增长而增大。因此，在今后冷却塔爆破工程中应避免过长的整体延期时间，从而减小塔体触地瞬间的偏转角，降低冷却塔爆破对周围建筑的振动危害。

4 结 论

(1)根据保留区人字柱的破坏形态，判定塔体倾倒时保留区形成的塑性铰为弱柱型塑性铰。当塔体偏转角度达到塑性铰转动极值后，保留区人字柱柱端发生压剪破坏并丧失承载能力，塔体开始下坐并伴随着塔体转动。在触地瞬间，塔体的转动角度随塑性铰转角极值的增大而增大。

(2)通过建立塔体数学模型，以最大线应变理论作为塔体触地破碎的力学依据，分析塔体触地瞬间塔壁截面应力状态，求得等效应力最大位置，比较该处应力与许用拉应力的大小，得出塔体在碰撞瞬间除切口处塔壁发生破碎外，塔壁上部也将发生解体破碎，这与冷却塔实际倒塌过程中破碎现象相吻合。当被爆物确定后，上部塔壁破碎截面在塔体母线上的坐标值随塑性铰转动极值的增大而减小。

(3)冷却塔连续倒塌过程诱发的塌落振动对周围建筑的危害最严重，而其振速峰值随塔体触地瞬间转动角度的增大而增大。鉴于此，对今后的冷却塔爆破拆除施工，在保证冷却塔正常倾覆条件的前提下，应避免爆破延期时间过长，从而控制触地瞬间塔体的转动角度，降低冷却塔爆破对周围建筑的危害。