4. 刚体动力学方程：牛顿-欧拉方程、平动动力学、转动动力学、惯性张量与质量特性

各位同学，今天我们来啃一块硬骨头——刚体动力学方程。说实话，这块内容我当年学的时候也觉得头大，但后来在项目中摔过几次跟头，才真正体会到它的分量。

刚体动力学，说白了就是研究飞行器怎么动、为什么这么动。你想想看，一架飞机在天上飞，既有整体的平移，又有绕重心的旋转，这两者还互相耦合。牛顿-欧拉方程就是描述这个过程的数学工具。

4.1 牛顿-欧拉方程：飞行的基本法则

牛顿-欧拉方程把刚体运动拆成两部分：平动和转动。平动用牛顿第二定律，转动用欧拉方程。我个人习惯把它们写在一起：

平动：F = m * a_cm
转动：M = I * α + ω × (I * ω)

这里F是合外力，m是质量，a_cm是质心加速度。M是合力矩，I是惯性张量，α是角加速度，ω是角速度。

嗯，这里要注意：转动方程里那个叉乘项ω × (I * ω)，很多人容易漏掉。我在项目中遇到过一位同事，仿真结果总是不对，查了两天才发现是忘了加这一项。这个叉乘项代表的是陀螺力矩，在高速旋转时影响非常大。

4.2 平动动力学：飞行器怎么走

平动动力学研究的是飞行器质心的运动轨迹。说白了，就是飞机往哪个方向飞、飞多快。

在惯性坐标系下，平动方程很简单：

m * dV/dt = F_gravity + F_aero + F_thrust

其中V是速度矢量，F_gravity是重力，F_aero是气动力，F_thrust是推力。

但实际仿真时，我们通常把方程投影到机体坐标系下。为什么？因为气动力和推力在机体坐标系下表达更直接。我曾经做过一个项目，把方程全放在惯性系里算，结果每次都要做坐标变换，代码又乱又容易出错。后来改成机体坐标系，清爽多了。

机体坐标系下的平动方程长这样：

m * (dV_b/dt + ω_b × V_b) = F_b

多出来的ω_b × V_b项，就是所谓的「牵连加速度」。你想想看，飞机在旋转，机体坐标系也在转，所以速度的变化率要加上这个修正项。

4.3 转动动力学：飞行器怎么转

转动动力学描述的是飞行器绕质心的旋转运动。这部分比平动复杂得多，因为涉及到惯性张量和角速度的耦合。

欧拉方程在机体坐标系下的形式：

I * dω/dt + ω × (I * ω) = M

这里I是惯性张量，是一个3×3的矩阵。ω是角速度矢量，M是外力矩。

为什么会有ω × (I * ω)这一项？我打个比方：你坐在转椅上，手里拿着两个哑铃。当你把哑铃收回来时，转速会变快——这就是角动量守恒。ω × (I * ω)项描述的就是这种角动量在旋转坐标系下的变化。

实际仿真时，我们通常把欧拉方程展开成三个标量方程：

Ixx * dp/dt = (Iyy - Izz) * q * r + Mx
Iyy * dq/dt = (Izz - Ixx) * r * p + My
Izz * dr/dt = (Ixx - Iyy) * p * q + Mz

其中p、q、r分别是绕x、y、z轴的角速度分量。Ixx、Iyy、Izz是主惯性矩。

4.4 惯性张量与质量特性：飞行器的「体质」

惯性张量，说白了就是描述飞行器质量分布的一个矩阵。它决定了飞行器在受到力矩时，会怎么转、转多快。

惯性张量的完整形式：

I = [Ixx  -Ixy  -Ixz
     -Ixy  Iyy  -Iyz
     -Ixz  -Iyz  Izz]

对角线元素Ixx、Iyy、Izz是主惯性矩，非对角线元素Ixy、Ixz、Iyz是惯性积。如果飞行器关于坐标平面对称，惯性积就为零。

质量特性包括：

• 总质量m：决定平动惯性
• 质心位置：气动力矩的力臂
• 惯性张量I：决定转动惯性

这些参数怎么来的？通常来自CAD模型的质量属性分析，或者实验测量。我在项目中常用SolidWorks或CATIA导出这些数据，但要注意：CAD模型和实际产品总有偏差，最好留点裕量。

4.5 知识体系总览

下面这张图是我自己整理的刚体动力学知识框架，方便大家理解各部分之间的关系：

这张图把刚体动力学拆成了三个模块。平动和转动通过角速度ω耦合在一起，而惯性张量则是转动部分的「灵魂参数」。我建议大家在写仿真代码时，也按照这个结构来组织代码，模块清晰，调试方便。

好了，刚体动力学方程就讲到这里。记住：牛顿-欧拉方程是基础，平动和转动要分开理解，惯性张量是转动特性的关键。下次我们聊姿态表示方法——四元数和欧拉角，那又是另一个有意思的话题了。

---

公众号:蓝海资料掘金营，微信deep3321