防波堤施工過程中地基變形及穩定性數值分析

周建鋒 耿瑞巖 丁頁嶺

(1.無錫市地震局，江蘇 無錫 214000； 2.連云港市地震數據處理中心，江蘇 連云港 222002； 3.江蘇省地震局，江蘇 南京 210000)

防波堤施工過程中地基變形及穩定性數值分析

周建鋒1耿瑞巖2丁頁嶺3

(1.無錫市地震局，江蘇 無錫 214000； 2.連云港市地震數據處理中心，江蘇 連云港 222002； 3.江蘇省地震局，江蘇 南京 210000)

以堤身填筑過程引起的地基沉降以及對邊坡穩定影響為研究對象，使用有限元軟件建立流固耦合數值分析模型，模擬了某防波堤工程堤身填筑過程，對沉降規律以及邊坡穩定進行了計算分析，為現場施工提供參考。

防波堤，地基沉降，地基水平位移，邊坡穩定性

1 研究背景

防波堤施工過程中會引起一系列的地基變形，同時，堤身穩定會在施工過程中發生變化，一些學者據此展開了一些研究。對于地基變形，通常分為在豎直方向上的沉降和水平方向的位移，例如劉銳等[1]通過研究監測地基沉降的沉降盤，對某防波堤工程進行動態的沉降監測，分析了各級加載階段的沉降變化規律；郭少龍等[2]根據現場實測數據，通過對地基沉降的觀測，分析了基于MIDAS/GTS軟件建立的地基加固效果研究的數值模型的可靠性；對于防波堤穩定性計算，之前研究者們常采用極限平衡法和強度折減法等，例如路衛衛[3]采用有限元強度折減法進行防波堤穩定計算，分析了孔隙水壓力及有效應力隨時間的變化以及對防波堤穩定性的影響；宣廬峻[4]結合國外某軟土地質條件下的防波堤工程案例,采用基于Biot固結理論的有限元法模擬整個防波堤建造過程的穩定與變形分析,并與常規計算方法的結果進行了比較。因此，建立有限元模型，對施工過程中地基變形和整體穩定性的計算，對保障施工安全有至關重要的作用。

2 有限元模型

2.1 模型建立

在荷載作用下，土體的超孔隙水壓力逐漸消散，內部含水緩慢滲出，體積逐漸減小，有效應力在一定范圍內逐漸增大，這種現象稱為土的固結。工程實踐中，軟土的固結沉降主要發生在超孔隙水壓力完全消散之前，但在固結完成(超靜水壓力完全消散)之后，軟土仍產生一定量緩慢的沉降變形，這就是土的蠕變作用。因此，完整的土體沉降應該是軟土固結沉降與軟土蠕變共同作用的結果。固結基本方程基于比奧(Biot)理論，滲流問題采用達西(Darcy)定理，基于小應變理論，假設土體骨架彈性變形[5,6]。根據太沙基(Terzaghi)原理，土體中的應力分為有效應力和孔隙壓力：

σ=σ′+m×(Psteady+Pexcess)

(1)

其中，σ=(σxx,σyy,σzz,σxy,σyz,σzx)T;m=(1,1,1,0,0,0)T,m包含單位正應力分量和零剪應力分量的矢量；σ為總應力矢量；σ′為有效應力；Pexcess為超孔隙水壓力；固結過程最終的穩態解表示為Psteady。

2.2 三維簡化

由于塑料排水板的形狀比較特殊，在理論分析過程中存在一定困難，一般將塑料排水板換算為等周長的圓截面排水體——等效砂井[7]。

塑料排水板在平面上按等邊三角形布置，其有效排水范圍為正六邊形，該有效范圍內的水通過位于其中的排水體排出。在進行固結分析時，將每個排水體的影響范圍轉換為一個等面積的圓[8]，以方便計算。

砂井地基的工作狀態是三維的。由于砂井布置比較密集，需劃分的單元數目很多，導致模型計算量非常大，降低了計算效率，因此把三維砂井轉換為二維來研究。轉換的前提是保證前后固結度與平均孔壓不變，首先增大砂井的間距，以減少單元數目，從而減少計算量；然后用沿著縱向連續分布的砂墻替換有一定間隔的砂井，即把原來的砂井地基轉變成了砂墻地基[9]。根據Barron理論，考慮井阻和涂抹作用的影響，只要調整滲透系數即可保證砂井與砂墻等效。

經過兩次簡化后，三維的塑料排水板將用一排排二維的砂墻地基代替，這極大提高了運算效率。

2.3 數值模型

考慮到防波堤整體長度達5 km左右，各斷面為典型的平面應變問題，因此在防波堤上取一典型斷面進行平面應變建模分析。所建立的典型數值分析模型如圖1所示。

根據施工進度安排，模擬了在設計低水位(-0.10 m)條件下研究斷面的以下幾個過程：施打塑料排水板、鋪設土工格柵、鋪設二片石、拋填堤心石至第一加載分界線，拋填堤心石和墊層石至第二加載分界線，拋填堤心石、墊層石和護面塊體至堤頂。

3 結果分析

圖2～圖4為計算結果云圖。按照施工過程，進行各個填筑階段的總體豎直沉降、水平位移、超孔隙水壓力計算，同時，對坡腳穩定性進行分析。經長期排水固結，達到最終沉降基本穩定時(最大超靜孔壓小于10 kPa)，累計沉降量約為1 040 mm，工后沉降184.3 mm。

3.1 第一階段堤身填筑

在排水板施工完成并鋪設土工格柵后，拋填0.5 m厚的二片石。二片石鋪設后，維持穩定一段時間使地基土體排水固結沉降，共歷時6 d，此間最大沉降為9.3 mm，沉降速率約為1.55 mm/d，最大超靜孔壓為5.28 kPa。

第一階段堤身填筑至標高-5.9 m，本級填筑高度2.1 m，累計填筑高度2.6 m，填筑后維持穩定40 d，使土體排水固結。此間本級沉降量為117.9 mm，最大累計沉降量約為127.2 mm，本級沉降速率為2.87 mm/d，累計平均沉降速率為2.71 mm/d，剛填筑至指定標高時超靜孔壓約為23.2 kPa，堤身坡腳穩定安全系數為6.0，固結沉降后最大超靜孔壓降低為17.4 kPa，堤身填筑第一階段累計水平位移分布如圖2所示。

3.2 第二階段堤身填筑

第二階段堤身填筑至標高-2.26 m，本級填筑高度3.64 m，累計填筑高度6.24 m，填筑后維持穩定100 d，使土體排水固結。此間本級沉降量為316.2 mm，最大累計沉降量約為443.3 mm，本級沉降速率為3.13 mm/d，累計平均沉降速率為2.99 mm/d，剛填筑至指定標高時超靜孔壓約為51.8 kPa，堤身坡腳穩定安全系數為2.59，固結沉降后最大超靜孔壓降低為28.7 kPa，堤身填筑第二階段累計水平位移分布如圖3所示。

3.3 第三階段堤身填筑

第三階段堤身填筑至標高+5 m，本級填筑高度7.26 m，累計填筑高度13.5 m，填筑后維持穩定使土體排水固結。按30 d完成本級填筑，完成時本級沉降量為412.4 mm，最大累計沉降量約855.7 mm，本級沉降速率13.7 mm/d，累計平均沉降速率4.81 mm/d，

剛填筑至指定標高時超靜孔壓約為52.6 kPa，堤身坡腳穩定安全系數1.35。經一年排水固結沉降后，最大超靜孔壓降低為36.5 kPa，最大累計沉降為922.1 mm，一年工后沉降為66.4 mm。堤身填筑第三階段累計水平位移分布如圖4所示。

4 結語

建立有限元流固耦合模型，對某防波堤施工過程進行了數值模擬計算，選取典型斷面進行了固結沉降、水平位移、超孔隙水壓力和整體穩定的計算和分析，主要結論如下：

1)隨著堤身的填筑，地基沉降逐漸增大，沉降速率也逐步增大，最大沉降速率出現在第三級填筑階段，計算斷面的最大沉降速率為1.37 cm/d，研究斷面的最大沉降為1.04 m。

2)由于施打塑料排水板，防波堤在填筑施工完成時固結度在80%左右，在施工完成一年后固結度在86%左右，達到最終穩定(最大超靜孔壓小于10 kPa)約需要2 548 d。

3)隨防波堤填筑高度的增大，堤身整體穩定性安全系數逐漸變小，在設計低水位時第三級填筑完成后堤身整體穩定安全系數最小，其值約為1.3，滿足規范要求(1.3～1.5)，防波堤整體穩定基本滿足要求，但在實際填筑過程中需密切關注堤身變形，控制填筑速度。

4)隨防波堤填筑高度的增大，堤身整體穩定性安全系數逐漸變小，在第三級填筑完成后堤身整體穩定安全系數最小，約為1.3，滿足規范要求。因此，防波堤整體穩定基本滿足要求，但在實際填筑過程中需密切關注堤身變形，控制填筑速度。

[1] 劉 銳,徐振坤,張秀勇，等.東營港北防波堤工程沉降位移分析[J].港口科技,2016(8):9-14.

[2] 郭少龍,蔣學煉,鹿 群，等.防波堤工程軟土地基排水固結處理的數值分析[J].土工基礎,2012,26(6)：88-91.

[3] 路衛衛.PLAXIS在防波堤穩定計算中的應用[J].水運工程,2014(4):40-42,62.

[4] 宣廬峻.國外某軟土地質防波堤工程的地基穩定分析[J].水運工程,2013(9):108-115.

[5] 趙 杰,程 翔,王桂萱.某大型填海造地工程人工島護岸結構固結沉降及蠕變數值分析[J].水運工程,2013(5):150-156.

[6] 王桂萱,易 劍,趙 杰.某大型填海造地工程的一維固結沉降數值分析[J].中國港灣建設,2013(6):10-14.

[7] 翟聚云,衛國祥,馬明江.土木工程專業巖土工程方向畢業設計選題研究[J].高等建筑教育,2014(3):107-109.

[8] 龍桂根,徐光輝,黃芝云,等.裸露山體困難立地綠化關鍵技術研究[J].江西林業科技,2011(3):27-32.

[9] 周冰峰,王培俊,蘇薇國，等.巖體平面應變模型試驗臺設計與有限元分析[J].工程機械,2014(3):37-40,52.

Numerical analysis of foundation deformation and stability in breakwater construction

Zhou Jianfeng1Geng Ruiyan2Ding Yeling3

(1.WuxiSeismologicalBureau,Wuxi214000,China; 2.LianyungangSeismicDataProcessCenter,Lianyungang222002,China; 3.JiangsuSeismologicalBureau,Nanjing210000,China)

Taking the foundation settlement caused by embankment construction process and the influence to slope stability as the research object, this dissertation simulates one certain breakwater construction processes and numerically analyses the foundation deformation and slope stability by means of finite element numerical analysis model of fluid structure interaction, so as to provide reference for site construction.

breakwater, foundation settlement, foundation water level, slope stability

1009-6825(2017)20-0062-02

2017-04-19

周建鋒(1978- )，男，工程師; 耿瑞巖(1978- )，男，助理工程師; 丁頁嶺(1981- )，女，工程師

TU441.6

A