4、发动机与推力模型：推力特性、油门响应、螺旋桨/喷气发动机简化模型、推力对飞行动力学的影响
发动机与推力模型是飞行控制系统建模中连接能量输入与运动输出的关键环节。在飞控仿真与调参工作中，推力模型不仅决定了飞机的配平点、爬升率和最大平飞速度，还直接影响纵向短周期与长周期模态的阻尼特性。本章将围绕推力特性的数学描述、油门通道的动态响应、螺旋桨与喷气发动机的简化建模方法，以及推力变化对飞行动力学的影响展开，帮助读者建立从油门指令到机体加速度的完整链路认知。
上图展示了发动机推力模型在飞控系统中的位置与内部结构。从油门指令到推力输出，需要经过伺服动态、推力产生和飞行动力学三个环节，同时速度与高度反馈会修正推力计算结果。下面将逐一展开各环节的建模方法与调参要点。
4.1 推力特性的数学描述
推力模型的输入通常是油门杆位置（或油门指令百分比），输出是沿机体纵轴方向的推力。在飞行控制系统中，推力特性需要描述为多个变量的函数：
T = f(δ_th, V, ρ, h, M)其中，δ_th为油门指令（0~1），V为空速，ρ为空气密度，h为高度，M为马赫数。对于亚音速飞行器，马赫数影响可以忽略，但密度和速度的影响必须保留。
推力特性的核心参数包括：
- 最大推力T_max：在海平面、标准大气、静止状态下的最大推力值，是推力模型的标称基准。
- 推力-油门线性度：理想情况下推力与油门呈线性关系，但实际发动机存在非线性，尤其在低油门区域推力响应迟缓。
- 速度影响因子：喷气发动机推力随速度增加而略有下降（冲压效应），螺旋桨发动机则表现为拉力随速度变化。
- 高度影响因子：空气密度下降导致推力减小，通常用密度比σ = ρ/ρ₀的幂函数表示。
在飞控仿真中，推力模型通常以查表或解析函数形式实现。查表法精度高但需要大量试验数据；解析法便于调参但精度有限。工程实践中常采用混合方式：用解析式描述趋势，用查表修正局部偏差。
4.2 油门响应动态建模
发动机从油门指令变化到推力输出并非瞬时完成，而是存在惯性延迟。这个延迟对飞控系统的影响不可忽视，尤其在需要快速推力响应的场景（如复飞、避障、编队保持）中。
油门响应动态通常用一阶惯性环节近似：
dT/dt = (T_cmd - T) / τ其中，T_cmd为稳态推力指令，τ为发动机时间常数。不同类型发动机的τ值差异显著：
| 发动机类型 | 时间常数τ（秒） | 典型应用 |
|---|---|---|
| 小型电动螺旋桨 | 0.1~0.3 | 多旋翼、小型固定翼 |
| 活塞式发动机 | 0.3~0.8 | 通航飞机、大型无人机 |
| 涡桨发动机 | 0.5~1.5 | 支线客机、运输机 |
| 涡扇/涡喷发动机 | 1.0~3.0 | 喷气式飞机、战斗机 |
在飞控调参时，如果忽略油门响应延迟，会导致仿真中推力变化过快，进而高估飞机的爬升率和加速能力。反之，如果τ设置过大，则可能使推力响应过慢，导致速度控制回路出现振荡。
4.3 螺旋桨发动机简化模型
螺旋桨发动机的推力产生机理与喷气发动机不同，其推力来源于桨叶旋转产生的拉力。简化建模时，通常采用以下两种方法：
方法一：功率固定型模型
假设发动机输出功率恒定，推力与速度成反比：
T = η_p · P / V其中，η_p为螺旋桨效率（通常取0.7~0.85），P为发动机轴功率，V为空速。该模型适用于速度变化范围不大的巡航状态，但在低速或高速极端工况下误差较大。
方法二：拉力系数模型
引入进距比J = V/(nD)，其中n为螺旋桨转速，D为桨径。拉力系数C_T = T/(ρn²D⁴)是J的函数，通过查表或拟合曲线获得：
T = C_T(J) · ρ · n² · D⁴该模型精度更高，但需要螺旋桨性能数据。对于飞控仿真，如果缺乏桨叶数据，可使用简化经验公式：
C_T = C_T0 · (1 - J/J_max)其中，C_T0为静态拉力系数，J_max为失速进距比。
螺旋桨发动机的油门响应除了发动机本身的时间常数外，还需考虑螺旋桨的转动惯量。在变距桨系统中，桨距角变化也会引入额外动态，但飞控建模中通常将其合并到一阶惯性环节中。
4.4 喷气发动机简化模型
喷气发动机的推力模型相对简单，因为其推力不直接依赖空速（在亚音速范围内），主要受油门和高度影响。常用简化模型为：
T = T_max · δ_th · (ρ/ρ₀)^0.7其中，T_max为海平面最大推力，δ_th为油门指令，ρ/ρ₀为密度比。该模型忽略了速度对推力的影响，适用于马赫数小于0.6的飞行包线。
对于更精确的建模，可加入速度修正项：
T = T_max · δ_th · (ρ/ρ₀)^0.7 · (1 - k·M)其中，k为速度衰减系数，典型值为0.1~0.3。该修正反映了高速时进气道总压恢复下降导致的推力损失。
喷气发动机的油门响应时间常数较大，且存在明显的延迟-惯性特性。在飞控仿真中，可用二阶系统近似：
d²T/dt² + 2ζω_n·dT/dt + ω_n²·T = ω_n²·T_cmd其中，ζ为阻尼比（典型值0.7~1.0），ω_n为自然频率（典型值0.5~2 rad/s）。二阶模型能更好地反映喷气发动机的加速-减速不对称性。
4.5 推力对飞行动力学的影响
推力不仅是飞机前进的动力，还通过多种机制影响飞行动力学特性。在飞控建模与调参中，必须考虑以下影响：
（1）纵向配平点移动
推力变化会改变飞机的俯仰力矩平衡。对于推力线低于重心的飞机（大多数常规布局），增加推力会产生低头力矩，需要升降舵配平来补偿。在飞控仿真中，推力变化引起的配平变化会影响升降舵的偏转范围，进而影响操纵裕度。
（2）短周期阻尼变化
推力对短周期模态的影响主要体现在速度阻尼项上。推力随速度增加而减小的发动机（如螺旋桨），会降低速度稳定性，使短周期阻尼减小；推力随速度增加而增加的发动机（如冲压发动机），则增强速度稳定性。
（3）推力-俯仰耦合
当推力线偏离重心时，推力变化会直接产生俯仰力矩，形成推力-俯仰耦合。这种耦合在飞控设计中需要特别注意，因为油门指令可能引起非预期的俯仰响应。在自动油门与自动俯仰控制同时工作时，耦合可能导致控制通道间相互干扰。
（4）非最小相位特性
对于推力线低于重心的飞机，增加推力会产生低头力矩，导致初始下沉，随后才加速爬升。这种“先下沉后上升”的响应具有非最小相位特性，在高度控制回路中可能引起初始反向响应，增加控制难度。
4.6 案例拆解：某小型无人机推力模型辨识与调参
下面通过一个实际案例，展示推力模型在飞控调参中的应用。
背景与约束：某型电动固定翼无人机，翼展2.4m，起飞重量8kg，采用前拉式螺旋桨布局，电机最大功率800W。飞控系统需要实现精确的空速保持功能，但试飞中发现空速控制存在持续振荡，振荡周期约3秒。
初始现象：在自动油门模式下，空速围绕目标值±3m/s振荡，油门指令在30%~70%之间周期性摆动。手动飞行时飞机响应正常，排除气动模型问题。
关键数据：
- 电机时间常数（厂家数据）：τ = 0.15s
- 螺旋桨直径：D = 0.4m
- 巡航速度：V = 20m/s
- 巡航推力需求：T ≈ 8N（估算）
- 飞控油门控制频率：50Hz
分析判断：初步怀疑推力模型中的时间常数设置过小，导致仿真中推力响应过快，而实际电机+螺旋桨系统存在更大的延迟。进一步分析发现，螺旋桨的转动惯量（J_prop ≈ 0.002 kg·m²）在加减速过程中会产生额外延迟，等效时间常数约为0.4s。此外，推力-速度曲线在巡航点附近斜率较大，导致速度反馈增益过高。
实际操作步骤：
- 将推力模型中的时间常数从0.15s调整为0.45s，包含电机和螺旋桨的综合惯性。
- 在推力模型中增加推力变化率限制：最大加速速率20N/s，最大减速速率15N/s。
- 将空速控制回路的比例增益降低30%，积分增益保持不变。
- 在仿真中验证调整后的模型，确认振荡消失。
结果对比：调整后，空速控制误差从±3m/s减小到±0.8m/s，油门指令变化频率明显降低，电机温度下降约15°C（因频繁加减速减少）。
风险与代价：降低比例增益后，空速恢复时间从2秒增加到4秒，在阵风扰动下空速偏差略有增大。后续通过增加前馈补偿（基于空速误差变化率的微分项）来弥补响应速度损失。
复盘经验：推力模型的时间常数不能仅依赖电机厂家数据，必须考虑螺旋桨惯量和调速器动态。在飞控调参前，建议先通过阶跃油门试验辨识实际推力响应特性。