3. AUV动力学建模的力学基础

要为AUV建立六自由度动力学模型,第一步不是急着写矩阵,而是先明确“用哪套力学语言描述运动”。工程上最常用的两条路径是牛顿-欧拉方程和拉格朗日方程。二者描述的是同一个物理世界,只是切入角度不同:前者从力与力矩的平衡出发,逐体分析;后者从能量出发,用广义坐标统一描述。理解它们的差异,才能判断在AUV建模中何时用哪一种、如何互相校验。

一、牛顿-欧拉方程:从力与力矩直接出发

牛顿第二定律给出平动关系,欧拉方程给出转动关系。对刚体而言,在惯性系下可写为:合力等于质量乘以质心加速度,合力矩等于角动量变化率。真正用于AUV建模时,通常把方程投影到艇体坐标系(原点取在艇体某固定点,常见为重心或浮心,x轴指向艇首,z轴指向艇底,y轴按右手系确定)。设艇体速度为 v=[u, v, w]T,角速度为 ω=[p, q, r]T,则刚体六自由度方程可整理成统一形式:

M_RB · ν̇ + C_RB(ν) · ν = τ_RB
ν = [u, v, w, p, q, r]^T

其中 M_RB 是刚体广义质量矩阵,包含质量m、惯性张量I及其耦合项;C_RB(ν) 是刚体科氏-向心项矩阵,由速度的二次项构成;τ_RB 是作用在艇体上的外力与力矩(螺旋桨推力、舵力、水动力、浮力与重力等)。牛顿-欧拉法的优点是物理图像清晰:每一项力都能对应到具体部件和具体作用点,便于分项建模和逐项排查。缺点是当坐标系选择变化或需要引入约束时,力矩项和耦合项容易写错符号。

二、拉格朗日方程:从能量统一描述

拉格朗日方程以广义坐标 q 为变量,通过动能T与势能V构造拉格朗日函数 L=T−V,方程为:

d/dt(∂L/∂q̇) − ∂L/∂q = Q

对AUV,常取位置与姿态角作为广义坐标,动能由平动动能与转动动能组成,势能主要来自重力与浮力的高度差。拉格朗日法的优势在于:系统能量是标量,坐标变换时不易出错;引入附加质量、附加惯性等水动力效应时,可以把它们并入动能项统一处理;对多体或带约束的系统扩展性好。缺点是中间推导量大,且外力(如螺旋桨推力、舵力)仍需通过广义力Q逐项映射,物理直观性不如牛顿-欧拉法。

三、两种方法在AUV建模中的适用性对比

实际建模中,二者并非互斥。常见做法是:用牛顿-欧拉法搭建刚体主体方程,因为AUV的推力、舵力、浮力作用点明确,便于分项;用拉格朗日法处理附加质量与附加惯性,或在验证阶段用能量法交叉校验。下面这张图对比了两条建模路径的输入、处理与输出关系,并给出汇合到统一六自由度形式的衔接点。

AUV 六自由度建模的两条力学路径 牛顿-欧拉法(力/力矩路径) 输入:艇体速度 ν、角速度 ω 外力:推力、舵力、水动力、浮力/重力 处理:逐体受力分析 投影到艇体坐标系,写力矩平衡 输出:M_RB、C_RB(ν)、τ_RB 物理直观,便于分项排查 注意:符号与坐标系约定 力矩耦合项易出错,需统一参照系 拉格朗日法(能量路径) 输入:广义坐标 q、q̇ 位置、姿态角及其导数 处理:构造 T、V 与 L=T−V 附加质量并入动能,势能含浮力高度差 输出:广义力 Q 与运动方程 标量能量推导,坐标变换稳健 注意:外力映射为广义力 推力/舵力需逐项投影,直观性较弱 统一六自由度一般形式 M·ν̇ + C(ν)·ν + D(ν)·ν + g(η) = τ 刚体主体方程 附加质量/校验 M = M_RB + M_A,C = C_RB + C_A

四、建立刚体六自由度一般形式

把刚体项与水动力项合并,AUV六自由度动力学一般形式通常写成:

M · ν̇ + C(ν) · ν + D(ν) · ν + g(η) = τ
M = M_RB + M_A
C(ν) = C_RB(ν) + C_A(ν)

这里 M 是广义质量矩阵,由刚体质量矩阵 M_RB 与附加质量矩阵 M_A 叠加;C(ν) 是科氏-向心矩阵,同样由刚体项与附加质量项组成;D(ν) 是阻尼项,g(η) 是恢复力项(重力与浮力),τ 是控制力与外力。需要强调:C(ν) 必须满足反对称性质(在无阻尼理想情况下),这是能量守恒的体现,也是检查矩阵是否写错的重要手段。

五、广义质量矩阵与科氏力项的具体构成

以艇体坐标系原点取在重心为例,刚体广义质量矩阵可写为分块形式:

M_RB = [ m·I3      −m·S(r_g) ]
       [ m·S(r_g)    I_b     ]

其中 m 为质量,I3 为三阶单位阵,S(·) 为反对称矩阵算子,r_g 为重心相对原点的位置,I_b 为惯性张量。若原点取在重心,则 r_g=0,耦合项消失,矩阵退化为对角块形式,这也是很多教材把原点设在重心的原因。附加质量矩阵 M_A 通常按六个自由度的附加质量与附加惯性系数填写,工程上常假设其对称,非对角项根据艇体对称性取舍。

科氏力项由 M 与速度构造。对分块矩阵,C(ν) 的一般构造规则是:对任意向量 a,定义 S(a) 为反对称矩阵,则 C(ν) 可由 M 的分块与 ν 组合得到,保证 νTC(ν)ν=0。实际编程时,建议不要手写展开式,而是用“分块矩阵+反对称算子”的方式生成,再用数值检查验证反对称性。例如取一组示例参数(仅作教学假设):m=50 kg,I_b=diag(4.2, 6.8, 7.5) kg·m²,r_g=0,附加质量取 M_A=diag(12, 18, 20, 1.5, 2.0, 2.2),代入后检查 C(ν) 是否满足 νTC(ν)ν≈0,若残差超过 1e-9 量级,说明构造有误。

建模小技巧:先用牛顿-欧拉法写出刚体部分,确认符号与坐标系一致;再用拉格朗日法独立推导一遍能量项,两者结果在无外力时应一致。若不一致,优先检查科氏项是否满足反对称性,以及附加质量是否被重复计入。

六、方法选择建议

对于以控制设计为目标的AUV建模,推荐“牛顿-欧拉为主、拉格朗日为辅”:主体方程用牛顿-欧拉法搭建,便于把推力、舵力、浮力逐项对应到物理部件;附加质量与恢复力用能量法处理,减少坐标变换错误;最后用反对称性和能量守恒做数值校验。这样既保留了工程直观性,又利用了拉格朗日法的稳健性,是AUV六自由度建模中较为实用的组合策略。