海上風(fēng)電單樁塔筒螺栓連接法蘭數(shù)值模擬研究

孫 智，李 令，曾 佳，庫 猛，李宗澤，3，金 磊，高 鋒，

(1.中國電建集團(tuán)華東勘測設(shè)計研究院有限公司，浙江 杭州 311122；2.浙江大學(xué)海洋學(xué)院，浙江 舟山 316021；3.浙江華東工程咨詢有限公司，浙江 杭州 311122)

0 引 言

我國沿海地區(qū)海上風(fēng)能資源豐富，海上風(fēng)電場建設(shè)已成為我國風(fēng)電開發(fā)的重要方式[1]。在海上風(fēng)電場中，塔架作為發(fā)電機(jī)組的支撐結(jié)構(gòu)，其安全性是整個發(fā)電系統(tǒng)正常運(yùn)行的基本保證。當(dāng)前，大型海上風(fēng)電場風(fēng)機(jī)塔架長度長、直徑大，一般由多段塔筒組成，通過法蘭-螺栓連接[2-3]。風(fēng)電場塔架倒塌事件時有發(fā)生，倒塔的原因除了極端強(qiáng)風(fēng)天氣造成的塔筒受荷過大以外，還與風(fēng)電設(shè)備安裝連接的高強(qiáng)螺栓強(qiáng)度不足、螺栓緊固施工安裝不當(dāng)?shù)扔嘘P(guān)[4]。螺栓對塔架安全至關(guān)重要，因此分析螺栓的受力特征及應(yīng)力規(guī)律具有重要的研究價值和工程應(yīng)用價值。

目前已有學(xué)者對風(fēng)機(jī)塔架的法蘭螺栓連接問題進(jìn)行研究。龔國偉[5]對大型風(fēng)力發(fā)電機(jī)組法蘭螺栓連接的設(shè)計計算方法進(jìn)行了研究；吉亮等[6]建立螺栓軸力計算模型，并對包含翹曲缺陷法蘭螺栓的疲勞特性開展研究；徐祉康等[7]基于Petersen方法對風(fēng)力發(fā)電機(jī)組塔筒M64大直徑螺栓法蘭進(jìn)行拉伸試驗(yàn)，并對其承載能力進(jìn)行了分析研究。塔筒在風(fēng)荷載作用下，法蘭螺栓受力特征的研究較少。

本文以大連某海上風(fēng)電場為工程背景，對塔筒所受風(fēng)荷載進(jìn)行理論計算，建立螺栓-法蘭數(shù)值模型，根據(jù)理論計算結(jié)果對模型施加荷載，分析螺栓的應(yīng)力特征，研究成果可為風(fēng)電場塔筒連接設(shè)計提供參考。

1 工程背景

海上風(fēng)電項目位于大連市莊河海域，場址中心離岸距離約35.2 km，規(guī)劃面積約55.8 km2(見圖1)，裝機(jī)規(guī)模為350 MW。

圖1 海上風(fēng)電場位置示意

工程所在海域內(nèi)地基土表層以粉砂為主，區(qū)域內(nèi)特殊性巖土主要為軟土，為淤泥質(zhì)粉質(zhì)粘土，具有天然含水量高，呈流變性，承載力低，壓縮性高，抗剪強(qiáng)度低，滲透性低等特點(diǎn)。該處風(fēng)速分布較為均勻，年平均風(fēng)速為4.2 m/s，濱海年平均風(fēng)速為4.7 m/s。

2 塔筒安裝簡介

該工程塔架由4段塔筒安裝組成，分別為上、中1、中2以及下塔筒，塔架組成示意如圖2所示。塔架安裝順序?yàn)椋喊惭b前準(zhǔn)備、底段塔筒安裝、中1及中2段塔筒安裝、上段塔筒安裝[8-9]。

圖2 海上風(fēng)電塔筒組成示意

塔架安裝準(zhǔn)備工作包括基礎(chǔ)頂部法蘭面清理，塔筒位置標(biāo)記，電氣設(shè)備以及安裝工具檢查，螺栓螺母以及墊圈的檢查等工作。

安裝底段塔筒時，工程上采用垂直起吊的方式。將塔筒吊具吊帶一端掛在主吊鉤上，另一端緊固在底塔上法蘭。緩緩起鉤至吊帶即將受力時停止起鉤，將底段塔筒與設(shè)備構(gòu)架在工裝上的固定裝置卸除。再次緩緩起鉤至底段塔筒離地100 mm靜置觀察1～2 min，同時在底段塔筒的下法蘭系上攬風(fēng)繩。塔筒抬起到基礎(chǔ)上部，快要下落到既定位置時，找準(zhǔn)兩個法蘭已標(biāo)記的位置，先使用3個螺栓進(jìn)行位置固定，隨后使塔筒緩慢下降。法蘭對齊，安裝剩余螺栓、螺母和墊圈。

3 塔架所受彎矩和剪力計算

根據(jù)日本土木協(xié)會規(guī)范《風(fēng)力發(fā)電設(shè)備塔架結(jié)構(gòu)設(shè)計指南及解說》[10]，對正常發(fā)電下風(fēng)機(jī)塔架風(fēng)荷載進(jìn)行計算。海上風(fēng)電機(jī)組塔架主要受到彎矩和剪力作用，在空間中，均可分解為3個方向的分量。塔架受荷具體方向如圖3所示。圖中，MXF、MYF、MZF分別為x、y、z方向受到的彎矩，kN·m；FXF、FYF、FZF分別為x、y、z方向受到的剪力，kN。

圖3 塔架受荷方向示意

該工程風(fēng)電機(jī)組各部件尺寸及質(zhì)量如表1所示。其中，風(fēng)電塔筒高96 m，風(fēng)電機(jī)組裝機(jī)容量6 MW，葉片數(shù)3個，風(fēng)輪直徑168 m，輪轂中心高度98.5 m。

表1 風(fēng)電機(jī)組各部件尺寸及質(zhì)量

風(fēng)機(jī)塔架的直徑往往隨塔架高度的不同而變化。為了更精確地得到風(fēng)機(jī)塔架不同高度處的剪力和彎矩大小的理論解，對于96 m高的塔架劃分出30個等距的小段，其中每3.2 m為一段，共計31個截面，計算每一個高度處截面的受力情況。平均剪力QDi和平均彎矩MDi的計算公式為

(1)

(2)

(3)

式中，QDi為對應(yīng)風(fēng)速UHi正常發(fā)電時的平均剪力，N，風(fēng)速為4 m/s時對應(yīng)風(fēng)速級別為1級；ρ為空氣密度，取ρ=1.293 kg/m3；CTi為對應(yīng)風(fēng)速UHi正常發(fā)電時的風(fēng)輪推力系數(shù)，m/s；R為風(fēng)力機(jī)的風(fēng)輪半徑，m；CDN為機(jī)艙的平均阻力系數(shù)；CDT為塔架的平均阻力系數(shù)；AN為機(jī)艙的可見面積，m2；Ht為相應(yīng)截面處風(fēng)機(jī)塔架高度，m；h為截面處的評估高度，m；dz為計算截面高度z處的塔架直徑，m；Uiz為對應(yīng)風(fēng)速級別UHi時截面高度z處的風(fēng)速，m/s；Hh為輪轂高度，m；α為風(fēng)速垂直方向的分布指數(shù)，與場地地面粗糙度有關(guān)，本文取0.1。根據(jù)式(1)～(3)，可以計算出在風(fēng)速為4 m/s時，塔架各處所受平均剪力和彎矩數(shù)值。

塔架各部風(fēng)速以及荷載，計算結(jié)果如圖4所示。由圖4可知，對于海上風(fēng)力發(fā)電塔架，隨著塔架高度的降低風(fēng)速隨之減小，所受剪切力相應(yīng)增大但變化并不明顯，而塔架所受的彎矩力隨著高度的減小逐漸增大且變化非常顯著，例如本例中塔架在96 m高度處受到的彎矩為404.061 kN·m，在13.88 m高度處所受彎矩達(dá)到13 787.763 kN·m。

圖4 塔架各處受力示意

4 螺栓數(shù)值分析

4.1 塔筒法蘭-螺栓數(shù)值模型

基于ABAQUS軟件建立塔筒法蘭-螺栓數(shù)值模型，如圖5所示。上下法蘭直徑6 m，由均勻分布的80個螺栓相連接。其中，法蘭材料屬性為理想彈性，螺栓材料屬性為理想彈塑性。塔筒法蘭與連接螺栓的屈服強(qiáng)度分別為345 MPa和1 040 MPa，其彈性模量E分別取211 GPa和206 GPa，泊松比按0.3取值。各部件之間的接觸設(shè)置為法向無摩擦，切向“罰”函數(shù)的接觸，摩擦系數(shù)為0.12[11-12]。

圖5 塔筒法蘭-螺栓數(shù)值模型

根據(jù)塔架受荷計算結(jié)果(見圖4)，對數(shù)值模型施加荷載，分析在塔架在正常發(fā)電時風(fēng)速為4 m/s的工況下，螺栓的受力特征及規(guī)律。在施加風(fēng)荷載前，對每個螺栓施加預(yù)緊力，預(yù)緊力值為550 kN，待螺栓上預(yù)緊力施加完畢，將螺栓變形后長度鎖定，隨后再施加風(fēng)荷載[13]。

4.2 結(jié)果與討論

以塔筒中1段與中2段連接處螺栓為例，通過圖4選取荷載數(shù)值，對數(shù)值模型施加荷載，分析螺栓受力特征。

圖6為塔筒正面迎風(fēng)螺栓在正常風(fēng)荷載工況下的應(yīng)力云圖。由圖6可知，在風(fēng)荷載作用下，塔筒迎風(fēng)側(cè)應(yīng)力大于背風(fēng)側(cè)，背風(fēng)側(cè)螺栓應(yīng)力值僅有約350 MPa，而迎風(fēng)側(cè)螺栓應(yīng)力值最高可達(dá)577 MPa，比背風(fēng)側(cè)高64.9%。

圖6 螺栓應(yīng)力分布云圖

圖7為螺栓迎風(fēng)側(cè)應(yīng)力曲線。由圖7可知，在風(fēng)荷載作用下，螺栓與墊圈接觸的位置出現(xiàn)應(yīng)力峰值，即螺栓的上下根部位置受力較大，且下根部受力要大于上根部。本文算例中，螺栓下根部應(yīng)力為577 MPa，上根部應(yīng)力為551 MPa，其余位置維持在540 MPa左右。因此，螺栓受到風(fēng)荷載時，危險位置為螺栓的上下根部，對螺栓進(jìn)行設(shè)計研究以及應(yīng)力分析時，應(yīng)著重分析螺栓的根部位置。

圖7 螺栓迎風(fēng)側(cè)應(yīng)力曲線

5 結(jié) 論

本文以大連市某海上風(fēng)電項目為工程背景，根據(jù)相關(guān)規(guī)范對海上風(fēng)力發(fā)電塔架在正常工況下所受風(fēng)荷載進(jìn)行理論計算，建立了塔筒法蘭-螺栓的數(shù)值模型，研究了螺栓的應(yīng)力特征及規(guī)律，得到以下結(jié)論：

(1)海上風(fēng)力發(fā)電塔架所受剪切力以及彎矩隨著塔架高度的降低而增大，剪切力變化并不明顯，而彎矩力的增加非常顯著。

(2)連接塔筒法蘭的螺栓在受到風(fēng)荷載時，迎風(fēng)側(cè)的應(yīng)力大于背風(fēng)側(cè)，本文算例中迎風(fēng)側(cè)應(yīng)力達(dá)到577 MPa，比背風(fēng)側(cè)的應(yīng)力高64.9%。

(3)螺栓迎風(fēng)側(cè)的應(yīng)力峰值出現(xiàn)在根部位置，其余位置應(yīng)力較小且平穩(wěn)，螺栓下根部受力最大，是螺栓最危險的位置。