3. ROV机械臂的静力学与动力学特性
ROV机械臂与陆用工业机械臂最本质的差别,不在于关节形式,而在于它始终工作在“水”这个介质里。水既提供了浮力,也带来了附加质量、阻尼和拖曳力;同时,机械臂各连杆的重力被浮力部分抵消,但抵消程度随姿态、排水体积和海水密度变化。因此,分析ROV机械臂的关节力矩,不能简单套用空气中的重力项,而必须把重力、浮力、水动力载荷分开建模,再叠加到关节力矩方程中。本节围绕这三个物理量展开,说明它们如何进入静力学与动力学模型,以及在实际操作中如何判断它们的影响。
一、坐标系与符号约定
为避免后续公式与图示产生歧义,先统一约定:以机械臂基座为原点建立基坐标系,z 轴竖直向上为正,重力加速度方向为 −z;关节力矩以驱动关节转动方向为正,即关节角增大方向为正。连杆 i 的质量记为 mi,其质心在基坐标系中的位置记为 ri,连杆排水体积记为 Vi。浮力方向恒为 +z,重力方向恒为 −z。全文与图示均采用这一约定,若切换到机械臂本体系,需明确说明变换关系。
二、水下重力补偿:等效重力与“水中重量”
在空气中,连杆 i 对关节产生的重力矩为 mi·g 与力臂的叉乘。进入水中后,连杆受到向上的浮力 Fb,i = ρ·g·Vi,其中 ρ 为海水密度(工程估算常取 1025 kg/m³)。于是连杆的“水中重量”为:
W'_i = m_i·g − ρ·g·V_i = g·(m_i − ρ·V_i)
定义浮力比 ki = ρ·Vi / mi,则 W'_i = mi·g·(1 − ki)。当 ki 接近 1 时,连杆接近中性浮力,重力矩被大幅抵消;当 ki 小于 1 时,连杆仍有“下沉”趋势。实际 ROV 机械臂的连杆多采用铝合金或钛合金加浮力材料,ki 通常在 0.6~0.95 之间,具体取决于是否填充合成泡沫。
需要强调的是,重力补偿不是把重力项删掉,而是用等效重力 W'_i 替换 mi·g。关节静力矩可写为:
τ_gravity,i = Σ (r_j × W'_j·ẑ) 的 z 分量(对转动关节)
其中求和覆盖关节 i 之后的所有连杆。由于浮力作用点通常取连杆排水体积的浮心,而重力作用点在质心,两者不重合时还会产生一个附加力偶,工程上称为“浮心力矩”。当浮心与质心偏差较大(例如连杆一端加装大体积浮力块)时,这一项不可忽略,需在静力学模型中单独计入。
三、水动力载荷:拖曳力、附加质量与阻尼
机械臂在水中运动时,连杆受到的水动力载荷主要包括三类:
- 拖曳力(阻力):Fd = 0.5·ρ·Cd·A·v·|v|,与速度平方相关,方向与运动方向相反。Cd 取决于连杆截面形状,圆柱形连杆典型值约 0.8~1.2,A 为迎流投影面积。
- 附加质量力:Fa = ρ·V·Ca·a,与加速度成正比,相当于给连杆“增加”了一部分质量。Ca 对细长圆柱约 0.9~1.0,对球形约 0.5。附加质量会显著改变动力学方程中的惯性矩阵。
- 粘性阻尼与涡激效应:低速时以线性阻尼为主,高速或特定流速下可能出现涡激振动,表现为周期性侧向力,对精细抓取作业不利。
把这三类载荷折算到关节上,得到关节的水动力力矩 τhydro。它与关节速度、加速度以及连杆姿态都相关,因此动力学方程可写为:
M(q)·q̈ + C(q,q̇)·q̇ + τ_gravity(q) + τ_hydro(q,q̇,q̈) = τ_joint
其中 M(q) 已包含附加质量贡献,τhydro 包含拖曳与阻尼项。与空气中相比,水中 M(q) 更大、C(q,q̇) 更复杂,且 τhydro 与速度强相关,导致机械臂在水中“更沉、更黏、更难急停”。
四、三类载荷对关节力矩的影响对比
为便于理解三者量级与作用方式,下表给出一个教学示例(参数为合理假设,非实测):假设某两连杆 ROV 机械臂,连杆 1 质量 6 kg、排水体积 4.5 L,连杆 2 质量 3 kg、排水体积 2.2 L,海水密度 1025 kg/m³,g 取 9.81 m/s²。
| 载荷类型 | 作用方向 | 与姿态/运动的关系 | 对关节力矩的典型影响 |
|---|---|---|---|
| 重力(空气中) | −z | 仅与姿态有关 | 基准项,随关节角周期变化 |
| 浮力 | +z | 与排水体积、姿态有关 | 抵消约 60%~95% 重力矩 |
| 拖曳力 | 与速度反向 | 与速度平方相关 | 高速运动时显著增大,低速可忽略 |
| 附加质量 | 与加速度反向 | 与加速度相关 | 增大惯性力矩,影响动态响应 |
以连杆 1 为例,其空气中重量约 58.9 N,浮力约 45.2 N,水中重量约 13.7 N,浮力比约 0.77。若连杆 1 质心到关节 1 的水平力臂为 0.25 m,则空气中重力矩约 14.7 N·m,水中仅约 3.4 N·m。可见浮力补偿使静力矩下降约 77%,这正是 ROV 机械臂能采用较小驱动器的原因之一。但若机械臂末端加装工具或样品,附加质量与拖曳力会迅速上升,此时动力学项可能超过静力矩项。
五、实际分析流程与判断依据
在工程实践中,可按以下顺序处理:第一步,测量或估算各连杆质量、质心位置、排水体积与浮心位置,计算浮力比 ki,得到等效重力 W′i;第二步,用等效重力替换重力项,建立静力学模型,计算典型姿态下的关节静力矩,判断驱动器是否满足要求;第三步,根据作业速度范围估算拖曳力与附加质量,判断动力学项是否超过静力矩项;第四步,在控制器中把 τgravity 作为前馈补偿项,把 τhydro 作为阻尼补偿或速度前馈项。
判断依据可以这样设定:若某关节在最大力臂姿态下的静力矩已接近驱动器额定力矩的 80%,说明浮力补偿不足或连杆配置偏重,应增加浮力材料或调整质心;若在快速运动时关节出现明显滞后或过冲,而静力矩余量充足,则说明拖曳力与附加质量占主导,应降低运动速度或提高控制器阻尼增益。对于精细抓取作业,还应关注低速下的粘性阻尼与涡激效应,必要时通过降低末端速度、缩短力臂来减小水动力扰动。
综上,ROV机械臂的关节力矩分析必须把重力、浮力、拖曳力、附加质量分开处理,再按姿态与运动状态叠加。只有明确各载荷的作用方向、量级与变化规律,才能在操作中合理设定速度、选择姿态、配置浮力材料,并对关节力矩做出准确判断。