Developable Quad Meshes and Contact Element Nets

一句话总结

本文提出一种作用于带顶点权重的四边形网格的局部可展性判据：给每个面配一个”接触元”（一个带权点加法向量），给每条边配一个”直纹向量”，由此得到”每个面一条可展性方程”，并可用标准优化求解，从而实现可展曲面的交互设计与放样（lofting）。

研究背景

• 领域现状：可展曲面（把平整、不可拉伸的板材弯折而成的形状）在建筑与制造中价值很高，几何建模研究已持续数十年。可展性有多种等价刻画——高斯曲率为零、存在到平面的等距映射、沿直纹切平面恒定、第二基本形式降秩等——每一种都催生了一类离散建模方法（基于样条、基于直纹四边形网格、基于正交测地网、基于三角网格局部铰链、基于高度场 Hessian 秩最小化等）。
• 核心痛点：以往方法大多要求网格的边对齐到曲面上的特殊曲线（如直纹、测地线），或者需要同时维护一张到平面的展开图。这带来两个麻烦：一是当直纹分布随形变而改变时难以处理；二是难以描述”没有沿直纹裁剪、带任意边界”的可展条带，放样这类经典问题因此不易求解。此前一直缺少一个作用于一般四边形网格的局部可展性判据。
• 本文 idea：借用”可展 \(\Leftrightarrow\) 第二基本形式降秩”这一微分几何刻画，并沿用棋盘格（checkerboard）思想，把它离散成”每个面一条方程”的局部条件。网格的边不必对齐任何特殊曲线，也不需要展开图，因此能无缝嵌入 remeshing、trimming、merging 等标准建模流程。

方法

整体思路：给任意四边形网格的每个顶点赋一个正权重，由此在每个面内嵌一个共面四边形，其所在平面即该面的切平面 \(\tau_f\)、法向 \(\boldsymbol{n}_f\)，中心为接触点 \(b_f\)——”接触元”就此定义。相邻面切平面的交线给出直纹向量；当一个面四周的平均直纹向量彼此平行时，该面即为离散可展。全部条件写成约束的平方和，用 Levenberg–Marquardt 优化。

flowchart LR
  A["带顶点权重的四边形网格"] --> B["每面内嵌共面四边形 → 接触元 (切平面/法向/接触点)"]
  B --> C["相邻面切平面求交 → 每条边的直纹向量"]
  C --> D["每个面的平均直纹向量"]
  D --> E["可展性条件: 平均直纹向量平行"]
  E --> F["优化 (约束+公平项) 投影到可展流形"]
  F --> G["加权 Doo-Sabin 细分 → 有理双二次样条曲面"]

关键设计一：由顶点权重定义接触元。 对每个面 \(f=v_0v_1v_2v_3\)，在每条边上取一个分点 \(m_{e_i}\)。一个引理指出：四个边分点共面当且仅当四段距离比之积为 1。而用顶点权重 \(w_i>0\) 写成 \(m_{e_i}=\dfrac{1}{w_i+w_{i+1}}(w_i v_i + w_{i+1} v_{i+1})\) 时，该乘积恒为 1，即任意正权重都保证边分点共面，反之任意共面边分点都可由权重表示。用齐次坐标 \(\hat v_i=(w_i v_i, w_i)\in\mathbb{R}^4\) 表示带权点后，边分点就是相邻齐次坐标的中点，接触点 \(b_f\) 是四顶点的加权质心。权重提供额外自由度，对粗网格的可展化效果尤其明显。

关键设计二：直纹向量与每面可展条件。 对相邻两面 \(f_i, f_{i+1}\)，其切平面交线 \(r_{e_i}=\tau_{f_i}\cap\tau_{f_{i+1}}\) 是离散直纹，方向由 \(\boldsymbol{r}_{\vec e_i}=\boldsymbol{n}_{f_i}\times\boldsymbol{n}_{f_{i+1}}\) 给出。绕一个面把对边的直纹向量做加权平均，得到两个平均直纹向量 \(\boldsymbol{r}_{\vec e13}, \boldsymbol{r}_{\vec e02}\)。离散可展性定义为：每个面上这两个平均直纹向量平行，即 \(\boldsymbol{r}_{\vec e13}\times\boldsymbol{r}_{\vec e02}=0\)。它同时蕴含了”每个面只有单一直纹方向”“该面对应的离散第二基本形式行列式为零 \(\det(\mathrm{II}_f)=0\)”等多条性质，是对光滑可展多种特征的同步离散化。平坦区域、拐点直纹、沿直纹排布的面等特殊情形都自然满足该条件。

关键设计三：优化框架与公平项。 变量含顶点、顶点权重、每面法向、每边直纹向量等；把上述几何关系写成一组多项式次数不超过 2 的约束（顶点、法向、平均边向量、直纹、可展 \(E_\mathrm{dev}\) 等能量），再加上顶点/法向/直纹场的公平（fairness）项抑制抖动与几何奇异，目标函数为各项加权和。对法向场施加公平能显著防止如尖锐脊线这类几何奇异。迭代末尾把正则项权重清零，让 \(E_\mathrm{dev}\) 真正趋于零。设计工具在此基础上叠加：位置约束（交互拖拽）、滑移约束（曲面沿参考曲线/点云滑动，用切平面距离近似）、软等距约束（模拟不可拉伸或弹性材料的正则化）。

关键设计四：细分与多分辨率。 加权 Doo–Sabin 细分作用在齐次坐标 \(\hat v_i\) 上时，会插值接触元（保持面心与切平面不变），其极限是一张分片有理双二次样条曲面（水密），控制点即最初的带权顶点。这让方法能导出可转 NURBS 的光滑曲面，并支持多分辨率建模：用一张小的控制网格 \(M\)，令细分 \(k\) 次后的 \(S^k M\) 离散可展，从而以少量 Bézier 面片表示近可展样条曲面。此外可证明接触元网格能对可展曲面做到二阶（Taylor）逼近，配合投影变换（把抛物柱面映射为二次锥面）还能沿整条直纹二阶逼近。

实验结果

方法实现为 Rhinoceros3D 的 C++ 插件，用 oneMKL/pardiso 求解稀疏系统，交互设计时优化在后台持续运行。可展性主要通过高斯像是否退化成曲线、以及正齿面（orthotomic）是否退化成曲线来验证（作者指出直接数值计算高斯曲率误差过大、不可靠）。一个突出性质是可展判据在 remeshing 与投影/仿射变换下近乎不变——把一张可展网格用另一张公平网格重新表示，仍近似满足判据。下表取该”重网格化 / 投影不变性”实验，报告变换后残余可展能量 \(E_\mathrm{dev}\) 仍极小：

配置	面数 \(\lvert F\rvert\)	总 \(E_\mathrm{dev}\)	每面 \(E_\mathrm{dev}\)
原网格 (a)	734	\(5.5\times10^{-5}\)	\(7.4\times10^{-8}\)
重网格化 (b)	175	\(2.0\times10^{-4}\)	\(1.1\times10^{-6}\)
重网格化 (c)	900	\(1.9\times10^{-5}\)	\(2.1\times10^{-8}\)
投影变换 (d)	734	\(2.5\times10^{-5}\)	\(3.4\times10^{-8}\)

此外：抛物柱面/二次锥面等能被离散网格精确表示（残余能量低至 \(10^{-28}\sim10^{-30}\) 量级）；把先前工作（正交测地网方法）的结果再经本文优化，高斯像质量可进一步提升；交互放样、含奇异点与折痕的分片可展、多分辨率建模等均在文中给出算例，单次迭代耗时从毫秒到亚秒级，交互可用。

亮点与局限

• 亮点：
  • 首个作用于一般四边形网格的局部可展性判据，边无需对齐直纹或测地线，也无需维护展开图，天然兼容 remeshing / trimming / merging。
  • 判据在重网格化、投影/仿射变换下近乎不变，便于裁剪与拼接。
  • 可直接导出水密的有理双二次样条曲面（加权 Doo–Sabin 细分极限），利于转 NURBS、接入 Rhino 工作流。
  • 能优雅求解经典的可展放样问题，并可表达折痕、奇异点等设计元素。
• 局限：
  • 优化需要用户有意义地选取各能量项权重，对不熟悉可展曲面特性的用户门槛较高。
  • 目前不与”曲面用分片可展逼近”类工作比较逼近能力，因为那需要先做分片可展分解；这是为保留用户引导式设计而有意为之。
  • 放样问题本身病态：某些边界（如两条异面直线）在特定直纹走向下无解，关闭直纹公平项时优化会陷入局部极小、在边界产生奇异。

延伸思考

• 该工作把”接触元”从微分几何概念落成可优化的离散建模基元，作者认为可发展出更完整的接触元网格理论；沿此思路，后续在 Laguerre / Möbius 球几何中用球面片、圆边等推广离散微分几何的工作，与本文属同一研究脉络。
• “降秩第二基本形式”的离散化思想与高度场 Hessian 秩最小化、高斯像细化等可展逼近方法互补：本文强于交互式设计与形变，而后者强于自动的分片可展逼近，二者结合或能兼顾”自动分解 + 设计美学”。
• 判据的投影不变性来自齐次坐标线性表达，这提示：在带权/齐次框架下处理几何约束，可能让更多”投影不变”的几何性质（如共轭网、直纹族）获得同样简洁的离散优化形式。