王華, 李智遠(yuǎn)
(1.貴州省質(zhì)安交通工程監(jiān)控檢測中心有限責(zé)任公司, 貴州 貴陽 550000;2. 中南大學(xué) 土木工程學(xué)院, 湖南 長沙 410004)
預(yù)應(yīng)力錨索框架梁是巖土錨固工程中非常重要的分支,以其安全可靠、加固合理等優(yōu)點(diǎn)被廣泛應(yīng)用于邊坡治理。在預(yù)應(yīng)力錨索框架服役過程中,受運(yùn)營環(huán)境和管理水平影響,其預(yù)應(yīng)力會發(fā)生損失,且支護(hù)結(jié)構(gòu)會發(fā)生老化,使錨固效果降低,嚴(yán)重時會導(dǎo)致預(yù)應(yīng)力錨索框架失效而引起邊坡失穩(wěn)。預(yù)應(yīng)力錨索框架內(nèi)力計算及預(yù)應(yīng)力錨索框架加固邊坡穩(wěn)定性分析是當(dāng)前主要研究方向,如朱晗迓等討論了錨墩框架梁系統(tǒng)的受力機(jī)理,結(jié)合監(jiān)測數(shù)據(jù),研究了錨墩在錨索預(yù)應(yīng)力擴(kuò)散中的分擔(dān)作用,并利用Winkler彈性地基梁理論計算了梁體內(nèi)力;黃志毅、馬洪生結(jié)合離心模型試驗,驗證了初參數(shù)法彈性地基梁內(nèi)力求解方法對預(yù)應(yīng)力錨索框架梁的適用性,并分析了框架梁和地基的力學(xué)參數(shù)對梁體內(nèi)力的影響;丁秀美等使用FLAC3D快速拉格朗日差分程序?qū)︻A(yù)應(yīng)力錨索框架作用下坡體附加應(yīng)力分布進(jìn)行了初步研究;呂慶等通過現(xiàn)場試驗,研究了不同預(yù)應(yīng)力錨固結(jié)構(gòu)作用下不同彈性模量的裂隙巖質(zhì)邊坡壓縮變形規(guī)律及附加應(yīng)力,討論了預(yù)應(yīng)力錨索間距和框架梁高度對邊坡表面壓縮效應(yīng)的影響。調(diào)研結(jié)果表明山區(qū)開挖路塹時難以保證邊坡坡面的完全平整,致使預(yù)應(yīng)力錨索框架梁與坡面之間常存在接觸不密實(shí)的情況,不僅會導(dǎo)致框架梁受力狀況改變,還會影響框架梁本身的健康狀況。而目前對這一問題鮮有研究。為此,該文采用FLAC3D建立多工況數(shù)值模型,研究框架梁懸空對預(yù)應(yīng)力錨索框架內(nèi)力的影響。
邊坡防護(hù)工程中,由于施工工藝欠佳或坡體表面不平整,預(yù)應(yīng)力錨索框架與坡體往往不能完全接觸。同時由于巖土體是散體,在預(yù)應(yīng)力錨索框架營運(yùn)階段長時間經(jīng)歷雨水沖刷后,梁體底部的松散土體可能會被掏空,造成框架梁體懸空,減小了梁體與坡體的受力面,改變了梁體的支撐條件,使梁體內(nèi)力增大導(dǎo)致梁體破壞,也可能導(dǎo)致梁體整體下滑。圖1為貴州省某地區(qū)框架梁底破壞后梁體下滑情況。
圖1 框架梁底破壞后梁體下滑
通過實(shí)地調(diào)查,根據(jù)懸空位置不同將其分為懸臂端懸空、錨頭下方懸空和跨中懸空(見圖2)。
圖2 框架梁懸空示意圖
以貴州省常見公路邊坡形式為分析對象,采用FLAC3D建立數(shù)值模型(見圖3),模型尺寸與邊界條件見圖4。垂直于紙面向里為Y軸正方向,Y方向?qū)挾葹?.5 m,框架梁居中布置于坡表面,邊坡坡率為1∶1.0,坡高10 m,邊坡坡后角為零。固定底部X、Y、Z方向,左右側(cè)X方向及前后側(cè)Y方向位移。
圖3 公路邊坡數(shù)值模型
圖4 數(shù)值模型尺寸與邊界條件(單位:m)
為計算簡化,建模時作如下假定:1) 鑒于主要研究對象為均質(zhì)邊坡,巖土體采用Mohr-Coulomb本構(gòu)模型;2) 不考慮錨索框架自身重力及與坡面的摩擦力;3) 巖土體在自重作用下已固結(jié)完成。
參考《工程巖體分級標(biāo)準(zhǔn)》中Ⅳ、Ⅴ級巖體的力學(xué)參數(shù)范圍,設(shè)定巖土體密度為2 250 kg/m3,強(qiáng)度參數(shù)見表1。
表1 巖土體的強(qiáng)度參數(shù)
模擬錨索時,錨索單元的內(nèi)錨固段長取8m,自由段長取15 m,錨固角(錨索與水平面的交角)取25°,錨固段砂漿外圈直徑取13 cm。6股錨索規(guī)格相同,根據(jù)GB/T 5224-2014《預(yù)應(yīng)力混凝土用鋼絞線》,選用規(guī)格為:高強(qiáng)度低松弛標(biāo)準(zhǔn)型鋼鉸線,記為UPS15.20-1860,為由7根鋼絲捻制的標(biāo)準(zhǔn)型鋼絞線,強(qiáng)度級別為1 860 MPa,公稱直徑Dm=15.20 mm,截面參考面積Sm=140 mm2。根據(jù)規(guī)范設(shè)定的錨索單元參數(shù)見表2。
表2 錨索單元參數(shù)
采用內(nèi)置的梁單元模擬框架梁,模型中框架梁沿坡面放置,框架梁尺寸及布置見圖5。梁截面尺寸為50 cm×50 cm,采用對稱配筋。
圖5 框架梁尺寸(單位:cm)
梁單元的參數(shù)參考C30砼,考慮配筋對砼彈性模量的提升可達(dá)10%左右,鋼筋砼彈性模量取32.4 GPa??蚣芰毫簡卧獏?shù)見表3。
表3 框架梁梁單元參數(shù)
FLAC3D中,結(jié)構(gòu)單元之間、結(jié)構(gòu)單元與巖土體單元之間的相互作用通過設(shè)置節(jié)點(diǎn)連接來實(shí)現(xiàn)。在預(yù)應(yīng)力錨索框架梁支護(hù)結(jié)構(gòu)模擬中,將梁節(jié)點(diǎn)與坡面剛接實(shí)現(xiàn)二者同步變位,將錨索索頭節(jié)點(diǎn)與相應(yīng)梁單元節(jié)點(diǎn)剛性連接實(shí)現(xiàn)錨頭和梁體的協(xié)調(diào)變形。圖6為預(yù)應(yīng)力錨索框架加固邊坡的坡面正視圖,其中黑色圓點(diǎn)為梁單元的節(jié)點(diǎn)。坡面寬8.5 m,進(jìn)行17等分,坡面高14.14 m,進(jìn)行29等分,則每個巖土體網(wǎng)格近似為0.5 m×0.5 m正方形??蚣芰嚎v梁長14 m,等分為28個子單元,橫梁長8 m,等分為16個子單元,每個節(jié)點(diǎn)均位于巖土體網(wǎng)格中心。錨索鎖定后的預(yù)應(yīng)力值取500 kN。
圖6 框架梁與坡體耦合關(guān)系示意圖
通過刪除與梁單元節(jié)點(diǎn)連接的巖土體單元模擬對應(yīng)的梁單元節(jié)點(diǎn)脫離巖土體,針對縱梁懸臂端、橫梁懸臂端、錨頭下方、縱梁跨中和橫梁跨中懸空工況研究框架梁懸空對預(yù)應(yīng)力錨索框架內(nèi)力的影響。
各梁體懸空工況下框架梁內(nèi)力的求解過程如下:1) 建立預(yù)應(yīng)力錨索框架加固邊坡模型,預(yù)設(shè)錨索單元的軸向力,此時整個力學(xué)模型并未達(dá)到平衡狀態(tài),進(jìn)行求解,得到初始無懸空狀態(tài)下預(yù)應(yīng)力錨索框架的內(nèi)力。2) 保持其他參數(shù)不變,刪除各工況對應(yīng)的巖土體單元,重新進(jìn)行求解,得到懸空工況下預(yù)應(yīng)力錨索框架的內(nèi)力。
因為框架梁懸空的位置不同,對縱梁和橫梁內(nèi)力的影響規(guī)律不一樣,對縱梁和橫梁彎矩分別進(jìn)行監(jiān)測,定義右側(cè)縱梁為1#梁,下部橫梁為2#梁,中部橫梁為3#梁(見圖7)。
圖7 梁體的編號
(1) 錨頭下方懸空工況。6個錨頭處懸空狀況相同,以錨頭為圓心,半徑0.5、1.0、1.5、2.0 m內(nèi)梁體懸空,對應(yīng)模擬工況見圖8。由于篇幅有限,僅給出部分懸空半徑示意圖。沿半徑方向每刪除一個巖土體單元,懸空半徑增加0.5 m。
圖8 錨頭下方懸空工況數(shù)值模型
(2) 縱梁懸臂端懸空工況。4個縱梁懸臂端懸空狀況相同,以縱梁懸臂端端點(diǎn)為圓心,半徑0.5、1.0、1.5、2.0 m內(nèi)梁體懸空,部分模型見圖9。
圖9 縱梁懸臂端懸空工況數(shù)值模型
(3) 橫梁懸臂端懸空工況。6個橫梁懸臂端懸空狀況相同,以橫梁懸臂端端點(diǎn)為圓心,半徑0.5、1.0、1.5、2.0 m內(nèi)梁體懸空,部分模型見圖10。
圖10 橫梁懸臂端懸空工況數(shù)值模型
(4) 縱梁跨中懸空工況。4個縱梁跨中懸空狀況相同,以縱梁跨中為圓心,半徑0.5、1.0、1.5、2.0、2.5 m內(nèi)梁體懸空,部分模型見圖11。
圖11 縱梁跨中懸空工況數(shù)值模型
(5) 橫梁跨中懸空工況。3個橫梁跨中懸空狀況相同,以橫梁跨中為圓心,半徑0.5、1.0、1.5、2.0 m內(nèi)梁體懸空,部分模型見圖12。
圖12 橫梁跨中懸空工況數(shù)值模型
不同錨頭下方懸空半徑下梁體內(nèi)力見圖13,定義框架梁靠近坡側(cè)受拉彎矩為正。由圖13可知:錨頭下方懸空,縱梁錨頭下方截面、跨中截面、橫梁錨頭下方截面、跨中截面彎矩絕對值均增大。
圖13 不同錨頭下方懸空半徑下各梁彎矩
不同縱梁懸臂端懸空半徑下梁體內(nèi)力見圖14。由圖14可知:縱梁懸臂端懸空,縱梁兩端負(fù)彎矩變大,正彎矩變小,中部彎矩變化不明顯;下部橫梁錨頭下方和跨中彎矩絕對值增大;中部橫梁彎矩變化不明顯。
圖14 不同縱梁懸臂端懸空半徑下各梁彎矩
不同橫梁懸臂端懸空半徑下梁體內(nèi)力見圖15。
圖15 不同橫梁懸臂端懸空半徑下各梁彎矩
由圖15可知:橫梁懸臂端懸空,縱梁正、負(fù)彎矩均變大;橫梁錨頭下方彎矩絕對值變小,跨中彎矩絕對值變大。懸臂端懸空在忽略梁體自重的情況下可等效為懸臂梁長度變短,變化規(guī)律與使用Winkler彈性地基梁模型的研究結(jié)論一致,說明通過刪除與梁單元節(jié)點(diǎn)連接的巖土體單元模擬對應(yīng)的梁體脫離土體懸空的方法可行。
不同縱梁跨中懸空半徑下梁體彎矩見圖16。由圖16可知:縱梁跨中懸空,縱梁正、負(fù)彎矩絕對值均減小,且中間錨頭下方與跨中彎矩有變號趨勢;橫梁正、負(fù)彎矩絕對值均增大。
圖16 不同縱梁跨中懸空半徑下各梁彎矩
橫梁跨中懸空,彎矩變化與縱梁跨中懸空相似,橫梁正、負(fù)彎矩絕對值均減小,且跨中彎矩有變號趨勢;縱梁正、負(fù)彎矩絕對值均增大。因篇幅限制,圖表不再列出。
定義使任意梁形成簡支約束或使整個懸臂段懸空的工況懸空半徑為這類工況的最大懸空半徑,則錨頭下方懸空、橫梁跨中懸空、橫梁懸臂端懸空和縱梁懸臂端懸空所能達(dá)到的最大懸空半徑為2 m,縱梁跨中懸空所能達(dá)到的最大懸空半徑為2.5 m。定義i工況的實(shí)際懸空半徑ri與這類工況所能達(dá)到的最大懸空半徑rmax的比例為懸空比rr,即
rr=ri/rmax×100%
(1)
定義各懸空工況下框架梁最大彎矩絕對值Mi與無懸空工況下框架梁最大彎矩Mmax的比值為最大彎矩相對比例Mr,即
Mr=Mi/Mmax×100%
(2)
各工況下最大彎矩相對比例見圖17。由圖17可知:錨頭下方懸空造成的框架最大彎矩絕對值增長程度比懸臂端和跨中懸空時大,對梁體內(nèi)力的影響大,為最不利懸空位置。
基于FLAC3D建立數(shù)值模型,通過刪除與梁單元節(jié)點(diǎn)連接的巖土體單元模擬框架梁懸空,研究各部位懸空對框架梁內(nèi)力的影響。結(jié)論如下:
圖17 各工況下最大彎矩相對比例
(1) 錨頭下方懸空導(dǎo)致所有梁體的彎矩絕對值都增大,懸臂端和跨中懸空導(dǎo)致部分梁體彎矩絕對值增大、部分減小。
(2) 相對于懸臂和跨中懸空,錨頭下方懸空造成的框架最大彎矩絕對值增長程度較大,錨頭下方懸空對梁體內(nèi)力的影響大于懸臂懸空和跨中懸空,對梁體內(nèi)力的影響大,為最不利懸空位置。