張 旭， 劉安宇

(天津工業(yè)大學(xué) 機(jī)械工程學(xué)院， 天津 300387)

風(fēng)力發(fā)電已經(jīng)成為世界上最具發(fā)展前景的能源開發(fā)方式.目前，利用風(fēng)能最為行之有效的裝置是風(fēng)力發(fā)電機(jī)，其中最具重要性的部件即為葉片，因為其性能與風(fēng)力機(jī)運(yùn)行現(xiàn)狀及穩(wěn)定性存在直接關(guān)聯(lián).在風(fēng)力機(jī)的運(yùn)行過程中，葉片周圍空氣的流動對葉片產(chǎn)生的作用力促使葉片出現(xiàn)一定程度的形變，這一變化會導(dǎo)致周邊空氣的流向發(fā)生變化.為了得到精準(zhǔn)的葉片動力學(xué)響應(yīng)，首要要考慮葉片及流場間存在的作用，即必須考慮葉片與其周圍流場的耦合作用.

國外研究人員對流體和固體風(fēng)力發(fā)電機(jī)組進(jìn)行了大量研究.Tongchitpakdee等對偏航條件下的大型風(fēng)機(jī)進(jìn)行了數(shù)值模擬計算，并對其氣動性能進(jìn)行了系統(tǒng)分析[1]；Rachid等在定常情況及非定常情況下對NACA4415翼型的葉片彎曲及扭轉(zhuǎn)效應(yīng)進(jìn)行分析研究，延長了風(fēng)力機(jī)的使用壽命[2]；Kim等考慮了氣動彈性對葉片的影響，并對大型風(fēng)力機(jī)葉片性能進(jìn)行了預(yù)測[3].

自80年代中期至今，國內(nèi)研究人員對實體研究流固耦合問題與風(fēng)力發(fā)電機(jī)整機(jī)的流固耦合進(jìn)行了大量的研究.徐建源等利用ADAMS軟件建立了風(fēng)波聯(lián)合作用海上風(fēng)力機(jī)模型，并對其進(jìn)行流固耦合分析，得到了復(fù)雜工況下的風(fēng)力機(jī)葉片的動態(tài)響應(yīng)參數(shù)[4]；李少華等利用Fluent軟件并結(jié)合SSTk-ω模型對1.2 MW風(fēng)力機(jī)旋轉(zhuǎn)風(fēng)輪流場進(jìn)行模擬，但是并沒有考慮風(fēng)力機(jī)塔筒與流場間的相互作用對數(shù)據(jù)造成的影響[5]；胡丹梅等在水平軸風(fēng)力機(jī)模型不同尖速比條件下對風(fēng)輪下游流場速度進(jìn)行測量，獲得風(fēng)力機(jī)葉片尾跡的流場定量信息，由于該測量方法容易受外界因素的影響，因此，流動情況與實際情況有一定誤差[6-7]；田琳琳等結(jié)合制動盤理論與CFD方法，并通過Fluent軟件對9臺風(fēng)力機(jī)的尾流互相干擾情況進(jìn)行模擬，得到風(fēng)力機(jī)尾流相互影響的結(jié)果[8]；陳海萍等研究了葉片周圍流場對葉片結(jié)構(gòu)特性的影響[9]；潘萍萍等通過Fluent軟件對1.5 MW風(fēng)力發(fā)電機(jī)塔筒進(jìn)行流固耦合分析，得到了風(fēng)力機(jī)塔筒的動態(tài)評估方法[10].

本文應(yīng)用ANSYS軟件對改進(jìn)后的2 MW風(fēng)力機(jī)葉片模型進(jìn)行流固耦合分析，通過對傳統(tǒng)翼型NACA4412的改進(jìn)，應(yīng)用鋪層設(shè)計把葉梢、葉根、展向等厚度改成漸變式.首先得到風(fēng)力機(jī)葉片及風(fēng)輪的流場計算結(jié)果，再將風(fēng)力機(jī)流場的計算結(jié)果施加到固體計算中，得到風(fēng)力機(jī)的應(yīng)力分布，進(jìn)而得出應(yīng)力最大值及變形量最大位置，所得結(jié)果可以為風(fēng)力機(jī)的設(shè)計和安全運(yùn)行提供指導(dǎo)[11].

1 流固耦合方法

流固耦合問題從數(shù)據(jù)傳輸?shù)慕嵌葋砜矗溆嬎憧梢苑譃閮煞N：單向流固耦合和雙向流固耦合.單向流固耦合是指通過流體計算獲得的壓力、速度等數(shù)據(jù)可以通過連接界面導(dǎo)入固體結(jié)構(gòu)，或者將從固體計算獲得的網(wǎng)格位移轉(zhuǎn)移到流體的計算.

計算流體力學(xué)簡稱CFD，主要使用離散點上的一系列變量值代替時域和空間域中的一系列物理量.通過建立可以表示變量之間關(guān)系的代數(shù)方程，將基本方程離散化，通過CFD的數(shù)值模擬，可以顯示流動領(lǐng)域的運(yùn)動.在通過CFD分析的工程應(yīng)用中，可以深入研究耦合問題以實現(xiàn)實驗指南，達(dá)到節(jié)約成本的目的.流固耦合分析流程圖如圖1所示.

圖1 耦合系統(tǒng)流程示意圖Fig.1 Schematic diagram of coupling system

2 三維模型建立

本文利用NASYS APDL建立風(fēng)力機(jī)葉片的三維實體模型.采用傳統(tǒng)翼型NACA4412應(yīng)用復(fù)合材料建立葉片有限元模型如圖2a所示，該模型總長44.5 m，翼型最大弦長3.1 m，根部為直徑2.34 m的圓形，質(zhì)量為10.6 t.為了加強(qiáng)葉片穩(wěn)定性，在葉片中部添加兩個厚為1.6 cm的加強(qiáng)筋，葉根部也進(jìn)行了加厚，具體模型如圖2b所示.葉片的建模工作主要在傳統(tǒng)翼型NACA4412模型下進(jìn)行改進(jìn)，應(yīng)用鋪層設(shè)計把葉梢、葉根、展向等厚度改成漸變式，這樣不僅更加接近實體模型，在強(qiáng)度、剛度上也有所加強(qiáng)，風(fēng)輪模型如圖2c所示.

3 靜止單葉片流固耦合分析

3.1 流場計算分析

根據(jù)葉片結(jié)構(gòu)在ANSYS軟件中繪得實體模型，有針對地對流體實際范圍進(jìn)行了設(shè)定，同時進(jìn)行了網(wǎng)格劃分，不過因常規(guī)風(fēng)力機(jī)體積過大，因此在制定計算域的過程中需要將其設(shè)置為葉片的數(shù)倍以上.設(shè)置的流體域范圍為202 m×103 m×63 m，所劃分的網(wǎng)格數(shù)為4 601 159，計算域及計算域網(wǎng)格模型如圖3所示.

在進(jìn)行流場分析時應(yīng)先抑制葉片固體域，以實現(xiàn)對流體域的保護(hù)，設(shè)置流體控制溫度為25 ℃.依照入口實際條件對其速度進(jìn)行控制，并將這一設(shè)定速度稱為恒速，流體的方向則需要垂直于入口界面，結(jié)合實際情況設(shè)置流體速度大小為13 m/s，而流體出口處設(shè)定為自由流出，風(fēng)力機(jī)葉片邊界條件設(shè)置為光滑無滑移壁面.風(fēng)場流速圖如圖4a所示，葉片在迎風(fēng)面和背風(fēng)面的壓力分布如圖4b、c所示，在背風(fēng)面葉片受到較小的壓力，而迎風(fēng)面葉片受到壓力最大.

圖2 葉片及風(fēng)輪模型Fig.2 Models for blade and wind wheel

圖3 計算域及計算域網(wǎng)格模型(靜止?fàn)顟B(tài))Fig.3 Models for computational domain and computationaldomain grid (static state)

圖4 風(fēng)場流速及迎風(fēng)面、背風(fēng)面壓力分布Fig.4 Wind field velocity and pressure distributionof windward and leeward surfaces

3.2 葉片靜力分析

葉片翼面材料為玻璃鋼，該材料為線性材料，用MP命令定義其材料特性.其各層的線性材料特性為正交異性，具體參數(shù)為：展向模量42.6×109Pa，徑向模量16.5×109Pa，剪切模量5.5×109Pa，泊松比0.22，密度1 950 kg/m3.

通過Static Structural接收風(fēng)場的分析結(jié)果數(shù)據(jù)，并實現(xiàn)結(jié)構(gòu)靜力的相關(guān)研究.首先單獨保留葉片，確保流程部分處于制約狀態(tài).由于葉片實際網(wǎng)格劃分直接影響最后的研究結(jié)果，因此在劃分過程中要確保其規(guī)則性.在Static Structural模塊下選擇Engineering Data，定義材料參數(shù)密度、彈性模量、泊松比和剪切模量，導(dǎo)入固體模型，定義網(wǎng)格劃分參數(shù)劃分網(wǎng)格，設(shè)置的網(wǎng)格數(shù)為980 916，定義固定約束.劃分網(wǎng)格后的模型如圖5所示.

圖5 葉片網(wǎng)格模型Fig.5 Model for blade grid

4 旋轉(zhuǎn)狀態(tài)風(fēng)輪流固耦合分析

4.1 流場計算分析

將Fluent模塊的Solution與Static Structural模塊的Setup建立連接進(jìn)行單向流固耦合分析，導(dǎo)入風(fēng)輪分別建立旋轉(zhuǎn)流體域和外部空氣流體域，旋轉(zhuǎn)流體域半徑為48 m，高度為6 m.外部空氣流體域尺寸為606 m×296 m×296 m，設(shè)置網(wǎng)格劃分參數(shù)及膨脹層，計算域及計算域網(wǎng)格模型如圖6所示.

圖6 計算域及計算域網(wǎng)格模型(旋轉(zhuǎn)狀態(tài))Fig.6 Models for computational domain and computationaldomain grid (rotation state)

設(shè)置室溫為25 ℃，入口速度為13 m/s，定義旋轉(zhuǎn)軸為y軸，旋轉(zhuǎn)速度為1.894 49 rad/s，風(fēng)場流速圖如圖7a所示，通過Contour云圖顯示的壓力云圖如圖7b所示，可以看出壓力從葉根到葉尖逐漸減少，葉根最大，而在葉尖處最小.

圖7 風(fēng)場流速及風(fēng)輪壓力分布Fig.7 Wind field velocity and pressuredistribution of wind wheel

4.2 風(fēng)輪結(jié)構(gòu)靜力分析

在Static Structural模塊下選擇Engineering Data進(jìn)行材料參數(shù)定義：密度為1 950 kg/m3，彈性模量為42.6×109Pa，泊松比為0.22，剪切模量為5.5×109Pa，導(dǎo)入固體模型，劃分網(wǎng)格，設(shè)置葉輪網(wǎng)格數(shù)為1 393 827，定義固定約束.劃分網(wǎng)格后的模型如圖8所示.

圖8 風(fēng)輪網(wǎng)格模型Fig.8 Model for wind wheel grid

4.3 結(jié)果輸出

導(dǎo)入流體分析結(jié)果即壓力，選擇導(dǎo)入壓力的面為葉片外表面，得到應(yīng)變分布如圖9所示.由圖9可知，在額定速度下風(fēng)力機(jī)運(yùn)行時，風(fēng)力機(jī)結(jié)構(gòu)受風(fēng)場的影響出現(xiàn)變形，變形的趨勢與在風(fēng)場中所受到的葉片表面壓力分布是基本吻合的，此外在形變量最大的葉尖處，相對于葉根處的約束部分最遠(yuǎn)，并且此處葉片材料是最薄的，與文獻(xiàn)[11]實驗結(jié)果基本吻合.

圖9 風(fēng)輪應(yīng)變分布Fig.9 Stain distribution of wind wheel

5 結(jié) 論

本文合理應(yīng)用ANSYS Workbench實現(xiàn)風(fēng)力系統(tǒng)中風(fēng)輪及葉片的單向流固耦合分析，得到兩者在運(yùn)行數(shù)據(jù)上存在的關(guān)聯(lián).實驗中得到的數(shù)據(jù)為進(jìn)行葉片疲勞及共振頻率測試提供了一定的參考，為之后的研究提供基礎(chǔ).在實際運(yùn)行工作中，葉片的固有頻率為耦合之后整個系統(tǒng)的固有頻率，這個固有頻率包括了風(fēng)力機(jī)塔架在系統(tǒng)中對葉片的影響，系統(tǒng)中旋轉(zhuǎn)離心力的作用，風(fēng)場在系統(tǒng)中對葉片的影響，系統(tǒng)中重力作用等一系列的作用，這些影響與單個葉片在靜止?fàn)顟B(tài)下的模態(tài)頻率有一定的區(qū)別，但一般而言頻率的變化很小.在今后的葉片優(yōu)化設(shè)計過程中，同樣需要關(guān)注該部分問題，以便得到的模型與實際運(yùn)作更為接近，減小并優(yōu)化設(shè)計周期.