羅如登 鮑志斌 王正陽

①中南大學(xué)土木工程學(xué)院(湖南長沙，410000)

②中交二航局建筑科技有限公司(湖北武漢，430000)

引言

工程爆破廣泛應(yīng)用于山區(qū)土木工程基礎(chǔ)開挖項目中，為推動我國基礎(chǔ)設(shè)施建設(shè)發(fā)揮了非常重要的作用。在鉆爆法施工過程中，爆破安全一直是工程爆破行業(yè)關(guān)注的重要問題。采用鉆爆法開挖時，經(jīng)常會面臨對爆破點鄰近目標(biāo)(如既有橋梁)的保護。控制爆破地震波對鄰近橋梁結(jié)構(gòu)安全性的影響，同時又確保在建項目的高效爆破施工，是爆破施工單位、行政主管部門、監(jiān)理單位等多方共同面臨的技術(shù)難題[1-7]。

丁茂瑞[8]利用二維數(shù)值分析研究了新建隧道爆破開挖對鄰近既有橋梁的安全性動力影響。黃赫烜等[9]通過監(jiān)測數(shù)據(jù)，研究了隧道爆破和鉆孔灌注樁施工對鄰近高鐵橋梁的影響，并結(jié)合有限元模型分析了橋梁的安全性和可靠性。孟靈鑫等[10]基于FLAC3D用瞬態(tài)動力學(xué)方法，分析了爆破振動對鄰近橋面及橋墩的影響。李杰[11]在水下爆破中，對鄰近橋梁振動危害進行了實時監(jiān)測及數(shù)據(jù)分析。

近年來，隨著高精度雷管、數(shù)碼電子雷管等新型爆破器材的問世，通過微差爆破技術(shù)可改善爆破效率、提高爆破質(zhì)量，同時可有效減少爆破地震波、沖擊波等的危害效應(yīng)[1-3，12]。

貴黃高速公路TJ13 標(biāo)段路基邊坡，爆區(qū)之間有一既有高鐵橋梁需要保護。為確保爆破安全，需結(jié)合鄰近保護目標(biāo)的振動危害效應(yīng)監(jiān)測信號及數(shù)值模擬分析，開展微差爆破方案優(yōu)化研究。

1 現(xiàn)場試驗設(shè)計

貴黃高速公路TJ13 標(biāo)段路基邊坡長260 m，爆區(qū)寬20 m，開挖高度16～20 m，中硬巖，總石方量約20 000 m3。爆區(qū)I和爆區(qū)II之間的既有高鐵橋梁需被重點保護。其中，七號墩距爆區(qū)最近，離爆區(qū)I北面340 m，爆區(qū)II南面582 m，如圖1 所示。

試驗方案設(shè)計如下:

1) 爆區(qū)I距離既有橋梁更近?，F(xiàn)場試驗時，采用TC-4850 爆破測振儀，在爆區(qū)I微差爆破時對鄰近高鐵橋梁的振動響應(yīng)進行實時數(shù)據(jù)采集。

2) 將爆區(qū)I的邊坡開挖成高1.5 m、寬1.5 m的臺階，在邊坡中部鉆3 個炮孔，孔徑100 mm、孔深3.0 m、孔距3.0 m，單耗取0.35 kg/m3，逐孔起爆，孔間微差時間初步設(shè)定為6 ms。

3) 測點布置在七號墩以及離爆源中心直線距離最短的第六跨梁，測點編號如圖2 所示。迎爆側(cè)，從七號墩底開始間隔2 m 布置1＃～5＃測點，從第六跨梁跨中梁體開始間隔2 m 布置6＃～10＃測點，橋墩上1＃測點、梁體上7＃測點距爆源最近。

2 數(shù)值模型建立

2.1 計算模型

基于ANSYS/ LS-DYNA軟件，建立爆區(qū)I邊坡和鄰近既有橋梁的1∶1 數(shù)值計算模型。采用ALE算法，Solid164 單元，長340 m，寬260 m，按實際尺寸建立3 個炮孔，微差起爆時間6 ms，總計算時長200 ms。橋梁結(jié)構(gòu)建模時進行簡化，跨徑30 m，梁和墩模型材料分別為C55 和C50 混凝土?？傮w建模如圖3 所示。

2.2 材料模型及參數(shù)

巖石采用彈塑性動力模型(Mat_Plastic_Kinematic)。該模型能精確地模擬巖石材料的動力特性，其屈服條件為[13]:

式中:sij為柯西應(yīng)力張量；aij為屈服面中心點處的應(yīng)力偏張量；F和C為輸入常數(shù)；σy為屈服應(yīng)力；σ0和β分別為初始屈服應(yīng)力和硬化參數(shù)；為應(yīng)變率；Ef為塑性硬化模量，Ef=EtE/(E -Et)；E為彈性模量；Et為切線模量；為有效塑性應(yīng)變。

根據(jù)地質(zhì)資料，巖石力學(xué)參數(shù)如表1 所示。

表1 巖石參數(shù)Tab.1 Parameters of rock

梁和墩分別采用標(biāo)號為C55 和C50 的混凝土，混凝土材料參數(shù)見表2。

表2 混凝土參數(shù)Tab.2 Parameters of concrete

炸藥采用高能炸藥材料模型(Mat_High_Explosive_Burn)，結(jié)合JWL 狀態(tài)方程[13]，對炸藥起爆點、起爆時間加以控制，實現(xiàn)對炸藥爆炸過程的模擬。

式中:p為爆轟壓力；V是相對體積；E是單位體積內(nèi)能；ω、A、B、R1、R2為材料常數(shù)。

具體參數(shù)如表3 所示。

表3 炸藥材料及狀態(tài)方程參數(shù)Tab.3 Parameters of explosive material and equation of state

2.3 邊界條件

在有限元模型中，需要采用有限域來模擬爆破地震波作用下炸藥爆炸產(chǎn)生的鄰近高架橋的動力響應(yīng)規(guī)律。在ANSYS 建模過程中，需要在范圍有限的模型邊界上施加無反射邊界條件，消除爆破地震波傳遞到模型邊界時發(fā)生的反射、折射現(xiàn)象，確保計算結(jié)果的準(zhǔn)確性。在計算模型中，除上部邊界面外，其余模型面均施加無反射邊界。

2.4 數(shù)值計算方案驗證

基于上述試驗和數(shù)值模擬方案，將逐孔起爆微差時間6 ms 時，實際采集到的橋墩和跨梁上各測點的爆破振速峰值與數(shù)值模擬計算所得各測點的振速峰值進行對比，如圖4 所示?？傮w吻合情況良好，最大偏差不超過10%，驗證了上述所采用的數(shù)值模擬方案基本科學(xué)合理。

3 數(shù)值模擬分析

采用控制變量法進行數(shù)值模擬及對比分析。保持測點位置不變，研究微差時間對鄰近橋梁同一測點處振速響應(yīng)的影響規(guī)律；保持微差時間不變，研究測點所處空間位置對測點處振速響應(yīng)的影響規(guī)律。

在上述數(shù)值模擬方案基礎(chǔ)上，增加樣本量以提升數(shù)值模擬及數(shù)據(jù)分析的可靠性。分別增設(shè)0、5、10、15、20、25 ms 共6 組逐孔起爆微差時間，并在六號墩相對于七號墩同樣位置處增設(shè)11?！?5＃共5 個測點，在第六跨梁遠(yuǎn)離爆源的另半跨梁體上對稱增設(shè)16?！?0＃共5 個測點，共計20 個測點，測點編號及位置如圖5 所示。分別計算得到20 個測點在6組微差時間下的爆破振速時程曲線。提取這120 組振速時程曲線的峰值，進行不同微差起爆時差條件下鄰近高鐵橋梁的爆破振動響應(yīng)分析。

3.1 微差時間對橋梁振速響應(yīng)的影響規(guī)律

將1＃測點分別在6組微差時間下的振速響應(yīng)峰值繪成散點形式，并用樣條曲線連接。對2＃～10＃測點進行同樣操作，結(jié)果如圖6(a)所示；對11?！? 0＃測點進行同樣操作，結(jié)果如圖6(b)所示。10條不同顏色的曲線，分別代表10 個不同測點處的綜合振速峰值隨不同微差時間的變化規(guī)律。

由圖6 分析可知:

1) 1?！?0＃測點和11＃～20＃測點分布規(guī)律相近，故振速響應(yīng)隨微差時間的變化規(guī)律大致相同。

2) 綜合測點處振速峰值隨微差時間變化規(guī)律有兩種類型:隨微差時間的增加而快速下降，隨后微差時間繼續(xù)增加時趨于平緩，成波浪狀微微起伏，此規(guī)律對應(yīng)的是1＃、3＃、4＃、6＃、10＃測點和11＃、13＃、14＃、16＃、20＃測點，其中，1＃、3＃、4＃、11＃、13＃、14＃測點分布在橋墩上，6＃、10＃、16＃、20＃測點分布在跨梁上；沒有隨微差時間增加而快速下降的階段，在微差時間增加的全程均緩慢上升或下降，成波浪狀微微起伏，對應(yīng)的是2＃、5＃、7＃、8＃、9＃測點和12＃、15＃、17＃、18＃、19＃測點，其中，2＃、5＃、12＃、15＃測點分布在橋墩上，7＃、8＃、9＃、17＃、18＃、19＃測點分布在跨梁上。

進一步分析可知:

1) 微差時間對鄰近橋梁的振速響應(yīng)有顯著影響，但只集中在一個區(qū)間范圍內(nèi)，超出此區(qū)間，則微差時間對振速響應(yīng)影響不顯著。炮孔之間微差時間若過小，則接近于同時起爆；若過大，則接近于逐孔起爆。數(shù)值模擬結(jié)果及試驗測試數(shù)據(jù)分析表明，炮孔微差時間很小或很大情況下，改變微差時間來控制爆破振動響應(yīng)的效果均不明顯。本工程背景下對鄰近橋梁振速響應(yīng)有顯著影響的微差時間取值范圍是0～5 ms，而5～10 ms 則是一個下降到上升的過渡階段，這里存在一個極小值。

2) 針對不同位置處的測點，對其振速峰值有顯著影響的微差時間區(qū)間范圍不同。

3)隨著微差時間的增加，結(jié)構(gòu)振速響應(yīng)的變化并非成單調(diào)遞增或遞減趨勢，而是反復(fù)增減，故總存在一個振速響應(yīng)的極小值。

3.2 空間位置對橋梁振速響應(yīng)的影響規(guī)律

以微差時間為控制變量，研究測點所處空間位置對測點處振速響應(yīng)的影響規(guī)律。列出在6 組微差時間下鄰近橋梁七號墩上1＃～5＃測點振速響應(yīng)的空間分布規(guī)律，如圖7(a)所示；六號墩上11＃～15＃測點振速響應(yīng)的空間分布規(guī)律如圖7(b)所示。第六跨梁上6?！?0＃測點振速響應(yīng)的空間分布規(guī)律如圖8(a)所示；16?！?0＃測點振速響應(yīng)的空間分布規(guī)律如圖8(b)所示。7 種不同顏色的曲線中，前6 條分別代表在6 組不同的微差時間下各測點的振速分布規(guī)律，第7 條曲線代表地表振速峰值隨爆心距呈指數(shù)衰減的規(guī)律曲線[14]:

式中:v為爆破振動質(zhì)點最大振速，cm/s；Q為炸藥量，kg，齊爆時取總藥量，延時起爆時取單段最大藥量；R為爆源的距離，m；d為比例距離；K和α分別為爆源至測點間地形、地質(zhì)條件的相關(guān)系數(shù)和衰減指數(shù)。

由數(shù)據(jù)分析可知:

1) 從第7 條曲線分析可知，1?！?＃測點的振速響應(yīng)隨微差時間的變化規(guī)律大致和11＃～15＃測點的振速響應(yīng)隨微差時間的變化規(guī)律相同；6?！?0＃測點的振速響應(yīng)隨微差時間的變化規(guī)律大致和16＃～20＃測點的振速響應(yīng)隨微差時間的變化規(guī)律相同；再次印證了所采用數(shù)值模擬方案的正確性。

2) 跨梁上測點振速起初隨測點離爆源的距離增加而減小，但最遠(yuǎn)的靠近橋墩頂處的跨梁上測點的振速反而增大，存在一定程度的高程放大效應(yīng)。

由此可知，鄰近高鐵橋梁受爆破作用后的振速響應(yīng)不僅和爆破作用大小有關(guān)，也和結(jié)構(gòu)所受約束條件有關(guān)，在橋墩底部和主梁靠近墩頂處這些端部固結(jié)位置，會出現(xiàn)離爆源較遠(yuǎn)的地方振速響應(yīng)反而增大的情況，在爆破安全校核過程中應(yīng)予以重視，確保爆破安全。

4 結(jié)論

1) 微差時間對路基邊坡爆破時鄰近橋梁的振速響應(yīng)有明顯影響，規(guī)律體現(xiàn)為非線性、非單調(diào)性。微差時間對鄰近橋梁的作用效應(yīng)不能線性疊加，在有限區(qū)間內(nèi)存在一個最佳的微差時間，使鄰近橋梁的爆破振動響應(yīng)最小，并可通過現(xiàn)場爆破振動監(jiān)測及數(shù)值模擬方法進行確定。

2) 對于橋梁的墩柱部分，爆破振速響應(yīng)與爆源的距離成反比。但橋梁結(jié)構(gòu)的梁體部分，在靠近墩柱處的振速則并非隨離爆源距離的增加而遞減，爆破安全校核環(huán)節(jié)應(yīng)引起重視。