陳光明， 張淑華 ， 袁 宇， 劉春陽， 孫興毅

(河海大學(xué) 港口海岸與近海工程學(xué)院，江蘇 南京210098)

樁基平臺(tái)具有抵御波浪和地震的能力，廣泛應(yīng)用于海洋石油開采。深水樁基平臺(tái)在地震作用下會(huì)引起結(jié)構(gòu)周圍水體的輻射波浪運(yùn)動(dòng)，水體波動(dòng)將對(duì)樁基平臺(tái)水下部分產(chǎn)生動(dòng)水壓力。該動(dòng)水壓力不僅會(huì)改變結(jié)構(gòu)的動(dòng)力特性，還會(huì)影響結(jié)構(gòu)的動(dòng)力響應(yīng)［1-2］。因此，動(dòng)水壓力問題是一個(gè)十分復(fù)雜的問題。

目前國內(nèi)外學(xué)者對(duì)深水中樁柱結(jié)構(gòu)動(dòng)水壓力方面做了不少研究。羅傳信等［3］進(jìn)行波浪與地震響應(yīng)研究，分析得出動(dòng)水阻尼對(duì)大波高、長周期波浪荷載下平臺(tái)響應(yīng)影響較大，在進(jìn)行動(dòng)力分析時(shí)應(yīng)加以考慮。高學(xué)奎［4］采用附加質(zhì)量的方法，結(jié)合有限元方法來分析深水橋墩的地震響應(yīng)，認(rèn)為地震動(dòng)水壓力對(duì)深水橋墩的影響是顯著的。劉洋［5］利用ADINA 軟件建立導(dǎo)管架海洋平臺(tái)有限元模型，分析表明考慮附加水質(zhì)量和動(dòng)水阻尼能夠減小平臺(tái)的頻率，并且驗(yàn)證得到阻尼對(duì)平臺(tái)的頻率影響很小。李富榮等［6］在Morison 方程的基礎(chǔ)上，以單柱式橋墩為研究對(duì)象，用附加水質(zhì)量法考慮動(dòng)水壓力對(duì)橋墩的影響，結(jié)果表明考慮動(dòng)水壓力效應(yīng)是有必要的。前人研究主要集中水體阻尼及水體附加質(zhì)量引起的水體慣性力對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的影響，很少考慮樁土非線性作用。文中根據(jù)動(dòng)力時(shí)程分析法的基本理論，基于Ansys 有限元軟件，建立樁基平臺(tái)-土三維模型，分析附加水質(zhì)量和動(dòng)水阻尼對(duì)樁基響應(yīng)的影響規(guī)律，并比較樁基空間布置及傾斜方位對(duì)其內(nèi)力的影響規(guī)律。

1 高樁承臺(tái)有限元模型和荷載環(huán)境

1.1 高樁承臺(tái)工程概況

針對(duì)某一工程為研究對(duì)象，承臺(tái)厚度2 m，采用9 根3 ∶1 斜樁，泥面處高程為0 m，樁底高程－42 m。高樁承臺(tái)布置見圖1，樁基布置見圖2 和表1，有限元模型見圖3，土層物理性質(zhì)參數(shù)如表2 所示。

圖1 高樁承臺(tái)布置Fig.1 Location of the elevated pile cap

圖2 樁位置布置Fig.2 Location of the pile

表1 樁的空間位置布置Tab.1 Spatial position of the pile arrangement

圖3 有限元模型Fig.3 Finite element model

表2 土層物理性質(zhì)參數(shù)Tab.2 Parameters of the soil physical property

1.2 高樁承臺(tái)有限元模型

承臺(tái)結(jié)構(gòu)材料為C40 混凝土，采用solid95 單元，單元尺寸為1 m。樁基材料為Q235B 鋼管樁，直徑為0.5 m，壁厚18 mm，采用beam188 單元，單元尺寸為1 m。

1.3 荷載環(huán)境

1.3.1 動(dòng)水壓力荷載 以Morison 方程為基礎(chǔ)，把水等效為附加質(zhì)量來考慮水對(duì)結(jié)構(gòu)的影響，忽略樁基對(duì)海水運(yùn)動(dòng)的影響，認(rèn)為海水對(duì)樁基的作用由未受擾動(dòng)的加速度場(chǎng)和速度場(chǎng)引起的沿水運(yùn)動(dòng)方向作用于樁基上的慣性力和阻尼力所引起，采用簡化的Morison 方程表示地震動(dòng)引起的樁基動(dòng)水壓力:

其中，MW= (CM－1)ρV 為動(dòng)水附加質(zhì)量系數(shù)，CW=0.5CDρAPσ為動(dòng)水附加阻尼系數(shù)。則節(jié)點(diǎn)i 處等效附加水質(zhì)量為

式中，海水密度ρ 取1.2 ×1.03kg/m3，R1為鋼管樁的外徑0.5 m，li為第i 單元有效長度為1 m，《海港水文規(guī)范》［7］對(duì)于圓柱體結(jié)構(gòu)取CM= 2.0，CD=1.2，σ = 0.5。

文中附加水質(zhì)量采用MASS21 單元模擬，動(dòng)水阻尼采用COMBIN14 模擬，樁土非線性作用采用COMBIN39 模擬，具體如圖4，5 所示。

1.3.2 地震荷載 文中選取前10 s 時(shí)間內(nèi)EL-Centro 波，且只考慮水平地震荷載作用。該區(qū)域地震設(shè)防烈度為8 度，地震波的調(diào)幅系數(shù)為0.2 g，需對(duì)加速度曲線上每個(gè)點(diǎn)乘以0.2 g 后施加于結(jié)構(gòu)上(見圖6)。

圖4 動(dòng)水壓力分析模型Fig.4 Hydrodynamic pressure analysis model

圖5 樁土非線性模型Fig.5 Nonlinear model of the pile-soil

圖6 EI-Centrol 地震波輸入Fig.6 Time-history curves of the EI-Centrol motion

2 樁土非線性作用和模態(tài)分析

埋入土壤的樁基部分所受土壤非線性作用力通過非線性彈簧單元COMBIN39 模擬。樁基和土體之間的阻尼系數(shù)為0.02。根據(jù)表2 中的內(nèi)摩擦角、不排水抗剪強(qiáng)度等參數(shù)，結(jié)合《海上固定平臺(tái)規(guī)劃、設(shè)計(jì)和建造的推薦作法工作應(yīng)力設(shè)計(jì)法》［8］，計(jì)算出樁土的橫向荷載-位移傳遞曲線(p－y 曲線)、軸向荷載－位移傳遞曲線(t－z 曲線)以及樁端荷載－位移傳遞曲線(q－z 曲線)，然后將荷載－位移傳遞曲線離散建立非線性彈簧單元實(shí)常數(shù)。沿樁基橫截面x，y 方向施加單元長度為1 的非線性彈簧，沿樁的軸線方向施加單元長度為1 的非線性彈簧，具體見圖4，5。

2.1 樁土側(cè)向荷載-位移傳遞曲線

根據(jù)室內(nèi)三軸試驗(yàn)和現(xiàn)場(chǎng)試樁樁周土極限土抗力一半時(shí)的樁撓曲變形y50= Aε50B，式中:A =2.5，B = 0.5 m，ε50= 0.02。

極限水平土抗力沿樁身分布，當(dāng)x ＜xr，pu=3Cu+ γx + JCux/B，x ＞ xr，pu= 9Cu，其中xr=6B/(γB/Cu+ J)，x 為從泥面算起的任意深度，J =0.5。

水下軟黏土在周期荷載作用下，當(dāng)y ＜3y50時(shí)，pu= 0.5(y/y50)1/3;當(dāng)y ≥3y50時(shí)，分為兩種情況:(1)x ≥xr時(shí)，p = 0.72pu;(2)x ＜xr時(shí)，當(dāng)y ≥15y50，p = 0.72pu(x/xr)，3y50＜ y ＜ 15y50，為 一直線。

2.2 樁土軸向荷載－位移曲線

對(duì)于打入黏性地基的樁，其樁身側(cè)摩阻力fs和樁端極限阻力qp的一般表達(dá)式:fs= αCu，qp= 9Cu，當(dāng)ψ ≤1.0，α = 0.5ψ－50，當(dāng)ψ ＞1.0，α = 0.5ψ－0.25，式中，ψ = Cu/P0，P0為相應(yīng)點(diǎn)的有限覆蓋土壓力。

2.3 模態(tài)分析

利用Ansys 做動(dòng)力時(shí)程分析時(shí)需要采用Alpha阻尼和Beta 阻尼來定義瑞利(Rayleigh)阻尼常數(shù)α和β。阻尼矩陣是在用這些常數(shù)乘以質(zhì)量矩陣［M］和剛度矩陣［K］后計(jì)算出來的。通常α 和β 值不是直接得到的，而是用振型阻尼比ξi計(jì)算出來的。ξi是某個(gè)振型i 的實(shí)際阻尼和臨界阻尼之比。根據(jù)《建筑抗震設(shè)計(jì)規(guī)范》［9］的規(guī)定:文中取ξ1= ξ2= 0.05。

其中，ωi為結(jié)構(gòu)的第1 階和第2 階自振頻率，如表3所示。

表3 結(jié)構(gòu)自振頻率Tab.3 Natural frequencies of the structures

3 有限元分析結(jié)構(gòu)與分析

3.1 動(dòng)水壓力對(duì)高樁承臺(tái)結(jié)果分析

為了分析高樁承臺(tái)動(dòng)水壓力的影響，分別計(jì)算考慮附加水質(zhì)量和動(dòng)水阻尼、考慮附加水質(zhì)量但不考慮動(dòng)水阻尼、不考慮附加水質(zhì)量也不考慮動(dòng)水阻尼時(shí)各樁基的位移峰值和樁頂彎矩峰值，計(jì)算結(jié)果如表4，5 所示。

表4 位移峰值Tab.4 Peak of the displacement

從表4 可知，沿地震波傳播方向的同一列樁基同時(shí)達(dá)到最大水平位移，且大小相等。附加水質(zhì)量使平臺(tái)位移峰值增大20% 左右，動(dòng)水阻尼使平臺(tái)位移峰值減?。?4% 左右，同時(shí)考慮附加水質(zhì)量和動(dòng)水阻尼時(shí)，樁1、樁4、樁7 號(hào)樁平臺(tái)位移峰值增大4.10%，樁2、樁5、樁8 號(hào)樁平臺(tái)位移峰值增大2.78%，樁3、樁6、樁9 號(hào)樁平臺(tái)位移峰值增大1.12%，且共同作用時(shí)不具有線性疊加效應(yīng)。

根據(jù)表5 中樁1 ～樁9y 方向和z 方向的峰值彎矩可知，附加水質(zhì)量和動(dòng)水阻尼對(duì)彎矩峰值的減小趨勢(shì)與樁空間布置存在一定的相關(guān)性。樁1 和樁2平行，樁1 和樁2 的增大或減小效應(yīng)相似;y 方向的峰值彎矩減小8.8% 左右;z 方向的峰值彎矩減小3.5% 左右;樁3 和樁1 沿地震加速度傳播方向上空間對(duì)稱，y，z 方向的峰值彎矩減小幅度呈現(xiàn)反對(duì)稱趨勢(shì)。樁4 和樁5 近似平行，y 方向的峰值彎矩減小幅度具有相似性。樁6、樁8、樁9 傾斜方向近似相同，y方向和z 方向的峰值彎矩減小11% 左右。樁7 與樁8沿地震加速度傳播方向上空間對(duì)稱，y 方向的峰值彎矩減小6% 左右，z 方向的峰值彎矩減小15%。綜上可得，當(dāng)樁基空間位置平行時(shí)，樁頂彎矩峰值存在相似的增大或減小趨勢(shì)，當(dāng)樁基之間近似空間對(duì)稱時(shí)，則出現(xiàn)相反的規(guī)律?？傊?，附加水質(zhì)量使彎矩峰值增大，且增大效應(yīng)不明顯，動(dòng)水阻尼使彎矩峰值減小，減小效應(yīng)明顯，兩者共同作用時(shí)具有線性疊加效應(yīng)。

表5 彎矩峰值Tab.5 Peak of the bending moment

3.2 考慮附加水質(zhì)量和動(dòng)水阻尼作用下各方位樁之間內(nèi)力分析

根據(jù)樁身彎矩峰值表，彎矩最大值發(fā)生在t =2.28 s 附近，此時(shí)加速度指向X 負(fù)方向，根據(jù)水平荷載作用方向和樁傾斜方向，此時(shí)樁1 ～樁5 為正斜樁，樁6 ～樁9 為負(fù)斜樁。

樁1 和樁2 空間平行，樁1 和樁3 關(guān)于樁2 所在的鉛直面空間對(duì)稱。圖7，8 表明，樁3 的彎矩峰值最小，樁1 的彎矩峰值最大，樁2 總體與樁1 的彎矩和位移相似，在樁頂附近有部分差異。樁7 和樁8 關(guān)于它們之間的鉛直面空間對(duì)稱，樁8 和樁9 空間近似平行。圖9，10 表明，樁7、樁8、樁9 的彎矩峰值圖基本重合，泥面以上樁9 位移峰最大，樁7 的位移峰值最小，泥面以下3 根樁位移峰值圖基本重合。圖7，9表明負(fù)斜樁和正斜樁承受方向相反，且負(fù)斜樁比正斜樁能承受更大的彎矩。圖11 表明沿地震加速度傳播方向上，樁頂位移相同。

圖7 1 ～3 號(hào)樁彎矩峰值Fig.7 Peak of the bending moment on 1 ～3 piles

圖8 1 ～3 號(hào)樁位移峰值Fig.8 Peak of the bending moment on 1 ～3 pile

圖9 7 ～9 號(hào)樁彎矩峰值Fig.9 Peak of moment on 7 ～9 piles

圖10 7 ～9 號(hào)位移峰值Fig.10 Peak of displacement on 7 ～9 piles

圖11 1，4，7 號(hào)樁位移峰值Fig.11 Peak of displacement on 1，4，7 piles

4 結(jié) 語

基于Ansys 有限元軟件，建立了樁基平臺(tái)-土系統(tǒng)的地震反應(yīng)分析三維有限元模型。通過p－y，t－z，q－z 曲線考慮樁土之間的非線性作用，重點(diǎn)比較了附加水質(zhì)量和動(dòng)水阻尼高樁承臺(tái)的位移和樁基內(nèi)力影響，并討論了樁基空間布置對(duì)結(jié)構(gòu)內(nèi)力影響，結(jié)果如下:

1)附加水質(zhì)量使彎矩產(chǎn)生小幅度的增大效應(yīng)，動(dòng)水阻尼使彎矩產(chǎn)生較大的減小效應(yīng)，且占主導(dǎo)地位，兩者共同作用時(shí)具有線性疊加效應(yīng)。

2)附加水質(zhì)量使位移產(chǎn)生增大效應(yīng)，動(dòng)水阻尼使彎矩產(chǎn)生較大的減小效應(yīng)，且附加水質(zhì)量產(chǎn)生的增大效應(yīng)大于動(dòng)水阻尼的減小效應(yīng)，兩者共同作用時(shí)不具有線性疊加效應(yīng)。

3)沿地震傳播方向上，同一列樁的樁頂位移基本相同。

4)附加水質(zhì)量和動(dòng)水阻尼對(duì)樁基彎矩峰值的增大或減小效應(yīng)存在一定的相關(guān)性，當(dāng)樁基空間位置平行時(shí)，樁頂彎矩峰值存在相似的增大或減小趨勢(shì)，當(dāng)樁基之間近似空間對(duì)稱時(shí)，則出現(xiàn)相反的規(guī)律。

5)在考慮樁土非線性和動(dòng)水壓力作用時(shí)，高樁承臺(tái)樁基正斜樁的承載力大于負(fù)斜樁的承載力。

［1］唐寰澄.世界跨海工程概況［C］//中國土木工程學(xué)會(huì)橋梁及結(jié)構(gòu)工程學(xué)會(huì)第十二屆年會(huì)論文集(上冊(cè)).廣州:中國土木工程學(xué)會(huì)，1996:93-103.

［2］J.F.威爾遜.海洋結(jié)構(gòu)動(dòng)力學(xué)［M］.楊國金，譯.北京:石油工業(yè)出版社，1991.

［3］羅傳信，劉春山，唐紅英.動(dòng)水阻尼對(duì)樁基平臺(tái)動(dòng)力響應(yīng)的影響問題［J］.海報(bào)，1989(11):794-800.

LUO Chuanxin，LIU Chunshan，TANG Hongying.Dynamic water damping effect on dynamic response of pile foundation platform problem［J］.Journal of Posters，1989(11):794-800.(in Chinese)

［4］高學(xué)奎，朱晞.地震動(dòng)水壓力對(duì)深水橋墩的影響［J］.北京交通大學(xué)學(xué)報(bào)，2006，30(1):55-58.

GAO Xuekui，ZHU Xi. The influence of the vibration water pressure of deep water piers［J］. Journal of Beijing Jiaotong University，2006，30(1):55-58.(in Chinese)

［5］劉洋.導(dǎo)管架海洋平臺(tái)海嘯荷載作用下的受力分析［D］.大慶:東北石油大學(xué)，2013.

［6］李富榮，陳國興，王志華.考慮動(dòng)水壓力影響的單柱式橋墩地震反應(yīng)分析［J］. 地震工程與地震工程震動(dòng)，2008，28(2):114-121.

LI Furong，CHEN Guoxing，WANG Zhihua.Considering the influences of dynamic water pressure single column pier earthquake response analysis［J］.Earthquake Engineering and Engineering Vibration，2008，28(2):114-121.(in Chinese)

［7］交通部第一航務(wù)工程勘察設(shè)計(jì)院.JTS145-2-2013 海港水文規(guī)范［S］.北京:人民交通出版社，2013.

［8］海洋石油工程專業(yè)標(biāo)準(zhǔn)化委員會(huì).SY T10030-2004 海上固定平臺(tái)規(guī)劃、設(shè)計(jì)和建造的推薦作法工作應(yīng)力設(shè)計(jì)法［S］.北京:中國石油出版社，2004.

［9］中華人民共和國住房和城鄉(xiāng)建設(shè)部.GB50011-2010 建筑抗震設(shè)計(jì)規(guī)范［S］.北京:中國建筑工業(yè)出版社，2010.