程聯(lián)偉

(中鐵第五勘察設計院集團有限公司,北京102600)

0 前 言

本文以西安某深基坑支護工程為背景,結合FLAC3D有限差分軟件,建立了空間三維有限差分分析模型,對預應力錨索復合土釘支護結構施工進行了三維動態(tài)數(shù)值模擬,分析總結了在基坑分步開挖過程中,深基坑側壁土壓力與變形;支護結構的受力分布、水平位移;地表面的沉降;坑底的隆起空間分布情況及變化規(guī)律。

1 工程概況

擬建工程建筑面積約為87 000 m2,基坑開挖深度為-16.8 m。根據(jù)鉆探揭露,場地地層在地基設計影響范圍內分為6層:第1層為填土(Q4ml);第2層為黃土(Q3eol);第3層為黃土(Q2eol);第4層為古土壤(Q3el);第5層為粉質粘土(Q2al);第6層為粉質粘土(Q2al)。各土層物理力學性質見表1。在勘察期間,實測場地地下水穩(wěn)定水位埋深7.30 m～9.80 m,屬潛水類型。

2 建立計算模型及計算分析

2.1 建立計算模型

依據(jù)Duncan與Goodman建議邊坡的計算范圍取邊坡坡腳以下邊坡高度H深、坡腳左右水平2H～3H的范圍。本文采用70 m×3.75 m×33.6 m的計算模型;土釘、錨索單元采用FLAC3D中的錨索結構單元,土釘及錨桿設計參數(shù)見表2;混凝土面層采用FLAC3D中的殼單元來模擬,用戶可以根據(jù)計算的精確程度定義三角形單元的大小、殼結構單元的厚度、材料屬性以及殼結構與土體之間的作用方式。通過對不同位置上的節(jié)點設置監(jiān)測歷史記錄,可以得到作用于噴射混凝土面層的土壓力變化過程;腰梁采用FLAC3D中的梁結構單元進行模擬。

表1 土層物理力學指標

表2 土釘與錨桿的設計參數(shù)

2.2 基坑開挖施工過程模擬

本文采用FLAC3D軟件對基坑開挖施工過程的模擬采用以下步驟:①基坑第一步開挖,設置第1排土釘,設置坑壁噴護面層,進行模型平衡計算。②基坑第二步開挖,設置第2排土釘,設置坑壁噴護面層,進行模型平衡計算。③分別進行第三、第四等各步依此類推,每步超挖深度為0.5 m,依次開挖支護至坑底。

2.3 結構水平位移特征

在基坑開挖設計過程中,基坑支護結構的水平位移的大小和范圍如何隨基坑開挖而變化是設計者最關心的問題,研究每個工序對支護結構的水平位移的影響變化是非常必要的。圖1為開挖結束后支護結構的水平位移等值線圖,圖2為在FLAC3D模型監(jiān)測到的基坑坑壁隨每步施工的水平位移隨深度變化分布曲線,圖3為采用FLAC3D進行模擬的結果與實測結果水平位移對比曲線。

圖1 基坑水平位移等值線圖

規(guī)律總結與結果分析[1-3]:①在基坑開挖支護前期,基坑變形主要是由于基坑開挖后應力不均,在應力重分布過程中引起基坑的初步變形,此時變形最大位置在開挖面以下的土體中。特別是前三步開挖,以此變形為主。隨開挖的進行,基坑變形的增大,此類變形逐漸失去主導地位。②基坑坑壁的水平變形,隨開挖進行逐漸增大。每一步開挖均引起一定的水平位移增值,并且每一步的開挖所引起的邊坡水平位移增量相差較小。同時,坑壁的最大水平位移增量常出現(xiàn)在上一步開挖面下方不遠處,由此可見,每次開挖影響最大的位置在臨近的臨空面附近。③距離坑壁較遠的位置土體的水平位移近似為零,遠離坑壁方向的土體的位移逐漸變小,靠近開挖面的土體水平位移較大?；幼畲笞冃伟l(fā)生在約距基坑底面基坑開挖深度處,在支護設計過程中應充分考慮這一特點,為更好控制基坑變形需增強此處支護結構的設計。④設置土釘及錨桿位置的坑壁水平位移明顯減小,坑壁水平位移沿深度變化同一高度的土體距土釘較近位置的水平位移較小,越遠離土釘?shù)奈恢?水平位移越變大。這一現(xiàn)象說明土釘對土體側向變形限制作用非常明顯。⑤由于土釘及錨桿對被加固土體的約束作用,能較好地限制基坑變形。經過土釘及錨桿加固后的復合土體內部的水平位移相對較小,表明土釘錨桿及土體一個彼此相互作用的有機整體,因此在結構分析過程中不應將各部分簡單的分開考慮。⑥土釘支護結構的柔性特點決定著支護結構抵抗側向位移較弱,每步的開挖深度雖然很小,由于側向及底部約束的解除,盡管坑壁仍然能夠保持自立穩(wěn)定,邊坡還是產生一定的水平位移。開挖深度除第一步為2.0 m外,其他各步均為1.5 m,卸載量變化不大,使得坑壁水平變形曲線形狀非常相似。

圖2 模型監(jiān)測水平位移變化曲線

采用FLAC3D分析結果與實測基坑水平位移存在差異的主要原因是:①由于基坑水平位移觀測基準點的位置位于基坑變形影響范圍內,使得所觀測的數(shù)據(jù)存在誤差;②由于基坑開挖前的降水使得坑壁土體水份流失,力學性質改變等影響因素為在FLAC3D計算時考慮;③采用網格進行區(qū)域劃分,采用網格節(jié)點上離散場值來代替區(qū)域內連續(xù)分布的場,用有限個網格節(jié)點來代替連續(xù)的求解區(qū)域帶來的誤差。

圖3 模擬與實測水平位移變化曲線

2.4 結構垂直位移特征

在此工程中,引起地面沉降的主要因素是地下降水。但是,在模擬過程中僅考慮降水結束后,由基坑開挖所引起的地面沉降大小及變化規(guī)律,如圖4所示。

圖4 基坑開挖垂直位移等值線圖

規(guī)律總結與結果分析:①隨著基坑的開挖,天然土層的原始平衡狀態(tài)被破壞,形成二次應力場[1]?；娱_挖影響范圍內應力重新分布,坑腳處應力集中分布?；拥酌嬗捎陂_挖卸載,地面土體出現(xiàn)一定量的回彈。新形成的坑壁有向下移動的運動趨勢,對卸載后的基坑底面土體產生擠壓作用,從而土體形成擠壓隆起現(xiàn)象。②基坑地表的豎直位移沿遠離坑壁方向呈曲線分布,每步開挖上一層土體的沉降量最大,與水平位移特征相同,每次開挖地表的變化曲線非常相似?；拥孛尕Q向變形較小,在土釘支護范圍內表現(xiàn)十分明顯,由此可見,土釘及錨桿的設置不但可以增加土體的整體強度,而且抑制了土體豎向變形的趨勢[4-5]。

2.5 土釘及錨桿受力特征

在模型分析中,土釘和錨桿均用cable單元進行模擬,假設結構單元的拉力在一個小的范圍內是均勻分布的,在此模型中每根土釘按照Δ z=1.0 m、每根錨桿按照Δz=0.955 m的原則等分為若干段,拉力大小相等。計算每個節(jié)點的應力在每個時間步長中的變化,并進行積分得到每根土釘?shù)膽Ψ植肌?/p>

規(guī)律總結與結果分析:①結構中各排土釘承受的拉力在土釘長度方向呈相似曲線分布。每一步開挖均會引起土釘拉力的一定增量,土釘每段的增量近似相同,曲線整體分布趨勢相近。各排錨桿的自由段所承受的拉力相同為水平直線分布,隨后呈現(xiàn)近似直線減小,此直線段部分很短,然后變?yōu)橥瓜蚯€減小為零。②對同一排土釘?shù)慕Y果分析表明,土釘中間靠近土釘與面層連接的一端承受的拉力最大,向兩端方向受力逐步遞減。各排土釘在靠近面層的一端,土釘拉力變化較小,均接近于一個定值;在土釘?shù)哪┒?土釘承受的拉力均遞減至接近于零,并且越接近土釘末端所承受的拉力減小速度越快。由此可知土體變形較大的地方常是承受拉力較大的地方,也是坑壁最危險的地方。③在基坑開挖過程中,每一排土釘或錨桿承受的拉力增量最大是在其下一步開挖和支護結構設置過程中產生的。因此,為確保結構安全施工順利進行,應該保證上一排土釘或錨桿的結構強度達到設計要求,達到一定的承載能力后,才能進行下一步施工步驟。④在復合型土釘支護結構中,上部坡頂及下部坡底的土釘承受拉力較小,其中坡底的土釘承受的拉力大于坡頂部位的土釘拉力;位于中間部位的錨桿承受較大的拉力,其中第1排錨桿承受拉力最大;位于錨桿中間兩排土釘承受的拉力同樣較小,最大拉力位于靠近土釘與面層連接的一端。因此,在設計時為了更有效的控制地面變形,建議適當增加第1排錨桿的長度;錨桿中間土釘受力較小,不能充分發(fā)揮作用,建議適當減小位于錨桿中間的土釘長度,滿足其局部穩(wěn)定要求即可;下部坡底部位的土釘受力特點,建議采用短而密的土釘布置方式。

綜上:由FLAC3D軟件模擬的單根土釘拉力變化規(guī)律、每步開挖所引起的各排土釘及錨桿的拉力增長大小、土釘承受最大拉力的位置等與理論分析及實測數(shù)據(jù)變化規(guī)律基本吻合,可作為研究復合型土釘支護的有效工具。

2.6 支護結構的應力應變特征

支護結構的應力特征主要包括:最大主應力σ1、第二主應力 σ2和最小主應力 σ3、剪切應變等變化特征。在各步開挖結束后,在基坑底部靠近坡腳位置,受力狀況比較復雜。因此,在復合型土釘支護設計過程中,應對此位置給予足夠重視。模型中應力符合屈服準則的塑性區(qū)范圍如圖5所示,在坡面及坡腳處均存在過去曾經屈服單元,從而進一步說明坡腳處應力的復雜性及支護設計的重要性。必要時應考慮對坡腳采取提前加固以確保結構安全。

圖5 應力符合屈服準則的塑性區(qū)范圍

3 結 論

本文采用FLAC3D巖土數(shù)值分析軟件對土釘+錨桿復合型土釘支護結構進行三維立體模擬分析研究,考慮了基坑分步開挖對支護結構的影響,從而分析和總結了支護結構中各因素隨開挖進行所產生變化規(guī)律[6-7]。

(1)坑壁水平位移的分布規(guī)律。在基坑開挖支護前期,基坑變形主要是在應力重分布過程中產生的初步變形。隨開挖進行水平變形逐漸增大,靠近開挖面的土體水平位移較大?；幼畲笞冃伟l(fā)生在約距基坑底面1/3基坑開挖深度處。土釘對土體側向變形限制作用非常明顯,土釘錨桿及土體是一個彼此相互作用的有機整體,因此在結構分析過程中不應將各部分簡單的分開考慮。

(2)地面沉降分布規(guī)律。隨著基坑的開挖,基坑開挖影響范圍內應力重新分布,坑腳處應力集中分布?；拥乇淼呢Q直位移沿遠離坑壁方向呈曲線分布,每步開挖上一層土體的沉降量最大?；拥酌嫱馏w出現(xiàn)一定量的回彈。新形成的坑壁有向下移動的運動趨勢,對卸載后的基坑底面土體產生擠壓作用,從而土體形成擠壓隆起現(xiàn)象。

(3)土釘及錨桿的軸拉力分布。各排土釘承受的拉力在土釘長度方向呈相似曲線分布。每一步開挖均會引起土釘拉力的一定增量,土釘每段的增量近似相同,曲線整體分布趨勢相近。土釘中間靠近土釘與面層連接的一端承受的拉力最大,向兩端方向受力逐步遞減。

(4)復合型土釘支護的應力應變特征。被加固土體沿潛在優(yōu)勢滑裂面上存在著滑移的趨勢,對基坑底部坡腳位置的土體產生擠壓作用,使得基坑底面隆起變形,在基坑底部靠近坡腳位置,受力狀況比較復雜,在坡面及坡腳處均存在過去曾經屈服單元。

[1] 王永焱,林在貫.黃土的結構與物理力學性質[M].北京:科學出版社,1990.

[2] 錢鴻縉,王繼唐,等編.濕陷性黃土地基[M].北京:中國建筑工業(yè)出版社,1985.

[3] 鄭建國,張?zhí)K民.濕陷性黃土在增濕時的強度特性[J].水文地質工程地質,1989,(2):22-25.

[4] 高永貴,劉宗洲.西安地區(qū)深基坑邊坡支護設計探討[J].西安建筑科技大學學報,1998,30(3):261-264.

[5] Schlosser F,Unterreiner P.Soil nailing in france:research and practice[J].Transportation Research Record,1991,1330:72-79.

[6] Schlosser F.Behavior and Design of Soil Nailing[C]//Proc.of Symposium on Recent Developments in Ground Improvement Techiques,Bangkok,1982:380-420.

[7] 李厚恩,秦四清.預應力錨索復合土釘支護的現(xiàn)場測試研究[J].工程地質學報,2008,16(3):393-400.