Z型鋼板樁碼頭轉(zhuǎn)角段三維數(shù)值模擬

郭浩霖，李春陽，王志斌

(中交第四航務(wù)工程勘察設(shè)計院有限公司，廣東 廣州 510290)

板樁碼頭具有結(jié)構(gòu)簡單、對復(fù)雜的地質(zhì)條件適應(yīng)性強、施工便捷、造價低等優(yōu)點[1]，得到廣泛應(yīng)用。對于挖入式港池而言，板樁碼頭存在向內(nèi)側(cè)或向外側(cè)彎折的轉(zhuǎn)角段，見圖1。

圖1 板樁碼頭轉(zhuǎn)角段

轉(zhuǎn)角段具有明顯的空間特征，目前國內(nèi)外針對板樁碼頭轉(zhuǎn)角段受力與變形的相關(guān)研究較少。板樁碼頭轉(zhuǎn)角段設(shè)計一般采用板樁碼頭標準段的計算成果。為了探討板樁碼頭轉(zhuǎn)角段的受力與變形特性，本文以廣州南沙港某Z型鋼板樁碼頭為例，運用國際大型通用巖土有限元軟件建立三維有限元模型，對板樁碼頭標準段、內(nèi)轉(zhuǎn)角段和外轉(zhuǎn)角段的受力與變形特性進行對比分析，為板樁碼頭轉(zhuǎn)角段的精細化設(shè)計提供參考。

1 工程概況

廣州南沙港某板樁碼頭頂面高程為5.6 m，港池前沿設(shè)計底高程為-5.7 m。前墻采用Z型鋼板樁，頂高程1.9 m、底高程-19 m。錨碇結(jié)構(gòu)采用現(xiàn)澆C40鋼筋混凝土錨碇墻，頂高程4.0 m、墻高4 m、墻厚0.6 m。板樁墻后回填中粗砂，錨碇墻前回填10～100 kg塊石，頂寬10 m。前墻與錨碇墻采用鋼拉桿相連，鋼拉桿間距2.1 m、安裝高程為1.5 m、直徑為70 mm。地基表層及下臥的淤泥和淤泥質(zhì)土采用真空聯(lián)合堆載預(yù)壓進行處理，處理范圍為碼頭前沿線向海側(cè)25 m至碼頭前沿線向陸側(cè)45 m。碼頭典型斷面見圖2。

圖2 碼頭典型斷面(高程：m)

2 模型建立

2.1 土體本構(gòu)模型及單元選取

土體本構(gòu)模型采用小應(yīng)變土體硬化模型(簡稱HSS模型)，該模型為彈塑性雙曲線模型，其在土體硬化模型的基礎(chǔ)上引入了小應(yīng)變屬性[2]，考慮了土的受荷歷史和剛度的應(yīng)變相關(guān)性，可模擬從小應(yīng)變到大應(yīng)變范圍內(nèi)土體的不同響應(yīng)[3]。

為真實模擬板樁-土的相互作用，用板單元模擬板樁結(jié)構(gòu)，用實體單元模擬土體、回填材料，在板樁與土體之間加入接觸界面單元，用點對點錨桿單元模擬拉桿，錨碇墻也采用實體單元。板樁碼頭不同結(jié)構(gòu)段的三維有限元模型見圖3。

圖3 板樁碼頭三維有限元模型

2.2 土體參數(shù)取值

反映土體特性的相關(guān)參數(shù)取值是巖土有限元計算的關(guān)鍵，常規(guī)地勘報告往往無法提供所有土體建模參數(shù)，需要參考相關(guān)的理論研究和工程經(jīng)驗綜合考慮確定。文獻[4-6]對HSS模型整套參數(shù)的取值進行了深入研究，為工程應(yīng)用提供了借鑒。本文土體強度根據(jù)地勘報告的三軸試驗取值，變形模量根據(jù)地勘報告的壓縮模量由經(jīng)驗公式計算確定，土體參數(shù)見表1。

表1 土體參數(shù)

2.3 模型邊界

為了減弱邊界效應(yīng)，模型邊界越大越好，但計算量相應(yīng)增加。綜合考慮模型的精度和計算效率，模型寬度方向(即模型X軸邊界)為碼頭前沿線向海側(cè)50 m至碼頭前沿線向陸側(cè)100 m，共150 m；模型長度方向(即模型Y軸邊界)標準段模型選取了10倍拉桿間距，即21 m，轉(zhuǎn)角段模型基于對稱性考慮，長度與寬度取值相同，為150 m；模型高度方向(即模型Z軸邊界)為碼頭頂高程至⑦1黏土-粉質(zhì)黏土底高程。

2.4 拉桿、錨碇墻交叉區(qū)處理

由于板樁碼頭內(nèi)轉(zhuǎn)角段兩側(cè)拉桿和錨碇墻在平面位置上相互交叉，為了保證兩個方向拉桿和錨碇墻的受力與變形互不干擾，在拉桿交叉處將拉桿沿豎向錯開布置，一側(cè)拉桿安裝高程取1.5 m(對應(yīng)錨碇墻的底高程為0 m)、另一側(cè)取0.5 m(對應(yīng)錨碇墻的底高程為-1 m)，在錨碇墻交叉處設(shè)置結(jié)構(gòu)縫處理。

2.5 施工階段

在軟件中定義施工階段如下：1)場地初始應(yīng)力平衡；2)場地開挖至拉桿和錨碇墻底高程；3)打設(shè)前墻鋼板樁，現(xiàn)澆錨碇墻，安裝拉桿；4)墻后回填，上部結(jié)構(gòu)施工；5)墻前港池開挖至設(shè)計底高程；6)施加碼頭前沿堆貨荷載、剩余水壓力、波吸力或系纜力等荷載。

3 計算結(jié)果及對比分析

3.1 前墻

分別提取標準段、內(nèi)轉(zhuǎn)角段和外轉(zhuǎn)角段前墻的彎矩和位移云圖，見圖4、5。

圖4 前墻彎矩云圖

圖5 前墻位移云圖

經(jīng)數(shù)據(jù)提取分析可知，3種結(jié)構(gòu)段前墻的彎矩和位移最大值基本一致，豎向分布基本相同，在-3.4 m處彎矩最大，在-4.8 m處位移最大。為更清晰地表示前墻的受力變形規(guī)律，分別提取3種結(jié)構(gòu)段前墻水平方向的彎矩和位移最大值，取相反值進行數(shù)據(jù)處理，見圖6、7。

圖6 前墻彎矩最大值相反值分布對比

圖7 前墻位移最大值相反值分布對比

從圖6、7可見，標準段前墻的彎矩和位移水平分布都很均勻，內(nèi)轉(zhuǎn)角段和外轉(zhuǎn)角段彎矩和位移在轉(zhuǎn)角處均最小，隨著與轉(zhuǎn)角處距離的增加而增大，彎矩在距轉(zhuǎn)角點約15 m處達到峰值并趨于平穩(wěn)，位移在距轉(zhuǎn)角點約20 m處達到峰值并趨于平穩(wěn)。

3.2 拉桿

分別提取標準段、內(nèi)轉(zhuǎn)角段和外轉(zhuǎn)角段的拉桿拉力進行數(shù)據(jù)處理，見圖8。

圖8 拉桿拉力分布對比

從圖8可見，3種結(jié)構(gòu)段的拉桿拉力最大值基本一致。標準段的拉桿拉力分布很均勻；內(nèi)轉(zhuǎn)角段在轉(zhuǎn)角處拉桿幾乎不受力，拉桿拉力隨著與轉(zhuǎn)角處距離的增加而增大，在距轉(zhuǎn)角點約20 m處達到峰值后又小幅回落，然后在距轉(zhuǎn)角點約40 m處再次升至峰值并趨于平穩(wěn)；外轉(zhuǎn)角段在轉(zhuǎn)角處拉桿拉力較小，約為最大值的1/6，隨著與轉(zhuǎn)角處距離的增加而增大，在距轉(zhuǎn)角點約20 m處達到峰值并趨于平穩(wěn)。

3.3 錨碇墻

提取標準段、內(nèi)轉(zhuǎn)角段和外轉(zhuǎn)角段錨碇墻的位移云圖，見圖9。

圖9 錨碇墻位移云圖

從圖9位移云圖可見，3種結(jié)構(gòu)段的錨碇墻位移最大值相差不大。分別提取3種結(jié)構(gòu)段前墻水平方向的位移最大值相反值進行數(shù)據(jù)處理，見圖10。

圖10 錨碇墻位移絕對值的最大值分布對比

從圖10可知，標準段錨碇墻水平方向的位移分布很均勻；內(nèi)轉(zhuǎn)角段和外轉(zhuǎn)角段在轉(zhuǎn)角處位移均最小；外轉(zhuǎn)角段錨碇墻的位移隨著與轉(zhuǎn)角處距離的增加而增大，在距轉(zhuǎn)角點約40 m處達到峰值并趨于平穩(wěn)；內(nèi)轉(zhuǎn)角段由于兩個方向的拉桿豎向相互錯開，錨碇墻的高度有差異，因此錨碇墻的位移分布較為復(fù)雜，底高程為0 m的錨碇墻距轉(zhuǎn)角點約15 m處位移最大，底高程為-1m的錨碇墻距轉(zhuǎn)角點約30 m處位移最大。

綜上所述，根據(jù)前墻、拉桿、錨碇墻的受力與變形特征，Z型鋼板樁碼頭的轉(zhuǎn)角段范圍為15～40 m，在該范圍內(nèi)前墻的彎矩和位移、拉桿拉力、錨碇墻位移均比標準段小，這與轉(zhuǎn)角處的邊界效應(yīng)以及拉桿、錨碇墻在該處交錯布置的空間效應(yīng)有關(guān)。

4 結(jié)論

1)Z型鋼板樁碼頭的轉(zhuǎn)角段范圍為15～40 m，在該范圍內(nèi)前墻的彎矩和位移、拉桿拉力、錨碇墻位移均比標準段小，按照標準段的計算成果進行轉(zhuǎn)角段設(shè)計偏于保守，可進一步優(yōu)化。

2)轉(zhuǎn)角段在轉(zhuǎn)角處前墻的彎矩和位移、拉桿拉力、錨碇墻位移均最小。

3)內(nèi)轉(zhuǎn)角段前墻的彎矩、位移隨著與轉(zhuǎn)角處距離的增加而增大，分別在距轉(zhuǎn)角點約15、20 m處達到峰值并趨于平穩(wěn)；拉桿拉力隨著與轉(zhuǎn)角處距離的增加而增大，在距轉(zhuǎn)角點約20 m處達到峰值后又小幅回落，然后在距轉(zhuǎn)角點約40 m處再次升至峰值并趨于平穩(wěn)；錨碇墻的位移分布較復(fù)雜，規(guī)律性不強。

4)外轉(zhuǎn)角段前墻的彎矩和位移分布特征同內(nèi)轉(zhuǎn)角段；拉桿拉力、錨碇墻位移隨著與轉(zhuǎn)角處距離的增加而增大，分別在距轉(zhuǎn)角點約20、40 m處達到峰值并趨于平穩(wěn)。