CN103761374A - 一种激光深熔焊热源模型的建模方法 - Google Patents

一种激光深熔焊热源模型的建模方法 Download PDF

Info

Publication number
CN103761374A
CN103761374A CN201410012392.3A CN201410012392A CN103761374A CN 103761374 A CN103761374 A CN 103761374A CN 201410012392 A CN201410012392 A CN 201410012392A CN 103761374 A CN103761374 A CN 103761374A
Authority
CN
China
Prior art keywords
heat source
laser
model
weld joint
source 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.)
Granted
Application number
CN201410012392.3A
Other languages
English (en)
Other versions
CN103761374B (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.)
Shanghai University of Engineering Science
Original Assignee
Shanghai University of Engineering Science
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 Shanghai University of Engineering Science filed Critical Shanghai University of Engineering Science
Priority to CN201410012392.3A priority Critical patent/CN103761374B/zh
Publication of CN103761374A publication Critical patent/CN103761374A/zh
Application granted granted Critical
Publication of CN103761374B publication Critical patent/CN103761374B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Laser Beam Processing (AREA)

Abstract

本发明公开了一种激光深熔焊热源模型的建模方法,建立T型接头实体模型,建模时综合考虑焊缝区材料属性,及网格顺序划分问题将焊缝区建为3个体,盖板焊缝区,空气膜焊缝区,芯板焊缝区;基于ANSYS,顺序划分网格,优先划分焊缝区;建立三维激光深熔焊热源模型,基于APDL,求解单元热值,热源模型控制方程为:
Figure DDA0000455613490000011
其中,x,y,z为笛卡尔坐标,QM为热源峰值,r1激光找在材料表面的有效半径r2,H1为热源形状参数,H0热源上部分高度,H为焊缝实际熔深,按次热源模型求解热值,存储到数组中,基于ANSYS软件,将求解热值,加载到有限元模型的相应单元,设置边界条件,计算求解,完成建模过程。

Description

一种激光深熔焊热源模型的建模方法
技术领域
本发明属于激光焊接技术领域,特别涉及激光焊接中深熔焊的热源模型的建模方法。
背景技术
激光焊接具有能量集中、深宽比大、热影响区小、变形小、适应能力强、焊接质量高等优点。因此在工业领域得到广泛应用,特别是在汽车、船舶、航空航天的轻量化过程中有着不可替代的优势。激光焊接可以获得较大的深宽比,从而可以实现T型接头的深熔焊。激光焊接涉及激光与材料的相互作用,小孔效应,等离子体效应,材料热物理属性,焊接熔池的传热和传质,等多变量、多耦合现象,***的实验研究费时、费力。现实中,由于激光的局部加热,在熔池及其附近激光与材料相互作用,数据不可能测量。所以,模拟是激光焊接技术中获取必要信息和预测焊接温度场分布的唯一有效方法。
建立合理的热源模型是焊接数值模拟结果准确、可靠的前提。对于激光深熔焊来说,激光热量集中,其功率密度高于106W/cm2,可以在极短时间内将金属汽化,形成小孔效应。而小孔类似光学的黑体,射入的激光经多次反射后全部被吸收,此后小孔就像焊接热源一样实现焊接过程。
发明内容
本发明旨在模拟激光深熔焊的小孔效应和熔池形状。
本发明的技术方案是,一种激光深熔焊热源模型的建模方法,包括以下步骤:
第一步,建立T型接头实体模型,建模时综合考虑焊缝区材料属性,及网格顺序划分问题将焊缝区建为3个体,盖板焊缝区,空气膜焊缝区,芯板焊缝区;
第二步,基于ANSYS,顺序划分网格,优先划分焊缝区;
第三步,建立三维激光深熔焊热源模型,基于APDL,求解单元热值,热源模型控制方程为:
Q Laser ( x , y , z ) = Q M exp [ - 9 ( x 2 + z 2 ) r 1 2 log ( H 1 y ) ] y ≥ H 0 q Laser ( x , y , z ) = Q M exp [ - 3 ( x 2 + z 2 ) r 2 2 ] H 0 > y ≥ H Q M = 9 ηQ π r 1 2 H 1 ( 1 - e - 3 ) H 0 = H 1 exp [ - 3 x 2 + z 2 r 2 2 ] - - - ( 1 )
其中,x,y,z为笛卡尔坐标,QM为热源峰值,r1激光找在材料表面的有效半径r2,H1为热源形状参数,H0热源上部分高度,H为焊缝实际熔深,按次热源模型求解热值,存储到数组中,
第四步,基于ANSYS软件,将求解热值,加载到有限元模型的相应单元,设置边界条件,计算求解,完成建模过程,
激光焊接形成的温度场控制方程为:
ρ ∂ ( c p T ) ∂ t = ∂ ∂ x ( κ ∂ T ∂ x ) + ∂ ∂ y ( κ ∂ T ∂ y ) + ∂ ∂ z ( κ ∂ T ∂ z ) + q Laser - - - ( 2 )
其中ρ是质量密度,Cp比热容,T是温度,t时间,κ是热传导系数,x,y,z是笛卡尔坐标,qLaser激光热源输入的能量,
本发明建立了全新的激光深熔焊热源模型,可以很好地模拟激光深熔焊温度场。基于APDL,建立实体模型,顺序划分网格,得到有限元模型;编写热源程序代码,计算热值存储于数组;设置边界条件,加载热值,计算求解。得到的T型接头激光深熔焊温度场与实验得到的结果吻合良好。与以前国内外所用较多热源进行了对比,表明本发明热源更适合激光深熔焊。且熔深越大,效果越明显。
附图说明
图1是本发明的新型激光焊接体热源模型
具体实施方式
本发明的具体实施建模步骤:
第一步,建立T型接头实体模型,建模时综合考虑焊缝区材料属性,及网格顺序划分问题将焊缝区建为3个体,盖板焊缝区,空气膜焊缝区,芯板焊缝区。
第二步,基于ANSYS,顺序划分网格。优先划分焊缝区。
第三步,建立三维激光深熔焊热源模型,基于APDL,求解单元热值。
热源模型控制方程为:
Q Laser ( x , y , z ) = Q M exp [ - 9 ( x 2 + z 2 ) r 1 2 log ( H 1 y ) ] y ≥ H 0 q Laser ( x , y , z ) = Q M exp [ - 3 ( x 2 + z 2 ) r 2 2 ] H 0 > y ≥ H Q M = 9 ηQ π r 1 2 H 1 ( 1 - e - 3 ) H 0 = H 1 exp [ - 3 x 2 + z 2 r 2 2 ] - - - ( 1 )
其中,x,y,z为笛卡尔坐标,QM为热源峰值,r1激光找在材料表面的有效半径r2,H1为热源形状参数,H0热源上部分高度,H为焊缝实际熔深,按次热源模型求解热值,存储到数组中。
第四步,基于ANSYS软件,将求解热值,加载到有限元模型的相应单元,设置边界条件,计算求解,完成建模过程。
激光焊接形成的温度场控制方程为:
ρ ∂ ( c p T ) ∂ t = ∂ ∂ x ( κ ∂ T ∂ x ) + ∂ ∂ y ( κ ∂ T ∂ y ) + ∂ ∂ z ( κ ∂ T ∂ z ) + q Laser - - - ( 2 )
其中ρ是质量密度,Cp比热容,T是温度,t时间,κ是热传导系数,x,y,z是笛卡尔坐标,qLaser激光热源输入的能量。如图1所示,是本发明的新型激光焊接体热源模型。
本发明的实验验证过程,选用的焊接工艺参数为激光功率Q=4500W,焊接速度V=25mm/s,激光的吸收效率为η=0.9,热源模型在材料表面的有效光斑半径r0=0.7mm,热源模型高度H=5.3mm,保护气体选用純氩,20L/min。
使用材料为304L不锈钢,其熔点为1450℃,沸点为2800℃。
使用IPG  YLS-5000型激光器,KUKA  KR60HA六轴机器人,实施焊接。
基于ANSYS二次开发技术,实现参数化建模,选用典型的solid70热实体单元,采用非均匀划分技术划分网格;应用数组存储和矢量计算,实现参数化加载求解。
最后,应用国内外常用热源模型,使用同一台计算机,采用同样的算法,模拟T型接头激光深熔焊。
本发明与国内外常用热源相比可获得更好的熔宽、熔深,熔池形状与实验焊缝更吻合。

Claims (1)

1.一种激光深熔焊热源模型的建模方法,其特征在于,包括以下步骤: 
第一步,建立T型接头实体模型,建模时综合考虑焊缝区材料属性,及网格顺序划分问题将焊缝区建为3个体,盖板焊缝区,空气膜焊缝区,芯板焊缝区; 
第二步,基于ANSYS,顺序划分网格,优先划分焊缝区; 
第三步,建立三维激光深熔焊热源模型,基于APDL,求解单元热值,热源模型控制方程为: 
其中,x,y,z为笛卡尔坐标,QM为热源峰值,r1激光找在材料表面的有效半径r2,H1为热源形状参数,H0热源上部分高度,H为焊缝实际熔深,按次热源模型求解热值,存储到数组中, 
第四步,基于ANSYS软件,将求解热值,加载到有限元模型的相应单元,设置边界条件,计算求解,完成建模过程, 
激光焊接形成的温度场控制方程为: 
Figure FDA0000455613460000012
其中ρ是质量密度,Cp比热容,T是温度,t时间,κ是热传导系数,x,y,z是笛卡尔坐标,qLaser激光热源输入的能量。 
CN201410012392.3A 2014-01-10 2014-01-10 一种激光深熔焊热源模型的建模方法 Expired - Fee Related CN103761374B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410012392.3A CN103761374B (zh) 2014-01-10 2014-01-10 一种激光深熔焊热源模型的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410012392.3A CN103761374B (zh) 2014-01-10 2014-01-10 一种激光深熔焊热源模型的建模方法

Publications (2)

Publication Number Publication Date
CN103761374A true CN103761374A (zh) 2014-04-30
CN103761374B CN103761374B (zh) 2018-07-24

Family

ID=50528610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410012392.3A Expired - Fee Related CN103761374B (zh) 2014-01-10 2014-01-10 一种激光深熔焊热源模型的建模方法

Country Status (1)

Country Link
CN (1) CN103761374B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105550434A (zh) * 2015-12-10 2016-05-04 南车株洲电力机车有限公司 一种机车车体轻量化优化方法
CN105718690A (zh) * 2016-01-26 2016-06-29 南京航空航天大学 基于时间和空间主动跟踪的激光3d打印熔池凝固行为数值模拟方法
CN106670666A (zh) * 2016-12-14 2017-05-17 南京航空航天大学 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法
CN107649804A (zh) * 2017-10-17 2018-02-02 华中科技大学鄂州工业技术研究院 一种增材制造熔深在线检测和控制***
CN109522610A (zh) * 2018-10-25 2019-03-26 南京航空航天大学 一种焊缝区混合晶区有限元建模方法
CN109933007A (zh) * 2019-01-30 2019-06-25 中国科学院金属研究所 一种激光-mig复合焊接热源模型的建立方法
CN112496614A (zh) * 2020-12-17 2021-03-16 上海工程技术大学 一种深熔焊的有限元仿真热源模型获取方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
杜汉斌等: "激光穿透焊温度场及流动场的数值模拟", 《焊接学报》 *
毕艳霞: "T型接头焊接温度场与应力场的数值模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑(月刊)》 *
熊智军 等: "激光深熔焊接热源模型的研究进展", 《焊接》 *
胥国祥 等: "铝合金T型接头激光+GMAW复合热源焊温度场的有限元分析", 《金属学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105550434A (zh) * 2015-12-10 2016-05-04 南车株洲电力机车有限公司 一种机车车体轻量化优化方法
CN105550434B (zh) * 2015-12-10 2018-08-03 南车株洲电力机车有限公司 一种机车车体轻量化优化方法
CN105718690A (zh) * 2016-01-26 2016-06-29 南京航空航天大学 基于时间和空间主动跟踪的激光3d打印熔池凝固行为数值模拟方法
CN106670666A (zh) * 2016-12-14 2017-05-17 南京航空航天大学 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法
CN107649804A (zh) * 2017-10-17 2018-02-02 华中科技大学鄂州工业技术研究院 一种增材制造熔深在线检测和控制***
CN107649804B (zh) * 2017-10-17 2023-06-23 华中科技大学鄂州工业技术研究院 一种增材制造熔深在线检测和控制***
CN109522610A (zh) * 2018-10-25 2019-03-26 南京航空航天大学 一种焊缝区混合晶区有限元建模方法
CN109933007A (zh) * 2019-01-30 2019-06-25 中国科学院金属研究所 一种激光-mig复合焊接热源模型的建立方法
CN112496614A (zh) * 2020-12-17 2021-03-16 上海工程技术大学 一种深熔焊的有限元仿真热源模型获取方法
CN112496614B (zh) * 2020-12-17 2022-09-13 上海工程技术大学 一种深熔焊的有限元仿真热源模型获取方法

Also Published As

Publication number Publication date
CN103761374B (zh) 2018-07-24

Similar Documents

Publication Publication Date Title
CN103761374A (zh) 一种激光深熔焊热源模型的建模方法
Ai et al. The prediction of the whole weld in fiber laser keyhole welding based on numerical simulation
CN109513924B (zh) 一种激光选区熔化过程中表面粗糙度控制方法
Tan et al. Numerical modeling of transport phenomena and dendritic growth in laser spot conduction welding of 304 stainless steel
CN112149330B (zh) 一种风电塔筒油封平台焊接残余应力预测、焊接工艺优化方法
Li et al. An evolutionary keyhole-mode heat transfer model in continuous plasma arc welding
Tsirkas et al. Evaluation of distortions in laser welded shipbuilding parts using local-global finite element approach
Yang et al. Identifying optimal process parameters in deep penetration laser welding by adopting Hierarchical-Kriging model
CN109359337B (zh) 一种可变光斑与粉末颗粒尺度耦合作用的仿真方法
Satyanarayana et al. Numerical simulation of the processes of formation of a welded joint with a pulsed Nd: YAG laser welding of ZR–1% NB alloy
CN103279629A (zh) 激光深熔点焊的能量利用率数值计算方法
CN103279630B (zh) 激光深熔点焊匙孔动态成形数值模拟方法
Chen et al. Numerical analysis of the effects of weld parameters on distortions and residual stresses in butt welded steel plates
Akella et al. Heat flux for welding processes: model for laser weld
CN117236139B (zh) 一种风电塔筒焊接残余应力预测方法
Nascimento et al. A literature review in heat source thermal modeling applied to welding and similar processes
Ku et al. 3D residual stress simulation of an excavate and weld repair mockup
Niu et al. Study on stress asymmetry of pipeline welded joint in mountainous areas with large slopes
Bin-Mansoor et al. Laser pulse heating of steel surface: consideration of phase-change process
Mandal et al. Fusion zone and HAZ prediction through 3-D simulation of welding thermal cycle
Farahani et al. Thermal energy storage control using phase change materials in a rectangular energy storage chamber with metal foam gradient and magnetic field
Zhang et al. Optimization of welding scheme for side beam of low pressure cylinder frame base on local-global mapping model
Narayanareddy et al. Finite element modeling of TIG welding for 316L stainless steel plate using Sysweld
Yamashita et al. Numerical simulation of laser welding processes with CIP finite volume method
CN108038288A (zh) 一种纳米核壳粒子温度场的预测方法及***

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Zhang Peilei

Inventor after: Chu Zhentao

Inventor after: Yang Gansheng

Inventor after: Ma Kai

Inventor after: Zhang Wenrui

Inventor after: Lu Yunlong

Inventor after: Shi Haichuan

Inventor after: Yu Zhishui

Inventor before: Chu Zhentao

Inventor before: Yang Gansheng

Inventor before: Ma Kai

Inventor before: Zhang Wenrui

Inventor before: Lu Yunlong

Inventor before: Shi Haichuan

Inventor before: Yu Zhishui

Inventor before: Zhang Peilei

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180724

CF01 Termination of patent right due to non-payment of annual fee