[论文解读] Density of triangulated ternary disc packings

Thomas Fernique, Daria Pchelina|arXiv (Cornell University)|Nov 5, 2022

Mechanical Behavior of Composites被引用 1

一句话总结

本文研究了使用三种不同圆盘尺寸（三元圆盘堆积）在平面上的最高密度堆积，重点分析了每个孔洞均由三个相互相切的圆盘围成的三角剖分构型。通过借鉴开普勒猜想证明中使用的局部密度重分配方法，作者对164组半径对 (r, s)（其中 r > s > 0）进行了分类，识别出16组三角剖分堆积达到最大密度，15组两尺寸三角剖分堆积为最优，40组非三角剖分堆积更密集，另有93组因结构或计算限制导致该方法失效。

ABSTRACT

We consider ternary disc packings of the plane, i.e. the packings using discs of three different radii. Packings in which each "hole" is bounded by three pairwise tangent discs are called triangulated. Connelly conjectured that when such packings exist, one of them maximizes the proportion of the covered surface: this holds for unary and binary disc packings. For ternary packings, there are 164 pairs (r, s), 1 > r > s, allowing triangulated packings by discs of radii 1, r and s. In this paper, we enhance existing methods of dealing with maximal-density packings in order to study ternary triangulated packings. We prove that the conjecture holds for 31 triplets of disc radii and disprove it for 40 other triplets. Finally, we classify the remaining cases where our methods are not applicable. Our approach is based on the ideas present in the Hales' proof of the Kepler conjecture. Notably, our proof features local density redistribution based on computer search and interval arithmetic.

研究动机与目标

• 确定使用三种不同半径的三元圆盘堆积中，哪些能在平面上实现最高的可能密度。
• 评估对于给定的半径对，三角剖分堆积（每个孔洞由三个相切圆盘围成）是否能最大化密度。
• 将证明开普勒猜想时使用的局部密度重分配方法扩展至更复杂的三元圆盘堆积情形。
• 根据三角剖分、两尺寸或非三角剖分堆积是否实现最大密度，将164组可能的半径对 (r, s)（其中 1 > r > s > 0）进行分类。

提出的方法

• 将开普勒猜想证明中使用的局部化技术适配至圆盘中心的Delaunay三角剖分，实现密度的局部重分配。
• 定义顶点势能和星形约束，以建模局部密度贡献，确保在三角剖分构型中的一致性。
• 使用线性规划与区间分析，验证在由几何与拓扑条件导出的约束下，势能分布的可行性。
• 构建一个在 R^9 中的9维多面体，以表示紧致顶点势能与质量参数 (m1, mr, ms) 的有效组合，通过检查其非空性来验证方法的有效性。
• 应用系统性验证流程，识别违反局部密度不等式的异常三角形，通过迭代添加约束条件来优化势能解。
• 根据多面体是否为空、是否存在解，或因结构复杂性或计算限制导致方法失效，对各类情况进行分类。

实验结果

研究问题

• RQ1 在164组半径对 (r, s)（其中 1 > r > s > 0）中，有多少组的三元圆盘堆积最大密度由三角剖分构型实现？
• RQ2 在哪些情况下，最大密度并非由三角剖分三元圆盘堆积实现，而是由仅使用两种圆盘尺寸的三角剖分二元堆积实现？
• RQ3 对于哪些半径对，非三角剖分堆积的密度严格高于任何三角剖分堆积？
• RQ4 为何标准密度重分配技术在某些半径对上失效？结构或计算上的何种障碍导致无法得出明确分类？
• RQ5 该方法能否扩展至处理具有多个星形层或非唯一局部构型的情形？

主要发现

• 对于16组半径对 (r, s)，最大密度由三角剖分三元圆盘堆积实现。
• 另有15组半径对，最大密度由三角剖分二元圆盘堆积（仅使用三种尺寸中的两种）实现，而非完整的三元构型。
• 对于40组半径对，发现非三角剖分堆积的密度严格高于任何三角剖分堆积，表明在这些情况下三角剖分并非最优。
• 在52组情况下，可行势能值的多面体被判定为空，导致方法因几何或代数不相容而无法适用。
• 四组情况（45, 62, 75, 144）仍未解决：多面体非空，但不存在任何有效的势能构型能同时满足所有三角形的局部密度不等式，表明需要采用更非局部的方法。
• 作者识别出22组情况，在饱和三角剖分堆积中，圆盘存在多个可能的星形层，这使得标准势能方法失效，需采用更复杂的重分配策略。