本发明专利技术公开了一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法，针对现有焊缝形貌预测方法无法考虑熔池与羽辉耦合行为进而导致万瓦级激光焊接焊缝形貌预测精度低的问题，本发明专利技术采用基于压力的可压缩两相流数值计算方法求解熔池与羽辉耦合行为，从而实现万瓦级激光焊接焊缝形貌的高精度预测。首先，输入初始时刻的焊缝形貌函数和焊接参数；其次，采用基于压力的可压缩两相流数值计算方法获取下一时刻的焊缝形貌函数；进而，绘制下一时刻焊缝形貌函数并提取焊缝形貌及焊缝形貌特征。与现有焊缝形貌预测方法相比，本发明专利技术可准确计算万瓦级激光焊接中熔池与羽辉耦合行为，算法简单容易实施，计算效率高，物理守恒性好，可实现对万瓦级激光焊接焊缝形貌的高精度预测。对万瓦级激光焊接焊缝形貌的高精度预测。对万瓦级激光焊接焊缝形貌的高精度预测。
【技术实现步骤摘要】
一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法
[0001]本专利技术涉及一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法，其属于焊缝形貌预测方法领域。
技术介绍
[0002]焊缝形貌是决定万瓦级激光焊接成形质量和成形效率的关键因素。准确地预测焊缝形貌，可以有效地减少万瓦级激光焊接工艺实验，降低实验成本，提高企业研发和生产效率。
[0003]目前焊缝形貌预测主要采用数值模拟方法，其原理是利用数值算法求解焊缝形成过程中的传热、流动及自由表面演化等方程，进而实现对焊缝形貌变化过程的有效预测。然而在万瓦级激光焊接过程中，羽辉会产生高速的流动(马赫数接近甚至大于1)形成可压缩特性。这种特性会使羽辉压力场产生显著变化(最大压力差超过1个大气压)，进而强烈改变焊缝形貌的变化过程。
[0004]现有焊缝形貌预测方法如不可压缩解法(Ai Y,Jiang P,Shao X,et al.The prediction of the whole weld in fiber laser keyhole welding based on numerical simulation[J].Applied Thermal Engineering,2016,113:980
‑
993)无法准确描述这种可压缩特性。而基于密度的(Density
‑
based)可压缩解法(Amara E H,Fabbro R,Bendib.Modeling of the compressible vapor flow induced in a keyhole during laser welding[J].Journal of Applied Physics,2003,93:4289
‑
4296)在求解熔池低速流动时(马赫数小于0.3)又会产生压力场振荡，进而给数值计算的实施带来困难。有学者提出不可压缩+可压缩的组合方法(Pang S,Chen X,Li W,et al.Efficient multiple time scale method for modeling compressible vapor plume dynamics inside transient keyhole during fiber laser welding[J].Optics&Laser Technology,2016,77:203
‑
214)，此方法需要在熔池和羽辉相界面处构造不同的数值通量和人工边界条件，会带来物理场不守恒和预测精度低等问题。上述难点制约了焊缝形貌预测方法在万瓦级激光焊接中的应用，因此迫切需要发展一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法。
技术实现思路
[0005]为解决万瓦级激光焊接中焊缝形貌预测精度低的问题，本专利技术提出了一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法。
[0006]本专利技术采用如下技术方案：一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法，包括如下步骤：
[0007]S1,输入初始时刻的焊缝形貌函数H
n
：
[0008][0009]其中，公式(1)中，Δx
min
为最小网格尺寸，φ
n
为初始时刻工件截面内各点到焊缝
的距离，tanh为双曲正切函数；
[0010]S2,确定初始时刻的焊接参数，各焊接参数表达如下：
[0011][0012]S3,获取下一时刻的焊缝形貌函数H
n+1
，获取过程如下：
[0013]S3.1,计算初始时刻的动量源S
n
：
[0014][0015]其中，公式(2)中，分别代表由惯性、表面力和体积力引起的动量源，为梯度算子，
·
为散度算子，上标
T
代表矩阵转置，g为重力加速度常数，k
b
为玻尔兹曼常数；
[0016]S3.2,计算初始时刻的能量源Q
n
：
[0017][0018]其中，公式(3)中，r为各点到激光光斑中心的距离；
[0019]S3.3,计算对流后的密度ρ
*
、速度u
*
和温度T
*
：
[0020][0021]其中，公式(4)中，Δt为计算时间步长；
[0022]S3.4,计算下一时刻的压力p
n+1
：
[0023][0024]其中，公式(5)中，R为理想气体常数；
[0025]S3.5,计算下一时刻的密度ρ
n+1
、速度u
n+1
、温度T
n+1
：
[0026][0027]S3.6,计算下一时刻的焊缝形貌函数H
n+1
：
[0028][0029]S4,绘制下一时刻焊缝形貌函数H
n+1
的分布云图；
[0030]S5,提取下一时刻的焊缝形貌；
[0031]S6,输出下一时刻的焊缝形貌特征。
[0032]与现有方法相比，本专利技术的优点如下：
[0033](1)与不可压缩解法相比，本方法考虑了羽辉可压缩性带来的压力场变化，可准确计算出熔池与羽辉的耦合行为，从而实现万瓦级激光焊接焊缝形貌的高精度预测；
[0034](2)与可压缩解法相比，本方法将压力作为控制方程的直接解，压力计算精度高，计算稳定性好；
[0035](3)与不可压缩+可压缩的组合方法相比，本方法采用一套控制方程求解熔池和羽辉两种流体，具有算法简单容易实施、计算效率高、物理守恒性好等优点。
附图说明
[0036]图1是本专利技术提出的适用于万瓦级激光焊接的高精度焊缝形貌预测方法的流程过程图。
[0037]图2是初始时刻的焊缝形貌函数定义的示意图。
[0038]图3是下一时刻焊缝形貌函数分布云图的示意图。
[0039]图4是下一时刻的焊缝形貌示意图。
[0040]图5是下一时刻的焊缝形貌特征示意图。
具体实施方式
[0041]下面结合附图和实施例对本专利技术进行进一步阐述。
[0042]针对现有焊缝形貌预测方法无法考虑熔池与羽辉耦合行为进而导致万瓦级激光焊接焊缝形貌预测精度低的问题，本专利技术采用基于压力的可压缩两相流数值计算方法求解熔池与羽辉耦合行为，从而实现万瓦级激光焊接焊缝形貌的高精度预测。
[0043]结合图1，本专利技术所述的一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法，包括如下步骤：
[0044]第一步：输入初始时刻的焊缝形貌函数
[0045]本专利技术以双曲正切函数表示焊缝形貌H
n
，焊缝形貌定义为双曲正切函数的0.5等值面。初始时刻的焊缝形貌函数H
n
见公式(1)：
[0046][0047]其中，Δx
min
为最小网格尺寸，本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法，其特征在于，包括如下步骤：S1,输入初始时刻的焊缝形貌函数H
n
：其中，公式(1)中，Δx
min
为最小网格尺寸，φ
n
为初始时刻工件截面内各点到焊缝的距离，tanh为双曲正切函数；S2,确定初始时刻的焊接参数，各焊接参数表达如下：S3,获取下一时刻的焊缝形貌函数H
n+1
，获取过程如下：S3.1,计算初始时刻的动量源S
n
：其中，公式(2)中，分别代表由惯性、表面力和体积力引起的动量源，为梯度算子，为散度算子，上标
T
代表矩阵转置，g为重力加速度常数，k
b
为玻尔兹曼常数；S3.2,计算初始时刻的能量源Q
n
：
其中，公式(3)中，r为各点到激光光斑中心的距离；S3.3,计算对流后的密度ρ
*
、速度u
*
和温度T
*
：其中，公式(4)中，Δt为计算时间步长；S3.4,计算下一时刻的压力p
n+1
：其中，公式(5)中，R为理...
【专利技术属性】
技术研发人员：钱龙根，秦增蓉，朱宇，
申请(专利权)人：常熟理工学院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。