鋼管混凝土格構(gòu)式新型風(fēng)力發(fā)電機(jī)塔架風(fēng)振響應(yīng)

李 斌，文昊天，高春彥

（內(nèi)蒙古科技大學(xué) 土木工程學(xué)院土木工程系，內(nèi)蒙古 包頭 014010）

風(fēng)力發(fā)電單機(jī)容量的不斷增大使得應(yīng)用于大型風(fēng)電機(jī)組的傳統(tǒng)塔架結(jié)構(gòu)安全性越來(lái)越低[1]，風(fēng)電事故時(shí)有發(fā)生。當(dāng)風(fēng)力發(fā)電機(jī)塔筒超過(guò)一定高度后，其結(jié)構(gòu)的力學(xué)性能下降非常明顯，主要表現(xiàn)在：塔頂在水平荷載作用下變形過(guò)大，塔身與上部風(fēng)機(jī)產(chǎn)生共振，塔身容易發(fā)生疲勞破壞[2]。塔架是風(fēng)力發(fā)電機(jī)組的主要承載結(jié)構(gòu)，不但要承受機(jī)艙重量、風(fēng)輪作用力以及風(fēng)對(duì)塔架的彎、剪、扭等靜力作用，還要承受振動(dòng)作用力[3]。而且鋼制高聳結(jié)構(gòu)對(duì)局部缺陷十分敏感，其臨界承載力只能達(dá)到理論值的1/5～1/3，很容易發(fā)生失穩(wěn)破壞[4]。塔架是為風(fēng)力發(fā)電機(jī)組正常工作服務(wù)的結(jié)構(gòu)，其發(fā)展水平應(yīng)與風(fēng)力發(fā)電機(jī)組的發(fā)展相協(xié)調(diào)。為了滿足現(xiàn)代風(fēng)力發(fā)電機(jī)的工作需求，必須要保證塔架在風(fēng)荷載作用下的安全。

根據(jù)各國(guó)風(fēng)機(jī)受損的歷史統(tǒng)計(jì)數(shù)據(jù)，塔架是風(fēng)機(jī)受損部件中遭受破壞比率較大的部件之一，占到19%[5]。從功能性角度來(lái)看，塔架是風(fēng)電機(jī)組中重要的支撐結(jié)構(gòu)，為風(fēng)力發(fā)電機(jī)提供必要的工作高度；從荷載角度來(lái)看，塔架的重量占到風(fēng)機(jī)總重的1/2左右；從經(jīng)濟(jì)性角度來(lái)看，塔架的成本大約為整個(gè)風(fēng)機(jī)制造總成本的15%左右，而且塔架倒塌必定造成風(fēng)力發(fā)電機(jī)的損傷[6]。因此，鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架在風(fēng)荷載作用下動(dòng)力響應(yīng)的研究對(duì)于提高塔架的設(shè)計(jì)水準(zhǔn)、保證風(fēng)機(jī)的運(yùn)行安全和整機(jī)性能都有重要意義。

1 風(fēng)速時(shí)程模擬

對(duì)本文所采用的鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架不同高度處的13個(gè)計(jì)算點(diǎn)進(jìn)行風(fēng)速時(shí)程模擬，計(jì)算點(diǎn)位置選取如圖1所示。

圖1 塔架橫截面與計(jì)算點(diǎn)位置示意圖

其中13號(hào)節(jié)點(diǎn)為風(fēng)輪中心位置。選取我國(guó)規(guī)范所采用的Davenport脈動(dòng)水平風(fēng)速譜以及Shiotami的簡(jiǎn)化空間相關(guān)性系數(shù)表達(dá)式。為了得到較為準(zhǔn)確的脈動(dòng)風(fēng)速時(shí)程曲線，采用諧波合成法對(duì)塔架工作地區(qū)的脈動(dòng)風(fēng)進(jìn)行風(fēng)速時(shí)程模擬；為了得到更準(zhǔn)確的平均風(fēng)速，以對(duì)數(shù)律求解平均風(fēng)速廓線。具體模擬參數(shù)如表1所示。

表1 風(fēng)速時(shí)程模擬參數(shù)

借助MATLAB軟件編寫(xiě)程序，最終生成各計(jì)算點(diǎn)的風(fēng)速時(shí)程曲線。

2 風(fēng)力發(fā)電機(jī)塔架有限元模型建立

本文所采用的2 MW風(fēng)力發(fā)電機(jī)塔架高度為60 m，風(fēng)機(jī)轉(zhuǎn)子中心高度為61.5 m。塔架分12層，層高有4 m、5 m和6 m三種。應(yīng)用ABAQUS軟件建立鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架有限元模型，鋼管模型采用殼單元，混凝土模型采用實(shí)體單元。將鋼管與混凝土的接觸定義為剛體約束，忽略鋼管與混凝土在受力過(guò)程中的局部滑移。將所有焊接節(jié)點(diǎn)設(shè)置為綁定約束，另外忽略風(fēng)力發(fā)電機(jī)塔架的基礎(chǔ)與地基土之間復(fù)雜的相互作用，將之視為剛性連接。為了提高網(wǎng)格質(zhì)量，對(duì)裝配件整體布種，對(duì)不同單元分別指派網(wǎng)格控制屬性。在分析步中設(shè)置頻率分析步，提取塔架的各階頻率和振型，用于后續(xù)風(fēng)振響應(yīng)計(jì)算。

圖2 風(fēng)力發(fā)電機(jī)塔架有限元模型

3 風(fēng)力發(fā)電機(jī)塔架風(fēng)致響應(yīng)分析

由模擬得到的各個(gè)計(jì)算點(diǎn)風(fēng)速時(shí)程曲線求解對(duì)應(yīng)的風(fēng)荷載時(shí)程曲線，利用ABAQUS有限元軟件將風(fēng)荷載以節(jié)點(diǎn)力的形式作用于塔架模型對(duì)應(yīng)的計(jì)算點(diǎn)，分析得到鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架各計(jì)算點(diǎn)的位移、速度和加速度響應(yīng)曲線。

3.1 塔架風(fēng)輪中心位置處的風(fēng)振響應(yīng)

提取風(fēng)輪中心位置處的位移、速度和加速度響應(yīng)時(shí)程曲線分別如圖3、圖4、圖5所示。

圖3 13號(hào)節(jié)點(diǎn)位移響應(yīng)時(shí)程曲線

圖4 13號(hào)節(jié)點(diǎn)速度響應(yīng)時(shí)程曲線

圖5 13號(hào)節(jié)點(diǎn)加速度響應(yīng)時(shí)程曲線

從圖3可以看出，塔架風(fēng)輪中心位置處的順風(fēng)向位移在0.07 m～0.23 m之間，位移響應(yīng)峰值為0.23 m，平衡位置與初始位置的水平距離約為0.17 m，在模擬的時(shí)段內(nèi)13號(hào)節(jié)點(diǎn)位置處的位移響應(yīng)表現(xiàn)出很強(qiáng)的隨機(jī)性。

由圖4可知13號(hào)節(jié)點(diǎn)的速度響應(yīng)在±0.06 m/s之間，大部分時(shí)間處于±0.02 m/s之間，速度峰值響應(yīng)為0.06 m/s。

由圖5可知，風(fēng)輪中心位置處加速度峰值響應(yīng)為0.11 m/s2，同樣也表現(xiàn)出極強(qiáng)的隨機(jī)性和波動(dòng)性。

3.2 塔架整體風(fēng)振響應(yīng)

從表2可以看出隨著高度的增加，結(jié)構(gòu)的位移響應(yīng)均值、位移響應(yīng)極值以及加速度響應(yīng)均方差都在增加。

位移響應(yīng)極值與位移響應(yīng)均值的差距逐漸增大，脈動(dòng)風(fēng)對(duì)結(jié)構(gòu)位移響應(yīng)的作用是隨高度不斷增加的。在61.5 m高度處脈動(dòng)風(fēng)位移響應(yīng)分量已接近總風(fēng)響應(yīng)的30%。

圖6、圖7所示分別為鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架各節(jié)點(diǎn)位移的響應(yīng)均值和加速度響應(yīng)的均方差隨塔架高度變化的關(guān)系。

從圖中可以看出位移響應(yīng)均值和加速度響應(yīng)均方差隨高度的增加逐漸變大，且變化速度越來(lái)越快。這主要是因?yàn)殡S著高度的增加結(jié)構(gòu)的橫截面面積減小、剛度下降而荷載反而增大所引起的。

圖6 塔架位移響應(yīng)均值

圖7 塔架加速度響應(yīng)均方差

為了研究不同風(fēng)速工況下鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架的風(fēng)振響應(yīng)，在不改變其他參數(shù)的條件條件下，分別取vˉ10=10 m/s、20 m/s、25 m/s、30 m/s對(duì)塔架的8號(hào)節(jié)點(diǎn)和13號(hào)節(jié)點(diǎn)進(jìn)行風(fēng)振響應(yīng)計(jì)算。

圖8中就風(fēng)力發(fā)電機(jī)塔架8號(hào)節(jié)點(diǎn)和13號(hào)節(jié)點(diǎn)在不同風(fēng)速工況下的加速度響應(yīng)均方差進(jìn)行了比較。從圖中可以看出節(jié)點(diǎn)的風(fēng)振響應(yīng)隨著基本風(fēng)速的增加單調(diào)增加。節(jié)點(diǎn)加速度響應(yīng)均方差在10 m/s～20 m/s的風(fēng)速范圍內(nèi)增長(zhǎng)較為緩慢，在20 m/s～25 m/s的風(fēng)速范圍內(nèi)增速明顯增大，之后又有所下降。相比之下，高度越高，風(fēng)振響應(yīng)對(duì)風(fēng)速的變化越敏感。

圖8 不同風(fēng)速工況下節(jié)點(diǎn)加速度響應(yīng)均方差

不同結(jié)構(gòu)形式的阻尼比取值不同，阻尼比會(huì)影響動(dòng)力時(shí)程計(jì)算的結(jié)果。為了分析不同阻尼比對(duì)塔架風(fēng)振響應(yīng)的影響，取基本風(fēng)速其他參數(shù)保持不變，分別取阻尼比ζ=0.02、0.035、0.05，計(jì)算塔架8號(hào)節(jié)點(diǎn)與13號(hào)節(jié)點(diǎn)的加速度響應(yīng)時(shí)程。

由圖9可以看出，節(jié)點(diǎn)的風(fēng)振響應(yīng)隨著阻尼比的增大而減小。在鋼管中灌入混凝土可以提高結(jié)構(gòu)的阻尼比，在一定程度上達(dá)到減振的目的。

圖9 不同阻尼比時(shí)節(jié)點(diǎn)加速度響應(yīng)均方差

4 塔架順風(fēng)向風(fēng)振系數(shù)

4.1 塔架結(jié)構(gòu)分段

本文研究的鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架為高聳結(jié)構(gòu)，為保證準(zhǔn)確性，對(duì)風(fēng)振系數(shù)的分析和計(jì)算必須分段進(jìn)行。

如圖10所示，將鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架沿高度方向分為6段。

各塔段結(jié)構(gòu)參數(shù)如表2所示，其中迎風(fēng)面寬度為各塔段二分之一高度處的迎風(fēng)面寬度。

4.2 不同方法中的風(fēng)振系數(shù)計(jì)算

計(jì)算所選參數(shù)見(jiàn)表3至表7。

圖10 塔架分段圖

4.3 結(jié)果分析

將上一節(jié)中按照四種不同方法計(jì)算得到的風(fēng)振系數(shù)進(jìn)行比較，見(jiàn)圖11。

圖11 不同方法計(jì)算所得的風(fēng)振系數(shù)對(duì)比

從圖中可以看出，基于隨機(jī)振動(dòng)理論按Davenport譜和按時(shí)程分析結(jié)果所得到的風(fēng)振系數(shù)基本吻合。在按隨機(jī)振動(dòng)理論計(jì)算塔架響應(yīng)時(shí)只取了第1階振型，這表明基于隨機(jī)振動(dòng)理論計(jì)算鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架風(fēng)振系數(shù)只取第1階振型進(jìn)行計(jì)算是可行的。因?yàn)橹蝗×似鹂刂谱饔玫牡?階振型，隨機(jī)振動(dòng)理論的計(jì)算值較數(shù)值模擬結(jié)果略微偏低。

表3 鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架各塔段結(jié)構(gòu)參數(shù)

表4 基于Davenport風(fēng)譜計(jì)算所得的風(fēng)電塔架風(fēng)振系數(shù)

表5 基于《建筑結(jié)構(gòu)荷載規(guī)范》（GB50009-2012）計(jì)算所得的風(fēng)電塔架風(fēng)振系數(shù)

表6 基于《高聳結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB50135-2006）計(jì)算所得的風(fēng)電塔架風(fēng)振系數(shù)

表7 基于數(shù)值分析結(jié)果計(jì)算所得的風(fēng)振系數(shù)

按照《高聳結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB50135-2006）計(jì)算所得的風(fēng)振系數(shù)與數(shù)值模擬的結(jié)果相比較，在50 m以下尚可，超過(guò)50 m以后出現(xiàn)了偏低的情況，這對(duì)結(jié)構(gòu)的安全性是不利的。現(xiàn)行《高聳結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB50135-2006）計(jì)算參數(shù)少，計(jì)算便利，但參數(shù)的物理意義不夠明確，且針對(duì)不同結(jié)構(gòu)形式的參數(shù)選用描述模糊、不夠全面，在確定系數(shù)時(shí)往往需要自行插值計(jì)算，受主觀影響較大。

按照《建筑結(jié)構(gòu)荷載規(guī)范》（GB50009-2012）計(jì)算所得的風(fēng)振系數(shù)與數(shù)值模擬的結(jié)果吻合度最差。新版《建筑結(jié)構(gòu)荷載規(guī)范》（GB50009-2012）關(guān)于風(fēng)振系數(shù)計(jì)算的參數(shù)設(shè)置較為詳細(xì)，每個(gè)參數(shù)的物理意義清晰，但對(duì)于鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架這種結(jié)構(gòu)形式來(lái)說(shuō)計(jì)算結(jié)果過(guò)于保守，并不十分適用。

5 結(jié)語(yǔ)

（1）在基本風(fēng)速為29.94 m/s時(shí)，鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架頂點(diǎn)的脈動(dòng)風(fēng)位移響應(yīng)分量接近總位移響應(yīng)的30%，風(fēng)振響應(yīng)顯著，脈動(dòng)風(fēng)對(duì)結(jié)構(gòu)響應(yīng)影響很大。

（2）結(jié)構(gòu)的風(fēng)振響應(yīng)與風(fēng)速和結(jié)構(gòu)的阻尼比相關(guān)，它隨著風(fēng)速的增加而增加，隨著結(jié)構(gòu)阻尼比的增加而減小。鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架結(jié)構(gòu)阻尼比較純鋼結(jié)構(gòu)塔架大，具有更好的抗風(fēng)振能力。

（3）現(xiàn)行《建筑結(jié)構(gòu)荷載規(guī)范》（GB50009-2012）和《高聳結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB50135-2006）提供的風(fēng)振系數(shù)計(jì)算方法對(duì)于新型鋼管混凝土格構(gòu)式風(fēng)力發(fā)電機(jī)塔架并不適用，誤差較大。