Waste-to-Value: Reutilized Material Maximization for Additive and Subtractive Hybrid Remanufacturing

一句话总结

面向”增减材混合再制造”（ASHRM），这篇工作首次提出一个能严格保证无碰撞约束、适用于一般三维形状的计算框架：给定初始旧件与目标件，通过引入裁剪平面 + 两阶段优化，最大化可复用的旧件体积，把”增材从零打印目标件”转化为”尽量沿用旧料、只补差量”，在 12 个案例上平均节省 69% 材料。

研究背景

再制造（remanufacturing）把磨损、报废或旧型号的零件恢复成近新甚至更优的状态，核心目标是最大化复用材料，尤其对航空航天、模具等高价值材料场景意义重大。增减材混合制造（ASHM）同时具备加料（增材 AM）与去料（减材 SM）能力，能处理复杂形状，正越来越多地用于再制造（ASHRM）。

ASHRM 流程规划的关键任务，是在保证可制造的前提下确定三类体积并最大化复用体积：

• 减材体积 \(V_S\)：从初始件上被切除的部分；
• 增材体积 \(V_A\)：需要补打的部分；
• 复用体积 \(V_R\)：保留下来直接沿用的旧料。

理想情况下，初始件 \(M\) 与目标件 \(M^*\) 做布尔运算，交集 \(V_R = M \cap M^*\) 即为理想（最大）复用体积，\(V_S = M - V_R\)，\(V_A = M^* - V_R\)。

难点在于：对复杂模型，这种理想切分往往违反制造过程中的无碰撞约束，而现有 ASHRM 规划普遍忽视这些约束，导致实际加工时刀具/喷嘴无法到达而失败。约束的复杂性来自两方面：

• 带动态约束的体积优化：可表述为寻找体积密度函数以最大化复用体积的广义变分问题；由于已去除区域会退出后续约束，约束项随优化过程动态演化。
• 异质的无碰撞约束：局部约束（喷嘴出口与初始打印基面的碰撞）与全局约束（整个刀具/喷嘴组件与动态存在材料的碰撞）性质迥异，显著放大搜索空间。

已有工作要么面向 CAD 特征、无法用于自由曲面，要么只处理简单高度场曲面，或只考虑全局碰撞、用贪心法处理简单几何。本文是首个在一般形状上严格考虑无碰撞约束的 ASHRM 规划方法。

方法

给定初始模型 \(M\) 与目标模型 \(M^*\)，算法先确定二者相对摆放，再生成一套可制造的再制造流程规划，使复用体积最大。采用 3+2 轴定位式制造策略（加工前固定旋转轴，过程中只动直线轴），以降低规划复杂度并提升稳定性。

两类无碰撞约束

• 可达性约束（全局）：加工刀具或打印喷嘴与动态变化的存在材料之间不能碰撞（增材可达 + 减材可达统称可达性）。
• 平凸约束（局部，plano-convex）：在 3+2 定位增材下，喷嘴必须打印在分片平坦且凸的表面上，才能避免喷嘴出口与初始打印基面 \(B\)（即 \(V_R\) 与 \(V_A\) 的接触面）碰撞。

当初始布尔切分得到的打印基面违反平凸约束时，引入一组裁剪平面把理想复用体积切开，使新的打印基面满足平凸约束。

问题重述：从变分到组合优化

核心思路是用裁剪平面 \(\mathcal{C}\) 替代连续体积密度函数，把广义变分问题降维成组合优化：

\[\arg\max_{v(x)} \int_{\Omega} v(x)\, dx \quad \Longrightarrow \quad \arg\max_{\mathcal{C}} V_R(\mathcal{C})\]

约束为 \(g_{plano\text{-}convex}(V_A)=0\) 与 \(g_{accessibility}(V_S, V_A)=0\)，其中 \(V_R(\mathcal{C})\) 是被 \(\mathcal{C}\) 裁剪后的复用体积。整体是一个两阶段算法。

第一阶段：复用体积初始化（只管平凸约束）

目标是找一组既满足平凸约束、又尽量大的复用体积的裁剪平面。

• 生成候选裁剪平面：对初始打印基面的每个连通分量 \(b_i\)，采样 \(N\) 个来自上半球的法向 \(n_{i,j}\)，取 \(b_i\) 上沿该法向到加工底座投影距离最小的点 \(p_{i,j}\)，构成候选平面 \(c_{i,j} = \langle n_{i,j}, p_{i,j}\rangle\)。每个候选平面切割理想复用体积，得到新的打印基面。
• 最优组合选择——加权集合覆盖（WSCP）：某些裁剪平面能同时解决多个连通分量。为每个 \(c_{i,j}\) 赋集合 \(s_{i,j}\)（若 \(b_k\) 与加工底座位于 \(c_{i,j}\) 异侧则将 \(b_k\) 加入），并赋权重 \(W_{i,j} = 1 - Volume(V_R^{i,j}) / Volume(V_R)\)。通过整数线性规划求解：

\[\min \sum_{i=1}^{N}\sum_{j=1}^{M} W_{i,j}\cdot Z_{i,j}\]

约束为 \(Z_{i,j}\in\{0,1\}\)，且每个初始打印基面 \(b_k\) 至少被一个选中的裁剪平面覆盖。用 GLPK 求解器得到最优组合，从而得到满足平凸约束的复用体积 \(V_R^{\mathcal{C}}\)。

第二阶段：复用体积最大化（引入可达性约束）

基于第一阶段定义 SM-then-AM 体积（先切后补的部分）：

\[V_T^{\mathcal{C}} = V_R - V_R^{\mathcal{C}}, \quad V_A^{\mathcal{C}} = V_A + V_T^{\mathcal{C}}, \quad V_S^{\mathcal{C}} = V_S + V_T^{\mathcal{C}}\]

• 基于图的可达性分析：构建碰撞依赖图 \(G_{cd}\)（预处理中通过几何碰撞检测建立），在优化过程中只需从图上删除 \(V_T^{\mathcal{C}}\) 内的节点及其边，再检查是否还有不可达节点，即可验证可达性，完全避免优化循环中的重复几何碰撞检测。
• 可行复用体积变换：迭代地沿法向 \(n_i\) 平移每个裁剪平面（步长取该法向到加工底座最短投影距离的 30%），每次迭代后用 \(G_{cd}\) 检查可达性，直到得到可行复用体积。
• 全局平面转局部平面：无限大的全局裁剪平面限制了精细优化，于是把每个全局平面拆成一组局部裁剪平面 \(\mathcal{C}_L = \{c_1^L, c_2^L, \dots\}\)，每个局部平面只裁剪 \(V_R\) 的局部区域（表示为 3D 平面上的四边形）。
• 蒙特卡洛树搜索（MCTS）：用 MCTS 探索局部裁剪平面位姿的巨大搜索空间。树的每一层对应一个局部平面 \(c_i^L\)，层数等于 \(|\mathcal{C}_L|\)；节点状态为该平面的法向与原点，节点值记录奖励与被选次数。经选择、扩展、模拟、回传四阶段迭代，配合 UCB 函数引导搜索优先选中能产生更大 \(V_R^{\mathcal{C}_L}\) 且满足约束的位姿。达到迭代阈值后，取满足可达性与平凸约束、且 \(V_R^{\mathcal{C}_L}\) 最大的解。

摆放确定（预处理）

\(M^*\) 相对 \(M\) 的摆放直接影响可制造复用体积，而单纯最大化 \(V_R\) 并不能保证最优（约束可能让大 \(V_R\) 对应更小的 \(V_R^{\mathcal{C}_L}\)）。采用解耦旋转与平移的由粗到细策略：粗阶段用第一评估函数（结合 \(V_R\) 体积与鼓励”直立而非平躺”的项，使打印基面更高）打分，选前五名旋转；细阶段对这五个固定旋转微调平移，用第二评估函数（进一步考虑打印基面相对高度与 \(V_R\) 高度）打分，取最高者。

可制造块分解（后处理）

• 减材块分解：\(V_S\) 表面已满足全局可达，用 UG-NX 沿 \(\pm X/\pm Y\) 四个方向生成定位粗加工块与精加工刀路，剩余不可达体积再计算额外可达方向分解。
• 增材块分解：依次用局部裁剪平面切 \(V_R^{\mathcal{C}_L}\) 得增材块，每个局部平面作为打印基面进一步分解成自支撑块；若块间仍有碰撞，沿打印方向切片、做碰撞检测、建有向图记录依赖，拓扑排序后进一步分解成无碰撞块。
• 刀路生成：反转块序列，连同方向输入 UG-NX 生成 SM/AM 刀路。

实验结果

C++ 实现，用 CGAL、Libigl、Eigen 做几何处理；碰撞检测 AM/SM 各采样 100 方向，第一阶段评估 100 个法向，第二阶段 MCTS 迭代 500 次，自支撑阈值 50°。

• 计算结果：在 12 个模型（涵盖模型修复、缩放、备件再制造、类型变换四类应用）上验证，平均节省 69% 材料。当 \(V_R\) 刀具/喷嘴均可达、打印基面平坦且位于顶部时，最终复用体积接近 \(V_R\)（如 Bunny-Cow，\(R_F=0.73\)，\(R_V=0.74\)）；Bearing 案例最终复用比 \(R_F=0.90\)，节省最高达 90%。平均耗时约 8 分钟，MCTS（第二阶段）是主要瓶颈。
• 大规模随机测试：从 Thingi10K 与经典图形模型库随机配对做初始/目标件，展示 50 个成功案例；几何与尺寸相近时复用率高，差异大时因难找好摆放且受可达性约束，\(R_V\) 到 \(R_F\) 下降明显。
• 刀具尺寸影响：更长的刀具更易可达，可得到更大的 \(V_R^{\mathcal{C}_L}\)，说明算法能适应不同刀具尺寸。
• 策略对比：相同时间下，两阶段策略优于自顶向下（全局平面搜索低效）与自底向上（单平面搜索空间巨大）；WSCP 用整数线性规划求最优解，优于单平面法与贪心法；带 UCB 的 MCTS 明显优于随机搜索（后者大量位姿违反约束）；摆放策略上，由粗到细解耦优于更密采样的联合优化——后者在 Hand-Hand 上因缺”直立”启发式选了倾斜摆放导致零复用，且本文方法计算更快。
• 摆放消融：只看体积会让大量 \(V_R\) 被迫转成 SM-then-AM 体积；只看平凸约束又使 \(V_R\) 偏小；本文用 \(w_1=0.4, w_2=0.6\) 权衡二者，得到兼具较大 \(V_R\) 与最大 \(V_R^{\mathcal{C}}\)。
• 仿真与物理验证：用 UG-NX 生成粗/精加工与增材刀路，全部案例通过 NX 的碰撞检测仿真。受限于昂贵金属 5 轴 ASHM 机的可得性，在改装的桌面 5 轴 CNC 上用 PLA 做物理验证，成功再制造 Bearing 与 Sculpture，相比纯增材分别省料 90% 与 69%。作为首个对一般模型做物理验证的工作，结果存在块间接缝、增材表面质量差等缺陷（部分源于设备改装导致旋转轴错位与挤出不稳），可通过 AM 后追加 SM 精加工改善。
• 制造策略对比：纯 AM 浪费约 9 倍材料；纯 SM 不额外耗料但留下缺失区域；分件打印再胶接引入人工、定位不准且更费时；本文 ASHRM 综合最优。

亮点与局限

亮点：

• 首个严格考虑无碰撞约束、适用于一般形状的 ASHRM 规划框架，并做了物理验证。
• 降维重述：用裁剪平面把病态的体积变分问题转化为可高效求解的组合优化。
• 约束分治：局部平凸约束用裁剪平面 + WSCP 处理，全局可达性用碰撞依赖图 \(G_{cd}\) 处理，避免优化循环中反复做几何碰撞检测。
• 全局转局部平面 + MCTS：局部裁剪带来更精细的复用体积提升，MCTS + UCB 高效探索巨大位姿空间。

局限：

• 高度重叠表面：\(M\) 与 \(M^*\) 仅差一个孔洞等高相似情形下，布尔运算易出错、边界定义模糊，连鲁棒布尔方法也难解决。
• 未优化制造时间：算法只优化省料，不优化块数或刀路时间。
• 可扩展性：运行时间随三角面片数近似线性增长，超高面数模型可能需数小时，瓶颈在布尔运算等密集几何操作；且仅适用于水密实体模型，薄壳模型的混合制造仍待探索。

延伸思考

• 当前是单次 SM-then-AM，未来引入多段 SM-AM 交替序列有望进一步扩大复用体积，但会大幅增加规划复杂度与约束耦合。
• “把连续变分问题约束到裁剪平面这一低维、组合式参数空间”的思路，本质是用可枚举的几何基元换取可求解性，这类降维参数化在其他受制造约束的几何优化中或可迁移。
• 可达性用碰撞依赖图做增量更新、避免重复几何检测，是把几何约束”预编译成图查询”的典型做法，对交互式或大规模制造规划的加速很有启发。
• 省料与制造时间/表面质量之间存在明显权衡，把制造时间目标纳入优化（如最小化块数）是走向实用的自然方向。