激光粉末床熔融（L-PBF）增材制造中通过全零件热历史控制实现Alloy 718性能优化

《Additive Manufacturing》：Property optimization through full-part thermal history control in laser powder bed fusion additive manufacturing

为解决L-PBF制造的Alloy 718需额外时效热处理致成本高、周期长的问题，研究者开展全零件热历史控制研究，通过动态调控激光功率和基板温度，在实现原位时效（in-situ aging）的同时，使零件平均硬度达439 HV且均匀性显著提升，为免后处理高性能构件制造提供新路径。

　　

在金属增材制造（Additive Manufacturing, AM）尤其是激光粉末床熔融（Laser Powder Bed Fusion, L-PBF）技术飞速发展的今天，航空航天与能源领域对高性能金属构件的需求日益增加。Alloy 718（又称Inconel 718）作为一种沉淀硬化型镍基高温合金，因其优异的高温强度、抗蠕变和耐腐蚀性能，成为L-PBF制造关键部件的首选材料之一。然而，传统L-PBF工艺流程存在一个显著的痛点：打印出的零件（as-fabricated）通常需要进行耗时的后处理热处理（如固溶和时效），才能析出足够的γ′（Ni₃(Al, Ti)）和γ″（Ni₃Nb）强化相以达到所需硬度与强度。这不仅延长了交付周期，增加了能耗与成本，且热处理过程中的热循环可能引发晶粒粗化、变形甚至开裂，削弱了L-PBF技术本身具备的“近净成形”优势。

有没有可能让“热处理”在打印过程中就完成？事实上，L-PBF工艺中层层沉积带来的热循环（thermal cycling）本身就能使零件经历温度变化。如果能精准控制这一热历史（thermal history），使零件在制造过程中就达到理想的时效温度区间，便可实现“原位时效（in-situ aging）”，从而省去后续热处理步骤。但这并非易事，因为零件几何形状（如带有悬垂结构的部件）会导致热量积累（heat accumulation）不均，使得不同部位的热历史截然不同，最终导致零件性能（如硬度）空间分布极不均匀。为此，来自卡内基梅隆大学（Carnegie Mellon University）等机构的研究人员提出了一种创新的“全零件热历史控制”框架，通过动态调控工艺参数，在L-PBF过程中直接优化Alloy 718零件的硬度分布，该研究成果发表在《Additive Manufacturing》期刊上。

为达成上述目标，研究者综合运用了多项关键技术方法：首先，利用Gleeble热模拟试验机开展等温时效实验，结合显微硬度测试绘制时间-温度-硬度图谱，并拟合改进的Johnson-Mehl-Avrami-Kolmogorov (JMAK) 模型以描述非等温硬化动力学；其次，建立基于轴对称集总层（axisymmetric lumped-layer）的有限体积热传导模型，以较低计算成本模拟打印全过程的温度场演变；最后，将热模型与硬化动力学模型耦合，采用增广拉格朗日微分动态规划（Augmented Lagrangian Differential Dynamic Programming）算法，以目标硬度分布为约束，逆向求解随时间变化的激光功率与基板温度（baseplate temperature）轨迹，并在Aconity MIDI及EOS M 290设备上进行实验验证。

1. Introduction（引言）

研究指出传统L-PBF流程多采用恒定工艺参数，导致几何依赖的热积累引发性能不均。通过动态调节激光功率、扫描速度、基板温度等参数来控制最终性能（Property-driven process planning），是解决这一逆问题的关键。本文首次实现了针对近净形零件（near-net shape part）全区域硬度的控制。

2. Materials & methods（材料与方法）

2.1. Alloy 718 powder：实验采用VDM Metals与EOS两款Alloy 718粉末，成分符合标准，使用前经筛分与烘烤处理。

2.2. Selection of axisymmetric geometry for optimization：选用倒圆锥（inverted cone）作为研究对象，因其轴对称特性可简化热模型为二维计算，且其悬垂结构天然导致显著的热积累与硬度梯度，便于验证优化效果。

2.3. Fabrication using constant process parameters：采用恒定参数（激光功率285 W，基板温度150°C等）打印倒圆锥与倒金字塔，作为对照组。

2.4. Sample preparation, hardness testing, and microscopy：沿零件中心线剖分样品，采用维氏显微硬度计（1 kg力）测绘硬度分布，并利用SEM及STEM-EDS分析析出相。

2.5. Generation of time-temperature-hardness maps：利用Gleeble进行不同温度（450-800°C）与时间（0.05-5 h）的等温处理，建立硬度数据集，并通过有理多项式校准温度分布。

2.6. Calibration of hardening kinetics model：将硬度转换为等效析出相体积分数（volume fraction transformed），利用迭代拟合标定非等温JMAK模型参数（预指数因子A，激活能E_a，Avrami指数n）。

2.7. Thermal dynamics integrated with hardening kinetics：构建轴对称有限体积热模型，计算域离散为1.25 mm × 1.25 mm的体素（voxels），区分实体与粉末热物性，耦合JMAK模型计算每一步的硬化状态。

2.8. State-space representation for optimization and numerical solution：将热-硬化耦合系统表述为状态空间方程，采用矩阵指数法求解，并利用优化算法求解使目标区域硬度均匀分布于450 HV的时变激光功率与基板温度轨迹。

2.9. Fabrication using optimal parameters：依据优化轨迹，将倒圆锥切分为不同区段分别赋予对应激光功率，配合分段变化的基板温度，完成优化样品的打印。

3. Results & discussion（结果与讨论）

3.1. Initial observation and validation of in-situ precipitation hardening：恒定参数打印的倒金字塔零件因热积累，底部高温区硬度较顶部低混区提升约40%，STEM-EDS证实了γ′相的析出，验证了L-PBF过程中发生原位时效的可能性。

3.2. Characterization of hardening kinetics：绘制的等温时效曲线显示，峰值硬度约475 HV出现在725°C保温5 h左右；高于此温度因δ相（Ni₃Nb）析出导致过时效软化。显微组织观察揭示了枝晶间成分偏析对γ′/γ″析出及δ相转变的影响。

3.3. Observation of precipitation within compositionally inhomogeneous microstructures：L-PBF快淬组织存在胞/枝晶偏析，Nb富集的胞壁促进γ″早期析出，但也会加速有害δ相形成，单一JMAK模型虽简化了这些微观差异，但足以捕捉整体硬化趋势以指导工艺控制。

3.4. Prediction of in-situ hardening during fabrication：拟合的JMAK参数（n≈0.61，接近位错主导析出的理论值0.67）能准确预测恒定参数下倒圆锥零件的硬度分布（模拟vs实验误差在置信区间内）。

3.5. Successful optimization of parameter trajectories to maximize hardness：

• •
  3.5.1. Model validation using constant parameter experiments：恒定参数下倒圆锥平均硬度374 ± 41.3 HV，模拟值356 ± 46 HV，模型可靠性得到验证。
• •
  3.5.2. Calculation and experimental validation of optimal parameter trajectories：优化后轨迹显示，在热积累剧烈区降低基板温度以散热，在需更多热输入区提高基板温度或激光功率。实验打印零件平均硬度达439 ± 29 HV（目标450 HV），均方根误差33 HV，且在8%内可重复，虽总打印时间因缓冷增加约3 h，但远小于传统直接时效处理的18 h。
• •
  3.5.3. Quality and accuracy of simulated optimal solution：硬度未完全均匀源于几何体素数量权重差异及当前控制变量局限（层内功率不可变），未来可通过层内调制或改进成本函数解决。
• •
  3.5.4. Implications of simplified thermal and hardening dynamics：集总层轴对称模型虽忽略熔池细节与非轴对称特征，但在约束范围内计算高效，未来可朝三维自适应网格与非线性存储优化拓展。
• •
  3.5.5. Optimal solution dependence on initial input trajectories：基于物理意义的初始轨迹（介于未硬化与全硬化之间）有助于梯度优化算法逃离平坦区并收敛至最优解。
  3.6. Observation of strain-age cracking and oxidation in high temperature regions：恒定参数高热区因硬化与残余应力出现应变时效开裂（strain-age cracking），而优化样品因整体维持较高温度降低了热梯度与应力（偏转量减少约20%），避免了开裂，但表面氧化变色加重，需在惰性气体严格管控下权衡。

4. Conclusions（结论）

本研究成功演示了一套集成实验与建模的优化框架，通过全零件热历史控制指导L-PBF制造Alloy 718的硬度分布。主要结论包括：通过Gleeble实验标定了对L-PBF条件适用的JMAK硬化动力学模型；建立了结合集总层热模型与硬化动力学的预测工具，可毫米级分辨率模拟全零件硬度；利用数值优化求解出时变激光功率与基板温度轨迹，使倒圆锥零件在无需后时效热处理下获得439 ± 29 HV的均匀硬度，重复性在8%以内。这一将热历史、材料动力学与数值优化结合的“属性驱动”框架，不仅证明了在打印过程中集成热处理的可行性，也为未来在增材制造中直接控制其它基于热历史的材料属性（如残余应力、晶粒尺寸）提供了基础，有望在资源受限或需快速响应的场景下实现高性能金属构件的免后处理制造。