4. 发动机推力模型
发动机推力模型，说白了就是回答一个问题：这发动机到底能推多大力？
我在做飞行器设计时，最头疼的就是推力建模。你想想看，发动机推力不是个固定值，它随着飞行高度、速度、油门开度一直在变。搞不定这个，后面的轨迹仿真全是白搭。
这一章，我带你搭建两种最常见的推力模型：航空发动机（涡扇/涡喷）和火箭发动机。咱们不搞复杂的CFD，用工程上够用的简化模型。
4.1 航空发动机推力特性
先说涡扇和涡喷。这两种发动机的原理我就不展开了，咱们直接看影响推力的两个核心因素：高度和马赫数。
4.1.1 高度影响
高度越高，空气越稀薄。发动机需要吸入空气来燃烧，空气少了，推力自然下降。
我习惯用一个简单的高度因子来修正海平面推力：
# 高度因子计算
def altitude_factor(h):
"""
h: 高度，单位米
返回：0~1之间的高度因子
"""
# 对流层内（0~11000m）按大气密度比例
if h <= 11000:
rho_ratio = (1 - h/44330.0) ** 4.256
else:
# 平流层简化处理
rho_ratio = 0.297 * np.exp((11000 - h)/6340)
return max(0.1, rho_ratio) # 限制最小值，防止除零我的经验：这个模型在0~15000米范围内误差在5%以内。超过这个范围，建议查标准大气表。
4.1.2 马赫数影响
马赫数对推力的影响比较复杂。低速时，随着速度增加，进气量增大，推力反而上升。但到了跨音速区，激波损失开始显现，推力会下降。
我整理了一个经验公式，在多个项目中验证过：
def mach_factor(Ma):
"""
Ma: 马赫数
返回：马赫数因子
"""
# 典型涡扇发动机的马赫数特性
if Ma < 0.8:
factor = 1.0 + 0.3 * Ma
elif Ma < 1.2:
factor = 1.24 - 0.3 * (Ma - 0.8)
else:
factor = 1.12 - 0.1 * (Ma - 1.2)
return max(0.5, factor)注意：这个公式适用于典型的亚音速涡扇发动机。如果是超音速涡喷，曲线会完全不同。我曾经在某个项目中直接套用，结果仿真结果完全不对——后来才发现发动机类型搞错了。
4.1.3 综合推力模型
把高度因子和马赫数因子乘起来，再乘以海平面最大推力，就得到了当前状态下的推力：
def turbofan_thrust(h, Ma, throttle=1.0):
"""
涡扇发动机推力模型
h: 高度(m)
Ma: 马赫数
throttle: 油门开度(0~1)
返回：推力(N)
"""
# 海平面最大推力（示例值）
T_max_sl = 120000 # 120kN
# 计算因子
f_h = altitude_factor(h)
f_m = mach_factor(Ma)
# 油门修正（非线性）
f_throttle = throttle ** 1.2
# 最终推力
T = T_max_sl * f_h * f_m * f_throttle
return T4.2 火箭发动机推力模型
火箭发动机就简单多了。它自带氧化剂，不依赖大气。所以高度影响主要是喷管膨胀的问题。
4.2.1 真空推力与海平面推力
火箭发动机有两个关键参数：
- 海平面推力：在地面测试时测得的推力
- 真空推力：在真空中能产生的推力
两者之间的关系：
def rocket_thrust(h, T_vac, T_sl):
"""
火箭发动机推力模型
h: 高度(m)
T_vac: 真空推力(N)
T_sl: 海平面推力(N)
返回：当前高度推力(N)
"""
# 大气压力随高度变化
p0 = 101325 # 海平面大气压(Pa)
if h <= 11000:
p = p0 * (1 - h/44330.0) ** 5.256
else:
p = p0 * 0.223 * np.exp((11000 - h)/6340)
# 推力随高度线性变化
T = T_vac - (T_vac - T_sl) * (p / p0)
return T核心要点：火箭发动机的推力随高度增加而增大。因为外部大气压降低，喷管出口压力与外部压力的压差增大，推力自然就上去了。
4.2.2 比冲模型
比冲是衡量发动机效率的关键指标。对于火箭发动机，比冲也随高度变化：
def specific_impulse(h, Isp_vac, Isp_sl):
"""
比冲随高度变化模型
"""
p0 = 101325
if h <= 11000:
p = p0 * (1 - h/44330.0) ** 5.256
else:
p = p0 * 0.223 * np.exp((11000 - h)/6340)
Isp = Isp_vac - (Isp_vac - Isp_sl) * (p / p0)
return Isp4.3 两种发动机的对比
| 特性 | 涡扇/涡喷发动机 | 火箭发动机 |
|---|---|---|
| 依赖大气 | 是 | 否 |
| 推力随高度变化 | 下降 | 上升 |
| 推力随速度变化 | 复杂（先升后降） | 基本不变 |
| 比冲范围 | 3000~6000s | 250~450s |
| 适用场景 | 大气层内飞行 | 大气层内外 |
我的建议：做飞行器仿真时，先把发动机类型搞清楚。如果是导弹，大概率用火箭发动机；如果是飞机，用涡扇/涡喷。别搞混了——我见过有人用火箭发动机模型去仿真客机，结果推力随高度增加，完全反了。
4.4 本章知识体系
下面这张图帮你理清本章的核心逻辑：
4.5 工程实践建议
最后，给你几个我在项目中踩过的坑：
- 别忘了油门响应延迟。发动机不是瞬间就能输出你想要的推力，有个动态过程。我一般加一个一阶惯性环节来模拟。
- 高度因子别用线性近似。大气密度不是线性变化的，用指数模型更准。
- 马赫数因子要实测数据校准。经验公式只能给个大概，真要精确仿真，得找发动机厂家要数据。
我曾经踩过的坑：在某次飞行器性能评估中，我直接用海平面推力算全弹道，结果仿真出来的爬升率比实际高了一倍。后来才发现，到了高空推力已经掉到海平面的40%了。从那以后，我每次做仿真都会先检查推力模型的高度修正。
好了，发动机推力模型就讲到这里。代码可以直接拿去用，但记得根据你的具体发动机参数调整一下系数。
公众号:蓝海资料掘金营，微信deep3321