吳樂文, 趙奎
(1.廣西華錫集團股份有限公司,廣西 柳州 545006;2.江西理工大學(xué)資源與環(huán)境工程學(xué)院,江西 贛州341000)
許多國內(nèi)外學(xué)者針對控制著工程巖體穩(wěn)定性的巖體結(jié)構(gòu)面力學(xué)特性這一熱點問題進行了持續(xù)性研究.眾所周知,結(jié)構(gòu)面的存在對巖體力學(xué)性質(zhì)有著決定性影響,因此,研究巖體結(jié)構(gòu)面力學(xué)性質(zhì)是非常有意義的.
孫廣忠教授在文獻[1]中歸納了結(jié)構(gòu)面的4種起伏形態(tài)分別為平直結(jié)構(gòu)面、臺階狀結(jié)構(gòu)面、鋸齒狀結(jié)構(gòu)面、波狀結(jié)構(gòu)面,4種結(jié)構(gòu)面的力學(xué)效應(yīng)是不同的.除了鋸齒狀結(jié)構(gòu)面,其他3種結(jié)構(gòu)面力學(xué)特性的研究已較多,一些學(xué)者對巖體結(jié)構(gòu)面錨固特性采用室內(nèi)試驗或理論分析進行了相關(guān)研究[2-4].但是對結(jié)構(gòu)面注漿加固前后力學(xué)特性的研究還比較少,尤其是采用數(shù)值分析方法的研究更少.因此,本文通過3種類型和2種臺階高度結(jié)構(gòu)面剪切模擬試驗,對結(jié)構(gòu)面注漿前后抗剪特性進行了分析研究.該研究對揭示其加固機制具有重要的意義.
通過ANSYS有限元分析軟件建立數(shù)值分析模型,再通過轉(zhuǎn)換導(dǎo)入FLAC3D軟件中進行分析和計算[5-12].為盡可能得到巖石結(jié)構(gòu)面剪切特性的規(guī)律,模型尺寸:長×寬×高為 100 mm×100 mm×100 mm.為了模擬結(jié)構(gòu)面,在模型中間設(shè)置厚度為2 mm的結(jié)構(gòu)面模擬其注漿前后效果.設(shè)置的結(jié)構(gòu)面類型包括1個臺階、2個臺階以及組合(1個臺階與1個鋸齒)3種類型,臺階寬度為20 mm,每種類型包括2種臺階高度,分別為10 mm、6 mm,則3種類型結(jié)構(gòu)面模型高寬比分別為0.5和0.3.如圖1~圖3所示分別為3種結(jié)構(gòu)面模型示意圖.
圖1 1個臺階結(jié)構(gòu)面模型
圖2 2個臺階結(jié)構(gòu)面模型
圖3 組合結(jié)構(gòu)面模型
邊界條件設(shè)置為底部固定約束,下半部分模型四周約束水平方向位移,上半部分頂部為自由邊界.
計算采用Mohr-Coulomb準(zhǔn)則,在模型頂部施加法向荷載,上半部分巖石施加水平位移荷載,荷載大小2×10-4mm/時步.模擬試驗中只考慮法向荷載,因為法向荷載遠遠大于巖石自重,所以計算略去巖石自重[13-15].如表1、表2所示,巖石和結(jié)構(gòu)面參數(shù)是根據(jù)RMT-150C型巖石力學(xué)試驗系統(tǒng)對巖石試塊進行室內(nèi)單軸和剪切試驗的結(jié)果選取.
表1 巖石參數(shù)選取
表2 結(jié)構(gòu)面參數(shù)選取
通過數(shù)值計算分析討論臺階狀結(jié)構(gòu)面剪切破壞形式以及規(guī)律,探討臺階數(shù)量及結(jié)構(gòu)面形態(tài)與抗剪強度關(guān)系、施加法向應(yīng)力大小與抗剪強度關(guān)系.由于試驗結(jié)果圖太多,此處只列出圖4~圖6,分別為試驗中高寬比0.5和0.3時,結(jié)構(gòu)面最大、最小主應(yīng)力云圖以及組合剖面塑性云圖.本文中所有圖示除非特別說明,剪切方向均為從左往右.
通過分析圖4、圖5、圖6等1個臺階、2個臺階、組合結(jié)構(gòu)面模型剖面的最大主應(yīng)力、最小主應(yīng)力以及塑性區(qū)云圖得到下述研究結(jié)果.
1)拉應(yīng)力區(qū)域皆形成于3種結(jié)構(gòu)面的臺階及鋸齒下部,拉應(yīng)力最大處位于底部.
2)結(jié)構(gòu)面上半部分巖石在剪切作用下有往上抬起效應(yīng),導(dǎo)致臺階右邊下部巖石受到明顯壓應(yīng)力作用,臺階右邊底角轉(zhuǎn)折處形成壓應(yīng)力集中區(qū)域.
3)組合類型巖石在剪切作用下鋸齒下方和鋸齒頂部均表現(xiàn)為拉應(yīng)力.而壓應(yīng)力較集中區(qū)域為鋸齒正下方與臺階底部,且在鋸齒右面底部達到最大壓應(yīng)力,說明剪切過程上部巖石有往上抬起效應(yīng).
圖4 1個臺階剖面最大主應(yīng)力云圖
圖5 2個臺階剖面最小主應(yīng)力云圖
圖6 組合剖面塑性區(qū)云圖
4)由塑性云圖6可知,因巖石抗拉強度遠遠小于巖石抗壓強度,導(dǎo)致臺階底部在剪切應(yīng)力和法向應(yīng)力的作用下最先表現(xiàn)為拉伸破壞.組合類型模型中鋸齒和臺階的底部都已開始出現(xiàn)拉伸破壞,且最大拉應(yīng)力區(qū)域與破壞區(qū)域一致,臺階最終因張剪滑移而破壞.
為了對比分析臺階高寬比不同的情況下3種結(jié)構(gòu)面剪切破壞規(guī)律,進行了臺階高寬比為0.3時結(jié)構(gòu)面剪切試驗.
通過分析臺階高寬比為0.3的1個臺階結(jié)構(gòu)面模型剖面的最大主應(yīng)力、最小主應(yīng)力以及塑性區(qū)云圖可知:
1)臺階高度降低后,臺階高寬比為0.3時結(jié)構(gòu)面主應(yīng)力云圖與臺階高寬比為0.5時的應(yīng)力云圖規(guī)律是一樣的.臺階底部受到拉應(yīng)力作用,拉應(yīng)力集中區(qū)域出現(xiàn)在臺階右下角.
2)沿順剪方向,臺階只在左下部分出現(xiàn)斜向上的拉伸破壞痕跡,表明當(dāng)臺階高寬比較小時,臺階一般因剪切作用產(chǎn)生壓切破壞,且對于2個臺階結(jié)構(gòu)面模型壓切破壞同時發(fā)生于2個臺階.因此,當(dāng)臺階高寬比較小時,臺階因剪切作用而產(chǎn)生壓切破壞.破壞區(qū)域分布與文獻[1]分析一致.
3)組合結(jié)構(gòu)面模型中拉應(yīng)力集中區(qū)域位于臺階和鋸齒的底部,臺階底部受彎矩作用明顯,由此可見:高寬比越大,臺階狀結(jié)構(gòu)面張性破壞效應(yīng)越明顯.鋸齒底部的拉應(yīng)力明顯小于臺階底部的拉應(yīng)力;鋸齒高度降低后,鋸齒下部的拉應(yīng)力集中效應(yīng)降低了,由此可見:起伏角α越大,結(jié)構(gòu)面張性破壞效應(yīng)越明顯.
4)2個臺階結(jié)構(gòu)面與1個臺階結(jié)構(gòu)面模型類似,沿受剪方向,出現(xiàn)于臺階左下角的拉伸破壞區(qū)逐步向臺階頂端擴展,最終導(dǎo)致臺階形成壓切破壞.鋸齒下方拉伸破壞區(qū)在左下與右下部位明顯.
通過對3種類型結(jié)構(gòu)面注漿前后的試驗過程監(jiān)測得到了法向應(yīng)力為 1 MPa、2 MPa、3 MPa、4 MPa、5 MPa下峰值剪切應(yīng)力的大小,如表3所示.
表3 不同法向應(yīng)力對應(yīng)峰值剪切應(yīng)力/MPa
通過數(shù)值模擬驗證了法向應(yīng)力與結(jié)構(gòu)面抗剪強度的關(guān)系特性.通過數(shù)值計算結(jié)果可知:抗剪強度和法向應(yīng)力擬合相關(guān)系數(shù)近似為1.0.結(jié)構(gòu)面抗剪強度隨法向應(yīng)力增大而呈線性增長.
依據(jù)文獻[4]的研究可知,數(shù)值模擬中等效參數(shù)方法可以模擬結(jié)構(gòu)面注漿效果,因此,通過改變中間層參數(shù)模擬注漿效果,結(jié)果如圖7~圖9所示.
圖7所示為1個臺階結(jié)構(gòu)面注漿與未注漿結(jié)構(gòu)面抗剪強度對比條形圖,注漿后結(jié)構(gòu)面抗剪強度在5個等級法向應(yīng)力下提高幅度明顯,分別達到14.11%、10.49%、8.60%、8.64%、8.70%.注漿后強度平均提高了10.11%.
圖8所示為2個臺階結(jié)構(gòu)面注漿與未注漿結(jié)構(gòu)面抗剪強度對比條形圖,注漿后結(jié)構(gòu)面抗剪強度在5個等級法向應(yīng)力下提高幅度明顯,分別為11.78%、8.33%、13.24%、6.68%、5.93%.注漿后強度平均提高了9.19%.
圖9所示為組合結(jié)構(gòu)面注漿與未注漿抗剪強度對比條形圖,注漿后結(jié)構(gòu)面抗剪強度在5個等級法向應(yīng)力下提高幅度明顯,分別為24.09%、7.03%、18.80%、16.20%、15.70%.注漿后強度平均提高了16.35%.
圖7 1個臺階結(jié)構(gòu)面注漿與未注漿強度對比
圖8 2個臺階結(jié)構(gòu)面注漿與未注漿強度對比
圖9 組合結(jié)構(gòu)面注漿與未注漿強度對比
通過對圖7~圖9的分析,進一步確定巖石結(jié)構(gòu)面注漿效果的等效參數(shù)模擬方法是可行的.
表4為5個等級法向應(yīng)力下1個臺階模型和組合結(jié)構(gòu)面模型抗剪強度大小關(guān)系.
表4 不同法向應(yīng)力2種類型結(jié)構(gòu)面抗剪強度比較/MPa
由表4可以看出,組合結(jié)構(gòu)面模型相對于1個臺階結(jié)構(gòu)面在增加一個鋸齒的情況下,結(jié)構(gòu)面抗剪強度并沒有明顯的提高,說明大爬坡角對結(jié)構(gòu)面抗剪強度大小起控制作用,定義其為控制性爬坡角;小爬坡角對結(jié)構(gòu)面抗剪強度大小不起控制性作用,定義其為非控制性爬坡角.
由表4中結(jié)構(gòu)面抗剪強度增長率可以看出,組合情況下臺階(即爬坡角為90°)結(jié)構(gòu)面模型對結(jié)構(gòu)面抗剪強度起決定性作用,鋸齒對抗剪強度影響較小,平均增長率僅為2.57%.
由此,結(jié)合前文分析,可以提出并得到組合情況下巖石結(jié)構(gòu)面抗剪強度修正公式:
式(1)中:ki為非控制性爬坡角因子,以本文中鋸齒狀結(jié)構(gòu)面模型爬坡角為30°情況下,其非控制性爬坡角因子ki=0.025 7.
在開展室內(nèi)單軸、剪切試驗的基礎(chǔ)上,開展了結(jié)構(gòu)面數(shù)值模擬試驗.通過數(shù)值計算,分析了3種類型結(jié)構(gòu)面抗剪力學(xué)特性,得出結(jié)論如下:
1)在剪切試驗中,高寬比較大時,結(jié)構(gòu)面一般發(fā)生剪斷破壞;高寬比較小時,結(jié)構(gòu)面一般發(fā)生壓切破壞.
2)結(jié)構(gòu)面抗剪強度隨著法向應(yīng)力的增加近似呈線性增長.
3)通過設(shè)置中間層,采用等效參數(shù)模擬方法進行結(jié)構(gòu)面注漿模擬試驗是可行的.數(shù)值模擬結(jié)果表明,注漿后結(jié)構(gòu)面抗剪強度得到顯著提高,從而在微觀上揭示了采用注漿加固巖石結(jié)構(gòu)面的可行性.
4)以數(shù)值試驗為依據(jù)提出了非控制性爬坡角因子ki以及組合結(jié)構(gòu)面抗剪強度修正公式.計算出本文中組合結(jié)構(gòu)面的非控制性爬坡角因子為0.025 7.因此,可以推斷出同一巖體中結(jié)構(gòu)面抗剪強度大小主要由其中較大的爬坡角控制.
[1]杜時貴.巖體結(jié)構(gòu)面的工程性質(zhì)[M].北京:地震出版社,1999.
[2]白紅杰,陸文,肖正學(xué),等.巖石結(jié)構(gòu)面抗剪強度方向性模擬研究[J].礦業(yè)研究與開發(fā),2007,10(5):61-63.
[3]林杭,曹平,周正義.FLAC3D模擬全長注漿錨桿的作用效果[J].巖土力學(xué),2005,26(增刊 1):167-170.
[4]程東幸,潘煒,劉大安,等.錨固節(jié)理巖體等效力學(xué)參數(shù)三維離散元模擬[J].巖土力學(xué),2006,27(12):2128-2132.
[5]王中文,方建勤.巖體結(jié)構(gòu)面錨固剪切特性的數(shù)值分析[J].煤炭學(xué)報,2010,35(5):729-733.
[6]韓立軍,宗義江,韓貴雷,等.巖石結(jié)構(gòu)面注漿加固抗剪特性試驗研究[J].巖土力學(xué),2011,32(9):2570-2576.
[7]楊明.巖體結(jié)構(gòu)面力學(xué)特性及其錨固效應(yīng)的數(shù)值計算研究[D].長沙:中南大學(xué),2009.
[8]劉自由,江學(xué)良,林杭.軟弱結(jié)構(gòu)面加樁特性的數(shù)值分析[J].中南大學(xué)學(xué)報(自然科學(xué)版):2011,42(5):1461-1466.
[9]李慶勇,遠方,童中明.混凝土剪切破壞數(shù)值模擬[J].低溫建筑技術(shù),2011(5):51-53.
[10]周蓮君,彭振斌,何忠明,等.結(jié)構(gòu)面剪切特性的試驗與數(shù)值模擬分析[J].科技導(dǎo)報,2009,27(4):31-35.
[11]李大勇,潘軍剛.土的直接剪切試驗三維數(shù)值模擬研究[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2008,27(6):16-20.
[12]胡黎明,馬杰,張丙?。奔粼囼炛薪佑|面漸進破壞的數(shù)值模擬[J].清華大學(xué)學(xué)報(自然科學(xué)版),2008,48(6):943-946.
[13]吳樂文.巖石結(jié)構(gòu)面注漿前后力學(xué)特性研究[D].贛州:江西理工大學(xué),2012.
[14]楊松林,徐衛(wèi)亞,朱煥春.錨桿在節(jié)理中的加固作用[J].巖土力學(xué),2002,23(5):604-607.
[15]金解放,鐘海兵,吳越等.靜載荷與循環(huán)沖擊作用下巖石損傷變量定義方法的選擇[J].有色金屬科學(xué)與工程,2013,4(4):85-90.