★ FEATURED ARTICLE 做激光熔覆仿真的朋友应该都有过这种经历单道熔覆的模型调通了温度场看着挺像回事但只要一把手伸到“多层多道”各种问题就全冒出来了——层与层之间热量来不及散、第一道还没冷却第二道又压上来熔池形状乱跳算着算着直接不收敛就算算完了应力结果也总让人觉得不可信。这篇文章主要聊我基于COMSOL做多层多道激光熔覆模型的完整思路包括物理场怎么选、生死单元怎么处理、移动热源怎么加载、热应力怎么耦合以及做“仿真视频”时怎么把结果展示得既准确又直观。适合正在做激光熔覆、增材制造类仿真的研究生和工程师参考尤其是被多层多道“多道搭接”和“逐层堆积”搞到头皮发麻的那批人应该能少走不少弯路。1. 多层多道到底难在哪从单道焊斑到堆积体的物理跃迁先聊一个很多人没太当回事、但后面处处踩坑的问题多层多道不是“单道复制粘贴”它在物理上已经变了一个东西。单独算一道激光熔覆激光扫过去熔池形成冷却基体只是一个有限厚度的受热板边界条件相对简单。到了多层多道情况完全不一样了。1.1 热量累积效应每一道都站在前一道的“余温”上第一道熔覆完成后基体表面温度并不会立刻降到室温。实际工况里下一道往往在十几秒甚至几秒后就开始了。多道搭接时激光扫描路径是蛇形来回的前一道刚扫完的位置可能还保持一两百摄氏度甚至更高后一道马上又来加热。这样造成的直接后果是熔池的最高温度逐道升高熔池尺寸逐道扩大冷却速度却逐道变慢。这种热量累积在仿真里如果不处理就会出现一个很尴尬的画面第一道温度场看着舒服第二道还行第三道开始熔池温度逼近甚至超过材料沸点然后求解器开始叫唤或者结果里出现一片花花绿绿的“高温噪音”。所以要诀是仿真模型里的每一道都必须真实继承上一道结束时的温度场作为初始条件而不是每一道都从“室温开始”。这是多层多道与单道最大的物理差异也是最容易让新手模型“失真”的地方。1.2 熔池流动与固液界面忽略什么、保留什么COMSOL做激光熔覆理论上最完整的是共轭传热加层流把熔池里的Marangoni对流做出来加上固液相变、辐射、蒸发潜热一大堆。一个真实的单道熔池模型网格动辄上百万单核跑一个晚上不稀奇。到了多层多道如果每一道都开完整流动算到半年后见吧。个人经验是多层多道场景下先放弃熔池流场细节把核心放在“逐层热量输入—散热—固液相变—残余应力”这条主线上。因为多层多道最关心的往往是搭接率、熔覆层形貌、热影响区范围、残余应力分布和变形趋势而不是某个瞬间熔池里涡流长什么样。用“有效导热系数”去近似熔池对流带来的热量均化是完全够用的工程简化。1.3 生死单元的引入堆积过程的离散化表达多层多道的几何不是一开始就全部存在的。激光没到的地方粉末还没铺上去那一层熔覆层物理上还不存在。如果一开始全建模你会发现基体还没被加热上面的“空中楼阁”却已经在导热了热量直接穿过还没形成的位置跑到前面去结果完全是错的。生死单元是解决这个问题的标准手段。简单理解就是几何模型先建全但在求解开始之前把尚未熔覆的单元全部“杀死”——不让它们参与导热计算、不承受热载荷直到激光扫描到对应位置时再“激活”它们使其参与后续计算。这有点像工地先搭好所有楼层的脚手架但哪层封顶哪层才开始住人。每一层熔覆层单元被激活的时间必须和激光扫描到该位置的时间严格对应起来。2. COMSOL几何装配与网格控制层次激活和生死单元的落地思路很多人卡在“生死单元”这一步倒不是不理解原理而是不知道在COMSOL里具体怎么落。这里分享一套我验证过多轮的实操思路适合做多层多道的层状搭接模型。2.1 几何分层建模把每道熔覆层切成一格格“积木”先在CAD里或COMSOL的几何节点里把每一道熔覆层单独画成长方体或弧形截面的独立域。比如做三层五道那就是三个层、五个条状域共15个独立域。这些域在几何关系上是“相邻”但不是“联合体”默认的连接关系最好使用Form Assembly装配体让每个域的网格可以独立生成和控制。为什么要用装配体而不是联合体因为后续做“逐道激活”时你会需要把某一道的域单独选中、单独赋材料、单独控制网格。用联合体的话域之间会被强制合并界面后期选择变得非常别扭。装配体虽然会让界面处多一层内部边界但在热分析里影响很小完全可以接受。这一步做完你会得到一个“齿梳状”的几何。如果觉得麻烦也可以用Work Plane在基体上方拉伸出面、再分割成条状域思路一样。2.2 未激活域的“准真空”处理弱导热而不是直接删除在COMSOL的标准传热接口里没有现成的“生死单元”勾选框。最工程化的替代手段是为每道熔覆层设置一个“激活状态”变量用解析阶跃函数Step函数控制它的导热系数和比热容。具体做法是未激活前把该域的导热系数设为一个极小的值比如实际值的1e-6密度设为实际值但比热容设到一个不产生热量缓冲的极小值等效于“这个位置不存在”。激光扫到该位置时Step函数在几个时间步内从0切换到1这一道的材料参数恢复到真实值单元开始正常吸热、导热、散热。有同学问那我不激活的时候把导热系数设为0不就行了千万别。导热系数为0会让传热方程在该区域退化矩阵条件数急剧恶化不收敛概率大增。设成1e-6的弱导热既保证“热量不穿过去”又维持了方程的病态程度可控实测非常稳。2.3 网格无关性与扫描路径附近的加密策略多层多道仿真的网格量一定要省着用。激光熔覆的局部热源非常集中熔池半径也就一两毫米可基体可能是一块几十毫米甚至上百毫米的板。如果全用细网格自由度轻松破百万工作站都扛不住。我通常的做法是“三级网格”激光扫描路径正下方的熔覆层和基体表层网格尺寸压到熔池半径的四分之一到五分之一约0.2~0.5mm基体中下层、远离扫描线的两侧网格尺寸放到5~10mm中间过渡区用自由四面体或者映射网格做渐进过渡避免尺寸跳变过大。网格无关性验证也不能省。我一般是固定热源参数不变把熔覆层附近网格从1.0mm开始往下加密对比同一时刻熔池最高温度和熔池宽度当温度差异小于2%、熔池宽度变化小于3%时认为网格基本无关。多层多道模型不必每次都做全套无关性验证但第一层第一道的模型必须要做一次后续在相同热源下沿用这个网格密度即可。3. 移动热源与熔池热行为的建模细节生死单元解决了“材料什么时候出现”的问题但热量怎么加、加在哪里、怎么移动是另一个大坑。移动热源的建模质量直接决定温度场的可信度。3.1 热源模型选型高斯面热源与双椭球体热源的取舍激光熔覆最常用的热源有两类高斯面热源Gaussian surface heat flux和双椭球体热源Double ellipsoid volumetric heat source。单道快速粗算可以用高斯面热源表达式简单参数也少。但多层多道场景里我更推荐双椭球体热源因为激光能量不是只作用在表面粉末层吸收能量后会把热量带进熔覆层内部体热源能更真实地反映这个“体积加热”过程。而且双椭球前后半轴长度可以分别设置能反映激光扫描时熔池前后不对称的实际形状——前半个椭球短前方材料未预热后半个椭球长后方有拖尾热影响。双椭球热源在COMSOL里实现时常见做法是用一个“局部坐标系”跟随激光位置移动在热源表达式中用空间坐标减去激光当前位置把热源中心始终固定在激光所在坐标处。表达式大致是Q(x,y,z,t) (6√3·η·P)/(a·b·c·π√π) · exp(-3(x-x0(t))^2/a^2) · exp(-3(y-y0(t))^2/b^2) · exp(-3(z-z0(t))^2/c^2)其中x0(t)、y0(t)是激光在扫描平面内的瞬时坐标用参数化扫描或瞬态求解的时间表达式控制。照搬公式之前有一个更基础的问题要先想清楚熔覆层在空间中是不断“新增”的。如果单纯把热源加载在整个域的方程里未激活区域会因为材料被“弱化”而不产生响应但热源项本身还是会施加在那些区域上。所以一个更稳妥的操作是在热源表达式里乘上该位置的“激活状态值”让热源只作用于已经激活的熔覆层域不会白费力气去加热一片“真空”。3.2 热物性参数与相变潜热的等效处理多层多道仿真对材料的温度相关物性非常敏感。钢材从室温到熔点的热导率和比热容变化很大如果全套用室温值熔池温度会被高估百分之二三十都不稀奇。我建议至少在材料节点里定义热导率随温度变化至少定义20°C、200°C、400°C、600°C、800°C、熔点附近、超过熔点几个点比热容随温度变化同样多点定义熔融相变潜热用“表观比热容法”处理即在一个小的温度区间内把潜热折算为比热容的峰值叠加进去。表观比热容法实现起来非常方便在材料定义里把一个宽度为ΔT的相变区间内的比热容设为基础比热 L/ΔTL是熔化潜热ΔT取20~50K。这样既不会引入额外的“相变接口”增加计算负担又能在温度曲线里看到明显的“平台期”非常实用。实际仿真中很容易忽略的是辐射和表面对流换热。激光熔覆的冷却过程中熔池和高温区域的对流散热与辐射散热占比非常大。在传热物理场里封闭边界默认是绝热的必须手动给所有暴露在空气里的边界加上“对流热通量”和“表面对环境辐射”两个条件。这里注意不要写成“表面对表面辐射”那个会引入视角因子计算成本高且多道场景下没必要。3.3 扫描路径与层间冷却时间时间步长的设定逻辑多层多道的扫描路径通常是蛇形往返每一道与上一道间隔几毫米层与层之间可能还有冷却等待时间。这些在模型里都要按真实工艺写清楚。时间步长的设定有个朴素的规则每一时间步激光移动的距离不要超过最小网格尺寸。比如激光扫描速度是5mm/s网格是0.25mm那时间步长建议不超过0.05s。再结合数值稳定性条件我一般先按0.1s步长粗跑一遍看温度场有没有锯齿状跳跃如果有就把步长减半。算完一道后风速冷却或者空冷等待阶段可以放大时间步长到0.5~1s这样能省不少计算时间但不会牺牲熔池细节。这里有个小技巧COMSOL的瞬态求解器可以使用“自适应时间步长”但峰值的捕捉和步长的稳定性很多时候不如固定小步长来得可靠。多层多道模型我建议先用固定步长跑通一道再切换自适应步长提升计算效率两步配合。4. 从热分析到结构分析热应力、变形与裂纹风险评估温度场只是第一步多层多道最终关心的往往是熔覆完成后能不能用——也就是残余应力、变形和裂纹风险。COMSOL把温度场结果作为载荷再做结构力学分析这个过程有一个非常容易被忽略的坑必须单独讲。4.1 单向耦合还是双向耦合严格来说熔覆过程存在热力双向耦合应力变形反过来会影响接触和热通量和变形热但这在激光熔覆里非常小。工程计算几乎都采用单向耦合——先求瞬态温度场再把结果导入结构分析求出应力和变形。这是业界标准的简化放心用。在COMSOL里一种常见做法是同时建立“固体传热”和“固体力学”两个物理场在固体力学里把温度变量T直接作为体载荷引用一起瞬态求解。这样省去了数据传递但计算量会翻倍。另一种做法是分两步先算传热保存瞬态结果再建立第二个研究类型比如“瞬态结构”用“固体传热”的结果作为温度载荷。多层多道模型由于自由度已经很大我倾向于第二种至少内存压力小一些。4.2 参考温度与材料高温性能的坑结构分析引用温度载荷后材料必须定义热膨胀系数、弹性模量、屈服强度随温度的变化这个千万不能省。尤其是屈服强度很多材料在800°C以上只剩室温屈服强度的百分之二三十。如果从头到尾用一个固定屈服强度冷却时的塑性区会被严重低估残余应力结果会显得“太乐观”。还有一个非常容易忽略的点热膨胀的参考温度。在COMSOL的固体力学设置里热膨胀默认的参考应变温度是“零应变参考温度”。如果忘了设置材料初始时刻就会因为从0°C膨胀到室温而产生热应力结果里会凭空多出一大块应力怎么看怎么离谱。正确做法是结构分析里的“初始条件”要和传热分析的初始温度一致。比如传热初始场是室温20°C结构力学分析的初始热应变参考温度就应该设成20°C这样初始时刻没有任何热应变后面所有的应力都来自温度变化。4.3 逐层应力历史的累积问题多层多道的残余应力不是“最后一道冷却完”的那一个瞬间定的。每一道熔覆、冷却都会在前一道基础上叠加新的应力。如果仿真的结构分析是从头到尾瞬态跟踪的那天然会带上历史累积但代价是计算量巨大。工程上有个折中办法把每一道结束后的应力场作为下一道的初始应力条件逐道“接力”计算。这样重心只放在每道结束时刻与下一道开始时刻的衔接上中间几十秒的空冷期不细化计算。我实测下来和前一种全程跟踪的结果差异在5%以内效率却能提升好几倍。判断开裂风险时我习惯同时看两个量冷却到室温后的最大主应力判断脆性开裂风险和Von Mises应力判断塑性流动和整体失稳风险。另外熔覆层与基体交界处的应力集中明显可以额外切一个截面图观察该处的应力梯度变化。如果看到交界处应力超过材料高温屈服强度那实际工艺大概率会在这个位置起裂。5. 多层多道仿真的高频坑与排查路线多层多道仿真做得久了会遇到很多奇奇怪怪的错误和反直觉的现象。这里把我踩过的、以及帮别人排查过的最常见的几类问题整理出来直接说结论和排查路径。5.1 求解发散先看热源功率加载曲线再看步长最常见的发散原因有两个一是初始时刻热源突然全功率加载相当于冷的材料被电焊枪怼上去时间步长稍大数值瞬间爆炸二是未激活区域的弱导热系数设置得太小导致矩阵几乎奇异。排查路线我是这样走的先把热源功率在0.1s内从0线性上升到目标值这一招能解决一半以上的发散然后把时间步长减小到原来四分之一看是否稳定最后检查弱导热系数的量级最小不建议低于实际值的1e-5否则宁可多承担一点点寄生热量也要保证矩阵良态。我见过最离谱的一次发散花了半天时间才发现是热源坐标系的方向反了激光在“倒着走”。所以排查发散时第一步先看热源运动轨迹的动画确认扫描路径和工序一致别急着调求解器参数。5.2 熔池温度高得离谱吸收率、热源体积和网格尺寸三连查如果熔池最高温度超过材料沸点一两千度通常不是物理问题而是模型参数问题。激光吸收率是最先要查的不锈钢粉末对光纤激光的吸收率一般在0.3~0.5之间不是0.9。很多教程为了温度好看会把吸收率调得很高实际会失真。其次是热源椭球的体积热源能量密度是功率除以椭球体积如果椭球尺寸设小了能量密度自然爆表。最后看网格网格太粗的情况下热源作用区域被平均化峰值温度会被低估但这个和“温度爆表”往往是反方向——网格太粗更常见的是熔池温度不足。不过也有一种特殊情况网格太粗时如果时间步长也大数值振荡会在个别节点上打出“尖刺温度”这个用5.1的排查思路就能定位。5.3 冷却阶段温度锯齿空间步长和时间步长不匹配冷却阶段温度曲线出现锯齿状波动几乎都是网格局部尺寸和时间步长不匹配造成的。特别是在熔覆层和基体交界处如果网格尺寸跳跃太大冷却收缩时会造成局部数值振荡。处理方法很简单把过渡区网格的尺寸变更加平滑或者把冷却阶段的时间步长缩小到熔覆阶段的一半。如果锯齿仍然存在可以在后处理时对结果做“平滑平均”但那只掩盖问题最好还是回到网格适配。5.4 应力结果“悬空”的应力云图应力结果出来后经常会看到熔覆层表面一片应力云图但仔细一看基体底部完全没有应力。这不一定是错的——如果基体足够厚底部约束较弱应力主要集中在表面热影响区完全合理。但如果应力云图里有明显的“飘浮感”比如某一道的应力场和相邻道不连续多数是因为逐道接力计算时初始应力没有继承好。检查一下结构分析的初始条件是不是包含了前一道的应力场而不是默认从零应力开始。6. 视频呈现与模型配套让仿真结果可信又看得懂标题里提到“视频与模型”这其实是个挺值得展开的点。COMSOL仿真结果如果只是给导师看或者写论文静态图就够了。但如果要演示工艺过程、输出给现场工程师看一段流畅的仿真视频比十张云图都管用。视频做得好不好牵涉到一个问题你能不能让观看者一眼看懂每道工序在干什么。6.1 视频的核心帧设计熔池、热影响区和应力演变三视图我做仿真视频时一般会固定三个视角俯视图显示激光热源移动、熔池拖尾、搭接带的过渡剖切侧视图显示熔覆层与基体交界处的热影响区深度变化一个旋转视角显示整个基体的应力和变形演变。关键是这三个视角务必使用同一套色标温度色标、应力色标。很多演示视频的问题在于每个视角换了色标观众根本没法对比。COMSOL导出动画时可以把三个视图放到同一个布局里同步导出导出分辨率设置为1920x1080帧率控制在每秒10~15帧就够了太高会让文件大得没法打开。6.2 时间轴标注把工序时间融进画面多层多道的视频如果没有任何时间轴标注观众会看到一道一道光在板子上跑来跑去但完全不知道对应实际工艺中的第几道、什么时刻。我通常会在模型三维视图里添加一个文本显示实时更新当前时间、当前道次和当前激光坐标这样画面信息量立刻就不一样了。COMSOL里可以通过“全局变量探针”把当前时间、温度最大值等变量放上视图再用“显示坐标”控制文本位置。这个小细节花不了十分钟但对视频质感和可信度提升非常明显。6.3 模型文件的组织让别人拿到就能跑如果你的模型要分享给同行或学生复现文件组织就很重要。COMSOL的.mph文件是有版本兼容性的分享前一定要注明软件版本号和组件版本号比如“COMSOL 6.2 传热模块结构力学模块”。另外建议把全局定义参数、材料参数、事件时间表都放在“全局参数”里集中管理导出一份Excel表格附在模型文件旁边。最实用的一点在后处理里把“第一道结束”“第二道结束”“冷却到室温”几个关键时间节点的截图事先导出成PNG放到一个说明文档里。这样别人打开模型后即使不重新求解也能通过这套截图快速判断计算结果是否正常。6.4 一个“以帧代全”的验证思路最后分享一个验证极有帮助的习惯在模型里放置几个固定探针坐标对应真实工艺里可能布置热电偶的位置例如熔覆层中心点沿深度方向每1mm取一个点。仿真算完后把所有探针的温度-时间曲线导出来和单道实验的测温曲线做对比。曲线形态和峰值温度的误差在10%以内这个模型的可靠性就基本站得住。如果误差明显优先回头检查热源吸收率和边界散热系数而不是去调网格。我曾经用这个方法帮某项目排查出边界对流系数被高估了一倍的问题修正后温度曲线和实测几乎重合。这个习惯建议一开始就养起来而不是等到模型算完后再补。7. 最后再分享一个小技巧善用“逐道重置”而不是“一道重置”多层多道仿真里还有一个很多人不太注意的问题每算完一道要不要把变量重置成初始状态。如果你把整个求解过程设置为一个连续瞬态让激光走完第一道、第二道……直到最后一道那所有道都会在同一个时间轴上求解一道叠一道热量累积自然出现。但如果你偷懒把“每道”设置成独立的研究步骤研究之间又不共享“存储解”那第二道默认从室温重新开始得到的结果根本就不是多层多道而是把单道重复算了几遍。正确的做法是在一个“瞬态研究”里定义事件时间序列把扫描、冷却、换道这些动作全部按时间排进同一个时间表。激光功率的启停、生死单元的激活、每道扫描坐标的切换全部由时间变量来控制。COMSOL的“事件接口”可以很方便地在特定时间点触发状态变化比如一道结束立即关闭热源、开启空冷下一道开始时再打开热源。把这个时间流程理清了多层多道仿真的“魂”就抓住了后面一切结果都顺畅了。激光熔覆的多层多道仿真说难也难说通也通。难在物理机制多、工序环节杂、计算资源有限通在其实核心就这么几件事——热量累积的正确继承、生死单元的准确切换、移动热源的稳定加载、应力历史的逐层累积。把这四条主线吃透再配合一套可靠的网格策略和排查路线模型跑出来的结果基本就能和真实工艺对上了。希望这篇稿子能帮正在跟熔覆模型较劲的朋友们省下几个熬夜的晚上。
阅读完成 · 觉得有帮助？
