单摆模拟器

单摆模拟器的物理模型中，质量 \( m \) 的质点用长度 \( L \) 的轻绳悬吊，在重力加速度 \( g \) 下运动。对于摆角 \( \theta \)，考虑阻尼系数 \( b \)，运动方程为 \( \frac{d^2\theta}{dt^2} = -\frac{g}{L} \sin\theta - b\frac{d\theta}{dt} \)。本模拟器不采用线性近似 \( \sin\theta \approx \theta \)，而是直接保留非线性项，从而能够重现振幅相关的周期变化和混沌行为。采用4阶龙格-库塔法（RK4）进行数值积分，足够小的时间步长 \( \Delta t \) 确保高精度遵守能量守恒定律。引入角速度 \( \omega = d\theta/dt \)，状态方程变为 \( \frac{d\theta}{dt} = \omega \)、\( \frac{d\omega}{dt} = -\frac{g}{L}\sin\theta - b\omega \)，RK4依次迭代更新。通过"摆动动画"、"角变位波形"、"相空间（\( \theta-\omega \) 平面）"三个可视化窗口，可直观观察阻尼和初始条件对非线性振动的影响。

长度为 $L$ 的单摆在振幅较小时可视为简谐振动，其周期由下式给出。

$T = 2\pi\sqrt{\dfrac{L}{g}}$

值得注意的是，周期既不依赖振幅也不依赖质量，仅由长度 $L$ 与重力加速度 $g$ 决定（等时性）。这正是摆钟的原理。由于周期与长度的平方根成正比，将 $L$ 增大到4倍时周期变为2倍。反之，可由周期与长度测定 $g$。

大振幅修正：等时性是在微小振幅（$\sin\theta\approx\theta$）下的近似。当振幅 $\theta_0$ 增大时周期略微变长，修正式为 $T \approx 2\pi\sqrt{\dfrac{L}{g}}\left(1+\dfrac{\theta_0^2}{16}+\cdots\right)$。

阻尼：存在空气阻力或轴摩擦时，振幅逐渐衰减（阻尼振动）。在相空间（角度-角速度平面）中，无阻尼的摆为闭合轨道，有阻尼时则向内卷曲成螺旋，可直观把握运动情形。在本模拟器中可改变振幅、长度与阻尼，观察周期与轨道的变化。

工业实际案例
汽车悬挂设计中，减振器原理类似本模拟器的阻尼模型。丰田、日产等车企通过非线性振动分析优化减振参数，在乘坐舒适和操纵稳定间取得平衡。建筑业则在超高层建筑（如"阿倍野哈鲁卡斯"）制振装置中应用摆的原理，调整阻尼系数和周期以抑制风致摇晃。

科研·教育应用
大学物理实验和机械工程课程中，用本模拟器直观学习非线性现象（混沌、能量耗散）。相比纯理论公式，学生通过调参自主探索"为什么振幅衰减"、"位相平面为何螺旋收缩"等问题，学习效果显著。东京大学、东北大学等顶级院校已引入类似工具。

CAE分析的前置步骤
本模拟器是正式CAE软件（如ANSYS、Abaqus）前的概念验证阶段。实务中先用本工具确定摆的基本参数（固有频率、阻尼比），再作为边界条件输入详细3D有限元分析。例如机器人关节振动抑制，先在本模拟器调整阻尼系数，再对实机CAE模型进行扭矩优化，可降低试原型20～30%的成本。

"单摆运动总是正弦波形"的想法是不对的。初始角度超过30°时非线性效应明显，周期变长，波形失真。本模拟器通过RK4求解非线性微分方程，能实时展示这种偏离。

"阻尼系数越大越快停止"听起来合理，但实际上过度阻尼会导致无振荡返回，即振幅快速衰减但不再摇晃。从相空间（θ-ω）的轨迹看，临界阻尼和过阻尼状态的差异一目了然。

"重力加速度变化只是让运动快慢改变"的理解太简化了。重力是复原力的来源，小g区域的非线性程度相对增强。月球低重力下，阻尼影响相对更明显。