第4章：刚体运动学基础：线速度、角速度、牛顿-欧拉方程在飞行器上的应用

各位同学，大家好。今天我们聊一个非常核心的话题——刚体运动学。说白了，就是研究飞机这个“硬邦邦”的东西，在空中是怎么动起来的。

我刚开始做飞行仿真的时候，总觉得这玩意儿太理论了，不就是个速度加个角速度嘛。直到有一次，我负责的六自由度模型在俯仰机动时，算出来的过载和试飞数据差了将近15%。查了三天，最后发现是角速度的耦合项没加进去。嗯，从那以后，我再也不敢小看这些“基础”了。

4.1 刚体假设与飞行器建模

首先，我们得明确一个前提：刚体假设。

什么叫刚体？就是物体内部任意两点之间的距离，在运动过程中保持不变。说白了，飞机不会像橡皮泥一样变形。

你可能会问：“飞机机翼在飞行中明明会弯曲啊？”没错，真实飞机确实有弹性变形。但在初步的飞行动力学分析中，我们忽略这些弹性效应。为什么？因为弹性模态的频率通常远高于刚体模态，对短周期运动影响不大。我个人习惯在初步设计阶段，先用刚体模型把控制律调通，再考虑弹性修正。

4.2 线速度与角速度：两个“动”的区别

飞行器的运动，说白了就两种：平动和转动。

• 平动：质心在空间中的移动，用线速度 V 描述。
• 转动：飞行器绕质心的旋转，用角速度 ω 描述。

这里有个容易混淆的点：线速度是矢量，有大小有方向。角速度也是矢量，但它的方向是沿着旋转轴的，用右手定则判断。

我记得有一次带新人做项目，他问我：“老师，飞机的滚转角速度 p 是正的，为什么机翼的升力反而减小了？”我告诉他，角速度本身不直接产生力，它产生的是运动学耦合。你想想看，滚转时左右机翼的相对气流速度不一样，这才是升力变化的原因。

4.3 坐标系与运动学关系

要描述飞行器的运动，我们得先选好“尺子”——坐标系。常用的有三个：

坐标系	原点	轴定义	用途
地面坐标系	地面某点	X北、Y东、Z地	导航、轨迹
机体坐标系	飞行器质心	X机头、Y右翼、Z向下	气动力、力矩
气流坐标系	飞行器质心	X沿来流、Z在对称面	升力、阻力

运动学关系，说白了就是不同坐标系之间的转换。比如，地面系下的速度，怎么变成机体系下的速度？

这里有一个经典的公式：

V_body = R * V_earth

其中 R 是旋转矩阵，由欧拉角（滚转 φ、俯仰 θ、偏航 ψ）构成。我建议你把这个矩阵背下来，或者至少知道怎么推导。我在做飞控代码时，这个矩阵几乎每天都要用。

4.4 牛顿-欧拉方程：飞行的“物理引擎”

好了，有了运动学，我们还需要动力学——也就是力和运动的关系。这就是牛顿-欧拉方程。

牛顿第二定律用于平动：

F = m * (dV/dt + ω × V)

欧拉方程用于转动：

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

这里有个坑，我当年就踩过。你看第二个方程，有个 ω × (I * ω) 项，这叫陀螺力矩项。很多人写代码时容易漏掉它。我曾经在仿真一个高速滚转的导弹时，发现滚转通道发散得特别快。查了半天，原来是忘了加这一项。加上之后，模型就稳定了。

4.5 在飞行器上的具体应用

我们把这些方程用到飞行器上，会得到什么？

首先，把力 F 分解为：气动力、推力、重力。力矩 M 分解为：气动力矩、推力力矩。

然后，代入牛顿-欧拉方程，就得到了飞行器的六自由度运动方程。这六个方程，是飞行仿真、飞控设计、稳定性分析的基础。

我个人习惯，在写仿真代码时，会把方程写成状态空间的形式：

状态向量: [u, v, w, p, q, r, φ, θ, ψ, x, y, z]
输入向量: [δe, δa, δr, δT]

其中 u,v,w 是机体系下的线速度分量，p,q,r 是角速度分量。这样写的好处是，可以直接用数值积分方法（比如龙格-库塔法）进行仿真。

4.6 知识体系总览

为了让大家更直观地理解本章的知识结构，我画了一张图：

这张图把本章的知识脉络串起来了。从刚体假设出发，到两种运动形式，再到坐标系和欧拉角，最后落到牛顿-欧拉方程和实际应用。你写代码的时候，可以照着这个结构来组织你的仿真程序。

4.7 一个小例子：滚转时的耦合效应

最后，我给大家留一个思考题。假设飞机在做纯滚转运动（p ≠ 0，q = r = 0），初始时迎角为 5°。请问：滚转过程中，飞机的俯仰角速度 q 会不会发生变化？

答案是：会。因为滚转时，升力矢量在机体坐标系下发生了旋转，产生了俯仰方向的力矩。这就是运动学耦合的一个典型例子。

好了，这一章的内容就到这里。刚体运动学是飞行器动力学的基础，看似简单，但细节很多。希望大家在写代码或做仿真时，能时刻想着这些“基础”的东西，它们往往决定了你的模型准不准。

---

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