《Digital Discovery》：A graph-based approach to selection of feasible compositions for compositionally graded alloy design
编辑推荐：
成分梯度合金（Compositionally Graded Alloys, CGAs）为需在空间变化条件下服役的工程构件提供了一种有前景的解决方案。然而，由于潜在的化学不相容性、力学失配或热膨胀差异可能危及可制造性或部件性能，沿梯度选择合适成分面临挑战。本研究
成分梯度合金（Compositionally Graded Alloys, CGAs）为需在空间变化条件下服役的工程构件提供了一种有前景的解决方案。然而，由于潜在的化学不相容性、力学失配或热膨胀差异可能危及可制造性或部件性能，沿梯度选择合适成分面临挑战。本研究在基于图论的CGA设计计算方法的最新进展基础上，解决了设计空间中存在多个孤立可行成分区域时如何进行操作的问题。采用五元Nb–Cr–V–W–Zr体系，利用CALPHAD（PHAse图计算）模拟、混合律模型和经验估算相结合的方法计算了热学性质和力学性质。所得数据嵌入标号属性图（Labeled Property Graph, LPG）结构，并通过约束条件进行筛选，这些约束降低了凝固相关缺陷的可能性、促进了梯度上的相稳定性，并确保了增材制造（Additive Manufacturing, AM）的工艺适用性。一组约束条件产生了四个孤立的子图，包含多组相连成分，研究人员进一步利用关键性能指标评估了这些子图在CGA设计和制造背景下的各自优势与不足。总体而言，本工作建立了子图分析与选择方法，为将设计空间缩小至一个符合性能目标和制造约束的可遍历成分组提供了指导，适用于任何CGA设计场景。
本研究由Arroyave研究团队开展，发表在《Digital Discovery》期刊上，核心贡献在于提出了一套系统化的数据驱动框架，用于在高维合金设计空间中分析和选择适用于成分梯度合金设计的子图。
**研究背景与动机**
单质合金长期以来作为工程材料的基础，提供均匀的中间性能，但无法满足经历空间变化服役条件的构件对局部性能定制的需求。现代构件常面临复杂的热梯度、力学载荷或腐蚀环境，单一均匀合金难以兼顾。成分梯度合金通过在同一构件中实现化学成分的渐进变化，使材料在不同区域展现差异化的微观组织和性能，避免了突变界面引起的力学或热学不相容问题。近年来，增材制造技术的进步使具有平滑成分过渡的梯度合金制造成为可能。然而，CGA设计面临重大挑战：力学、热学和化学性质的相互作用复杂，改善某一性质往往以牺牲另一性质为代价；合金系统间的相容性至关重要，热膨胀系数失配可能导致残余应力、开裂或界面剥离。传统路径规划算法如快速探索随机树（Rapidly-Exploring Random Trees, RRT
*）虽已被应用于CGA研究，但在过约束设计空间中计算昂贵且难以解释。这促使研究人员寻求更高效、可扩展的框架来处理约束驱动的CGA设计环境。标号属性图（LPG）作为一种图结构，由节点和边组成且均可存储键值对属性，能够直接表征成分值、工艺参数或性能指标等描述性属性，形成更紧凑、信息丰富的结构以支持复杂关系的高效查询。Allen等人此前已证明LPG在金属增材制造中的效用，本研究在此基础上扩展，引入子图分析技术以识别和选择适用于CGA设计的可行成分区域。
**主要技术方法**
研究采用的技术路线包括：（1）基于CALPHA框架，使用Thermo-Calc的TC-PYTHON API开展热力学模拟，获取相稳定性数据、固相线/液相线温度及平衡相分数；采用混合律估计Pugh比，运用Maresca-Curtin模型估算1000°C下的屈服强度，采用Nabarro-Herring模型评估蠕变速率，并基于Scheil模拟计算Kou热裂判据。（2）将数据嵌入标号属性图结构，节点代表合金成分及其属性，边连接成分相似的节点；通过约束过滤（凝固温度范围≤50°C、平衡BCC相分数≥0.999、Scheil BCC相分数≥0.999）提取满足条件的子图。（3）采用均匀流形近似与投影（Uniform Manifold Approximation and Projection, UMAP）进行降维可视化，结合核密度估计（Kernel Density Estimation, KDE）验证聚类模式的真实性，使用NetworkX进行图操作与子图提取。样本数据来源为五元Nb–Cr–V–W–Zr体系的约10,000个独特成分，元素选择基于其高温合金研究相关性：Nb稳定BCC相，V提供韧性和固溶强化，W提供高熔点和抗蠕变性，Cr提供抗氧化性，Zr提供耐腐蚀性。
**研究结果**
"凝固范围约束"：将凝固温度范围限制在≤50°C后，10,274个节点缩减为352个，形成5个子图（其中4个有效）。子图0包含293个节点，提供多条可遍历路径；其余子图节点数分别为37、14、8。UMAP投影显示各子图的成分分布特征，密度编码揭示了特性变化。
"平衡BCC相分数约束"：在12个温度点（25°C至2750°C）要求平衡BCC相分数≥0.999，每个温度点均产生单一子图，显示完全连通的成分空间。1000°C时可在Cr富集区至V富集区再至Nb和Zr富集区规划多步路径；2500°C时则可设计更短的高温耐受路径。节点数量随温度变化呈峰形，约1700°C时达到峰值，反映熵效应对BCC相稳定性的促进作用。
"Scheil BCC相分数约束"：这是最具区分度的约束，产生四个明显分离的子图，证明了单独使用UMAP无法识别可行空间中间断区域、必须借助图论表示的关键价值。子图0（118个节点）以V、Cr、Nb为主要成分，W含量低；子图1富V、Nb；子图2富W；子图3富Nb、含少量Zr和W。
**各子图材料性能评估**
"Kou热裂判据评估"：子图0的Kou值最低，热裂敏感性最小；子图2和3的Kou值较高，制造过程中裂纹形成风险更大，但各子图内部均存在可遍历的低Kou值区域。
"平衡密度分析"：子图0密度最低且分布窄，适合轻量化应用；子图1密度约8.5 g·cm
?3；子图2和3密度范围宽，有利于设计梯度密度分布。
"热膨胀系数评估"：四个子图呈现阶梯式CTE值分布，各子图内部CTE变化范围相似，内部CTE失配不构成设计顾虑，但与其他子图的匹配需关注。
"力学性能分析"：子图3屈服强度最低（峰值约100 MPa），与成分复杂度导致的固溶强化效应一致；子图1和2屈服强度最高但分布宽，路径选择灵活性大。蠕变性能方面，各子图范围相似，子图2蠕变速率分布最宽、最大值最高，但整体最低值最优。Pugh比方面，子图2值低且分布窄（中心约2.25），显示潜在脆性；子图3全部高于延性-脆性阈值1.75，延性最佳。
"性能权衡综合分析"：雷达图和箱线图分析显示，子图0综合平衡，低密度、中等热裂敏感性、高CTE、低屈服强度；子图1热裂抗力好、CTE有利但密度较高；子图2高强度、热学性能优，但蠕变和密度变异大；子图3低密度、抗蠕变好，但机械强度弱。最终选择取决于应用优先级和路径规划策略。
**讨论与结论总结**
本研究的工作流程支持分阶段设计：先过滤绝对需求，再优先选择感兴趣的性能指标进行子图筛选。Scheil过滤产生的四个子图提供了充分的案例研究素材，展示了CGA设计场景下需考虑的关键因素。研究验证了所提出的图框架可用于识别满足可制造性和相稳定性约束的设计空间中的成分连通区域。将候选成分限制在单一可遍历子图内，即可设计连续的成分路径，预测其在平衡和非平衡凝固条件下均保持完全BCC相或具有窄凝固范围。该框架提供了从高维计算筛选到物理可实现的梯度材料设计的实用路径。
研究结论指出，通过整合基于CALPHAD的热力学建模、力学性能估算和基于图的数据结构，该框架支持对CGA设计和增材制造适用的子图选项进行高效探索。三种过滤约束的应用揭示了不同约束如何塑造设计空间，特别是Scheil BCC约束产生了四个可区分的可遍历子图。通过Kou热裂判据、密度、CTE、屈服强度、蠕变速率和Pugh比等关键性能的分析，详细展示了各子图的多指标表现及竞争性能间的权衡。子图0在整体上最为平衡，具有一致的低密度和中等热裂敏感性，但低屈服强度和高CTE是其局限。未来工作可将该图框架与闭环或自主实验平台集成，利用实测性质实时更新节点属性和约束评估，实现可行子图和成分路径的动态优化。
