[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN111046535A - 一种激光加工热分布计算方法 - Google Patents

一种激光加工热分布计算方法 Download PDF

Info

Publication number
CN111046535A
CN111046535A CN201911166561.8A CN201911166561A CN111046535A CN 111046535 A CN111046535 A CN 111046535A CN 201911166561 A CN201911166561 A CN 201911166561A CN 111046535 A CN111046535 A CN 111046535A
Authority
CN
China
Prior art keywords
substrate
laser
thermal conductivity
laser beam
heat
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
CN201911166561.8A
Other languages
English (en)
Other versions
CN111046535B (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.)
Jinan University
Original Assignee
Jinan University
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 Jinan University filed Critical Jinan University
Priority to CN201911166561.8A priority Critical patent/CN111046535B/zh
Publication of CN111046535A publication Critical patent/CN111046535A/zh
Application granted granted Critical
Publication of CN111046535B publication Critical patent/CN111046535B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明涉及激光加工技术领域,具体公开了一种激光加工热分布计算方法,包括(1)选定基板并定义基板的物性参数;(2)设定基板的初始温度,在基板上建立三维坐标系,确定基板的初始条件;(3)定义激光光束的物理参数,设定激光光束的能量分布以高斯分布热流密度的形式加载;(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,确定激光光束照射的外界条件等步骤。本发明的激光加工热分布计算方法,采用激光光束以高斯分布热流密度的形式加载在基板上,建立了三维的热导模型,推算出基板的热分布情况,相比于现有技术中积分平均值的计算方式,大大减少了模型误差,能够快速计算出更加精确的基板热源分布情况。

Description

一种激光加工热分布计算方法
技术领域
本发明涉及激光加工技术领域,尤其涉及一种激光加工热分布计算方法。
背景技术
根据激光束与材料相互作用的机理,大体可将激光加工分为激光热加工和光化学反应加工两类,激光热加工是指利用激光束投射到材料表面产生的热效应来完成加工过程,包括激光焊接、激光雕刻切割、表面改性、激光镭射打标、激光钻孔和微加工等。
人们为了简化激光照射在材料表面的热源表达式,通常采用了热源照射区域上热分布的积分平均值取代这一区域的真正热源分布(参见2016年虞钢等编的《激光先进制造技术及其应用》,国防工业出版社),这样就会产生模型误差,从而无法计算出精准的热源分布情况。
发明内容
针对现有技术中的技术问题,本发明提供一种激光加工热分布计算方法。
一种激光加工热分布计算方法,包括以下步骤:
(1)选定基板并定义基板的物性参数,基板的物性参数至少包括厚度X、表面半径R、密度ρ、比热容c、热导率κ和熔点K;
(2)设定基板的初始温度为Ta,并在基板上建立三维坐标系,且基板的初始条件为:T(x,y,z)|t=0=Ta
(3)定义激光光束的物理参数,包括激光功率P、光束半径rb,设定激光光束的能量分布以高斯分布热流密度的形式加载,满足:
Figure BDA0002287605580000011
其中α为基板的吸收系数;
(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:
Figure BDA0002287605580000021
s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:
Figure BDA0002287605580000022
其中,热导率k=D·ρc,D为扩散率,Ω为激光照射区域范围,hc为基板表面散热系数;
(5)确定三维热导模型,并对三维热导模型进行转换计算;三维热导模型满足:
Figure BDA0002287605580000023
其中,t为激光光束照射时间;
(6)调用Matlab软件进行数值实验。
进一步的,步骤(5)中对三维热导模型进行转换计算包括:
(501)展开计算为:
Figure BDA0002287605580000024
(502)将三维热导模型转换为圆柱坐标形式模型:
Figure BDA0002287605580000031
进一步的,步骤(6)包括:
(601)对实验参数进行赋值,实验参数包括基板的材质、尺寸、物性参数、吸收系数α、初始温度Ta以及激光光束的物理参数;
(602)将数值填入Matlab软件中进行计算并输出仿真结果。
进一步的,步骤(601)中,还包括设定激光照射时间t,以及时间步长Δt。
进一步的,步骤(601)中基板的材质为316不锈钢,基板的厚度X为0.3mm,基板的表面半径R为0.03mm,基板的密度ρ为8×103kg/m3,基板的比热容c为500J/(kg·K),基板的热导率κ为21.5W/(m·K),基板的熔点K为1673K,基板的吸收系数α为1,基板的初始温度Ta为20℃;
激光功率P为200W,光束半径rb为0.5mm。
进一步的,激光照射时间t为0.001s,时间步长Δt为0.0001s。
本发明实施例的激光加工热分布计算方法,采用激光光束以高斯分布热流密度的形式加载在基板上,建立了三维的热导模型,再根据激光光束传播的特性引入圆柱坐标,将三维问题转化为二维问题,推算出基板的热分布情况,相比于现有技术中积分平均值的计算方式,大大减少了模型误差,能够快速计算出更加精确的基板热源分布情况。
附图说明
为了更清楚的说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见的,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
图1为本发明实施例的激光加工热分布计算方法的步骤流程图;
图2为本发明实施例的激光加工热分布计算方法的三维坐标示意图;
图3为本发明实施例的激光加工热分布计算方法的实验仿真图。
具体实施方式
下面将结合本发明中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通的技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例,都属于本发明的保护范围。
本发明提供一种激光加工热分布计算方法,如图1所示,包括以下步骤:
步骤S101:选定基板并定义基板的物性参数,基板的物性参数至少包括厚度X、表面半径R、密度ρ、比热容c、热导率κ和熔点K。
密度ρ、比热容c、热导率κ和熔点K均与基板所选的材质有关,本实施例中的基板可根据不同需求选用不同材质。
步骤S102:设定基板的初始温度为Ta,并在基板上建立三维坐标系,确定基板的初始条件为:T(x,y,z)|t=0=Ta
x,y,z分别表示三维坐标系的x轴、y轴、z轴,具体建立如图2所示。
步骤S103:定义激光光束的物理参数,包括激光功率P、光束半径rb,设定激光光束的能量分布以高斯分布热流密度的形式加载,满足:
Figure BDA0002287605580000041
Figure BDA0002287605580000042
其中α为基板的吸收系数。
本发明实施例将激光光束的能量分布设定为高斯分布,即满足:
Figure BDA0002287605580000043
所以激光光束的能量分布以高斯分布热流密度的形式加载,满足
Figure BDA0002287605580000051
步骤S104:将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:
Figure BDA0002287605580000052
s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:
Figure BDA0002287605580000053
其中,热导率κ=D·ρc,D为基板的扩散率,Ω为激光照射区域范围,hc为基板表面散热系数。
本步骤确定出激光光束加载的边界条件以及外界条件,分别为
Figure BDA0002287605580000054
Figure BDA0002287605580000055
s∈Ω,
Figure BDA0002287605580000056
步骤S105:确定三维热导模型,并对三维热导模型进行转换计算;三维热导模型满足:
Figure BDA0002287605580000057
其中,t为激光光束照射时间。
结合步骤S101至S104,采用所有限定条件确定出三维热导模型,对该三维热导模型进行转换计算,具体为:
步骤S1051:展开计算为:
Figure BDA0002287605580000058
步骤S1052:将三维热导模型转换为圆柱坐标形式模型:
Figure BDA0002287605580000061
步骤S106:调用Matlab软件进行数值实验。
具体的,步骤S106包括:
步骤S1061:对实验参数进行赋值,实验参数包括基板的材质、尺寸、物性参数、吸收系数α、初始温度Ta以及激光光束的物理参数;
步骤S1062:将数值填入Matlab软件中进行计算并输出仿真结果。
具体的,本实施例的步骤S1061中还包括设定激光照射时间t,以及时间步长Δt。
具体的,本实施例的步骤S1061中,对实验参数进行赋值可采取:基板的材质为316不锈钢,基板的厚度X为0.3mm,基板的表面半径R为0.03mm,基板的密度ρ为8×103kg/m3,基板的比热容c为500J/(kg·K),基板的热导率κ为21.5W/(m·K),基板的熔点K为1673K,基板的吸收系数α为1,基板的初始温度Ta为20℃;激光功率P为200W,光束半径rb为0.5mm,激光照射时间t为0.001s,时间步长Δt为0.0001s。
将以上数值填入Matlab软件中进行计算并输出仿真结果,如图3所示,为本发明实施例的基板经过激光光束照射后的温度分布情况,由该实验结果可知,316不锈钢材质的基板只有在激光作用的区域内,即半径r为0的附近区域,温度才会有明显地变化(紫色部分),而其他离激光作用区域相对较远的区域(蓝色部分)的温度则基本不变,保持为初始温度,分析此实验结果的原因,这种现象主要是由于基板的热导率较低导致的,从而导致了基板表面和基板背面的温差较大。
本发明实施例的激光加工热分布计算方法,采用激光光束以高斯分布热流密度的形式加载在基板上,建立了三维的热导模型,再根据激光光束传播的特性引入圆柱坐标,将三维问题转化为二维问题,推算出基板的热分布情况,相比于现有技术中积分平均值的计算方式,大大减少了模型误差,能够快速计算出更加精确的基板热源分布情况。
以上借助具体实施例对本发明做了进一步描述,但是应该理解的是,这里具体的描述,不应理解为对本发明的实质和范围的限定,本领域内的普通技术人员在阅读本说明书后对上述实施例做出的各种修改,都属于本发明所保护的范围。

Claims (6)

1.一种激光加工热分布计算方法,其特征在于,包括以下步骤:
(1)选定基板并定义基板的物性参数,所述基板的物性参数至少包括厚度X、表面半径R、密度ρ、比热容c、热导率κ和熔点K;
(2)设定所述基板的初始温度为Ta,并在所述基板上建立三维坐标系,且所述基板的初始条件为:T(x,y,z)|t=0=Ta
(3)定义激光光束的物理参数,包括激光功率P、光束半径rb,设定激光光束的能量分布以高斯分布热流密度的形式加载,满足:
Figure FDA0002287605570000011
其中α为基板的吸收系数;
(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:
Figure FDA0002287605570000012
s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:
Figure FDA0002287605570000013
其中,热导率κ=D·ρc,D为扩散率,Ω为激光照射区域范围,hc为基板表面散热系数;
(5)确定三维热导模型,并对所述三维热导模型进行转换计算;所述三维热导模型满足:
Figure FDA0002287605570000014
其中,t为激光光束照射时间;
(6)调用Matlab软件进行数值实验。
2.如权利要求1所述的一种激光加工热分布计算方法,其特征在于,所述步骤(5)中对所述三维热导模型进行转换计算包括:
(501)展开计算为:
Figure FDA0002287605570000021
(502)将三维热导模型转换为圆柱坐标形式模型:
Figure FDA0002287605570000022
3.如权利要求2所述的一种激光加工热分布计算方法,其特征在于,所述步骤(6)包括:
(601)对实验参数进行赋值,所述实验参数包括所述基板的材质、尺寸、物性参数、吸收系数α、初始温度Ta以及所述激光光束的物理参数;
(602)将数值填入Matlab软件中进行计算并输出仿真结果。
4.如权利要求3所述的一种激光加工热分布计算方法,其特征在于,所述步骤(601)中,还包括设定激光照射时间t,以及时间步长Δt。
5.如权利要求3所述的一种激光加工热分布计算方法,其特征在于,所述步骤(601)中所述基板的材质为316不锈钢,所述基板的厚度X为0.3mm,所述基板的表面半径R为0.03mm,所述基板的密度ρ为8×103kg/m3,所述基板的比热容c为500J/(kg·K),所述基板的热导率k为21.5W/(m·K),所述基板的熔点K为1673K,所述基板的吸收系数α为1,所述基板的初始温度Ta为20℃;
所述激光功率P为200W,光束半径rb为0.5mm。
6.如权利要求4所述的一种激光加工热分布计算方法,其特征在于,所述激光照射时间t为0.001s,所述时间步长Δt为0.0001s。
CN201911166561.8A 2019-11-25 2019-11-25 一种激光加工热分布计算方法 Active CN111046535B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911166561.8A CN111046535B (zh) 2019-11-25 2019-11-25 一种激光加工热分布计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911166561.8A CN111046535B (zh) 2019-11-25 2019-11-25 一种激光加工热分布计算方法

Publications (2)

Publication Number Publication Date
CN111046535A true CN111046535A (zh) 2020-04-21
CN111046535B CN111046535B (zh) 2022-12-09

Family

ID=70233333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911166561.8A Active CN111046535B (zh) 2019-11-25 2019-11-25 一种激光加工热分布计算方法

Country Status (1)

Country Link
CN (1) CN111046535B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112032613A (zh) * 2020-08-11 2020-12-04 南京嘉视信电子有限公司 一种能够快速散热的隧道灯及其设计方法
CN113369694A (zh) * 2021-04-25 2021-09-10 西安交通大学 一种双光束耦合激光增材成形方法及装置
CN114248022A (zh) * 2021-11-19 2022-03-29 深圳市裕展精密科技有限公司 用于板材的激光打孔方法、激光控制系统及存储介质
CN115803200A (zh) * 2020-06-24 2023-03-14 西门子股份公司 预热用于增材制造的构建板

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101788459A (zh) * 2010-02-08 2010-07-28 暨南大学 一种准连续方式的分光波长组合方法
US20100292965A1 (en) * 2009-05-18 2010-11-18 Bekir Sami Yilbas Method of modeling phase changes due to laser pulse heating
CN102955881A (zh) * 2012-10-30 2013-03-06 温州大学 一种计算集成电路芯片焊点热疲劳失效概率的方法
DE102012015167A1 (de) * 2012-07-31 2013-03-21 Daimler Ag Verfahren zur Fehlererkennung in Bauteilen
EP2754524A1 (de) * 2013-01-15 2014-07-16 Corning Laser Technologies GmbH Verfahren und Vorrichtung zum laserbasierten Bearbeiten von flächigen Substraten unter Verwendung einer Laserstrahlbrennlinie
CN109530696A (zh) * 2018-12-28 2019-03-29 天津镭明激光科技有限公司 一种基板作为零件一部分的激光选区熔化成形方法
CN109883443A (zh) * 2019-02-19 2019-06-14 北京工业大学 一种线结构光传感器空间姿态标定方法
CN110057444A (zh) * 2019-05-29 2019-07-26 中国人民解放军32027部队 一种基于热传导逆问题的激光强度反演方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100292965A1 (en) * 2009-05-18 2010-11-18 Bekir Sami Yilbas Method of modeling phase changes due to laser pulse heating
CN101788459A (zh) * 2010-02-08 2010-07-28 暨南大学 一种准连续方式的分光波长组合方法
DE102012015167A1 (de) * 2012-07-31 2013-03-21 Daimler Ag Verfahren zur Fehlererkennung in Bauteilen
CN102955881A (zh) * 2012-10-30 2013-03-06 温州大学 一种计算集成电路芯片焊点热疲劳失效概率的方法
EP2754524A1 (de) * 2013-01-15 2014-07-16 Corning Laser Technologies GmbH Verfahren und Vorrichtung zum laserbasierten Bearbeiten von flächigen Substraten unter Verwendung einer Laserstrahlbrennlinie
CN109530696A (zh) * 2018-12-28 2019-03-29 天津镭明激光科技有限公司 一种基板作为零件一部分的激光选区熔化成形方法
CN109883443A (zh) * 2019-02-19 2019-06-14 北京工业大学 一种线结构光传感器空间姿态标定方法
CN110057444A (zh) * 2019-05-29 2019-07-26 中国人民解放军32027部队 一种基于热传导逆问题的激光强度反演方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHENGHAO YUAN ET AL.: "A three-dimensional unconditionally stable ADI-FDTD method in the cylindrical coordinate system", 《 IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》 *
MIHAIOANE ET AL.: "Mathematical modeling of two-photon thermal fields in laser–solid interaction", 《OPTICS & LASER TECHNOLOGY》 *
丁宇等: "连续波1064nm激光辐照铝金属的热分布研究", 《光电技术应用》 *
朱大伟等: "三维柱坐标ADI-FDTD算法及其CPML实现", 《微波学报》 *
李艳丽等: "激光工艺参数对45钢淬透性能的影响", 《新技术新工艺》 *
王金艳: "空心光泵浦半导体激光器的热特性研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115803200A (zh) * 2020-06-24 2023-03-14 西门子股份公司 预热用于增材制造的构建板
CN112032613A (zh) * 2020-08-11 2020-12-04 南京嘉视信电子有限公司 一种能够快速散热的隧道灯及其设计方法
CN113369694A (zh) * 2021-04-25 2021-09-10 西安交通大学 一种双光束耦合激光增材成形方法及装置
CN113369694B (zh) * 2021-04-25 2022-06-07 西安交通大学 一种双光束耦合激光增材成形方法及装置
CN114248022A (zh) * 2021-11-19 2022-03-29 深圳市裕展精密科技有限公司 用于板材的激光打孔方法、激光控制系统及存储介质
CN114248022B (zh) * 2021-11-19 2023-11-10 富联裕展科技(深圳)有限公司 用于板材的激光打孔方法、激光控制系统及存储介质

Also Published As

Publication number Publication date
CN111046535B (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
CN111046535A (zh) 一种激光加工热分布计算方法
CN104148636B (zh) 一种控制金属零件增材制造热变形的扫描路径生成方法
CN108889948B (zh) 一种用于薄壁件增材制造的分区扫描方法
CN111872395B (zh) 镍基合金单晶或定向凝固零件的制备方法
CN112380752B (zh) 利用预测金属薄板焊接热处理数值改进其焊接工艺的方法
US20170364058A1 (en) Controlled thin wall thickness of heat exchangers through modeling of additive manufacturing process
JP2014508668A (ja) 三次元体の製造方法
CN113414407B (zh) 一种镍基合金增材制造方法和镍基合金零件
WO2019049981A1 (ja) 積層造形物の解析方法及び積層造形物の解析装置、並びに積層造形物の製造方法及び積層造形物の製造装置
CN108971267B (zh) 一种微变形板材激光矫平方法
CN105598448A (zh) 一种金属材料激光3d打印原位预热温度的控制方法
JP2017211977A (ja) データ作成装置、3次元積層システム、設定方法及びプログラム
CN105127426A (zh) 一种三维结构组合烧结加工方法
CN117798382A (zh) 基于有限元及轨迹追踪的增材制造场量宏观特性预测方法
Gautam et al. Laser forming of mild steel sheets using different surface coatings
US20180099358A1 (en) Metallic Sleeve For Reducing Distortion In Additive Manufacturing
JP4739507B2 (ja) 3次元造形装置
CN113523302B (zh) 选区激光熔化成型镁合金烧损的抑制方法
Sun et al. Melting cell based compensated design method for improving dimensional accuracy of additively manufactured thin channels
CN107552785A (zh) 一种3d打印、热处理一体化加工方法
JP5343052B2 (ja) 構造物解析方法、プログラムおよび解析装置
CN205834237U (zh) 带大颗粒氧化物回收导向装置的3d打印设备
CN105799175A (zh) 一种维持工作缸整体温度均衡的装置及其控制方法
CN113423525B (zh) 用于在增材制造中加热基材的方法
CN112699592A (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
GR01 Patent grant
GR01 Patent grant