螺旋桨三维建模与水动力数值分析

0.5022 1.163 0.4296 1.018 0.3873 0.9375 0.2996 0.7592
计算所得到的敞水性能曲线及与试验计算结果的 比较如图 4 所示.图中，螺旋桨水动力参数的计算结果
与试验结果能够较好的吻合，KT，KQ的误差最大不超 过 3.13%和 6.5%.同时我们也发现随着螺旋桨载荷的 增加（进速系数减小），计算值与试验值偏差有增加的 趋势，这种现象在使用RANS数值模拟中普遍存在[4]. 存在这些差异的原因主要是试验环境中的管道壁面和 桨毂轮廓的影响，进流速度的不均匀，这些在CFD中 模拟不能完全吻合.与试验结果相比,计算所得的KT-J 线斜率绝对值在 0.5≤J≤1.1 中稍偏大，而KQ-J线计算 结果与试验数据相比在小进速下误差稍大，其中最大 偏差在 6.3%左右，在J＝1.55 附近两者相交，随后计 算值比试验值略小，偏差在 5.48%.小进速下KQ-J误差 大的原因有两种可能：1）在数值模拟中没有使用空化 模型，空泡现象的存在可能对螺旋桨性能有影响；2） 进速系数小的情况下螺旋桨周围流场没有完全发展成 湍流.对于η0-J（η0为旋转效率）曲线，当J≤0.8 时计 算值与试验值几乎完全吻合，之后两者之间的偏差有 所增加，为负偏差，在J＝1.6 时的偏差最大达到 7.41%.
关键词：船舶；船舶推进装置；螺旋桨；三维建模；计算流体力学方法 中图分类号：U664.33 文献标识码：A 文章编号：1000-6982 (2008) 06-0023-04
3D modeling and numerical analysis for hydrodynamic force of propeller
OXYZ平行. O′ X1Y1Z1坐标系的规定如图 1（b）中所 示. O′ X1Y1Z1坐标系可以通过一次旋转与 O′X ′Y ′Z ′ 坐 标系重合，由此可得：
⎛ ⎜ ⎜
X1 Y1
⎞ ⎟ ⎟
=
⎛ ⎜ ⎜
1 0
0 cosϕ
0 ⎞⎛ X′⎞
sin
ϕ
⎟ ⎟
⎜ ⎜
Y
′
⎟ ⎟
⎜⎝ Z1 ⎟⎠ ⎜⎝ 0 − sinϕ cosϕ ⎟⎠ ⎜⎝ Z ′ ⎟⎠
总的来说计算结果和试验误差在工程所允许的范 围内，利用 CFD 技术可以对螺旋桨的敞水性能做出较 为可靠的预报，如何进一步提高计算精度值得深究.
KT, KQ, η0
1.3
1.2
试验
1.1
1.0
10KQ
数值
0.9
0.8
0.7 0.6
η0
0.5
0.4
0.3
KT
0.2
0.1
0 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6
3 数值计算结果的考察分析
3.1 敞水性能曲线计算结果与试验对比
本次数值计算选用的螺旋桨直径 0.25m，其转速n 为 900r/min，其对应的雷诺数为 9.33×105，在这里通 过固定转速，改变进速的办法来改变进速系数J.进速
系数J（J=VA/nD）的变化范围为 0.5~1.6，其中，VA为 水的进速，m/s；n为螺旋桨的旋转速度，r/min.
图 2 螺旋桨的三维实体模型
2 数值计算过程
2.1 计算域与网格划分
鉴于来流的均匀和螺旋桨几何上的周期性，为进
行敞水试验只需取单个桨叶所在的单通道作为计算域
即可分析，由此节省计算时间提高计算效率.单桨叶通
道内外边界定在桨毂和直径为螺旋桨直径数倍的圆柱
体表面.进出口面为 90°的扇面.在近壁面区，流动情况
0 引言 获得螺旋桨水动力性能数据主要通过实验和计算
的方法.目前广泛采用的实验方法是螺旋桨模型敞水 试验，它是检定和分析螺旋桨性能较为简单的方法. 近年来随着实验条件的进步，出现了激光多普勒仪 （LDV），它的不干扰流场和高精度的特性被作为一种 先进的流场测量手段，可是这需要非常精密的试验设 备，同时，每个桨模的制作周期长，费用昂贵，不具 有快速预报螺旋桨敞水性能的能力.从 1912 年茹可夫 斯基提出螺旋桨旋涡理论以来，螺旋桨计算方法发展 已有数十年，基于势流理论而建立升力面理论，面元 法等[1]已有相当水平并广泛应用到螺旋桨的设计与预 报当中.然而，这些方法的计算过程异常烦琐，需要事 先掌握许多预备知识.N-S方程考虑了流体的粘性影
表 2 给出使用k-ε 模型在J=0.6～0.9 时三种不同的 网格方案下的KT，KQ（分别为推力系数和扭矩系数） 的计算结果.三种网格形式对KT，KQ的计算结果的影 响不是很明显，不超过 1.3%，故使用“Base”方案来 分析以下的计算结果.
表 2 不同网格形式下的KT，KQ结果比较
J
Base Extended Fine
YAO Zhen-qiu, GAO Hui, YANG Chun-lei
(College of Naval Architecture & Ocean Engineering, Jiangsu University of Science and Technology, Zhengjiang 212003, China)
等数据绘制正投影图﹑侧投影图﹑伸张轮廓图；而对
于三维桨建模来说，则是通过螺旋桨基本参数和各叶
切面二维形状尺寸建立与螺旋桨曲面所有型值点空间
坐标的关系式.图 1 中OH是基线.θ为纵斜角，ϕ为螺距
角.全局坐标系OXYZ的OXY平面与螺旋桨轮毂端面平
行 O′ 为基线与圆柱面的交点，坐标系 O′X ′Y ′Z ′ 与
J=0.6 J=0.8 J=0.9 J=1.1
KT 10KQ
KT 10KQ
KT 10KQ
KT 10KQ
0.5024 1.167 0.4283 1.016 0.3870 0.9350 0.2997 0.7582
0.5023 1.167 0.4237 1.013 0.3868 0.9386 0.2998 0.7581
响，对于因粘性而产生的桨叶表面边界层的生成、发 展、分离桨叶梢涡的形成都有可能预报[2]，但早期的 解决方法需要大量的计算，制约了它的应用.现在随着 计算机的迅猛发展，已经具备了几年以前还难以想象 的计算能力，使得一些比较成熟的船舶性能计算方法， 开始在实践中发挥重要的作用.
本文利用螺旋桨基本参数建立螺旋桨三维模型并 采用 Fluent 软件对其敞水性能和流场进行数值仿真计 算，给出了敞水性征曲线，并就计算结果与试验测量 值作了比较和分析.
（1）
由式（1）得
⎛ X′⎞ ⎛1
⎜ ⎜
Y
′
⎟ ⎟
=
⎜ ⎜
0
0 cosϕ
−
0 sin
ϕ
⎞ ⎟ ⎟
⋅
⎛ ⎜ ⎜
X1 Y1
⎞ ⎟ ⎟
⎜⎝ Z ′ ⎟⎠ ⎜⎝ 0 sinϕ cosϕ ⎟⎠ ⎜⎝ Z1 ⎟⎠
（2）
由图1（a）可得到坐标变换关系：
⎛ ⎜ ⎜
X Y
⎞ ⎟ ⎟
=
⎛ ⎜ ⎜
X
′
+
Ri Y
cos ′
Ψ
⎞ ⎟ ⎟
⎜⎝ Z ⎟⎠ ⎝⎜ Z ′ − Ri tanθ ⎠⎟
（3）
其中， Ψ = Y ′ R
由式（2）~式（3）整理可得：
ቤተ መጻሕፍቲ ባይዱ
⎛
⎛ ⎜ ⎜⎜⎝
X Y Z
⎞ ⎟ ⎟⎟⎠
=
⎜ ⎜ ⎜ ⎜ ⎜
Ri
cos
Y1
cosϕ
− Ri
Z1
sin
ϕ
Ri
sin
Y1
cosϕ
− Ri
Z1
sin
ϕ
⎞ ⎟ ⎟ ⎟ ⎟ ⎟
⎜⎝Y1 sinϕ + Z1 cosϕ − Ri tanθ ⎟⎠
1 螺旋桨三维模型的建立 1.1 螺旋桨曲面型值的坐标变换方法
传统的螺旋桨手工作图法的几何视图表达方式是
收稿日期：2007-11-05；修回日期：2008-02-22 作者简介：姚震球（1966-），男，副教授，博士，主要从事船舶设计研究.
— 23 —
万方数据
根据螺旋桨提供的桨叶轮廓尺寸表和桨叶切面尺寸表
变化很大，特别是粘性底层，几乎是层流，湍流应力
几乎不起作用，所以采用低 Re 数 k-ε 模型求解粘性影
响比较明显的区域，这就要求在壁面区划分比较密的
网格，越靠近壁面，网格越细.沿桨叶面以及桨毂表面
的法向方向按一定的比例伸出 6 层边界层网格，根据
经验第一层高度定义为 0.001D，其它区域都使用四面
体网格.为了检验出口面距离以及不同网格形式对计
算结果的影响，建立三种方案，如表 1 所示，面网格
划分见图 3.
表 1 网格划分方案
标记 出口到原 螺旋桨面网格形式 网格
点距离
数目
导边、随边、叶梢、
叶根部附近区域采
Base
0.72D 用长度为 0.0025D
492992
的三角形面网格，
其余长度为 0.005D
（4）
式（4）即为曲面型值点的局部坐标到全局坐标的
转换公式.
Z
θO Y
XM
Z1 Z′
Y1
Y′ O′ X1
X′
ϕ
ψ Ri
B AH (a)
Z1
M
Z′ ϕ Y1
O′ Y′
ϕ
A′ B′
(b)
图 1 螺旋桨投影原理图
1.2 螺旋桨曲面型值点坐标进行计算 本文研究的螺旋桨为四叶的MAU桨，其盘面比
— 24 —
0.55，螺距比 1.6，毂径比 0.18，纵倾角 10°.根据各叶 切面几何参数，由式（4），通过计算程序就可以方便 地得到螺旋桨所有曲面型值点的坐标.由于桨叶轮廓 尺寸表和桨叶切面尺寸表只提供了 0.2R~0.95R的型 值，所以必须对桨叶叶根处插值进行修正.为了计算方 便需要将桨作一些处理：桨毂用无限长的圆柱体来代 替，导边随边与桨毂表面用光顺的曲线连接.在 pro/ENGINEER中得到光滑曲面[3]，最终生成的螺旋桨 实体模型如图 2 所示，其中，Z轴与螺旋桨的旋转轴方 向一致，以船艏方向为正向，Y轴与桨叶参考线一致， X轴服从右手法则.
Abstract: The coordinate transformation formula for transforming the local to the global coordinate is deduced and introduced the 3D modeling at the same time. Dispose the propeller locally for the convenience of numerical analysis and provide the program. Further more, the CFD method is used to simulate numerically hydrodynamic characteristics of MAU propeller in order to have the capability of quick prediction of propeller open water performance, then compared with the measurement. The relevant investigation on the open water behavior of this propeller is also inspected. Key words: ship; marine propulsion; propeller; 3D modeling; CFD
Extended
2D
同上
647512
长度为 0.0025D 的
Fine
2D
三角形面网格
1635235
万方数据
(a) Base(Extended)
(b) Fine
图 3 螺旋桨面网格
2.2 边界条件 速度进口边界给定均匀来流的各个速度分量；出
口边界给定表压为 0（相对于工作压强）；其上远扬外 边界同样设为速度进口边界；叶片和桨毂定义为固体 壁面，无滑移条件；整个桨模的中心与坐标原点重合； 单桨叶通道的法线夹角为 90°的两个周向侧面设置为 旋转周期性边界.计算域流体则按单旋转坐标系，设置 绕 z 轴以角速度 n 旋转.
SHIP ENGINEERING Vol.30 No.6 2008
船舶工程 总第 30 卷，2008 年第 6 期
螺旋桨三维建模与水动力数值分析
姚震球，高 慧，杨春蕾
（江苏科技大学 船舶与海洋工程学院，镇江 212003）
摘 要：推导了螺旋桨叶切面局部坐标系到全局坐标系的坐标转换公式，给出了三维实体建模过程. 为数值计算方便，对桨模做了一些局部处理.然后运用计算流体力学方法（CFD）对螺旋桨的水动力特性 进行数值模拟，以尽快形成螺旋桨敞水性能 CFD 计算的快速预报的能力.文章以 MAU 型桨作为研究对 象，给出其敞水性能的数值计算结果并与试验值做了比较，获得良好的结果.同时还对该桨型周围流场进 行了一些考察.