邊坡爆破對既有高邊坡的振動影響分析

2018-05-08 11:43:37■王華

福建交通科技 2018年2期

■王華

（福建省交通規劃設計院，福州 350004）

1 引言

隨著公路工程不斷建設發展，新建邊坡臨近既有邊坡的情況越來越普遍，因此分析新建邊坡建設的影響并提出合理的應對措施，是此類邊坡工程建設過程中亟需解決的問題。已有文獻更多地集中在研究邊坡開挖后的長期穩定性問題[1,2]。但是，新建邊坡爆破施工易造成既有邊坡巖體原有的節理、裂隙進一步擴展，嚴重時還可能造成巖體沿已有結構面滑移，直接影響邊坡的整體穩定性[3,4]。因此,分析爆破振動效應對邊坡的影響顯得尤為重要。國內有研究者利用有限元法對爆破荷載作用進行數值模擬分析，得到了較為準確的結果，這說明了利用數值模擬方法解決該類問題的可靠性[5,6]。

本文以福建省漳州市某公路工程新建邊坡作為研究對象，對臨近高速公路匝道高邊坡段落范圍內爆破施工過程進行三維動力數值模擬分析，確定了爆破荷載在不同位置對既有高邊坡的振動影響，制定了安全合理的爆破方案，并對同類工程具有指導意義。

2 爆破振動效應判據標準

我國《爆破安全規程》（GB6722-2014）[9]中，提出了永久性匝道高邊坡各頻率范圍內的安全允許爆破振速如表1所示。根據依托工程情況，采用淺孔爆破施工方案，因此取振動速度的控制值為10cm/s。

3 爆破振動影響數值模擬分析

3.1 工程概況

福建省漳州市某公路工程新建邊坡工程與既有高速公路匝道邊坡相鄰，其交叉位置相互關系具體如圖1所示。

表1 爆破振動安全允許標準

新建公路邊坡最大邊坡高度為29.6m，共4階，爆破開挖區域基巖主要以中風化花崗巖和散體狀強風化花崗巖為主。既有匝道高邊坡的最大邊坡高度為35.1m，共4階，以散體狀強風化花崗巖～碎塊狀強風化花崗巖為主，采用M7.5漿砌片石拱型骨架防護，典型橫斷面如圖2所示。

由于依托工程新建邊坡與既有邊坡的空間位置關系復雜，且爆破振動荷載作用過程對邊坡動力響應有較大影響，因此本文采用基于時程分析的三維動力有限元法進行數值模擬。

3.2 爆破計算模型及參數

基于MIDASGTS NX有限元分析軟件，建立三維有限元模型，并施加爆破時程荷載，分析既有匝道高邊坡的振速和應力響應。通過調整爆破位置，從而確定滿足安全爆破振動速度的爆破距離。

圖1 交叉位置相互關系

圖2 既有匝道高邊坡橫斷面圖

計算模型縱向長313m，橫向長598m，高度方向長60m，采用三維有限元模型，共劃分23468個單元，31810個節點，巖土采用十結點四面體高階實體單元模擬，如圖3～4所示。

圖3 邊坡開挖前計算模型

圖4 邊坡開挖后計算模型

3.2.1 本構模型與邊界條件

巖土具有復雜的變形特征，且隨加載條件的變化而發生變化，本構模型雖然有很多種，但實際上并不存在一種本構模型能模擬實際土體的所有復雜變形特征。一個合理的本構模型應該滿足參數易確定性，理論上的嚴格性和計算機實現的可能性。本文研究的邊坡巖體主要由全風化花崗巖和散體狀強風化花崗巖組成，可視為各向同性體，計算時采用Mohr-Coulomb強度準則。

由于巖體的動彈性模量Ed都普遍大于靜彈性模量Eme，可利用巖體完整性系數Kv得到動彈性模量和靜彈性模量的換算關系。按地質縱斷圖中的描述，取Kv為0.35～0.15，取j為0.2。因此，選定爆破動力分析計算參數如表2所示。

表2 爆破動力分析計算參數

對于動力分析，建立一般的邊界條件會由于波的反射作用而產生很大的誤差。計算時采用粘彈性邊界既可以約束動力問題中的零頻分量，能夠模擬人工邊界外半無限介質的彈性恢復性能，具有良好的穩定性和較高的精度。因此邊界條件中除地表為自由邊界外，其余邊界面均施加了法向曲面彈簧。

3.2.2 爆破等效荷載

淺孔爆破施工方案根據開挖深度分一個或兩個臺階進行爆破，邊坡采用光面爆破。本次計算擬定炮孔直徑45mm，深度H=4m，底板抵抗線Wp=1.45m，超鉆深度h=0.2m，孔距 a=2m，排距 b=2m；炸藥單耗 q=0.45kg/m3，炮孔方向：中間主炮孔取垂直孔，邊坡光面孔與邊坡坡率相同。考慮到爆破振速跟炸藥量密切相關，為控制振速，單次爆破按縱向小臺階爆破，單次縱向爆破4m，共3個爆孔，單孔裝藥量5kg，共15kg。

爆破荷載采用美國National Highway Institute里提及的公式，每1kg炸藥的爆破壓力如下：

式中：Pdet為爆破壓力，ρ為炸藥比重（取為1g/cm3），Ve為爆破速度（取為300000cm/s），PB為作用在孔壁面上最大爆破壓力，dc為火藥直徑，dh為孔眼直徑。

根據公式（2）求得每1kg炸藥的最大爆破壓力PB為20.9MPa。為控制振速，考慮施工時按縱向小臺階爆破，單次縱向爆破4m，共3個爆孔，單孔裝藥量5kg，共15kg。本次計算取單孔爆破壓力為104.5MPa，考慮三個孔同時起爆。

由于爆破壓力是隨時間變化的，計算中采用以下動壓力時程公式：

式中：B=16338，為荷載常量。

3.2.3 計算工況

為了分析在邊坡不同位置起爆對既有高邊坡的影響，本文設計了6種工況，各工況的起爆位置如圖5所示。表3列出各工況爆破點距新建邊坡與既有邊坡交點的距離。

圖5 爆破荷載施加位置

表3 爆破荷載施加位置

3.3 數值模擬結果與分析

3.3.1 爆破振動速度結果

計算得到的工況1（爆破點距離既有匝道高邊坡最近）既有邊坡的最大爆破振動合速度為155.7cm/s，位于爆破點附近新建邊坡和既有匝道高邊坡的交點處，發生于0.043s時刻。

繪制出不同工況下既有匝道高邊坡最大振動速度隨距離變化曲線分別如圖6、表4所示。經內插可知起爆點距新建邊坡與既有邊坡交點23m以外時，既有邊坡的最大振動速度小于10cm/s，滿足《爆破安全規程》允許振動速度要求。

圖6 各工況既有匝道高邊坡最大振動速度變化曲線

表4 各工況既有匝道高邊坡最大振動速度統計表

3.3.2 應力結果

圖7～8為計算工況1爆破后，既有匝道高邊坡表面監測點的第一主應力（最大拉應力）和最大剪應力結果。當爆破時刻為0.043s時，邊坡表面最大拉應力達到最大值，為0.64MPa，最大剪應力為 0.32MPa，已接近或超過M7.5砂漿的抗拉和抗剪強度0.54MPa，說明在工況1爆破作業下，存在漿砌片石開裂脫落、掉塊的危險。