馬哲超，盛佳珺

（中交第三航務(wù)工程勘察設(shè)計(jì)院有限公司，上海200032）

?

桶式結(jié)構(gòu)穩(wěn)定性三維有限元數(shù)值分析

馬哲超，盛佳珺

（中交第三航務(wù)工程勘察設(shè)計(jì)院有限公司，上海200032）

使用三維有限元數(shù)值分析方法，結(jié)合土體彈塑性模型，計(jì)算分析桶式結(jié)構(gòu)在填土和使用荷載作用下，結(jié)構(gòu)的位移和穩(wěn)定性，通過對(duì)比不同工況和不同參數(shù)的計(jì)算結(jié)果，得出影響桶式結(jié)構(gòu)穩(wěn)定性的關(guān)鍵因素，即桶式結(jié)構(gòu)在13 m以下的填土高度，結(jié)構(gòu)的穩(wěn)定性在施工期有較大的安全儲(chǔ)備，而使用期使用荷載的大小、與結(jié)構(gòu)的距離是影響整體穩(wěn)定性的主要因素。

桶式基礎(chǔ)；軟土地基；位移；穩(wěn)定性

0　引言

對(duì)于桶式基礎(chǔ)結(jié)構(gòu)，目前設(shè)計(jì)上沒有專門的規(guī)范，因此一些技術(shù)參數(shù)必須通過深入系統(tǒng)的研究來獲得。桶式結(jié)構(gòu)設(shè)計(jì)計(jì)算核心是荷載作用下結(jié)構(gòu)的整體變形與受力分析，即結(jié)構(gòu)與地基的共同作用問題。要搞清楚結(jié)構(gòu)的變形和整體穩(wěn)定性，最重要的是對(duì)地基土的變形與受力情況的正確反映，而這對(duì)于數(shù)值分析往往是最困難的[1-4]。因此，本文采用彈性模型定性分析桶式基礎(chǔ)結(jié)構(gòu)的穩(wěn)定性，結(jié)合現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)反演分析結(jié)構(gòu)變形，從中尋找設(shè)計(jì)的控制參數(shù)，為此新型桶式基礎(chǔ)結(jié)構(gòu)設(shè)計(jì)、優(yōu)化提供可靠的技術(shù)依據(jù)。

1　計(jì)算原理

徐圩港區(qū)位于連云港區(qū)南翼，埒子口以西至小丁港之間海岸，隸屬連云區(qū)，目前為開敞海岸，防波堤工程建設(shè)是徐圩港區(qū)起步建設(shè)的前提條件。港區(qū)設(shè)計(jì)高水位5.41 m，設(shè)計(jì)低水位0.47 m，極端高水位6.56 m，極端低水位-0.68 m。

結(jié)構(gòu)設(shè)計(jì)方案由下部鋼筋混凝土橢圓腔體結(jié)構(gòu)件和上部擋土墻組成。每1組由1個(gè)基礎(chǔ)桶結(jié)構(gòu)件和上部擋土墻結(jié)構(gòu)組成；基礎(chǔ)桶體呈橢圓形，桶內(nèi)通過隔板劃分隔倉；上部擋土結(jié)構(gòu)坐落在基礎(chǔ)桶底板上，通過底板上的圈梁連接。見圖1。

本文整體穩(wěn)定分析采用有限元強(qiáng)度折減法，其原理是通過不斷減低邊坡巖土體抗剪強(qiáng)度參數(shù)直至達(dá)到極限破壞狀態(tài)為止，程序自動(dòng)根據(jù)彈塑性有限元計(jì)算結(jié)果得到滑動(dòng)破壞面的巖土體抗剪強(qiáng)度參數(shù)的降低倍數(shù)，即為邊坡強(qiáng)度安全系數(shù)。

對(duì)于采用摩爾-庫倫強(qiáng)度準(zhǔn)則模擬的材料，則有：

這種強(qiáng)度折減安全系數(shù)的定義與邊坡穩(wěn)定分析的極限平衡安全系數(shù)的定義形式是一致的。

2　桶式基礎(chǔ)防波堤三維有限元框架的建立

2.1基本假設(shè)

桶式結(jié)構(gòu)計(jì)算屬于典型的土與結(jié)構(gòu)相互作用問題，需要采用有限元增量法計(jì)算。在建立有限元分析模型時(shí)做如下假設(shè)：

1）不考慮施工對(duì)土體性質(zhì)的影響；

2）同一層土體視為均質(zhì)各向同性材料；

3）數(shù)值計(jì)算不考慮地震等動(dòng)力荷載的影響；4）數(shù)值計(jì)算沒有考慮土體的流變特性。

2.2計(jì)算參數(shù)的確定

數(shù)值模擬中土體采用摩爾-庫倫彈塑性模型；上層擋土墻結(jié)構(gòu)，下層桶頂蓋板以及下層桶壁采用各向同性彈性材料的板單元來模擬。為了模擬結(jié)構(gòu)與土之間的相互作用及荷載傳遞特性，在結(jié)構(gòu)表面設(shè)置接觸面單元，接觸面單元參數(shù)根據(jù)鄰近土層參數(shù)取值。

對(duì)于任何數(shù)值分析，準(zhǔn)確的計(jì)算參數(shù)是保證計(jì)算合理性的基礎(chǔ)。為了進(jìn)行連云港港徐圩港區(qū)桶式基礎(chǔ)結(jié)構(gòu)與土體相互作用分析，進(jìn)行了大量的室內(nèi)與現(xiàn)場(chǎng)土工試驗(yàn)，得到了計(jì)算所需的各種參數(shù)見表1。桶式結(jié)構(gòu)主要參數(shù)為彈性模量3.25× 107MPa，水中重度14 kN/m3，泊松比0.2。

表1　土體主要參數(shù)表Table 1　Main parameters of soil

地基土體的初始側(cè)壓力系數(shù)K0是一個(gè)非常重要的參數(shù)，它一方面決定了桶式基礎(chǔ)結(jié)構(gòu)下桶與地基之間摩擦力的大小，另一方面作為土體的水平抗力對(duì)結(jié)構(gòu)的變形與整體穩(wěn)定性起決定性的作用。本計(jì)算過程中，K0通過下式計(jì)算：

式中：ψ為土體的內(nèi)摩擦角。對(duì)于淤泥層，ψ=2.3°，K0=0.96；對(duì)于粉質(zhì)黏土層，ψ=11.5°，K0=0.8。

地基土體與防波堤上下桶結(jié)構(gòu)接觸面的摩擦系數(shù)通過離心模型試驗(yàn)獲得，根據(jù)離心模型試驗(yàn)結(jié)果，得到結(jié)構(gòu)與土體之間摩擦系數(shù)的范圍。離心試驗(yàn)得到的摩擦系數(shù)為0.12～0.14，計(jì)算采用的摩擦系數(shù)取0.115。

2.3計(jì)算邊界的選取

為消除邊界效應(yīng)，選取如圖2所示的3個(gè)桶式結(jié)構(gòu)段作為研究對(duì)象。x軸方向?yàn)轳g岸軸線垂直方向，x軸正向指向海側(cè)，從駁岸結(jié)構(gòu)兩側(cè)邊緣向海側(cè)和陸側(cè)各延伸120 m，跨度為240 m；y軸方向?yàn)轳g岸結(jié)構(gòu)軸線方向，y方向?qū)挾?3 m；z方向?yàn)樨Q直方向，模型底部高程取-60 m。模型底部施加3個(gè)方向約束，周圍施加法向約束，水位取平均水位+2.94 m，吹填淤泥至+8.5 m。

圖2　桶式基礎(chǔ)三維有限元分析模型Fig.2　3D finite element model of the bucket-based foundation

2.4計(jì)算工況

本文重點(diǎn)研究桶式基礎(chǔ)結(jié)構(gòu)在施工期和使用期的變形與整體穩(wěn)定性，具體計(jì)算工況及荷載為：施工期荷載為+8.5 m吹填淤泥+施工荷載0.5 kPa，使用期荷載為+7.5 m吹填淤泥（已處理）+堆場(chǎng)使用荷載2～5 t+碼頭線荷載300 kN/m。

3　桶式基礎(chǔ)防波堤三維有限元計(jì)算結(jié)果

3.1位移計(jì)算結(jié)果

前文提到徐圩港區(qū)已進(jìn)行桶式基礎(chǔ)建設(shè)，并且已有完整的位移觀測(cè)資料，所以可根據(jù)現(xiàn)場(chǎng)實(shí)際工況計(jì)算出桶體位移（圖3），再與實(shí)測(cè)位移相比較得出數(shù)值計(jì)算位移修正系數(shù)，進(jìn)而通過修正系數(shù)校對(duì)計(jì)算位移，得出桶體各個(gè)工況下的位移。

現(xiàn)場(chǎng)監(jiān)測(cè)水平位移是1.92 cm，數(shù)值計(jì)算結(jié)果11.06 cm，修正系數(shù)為m=0.181。

根據(jù)以上計(jì)算結(jié)果匯總為表2。

圖3　實(shí)際工況位移（陸域回填4 m厚石頭）Fig.3　Displacement of the bucket-based foundation under actual working condition(backfill 4 m stones)

表2　水平位移匯總表Table2　Summary of horizontal displacement

從水平位移計(jì)算結(jié)果來看，三維有限元彈塑性模型計(jì)算軟土水平位移較實(shí)測(cè)偏大，需通過反演得出的地區(qū)土體水平位移修正系數(shù)進(jìn)行修正。修正后的結(jié)果顯示，桶式基礎(chǔ)結(jié)構(gòu)無論在施工期和使用期結(jié)構(gòu)位移都在安全范圍，結(jié)構(gòu)位移大部分在施工期已發(fā)生，除非使用期結(jié)構(gòu)附近陸域部分存在較大使用荷載，否則不會(huì)發(fā)生較大水平位移。同時(shí)從位移云圖中可知，桶式基礎(chǔ)結(jié)構(gòu)最大水平位移發(fā)生在擋土墻底部和下桶底面，位移呈帶狀分布，這說明桶式基礎(chǔ)結(jié)構(gòu)在底部插入工程性質(zhì)較好的土層前提下，在吹填荷載作用下最可能發(fā)生的是大滑弧破壞，不會(huì)沿著桶底面破壞，間接證明了結(jié)構(gòu)平面尺度的合理性。

3.2承載力系數(shù)計(jì)算結(jié)果

圖4證實(shí)了上文的結(jié)論，即桶式基礎(chǔ)結(jié)構(gòu)在底部插入工程性質(zhì)較好的土層前提下，在吹填荷載作用下最可能發(fā)生的是大滑弧破壞，不會(huì)沿著桶底面破壞。

各階段承載力系數(shù)曲線如圖5所示。

由施工期承載力系數(shù)曲線可知，桶式基礎(chǔ)結(jié)構(gòu)施工期整體穩(wěn)定最小承載力系數(shù)為1.438，故桶式基礎(chǔ)結(jié)構(gòu)安裝完成后，可以一次性吹填成陸，不需要分級(jí)吹填。因此，采用該結(jié)構(gòu)可以大大縮短整個(gè)工程工期。

圖4　位移增量云圖（滑裂面）Fig.4　Incremental displacement cloud(sliding surface)

圖5　承載力系數(shù)曲線Fig.5　Bearing capacity factor curve

由使用期5 t荷載的承載力系數(shù)曲線和使用期5 t荷載淤泥層強(qiáng)度提高的承載力系數(shù)曲線對(duì)比可知，在使用期后方陸域5 t使用荷載工況下，桶式基礎(chǔ)結(jié)構(gòu)整體穩(wěn)定承載力系數(shù)為1.16，原因?yàn)橛倌鄬映休d強(qiáng)度不足發(fā)生破壞，導(dǎo)致承載力系數(shù)偏小，提高淤泥層強(qiáng)度后承載力系數(shù)大幅提高，說明桶式駁岸結(jié)構(gòu)使用期后方陸域若有較高使用荷載，必須對(duì)淤泥層進(jìn)行加固，提高其承載能力。

由前50 m 2 t荷載，后50 m 5 t荷載的承載力系數(shù)曲線可知，桶式駁岸結(jié)構(gòu)在使用期前50 m范圍均載2 t，50～100 m范圍均載5 t工況下，整體穩(wěn)定承載力系數(shù)為1.423。對(duì)比使用期5 t荷載的承載力系數(shù)曲線可以分析出，桶后場(chǎng)地堆載在桶體附近減載可以提高承載力系數(shù)。

由使用期2 t荷載的承載力系數(shù)曲線可知，桶式基礎(chǔ)結(jié)構(gòu)在使用期均載2 t工況下，整體穩(wěn)定承載力系數(shù)為1.496。對(duì)比5 t荷載和范圍減載的承載力系數(shù)曲線可以分析出，減載范圍在50 m以內(nèi)效果明顯，超過50 m對(duì)承載力系數(shù)的提高有限。

4　結(jié)語

結(jié)合已建成的桶式基礎(chǔ)防波堤工程，通過軟土彈塑性本構(gòu)模型，使用三維有限元數(shù)值分析方法，計(jì)算分析了桶式基礎(chǔ)結(jié)構(gòu)的位移和穩(wěn)定性。研究結(jié)論如下：

1）桶式基礎(chǔ)結(jié)構(gòu)在底部插入工程性質(zhì)較好的土層前提下，在吹填荷載作用下最可能發(fā)生的是大滑弧破壞，不會(huì)沿著桶底面破壞。

2）桶式基礎(chǔ)結(jié)構(gòu)在施工期有較大安全儲(chǔ)備，基礎(chǔ)結(jié)構(gòu)安裝完成后，可以一次性吹填成陸，不需要分級(jí)吹填。

3）桶式基礎(chǔ)結(jié)構(gòu)，在使用期的整體穩(wěn)定性，主要取決于使用荷載大小和荷載與結(jié)構(gòu)的距離。為保證使用安全，可采取結(jié)構(gòu)一定范圍內(nèi)限載等措施，并且如有較大使用荷載，必須相應(yīng)提高土體淤泥層的強(qiáng)度。

4）三維有限元法結(jié)合土體彈塑性模型，計(jì)算土-結(jié)構(gòu)作用水平位移時(shí)，結(jié)果偏大，需使用反演得出的土體水平位移修正系數(shù)進(jìn)行修正。

[1]中交第三航務(wù)工程勘察設(shè)計(jì)院有限公司.連云港港徐圩港區(qū)直立式結(jié)構(gòu)東防波堤工程初步設(shè)計(jì)[R].2012. CCCC Third Harbor Consultants Co.,Ltd.Preliminary design of the east up-right breakwater project in Xuwei,Lianyungang Port[R]. 2012.

[2]天津大學(xué).連云港港徐圩港區(qū)防波堤工程桶型基礎(chǔ)結(jié)構(gòu)斷面波浪模型試驗(yàn)[R].2011. Tianjin University.Wave model test of bucket-based structure sections in the breakwater project in Xuwei,Lianyungang Port[R]. 2011.

[3]南京水利科學(xué)研究院.連云港港徐圩港區(qū)防波堤工程離心模型試驗(yàn)[R].2012. Nanjing Hydraulic Research Institute.Centrifugal modeling test of breakwater project in Xuwei,Lianyungang Port[R].2012.

[4]南京水利科學(xué)研究院.新型桶式基礎(chǔ)防波堤與地基動(dòng)力相互作用數(shù)值模擬研究[R].2012. Nanjing Hydraulic Research Institute.Numerical simulation of the dynamic interactions between bucket-based breakwater and the ground[R].2012.

3D FEM numerical analysis on stability of bucket-based structure

MA Zhe-chao,SHENG Jia-jun
（CCCC Third Harbor Consultants Co.,Ltd.,Shanghai 200032,China）

Combined with the soil elastic-plastic constitutive model，we analyzed the displacement and stability of the bucket-based foundation under the filling and loading by use the 3D finite element model.By comparing the results of different conditions and different parameters,we obtained the key factors affecting the stability of bucket based foundation.It is showed that during the construction period,the stability of bucket foundation structure with fill load less than 13 m in height is quite well,while during the service period,the load and the distance where loads are exerted to the structure are the key factors affecting the stability of bucket foundation structure.

bucket based foundation;soft soil ground;displacement;stability

U655.2

A

2095-7874（2016）03-0012-04

10.7640/zggwjs201603003

2016-01-13

江蘇省科技支撐計(jì)劃項(xiàng)目（BE2013663）；江蘇省交通運(yùn)輸科技項(xiàng)目（2013Y20）

馬哲超（1986—），男，遼寧營口人，碩士，助理工程師，巖土工程專業(yè)。E-mail：mazc@theidi.com