📚 课程目录第4章船舶流体动力与参数获取

4. 船舶流体动力与参数获取

在建立船舶操纵运动数学模型时，船体受到的水动力是决定仿真可信度的核心。工程上通常把船体水动力分解为惯性项与阻尼项：惯性项由附加质量描述，阻尼项由阻尼系数描述；再叠加舵力与螺旋桨推力，构成完整的六自由度（工程仿真中常简化为纵荡、横荡、艏摇三自由度）运动方程。本节围绕这四类参数——附加质量、阻尼系数、舵力模型、螺旋桨推力模型——说明它们的物理含义、获取途径与典型估算方法，并给出可用于MATLAB仿真的参数组织方式。

坐标系与符号约定

为避免正负号混乱，先统一约定：采用船体坐标系，原点取在船舶重心（或设计水线面中心），x轴指向船首，y轴指向右舷，z轴指向龙骨下方；纵荡速度u、横荡速度v、艏摇角速度r均以该坐标系为正方向。舵角δ以右舵为正，螺旋桨转速n以产生正推力（前进方向）为正。后文所有公式与图示均遵循这一约定。

附加质量：惯性项的来源

船体在流体中加速时，必须同时带动周围流体一起加速，这部分“被带动的流体质量”即附加质量。三自由度模型中主要关注三个量：纵荡附加质量 \(m_x\)、横荡附加质量 \(m_y\) 与艏摇附加惯量 \(J_z\)。它们与船体质量同量纲，但数值差异很大：细长船体的 \(m_x\) 通常只有船体质量的百分之几到十几，而 \(m_y\) 可达船体质量的0.5～1.5倍，艏摇附加惯量则常与船体自身转动惯量同量级。

获取附加质量有三条途径。其一是势流理论计算，如基于细长体理论的切片法，适合方案设计阶段；其二是经验公式，例如以船长L、吃水T、方形系数Cb为自变量的回归式；其三是从CFD或约束模型试验的加速度响应中辨识。教学仿真中可采用如下量级估算（示例假设，非实测）：

% 典型内河货船示例参数（假设值，仅用于教学演示）
L = 80; B = 13; T = 4.2; Cb = 0.82;
rho = 1025;            % 海水密度 kg/m^3
m   = rho * L*B*T*Cb;  % 船体质量
mx  = 0.08 * m;        % 纵荡附加质量
my  = 0.90 * m;        % 横荡附加质量
Jz  = 0.05 * m * L^2;  % 艏摇附加惯量（量级估算）

需要强调，这些系数是量级参考，实际取值应结合船型资料或辨识结果修正；在仿真中把它们单独定义为参数变量，便于后续灵敏度分析。

阻尼系数：速度相关的阻力

阻尼项描述船体运动受到的水动力阻力，通常写成速度的多项式。纵荡方向以二次阻力为主：\(X = -X_{uu}u|u|\)；横荡与艏摇方向在小扰动下可线性化，用 \(Y_v\)、\(Y_r\)、\(N_v\)、\(N_r\) 四个水动力导数描述，它们分别表示单位横荡速度或艏摇角速度引起的横向力与艏摇力矩。工程上常把阻尼写成“线性+非线性”组合，兼顾低速操纵与高速直航两种工况。

阻尼系数的获取以约束模型试验（PMM、平面运动机构）和CFD计算为主，经验回归式（如基于船型参数的藤井公式）可作为初值。下面给出一个可直接嵌入仿真脚本的阻尼参数组织示例：

% 阻尼与水动力导数（示例假设值）
Xuu = 0.5 * rho * 0.8 * B*T;   % 纵荡二次阻力系数
Yv  = -0.35 * m;               % 横荡阻尼力对 v 的导数
Yr  = -0.10 * m * L;           % 横荡阻尼力对 r 的导数
Nv  = -0.08 * m * L;           % 艏摇力矩对 v 的导数
Nr  = -0.06 * m * L^2;         % 艏摇力矩对 r 的导数

注意 \(Y_v\)、\(N_r\) 等导数本身为负值，表示阻尼作用方向与运动方向相反；在代码中保留负号，可避免在运动方程中重复取负导致符号错误。

舵力模型与螺旋桨推力模型

舵力通常用机翼升力公式描述：法向力 \(F_N = \frac{1}{2}\rho A_R V_R^2 C_L(\alpha)\)，其中 \(A_R\) 为舵面积，\(V_R\) 为舵处有效来流速度，攻角 \(\alpha\) 由舵角δ与船体横向速度、艏摇引起的局部漂角合成。小舵角下 \(C_L\) 近似线性，大舵角时需考虑失速饱和。舵力对重心的力矩还需乘以力臂 \(x_R\)（舵到重心的纵向距离）。

螺旋桨推力常用四象限模型或简化公式 \(T = \rho n^2 D^4 K_T(J)\)，其中进速比 \(J = u(1-w)/(nD)\)，w为伴流系数，D为桨径，\(K_T\) 由敞水特性曲线插值得到。仿真中若只关注操纵性，可用二次近似 \(T = K_1 n^2 - K_2 n u\) 简化，但需保证在常用航速区间内与敞水曲线吻合。

典型船舶参数估算流程与结果判读

把上述参数组织成MATLAB可用的结构体，是本节最实用的交付物。建议按“船型几何—质量惯性—水动力导数—执行机构”四层组织，每层保留来源标注，便于追溯。下面给出一个可直接运行的参数装配脚本框架：

ship.L = 80; ship.B = 13; ship.T = 4.2; ship.Cb = 0.82;
ship.rho = 1025;
ship.m  = ship.rho*ship.L*ship.B*ship.T*ship.Cb;
ship.mx = 0.08*ship.m;  ship.my = 0.90*ship.m;
ship.Jz = 0.05*ship.m*ship.L^2;
ship.Xuu = 0.5*ship.rho*0.8*ship.B*ship.T;
ship.Yv = -0.35*ship.m; ship.Yr = -0.10*ship.m*ship.L;
ship.Nv = -0.08*ship.m*ship.L; ship.Nr = -0.06*ship.m*ship.L^2;
ship.AR = 6.5; ship.xR = -0.5*ship.L;   % 舵面积与力臂
ship.D  = 2.6; ship.w  = 0.25;          % 桨径与伴流系数

参数装配完成后，应做两项基本校验。第一是量纲与量级校验：附加质量与船体质量同量纲，阻尼导数应使横荡、艏摇运动在无控条件下呈衰减趋势，若仿真中出现横荡速度持续发散，多半是 \(Y_v\) 符号取反或量级过小。第二是稳态校验：给定固定舵角与转速，仿真至稳态后，检查回转直径与航速是否落在该船型的合理区间（内河货船满舵回转直径通常为3～5倍船长），以此判断舵力与推力参数是否匹配。若回转直径明显偏大，优先检查舵面积与 \(C_L\) 斜率；若稳态航速偏低，则检查 \(K_T\) 插值区间与伴流系数。

建议在参数结构体中为每个系数增加 source 字段（如 'empirical'、'CFD'、'PMM'），仿真报告中一并输出，便于评审时区分假设值与辨识值。

最后需要提醒：本节给出的数值均为教学示例假设，不能替代实船资料。实际项目中，附加质量与阻尼导数应优先采用船模试验或CFD结果，经验公式仅用于初值；舵力与推力模型则需用试航数据（如Z形试验、回转试验）反向校核。只有参数来源清晰、符号约定统一、稳态结果可判读，后续的航向保持与轨迹跟踪仿真才具备工程参考价值。