李久偉 路建華 周華僑 黃賜榮

（中國(guó)電建集團(tuán)華東勘測(cè)設(shè)計(jì)研究院，浙江 杭州 400045）

0 引言

風(fēng)能是一種分布廣泛、儲(chǔ)量豐富的新能源，大力發(fā)展風(fēng)電產(chǎn)業(yè)對(duì)調(diào)整能源結(jié)構(gòu)、推進(jìn)能源生產(chǎn)具有重要意義。根據(jù)國(guó)際可再生能源署統(tǒng)計(jì)，截至2021 年底，風(fēng)能發(fā)電在新能源發(fā)電量的占比已超過(guò)25%，同比漲幅高達(dá)13%。在產(chǎn)業(yè)政策引導(dǎo)和市場(chǎng)需求驅(qū)動(dòng)的雙重作用下，我國(guó)陸上風(fēng)電裝機(jī)容量已躍居世界第一位，成為全球陸上風(fēng)電開(kāi)發(fā)與利用的主要推動(dòng)力。

在電價(jià)補(bǔ)貼持續(xù)退坡、平價(jià)上網(wǎng)的行業(yè)背景下，發(fā)展大單機(jī)容量、大葉輪直徑和高塔筒的大型高空化風(fēng)電機(jī)組是降低度電成本、提高收益率的有效技術(shù)手段，已成為陸上風(fēng)電主流發(fā)展方向。

目前商業(yè)化應(yīng)用程度最廣的風(fēng)電塔架是柔性全鋼塔架，其具有標(biāo)準(zhǔn)化程度高、產(chǎn)業(yè)鏈成熟、施工快等優(yōu)點(diǎn)。為進(jìn)一步實(shí)現(xiàn)發(fā)電量的最大化提升，柔塔的整體剛度以及結(jié)構(gòu)固有頻率隨輪轂高度的增加而下降，使柔塔渦激共振所需的風(fēng)速進(jìn)一步降低，更容易觸發(fā)渦激振動(dòng)，增加倒塔事故風(fēng)險(xiǎn)，因此如何保障柔塔的安全性就成了行業(yè)亟待解決的問(wèn)題。該文從空氣動(dòng)力學(xué)的原理出發(fā)，分析了風(fēng)電塔架渦激振動(dòng)的原理，并結(jié)合理論計(jì)算等多種研究手段，給出了渦激振動(dòng)響應(yīng)的不同計(jì)算方式，最后提出了不同類(lèi)型的防治措施。

1 塔筒的渦振原理

1.1 起振風(fēng)速及鎖定現(xiàn)象

當(dāng)氣流流過(guò)塔筒的表面時(shí)，在流體黏性和逆壓梯度的作用下，流體在塔筒尾部?jī)蓚?cè)會(huì)交替生成旋向相反、周期性脫落的旋渦，隨著風(fēng)速增大，渦脫速度和頻率也會(huì)相應(yīng)增大，當(dāng)渦脫頻率增大到與塔筒固有頻率接近時(shí)，塔筒就會(huì)發(fā)生橫風(fēng)向的振動(dòng)現(xiàn)象，即渦激振動(dòng)。

塔筒發(fā)生渦振的臨界風(fēng)速被稱(chēng)為起振風(fēng)速vcr，該風(fēng)速與結(jié)構(gòu)的自振頻率，塔筒直徑等結(jié)構(gòu)參數(shù)有關(guān)，該風(fēng)速的計(jì)算見(jiàn)式（1）。

式中：為塔筒的直徑，f為塔筒的第階自振頻率，S為斯特勞哈爾數(shù)。

當(dāng)渦脫頻率與結(jié)構(gòu)固有頻率接近時(shí)，渦脫頻率便會(huì)鎖定在結(jié)構(gòu)的固有頻率上，如圖1 所示，當(dāng)風(fēng)速大于起振風(fēng)速時(shí)，塔筒振動(dòng)又會(huì)反過(guò)來(lái)改變流場(chǎng)，此時(shí)渦脫頻率不會(huì)隨風(fēng)速的增大而增大，在一定風(fēng)速區(qū)間內(nèi)，渦脫頻率會(huì)持續(xù)與結(jié)構(gòu)自振頻率保持一致，結(jié)構(gòu)會(huì)持續(xù)發(fā)生渦激振動(dòng)，該現(xiàn)象被稱(chēng)為渦激振動(dòng)的“鎖定現(xiàn)象”。

圖1 風(fēng)電塔筒渦激振動(dòng)

1.2 渦激共振區(qū)高度

由于沿塔筒高度方向的風(fēng)速有差異，因此當(dāng)風(fēng)剖面處于鎖定風(fēng)速區(qū)間內(nèi)時(shí)，塔筒才會(huì)發(fā)生渦振。嚴(yán)格來(lái)說(shuō)，塔筒發(fā)生渦激振動(dòng)是由于部分高度區(qū)域的渦脫頻率與結(jié)構(gòu)自振頻率相近，因此導(dǎo)致整個(gè)塔筒發(fā)生共振現(xiàn)象，該高度區(qū)間被稱(chēng)為“渦振共振區(qū)”。

由起振風(fēng)速公式（1）可知，起振風(fēng)速與塔筒自振頻率成正比，當(dāng)塔筒結(jié)構(gòu)越柔，自振頻率越小，其發(fā)生渦激振動(dòng)的臨界風(fēng)速也會(huì)相應(yīng)更小，更可能發(fā)生渦激振動(dòng)。而混凝土結(jié)構(gòu)塔筒由于其剛度較大，自振頻率也會(huì)更大，因此在同等高度下混凝土塔筒不容易發(fā)生渦振。

2 塔筒渦振響應(yīng)計(jì)算

在渦激振動(dòng)的相關(guān)研究中，塔筒的自振頻率、起振風(fēng)速對(duì)塔筒的安全性至關(guān)重要。在渦振領(lǐng)域另一個(gè)主要研究方向?yàn)闇u振的響應(yīng)計(jì)算，常見(jiàn)的渦振響應(yīng)計(jì)算方法有理論計(jì)算、風(fēng)洞模型試驗(yàn)以及數(shù)值模擬，該文將介紹三種渦振響應(yīng)的計(jì)算方法，并對(duì)計(jì)算過(guò)程中的重難點(diǎn)進(jìn)行介紹。

2.1 理論計(jì)算方法

塔筒在風(fēng)荷載作用下，會(huì)發(fā)生橫風(fēng)向的渦激振動(dòng)，針對(duì)塔筒建立運(yùn)動(dòng)方程，其中、、分別為順風(fēng)向、橫風(fēng)向、高度方向。塔筒的渦振位移由式（2）所控制。

式中：（，）為高度處塔筒的無(wú)量綱位移，即高度處塔筒的位移與直徑之比；為塔筒結(jié)構(gòu)阻尼比；λ為質(zhì)量比，λ=/2m，其中為空氣密度，為高度處塔筒的單位長(zhǎng)度質(zhì)量；U為無(wú)量綱風(fēng)速，U=/（D），為結(jié)構(gòu)圓頻率；為尾流振子角位移，其定義如圖2 所示；為無(wú)量綱風(fēng)速U與無(wú)量綱起振風(fēng)速U，=U/U；為與的比值，其中為尾流振子扭轉(zhuǎn)中心到質(zhì)心的距離，如圖2 所示；為尾流振子的阻尼比，其通過(guò)=/（2 2 π），其中為與Magnus 效應(yīng)有關(guān)的參數(shù)。

圖2 塔筒尾流振子示意圖[5]

求解方程組式（2），可推導(dǎo)出塔筒在z 高度處的有效渦激力系數(shù)幅值C。

式中：（）為高度處的風(fēng)速，（）為高度處的塔筒直徑，C（）為高度處的有效渦激力系數(shù)。

采用伽遼金方法，可以將上述運(yùn)動(dòng)方程（2）的解分解為伽遼金函數(shù)與廣義位移時(shí)程的乘積，見(jiàn)式（5）。

式中：為伽遼金函數(shù)的階數(shù)，φ為方向的振型函數(shù)，q為方向的廣義位移。

將式（4）、式（5）代入式（2），并乘以φ后沿高度方向積分，可求得塔筒的渦激振動(dòng)方程為式（6）。

由式（3）、式（6）可以看出有效渦激力系數(shù)幅值C和位移幅值相關(guān)，因此對(duì)某一風(fēng)速下的塔筒渦激振動(dòng)位移，可以通過(guò)迭代法進(jìn)行求解。

2.2 風(fēng)洞試驗(yàn)測(cè)量

對(duì)于渦激振動(dòng)的響應(yīng)，塔筒的結(jié)構(gòu)參數(shù)、塔筒結(jié)構(gòu)的局部剛度變化、風(fēng)場(chǎng)的湍流特性等都會(huì)在極大程度上影響渦激振動(dòng)的響應(yīng)，因此對(duì)該文的研究對(duì)象風(fēng)電機(jī)組塔筒，只能進(jìn)行縮尺試驗(yàn)。為了得到塔筒的渦振響應(yīng)以及起振風(fēng)速，可以制作塔筒的氣彈模型進(jìn)行風(fēng)洞試驗(yàn)，來(lái)獲得渦激振動(dòng)的位移響應(yīng)以及起振風(fēng)速。塔筒的氣彈模型設(shè)計(jì)除滿(mǎn)足外形幾何相似、剛度和質(zhì)量分布相似外，還需要慣性參數(shù)、彈性參數(shù)、重力參數(shù)、黏性參數(shù)、斯特羅哈數(shù)、結(jié)構(gòu)衰減率等的一致性。

對(duì)于氣彈模型風(fēng)洞試驗(yàn)，塔筒的局部連接等很難滿(mǎn)足縮尺比，一般可通過(guò)模型制作材料來(lái)調(diào)節(jié)模型頻率來(lái)滿(mǎn)足縮尺相似關(guān)系。風(fēng)洞試驗(yàn)中塔筒順風(fēng)向以及橫風(fēng)向的風(fēng)致加速度可以采用加速度傳感器來(lái)監(jiān)測(cè)獲取，風(fēng)洞實(shí)際風(fēng)速由風(fēng)洞中的風(fēng)速測(cè)速儀實(shí)測(cè)得到。由于塔筒風(fēng)洞試驗(yàn)時(shí)風(fēng)速較大，無(wú)接觸式圖像識(shí)別技術(shù)可以避免對(duì)風(fēng)洞內(nèi)部的流場(chǎng)產(chǎn)生干擾，可以獲取可靠的高精度位移時(shí)程數(shù)據(jù)，塔筒模型的渦振位移同樣可以采用無(wú)接觸式圖像識(shí)別技術(shù)進(jìn)行采集。近些年來(lái)，無(wú)接觸式圖像識(shí)別技術(shù)憑借無(wú)接觸、精度高、干擾小、成本低的優(yōu)點(diǎn)，已經(jīng)成為位移測(cè)量的主要方式和手段，未來(lái)無(wú)接觸式圖像識(shí)別技術(shù)具有廣泛的應(yīng)用前景。

2.3 渦激振動(dòng)數(shù)值模擬

目前，計(jì)算流體動(dòng)力學(xué)（CFD）已逐步成為模擬塔筒周邊的漩渦脫落情況的有效方法之一。塔筒的渦振也可以利用CFD 技術(shù)進(jìn)行計(jì)算，空氣和塔筒的相互作用是流固耦合作用，塔筒受到氣動(dòng)力的作用發(fā)生振動(dòng)，而振動(dòng)反過(guò)來(lái)會(huì)對(duì)流場(chǎng)產(chǎn)生影響。為進(jìn)一步分析塔筒的繞流特性，需要開(kāi)展塔筒流固耦合計(jì)算。

流固耦合的計(jì)算步驟主要有網(wǎng)格劃分、計(jì)算設(shè)置、LES模型流場(chǎng)求解、提取塔筒表面氣動(dòng)力、Newmark-β 法迭代求解。流體與固體的求解過(guò)程按上述順序求解，流場(chǎng)計(jì)算所得氣動(dòng)力作為固體求解的輸入條件，而固體求解所得位移、速度等響應(yīng)將傳輸給流體域網(wǎng)格，通過(guò)網(wǎng)格的移動(dòng)實(shí)現(xiàn)固體位移變化，依次迭代直至完成計(jì)算。

對(duì)塔筒的數(shù)值模擬，為了對(duì)比試驗(yàn)縮尺模型計(jì)算結(jié)果，基于CFD 技術(shù)做縮尺模型模擬，測(cè)量塔筒頂部的渦振位移以及加速度。通過(guò)風(fēng)洞試驗(yàn)與數(shù)值模擬得到的位移、加速度對(duì)比來(lái)驗(yàn)證數(shù)值模擬結(jié)果的準(zhǔn)確性。

3 塔筒渦振防治措施研究

塔筒避免渦激振動(dòng)的關(guān)鍵在于如何避免與結(jié)構(gòu)固有頻率接近，因此該文從塔筒結(jié)構(gòu)形式、氣動(dòng)力干擾以及阻尼耗能的角度出發(fā)，分別給出不同的渦振減振思路。

3.1 結(jié)構(gòu)式防治措施

通過(guò)該文第1 節(jié)的相關(guān)研究可以看出，相比鋼混組合式塔筒，全鋼塔筒剛度較小，更容易發(fā)生渦激振動(dòng)，且隨著塔筒高度的高度達(dá)到160M 級(jí)，鋼塔的渦振發(fā)生率也在大幅度提高。在同等高度下，混凝土塔筒不容易發(fā)生渦振，因此為了減小塔筒結(jié)構(gòu)發(fā)生渦振的可能性，混凝土塔筒更加具有優(yōu)勢(shì)，而在一定程度上，改變結(jié)構(gòu)的形式會(huì)增加構(gòu)設(shè)計(jì)成本。但在塔架高度和葉輪直徑的雙重增加下，鋼混塔筒可能隨著高度的快速增大而面臨渦激振動(dòng)的相關(guān)困擾，因此從流場(chǎng)干擾、阻尼耗能的角度考慮，研發(fā)相關(guān)的渦振減振裝置非常迫切。

3.2 氣動(dòng)式防治措施

改變塔筒周?chē)鲌?chǎng)也是減小渦激振動(dòng)的有效措施。從流場(chǎng)的角度來(lái)看，可以通過(guò)不同的擾流裝置改變結(jié)構(gòu)的尾流，阻止規(guī)律渦脫的產(chǎn)生，進(jìn)而實(shí)現(xiàn)渦振的抑制。對(duì)高聳塔筒結(jié)構(gòu)，目前常見(jiàn)的氣動(dòng)式減振措施僅有擾流線，未來(lái)可以考慮布置擾流板、穿孔罩等多種減振構(gòu)件。為了達(dá)到最佳的減振效果，可以多種研究手段對(duì)不同氣動(dòng)式的防治措施進(jìn)行對(duì)比，并針對(duì)各防治裝置的型式、結(jié)構(gòu)參數(shù)和材質(zhì)等參數(shù)對(duì)防治裝置減振效果的影響，給出最優(yōu)的布置方案。

除此之外，相關(guān)研究表明，通過(guò)圓柱表面布置吸氣孔改變局部流場(chǎng)的分布來(lái)進(jìn)行減振具有非常好的效果，該方式對(duì)圓柱渦振的抑振效果可達(dá)到70%。

3.3 阻尼式防治措施

阻尼減振技術(shù)目前在高層建筑、橋梁等的抗風(fēng)抗震領(lǐng)域廣泛應(yīng)用，具有可靠高效、安裝方便等優(yōu)點(diǎn)。在高聳結(jié)構(gòu)中，可以利用TMD 進(jìn)行減振。如圖3 所示，TMD 系統(tǒng)的基本減振原理為：當(dāng)塔筒渦振時(shí)，TMD 就會(huì)產(chǎn)生一個(gè)與渦振方向相反的慣性力作用在結(jié)構(gòu)上，使塔筒的振動(dòng)響應(yīng)衰減。TMD 有它自身的振動(dòng)頻率和阻尼，通過(guò)改變質(zhì)量或剛度調(diào)整TMD 的自振頻率，使其接近主結(jié)構(gòu)的基本頻率。其基本減振原理為當(dāng)塔筒振動(dòng)時(shí)，TMD 就會(huì)產(chǎn)生一個(gè)與振動(dòng)方向相反的慣性力作用在結(jié)構(gòu)上，使塔筒的振動(dòng)響應(yīng)衰減。相關(guān)研究表明，當(dāng)TMD 系統(tǒng)的自振頻率與塔筒的渦振頻率接近時(shí)，可以達(dá)到最佳減振效果。

如圖3 所示，TMD 系統(tǒng)的基本減振原理如下：當(dāng)塔筒渦振時(shí)，TMD 就會(huì)產(chǎn)生一個(gè)與渦振方向相反的慣性力作用在結(jié)構(gòu)上，使塔筒的振動(dòng)響應(yīng)衰減。TMD 可以安裝在塔筒頂部區(qū)域，在TMD 設(shè)計(jì)之前，要根據(jù)減振要求確定需要減振的響應(yīng)量以及塔筒結(jié)構(gòu)的主要模態(tài)階數(shù)。如果塔筒的渦振響應(yīng)由第一階模態(tài)控制，那么TMD 的優(yōu)化參數(shù)應(yīng)該按照第一階模態(tài)確定，靠近塔筒的頂部位置可根據(jù)TMD 優(yōu)化理論進(jìn)行計(jì)算，確定最優(yōu)設(shè)置時(shí)TMD 的頻率比、阻尼比、質(zhì)量比參數(shù)情況。

圖3 塔筒加裝TMD 的力學(xué)模型

4 結(jié)語(yǔ)

隨著陸上風(fēng)電的建設(shè)和“3060 目標(biāo)”的不斷深化，我國(guó)的風(fēng)電行業(yè)目前正迎來(lái)快速發(fā)展的階段。塔架作為風(fēng)力發(fā)電機(jī)組的重要支撐結(jié)構(gòu)，是保障機(jī)組安全運(yùn)行的核心部件，塔架的安全性和可靠性對(duì)機(jī)組的持續(xù)發(fā)電至關(guān)重要。在陸上風(fēng)電機(jī)組大型化的發(fā)展趨勢(shì)下，傳統(tǒng)柔塔所存在的渦激振動(dòng)問(wèn)題嚴(yán)重威脅風(fēng)電機(jī)組的穩(wěn)定運(yùn)行，導(dǎo)致其在向高空化發(fā)展中遇到一定技術(shù)瓶頸。在塔架高度和葉輪直徑的雙重增加下，渦振愈發(fā)頻繁，亟須解決。為解決上述技術(shù)難題，該文通過(guò)深入研究渦振的機(jī)理，提出了多種塔筒渦振響應(yīng)的計(jì)算方法，并從塔筒結(jié)構(gòu)形式、氣動(dòng)力干擾、阻尼耗能的角度出發(fā)，分別給出了渦振減振原理以及減振設(shè)計(jì)思路，對(duì)解決超高聳風(fēng)電機(jī)組塔架的渦振問(wèn)題具有一定的指導(dǎo)意義。