一维碰撞的动量能量分析
在经典力学中，碰撞是一个极其重要的物理过程。无论是微观世界的粒子散射，还是宏观世界的交通事故分析，动量守恒与能量守恒都是揭示碰撞前后状态变化的黄金法则。本文将深入探讨一维碰撞（即物体在同一直线上的碰撞）的数学推导与物理本质，重点分析完全弹性碰撞与完全非弹性碰撞。
在一维碰撞问题中，我们通常考察一个封闭系统（不受外力或外力合力为零）。设两个质点的质量分别为 $m_1$ 和 $m_2$，碰撞前的初始速度分别为 $u_1$ 和 $u_2$，碰撞后的末速度分别为 $v_1$ 和 $v_2$。
1. 动量守恒定律
根据牛顿第三定律，系统在碰撞过程中所受的合外力为零，因此系统的总动量保持不变：
$$m_1 u_1 + m_2 u_2 = m_1 v_1 + m_2 v_2$$
动量守恒适用于所有类型的碰撞（无论是弹性还是非弹性）。
2. 机械能守恒与恢复系数
虽然动量总是守恒的，但机械能却不一定。在碰撞过程中，部分机械能可能会转化为热能、声能或内能。为了定量描述碰撞的能量损失，我们引入了恢复系数（Coefficient of Restitution, 简记为 $e$）：
$$e = \frac{v_2 - v_1}{u_1 - u_2}$$
恢复系数 $e$ 的取值范围为 $0 \le e \le 1$：
- 当 $e = 1$ 时，为完全弹性碰撞，系统动能无损失。
- $0 < e < 1$ 时，为非弹性碰撞，部分动能转化为其他形式的能量。
- 当 $e = 0$ 时，为完全非弹性碰撞，碰撞后两物体粘合在一起，以相同的速度运动。
完全弹性碰撞的数学推导
在完全弹性碰撞中，不仅动量守恒，系统的总动能也守恒：
$$\frac{1}{2} m_1 u_1^2 + \frac{1}{2} m_2 u_2^2 = \frac{1}{2} m_1 v_1^2 + \frac{1}{2} m_2 v_2^2$$
结合动量守恒方程，我们可以联立求解未知数 $v_1$ 和 $v_2$。经过数学化简，可以得到经典的一维完全弹性碰撞速度公式：
$$v_1 = \frac{(m_1 - m_2)u_1 + 2m_2 u_2}{m_1 + m_2}$$
$$v_2 = \frac{(m_2 - m_1)u_2 + 2m_1 u_1}{m_1 + m_2}$$
特殊情况分析
- 两物体质量相等 ($m_1 = m_2$)：
代入公式可得 $v_1 = u_2$，$v_2 = u_1$。这表明两个质量相等的物体发生完全弹性碰撞后，彼此交换了速度。例如，牛顿摆（Newton's cradle）就是这一原理的完美体现。 - 被撞物体静止 ($u_2 = 0$) 且质量极小 ($m_2 \ll m_1$)：
此时 $v_1 \approx u_1$，$v_2 \approx 2u_1$。重球撞击轻球时，轻球会以接近入射球两倍的速度弹开。
完全非弹性碰撞与能量损失
在完全非弹性碰撞中，$e = 0$，这意味着碰撞后两物体融为一体，具有共同的速度 $v$：
$$v = \frac{m_1 u_1 + m_2 u_2}{m_1 + m_2}$$
此时的动能损失 $\Delta E$ 可以通过碰撞前的总动能减去碰撞后的总动能来计算：
$$\Delta E = \left( \frac{1}{2} m_1 u_1^2 + \frac{1}{2} m_2 u_2^2 \right) - \frac{1}{2}(m_1 + m_2)v^2$$
将共同速度 $v$ 代入上式并化简，可得能量损失的便捷计算公式：
$$\Delta E = \frac{1}{2} \frac{m_1 m_2}{m_1 + m_2} (u_1 - u_2)^2$$
从公式中可以看出，能量损失 $\Delta E$ 永远大于或等于零，且只取决于两物体的相对速度平方以及折合质量。这部分能量通常转化为热能或形变能。
编程实践：模拟一维弹性碰撞
为了更好地理解上述物理过程，我们可以用 Python 编写一个简单的一维完全弹性碰撞模拟函数。
def elastic_collision(m1, u1, m2, u2):
"""
计算一维完全弹性碰撞后的末速度
:param m1: 物体1质量
:param u1: 物体1初速度
:param m2: 物体2质量
:param u2: 物体2初速度
:return: (v1, v2) 碰撞后的末速度
"""
v1 = ((m1 - m2) * u1 + 2 * m2 * u2) / (m1 + m2)
v2 = ((m2 - m1) * u2 + 2 * m1 * u1) / (m1 + m2)
return v1, v2
# 示例：质量为 2kg 的物体以 5m/s 撞击静止的 3kg 物体
m_1, u_1 = 2.0, 5.0
m_2, u_2 = 3.0, 0.0
v_1, v_2 = elastic_collision(m_1, u_1, m_2, u_2)
print(f"碰撞后物体1速度: {v_1:.2f} m/s")
print(f"碰撞后物体2速度: {v_2:.2f} m/s")
运行该代码，我们可以准确预测不同质量和初速度配置下的碰撞结果，这在游戏物理引擎开发和工程碰撞模拟中具有广泛的应用价值。