苗滿(mǎn)勝，劉 軍

(河南高速公路發(fā)展有限責(zé)任公司， 河南 鄭州 450000)

土工格室是20世紀(jì)80年代初在國(guó)際上出現(xiàn)的一種新型土工合成材料，其材質(zhì)輕、耐磨損、伸縮自如、施工速度快并提供較高的側(cè)向限制，廣泛應(yīng)用于土木工程領(lǐng)域。由于其結(jié)構(gòu)輕、施工方便、造價(jià)低廉，能夠滿(mǎn)足綠化墻面的要求，在公路邊坡防護(hù)中得到了廣泛的應(yīng)用[1-7]。一些學(xué)者針對(duì)土工格室柔性擋墻這種新型邊坡支護(hù)結(jié)構(gòu)展開(kāi)了相關(guān)的研究，取得了一些成果。謝永利等[8]運(yùn)用有限元分析軟件MARC對(duì)土工格室柔性擋墻的工程性狀進(jìn)行了彈塑性數(shù)值仿真分析。宋飛等[9]運(yùn)用PLAXIS軟件研究了在不同高寬比、坡度及路基表面荷載下對(duì)土工格室擋墻變性性狀的影響。馬新巖等[10]采用數(shù)值模擬和結(jié)構(gòu)計(jì)算相結(jié)合的方法，結(jié)合實(shí)際工程，對(duì)路塹土工格室柔性擋墻進(jìn)行實(shí)例分析，并根據(jù)施工經(jīng)驗(yàn)總結(jié)柔性擋墻施工工序。王根等[11]以瑞鷹高速土工格室柔性擋墻工程為背景，采用ADINA有限元數(shù)值仿真分析，考慮墻體與路基土的耦合關(guān)系，研究了土工格室多級(jí)柔性擋墻的受力機(jī)理，對(duì)比多級(jí)擋墻在多種工況下的受力情況。楊利等[12]分析了土工格室加筋粗粒土的效果分析。劉軍等[13]結(jié)合工程分析了土工格室柔性擋墻的敏感因素。而在土工格室擋墻設(shè)計(jì)過(guò)程中，土工格室柔性擋墻在不同設(shè)計(jì)參數(shù)下的受力及變形規(guī)律仍不清楚，因此土壓力計(jì)算方法缺少合理的理論基礎(chǔ)，理論研究仍然落后于工程實(shí)踐，需要進(jìn)行深入的研究。

本文采用巖土工程有限差分軟件FLAC3D計(jì)算分析了設(shè)計(jì)參數(shù)下土工格室柔性擋墻的受力及變形特征，研究了厚度、坡比以及彈性模量對(duì)于柔性擋墻受力及變形性狀的影響。為柔性擋墻土壓力計(jì)算方法的提出提供了理論基礎(chǔ)，對(duì)于工程設(shè)計(jì)具有一定的參考意義。

1 模型的建立及相關(guān)參數(shù)

1.1 模型的建立

土工格室柔性擋墻數(shù)值計(jì)算模型包括幾何模型與材料模型。本文根據(jù)工程實(shí)際問(wèn)題和研究目的需要建立柔性擋墻計(jì)算模型。計(jì)算模型是由擋墻墻身、基礎(chǔ)、墻后填土三部分組成。經(jīng)過(guò)查閱參考資料[13]以及試算，確定各部分的模型尺寸為：擋墻長(zhǎng)度15 m，高度為6 m，厚度為3 m，擋墻設(shè)計(jì)坡率為1∶0.5；基礎(chǔ)寬度為4 m，厚度為1 m，墻后土體寬度取15 m。其三維實(shí)體幾何模型與邊界示意圖分別如圖1、圖2所示。由于該模型采用整體擋墻的一部分進(jìn)行分析，所以需要確定擋墻模型邊界條件。擋土墻前立面和后立面邊界施加X(jué)方向位移約束為零，擋土墻前后兩兩個(gè)側(cè)面邊界施加Z方向位移約束為零，柔性擋墻底面邊界在X、Y、Z三個(gè)方向施加位移約束為零。

1.2 網(wǎng)格劃分

FLAC3D數(shù)值分析過(guò)程中，單元?jiǎng)澐衷蕉?，?jì)算結(jié)果越精細(xì)，但是單元?jiǎng)澐诌^(guò)多，會(huì)增加計(jì)算機(jī)負(fù)擔(dān)，導(dǎo)致計(jì)算過(guò)程緩慢。綜合考慮各種因素，本次計(jì)算單元網(wǎng)格劃分采用“掃略網(wǎng)格劃分”，將擋墻墻體作為重點(diǎn)研究對(duì)象，進(jìn)行網(wǎng)格加密處理；基礎(chǔ)部分和土體部分作為非研究重點(diǎn)，網(wǎng)格劃分時(shí)相對(duì)疏一些。具體單元網(wǎng)格劃分如圖3。整個(gè)模型共17 904個(gè)節(jié)點(diǎn)，15 585個(gè)單元。

圖1 幾何模型示意圖(單位：m)

圖2 邊界示意圖

圖3網(wǎng)格劃分立體圖

1.3 模型參數(shù)的選取

在土工格室柔性擋墻中，土體本構(gòu)模型Mohr-Coulomb理想彈塑性模型，該理論比較完善且簡(jiǎn)單實(shí)用，被廣泛應(yīng)用于巖土工程和理論實(shí)踐中[14]。FLAC3D數(shù)值分析軟件中的Mohr-Coulomb理想彈塑性模型需要輸入5個(gè)參數(shù)，分別為黏聚力c、內(nèi)摩擦角φ、密度ρ、剪切模量G和體積模量K。其中黏聚力c、內(nèi)摩擦角φ和密度ρ可以通過(guò)土樣的常規(guī)試驗(yàn)來(lái)獲取，而剪切模量G和體積模量K可根據(jù)彈性模量E和泊松比μ求得，其計(jì)算方法如式(1)所示。

(1)

參考工程實(shí)際勘察報(bào)告及課題組試驗(yàn)結(jié)果對(duì)柔性擋墻強(qiáng)度參數(shù)取值，擋墻各部分具體參數(shù)見(jiàn)表1。

表1 模型材料計(jì)算參數(shù)

1.4 接觸面優(yōu)化

本文對(duì)擋墻墻體、地基基礎(chǔ)以及墻后填土相連接的結(jié)構(gòu)面采用接觸面形式處理，即采用FLAC3D軟件中自帶的Interface接觸單元來(lái)模擬[15-16]。Interface接觸單元厚度為零，但具有方向性，是通過(guò)將三角形面積單元分配到各節(jié)點(diǎn)上來(lái)表示的，在設(shè)置接觸面單元時(shí)需要Interface單元平鋪到事先確定的目標(biāo)面上，并與另一個(gè)接觸面連接在一起。

接觸面上受到的切向力小于最大切向力時(shí)處于彈性階段，可通過(guò)式(2)和式(3)來(lái)確定其某一時(shí)刻的法向力與切向力。

Fn(t+Δt)=KnμnA+σnA

(2)

(3)

當(dāng)|Fs|=Fsmax，接觸面滑動(dòng)時(shí)，切向力保持不變，進(jìn)入塑性變形階段，法向應(yīng)力隨著位移增大而增大，可按式(4)計(jì)算。

(4)

2 工程實(shí)例計(jì)算

2.1 初始應(yīng)力平衡

土工格室柔性擋墻建造之前，施工場(chǎng)地地表位移為零，但是土體內(nèi)的應(yīng)力卻是一直存在的，這種沒(méi)有位移只有應(yīng)力存在的狀態(tài)即地應(yīng)力平衡狀態(tài)。數(shù)值模擬過(guò)程中，初始應(yīng)力平衡是為了使數(shù)值模型獲得一個(gè)存在初始應(yīng)力但初始應(yīng)變?yōu)榱愕臓顟B(tài)。考慮到柔性擋墻施工場(chǎng)地中自重應(yīng)力是主要初始應(yīng)力，顯然其自重應(yīng)力是外力，地應(yīng)力場(chǎng)是將自重應(yīng)力下的場(chǎng)地內(nèi)力提取出來(lái)施加于數(shù)值模型上，開(kāi)始工程數(shù)值計(jì)算時(shí)再施加自重應(yīng)力，使得場(chǎng)地內(nèi)的內(nèi)力與外力平衡從而獲得較精確的、沒(méi)有受到人為因素干擾情況下的數(shù)值模型初始狀態(tài)。

邊坡坡體在自重應(yīng)力下的水平有效應(yīng)力與水平位移見(jiàn)圖4和圖5。

2.2 水平土壓力及變形對(duì)比分析

基于前文分析成果，并將數(shù)值模擬計(jì)算結(jié)果同現(xiàn)場(chǎng)測(cè)試和模型試驗(yàn)進(jìn)行對(duì)比。土壓力數(shù)值模擬結(jié)果與對(duì)比曲線(xiàn)見(jiàn)圖6、圖7。水平位移數(shù)值模擬與曲線(xiàn)見(jiàn)圖8、圖9。

圖4 重力作用下的初始水平有效應(yīng)力分布

圖5 重力作用下的初始水平位移分布

圖6 柔性擋墻墻體水平受力云圖

圖7 柔性擋墻水平土壓力數(shù)值模擬計(jì)算結(jié)果

圖8柔性擋墻墻體水平位移云圖

從圖7中可以看處，柔性擋墻水平土壓力數(shù)值計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)值存在一定差距，數(shù)值計(jì)算結(jié)果較實(shí)際偏小。數(shù)值計(jì)算解得最大土壓力值為55.32 kPa，位于擋墻高度2.5 m位置處；而現(xiàn)場(chǎng)實(shí)測(cè)土壓力值最大為57.62 kPa，且位于墻高1/2H處(擋墻高度為3 m)。這主要是由于數(shù)值分析在理想狀態(tài)下建立計(jì)算模型和設(shè)定邊界條件，無(wú)法完全準(zhǔn)確模擬擋墻具體和施工機(jī)械影響，所以現(xiàn)場(chǎng)實(shí)際測(cè)量存在一定差異，而且現(xiàn)場(chǎng)實(shí)際測(cè)量除了容易受周邊環(huán)境影響，還會(huì)出現(xiàn)一些不可避免的人為操作誤差。

從圖9中可以看到，柔性擋墻墻體水平位移數(shù)值模擬計(jì)算結(jié)果曲線(xiàn)與現(xiàn)場(chǎng)測(cè)試和模型試驗(yàn)基本一致，均以墻高1/2H為界，分界以上水平位移隨墻高減小，分界以下沿墻高增大。數(shù)值模擬算得最大墻體位移為4.52 cm，比現(xiàn)場(chǎng)測(cè)試結(jié)果要小15.04%，比模型試驗(yàn)要大2.03%。則說(shuō)明本文數(shù)值模擬計(jì)算結(jié)果與現(xiàn)場(chǎng)測(cè)試和模型試驗(yàn)結(jié)果在數(shù)值上差距不大，且曲線(xiàn)發(fā)展趨勢(shì)也基本相似，因此可以認(rèn)為本文數(shù)值分析在墻體水平位移方面比較接近實(shí)際。

圖9柔性擋墻水平位移數(shù)值模擬計(jì)算結(jié)果

3 柔性擋墻設(shè)計(jì)參數(shù)影響分析

根據(jù)前文分析可知，了解各個(gè)擋墻設(shè)計(jì)參數(shù)對(duì)于土工格室柔性擋墻的受力及變形的影響程度，對(duì)于擋墻設(shè)計(jì)及施工有著重要意義。本文選取擋墻厚度、擋墻坡比、擋墻彈性模量三個(gè)參數(shù)設(shè)立3種工況展開(kāi)分析，具體如表2所示。以期把握此類(lèi)參數(shù)對(duì)于擋墻受力與變形的影響，從而對(duì)擋墻設(shè)計(jì)進(jìn)行優(yōu)化，對(duì)擋墻土壓力計(jì)算提供依據(jù)。

表2 數(shù)值模擬各工況設(shè)計(jì)表

3.1 擋墻厚度影響分析

圖10為墻高為6 m，坡比為1∶0.5的擋柔性墻，墻身厚度在2 m到5 m范圍內(nèi)變化時(shí)的墻背水平應(yīng)力及位移。

圖10柔性擋墻墻體厚度變化影響圖

從圖10(a)中可以看出，高為6 m的柔性擋墻墻體厚度不同時(shí)，土壓力沿?fù)鯄Ω叨确植疾町愝^大。當(dāng)墻體厚度為2 m時(shí)，擋墻水平土壓力整體較小，但是在墻高為3 m位置處，土壓力急劇增大，達(dá)到69.73 kPa；當(dāng)墻體厚度為3 m時(shí)土壓力分布較均勻，在墻高2.5 m位置處達(dá)到最大，最大為55.32 kPa；當(dāng)墻體厚度為4 m和5 m時(shí)，擋墻水平土壓力值相差不大，均呈不規(guī)則的線(xiàn)性分布，且隨擋墻高度降低逐漸變大。

從圖10(b)中可以看出，不同擋墻厚度下墻背水平位移曲線(xiàn)與土壓力曲線(xiàn)具有一定相似性。墻厚為2 m時(shí)擋墻變形較大，尤其在墻高為3 m位置處達(dá)到6.34 cm，變形量大于擋墻厚度3%，說(shuō)明2 m的擋墻厚度并不能滿(mǎn)足該擋墻安全設(shè)計(jì)要求。墻體厚度為3 m、4 m和5 m時(shí)，擋墻墻體水平位移相差不大，這也說(shuō)明在滿(mǎn)足擋墻安全設(shè)計(jì)要求基礎(chǔ)上，盲目增加擋墻厚度對(duì)于擋墻變形的控制有限，反而只會(huì)增大工程造價(jià)。

3.2 擋墻坡比影響分析

圖11為墻高為6 m、墻厚為3 m的土工格室柔性擋墻，在變化擋墻坡比情況下的墻背土壓力及水平位移曲線(xiàn)。

圖11柔性擋墻坡比變化影響圖

從圖11(a)可以看出，柔性擋墻坡度越陡，擋墻水平土壓力越大。擋墻坡度為1∶0.75和1∶0.50時(shí)，擋墻水平土壓力分布曲線(xiàn)相似，均在墻高2.5 m處分界，分界以上部分，土壓力沿墻高線(xiàn)性減小，分界以下部分恰好相反，分界點(diǎn)處土壓力值最大，分別為55.32 kPa和49.96 kPa；柔性擋墻坡度為1∶0.25時(shí)，土壓力沿墻高呈線(xiàn)性分布，在擋墻底部達(dá)到最大為60.11 kPa。

結(jié)合圖11(a)和圖11(b)可以看出，擋墻坡比在1∶0.75到1∶0.50之間時(shí)，坡比對(duì)于擋墻水平土壓力及墻體位移影響不大，說(shuō)明在擋墻設(shè)計(jì)過(guò)程中，可適當(dāng)?shù)卦谄卤?∶0.75到1∶0.50之間提高擋墻坡度，從而減少擋墻占地面積；但是當(dāng)擋墻設(shè)計(jì)坡比大于1∶0.50時(shí)，要想繼續(xù)增大擋墻坡來(lái)減少占地面積度則需要謹(jǐn)慎，因?yàn)槠卤葹?∶0.25的柔性擋墻墻體水平位移較大，且沿?fù)鯄Ω叨瘸示€(xiàn)性增大，最大位移點(diǎn)達(dá)到5.57 cm，說(shuō)明該柔性擋墻設(shè)計(jì)較陡時(shí)，其受力形態(tài)和位移變形規(guī)律產(chǎn)生了較大變化，在施工過(guò)程中及施工完成后，擋墻可能還發(fā)生了繞墻頂轉(zhuǎn)動(dòng)的位移。

3.3 擋墻彈性模量影響分析

圖12為墻體彈性模量對(duì)墻高為6 m、坡比1∶0.50、擋墻厚度為3 m的土工格室柔性擋墻水平方向土壓力及墻體位移的影響。圖中模量是指墻體彈性模量，取值范圍為5×106Pa～5×1010Pa。

圖12柔性擋墻墻身彈性模量對(duì)土壓力及變形的影響

從圖12(a)中可以看出，擋墻水平土壓力隨著擋墻墻體彈性模量的增大而增大。墻體彈性模量在5×106Pa到5×108Pa之間時(shí)，擋墻水平土壓力近似呈鼓形分布，只是最大土壓力值所在擋墻位置不同，墻體彈性模量為5×106Pa和5×107Pa時(shí)，最大土壓力值位于墻高2.5 m處，分別為43.21 kPa和55.32 kPa；墻體彈性模量為5×108Pa時(shí)，最大土壓力位于墻高3 m處為61.32 kPa。墻體彈性模量為5×109Pa和5×1010Pa時(shí)，墻體材料更接近剛性材料，擋墻土壓力曲線(xiàn)也與剛性擋墻相似，在墻高1/3H附近處土壓力值最大，分別為71.19 kPa和84.76 kPa。

從圖12(b)中可以看出，柔性擋墻墻體水平位移隨擋墻墻體彈性模量的增大而降低。墻體彈性模量在5×105Pa～5×107Pa之間，擋墻墻體位移呈鼓形分布，最大位移點(diǎn)位于擋墻中部達(dá)到9.21 cm，是擋墻厚度的3.07%，不滿(mǎn)足設(shè)計(jì)安全要求，說(shuō)明擋墻材料彈性模量過(guò)低，會(huì)導(dǎo)致?lián)鯄Ξa(chǎn)生較大變形，因此擋墻施工前一定要進(jìn)行對(duì)格室材料進(jìn)行檢查驗(yàn)收；墻體彈性模量在5×106Pa～5×107Pa之間時(shí)，擋墻水平位移規(guī)律基本一致，呈明顯非線(xiàn)性分布，在墻高1/2H附近處達(dá)到最大；墻體彈性模量大于5×108Pa時(shí)，擋墻墻體水平位移曲線(xiàn)發(fā)生了很大變化，隨擋墻高度線(xiàn)性增加，在墻頂位置達(dá)到最大。

4 結(jié) 論

本文應(yīng)用FLAC3D軟件對(duì)土工格室柔性擋墻進(jìn)行了數(shù)值模擬研究，分析了不同設(shè)計(jì)參數(shù)對(duì)于6 m高的路肩式土工格室柔性擋墻的受力及變形的影響，并得出了如下結(jié)論：

(1) 6 m高的柔性擋墻墻體厚度為2 m～3 m時(shí)，擋墻水平土壓力和墻體位移均呈鼓形分布，最大值位于墻高1/2H附近處；墻體厚度大于3 m時(shí)，墻體水平位移變小，墻背土壓力及位移曲線(xiàn)也逐漸過(guò)渡為線(xiàn)性模式。

(2) 1∶0.50和1∶0.75兩種坡率的柔性擋墻其擋墻水平土壓力及墻體位移曲線(xiàn)相似，數(shù)值相差不大，且都以墻高1/2H分界；當(dāng)柔性擋墻坡比為1∶0.25時(shí)，擋墻水平土壓力值及位移值，都出現(xiàn)了不同程度的增加，且位移曲線(xiàn)也近似變?yōu)橹本€(xiàn)，在墻頂位置達(dá)到最大，為5.57 cm。

(3) 6 m高的柔性擋墻水平土壓力隨著墻體彈性模量的增大而增大，水平位移恰好相反，隨墻體彈性模量增大而減小。當(dāng)體彈性模量在5×106Pa～5×108Pa之間時(shí)，擋墻土壓力及墻體水平位移呈鼓形分布，在墻高1/2H附近處達(dá)到最大；當(dāng)墻體彈性模量在5×109Pa～5×1010Pa之間時(shí)，水平土壓力及位移曲線(xiàn)發(fā)生明顯轉(zhuǎn)變，位移曲線(xiàn)大致呈線(xiàn)性分布，土壓力也變得與剛性擋墻所類(lèi)似，在墻高1/3H處達(dá)到最大。