2. 过渡段动力学建模:多旋翼与固定翼气动耦合
好,咱们进入正题。过渡段动力学建模,说白了就是搞清楚飞机在从多旋翼模式切换到固定翼模式时,到底受了哪些力、这些力怎么互相影响。我刚开始做VTOL时,总觉得过渡段就是“旋翼慢慢停、机翼慢慢飞”,结果第一次试飞差点栽跟头——飞机在中间状态突然抬头,差点翻过去。嗯,从那以后我就老老实实把动力学模型先吃透了。
2.1 多旋翼与固定翼气动耦合
先问一个问题:过渡段里,旋翼和机翼的气流是各干各的吗?
当然不是。旋翼产生的下洗气流会直接打在机翼上,改变机翼的有效迎角。反过来,前飞速度带来的来流也会改变旋翼的拉力方向。这就是气动耦合——两个系统互相“捣乱”。
我个人习惯把耦合分成两类:
- 旋翼对机翼的干扰:旋翼下洗气流增加机翼局部风速,导致升力分布不均。我在项目中遇到过,机翼内侧先失速,就是因为下洗太强。
- 机翼对旋翼的干扰:前飞时,机翼上表面低压区会“吸”旋翼尾流,改变旋翼入流分布。说白了,旋翼的拉力效率会下降。
核心结论:过渡段中,气动耦合不可忽略。忽略它,你的控制律大概率会发散。
2.2 过渡段受力分析
咱们把飞机当成一个刚体,看看它受哪些力。你想想看,这时候飞机既不是纯旋翼机,也不是纯固定翼,所以受力很杂。
| 力的类型 | 来源 | 特点 |
|---|---|---|
| 重力 | 地球引力 | 始终向下,大小不变 |
| 旋翼拉力 | 多旋翼电机+桨 | 方向可调,随转速变化 |
| 气动升力 | 机翼+机身 | 随空速和迎角变化 |
| 气动阻力 | 全机 | 与空速平方成正比 |
| 耦合干扰力 | 旋翼-机翼相互作用 | 非线性,难建模 |
我记得有一次做风洞实验,发现耦合干扰力在过渡段中间区域(空速15-20m/s)最大,能占到总升力的15%。这个量级,控制上必须考虑。
2.3 关键状态变量定义
建模之前,先定好状态变量。我建议用下面这组,工程上最常用:
- 空速 V:飞机相对空气的速度,单位 m/s。过渡段核心变量。
- 迎角 α:机翼弦线与来流夹角。注意,旋翼下洗会改变局部迎角。
- 俯仰角 θ:机体轴与水平面夹角。控制律里常用。
- 旋翼总拉力 T:所有旋翼拉力的矢量和。
- 过渡进度因子 σ:0表示纯多旋翼模式,1表示纯固定翼模式。这是我个人习惯引入的变量,方便做增益调度。
小技巧:σ 可以定义为空速的函数,比如 σ = V / V_transition,其中 V_transition 是过渡完成速度。这样控制律里直接查表就行。
2.4 简化模型推导
好,现在咱们推导一个工程可用的简化模型。为什么要简化?因为完整的气动耦合模型是偏微分方程组,飞控芯片算不过来。
假设条件:
- 忽略旋翼对机翼的诱导阻力(工程上影响较小)
- 假设旋翼拉力方向始终垂直于机体(近似成立)
- 忽略侧向力和横航向耦合(先做纵向)
那么,纵向动力学方程可以写成:
m * dV/dt = T * cos(α+θ) - D - m*g*sin(θ)
m * V * dγ/dt = T * sin(α+θ) + L - m*g*cos(θ)
I_y * d²θ/dt² = M_aero + M_rotor
其中:
- D 为阻力,L 为升力,M_aero 为气动俯仰力矩
- M_rotor 为旋翼产生的俯仰力矩
- γ 为航迹角
这里要注意,升力 L 不仅来自机翼,还包含旋翼下洗的贡献。我一般用修正系数 k_couple 来近似:
L = 0.5 * ρ * V² * S * CL(α) + k_couple * T
k_couple 的值可以通过 CFD 或风洞标定得到。我在项目中取 0.1~0.3,具体看布局。
避坑指南:我曾经在简化模型里直接忽略了 M_rotor,结果仿真时飞机在过渡段末端出现高频振荡。后来发现旋翼转速变化产生的陀螺力矩在低速时不可忽略。所以,M_rotor 项至少保留一个一阶近似。
2.5 知识体系总览
下面这张图是我自己画的,把本章的核心逻辑串起来了。你一看就明白:
这张图从左到右、从上到下,就是咱们建模的完整流程。先搞清楚三个输入,然后做受力分析,再定义状态变量,最后推导出简化模型。每一步都有坑,但按这个框架走,至少不会漏东西。
公众号:蓝海资料掘金营,微信deep3321