第4章刚体动力学基础

4. 刚体动力学基础：牛顿-欧拉方程、动量与动量矩定理、转动惯量

各位同学，今天我们来聊聊刚体动力学。说实话，这是飞行器建模里最绕不开的一块硬骨头。我当年刚入行时，总觉得这部分理论太抽象，直到第一次做六自由度仿真，模型在天上乱飞，我才意识到——嗯，基础不牢，地动山摇。

刚体动力学，说白了就是研究飞行器这个“硬家伙”在力和力矩作用下怎么动。我们不考虑它变形，只关心它作为一个整体怎么飞、怎么转。今天我会把牛顿-欧拉方程、动量与动量矩定理、转动惯量这三个核心概念串起来讲，结合我这些年踩过的坑，帮你把这块地基打牢。

我个人习惯把牛顿-欧拉方程看作飞行器动力学的“宪法”。它分两部分：

写成数学形式就是：

牛顿方程：F = m * a_c
欧拉方程：M = I * α + ω × (I * ω)

这里有个细节我特别想强调——坐标系的选择。我在项目中遇到过好几次，新手直接把牛顿方程写在机体坐标系里，结果加速度项算出来全是错的。为什么？因为牛顿方程只在惯性系下成立。你想想看，飞行器自己在转，你在它身上测加速度，那能准吗？

核心要点：牛顿方程必须在惯性坐标系（如地心惯性系）下写，欧拉方程通常在机体坐标系下写更方便。这是刚体动力学建模的第一条铁律。

实际工程中，我们通常的做法是：

这样做的好处是转动惯量矩阵在机体坐标系下是常数，省去了实时更新的麻烦。我曾经见过一个团队把转动惯量写成时变矩阵，仿真步长稍微大一点就发散——其实完全没必要。

有时候直接用牛顿-欧拉方程求解挺麻烦的，尤其是系统里有外力不好建模的时候。这时候，动量定理和动量矩定理就派上用场了。

动量定理：

dP/dt = F_ext
其中 P = m * v_c（动量 = 质量 × 质心速度）

动量矩定理：

dH/dt = M_ext
其中 H = I * ω（动量矩 = 转动惯量 × 角速度）

说白了，这两个定理就是牛顿-欧拉方程的积分形式。但为什么还要单独拿出来讲？因为有些场景下用它们更方便。

举个例子。我在做某型无人机碰撞分析时，撞击力很难精确测量。但如果用动量定理，只需要知道碰撞前后的速度变化，就能反推出冲量大小。这就是工程上的“避重就轻”——算不了力，就算动量变化。

实用技巧：当系统受到冲击载荷或短时强作用力时，优先考虑用动量/动量矩定理。我曾经用这个方法在半天内解决了某型号的着陆冲击分析，而用传统方法至少要算一周。

这里要特别提醒一点：动量矩定理中的矩点选择很关键。如果选质心，形式最简单；如果选其他点，会多出惯性力项。我建议初学者一律选质心作为矩点，等熟练了再考虑其他选择。

转动惯量，说白了就是衡量飞行器抵抗转动的能力。质量越大、质量分布离转轴越远，转动惯量就越大。

对于飞行器来说，我们通常关心三个轴的转动惯量：

除此之外，还有惯性积 I_xy、I_xz、I_yz。它们描述的是质量分布不对称带来的耦合效应。如果飞行器完全对称，惯性积为零。但现实中的飞行器多少都有点不对称——比如单发螺旋桨飞机的螺旋桨旋转就会产生不对称力矩。

避坑指南：我曾经在建模时忽略了惯性积，结果仿真出来的横航向响应跟试飞数据对不上。查了三天才发现是惯性积没加。对于非对称布局的飞行器（如直升机、倾转旋翼机），惯性积绝对不能忽略。

转动惯量的计算有两种方法：

我个人建议：理论计算和实验测量都要做。理论计算帮你理解质量分布的影响，实验测量给你真实数据。两者对得上，说明模型靠谱；对不上，说明有地方漏了。

下面这张图是我自己总结的刚体动力学知识框架，你可以把它当作学习路线图：

从这张图可以看得很清楚：牛顿-欧拉方程是顶层总纲，往下分支出平动和转动两套理论，再往下是具体的定理和参数，最后落到工程应用上。我建议你学习时也按这个顺序来，先搭框架，再填细节。

今天的内容就到这里。总结三个你必须记住的点：

最后说句掏心窝的话：刚体动力学这部分，光看书是学不会的。我建议你找一个小型飞行器的参数，自己动手建个六自由度模型，跑一遍仿真。遇到问题再回来翻书，印象会深得多。下一章我们会讲运动学，把位置、姿态、角速度这些量之间的关系理清楚——嗯，到时候见。