高能束焊数值模拟中热源模型参数的确定研究.pdf

高能束焊数值模拟中热源模型参数的确定 清华大学机械工程系 王煜赵海燕吴照张建强刘点兵 摘要在高能束焊的数值模拟中，热源模型参数的选择是否合适，决定了计算结果的准确与否． 以往一般借助经验，通过反复试算决定模型参数，造成计算效率低下．本文提出了一种用解析法计算 高能柬焊双椭球热源模型参数的方法，并以实例计算验证了方法的可行性。应用此方法可以方便的确 定实际计算所需的热源参数，简化试算过程，提高计算效率． 关键词：高能束流焊 双椭球热源模型解析法 0前言 高能束流焊接是一种重要的焊接方法，其焊接的功率密度(Power Density)可达到 105 w，cm2以上，属于高功率密度的热源。束流由单一的离子、电子、光子或二种以 上的粒子组合而成，如等离子弧焊、电子束焊、激光束焊及复合热源焊(激光束+Arc 焊后成品热影响区小，焊缝质量高，所以高能束焊得到了日益广泛的应用。为了提高生 产效率和产品质量，对焊接过程进行精确控制是必要的。这需要对焊接过程进行深入详 细的研究。值由于金属熔化凝固时温度高，而且又是一个动态过程，因此采用实验方法 进行实时测量十分困难。在这方面，数值模拟为我们提供了重要的研究手段。 对焊接热过程的准确模拟是确保热应力变形分析可靠性的重要前提，为了准确计 算，针对不同的加热过程，应采用不同的热输入模式(热源模型)进行计算。由于高能 束流的能量密度高，加热范围小，加热，冷却速度快和存在锁孔效应等特点，用于普通熔 化焊计算的Gauss面热源模型不适于描述这一过程，而采用能够反映沿深度方向加热特 点的双椭球形体热源模型进行模拟。 1双椭球热源模型‘21 双椭球模型所描述的热流输入分布在一定的体积内。模型考虑了焊接中的“挖掘” “搅拌”作用，适用于描述高能柬焊如电子束焊、激光焊等深宽比较大的深熔焊接过程 的热输入。如图1所示： 图1 双椭球体热源模型 模型内部热流输入分布为： 椭球沿Y轴前半部分 l lO (1) g(x，y，z，，)：．6ff3f／Q_．-3x2／a2e-322Ib2e-3㈨{r．，v)12／砰 abc，r4tc 2P31y+v(r-t)1饥i 己 黼沿乙嚣篙竺蜒e_3．d／,,2C_322／h 式中： abcrc4z a，b，cl，C2：热源形状参数lm’I Q：热输入功率【Wl “，fr：前后部分的能量分布参数 模型分成前后长度不同的两部分，是为了能更好的反映出焊接过程中移动热源的 前端和后端不同的温度梯度分布(前端较陡，后端较缓)。在一般计算时，作为简化处 理，取cI=c2=c。热源可由统一方程描述为： r一，)12_： g(x，．y，z，，)：粤P一3—2“2P-3=／hc-31y+、，I (3) 由于热流分布在椭球形体积内，能够反映出高能束流沿深度方向对焊件进行加热的特点， 因此对于焊接温度场可以进行更为准确的模拟。但双椭球模型中四个形状参数a、b、c．、C2 一般由经验给定，没有定量化的公式描述。这对于实际应用是非常不便的。在计算温度场时， 评定计算结果是否合理的判据之一是能否获得较准确的焊缝熔池形状。对于功率同样大小的 热源，当功率密度分布不同时，所得的熔池形状会有很大差异。在使用双椭球模型进行计算 时，常常需要经过多次试算，观察当焊接进入稳态时的熔池尺寸，根据情况调整形状参数， 直至能够获得合理形状尺寸的熔池。此时才能将调整好的热源模型应用于实际的温度场或热 力耦合计算中。而进行有限元运算通常是一项非常耗时的计算过程，而为获得可用热源模型 所进行的反复试算会消耗大量的时问，造成工作效率低下