孫開旗 肖玉德

（安徽交通職業(yè)技術學院，安徽 合肥 230051）

0 引言

洪水沖毀橋梁事故在2011 年，中國北部的巴彥塔拉河鐵路橋被洪水沖毀；2020 年7 月有400 年歷史的安徽的屯溪老大橋被沖毀；2020 年8 月陜西洛南暴雨引發(fā)洪水沖垮26座橋梁；2020 年9 月韓國江原道平昌郡珍富面松亭橋被洪水沖垮。因此洪水對橋梁的危害性較大，也是橋梁沖毀的罪魁禍首之一，需要對洪水和橋墩之間的相互作用進行分析，才可以發(fā)揮橋墩建設的作用。

橋墩結構包括方形墩、片狀墩、圓頭墩以及圓柱墩等結構形式，而在中、小橋梁中，串列雙圓形柱橋墩廣泛應用于橋梁下部結構橋墩中。水流在流動過程中，過流斷面會出現(xiàn)收縮的現(xiàn)象，此時柱體周圍的水流速度會隨著變化而發(fā)生變化，不斷增加速度后，黏性力會出現(xiàn)問題，此時也形成了圓柱繞流的現(xiàn)象。柱體的出現(xiàn)造成了雍水的現(xiàn)象，因此水流自身也會對柱體造成作用力，此過程也是導致作用力復雜的主要因素，同時也會對流體產(chǎn)生較大的干擾，造成不可忽視的影響。

該文基于CFD 的方法，采用TM_LES 湍流模型，對不同間距的串列雙圓形橋墩進行數(shù)值模擬，分析了氣動特性系數(shù)等特征參數(shù)及水場流動特性，以期對實際橋梁工程中水中橋墩的間距布置提供參考。

1 氣動系數(shù)

水流及風本質上均是流體，水流對橋墩的氣動系數(shù)，可以參照風對建筑物的氣動系數(shù)定義。由平均風作用引起的靜荷載稱為靜力風荷載[2]。在橫向風作用下，對建筑物可以用靜力三分力表示，相應的氣動力系數(shù)在風軸上稱為阻力系數(shù)、升力系數(shù)和扭轉力矩系數(shù)。

式中：FD、FL、FT分別為阻力、升力和扭轉力矩；CD、CL、CM分別為阻力系數(shù)、升力系數(shù)和扭轉力矩系數(shù)；ρ為空氣密度；U為來流風速；D為建筑物特征高度；B為建筑物特征寬度。

2 工程實例

2.1 計算模型

計算模型圖如圖1 所示，上、下游圓形橋墩直徑均為D，兩個橋墩中心距為L，且2 個圓形橋墩的幾何中心在同一水平線上，考慮橋墩不同間距L/D=2.0、2.3、2.5、2.7、3.0、3.5、4.0、4.5、5.0m。

圖1 計算模型圖

該文入口邊界為速度入口，來流速度為U=10 m/s，出口邊界為壓力出口，上下邊界為對稱滑移邊界，兩個圓形橋墩為無滑移壁面邊界條件。流體為25℃水，ρ取1.0kg/m3，動黏性系數(shù)為0.001003，湍流強度2.1%。

計算工具采用同濟大學劉十一博士開發(fā)的“小牛CFD”程序，該程序的特點如下：首先是校該文件的優(yōu)勢。這類校本為自主研發(fā)校本，需要引擎的支持，在使用前期，需要明確腳本的語言，盡可能地保證語法簡明扼要，提高整體功能性。在語法方面需要按照C++/C#的方式進行處理，如果涉及矩陣操作問題，就需要將其與Matlab 進行聯(lián)合，此時才可編寫較為復雜的程序，軟件內部也可提供較為復雜的模板形式。其次，針對軟件內部具體情況，對網(wǎng)格進行劃分，保證建模的靈活性需求。在輸入幾何模型后，可以自動生產(chǎn)網(wǎng)絡，用戶需要按照網(wǎng)格的引導針對不同的區(qū)域設置網(wǎng)格密度，建模后還可形成曲線或者直線對創(chuàng)建區(qū)域進行設定，對創(chuàng)建區(qū)域需要應用布爾運算進行操作。再次，其中還包括高效求解器。這類設備可以結合實際需求，采用有限體積法形成最快的線性方程組求解方式，以滿足高校并行和流暢性方面的需求，最終滿足精度。除此之外，還擅長使用動網(wǎng)格、流固耦合計算方式。這類計算方式可以滿足獨特性網(wǎng)格模型需求，并且使網(wǎng)格變形最大化。獨創(chuàng)的求解方式在計算方面整體效率會提升。如果選擇使用流固耦合計算方式，在計算的過程中需要使用最新的流固耦合算法，針對腳本中出現(xiàn)的各類語言進行定義，保證功能性需求，還可以實現(xiàn)復雜化的系統(tǒng)模擬。最后，完善處理功能。此時可以選擇使用邊計算邊顯示的方式，還可針對相關數(shù)據(jù)進行錄像。自定義圖樣可以直接顯示物體的壓強，后續(xù)使用“速度矢量紋理”的對流線進行顯示，使其具備直觀性特點；為了滿足速度需要，還可以采用Direct 3D 顯卡加速技術。

2.2 網(wǎng)格劃分

網(wǎng)格劃分是較為重要的一項環(huán)節(jié)，可以適應特征變化，并且還可滿足加密需求，無論梯度較小或者較大都可形成均勻過渡的需要[3-4]。

傳統(tǒng)迭代求解方法（如雅可比法、共軛梯度法）能快速抹平相鄰節(jié)點的殘差，但是對大尺度的殘差抹平效率很低；以SOR（超松弛）迭代法為例，隨著殘差梯度減小，殘差收斂速度越來越慢。對長周期殘差，迭代效率非常低；求解大規(guī)模線性方程組的效率很低。該工具采用代數(shù)多重網(wǎng)格迭代法，針對不同的尺度同時建立線性方程組，交替執(zhí)行“殘差抹平”操作。這類方式的優(yōu)點是不同的周期也可以達到快速地抹平的需求。

代數(shù)多重網(wǎng)格沒有依賴性，可以結合矩陣結構生成多重網(wǎng)格。小尺度相鄰的頂點可以進行合并，還可切割成為2 個區(qū)域，保證頂點的相似相等，減少切面的邊數(shù)。

代數(shù)多重計算的流程如下：首先針對流體時間獲得原始線性方程組；其次，進行代數(shù)網(wǎng)格的自動劃分，結合不同尺寸生成線性方程組；最后交替進行“殘差抹平”迭代計算。

程序自動進行的網(wǎng)格劃分圖如圖2 所示。

圖2 網(wǎng)格劃分圖

2.3 計算工況

考慮雙圓形橋墩不同間距工況：L/D=2.0、2.3、2.5、2.7、3.0、3.5、4.0、4.5、5.0，共9 個不同間距變化工況。為了結果對比，也對單圓形橋墩在相同條件下進行了數(shù)值模擬。

3 數(shù)值模擬結果

通過數(shù)值模擬，得到不同間距工況下阻力系數(shù)見表1、圖3。

圖3 不同間距工況下阻力系數(shù)圖

表1 不同間距工況下阻力系數(shù)表

對串列雙圓柱橋墩在不同間距時的數(shù)值模擬結果進行分析，可得到以下結論。

對上游圓柱來說，阻力系數(shù)整體呈上升趨勢，并最終趨向于單圓柱情況。在2.0 ≤L/D≤3.0 時呈增長趨勢，數(shù)值變化趨勢較快；在3.0 ≤L/D≤5.0 時，數(shù)值變化趨勢不大，隨著間距比的增大，阻力系數(shù)趨向于單圓柱情況。

對下游圓柱來說，阻力系數(shù)整體呈下降趨勢，并最終趨向于單圓柱情況。在2.0 ≤L/D≤3.0 時呈逐漸降低趨勢，數(shù)值變化趨勢也較快；在3.0 ≤L/D≤5.0 時，總體數(shù)值也呈降低趨勢，但變化較平緩，隨著間距比的加大，阻力系數(shù)也趨向于單圓柱情況。

間距比在2.5 ≤L/D≤3.0，阻力系數(shù)發(fā)生了明顯的跳躍，說明氣動力跳躍臨界間距范圍在2.5 ≤L/D≤3.0。

由上文可知，間距比在2.5 ≤L/D≤3.0，阻力系數(shù)發(fā)生明顯地跳躍，研究當L/D=2.7 時，不同的水流速度（U=5、10、15、20、25、30、35、40 m/s）對阻力系數(shù)的影響，見表2。

表2 L/D=2.7 時不同水流速度下阻力系數(shù)表

通過對串列雙圓柱橋墩在L/D=2.7 時不同水流速度下數(shù)值模擬結果，由圖4 可知：對上游圓柱來說，當水流速度為5.0 ≤U≤15.0 m/s 時，阻力系數(shù)有所波動；當水流速度U＞15.0 m/s 時，阻力系數(shù)變化不大，并最終趨向于0.60。對下游圓柱來說，阻力系數(shù)整體呈下降趨勢，水流速度在5.0 ≤U≤15.0 m/s 時，阻力系數(shù)波動較大；水流速度U＞15.0 m/s，阻力系數(shù)變化較平緩，并最終也趨向于0.60。

圖4 L/D=2.7 時不同水流速度下阻力系數(shù)

上游圓柱阻力系數(shù)受水流速度影響較小，而下游圓柱阻力系數(shù)受水流速度影響相對較大，當水流速度達到40 m/s 以上時，兩者阻力系數(shù)趨于穩(wěn)定值0.60。

4 流場特性

通過數(shù)值模擬計算，得到不同間距下的雙圓柱橋墩速度流線圖，限于篇幅原因，該文只列出間距比為2.0、2.5、3.0、3.5、4.0、5.0 的速度流線圖，如圖5 所示。從圖5 可知，隨著間距比增大，上、下游圓柱橋墩都出現(xiàn)會周期性脫落的漩渦，最后形成湍流，此狀態(tài)接近單圓柱的流場形態(tài)。

圖5 速度流線圖

5 結論

該文基于CFD 的方法，采用TM_LES 湍流模型，對不同間距的串列雙圓形橋墩及串列雙圓柱橋墩在L/D=2.7 時不同水流速度下進行數(shù)值模擬，分析氣動特性系數(shù)等特征參數(shù)及水場流動特性，得到以下結論，為實際橋梁工程中水中橋墩的間距布置提供參考：1）串列雙圓柱存在針對氣動系數(shù)的臨界間距，間距比在2.5 ≤L/D≤3.0 時，阻力系數(shù)發(fā)生了明顯的跳躍，說明臨界間距范圍在2.5 ≤L/D≤3.0 內。2）對上游圓柱，阻力系數(shù)整體呈上升趨勢，并最終趨向于單圓柱情況。在2.0 ≤L/D≤3.0 時呈增長趨勢，數(shù)值變化趨勢較快；在3.0 ≤L/D≤5.0 時，數(shù)值變化趨勢不大，隨著間距比的增大，阻力系數(shù)趨向于單圓柱情況。

對下游圓柱，阻力系數(shù)整體呈下降趨勢，并最終趨向于單圓柱情況。在2.0 ≤L/D≤3.0 時呈逐漸降低趨勢，數(shù)值變化趨勢也較快；在3.0 ≤L/D≤5.0 時，總體數(shù)值也呈降低趨勢，但變化較平緩，隨著間距比的加大，阻力系數(shù)也趨向于單圓柱情況。3）在臨界間距范圍內，上游圓柱阻力系數(shù)受水流速度影響較小，而下游圓柱阻力系數(shù)受水流速度影響相對較大，當水流速度達到40 m/s 以上時，兩者阻力系數(shù)趨于穩(wěn)定值0.60。4）隨著間距比增大，上、下游圓柱橋墩均出現(xiàn)會周期性脫落的漩渦，最后形成湍流，該狀態(tài)接近單圓柱的流場形態(tài)。