高文樂(lè)，劉志成，李坤鵬，史卓金

(1.山東科技大學(xué) 土木工程與建筑學(xué)院，青島 266590;2.中鐵十四局建筑分公司，濟(jì)南 250014)

隨著國(guó)家對(duì)環(huán)境保護(hù)的重視，一些高污染、高耗能的火力發(fā)電被關(guān)?；蛘哌M(jìn)行升級(jí)改造，許多廢棄的高聳煙囪等構(gòu)筑物需要拆除。往往受場(chǎng)地大小的限制，單向低位爆破切口的爆破形式已經(jīng)不能滿足倒塌周圍空間的要求，因此折疊爆破拆除的形式應(yīng)用越來(lái)越多。為了確保煙囪爆破拆除的經(jīng)濟(jì)和安全，許多學(xué)者通過(guò)數(shù)值模擬對(duì)煙囪的爆破拆除的過(guò)程進(jìn)行了計(jì)算。其中薛江波對(duì)磚煙囪破壞分析之后，通過(guò)分離式建模進(jìn)行數(shù)值模擬的方式對(duì)鋼筋混凝土煙囪單向倒塌爆破拆除進(jìn)行計(jì)算，數(shù)值模擬得出的數(shù)據(jù)與在爆破現(xiàn)場(chǎng)測(cè)到的數(shù)據(jù)高度一致[1]。Munjiza和Friedhelm采集煙囪爆破過(guò)程中爆堆范圍以及煙囪整體觸地速度的數(shù)據(jù)，與數(shù)值模擬進(jìn)行對(duì)比分析，試驗(yàn)證明數(shù)值模擬的模型和前處理數(shù)據(jù)可以用到工程實(shí)際中[2,3]。武哲分析了處于復(fù)雜環(huán)境中的一棟90 m高框架剪力墻結(jié)構(gòu)樓房采用三個(gè)爆破切口進(jìn)行爆破拆除的技術(shù)，最終得到了理想的爆破效果[4]。任江通過(guò)對(duì)煙囪的結(jié)構(gòu)分析，考慮煙道和灰斗板對(duì)煙囪倒塌的不利影響，通過(guò)封堵現(xiàn)有煙道，應(yīng)用復(fù)式梯形的切口形式，開(kāi)設(shè)定向窗和導(dǎo)向窗對(duì)爆破方向精準(zhǔn)定位等，最終實(shí)現(xiàn)了煙囪按照預(yù)定方向倒塌[5]。葉家明通過(guò)應(yīng)力條件計(jì)算出水塔爆破缺口的圓心角大小和剛度條件計(jì)算出水塔爆破缺口高度，并基于ANSYS/LS-DYNA軟件，對(duì)水塔的定向倒塌進(jìn)行數(shù)值模擬，數(shù)值模擬分析也驗(yàn)證了水塔不能倒塌的結(jié)果[6]。采用基于LS-DYNA聯(lián)合Hypermesh的方法對(duì)壽光巨能熱電廠100m高的煙囪進(jìn)行了數(shù)值模擬，并進(jìn)一步利用數(shù)值模擬的方法對(duì)上下爆破切口的最佳起爆時(shí)差進(jìn)行研究，對(duì)相似工程爆破設(shè)計(jì)和施工具有較大參考價(jià)值。

1 工程概述

壽光市侯鎮(zhèn)南部巨能熱電廠，是一座廢棄電廠，內(nèi)部發(fā)電設(shè)備已經(jīng)拆除。其院內(nèi)有一高為100 m的鋼筋混凝土煙囪需進(jìn)行爆破拆除。周圍環(huán)境為：煙囪東距3層廠房50 m，南距魯麗熱電廠120 m、本廠平房80 m，西、北距最近廠房70 m。如圖1。

圖 1 爆破周圍環(huán)境示意圖(單位：m)Fig. 1 Surrounding environment of blasting(unit:m)

煙囪頂部排煙口直徑3.5 m，壁厚160 mm，煙囪底部外直徑8.3 m、壁厚360 mm?；炷翢焽鑳?nèi)部設(shè)置有耐火磚內(nèi)襯和隔熱層，其中隔熱層厚50 mm，內(nèi)襯在5 m以下的厚度為230 mm，5 m以上為120 mm；豎向配筋為0～15 m為φ22、15～50 m為φ20、50～80 m為φ16，底部為雙層鋼筋網(wǎng)片結(jié)構(gòu)，上部為單層鋼筋網(wǎng)片結(jié)構(gòu)。煙囪底部基座為直徑15 m，深1.8 m鋼筋混凝土結(jié)構(gòu)[7]。

2 爆破方案及爆破參數(shù)

2.1 爆破拆除方案的選擇

煙囪的爆破拆除主要有定向倒塌、原地坍塌和折疊倒塌三種方式，折疊倒塌又分為異向和同向折疊倒塌[8]。待拆煙囪在煙囪向下的沖擊力不夠的情況下，原地坍塌方案很難對(duì)煙囪的倒塌方向進(jìn)行控制；同時(shí)，周圍距離明顯無(wú)法滿足高度100 m煙囪的定向倒塌爆破；根據(jù)待拆煙囪的周圍環(huán)境以及對(duì)煙囪結(jié)構(gòu)的分析，選擇折疊倒塌爆破方案。由于作業(yè)場(chǎng)地大小的限制，爆破拆除時(shí)不僅在煙囪底部設(shè)置爆破切口，還要在其中部的適當(dāng)位置炸開(kāi)一個(gè)或一個(gè)以上的切口，使煙囪朝向兩個(gè)或兩個(gè)以上的方向進(jìn)行同向或反向分段折疊倒塌。在正南正北方向上距離最近建筑物的距離大于70 m，滿足倒塌范圍需求，有足夠的安全距離，故選用雙向折疊倒塌爆破拆除[9]。即充分利用中間平臺(tái)，分+0～+50 m和50～100 m兩段，分別設(shè)計(jì)上下兩個(gè)爆破切口，上部向正北、下部向正南方向倒塌的雙向折疊倒塌。

2.2 爆破切口參數(shù)

(1)切口位置

根據(jù)充分的研究和分析得出，煙囪的下部爆破切口設(shè)置為0.5 m以上，此位置外徑為8.3 m、壁厚0.36 m，重量為1200.98 t，重心高度為23.6 m。上部切口設(shè)置在煙囪50.5 m處，此位置外徑為5.6 m、壁厚0.24 m，重量為322.37 t，重心高度為18.8 m。

(2)切口高度

切口高度是保證定向倒塌的一個(gè)重要參數(shù)。一般情況下[1]，爆破切口的高度通常取為h=(1.5～2.0)δ。為確保爆破時(shí)煙囪按預(yù)定方向準(zhǔn)確倒塌和減小爆堆范圍，實(shí)際工程中切口高度通常取大于切口位置的壁厚5倍以上，因此得出其下部切口h=5.0×0.36=1.8 m，取3 m；煙囪上部切口取h=5×0.24=1.2 m，考慮上切口打眼施工困難，實(shí)際取1 m。

(3)切口長(zhǎng)度

煙囪倒塌方向和爆堆范圍較大地受到切口長(zhǎng)度的影響。為取最佳的切口長(zhǎng)度，應(yīng)根據(jù)式(1)進(jìn)行取值

(1)

式中：L為爆破切口長(zhǎng)度；s為煙囪爆破部位的外周長(zhǎng)。

煙囪的爆破切口圓心角取α=210°，下部的切口長(zhǎng)度取14 m，上部切口長(zhǎng)度取10 m。

(4)上下切口的起爆時(shí)差

折疊爆破達(dá)到最佳效果應(yīng)當(dāng)滿足：保證下切口起爆時(shí)，防止上、下段筒體間出現(xiàn)傾倒方向偏差。

通過(guò)計(jì)算切口高度可知，在切口閉合瞬間，上段筒體重心在切口水平面的投影到塔心的距離大于或等于該處半徑[10]，因此，起爆時(shí)差為

(2)

式中：φ為煙囪與Y軸的夾角；J1為轉(zhuǎn)動(dòng)慣量；L1為煙囪計(jì)算長(zhǎng)度；Δt為經(jīng)驗(yàn)值，一般取100～500 ms(筒體材料強(qiáng)度高時(shí)取小值)。

目前，國(guó)內(nèi)上、下爆破切口的起爆時(shí)差選1.4～7.0 s，根據(jù)理論和工程經(jīng)驗(yàn)分析得出[7]，實(shí)際爆破采用半秒6段雷管，即在上部爆破煙囪形成5°的傾斜時(shí)下部切口開(kāi)始起爆[11]，上部與下部的起爆時(shí)差約為2.5 s。

3 煙囪傾倒過(guò)程的動(dòng)能分析

煙囪折疊倒塌過(guò)程中形成了鉸接雙連桿折疊下落的現(xiàn)象，為雙自由度動(dòng)力學(xué)的問(wèn)題，采用拉格朗日方程建立動(dòng)力學(xué)模型[12]。見(jiàn)圖2。

圖 2 煙囪折疊倒塌過(guò)程圖Fig. 2 Schematic diagram of the collapse of the chimney

(1)上段筒體的幾何學(xué)關(guān)系

x1=h2·cosθ2-L1·sinθ1

(3)

Z1=h2·cosθ2+L1·cosθ1

(4)

式中：X1為煙囪上段質(zhì)心水平位移；Z1為煙囪上段質(zhì)心垂直位移；h1為上段筒體長(zhǎng)度；h2為下段筒體長(zhǎng)度；L1為上段筒體質(zhì)心到上部切口的距離；L2為下段筒體質(zhì)心到下部切口的距離；θ1為上段筒體轉(zhuǎn)角；θ2為下段筒體轉(zhuǎn)角。

將X1，Z1進(jìn)行微分，可得出煙囪上段質(zhì)心水平向及垂直向速度分量，從而得出上段平動(dòng)部分動(dòng)能

(5)

(6)

(7)

m1為上段筒體的質(zhì)量，而上段繞其重心轉(zhuǎn)動(dòng)的轉(zhuǎn)動(dòng)部分動(dòng)能為

(8)

J1c為上段繞其質(zhì)心的轉(zhuǎn)動(dòng)慣量。

由式(7)和式(8)兩式可得上段的總動(dòng)能為

(9)

(2)下段筒體的動(dòng)能為

(10)

J2為下段筒體繞其質(zhì)心的轉(zhuǎn)動(dòng)慣量。

(3)從而得到運(yùn)動(dòng)過(guò)程中系統(tǒng)的總動(dòng)能為

T=T1+T2

(11)

(12)

4 計(jì)算模型

采用SolidWorks建立精準(zhǔn)三維立體模型，之后導(dǎo)入到Hypermesh軟件中對(duì)其進(jìn)行細(xì)化處理。利用Hypermesh軟件強(qiáng)大的網(wǎng)格劃分功能，對(duì)模型網(wǎng)格進(jìn)行加密劃分。考慮到煙囪模型關(guān)鍵地方的約束較多，后期需要添加大量的約束，采用直接法建立鋼筋混凝土共節(jié)點(diǎn)分離式模型[13]，將鋼筋節(jié)點(diǎn)與混凝土節(jié)點(diǎn)直接產(chǎn)出單元體。模型材料參數(shù)設(shè)定中，鋼筋設(shè)置為BEAM161，混凝土設(shè)置為SOLID164 實(shí)體單元。據(jù)鋼筋混凝土結(jié)構(gòu)的力學(xué)特征研究得出，將混凝土模型設(shè)置為* MAT_BRITTLE_DAMAGE ，混凝土的彈性模量為20 GPa，泊松比為0.2，拉伸極限強(qiáng)度為5.0 MPa，剪切極限強(qiáng)度為6.0 MPa，鋼筋模型設(shè)置為*MAT _ PLASTIC _KINEMATIC ，其密度為7825 kg /m3，彈性模量為200 GPa，泊松比為0.25，屈服強(qiáng)度為335 MPa，爆破切口失效模型設(shè)置為* MAT_ADD_EROSION，將鋼筋的抗拉失效應(yīng)變定義為0.15。

通過(guò)修改關(guān)鍵字中的CONTROL_TIMESTEP中的DT2MS控制計(jì)算時(shí)間，當(dāng)DT2MS>0時(shí)，模型中離散的所有小單元都會(huì)以軟件內(nèi)部設(shè)定好的時(shí)間步長(zhǎng)進(jìn)行計(jì)算，全部的小單元都會(huì)被質(zhì)量縮放。當(dāng)DT2MS<0時(shí)，并不是所有的小單元都會(huì)質(zhì)量縮放了，只有時(shí)間步長(zhǎng)小于DT2MS的單元才會(huì)被質(zhì)量縮放[14]。

5 上下切口最佳延時(shí)時(shí)差研究

前面已經(jīng)對(duì)該爆破過(guò)程進(jìn)行了數(shù)值模擬，模擬結(jié)果與實(shí)際爆破結(jié)果高度一致。在此基礎(chǔ)上對(duì)其進(jìn)一步研究，分別設(shè)置上下切口的起爆時(shí)差為1s、1.5 s、2 s、2.5 s、3 s，通過(guò)爆堆最遠(yuǎn)距離、煙囪倒地破碎程度和煙囪觸地速度這三個(gè)變量對(duì)起爆時(shí)差進(jìn)行優(yōu)化。

(1)倒塌過(guò)程的對(duì)比(見(jiàn)圖3)

圖 3 不同起爆時(shí)差下的倒塌過(guò)程模擬效果比較Fig. 3 Comparison of simulation effects of collapse process with different delay time

通過(guò)對(duì)比分析圖3，可以看出隨著上下切口起爆時(shí)差的增長(zhǎng)，煙囪上半段的倒塌角度越小，上半段觸地的時(shí)間越晚。分析倒塌過(guò)程圖可以看出上部爆破切口部分的閉合受到下部爆破切口的影響，當(dāng)上下部爆破切口的起爆時(shí)差相差較小時(shí)，煙囪上下兩段的倒塌方向相反。煙囪下半段倒塌會(huì)通過(guò)爆破切口保留部分的煙囪向反方向拉動(dòng)上半部分，導(dǎo)致上半部分落地速度較快。當(dāng)起爆時(shí)差為3 s時(shí)，上部和下部爆破切口起爆時(shí)差相差較大，上部煙囪爆破切口部分的鋼筋混凝土已經(jīng)被拉壞，下半段煙囪對(duì)上半段的倒塌速度影響較小，上部分煙囪斷裂后下部分煙囪才開(kāi)始傾斜，是兩個(gè)獨(dú)立的爆破過(guò)程。上下部煙囪爆破倒塌沒(méi)有銜接起來(lái)，會(huì)導(dǎo)致爆破倒塌方向的不確定性增大[15]。

(2)爆堆倒塌最遠(yuǎn)距離的對(duì)比

通過(guò)分析圖4不同起爆時(shí)差下的爆堆范圍及倒塌最遠(yuǎn)距離得出不同起爆時(shí)差下的煙囪觸地后基本破碎，只有起爆時(shí)差為2.5 s時(shí)的倒塌最遠(yuǎn)距離最小為50 m，其他起爆時(shí)差的倒塌最遠(yuǎn)距離均為55 m左右。因此起爆時(shí)差為2.5 s時(shí)的倒塌范圍是最小的，即為最安全的爆堆范圍。

(3)煙囪觸地速度的對(duì)比

通過(guò)圖5不同起爆時(shí)差下的煙囪觸地速度可以分析得出：最大觸地速度發(fā)生時(shí)刻隨著起爆時(shí)差的增大而越來(lái)越大，起爆時(shí)差為1 s、1.5 s、2 s、2.5 s、3 s下的最大觸地速度出現(xiàn)時(shí)刻分別為6.5 s、6.9s、7.5 s、7.9 s、8.2 s，同時(shí)證明圖3中得出的隨著上下切口起爆時(shí)差的增大觸地時(shí)間越來(lái)越晚的結(jié)論正確。起爆時(shí)差為1 s的煙囪整體最大觸地速度為16 m/s，其余起爆時(shí)差的最大整體觸地速度隨起爆時(shí)差增大而逐漸減小。

6 數(shù)值模擬結(jié)果及分析

分析各項(xiàng)指標(biāo)，雙向折疊倒塌爆破的上下切口的最優(yōu)起爆時(shí)差為2.5 s，和爆破設(shè)計(jì)的公式計(jì)算得出的上下爆破切口的起爆間隔時(shí)差相同。圖6顯示了煙囪倒塌與工程實(shí)際的比較，數(shù)值模擬結(jié)果精確仿真了工程實(shí)際的倒塌過(guò)程。煙囪倒塌的方向也按照爆破設(shè)計(jì)的方向倒塌，折疊爆破實(shí)現(xiàn)了減小爆堆范圍的要求。0 s時(shí)上部爆破切口起爆，在重力作用下開(kāi)始向著爆破切口方向傾斜。2.5 s時(shí)下部爆破切口起爆，此時(shí)上部煙囪大概傾斜了5°，減小了煙囪出現(xiàn)前沖現(xiàn)象的概率。4.2 s時(shí)上下段的煙囪傾斜角度不斷加大，傾斜速度不斷加快。6.6 s時(shí)近乎自由落體式的速度觸地破壞，實(shí)際的煙囪倒塌過(guò)程與數(shù)值模擬計(jì)算相同時(shí)刻的煙囪倒塌的過(guò)程基本一致。

圖 5 不同起爆時(shí)差下的煙囪整體觸地速度Fig. 5 The overall velocity of the chimney touching the ground with different delay time

圖 6 實(shí)際爆破效果與模擬效果的比較Fig. 6 Comparison of actual blasting effect and simulation results

結(jié)合數(shù)值模擬的結(jié)果可以將煙囪的倒塌過(guò)程分為三個(gè)階段：第一階段為上切口的形成至保留區(qū)出現(xiàn)裂縫的階段；鋼筋和混凝土共同承擔(dān)上部煙囪的壓力作用、上部切口形成而產(chǎn)生彎矩，出現(xiàn)拉應(yīng)力作用。此時(shí)處于彈性階段，應(yīng)力與應(yīng)變呈現(xiàn)線性正比例關(guān)系。隨著彎矩的增大，爆破切口保留處的應(yīng)變也越來(lái)越大，由于混凝土抗拉強(qiáng)度遠(yuǎn)小于抗壓強(qiáng)度，混凝土在拉應(yīng)力作用下開(kāi)始出現(xiàn)塑性變形，但鋼筋彈性模量較大，其仍然屬于彈性階段[16]。當(dāng)上部切口保留處的拉應(yīng)力大于混凝土所能承受的最大拉應(yīng)力時(shí)，爆破切口保留處開(kāi)始出現(xiàn)裂縫。裂縫的出現(xiàn)標(biāo)志進(jìn)入第二階段，隨著彎矩和變形的增加混凝土逐漸退出工作，在受拉鋼筋達(dá)到屈服之前，所有的拉應(yīng)力由鋼筋承擔(dān)。接下來(lái)進(jìn)入破壞階段，為第三階段，鋼筋達(dá)到屈服值以后進(jìn)入強(qiáng)化和頸縮階段，此時(shí)裂縫快速變大、上部煙囪傾斜角逐漸增大直到從爆破切口處斷開(kāi)觸地破壞，下部切口爆破后也會(huì)產(chǎn)生和上部切口相同的三個(gè)階段[17]。

7 結(jié)論

(1)基于LS-DYNA聯(lián)合Hypermesh的數(shù)值模擬方法對(duì)100 m高鋼筋混凝土煙囪的雙向折疊爆破進(jìn)行了仿真模擬，數(shù)值模擬的結(jié)果與工程現(xiàn)場(chǎng)數(shù)據(jù)高度一致。

(2)分析了煙囪倒塌的觸地快慢受到起爆時(shí)差的影響，得出起爆時(shí)差越大，觸地時(shí)間就越晚的結(jié)論。同時(shí)起爆時(shí)差越大，上部爆破切口倒塌方向的爆堆就越大。爆破時(shí)差為2.5 s時(shí)的爆堆堆積范圍是最小的。

(3)分析了煙囪的觸地速度，得出整體觸地速度越快對(duì)地面的振動(dòng)也就越大，越不安全的結(jié)論。分析得出起爆時(shí)差為1 s時(shí)的煙囪整體觸地速度最大為16 m/s，其余起爆時(shí)差下煙囪整體觸地最大速度隨起爆時(shí)差增大而逐漸減小。也就意味著起爆時(shí)差越大對(duì)地面的振動(dòng)越小，越有利于周邊建筑物的安全。

(4)綜合分析模擬結(jié)果，爆破時(shí)差為2.5 s時(shí)爆堆的范圍最小，煙囪整體倒塌觸地速度也相對(duì)較小為13.8m/s，爆堆破碎程度與時(shí)差影響關(guān)系不明顯，所以2.5 s是最佳的起爆時(shí)差。從而也證明了爆破時(shí)差計(jì)算公式是正確可行的，對(duì)類似的爆破工程具有一定的指導(dǎo)意義。