CN109933007A - 一种激光-mig复合焊接热源模型的建立方法 - Google Patents
一种激光-mig复合焊接热源模型的建立方法 Download PDFInfo
- Publication number
- CN109933007A CN109933007A CN201910090055.9A CN201910090055A CN109933007A CN 109933007 A CN109933007 A CN 109933007A CN 201910090055 A CN201910090055 A CN 201910090055A CN 109933007 A CN109933007 A CN 109933007A
- Authority
- CN
- China
- Prior art keywords
- heat source
- laser
- welding
- mig
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Abstract
本发明公开了一种激光‑MIG复合焊接热源模型的建立方法,属于激光焊接技术领域。该方法首先建立对接接头激光‑MIG复合焊接三维有限元网格模型,然后基于FORTRAN语言,开发复合热源模型,采用有限元模拟软件计算求解控制方程,进行温度场计算,即完成了激光‑MIG复合焊接热源模型的建立。利用本模型可获取不同焊接工艺参数下的温度场,对准确认识金属材料激光‑MIG复合焊接温度场的分布规律、控制焊接质量和使用性能提供技术支持。
Description
技术领域
本发明涉及激光焊接技术领域,具体涉及一种激光-MIG复合焊接热源模型的建立方法。
背景技术
激光焊接作为一种精密、高效、快速的高能束焊接方法,在汽车、船舶、航空航天等领域得到了越来越广泛的应用。但常规的激光焊接存在着激光能量利用率低、装配精度要求高等问题。为了解决这些问题,激光-熔化极惰性气体保护焊(Metal Inert GasWelding,以下简称MIG)复合焊接技术受到了广泛关注。激光-MIG复合焊将激光焊与MIG焊的优点进行了有机结合,不仅提高了激光能量利用率和焊接过程的稳定性,改善焊缝质量,而且降低了对焊接件装配精度的要求,同时可弥补MIG焊存在的焊速慢、效率低、变形大等不足。目前,激光-MIG焊接技术已应用于德国大众、Aker Warnow Werf船厂等汽车和船舶的制造中。由于工业结构件往往形状复杂、焊接道次多,制造工艺不成熟,如果采用物理实验或者依靠经验的方法需要做多次试验,周期长,成本高,严重影响结构件的生产制造。采用数值模拟的方法来代替前期不必要的物理实验,不仅可以节约大量的成本和人力物力资源,而且可以缩短产品研制周期,提高产品竞争力。
建立合理的热源模型是焊接过程数值模拟结果准确、可靠的前提。因此,激光-MIG复合焊接过程模拟首要解决的问题是建立合理的激光-MIG复合焊接热源模型。在激光和MIG焊热源的双重作用下,焊缝截面形貌呈“倒喇叭”型,这种形貌难以采用有限元软件中自带的高斯面热源模型和双椭球体热源模型来描述。所以,要实现激光-MIG复合焊接过程模拟,必须建立合适的焊接热源模型,以便优化焊接工艺。
发明内容
本发明的目的在于提供一种激光-MIG复合焊接热源模型的建立方法,以解决现有有限元软件自带的热源模型无法描述“倒喇叭”型焊缝截面形貌的问题。
为实现上述目的,本发明所采用的技术方案如下:
一种激光-MIG复合焊接热源模型的建立方法,该方法包括如下步骤:
第一步,建立对接接头激光-MIG复合焊接三维有限元网格模型:在三维造型软件中建立几何模型,利用网格划分软件或者有限元模拟软件对工件进行网格划分;网格划分时在焊缝及其周围区域进行网格细化,以减少误差。
第二步,建立三维激光-MIG复合焊热源模型,基于FORTRAN语言,开发组合热源模型;其中:
工件的上半部分采用双椭球体热源模型,控制方程如公式(1)~(3):
qV1=qf+qr (3);
公式(1)~(3)中:QV1为双椭球体热源占有的能量,af、ar、b、c为双椭球体热源形状参数,H1为双椭球体热源作用高度;qf和qr分别为前、后半椭球体内的热流密度分布;qV1为双椭球体热源热流密度分布;
工件的下半部分采用圆柱体热源模型,控制方程如公式(4):
公式(4)中:QV2为圆柱体热源占有的能量,r0为高斯面热源有效加热半径,H2为圆柱体热源作用高度,其值为H2=d-H1,d为工件厚度;qV2为圆柱体热源热流密度分布;
对于整个激光-MIG复合焊热源来说,总功率Qtotal等于各组合部分(工件的上半部分和下半部分)热源功率之和,如公式(5);
公式(5)中:QV1=χ1Qtotal,QV2=χ2Qtotal,χ1+χ2=1,χ1和χ2分别为上、下部分热源的功率分配系数;
第三步,基于有限元模拟软件求解控制方程,进行温度场计算:
将第二步所建立的激光-MIG复合焊热源模型加载到有限元网格模型上,设置热边界条件,计算求解,完成热源模型建立,获得焊接温度场和焊缝截面形貌。
激光-MIG复合焊接形成的温度场控制方程如公式(6);
公式(6)中:T为温度,t为时间,ρ为材料的密度,cp为材料的比热容,k为材料的导热系数,x、y、z为坐标系中的三向坐标,qV1和qV2为激光-MIG复合热源输入的能量。
本发明的有益效果如下:
本发明提供一种激光-MIG复合焊接热源模型的建立方法,针对对接接头的激光-MIG复合焊接过程,建立了有限元网格模型,采用FORTRAN语言编写复合热源模型程序代码,设置合理的初始条件和边界条件,利用有限元模拟软件进行计算求解。得到的对接接头截面形貌与实验得到的结果吻合良好,表明本发明热源模型适合激光-MIG复合焊接。
附图说明
图1为本发明建立激光-MIG复合焊接体热源模型示意图。
图2为实际焊缝与本发明热源模型模拟的焊缝截面形貌对比。
具体实施方式
以下结合附图详述本发明。
本发明为激光-MIG复合焊接热源模型的建立方法,如图1所示,具体按以下步骤实现:
第一步,建立对接接头激光-MIG复合焊接三维有限元网格模型:在三维造型软件中建立几何模型,利用网格划分软件或者有限元模拟软件对工件进行网格划分,优先划分焊缝及其周围区域;
第二步,建立三维激光-MIG复合焊热源模型,基于FORTRAN语言,开发组合热源模型;其中:
工件的上半部分采用双椭球体热源模型,控制方程如公式(1)~(3);
qV1=qf+qr (3);
其中,QV1为双椭球体热源占有的能量,af、ar、b、c为双椭球体热源形状参数,H1为双椭球体热源作用高度,qf和qr分别为前、后半椭球体内的热流密度分布,qV1为双椭球体热源热流密度分布。
工件的下半部分采用圆柱体热源模型,控制方程如公式(4):
其中,QV2为圆柱体热源占有的能量,r0为高斯面热源有效加热半径,H2为圆柱体热源作用高度,其值为H2=d-H1,d为工件厚度;qV2为圆柱体热源热流密度分布。
对于整个激光-MIG复合焊热源来说,总功率Qtotal等于各组合部分热源功率之和,如公式(5):
公式(5)中,QV1=χ1Qtotal,QV2=χ2Qtotal,χ1+χ2=1,χ1和χ2分别为上、下部分热源的功率分配系数。
第三步,基于有限元模拟软件求解控制方程,进行温度场计算:
将所建立的激光-MIG复合焊热源模型加载到有限元网格模型上,设置热边界条件,计算求解,完成热源模型建立,获得焊接温度场和焊缝截面形貌。
激光-MIG复合焊接形成的温度场控制方程如公式(6):
其中,T为温度,t为时间,ρ为材料的密度,cp为材料的比热容,k为材料的导热系数,x、y、z为坐标系中的三向坐标,qV1和qV2为激光-MIG复合热源输入的能量。
实施例1:
本实施例以5083-H111铝合金对接接头的激光-MIG复合焊接过程为例进行模拟计算,验证本发明的有益效果。
激光-MIG复合焊接实验使用的母材为5083-H111铝合金,试板尺寸为100mm×60mm×6mm,焊材牌号ER5087,焊丝直径为1.2mm。焊前,依次采用丙酮、7wt.%氢氧化钠水溶液和28wt.%硝酸溶液对5083-H111铝合金试板表面进行清洗。使用HL4406D型YAG激光器和Panasonic YD-500AG2电源实施焊接,焊接工艺参数为激光功率P=4kW,焊接速度v=12mm/s,MIG焊的电流80A,电压为20V。焊后,采用线切割切取焊缝中部,经研磨抛光后,采用凯勒试剂腐蚀,利用体式显微镜获得焊缝截面形貌。
第一步:基于Visual Mesh软件进行实体建模和网格划分,靠近焊缝的区域采用单元较小的加密型网格,离焊缝越远,网格尺寸越大,共29980个单元,25437个节点。
第二步:加载本发明所提供的热源模型,相关参数如下表:
第三步:设置工件地面与工作台之间存在传导换热,其他表面的热边界条件为热对流和热辐射。采用Sysweld软件求解焊接热过程,提取焊缝截面形貌。
为了验证模型计算结果的准确性,将计算结果与实际焊缝进行比较,如图2所示,通过比较可以看出,本发明所提供的热源模型得到了与实际焊缝相似度很高的熔池形貌,且尺寸吻合较好。
Claims (4)
1.一种激光-MIG复合焊热源模型的建立方法,其特征在于:该方法包括如下步骤:
第一步,建立激光-MIG复合焊接三维有限元网格模型:在三维造型软件中建立几何模型,利用网格划分软件或者有限元模拟软件对工件进行网格划分;
第二步,建立三维激光-MIG复合焊热源模型,基于FORTRAN语言,开发组合热源模型;其中:
工件的上半部分采用双椭球体热源模型,控制方程如公式(1)~(3):
qV1=qf+qr (3);
公式(1)~(3)中:QV1为双椭球体热源占有的能量,af、ar、b、c为双椭球体热源形状参数,H1为双椭球体热源作用高度;qf和qr分别为前、后半椭球体内的热流密度分布;qV1为双椭球体热源热流密度分布;
工件的下半部分采用圆柱体热源模型,控制方程如公式(4):
公式(4)中:QV2为圆柱体热源占有的能量,r0为高斯面热源有效加热半径,H2为圆柱体热源作用高度,其值为H2=d-H1,d为工件厚度;qV2为圆柱体热源热流密度分布;
对于整个激光-MIG复合焊热源来说,总功率Qtotal等于各组合部分热源功率之和,如公式(5);
公式(5)中:QV1=χ1Qtotal,QV2=χ2Qtotal,χ1+χ2=1,χ1和χ2分别为上、下部分热源的功率分配系数;
第三步,基于有限元模拟软件求解控制方程,进行温度场计算:
将第二步所建立的激光-MIG复合焊热源模型加载到有限元网格模型上,设置热边界条件,计算求解,完成热源模型建立,获得焊接温度场和焊缝截面形貌。
2.根据权利要求1所述的激光-MIG复合焊热源模型的建立方法,其特征在于:第一步中,三维有限元网格模型是针对对接接头的激光-MIG复合焊接过程建立。
3.根据权利要求1所述的激光-MIG复合焊热源模型的建立方法,其特征在于:第一步中,网格划分时在焊缝及其周围区域进行网格细化,以减少误差。
4.根据权利要求1所述的激光-MIG复合焊热源模型的建立方法,其特征在于:第三步中,激光-MIG复合焊接形成的温度场控制方程如公式(6):
公式(6)中:T为温度,t为时间,ρ为材料的密度,cp为材料的比热容,k为材料的导热系数,x、y、z为坐标系中的三向坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910090055.9A CN109933007A (zh) | 2019-01-30 | 2019-01-30 | 一种激光-mig复合焊接热源模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910090055.9A CN109933007A (zh) | 2019-01-30 | 2019-01-30 | 一种激光-mig复合焊接热源模型的建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109933007A true CN109933007A (zh) | 2019-06-25 |
Family
ID=66985412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910090055.9A Pending CN109933007A (zh) | 2019-01-30 | 2019-01-30 | 一种激光-mig复合焊接热源模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109933007A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110640317A (zh) * | 2019-09-19 | 2020-01-03 | 中国科学院兰州化学物理研究所 | 一种镍钛形状记忆合金的激光焊接方法 |
CN113305435A (zh) * | 2021-05-25 | 2021-08-27 | 华中科技大学 | 激光搅拌焊接中气孔抑制的工艺参数优化方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008200691A (ja) * | 2007-02-19 | 2008-09-04 | Daihen Corp | アーク溶接のビード形状シミュレーション装置 |
CN102708237A (zh) * | 2012-04-27 | 2012-10-03 | 上海飞机制造有限公司 | T型接头双激光束双侧同步焊接熔池流场的模拟方法 |
CN103049623A (zh) * | 2013-01-18 | 2013-04-17 | 哈尔滨工业大学 | 一种激光焊接热源模型的建立方法 |
CN103366044A (zh) * | 2013-06-04 | 2013-10-23 | 广东电网公司电力科学研究院 | 一种基于ansys预测cmt焊缝形状尺寸的方法 |
CN103761374A (zh) * | 2014-01-10 | 2014-04-30 | 上海工程技术大学 | 一种激光深熔焊热源模型的建模方法 |
CN105740577A (zh) * | 2016-03-08 | 2016-07-06 | 南京航空航天大学 | 一种Invar钢PMIG摆动焊温度场及变形模拟方法 |
CN106670666A (zh) * | 2016-12-14 | 2017-05-17 | 南京航空航天大学 | 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法 |
-
2019
- 2019-01-30 CN CN201910090055.9A patent/CN109933007A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008200691A (ja) * | 2007-02-19 | 2008-09-04 | Daihen Corp | アーク溶接のビード形状シミュレーション装置 |
CN102708237A (zh) * | 2012-04-27 | 2012-10-03 | 上海飞机制造有限公司 | T型接头双激光束双侧同步焊接熔池流场的模拟方法 |
CN103049623A (zh) * | 2013-01-18 | 2013-04-17 | 哈尔滨工业大学 | 一种激光焊接热源模型的建立方法 |
CN103366044A (zh) * | 2013-06-04 | 2013-10-23 | 广东电网公司电力科学研究院 | 一种基于ansys预测cmt焊缝形状尺寸的方法 |
CN103761374A (zh) * | 2014-01-10 | 2014-04-30 | 上海工程技术大学 | 一种激光深熔焊热源模型的建模方法 |
CN105740577A (zh) * | 2016-03-08 | 2016-07-06 | 南京航空航天大学 | 一种Invar钢PMIG摆动焊温度场及变形模拟方法 |
CN106670666A (zh) * | 2016-12-14 | 2017-05-17 | 南京航空航天大学 | 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110640317A (zh) * | 2019-09-19 | 2020-01-03 | 中国科学院兰州化学物理研究所 | 一种镍钛形状记忆合金的激光焊接方法 |
CN113305435A (zh) * | 2021-05-25 | 2021-08-27 | 华中科技大学 | 激光搅拌焊接中气孔抑制的工艺参数优化方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112149330B (zh) | 一种风电塔筒油封平台焊接残余应力预测、焊接工艺优化方法 | |
CN106513679B (zh) | 一种适用于大型金属零件的激光增材制造方法 | |
CN103246774B (zh) | 数值模拟p92钢管焊接热影响区宽度的方法 | |
CN105740577A (zh) | 一种Invar钢PMIG摆动焊温度场及变形模拟方法 | |
CN103972628B (zh) | 自定位阶梯扭波导及其真空钎焊工艺 | |
CN111310296B (zh) | 一种应用于复杂曲面的变密度金属板材点阵芯体设计方法 | |
CN109933007A (zh) | 一种激光-mig复合焊接热源模型的建立方法 | |
CN105138772A (zh) | 变截面构件电子束焊的有限元模拟方法 | |
Xu et al. | Finite element analysis of residual stress in hybrid laser-arc welding for butt joint of 12 mm-thick steel plate | |
CN106670666A (zh) | 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法 | |
CN103522014A (zh) | 激光组合焊接厚壁金刚石钻头制造方法 | |
CN111695256B (zh) | 一种基于能量分配系数的激光电弧复合热源的建模方法 | |
CN111680442B (zh) | 一种基于坡口压缩系数的激光电弧热源模型的建模方法 | |
CN112883623A (zh) | 一种Fe-36%Ni合金PMIG摆动焊热-机械耦合有限元建模方法 | |
CN112016228B (zh) | 一种基于水冷压缩系数的水下焊接热源模型的建模方法 | |
CN104985326A (zh) | 一种T型接头双侧激光-InFocus电弧复合焊接方法 | |
CN109885946B (zh) | 一种确定复合热源的能量分配的方法及焊接模拟方法 | |
CN104438480B (zh) | 一种采用工程拼焊板实现复杂零件渐进成形的加工方法 | |
CN112329309B (zh) | 环状熔核电阻点焊数值模拟方法 | |
Tian et al. | Sensitivity analysis for process parameters influencing electric arc cutting | |
Zhao et al. | Overlapping optimization of hybrid deposited and micro-rolling additive manufacturing | |
CN112317960A (zh) | 一种基于icme的激光焊接全工艺方法 | |
CN113033039B (zh) | 单丝气电立焊焊接热源模型的建立方法 | |
CN111914383A (zh) | 一种基于壳-实体耦合单元的t型结构激光焊接模拟方法 | |
CN217253493U (zh) | 一种真空钎焊分级焊接工装 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190625 |