CN111046535A - 一种激光加工热分布计算方法 - Google Patents
一种激光加工热分布计算方法 Download PDFInfo
- 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
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 39
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- 239000000758 substrate Substances 0.000 claims abstract description 62
- 230000004907 flux Effects 0.000 claims abstract description 7
- 238000010521 absorption reaction Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 8
- 238000002844 melting Methods 0.000 claims description 7
- 230000008018 melting Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 5
- 229910000619 316 stainless steel Inorganic materials 0.000 claims description 4
- 238000000034 method Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 238000003754 machining Methods 0.000 claims 3
- 230000001678 irradiating effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000003698 laser cutting Methods 0.000 description 1
- 238000010147 laser engraving Methods 0.000 description 1
- 238000010330 laser marking Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000006552 photochemical reaction Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
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;
(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:其中,热导率k=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),基板的熔点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所示。
步骤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),基板的熔点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;
(4)将激光光束垂直照射在基板表面,并移动激光光束通过表面热源边界条件进行加载,满足:s∈Ω;激光照射区域Ω外的边界与空气接触,激光光束照射的外界条件为:其中,热导率κ=D·ρc,D为扩散率,Ω为激光照射区域范围,hc为基板表面散热系数;
(5)确定三维热导模型,并对所述三维热导模型进行转换计算;所述三维热导模型满足:
其中,t为激光光束照射时间;
(6)调用Matlab软件进行数值实验。
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。
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)
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)
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部队 | 一种基于热传导逆问题的激光强度反演方法 |
-
2019
- 2019-11-25 CN CN201911166561.8A patent/CN111046535B/zh active Active
Patent Citations (8)
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)
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)
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 |