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

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

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

0 引 言

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

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

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

1 工程背景

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

圖1 海上風(fēng)電場(chǎ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 塔筒安裝簡(jiǎn)介

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

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

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

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

3 塔架所受彎矩和剪力計(jì)算

根據(jù)日本土木協(xié)會(huì)規(guī)范《風(fēng)力發(fā)電設(shè)備塔架結(jié)構(gòu)設(shè)計(jì)指南及解說(shuō)》[10]，對(duì)正常發(fā)電下風(fēng)機(jī)塔架風(fēng)荷載進(jìn)行計(jì)算。海上風(fēng)電機(jī)組塔架主要受到彎矩和剪力作用，在空間中，均可分解為3個(gè)方向的分量。塔架受荷具體方向如圖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個(gè)，風(fēng)輪直徑168 m，輪轂中心高度98.5 m。

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

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

(1)

(2)

(3)

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

塔架各部風(fēng)速以及荷載，計(jì)算結(jié)果如圖4所示。由圖4可知，對(duì)于海上風(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個(gè)螺栓相連接。其中，法蘭材料屬性為理想彈性，螺栓材料屬性為理想彈塑性。塔筒法蘭與連接螺栓的屈服強(qiáng)度分別為345 MPa和1 040 MPa，其彈性模量E分別取211 GPa和206 GPa，泊松比按0.3取值。各部件之間的接觸設(shè)置為法向無(wú)摩擦，切向“罰”函數(shù)的接觸，摩擦系數(shù)為0.12[11-12]。

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

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

4.2 結(jié)果與討論

以塔筒中1段與中2段連接處螺栓為例，通過(guò)圖4選取荷載數(shù)值，對(duì)數(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)力曲線(xiàn)。由圖7可知，在風(fēng)荷載作用下，螺栓與墊圈接觸的位置出現(xiàn)應(yīng)力峰值，即螺栓的上下根部位置受力較大，且下根部受力要大于上根部。本文算例中，螺栓下根部應(yīng)力為577 MPa，上根部應(yīng)力為551 MPa，其余位置維持在540 MPa左右。因此，螺栓受到風(fēng)荷載時(shí)，危險(xiǎn)位置為螺栓的上下根部，對(duì)螺栓進(jìn)行設(shè)計(jì)研究以及應(yīng)力分析時(shí)，應(yīng)著重分析螺栓的根部位置。

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

5 結(jié) 論

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

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

(2)連接塔筒法蘭的螺栓在受到風(fēng)荷載時(shí)，迎風(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)，螺栓下根部受力最大，是螺栓最危險(xiǎn)的位置。