朱 亮，陳 明,2，盧文波,2，嚴(yán) 鵬,2

(1.武漢大學(xué)水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072； 2.武漢大學(xué)水工巖石力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072)

?

輪廓爆破下柱狀節(jié)理巖體開裂過(guò)程的數(shù)值模擬

朱 亮1，陳 明1,2，盧文波1,2，嚴(yán) 鵬1,2

(1.武漢大學(xué)水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072； 2.武漢大學(xué)水工巖石力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072)

以白鶴灘水電站柱狀節(jié)理玄武巖為研究對(duì)象，根據(jù)現(xiàn)場(chǎng)的柱狀節(jié)理分布統(tǒng)計(jì)資料生成數(shù)值模型，模擬邊坡開挖輪廓爆破下柱狀節(jié)理巖體的開裂過(guò)程。結(jié)果表明：輪廓爆破時(shí)，柱狀節(jié)理面對(duì)巖體裂紋的發(fā)展起著導(dǎo)向和控制作用；爆炸荷載作用下炮孔周圍的節(jié)理面首先張開，然后裂縫沿著炮孔連心線附近的節(jié)理面擴(kuò)展，最后形成貫通的由大量柱狀節(jié)理面連接而成的裂縫；預(yù)裂爆破條件下，預(yù)裂縫寬度先增大后減少，最后達(dá)到穩(wěn)定的開度，而且預(yù)裂爆破在柱狀節(jié)理巖體中產(chǎn)生的裂紋范圍比光面爆破產(chǎn)生的大。

爆炸力學(xué)；開裂過(guò)程；輪廓爆破；巖體；柱狀節(jié)理

中國(guó)西南已建、在建和將建的多個(gè)水電站中，均揭露出大量二疊系峨眉山玄武巖組。如四川樂(lè)山境內(nèi)的銅街子水電站和金沙江中下游的溪洛渡和白鶴灘水電站，均存在柱狀節(jié)理玄武巖，其中以白鶴灘水電站最為典型，其邊坡、壩基和地下廠房勘探平洞均揭露大量密集發(fā)育的細(xì)長(zhǎng)狀柱狀節(jié)理的玄武巖組。

柱狀節(jié)理玄武巖一般較為破碎，巖石質(zhì)量指標(biāo)(RQD)值較低，完整性較差，是西南地區(qū)一種比較典型的非貫通節(jié)理巖體[1]。針對(duì)柱狀節(jié)理巖體特殊的構(gòu)造和對(duì)工程產(chǎn)生的影響，許多人做了研究。R.Mallet[2]在1875年就對(duì)玄武巖柱狀節(jié)理特征和形成機(jī)制的進(jìn)行了研究。劉海寧等[3]開展的柱狀節(jié)理玄武巖體的脆性材料模型真三軸物理模擬試驗(yàn)表明其強(qiáng)度變形的各向異性特征顯著。石安池等[4]系統(tǒng)分析了白鶴灘柱狀節(jié)理玄武巖的基本力學(xué)特性和不同試驗(yàn)加載條件下的巖體變形機(jī)制。徐衛(wèi)亞等[5-6]認(rèn)為RMR法和GSI法能較好的反映白鶴灘柱狀節(jié)理巖體的工程力學(xué)特性，采用各向異性彈塑性模型計(jì)算更能反映柱狀節(jié)理的特殊性。朱珍德等[7]認(rèn)為Cosserat各向異性本構(gòu)模型比各向同性本構(gòu)模型更能反映工程實(shí)際，對(duì)工程優(yōu)化設(shè)計(jì)有著重要的意義?？梢?，研究者對(duì)柱狀節(jié)理玄武巖的成因、力學(xué)特性、本構(gòu)關(guān)系及靜力破壞特性等方面的研究均取得了很大進(jìn)展，但其爆破開挖動(dòng)力擾動(dòng)荷載下的破壞特性等方面有待進(jìn)一步研究。王志亮等[8]采用LS-DYNA和UDEC相結(jié)合的方法，研究了爆炸荷載作用下節(jié)理巖體中裂紋的發(fā)展過(guò)程。結(jié)果顯示UDEC中自帶的泰森多邊形節(jié)理生成器在構(gòu)造節(jié)理面分布及體現(xiàn)爆炸荷載作用下裂紋擴(kuò)展具有很大的優(yōu)勢(shì)。

節(jié)理巖體在爆炸荷載作用下的破壞情況與節(jié)理的存在有很大關(guān)系[9]。UDEC程序在研究節(jié)理巖體方面的方法簡(jiǎn)單、結(jié)果真實(shí)，因此本文中采用UDEC模擬研究輪廓爆破下柱狀節(jié)理巖體的開裂過(guò)程。

1 柱狀節(jié)理巖體裂紋發(fā)展機(jī)理分析

柱狀節(jié)理玄武巖是一種典型的節(jié)理巖體，巖石強(qiáng)度大，節(jié)理面發(fā)育且強(qiáng)度低。爆炸荷載作用下，炮孔壁附近的節(jié)理面容易發(fā)生破壞產(chǎn)生裂紋。根據(jù)成艷[10]對(duì)長(zhǎng)白山第四紀(jì)玄武巖隧道光面爆破的研究，爆炸荷載作用下，柱狀節(jié)理巖體中的裂紋的發(fā)展?jié)M足：

(1)

式中：K為應(yīng)變能密度因子，θ為裂紋發(fā)展的方向與已經(jīng)存在的節(jié)理裂隙之間的夾角，δ為起裂角。在爆破作用下產(chǎn)生的裂紋主要屬于張開型裂紋[11]，裂紋尖端的應(yīng)變能密度因子為

(2)

式中：K1為裂紋尖端的應(yīng)力強(qiáng)度因子；E為巖石的彈性模量；μ=0.2，為巖石的泊松比。式兩端同時(shí)對(duì)θ求導(dǎo)，結(jié)合式(1)，可得：

(3)

(4)

從式(3)～(4)可以看出在爆炸荷載作用下，節(jié)理巖體中產(chǎn)生的裂紋會(huì)沿著θ=0的方向發(fā)展，也就是沿著巖體中原有的節(jié)理裂隙的方向發(fā)展。因此，當(dāng)炮孔壁附近的節(jié)理面破壞產(chǎn)生裂紋時(shí)，柱狀節(jié)理裂隙對(duì)之后裂紋的發(fā)展起著控制和導(dǎo)向的作用。

2 柱狀節(jié)理巖體開裂模擬的數(shù)值模型

2.1 離散元模型

參照文獻(xiàn)[5]的分類方法，可將白鶴灘柱狀節(jié)理巖體分為3類：(1) 非規(guī)則小型，柱狀節(jié)理呈細(xì)長(zhǎng)狀密集發(fā)育，直徑為13～25 cm；(2) 非規(guī)則中型，柱狀節(jié)理發(fā)育不規(guī)則，未切割成完整的柱體，直徑為25～40 cm，其內(nèi)微裂隙較發(fā)育，但相互咬合，未完全切斷；(3)非規(guī)則大型，柱狀節(jié)理發(fā)育不規(guī)則，柱體粗大，直徑0.5～2.5 m，切割不完全，嵌合緊密。本文中以大量出現(xiàn)的第2類巖體為研究對(duì)象，采用Voronoi隨機(jī)節(jié)理生成器，在10 m×8 m的區(qū)域內(nèi)，設(shè)置節(jié)理平均邊長(zhǎng)0.25 m，令迭代次數(shù)小于100來(lái)控制柱面形狀，使其以四邊形和五邊形為主，模型如圖1和2所示。為了利于能量的耗散，模型下側(cè)邊界和左右邊界均設(shè)為無(wú)反射邊界。光面爆破布置3個(gè)炮孔，光爆孔的抵抗線為1.5 m，孔徑為90 mm，孔距1 m，爆破藥卷直徑25 mm。預(yù)裂爆破布置3個(gè)炮孔，預(yù)裂孔的抵抗線為5.5 m，孔徑為90 mm，孔距1 m，爆破藥卷直徑32 mm。

圖1 光面爆破數(shù)值模擬Fig.1 Numerical simulation model of smooth blasting

圖2 欲裂爆破數(shù)值模擬Fig.2 Numerical simulation model of pre-splitting blasting

輪廓爆破不可避免的會(huì)對(duì)保留巖體造成損傷，巖體中產(chǎn)生裂紋。通過(guò)監(jiān)測(cè)節(jié)理面兩側(cè)柱體面向臨空面(y方向)的速度和位移，揭示裂紋的發(fā)展變化過(guò)程。光面爆破監(jiān)測(cè)點(diǎn)位置如圖1所示，A位于炮孔下側(cè)0.03 m，B位于炮孔連心線下側(cè)0.07 m。預(yù)裂爆破過(guò)程中，通過(guò)監(jiān)測(cè)炮孔連心線上節(jié)理面兩側(cè)柱體的速度和位移，表明預(yù)裂縫發(fā)展過(guò)程。監(jiān)測(cè)點(diǎn)位置如圖2所示，a位于炮孔連心線上側(cè)0.01 m，b位于炮孔連心線下側(cè)0.02 m。

2.2 計(jì)算參數(shù)

2.2.1 爆破荷載

UDEC計(jì)算中爆炸波輸入通常被簡(jiǎn)化為三角脈沖，參數(shù)包括峰值、沖擊時(shí)間和脈沖寬度。輪廓爆破采用不耦合裝藥，這能夠有效的降低作用于炮孔壁上的峰值壓力，且爆生氣體充滿炮孔時(shí)的壓力通常小于臨界壓力。此時(shí)孔壁的初始沖擊波峰值壓力為[12]

(5)

圖3 爆破荷載曲線Fig.3 The history of blasting load

式中：ρ=1.1 g/cm3，為炸藥密度；D=4 000 m/s，為炸藥爆速；γ=3，為炸藥等熵指數(shù)；ν=1.4，為絕熱指數(shù)；dc和db分別為裝藥和炮孔的直徑；pk=200 MPa，為爆生氣體膨脹過(guò)程中的臨界壓力。計(jì)算知光面爆炸荷載峰值為17 MPa，預(yù)裂爆破荷載峰值為34 MPa。設(shè)荷載上升時(shí)間為2 ms，下降時(shí)間為6 ms，總歷程為8 ms。三角形荷載(σ)曲線見圖3。

2.2.2 巖體物理力學(xué)參數(shù)

柱狀節(jié)理構(gòu)造是玄武巖中最為常見的一種規(guī)則柱體形態(tài)的原生張性破裂構(gòu)造，UDEC中將柱狀節(jié)理巖體視為完整的巖石被巖體中的不連續(xù)面分割而成。巖石采用UDEC程序中的Mohr-Coulomb模型，節(jié)理面滿足庫(kù)侖滑移準(zhǔn)則。數(shù)值模型采用文獻(xiàn)[13-14]中的柱狀節(jié)理巖體巖石以及節(jié)理面的有關(guān)力學(xué)參數(shù)，如下表1和表2所示，表1中：pR為內(nèi)聚力，α為摩擦角，σs為抗剪強(qiáng)度；表2中：Kn為法向剛度，Kt為切向剛度，αi為內(nèi)摩擦角，ρJ為粘聚力，σt為抗拉強(qiáng)度。

表1 巖石力學(xué)參數(shù)

表2 節(jié)理面力學(xué)參數(shù)

3 模擬結(jié)果分析

3.1 巖體新增裂紋特性

圖4和圖5分別給出了光面爆破和預(yù)裂爆破的裂紋發(fā)展過(guò)程。光面爆破采用不耦合裝藥，作用于爆孔壁上的沖擊荷載峰值壓力低于柱狀節(jié)理巖體中巖石的動(dòng)極限抗壓強(qiáng)度，炮孔壁附近巖體沒(méi)有被“壓碎”。柱狀節(jié)理巖體的柱體緊密嵌合成為一個(gè)整體，節(jié)理面的抗拉強(qiáng)度很小。沖擊荷載作用下，巖體中產(chǎn)生徑向壓應(yīng)力和環(huán)向拉應(yīng)力。拉應(yīng)力大于節(jié)理面的抗拉強(qiáng)度，與連心線夾角小的節(jié)理面張開，如圖4(a)所示，炮孔壁周圍產(chǎn)生沿著柱狀節(jié)理面的張拉裂紋(紅色)。同時(shí)在壓縮應(yīng)力作用下與連心線夾角大的節(jié)理面發(fā)生剪切屈服，出現(xiàn)剪切裂紋(藍(lán)色)。裂紋優(yōu)先沿著已有裂紋的方向發(fā)展，同時(shí)相鄰炮孔具有導(dǎo)向作用，形成沿著連心線的貫通裂紋。裂紋繼續(xù)向自由面擴(kuò)張，如圖4(c)所示，此時(shí)巖體中已產(chǎn)生了大量裂隙，并能在模型上部臨空面上看到輕微鼓包現(xiàn)象。隨著保留巖體中產(chǎn)生通向自由面的貫通裂紋，爆生壓力作用結(jié)束。最后巖體在慣性作用下進(jìn)一步發(fā)生破壞和拋擲。

圖4 光面爆破裂紋發(fā)展過(guò)程Fig.4 Developing process of cracks in smooth blasting

圖5 預(yù)裂爆破裂紋發(fā)展過(guò)程Fig.5 Developing process of cracks in presplitting blasting

對(duì)比圖4和圖5，可以看出預(yù)裂爆破和光面爆破下柱狀節(jié)理巖體中初始裂紋的產(chǎn)生過(guò)程和形狀都比較相似。從2種工況的裂紋分布來(lái)看，裂紋范圍不均勻。爆破過(guò)程中，裂紋沿著已有的柱狀節(jié)理面發(fā)展。由于柱狀節(jié)理面分布隨機(jī)，有利于裂紋產(chǎn)生的方向，巖體的破壞程度大，節(jié)理分布對(duì)巖體產(chǎn)生裂紋的范圍有很大影響。

圖6 輪廓爆破細(xì)部圖Fig.6 Detail drawing of contour blasting

圖6給出了2種爆破的細(xì)部圖?？梢钥闯?，2種工況保留巖體的損傷程度不同。光面爆破下保留巖體中裂紋的最大范圍為0.31 m，預(yù)裂爆破下保留巖體的裂紋范圍為0.61 m，說(shuō)明預(yù)裂爆破對(duì)保留巖體產(chǎn)生的損傷比光面爆破大。由于預(yù)裂爆破的抵抗線以及爆炸荷載峰值壓力比光面爆破的大，因此預(yù)裂對(duì)保留巖體的損傷程度大于光面對(duì)保留巖體的損傷程度。

從圖6可以看出，光面和預(yù)裂縫的形成與柱狀節(jié)理面關(guān)系密切。相鄰炮孔同時(shí)起爆后，徑向裂隙不再沿著炮孔連心線的方向貫通，而是沿著連心線附近的節(jié)理面產(chǎn)生并貫通。爆炸荷載作用下，炮孔附近的節(jié)理面張開，產(chǎn)生裂紋。當(dāng)新裂紋擴(kuò)展到柱狀節(jié)理面時(shí)，其擴(kuò)展路徑嚴(yán)重受到其影響，裂紋擴(kuò)展歸并到柱狀節(jié)理面上。同時(shí)其他方向的裂紋擴(kuò)展受到抑制，巖體裂紋沿著節(jié)理面形成貫通，最終形成由柱狀節(jié)理面連接而成的光面和預(yù)裂縫?？梢?，當(dāng)巖體中的節(jié)理面破壞后，巖體的破壞將一直沿著對(duì)巖體的破壞起控制性作用的柱狀節(jié)理面發(fā)生，柱狀節(jié)理的存在對(duì)巖體的破壞具有導(dǎo)向作用。

3.2 裂紋發(fā)展過(guò)程

由監(jiān)測(cè)結(jié)果可以得出節(jié)理面兩側(cè)柱體的相對(duì)速度和位移，由此說(shuō)明爆炸荷載作用下巖體中裂紋的發(fā)展?fàn)顩r。圖7～8給出了2種爆破的位移(r)和速度(u)曲線。

圖7 光面爆破裂紋發(fā)展曲線Fig.7 Developing curves of the crack in smooth blasting

圖8 預(yù)裂爆破預(yù)裂縫發(fā)展曲線Fig.8 Developing curves of precracking in presplitting blasting

圖7(a)表明伴隨爆炸荷載的作用，玄武巖柱體之間逐漸產(chǎn)生相對(duì)位移，節(jié)理面之間產(chǎn)生滑移及張開變形。爆炸荷載的加載階段巖體存在負(fù)向的壓縮位移，與文獻(xiàn)[15]的理論研究及數(shù)值模擬一致。加載過(guò)程結(jié)束后，裂紋發(fā)展逐漸停止。爆炸荷載作用下，巖體中產(chǎn)生強(qiáng)烈的爆炸壓縮波，保留巖體被壓縮。應(yīng)力波過(guò)后，壓縮應(yīng)變回彈。節(jié)理面上側(cè)柱體A獲得沖擊速度并發(fā)生剛體位移，節(jié)理面張開。隨著柱體與周圍柱體之間的相互作用損耗大量能量，相對(duì)速度逐漸減少為零，裂紋開度達(dá)到穩(wěn)定(圖7(b))。

圖8監(jiān)測(cè)結(jié)果表明，預(yù)裂爆破過(guò)程中，柱狀節(jié)理面開度達(dá)到最大后緩慢減小，最后形成穩(wěn)定的預(yù)裂縫。預(yù)裂孔位于巖體深處，待開挖巖體和保留巖體在爆炸壓縮波作用下，向相反方向移動(dòng)。連心線上節(jié)理面張開，同時(shí)巖體內(nèi)部積蓄壓縮勢(shì)能。應(yīng)力波過(guò)后，壓縮應(yīng)變回彈，預(yù)裂縫逐漸減小，最后達(dá)到穩(wěn)定。

4 結(jié) 論

采用離散元UDEC數(shù)值模擬技術(shù)，對(duì)輪廓爆下柱狀節(jié)理巖體的開裂過(guò)程進(jìn)行了分析，得到以下初步結(jié)論。

(1) 輪廓爆破時(shí)，由于柱狀節(jié)理巖體的節(jié)理強(qiáng)度較低，沿著節(jié)理面產(chǎn)生張拉或剪切破壞，而柱體強(qiáng)度高，一般不產(chǎn)生貫穿柱體的裂紋。

(2) 光面爆破和預(yù)裂爆破的成縫過(guò)程基本一致；炮孔壁附近的柱狀節(jié)理面首先張開，裂紋沿著節(jié)理面貫通，形成由柱狀節(jié)理面連接而成的輪廓面；柱狀節(jié)理面對(duì)裂縫的發(fā)展起著控制和導(dǎo)向作用。

(3) 光面爆破時(shí)，保留巖體在爆炸應(yīng)力波的作用下被壓縮，作用過(guò)程結(jié)束后，巖體應(yīng)變回彈，節(jié)理面張開產(chǎn)生裂紋；預(yù)裂爆破中，預(yù)裂縫的形成是一個(gè)先張開后逐漸縮小穩(wěn)定的過(guò)程；光面爆破下保留巖體的裂紋范圍比預(yù)裂爆破的小。

[1] 張倬元,王士天,王蘭生.工程地質(zhì)分析原理[M].北京:地質(zhì)出版社,1994.

[2] Mallet R. On the origin and mechanism of production of the prismatic (or columnar) structure of basalt[J]. Philosophical Magazine, 1875,4(50):122-135,201-226.

[3] 劉海寧,王俊梅,王思敬.白鶴灘柱狀節(jié)理巖體真三軸模型試驗(yàn)研究[J].巖土力學(xué),2010,31(1):163-171. Liu Hai-ning, Wang Jun-mei, Wang Si-jing. Experimental research of columnar jointed basalt with true apparatus at Baihetan Hydropower Station[J]. Rock and Soil Mechanics, 2010,31(1):163-171.

[4] 石安池,唐鳴發(fā),周其健.金沙江白鶴灘水電站柱狀節(jié)理玄武巖巖體變形特性研究[J].巖石力學(xué)與工程學(xué)報(bào),2008,27(10):2079-2086. Shi An-chi, Tang Ming-fa, Zhou Qi-jian. Research of deformation characteristics of columnar jointed basalt at Baihetan Hydropower Station[J]. Chinese Journal of Rock Mechanics and Engineering, 2008,27(10):2079-2086.

[5] 徐衛(wèi)亞,鄭文棠,石安池.水利工程中的柱狀節(jié)理巖體分類及質(zhì)量評(píng)價(jià)[J].水利學(xué)報(bào),2011,42(3):262-279. Xu Wei-ya, Zheng Wen-tang, Shi An-chi. Classification and quality assessment of irregular columnar jointed basaltic rock mass for hydraulic engineering[J]. Journal of Hydraulic Engineering, 2011,42(3):262-279.

[6] 徐衛(wèi)亞,鄭文棠,寧宇,等.柱狀節(jié)理壩基巖體三維各向異性數(shù)值分析[J].巖土力學(xué),2010,31(3):949-955. Xu Wei-ya, Zheng Wen-tang, Ning Yu, et al. 3D anisotropic numerical analysis of rock mass with columnar for dam foundation[J]. Rock and Soil Mechanics, 2010,31(3): 949-955.

[7] 朱珍德,秦天昊,王士宏.基于Cosserat理論的柱狀節(jié)理巖體各向異性本構(gòu)模型研究[J].巖石力學(xué)與工程學(xué)報(bào),2010,29(2):4068-4076. Zhu Zhen-de, Qin Tian-hao, Wang Shi-hong, et al. Study of anisotropic constitutive model of columnar jointed rock masses based on cosserat theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2010,29(2):4068-4076.

[8] Wang Z L, Konietky H. Modeling of blast-induced fractures in jointed rock masses[J]. Engineering Fracture Mechanics, 2009,76:1945-1955.

[9] 劉紅巖,秦四清,楊軍.爆炸荷載下巖石破壞的數(shù)值流形方法模擬[J].爆炸與沖擊,2007,27(1):50-56. Liu Hong-yan,Qin Si-qing, Yang Jun. Simulation of rock failure by numerical manifold method under blasting load[J]. Explosion and Sock Waves, 2007,27(1):50-56.

[10] 成艷.長(zhǎng)白山第四紀(jì)玄武巖隧洞光面爆破研究[D].吉林:吉林大學(xué)建設(shè)工程學(xué)院,2012.

[11] 齊景嶽,劉正雄.隧道現(xiàn)代爆破技術(shù)[M].北京:中國(guó)鐵道出版社,1995.

[12] Ding D, Zhu C. Estimaing amount of explosive for fracture plane control blasting with notched boreholes[J]. Transactions of NFsoc, 1993,3(2):18-21.

[13] 邴鳳山.中國(guó)水力發(fā)電年鑒2010[M].北京:中國(guó)電力出版社,2011.

[14] 朱道建,楊林德,蔡永昌.柱狀節(jié)理巖體壓縮破壞過(guò)程模擬及機(jī)制分析[J].巖石力學(xué)與工程學(xué)報(bào),2009,28(4):716-724. Zhu Dao-jian, Yang Lin-de, Cai Yong-chang. Simulation of compressive failure of columnar jointed rock mass and its failure mechanism analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2009,28(4):716-724.

[15] 金李,盧文波,陳明,等.節(jié)理巖體的爆破松動(dòng)機(jī)理[J].爆炸與沖擊,2009,29(5):274-280. Jin Li, Lu Wen-bo, Chen Ming, et al. Mechanism of jointed rock loosing under blasting load[J]. Explosion and Shock Waves, 2009,29(5):274-280.

(責(zé)任編輯 王小飛)

Numerical simulation on the cracking process of columnar jointed rock masses in contour blasting

Zhu Liang1, Chen Ming1,2, Lu Wen-bo1,2, Yan Peng1,2

(1.StateKeyLaboratoryofWaterResourcesandHydropowerEngineeringScience,WuhanUniversity,Wuhan430072,Hubei,China; 2.KeyLaboratoryofRockMechanicsinHydraulicStructuralEngineering,MinistryofEducation,WuhanUniversity,Wuhan430072,Hubei,China)

With Baihetan columnar joints basalt as the research object, numerical model is generated according to the local distribution statistics. Based on numerical analysis, the cracking process of columnar jointed rock masses with the slope excavated by contour blasting is simulated. The results show that the columnar joints control and lead the generation of cracking in contour blasting. The joints around the hole are firstly opened by the explosion load. Then they extend along the joints near the defiled line of holes. Finally the cracks formed by a large number of columnar joints are connected. The process of the width of pre-splitting crack is that it reduces to stabile opening after increases under the condition of pre-splitting blasting. And the range of cracks in pre-splitting blasting is larger than these in smooth blasting.

mechanics of explosion; cracking process; contour blasting; rock; columnar joints

10.11883/1001-1455(2015)04-0555-06

2013-01-10；

2014-05-06

國(guó)家自然科學(xué)基金項(xiàng)目(51279146)；國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2011CB013501)

朱 亮(1990- )，男，碩士研究生； 通訊作者： 陳 明，whuchm@yahoo.com.cn。

O383 國(guó)標(biāo)學(xué)科代碼： 13035

A