仿鸮翼型風(fēng)力機(jī)葉片氣動(dòng)彈性變形分析

摘要： 為抑制風(fēng)力機(jī)葉片的氣動(dòng)彈性變形，提出一種用于大型風(fēng)力機(jī)的仿鸮類翼型改型設(shè)計(jì)方案.通過(guò)單向流固耦合方法分析了仿生翼型對(duì)葉片氣動(dòng)彈性變形的影響，并通過(guò)模態(tài)分析和共振響應(yīng)闡述了仿生翼型抑制葉片氣動(dòng)彈性變形的機(jī)理.結(jié)果發(fā)現(xiàn)：與模型葉片相比，仿生翼型葉片氣動(dòng)彈性變形量降低了11.05%.仿生翼型使葉片壓力面應(yīng)力分布趨勢(shì)發(fā)生改變，使葉片壓力面的最大應(yīng)力值由葉中轉(zhuǎn)移至前緣.仿生翼型使葉片上下表面的壓力差提高，且吸力面的壓力值約為原來(lái)的2倍，使葉片提高了在揮舞方向的抗變形能力.同時(shí)，葉片的氣動(dòng)性能有所提高.仿生翼型使葉片的一階至六階固有頻率升高；仿生翼型使葉片在揮舞方向上的共振變形量減小了89.23%，葉片的共振速度和加速度幅值也降低了.故仿生翼型能夠有效抑制葉片的氣動(dòng)彈性變形，其設(shè)計(jì)方法及結(jié)論可為抑制大型風(fēng)力機(jī)葉片的氣動(dòng)彈性變形研究提供一種新思路.

關(guān)鍵詞： 風(fēng)力機(jī)葉片；氣動(dòng)彈性變形；仿鸮翼型；單向流固耦合；固有頻率

中圖分類號(hào)： S277.9 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)： 1674-8530（2024）04-0365-08

DOI：10.3969/j.issn.1674-8530.23.0066

陳坤，趙培堯，馮文慧，等.仿鸮翼型風(fēng)力機(jī)葉片氣動(dòng)彈性變形分析［J］.排灌機(jī)械工程學(xué)報(bào)，2024，42（4）：365-372.

CHEN Kun， ZHAO Peiyao， FENG Wenhui， et al. Analysis of aeroelastic deformation of bionic airfoil wind turbine blades［J］.Journal of drainage and irrigation machinery engineering（JDIME），2024，42（4）：365-372.（in Chinese）

Analysis of aeroelastic deformation of bionic airfoil wind turbine blades

CHEN Kun*， ZHAO Peiyao， FENG Wenhui， GUI Hongliang， HAO Zhenhua

（School of Mechanical Engineering， Xinjiang University， Urumqi， Xinjiang 830017， China）

Abstract： In order to suppress the aeroelastic deformation of wind turbine blades， a modified design of the bionic owl′s airfoil for large wind turbines was proposed. The effect of the bionic airfoil on the aeroe-lastic deformation of the blade was analyzed by the unidirectional fluid-structure coupling method， and the mechanism of the bionic airfoil suppressing the aeroelastic deformation of the blade was elaborated by modal analysis and resonance response. The results find that there is an 11.05% reduction in the aeroelastic deformation of the bionic airfoil blade compared to the model blade. The bionic airfoil makes the blade pressure surface stress distribution trend change， and makes the maximum stress value of the blade pressure surface transfer from the middle of the blade to the leading edge. The bionic airfoil makes the pressure difference between the upper and lower surface of the blade increase， and the pressure value of the suction surface is about twice the original， so that the blade in the flapping direction to improve the resistance to deformation. At the same time， the aerodynamic performance of the blades has been improved. The bionic airfoil pattern increases the first to sixth order intrinsic frequency of the blade， the bionic airfoil pattern causes an 89.23% reduction in the resonant deformation of the blade in the flapping direction， and also reduces the resonant velocity and acceleration amplitude of the blade. Therefore， the bionic airfoil can effectively suppress the aeroelastic deformation of the blade， and its design method and conclusion can provide a new idea for the research of suppressing aeroelastic deformation of large wind turbine blade.

Key words： wind turbine blade;aeroelastic deformation;bionic owl′s airfoil;unidirectional fluid-structure coupling;inherent frequency

當(dāng)今社會(huì)正面臨著能源短缺和環(huán)境惡化的問(wèn)題，急需大量的綠色清潔能源［1］.風(fēng)能作為綠色清潔能源的重要組成部分，它在實(shí)現(xiàn)全球能源轉(zhuǎn)型及碳中和目標(biāo)的過(guò)程中發(fā)揮著重要作用［2］.葉片是風(fēng)力機(jī)中唯一的捕風(fēng)部件，早期的葉片寬度大、長(zhǎng)度小，且運(yùn)行速度小.隨著風(fēng)電產(chǎn)業(yè)的發(fā)展，現(xiàn)在應(yīng)用的葉片不斷細(xì)長(zhǎng)化，且運(yùn)行速度加快，葉片的氣動(dòng)彈性變形問(wèn)題愈發(fā)顯著.同時(shí)，葉片斷裂的現(xiàn)象逐漸增多，導(dǎo)致了維修葉片的費(fèi)用急劇增大，因此，抑制風(fēng)力機(jī)葉片的氣動(dòng)彈性變形問(wèn)題逐漸成為研究熱點(diǎn)［3-4］.

單麗君等［5］采用多目標(biāo)遺傳算法對(duì)3 MW風(fēng)力機(jī)葉片進(jìn)行優(yōu)化，發(fā)現(xiàn)葉片的固有頻率得到提高，葉片展向撓度有所減小.張新虹等［6］探究了阻尼對(duì)葉片變形抑制的影響，發(fā)現(xiàn)存在最佳的阻尼值，當(dāng)超過(guò)該值后，葉片的變形反而會(huì)增大.楊瑞等［7］研究了阻尼層厚度對(duì)葉片變形的影響.SUN等［8］對(duì)比研究了穿孔阻尼葉片與標(biāo)準(zhǔn)葉片，發(fā)現(xiàn)葉片的結(jié)構(gòu)阻尼和結(jié)構(gòu)剛度都有所提高.周震霆等［9］研究了氣動(dòng)阻尼對(duì)葉片氣動(dòng)彈性振動(dòng)的影響，發(fā)現(xiàn)葉片達(dá)到穩(wěn)態(tài)的時(shí)間縮短了，葉片展向變形也得到了抑制.KHALID等［10］發(fā)現(xiàn)與傳統(tǒng)的閉孔結(jié)構(gòu)（如蜂窩）相比，增氧結(jié)構(gòu)的阻尼比更高，葉片的振動(dòng)幅度更小.CHANG等［11］研究發(fā)現(xiàn)Gurney襟翼結(jié)構(gòu)能夠減小葉片氣動(dòng)彈性振動(dòng)的振幅和振動(dòng)范圍，提高葉片的氣動(dòng)彈性穩(wěn)定性.LI等［12］提出采用單向電纜擺錘阻尼器來(lái)緩解大型風(fēng)力機(jī)葉片在擺動(dòng)方向的振動(dòng)響應(yīng).DAS等［13］將形狀記憶合金（SMA）沿葉片展向放置在內(nèi)部，用于葉片的加固.研究發(fā)現(xiàn)葉片在擺動(dòng)和揮舞方向的變形明顯減少.JIANG等［14］對(duì)三維角互鎖編織復(fù)合材料（3DAWC）懸臂梁的自由振動(dòng)進(jìn)行了研究，發(fā)現(xiàn)與傳統(tǒng)的雙軸經(jīng)編復(fù)合材料（BWKC）相比，3DAWC具有更好的動(dòng)態(tài)性能.

隨著仿生學(xué)的興起，學(xué)者們?cè)噲D從生物優(yōu)異的特征上尋求抑制風(fēng)力機(jī)葉片氣動(dòng)彈性變形的方法.MENG等［15］受竹子彎曲變形后可迅速恢復(fù)的啟發(fā)，研究了竹壁層結(jié)構(gòu)對(duì)葉片氣動(dòng)彈性變形的影響，發(fā)現(xiàn)仿竹壁層結(jié)構(gòu)能夠通過(guò)散熱的方式有效地消耗葉片產(chǎn)生振動(dòng)的能量，從而抑制葉片振動(dòng)的位移、速度和加速度.同年，孟杰等［16］根據(jù)毛竹竹節(jié)規(guī)律性的分布和隔板支撐竹竿的結(jié)構(gòu)設(shè)計(jì)了環(huán)形剪切腹板，發(fā)現(xiàn)環(huán)形剪切腹板可以顯著提高葉片根部的抗彎性能.張立等［17］研究了仿葉脈結(jié)構(gòu)的腹板，發(fā)現(xiàn)可以提高葉片的抗共振性能，且對(duì)葉片的振動(dòng)變形有抑制效果.賈文雋等［18］研究了海鷗翅翼的翼型和結(jié)構(gòu)，發(fā)現(xiàn)這2種結(jié)構(gòu)具有抑制機(jī)翼氣動(dòng)彈性的效果.根據(jù)文獻(xiàn)研究可知，目前基于仿生學(xué)對(duì)抑制葉片氣動(dòng)彈性變形的研究比較熱門，但是從葉片翼型設(shè)計(jì)角度研究抑制葉片氣動(dòng)彈性變形的文獻(xiàn)較少.

文中基于鸮翅膀展向截面翼型，提出一種適用于大型風(fēng)力機(jī)葉片設(shè)計(jì)的仿生翼型.將仿生翼型應(yīng)用于葉片上，通過(guò)數(shù)值計(jì)算，探究仿生翼型對(duì)風(fēng)力機(jī)葉片氣動(dòng)彈性變形的影響，并闡述仿生翼型抑制葉片氣動(dòng)彈性變形的機(jī)理.

1 葉片模型

1.1 5 MW風(fēng)力機(jī)葉片模型

選用美國(guó)可再生能源實(shí)驗(yàn)室研發(fā)的5 MW風(fēng)力機(jī)葉片作為研究對(duì)象.葉片由5種DU翼型和1種NACA翼型組成［19］，如圖1所示.5 MW風(fēng)力機(jī)的風(fēng)輪直徑為126 m，其中葉片長(zhǎng)61.5 m，輪轂直徑為3 m.風(fēng)力機(jī)的額定風(fēng)速為11.4 m/s，額定轉(zhuǎn)速為12.1 r/min，葉尖速比為7.在葉片模型中，葉片展向截面翼型的50%弦長(zhǎng)處為5 mm厚度的平板腹板.

1.2 相似性設(shè)計(jì)

由于5 MW風(fēng)力機(jī)葉片尺寸過(guò)大，在能精確地模擬風(fēng)力機(jī)葉片流場(chǎng)的要求下，其所需的計(jì)算資源過(guò)于龐大.同時(shí)，5 MW風(fēng)力機(jī)葉片的試驗(yàn)研究也難以開(kāi)展.因此，基于相似原理對(duì)5 MW風(fēng)力機(jī)葉片進(jìn)行相似化設(shè)計(jì).

根據(jù)相似原理，風(fēng)力機(jī)縮比模型和原型相似應(yīng)同時(shí)滿足幾何相似、運(yùn)動(dòng)相似和動(dòng)力相似.動(dòng)力相似則需要模型和原型在流動(dòng)中同時(shí)滿足雷諾數(shù)相似、歐拉數(shù)相似、馬赫數(shù)相似、弗勞德數(shù)相似、斯特勞哈數(shù)相似等.根據(jù)楊洪磊等［20］對(duì)風(fēng)力機(jī)葉片相似的相關(guān)研究，發(fā)現(xiàn)原型與模型之間同時(shí)滿足幾何相似、運(yùn)動(dòng)相似和動(dòng)力相似是困難的，研究目標(biāo)不同，可以適當(dāng)調(diào)整相似參數(shù)，對(duì)控制參數(shù)有所側(cè)重，確保主要因素的相似性.文中優(yōu)先滿足幾何相似、運(yùn)動(dòng)相似和雷諾數(shù)相似準(zhǔn)則，確保模型和原型在相似的流場(chǎng)中，功率和動(dòng)量矩系數(shù)均相同.風(fēng)力機(jī)模型設(shè)計(jì)是嚴(yán)格按照相似律的規(guī)定進(jìn)行的.

風(fēng)力機(jī)葉片長(zhǎng)度L、葉輪轉(zhuǎn)速Ω、轉(zhuǎn)矩M、功率P的相似關(guān)系分別為L(zhǎng)1L2=50，Ω1Ω2=50-1，M1M2=503，P1P2=502，

其中比例因子采用50.根據(jù)相似率縮比得到模型風(fēng)力機(jī)的相關(guān)參數(shù)：風(fēng)輪直徑為2.52 m，葉片長(zhǎng)1.23 m，額定功率為2×10-3 MW，額定轉(zhuǎn)速為605 r/min，額定風(fēng)速和葉尖速比與原型相同.

1.3 仿生翼型設(shè)計(jì)

首先，采用逆向建模得到長(zhǎng)耳鸮翅膀模型，并在建模軟件SolidWorks中進(jìn)行切割處理，獲得了鸮翅膀展向每10%的截面翼型，并進(jìn)行多項(xiàng)式擬合和光順化處理.

通過(guò)對(duì)比鸮翅膀展向每10%的截面翼型，發(fā)現(xiàn)鸮翼展向50%截面翼型具有較高的升阻比，適合作為葉片設(shè)計(jì)的基礎(chǔ)翼型.但不能直接應(yīng)用于模型葉片中，模型葉片的組成翼型的相對(duì)厚度大，能保證葉片運(yùn)行時(shí)的剛度和不易發(fā)生過(guò)度彎曲現(xiàn)象.

為保證葉片具有較為安全的剛度，需將仿生翼型與原型葉片上的NACA64翼型的相對(duì)厚度保持一致.首先，將2種翼型的弦長(zhǎng)保持為0.1 m且水平放置.接著，通過(guò)觀察兩翼型上下弦線的結(jié)構(gòu)特點(diǎn)，發(fā)現(xiàn)翼型上弦線的相對(duì)厚度相同，但最大厚度位置不同；鸮翼50%截面翼型的下弦線為上凹，NACA64翼型下弦線為下凸，兩者的相對(duì)厚度相差達(dá)到最大.因此，保留鸮翼50%截面翼型上弦線和NACA64翼型下弦線，采用坐標(biāo)組合的方式，通過(guò)多項(xiàng)式擬合和光順化處理后，建立了仿生翼型，如圖2所示，圖中橫坐標(biāo)X為翼型橫向坐標(biāo)，縱坐標(biāo)Y為翼型縱向坐標(biāo).

將仿生翼型應(yīng)用在模型葉片展向70%～100%處，保持重構(gòu)后的截面翼型和原截面翼型的弦長(zhǎng)和槳距角一致，建立了仿生葉片模型，如圖3所示.

2 數(shù)值模擬

2.1 流場(chǎng)數(shù)值模擬

基于k-ω SST湍流模型對(duì)模型風(fēng)力機(jī)葉片和仿生葉片進(jìn)行了流場(chǎng)模擬.采用圓柱體計(jì)算域，計(jì)算域由靜止域、交界面和旋轉(zhuǎn)域組成，交界面是靜止域和旋轉(zhuǎn)域共有的面，用于傳遞信息.計(jì)算域采用非結(jié)構(gòu)網(wǎng)格劃分，并對(duì)風(fēng)力機(jī)表面及附近網(wǎng)格進(jìn)行了加密處理.基于k-ω SST湍流模型對(duì)近壁面y+＜1的要求，確定了葉片表面邊界層的第1層高度為0.316×10-4 m.通過(guò)網(wǎng)格無(wú)關(guān)性驗(yàn)證，如圖4所示，n為網(wǎng)格的總體數(shù)量，M為在對(duì)應(yīng)網(wǎng)格數(shù)量下測(cè)得的風(fēng)力機(jī)轉(zhuǎn)矩.最終確定了模型風(fēng)力機(jī)流場(chǎng)的網(wǎng)格數(shù)量為1 172萬(wàn)，仿生風(fēng)力機(jī)流場(chǎng)的網(wǎng)格數(shù)量為1 192萬(wàn).圖5為流場(chǎng)中葉片周圍網(wǎng)格的細(xì)節(jié).

模型風(fēng)力機(jī)葉片與仿生葉片的邊界條件相同，均采用速度入口和壓力出口［21］.入口的速度設(shè)為風(fēng)力機(jī)額定風(fēng)速11.4 m/s，出口的壓力設(shè)為大氣壓.采用滑移網(wǎng)格解決了葉片的旋轉(zhuǎn)問(wèn)題，風(fēng)力機(jī)轉(zhuǎn)速為605 r/min.

穩(wěn)態(tài)流場(chǎng)采用k-ω SST湍流模型，算法采用SIMPLEC算法，離散格式為二階迎風(fēng)格式，求解過(guò)程中，殘差設(shè)為1×10-4.流場(chǎng)模擬過(guò)程中控制方程采用納維-斯托克斯（N-S）方程，其表達(dá)式為

Qt+（Fi+Fv）x+（Gi+Gv）y+（Hi+Hv）z=0，（1）

式中：守恒變量向量Q=（ρ，ρu，ρv，ρw，ρE）T；Fi，Gi，Hi為對(duì)流項(xiàng)通量向量；Fv，Gv，Hv為黏性項(xiàng)通量向量.

2.2 結(jié)構(gòu)場(chǎng)數(shù)值模擬

風(fēng)力機(jī)葉片大多采用玻璃纖維和碳纖維等增強(qiáng)材料，其中E-玻璃纖維在強(qiáng)度、延展性、絕緣性、耐熱性和耐濕性方面具有優(yōu)良的性能，使用較為普遍.因此，文中設(shè)置葉片的材料為E-玻璃纖維，材料密度為2.54 g/cm3，楊氏模量為7.3×104 MPa，泊松比為0.22，材料的許用應(yīng)力為686 MPa.

風(fēng)力機(jī)葉片采用非結(jié)構(gòu)網(wǎng)格進(jìn)行劃分，對(duì)葉片前緣和尾緣處進(jìn)行細(xì)化網(wǎng)格.對(duì)風(fēng)力機(jī)輪轂處進(jìn)行固定約束，并設(shè)置風(fēng)力機(jī)轉(zhuǎn)速為605 r/min.將流場(chǎng)中葉片表面的壓力信息傳遞至結(jié)構(gòu)場(chǎng)葉片表面進(jìn)行分析.結(jié)構(gòu)場(chǎng)模擬過(guò)程中控制方程是基于牛頓第二定律獲得的，不考慮能量方程，其表達(dá)式為

式中：ρs為結(jié)構(gòu)密度；d··s為結(jié)構(gòu)域的加速度矢量；σs為柯西應(yīng)力張量； fs為體積力矢量.

在流固耦合交界面處，流場(chǎng)與結(jié)構(gòu)場(chǎng)的應(yīng)力應(yīng)變等信息應(yīng)對(duì)應(yīng)相等，其方程為

τf·nf=τs·ns，

df=ds，（3）

式中：τf，τs分別為流場(chǎng)和結(jié)構(gòu)場(chǎng)的應(yīng)力；nf，ns分別為流場(chǎng)和結(jié)構(gòu)場(chǎng)的單位向量；df，ds分別為流場(chǎng)和結(jié)構(gòu)場(chǎng)的位移.

3 結(jié)果分析

3.1 數(shù)值計(jì)算模型有效性驗(yàn)證

通過(guò)在流場(chǎng)中監(jiān)測(cè)風(fēng)力機(jī)的轉(zhuǎn)矩，選取風(fēng)力機(jī)的輸出功率和功率系數(shù)作為驗(yàn)證數(shù)值計(jì)算模型有效性的流場(chǎng)評(píng)價(jià)指標(biāo).

通過(guò)流場(chǎng)計(jì)算，與1∶50相似化處理的模型葉片輸出功率值相比，數(shù)值計(jì)算模型的功率損失在1.50%以內(nèi).將數(shù)據(jù)與美國(guó)可再生能源重點(diǎn)實(shí)驗(yàn)室公布的5 MW葉片功率系數(shù)數(shù)據(jù)進(jìn)行對(duì)比，數(shù)值計(jì)算模型的功率系數(shù)誤差e也在1.50%以內(nèi)，如表1所示.從風(fēng)力機(jī)的輸出功率P和功率系數(shù)CT上驗(yàn)證了數(shù)值計(jì)算模型的有效性.

同時(shí)，以KYLE等［22］對(duì)額定工況下NREL 5 MW葉片展向截面壓力系數(shù)的研究為參考標(biāo)準(zhǔn)，進(jìn)行數(shù)值計(jì)算模型展向截面壓力系數(shù)的對(duì)比驗(yàn)證.圖6為文獻(xiàn)中和文中數(shù)值計(jì)算獲得的葉片展向80%處截面壓力系數(shù)分布曲線.X/C為相對(duì)弦向位置，X為翼型表面點(diǎn)的弦向坐標(biāo)，C為翼型弦長(zhǎng).從圖中可以發(fā)現(xiàn)，數(shù)值計(jì)算的結(jié)果與文獻(xiàn)中壓力系數(shù)分布之間的誤差較小，具有良好的擬合度.從葉片展向截面壓力系數(shù)分布上驗(yàn)證了數(shù)值計(jì)算模型的有效性.

3.2 氣動(dòng)彈性變形分析

圖7為模型葉片和仿生葉片在額定工況下的變形量及其變化幅度隨葉片展向的分布，圖中r為點(diǎn)到輪轂中心的位移，R為風(fēng)輪半徑，r/R為點(diǎn)在葉片展向的相對(duì)位置，y為葉片展向的撓度，SymbolDA@y為仿生葉片與模型葉片展向的氣動(dòng)彈性變形量差值.由圖7可以看出，靠近葉根處，仿生葉片與模型葉片的變形相差較小.總體上，仿生葉片的變形量比模型葉片小，且與模型葉片的變形量差值沿展向逐漸增大.與模型葉片相比，仿生葉片葉尖處的氣動(dòng)彈性變形量被抑制了11.05%.

圖8為模型葉片和仿生葉片不同振動(dòng)形式下的變形量及其變化隨葉片展向的分布.可看出揮舞振動(dòng)是葉片變形的主要形式，葉片在擺動(dòng)方向的變形量約為揮舞方向的14.70%.與模型葉片相比，仿生葉片葉尖處揮舞方向的變形量被抑制了11.18%，擺振方向變形量被抑制了9.30%.

3.3 葉片應(yīng)力應(yīng)變分析

圖9為模型葉片和仿生葉片在額定風(fēng)速下壓力面的壓力分布.從圖中可以看出，模型葉片壓力面上葉中部的應(yīng)力值σ最大，葉尖和葉根處應(yīng)力值最小，總體上呈現(xiàn)為由葉中向葉片兩端輻射狀降低的應(yīng)力分布趨勢(shì).模型葉片葉中部承受的力最大，因此，葉片在葉中處發(fā)生斷裂或失效的概率更大.與模型葉片相比，仿生葉片表面的應(yīng)力分布發(fā)生了改變.靠近葉根的葉片前緣處應(yīng)力值最大，在葉片的翼型段，呈現(xiàn)為由葉片前緣至尾緣逐漸減小的應(yīng)力分布趨勢(shì).究其原因主要是葉片展向70%Symbol～A@100%處的翼型增加了鸮翼的結(jié)構(gòu)特征.總體上，模型葉片和仿生葉片產(chǎn)生形變后的最大應(yīng)力值均遠(yuǎn)小于材料的許用應(yīng)力，在額定風(fēng)況下可安全運(yùn)行.

圖10為模型葉片和仿生葉片在額定風(fēng)速下壓力面的應(yīng)變?chǔ)欧植?從圖中可看出風(fēng)力機(jī)葉片表面的應(yīng)變分布趨勢(shì)與其應(yīng)力分布趨勢(shì)相同.模型葉片壓力面呈現(xiàn)為葉中向兩端輻射狀降低的應(yīng)變分布趨勢(shì)，仿生葉片壓力面在翼型段呈現(xiàn)為由葉片前緣至尾緣不斷降低的應(yīng)變分布趨勢(shì).

3.4 葉片表面壓力分布

圖11和圖12分別為模型葉片和仿生葉片壓力面和吸力面的壓力p分布云圖.從圖中可以發(fā)現(xiàn)，模型葉片和仿生葉片在壓力面上的壓力分布是大致相同的.

同時(shí)，與模型葉片吸力面壓力分布相比，仿生翼型改變了葉片仿生段的壓力分布.總體上，仿生葉片吸力面所受到的負(fù)壓值低于模型葉片.

3.5 葉片展向截面壓力系數(shù)分析

圖13為模型葉片和仿生翼型展向80%和90%截面處的壓力系數(shù)Cp的分布情況.

由圖13可以看出，在80%截面處，仿生葉片壓力面在0.05＜X/C＜0.25和X/C＞0.75的范圍內(nèi)的上下表面壓力差明顯高于模型葉片.仿生翼型較NACA64翼型更能提高葉片的氣動(dòng)性能.在額定風(fēng)況下，翼型壓力面的正壓及吸力面負(fù)壓的絕對(duì)值在靠近前緣時(shí)逐漸增大，前緣處存在最大壓差，在尾緣處壓力最小.仿生葉片與模型葉片壓力面受到的壓力是幾乎一致的，但仿生葉片在吸力面上的壓力總體上是高于模型葉片的.因此，當(dāng)葉片發(fā)生揮舞方向的變形時(shí)，仿生葉片受到的阻力作用是高于模型葉片的，其變形量也會(huì)因此減小.

3.6 模態(tài)振型和固有頻率分析

圖14和圖15分別為模型葉片和仿生葉片的前六階模態(tài).從圖中可看出，葉片的一階振型為揮舞振動(dòng)變形，二階三階振型為揮舞振動(dòng)和擺振耦合形式的變形.四階振型為擺振變形，五階六階振型為擺振和揮舞振動(dòng)耦合形式的變形.從葉片的變形量s來(lái)看，仿生葉片的一階模態(tài)和六階模態(tài)的變形量相較于模型葉片來(lái)說(shuō)有所減小，仿生葉片的六階固有頻率均大于模型葉片.這是因?yàn)椋律硇透淖兞巳~片的形狀，改變了葉片的剛度分布.同時(shí)，仿生翼型使葉片的固有頻率提高了，可在一定程度上防止葉片發(fā)生共振.

3.7 共振響應(yīng)分析

圖16為模型葉片和仿生葉片在葉尖角點(diǎn)揮舞方向的振動(dòng)位移y1隨頻率f變化的曲線.

由圖16可看出，模型葉片和仿生葉片均在34 Hz左右出現(xiàn)了位移急劇變化.從葉片的模態(tài)分析中可知，34 Hz左右對(duì)應(yīng)的是模型葉片和仿生葉片的前三階模態(tài)頻率.在揮舞共振頻率下，與模型葉片相比，仿生葉片葉尖角點(diǎn)的振動(dòng)位移減少了89.23%.

3.8 輸出功率分析

在額定工況中，對(duì)模型風(fēng)力機(jī)和仿生風(fēng)力機(jī)的轉(zhuǎn)矩進(jìn)行監(jiān)測(cè).當(dāng)風(fēng)力機(jī)運(yùn)行穩(wěn)定后，模型葉片的轉(zhuǎn)矩值為31.107 9 N·m，仿生葉片的轉(zhuǎn)矩值為32.135 0 N·m.由轉(zhuǎn)矩值進(jìn)而獲得模型葉片的輸出功率為1 970.86 W，仿生葉片的輸出功率為2 035.93 W.與模型葉片相比，仿生翼型可使風(fēng)力機(jī)的輸出功率提高3.3%.可見(jiàn)，仿生翼型對(duì)風(fēng)力機(jī)的氣動(dòng)性能具有提升的效果.

4 結(jié) 論

1） 通過(guò)對(duì)風(fēng)力機(jī)葉片氣動(dòng)彈性變形的分析可知，在額定工況下，仿生翼型可有效減少風(fēng)力機(jī)葉片的氣動(dòng)彈性位移量約11.05%，其中在揮舞方向的變形量可有效降低11.18%.

2） 通過(guò)對(duì)葉片表面應(yīng)力應(yīng)變和截面壓力的分析可知，仿生翼型使葉片壓力面由葉中向兩端輻射狀降低的應(yīng)力分布轉(zhuǎn)變?yōu)橛汕熬壷廖簿壊粩嘟档偷膽?yīng)力分布.仿生翼型的上下表面壓力差高于NACA64翼型，吸力面上總體的壓力也高于NACA64翼型.

3） 通過(guò)對(duì)葉片固有頻率和共振響應(yīng)的分析可知，揮舞振動(dòng)是葉片主要的振動(dòng)形式，其次為擺振.仿生翼型對(duì)葉片的固有頻率有所提高.仿生翼型可有效降低葉片發(fā)生一階揮舞共振的位移、速度和加速度.

參考文獻(xiàn)（References）

［1］ 黃顯峰，周引航，張啟凡，等.基于云模型的水光互補(bǔ)清潔能源基地容量配置方案優(yōu)選［J］.水利水電科技進(jìn)展，2024，44（1）：44-51.

ＨUANG Xianfeng， ZHOU Yinhang， ZHANG Qifan， et al. Capacity allocation scheme optimization based on cloud model for a hydro-photovoltaic complementary clean energy base［J］. Advances in science and techno-logy of water resources， 2024，44（1）：44-51.（in Chinese）

［2］ 陳中亞，吳煒，何俊尉，等.風(fēng)電機(jī)組法蘭用摩擦墊片的摩擦系數(shù)實(shí)驗(yàn)研究［J］. 機(jī)電工程，2022，39（9）： 1319-1324.

CHEN Zhongya，WU Wei，HE Junwei，et al. Experi-mental study on friction coefficient of friction gaskets for the flanges of wind turbine［J］. Journal of mechanical ＆ electrical engineering， 2022，39（9）： 1319-1324. （in Chinese）

［3］ 楊從新， 張宇婷， 岳念西.IDDES 方法模擬風(fēng)切變下大型水平軸風(fēng)力機(jī)流場(chǎng)特性［J］. 西華大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2024， 43 （1）： 78-86.

YANG Congxin，ZHANG Yuting，YUE Nianxi. The flow field characteristics of large horizontal axis wind turbine in shear inflow based on IDDES method［J］. Journal of Xihua University（natural science edition）， 2024， 43 （1）： 78-86. （in Chinese）

［4］ 楊從新，劉文杰，李壽圖，等.山區(qū)湍流特征及其對(duì)風(fēng)力機(jī)功率的影響［J］. 西華大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2023， 42（2）：28-38.

YANG Congxin， LIU Wenjie， LI Shoutu， et al. Characteristics of turbulence in mountainous areas and its effect on the power of wind turbines［J］. Journal of Xihua University（natural science edition），2024， 43（2）： 28-38. （in Chinese）

［5］ 單麗君，柳敬元. 3 MW風(fēng)力機(jī)葉片振動(dòng)特性分析［J］. 可再生能源，2019，37（7）：1067-1072.

SHAN Lijun， LIU Jingyuan. Vibration analysis of 3 MW wind turbine blade ［J］. Renewable energy resources， 2019， 37（7）： 1067-1072. （in Chinese）

［6］ 張新虹，潘宏剛，洪志同，等. 阻尼葉片振動(dòng)特性試驗(yàn)研究［J］. 沈陽(yáng)工程學(xué)院學(xué)報(bào)（自然科學(xué)版），2020，16（2）：21-24.

ZHANG Xinhong， PAN Honggang， HONG Zhitong， et al. Experimental study on vibration characteristics of damped blades［J］. Journal of Shenyang Institute of Engineering（natural science）， 2020， 16（2）：21-24. （in Chinese）

［7］ 楊瑞，楊偉，陳志龍. 大型風(fēng)力機(jī)葉片阻尼層厚度分析與抑顫研究［J］. 蘭州理工大學(xué)學(xué)報(bào)，2021，47（4）：45-51.

YANG Rui， YANG Wei， CHEN Zhilong. Analysis of thickness of damping layer of large wind turbine blades and study on flutter suppression ［J］. Journal of Lanzhou University of Technology， 2021， 47（4）：45-51. （in Chinese）

［8］ SUN D G， GUO J J， SONG Y， et al. Flutter stability analysis of a perforated damping blade for large wind turbines［J］." Journal of sandwich structures and materials， 2019， 21（3）： 973-989.

［9］ 周震霆，賀星，劉永葆. 氣動(dòng)阻尼對(duì)裂紋葉片振動(dòng)特性的影響研究［J］. 燃?xì)廨啓C(jī)技術(shù)，2022，35（3）：42-48.

ZHOU Zhenting， HE Xing， LIU Yongbao. Influence of aerodynamic damping on vibration characteristics of cracked blade ［J］. Gas turbine technology，2022，35（3）：42-48. （in Chinese）

［10］ KHALID S A， KHAN A M， SHAH O R. A numerical study into the use of auxectic structures for structural damping in composite sandwich core panels for wind turbine blades［J］. Journal of energy resources technology， 2022，144（3）： 031301.

［11］ CHANG L， YU Y J，LIU T R，et al. Aeroelastic vibration analysis of wind turbine blade with Gurney flap［J］. Journal of mechanical engineering science， 2021， 235（20）： 4913-4925.

［12］ LI W R， QIAO L， LI G， et al. Vibration control of large wind turbine blades with unidirectional cable pendulum damper［J］. International journal of structural stability and dynamics， 2022，23（7）： 2350082.

［13］ DAS S， SAJEER M M， CHAKRABORTY A， et al. Shape memory alloy-based centrifugal stiffening for response reduction of horizontal axis wind turbine blade［J］. Structural control and health monitoring， 2021，28（3）：e2669.

［14］ JIANG X， JIANG Y， LU S， et al. Free vibration of three-dimensional angle-interlock woven composite multilayered cantilever beam［J］. Journal of vibration engineering amp; technologies， 2023，11：3387-3398.

［15］ MENG J， SUN D G. Research on vibration suppression of wind turbine blade with a multi-layer porous damping structure based on bamboo wall microstructure［J］. Journal of the Brazilian Society of Mechanical Sciences and Engineering， 2021， 43：153.

［16］ 孟杰，張喜清，孫大剛.帶有仿生環(huán)形剪切腹板風(fēng)機(jī)葉片抗彎強(qiáng)度［J］.科學(xué)技術(shù)與工程，2021，21（31）：13354-13360.

MENG Jie， ZHANG Xiqing， SUN Dagang. The effects of ring shear webs on the bending strength of wind turbine blade ［J］. Science technology and engineering，2021，21（31）：13354-13360. （in Chinese）

［17］ 張立，劉宇航，李春，等. 基于仿生設(shè)計(jì)的風(fēng)力機(jī)葉片腹板力學(xué)性能分析［J］. 熱能動(dòng)力工程， 2019， 34（9）：141-147.

ZHANG Li， LIU Yuhang， LI Chun， et al. Mechanical performance analysis of wind turbine blade web based on bionic design ［J］. Journal of engineering for thermal energy and power， 2019， 34（9）：141-147. （in Chinese）

［18］ 賈文雋，華欣，史繼拓，等.基于流固耦合的仿生海鷗機(jī)翼氣動(dòng)彈性分析［J］.兵器裝備工程學(xué)報(bào)，2020，41（7）：106-111.

JIA Wenjun， HUA Xin， SHI Jituo， et al. Aeroelastic analysis of bionic seagull wing based on fluid-structure coupling ［J］. Journal of ordnance equipment enginee-ring，2020，41（7）：106-111.（in Chinese）

［19］ 張浩. L小翼對(duì)5 MW水平軸風(fēng)力機(jī)氣動(dòng)性能的影響［D］.蘭州：蘭州理工大學(xué)，2019.

［20］ 楊洪磊，張明明，徐建中. 海上風(fēng)電葉片智能控制實(shí)驗(yàn)系統(tǒng)研制及控制效果分析［J］. 工程熱物理學(xué)報(bào)，2021，42（3）：611-618.

YANG Honglei，ZHANG Mingming，XU Jianzhong.Development of smart control experiment system for offshore wind turbine blades and analysis of control effect［J］. Journal of engineering thermophysics， 2021， 42（3）： 611-618.（in Chinese）

［21］ 姚偉偉. 基于鸮翼特征的仿生耦合風(fēng)力機(jī)葉片氣動(dòng)性能研究［D］. 烏魯木齊：新疆大學(xué)， 2021.

［22］ KYLE R， WANG F， FORBES B. The effect of a leading edge erosion shield on the aerodynamic performance of a wind turbine blade［J］.Wind energy，2020，23（4）：953-966.

（責(zé)任編輯 朱漪云）

收稿日期： 2023-04-17； 修回日期： 2023-08-17； 網(wǎng)絡(luò)出版時(shí)間： 2024-04-11

網(wǎng)絡(luò)出版地址： https：//link.cnki.net/urlid/32.1814.TH.20240408.1416.002

基金項(xiàng)目： 新疆維吾爾自治區(qū)自然科學(xué)基金資助項(xiàng)目（2022D01C33）；新疆維吾爾自治區(qū)研究生科研創(chuàng)新項(xiàng)目（XJ2022G049）；國(guó)家自然科學(xué)基金資助項(xiàng)目（51765062）

第一作者簡(jiǎn)介： 陳坤（1984—），女，河南太康人，副教授，博士（通信作者，chenkun_34@163.com），主要從事流體機(jī)械仿生技術(shù)的應(yīng)用研究.

第二作者簡(jiǎn)介： 趙培堯（1996—），男，河北邯鄲人，碩士研究生（zhaopeiyao_91@163.com），主要從事風(fēng)力機(jī)結(jié)構(gòu)仿生的應(yīng)用研究.