劉鵬飛，范俊奇，郭佳奇，朱斌忠

（1. 河南理工大學(xué)土木工程學(xué)院，河南 焦作 454003；2. 煤炭安全生產(chǎn)與清潔高效利用省部共建協(xié)同創(chuàng)新中心，河南 焦作 454003；3. 軍事科學(xué)院國防工程研究院，河南 洛陽 471023）

花崗巖是大陸地殼的主要組成部分，是一種巖漿在地表以下凝結(jié)形成的巖漿巖，屬于深層侵入巖。同時(shí)花崗巖也是一種典型的硬脆巖，具有強(qiáng)度高、脆性大、硬度高等特點(diǎn)。若深部地下工程處于花崗巖地層中，在施工過程中巖爆、板裂、層裂、片幫等工程地質(zhì)災(zāi)害頻發(fā)，將對工程的進(jìn)度、成本以及人員的安全造成極大的影響和危害[1-2]。因此，正確認(rèn)識(shí)硬脆巖在不同力學(xué)狀態(tài)下的變形破壞規(guī)律對預(yù)測和防治巖爆等地質(zhì)災(zāi)害具有積極作用。此外，研究不同應(yīng)力狀態(tài)下巖體的力學(xué)性質(zhì)對正確認(rèn)識(shí)巖體破壞機(jī)制、準(zhǔn)確分析圍巖穩(wěn)定性以及合理制定開挖方案和支護(hù)措施有著重要的意義[3]。

國內(nèi)外學(xué)者對硬巖在不同力學(xué)環(huán)境下的巖石力學(xué)性質(zhì)做了大量的研究工作，并取得了許多成果。Martin 等[4]利用剛性試驗(yàn)機(jī)對大量脆性巖石進(jìn)行試驗(yàn)，得到了硬巖完整的應(yīng)力-應(yīng)變曲線和破壞形態(tài)；Zong 等[5]為研究圍壓對硬巖的巖石力學(xué)性質(zhì)及損傷演化的影響規(guī)律，對砂巖進(jìn)行了大量三軸試驗(yàn)，通過分析不同應(yīng)力狀態(tài)下砂巖的應(yīng)力-應(yīng)變特征和強(qiáng)度變形特征，得出了峰值強(qiáng)度、彈性模量和變形模量隨圍壓呈線性增長的規(guī)律；蒲超等[6]通過對砂巖進(jìn)行三軸壓縮試驗(yàn)，定量分析了圍壓對硬巖強(qiáng)度和變形特征的影響規(guī)律，并基于Mohr-Coulomb 強(qiáng)度準(zhǔn)則求解出砂巖的巖石力學(xué)參數(shù)。上述研究對硬巖巖石力學(xué)的發(fā)展起到了極大的促進(jìn)作用，但受載巖體破壞的本質(zhì)是能量積聚和釋放，因此研究不同力學(xué)環(huán)境下巖石破壞過程中的能量演化規(guī)律對進(jìn)一步認(rèn)識(shí)巖石力學(xué)性質(zhì)和揭示受載巖石破壞機(jī)理具有重要的意義。為此，眾多學(xué)者對硬巖在不同力學(xué)環(huán)境下變形破壞過程中的能量演化機(jī)制開展了深入和系統(tǒng)的研究。Mikhalyuk 等[7]在不同侵入巖實(shí)驗(yàn)的基礎(chǔ)上，建立了動(dòng)態(tài)變形能量隨加載強(qiáng)度和加載條件變化的一般機(jī)理，并比較了準(zhǔn)彈性和彈塑性變形區(qū)動(dòng)態(tài)加載能耗散的定性數(shù)據(jù)；謝和平等[8]通過大量單軸、三軸試驗(yàn)，分析了巖石在不同力學(xué)狀態(tài)下破壞過程中的能量演化機(jī)制和能量對巖石強(qiáng)度的影響規(guī)律，并建立了變形破壞過程中巖石能量耗散、釋放與巖石強(qiáng)度和整體破壞情況的內(nèi)在聯(lián)系；田勇等[9]通過對灰?guī)r巖樣進(jìn)行三軸壓縮試驗(yàn)，揭示了巖石在壓縮過程中各階段的能量轉(zhuǎn)化方式；于水生等[10]通過三軸壓縮試驗(yàn)分析了花崗巖能量特征與應(yīng)力、應(yīng)變及圍壓之間的關(guān)系；Li 等[11]對動(dòng)荷載作用下巖石破裂過程中的能量特征進(jìn)行了定量分析，發(fā)現(xiàn)能量與細(xì)觀形貌特征之間存在對應(yīng)關(guān)系，即破壞過程中能量吸收越多，巖石碎裂程度越明顯。

本研究將通過室內(nèi)三軸壓縮試驗(yàn)與單軸壓縮試驗(yàn)對陜西省寶雞市某隧道花崗巖的力學(xué)性質(zhì)與受載破壞過程中的能量演化規(guī)律進(jìn)行分析，探明不同圍壓下花崗巖在受載破壞過程中的能量演化過程，探索花崗巖的能量演化規(guī)律與力學(xué)破壞機(jī)制之間的關(guān)系，從而揭示花崗巖在加載破壞時(shí)的能量演化規(guī)律，為深部地下工程巖體力學(xué)參數(shù)的選取、工程設(shè)計(jì)計(jì)算以及開挖、支護(hù)方案的選取提供理論依據(jù)。

1 試驗(yàn)設(shè)備和方法

1.1 試樣及試驗(yàn)設(shè)備

本研究的巖樣取自陜西省寶雞市某隧道，巖樣主要成分為微風(fēng)化中粗粒黑白云母花崗巖，巖質(zhì)較硬，巖體較完整且表面光滑，無明顯缺陷。為了減少試驗(yàn)結(jié)果的離散性，盡量采用從一個(gè)大巖塊上密集套鉆的方法獲取巖石試件，鉆孔取芯后將試件裝入巖芯箱，運(yùn)送至實(shí)驗(yàn)室。嚴(yán)格按照標(biāo)準(zhǔn)，將巖樣加工成直徑50 mm、高100 mm 的圓柱體，且試件兩端面的不平整度不大于0.5 mm，試件的直徑和高度誤差不大于0.3 mm，兩端面垂直于試件軸線且最大偏差角不大于0.25°，巖石縱波波速在4 km/s 左右，平均容重為2.41 g/cm3。

如圖1 所示，本試驗(yàn)是在RMT-150B 巖石力學(xué)試驗(yàn)系統(tǒng)上進(jìn)行的，該系統(tǒng)最高可以施加50 MPa的圍壓，最大軸向力1 000 kN，圍壓加載速率范圍為0.001～1 MPa/s，最大壓縮變形量為20 mm，可控制變形率范圍為0.001～1 mm/s，可控制加載率范圍為0.01～100 kN/s。

圖1 RMT-150B 巖石力學(xué)試驗(yàn)系統(tǒng)Fig. 1 RMT-150B rock mechanics test system

1.2 試驗(yàn)方法

三軸壓縮試驗(yàn)是在RMT-150B 巖石力學(xué)試驗(yàn)系統(tǒng)上進(jìn)行的。在加載之前先將巖樣用乳膠套包裹好，然后在花崗巖巖樣的上下兩端分別墊上直徑50 mm、高度25 mm 的剛性墊塊，最后將巖樣放入壓力缸中進(jìn)行加載。以靜水壓力狀態(tài)加載巖樣到預(yù)定的圍壓值，然后以恒定位移速率施加軸壓直至巖樣被破壞，設(shè)置0、5、10、15 和20 MPa 共5 個(gè)圍壓水平，每組應(yīng)變率下有3 個(gè)巖樣，共15 個(gè)巖樣，每組各取一個(gè)典型巖樣進(jìn)行分析。

2 不同圍壓下花崗巖力學(xué)性質(zhì)試驗(yàn)結(jié)果及分析

2.1 不同圍壓下花崗巖的應(yīng)力-應(yīng)變曲線

不同圍壓下花崗巖的應(yīng)力-應(yīng)變關(guān)系曲線如圖2 所示，可以看出5 種不同圍壓下巖石的應(yīng)力-應(yīng)變曲線特征大致相同，分為4 個(gè)階段：壓密階段、線彈性變形階段、塑性變形階段和破壞階段。圍壓不同時(shí)，花崗巖峰前階段的應(yīng)力-應(yīng)變曲線變化趨勢基本相同。本次試驗(yàn)采用的花崗巖具有明顯的硬脆性，因此在三軸加載峰前階段，隨著軸向荷載增大，應(yīng)力-應(yīng)變曲線呈直線或類似直線直至破壞，表現(xiàn)出較好的彈性特性。從圖2 可以看出，圍壓越大，應(yīng)力-應(yīng)變曲線峰前階段的彈性模量越大，但其增幅不大，說明在0～20 MPa 的高圍壓下花崗巖表現(xiàn)出一定的壓硬性。在低圍壓情況下，由于巖石材料的非均質(zhì)性，強(qiáng)度低的材料在巖石加載屈服階段先達(dá)到其極限承載力而產(chǎn)生塑性變形，而強(qiáng)度較高的材料達(dá)到其應(yīng)力峰值時(shí)還未發(fā)生破壞。隨著材料變形以及損傷完全發(fā)育，材料的軸向承載能力逐漸減弱，而巖樣將在已經(jīng)發(fā)生塑性變形的區(qū)域產(chǎn)生應(yīng)力集中，最終導(dǎo)致變形局部化，巖樣發(fā)生脆性破壞。在圍壓較高的情況下，同樣強(qiáng)度較低的材料先達(dá)到其極限應(yīng)力而產(chǎn)生塑性變形，但由于圍壓效應(yīng)，巖石材料的強(qiáng)度大幅度增強(qiáng)，因此在峰值應(yīng)力附近有較強(qiáng)的承載能力，所以峰值階段的巖樣還有一定的承載能力；但隨著變形逐漸增加，巖樣內(nèi)部損傷不斷發(fā)育而導(dǎo)致巖樣承載能力逐漸降低，最終巖樣達(dá)到其承載極限而發(fā)生屈服破壞，巖樣表現(xiàn)為延塑性破壞。在高圍壓情況下，峰值階段花崗巖出現(xiàn)屈服平臺(tái)，塑性變形明顯增加，與Yang 等[1]的試驗(yàn)結(jié)果相同。此外，隨著圍壓增大，花崗巖的峰值應(yīng)力和峰值應(yīng)變均逐漸增大，在峰后階段，圍壓越大，應(yīng)力衰減速率也越快。

圖2 不同圍壓下花崗巖的應(yīng)力-應(yīng)變曲線Fig. 2 Stress-strain curves of granite under different confining pressures

2.2 不同圍壓下花崗巖峰值的應(yīng)力-應(yīng)變特性

表1 顯示了不同圍壓下花崗巖的峰值應(yīng)力和峰值應(yīng)變。圍壓為5 MPa 時(shí)，花崗巖的峰值應(yīng)力比單軸壓縮時(shí)的峰值應(yīng)力增加70.92%，峰值應(yīng)變增加17.48%；當(dāng)圍壓從5 MPa 增加到10 MPa 時(shí)，巖樣的抗壓強(qiáng)度增加了約21.54%，峰值應(yīng)變增加了10.72%；當(dāng)圍壓從10 MPa 增加到15 MPa 時(shí)，花崗巖的抗壓強(qiáng)度增加了14.39%，峰值應(yīng)變增加了14.43%；當(dāng)圍壓從15 MPa 增加到20 MPa 時(shí)，巖樣的峰值應(yīng)力增加了16.43%，峰值應(yīng)變增加了13.21%。圖3 為不同圍壓下峰值應(yīng)力、應(yīng)變與圍壓的關(guān)系，可以看出峰值應(yīng)力、峰值應(yīng)變和圍壓具有較強(qiáng)的線性關(guān)系?；贛ohr-Coulomb 強(qiáng)度準(zhǔn)則，對試驗(yàn)結(jié)果擬合得到圍壓-峰值應(yīng)力曲線，如圖3 所示。

式中： σ1為最大主應(yīng)力， σ3為最小主應(yīng)力； σc為花崗巖的理論單軸抗壓強(qiáng)度； ξ為圍壓對花崗巖承載能力 的影響系數(shù)。

式中： φ為巖石材料的內(nèi)摩擦角，c 為巖石材料的黏聚力。通過式（1）、式（2）和式（3）計(jì)算得到花崗巖巖樣 的內(nèi)摩擦角 φ為57.629°，黏聚力c 為23.548 MPa。

表1 不同圍壓下花崗巖的力學(xué)性質(zhì)參數(shù)Table 1 Mechanical properties of granite under different confining pressures

2.3 不同圍壓下花崗巖的破壞模式

圖3 圍壓與峰值應(yīng)力、應(yīng)變的關(guān)系Fig. 3 Relationship of peak stress and strain with confining pressure

不同圍壓下花崗巖試件的典型加載破裂特征如圖4 所示。在加載破壞過程中巖樣呈現(xiàn)微鼓狀，產(chǎn)生了明顯沿橫向擴(kuò)容現(xiàn)象；在破壞時(shí)，花崗巖試件發(fā)出清脆的爆裂聲，體現(xiàn)了巖樣的脆性破壞特征。

從圖4 可以看出，花崗巖的破壞形式與圍壓有顯著的相關(guān)性，在三軸壓縮時(shí)花崗巖試樣的主要破壞形式為劈裂破壞和剪切破壞，且圍壓越低，劈裂破壞越完全。當(dāng)圍壓為5 MPa 時(shí)，如圖4（a）所示，巖樣表面只出現(xiàn)兩條主劈裂裂紋；當(dāng)圍壓為10 MPa 時(shí)，巖樣表面可見兩條平行裂紋，且?guī)r樣表現(xiàn)為張剪破壞，其端部損傷嚴(yán)重；當(dāng)圍壓為15 MPa 時(shí)，巖樣表現(xiàn)為典型的Y 型共軛剪切；當(dāng)圍壓為20 MPa 時(shí)，花崗巖主要表現(xiàn)為剪切破壞，在巖樣的上端、剪切裂紋的起始部分出現(xiàn)張性的次生裂紋，巖樣底端破碎較嚴(yán)重，并出現(xiàn)巖粉。由于圍壓較高時(shí)，巖樣內(nèi)部的礦物顆粒相互黏結(jié)、嵌合較緊密，導(dǎo)致巖樣的黏聚力提高，因此不容易發(fā)生斷裂。

圖4 不同圍壓下花崗巖巖樣的破壞特征Fig. 4 Failure characteristics of granite samples under different confining pressures

3 不同圍壓下花崗巖的能量演化規(guī)律

巖石加載破壞是一個(gè)伴隨著能量吸收、儲(chǔ)存、耗散以及釋放的過程。如不考慮外部溫度變化產(chǎn)生的熱能，外力對巖石系統(tǒng)做功一部分以彈性應(yīng)變能的形式儲(chǔ)存，另一部分則以塑性變形能、損傷能等耗散能的形式耗散；當(dāng)巖石從外部系統(tǒng)吸收并儲(chǔ)存的彈性應(yīng)變達(dá)到其儲(chǔ)能極限時(shí)，巖石發(fā)生破壞，儲(chǔ)存的能量一部分轉(zhuǎn)化為表面能用于巖石破裂損傷，而多余的能量以巖塊動(dòng)能、熱能、各種輻射能等形式向外界釋放。研究巖石在加載下變形破壞過程的應(yīng)變能、彈性應(yīng)變能和耗散能的演化規(guī)律，從能量角度揭示巖石的力學(xué)特性和損傷機(jī)制，對于認(rèn)識(shí)巖石的損傷演化具有重要意義。

3.1 能量分析方法

當(dāng)試驗(yàn)機(jī)對巖樣做功產(chǎn)生變形時(shí)，假定巖石系統(tǒng)為封閉系統(tǒng)，與外界無熱量交換，由熱力學(xué)第一定律可知，巖石系統(tǒng)從外界吸收的能量可轉(zhuǎn)化為可釋放的彈性應(yīng)變能和耗散能[12]

式中：W 為試驗(yàn)機(jī)對巖樣做的總功，U 為總能量，Ue為彈性應(yīng)變能，Ud為耗散應(yīng)變能。圖5 展示了巖樣加載破壞過程中彈性應(yīng)變能Ue和耗散應(yīng)變能Ud的關(guān)系，其中： σi為主應(yīng)力， εi為主應(yīng)變，Ei為巖石卸荷彈 性模量。本試驗(yàn)為常規(guī)三軸壓縮試驗(yàn)，即 σ2= σ3，總能量

式中： σ2為中間主應(yīng)力； ε1、 ε2、 ε3分 別為最大主應(yīng)力、中間主應(yīng)力和最小主應(yīng)力方向的應(yīng)變； μ為巖石的 泊松比；Ei為巖石卸荷彈性模量，計(jì)算時(shí)可用初始彈性模量E0代替。

3.2 不同圍壓下的能量演化過程

圖5 巖石壓縮過程中Ud 和Ue 的關(guān)系Fig. 5 Relationship between Ud and Ue in rock failure

圖6 給出了每立方米花崗巖在不同圍壓下受載破壞過程中的能量演化特征，其中， σci和 σcd分別為起裂應(yīng)力和擴(kuò)容應(yīng)力。根據(jù)不同圍壓下花崗巖的能量變化規(guī)律，將其應(yīng)力-應(yīng)變曲線分為初始?jí)好茈A段、彈性階段、微裂紋穩(wěn)定和加速擴(kuò)展階段以及峰后破壞階段[4]。

圖6 不同圍壓下花崗巖的能量演化特征Fig. 6 Energy evolution characteristics of granite under different confining pressures

在初始?jí)好茈A段，巖樣從外界吸收的能量基本轉(zhuǎn)化為使巖石內(nèi)部微裂隙發(fā)生閉合和摩擦滑移的耗散能Ud，巖樣內(nèi)部基本沒有儲(chǔ)存彈性應(yīng)變能。在線彈性變形階段，巖樣的原生裂紋已經(jīng)被壓密，巖石從外界吸收的能量主要以彈性應(yīng)變能的形式儲(chǔ)存在巖樣中，巖石內(nèi)部儲(chǔ)存的彈性應(yīng)變能隨著彈性變形的發(fā)生而增大，但儲(chǔ)能速率恒定，彈性應(yīng)變能呈直線型。在微裂紋穩(wěn)定與加速擴(kuò)展階段，巖樣內(nèi)部萌生大量的微裂紋，隨著裂紋的發(fā)展、匯合，巖樣變形增加。隨著軸向應(yīng)力不斷增加，巖樣內(nèi)部微裂紋穩(wěn)定擴(kuò)展，此時(shí)耗散能開始增加，彈性應(yīng)變能繼續(xù)保持增長但增長的速率平緩。隨著軸向荷載繼續(xù)增大，巖樣內(nèi)部微裂紋開始加速擴(kuò)展，此時(shí)耗散能的增長速率逐漸變大，而彈性能增長速率急劇減小，說明巖石由于能量耗散的作用，結(jié)構(gòu)發(fā)生了大的變化。在此階段圍壓對加載過程中能量演化規(guī)律的影響顯著，圍壓越大，巖樣在此階段儲(chǔ)存的彈性應(yīng)變能越大，產(chǎn)生的耗散能與總能量的比越小，表明在圍壓高的情況下圍壓效應(yīng)會(huì)使巖樣內(nèi)部的損傷發(fā)育不完全。處于破壞階段時(shí)，由于施加了軸向荷載，巖樣內(nèi)部微裂紋匯合、貫通，形成宏觀破裂面，彈性應(yīng)變能轉(zhuǎn)化為耗散能，然后以表面能、動(dòng)能等形式釋放，彈性能密度急劇大幅減小，耗散能密度大幅增加，隨著破壞過程的進(jìn)行，彈性應(yīng)變能處于較低的水平，最終因整體構(gòu)架強(qiáng)度喪失而導(dǎo)致巖石破壞。

從圖6 可以看出，在微裂紋穩(wěn)定擴(kuò)展階段巖樣儲(chǔ)存彈性應(yīng)變能的速率較均勻，彈性能表現(xiàn)為直線形式；而在巖樣進(jìn)入微裂紋加速擴(kuò)展階段后，花崗巖彈性應(yīng)變能儲(chǔ)存的速率大幅減小。因此根據(jù)花崗巖在微裂紋穩(wěn)定擴(kuò)展階段的線性儲(chǔ)能規(guī)律可以得出巖樣的應(yīng)力閾值 σci和 σcd。由表2 和圖7 可知，圍壓 越大，起裂應(yīng)力、擴(kuò)容應(yīng)力越大，巖樣起裂點(diǎn)處與擴(kuò)容點(diǎn)處的各能量也越大。

表2 不同圍壓下的應(yīng)力閾值及能量Table 2 Stress threshold and energy under different confining pressures

圖7 不同圍壓下的應(yīng)力閾值及能量Fig. 7 Stress threshold and energy under different confining pressures

3.3 不同圍壓下的峰值能量特征

深部地下工程在高地應(yīng)力環(huán)境下施工必然會(huì)導(dǎo)致彈性應(yīng)變能釋放和耗散能耗散。彈性應(yīng)變能的轉(zhuǎn)化和耗散均會(huì)使巖石產(chǎn)生損傷而強(qiáng)度降低，彈性應(yīng)變能釋放導(dǎo)致了巖石破壞，剩余的能量轉(zhuǎn)化為巖石的動(dòng)能，最終導(dǎo)致地下工程在施工過程中可能會(huì)引發(fā)一系列地質(zhì)災(zāi)害，如大變形、巖爆等，因此需要對巖石在不同力學(xué)環(huán)境下的能量轉(zhuǎn)化機(jī)制進(jìn)行深入和系統(tǒng)的研究。圖8 和表3 為不同圍壓下每立方米巖樣的總能量、彈性應(yīng)變能和耗散能密度。

由圖8 和表3 可知，總能量、彈性應(yīng)變能和耗散能隨著圍壓的增大而增大，且具有很強(qiáng)的正線性，圖8 還顯示了對試驗(yàn)結(jié)果進(jìn)行線性擬合得到的圍壓與各能量之間的關(guān)系。低圍壓情況下，巖石的承載能力與峰值應(yīng)變都較小，巖樣在破壞前所吸收、儲(chǔ)存的能量相應(yīng)也比較小。從表3中可知，當(dāng)圍壓較大時(shí)，峰值點(diǎn)處的耗散能占總能量的33.5%，彈性應(yīng)變能的密度達(dá)到1.099 MJ/m3。由2.1 節(jié)可知，在高圍壓情況下巖樣還具有較高的承載能力，在荷載作用下發(fā)生變形、損傷并逐漸發(fā)育，此時(shí)巖樣從外界系統(tǒng)吸收并儲(chǔ)存了大量的彈性應(yīng)變能：一部分轉(zhuǎn)化成耗散能，在耗散能的作用下材料顆粒間的黏結(jié)力不斷喪失，微裂隙不斷發(fā)育，巖石強(qiáng)度不斷喪失，因此巖石釋放的耗散能占總能量的比例較大，此時(shí)巖石表現(xiàn)出延塑性特征；另一部分則以彈性應(yīng)變能的形式儲(chǔ)存在巖石內(nèi)部，當(dāng)巖石破壞時(shí)轉(zhuǎn)化為耗散能，并以動(dòng)能、熱能和輻射能等形式釋放。在實(shí)際工程中，在高地應(yīng)力區(qū)域內(nèi)施工極有可能產(chǎn)生巖塊彈射并發(fā)生巖爆，因此需要采用應(yīng)力釋放施工措施降低巖爆發(fā)生的概率。

圖8 圍壓與峰值點(diǎn)處能量的關(guān)系Fig. 8 The relationship between confining pressure and energy at peak point

表3 不同圍壓下峰值點(diǎn)處的能量Table 3 Energy at peak points under different confining pressures

4 基于能量耗散的巖樣損傷模型

巖樣的破壞實(shí)質(zhì)上是在能量耗散作用下發(fā)生損傷并產(chǎn)生失穩(wěn)的一種現(xiàn)象，因巖樣加載過程中能量轉(zhuǎn)化規(guī)律對巖石的變形破壞影響顯著，其變形破壞過程中的能量演化規(guī)律能夠真實(shí)反映巖體破壞的本質(zhì)特征，可將巖體單元的能量損傷量定義為[8]

式中：Ud為巖樣某一時(shí)刻所產(chǎn)生的耗散能，Udmax為巖石完全破壞時(shí)需要的耗散能。根據(jù)Kachanov 對損傷變量的定義[13-14]，得到基于能量轉(zhuǎn)化原理的巖石損傷變量

式中： σp為巖石的峰值強(qiáng)度， σr為巖石的殘余強(qiáng)度。由式（1）和式（8）得到考慮圍壓效應(yīng)的巖石損傷變 量

根據(jù)彈塑性力學(xué)以及相關(guān)的研究成果，巖石在三軸應(yīng)力狀態(tài)下的本構(gòu)關(guān)系

式中：E 為彈性模量?？紤]能量耗散對巖石材料的影響，根據(jù)式（8）和式（10）得到基于耗散能的損傷演化模型

根據(jù)式（9）得出加載過程中花崗巖在不同圍壓下?lián)p傷變量的演化規(guī)律，如圖9 所示。在加載初期的壓密階段，巖樣內(nèi)部初始裂隙被壓密導(dǎo)致部分能量耗散，此時(shí)巖樣應(yīng)力水平較低，因此損傷發(fā)育緩慢；隨著荷載的增加，巖石內(nèi)部微裂紋開始穩(wěn)定擴(kuò)展，而新裂紋的萌生與擴(kuò)展需要消耗大量能量，此時(shí)巖石的損傷變量持續(xù)增大；隨著軸向荷載的增大，巖樣內(nèi)部的微裂紋開始加速擴(kuò)展、匯合，最終形成宏觀裂紋，巖樣能量急劇耗散，應(yīng)力水平較高，巖樣的損傷變量在該階段快速增大直至巖樣破壞。

如圖9 所示，損傷突變點(diǎn)均在峰值應(yīng)力附近，但圍壓較小時(shí)，在損傷突變點(diǎn)前巖樣損傷表現(xiàn)出明顯的突變型，表明此時(shí)巖石內(nèi)部損傷發(fā)育未完全就發(fā)生了破壞，使得大量的裂紋萌生在巖樣破壞前后，此時(shí)易發(fā)生劈裂破壞；當(dāng)圍壓較高時(shí)，損傷突變點(diǎn)之前的損傷緩慢增加，圍壓效應(yīng)使得巖石材料的強(qiáng)度大幅提高，而巖石在較高的應(yīng)力下產(chǎn)生變形，巖樣內(nèi)部損傷不斷發(fā)育，最終達(dá)到其承載極限而屈服破壞，巖樣最終表現(xiàn)為剪切破壞。

圖9 不同圍壓下的花崗巖損傷曲線Fig. 9 Damage curves of granite under different confining pressures

5 結(jié) 論

基于花崗巖試樣的三軸壓縮試驗(yàn)結(jié)果，研究了不同圍壓下花崗巖巖石的力學(xué)性質(zhì)及能量演化特征，得到如下結(jié)論。

（1）本試驗(yàn)所用的花崗巖具有明顯的脆性，圍壓對花崗巖的抗壓強(qiáng)度、峰值應(yīng)變影響顯著，且?guī)r樣的峰值應(yīng)力與峰值應(yīng)變和圍壓具有較強(qiáng)的正線性關(guān)系，基于Mohr-Coulomb 強(qiáng)度準(zhǔn)則，求出花崗巖的黏聚力為23.548 MPa，內(nèi)摩擦角為57.629°。

（2）巖樣的破壞模式與圍壓有較強(qiáng)相關(guān)性，當(dāng)圍壓從5 MPa 增大到20 MPa 時(shí)，巖石表現(xiàn)為典型劈裂破壞、張剪破壞、Y 型共軛剪切、典型剪切破壞，這是由圍壓效應(yīng)對巖樣摩擦強(qiáng)度與黏聚強(qiáng)度的影響導(dǎo)致的。

（3）根據(jù)花崗巖線性儲(chǔ)能規(guī)律，得出花崗巖的應(yīng)力閾值與閾值點(diǎn)處能量，圍壓越大，起裂應(yīng)力、擴(kuò)容應(yīng)力越大，且?guī)r樣起裂點(diǎn)處與擴(kuò)容點(diǎn)處的各能量也越大；巖石在加載破壞過程中吸收的總能量、儲(chǔ)存的彈性應(yīng)變能以及耗散能都隨著圍壓的增大而線性增大。

（4）基于花崗巖在加載破壞過程中的能量耗散機(jī)制，提出了巖石損傷演化模型，得到了花崗巖的損傷變量D 在不同圍壓下加載破壞過程中的演化規(guī)律，發(fā)現(xiàn)在低圍壓情況下巖石損傷具有突變性，高圍壓情況下在損傷突變點(diǎn)之前損傷緩慢增加。