何海建，楊 揚，孟令銳，晁貫良，董姝言

(許昌許繼風(fēng)電科技有限公司，河南 許昌 461000)

0 引 言

由于國家政策的大力支持，近幾年風(fēng)力發(fā)電機組的國產(chǎn)化程度逐漸提高[1]。為了保證風(fēng)電機組設(shè)備的可靠性，需要對風(fēng)機中各部件的強度等性能進行分析。

塔筒作為風(fēng)力發(fā)電機的一個重要的組成部分，它的結(jié)構(gòu)穩(wěn)定性直接影響著風(fēng)力發(fā)電機的運行安全[2]。塔筒各段之間通過法蘭進行連接，風(fēng)力發(fā)電機組運行時，作用于葉輪和機艙的載荷通過主機架傳遞到塔筒上，同時由于風(fēng)速和風(fēng)向的不穩(wěn)定性，塔筒法蘭還要承受動態(tài)載荷，在長期載荷作用下，塔筒法蘭極易遭到破壞，如果連接法蘭出現(xiàn)問題，不僅會對風(fēng)力發(fā)電的性能造成影響，甚至?xí)l(fā)安全事故，因此為了保證風(fēng)力發(fā)電機組的正常運行，需要對風(fēng)機塔筒法蘭的極限及疲勞強度進行校核[3]。塔筒法蘭的破壞方式主要為極限載荷工況下的局部塑性變形或破壞以及隨機載荷作用下的疲勞斷裂[4]。

風(fēng)機塔筒法蘭強度計算時最常用的是有限元分析法，晁貫良[5]采用ANSYS軟件建立了塔頂法蘭連接系統(tǒng)的有限元模型，利用子模型分析技術(shù)及工況疊加計算方法對塔頂法蘭進行了極限強度校核，但計算模型中忽略了塔頂連接螺栓對法蘭極限強度的影響；杜靜[6]利用MSC.Marc/Mentat軟件建立了風(fēng)力機法蘭有限元分析模型，運用雨流計數(shù)法與Palmgren-Miner線性累積損傷理論相結(jié)合的疲勞計算方法，對塔頂法蘭疲勞強度進行了研究，但未考慮不同方向上的剪應(yīng)力對塔筒法蘭疲勞損傷產(chǎn)生的影響。

本研究以某大型風(fēng)力發(fā)電機為例，建立塔筒法蘭連接系統(tǒng)的有限元模型，計算中將考慮法蘭連接螺栓對法蘭強度的影響，同時提出將臨界平面算法與剪應(yīng)力算法相結(jié)合，進行塔筒法蘭疲勞損傷計算。

1 塔筒法蘭的幾何模型

塔筒法蘭幾何模型如圖1所示。

圖1 塔筒法蘭幾何模型

本研究采用Solidworks軟件進行塔筒法蘭連接系統(tǒng)創(chuàng)建，為了考慮連接螺栓對塔筒法蘭計算的影響，法蘭幾何模型中保留了連接螺栓孔。同時為了提高計算效率，模型中簡化了對計算結(jié)果影響較小的塔筒平臺、爬梯等附屬設(shè)備。

2 塔筒法蘭有限元模型及載荷坐標系

為了準確模擬塔筒法蘭的真實受力狀態(tài)，本研究在有限元模型中創(chuàng)建了與塔筒法蘭相鄰的部件，主要包括上段塔筒、連接螺栓、墊圈以及下段塔筒。

本研究將風(fēng)機塔筒法蘭的幾何裝配模型導(dǎo)入有限元分析軟件ANSYS中進行網(wǎng)格劃分，有限元模型剖視圖如圖2所示。

圖2 塔筒法蘭連接有限元模型

上段塔筒、連接法蘭、墊圈以及下段塔筒均采用高階六面體單元進行網(wǎng)絡(luò)劃分，為了減小結(jié)構(gòu)的關(guān)鍵部位對計算結(jié)果的影響，本研究對上、下法蘭圓角位置的網(wǎng)格進行了適度加密[7]，同時為了提高塔筒法蘭疲勞強度計算效率以及采用臨界平面方法進行疲勞計算，在上、下法蘭的圓角位置創(chuàng)建了一層厚度較薄的殼單元，殼單元與下覆實體單元之間通過共節(jié)點的方式進行連接。此外，連接螺栓采用梁單元進行模擬，單元截面積與螺栓應(yīng)力面積相同。上、下法蘭之間通過標準摩擦接觸方式進行連接，摩擦系數(shù)為0.2。塔筒法蘭有限元模型共包含167 040個單元，722 305個節(jié)點。

為了準確模擬載荷的傳遞，本研究在塔筒上、下法蘭接合面中心位置建立節(jié)點并施加外載荷，并將該中心節(jié)點通過MPC剛性梁單元與上段塔筒頂部位置進行連接，以便模擬載荷的傳遞，根據(jù)塔筒實際的受力狀態(tài)，需要對下段塔筒底部端面所有節(jié)點的平動自由度進行約束。

塔筒法蘭各部件所用材料均為低合金高強度結(jié)構(gòu)鋼，其彈性模量為2.1×1011Pa，泊松比為0.3，密度為7 850 kg/m3。塔筒法蘭由Q345鍛造而成，其最大壁厚為150 mm，根據(jù)《GB/T1591-2008低合金高強度結(jié)構(gòu)鋼》塔筒法蘭抗拉強度為470 MPa，屈服強度為285 MPa。

風(fēng)力發(fā)電機組運行時，葉片產(chǎn)生的氣動載荷以及由于風(fēng)輪旋轉(zhuǎn)和機艙對風(fēng)引起的離心力、慣性力通過主機架傳遞到塔筒上，這些載荷和塔筒自身的重力構(gòu)成了塔筒法蘭載荷[8]。塔筒法蘭所承受的載荷通過GH Bladed軟件根據(jù)GL規(guī)范計算仿真得到，載荷方向與GL規(guī)定的塔筒載荷坐標系方向一致，坐標系如圖3所示。坐標系原點位于塔筒中心位置，XF沿水平方向并指向下風(fēng)向，ZF豎直向上，YF位于水平方向并與XF和ZF組成右手坐標系。

圖3 塔筒載荷坐標系

3 塔筒法蘭靜強度分析

以某MW級風(fēng)力發(fā)電機組為例，其塔筒法蘭中心位置的極限工況載荷如表1所示。

塔筒法蘭靜強度計算分兩個載荷步進行，第一個載荷步，對連接螺栓組施加預(yù)緊力載荷；第二個載荷步，鎖定第一步施加的螺栓預(yù)緊力，并施加表1中的極限工況載荷[9]。由于上、下法蘭之間采用標準摩擦接觸方式進行連接，有限元模型采用隱式非線性牛頓-拉斐遜迭代算法進行求解。根據(jù)GL規(guī)范，為了保證載荷與材料的安全設(shè)計值，塔筒法蘭靜強度安全系數(shù)應(yīng)通過下式進行計算：

表1 塔筒法蘭極限工況載荷表

(1)

式中：γf—載荷安全系數(shù)，取1.0；γm—材料安全系數(shù)，取1.1；σeqv，max—塔筒法蘭最大等效應(yīng)力，MPa；Re—塔筒法蘭屈服強度，MPa。

塔筒法蘭危險工況下等效應(yīng)力云圖如圖4所示。

圖4 塔筒法蘭等效應(yīng)力云圖

從該云圖可以看出：塔筒法蘭靜強度危險位置發(fā)生在圓角過渡處，最大應(yīng)力值為235.4 MPa，對應(yīng)靜強度安全系數(shù)為1.1，大于1，因此，塔筒法蘭在極限工況下運行能夠滿足靜強度設(shè)計要求。

4 塔筒法蘭疲勞強度分析

工程實踐表明：部件的初始裂紋通常首先出現(xiàn)在應(yīng)力較大位置，對于塔筒法蘭，倒圓角區(qū)域為其應(yīng)力集中位置，為了提高計算效率，本研究僅計算塔筒法蘭圓角位置的疲勞強度。塔筒法蘭疲勞強度計算流程如圖5所示。

由于塔筒法蘭疲勞強度主要受到My載荷分量的影響，本研究首先在塔筒法蘭接合面中心位置施加正、負兩個方向上的疲勞極限載荷(+Mymax、-Mymin)，并在這兩個疲勞極限工況下，對風(fēng)機塔筒法蘭的有限元模型進行非線性求解。為了考慮法蘭連接螺栓的預(yù)緊力效應(yīng)，本研究對這兩個疲勞極限工況進行非線性求解時，分兩個載荷步進行計算，第一個載荷步施加連接螺栓預(yù)緊力載荷，第二個載荷步施加疲勞極限外載荷，施加疲勞極限外載荷時分為5個子步進行，根據(jù)各載荷步的計算結(jié)果，能夠得到當疲勞載荷由負向極值變化到正向極值時，塔筒法蘭倒圓位置載荷與各臨界平面剪應(yīng)力的關(guān)系曲線。以法蘭某個熱點位置為例，其載荷與危險臨界平面剪應(yīng)力關(guān)系曲線如圖6所示。

圖5 塔筒法蘭疲勞強度計算流程

圖6 法蘭某熱點位置載荷與危險臨界平面剪應(yīng)力關(guān)系曲線

其次，根據(jù)載荷與臨界平面剪應(yīng)力關(guān)系曲線以及疲勞時序載荷譜，通過GH-Bladed軟件可計算得到法蘭圓角位置各臨界平面下的剪應(yīng)力譜。最后，本研究將法蘭各臨界平面下的剪應(yīng)力譜、應(yīng)力譜循環(huán)次數(shù)以及法蘭S-N曲線數(shù)據(jù)輸入至疲勞分析軟件，基于Miner線性累積損傷理論[10]，采用Goodman平均應(yīng)力修正方法計算塔筒法蘭疲勞損傷及安全系數(shù)。

根據(jù)GL規(guī)范，塔筒法蘭作為鍛造件，其S-N曲線可由相關(guān)材料參數(shù)擬合得到，用于擬合塔筒法蘭S-N曲線的輸入?yún)?shù)如表2所示。

表2 塔筒法蘭相關(guān)輸入?yún)?shù)

由于在疲勞分析軟件中使用了Goodman平均應(yīng)力修正方法，表2中平均應(yīng)力修正系數(shù)Fm=1，即S-N曲線中未考慮平均應(yīng)力的影響。根據(jù)表2中的輸入?yún)?shù)，可得到塔筒法蘭S-N曲線如圖7所示。

圖7 塔筒法蘭S-N曲線

計算結(jié)果顯示：在疲勞載荷作用下，塔筒法蘭圓角位置的疲勞強度安全系數(shù)為5.163，因此，在20年的設(shè)計壽命周期內(nèi)，該機組塔筒法蘭能夠滿足疲勞強度設(shè)計要求。

5 結(jié)束語

本研究介紹了MW級風(fēng)力發(fā)電機組塔筒法蘭極限強度及疲勞強度計算方法。提出了進行塔筒法蘭疲勞計算時采用臨界平面算法與剪應(yīng)力算法相結(jié)合的疲勞損傷計算方法，該方法能夠更為準確地考慮到塔筒法蘭的真實疲勞應(yīng)力歷程。

計算結(jié)果表明：塔筒法蘭的極限強度及疲勞強度安全系數(shù)均大于1，能夠滿足強度設(shè)計要求，同時塔筒法蘭的危險位置與工程實際較為吻合，說明提出的計算方法是可靠實用的，能夠為風(fēng)力發(fā)電機組塔筒法蘭的設(shè)計提供必要的依據(jù)，具有一定的工程實用價值。

[1] 周姝燦，唐俊熙.風(fēng)電、儲能參與系統(tǒng)調(diào)頻的協(xié)調(diào)控制策略研究[J].機電工程技術(shù)，2017，46(11)：31-35.

[2] 關(guān)媛媛，劉大江，邢 艷.L型法蘭盤在風(fēng)力發(fā)電機上的應(yīng)用[J].硅谷，2014(6):77-78.

[3] 張 榮，何玉林，牛興海，等.MW級風(fēng)力發(fā)電機塔頂法蘭連接系統(tǒng)有限元分析[J].制造業(yè)自動化，2012，34(23):112-114.

[4] Germanischer Lloyd. Rules and guidelines industrial services: guideline for the certification of wind turbines[M]. Hamburg: Germanischer Lloyd，2010.

[5] 晁貫良，蘇鳳宇，董姝言，等.MW級風(fēng)力發(fā)電機塔頂法蘭有限元計算[J].機械與電子，2016，34(1):46-50

[6] 杜 靜，牛興海，何玉林，等.大功率風(fēng)力發(fā)電機組塔頂法蘭疲勞分析[J].機械科學(xué)與技術(shù)，2012，31(10):53-57

[7] 杜 靜，?；塾?，石秉楠.風(fēng)力發(fā)電機塔筒頂部法蘭的有限元分析[J].現(xiàn)代制造工程，2014(4):54-58.

[8] 姚興佳，劉光德，邢作霞，等.大型變速風(fēng)力發(fā)電機組總體設(shè)計中的幾個問題探討[J].沈陽工業(yè)大學(xué)學(xué)報，2006，28(2):5-11.

[9] 余 娟，余明友，金建偉，等.螺栓連接中預(yù)緊力的有限元分析[J].軟件導(dǎo)刊，2013(9):64-66.

[10] 濮良貴，紀名剛.機械設(shè)計[M].北京：高等教育出版社，2001.