簡支混凝土梁橋Pushover分析簡化方法1

王 草包琦瑋李全旺

1）清華大學土木工程系，北京 100084

2）北京市市政設計研究總院，北京 100035

簡支混凝土梁橋Pushover分析簡化方法1

王 草1)包琦瑋2)李全旺1)

1）清華大學土木工程系，北京 100084

2）北京市市政設計研究總院，北京 100035

靜力Pushover分析方法是研究結構在地震作用下力學性能的非線性分析方法，目前已經得到了廣泛的運用。本文主要討論了對于簡支梁橋的簡化Pushover分析方法。對于非高墩的規則橋梁，結構的重量主要來自于上部結構，因而可以把分析模型簡化為單自由度體系。當給定了墩柱潛在塑性鉸區的彎矩-曲率關系曲線后，可以通過簡單的運算得到其Pushover曲線，從而避免了有限元建模的麻煩。該方法對于工程運用有著一定的參考價值。最后，結合算例表明了該簡化方法的有效性。

簡支梁橋 Pushover分析 抗震評估 簡化方法

引言

在大震作用下，傳統的基于承載力的抗震設計不足以預估結構的變形和結構在大震作用下的實際行為，因此結構基于位移的抗震設計越來越受到重視（毛建猛，2008）。

作為計算結構在地震作用下非線性反應的簡便方法，Pushover方法由于原理簡單、計算方便且精度滿足工程需求，已得到越來越廣泛的研究（秦泗鳳，2008；馬千里等，2008），并已經被國內外的一些規范所采用（中華人民共和國國家標準，2001；FEMA 273，1996）。

在對規則橋梁進行Pushover分析時，當結構的自重主要來自于上部結構時，可以將體系簡化為單自由度體系，這樣，Pushover分析過程便可以得到進一步簡化。

本文首先介紹了靜力Pushover分析方法的原理，然后提出了對于簡支梁橋的Pushover分析簡化方法，最后結合算例說明了該方法的有效性。

1 靜力Pushover分析方法

對于n維多自由度體系，若已知其質量矩陣、阻尼矩陣和剛度矩陣分別為M，C和K，則可以基于D' Alembert原理建立其運動方程（劉晶波等，2005）：

對于一個n維線性空間，n個線性無關的向量可以組成一組基；該空間內的其他任意向量可以用這組基的線性組合來表示。容易證明，對于n個自由度體系而言，其n階主振型{φ1,φ2,…φn}便是其形狀向量的一組基。這樣，式（1）的解（體系的變形形狀向量）可以用φ1，φ2,…φn的線性組合表示：

式中，αi(t)為時變廣義坐標。

在如此的線性表示下，并假定Rayleigh阻尼，則可以利用主振型向量對M和K的正交性，對式（1）進行解耦：

式中，ξj和ωj分別為第j階的主振型對應的廣義阻尼比和頻率；Mj為第j階廣義質量，j=1，2，3，… n。

現考慮由第j振型引起的第i個質點所受的地震力Fij。令，由式（3）可得：

式中，Sa( j )為第j振型對應的加速度譜（Gupta等2000）。

進一步可得：

現對上述n–DOF作如下2點假定：①體系的地震響應僅由基本振型決定；②位移向量x由控制點處的位移xn和結構的基本振型形狀向量φ1表示，即x=xnφ1（認為在地震動過程中，結構的變形形狀保持不變）。

由上述2點假定，令式（5）中的j=1，則：

另一方面，若xn是向量x中最大的元素，則：

由式（7）可得：

對結構施加逐漸增大的水平力F時，可以得到其“水平力–端頂位移”曲線，稱之為Pushover曲線；基于式（6）和式（8）可以將該曲線轉化為“位移–加速度”譜曲線，稱之為“能力譜曲線”。如圖1所示，將能力譜曲線和設計反應譜（需求譜）曲線繪制在同一坐標系中；隨著水平力F的增加（沿著能力譜曲線箭頭走勢），若這2條曲線有交點，則該點所對應的水平位移值xlim即為結構在地震力作用下的最大水平位移響應。然后，若xlim不大于給定的位移閾值，則認為結構在此地震作用下安全；否則，則結構在此地震強度下不安全（Chopra等，2000）。這便是Pushover分析方法的思路。關于xlim的求解方法在眾多文獻中均有介紹，在此不再贅述（毛建猛，2008；馬千里等，2008；Chopra等，2000）。然而需要注意的是，Pushover分析方法是建立在兩個簡化條件的基礎上，而這樣的假定在簡化問題的同時也會帶來一定的計算誤差。由此，眾多學者在研究中都給出了Pushover分析的改進方法（Chopra等，2002；Fajfar等，1996）。由于本文的主旨是給出Pushover分析的一個簡化方法，而非討論該法的精確程度，故而這些改進算法在下文中不予涉及。

圖1 Pushover分析方法圖示Fig. 1 Illustration for the Pushover analysis method

2 對于簡支梁橋的Pushover簡化分析

圖2 獨立墩柱的墩頂水平位移圖Fig. 2 Illustration for the horizontal displacement of a single abutment

基于Pushover分析方法求解簡支梁橋橋墩的頂部位移時，注意到其較之普通高層建筑結構，有著如下特殊性質：①主要質量集中在頂部，而普通高層建筑的質量沿高度方向分散分布；②在墩頂施加水平作用力，可將墩柱看作一個彈（塑）性桿件來研究其變形問題。

考慮如圖2所示的墩柱。當橋墩受到的端頂水平推力較小時，沿墩柱高度方向的彎矩呈線性分布，墩頂處產生的水平位移Δtot在彈性范圍內。當水平推力逐漸增加，于墩底處的彎矩值達到截面的屈服彎矩My時，會在墩底處形成塑性鉸。

如圖3所示，用彈簧模型來模擬塑性鉸的力學特征（秦泗鳳，2008），其中，Lp為等效塑性鉸區的長度。將墩柱的塑性變形等效為Lp長度內的塑性段變形Δpla，則有：

式中，Δela為墩頂處的彈性變形。

對于塑性鉸區的截面，為了簡化問題，通常可以將其彎矩對曲率的曲線（M-φ曲線）雙線性化（中華人民共和國行業標準，2008）。如圖4所示，其中φy為截面的屈服曲率；φu為截面的極限曲率。

圖3 塑性鉸彈簧模型Fig. 3 Spring model for the plastic hinge

圖4 彎矩-曲率的雙折線化模型Fig. 4 Bilinear model for the moment–curvature curve

由圖3則有：

考慮墩頂處的彈性變形量Δela，設曲率沿墩高方向呈線性分布，則有：

將式（10）和式（11b）代入式（9）中，則可以得到在塑性范圍內墩柱的水平位移和塑性鉸區內截面塑性曲率之間的對應關系：

對于多數簡支梁橋或者連續梁橋，認為其結構的主要重量來自于上部結構。在此簡化條件下，對圖2中的墩柱進行受力分析，同時考慮P–Δ效應，則有：式中，M為墩柱底的彎矩；F為墩柱受到的水平作用力；L為墩高；P為上部結構產生的豎向荷載（等效集中作用于墩頂）。

另一方面，如圖4所示，對于塑性鉸區的截面，存在著如下近似的彎矩-曲率關系：

這樣，可以聯立式（12）、（13）和（14）來建立墩柱頂的水平作用力和位移之間的對應關系（即Pushover曲線）：

按前節介紹的Pushover分析步驟，需要將式（15）中的Pushover曲線轉化為能力譜曲線。對于這樣的簡化情況，有：

式中，g為重力加速度。

式（15）和式（16）建立了獨立混凝土墩柱的能力譜曲線，其表達的形式簡單，而且避免了有限元軟件建模分析的麻煩。

如果給出了墩柱精確的M-φ曲線，則可以按照如下方法來得到Pushover曲線：給定曲率值φ，然后由式（12）計算得到Δ（?），再由M-φ曲線給出彎矩M，進而由式（13）解出F。從這里也可以看到，Pushover曲線可以由M-φ曲線通過線性變換得到。

如果考慮墩柱頂的支座剛度，如圖5所示，則可按如下方法來得到Pushover曲線。

更為一般地，如果某一聯橋有多個墩柱，則可以先分別得到每個墩柱的Pushover曲線，然后將這些曲線線性疊加，即可得到整聯的F-Δ曲線。

對于規則圓形或矩形的墩柱截面，由規范可以直接算得M-φ曲線的關鍵參數，從而方便地得到Pushover曲線（中華人民共和國行業標準，2008）；而對于復雜、不規則的截面，可以通過迭代算法（葉列平，2005）給出其M-φ曲線。

3 算例分析

算例引用康釗（2005）對官園橋5號軸橋墩進行順橋向Pushover分析的結果，以該文獻中的計算結果與本文方法給出的結果進行對比。

考慮的橋墩為門式雙墩柱–蓋梁結構，墩柱尺寸如圖6所示，柱高為5m。抗震設防烈度為8度，設計基本地震加速度為0.20g，場地類別為II類。考慮如圖7所示的模型，其中用彈簧模擬支座，其水平抗推剛度為K=45481kN/m。

塑性鉸長度Lp（cm）由規范給出（中華人民共和國行業標準，2008）。

圖5 考慮支座水平位移的Pushover分析Fig. 5 Pushover analysis considering the horizontal displacement of bearing

圖6 墩柱的截面尺寸Fig. 6 Section size of the Guanyuanqiao abutment at Beijing

圖7 Pushover計算模型Fig. 7 Pushover analysis model

其中，b為矩形截面的短邊；H為塑性鉸截面到反彎點的距離（cm）；fy為縱向鋼筋抗拉強度標準值（MPa）；ds為縱向鋼筋的直徑（cm）。

康釗（2005）給出了墩柱潛在塑性鉸區截面的彎矩-曲率曲線，如圖8所示。由此M-φ曲線，結合本文的方法，則可以得到Pushover曲線，如圖9所示。

圖8 單個墩柱的彎矩-曲率曲線Fig. 8 Moment-curvature curve for single abutment

圖9 Pushover分析結果Fig. 9 Result of the Pushover analysis

為了對結構進行抗震性能分析，將Pushover曲線轉化為能力譜曲線，并與需求譜曲線畫在同一坐標系中，如圖10所示。

應注意圖10中采用的是彈塑性反應譜。為了考慮結構的彈塑性，Vidic等（1994）提出用折減系數R(μ)來對彈性反應譜進行修正：

其中，μ為位移延性系數（秦泗鳳，2008），由式（9）和式（11）可得：

R(μ)與μ的關系可表示為：

其中，Tg為結構的特征周期；c1，c2，cR和cT由結構的滯回性質確定，取值見表1。

表1 取ξ = 0.05時，式（21）中的參數取值Table 1 Values for parameters in equation (21) when ξ=0.05

在圖10中從康釗（2005）的結果和本文方法給出的結果均能看到，在8度設防烈度E2地震作用下，支座頂的水平位移為60mm，墩頂的水平位移約為30mm，為墩柱高度的1/167，滿足變形控制要求。同時，從圖10中也可以看到，用本文方法得到的近似Pushover分析方法與非簡化情況相比，二者給出的結果符合的很好。

圖10 確定結構的性能點Fig. 10 Determination of the performance point

4 結論

本文主要討論了簡支梁橋的簡化Pushover分析方法。對于非高墩橋梁，由于結構重力主要來自于上部結構，故而可以將之簡化為單自由度體系。在給定截面的彎矩-曲率曲線后，則通過簡單的運算即可得到結構Pushover曲線，從而省去了有限元建模的麻煩。算例分析結果表明，該方法的精度滿足工程要求，從而對工程運用有一定的參考價值。

康釗，2005．北京市在用市政橋梁基于位移的抗震安全評估（碩士學位論文）．清華大學.

劉晶波，杜修力，2005．結構動力學．北京：機械工業出版社．

馬千里，葉列平，陸新征，2008．MPA法與Pushover法的準確性對比．華南理工大學學報（自然科學版），36（11）：121—128.

毛建猛，2008．Pushover分析方法的改進研究（博士學位論文）．中國地震局工程力學研究所.

秦泗鳳，2008．橋梁抗震性能評價的靜力非線性分析方法研究（博士學位論文）．大連理工大學.

葉列平，2005．混凝土結構（上冊）．北京：清華大學出版社.

中華人民共和國行業標準，2008．公路橋梁抗震設計細則（JTG/T B02-01-2008）．北京：人民交通出版社.

中華人民共和國國家標準，2001．建筑結構抗震設計規范（GB 50011-2001）．北京：中國建筑工業出版社.

Chopra A.K.，Goel R.K.，2000．Evaluation of NSP to estimate seismic deformation：SDF systems．Journal of Structural Engineering，126（4）：482—490.

Chopra A.K.，Goel R.K.，2002．A modal pushover analysis procedure for estimating seismic demands for buildings．Earthquake Engineering & Structural Dynamics，31（3）：561—582.

Fajfar P.，Ga?per?i? P.，1996．The N2 method for the seismic damage analysis of RC buildings．Earthquake Engineering & Structural Dynamics，25（1）：31—46.

FEMA 273，1996．NEHRP commentary on the guidelines for the rehabilitation of buildings．Washington，D.C.：Federal Emergency Management Agency.

Gupta B.，Kunnath S.K.，2000．Adaptive spectra-based pushover procedure for seismic evaluation of structures．Earthquake Spectra，16（2）：367—392.

Vidic T.，Fajfar P.，Fischinger M.，1994．Consistent inelastic design spectra：strength and displacement．Earthquake Engineering & Structural Dynamics，23（5）：507—521.

Simplified Pushover Analysis Method for Simply Supported Concrete Bridges

Wang Cao1), Bao Qiwei2)and Li Quanwang1)

1) Department of Civil Engineering, Tsinghua University, Beijing 100084, China
2) Beijing General Municipal Engineering Design and Research Institute, Beijing 100035, China

The Nonlinear Static Pushover analysis is a method utilized for description of the structural seismic performance under strong earthquakes, and it has been widely used by many researchers. This paper mainly discusses the simplified Pushover method for simply supported concrete bridges. It is noted that for regular bridges with short piers, the structural weight comes mainly from the superstructure, with which the analytical model can be simplified as a single degree system. The Pushover curve can be obtained with simple algebra given the moment –curvature relationship of the section in the potential hinge area, indicating that one does not need to establish a finite element model to get the Pushover curve. Thus, the proposed method can be a valuable reference for practice. Finally, the accuracy of the simplified method is presented with an illustrative example.

Simply supported concrete bridges; Pushover analysis; Seismic safety assessment; Simplified method

王草，包琦瑋，李全旺，2014．簡支混凝土梁橋Pushover分析簡化方法．震災防御技術，9（3）：439—446．

10.11899/zzfy20140310

清華大學自主科研基金（20121087910）

2014-03-11

王草，男，生于1993年。碩士研究生。主要從事橋梁安全評估工作。Email：mornsun17@gmail.com