葉禮裕 王 超 黃 勝 常 欣

（哈爾濱工程大學船舶工程學院 哈爾濱 150001）

0 引 言

螺旋槳幾何造型復雜，包含螺距、縱傾、側斜、拱度、弦長、槳葉剖面等幾何參數(shù)，這些幾何參數(shù)相互聯(lián)系和制約，影響著整個螺旋槳的水動力、噪聲、空泡、振動等性能.在螺旋槳的理論設計中，螺旋槳的幾何參數(shù)往往是給定的，但是給定的螺旋槳幾何參數(shù)得到的螺旋槳性能并不能滿足實際要求，需要進行參數(shù)的優(yōu)化設計得到性能更佳的槳.由于螺旋槳的各個性能的相互沖突，要想得到螺旋槳參數(shù)使得各方面性能達到最優(yōu)是很難的，改善螺旋槳的一些性能就有可能使其他性能的惡化.因此，如何權衡螺旋槳幾何參數(shù)得到接近理想的槳一直是研究人員們所關注的.

粒子群優(yōu)化（particle swarm optimization，PSO）是 Kennedy和 Eberhart［1］于1995年提出的一種優(yōu)化算法.PSO具有較強的全局收斂能力和魯棒性，非常適合于在復雜環(huán)境中對優(yōu)化問題的求解.在解決多目標優(yōu)化問題時，維持著一個較大數(shù)量的種群，高效獲得多個Pareto解集，成為當前研究的熱點問題［2］.采用多目標優(yōu)化可以得到一些目標得到優(yōu)化同時其他目標不惡化的情況，也就是Pareto解集或非支配解集［3］.當?shù)玫絇areto解集或非支配解集之后，決策者們可以根據(jù)實際情況來選擇合理解.在船舶設計領域，商船主要關心的是經(jīng)濟性與節(jié)能，希望得到的是效率最優(yōu)且振動、噪聲等性能不惡化的螺旋槳；軍艦主要關心的是快速性與隱身性，希望得到的是推力與噪聲最優(yōu)且效率減小的螺旋槳.

本文以節(jié)能為目的，選取了DTRC4382型槳為母型，基于螺旋槳水動力性能的面元理論預報程序，采用粒子群算法，螺旋槳的各個幾何參數(shù)用B樣條曲線來擬合，以提高螺旋槳的推力系數(shù)和敞水效率為目標進行螺旋槳優(yōu)化設計，選取Pareto前沿上3個點為方案槳，分析了原槳和方案槳的水動力性能的不同特點.

1 螺旋槳的面元法理論

與螺旋槳的升力線理論和升力面理論不同，螺旋槳的面元法不對物體形狀做任何假設，并在真實物面上滿足邊界條件，能比較精確預報螺旋槳水動力性能，使其在近幾年來得到了廣泛的應用［4］.

假定螺旋槳在速度V0的來流，以角速度ω轉動.利用格林公式將描述不可壓、無粘、無旋的勢流問題的拉普拉斯方程轉化為物體邊界上的積分方程，從而將物體繞流問題的求解轉化為任意物面上未知節(jié)點強度的求解［5］.

式中：Δφ為通過尾渦面的速度勢跳躍，可表示為

在螺旋槳定常問題中，尾渦面的速度勢跳躍Δφ在同一半徑處是常量，由法向偶極子分布與渦環(huán)分布的等價關系可知，Δφ即為尾渦強度.可在升力體的尾緣處滿足庫塔條件以確定該尾渦強度.庫塔條件有多種形式，這里采用壓力庫塔條件，要求在升力體的尾緣處上下表面的壓力差為零，即

2 多目標粒子群算法

2.1 算法的數(shù)學模型

多目標優(yōu)化往往是在一組可行解中選擇Pareto解集的過程，通常多目標優(yōu)化問題包括：n個決策變量，m個目標函數(shù)，k個約束條件，可將其數(shù)學模型表示為

式中：x＝ （x1，x2，…，xn）∈X 代表決策向量.

2.2 算法關鍵問題的處理

盡管基于PSO算法在處理單目標問題時取得了較大的成功，但是并不能直接用于處理多目標問題［6］.目前，MOPSO算法主要采用 NSGA－II非支配排序策略，采用變異操作來保持種群的多樣性，選擇最優(yōu)粒子可采用隨機選擇、輪盤賭或非支配解的疏密等方法［7－9］.

在處理 MOPSO關鍵問題上，本文采用了NSGA－II非支配排序策略來構成Pareto最優(yōu)解集，采用動態(tài)加權法來選擇最優(yōu)粒子來提高Pareto解的多樣性，采用慣性權重來提高全局尋優(yōu)能力.針對種群飛行速度減弱時，采用變異操作來引導粒子跳出局部最優(yōu)解.

NSGA－II中非支配排序思想是目前MOPSO形成Pareto最優(yōu)解的主流方法，操作方法如下，首先將當前支配解集Pt和當前的子代Qt合并成種群Rt，對Rt進行非劣排序和擁擠距離排序得到新的種群Pt＋1.對已形成的Pareto解集，采用動態(tài)加權法來選擇最優(yōu)粒子，按下式動態(tài)計算Pareto解集中各粒子適應度并將當前動態(tài)適應度最大的粒子作為全局最優(yōu)粒子［10］.

粒子群算法模型可被抽象成每個個體由當前位置和速度組成.在每次迭代過程中，為了克服早熟收斂，提高PSO算法開拓能力和保持收斂速度，槳自適應慣性權重和變異操作引入粒子群算法中，當前粒子的速度和速度可按下式更新.

式中：xi為粒子位置；vi為粒子飛行速度；w為慣性權重；r1，r2為［0，1］范圍變化的隨機數(shù)；pi為個體最優(yōu)粒子的位置；pg為全局最優(yōu)粒子位置.k為壓縮因子，表示為對粒子的飛行速度進行約束.

加速因子隨迭代次數(shù)增加而增加，在后期有利于提高全局搜索能力，特別對多峰問題有利.α為加速因子，按下式更新

慣性權重w是衡量上一代的進化速度對現(xiàn)有速度的影響程度，對粒子群算法的收斂性能影響很大.大的慣性權重得到的粒子速度大，有利于擴大搜索空間.小的慣性權重得到的粒子速度小，有利于局部搜索.可按下式調(diào)整.

式中：w0∈ （0，1），w1＞ w0均 為 常 數(shù)；r為 在（0，1）間的隨機數(shù).

本文選擇了4種具有代表性的多目標優(yōu)化測試函數(shù)ZDT1～ZDT4.算法參數(shù)選取：種群規(guī)模N＝100，迭代次數(shù)Nt＝250，算法運行100次，w0＝0.3，w1＝0.8.圖1給出了該算法在測試函數(shù)上所得到的Pareto解.由圖1可知，該算法能較好地保持所得最優(yōu)解在Pareto前端上的收斂性和分布均勻性.

圖1 測試函數(shù)所得Pareto前沿

圖2 推力系數(shù)與敞水效率的散點圖

3 多目標螺旋槳參數(shù)優(yōu)化

本文優(yōu)化設計是以提高推進效率為目標進行的優(yōu)化設計.在明確設計任務的情況下，結合傳統(tǒng)設計確定各參數(shù)的取值范圍，評估各參數(shù)范圍的可行性，提前減少不合理的參數(shù)情況，提高效率.由于設計參數(shù)及其取值范圍的多樣性，設計空間理論上是無窮大的.實際操作中，設計空間的規(guī)模，即方案個數(shù)要適當；方案較少不具有代表性，優(yōu)化效率將會很低，而方案過多則會造成計算成本的大大增加，將造成資源的浪費，也不利于設計效率和設計質(zhì)量的提高.

3.1 優(yōu)化模型和優(yōu)化目標函數(shù)

在設計進速J＝0.899下，以DTR4382槳為母型，目標函數(shù)取螺旋槳推力系數(shù)和敞水效率，但未對其他性能加以限制，目的是為了獲得更多數(shù)據(jù)信息.在實際設計需求中，對螺旋槳其他性能的限定是必不可少的.螺旋槳的螺距分布、側斜分布、縱傾分布、拱度分布和弦長分布可采用B樣條曲線的參數(shù)化表達方式，并根據(jù)這些幾何參數(shù)數(shù)值大小對控制點的變化范圍加限制，以提高搜索效率.對目標槳進行搜索時，將種群數(shù)目設為20，迭代次數(shù)15次，Pareto解的個數(shù)設為20個，對權重的分配與上文的測試函數(shù)選取的相同.

優(yōu)化前的推力系數(shù)和敞水效率的坐標點記為original，而優(yōu)化后的坐標點組成了Pareto前沿，將Pareto前沿解記為Pareto.圖2給出了推力系數(shù)與敞水效率的散點圖.由圖可知，Pareto解的個數(shù)較多且分布較均勻，能夠構成Pareto前端線.original點離Pareto前沿有一定距離.可見，推力系數(shù)或者敞水效率都得到了較大程度的優(yōu)化.該圖也反應了推力系數(shù)和敞水效率是相互矛盾的，推力系數(shù)最優(yōu)必須以降低敞水效率為代價，而敞水效率最優(yōu)就必須以降低推力系數(shù)為代價.

表1 原槳和方案槳的推進系數(shù)和敞水效率

Pareto前沿的每個坐標點代表一個設計方案，船舶設計人員可以根據(jù)不同的船舶需求選擇合適的方案.商船主要關心的是節(jié)能與經(jīng)濟性，可從效率較高的那端選擇.軍艦主要關心的是快速性，可從推力系數(shù)較高的那端選擇.對設計方案進行選擇時，也應對方案進行可行性評估，包括螺旋槳幾何外形光順程度、空泡性能及強度等.為了分析的方便，本文從Pareto前沿上只選擇了3個方案進行分析.3個方案所對應的推進系數(shù)和敞水效率見表1，以下簡稱方案槳.其中，與原槳相比，方案1的敞水效率較高且推進效率不減?。环桨?較為折中，推力系數(shù)與敞水效率都得到相同程度的優(yōu)化；方案3的推進系數(shù)較高且敞水效率不減小.三個方案槳選擇于Pareto前沿的兩端與中部，能夠反映其他方案的特點.

表2給出了原槳和3個方案槳在不同進速下的推力系數(shù)kt和敞水效率ηo的對比.由表可知，方案槳的推力系數(shù)kt和敞水效率ηo都比原槳提高了.在不同進速，方案1能保證敞水效率比其他槳高；方案2比較折中，能保證推力系數(shù)優(yōu)于方案1且敞水效率優(yōu)于方案3；方案3也能保證推進效率比其他槳高.另外，當進速系數(shù)較低時，方案槳的推力系數(shù)明顯比原槳大，而敞水效率相差不大.隨著進速系數(shù)的增大，方案槳的敞水效率較原槳大得多，而推力系數(shù)與原槳越接近.

表2 不同進速下敞水效率與推力系數(shù)對比

3.2 優(yōu)化前后螺旋槳幾何參數(shù)和外形

圖3～5分別給出了原槳和3個方案槳的弦長、螺距及拱度分布，與原槳相比，方案槳的這3個幾何參數(shù)的分布曲線均有改變且能保證光順.3個方案槳的弦長分布變化趨勢相同，不同半徑處的弦長均減小，將導致整個槳葉的盤面比減小，使得螺旋槳剖面上的載荷越大，其吸力面的壓力值越小，壓力面的壓力值越大，該變化將會對螺旋槳的空泡性能產(chǎn)生不利的影響，但是對螺旋槳敞水效率的提高是有利.對于螺距，方案槳的螺距分布比較光順，槳葉葉根附近的螺距減小，而靠近葉梢附近的螺距增加，由螺旋槳環(huán)流理論可以知道，這種變化趨勢將使得螺旋槳負荷徑向分布的最大值位置偏向葉梢，槳葉的效率越高，但葉梢處的空泡性能會變得很差.螺旋槳葉梢附近螺距的增加將增大其葉剖面攻角，進而增大螺旋槳的推力系數(shù)和扭矩系數(shù)，螺旋槳的敞水效率也會明顯增加.對于拱度，方案槳的拱度分布變化較大，3個方案槳在不同半徑處的拱度都增加，而方案3增加最多，對應推力系數(shù)最大，由此推斷拱度與推力系數(shù)有較大的關系；另外，拱度決定葉剖面的負荷沿弦向的分布，拱度增大能避免導邊附近出現(xiàn)負壓峰，推遲空泡的發(fā)生，但對整個槳葉的強度不利.

圖3 弦長分布

圖6 ～7分別給出了原槳和3個方案槳的側斜和縱傾分布.與原槳相比，3個方案槳的側斜分布變化趨勢相同，葉根附近的側斜減小而葉梢附近的側斜增大，槳葉葉梢附近的側斜增加，螺旋槳的水動力性能會提高，但方案槳的側斜分布變化并不明顯，主要是因為側斜對螺旋槳水動力性能的影響不大.對于縱傾，方案槳的槳葉葉根處縱傾減小而葉梢附近的縱傾增加.隨著縱傾的增加，螺旋槳的推力系數(shù)及扭矩系數(shù)會提高，但變動幅度很小，相比其他參數(shù)的影響，縱傾對螺旋槳水動力性能的影響相對較小，而葉梢處的縱傾增大將改善槳葉梢部的三維繞流，進而影響槳葉表面的壓力分布，改變梢渦空泡的初生特性.

圖4 螺距分布

圖5 厚度分布

在進行螺旋槳單參數(shù)的優(yōu)化設計過程中，其他參數(shù)不變，僅改變單個參數(shù)，會造成螺旋槳變化不夠協(xié)調(diào)造成螺旋槳外形的畸變.在螺旋槳多參數(shù)優(yōu)化中，各個參數(shù)之間不能很好的耦合，也會造成螺旋槳外形不光順.螺旋槳的幾何參數(shù)采用怎樣的參數(shù)化表達方式和怎樣進行變形控制，都將影響優(yōu)化后的螺旋槳各幾何參數(shù)的耦合程度.人們希望看到的是各組參數(shù)的變化能耦合在一起，以達到協(xié)同影響螺旋槳性能的目的.圖8給出了原槳和方案槳的模型.由圖8可知，在幾何參數(shù)變化后，3個方案槳都能保持較光順的幾何形狀.可見，本文采用B樣條曲線擬合方式對進行螺旋槳的多參數(shù)優(yōu)化設計，能使得優(yōu)化槳各個參數(shù)之間能夠很好地耦合.與原槳相比，3個方案槳比較直觀的變化是槳葉面積變小，葉梢處的側斜增加.

圖6 側斜分布

圖7 縱傾分布

圖8 原槳與方案槳的模型

4 結 論

1）在設計進速下，優(yōu)化后得到的Pareto解的個數(shù)較多且分布均勻，能夠形成Pareto前沿線，每個Pareto解都可作為方案槳，優(yōu)化后螺旋槳的推力系數(shù)和扭矩系數(shù)都能得到不同程度的提高.

2）在不同進速下，所選的3個方案槳的推力系數(shù)和敞水效率都能比原槳高.方案1能保證敞水效率比其他槳高，方案2能保證推力系數(shù)和敞水效率得到相同程度提高，方案3也能保證推進效率比其他槳高.

3）所選的3個方案槳的各個幾何參數(shù)都出現(xiàn)了變化，從而引起了螺旋槳幾何外形的改變.但是，螺旋槳的幾何參數(shù)采用B樣條曲線表達以及對控制點的變化范圍加以限制，各個參數(shù)之間能夠很好地耦合，能夠形成較光順的螺旋槳幾何外形.

為了獲得更多數(shù)據(jù)信息，未對螺旋槳其他性能加以限制，只對推力系數(shù)和敞水效率的進行尋優(yōu)，優(yōu)化能達到節(jié)能的目的，可能帶來其他性能的惡化.本文為后續(xù)對螺旋槳各方面性能的統(tǒng)籌優(yōu)化設計，以實現(xiàn)螺旋槳的全面優(yōu)化打下基礎.

［1］EBERHART R，KENNEDY J.A new optimizer using particle swarm theory［C］.Proc of the 6th Int’l Symposium on Micro Machine and Human Science.Piscataway，NJ：IEEE Service Center，1995：39－43.

［2］SIERRA M R，COELLO C A C.Multi－objective particle swarm optimizers：A survey of the state－ofthe－art［J］.Int J of Computational Intelligence Research，2006，2（3）：287－308.

［3］TAHER N，MOHAMMAD R N，RASOUL A A，et al.A multi－objective fuzzy adaptive pso algorithm for location of automatic voltage regulators in radial distribution networks［J］.International Journal of Control，Automation，and Systems，2012，10（4）：772－777.

［4］蘇玉民，黃 勝.船舶螺旋槳理論［M］.哈爾濱：哈爾濱工程大學出版社，2003.

［5］王 超.螺旋槳水動力性能、空泡及噪聲性能的數(shù)值預報［D］.哈爾濱：哈爾濱工程大學，2010.

［6］侯 磊.基于多目標粒子群算法的船舶主尺度優(yōu)化設計研究［J］.船舶力學，2011，15（7）：784－790.

［7］KAVEH A，LAKNEJADI K.A novel hybrid charge system search and particle swarm optimization method for multi－objective optimization［J］.Expert Systems with Application，2011，38（12）：475－488.

［8］REZA T M，MOZHGAN A，AZAR S B，et al.Solving a multi－objective job shop scheduling problem with sequence－dependent setup times by a Pareto archive PSO combined with genetic operators and VNS［J］.The International Journal of Advanced Manufacturing Technology，2011，53（5）：733－750.

［9］RENAN S M，MAURO R，VLADIMIRO M，et al.Multi－objective evolutionary particle swarm optimization in the assessment of the impact of distributed generation［J］.Electric Power Systems Research，2012，89：100－108.

［10］陳民鈾，張聰譽，羅辭勇，等.自適應進化多目標粒子群優(yōu)化算法［J］.控制與決策，2009，24（12）：1851－1855.