王云飛,鄭曉娟,焦華喆,陳峰賓,趙洪波

(1.河南理工大學(xué)土木工程學(xué)院，河南 焦作 454000；2.焦作師范高等?？茖W(xué)校管理學(xué)院，河南 焦作 454000)

花崗巖破壞過程能量演化機制與能量屈服準(zhǔn)則*

王云飛1,鄭曉娟2,焦華喆1,陳峰賓1,趙洪波1

(1.河南理工大學(xué)土木工程學(xué)院，河南 焦作 454000；2.焦作師范高等?？茖W(xué)校管理學(xué)院，河南 焦作 454000)

為了明確巖石破壞的能量演化特性，結(jié)合單軸實驗和顆粒流程序獲得花崗巖的細(xì)觀力學(xué)參數(shù)，進行不同應(yīng)力狀態(tài)的花崗巖實驗，研究不同圍壓下花崗巖破壞過程的能量演化機理并推導(dǎo)能量屈服準(zhǔn)則。獲得以下主要結(jié)論：花崗巖破壞過程中低圍壓下內(nèi)部損傷出現(xiàn)較早而高圍壓較晚，表明低圍壓花崗巖內(nèi)部損傷是漸進發(fā)展過程，而高圍壓下內(nèi)部損傷一旦出現(xiàn)便快速發(fā)展破壞；高圍壓花崗巖峰值前一定應(yīng)變范圍彈性應(yīng)變能基本保持不變，吸收的能量全部轉(zhuǎn)化為耗散能，表明高圍壓破壞時花崗巖內(nèi)部損傷程度嚴(yán)重；彈性應(yīng)變能經(jīng)歷不斷積累并達到彈性儲能極限而后減小的變化過程，而彈性儲能極限與圍壓之間存在線性變化規(guī)律，因此高圍壓下巖體開挖卸荷時極易誘發(fā)大量彈性應(yīng)變能的急劇釋放，引起圍巖失穩(wěn)甚至發(fā)生巖爆；花崗巖峰值破壞時的能量比與圍壓無關(guān)，為一定值；基于能量原理導(dǎo)出了能量屈服準(zhǔn)則，該準(zhǔn)則包含巖性參數(shù)和所有主應(yīng)力，能夠綜合反映巖石破壞影響因素。

固體力學(xué)；能量屈服準(zhǔn)則；顆粒流程序；花崗巖；能量演化機理

外力對巖石做功使它發(fā)生變形甚至破壞，是個能量轉(zhuǎn)化過程，巖石受載破壞是能量驅(qū)動的結(jié)果，因而研究巖石破壞過程的能量轉(zhuǎn)化機制具有重要意義。在這個方面，已開展了大量的研究：巖石整體破壞與能量耗散和釋放的內(nèi)在聯(lián)系[1-2]，內(nèi)部損傷與宏觀破壞的能量解釋[3]，卸荷引起巖爆的能量分析[4-5]，加卸載速率[6]、圍壓[7]對巖石破壞中能量變化的影響，采用顆粒流軟件對巖土類材料的研究[8-9]。在屈服準(zhǔn)則方面，有深部巖體強度準(zhǔn)則[10]，Mohr-Coulomb 強度理論及其發(fā)展[11]，從能量轉(zhuǎn)化角度建立的能量屈服準(zhǔn)則[12-13]。

由上可見，對巖石破壞的能量轉(zhuǎn)化研究較多，但對花崗巖全應(yīng)力應(yīng)變曲線的不同變形階段能量轉(zhuǎn)化機理研究很少，而花崗巖又是大型地下工程、礦山巷道施工面臨的普遍巖體。因此，在本文中，對花崗巖在不同圍壓下各變形階段的能量轉(zhuǎn)化進行詳細(xì)研究，并從彈性儲能極限的新角度建立能量屈服準(zhǔn)則，擬為分析地下巖體工程的穩(wěn)定性提供參考。

1 顆粒流理論與模型

P.A.Cundall等[14]引入分子動力學(xué)思想創(chuàng)建的顆粒流理論，著重從細(xì)觀力學(xué)角度分析材料的破壞特性，能夠研究線彈性直至破壞的大變形過程。它主要用于巖石材料的研究，無需定義材料的本構(gòu)模型，只需設(shè)定細(xì)觀力學(xué)參數(shù)和顆粒間的黏結(jié)模型，便可表達復(fù)雜的非線性應(yīng)力應(yīng)變關(guān)系。

顆粒流程序(particle flow code, PFC)采用平行黏結(jié)模型表征巖石類材料，平行黏結(jié)發(fā)生在半徑為R的圓形范圍內(nèi)，如圖1所示。A和B表示兩個接觸顆粒，OA、OB和OC分別表示A、B顆粒和平行黏結(jié)的中心，R和L分別為平行黏結(jié)半徑和厚度，F(xiàn)和M為平行黏結(jié)所承受的力和彎矩，F(xiàn)分解為法向和切向分量Fn、Fs，M分解為法向和切向分量Mn、Ms，A為平行黏結(jié)截面面積，I和J分別為平行黏結(jié)截面慣性矩和極慣性矩。

圖1 平行黏結(jié)模型Fig.1 Parallel bond model

作用在平行黏結(jié)上的法向應(yīng)力σ和切向應(yīng)力τ分別為[15]：

(1)

(2)

試樣受載后微小顆粒間受力增加，當(dāng)達到其黏結(jié)強度便發(fā)生斷裂形成內(nèi)部損傷，損傷積累便形成宏觀裂紋導(dǎo)致巖石破壞，這與實際巖石的破壞機制相符。

2 花崗巖特性與細(xì)觀力學(xué)參數(shù)

花崗巖密度為2.656 g/cm3，將它加工為?50 mm×100 mm的試樣進行單軸實驗，測得花崗巖的單軸抗壓強度為149.67 MPa，彈性模量為67.43 GPa?；◢弾rPFC試樣由服從正態(tài)分布的大小顆粒及其之間的黏結(jié)組成，最大和最小粒徑之比為1.66，這與巖石由微小礦物顆粒黏結(jié)組成的結(jié)構(gòu)相似。PFC通過設(shè)置細(xì)觀力學(xué)參數(shù)來表征巖石的宏觀力學(xué)行為，由室內(nèi)單軸實驗曲線反演獲得花崗巖的細(xì)觀力學(xué)參數(shù)，顆粒間的摩擦因素為0.5，顆粒彈性模量Ec=67.96 GPa，顆粒法向與切向剛度比kn/ks=1.8，平行黏結(jié)半徑系數(shù)λ=1.0，平行黏結(jié)彈性模量Ec=67.96 GPa，平行黏結(jié)法向與切向剛度比kn/ks=1.8，平行黏結(jié)法向強度的平均值σn,mean=98 MPa，標(biāo)準(zhǔn)差σn,dev=24.5 MPa，平行黏結(jié)切向強度的平均值τn,mean=294 MPa，標(biāo)準(zhǔn)差τn,dev=73.5 MPa。

3 能量演化機制

3.1 能量原理

工程修建中引起的加卸載活動，改變了巖體的受力條件。在巖體所受應(yīng)力改變過程中外力作功，其中一部分以彈性應(yīng)變能的形式存儲于巖體中，另一部分為巖體內(nèi)部損傷所消耗。

圖2 耗散能與可釋放應(yīng)變能的關(guān)系Fig.2 Relationship between dissipation strain energy and releasable strain energy

取單位巖體進行破壞過程中的能量轉(zhuǎn)化分析，假設(shè)這個過程中無熱交換和聲發(fā)射等能量損失，則外力對單位巖體的總輸入能量U可表示為[1]：

U=Ud+Ue

(3)

式中：Ud為單位巖石耗散能，Ue為巖石彈性變形儲存的能量，即彈性應(yīng)變能。

圖2為耗散能Ud與彈性應(yīng)變能Ue的關(guān)系示意圖，其中陰影部分面積為由于變形而存儲于巖石中的彈性應(yīng)變能Ue，卸載模量Eu與應(yīng)力應(yīng)變曲線及水平軸圍成的面積為巖石內(nèi)部損傷形成所需的耗散能Ud。

在三向應(yīng)力狀態(tài)下，巖石單元的能量計算公式[16]可表示為：

(4)

(5)

(6)

(7)

式中：σi(i=1,2,3)為主應(yīng)力，εi(i=1,2,3)為主應(yīng)變，規(guī)定壓應(yīng)變?yōu)檎瑒t巖石受壓側(cè)向膨脹ε2、ε3為負(fù)值。εei(i=1,2,3)為相應(yīng)的彈性應(yīng)變，E0為初始彈模，Eu和νu分別為卸載彈模和泊松比。尤明慶等[17]、余賢斌等[18]指出，可用E0近似代替Eu，為了簡化計算本文中采用了：Eu=E0，νu=ν。

3.2 能量轉(zhuǎn)化分析

外力對巖石作功使它變形，在巖石內(nèi)部存儲彈性應(yīng)變能，同時外力作用下巖石內(nèi)部損傷發(fā)展也消耗能量。當(dāng)巖石積蓄的彈性能達到極限時，便釋放轉(zhuǎn)變?yōu)楹纳⒛芷茐膸r石，巖石破壞是能量驅(qū)動的結(jié)果。

圖3為花崗巖破壞的彈性能和耗散能轉(zhuǎn)化機制?？梢姡?/p>

(1)花崗巖在壓縮的彈性階段，吸收的總能量基本都以可釋放彈性應(yīng)變能的形式存儲，原因在于，彈性變形階段花崗巖內(nèi)部損傷極其微少，耗散能基本為零。

(2)隨著外力作用進入屈服階段，雖然花崗巖吸收的總能量仍以可釋放彈性應(yīng)變能的形式存儲，但耗散能也在不斷增加，內(nèi)部損傷不斷增多，特別在峰值點附近耗散能增速顯著，峰值附近內(nèi)部損傷急劇增多。在這個階段，彈性應(yīng)變能的增速不斷減小，而耗散能的增速不斷增大。

(3)峰后階段，外力作功巖石繼續(xù)吸收能量，但彈性應(yīng)變能在峰后快速釋放，相應(yīng)的耗散能急劇增加。釋放的彈性應(yīng)變能轉(zhuǎn)化為巖石內(nèi)部損傷和裂隙擴展所需的表面能。在這個階段，外力作功加彈性能的釋放，促使巖石內(nèi)部裂隙的快速擴展貫通，導(dǎo)致巖石失穩(wěn)破壞。

圖3 花崗巖破壞過程能量轉(zhuǎn)化機制Fig.3 Energy conversion mechanism of granite in failure process

由巖石吸收總能量、彈性應(yīng)變能和耗散能的轉(zhuǎn)化可見，巖石吸收的總能量在彈性階段增速不斷增大，在屈服階段增速逐漸減小，在峰后增速基本恒定。彈性應(yīng)變能經(jīng)歷先增加后減小的變化趨勢，在應(yīng)力峰值處達到最大儲能極限。耗散能從屈服點附近開始出現(xiàn)，增速逐漸增大，在峰值處增速達到最大，在殘余強度階段增速基本與總吸收能增速相同。

對比不同圍壓下花崗巖能量轉(zhuǎn)化可見，隨著圍壓的增大，花崗巖耗散能出現(xiàn)對應(yīng)的軸向應(yīng)變不斷增加，即外力作用下在單軸和低圍壓時花崗巖內(nèi)部損傷出現(xiàn)較早(見圖3(a))，高圍壓時內(nèi)部損傷出現(xiàn)較晚(見圖3(d))，這也表明單軸和低圍壓下花崗巖內(nèi)部損傷是漸進發(fā)展過程，而高圍壓下內(nèi)部損傷一旦出現(xiàn)便快速發(fā)展。隨著圍壓的增加，屈服階段彈性應(yīng)變能的增速逐漸減小，高圍壓下峰值前一定應(yīng)變范圍內(nèi)彈性應(yīng)變能基本不再增加，表明高圍壓下峰值附近區(qū)間花崗巖吸收的總能量都轉(zhuǎn)化為耗散能用于內(nèi)部損傷和裂隙發(fā)展所需的表面能，進一步說明花崗巖在高圍壓下峰值破壞時內(nèi)部損傷程度較高。彈性應(yīng)變能達到極值不能繼續(xù)增大，便進入峰后釋放階段，單軸和低圍壓下彈性應(yīng)變能急劇釋放對應(yīng)巖石的突然失穩(wěn)破壞(見圖3(a))；高圍壓下由于圍壓的約束作用能量釋放困難，更沒有突然釋放的條件，因而彈性應(yīng)變能減小到某個值后，在殘余強度階段仍保持較大值，對應(yīng)于花崗巖的延性破壞(見圖3(d))。

圖3還顯示，在花崗巖變形破壞過程中，彈性應(yīng)變能經(jīng)歷從最初的不斷積蓄增加到峰后釋放減小期間存在儲能最大值，稱為彈性儲能極限。圖4為不同圍壓下的花崗巖彈性儲能極限。可以看出，彈性儲能極限與圍壓存在良好的線性關(guān)系，單軸時極限彈性應(yīng)變能為0.154 MJ/m3，50 MPa圍壓時為0.647 MJ/m3，100 MPa圍壓時為1.281 MJ/m3。

由圖4可見，高圍壓下花崗巖破壞時積蓄的彈性應(yīng)變能較大，又結(jié)合圖3可知，高圍壓下花崗巖破壞時的內(nèi)部損傷嚴(yán)重，因此，處于較高三向壓力狀態(tài)的花崗巖體在地下工程修建開挖卸荷時，圍壓突然降低(巷道周邊接近單軸應(yīng)力狀態(tài))極易誘發(fā)積蓄在巖體中大量彈性能的急劇釋放，導(dǎo)致圍巖失穩(wěn)破壞甚至發(fā)生巖爆。

由圖5可見，花崗巖在峰值破壞時，彈性應(yīng)變能與吸收總能量的比Ue/U和耗散能與吸收總能量的比Ud/U，與圍壓基本沒有關(guān)系。在花崗巖破壞時，彈性應(yīng)變能與吸收總能量的比為0.75，相應(yīng)的耗散能與吸收總能量的比為0.25。

圖4 花崗巖彈性儲能極限與圍壓的關(guān)系Fig.4 Relationship between granite’s maximum elastic strain energy and confining pressures

圖5 花崗巖峰值能量比與圍壓關(guān)系Fig.5 Relationship between granite’s peak energy ratio and confining pressures

4 能量屈服準(zhǔn)則

以上表明，花崗巖單位巖體彈性儲能極限Ue,max與圍壓存在線性變化規(guī)律，即：

Ue,max=U0+kσ3

(8)

式中：U0為花崗巖單軸破壞的儲能極限，k為與巖性相關(guān)的常數(shù)。

由圖3可見，不同圍壓下花崗巖單位巖體彈性應(yīng)變能隨應(yīng)變增加不斷增大，當(dāng)達到儲能極限時，可釋放彈性應(yīng)變能導(dǎo)致巖體破壞。因而，巖體單元破壞時有：

(9)

將式(8)代入式(9)，并用初始彈性模量E0替代卸荷彈性模量Eu，得：

(10)

由式(7),花崗巖的單軸壓縮儲能極限U0為：

(11)

將式(11)代入式(10)，并經(jīng)化簡可得：

(12)

令2E0k=K，則有能量屈服準(zhǔn)則：

(13)

圖6 花崗巖能量屈服準(zhǔn)則Fig.6 Energy yield criterion of granite

式中：K為與巖性相關(guān)的參數(shù)。通過實驗獲得巖性參數(shù)K，便可利用式(13)，判定不同應(yīng)力狀態(tài)下巖體工程的穩(wěn)定性。

為了驗證能量屈服準(zhǔn)則的計算精度，采用不同圍壓強度的線性擬合，給出能量屈服準(zhǔn)則的計算強度：σc=161.06 MPa，K=1.618 5 GPa。

圖6為能量屈服準(zhǔn)則計算結(jié)果，計算值與實驗結(jié)果很接近，最大相對誤差為9.89%。由式(13)可見，基于能量原理建立的能量屈服準(zhǔn)則物理意義明確，即花崗巖達到最大儲能極限巖石發(fā)生破壞，且能量準(zhǔn)則包括彈性模量、泊松比和所有主應(yīng)力，能夠從本質(zhì)上反映巖石破壞的綜合影響因素。

5 結(jié) 論

(1)花崗巖低圍壓時內(nèi)部損傷出現(xiàn)較早，高圍壓下內(nèi)部損傷出現(xiàn)較晚，表明花崗巖在低圍壓下?lián)p傷是漸進發(fā)展過程，而高圍壓下內(nèi)部損傷一旦出現(xiàn)便快速發(fā)展。

(2)花崗巖在高圍壓下峰值之前一定應(yīng)變范圍內(nèi)彈性應(yīng)變能基本不再增加，吸收的能量全部轉(zhuǎn)化為耗散能，說明花崗巖高圍壓下峰值破壞時內(nèi)部損傷程度較高。

(3)破壞過程中，彈性應(yīng)變能經(jīng)歷不斷積累增大到儲能極限而后減小的過程，其彈性儲能極限與圍壓之間存在良好的線性關(guān)系。

(4)高圍壓下花崗巖積蓄的彈性應(yīng)變能較大，因而工程開采卸圍壓時，極易誘發(fā)大量彈性應(yīng)變能的急劇釋放，導(dǎo)致圍巖失穩(wěn)破壞甚至發(fā)生巖爆。

(5)花崗巖在峰值破壞時，彈性應(yīng)變能與吸收總能量的比為0.75，耗散能與吸收總能量的比為0.25，與圍壓沒有關(guān)系。

(6)基于巖石破壞的能量演化機制導(dǎo)出能量屈服準(zhǔn)則，該準(zhǔn)則物理意義明確，包含巖性參數(shù)和所有主應(yīng)力，能夠從本質(zhì)上反映巖石破壞的綜合影響因素。

[1] 謝和平，鞠楊，黎立云.基于能量耗散與釋放原理的巖石強度與整體破壞準(zhǔn)則[J].巖石力學(xué)與工程學(xué)報，2005，24(17)：3003-3010. Xie Heping, Ju Yang, Li Liyun. Criteria for strength and structural failure of rocks based on energy dissipation and energy release principles[J]. Chinese Journal of Rock Mechanics and Engineering, 2005,24(17):3003-3010.

[2] 謝和平,彭瑞東,鞠楊.巖石變形破壞過程中的能量耗散分析[J].巖石力學(xué)與工程學(xué)報,2004,23(21):3565-3570. Xie Heping, Peng Ruidong, Ju Yang. Energy dissipation of rock deformation and fracture[J]. Chinese Journal of Rock Mechanics and Engineering, 2004,23(21):3565-3570.

[3] 黃達，黃潤秋，張永興.粗晶大理巖單軸壓縮力學(xué)特性的靜態(tài)加載速率效應(yīng)及能量機制試驗研究[J].巖石力學(xué)與工程學(xué)報，2012，31(2)：245-255. Huang Da, Huang Runqiu, Zhang Yongxing. Experimental investigations on static loading rate effects on mechanical properties and energy mechanism of coarse crystal grain marble under uniaxial compression[J]. Chinese Journal of Rock Mechanics and Engineering, 2012,31(2):245-255.

[4] 陳衛(wèi)忠，呂森鵬，郭小紅，等.基于能量原理的卸圍壓試驗與巖爆判據(jù)研究[J].巖石力學(xué)與工程學(xué)報，2009，28(8):1530-1540. Chen Weizhong, Lu Senpeng, Guo Xiaohong, et al. Research on unloading confining pressure tests and rockburst criterion based on energy theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2009,28(8):1530-1540.

[5] Sanchidrian J A, Segarra P, López L M. Energy components in rock blasting[J]. Rock Mechanics and Mining Sciences, 2007,44(1):130-147.

[6] 黎立云，謝和平，鞠楊，等.巖石可釋放應(yīng)變能及耗散能的實驗研究[J].工程力學(xué)，2011，28(3):35-40. Li Liyun, Xie Heping, Ju Yang, et al. Experimental investigations of releasable energy and dissipative energy within rock[J]. Engineering Mechanics, 2011,28(3):35-40.

[7] 尤明慶，華安增.巖石試樣破壞過程的能量分析[J].巖石力學(xué)與工程學(xué)報，2002，21(6)：778-781. You Mingqing, Hua Anzeng. Energy analysis on failure process of rock specimens[J]. Chinese Journal of Rock Mechanics and Engineering, 2002,21(6):778-781.

[8] Zhou Yu, Wu Shunchuan, Gao Yongtao, et al. Macro and meso analysis of jointed rock mass triaxial compression test by using equivalent rock mass (ERM) technique[J]. Journal of Central South University, 2014,21(3):1125-1135.

[9] 羅勇,龔曉南,連峰.三維離散顆粒單元模擬無黏性土的工程力學(xué)性質(zhì)[J].巖土工程學(xué)報,2008,30(2):292-297. Luo Yong, Gong Xiaonan, Lian Feng. Simulation of mechanical behaviors of granular materials by three-dimensional discrete element method based on particle flow code[J]. Chinese Journal of Geotechnical Engineering, 2008,30(2):292-297.

[10] 周小平,錢七虎,楊海清.深部巖體強度準(zhǔn)則[J].巖石力學(xué)與工程學(xué)報,2008,27(1):117-123. Zhou Xiaoping, Qian Qihu, Yang Haiqing. Strength criteria of deep rock mass[J]. Chinese Journal of Rock Mechanics and Engineering, 2008,27(1):117-123.

[11] 俞茂宏,昝月穩(wěn),范文,等.20 世紀(jì)巖石強度理論的發(fā)展:紀(jì)念 Mohr-Coulomb 強度理論100周年[J].巖石力學(xué)與工程學(xué)報,2000,19(5):545-550. Yu Maohong, Zan Yuewen, Fan Wen, et al. Advances in strength theory of rock in 20 century: 100 years inmemory of the Mohr-Coulomb strength theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2000,19(5):545-550.

[12] 高紅,鄭穎人,馮夏庭.巖土材料能量屈服準(zhǔn)則研究[J].巖石力學(xué)與工程學(xué)報,2007,26(12):2437-2443. Gao Hong, Zheng Yingren, Feng Xiating. Study on energy yield criterion of geomaterials[J]. Chinese Journal of Rock Mechanics and Engineering, 2007,26(12):2437-2443.

[13] 周輝,李震,楊艷霜,等.巖石統(tǒng)一能量屈服準(zhǔn)則[J].巖石力學(xué)與工程學(xué)報,2013,32(11):2170-2185. Zhou Hui, Li Zhen, Yang Yanshuang, et al. Unified energy yield criterion of rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2013,32(11):2170-2185.

[14] Cundall P A, Strack O D. A discrete numerical model for granula assemblies[J]. Géotechnique, 1979,29(1):47-65.

[15] Itasca Consulting Group. PFC3D: Particle flow code in 3 dimensions[R]. Minneapolis, USA: Itasca Consulting Group, 2008.

[16] Solecki R, Conant R J. Advanced mechanics of materials[M]. London: Oxford University Press, 2003.

[17] 尤明慶,蘇承東.大理巖試樣循環(huán)加載強化作用的試驗研究[J].固體力學(xué)學(xué)報,2009,29(1):66-72. You Mingqing, Su Chengdong. Experimental study on strengthening of marble specimen in cyclic loading of uniaxial or pseudo-triaxial compression[J]. Chinese Journal of Solid Mechanics, 2009,29(1):66-72.

[18] 余賢斌,謝強,李心一,等.巖石直接拉伸與壓縮變形的循環(huán)加載實驗與雙模量本構(gòu)模型[J].巖土工程學(xué)報,2005,27(9):988-993. Yu Xianbin, Xie Qiang, Li Xinyi, et al. Cycle loading tests of rock samples under direct tension and compression and bi-modular constitutive model[J]. Chinese Journal of Geotechnical Engineering, 2005,27(9):988-993.

(責(zé)任編輯 丁 峰)

Energy evolution mechanism and energy yield criterion in granite’s failure process

Wang Yunfei1, Zheng Xiaojuan2, Jiao Huazhe1,Cheng Fengbin1, Zhao Hongbo1

(1.SchoolofCivilEngineering,HenanPolytechnicUniversity,Jiaozuo454000,Henan,China;2.ManagementSchool,JiaozuoTeachersCollege,Jiaozuo454000,Henan,China)

To understand the energy evolution mechanism in the rock failure process, this paper firstly obtained the meso-mechanical parameters of granite using uniaxial compression experiments and particle flow codes, then tested the granite under different confining pressures and finally analyzed its energy evolution mechanism in the failure process and deduced its energy yield criterion. The main results are as follows: The internal damage of granite in the failure process occurs earlier under lower confining pressures while later under higher confining pressures, which shows that the internal damage under lower confining pressures is a progressive development process but under higher confining pressures the internal damage rapidly develops into failure once it occurs. The granite’s elastic strain energy remains constant in a certain strain range before the peak under higher confining pressures,and the overall energy absorbed transforms into dissipation energy, which shows that the granite internal damage under higher confining pressures is more severe. The elastic strain energy increases and reaches the elastic strain energy limit and then decreases. There exists a linear relationship between the elastic strain energy limit and the confining pressure, therefore rock excavation under high confining pressures is likely to induce a rapid release of a large amount of elastic strain energy which causes the surrounding rock to become unstable and even to burst. The energy ratio at the granite’s peak failure is a definite value and independent of the confining pressure. The energy yield criterion is derived based on the principle of energy. It includes lithology parameters and all principal stresses and can reflect the comprehensive factors influencing the rock failure.

solid mechanics; energy yield criterion; particle flow code; granite; energy evolution mechanism

10.11883/1001-1455(2016)06-0876-07

2015-03-04； < class="emphasis_bold">修回日期：2015-05-20

2015-05-20

國家自然科學(xué)基金項目 (51104057)；煤炭工業(yè)協(xié)會科學(xué)技術(shù)研究計劃項目(MTKJ2013-338)； 河南省教育廳重點項目(13A440323)；河南省高?？萍紕?chuàng)新團隊支持計劃項目(15IRTSTHN029)

王云飛(1978— )，男，博士,wyf_ustb@126.com。

O346.4;TU452 <國標(biāo)學(xué)科代碼：1303540 class="emphasis_bold"> 國標(biāo)學(xué)科代碼：1303540 文獻標(biāo)志碼：A國標(biāo)學(xué)科代碼：1303540

A