鋼管樁碼頭地震易損性分析的單自由度模型*

高樹飛，馮云芬

(聊城大學(xué)建筑工程學(xué)院，山東 聊城 252000)

鋼管樁由于水平承載力高，而且彈性和韌性較好，在全直樁碼頭中得到廣泛應(yīng)用。雖然鋼管樁碼頭的承載力和延性相比同等尺度的混凝土樁碼頭要高，但是在強(qiáng)震作用下仍會發(fā)生較大破壞。例如，1995年的日本Hyogoken-Nanbu地震造成神戶港Takahama碼頭的鋼管樁發(fā)生較大破壞，鋼管樁在樁帽和入土部分處均出現(xiàn)較大的彎曲變形，致使碼頭運(yùn)行中斷[1]。因此，對鋼管樁碼頭進(jìn)行合理的抗震設(shè)計(jì)和評估，對于港口的防災(zāi)和減災(zāi)具有十分重要的意義。傳統(tǒng)的高樁碼頭的抗震設(shè)計(jì)方法[2]，無論是基于力的方法，還是基于位移的方法，均是考慮確定性評估碼頭的抗震性能，無法考慮不確定性(包括地震動不確定性和結(jié)構(gòu)不確定性)對碼頭抗震性能的影響，因此近年來從概率角度分析評估碼頭抗震性能正逐漸引起研究人員的關(guān)注，而地震易損性分析方法則是最常用的概率評估方法。該方法通過計(jì)算結(jié)構(gòu)在不同地震動強(qiáng)度下發(fā)生不同破壞狀態(tài)的條件概率，從宏觀角度描述地震動強(qiáng)度對結(jié)構(gòu)破壞的影響，其分析結(jié)果可應(yīng)用于地震引起的災(zāi)害損失評估和抗震加固改造決策，而這一條件概率即為結(jié)構(gòu)易損性。

對于高樁碼頭的地震易損性研究，已有很多學(xué)者開展了研究，但基本上都是國外學(xué)者[3-7]開展的研究，國內(nèi)目前只有馮云芬等[8]開展過相關(guān)研究，而且上述研究均針對混凝土樁碼頭。另外，在碼頭地震易損性分析中，通常須采用非線性動力時(shí)程分析方法對碼頭有限元模型進(jìn)行分析，以確定碼頭在不同地震動強(qiáng)度下的地震反應(yīng)值，再將反應(yīng)值與碼頭結(jié)構(gòu)的不同破壞狀態(tài)的定義值進(jìn)行對比，則可得到結(jié)構(gòu)在不同地震動強(qiáng)度水平下發(fā)生不同破壞狀態(tài)的概率，然后將破壞的概率擬合成易損性曲線[8]。由于碼頭結(jié)構(gòu)的破壞狀態(tài)一般根據(jù)規(guī)范確定，因此合理準(zhǔn)確地確定碼頭的地震反應(yīng)值是易損性分析的關(guān)鍵。在上述研究中，對于碼頭的時(shí)程分析均是直接針對實(shí)際結(jié)構(gòu)開展的，由于復(fù)雜的樁-土相互作用，且易損性分析須采用大量的地震波以考慮地震動的不確定性，故使得易損性分析的計(jì)算量非常大，不便于工程應(yīng)用。目前，在工程結(jié)構(gòu)的非線性時(shí)程分析中，為降低計(jì)算的復(fù)雜性，將原結(jié)構(gòu)簡化為單自由度體系已是一種通行的做法，針對不同結(jié)構(gòu)形式，已有學(xué)者提出了多種單自由度模型[9-11]，但尚未對鋼管樁碼頭易損性分析開展相關(guān)研究。

建立等效單自由度模型的關(guān)鍵是確立合理的恢復(fù)力模型，本文基于光滑型骨架曲線和Masing準(zhǔn)則建立鋼管樁碼頭的恢復(fù)力模型，其中骨架曲線可由碼頭的推覆分析確定；為驗(yàn)證本文提出的單自由度模型在碼頭地震易損性分析中的合理性和有效性，對一個(gè)鋼管樁碼頭案例進(jìn)行易損性分析，比較原型結(jié)構(gòu)和單自由度模型的易損性曲線的差異。

1 單自由度模型的構(gòu)建

考慮到鋼管樁是鋼構(gòu)件，而鋼構(gòu)件本身并不具有負(fù)剛度特性，但由于P-Δ效應(yīng)的影響，其骨架曲線也可能出現(xiàn)負(fù)剛度現(xiàn)象。考慮到碼頭的水平承載力主要由陸側(cè)樁提供(陸側(cè)樁的自由長度小)，通?？珊雎訮-Δ效應(yīng)的影響，故本文不考慮骨架曲線的負(fù)剛度現(xiàn)象。結(jié)合文獻(xiàn)[10]和[12]對鋼管樁碼頭恢復(fù)力特性進(jìn)行研究，本文采用式(1)描述的曲線型骨架線，見圖1。

(1)

式中：F為鋼管樁碼頭水平承載力；Δ為水平位移；Δ1和F1分別為首個(gè)樁塑性鉸形成時(shí)的上部結(jié)構(gòu)水平位移和水平力；α和β為參數(shù)，可按表1確定。文獻(xiàn)[12]在確定α和β時(shí)考慮了樁-土相互作用。對于Δ1和F1，可通過對碼頭進(jìn)行推覆分析予以確定。對于滯回規(guī)則，同樣參考文獻(xiàn)[10]和[12]，采用Masing準(zhǔn)則。

注：Δm為最大水平位移；Fm為最大水平位移對應(yīng)的水平力。

圖1恢復(fù)力模型

表1 參數(shù)α和β的值

注:Su為黏土的不排水抗剪強(qiáng)度，對Su的中間值可通過線性插值予以確定α和β。

在確定好單自由度模型的恢復(fù)力模型后，對于單自由度模型的質(zhì)量m，可以按照式(2)計(jì)算：

m=KT2/(4π2)

(2)

式中：K為碼頭彈性剛度，可取為碼頭Pushover曲線的原點(diǎn)到首個(gè)樁塑性鉸形成點(diǎn)連線的斜率；T為碼頭的自振周期，可對實(shí)際碼頭進(jìn)行振型分析后確定，取為振型參與系數(shù)最大(碼頭橫向反應(yīng))的振型對應(yīng)的周期。對于單自由度模型的阻尼比，一般參考規(guī)范[13]，取為5%。

2 易損性分析

參考文獻(xiàn)[8]，選用位移作為衡量碼頭抗震性能的指標(biāo)，相應(yīng)的易損性函數(shù)可用下式表示：

(3)

式中：C和D分別為鋼管樁碼頭的位移能力和位移需求；IM為地震動的強(qiáng)度參數(shù)；Φ(·)為標(biāo)準(zhǔn)正態(tài)分布函數(shù)；μD為鋼管樁碼頭在不同破壞極限狀態(tài)下的位移能力均值；μD為碼頭在相應(yīng)的IM下的位移需求均值；βC和βD分別為碼頭位移能力和位移需求的對數(shù)標(biāo)準(zhǔn)差。另外，當(dāng)位移能力μC是由Pushover分析確定時(shí)，βC的值可為0.3。

在確定碼頭易損性函數(shù)的方法中，云圖法的計(jì)算量較小，故本文考慮采用該方法確定碼頭易損性函數(shù)。采用云圖法確定碼頭位移需求的概率模型時(shí)，假定位移需求的概率分布為對數(shù)正態(tài)，且D和IM之間滿足以下的指數(shù)關(guān)系：

D(IM)=a(IM)bε

(4)

式中：a和b為參數(shù)，可通過回歸分析計(jì)算；ε為均值是1的對數(shù)正態(tài)隨機(jī)變量，同時(shí)其對數(shù)標(biāo)準(zhǔn)差為βD。對式(4)的等號兩側(cè)分別取對數(shù)可得：

我早已設(shè)計(jì)好，這次一定要雙手捧花，單膝跪地，紳士般求婚，“嫁給我吧?！比缓髶肀В瑹崃业奈?，吻個(gè)死去活來。

ln[D(IM)]=lna+bln(IM)+lnε

(5)

(6)

(7)

式中：Di為第i次非線性時(shí)程分析得到的鋼管樁碼頭位移需求；n為計(jì)算中采用的地震波數(shù)目。由此，利用式(6)和(7)確定μD和βD后，將數(shù)值代入式(3)，則可得到鋼管樁碼頭的易損性函數(shù)。

對于鋼管樁碼頭的破壞極限狀態(tài)，本文考慮文獻(xiàn)[14]定義3個(gè)破壞狀態(tài)，分別為破壞狀態(tài)I-最小破壞、II-可控且可修復(fù)的破壞和III-可保障生命安全的破壞，并用鋼管樁的塑性鉸材料應(yīng)變εs表示，見表2，其中地基深處樁塑性鉸是指塑性鉸所在位置的深度大于10倍樁徑。

表2 鋼管樁碼頭破壞狀態(tài)

3 案例分析

3.1 案例概述

鋼管樁碼頭斷面見圖2，排架間距為6.3 m，碼頭分段長度為47.1 m。鋼管樁直徑為800 mm，壁厚是10 mm，鋼材為Q345。岸坡土體為單一的砂土，其內(nèi)摩擦角為30°。

圖2 鋼管樁碼頭斷面(高程：m；尺寸：mm)

3.2 地震波

地震波從美國太平洋地震工程研究中心的強(qiáng)震數(shù)據(jù)庫中選取，共50條。限于篇幅，地震波的相關(guān)信息從略。

3.3 易損性分析結(jié)果

采用SAP2000建立碼頭有限元數(shù)值分析模型，相關(guān)建模方法可參考文獻(xiàn)[12]。首先進(jìn)行碼頭推覆分析，可得碼頭的荷載-變形曲線(Pushover曲線)，見圖3，圖中給出了首個(gè)鋼管樁塑性鉸出現(xiàn)時(shí)對應(yīng)的位置，并給出了不同破壞狀態(tài)對應(yīng)的位置。不同破壞狀態(tài)I、II和III下的μC分別為8.7、13.2、17.4 cm。

圖3 Pushover曲線

利用推覆分析得到的首個(gè)樁塑性鉸形成時(shí)的碼頭水平位移Δ1和水平力F1，可由式(1)確定等效單自由度模型的骨架曲線，再結(jié)合Masing準(zhǔn)則和式(2)即可建立碼頭的單自由度模型。采用云圖法確定碼頭易損性時(shí)需要計(jì)算每一條地震波下實(shí)際碼頭和單自由度模型的位移需求，但是需要指出的是，本文中的單自由度模型只是用于確定碼頭的橫向位移反應(yīng)，而計(jì)算碼頭總位移需求則需要考慮碼頭扭轉(zhuǎn)和雙向水平地震動的影響，本文采用文獻(xiàn)[14]的動力放大系數(shù)DMF來考慮這一影響，即在計(jì)算確定的碼頭橫向位移需求的基礎(chǔ)上乘以DMF，進(jìn)而得到碼頭的總位移需求D。DMF的計(jì)算公式如下：

(8)

式中：e0為碼頭的初始偏心距；Ll是碼頭分段長度。對于本案例，經(jīng)計(jì)算可得e0=6.40 m、Ll=47.1 m，代入式(8)后可得DMF=1.498。

通過非線性時(shí)程分析確定碼頭的總位移需求后，考慮地震動強(qiáng)度參數(shù)IM選用峰值地面加速度PGA，可得位移需求D和PGA之間的關(guān)系，見圖4，再根據(jù)式(6)和式(7)則可確定μD和βD。對于實(shí)際碼頭和對應(yīng)的單自由度模型，μD可分別由式(9)和(10)計(jì)算，βD分別為0.371 5、0.341 5 cm。

lnμD=0.696 8ln(PGA)+3.081 8

(9)

lnμD=0.613 0ln(PGA)+2.992 8

(10)

圖4 位移需求D的回歸分析

圖5 易損性曲線

將式(9)和(10)以及計(jì)算得到的βD代入式(3)，可得實(shí)際碼頭和單自由度模型的易損性函數(shù)，將函數(shù)用幾何圖形表示出來，即可得易損性曲線，見圖5。由圖5不難看出，隨著PGA的增大，碼頭結(jié)構(gòu)發(fā)生破壞的概率逐漸增大，而且發(fā)生破壞狀態(tài)III的概率低于破壞狀態(tài)I和II，說明碼頭發(fā)生不可修復(fù)的破壞的概率較低；總體上實(shí)際碼頭和單自由度模型的易損性曲線相差不大，特別是在破壞狀態(tài)I的情況下；在峰值地面加速度PGA小于0.5g的情況下，二者的差別很小，但隨后二者的差距變大?？紤]到二者易損性曲線的差別很小，因此在對實(shí)際鋼管樁碼頭進(jìn)行易損性分析時(shí)，可以按本文的方法構(gòu)建單自由度模型以代替原結(jié)構(gòu)進(jìn)行分析，以減小計(jì)算量和計(jì)算難度；另外，考慮到單自由度模型和實(shí)際碼頭易損性曲線在PGA大于0.5g的情況下有一定差別，故在此情況下的單自由度模型計(jì)算結(jié)果應(yīng)慎重使用。

4 結(jié)論

1) 采用曲線型骨架線和Masing準(zhǔn)則構(gòu)建的單自由度模型可代替實(shí)際鋼管樁碼頭進(jìn)行易損性分析，其易損性曲線與實(shí)際碼頭的曲線存在一定差別，但差別不大。

2) 采用單自由度模型進(jìn)行鋼管樁碼頭易損性分析可以極大地降低計(jì)算量，但仍然無法避免進(jìn)行大量的非線性時(shí)程分析，在下一步研究中可考慮對單自由度模型進(jìn)行參數(shù)分析，研究易損性曲線和相關(guān)參數(shù)的關(guān)系，進(jìn)而建立標(biāo)準(zhǔn)化的碼頭易損性曲線供碼頭抗震評估使用。