拉筋帶對地震作用下土工格室擋墻力學性能的影響

靳飛飛,宋 飛,石 磊,朱 婕

(長安大學 公路學院,西安 710064)

中國地震頻發,地震觸發的邊坡滑坡造成了大量的人員傷亡和經濟損失[1]。土工格室擋墻結構因其結構巧、施工易、造價低等優勢,在公路、鐵路、機場、電力和冶金礦山等工程領域中得到了廣泛的應用[2]。通過土工格室擋墻對邊坡進行支護其抗震性能得到顯著提高[3-4],然而,目前關于土工格室擋墻支護邊坡的動力問題研究較少,地震作用下土工格室擋墻支護邊坡的穩定性研究遠落后于工程實踐,因此對其進行研究十分必要且迫切。

目前,前人對靜力條件下土工格室擋墻支護邊坡的穩定性研究已經取得豐碩的成果[5-7]。文獻[8]采用梁模型來模擬土工格室的柔性板基礎,對加固邊坡的土工格室層不同位置,不同層數以及加固層之間的垂直間距、長度、厚度、土工格室的楊氏模量進行了研究;文獻[9]通過室內模型試驗研究了土工格室參數對擋墻支護邊坡穩定性的影響;文獻[10]采用有限元軟件對不同高寬比和不同土體強度下的土工格室擋墻的破壞模式進行了研究,給出了不同工況下土工格室擋墻破壞模式的數學表達式;文獻[11]基于極限平衡法,提出了一種分析土工格室加筋邊坡穩定性的方法;文獻[12]采用極限平衡分析方法,提出了路堤臨界高度的計算模型,并討論了格室參數對路堤臨界高度的影響規律。

迄今為止,人們對于地震荷載作用下土工格室擋墻支護邊坡的穩定性研究卻鮮有涉及。鑒于此,本文利用FLAC3D數值軟件建立土工格室擋墻支護邊坡數值模型[13-14],系統分析了不同位置和不同數量的拉筋帶對邊坡動力響應的影響,以期為土工格室擋墻結構的抗震設計提供可靠的依據。

1 格室約束圍壓理論解

土工格室內部充填料在外荷載σ1作用下將產生軸向壓縮和側向膨脹變形,引起格室產生拉伸應變和拉應力,而格室的拉應力對于格室內充填料產生側向約束力,如圖1所示,這種側向約束力稱為約束圍壓[15]。

格室對于充填料的約束圍壓可按照式(1)進行計算[15],即

(1)

(2)

(3)

式中:σg為格室約束圍壓;Mt為格室條帶剛度;εc為格室周向應變;Dε為格室在周向應變εc時等效直徑;D0為格室初始等效直徑;ε1為格室軸向應變。

2 邊坡模型建立

2.1 數值模型建立

參考文獻[13-14]建立土工格室擋墻支護邊坡數值模型,如圖2所示。該模型由基底和邊坡兩部分組成,基底長×寬×高=7.0 m×0.4 m×0.2 m,邊坡由格室擋墻和墻后坡體兩部分組成,擋墻由14層土工格室及格室內充填料組成,格室內充填料為礫石,基底和坡體填料均為砂土,每個格室單元的高度和等效直徑均為20 cm,格室厚度為0.12 cm,邊坡高為2.8 m,厚為0.4 m,坡比為1∶0.5。由于FALC3D中前處理建模具有一定的局限性,為了體現土工格室單元的曲面效應,本文利用Rhino軟件建立土工格室單元,應用Rhino To Flac3D接口程序將土工格室單元導入到FALC3D中進行動力計算。

圖2 邊坡數值模型(m)

邊坡數值模型共1 338個單元,2 235個節點,土工格室與坡體、填料之間均設置接觸面,該接觸面服從Mohr-Coulomb準則。文獻[13]通過對土工格室條帶和節點進行拉伸試驗,測得節點強度和條帶屈服強度接近?；诒疚臄抵的M提取的格室應變計算得到格室拉力均小于格室條帶屈服強度和節點強度。本文數值模擬中礫石和砂土采用彈塑性本構關系,屈服準則均采用適用于巖土體力學行為的Mohr-Coulomb準則[16-20],土工格室采用彈塑性本構關系。模型中材料物理力學參數[13-14]如表1、2所示。

表2 接觸面力學參數

2.2 邊界條件、阻尼和地震波選取

為了避免地震波傳播時產生的反射及能量耗散,本次數值模型的四周及底部邊界均設置為自由場邊界。阻尼選取局部阻尼,阻尼系數為0.158。

選取日本神戶海洋氣象臺觀測點記錄的近場脈沖地震波Kobe波作為輸入地震波[13-14],XZ雙向耦合施加于模型底部,Z向加速度為X向加速度的0.38倍,其中X向地震波加速度峰值為1.193g,卓越頻率為1.9 Hz,主頻為1～3.8 Hz;Z向地震波加速度峰值為0.463g,卓越頻率為1.18 Hz,主頻為0.9～5 Hz,持時為20 s,如圖3所示。

圖3 Kobe波加速度時程曲線

2.3 監測點布置

為了便于分析,對不同位置的土工格室進行定義,將臨空面處的格室定義為墻面格室,將墻體內的格室定義為墻體格室。

在墻面自坡腳至坡頂布置14個位移監測點(編號依次為D1～D14);在墻面格室自坡腳至坡頂布置14個加速度監測點(編號依次為A1～A14);在墻面格室的內壁布置應變監測點,每個格室單元前后對稱位置各布設一個監測點,計算時取二者的平均值作為最終應變值(編號依次為Y1～Y14,Y1′～Y14′);在坡頂自坡肩向坡體深處共布置10個沉降監測點(編號依次為C1～C10),各監測點布置如圖4所示。

圖4 監測點布置

3 計算結果分析

不同位置的拉筋帶,自墻面向坡體方向,均為9個格室單元長度(包括墻體3個格室單元)。為了研究拉筋帶的力學特性,在9個格室單元內均布置應變監測點和加速度監測點,其中應變監測點在每個格室單元前后對稱位置各布設一個,計算時取二者的平均值作為最終應變值,拉筋帶及其監測點布置如圖5所示。

圖5 拉筋帶及其監測點布置示意圖

3.1 數值模擬結果驗證

為了驗證數值模型的準確性,建立與文獻[13-14]中的5個試驗模型對應的數值模型,數值模型中的參數設置、試驗條件、試驗要素、監測點布設、監測指標以及地震波選取均按照試驗[13-14]提供的設置。圖6和圖7分別為第1次和第3次振動的數值計算結果與振動臺試驗結果對比,通過對比可以看出,第1次振動和第3次振動的數值模擬計算結果與振動臺試驗結果及變化規律相似,說明數值計算與試驗結果相吻合,數值模擬計算結果是準確可靠的。

圖6 第1次振動數值模擬與試驗結果對比

圖7 第3次振動數值模擬與試驗結果對比

3.2 拉筋帶位置對邊坡動力響應的影響

根據文獻[21]的研究,本文選擇0、H/3、H/2、3H/4、H(H為邊坡高度)這5個位置作為邊坡拉筋帶的布設位置。以圖2所示模型為基本模型,首先布設兩層拉筋帶,通過計算確定兩層拉筋帶布設最優位置;基于兩層拉筋帶最優位置,進行3層拉筋帶的布設,最后,進行4層拉筋帶的布設。

3.2.1 拉筋帶位置對格室約束圍壓的影響

圖8(a)為布設兩層拉筋帶墻面格室約束圍壓沿高程分布。由圖可知,當在(0,H/3)、(0,H/2)和(0,3H/4)布設拉筋帶時,格室約束圍壓的最大值依次出現在H/3、H/2和3H/4處,說明在坡體中布設拉筋帶對該位置附近的坡體起到了有效的錨固作用,而在(0,H)布設拉筋帶時,格室約束圍壓的最大值出現在約H/3處,這是由于土壓力合力作用點位置大約在距離墻底約H/3處,且坡頂處拉筋帶緊靠坡頂臨空面,在地震作用下,格室內填料將向臨空面運動,格室約束效果減弱;當在(H/3,H/2)和(H/3,3H/4)布設拉筋帶時,二者中下部的格室約束圍壓值基本接近,中上部的格室約束圍壓差值增大,且在(H/3,H)布設時,格室約束圍壓達到所有布設方案的最大值,說明在(H/3,H)處布設拉筋帶對于墻體的整體穩定最為有利;當在(H/2,3H/4)布設拉筋帶時,與(H/3,3H/4)相比,中上部格室約束圍壓二者基本接近,但中下部前者遠小于后者,說明在H/3處布設拉筋帶對于提高擋墻穩定具有重要作用;當在(3H/4,H)布設拉筋帶時,與(0,H)相比,坡腳附近格室約束圍壓前者小于后者,中上部卻相反。

圖8(b)為布設3層拉筋帶時墻面格室約束圍壓沿高程分布。由圖可知,當在(0,H/3,H)布設拉筋帶時格室約束圍壓最大,最大值相較于(H/3,H),格室約束圍壓由97.61 kPa增加至105.71 kPa,增幅達到8.29%,說明在坡腳處布設拉筋帶對于提高墻體穩定同樣具有重要作用;當在(H/3,H/2,H)和(H/3,3H/4,H)布設拉筋帶時,中下部格室約束圍壓前者大于后者,而中上部卻相反?？梢?布設拉筋帶雖然在一定程度上能夠提高邊坡穩定,但其布設位置不合理時將不能充分發揮格室力學性能,造成材料浪費。此外,3種方案下格室約束圍壓的最大值均出現在H/3處,說明在H/3處布設拉筋帶是必要且合理的。

圖8(c)為布設4層拉筋帶時墻面格室約束圍壓沿高程分布。由圖可知,在(0,H/3,3H/4,H)布設拉筋帶時,格室約束圍壓大于(0,H/3,3H/4,H),對應的最大值為111.39 kPa,相較于(0,H/3,H),增幅達到5.38%,可見,在邊坡中下部和中上部合理布設拉筋帶,能夠顯著增強邊坡抗震性能。

綜上所述,在實際工程中應優先考慮在0,H/3,3H/4和H這4個位置布設拉筋帶。

3.2.2 拉筋帶位置對擋墻位移的影響

圖9(a)為布設兩層拉筋帶時墻面永久水平位移沿高程變化。由圖可知,當在(0,H/3)和(0,H/2)布設拉筋帶時,位移沿高程呈非線性增加,在坡頂處位移達到最大值;而在(0,3H/4)和(0,H)布設拉筋帶時,位移沿高程呈“增加-減小”的兩段形態,最大位移出現在H/3處,說明拉筋帶布設位置上移能夠對坡頂附近的位移起到有效遏制作用;當在(H/3,H/2)布設拉筋帶時,位移最大值出現在坡頂處,而在(H/3,3H/4)和(H/3,H)布設拉筋帶時,位移最大值又出現在H/3處,由此說明,僅在邊坡中下部布設拉筋帶時最大位移出現在坡頂附近,而僅在邊坡中上部布設拉筋帶時最大位移出現在H/3處;當在(H/2,3H/4)和(H/2,H)布設拉筋帶時,邊坡中下部位移值二者基本接近,但中上部位移逐漸減小,而在(3H/4,H)布設拉筋帶時,坡腳附近的位移大于坡頂附近。

圖9(b)為布設3層拉筋帶時墻面永久水平位移沿高程變化。由圖可知,當在(0,H/3,H)布設拉筋帶時,墻面永久水平位移值最小,其位移值相較于(H/3,H),由19.78 cm減小至18.65 cm,減幅達到5.68%;而在(H/3,H/2,H)和(H/3,3H/4,H)布設拉筋帶時,二者中上部位移逐漸減小,中下部位移基本不變。

圖9(c)為布設4層拉筋帶時墻面永久水平位移沿墻高變化。由圖可知,當在(0,H/3,3H/4,H)布設拉筋帶時位移最小,最小值為18.06 cm,相較于(0,H/3,H),減幅達到5.16%;而在(0,H/3,H/2,H)布設拉筋帶時,其位移值小于(0,H/3,3H/4,H),坡腳附近二者位移差值量較小,但在坡頂附近位移差值量較大,可見保持邊坡中下部拉筋帶位置和數量不變的情況下,在3H/4處布設拉筋帶更能有效提高邊坡的整體穩定。

綜上所述,在坡腳和坡頂附近布設拉筋帶更有助于提高邊坡的抗震性,這是由于邊坡破壞時,在坡腳發生剪切破壞,坡頂發生張拉破壞,當在0和H/3布設拉筋帶時可阻止坡腳處剪切破壞的發生,在3H/4和H布設拉筋帶時可阻止坡頂張拉破壞的發生,因此在這4個位置布設拉筋帶是合理的。

3.2.3 拉筋帶位置對擋墻加速度的影響

圖10(a)為布設兩層拉筋帶時墻面水平峰值加速度沿高程變化。由圖可知,峰值加速度沿墻高表現為“增加-減小-增加”的三段形態,在坡頂處峰值加速度達到最大值。具體來看,當在(0,H/3)、(0,H/2)、(0,3H/4)和(0,H)布設拉筋帶時,邊坡中下部峰值加速度逐漸增大,而中上部逐漸減小,說明隨著拉筋帶布設位置的上移,坡頂附近加速度放大效應減弱;當在(H/3,H/2)、(H/3,3H/4)和(H/3,H)布設拉筋帶時,峰值加速度逐漸減小,且在(H/3,H)處布設時峰值加速度達到所有方案的最小值,說明在該處布設拉筋帶加速度放大效應最弱;當在(H/2,3H/4)和(H/2,H)布設拉筋帶時,在邊坡中下部二者峰值加速度接近,中上部逐漸減小,而在(3H/4,H)布設拉筋帶時,坡腳附近峰值加速度大于坡頂附近,說明在坡頂處布設拉筋帶,對于減小坡頂附近的加速度放大效應更為有利。

圖10 拉筋帶位置對加速度的影響

圖10(b)為布設3層拉筋帶時墻面水平峰值加速度沿墻高變化。由圖可知,當在(0,H/3,H)布設拉筋帶時,峰值加速度最小,相較于(H/3,H)峰值加速度的最大值由17.05 m/s2減小至16.07 m/s2,減幅達到5.77%;而在(H/3,H/2,H)和(H/3,3H/4,H)布設拉筋帶時,邊坡中下部二者峰值加速度接近,而上部逐漸減小,說明在保持邊坡中下部布設位置不變的情況下,拉筋帶位置越靠近坡頂,加速度放大效應越小。

圖10(c)為布設4層拉筋帶時墻面水平峰值加速度沿墻高變化。由圖可知,當在(0,H/3,3H/4,H)布設拉筋帶時,峰值加速度最小,峰值加速度最大值為15.53 m/s2,相較于(0,H/3,H),減幅達到3.36%;而在(0,H/3,H/2,H)布設拉筋帶時,在坡腳處二者數值相近,在坡頂處二者差值量增大。

綜上所述,在0、H/3、3H/4和H處布設時,峰值加速度放大效應最小,邊坡最為穩定。

3.2.4 拉筋帶位置對坡頂沉降的影響

圖11為不同拉筋帶布設時坡頂沉降分布。由圖可知,不同布設方案下,坡頂沉降自坡肩向坡體延伸,均表現為先增加后減小的“V”型分布。具體來看,由圖11(a)可知,當在(0,H/3)、(0,H/2)、(0,3H/4)和(0,H)處布設拉筋帶時,坡頂沉降值在(0,H)較小,最小值為18.41 cm;當在(H/3,H/2)、(H/3,3H/4)和(H/3,H)布設拉筋帶時,坡頂沉降值在(H/3,H)較小,最小值為17.08 cm;當在(H/2,3H/4)和(H/2,H)布設拉筋帶時,坡頂沉降值在(H/2,H)較小,最小值為17.59 cm,可見,在(H/3,H)布設拉筋帶坡頂沉降最小,邊坡最為穩定。由圖11(b)可知,在(0,H/3,H)布設拉筋帶時沉降最小,最小值為16.01 cm,相較于(H/3,H),減幅達到6.26%;而在(H/3,H/2,H)和(H/3,3H/4,H)布設拉筋帶時,二者差值不大。由圖11(c)可知,在(0,H/3,3H/4,H)布設拉筋帶時沉降最小,最小值為15.28 cm,相較于(0,H/3,H),減幅達到4.51%。

圖11 拉筋帶位置對坡頂沉降的影響

綜上所述,在0,H/3,3H/4,H處布設拉筋帶對于減小坡頂沉降具有顯著作用。

3.3 拉筋帶數量對邊坡動力響應的影響

基于上述分析,選取兩層拉筋帶、3層拉筋帶和4層拉筋帶布設的最優位置,以此來研究不同數量拉筋帶對于邊坡動力響應的影響,選取工況如表3所示。

表3 拉筋帶數量計算工況

3.3.1 拉筋帶數量對邊坡滑裂面的影響

圖12展示了計算終態時各工況的塑性云圖。由圖可知,在未布設拉筋帶時(工況1),邊坡的塑性區較小,主要出現在坡面附近,邊坡滑裂面近似于圓弧形;當在H/3和H處布設拉筋帶時(工況2),坡腳和坡頂附近的邊坡塑性區向坡體發展較為明顯,塑性區域增大,邊坡滑裂面向坡體轉移,但H/3～H高程之間的塑性區并未得到明顯的發展;當在0、H/3和H處布設拉筋帶時(工況3),坡腳附近的塑性區進一步向坡體延伸,且H/3～H之間的塑性區也有進一步向坡體發展的趨勢,但坡頂處的塑性區未得到明顯發展,滑裂面位置僅在坡腳處向坡體轉移;當在0、H/3、3H/4和H處布設拉筋帶時(工況4),邊坡塑性區整體向坡體明顯發展,塑性區域增大,滑裂面整體轉移至拉筋帶末端,邊坡穩定性得到提高。

圖12 各工況下邊坡塑性云圖

綜合而言,在選擇拉筋帶布設位置時,應優先考慮坡腳,坡頂以及H/3～H高程之間的塑性區發展情況,從而確定最優布設位置。

3.3.2 拉筋帶數量對格室約束圍壓的影響

圖13為不同工況下墻面格室約束圍壓沿高程分布。由圖可知,不同工況下墻面格室約束圍壓沿高程均呈“增加-減小”的兩段形態,在H/3處達到最大值,由于底層格室內充填料受到上部格室及其內充填料的壓重作用較大,對于格室內充填料的變形具有一定的限制作用,而邊坡中部受到上部格室及其充填料的壓重作用減小,且H/3處土壓力最大,沿墻高,格室壓重作用減弱,頂層格室內部充填料在地震力作用下,向坡頂臨空面運動,格室約束圍壓較小。另外,工況1中,格室約束圍壓峰值為88.92 kPa;工況2中,格室約束圍壓峰值為97.61 kPa,相較于工況1,增幅為9.77%,且H/3～H高程格室約束圍壓增幅大于0～H/3高程;工況3中,格室約束圍壓峰值為105.71 kPa,相較于工況2,增幅為8.29%,且0～H/3高程和H/3～H高程格室約束圍壓增幅基本一致;工況4中,格室約束圍壓峰值為108.39 kPa,與工況3相比,增幅為2.54%,且0～H/3高程格室約束圍壓增幅較小,而H/3～H高程增幅較為明顯,說明在3H/4處設置拉筋帶對于增強邊坡上部穩定性具有重要作用。

圖13 不同工況墻面格室約束圍壓沿高程變化

圖14展示了不同工況下拉筋帶格室約束圍壓沿水平方向變化。由圖可知,格室約束圍壓自墻面向坡體延伸呈先增大后減小的“人”字形分布,最大值均出現在擋墻與坡面交界處,由于在地震作用下,與坡面交界處的拉筋帶格室單元受到拉筋帶上下坡體的推壓作用較大,該處格室單元需要發揮更大的約束作用以低抗坡體的推壓作用,此處格室約束圍壓最大;向坡體延伸,拉筋帶格室單元受到坡體的壓重作用較大,起到主要作用,因此越向坡體格室約束圍壓越小。此外,不同工況下墻面格室約束圍壓小于墻體格室約束圍壓,由于墻體格室受到周邊格室及坡體填料的相互擠壓作用較大,限制了格室的拉伸變形,而墻面格室緊靠臨空面,在地震力作用下,受到周邊格室擠壓作用時,將向臨空面方向產生變形,因此墻體格室約束圍壓大于墻面格室約束圍壓。不同工況下格室約束圍壓由大到小依次為H/3>3H/4>H>0,說明H/3處格室受力最大。

圖14 拉筋帶處格室約束圍壓沿水平方向變化

綜合而言,拉筋帶布設位置合理時,隨著拉筋帶數量的增加邊坡的整體穩定性得到顯著提高。但拉筋帶數量布設較多時,施工難度加大,工時較長。

3.3.3 拉筋帶數量對擋墻位移的影響

圖15展示了不同工況下墻面永久水平位移沿高程變化。由圖可知,不同工況下位移沿高程均表現出明顯的非線性,在H/3處位移達到最大值。此外,工況1中,最大位移為22.02 cm,墻頂位移小于墻腳位移;工況2中,最大位移為19.93 cm,相較于工況1,減幅達到9.46%,且墻腳處位移由19.56 cm減小至8.11 cm,減幅達到58.56%;工況3中,最大位移為19.93 cm,相較于工況2,減幅達到6.44%,且0～H/3減幅最大,H/3～H高程減幅較小;工況4中,最大位移為17.51 cm,相較于工況3,0～H/3減幅較小,H/3～H高程減幅較大。4種工況下位移由大到小依次為工況1、工況2、工況3、工況4,即隨著拉筋帶數量的增加,墻面永久水平位移逐漸減小,且減幅趨弱。

圖15 不同工況下墻面永久水平位移沿高程變化

圖16展示了不同工況下邊坡水平位移云圖。由圖可知,不同工況下邊坡水平位移最大值均出現在H/3處。未布設拉筋帶時(工況1),坡面和坡體均出現了位移,說明擋墻結構抗震性能較弱;當在H/3和H布設拉筋帶時(工況2),位移僅出現在墻體及拉筋帶錨固區域,坡腳處水平位移得到明顯減小;當在0、H/3和H布設拉筋帶時(工況3),墻腳水平位移進一步減小,H/3～H高程之間的位移等勢線向坡體發展;當在0、H/3、3H/4和H處布設拉筋帶時(工況4),墻腳附近水平位移未得到明顯改變,H/3～3H/4高程之間的位移等勢線繼續向坡體發展。綜合而言,拉筋帶數量增加,墻面永久水平位移逐漸減小,邊坡抗震性得到顯著提高。

圖16 不同工況下邊坡水平位移云圖

3.3.4 拉筋帶數量對擋墻加速度的影響

圖17展示了墻面水平峰值加速度沿高程分布。由圖可知,不同工況下峰值加速度沿高程分布大體相同,均表現為“增加-衰減-增加”的三段形態,這與文獻[22]研究結果一致。這是由于在墻腳附近(H<50 cm),擋墻結構對地震波的衰減效應并不明顯,臨空面放大效應仍占據主導因素,從而出現了加速度隨高程而增加的現象;墻高到達一定高度時(H=50 cm左右),隨著高程增加,土工格室及內部充填料對地震波的阻尼增大,濾波作用增強,入射波經過格室內充填料的吸收濾波后,阻尼衰減效應明顯增高,從而超過臨空面放大效應而占據主導地位,此時加速度隨高程增加而減小;當高程增加至H=110 cm時,沿高程至墻頂,由于墻頂臨空面的存在,入射波在臨空面處經過多次反射、折射和疊加,臨空面放大效應超過阻尼衰減效應再次占據主導因素,從而加速度沿高程的增加而增加。另外,4種工況下,峰值加速度由大到小依次為工況1、工況2、工況3、工況4,即隨著加筋帶數量的增加,峰值加速度逐漸減小。

圖17 墻面水平峰值加速度沿高程變化

圖18為不同工況下拉筋帶處峰值加速度沿水平分布。由圖可知,不同工況下峰值加速度自墻面向坡體延伸,均表現為先增大后減小的“人”字形分布,最大值出現在墻體與坡面交界處。這是由于地震波沿水平方向穿過土工格室單元時,土工格室及其內部充填料將會對地震波進行反射和透射,從而造成地震能量的衰減,出現了墻面峰值加速度小于墻體峰值加速度;而位于坡體內部的拉筋帶,其四周被坡體材料所包裹,地震波沿水平方向傳播至坡面時,坡面放大效應占據主導位置,出現了坡面峰值加速度大于坡體峰值加速度。不同位置的拉筋帶,峰值加速度由高到低依次為H>3H/4>H/3>0,說明在坡頂處加速度放大較為顯著。綜上所述,隨著拉筋帶數量的增加,峰值加速度逐漸減小;不同位置的拉筋帶,坡頂處峰值加速度最大。

圖18 不同工況下拉筋帶加速度沿水平方向變化

3.3.5 拉筋帶數量對坡頂沉降的影響

圖19為不同工況坡頂沉降分布。由圖可知,不同工況下,自坡肩向坡體深處延伸,沉降均呈先增大后減小的“V”型分布,在墻體與坡體交界處沉降最大;由于格室的約束作用以及下層格室擋墻的阻礙作用,坡肩處沉降最小;向坡體延伸,與擋墻接觸處的坡體填料隨墻體位移最先向下滑落,使得坡體頂部與頂層格室之間產生空腔,在地震荷載作用下,頂層格室及其內充填料整體下沉,并達到最大值;繼續向坡體深處延伸,坡體填料滑落趨勢漸緩,沉降減小。4種工況下,沉降由大到小對應的沉降值依次為18.28、17.08、16.01、15.29 cm,減幅依次為6.56%、6.28%、4.51%?？梢?沉降隨拉筋帶數量的增加而降低,但減幅逐漸減小。

圖19 不同工況坡頂沉降分布

3.4 土工格室擋墻抗震性能分析

土工格室擋墻是將土工格室和充填料組成的結構層按照一定的坡度層層疊加組合而成。格室擋墻的格室本身、充填料與格室之間、層與層之間,具備可靠的抗震構造措施[23]。

土工格室擋墻具有足夠的豎向承載能力,土工格室對充填料有著較好的約束作用,使得擋墻結構整體剛度較大,結構穩定性較好,在地震時能夠支撐上部結構的所有重量和使用載荷;土工格室擋墻結構具有可變的水平剛度,地震力較小時格室層間摩阻力以及格室與充填料的相互作用具有足夠的剛度,以抵抗邊坡的破壞,地震力較大時,擋墻水平剛度較小,格室層間產生水平滑動,隨地震動發生協調變形;格室擋墻結構具有一定的水平彈性恢復力,地震中格室層間的咬合摩阻力作用,使擋墻結構具有自動復位功能,地震力結束后,擋墻中下部發生撓曲變形,擋墻整體結構未發生坍塌等工程事故;土工格室擋墻結構具有足夠的阻尼,土工格室及其充填料對地震能量具有顯著耗能作用。

4 結 論

本文利用FLAC3D軟件,建立土工格室擋墻支護邊坡數值模型,系統分析了拉筋帶位置和拉筋帶數量對邊坡動力響應的影響,得到主要結論如下:

1)不同方案下,格室約束圍壓沿墻高均表現為“增加-減小”的兩段形態;拉筋帶布設位置合理時,位移最大值出現在H/3處;水平峰值加速度沿墻高呈“增加-減小-增加”的三段分布,坡頂處加速度最大;沉降自坡肩向坡體延伸呈先增加后減小的“V”形分布,在墻體與坡面交界處沉降最大。

2)布設兩層拉筋帶時,在(H/3,H)處為最優位置;布設3層拉筋帶時,在(0,H/3,H)處為最優位置;布設4層拉筋帶時,在(0,H/3,3H/4,H)處為最優位置。

3)不同工況下,拉筋帶處格室約束圍壓和水平峰值加速度自墻面向坡體延伸均表現為先增大后減小的“人”字形分布,最大值均出現在擋墻和坡面交界處;且墻面加速度小于墻體加速度,土工格室及其填料對于加速度的傳播具有一定衰減作用。

4)不同工況下,拉筋帶位置處格室約束圍壓最大值出現在H/3處,峰值加速度最大值出現在H處;不同布設方案下,在H/3和H處布設拉筋帶時,邊坡位移和沉降最小,實際工程中應優先考慮在H/3和H處布設。