我一直觉得“垂直发射火箭的最大高度”是航天动力学里最能暴露物理本质、也最容易把人绕晕的入门题。火箭想要飞得更高要么老老实实忍受重力带来的速度损失要么推力不足、结构失控还没离开发射台就让一切归零——教材当然不会教“原地爆炸”的部分它只会把问题揉进一个编号类似“例题11.1”的题目里单级火箭、垂直发射、给定喷气速度和质量流率求最大高度顺便让你体会什么叫重力损失。这篇文章我就把这类题从运动方程到数值验证完整走一遍既讲清最大高度的推导全过程也把重力损失单独拎出来做一次直观的对比实验。这道题到底能解决什么问题往小说它是理论力学和航天器动力学课程里的经典入门题往大说它决定了你以后看待运载火箭的第一直觉为什么火箭不早点关机为什么发射后要转弯为什么推力不能太小如果你正在学理论力学、空气动力学或航天器动力学或者你在玩模型火箭想知道“发动机参数到底怎么影响最终高度”这篇就是照着教材例题做的完整实战拆解。我自己做这类题的习惯是先建方程再求解析解然后用数值积分反算一遍两边互相验证。这样推导里哪怕有一个符号错误都会在对照中原形毕露。下面正式开始。1. 先把例题的物理直觉建立起来1.1 为什么垂直发射这么“吃亏”垂直发射是所有飞行剖面里重力损失最大的一种。因为整个有动力段推力方向都和重力方向完全相反重力每秒都在从发动机给出的加速度里硬性扣除约9.8m/s²。如果发动机的推重比只有1.1那它可以用来“加速”的净余量其实只有0.1g绝大多数推力都用在和重力拔河上。生活里可以这样理解你拎着一袋20斤的大米爬五楼腿确实在发力但那部分力气并没有变成你水平冲刺的速度它只是用来“对抗重力”。垂直发射就是把这种“抗重力开销”放到最大每升一米都要支付一份势能每一秒都要付出一次速度账面损失。这个现象就是题眼重力损失一般用速度单位m/s度量记作 Δv_gravity。它并不是什么玄学概念而是火箭方程给出的理想速度增量里被重力这个“常驻税种”划走的一笔开销。等算完整道题你会看到它的量级甚至比很多人想象得要大得多。1.2 例题11.1的完整原始设定在某本经典航天动力学教材里这道题以极其朴素的方式出现单级火箭垂直发射忽略空气阻力忽略地球自转重力加速度取常数。我给出一个数值化之后的版本方便复现和检验物理量符号数值初始总质量m05000 kg燃料质量m_fuel4000 kg烧完后结构质量m_e1000 kg有效喷气速度ve2500 m/s质量流率ṁ50 kg/s重力加速度g9.8 m/s²由这些参数可以立刻算出三件事燃烧时间 tb 4000 / 50 80 秒发动机推力 T ṁ·ve 50 × 2500 125000 N初始推重比 T/(m0·g) 125000 / (5000 × 9.8) ≈ 2.55。也就是说这枚理想火箭在地面时推力是自身重量的2.55倍能顺利离地。为什么敢做这么多“理想化假设”因为这道题的目的不是预测真实火箭的轨迹而是把重力损失这个核心因素隔离出来。空气阻力会影响低空段的精确数值地球自转会影响最终入轨几何但如果一上来就把它们全部塞进方程你反而看不清“纯重力到底从你的速度预算里偷走了多少”。例题的价值就是先做减法把次要因素剥干净让主要矛盾暴露在聚光灯下。1.3 力学模型一个质量不断变小的质点火箭的运动方程本质上还是牛顿第二定律只不过质量是时间的函数。以竖直向上为正方向作用在“火箭壳体加剩余燃料”这个控制体上的外力包括发动机推力 T 和重力 mg于是T - mg m·a把 T ṁ·ve 和 m(t) m0 - ṁ·t 代进去得到瞬时加速度a(t) (ṁ·ve) / (m0 - ṁ·t) - g这一步很多人会踩坑他们想直接套高中物理的匀加速公式。但请注意这里的加速度是随时间变化的。原因有两个一是推力恒定而质量不断减少相同推力产生的加速度会越来越大二是重力项始终存在像一个固定税率的“常驻扣款”。也就是说火箭的推力虽然不变但“越来越轻的自己推自己”当然越飞越猛。一个重要的自查点初始加速度必须为正。代入数值初始 a(0) 125000 / 5000 - 9.8 15.2 m/s²说明火箭能离地。如果初始推重比小于1算出来的 a(0) 就会是负的那这枚火箭会先在发射台上“寸步难行”燃料白白烧完也飞不起来——用标题的话说想飞得更高至少先把重力和推力的胜负关系搞清楚。2. 最大高度的完整推导从运动方程到623.6km2.1 动力段速度与高度都要靠积分把加速度对时间积分一次得到动力段速度v(t) ve·ln(m0 / m(t)) - g·t其中 m(t) m0 - ṁ·t。这个公式的第一项正是齐奥尔科夫斯基火箭方程给出的理想速度增量第二项就是累计的重力损失。你可以这么看火箭方程告诉你“如果不被重力拖后腿你能取得多少速度”而重力项告诉你“现实里每秒钟要交多少税”。再积分一次得到高度。这一步最容易出错因为被积函数里有对数项。我直接给出换元法的中间结果对 h(t) ∫v(t)dt令 u m0 - ṁ·τ整理后得到h(t) ve·[ t - (m(t)/ṁ)·ln(m0/m(t)) ] - (1/2)·g·t²这个公式看起来不友好但每个部分都有物理意义。方括号代表“理想速度随时间的积分结果”后面的半重力项则是重力在位移上的累计后果。代入 t 80sm(80) 1000kg可以得到v1 2500·ln(5000/1000) - 9.8×80 3239.6 m/sh1 2500×(80 - (1000/50)×ln5) - 0.5×9.8×6400 88168 m ≈ 88.2 km也就是说发动机工作80秒后火箭飞了88.2公里高速度约3.24km/s。它此时还在快速上升远没到最高点。要注意的是有动力段结束时加速度已经拉到 a(80) (125000 - 1000×9.8) / 1000 115.2 m/s²约11.8g。这种过载在真实载人火箭里是绝对不允许的但教材例题里它只是理想化模型的副产品。你真正该关心的是关机点速度而不是把过载也当作可接受的设计指标。2.2 滑行段把剩余速度换成高度燃料烧完后发动机不再提供推力火箭只剩一个1000kg的“空壳”做惯性滑行。此时质量恒定受力只有重力运动变成标准的匀减速直线运动。设滑行段上升高度为 h2最高点速度必须降为0。由能量守恒或者直接套匀减速公式(1/2)·m·v1² m·g·h2于是 h2 v1² / (2g)。代入 v1 3239.6 m/sh2 3239.6² / (2×9.8) ≈ 535.5 km最大高度为H h1 h2 88.2 535.5 623.7 km ≈ 624 km很多初学朋友看到这个数字会愣住燃料才烧80秒为什么还能再爬升535公里答案在于滑行段几乎是在“真空边缘以上”飞行。虽然这里把大气阻力忽略了但物理直觉是对的高空空气稀薄火箭以3.2km/s的速度冲上去惯性足以把它送到很远的高度。教材故意选这个题目就是让你建立一个观念最大高度出现在速度降为零的时刻而不是发动机停机的时刻。2.3 把结果放进数量级坐标系里623.7公里是个什么概念国际空间站轨道高度约400公里地球半径约6371公里。所以这枚例题火箭的亚轨道顶点已经越过了近地轨道高度甚至接近很多低轨卫星的运行高度。当然它没有横向速度不会入轨只是会落回来。这时候值得做一个反事实对比如果完全没有重力损失理想末速度会是 v_ideal ve·ln(m0/m_e) 2500×ln5 ≈ 4023.6 m/s。用同样的重力场做滑行分析对应的最大高度将是h1_ideal v_ideal²/(2g) 119.5 km 825.9 km 945.4 km也就是说重力损失让这枚例题火箭的最大高度从接近945公里降到624公里损失了约321公里。这个差距不是“修正量”而是主导量。这恰恰是这类例题最想让你看见的事重力损失不是小账而是大账。3. 重力损失解析到底损失了什么东西3.1 核心公式Δv_g g·tb从一个最简单的动作开始。竖直飞行中重力对加速度的贡献每时每刻都是 -g只要发动机在燃烧这9.8m/s²的“税”就一直扣。累计 tb 秒速度损失就是Δv_g g·tb代入本题Δv_g 9.8×80 784 m/s。这就是为什么说重力损失和“烧了多久”强相关而和喷气速度、质量比并没有直接关系。同一枚火箭如果把质量流率调大一倍燃烧时间会从160秒缩到80秒重力损失直接从1568m/s减半到784m/s。这个结论对工程选型极其有用想降低重力损失最直接的办法是缩短燃烧时间而不是一门心思提高比冲。当然真实任务里推力方向会随姿态变化g也会随高度变化。更严谨的写法是把它扩展成矢量积分 Δv_g ∫(g·sinγ)dt其中γ是速度矢量与当地水平面的夹角。垂直发射时γ 90°sinγ 1就退化成本题的一维形式。这也是“垂直发射重力损失最大”的原因sinγ取到了最大值。3.2 快烧还是慢烧一次控制变量的对比实验为了把重力损失的分量看得更清楚我把例题改成方案B保持初始质量、燃料质量、喷气速度都不变只把质量流率从50kg/s降到25kg/s。于是理想速度增量不变但燃烧时间变成160秒初始推重比只剩1.28。参数方案A快速燃烧方案B慢速燃烧质量流率50 kg/s25 kg/s燃烧时间80 s160 s初始推重比2.551.28理想速度增量4023.6 m/s4023.6 m/s重力损失784 m/s1568 m/s动力段末速度3239.6 m/s2455.6 m/s动力段末高度88.2 km113.6 km滑行段上升高度535.5 km307.6 km最终最大高度623.7 km421.2 km这个表的信息量很大。慢速燃烧方案前期推力小、加速度低在动力段里“磨蹭”的时间更长因此动力段高度反而更高但它最终总高度却少了200多公里。原因就是关机点速度差太多一个带着3.2km/s关机一个只带2.5km/s关机。速度差距在滑行高度上被平方放大于是重力损失在总高度上体现出压倒性影响。这也解释了为什么运载火箭的第一级几乎都是“大力出奇迹”的设计宁可发动机重一点、结构粗一点也要在尽量短的时间内把火箭推离稠密大气用较小的重力损失换取更高的净速度和弹道收益。3.3 真实火箭为什么不是一直垂直加速既然垂直发射重力损失最大为什么火箭不从一开始就水平加速答案有两层。低空大气太稠密水平速度会产生巨大动压和热流必须先垂直冲过稠密区离开稠密大气后火箭也不能立刻大角度转弯而是采用“重力转弯”程序——推力方向始终跟随速度方向让重力负责慢慢把速度矢量从竖直转向水平。重力转弯的本质是让重力从“纯敌人”变成“教练”它一方面继续造成重力损失另一方面不断改变速度方向使推进剂产生的速度增量越来越多的朝水平方向注入。最终近地轨道任务的重力损失仍然有几公里每秒但相比纯垂直发射已经小很多。例题里的垂直发射故意把这种工程复杂性全部剥掉只剩一个清洁的物理事实重力和推力之间最简单的拔河。4. 用数值积分验证解析解别让推导变成自说自话4.1 用Python做一维弹道积分解析推导再漂亮也值得用一个数值积分来“对账”。我习惯把运动方程直接交给欧拉法虽然它精度不如RK4但对于这种光滑的一维问题只要步长取足够小已经足够验证解析结果。dt 0.01 h 0.0 # 高度单位m v 0.0 # 速度单位m/s m 5000.0 # 当前质量单位kg m_empty 1000.0 mdot 50.0 ve 2500.0 g 9.8 while v 0: if m m_empty: thrust mdot * ve a (thrust - m * g) / m m max(m_empty, m - mdot * dt) else: a -g v_new v a * dt h (v v_new) / 2 * dt # 梯形积分降低位移误差 v v_new if v 0: break print(最大高度约, h / 1000, km)这段代码的核心逻辑很简单每个时间步重新计算当前质量和当前加速度先更新速度再用梯形公式更新高度燃料耗尽后加速度只剩 -g当速度降到0循环停止此刻的高度就是最大高度。我故意没有在动力段结束时停止循环因为最高点永远在关机之后。4.2 数值结果和解析结果对账按上面代码用 dt0.01s 跑一遍最大高度会落在623公里附近和解析解623.7km的差距通常只有几公里量级并且步长越小越接近解析值。如果不做梯形积分只用朴素欧拉 h v*dt误差会比梯形版本明显大。这本身也是数值方法的一个教训位移积分尽量用梯形或更高阶方法不然算这种高空滑行问题会出现肉眼可见的偏差。更有意思的是打印“达到最高点的时间”。解析上滑行段时间 t_coast v1/g 3239.6/9.8 ≈ 330.6秒加上动力段80秒总时间约410.6秒。这意味着在发动机熄火后火箭还要继续向上漂5分半钟才抵达顶点。许多初算者拿到623km这个数字时总以为它是烧完燃料那一瞬间的高度——错。火箭80秒烧完燃料时只有88km真正的顶点在约6分50秒时才出现。所以看遥测数据时“关机点速度”比“关机点高度”更能判断最终弹道高度。4.3 什么时候该相信这个例题解析模型做了大量简化我把边界条件讲透。忽略空气阻力对低空段影响极大真实火箭在30km以下会受到明显动压和阻力最大高度会显著缩水重力加速度取常数而624km高度上实际重力已经降到约8.1m/s²滑行段后半段会比模型“更省力”但影响量级在10%以内地球自转给火箭的初始东向速度在垂直发射模型里用不上但如果算入轨它是一项免费的横向速度增量。所以例题的价值不在精确预报而在于隔离出重力损失这一单一因素让你看清量级。做工程估算时解析解当初步骨架数值积分加气动模型当修正两者配合才是正确姿势。5. 常见误区与排查实录5.1 最容易被绕进去的五个坑常见误区正确理解最大高度出现在燃料烧完瞬间最高点出现在关机后的滑行段速度减到0才到顶点直接用 h½at² 算动力段高度加速度随质量变化必须积分或逐段累加重力损失 能量白白浪费部分转化为势能但在速度预算口径里确实是损失推重比越大一定越好理想例题里是真实工程受结构过载、发动机增重约束忽略初始推重比是否大于1推重比≤1时初始加速度≤0火箭连离地都做不到这五条我都在实际教学和自算中见过。第一条最经典拿到623km后第一反应往往是“燃料烧完才88km后面的高度哪来的”。答案就在滑行段。第二条是计算习惯问题很多初学朋友惯性套用匀变速公式结果动力段高度差了近30km再一平方放大总高度千差万别。第三条值得多说一句。“重力损失”损失的是火箭方程里的速度预算不是热力学意义上的能量。被重力抵消的那部分推力确实作为势能存进了火箭的高度里——这就是为什么火箭能到达高位置。但在判断能否入轨、能飞多高的“速度指标”上这笔开销无法参与最终速度所以工程上仍然叫损失。它不是无用功而是预算被划走。5.2 快速自查你的推导是否靠谱你算完后可以用三个办法快速检查量纲检查任何公式中的每一项都必须是m/s或m不能出现m²/s之类混搭项。极限检查如果质量流率趋近于0燃烧时间趋近无穷实际末速度应该趋近一个很糟糕的值甚至无法起飞如果喷气速度趋近无穷重力损失占比趋近0。能量复核滑行段最高点处½·v1²应该恰好等于g·h2把这个代入算一遍如果两边不相等说明前面的速度或高度公式有问题。还有一个工程习惯在跑数值积分时把动力段结束时刻的v、h、m打印出来和解析公式算出的v1、h1对比。如果这三项都对得上后面滑行段基本不会错。我每次推导的最后一公里都是靠这种“打印中间状态”的习惯来兜底。6. 最后再聊点个人经验我第一次算这类题时也觉得“既然有空气阻力那重力损失这事肯定没那么严重”。直到我把方案A和方案B放一起算完才意识到重力损失才是那个能把射高砍掉三分之一的大头空气阻力在低空段的贡献反而更像是叠加修正。后来我帮人估算探空火箭和模型火箭的最大高度起始步骤永远是先评估发动机推重比再算燃烧时间然后立刻估计重力损失的量级。判断一枚火箭能飞多高与其去精算阻力系数不如先把重力损失这笔账记牢。如果你也打算拿这个流程去估算自己的模型火箭我建议在发射前做一次地面静推力测试至少测出平均推力和真实燃烧时间然后代入这套框架算一遍。你会发现“为什么同样的发动机有人装出2km射高有人只能上300m”——推力一致的情况下多出来的那部分火箭重量要么把初始推重比压到1附近要么把燃烧时间拖长到重力损失爆炸。这道例题教给我的说到底就一句话想飞得更高你必须先和重力把账算清楚清楚到既不幻想绕开它也不用蛮力挑战它。