武博強(qiáng),劉 青,楊德宏,王 勇,蘇衛(wèi)衛(wèi)
(1.中交第一公路勘察設(shè)計研究院有限公司,陜西 西安 710065;2. 中鐵第一勘察設(shè)計院集團(tuán)有限公司,陜西 西安 710043)
公路土質(zhì)路塹邊坡通常根據(jù)規(guī)范中規(guī)定的坡率法進(jìn)行相應(yīng)的路基邊坡設(shè)計,但有些項目路塹高邊坡由于占地問題而導(dǎo)致區(qū)域放坡受限,即使一級邊坡采用擋土墻收坡后仍不能滿足占地要求,只能采用陡坡率,從而不能保證路塹邊坡的穩(wěn)定性,因此需要進(jìn)行邊坡預(yù)加固。邊坡預(yù)加固方法中,注漿預(yù)加固可以保證土質(zhì)邊坡在較陡設(shè)計坡率下仍可以滿足工程需要,是解決公路路塹邊坡放坡受限的有效方法之一,但國內(nèi)對預(yù)加固后土質(zhì)邊坡穩(wěn)定性研究較少,分析方法處于摸索階段,至今沒有得到統(tǒng)一,目前國內(nèi)外研究學(xué)者普遍采用的方法皆可歸類為抗剪強(qiáng)度加固區(qū)域置換法。史秀志等[1]基于Flac3D建立均質(zhì)假三維邊坡模型討論了在不同坡角下黏聚力和內(nèi)摩擦角對均質(zhì)邊坡穩(wěn)定性的影響程度,但實際工程中不可能將水泥土的抗剪強(qiáng)度應(yīng)用于整個邊坡模型,且水泥土樁的抗剪強(qiáng)度與樁間土體抗剪強(qiáng)度的置換關(guān)系沒有明確。蔡先慶等[2]同樣基于Flac3D建立模型對邊坡抗滑樁前土體進(jìn)行注漿加固,保證抗滑樁樁頂位移量滿足要求,但注漿區(qū)域的選取與加固區(qū)域土體綜合參數(shù)的選取方法均未說明來由。根據(jù)抗剪強(qiáng)度土體置換原理,該方法其實并不適用于有限元及拉格朗日元模型,加固土體與原土體的抗剪強(qiáng)度、剪切模量等參數(shù)相差較大,因此該方法是一種近似方法,但此方法獲得的穩(wěn)定系數(shù)的準(zhǔn)確度并沒有學(xué)者進(jìn)行研究。為探究該方法原理及準(zhǔn)確度,以內(nèi)蒙古經(jīng)烏高速公路風(fēng)積沙路塹邊坡為例[3],基于實際注漿試驗數(shù)據(jù)分析了加固土體與原土體抗剪強(qiáng)度的關(guān)系,給出了抗剪強(qiáng)度加固區(qū)域置換法的正確分析步驟,總結(jié)了不同加固區(qū)域與邊坡穩(wěn)定性的關(guān)系,討論了高邊坡穩(wěn)定系數(shù)產(chǎn)生偏差的原因,并在此方法的基礎(chǔ)上提出了更貼近實際的注漿樁體模型分析法,基于現(xiàn)場注漿試驗數(shù)據(jù)建立實體邊坡加固模型,討論兩種方法的優(yōu)缺點,旨在為國內(nèi)外研究人員和設(shè)計者在公路土質(zhì)邊坡預(yù)加固穩(wěn)定性分析方面提供一些新的啟示。
抗剪強(qiáng)度加固區(qū)域置換法的原理是將加固土樁的抗剪強(qiáng)度通過樁截面面積與樁間加固面積的關(guān)系置換為加固區(qū)域土體的綜合抗剪強(qiáng)度,采用該抗剪強(qiáng)度進(jìn)行邊坡穩(wěn)定性分析,相關(guān)公式如下[4]:
τps=nτp+(1-n)τs,
(1)
cps=ncp+(1-n)cs,
(2)
tanφps=ntanφp+(1-n)tanφs,
(3)
式中,τps為加固土體抗剪強(qiáng)度;τs為土體抗剪強(qiáng)度;τp為樁體抗剪強(qiáng)度;n為置換率;cps為加固土體黏聚力;cp為樁體黏聚力;cs為土體黏聚力;φps為加固土體內(nèi)摩擦角;φp為樁體內(nèi)摩擦角;φs為土體內(nèi)摩擦角。
根據(jù)上述公式(1)~(3),加固土體的抗剪強(qiáng)度和土體的抗剪強(qiáng)度關(guān)系可簡化為過原點且斜率為k的直線公式(4)~(5),根據(jù)前人的研究成果[5-6]及內(nèi)蒙古經(jīng)烏高速公路項目的現(xiàn)場注漿試驗情況,注漿后土顆粒間隙被一定的注漿壓力沖破,形成多個貫通的注漿通道。顆粒間空隙被水泥漿液填充,土體間聯(lián)結(jié)增強(qiáng)。注漿有效半徑之外的土體由于漿液的壓密作用孔隙率也有一定程度的減小,顆粒從“指尖對指尖”變化為“指尖對指縫”。由于顆粒間的嵌擠作用,黏聚力增加明顯,從公式(4)~(5)可以看出,Kc的提升幅度遠(yuǎn)大于Kφ。
cps=KCcs,
(4)
tanφps=Kφtanφs,
(5)
式中,Kc為黏聚力提高系數(shù);Kφ為內(nèi)摩擦角正切值提高系數(shù)。
抗剪強(qiáng)度加固區(qū)域置換法原理簡單明了,但該加固區(qū)域土體采用綜合抗剪強(qiáng)度參數(shù)本身就是一種近似的方法,密度、泊松比等根據(jù)土體自身確定的參數(shù)均不應(yīng)通過面積置換綜合取值,因此更適用于指標(biāo)單一的極限平衡法,且方法的應(yīng)用準(zhǔn)確度未知。
風(fēng)積沙僅存在表觀黏聚力,無真實黏聚力,由于占地緊張,路段采用1∶1,1∶1.25,1∶1.25坡率8 m 分級逐級放坡,并設(shè)置2 m寬度臺階,風(fēng)積沙路塹邊坡計算模型如圖1所示,相關(guān)邊坡巖土體參數(shù)見表1。
圖1 路塹邊坡計算模型Fig.1 Calculation model of cut slope
表1 風(fēng)積沙土體計算參數(shù)Tab.1 Calculation parameters of aeolian sandy soil
研究區(qū)路塹邊坡屬于地質(zhì)條件簡單且破壞后危害較輕的邊坡,風(fēng)積沙土體特性為非黏性土,其可能發(fā)生的滑動區(qū)域為折線形區(qū)域,采用簡布法進(jìn)行計算,根據(jù)可能的滑動區(qū)域確定注漿預(yù)加固區(qū)域后,邊坡土體組成發(fā)生了改變,預(yù)加固后邊坡加固區(qū)為黏性土體,未加固區(qū)為非黏性土體,應(yīng)按照沿加固區(qū)內(nèi)部以及交界面圓弧形滑動和沿加固區(qū)外折線形滑動兩種情況分別進(jìn)行分析,圓弧形滑動采用Bishop法,折線形滑動采用隱式解法進(jìn)行分析,并采用簡布法進(jìn)行校核[7-12]。
路塹邊坡穩(wěn)定性分析應(yīng)對可能發(fā)生的局部滑動及整體滑動進(jìn)行逐級分析,邊坡模型穩(wěn)定性分析結(jié)果由表2可知,一級坡單獨發(fā)生局部滑動的可能性最大,一級坡牽引二級坡發(fā)生局部滑動的可能性次之,整體滑動時的邊坡穩(wěn)定系數(shù)為Fs,屬于不穩(wěn)定狀態(tài),二級坡和三級坡發(fā)生局部滑動的可能性相對較小,屬于欠穩(wěn)定狀態(tài)。因此應(yīng)根據(jù)一級坡、一級坡牽引二級坡滑動以及邊坡整體滑動的潛在滑動面綜合確定注漿區(qū)域與注漿深度。
表2 邊坡模型穩(wěn)定性分析結(jié)果Tab.2 Analysis result of stability of slope model
在注漿試驗中,注漿有效加固半徑為0.25 m,樁間距為1.5 m,注漿提升速度50 cm/min,噴漿泵量35~50 L/min,正方形分布,樁體黏聚力為1.1 MPa,內(nèi)摩擦角為42°,據(jù)此得到加固區(qū)域置換率為0.087 2,采用抗剪強(qiáng)度加固區(qū)域置換法后的綜合黏聚力為96.8 kPa,綜合內(nèi)摩擦角為34.8°。為保證施工質(zhì)量,注漿深度不小于潛在滑動面以下0.5 m,一級坡局部滑動、一級坡牽引二級坡滑動以及邊坡整體滑動的注漿區(qū)域如圖2~圖4所示。
圖2 一級坡局部滑動注漿區(qū)域Fig.2 Grouting area of local sliding of first grade slope
圖3 一級坡牽引二級坡滑動注漿區(qū)域Fig.3 Grouting area of sliding of secondary slope draged by primary slope
圖4 整體滑動注漿區(qū)域Fig.4 Grouting area of integral sliding
確定3種注漿區(qū)域后,假設(shè)邊坡潛在最危險滑動面位于加固區(qū)內(nèi),將發(fā)生圓弧形滑動破壞,但計算結(jié)果顯示潛在最危險圓弧滑動面均位于非加固區(qū)內(nèi),一級坡單獨預(yù)加固時最危險圓弧滑動面位于三級坡區(qū)域,一級坡與二級坡均預(yù)加固時最危險圓弧滑動面與一級坡單獨預(yù)加固時一致,仍位于三級坡區(qū)域(圖5),邊坡整體預(yù)加固時最危險圓弧滑動面位于加固區(qū)域以外(圖6),這說明由于注漿區(qū)域深度較淺,加固土體由于黏聚力的大幅度增加使得抗滑力也大幅度增加,在加固區(qū)域中產(chǎn)生圓弧滑動的可能性小,應(yīng)以加固區(qū)域外產(chǎn)生的折線形滑動破壞分析為準(zhǔn)。
圖5 一、二級邊坡預(yù)加固后潛在最危險圓弧滑動面Fig.5 Potential most dangerous arc sliding surface after pre-reinforcement of primary and secondary slopes
圖6 邊坡整體預(yù)加固后潛在最危險圓弧滑動面Fig.6 Potential most dangerous arc sliding surface after integral pre-reinforcement of slope
折線形滑動穩(wěn)定性分析結(jié)果由表3可知,一級邊坡與二級邊坡預(yù)加固時的最危險滑動位置與僅一級邊坡預(yù)加固時的滑動位置相同,均為三級坡局部折線形滑動?;瑒游恢门c圖5中的圓弧滑動面位置相近,最小穩(wěn)定系數(shù)均為1.01,與未加固時的三級坡局部滑動穩(wěn)定系數(shù)相同。這說明未加固時局部滑動的穩(wěn)定系數(shù)雖然最小,但針對局部區(qū)域預(yù)加固后僅對加固區(qū)域有效果,先前安全系數(shù)較高的區(qū)域反而發(fā)生局部滑動,因此應(yīng)針對整體滑動面進(jìn)行預(yù)加固。邊坡整體預(yù)加固后的最危險滑動面位置相較圖6中的圓弧滑動面更貼近坡面(圖7),為整體式淺層折線滑動,最小穩(wěn)定系數(shù)為1.185,屬于基本穩(wěn)定狀態(tài),但未滿足設(shè)計規(guī)范的要求,且在整體邊坡加固分析中,各級邊坡產(chǎn)生局部滑動的穩(wěn)定系數(shù)均較大。這說明由于加固區(qū)域的抗剪強(qiáng)度大,局部剪出的可能性均較小。
表3 邊坡預(yù)加固模型穩(wěn)定性分析結(jié)果Tab.3 Stability analysis result of slope pre-reinforcement model
圖7 邊坡整體預(yù)加固后潛在最危險滑動面對比Fig.7 Comparison of potential most dangerous sliding faces after integral pre-reinforcement of slope
由于水泥土樁的抗剪強(qiáng)度是風(fēng)積沙土體抗剪強(qiáng)度的1 100倍,在樁間距處于一個區(qū)間的時候,即使樁間距接近該區(qū)間的上限,樁間土體由于擠土作用形成土拱,邊坡發(fā)生滑動的可能性仍較小,但樁間距的擴(kuò)大使得強(qiáng)度置換法的置換率將變小,綜合抗剪強(qiáng)度的黏聚力和內(nèi)摩擦角均有不同程度的下降。所以根據(jù)強(qiáng)度置換法的原理及計算結(jié)果,在注漿參數(shù)確定的前提下,該邊坡注漿后穩(wěn)定系數(shù)不會比置換法得到的穩(wěn)定系數(shù)小。置換法是一種相對保守的方法,但破壞方式及最危險滑動面類似,均為注漿區(qū)域以外的折線形滑動破壞?;诖?,筆者提出一個更貼近實際情況的方法進(jìn)行對比論證,即建立注漿樁體模型法。采用柱體單元建立邊坡加固模型,基于強(qiáng)度折減法及位移剪應(yīng)變增量等多方面分析邊坡注漿預(yù)加固后的實際穩(wěn)定性。注漿樁體模型法所依據(jù)的注漿區(qū)域及土體注漿等參數(shù)完全同上文,加固土體與土體接觸面參數(shù)根據(jù)式(6)取值[13](表4)。
表4 注漿樁體參數(shù)Tab.4 Parameters of grouting pile
(6)
式中,K為體積模量;G為剪切模量;Δzmin為接觸面法向連接區(qū)域上的最小尺寸,根據(jù)模型網(wǎng)格尺寸確定。(該公式引自文獻(xiàn)[13])
據(jù)此建立注漿預(yù)加固模型,注漿樁體為圓柱實體模型(圖8),考慮國內(nèi)目前普遍采用平臺注漿法,注漿設(shè)計頂面應(yīng)位于設(shè)計坡面線以上0.5 m,便于刷坡時鑿除停漿后浮漿部分。采用Rhino軟件建立預(yù)加固邊坡模型后通過Griddle插件輸出至Flac3D軟件進(jìn)行分析,由于注漿樁體呈正方形分布,樁間距為1.5 m,坡長方向均呈規(guī)律性分布,為減少計算時長,模型y方向(坡長方向)延伸長度為1.5 m,即樁體直徑0.5 m,兩側(cè)各0.5 m土體。具體建立步驟為:(1)在邊坡注漿區(qū)域建立樁體模型;(2)邊坡二維平面擠出邊坡三維模型;(3)將樁體與邊坡重合部分進(jìn)行布爾差集運算,去除重復(fù)部分;(4)生成網(wǎng)格模型(圖9)。
圖8 加固土樁實體網(wǎng)格模型Fig.8 Solid grid model of reinforced soil pile
圖9 預(yù)加固邊坡實體網(wǎng)格模型Fig.9 Solid grid model of pre-reinforced slope
計算分析得到基于強(qiáng)度折減法[14-20]的邊坡穩(wěn)定系數(shù)Fs=1.221>1.20(圖10),屬于穩(wěn)定狀態(tài),滿足設(shè)計規(guī)范的要求,相較強(qiáng)度置換法穩(wěn)定系數(shù)增大30.4%。根據(jù)最大剪應(yīng)變增量圖顯示(圖11),最危險滑動面與強(qiáng)度置換法基本相同,破壞時邊坡的x向最大位移為6.68 m,破壞方式是加固區(qū)域以下折線型滑動,與前文的宏觀判斷基本一致。所以強(qiáng)度置換法按照潛在滑裂面位置確定注漿深度是正確的,但穩(wěn)定性分析結(jié)果偏小,且較注漿樁體模型法誤差較大。兩種方法的對比如表5所示。
圖10 預(yù)加固邊坡x向位移云圖(單位:m)Fig.10 Nephogram of pre-reinforced slope in x direction(unit:m)
圖11 預(yù)加固邊坡最大剪應(yīng)變增量云圖Fig.11 Nephogram of maximum shear strain increment of pre-reinforced slope
表5 預(yù)加固方法優(yōu)缺點Tab.5 Advantages and disadvantages of pre-reinforcement method
(1)通過實例分析,抗剪強(qiáng)度加固區(qū)域置換法抗剪強(qiáng)度主要提升在黏聚力,內(nèi)摩擦角變化較小,且該方法更適用于極限平衡法,從置換法的原理得出的邊坡穩(wěn)定性結(jié)果偏小,導(dǎo)致在后續(xù)設(shè)計中注漿深度普遍較大,實際穩(wěn)定性結(jié)果較為保守。
(2)通過對深挖路塹邊坡可能的局部滑動以整體滑動發(fā)現(xiàn),雖然局部滑動的可能性最大,但確定注漿區(qū)域應(yīng)通過整體潛在滑動面確定,對于非黏性土,加固前后的滑動破壞模式均為折線形滑動。
(3)注漿樁體模型法很好地彌補(bǔ)了抗剪強(qiáng)度加固區(qū)域置換法在計算準(zhǔn)確度上的不足,采用與實際相符的實體建模,真實地反映了邊坡預(yù)加固后的穩(wěn)定性。邊坡潛在破壞面與強(qiáng)度置換法在滑動面位置與形狀的判斷結(jié)果上相似,計算準(zhǔn)確度較強(qiáng)度置換法提升了30.4%。注漿樁體模型法在使用前應(yīng)進(jìn)行注漿試驗確定準(zhǔn)確的注漿參數(shù)并通過加固體試樣相關(guān)試驗確定接觸面等參數(shù)。此方法不僅適用于風(fēng)積沙等非黏性土體邊坡,也適用于黏性土邊坡的加固,可根據(jù)不同土體的種類、孔隙率、含水量等調(diào)整加固樁體的類型與注漿參數(shù)。
(4)根據(jù)抗剪強(qiáng)度加固區(qū)域置換法與注漿樁體模型法的優(yōu)缺點比較,推薦在工可及初步設(shè)計階段采用抗剪強(qiáng)度加固區(qū)域置換法進(jìn)行邊坡預(yù)加固穩(wěn)定性分析,便于快速定量分析設(shè)計,保證工程量預(yù)算,在施工圖設(shè)計階段采用注漿樁體模型法進(jìn)行精細(xì)化設(shè)計。