多次地震動下鐵路圓端空心墩的位移響應研究

王 猛

(西南交通大學土木工程學院橋梁系，四川成都 610031)

隨著國內抗震規范的修訂，基于結構性能的抗震設計理念也越來越深入人心[1-3]。相關資料[4-6]表明，由于區域、地質結構、應力場強度等方面的不同，在地震類型中，主余震型占60 %，震群型占25 %，孤立型地震序列僅占15 %。整體上，結構在使用期間遭遇的地震動次數往往較多。有工程實例表明，個別結構在主震中并未發生倒塌現象，卻在隨后的余震中損壞[7-8]。因此有很多學者致力于研究多次地震動作用對結構的影響，如：Mhin[9]于20世紀80年代，通過研究非線性單自由度體系對主余震的響應情況進行研究，認為主余震要求結構有更高的延性性能；Amadio.C等[10]學者同樣通過單自由度體系來研究多次地震動下的結構損傷情況，認為多次地震會加劇結構的損傷；Sunasaka[11]等學者對城市高架橋進行研究，認為后續地震會加劇結構損傷，進而導致結構破壞。

國內方面，張沛洲、康謹之等[12]學者以單自由度為例分析主余震對結構累積損傷及抗震性能的影響；于曉輝、呂大剛等[13]學者通過對比計算結構在考慮與不考慮多次地震動下的損傷差值，來研究其對結構的損傷影響；李洪泉等[14]學者通過試驗研究分析結構在地震作用下的累計損傷情況；張煜敏[15]等學者以連續梁為例研究多次地震序列對橋梁結構產生的潛在震害。

以上學者研究均表明，在多次地震作用下，橋梁結構響應往往更為劇烈，導致的損傷也隨之增加，在主震作用下未倒塌的結構在隨后的余震中存在倒塌風險。但以上研究內容多數以單自由度為主，較少涉及多自由度結構，對山區鐵路橋墩常用結構形式圓端空心墩也并未提及，再綜合現行鐵路抗震設計規范[16]并未考慮多次地震動(不同于頻遇地震)這一現象對橋梁結構的影響，作者以秋末河大橋為工程背景，基于抗震分析軟件OpenSees對鐵路圓端空心墩進行傳統及累計增量動力分析(IDA)，對比同一結構在單一地震及多次地震動下位移響應情況，來研究多次地震對鐵路空心高墩的位移影響規律，進而給出相關抗震分析計算建議。

1 工程概況

秋末河大橋位于成都至蘭州鐵路段內，為一段曲線梁橋，前連成都，后接蘭州，其復雜的水文地質條件給設計、施工帶來諸多不便。該橋橋墩普遍較高，在30.5～42.5m之間，設計時采用變截面圓端形式，外部變化比例為40∶1，內部空心部分變化比例為70∶1，橋墩上下端部為變截面圓端實心墩，端部位置設置倒角，橋墩大樣如圖1所示，主要截面大樣如圖2所示。

圖1 橋墩構造

圖2 橋墩橫截面

2 模型建立

2.1 Midas全橋及單墩模型

采用橋梁通過有限元軟件Midas/civil建立秋末河大橋全橋有限元模型(圖3)。上部主梁、下部橋墩、樁基礎采用梁單元模擬，其中樁—土相互作用通過在樁基礎上施加等效節點彈簧來模擬，彈簧剛度采用“m”法進行計算。

圖3 Midas/Civil有限元模型

由于秋末河大橋橋墩截面形式相同，均為變截面圓端橋墩，全橋模型分析計算較為繁瑣，為提高計算效率，取一個具有代表性的橋墩進行研究，所選取的橋墩設計參數如表1所示。通過全橋有限元模型計算出主梁作用于單墩墩頂的豎向力，隨后建立Midas/Civil單墩有限元模型，計算其動力特性，累積質量參與系數隨振型變化情況如圖4所示。

圖4 累積質量參與系數分布

表1 結構設計參數

從圖4的變化情況來看，當順橋向累積質量參與系數達到90 %以上時，振型數量為17，并非單一振型控制結構，該圓端空心墩結構較柔，剛度相對較小，和以往研究的低階振型為主的結構有所不同。

2.2 OpenSees單墩模型

由于OpenSees處理抗震過程中涉及到的非線性能力較強，因而廣泛應用于結構抗震分析中。本次模擬內容包括以下幾方面：材料的本構關系、單元類型的選取、結構單元的離散、結果輸出。模型參數見表1所示，本模型所涉及材料為鋼筋和混凝土。對于鋼筋采用OpenSees中自帶的Steel01理想彈塑性材料；對于混凝土材料，將其分為強約束區、弱約束區與非約束區，分別采取不同參數、相同本構類型的Concrete02材料進行模擬，其本構曲線如圖5所示。由于本橋墩高階振型較為明顯，因此橋墩單元數量取為30。建立單墩結構后，將全橋模型中計算出的豎向力作用于單墩墩頂，并將其轉化為質量進行數值模擬。

圖5 混凝土本構曲線

2.3 地震波的選取

當地地質資料表明，秋末河大橋地震動反應譜特征周期為0.45s，結合本文所研究內容為多次地震作用對鐵路圓端空心墩的位移響應影響規律，所選取地震波均為實際記錄地震波，分別為ElCentr波、唐山遷安波、汶川綿竹波。

3 計算分析

3.1 累積增量動力分析

傳統增量動力分析(IDA)[17]是指，對結構作用一系列相同變化規律、不同峰值加速度的地震動，通過分析結構的響應變化規律，來研究結構的抗震性能，通常選用某個特定的指標進行評價，如位移延性、曲率延性、位移角等。該方法首先由Bertero[18]提出，隨后被美國相關規范所采納[19-20]，該方法發展至今，已逐漸被人們作為一種評價結構整體抗側傾能力的常用手段。

雖然該方法被廣泛接受，但其對考慮結構在多次地震動下的響應情況尚有欠缺。對此，作者提出累積增量動力分析方法，來考慮結構在多次地震動下的響應。具體加載方式如圖6所示，在同一次分析工況中，同時作用一系列相同規律、不同峰值的地震動，為了使結構充分振動，經試算后，在每個地震動后增加30s的間隔期。

圖6 累積增量動力分析加載方式

3.2 計算結果

結構在傳統增量動力分析下的位移情況在后續對比數據中給出；結構在累積增量動力分析下的位移時程曲線如圖7～圖9所示。將同一地震波下不同分析方法得到的位移結果進行對比，如圖10～圖12所示，兩種分析方法作用下，絕對最大位移比值、殘余位移比值情況見表2～表4所示。

圖7 El Centro地震波累積IDA結果

3.3 結果分析

從傳統IDA計算結果上看，結構絕對最大位移及殘余位移隨著峰值加速度的增加而增加，其原因較為明顯，對于最大位移響應而言，地震波加速度增加，結構的慣性力也隨之增加，結構最大位移響應也隨之增加；對于殘余位移而言，其變化情況需分兩種情況分析：當結構處于彈性范圍內時，殘余位移隨著加速度變大而變大(近似線性增加)，當結構處于彈塑性范圍內時，殘余位移隨著加速度的變大而變大，且速度較彈性明顯，其原因在于結構進入彈塑性后，結構會產生相應新的塑性位移。

表2 El Centro地震波作用下結構位移

表3 唐山遷安地震波作用下結構位移

表4 汶川綿竹地震波作用下結構位移

圖8 唐山遷安地震波累積IDA結果

圖9 汶川綿竹地震波累積IDA結果

圖10 位移對比—El Centro

圖11 位移對比—唐山遷安

圖12 位移對比—汶川綿竹

從累積IDA計算結果上看，結構絕對最大位移及殘余位移同樣隨著峰值加速度的增加而增加。對于累積分析，結構的最大位移及殘余位移較傳統分析更大，原因在于：在進行累積IDA計算時，較高地震波工況的開始狀態為較低地震波工況的結束狀態，繼承了上一工況后的結構狀態。若上一工況結束狀態存在殘余位移，則下一工況的殘余位移會更大；若上一工況結束時存在損傷、材料剛度退化等不利結構受力的狀態，則下一工況結構會發生更大的損傷等現象，故而結構的位移響應會更大。

將兩種分析方法再次進行對比，從計算對比圖(圖10～圖12)、對比表(表2～表7)所示內容，可以發現，累積IDA作用下的結構位移響應存在“跑偏”現象，即位移向某一方向增加；同時，對于結構最大位移響應而言，當結構處于彈性狀態，其最大位移響應的誤差在5 %以內，在工程可接受誤差范圍內；當結構處于彈塑性狀態，最大位移響應差異較大。對于結構殘余位移而言，無論結構處于彈性還是彈塑性狀態，其差異均很大，不容忽視。

表5 結構位移對比—El Centro

表6 結構位移對比—唐山遷安

表7 結構位移對比—汶川綿竹

通過以上對比分析可知，多次地震動對鐵路圓端空心墩的位移響應有較大的影響，主要原因在于前一地震工況產生的結構剛度退化、材料損傷及殘余位移對后一地震工況有較大的影響。本文涉及的是圓端空心墩，結構高階振型較為明顯，不再是低階振型控制結構響應的結構，希望以后再進行相關結構的抗震位移性能分析時，能夠參考本文的相關分析結果。

4 結論

(1)本文以秋末河大橋為工程背景，通過Midas/civil建立全橋模型計算后簡化出單墩模型，通過OpenSees對圓端空心墩進行地震作用下位移性能分析。

(2)該圓端空心墩結構較柔，高階振型較為明顯，并非低階振型控制的結構，在進行單一地震動計算時，即進行傳統IDA，結構最大位移及殘余位移隨峰值加速度的增加而增加。

(3)在考慮多次地震動對結構影響情況下，即進行累計IDA計算，結構的位移響應與不考慮該影響時位移響應有較大區別，當結構處于彈性時，最大位移響應差異不大，誤差在5 %以內，為工程可接受誤差；當結構處于彈塑性時，最大位移響應差異較大，不可忽略；且隨著彈塑性的發展，多次地震動對最大位移的影響也更為明顯。

(4)多次地震動對圓端空心墩結構的殘余位移影響很大，即使結構處于彈性狀態，其差異也較大，不容忽視；隨著結構彈塑性的發展，多次地震動對其影響(即兩種情況下的殘余位移比值)也趨于穩定。