CN110889213B - 一种基于切片分割加载体的热源焊接模拟方法 - Google Patents

一种基于切片分割加载体的热源焊接模拟方法 Download PDF

Info

Publication number
CN110889213B
CN110889213B CN201911128606.2A CN201911128606A CN110889213B CN 110889213 B CN110889213 B CN 110889213B CN 201911128606 A CN201911128606 A CN 201911128606A CN 110889213 B CN110889213 B CN 110889213B
Authority
CN
China
Prior art keywords
welding
model
loaded
generation rate
molten pool
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.)
Active
Application number
CN201911128606.2A
Other languages
English (en)
Other versions
CN110889213A (zh
Inventor
迟露鑫
麻永林
邢淑清
陈重毅
梁仕发
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Technology
Original Assignee
Chongqing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Technology filed Critical Chongqing University of Technology
Priority to CN201911128606.2A priority Critical patent/CN110889213B/zh
Publication of CN110889213A publication Critical patent/CN110889213A/zh
Application granted granted Critical
Publication of CN110889213B publication Critical patent/CN110889213B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Abstract

本发明公开了一种基于切片分割加载体的热源焊接模拟方法，包括：获取目标焊接件模型；基于目标焊接件模型建立焊缝模型；将焊缝模型沿熔池深度方向及熔池长度方向划分为多个待加载体；对每个待加载体分别加载体生热率热源模型，下层待加载体的体生热率大于或等于上层待加载体的体生热率；提取体生热率热源模型加载后的目标焊接件模型模拟数据。与现有技术相比，本发明在总热输入不变的前提下，施加不同比例的体生热率，即下方较高，上方较低，通过焊材上下表明温度分布的调节，实现焊缝、热影响区较为精确的模拟计算，提高了对SAPH440钢材料的焊接模拟结果的准确性。

Description

一种基于切片分割加载体的热源焊接模拟方法
技术领域
本发明涉及实验仿真领域，具体涉及一种基于切片分割加载体的热源焊接模拟方法。
背景技术
研究学者通过对焊接过程数值模拟的研究，希望在创建焊接温度场模型时，能够客观地反映焊接过程，尽可能地考虑进所有影响焊接过程的因素，以便能够准确预测焊接温度场分布。
然后现有技术中采用普通体生热加载的方法，并不能得到准确的模拟结果，例如，在针对SAPH440钢的焊接模拟中，普通体生热加载在模拟熔池温度分布存在一定的局限性，受到焊接方法、焊接材料、焊接参数、焊件厚度等影响，使得热源分布不能与实际热源吻合，模拟结果存在较大误差，而焊接温度场分布直接影响接头组织性能，更决定接头残余应力大小。因此，如何采用基于实验数据，修正现有的热源模型，使得热源模型能够针对SAPH440钢脉冲弧焊准确模拟不同工艺的焊接温度场分布，成为了本领域技术人员急需解决的问题。
发明内容
针对现有技术存在的上述不足，本发明要解决的技术问题是：如何提高焊接模拟结果的准确性。
为解决上述技术问题，本发明采用了如下的技术方案：
一种基于切片分割加载体的热源焊接模拟方法，包括：
获取目标焊接件模型；
基于目标焊接件模型建立焊缝模型；
将焊缝模型沿熔池深度方向及熔池长度方向划分为多个待加载体；
对每个待加载体分别加载体生热率热源模型，下层待加载体的体生热率大于或等于上层待加载体的体生热率；
提取体生热率热源模型加载后的目标焊接件模型模拟数据。
优选地，体生热率采用以下公式计算：
Figure BDA0002277634620000011
其中，HGEN为体生热率；η为焊接热效率；U为电弧电压；I为焊接电流；v为焊速；S为加载的体单元面积；t为加载时间。
优选地，所述焊缝模型中，熔池深度方向与熔池长度方向的尺寸比为10：4。
优选地，将焊缝模型沿熔池深度方向按照截面面积3：3：4的比例，沿熔池长度方向按照截面面积2：1：2的比例切割分层得到待加载体。
优选地，沿熔池深度方向的待加载体的体生热率比为3：3：4，沿熔池长度方向的待加载体的体生热率比为1：2：1。
综上所述，本发明公开了一种基于切片分割加载体的热源焊接模拟方法，包括：获取目标焊接件模型；基于目标焊接件模型建立焊缝模型；将焊缝模型沿熔池深度方向及熔池长度方向划分为多个待加载体；对每个待加载体分别加载体生热率热源模型，下层待加载体的体生热率大于或等于上层待加载体的体生热率；提取体生热率热源模型加载后的目标焊接件模型模拟数据。与现有技术相比，本发明在总热输入不变的前提下，施加不同比例的体生热率，即下方较高，上方较低，通过焊材上下表明温度分布的调节，实现焊缝、热影响区较为精确的模拟计算，提高了对SAPH440钢材料的焊接模拟结果的准确性。
附图说明
图1为本发明公开的一种基于切片分割加载体的热源焊接模拟方法的一种具体实施方式的流程图；
图2为本发明实施例中示例目标焊接件模型的网格划分示意图；
图3为本发明实施例中示例焊缝模型的分割示意图；
图4为本发明中按照不同方式沿熔池深度划分待加载体之后的模拟数据对比图；
图5为本发明中按照不同方式沿熔池长度划分待加载体之后的模拟数据对比图；
图6为现实中的焊接接头示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步的详细描述说明。
如图1所示，本发明公开了一种基于切片分割加载体的热源焊接模拟方法，包括：
获取目标焊接件模型；
基于目标焊接件模型建立焊缝模型；
将焊缝模型沿熔池深度方向及熔池长度方向划分为多个待加载体；
对每个待加载体分别加载体生热率热源模型，下层待加载体的体生热率大于或等于上层待加载体的体生热率；
提取体生热率热源模型加载后的目标焊接件模型模拟数据。
现有技术中，在进行焊接模拟时，对于焊缝区域，单位体积的热输入相同，有限元模型下层体积相对于上层体积变小，然而上下层有限元模型在计算过程中却有同样的接触面积向两侧母材导热，使得焊缝下方温度偏低，造成焊缝熔深温度梯度过大，熔深较短，与实际不符。特别是网格的疏密代表着计算的精确度，因为焊接弧柱直径实际3mm～4mm，沿着熔深方向加载一层单元，造成热源在熔深的深度不够，电弧造成熔池后托性，沿着焊接方向上熔池长度往往与实际热源宽度不符。因此，本发明焊缝模型沿熔池深度方向及熔池长度方向划分为多个待加载体，并且根据深度的不同施加不同的体生热率，从而使得在总热输入不变的前提下，焊缝模型中下方的待加载体的体生热率较高，上方的待加载体的体生热率较低，通过焊材上下表明温度分布的调节，实现焊缝、热影响区较为精确的模拟计算，提高了对SAPH440G钢材料的焊接模拟结果的准确性。
具体实施时，体生热率采用以下公式计算：
Figure BDA0002277634620000031
其中，HGEN为体生热率；η为焊接热效率；U为电弧电压；I为焊接电流；v为焊速；S为加载的体单元面积；t为加载时间。
在本发明中，η的取值可为0.85，t的取值可为0.1s。初始条件：Tt＝0＝25℃，边界条件：对流换热系数α＝10.8W/(m2·℃)。
具体实施时，所述焊缝模型中，熔池深度方向与熔池长度方向的尺寸比为10：4。
如图6所示，通过实验发现，焊接接头尺寸通常为熔宽10mm,熔深4mm,熔池长度19mm，因此取熔池深度方向与熔池长度方向的尺寸比为10：4。
如图3所示，具体实施时，将焊缝模型沿熔池深度方向按照截面面积3：3：4的比例、沿熔池长度方向按照截面面积2：1：2的比例切割分层得到待加载体。
具体实施时，沿熔池深度方向的待加载体的体生热率比为3：3：4，沿熔池长度方向的待加载体的体生热率比为1：2：1。
在本发明中，以大小为100mm×240mm×4mm的目标焊接件为例，焊接时为避免热源尺寸的失真，沿着熔深方向焊缝及粗晶区处的网格为0.1mm，远离焊缝区域的网格粗大，焊件有限元模型如图2所示。在模拟实际焊接过程，利用APDL语言中的*DO循环代码模拟单元的激活及热源的移动。
对有限元模型每秒填缝金属熔深方向加载热源以单层(方案一)、分割两层(方案二)、分割三层(方案三)进行对比。
焊缝加载热源按传统均匀分布形式(方案一)，熔池最高温度可达2219℃，但是，在焊缝底部温度未达到1500℃，熔池深度方向温度梯度过大，模型没有熔透。
焊缝热源加载沿着熔深方向均分两层(方案二)，总热输入按1:2的比例分配能量，温度沿熔深温度梯度较缓，焊缝底部的温度明显高于1500℃，但熔化金属范围较窄。
焊缝热源加载沿着熔深方向为三层(方案三，即本发明的方案)，总热输入按3:3:4的比例施加能量在上、中、下三层，在总热输入不变的前提下，总体温度梯度较小，在熔深方向上温度梯度均匀。焊缝底部的熔化金属较宽范围内温度在1500℃以上，基本熔透。改变熔池截面的热输入可以有效地控制熔深及熔池的温度分布。如附图4所示。
在沿着熔池深度方向将热输入分3层的基础上，将熔池长度方向热源载荷将单位时间填充长度切割为不同片数，采用以下三个方案进行模拟：不切割单片加载(方案一)，切割2片加载，且比例为1：1.5(方案二)，切割3片体加载，以比例为1：2：1加载热源(方案三，即本发明的方案)。
沿着焊接方向熔池长度方向加载单片单元(方案一)，相当于加载长度为每秒焊接速度，加载体积较大，温度明显低于实际焊接温度，熔长由于模拟时沿焊接方向加载单元过少，小于焊接时中熔池熔长，熔池宽度较小，熔池前方等温线梯度较大，误差较大。
沿着焊接方向熔池长度方向加载2片单元(方案二)，按着1：1的形式加载，整体热输入较单片加载单元有所提高，熔长随着加载单元的增加变长，热源模型尺寸超过熔池宽度尺寸，然而，熔池长度后托性明显不足。
沿着焊接方向熔池长度方向加载3片单元(方案三)，接近实际电弧弧柱直径，焊缝的最高温度区间明显提升，热输入所加载的体积与实际热加载区较为接近。熔长为6.3mm，熔宽为4.2mm。熔池宽度没有发生改变，熔长大幅度增大，更加接近实际。如图5所示。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管通过参照本发明的优选实施例已经对本发明进行了描述，但本领域的普通技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离所附权利要求书所限定的本发明的精神和范围。

Claims (3)

1.一种基于切片分割加载体的热源焊接模拟方法，其特征在于，包括：
获取目标焊接件模型；
基于目标焊接件模型建立焊缝模型；
将焊缝模型沿熔池深度方向及熔池长度方向划分为多个待加载体；
对每个待加载体分别加载体生热率热源模型，下层待加载体的体生热率大于或等于上层待加载体的体生热率；
提取体生热率热源模型加载后的目标焊接件模型模拟数据；
其中，将焊缝模型沿熔池深度方向按照截面面积3：3：4的比例、沿熔池长度方向按照截面面积2：1：2的比例切割分层得到待加载体；沿熔池深度方向的待加载体的体生热率比为3：3：4，沿熔池长度方向的待加载体的体生热率比为1：2：1。
2.如权利要求1所述的基于切片分割加载体的热源焊接模拟方法，其特征在于，体生热率采用以下公式计算：
Figure FDA0004149468330000011
其中，HGEN为体生热率；η为焊接热效率；U为电弧电压；I为焊接电流；v为焊速；S为加载的体单元面积；t为加载时间。
3.如权利要求1所述的基于切片分割加载体的热源焊接模拟方法，其特征在于，所述焊缝模型中，熔池深度方向与熔池长度方向的尺寸比为10：4。
CN201911128606.2A 2019-11-18 2019-11-18 一种基于切片分割加载体的热源焊接模拟方法 Active CN110889213B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911128606.2A CN110889213B (zh) 2019-11-18 2019-11-18 一种基于切片分割加载体的热源焊接模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911128606.2A CN110889213B (zh) 2019-11-18 2019-11-18 一种基于切片分割加载体的热源焊接模拟方法

Publications (2)

Publication Number Publication Date
CN110889213A CN110889213A (zh) 2020-03-17
CN110889213B true CN110889213B (zh) 2023-05-02

Family

ID=69747852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911128606.2A Active CN110889213B (zh) 2019-11-18 2019-11-18 一种基于切片分割加载体的热源焊接模拟方法

Country Status (1)

Country Link
CN (1) CN110889213B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109885946A (zh) * 2019-02-26 2019-06-14 大族激光科技产业集团股份有限公司 一种确定复合热源的能量分配的方法及焊接模拟方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3899008B2 (ja) * 2002-10-08 2007-03-28 新日本製鐵株式会社 突合せ溶接継手の疲労強度向上方法
US20090266870A1 (en) * 2008-04-23 2009-10-29 The Boeing Company Joined composite structures with a graded coefficient of thermal expansion for extreme environment applications
US8498728B2 (en) * 2009-10-19 2013-07-30 Edison Welding Institute, Inc. Remote high-performance modeling system for material joining and material forming
JP5504855B2 (ja) * 2009-12-04 2014-05-28 株式会社Ｉｈｉ レーザ溶接方法
JP2011198684A (ja) * 2010-03-23 2011-10-06 Sumitomo Wiring Syst Ltd コネクタ及びその製造方法
WO2013014481A1 (fr) * 2011-07-26 2013-01-31 Arcelormittal Investigación Y Desarrollo Sl Pièce d'acier soudée préalablement mise en forme à chaud à très haute résistance mécanique et procédé de fabrication
CN103801795B (zh) * 2013-12-12 2017-05-03 内蒙古科技大学 一种超细晶钢ss400中厚板焊接方法
CN106132622B (zh) * 2013-12-27 2018-08-31 杰富意钢铁株式会社 电阻点焊方法
CN104007182B (zh) * 2014-06-12 2016-04-27 重庆理工大学 一种利用声发射信号频谱定量检测电阻点焊熔核形核质量的方法
EP3009220B1 (de) * 2014-10-17 2018-12-19 Rehm Thermal Systems GmbH Rolle-zu-Rolle-Fertigungsanlage und -verfahren für verkettete kontinuierliche und diskontinuierliche Verarbeitungsprozesse
CN105234599B (zh) * 2015-10-20 2018-06-12 沈阳富创精密设备有限公司 焊接温度场控制系统及方法
CN108022236B (zh) * 2017-11-28 2019-07-09 华中科技大学 基于焊缝图像提取拟合的焊接应力变形预测方法及系统
CN109299554A (zh) * 2018-09-30 2019-02-01 西南石油大学 一种激光-电弧复合焊接热源模型的建立方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109885946A (zh) * 2019-02-26 2019-06-14 大族激光科技产业集团股份有限公司 一种确定复合热源的能量分配的方法及焊接模拟方法

Also Published As

Publication number Publication date
CN110889213A (zh) 2020-03-17

Similar Documents

Publication Publication Date Title
CN108399307A (zh) 一种激光3d打印有限元模拟方法
Montevecchi et al. Finite element modelling of wire-arc-additive-manufacturing process
CN109777942B (zh) 超大型压力容器局部焊后热处理焊接应力与变形调控方法
JP5712306B2 (ja) 三次元体の製造方法
CN102693336B (zh) 一种预测大型管道焊接热循环参数的方法
KR101996933B1 (ko) 적층 생성된 구조물에서 응력들 및 형상 편차들을 결정하기 위한 방법, 컴퓨터-판독가능 데이터 캐리어, 컴퓨터 프로그램, 및 시뮬레이터
CN103246774B (zh) 数值模拟p92钢管焊接热影响区宽度的方法
CN112149330B (zh) 一种风电塔筒油封平台焊接残余应力预测、焊接工艺优化方法
CN112380752B (zh) 利用预测金属薄板焊接热处理数值改进其焊接工艺的方法
CN105138772B (zh) 变截面构件电子束焊的有限元模拟方法
Jithin et al. FE modeling for single spark in EDM considering plasma flushing efficiency
CN110705156A (zh) 一种多层多道焊仿真方法、装置、设备及存储介质
CN105598448A (zh) 一种金属材料激光3d打印原位预热温度的控制方法
CN110889213B (zh) 一种基于切片分割加载体的热源焊接模拟方法
Cheepu et al. Numerical simulations of the effect of heat input on microstructural growth for MIG-based wire arc additive manufacturing of Inconel 718
SG173281A1 (en) A method of restoring a metallic component
CN109735843B (zh) 一种增加激光熔覆高硬度合金层厚度的工艺方法及其激光熔覆修复的产品
CN105799175B (zh) 一种维持工作缸整体温度均衡的装置及其控制方法
JP2019084723A (ja) 三次元造形方法及び三次元造形装置、及びこれらによって造形された立体形状物
Chua et al. Investigation of penetration depth and efficiency of applied heat flux in a directed energy deposition process with feeding of Ti-6Al-4V wires
CN112883623A (zh) 一种Fe-36%Ni合金PMIG摆动焊热-机械耦合有限元建模方法
CN107138727A (zh) 一种具有点阵冷却结构的扇形封严块制备方法
CN110866359B (zh) 一种基于修正双椭球热源模型的焊接模拟方法
CN109885946B (zh) 一种确定复合热源的能量分配的方法及焊接模拟方法
Zubairuddin et al. Finite element simulation of weld bead geometry and temperature distribution during GTA welding of modified 9Cr-1Mo steel and experimental validation

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
GR01 Patent grant
GR01 Patent grant