4. 材料本构与相变模型

在增材制造热历史与变形仿真中，热源模型决定“输入多少热”，而材料本构决定“材料如何响应热”。同一套温度场，若本构选择不同，残余应力与翘曲的预测结果可能相差数倍。本节围绕三条主线展开：常用本构的适用边界、固态相变对热应变与力学响应的耦合影响，以及材料参数的获取途径。

一、三类本构的物理含义与适用边界

金属增材制造中，材料在熔池附近经历熔化—凝固—快速冷却，在远离熔池区域则经历多次热循环。不同位置、不同时刻的力学行为差异极大，因此需要按温度区间和加载速率选择本构。

弹塑性本构是最基础的一类。它假设总应变由弹性应变与塑性应变叠加，屈服后应力由硬化曲线决定。常用形式为等向硬化加随动硬化组合，以描述循环加载下的包辛格效应。在增材制造中，弹塑性模型适合描述固态、温度较低、加载速率不敏感的区域，例如基板附近和已沉积层的低温段。其关键参数包括弹性模量、泊松比、屈服强度、硬化模量，且这些参数均随温度变化。

粘塑性本构引入应变率相关性，屈服面随等效塑性应变率变化。典型形式为Perzyna型或Chaboche粘塑性模型。在激光粉末床熔融中，熔池附近冷却速率可达10⁴～10⁶ K/s，材料在高温下表现出明显的率敏感性，粘塑性模型能更合理地描述快速冷却过程中的应力松弛与塑性流动。其参数除弹塑性参数外，还需粘性参数、率敏感指数等。

蠕变本构描述恒定应力下应变随时间累积的行为，常用Norton-Bailey幂律形式。在增材制造中，蠕变主要影响高温停留阶段，例如预热基板、层间保温或电子束熔融的高温成形过程。若忽略蠕变，长时间高温区的应力预测会偏高，导致残余应力与变形评估失真。

二、固态相变对热应变与力学响应的影响

固态相变是增材制造仿真中容易被忽略却影响显著的因素。以Ti-6Al-4V为例，冷却过程中β相向α相转变，伴随体积变化和相变塑性。相变对仿真的影响体现在三个方面。

第一，热应变不再仅由热膨胀系数决定。相变引起的体积变化会叠加到热应变上。若仅用单一热膨胀系数，冷却段的应变预测会偏离实际。正确做法是将热应变分解为热膨胀项与相变项，相变项由相分数演化决定。

第二，相变塑性。在相变过程中，即使应力低于屈服强度，材料也会产生不可逆变形，称为相变塑性。其大小与相变体积变化、应力状态和相变进度相关。对于Ti-6Al-4V和部分钢种，相变塑性对最终残余应力的贡献不可忽略。

第三，力学性能随相分数变化。弹性模量、屈服强度等参数是各相性能的混合结果。若采用单一相参数，高温段的应力响应会失真。常用混合法则为线性混合或Voigt-Reuss混合，需根据材料体系选择。

三、常用金属材料参数获取途径

本构模型的精度最终取决于参数质量。增材制造常用金属包括Ti-6Al-4V、316L不锈钢、IN718镍基合金、AlSi10Mg铝合金等。参数获取通常有四条途径，各有适用场景。

文献与手册是最快捷的起点。ASM手册、材料供应商数据表可提供室温至中温段的弹性模量、屈服强度、热膨胀系数。但需注意：增材制造材料的微观组织与铸锻材不同，直接套用铸锻参数会引入系统误差。建议优先选用同材料、同工艺的文献数据。

热力学计算软件如JMatPro、Thermo-Calc可计算平衡或亚稳相图、相分数随温度变化、热膨胀系数、弹性模量等。对于相变耦合仿真，这类工具能提供相分数演化曲线，是连接相变与本构的关键输入。但计算结果是基于热力学平衡假设，对快速冷却的非平衡过程需谨慎修正。

热模拟实验是获取高温力学参数最直接的方法。Gleeble热模拟试验机可对试样施加与实际增材制造相近的热循环，同时测量应力—应变响应。通过不同温度、不同应变率下的压缩或拉伸试验，可拟合弹塑性、粘塑性参数。蠕变参数则需在恒定温度和载荷下进行长时间试验。这类实验成本较高，但数据最贴近实际工艺。

反演标定是当实验数据不足时的补充手段。以实测的残余应力或变形为目标，通过优化算法调整本构参数，使仿真结果逼近实测值。常用方法包括响应面法、遗传算法、贝叶斯标定等。反演标定需注意参数可辨识性，避免多参数同时拟合导致非唯一解。

四、实操中的参数组织与验证

在实际仿真项目中，建议将材料参数组织为温度相关的表格形式，每个温度点对应一组完整参数。以Ti-6Al-4V为例，可设置20℃、200℃、400℃、600℃、800℃、1000℃等温度节点，中间温度线性插值。相变参数单独定义相分数随温度的变化曲线，并指定相变开始与结束温度。

# 示例：Ti-6Al-4V 弹塑性参数表（示例数据，非实测）
# T(℃)   E(GPa)   ν      σy(MPa)   H(GPa)
  20      110     0.34    880       12
  200     95      0.34    720       10
  400     78      0.35    550       8
  600     58      0.36    380       5
  800     35      0.38    180       2
  1000    15      0.40    50        0.5

# 相变参数示例
# β→α 相变开始温度：980℃
# β→α 相变结束温度：800℃
# 相变体积变化：+0.5%（示例值）
# 相变塑性系数：5e-5 /MPa（示例值）

参数设置完成后，需通过简单验证算例检查合理性。例如单道单层沉积的基板翘曲量、冷却后残余应力分布趋势。若翘曲方向与量级明显偏离文献或经验，应优先检查热膨胀系数、屈服强度随温度的变化是否合理，以及相变参数是否被正确激活。验证通过后，再将本构模型应用于多道多层的大规模仿真。

综上，本构模型的选择与参数获取是增材制造热力耦合仿真的核心环节。弹塑性、粘塑性、蠕变三类本构各有适用区间，固态相变则通过相变应变、相变塑性和相分数混合性能三条路径影响力学响应。参数获取应优先采用同工艺实验数据，辅以热力学计算和反演标定，并通过简单算例验证后再投入大规模仿真。