基于模式函数和变分法的螺旋桨最佳环量计算方法

2019-03-30韩用波周伟新董郑庆董世汤

舰船科学技术 2019年3期

韩用波，周伟新，董郑庆，冯君，董世汤

(1. 中国船舶科学研究中心，江苏无锡 214082；2. 中船重工（上海）节能技术发展有限公司，上海 200011)

0 引言

螺旋桨能够将70%左右的机械动力转化为推进力，成为当前各类船舶上应用最为广泛的推进器形式。对于民船螺旋桨，特别是散货船、油船和集装箱船等三大主流船型，推进效率是这类船舶螺旋桨设计时最优先考虑的指标。在螺旋桨主参数确定后，螺旋桨的效率主要由其径向环量分布决定。开展螺旋桨最佳径向环量的研究可以从机理上给出提高推进器效率的方向。因此以效率为目标的螺旋桨最佳环量计算一直是螺旋桨设计研究的一个重要内容。

螺旋桨最佳环量的计算方法是伴随着升力线理论的发展而提出的。早期，Betz发展螺旋桨升力线理论时，根据在效率较高处增加环量，在效率较低处减小环量的思想，Betz得出了敞水条件下最佳环量分布条件。Lerbs[1]首次提出了求解最佳环量分布的变分法，建立了优化方程式，限于当时的计算条件，未能在计算机上实现。但Lerbs提出的求解最佳环量分布的变分法思想，逐渐成为求解敞开螺旋桨、导管桨、对转桨等各类组合推进桨叶最佳环量的主流方法，一直沿用至今。近年来，国外关于螺旋桨最佳环量计算方面的研究工作，主要集中在麻省理工学院Kerwin教授及其学生身上[2 – 4]。国内苏玉民[5]、曹庆明[6]、孙文愈[7]等也采用变分法开展了船用推进器最佳环量计算的研究工作。

采用传统变分法进行最佳环量分布数值求解时，需要对升力线模型进行离散，将连续的环量分布离散为在涡格内均布的常量，连续的自由涡片离散为从涡格边界泄出的集中涡。根据Kerwin[8]的研究成果，在使用离散化升力线模型时，在叶梢和叶根2个自由端存在1/2阶的奇性，因此需要对2个端部的网格划分和控制点布置进行特殊处理。本文给出一种求解螺旋桨最佳环量的改良方法，将径向环量分布用正弦级数表达出来，通过求解正弦级数模式函数的各项系数，从而得到螺旋桨最佳环量分布形式。

1 基本理论

1.1 正弦级数模式函数表达螺旋桨径向环量分布

1.2 螺旋桨诱导速度的计算

诱导因子仍按Lerbs的定义不变，则诱导速度的计算式如下：

由式（1）可得：

可以将式（4）写成如下形式：

其中：

1.3 功率系数和推力系数的计算

将式（1）中的环量分布改写成如下形式：

功率系数

令

改写

又令

则功率系数计算式可与成以下形式：

同理，推力系数可以写成如下形式：

其中：

1.4 变分法求解最佳环量

求解最佳环量分布的条件是：推力系数须滿足指定值而功率系数为最小，或功率系数满足指定值而推力系数最大，这一问题可转化为多元函数求条件极值的问题。以指定推力系数求的极小值为例，引入Lagrange乘子，并构建辅助函数，则求极值条件的方程式为：

其中：

联立式（16）~式（18）可得：

整理后可写成：

但是，这一轮的过程是建立在某一给定的水动力螺距角分布的基础上的，故需有迭代过程。在第一轮求解时可以先用进角分布代替水动力螺距角作为初值。在求出表达环量分布的正弦级数各项系数后，重新计算诱导速度和水动力螺距角分布，用新的水动力螺距角分布再重新计算得出新的最佳环量分布，一直迭代到和诱导速度收敛。至此方求出表达最佳环量的正弦级数各项系数的最终值后，代入式（1）最佳环量分布得解。