真空預壓排水固結法有限元模擬及分析

■吳炎輝

（東山縣海通建設工程有限公司，漳州 363400）

1 前言

隨著我國經濟的快速發展，沿海地區的圍海造地工程不斷增加。但由于沿海地區分布著大量的軟土，其含水量高、抗剪強度低以及滲透系數小等特點使其成為在該區域建設工程時的一大難題。為使工程項目達到預期的安全標準和進度要求，就需要對軟土進行人工處理，真空預壓排水固結法是軟土加固處理中重要的方法之一。

真空預壓法自1952 年首次提出至今已有相當大的發展，在沿海地區的應用迅速增多，為了保證建設質量和建筑物的正常使用，對真空預壓法進行更深入的研究十分重要。

2 有限元模型的建立

有限元數值計算方法是目前運用最廣泛的數值計算方法之一，其中大型有限元分析軟件ABAQUS 具有很好的處理非線性和復雜邊界條件的能力，因此本文擬用ABAQUS 對真空預壓過程進行有限元建模和分析，并對模型中的參數進行敏感度分析。

2.1 真空預壓法簡介

如圖1 所示，真空預壓法是通過在軟土地基內部打入豎向排水系統，同時在軟土表面鋪設砂墊層作為水平排水系統，用不透氣的密封膜將表面密封使其與大氣隔絕，然后利用抽真空裝置對排水系統進行抽真空，使排水系統內形成一定的真空度并逐漸擴展至加固區土體內部并產生孔隙水壓力差，從而讓孔隙水向排水系統流動并最終排出，最終達到加固土體的目的。

圖1 真空預壓排水系統

2.2 砂井與砂墻的等效

砂井地基嚴格意義上是屬于三維固結問題，應該采用三維有限元模型進行計算，但由于三維有限元分析計算量太大不利于應用，因此將三維問題轉換為二維問題進行計算分析的方法被廣泛采用。

如圖2 所示，將三維和二維問題進行等效的原則是保持等效前后不同時刻具有相同的固結度，其主要方法是通過修改滲透系數達到的，等效后與等效前的滲透系數之比為：

其中：khp為二維模型中土體的等效滲透系數；n 為涂抹區直徑與砂井直徑之比；s 為砂井有效排水直徑與砂井直徑之比。

圖2 三維與二維問題等效示意圖

2.3 模型參數的選取

該模型在大型通用有限元軟件ABAQUS 里建立，采用平面應變有限元模型，其尺寸如圖3 所示，由于的二維模型是左右對稱的，為了便于分析，該模型取整個二維邊界的右半邊進行建模，土體高10m，寬15m，加固區寬度5m，豎向排水板深度5m，模型左右兩側采用限制水平位移的邊界條件以模擬對稱及土體無限延伸的情況，底部采用限制水平和豎直位移的邊界條件用以模擬半空間無限體的條件，固結時間設置為約694h。土體及排水板參數如表1 所示：

圖3 有限元模型尺寸及邊界條件示意圖

表1 有限元模型計算參數

3 有限元模型計算結果分析

3.1 位移計算結果分析

井阻為真空度在排水板內傳遞的衰減率，該模型采用的井阻參數為16kPa/m，其排水板內真空度分布圖如圖4 所示。

圖4 排水板真空度分布圖

圖5 和圖6 分別是預壓結束時得到的水平位移和豎向位移云圖，從圖中可以看出，土體沉降主要集中在加固區中心，并隨著與中心距離的增大而逐漸減少，其中心點沉降約為40mm，加固區邊緣處的沉降約為35mm，而距離中心10m 和15m 處的沉降分別為10mm 和8mm，不足中心點沉降的1/3。

水平位移方面，加固區中心點處的水平位移幾乎為0，加固區邊緣的水平位移約為18mm，而最大水平位移發生在距離加固中心約6.5m 處，其水平位移約為22cm，方向指向加固區。

圖5 預壓結束時豎向位移云圖

圖6 預壓結束時水平向位移云圖

不同深度處橫向水平位移分布情況如圖7 所示，雖然位移量有所差異，但不同深度處的最大水平位移都在離中心約6.5m 的位置。

圖7 不同深度處水平向位移的橫向分布

3.2 超孔隙水壓力計算結果分析

超孔隙水壓力在模型中的分布情況如圖8 所示。從圖中可以看出，真空荷載使土體內產生負超孔隙水壓力，并形成一種以加固區為中心向外逐漸擴散的形式。

圖8 預壓結束時超孔隙水壓力云圖

圖9 所示的為距離加固區不同距離的超孔隙水壓力隨深度的分布。從圖中可以看出，越遠離加固區，其超孔隙水壓力的下降得越快，如圖10 所示，土體表面、2m 深和4m 深處超孔隙水壓力隨水平距離的變化率分別為14.5kPa/m、8.9kPa/m 和1.8kPa/m。

圖9 超孔隙水壓力隨深度的分布

圖10 超孔隙水壓力沿水平方向的分布

3.3 土體固結效率對井阻的敏感度分析

井阻大小受排水板材料、抽真空器材及加固區土體情況等諸多因素的影響。圖11、圖12 和圖13 為不同井阻情況下加固區中心處土體不同深度的沉降情況。

圖11 加固區中心土體表面沉降值

圖12 加固區中心3m 深處沉降值

圖13 加固區中心5m 深處沉降值

在加固區中心點處，不同深度的沉降差隨井阻大小的變化而改變。隨著深度的增加，沉降差逐漸增大，如圖14所示。井阻16kPa/m 與5kPa/m 的沉降差在表面可達到0.041m，其斜率為0.0037，在3m 和5m 深度處沉降差分別為0.033m 和0.019m，其斜率分別為0.003 和0.0017，可以看出，加固區較深位置的沉降量受井阻影響較小。

圖14 中心點不同深度處沉降差

4 結論

（1）土體最大豎向位移處于加固區中心，在加固區以外的地方，其豎向位移大幅減小。而最大水平位移發生在靠近加固區邊緣的地方。

（2）加固區土體內的真空度可以維持在與排水板內真空度相當的水平，在遠離加固區的地方，由于沒有采用密封措施，其真空度則難以保持在一個較高的水平。

（3）井阻對土體固結效率有顯著影響，如果能使井阻維持在一個較小的數值，對土體效果能起到非常積極的作用。

[1]費康,張建偉.ABAQUS 在巖土工程中的應用[M].北京市:水利水電出版社,2010.

[2]Li B,Wu S F,Chu J,et al.Evaluation of Two Vacuum Preloading Techniques Using Model Tests[C].ASCE,2011.636-645.

[3]陳環.真空預壓法機理研究十年[J].港口工程,1991,(4):17-25.

[4]Indraratna B,Sathananthan I,Rujikiatkamjorn C,et al.Analytical and numerical modeling of soft soil stabilized by prefabricated vertical drains incorporating vacuum preloading[J].International Journal of Geomechanics,2005,5(2):114-124.

[5]張澤鵬,李約俊,馮淦清,等.塑料排水板在真空預壓加固軟基中的作用[J].廣州大學學報(自然科學版),2002,1(2):68-71.

[6]沈珠江,陸舜英.軟土地基真空排水預壓的固結變形分析[J].巖土工程學報,1986,8(3):7-15.