王貴春，曹宗恒

(鄭州大學(xué) 土木工程學(xué)院 河南 鄭州 450001)

0 引言

渦激振動是由結(jié)構(gòu)尾流區(qū)漩渦脫落引起的結(jié)構(gòu)周期性振動，具有強(qiáng)迫、自激及限幅的特點(diǎn)。斜拉索在流體作用下都可能發(fā)生渦激振動，從而導(dǎo)致結(jié)構(gòu)疲勞破壞，影響結(jié)構(gòu)正常使用。所以研究斜拉索渦激振動響應(yīng)特性具有重要意義。

Feng[1]在1968年進(jìn)行風(fēng)洞試驗(yàn)，研究了單圓柱體渦激振動響應(yīng)。隨著計(jì)算機(jī)技術(shù)的迅速發(fā)展，很多學(xué)者采用計(jì)算流體動力學(xué)(compute fluid dynamics,CFD)方法研究圓柱體渦激振動響應(yīng)特性。李廣望[2]采用任意Lagrange-Euler方法計(jì)算了圓柱體的渦激振動響應(yīng)。徐楓[3]采用動網(wǎng)格技術(shù)對不同截面形式的柱體進(jìn)行了數(shù)值模擬，但是柱體運(yùn)動引起的計(jì)算網(wǎng)格畸變、纏繞等問題不能得到有效解決。劉志文等[4]采用“剛性區(qū)域+動網(wǎng)格區(qū)域+靜止網(wǎng)格區(qū)域”的方式建立網(wǎng)格，對矩形斷面橋梁渦激振動進(jìn)行計(jì)算，計(jì)算過程中雖能保證壁面附近計(jì)算網(wǎng)格的質(zhì)量，但是運(yùn)動網(wǎng)格劃分方式較復(fù)雜且不同劃分方式網(wǎng)格間的連接容易出現(xiàn)問題，尤其是當(dāng)網(wǎng)格出現(xiàn)變形時(shí)，動網(wǎng)格區(qū)域計(jì)算網(wǎng)格質(zhì)量大幅下降，容易影響計(jì)算結(jié)果的精度。劉軍[5]采用Fluent對高雷諾數(shù)下輸電導(dǎo)線渦激振動進(jìn)行了分析，并詳細(xì)分析了高質(zhì)量比和低質(zhì)量比參數(shù)對渦激振動的影響，其研究結(jié)果表明：低質(zhì)量比渦激振動振幅要大于高質(zhì)量比渦激振動振幅。曲寶旭[6]采用重疊網(wǎng)格法對具有彈性支撐不同自由度下輸電導(dǎo)線的渦激振動進(jìn)行了分析，結(jié)果表明：不同自由度對渦激振動橫風(fēng)向振動具有一定影響。

當(dāng)前對圓柱渦激振動數(shù)值模擬主要是采用動網(wǎng)格技術(shù)，采用重疊網(wǎng)格法的相關(guān)報(bào)道相對較少。本文采用重疊網(wǎng)格技術(shù)建立網(wǎng)格，應(yīng)用Newmark-β法求解結(jié)構(gòu)振動方程，并且編制UDF(user defined functions)，嵌入到Fluent中，通過宏命令DEFINE_ZONE_MOTION傳遞運(yùn)動速度，以更新運(yùn)動網(wǎng)格位置，對斜拉索結(jié)構(gòu)渦激振動響應(yīng)進(jìn)行數(shù)值模擬，從而解決結(jié)構(gòu)振動過程中可能產(chǎn)生的網(wǎng)格畸變和負(fù)體積問題，實(shí)現(xiàn)流固耦合計(jì)算。

1 理論基礎(chǔ)

1.1 流體控制方程

不可壓縮流體N-S方程和連續(xù)性方程分別為

(1)

(2)

式中：υ=μ/ρ為分子黏度，μ為動力黏度；fi為單位質(zhì)量流體所受到的體積力(質(zhì)量力)；p為流體介質(zhì)微元體上的壓力；ρ表示流體密度；ui表示x、y、z方向的流體速度分量。

1.2 結(jié)構(gòu)振動控制方程

以文獻(xiàn)[7]中斜拉索實(shí)驗(yàn)?zāi)Ｐ蜑楸尘斑M(jìn)行斜拉索渦激振動分析，斜拉索直徑D設(shè)置為0.066 m。如圖1所示，斜拉索可簡化為兩自由度質(zhì)量-彈簧-阻尼系統(tǒng)，振動方程表示為

圖1 兩自由度質(zhì)量-彈簧-阻尼體系

(3)

2 數(shù)值模型建立

在劃分網(wǎng)格之前，要確定計(jì)算域的大小，合理的計(jì)算域大小可以使數(shù)值計(jì)算效率提高，也能保持計(jì)算結(jié)果的準(zhǔn)確性。通常可以根據(jù)阻塞率的大小來確定計(jì)算域大小，在CFD模型中阻塞率一般控制在5%以內(nèi)。二維流場計(jì)算域大小為24D×50D，其中D為圓柱的直徑，如圖2所示。定義圓柱截面形心所在位置為計(jì)算域坐標(biāo)原點(diǎn)，設(shè)定重疊區(qū)域直徑大小為3D，坐標(biāo)原點(diǎn)到入口邊界(INLET)的水平距離為15D，與出口邊界(OUTLET)的水平距離為30D。上下兩側(cè)均為對稱邊界(SYMMETRY)，其到坐標(biāo)原點(diǎn)的垂直距離均為12D。拉索表面為無滑移壁面邊界(NO-SLIPWALL)。在進(jìn)行拉索渦激振動計(jì)算時(shí)，湍流模采用剪切應(yīng)力輸運(yùn)模型(shear stress transition，SST)k-ω模型，湍流度設(shè)置為5%。該模型能較好地處理壁面附近的流動得到可靠的結(jié)果，適用于存在逆壓梯度和流動分離的圓柱繞流場計(jì)算[8]。

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

2.1 重疊網(wǎng)格劃分及邊界條件

重疊網(wǎng)格由背景網(wǎng)格和前景網(wǎng)格組成，也即流體計(jì)算域根據(jù)計(jì)算需要被劃分為多個(gè)部件區(qū)域。首先在各區(qū)域單獨(dú)創(chuàng)建網(wǎng)格，然后通過插值計(jì)算實(shí)現(xiàn)相鄰計(jì)算域間的流場數(shù)據(jù)交換。與傳統(tǒng)動網(wǎng)格計(jì)算相比，重疊網(wǎng)格在計(jì)算過程中，不會產(chǎn)生網(wǎng)格纏繞、畸變和負(fù)體積等問題，能始終保持計(jì)算網(wǎng)格的質(zhì)量，從而保證計(jì)算結(jié)果的準(zhǔn)確性。

采用CFD軟件Fluent中的重疊網(wǎng)格(overset mesh)技術(shù)，進(jìn)行圓柱體兩自由度渦激振動響應(yīng)計(jì)算。重疊網(wǎng)格技術(shù)的實(shí)施，包括兩個(gè)主要步驟：1)在不同網(wǎng)格間建立通信關(guān)系；2)利用從步驟1中獲得的信息，進(jìn)行數(shù)值求解。

2.2 網(wǎng)格獨(dú)立性與時(shí)間無關(guān)性檢驗(yàn)

表1 網(wǎng)格方案詳細(xì)參數(shù)設(shè)置

從圖3(a)可以看出在不同網(wǎng)格方案設(shè)置下，拉索氣動響應(yīng)隨著時(shí)間步的減小都能達(dá)到一個(gè)穩(wěn)定的值，不同的是網(wǎng)格數(shù)量越大以及寬高比越小，也就是網(wǎng)格的整體質(zhì)量越好(A1,A2網(wǎng)格)，得到的阻力系數(shù)均值較A3網(wǎng)格方案得到的結(jié)果稍大。從圖3(b)中可以看出，在不同網(wǎng)格方案下，采用A1和A2網(wǎng)格得到的St數(shù)在整體趨勢上更加吻合，結(jié)果也比較接近?？梢缘贸鼋Y(jié)論，采用A1和A2網(wǎng)格方案，得到St數(shù)與其他學(xué)者的結(jié)果比較一致，如表2所示。可以認(rèn)為A2、A3網(wǎng)格方案已通過網(wǎng)格獨(dú)立性以及時(shí)間無關(guān)性檢驗(yàn)，繼續(xù)對網(wǎng)格進(jìn)行加密意義不大。綜上所述，在考慮計(jì)算精確度和計(jì)算效率下，選擇A2網(wǎng)格方案和時(shí)間步大小為0.002 s進(jìn)行數(shù)值模擬。

表2 Re=3 900時(shí)St數(shù)對比

圖3 阻力系數(shù)均值與St數(shù)隨時(shí)間步大小的變化

3 斜拉索渦激振動響應(yīng)特性分析

拉索在實(shí)際服役時(shí)，渦激振動起振的折減風(fēng)速為Ur=U/(fn·D)，通常小于10，所以選取的計(jì)算折減風(fēng)速Ur范圍為3.5～7.5，相應(yīng)的雷諾數(shù)Re位于6 500～14 000之間。

3.1 St數(shù)及渦脫頻率分析

圖4(a)中可以看出，在折減風(fēng)速Ur=3.5～4.64之間時(shí)，St數(shù)隨折減風(fēng)速變化呈線性下降的趨勢；當(dāng)Ur=4.727時(shí)，St=0.203；在Ur=4.64～4.94之間時(shí)，隨著Ur的增大St呈倒“V”字形變化，即St數(shù)先突然增大，又急劇減小，此時(shí)對應(yīng)拉索發(fā)生渦激振動的折減風(fēng)速區(qū)間，可以推斷這個(gè)區(qū)間內(nèi)流場變化比較劇烈；當(dāng)Ur>4.94時(shí)，St數(shù)又突然增大之后呈線性減小趨勢。如圖4(b)所示，隨著折減風(fēng)速增大，由于拉索發(fā)生了頻率鎖定現(xiàn)象(lock-in)，所以在某一風(fēng)速區(qū)間內(nèi)拉索渦脫頻率并沒有線性增大；當(dāng)渦脫頻率約為3.494 Hz時(shí)，對應(yīng)的折減風(fēng)速為4.94；當(dāng)Ur>4.59時(shí)，拉索渦脫頻率沒有隨著折減風(fēng)速的增大而增大，而是與在靠近拉索的一階固有頻率，所以，可以認(rèn)為此處為鎖定區(qū)間的起點(diǎn)；當(dāng)Ur=5.04時(shí)，拉索渦脫頻率與結(jié)構(gòu)固有頻率解鎖，渦脫頻率隨折減風(fēng)速的增加開始呈線性變化。從圖4(b)中可以看出，拉索發(fā)生鎖定的區(qū)間比較窄，也就說，拉索渦激振動對風(fēng)速的變化非常敏感。

圖4 不同折減風(fēng)速下渦激振動St數(shù)與升力系數(shù)頻率分析結(jié)果

3.2 振幅響應(yīng)分析

圖5給出了不同折減風(fēng)速下拉索橫風(fēng)向無量綱位移時(shí)程曲線。從圖5中可以看出，在高質(zhì)量比下拉索在折減風(fēng)速較低時(shí)，位移呈現(xiàn)微幅簡諧振動，拉索的振動幅度也很小。隨著折減風(fēng)速的增加拉索振動幅度逐漸增大，并且出現(xiàn)了“拍振”現(xiàn)象。橫向位移在折減風(fēng)速Ur=4.67～4.94之間顯著增加。

圖5 典型折減風(fēng)速下拉索橫向振動位移時(shí)程曲線

如圖6(a)所示，可以看出，在鎖定區(qū)間拉索橫風(fēng)向振動頻率接近拉索固有頻率。在折減風(fēng)速Ur=4.64～4.94之間時(shí)，漩渦脫落頻率與拉索固有頻率值非常接近，即拉索發(fā)生了“鎖定”，拉索作大幅振動，并且在折減風(fēng)速Ur=4.94時(shí)，拉索達(dá)到最大振幅0.20D左右(如圖6(b)所示)。

圖6 不同折減風(fēng)速下拉索橫向主振頻率與最大位移

3.3 拉索渦激振動軌跡分析

圖7(a)～(g)給出了不同折減風(fēng)速下拉索渦激振動的軌跡圖。拉索渦激振動是以空氣為介質(zhì)進(jìn)行計(jì)算，拉索的質(zhì)量比m*=159.4，此時(shí)拉索應(yīng)屬于高質(zhì)量比下渦激振動。高質(zhì)量比情況下，渦激振動的運(yùn)動軌跡形式更加復(fù)雜。在較低或者較高折減風(fēng)速下，拉索振動軌跡主要呈橢圓形式，同時(shí)運(yùn)動的方向和位置也在不斷發(fā)生變化。在接近結(jié)構(gòu)固有頻率時(shí)，拉索在鎖定區(qū)內(nèi)也即流固耦合作用較強(qiáng)時(shí)，拉索運(yùn)動軌跡大致呈“8”字形。低質(zhì)量比圓柱渦激振動軌跡一般呈單一的“8”字形，高質(zhì)量比結(jié)構(gòu)渦激振動的運(yùn)動形式更加豐富，在高質(zhì)量比參數(shù)下拉索不再是單一的“8”字形運(yùn)動，這與文獻(xiàn)[6]的結(jié)論一致。

圖7 拉索渦激振動質(zhì)心運(yùn)動軌跡

4 質(zhì)量比對渦激振動影響分析

本節(jié)設(shè)定質(zhì)量比m*=7.9作為低質(zhì)量比情況，阻尼比保持ζ=0.546%不變(相對應(yīng)的質(zhì)量阻尼m*ζ=0.041 5)，對拉索進(jìn)行數(shù)值模擬。

從圖8中可以看到，在質(zhì)量比為159.4時(shí)，升阻力系數(shù)的均值是比較接近的，也就是說，在高質(zhì)量比參數(shù)下，流場的變化并不是很強(qiáng)烈。對比不同質(zhì)量比下拉索升阻力系數(shù)均值可以發(fā)現(xiàn)，低質(zhì)量比參數(shù)下拉索升阻力系數(shù)均值變化的范圍更廣，橫風(fēng)向振動大小對升阻力大小有較大的影響，在結(jié)構(gòu)振幅較大時(shí)，其升阻力也有一個(gè)非常顯著的增大過程，但隨著振幅的降低，其升阻力又重新降低，直到保持一個(gè)穩(wěn)定的值?？梢酝茢?，低質(zhì)量比參數(shù)下流場與結(jié)構(gòu)之間的相互作用也更加顯著。

圖8 不同折減風(fēng)速及質(zhì)量比下升阻力系數(shù)均值對比

如圖9、10所示，相較于高質(zhì)量比情況，拉索渦激振動的鎖定區(qū)間也有所變大，由鎖定區(qū)過渡到非鎖定區(qū)較慢。當(dāng)頻率比為1時(shí)，拉索達(dá)到最大振幅為0.82D，為高質(zhì)量比參數(shù)下最大振幅的4.1倍。在設(shè)置計(jì)算的有限工況內(nèi)，渦激振動處于鎖定區(qū)間時(shí)的鎖定頻率比在1.02左右，也就是說，鎖定頻率發(fā)生了小幅度的偏移。當(dāng)頻率比達(dá)2.318時(shí)，渦激振動振幅顯著減小至0.05D，此時(shí)渦激振動已脫離鎖定區(qū)間。

圖9 不同折減風(fēng)速下橫風(fēng)向最大位移變化

對比圖8～10可以推斷：質(zhì)量比越小，結(jié)構(gòu)發(fā)生渦激振動時(shí)鎖定區(qū)間越寬，最大振幅也越大。隨著質(zhì)量比減小，發(fā)生鎖定時(shí)的頻率比也會有一定的偏移。

圖10 不同折減風(fēng)速下頻率比變化

5 結(jié)論

通過對斜拉索進(jìn)行渦激振動數(shù)值模擬，得出了4個(gè)結(jié)論。1)拉索發(fā)生渦激振動時(shí)，阻力系數(shù)均值首先降低至最小值，而后又顯著增加到最大值。之后隨著折減風(fēng)速的繼續(xù)增大，阻力系數(shù)均值呈不斷減小趨勢。升力系數(shù)均方根值隨折減速度的增大(Ur=3.5～4.51)基本保持不變，之后升力系數(shù)突然降低，又在Ur=4.59～4.94之間保持顯著增大的趨勢，直至最大值，之后隨折減風(fēng)速的增加呈減小趨勢；2)對于高質(zhì)量比渦激振動，其振幅雖然較低，但由于振動產(chǎn)生的疲勞破壞要特別注意。只有在接近拉索一階固有頻率的范圍內(nèi)，拉索振幅顯著增大，達(dá)到了0.2D，之后雖折減風(fēng)速進(jìn)一步增加，但振幅迅速減小，拉索振幅的變化對于折減風(fēng)速的變化非常敏感；3)對拉索渦激振動的運(yùn)動軌跡進(jìn)行了詳細(xì)的分析，拉索屬于高質(zhì)量比渦激振動，其運(yùn)動軌跡更加豐富，只有在接近發(fā)生渦激振動的范圍內(nèi)，拉索大致呈“8”字形運(yùn)動軌跡，這與低質(zhì)量比圓柱的渦激振動軌跡有顯著差異；4)質(zhì)量比會對拉索渦激振動產(chǎn)生一定的影響，是影響渦激振動的重要因素。通過對低質(zhì)量比渦激振動分析，可以得出結(jié)論：質(zhì)量比越小，振動鎖定區(qū)間越大，拉索在鎖定區(qū)間內(nèi)振幅也越大。