一、小振幅下的等时性推导 #

咱们直接从钟摆的运动微分方程入手。设摆长为$L$，摆球质量为$m$，摆角为$\theta$（以竖直向下为平衡位置，逆时针为正）。根据转动定律，力矩等于转动惯量乘以角加速度：

$$
I \ddot{\theta} = -\tau
$$

其中转动惯量$I = mL^2$，重力产生的回复力矩$\tau = mgL \sin\theta$（负号表示力矩与摆角方向相反）。代入后化简得到核心方程：

$$
\ddot{\theta} + \frac{g}{L} \sin\theta = 0
$$

这是个非线性微分方程，但当振幅$\theta_0 < 20^\circ$时，$\theta$很小，我们可以对$\sin\theta$做泰勒展开：

$$
\sin\theta = \theta - \frac{\theta^3}{6} + \frac{\theta^5}{120} - \dots
$$

小振幅下，$\theta^3$及更高阶项的影响可以忽略，近似为$\sin\theta \approx \theta$，此时方程退化为线性简谐振动方程：

$$
\ddot{\theta} + \omega^2 \theta = 0
$$

其中$\omega = \sqrt{\frac{g}{L}}$是角频率。这个方程的解为$\theta(t) = \theta_0 \cos(\omega t + \phi)$，对应的周期为：

$$
T = \frac{2\pi}{\omega} = 2\pi\sqrt{\frac{L}{g}}
$$

你看，这个周期表达式里完全没有出现振幅$\theta_0$——这就从数学上严格证明了小振幅下钟摆的等时性。

二、大振幅下的周期偏离（破坏等时性的数学根源） #

当振幅$\theta_0$增大时，$\sin\theta \approx \theta$的近似不再成立，我们必须求解原非线性方程。这里用能量守恒+椭圆积分的方法来推导周期：

1. 能量守恒建立积分式 #

摆球的机械能守恒：动能+势能=常数。当摆到最大振幅$\theta_0$时，速度为0，势能最大；平衡位置时势能最小，动能最大。据此列出等式：

$$
\frac{1}{2}mL^2 \dot{\theta}^2 + mgL(1 - \cos\theta) = mgL(1 - \cos\theta_0)
$$

整理后得到角速度与摆角的关系：

$$
\dot{\theta} = \sqrt{\frac{2g}{L} (\cos\theta - \cos\theta_0)}
$$

周期$T$是摆球完成一次全摆动的时间，也就是4倍的从$\theta=0$到$\theta=\theta_0$的单方向运动时间，因此：

$$
T = 4 \int_{0}^{\theta_0} \frac{d\theta}{\sqrt{\frac{2g}{L} (\cos\theta - \cos\theta_0)}}
$$

2. 变量替换为椭圆积分标准形式 #

利用三角恒等式$\cos\theta - \cos\theta_0 = 2\left[\sin^2\left(\frac{\theta_0}{2}\right) - \sin^2\left(\frac{\theta}{2}\right)\right]$，令$k = \sin\left(\frac{\theta_0}{2}\right)$，并做变量替换$\sin\left(\frac{\theta}{2}\right) = k \sin\phi$（$\phi$从0到$\frac{\pi}{2}$对应$\theta$从0到$\theta_0$）。

对$\theta$求导可得：

$$
\frac{d\theta}{d\phi} = \frac{2k \cos\phi}{\cos\left(\frac{\theta}{2}\right)} = \frac{2k \cos\phi}{\sqrt{1 - k^2 \sin^2\phi}}
$$

代入周期积分式，化简后得到：

$$
T = 4\sqrt{\frac{L}{g}} \int_{0}^{\frac{\pi}{2}} \frac{d\phi}{\sqrt{1 - k^2 \sin^2\phi}}
$$

这个积分就是第一类完全椭圆积分，记为$K(k)$，即$K(k) = \int_{0}^{\frac{\pi}{2}} \frac{d\phi}{\sqrt{1 - k^2 \sin^2\phi}}$。

3. 幂级数展开看周期与振幅的关系 #

将$K(k)$展开为幂级数（$k < 1$，对应$\theta_0 < 180^\circ$，这是摆能做完整摆动的前提）：

$$
K(k) = \frac{\pi}{2} \left[1 + \left(\frac{1}{2}\right)^2 k^2 + \left(\frac{1 \cdot 3}{2 \cdot 4}\right)^2 k^4 + \left(\frac{1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6}\right)^2 k^6 + \dots\right]
$$

把$k = \sin\left(\frac{\theta_0}{2}\right)$代入，当$\theta_0$很小时，$\sin\left(\frac{\theta_0}{2}\right) \approx \frac{\theta_0}{2}$（弧度制），所以$k^2 \approx \frac{\theta_0^2}{4}$，代入周期表达式：

$$
T = 2\pi\sqrt{\frac{L}{g}} \left[1 + \frac{1}{16}\theta_0^2 + \frac{11}{3072}\theta_0^4 + \dots\right]
$$

这里就能清晰看到大振幅对周期的影响：

• 当$\theta_0$很小时（比如$<20^{\circ$，约0.35弧度），$\theta_0}2 \approx 0.12$，$\frac{1}{16}\theta_0^2 \approx 0.0075$，周期仅比小振幅近似值大0.75%，几乎可以忽略，这就是我们说的“等时性”；
• 当$\theta_0$增大时，$\theta_0^2$及更高阶项的占比越来越大，周期$T$会随着振幅$\theta_0$的增大而单调增加——这就是大振幅破坏等时性的数学根源：非线性项$\sin\theta$的高阶项无法忽略，导致周期与振幅直接相关。

比如当$\theta_0 = 90^\circ$（$\pi/2$弧度）时，代入计算可得周期约为小振幅周期的1.18倍，偏差已经非常明显了。

内容的提问来源于stack exchange，提问作者user