胡立華，李德建

(中南大學(xué)土木工程學(xué)院，湖南長(zhǎng)沙410075)

風(fēng)對(duì)結(jié)構(gòu)的振動(dòng)反應(yīng)分為發(fā)散振動(dòng)和限幅振動(dòng)，顫振和馳振屬于發(fā)散振動(dòng)，渦激振和抖振屬于限幅振動(dòng)。發(fā)散振動(dòng)具有危險(xiǎn)性，是一種急性病，大跨度橋梁在施工中結(jié)構(gòu)剛度小，容易在短時(shí)間內(nèi)受風(fēng)致振動(dòng)影響，因此必須對(duì)施工階段顫振穩(wěn)定性進(jìn)行研究，大跨度橋梁的抗風(fēng)性能首先取決于其動(dòng)力特性［1］。因此，必須根據(jù)懸臂施工的剛度對(duì)結(jié)構(gòu)的動(dòng)力性能作出評(píng)價(jià)和分析，以保證施工的順利進(jìn)行。以酉水大橋(80+145+80)m斜交高墩大跨預(yù)應(yīng)力混凝土連續(xù)箱梁橋?yàn)楣こ瘫尘?，?duì)最高橋墩施工階段、最大懸臂施工階段、成橋階段3個(gè)階段的動(dòng)力特性進(jìn)行分析，依此分析施工階段的顫振穩(wěn)定性，得出一些結(jié)論。

1 結(jié)構(gòu)動(dòng)力特性分析方法

1.1 結(jié)構(gòu)動(dòng)力特性有限元分析

利用Midas－civil對(duì)橋梁的動(dòng)力特性進(jìn)行分析。結(jié)構(gòu)動(dòng)力特性主要與結(jié)構(gòu)的質(zhì)量分布和剛度特性(即模型的邊界條件、材料的彈性模量和截面的幾何特性)有關(guān)，故必須精確模擬分析橋梁結(jié)構(gòu)的質(zhì)量和剛度特性。在成橋階段，橋梁鋪裝等二期恒載也必須轉(zhuǎn)化為均布質(zhì)量分布于主梁?jiǎn)卧?。酉水大橋?dòng)力有限元模型，在軟件Midas－civil中采用空間梁?jiǎn)卧M，每個(gè)節(jié)點(diǎn)6個(gè)自由度，對(duì)應(yīng)3個(gè)線自由度和3個(gè)角自由度，同時(shí)考慮梁的剪切變形。

結(jié)構(gòu)的動(dòng)力方程為:

式中:［M］，［C］和［K］分別為結(jié)構(gòu)的質(zhì)量矩陣、阻尼矩陣和剛度矩陣;{ˉδ}，{ˉδ}和{δ}分別為結(jié)構(gòu)的加速度、速度和位移向量;{F}為作用在橋梁空間梁?jiǎn)卧系牧ο蛄俊?/p>

求結(jié)構(gòu)自振特性時(shí)，一般不考慮阻尼的影響。令［C］=0，{F}=0，則得到其無阻尼自振方程為:

式(2)具有非零解，則:

即式(3)的特征方程為:

對(duì)于式(4)特征值的求解方法較多，常有的有多重Ritz向量法、Lanczos向量迭代法、子空間迭代法、Rayleigh－Ritz法和逆迭代法。從結(jié)構(gòu)分析的角度來說，并不是對(duì)所有振型和頻率都考慮，一般只需求出較低的幾階即可，這樣不僅節(jié)省振型存儲(chǔ)空間，而且縮短計(jì)算時(shí)間。本文特征值的求解采用子空間迭代法進(jìn)行，這種方法能充分利用［M］和［K］2個(gè)矩陣的稀疏帶狀性質(zhì)，從而提高求解效率。經(jīng)驗(yàn)表明，子空間迭代法是目前求解大型結(jié)構(gòu)自振頻率和振型的最有效方法之一［2－3］。

本文使用子空間迭代法計(jì)算橋梁結(jié)構(gòu)前10階頻率和振型，并給出前5階頻率和振型［4］。分析控制最大迭代次數(shù)設(shè)為20，迭代收斂誤差設(shè)定為le－10。

1.2 結(jié)構(gòu)動(dòng)力特性現(xiàn)場(chǎng)測(cè)試

為對(duì)比結(jié)構(gòu)動(dòng)力特性有限元分析方法的正確性，同時(shí)對(duì)結(jié)構(gòu)的動(dòng)力特性進(jìn)行現(xiàn)場(chǎng)測(cè)試。試驗(yàn)采用由拾振器(位移、速度和加速度傳感器)、數(shù)據(jù)采集儀、電荷放大器及動(dòng)力分析儀組成的動(dòng)力測(cè)試系統(tǒng)［5］，本結(jié)構(gòu)動(dòng)力性能現(xiàn)場(chǎng)試驗(yàn)采用德國(guó)生產(chǎn)的國(guó)際公認(rèn)的 HBM動(dòng)力測(cè)試與分析系統(tǒng)完成。HBM動(dòng)力測(cè)試與分析系統(tǒng)如圖1所示。

圖1 振動(dòng)測(cè)試、分析示意圖Fig.1 Vibration test and analysis diagram

當(dāng)橋跨結(jié)構(gòu)無車輛通過時(shí)，記錄結(jié)構(gòu)在環(huán)境隨機(jī)激勵(lì)下的加速度響應(yīng)，通過頻譜分析可獲得結(jié)構(gòu)的自振特性。

2 酉水大橋斜交高墩大跨度連續(xù)箱梁橋動(dòng)力特性分析

2.1 工程概況

酉水大橋主橋?yàn)橐蛔苯桓叨沾罂缍葢冶蹪仓A(yù)應(yīng)力混凝土連續(xù)箱梁橋，跨徑布置為80 m+ 145 m+80 m。主墩為5號(hào)和6號(hào)墩，墩身采用雙肢矩形空心墩，順?biāo)鞑贾?，主橋墩軸線與橋梁中心線成65°夾角。其中，5號(hào)墩高61 m，在墩底往上20 m處設(shè)置一個(gè)變截面段，橫向放坡40∶1，縱向放坡80∶1，20 m以上部分為等截面，空心墩壁厚125 cm。橋梁布置如圖2所示。

圖2 斜交高墩構(gòu)造圖Fig.2 High skew pier structure

2.2 各施工階段結(jié)構(gòu)自振測(cè)試

為進(jìn)行斜交高墩大跨度連續(xù)箱梁橋在最高墩施工階段、最大懸臂施工階段及成橋階段的動(dòng)力特性研究，在橋梁邊跨跨中以及墩身橫向、縱向設(shè)置加速度計(jì)測(cè)試橋梁豎向、橫向、縱向在動(dòng)力激振作用下的加速度值，加速度布置及編號(hào)如圖3和圖4所示。

圖3 邊跨跨中截面加速度傳感器布置圖Fig.3 Side span midspan section acceleration sensor layout

圖4 橋梁5號(hào)橋墩截面加速度傳感器布置圖ig.4 Section of the bridge No5 pier acceleration sensor layout

通過HBM動(dòng)態(tài)數(shù)據(jù)測(cè)試儀采集各施工階段結(jié)構(gòu)各方向脈動(dòng)加速度波，應(yīng)用HBM動(dòng)態(tài)數(shù)據(jù)分析儀分析各結(jié)構(gòu)加速度頻譜，從而得到各測(cè)試點(diǎn)方向的結(jié)構(gòu)自振頻率。

2.3 最高墩施工階段結(jié)構(gòu)動(dòng)力特性

通過Midas－Civil有限元分析軟件建立5號(hào)墩最高墩模型，將自重荷載轉(zhuǎn)化為質(zhì)量，計(jì)算最高墩階段動(dòng)力特性，給出前5階計(jì)算結(jié)果如表1和圖5所示。

表1 最高墩動(dòng)力特性計(jì)算結(jié)果Table 1 Dynamic characteristic calculation results of themaximum pier

從表1可知，在最高墩施工階段，橋墩動(dòng)力特性主要表現(xiàn)為橋墩縱漂。

通過HBM動(dòng)態(tài)數(shù)據(jù)采集儀，在最高橋墩階段，測(cè)點(diǎn)2和3的實(shí)測(cè)加速度脈動(dòng)波型如圖5所示。

圖5 最高橋墩各測(cè)點(diǎn)加速度頻譜Fig.5 Maximum pier acceleration spectrum of each measuring point

最高墩階段結(jié)構(gòu)振動(dòng)頻率測(cè)試結(jié)果與理論計(jì)算值對(duì)比如表2所示。

表2 最高墩階段振動(dòng)測(cè)試頻率與理論頻率比較Table 2 Vibration test frequency and theoretical frequency comparison of themaximum pier

實(shí)測(cè)橋墩橫向、縱向一階自振頻率均比理論計(jì)算值大，說明最高墩橫向、豎向剛度良好。

2.4 最大懸臂階段動(dòng)力特性分析

同理，通過Midas－Civil有限元分析軟件建立5號(hào)墩最大懸臂模型，將自重荷載轉(zhuǎn)化為質(zhì)量，計(jì)算最大懸臂階段動(dòng)力特性，給出前5階計(jì)算結(jié)果如表3和圖6所示。

表3 最大懸臂階段動(dòng)力特性計(jì)算結(jié)果Table 3 Dynamic characteristic calculation results of the largest cantilever

從表3可知，在最大懸臂施工階段，橋梁T構(gòu)的動(dòng)力特性主要表現(xiàn)為橋墩縱漂。

通過HBM動(dòng)態(tài)數(shù)據(jù)采集儀，在最大懸臂施工階段，測(cè)點(diǎn)1，2和3的實(shí)測(cè)加速度脈動(dòng)波型如圖6所示。

圖6 最大懸臂階段各測(cè)點(diǎn)加速度頻譜Fig.6 The largest cantilever acceleration spectrum of each measuring point

最大懸臂階段結(jié)構(gòu)振動(dòng)頻率測(cè)試結(jié)果與理論計(jì)算值對(duì)比見表4。

表4 最大懸臂階段振動(dòng)測(cè)試頻率與理論頻率比較Table 4 Vibration test frequency and theoretical frequency comparison of the largest cantilever

從表4可知，最大懸臂結(jié)構(gòu)3測(cè)點(diǎn)實(shí)測(cè)結(jié)構(gòu)頻率比理論計(jì)算值大，說明最大懸臂階段結(jié)構(gòu)橫向、縱向豎向剛度良好。

2.5 成橋階段結(jié)構(gòu)動(dòng)力特性分析

同理，通過Midas－Civil有限元分析軟件建立酉水大橋成橋模型，主橋左邊引橋連接，右邊與橋臺(tái)連接，將二期恒載轉(zhuǎn)化為質(zhì)量，計(jì)算酉水大橋成橋階段動(dòng)力特性，給出前5階計(jì)算結(jié)果如表5和圖7所示。

表5 成橋階段動(dòng)力特性計(jì)算結(jié)果Table 5 Dynamic characteristic calculation results of the bridge

由表5可知，在成橋階段，橋梁結(jié)構(gòu)動(dòng)力特性主要表現(xiàn)為橋墩縱漂。

同理，在成橋階段，通過HBM動(dòng)態(tài)數(shù)據(jù)采集儀，在最大懸臂施工階段，測(cè)點(diǎn)1，2和3的實(shí)測(cè)加速度脈動(dòng)波型如圖7所示。

成橋階段結(jié)構(gòu)振動(dòng)頻率測(cè)試結(jié)果與理論計(jì)算值對(duì)比見表6。

表6 成橋階段振動(dòng)測(cè)試頻率與理論頻率比較Table 6 Vibration test frequency and theoretical frequency comparison of the bridge

從表6可知，主橋成橋階段實(shí)測(cè)自振頻率比理論計(jì)算值大，說明成橋階段結(jié)構(gòu)縱橋向、橫橋向、豎向剛度良好。

圖7 成橋階段各測(cè)點(diǎn)加速度頻譜ig.7 The bridge acceleration spectrum of eachmeasuring point

3 酉水大橋斜交高墩大跨度連續(xù)箱梁橋顫振研究

超過一定的風(fēng)速時(shí)，橋梁結(jié)構(gòu)吸收的能力將大于阻尼消耗的能量，該風(fēng)速成為臨界風(fēng)速Vcr，由于空氣動(dòng)力產(chǎn)生彎曲和扭轉(zhuǎn)的耦合振動(dòng)，振幅將無限增加至橋梁破壞［6］。對(duì)于顫振自激型發(fā)散振動(dòng)，應(yīng)該使其臨界風(fēng)速提高到高于橋梁在設(shè)計(jì)使用期內(nèi)橋位處可能出現(xiàn)的最大風(fēng)速［7］。

3.1 顫振檢驗(yàn)風(fēng)速

酉水大橋地處V型河谷地帶，橋面距離地面的高度為87.16 m，根據(jù)《公路橋梁抗風(fēng)設(shè)計(jì)規(guī)范》［8］，吉首地區(qū)100 a重現(xiàn)期的10 min平均年最大風(fēng)速為24.2 m/s，故得到酉水大橋主梁100 a重現(xiàn)期的設(shè)計(jì)風(fēng)速為:

考慮施工階段設(shè)計(jì)風(fēng)速的重現(xiàn)期為10 a，則施工階段的設(shè)計(jì)風(fēng)速為:

顫振檢驗(yàn)風(fēng)速為:

3.2 橋梁結(jié)構(gòu)顫振臨界分析

根據(jù)《公路橋梁抗風(fēng)設(shè)計(jì)規(guī)范》可知，對(duì)于主跨跨徑小于300 m的橋梁，如果其主梁斷面的寬高比B/H＜4，則其顫振臨界風(fēng)速的計(jì)算公式為:

或者

式中:B為主梁斷面寬;H為主梁斷面高度;ft為橋梁結(jié)構(gòu)扭轉(zhuǎn)頻率。

顫振臨界風(fēng)速取以上兩式計(jì)算結(jié)果的較小者。

當(dāng)橋面較寬，扭轉(zhuǎn)基頻較高，大跨度橋梁的抗風(fēng)穩(wěn)定性就越好;反之，當(dāng)橋面較窄，扭轉(zhuǎn)基頻較低，大跨度橋梁的抗風(fēng)穩(wěn)定性就越差［9－11］。本橋梁結(jié)構(gòu)橋面寬度為12 m不變，因此選擇橋梁結(jié)構(gòu)剛度最小，扭轉(zhuǎn)頻率最低時(shí)作顫振穩(wěn)定分析。

從上部分橋梁各施工階段結(jié)構(gòu)動(dòng)力特性分析可知，橋梁結(jié)構(gòu)在最大懸臂施工階段結(jié)構(gòu)剛度小，扭轉(zhuǎn)頻率最低，在風(fēng)荷載作用下容易發(fā)生顫振。

根據(jù)最大懸臂施工階段結(jié)構(gòu)動(dòng)力特性分析結(jié)果，計(jì)算酉水大橋在最大懸臂階段扭轉(zhuǎn)基頻ft= 1.131 Hz，取B=12 m，H=4 m。因此，臨界顫振風(fēng)速為:

式(8)中:Vcr=67.86 m/s

式(9)中:Vcr=54.288 m/s

故取Vcr=54.288 m/s

由以上的計(jì)算結(jié)果可知，Vcr＞［Vcr］=49.945 m/s，顫振臨界風(fēng)速大于顫振檢驗(yàn)風(fēng)速，因此橋梁結(jié)構(gòu)在自然風(fēng)速下不會(huì)發(fā)生顫振失穩(wěn)現(xiàn)象。

4 結(jié)論

(1)橋梁在最高橋墩施工階段、最大懸臂施工階段和成橋階段3個(gè)階段，橋梁動(dòng)力特性表現(xiàn)為主橋墩縱漂，在最大懸臂施工階段，懸臂長(zhǎng)度大，結(jié)構(gòu)剛度小，橋墩縱漂明顯。

(2)橋梁在最高墩施工階段、最大懸臂階段和成橋階段3個(gè)階段，動(dòng)力特性實(shí)測(cè)自振頻率比理論值大，說明橋梁在縱橋向、橫橋向、豎向剛度良好。

(3)橋梁結(jié)構(gòu)在各施工階段扭轉(zhuǎn)剛度不大，在最大懸臂施工階段扭轉(zhuǎn)剛度最小，扭轉(zhuǎn)頻率最低，為自然風(fēng)顫振最不利階段，其顫振臨界風(fēng)速為54.288 m/s，大于顫振檢驗(yàn)風(fēng)速49.945 m/s，橋梁結(jié)構(gòu)顫振穩(wěn)定。

［1］顏娟，張哲，黃才良.廣東金馬大橋動(dòng)力分析［J］.公路交通科技，2003，20(1):71－73.YAN Juan，ZHANG Zhe，HUANG Cailiang.Dynamic analysis for guangdong jinma bridge［J］.Journal of Highway and Transportation Research and Development，2003，20(1):71－73.

［2］楊永賢，王福敏，周筠莉.城市輕軌連續(xù)剛構(gòu)橋動(dòng)力特性分析［J］.公路交通技術(shù)，2007，(2):94－97.YANG Yongxian，WANG Fumin，ZHOU Yunli.Analysis of dynamic behavior of continuous rigid frame bridge for urban light rail［J］.Technology of Highway and Transport，2007，(2):94－97.

［3］豐碩，項(xiàng)貽強(qiáng)，汪勁豐.大跨徑連續(xù)剛構(gòu)橋的動(dòng)力性能及地震響應(yīng)分析［J］.中南公路工程，2002，30(4):77－81.FENG Shuo，XIANG Yiqiang，WANG Jinfeng.Analysis on dynamic characteristics and seismic responses for long－span continuous rigid－framed bridges［J］.Central South Highway Engineering，2002，30(4):77－81.

［4］翟建平.肇慶西江特大橋動(dòng)力特性分析［J］.鐵道科學(xué)與工程學(xué)報(bào)，2011，8(3):13－17.ZHAIJianping.Analysis on the dynamic properties of zhaoqing xijiang bridge［J］.Journal of Railway Science and Engineering，2011，8(3):13－17.

［5］王解軍，王明明，梁錦鋒.變剛度預(yù)應(yīng)力簡(jiǎn)支梁橋固有振動(dòng)分析［J］.中南公路工程，2005，30(3):62－64.WANG Jiejun，WANGMingming，LIANG Jinfeng.Natural vibration analysis for a simply supported girder bridge with varying rigidity and prestressing force［J］.Central South Highway Engineering，2005，30(3):62－64.

［6］張迪.龍華特大橋最大懸臂狀態(tài)抗風(fēng)分析［D］.哈爾濱:東北林業(yè)大學(xué)，2007.ZHAGN Di.Analysis of wind resistant of longhua bridge in the longest cantilever state［D］.Harbin:Dongbei Forestry University，2007.

［7］孫娜.橋梁抗風(fēng)控制及設(shè)計(jì)探討［J］.交通標(biāo)準(zhǔn)化，2013 (8):142－144.SUN Na.Control and design of bridge wind－resistant［J］.Communications Standardization，2013(8):142－144.

［8］JTG/TD63－01－200，公路橋梁抗風(fēng)設(shè)計(jì)規(guī)范［M］.JTG/TD63－01－2004，Design specification for highway bridges［M］.

［9］李國(guó)豪.橋梁結(jié)構(gòu)穩(wěn)定與振動(dòng)［M］.北京:中國(guó)鐵道出版社，1992.LI Guohao.Stability and vibration of bridge structure［M］.Beijing:China Railway Press，1992.

［10］張金武，李龍安.公路橋梁抗風(fēng)設(shè)防標(biāo)準(zhǔn)的初探［J］.橋梁建設(shè)，2002(3):29－32.ZHANG Jinwu，LI Longan.Discussion and study about anti－wind design criteria for highway bridges［J］.Bridge Construction，2002(3):29－32.

［11］王愛勤，李龍安.公路大跨度橋梁抗風(fēng)性能評(píng)價(jià)方法［J］.長(zhǎng)安大學(xué)學(xué)報(bào)，2006，26(5):58－62.WANG Aiqin，LI Longan.Evaluation method of windresistance capability of large－span bridge on highway［J］.Journal of Chang＇an University，2006，26(5):58－62.