4. 船舶阻力与推进需求分析

在为一条目标船选配螺旋桨之前,必须先回答一个基本问题:船在给定航速下究竟需要多大的推力?这个问题的答案来自船体阻力曲线,而螺旋桨的任务就是把电机输出的功率转换成恰好等于该阻力的推力。本节围绕“阻力—有效功率—推力需求—螺旋桨工作点”这条主线,说明如何用工程估算方法把船体需求换算到螺旋桨轴上,并初步确定目标船型的螺旋桨工作点。

一、船体阻力曲线的构成与读法

船舶总阻力通常分解为摩擦阻力、剩余阻力(含兴波阻力与粘压阻力)以及附体阻力、空气阻力等附加项。工程上常用无量纲的阻力系数表达:

R_t = 0.5 · ρ · V² · S · C_t
C_t = C_f + C_r + ΔC_a
C_f = 0.075 / (lg Re − 2)²      (ITTC-1957 摩擦阻力公式)
Re = V · L / ν

其中 ρ 为水的密度(海水约 1025 kg/m³),V 为航速(m/s),S 为湿表面积,C_f 为摩擦阻力系数,C_r 为剩余阻力系数,ΔC_a 为粗糙度与附体修正。对一条给定船,S、L 是几何常量,C_r 随傅汝德数 Fr = V/√(gL) 变化,因此总阻力曲线 R_t(V) 是一条随航速上升而加速增长的曲线——低速段近似平方增长,进入兴波显著区后增长更快。

实际估算时,通常以母型船或系列图谱(如泰勒系列、BSRA 系列)的 C_r–Fr 曲线为输入,逐点计算 3~5 个航速下的 R_t,再拟合成曲线。下面给出一个教学示例(假设值,非实测):某内河货船 L=68 m,B=12.6 m,T=2.8 m,S≈620 m²,设计航速 12 kn(6.17 m/s),海水密度取 1025 kg/m³,ν=1.19×10⁻⁶ m²/s。

航速 V (kn)V (m/s)Re (×10⁸)C_f (×10⁻³)C_r (×10⁻³)R_t (kN)
84.122.351.830.6238.5
105.142.941.760.9562.0
126.173.531.711.4898.6
147.204.111.672.35158.2

可以看到,航速从 12 kn 提高到 14 kn,阻力从约 98.6 kN 跃升到 158.2 kN,增幅超过 60%,而航速只提高 16.7%。这正是“阻力曲线陡升段”的典型特征,也是螺旋桨与电机匹配时最容易出问题的区域。

二、有效功率与推力需求的换算

有效功率 P_E 定义为拖曳船体所需功率:

P_E = R_t · V        (单位:W,V 用 m/s)

以 12 kn 工况为例,P_E = 98.6 kN × 6.17 m/s ≈ 608 kW。但 P_E 只是“船体需要的功率”,螺旋桨轴端需要提供的功率还要除以推进效率:

P_D = P_E / η_D
η_D = η_O · η_H · η_R
T = R_t / (1 − t)

其中 η_D 为总推进效率,η_O 为敞水效率,η_H 为船身效率,η_R 为相对旋转效率;t 为推力减额分数。对上述内河货船,取 t=0.18、w=0.22(伴流分数),则船身效率 η_H=(1−t)/(1−w)=0.82/0.78≈1.05;若 η_R≈1.0、η_O≈0.55,则 η_D≈0.58。于是轴功率 P_D≈608/0.58≈1048 kW,所需推力 T=98.6/0.82≈120.2 kN。

这一步换算的意义在于:船体阻力曲线给出的是“需求”,螺旋桨推力曲线给出的是“供给”,两者必须在同一航速下相交。注意推力 T 与阻力 R_t 之间差一个推力减额系数,不能直接画在同一坐标上比较,必须先换算到同一基准。

三、螺旋桨工作点的确定

螺旋桨的推力与扭矩随进速系数 J = V_A/(nD) 变化,其中 V_A = V(1−w) 为进速,n 为转速,D 为直径。对给定螺旋桨,其敞水特性 T=K_T·ρn²D⁴、Q=K_Q·ρn²D⁵ 是已知曲线。把船体需求换算成“螺旋桨必须提供的推力”后,工作点就是推力需求曲线与螺旋桨推力曲线的交点,同时还要满足扭矩平衡:

Q_P = P_D / (2πn)
K_Q = Q_P / (ρ n² D⁵)

工程上更常用的做法是绘制“螺旋桨负载线”:以 n 为横坐标,画出螺旋桨所需扭矩 Q_P(n) 与电机可输出扭矩 T_m(n) 两条曲线,交点即稳态工作点。对上述示例,若选 D=1.6 m、四叶、盘面比 0.55 的定距桨,在 12 kn 时取 n=280 r/min(4.67 r/s),V_A=6.17×0.78≈4.81 m/s,J=4.81/(4.67×1.6)≈0.644。查该桨敞水图谱得 K_T≈0.20、K_Q≈0.032,则:

T = 0.20 × 1025 × 4.67² × 1.6⁴ ≈ 120.6 kN
Q = 0.032 × 1025 × 4.67² × 1.6⁵ ≈ 3.06 kN·m
P_D = 2π × 4.67 × 3.06 ≈ 89.8 kW?

此处出现明显矛盾:按推力算出的 P_D 应约 1048 kW,而按扭矩算只有约 90 kW。这说明所取 K_Q 或转速不自洽——实际计算中必须用同一组 (J, K_T, K_Q) 同时校核推力和扭矩,若不一致则调整 n 或 D 重新迭代。这正是工作点估算中最容易出错、也最需要反复校核的环节。正确的做法是:先由推力需求反算所需 n,再由该 n 下的 K_Q 求扭矩与功率,最后与电机额定点比较。

目标船阻力曲线与螺旋桨工作点估算 航速 V (kn) 阻力 R_t (kN) 0 8 10 12 14 0 50 100 150 200 R_t(V) 总阻力曲线 T = R_t/(1−t) 推力需求 设计工作点 12 kn T≈120 kN,P_E≈608 kW 阻力陡升区 12→14 kn 阻力 +60% 电机易过载 船体阻力 R_t C_f + C_r + ΔC_a P_E = R_t·V 推力需求换算 T = R_t/(1−t) P_D = P_E/η_D 螺旋桨工作点 J = V_A/(nD) K_T、K_Q 校核

四、估算结果如何判断与使用

得到工作点后,需要做三项判断。第一,功率裕度:电机额定功率应大于 P_D 并留 10%~15% 裕度,以覆盖污底、风浪、浅水等阻力增值。第二,转速匹配:工作点转速应落在电机恒扭矩区与恒功率区的合理位置,避免长期在弱磁高速区运行导致扭矩不足。第三,重载校核:把阻力曲线整体上移 15%~20% 模拟污底工况,检查螺旋桨负载线是否越过电机扭矩限值——若越过,说明该桨偏大或减速比偏小,需要减小 D 或调整传动比。

估算阶段的所有系数(t、w、η_O、ΔC_a)都应注明来源:来自母型船、系列图谱还是经验取值。同一份计算书中系数口径必须统一,否则推力与扭矩校核会出现本节示例中那种“功率对不上”的矛盾。

最后要强调,阻力曲线与工作点估算只是匹配设计的起点。实船试航时应用轴功率仪和转速传感器记录实际 P_D–n 关系,与估算曲线对比;若实测功率明显高于估算值,优先排查污底、螺旋桨加工偏差和轴系损失,而不是直接加大电机功率。把估算、选型、试航验证串成闭环,螺旋桨与电机的匹配才算真正落地。