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

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

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

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

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

1 格室約束圍壓理論解

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

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

(1)

(2)

(3)

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

2 邊坡模型建立

2.1 數(shù)值模型建立

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

圖2 邊坡數(shù)值模型(m)

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

表2 接觸面力學參數(shù)

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

為了避免地震波傳播時產(chǎn)生的反射及能量耗散,本次數(shù)值模型的四周及底部邊界均設(shè)置為自由場邊界。阻尼選取局部阻尼,阻尼系數(shù)為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 監(jiān)測點布置

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

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

圖4 監(jiān)測點布置

3 計算結(jié)果分析

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

圖5 拉筋帶及其監(jiān)測點布置示意圖

3.1 數(shù)值模擬結(jié)果驗證

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

圖6 第1次振動數(shù)值模擬與試驗結(jié)果對比

圖7 第3次振動數(shù)值模擬與試驗結(jié)果對比

3.2 拉筋帶位置對邊坡動力響應(yīng)的影響

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

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

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

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

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

綜上所述,在實際工程中應(yīng)優(yōu)先考慮在0,H/3,3H/4和H這4個位置布設(shè)拉筋帶。

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

圖9(a)為布設(shè)兩層拉筋帶時墻面永久水平位移沿高程變化。由圖可知,當在(0,H/3)和(0,H/2)布設(shè)拉筋帶時,位移沿高程呈非線性增加,在坡頂處位移達到最大值;而在(0,3H/4)和(0,H)布設(shè)拉筋帶時,位移沿高程呈“增加-減小”的兩段形態(tài),最大位移出現(xiàn)在H/3處,說明拉筋帶布設(shè)位置上移能夠?qū)ζ马敻浇奈灰破鸬接行Ф糁谱饔?當在(H/3,H/2)布設(shè)拉筋帶時,位移最大值出現(xiàn)在坡頂處,而在(H/3,3H/4)和(H/3,H)布設(shè)拉筋帶時,位移最大值又出現(xiàn)在H/3處,由此說明,僅在邊坡中下部布設(shè)拉筋帶時最大位移出現(xiàn)在坡頂附近,而僅在邊坡中上部布設(shè)拉筋帶時最大位移出現(xiàn)在H/3處;當在(H/2,3H/4)和(H/2,H)布設(shè)拉筋帶時,邊坡中下部位移值二者基本接近,但中上部位移逐漸減小,而在(3H/4,H)布設(shè)拉筋帶時,坡腳附近的位移大于坡頂附近。

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

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

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

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

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

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

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

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

綜上所述,在0、H/3、3H/4和H處布設(shè)時,峰值加速度放大效應(yīng)最小,邊坡最為穩(wěn)定。

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

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

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

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

3.3 拉筋帶數(shù)量對邊坡動力響應(yīng)的影響

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

表3 拉筋帶數(shù)量計算工況

3.3.1 拉筋帶數(shù)量對邊坡滑裂面的影響

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

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

綜合而言,在選擇拉筋帶布設(shè)位置時,應(yīng)優(yōu)先考慮坡腳,坡頂以及H/3～H高程之間的塑性區(qū)發(fā)展情況,從而確定最優(yōu)布設(shè)位置。

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

圖13為不同工況下墻面格室約束圍壓沿高程分布。由圖可知,不同工況下墻面格室約束圍壓沿高程均呈“增加-減小”的兩段形態(tài),在H/3處達到最大值,由于底層格室內(nèi)充填料受到上部格室及其內(nèi)充填料的壓重作用較大,對于格室內(nèi)充填料的變形具有一定的限制作用,而邊坡中部受到上部格室及其充填料的壓重作用減小,且H/3處土壓力最大,沿墻高,格室壓重作用減弱,頂層格室內(nèi)部充填料在地震力作用下,向坡頂臨空面運動,格室約束圍壓較小。另外,工況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處設(shè)置拉筋帶對于增強邊坡上部穩(wěn)定性具有重要作用。

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

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

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

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

3.3.3 拉筋帶數(shù)量對擋墻位移的影響

圖15展示了不同工況下墻面永久水平位移沿高程變化。由圖可知,不同工況下位移沿高程均表現(xiàn)出明顯的非線性,在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,即隨著拉筋帶數(shù)量的增加,墻面永久水平位移逐漸減小,且減幅趨弱。

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

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

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

3.3.4 拉筋帶數(shù)量對擋墻加速度的影響

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

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

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

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

3.3.5 拉筋帶數(shù)量對坡頂沉降的影響

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

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

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

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

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

4 結(jié) 論

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

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

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

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

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