楊吉新 劉 暢 黎建華 孫亭亭

(武漢理工大學(xué)交通學(xué)院1) 武漢 430063) (中國(guó)市政工程中南設(shè)計(jì)研究總院有限公司2) 武漢 430000)(中交第二航務(wù)工程勘察設(shè)計(jì)院有限公司3) 武漢 430060)

0 引 言

斜拉橋因其美觀的造型及較好的跨越能力在現(xiàn)代橋梁工程中的應(yīng)用相當(dāng)廣泛，但是大多數(shù)斜拉橋平面采用直線布置，當(dāng)?shù)匦蔚孛矎?fù)雜，橋墩位置受到限制等情況時(shí)，直線布置無(wú)法滿足設(shè)計(jì)需求，此時(shí)曲線布置成為了解決這一矛盾的方案之一[1].目前國(guó)內(nèi)外曲線斜拉橋的設(shè)計(jì)及工程案例較為少見(jiàn)，國(guó)內(nèi)外學(xué)者對(duì)于曲線斜拉橋的理論探討也不多.世界上首座S形曲線斜拉橋修建于日本首都高速公路Katoushika—Edogawa線上，橋梁全長(zhǎng)455 m，橋型布置形式為四跨曲線連續(xù)箱梁及兩座主塔、密索結(jié)構(gòu)[2].

大量的實(shí)驗(yàn)表明:橋梁靜風(fēng)失穩(wěn)是造成橋梁傾覆結(jié)構(gòu)破壞的主要原因之一.靜風(fēng)失穩(wěn)作為橋梁風(fēng)荷載破壞中最危險(xiǎn)的破壞形式之一，對(duì)比風(fēng)載的其他作用，具有可預(yù)見(jiàn)性低、危害性大的特征[3].而具有特殊線形的斜拉橋由于結(jié)構(gòu)受力的不同，其靜風(fēng)失穩(wěn)過(guò)程也與直線斜拉橋有一定差異，為了指導(dǎo)曲線斜拉橋抗風(fēng)設(shè)計(jì)，保證其結(jié)構(gòu)穩(wěn)定性，研究曲線斜拉橋的靜風(fēng)失穩(wěn)無(wú)疑對(duì)豐富曲線斜拉橋的理論研究具有極其重要的意義[4].

1 工程背景

大連南部濱海大道西延伸線工程由主線工程和支線工程組成，主線工程起點(diǎn)為南部濱海大道星海灣跨海大橋西引橋登陸點(diǎn)，朝凌水灣項(xiàng)目規(guī)劃道路西行.支線工程以七賢東路和主線交匯處為起點(diǎn)，沿七賢東路向北，通過(guò)新建橋梁跨越黃浦路，并與學(xué)子街相連，繼續(xù)向北與紅綾路平交.本曲線斜拉橋建于南部濱海大道西延伸線工程(四標(biāo)段)，位于遼寧省大連市高新園區(qū)，緊鄰小平島與大連軟件園區(qū).主線斜拉橋平面為S形布置，采用雙塔單索面結(jié)構(gòu)，跨徑為50 m+96 m+192 m+70 m=408 m.橋型布置圖見(jiàn)圖1.

圖1 橋型布置圖

2 理論依據(jù)

風(fēng)遇到障礙會(huì)導(dǎo)致流場(chǎng)發(fā)生改變，所以風(fēng)在經(jīng)過(guò)主梁時(shí)流場(chǎng)會(huì)重分布從而產(chǎn)生風(fēng)壓對(duì)主梁施加風(fēng)荷載[5].主梁在風(fēng)載下產(chǎn)生變形，發(fā)生變形由于空間姿態(tài)的變化將對(duì)流場(chǎng)產(chǎn)生影響，又會(huì)繼續(xù)引發(fā)結(jié)構(gòu)的扭轉(zhuǎn)變形，故斜拉橋的靜風(fēng)失穩(wěn)過(guò)程為非線性變化過(guò)程.對(duì)其模擬通常采用逐段增加風(fēng)速進(jìn)行迭代求解，求得結(jié)果.具體手段為通過(guò)擬定初始風(fēng)速與風(fēng)速增量計(jì)算主梁受到的靜風(fēng)荷載，求得結(jié)構(gòu)位移、扭轉(zhuǎn)角后進(jìn)行內(nèi)部迭代，內(nèi)插求得三分力系數(shù)后疊加風(fēng)速增量重新計(jì)算靜風(fēng)荷載，進(jìn)行外部迭代，最終判斷靜風(fēng)失穩(wěn)臨界風(fēng)速[6].

斜拉橋的靜風(fēng)失穩(wěn)模型主要有二維模型和三維模型兩種.二維模型將橋梁扭轉(zhuǎn)失穩(wěn)做了大量簡(jiǎn)化，未考慮到實(shí)際情況中的彎扭耦合失穩(wěn)的組合作用，尤其曲線斜拉橋由于其線型、結(jié)構(gòu)受力的特殊性，單純將荷載和氣動(dòng)常數(shù)的計(jì)算簡(jiǎn)化為線性過(guò)程無(wú)法準(zhǔn)確模擬工程實(shí)際[7].綜上，二維模型過(guò)于簡(jiǎn)化了眾多影響因素會(huì)導(dǎo)致計(jì)算結(jié)果相對(duì)不精確甚至無(wú)法求解臨界風(fēng)速，為更精確模擬研究曲線斜拉橋的抗風(fēng)性能，采用更加符合實(shí)際簡(jiǎn)化較少的非線性三維靜風(fēng)失穩(wěn)模型[8]，平衡方程為

K(δ)Δδ=ΔP(δ)

(1)

式中：K(δ)為結(jié)構(gòu)的切線剛度矩陣，它可包含結(jié)構(gòu)幾何剛度、線彈性剛度,以及塑性剛度矩陣；Δδ,ΔP(δ)分別為結(jié)構(gòu)位移增量及風(fēng)荷載增量向量，迭代的過(guò)程中要不斷的修正δ以滿足范數(shù)的要求.

3 靜三分力系數(shù)求解

研究橋梁抗風(fēng)性能，靜三分力系數(shù)是及其重要的研究指標(biāo)，采用常規(guī)解析法求解時(shí)由于主梁橫斷面的不同，自然風(fēng)變化無(wú)常等特性不易求出準(zhǔn)確解.因此，大多數(shù)工程實(shí)際中對(duì)橋梁結(jié)構(gòu)的三分力系數(shù)進(jìn)行求解一般采用規(guī)范要求的風(fēng)洞試驗(yàn)[9].風(fēng)洞試驗(yàn)準(zhǔn)備工作及操作周期較為漫長(zhǎng)，同時(shí)試驗(yàn)費(fèi)用需求較大，故在不具備風(fēng)洞試驗(yàn)條件的情況下，采用數(shù)值模擬的計(jì)算流體力學(xué)(CFD)法求解靜三分力系數(shù)成為首選[10].為檢驗(yàn)數(shù)值模擬方法的正確性，選取大連理工大學(xué)實(shí)驗(yàn)室中某類(lèi)流線型主梁斷面風(fēng)洞試驗(yàn)作為樣本數(shù)據(jù)[11]，對(duì)比試驗(yàn)數(shù)據(jù)與數(shù)值模型計(jì)算值見(jiàn)圖2.由圖2可知，模擬值的誤差均在10%以內(nèi).通過(guò)CFD法求得的計(jì)算值與風(fēng)洞試驗(yàn)數(shù)據(jù)相差無(wú)幾，說(shuō)明數(shù)值模擬手段對(duì)計(jì)算主梁三分力系數(shù)具有較高的精度.

圖2 阻力和升力系數(shù)計(jì)算與試驗(yàn)對(duì)比圖

采用WORKBENCH建立數(shù)值計(jì)算模型并采用ICEM網(wǎng)格劃分模塊對(duì)已建立模型劃分網(wǎng)格，主梁周?chē)捎肙-block進(jìn)行網(wǎng)格劃分[12]，主梁邊界層網(wǎng)格采用拋物線式加密，最接近邊界層的網(wǎng)格最小尺寸為1×10-4m.模型計(jì)算域入口邊界距離主梁為1.5倍主梁寬，上下邊界距離主梁距離均為1.5倍主梁寬，出口邊界距離主梁3倍主梁寬.通過(guò)旋轉(zhuǎn)主梁一定角度來(lái)模擬各類(lèi)風(fēng)攻角，模型全局網(wǎng)格劃分以及局部網(wǎng)格見(jiàn)圖3.

圖3 模型網(wǎng)格劃分

在FLUENT中分別對(duì)主梁斷面在-5°,-3°,0°+3°,+5°五類(lèi)風(fēng)攻角作用的數(shù)值模型進(jìn)行計(jì)算[13]，選取風(fēng)場(chǎng)的入口邊界平均風(fēng)速為10 m/s，上下邊界采用對(duì)稱(chēng)邊界模擬，出口定義為自由流動(dòng)，計(jì)算過(guò)程中對(duì)靜三分力系數(shù)實(shí)時(shí)監(jiān)控，計(jì)算并繪制不同風(fēng)攻角作用下模型風(fēng)速以及壓力云圖見(jiàn)圖4～5.

圖4 -5°,0°,5°風(fēng)攻角模型風(fēng)速云圖

圖5 -5°,0°,5°風(fēng)攻角模型壓力云圖

由圖4～5可知，風(fēng)攻角的變化導(dǎo)致主梁周?chē)L(fēng)速、風(fēng)壓發(fā)生改變，風(fēng)攻角為±5°時(shí)，主梁周?chē)娘L(fēng)速、風(fēng)壓達(dá)到最大.隨著風(fēng)攻角的增加，主梁周?chē)鲌?chǎng)發(fā)生改變，風(fēng)速在斷面端點(diǎn)處改變方向?qū)е嘛L(fēng)速增大，風(fēng)壓重新分布，主梁受到的阻力增大.同時(shí)風(fēng)壓集中在主梁迎風(fēng)側(cè)，風(fēng)攻角為-5°時(shí)，主梁下端呈現(xiàn)負(fù)壓，風(fēng)攻角逐漸增大，下端負(fù)壓區(qū)域面積減小，上端由正壓變?yōu)樨?fù)壓，主梁受到的升力和扭矩增大.因此,風(fēng)攻角的增加會(huì)增加橋梁風(fēng)載，尤其對(duì)于曲線橋梁，自身受索力作用影響，抗扭能力較弱，風(fēng)攻角的增加會(huì)使得曲線梁橋更易出現(xiàn)失穩(wěn).

繪制主梁橫斷面靜三分力系數(shù)變化圖見(jiàn)圖6，可以看出風(fēng)攻角為-5°時(shí)Cl,Cm值即升力系數(shù)與扭矩系數(shù)均最小.隨著風(fēng)攻角從-5°～+5°增加，兩者逐漸增大.當(dāng)風(fēng)攻角為0°時(shí)Cd值即阻力系數(shù)最小，隨著風(fēng)攻角絕對(duì)值的增加，主梁迎風(fēng)面積加大，阻力系數(shù)值隨之增大.

圖6 主梁靜三分力系數(shù)

為了分析不同風(fēng)攻角情況下主梁周?chē)膲毫Ψ植甲兓?，?duì)主梁節(jié)點(diǎn)進(jìn)行編號(hào)[14]，指定主梁迎風(fēng)側(cè)下表面端點(diǎn)為起始點(diǎn)，逆時(shí)針沿主梁輪廓線依次對(duì)主梁節(jié)點(diǎn)編號(hào)，根據(jù)數(shù)值模型中提取的不同風(fēng)攻角作用下主梁各周向節(jié)點(diǎn)壓力值繪制折線圖,見(jiàn)圖7.

圖7 5類(lèi)風(fēng)攻角作用下主梁周向節(jié)點(diǎn)風(fēng)壓圖

由圖7可知，隨著風(fēng)攻角的增加，圖中兩側(cè)數(shù)據(jù)值即主梁迎風(fēng)側(cè)的節(jié)點(diǎn)風(fēng)壓變化均逐漸增大，中間數(shù)據(jù)及主梁背風(fēng)測(cè)風(fēng)壓波動(dòng)較小.迎風(fēng)側(cè)主梁上表面負(fù)壓值出現(xiàn)較大值，導(dǎo)致主梁轉(zhuǎn)矩大幅度增加.曲線斜拉橋在風(fēng)攻角為±5°時(shí)迎風(fēng)側(cè)主梁負(fù)壓出現(xiàn)極大波動(dòng)，梁體會(huì)出現(xiàn)較明顯的扭轉(zhuǎn).

4 曲線斜拉橋靜風(fēng)失穩(wěn)分析

在ANSYS中建立曲線斜拉橋有限元模型并劃分網(wǎng)格見(jiàn)圖8.

圖8 全橋有限元模型

將計(jì)算求得的靜三分力系數(shù)編輯為數(shù)據(jù)文件，采用ANSYS中程序設(shè)計(jì)語(yǔ)言編程，設(shè)置自重作用下為初始狀態(tài)，求解全橋的靜風(fēng)穩(wěn)定性[15-17].通過(guò)風(fēng)速增量逐級(jí)迭代求解非線性靜風(fēng)失穩(wěn)臨界風(fēng)速，風(fēng)速臨界點(diǎn)為迭代開(kāi)始發(fā)散時(shí).通過(guò)數(shù)據(jù)分析求得當(dāng)風(fēng)速疊加到226 m/s時(shí)程序內(nèi)部迭代開(kāi)始發(fā)散，故該曲線斜拉橋的失穩(wěn)臨界風(fēng)速值為226 m/s.選擇主跨跨中節(jié)點(diǎn)，繪制其豎向變形增量迭代曲線以及扭轉(zhuǎn)角增量迭代曲線見(jiàn)圖9.

圖9 主跨跨中處節(jié)點(diǎn)豎向變形和扭轉(zhuǎn)角

由圖9可知，橋梁結(jié)構(gòu)的豎向位移、扭轉(zhuǎn)角度呈現(xiàn)非線性變化，曲線斜拉橋跨中豎向位移、扭轉(zhuǎn)角前期變化幅度不大，當(dāng)風(fēng)速達(dá)到180 m/s時(shí)變形出現(xiàn)極速增長(zhǎng)，在靜風(fēng)失穩(wěn)臨界值226 m/s附近，位移、扭轉(zhuǎn)變形開(kāi)始發(fā)散，說(shuō)明此時(shí)結(jié)構(gòu)已經(jīng)失穩(wěn)，圖中展示的豎向位移、扭轉(zhuǎn)角的大幅度增加，表明結(jié)構(gòu)失穩(wěn)形態(tài)呈現(xiàn)為彎扭耦合變形失穩(wěn).整個(gè)失穩(wěn)過(guò)程中，風(fēng)速逐級(jí)增加風(fēng)攻角不斷變化，導(dǎo)致結(jié)構(gòu)的空間狀態(tài)不斷發(fā)生變化，每次變化后的橋梁風(fēng)攻角不止包含起始風(fēng)攻角，需加上結(jié)構(gòu)變形產(chǎn)生的扭轉(zhuǎn)角.計(jì)算結(jié)果表明，角度、空間狀態(tài)的不斷變化導(dǎo)致結(jié)構(gòu)的靜三分力系數(shù)發(fā)生變化，從而結(jié)構(gòu)所受靜風(fēng)荷載隨靜三分力系數(shù)改變而改變，主梁整體變形不遵從線性變化規(guī)律，如上圖所示呈現(xiàn)非線性變化.曲線斜拉橋由于空間結(jié)構(gòu)的復(fù)雜性，橋梁會(huì)出現(xiàn)較大的扭轉(zhuǎn)變形，對(duì)風(fēng)荷載的作用比較敏感，相比于直線斜拉橋，扭轉(zhuǎn)變形更大.

5 結(jié) 論

1) 風(fēng)攻角的變化導(dǎo)致主梁周?chē)L(fēng)速、風(fēng)壓發(fā)生改變，風(fēng)攻角為±5°時(shí)，主梁周?chē)娘L(fēng)速、風(fēng)壓達(dá)到最大.風(fēng)壓集中在主梁迎風(fēng)側(cè)，風(fēng)攻角為-5°時(shí)，主梁下端呈現(xiàn)負(fù)壓，風(fēng)攻角逐漸增大，下端負(fù)壓區(qū)域面積減小，上端由正壓變?yōu)樨?fù)壓，主梁受到的升力和扭矩增大.因此風(fēng)攻角的增加會(huì)增加橋梁風(fēng)載，尤其對(duì)于曲線橋梁，自身受索力作用影響，抗扭能力較弱，風(fēng)攻角的增加會(huì)使得曲線梁橋更易出現(xiàn)失穩(wěn).

2) 隨著風(fēng)攻角的增加，升力系數(shù)與轉(zhuǎn)矩系數(shù)呈現(xiàn)增長(zhǎng)趨勢(shì)；阻力系數(shù)在風(fēng)攻角為0°時(shí)最小，隨著風(fēng)攻角絕對(duì)值的增加，主梁迎風(fēng)面積加大，阻力系數(shù)值隨之增大.

3) 隨著風(fēng)攻角的增加，迎風(fēng)側(cè)主梁風(fēng)壓變化均逐漸增大，中間數(shù)據(jù)值以及背風(fēng)側(cè)主梁風(fēng)壓波動(dòng)較小.迎風(fēng)側(cè)主梁上表面負(fù)壓值出現(xiàn)較大值，導(dǎo)致主梁扭矩大幅度增加.曲線斜拉橋在風(fēng)攻角較大時(shí)迎風(fēng)側(cè)主梁負(fù)壓出現(xiàn)極大波動(dòng)，梁體會(huì)出現(xiàn)較明顯的扭轉(zhuǎn)現(xiàn)象.

4) 曲線斜拉橋跨中位移在豎向上前期變化幅度不大，扭轉(zhuǎn)角度變化平穩(wěn).當(dāng)風(fēng)速超過(guò)180 m/s后出現(xiàn)極速增加，在臨界值226 m/s附近結(jié)構(gòu)已經(jīng)失穩(wěn)，且失穩(wěn)形態(tài)呈現(xiàn)為彎扭耦合式.曲線斜拉橋由于空間結(jié)構(gòu)的復(fù)雜性，橋梁會(huì)出現(xiàn)較大的扭轉(zhuǎn)變形，對(duì)風(fēng)荷載的作用比較敏感，相比于直線斜拉橋，扭轉(zhuǎn)變形更大.