賈 旭
(湖南省建工集團(tuán)有限公司,湖南 長(zhǎng)沙 410000)
根據(jù)臨桂新區(qū)工程設(shè)計(jì),按照防洪標(biāo)準(zhǔn),太平河河道需滿足100年一遇洪水標(biāo)準(zhǔn),原太平河河道寬度將由20 m增加至不小于40 m。同時(shí),為減少新修河道占用周邊農(nóng)田,并與周邊環(huán)境相協(xié)調(diào),新設(shè)計(jì)河道岸線在下游處曲率逐漸變大,與渡槽走向斜交。同時(shí),原渡槽下方,當(dāng)?shù)卮迕褚蓝战ㄓ型ㄐ泻?jiǎn)易板橋需同時(shí)進(jìn)行拆除,為不影響后續(xù)行人通行,設(shè)計(jì)方案中確定將渡槽與橋梁同時(shí)改建,故而,太平河渡槽改建包含渡槽及人行橋兩部分內(nèi)容。重建渡槽位于新設(shè)計(jì)河道彎段,設(shè)計(jì)時(shí)僅考慮槽墩與渡槽中心線垂直,未考慮河道內(nèi)水流方向,為確定槽墩與水流是否順直,本文針對(duì)河道水流流場(chǎng)進(jìn)行模擬分析。
太平河位于桂林市臨桂新區(qū)境內(nèi),太平河渡槽位于河道中部,原渡槽跨度與老河道寬度基本相同,長(zhǎng)度約20 m,渡槽為簡(jiǎn)支梁式,共3跨,4個(gè)槽墩,河道中間槽墩為單矩形柱式,采用漿砌石砌筑,槽身寬度約0.6 m,原渡槽使用年限較久,槽身及槽墩結(jié)構(gòu)已有多處漏水、脫落,加之跨度不滿足新河道寬度,故新渡槽改建方案確定將老渡槽進(jìn)行拆除,并在原址處修建新渡槽。
照前述渡槽改建內(nèi)容,設(shè)計(jì)單位結(jié)合太平河河道設(shè)計(jì)岸線確定了渡槽設(shè)計(jì)初步方案,新建渡槽設(shè)計(jì)過(guò)水流量與原渡槽相同,同樣采用簡(jiǎn)支梁式渡槽,每跨9.6 m,設(shè)4個(gè)槽墩及兩個(gè)邊墩,槽墩由槽身及人行橋共用,兩者并列平行布置。渡槽外空尺寸為0.9 m×1.1 m(寬×高),人行橋橋面寬2.3 m,河道內(nèi)槽墩為適應(yīng)水流采用圓端形墻狀墩身,槽墩尺寸為4.4 m×0.6 m。
槽墩原設(shè)計(jì)初衷是為適應(yīng)水流要求,盡可能減少槽墩對(duì)水流影響,降低水流對(duì)槽墩的沖擊,提高槽墩的穩(wěn)定性。然而,由于渡槽需要在原址上重建,該部位位于新設(shè)計(jì)河道彎段,設(shè)計(jì)時(shí)僅考慮槽墩與渡槽中心線垂直,未考慮河道內(nèi)水流方向,為確定槽墩與水流是否順直,需先對(duì)河道水流流場(chǎng)進(jìn)行模擬。本文選用程凱等二維模型的計(jì)算與實(shí)測(cè)數(shù)據(jù)對(duì)比結(jié)果[1],采用二維模型平面模型對(duì)彎曲復(fù)式河道進(jìn)行模擬計(jì)算。
根據(jù)初步設(shè)計(jì)水文分析資料,太平河河道入口20年一遇洪水流量為121 m3/s,按照河道斷面形式計(jì)算入口處平均流速0.67 m/s,同時(shí)算得河道水力半徑為2.67 m。
基于Wilson等利用激光多譜勒測(cè)速儀(LDV)測(cè)量室內(nèi)動(dòng)床河灣模型中的二次環(huán)流[2],并將測(cè)量結(jié)果與標(biāo)準(zhǔn)k-ξ雙方程模型計(jì)算結(jié)果比較。表明在采用合理網(wǎng)格布置和差分格式的前提下,標(biāo)準(zhǔn)k-ξ紊流模型能夠較好地預(yù)測(cè)彎曲河道中的主流和二次流流速。
直角坐標(biāo)系下,標(biāo)準(zhǔn)k-ξ模型的湍動(dòng)能及其比耗散率輸運(yùn)方程為:
式中:Gk為由層流速度梯度產(chǎn)生的湍流動(dòng)能;Gω為由ω方程產(chǎn)生的湍流動(dòng)能;Tk和Tω為k和ω的擴(kuò)散率;Yk和Yω為由擴(kuò)散而產(chǎn)生的湍流。
FLuent軟件中提供了除標(biāo)準(zhǔn)k-ξ模型外的其他兩種k-ξ模型,包括RNG模型以及Realizable模型,針對(duì)雷諾數(shù)較高的情況(Re>105),結(jié)合標(biāo)準(zhǔn)壁面函數(shù),采用Realizable k-ξ模型,數(shù)值計(jì)算精度較高。
通過(guò)設(shè)置監(jiān)測(cè)點(diǎn)獲得無(wú)槽墩情況下槽墩上游面中心所在位置處流速及速度方向,監(jiān)測(cè)點(diǎn)由左岸至右岸分別為P_1、P_2、P_3、P_4。數(shù)值計(jì)算后速度場(chǎng)及監(jiān)測(cè)點(diǎn)流速及速度方向統(tǒng)計(jì),見(jiàn)圖 1、表 1。
圖1 河道流場(chǎng)圖
表1 監(jiān)測(cè)點(diǎn)數(shù)據(jù)統(tǒng)計(jì)表
根據(jù)以上數(shù)值模擬結(jié)果,在彎段修建渡槽,采用圓端形槽墩,墩柱軸線與流速方向存在一定角度,對(duì)水流條件影響較大。
為緩解上述問(wèn)題,工程提出新的渡槽設(shè)計(jì),即采用對(duì)不均勻流場(chǎng)適應(yīng)性較強(qiáng)的單圓柱形槽墩,槽墩直徑為0.8 m。
新建渡槽設(shè)計(jì)過(guò)水流量與原渡槽相同,同樣采用簡(jiǎn)支梁式渡槽,每跨9.6 m,設(shè)4個(gè)槽墩及兩個(gè)邊墩,槽墩由槽身及人行橋共用,兩者并列平行布置。渡槽外空尺寸為0.9 m×1.1 m(寬×高),人行橋橋面寬2.3 m,河道內(nèi)槽墩為適應(yīng)水流采用圓端形墻狀墩身,槽墩尺寸為4.4 m×0.6 m。
為評(píng)價(jià)新方案單圓柱形槽墩與原設(shè)計(jì)方案對(duì)水流的影響,本文采用與河道流場(chǎng)分析相同的Fluent軟件進(jìn)行數(shù)值模擬。
根據(jù)李彬,孫東坡等研究[3],橋梁與水流斜交且橋墩為圓端形時(shí),橋墩阻水寬度計(jì)算采用下式:
式中:B為橋墩阻水寬度,L為橋軸線法線方向的橋墩寬度,b為橋軸線方向的橋墩寬度,α為橋墩長(zhǎng)軸軸線與水流方向間夾角。
計(jì)算結(jié)果見(jiàn)表2。
表2 槽墩特征長(zhǎng)度
分別建立同條件下各種槽墩的入流模型,各槽墩尺寸不變,圓端形槽墩與水流方向角度取前述確定角度,計(jì)算區(qū)域以槽墩直徑和阻水寬度中的大值D為特征長(zhǎng)度,計(jì)算域的長(zhǎng)寬分別為30D和20D,槽墩中心距上游及下游面分別為10D和20D,距離左右側(cè)均為10D。
數(shù)值模擬時(shí),邊界條件選用速度入流和自由出流;流場(chǎng)的上下邊界采用對(duì)稱邊界條件,以減少邊界對(duì)水流的影響;槽墩設(shè)定為無(wú)滑移壁面條件,入流速度及與槽墩的夾角均采用前述計(jì)算結(jié)果。
求解器為基于壓力的分離式求解器,采用二階隱式時(shí)間格式的非定常流動(dòng)進(jìn)行計(jì)算,壓力與速度的耦合計(jì)算采用PISO算法,壓力的離散為二階精度離散,動(dòng)量方程采用二階迎風(fēng)格式,計(jì)算中壓力、密度、體積和動(dòng)量項(xiàng)的欠松弛因子分別為0.7、1、1、0.7。
阻力是流體對(duì)槽墩的沿流動(dòng)方向的壓力,即為槽墩對(duì)水流的阻滯作用力。根據(jù)阻力系數(shù)公式,計(jì)算槽墩阻力系數(shù),見(jiàn)下式:
式中:Fd為流體作用于單位長(zhǎng)度柱體上沿順流方向的分力。
結(jié)果見(jiàn)表3。
表3 各類型槽墩阻力系數(shù)
經(jīng)過(guò)數(shù)值模擬后的速度流線見(jiàn)圖2~圖6。
圖2 1#槽墩速度矢量圖
圖3 2#槽墩速度矢量圖
圖4 3#槽墩速度矢量圖
圖5 4#槽墩速度矢量圖
圖6 圓柱形槽墩速度矢量圖
1)通過(guò)圖2~圖6可以看出,1#~4#圓端形槽墩偏轉(zhuǎn)一側(cè)出現(xiàn)漩渦,漩渦的位置與流速大小及軸線偏轉(zhuǎn)角度大小均相關(guān),隨著槽墩偏轉(zhuǎn)角度增大,漩渦位置往上游方向偏移,流速越大漩渦面積越大,而新方案中圓柱形槽墩漩渦均位于槽墩下游位置,與槽墩相對(duì)位置固定。
2)通過(guò)表2槽墩阻力系數(shù)表統(tǒng)計(jì)可看出,隨著墻狀槽墩偏移角度的增大,槽墩對(duì)水流阻力系數(shù)呈增大趨勢(shì),而新方案中圓柱型槽墩阻力系數(shù)遠(yuǎn)小于墻狀槽墩,對(duì)水流的阻滯較小,故新方案中圓柱形槽墩水流條件較好。
通過(guò)方案比選,可見(jiàn)采用單圓柱形槽墩在結(jié)構(gòu)穩(wěn)定性及對(duì)水流影響方面均優(yōu)于原方案設(shè)計(jì)中圓端形槽墩。同時(shí),圓柱形槽墩對(duì)彎段水流適應(yīng)性強(qiáng),槽墩基礎(chǔ)設(shè)計(jì)具有更加明顯的針對(duì)性。因此,在河道岸線多變、水流條件復(fù)雜的河道,跨河建筑物推薦采用單圓柱形槽墩或橋墩。