張 晗，張洪福

(東北林業(yè)大學，黑龍江 哈爾濱 150040)

1 水平軸風力發(fā)電機模型

為了控制風致振動，本文將小型水平軸風力發(fā)電機模型放置在橋梁節(jié)段模型上面。風機模型直徑為60 mm，模型的縮尺比為1∶40。每個風力機模型都由一個多葉片轉(zhuǎn)子組成，該轉(zhuǎn)子連接到一個小型交流電(AC)發(fā)電機，葉片連接在小型交流電發(fā)電機上，風力機的重量包括發(fā)電機，葉片和支撐桿等，一系列總重量不超過14 g。風機最大能量變換效率體現(xiàn)在葉尖速比最優(yōu)的情況下，葉尖速比定義公式如下

(1)

式中:ω為旋轉(zhuǎn)角速度，rad·s-1；R為葉片的半徑，m；V為額定風速，m·s-1；n為風輪轉(zhuǎn)速，r·min-1。λ一般與葉片數(shù)量、葉片翼型、類型有關。本實驗選取的葉片葉尖速比為2～3左右，具有較高的功率系數(shù)Cp=0.5左右。

2 風力機尾流特征

本文通過數(shù)值模擬探究單個風力機順流向渦的空間形態(tài)。由于風機的轉(zhuǎn)子對模擬結果影響較大，因此只對風機轉(zhuǎn)子進行數(shù)值模擬。

2.1 計算域的選取

計算模型中的風機采用1∶40的縮尺模型，利用ICEM軟件建模，計算域大小如圖2所示。計算域分為內(nèi)流場和外流場兩部分，內(nèi)流場為風輪轉(zhuǎn)動區(qū)域，外流場為風輪前面的來流區(qū)和風輪后面的尾流區(qū)，為靜止區(qū)域。計算域整體為矩形形狀，內(nèi)流場與外流場之間設置Interface面來保證兩個流場之間的計算數(shù)據(jù)相對流通。整個計算域展向長度為24D(D為風輪直徑)，計算域高度為8D，風機旋轉(zhuǎn)中心到速度入口的距離為5D，風機旋轉(zhuǎn)中心到壓力出口的距離為19D。

2.2 網(wǎng)格劃分

網(wǎng)格劃分后的總數(shù)量為463萬，對風機和內(nèi)流場采用局部加密的方式，風力機表面網(wǎng)格數(shù)量為1.6萬。網(wǎng)格質(zhì)量最小為0.225，最大為1，總體質(zhì)量較好。流域網(wǎng)格采用混合網(wǎng)格形式，內(nèi)流場由于風機葉片旋轉(zhuǎn)角度較為復雜采用非結構化網(wǎng)格形式，外流場為結構化網(wǎng)格。

2.3 邊界條件設置與求解策略

本文采用分離求解器基于壓力型進行求解，求解類型為瞬態(tài)?；陲L力機流場的湍流特性，選擇SSTk-w湍流模型；計算域條件設置內(nèi)流場旋轉(zhuǎn)為滑移網(wǎng)格。邊界條件中，外流場入口速度為7 m·s-1，風機設置為Moving Wall，與內(nèi)流場旋轉(zhuǎn)域同步旋轉(zhuǎn)，轉(zhuǎn)速均為523.8 r·s-1且無滑移。求解方法選擇壓力速度耦合算法SIMPLE算法；為提高計算準確性，采用二階迎風格式。對于殘差收斂設定為10-5，時間步長為0.000 1 s，監(jiān)控風機的升力系數(shù)、阻力系數(shù)和力矩系數(shù)的變化，選擇初始化之后，進行迭代計算。

2.4 數(shù)值模擬結果

風力機在來流風作用下旋轉(zhuǎn)會產(chǎn)生大尺度的順流向渦，該順流向渦以螺旋狀向下游漂移，多個文獻表明該類型渦結構可抑制展向渦的形成與發(fā)展，從而抑制結構的渦激振動。

3 橋梁節(jié)段模型風洞試驗參數(shù)簡介

本次風洞試驗在哈爾濱工業(yè)大學浪槽與風洞聯(lián)合實驗室內(nèi)進行，實驗室是一座閉口回流式風洞，由兩個試驗段組成。本次風洞試驗均是在小風洞實驗段完成，實驗風速范圍從3～50 m·s-1，為連續(xù)可調(diào)。該矩形工作試驗段寬4 m，高3 m，長25 m。在試驗過程中，自由流的湍流強度小于0.46%，自由流的不均勻度小于1%，平均氣流偏角小于0.5°。

4 尾流試驗結果與分析

尾流區(qū)的漩渦脫落是引起橋梁渦激振動的主要原因，因此研究尾流區(qū)的速度變化有著重要的意義。本次實驗分別以裸梁和成橋狀態(tài)下增設欄桿作為基礎工況進行對比。實驗中的風速均采用折減風速U/fB，這里U為來流風速，f為模型的振動頻率，B為模型寬度。風力機通過交錯布置放置在橋梁節(jié)段模型上。其中h是小型風力機的高度，L是風力機展向的間距。風力機沿展向方向間隔為2H、3H(H為橋梁模型的高度，H=0.1m)，旋轉(zhuǎn)軸的高度分別為h/H=0.28、0.38。

選取風力機旋轉(zhuǎn)軸高度為h/H=0.28，展向間隔為3H的實驗工況研究尾流風速的變化。為了與模型懸掛狀態(tài)下的尾流數(shù)據(jù)進行對比，本次實驗選取彈簧懸掛系統(tǒng)狀態(tài)下的折減風速為U/fvB=1.328和U/fαB=1.430。監(jiān)測點所用儀器為熱線風速儀和二維熱線探頭，采樣頻率為1 000 Hz。針對以上兩種風速，分別對橋梁節(jié)段模型順風向(U)進行監(jiān)測分析。

無控和安裝小型水平軸風力發(fā)電機模型的狀態(tài)下監(jiān)測點位置圖。監(jiān)測點分布在直線z/H=0以及z/H=1.5，y/H=0，這里z/H=0，y/H=0對應安裝的風機模型的位置，z/H=1.5，y/H=0對應兩個風力機模型的中間位置。沿著x方向上共有六個監(jiān)測點x=4H、8H、12H、16H、20H、24H。y軸為無量綱順流向脈動風速為u′/U0，u′為順流向脈動風速，U0為來流風速。

無控狀態(tài)下和有控狀態(tài)尾流區(qū)順流向脈動風速值隨著順流向距離的增大而逐漸衰減，展向渦能量逐漸耗散。豎向渦振情況，在近尾流區(qū)無控要比有控的脈動風速值低；在有控的工況下，順流向風速脈動值要衰減的快，遠尾流區(qū)有控比無控狀態(tài)下脈動風速值低，說明在有風機的情況下旋轉(zhuǎn)產(chǎn)生的順流向渦能夠抑制展向渦的形成與脫落，加快順流向脈動風速值的降低。扭轉(zhuǎn)渦振情況，在近尾流區(qū)有控和無控相比，順流向脈動風速值要偏低，在遠尾流區(qū)有控和無控順流向脈動風速值相差不是太大。

尾流測試結果：功率譜密度

功率譜可以反映主梁漩渦脫落的頻率與強度。一般而言，橋梁旋渦脫落若能夠被抑制，則渦激共振現(xiàn)象也可被抑制。本文中，選取尾流區(qū)監(jiān)測位置位于直線z/H=0和z/H=1.5，y/H=0，x/H=20上，展向間距為3H，旋轉(zhuǎn)軸高度為h/H=0.28的風機布置在橋梁節(jié)段模型的迎風端和背風端。

無控狀態(tài)時，在折減風速U/fvB=1.328處出現(xiàn)強烈的豎向渦激共振，展向渦能量非常大。在折減風速U/fαB=1.430出現(xiàn)強烈的扭轉(zhuǎn)渦激共振。當布置小型水平軸風力發(fā)電機時，峰值能夠完全消失，表明安裝小型水平軸風力發(fā)電機對橋梁節(jié)段模型渦激共振有明顯的控制效果。

5 結 論

水平軸風力機是一種有效收集風能的裝置，其葉片旋轉(zhuǎn)可以產(chǎn)生順流向渦結構。主梁展向渦是導致渦振產(chǎn)生的主要原因，本文采以多葉片水平軸風力機為研究對象，采用節(jié)段模型風洞試驗探究了風力機尾渦對典型大跨橋梁主梁展向渦的抑制作用，采用數(shù)值模擬觀察了風力機尾流形態(tài)。數(shù)值模擬結果表明，風力機的尾流為順流向螺旋狀渦結構。試驗結果表明在豎向渦振發(fā)生風速范圍內(nèi)，近尾流區(qū)受葉片旋轉(zhuǎn)影響，風力機增加了風速的脈動，但在遠尾流區(qū)降低了風速脈動；在扭轉(zhuǎn)渦振發(fā)生范圍內(nèi)，遠近尾流區(qū)的脈動均被減小；風速頻譜分析結果表明風力機可以有效抑制主梁周期性渦脫，從而減弱渦激振動。