基于ABAQUS的真空預壓插板深度的研究

張澤波，朱 平，左 軍

（1.天津大學 建筑工程學院，天津 300072；2.天津市北洋水運水利勘察設計研究院有限公司，天津300452）

基于ABAQUS的真空預壓插板深度的研究

張澤波1，朱 平1，左 軍2

（1.天津大學 建筑工程學院，天津 300072；2.天津市北洋水運水利勘察設計研究院有限公司，天津300452）

通過 ABAQUS大型有限元軟件模擬了實際工程中的部分加固過程，驗證了此類工程數值模擬的可行性。ABAQUS有限元軟件計算的地基沉降規律和沉降值與實測結果非常接近，其中不一致部分為簡化計算引起的誤差。通過模擬不同插板深度的加固，對比了對應深度的加固效果，分析了底層的固結沉降效果。得出真空預壓插板深度的變化會對工程加固效果產生影響，其中對底層土體影響較為顯著；排水板打設穿過目標軟基底部一定深度可有效提高底層土體固結度，進而提高軟基整體加固效果。

真空預壓；插板深度；ABAQUS；固結度

引 言

天津濱海新區發展速度逐年加快，隨著其快速發展，也出現了很多的現實問題，其中很重要的一個問題就是土地缺乏。為了緩解這一問題，此地區進行了大面積的圍海造陸工程。天津濱海新區的圍海造陸工程所用吹填土來自港池航道的疏浚淤泥或淺海底部的泥沙，這些泥沙經過水力作用吹填到指定區域形成陸域場地。由于這些泥沙的含水率很高，其范圍一般為 80 %～200 %之間［1］，因此其承載力極地，強度指標很不理想。為了能快速有效地加固此類軟土地基，一般的加固方法是無法完成的，而真空預壓法對于此類地基的處理效果能達到設計要求。

真空預壓法是利用抽真空設備產生的真空負壓作為預壓荷載使土中水在壓差下排出從而使軟土地基得到加固的一種處理方法。真空預壓系統由3部分組成，包括：真空壓力系統、密封系統和排水系統。排水系統在整個真空預壓軟基處理工程中的變化會對處理效果產生顯著影響，值得學者和技術人員進行深入研究。

近年來很多學者對于真空預壓法進行了相關研究，并取得了一定成果。施建勇、雷國輝等［2］研究了軟土的壓縮特性，提出了在真空預壓處理過程中軟基沉降的計算方法；陳平山、莫海鴻［3］應用等效原理，將塑料排水板等效成砂井地基，對真空預壓軟基處理中的固結度進行了研究，提出了更有效的簡化方法；閆澍旺等［4，5］提出了真空預壓軟基加固中固結度的估算方法，并對真空預壓有效加固深度進行了研究；Chai J.C.等［6］通過實驗與數值模擬的研究，建立了軟基處理中真空預壓法最優插板深度的解析方程，此方程決定了軟基的最大整體固結度和相應的沉降量；Rujikiatkamjion C.和Indraratna B.［7］經過實驗論證和數值計算提出了真空預壓中軟基固結度的解析解，對土體的豎向固結和水平向固結都進行了研究。以往對真空預壓法的研究都是針對其整體的固結度或沉降量進行的研究，對單獨各層或軟基底層的固結沉降研究很少。本文先通過ABAQUS大型有限元軟件模擬了實際工程中的部分加固過程，驗證了此類工程數值模擬的可行性，然后模擬了不同插板深度的加固，對比了對應深度的加固效果，著重分析了底層的固結沉降效果，得出了提高處理效果、增加數據可靠度的工程插板方式，可為類似軟基加固工程提供借鑒。

1 工程實例有限元分析與驗證

1.1 工程概況

天津臨港經濟區南部區域緯五路一期工程地基處理工程位于臨港中港池以北，地基處理總面積57.3萬m2，其中道路加固區面積約為49萬m2，共分為22個處理區，其中1～4區和6～22區的地基處理方式為真空預壓法。根據設計圖紙要求，1～4區在原泥面鋪設一層編織布，兩層竹笆和一層短纖土工布，吹填砂性土。吹填后場地應整平，打設塑料排水板前后均應測量砂墊層表面高程，以便推算打板期間沉降量，排水板打設底高程為-12.5 m。6～22區淺層真空預壓在原泥面鋪設一層編織布，一層短纖土工布，作為工作墊層；深層真空預壓在已鋪設完畢砂墊層的淺層預處理面打設排水板、鋪膜、吹砂等。本文選取3區作為研究對象，對3區用ABAQUS有限元軟件來模擬分析。

3區面積為90 m×234.5 m，水平排水系統的砂墊層為1 400mm厚的粉細砂，排水濾管布置于其中，豎向排水系統為正方形布置的塑料排水板，其間距為850mm。工程應滿足設計提出的以下卸載要求：1）按實測沉降曲線推算的固結度大于85 %；2）連續5天實測地表平均沉降速率不大于3.0mm/d。有限元分析對應的場地模型示意如圖1所示。

圖1 場地模型示意

3區為人工吹填土，加固區土層分布情況為：①淤泥質土層，呈灰色，高壓縮性；②粉質粘土，呈灰色，中間混有粉土斑；③淤泥質粘土，呈褐灰色，壓縮性較高；④淤泥質粘土，呈灰色，飽和，高壓縮性；⑤淤泥質粘土和粉質粘土疊加，呈灰色，飽和度較高，高壓縮性。在此軟基中埋入7個磁環，用來測量各土層沉降量。各土層的對應的土性參數如表1所示。

表1 各層土性參數

1.2 固結理論分析

真空預壓處理軟土地基的過程實際上就是軟基在真空負壓荷載下的排水固結過程。而土體的固結沉降及其孔隙壓力消散的計算需要一定的固結理論的支撐。太沙基一維固結理論忽略了變形綜合條件對固結總壓力的影響，此理論是在固結變形只發生在豎直方向的假定下成立的，因此這個理論不符合真空預壓的真實固結過程。巴隆的軸對稱固結理論由于計算太過復雜，計算所用到的各種參數很難獲取，不容易推廣到實際工程中去應用。比奧固結理論較真實的反映了土體固結變形與孔隙壓力的相互關系，同時考慮了垂直沉降和水平位移，其所用到的各種計算參數條件也較容易獲取，因此此理論在實際工程中比較實用［8］。根據真空預壓原理和平面應變比奧固結理論，比奧固結微分方程如下：

1.3 排水板等效原理

現今工程所使用的塑料排水板與傳統的砂井相比，雖然其材料不同，但是在理論計算中其加固原理和理論方法是相同的，因此塑料排水板在模擬計算時可以等效為砂井。又由于在實際的真空預壓工程中軟土中的排水板布置密集，可以把排水板在縱向的布置形式看成是一個平面，將三維固結問題轉化為較為簡單的平面固結問題，這樣就大大簡化了網格，減小了計算量。將三維的排水板轉化為平面形式即二維排水砂墻，其滲透系數調整公式如下所示［9］：

其中：kxa、kza分別為軟土地基實測的水平和豎直滲透系數；kxp、kzp分別為有限元計算時所用滲透系數；Dx、Dz分別為滲透系數調整系數；μ為土的泊松比；n為砂井有效排水半徑和砂井半徑之比。

1.4 ABAQUS有限元模型

由于3區場地兩邊的非加固區地形相似，又處于同一處理區，因此土質參數相近。3區場地按中心軸在橫向呈對稱狀，其橫向寬度為90 m，因此按中心線取 45 m作為有限元模型區域的加固區。又由于真空預壓加固過程中對加固區以外較遠處的沉降和水平位移影響較小，因此有限元模型的影響區只去45 m進行計算，也就是從加固區邊緣向外延伸45 m。有限元計算模型的網格劃分如圖2所示。

圖2 有限元計算模型網格劃分示意

由于加固區在鋪設砂墊層以后的地面高程為5 m，塑料排水板打設深度底部高程為-12.5 m，考慮到真空預壓可以影響到排水板底部以下的土體，再結合土質條件，有限元模型的計算深度取30 m。邊界條件的選擇為：模型底部為雙向約束；兩個側邊界為水平約束；表面無約束，為自由變形面；影響區表面孔壓邊界條件為 0；加固區排水板單元中的孔壓為線性施加，頂部真空度為85 kPa，沿程損失為3 kPa/m。土體采用摩爾—庫倫本構模型。荷載步為：真空壓力為瞬時一次施加，恒載時間為工程記錄的實際的恒載時間。網格劃分采用的單元類型為8節點的平面應變孔壓單元，四邊形結構化網格。由于加固區需要重點分析，網格劃分較密，影響區網格較為稀疏。

1.5 結果分析

有限元模型計算結果中的地基沉降位移如圖3所示。

圖3 地基沉降位移

由圖3可以看出，加固區沉降規律為：加固區中心以及向外的一定的范圍內軟基沉降量較大；加固區表層土體的沉降量較大，隨著深度的增加沉降量越來越小；真空預壓可以影響到塑料排水板一下一定深度范圍內的土體，但是其沉降量很微小。真空預壓可以影響到加固區以外的土體，其水平方向有向加固區收縮的趨勢，在一定的范圍內也會發生沉降。這些規律都與真空預壓理論相一致。

圖4為加固區中心土體的有限元模型計算沉降曲線和工程實測沉降曲線對比圖，此處的沉降值為真空加載下的地基沉降值，不包括插板期間的沉降。我們可以看出有限元計算的地基沉降規律和沉降值與實測值非常接近，其中不一致部分為簡化計算引起的誤差。有限元可很有效的模擬實際的真空預壓處理工程。

圖4 計算沉降曲線和工程實測沉降曲線對比

圖5 磁環沉降數值

圖5為實測加固中心區各土層所安裝磁環的沉降量數值圖，從圖中可以看出表層軟土的沉降量最大，底層的沉降量很小，其固結度較低。從前面有限元模型計算沉降曲線和工程實測沉降曲線對比圖可知數值計算值與實測值非常接近，因此我們可以用數值計算值來分析和預測實際的處理工程。實際工程的各土層的分層沉降量也可用數值計算值來模擬。

2 不同插板深度土體分層固結度分析

在第二節有限元模型分析結果的基礎上，本節進行了不同插板深度的模型分析。第二節所述軟基加固工程的實際加固深度的土底高程為-13 m，而其所打設塑料排水板底部高程為-12.5 m，這是因為真空預壓的影響深度可以達到排水板以下一定范圍內。但是各土層的分層沉降量相差很大，底部土層的沉降很小，因此其固結度與土體的整體固結度相差很大。

在第二節有限元模型的基礎上，變換排水板深度建立三個模型A、B、C。對應排水板底部高程分別為-12.5 m、-13.0 m、-13.5 m。也即A的排水板打設深度在所加固軟土底部以上0.5 m，B的排水板正好打到所處理軟土底部，C的排水板打穿所加固軟土層到其底部以下0.5 m。整理模型計算結果，根據A、B、C三個模型所得的各土層分層沉降量，然后計算A、B、C各土層的卸載時的固結度。其卸載時的固結度數值如圖6所示。

圖6 各土層固結度數值

由圖6可以看出隨著插板深度的增加上不和中部土層的固結度變化不明顯，而底部土層的固結度有明顯的增大。這是由于雖然真空預壓的影響深度可以達到排水板一下一定的范圍，但是由于所處理地基為淤泥質粘土，其滲透系數很小，在排水板未到達的深度土中孔隙水無法有效的排出，致使其沉降很小，固結度較低。在實際工程中容易出現底層軟基未能達到整體處理要求，其強度指標的提高不明顯。當排水板打穿目標軟基時，底層土中孔隙水能更有效的排出，固結度接近軟基整體的處理要求，處理效果更好。我們還可以推理得出，在達到同樣的處理效果的前提下，當排水板打設深度穿透目標軟基底部時，處理工期會縮短。

3 結 論

1）ABAQUS有限元軟件計算的地基沉降規律和沉降值與實測結果非常接近，其中不一致部分為簡化計算引起的誤差。有限元可很有效的模擬實際的真空預壓處理工程。也可對實際工程進行預測。

2）真空預壓可以影響到加固區以外的土體，其水平方向有向加固區收縮的趨勢，在一定的范圍內也會發生沉降。

3）真空預壓插板深度的變化會對工程加固效果產生影響，其中對底層土體影響較為顯著。排水板打設穿過目標軟基底部一定深度可有效提高底層土體固結度，進而提高軟基整體加固效果。

［1］董志良，張功新，周琦，等.天津濱海新區吹填造陸淺層超軟土加固技術研發及應用［J］.巖石力學與工程學報，2011，30（5）：1073-1080.

［2］施建勇，雷國輝，艾英缽，等.關于真空預壓沉降計算的研究［J］.巖石力學，2006，27（3）：365-368.

［3］陳平山，莫海鴻，等.真空預壓加固地基等效固結度簡化計算［J］.華南理工大學學報（自然科學版），2008，36（1）：139-144.

［4］Chu J，Yan S W.Estimation of degree of consolidation for vacuum preloading projects［J］.Int Journal of Geomechanics，ASCE，2005，5（2）：158-165.

［5］Chu J，Yan S W.Effective depth of vacuum preloading［J］.Lowland Technology International，2006，8（2）：1-8.

［6］Chai J C，Carter J P，Hayashi S.Vacuum consolidation and its combination with embankment loading［J］.Canadian Geotechnical Journal，2006，43（5）：985-996.

［7］Rujikiatkamjorn C，Indraratna B.Analytical solutions and design curves for vacuum-assisted consolidation with both vertical and horizontal drainage［J］.Canadian Geotechnical Journal，2007，44：188-200.

［8］孫立強.超軟吹填土地基真空預壓理論及模型試驗的研究［D］.天津：天津大學，2010.

［9］趙維炳，陳永輝，龔友平.平面應變有限元分析中砂井的處理方法［J］.水利學報，2004，（01）：16-20.

Research on Depth of Vacuum Preloading Plug-board Based on ABAQUS

Zhang Zebo1，Zhu Ping1，Zuo Jun2
（1.Institute of Civil Engineering，Tianjin University，Tianjin 300072，China； 2.Tianjin Beiyang Water Transport and Water Conservancy Survey and Design Institute Co.，Ltd.，Tianjin 300452，China）

Large-scale finite element software ABAQUS has been used to simulate partial reinforcing process of real projects，which proves the feasibility of numerical simulation.The settlement regularity and sedimentation value obtained by using ABAQUS finite element software basically accord with the measured results，and the inconformity of data comes from the deviation caused by simplified calculation.By simulating the reinforcement process based on different depths of plug-board，the reinforcement effects are compared each other，and the consolidation settlement effects of bottom layer are analyzed.The research results show that the depth variation of vacuum preloading plug-board affects the reinforcing result，especially the soil body at bottom.To drive drain board into soft base to a certain extent can effectively improve the consolidation degree of soil body at bottom，and enhance the integral reinforcing effect of soft base.

vacuum preloading； depth of plug-board； ABAQUS； degree of consolidation

TU472.3+3

1004-9592（2016）05-0097-05

10.16403/j.cnki.ggjs20160524

2015-12-06

張澤波（1990-），男，碩士研究生，研究方向為軟土地基加固處理。