肖 旺,蘇永華,方硯兵
(湖南大學(xué) 土木工程學(xué)院,湖南 長沙 410082)
?
考慮峰后特性的隧道圍巖錨固力學(xué)效應(yīng)
肖旺,蘇永華,方硯兵
(湖南大學(xué)土木工程學(xué)院,湖南長沙410082)
圍繞處于峰后狀態(tài)的隧道圍巖在全長黏結(jié)錨桿支護(hù)下的力學(xué)機(jī)制進(jìn)行研究。按力學(xué)特征,將在次生應(yīng)力作用下的峰后狀態(tài)圍巖由隧道臨空面往深處依次劃分為殘余強(qiáng)度區(qū)、塑性軟化區(qū)及彈性變形區(qū)?;阱^桿黏結(jié)抗滑移變形支護(hù)抗力產(chǎn)生原理,導(dǎo)得錨桿軸向應(yīng)力及抗剪強(qiáng)度方程。將全長黏結(jié)錨桿通過巖石-錨固面對圍巖提供的支護(hù)力等效為圍巖單元徑向體積力,建立了考慮錨桿支護(hù)效應(yīng)的微元體平衡方程;將該平衡方程分別與考慮峰后體積擴(kuò)容及軟化特征的殘余帶、軟化帶塑性相容方程聯(lián)合,推出了隧道峰后圍巖殘余區(qū)、軟化區(qū)在注漿錨桿支護(hù)下位移、應(yīng)力及范圍的計算公式?;诶碚撗芯康膶?shí)例分析,進(jìn)一步揭示了注漿錨桿與圍巖內(nèi)部錨固效應(yīng),錨桿端部預(yù)緊力對峰后圍巖軟化、擴(kuò)容性質(zhì)的改善作用和對圍巖各區(qū)范圍發(fā)展及隧道周邊變形的控制作用。
隧道工程;隧道圍巖;解析法;軟巖擴(kuò)容;錨固效應(yīng);應(yīng)變軟化;界面剪應(yīng)力;預(yù)緊力
在隧道開挖過程中應(yīng)力集中如果超過巖石的峰值強(qiáng)度,會導(dǎo)致巖石峰后軟化和峰后體積擴(kuò)容[1],前者導(dǎo)致圍巖承載能力下降,而后者是導(dǎo)致圍巖大變形的主要原因。圓形隧洞圍巖的變形與壓力的彈塑性力學(xué)分析最早是由Fenner提出來的,后來Kastner[2]作了重要修正,但他們都假設(shè)隧洞圍巖為理想彈塑性介質(zhì),且認(rèn)為巖體破壞后無體積變化,這與實(shí)際情況具有一定偏差,特別是軟弱或破碎巖體隧道的開挖必將導(dǎo)致圍巖應(yīng)變軟化和非線性體積膨脹,產(chǎn)生擴(kuò)容變形。研究結(jié)果表明考慮圍巖應(yīng)變軟化和塑性擴(kuò)容更接近實(shí)際[3-11]。蘇永華[3]、孫闖[4]基于深部圍巖在開挖卸載應(yīng)力擾動之后其力學(xué)性能弱化的現(xiàn)實(shí),認(rèn)為其應(yīng)力-應(yīng)變關(guān)系服從峰后軟化的節(jié)理模型。袁文伯[5]、馬念杰[6]考慮巖體的應(yīng)變軟化和塑性流動特性,以巖體的應(yīng)力-應(yīng)變峰后破壞段作為巖體的塑性軟化強(qiáng)度,建立軟化模型,得到了更符合實(shí)際、適用性更廣泛的巷道圍巖塑性區(qū)半徑和應(yīng)力的一般解。付國彬[7]、姚國圣[8]、Brown E T和Hoek E[9]同時考慮了圍巖應(yīng)變軟化和擴(kuò)容特性(體積膨脹),基于應(yīng)力-應(yīng)變?nèi)€段模型,得到了隧道圍巖彈塑性解答;范文[10]、Wang[11]基于非關(guān)聯(lián)彈塑性準(zhǔn)則,引入剪脹角模擬軟巖的塑性擴(kuò)容特性,研究了圍巖擴(kuò)容的變形、應(yīng)力解析解。以上成果顯著,較接近工程實(shí)際。
錨桿在巖體中的作用機(jī)制一直是國內(nèi)外研究的主要問題,并取得了大量研究成果[12-15]。文競舟[13]、C.Li和B.Stillborg[14]、Y.Cai[15]等在不考慮巖體塑性軟化和體積擴(kuò)容特性的條件下研究了注漿錨桿與圍巖的相互作用及圍巖的變形、應(yīng)力解析。
研究表明,在隧道圍巖變形的4個力學(xué)分區(qū)[1]中,塑性強(qiáng)化區(qū)處于峰前變形階段,是圍巖承載的主體;塑性軟化區(qū)和塑性殘余區(qū)處于峰后變形階段,巖體強(qiáng)度已經(jīng)損傷,是實(shí)施支護(hù)的主要對象。錨桿支護(hù)主要通過外支撐和內(nèi)加固的形式作用于圍巖。外支撐提供徑向抗力施加圍壓,圍壓的增大會使擴(kuò)容量隨之減弱,因此可以減少體積擴(kuò)容。內(nèi)加固一方面通過提高殘余區(qū)力學(xué)參數(shù),降低圍巖軟化模量;另一方面將錨桿界面的剪應(yīng)力以體積力[13]的形式引入圓形隧道圍巖中。本文在前人研究的基礎(chǔ)上考慮峰后圍巖變形的應(yīng)變軟化及塑性剪脹擴(kuò)容特性,基于Mohr-Coulomb強(qiáng)度準(zhǔn)則和非關(guān)聯(lián)彈塑性準(zhǔn)則,并結(jié)合錨桿的加固機(jī)理[16],對圓形洞室峰后圍巖在錨桿支護(hù)作用下進(jìn)行彈塑性分析,得到了更為全面的應(yīng)力及變形解答。
1.1基本力學(xué)模型
洞室開挖后力學(xué)模型如圖1所示:①洞室開挖半徑為r0,軟化區(qū)半徑Rp和殘余區(qū)半徑Rb;②錨桿有效錨固長度為L,錨固半徑L0=L+r0;③圍巖原巖應(yīng)力場為σz;④把所有作用在圍巖表面上的支護(hù)阻力pi稱為外部支撐,把錨桿與圍巖相互黏結(jié)的作用力稱為內(nèi)部加固,通過錨桿端部預(yù)緊力傳遞。
圖1 圍巖彈塑性變形區(qū)域圖Fig.1 Elastic-plastic deformation zone of surrounding rock
σr,σθ分別為隧道圍巖的徑向、切向應(yīng)力;εr,εθ分別為隧道圍巖的徑向、切向應(yīng)變。文中上角標(biāo)為“e”,“p”,“b”分別為彈性區(qū)、塑性軟化區(qū)和塑性殘余區(qū)的量。
基本假定:(1)認(rèn)為巖體為各向同性、均質(zhì)連續(xù)介質(zhì),無限長圓形隧道可按軸對稱平面應(yīng)變問題處理;(2)錨桿和注漿錨固劑均處于彈性狀態(tài),錨桿與圍巖完全黏結(jié),不產(chǎn)生滑移。認(rèn)為注漿體和巖體力學(xué)性質(zhì)一致,只考慮錨桿與圍巖的相互作用;(3)在塑性軟化區(qū),認(rèn)為巖體強(qiáng)度隨變形發(fā)展而衰減主要是由于內(nèi)聚力變化的結(jié)果[1,5],因此只考慮內(nèi)聚力的軟化和強(qiáng)化,不考慮內(nèi)摩擦角的軟化和強(qiáng)化[7-9]。
1.2圍巖彈塑性軟化與殘余模型
根據(jù)巖體全應(yīng)力-應(yīng)變曲線,如圖2(a)所示,其變形可分為彈性區(qū)、塑性強(qiáng)化區(qū)、塑性軟化區(qū)和塑性殘余區(qū),為了簡化,把彈性區(qū)和塑性強(qiáng)化區(qū)看成峰前彈性變形階段[1]。
圖2 圍巖峰后軟化模型Fig.2 Surrounding rock post-peak softening model
巖石峰后應(yīng)變軟化定律簡化為理想殘余塑性模型[7-11],見圖2(b);由研究可知[1,7-9]凝聚力c會隨著塑性應(yīng)變的增加而逐漸減小至殘余值。假設(shè)c隨塑性切應(yīng)變線性軟化[1,7-9],見圖2(c)。
根據(jù)圖2(b),采用線性Mohr-Coulomb準(zhǔn)則[1,7-10]:
(1)
式中,Kp=(1+sinφ)/(1-sinφ),φ為巖體內(nèi)摩擦角;σc為巖體峰值強(qiáng)度,σc=2c0cosφ/(1-sinφ)。
在塑性軟化區(qū),式(1)變?yōu)椋?/p>
(2)
(3)
令軟化系數(shù)k=Mc,在圍巖殘余區(qū),根據(jù)式(1)有:
(4)
1.3圍巖彈塑性擴(kuò)容模型
考慮巖體峰后塑性軟化區(qū)和殘余區(qū)巖體發(fā)生擴(kuò)容,擴(kuò)容系數(shù)與應(yīng)變關(guān)系[1,7-9]如圖3所示。
圖3 圍巖峰后擴(kuò)容模型Fig.3 Surrounding rock post-peak dilatancy model
考慮巖體擴(kuò)容的非關(guān)聯(lián)流動法則[7-11],在塑性軟化區(qū)有:
(5)
在塑性殘余區(qū)有:
(6)
1.4錨桿力學(xué)分析
1.4.1錨桿作用機(jī)理
在隧洞圍巖中施加錨桿,錨桿對圍巖的錨固作用體現(xiàn)在徑向和切向的錨固力。徑向錨固力對圍巖施加圍壓,圍壓增大使塑性擴(kuò)容降低,穩(wěn)定性增強(qiáng)。錨桿貫穿隧道圍巖的軟弱面,切向錨固力改善軟弱面的力學(xué)性質(zhì),使圍巖抗剪強(qiáng)度參數(shù)c,φ提高,殘余強(qiáng)度增大,進(jìn)而降低圍巖軟化模量。
根據(jù)文獻(xiàn)[16]錨桿作用機(jī)理,在圍巖殘余區(qū)
(7)
式中,cm為施加錨桿后殘余區(qū)黏聚力增大值;σs為錨桿屈服強(qiáng)度;ds為錨桿直徑;Sc,Sl分別為錨桿沿隧道縱向和橫向的布置間距。
施加錨桿后的巖體強(qiáng)度為:
(8)
施加錨桿后圍巖軟化模量為:
(9)
1.4.2全長黏結(jié)錨桿界面剪應(yīng)力模型
沿錨桿徑向的剪應(yīng)力是由于洞室開挖,圍巖變形引起的。通過大量的拉拔試驗[14-15]表明,圍巖變形產(chǎn)生剪應(yīng)力的主要原因是錨桿和圍巖之間的相對位移。從測試結(jié)果[17-18]可以看出,當(dāng)接觸面上沒有滑移時,即錨桿與巖體完全黏結(jié),錨桿與圍巖界面剪應(yīng)力與剪切位移呈線性增加關(guān)系。圖4為錨桿受力圖。
圖4 錨桿受力圖Fig.4 Stress analysis of bolt
如圖4所示,取錨桿軸向微元段,荷載以圖示向右為正,τ(r)為錨桿界面剪應(yīng)力分布,P(r)為錨桿軸力。假定錨固體表面剪應(yīng)力與剪切位移呈線彈性關(guān)系[17-18],則有:
(10)
式中,r為錨桿沿隧道半徑方向的坐標(biāo);ε(r)為錨桿的軸向正應(yīng)變;K為錨固體與錨固層界面的拉拔剪切剛度[14-15,17-18]。
(11)
式中,Er,Eg分別為圍巖和灌漿的彈性模量;νr,νg分別為圍巖和灌漿的泊松比。
剪應(yīng)力在錨桿長度上的積分即為軸力分布:
(12)
軸向拉應(yīng)變
(13)
式中Es為錨桿的彈性模量。
聯(lián)立式(10)、(12)和(13)有
(14)
利用邊界條件:P(r)|r=r0=P0;P(r)|r=L0=0解得:
(15)
(16)
式中,D1=-P0emr0/(e2mL0-e2mr0);D2=P0em(r0+2L0)/ (e2mL0-e2mr0),P0為錨桿端部軸力,即端部預(yù)緊力,m2=4K/(Esds)。
1.5體積力計算模型
把錨桿對圍巖的支護(hù)作用力以等效徑向體積力f(r)作用于圓形隧道圍巖,假設(shè)錨桿沿隧道斷面對稱分布[13],即可分析僅含有單根錨桿的圍巖楔形單元體,如圖5所示。
圖5 錨桿與圍巖楔形單元Fig.5 Wedge element of bolt and surrounding rock
如圖5所示,取錨桿微段dr界面上的合力dQ為
(17)
此微段的體積dV為:
dV=Scrαdr,
(18)
(19)
式中Sc,α分別為錨桿沿隧道縱向間距和環(huán)向夾角。
將剪應(yīng)力表達(dá)式(16)代入式(19)即可得到微元體積力f(r)的表達(dá)式。
軸對稱錨固區(qū)內(nèi)圍巖應(yīng)滿足平衡微分方程:
(20)
幾何方程:
(21)
2.1洞室圍巖應(yīng)力、位移的分析
2.1.1彈性區(qū)(Rp≤r<+∞)
(1)彈性區(qū)應(yīng)力[1-2]
(22)
式中,σRp為彈塑性交界處的徑向應(yīng)力,r=Rp時,σRp=(2σz-σc)/(1+Kp)。
(2)彈性區(qū)變形與位移[1-2]
(23)
式中,C=(σzsinφ+ccosφ)(1+μ)/E,E為巖石彈性模量,μ為泊松比。
2.1.2塑性軟化區(qū)(Rb≤r