梁 巖, 張卓航， 班亞云， 王 艷

(鄭州大學 土木工程學院，河南 鄭州 450001)

公路橋梁是交通工程的重要組成部分，近年來國內外橋梁抗震性能評價、加固改造實踐和震害經驗表明，對公路橋梁進行抗震性能評價，并對不滿足評價要求的橋梁采取適當的抗震加固改造對策，是減輕地震災害的重要途徑[1]。橋梁地震易損性分析是橋梁抗震性能評價的重要依據[2]，可預測在地震作用下，結構達到或超越某一損傷指標的概率。易損性分析方法在考慮地震動隨機性及結構損傷不確定性的同時，能夠有效地將橋梁在地震中的失效概率和橋梁的承載能力建立起某種概率上的聯系。戶東陽[3]基于往年收集的渝昆高鐵沿線實測地震動記錄資料，研究了不同形式支座的抗震性能，以及橋墩的地震易損性。王君[4]、陳偉等[5]均以國內典型的高速鐵路橋梁為工程背景，采用地震時程分析方法(IDA)對橋梁的地震易損性進行分析。Basoez等[6]同樣根據收集的震害資料，依據邏輯回歸法，分析地震作用下橋梁構件的易損性。

隨著國內高烈度區(qū)公路橋梁的建造，多跨連續(xù)梁-剛構橋逐漸增多，研究該橋型地震易損性對其抗震性能評估有重要意義。而國內外對橋梁地震易損性的研究多集中于橋墩，對于耐久性損傷下橋墩抗震性能的研究僅針對帶支座的橋墩，未考慮墩頂約束作用，未重視服役期內材料性能的退化，且忽略了橋梁系統(tǒng)的損傷概率與單一構件的損傷概率之間的差異性，僅以構件的損傷概率分析地震作用下橋梁整體的損傷，導致結構存在安全隱患。

本文基于地震易損性分析理論對某多跨連續(xù)梁-剛構橋梁構件(橋墩、支座、橋臺)進行易損性研究，分析橋梁構件在地震作用下的失效概率，并采用界限法評估整個橋梁系統(tǒng)的抗震性能。

1 非線性有限元模型建立

OpenSEES可對結構進行模態(tài)分析、動力非線性分析以及結構體系在地震作用下的可靠性分析[7]?；贠penSEES建立某連續(xù)梁-剛構橋非線性有限元模型。該橋長6×60 m，橋型如圖1所示。

圖1 算例橋梁計算簡圖 (m)Figure 1 Diagram of the example bridge (m)

主梁等高(3.35 m)，采用C50混凝土。橋面寬度為12.2 m。橋墩1#、4#及5#為分離式菱形實體墩，采用C40混凝土；2#、3#墩采用C50混凝土，為墩梁固接橋墩。縱筋為32 mm帶肋鋼筋，屈服強度335 MPa；箍筋為16 mm光圓鋼筋，屈服強度235 MPa。該橋實際工程所處位置墩底主要為弱風化巖，Ⅱ類場地。所處區(qū)域7度抗震設防，地震分組為第一組。各構件在不同損傷狀態(tài)下用以評估的損傷指標Sc和不確定性βc見表1。

表1 橋梁不同構件的損傷指標及不確定性Table 1 Damage index and uncertainty of different bridge components

采用OpenSEES材料庫中的Concrete 01本構和Steel 02本構分別對混凝土材料和鋼筋材料進行模擬。橋墩非線性特征采用Nonlinear Beam Column模擬，橋梁上部結構采用Elastic Beam Column模擬，橋臺球型支座采用非線性連接單元模擬及ZeroLength Element模擬。依據本課題組前期所做鋼筋混凝土柱試驗[8]建立模型，采用相同加載方式加載，對比模型模擬結果與實際試驗結果如圖2所示。試驗結果與計算結果較為吻合，驗證了建模方法及建模參數的可靠性。

圖2 滯回曲線對比Figure 2 Hysteretic curve comparison

根據算例橋梁所處區(qū)域的地質條件，基于抗震設計規(guī)范[9]生成目標反應譜。依照所生成的反應譜，從PEER強震庫中合理選取地震強度等級均為7級的地震動記錄，用作非線性動力分析的原始地震動。所生成的目標反應譜以及整合后的均值反應譜如圖3所示。

圖3 地震動反應譜 (ξ=0.05)Figure 3 Seismic response spectrum (ξ=0.05)

地震反應譜加速度SA對地震動擬合效果較好[10]。選取SA作為地震動強度指標，以橋梁第一階周期1.012 s所對應的實測SA為基準調整，使其均勻分布于0.1g～1.0g，形成10組共150條不同等級的地震數據，對有限元模型進行非線性動力時程分析。IDA分析中地震動的調幅公式：

ag(i)(t)=kia(t)。

(1)

式中：ag(i)(t)為第i次調幅后地震動；a(t)為原始地震動；ki為調幅系數。

2 關鍵構件易損性分析

分析橋梁各構件在模擬地震動下的易損性，得到在不同強度等級的地震作用下橋臺、支座、橋墩達到各損傷指標狀態(tài)的超越概率。

根據易損性分析理論，對模型輸入所選取的150條地震動并進行非線性動力分析，采用IDA分析方法，調幅地震動，得出計算樣本數據點，用對數線性回歸分析方法整合數據點，得出各構件在地震作用下需求響應值與反應譜SA的對數擬合函數。結合表1給出各構件的損傷指標參數，根據式(2)計算構件失效概率[11]：

(2)

式中：Pf為失效概率；IM為地震動強度參數；βd為結構在地震作用下最大反應的對數標準差，參見式(3)取值[12]：

(3)

分別計算各構件(橋墩、支座、橋臺)在不同強度地震作用下，所產生輕微損傷、中等損傷、嚴重損傷以及完全破壞下的超越概率并繪制曲線，對比分析其易損性曲線，評估橋梁各構件的抗震性能。

2.1 橋墩易損性分析

取各橋墩的150個不同強度地震(0.1g～1.0g)作用下的響應值和對應的地震動反應譜加速度SA，采用對數線性回歸方法，得到各橋墩的地震響應值與SA的對數回歸分析結果，將結果代入式(2)，從而得到各橋墩的易損性函數式。

為計算橋墩的地震易損性損傷指標，通過彎矩-曲率分析，計算橋墩鋼筋屈服時曲率值并計算其等效屈服曲率φ′y及混凝土應變εc=0.004 時的曲率φc4，將曲率值代入并得到橋墩不同損傷狀態(tài)時的損傷指標界限值μcy1、μcy、μc4、μcmax：

(4)

(5)

(6)

(7)

μcmax=μc4+3。

(8)

式中：Δcy為橋墩屈服時墩頂相對位移；Δc4為混凝土截面最外側應變εc達到0.004時墩頂相對位移；Lp為等效塑性鉸長度；L為懸臂橋墩墩高。

所得橋墩損傷狀態(tài)指標如表2所示。1#、3#和5#橋墩易損性曲線如圖4所示。

表2 橋墩損傷狀態(tài)指標Table 2 Damage status index of pier

圖4 橋墩地震易損性曲線Figure 4 Seismic fragility curve of the pier

在不同損傷狀態(tài)下，各橋墩的超越概率隨SA的增大明顯升高；對比各橋墩發(fā)生在輕微損傷、中等損傷、嚴重損傷、完全破壞的易損性曲線，前3種較為接近，且隨著地震強度增大其超越概率在0.2～0.7之間上升較快，0.7之后趨于平穩(wěn)；同一地震強度下，橋墩出現完全破壞的概率明顯低于其他損傷的概率。

為分析橋墩類型對自身易損性的影響，各橋墩不同損傷狀態(tài)下的超越概率曲線見圖5。

圖5 橋墩在不同損傷狀態(tài)下的地震易損性曲線Figure 5 Seismic fragility curves of the pier in different damage states

對于1#與4#橋墩，譜加速度SA=0.6g時，輕微損傷、中等損傷以及嚴重損傷狀態(tài)的超越概率大于80%，完全破壞概率較小，約為20%；5#橋墩在譜加速度SA=0.6g發(fā)生輕微損傷、中等損傷、嚴重損傷以及完全破壞的概率分別為60%、45%、34%、1%；2#和3#橋墩在譜加速度SA=0.6g出現輕微損傷、中等損傷、嚴重損傷概率小于9%，完全破壞概率較小。

由于橋墩延性變形能力與橋墩高度呈現正相關，且1#、4#及5#橋墩混凝土采用C40，比2#與3#橋墩C50強度低，同時由于上部梁體對剛構墩的固結作用，在相同損傷狀態(tài)下，1#與4#橋墩出現損傷概率較大，5#橋墩次之，2#和3#橋墩的概率最小，隨著墩高增加，橋梁變形能力增大，其地震易損性逐漸減小。

由于1#和4#墩高度接近且均為支座墩，易損性曲線基本相同；2#和3#墩高度接近且均為墩梁固結，易損性曲線基本相同。橋墩易損性與材料強度、墩高、墩梁連接方式及墩頂荷載等有關。

2.2 支座、橋臺易損性分析

易損性分析方法與橋墩相同。支座位移損傷指標和橋臺位移損傷指標如表3所示。支座易損性曲線未單獨繪制，將與其他構件共同分析。左、右橋臺易損性曲線見圖6。

表3 支座和橋臺損傷狀態(tài)指標Table 3 Damage status index of support and abutment

圖6 橋臺地震易損性曲線Figure 6 Seismic fragility curves of the abutment

橋臺不同損傷狀態(tài)下的易損性曲線斜率較大，在SA=0.4g時，發(fā)生輕微和中等損傷狀態(tài)的概率均高于90%，在SA=0.6g時，嚴重損傷和完全破壞的超越概率分別為69%、30%；左右橋臺發(fā)生同一損傷的易損性曲線無明顯差別。

3 耐久性損傷橋墩時變地震易損性分析

本文基于橋址位置并參考《混凝土結構耐久性評定標準》(CECS 220—2007)[13]，以5#橋墩為例，分析近海環(huán)境下其氯離子侵蝕致使鋼筋銹蝕的變化規(guī)律，得到不同服役期鋼筋各項參數值，見表4。

表4 不同服役期的鋼筋特征值Table 4 Characteristic values of reinforcement in different service periods

為研究不同服役期鋼筋的銹蝕率對橋墩地震響應的影響，首先計算橋墩在0、20、40、60、80、100 a服役期力學性能，基于已建立的OpenSEES非線性有限元模型，對比分析不同服役期內易損性曲線如圖7所示。

圖7 橋墩時變易損性曲線Figure 7 Time varing fragility curve of pier

結果顯示，在同強度地震作用下，5#橋墩在4種狀態(tài)的易損性隨著服役期增加而逐漸升高。相比于0 a鋼筋未銹蝕的橋墩構件，20、40、60、80、100 a服役期的超越概率在中等損傷狀態(tài)下，分別增大3.1%、11.3%、20.9%、29.1%、34.1%；嚴重損傷狀態(tài)下，分別增大3.9%、13.5%、25.1%、33.4%、35.4%。

4 橋梁系統(tǒng)易損性分析

橋梁結構理論中串聯體系和并聯系統(tǒng)的失效概率[14]可用以下公式計算：

Psys=P(F1∪F2∪F3∪…∪Fn)

(9)

圖8 橋梁系統(tǒng)在不同損傷狀態(tài)下的地震易損性曲線Figure 8 Seismic fragility curves of the bridge system in different damage states

(10)

式中：Psys為整個系統(tǒng)失效的概率；Fi(i=1,2,…,n)表示結構第i個構件失效。

根據各個構件的易損性函數從而建立整個工程結構系統(tǒng)的易損性函數，即界限估計法。一階界限法通過假定橋梁構件失效模式的正相關性 (ρ=0和ρ=1)，評估橋梁系統(tǒng)的最大與最小失效概率，從而確定系統(tǒng)失效概率的上限值和下限值，可表示為[15]

(11)

式中：Pi為第i個構件發(fā)生破壞的概率；m為構件數目。

一階界限法上限值的計算公式為

(12)

根據橋墩、支座和橋臺的地震易損性分析結果，橋梁系統(tǒng)在各損傷狀態(tài)下的易損性曲線界限值如圖8所示，匯總橋梁各構件以及橋梁系統(tǒng)易損性曲線于圖9。

圖9 橋梁系統(tǒng)與各構件在不同損傷狀態(tài)下的地震易損性曲線Figure 9 Seismic fragility curves of the bridge system and each component in different damage states

由于橋梁各構件的最大超越概率決定了系統(tǒng)易損性的下限值，相比于單個構件，橋梁系統(tǒng)超越概率更大，SA=0.6g時發(fā)生嚴重損傷的超越概率已經達到95%以上。若忽略各主要構件對全橋系統(tǒng)抗震性能的影響，以單個構件的易損性來評估橋梁系統(tǒng)的易損性，會高估橋梁抗震性能。

橋梁不同構件的易損性存在較大差異。對于本文實際工程，2#和3#支座最易發(fā)生輕微破壞，其次是橋臺，2#和3#剛構墩發(fā)生輕微破壞的概率最小；中等損傷狀態(tài)下，兩個橋臺易損性最大，其次為1#和4#橋墩，對于2#和3#橋墩以及4#支座最不易發(fā)生中等損傷；1#和4#橋墩最易發(fā)生嚴重損傷，其次是兩個橋臺，SA=0.6g時橋臺發(fā)生嚴重損傷的超越概率為60%，2#、3#橋墩以及4#支座嚴重損傷概率較小，SA=0.6g時，1#、5#支座較2#、4#支座發(fā)生嚴重損傷的超越概率增大25%；兩個橋臺和5#支座完全破壞的概率較大，2#和3#橋墩發(fā)生完全破壞的概率較小。整體而言，該橋梁兩個橋臺相對薄弱，2#、3#橋墩以及4#支座在相同損傷狀態(tài)下的超越概率低于其他構件。

5 結論

(1) 橋墩的抗震性能與材料強度、墩高、墩梁連接方式及墩頂荷載等有關。由于橋墩高度與延性變形能力呈正相關，相同損傷狀態(tài)下，低墩的抗震性能較差。在同座橋梁橋墩高度不一時，應提高低橋墩的抗震性能。

(2) 墩頂支座的易損性受其延性及相鄰支座的影響，由于橋墩與上部結構之間存在位移傳遞效應，橋臺處支座相比于橋墩處支座，其易損性更大。SA=0.6g時，1#、5#橋臺處支座發(fā)生嚴重損傷的超越概率相較于2#、3#橋墩處支座明顯升高25%。

(3) 相同強度地震作用下，橋梁不同構件易損性差異較大。2#和3#剛構墩發(fā)生嚴重損傷的概率最小，對于高墩多跨連續(xù)剛構橋，橋臺易損性較大，而2#、3#高墩及4#支座易損性較小。全壽命周期內耐久性損傷后，橋梁關鍵抗震構件橋墩輕微損傷、中等損傷及嚴重損傷的超越概率顯著增大，完全破壞的超越概率變化較小。

(4) 與單個構件相比，橋梁系統(tǒng)更易發(fā)生損傷。忽略各構件的相關性及其對全橋系統(tǒng)抗震性能的影響，采用單個構件易損性來評估橋梁易損性會高估橋梁抗震性能。