终端速度(当 $dv/dt = 0$ 时):
$$v_t = \sqrt{\frac{2mg}{C_d \rho A}}$$真空落体($C_d=0$):$v = v_0 + gt$,$h = h_0 - v_0 t - \frac{1}{2}gt^2$
数值积分:用Euler法(Δt = 0.01s)计算全部数据点后绘图
改变质量、空气阻力系数和截面积,对自由落体进行时间积分。比较真空落体与空气中落体,体验终端速度的形成过程。
终端速度(当 $dv/dt = 0$ 时):
$$v_t = \sqrt{\frac{2mg}{C_d \rho A}}$$真空落体($C_d=0$):$v = v_0 + gt$,$h = h_0 - v_0 t - \frac{1}{2}gt^2$
数值积分:用Euler法(Δt = 0.01s)计算全部数据点后绘图
这个模拟器最核心的控制方程,描述了物体在重力与空气阻力共同作用下的运动。空气阻力与速度的平方成正比,这是高速运动时的典型模型。
$$m\frac{dv}{dt}= mg - \frac{1}{2}C_d \rho A v^2$$m:物体质量 (kg)
dv/dt:加速度 (m/s²)
g:重力加速度 (≈9.8 m/s²)
Cd:阻力系数,取决于物体形状
ρ:空气密度 (kg/m³)
A:物体在运动方向上的投影截面积 (m²)
v:瞬时速度 (m/s)
等式右边第一项是向下的重力,第二项是向上的空气阻力。
当加速度为零(dv/dt=0)时,重力与阻力平衡,速度达到稳定值,即终端速度。这是上面控制方程的一个特解。
$$v_t = \sqrt{\frac{2mg}{C_d \rho A}}$$vt:终端速度 (m/s)。公式表明,终端速度与质量的平方根成正比,与阻力系数、空气密度和截面积的平方根成反比。这就是为什么重而流线型的物体(如子弹)终端速度极高,而轻且迎风面积大的物体(如蒲公英种子)终端速度很低。
汽车工程与空气动力学:在CFD(计算流体力学)仿真中,工程师会计算不同车体设计的阻力系数Cd。这个仿真器中的原理与之相通,通过调整Cd和截面积A,可以直观理解这些参数如何影响物体在流体中的最终运动状态,用于初步验证设计。
包装跌落测试与产品安全:在CAE软件(如Abaqus, LS-DYNA)中进行跌落冲击分析时,需要给产品设定一个合理的初始撞击速度。这个速度往往就是考虑空气阻力后的终端速度,而非真空自由落体速度,这样模拟结果才更真实。
航空航天与跳伞运动:航天器返回舱或跳伞运动员的下落轨迹规划,严重依赖于终端速度的计算。不同高度下空气密度ρ的变化,会直接影响终端速度,从而决定开伞时机、着陆冲击力等关键因素。
运动器材设计与环境科学:高尔夫球表面的凹坑、滑雪运动员的姿势都是为了优化阻力系数Cd和截面积A,从而改变空气阻力,影响运动轨迹或终端速度。同样,研究雨滴、沙尘的沉降速度也基于此模型。
在使用本模拟器时,有几个初学者尤其容易陷入的误区。首先,人们常认为“质量变为两倍则终端速度也会变为两倍”,但这是错误的。观察终端速度公式 $v_t = \sqrt{2mg / (C_d \rho A)}$ 可以发现,速度与质量的平方根 $\sqrt{m}$ 成正比。也就是说,即使质量变为4倍,终端速度也只会变为2倍。例如,比较100kg与400kg的物体,速度差仅为2倍。这种非线性关系有违直觉,需要特别注意。
其次,参数的单位和尺度应基于现实设定。请不要随意将质量设为“1”、截面积设为“1”就运行模拟,然后仅得出“哦”这样的结论就结束。若模拟雨滴,质量应约为0.001kg(1克),直径为数毫米。若是跳伞者,质量70kg、截面积根据姿势在0.5~1.0 m²左右较为合理。若将工具的滑块调至极端值,则可能得到不现实的结果(例如时速数千公里的终端速度)。实际工程中使用的数据,应始终考虑现实的数量级。
最后,请勿忘记本模型基于“恒定密度的静止流体中”这一理想化条件。实际上,高度降低时空气密度 $\rho$ 会增加;物体高速下落时周围气流可能变得紊乱,导致阻力系数 $C_d$ 本身随速度变化。此外,若物体发生旋转或以不稳定姿态下落,运动将更为复杂。本模拟器旨在帮助学习基础知识,并未包含这些高阶现象。首先扎实理解基础模型,进而认识其局限性,是通往下一步的捷径。