小编亲测 | 复现Li & Liou ：DED激光熔覆

论文的研究对象很具体：在Ti6Al4V钛合金基板上加工两种V形槽——一种是侧壁45°的三角形槽，一种是侧壁90°的矩形槽。然后用DED激光工艺回填，对比修复质量。实验发现，三角形槽修复后界面结合良好，拉伸试样断裂在沉积区而非界面；矩形槽却因为激光打不到垂直侧壁，底部出现了大孔隙，直接判废，小编复现的是V形槽，文末有资源如下。

一、论文原文在讲什么

论文的实验设置很清晰：基板尺寸60×19×10 mm，V形槽位于长度中点，深3mm。激光功率650W，扫描速度210 mm/min，层厚0.5mm，光斑直径1mm。三角形槽采用1+2+3+4+5+6的道次分布，共21道；矩形槽每层3道，共18道。

图1 论文Fig.1：两种V形槽修复几何与扫描策略对比

实验结果很直观：三角形槽修复后的横截面（论文Fig.3c）显示沉积材料与基板结合紧密，没有分层或孔隙；矩形槽（Fig.3d）却在侧壁底部出现了明显的大孔隙——激光束无法有效照射到90°垂直侧壁，粉末熔化了却没能跟基板熔合在一起。拉伸测试也印证了这一点：三角形槽修复件的平均抗拉强度1037 MPa，试样断裂在沉积区，说明界面强度比沉积区还高。

仿真部分，论文用ANSYS Mechanical APDL建立了顺序耦合热-力模型。先做瞬态温度场，再把温度场作为体载荷做弹塑性应力分析。沉积单元沿扫描方向依次激活（单元生死法），用基板侧表面的热电偶TC1标定模型。温度场验证通过后，再预测残余应力和变形。

二、COMSOL复现现场：从ANSYS到COMSOL的转译

论文用的是ANSYS APDL，小编把它转到了COMSOL 6.4。两者在原理上完全对应：ANSYS的"单元生死"对应COMSOL的Quiet/Active属性缩放；ANSYS的均匀体热源对应COMSOL的域热源节点；ANSYS的温度相关对流系数表，对应COMSOL的插值函数。

全局参数长这样（可直接抄）：

L_sub=60[mm], W_sub=19[mm], t_sub=10[mm]
d_groove=3[mm], x_groove=30[mm]
P_las=650[W], eta_abs=0.3, vs=210[mm/min]
r_las=0.5[mm], t_lay=0.5[mm], w_trk=1[mm]
N_trk=21, t_trk=W_sub/vs=5.4286[s]
T_amb=25[degC], eps_r=0.3, h_clamp=100[W/(m^2*K)]
alpha_q=1e-4, tsm=0.03[s]

几何建模是这道题的难点。基板是一个长方体，中间用布尔差集挖掉V形槽。然后槽内要建21个独立的沉积道棱柱——每个道一个独立的域，这样才能逐道激活。三角形槽的截面是倒三角，顶点在(30, 7)，顶边宽6mm。第L层有L条道，每条道宽1mm、高0.5mm，最外侧道的下顶点自动落在45°槽壁上。

热源表达式是论文式(3)的COMSOL转写——均匀体热源：

Q0 = eta_abs*P_las/(pi*r_las^2*t_lay) = 4.97e11 [W/m^3]
Q_src = Q0 * (k_cur>=1) * flc2hs((r_las-r_dist)/dsm,1)
      * flc2hs((z-zlow)/dsm,1) * flc2hs((zlow+t_lay-z)/dsm,1)

三个flc2hs因子分别是：沉积期开关、激光圆盘掩膜、当前层z区间窗。任一位置只有同时满足"属于当前道、在激光圆盘内、在当前层高度"时才获得热源。

激活策略跟之前Nain论文类似，但这里是逐道激活而非逐层激活。21个沉积域各有一个激活变量actk，在激光开始扫描第k道的时刻从alpha_q·k_int平滑升到k_int。只缩放热导率k，不缩放比热Cp和密度ρ——未激活单元保留热容，只切断导热，避免激活瞬间能量凭空注入造成温度跳变。

边界条件方面，论文把对流和辐射合并成一个用户定义热通量。COMSOL里如果分别用"热通量(对流)"+"表面对环境辐射"，辐射节点的T⁴公式在负温度下会反转抽热，导致发散。所以手册里写了一个保护版通量：

q_prot = h_conv(T)*(T_amb-T) + eps_r*5.67e-8*((T_amb+273.15)^4 - max(T,293.15)^4)

关键技巧：辐射项自变量用max(T,293.15)钳位，防止T<0时T⁴反号。全部插值函数外推方式设"常数"，避免低温区线性外推出负热导率。求解器只用PARDISO直接求解器，关闭迭代求解器——这是防止负温度链式发散的三件套。

三、真实跑出来的结果

下面三张图是实际复现的截图，不是示意图，是COMSOL结果窗口里直接截的。

图2 TC1热历史（基板侧表面，距槽中心4mm）—— 真实复现截图

TC1温度曲线长这样：从25°C起步，随着21道激光依次扫过，出现了21个清晰的升温峰。峰值温度从第一道约850°C逐渐爬升，中间几道因热积累效应峰值更高，最后几道稳定在800°C左右。114秒沉积结束后进入冷却段，温度缓慢下降。这跟论文Fig.8b的趋势完全对得上——仿真红线与实验黑线几乎重合，说明热源模型和边界条件标定是靠谱的。

图3 模型树、热源体积分验证与几何预览 —— 真实复现截图

这张图展示了模型树结构和热源验证。在"体积分1"节点里，对沉积域的Q_src做体积分，t=103s时结果为195.16 W。激光功率650W，吸收率0.3，有效输入功率就是195W——体热源积分值跟理论值对上了，说明热源表达式没有量纲错误，几何选择也没漏域。右边是三维几何预览，基板、V形槽和21条沉积道一目了然。

图4 t=90s时的三维温度场 —— 真实复现截图

90秒时，第16道左右正在沉积。温度场显示激光所在区域温度最高（红黄色，超过2000°C），热量向周围基板和已沉积层扩散。可以看到V形槽底部已经被前几道填满，当前激光在槽的中上部移动。论文Fig.7显示第1道末端峰值1854°C，第16道2029°C，第21道2284°C——的复现结果在同一量级，热积累趋势一致。

图5 论文Fig.8：TC1热历史仿真与实验对比

四、复现这篇论文，到底能用在哪？

这篇论文的仿真模型跑通了，但真正的价值不止于"复现"。我琢磨了一下，至少有四个方向可以直接落地：

第一，修复槽形优化设计。论文最核心的发现是：45°三角形槽能修，90°矩形槽修不了。但45°就是最优吗？如果换成30°或60°呢？仿真模型可以帮你扫一遍侧壁角度-修复质量的参数空间，找到既能保证侧壁熔合、又尽量少切削基板材料的"甜点角度"。这比每次上机试切省下的时间和材料费，够跑几百次仿真了。

第二，残余应力与变形预测。论文的顺序耦合热-力模型预测了：三角形槽修复后基板顶面最大纵向拉应力1525 MPa，最大位移0.08 mm；矩形槽位移更大（0.122 mm），且局部等效应力超过屈服强度，可能诱发裂纹。这些数字直接决定了修复件能不能直接装机使用，还是需要后续去应力退火或矫形。

第三，工艺参数快速标定。论文用TC1热电偶标定了模型，确认吸收率0.3是合理的。如果你的实验换了激光器（比如从YAG换成光纤激光）、换了粉末粒度、或者换了保护气成分，吸收率会变，但标定方法不变——先假设一个值，跑仿真，对比热电偶曲线，迭代到吻合。这个闭环流程是仿真工程师的看家本领。

第四，当作COMSOL进阶教材。这个模型比Nain的薄壁沉积更复杂：V形槽几何、21个独立域的逐道激活、zigzag扫描策略、温度相关的对流系数、负温度防护三件套。如果你已经跑通了薄壁沉积，拿这个模型进阶，能学到更多工程化的调试技巧。

说到底，仿真不是目的，是手段。论文里的温度场告诉你"激光能不能熔到侧壁"，应力场告诉你"修完会不会变形开裂"。这两个问题搞清楚了，修复工艺才算真正从"经验试错"走向了"定量设计"。

写在最后：三个提醒

如果你也想复现这篇论文，有三件事我踩过坑，提前告诉你：

第一，负温度防护三件套缺一不可。这篇论文的对流系数h(T)在25°C时只有12 W/(m²·K)，如果仿真中某个点温度下冲到低值，内置辐射节点的T⁴公式会反号，把该点当成"热源"而不是"散热"，温度越算越低，最终发散。解决办法：禁用内置辐射节点，改用钳位通量q_prot；插值函数外推设"常数"；只用PARDISO直接求解器。这三件套缺一个，模型都可能跑崩。

第二，只缩放k、不缩放Cp和ρ。有些教程教你把Quiet/Active单元的所有物性都缩到1e-4倍，但这会导致激活瞬间能量凭空注入——未激活单元本来温度是25°C，一激活突然变成真实材料，比热和密度也瞬间跳变，温度场会出现诡异的尖峰。正确的做法是只缩导热率k，保留热容Cp和密度ρ，让未激活单元"存得住热、传不出去"。

第三，热源体积分必须自检。Q0 = eta_abs*P/(pi*r²*t_lay) 的量纲是W/m³，但你要确认这个体热源施加在正确的几何域上。在COMSOL里加一个"体积分"派生值，对沉积域的Q_src积分，结果应该约等于eta_abs*P（195W）。如果对不上，要么是表达式写错了，要么是几何选择漏了域——这是检验热源是否生效的最快办法。

论文上的公式再漂亮，也只有你自己在COMSOL里点下"计算"、看到21个温度峰值依次蹦出来的那一刻，才算真正读懂了。从"能跑"到"能跑对"，中间隔着一百个报错和一千次调试——但这就是仿真的日常。

几个思考

1. 这篇论文的热源是均匀体热源，不是双椭球。Li & Liou用的是Q = αP/(πr²t)，一个平顶圆柱体热源。这跟Nain论文的Goldak双椭球完全不同——双椭球前后不对称，能模拟熔池拖尾；均匀热源简单，但假设激光能量在光斑内均匀分布。如果你的激光是高斯分布，均匀热源会高估中心温度、低估边缘温度。你有没有想过，自己的激光到底是什么分布？

2. 单元生死法 vs Quiet/Active法，哪个更适合修复？论文用ANSYS的单元生死法，小编在COMSOL里用属性缩放法。生死法需要方程重编号，实现复杂但计算快；缩放法简单但自由度不会减少。对于21个域的修复模型，两种方法差别不大；但如果修复区域是复杂曲面、几百个域，缩放法的求解时间会不会拖垮你？

3. 热-力耦合的顺序耦合够吗？论文先做热分析，再把温度场作为载荷做力学分析，这是顺序耦合。但DED过程中材料经历熔化-凝固-相变，高温下的应力松弛和固态相变（Ti6Al4V的α→β转变）都没考虑进去。如果修复的是航空发动机叶片，残余应力预测差10%可能就是安全和不安全的区别。什么时候必须用完全耦合？顺序耦合的误差边界在哪？

4. 从薄壁沉积到V形槽修复，建模思路变了什么？Nain的薄壁是20层规整长方体，每层几何一样；Li的修复是21个不规则梯形/三角形道，每道尺寸不同。前者可以用Array阵列复 制，后者必须逐道画截面。这提醒：增材仿真的建模成本，很大程度上取决于几何复杂度。有没有可能用参数化几何（COMSOL的LiveLink for MATLAB或CAD导入）来自动生成修复路径？

这些问题没有标准答案。但正是这些问题，让一篇复现文章不只是"我跑出来了"，而是"我跑完之后，开始想下一步该怎么走"。

参考标准与文献
[1] Li L, Zhang X, Liou F. (2021). Experimental and Numerical Investigation in Directed Energy Deposition for Component Repair. Materials 2021, 14, 1409. doi:10.3390/ma14061409
[2] Zhang X, et al. (2019). Experimental characterization of a direct metal deposited cobalt-based alloy on tool steel for component repair. JOM 71, 946–955.
[3] Graf B, Gumenyuk A, Rethmeier M. (2012). Laser metal deposition as repair technology for stainless steel and titanium alloys. Phys. Procedia 39, 376–381.
[4] COMSOL Multiphysics 6.4 Documentation: Heat Transfer Module User's Guide.