不同連接屬性下低樁承臺重力式碼頭動力特性

周 丹，王多銀，趙大淮，段倫良

(重慶交通大學 河海學院，重慶400074)

隨著庫區水位的日漸升高，地基條件良好的碼頭岸線越來越少，部分碼頭不得不修建在極端地基上[1]。重力式碼頭因施工相對簡單、成本相對低廉、適應性強等優點而得到廣泛應用[2]。樁-土-承臺協同發揮作用過程中，樁側土體出現的成拱效應有效改善地基承載力[3]，因此在深水區橋梁碼頭等工程中應用十分廣泛[4-5]。低樁承臺重力式碼頭兼有樁基和重力式碼頭的優點，能很好地適用于內河軟土地區碼頭的建設。近年來國內外歷次地震災害調查結果[6-7]表明：碼頭在地震作用下極易遭到破壞，從而造成巨大經濟損失。因此研究低樁承臺重力式碼頭在地震作用下的動力響應的重要性不言而喻。

重力式碼頭與基礎之間以摩擦接觸的方式進行連接，建筑和橋梁中樁基-承臺與上部結構則采用剛性連接，即錨固連接。為便于描述，將重力式碼頭稱為上部結構，將樁-承臺結構稱為下部結構。地震作用下，低樁承臺重力式碼頭上下部結構連接節點的工作形態直接影響碼頭整體結構安全性能。目前還沒有低樁承臺重力式碼頭上下部結構中采用錨固連接的相關研究，本文在低樁承臺重力式碼頭上下部結構之間分別采用錨固連接和摩擦接觸連接(圖1)，利用ABAQUS有限元軟件對兩種連接方式形成的低樁承臺重力式碼頭整體結構建立有限元模型，在地震荷載作用下進行動力時程分析，來探討不同連接屬性對結構動力響應的影響，為今后該種結構的抗震設計提供參考。

圖1 錨固連接與摩擦接觸連接

1 分析方法與求解

1.1 動力平衡方程的建立

根據平衡條件，建立低樁承臺重力式碼頭結構在地震作用下的平衡方程[8]：

(1)

1.2 動力平衡方程的求解

求解動力平衡方程的基本方法有：時域分析法、頻域分析法和振型疊加法，在數值分析中常采用振型疊加法和直接積分法，本文采用基于New-mark法的隱式積分法[9]：

(2)

(3)

(4)

(5)

(6)

(7)

(8)

2 有限元模型的建立

2.1 工程概況

某500噸級的低樁-承臺塊體階梯型結構重力式碼頭，所處場地為8度地震烈度帶，主要建筑物等級為一級，該碼頭設計高水位為6.245 m，設計低水位為4.300 m；碼頭所處位置地勢陡峭，其地基土土層由下至上分別為頁巖、黏土和回填土。取沿河向長度3.5 m進行研究，承臺底寬3.5 m、高1.8 m，胸墻高2 m，采用兩樁承臺重力式碼頭結構，樁徑D為500 mm。結構斷面見圖2。

圖2 結構斷面(單位：mm)

2.2 有限元模型及相關參數

根據實際工程建立同尺寸的三維低樁承臺重力式碼頭結構-土有限元模型，下部為地基土，樁底嵌入頁巖，塊體碼頭后為回填土。本文取樁軸外半徑12D以內的土體、取樁端以下厚15D的地基土進行研究[10]。模型設計的材料相關參數見表1。

表1 模型參數

模型中重力式碼頭采用鋼筋混凝土結構，樁基采用鋼筋混凝土結構，按構造配筋。計算模型中混凝土采用塑性損傷模型[11]，地基土與回填土選用Mohr-Coulom塑性模型[12-13]，鋼筋采用線彈性模型[14]。結構與土之間的接觸部位存在滑移或脫開等現象，屬于非線性問題，在ABAQUS中通過“接觸對” 實現；鋼筋與混凝土分別采用不同的單元分離建模，之后進行自由度耦合。

模型荷載中水上土體按照天然密度考慮，水下土體按照飽和密度考慮，水壓力按靜水壓力考慮，地震荷載加速度時程曲線見圖3。地基在未受到擾動時就存在天然應力，對分析過程產生了影響，因此本文采用自動平衡，即在step模塊中創建geostatic分析步進行初始應力平衡。考慮樁-土-承臺共同受力[15]，建立碼頭結構與土體有限元模型，見圖4。

圖3 0.2g EL-Centro波加速度時程曲線

圖4 碼頭-土體有限元模型

3 不同上下部結構連接屬性對結構動力響應的影響

3.1 結構水平動位移響應

通過ABAQUS精細化建模，得到兩種連接屬性對結構水平動位移影響規律，見圖5及表2。由圖5可知，在地震作用下，兩種連接屬性所對應的上下部結構水平動位移的變化趨勢基本一致，且上部結構的水平動位移始終大于下部結構；在地震過程中，采用摩擦接觸和錨固連接所對應結構水平動位移的最大值分別為-5.1、-3.6 mm，地震結束時，兩種連接屬性所對應結構的殘余位移分別為-1.7、-1.2 mm。

由表2可知，在地震荷載作用下，采用摩擦接觸連接的低樁承臺重力式碼頭上下部結構水平動位移之差的最大值、重力式碼頭頂部與底部水平動位移之差的最大值及在地震結束時的殘余水平動位移之差均比采用錨固連接時大。

圖5 上下部結構水平動位移最大值時程曲線

表2 上下部結構水平位移差值

3.2 結構水平加速度響應

兩種連接屬性對結構水平加速度的影響規律見圖6及表3。

圖6 上下部結構水平峰值加速度時程曲線

表3 上下部結構峰值加速度絕對值之差

由圖6可知，在地震作用下，兩種連接屬性所對應的上下部結構的水平加速度隨時間的變化趨勢基本一致，且上部結構水平加速度絕對值均大于下部結構；由于地震響應的遲滯性，上下部結構達到峰值水平加速度的時間點均稍晚于輸入地震荷載峰值加速度的時間點；在地震結束時，上下部結構的水平加速度都趨于0。

由表3可知，摩擦接觸連接的結構達到峰值加速度的時間點晚于采用錨固連接結構；摩擦接觸連接所對應的上下部結構水平加速度之差的最大值比采用錨固連接結構大；錨固連接所對應的上部結構峰值加速度占輸入地震荷載加速度峰值的比率比摩擦接觸連接結構大。這都是由于采用錨固連接結構的整體性要好于摩擦接觸連接的結構。

3.3 結構應力響應

兩種連接方式對低樁承臺重力式碼頭結構應力響應的影響規律見圖7～11及表4。

圖7 結構應力云圖

圖8 碼頭底面拉應力云圖

圖9 碼頭底面壓應力云圖

圖10 承臺頂面拉應力云圖

圖11 承臺頂面壓應力云圖

表4 上下部結構峰值應力及其比值

由圖7可知，兩種連接方式所對應結構整體應力都隨高度的增加而減小；下部結構的應力都大于上部結構，樁身應力分布不均勻且樁底部應力最大。由圖8、9可知，在地震荷載作用下，兩種連接屬性所對應的重力式碼頭底面的拉應力均沿河流方向呈帶狀分布，且錨固連接重力式碼頭底面的應力最大值較摩擦接觸結構大；錨固連接所對應的重力式碼頭底面的壓應力分布出現一定的集中現象，而摩擦接觸所對應的重力式碼頭底面的壓應力分布較分散。由圖10、11可知，摩擦接觸所對應的結構承臺頂面的峰值拉壓應力均小于錨固連接所對應的結構，這也是由于摩擦接觸能在一定程度上緩解應力集中所致。

由表4可知，由于錨固連接應力傳遞效率明顯好于摩擦接觸連接屬性，使得錨固連接所對應結構的最大、最小拉(壓)應力均大于摩擦接觸連接屬性所對應的結構，且錨固連接所對應結構的最大與最小拉應力之比及最大與最小壓應力之比均大于摩擦接觸連接屬性所對應的結構。

3.4 結構動土壓力響應

兩種連接屬性對結構動土壓力響應影響見圖12、13及表5。

圖12 結構動土壓力沿高度變化曲線

圖13 上下部結構動土壓力最大值時程曲線

表5 上下部結構動土壓力及其比值

由圖12可知，兩種連接屬性所對應的結構沿高度方向的動土壓力均呈先振蕩減小再線性減小趨勢；錨固連接所對應的結構沿結構高度方向的動土壓力始終小于摩擦接觸所對應的結構。由圖13可知，兩種連接屬性所對應上下部結構動土壓力隨時間的變化趨勢基本一致，下部結構的動土壓力均大于上部結構。

由表5可知，兩種連接屬性所對應的下部結構的峰值動土壓力均大于上部結構，摩擦接觸所對應的上下部結構最大動土壓力值、最大與最小動土壓力之比均大于錨固連接結構，這說明，摩擦接觸連接所對應的結構在地震荷載作用下更容易發生應力集中現象。

4 結論

1)兩種連接屬性所對應結構的水平動位移，其上下部結構變化趨勢基本一致，都隨結構高度增加逐漸增大，且上部結構始終大于下部結構；在地震結束時，錨固連接所對應結構的殘余位移比摩擦接觸結構小。

2)兩種連接屬性所對應結構的水平加速度絕對值，都隨高度的增加而增大，上部結構均大于下部結構；摩擦接觸所對應的結構達到峰值加速度的時間點明顯滯后于錨固連接結構，摩擦接觸所對應結構的上下部結構水平峰值加速度之差大于采用錨固連接屬性的結構；錨固連接所對應的上部結構峰值加速度占輸入地震荷載加速度峰值的比率大于采用摩擦接觸連接的結構。

3)錨固連接所對應的重力式碼頭底面的峰值拉、壓應力及承臺頂面的峰值拉、壓應力均大于摩擦接觸連接所對應的結構；錨固連接所對應結構的最大與最小拉(壓)應力均小于摩擦接觸連接屬性所對應的結構。

4)兩種連接屬性所對應上下部結構的動土壓力沿結構高度均呈先振蕩減小再線性減小趨勢，錨固連接所對應的結構沿結構高度方向的動土壓力始終小于摩擦接觸所對應的結構；摩擦接觸所對應的上下部結構最大動土壓力值、最大與最小動土壓力之比均大于錨固連接結構。