武 璇，仉永超，靳良真，李廣輝

(濰柴重機(jī)股份有限公司，山東 濰坊 261108)

0 引言

基于環(huán)流理論設(shè)計(jì)螺旋槳的過(guò)程中，要使螺旋槳的性能達(dá)到最佳狀態(tài)，應(yīng)合理選擇螺旋槳各參數(shù)的分布形式。影響螺旋槳性能的參數(shù)主要有槳轂形狀、盤面比、側(cè)斜和縱傾等，在螺旋槳的優(yōu)化設(shè)計(jì)中必須首先知道這些參數(shù)對(duì)螺旋槳水動(dòng)力性能和空泡性能的影響。

王超等基于速度勢(shì)面元法研究了側(cè)斜與縱傾對(duì)螺旋槳水動(dòng)力性能的影響。胡健等基于CFD方法，探討了縱傾螺旋槳敞水性能和空化性能的影響。王詩(shī)洋等基于速度勢(shì)面元法，研究了螺旋槳幾何參數(shù)對(duì)其水動(dòng)力性能的影響。以上文獻(xiàn)研究了螺旋槳幾何參數(shù)，并取得一系列成果，但對(duì)葉切面參數(shù)對(duì)螺旋槳水動(dòng)力性能的影響研究較少。為此，本文以KP505螺旋槳為原型，借助CFD軟件對(duì)螺旋槳的水動(dòng)力性能進(jìn)行數(shù)值模擬，并探討不同螺距、拱度的徑向分布對(duì)螺旋槳水動(dòng)力性能的影響，總結(jié)各性能參數(shù)的變化規(guī)律。

1 控制方程

在螺旋槳水動(dòng)力性能數(shù)值模擬過(guò)程中，流體均視為不可壓縮流體，因此基本控制方程包括連續(xù)性方程和動(dòng)量守恒方程，如下所示：

(1)

(2)

(3)

(4)

式中：

μ

為湍流粘度；

G

為由平均速度梯度引起的湍動(dòng)能

k

的生成項(xiàng)；

C

、

σ

、

σ

為常數(shù)，分別等于1.9、1.0、1.2。

2 幾何模型建立

本文以KP505螺旋槳為研究對(duì)象，其模型主要參數(shù)如下：直徑0.25 m，葉數(shù)5，盤面比0.8，側(cè)斜角32°，螺距比0.996 7，翼型NACA66。

根據(jù)文獻(xiàn)[7]提供的螺旋槳的幾何參數(shù)，且已知NACA66翼型的拱度、厚度沿弦長(zhǎng)的分布形式，便可以通過(guò)計(jì)算得到各切面二維坐標(biāo)；然后根據(jù)坐標(biāo)轉(zhuǎn)換公式，通過(guò)Fortran自編代碼計(jì)算槳葉離散點(diǎn)三維坐標(biāo)(其中：

x

軸指向螺旋槳后方，

y

軸垂直向上，

z

軸按照右手定則)；最后將槳葉的三維坐標(biāo)點(diǎn)導(dǎo)入U(xiǎn)G軟件中，通過(guò)三維坐標(biāo)點(diǎn)生成螺旋槳曲面，并進(jìn)行曲面縫合、實(shí)體陣列生成實(shí)體三維模型?？紤]槳轂對(duì)推進(jìn)器水動(dòng)力性能的影響較小，本文采用假定的槳轂進(jìn)行三維建模。在三維軟件中建立的KP505螺旋槳空間坐標(biāo)點(diǎn)及三維幾何模型見(jiàn)圖1。

圖1 KP505螺旋槳空間坐標(biāo)點(diǎn)和三維幾何模型

3 螺旋槳水動(dòng)力性能數(shù)值模擬

3.1 計(jì)算域與邊界條件設(shè)置

本文采用STAR-CCM+軟件對(duì)螺旋槳水動(dòng)力性能進(jìn)行數(shù)值模擬時(shí)，流場(chǎng)計(jì)算域的大小劃分見(jiàn)圖2。將計(jì)算域劃分為靜止域和旋轉(zhuǎn)域，且均為圓柱體。為了更好地模擬螺旋槳復(fù)雜的尾流流動(dòng)，本文設(shè)置速度入口與槳盤面的距離為5

D

(

D

為螺旋槳直徑)，壓力出口與槳盤面距離為13

D

，靜止區(qū)域半徑為3

.

5

D

，旋轉(zhuǎn)區(qū)域的軸向距離為0

.

6

D

，旋轉(zhuǎn)區(qū)域直徑為1

.

2

D

。對(duì)該區(qū)域采用動(dòng)參考系模型(Moving Reference Frame，MRF)方法設(shè)置旋轉(zhuǎn)坐標(biāo)系，實(shí)現(xiàn)螺旋槳的旋轉(zhuǎn)運(yùn)動(dòng)。

圖2 計(jì)算域劃分示意圖

邊界條件的設(shè)置：計(jì)算域軸向上游邊界設(shè)為速度入口，以給定進(jìn)口邊界上各節(jié)點(diǎn)的速度值；軸向下游邊界設(shè)為壓力出口；槳葉、槳轂表面均設(shè)為無(wú)滑移壁面；靜止域周向邊界設(shè)為對(duì)稱平面；靜止域和旋轉(zhuǎn)域相互接觸界面均設(shè)置為交界面以傳遞數(shù)據(jù)信息，從而保證流動(dòng)的連續(xù)性。

3.2 網(wǎng)格劃分與無(wú)關(guān)性驗(yàn)證

圖3 計(jì)算域軸向網(wǎng)格劃分示意圖

表1 3種不同密度網(wǎng)格數(shù)量

不同網(wǎng)格密度進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證時(shí)，螺旋槳轉(zhuǎn)速取

n

=20 rad/s，且在計(jì)算中保持轉(zhuǎn)速不變，通過(guò)改變?nèi)肟谒俣戎祦?lái)得到不同的進(jìn)速系數(shù)。不同網(wǎng)格密度的螺旋槳水動(dòng)力性能計(jì)算結(jié)果與試驗(yàn)值誤差見(jiàn)表2(KP505敞水試驗(yàn)數(shù)據(jù)來(lái)自文獻(xiàn)[9])。由表2可知：在進(jìn)速系數(shù)

J

=0.3～0.7工況下，3種網(wǎng)格密度計(jì)算所得到的結(jié)果均與試驗(yàn)值吻合較好，誤差均在5%以內(nèi)，僅粗密度網(wǎng)格所對(duì)應(yīng)的計(jì)算結(jié)果與試驗(yàn)值誤差稍微偏大；而中密度與細(xì)密度網(wǎng)格所對(duì)應(yīng)的計(jì)算結(jié)果與試驗(yàn)值誤差并無(wú)明顯減小，反而出現(xiàn)波動(dòng)現(xiàn)象。在

J

=0.9工況下，3種網(wǎng)格密度計(jì)算所得到的值均出現(xiàn)誤差較大現(xiàn)象。這是因?yàn)樵诼菪龢ㄋ阅苡?jì)算過(guò)程中，所采用的湍流模型均基于湍流充分發(fā)展的前提下，而對(duì)于螺旋槳槳模敞水試驗(yàn)中存在層流區(qū)域和過(guò)渡區(qū)域，因此采用充分發(fā)展的湍流模型會(huì)導(dǎo)致螺旋槳敞水?dāng)?shù)值模擬誤差偏大。但整體而言，隨網(wǎng)格密度的增加，其數(shù)值計(jì)算精度并未明顯改善。由于計(jì)算資源的限制，本文選用中密度網(wǎng)格進(jìn)行后續(xù)計(jì)算研究。

表2 不同網(wǎng)格密度下KP505螺旋槳水動(dòng)力性能計(jì)算結(jié)果

3.3 不同螺距與拱度的徑向分布

為討論螺距、拱度分別對(duì)螺旋槳水動(dòng)力性能的影響，在其他參數(shù)不變的情況下，只改變螺旋槳螺距、拱度的徑向分布。為了便于分析與比較，在原母型槳螺距、拱度的基礎(chǔ)上，分別選取了0.9倍和1.1倍徑向螺距分布0.6倍和1.4倍徑向拱度分布。螺距、拱度沿徑向的分布形式見(jiàn)圖4。

r—螺旋槳任意處的半徑；R—螺旋槳半徑。

3.4 螺距與拱度對(duì)螺旋槳水動(dòng)力性能影響分析

J

=0.1～1.1時(shí)不同螺距下對(duì)應(yīng)的敞水曲線見(jiàn)圖5。從計(jì)算的敞水曲線結(jié)果可以看出，

K

、

K

和

η

隨

J

的變化規(guī)律一致且與常規(guī)圖譜槳的敞水圖譜隨螺距的變化規(guī)律一致。隨著螺旋槳的螺距增大，其推力和轉(zhuǎn)矩系數(shù)均明顯增大。另外，在螺旋槳盤面比不變的情況下，增加螺距會(huì)導(dǎo)致螺旋槳單位面積的平均推力增大，使得葉元體最大減壓系數(shù)增大且更易產(chǎn)生空泡現(xiàn)象。而對(duì)于敞水效率而言，在中、低進(jìn)速系數(shù)下，螺距越小，效率越高，其最高效率點(diǎn)隨螺距的增大所對(duì)應(yīng)的

J

也隨之增大。

圖5 不同螺距下螺旋槳敞水曲線對(duì)比

J

=0.1～0.9時(shí)不同拱度所對(duì)應(yīng)的敞水曲線見(jiàn)圖6。從計(jì)算的敞水曲線結(jié)果可以看出：

K

、

K

和

η

隨進(jìn)速系數(shù)

J

的變化規(guī)律一致。另外，隨著螺旋槳的拱度增大，其推力和轉(zhuǎn)矩系數(shù)均明顯增大；而對(duì)于敞水效率而言，在最高效率點(diǎn)之前，效率隨著拱度的增加而降低，其效率值變化幅度較小。

圖6 不同拱度下螺旋槳敞水曲線對(duì)比

從

J

=0.4時(shí)不同拱度下槳葉表面壓力分布可知：隨著葉切面拱度的增大，導(dǎo)致螺旋槳槳葉表面壓力特性的變化，其葉面的壓力沿弦長(zhǎng)方向的變化較為明顯，葉面、葉背的高壓區(qū)域由隨邊向?qū)н呥^(guò)渡，導(dǎo)致螺旋槳葉面、葉背的壓力差增大，所以導(dǎo)致螺旋槳推力的增大。另一方面，隨著葉切面拱度的增大，槳葉切面上的水流速度過(guò)快,葉背壓力明顯下降,使得葉背更易達(dá)到水的汽化壓力臨界值，進(jìn)而導(dǎo)致葉背區(qū)域的空泡出現(xiàn)更早且面積更大；槳單位面積上承受壓力過(guò)大引起的空泡，除了會(huì)產(chǎn)生螺旋槳的空泡剝蝕外,同時(shí)會(huì)誘導(dǎo)槳產(chǎn)生空化噪聲。所以在螺旋槳設(shè)計(jì)優(yōu)化時(shí)應(yīng)選擇合適的拱度比，以設(shè)計(jì)出性能最佳的螺旋槳。

4 結(jié)論

(1)通過(guò)網(wǎng)格無(wú)關(guān)性驗(yàn)證，得到計(jì)算誤差基本在5%以下，驗(yàn)證了本文所建立的幾何模型、數(shù)值模型對(duì)數(shù)值求解的準(zhǔn)確性。

(2)盤面比不變時(shí)，隨螺距、拱度的增大，螺旋槳推力顯著增大，其單位面積上的平均推力增大，導(dǎo)致葉切面的最大減壓系數(shù)增大且更易產(chǎn)生空泡。

(3)在中、低進(jìn)速系數(shù)下，螺距越小，敞水效率越高；而隨拱度的增大，敞水效率無(wú)明顯變化。