楊海洋，鐘鐵毅，夏 禾

（北京交通大學(xué)土木建筑工程學(xué)院，北京 100044）

鐵路懸索橋縱向非一致激勵地震響應(yīng)分析

楊海洋，鐘鐵毅，夏 禾

（北京交通大學(xué)土木建筑工程學(xué)院，北京 100044）

以某三跨鐵路懸索橋為例，基于大型通用有限元軟件ANSYS，采用行波輸入與多點激勵兩種非一致地震激勵進行鐵路懸索橋地震響應(yīng)計算分析。研究獲得大跨度懸索橋在非一致激勵輸入下響應(yīng)變化規(guī)律，并通過分析橋梁各響應(yīng)控制振型變化獲得非一致激勵對橋梁地震響應(yīng)影響機理，從而揭示兩種非一致激勵考慮方式對橋梁響應(yīng)影響的根本區(qū)別，為懸索橋非一致激勵下地震響應(yīng)計算分析提供重要參考依據(jù)。

懸索橋；非一致激勵；控制振型；行波效應(yīng)；多點激勵

大跨度橋梁的抗震性能已受到高度重視，并有大量研究［1－4］。對其分析過程中關(guān)鍵為確定地震動輸入。地震動本質(zhì)為空間變化，由于橋梁跨度大、基底延伸長，地基各點振幅、相位差別較大，考慮地震動多點非一致輸入極為重要［5］。

Abdel-Ghaffar等［6］利用實際地震記錄分別在時、頻域內(nèi)對金門大橋進行地震反應(yīng)分析。結(jié)果表明，一致輸入反應(yīng)分析并非代表最不利情況，行波效應(yīng)對懸索橋的地震反應(yīng)影響顯著。Nakamura等［7］用多點激勵反應(yīng)譜方法對金門大橋進行分析，認為受其柔性影響，以動力反應(yīng)分量為主，行波效應(yīng)及多點激勵對動力反應(yīng)分量影響顯著，且多點激勵會致橋塔有小響應(yīng)。Dumanoglu等［8］對三座歐洲大跨度箱梁懸索橋行波效應(yīng)比較分析發(fā)現(xiàn)，在懸索橋的某些部位行波效應(yīng)會增大地震反應(yīng)。Harichandran等［9］認為對懸索橋進行一致輸入分析會高估某些部位反應(yīng)，低估其它部位反應(yīng)。胡世德等［10］通過對江陰長江大橋地震反應(yīng)分析，認為行波效應(yīng)會使結(jié)構(gòu)反應(yīng)增大，而多點激振影響較小。郭發(fā)忠［11］通過對某三跨懸索橋的地震響應(yīng)分析，認為由于多支承激勵引起的柔性結(jié)構(gòu)地震反應(yīng)中擬靜力部分占比較小，總體上行波效應(yīng)對結(jié)構(gòu)內(nèi)力反應(yīng)影響不大。有研究［12－14］只考慮行波效應(yīng)或直接生成人工地震波，雖用兩種方案進行計算，但均未能揭示兩種非一致激勵考慮方式對橋梁響應(yīng)影響的根本區(qū)別。只考慮行波效應(yīng)會高估各點間相關(guān)性，直接生成人工地震波而不進行譜擬合則會致生成的人工地震波反應(yīng)譜值與實際地震波不符，如長周期段譜值偏小等，且各點間地震波反應(yīng)譜值差異亦會造成橋梁地震響應(yīng)變化。

本文以某三跨鐵路懸索橋為例，基于大型通用有限元軟件ANSYS采用行波輸入與多點激勵兩種非一致地震激勵對鐵路懸索橋進行地震響應(yīng)計算分析，獲得該橋的響應(yīng)變化規(guī)律。對歸一化后橋梁地震響應(yīng)時程進行傅里葉變換并分析譜值，獲得橋梁各響應(yīng)控制振型變化，給出非一致激勵對橋梁地震響應(yīng)影響機理，揭示兩種非一致激勵考慮方式對橋梁響應(yīng)產(chǎn)生影響的根本區(qū)別。

1 多點人工地震波生成

將各點隨機相位角項φmk設(shè)為一致，即與各點位置無關(guān)φk，便得局部場地收斂的多點地震波［15］。進行反應(yīng)譜擬合時，本文采用考慮相位譜的人造地震動反應(yīng)譜擬合方法［16］，可得合理的加速度時程曲線，但會在模擬所得地震地面運動中引入長周期分量，造成直接積分的速度、位移時程曲線嚴重失真。本文通過時域內(nèi)基線調(diào)整及頻域內(nèi)過濾掉不合理長周期分量方法達到零線調(diào)整目的。本文加速度時程與某組地震波位移時程對比見圖1。由圖1看出，考慮多點相關(guān)性多點激勵輸入的各支承點位移時程較接近，無明顯延后現(xiàn)象。

圖1 本文生成的加速度時程及各點位移時程示例Fig．1 Acceleration time history and disp lacement time histories

2 計算模型及自振特性分析

本文以擬建某特大鐵路懸索橋為計算分析模型，用ANSYS有限元分析軟件建模。橋梁為地錨式懸索橋，跨度98 m＋660 m＋98 m，中跨設(shè)置吊桿，在跨中設(shè)置三對柔性中央扣，邊跨無吊桿，主塔為混凝土，高196 m、156 m，主梁用上承式鋼桁梁，橋面寬24 m，桁高12 m。主桁與橋塔橫梁相交處只設(shè)豎向、橫橋向支座，順橋向無約束，主梁兩端只約束豎向、橫橋向位移。懸索橋立面圖及各支承點編號見圖2。

圖2 懸索橋立面圖Fig．2 Elevation of the suspension bridge

為簡化計算、明確主梁內(nèi)力變化，橋梁有限元模型用魚骨式建模，見圖3。橋梁前10階振型及描述見表1。橋梁前10階振型見圖4。設(shè)置柔性中央扣后橋梁縱飄振型出現(xiàn)被延后，橋梁1階振型為主梁對稱橫彎。橋梁周期較長，表明該橋?qū)偃嵝越Y(jié)構(gòu)。應(yīng)對橋梁內(nèi)力、位移進行分析研究。

圖3 ANSYS有限元模型Fig．3 ANSYS finite elementmodel of full-bridge

表1 前10階振型及描述Tab.1 Top 10 vibration modes of the bridge

圖4 前10階振型圖Fig．4 Top 10 vibration shape diagram

3 非一致激勵下橋梁地震響應(yīng)分析

本文將系統(tǒng)對大跨度懸索橋在兩種不同非一致激勵方式下地震響應(yīng)進行分析。為反映不同輸入對橋梁地震響應(yīng)統(tǒng)計意義的變化規(guī)律，生成5組多點地震波進行計算，每組地震波含3種波速，分別為500 m／s、1000 m／s及2 000 m／s，每種波速下均含懸索橋的6個支承點地震波。計算一致激勵時用多點輸入的第1點地震波計算，行波輸入則對第1點地震波做延時處理，以保證幾種輸入方式的可比性。采用常遇地震順橋向加載，按50年超越概率63%計算，加速度峰值為0．79 m／s2，計算中考慮幾何非線性影響，并設(shè)地震波從左向右傳播。計算結(jié)果均為用5組地震波分別計算平均值，且均為與自重荷載組合結(jié)果。分析主梁位移時為相對主梁左端相對位移，分析橋塔位移時則為相對各塔底相對位移。

3.1 主纜及吊桿

5組地震波作用、不同輸入方式下主纜及吊桿應(yīng)力平均值見圖5。由表2看出，隨波速變化不同非一致激勵方式下地震響應(yīng)均有不同。行波輸入下隨波速減小主纜及吊桿最大應(yīng)力均有增大，但幅度較小，主纜應(yīng)力最大增大4．82%，吊桿最大應(yīng)力最大增加3．12%；多點激勵作用下主纜及吊桿應(yīng)力與一致輸入基本一致，幅值變化約1%?？紤]本文加載地震動峰值較小，地震造成的響應(yīng)與自重下內(nèi)力相比較小，按罕遇地震加載時各響應(yīng)浮動百分比應(yīng)增大。多種工況作用下吊桿最小應(yīng)力為0值時出現(xiàn)于跨中正中部位，表明地震作用下該位置吊桿可能出現(xiàn)松弛，設(shè)計時應(yīng)采取必要構(gòu)造措施。不同輸入方式的吊桿應(yīng)力隨吊桿位置分布見圖5。由圖5看出，兩種非一致輸入方式的吊桿應(yīng)力分布與一致激勵基本一致，僅應(yīng)力幅值有區(qū)別。

表2 主纜及吊桿應(yīng)力Tab.2 Stress of main cable and suspender

3.2 主梁

不同輸入下的主梁最大地震響應(yīng)見表3。由表3看出，行波輸入下，隨波速的減小主梁最大豎向位移、彎矩及剪力均呈增長趨勢，波速500m／s時增幅分別為403．03%、10．88%及4．37%；而梁端順橋向位移及主梁最大軸力變化規(guī)律不明顯，除波速1 000 m／s有所減小外基本與一致激勵持平，且梁端順橋向位移有減小趨勢；多點激勵作用下僅主梁豎向位移有所增長，波速1 000 m／s時增長100%，但總體隨波速變化規(guī)律不明顯；而梁端位移、主梁最大軸力、最大豎向彎矩及剪力亦與一致激勵基本持平。

為進一步研究兩種不同非一致輸入下主梁的響應(yīng)變化規(guī)律，主梁豎向位移、主梁軸力、主梁豎向彎矩及主梁豎向剪力沿梁長包絡(luò)見圖6～圖9。由4圖看出，不同輸入下主梁響應(yīng)發(fā)生變化區(qū)域均為中跨，兩側(cè)邊跨響應(yīng)基本無影響。行波輸入下主梁跨中豎向位移顯著增大，包絡(luò)形狀由一致激勵的中跨1／4跨、3／4跨雙峰變?yōu)?／2跨單峰，跨中彎矩明顯增大，剪力變化不明顯。多點激勵輸入下主梁跨中豎向位移有一定程度增長，但包絡(luò)仍為雙峰形式，主梁豎向彎矩、剪力包絡(luò)無明顯變化。兩種輸入下主梁軸力分布均與一致激勵相同，區(qū)別僅在于峰值大小。因此，兩種輸入方式對主梁響應(yīng)有較大不同，考慮非一致激勵影響時應(yīng)謹慎選則。

表3 主梁最大地震響應(yīng)Tab.3 Maximum of seism ic responses ofmain girdle

圖5 吊桿應(yīng)力分布Fig．5 Distribution of suspension stress

圖6 主梁豎向位移包絡(luò)Fig．6 Vertical displacement envelope ofmain girdle

圖7 主梁軸力包絡(luò)Fig．7 Axial force envelope ofmain gird le

圖8 主梁豎向彎矩包絡(luò)Fig．8 Verticalmoment envelope ofmain girdle

圖9 主梁豎向剪力包絡(luò)Fig．9 Vertical shearing force envelope ofmain girdle

圖10 主梁端部順橋向位移傅里葉幅值譜Fig．10 Amplitude spectra of longitudinal displacement ofmain girdle

為進一步分析橋梁在非一致激勵作用下響應(yīng)變化機理，本文對主梁各位置響應(yīng)時程歸一化的傅里葉幅值譜計算分析。圖10為兩種非一致激勵對主梁響應(yīng)影響機理區(qū)別的端部順橋向位移幅值譜。由圖10看出，行波輸入下各響應(yīng)幅值譜在頻率點0．253 Hz處譜值均出現(xiàn)不同程度增長，且隨波速減小而加深，而多點激勵作用下未見類似情況。0．253 Hz處為主梁一階豎彎振型，其形狀與行波500m／s的主梁豎向位移包絡(luò)一致。

3.3 橋塔

不同輸入方式下兩側(cè)橋塔最大順橋向地震響應(yīng)見表4。由表4看出，行波輸入下橋塔最大順橋向位移與塔底軸力均有所增加，其中矮塔順橋向位移增加107．28%，塔底軸力增加2．87%，而塔底彎矩及剪力幅值則略有減小，最大減幅分別達27．21%及26．48%，但隨波速增大各項響應(yīng)均越接近一致激勵結(jié)果；多點激勵作用下各項地震響應(yīng)均基本與一致激勵持平，無明顯變化規(guī)律。

表4 橋塔最大順橋向地震響應(yīng)Tab.4 Maxim um of longitudinal seism ic responses of towers

圖11～圖14為橋塔順橋向位移、軸力、順橋向彎矩及順橋向剪力包絡(luò)，由于包絡(luò)線基本對稱，僅給出一側(cè)。由4圖看出，不同輸入下，橋塔順橋向位移包絡(luò)發(fā)生明顯變化：一致輸入時橋塔最大順橋向位移發(fā)生于塔梁固結(jié)處與塔頂間橋塔中部，考慮非一致激勵后塔頂順橋向位移均不同程度增大，矮塔更明顯。橋塔軸力包絡(luò)無明顯區(qū)別，僅幅值有所區(qū)別。橋塔順橋向彎矩（剪力）包絡(luò)在波速500 m／s時變化較大，并仍以矮塔明顯。此矮塔支承點輸入地震波延時較長，且矮塔自振周期較高塔短的緣故。

圖11 橋塔順橋向位移包絡(luò)Fig．11 Longitudinal displacement envelope of tower

圖12 橋塔軸力包絡(luò)Fig．12 Axial force envelope of tower

圖13 橋塔順橋向彎矩包絡(luò)Fig．13 Longitudinalmoment envelope of tower

圖14 橋塔順橋向剪力包絡(luò)Fig．14 Longitudinal shearing force envelope of towe

對橋塔各項地震響應(yīng)傅里葉幅值譜進行分析。矮塔塔底順橋向彎矩幅值譜見圖15。由圖15看出，該結(jié)果與主梁分析結(jié)果一致，橋塔地震響應(yīng)幅值譜仍為行波輸入下0．253 Hz處幅值增加，而多點激勵作用下則無明顯變化。結(jié)果表明，行波輸入對大跨度懸索橋影響主要表現(xiàn)為對一階豎向振型激發(fā)，且能使該振型成為部分響應(yīng)的控制振型，多點激勵對主要振型無顯著影響。其原因為行波輸入下各點間相對位移較大，擬靜力位移導(dǎo)致橋梁一階振型激發(fā)，多點激勵因各點間位移時程相差較小，故無此現(xiàn)象產(chǎn)生。

圖15 矮塔塔底順橋向彎矩傅里葉幅值譜Fig．15 Amplitude spectra of longitudinalmoment of tower

4 結(jié) 論

本文結(jié)合某大跨度鐵路懸索橋，基于大型有限元通用軟件ANSYS，進行兩種非一致激勵下地震動響應(yīng)分析，結(jié)論如下：

（1）行波輸入下橋梁主纜及吊桿應(yīng)力隨波速減小而增加，而多點激勵作用下，主纜及吊桿應(yīng)力與一致輸入基本相同。兩種輸入方式下，吊桿應(yīng)力分布與一致激勵無明顯區(qū)別。

（2）兩種非一致激勵作用時主梁豎向位移較一致激勵均有所增大，行波輸入時增大更顯著，且隨波速減小增幅增大，并伴隨主梁豎向彎矩、位移包絡(luò)形狀改變；多點輸入作用時主梁豎向位移隨波速變化規(guī)律不明顯，豎向彎矩、位移包絡(luò)無明顯改變。兩種非一致激勵作用時主梁豎向剪力峰值及包絡(luò)變化不明顯。

（3）行波輸入時橋塔最大順橋向位移與塔底軸力均有所增大。而塔底彎矩及剪力幅值則略有減小，但波速越快各項響應(yīng)越接近一致激勵結(jié)果；多點激勵作用下各項地震響應(yīng)均與一致激勵基本持平，無明顯變化規(guī)律。兩種非一致激勵作用時塔頂位移均顯著增大。

（4）行波輸入對大跨度懸索橋影響主要表現(xiàn)為對一階豎向振型激發(fā)，而多點激勵對主要控制振型無顯著影響。各點間相對位移較大，擬靜力位移致橋梁一階振型激發(fā)，多點激勵由于各點間位移時程相差較小。

（5）不同輸入下懸索橋地震響應(yīng)的主要控制振型未發(fā)生改變時其響應(yīng)幅值變化不會太大。進行懸索橋地震響應(yīng)分析前應(yīng)先分析各響應(yīng)主要控制振型，以便選擇合理的地震動輸入方式。

（6）本文計算結(jié)果僅考慮常遇地震，所得橋梁地震響應(yīng)變化幅度僅作規(guī)律研究。罕遇地震作用下響應(yīng)變化幅度會更明顯，設(shè)計中應(yīng)注意。

［1］李吉濤，楊慶山，劉陽冰．多點地震激勵下大跨連續(xù)鋼構(gòu)橋易損性分析［J］．振動與沖擊，2013，32（5）：75－80．

LIJi-tao，YANG Qing-shan，LIU Yang-bing．Fragility analysis of long span continuous rigid frame bridge undermulti-support excitations［J］．Journal of Vibration and Shock，2013，32（5）：75－80．

［2］柳國環(huán)，李宏男，田利．九江長江大橋在多點多維地震激勵下的反應(yīng)分析［J］．振動與沖擊，2009，28（9）：224－230．

LIU Guo-huan，LIHong-nan，TIAN Li．Response analysis of Jiujiang Yangtze river highway bridge under spatially variable earthquake ground motions［J］．Journal of Vibration and Shock，2009，28（9）：224－230．

［3］Lin JH，Zhang Y H，Li Q S，et al．Seism ic spatial effects for long-span bridges using the pseudo excitation method［J］Engineering Structures，2004，26（9）：1207－1216．

［4］劉洪兵，朱晞．大跨度斜拉橋多支承激勵地震響應(yīng)分析［J］．土木工程學(xué)報，200l，34（6）：38－44．

LIU Hong-bin，ZHU Xi．Seismic response analysis of longspan cable-stayed bridges under multi-support excitations［J］．China Civil Engineering Journal，2001，34（6）：38－44．

［5］苗家武，胡世德，范立礎(chǔ)．大型橋梁多點激勵效應(yīng)的研究現(xiàn)狀與發(fā)展［J］．同濟大學(xué)學(xué)報，1999，27（2）：189－123．

MIAO Jia-wu，HU Shi-de，F(xiàn)AN Li-chu．State of art and development of study on seismic performance of long span bridges under multiple-support excitations［J］．Journal of Tongji University，1999，27（2）：189－123．

［6］Abdel-Ghaffar A M，Rubin L I．Suspension bridge response to multiple-support excitations［J］．Journal of Engineering Mechanics，ASCE，1982，108（2）：419－435．

［7］Nakamura Y，Kiureghian A D，Liu D V．Multi-support response analysis of the golden gate bridge［R］．University of California atBerkely，Report No．UCB／EERC－93／05，1993．

［8］Dumanoglu A A，Brownjohn JMW，Severn R T．Seismic analysis of the fatih sultan mehmet suspension bridge［J］．Earthquake Engineering and Strutural Dynamics，1992，21：881－906．

［9］Harichandran R S，Hawwari A，Sweidan B N．Response of long-span bridges to spacially varying ground motion［J］．Journal of Structural Engineering ASCE，1996，122（5）：476－484．

［10］胡世德，范立礎(chǔ)．江陰長江公路大橋縱向地震反應(yīng)分析［J］．同濟大學(xué)學(xué)報，1994，22（4）：433－438．

HU Shi-de，F(xiàn)AN Li-chu．The longitudinal earthquake response analysis for the Jiangyin Yangtze river bridge［J］．Journal of Tongji University，1994，22（4）：433－438．

［11］郭發(fā)忠．大跨度懸索橋行波激勵地震反應(yīng)分析［J］．公路交通技術(shù)，2011（3）：60－62．

GUO Fa-zhong．Analysis for seisimic response of traveling wave excitation of large-span suspension bridges［J］．Technology of Highway and Transport，2011（3）：60－62．

［12］張卉，王志清，彭元誠，等，四渡河大跨懸索橋空間地震響應(yīng)分析［J］．地震工程與振動，2007，27（3）：70－76．

ZHANG Hui，WANG Zhi-qing，PENG Yuan-cheng，et al．Spatial seismic response analysis of long-span suspension bridge across Sidu river［J］．Earthquake Engineering and Engineering Vibration，2007，27（3）：70－76．

［13］孫永濤．大跨懸索橋在多點激勵作用下的地震反應(yīng)分析［D］．哈爾濱：哈爾濱工業(yè)大學(xué)，2009．

［14］陳海波．大跨度懸索橋的地震反應(yīng)分析［D］．合肥：合肥工業(yè)大學(xué)，2010．

［15］董汝博，周晶，馮新．一種考慮局部場地收斂性的多點點地震動合成方法［J］．振動與沖擊，2007，26（4）：5－9．

DONG Ru-bo，ZHOU Jing，F(xiàn)ENG Xin．A local convergent method for simulating multi-point earthquake ground motion［J］．Journal of Vibration and Shock，2007，26（4）：5－9．

［16］胡聿賢，何訓(xùn)．考慮相位譜的人造地震動反應(yīng)譜擬合［J］．地震工程與工程振動，1986，6（2）：37－51．

HU Yu-xian，HE Xun．Phase angle consideration in generating response spectrum-compatible groundmotion［J］．Journal of Earthquake Engineering and Engineering Vibration，1986，6（2）：37－51．

Seism ic responses analysis of a railway suspension bridge under longitudinal non-uniform excitations

YANGHai-yang，ZHONG Tie-yi，XIAHe
（School of Civil Engineering，Beijing Jiaotong University，Beijing 100044，China）

Based on the general finite element software ANSYS，seismic responses under non-uniform excitations of a long-span railway suspension bridge were studied．Two different inputmethods，i．e．，traveling seismic wave excitation and multi-support excitation，were used to analyze the seismic responses of the bridge，and then the response regulations were discussed．Analysing the differences between the response control modes，the influence mechanisms of the nonuniform excitations were discussed，thus the fundamental differences between effects of the two input methods were revealed．The results can provide an important basis for the seismic analysis of suspension bridges．

suspension bridge；non-uniform excitation；controlmode；travelling wave effect；multi-support excitation

U24

：A

10．13465／j．cnki．jvs．2014．22．028

鐵道部科技研究開發(fā)計劃（2010G004－I）；國家重點研究計劃“973”項目（2013CBO36203）

2013－11－08 修改稿收到日期：2014－01－15

楊海洋男，博士生，1991年6月生

鐘鐵毅男，博士后，副教授，1963年4月生