盧高明 李元輝 張希巍 金長宇

(1.深部金屬礦山安全開采教育部重點(diǎn)實驗室，遼寧 沈陽110819;2.東北大學(xué)資源與土木工程學(xué)院，遼寧 沈陽110819)

大量工程實踐表明，地下工程巖體多處于復(fù)雜的應(yīng)力狀態(tài)，在工程巖體施工過程中，巖石通常處于反復(fù)加卸載過程，巖石的強(qiáng)度、變形及斷裂損傷力學(xué)特性與所受的應(yīng)力狀態(tài)以及承受的荷載形式密切相關(guān)［1-3］。因此研究不同應(yīng)力狀態(tài)循環(huán)加卸載條件下巖石的強(qiáng)度、變形和損傷力學(xué)特性具有重要的理論意義和工程價值。國內(nèi)外許多學(xué)者對循環(huán)載荷破壞過程中巖石的變形破壞特征、力學(xué)特性、破壞機(jī)制及斷裂損傷力學(xué)特性進(jìn)行了相關(guān)的研究［4-7］。E. Eberhardt等［8］對大量花崗巖開展了循環(huán)加卸載試驗，研究了循環(huán)荷載條件下脆性巖石微裂紋擴(kuò)展條件、損傷力學(xué)特性及其斷裂準(zhǔn)則。J. Xiao 等［9］通過花崗巖循環(huán)加卸載試驗對多種損傷變量優(yōu)缺點(diǎn)進(jìn)行了評價，結(jié)果表明這些損傷變量能夠真實體現(xiàn)巖石疲勞損傷的演化過程。循環(huán)加、卸載試驗對巖石在壓縮條件下的變形和損傷特性的定量分析提供了有效技術(shù)手段。周家文等［10］以向家壩砂巖單軸循環(huán)加卸載室內(nèi)力學(xué)試驗結(jié)果為基礎(chǔ)，對脆性巖石峰值強(qiáng)度、彈性常數(shù)變化規(guī)律及其損傷力學(xué)特性進(jìn)行了研究，給出一種根據(jù)應(yīng)力－應(yīng)變曲線計算損傷變量的方法，并驗證了該方法與聲發(fā)射測試結(jié)果較為一致。此外，C. D. Martin 等［11］對花崗巖進(jìn)行了不同應(yīng)力狀態(tài)下的循環(huán)荷載試驗，初步確立了損傷變量與內(nèi)摩擦角、黏聚力的關(guān)系。E.Eberhardt 等［12］根據(jù)C.D.Martin 等人的研究成果，又對花崗巖峰前損傷對其峰值強(qiáng)度和變形的影響進(jìn)行了研究。

上述學(xué)者通過循環(huán)加卸載試驗對脆性巖石的強(qiáng)度、變形、彈性常量及損傷特性等進(jìn)行了定量化的研究，并取得了較為顯著的成果，研究不同應(yīng)力狀態(tài)和不同荷載承受方式下巖石的斷裂損傷機(jī)制具有重要意義。本研究在總結(jié)前人研究成果的基礎(chǔ)上，對馬城鐵礦輝綠巖開展了4 種不同應(yīng)力路徑巖石力學(xué)試驗，并提出一種巖石循環(huán)加卸載過程中相對應(yīng)變損傷參數(shù)，結(jié)合現(xiàn)有損傷力學(xué)理論，對比分析了單軸和三軸條件下絕對應(yīng)變損傷參數(shù)和相對應(yīng)變損傷參數(shù)的變化規(guī)律，對循環(huán)荷載條件下脆性巖石的損傷破壞演化機(jī)制進(jìn)行了研究。

1 巖石試件、試驗儀器及試驗方法

1.1 試件的制備與選擇

試驗試樣為馬城鐵礦輝綠巖，按照國際巖石力學(xué)建議的方法，把取得的巖石試樣加工成直徑為50 mm，高為100 mm 的圓柱型標(biāo)準(zhǔn)試件，其加工精度符合國際巖石力學(xué)學(xué)會建議的實驗室規(guī)范。在自然風(fēng)干狀態(tài)下對試件進(jìn)行波速和密度測試，測試結(jié)果如圖1 所示。測得的P 波波速為4 000 ～6 500 m/s，密度為2.65 ～2.95 g/cm3，為了減小試驗結(jié)果的離散性，試驗所選試件主要集中在圖1 所示橢圓形區(qū)域內(nèi)。

1.2 試驗儀器

試驗所采用的儀器是Rockman207 硬巖三軸試驗系統(tǒng)［13］。該試驗系統(tǒng)包括常規(guī)軸向加載系統(tǒng)、流變軸向加載系統(tǒng)和圍壓加載系統(tǒng)，能夠完成復(fù)雜路徑下的剪切和流變試驗。試驗系統(tǒng)具有3 個獨(dú)立的LVDT測量和反饋通道，以軸向或徑向LVDT信號作為反饋控制。在常規(guī)三軸試驗?zāi)Ｊ较?，高剛性的試驗機(jī)支架、高性能的伺服閥和準(zhǔn)確的徑向變形測量，由計算機(jī)通過變形和負(fù)荷準(zhǔn)確地控制試驗的全過程，為獲得巖石試件的破壞全過程應(yīng)力－應(yīng)變曲線提供了條件，并提供可信度較高的試驗數(shù)據(jù)。

圖1 輝綠巖P 波波速－密度關(guān)系Fig.1 Relationship between P-wave velocity and density of diabase

1.3 試驗方法

(1)單軸壓縮變形試驗，為確定馬城鐵礦輝綠巖單軸抗壓強(qiáng)度和基本巖石力學(xué)參數(shù)，為后續(xù)設(shè)計三軸試驗應(yīng)力路徑提供依據(jù)。軸壓加載速率1 MPa/s，環(huán)向變形0.03 mm/min。

(2)三軸壓縮變形試驗，研究不同圍壓下輝綠巖全破壞過程力學(xué)行為。圍壓加載速率為2 MPa/min，軸壓加載速率為1 MPa/s，環(huán)向變形控為0.04 mm/min。

(3)單軸、三軸不斷增大荷載的循環(huán)加卸載試驗，每次循環(huán)施加的最大荷載比前一次循環(huán)施加的最大荷載為大，分別約為平均強(qiáng)度的20%、40%、60%、80%進(jìn)行卸載，每次卸載至0 MPa，最后加載至試件破壞。

2 試驗結(jié)果和分析

2.1 單軸壓縮應(yīng)力－應(yīng)變曲線

圖2 為一典型單軸壓縮全過程應(yīng)力－應(yīng)變關(guān)系曲線。從圖2 可以清晰地看到巖石變形的4 個階段(即裂隙的壓密、穩(wěn)定發(fā)展、不穩(wěn)定發(fā)展及破壞后階段)，其破壞形式屬于典型的II 型破壞。在壓密階段，曲線呈明顯非線性，試件體積壓縮顯著，而徑向變形較小，說明該巖石內(nèi)部節(jié)理裂隙發(fā)展較充分。屈服點(diǎn)以后，微破裂的發(fā)展發(fā)生了質(zhì)的變化，試件由體積壓縮轉(zhuǎn)為擴(kuò)容，且擴(kuò)容出現(xiàn)在峰值強(qiáng)度前。達(dá)到峰值強(qiáng)度后，試件表面出現(xiàn)明顯微裂紋，且微裂紋的迅速發(fā)展，在此階段，徑向變形迅速增加，體積變形擴(kuò)容迅速，軸向變形有減小趨勢。破壞后的試件仍保持整體狀態(tài)，試件大部分呈典型單斜面剪切破壞，表現(xiàn)出明顯的脆性破壞特征。

圖2 輝綠巖單軸壓縮應(yīng)力－應(yīng)變關(guān)系曲線Fig.2 Stress-strain curves of diabase under uniaxial compression

2.2 循環(huán)加卸載強(qiáng)度分析

馬城鐵礦輝綠巖試件不同圍壓下(0、10、20、30 MPa)單調(diào)壓縮和循環(huán)荷載峰值強(qiáng)度對比關(guān)系如圖3所示。由圖3 可知，在單軸條件下，循環(huán)加卸載平均峰值強(qiáng)度比單調(diào)荷載低10% ～20%，出現(xiàn)單軸循環(huán)加卸載脆性巖石峰值強(qiáng)度的“弱化”現(xiàn)象。這可能跟巖石內(nèi)部微裂隙的斷裂擴(kuò)展力學(xué)特性有關(guān)，在達(dá)到峰值強(qiáng)度之前，局部孔隙的閉合、裂隙的萌生和擴(kuò)展就已經(jīng)影響了巖石力學(xué)性能，巖石內(nèi)部損傷的積累程度隨循環(huán)次數(shù)的增加而增加，導(dǎo)致峰值強(qiáng)度降低。三軸條件下，峰值強(qiáng)度出現(xiàn)離散性，這可能是由于圍壓作用約束了循環(huán)加卸載過程中巖石內(nèi)部微裂隙的斷裂擴(kuò)展，抑制了巖石內(nèi)部損傷的積累程度。對于脆性巖石循環(huán)加卸載的峰值強(qiáng)度來說，與循環(huán)加卸載的次數(shù)、加載方式、巖石內(nèi)部微觀結(jié)構(gòu)以及圍壓條件等都有密切關(guān)系。

圖3 輝綠巖峰值強(qiáng)度與圍壓的關(guān)系Fig.3 Relationship between peak strength and confining pressure of diabase

2.3 循環(huán)加卸載應(yīng)力－應(yīng)變曲線分析

圖4 為單軸和三軸條件下馬城鐵礦輝綠巖試件不斷增大荷載的循環(huán)加卸載應(yīng)力－應(yīng)變曲線。

圖4 單軸和三軸狀態(tài)下輝綠巖應(yīng)力－應(yīng)變曲線Fig.4 Stress-strain curves of diabase under uniaxial and triaxial compression

由圖4 可知，與單軸循環(huán)加卸載相比，在有圍壓條件下峰值強(qiáng)度和彈性極限顯著增加，峰值前的變形也顯著增大。不管是單軸還是三軸條件下，卸載曲線都不與加載曲線重合，而是低于加載曲線，應(yīng)變也不能恢復(fù)到其下限值。加載曲線和卸載曲線形成塑性滯回環(huán)［14］，隨著循環(huán)次數(shù)的增加，塑性滯回環(huán)的面積逐漸增大，并向應(yīng)變增大的方向“移動”。單軸條件下的應(yīng)力－應(yīng)變曲線呈現(xiàn)明顯的內(nèi)凹現(xiàn)象，體現(xiàn)在切線模量上，卸載過程切線模量逐漸減小，重新加載過程切線模量先增加后減小。三軸條件下該曲線的內(nèi)凹現(xiàn)象得到了“弱化”，這是由于圍壓作用改變了巖石內(nèi)部的微裂紋斷裂擴(kuò)展機(jī)制，圍壓作用使得內(nèi)部微裂紋的擴(kuò)展需要在更大的應(yīng)力條件下才能發(fā)生，這也解釋了圍壓升高導(dǎo)致峰值強(qiáng)度增大的原因。

2.4 循環(huán)加卸載彈性常量分析

在循環(huán)加卸載應(yīng)力－應(yīng)變曲線中，加載曲線直線段大致與卸載曲線的割線平行，一般可將卸載曲線的割線斜率作為彈性模量，按照文獻(xiàn)［15］給出的方法分別計算單軸和三軸條件下循環(huán)加卸載過程中的彈性常量，如圖5 所示。

從圖5(a)可知，單軸和三軸循環(huán)加卸載過程得到的泊松比都呈現(xiàn)逐漸增加的趨勢，這說明隨著循環(huán)次數(shù)的增加，徑向應(yīng)變的增大速度大于軸向應(yīng)變，同時也說明巖石內(nèi)部的微裂隙隨加卸載的進(jìn)行不斷發(fā)生斷裂擴(kuò)展。從單軸和三軸泊松比對比來看，第一次循環(huán)差距較大，隨循環(huán)次數(shù)的增加差距逐漸減小并超過三軸泊松比值，這說明單軸條件下徑向應(yīng)變增大速度大于三軸條件，其實質(zhì)為圍壓作用抑制了巖石內(nèi)部微裂隙的斷裂擴(kuò)展條件，對徑向應(yīng)變影響作用較大。從圖5(b)可知，單軸和三軸循環(huán)加卸載過程得到的彈性模量變化趨勢具有較大差異性，單軸彈性模量隨循環(huán)次數(shù)的增加呈現(xiàn)穩(wěn)定增長趨勢，而三軸彈性模量則隨循環(huán)次數(shù)的增加而逐漸減小。通過對單軸和三軸條件下泊松比和彈性模量的對比分析可知，圍壓對循環(huán)加卸載過程中脆性巖石的彈性常量有較大影響作用。

圖5 單軸和三軸條件下循環(huán)加卸載彈性常量Fig.5 Elastic constants of cyclic loading and unloading under uniaxial and triaxial conditions

2.5 循環(huán)加卸載應(yīng)變損傷參數(shù)分析

通過以上對馬城鐵礦輝綠巖力學(xué)特性的試驗研究可知，脆性巖石循環(huán)加卸載過程其實是損傷［16-17］破壞逐漸累積的過程，按照E. Eberhardt 等［8，12］關(guān)于循環(huán)加卸載過程損傷參數(shù)的定義，其計算公式為

式中，wax，wrad和wvol分別為軸向、徑向和體積應(yīng)變損傷參數(shù);εax(i)，εrad(i)和εvol(i)分別為每次循環(huán)中軸向、徑向和體積永久應(yīng)變;n 為循環(huán)加卸載次數(shù)。

式(1)、式(2)、式(3)反映了循環(huán)加卸載過程每次循環(huán)中永久應(yīng)變與總永久應(yīng)變的關(guān)系，即每次循環(huán)中巖石損傷與積累的損傷和的關(guān)系，但忽略了彈性應(yīng)變的作用，由于循環(huán)加卸載過程中的彈性應(yīng)變對巖石的損傷破壞力學(xué)特性也起著重要作用，故在以上應(yīng)變損傷參數(shù)基礎(chǔ)上，補(bǔ)充一種新的應(yīng)變損傷參數(shù)，結(jié)合原有應(yīng)變損傷參數(shù)，能夠更好地反映巖石循環(huán)加卸載過程中的損傷演化規(guī)律，計算公式為

式中，ε'ax(i)，ε'rad(i)和ε'vol(i)分別為每次循環(huán)中軸向、徑向和體積彈性應(yīng)變。其他文獻(xiàn)中稱第1 種損傷參數(shù)為絕對應(yīng)變損傷參數(shù)，為了區(qū)分2 種應(yīng)變損傷參數(shù)，不妨稱第2 種損傷參數(shù)為相對應(yīng)變損傷參數(shù)。2 種應(yīng)變損傷參數(shù)計算結(jié)果如圖6 和圖7 所示。

圖6 單軸條件下2 種損傷參數(shù)計算結(jié)果Fig.6 Calculation results of two damage parameters under uniaxial conditions

由圖6(a)和圖7(a )可知，單軸和三軸條件下，絕對應(yīng)變損傷參數(shù)變化趨勢基本相同，軸向、徑向和體積絕對應(yīng)變損傷參數(shù)均隨循環(huán)次數(shù)的增加而不斷增大，說明絕對應(yīng)變損傷參數(shù)受圍壓的影響作用不大，并且隨著循環(huán)加卸載過程的進(jìn)行，巖石內(nèi)損傷累積程度逐漸增加。第1次循環(huán)具有較大的應(yīng)變損傷參數(shù)值，隨后每次循環(huán)都有不同程度的增加，說明脆性巖石內(nèi)部微裂隙等缺陷在被壓密和釋放的過程中，有一部分缺陷在卸載完成之后并不能完全被釋放;第1 次加卸載之后，微裂隙的調(diào)整已經(jīng)完成，繼續(xù)進(jìn)行循環(huán)加卸載時，每次循環(huán)都有新的斷裂擴(kuò)展產(chǎn)生，導(dǎo)致巖石內(nèi)部損傷逐漸累積。

圖7 三軸條件下2 種損傷參數(shù)計算結(jié)果Fig.7 Calculation results of two damage parameters under triaxial conditions

由圖6(b)和圖7(b)可知，單軸和三軸條件下，軸向和體積相對應(yīng)變損傷參數(shù)隨循環(huán)次數(shù)的增加而減小，根據(jù)絕對應(yīng)變損傷參數(shù)的變化趨勢已經(jīng)知道，巖石內(nèi)部的損傷積累程度隨循環(huán)次數(shù)的增加而不斷增大，但從相對應(yīng)變損傷參數(shù)計算結(jié)果來看，軸向和體積永久應(yīng)變與總應(yīng)變的比值隨循環(huán)次數(shù)的增加而逐漸降低，這說明彈性應(yīng)變比永久應(yīng)變增長要快。由于巖石內(nèi)部的初始損傷以及加載作用形成的次生損傷，在第1 次循環(huán)時產(chǎn)生了較大的相對應(yīng)變損傷參數(shù)值，但不超過0.3，說明該巖石在破壞之前產(chǎn)生的大部分應(yīng)變?yōu)閺椥詰?yīng)變，并且彈性應(yīng)變隨循環(huán)次數(shù)的增加而不斷增大。從圖6 和圖7 對比結(jié)果可以看出，三軸條件下2 種應(yīng)變損傷參數(shù)值都略低于單軸條件下，說明圍壓對巖石內(nèi)部損傷的積累程度具有一定的約束作用。

3 結(jié) 論

(1)根據(jù)單軸和三軸壓縮變形及循環(huán)加卸載試驗結(jié)果，馬城鐵礦輝綠巖屬于脆性硬巖，破壞形式屬于典型的脆性破壞;單軸條件下循環(huán)加卸載峰值強(qiáng)度比單調(diào)荷載低10% ～20%，三軸條件下峰值強(qiáng)度出現(xiàn)離散性，這與巖石的巖性、加載方式以及循環(huán)次數(shù)等因素相關(guān)。

(2)單軸條件下循環(huán)加卸載過程中應(yīng)力－應(yīng)變曲線呈現(xiàn)明顯內(nèi)凹現(xiàn)象，三軸條件下內(nèi)凹現(xiàn)象得到了“弱化”，體現(xiàn)在切線模量上，這與巖石內(nèi)部微裂隙的閉合、張開、擴(kuò)展等斷裂損傷機(jī)制具有密切聯(lián)系;對于循環(huán)加卸載過程中的彈性常量，圍壓對泊松比的影響不大，對彈性模量有較大影響作用。

(3)本研究給出一種相對應(yīng)變損傷參數(shù)，通過對單軸和三軸條件下2 種應(yīng)變損傷參數(shù)的對比分析可知，絕對應(yīng)變損傷參數(shù)隨循環(huán)次數(shù)的增加而不斷增大，幾乎不受圍壓作用的影響，而相對應(yīng)變損傷參數(shù)受到圍壓的約束作用。

［1］ Eberhardt E，Stead D，Stimpson B，et al. The effect of neighbouring cracks on elliptical crack initiation and propagation in uniaxial and triaxial stress fields［J］. Engineering Fracture Mechanics，1998，59(2):103-115.

［2］ Cai M，Kaiser P，Tasaka Y，et al. Generalized crack initiation and crack damage stress thresholds of brittle rock masses near underground excavations［J］.International Journal of Rock Mechanics and Mining Sciences，2004，41(5):833-847.

［3］ 馮春林，畢忠偉，吳永豐，等. 循環(huán)荷載下巖石臨界強(qiáng)度的測定［J］.金屬礦山，2011(2):13-16.

Feng Chunlin，Bi Zhongwei，Wu Yongfeng，et al. Determination of critical strength of rocks subjected to cycle loading［J］.Metal Mine，2011(2):13-16.

［4］ 章清敘，葛修潤，黃 銘，等. 周期荷載作用下紅砂巖三軸疲勞變形特性試驗研究［J］.巖石力學(xué)與工程學(xué)報，2006，25(3):473-478.

Zhang Qingxu，Ge Xiurun，Huang Ming，et al. Testing study on fatigue deformation law of red-sandstone under triaxial compression with cyclic loading［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(3):473-478.

［5］ Xue L. Damage accumulation and fracture initiation in uncracked ductile solids subject to triaxial loading［J］.International Journal of Solids and Structures，2007，44(16):5163-5181.

［6］ Song H，Zhang H，Kang Y，et al. Damage evolution study of sandstone by cyclic uniaxial test and digital image correlation［J］. Tectonophysics，2013，608:1343-1348.

［7］ 葛修潤，蔣 宇，盧允德，等. 周期荷載作用下巖石疲勞變形特性試驗研究［J］. 巖石力學(xué)與工程學(xué)報，2003，22(10):1581-1585.

Ge Xiurun，Jiang Yu，Lu Yunde，et al. Testing study on fatigue deformation law of rock under cyclic loading［J］. Chinese Journal of Rock Mechanics and Engineering，2003，22(10):1581-1585.

［8］ Eberhardt E，Stead D，Stimpson B，et al. Changes in acoustic event properties with progressive fracture damage［J］.International Journal of Rock Mechanics and Mining Sciences，1997，34(3):71-71.

［9］ Xiao J，Ding D，Jiang F，et al.Fatigue damage variable and evolution of rock subjected to cyclic loading［J］.International Journal of Rock Mechanics and Mining Sciences，2010，47(3):461-468.

［10］ 周家文，楊興國，符文熹，等.脆性巖石單軸循環(huán)加卸載試驗及斷裂損傷力學(xué)特性研究［J］. 巖石力學(xué)與工程學(xué)報，2010，29(6):1172-1183.

Zhou Jiawen，Yang Xingguo，F(xiàn)u Wenxi，et al.Experimental test and fracture damage mechanical characteristics of brittle rock under uniaxial cyclic loading and unloading conditions［J］.Chinese Journal of Rock Mechanics and Engineering，2010，29(6):1172-1183.

［11］ Martin C D，Chandler N A.The progressive fracture of lac du bonnet granite［J］. International Journal of Rock Mechanics and Mining Sciences and Geomechanical Abstracts，1995，31(6):643-659.

［12］ Eberhardt E，Stead D，Stimpson B.Quantifying progressive pre-peak brittle fracture damage in rock during uniaxial compression［J］.International Journal of Rock Mechanics and Mining Sciences，1999，36(3):361-380.

［13］ 張希巍，楊成祥，張 柬，等.紅透山銅礦深部片麻巖力學(xué)行為試驗研究［J］.巖石力學(xué)與工程學(xué)報，2013，32(S2):3228-3237.

Zhang Xiwei，Yang Chengxiang，Zhang Jian，et al. Experimental study of mechanical behavior of deep gneiss in Hongtoushan copper Mine［J］. Chinese Journal of Rock Mechanics and Engineering，2013，32(S2):3228-3237.

［14］ 肖福坤，申志亮，劉 剛，等.循環(huán)加卸載中滯回環(huán)與彈塑性應(yīng)變能關(guān)系研究［J］. 巖石力學(xué)與工程學(xué)報，2014，33(9):1791-1797.

Xiao Fukun，Shen Zhiliang，Liu Gang，et al. Relationship between hysteresis loop and elastoplastic strain energy during cyclic loading and unloading［J］. Chinese Journal of Rock Mechanics and Engineering，2014，33(9):1791-1797.

［15］ 蔡美峰，何滿潮，劉東燕.巖石力學(xué)與工程［M］.北京:科學(xué)出版社，2002:50-58.

Cai Meifeng，He Manchao，Liu Dongyan. Rock Mechanics and Engineering［M］.Beijing:Science Press，2002:50-58.

［16］ 夏 冬，楊天鴻，常 宏.含水煤巖損傷破壞過程中聲發(fā)射特征的研究［J］.金屬礦山，2014(6):1-9.

Xia Dong，Yang Tianhong，Chang Hong.AE properties of moisturecontaining coal and rock during the failure progress［J］. Metal Mine，2014(6):1-9.

［17］ 李旭絨，紀(jì)洪廣，王 軍，等.混凝土井壁受鹽害侵蝕后的強(qiáng)度損傷試驗研究［J］.金屬礦山，2014(10):157-160.

Li Xurong，Ji Hongguang，Wang Jun，et al.Test research on strength damage of concrete shaft lining under salt disaster corrosion［J］.Metal Mine，2014(10):157-160.