摘 要：采用用戶自定義函數(shù)（UDF）中三大宏編程更新墩臺(tái)的質(zhì)心位置和振動(dòng)狀態(tài)，通過(guò)建立墩臺(tái)渦激振動(dòng)數(shù)值模型并結(jié)合嵌套網(wǎng)格技術(shù)，對(duì)有無(wú)前置隔板墩臺(tái)的水動(dòng)力特性進(jìn)行分析，結(jié)果表明：設(shè)置前置隔板后，約化速度為6～7時(shí)墩臺(tái)的回流區(qū)順流向長(zhǎng)度和渦流長(zhǎng)度隨著約化速度增大顯著減小，而無(wú)前置隔板時(shí)墩臺(tái)表現(xiàn)出相反趨勢(shì)，前置隔板有效抑制了尾流擾動(dòng)。約化速度為7～9時(shí)，有前置隔板墩臺(tái)的回流區(qū)順流向長(zhǎng)度和渦流長(zhǎng)度比無(wú)前置隔板墩臺(tái)的小，前置隔板減弱了渦激振動(dòng)對(duì)墩臺(tái)的影響。前置隔板有效延遲了旋渦脫落時(shí)間，墩臺(tái)上下表面壓力分布更為均勻，減弱了墩臺(tái)的渦激振動(dòng)響應(yīng)。

關(guān)鍵詞：UDF；前置隔板；墩臺(tái)；渦激振動(dòng)抑制；嵌套網(wǎng)格

中圖分類號(hào)：O357.5+2 文獻(xiàn)標(biāo)志碼：A doi：10.3969/j.issn.1000-1379.2025.02.024

引用格式：施佳科，李東風(fēng)，陳姿妙，等.基于UDF的前置隔板墩臺(tái)渦激振動(dòng)抑制數(shù)值分析[J].人民黃河，2025，47（2）：156-160.

基金項(xiàng)目：浙江省自然科學(xué)基金資助項(xiàng)目（LZJWZ22C030001，LZJWZ22E090004）；浙江省水利廳水利科技項(xiàng)目（RB2115）；國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2016YFC0402502）；國(guó)家自然科學(xué)基金資助項(xiàng)目（51979249）

NumericalAnalysisofVortexInducedVibrationSuppression ofFrontPartitionPierBasedonUDF

SHIJiake1，2，LIDongfeng2，CHENZimiao2，WULuobin2，F(xiàn)ANGZhihao1，3，BAIFuqing2，HUJianyong2

（1.SchoolofNavalArchitectureandMaritime，ZhejiangOceanUniversity，Zhoushan316022，China；2.SchoolofHydraulicEngineering，ZhejiangUniversityofWaterResourcesandElectricPower，Hangzhou310018，China；3.ShanghaiZhongshenIndustryCo.，Ltd.，Shanghai200131，China）

Abstract：UserDefinedFunctions（UDF）withthreemajormacroprogrammingwasusedtoupdatethecentroidpositionandvibrationstateof thepier，anumericalmodelofvortexinducedvibrationofthepierwasestablishedandcombinedwithnestedgridtechnologytoanalyzethe hydrodynamiccharacteristicsofthepierwithandwithoutprepartitionplates.Theresultsshowthataftersettingupthefrontpartition，the downstreamlengthandvortexlengthofthebackflowzoneofthepiersignificantlydecreasewiththeincreaseofthereducedvelocitywhenthe reducedvelocityis6-7，whilethepierwithoutaprebaffleshowstheoppositetrend.Theprebaffleeffectivelysuppressesthewakedisturb？ ance.Whenthereducedspeedisbetween7and9，thelengthofthedownstreamandvortexintherecirculationzoneofthepierwithafront baffleissmallerthanthatofthepierwithoutafrontbaffle.Thefrontbafflereducestheimpactofvortexinducedvibrationonthepier.The frontpartitioneffectivelydelaysthevortexsheddingtime，andthepressuredistributionontheupperandlowersurfacesofthepierismoreu？ niform，weakeningthevortexinducedvibrationresponseofthepier.

Keywords：UDF；frontpartition；pier；vortexinducedvibrationsuppression；nestinggrid

在水工建筑物設(shè)計(jì)中，河道、溢洪道和消力池等水道中的閘墩、分隔墩、消力墩等多以墩臺(tái)形式呈現(xiàn)[1-2]，水流繞過(guò)墩臺(tái)后產(chǎn)生渦激振動(dòng)，會(huì)引起墩臺(tái)結(jié)構(gòu)破壞[3-4]。學(xué)者們對(duì)海洋工程、機(jī)械工程和航空航天工程等領(lǐng)域的流體流動(dòng)旋渦抑制研究較為深入[5-7]，但在水利工程領(lǐng)域的研究不多[8]。流體流動(dòng)旋渦抑制分為主動(dòng)抑制和被動(dòng)抑制，其中主動(dòng)抑制需要對(duì)水流運(yùn)動(dòng)系統(tǒng)額外輸入能量，雖然可以抑制旋渦，但額外增加了能源動(dòng)力費(fèi)用，對(duì)水工建筑物的運(yùn)行和管理也會(huì)帶來(lái)一定困難。相較而言，被動(dòng)抑制只需要優(yōu)化設(shè)計(jì)墩臺(tái)的幾何外形，或添加隔板等附屬擾流裝置，就可以抑制流體分離和渦激振動(dòng)[1-2]。Zdravkov？ ich[9]對(duì)各種流體流動(dòng)旋渦被動(dòng)抑制技術(shù)進(jìn)行了全面總結(jié)；Gao等[10]對(duì)設(shè)有不同長(zhǎng)度前置隔板的圓柱體進(jìn)行了風(fēng)洞實(shí)驗(yàn)，發(fā)現(xiàn)前置隔板長(zhǎng)度等于圓柱體直徑時(shí)，與無(wú)前置隔板相比，其平均阻力系數(shù)、均方根升力系數(shù)分別減小36%、63.6%。然而，以往研究尚未揭示抑制水工建筑物圓形墩臺(tái)渦激振動(dòng)的機(jī)理。

本文以Khalak等[11]開(kāi)展的經(jīng)典實(shí)驗(yàn)為基礎(chǔ)，構(gòu)建墩臺(tái)渦激振動(dòng)數(shù)值模型，設(shè)置計(jì)算域邊界條件，應(yīng)用UDF方法求解墩臺(tái)雙向流固耦合問(wèn)題，通過(guò)將相同幾何動(dòng)力參數(shù)的無(wú)前置隔板墩臺(tái)渦激振動(dòng)數(shù)值模型模擬結(jié)果與Khalak等[11]得出的經(jīng)典實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證數(shù)值模型的精確性，分析有無(wú)前置隔板墩臺(tái)尾跡變化情況、旋渦脫落情況、瞬時(shí)壓力系數(shù)分布，進(jìn)而揭示墩臺(tái)渦激振動(dòng)抑制機(jī)理。

1 墩臺(tái)渦激振動(dòng)數(shù)值模型構(gòu)建和驗(yàn)證

1.1 流場(chǎng)控制方程

1.3 模型定解條件

建立墩臺(tái)渦激振動(dòng)數(shù)值模型（見(jiàn)圖1），計(jì)算域尺寸為20D×40D（其中D為墩臺(tái)直徑，取值0.0381m）。設(shè)置計(jì)算域邊界條件與Khalak等[11]開(kāi)展的經(jīng)典實(shí)驗(yàn)邊界條件相同，即水流流入計(jì)算域的入口邊界為速度進(jìn)口邊界，出口邊界為壓力出口邊界，水流與墩臺(tái)的接觸面為固壁邊界（無(wú)滑移壁面邊界），計(jì)算域的頂部和底部均為對(duì)稱邊界。將來(lái)流方向設(shè)置為x方向，垂直來(lái)流方向設(shè)置為y軸方向。把結(jié)構(gòu)初始速度和位移均為0作為初始條件。

墩臺(tái)中心位于距離速度進(jìn)口邊界和兩個(gè)對(duì)稱邊界10D處，前置隔板長(zhǎng)度、厚度分別取D、0.06D[12]。實(shí)際來(lái)流速度U∞為0.0762～0.1524m/s時(shí)，對(duì)有無(wú)前置隔板的墩臺(tái)渦激振動(dòng)進(jìn)行數(shù)值分析。上述實(shí)際來(lái)流速度范圍對(duì)應(yīng)約化速度Ur范圍為5～10，其表示結(jié)構(gòu)振動(dòng)一個(gè)周期內(nèi)流體路徑與結(jié)構(gòu)特征尺寸的比值，是速度的無(wú)量綱參數(shù)。

1.4 嵌套網(wǎng)格劃分與插值處理

為分析隔板與墩臺(tái)周?chē)乃w流動(dòng)特征，采用嵌套網(wǎng)格技術(shù)將整個(gè)計(jì)算域劃分成背景網(wǎng)格（網(wǎng)格數(shù)量為41936個(gè)），將與墩臺(tái)同心半徑為3D的圓作為伴隨運(yùn)動(dòng)區(qū)域。對(duì)網(wǎng)格進(jìn)行插值處理，組分網(wǎng)格最外層邊界被定義為Overset邊界，以便其與背景網(wǎng)格之間進(jìn)行流場(chǎng)信息傳輸，插值前后網(wǎng)格對(duì)比見(jiàn)圖2。

1.5 UDF實(shí)現(xiàn)原理

在Fluent中每個(gè)剛體都有定義的質(zhì)心。DEFINE_CG_MOTION是用戶自定義函數(shù)的一種宏，用于定義剛體的質(zhì)心運(yùn)動(dòng)。函數(shù)控制剛體質(zhì)心的運(yùn)動(dòng)方式，通過(guò)計(jì)算線速度和角速度，間接更新剛體位置，描述剛體在流場(chǎng)中的運(yùn)動(dòng)規(guī)律。結(jié)構(gòu)振動(dòng)會(huì)導(dǎo)致流體邊界隨時(shí)間變化，UDF可定義邊界隨振動(dòng)結(jié)構(gòu)的運(yùn)動(dòng)規(guī)律。DEFINE_GRID_MOTION是一種用于定義網(wǎng)格動(dòng)態(tài)運(yùn)動(dòng)的宏，通過(guò)直接修改網(wǎng)格中的節(jié)點(diǎn)坐標(biāo)，實(shí)現(xiàn)復(fù)雜的邊界運(yùn)動(dòng)和流固耦合，靈活模擬動(dòng)態(tài)網(wǎng)格。UDF還可以自定義輸出，提取研究所需數(shù)據(jù)。DEFINE_EXECUTE_AT_END是一種執(zhí)行用戶自定義計(jì)算任務(wù)的宏。在每個(gè)仿真時(shí)間步結(jié)束時(shí)，采用Fluent調(diào)用龍哥庫(kù)塔四階方程函數(shù)，以便更新墩臺(tái)的振動(dòng)狀態(tài)（位移和速度）。

1.6 雙向流固耦合求解

雙向流固耦合求解過(guò)程見(jiàn)圖3。

1.7 模型預(yù)測(cè)精度驗(yàn)證

對(duì)墩臺(tái)結(jié)構(gòu)參數(shù)進(jìn)行無(wú)量綱設(shè)定：質(zhì)量比為10，阻尼比為0.00542，固有頻率為0.4Hz。通過(guò)將相同幾何動(dòng)力參數(shù)的無(wú)前置隔板墩臺(tái)渦激振動(dòng)數(shù)值模型預(yù)測(cè)結(jié)果與Khalak等[11]得出的經(jīng)典實(shí)驗(yàn)結(jié)果進(jìn)行比較，驗(yàn)證數(shù)值模型的預(yù)測(cè)精度。選擇最大振幅比與振動(dòng)頻率進(jìn)行對(duì)比，結(jié)果見(jiàn)圖4[圖（b）中fn為結(jié)構(gòu)固有頻率，fb為結(jié)構(gòu)實(shí)際振動(dòng)頻率，當(dāng)兩者一致時(shí)會(huì)產(chǎn)生共振]，預(yù)測(cè)結(jié)果與經(jīng)典實(shí)驗(yàn)結(jié)果（定量和變化趨勢(shì)）較一致，說(shuō)明數(shù)值模型具有較高的預(yù)測(cè)精度。

2 前置隔板墩臺(tái)渦激振動(dòng)抑制機(jī)理分析

2.1 墩臺(tái)尾跡變化

圖5為來(lái)流約化速度Ur=7時(shí)有無(wú)前置隔板墩臺(tái)的速度均方根（Urms）等值線，其可反映墩臺(tái)尾跡旋渦的強(qiáng)度和范圍。圖5中回流區(qū)順流向長(zhǎng)度Lr代表時(shí)均順流向速度最小值所處位置與墩臺(tái)中心的距離，對(duì)稱分布的2個(gè)Urms峰值代表剪切層旋滾位置，渦流長(zhǎng)度Lv代表墩臺(tái)中心和剪切層旋滾位置之間的流向距離，尾流寬度W代表2個(gè)Urms峰值之間的縱向距離。分析圖5可知，速度均方根等值線在尾跡軸上呈現(xiàn)出良好的對(duì)稱性。墩臺(tái)的存在以及剪切層的分離，使尾跡區(qū)產(chǎn)生較為明顯的低流速區(qū)域?；亓鲄^(qū)末端與尾跡軸上速度最大值所處位置密切對(duì)應(yīng)。有前置隔板墩臺(tái)的Lr、Lv值分別比無(wú)前置隔板墩臺(tái)的降低了3.0%、2.2%。從x=0.5D、y=0.2D處剪切層附近的等值線可以看出，有前置隔板墩臺(tái)的水體流速梯度較小，抑制了分離流的擴(kuò)散效應(yīng)，剪切層之間的距離更近。

有無(wú)前置隔板墩臺(tái)的尾跡參數(shù)見(jiàn)表1。Ur為6～7（實(shí)際來(lái)流速度為0.09144～0.10668m/s）時(shí)，有前置隔板墩臺(tái)的Lv和Lr值隨著約化速度增大顯著減小，而無(wú)前置隔板墩臺(tái)的Lv、Lr值隨著約化速度增大顯著增大，說(shuō)明在此速度區(qū)間內(nèi)前置隔板有效抑制了尾跡的擾動(dòng)。Ur為7～9時(shí)，有前置隔板墩臺(tái)的Lv、Lr值均小于無(wú)前置隔板墩臺(tái)的?；亓鲄^(qū)較短時(shí)，渦流強(qiáng)度低，渦激振動(dòng)幅度小。渦流長(zhǎng)度較短時(shí)，渦流作用時(shí)間不夠持久，產(chǎn)生的振動(dòng)幅度較小。因此，設(shè)置前置隔板減弱了渦激振動(dòng)對(duì)墩臺(tái)的影響。

2.2 墩臺(tái)周期性旋渦脫落

圖6為無(wú)前置隔板墩臺(tái)的渦量等值線，可反映墩臺(tái)尾跡旋渦的周期性脫落情況。y/D<0的渦量等值線反映0～T/2時(shí)段的旋渦脫落情況（T為旋渦脫落周期），y/D>0的渦量等值線反映T/2～T時(shí)段的旋渦脫落情況。分析圖6可知，墩臺(tái)兩側(cè)剪切層中顯示出較高的渦流強(qiáng)度，剪切層不斷向下游移動(dòng)。t=0時(shí)刻的渦量等值線顯示了從底部剪切層開(kāi)始的逆時(shí)針旋渦，隨著時(shí)間延長(zhǎng)，旋渦范圍增大。0～T/2時(shí)段，渦脫后渦量逐漸減小，同時(shí)向下游移動(dòng)并擴(kuò)展，渦量等值線在渦流中分布較為均勻，渦流渦量耗散較慢。當(dāng)時(shí)間超過(guò)T/2時(shí)，來(lái)自上部剪切層的渦流出現(xiàn)夾帶y/D=0以下旋轉(zhuǎn)方向相反的渦流趨勢(shì)。當(dāng)時(shí)間超過(guò)3T/4時(shí)，即將脫離的渦旋有被墩臺(tái)尾跡另一側(cè)水流切斷趨勢(shì)。

有前置隔板墩臺(tái)的渦量等值線見(jiàn)圖7。分析可知，與無(wú)前置隔板墩臺(tái)相比，0～T/2時(shí)段有前置隔板墩臺(tái)時(shí)其尾部旋渦在發(fā)展時(shí)受流體干擾較持久，導(dǎo)致脫落時(shí)尺寸較小，渦結(jié)構(gòu)受到明顯抑制，表明前置隔板延緩了流體在墩臺(tái)后部的旋渦脫落行為。渦量等值線在渦流中分布較不均勻，前置隔板削弱渦量較為明顯。因此，設(shè)置前置隔板墩臺(tái)增加了旋渦脫落前流體的動(dòng)能耗散，減少了流場(chǎng)對(duì)墩臺(tái)的作用力。

2.3 瞬時(shí)壓力系數(shù)分布

圖8為有無(wú)前置隔板墩臺(tái)表面的瞬時(shí)壓力系數(shù)Cp分布情況（圖中x/D為-1.5～-0.5時(shí)對(duì)應(yīng)于順流向前置隔板結(jié)構(gòu)，x/D為-0.5～0.5時(shí)對(duì)應(yīng)于順流向裸墩臺(tái)結(jié)構(gòu)）。對(duì)于無(wú)前置隔板墩臺(tái)，約化速度為5時(shí)Cp分布呈現(xiàn)傾斜的“8”字形，墩臺(tái)上下表面Cp差值顯著，導(dǎo)致振動(dòng)響應(yīng)增強(qiáng)；約化速度增大到8和10時(shí)，墩臺(tái)上下表面Cp分布較一致，從而減弱了振動(dòng)響應(yīng)。對(duì)于有前置隔板墩臺(tái)，約化速度為8、10時(shí)x/D為-1.5～-1.0（隔板前半部分）處墩臺(tái)上表面Cp大于下表面Cp，x/D為-1.0～-0.5（隔板后半部分）處墩臺(tái)下表面Cp大于上表面Cp，兩種方向相反的水流力作用在有前置隔板墩臺(tái)表面，對(duì)總水流升力幾乎沒(méi)有貢獻(xiàn)。x/D為-1.5～0.5（從隔板起始位置到隔板尾部）處，隔板上的Cp差呈現(xiàn)先減小后增大再減小趨勢(shì)。原因是前置隔板上的水體流動(dòng)類似于在細(xì)長(zhǎng)鈍體上的流動(dòng)，在隔板前緣與墩臺(tái)交界處存在明顯的回流區(qū)，使墩臺(tái)前端速度波動(dòng)，導(dǎo)致前置隔板上存在過(guò)渡流[13]。約化速度為10時(shí)，有前置隔板墩臺(tái)上下表面壓力分布均勻，沒(méi)有明顯的壓力差，導(dǎo)致升力幾乎為0，從而抑制了渦激振動(dòng)。

3 結(jié)論

1）通過(guò)建立渦激振動(dòng)數(shù)值模型，采用嵌套網(wǎng)格技術(shù)實(shí)現(xiàn)了流場(chǎng)信息傳輸，基于UDF三大宏編程實(shí)現(xiàn)網(wǎng)格動(dòng)態(tài)運(yùn)動(dòng)以及更新墩臺(tái)的質(zhì)心位置和振動(dòng)狀態(tài)。通過(guò)將模型預(yù)測(cè)結(jié)果與經(jīng)典試驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證了模型具有較高的預(yù)測(cè)精度。

2）與無(wú)前置隔板墩臺(tái)相比，有前置隔板墩臺(tái)改善了尾部的水動(dòng)力特性。約化速度為6～7時(shí)，前置隔板有效抑制了尾流擾動(dòng)。約化速度為7～9時(shí)，有前置隔板時(shí)墩臺(tái)回流區(qū)順流向長(zhǎng)度和渦流長(zhǎng)度比無(wú)前置隔板墩臺(tái)時(shí)小，前置隔板減弱了渦激振動(dòng)對(duì)墩臺(tái)的影響。

3）前置隔板延遲了旋渦脫落發(fā)生時(shí)間，墩臺(tái)上下表面壓力分布較為均勻，導(dǎo)致升力幾乎為0，有效減少了渦激振動(dòng)帶來(lái)的不利影響。

參考文獻(xiàn)：

[1] 閆秀秀，尹進(jìn)步，張曜，等.導(dǎo)流洞復(fù)合式出口消能體型的研究與應(yīng)用[J].人民黃河，2024，46（4）：157-162.

[2] 趙雪萍，李璐，趙玉良，等.溝水坡水庫(kù)溢洪道彎道水流及消能實(shí)驗(yàn)研究[J].人民黃河，2022，44（5）：128-131，146.

[3] 吳新，鄒志風(fēng).旋轉(zhuǎn)方式對(duì)并列雙圓柱繞流特性的影響[J].人民黃河，2017，39（11）：71-75.

[4] 黃友明.水力學(xué)的基本理論及其在土木工程中的應(yīng)用：評(píng)《水力學(xué)（第三版）》[J].人民黃河，2024，46（4）：163.

[5] 孔令輝，劉東，張雷，等.水力機(jī)械抗空蝕涂層研究進(jìn)展[J].人民黃河，2023，45（6）：157-162.

[6] 楊卓媛，夏軍強(qiáng)，宋紅波，等.大型水庫(kù)運(yùn)用對(duì)河口尾閭河道出汊的影響研究[J].人民黃河，2024，46（4）：37-42.

[7] 梁艷潔，趙正偉.黃河口漁洼河段河勢(shì)及治理方案研究[J].人民黃河，2023，45（11）：65-69.

[8] 郭博文，高玉琴，范冰，等.基于LES的水工鋼閘門(mén)流激振動(dòng)數(shù)值模擬研究[J].人民黃河，2024，46（3）：126-131.

[9] ZDRAVKOVICHMM.ReviewandClassificationofVarious AerodynamicandHydrodynamicMeansforSuppressing VortexShedding[J].JournalofWindEngineeringandIndus？ trialAerodynamics，1981，7（2）：145-189.

[10] GAODL，CHENGB，HUANGYW，etal.FlowCharacter？ isticsofaFixedCircularCylinderwithanUpstreamSplitter Plate：OnthePlate？LengthSensitivity[J].Experimental ThermalandFluidScience，2020，117：110135.

[11] KHALAKA，WILLIAMSONHKC.Motions，F(xiàn)orcesandMode TransitionsinVortex？InducedVibrationsatLowMass？ Damping[J].JFluidsStruct，1999，13（7-8）：813-851.

[12] WANGEH，ZHAOSH，XUWH，etal.EffectofSplitter PlateLengthonFIVofCircularCylinder[J].International JournalofMechanicalSciences，2023，254：108413.

[13] CHUTKEYK，SURIYANARAYANANP，VENKATAKRISHNAN L.NearWakeFieldofCircularCylinderwithaForward SplitterPlate[J].JournalofWindEngineeringIndustrial Aerodynamics，2018，173：28-38.

【責(zé)任編輯 栗 銘】