CN111046535B - 一种激光加工热分布计算方法 - Google Patents
一种激光加工热分布计算方法 Download PDFInfo
- Publication number
- CN111046535B CN111046535B CN201911166561.8A CN201911166561A CN111046535B CN 111046535 B CN111046535 B CN 111046535B CN 201911166561 A CN201911166561 A CN 201911166561A CN 111046535 B CN111046535 B CN 111046535B
- Authority
- CN
- China
- Prior art keywords
- substrate
- laser
- laser beam
- thermal conductivity
- 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.)
- Active
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明涉及激光加工技术领域,具体公开了一种激光加工热分布计算方法,包括(1)选定基板并定义基板的物性参数;(2)设定基板的初始温度,在基板上建立三维坐标系,确定基板的初始条件;(3)定义激光光束的物理参数,设定激光光束的能量分布以高斯分布热流密度的形式加载;(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,确定激光光束照射的外界条件等步骤。本发明的激光加工热分布计算方法,采用激光光束以高斯分布热流密度的形式加载在基板上,建立了三维的热导模型,推算出基板的热分布情况,相比于现有技术中积分平均值的计算方式,大大减少了模型误差,能够快速计算出更加精确的基板热源分布情况。
Description
技术领域
本发明涉及激光加工技术领域,尤其涉及一种激光加工热分布计算方法。
背景技术
根据激光束与材料相互作用的机理,大体可将激光加工分为激光热加工和光化学反应加工两类,激光热加工是指利用激光束投射到材料表面产生的热效应来完成加工过程,包括激光焊接、激光雕刻切割、表面改性、激光镭射打标、激光钻孔和微加工等。
人们为了简化激光照射在材料表面的热源表达式,通常采用了热源照射区域上热分布的积分平均值取代这一区域的真正热源分布(参见2016年虞钢等编的《激光先进制造技术及其应用》,国防工业出版社),这样就会产生模型误差,从而无法计算出精准的热源分布情况。
发明内容
针对现有技术中的技术问题,本发明提供一种激光加工热分布计算方法。
一种激光加工热分布计算方法,包括以下步骤:
(1)选定基板并定义基板的物性参数,基板的物性参数至少包括厚度X、表面半径R、密度ρ、比热容c、热导率κ和熔点K1;
(2)设定基板的初始温度为Ta,并在基板上建立三维坐标系,且基板的初始条件为:T(x,y,z)|t=0=Ta;
(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:其中,热导率κ=D·ρc,D为扩散率,Ω为激光照射区域范围,hc为基板表面散热系数;
(5)确定三维热导模型,并对三维热导模型进行转换计算;三维热导模型满足:
其中,t为激光光束照射时间;
(6)调用Matlab软件进行数值实验。
进一步的,步骤(5)中对三维热导模型进行转换计算包括:
(501)展开计算为:
(502)将三维热导模型转换为圆柱坐标形式模型:
进一步的,步骤(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),基板的熔点K1为1673K,基板的吸收系数α为1,基板的初始温度Ta为20℃;
激光功率P为200W,光束半径rb为0.5mm。
进一步的,激光照射时间t为0.001秒,时间步长Δt为0.0001秒。
本发明实施例的激光加工热分布计算方法,采用激光光束以高斯分布热流密度的形式加载在基板上,建立了三维的热导模型,再根据激光光束传播的特性引入圆柱坐标,将三维问题转化为二维问题,推算出基板的热分布情况,相比于现有技术中积分平均值的计算方式,大大减少了模型误差,能够快速计算出更加精确的基板热源分布情况。
附图说明
为了更清楚的说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见的,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
图1为本发明实施例的激光加工热分布计算方法的步骤流程图;
图2为本发明实施例的激光加工热分布计算方法的三维坐标示意图;
图3为本发明实施例的激光加工热分布计算方法的实验仿真图。
具体实施方式
下面将结合本发明中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通的技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例,都属于本发明的保护范围。
本发明提供一种激光加工热分布计算方法,如图1所示,包括以下步骤:
步骤S101:选定基板并定义基板的物性参数,基板的物性参数至少包括厚度X、表面半径R、密度ρ、比热容c、热导率κ和熔点K1。
密度ρ、比热容c、热导率κ和熔点K1均与基板所选的材质有关,本实施例中的基板可根据不同需求选用不同材质。
步骤S102:设定基板的初始温度为Ta,并在基板上建立三维坐标系,确定基板的初始条件为:T(x,y,z)|t=0=Ta。
x,y,z分别表示三维坐标系的x轴、y轴、z轴,具体建立如图2所示。
步骤S104:将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:其中,热导率κ=D·ρc,D为基板的扩散率,Ω为激光照射区域范围,hc为基板表面散热系数。
步骤S105:确定三维热导模型,并对三维热导模型进行转换计算;三维热导模型满足:
其中,t为激光光束照射时间。
结合步骤S101至S104,采用所有限定条件确定出三维热导模型,对该三维热导模型进行转换计算,具体为:
步骤S1051:展开计算为:
步骤S1052:将三维热导模型转换为圆柱坐标形式模型:
步骤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),基板的熔点K1为1673K,基板的吸收系数α为1,基板的初始温度Ta为20℃;激光功率P为200W,光束半径rb为0.5mm,激光照射时间t为0.001秒,时间步长Δt为0.0001秒。
将以上数值填入Matlab软件中进行计算并输出仿真结果,如图3所示,为本发明实施例的基板经过激光光束照射后的温度分布情况,由该实验结果可知,316不锈钢材质的基板只有在激光作用的区域内,即半径r为0的附近区域,温度才会有明显地变化(紫色部分),而其他离激光作用区域相对较远的区域(蓝色部分)的温度则基本不变,保持为初始温度,分析此实验结果的原因,这种现象主要是由于基板的热导率较低导致的,从而导致了基板表面和基板背面的温差较大。
本发明实施例的激光加工热分布计算方法,采用激光光束以高斯分布热流密度的形式加载在基板上,建立了三维的热导模型,再根据激光光束传播的特性引入圆柱坐标,将三维问题转化为二维问题,推算出基板的热分布情况,相比于现有技术中积分平均值的计算方式,大大减少了模型误差,能够快速计算出更加精确的基板热源分布情况。
以上借助具体实施例对本发明做了进一步描述,但是应该理解的是,这里具体的描述,不应理解为对本发明的实质和范围的限定,本领域内的普通技术人员在阅读本说明书后对上述实施例做出的各种修改,都属于本发明所保护的范围。
Claims (5)
1.一种激光加工热分布计算方法,其特征在于,包括以下步骤:
(1)选定基板并定义基板的物性参数,所述基板的物性参数至少包括厚度X、表面半径R、密度ρ、比热容c、热导率κ和熔点K1;
(2)设定所述基板的初始温度为Ta,并在所述基板上建立三维坐标系,且所述基板的初始条件为:T(x,y,z)|t=0=Ta;
(3)定义激光光束的物理参数,包括激光功率P、光束半径rb,设定激光光束的能量分布以高斯分布热流密度的形式加载,满足:其中α为基板的吸收系数;其中,q(r)表示激光光束的能量分布,r表示基板上的点到坐标原点的距离;
(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:其中,热导率κ=D·ρc,D为扩散率,Ω为激光照射区域范围,hc为基板表面散热系数;
(5)确定三维热导模型,并对所述三维热导模型进行转换计算;所述三维热导模型满足:
其中,t为激光光束照射时间;
(6)调用Matlab软件进行数值实验;
所述步骤(5)中对所述三维热导模型进行转换计算包括:
(501)展开计算为:
(502)将三维热导模型转换为圆柱坐标形式模型:
2.如权利要求1所述的一种激光加工热分布计算方法,其特征在于,所述步骤(6)包括:
(601)对实验参数进行赋值,所述实验参数包括所述基板的材质、尺寸、物性参数、吸收系数α、初始温度Ta以及所述激光光束的物理参数;
(602)将数值填入Matlab软件中进行计算并输出仿真结果。
3.如权利要求2所述的一种激光加工热分布计算方法,其特征在于,所述步骤(601)中,还包括设定激光照射时间t,以及时间步长Δt。
4.如权利要求3所述的一种激光加工热分布计算方法,其特征在于,所述步骤(601)中所述基板的材质为316不锈钢,所述基板的厚度X为0.3mm,所述基板的表面半径R为0.03mm,所述基板的密度ρ为8×103kg/m3,所述基板的比热容c为500J/(kg·K),所述基板的热导率κ为21.5W/(m·K),所述基板的熔点K1为1673K,所述基板的吸收系数α为1,所述基板的初始温度Ta为20℃;
所述激光功率P为200W,光束半径rb为0.5mm。
5.如权利要求3所述的一种激光加工热分布计算方法,其特征在于,所述激光照射时间t为0.001秒,所述时间步长Δt为0.0001秒。
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 CN111046535A (zh) | 2020-04-21 |
CN111046535B true 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) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112032613A (zh) * | 2020-08-11 | 2020-12-04 | 南京嘉视信电子有限公司 | 一种能够快速散热的隧道灯及其设计方法 |
CN113369694B (zh) * | 2021-04-25 | 2022-06-07 | 西安交通大学 | 一种双光束耦合激光增材成形方法及装置 |
CN114248022B (zh) * | 2021-11-19 | 2023-11-10 | 富联裕展科技(深圳)有限公司 | 用于板材的激光打孔方法、激光控制***及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101788459A (zh) * | 2010-02-08 | 2010-07-28 | 暨南大学 | 一种准连续方式的分光波长组合方法 |
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部队 | 一种基于热传导逆问题的激光强度反演方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8160846B2 (en) * | 2009-05-18 | 2012-04-17 | King Fahd University Of Petroleum & Minerals | Method of modeling phase changes due to laser pulse heating |
-
2019
- 2019-11-25 CN CN201911166561.8A patent/CN111046535B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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)
Title |
---|
A three-dimensional unconditionally stable ADI-FDTD method in the cylindrical coordinate system;Chenghao Yuan et al.;《 IEEE Transactions on Microwave Theory and Techniques》;20021210;第50卷(第10期);第2401-2405页 * |
Mathematical modeling of two-photon thermal fields in laser–solid interaction;MihaiOane et al.;《Optics & Laser Technology》;20040430;第36卷(第3期);第219-222页 * |
三维柱坐标ADI-FDTD算法及其CPML实现;朱大伟等;《微波学报》;20191024;第35卷(第5期);第8-12,39页 * |
激光工艺参数对45钢淬透性能的影响;李艳丽等;《新技术新工艺》;20110125(第1期);第50-51页 * |
空心光泵浦半导体激光器的热特性研究;王金艳;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20140815;第I135-132页 * |
连续波1064nm激光辐照铝金属的热分布研究;丁宇等;《光电技术应用》;20161215;第31卷(第6期);第28-32页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111046535A (zh) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111046535B (zh) | 一种激光加工热分布计算方法 | |
JP5712306B2 (ja) | 三次元体の製造方法 | |
EP3608040B1 (en) | Powder-bed additive manufacturing device | |
Pinkerton et al. | Modelling the geometry of a moving laser melt pool and deposition track via energy and mass balances | |
US10584923B2 (en) | Systems and methods for heat exchanger tubes having internal flow features | |
CN105598448B (zh) | 一种金属材料激光3d打印原位预热温度的控制方法 | |
CN112380752B (zh) | 利用预测金属薄板焊接热处理数值改进其焊接工艺的方法 | |
CN104190930A (zh) | 一种同质功能梯度材料及结构的激光增材制造方法 | |
CN112191849A (zh) | 基于温度分布的梯度多孔散热装置设计及增材制造方法 | |
CN112149330B (zh) | 一种风电塔筒油封平台焊接残余应力预测、焊接工艺优化方法 | |
CN105718690A (zh) | 基于时间和空间主动跟踪的激光3d打印熔池凝固行为数值模拟方法 | |
US20180099358A1 (en) | Metallic Sleeve For Reducing Distortion In Additive Manufacturing | |
JP2003321704A (ja) | 積層造形法およびそれに用いる積層造形装置 | |
CN110545940A (zh) | 基于粉末床增材制造工件的方法、为该前述方法建立校正参数的方法和用于该后述方法的计算机程序产品 | |
CN218873735U (zh) | 一种金属激光3d打印熔化成形装置 | |
Sun et al. | Melting cell based compensated design method for improving dimensional accuracy of additively manufactured thin channels | |
Singh et al. | Modelling of dimensional accuracy in precision investment casting using Buckingham’s Pi approach | |
CN108607990A (zh) | 一种高表面质量的快速增材制造方法 | |
CN115041710A (zh) | 一种多能束增材制造的三维温度场控制装置 | |
JP2007075885A (ja) | レーザ曲げ加工予測装置、レーザ曲げ加工予測プログラム、レーザ曲げ加工予測方法 | |
JP2019000905A (ja) | 角度付きろう付け継手を有する複合構成要素、クーポンろう付け方法および関連する記憶媒体 | |
JP6799083B2 (ja) | 熱伝達装置、ターボ機械ケーシング、および関連の記憶媒体 | |
JP7358821B2 (ja) | 積層造形装置および積層造形方法 | |
EP3705211A1 (en) | Method for heating a base material in additive manufacturing | |
US20190085719A1 (en) | Anti-deflection feature for additively manufactured thin metal parts and method of additively manufacturing thin metal parts |
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 |