單向靜載下柱狀藥包爆破裂紋擴(kuò)展規(guī)律及機(jī)制

彭建宇，李元輝，張鳳鵬，范光華

(東北大學(xué) 深部金屬礦山安全開采教育部重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110819)

單向靜載下柱狀藥包爆破裂紋擴(kuò)展規(guī)律及機(jī)制

彭建宇，李元輝，張鳳鵬，范光華

(東北大學(xué) 深部金屬礦山安全開采教育部重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110819)

將巖體內(nèi)柱狀藥包爆破簡(jiǎn)化為平面應(yīng)變問(wèn)題，采用ANSYS/LS-DYNA軟件進(jìn)行數(shù)值計(jì)算，研究單向靜應(yīng)力對(duì)爆破裂紋的影響規(guī)律及機(jī)制。結(jié)果表明：無(wú)靜應(yīng)力時(shí)，爆破裂紋區(qū)形狀大致為圓形；而施加單向靜應(yīng)力后，裂紋區(qū)形狀演化成橢圓形，橢圓的長(zhǎng)軸方向與靜應(yīng)力方向相同，且隨著靜載荷增加，裂紋區(qū)長(zhǎng)軸基本保持不變，短軸長(zhǎng)度降低。這是由于隨靜載荷變化，沿著加載方向單元的最大主應(yīng)變基本不變；而垂直加載方向單元的最大主應(yīng)變隨著靜載荷增加而減小，從而導(dǎo)致該方向裂紋長(zhǎng)度降低。

單向靜應(yīng)力；巖石爆破；ANSYS/LS-DYNA；裂紋擴(kuò)展

隨著淺部資源的消耗，越來(lái)越多的礦山進(jìn)入深部開采狀態(tài)[1]。深部礦巖在爆破開挖過(guò)程中不僅受爆炸產(chǎn)生的沖擊載荷的作用，還會(huì)受到地應(yīng)力的影響，巖石的破壞是初始地應(yīng)力和爆炸載荷共同作用的結(jié)果。因此，開展地應(yīng)力下的爆破裂紋擴(kuò)展規(guī)律研究具有重要意義。

國(guó)內(nèi)外學(xué)者采用理論分析、室內(nèi)試驗(yàn)及數(shù)值模擬等多種手段研究初始應(yīng)力下巖石的爆破問(wèn)題。理論方面，戴俊等[2]以柱裝藥爆破漏斗理論為基礎(chǔ)，對(duì)高地應(yīng)力巷道崩落爆破標(biāo)準(zhǔn)漏斗的形成進(jìn)行了研究；陳明等[3]采用摩爾庫(kù)侖準(zhǔn)則及最大拉應(yīng)力準(zhǔn)則，研究了地應(yīng)力對(duì)爆生裂隙區(qū)比例半徑的影響。實(shí)驗(yàn)研究方面，肖正學(xué)等[4]的實(shí)驗(yàn)結(jié)果表明初始應(yīng)力場(chǎng)能改變爆破應(yīng)力波的傳播規(guī)律，對(duì)裂紋擴(kuò)展具有導(dǎo)向作用。楊立云等[5]設(shè)計(jì)了動(dòng)靜組合加載系統(tǒng)，結(jié)合焦散線實(shí)驗(yàn)手段，研究靜應(yīng)力場(chǎng)對(duì)爆生裂紋擴(kuò)展規(guī)律的影響。徐穎等[6]開展了高地應(yīng)力條件下深部圍巖爆破開挖三維相似模型試驗(yàn)。此外，由于數(shù)值模擬方法具有經(jīng)濟(jì)、快速、易獲取數(shù)據(jù)指標(biāo)等優(yōu)勢(shì)。Doze F. V., et al[7]用數(shù)值模擬的手段研究靜應(yīng)力下的爆破過(guò)程，結(jié)果顯示靜載荷方向?yàn)楸屏鸭y擴(kuò)展的主導(dǎo)方向。楊棟等[8]采用有限差分程序FLAC3D研究了爆破荷載與地應(yīng)力動(dòng)態(tài)卸荷復(fù)合作用下隧道圍巖損傷分布。魏晨慧等[9]模擬了地應(yīng)力條件下切縫藥包爆破裂紋擴(kuò)展過(guò)程，結(jié)果表明裂紋的擴(kuò)展方向受切縫角度和最大地應(yīng)力方向兩個(gè)條件影響。然而，已有的靜應(yīng)力下爆破過(guò)程模擬工作多局限在裂紋擴(kuò)展形態(tài)方面，缺少相關(guān)機(jī)理分析。本文采用有限元軟件ANSYS/LS-DYNA[10-11]模擬單向靜載荷下巖石爆破效果，并提取模型關(guān)鍵位置應(yīng)力、應(yīng)變，研究靜載荷下巖石的動(dòng)態(tài)破壞機(jī)制，為深部巖體工程開挖設(shè)計(jì)提供一定的參考。

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

1.1 數(shù)值計(jì)算模型

圖1 單軸靜應(yīng)力下柱狀藥包爆破模型

柱狀藥包爆破時(shí)，藥包長(zhǎng)度遠(yuǎn)大于藥包直徑，長(zhǎng)度方向各截面的破壞情況類似，可看作平面應(yīng)變模型。本文只取垂直于炸藥長(zhǎng)度方向的一個(gè)截面進(jìn)行建模，計(jì)算模型如圖1所示。模型尺寸為300mm×300mm，藥卷直徑為6mm，位于模型幾何中心。模型使用單層網(wǎng)格，在網(wǎng)格厚度方向施加位移約束。本文僅開展單向靜載荷下的爆破模擬研究，為了便于分析，建立圖1所示的坐標(biāo)系。取炮孔中心為坐標(biāo)原點(diǎn)，靜應(yīng)力加載方向?yàn)閤方向，無(wú)靜載荷方向?yàn)閥方向。施加的初始?jí)簯?yīng)力范圍為0～60 MPa。

1.2 巖體材料模型及參數(shù)

(1)

(2)

(3)

(4)

巖石隨動(dòng)硬化的Mises屈服條件滿足式(5)。

(5)

數(shù)值模擬時(shí)，采用MAT_ADD_EROSION關(guān)鍵字定義巖體破壞準(zhǔn)則，該關(guān)鍵字可以確定多種破壞標(biāo)準(zhǔn)。當(dāng)巖石材料單元達(dá)到關(guān)鍵字中設(shè)定的強(qiáng)度準(zhǔn)則就會(huì)失效刪除，在巖體中出現(xiàn)裂縫，模擬爆破裂紋形成過(guò)程。模擬采用的巖石力學(xué)參數(shù)見表1。

1.3 炸藥材料參數(shù)

LS-DYNA程序中，用JWL狀態(tài)方程描述炸藥爆炸產(chǎn)生沖擊壓力的過(guò)程，JWL方程中壓力和比容的關(guān)系見式(6)。

(6)

式中：A、B、R1、R2、ω為炸藥JWL狀態(tài)方程參數(shù)；PCJ為初始?jí)毫Γ籚為相對(duì)體積；E0為初始比內(nèi)能。

模擬采用的炸藥參數(shù)見表2。

表1 巖石材料參數(shù)

表2 炸藥材料參數(shù)

2 計(jì)算結(jié)果分析

2.1 單向靜應(yīng)力下爆破破壞范圍

圖2為不同靜應(yīng)力條件下巖石爆破后的最終破壞形態(tài)?？梢钥闯銎茐膮^(qū)域由粉碎區(qū)及裂隙區(qū)兩部分組成。爆破形成的粉碎區(qū)大小受靜應(yīng)力影響較小，基本不發(fā)生變化。但由圖2(a)可以看出，無(wú)靜應(yīng)力時(shí)，巖石破壞裂隙區(qū)對(duì)稱分布，大致為圓形；施加單向靜應(yīng)力后，裂隙區(qū)變?yōu)闄E圓形(圖2(b)～(d))。

圖3為爆破后裂紋長(zhǎng)度隨靜應(yīng)力的變化曲線。隨著靜應(yīng)力的增加，x方向裂紋長(zhǎng)度變化不大，而y方向裂紋長(zhǎng)度逐漸減小。因此，靜應(yīng)力不影響自身方向的裂隙區(qū)大小，但對(duì)垂直于初始應(yīng)力方向上的裂紋形成起到明顯的抑制作用。

圖2 不同靜應(yīng)力下的爆破破壞區(qū)域

圖3 裂紋長(zhǎng)度隨靜應(yīng)力的變化

2.2 初始應(yīng)力對(duì)爆破效果的影響機(jī)制

圖4為不同靜應(yīng)力條件下炮孔附近同一單元的壓力時(shí)程曲線?？梢钥闯?，爆炸近區(qū)單元的應(yīng)力很大，不同靜應(yīng)力下爆破后單元的壓力相差很小，曲線基本重合在一起。說(shuō)明在近區(qū)爆炸沖擊波占主導(dǎo)作用，靜應(yīng)力對(duì)爆炸應(yīng)力波影響較小。圖5不同靜應(yīng)力條件下爆破過(guò)程中遠(yuǎn)區(qū)某一單元的壓力時(shí)程曲線。在爆炸遠(yuǎn)區(qū)，爆炸應(yīng)力峰值較爆破近區(qū)衰減明顯，不同靜應(yīng)力下單元壓力相差明顯，說(shuō)明靜應(yīng)力對(duì)爆破遠(yuǎn)區(qū)應(yīng)力場(chǎng)影響明顯。此外，由圖4及圖5還可以看出，初始應(yīng)力不影響應(yīng)力波的到達(dá)時(shí)間。

圖4 爆炸近區(qū)單元壓力時(shí)程曲線

圖5 爆炸遠(yuǎn)區(qū)單元壓力時(shí)程曲線

以炮孔中心為原點(diǎn)，取模型上與x方向成45°角的直線L1作為分析路徑(圖1)。L1上單元的峰值壓力和峰值有效應(yīng)力隨距離衰減曲線分別如圖6及圖7所示。由圖6、圖7可以看出，無(wú)靜應(yīng)力和存在x方向靜應(yīng)力兩種條件下的爆破，單元峰值壓力和峰值有效應(yīng)力衰減規(guī)律基本一致，說(shuō)明靜應(yīng)力不改變爆炸應(yīng)力波傳播過(guò)程中的衰減規(guī)律。

取模型上半徑為140mm的1/4圓為分析線L2，選取L2上任意一點(diǎn)與炮孔中心連線為直線L，L與y軸的夾角為θ(圖1)。不同靜應(yīng)力下分析線L2上單元最大主應(yīng)變隨角度θ的變化曲線如圖8所示。從圖8中可以看出，沒(méi)有靜應(yīng)力時(shí)，L2上各方向最大主應(yīng)變數(shù)值較為接近，且0°和90°方向上單元的最大主應(yīng)變大小相同，因此巖石破壞區(qū)近似為圓形。存在靜應(yīng)力時(shí)，隨著θ變大，單元最大主應(yīng)變逐漸增大；靜應(yīng)力越大，0°方向(垂直靜應(yīng)力方向)單元最大主應(yīng)變?cè)叫。?0°(平行靜應(yīng)力方向)單元的最大主應(yīng)變基本不變。因此，x方向靜應(yīng)力越大，y方向的裂紋受到抑制作用越明顯，而x方向的裂紋擴(kuò)展基本不受影響。這很好地解釋了圖3裂紋長(zhǎng)度隨靜應(yīng)力的變化規(guī)律。

圖6 L1上單元峰值壓力衰減曲線

圖7 L1上單元峰值有效應(yīng)力衰減曲線

圖8 L2上單元最大主應(yīng)變隨角度變化曲線

3 結(jié)論

1)靜應(yīng)力改變巖石破壞區(qū)域形態(tài)。無(wú)靜應(yīng)力作用下，巖石在柱狀藥包作用下的裂紋區(qū)域近似為圓形。施加單向靜應(yīng)力后，巖石破壞區(qū)域近似為橢圓形，且靜應(yīng)力方向?yàn)閹r體爆破裂紋擴(kuò)展的主方向；隨著靜應(yīng)力增加，破壞區(qū)短軸長(zhǎng)度減小，長(zhǎng)軸長(zhǎng)度相對(duì)變化不大。

2)靜應(yīng)力不影響巖體內(nèi)應(yīng)力波衰減規(guī)律，但是能改變巖體不同位置的應(yīng)變分布。無(wú)論是否存在初始靜應(yīng)力，爆破過(guò)程中單元峰值壓力和有效應(yīng)力均以指數(shù)形式衰減。此外，靜應(yīng)力方向的最大主應(yīng)變基本不隨靜載荷變化而變化，垂直靜應(yīng)力方向的最大主應(yīng)變隨著靜載荷增加而減小。因此初始靜載荷能夠改變爆破載荷下巖體內(nèi)的應(yīng)變分布規(guī)律，進(jìn)而影響裂紋擴(kuò)展過(guò)程。

[1] 劉東，郭進(jìn)平. 深部高應(yīng)力條件下無(wú)底柱分段崩落采礦法的應(yīng)用[J]. 中國(guó)礦業(yè)，2015，24(12): 109-111.

[2] 戴俊，錢七虎. 高地應(yīng)力條件下的巷道崩落爆破參數(shù)[J]. 爆炸與沖擊，2007，27(3): 272-277.

[3] 陳明，盧文波，周創(chuàng)兵，等. 初始地應(yīng)力對(duì)隧洞開挖爆生裂隙區(qū)的影響研究[J]. 巖土力學(xué)，2009，30(8): 2254-2258.

[4] 肖正學(xué)，張志呈，李端明. 初始應(yīng)力場(chǎng)對(duì)爆破效果的影響[J]. 煤炭學(xué)報(bào)，1996，21(5): 497-501.

[5] 楊立云，楊仁樹，許鵬，等. 初始?jí)簯?yīng)力場(chǎng)對(duì)爆生裂紋行為演化效應(yīng)的實(shí)驗(yàn)研究[J]. 煤炭學(xué)報(bào)，2013，38(3): 404-410.

[6] 徐穎，袁璞. 爆炸荷載下深部圍巖分區(qū)破裂模型試驗(yàn)研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2015，34: 3844-3851.

[7] Donze F V, Bouchez J, Magnier S A. Modeling fractures in rock blasting[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(8)： 1153-1163.

[8] 魏晨慧，朱萬(wàn)成，白羽，等. 不同地應(yīng)力條件下切縫藥包爆破的數(shù)值模擬[J]. 爆炸與沖擊，2016，36(2): 161-169.

[9] 楊棟，李海波，夏祥，等. 高地應(yīng)力條件下爆破開挖誘發(fā)圍巖損傷的特性研究[J]. 巖土力學(xué)，2014，35(4): 1110-1122.

[10] LSTC. LS-DYNA Keyword User's Manual[M]. California: Livemore Software Technology Corporation, 2001.

[11] 葉海旺，王亮，張迎吉，等. 金堆城南露天數(shù)碼雷管減震爆破試驗(yàn)研究[J]. 中國(guó)礦業(yè)，2014，23(12): 136-139.

Crack propagation rule and mechanism of cylindrical blasting under uniaxial static load

PENG Jianyu，LI Yuanhui，ZHANG Fengpeng，F(xiàn)AN Guanghua

(Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines，Northeastern University，Shenyang 110819，China)

The cylindrical blasting in rock mass was simplified as plane strain problems and the effect of uniaxial static stress on blasting crack was studied by ANSYS/LS-DYNA software. The results show that the shape of the crack area is roughly circular when no static stress is applied, while the crack area is elliptic after the static stress is applied. The long axis direction of the ellipse is the same as the static stress direction. As the static load increases, the long axis of the crack region remains unchanged, and the length of the minor axis decreases. This is because the maximum principal strain along the loading direction is basically unchanged with increasing static load; and the maximum principal strain of element in the vertical load direction decreases, resulting in the crack length decreases.

uniaxial static stress; rock blasting; ANSYS/LS-DYNA; crack propagation

2016-10-18

國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目資助(編號(hào)：2016YFC0801605；2016YFC0600706)；國(guó)家自然科學(xué)基金項(xiàng)目資助(編號(hào)：51674061)

彭建宇(1987-)，男，博士。

O383.1；TD853

A

1004-4051(2017)01-0088-04