首页/文章/ 详情

PMS | 综述: 面向适氢结构功能一体化合金设计的多尺度计算

5月前浏览831





导读

尽管先进适氢金属材料(AHMM)在抗氢、储氢、氢敏已取得诸多进展,但承载结构合金的氢脆问题、功能储氢合金的能力和循环充放限制仍然存在明显挑战,且受服役环境影响显著;现有计算模拟方法虽阐明了氢与晶体缺陷的基础作用,却缺乏对合金材料与氢相互作用过程中力学、化学、微观组织耦合行为的系统研究,难以支撑工程应用中的材料选型与寿命评估。本文聚焦合金材料适氢性能预测的核心难题,将AHMM分为氢暴露结构金属(HESMs)、氢结合结构金属(HISMs)、氢势功能金属(HPFMs)三类,梳理了多尺度计算设计体系,提出需根据不同材料类别差异化组合模拟方法,并指出未来需完善材料数据库、强化模型验证与实验协同,推动耐氢结构功能合金计算设计走向工程实用化,为AHMM研发提供理论参考。


文章封面 部分截图

   


部分内容  图文解读

图1清晰呈现本研究的核心框架及三类氢相关金属材料的具体应用场景:其中a为研究整体流程图,直观展示了从材料设计、计算模拟到性能验证的完整逻辑;b为HESMs的典型应用实例,如航空涡轮发动机等需在含氢环境中稳定运行的设备;c为HISMs的应用场景,包括金属氢化物、泡沫铜储罐等储氢设备;d为HPFMs的应用示例,如用于储氢系统的金属-石墨交替双极板,清晰展现各类材料的实际应用场景与功能定位。

   

图1 氢能应用中金属材料的多尺度表征涉及多个领域:(a) 本文概述不同计算设计领域中金属与氢的相互作用及储氢行为;(b) 航空涡轮发动机为 HESMs;(c) 金属氢化物与泡沫铜储罐为HISMs;(d) 金属 - 石墨交替双极板及其截面为 HPFMs。

微观组织会显著影响氢的吸附、扩散与释放,通过晶粒细化、相比例调控和缺陷调控,可兼顾储氢能力与抗氢脆性能。表 1 系统总结了三类AHMM的氢作用机理、关键性能及适配计算方法。

表1 AHMM与氢的相互作用总结

   

图 2a 示意了 Mg₂Ni 合金放氢过程的结构–功能关联与氢演变,金属氢化物(如 MgH₂、LaNi₅H₆)作为储氢功能材料依靠可逆反应储放氢,其储氢系统为集成配套组件的工程装置,二者设计重点不同;微观缺陷调控储氢表面反应,镁合金特殊晶界可降低氢解离能垒,晶界调控是优化功能金属表面性能的有效手段,此外可结合氢作用、电荷特征分析各类先进临氢金属材料,并通过多尺度耦合计算研究多晶材料氢输运行为(图 2b、2c)。

   

图 2:(a) 镁镍功能金属放氢 7000 秒后的储氢行为模拟计算结果;(b) 中锰钢马氏体与奥氏体结构的氢原子分布、电荷重排及差分电荷密度第一性原理计算结果;(c) 结合相场、密度泛函与分子动力学方法,整合原子尺度数据,完成多晶二氧化钛氢输运的多尺度模拟。

密度泛函理论(DFT)作为氢能材料设计的核心手段,可解析化学储氢合金的键合与表面相互作用,针对镁基、钯基等合金揭示氢相关特性并指导材料开发,第一性原理计算还能关联电子结构与宏观性能、计算合金相生成焓、表征氢在金属表界面的吸附行为,依托金属氢化物稳定性与成分结构的关联优化储氢性能,适配三类临氢金属材料研发需求,为多尺度模拟提供基础参数,同时可搭建量子尺度与宏观性能的关联,结合相关理论参数预测合金氢相关动态指标,文中图 3a-e 分别展示了其相关作用机理、优化效果及计算结果。

   

图 3 基于 DFT 可视化不同场景下的氢 - 金属相互作用:(a) 金属表面氢吸附电荷分布;(b) 氢原子重叠的电子密度,体现共价键特征;(c) 电子局域化函数,表征键合与电子分布;(d) 储氢功能金属的优化晶体结构;(e) 三元合金焓值曲面的内能与生成焓计算。

探究氢的输运及其对先进临氢金属材料的影响,对储氢技术发展和氢脆抑制至关重要;分子动力学(MD )与动力学蒙特卡洛( KMC) 是互补的核心模拟方法,MD 擅长解析原子尺度氢的扩散、缺陷作用及氢脆机制,获取关键参数但计算受限,KMC 依托相关参数模拟长时尺度氢行为、预测材料长期演化,二者耦合实现跨尺度仿真,广泛应用于各类临氢材料研究,可解析多种氢致相关问题,模拟精度受原子作用势等限制,图4 展示了二者的相关模拟应用。

   

图 4 为不同材料体系氢作用的多尺度模拟结果:(a) KMC 模型构建 NiO 与 YSZ 三维晶粒组织,作为氢还原与烧结模拟初始结构;(b) MD 模拟 430 不锈钢中氢原子与螺位错的时效演化作用;(c) 结合 KMC 分析氢渗入后纯铁单晶拉伸过程的裂纹扩展与应力分布。

晶体塑性模型(CPM)可有效表征多晶金属在氢作用下的非均匀变形与氢脆行为,纳入氢相关本构规律后,能预测局部应力集中、裂纹萌生与氢致损伤;图 5 结合实例直观呈现:氢会诱发 316L 不锈钢滑移局部化,可提高材料滑移系剪切应变率,降低合金界面能并加速裂纹扩展,同时显著影响镍、钯单晶的强度与加工硬化性能,为解析面心立方金属氢脆机理提供重要模拟手段。

   

图 5 基于CPM分析面心立方金属的氢致作用:(a) 结合金相、电子背散射衍射,得到 316L 不锈钢滑移系塑性应变与晶体取向分布;(b) 考虑氢效应,修正特定滑移系的剪切应变率本构关系;(c) 氢降低界面能,诱发 Al-Zn-Si 合金裂纹扩展,展示不同应变下的三维晶粒结构与取向变化;(d) 对比氢对镍、钯单晶剪切应力、拉伸应力及加工硬化率的影响。

人工智能结合 DFT、MD 等计算方法,可助力结构材料与储氢功能金属的研发。如图 6a 所示,XPEAK智能模拟平台可减少实验试错,定向优化材料动力学性能、指导理性设计。人工智能结合实验验证,是加速固态储氢材料实用化研发的有效途径。图 6b 展示了神经网络、支持向量机、决策树、线性回归等算法的单独应用,及其与各类计算方法的耦合联用模式。

   

图 6:(a) 基于 XPEAK 统计数据库元素组成,标注各元素出现频次;(b) 机器学习算法应用于氢能材料多类研究方向,结合仿真计算与材料性能,助力固态储氢体系研发。

多尺度耦合仿真方法将量子力学、随机模拟、热力学与连续介质力学相结合,统一应用于钢、锆基、镍基合金体系(图 7)。原子模拟提供微观机理参数,DFT 可计算电子结构能量、氢 - 缺陷结合能、固溶热力学与氢化物生成焓,并揭示氢扩散中的量子隧穿效应;MD 则表征氢与缺陷的动态作用、短程扩散及氢化物早期形核过程。DFT 与 MD 耦合研究还阐明了空位 - 氢作用机制,明确氢含量对空位形成能及缺陷形貌的调控规律。

   

图 7 多模型耦合仿真预测金属氢行为:(a) 相图计算结合 DFT,分析双相不锈钢钝化膜电子特性与氢脆关联;(b) 元素掺杂调控结构金属,优化氢渗透性能;(c) CPM 与 PFM 联用,研究充氢镍合金的断裂演化;(d) PFM 耦合 CPM,模拟锆合金氢致损伤演变;(e) DFT 结合有限元,计算镍基结构金属氢致结合能衰减。

图8系统展示AHMMs多尺度仿真全链条成果与现存挑战:a、b、c 分别直观呈现氢致晶格演变、氢化物相变规律及晶界氢敏感特性,验证了多尺度耦合方法在微观机理解析上的优势;d 揭示跨尺度计算误差、算力受限等现实短板;e 明确指出整合多种计算模拟策略与实验校验的一体化跨尺度框架,是提升材料预测精度、推动储氢合金定向设计的核心发展方向。

   

图 8 多方法耦合计算模拟结果:(a) 高性能储氢功能金属材料的晶格氢化改性机制;(b) 金属氢化物核壳结构形成过程计算;(c) 不同倾斜角度下稳定晶界构型,体心立方配位原子为蓝色,其余为白色;(d) 相图计算耦合相场、晶体塑性模拟,展示氢在高熵结构金属中的气泡效应;(e) 融合密度泛函理论、量子机器学习原子势、蒙特卡洛 / 分子动力学与实验验证的跨尺度计算设计框架。

多尺度耦合模拟是对真实物理体系的必要简化与近似表征,过程中必然存在信息损耗。如图9所示,通过不确定性量化、优化采样算法及创新耦合策略,精准识别、定量评估并合理降低这类信息损失,与单一计算方法的迭代优化同等关键。

   

图9 AHMM设计用多尺度计算方法的关键对比

图10完整展示了当前氢能金属材料研究的现存瓶颈与优化方向:a 至 e 依次反映合金成分效应、纳米氢分布、位错作用、光谱干扰及实验表征局限等理论模拟与实验验证的突出问题;f 提出微观结构调控与结构功能一体化设计策略,通过优化组织形貌协同兼顾氢吸附存储与塑性变形能力,结合多尺度仿真、先进实验技术与智能算法,可有效突破现有研究短板,助力高性能抗脆、长寿命氢能金属材料开发。

   

图10 AHMMs的多尺度表征策略:(a) 三维原子探针重构图(C 蓝、D 红)及 Fe、C、D 元素分布,揭示氢相关研究必要性;(b) 高分辨 TEM 显示室温下氢占据析出相及相界行为;(c) 电子通道衬度成像观测充氢后 AHMM 中位错运动;(d) 二次离子质谱表征金属污染及充氢前后能谱变化;(e) 电子背散射衍射与核平均取向差图,展示氢致退火态、5%/15% 预应变金属失效差异;(f) 耦合模拟方法需重点研究的氢作用相关未来方向

适氢合金的精准设计,核心难点在于多尺度计算参数的逐级实验验证,而非仅考核最终性能。可靠的多尺度模型,其原子能级能量、介观动力学系数、宏观本构关系等基础参数,都必须依托实验数据校准。图 11 将各类仿真输出与对应实验表征手段一一关联,实现模型校准直观可追溯,全面提升材料预测设计体系的可靠性。

   

图11 计算模拟结果与对应实验表征技术关联路线图

图12展示了仿真–实验–算法闭环迭代的材料自主优化设计框架,结合计算模型技术成熟度分级标准,以实验数据持续校准多尺度模拟、机器学习补齐研究盲区,打通理论研究与工程应用的壁垒,推动临氢合金的高效、可靠定向研发。

   

图12 面向 AHMM 设计融合先进表征技术未来方向:(a) 同步辐射 X 射线衍射表征晶格应变与织构演变;(b) 原位中子衍射研究充氢对堆垛层错的影响



文章信息

简而言之,最优计算策略取决于材料类别与服役挑战。对于HESMs,氢脆(HE)可以通过密度泛函理论(DFT)/分子动力学(MD)进行氢-缺陷定量表征,并借助晶体塑性模型(CPM)/相场法(PFM)/有限元法(FEM)建立失效模型,优先开展缺陷工程调控。同时,化学-力学计算流与基准数据集亦为裂纹预测验证所必需。对于HISMs,基于FEM/CPM的应力分析必须与长期扩散及界面稳定性模型(如动力学蒙特卡洛(KMC)、相图计算(CALPHAD和PFM)耦合,聚焦应力-扩散耦合行为。面向物理储氢,需建立融合氢传输、界面稳定性与循环损伤的环境感知本构模型,并与高压、低温条件下的实验相关联。对于HPFMs,通过DFT/机器学习(ML)进行电子结构筛选以驱动材料发现,同时采用MD/KMC优化传输动力学。研究重心转向通过氢化物形核与生长的介观尺度模型实现功能耐久性,兼顾热力学与微观组织架构设计,以缓解化学储氢中的退化问题。尽管存在通用的计算模拟方法集,但其集成方式必须针对具体材料类别进行定制。


文章信息


主要作者情况:

吴仁豪(共同第一作者,通讯作者), Zaigham Saeed Toor(共同第一作者), Hyoung Seop Kim(通讯作者)


通讯作者单位:

1.韩国浦项科技大学金属与环保材料研究所

2.日本东北大学材料科学高等研究所


文章标题:

Multiscale Computational Design of Hydrogen-Compatible Integrated Structural–Functional Alloys


期刊名称:

Progress in Materials Science


文章链接:

https://doi.org/10.1016/j.pmatsci.2026.101727

来源:我的博士日记
断裂化学通用航空电子增材裂纹理论材料分子动力学多尺度人工智能曲面
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-05-12
最近编辑:5月前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 133文章 155课程 0
点赞
收藏
作者推荐

任意多晶微观结构生成,GUI操作,模型直接下载

在金属材料、陶瓷及复合材料的微观力学研究中,构建一个符合统计学特征的多晶代表性体积单元(RVE)往往是科研工作的第一步。然而,传统的建模方法往往面临重重困难:使用商业软件手动分割效率低下;利用专业建模软件(如 Neper)虽然强大,但命令行操作和复杂的参数配置让许多初学者望而却步;而自编程序生成 Voronoi 镶嵌模型,又难以精准控制晶粒尺寸分布和形状统计特征。有没有一种工具,既能保证模型的科学性,又能像“点外卖”一样简单快捷?今天,我们要向大家强烈推荐一个在线神器——Synthetmic。【核心介绍:什么是 Synthetmic?】Synthetmic 是由赫瑞-瓦特大学(Heriot-Watt University)的 David Bourne 博士开发的一款基于 R 语言 Shiny 框架的在线 GUI 工具。它的核心使命是通过数学算法,快速生成具有特定几何特性的合成多晶微观结构。该工具不仅支持传统的 Voronoi 镶嵌,更引入了功能强大的 Laguerre 镶嵌(权重 Voronoi)算法。这意味着你不再受限于匀称的晶粒,而是可以生成具有特定体积分布、更接近真实金属组织的复杂模型。网站地址:https://david-bourne.shinyapps.io/synthetmic-gui/【功能亮点:为什么它值得收藏?】零门槛,全在线操作: 无需安装任何环境,打开浏览器即可完成从参数配置到模型生成的全过程。高度可定制的统计控制: * 晶粒数量: 自由设定生成 10 到 1000+ 个晶粒。尺寸分布: 支持对晶粒体积的对数正态分布(Log-normal)进行精准控制,模拟不同加工状态下的组织。空间排布: 通过调整点过程参数,控制晶粒的密集程度与均匀性。实时可视化预览: 网页右侧提供 3D 实时渲染,调整左侧参数后,模型形态即刻更新,真正实现“所见即所得”。多格式导出: 生成的模型支持导出为坐标数据、拓扑连接信息等,方便后续导入 ABAQUS、ANSYS 或自编的有限元/晶体塑性(CPFEM)程序中。【操作流程:三步搞定】第一步:设定全局参数。 在左侧面板选择晶粒总数及 RVE 尺寸。第二步:精修几何特征。 调整权重系数(Weights)和偏度,生成不规则或特定分布的晶粒形状。第三步:导出与应用。 预览满意后,点击下载按钮获取几何模型文件。五大基础案测试例如下:一:2D单相对数正态分布1000个晶粒模型(尺寸1*1)(0.1s)二:2D双相对数正态分布500个晶粒模型(0.1s)(体积分数0.2:0.8,晶粒个数450:50)自动生成对应的统计信息:点击下载可以直接获取所有相关的晶粒信息三:随机的500个晶粒的3D单相模型(0.1s)四:随机的1000个晶粒的3D双相模型(0.5s)二次相分布于晶界五:随机的2000个晶粒的周期性单相模型(2s) 同时内置丰富的可视化,比如调整透明度,切片显示,修改颜色等等来源:我的博士日记

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈