一维势能景观

只看势能曲线，能预测粒子怎么运动吗？ 大部分能。在势能图上，按总能量画一条水平线：曲线低于这条线的地方，粒子才去得了，两者相交处就是转折点。谷底弯得多厉害，只决定小振动频率。谷底要是像四次方曲线那样平，弯曲为零，算出的小振动频率也是零，可粒子照样被困在里面，是稳定的。

改变势阱形状，拖动质点后释放。观察力和速度怎样变化，再从能量图与相图读出转折点、允许区和周期。

质点位置 $0$

速度 $0$

机械能 $0$

完整记录的振动周期 $0$

选中时刻 $0$

保守力 $0$

黏性阻力 $0$

势能 $0$

动能 $0$

耗散能量 $0$

局部小振动角频率 $0$

选定极小值 $0$

初始能量左边界 $0$

初始能量右边界 $0$

保守能量积分周期 $0$

全记录最大能量平衡误差 $0$

半步长位置差 $0$

完整记录的鞍点穿越次数 $0$

周期积分分辨率差 $0$

采用 SI 单位，向右为正。质点在给定的多项式力场中沿水平约束运动；导轨与球体是示意，不模拟滚动或力场执行机构。能量曲线是图表，不是真实轨道。灰色相曲线按释放能量计算，青色为数值轨迹。每次试验共用一个长度比例，力和速度各用固定比例。每步检查机械能与耗散热的总账；半步长位置比较最多覆盖前十秒。

能量图预测整条轨迹

背景势能图告诉你质点在每个位置有多少势能。曲线的斜率给出力，平坦的点就是平衡位置。在总能量处画一条水平线，曲线低于这条线的位置，就是经典质点能到的地方。这条曲线不是轨道：质点沿横轴运动，纵轴表示能量。

为什么值得理解OpenStax 大学物理第 8.4 节用势能图判断平衡稳不稳定、质点在哪里转折。这个实验更进一步：可以改势阱形状，在相空间里看运动，用积分算周期，还能看数值积分的能量漂了多少。

核心问题: 势阱底部不再是抛物线形状时，哪些预测还成立？
一句话直觉: 质点在哪里转折，要用完整的能量方程算；小振动频率只看平衡点附近弯得多厉害。

核心数学模型

可以改形状的势阱

U(x)=\frac{kx^2}{2}+\frac{qx^4}{4}+fx

二次和四次系数都不取负值，只要有一个为正，势阱就能困住质点。倾斜项把平衡点推向一侧。

力与稳定性

\begin{aligned}F(x)&=-U\prime(x),\\ U\prime(x_*)&=0,\\ \omega_0^2&=\frac{U\prime\prime(x_*)}{m}\end{aligned}

极小值处曲率为正，小振动就是简谐振动；曲率为零时，要看更高阶的项才能判断稳不稳定。

允许区域与转折点

\begin{aligned}E&=\frac12mv^2+U(x),\\ U(x)&\le E,\\ U(x_\pm)&=E\end{aligned}

质点在转折点速度为零。从平衡点静止释放时，允许区域缩成一个点。

有限振幅下的周期

\begin{aligned}T(E)&=\sqrt{2m}\\ &\quad\times\int_{x_-}^{x_+}\frac{dx}{\sqrt{E-U(x)}}\end{aligned}

这个积分独立于模拟单独计算。被积函数在转折点发散，余弦代换能去掉这种简单奇点；读数还显示分辨率减半后结果差多少。

能量漂移与小振动参考

\begin{aligned}\delta E(t)&=E(t)-E(0),\\ x_{\mathrm{lin}}(t)&=x_*\\ &\quad+(x_0-x_*)\cos\omega_0t\\ &\quad+\frac{v_0}{\omega_0}\sin\omega_0t\end{aligned}

参考曲线只在平衡点附近成立。小振动频率为零时，参考变成匀速运动。这只是线性近似，并不说明真实的极小值是中性稳定的。

常见难点

能量不是高度

容易误解质点会顺着曲线往下滑。

正确理解纵轴是能量。质点实际上沿横轴来回运动。

力取决于斜率

容易误解势能最大处的力最大。

正确理解力看的是该点的斜率，不是高度。势阱倾斜后，平衡点处的势能甚至可以不为零。

曲率为零，要看高阶项

容易误解小振动频率为零就说明不稳定。

正确理解正的四次项照样能困住质点。质点仍然来回振动，周期有限，并随振幅变化。

动手实验

01
先检验简谐势阱
先预测简谐周期，再选“简谐运动”，展开“周期、受力与数值检查”，比较两项周期读数。
你会看到： 在两个转折点之间积分得到的周期，与简谐振动的解析周期一致，也与轨迹上反复转折的间隔一致。
02
改变势阱形状
先预测释放点处力的方向，再改变“二次刚度”和“四次束缚”，展开数值检查，对照实际受力与局部线性预测。
你会看到： 释放能量和转折点一起变化，控件上显示的是实际系数。
03
让极小值变平
选择“平坦阱底”，先预测释放位移减半后周期如何变化，再把质点拖近极小值。
你会看到： 小振动频率一直是零，实际周期却随释放幅度变小而变长。
04
移动平衡点
选择“移动极小值”，在相图中横向和纵向拖动圆点，比较改动前后的能量边界。
你会看到： 平衡点移动了。改变初速度，允许区间跟着变，势能曲线本身不变。