[论文解读] Gravitational potential and energy of homogeneous rectangular parallelepiped

Z. F. Seidov, P. I. Skvirsky|arXiv (Cornell University)|Feb 28, 2000

Dynamics and Control of Mechanical Systems参考文献 2被引用 7

一句话总结

本文通过符号积分，推导出均匀长方体（RP）的引力势和能量的精确解析表达式。利用Mathematica求解，得到初等函数形式的闭式解，可在体内外任意点实现精确计算，关键结果包括立方体中心和顶点处的势能，以及立方体、细杆和薄板的精确势能公式。

ABSTRACT

Gravitational potential and gravitational energy are presented in analytical form for homogeneous rectangular parallelepiped.

研究动机与目标

• 通过初等函数推导均匀长方体（RP）引力势的精确解析表达式，以闭式形式表示。
• 通过在体积分区内对势能进行积分，计算RP的引力势能，得到精确公式。
• 研究立方体、细杆和薄板等极限情况，以验证一般解并探讨其标度行为。
• 比较不同形状（包括椭球和RP）的WUM比值（势能 /（中心势能 × 质量）），揭示普遍趋势。

提出的方法

• 将牛顿引力势表述为对RP体积的三重积分，积分限固定。
• 使用符号计算（Mathematica）逐步求解积分：先对z积分，再对y积分，最后对x积分，得到对称的最终表达式。
• 采用函数 $ v(x',y',z') = y'z'\ln[x' + \sqrt{x'^2 + y'^2 + z'^2}] - \frac{x'^2}{2}\arctan\left(\frac{y'z'}{x'\sqrt{x'^2 + y'^2 + z'^2}}\right) $ 作为势能的核心组成部分。
• 通过在物体体积内对势能进行三重积分推导势能，利用对称性简化，并通过极限情况验证。
• 应用标度不变性，将结果表示为无量纲形式，以便在不同形状间进行普遍比较。
• 通过已知极限（线段、矩形、椭球）验证结果，并将WUM比值与椭球的已知值 $ 2/5 $ 进行比较。

实验结果

研究问题

• RQ1 在空间中任意点，均匀长方体的引力势的精确解析表达式是什么？
• RQ2 RP的引力势能如何依赖于其尺寸和密度？能否以闭式形式表达？
• RQ3 对于均匀立方体，其引力势在关键对称点（中心、顶点、面心）的值是多少？
• RQ4 势能与中心势能和质量乘积之比（WUM）在不同RP几何构型下如何变化？与椭球情况相比有何差异？
• RQ5 推导出的RP公式能否在极限情况下退化为已知的细杆和薄板结果？

主要发现

• 均匀长方体在中心处的引力势为 $ U(0,0,0) = 24\ln\left(\frac{1+\sqrt{3}}{\sqrt{2}}\right) - 2\pi \approx 9.52017 $，单位为 $ a^2 G \rho $。
• 立方体顶点处的势能恰好为其中心势能的一半：$ U(1,1,1) = \frac{1}{2} U(0,0,0) $。
• 边长为 $ 2a $ 的均匀立方体的势能为 $ W_C = 30.117 G \rho^2 a^5 $，由通解公式（31）导出。
• 对于细长棒（$ a \gg b = c $），势能按 $ W_{\text{stick}} = \frac{32}{3} G \rho^2 a b^4 \ln\left(\frac{a}{b}\right) $ 标度。
• 对于细长方形板（$ a \to 0 $，$ b = c $），势能为 $ W_S = 47.5714 G \rho^2 b^3 a^2 $，其WUM比值约为0.421673。
• WUM比值（势能 /（中心势能 × 质量））在立方体处达到约0.39544的最小值，趋近于细杆的0.5，趋近于薄板的约0.421673，表现出非单调行为，与椭球的 $ WUM = 2/5 = 0.4 $ 不同。