張新榮,程 振,劉成浩,李 鵬
(中水東北勘測(cè)設(shè)計(jì)研究有限責(zé)任公司,吉林長(zhǎng)春 130021)
風(fēng)能作為一種分布廣泛、無污染的可再生能源,越來越受到人們的重視。風(fēng)力發(fā)電是可再生能源技術(shù)中成本降低最快的發(fā)電技術(shù)之一,成為我國可持續(xù)發(fā)展戰(zhàn)略的重要組成部分。大型風(fēng)力發(fā)電機(jī)組大都地處寬闊邊遠(yuǎn)地域,如近海和戈壁灘、草原等三北地區(qū),受高溫、高寒、高濕、大風(fēng)沙惡劣的自然環(huán)境影響,其機(jī)組設(shè)備處于100 m以上的高空,會(huì)產(chǎn)生各種各樣復(fù)雜的故障[1]。而各種故障中,大部分故障的產(chǎn)生于機(jī)組主要結(jié)構(gòu)的振動(dòng)有關(guān)。對(duì)風(fēng)力發(fā)電機(jī)組實(shí)現(xiàn)振動(dòng)檢測(cè),對(duì)其運(yùn)行狀態(tài)進(jìn)行定量評(píng)估,可有計(jì)劃地進(jìn)行設(shè)備維護(hù),降低故障發(fā)生率,從而提高風(fēng)力發(fā)電機(jī)組可靠性,提高風(fēng)能利用率。振動(dòng)檢測(cè)及故障診斷的最重要依據(jù)是分析結(jié)構(gòu)的固有頻率、特征頻率及振動(dòng)幅值,即機(jī)組的振動(dòng)特性分析。
文章分析了風(fēng)力發(fā)電機(jī)組振動(dòng)模式,簡(jiǎn)化建立了機(jī)組軸向振動(dòng)數(shù)學(xué)模型,推導(dǎo)模型的振動(dòng)方程,依據(jù)動(dòng)力學(xué)基礎(chǔ)知識(shí)及機(jī)組固有結(jié)構(gòu)特性計(jì)算其前4階固有頻率的解析解。
風(fēng)力發(fā)電機(jī)組主要振動(dòng)體可分為:葉片、塔筒及機(jī)艙,其中葉片及塔筒為彈性振動(dòng)體,而機(jī)艙的剛度較大,主要以質(zhì)量慣性參與振動(dòng)。大型風(fēng)電機(jī)組整機(jī)振動(dòng)有3種形式:軸向振動(dòng)、切向振動(dòng)及扭轉(zhuǎn)振動(dòng),如圖1所示,下面將針對(duì)每種振型做振動(dòng)特性分析。
圖1 風(fēng)力發(fā)電機(jī)組整機(jī)振動(dòng)形式
圖2 風(fēng)力發(fā)電機(jī)組質(zhì)量不平衡引起切向振動(dòng)
從軸向觀察風(fēng)力發(fā)電機(jī)組,由安裝或初始制造誤差、結(jié)冰或積灰引起的質(zhì)量不平衡,均會(huì)導(dǎo)致風(fēng)機(jī)發(fā)生切向振動(dòng),如圖2所示。質(zhì)量不平衡故障的產(chǎn)生,究其根本原因,是風(fēng)在機(jī)組每個(gè)葉片上的合力不一致,即升力、阻力不一致,如圖3所示。
圖3 風(fēng)力發(fā)電機(jī)組質(zhì)量不平衡引起切向振動(dòng)
圖4 風(fēng)力發(fā)電機(jī)組不平衡量引起扭轉(zhuǎn)振動(dòng)
俯視觀察風(fēng)力發(fā)電機(jī)組,如圖4所示。假設(shè)不平衡量在葉輪上某點(diǎn)處,機(jī)組運(yùn)行產(chǎn)生離心力的水平分量為F,且輪轂與塔筒的中心距離L,這時(shí)會(huì)產(chǎn)生繞塔筒的扭矩T,從而激起風(fēng)力發(fā)電機(jī)組扭轉(zhuǎn)振動(dòng)。
由圖2,3可以看出,質(zhì)量不平衡力只存在于葉輪面內(nèi),不會(huì)引起風(fēng)力發(fā)電機(jī)組的軸向振動(dòng),軸向振動(dòng)是由氣動(dòng)不平衡引起的。大氣在靠近地表運(yùn)動(dòng)產(chǎn)生邊界層,厚度在10 m至數(shù)百米,如圖5所示。風(fēng)力發(fā)電機(jī)組運(yùn)行于氣流邊界層中,風(fēng)速在垂直地表方向?yàn)榈湫偷膶訝罘植肌?/p>
在大型風(fēng)電場(chǎng)建設(shè)及運(yùn)行中,由于安裝或初始制造誤差造成的切向振動(dòng)極少出現(xiàn),由寒冷地區(qū)葉片結(jié)冰等情況引起的扭轉(zhuǎn)振動(dòng)偶有發(fā)生,而大氣氣流引起的軸向振動(dòng)更為常見。本文針風(fēng)力發(fā)電機(jī)組軸向振動(dòng)進(jìn)行理論分析,通過求解微分方程獲得機(jī)組軸向振動(dòng)固有頻率的解析解。
1)振動(dòng)方程
在動(dòng)力分析中視結(jié)構(gòu)為單自由度體系,并假設(shè)塔筒為等厚鋼管制成,其運(yùn)動(dòng)方程為:
圖5 軸向振動(dòng)模型及受力分析圖
假定,風(fēng)載是平穩(wěn)隨機(jī)過程且具有固定力譜密度CF(ω),響應(yīng)譜SS(ω)對(duì)式(1)進(jìn)行傅氏變換,則有:
2)固有頻率
結(jié)構(gòu)簡(jiǎn)化圖見圖6,塔頂機(jī)艙中心設(shè)為坐標(biāo)原點(diǎn)。
不考慮內(nèi)部及外部空氣摩擦力因素,則有:
圖6 軸向振動(dòng)模型結(jié)構(gòu)簡(jiǎn)化圖
若傳遞函數(shù)L(m)分母為零,即ω=K,則響應(yīng)譜SS(ω)的值為無窮大,因此風(fēng)電機(jī)組軸向振動(dòng)的一階固有頻率為:
由結(jié)構(gòu)力學(xué)可知,距離坐標(biāo)原點(diǎn)x處的彎矩為:
積分一次得到轉(zhuǎn)角:
式中:E為彈性模數(shù),E=2.06×1011Pa;I為慣性矩;g為重力加速度;G為當(dāng)量重,G=GD+GT/3,GD為塔頂重量,GT為塔筒重量;L為塔頂?shù)降孛婢嚯x。
本文分析了風(fēng)力發(fā)電機(jī)組整機(jī)振動(dòng)3種不同振動(dòng)模式,簡(jiǎn)化建立了機(jī)組軸向振動(dòng)數(shù)學(xué)模型,依據(jù)動(dòng)力學(xué)基礎(chǔ)知識(shí)簡(jiǎn)化推導(dǎo),并最終建立模型的振動(dòng)方程。機(jī)組振動(dòng)可簡(jiǎn)化為單自由度振動(dòng)系統(tǒng)模型,求解振動(dòng)方程可得到軸向振動(dòng)固有頻率特性。
[1]郭新生.風(fēng)能利用技術(shù)[M].北京:化學(xué)工業(yè)出版社,2007.
[2]芮曉明,柳亦兵,馬志勇.風(fēng)力發(fā)電機(jī)組設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,2010.
[3]陳云程,陳孝耀,朱成名.風(fēng)力機(jī)設(shè)計(jì)與應(yīng)用[M].上海:科學(xué)技術(shù)出版社,1990.
[4]姚興佳,張雅楠,郭慶鼎.風(fēng)力發(fā)電機(jī)組的多模變槳距控制[J].沈陽工業(yè)大學(xué)學(xué)報(bào),2009(2):28—31.
[5]屈維德.機(jī)械振動(dòng)手冊(cè)[M].北京:機(jī)械工業(yè)出版社,1992.
[6]勒古里雷斯(LeGourieres,D).風(fēng)力機(jī)的理論與設(shè)計(jì),1990.
[7]吳運(yùn)東.ZDF13-20kW風(fēng)力發(fā)電機(jī)風(fēng)振計(jì)算[J].太陽能學(xué)報(bào),1988(1):9—12.
[8]林起孝.FD6BH型3×4kW風(fēng)力發(fā)電機(jī)組塔架計(jì)算[J].福建機(jī)械,1984(1).
[9]Davenport,A.G.,Thespectrumofhorizontalgustiness near the ground in high winds.Quart.J.Royal Mctcoral.Soc,1961,187:194-211.