李建華,魏克湘,毛文貴,王高升

(1.湖南工程學(xué)院 機(jī)械工程學(xué)院，湘潭 411104；2.湖南工程學(xué)院 湖南省風(fēng)電裝備與電能變換協(xié)同創(chuàng)新中心，湘潭 411104)

小型風(fēng)力發(fā)電機(jī)塔架分析及優(yōu)化設(shè)計(jì)*

李建華1,2,魏克湘1,毛文貴1,王高升1

(1.湖南工程學(xué)院 機(jī)械工程學(xué)院，湘潭 411104；2.湖南工程學(xué)院 湖南省風(fēng)電裝備與電能變換協(xié)同創(chuàng)新中心，湘潭 411104)

塔架是風(fēng)力發(fā)電機(jī)的重要受力組件，塔架的截面形狀、塔架高度和最小直徑對(duì)塔架的力學(xué)性能影響顯著.提出基于遺傳算法優(yōu)化塔架的厚度和直徑，避免因發(fā)電機(jī)的激勵(lì)引起的共振、使塔架的質(zhì)量最輕，降低塔架的制造成本.為風(fēng)力發(fā)電機(jī)塔架的設(shè)計(jì)制造提供了理論依據(jù).

小型風(fēng)力發(fā)電機(jī)；塔架；遺傳算法；優(yōu)化

0 前言

天然氣、石油等傳統(tǒng)能量的減少和環(huán)境污染嚴(yán)重，風(fēng)能作為取之不盡的清潔能源已備受關(guān)注，受到高度重視.風(fēng)力發(fā)電機(jī)將風(fēng)能轉(zhuǎn)換為電能，為了獲得良好的風(fēng)能資源，風(fēng)力發(fā)電機(jī)置于高的塔架上.塔架是風(fēng)力機(jī)的重要組件，是風(fēng)機(jī)的主要受力部件之一.現(xiàn)代風(fēng)力發(fā)電機(jī)塔架采用薄殼型筒狀高聳結(jié)構(gòu)，底端固定，頂端自由且承受機(jī)艙和風(fēng)輪的重力和氣動(dòng)載荷，容易發(fā)生振動(dòng)和失穩(wěn).隨著機(jī)組的單機(jī)容量不斷增大，塔架重量和外型尺寸隨之增大，高度也隨之增加，作用在塔架上的重力載荷和氣動(dòng)載荷交變性和時(shí)變性更加顯著.同時(shí)，漿葉的掃風(fēng)范圍大，風(fēng)載的情況變化復(fù)雜.近年來(lái)，塔筒在風(fēng)力發(fā)電機(jī)組運(yùn)行過(guò)程中倒塌或失穩(wěn)的事故時(shí)有發(fā)生[1-4].因此，塔架的力學(xué)分析對(duì)于研究風(fēng)力發(fā)電機(jī)的安全性能具有重要意義.

目前，對(duì)風(fēng)力發(fā)電機(jī)的研究大都集中在葉片、發(fā)動(dòng)機(jī)性能等方面，對(duì)于塔架的研究很少.本文根據(jù)風(fēng)力發(fā)電機(jī)的實(shí)際受力對(duì)塔架進(jìn)行設(shè)計(jì)和優(yōu)化，保證塔架的剛度和穩(wěn)定性,避免塔架的最低自然頻率不會(huì)被葉片的轉(zhuǎn)動(dòng)頻率顯著影響等前提下,減小塔架重量,進(jìn)而降低成本.為風(fēng)力發(fā)電機(jī)塔架的設(shè)計(jì)制造提供了理論依據(jù).

1 風(fēng)力發(fā)電機(jī)塔架的計(jì)算模型

在文獻(xiàn)[5]研究中，通過(guò)靜力分析得知，塔的最大應(yīng)力出現(xiàn)在根部，隨著塔架高度的增加，應(yīng)力逐漸減少，得出了“塔架的截面按照變截面設(shè)計(jì)是正確的”結(jié)論.因此，本文中風(fēng)力發(fā)電機(jī)塔架基本結(jié)構(gòu)仍采用變截面，如圖1所示.本設(shè)計(jì)主要是確定塔架截面形狀、各部分鋼材的厚度和塔架頂部和底部的直徑，在設(shè)計(jì)底座之前還應(yīng)計(jì)算水平受力和底部?jī)A覆力矩.滿足最低制造和運(yùn)輸成本的最小塔架質(zhì)量將通過(guò)優(yōu)化問(wèn)題解決.

圖1 塔架結(jié)構(gòu)示意圖

風(fēng)力發(fā)電機(jī)組在運(yùn)行時(shí)，涉及的IEC負(fù)荷包括：A風(fēng)輪靜止時(shí)的待機(jī)負(fù)荷，B正常運(yùn)行中的疲勞負(fù)荷，C偏航負(fù)荷三種形式.塔架受力狀況極為復(fù)雜，塔架負(fù)荷主要包括水平拉力和底部?jī)A覆力矩.由極端風(fēng)速導(dǎo)致的塔架承受的水平力通常大于最大推力Tmax,因此，選取塔架形狀時(shí)需要考慮減小阻力.

(1)水平拉力

塔架上垂直高度方向單位高度上受到的水平拉力：

(1)

式中：d(y)表示直徑；Cd為拉力系數(shù)；U(y)為平均風(fēng)速，由公式(2)計(jì)算；ρ為風(fēng)的密度.

(2)

式中：hr為參考高度，一般取10 m;n是一個(gè)依賴于地表粗糙度的指數(shù)項(xiàng).

文獻(xiàn)[6]給出了n和粗糙度z之間的關(guān)系：

+0.24

(3)

(2) 底部?jī)A覆力矩

(4)

對(duì)于線性錐形塔架，忽略因滑動(dòng)連接導(dǎo)致的偏移量，有：

(5)

式中：M0由IEC負(fù)荷形式確定；dh為塔架底座直徑；d0為塔頂最小直徑；h為塔架高度.

設(shè)計(jì)塔架除了要滿足基本結(jié)構(gòu)的強(qiáng)度和剛度外，還要考慮塔架偏向和自然頻率，塔架偏向即風(fēng)力機(jī)和塔架的最大偏移，這個(gè)量必須取得很小，以保證塔架結(jié)構(gòu)形狀在負(fù)荷作用下不變的假設(shè)，從而不導(dǎo)致負(fù)荷變化.自然頻率的計(jì)算是用來(lái)判斷葉片經(jīng)過(guò)的頻率與自然頻率是否吻合，以致產(chǎn)生共振.

(3)塔架的優(yōu)化

塔架的質(zhì)量與制造、運(yùn)輸和地基成本密切相關(guān)，需要在設(shè)計(jì)中最小化.本文基于遺傳算法進(jìn)化策略進(jìn)行優(yōu)化塔架底座直徑、塔頂最小直徑和各段厚度.優(yōu)化過(guò)程中將塔架質(zhì)量的最小化、最大容量因數(shù)在安全范圍內(nèi)最大化及塔架和風(fēng)力發(fā)電機(jī)自然頻率的最大化三個(gè)目標(biāo)結(jié)合起來(lái)，適應(yīng)度目標(biāo)函數(shù)如下：

(6)

式中：am，aCF代表塔架質(zhì)量和最大容量因數(shù)優(yōu)化中所占比份.mt為塔架質(zhì)量；CF為最大容量因數(shù)；f為塔架和風(fēng)力發(fā)電機(jī)自然頻率.ci代表塔架底座直徑、塔頂最小直徑和各段厚度各優(yōu)化對(duì)象.

2 塔架的設(shè)計(jì)

2.1 塔架的截面形狀

本文以小型風(fēng)力發(fā)電機(jī)為例，塔架材料和性能要求如表1所示.取直徑為290 mm，厚度為3.9 mm，高度為18 m 時(shí)對(duì)圓形截面和八邊形兩種截面形狀塔架進(jìn)行分析，如表2所示，面對(duì)同樣的環(huán)境，八邊形截面形狀塔架要比圓形截面塔架質(zhì)量輕，風(fēng)力機(jī)偏移小.

表1 塔架材料和性能要求

表2 兩種截面形狀的塔架性能

2.2 三種IEC負(fù)荷下塔架的性能分析

針對(duì)八邊形截面形狀塔架，根據(jù)運(yùn)行環(huán)境和已知風(fēng)速等參數(shù)，參考文獻(xiàn)[7]可獲得風(fēng)輪靜止時(shí)的待機(jī)負(fù)荷A類型下，水平推力為2200 N，力矩為0 N；正常運(yùn)行中的疲勞負(fù)荷B類型中平均推力為1143 N，作用為風(fēng)力發(fā)電機(jī)軸承上的軸承彎矩為526 N·m；偏航負(fù)荷C類型下水平推力為1143 N，與其對(duì)應(yīng)的回轉(zhuǎn)力矩為3097 N·m.對(duì)以上三種類型負(fù)荷工況下，分析塔架水平拉力、底部?jī)A覆力矩、風(fēng)力機(jī)偏移和自然頻率隨塔架高度和最小直徑的變化情況如圖2、圖3所示，由圖2可知塔架水平拉力、底部?jī)A覆力矩、風(fēng)力機(jī)偏移、自然頻率和最大容量因數(shù)隨著高度的增加而增加，自然頻率隨著高度的增加而減小.塔架設(shè)計(jì)要防止共振現(xiàn)象的發(fā)生.由圖3可知最小直徑對(duì)最優(yōu)塔架的質(zhì)量和最大容量因數(shù)的分布有重要影響.減小最小直徑會(huì)導(dǎo)致最大容量因數(shù)增加.同時(shí)塔頂?shù)脑O(shè)計(jì)會(huì)影響底座的彎曲力矩.實(shí)際設(shè)計(jì)中最小直徑的選取由塔架法蘭和偏航機(jī)構(gòu)的設(shè)計(jì)確定.

圖2 塔架的性能與塔架基礎(chǔ)高度關(guān)系

圖3 塔架的性能與塔架最小直徑關(guān)系

3 塔架的優(yōu)化

采用遺傳算法[8]對(duì)八邊形截面形狀塔架進(jìn)行如公式(6)所示的單目標(biāo)優(yōu)化，設(shè)計(jì)變量的范圍如下：塔架的直徑范圍為[0.17 m,0.5 m]，塔架的三段厚度范圍為[2 mm,50 mm].遺傳算法采用的參數(shù)設(shè)置如下：種群規(guī)模為2000，交叉變異概率為0.1，個(gè)體數(shù)為1000，塔架質(zhì)量和最大容量因數(shù)優(yōu)化中所占比份各為0.5，終止代數(shù)為200.通過(guò)優(yōu)化獲得優(yōu)化結(jié)果如表3所示.比優(yōu)化前最小塔架質(zhì)量小6.24%.

表3 塔架優(yōu)化結(jié)果

4 結(jié) 論

本文針對(duì)塔架的形狀、塔架基礎(chǔ)高度和最小直徑對(duì)塔架的水平受力、水平拉力、底部?jī)A覆力矩、風(fēng)力機(jī)偏移、自然頻率和最大容量因數(shù)等性能的影響.基于遺傳算法對(duì)小型風(fēng)力發(fā)電機(jī)塔架進(jìn)行優(yōu)化，通過(guò)優(yōu)化塔架的三段壁厚和直徑改善了塔架的固有頻率，減輕了塔架的質(zhì)量，為塔架的設(shè)計(jì)和制造提供理論參考.

[1] 金映麗,谷繼孟,馬禮鵬. 大型風(fēng)力發(fā)電機(jī)組塔架的靜強(qiáng)度及模態(tài)分析[J]. 機(jī)械工程與自動(dòng)化,2017,(1): 63-65.

[2] 張 羽, 蔡 新, 高 強(qiáng),等. 風(fēng)力機(jī)塔架結(jié)構(gòu)研究概述[J]. 工程設(shè)計(jì)學(xué)報(bào),2016,23(2): 108-115.

[3] 尹堯杰, 廖明夫, 呂 品. 電壓跌落對(duì)風(fēng)力機(jī)塔架振動(dòng)影響分析[J]. 機(jī)械科學(xué)與技術(shù),2016, 35(4): 501-506.

[4] 郭 鵬, 徐 明, 白 楠, 等. 基于SCADA運(yùn)行數(shù)據(jù)的風(fēng)電機(jī)組塔架振動(dòng)建模與監(jiān)測(cè)[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2013,33(5):128-135.

[5] 趙吉文，劉永斌，陳軍寧，等. 基于混沌搜索的600KW風(fēng)力機(jī)塔架結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J]. 系統(tǒng)仿真學(xué)報(bào),2007, 19(11):2425-2428.

[6] Manwell JF, Mcgowan JG, Rogers AL. Wind Energy Explained-theory[M]. Design and Application.Wiley, New York,2002.

[7] 郭 靖, 丁勇鋼, 汪正軍,等.小型風(fēng)力機(jī)：分析、設(shè)計(jì)與應(yīng)用[M]. 北京：機(jī)械工業(yè)出版社, 2013.

[8] Xin Liua, Zhiyong Zhang. Optimization of Astronaut Landing Position Based on Micro Multi-objective Genetic Algorithms[J]. Aerospace Science and Technology, 2013, 29(1): 321-329.

[9] 李建華，魏克湘，王高升，等.小型水平軸風(fēng)力機(jī)葉片氣動(dòng)性能優(yōu)化[J].湖南工程學(xué)院學(xué)報(bào)(自科版)，2016，26(3)：24-27

AnalysisandOptimizationDesignForaSmallWindTurbineTower

LI Jian-hua1, 2, WEI Ke-xiang1, MAO Wen-gui1, WAN Gao-sheng1

(1.College of Mechanical Engineering, Hunan Institute of Engineering, Xiangtan 411104, China;2. Hunan Province Cooperative Innovation Center for Wind Power Equipment and Energy Conversion, Hunan Institute of Engineering, Xiangtan 411104, China)

The tower is an important force component of the wind turbine. The cross section shape, tower height and minimum diameter of the tower have significant influence on the mechanical properties of the tower. It is proposed to optimize the thickness and diameter of the tower based on genetic algorithm, avoid the resonance caused by the excitation of the generator, make the tower the lightest quality and reduce the tower manufacturing cost, which provides a theoretical basis for the design and manufacture of wind turbine tower.

small wind turbine； tower； genetic algorithm； optimize

2017-04-13

國(guó)家自然科學(xué)基金項(xiàng)目(11472103)；湘潭市科技計(jì)劃項(xiàng)目(FJ20164001)；湖南省自然科學(xué)基金資助項(xiàng)目(2016JJ6026).

李建華(1976-),男，碩士，講師，研究方向：機(jī)械設(shè)計(jì)及理論.

TK83

A

1671-119X(2017)04-0026-04