《METALLURGICAL AND MATERIALS TRANSACTIONS A-PHYSICAL METALLURGY AND MATERIALS SCIENCE》：A Precipitation Model for Simulation of Mechanical Properties of Al–Si–Mg–Cu Foundry Alloys
编辑推荐：
针对Al–Si–Mg–Cu铸造铝合金，研究人员开发了一种析出模型（precipitation model）。该类合金的时效响应强烈受先前加工条件的影响，包括铸造参数以及后续的固溶处理（solution heat treatment, SHT）。该模型首先预测铸
针对Al–Si–Mg–Cu铸造铝合金，研究人员开发了一种析出模型（precipitation model）。该类合金的时效响应强烈受先前加工条件的影响，包括铸造参数以及后续的固溶处理（solution heat treatment, SHT）。该模型首先预测铸态（as-cast, AC）微观组织，特别是初生颗粒的体积分数与尺寸，以及在不进行固溶处理时Mg和Si的固溶度水平。当包含固溶处理时，模型模拟关键微观组织参数的演变以及整个固溶处理过程中Mg和Si在固溶体中浓度的相应变化。这些预测结果确定了进入最终人工时效热处理时的微观组织条件。时效热处理过程中的析出行为通过一个先前为变形合金（wrought alloys）开发的模型进行模拟，并进行了少量修改，以适应铸造合金体系的特征。模型预测结果与在170 °C下不同时效时间获得的拉伸性能测量值进行了对比，测试样品包括水冷或空冷条件下铸造后直接时效的样品，以及时效前经过中间固溶处理的样品。模拟结果与实验数据吻合良好。特别是，该模型成功再现了经固溶处理后相对于铸造后直接时效样品所观察到的屈服应力增加现象，该效应归因于固溶处理过程中含Mg颗粒的溶解以及由此导致的固溶体中可用Mg含量的增加。
铸造Al–Si–Mg–Cu合金力学性能模拟：从铸态组织到人工时效的集成析出模型解读
一、研究背景与问题提出
Al–Si–Mg–Cu铸造铝合金因其低密度、优良的铸造性能和可靠的力学性能，被广泛应用于汽车轮毂、悬挂部件、副车架及动力总成壳体等结构件中。这类合金兼具优异的流动性和抗铸造缺陷能力，特别适用于高压压铸（HPDC）、砂型铸造和金属型铸造等复杂几何形状的成形工艺。随着轻量化高性能铸件需求的日益增长，建立能够关联合金成分、加工条件与最终力学性能的预测模型变得愈发迫切。
然而，铸造合金的时效响应强烈依赖于此前的加工历史，包括铸造参数和后续固溶处理（SHT）条件。铸态组织中初生颗粒的分布、固溶体中合金元素的浓度水平，以及固溶处理过程中颗粒的溶解与浓度梯度的均匀化，都会直接决定最终人工时效阶段强化析出物的形成与演变。因此，准确预测时效强化效果的前提，是对铸造和固溶处理阶段微观组织演变进行可靠的定量描述。
此前，Myhr和Grong等人已针对6xxx系变形铝合金开发了名为NaMo的集成析出、屈服应力与加工硬化模型，并成功应用于预测Al–Mg–Si合金的析出结构和抗拉强度。但该模型直接应用于铸造合金时存在系统性偏差——模拟的析出动力学过慢，导致预测的峰值强度出现时间晚于实验观察。此外，铸造合金中Si含量远高于变形合金，而Mg含量通常较低，这种成分差异对空位保留和空位辅助扩散行为产生显著影响，进而改变时效析出动力学。
本研究旨在将NaMo模型扩展至铸造铝合金体系，建立覆盖铸造、固溶处理和人工时效全流程的集成模拟框架，以统一的方法预测合金从铸态到最终时效态的微观组织演变和力学性能。
二、研究内容概述
研究人员开发了一个完整的铸造Al–Si–Mg–Cu合金析出模型，该模型首先预测铸态微观组织，包括初生颗粒的体积分数和尺寸，以及Mg和Si的固溶度水平；当包含固溶处理时，模型模拟关键微观组织参数的演变和固溶体浓度的变化；最终通过改进的NaMo析出模型模拟人工时效过程中的析出行为。模型预测与实验数据进行了系统对比，验证了其预测铸造后直接时效和经固溶处理后时效样品屈服应力的能力。
三、关键技术方法
研究人员采用了一系列简化的解析表达式和数值模拟工具构建集成模型。铸态Si分布采用二元Scheil凝固模型估算，Mg固溶度通过Alstruc计算工具的系统模拟结果拟合得到经验关系式。固溶处理过程中Si和Mg浓度梯度的均匀化采用基于误差函数的扩散解析解描述，Si颗粒粗化采用扩散控制的Ostwald熟化动力学方程。人工时效阶段的析出行为基于Kampmann–Wagner数值框架的NaMo模型，通过引入动力学缩放因子ω=2.8修正空位辅助扩散效应，并单独建立冷却过程中粗大颗粒析出的形核律。屈服强度模型综合了固溶强化、位错强化、析出强化和相界强化四种贡献。实验验证基于Hydro研发中心的铸造试验，涉及三种不同Cu含量的合金，涵盖水冷和空冷两种冷却条件。
四、研究结果
4.1 铸态微观组织预测
模型采用Scheil凝固假设估算共晶分数和Si的固溶浓度，建立了铸态Si颗粒体积分数和平均半径与二次枝晶臂间距（SDAS）的关系。对于Mg，通过Alstruc模拟建立了铸态Mg固溶浓度与初始Mg含量和SDAS的经验关系式，发现Mg固溶度随SDAS增大而增加，这归因于较大SDAS对应更长的局部凝固时间，促进了Mg向枝晶基体的反向扩散。
4.2 固溶处理过程中的浓度均匀化
固溶处理期间，Si浓度梯度的均匀化通过扩散方程解析解描述，以SDAS作为特征扩散距离。Mg浓度的增加通过类似的误差函数解描述，其平衡浓度采用Arrhenius型表达式拟合Alstruc模拟结果获得。模型假设含Mg颗粒间的平均距离等于SDAS。
4.3 Si颗粒粗化动力学
固溶处理过程中Si颗粒的粗化采用扩散控制熟化动力学描述，粗化速率常数k?根据文献数据确定为2.3×10?? m3K/s。
4.4 时效析出模型的修正与验证
针对铸造合金与变形合金的动力学差异，研究人员引入动力学缩放因子ω=2.8，将形核和长大/溶解速率方程统一缩放，等效于引入有效时效时间t_eff=ωt。这一修正的物理依据是：高Si含量铸造合金在淬火后保留更多空位，Si-rich原子构型稳定了空位，从而加速了空位辅助扩散。冷却过程中粗大颗粒的析出通过单独的形核律描述，形核参数根据水冷和空冷样品的强度差异校准。
4.5 屈服强度模型
屈服强度由纯铝基体强度、固溶强化、位错强化、析出强化和相界强化线性叠加构成。位错强化考虑了冷却过程中因热膨胀失配产生的几何必要位错（GNDs）和塑性变形产生的GNDs。析出强化通过NaMo模型中的粒子尺寸分布直接耦合计算，初生颗粒的强化贡献采用上限估计。相界强化采用Hall–Petch型关系，以SDAS作为特征尺度。
五、讨论与结论
模型预测与实验数据在多种条件下均表现出良好的一致性。特别是，模型成功再现了经固溶处理后样品相对于铸造后直接时效样品屈服应力增加的现象，该效应归因于固溶处理过程中含Mg颗粒的溶解以及由此导致的固溶体中可用Mg含量的增加。这一结果表明，所建立的集成模拟框架能够有效捕捉铸造、固溶处理和人工时效全流程中微观组织演变与力学性能的关键关联，为铸造铝合金的工艺优化和性能预测提供了可靠的计算工具。
