李艷龍陳云娟高濤敬藝馬永安

(1.山東建筑大學(xué)土木工程學(xué)院，山東 濟(jì)南 250101；2.山東建筑大學(xué)建筑結(jié)構(gòu)加固改造與地下空間工程教育部重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250101)

0 引言

自然界中，由于長期的各種地質(zhì)作用，巖體內(nèi)部將會(huì)出現(xiàn)各種各樣的節(jié)理、裂隙、孔隙、斷層等弱結(jié)構(gòu)面，其相互交叉切割，形成了更加復(fù)雜的組合弱結(jié)構(gòu)面。在工程擾動(dòng)作用下，這些面將嚴(yán)重影響巖體的力學(xué)行為、破壞變形特性、損傷變化規(guī)律以及工程進(jìn)度等，研究含弱結(jié)構(gòu)面巖體的破壞規(guī)律對(duì)于保證工程順利施工具有重要的理論價(jià)值。在一些地質(zhì)工程中，裂隙的分布形態(tài)對(duì)巖體的破壞行為和力學(xué)行為具有很大影響[1-3]。因此，研究不同角度節(jié)理巖體的圍壓卸壓損傷規(guī)律，揭示不同角度裂隙擴(kuò)展貫通機(jī)制，對(duì)掌握隧道開挖工程破壞和失穩(wěn)具有重大意義。

學(xué)者們對(duì)裂隙巖體的力學(xué)性質(zhì)及損傷演化規(guī)律開展了許多研究。孫旭曙等[4]通過不同角度單一預(yù)制貫通節(jié)理巖體試件并開展常規(guī)三軸壓縮試驗(yàn)，系統(tǒng)分析了巖體的彈性模量、黏聚力、抗壓強(qiáng)度等與節(jié)理傾角之間的內(nèi)在關(guān)系；李勇等[5]基于室內(nèi)試驗(yàn)與PFC2D模擬軟件，研究了在不同裂隙傾角和單軸壓縮作用下的貫通機(jī)理，并從應(yīng)力場和位移場分析了平行雙裂隙的裂紋擴(kuò)展演化規(guī)律；李露露等[6]研究了不同角度裂隙條件下三叉裂隙類巖石試樣在單軸壓縮條件下的破壞規(guī)律，從宏觀和細(xì)觀的角度分析了交叉裂隙角度的破裂裂紋演化擴(kuò)展規(guī)律；王艷磊等[7]通過室內(nèi)試驗(yàn)系統(tǒng)分析了單軸壓縮下不同裂隙數(shù)量和不同裂隙開度條件下對(duì)低強(qiáng)度巖石應(yīng)力-應(yīng)變曲線和破壞模式的影響；汪杰等[8]提出了節(jié)理傾角對(duì)巖體初始損傷的影響，運(yùn)用損傷力學(xué)理論，建立了節(jié)理巖體損傷演化及損傷本構(gòu)模型；張波等[9-11]通過單軸壓縮試驗(yàn)對(duì)含有充填材料的節(jié)理巖體以及自制交叉裂隙相似材料試件進(jìn)行了科學(xué)探索，研究了裂隙充填與否對(duì)節(jié)理巖體力學(xué)行為的影響，并利用有限元軟件深入研究了試件的斷裂損傷機(jī)理。在裂隙巖體破裂擴(kuò)展過程中，從能量方面對(duì)巖體進(jìn)行監(jiān)測(cè)同樣是一個(gè)重要依據(jù)[12-13]。王桂林等[14]基于室內(nèi)單軸試驗(yàn)和巖石能量理論，研究不同傾角節(jié)理砂巖能量演化機(jī)制，并將節(jié)理砂巖能量演化過程劃分為5段；PENG等[15]基于常規(guī)三軸壓縮試驗(yàn)，研究了能量轉(zhuǎn)化與試樣破壞之間的關(guān)系，并建立了其在峰前剛度退化的損傷演化模型。蒲超等[16]對(duì)千枚巖進(jìn)行不同圍壓下的三軸壓縮試驗(yàn)，研究圍壓對(duì)千枚巖變形破壞特征和能量演化特征的影響。

上述針對(duì)節(jié)理巖體試驗(yàn)研究主要采用單軸、雙軸或常規(guī)三軸試驗(yàn)。文章在此基礎(chǔ)上采用數(shù)值模擬軟件FLAC3D，基于單節(jié)理砂巖試樣模擬含節(jié)理隧道開挖，再結(jié)合試樣塑性區(qū)的發(fā)展?fàn)顩r全過程記錄并分析了裂隙破壞情況，并對(duì)比分析各工況作用下不同角度節(jié)理試樣的裂隙擴(kuò)展行為，分析其裂隙發(fā)展規(guī)律和特征，再根據(jù)能量耗散原理，進(jìn)一步研究了單節(jié)理砂巖力學(xué)性質(zhì)和破壞模式，為實(shí)際隧道工程建設(shè)提供了理論依據(jù)和技術(shù)參考。

1 數(shù)值模型建立

為研究不同角度節(jié)理對(duì)巖體裂隙擴(kuò)展規(guī)律影響，采用摩爾—庫倫模型，模擬含節(jié)理砂巖試樣在真三軸三向六面受力狀態(tài)下的單面卸荷試驗(yàn)，研究不同角度節(jié)理對(duì)巖體試樣的破裂規(guī)律與特征。以30°節(jié)理砂巖試件為例，模型示意圖如圖1所示，其中模型尺寸為70 mm×45 mm×140 mm(長×寬×高)的長方體試樣，節(jié)理傾角(節(jié)理與水平面的夾角)α分別為15°、30°、45°、60°、75°共5種工況。節(jié)理長度為15 mm，其位置設(shè)置在模型正中間。預(yù)制節(jié)理位置及距離示意圖如圖2所示[17]。砂巖的模擬參數(shù)見表1。

表1 砂巖FLAC3D模擬參數(shù)表

圖1 30°節(jié)理砂巖試樣模型圖

圖2 預(yù)制節(jié)理位置及距離示意圖/mm

2 節(jié)理砂巖裂隙擴(kuò)展分析

為了探索單節(jié)理砂巖試樣卸荷過程中試樣變形破壞的本質(zhì)特征，采用FLAC3D軟件對(duì)試樣進(jìn)行相同軸壓作用下的簡化三軸單面卸荷模擬試驗(yàn)。數(shù)值模擬時(shí)，正值應(yīng)力區(qū)域?yàn)槭芾瓍^(qū)，負(fù)值應(yīng)力區(qū)域?yàn)槭軌簠^(qū)。選取試驗(yàn)?zāi)M方案為:挖掘含有裂隙的隧道，實(shí)際上是周圍巖體一側(cè)圍壓消失，且軸壓不變的過程，故模擬試驗(yàn)分別施加低圍壓5、8和11 MPa，同時(shí)施加軸壓45 MPa，運(yùn)行達(dá)到平衡；卸去不含節(jié)理一側(cè)的圍壓σ3，其它側(cè)圍壓保持不變，運(yùn)行達(dá)到平衡。5種不同節(jié)理傾角的試樣試驗(yàn)均采用此方案。

2.1 相同圍壓、軸壓卸荷前裂隙擴(kuò)展分析

在相同圍壓、軸壓作用下，對(duì)試樣進(jìn)行簡化三軸數(shù)值模擬試驗(yàn)。5種不同角度單節(jié)理試樣在相同圍壓、軸壓作用下的模擬破壞(圍壓為5 MPa)如圖3所示。

圖3 不同角度單節(jié)理砂巖試樣未卸荷(圍壓5 MPa)破壞圖

5種工況試樣在相同圍壓、軸壓作用后，不同角度單節(jié)理試樣的裂隙發(fā)育大致相同，主要在預(yù)制節(jié)理面尖部產(chǎn)生剪切翼裂紋，為張拉破壞。

不同圍壓作用下單節(jié)理試樣塑性區(qū)體積與節(jié)理角度關(guān)系圖如圖4所示。在相同圍壓、軸壓作用下，節(jié)理傾角由15°增加到75°，對(duì)試樣的破壞影響規(guī)律基本一致，塑性區(qū)體積基本上呈先上升后下降的趨勢(shì)。節(jié)理傾角為30°～60°的試樣破壞最大，15°試樣和75°試樣相對(duì)于其他角度試樣破壞體積較小。

圖4 不同角度單節(jié)理試樣塑性區(qū)體積與節(jié)理角度關(guān)系圖

2.2 相同圍壓、軸壓卸荷后裂隙擴(kuò)展分析

在相同圍壓、軸壓作用下，對(duì)試樣進(jìn)行單面卸荷數(shù)值模擬。5種工況試樣主要生成翼裂紋與次生裂紋，最初破壞產(chǎn)生于預(yù)制節(jié)理面尖部的翼裂紋，并沿節(jié)理的方向進(jìn)行擴(kuò)展，主要為張拉破壞，然后出現(xiàn)次生裂紋進(jìn)而導(dǎo)致試樣發(fā)生剪切破壞。5種不同角度單節(jié)理試樣單面卸荷后模擬塑性破壞圖和最大主應(yīng)力圖(圍壓5 MPa)分別如圖5、6所示。

圖5 不同角度單節(jié)理砂巖試樣卸荷(圍壓5 MPa)塑性破壞圖

15°和60°節(jié)理試樣，在卸荷平衡過程中，節(jié)理面受壓導(dǎo)致沿預(yù)制裂隙面發(fā)生多處破壞，試樣破壞主要為節(jié)理尖端產(chǎn)生的翼裂紋與位于預(yù)制節(jié)理左上部位的次生裂隙，以及卸壓面部位的張拉破壞。與 60°試樣相比，15°試樣次生裂隙發(fā)育不明顯。

圖6 不同角度單節(jié)理砂巖試樣卸荷(圍壓5 MPa)最大應(yīng)力圖

30°和45°節(jié)理試樣，整體破壞嚴(yán)重，整個(gè)試樣貫通破壞，在預(yù)制節(jié)理面中與15°、60°和75°試樣相比，其翼裂紋和次生裂紋更加明顯，發(fā)育程度更加完整，這也說明了30°和45°試樣的抗壓強(qiáng)度較低一些。30°和45°模擬試樣相比其他角度試樣存在大量的抗拉破壞單元，主要集中在節(jié)理尖端部位，但整體上看還是以剪切破壞為主；30°和45°試樣裂紋的擴(kuò)展方向與最大主應(yīng)力的方向平行，這說明了張拉應(yīng)力的增長對(duì)裂紋的擴(kuò)展方向和速度有一定的影響，而且兩者單節(jié)理砂巖試樣平衡時(shí)的破壞形式都大致呈y形。

75°試樣，整體破壞最小，僅在預(yù)制節(jié)理尖部與卸荷面存在破壞單元，其破壞類型以剪切破壞為主，這說明75°試樣比其他角度節(jié)理試樣的抗壓能力強(qiáng)，不易被破壞。其最大應(yīng)力相比其他節(jié)理角度試樣較小，張拉應(yīng)力主要分布在節(jié)理左側(cè)，卸荷面右側(cè)。

2.3 不同圍壓裂隙擴(kuò)展分析

在不同圍壓、相同軸壓作用下，對(duì)試樣僅一面進(jìn)行卸壓模擬。采用控制變量的方式，對(duì)相同節(jié)理角度的單節(jié)理試樣，在恒定軸壓σ1作用下，設(shè)定不同圍壓σ3進(jìn)行模擬。30°單節(jié)理砂巖試樣分別在圍壓5、8、11 MPa作用下單面卸荷模擬塑性圖和最大應(yīng)力分別如圖7、8所示。

圖7 不同圍壓作用下30°單節(jié)理砂巖試樣模擬塑性圖

在相同軸壓作用下，加載初始值分別為5、8和11 MPa的作用過程中，試樣將會(huì)產(chǎn)生少量破壞或者不破壞；在進(jìn)行單面卸荷(其他條件均不變)過程中，其他三面圍壓限制了試樣在軸壓作用下側(cè)面方向的擴(kuò)張，使其向臨空面發(fā)展；在恒定軸壓作用下，卸荷面試樣從受壓轉(zhuǎn)換成受拉，再逐漸到受拉破壞，主要是在卸荷面和節(jié)理面附近產(chǎn)生貫通試樣的受拉剪切破壞。在卸荷破壞前，圍壓為11 MPa時(shí)的試樣比圍壓為5 MPa時(shí)的試樣承受能力較強(qiáng)，破壞范圍更小，積累能量更多，消耗能量更?。幌啾葒鷫簽?1 MPa時(shí)的試樣，圍壓為5 MPa時(shí)的試樣在承受相同的軸壓作用下，其抗壓和抗拉能力較小，塑性區(qū)破壞的范圍更大一些。在卸荷破壞后，圍壓為11 MPa時(shí)的試樣較其他兩種情況的試樣塑性區(qū)破壞范圍較小，主要是因?yàn)槠渌葱秶鷫旱?個(gè)加載面的圍壓相對(duì)較大，能夠增大單節(jié)理砂巖試樣的抗壓、抗拉能力以及整體穩(wěn)定能力，使得單節(jié)理試樣破壞范圍更小，能夠盡快進(jìn)入平衡狀態(tài)。

圖8 不同圍壓作用下30°單節(jié)理砂巖試樣模擬最大應(yīng)力圖

3 能量耗散規(guī)律研究

3.1 能量計(jì)算原理

巖石在加卸載破壞過程中，能量轉(zhuǎn)化的變形特征可以分為可逆變形和不可逆變形，前者可以產(chǎn)生彈性應(yīng)變能，而后者則是以塑性變形、損傷、摩擦等方式而產(chǎn)生耗散能。巖石的能量轉(zhuǎn)化與變形破壞密不可分，應(yīng)用數(shù)值模擬軟件FLAC3D對(duì)巖石破壞全過程分析時(shí)，輸入總能量主要轉(zhuǎn)化為彈性應(yīng)變能和耗散能對(duì)模擬試樣產(chǎn)生影響，且從熱力學(xué)第一定律角度出發(fā)，忽略模擬系統(tǒng)與外界發(fā)生的熱交換現(xiàn)象，那么單節(jié)理巖體變形過程中輸入總能量將全部用在儲(chǔ)存于巖體內(nèi)部的彈性應(yīng)變能以及巖體內(nèi)部耗散的能量，其計(jì)算公式由式(1)和(2)表示為

式中W為外界輸入能，J；We為儲(chǔ)存于巖石內(nèi)可釋放的彈性應(yīng)變能，J；Wd為耗散能，主要消耗于巖石內(nèi)部損傷和塑性變形，彈性應(yīng)變能分解為偏差分量和體積分量，J；V為子單元體積，m3；G為剪切模量，Pa；K為體積模量，Pa；為偏應(yīng)力，Pa；為平均應(yīng)力，Pa。

3.2 單節(jié)理砂巖試樣單面卸荷能量分析

根據(jù)上述能量計(jì)算原理，分析單節(jié)理砂巖能量演化過程，揭示不同傾角的能量積累與耗散特征。通過FLAC3D內(nèi)置的fish語言，可跟蹤分析單節(jié)理砂巖試樣單面卸荷過程中內(nèi)部能量積累與消耗情況，實(shí)時(shí)監(jiān)測(cè)內(nèi)部耗散能及應(yīng)變能，即可分析試樣在圍壓卸壓過程中能量的變化規(guī)律。峰值前后能量對(duì)比見表2，表中差值耗散能等于峰后耗散能減去峰前耗散能，差值塑性區(qū)體積等于峰值后塑性區(qū)體積減去峰值前塑性區(qū)體積。峰值前后耗能及塑性區(qū)體積差對(duì)比線如圖9所示。

表2 試樣在圍壓卸壓過程中峰值前后能量對(duì)比表

圖9 不同角度節(jié)理峰值前后耗散能差及塑性區(qū)體積差對(duì)比曲線圖

根據(jù)表2和圖9的分析可知，在相同圍壓、軸壓作用下，節(jié)理傾角在15°～90°范圍內(nèi)，峰值前后差值塑性區(qū)總體積基本上呈現(xiàn)出先增加后減少的趨勢(shì)，曲線近似成倒U形，且相同低圍壓、軸壓作用下不同角度節(jié)理試樣差值耗散能與差值塑性區(qū)體積演化規(guī)律基本相似，差值耗散能與差值塑性區(qū)體積呈正相關(guān)。

不同角度的節(jié)理巖體試樣在相同低圍壓、軸壓條件下，結(jié)合圖5中節(jié)理試樣卸壓破壞塑性圖可知，30°和45°試樣的抗壓強(qiáng)度較低，抵抗變形能力較弱，破壞區(qū)域較為嚴(yán)重，導(dǎo)致卸壓后耗散能增長較大；15°、60°和75°試樣承受能力較強(qiáng)，破壞較少，耗散能增長較小，這也是不同角度節(jié)理試樣卸壓前后耗散能差相差幾十倍的原因。

4 含有多裂隙的隧道開挖數(shù)值模擬

建立含不同角度節(jié)理的隧道數(shù)值模型，研究相同工況作用下，隧道開挖前后不同角度節(jié)理對(duì)隧道開挖時(shí)的影響?？紤]到邊界效應(yīng)的影響，隧道模型邊界設(shè)置為隧道開挖尺寸的3.5倍，這樣對(duì)隧道模型開挖后的分析更具有準(zhǔn)確性。建立的隧道模型尺寸為60 m×64 m×60 m(長×寬×高)，隧道口尺寸為8 m×10 m(寬×高)。

三維數(shù)值計(jì)算模型采用六面體單元進(jìn)行網(wǎng)格劃分，整個(gè)模型共有306 000個(gè)實(shí)體單元，309 240個(gè)網(wǎng)格節(jié)點(diǎn)，所有單元采用摩爾-庫倫模型。模型邊界條件為:前后兩側(cè)邊界施加y方向位移約束，左右兩側(cè)邊界施加x方向位移約束，底部邊界受z方向位移約束，上部地表為自由邊界。模型隧道口兩側(cè)分布有12道節(jié)理，其傾角分別為15°、30°、45°、60°和75°，每道節(jié)理尺寸為4 m×1.5 m(長×寬)。隧道開挖數(shù)值模型的塑性圖和最大主應(yīng)力圖分別如圖10、11所示，其中節(jié)理均勻分布在模擬隧道口兩側(cè)，并且角度一致，完全相同。

圖10 不同節(jié)理傾角塑性圖

圖11 不同節(jié)理傾角最大應(yīng)力圖

根據(jù)圖10、11分析可知，破壞塑性區(qū)主要分布在隧道口兩側(cè)和節(jié)理處，破壞情況基本上是剪切破壞。節(jié)理與隧道之間的巖橋是拉剪破壞，破壞面是從隧道口底部逐漸向上開裂破壞直到與節(jié)理連通。與隧道口較近的第一排節(jié)理，節(jié)理之間的區(qū)域均為剪切破壞，最大主應(yīng)力為正值，故節(jié)理與節(jié)理之間的巖橋同樣是拉剪破壞；與隧道口較遠(yuǎn)的第二排節(jié)理，主要是在節(jié)理尖端部位發(fā)生剪切破壞，破壞較小，對(duì)隧道開挖的影響也較小。

相同低圍壓、軸壓作用下，節(jié)理隧道開挖前后能量對(duì)比見表3，開挖前后差值耗散能與差值塑性區(qū)對(duì)比曲線如圖12所示。開挖前后達(dá)到平衡時(shí)得到的差值塑性區(qū)體積和差值耗散能均隨著角度的增大呈現(xiàn)先增大后減小的規(guī)律，曲線近似成倒U形，與上述單節(jié)理砂巖的差值塑性區(qū)體積和差值耗散能所得到的規(guī)律相似。

圖12 節(jié)理隧道開挖前后差值耗散能與差值塑性區(qū)對(duì)比曲線圖

表3 節(jié)理隧道開挖前后能量表

5 結(jié)論

為研究不同角度節(jié)理對(duì)巖石物理性質(zhì)的影響程度，開展了含單節(jié)理巖體低圍壓卸壓作用下的模擬試驗(yàn)，并通過能量角度對(duì)其進(jìn)行分析與討論，得到以下結(jié)論:

(1)在低圍壓作用下單節(jié)理試樣卸荷前后峰值時(shí)，相比15°和75°節(jié)理試樣，30°～60°之間的節(jié)理試樣發(fā)育程度較高，且破壞形式大致呈y形。

(2)在相同恒定軸壓作用下，在低圍壓范圍內(nèi)進(jìn)行第三主應(yīng)力σ3單面卸荷過程中，隨著圍壓的增大，含單節(jié)理巖體試樣破壞更不易，超過圍壓一定范圍，其破壞更加嚴(yán)重。

(3)相同低圍壓、軸壓作用下不同角度節(jié)理試樣差值耗散能隨著角度的增大呈現(xiàn)先增大后減小的規(guī)律，曲線近似成倒U形，差值耗散能與差值塑性區(qū)體積演化規(guī)律基本相似，并且開挖過程中兩演化規(guī)律與數(shù)值模擬規(guī)律較為吻合。