張 科，劉享華，楊紅宣，范文臣

(1.昆明理工大學(xué)電力工程學(xué)院，云南 昆明 650500；2.昆明理工大學(xué)建筑工程學(xué)院，云南 昆明 650500；3.中南大學(xué)資源與安全工程學(xué)院，湖南 長沙 410083)

由于漫長的地質(zhì)構(gòu)造作用，巖體內(nèi)部除了存在裂隙缺陷，也孕育了孔洞缺陷[1]。這些地質(zhì)缺陷是影響巖體強(qiáng)度和變形破壞特性的主要因素，也是影響工程安全的重要地質(zhì)條件。因此，開展含地質(zhì)缺陷的巖體力學(xué)特性試驗(yàn)研究，對(duì)于確保巖體工程的穩(wěn)定具有重要的理論意義。

近幾十年來，國內(nèi)外學(xué)者們對(duì)含裂隙缺陷的巖體力學(xué)特性進(jìn)行了大量的試驗(yàn)、理論以及數(shù)值模擬研究，取得了豐碩的成果[2-5]。而對(duì)含孔洞缺陷的巖體力學(xué)特性方面的研究相對(duì)較少，Lajtai[6]最早采用石膏模型材料試驗(yàn)研究壓縮作用下含孔洞試件的破壞過程。胡盛斌等[7]采用水泥砂漿材料試驗(yàn)研究了單軸循環(huán)荷載作用下含孔洞試件的破壞特征。宋義敏等[8]采用石膏模型材料試驗(yàn)研究了含孔洞試件拉斷裂和遠(yuǎn)場(chǎng)斷裂破壞、V型坑破壞和分區(qū)破裂化3種破壞形式的形成機(jī)制。李地元等[9]研究了單軸壓縮條件下含雙側(cè)預(yù)制孔洞的花崗巖巖板力學(xué)特性和變形破裂演化規(guī)律。Chen等[10]對(duì)大比尺含圓形孔洞模型進(jìn)行了雙向加載試驗(yàn)。朱譚譚等[11]對(duì)含雙圓形孔洞的砂巖試樣進(jìn)行了單軸壓縮試驗(yàn)?？锥春土严对趲缀涡螒B(tài)上存在著差異性，所以在受力條件下將表現(xiàn)出不同的力學(xué)特性。研究表明，壓縮應(yīng)力作用下含裂隙巖石類脆性材料的破裂演化過程中產(chǎn)生翼型裂紋和次生裂紋[2-5]，翼型裂紋屬于拉裂紋，在裂隙尖端或靠近尖端的位置萌生；次生裂紋屬于剪切裂紋，出現(xiàn)在裂紋尖端。而含孔洞巖石類脆性材料在壓縮加載過程中觀察到3種裂紋模式[8]：初始拉裂紋，形成于孔頂和孔底拉應(yīng)力集中區(qū)；剪切裂紋，形成于孔壁左右兩側(cè)壓應(yīng)力集中區(qū)；遠(yuǎn)場(chǎng)裂紋，形成于孔洞周邊的遠(yuǎn)處。巖體工程中常采用灌漿充填技術(shù)處理地質(zhì)缺陷，加固巖體[12-16]?？锥慈毕莺土严度毕萃瑫r(shí)出現(xiàn)在巖體中，朱譚譚等[17]對(duì)孔洞-裂隙組合型缺陷巖體的力學(xué)特性進(jìn)行了試驗(yàn)研究，但灌漿體對(duì)含孔洞裂隙巖體力學(xué)特性的影響研究卻鮮有報(bào)道。

本文采用類巖石模型材料制備含孔洞裂隙模型試件，單軸壓縮試驗(yàn)對(duì)比孔洞不充填和完全灌漿充填2種工況，使用RFPA2D軟件進(jìn)行數(shù)值模擬，研究不同裂隙傾角條件下灌漿體對(duì)含孔洞裂隙巖體應(yīng)力環(huán)境、破裂模式、力學(xué)特性參數(shù)的影響，為理解孔洞灌漿體的加固作用機(jī)制提供了一種途徑。

1 試件制備與試驗(yàn)方案

1.1 試件制備

若直接采用天然巖石進(jìn)行試驗(yàn)，試件加工困難且無法保證一致性，因此本試驗(yàn)制備類巖石模型材料，這類材料具有與巖石相似的結(jié)構(gòu)特征和破壞特征，同時(shí)還具有力學(xué)特性穩(wěn)定、制作方便以及成本低廉等優(yōu)點(diǎn)。類巖石模型材料配比參考Zhang等的研究[18]，配比為白色硅酸鹽水泥∶細(xì)砂∶水=2∶2∶1 (體積比)。試件尺寸為150 mm (長)×150 mm (寬)×30 mm (厚)。在巖溶地區(qū)隧道工程的修建過程中難免遇到巖溶問題，如果未能及時(shí)處理，容易造成塌方事故，給施工和運(yùn)行造成重大安全隱患。通過地質(zhì)勘察發(fā)現(xiàn)圓形、類圓形是天然溶洞主要的形狀之一[19]。因此，在模型材料中心位置添加圓形預(yù)制孔洞，以模擬天然溶洞，孔洞半徑為25 mm(圖1a)。在模型材料初凝前插入云母片模擬閉合型裂隙，裂隙長度為10 mm，裂隙延長線與孔洞圓心重合。設(shè)計(jì)3種不同的裂隙傾角β，定義為裂隙與水平方向的夾角，分別設(shè)置為0°，30°和90°，分別表征水平、傾斜和豎直裂隙(圖1a)。試件在標(biāo)準(zhǔn)鋼模具中澆注成型，澆筑完24 h后拆模，然后將其放入混凝土養(yǎng)護(hù)箱按標(biāo)準(zhǔn)養(yǎng)護(hù)方式養(yǎng)護(hù)28 d。模型材料和裂隙的基本物理力學(xué)參數(shù)與文獻(xiàn)[18]相同。

試驗(yàn)對(duì)比孔洞不充填和完全灌漿充填2種工況。而孔邊裂隙為閉合型，無法對(duì)其進(jìn)行灌漿處理，所以試驗(yàn)未考慮裂隙灌漿充填這種工況。對(duì)于孔洞完全灌漿充填工況，相應(yīng)的試件制備是在養(yǎng)護(hù)完成的含孔洞裂隙物理模型試件基礎(chǔ)上，按白色硅酸鹽水泥∶細(xì)砂=1∶2 (體積比)以及適量水配制低強(qiáng)度水泥砂漿，用于模擬灌漿體，將其完全充填孔洞(圖1b)，然后再養(yǎng)護(hù)72 h。養(yǎng)護(hù)完成后，將試件上下表面磨平，以便于在試驗(yàn)中觀察裂紋擴(kuò)展過程，物理模型試件如圖2所示。

1.2 試驗(yàn)裝置和方法

單軸壓縮實(shí)驗(yàn)在RYL-600型微機(jī)控制巖石剪切流變儀上完成(圖3)。試驗(yàn)加載過程中采用位移控制，加載速率設(shè)置為0.6 mm/min，試驗(yàn)加載持續(xù)到試件完全破壞。采用配套的控制記錄軟件記錄軸向位移和荷載數(shù)據(jù)，采用高清攝像機(jī)記錄裂紋起裂、擴(kuò)展、貫通和失效等破壞過程。

圖3 試驗(yàn)裝置Fig.3 Experimental setup

圖4 數(shù)值計(jì)算模型Fig.4 Numerical model

2 數(shù)值模型

RFPA2D是基于有限元和統(tǒng)計(jì)損傷理論的數(shù)值模擬軟件，可以真實(shí)地模擬巖石類脆性材料的破壞過程。相應(yīng)的數(shù)值模型試件尺寸為長×寬=150 mm×150 mm，模型網(wǎng)格劃分為300×300=90 000個(gè)細(xì)觀基元(圖4)。破壞判據(jù)采用修正的Morh-Coulomb準(zhǔn)則和拉破壞準(zhǔn)則。模型試件上部施加法向位移控制條件，下部設(shè)置為法向位移約束。

3 試驗(yàn)及模擬結(jié)果分析

3.1 孔洞灌漿體對(duì)應(yīng)力環(huán)境的影響

采用RFPA2D軟件計(jì)算數(shù)值模型試件的應(yīng)力分布，以軸向加載位移量為0.17 mm為例，不同裂隙傾角的試件最大和最小主應(yīng)力云圖如圖5所示，各試件最大壓應(yīng)力和最大拉應(yīng)力的計(jì)算結(jié)果見表1?？梢钥闯觯锥床怀涮詈涂锥赐耆酀{充填試件的主應(yīng)力集中區(qū)域大體一致：“X”型最大壓應(yīng)力帶交于孔洞左右兩側(cè)，且孔洞和裂隙周邊應(yīng)力更大，最大拉應(yīng)力則集中于孔頂和孔底。相比于孔洞不充填試件，完全灌漿充填試件的灌漿體減輕了應(yīng)力集中，表現(xiàn)為最大壓應(yīng)力和最大拉應(yīng)力下降，降幅分別為7.6%～9.8%、0.4%～9.8%。

圖5 不同裂隙傾角β下的主應(yīng)力云圖(MPa)Fig.5 Contours of the principal stresses of different flaw angles β (MPa)

3.2 破裂模式

3.2.1 裂隙傾角對(duì)破裂模式的影響

單軸壓縮作用下不同裂隙傾角的含孔洞裂隙物理模型和數(shù)值模型試件(孔洞不充填和完全灌漿充填)宏觀破裂模式如圖6所示。

(1)水平裂隙(β=0°)：RFPA2D數(shù)值模擬表明，在加載初期，孔頂和孔底產(chǎn)生拉應(yīng)力集中現(xiàn)象(圖5)。這些裂紋立即擴(kuò)展到一定長度，由于拉應(yīng)力得以釋放，繼續(xù)加載后一直保持穩(wěn)定，并未導(dǎo)致試件貫通破壞。隨后，數(shù)值模擬發(fā)現(xiàn)孔壁左右兩側(cè)壓應(yīng)力集中(圖5)，導(dǎo)致發(fā)生碎屑崩落的現(xiàn)象，類似于地下工程開挖過程中發(fā)生的巖爆現(xiàn)象?？锥醋笊蟼?cè)稍遠(yuǎn)處形成遠(yuǎn)場(chǎng)裂紋；預(yù)制裂隙右端形成剪切裂紋沿試件左下端發(fā)展。隨著施加荷載的進(jìn)一步增加，宏觀裂紋搭接貫通，最終導(dǎo)致試件的失穩(wěn)破壞(圖6)，此時(shí)數(shù)值模擬結(jié)果顯示應(yīng)力集中現(xiàn)象不明顯。

表1 應(yīng)力計(jì)算結(jié)果

(2)傾斜裂隙(β=30°)：在加載初期，除了在孔頂和孔底產(chǎn)生初始拉裂紋，孔邊預(yù)制裂隙尖端也觀察到拉裂紋。隨著施加荷載的增加，孔壁左右兩側(cè)碎屑崩落，孔洞左下側(cè)稍遠(yuǎn)處形成遠(yuǎn)場(chǎng)裂紋；而預(yù)制裂隙尖端產(chǎn)生剪切裂紋，向右上端延伸和擴(kuò)展。繼續(xù)加載，遠(yuǎn)場(chǎng)裂紋、預(yù)制裂隙和剪切裂紋搭接貫通，使得試件產(chǎn)生最終破壞(圖6)，此后數(shù)值模擬結(jié)果表明應(yīng)力趨于平均分布。

圖6 不同裂隙傾角β下的試件單軸壓縮破裂模式Fig.6 Uniaxial compression failure modes of specimens under different flaw angles β

(3)豎直裂隙(β=90°)：孔洞左上側(cè)和右下側(cè)形成的遠(yuǎn)場(chǎng)裂紋導(dǎo)致試件的失穩(wěn)破壞(圖6)。

RFPA2D數(shù)值模擬所再現(xiàn)的破壞過程與試驗(yàn)結(jié)果大體吻合。本次試驗(yàn)中水平和傾斜裂隙試件的破裂模式與朱譚譚等[16]的含孔洞張開型裂隙砂巖試驗(yàn)結(jié)果基本一致。由于裂隙開合狀態(tài)以及試件材料的不同，導(dǎo)致豎直裂隙試件的破裂模式有所不同，探究這兩因素的影響效應(yīng)是今后研究的方向。

3.2.2 孔洞灌漿體對(duì)破裂模式的影響

孔洞完全灌漿充填試件的裂紋擴(kuò)展特征與孔洞不充填試件類似，這是由于灌漿體強(qiáng)度遠(yuǎn)低于圍巖強(qiáng)度，雖然使得應(yīng)力環(huán)境有所改善，但仍不足以改變最終破裂模式。水平和豎直裂隙試件(β=0°和90°)最終的破裂呈現(xiàn)為遠(yuǎn)場(chǎng)裂紋剪斷巖橋，均未沿預(yù)制裂隙發(fā)生破壞；而沿遠(yuǎn)場(chǎng)裂紋、預(yù)制裂隙和次生剪切裂紋的剪切滑移導(dǎo)致傾斜裂隙試件的最終破壞。

3.3 力學(xué)特性參數(shù)

3.3.1 裂隙傾角對(duì)力學(xué)特性參數(shù)的影響

單軸壓縮作用下不同裂隙傾角的含孔洞裂隙試件(孔洞不充填和完全灌漿充填)力學(xué)特性參數(shù)見圖7。對(duì)于含孔洞裂隙試件，預(yù)制裂隙傾角顯著影響試件的破裂模式，而宏觀強(qiáng)度參數(shù)與巖體中裂紋起裂、擴(kuò)展和貫通過程緊密相關(guān)。傾斜裂隙試件破壞主要受預(yù)制裂隙控制，破壞模式為沿裂隙面滑移破壞，導(dǎo)致試件峰值強(qiáng)度最低，試驗(yàn)和模擬結(jié)果分別為7.20和7.56 MPa。遠(yuǎn)場(chǎng)裂紋是水平和豎直裂隙這兩種結(jié)構(gòu)型式試件破壞的重要表現(xiàn)形式，均呈現(xiàn)為穿切巖橋剪切破壞，造成其峰值強(qiáng)度高于傾斜裂隙試件。

圖7 灌漿體對(duì)力學(xué)特征參數(shù)的影響Fig.7 Influence of grouting body on mechanical characteristic parameters

3.3.2 孔洞灌漿體對(duì)力學(xué)特性參數(shù)的影響

(1)對(duì)于峰值強(qiáng)度變化特征，蘇海健等[13]、尹乾等[15]研究發(fā)現(xiàn)低強(qiáng)度充填物對(duì)裂隙巖體峰值強(qiáng)度的增強(qiáng)效應(yīng)并不明顯，與不充填相比僅增大1.49%～ 3.96%；而本次試驗(yàn)采用的孔洞灌漿體也屬于低強(qiáng)度材料，但相比于孔洞不充填試件，孔洞完全灌漿充填的物理模型和數(shù)值模型試件的峰值強(qiáng)度增幅分別為5.45%～23.33%、4.18%～14.29%，加固效果明顯。另外，還發(fā)現(xiàn)傾斜裂隙試件的峰值強(qiáng)度增幅最大。

(2)對(duì)于殘余強(qiáng)度變化特征，孔洞完全灌漿充填試件破壞后并沒有完全失去承載力，除了巖石原生和次生裂隙內(nèi)部之間摩擦力，灌漿體在受力條件下表現(xiàn)出彈塑性特征，所以還能承受荷載，導(dǎo)致殘余強(qiáng)度高于孔洞不充填試件。孔洞完全灌漿充填的物理模型和數(shù)值模型試件殘余強(qiáng)度的增幅分別達(dá)到31.68%～161.71%、22.54%～73.85%。在這一點(diǎn)上，孔洞灌漿體加固效果優(yōu)于裂隙充填工況，蘇海健等[13]、尹乾等[15]研究發(fā)現(xiàn)殘余強(qiáng)度有所增大。這是因?yàn)槌涮钗镎急炔淮?，所以影響輕微。

(3) 對(duì)于變形變化特征，孔洞完全灌漿充填情況下物理模型試件的彈性模量比孔洞不充填試件減小了1.47%～20.26%，而數(shù)值模型試件的彈性模量基本保持不變。

3.4 孔洞灌漿體的加固作用機(jī)制

綜合上述試驗(yàn)和模擬結(jié)果，孔洞灌漿體的加固作用可能存在以下3種機(jī)制：

(1)當(dāng)孔洞灌漿處理后，灌漿體與含孔洞裂隙巖體形成了一種復(fù)合結(jié)構(gòu)體，盡管灌漿體強(qiáng)度不高，但巖體受荷變形壓縮灌漿體，灌漿體又對(duì)圍巖產(chǎn)生反作用力，灌漿體與含孔洞裂隙巖體之間形成了力學(xué)上相互作用，改善了孔洞周邊的應(yīng)力環(huán)境，增強(qiáng)了圍巖的自支撐能力。

(2)另一方面，灌漿體和圍巖交界面會(huì)產(chǎn)生局部支護(hù)力，尤其在孔邊裂隙附近，而已有研究證明小的表面荷載對(duì)摩擦型介質(zhì)中屈服區(qū)范圍可能會(huì)產(chǎn)生重大影響[20]。

(3) 另外，圍巖破壞并不意味著承載能力的消失，灌漿體在受力條件下表現(xiàn)出彈塑性特征，改善了含孔洞裂隙巖體的峰后特征，破壞后的試件還有較高的殘余強(qiáng)度。

4 結(jié)論

(1)孔洞灌漿體能夠改善孔洞周圍的應(yīng)力環(huán)境，應(yīng)力分析表明灌漿體與含孔洞裂隙巖體之間形成了力學(xué)上的相互作用，灌漿充填后的試件最大壓應(yīng)力和最大拉應(yīng)力均小于孔洞不充填試件，降幅分別為7.6%～9.8%和0.4%～9.8%，使得孔洞上下兩端初始拉裂紋以及孔壁碎屑崩落引起的破壞程度降低。

(2)孔洞完全灌漿充填的物理模型和數(shù)值模型試件最終破裂模式與孔洞不充填試件類似，傾斜裂隙試件破壞主要受預(yù)制裂隙控制，而水平和豎直裂隙試件的破裂模式為巖橋剪斷。

(3)孔洞完全灌漿充填的物理模型和數(shù)值模型試件峰值強(qiáng)度和殘余強(qiáng)度均顯著高于孔洞不充填試件，峰值強(qiáng)度增幅分別為5.45%～23.33%、4.18%～ 14.29%，殘余強(qiáng)度增幅分別為31.68%～ 161.71%、22.54%～73.85%。灌漿充填后的力學(xué)參數(shù)變化幅度與裂隙傾角有關(guān)，傾斜裂隙試件的加固效果最好。

(4)本次研究初步探究了單軸壓縮作用下充填試件的力學(xué)特性，下一步將開展雙軸和三軸壓縮作用下不同孔洞形狀的相關(guān)試驗(yàn)和數(shù)值模擬研究工作。