CN107330181A - 预测激光焊接变形量的实现方法 - Google Patents
预测激光焊接变形量的实现方法 Download PDFInfo
- Publication number
- CN107330181A CN107330181A CN201710505266.5A CN201710505266A CN107330181A CN 107330181 A CN107330181 A CN 107330181A CN 201710505266 A CN201710505266 A CN 201710505266A CN 107330181 A CN107330181 A CN 107330181A
- Authority
- CN
- China
- Prior art keywords
- source
- heat source
- gauss
- molten bath
- workpiece
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
- B23K37/00—Auxiliary devices or processes, not specially adapted to a procedure covered by only one of the preceding main groups
Abstract
一种预测激光焊接变形量的实现方法,通过采集焊接后熔池区域的几何形貌并建立有限元模型以及由高斯面热源和高斯体热源组成的热源三维扫描模型,然后拟合热源参数并识别非完全耦合过程的耦合参数,再建立生产中的待连工件的有限元模型,并载入热源三维扫描模型和耦合参数,进行有限元仿真计算,最终获得待连工件的预测变形量。本发明提高了温度场计算的准确性,提高了变形预测的准确性,大大减少了焊接变形预测的时间。
Description
技术领域
本发明涉及的是一种激光焊接领域的技术,具体是一种预测激光焊接变形量的实现方法。
背景技术
激光焊接过程逐渐开始广泛的应用于汽车车身板件的连接过程中,其具有自动化程度高、 焊接速度快、焊接熔深深、接合强度高等优点,但也伴随着焊接变形、残余应力等问题。激光 焊接引起的变形往往超出车身板件的装配公差,对装配精度有很大的影响。
发明内容
本发明针对现有技术无法通过定义焊丝网格轮廓控制焊缝熔池区形貌、无法准确拟合焊 缝处熔池区形貌,焊接变形仿真效果较差的缺陷,提出一种预测激光焊接变形量的实现方法, 提高了温度场计算的准确性,提高了变形预测的准确性,大大减少了焊接变形预测的时间。
本发明是通过以下技术方案实现的:
本发明通过采集焊接后熔池区域的几何形貌并建立有限元模型以及由高斯面热源和高 斯体热源组成的热源三维扫描模型,然后拟合热源参数并识别非完全耦合过程的耦合参数,再 建立生产中的待连工件的有限元模型,并载入热源三维扫描模型和耦合参数,进行有限元仿真 计算,最终获得待连工件的预测变形量。
所述方法具体包括以下步骤:
1)对试验工件进行激光焊接,并测量焊接后的工件变形量;
2)对焊缝区域进行金相显微观察,获得熔池区域的几何形貌;
3)建立试验工件的有限元模型;
4)建立由高斯面热源和高斯体热源组成的热源三维扫描模型,拟合热源参数;
5)识别非完全耦合过程的耦合参数;
6)建立生产中的待连工件的有限元模型,并载入热源三维扫描模型和耦合参数,进行有 限元仿真计算,获得待连工件的预测变形量。
所述的步骤2)具体包括以下步骤:
2.1)将焊缝区域切割为2×2cm方块;
2.2)对方块进行镶嵌、打磨、抛光和腐蚀;
2.3)对方块进行金相显微观察获得熔池区域的几何形貌。
所述的步骤4)具体包括以下步骤:
4.1)生成由高斯面热源和高斯体热源组成的热源三维扫描模型;
4.2)设定一组热源参数,进行试验工件的仿真焊接;
4.3)对比仿真结果中的熔池区域的几何形貌和金相中的熔池区域的几何形貌的关键特征 (如熔池边界形貌、熔池面积),关键特征符合则将该组热源参数作为最终值,否则回到步骤4.2)。
所述的热源参数包括:吸收率η,功率分配系数μ,面热源半径R0,体热源上面半径ru, 体热源下面半径r1,体热源深度D。
所述的高斯面热源的能量分布为
所述的高斯体热源的能量分布为
所述的步骤5)具体包括以下步骤:
5.1)设定一组非完全耦合过程的耦合参数;
5.2)进行试验工件的仿真焊接;
5.3)仿真结果中的试验工件变形量与实际变形量之差小于阈值则该组耦合参数作为最终 值,否则回到步骤5.1)。
所述的耦合参数包括温度场计算步长ttemp和温度场个数n。
所述的焊接为激光焊接。
附图说明
图1为本发明方法流程示意图;
图2为试验工件结构图;
图3为试验工件变形量示意图;
图4为试验工件的有限元模型示意图;
图5为几何形貌对比示意图;
图中:(A)为金相结果的几何形貌图;(B)为仿真结果的几何形貌图;
图6为热源三维扫描模型结构示意图;
图7为实施例中待连工件的预测变形量有限元仿真示意图。
具体实施方式
如图1所示,本实施例包括以下步骤:
1)对试验工件进行焊接,并测量焊接后的工件变形量。
如图2所示,所述的试验工件为汽车车身中的B柱区域板件,即采用激光焊接方法来焊 接B柱内板和B柱加强版。B柱加强版和B柱内板为试验工件。对试验工件进行光学扫描,获 得试验工件的几何外形,用于建立有试验工件的限元模型。
如图3所示,所述的试验工件完成焊接后,进行三维光学扫描,获得焊接后的试验工件 的变形量。
2)对焊缝区域进行金相显微观察,获得熔池区域的几何形貌。
2.1)将焊缝区域切割为2×2cm方块。
2.2)对方块进行镶嵌、打磨、抛光和腐蚀。
2.3)对方块进行金相显微观察获得熔池区域的几何形貌。
3)建立试验工件的有限元模型。如图4所示,采用Simufact Welding软件建立的试验工 件的有限元模型,包括三维体网格生成的板件网格模型、焊缝路径设置以及夹具位置设置。
4)建立由高斯面热源和高斯体热源组成的热源三维扫描模型,拟合热源参数。
4.1)生成由高斯面热源和高斯体热源组成的热源三维扫描模型。
4.2)设定一组热源参数,进行试验工件的仿真焊接。
4.3)如图5所示,对比仿真结果中的熔池区域的几何形貌和金相中的熔池区域的几何形 貌的关键特征(熔池边界形貌、熔池面积),关键特征符合则将该组热源参数作为最终值,否则 回到步骤4.2)。热源参数最终值如表1所示。
表1热源参数数值
如图6所示,所述的热源三维扫描模型包括高斯面热源和高斯体热源。热源参数包括:吸收率η,功率分配系数μ,面热源半径R0,ru体热源上面半径,r1体热源下面半径,D体热源深度。对于输入功率为Q,吸收率η的激光热源,高斯面热源能量分布为高斯体热源能量分布为 其中:QVol=μηQ,QSurf=(1-μ)ηQ,r坐 标点的径向坐标,r0为深度d以下的体热源半径,α为激光入射的倾 斜角度。
5)识别非完全耦合过程的耦合参数。
5.1)设定非完全耦合过程的耦合参数,耦合参数包括温度场计算步长ttemp和温度场个数n。
5.2)进行试验工件的仿真焊接。
5.3)仿真结果中的试验工件变形量与实际变形量之差小于阈值则该组耦合参数作为最终 值,否则回到步骤5.1)。得到温度场计算步长ttemp和温度场个数n的最终值为0.02s和5。
6)建立生产中的待连工件的有限元模型,并载入热源三维扫描模型和耦合参数,进行有 限元仿真计算,获得待连工件的预测变形量,如图7所示。仿真结果与扫描结果对比如表2所 示,最大误差不超过实测误差的10%。
表2仿真结果与扫描结果对比
变形量 | 1 | 2 | 3 | 4 | 5 | 6 |
扫描(mm) | 1.01 | 1.15 | 1.25 | 1.23 | 0.89 | 0.21 |
仿真(mm) | 0.94 | 1.12 | 1.19 | 1.13 | 0.85 | 0.20 |
误差(%) | -6.9 | -2.6 | -4.8 | -8.1 | -4.5 | -4.8 |
与现有技术相比,本发明通过焊接试验及相应焊接件的金相试验与光学扫描,即可得到 预测程度较高的仿真模型,无需进行焊接过程的温度曲线测量与残余应力测量。仿真过程采用 非完全耦合的过程,相比相同热分析步长的完全耦合过程,仿真精度将近,仿真时长缩短为36%, 如表3所示。
表3非完全耦合方法与完全耦合方法结果对比
实例 | 非完全耦合 | 完全耦合 | 扫描结果 |
热分析步长(s) | 0.02s | 0.02s | \ |
耦合分析步长(s) | 0.1s | 0.02s | \ |
计算时间 | 0.36T | T | \ |
测点最大变形量(mm) | 1.19 | 1.21 | 1.25 |
相对误差(%) | -4.8 | -3.2 | 0 |
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式 对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围 内的各个实现方案均受本发明之约束。
Claims (8)
1.一种预测激光焊接变形量的实现方法,其特征在于,通过采集焊接后熔池区域的几何形貌并建立有限元模型以及由高斯面热源和高斯体热源组成的热源三维扫描模型,然后拟合热源参数并识别非完全耦合过程的耦合参数,再建立生产中的待连工件的有限元模型,并载入热源三维扫描模型和耦合参数,进行有限元仿真计算,最终获得待连工件的预测变形量。
2.根据权利要求1所述的方法,其特征是,所述方法具体包括以下步骤:
1)对试验工件进行焊接,并测量焊接后的工件变形量;
2)对焊缝区域进行金相显微观察,获得熔池区域的几何形貌;
3)建立试验工件的有限元模型;
4)建立由高斯面热源和高斯体热源组成的热源三维扫描模型,拟合热源参数;
5)识别非完全耦合过程的耦合参数;
6)建立生产中的待连工件的有限元模型,并载入热源三维扫描模型和耦合参数,进行有限元仿真计算,获得待连工件的预测变形量。
3.根据权利要求2所述的方法,其特征是,所述的步骤2)具体包括以下步骤:
2.1)将焊缝区域切割为2×2cm方块;
2.2)对方块进行镶嵌、打磨、抛光和腐蚀;
2.3)对方块进行金相显微观察获得熔池区域的几何形貌。
4.根据权利要求2所述的方法,其特征是,所述的步骤4)具体包括以下步骤:
4.1)生成由高斯面热源和高斯体热源组成的热源三维扫描模型;
4.2)设定一组热源参数,进行试验工件的仿真焊接;
4.3)对比仿真结果中的熔池区域的几何形貌和金相中的熔池区域的几何形貌的关键特征,关键特征符合则将该组热源参数作为最终值,否则回到步骤4.2);
所述的热源参数包括:吸收率η,功率分配系数μ,面热源半径R0,体热源上面半径ru,体热源下面半径r1,体热源深度D。
5.根据权利要求4所述的方法,其特征是,所述的关键特征包括:熔池边界形貌和熔池面积。
6.根据权利要求4所述的方法,其特征是,所述的高斯面热源的能量分布为
所述的高斯体热源的能量分布为
7.根据权利要求2所述的方法,其特征是,所述的步骤5)具体包括以下步骤:
5.1)设定一组非完全耦合过程的耦合参数;
5.2)进行试验工件的仿真焊接;
5.3)仿真结果中的试验工件变形量与实际变形量之差小于阈值则该组耦合参数作为最终值,否则回到步骤5.1)。
8.根据权利要求7所述的方法,其特征是,所述的耦合参数包括温度场计算步长ttemp和温度场个数n。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710505266.5A CN107330181A (zh) | 2017-06-28 | 2017-06-28 | 预测激光焊接变形量的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710505266.5A CN107330181A (zh) | 2017-06-28 | 2017-06-28 | 预测激光焊接变形量的实现方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107330181A true CN107330181A (zh) | 2017-11-07 |
Family
ID=60198376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710505266.5A Pending CN107330181A (zh) | 2017-06-28 | 2017-06-28 | 预测激光焊接变形量的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107330181A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108399280A (zh) * | 2018-01-30 | 2018-08-14 | 南京理工大学 | 一种预测激光选区熔化成型件变形的有限元仿真方法 |
CN108971823A (zh) * | 2018-10-10 | 2018-12-11 | 桂林航天工业学院 | 一种基于焊接工件三维重构的焊接装置及方法 |
CN110355472A (zh) * | 2019-05-10 | 2019-10-22 | 南京航空航天大学 | 一种基于约束力预测薄板拼焊变形的方法 |
CN111805081A (zh) * | 2020-05-20 | 2020-10-23 | 中车青岛四方机车车辆股份有限公司 | 抑制铝合金焊接热裂纹的方法及焊接装置 |
CN111914383A (zh) * | 2019-05-10 | 2020-11-10 | 南京航空航天大学 | 一种基于壳-实体耦合单元的t型结构激光焊接模拟方法 |
CN112016228A (zh) * | 2020-08-14 | 2020-12-01 | 河海大学常州校区 | 一种基于水冷压缩系数的水下焊接热源模型的建模方法 |
CN113191037A (zh) * | 2021-03-23 | 2021-07-30 | 深圳市信维通信股份有限公司 | 一种基于形位公差的插损预测方法及终端 |
CN113523656A (zh) * | 2021-07-05 | 2021-10-22 | 江苏徐工工程机械研究院有限公司 | 一种结构件焊接反变形设计方法 |
CN114245764A (zh) * | 2019-08-09 | 2022-03-25 | 株式会社神户制钢所 | 层叠造型物的层叠计划方法、层叠造型物的制造方法以及制造装置 |
CN117252000A (zh) * | 2023-09-15 | 2023-12-19 | 南昌航空大学 | 电子束倾斜焊接的热源模型构建方法及系统、仿真方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102152016A (zh) * | 2010-02-03 | 2011-08-17 | 株式会社日立制作所 | 焊接变形分析方法 |
CN103049623A (zh) * | 2013-01-18 | 2013-04-17 | 哈尔滨工业大学 | 一种激光焊接热源模型的建立方法 |
CN104298817A (zh) * | 2014-09-26 | 2015-01-21 | 中国船舶重工集团公司第七二五研究所 | 一种用于模拟高能束焊接的双椭圆指数衰减体热源模型 |
-
2017
- 2017-06-28 CN CN201710505266.5A patent/CN107330181A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102152016A (zh) * | 2010-02-03 | 2011-08-17 | 株式会社日立制作所 | 焊接变形分析方法 |
CN103049623A (zh) * | 2013-01-18 | 2013-04-17 | 哈尔滨工业大学 | 一种激光焊接热源模型的建立方法 |
CN104298817A (zh) * | 2014-09-26 | 2015-01-21 | 中国船舶重工集团公司第七二五研究所 | 一种用于模拟高能束焊接的双椭圆指数衰减体热源模型 |
Non-Patent Citations (3)
Title |
---|
卢宇峰等: "激光焊接圆锥体热源模型及参数研究", 《第十六次全国焊接学术会议》 * |
宋燕利等: "基于BP和GA的激光焊接热源模型参数优化", 《塑性工程学报》 * |
陈细亮: "高温合金环形件电子束焊接温度场和变形的数值模拟", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108399280B (zh) * | 2018-01-30 | 2020-05-05 | 南京理工大学 | 一种预测激光选区熔化成型件变形的有限元仿真方法 |
CN108399280A (zh) * | 2018-01-30 | 2018-08-14 | 南京理工大学 | 一种预测激光选区熔化成型件变形的有限元仿真方法 |
CN108971823A (zh) * | 2018-10-10 | 2018-12-11 | 桂林航天工业学院 | 一种基于焊接工件三维重构的焊接装置及方法 |
CN111914383A (zh) * | 2019-05-10 | 2020-11-10 | 南京航空航天大学 | 一种基于壳-实体耦合单元的t型结构激光焊接模拟方法 |
CN110355472A (zh) * | 2019-05-10 | 2019-10-22 | 南京航空航天大学 | 一种基于约束力预测薄板拼焊变形的方法 |
CN114245764A (zh) * | 2019-08-09 | 2022-03-25 | 株式会社神户制钢所 | 层叠造型物的层叠计划方法、层叠造型物的制造方法以及制造装置 |
CN114245764B (zh) * | 2019-08-09 | 2023-08-29 | 株式会社神户制钢所 | 层叠造型物的层叠计划方法、层叠造型物的制造方法以及制造装置 |
CN111805081A (zh) * | 2020-05-20 | 2020-10-23 | 中车青岛四方机车车辆股份有限公司 | 抑制铝合金焊接热裂纹的方法及焊接装置 |
CN112016228A (zh) * | 2020-08-14 | 2020-12-01 | 河海大学常州校区 | 一种基于水冷压缩系数的水下焊接热源模型的建模方法 |
CN113191037A (zh) * | 2021-03-23 | 2021-07-30 | 深圳市信维通信股份有限公司 | 一种基于形位公差的插损预测方法及终端 |
CN113523656A (zh) * | 2021-07-05 | 2021-10-22 | 江苏徐工工程机械研究院有限公司 | 一种结构件焊接反变形设计方法 |
CN113523656B (zh) * | 2021-07-05 | 2023-01-17 | 江苏徐工工程机械研究院有限公司 | 一种结构件焊接反变形设计方法 |
CN117252000A (zh) * | 2023-09-15 | 2023-12-19 | 南昌航空大学 | 电子束倾斜焊接的热源模型构建方法及系统、仿真方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107330181A (zh) | 预测激光焊接变形量的实现方法 | |
CN204747769U (zh) | 一种激光在线测量加工检测装置 | |
CN104808268B (zh) | 多焦点透镜及其加工设备、加工方法 | |
CN110328848A (zh) | 一种激光修复方法及装置 | |
CN104674210A (zh) | 一种工件激光自动化修复方法 | |
CN106312323B (zh) | 一种变形叶片的激光喷丸校形方法及装置 | |
WO2018032771A1 (zh) | 一种用于航空发动机支架的激光喷丸校形质量控制方法 | |
CN103323526A (zh) | 一种基于超声相控阵和支持向量机的焊缝缺陷检测识别方法 | |
CN107301261A (zh) | 基于comsol温度模型模拟计算激光加工与焊接过程中温度场的方法 | |
CN103488832A (zh) | 一种复杂曲面零件破损区域的几何修复方法 | |
CN101219528A (zh) | 燃气涡轮机轮叶和叶片的水射流剥离和再定轮廓 | |
CN106312299B (zh) | 航空发动机支架激光喷丸校形形状精度在线控制的方法与装置 | |
CN106270005B (zh) | 一种叶片激光喷丸校形的方法与装置 | |
WO2021212894A1 (zh) | 未熔合缺陷的无损检测方法、检测标准件及其制造方法 | |
CN112380677B (zh) | 基于激光清洗剥离碳纤维材料脱模剂中温度场的模拟方法 | |
CN109813758B (zh) | 基于红外检测图像信号的层间脱粘缺陷复杂边界轮廓定量重构方法 | |
CN112404883B (zh) | 一种实时精准液态微型锻造增材再制造方法及装置 | |
CN104807818B (zh) | 三维静态及动态显微检测系统及方法 | |
CN111767666B (zh) | 一种汽车零部件激光焊连接的cae仿真模拟方法 | |
CN110320236B (zh) | 大型风力机叶片内部缺陷深度的红外测量方法 | |
CN110548876A (zh) | 一种铺粉式再制造装置及方法 | |
CN109557675A (zh) | 基于非球面镜像差效应的长焦深激光光束均匀化光学系统 | |
CN106910694A (zh) | 一种自动检测集成电路器件翘曲度的方法及装置 | |
CN112881533A (zh) | 一种超声检测试块及其制备方法 | |
CN114406476A (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 |
Application publication date: 20171107 |
|
RJ01 | Rejection of invention patent application after publication |