3、动力学建模入门：牛顿-欧拉方程、漂浮体的质量与惯量、外力与力矩分析

各位同学，欢迎来到动力学建模的入门课。说实话，很多搞控制的人一听到「建模」两个字就头疼。我当年也一样，总觉得这是搞力学的人的事。但后来我发现，不懂建模，控制算法就是空中楼阁。你想想看，连被控对象长什么样都不知道，你怎么设计控制器？

这一节，我们就来啃下这块硬骨头。别怕，我会用最接地气的方式讲清楚。

3.1 牛顿-欧拉方程：漂浮体的运动圣经

漂浮体在太空或水下运动，本质上就是刚体运动。描述这种运动，最经典的工具就是牛顿-欧拉方程。说白了，牛顿管平动，欧拉管转动。

平动部分（牛顿第二定律）：

m * a = F_ext

其中 m 是质量，a 是加速度，F_ext 是合外力。这个公式大家高中就学过，但在漂浮体上，要注意一点：所有量都要在惯性系下表达。我在做卫星姿态控制时，就吃过这个亏——把体坐标系下的力直接代进去算，结果全错了。

转动部分（欧拉方程）：

I * ω_dot + ω × (I * ω) = τ_ext

这里 I 是惯量张量，ω 是角速度，τ_ext 是外力矩。注意那个叉乘项 ω × (I * ω)，这是漂浮体特有的「陀螺耦合项」。为什么会有它？因为漂浮体在旋转时，不同轴的转动会互相影响。我见过不少新手直接忽略这一项，结果仿真出来的姿态乱飘。

3.2 质量与惯量：漂浮体的「体重」和「身材」

质量 m 很简单，就是物体所含物质的多少。但在漂浮体控制中，我们更关心的是惯量。惯量说白了就是物体抵抗转动的能力。质量越大越难推，惯量越大越难转。

惯量张量 I 是一个 3×3 的矩阵：

I = | Ixx  Ixy  Ixz |
    | Iyx  Iyy  Iyz |
    | Izx  Izy  Izz |

对角线元素 Ixx、Iyy、Izz 是主惯量，非对角线元素是惯性积。如果物体质量分布对称，惯性积为零。我建议你们在实际项目中，尽量把传感器和执行器对称安装，这样惯量矩阵就是对角阵，控制会简单很多。

惯量的计算：

• 简单形状：球体 I = (2/5)mR²，圆柱 I = (1/2)mR²（绕轴线）
• 复杂形状：用 CAD 软件或实验测量
• 变惯量：燃料消耗、机械臂运动都会改变惯量

3.3 外力与力矩分析：谁在推你？

漂浮体受到的外力和力矩，主要来自这几个方面：

类型	来源	特点
重力	地球引力	均匀场，作用在质心
推力	推进器/发动机	可控，但有限幅
扰动力	环境干扰	不可控，需要补偿
科里奥利力	旋转参考系	速度相关，容易忽略

力矩分析要特别注意：

• 推力不通过质心时，会产生力矩
• 重力在质心处不产生力矩，但重心偏移时会产生恢复力矩
• 流体阻力会产生与速度平方成正比的阻尼力矩

3.4 知识体系总览

下面这张图，是我自己总结的动力学建模核心逻辑。你看一遍，就能把整个章节串起来。

这张图把整个章节串起来了。你从牛顿-欧拉方程出发，左边走平动路线，右边走转动路线。平动需要质量，转动需要惯量。然后两边分别加上外力和外力矩，最后汇合到动力学模型，用于控制算法设计。嗯，就是这么个逻辑。

3.5 小结与实操建议

这一节的内容，说白了就是三件事：

1. 记住牛顿-欧拉方程的形式，尤其是那个叉乘项
2. 搞清楚质量和惯量，惯量张量一定要测量或计算准确
3. 列出所有外力和力矩，别漏掉科里奥利力和扰动力

我建议你找个简单的漂浮体模型（比如一个立方体），手动推导一遍它的动力学方程。写出来，然后仿真看看。我第一次做的时候，花了整整一个周末，但做完之后，对建模的理解就完全不一样了。