参数
操作
预设方案
显示设置
显示网格
粒子轨迹
力向量
统计
0
引力源数
0
粒子数
0
总动能
60
FPS
点击画布放置引力源，体验时空橡皮膜类比下网格的弯曲变形。观察测试粒子描绘出的轨道，直观感受广义相对论的物理图像。
两个质量 M 与 m 之间的引力：
本模拟器在每个时间步计算加速度 a = G·M/r²，并通过Verlet积分更新粒子的位置和速度。引入软化参数ε以防止r→0时的数值发散。
网格位移量：
爱因斯坦用"重物压弯橡皮膜"来比喻广义相对论中质量弯曲时空的概念。本模拟器的二维网格变形正是这一类比的直观呈现。真实时空是四维的，但这种可视化对于直觉理解非常有效。
逃逸速度：
粒子初速度低于逃逸速度时做椭圆轨道运动，高于时做双曲线轨道飞离。在CAE领域，轨道计算技术被应用于火箭轨道设计和人造卫星摄动分析。
引力势满足泊松方程 ∇²φ = 4πGρ，在数学上与FEM中热传导、静电场、结构势能等问题同构。网格变形可视化与网格变形分析及拓扑优化采用相同的概念。N体数值积分也是MD（分子动力学）模拟的基础。
模拟的核心是牛顿万有引力定律，计算测试粒子受到的引力加速度：
$$ \vec{a}= - \sum_i \frac{G M_i}{(\|\vec{r}- \vec{r}_i\| + \epsilon)^2}\cdot \frac{\vec{r}- \vec{r}_i}{\|\vec{r}- \vec{r}_i\|}$$这里，$\vec{a}$ 是测试粒子的加速度，$G$ 是引力常数，$M_i$ 是第 $i$ 个引力源的质量，$\vec{r}$ 和 $\vec{r}_i$ 分别是粒子和引力源的位置矢量，$\epsilon$ 就是软化参数，用于防止距离 $\|\vec{r}- \vec{r}_i\|$ 趋近于零时的数值发散。
为了计算粒子轨迹，模拟器采用了Verlet数值积分方法，这是一种在分子动力学和游戏物理中常用的稳定算法：
$$ \vec{r}(t + \Delta t) = 2\vec{r}(t) - \vec{r}(t - \Delta t) + \vec{a}(t) \Delta t^2 $$其中 $\vec{r}(t)$ 是当前时刻的位置，$\vec{r}(t - \Delta t)$ 是上一时刻的位置，$\Delta t$ 是时间步长。这个公式直接用加速度更新位置，无需显式计算速度，具有良好的能量守恒特性。
航天器轨道设计与导航：在规划月球探测器或火星车的轨道时，工程师需要精确计算天体（地球、月球、太阳）引力造成的复杂时空弯曲效应。类似的引力计算和轨道积分是任务成功的基石。
引力透镜效应研究：在天文学中，遥远星系发出的光在经过大质量天体（如星系团）附近时，其路径会因为时空弯曲而发生偏折，形成多个像或光环。模拟质量分布对光路的影响，是分析引力透镜现象的关键。
CAE软件中的多体动力学分析：在机械系统仿真中，当分析卫星太阳翼的展开、空间机械臂的运动时，需要考虑部件间微小的引力相互作用及其长期累积效应，其控制方程与模拟器中的粒子运动类似。
游戏与科幻视觉特效：为了在电影或游戏中真实渲染黑洞等极端天体周围的景象，特效团队需要基于广义相对论计算光线在弯曲时空中的路径，模拟器中的“引力井”概念是理解这一视觉扭曲的直观起点。
首先需要明确，本模拟器是一个二维类比。实际的时空弯曲发生在三维空间加时间的四维维度中，这个“橡胶膜”意象仅作为理解辅助。尤其要把握本质：粒子“下落”并非朝向膜的“下方”，而是在扭曲几何结构中行进的结果。
其次，关于“软化”参数设置。若将其设为0.1或0.01等较小值，粒子易与引力源剧烈碰撞导致计算发散。实际应用中需根据模拟尺度调整：例如星系碰撞模拟中恒星间距极大，通常取较大软化值以保持计算稳定；反之对小尺度系统进行精密计算时需减小该值，但需同步缩小时间步长。本模拟器中“粒子初速度”较大时轨道失稳，正是因固定时间步长导致的现象，这也是实际工程中需要注意的要点。
最后，“Verlet积分法”的特性。该方法虽具有良好能量守恒性，但速度不直接计算（通过位置差分求得），因此后续添加速度相关力（如空气阻力）时需要特别处理。本模拟器仅涉及纯重力场故能直接使用，但在扩展应用时需将此特性纳入考量。
本模拟器核心的“质点在引力场中的运动计算”，其实与分子动力学（MD）的数值计算方法高度相似。MD通过计算原子间作用力（如Lennard-Jones势），同样采用Verlet积分追踪原子轨迹。尽管引力遵循 $1/r^2$ 规律而分子间力遵循 $1/r^{12}$ 等不同法则，但防止计算发散的截断法与软化思路是相通的。
另一重要应用在于机械臂与操作器的轨迹规划。若将引力势阱视为“目标位置势场”，粒子视为“机器人末端执行器”，并设置障碍物排斥势场，即可通过相同力学计算生成机器人自然趋近目标且规避障碍的路径。这正是一种经典而有效的人工势场法。
此外，流体模拟中的粒子法（如SPH法）在计算邻近粒子相互作用时，也必须采用类似“软化”处理的平滑化技术（核函数）。这种通过力场平滑化防止计算发散的工程基础技术，已跨越学科边界得到广泛应用。
第一步应夯实数值计算与经典力学基础。理解本模拟器的Verlet积分法时，可对比欧拉法、龙格-库塔法等其他数值积分方法。例如在相同初始条件下使用欧拉法，能直观体验能量不守恒导致的轨道膨胀等“发散”现象。推荐《欧拉法进阶：数值计算基础》这类入门教材。
数学背景方面，微分方程（特别是常微分方程初值问题）与向量分析至关重要。粒子运动本质是牛顿运动方程 $m\ddot{\vec{r}} = \vec{F}$ 这类二阶常微分方程，计算机模拟正是对其求解。同时，学习将力 $\vec{F}$ 表示为势能 $U$ 的梯度（$-\nabla U$）的方法，能深化对守恒力的理解——如同球体沿势能谷（引力势阱）滚落的意象。
熟悉本工具后，建议探索“多体问题”与“相对论效应”。当前模拟接近引力源固定的“限制性三体问题”，若多个引力源动态相互作用（如双星系统行星轨道），则会出现混沌轨迹。而要理解水星近日点进动等牛顿力学无法解释的细微现象，则需要进阶到后牛顿近似等广义相对论领域。通过本模拟器反复调整参数培养的“物理直觉”，将成为探索这些高阶课题的基石。