› 牛顿摆 返回
碰撞物理模拟器

牛顿摆 — 弹性碰撞与动量守恒

拉起 k 个球后松手,另一侧会有 k 个球以相同速度飞出。可视化动量柱与动能柱在碰撞前后守恒的过程,并直观体会降低恢复系数 e 时运动衰减的原因。

参数设置
5 个
1 个
1.00
e=1 为理想弹性,e<1 时会衰减(实物钢球约为 0.95~0.99)。
45 °
预设

也可以拖动最左边的球把它拉起,松开后碰撞连锁开始。

实时数值
0
运动中的球
0.00
总动量 p
0.00
总动能
0.0%
能量损失(相对初始值)
—
周期 T [s]
实时牛顿摆
速度矢量 轨迹 碰撞
动量守恒 + 动能守恒 → 进 k 个就出 k 个。 等质量弹性碰撞中,同时满足两个守恒定律的解只有“相同数量的球以相同速度从另一侧飞出”这一个。
动量与动能(碰撞前后守恒)
总动量 p
0.00
总动能
0.00
虚线=初始值的基准线。e=1 时,球每次在最低点碰撞时,p 和 KE 都会回到这条线(碰撞前后守恒)。e<1 时,每次碰撞后 KE 柱能达到的高度都会下降,逐渐衰减。
计算结果
0.00
动能 J/kg
0.00
势能 J/kg
0.00
总能量 J/kg
—
周期 T (s)
与 CAE 的关系 它是碰撞仿真(LS-DYNA、Abaqus/Explicit)的基础模型。恢复系数的建模是缓冲吸能设计和碰撞安全分析的关键,球磨机和冲击成形也受同样的接触力学支配。
理论与主要公式
  • 单摆运动: $\ddot{\theta}= -\dfrac{g}{L}\sin\theta$(半隐式欧拉积分)
  • 弹性碰撞(等质量): $v_1' = \tfrac{1-e}{2}v_1 + \tfrac{1+e}{2}v_2,\; v_2' = \tfrac{1+e}{2}v_1 + \tfrac{1-e}{2}v_2$
  • 守恒定律: $\sum m_i v_i = \text{const}$(动量),$e=1$ 时 $\sum \tfrac12 m_i v_i^2 = \text{const}$(动能)
  • 单摆周期: $T = 2\pi\sqrt{\dfrac{L}{g}}$
  • 每次碰撞的动能损失: $\propto (1-e^2)$

验证:N=5 时拉起 k=1 只有 1 个(其余静止),k=2 时 2 个,k=3 时 3 个从另一侧飞出(e=1 时 p 和 KE 完全守恒)。

什么是牛顿摆

🙋
牛顿摆就是那种桌面玩具吧:把一端的一个球拉起再松手,另一侧就只有一个球飞出去。为什么中间的球不动,只有两端的球在动?
🎓
大致来说,是因为动量和动能都守恒。质量相等的球发生完全弹性碰撞(恢复系数 e=1)时,撞上去的球停下来,它的动量和能量“原封不动”地传给下一个球。在模拟器里把“恢复系数”滑块保持在 1,拉起一个球试试,你会清楚地看到每次碰撞时动量柱和动能柱都回到基准线(=碰撞前后守恒)。
🙋
原来如此!那如果把恢复系数调到小于 1 会怎样?这个“e”表示什么?
🎓
实际中 e=1 的完全弹性碰撞几乎不存在,多为 0.95~0.99。e 是碰撞后与碰撞前相对速度之比,小于 1 就是非弹性碰撞。也就是说,每次碰撞都有一部分动能变成热、声音和变形而损失掉。在模拟器里把 e 降到 0.8 左右试试,应该能观察到:碰撞前后动量不变(p 守恒),但每次碰撞后动能柱能达到的高度都在下降,最终全部停下来。
🙋
明白了!把“拉起的球数 k”设为 2,另一侧真的飞出 2 个球。这也是由守恒定律决定的吗?
🎓
没错。进 k 个就出 k 个——这是同时满足动量守恒和动能守恒的唯一解。假设“只有 1 个球以 k 倍速度飞出”,动量对得上但动能对不上;反过来假设“所有球都以相同的慢速运动”,能量对得上但动量对不上。能同时满足两者的只有“相同数量的球以相同速度飞出”的情况。在 CAE 中做振动和冲击分析时,这些守恒定律也是最基本的基础。

常见问题

把球拖到足够的高度再松开鼠标。拉起角度太小时,碰撞连锁会变弱。另外,恢复系数小于 1 时能量会衰减,连锁容易中断。请把恢复系数设为 1.0 再确认。
恢复系数 e 决定碰撞后的速度比:e=1 为完全弹性碰撞,e=0 为完全非弹性碰撞(粘在一起)。在 CAE 中用于缓冲吸能材料的设计(减小 e 以吸收能量)和碰撞安全部件的性能评价。
本工具把绳长固定为 L=1.0 m,重力加速度固定为 g=9.81 m/s²。周期为 T=2π√(L/g):L 越长周期越长,碰撞间隔变大;g 越大周期越短。无论哪种情况,“进 k 个就出 k 个”的碰撞结果都不会改变。
本工具简化显示了等质量的弹性碰撞,而 LS-DYNA 和 Abaqus/Explicit 用有限元法分析复杂形状和材料之间的接触。恢复系数的概念和动量守恒定律是共通的,可以直观理解冲击分析的基本概念。

现实世界中的应用

碰撞安全设计(汽车、铁路):在车辆碰撞仿真(LS-DYNA 等)中,恢复系数 e 是模拟缓冲吸能材料特性的重要参数。e 的值决定吸收的能量(即对乘员的冲击)。

制造工艺(球磨机、冲压):在粉碎机和金属加工的冲击工艺中,介质(球或锤)之间的碰撞左右着产品质量。碰撞时的能量传递效率正是用本工具所学的物理来分析的。

体育用品开发:高尔夫球杆头与球、网球拍与球等的碰撞中,恢复系数(COR)是性能的关键。有些器材在规则中规定了上限值。

航天(对接、采样返回):航天器对接和小行星探测器采样时的碰撞以极低的相对速度进行。利用非弹性碰撞(e<1)耗散动能的设计很重要。

常见误解与注意事项

首先,你是不是以为“因为有 5 个球,拉起一端的 1 个,另一侧也只会动 1 个”? 其实这只是在初始条件完全对称、并且是完全弹性碰撞(e=1)这种理想情况下才成立。例如一开始同时拉起 2 个球再松开,另一侧就会飞出 2 个球。这是由动量和能量守恒必然决定的。请在模拟器中确认。其次,是参数设置上的陷阱。 把“拉起角”调大后,单摆的周期会比小振幅值 2π√(L/g) 更长(45° 时长约 4%,60° 时长约 7%)。显示的周期是小振幅的近似值,因此要记住,角度大时它与实际的碰撞间隔不同。最后,是“恢复系数 e 只由材料决定”这一误解。 实际上它还会随碰撞速度、温度和表面状态而变化。在 CAE 碰撞分析中,常常不把 e 当作常数,而是使用依赖速度的模型。观察本工具中改变 e 时的行为变化时,诀窍是不要只看成单纯的“衰减”,而要与“能量损失的比例”联系起来。

使用指南

  1. 在 2~9 个范围内设置“球的数量 N”,改变牛顿摆的规模
  2. 把“拉起的球数 k”设为 1~(N−1) 个,确认进 k 个就从另一侧出 k 个
  3. 在 0.60~1.00 范围内改变恢复系数 e,观察完全弹性碰撞(e=1.0)与非弹性碰撞的差别(e=1 时每次碰撞动量柱和动能柱都回到基准线,能量损失保持 0%)
  4. 调整拉起角 θ 和播放速度(0.25×/1×/2×),以慢放/快进观察碰撞连锁。等质量弹性碰撞中动量守恒 Σm·v=常数 始终成立

具体计算示例

本工具的绳长固定为 L=1.0 m,g=9.81 m/s²,小振幅周期为 T=2π√(L/g)≈2.01 秒(与质量无关;拉起角 45° 时实际周期约长 4%)。每个球按 1 kg 处理,能量以 J/kg 显示。把球拉起 45° 时高度为 h=L(1−cos45°)≈0.293 m,势能为 g·h≈2.87 J/kg,在最低点变成同样大小的动能,速度为 √(2gh)≈2.40 m/s。k=2 时总能量加倍,约为 5.75 J/kg。恢复系数 e<1 时,撞击静止球的一次碰撞会损失动能的 (1−e²)/2(e=0.95 时约 4.9%)。5 个球排成一排时,每次摆动中碰撞会连锁发生,因此能量下降得更快。

实务注意事项