風(fēng)荷載作用下風(fēng)電葉片的受力研究

劉 霞(河北大學(xué)建筑工程學(xué)院，河北保定 071000)

風(fēng)荷載作用下風(fēng)電葉片的受力研究

劉 霞
(河北大學(xué)建筑工程學(xué)院，河北保定 071000)

文章對(duì)在風(fēng)荷載作用下的風(fēng)電葉片進(jìn)行損傷識(shí)別研究，在 ANSYS 中建立有限元模型，對(duì)結(jié)構(gòu)進(jìn)行模態(tài)分析，提取前十階頻率和模態(tài)振型，得出風(fēng)電葉片的振動(dòng)特性。采用振型分解反應(yīng)譜法和譜分析技術(shù)計(jì)算風(fēng)荷載作用下的風(fēng)電葉片的動(dòng)力響應(yīng)，為以后風(fēng)電葉片的實(shí)時(shí)監(jiān)測(cè)做準(zhǔn)備。

風(fēng)電葉片; 模態(tài)分析; 譜分析; ANSYS

風(fēng)能是一種清潔、可再生的綠色能源，具有生態(tài)、環(huán)保、對(duì)環(huán)境無(wú)污染、對(duì)生態(tài)無(wú)破壞等優(yōu)點(diǎn)，對(duì)于人類社會(huì)可持續(xù)發(fā)展具有重要意義。風(fēng)電葉片是捕獲風(fēng)能的最主要部件，也是風(fēng)力發(fā)電機(jī)的關(guān)鍵部件，其質(zhì)量可靠性是保證機(jī)組正常穩(wěn)定運(yùn)行的決定因素。葉片也是風(fēng)力發(fā)電機(jī)組中受力最復(fù)雜的部件，各種激勵(lì)力幾乎都是通過(guò)葉片傳遞出去的，無(wú)論是地球附面層形成的風(fēng)的不均勻流，還是重力影響以及陣風(fēng)等因素，都是作用在葉片上。葉片的根部是葉片結(jié)構(gòu)的關(guān)鍵截面，因?yàn)樵谌~根處的載荷最大，而葉根連接處主要受復(fù)合材料的剪切、擠壓或膠接剪切強(qiáng)度控制，復(fù)合材料的這些強(qiáng)度均低于其抗拉強(qiáng)度，因此葉根是葉片最危險(xiǎn)的部位之一。

本文采用大型通用有限元軟件ANSYS建立三維模型，采用振型分解反應(yīng)譜法計(jì)算風(fēng)荷載作用下的風(fēng)電葉片的動(dòng)力響應(yīng)，為以后風(fēng)電葉片的實(shí)時(shí)監(jiān)測(cè)做準(zhǔn)備。

1 脈動(dòng)風(fēng)速譜[1]

許多風(fēng)工程專家對(duì)水平陣風(fēng)功率譜進(jìn)行了研究，得到了不同形式的風(fēng)速譜表達(dá)式，其中應(yīng)用較為廣泛的是Davenport提出的脈動(dòng)風(fēng)速譜，我國(guó)風(fēng)譜一般取達(dá)文波特風(fēng)譜。Davenport根據(jù)世界上不同地點(diǎn)、不同高度實(shí)測(cè)得到的90多次強(qiáng)風(fēng)記錄，并假定水平陣風(fēng)譜中的湍流積分尺度L沿高度不變，取常數(shù)1 200 m，并取脈動(dòng)風(fēng)速譜為離不同地面高度實(shí)測(cè)值的平均，建立了縱向脈動(dòng)風(fēng)速譜經(jīng)驗(yàn)的數(shù)學(xué)表達(dá)式如下：

式中：SV(n)為脈動(dòng)風(fēng)速功率譜；k為地面粗糙度系數(shù)；n為脈動(dòng)風(fēng)頻率；V1C為標(biāo)準(zhǔn)高度為10 m處的平均風(fēng)速。

2 攝型分解反應(yīng)譜法[4]

振型分解反應(yīng)譜法的流程為模態(tài)分析、譜分析、模態(tài)擴(kuò)展和合并模態(tài)。它通過(guò)理想簡(jiǎn)化的單質(zhì)點(diǎn)體系的最大風(fēng)荷載反應(yīng)來(lái)描述風(fēng)荷載的特性，考慮了當(dāng)時(shí)風(fēng)荷載運(yùn)動(dòng)特性與結(jié)構(gòu)自身的動(dòng)力特性之間的關(guān)系，通過(guò)反應(yīng)譜來(lái)計(jì)算結(jié)構(gòu)動(dòng)力特性(自振周期、振型)所產(chǎn)生的共振效應(yīng)，即將模態(tài)分析結(jié)果與一個(gè)已知的風(fēng)荷載反應(yīng)譜耦合起來(lái)，然后通過(guò)計(jì)算確定結(jié)構(gòu)在風(fēng)載荷作用下的動(dòng)力響應(yīng)。在采用振型分解反應(yīng)譜法對(duì)風(fēng)電葉片進(jìn)行風(fēng)譜效應(yīng)分析時(shí)，風(fēng)電葉片經(jīng)離散化后，可對(duì)運(yùn)動(dòng)狀態(tài)中各節(jié)點(diǎn)按下式建立動(dòng)力平衡微分方程:

振型分解反應(yīng)譜法求解風(fēng)荷載時(shí)，在風(fēng)電葉片模態(tài)分析的基礎(chǔ)上，按照功率譜密度分析方法，取對(duì)風(fēng)電葉片影響最大的Y方向風(fēng)速激勵(lì)進(jìn)行分析。風(fēng)電葉片的阻尼比取為 0.2。將用MATLAB模擬出的Davenport風(fēng)荷載反應(yīng)譜帶入 ANSYS譜分析過(guò)程，采用 SRSS法合并模態(tài)，則可得到風(fēng)電葉片的風(fēng)速譜譜分析結(jié)果。

高層建筑和高聳結(jié)構(gòu)等,一般認(rèn)為這些條件是可以得到滿足的,因而采用SRSS方法可以給出滿足工程需要的結(jié)果。SRSS法相對(duì)于CQC法要簡(jiǎn)潔方便,并且計(jì)算效率要高。正是由于這一點(diǎn),該法長(zhǎng)期以來(lái)一直為許多國(guó)家的抗震和抗風(fēng)設(shè)計(jì)規(guī)范所采用。我國(guó)規(guī)范就是世界上最早采用SRSS法的國(guó)家之一。

此為平方和開平方(SRSS)方法。

3 風(fēng)電葉片實(shí)例

葉片是風(fēng)力機(jī)的關(guān)鍵部件之一,葉片的結(jié)構(gòu)和強(qiáng)度對(duì)風(fēng)力機(jī)的可靠性起重要作用。由于玻璃鋼/復(fù)合材料具有重量輕、耐腐蝕、剛度可設(shè)計(jì)等一系列優(yōu)點(diǎn),國(guó)內(nèi)外已普遍采用玻璃鋼/復(fù)合材料的葉片。

現(xiàn)代風(fēng)力機(jī)通常是采用三葉片的上風(fēng)或下風(fēng)結(jié)構(gòu)。風(fēng)力機(jī)葉片翼型形狀是由飛機(jī)機(jī)翼發(fā)展而來(lái)的，葉片通常由翼型系列組成。常用的翼型有NACA44xx系列等航空翼型。根據(jù)不同的設(shè)計(jì)需要選取翼型。

葉片主要設(shè)計(jì)參數(shù)如表1所示，在Profili軟件中選擇合適的翼型[2]。

表1 風(fēng)機(jī)葉片的主要設(shè)計(jì)參數(shù)

對(duì)于葉片翼展方向鋪層層數(shù)不相同的部位要對(duì)其定義相應(yīng)的材料,也就是說(shuō)在建模時(shí)用多少個(gè)部分拼接成葉面,理論上就要建立多少個(gè)鋪層材料。葉片鋪層數(shù)據(jù)如表2所示。

表2 葉片鋪層數(shù)據(jù)

葉片的鋪層情況非常復(fù)雜，沿葉片縱向分布不均，沿葉片橫向也是不均勻的，這就給單元屬性的賦值帶來(lái)了很大的困難。故采用分段分塊定義層合板模型，然后再分段分塊將模型賦給單元，這樣最大限度地模擬了葉片鋪層的真實(shí)情況。

三維風(fēng)電葉片模型采用通用有限元軟件ANSYS建立，用shell99單元模擬三維不同鋪層，模型如圖1所示。利用ANSYS模態(tài)分析，分別計(jì)算出風(fēng)電葉片葉根部位損傷5 %、10 %、20 %的模態(tài)變化。

圖1 風(fēng)電葉片模型

4 模態(tài)分析與譜分析結(jié)果

通過(guò)不同損傷程度下的模態(tài)分析(表3)可知，在葉片損傷不大的情況下，頻率的變化不明顯。這就說(shuō)明，單單依靠損傷前后模態(tài)的分析，對(duì)葉片的監(jiān)測(cè)還是欠缺的。所以，本文就將模態(tài)分析與譜分析相結(jié)合起來(lái)，更加能夠明確的監(jiān)測(cè)到葉片的損傷位置。

在模態(tài)分析之后，我們引入Davenport 風(fēng)速譜，其功率譜與時(shí)程曲線見圖2、圖3。本文選取 10 頻率及對(duì)應(yīng)下的 Davenport 風(fēng)速譜值如表4所示。

在ANSYS譜分析過(guò)后，我們可以分別得到不同節(jié)點(diǎn)的位移響應(yīng)譜與速度響應(yīng)譜，通過(guò)損傷前后響應(yīng)譜的不同對(duì)比，可以證明出葉片的損傷程度，損傷的越明顯，響應(yīng)譜的差別越大。圖4、圖5是損傷前后的對(duì)比位移響應(yīng)譜。由圖4、圖5可知，損傷后比損傷前的位移值大。

表3 損傷前后頻率變化

圖2 Davenport脈動(dòng)風(fēng)速功率譜

圖3 Davenport風(fēng)速時(shí)程曲線

表4 譜值-頻率對(duì)應(yīng)

圖4 損傷前的位移功率譜

圖5 損傷后的位移功率譜

圖6、圖7是損傷前后的對(duì)比速度響應(yīng)譜。由圖6、圖7可知，損傷后比損傷前的速度變化明顯。

當(dāng)在風(fēng)電葉片受損嚴(yán)重時(shí)，通過(guò)損傷前后模態(tài)分析對(duì)比，可以分析確定風(fēng)電葉片的損傷。當(dāng)風(fēng)電葉片損傷不明顯時(shí)，可以在利用譜分析的理論，通過(guò)分析某點(diǎn)的位移響應(yīng)譜和速度響應(yīng)譜，可以得到損傷前后不同的響應(yīng)譜的結(jié)果，為以后的損傷識(shí)別提供幫助。

5 結(jié)論

本文通過(guò)MATLAB模擬脈動(dòng)風(fēng)速功率譜，并利用ANSYS成功而用振型分解反應(yīng)譜法，將風(fēng)速功率譜加載到風(fēng)電葉片當(dāng)中，并對(duì)風(fēng)電葉片進(jìn)行模態(tài)分析、譜分析、模態(tài)擴(kuò)展和合并模態(tài)，以此方便知道風(fēng)電葉片的應(yīng)力分布以及各階模態(tài)的分布，并且當(dāng)風(fēng)荷載加載到風(fēng)電葉片上時(shí)，可得到實(shí)時(shí)的葉片的位移、速度以及加速度響應(yīng)譜。為以后風(fēng)電葉片的實(shí)時(shí)監(jiān)測(cè)以及在強(qiáng)風(fēng)荷載下風(fēng)電葉片的損傷提供依據(jù)。

圖6 損傷前的速度功率譜

圖7 損傷后的速度功率譜

[1] 王博，風(fēng)荷載作用下的高聳塔架結(jié)構(gòu)的動(dòng)力響應(yīng)[D].蘭州理工大學(xué)，2014.

[2] 林海晨，風(fēng)力機(jī)葉片的有限元建模[J]，綿陽(yáng)師范學(xué)院學(xué)報(bào)，2007,26(8):43-47.

[3] 張春麗，黃爭(zhēng)鳴，董國(guó)華.基于非線性本構(gòu)關(guān)系的復(fù)合材料風(fēng)機(jī)葉片有限元極限分析與設(shè)計(jì)[D].同濟(jì)大學(xué)，2007.

[4] 孔劍.基于ANSYS的鋼閘門地震反應(yīng)譜分析[D].河海大學(xué)，2014.

劉霞(1990～)，女，碩士，研究方向?yàn)榻Y(jié)構(gòu)健康監(jiān)測(cè)與損傷識(shí)別。

TU311.3

A

[定稿日期]2016-08-03