从焊接工艺和模拟仿真两头打交道久了,我越来越觉得,电弧和熔池根本没法分开谈。无论是TIG电弧熔池一体化,还是MIG电弧熔滴蒸汽一体化,本质都是同一件事:把电弧、熔池、熔滴、金属蒸气当成一个相互咬合的耦合系统来理解和控制。现场遇到的咬边、飞溅大、熔深不足、气孔密集,追根溯源,问题全都出在电弧和熔池/熔滴之间那点"纠缠"上。
早期做模拟的人习惯把电弧当成一个固定热源直接砸在工件上,算熔池的流场和温度场;或者反过来,把熔池表面当成一个给定温度的边界,去解电弧的温度分布。分开算在当时是计算资源限制下的务实选择,可一旦推到深水区——脉冲焊参数优化、低飞溅控制、短路过渡临界判断、增材制造中的熔池闭环调节——就发现完全不够用了。原因不复杂:电弧的热流密度和压力分布是由熔池表面的温度决定的,而熔池的流动形态又反过来被电弧的热、力、电作用塑造,这是典型的强耦合问题。老办法用了这么多年,模拟结果和实测总是差一口"气",差的就是这口耦合的气。
这篇文章我把这两块掰开揉碎讲清楚:TIG焊里电弧和熔池怎么在一个体系里算,MIG焊里熔滴过渡和金属蒸气怎么改变电弧行为,以及在实际建模中哪些步骤最容易被卡住。适合正在做焊接数值模拟、焊接工艺窗口开发或者想搞懂电弧物理本质的人参考。
1. 先从现场聊起:为什么要把电弧和熔池捆在一根绳上
焊接电弧本质上是一个被电离的气体导体,电流从钨极(TIG)或焊丝(MIG)流过电弧,再流进工件。这个"气体导体"的温度在弧心可以到一万多开尔文,而熔池金属表面温度只有两千开尔文左右。两者之间隔着一个厚度只有几毫米、温度梯度却高达每毫米几千开尔文的区域。这个区域里的热交换、动量交换、电流传导,决定了熔池的形状、熔深、飞溅倾向,也决定了电弧本身的稳定性。
老的分开建模思路有一个隐含假设:电弧给熔池的热流密度分布是预先给定的"标准曲线",熔池表面的变形和温度变化不会反过来改变这条曲线。实事完全不是这样。熔池表面一旦凹陷,电弧弧长随之改变,阴极区和阳极区的边界条件跟着变,热流密度分布就被重新洗牌了。TIG焊中常见的"电弧挺度"现象——熔池变深时电弧明显变亮变高——就是这种耦合关系的直接体现。
1.1 一体化模型到底"一体"在哪里
所谓TIG电弧熔池一体化,简单说就是把钨极尖端、电弧等离子体区、熔池液态金属区、工件固体区放在同一个计算域里求解,不再单独给电弧和熔池分别设定"人为耦合边界"。这样做的第一个好处是,电弧在熔池表面的热输入、电流密度分布、剪切力分布全部由计算自动得出,而不是靠经验公式拍脑袋。第二个好处是,熔池表面的自由变形可以在同一个计算框架内实现,不需要人为假设"熔池是平的"。
对MIG电弧熔滴蒸汽一体化来说,"一体"的含义更进一层。MIG焊多了一个焊丝熔化后的熔滴过渡过程,熔滴在焊丝末端形成、长大、脱离、飞行、进入熔池,整个过程与电弧等离子体、金属蒸气之间实时相互作用。金属蒸气的浓度分布直接改变电弧的电导率和温度场,电弧温度场又改变熔滴受到的电磁力和等离子体流力。要把这一整套东西一起解出来,就需要把气相、液相、固相、金属蒸气浓度当成一个统一的多物理场体系。
这种"一体化"在物理上更接近真实,在数学上则意味着控制方程更加复杂。我曾对比过老式解耦方法和一体化模型对同一组TIG参数的预测结果,最大差异出现在熔池自由表面的凹陷深度上:解耦模型预测的凹陷深度比实测浅20%左右,而一体化模型和高速摄像测量值基本吻合。这个差异就是因为解耦模型没法准确捕捉熔池表面对电弧压力分布的反馈。
1.2 强耦合机制:热、力、电的三方博弈
电弧和熔池之间的交互可以归纳为热、力、电三个通道的博弈,任何一个通道的变化都会牵动另外两个。
热通道:电弧把焦耳热和电子携带的凝结热输入熔池,熔池表面温度升高后产生蒸发,蒸发带走热量,金属蒸气又反过来改变电弧的热导和辐射特性。热通道的耦合直接决定熔深和熔宽。
力通道:电弧对熔池表面施加压力(电弧静压力)、剪切力(等离子体流沿熔池表面切向流动产生)以及熔池金属蒸发产生的反作用力。这几个力决定熔池表面怎么凹陷、熔池内部怎么流动。TIG焊中电弧压力驱动熔池表面下凹,形成"熔池倒流"现象,直接关联成形质量。
电通道:电流从电弧进入熔池,再穿过工件。熔池表面的温度、成分、几何形状会影响电流密度分布,电流密度分布又决定焦耳热的落点和洛伦兹力的分布。MIG焊中,熔滴在洛伦兹力作用下被"收腰"加速,这是射流过渡形成的核心机制之一。
这三条通道不是各管各的,而是实时纠缠。模拟中最容易犯的错误就是只耦合热和力,忽略电通道的反馈,导致熔滴脱离频率算不准。
2. TIG焊:电弧加热熔池,熔池反过来捏住电弧
TIG焊的一体化建模,是所有焊接模拟里相对"友好"的起点,因为不涉及熔滴过渡,几何结构相对简单,焊丝和钨极位置稳定。但简单的几何不代表物理简单,至少有三个地方需要认真处理。
2.1 电弧-熔池界面的物理边界
TIG电弧和熔池的界面,是整个模型的灵魂位置。这个界面不是一条固定边界,而是一个可以上下变形、温度极高、有蒸发和质量交换的自由表面。
在实际模拟中,这个界面上的热交换分三块:电弧输入给熔池的热、熔池表面通过辐射和对流散掉的热、熔池表面蒸发带走的热。电弧输入热又分为两部分:弧柱传导和对流过来的热,以及电子在阳极附近凝结时释放的凝结热。阳极压降效应让电子在接近阳极表面时被加速获得额外动能,碰撞到熔池表面后动能转化为热能。这一部分的功率密度不可小觑,在高电流密度条件下,阳极凝结热可以占总输入热量的三成以上。
界面上的力平衡同样复杂。电弧压力向下压熔池表面,熔池内部的浮力和表面张力向上撑住,重力向下施加。三力平衡的微小偏移就会造成熔池液面高低变化。模拟时通常用压力平衡方程来确定自由表面的位置,而不是人为给定熔池形状。我做过一个对比案例:在200A电流下,如果把熔池表面强制设为平面,计算出的熔深比实测低了将近15%,只有让表面自由变形,熔深才接近实测值。
界面附近还有强烈的蒸发。钨极TIG焊熔池温度通常在1800到2200开尔文范围,这时低碳钢的蒸发通量已经可观,蒸发的铁原子会进入电弧边缘区域,改变电弧边缘的电导率分布,影响着电弧径向的热流分布。别小看这个效应,很多高速摄像里看到的电弧边缘亮度波动,源头就是蒸发。
2.2 一体化模型怎么预测熔深和熔池形貌
熔深预测是TIG模拟的金标准,算得准熔深,才有资格谈工艺优化。一体化模型预测熔深的关键,在于准确处理电流密度分布和熔池内部的流动。
熔池内部不是静止的液体。焊接电流沿熔池表面进入液态金属后,在电弧中心区域电流密度最大,电流线从表面中心向下流散,这个发散过程产生朝下的洛伦兹力,驱动液态金属沿中心向下流动,然后沿熔池底部向边缘翻涌,形成经典的"外向环流"。这个环流把高温金属不断送到熔池底部,熔池也跟着变深。如果换成表面活性元素(比如硫含量高的钢),熔池表面张力温度系数反转,液流方向反向,熔池就会变宽变浅。这类"小元素改大形貌"的效应,只有在一体化模型里才会自动浮现,因为表面张力边界是直接耦合在自由表面方程里的。
我做过的典型算例是304不锈钢TIG焊,电流从120A扫到250A,每组都记录熔深、熔宽和熔池表面凹陷量。模拟结果在低速区和实验中段吻合良好,误差在±8%以内。熔池表面的中心凹陷深度在120A时约0.6毫米,250A时达到1.4毫米左右,这个量级与实测吻合。凹陷深度直接影响电弧的等效弧长,弧长又改变电弧电压和热输入,这个闭环只有一体化模型能捕捉。
2.3 参数变化下熔池流动的"指挥棒"
在TIG一体化模型里,有几个参数对整个体系影响最大,调试时需要优先关注。
焊接电流:电流不只是加热功率,它决定洛伦兹力的大小,洛伦兹力又是熔池环流的主要驱动力。电流升高时,熔池环流速度显著增大,从几十厘米每秒可以跳到每秒数米量级,熔池形态也随之变化。
弧长:弧长改变的是电弧电压和热流分布。弧长太短,电弧被熔池表面的蒸发气体强烈压窄,热流集中,容易咬边;弧长太长,热流分散,熔深浅。一体化模型里弧长是几何输入参数,但真实影响是通过计算得到的,不是经验修正。
保护气体:氩气、氦气、混合气的电离特性、导热系数不同,直接影响电弧温度和热流分布的形态。氦气电弧温度高、弧柱宽,适合传热需求大的工况。模拟时要注意气体物性随温度和成分的变化,不能取常温常数。
钨极几何:钨极尖端角度和直径影响电流密度初始分布,尖端角度越大电弧收缩越剧烈。这些细节模拟得很到位时,连钨极烧损对电弧形状的间接影响都可以纳入研究范围。
3. MIG焊:熔滴过渡和金属蒸气才是真正的主角
MIG焊的一体化模型比TIG难一个量级,因为计算域里多了一个动态变化的熔滴,熔滴的形态、位置、速度每时每刻都在变,而且和电弧等离子体、金属蒸气浓度场互相作用。
3.1 熔滴过渡模式与电弧的协同演化
MIG焊经典的过渡模式划分——短路过渡、滴状过渡、射流过渡、旋转射流过渡——本质上对应不同的电弧与熔滴作用平衡状态。
滴状过渡阶段,熔滴尺寸较大,脱离频率低,熔滴在焊丝末端长大到一定尺寸后,重力超过表面张力支持力才脱落。这个阶段电弧相对稳定,熔滴下落的路径由重力主导。模拟中需要捕捉熔滴颈部收缩、发生缩颈时电流集中加热导致颈部温度急剧升高、最终熔断的过程,这个过程如果转成动画看,是非常直观的"电-热-力耦合"课。
射流过渡阶段,电弧在熔滴下方呈锥形包住熔滴,熔滴被电磁力向下挤压,顶部不断有液体被剥落成细小颗粒。这个阶段熔滴尺寸小、频率高,模拟需要非常细致的网格分辨能力才能捕捉熔滴的周期性形成和脱落。射流过渡的临界电流不是固定值,它与焊丝直径、保护气体成分、焊丝表面状态都有关系,模拟的价值在于可以针对具体材料和焊丝规格找到这个临界点,而不必做大量实验扫描。
短路过渡阶段难度最大,熔滴接触熔池后形成液桥,电流通过液桥产生电阻热,液桥温度迅速升高并爆破,产生飞溅。这个阶段熔滴、熔池、电弧三者位置关系发生剧烈变化,液桥爆破的瞬态过程尺度在毫秒级,网格和时间步长都要足够细才能捕捉。
3.2 金属蒸气:看不见的"第三只手"
MIG焊的金属蒸气效应,是近十年电弧物理模拟里讨论最多、影响也最大的现象。焊丝和熔池蒸发出的铁蒸气、锰蒸气、铬蒸气等,进入电弧后,会显著改变电弧特性。
金属蒸气的电离能比氩气电离能低很多,这意味着在同样的电离条件下,金属蒸气更容易电离,电弧中心区域电导率会大幅上升。电导率上升的后果是电弧中心电流通道变窄,电流密度集中,于是电弧表现为"收缩"。很多人直观认为蒸气多会让电弧变亮变热,实际上弧心温度可能因为大电流通道集中在极窄区域反而有所下降,但电流密度和热流密度却升高了。
蒸气还有一个重要的动力学效应:蒸发出来的金属原子在熔滴和熔池表面形成的蒸气反冲压力,对熔滴表面产生向外的推力,这个力可以在一定程度上"托住"熔滴,推迟熔滴脱落。早期做MIG模拟忽略蒸气反冲力,算出的熔滴尺寸总是偏小、过渡频率总是偏大,原因就在这里。
我在模拟低碳钢MIG焊时,分别做了带金属蒸气模型和不带金属蒸气模型的对比。不带蒸气模型时,200A电流下电弧温度场中心温度稳定在约22000开尔文;加上蒸气模型后,因蒸气浓度增高导致的辐射冷却和电导率变化,中心温度下降到约14000开尔文。这个差异直接决定了对熔滴加热程度的判断。从实测电弧光谱反推的温度来看,带蒸气模型的结果明显更接近真实。
3.3 飞溅到底是怎么来的
MIG焊的飞溅问题是工艺现场最头疼的问题之一,一体化模型能直接揭示飞溅的机理。飞溅主要有两类来源。
第一类是短路过渡的液桥爆破飞溅。当熔滴与熔池接触形成液桥后,液桥截面突然变窄,电流密度飙升,液桥在极短时间内被电阻热加热到汽化点,金属剧烈蒸发产生的高压将周围液态金属炸飞。模拟中可以在液桥直径最小处捕捉到温度急剧升高、压力波向外扩张的瞬态过程,这个位置和时间点,就是飞溅的"起爆点"。减少这类飞溅的关键是控制液桥爆破前的电流上升速率,这也是很多波形控制焊机设计控制策略的理论依据。
第二类是射流过渡中熔滴在电弧力作用下被撕碎时产生的小颗粒飞溅。当熔滴受到过大的等离子体流力和电磁力时,液态金属表面会形成细微凸起,凸起被气流撕裂成小颗粒高速飞出,形成飞溅。模拟中表现为熔滴边缘周期性产生小液滴并向外偏转。这类飞溅与保护气体流量、电流波形参数密切相关。
把飞溅机理吃透之后,再去调整焊接参数就是有的放矢了。比如短路过渡中以较低的电流增长率起步,让液桥有时间慢慢缩颈而不是瞬间爆断;射流过渡中避免电流过冲,控制电磁力的峰值。这些策略不是拍脑袋想出来的,全是模拟结果推导出来的。
4. 数值建模的完整流程(可直接抄作业)
说完了物理图景,进入实操环节。下面这套流程,是我在若干模拟项目里反复打磨出来的标准路线,按这个顺序走,踩坑概率会小很多。
4.1 控制方程与源项处理
一体化模型的核心方程组是典型的磁流体力学(MHD)加自由表面流动问题。气流区域和熔池区域共用一套连续介质假设下的Navier-Stokes方程组,但物性参数随温度和成分剧烈变化。
质量守恒方程、动量方程、能量方程、电流连续方程是四大支柱。动量方程里需要加入洛伦兹力源项,公式就是电流密度矢量叉乘磁感应强度矢量。能量方程里需要加入焦耳热源项、辐射损失项和电子凝结热边界源项。电流连续方程求解出电标量势,再从电场强度计算电流密度,磁场的获取通常通过磁矢势法求解。
这里有个关键细节:辐射损失不能简单用一个常系数近似。电弧温度从5000到20000开尔文范围内,辐射损失系数随温度变化非常剧烈。工程常用的做法是查表插值,或者用净辐射系数模型,这个系数在不同气体成分(纯氩还是氩气加蒸气)下差异很大。如果偷懒用常数辐射系数,电弧中心温度就会明显偏高,后续熔池热量输入的准确性就无从谈起。
自由表面的处理,TIG模型和MIG模型也要分开考虑。TIG熔池表面变形幅度相对温和,用Level-Set方法或者VOF方法都可以稳定求解。MIG熔滴过渡时表面拓扑变化剧烈,液桥的拉伸、断裂、合并,只能用VOF方法配合几何重构(Geo-Reconstruct)才能稳健处理。用Level-Set勉强做MIG熔滴滴落,大概率会出现界面质量不守恒的尴尬。
4.2 自由表面、蒸发和边界条件设置
边界条件是一体化模型最容易出错的地方,也最考验经验。我把关键边界列成下面这张表:
| 边界位置 | 主要设置 | 容易踩的坑 |
|---|---|---|
| 钨极/焊丝表面 | 电流入口边界,给定电流密度,热边界可设对流冷却 | 电流密度分布给定方式不同,弧根位置会漂移 |
| 电弧外围 | 压力出口或开放边界,保护气体入口按流量给定 | 出口边界太近会约束电弧膨胀,导致温度场失真 |
| 熔池表面 | 自由表面,耦合热流、剪切力、蒸发通量 | 蒸发通量系数太大会导致熔池表面温度异常降低 |
| 工件底面和侧面 | 热对流和辐射边界,底面通常固定为环境温度或给定换热系数 | 底面换热系数对熔深结果的影响比很多人想的敏感 |
| 对称面 | 轴对称模型时的对称边界 | 三维模型没有这个问题,但计算量翻倍 |
蒸发通量的处理,推荐用Langmuir蒸发方程结合环境压力的反压力修正。实际蒸发的难点是蒸发系数(蒸发概率)并不等于理论值,需要根据实验校准,通常的取值区间在0.01到1之间。我在低碳钢熔池模拟中校准过几组,发现系数取0.05附近时,熔池表面温度和蒸发冷却效果与热电偶测量最接近。这个系数的敏感性很高,是模拟结果不确定性的主要来源之一,之前提过,它值得你花时间反复标定。
4.3 网格划分与求解策略
网格划分直接决定一体化模型的成败。电弧弧柱区域温度梯度和速度梯度都大,网格要足够密;熔池区域自由表面附近必须要有精细边界层网格来分辨流体剪切作用;工件内部则可以用较粗的网格节省算力。
我的经验值供参考:TIG轴对称模型,电弧区域最小网格尺度取0.1毫米左右,熔池表面沿向网格细化到0.05毫米;MIG三维模型,焊丝和熔滴附近最小网格尺度取0.05毫米,总网格量通常在几十万到几百万之间。如果网格尺度超过0.3毫米,电弧边缘的温度场和电流密度分布就会明显失真,熔滴缩颈位置的捕捉也会出问题。
求解策略上,强烈建议先稳态后瞬态。第一步,固定熔池表面为平面,只解电弧温度场和流场,得到初始电弧状态;第二步,开启自由表面变形,让熔池参与耦合;第三步,电流从小值逐步加载到大值,避免一次性大电流导致初期振荡发散。MIG模型还需要额外处理焊丝进给:通常用质量入口边界条件模拟焊丝连续送入,在焊丝末端累积的金属形成熔滴。瞬态计算的时间步长从0.1到10微秒之间按需调整,液桥断裂前后要加密到微秒量级。
5. 典型模拟结果怎么读、怎么用
模拟算完不等于工作结束,读懂结果、把模拟结论转化成工艺语言,才是投入产出的关键。
5.1 从温度场看过渡模式转变
MIG模型的温度场云图里,最容易识别的是熔滴过渡模式转变。滴状过渡阶段,电弧呈现典型的钟形,熔滴宿主在焊丝底端缓慢长大,温度等值线在熔滴周围形成一层明显的"热包裹"。射流过渡阶段,电弧顶部收窄、底部扩展成锥形,熔滴顶部不断被剥落成密集细颗粒,温度场呈现明显的高频脉动。
我常用一个指标来定量判断过渡模式:熔滴脱离频率与熔滴直径的乘积。滴状过渡时这个值稳定在0.3到0.5之间,射流过渡时跳到2以上。模拟结果里看这个指标的变化趋势,再结合高速摄像验证,就能很可靠地区分过渡模式转变点。这个转变点对应的电流值,就是实际焊接工艺窗口的下限参考。
5.2 蒸气浓度与电弧压力分布的关系
金属蒸气的浓度场分布,在结果后处理中经常被忽视,但它实际上解释了电弧行为的很多"异常"现象。比如模拟中发现电弧电压低于经验公式预测值,检查一下蒸气浓度场,往往能看到弧心区域的蒸气浓度已经超过10%,因为蒸气电离能低,同样电流下维持电弧所需电压自然降低。
电弧压力分布直接影响熔池表面凹陷和飞溅倾向。读取模拟结果里熔池表面中心的电弧压力峰值,正常TIG在200A电流下大概在200到500帕之间。如果算出来压力峰值异常偏高,就要检查是不是忽略了蒸气反冲力或者边界条件有误。MIG焊中电弧压力分布不再对称,熔滴下落一侧压力明显更大,这种不对称是熔滴偏转和飞溅的重要驱动因素。
5.3 工艺参数窗口的快速筛选
一体化模拟最重要的工程价值,就是代替部分实验来做参数窗口的快筛。传统方法每调整一个参数组就要做焊接试验加金相检测,一轮下来至少半天。模拟虽然单组算得慢,但可以并行跑多组参数,而且可以得到温度、流场、蒸汽浓度等实验中难以直接测得的中间量。
我在实际项目中做过这样一件事:针对某高强钢MIG焊的飞溅率问题,用模拟扫了电流、电压、干伸长三因素共27组参数组合,从中筛出飞溅率预估最低的5组再去实验验证。实验结果显示,模拟预估的飞溅率排序与实验实测排序高度一致,前3组里有两组达到了现场工艺要求。这比盲目做实验试错效率高太多。
做这类参数扫掠时,建议每组模拟保持标准焊接位置(比如平焊),只改变目标参数。否则位置因素和参数因素混在一起,结果没法归因。
6. 我踩过的坑和排查技巧
一体化模拟经常出现看似正常实则偏离物理的结果。分享几个我踩过多次的坑和处理思路。
6.1 收敛困难的排查顺序
模型发散或振荡时,别急着调小松弛因子。先按顺序排查以下三类问题。
第一,边界条件是不是自洽。电流入口面积、电流密度值、总电流三者的乘积关系对不对。如果焊丝直径是1.2毫米,电流密度按200A输入,入口半径稍微大了0.1毫米,总电流偏差就超过10%,整个电磁力场就开始偏离物理。
第二,物性参数表在温度区间内是否连续且有物理意义。高温下材料热导率和比热的插值外推,是最容易在计算中产生虚假能量峰的地方。我吃过一次大亏:某钢材热导率在3000开尔文以上外推时出现了不合理的下降趋势,模拟熔池深度的偏差高达30%,一直查不出来,最后逐行核对物性表才发现。
第三,网格质量在自由表面附近是否足够。界面附近网格纵横比过大会导致压力场震荡,这在VOF类方法里尤其常见。遇到振荡先用瞬间的动态压力分布图定位振荡区域,看是不是集中在熔池表面附近。
6.2 界面模糊、蒸发敏感、材料属性的大坑
界面模糊是VOF方法的经典难题。熔滴表面被"糊"出过渡层,厚度超过两个网格,导致界面处表面张力计算不准确,熔滴脱离时间偏移。解决思路主要有:提高网格分辨率、使用界面压缩格式(如界面压缩项)、以及在重构时选择更高阶的几何重构格式。这三个手段组合使用,通常可以把界面模糊控制在1到2个网格尺度以内。
蒸发模型的敏感性前面已经反复强调过。这里再补一个排查技巧:如果你发现模拟中熔池表面温度一直上不去,比如总在1850开尔文上下徘徊,而实测熔池表面可以达到2000开尔文以上,大概率就是蒸发系数设置过大,蒸发冷却吸走了太多热量。反过来,如果熔池表面温度明显高于实测,且电弧亮度异常集中,可能是蒸发系数太小或者完全没开蒸发模型。蒸发系数这个参数值得你专门做一组标定实验,用不同电流下的实测熔池表面温度曲线去反推最合理的取值。
材料属性方面,还要注意液态金属的表面张力温度系数。这个参数的正负号直接改变熔池环流方向。低碳钢中硫含量从0.008%到0.05%的变化,可以反转表面张力温度系数的符号,导致熔池从纵深向变为宽浅向。模拟时如果对材料成分没有把握,应当做表面张力温度系数的区间敏感性分析。
6.3 几个小技巧与建议
最后分享几个实战技巧。
第一个技巧关于时间步长的自适应。熔滴过渡的瞬态尺度在毫秒级,但熔池热传导的时间尺度在秒级。如果全程用小步长,算到熔池稳定要等太久。我习惯采用双时间尺度策略:液桥断裂、熔滴分离等剧烈事件前后使用微秒级步长,熔池缓慢演化阶段使用毫秒级步长。具体判断标准是界面最大移动速度乘以时间步长不超过网格最小尺度的四分之一。
第二个技巧是重视后处理中的守恒性检查。每算完一组,都检查一下总电流是否等于边界输入的电流值,总能量收支的残差是否在合理范围。守恒性问题往往是建模疏漏的早期信号,比看温度场云图更可靠。
第三个技巧有关模型降阶。完全一体化模型的计算成本仍然偏高,不适合做大规模参数寻优。我常用的策略是先用一体化模型建立若干典型工况的基准数据库,再基于基准数据训练一个简化模型——比如用响应面或者神经网络来建立熔深、飞溅率与关键工艺参数的关系。这样既保留了一体化模型的物理精度,又能满足工艺现场快速响应的需求。
我个人在实际操作中的体会是,一体化模拟的另一个隐性收益在于帮助建立"耦合直觉"。当你反复观察过熔滴脱落瞬间的温度场重分布、熔池表面凹陷对电弧压力的影响、蒸气浓度变化对电弧电压的调制之后,处理现场问题时会多一个思维维度。飞溅大了先想液桥爆破时机,熔深浅了先看熔池环流方向,这些判断在纸上谈兵时很难建立,但模型会逼着你把因果关系捋清楚。后面如果再遇到电弧不稳或者成形缺陷,顺着热、力、电三条线去分析,很多问题会清晰很多。