風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法

2016-12-21  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

隨著風(fēng)電機(jī)組單機(jī)容量的不斷增大,風(fēng)電機(jī)組葉片尺寸也越來越長。對(duì)于超長風(fēng)電機(jī)組葉片的氣動(dòng)計(jì)算及結(jié)構(gòu)計(jì)算,常規(guī)的葉片剛性體假設(shè)所引起的誤差越來越大,因此必須考慮葉片在風(fēng)載條件下結(jié)構(gòu)變形對(duì)氣動(dòng)力的影響。流固耦合(fluid-structure interaction, FSI)是研究風(fēng)電機(jī)組葉片氣彈特性及結(jié)構(gòu)特性的重要手段,國內(nèi)外通過流固耦合對(duì)葉片的結(jié)構(gòu)響應(yīng)、流場變化、葉片載荷分布等特性進(jìn)行了一系列研究。流固耦合按照求解方法可分為完全耦合和交替求解。李媛等以MPCCI為數(shù)據(jù)交換平臺(tái)連接商用CFD,CSD軟件,模擬考慮風(fēng)切變時(shí)輪轂高度額定風(fēng)速下的風(fēng)輪全三維流固耦合,分析了風(fēng)剪切對(duì)葉片繞流、變形以及耦合作用對(duì)載荷分布的影響規(guī)律。胡丹梅等采用k-SST紊流模型和滑移網(wǎng)格技術(shù),對(duì)海上風(fēng)電機(jī)組葉片進(jìn)行了流固藕合計(jì)算分析,發(fā)現(xiàn)葉片流固耦合作用使葉片氣動(dòng)攻角、扭矩增大。Zahle等采用不可壓縮重疊網(wǎng)格技術(shù),分析了水平軸風(fēng)電機(jī)組風(fēng)輪與塔架間的耦合作用。Hsu等通過對(duì)5 MW三葉片水平軸風(fēng)電機(jī)組風(fēng)輪在有無塔架情況下進(jìn)行了流固耦合研究,發(fā)現(xiàn)葉片轉(zhuǎn)過塔架時(shí)由于塔影效應(yīng)使單個(gè)葉片氣動(dòng)扭矩降低10%-12%。 Zhang等通過對(duì)比研究,發(fā)現(xiàn)風(fēng)切變對(duì)風(fēng)電機(jī)組葉片變形及應(yīng)力的影響明顯大于流固耦合所產(chǎn)生的作用。任年鑫等基于三維黏性不可壓縮Navier-Stokes方程和重整化群k-ε(RNG)湍流模型,數(shù)值模擬美國可再生能源實(shí)驗(yàn)室(NREL)5 MW海上風(fēng)電機(jī)組的氣動(dòng)性能,并研究了浮式平臺(tái)不同運(yùn)動(dòng)幅值及運(yùn)動(dòng)周期對(duì)風(fēng)電機(jī)組葉片氣動(dòng)性能的影響規(guī)律。陳海萍等選用k-ε湍流模型,將風(fēng)作為黏性不卡壓縮流體,計(jì)算出流固耦合作用下風(fēng)電機(jī)組葉片的應(yīng)變分布。李德源等采用計(jì)算多體系統(tǒng)動(dòng)力學(xué)理論,基于R-W(Roberson-Wittenburg)建模方法,建立柔性葉片的氣彈耦合方程,結(jié)合相應(yīng)的風(fēng)電機(jī)組氣動(dòng)載荷分析模塊,實(shí)現(xiàn)了葉片的氣彈耦合分析。王旭東等將葉片的啟動(dòng)與結(jié)構(gòu)力學(xué)模型進(jìn)行耦合,研究葉片變形對(duì)來流的影響,提出一種旋轉(zhuǎn)葉片結(jié)構(gòu)性能分析方法。Dong等應(yīng)用CFD-CSD禍合方法對(duì)風(fēng)電機(jī)組葉片的載荷及氣彈特性進(jìn)行研究,發(fā)現(xiàn)葉片氣彈變形顯著降低氣動(dòng)載荷。MEXICO (model experiment incontrolled conditions)實(shí)驗(yàn)是由歐盟資助的一項(xiàng)大型風(fēng)電機(jī)組實(shí)驗(yàn)。梁明軒通過對(duì)葉片模態(tài)分析得到葉片各階模態(tài),發(fā)現(xiàn)葉片扭轉(zhuǎn)、揮舞、擺振三者之間的耦合,并對(duì)葉片氣動(dòng)彈性問題進(jìn)行研究,建立葉片顫振方程,確定了葉片顫振發(fā)散條件。楊華等應(yīng)用風(fēng)電機(jī)組旋轉(zhuǎn)葉片表面壓力的三維PIV(particle image velocimetry)流場測量方法,測試得到無偏航工況下葉片上5個(gè)測試斷面的表面壓力系數(shù)分布,并計(jì)算了該工況下葉輪所受的軸向力和切向力。潘旭運(yùn)用Ansys、CFX基于Ansys Workbench軟件平臺(tái)對(duì)某MW級(jí)水平軸風(fēng)電機(jī)組葉片進(jìn)行流固耦合數(shù)值模擬,分析了藕合過程中葉片的變形及應(yīng)力變化。陸洋等利用Newmark數(shù)值積分法獲得葉片氣彈響應(yīng)的穩(wěn)態(tài)周期解,分別以NREL Phase VI非定常空氣動(dòng)力學(xué)實(shí)驗(yàn)及其公開的1.5 MW風(fēng)電機(jī)組葉片為算例計(jì)算了有/無預(yù)彎葉片的氣彈響應(yīng)。以上研究均通過建立風(fēng)輪旋轉(zhuǎn)的計(jì)算模型對(duì)風(fēng)輪流固禍合進(jìn)行研究,網(wǎng)格劃分困難,涉及到移動(dòng)域與靜止域網(wǎng)格匹配等復(fù)雜過程,計(jì)算成本高。為此,本文采用風(fēng)切變形式模擬風(fēng)輪旋轉(zhuǎn)及來流風(fēng)速的綜合效應(yīng),研究風(fēng)電機(jī)組葉片流固耦合數(shù)值模擬的一種新方法。


 1流體計(jì)算模型


1.1流體域求解方法


流場中的控制體既不能產(chǎn)生質(zhì)量也不能消滅質(zhì)量,因此它滿足質(zhì)量守恒定律。該定律在流場中的數(shù)學(xué)表達(dá)式為連續(xù)性方程

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys圖片1

式中,p,t,V分別為流體密度、時(shí)間和流體速度。

動(dòng)量守恒定律可表述為:控制體中流體的動(dòng)量對(duì)時(shí)間的變化率等于外界作用在該微元體上的各種力之和。其在流場中的數(shù)學(xué)表達(dá)式為運(yùn)動(dòng)方程

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys圖片2

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片3


能量守恒定律是流動(dòng)系統(tǒng)必須滿足的基本定理,在不考慮溫度變化條件下,其數(shù)學(xué)表達(dá)式為


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片4


式中,F為質(zhì)量力,P 為表面壓強(qiáng)。


1.2葉素速度


風(fēng)電機(jī)組風(fēng)輪運(yùn)行過程中實(shí)際參考風(fēng)速為來流風(fēng)速和旋轉(zhuǎn)風(fēng)速的共同作用,通過改變?nèi)~片不同葉展處的風(fēng)速可實(shí)現(xiàn)正常來流風(fēng)速和旋轉(zhuǎn)風(fēng)速的疊加效果。不同葉展處的參考風(fēng)速Vref的計(jì)算公式為

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片5

式中,Vo為自然來流風(fēng)速;w 為風(fēng)輪旋轉(zhuǎn)角速度;r 為葉展位置的旋轉(zhuǎn)半徑;a為軸向誘導(dǎo)系數(shù);a’為切向誘導(dǎo)系數(shù)。

葉素速度矢量的三角形關(guān)系如圖1所示,圖中α、β、?分


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片6
風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys分析案例圖片7


 2結(jié)構(gòu)計(jì)算模型

2.1結(jié)構(gòu)計(jì)算方法


風(fēng)電機(jī)組葉片為多自由度有阻尼的非線性結(jié)構(gòu),其結(jié)構(gòu)動(dòng)力學(xué)運(yùn)動(dòng)方程為

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys分析案例圖片8

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys結(jié)果圖片9

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys結(jié)果圖片10


2.2葉片結(jié)構(gòu)模型


為了用靜止葉片模擬葉片旋轉(zhuǎn)效應(yīng),只通過將葉素處的速度按照速度矢量關(guān)系修正還不能達(dá)到應(yīng)有效果,需要對(duì)葉片不同截面處翼型的氣動(dòng)扭角進(jìn)行修正。如圖1所示,以旋轉(zhuǎn)導(dǎo)致的風(fēng)速方向作為參考方向,將該截面處的翼型進(jìn)行順時(shí)針旋轉(zhuǎn),旋轉(zhuǎn)角θ等于人流角?。新翼型各點(diǎn)坐標(biāo)相對(duì)于原坐標(biāo)的計(jì)算公式為


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys結(jié)果圖片11


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys結(jié)果圖片12

根據(jù)扭角修正方法,對(duì)某5 MW風(fēng)電機(jī)組葉片各截面處的翼型扭角進(jìn)行修正,修正后建立整個(gè)葉片的結(jié)構(gòu)模型,如圖2所示。


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片13


3邊界條件


耦合系統(tǒng)以MPCCI為數(shù)據(jù)交換平臺(tái),采用弱耦合方式,數(shù)據(jù)傳遞過程中固體與流體分別求解,只在禍合界面上進(jìn)行數(shù)據(jù)交換,禍合過程如圖3所示。流固藕合交界面上分別滿足邊界位移及應(yīng)力的連續(xù)性,即


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片14


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)的效果圖片15


4算例分析


4.1模型基本參數(shù)

以某年平均風(fēng)速8 m/s 的二類風(fēng)電場中5 MW風(fēng)電機(jī)組葉片為研究對(duì)象,該機(jī)組風(fēng)輪直徑127 m,葉片長62 m,額定風(fēng)速12.2 m/s,風(fēng)輪額定轉(zhuǎn)速13.5 r/min。在額定風(fēng)速及額定轉(zhuǎn)速條件下,根據(jù)上述方法對(duì)葉片進(jìn)行計(jì)算,得到各翼展處的參考風(fēng)速及修正角度如表1所示。


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)課程圖片16

  
根據(jù)表1的參考風(fēng)速及修正角度,進(jìn)行風(fēng)電機(jī)組葉片結(jié)構(gòu)建模,得到修正前后葉片的結(jié)構(gòu)模型(圖2)。對(duì)修正后的結(jié)構(gòu)模型進(jìn)行有限元離散化,得到流場計(jì)算的全局網(wǎng)格及局部網(wǎng)格單元如圖4所示。計(jì)算節(jié)點(diǎn)共計(jì)30萬,網(wǎng)格單元共計(jì)86萬。

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys培訓(xùn)課程圖片17



采用等效質(zhì)量、剛度分布的方法建立風(fēng)電機(jī)組葉片殼結(jié)構(gòu)模型。葉根段厚度為0.12 m,葉尖段為0.012 m,中間各段厚度沿展向線性變化,網(wǎng)格采用三角形殼單元S3R。節(jié)點(diǎn)總數(shù)5000,單元總數(shù)10000。用于結(jié)構(gòu)計(jì)算的葉片有限元網(wǎng)格單元如圖5所示。


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys workbanch圖片18


4.2結(jié)果分析

風(fēng)輪旋轉(zhuǎn)模擬模型與均勻來流初始模型的風(fēng)電機(jī)組葉片表面壓力計(jì)算結(jié)果對(duì)比如圖6所示。兩種模型的葉片表面壓力呈現(xiàn)相同的變化規(guī)律,由葉根到葉尖逐漸減小;靠近葉尖處同一翼型截面上壓力沿前緣到尾緣先迅速減小到最小值,后逐漸增大。數(shù)值上,旋轉(zhuǎn)模擬模型葉片的表面壓力小于均勻來流初始模型,整體誤差小于5%,驗(yàn)證了風(fēng)輪旋轉(zhuǎn)模擬模型的可靠性。


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys workbanch圖片19

風(fēng)輪旋轉(zhuǎn)模擬模型與均勻來流初始模型的風(fēng)電機(jī)組葉片40 m展向位置翼型周圍氣壓、流速的計(jì)算結(jié)果如圖7及圖8所示

兩種模型葉片40 m展向位置翼型周圍氣壓分布,與MEXICO實(shí)驗(yàn)60%半徑位置測試壓力數(shù)據(jù)呈現(xiàn)相同的變化規(guī)律,翼型升力側(cè)壓力小于壓力側(cè)壓力;翼型壓力側(cè)表面壓力從前緣到尾緣呈現(xiàn)先減小后增大的變化趨勢,最大壓力位于前緣點(diǎn)靠下位置,最小壓力位于升力側(cè)30%弦長處。

兩種模型葉片40 m展向位置翼型周圍空氣流速分布規(guī)律相同,翼型升力側(cè)流速大于壓力側(cè)流速;翼型壓力側(cè)附近流速從前緣到尾緣先增大后減小,最小流速位于前緣點(diǎn)靠下位置,最大流速位于升力側(cè)約25%弦長處。

風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys workbanch圖片20
風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys圖片21


兩種模型葉片40 m展向位置翼型周圍氣壓及流速計(jì)算結(jié)果的整體誤差在3%以內(nèi)。

風(fēng)輪旋轉(zhuǎn)模擬模型與均勻來流初始模型的風(fēng)電機(jī)組葉片結(jié)構(gòu)位移計(jì)算結(jié)果對(duì)比如圖9所示。兩種模型葉片的結(jié)構(gòu)位移從葉根到葉尖逐漸增大,在葉尖達(dá)到最大值。旋轉(zhuǎn)模擬模型葉片的最大結(jié)構(gòu)位移計(jì)算結(jié)果為5.98 m,均勻來流初始模型葉片的最大結(jié)構(gòu)位移計(jì)算結(jié)果為5.74 m,誤差小于3%。


風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法ansys圖片22


結(jié)論

以風(fēng)切變形式模擬風(fēng)輪旋轉(zhuǎn)與來流風(fēng)速綜合效應(yīng),對(duì)風(fēng)電機(jī)組葉片各展向位置翼型的扭角進(jìn)行修正,提出一種用于風(fēng)電機(jī)組葉片流固禍合數(shù)值模擬的風(fēng)輪旋轉(zhuǎn)模擬模型。在輪轂中心處風(fēng)速為額定風(fēng)速條件下,應(yīng)用風(fēng)輪旋轉(zhuǎn)模擬模型對(duì)風(fēng)電機(jī)組葉片進(jìn)行三維流固耦合數(shù)值計(jì)算,并與均勻來流初始模型計(jì)算結(jié)果及文獻(xiàn)實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比分析。風(fēng)輪旋轉(zhuǎn)模擬模型葉片40 m展向位置翼型周圍氣壓的仿真計(jì)算結(jié)果與MEXICO實(shí)驗(yàn)測試結(jié)果的分布規(guī)律相同。兩種模型葉片的表面壓力、翼型周圍氣壓及流速計(jì)算結(jié)果的分布規(guī)律相同,風(fēng)輪旋轉(zhuǎn)模擬模型得到的葉片面壓力小于均勻來流初始模型的計(jì)算結(jié)果,整體誤差小于5%;葉片40 m展向位置翼型周圍氣壓及流速計(jì)算結(jié)果,以及葉片最大結(jié)構(gòu)位移計(jì)算結(jié)果的誤差均小于3%。

作者:王偉龍,田德,鄧英,林俊杰



開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:風(fēng)電機(jī)組葉片流固耦合的數(shù)值模擬方法 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識(shí) ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn) 

編輯
在線報(bào)名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com