(1.河北工程大學(xué)水利水電學(xué)院，河北 邯鄲 056038；2.河北工程大學(xué) 河北省智慧水利重點實驗室，河北 邯鄲 056038；3.南水北調(diào)中線干線工程建設(shè)管理局河北分局，河北 石家莊 050035)

橋梁建設(shè)是常見的涉河建設(shè)內(nèi)容，橋梁是兩岸物質(zhì)運(yùn)輸、人員流動的重要載體，同時橋梁建設(shè)對河流的擾動影響著河流的度汛安全[1-2]。我國《防洪法》第二十七條中規(guī)定涉河項目建設(shè)應(yīng)符合防洪標(biāo)準(zhǔn)，不得危害堤防安全、影響河勢穩(wěn)定、妨礙行洪暢通等，在實際工程中橋梁橋墩繞流在多方面對河流的防洪安全有著復(fù)雜的影響[3-4]。在發(fā)生歷時洪水時，涉河橋梁對所在河段的防洪安全影響主要包括橋前壅水與水流經(jīng)過橋梁時橋墩繞流改變原有河勢，主要體現(xiàn)在水力因素與泥沙因素[5]。就橋前壅水而言，國內(nèi)外已取得了諸多計算及實際工程經(jīng)驗，在單圓柱與雙圓柱橋墩、橋河正交與斜交，單線、雙線以及多線橋梁布置工程條件下，對其進(jìn)行試驗研究并取得符合相應(yīng)工程條件的計算公式及計算方法[6-11]。在橋墩繞流方面，由橋墩對水流的擾動改變河流原有的河勢，對下游防洪安全產(chǎn)生諸多不可控因素，當(dāng)橋河斜交時相對使下游橋孔的平均流速比同跨類型正交橋增大10%～15%，在計算橋墩基礎(chǔ)埋深時應(yīng)將設(shè)計流速增大10%～15%[12]。同時,由于相鄰橋墩間的壓縮作用使水流向一側(cè)偏轉(zhuǎn)，對河流行洪不利。在斜交角度小于60°時，李付軍等[13]建立的偏角公式與試驗資料吻合較好。在一組橋墩含多個橋墩串聯(lián)時，對水流的導(dǎo)流作用及合適布置形式尚未有明確結(jié)論。對于一般的實際工程水力學(xué)問題，通常采用實體物理模型試驗與計算機(jī)仿真數(shù)值計算兩種方法。實體物理模型的傳統(tǒng)水力學(xué)試驗，其受場地、人力、造價等因素影響，具有不便性和研究周期較長等缺點。數(shù)值方法具有不受場地影響，更加方便快捷等優(yōu)勢。經(jīng)過前人不斷完善流體力學(xué)及計算水力學(xué)的基本理論，隨著計算機(jī)計算性能的提高，以及大量的實際工程經(jīng)驗，得到了很多適用于水力學(xué)湍流計算模型及實際工程的計算軟件。MIKE21 HD中水動力學(xué)模塊適用于河流水力學(xué)計算[14]。HEC-RAS常用于工程水面線的計算，鄧顯羽等[15]運(yùn)用該軟件對密江特大橋進(jìn)行了防洪評價。FLUENT 是目前應(yīng)用比較廣泛的流體力學(xué)商業(yè)計算軟件。何國建等[16]基于FLUENT軟件運(yùn)用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行了橋墩群河道水流的三維數(shù)值分析。本文基于FLUENT軟件，應(yīng)用標(biāo)準(zhǔn)k-ε湍流模型模擬了串聯(lián)三圓柱形橋墩在順直來流條件下的繞流流場，分析了單組串聯(lián)三圓柱橋墩在橋墩串聯(lián)軸線與水流夾角和不同橋墩間距不同排列組合作用下，對水流偏導(dǎo)作用的機(jī)理及由于繞流尾跡偏轉(zhuǎn)致使橋墩下游漩渦阻力分布不均對水流偏導(dǎo)的作用程度。

1 數(shù)值方法

1.1 計算工況

為了探究串聯(lián)圓柱形橋墩對水流的導(dǎo)流作用及不同布置條件下對水流的導(dǎo)流程度，取串聯(lián)圓柱橋墩串聯(lián)軸線與水流夾角和相鄰圓柱橋墩間距兩個控制因素，進(jìn)行排列組合進(jìn)行多種工況模擬分析。在初始流速分別為0.5m/s和1.0m/s兩種來流條件下,對圓柱圓心間距L與直徑D比分別為2.0、2.5、3.0、3.5、4.0，橋墩串聯(lián)軸線與夾角分別為5°、10°、15°、20°、25°、30°組合工況進(jìn)行模擬研究。每組流量30種工況，共計60組。具體尺寸以及角度組合的模擬工況見表1。

表1 模擬計算工況

1.2 模型建立

為了建立更加有效合理的計算區(qū)域模型，整個計算域模型以單個繞流圓柱橋墩直徑為單位，直徑D為2m。模擬水道長300D，寬20D，繞流橋墩直徑為1D，其中間圓柱中心距左側(cè)進(jìn)口100D，串聯(lián)三圓柱形橋墩布置以中間圓柱為中心并按中心對稱排列(見圖1)。

圖1 二維計算水槽模型

1.3 網(wǎng)格劃分

采用三角形網(wǎng)格對整個計算域進(jìn)行劃分，為提高計算精度同時兼顧計算效率，圓柱向四周由密變疏，網(wǎng)格尺寸增長率為1.05，見圖2。同時為了進(jìn)一步提高圓柱繞流邊界層的計算精度和仿真真實性，對圓柱周圍重新劃分網(wǎng)格，采用矩形網(wǎng)格，由邊界向四周的網(wǎng)格尺寸增長率為1.1(見圖3)。

圖2 整體計算網(wǎng)格

圖3 圓柱邊界局部網(wǎng)格

1.4 計算湍流模型

在計算流體力學(xué)方面，經(jīng)過前人的不斷探索改進(jìn)已有不少成熟的湍流計算模型，本文采用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行模擬計算。在標(biāo)準(zhǔn)k-ε模型中輸運(yùn)方程為

Gk+Gb-ρε-YM+SK

(1)

(2)

湍流耗散率計算如下:

式中：Gk為平均速度的湍動能生成率；Gb為浮力產(chǎn)生的湍動能；C1ε、C2ε、C3ε為常數(shù)；σk、σε為有關(guān)湍流的普朗特數(shù)。

2 模型驗證

2.1 繞流水力特性驗證

繞流阻力系數(shù)對比見表2。將雷諾數(shù)Re=2×104對單圓柱繞流阻力系數(shù)與前人的研究成果進(jìn)行比較，對比結(jié)果表明本文數(shù)值計算結(jié)果與文獻(xiàn)[17]、[18]數(shù)據(jù)吻合很好。

表2 繞流阻力系數(shù)對比

2.2 繞流流場結(jié)構(gòu)驗證

對比文獻(xiàn)[19]通過物理水槽模型試驗對單圓柱橋墩繞流所得結(jié)果，本文所得橋墩繞流流場結(jié)構(gòu)(見圖4)與物理試驗結(jié)果流場結(jié)構(gòu)在尾跡擺動、渦體分布規(guī)律上相似。

綜上，本文建立的數(shù)學(xué)方法得到了合理的數(shù)值解，滿足計算求解需要。

圖4 單圓柱繞流尾跡

3 結(jié)果分析

3.1 繞流局部流場結(jié)構(gòu)

當(dāng)液體流經(jīng)圓柱體表面時，所受的阻力有兩種：一是作用在圓柱表面的摩擦阻力，二是漩渦發(fā)生引起的漩渦阻力。串聯(lián)多圓柱串聯(lián)軸線與水流方向成一定夾角時，由于串聯(lián)圓柱形成類似導(dǎo)流墻，對水流具有偏折作用。由圖5可以看出，在圓柱橋墩間，因橋墩間對水流相互壓縮作用，圓柱繞流尾跡向左側(cè)偏轉(zhuǎn)，致使因漩渦產(chǎn)生的漩渦阻力主要集中左側(cè)，在串聯(lián)圓柱左側(cè)形成高阻低流速區(qū)，過流能力驟減，河流主要過流向另一側(cè)集中。因此，圓柱對水流的偏導(dǎo)作用主要由兩方面組成：一是串聯(lián)多圓柱成一定角度時形成類似導(dǎo)流墻對水流的折射作用；二是前方圓柱后的回流區(qū)受后方圓柱的壓縮作用偏向一側(cè)，使漩渦阻力向一側(cè)集中，增大動能損失減小過流流量。

圖5 串聯(lián)三圓柱繞流局部流速矢量

3.2 尾跡偏轉(zhuǎn)

在順直水流條件下，單圓柱的繞流漩渦阻力在圓柱后部，圓柱后部的脫落渦對稱分布或左右交替對稱。在串聯(lián)多圓柱時，由于圓柱間對繞流邊界水流的壓縮作用，串聯(lián)軸線與水流成一定夾角會使圓柱后部的尾跡回流區(qū)發(fā)生一定的偏轉(zhuǎn)，漩渦阻力向一側(cè)集中，如圖5為漩渦阻力區(qū)，水流因漩渦阻力造成的動能損失增大，使過流能力減小。由于橋墩間的壓縮作用對水流的壓縮影響，圓柱繞流分離點左右呈不對稱分布，水流受偏導(dǎo)作用使圓柱后尾跡回流區(qū)偏向一側(cè)。定義左右分離點連線的法向方向與水流方向夾角為繞流尾跡偏轉(zhuǎn)角度β。

3.2.1 第一根橋墩繞流尾跡偏角β1

水流過第一橋墩向兩側(cè)均勻分流后，繞流尾跡受橋墩間壓縮影響，繞流尾跡向一側(cè)偏轉(zhuǎn)。繞流偏轉(zhuǎn)角度與橋墩串聯(lián)軸線與水流方向夾角呈正相關(guān)。在橋墩串聯(lián)軸線與水流方向夾角一定時，橋墩間距在2.0D～4.0D范圍內(nèi)，橋墩間距越小壓縮作用越明顯，偏角越大。根據(jù)模擬結(jié)果(見圖6)，在橋墩間距為4.0D時，偏角不超過2.5°。橋墩串聯(lián)軸線與水流方向夾角為5°時偏角不超過2°。此時繞流尾跡偏轉(zhuǎn)較小，水流受漩渦阻力分布不均引起的偏導(dǎo)作用較小。

圖6 第一根橋墩繞流尾跡偏角

3.2.2 第二根橋墩繞流尾跡偏角β2

相較于第一個橋墩，中間橋墩繞流受前后橋墩的兩方向壓縮作用，繞流尾跡偏角相較于前面橋墩繞流尾跡偏角明顯增大。由圖7可以看出，當(dāng)橋墩間距一定時，在橋墩串聯(lián)軸線與水流方向夾角在5°～16°范圍內(nèi)，繞流尾跡偏角隨橋墩串聯(lián)軸線與水流方向夾角正相關(guān)且增長迅速。在16°～30°范圍內(nèi)，橋墩間距為2.0D時，尾跡偏轉(zhuǎn)角度緩慢增長；橋墩間距為2.5D～4.0D時，串聯(lián)軸線與水流方向夾角在16°左右達(dá)到最大值，尾跡偏轉(zhuǎn)角度開始緩慢減小。另外，在夾角為5°時，尾跡偏角非常小，在5°～10°之間偏轉(zhuǎn)角增加明顯，在10°～15°時增長最為迅速。

圖7 第二根橋墩繞流尾跡偏角

當(dāng)串聯(lián)軸線與水流夾角一定時，在5°～13°范圍內(nèi)，繞流尾跡偏角與橋墩間距呈正相關(guān)；在13°～30°時，繞流尾跡偏角與橋墩間距呈負(fù)相關(guān)；夾角為13°時，在任一橋墩間距下，繞流尾跡偏角相等。

4 結(jié) 論

通過數(shù)值模擬進(jìn)行串聯(lián)三圓柱在橋墩間距2.0D～4.0D范圍內(nèi)，橋墩串聯(lián)軸線與水流夾角對水流偏導(dǎo)作用分析得出如下結(jié)論：

a.串聯(lián)多圓柱橋墩對水流的偏導(dǎo)作用主要表現(xiàn)在：串聯(lián)橋墩形成的類似導(dǎo)流墻對水流的反射作用和橋墩繞流尾跡受壓縮作用向一側(cè)偏轉(zhuǎn)使漩渦阻力向該側(cè)集中，致使水流動能損失增大，過流能力降低。

b.第一個橋墩尾跡偏角與橋墩串聯(lián)線和水流夾角呈正相關(guān)，與橋墩間距負(fù)相關(guān)。

c.第二個橋墩尾跡偏角在橋墩間距2.5D～4.0D范圍內(nèi)，偏角隨橋墩串聯(lián)線和水流夾角先增大后減小。在橋墩串聯(lián)線和水流夾角在5°～13°范圍內(nèi)偏角與橋墩間距正相關(guān)，在13°～30°范圍內(nèi)負(fù)相關(guān)。

本文利用FLUENT軟件對串列三圓柱橋墩在不同橋墩間距以及不同串列軸線與水流方向夾角組合工況進(jìn)行模擬研究。根據(jù)模擬結(jié)果，分析了串列三圓柱橋墩對水流的偏導(dǎo)作用，為涉河橋梁工程規(guī)劃、設(shè)計、建設(shè)以及河道防洪提供參考。