CN114273671A - 一种双光束激光粉末床熔融模拟仿真方法 - Google Patents
一种双光束激光粉末床熔融模拟仿真方法 Download PDFInfo
- Publication number
- CN114273671A CN114273671A CN202111519158.6A CN202111519158A CN114273671A CN 114273671 A CN114273671 A CN 114273671A CN 202111519158 A CN202111519158 A CN 202111519158A CN 114273671 A CN114273671 A CN 114273671A
- Authority
- CN
- China
- Prior art keywords
- laser
- metal
- powder bed
- ray
- temperature
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 239000000843 powder Substances 0.000 title claims abstract description 51
- 238000004088 simulation Methods 0.000 title claims abstract description 23
- 230000004927 fusion Effects 0.000 title claims abstract description 16
- 238000007639 printing Methods 0.000 claims abstract description 36
- 238000010521 absorption reaction Methods 0.000 claims abstract description 34
- 230000008569 process Effects 0.000 claims abstract description 33
- 238000012544 monitoring process Methods 0.000 claims abstract description 22
- 230000001052 transient effect Effects 0.000 claims abstract description 20
- 239000002184 metal Substances 0.000 claims description 42
- 230000008859 change Effects 0.000 claims description 19
- 238000002844 melting Methods 0.000 claims description 15
- 230000008018 melting Effects 0.000 claims description 15
- 239000012071 phase Substances 0.000 claims description 15
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 239000007769 metal material Substances 0.000 claims description 9
- 230000008021 deposition Effects 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 7
- 229910001338 liquidmetal Inorganic materials 0.000 claims description 7
- 238000013459 approach Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000004134 energy conservation Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 5
- 238000010309 melting process Methods 0.000 claims description 5
- 239000007790 solid phase Substances 0.000 claims description 5
- 125000004432 carbon atom Chemical group C* 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000009835 boiling Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 3
- 238000001704 evaporation Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000007791 liquid phase Substances 0.000 claims description 3
- 235000011837 pasties Nutrition 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000002310 reflectometry Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000007499 fusion processing Methods 0.000 abstract 1
- 238000005457 optimization Methods 0.000 abstract 1
- 238000003908 quality control method Methods 0.000 abstract 1
- 230000007547 defect Effects 0.000 description 4
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007648 laser printing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Landscapes
- Powder Metallurgy (AREA)
Abstract
本发明公开了一种双光束激光粉末床熔融模拟仿真方法，包括步骤一：建立双光束激光粉末床熔融打印过程模型；步骤二：根据步骤一建立的模型，控制激光运动，并计算、监测激光能量的吸收情况；步骤三：根据步骤一建立的模型，计算、监测构件温度场、流场的瞬态特征及空间分布；步骤四：根据步骤一建立的模型，计算、监测熔池和沉积道的瞬态特征及空间分布。利用本发明的模拟仿真方法，可以对多激光束粉末床熔融过程进行虚拟打印及工艺过程优化，有望为多激光束粉末床熔融打印成形与质量控制提供技术依据。
Description
技术领域
本发明属于激光增材制造领域，具体涉及一种多激光打印方法，尤其涉及一种双激光粉末床熔融过程的模拟仿真方法。
背景技术
激光粉末床熔融(LPBF)是一种高精度、高性能的金属零件增材制造技术。该技术非常适用于薄壁、复杂内腔和内部流道等复杂形状部件的整体制造，具有广泛的应用前景。为了满足航空航天和能源领域大型零件高效制造的迫切需求，激光粉末床熔融正在向多激光束制造方向发展。多激光粉末床熔融(ML-PBF)是一种更复杂的增材制造技术，包含多个单激光束模块以及扫描区域的可控拼接。
该技术可以极大地提高打印尺寸和生产率，同时保持其优越的精度特性。然而，多激光粉末床熔融过程的重叠区域容易出现成形异常、内部孔隙等缺陷。与单激光粉末床熔融制备的样品相比，多激光粉末床熔融制备的样品中重叠区域的表面精度、相对密度甚至力学性能在一定程度上可能有所下降。仅依托实验试错法探索费事费力，且难以对打印过程中的激光-物质作用机制以及温度场、流场及沉积几何特征演化过程进行全面研究。因此，需要结合高置信度多物理场耦合模拟仿真，对多激光粉末床熔融关键物理现象、成形过程与打印异常结果进行高分辨率剖析。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足，采用一种双光束激光粉末床熔融模拟仿真方法，对熔池、沉积道进行高时间分辨率与高空间分辨率监测，包括温度分布、内部流体运动与沉积轮廓形貌情况。
为实现上述技术目的，本发明采用如下方案：
一种双光束激光粉末床熔融模拟仿真方法，包括：
步骤一：建立双光束激光粉末床熔融打印过程模型；
步骤二：根据步骤一建立的模型，控制激光运动，并计算、监测激光能量的吸收情况；
步骤三：根据步骤一建立的模型及步骤二对激光运动的控制与激光能量吸收情况计算结果，进一步计算、监测构件温度场、流场的瞬态特征及空间分布；
步骤四：根据步骤一至步骤三双光束激光粉末床熔融过程建模、激光能量吸收率与温度场、流场模拟结果，进一步计算、监测熔池和沉积道的瞬态特征及空间分布。
具体地，所述步骤一具体为：
耦合求解连续性方程、质量守恒方程、能量守恒方程和VOF相方程，其中连续性方程和动量守恒方程：
上式中，
为流速，ρ为密度，μ为动态粘度，P为压力,
为重力加速度，
为网格单元体心的位置矢量，
为流体流动时所受到的力，分别是马兰戈尼力、蒸汽反冲压力、表面张力以及糊状区域阻尼力，具体表述如下:
上式中
为表面张力温度系数，T为温度，
为单位法向矢量，ρ1和ρ2分别为固态金属与气体的密度，
为界面项，引入该变量是为了将界面力转换为体积力，P0为压力，Hv为汽化潜热，M为摩尔质量，Tb为金属材料沸点，R为气体常数，γ为金属的表面张力系数，k为金属/气体界面的曲率，Kc为渗透系数为渗透系数，f为熔融金属的体积分数，为避免温度接近固相温度时零作为除数而引入一个极小值e0，fl由以下方程给出：
上式(7)中Tl和Ts分别是材料液相线温度和固相线温度，
能量守恒方程：
Qh＝Qh1+Qh2 (9)
上式中Cp为比热容，Qh1与Qh2定义为1号激光与2号激光的热源，两束激光作用在金属球体粉末项表面，Qh为两束激光热输入，Ql是两束激光在打印过程中辐射、对流和蒸发造成的热损失，即：
上式中辐射热损失、对流热损失、蒸发热损失分别表示为：
qc＝-hc(T-Tref) (12)
上式中σb为Stefan-Boltzmann常数，ε是发射率，hc是对流传热系数，Tref为环境温度，源项SLatent用于计算固液相转变过程中与潜热有关的能量变化，由下式定义：
SLatent＝SpT+Sc (14)
其中Hf为相变潜热，fl′为前一时间步的迭代值，温度相关的液态金属体积分数函数F及F-1定义为：
F-1＝f(Tl-Ts)+Ts (18)
制造过程中金属材料气液界面会随时间和空间的变化而变化，求解VOF相方程计算自由界面的动态形貌：
定义α为金属相的体积分数，若α＝1，则此处内部为金属区域，若α＝0，则单元内部为气体区域，若0＜α＜1即为气体与金属混合区域，在金属/气体界面上，热物性参数由加权函数给出：
ρ＝αρ1(1-α)ρ2 (20)
k＝αk1+(1-α)k2 (21)
Cp＝αCp1+(1-α)Cp2 (22)
μ＝αμ1+(1-α)μ2 (23)
其中k1和k2分别为金属与气体的热导率，Cp1和Cp2分别为金属与气体的比热容，μ1和μ2分别为金属与气体的粘度。
具体地，所述步骤二具体为：
两个激光束相互靠近，直到相遇，之后这两束激光束继续沿各自的扫描方向移动，采用射线追踪法计算激光束的实时能量吸收与反射情况，
射线追踪法将式(9)中的两束激光热源分别离散为若干个子光束，即射线，变量包括射线的初始坐标
功率Pi、方向向量
阳反射频率fi，射线的初始坐标
满足：
上式(24)中x0和y0分别为各光束中心的x，y坐标，zi初始化为计算域的z方向最大坐标，激光功率初始化为：
其中f为能量分布分数，激光束的初始方向为竖直向下，因此激光射线的方向向量
先定义为(0，0，-1)，
为单位向量，反射频率fi初始化为0，激光照射到粉体表面，
会发生相应的变化，任意表面网格与射线之间的垂直距离为：
其中
即射线坐标指向某单元格坐标的向量，其中(xm，ym，zm)为任意表面单元格的坐标，
与入射射线的角度θ1为：
当某入射激光射线照射金属某表面位置时，其部分功率被粉体吸收，其余被反射，被反射的射线向量更新为
反射射线的方向：
反射射线的频率：
反射射线的功率：
其中，
为界面法向量，ξR为自由表面的反射率，其依赖于入射射线和表面法向量之间的角度θ2，即：
上式中，nR和kR分别表示金属材料在激光能量吸收过程中光线折射率的实部和虚部，相应的表面网格的热源为：
各束激光吸收率定义为吸收功率与总光束功率的比值，即：
其中，A为激光吸收率，qlaser，m为第m个单元的计算热源，Ncell为数值网格的总数，Pi为第i条激光的功率，Nlaser为激光束的总数。
具体地，所述步骤二中：两束激光扫描路径在或不在同一条线上，并且方向相同或相反。
具体地，所述步骤二中：在打印过程中两个激光束的参数相同或不同。
具体地，所述步骤二中：采用射线追踪法计算激光束多重反射与能量吸收率，并监测从打印开始到打印结束全过程的激光能量吸收情况。
具体地，所述步骤三中：根据步骤二计算得出的激光能量吸收情况，计算并监测从打印开始到结束全过程中的温度场变化过程及其空间分布特征。
具体地，所述步骤三中：对双光束激光粉末床熔池模型求解，获得熔池液态金属的流动速度场，计算并监测出打印过程中流场的瞬态特征及空间分布。
具体地，所述步骤四中：采用VOF方法捕获液/气相的自由界面，针对每一时间步，依据粉体与基材的固相线温度及VOF界面特征，提取熔池三维轮廓数据，分析并监测打印过程中熔池的形貌变化。
具体地，所述步骤四中：基于凝固后熔池液态金属形成的沉积形貌，计算出沉积道特征，并实现沉积道特征的演化过程监测。
相比于现有技术，本发明具有以下有益效果：
(1)本发明充分利用先进计算机模拟、仿真方法，深入研究双光束激光粉末床熔融打印过程，揭示打印过程中熔池、沉积道、温度场和流场的瞬态特征及空间分布规律，有望为多激光粉末床融合实现控形一致性和统一性提供科学依据；
(2)本发明可以对熔池、沉积道进行高时间分辨率与空间分辨率监测，包括温度范围、轮廓形貌和内部流体运动情况，可用于深入分析打印缺陷及形成机理，控制缺陷提升打印质量。
附图说明
图1为双光束激光粉末床熔融相向打印瞬态过程模型；
图2为1号激光能量吸收率瞬态变化图；
图3为2号激光能量吸收率瞬态变化图；
图4为打印过程中温度场随时间的变化；
图5为打印过程中流场随时间的变化；
图6为打印结束后沉积道三维模型；
图7为打印过程中熔池随时间的变化；
图8为打印过程中熔池尺寸的瞬态变化图；
图9为打印过程中熔池整个表面积的瞬态变化图；
图10为打印过程中熔池整个体积的瞬态变化图。
具体实施方式
以下结合附图对本发明的实施例作进一步详细描述。
一种双光束激光粉末床熔融模拟仿真方法，包括：
步骤一：建立双光束激光粉末床熔融瞬态过程中相向打印模型。
本发明的实例采用的相关参数为：激光功率P＝500W，重叠区长度d＝1mm，激光扫描速度v＝1000mm/s。本发明基于开源计算流体力学库OpenFOAM开发了一种热流体模型以模拟仿真激光粉末床熔融制造过程中熔池的运动情况。耦合求解连续性方程、质量守恒方程、能量守恒方程和VOF相方程。其中连续性方程和动量守恒方程：
上式中，
为流速，ρ为密度，μ为动态粘度，P为压力，
为重力加速度，
为网格单元体心的位置矢量。
为流体流动时所受到的力，分别是马兰戈尼力、蒸汽反冲压力、表面张力以及糊状区域阻尼力，具体表述如下：
上式中
为表面张力温度系数，T为温度，
为单位法向矢量，ρ1和ρ2分别为固态金属与气体的密度，
为界面项，引入该变量是为了将界面力转换为体积力。P0为压力，Hv为汽化潜热，M为摩尔质量，Tc为金属材料沸点，R为气体常数。γ为金属的表面张力系数，k为金属/气体界面的曲率。Kc为渗透系数，f为熔融金属的体积分数，为避免温度接近固相温度时零作为除数而引入一个极小值e0。fl由以下方程给出：
上式中Tl和Ts分别是材料液相线温度和固相线温度。
能量守恒方程：
Qh＝Qh1+Qh2 (9)
上式中Cp为比热容，Qh1与Qh2定义为1号激光与2号激光的热源，两束激光作用在金属球体粉末项表面，Qh为两束激光热输入，Ql是两束激光在打印过程中辐射、对流和蒸发造成的热损失，即：
上式中辐射热损失、对流热损失、蒸发热损失分别表示为：
qc＝-hc(T--Tref) (12)
上式中σb为Stefan-Boltzmann常数，ε是发射率，hc是对流传热系数，Tref为环境温度。源项SLatent用于计算固液相转变过程中与潜热有关的能量变化，由下式定义：
SLatent＝SpT+Sc (14)
其中Hf为相变潜热，fl′为前一时间步的迭代值。温度相关的液态金属体积分数函数F及F-1定义为：
F-1＝f(Tl-Ts)+Ts (18)
制造过程中金属材料气液界面会随时间和空间的变化而变化，求解VOF相方程计算自由界面的动态形貌：
定义α为金属相的体积分数，若α＝1，则此处内部为金属区域，若α＝0，则单元内部为气体区域，若0＜α＜1即为气体与金属混合区域。在金属/气体界面上，热物性参数由加权函数给出：
ρ＝αρ1(1-α)ρ2 (20)
k＝αk1+(1-α)k2 (21)
Cp＝αCp1+(1-α)Cp2 (22)
μ＝αμ1+(1-α)μ2 (23)
上式中k1和k2分别为金属与气体的热导率，Cp1和Cp2分别为金属与气体的比热容，μ1和μ2分别为金属与气体的粘度。
步骤二：根据步骤一建立的模型，控制激光运动，并计算、监测激光能量的吸收情况。
如图1所示，两个激光束相互靠近，直到相遇，之后这两束激光束继续沿各自的扫描方向移动。本发明模拟了材料与热源之间的实时能量吸收与反射情况，采用射线追踪法计算激光束的实时能量吸收与反射情况。
射线追踪法将式(9)中的两束激光热源分别离散为若干个子光束(n个)，即射线，主要变量包括射线的初始坐标
功率Pi、方向向量
和反射频率fi。射线的初始坐标
满足：
上式中x0和y0分别为各光束中心的x，y坐标，zi初始化为计算域的z方向最大坐标。激光功率初始化为：
上式中f为能量分布分数，激光束的初始方向为竖直向下，因此激光射线的方向向量
先定义为(0，0，-1)。
为单位向量，反射频率fi初始化为0。激光照射到粉体表面，
会发生相应的变化。任意表面网格与射线之间的垂直距离为：
上式中
即射线坐标指向某单元格坐标的向量。其中(xm，ym，zm)为任意表面单元格的坐标。
与入射射线的角度θ1为：
当某入射激光射线照射金属某表面位置时，其部分功率被粉体吸收，其余被反射。被反射的射线向量更新为
反射射线的方向：
反射射线的频率：
反射射线的功率：
上式中，
为界面法向量，ξR为自由表面的反射率，其依赖于入射射线和表面法向量之间的角度θ2(θ2应该小于90°)，即：
上式中，nR和kR分别表示金属材料在激光能量吸收过程中光线折射率的实部和虚部。同时，相应的表面网格的热源为：
各束激光吸收率定义为吸收功率与总光束功率的比值，即：
上式中，A为激光吸收率，qlaser，m为第m个单元的计算热源，Ncell为数值网格的总数，Pi为第i条激光的功率，Nlaser为激光束的总数。如图2和图3所示，通过两个激光束相向运动模拟仿真结果，计算出在此过程中粉体对激光能量的吸收情况。
步骤三：根据步骤一建立的模型及步骤二对激光运动的控制与激光能量吸收情况计算结果，进一步计算、监测构件温度场、流场的瞬态特征及空间分布。
如图4和图5所示，模拟了两束激光在打印过程中接近、交汇、远离的过程，通过耦合求解方程(1)，(2)，(8)，(19)，获得构件温度场及熔池内部熔融金属的流场的变化情况。
步骤四：根据步骤一至步骤三双光束激光粉末床熔融过程建模、激光能量吸收率与温度场、流场模拟结果，进一步计算、监测熔池和沉积道的瞬态特征及空间分布。
如图7所示，在步骤一建立双光束激光粉末床熔融打印过程模型的基础上得出了熔池瞬态特征。如图8至图10所示，依据所用金属固相线温度和VOF界面特征，提取熔池三维轮廓数据，包括熔池深度、宽度、长度和高度。在此基础上，获得沉积道的瞬态特征及空间分布(如图6所示)，即左侧、中间(重叠区域)和右侧的最大高度和宽度。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。
Claims (10)
1.一种双光束激光粉末床熔融模拟仿真方法，其特征在于，包括：
步骤一：建立双光束激光粉末床熔融打印过程模型；
步骤二：根据步骤一建立的模型，控制激光运动，并计算、监测激光能量的吸收情况；
步骤三：根据步骤一建立的模型及步骤二对激光运动的控制与激光能量吸收情况计算结果，进一步计算、监测构件温度场、流场的瞬态特征及空间分布；
步骤四：根据步骤一至步骤三双光束激光粉末床熔融过程建模、激光能量吸收率与温度场、流场模拟结果，进一步计算、监测熔池和沉积道的瞬态特征及空间分布。
2.根据权利要求1所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤一具体为：
耦合求解连续性方程、质量守恒方程、能量守恒方程和VOF相方程，其中连续性方程和动量守恒方程：
上式中，
为流速，ρ为密度，μ为动态粘度，P为压力，
为重力加速度，
为网格单元体心的位置矢量，
为流体流动时所受到的力，分别是马兰戈尼力、蒸汽反冲压力、表面张力以及糊状区域阻尼力，具体表述如下：
上式中
为表面张力温度系数，T为温度，
为单位法向矢量，ρ1和ρ2分别为固态金属与气体的密度，
为界面项，引入该变量是为了将界面力转换为体积力，P0为压力，Hv为汽化潜热，M为摩尔质量，Tb为金属材料沸点，R为气体常数，γ为金属的表面张力系数，k为金属/气体界面的曲率，Kc为渗透系数为渗透系数，f为熔融金属的体积分数，为避免温度接近固相温度时零作为除数而引入一个极小值e0，fl由以下方程给出：
上式(7)中Tl和Ts分别是材料液相线温度和固相线温度，
能量守恒方程：
Qh＝Qh1+Qh2 (9)
上式中Cp为比热容，Qh1与Qh2定义为1号激光与2号激光的热源，两束激光作用在金属球体粉末项表面，Qh为两束激光热输入，Ql是两束激光在打印过程中辐射、对流和蒸发造成的热损失，即：
上式中辐射热损失、对流热损失、蒸发热损失分别表示为：
qc＝-hc(T-Tref) (12)
上式中σb为Stefan-Boltzmann常数，ε是发射率，hc是对流传热系数，Tref为环境温度，源项SLatent用于计算固液相转变过程中与潜热有关的能量变化，由下式定义：
SLatent＝SpT+Sc (14)
其中Hf为相变潜热，fl′为前一时间步的迭代值，温度相关的液态金属体积分数函数F及F-1定义为：
F-1＝f(Tl-Ts)+Ts (18)
制造过程中金属材料气液界面会随时间和空间的变化而变化，求解VOF相方程计算自由界面的动态形貌：
定义α为金属相的体积分数，若α＝1，则此处内部为金属区域，若α＝0，则单元内部为气体区域，若0＜α＜1即为气体与金属混合区域，在金属/气体界面上，热物性参数由加权函数给出：
ρ＝αρ1(1-α)ρ2 (20)
k＝αk1+(1-α)k2 (21)
Cp＝αCp1+(1-α)Cp2 (22)
μ＝αμ1+(1-α)μ2 (23)
其中k1和k2分别为金属与气体的热导率，Cp1和Cp2分别为金属与气体的比热容，μ1和μ2分别为金属与气体的粘度。
3.根据权利要求2所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤二具体为：
两个激光束相互靠近，直到相遇，之后这两束激光束继续沿各自的扫描方向移动，采用射线追踪法计算激光束的实时能量吸收与反射情况，
射线追踪法将式(9)中的两束激光热源分别离散为若干个子光束，即射线，变量包括射线的初始坐标
功率Pi、方向向量
和反射频率fi，射线的初始坐标
满足：
上式(24)中x0和y0分别为各光束中心的x，y坐标，zi初始化为计算域的z方向最大坐标，激光功率初始化为：
其中f为能量分布分数，激光束的初始方向为竖直向下，因此激光射线的方向向量
先定义为(0，0，-1)，
为单位向量，反射频率fi初始化为0，激光照射到粉体表面，
会发生相应的变化，任意表面网格与射线之间的垂直距离为：
其中
即射线坐标指向某单元格坐标的向量，其中(xm，ym，zm)为任意表面单元格的坐标，
与入射射线的角度θ1为：
当某入射激光射线照射金属某表面位置时，其部分功率被粉体吸收，其余被反射，被反射的射线向量更新为
反射射线的方向：
反射射线的频率：
反射射线的功率：
其中，
为界面法向量，ξR为自由表面的反射率，其依赖于入射射线和表面法向量之间的角度θ2，即：
上式中，nR和kR分别表示金属材料在激光能量吸收过程中光线折射率的实部和虚部，相应的表面网格的热源为：
各束激光吸收率定义为吸收功率与总光束功率的比值，即：
其中，A为激光吸收率，qlaser，m为第m个单元的计算热源，Ncell为数值网格的总数，Pi为第i条激光的功率，Nlaser为激光束的总数。
4.根据权利要求3所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤二中：两束激光扫描路径在或不在同一条线上，并且方向相同或相反。
5.根据权利要求3所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤二中：在打印过程中两个激光束的参数相同或不同。
6.根据权利要求3所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤二中：采用射线追踪法计算激光束多重反射与能量吸收率，并监测从打印开始到打印结束全过程的激光能量吸收情况。
7.根据权利要求3所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤三中：根据步骤二计算得出的激光能量吸收情况，计算并监测从打印开始到结束全过程中的温度场变化过程及其空间分布特征。
8.根据权利要求3所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤三中：对双光束激光粉末床熔池模型求解，获得熔池液态金属的流动速度场，计算并监测出打印过程中流场的瞬态特征及空间分布。
9.根据权利要求3所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤四中：采用VOF方法捕获液/气相的自由界面，针对每一时间步，依据粉体与基材的固相线温度及VOF界面特征，提取熔池三维轮廓数据，分析并监测打印过程中熔池的形貌变化。
10.根据权利要求9所述的双光束激光粉末床熔融模拟仿真方法，其特征在于，所述步骤四中：基于凝固后熔池液态金属形成的沉积形貌，计算出沉积道特征，并实现沉积道特征的演化过程监测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111519158.6A CN114273671A (zh) | 2021-12-13 | 2021-12-13 | 一种双光束激光粉末床熔融模拟仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111519158.6A CN114273671A (zh) | 2021-12-13 | 2021-12-13 | 一种双光束激光粉末床熔融模拟仿真方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114273671A true CN114273671A (zh) | 2022-04-05 |
Family
ID=80871862
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111519158.6A Pending CN114273671A (zh) | 2021-12-13 | 2021-12-13 | 一种双光束激光粉末床熔融模拟仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114273671A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117428210A (zh) * | 2023-12-20 | 2024-01-23 | 中国商用飞机有限责任公司 | 多激光选区熔融搭接方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007147221A1 (en) * | 2006-06-20 | 2007-12-27 | Katholieke Universiteit Leuven | Procedure and apparatus for in-situ monitoring and feedback control of selective laser powder processing |
| WO2018140019A1 (en) * | 2017-01-26 | 2018-08-02 | Siemens Product Lifecycle Management Software Inc. | System and method for adaptive domain reduction for thermo-structural simulation of additive manufacturing process |
| CN109359337A (zh) * | 2018-09-14 | 2019-02-19 | 华中科技大学 | 一种可变光斑与粉末颗粒尺度耦合作用的仿真方法 |
| CN111112621A (zh) * | 2020-01-22 | 2020-05-08 | 南京理工大学 | 一种激光定向能量沉积熔池形貌尺寸预测与监控方法 |
| CN111283192A (zh) * | 2020-01-22 | 2020-06-16 | 南京理工大学 | 一种激光粉末床熔融增材制造熔池监测与孔隙控制方法 |
| US20200341452A1 (en) * | 2019-04-23 | 2020-10-29 | Dassault Systems Simulia Corp | Machine learning with fast feature generation for selective laser melting print parameter optimization |
| WO2021015782A1 (en) * | 2019-07-25 | 2021-01-28 | Siemens Industry Software Inc. | Image-based defect detections in additive manufacturing |
| CN113239647A (zh) * | 2021-06-01 | 2021-08-10 | 南京航空航天大学 | 一种三维脉冲激光深熔焊动力学仿真数值建模方法 |
| CN113779793A (zh) * | 2021-09-10 | 2021-12-10 | 电子科技大学 | 一种基于射线追踪的激光选区熔化的热源建模方法 |
-
2021
- 2021-12-13 CN CN202111519158.6A patent/CN114273671A/zh active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007147221A1 (en) * | 2006-06-20 | 2007-12-27 | Katholieke Universiteit Leuven | Procedure and apparatus for in-situ monitoring and feedback control of selective laser powder processing |
| WO2018140019A1 (en) * | 2017-01-26 | 2018-08-02 | Siemens Product Lifecycle Management Software Inc. | System and method for adaptive domain reduction for thermo-structural simulation of additive manufacturing process |
| CN109359337A (zh) * | 2018-09-14 | 2019-02-19 | 华中科技大学 | 一种可变光斑与粉末颗粒尺度耦合作用的仿真方法 |
| US20200341452A1 (en) * | 2019-04-23 | 2020-10-29 | Dassault Systems Simulia Corp | Machine learning with fast feature generation for selective laser melting print parameter optimization |
| WO2021015782A1 (en) * | 2019-07-25 | 2021-01-28 | Siemens Industry Software Inc. | Image-based defect detections in additive manufacturing |
| CN111112621A (zh) * | 2020-01-22 | 2020-05-08 | 南京理工大学 | 一种激光定向能量沉积熔池形貌尺寸预测与监控方法 |
| CN111283192A (zh) * | 2020-01-22 | 2020-06-16 | 南京理工大学 | 一种激光粉末床熔融增材制造熔池监测与孔隙控制方法 |
| CN113239647A (zh) * | 2021-06-01 | 2021-08-10 | 南京航空航天大学 | 一种三维脉冲激光深熔焊动力学仿真数值建模方法 |
| CN113779793A (zh) * | 2021-09-10 | 2021-12-10 | 电子科技大学 | 一种基于射线追踪的激光选区熔化的热源建模方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117428210A (zh) * | 2023-12-20 | 2024-01-23 | 中国商用飞机有限责任公司 | 多激光选区熔融搭接方法 |
| CN117428210B (zh) * | 2023-12-20 | 2024-03-08 | 中国商用飞机有限责任公司 | 多激光选区熔融搭接方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | A new ray tracing heat source model for mesoscale CFD simulation of selective laser melting (SLM) | |
| Khamidullin et al. | Modeling of the effect of powder parameters on laser cladding using coaxial nozzle | |
| Cook et al. | Simulation of melt pool behaviour during additive manufacturing: Underlying physics and progress | |
| Fürstenau et al. | Generating virtual process maps of SLM using powder-scale SPH simulations | |
| He et al. | Transport phenomena during direct metal deposition | |
| CN111283192B (zh) | 一种激光粉末床熔融增材制造熔池监测与孔隙控制方法 | |
| CN111112621B (zh) | 一种激光定向能量沉积熔池形貌尺寸预测与监控方法 | |
| CN114273671A (zh) | 一种双光束激光粉末床熔融模拟仿真方法 | |
| Wang et al. | Characterization of laser-powder interaction and particle transport phenomena during laser direct deposition of W‒Cu composite | |
| Gu et al. | The transient heat conduction MPM and GIMP applied to isotropic materials | |
| CN112199881A (zh) | 一种直接金属沉积增材模拟方法与*** | |
| Minaie et al. | Analysis of flow patterns and solidification phenomena in the die casting process | |
| Han et al. | Study on a multifield coupling mechanism and a numerical simulation method of a pulsed laser deposition process from a disk laser | |
| Magana et al. | Multiphysics modeling of thermal behavior, melt pool geometry, and surface topology during laser additive manufacturing | |
| Ebrahimi et al. | Revealing the effects of laser beam shaping on melt pool behaviour in conduction-mode laser melting | |
| Piscopo et al. | Mesoscale modelling of laser powder-based directed energy deposition process | |
| Guan et al. | Manipulation of three-dimensional Richtmyer-Meshkov instability by initial interfacial principal curvatures | |
| Jhang et al. | Systematic modeling approach for analyzing the powder flow and powder energy absorptivity in direct energy deposition system | |
| Wiederkehr et al. | Acquisition and optimization of three-dimensional spray footprint profiles for coating simulations | |
| Zhan et al. | A novel finite element model for simulating residual stress in laser melting deposition | |
| Jia et al. | A comprehensive comparison of modeling strategies and simulation techniques applied in powder-based metallic additive manufacturing processes | |
| Watari et al. | Two-fluid model to simulate metal powder bed fusion additive manufacturing | |
| Duan et al. | Numerical investigation of gas entrapment in metal additive manufacturing using a quasi-sharp-interface particle method | |
| Ji | Thermochemical non-equilibrium effects of shock reflection hysteresis in steady supersonic flows | |
| Dorussen et al. | An efficient ray tracing methodology for the numerical analysis of powder bed additive manufacturing processes |
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 |