2. 飞行器动力学基础:刚体运动方程、六自由度模型、质量变化对动力学的影响

各位同学,欢迎来到《自适应控制》的第二讲。今天咱们聊聊飞行器动力学的基础。说实话,这部分内容我当年学的时候也觉得有点枯燥,一堆公式和坐标系。但后来在实际项目中吃过亏,才明白这些基础有多重要。

你想想看,一架飞机在空中飞,它到底是怎么运动的?我们怎么用数学去描述它?更关键的是,当它投下货物、消耗燃油,质量一直在变,我们的控制律还能不能稳住?嗯,这就是今天要解决的核心问题。

2.1 刚体运动方程:从牛顿到欧拉

首先,我们把飞行器简化成一个刚体。什么叫刚体?就是忽略弹性形变,认为飞机内部各个点的相对位置不变。当然,真实飞机肯定有弹性,但那是高阶话题了,咱们先把刚体搞明白。

刚体运动可以拆成两部分:质心的平动绕质心的转动。平动用牛顿第二定律,转动用欧拉方程。我个人习惯把这两个方程并排写在笔记本上,每次建模前先看一眼,提醒自己别漏项。

平动方程(牛顿第二定律):

m * dV/dt = F_ext

其中 m 是质量,V 是速度矢量,F_ext 是合外力(重力、推力、气动力)。

转动方程(欧拉方程):

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

其中 I 是惯性张量,ω 是角速度矢量,M_ext 是合外力矩。

这里有个坑,我刚开始做仿真时经常掉进去——坐标系的选择。平动方程通常在地面惯性系下写,转动方程在机体坐标系下写。两个坐标系之间的转换,稍不注意就会出错。

我的小技巧:写代码前,先把坐标系画在纸上,标清楚每个矢量的分量是在哪个系下表达的。别嫌麻烦,这一步能省你后面三天调试时间。

2.2 六自由度模型:三个平动 + 三个转动

六自由度(6-DOF)模型,说白了就是描述飞行器在空间中的全部运动状态。哪六个?三个位置自由度(x, y, z)和三个姿态自由度(滚转 φ、俯仰 θ、偏航 ψ)。

为什么是六个?因为刚体在三维空间里,能做的运动无非就是前后左右上下移动,加上绕三个轴的旋转。你想想看,一架飞机不可能像蛇一样扭来扭去,对吧?

完整的六自由度模型通常写成状态空间形式。我个人习惯把状态向量写成这样:

状态向量 X = [x, y, z, φ, θ, ψ, u, v, w, p, q, r]^T

其中:

  • (x, y, z) —— 位置(地面坐标系)
  • (φ, θ, ψ) —— 姿态角(欧拉角)
  • (u, v, w) —— 速度在机体轴的分量
  • (p, q, r) —— 角速度在机体轴的分量

这里要注意,欧拉角在俯仰角 θ = ±90° 时会出现万向锁问题。我在做某型无人机仿真时就遇到过,俯仰拉到90度,偏航和滚转突然耦合在一起,控制律直接炸了。后来改用四元数才解决。所以,如果你做全姿态飞行(比如特技、导弹),建议用四元数代替欧拉角。

2.3 质量变化对动力学的影响

好,重点来了。前面讲的都是质量不变的情况。但实际飞行中,质量怎么可能不变?燃油消耗、货物投放、乘客走动……这些都是质量变化源。

质量变化会从两个层面影响动力学:

  1. 直接影响平动方程:m 在变,所以 d(mV)/dt 不能简单写成 m * dV/dt,而要考虑 dm/dt 项。
  2. 影响转动惯量和质心位置:质量分布变了,惯性张量 I 和质心位置都会变,进而影响欧拉方程。

咱们先看平动方程。考虑质量变化后,正确的形式是:

d(mV)/dt = F_ext
→ m * dV/dt + V * dm/dt = F_ext

看到了吗?多了一项 V * dm/dt。这一项在质量快速变化时(比如火箭级间分离、空投物资)绝对不能忽略。我曾经在做一个空投仿真时,忽略了这一项,结果飞机在投放瞬间出现了一个诡异的加速度跳变,查了两天才找到原因。

再看转动部分。质量变化会导致惯性张量 I 随时间变化。欧拉方程变成:

d(I * ω)/dt + ω × (I * ω) = M_ext
→ I * dω/dt + dI/dt * ω + ω × (I * ω) = M_ext

多出来的 dI/dt * ω 项,就是质量分布变化带来的附加力矩。比如燃油从机翼油箱向机身油箱转移,飞机的滚转惯量会减小,这时候如果不补偿,滚转控制就会变得过于灵敏。

注意:质量变化不仅影响动力学方程,还会影响控制系统的稳定性。传统的 PID 控制器在质量变化超过 20% 时,性能往往会显著下降。这就是为什么我们需要自适应控制——它能在线调整控制器参数,适应这些变化。

2.4 知识体系总览

为了让大家对本章内容有个整体把握,我画了一张图。这张图把刚体运动、六自由度模型、质量变化影响串在了一起,你看完应该能明白它们之间的逻辑关系。

飞行器动力学基础:知识体系 刚体运动方程 平动方程(牛顿第二定律) 转动方程(欧拉方程) 六自由度模型 (6-DOF) 3个平动自由度 (x,y,z) 3个转动自由度 (φ,θ,ψ) 质量变化扰动 影响 m、I、质心位置

从这张图可以看得很清楚:刚体运动方程是基础,它衍生出平动和转动两个分支,两者合在一起构成了六自由度模型。而质量变化扰动,就像一只看不见的手,同时影响着平动方程中的质量 m 和转动方程中的惯性张量 I。

2.5 实际项目中的体会

最后,我想分享一点个人经验。我在做某型运输机的自适应控制器时,发现质量变化带来的最大问题不是参数变化本身,而是你很难准确测量当前的质量和质心位置。燃油消耗还能通过流量计估算,但货物装载情况呢?乘客分布呢?这些信息往往是不精确的。

所以,在实际工程中,我们通常不会直接测量质量,而是通过在线辨识的方法去估计。比如用递推最小二乘法(RLS)实时估计质量参数,然后反馈给自适应控制器。这部分内容,我们会在后面的章节详细展开。

好了,今天的内容就到这里。记住一句话:飞行器动力学是自适应控制的地基,地基不稳,楼盖得再高也得塌。下一章,我们会进入自适应控制的核心——模型参考自适应控制(MRAC),到时候你会看到,今天学的动力学方程会直接出现在参考模型里。

本章要点回顾:

  • 刚体运动 = 质心平动 + 绕质心转动
  • 六自由度模型包含 3 个位置 + 3 个姿态
  • 质量变化引入 dm/dt 和 dI/dt 附加项
  • 实际工程中质量参数常通过在线辨识获取

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