基于箱型梁骨架式風(fēng)力機葉片性能研究

姜　艷，王　瑋，霍尹才，孫曉明，韓　梅

(1.黑龍江大學(xué) 水利電力學(xué)院，哈爾濱 150080；2.哈爾濱工程大學(xué) 船舶工程學(xué)院，哈爾濱 150001)

?

基于箱型梁骨架式風(fēng)力機葉片性能研究

姜艷1，王瑋2,*，霍尹才1，孫曉明1，韓梅1

(1.黑龍江大學(xué) 水利電力學(xué)院，哈爾濱 150080；2.哈爾濱工程大學(xué) 船舶工程學(xué)院，哈爾濱 150001)

摘要：葉片是風(fēng)力發(fā)電機組的關(guān)鍵部件，風(fēng)力機葉片所使用的材料將直接影響風(fēng)力機的性能。在探索新型碳纖維復(fù)合材料葉片的基礎(chǔ)上，通過箱型梁骨架式葉片有限元分析，研究復(fù)合材料風(fēng)力機葉片在載荷作用下的變形及其力學(xué)特性，進而揭示碳纖維復(fù)合材料的鋪層角度對變形和力學(xué)特性的影響規(guī)律，以期為新材料在風(fēng)力機葉片上的應(yīng)用提供技術(shù)支持和理論依據(jù)，對于提升我國自主設(shè)計和制造長壽命大型風(fēng)力機葉片的能力具有重要的科學(xué)價值和工程應(yīng)用意義。

關(guān)鍵詞：風(fēng)力機葉片；碳纖維復(fù)合材料；性能；有限元分析

0引言

風(fēng)能是可再生能源中最具活力的一種，己成為發(fā)展最快的能源[1]，受到世界各國的關(guān)注[2]。因此對高性能風(fēng)力機及其部件的研發(fā)十分迫切[3]。然而就風(fēng)力發(fā)電機組而言，葉片是其重要組成部件，風(fēng)力機葉片所使用的材料將直接影響風(fēng)力機的性能，且風(fēng)力機葉片性能的好壞將會直接影響到整個風(fēng)力機組運行的效率和穩(wěn)定性[4]。風(fēng)力機葉片受力比較復(fù)雜，再加上外形比較細長，屬于容易發(fā)生形變的柔性彈性體[5]。因此，研究葉片其性能和壽命在極端外界環(huán)境下的變化也逐漸成為研發(fā)人員關(guān)注的焦點。目前國內(nèi)外所采用的風(fēng)力機的葉片為玻纖增強聚脂樹脂、玻纖增強環(huán)氧樹脂[6-10]，但由于受極端溫度的影響，致使葉片結(jié)構(gòu)產(chǎn)生反復(fù)變形，從而影響風(fēng)力機的氣動性能與結(jié)構(gòu)性能，并且隨著風(fēng)力機組的大型化量產(chǎn)，這種形式的葉片將無法滿足要求，必須開發(fā)更為先進的材料，具備輕質(zhì)、高強以及剛性好等性能。為此，本文在探索新型碳纖維復(fù)合材料葉片的基礎(chǔ)上，通過箱型梁骨架式葉片有限元分析，研究復(fù)合材料風(fēng)力機葉片在載荷作用的變形及其力學(xué)特性，進而揭示碳纖維復(fù)合材料的鋪層角度對變形和力學(xué)特性的影響規(guī)律，以期為新材料在風(fēng)力機葉片上的應(yīng)用提供技術(shù)支持和理論依據(jù)，對于提升我國自主設(shè)計和制造長壽命大型風(fēng)力機葉片的能力具有重要的科學(xué)價值和工程應(yīng)用意義。

黑龍江省風(fēng)電產(chǎn)業(yè)將在“一帶一路”建設(shè)戰(zhàn)略的重要指引下，風(fēng)能行業(yè)走出國門將打開一個全新的市場，風(fēng)力發(fā)電作為成熟的代表性新能源產(chǎn)業(yè)，也將成為能源基礎(chǔ)設(shè)施建設(shè)的一大載體，故研究高寒地區(qū)氣候?qū)︼L(fēng)力發(fā)電機葉片性能的影響具有一定的實際應(yīng)用價值。

1風(fēng)力發(fā)電機葉片

1.1葉片有限元模型

風(fēng)力機葉片的幾何模型較為復(fù)雜?；贏NSYS Workbench中的DM建模模塊，本文采取使用自上而下的方式構(gòu)建葉片模型，首先建立葉片實體三維模型，在 DM 中建立葉片曲面參數(shù)模型，并將葉片建模前先將其分段，然后確定好截面數(shù)后建立每個截面上的關(guān)鍵點，再將相應(yīng)的點連接成線，最后將相應(yīng)的線連成面，進而實現(xiàn)葉片的三維實體建模。

1.2箱型梁骨架式葉片有限元分析

箱型梁骨架式葉片由復(fù)合材料蒙皮、箱型梁、隔板、金屬軸、法蘭盤5部分組成(圖1)，對葉片各構(gòu)件沿展長方向厚度進行了分區(qū)設(shè)計(表1)。此外，葉片分析時仍采用Solid187實體單元和Shell181殼單元進行網(wǎng)格劃分，共劃分了89 419個節(jié)點和78 553個單元。

蒙皮角度/(°)09045-45-4545900①區(qū)700～1800mm層數(shù)厚度/mm4422111111112244②區(qū)1800～3800mm層數(shù)厚度/mm53.511111111111153.5③區(qū)3800～6000mm層數(shù)厚度/mm53.511111111111153.5

注：沿葉片展長方向為0°,沿葉片環(huán)向為90°。

(b) 隔板分區(qū)設(shè)計

(c)箱型梁分區(qū)設(shè)計

(d)金屬軸分區(qū)設(shè)計

2計算結(jié)果分析

2.1蒙皮計算結(jié)果

通過對蒙皮應(yīng)力云圖和坐標位置-應(yīng)力關(guān)系的分析可知(圖2和圖3)，葉片最大應(yīng)力為198.84 MPa，出現(xiàn)在距葉根6 000 mm的葉尖處，蒙皮應(yīng)力較大區(qū)集中在①區(qū)，不會發(fā)生葉片根部的斷裂問題，使葉片與主體部分連接處具有很好的穩(wěn)定性能。

圖2　蒙皮應(yīng)力云圖Fig.2　Skin stress

圖3　蒙皮坐標位置-應(yīng)力關(guān)系Fig.3　The relation of position-stress of skin

2.2箱型梁計算結(jié)果

通過箱型梁上下翼緣及腹板應(yīng)力云圖和坐標位置-應(yīng)力關(guān)系的分析可知(圖4～圖7)：翼緣應(yīng)力最大值為258.4 MPa，出現(xiàn)在葉片第55 702個節(jié)點上，距葉根1 200 mm處，箱型梁翼緣應(yīng)力總體較大，說明其承受了較大部分的彎矩作用。腹板最大應(yīng)力為147.3 MPa出現(xiàn)在葉片第2 342個節(jié)點上距葉根3 800 mm處，該處應(yīng)力較大，是因為此處腹板厚度發(fā)生了梯度性變化，且腹板受力較均勻。由整體看，腹板、翼緣兩者共同受力，協(xié)調(diào)變形。

圖4　箱梁翼緣應(yīng)力云圖 Fig.4　Flange stress

圖5　箱梁翼緣坐標位置-應(yīng)力關(guān)系 Fig.5　The relation of position-stress of flange

圖6　箱梁腹板應(yīng)力云圖 Fig.6　Web stress

圖7　箱梁腹板坐標位置-應(yīng)力關(guān)系Fig.7　The relation of position-stress of web

2.3隔板計算結(jié)果

借助于隔板應(yīng)力云圖和隔板坐標位置-應(yīng)力關(guān)系分析可知(圖8和圖9)：該處應(yīng)力最大值為378 MPa出現(xiàn)在葉片第82 449個節(jié)點處(隔板與軸接觸處)，其值小于材料設(shè)計值，設(shè)計有效。

圖8　隔板應(yīng)力云圖 Fig.8　sect1ition stress

圖9　隔板坐標位置-應(yīng)力關(guān)系Fig.9　The relation of position-stress of sect1ition

2.4金屬軸計算結(jié)果

金屬軸應(yīng)力云圖見圖10，通過分析可知，金屬軸最大應(yīng)力為345.54 MPa，出現(xiàn)在葉片第2 311個節(jié)點上，距葉根955 mm處。通過對坐標位置-應(yīng)力關(guān)系圖分析(圖11)，可見在每一區(qū)段厚度改變處，應(yīng)力都有相應(yīng)的突變增大，故在分段處應(yīng)做加強處理。

圖10　金屬軸應(yīng)力云圖 Fig.10　Metal shaft stress

圖11　金屬軸坐標位置-應(yīng)力關(guān)系Fig.11　The relation of position-stress of metal shaft

相對于箱形梁骨架式風(fēng)力機葉片，全復(fù)合材料式葉片也進行了分析，計算得到該方案下葉片結(jié)構(gòu)的整體應(yīng)力見圖12、圖13。由有限元分析結(jié)果可知葉片最大應(yīng)力為187.37 MPa，出現(xiàn)在葉片蒙皮葉根過渡段第20 676個節(jié)點處，距葉根為972.45 mm，與結(jié)構(gòu)鋼葉片相比最大應(yīng)力降低了46%，說明全復(fù)合材料葉片能有效避免應(yīng)力集中現(xiàn)象的發(fā)生。從葉片整體應(yīng)力、應(yīng)變水平來看，應(yīng)力、應(yīng)變主要集中在葉根過渡和靠近葉根過渡區(qū)域，而葉尖區(qū)域應(yīng)力、應(yīng)變較小，同結(jié)構(gòu)鋼葉片受力相符，但全復(fù)合材料葉片整體受力比較均勻。對于葉片內(nèi)腹板，主要承受剪力作用，受力較為均勻，其最大應(yīng)力為40.51 MPa，出現(xiàn)在腹板第3 449個節(jié)點處，距葉根為1 213.8 mm，應(yīng)力、應(yīng)變水平相對較低。同時腹板內(nèi)側(cè)蒙皮受力較大，而腹板外側(cè)蒙皮受力較小。但是全復(fù)合材料葉片直接制作成本在2 969左右，其制作周期約為213.5 h。而箱型梁骨架式直接制作成本在1 884左右，其制作周期約為111.9 h。

圖12　復(fù)合材料葉片應(yīng)力云圖 Fig. 12　Composite blade stress

圖13　蒙皮應(yīng)力云圖Fig.13　Skin stress

在綜合考慮葉片受力性能、材料利用率、葉片結(jié)構(gòu)自重、造價等多種因素的基礎(chǔ)上，采用選取玻璃纖維增強樹脂基復(fù)合材料(箱型梁骨架式)葉片在滿足強度、剛度、制作工藝等條件下，材料利用率較高，自重、成本相對較小。

3結(jié)論

本文在探索新型碳纖維復(fù)合材料葉片的基礎(chǔ)上，通過箱型梁骨架式葉片有限元分析，研究復(fù)合材料風(fēng)力機葉片在載荷作用的變形及其力學(xué)特性，進而揭示碳纖維復(fù)合材料的鋪層角度對變形和力學(xué)特性的影響規(guī)律，以其為新材料在風(fēng)力機葉片上的應(yīng)用提供技術(shù)支持和理論依據(jù)。

此外，為了進一步揭示箱型梁骨架式風(fēng)力機葉片的好處，對全復(fù)合材料式葉片也進行了分析，其葉片最大應(yīng)力比箱型梁式葉片要低，但是結(jié)構(gòu)造價和制作周期要長，因此綜合各方面建議選擇箱型梁骨架式風(fēng)力機葉片作為風(fēng)力機葉片的首選。

參考文獻：

[1]張承東. 風(fēng)力機葉片的動力學(xué)特性分析及分形特征研究[D]. 天津:天津工業(yè)大學(xué), 2007.

[2]董永棋. 國外風(fēng)力發(fā)電機FRP葉片近況[J]. 纖維復(fù)合材料, 2001.9 (3):46-49.

[3]田虎楠, 王作棠, 孫文磊,等. 極端高溫下風(fēng)力機葉片的性能分析[J]. 流體機械, 2013, 41(9): 37-39.

[4]周茜茜. MW級風(fēng)力機葉片氣動特性及流固禍合特性研究[D]. 長春:吉林大學(xué), 2015.

[5]孫麗萍, 王昊,丁嬌嬌. 風(fēng)力發(fā)電機葉片的氣動彈性及顫振研究綜述[J]. 液壓與氣動, 2012, 10:1-5.

[6]廖猜猜, 徐建中, 席光. 基于一種改進的 PSO 算法在風(fēng)力機葉片優(yōu)化中的應(yīng)用[J]. 工程熱物理學(xué)報, 2008, 29:773-776.

[7]楊從新, 張釗. 基于改進的粒子群算法的風(fēng)力機葉片優(yōu)化設(shè)計[J]. 新技術(shù)新工藝, 2013, 10:49-52.

[8]劉雄, 陳嚴, 葉枝全. 遺傳算法在風(fēng)力機風(fēng)輪葉片優(yōu)化設(shè)計中的應(yīng)用[J]. 太陽能學(xué)報, 2006, 2:180-185

[9]Lu X, Balendra R. Temperature-related errors on aerofoil section of turbine blade[J]. Journal of Materials Processing Technology, 2001, 115(2): 240-244.

[10] Cheng Lu, Zhang Liwen，MU Zhengjun,et al. 3D FEM simulation of multi-stage forging process of a gas turbine compressor blade [J]. Journal of Materials Processing Technology, 2009, 198(2):463-470.

Performance of box girder skeleton of wind turbine blade

JIANG Yan1, WANG Wei2,*, HUO Yin-Cai1, SUN Xiao-Ming1, HAN Mei1

(1.SchoolofHydraulicandElectricPower,HeilongjiangUniversity,Harbin150080,China; 2.CollegeofShipbuildingEngineering,HarbinEngineeringUniversity,Harbin150001,China)

Abstract：The blade is a key component of the wind turbine. The wind turbine blade materials will directly affect the performance of the wind turbine. In the present paper, a new carbon fiber composite blade was studied, through finite element analysis of frame beam type blade, the mechanical function and deformation under the load has been investigated. The ply angle of each layer has also been studied. These may provide technical and theoretical basis for improving chinese independent design and manufacturing capacity of long-lived large wind turbine blades.

Key words：blade； carbon fiber composite; performance; finite element analysis

DOI:10.13524/j.2095-008x.2016.01.005

收稿日期：2015-12-01

基金項目：黑龍江省教育廳科學(xué)技術(shù)研究項目(12531527)；國家自然科學(xué)基金資助項目(51209147)

作者簡介：姜艷(1966-)，女，黑龍江哈爾濱人，教授，碩士，研究方向：水工結(jié)構(gòu)，E-mail:jiangyan@hlju.edu.cn;*通訊作者:王瑋(1979-)，女，黑龍江大慶人，講師，博士，碩士研究生導(dǎo)師，研究方向：海洋結(jié)構(gòu)物力學(xué)性能，E-mail：mimal@126.com。

中圖分類號：TK83

文獻標志碼：A

文章編號：2095-008X(2016)01-0023-06