One String to Pull Them All: Fast Assembly of Curved Structures from Flat Auxetic Linkages

问题与目标

本文研究可展开结构（deployable structures）的设计与制造：如何把一个自由曲面目标形状，分解为可平铺制造的模块，并且能够以极简的方式从平面状态展开成三维曲面。已有可展开结构方法在可逆性、可扩展性、以及对壳状双曲率几何的可靠支撑上都存在缺陷——有的方法无法逼近平均曲率变号（含负平均曲率）的几何，有的需要专用设备驱动、或存在卡死（jamming）风险，还有的需要按精确顺序手工组装大量离散片。

作者提出一个计算方法：把目标曲面用刚性四边形棱台瓷砖（quadrilateral frustum tiles）逼近，构成一个平面机械联动机构（flat mechanical linkage），瓷砖之间在顶点处用转动铰链连接。仅通过一根绳子的单次拉动（single string pull），平面联动机构就能平滑地组装成目标三维结构；松开绳子则回到平面，实现可逆驱动。整个机构只需单向拉动，因此与尺度无关。

核心贡献

• 把用户给定曲面模块化为尺度渐变的四边形棱台，并计算尊重几何曲率的铰链位置，得到平面机械联动机构；
• 给出一个一般性分析（定理）：只要闭合边界处的所有间隙，就足以将全部联动瓷砖引导到目标三维形态；
• 识别组装所需的最小起吊点集合（minimum set of lift points），并优化绳路以最小化通道摩擦，从而降低平面到组装过程中的总阻力；
• 支持大尺度范围上的可逆驱动，用人手、机械臂、建筑吊车驱动的多个原型验证。

方法

整体流水线：给定自由曲面 → 优化编码几何的四边形瓷砖布局 → 平铺瓷砖为无重叠布局并计算铰链 → 选择最小起吊点集合 → 计算低摩擦闭合绳路 → 制造平面联动机构并单次拉绳组装。

表面有理化（Surface Rationalization）

单位瓷砖是由八个顶点定义的四边形棱台，顶底面平行，含 \(h\) 厚度与相邻四边形间共面上的二面角 \(\theta\)。四片瓷砖围成一个单元胞，中间的空隙称为 gap。

初始化借鉴之前的 auxetic linkage 工作，用共形映射建立三维联动机构与二维平面之间的对应。从目标曲面 \(S\) 出发，用 Boundary First Flattening 展平到平面域 \(\Omega \subset \mathbb{R}^2\)，得到共形映射 \(c: S \to \Omega\)。在 \(\Omega\) 上覆盖规则四边形网格并裁掉边界外的四边形，得到平面顶面网格 \(M_{2D}\)；再用逆映射 \(c^{-1}\) 抬升得到三维顶面网格 \(M_{3D}\)。当共形尺度因子（CSF）超过 2 时，说明四边形 auxetic 图案无法局部扩张到足够逼近目标，需要引入分割（split）。

组装配置优化：共形抬升不保持四边形面的平面性，需要在保持贴近目标曲面的前提下优化。基于投影法（projection-based）最小化组合能量：

\[E_{Assembled}(v) = \omega_1 E_{Planar}(v) + \omega_2 E_{Square}(v) + \omega_3 E_{Surface}(v)\]

其中平面性能量为

\[E_{Planar} = \sum_{q \in Q(M_{3D})} \| V_q - P_P(V_q) \|_F^2\]

将每个四边形投影到其最佳拟合平面。表面贴合能量把顶点拉向目标曲面上的最近点 \(P_S\)；方正性能量由边长项 \(E_{Length}\) 和形状项 \(E_{Shape}\) 组成，用于缩小四边形尺寸差异（避免 auxetic 联动机构开合阻力过大）。目标边长取相邻两四边形平均边长的均值：\(L_{Target}(e_{ij}) = \frac{\bar{L}_a + \bar{L}_b}{2}\)。优化后挤出（extrude）得到带厚度的三维瓷砖网格 \(T_{3D}\)。

平面配置优化：从裁剪网格 \(M_{2D}\) 出发，优化其顶点以匹配 \(K_{3D}\) 的边长，得平面联动机构 \(K_{2D}\)，再挤出成平面瓷砖 \(T_{2D}\)。目标能量：

\[E_{Flat}(v) = \omega_1 E_{Edge}(v) + \omega_2 E_{Collision}(v) + \omega_3 E_{Fab}(v)\]

分别度量边长错配、平面配置中四边形重叠（沿用非穿透约束），以及可制造间隙角约束。用户可设最小间隙角 \(\theta_{min}\)，最大角固定为 \(90^\circ\)。

铰链放置策略

铰链要放在每两个相连瓷砖之间。当挤出使二面角 \(\theta > 90^\circ\) 时，若所有铰链都放在顶面，会在平面布局产生物理无效的重叠区。作者给出两种方案：其一是铰链一律放顶面，把重叠区从其中一片瓷砖上切掉，制造简单但削弱受压承载能力与内部美观；其二是根据关节处沿开合二面角方向的法曲率符号，决定铰链放顶面还是底面（曲率为负则二面角超过 \(90^\circ\)，铰链放底面），保证侧面完全接触。第二种方案会引发相邻瓷砖的连锁重叠，需要额外的全局重定位优化：

\[E_{Hinge}(v) = \omega_1 E_{Rigid}(v) + \omega_2 E_{Collision}(v) + \omega_3 E_{Conn}(v)\]

其中刚性约束按 Maxwell 框架对每片瓷砖施加 18 条边约束（4 顶边、4 底边、4 侧边、6 面对角线）；连接能量

\[E_{Conn} = \sum_{u \in H(T_{2D})} \| x_u - P_{Conn}(x_u) \|_2^2\]

惩罚铰链顶点对的分离，使成对铰链重合形成共享枢轴。

网格分割（Mesh Splitting）

有界共形尺度因子 \(\lambda\) 约束在 1 与常数 \(\sigma > 1\) 之间，\(\sigma\) 取决于联动拓扑：三角形 auxetic 最大面积扩张 \(\sigma = 4\)，四边形 auxetic 为 \(\sigma = 2\)。当尺度因子大于 2 时引入分割：沿网格方向、穿过最高高斯曲率点，层次化地二分网格，直到每部分满足 \(\sigma \le 2\)。完整分割保持联动机构内规则连接，不引入奇异点或不规则顶点，并与 auxetic 网格的栅格平行，便于穿绳且保持相同驱动机制。

绳路布置与驱动

边界约束（关键观察）：一旦边界瓷砖就位到目标配置，整个结构的几何就被唯一确定。据此先把绳子沿边界瓷砖穿过。定理表述为：由八顶点四边形棱台刚性瓷砖构成、每片在单一角点与邻居相连、且每个内部空隙为四顶点菱形的联动机构，在被一根连续边界绳围绕并张紧时是等静定的（isostatic），只允许六个平凡刚体运动，无内部无穷小挠曲。因此边界瓷砖收缩到目标配置时，内部联动逐步闭合并在所有间隙闭合后锁定。

最小起吊点选择：把每个 gap 建模为离散能量场，用重力势能（GPE）量化承载重要性：

\[g_i = \sum_{t_k \in T(G_i)} \frac{1}{4} m_k g (z_k - z_{min})\]

用离散 Morse–Smale 分割识别真正的能量峰及其盆地，构建 gap 之间的无向图（共享瓷砖则连边），边权取 \(w(u,v) = \min\{g_u, g_v\}\)，用 Kruskal 算法求最大生成树（MST）。两峰间唯一路径上的最弱边定义能垒：

\[barrier(i,j) = \min_{(u,v) \in Path_T(i,j)} w(u,v)\]

耦合系数 \(c(i,j) = \frac{barrier(i,j)}{g_j}\)（\(g_i \ge g_j\)）：\(c \approx 1\) 表示强耦合，\(c \approx 0\) 表示两峰几乎独立。据此构建有向无环图（DAG），用耦合阈值 \(\tau\) 判定是否耦合，连通分量即强耦合峰簇，每簇最高能量峰（入度为零节点）作为该簇起吊点，簇数即最小起吊点数。用阈值扫描策略（起始 \(\tau_0 = 0.8\)，失败则每次加 0.1 重建并重跑仿真）缓解阈值敏感性。该方法还能自动融合低分辨率网格造成的虚假重复峰。

最小摩擦绳路：用 Capstan 方程建模通道能量，绳以张力 \(T_1\) 进入、累计缠绕角 \(\theta_{Total}\)、通道壁摩擦系数 \(\mu_c\)：

\[E_{Channel} = T_1 (e^{\mu_c \theta_{Total}} - 1)\]

把平面配置表示为无向图 \(G = (V, E)\)，每个节点是一个 gap，共享顶点则连边，形成栅格状结构。每个 gap 节点带整数标签：0 与 1 分别表示竖直/水平 gap，\(-1\) 表示允许进入内部起吊点 gap 的虚拟边界入口，\(-2\) 表示网格分割产生的分割边界 gap。目标是找一条有效闭合走路 \(R\)，最小化累计转角 \(\min_R \theta_{Total}\)（转角在平面几何位置上按前驱-自身-后继质心夹角度量），并满足：最小化交叉（每个 gap 只被一条绳穿过，从边界进入内部时应从 \(x_b = -1\) 的 gap 进入），以及连接分割（可经过 \(x_b = -2\) 的边界节点但不能由其进入内部，以保持各分割部件连通）。

仿真与实现：用 Projective Dynamics 从平面到组装做仿真评估。将绳的驱动抽象为两种几何动作——把相邻侧面中点重合以闭合 gap 的 snap，以及在出口处把对应 gap 抬到指定高度的 lift，编码为位置约束以加速仿真。仿真能量：

\[E_{Simulation}(v) = \omega_1 E_{Rigid}(v) + \omega_2 E_{Collision}(v) + \omega_3 E_{Actuation}(v)\]

工具基于开源库 ShapeOp 与 libigl。约 10k 三角形网格：组装配置优化约 1.2 分钟，平面配置优化约 0.6 分钟（Mac M1 Pro，八核 3.2 GHz），驱动仿真额外约 5.5–6 分钟。给出各能量权重供复现，例如 \(E_{Assembled}\) 取 \(\omega_1 = 10000, \omega_2 = 10, \omega_3 = 0.1\)（高平面性权重保证可制造，低表面匹配权重允许适度偏离目标）。

制造与结果

在多种双曲率曲面上评估（如 ThreeShell、Seashell、Peanut、Saddle、Cloud、Lilium），并给出物理制造原型。

制造方式：auxetic 机构包含刚性瓷砖与柔性铰链两部分。给出两种技术：其一是双材料 3D 打印，柔性 TPU 做铰链、PLA 做瓷砖（TPU 关节 0.6mm 厚提供足够柔性与抗撕裂、0.8mm 宽避免过刚导致塑性变形，TPU 与 PLA 间用互锁结构增强粘合）；其二是激光切割 EPDM 橡胶铰链粘到 PLA 瓷砖上（旋转阻力更低但需额外手工粘合）。体现质量与制造便利之间的权衡。

铰链放置消融：顶面铰链会使凹特征区域下垂（因平面布局中切掉了重叠），双铰链方法在这些区域自动把铰链放底面，实现瓷砖侧面完全接触、防止下垂。

小尺度原型：用 Ultimaker S5 双材料 3D 打印，组装配置平均尺寸约 \(164 \times 121 \times 56\) mm、平均质量约 135 g。含分割的结构在分割区侧面嵌入磁体以精确对齐。

大尺度原型：制造了一把真实尺寸、可用的椅子，组装后约 3 英尺高，用激光切割胶合板箱体、橡胶条做铰链、涂层不锈钢缆穿过瓷砖内胶合的塑料通道。椅子能平滑展开、松张后保持稳定、可再压平折叠成紧凑形态。

平面性误差：以顶点到最佳拟合平面最大距离量化。3D 打印 Lilium（\(211 \times 212 \times 169\) mm，非锥形网格）经优化，顶面从 0.0074mm 降到 0.00036mm、底面从 0.1729mm 降到 0.00038mm、接触面从 0.587mm 降到 0.00039mm，通常降低 100 到 1000 倍。

为何先优化组装再优化平面：共形映射配有界尺度因子保证平面嵌入始终存在，平面配置边长总能匹配已优化好的组装配置，故不必联合优化。方法优先保证四边形 gap 而非精确曲线对齐，得到便于穿绳与组装的栅格状图案；用共形映射初始化以避免五价奇异点（quad 联动始终需要四顶点）。

应用案例

• 消费品：把椅子的人体工学曲线编码进模块化平面态，平面制造、平铺运输，用户单次拉绳即完成组装，也可折叠紧凑存放；
• 个性化医疗器械：夹板与姿态矫正器，按患者定制、平放存储、需要时拉绳组装；
• 快速机器人组装：机器人只需检测并拉动一根绳，用人形机器人组装含 193 个模块的 Lilium 原型，比逐块检测与运动规划简单得多；
• 空间探索与建筑结构：平面制造的模块化太空栖息舱可紧凑存储运输、机器人快速展开/拆卸/重组；建筑结构可做成中空框架，玻璃面板装在顶面，用少量吊车快速展开。

局限与未来工作

• 当前实现聚焦盘状拓扑（disk topology）输入曲面，推广到其他拓扑（如球状）主要受共形映射难以始终对齐四边形图案的限制；
• 非常复杂的几何（如 Stanford Bunny）需大量分割才能把尺度因子压到阈值以下，如何最优放置大量分割是独立的复杂研究问题；
• 理论保证边界约束下组装配置刚性，实测所有原型松张后仍能自支撑，但未经严格测试或在优化中强制，可作为单绳驱动自支撑结构的方向；
• 建筑级尺度部署的硬件规格（缆径、铰链强度）等工程难题留待未来；
• 当前仿真未纳入风载、热应力等环境载荷。

结论

方法把 auxetic 联动机构与驱动绳的交互仅通过瓷砖连接性抽象出来，使其对尺度、材料、制造技术不敏感；流水线统一了几何驱动的四边形铺砌、无碰撞展平、起吊点最小化与摩擦感知绳路布置，产出可展开、刚性、可逆的结构，无需专用设备，建立了一个高效、可扩展、可逆的组装范式。