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

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

拖动小球提起后松手,观察弹性碰撞链如何传递动量。调节重力加速度、绳长和恢复系数,直观感受能量的流动与转化。

参数设置
重力加速度 g
m/s²
绳长 L
m
恢复系数 e
速度倍率
×
预设场景
显示选项

拖动最左侧的球向上提起,松手即可开始。

计算结果
0.00
动能 J
0.00
势能 J
0.00
总能量 J
周期 T (s)
Cradle
能量时程曲线
理论与主要公式
  • 摆动方程: $\ddot{\theta}= -\dfrac{g}{L}\sin\theta$(Verlet 积分)
  • 弹性碰撞(等质量): $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$
  • 单摆周期: $T = 2\pi\sqrt{\dfrac{L}{g}}$
  • 动量守恒: $\sum m_i v_i = \text{const}$

什么是牛顿摆

🙋
为什么我拉起一个球,松开后,总是对面弹出一个球,而不是两个球各飞一半速度呢?
🎓
简单来说,这是为了同时满足“动量守恒”和“动能守恒”这两个铁律。在实际工程中,比如两个完全相同的钢球正面碰撞,它们不会“分享”速度。你可以在这个模拟器里试试:先拉起最左边一个小球,然后松开,观察碰撞后最右边小球的速度。试着把“恢复系数”滑块拉到最右边(e=1),这就是完美的弹性碰撞,你会看到能量几乎没有损失。
🙋
诶,真的吗?那如果我拉起两个球再松开,是不是对面也会弹出两个球?
🎓
没错!这正是牛顿摆最神奇的地方。拉起几个球,对面就会弹出几个球,这是质量和能量守恒的必然结果。你可以在模拟器里验证一下:同时拉起左边两个小球,然后松开。你会发现中间的小球好像“隐形”了一样,动量直接传递给了最右边的两个球。试着把“重力加速度”调小一点,比如模拟月球环境,你会发现整个碰撞过程变慢了,但传递的规律不变。
🙋
那如果我把“恢复系数”调小一点会怎么样?现实中不可能有e=1的完美碰撞吧?
🎓
你的观察很敏锐!工程现场常见的是非完全弹性碰撞,恢复系数e通常小于1,比如实际金属球在0.95到0.99之间。e<1意味着碰撞会损失一部分动能,转化为热量或变形能。你现在就可以在模拟器里试试:把恢复系数滑块从1慢慢往左拖,比如拖到0.8,然后拉起一个球。你会看到碰撞后,对面球摆起的高度明显降低了,而且整个系统摆动几下后就慢慢停下来了,这就是能量损耗的直观体现。

物理模型与关键公式

小球的摆动遵循单摆运动规律,其角加速度由重力沿切线方向的分力提供:

$$\ddot{\theta}= -\dfrac{g}{L}\sin\theta$$

其中,$\ddot{\theta}$是角加速度,$g$是重力加速度(可调节),$L$是绳长(可调节),$\theta$是摆角。模拟器使用Verlet积分法对这个方程进行数值求解,从而计算出小球每一帧的精确位置。

当两个等质量小球发生对心碰撞时,其碰撞后的速度由恢复系数$e$决定:

$$v_1' = e \cdot v_2, \quad v_2' = e \cdot v_1$$

这里$v_1, v_2$是碰撞前的速度,$v_1', v_2'$是碰撞后的速度。当$e=1$时为完全弹性碰撞,满足动量和动能双重守恒;$e\lt 1$时为非弹性碰撞,部分动能转化为内能。

现实世界中的应用

汽车碰撞安全分析:在LS-DYNA或Abaqus/Explicit等CAE软件中,恢复系数是定义接触界面的关键参数。工程师通过调整它来模拟车辆碰撞时保险杠、车架的吸能效果,从而优化安全设计。

球磨机与粉碎工艺:在矿业和材料加工中,球磨机利用钢球从高处落下撞击物料。模拟钢球间的碰撞和能量传递,对于设计高效节能的粉碎设备至关重要。

精密冲压与成型:在金属冲压过程中,模具与板料的接触碰撞类似于非弹性碰撞。通过控制“恢复”行为(即回弹),可以预测和补偿零件的最终形状,提高制造精度。

运动器材与生物力学:分析网球拍击球、高尔夫球杆击球,甚至人体关节的冲击(如跑步时脚掌落地),其核心物理模型都离不开碰撞理论和动量守恒定律。

常见误解与注意事项

首先,你是否误以为“因为有5个小球,所以拉动一端的一个球时,另一侧也只会弹出一个球”? 实际上,这仅在初始条件完全对称且为完全弹性碰撞(e=1)的理想情况下才会发生。例如,如果最初同时拉动并释放两个球,则另一侧会弹出两个球。这是由动量守恒和能量守恒定律必然决定的。你可以在模拟器中尝试验证。其次,参数设置的陷阱。 如果将“重力”设为0,小球碰撞后将沿直线匀速运动而不再返回。此时系统已不再是摆锤模型。在实际工程中,将参数设为极端值导致模型失效是典型问题。最后,关于“恢复系数e仅由材料决定”的误解。 实际上e值会受碰撞速度、温度及表面状态影响。在CAE碰撞分析中,常将e视为与速度相关的变量而非常数。使用本工具时,观察改变e值引发的行为变化时,建议将其与“能量损失比例”而非简单的“衰减”关联起来。

使用指南

  1. 在gravNum字段输入重力加速度值(9.8 m/s²为地球标准值,可调至1.6 m/s²模拟月球环境),gravParam选择单位
  2. 通过lenNCNum设置摆球绳长(范围0.5-2.0 m),影响摆动周期T=2π√(L/g)的计算结果
  3. 调节restNum恢复系数(0-1之间),0代表完全非弹性碰撞,1代表完全弹性碰撞,观察碰撞后球体反弹高度变化
  4. 启动模拟后,实时监测动能、势能、总能量数值变化及周期参数

具体计算示例

设置钢球质量m=0.5 kg、绳长L=1.2 m、g=9.8 m/s²、恢复系数e=0.95。初始释放高度h=1.2 m时,势能Ep=0.5×9.8×1.2=5.88 J。碰撞前速度v=√(2gh)=4.85 m/s,碰撞后速度v'=0.95×4.85=4.61 m/s。周期T=2π√(1.2/9.8)=2.20 s。每次碰撞后总能量衰减约(1-0.95²)×100%=9.75%。

实务注意事项

  1. 当恢复系数e<0.8时,摆球在5-8个周期内能量衰减至初值的10%以下,适用于模拟真实材料碰撞(橡胶e≈0.6,钢球e≈0.92)
  2. 绳长精确度影响周期计算:1 cm误差导致周期偏差约0.4%,模拟工程应用时需验证测量数据
  3. 重力加速度变化会改变动能与势能的相互转化速率,月球环境(g=1.62 m/s²)下周期约为地球的2.47倍
  4. 恢复系数与材料温度相关,低温钢球e值可增加0.03-0.05,高温环境需重新校准参数