参数设置
预设场景
显示选项
拖动最左侧的球向上提起，松手即可开始。
📐 物理模型
- 摆动方程： $\ddot{\theta}= -\dfrac{g}{L}\sin\theta$（Verlet 积分）
- 弹性碰撞（等质量）： $v_1' = e \cdot v_2,\; v_2' = e \cdot v_1$
- 单摆周期： $T = 2\pi\sqrt{\dfrac{L}{g}}$
- 动量守恒： $\sum m_i v_i = \text{const}$
CAE 应用：
这是碰撞仿真（LS-DYNA、Abaqus/Explicit）的基础模型。恢复系数的建模在冲击吸能材料设计和汽车碰撞安全分析中不可或缺。球磨机和冲压工艺中的接触力学与此相同。
什么是牛顿摆
🧑🎓
为什么我拉起一个球，松开后，总是对面弹出一个球，而不是两个球各飞一半速度呢？
🎓
简单来说，这是为了同时满足“动量守恒”和“动能守恒”这两个铁律。在实际工程中，比如两个完全相同的钢球正面碰撞，它们不会“分享”速度。你可以在这个模拟器里试试：先拉起最左边一个小球，然后松开，观察碰撞后最右边小球的速度。试着把“恢复系数”滑块拉到最右边（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<1$时为非弹性碰撞，部分动能转化为内能。
现实世界中的应用
汽车碰撞安全分析：在LS-DYNA或Abaqus/Explicit等CAE软件中，恢复系数是定义接触界面的关键参数。工程师通过调整它来模拟车辆碰撞时保险杠、车架的吸能效果，从而优化安全设计。
球磨机与粉碎工艺：在矿业和材料加工中，球磨机利用钢球从高处落下撞击物料。模拟钢球间的碰撞和能量传递，对于设计高效节能的粉碎设备至关重要。
精密冲压与成型：在金属冲压过程中，模具与板料的接触碰撞类似于非弹性碰撞。通过控制“恢复”行为（即回弹），可以预测和补偿零件的最终形状，提高制造精度。
运动器材与生物力学：分析网球拍击球、高尔夫球杆击球，甚至人体关节的冲击（如跑步时脚掌落地），其核心物理模型都离不开碰撞理论和动量守恒定律。
常见误解与注意事项
首先，你是否误以为“因为有5个小球，所以拉动一端的一个球时，另一侧也只会弹出一个球”？ 实际上，这仅在初始条件完全对称且为完全弹性碰撞（e=1）的理想情况下才会发生。例如，如果最初同时拉动并释放两个球，则另一侧会弹出两个球。这是由动量守恒和能量守恒定律必然决定的。你可以在模拟器中尝试验证。其次，参数设置的陷阱。 如果将“重力”设为0，小球碰撞后将沿直线匀速运动而不再返回。此时系统已不再是摆锤模型。在实际工程中，将参数设为极端值导致模型失效是典型问题。最后，关于“恢复系数e仅由材料决定”的误解。 实际上e值会受碰撞速度、温度及表面状态影响。在CAE碰撞分析中，常将e视为与速度相关的变量而非常数。使用本工具时，观察改变e值引发的行为变化时，建议将其与“能量损失比例”而非简单的“衰减”关联起来。
相关工程领域
这个简易模拟器背后的计算原理正是多体动力学领域的基础。例如汽车悬架系统或机器人手臂关节等由铰链连接、相互干涉的多刚体系统运动分析。此外，小球碰撞计算是离散元法的核心技术，常用于模拟矿石破碎流程、谷仓内颗粒流动及药品混合状态。而用于求解摆锤运动方程的“Verlet积分法”，同样是分子动力学中原子轨迹计算的基础算法。换言之，从玩具模拟到纳米级物质模拟，其核心思想一脉相承。碰撞安全设计（碰撞仿真）虽需处理复杂“面接触”分析而非此处“点”碰撞，但其出发点正源于此。
进阶学习建议
第一步建议通过公式追踪“守恒定律”。本工具中质量均相同，但可尝试推导质量不同两球（如m1=1kg, m2=3kg）的正面弹性碰撞，联立动量守恒与动能守恒方程求解。若能推导碰撞后速度，理解将大幅深化。其次，建议学习摆锤微分方程的数学背景。$$ \ddot{\theta}= -\dfrac{g}{L}\sin\theta $$ 该式在小角度条件下（sinθ≒θ）可近似为“简谐振动”方程，这是线性化的第一步。通过模拟器改变“摆线长度”时，观察周期变化并与近似公式 $$ T = 2\pi \sqrt{L/g} $$ 对比，是很好的练习。若想进一步深入，推荐研究“角动量守恒”与“非对心碰撞（斜碰撞）”。台球碰撞或球体斜碰墙壁等现象更贴近实际情况。通过本工具掌握基础后，可尝试挑战能在二维平面移动小球的更高级模拟器。