☰
Flow-3D激光送粉熔敷仿真:从粒子追踪到熔池分析实战指南
2026/10/2 22:27:46 网站建设 项目流程

刚开始接触"Flow3D 激光送粉增材制造"这个方向的时候,我其实是有些发怵的。激光熔敷本身是高温物理过程,粉末在气流里飞、熔池在基板上淌,还要搭积木一样一层层堆出零件形状,单靠氢脆试错改工艺参数,日子太难过了。做仿真就是要提前"看见"粉末怎么落、熔池怎么长、粒子怎么铺,把看不见的高温过程搬到屏幕上反复试。我自己踩了无数坑之后,才慢慢把 Flow-3D 里激光送粉熔敷这一套流程跑通,这里把最核心的东西掰开揉碎讲一遍。

需要提醒的是,这里说的"FDM 激光熔敷"不是 3D 打印常说的 FDM 熔融沉积,而是激光熔敷(Laser Cladding / Laser Metal Deposition)方向。有些资料里会把 Laser Metal Deposition 缩写写成 LMD,也有写成 DMD 的,实际指的都是同轴送粉激光熔敷这个增材制造工艺。Flow-3D 在这类场景里主要解决的是"粉末粒子群在气流中的运动轨迹"和"熔池流体流动与凝固"两个问题。标题里提到的"粒子"追踪和"粒子子直..."这类操作,对应到软件里其实就是两类东西:一是离散粒子模型里的粉末颗粒注入,二是用无质量粒子做流场示踪。这两个功能是后期分析送粉利用率、熔池流型最重要的抓手。

1. 激光熔敷的物理本质:熔池、粉末流与那个"粒子"难题

很多第一次接触这个方向的朋友,会觉得激光熔敷仿真就是"给一块板加个热源,再把粉末撒上去"。真开工之后才发现,这个工艺里至少有三组物理场在同时打架:载气携粉的湍流射流场、激光辐照下的熔池流场(含固液相变和表面张力驱动对流)、粉末穿越激光束时自身被预热甚至部分熔化的过程。

1.1 粉末粒子为什么是核心问题

熔敷质量的好坏,很大程度不取决于激光给了多少能量,而是取决于粉末粒子到底落在哪、落进去多少。同轴送粉的喷嘴设计,通常让粉末流汇聚在激光光斑正下方的一个小焦面上。粉斑直径、焦点高度、送粉角度这三个量,决定了粉末是"精准落进熔池被完全熔化"还是"飘散在热影响区变成废料"。

我见过太多人一开始把注意力全放在激光功率上,结果粉末利用率低得吓人,基板上形成一圈一圈的未熔化颗粒。实际上,粉末颗粒一旦进入激光辐照区,会被加热到很高温度,部分颗粒在落到熔池之前就已经熔化甚至汽化。如果粉末流汇聚点高于熔池表面,粉末会在半空"开花";如果汇聚点太低,粒子又会被熔池上方的金属蒸汽和热气吹飞。仿真里如果不把粒子追踪做好,这些现象根本看不见。

1.2 Flow-3D 的求解思路与优势

Flow-3D 在激光送粉熔敷仿真上的核心逻辑,是用有限体积法求解纳维-斯托克斯方程,同时用 VOF 方法追踪空气与金属熔体之间的自由表面,再用 FAVOR 方法处理基板和喷嘴这类固体几何边界。

跟其他通用 CFD 软件相比,它的优势主要体现在三个地方:

  • 自由表面处理非常成熟。熔池表面不是平的,激光作用区中心凹陷、边缘隆起,还有金属蒸汽反冲压和表面张力综合塑造曲面。VOF 方法天然适合这类问题。
  • 自带完整的凝固/熔化模型,可以设置固相线、液相线、潜热,熔池凝固时糊状区处理比较稳定。
  • 粒子模块灵活。既可以在入口注入粉末粒子,追踪其在气体中的运动、与壁面的碰撞,也可以在流体域内注入无质量粒子来示踪熔池流动。

需要强调的是,Flow-3D 并不是全能的。它对温度场的求解在焊道搭接这种长时间多道熔敷场景里,计算量会非常大。实际工程上我会用"局部熔敷+周期性边界"的思路,先算单道,确认工艺窗口,再做多道拼接的趋势分析。

1.3 一个典型的同轴送粉熔敷问题长什么样

我们以最常见的平板单道熔敷为例。基板是低碳钢,或者钛合金,表面有一层同轴送粉喷嘴,激光沿 X 方向移动,粉末由 Ar 气载送。仿真要回答的问题很明确:给定激光功率 P、扫描速度 v、送粉量 m、喷嘴倾角与保护气流量,熔敷层宽度、高度、稀释率各是多少?粉末利用率有多高?

这些问题在实验里要切样、腐蚀、看金相,还要用高速相机拍粒子轨迹。而在 Flow-3D 里,我们只需要在熔池凝固后提取沉积层的几何轮廓,再统计捕获粒子的数量,就能得到一组可比对的数据。把这组数据在不同工艺参数下扫一遍,就建立了"仿真工艺窗口",后续上机床试制时可以直接跳过一大批盲试。

2. Flow-3D 建模仿真的关键设置:从 CAD 模型到粉嘴模型

很多新手卡在第一步的"设置"上,不是因为软件难学,而是因为没有想清楚哪些物理过程要写进去、哪些可以砍掉。我自己的原则是:第一版模型尽量简化,只保留三类要素——基板与熔池、激光热源、粉末粒子注入。喷嘴气流如果暂时不研究粉斑汇聚半径,可以先不建,直接给定粒子入射速度、方向和分布范围。

2.1 几何模型与网格划分

先说几何。基板通常是一个长方体,尺寸不需要太大,能覆盖熔池影响区即可。比如单道熔敷,基板取 20 mm × 10 mm × 5 mm 足够,太大的话计算域会浪费大量网格。喷嘴如果要做同轴气流研究,可以把喷嘴简化为同轴圆环出口,出口处给环形速度边界。

网格划分是 Flow-3D 仿真的重中之重。激光光斑直径一般 1~3 mm,粉末粒径 50~150 μm,熔池深度 0.5~2 mm。这意味着整个计算域里的网格尺度必须同时分辨光斑能量分布和粒子与熔池的交互。

我自己常用的配置是:在激光光斑和熔池区域设置嵌套加密块,最小网格边长为粉末粒径的 1.5~2 倍。比如粒径 70 μm,网格边长取 100~150 μm。如果用均匀网格,整个计算域网格数会爆炸,但只对局部加密区做细网格,就能兼顾精度和计算时间。

网格尺寸和粒子尺寸的关系还有一个容易忽视的点:如果粒子直径只有网格尺寸的十分之一,粒子与流体的相互作用会被严重低估。在做"粉末撞击熔池液面后是否被捕获"这类问题时,细网格和粒子捕捉条件必须同时满足,否则粒子会穿透熔池液面直接飞出计算域。

2.2 物理模型选择

打开 Physics 选项卡,要勾选的模型包括:

  • VOF 模型(默认开启,用于追踪金属熔体与空气的自由表面)
  • 凝固/熔化模型(Solidification Model),设置材料的固相线温度、液相线温度、熔化潜热
  • 传热模型,包括对流和辐射。高温熔池表面辐射散热不可忽略,通常设置表面辐射率 0.3~0.5
  • 粘性流模型,层流即可。熔池尺寸小,流速虽然受马兰戈尼力驱动,但一般还是层流占主导
  • 粒子模型,选择 Particle Injection / Entrained Particles,按粉末实际物理属性设置密度、直径、初始速度和数量的分布

在流体属性 (Fluid Properties) 里,金属熔体需要给出随温度变化的密度、粘度、导热系数和比热。这里要说一个坑:很多人图省事,设置恒定的粘度,结果熔池里的马兰戈尼对流完全失真。粘度随温度的变化直接决定了表面张力梯度驱动的熔池内环流形态,这个环流又决定了熔池深度和稀释率。我是把纯铁或钛合金的随温度物性数据整理成表格,再录入软件。

2.3 激光热源怎么给最接近真实

Flow-3D 的激光热源一般用高斯分布来近似。对于连续激光熔敷,可以采用表面热源或者体热源模型。如果是深熔焊那种匙孔模式,体热源更合适;但送粉熔敷通常热输入密度相对低,属于热传导模式,表面热源就够用。

高斯面热源的表达式大家可能很熟:在光斑半径 r0 内,热流密度按 exp(-2r²/r0²) 分布。关键是要乘以材料对激光的吸收率。很多人拿真实激光功率直接输入软件,算出来的熔池深度比实验深一大截,就是因为没有考虑吸收率。对于波长 1 μm 左右的激光,钢表面吸收率大概 0.3~0.45,钛合金更低一些。如果表面有金属粉末覆盖,吸收率又会变化。这个系数最好用"仿真-实验对照"反推一次,之后再做趋势研究就准多了。

激光热源在软件里通常以专门的 Boundary Condition 或 Volume Heat Source 形式施加。注意扫描移动方向的拖尾效应——熔池尾部温度梯度比前部缓,热源不能简单静止。需要把热源位置写成随时间移动的函数,或者用网格块移动 + 热源跟随的方式实现。

2.4 粉末粒子的注入与运动

粉末粒子设置是本文另一个重点。Flow-3D 中注入粒子时,有几个参数直接影响仿真可信度:

  • 质量流率。由送粉量决定,比如 8 g/min,要换算成每秒注入的粒子总质量。
  • 粒径分布。实际雾化粉末不是单一粒径,而是服从 Rosin-Rammler 分布。可以把粒子注入设为多组直径,比如 50 μm、75 μm、100 μm,按质量占比分配。
  • 初始速度和方向。从喷嘴出口射入时,速度由载气流量和喷嘴结构决定。如果没有建喷嘴模型,就需要估算粉末出口速度。常见 Ar 载气流速 5~15 m/s。
  • 粒子与熔池的相互作用。激光熔敷里,粒子碰到熔池后应当被捕集、熔化并入流体域。Flow-3D 粒子模型中有多种边界条件,可以把熔池自由表面设置为"粒子捕集"(Stick/Capture),或者设置温度阈值,粒子到达熔池表面且温度高于液相线时并入流体。

我实际跑下来发现,粒子捕集阈值设置对结果影响很大。如果条件太宽松,粉末会被"吸"进熔池,利用率虚高;如果条件太严,又会出现粒子落进熔池表面又被反弹出来的怪现象。一般来说,粒子与熔池接触时只要速度法向分量低于某个阈值、且熔池表面位置在液相线以上,就可以判为捕获。

3. 送粉量与激光参数耦合:仿真中的粒子追踪与熔池定量化分析

模型跑通之后,重点就转移到了"怎么从仿真里提取有用的工艺结论"。我自己的习惯是,把仿真结果拆成三条线来读:粉末利用率、熔池几何尺寸(宽度、深度、高度)、温度场时间演化。

3.1 从粒子追踪结果里提取粉末利用率

Flow-3D 后处理里可以按粒子 ID 追踪每个颗粒的完整运动轨迹。把进入计算域的总粒子质量,和最终留在沉积层中的粒子质量做比值,就是粉末利用率。

这个数字在送粉熔敷里非常敏感。我举个例子:固定激光功率 2000 W、扫描速度 10 mm/s,送粉量从 5 g/min 提高到 15 g/min,粉末利用率会先升后降。低送粉量时所有粉末都有足够能量熔化,但沉积层薄;高送粉量时,大量粉末来不及吸收足够热量,部分粒子未熔化就被后续粉末"埋"进熔池,形成夹杂。

粒子轨迹图还有一个用法是检查喷嘴角度。我调过一个 45° 倾角喷嘴的案例,仿真显示下方粉末流受基板阻挡弹开,粉末汇聚点向扫描方向后方偏移,导致熔池前后部接受粉末不均匀。后来把喷嘴角度改成 55°,粉末汇聚点与激光光斑中心重新对齐,沉积层表面鱼鳞纹才均匀。

3.2 熔池几何尺寸与稀释率的定量提取

熔敷质量的两大指标是稀释率和接触角。稀释率定义是熔池内基板熔化深度 /(熔化深度 + 沉积层高度)× 100%。稀释率太低,涂层与基板结合不好;太高,合金元素被过度稀释,熔敷层性能下降。

在 Flow-3D 后处理里,沿着熔敷层中心取一个横截面,冻结熔池边界,直接就能读到熔深和沉积高度。我通常会用软件里的"探测曲线"功能,或者干脆把温度场等值线(固相线温度)导出到文本文件,用 Python 或者 MATLAB 算面积。这个比在界面里目测要准得多,也方便批量处理多组参数。

3.3 参数耦合的"工艺窗口"仿真方法论

单道熔敷的工艺参数主要有四个:激光功率、扫描速度、送粉量、光斑直径。它们不是独立起作用的。比如激光功率和扫描速度之比就是线能量密度,决定了熔池温度场;送粉量和线能量密度的比值决定了"单位质量粉末获得的能量",这个值是判断未熔缺陷的最佳指标。

我一般会做一组正交仿真矩阵,参数范围和取值如下表:

参数取值说明
激光功率 P1200/1600/2000/2400 W对应不同热输入
扫描速度 v6/10/14 mm/s影响熔池驻留时间
送粉量 m5/10/15 g/min影响沉积效率与缺陷
光斑直径 d2/3 mm影响能量密度分布

每组仿真算完之后,把稀释率、粉末利用率、沉积层高度整理成表格,找出"稀释率在 10%~25% 区间、无未熔颗粒"的参数组合,就得到了推荐的工艺窗口。

3.4 一个可复现的典型算例参数

为了让没有头绪的朋友能直接开跑,我把自己验证过的一个基准案例参数列出来。基板为 316L 不锈钢,粉末也用 316L,粒径 75 μm,激光功率 1500 W,光斑直径 2 mm,扫描速度 8 mm/s,送粉量 8 g/min,载气流速 6 m/s,喷嘴距基板高度 12 mm。

计算域 20 mm × 8 mm × 5 mm,最小网格 0.15 mm,网格总数约 350 万。粒子注入 10000 个/秒,总模拟时间 1.2 s。这个配置在普通 8 核工作站上大约要算 10~14 小时。结果里熔敷层高度约 0.8 mm,熔宽约 2.6 mm,稀释率 16%,粉末利用率 42%,和实验结果吻合度比较高。

4. 跑通案例之后的那些坑:网格、粒子离散相与温度场的务实排查技巧

仿真跑起来只是第一步,结果能不能用是另一回事。我前前后后在这套模型上踩过的坑可以列一份长长的清单,这里挑几个最典型的说透。

4.1 熔池凝固边界"跳变"与网格敏感性

以前用 Flow-3D 做激光熔敷,经常遇到一个诡异问题:熔池宽度每隔几个时间步就跳一下,曲线锯齿状。后来排查发现是网格太粗,固液相变界面在格子之间反复横跳。熔池固相线附近的温度梯度极陡,如果网格尺寸和熔池深度在同一量级,界面的数值分辨率就不够。

解决思路有几个层级:第一优选是加密熔池附近网格,让网格尺寸小于预期熔深的一半;第二是把凝固模型中糊状区参数调宽一些,让相变释放潜热的过程更平滑;第三是加强时间步限制,在熔池区域用更小的对流 CFL 数。我后来发现 Flow-3D 凝固模型里有个"mushy zone"参数设置,如果把它设太窄,数值求解稳定性会变差。适当放宽糊状区范围,不仅计算更稳,熔池轮廓也更接近实验。

4.2 粒子追踪"穿透"熔池表面:单向耦合的局限

接前面提到的粒子捕集问题,再展开细说。很多 CFD 软件默认粒子与流体是单向耦合,也就是说流体影响粒子轨迹,但粒子本身的质量和动量不反馈给流体。对熔敷粉末这种体积分数很低的场景,单向耦合大部分情况可行,但在"大量粒子持续冲击熔池液面"时,粒子群对液面的挤压效应不可忽略。

我遇到过一次非常离谱的结果:粉末粒子全部从熔池顶部穿过去,就像熔池根本不存在一样。最后发现原因有两个:一是粒子时间步长太大,粒子在单个时间步内移动距离超过几个网格,直接跳过了液面捕集层;二是粒子与流体之间的拖曳力计算在低网格分辨率下失真,粒子穿透了自由表面。

解决办法是把粒子子时间步长强制缩小,通常是 1/10 或 1/20 个流体时间步。同时把自由表面处的网格层数加密到 3 层以上,确保粒子穿越自由表面时一定有充足的空间触发捕集条件。

4.3 温度场初值对熔敷层形貌的影响

这个相对隐蔽。新用户在 Flow-3D 里通常设初始温度 25 ℃,但实际激光熔敷是连续过程,基板温度会逐渐升高。如果只仿真 1 秒钟内的单道熔敷,初始温度的影响不大;但若仿真 5 秒以上的连续扫描,基板尾部不断被预热,熔池体积会越来越大,沉积层也随之变宽。

我在做多道搭接仿真时,会把上一道熔敷结束后的温度场作为下一道的初始条件。Flow-3D 支持从一个 restart 文件延续计算,可以把整个扫描过程分成多段,每段的起始温度场自动继承上段结果。这样算出来的搭接区温度分布才真实,不然第二道熔敷的稀释率完全失真。

4.4 实验验证:仿真必须和一次实验"握手"

不管仿真做得多细,都不能完全替代实验。我的建议是,哪怕只做一次单道熔敷实验,也要把仿真模型的关键参数校核一遍。需要对比的物理量包括:熔高、熔宽、熔深、稀释率、粉末利用率。如果这几个量能和实验对得上,仿真模型才具备"工艺预测"的价值。

校准顺序也有讲究。先校激光吸收率,因为它的误差对熔池尺寸影响最大。调整办法是:固定所有参数,用不同吸收率跑三组仿真,把熔池深度与实验切片的熔深对比,线性插值。之后再校粉末捕集阈值,因为它直接影响沉积层高度和粉末利用率。最后才校边界对流换热系数。

5. 一些个人经验和后续可以继续做的事

仿真这个领域,入门容易,做深难。激光送粉熔敷尤其考验对物理过程的直觉——你得能想象出熔池表面那个凹陷在移动时怎么拉扯周围的液体,想象粉末粒子拖着高温尾迹冲进熔池时会不会被蒸汽喷出来。

我个人的习惯是,每跑一组工艺参数,都顺手把熔池内的流线、粒子轨迹、温度等值面各导出一张图存档。等参数多了之后回看这些图,往往能看出单一时刻看不出来的规律。比如有一次我发现,当激光功率超过某个阈值后,熔池内部出现了一个非常强劲的回流涡,这个涡把已经凝固的粉末重新卷进熔池,造成局部成分偏析。这类现象在实验金相里很难解释,但在仿真流线图里一目了然。

另外一个非常实用的小技巧是:用无质量粒子做熔池示踪。在熔池内部某个局部注入一批无质量粒子,让它们跟随熔体流动,后处理时按粒子轨迹着色,就能看到马兰戈尼对流的完整路径。这不仅帮助你理解熔池搅拌效果,还能直接验证稀释率是不是主要由对流贡献的。

后续要继续扩展的话,我建议往三个方向走:一是把粒子群算法和仿真结合,比如用粒子群优化(PSO)算法自动寻优激光功率、扫描速度、送粉量,避免手工试凑;二是做多道多层搭接仿真,研究热循环对残余应力的影响;三是加入保护气与金属蒸汽的组分扩散,研究激光熔敷过程中的烟尘和气孔形成机制。

老实说,我最初用 Flow-3D 做激光送粉熔敷仿真时,也是从"不知道怎么设粒子"开始的。后来把粉末利用率、熔池形貌、工艺窗口这一整套流程跑通,再回头看,仿真不是要取代实验,而是要把实验前的盲试时间压缩到最小,让你上机床那一刻更有底气。如果你正准备开这个方向,我的建议是:先别贪多,从单道、单参数基准仿真开始,跑通一个可信案例,再慢慢往上加复杂度。这条路,值得走。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询