上承式拱橋立柱在超速汽車下的動態(tài)效應(yīng)

吳隆 李海濤 王彥文

摘要：為了分析上承式拱橋立柱在移動荷載下的車一橋耦合效應(yīng)，在考慮拱肋的內(nèi)力的前提下通過理論分析，利用有限元方法建立起拱橋、車輛的模型，通過試驗監(jiān)測數(shù)據(jù)和數(shù)值分析結(jié)果分析汽車的超速現(xiàn)象對橋梁立柱應(yīng)力、位移及橋梁變形的影響。結(jié)果表明車輛的超速行駛對立柱的軸力產(chǎn)生明顯的影響，立柱在z方向剪力會產(chǎn)生應(yīng)力峰值，同時超速現(xiàn)象對立柱y方向剪力及橋梁的變形并沒有直接影響。

關(guān)鍵詞：橋梁工程;拱橋;立柱;有限元方法;拱肋內(nèi)力;超速

中圖分類號：U448.2

文獻標(biāo)志碼：A

文章編號：1672-1098（ 2020）04-0076-07

作者簡介：吳?。?996 -），男，安徽黃山人，在讀碩士，研究方向：橋梁工程。

拱橋作為一種特殊且常見的橋梁形式，由于其特殊形式及特殊受力形式，眾多學(xué)者對其在各種荷載作用下的動力響應(yīng)及損傷做了很多研究。Sher-van Ataei[1]等通過對動態(tài)放大系數(shù)（ DAF）的分析計算，研究了使用壽命較長的石拱橋在不同軸載列車作用下的動態(tài)響應(yīng)。M.R.Akbari[2]等用豎向簡諧荷載作用下的拱形梁對拱橋的非線性振動進行了建模和仿真，并通過Akbari-ganji method（ AGM）對其進行了分析，研究而拱形梁的非線性振動特性。Vicente N.Moreira[3]等運用塑性理論，分析了葡萄牙的5座建成年代較久的石拱橋承受高于所設(shè)想的動態(tài)荷載下的承載情況和可用性性能。H-Gonen[4]等研究了帕魯橋在恒載下的變形和應(yīng)力，同時將地震加速度作為荷載對該橋進行了數(shù)學(xué)建模及動力分析和計算。Abdulkadir Cnneyt Aydin[5]等研究了奧杜省的薩普戴爾橋在一定荷載作用下的受力性能和在不同位置受到荷載作用下的性能。

由于車輛性能的不斷改進及駕駛?cè)藛T的疏忽，車輛超速行駛現(xiàn)象已屢見不鮮?；诖似洮F(xiàn)象，本文對中國貴州省福泉市三江口拱橋考慮在受到相應(yīng)彎矩、軸力和剪力的作用下，在移動荷載作用下立柱的車一橋耦合效應(yīng)進行了理論分析、模擬和試驗，研究了立柱在超載和超速下的動態(tài)響應(yīng)。

1 工況描述

本橋位于中國貴州省福泉市三江口鎮(zhèn)，修建于1986年，主拱圈結(jié)構(gòu)為雙曲拱橋，橋面寬度0.82+7.0+0.82m，主跨徑55.0m，矢高6.85m，主拱圈由4根拱肋組成，拱肋截面高度70cm、寬度60cm，兩端分別有簡支梁混凝土引橋，設(shè)計時速30km/h。本橋在設(shè)計時主拱圈采用懸鏈線作為拱軸線，拱軸系數(shù)m=1.988。

2 車一橋耦合系統(tǒng)（VBCV）

目前，數(shù)值分析方法因為其較低的限制性和較低的花費被大量用于分析VBCV問題[6]，而且有文獻驗證了數(shù)值分析方法與試驗結(jié)果能較好地吻合[7-8]，對于車橋組合結(jié)構(gòu)，通常采用Newmark和Runge—Kutta法來正確的解決耦合振動[9]。

本文運用疊加方法解決VBCV問題，首先分別建立車輛和橋梁的方程，并通過邊界條件將兩個子系統(tǒng)進行組合疊加，在不同情況進行分析得到相應(yīng)的動態(tài)響應(yīng)。

2.1 車輛模型

采用經(jīng)典三軸車輛模型，該車總質(zhì)量30t，將車輛劃分成一個車體質(zhì)量，6個車輪及懸掛質(zhì)量。用彈簧一阻尼系統(tǒng)模擬車體一車輪和車輪一橋面之間的聯(lián)系，從而將車輛模型簡化成質(zhì)量、彈簧、阻尼系統(tǒng)所組成的三維有限元模型體系[10]。其中M為車體質(zhì)量;m1～m6為車輪及懸掛質(zhì)量;kt1～kt6ct1～ct6為懸掛的剛度和阻尼;ks1～ks6、cs1一cs6為輪胎的剛度和阻尼。同時還考慮了車輪豎向位移z1～z6、車體豎向位移zv、車體橫向及縱向轉(zhuǎn)動位移θv、φv，共9個自由度。

2.2 橋梁模型

對于本橋梁拱肋部分采用虛擬梁單元和實體單元拼接而成，左右兩側(cè)立柱采用梁單元，中間部分用梁單元模擬填料，橋面板用板單元模擬，根據(jù)施工圖建立三維有限元模型（見圖1），其中行車方向為x方向，圖示方向右為正;橫向為y方向，圖示方向左上方為正;豎向為z方向圖示方向上為正，圖1左下方符號處為原點。

2.3車一橋組合系統(tǒng)

1）車輛模擬

對于三軸車輛，建立起三軸車輛模型[11]并根據(jù)D' Alenmber原理，建立起動力學(xué)方程

由于在空腹時拱橋中，由于集中力的存在，使得拱的恒載壓力線為一條在集中力作用點處有轉(zhuǎn)折的曲線，而不是懸鏈線，甚至不是一條光滑的曲線[14]。所以在拱肋內(nèi)會產(chǎn)生彎矩和剪力。

通過查詢并計算在考慮彈性壓縮影響下拱肋任意截面的內(nèi)力，則任意截面的偏離彎矩為

3 試驗及數(shù)值分析

由于試驗條件的限制和花費，很難通過監(jiān)測得到完整的試驗數(shù)據(jù)，故部分結(jié)果需要通過數(shù)值分析的方法得到。

3.1測點布置

根據(jù)本文的目的，在試驗中布置測點對橋梁立柱受到的應(yīng)力及位移進行監(jiān)測，測點布置于立柱與拱肋相交處，如圖2和圖3所示。

3.2 試驗數(shù)據(jù)

在圖1中，選用每排的立柱及測點兩側(cè)實體單元作為研究對象，數(shù)值模擬中單元號如表1所示。監(jiān)測所得位移值及內(nèi)力值如表2和表3所示（均用數(shù)值模擬中單元號代替）。

3.3數(shù)值分析結(jié)果

對橋梁先進行整體分析，得到立柱的內(nèi)力及位移值如表4和表5所示。對比表2-表5數(shù)值可以發(fā)現(xiàn)，內(nèi)力及位移的數(shù)值模擬的結(jié)果相較于現(xiàn)場監(jiān)測相差并不是太大，故本文將將數(shù)值模擬所得到的結(jié)果作為試驗值進行分析。

1）動力特性

根據(jù)橋梁的材料、剛度、幾何尺寸等條件得到橋梁的前四階振型如圖4所示。

2）車輛勻速前進時立柱的動態(tài)響應(yīng)

由表3和表4可知，同一排3根立柱各值差異不大，現(xiàn)選用中間立柱作為研究對象。為了對比六立柱在汽車不同行駛速度下的動力響應(yīng)，設(shè)計時速分別為20km/h和30km/h的車輛勻速通過橋，測量6根立柱在兩種情況下的軸力、剪力及變形，由表格形式列出。

通過對比表6、表7可以發(fā)現(xiàn)，在未超過設(shè)計時速時，隨著車速的增大，軸力及z方向剪力最大值有少量增加，由于車輛對稱作用，y方向剪力幾乎等于零。在變形方面，在x方向有少量變化，y方向和z方向幾乎沒有變形，故在為超速時，速度對變形影響不大。

3.4 超速下車一橋耦合現(xiàn)象

由于現(xiàn)在車輛大多都是超速行駛，故本小節(jié)研究在超速情況下立柱的響應(yīng)。

設(shè)置汽車分別以40km/h、50km/h的速度勻速通過橋梁，對其進行模擬，可以得到六根立柱在超速下的內(nèi)力及變形，如圖5、圖6所示。

通過對比圖5、圖6發(fā)現(xiàn)，速度的增加對立柱所承受的軸應(yīng)力有著明顯的影響，軸應(yīng)力最大值及最小值均呈增長趨勢且變化幅度較大，20km/h時軸應(yīng)力最大值為31. 87kN/m2、最小值為-28. 8kN/m2;而50km/h時軸應(yīng)力最大值已達到59. 40kN/m2，最小值增長為-61. 97kN/m2。故在長期超速荷載的作用下，立柱的軸向性能會受到一定的損壞。

y方向剪應(yīng)力由于汽車荷載對稱作用變化并不明顯。而z方向剪應(yīng)力則呈現(xiàn)出先上升后下降的趨勢，由圖5（e）～（f）曲線變化可預(yù)測在車速達到35 km/h左右時z方向剪應(yīng)力可能會出現(xiàn)極值，此時立柱將承受較大的剪切力，可能會造成沿橋向剪切破壞。

4結(jié)論

（1）通過對比現(xiàn)場監(jiān)測結(jié)果及數(shù)值模擬結(jié)果可以發(fā)現(xiàn)二者值大致接近，故本模擬結(jié)果可以使用。

（2）在未超過設(shè)計時速時，隨著車速的增加，立柱軸力及z方向剪力有少量增加，y方向剪力幾乎不變，x方向有較小變形。

（3）在超速時，車輛速度的增長對立柱軸應(yīng)力有著顯著的影響。立柱z方向剪應(yīng)力在車輛超過設(shè)計時速后不久將達到峰值，故在長期車輛超速行駛下立柱有可能產(chǎn)生沿橋向剪切破壞。

（4）車速的變化對橋梁變形的影響并不是太大，故超速現(xiàn)象并不會直接影響到橋梁的變形損壞。

參考文獻：

[1] ATAEI S，MIRIS. Investigating dynamic amplifificationfactor of railway masonry arch bridges through dynamicload tests[J].Construction and Building Materials，2018 ，183：693 -705.

[2] AKBARI M R，NIMAFAR M，CANJI D D，et，a1.Investi-gation on non-linear vibration in arched beam for bridgesconstruction via AGM method[J].Applied Mathematicsand Computation， 2017， 298： 95 -110.

[3]MOREIRA V N，F(xiàn)ERNANDES J，MATOS J C，et，al.Reliability- based assessment of existing masonry archrailway bridges[J].Construction and Building Materi-als ，2016 ，115：544-554.

[4]CONEN H.DOCAN M，KARACASU M，et..Structuralfailures in refrofifit historical murat masonry arch rmdge[J]. Engineering Failure Analysis，2013， 35： 334-342.

[5]AYDIN A C，OZKAYA S G.The fifinite element analy-sis of collapse loads of single-spanned historic masonryarch bridges（ Ordu， Sarpdere Bridge）[J].EngineeringFailure Analysis， 2018， 84： 131 -138.

[6]CHARLES W J，DI S，YOZO F.Prediction of vehicle-induced local responses and application to a skewedgirder bridges[J].Engineering Structure， 2011 ，33：1088-1 097.

[7]ASHEBO D C，CHAN T H T，YU L Evaluation of dy-namic loads on a skew box girder continuous bridge partI： field test and modal analysis[J].Engineering struc-tures ，2007;29（6）：1 052-1 063.

[8] WANG T L，HUANG D， SHAHAWY M. Dynamic be-havior of slant-legged rigid-frame highway bridge[J].JStructure Engineering， 1994， 120（3）：885-902.

[9] KIM C W，KAWATANI M， KIM K B.Three-dimension-al dynamic analysis of vehicleinterection with roadway roughness[J].Computation structures， 2005， 83（ 19 -20）：1 627-1 645.

[10]鄧露，何維，俞揚，等，公路車一橋耦合振動的理論研究和應(yīng)用研究進展[J].中國公路學(xué)報，2018（ 31）：38 -54.

[11]《中國公路學(xué)報》編輯部，中國橋梁工程學(xué)術(shù)研究綜述. 2014[J].中國公路學(xué)報，2014，27（5）：1-96.

[12]HUMAR J L，KASHIF A H.Dynamic response analy-sis of slab type bridges[J].Joumal of Structures Engi-neering， 1995， 121（1）：48-62.

[13]LEITO P N， SILVA J G S，VELLASCO P C G D， et，a1. Composite（ steel-concrete） highway bridge fatigueassessment[ J]. Joumal of Constructional Steel Re-search， 2011， 67（1）：14-24.

[14]汪蓮，任偉新，趙青，等，橋梁工程[M].合肥：合肥工業(yè)大學(xué)出版社，2011：248-260.

[15]顧安邦，孫國柱，公路橋涵設(shè)計手冊（拱橋下冊）[S].北京：人民交通出版社，1996.

（責(zé)任編輯：丁寒）