程友良,薛占璞

(華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

風(fēng)力發(fā)電塔架振動(dòng)特性聯(lián)合仿真應(yīng)用研究

程友良,薛占璞

(華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

風(fēng)力機(jī)塔架是整個(gè)風(fēng)力發(fā)電機(jī)組的支撐部件，其承載能力決定風(fēng)力發(fā)電機(jī)組的安全性與正常運(yùn)行。針對(duì)風(fēng)力發(fā)電塔架振動(dòng)引起的風(fēng)力發(fā)電機(jī)組失穩(wěn)問(wèn)題，運(yùn)用聯(lián)合仿真技術(shù)，在Solidworks中建立塔架三維結(jié)構(gòu)動(dòng)力學(xué)模型。利用有限元分析軟件ANSYS對(duì)塔架的振動(dòng)特性及風(fēng)致響應(yīng)進(jìn)行分析，提取模態(tài)分析的前五階非零模態(tài)振型，對(duì)各階的位移與振型圖進(jìn)行分析。針對(duì)塔架在極端環(huán)境下動(dòng)力響應(yīng)變化情況，進(jìn)行脈動(dòng)風(fēng)速下塔架應(yīng)力變化分析。風(fēng)致響應(yīng)是衡量塔架正常運(yùn)行的重要指標(biāo)，應(yīng)力變化幅度決定塔架強(qiáng)度是否滿足強(qiáng)度要求，以防止塔架結(jié)構(gòu)斷裂。分析結(jié)果表明，塔架前三階振型為橫向彎曲，無(wú)明顯位移偏移量；第四階振型出現(xiàn)橫向彎曲和扭轉(zhuǎn)變形，振動(dòng)位移和應(yīng)力變化幅度較大；第五階振型出現(xiàn)扭轉(zhuǎn)。該結(jié)果為塔架的振動(dòng)特性提供了依據(jù)，從而避免了風(fēng)力發(fā)電機(jī)組失穩(wěn)現(xiàn)象的發(fā)生。

風(fēng)力發(fā)電塔架； 振動(dòng)特性； 位移； 有限元分析； 聯(lián)合仿真； 結(jié)構(gòu)動(dòng)力學(xué)

0 引言

隨著煤、石油等傳統(tǒng)能源日益減少及環(huán)境污染引發(fā)社會(huì)、生活等各類問(wèn)題，風(fēng)力發(fā)電作為一種清潔、無(wú)污染的可再生能源技術(shù)，越來(lái)越受到世界各國(guó)的重視，已成為全球應(yīng)對(duì)能源危機(jī)的重要舉措。風(fēng)力機(jī)單機(jī)容量從6 MW逐漸趨于10 MW，甚至15 MW超大容量。隨著單機(jī)容量的逐漸增大，塔架高度也隨之增加。塔架在風(fēng)力發(fā)電系統(tǒng)中除了承受風(fēng)力機(jī)及機(jī)艙重力之外，還承受風(fēng)切變、振動(dòng)等載荷作用，并且載荷的隨機(jī)性與交變性更為顯著[1-2]，塔架振動(dòng)特性是保證風(fēng)力發(fā)電系統(tǒng)穩(wěn)定運(yùn)行的關(guān)鍵。

近年來(lái)，許多學(xué)者對(duì)風(fēng)力機(jī)塔架進(jìn)行振動(dòng)特性分析，對(duì)錐筒式和分段式塔架的模態(tài)振型圖進(jìn)行比較，并根據(jù)各個(gè)模態(tài)分析結(jié)果，選擇一種最佳方案[2-3]。但是以往的塔架振動(dòng)特性分析多采用單一仿真技術(shù)，而較少采用聯(lián)合仿真分析技術(shù)。為準(zhǔn)確分析風(fēng)力機(jī)塔架振動(dòng)特性，本文利用聯(lián)合仿真技術(shù)，基于振動(dòng)特性基本理論，通過(guò)有限元分析軟件Solidworks及ANSYS對(duì)風(fēng)力機(jī)塔架的振動(dòng)特性及脈動(dòng)風(fēng)速下的響應(yīng)情況進(jìn)行分析，得到塔架在不同頻率下的振型圖及應(yīng)力變化分布，從而為塔架的應(yīng)用及參數(shù)在線監(jiān)測(cè)提供參考。

1 振動(dòng)特性基本理論

(1)

式中:M為塔架結(jié)構(gòu)的整體質(zhì)量矩陣；K為塔架結(jié)構(gòu)的整體剛度矩陣。

非零矢量φi的比例解為：

φi=[φ1iφ2iφ3iφ4i]T

(2)

將n個(gè)特征矢量φi按列排成n×n階陣，可得系統(tǒng)的特征矢量矩陣：

φ=[φ1φ2…φn]

(3)

此時(shí)特征矢量為模態(tài)矢量或模態(tài)矩陣，即為模態(tài)振型。

2 風(fēng)力機(jī)塔架結(jié)構(gòu)動(dòng)力學(xué)模型

利用Solidworks建立1.5 MW風(fēng)力機(jī)塔架三維結(jié)構(gòu)動(dòng)力學(xué)模型，底部直徑為4.2 m，頂部直徑為2.5 m，高度為70 m，葉片和機(jī)艙的質(zhì)量為78 300 kg，塔架材料為合金鋼，彈性模量E=2.0×1 011 Pa，泊松比μ=0.3，材料密度ρ=7.85×103kg/m3，塔架可分為兩段，法蘭直徑為4 m。為了保證風(fēng)力機(jī)塔架振動(dòng)特性分析與實(shí)際結(jié)構(gòu)相一致，增加風(fēng)力機(jī)及機(jī)艙，構(gòu)成風(fēng)力發(fā)電系統(tǒng)。塔架三維結(jié)構(gòu)動(dòng)力學(xué)模型如圖1所示。

圖1 塔架三維結(jié)構(gòu)動(dòng)力學(xué)模型

根據(jù)塔架的實(shí)際情況作出如下簡(jiǎn)化。

①省略塔架連接處螺栓及法蘭連接的配合情況，近似剛性連接;

②不考慮螺栓裝配的預(yù)應(yīng)力;

Dynamic Simulation of Impact of Shipping Services on Port

③簡(jiǎn)化對(duì)于風(fēng)力機(jī)及機(jī)艙不重要的裝配孔;

④地面與塔架的連接剛度無(wú)限大，塔架與地面連接進(jìn)行六自由度全約束。

根據(jù)結(jié)構(gòu)簡(jiǎn)化要求，在Solidworks中進(jìn)行模型的修改，格式保存成“.x-t”，并通過(guò)接口技術(shù)導(dǎo)入到ANSYS中。設(shè)置單元類型為六面體三維實(shí)體單元solid45，每個(gè)節(jié)點(diǎn)有三個(gè)移動(dòng)自由度[6-8]。根據(jù)簡(jiǎn)化的風(fēng)力塔架三維結(jié)構(gòu)動(dòng)力學(xué)模型，其計(jì)算域如圖2所示。

圖2 風(fēng)力機(jī)塔架計(jì)算域

計(jì)算域大小為長(zhǎng)40 m、寬30 m、高120 m，塔架模型位于計(jì)算域的中心位置。在網(wǎng)格劃分方面，塔架壁面的網(wǎng)格質(zhì)量起著至關(guān)重要的作用，塔架網(wǎng)格數(shù)量約為3.84×105。

3 塔架振動(dòng)特性分析

3.1 模態(tài)分析

在塔架結(jié)構(gòu)模態(tài)分析中，前五階非零模態(tài)振型是所有振型的基礎(chǔ)，與其對(duì)應(yīng)的模態(tài)頻率是振動(dòng)特性分析的重點(diǎn)[9-11]。振型圖表明：前三階振型中無(wú)明顯扭轉(zhuǎn)變形，第四階模態(tài)振型顯示徑向擠壓及扭轉(zhuǎn)變形。第一階、第二階振型為擺振，第三階振型為揮舞。根據(jù)振動(dòng)特性基本理論，在振動(dòng)過(guò)程中塔架的能量主要在第一階、第二階頻率處，塔架的主要振動(dòng)是擺振，而非扭振。

塔架前五階固有頻率、振型和總振幅如表1所示。由表1可知，前三階模態(tài)振型主要是橫向彎曲，變形量不大；第四階模態(tài)振型出現(xiàn)了縱向彎曲及扭轉(zhuǎn)，總振幅為0.115 m，變形最大的節(jié)點(diǎn)位于塔架的頂部處，屬于扭曲變形，應(yīng)避免以該頻率(f=3.488 Hz)工作，否則會(huì)產(chǎn)生由于振動(dòng)位移過(guò)大發(fā)生塔架結(jié)構(gòu)斷裂。隨著頻率的增加，塔架振型由橫向彎曲逐漸轉(zhuǎn)變?yōu)榕まD(zhuǎn)振動(dòng)。實(shí)際運(yùn)行時(shí)應(yīng)根據(jù)頻率及對(duì)應(yīng)的振型進(jìn)行監(jiān)測(cè)，避免由于振動(dòng)現(xiàn)象導(dǎo)致的破壞。

表1 塔架前五階固有頻率、振型和總振幅

3.2 脈動(dòng)風(fēng)載荷作用下動(dòng)力響應(yīng)分析

為了分析風(fēng)力機(jī)塔架在極端環(huán)境下的動(dòng)力響應(yīng)情況，對(duì)其應(yīng)力變化進(jìn)行分析，以保證塔架足夠的強(qiáng)度適應(yīng)風(fēng)力發(fā)電系統(tǒng)運(yùn)行。施加脈動(dòng)風(fēng)載荷，數(shù)值計(jì)算采用Simple算法，設(shè)置殘差小于10-6，對(duì)塔架周圍風(fēng)載荷滿足的連續(xù)性方程、動(dòng)量方程、能量方程等參數(shù)趨于漸近值時(shí)達(dá)到收斂。

塔架應(yīng)力變化曲線如圖3所示。

圖3 塔架應(yīng)力變化曲線

由圖3可知：最大應(yīng)力為45 kPa，對(duì)應(yīng)的頻率為2.39 Hz；最小應(yīng)力為30 kPa，對(duì)應(yīng)的頻率為3.81 Hz。塔架在30 s內(nèi)應(yīng)力變化幅度不大，屬于允許范圍內(nèi)，不易發(fā)生應(yīng)力集中以及振動(dòng)破壞。從對(duì)應(yīng)的頻率與振動(dòng)位移量可知，模態(tài)頻率有10%裕度，不易發(fā)生共振。

4 結(jié)束語(yǔ)

結(jié)合風(fēng)力機(jī)塔架結(jié)構(gòu)參數(shù)，利用聯(lián)合仿真技術(shù)，在Solidworks建立風(fēng)力機(jī)塔架模型，并通過(guò)ANSYS對(duì)其進(jìn)行模態(tài)分析及極端載荷下的應(yīng)力變化分析。

由模態(tài)分析可知塔架的固有頻率、振型及總振幅，振動(dòng)特性頻率主要集中在0.357～3.726 Hz,最大變形量為0.115 m，應(yīng)力變化平緩，不易發(fā)生應(yīng)力集中。

振動(dòng)特性結(jié)果為塔架的設(shè)計(jì)及運(yùn)行狀態(tài)的監(jiān)測(cè)提供了理論依據(jù)，塔架的動(dòng)力學(xué)分析是未來(lái)研究的重點(diǎn)。

[1] 張豐豪,何榕.結(jié)構(gòu)阻尼對(duì)風(fēng)力機(jī)塔架振動(dòng)特性的影響[J].太陽(yáng)能學(xué)報(bào),2015,36(10):2467-2473.

[2] 張偉.基于ANSYS的糊底機(jī)機(jī)架模態(tài)分析[J].自動(dòng)化儀表,2016,37(12):23-28.

[3] 趙榮珍,呂鋼.水平軸風(fēng)力發(fā)電機(jī)塔架的振動(dòng)模態(tài)分析[J].蘭州理工大學(xué)學(xué)報(bào),2009,35(2):33-36.

[4] 黃珊秋,陸萍.ZONDZ-40風(fēng)力機(jī)塔架的模態(tài)分析[J].太陽(yáng)能學(xué)報(bào),2001,22(2):153-156.

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

[6] 周長(zhǎng)城,胡文喜,熊文波.ANSYS基礎(chǔ)與典型范例[M].北京：電子工業(yè)出版社,2007.

[7] 柯世堂,王同光,曹九發(fā),等.考慮葉片旋轉(zhuǎn)和離心力效應(yīng)風(fēng)力機(jī)塔架風(fēng)振分析[J].太陽(yáng)能學(xué)報(bào),2015,36(1):33-40.

[8] 應(yīng)有,朱重喜,楊帆,等.大型風(fēng)電機(jī)組塔架主動(dòng)阻尼控制技術(shù)研究[J].太陽(yáng)能學(xué)報(bào),2015,36(1):54-60.

[9] 戴建鑫.風(fēng)力機(jī)塔架的有限元建模及靜動(dòng)態(tài)特性分析[D].蘭州:蘭州理工大學(xué),2011.

[10]王振宇,張彪,趙艷,等.臺(tái)風(fēng)作用下風(fēng)力機(jī)塔架振動(dòng)響應(yīng)研究[J].太陽(yáng)能學(xué)報(bào),2013,34(8):1434-1442.

[11]閆梁.風(fēng)力發(fā)電機(jī)塔架結(jié)構(gòu)風(fēng)壓數(shù)值模擬研究[D].包頭:內(nèi)蒙古科技大學(xué),2014.

Study on Application of Joint Simulation of Wind Power Generation Tower Vibration Performance

CHENG Youliang,XUE Zhanpu

(School of Energy Power and Mechanical Engineering,North China Electric Power University,Baoding 071003,China)

Wind turbine tower is the supporting component of the whole wind power generation unit; its bearing capacity determines the safety and normal operation of wind power unit.To solve the instability problem of wind power unit which is caused by the vibration of wind power generation tower,by using joint simulation technology,the three-dimensional structural dynamics model of tower is established in Solidworks.The vibration characteristics of the tower and the wind-induced response are analyzed by adopting the finite element analysis software ANSYS,the first five orders non-zero modal vibration types of modal analysis are extracted,and the displacement and vibration pattern of each order are analyzed.In accordance with the dynamic response variations of tower under extreme environment,the stress changes of tower under pulsed wind speed are analyzed.To avoid the fracture of tower,the wind-induced response is the important index to determine whether the strength of tower meets the requirement or not.The analysis results show that the first three orders vibration types of tower are lateral bending and no obvious shift; the fourth vibration type is bending and torsion deformation,its vibration displacement and the stress changes are large; torsion may occur in the fifth order vibration type.The results provide the basis for the vibration characteristics of the tower,so as to avoid the occurrence of instability phenomenon of wind power generation unit.

Wind power generation tower; Vibration characteristics; Displacement; Finite element analysis; Joint simulation; Structural dynamics

中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)基金資助項(xiàng)目(2016XS107)

程友良(1963—)，男，博士，教授，主要從事流體動(dòng)力學(xué)理論及其應(yīng)用、流體設(shè)備與節(jié)能、可再生能源理論及其應(yīng)用、清潔能源利用技術(shù)與設(shè)備等方向的研究。E-mail:shenghuo166@163.com。

TH-39;TP29

A

10.16086/j.cnki.issn1000-0380.201706006

修改稿收到日期：2017-03-10