王 波 張以良 熊 鷹

(海軍工程大學(xué)艦船工程系1) 武漢 430033) (92196部隊(duì)2) 青島 266000)

船用螺旋槳水動(dòng)力性能優(yōu)化設(shè)計(jì)*

王 波1)張以良2)熊 鷹1)

(海軍工程大學(xué)艦船工程系1)武漢 430033) (92196部隊(duì)2)青島 266000)

為尋找效率和空泡性能優(yōu)良的螺旋槳，利用螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)對(duì)某吊艙螺旋槳開(kāi)展迭代優(yōu)化，并對(duì)優(yōu)化結(jié)果進(jìn)行試驗(yàn)驗(yàn)證.螺旋槳各參數(shù)分布曲線(xiàn)采用貝塞爾樣條曲線(xiàn)擬合，以確保各參數(shù)在優(yōu)化過(guò)程中的光順性；應(yīng)用遺傳算法，結(jié)合螺旋槳面元法的定常水動(dòng)力性能預(yù)報(bào)工具，以提高螺旋槳效率及空泡性能為優(yōu)化目標(biāo)，對(duì)螺旋槳各參數(shù)分布進(jìn)行優(yōu)化及相關(guān)分析探討.結(jié)果表明，優(yōu)化槳較原槳效率略低，但空泡初生更晚，相同條件下空化面積更小，空泡性能得到有效改善.

螺旋槳；優(yōu)化；遺傳算法；定常面元法；貝塞爾樣條曲線(xiàn)

0 引 言

在螺旋槳設(shè)計(jì)中，眾多的幾何參數(shù)與螺旋槳各性能之間關(guān)系非常復(fù)雜，傳統(tǒng)設(shè)計(jì)方法中，很多參數(shù)往往是根據(jù)設(shè)計(jì)者的經(jīng)驗(yàn)進(jìn)行多次調(diào)整后形成最終設(shè)計(jì)方案.該“修改-驗(yàn)證-修改”的人工嘗試方法得到螺旋槳方案可能達(dá)到效率及空泡性能都較佳的效果，但并不一定是最佳效果；而且該方法需要大量的實(shí)際經(jīng)驗(yàn).近些年來(lái)，通過(guò)智能優(yōu)化算法對(duì)螺旋槳綜合性能進(jìn)行優(yōu)化研究越來(lái)越多[1-6].

本文利用螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)[7],以螺旋槳效率及空泡性能為優(yōu)化目標(biāo)，對(duì)某吊艙螺旋槳進(jìn)行優(yōu)化設(shè)計(jì)研究，并對(duì)優(yōu)化結(jié)果進(jìn)行模型試驗(yàn)驗(yàn)證.

1 螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)

1.1 簡(jiǎn)介

本螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)是通過(guò)調(diào)整原槳的幾何形狀，建立起優(yōu)化計(jì)算的樣本空間，再利用遺傳算法[8]配合面元法定常水動(dòng)力預(yù)報(bào)工具從中尋找最優(yōu)方案的方式實(shí)現(xiàn)螺旋槳的優(yōu)化.為保證優(yōu)化過(guò)程中螺旋槳各參數(shù)的光順性，螺旋參數(shù)化建模采用貝塞爾樣條曲線(xiàn)擬合[9].其優(yōu)化流程見(jiàn)圖1.

圖1 螺旋槳綜合性能優(yōu)化設(shè)計(jì)系統(tǒng)流程圖

1.2 優(yōu)化原理及方法

1) 優(yōu)化目標(biāo)的建立 以螺旋槳效率及空泡性能為優(yōu)化目標(biāo)，空泡性能以0.7R半徑處葉背最小壓力系數(shù)為判斷依據(jù).

(1)

式中：η0為螺旋槳效率；Cp為0.7R半徑處葉背最小壓力系數(shù).定義壓力系數(shù)為

(2)

式中：p0為葉背面某處的壓力，Pa；pv為水的氣化壓力，Pa；n為螺旋槳轉(zhuǎn)速，r/s；D為螺旋槳直徑，m.

利用權(quán)重系數(shù)法建立螺旋槳效率和空泡性能的多目標(biāo)函數(shù)，總的目標(biāo)函數(shù)可以表示為

f=w1f1+w2f2

(3)

式中：w1為螺旋槳效率權(quán)重系數(shù)；w2為空泡性能權(quán)重系數(shù).

2) 約束條件的建立 在優(yōu)化過(guò)程中，必須對(duì)優(yōu)化設(shè)計(jì)方案進(jìn)行評(píng)估分析及篩選，以保證最終得到的方案合理，符合實(shí)際要求.對(duì)不符合要求的方案的目標(biāo)函數(shù)實(shí)行懲罰.從以下3個(gè)方面進(jìn)行約束.

(1) 幾何約束.基于貝塞爾樣條曲線(xiàn)擬合實(shí)現(xiàn)螺旋槳參數(shù)化建模后，建立起了擬合曲線(xiàn)控制點(diǎn)與參數(shù)分布曲線(xiàn)的對(duì)應(yīng)關(guān)系，在優(yōu)化過(guò)程中，調(diào)整控制點(diǎn)的位置既可以改變參數(shù)分布曲線(xiàn)的形狀.在這一過(guò)程中，為避免曲線(xiàn)出現(xiàn)纏繞等不合理現(xiàn)象，對(duì)新生成的控制點(diǎn)位置進(jìn)行約束;(2) 水動(dòng)力約束.為保證優(yōu)化后的螺旋槳產(chǎn)生的推力符合要求，需對(duì)螺旋槳的推力系數(shù)進(jìn)行約束;(3) 強(qiáng)度約束.優(yōu)選過(guò)程中，為保證優(yōu)化后的螺旋槳滿(mǎn)足強(qiáng)度要求，對(duì)槳葉最大厚度進(jìn)行約束.

(4)

式中：VXi為擬合曲線(xiàn)控制點(diǎn)的X坐標(biāo)值；kt為螺旋槳推力系數(shù);σ1為推力系數(shù)誤差限；t為螺旋槳槳葉最大厚度，約束優(yōu)化方案的槳葉最大厚度不小于初始槳.

2 算例與分析

2.1 優(yōu)化對(duì)象

本文以某吊艙槳模型為母型，利用螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)進(jìn)行優(yōu)化.該母型槳主要參數(shù)見(jiàn)表1.

表1 槳的主要參數(shù)

2.2 優(yōu)化設(shè)置

優(yōu)化時(shí)保持螺旋槳直徑不變，將所有徑向參數(shù)(弦長(zhǎng)、螺距、縱傾、側(cè)斜、最大厚度及最大拱度)和剖面參數(shù)(弦向厚度分布和弦向拱度分布)作為優(yōu)化參數(shù).其中螺旋槳徑向參數(shù)分布采用三階貝塞爾曲線(xiàn)擬合[10-12]，剖面參數(shù)采用五階貝塞爾曲線(xiàn)擬合.限制條件設(shè)置：推力系數(shù)誤差限設(shè)置為±0.1%；強(qiáng)度約束條件設(shè)置槳葉最大厚度分布不小于原槳槳葉最大厚度分布.優(yōu)化目標(biāo)設(shè)置：效率和空泡性能.其中空泡性能以0.7R半徑處葉背最小壓力系數(shù)為判斷依據(jù)，采用權(quán)重系數(shù)法，權(quán)重系數(shù)分別取為0.5.遺傳算法設(shè)置：種群大小設(shè)置為30；優(yōu)化代數(shù)設(shè)置為500；交叉概率設(shè)置為0.8.

2.3 優(yōu)化結(jié)果

表2 槳優(yōu)化結(jié)果

表2表明，優(yōu)化槳較原槳推力系數(shù)保持不變，效率下降0.129%；而0.7R處葉背最小壓力系數(shù)較原槳提高了5.9%.可見(jiàn)，優(yōu)化槳效率略低，但具有更優(yōu)良的空泡性能.

圖2為是優(yōu)化前后螺旋槳徑向參數(shù)的變化.圖3為優(yōu)化前后剖面參數(shù)的變化情況.

圖2 優(yōu)化前后螺旋槳徑向參數(shù)分布的比較

圖3 優(yōu)化前后剖面參數(shù)分布的比較

圖2表明，弦長(zhǎng)分布、螺距分布及最大拱度分布變化較明顯，對(duì)螺旋槳性能變化較為敏感.弦長(zhǎng)在內(nèi)半徑有所減小，在外半徑有所增大；螺距在槳轂附近有所降低，在徑向其他位置均有所增加；縱傾在槳轂至0.7R半徑有所增加，但從0.7R半徑至葉梢有略微的減?。粋?cè)斜在槳轂至0.5半徑處有所增加，其他半徑處保持不變，側(cè)斜平衡度有所降低；最大厚度在槳轂附近有所增加；最大拱度分布在槳轂附近有所增加外，在徑向其他位置均是減小.

圖3表明，剖面最大厚度略微的往導(dǎo)邊方向靠近，而最大拱度略微的往隨邊方向靠近.

2.4 試驗(yàn)驗(yàn)證分析

本實(shí)驗(yàn)在海軍工程大學(xué)空泡水洞實(shí)驗(yàn)室完成.空泡水洞工作段封閉，長(zhǎng)2.6 m，截面為0.6 m×0.6 m帶圓角的矩形.上下水平段中心線(xiàn)的高度約為10 m，2個(gè)垂直段中心線(xiàn)間距為18 m，收縮比6.25∶1，調(diào)壓范圍0.05～2 bar，空泡數(shù)0.2～3.

試驗(yàn)?zāi)康氖潜容^優(yōu)化前后螺旋槳模型的空泡性能，以此驗(yàn)證優(yōu)化的合理性.原槳模型采用合金銅材質(zhì)，優(yōu)化模型采用鋁合金.模型見(jiàn)圖4.

圖4 試驗(yàn)槳模型

本次試驗(yàn)從2個(gè)方面來(lái)比較螺旋槳空泡性能的優(yōu)良性.(1) 比較優(yōu)化前后槳模的空泡初生的早晚；(2) 比較優(yōu)化前后螺旋槳在相同水速空化數(shù)和相同進(jìn)速系數(shù)下空化面積大小.

(5)

(6)

式中:V為水洞水流速度，m/s；D為槳模直徑，m；n為槳模轉(zhuǎn)數(shù)，r/s；ρ為水密度,kg/m3；p0為水洞中心壓力，Pa；pv為水汽化壓力，Pa.

1) 空泡初生速度對(duì)比 在常壓下，保持水洞水流速V=4.2 m/s，增加轉(zhuǎn)速，直至產(chǎn)生初生空泡，比較空泡初生的早晚.兩槳模產(chǎn)生初生空泡的工況見(jiàn)表3.表3表明，優(yōu)化槳空泡初生較原槳要更晚些.

表3 空泡初生試驗(yàn)實(shí)測(cè)工況

2) 相同水速空化數(shù)和進(jìn)速系數(shù)下空化面積大小對(duì)比 保持水流速度4.2 m/s，同時(shí)保持水洞在負(fù)0.76 bar，即保持相同的水速空化數(shù)σV=2.50.相同水速空化數(shù)下的試驗(yàn)觀(guān)測(cè)工況見(jiàn)表4.

表4 空泡對(duì)比試驗(yàn)實(shí)測(cè)工況

根據(jù)觀(guān)測(cè)到的結(jié)果，描繪出螺旋槳葉背空泡示意圖.通過(guò)觀(guān)察并對(duì)比不同工況下2螺旋槳的空化狀態(tài)，可以反映出優(yōu)化前后螺旋槳的空化效果.如圖5～圖8所示.

圖5 工況1下原槳與優(yōu)化槳試驗(yàn)對(duì)比

圖6 工況2下原槳與優(yōu)化槳試驗(yàn)對(duì)比

圖7 工況3下原槳與優(yōu)化槳試驗(yàn)對(duì)比

圖8 工況4下原槳與優(yōu)化槳試驗(yàn)對(duì)比

圖5～圖8表明，優(yōu)化槳在槳葉葉背上空泡面積比原槳的小.尤其是在高進(jìn)速系數(shù)的時(shí)候.

試驗(yàn)表明，優(yōu)化后的螺旋槳空泡初生較原槳更晚，同時(shí)在相同條件下優(yōu)化槳空化面積較原槳更小，有效地改善了優(yōu)化槳的空泡性能.

3 結(jié) 束 語(yǔ)

以提高螺旋槳效率及0.7R半徑處葉背最小壓力系數(shù)為優(yōu)化目標(biāo)，利用本單位自行編制的螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)對(duì)某吊艙模型槳各參數(shù)進(jìn)行迭代優(yōu)化.數(shù)值優(yōu)化結(jié)果表明，效率降低0.129%，0.7R半徑處葉背最小壓力系數(shù)提高5.9%.空泡試驗(yàn)對(duì)比驗(yàn)證，優(yōu)化槳較原槳空泡初生更晚；在相同條件下優(yōu)化槳的空化面積較原槳小，從而表明螺旋槳空泡性能得到有效改善.

[1]王國(guó)強(qiáng)，董世湯.船舶螺旋槳理論與應(yīng)用[M].哈爾濱：哈爾濱工程大學(xué)出版社.2005.

[2]CHEN Jenghorng， SHIH Yushan．Basic design of a series propeller with vibration consideration by genetic algorithm[J]．Marine Science and Technology，2007(12)：119-129.

[3]BERTETTA D，BRIZZOLARA S． CPP propeller cavitation and noise optimization at different pitches with panel code and validation by cavitation tunnel measurements[J]．Ocean Engineering，2012(53)：177-195．

[4]TAHERI R， MAZAHERI K．Hydrodynamic optimization of marine propeller using gradient and non-gradient based algoritnms[J]．Acta Polytechnica Hungarica，2013，10(3)：221-237．

[5]黃樹(shù)權(quán)．基于遺傳算法的螺旋槳性能優(yōu)化研究[D]．大連：大連理工大學(xué)，2009．

[6]趙 威，楊晨?。奥菪龢菥嗉肮岸鹊膬?yōu)化設(shè)計(jì)研究[J]．中國(guó)造船，2010，51(1)：1-8．

[7]張以良，熊 鷹.螺旋槳多參數(shù)協(xié)同優(yōu)化系統(tǒng)V1.0：中國(guó)，2014SR038756[P].2014.

[8]雷英杰，張善文．MATLAB遺傳算法工具箱及應(yīng)用[M]．西安：西安電子科技大學(xué)出版社，2005.

[9]朱心熊．自由曲線(xiàn)曲面造型技術(shù)[M]．北京：科學(xué)出版社，2000.

[10]張以良，熊 鷹，時(shí)立攀．一種螺旋槳參數(shù)化建模的方法[J]．艦船科學(xué)技術(shù)，2015，37(1)：34-38.

[11]姚震球，高 慧，楊春蕾，螺旋槳三維建模與水動(dòng)力數(shù)值分析[J]．船舶工程，2008(6)：55-58.

[12]王 超，黃 勝，解學(xué)參.基于CFD方法的螺旋槳水動(dòng)力性能預(yù)報(bào)[J]．海軍工程大學(xué)學(xué)報(bào)，2008(4)：48-52.

Design and Optimization of Hydrodynamic Performance of Marine Propeller

WANG Bo1)ZHANG Yi liang2)XIONG Ying1)

(DepartmentofNavalEngineering,NavalUniv.ofEngineering,Wuhan430033,China)1)(TheUnit92196thofPLA,Qingdao266000,China)2)

In order to seek for a propeller that both good at efficiency and cavitation performance, an iterative optimization design of a pod propeller was accomplished by using propeller comprehensive performance optimization design system-with independent intellectual property rights, and model tests were conducted to verify the optimization results. The propeller parameters distributions are described by Bezier-spline, so that the curve fairing can be ensured. Aiming at improving the efficiency and cavitation performance, genetic algorithm combined with a prediction method of steady hydrodynamic characteristics of propeller based on surface panel method is utilized to optimize and analyze propeller parameters distribution. The optimization results show that the efficiency of the optimized propeller is slightly lower, but the cavitation inception is later and the cavitation area is smaller at the same condition comparing to the original propeller, the cavitation performance have been improved.

propeller; optimization; genetic algorithm; steady surface panel method; Bezier-spline

2015-04-10

*海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室研究基金項(xiàng)目資助(批準(zhǔn)號(hào)：1106)

U644.33

10.3963/j.issn.2095-3844.2015.04.022

王 波(1990- ):男，碩士生，主要研究領(lǐng)域?yàn)榕灤黧w動(dòng)力性能