3. 温度场分布与热循环：热传导方程与边界条件、瞬态温度场模拟、热循环曲线特征分析
各位工程师朋友，咱们接着聊。上一节我们把热源模型和熔池行为掰开揉碎讲了一遍。这一节，咱们得把视野拉大一点——从熔池这个小局部，扩展到整个零件的温度场分布。
说白了，电子束增材制造的过程，就是一个反复加热、冷却的过程。你想想看，电子束在粉末床上扫来扫去，每一层、每一条轨迹，都在给零件注入热量。这些热量怎么传？怎么散？零件内部温度怎么变？这些搞不清楚，变形控制就是一句空话。
我个人习惯，做工艺仿真之前，先把热传导方程和边界条件理清楚。这是基本功，也是后面所有分析的地基。
核心观点：温度场分布决定了热应力分布，热应力分布决定了最终变形。控制变形，本质上就是控制温度场。
3.1 热传导方程与边界条件
先看热传导方程。三维瞬态热传导的控制方程长这样：
ρ·c·(∂T/∂t) = k·(∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z²) + Q其中：
- ρ — 材料密度 (kg/m³)
- c — 比热容 (J/(kg·K))
- k — 导热系数 (W/(m·K))
- T — 温度 (K)
- t — 时间 (s)
- Q — 内热源项 (W/m³)，也就是电子束输入的能量
这个方程看着简单，但实际用起来，坑不少。我在项目中遇到过最头疼的问题，就是材料参数随温度变化。你想想看，从室温到接近熔点，钛合金的导热系数能差好几倍。如果还用常数去算，那结果基本就是错的。
边界条件这块，我一般分三类来处理：
| 边界类型 | 数学表达 | 物理含义 |
|---|---|---|
| 第一类（Dirichlet） | T = T₀ | 给定边界温度，比如基板预热温度 |
| 第二类（Neumann） | -k·(∂T/∂n) = q | 给定热流密度，比如电子束热源输入 |
| 第三类（Robin） | -k·(∂T/∂n) = h·(T - T∞) | 对流换热，比如粉末床表面与腔室气体 |
嗯，这里要注意。电子束增材制造有个特殊之处——真空环境。真空下没有对流换热，辐射换热就成了主要散热方式。所以边界条件里，还得加上辐射项：
q_rad = ε·σ·(T⁴ - T∞⁴)ε是发射率，σ是斯特藩-玻尔兹曼常数。这个辐射项在高温区特别重要，温度越高，辐射散热占比越大。
个人经验：我建议在仿真中把辐射和对流分开处理。有些软件把两者合并成一个等效换热系数，虽然省事，但高温区误差很大。尤其是做大型零件仿真时，这个误差会被放大。
3.2 瞬态温度场模拟
搞清楚了方程和边界条件，接下来就是怎么算的问题。瞬态温度场模拟，说白了就是求解上面那个偏微分方程。但实际做起来，有几个关键点要把握好。
第一，时间步长的选择。电子束扫描速度很快，通常每秒几百毫米。如果时间步长太大，会漏掉熔池附近的温度峰值。如果太小，计算量又受不了。我个人的经验是：扫描阶段用0.1-0.5ms的时间步长，冷却阶段可以放大到0.1-1s。
第二，网格尺寸的权衡。熔池附近温度梯度极大，网格要细。远离熔池的地方，温度变化平缓，网格可以粗。我习惯用过渡网格，从熔池区域的0.1mm逐渐过渡到零件边缘的2-3mm。这样既保证精度，又不至于算到天荒地老。
第三，生死单元技术。这个在增材制造仿真里特别重要。粉末层的初始状态是松散的，导热系数很低。电子束扫过之后，粉末熔化凝固变成实体，导热系数瞬间提升。仿真中需要用生死单元来模拟这个过程——粉末层初始是"死"的，被扫到后才"活"过来。
下面是一个简单的瞬态温度场模拟流程：
1. 初始化：设置初始温度（通常为预热温度）
2. 激活当前层粉末单元
3. 施加电子束热源（移动热源模型）
4. 求解热传导方程，得到当前时间步的温度场
5. 更新材料参数（随温度变化）
6. 判断是否扫描完毕
- 否：移动到下一位置，重复步骤3-5
- 是：进入冷却阶段
7. 冷却阶段：移除热源，仅考虑散热边界
8. 输出温度场结果这个流程看着简单，但实际跑起来，一个大型零件可能要算好几天。我做过一个300mm长的钛合金零件，网格数量超过500万，在32核工作站上跑了整整72小时。所以，仿真策略很重要，不是所有细节都要算到极致。
避坑指南：我曾经犯过一个低级错误——忽略了粉末床的导热各向异性。粉末层在水平方向和垂直方向的导热系数差别很大，垂直方向因为颗粒接触点少，导热更差。如果当成各向同性处理，温度场会严重失真。
3.3 热循环曲线特征分析
温度场模拟做完，我们得到的是整个零件在任意时刻的温度分布。但光看云图还不够，我们需要提取关键点的温度随时间变化曲线——这就是热循环曲线。
热循环曲线长什么样？我画个示意图：
这张图里，我标出了几个关键特征：
- 峰值温度 Tₚ — 电子束扫过该点时达到的最高温度，通常接近甚至超过材料熔点
- 升温速率 — 从初始温度到峰值温度的斜率，电子束工艺的升温速率可达10⁵-10⁶ K/s
- 降温速率 — 峰值温度后的冷却斜率，决定了凝固组织和残余应力
- 热平衡温度 — 多次扫描后，零件整体温度趋于稳定的平台值
为什么热循环曲线这么重要？因为它直接决定了材料的微观组织和残余应力状态。
举个例子。我在做TC4钛合金零件时发现，如果热循环的峰值温度超过β相变温度（约995°C），冷却速率又足够快，就会形成细小的马氏体α'相。这种组织强度高但塑性差，容易开裂。反过来，如果峰值温度低于β相变温度，或者冷却速率慢，得到的是α+β双相组织，综合性能更好。
所以，控制热循环曲线，本质上就是在控制材料的最终性能。
3.4 多层多道热循环的叠加效应
单层单道的热循环好理解，但实际零件是几十层、几百层堆叠起来的。每一层扫描时，都会对已经凝固的下层产生一次"回火"效应。这种叠加效应，我称之为"热历史积累"。
热历史积累有几个典型特征：
- 峰值温度逐层降低 — 远离当前扫描层的区域，受热影响越来越小
- 热平衡温度逐渐升高 — 随着层数增加，零件整体温度不断爬升
- 冷却速率逐渐减慢 — 零件温度高了，散热驱动力变小了
我记得有一次做仿真，一个200层的零件，前50层热循环曲线差异很大，但到了100层以后，曲线基本稳定了。这说明零件进入了"准稳态"热循环阶段。这个阶段，每层的热条件几乎一样，工艺参数可以保持稳定。
关键结论：大尺寸零件的变形控制，重点在前几十层。这个阶段温度场变化剧烈，热应力积累最快。一旦进入准稳态，变形速率会明显放缓。
3.5 温度场仿真的实用建议
最后，我把自己这些年做温度场仿真的几点心得分享给大家：
| 问题 | 我的做法 | 原因 |
|---|---|---|
| 材料参数缺失 | 用JMatPro或CALPHAD计算 | 实测数据往往不全，计算数据比猜测靠谱 |
| 计算时间太长 | 采用自适应时间步长 | 加热阶段用小步长，冷却阶段用大步长 |
| 结果不收敛 | 检查网格质量，避免畸变单元 | 六面体网格比四面体更稳定 |
| 与实验对不上 | 校准热源参数和换热系数 | 仿真参数需要实验标定，不能照搬文献 |
嗯，这一节内容不少。热传导方程是理论基础，瞬态模拟是实践工具，热循环曲线是分析手段。三者环环相扣，缺一不可。你把这些搞透了，温度场这块就算拿下了。
一个小技巧：做温度场仿真时，我习惯先在局部区域做网格敏感性分析。用0.2mm、0.1mm、0.05mm三种网格尺寸分别算一个简单算例，看温度结果差异。如果0.1mm和0.05mm的结果相差不到5%，那就用0.1mm。这样既保证精度，又不浪费计算资源。
公众号:蓝海资料掘金营，微信deep3321