1. 课程导入与多尺度仿真全景

材料研发长期依赖“配方—工艺—性能”的试错循环,一次合金成分调整往往需要数周甚至数月的熔炼、热处理与力学测试。多尺度仿真的工程价值,在于把这一循环中的关键环节搬到计算机里:用电子尺度的计算判断相稳定性与键合强度,用原子尺度的模拟观察位错形核与界面迁移,用介观尺度的方法追踪晶粒长大与相变组织演化,再用连续介质尺度的有限元评估构件在真实载荷下的应力分布与失效位置。四个尺度并非彼此替代,而是通过参数传递与信息回填构成一条可追溯的“计算证据链”。

理解这条证据链,首先要建立尺度划分的直觉。电子尺度对应埃量级,处理的是电子密度、能带、总能量与原子间相互作用;原子尺度对应纳米量级,关注单个位错、裂纹尖端、晶界与点缺陷的行为;介观尺度跨越纳米到微米,处理的是位错群、晶粒、析出相与相场演化;宏观连续尺度则从微米延伸到米,用弹性、塑性、损伤本构描述工程构件的响应。尺度每放大一个量级,可模拟的时间与空间范围显著增加,但原子级细节被平均化,因此必须明确“哪些信息需要跨尺度传递、哪些可以粗粒化”。

材料多尺度仿真全景:尺度、方法、传递与工程应用 1 Å 1 nm 1 μm 1 mm 1 m 电子尺度 DFT / 第一性原理 电子密度、总能量 形成能、弹性常数 相稳定性、界面能 输出:能量/力学参数 原子尺度 MD / 分子动力学 位错形核与运动 裂纹扩展、界面迁移 扩散系数、层错能 输出:缺陷/界面行为 介观尺度 相场 / 位错动力学 晶粒长大、相变组织 析出相演化 位错群交互 输出:组织/本构参数 宏观连续尺度 FEM / 晶体塑性 应力应变场 损伤演化、失效位置 构件寿命评估 输出:工程性能预测 E, ν γ_sf, D σ_y, 组织 实验/宏观结果回填:修正势函数、界面参数与本构模型 典型工业应用场景与课程项目映射 高温合金叶片 γ/γ′相稳定性 位错攀移与蠕变 热障涂层应力 锂电正极材料 Li⁺迁移路径 相变与体积变化 颗粒开裂预测 汽车高强钢 马氏体相变 残余应力分布 碰撞吸能评估 半导体封装 界面分层与翘曲 热循环疲劳 焊点可靠性

方法谱系的选择取决于“要回答什么问题”。若关心合金中某相是否稳定、溶质原子倾向占据哪个位点,第一性原理是首选;若关心位错在应力下如何运动、裂纹沿哪个晶面扩展,分子动力学更直接;若关心热处理后晶粒尺寸分布与析出相体积分数,相场或蒙特卡洛方法更合适;若关心构件在循环载荷下的疲劳寿命,则需要晶体塑性有限元或宏观损伤模型。工程实践中,常见做法是“上层方法提供参数、下层方法消耗参数”:DFT 算出的层错能与界面能输入到 MD 势函数校验中,MD 得到的位错迁移率与界面迁移率输入到相场或位错动力学模型,介观模型输出的屈服强度与组织特征再输入到宏观有限元。反向路径同样重要——宏观试验测得的应力应变曲线可用于反推本构参数,再回填修正介观与原子尺度模型。

课程项目总览围绕一条完整的“计算证据链”展开。项目对象设定为一种假想的镍基高温合金,成分为 Ni-8Cr-5Co-2Mo-4Al-3Ti(质量分数,示例设定),目标是在 750 ℃ 下评估其蠕变抗力与组织稳定性。项目分为四个阶段:第一阶段用 DFT 计算 γ 与 γ′ 相的晶格常数、形成能与反相畴界能,为后续势函数提供基准;第二阶段用 MD 模拟 γ/γ′ 界面在 750 ℃ 下的迁移行为,提取界面迁移率与扩散系数;第三阶段用相场模型模拟时效过程中 γ′ 相的粗化动力学,得到不同时效时间下的颗粒尺寸分布;第四阶段将组织参数与蠕变本构输入晶体塑性有限元模型,计算单晶叶片模拟件在离心载荷下的蠕变应变与危险点位置。每个阶段都要求记录输入参数、计算设置、中间输出与判据,最终形成一份可追溯的仿真报告。

学习建议:本课程不要求你在第一遍就掌握所有方法的数学细节,但要求你始终清楚三件事——当前尺度在回答什么问题、输入参数从哪里来、输出结果将传递到哪个尺度。只要这条“传递链”不断,多尺度仿真就不会退化成四个孤立的计算练习。