脆性巖石蠕變裂紋成核宏細觀力學機理研究1)

李曉照 賈亞星 張騏爍 戚承志,?,

?(北京建筑大學土木與交通工程學院,北京 100044)

?(北京未來城市設計高精尖創(chuàng)新中心,北京 100044)

??(城市交通基礎設施建設北京國際科技合作基地,北京 100044)

引言

隨著我國經(jīng)濟不斷發(fā)展,深部地下工程建設需求與日俱增.深部高地應力環(huán)境導致了地下工程圍巖具有很強的時間效應,脆性巖石又是深部地下工程圍巖的主要組成,進而導致深部地下工程圍巖表現(xiàn)為明顯的脆性蠕變特性,其對深部地下工程圍巖微震及巖爆災害的評價有著重要的實際意義.

眾多學者通過蠕變試驗研究了脆性巖石的蠕變應變演化曲線,并分析了不同恒定軸向加載應力、分級軸向加載應力、分級圍壓卸載等對巖石蠕變應變演化曲線影響[1-2].也有學者研究了不同含水率、滲透壓及溫度對脆性巖石蠕變特性影響[3-6].在以上試驗的基礎上,有學者通過彈性、黏性及塑性模塊進行不同的串并聯(lián)組合,提出了巖石蠕變演化模型,并結(jié)合試驗對比分析校驗模型參數(shù),進而將模型應用于數(shù)值模擬軟件中,分析預測深部地下工程圍巖穩(wěn)定性[7-8].然而,裂紋變化行為嚴重影響固體材料變形力學性能[9-11],以上研究均是從宏觀角度研究了巖石的減速、穩(wěn)態(tài)及加速蠕變應變隨時間演化曲線,無法準確地分析巖石內(nèi)部局部裂紋擴展、貫通及成核現(xiàn)象對巖石蠕變應變演化曲線的影響.

有學者通過聲發(fā)射監(jiān)測技術(shù),對脆性巖石蠕變演化過程中,巖石內(nèi)部裂紋擴展情況進行了更加精確的分析[12-13],圖1(a)中,隨著時間演化,巖石蠕變過程中聲發(fā)射率是不斷變化的,每一次較大的聲發(fā)射率的突增,都對應著一次明顯的應變加速增大到減速增大的轉(zhuǎn)化過度,然而,在較小的聲發(fā)射率變化時,應變值則沒有明顯的改變;圖1(b) 中,累計聲發(fā)射率可以更加明顯地觀測到累計聲發(fā)射率演化趨勢與應變演化趨勢相似,每一次應變加速增大,對應著累計聲發(fā)射率的加速的增大,即巖石的應變率與聲發(fā)射率成正比例關系,此外裂紋速率、應變率及聲發(fā)射率三者之間存在相互正比例對應關系,本文的理論研究與結(jié)果討論中將給出詳細說明.因此通過上述結(jié)果可以看到,脆性巖石蠕變過程中,其蠕變演化曲線并不只是經(jīng)歷一次減速、穩(wěn)態(tài)及加速蠕變演化過程,而是中間可能存在多次應變率突變的情況.該現(xiàn)象與脆性巖石內(nèi)部微震事件的發(fā)生有著密切關系,每一次蠕變應變的加速增大,也就是聲發(fā)射率的突增現(xiàn)象,會導致一次微震現(xiàn)象,多次微震事件的發(fā)生,最終會導致巖石的整體斷裂破壞,甚至會產(chǎn)生嚴重的巖爆事件[14-15].

圖1 巖石蠕變過程中(a)聲發(fā)射信號率及(b)累計聲發(fā)射事件率時間演化關系[12]Fig.1 Evolutions of(a)rate of acoustic emission and(b)the number accumulated acoustic emission during creep of rocks[12]

在恒定荷載蠕變演化過程中,導致這種較大的聲發(fā)射率及應變率的突增與應變值的加速增大到減速增大的過度變化的原因,是由于脆性巖石內(nèi)部局部微裂紋成核損傷突變所導致,蠕變過程中較小的聲發(fā)射率出現(xiàn)只是巖石內(nèi)部裂紋擴展現(xiàn)象導致,此時并沒有局部裂紋成核現(xiàn)象出現(xiàn)[16].此外,基于巖石細觀裂紋擴展模型并結(jié)合裂紋亞臨界裂紋擴展法則,脆性巖石壓縮蠕變宏細觀力學模型被研究[17].然而,目前能夠解釋脆性巖石蠕變過程中,細觀裂紋成核損傷突變對巖石宏觀變形影響的宏細觀力學模型很少,因此本文將基于固體斷裂損傷力學模型,并結(jié)合巖石聲發(fā)射蠕變監(jiān)測試驗結(jié)果,研究分析細觀裂紋成核損傷突變影響下的巖石蠕變演化機理,進而為深部地下工程脆性蠕變過程中微震、巖爆事件的預測評價提供一定的理論支持.

1 理論模型

本文基于時間依賴的亞臨界裂紋擴展力學模型,引入裂紋成核損傷突變演化路徑函數(shù),并結(jié)合裂紋擴展、應變及生發(fā)射事件相關的損傷模型,提出了一種可以解釋脆性巖石裂紋成核損傷突變,對長期承受恒定應力狀態(tài)(軸向應力σ1,圍壓σ3=σ2)下的蠕變應變演化曲線影響的宏細觀力學模型(圖2).

其中圖2(b) 給出了脆性巖石壓縮蠕變過程中,內(nèi)部可能形成的5 個局部細觀裂紋成核現(xiàn)象,數(shù)字1 ～5 描述了裂紋成核順序,每一個局部細觀裂紋成核都會導致巖石的損傷出現(xiàn)一個突變情況,每一個損傷突變導致了巖石蠕變應變的突然加速增大,文獻[16] 詳細闡述了裂紋成核損傷突變影響下的應力?應變本構(gòu)機理.圖2(c) 給出了一種巖石內(nèi)部裂紋成核可能導致的損傷突變演化路徑,對于真實巖石內(nèi)部,由于細觀裂紋隨機分布特性,巖石內(nèi)部每一個裂紋成核導致的損傷突變大小各不相同,而且相鄰兩次裂紋成核損傷突變之間的時間間隔也是隨機的,因此,圖2(c) 通過?DCN1,?DCN2,?DCN3··· 來描述裂紋成核損傷突變值,?t1描述第一次裂紋成核的時間,?t2,?t3,?t4··· 來描述相鄰裂紋成核的間隔時間.

圖2 (a)巖石裂紋擴展模型;(b)局部裂紋成核模型;(c)裂紋成核損傷隨時間演化曲線Fig.2 (a)Rock model of crack growth;(b)nucleation model of localized cracks;(c)variation of damage from crack nucleation over time

圖1(a)中,壓縮應力作用下,脆性巖石內(nèi)部隨著初始裂紋(其尺寸為a,角度為φ) 沿著裂紋面滑動,導致其兩端新生成翼型裂紋不斷擴展(其長度為l),此翼型裂紋為張開型裂紋,其應力強度因子為[17-18]

其中,β 為常數(shù),μ 為初始裂紋摩擦系數(shù).模型中Fw為作用于初始裂紋面的楔力,其作用可以平衡巖石翼型擴展裂紋之間內(nèi)應力為翼型裂紋擴展極限長度,其物理意義是兩個相鄰翼型裂紋尖端相互貫通時的長度(圖2(a)),llim=[19],式中參數(shù)α=cos φ,巖石初始損傷Do=4πNV(αa)3/3,NV為巖石單位體積的微裂紋數(shù)量.模型計算中定義壓應力為負,結(jié)果分析圖中壓應力取正.圖1(a)模型所有初始裂紋尺寸或角度均假設為恒定值,可理解為真實巖石內(nèi)部隨機分布裂紋尺寸或角度的平均值,即通過一種平均化方法定性研究細觀裂紋影響下的巖石力學性能.眾多學者通過試驗定性地研究了預制微小裂紋尺寸及角度對巖石力學特性影響[20-25],但是微小預制裂紋尺度要遠大于與本文模型中提出細觀裂紋尺度,無法解釋細觀裂紋影響下的完整自然巖石力學特性.然而細觀尺度裂紋對巖石力學特性影響很難試驗研究,隨機裂紋尺寸及角度影響的巖石理論方程又很難建立,所以本文通過平均化假設的方法近似地研究細觀裂紋對巖石力學特性影響.

基于式(1)翼型裂紋尖端應力強度因子,脆性巖石的應力與裂紋擴展關系可表示為[17,19]

式(7) 描述了裂紋擴展硬化與裂紋擴展軟化曲線結(jié)果,其對應著巖石應變硬化與應變軟化的應力?應變曲線,存在一個裂紋啟裂應力值σ1ci(即l=0 時對應的軸向應力),和一個峰值強度值σ1peak,這兩個值為巖石蠕變應變演化結(jié)果計算的恒定軸向應力σ1值選取提供了依據(jù)(即σ1ci<σ1<σ1peak).

張開型亞臨界裂紋擴展速率與強度因子的關系式可表示為dl/dt=v(KI/KIC)n[26],引入式(1)應力強度因子表達式,圖1(a) 中脆性巖石內(nèi)部亞臨界裂紋擴展速率可表示為[17]

式中,v為特征裂紋速度,n為應力侵蝕指數(shù).根據(jù)式(8)解微分方程,可以求解出恒定應力狀態(tài)下裂紋長度隨時間演化結(jié)果l(t),其中微分方程中給定應力狀態(tài)的初始裂紋長度lo可通過式(7)計算得到.

軸向應變相關的巖石損傷可定義為[27]

式中,D為損傷因子,m和εo為材料常數(shù),ε1為軸向應變,γ 為泊松比,E為彈性模量,|·|為取絕對值.

裂紋相關的損傷可定義為[18]

聯(lián)立式(9)與式(10),軸向應變與裂紋長度關系可表示為

將翼型裂紋極限長度llim代入式(11),可以計算得到一個極限軸向應變ε1f,該值為巖石失效應變提供一定參考.

將式(8)計算得到的恒定應力作用下,亞臨界裂紋長度隨時間演化結(jié)果l(t)代入式(11),也就是用l(t)替換l,可以得到巖石的恒定荷載蠕變應變隨時間演化表達式

式(12) 很好地描述了巖石的蠕變應變演化曲線,然而,仍然無法描述脆性巖石內(nèi)部裂紋成核損傷突變對蠕變應變時間演化曲線影響.本文為了能夠通過理論模型描述裂紋成核損傷突變對蠕變演化曲線影響,引入了圖2(c)中裂紋成核損傷蠕變演化路徑,將該損傷函數(shù)引入式(12)則可以很好解釋裂紋成核損傷突變影響.然而,如何將裂紋成核損傷突變演化路徑函數(shù)引入式(12) 是一個關鍵問題.具體解決方法如下:通過式(12)可以看到,巖石蠕變應變時間演化結(jié)果與巖石初始損傷Do相關,針對一個巖石試樣來說,其初始損傷是恒定不變的,然而,為了描述裂紋成核損傷突變影響,將圖2(c) 提出的裂紋成核損傷蠕變時間演化路徑函數(shù)DCN(t)替換式(12)中的初始損傷Do,即每一次裂紋成核引起的巖石損傷突變值可以通過初始損傷突變近似地描述.則脆性巖石裂紋成核損傷突變影響下的蠕變時間演化表達式為

此外值得注意的是式(13)中裂紋長度時間演化參數(shù)l(t) 中同時也隱含著裂紋成核損傷突變演化函數(shù)DCN(t).具體用DCN(t)替換Do的原因:已知初始損傷Do=4πNV(αa)3/3,假設初始裂紋尺寸a及角度φ 保持不變,則在上述理論計算中每一次裂紋成核引起的損傷突變?DCN,導致了裂紋數(shù)量的NV的突然增大,這也近似地實現(xiàn)了裂紋成核損傷突變結(jié)果.實際巖石內(nèi)部裂紋成核損傷突變,不僅僅是裂紋數(shù)量的增多,裂紋尺寸、方向及位置也可能在不斷變化,由于本文模型局限性,無法真實模擬這些裂紋變化特征,所以只能通過改變裂紋數(shù)量近似等效的方法實現(xiàn)裂紋成核現(xiàn)象,為脆性巖石細觀裂紋成核與宏觀力學特性關系的理解提供一定參考.

此外,式(10) 給出了圖2(a) 脆性巖石細觀裂紋模型的損傷表達式,拆分變換式(10),巖石損傷方程也可以表示為

在式(14)中,巖石裂紋損傷演化與初始損傷Do密切相關,因此,通過控制初始損傷方法近似地實現(xiàn)裂紋成核損傷突變的方法是可行的.具體理論結(jié)果是否合理詳見下文結(jié)果分析.

此外,基于聲發(fā)射監(jiān)測數(shù)據(jù),巖石損傷為[28]

式中,Cd為巖石某一階段的累計聲發(fā)射數(shù),C0為巖石壓縮破壞全過程中累計聲發(fā)射數(shù),σr為殘余強度,σP為峰值強度.式(15) 可以很好地呈現(xiàn)蠕變試驗過程中,巖石內(nèi)部局部微裂紋成核損傷突變現(xiàn)象(見圖1).因此,可以通過巖石內(nèi)部裂紋擴展誘發(fā)的聲發(fā)射信號數(shù)據(jù)更加準確地判斷巖石內(nèi)部裂紋成核損傷突變演化路徑函數(shù)DCN(t),進而為通過式(13)更加真實地計算模擬巖石內(nèi)部裂紋成核規(guī)律影響下的巖石蠕變演化曲線結(jié)果提供重要依據(jù).

2 模型結(jié)果分析

2.1 模型合理性驗證

依據(jù)鹽巖力學性能[12],模型參數(shù)選取如下:初始損傷Do=0.06,初始裂紋摩擦系數(shù)、尺寸及角度分別為μ=0.6,a=2.9 mm,φ=45?.材料常數(shù)m=2,εo=0.1,β=0.32,斷裂韌度KIC=0.32 MPa·m1/2,彈性模量E=20 GPa,泊松比γ=0.3,特征裂紋速率v=0.16 m/s,應力侵蝕指數(shù)n=57,?t1=1440 h,?t2=4296 h,?DCN1=0.005,?DCN2=0.001 7.具體參數(shù)選取方法可參考文獻[17,19].

圖3 中給出了考慮裂紋成核損傷突變影響下的理論蠕變應變演化曲線,并通過與試驗結(jié)果[12](即圖1 試驗曲線)對比驗證模型合理性.巖石蠕變過程中分別在t1=1440 h,t2=4296 h 時,經(jīng)歷兩次由加速蠕變到減速蠕變的轉(zhuǎn)化,其分別對應兩次局部裂紋成核損傷突變,在模型中該兩次成核損傷分別賦值?DCN1=0.005,?DCN1=0.001 7,此外,巖石在最后產(chǎn)生加速蠕變變形,導致了巖石的最終破壞.

圖3 裂紋成核損傷突變下的理論與試驗[12] 蠕變演化曲線Fig.3 Theoretical and experimental[12] creep curves influenced by damage catastrophe from crack nucleation

此外,為進一步驗證模型合理性,圖4(a)給出了錦屏大理巖試樣蠕變演化過程中,歸一化裂紋損傷演化理論結(jié)果與數(shù)值試驗結(jié)果[29]對比分析曲線.由于本文模型局限性,理論與試驗結(jié)果數(shù)值上存在一定差異,因此采用了歸一化結(jié)果進行對比分析.模型參數(shù)如下:?DCN1=?DCN2=···?DCN9=0.000 5,?t1=?t2=···?t9=2.0×105s.數(shù)值試驗結(jié)果存在多次裂紋成核損傷突增現(xiàn)象,只是數(shù)值試驗更好地呈現(xiàn)了裂紋成核損傷時間、大小及數(shù)量的隨機特性,理論結(jié)果假設每一次裂紋成核損傷大小及時間間隔一定,不能很好地呈現(xiàn)隨機性,但是理論與數(shù)值試驗變化趨勢具有一定相似性.

圖4(b) 給出了地殼孕震階段的地殼應變演化曲線[30]與本文理論結(jié)果的歸一化應變演化對比曲線.由于本文采用實驗室尺度巖石模型,數(shù)值上難以與地殼尺度巖石達到一致,因此采用了歸一化對比分析.模型參數(shù)如下:KIC=1.61 MPa·m1/2,n=57,v=0.16 m/s,m=1,εo=1/68,β=0.276,φ=32?,a=3.1 mm,μ=0.51,Do=0.048,σ1=100 MPa,σ3=4 MPa,γ=0.2,?t1=4.8×105s,?t=3.6×105s,?t3=6.2×105s,?t4=?t5=?t8=1.4×106s,?t6=?t7=0.6×105s,?DCN1=?DCN3=?DCN5=?DCN6=?DCN7=5.0×10?4,?DCN2=1.0×10?3,?DCN4=1.0×10?3,?DCN8=2.5×10?3.地殼孕震階段,變形監(jiān)測區(qū)域的地應力可以近似理解為處于一個恒定應力狀態(tài),此時地殼變形可以近似理解為一個蠕變演化過程,由于地殼局部裂隙貫通成核損傷突變,導致了應變的突變,進而導致局部地震.成核損傷大小決定了應變增量大小,進而也影響地震震級.雖然數(shù)值上存在差異,但是演化趨勢相似,一定程度上解釋了裂紋成核損傷突變機理,為孕震階段應變演化評價提供了一定參考.

圖4 巖石蠕變裂紋成核(a)歸一化損傷,(b)應變演化理論與試驗[29-30] 對比曲線Fig.4 Comparison between theoretical and experimental[29-30] results of the evolution curves in(a)the normalized damage and(b)strain influenced by crack nucleation during creep

以上理論結(jié)果與試驗結(jié)果變化趨勢相似,然而在數(shù)值上存在一定差異,原因可能是:(1) 本文理論模型分別假設巖石內(nèi)部初始裂紋尺寸及角度均為相同的值,通過一種平均化方法近似地理論分析了隨機分布裂紋巖石力學性能;(2)模型中其他參數(shù)選取時存在一定誤差.

此外,不可能為了模擬蠕變結(jié)果,每一個理論模擬都要配合開展聲發(fā)射試驗,來確定裂紋成核突變發(fā)生時間、裂紋成核損傷突變大小,及裂紋成核數(shù)量,再反過來應用于理論模型中,實現(xiàn)模擬結(jié)果.對于相同類型的巖石,可通過監(jiān)測不同受力狀態(tài)的大量巖石試樣蠕變聲發(fā)射信號數(shù)據(jù)規(guī)律,近似地匯總統(tǒng)計出此類巖石不同應力狀態(tài)下的裂紋成核時間、成核損傷大小及成核數(shù)量變化規(guī)律,進而應用于本文模型中,為預測真實巖石蠕變裂紋成核現(xiàn)象提供一定支持.本文理論模型還很難模擬真實巖石隨機分布裂紋成核時間、成核損傷大小及成核數(shù)量的隨機特性,裂紋成核導致的損傷突變大小、裂紋成核發(fā)生時間以及裂紋成核數(shù)量影響下的脆性巖石蠕變演化,定性地研究了脆性巖石細觀裂紋蠕變成核機理.下面將對這些因素影響下的蠕變應變演化規(guī)律進行詳細討論分析.

2.2 裂紋成核損傷變量對蠕變演化影響

為了便于模型參數(shù)對比分析,下面討論中模型參數(shù)選取如下:Do=0.048,μ=0.51,a=3.1 mm,φ=45?,m=1,εo=1/68,β=0.32,KIC=1.61 MPa·m1/2,E=30 GPa,v=0.16 m/s,n=57,γ=0.2,σ1=90 MPa,σ3=2 MPa.

圖5 給出了當巖石整個蠕變演化過程中,只有一個裂紋成核損傷突變?DCN1情況下,裂紋成核損傷大小對蠕變過程中翼型裂紋擴展長度及裂紋速率隨時間演化曲線的影響結(jié)果.對于不同的裂紋成核損傷下,當發(fā)生裂紋成核損傷前,裂紋長度或裂紋速率隨時間演化曲線是一致的;當裂紋成核損傷發(fā)生時,隨著裂紋成核損傷突變值增大,在給定的時間點,裂紋長度或裂紋速率均增大,裂紋擴展至巖石最終失效長度時的時間越短.

圖5 中單個裂紋成核損傷大小對蠕變(a)裂紋長度及(b)裂紋速率演化曲線影響Fig.5 Effect of the size of individual crack nucleation damage on creep evolutions of(a)crack length and(b)crack velocity

圖6 給出了巖石蠕變演化過程中,裂紋成核損傷大小對宏觀軸向應變及應變率時間演化曲線的影響結(jié)果.該變化趨勢與上述裂紋長度及裂紋速率演化趨勢相似.在裂紋成核損傷突變發(fā)生后,在給定的時間,隨著裂紋成核損傷的增大,巖石的軸向應變及應變率均增大.隨著裂紋成核損傷突變的增大,巖石最終蠕變失效時間越短,圖7 中也清楚的給出了巖石最終蠕變失效時間隨裂紋成核損傷的變化趨勢.

圖6 單個裂紋成核損傷大小對蠕變(a)軸向應變及(b)應變率演化曲線影響Fig.6 Effect of the size of individual crack nucleation damage on evolutions of(a)axial strain and(b)strain rate during creep

圖7 單個裂紋成核損傷大小對巖石蠕變失效時間影響Fig.7 Effect of the size of individual crack nucleation damage on creep failure time of rocks

圖8 中給出了脆性巖石蠕變演化過程中,當只發(fā)生一次裂紋成核損傷突變時,裂紋成核損傷突變所發(fā)生的時間節(jié)點?t1對巖石蠕變裂紋擴展長度及裂紋速率的影響.在裂紋成核損傷突變發(fā)生后,隨著裂紋成核損傷時間點的增大,裂紋擴展至最終失效長度的時間越長,即裂紋成核越晚,巖石直至最終失效的時間越長,圖10 中也進一步給出了證明,巖石蠕變失效時間隨裂紋成核時間成線性增大.

圖8 單個裂紋成核損傷發(fā)生的時間對蠕變(a)裂紋長度及(b)裂紋速率演化曲線影響Fig.8 Effect of the happened time at individual crack nucleation damage on the evolutions of(a)crack length and(b)crack velocity during creep

圖 9給出了裂紋成核損傷時間節(jié)點對巖石蠕變軸向應變與應變率的影響結(jié)果.軸向應變及應變率演化趨勢與裂紋長度變化趨勢相似.

圖9 單個裂紋成核損傷發(fā)生時間對蠕變(a)軸向應變及(b)應變率演化曲線影響Fig.9 Effect of the happened time at individual crack nucleation damage on the evolutions of(a)axial strain and(b)strain rate during creep

圖10 單個裂紋成核損傷發(fā)生時間對巖石蠕變失效時間影響Fig.10 Effect of the happened time of individual crack nucleation damage on creep failure time of rocks

以上討論中,均是給出了一次裂紋成核損傷突變影響下,裂紋成核損傷大小及發(fā)生時間對巖石蠕變演化的影響.然而,對于真實巖石中,在整個長期蠕變演化過程中,由于巖石內(nèi)部裂紋隨機分布特性,導致巖石的非均勻力學性能,因此巖石內(nèi)部可能會隨機出現(xiàn)多次裂紋成核損傷突變現(xiàn)象.因此,下面將對多次裂紋成核損傷發(fā)生情況下的巖石蠕變裂紋及應變演化結(jié)果進行討論.

圖11(a)分別給出了發(fā)生2,3,4 次裂紋成核損傷突變影響下的裂紋長度隨時間演化曲線,每一次成核損傷值?DCN1,?DCN2,?DCN3,?DCN4均為0.005,第一次裂紋成核?t為0.6×108s,其他裂紋成核損傷之間時間間隔?t2,?t3,?t4均為0.2×108s.圖11(b)、圖12(a) 及圖12(b) 分別給出了圖11(a) 相對應的裂紋速率、軸向應變及應變率演化曲線.當每一次裂紋成核損傷突變值及相鄰裂紋成核損傷時間差一定時,隨著裂紋成核數(shù)量的增大,巖石蠕變失效時間減小,圖13 也總結(jié)歸納出了巖石蠕變失效時間隨裂紋成核數(shù)量變化規(guī)律.

圖11 裂紋成核數(shù)量對蠕變(a)裂紋長度及(b)裂紋速率演化曲線影響Fig.11 Effect of the number of crack nucleation on the evolutions of(a)crack length and(b)crack velocity during creep

圖12 裂紋成核數(shù)量對蠕變(a)軸向應變(b)應變率演化曲線影響Fig.12 Effect of the number of crack nucleation on the evolutions of(a)axial strain and(b)strain rate during creep

圖13 裂紋成核數(shù)量對巖石蠕變失效時間影響Fig.13 Effect of the number of crack nucleation on creep failure time of rocks

上述關于多次裂紋成核損傷對蠕變演化影響討論中,裂紋成核損傷大小及裂紋成核發(fā)生時間點均是被假設為定值,圖14 給出了巖石蠕變演化過程中,多次裂紋成核現(xiàn)象出現(xiàn)時,相鄰裂紋成核損傷大小及相鄰裂紋成核時間間隔不同的情況下,巖石內(nèi)部裂紋長度、裂紋速率、軸向應變及應變隨時間演化趨勢.此時,巖石蠕變過程中,裂紋長度、裂紋速率、軸向應變及應變率突變性更加接近呈現(xiàn)一種隨機特性.對于深部地下含有隨機分布裂紋的自然巖石,若要更加精確預測巖石非均勻性影響下的蠕變力學性能,可以通過圖1 與式(15) 中提到的聲發(fā)射監(jiān)測數(shù)據(jù),分析不同巖石裂紋擴展導致的聲發(fā)射率及累計聲發(fā)射變化值,預測裂紋成核損傷大小,并結(jié)合本文提出的理論模型方法.該模型的提出為深部地下工程圍巖長期穩(wěn)定性及短期施工階段微震及巖爆災害評價預測提供了重要理論支持.

圖14 多變量裂紋成核損傷對蠕變(a)裂紋長度及裂紋速率(b)軸向應變及應變率演化曲線影響Fig.14 Effect of the multiple variables of crack nucleation damage on the evolutions of(a)crack length and crack velocity,(b)axial strain and strain rate during creep

3 結(jié)論

基于固體斷裂損傷力學理論,推出了脆性巖石裂紋成核損傷突變影響下的蠕變演化機理模型.并分析討論了裂紋成核損傷大小、時間及數(shù)量對脆性巖石蠕變裂紋及應變演化特性影響.具體結(jié)論如下:

(1)巖石蠕變過程中,在同一時間發(fā)生一次裂紋成核損傷突變時,隨著裂紋成核損傷大小的增大,巖石裂紋長度、裂紋速率、軸向應變及應變率均增大,巖石最終蠕變失效時間減小.

(2)巖石蠕變過程中,發(fā)生一次裂紋成核損傷突變時,裂紋成核損傷突變發(fā)生的時間越早,巖石最終發(fā)生失效的時間越短,而且?guī)r石蠕變失效時間與裂紋成核時間成線性遞增關系.

(3)在巖石蠕變過程中,發(fā)生多次裂紋成核現(xiàn)象時,當巖石裂紋成核損傷大小及相鄰裂紋成核時間間隔相同的情況下,隨著裂紋成核數(shù)量的增大,巖石蠕變失效時間減小;當巖石裂紋成核損傷大小及裂紋成核時間間隔不同且隨機時,巖石蠕變變形演化更加接近含隨機分布裂紋的真實巖石非均勻性導致的蠕變演化規(guī)律,其也為深部地下工程微震及巖爆事件預測評價提供一定參考.