磁性摆锤混沌模拟器

摆锤的前端（质量m）受到重力、空气阻力和3块磁石的类库仑力。运动在二维平面内发生，遵循牛顿运动方程。

$$ m\ddot{\mathbf{r}}= -mg\hat{\mathbf{y}}- \gamma \dot{\mathbf{r}}+ \sum_{i=1}^{3}\frac{k(\mathbf{r}- \mathbf{r}_i)}{(||\mathbf{r}- \mathbf{r}_i||^2 + \epsilon^2)^{3/2}}$$

其中，$\mathbf{r}$ 是摆锤前端的位置矢量，$\mathbf{r}_i$ 是磁石i的位置，$g$ 是重力加速度，$\gamma$ 是阻尼系数（模拟器参数），$k$ 是磁力强度，$\epsilon$ 是正则化参数，防止计算发散。

这个二阶微分方程可以转化为关于位置和速度的一阶联立微分方程组，用RK4法数值求解。

$$ \frac{d}{dt}\begin{bmatrix}\mathbf{r}\\ \dot{\mathbf{r}}\end{bmatrix}= \begin{bmatrix}\dot{\mathbf{r}}\\ \frac{1}{m}\left( -mg\hat{\mathbf{y}}- \gamma \dot{\mathbf{r}}+ \sum_{i}\mathbf{F}_{mag,i}\right) \end{bmatrix} $$

这种形式叫"状态空间表示"，是计算机数值模拟的标准形式。右边力的计算是非线性的，这就是混沌行为的根源。

非线性动力学与混沌理论教学：初值敏感性、分形吸引盆等抽象概念能通过视觉直观理解，是最优的演示工具，常用于实验室教学。

复杂系统建模：气象预报、股市变动、生态个体数等许多非线性相互作用的系统，可以用这个简化模型理解其基本行为。

CAE非线性动力分析：汽车悬架、飞机着陆装置等具有多个稳定状态的机械系统，其动力特性模拟采用类似原理。对初值敏感的系统设计需特别注意。

艺术与设计：分形吸引盆生成的无穷复杂优美图案，启发了数字艺术和平面设计创意。

用这个模拟器时有几点要注意。首先是一个常见误解："初始位置改变一点点结果就大不同，所以计算必须超级精确"。事实是，初值敏感性确实是混沌的核心，但模拟本身的数值误差也会被当作"换了初值"。所以选择积分方法（这里用RK4）和时间步长(Δt)特别重要。比如把Δt从0.01秒改成0.1秒，轨迹会完全变样——这是"数值误差造成的混沌"而非真正的"决定论混沌"。CAE实务中也常出现这种误区：网格太粗、时间步太大，反而产生不符合物理的振动。

其次，参数搭配要平衡。阻尼系数γ接近0时能量守恒，轨迹复杂漂亮，但现实系统总有空气阻力，必然存在衰减。γ太大摆锤很快停止，混沌特性就消失了。比如质量m=1、磁力k=1000的系统，γ=0.1～0.5通常是混沌与衰减的好平衡点。CAE中阻尼系数设太离谱，模拟结果也会和试验数据不符。

最后，别觉得这个模拟"完全真实"。这里的磁力模型是简化版，没完全复现真实磁偶极子间的相互作用。参数ε是为了数值稳定性"加的辅助"。这种思路和CAE里用线性弹簧简化复杂现象是一样的：首先提取现象本质，而不是追求完美重现。