4. 增材制造材料热物理与力学性能

在增材制造中,残余应力并不是单纯由“热胀冷缩”造成的,而是材料在快速加热、快速冷却、反复热循环下,其热物理性能与力学性能随温度剧烈变化后共同作用的结果。理解这一节的关键,是把材料参数看作“温度的函数”,而不是一组常温常数。只有把热导率、比热、密度、热膨胀系数、弹性模量、屈服强度、相变潜热、蠕变参数等放到同一温度轴上观察,才能解释为什么同一材料在不同工艺参数下会呈现完全不同的应力与变形行为。

一、温度相关性能:决定温度场与应力场的基础

热物理性能决定热量如何进入、扩散和散失。热导率 k(T) 越高,熔池热量越容易向周围传导,温度梯度相对平缓,残余应力倾向降低;热导率低时,热量集中在熔池附近,梯度陡峭,应力集中更明显。比热 c_p(T) 与密度 ρ(T) 共同决定体积热容 ρc_p,它影响升温与降温速率。以常见金属为例,钛合金 TC4 在室温附近热导率约 7 W/(m·K),而铝合金 AlSi10Mg 约 150 W/(m·K),两者在激光粉末床熔融中的温度场形态差异显著,这也是铝件往往变形更分散、钛件局部应力更集中的原因之一。

力学性能同样随温度变化。弹性模量 E(T) 随温度升高而下降,屈服强度 σ_y(T) 在高温下显著降低。这意味着在熔池附近的高温区,材料几乎不能承担弹性应力,而远离熔池的低温区仍保持较高刚度,约束了高温区的收缩,从而在冷却后形成“近缝区受压、远区受拉”的典型应力分布。若把 E 和 σ_y 当作常数计算,往往会高估高温区应力、低估塑性应变积累。

二、相变:改变体积与应力路径的“隐藏变量”

相变对残余应力的影响常被低估。以钛合金为例,β→α 相变伴随体积变化和潜热释放,会改变局部冷却曲线;在钢类材料中,奥氏体向马氏体转变时体积膨胀,可在冷却后期引入压应力,甚至部分抵消热收缩拉应力。相变还使热膨胀系数在相变区间出现非线性,若仍用单一 α 值计算,热应变会明显失真。工程上处理这类问题,通常需要在材料模型中输入随温度变化的 α(T)、相变起始/结束温度、潜热以及相变塑性参数,而不是只给一个平均线膨胀系数。

三、蠕变与塑性行为:应力松弛与不可逆变形的来源

增材制造的热循环具有“高温停留时间短、反复多次”的特点。在高温阶段,即使时间很短,蠕变也会造成应力松弛,使部分弹性应变转为非弹性应变。蠕变通常用 Norton 律描述:ε̇_cr = A σ^n exp(-Q/RT),其中 A、n、Q 为材料参数。温度越高、应力越大,蠕变速率越快。塑性行为则决定材料何时屈服、如何硬化。若采用各向同性硬化模型,需给出不同温度下的屈服面与硬化曲线;若考虑循环加载,还应引入随动硬化以描述包辛格效应。二者共同决定冷却后残余应力的大小与方向。

温度相关材料性能 → 热-力响应 → 残余应力与变形 热物理性能 k(T), cp(T), ρ(T) • 热导率 k(T):决定温度梯度陡缓 • 体积热容 ρcp:决定升降温速率 • 热膨胀系数 α(T):热应变来源 • 相变潜热 L:改变冷却曲线 • 密度 ρ(T):影响熔池流动与传热 示例:TC4 k≈7,AlSi10Mg k≈150 W/(m·K) 力学性能 E(T), σy(T), 硬化 • 弹性模量 E(T):高温下降 • 屈服强度 σy(T):高温显著降低 • 硬化曲线:决定塑性应变积累 • 包辛格效应:循环加载需随动硬化 • 蠕变参数 A, n, Q:高温应力松弛 示例:高温区 E 下降 40%~60% 相变与微观组织演化 • β→α 相变:体积变化 + 潜热 • 奥氏体→马氏体:体积膨胀 • 相变塑性:应力下相变附加应变 • α(T) 非线性:不能取平均值 • 相变温度区间:影响应力反转 示例:马氏体相变可引入压应力 热-力耦合本构响应 热应变 + 弹性应变 + 塑性应变 + 蠕变应变 + 相变应变 残余应力分布与变形结果 近熔池区:高温低屈服 → 冷却后受压;远区:约束收缩 → 受拉 宏观表现:翘曲、开裂、层间剥离、尺寸超差 温度轴上的性能变化 室温 熔点 E(T) 下降 σy(T) 下降 α(T) 非线性

四、参数如何共同影响残余应力:一个可操作的判断链

在实际分析中,可以按以下顺序检查材料参数是否“用对”:第一,确认热导率与比热是否随温度输入,若只给常温值,温度场会失真;第二,确认热膨胀系数是否为瞬时值或割线值,并明确参考温度;第三,确认屈服强度与弹性模量是否覆盖从室温到熔点的范围;第四,确认是否考虑相变潜热与相变塑性;第五,确认高温蠕变是否被激活。以激光粉末床熔融 TC4 为例,若忽略蠕变,计算得到的残余应力往往偏高;若同时忽略 β→α 相变体积变化,应力方向甚至可能判断错误。

教学示例:在热-力耦合仿真中,可将材料参数按温度点离散输入,例如每 50 ℃ 一个数据点,并在相变区间加密到每 10 ℃。这样既能控制计算量,又能避免相变区间的非线性被平滑掉。

五、面向工程的简化处理建议

对于只需要判断变形趋势的快速评估,可以采用“等效温度相关”策略:保留 E(T)、σ_y(T)、α(T) 三条主曲线,将蠕变等效为高温区屈服强度折减,将相变等效为相变区间的附加体积应变。对于需要预测开裂或尺寸精度的场景,则必须引入蠕变与相变塑性,否则近缝区应力状态和塑性区宽度都会偏差较大。无论采用哪种简化,都应记录参考温度、坐标系与符号约定:本文约定拉应力为正、压应力为负,热膨胀以无应力参考温度为零点。

材料热物理与力学性能是增材制造残余应力分析的“输入地基”。温度相关性能决定温度场与热应变,相变改变体积与应力路径,蠕变与塑性决定应力能否松弛以及不可逆变形如何积累。把这些参数按温度轴组织起来,并与工艺热循环对应,才能从材料层面解释并控制残余应力与变形。