張新榮，程 振，劉成浩，李 鵬

（中水東北勘測設(shè)計研究有限責任公司，吉林長春 130021）

1 概述

風能作為一種分布廣泛、無污染的可再生能源，越來越受到人們的重視。風力發(fā)電是可再生能源技術(shù)中成本降低最快的發(fā)電技術(shù)之一，成為我國可持續(xù)發(fā)展戰(zhàn)略的重要組成部分。大型風力發(fā)電機組大都地處寬闊邊遠地域，如近海和戈壁灘、草原等三北地區(qū)，受高溫、高寒、高濕、大風沙惡劣的自然環(huán)境影響，其機組設(shè)備處于100 m以上的高空，會產(chǎn)生各種各樣復雜的故障[1]。而各種故障中，大部分故障的產(chǎn)生于機組主要結(jié)構(gòu)的振動有關(guān)。對風力發(fā)電機組實現(xiàn)振動檢測，對其運行狀態(tài)進行定量評估，可有計劃地進行設(shè)備維護，降低故障發(fā)生率，從而提高風力發(fā)電機組可靠性，提高風能利用率。振動檢測及故障診斷的最重要依據(jù)是分析結(jié)構(gòu)的固有頻率、特征頻率及振動幅值，即機組的振動特性分析。

文章分析了風力發(fā)電機組振動模式，簡化建立了機組軸向振動數(shù)學模型，推導模型的振動方程，依據(jù)動力學基礎(chǔ)知識及機組固有結(jié)構(gòu)特性計算其前4階固有頻率的解析解。

2 風電機組整機振動模式

風力發(fā)電機組主要振動體可分為：葉片、塔筒及機艙，其中葉片及塔筒為彈性振動體，而機艙的剛度較大，主要以質(zhì)量慣性參與振動。大型風電機組整機振動有3種形式：軸向振動、切向振動及扭轉(zhuǎn)振動，如圖1所示，下面將針對每種振型做振動特性分析。

2.1 切向振動

圖1 風力發(fā)電機組整機振動形式

圖2 風力發(fā)電機組質(zhì)量不平衡引起切向振動

從軸向觀察風力發(fā)電機組，由安裝或初始制造誤差、結(jié)冰或積灰引起的質(zhì)量不平衡，均會導致風機發(fā)生切向振動，如圖2所示。質(zhì)量不平衡故障的產(chǎn)生，究其根本原因，是風在機組每個葉片上的合力不一致，即升力、阻力不一致，如圖3所示。

2.2 扭轉(zhuǎn)振動

圖3 風力發(fā)電機組質(zhì)量不平衡引起切向振動

圖4 風力發(fā)電機組不平衡量引起扭轉(zhuǎn)振動

俯視觀察風力發(fā)電機組，如圖4所示。假設(shè)不平衡量在葉輪上某點處，機組運行產(chǎn)生離心力的水平分量為F，且輪轂與塔筒的中心距離L，這時會產(chǎn)生繞塔筒的扭矩T，從而激起風力發(fā)電機組扭轉(zhuǎn)振動。

2.3 軸向振動

由圖2，3可以看出，質(zhì)量不平衡力只存在于葉輪面內(nèi)，不會引起風力發(fā)電機組的軸向振動，軸向振動是由氣動不平衡引起的。大氣在靠近地表運動產(chǎn)生邊界層，厚度在10 m至數(shù)百米，如圖5所示。風力發(fā)電機組運行于氣流邊界層中，風速在垂直地表方向為典型的層狀分布。

在大型風電場建設(shè)及運行中，由于安裝或初始制造誤差造成的切向振動極少出現(xiàn)，由寒冷地區(qū)葉片結(jié)冰等情況引起的扭轉(zhuǎn)振動偶有發(fā)生，而大氣氣流引起的軸向振動更為常見。本文針風力發(fā)電機組軸向振動進行理論分析，通過求解微分方程獲得機組軸向振動固有頻率的解析解。

3 整機軸向振動固有頻率分析

1）振動方程

在動力分析中視結(jié)構(gòu)為單自由度體系，并假設(shè)塔筒為等厚鋼管制成，其運動方程為：

圖5 軸向振動模型及受力分析圖

假定，風載是平穩(wěn)隨機過程且具有固定力譜密度CF(ω)，響應(yīng)譜SS(ω)對式（1）進行傅氏變換，則有：

2）固有頻率

結(jié)構(gòu)簡化圖見圖6，塔頂機艙中心設(shè)為坐標原點。

不考慮內(nèi)部及外部空氣摩擦力因素，則有：

圖6 軸向振動模型結(jié)構(gòu)簡化圖

若傳遞函數(shù)L(m)分母為零，即ω=K，則響應(yīng)譜SS(ω)的值為無窮大，因此風電機組軸向振動的一階固有頻率為：

由結(jié)構(gòu)力學可知，距離坐標原點x處的彎矩為：

積分一次得到轉(zhuǎn)角：

式中：E為彈性模數(shù)，E=2.06×1011Pa；I為慣性矩；g為重力加速度；G為當量重，G=GD+GT/3，GD為塔頂重量，GT為塔筒重量；L為塔頂?shù)降孛婢嚯x。

5 結(jié)語

本文分析了風力發(fā)電機組整機振動3種不同振動模式，簡化建立了機組軸向振動數(shù)學模型，依據(jù)動力學基礎(chǔ)知識簡化推導，并最終建立模型的振動方程。機組振動可簡化為單自由度振動系統(tǒng)模型，求解振動方程可得到軸向振動固有頻率特性。

［1］郭新生.風能利用技術(shù)［M］.北京：化學工業(yè)出版社，2007.

［2］芮曉明，柳亦兵，馬志勇.風力發(fā)電機組設(shè)計［M］.北京：機械工業(yè)出版社，2010.

［3］陳云程，陳孝耀，朱成名.風力機設(shè)計與應(yīng)用［M］.上海：科學技術(shù)出版社，1990.

［4］姚興佳，張雅楠，郭慶鼎.風力發(fā)電機組的多模變槳距控制［J］.沈陽工業(yè)大學學報，2009（2）：28—31.

［5］屈維德.機械振動手冊［M］.北京：機械工業(yè)出版社，1992.

［6］勒古里雷斯（LeGourieres，D）.風力機的理論與設(shè)計，1990.

［7］吳運東.ZDF13-20kW風力發(fā)電機風振計算［J］.太陽能學報，1988（1）：9—12.

［8］林起孝.FD6BH型3×4kW風力發(fā)電機組塔架計算［J］.福建機械，1984（1）.

［9］Davenport,A.G.,Thespectrumofhorizontalgustiness near the ground in high winds.Quart.J.Royal Mctcoral.Soc，1961，187：194-211.