【摘要】近幾年來，國內(nèi)外學者對大跨徑橋梁的靜風穩(wěn)定性展開了廣泛的研究，相繼產(chǎn)生了一些分析方法。本文首先對現(xiàn)有的大跨徑橋梁靜風穩(wěn)定性分析方法進行了回顧與比較，并就其中的不足進行簡要討論，作為改進提出了一種新方法，采用這些方法對虎門大橋的靜風穩(wěn)定性進行了計算與比較。
關(guān)鍵詞 大跨徑橋梁 靜風穩(wěn)定性 非線性分析

一、引言
近幾年來，我國先后建成了主跨888m的虎門懸索橋和主跨1385m的江陰長江懸索橋以及主跨423m的南浦斜拉橋和602m的楊浦斜拉橋。交通部規(guī)劃將興建的沿海高等級公路干線上，還有五座大型跨海工程，它們自北向南依次跨越渤海海峽、長江口、杭州灣、珠江口伶仃洋以及瓊州海峽，這些跨海工程都需要建造大跨徑橋梁，為大跨徑橋梁提出了更高的要求。
但是，隨著橋梁跨徑的不斷增大，勢必會帶來一些新的問題。風荷載作用下大跨徑橋梁的靜力穩(wěn)定問題就是其例。早在1967年日本東京大學Hirai教授就在懸索橋的全橋模型風洞試驗中觀察到了靜力扭轉(zhuǎn)發(fā)散的現(xiàn)象，同濟大學風洞實驗室在對汕頭海灣二橋的風洞試驗中，也發(fā)現(xiàn)了斜拉橋由靜風引起的彎扭失穩(wěn)現(xiàn)象[1，2]。最近，Boonyapinyo， Miyata、謝旭等學者通過計算也表明了這種現(xiàn)象出現(xiàn)的可能性[1,3]。圖1為大跨徑橋梁結(jié)構(gòu)靜風失穩(wěn)過程示意圖。從圖中我們可以看出接線性方法進行大跨徑橋梁的靜風臨界風速的求解一般會過高地估計了橋梁的抗風能力，是偏于不安全的[4]。而采用僅考慮結(jié)構(gòu)幾何非線性分析方法計算出的先穩(wěn)風速也會偏高，只有全面考慮結(jié)構(gòu)幾何非線性與靜風荷載非線性的分析方法，才能比較真實地反映橋梁結(jié)構(gòu)失穩(wěn)的全過程。因此，有必要對大跨徑橋梁的非線性靜風穩(wěn)定性的分析方法進行全面的研究。

本文首先回顧和比較了目前幾種主要的靜風穩(wěn)定性分析方法，針對這幾種方法的不足，提出了一種新的求解大跨徑橋梁靜風穩(wěn)定性的方法--增量與內(nèi)外兩重迭代法，最后，通過一座實橋算例分析對這些方法進行了比較分析，得出了一些結(jié)論。

二、現(xiàn)有靜風穩(wěn)定性分析方法的回顧與比較
1．靜風穩(wěn)定性分析方法的回顧
早期的靜風穩(wěn)定性分析方法主要是線性方法，該方法根據(jù)結(jié)構(gòu)失穩(wěn)方式不同又可分為側(cè)傾失穩(wěn)和扭轉(zhuǎn)發(fā)散兩種[5]。下面以單跨懸索橋為例分別給出這兩種失穩(wěn)模態(tài)的計算公式。
（1）側(cè)傾失穩(wěn)臨界風速計算
Vlb＝KlbftB（ 1）
式中

其中，Cd為加勁梁阻力系數(shù)；Cl為加勁梁升力系數(shù)的斜率；ε為扭彎頻率比；Bc為主纜間距；h為加勁梁梁高。
(2)扭轉(zhuǎn)發(fā)散臨界風速計算
vtd＝KtdftB
式中

其中， m為單位橋長質(zhì)量（kg／m）；此為單位橋長質(zhì)量慣矩 （kg·平方米／m）； CMO為在攻角α=0度時升力矩系數(shù)CM的斜率，由風洞試驗測得；ft一般取一階對稱扭轉(zhuǎn)頻率（Hz）。
由于以上兩種公式均未考慮結(jié)構(gòu)非線性和靜風荷載非線性因素的影響，因此計算出的靜風臨界風速一般會明顯偏高，對結(jié)構(gòu)而言是偏于不安全的。
文獻［ 6」提出一種用于分析大跨徑懸索橋靜風穩(wěn)定性的實用方法--三角級數(shù)法，該方法綜合考慮了結(jié)構(gòu)幾何非線性和靜風荷載升力和升力矩共同作用的非線性影響。該方法是由計算某一風速下結(jié)構(gòu)的靜風響應和結(jié)構(gòu)臨界風速的計算兩部分組成。其中在計算結(jié)構(gòu)靜風響應時，將升力和升力矩曲線按分段直線擬會，將升力、升力矩、豎向位移和結(jié)構(gòu)扭轉(zhuǎn)角用一組三角級數(shù)表示，并分別代入到懸索橋豎向和扭轉(zhuǎn)平衡微分方程，聯(lián)立確定各級數(shù)項的待定系數(shù)，由于方程是非線性的，所以此項計算必須通過迭代方法完成。計算扭轉(zhuǎn)發(fā)散臨界風速時，在初始攻角下，先假定一初始風速V0。，再通過靜風響應計算得到一組相應的主纜索力和扭轉(zhuǎn)角，以此為新的初態(tài)增加一級風速重新計算，當前后兩次所計算的扭轉(zhuǎn)角的相對誤差超出允許值時，認為結(jié)構(gòu)出現(xiàn)扭轉(zhuǎn)發(fā)散，此時計算出的風速即為臨界風速。
文獻[4]采用增量法和迭代法相結(jié)合的方法進行大跨徑橋梁第二類靜風穩(wěn)定性有限元分析。其中增量法分為內(nèi)增量和外增量兩種，所謂內(nèi)增量，就是采用增量法進行結(jié)構(gòu)幾何非線性的求解，外增量就是不斷施加風速。而迭代法則是進行某一風速下結(jié)構(gòu)平衡狀態(tài)的求解。該方法雖然能較好地跟蹤結(jié)構(gòu)失穩(wěn)的全過程，但是由于它采用增量法進行結(jié)構(gòu)幾何非線性的計算，因此會出現(xiàn)計算誤差累計的問題。
2．靜風穩(wěn)定性分析方法的比較
以上分別介紹了三種不同的大跨徑橋梁靜風穩(wěn)定性分析方法，為了更好地認識和了解它們，本文對其作了比較分析，如表1所示。

三、對靜風穩(wěn)定性分析方法的改進
針對上述幾種靜風穩(wěn)定性分析方法的不足，本文在綜合考慮靜風荷載與結(jié)構(gòu)非線性影響的基礎(chǔ)上，采用增量與內(nèi)外兩重選代相結(jié)合的方法，實現(xiàn)了對大跨徑橋梁靜風穩(wěn)定性的精確求解。
作用在主梁單位長度的靜風荷載可分解為橫向風荷載PH、豎向風荷載PV和扭轉(zhuǎn)力矩M，具體表達式如下：

式中，CH（a），CV（a），CM（a）分別表示在有效攻角下主梁沿結(jié)構(gòu)坐標軸各方向的阻力、升力、升力短系數(shù)；所謂有效攻角，是指靜風初始攻角與靜風作用引起的主梁扭轉(zhuǎn)角之和。
按照桿系結(jié)構(gòu)空間穩(wěn)定理論，問題可歸結(jié)為求解如下形式的非線性方程：

式中，Ke和Kg分別為結(jié)構(gòu)的線彈性和幾何剛度矩陣；α為有效攻角。PH，PV，PM分別為體鈾下的風阻力、升力和升力矩；f為函數(shù)因子；上標G和W分別代表重力和風力。
從（4）式可知，不僅結(jié)構(gòu)的剛度是結(jié)構(gòu)變形的函數(shù)，而且右端項所表示的靜風荷載也是結(jié)構(gòu)變形的函數(shù)，為了求解該非線性方程，就必須采用迭代法。而為了跟蹤結(jié)構(gòu)變形的全過程，又必須采用增量法。為此，本文提出了采用增量與內(nèi)外兩重選代相結(jié)合的方法。增量法將風速按一定比例增加。而其中的內(nèi)層選代主要是進行結(jié)構(gòu)的非線性計算，而外層選代則是為了尋找結(jié)構(gòu)的某一風速下的平衡位置。該方法的具體實施步驟如下：
（l）假定一初始風速V0；
（2）計算在該風速下結(jié)構(gòu)所受的靜風荷載；
(3）采用 Newton- Rapson法求解（4）式，得到結(jié)構(gòu)位移U；
（4）從結(jié)構(gòu)位移U中提取單元扭轉(zhuǎn)角（為左右兩節(jié)點扭轉(zhuǎn)位移之和的平均值），重新計算結(jié)構(gòu)的靜風荷載；
（5）檢查三分力系數(shù)的歐幾里得范數(shù)是否小于允許值；
（6）如果小于允許值，則按預定步長增加風速，重復步驟（2）～（5）；否則，重復步驟(3)～(5)；
(7)如果在某一級風速下，出現(xiàn)選代不收斂，則恢復到上一級風速狀態(tài)，縮短步長，重新計算，直至相鄰兩次風速之差小于預定值為止。

四、算例分析與比較
分別采用上述幾種不同的方法對虎門大橋進行了靜風穩(wěn)定性分析。分析中采用的的結(jié)構(gòu)幾何和材料參數(shù)、各攻角下的靜力三分力系數(shù)均取自同濟大學土木工程防災國家重點實驗室--風洞實驗室的資料和試驗報告[7]。計算出的臨界風速列于表 2。圖 2為采用不同方法計算出的結(jié)構(gòu)主梁跨中斷面橫向位移隨風速的變化（因側(cè)傾失穩(wěn)和扭轉(zhuǎn)發(fā)散兩種方法無法跟蹤結(jié)構(gòu)失穩(wěn)的全過程，故未列出）。為了更好地說明本文所提的增量一內(nèi)外兩重迭的代法的穩(wěn)定性，分三種不同的初始風速加載方式考察了結(jié)構(gòu)靜風響應，其中方式一為風速從40m／s開始，每次增加 10m／s，直至110m／s；方式二為風速從70m／s開始，每次增加20m／s，直至 110m／s；方式三為一次性加至110m／s，計算結(jié)果如表3所示。

五、結(jié)論
通過對表1中所列的幾種主要大跨徑橋梁靜風穩(wěn)定性分析方法的回顧與比較，可以得出以下三點結(jié)論：
（1）采用線性方法來計算大跨徑橋梁的靜風穩(wěn)定性，一般會過高地估計結(jié)構(gòu)的抗靜風能力，偏于不安全，這與它們未考慮結(jié)構(gòu)幾何非線性和靜風荷載非線性有關(guān)。而采用三角級數(shù)法、增量迭代法和本文提出的改進方法--增量一兩重選代法，均能比較真實地反映結(jié)構(gòu)的實際抗靜風能力。
（2）線性方法是通過預先假定結(jié)構(gòu)的失穩(wěn)模式來求解結(jié)構(gòu)靜風失穩(wěn)風速，這將導致其無法跟蹤結(jié)構(gòu)失穩(wěn)的全過程，而其余三種方法則不同，它們均是通過風速的逐級加載來獲得結(jié)構(gòu)靜風；臨界風速，因此，它們完全能跟蹤結(jié)構(gòu)失穩(wěn)的全過程。
（3）增量迭代法是采用增量法進行結(jié)構(gòu)幾何非線性的求解，因此，勢必會造成計算誤差的積累，此外，如果要獲得結(jié)構(gòu)的靜風臨界風速，將只能進行風速的逐級加載，這勢必又會造成計算速度慢的缺點。而本文提出的改進法--增量一兩重迭代法無須進行風速的逐級加載就能獲得結(jié)構(gòu)靜風臨界風速，且有比較高的精度，這一點在算例中得到了充分的驗證。

參考文獻
［ l］V．Boonyapinyo,H.Yamada,T.Miyata,Wind-induced nonlinearl lateral-torsional buckling of cable-stayed bridges, Jouynal of Structural Engineering,ASCE,120(2);486-506，1994
[2]項海帆，林志興．橋梁抗風設(shè)計規(guī)范》的研究課題．結(jié)構(gòu)工程師（增刊），1998

［3］Nagai，M.，Xie X．，and etc， static and dynamic instablility analysis of 1400－Meter long-span cable－ stayed bridges，IABSE"Long－Span and High-rise Structures",1998
［4］《公路橋梁抗風設(shè)計指南）編寫組．公路橋梁抗風設(shè)計指南．北京：人民交通出版社，1996
［5］程進，肖汝誠，項海帆．大跨徑懸索橋靜風扭轉(zhuǎn)發(fā)散的級數(shù)解法．同濟大學學報，1992（2）：127-130