范佳豪, 劉慧杰, 姜會民, 劉小兵,2,3, 楊 群,2,3

(1. 石家莊鐵道大學 土木工程學院，河北 石家莊 050043；2. 河北省風工程和風能利用工程技術創(chuàng)新中心，河北 石家莊 050043；3. 石家莊鐵道大學 省部共建交通工程結(jié)構(gòu)力學行為與系統(tǒng)安全國家重點實驗室，河北 石家莊 050043)

0 引言

現(xiàn)代的橋梁不斷向著大跨和輕柔的方向發(fā)展，對風荷載愈發(fā)敏感[1-2]。扁平箱梁由于具有橫向抗彎、抗扭性能好及良好的顫振穩(wěn)定性[3-5]等優(yōu)點，被廣泛地應用于實際工程中，如國內(nèi)建設的江東大橋[6]、青島海灣大橋[7]，主梁都采用了扁平箱梁。掌握扁平箱梁的氣動力特性及其發(fā)生機理對于扁平箱梁橋的抗風設計及氣動優(yōu)化[8-9]具有重要的參考價值。

已有學者對扁平箱梁進行了一些研究。其中李永樂等[10]以瓊州海峽為背景，通過數(shù)值模擬對其主跨為3 500 m的超大跨度懸索橋中的扁平單箱梁進行了氣動力特性的研究。研究結(jié)果發(fā)現(xiàn),扁平單箱梁的三分力系數(shù)絕對值隨風攻角增大而增大，較大風攻角下的扁平單箱梁穩(wěn)定性能降低；而扁平箱梁的斯托羅哈數(shù)較小，表明該主梁型式具有較為良好的渦振性能。祝志文等[11]以丹麥大帶東橋主跨加勁梁為研究對象，對不同風攻角下的扁平箱梁進行了二維數(shù)值計算，并將二維和三維主梁的整體氣動力特性、表面壓力分布進行對比。結(jié)果發(fā)現(xiàn)，在典型風攻角下扁平箱梁具有單一斯托羅哈數(shù)，且均大于成橋的狀態(tài)。楊陽等[12]以寸灘長江大橋主橋的加勁梁為研究對象，采用風洞試驗和數(shù)值模擬相結(jié)合的方法對扁平鋼箱梁的氣動力特性和渦振性能進行了研究。結(jié)果表明，該寬體扁平鋼箱梁在不同風速下的靜力三分力系數(shù)十分接近，受雷諾數(shù)的影響較小，該梁對風攻角變化比較敏感。李春光等[13]通過節(jié)段模型試驗的方法，研究了寬高比為11.08的扁平鋼箱梁在不同條件下對其氣動力特性的影響。研究結(jié)果發(fā)現(xiàn)，隨著風攻角由負變正，扁平鋼箱梁斷面的扭轉(zhuǎn)氣動力特性偏于不安全，而且在懸挑臂導流板以及底板設豎直中央穩(wěn)定板不適于提高扁平鋼箱梁斷面的氣動力特性。

從國內(nèi)外學者的研究中可以發(fā)現(xiàn)，當前的研究多以某一工程實例為背景展開，其中研究的扁平箱梁的寬高比多為某一固定值。然而，隨著交通量的日益增長，橋面逐漸向著更寬幅的方向發(fā)展，當前研究尚不能預測更寬幅扁平箱梁的氣動力特性。因此，需要對扁平箱梁進行更為系統(tǒng)的研究，揭示截面寬高比對扁平箱梁氣動力特性的影響規(guī)律，以更好地指導實際工程。以某跨海大橋初步設計方案為工程背景，通過CFD數(shù)值模擬的方法研究了4種不同寬高比扁平箱梁在7個風攻角下的氣動力特性，并對其發(fā)生機理進行了分析。

1 數(shù)值計算方法與計算模型

1.1 數(shù)值計算方法

在流體力學中通過動量方程、連續(xù)性方程和能量方程來體現(xiàn)3個最基本的守恒定律(動量守恒定律、質(zhì)量守恒定律及能量守恒定律)。其中不可壓縮流動的連續(xù)性方程、動量守恒方程及能量方程分別為

(1)

(2)

(3)

式中，ρ為靜壓；τij為應力張量；gi、Fi分別為i方向上的重力體積力和外部體積力。

數(shù)值計算方法采用雷諾平均法。CFD數(shù)值分析中，合理選取湍流模型是準確模擬箱梁斷面周邊流場的基本前提[14]。SSTk-ω湍流模型綜合了標準k-ε模型和k-ω模型各自的優(yōu)缺點，使得SSTk-ω湍流模型能獲得優(yōu)于標準k-ε模型和k-ω模型的模擬結(jié)果。通常認為SSTk-ω湍流模型是RANS模型中求解表現(xiàn)較好的湍流模型之一，比較適合于流動分離區(qū)的求解[15]。所以數(shù)值計算模型采用SSTk-ω湍流模型。

1.2 數(shù)值計算模型

以某跨海大橋初步設計方案為工程背景，扁平箱梁節(jié)段模型尺寸如圖1所示，扁平箱梁模型寬度B=370 mm，高度H=64 mm，風攻角選用-6°～6°，風攻角方向沿風向斜向上為正，風速采用6 m/s。此扁平箱梁的寬高比為B/H=5.78，為了進一步研究不同寬高比對扁平箱梁氣動效應的影響，增加B/H=5.0、8.0、10.0 3個不同寬高比的扁平箱梁，然后對4個不同寬高比的扁平箱梁進行數(shù)值模擬計算。對于扁平箱梁的二維數(shù)值計算，首先應確定它的邊界條件和計算域，計算域選用矩形域進行計算，速度入口設置為velocity-inlet，出口設置為outflow，上、下邊界為symmetry，因?qū)ΨQ邊界距箱梁模型中心的尺寸、速度入口邊界對計算網(wǎng)格數(shù)量影響不大，故查閱文獻[16]、文獻[17]后將其確定為10B。經(jīng)過試算后可知，當尾流區(qū)長度大于20B時，計算結(jié)果不再發(fā)生變化，故尾流區(qū)長度確定為20B。圖2為扁平箱梁的計算域和邊界條件示意圖。通過試算當首層網(wǎng)格厚度為0.03 mm時可以保證無量綱Y+在1.0左右，故設置首層網(wǎng)格尺寸為0.03 mm。圖3給出了B/H=5.78扁平箱梁的局部網(wǎng)格示意圖。

圖1 節(jié)段模型尺寸圖(單位：mm)

圖2 計算域和邊界條件示意圖

圖3 扁平箱梁局部網(wǎng)格示意圖

1.3 參數(shù)定義及結(jié)果驗證

扁平箱梁的阻力系數(shù)CD、升力系數(shù)CL和扭矩系數(shù)CM的定義分別為

CD=2FD/ρU2B

(4)

CL=2FL/ρU2B

(5)

CM=2MT/ρU2B2

(6)

式中，F(xiàn)D、FL、MT分別為順風向阻力、橫風向升力和扭矩；ρ為空氣密度；U為風速；B為扁平箱梁寬度。三分力示意圖如圖4所示。

圖4 三分力示意圖

為了驗證數(shù)值計算結(jié)果的準確性，將數(shù)值計算結(jié)果與文獻[18]的風洞試驗結(jié)果(寬高比B/H=5.78)進行了對比，如圖5所示。從圖5可以看到，對于阻力系數(shù)而言，數(shù)值計算結(jié)果和風洞試驗結(jié)果略微有些差別，數(shù)值相差不大，但二者在變化規(guī)律上一致，均隨風攻角的增大呈現(xiàn)出先緩慢減小后逐漸增大的變化規(guī)律。對于扁平箱梁的升力系數(shù)和扭矩系數(shù)，數(shù)值計算結(jié)果與風洞試驗結(jié)果吻合較好。以上結(jié)果說明，數(shù)值計算結(jié)果具有一定的可靠性。

圖5 寬高比為5.78時三分力系數(shù)的對比

2 寬高比對扁平箱梁氣動力特性的影響規(guī)律

2.1 三分力系數(shù)

圖6為不同寬高比下扁平箱梁的三分力系數(shù)隨風攻角的變化曲線，從圖6可以看到扁平箱梁在不同寬高比下的阻力系數(shù)隨風攻角的變化規(guī)律基本一致，隨著風攻角由-6°增加到6°，扁平箱梁的阻力系數(shù)先緩慢減小后逐漸增大，在風攻角為6°時達到其最大值，最大值為寬高比為5.0的時候，其最大值約為0.14。同時可以看到，隨寬高比的增大，扁平箱梁在不同風攻角下的阻力系數(shù)均呈現(xiàn)出逐漸減小的規(guī)律。當風攻角為6°時，寬高比為5.0、5.78、8.0和10.0的扁平箱梁其阻力系數(shù)分別約為0.14、0.12、0.1和0.09。

圖6 不同寬高比扁平箱梁三分力系數(shù)隨風攻角的變化曲線

扁平箱梁在不同寬高比下的升力系數(shù)均隨風攻角的增加表現(xiàn)出了先增大后減小的變化規(guī)律。當B/H=5.0和5.78時，升力系數(shù)極大值發(fā)生在風攻角為2°時;當B/H=8.0和10.0時，升力系數(shù)極大值發(fā)生在風攻角為4°時。當風攻角為-6°～2°時，扁平箱梁在不同寬高比下的升力系數(shù)基本一致，隨著風攻角的增加，由-0.7左右逐漸增加至0.2附近。當風攻角為4°～6°時，扁平箱梁的升力系數(shù)隨寬高比的增大表現(xiàn)出了逐漸增大的規(guī)律。當B/H=5.0和5.78時，扭矩系數(shù)在風攻角達到0°時趨于平穩(wěn)，當B/H=8.0和10.0時，扭矩系數(shù)在風攻角達到2°時趨于平穩(wěn)。對比變扁平箱梁在不同寬高比下的扭矩系數(shù)可以看到，當風攻角為-6°～0°時，隨著寬高比的增加扭矩系數(shù)呈現(xiàn)出逐步減小的變化規(guī)律，當風攻角為4°～6°時，隨著寬高比的增加扭矩系數(shù)略有增大。當風攻角為-6°時，扭矩系數(shù)絕對值最大，寬高比為5.0、5.78、8.0和10.0的扁平箱梁其扭矩系數(shù)最大絕對值分別約為0.09、0.11、0.13和0.14。綜合以上結(jié)果不難發(fā)現(xiàn)，寬高比的增加會使扁平箱梁受到的阻力減小，但會使其受到更大的升力和扭矩。

2.2 風壓系數(shù)

扁平箱梁氣動力特性的發(fā)生機理可通過風壓分布進行分析[19-20]。為此，給出了不同寬高比下扁平箱梁表面的風壓分布，如圖7所示，需要說明的是，受篇幅限制，僅展示風攻角為±6°和0°。圖7的風壓系數(shù)定義為式(7)，其中P-P0為模型表面壓力與參考壓力的壓力差。為便于描述，將扁平箱梁的角點依次記作a、b、c、d、e、f，如圖4所示。記ab面、bc面及cd面分別為迎風側(cè)直腹板、上頂板及背風側(cè)直腹板，de面、ef面及fa面分別為背風側(cè)斜腹板、下底板及迎風側(cè)斜腹板。

圖7 不同寬高比扁平箱梁表面的平均風壓系數(shù)

Cp=(P-P0)/0.5ρU2

(7)

圖7為±6°和0°風攻角下不同寬高比下扁平箱梁表面平均風壓系數(shù)隨寬高比的變化曲線，從圖7可以看到扁平箱梁在不同風攻角下的平均風壓系數(shù)隨寬高比的變化規(guī)律基本一致。在-6°風攻角下，不同寬高比扁平箱梁受正壓的部分為ab面和部分bc面，其余部分均受負壓的影響。不同寬高比扁平箱梁ab面的平均風壓系數(shù)由角點a到角點b逐漸增大。扁平箱梁cd面和de面的平均風壓系數(shù)變化幅度比ef面和fa面要小。在0°風攻角下，不同寬高比下扁平箱梁的平均風壓系數(shù)僅有ab面為正值，其他面均為負值。在bc面上不同寬高比的負壓由角點b到角點c表現(xiàn)出逐漸減小的規(guī)律。還可以明顯看出在角點f附近受到的負壓最大，其值約為-0.9。在6°風攻角下，不同寬高比扁平箱梁受正壓的部分為ab面和靠近角點a的fa面，其余部分均受負壓的影響。不同寬高比扁平箱梁ab面的平均風壓系數(shù)由角點a到角點b逐漸減小。對于不同寬高比扁平箱梁的阻力系數(shù)，起主要影響作用的是扁平箱梁ab面、cd面、de面及fa面的風壓。ab面的正壓隨著風攻角的增大變化不大，cd面、de面及fa面的負壓隨著風攻角的增大逐漸減小，所以cd面、de面和fa面的風壓大小決定了箱梁阻力系數(shù)隨著寬高比的增大而逐步減小。對于不同寬高比扁平箱梁的升力系數(shù)，起主要影響作用的是扁平箱梁bc面和ef面的風壓，隨著風攻角的增大bc面的負壓越來越大，ef面的正壓越來越大，同時也受到扁平箱梁de面和fa面的影響。從圖7中可以看出扁平箱梁bc面和ef面的壓力差隨著寬高比的增大變化不明顯，這較好解釋了圖6(b)升力系數(shù)規(guī)律的一致性。

3 氣動力特性的流場機理分析

下面從流場角度來解釋不同寬高比扁平箱梁風壓系數(shù)的變化規(guī)律。圖8為不同寬高比-6°風攻角下的時均流線圖，從圖8可以看到，不同寬高比扁平箱梁所形成旋渦的位置基本相同，均出現(xiàn)在扁平箱梁的cd面、fa面和de面附近，但旋渦大小不一。當寬高比為5.0時，扁平箱梁cd面、fa面和de面均出現(xiàn)了旋渦，對扁平箱梁產(chǎn)生負壓作用；隨寬高比的增大，扁平箱梁cd面形成的旋渦越來越小，致使扁平箱梁所受到的負風壓逐漸減小，這驗證了圖7(a)中cd面和de面的平均風壓系數(shù)絕對值隨寬高比的增大而減小的變化規(guī)律。

圖8 -6°不同寬高比箱梁時均流線圖

圖9為不同寬高比0°風攻角下的時均流線圖，從圖9可以看到,靠近角點b附近的旋渦中心隨寬高比的增大，旋渦中心在逐漸向角點b移動。不難發(fā)現(xiàn)的是當旋渦隨寬高比增大逐漸變小的時候，cd面附近仍有3個旋渦。當寬高比為5.0時，旋渦對扁平箱梁產(chǎn)生負壓作用。隨寬高比的增大，扁平箱梁cd面形成的旋渦越來越小，致使cd面和de面的平均風壓系數(shù)絕對值隨寬高比的增大而減小。

圖9 0°不同寬高比箱梁時均流線圖

圖10為不同寬高比6°風攻角下的時均流線圖，從圖10可以看到，不同的寬高比所形成旋渦的位置基本相同。當寬高比為5.0和5.78時，扁平箱梁bc面、cd面均出現(xiàn)了旋渦。隨寬高比的增大，扁平箱梁bc面形成的旋渦無明顯變化，所以扁平箱梁bc面的風壓變化表現(xiàn)不明顯。但扁平箱梁cd面形成的旋渦隨寬高比的增大越來越小，當寬高比為8.0時，扁平箱梁cd面不再形成旋渦，這致使cd面的負壓隨寬高比的增大而減小，所以扁平箱梁cd面的平均風壓系數(shù)絕對值隨寬高比的增大而減小。

圖10 6°不同寬高比箱梁時均流線圖

4 結(jié)論

基于數(shù)值模擬研究了-6°～6°風攻角范圍內(nèi)不同寬高比扁平箱梁的氣動力特性，得到了以下結(jié)論：

(1)扁平箱梁的阻力系數(shù)隨風攻角由-6°增大到6°呈現(xiàn)出先緩慢減小后逐漸增大的變化規(guī)律；升力系數(shù)呈現(xiàn)出先增大后減小的變化規(guī)律；扭矩系數(shù)呈現(xiàn)出先增大后趨于平穩(wěn)的變化規(guī)律。寬高比的增加會使扁平箱梁受到的阻力減小，但會使其受到更大的升力和扭矩。

(2)不同的寬高比扁平箱梁所受的風壓不同，扁平箱梁迎風側(cè)直腹板所受風壓基本不隨寬高比發(fā)生變化，迎風側(cè)斜腹板所受風壓隨寬高比的變化規(guī)律與風攻角有關，其他面所受風壓隨寬高比的增大呈現(xiàn)出逐漸增大的規(guī)律。不同的風攻角下扁平箱梁所受的風壓也不同，隨著風攻角的增大，上頂板所受負壓越來越強，迎風側(cè)斜腹板所受負壓越來越弱。

(3)旋渦主要存在于扁平箱梁上頂板和尾流處，負攻角下在箱梁的迎風側(cè)斜腹板處也有旋渦出現(xiàn)。不同寬高比下扁平箱梁所形成旋渦的位置基本相同，但大小和強度不同。箱梁上頂板所形成的旋渦對箱梁升力系數(shù)的大小起決定性作用。