王永鼎，陳小鋒

（上海海洋大學 工程學院，上海 201306）

0 引言

海流機是海洋能開發(fā)利用的獲能機構(gòu)，其形式主要分為水平軸海流機和垂直軸海流機。海流機葉片是捕獲海流能的關(guān)鍵部件。水平軸海流機葉片是復雜的三維曲面體。海流機實際工作在近海岸的海底，由于海底面的存在，海流會形成呈梯度分布的邊界層速度，對海流機葉片表面負載產(chǎn)生一定的影響。海洋表面周期性波浪運動，也會影響海流機獲能效率。正常旋轉(zhuǎn)過程中，海流機葉片上可能會發(fā)生空化現(xiàn)象，導致葉片發(fā)生空蝕，縮短葉片的使用壽命[2]。因此，通過數(shù)值計算預測葉片的水動力性能具有重要意義。

海流機葉片的設計類似于風力機，但海流不同于氣流，海流作用在葉片上的載荷更大，工況更加復雜多變。葉片會發(fā)生繞展向的扭轉(zhuǎn)，軸向的揮舞和旋轉(zhuǎn)平面內(nèi)的擺振。浙江大學徐學濤等在文獻[3]中對水平軸潮流能發(fā)電裝置的葉片基于特定海況進行系統(tǒng)設計和校核。文獻[4]和文獻[5]分別對水平軸風力機葉片進行參數(shù)化建模，計算出葉片各截面處極限載荷，通過動力學分析得到葉片振動特性。文獻[6]基于葉素-動量理論，建立了水平軸海流機水動力學設計和校核模型。本文根據(jù)特定海況，參照風力機相關(guān)理論，設計海流機葉片實體外形并通過數(shù)值計算分析其模態(tài)和振型。

1 葉片水動力外形設計

1.1 海況和葉輪基本參數(shù)的確定

2000年前后，杭州灣跨海大橋橋址潮流調(diào)查時測得數(shù)據(jù)顯示，多處潮流流速達4～5m/s，最大流速可達5m/s以上[9]，此處取額定流速2.5m/s。假設海流機所處水深較深，不受海面波浪的影響，并認為海流流速在整個葉片掃略面內(nèi)均勻分布，不存在速度梯度。葉輪直徑D根據(jù)功率來計算，額定功率計算公式[8]：

額定功率P取100kW，海水密度ρ約為1027kg.m-3，額定流速V為2.5m/s，功率系數(shù)Cp預估0.4，葉輪半徑：

輪轂直徑為葉輪直徑0.1D。輪轂迎流端選取合理的流線型外形，結(jié)構(gòu)強度也要有保證。兩種介質(zhì)物理性質(zhì)對比如表1所示。

表1 海水和空氣對比

海流機葉輪參數(shù)整理如表2所示。

表2 葉輪參數(shù)

1.2 參數(shù)化建模

利用MATLAB聯(lián)合SolidWorks對葉片進行參數(shù)化建模。設計葉片參數(shù)如表3所示。

表3 葉片幾何參數(shù)

通過葉素-動量（BEM）理論[10]和Wilson設計法[1]計算葉素軸向和周向誘導因子a和b，然后就可以求得葉素的弦長和扭角等水動力外形參數(shù)。a和b需經(jīng)過反復迭代計算得到。算法流程如圖1示所。

圖1 Wilson設計方法

在MATLAB中采用最小二乘法將計算所得弦長和扭角進行曲線擬合，翼型截面沿展向半徑r作為自變量，其中弦長和扭角的擬合多項式為：

圖2和圖3分別為17個截面弦長和扭角優(yōu)化值和擬合曲線：

圖2 扭角擬合值和優(yōu)化計算值

圖3 弦長擬合值和優(yōu)化計算值

通過Profili軟件獲取NACA4415翼型，將原始翼型坐標（x0，y0）先移動到旋轉(zhuǎn)中心（0.25C，0），其中C為弦長。然后將各個截面的翼型按照圖3規(guī)律變換弦長，再各自繞旋轉(zhuǎn)中心旋轉(zhuǎn)對應扭角。

坐標變換公式如式(5)所示，變換前坐標（x0，y0）順時針繞點(a，b)旋轉(zhuǎn)β（單位：rad）后得到變換后坐標（x，y）。

海流機葉片輪廓如圖4所示。

圖4 水平軸海流機葉片

2 網(wǎng)格無關(guān)性驗證

與實驗驗證相比，數(shù)值計算周期短，成本低，但數(shù)值計算結(jié)果的可靠性受很多因素的影響，尤其是網(wǎng)格劃分質(zhì)量對計算結(jié)果的精度影響很大。劃分網(wǎng)格的數(shù)量，關(guān)鍵部位的網(wǎng)格密度，網(wǎng)格類型的選用會使計算結(jié)果的差異很大[11]。水平軸海流機葉片是復雜三維曲面，因此選用適應度較好的四面體網(wǎng)格。由于過于細密的網(wǎng)格對硬件要求較高，且對于計算結(jié)果的影響較小，在此選用中等密度網(wǎng)格。網(wǎng)格數(shù)量的驗證如表4所示。

表4 網(wǎng)格無關(guān)性驗證

由表4可知，網(wǎng)格數(shù)達到5萬左右時，結(jié)構(gòu)靜力學分析時的最大變形量前后誤差很小。因此網(wǎng)格數(shù)取5萬左右即可滿足計算精度。

3 靜力分析

靜力分析是通過數(shù)值計算獲取葉片在不同工況下的變形情況，從而預測葉片疲勞破壞發(fā)生的位置。葉片選用材料的物理特性如表5所示。

表5 葉片材料參數(shù)

靜力分析中分別設定如表6所示的三種工況，工況1為額定工況；工況2表示海流流速較快，葉輪轉(zhuǎn)速較高；工況3表示海流流速過大導致海流機停車。

表6 三種常見工況

計算結(jié)果如表7所示。

表7 靜力分析結(jié)果

由于篇幅有限，僅給出額定工況下的總最大變形圖5，最大等效應變圖6，最大等效應力圖7。

圖5 總最大變形

圖6 等效彈性應變

圖7 最大等效應力

通過對比三種工況發(fā)現(xiàn)，流速越大，葉片表面受到的壓力載荷也越大，葉片沿軸向變形逐漸增大，最大位移都發(fā)生在葉尖處。最大應力和最大應變發(fā)生在距葉尖1/3處，說明此處容易疲勞破壞，長時間工作后容易造成葉片折斷。

4 模態(tài)分析

葉片在海流中工作，同時受到水流沖擊，地磁場對葉片產(chǎn)生的重力和旋轉(zhuǎn)時的離心力作用。海流機葉片展向長，弦向短，是易發(fā)生振動的細長彈性體[3]。相比于氣流，海流流速穩(wěn)定，密度大，流向變化周期長，海流作用在葉片上的載荷比氣流要大，而且這些載荷同時具有交變和隨機性[7]，在設計工況下，海流機葉片發(fā)生振動不可避免。為避免葉片因共振而造成材料的疲勞破壞，需對其模態(tài)分析，了解葉片的振動特性。

海流機葉片的葉根是從橢圓形根部截面通過過渡表面放樣到翼型輪廓部分的放樣體。在進行單葉片模態(tài)分析時，可以將葉根截面的六個自由度全部約束，另一端自由。葉片前十階模態(tài)計算結(jié)果如表8所示。

表8 固有頻率和振型

根據(jù)振動理論，物體在發(fā)生振動的過程中，能量一般主要集中在第一和第二階模態(tài)[8]。因此這里僅給出海流機葉片的一階和二階模態(tài)的總變形云圖8和圖9。

圖8 一階揮舞

圖9 二階擺振和揮舞

葉輪在額定轉(zhuǎn)速下工作，作用在葉片上的外界載荷包括海流作用力、葉片自身重力和旋轉(zhuǎn)離心力。葉片以Y軸作為轉(zhuǎn)軸，Z軸作為展向，X軸在葉輪掃略的平面內(nèi)。葉片在旋轉(zhuǎn)過程中的固有頻率是海流機整體機械結(jié)構(gòu)耦合后的固有頻率，與靜止時測得的單葉片的固有頻率稍有偏差，但一般忽略這種很小的頻率變化，故可認為單個葉片的固有頻率就是其運行時的固有頻率[8]。

海流機在海水里主要受海流周期性變化和轉(zhuǎn)輪轉(zhuǎn)動的激勵[2]。海流的運動具有周期性，半日潮的潮流流速一天變化兩次，每次周期約為8小時，由此推導出海流激勵為：

海流流速一般為0.5～5m/s之間，根據(jù)上式計算葉輪轉(zhuǎn)速區(qū)間，進而可得葉輪轉(zhuǎn)動的激勵頻率：

根據(jù)振動原理，葉片發(fā)生共振條件是外界激勵與葉片的固有頻率相同或接近。由上計算可知，海流激勵和葉輪的激勵范圍分別為3.47×10-5Hz和0.139Hz～1.39Hz，由表8可知，上述兩種頻率不但不在海流機葉片前十階振動的頻率范圍內(nèi)，而且兩者的范圍相差甚遠，故不會發(fā)生共振。同時由表8中海流機葉片的第一、二階模態(tài)振型可知，揮舞和擺振是其主要的振動形式。葉片的振動節(jié)點處在葉片距葉根三分之二處，可見此位置易發(fā)生疲勞破壞。

5 結(jié)語

1）運用Wilson設計原理，以最大能量利用率為目標函數(shù)，通過MATLAB和SolidWorks聯(lián)合仿真，對水平軸海流機葉片進行參數(shù)化建模，得到海流機葉片實體模型。

2）利用ANSYS對葉片進行靜力分析發(fā)現(xiàn)，額定工況下選用結(jié)構(gòu)鋼材料，葉片總的最大變形量為42.747mm，最大的等效應變?yōu)?.067mm/mm，最大的等效應力為13268Mpa。

3）利用ANSYS對葉片進行模態(tài)分析發(fā)現(xiàn)，設計工況下海流機葉片的一階固有頻率為1318Hz，其主要的振動型式為軸向揮舞；二階固有頻率4164.3Hz，主要的振動型式為二階軸向揮舞和周向擺振。海流激勵和葉輪的激勵范圍分別為3.47×10-5Hz和0.139Hz～1.39Hz，兩者相差甚遠，根據(jù)振動理論，不會發(fā)生共振。葉片振動節(jié)點主要位于距葉跟三分之二處，工程應用中應加強這些部位的結(jié)構(gòu)強度。