張鳳鵬，彭建宇，李  瑩，李元輝

（東北大學深部金屬礦山安全開采教育部重點實驗室，遼寧 沈陽  110819）

摘  要：將無限大巖體內柱狀藥包爆破簡化為平面應變問題，基于雙線性隨動硬化屈服理論，采用ANSYS/LS-DYNA對爆破過程進行模擬，研究垂直炮孔方向的雙向地應力場對裂紋擴展規律的影響。數值結果表明：雙向等壓下，裂紋區形狀為圓形，其面積隨地應力增大而非線性減小，但裂紋密度增大，給出了裂紋區面積與壓應力之間的定量關系；而雙向不等壓情況下，裂紋區形狀為近似橢圓形，長軸出現在較大壓應力作用方向，短軸出現在較小壓應力作用方向，在較小壓應力不變的情況下，隨著較大壓應力增大，長軸先增大后減小，短軸減小，裂紋區面積先增加后減小，呈非單調變化。

關鍵詞：數值模擬；地應力；爆破；巖體；ANSYS/LS-DYNA；裂紋擴展

Simulation Study of The Effect of In-situ Stress on Rock Blasting

Zhang Fengpeng, Peng Jianyu, Li Ying, Li Yuanhui

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

Abstract：Cylindrical charge blasting in infinite rock was simplified to plane strain problem, based on bilinear kinematic hardening yield theory, the blasting process was simulated by ANSYS/LS-DYNA under different in-situ stresses vertical the axis of blasthole. Numerical results show that: Fracture zones are circles when the two direction compressive in-situ stresses are equal, cracks area decrease with the increasing of in-situ stress, but the crack density increases, and the quantitative relationship between compressive stress and crack zone area was given; Fracture zones are ellipses when the two direction compressive in-situ stresses are not equal, the long axis of fracture zone coincides with the larger compressive principal stress and the short axis of fracture zone coincides with the smaller compressive principal stress. Under the condition of invariable in smaller compressive stress, the length of long axis and the area of cracks decrease after increasing with increasing of the larger compressive stress, and the length of short axis decreases always.

Key words: numerical simulation; in-situ stress; blasting; Rock mass; ANSYS/LS-DYNA; crack evolution

0  引  言

隨著人們對能源需求量的增加和礦山開采強度的不斷增大，淺部資源日益減少，國內外礦山都相繼進入深部資源的開采狀態。深部巖體處于較高的地應力環境中，在深部巖體中進行開采爆破時，巖體處于高地應力和爆炸沖擊波的共同作用，兩種載荷相互影響，導致巖體爆破破壞過程比較復雜。目前，高應力條件下巖石的爆破日益受到國內外學者的重視^[1-3]，國內外學者從實驗、理論分析和數值模擬等方面對高應力深部巖體爆破過程進行研究。肖正學等^[4]通過對存在不同初始應力場的幾種材料的室內試驗，認為初始應力場的存在可改變爆轟波的傳播規律以及對裂紋發展起著導向作用。謝源^[5]對高應力條件下巖石爆破裂紋擴展規律的模擬進行了試驗研究，研究表明介質爆破裂紋的方向及大小與附加的主應力有關。白羽等^[6]研究了地應力對雙孔爆破的影響。戴俊等^[7]以柱狀裝藥爆破漏斗理論為基礎，對高地應力巷道崩落爆破標準漏斗形成進行研究。劉艷等^[8]通過隱式-顯式連續求解進行了地應力作用下巖體爆炸的數值模擬。

盡管國內外學者對高應力條件下巖石爆破進行了大量研究，但多局限在定性分析層面上，相關定量研究成果仍然比較少。本文將采用ANSYS/LS-DYNA對不同地應力條件下巖石爆破破壞特征進行數值模擬，研究破壞區范圍和靜應力之間的定量關系。

1  計算模型

1.1 炸藥爆炸狀態方程

本文模擬高能炸藥在花崗巖中的爆破過程。炸藥起爆后，體積膨脹，炸藥單元體內一點的壓力采用JWL狀態方程模擬：

式中，p為爆炸壓力（Pa）,F為炸藥化學能釋放率,D為炸藥爆速（m/s）,A_emax和v_e分別為炸藥最大橫截面積和體積,t，t₁分別為當前時間和炸藥內一點的起爆時間（s）,p_eos為炸藥的爆轟壓（Pa）,V為相對體積,E_V為內能參數（Pa）,A、B、R₁、R₂、ω為常數。具體炸藥參數見表1。

1.2 巖體材料模型

本文采用雙線性隨動硬化模型（*MAT_PLASTIC_KINEMATIC）作為巖體的彈塑性屈服模型^[9]，屈服應力與應變率的關系為：

式中，為初始屈服應力（Pa）,為應變率（s^-1）,C和P為Cowper-Symonds應變率參數，取C=2.5和P=4.0,β為硬化參數,E_P為塑性硬化模量（Pa）,E₀為楊氏模量（Pa）,E_tan為切線模量（Pa）,為巖體有效塑性應變，由下式定義：

式中，t為發生塑性應變累計時間（s）,（i, j=1, 2, 3）為巖體塑性應變偏量分量。

爆破過程中粉碎區采用Mises屈服破壞準則，而裂紋區采用拉伸破壞準則^[10]：

式中，為巖體中任一點的Mises有效應力（Pa）,σ_t為爆破產生的拉應力（Pa）；σ_ij（i, j=1, 2, 3）為巖體單元的應力分量；σ_cd和σ_td分別為巖體的單軸動態抗壓強度和抗拉強度（Pa），與靜態強度的關系近似取為：

                                                                                                                                （10）

式中，σ_c和σ_t0分別為巖體的單軸靜態抗壓強度和抗拉強度（Pa）。

本文以花崗巖為研究對象，力學參數見表2。

2  數值計算與分析

2.1 數值計算模型

本文模擬柱狀藥包在無限大巖體內的爆破過程，簡化為平面應變問題進行分析。模擬中取炮孔直徑d=38 mm，計算模型尺寸為2 m×2 m，炸藥位于模型的幾何中心，所建立的計算模型如圖1所示。由于該模型只是無限巖體的一部分，出現了人為的邊界，為了消除人為邊界處的反射波對結構動力響應的影響，在模型四周邊界添加無反射邊界條件，以模擬無限巖體。

2.2 初始應力場的加載方案

在ANSYS中進行隱式分析形成初始應力場，為了對比分析不同附加初始應力情況下巖體爆破特征，設計下列數值分析方案（見表3）：

表中n=s _y/s _x。

2.3 計算結果及分析

（1）雙向等壓地應力場對巖體爆破影響

圖2給出了采用等值雙向加載時不同地應力下的裂紋分布圖，結果表明，當地應力由0增加到20 MPa時，粉碎區的形狀沒有變化，面積略有減小，但變化不顯著。而裂紋在粉碎區外以裝藥中心為圓心呈放射狀分布，裂紋破壞區域為圓形，面積隨地應力增大明顯減小，因此，本文重點討論地應力對裂紋擴展的影響。

圖2  雙向等值地應力作用下的裂紋分布

Fig.2  Crack distribution of the rock model under equal in-situ stresses

表4列出了不同地應力下裂紋區直徑L的具體結果。可見隨著地應力增大，裂紋擴展直徑減小，由地應力為0 MPa時的52.6 d，減小到20 MPa時的22.1 d，但分布規律基本相同。由于裂紋區以環向拉伸破壞為主，隨著與爆腔距離的增加，環向拉應力逐漸減小，此時壓地應力抑制了該區域裂紋的擴展，導致裂紋區面積減小，裂紋密度增加，使該區域巖體破碎更充分，而裂紋區外巖體更完整，減小了爆破損傷區范圍。

圖3給出了裂紋區面積隨地應力變化的規律，其中縱軸A為每種工況下裂紋區的面積， ，采用炮孔橫截面積A₀進行無量綱化，；橫軸為地應力，s=s_x=s_y，采用巖石的抗拉強度s_t0（見表2）進行無量綱化。可見，隨著地應力增加，裂紋區的面積減小速度明顯放緩，有逐漸收斂的趨勢。將圖中曲線擬合為：                            

            （11）

也就是說對于本文研究的花崗巖，在使用文中炸藥爆破時，裂紋區面積與地應力存在式（11）的量化關系。當地應力不斷增大時，裂紋區面積收斂到402.2倍炮孔橫截面積，裂紋區半徑約為20倍炮孔半徑。

（2）雙向不等壓地應力對巖體爆破影響

實際工程中，通常垂直炮孔平面內的兩個主應力是不等的，為此，本文進一步探討雙向不等壓地應力對爆破破巖效果的影響。圖4給出了s_x=10 MPa情況下，s_y分別取10、20、30和40 MPa四種工況下的裂紋分布。結果表明：在雙向不等壓地應力作用下，不但裂紋區面積發生變化，形狀也發生明顯改變，由等壓時的圓形轉化為不等壓工況下的近似橢圓形，長軸出現在較大壓應力作用方向，長度為L_a，短軸出現在較小壓應力作用方向，長度為L_b，隨較大壓應力增加，長軸先增大后減小，短軸減小（見圖5），表明最大壓應力方向是破巖的主導方向。

圖6給出了裂紋區面積隨較大主應力的變化規律，圖中可見在雙向不等壓地應力作用下，裂紋區面積變化過程比較復雜，在一個主應力不變的情況下，隨另一主應力增加而先增加，后減小，呈非單調變化。

3  結  論

通過對不同地應力場下巖體爆破破壞過程進行模擬研究，得出如下結論：

（1）雙向等壓下，裂紋區形狀為圓形，其面積隨地應力增大而減小，但裂紋密度增大。對于無限大花崗巖巖體內柱狀藥包爆破問題，垂直藥柱平面內巖體破壞區面積與壓應力之間的定量關系為：

（2）雙向不等壓情況下，裂紋區形狀為近似橢圓形，長軸出現在較大壓應力作用方向，短軸出現在較小壓應力作用方向，在較小壓應力不變的情況下，隨較大壓應力增加，長軸先增大

后減小，短軸單調減小，裂紋區面積先增加，后減小，最大壓應力方向是巖體主破裂方向。

（3）地應力使爆破裂紋分布更集中，局部裂紋密度增大，使該區域巖體破碎更充分，而裂紋區外巖體更完整，減小了爆破損傷區范圍。

參考文獻:

[1]  Feng X T, Hudson J.Rock engineering design[M].Leiden: CRC Press/Balkema, 2011.

[2]  Petr Konicek, Kamil Soucek, Lubomir Stas, et al.Long-hole destress blasting for rockburst control during deep underground coal mining[J].International Journal of Rock Mechanics & Mining Sciences, 2013, 61(07): 141-153.

[3]  D.H.Diering.Tunnels under pressure in an ultra-deep Witwatersrand gold mine [J].The Journal of the South African Institute of Mining and Metallurgy, 2000, 100(09): 319-324.

[4]  肖正學,張志呈,李端明.初始應力場對爆破效果的影響[J].煤炭學報,1996,21(05):497-501.

[5]  謝源.高應力條件下巖石爆破裂紋擴展規律的模擬研究[J].湖南有色金屬,2002,18(08):1-3.

[6]  白羽,朱萬成,魏晨慧,等.不同地應力條件下雙孔爆破的數值模擬[J].巖土力學,2013,34(S1):466-471.

[7]  戴俊,錢七虎.高地應力條件下的巷道崩落爆破參數[J],爆炸與沖擊,2007,27(03):272-277.

[8]  劉艷,許金余.地應力場下巖體爆體的數值模擬[J].巖土力學,2007,28(11):2485-2488.

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

[10]  夏祥,李海波,李俊如,等.巖體爆生裂紋的數值模擬[J].巖土力學,2006,27(11):1987-1991.