首页/文章/ 详情

超声能场对激光增材制造熔池凝固过程的影响机制研究

1月前浏览837
这类软件好像有COMSOL,还有SYSWELD。貌似abaqus、ansys也有,不过我知识面比较窄还没关注。

选择用COMSOL Multiphysics软件对超声能场处理下激光熔丝增材制造熔池凝固过程进行仿真模拟。模拟主要包括两部分内容: 激光熔丝增材制造 A517Q 合金钢温度场仿真,探索成形过程中温度场的变化情况;对 不同超声能场施加方式下熔池凝固过程中熔池声压场和温度场进行模拟,分析超声能场对激光熔丝增材制造熔池凝固过程的影响机制。

超声能场辅助激光熔丝增材制造有限元模型的建立

为了进行有限元数值模拟,首先建立超声能场辅助激光熔丝增材制造相关模型,其中网格划分非常关键,决定了仿真模拟的计算精度和时间。与此同时,材料参数的正确选择对仿真结果的准确性也有较大影响,边界条件、热源模型和超声能场接触方式选择越贴近实际实验情况,模拟结果越准确。所以首先对几何模型建立、网格划分、材料 参数选择、热源模型选取和边界条件设置等过程进行说明。

几何模型及网格划分

首先建立增材制造成形过程几何模型,由于对增材制造过程模拟只是为了给后续熔池凝固过程研究提供熔池初始温度,所以沉积层截面形状对模拟结果几乎不存在影响,考虑到建模和有限元计算的便利性,将沉积层截面形状设置为矩形,如图(a)所示。 在对求解精度和求解时间进行综合考虑后,对网格划分选择了正常划分后局部加密,网格的疏密情况根据实际情况而定。整个沉积层以及与沉积层接触的基板部分对网格进行加密,其余部分网格尺寸较大,节省求解时间。完整网格包含49447个域单元、6766个边界单元和614个边单元,网格最小尺寸为0.04mm,如图(b)所示。
(a)模型;(b)网格划分
为了研究超声能场对熔池凝固过程的影响,重新进行了建模和网格划分。为了节省计算时间,对模型进行了简化,使用二维模型来模拟超声能场对熔池凝固过程的影响。 根据上超声设备和下超声设备实际作用位置,在模型中加入小长方体作为超声能场作用点,如图(a)所示。在保证计算精度的情况下,为节约模型计算时间,对整个模型进行分区域网格划分。首先对熔池部分进行网格划分,选择流体动力学作为网格划分标准,网格类型为自由三角形,最小网格尺寸为0.35μm,由于该部分为重点模拟区域, 所以对该区域网格进行细化处理,限制最大网格尺寸为0.1mm。对液态熔池与已凝固沉积层以及基板的边界层进行了重新定义,定义层数为2层,层厚根据调节因子进行调节, 调节因子限制为5,同样限制边界层最大网格尺寸为0.1mm。模型其他部分按照极细化 进行划分,网格类型选用自由三角形,最小网格尺寸为3μm,最终完整网格包含10673个域单元和577个边界单元。
(a)模型;(b)模型局部放大;(c)网格局部放大

材料参数

激光熔丝增材制造过程中,由于熔池温度在很短时间内存在剧烈变化,而材料的一些参数随着温度变化而变化,所以选择合理的材料参数非常重要。在 COMSOL Multiphysics 材料库中并没有本文所用的A517Q 合金钢,JMatPro软件可以通过输入材料成分才对材料的各种参数进行计算,所以使用JMatPro软件对A517Q合金钢密度、 杨氏模量、热导率和动力粘度等参数进行了计算,如图所示。其余参数保持恒定不 变,不随温度变化而变化,具体数值如表所示。
表A517Q 合金钢模拟所需材料参数
A517Q 合金钢材料参数随温度变化曲线:
(a)密度;(b)杨氏模量;(c)热导率;(d)动力 粘度

热源模型

在增材制造过程中,激光热源分布为中间多边缘少,所以选择平面高斯热源, 该热源分布函数与实际激光热源分布近乎相同,使用这种热源模型可以获得较为准确的 熔池温度场。其热流密度随着坐标与热源中心的距离呈现指数型分布,其方程式如式:

边界条件

使用 COMSOL Multiphysics 软件的全耦合方式进行超声能场对熔池凝固过程的模拟,其中包括固体和流体传热模块、层流模块以及固体力学模块。在进行温度场模拟时,主要对包括沉积层初始温度、热传导、热对流、热辐射等的边界条件进行设置。 通过激光熔丝增材制造单道单层温度场的模拟结果,确定熔池初始温度设置为2500K, 模型其他部分初始温度皆为293.15K,以此来模拟熔池凝固过程。热传导的由以下方程表示:
由于模拟过程中存在固-液相变,相变过程中材料会吸收或者释放热量,产生相变潜热,对凝固过程温度也有较大影响,在边界条件中用相变材料进行定义,给出了具体固相温度、液相温度和相变潜热具体数值,如上表所示。
在进行超声能场对熔池声压场模拟过程中,主要包括以下界面:上超声施加点,其加载载荷为正弦位移载荷A=A0sin2πft。其中超声频率f为20kHz,超声输出振幅A0可调整,最大为20μm。下超声施加点,其加载载荷也为正弦位移载荷A=A1sin2πft。其中 超声频率f为20kHz,超声输出振幅A1可调整,最大为16μm。基板底面与成形平台接 触面设置为低反射界面,模拟超声波传播到该位置时不再进行反射。固-液耦合界面,用来模拟超声波在固液界面的传播情况,进行耦合计算。设置声软化边界,因为超声从金 属熔体传播到空气中,金属熔体的声特性阻抗率大于空气,所以设置声软化边界来模拟熔池和空气的接触面。

超声能场不同参数对熔池凝固过程的影响

在研究超声能场不同参数对熔池凝固过程的影响机制之前,首先对激光熔丝增材制造单道单层温度场进行了模拟,以此来确定后续模拟时熔池的初始温度,模拟结果如图所示。从图(a)可以看出,激光光斑中心及其附近温度很高,但由于热源向前移 动速度很快,热源移动后,材料的温度降低很快,在很短的时间内,材料发生凝固,由液态变为固态。结合图(b)的温度曲线以及A517Q合金钢的固相温度和液相温度, 最终选择2500K作为熔池初始温度。
激光熔丝增材制造单道单层温度场模拟:
(a)三维温度场;(b)熔池内一点的温度-时间曲 线
因为实际成形过程中涉及到大量的物理化学变化,且超声能场对熔池凝固过程的影响是动态过程,本实验超声能场的周期为50 μs,熔池凝固过程虽然时间较短,但也以秒为单位,所以超声能场与熔池凝固过程的耦合较为复杂,计算步长小,计算量大。因此,为了降低数值计算的复杂程度,提高模拟的可行性并减少求解时间,在模拟过程中做出如下简化和假设:
(1)使用二维模型代替三维模型,在保证模拟结果准确度的基础上,大大缩短了数值计算时间,提高了计算效率。
(2)忽略增材制造过程中熔池运动行为,不考虑超声能场的引入对宏观温度上产生的影响。
(3)假设熔池形状保持不变,环境温度为293.15K,大气压力为1.01×105 Pa,忽略成形过程中保护气的影响。

上超声不同振幅对沉积层熔池声压场和温度场的影响

在激光熔丝增材制造单道单层温度场模拟的基础上,开展了超声能场对熔池凝固过程影响规律的研究,首先模拟了超声不同振幅对熔池凝固过程声压场的影响。模拟时间步长为1.0μs,从熔池温度2500K开始计算,直至整个熔池完全凝固为止。超声能场在熔池中传播时,熔池内部会出现拉应力和压应力交替出现的现象,熔池受到拉应力时会形成负压,当负压大于熔池内分子之间的结合力时,熔池内会生成大量气泡,即空化气泡。此时最大负压的峰值称为空化阈值,也就是说想要形成空化泡,熔池内的负压必须大于等于空化阈值。随后熔池受到压应力时会形成正压,空化泡在正压的作用下会闭合破碎,在空化泡的破碎时会产生强烈的对流和冲击波,导致局部范围内产生较大的温度梯度和压力梯度。这会引起枝晶破碎,增加过冷度,从而提高熔池中的形核率,进而细化晶粒。所以通过声压场的模拟来研究本实验条件下超声能场的空化效应作用情况。
因为超声波会在传播的过程中逐渐衰减,即超声实际振幅会随着距离的增加不断减小,所以为了研究不同超声输出振幅下熔池处超声能量的大小,对不同输出振幅下熔池处实际振幅进行了计算,如图所示。从图中可以看出,随着超声输出振幅的增加, 熔池处实际振幅大幅度提高,从2.11μm分别提高到3.56μm和5.26μm。说明超声的输出振幅影响较大,直接决定了熔池处实际振幅的大小。
不同超声输出振幅下熔池处实际振幅随时间变化曲线:
(a)输出振幅8 μm;(b)输出振幅 14 μm;(c)输出振幅20 μm
超声波产生的能量可以分为两部分,分别是声传播能和空化能。声传播能主要在介质中传播消耗,最后转化为成内能,空化能主要被空化气泡吸收并转化为其他形式的能量。空化气泡单位时间和空间内释放的能量称为空化强度,可以通过声强来反应空化强度的大小。熔池处声强的大小可以使用熔池处的实际振幅来计算获得,声强I计算为公式所示:
根据图的超声不同输出振幅下熔池处实际振幅大小,计算了熔池处声强,如图所示。对于轻金属熔体来说,发生超声空化对应的声强临界值约为100W/cm^2,在图中用黑色虚线进行标记。从图(a)可以看出,当超声输出振幅为8μm时,声强仅 为 95.08 W/cm^2,小于发生超声空化所需要的声强临界值,判断在熔池内超声空化现象 很弱,或者没有发生超声空化。当超声输出振幅增加到14μm时,声强提高到了288.52W/cm^2(b), 超过了超声空化所需要的声强阈值,判断熔池内发生了超声空化效 应。当超声输出振幅为20μm时,声强进一步提高,达到了580.33W/cm^2( c), 远大于超声空化的声强阈值,判断熔池内部超声空化效应较强且更为充分。
不同超声输出振幅下熔池处声强随时间变化曲线:
(a)输出振幅8 μm;(b)输出振幅 14 μm;(c)输出振幅 20 μm
除了声强可以反应超声的空化强度外,超声的声压也可以对空化效应的强弱进行定 量的分析,所以对超声空化的声压阈值进行了计算,其公式由表示,经过计算,A517Q合金钢的空化阈值PB为0.22MPa。
为了对超声能场施加后熔池内部声压的变化有更直观的了解,对超声不同输出振幅下熔池内部声压稳定时的声压云图进行了绘制,如图所示。图中空白区域代表该部分的声压大于超声空化的声压阈值,也就是说该部分可以发生超声空化效应。从图中可以看出,熔池内部声压呈现出从底部到顶部逐渐降低的趋势,且随着超声输出振幅的不断增大,可以发生超声空化的范围也不断增大。熔池内部最大声压值也不断增大,当超声输出振幅为8μm时,最大声压为0.22MPa,当超声输出振幅增加到14 μm时,最大 声压提高到0.41MPa,当超声输出振幅增加到20μm时,最大声压为0.52MPa,提高了约2.4倍。
不同超声输出振幅下熔池内部声压云图:
(a)输出振幅8 μm;(b)输出振幅14 μm;(c)输 出振幅20 μm
为了对不同超声输出振幅下熔池内部声压进行定量统计,取熔池内部十个超声周期内平均声压的变化来判断超声空化效应的发生情况。想要形成空化泡,熔池内的负压必须大于等于空化阈值(0.22 MPa)。从图可以看出,当超声输出振幅为8μm时,其内部声压变化趋于稳定后,负压为-0.1475MPa(a),没有达到超声空化的声压阈值,结合其声强值及声压云图,判断该参数下熔池内部超声空化效应较弱,对熔池影响 不大。当超声输出振幅为14μm时,声压变化趋于稳定后,负压为-0.3082MPa(b),超过了超声空化的声压阈值,熔池内部有空化泡生成,其正压为0.4983 MPa,可以 使产生的空化泡破裂,结合其声强值及声压云图,判断该参数下熔池内接近三分之一的区域都发生了超声空化效应。当超声输出振幅继续增大到20μm时,负压为-0.4591MPa (c),相比于超声输出振幅8μm时,其负压值提高了约2倍,且正压为0.6721MPa,更容易使空化泡发生破裂。综合以上分析,当超声输出振幅为20μm时,熔池内部超声空化效应更易发生且超声空化范围更大。
不同超声输出振幅下熔池内平均声压随时间变化曲线:
(a)输出振幅8 μm;(b)输出振幅 14 μm;(c)输出振幅20 μm
超声的声流效应对熔池的作用主要体现在对熔池温度场的影响上。超声波在熔体中传播时,超声能量会在熔体介质的粘度和声吸收作用下发生衰减,这促使声压随着超声波的传播,形成了一定的声压梯度,声压梯度的形成会使熔体内部发生快速流动,这种现象称为声流效应。声流效应使熔池内部热对流的速度可以提高5倍以上。熔池内部流体流动速度的提高可以对熔池存在一定的搅拌作用,起到破坏液/固界面,促进散热和加速颗粒分散等作用。
对超声不同输出振幅下熔池温度场进行了研究,熔池不同凝固时刻的温度场云图如图所示。从图中可以看出,在相同凝固时刻,超声不同输出振幅处理下,熔池内温度场有较大变化。当超声作用时间为0.01 s时,由于超声作用时间较短,声流效应对熔池温度场产生的变化较小,但熔池内最低温度明显随着超声输出振幅的增加而不断降低。 当超声作用时间提高到0.03s时,不同超声输出振幅的温度场出现明显变化。超声输出 振幅为8μm时,熔池内最高温度为2460K(a),超声输出振幅增加到14μm时, 熔池内最高温度降低到2420K(b),超声输出振幅为20μm时,熔池内最高温度降低到2350K(c)。当超声作用时间提高到0.05s时,随着超声输出振幅的提高, 熔池内温度降低幅度更大,当超声输出振幅为20μm时,熔池内最高温度仅为1780K, 接近完全凝固温度。
不同超声输出振幅下熔池凝固过程温度场:
(a)输出振幅8 μm;(b)输出振幅14 μm; (c)输出振幅20 μm
为了更直观的观察熔池凝固过程中固-液相变,绘制了超声不同输出振幅下固-液相变云图,如图所示。其中红色代表液相,蓝色代表固相,从图中可以看出,随着超声输出振幅的增大,熔池凝固速度加快,固相范围大幅度增加。在超声输出振幅20μm 作用0.05 s后,熔池超过一半已经转变为固相。
不同超声输出振幅下熔池凝固过程固-液相变:
(a)输出振幅8 μm;(b)输出振幅14 μm; (c)输出振幅20 μm
为了便于进行对比,分析不同超声输出振幅下熔池完全凝固所需时间,对熔池凝固过程中不同时刻的最高温度进行了统计,如图所示。从图中可以看出,当超声输出振幅增加时,同一时间,熔池内最高温度不断减小。当超声输出振幅为8μm时,熔池 完全凝固需要0.13s,超声输出振幅增大到14μm时,熔池完全凝固时间降低为0.08s, 超声输出振幅为20μm时,完全凝固时间仅为0.06s,熔池完全凝固时间大幅度降低。 这也从侧面证明,经过上超声处理后,熔池的冷却速度增大,冷却速度增大会导致过冷度增大,过冷度的增大会提高熔池的形核率,进而达成细化晶粒的目的。
不同超声输出振幅下熔池内最高温度随时间变化曲线:
(a)超声输出振幅8 μm;(b)超声 输出振幅14 μm;(c)超声输出振幅20 μm

下超声不同功率对沉积层熔池声压场和温度场的影响

为了研究下超声不同功率对熔池声压场和温度场的影响,首先模拟了不同功率下熔 池处的实际振幅,如图所示。从图中可以看出,随着超声功率的提高,熔池处实际振幅大幅度提高。当超声功率为200W时,熔池处实际振幅为4.25μm(a);当 超声功率增大到600W时,振幅提高为6.33μm(b);超声功率为1000W时, 熔池处实际振幅高达10.46μm(c),为总振幅(16μm)的66.2%,超声衰减较小。熔池处实际振幅更大,熔池受到超声能量的干扰更大。
不同超声功率下熔池处实际振幅随时间变化曲线:
(a)超声功率200 W;(b)超声功率600 W;(c)超声功率1000 W
根据上图的超声不同功率下熔池处实际振幅大小,计算了熔池处声强,如图所示。从下图中可以明显看出,下超声不同功率作用时,声强都远远超过了发生超声空化 所需要的声强阈值(100W/cm^2),即下超声作用时,熔池内部都存在超声空化效应。随 着超声功率的增加,声强不断增大,从419.67W/cm^2(200W),分别提高到933.77W/cm^2 (600 W)和2622.95W/cm^2(1000W)。下超声作用时,熔池声强值较大,说明其超声空化效应较易发生。
不同超声功率下熔池处声强随时间变化曲线:
(a)超声功率200 W;(b)超声功率600 W;(c)超声功率1000 W
对超声不同功率下熔池内部声压稳定时的声压云图进行了绘制,如图所示。图中空白区域可以发生超声空化效应。从图中可以看出,在下超声作用下,熔池内部声压分布并不均匀,在熔池的底部呈现出极大的负压值,熔池顶部呈现出较大的正压值,且随着超声功率的增大,熔池内正压值和负压值的绝对值都逐渐增大。当超声功率为200W时,最大负压值为-0.69MPa(a);超声功率增大到600W时,最大负压值为-0.92MPa(b);超声功率为1000W时,最大负压值进一步增大为-1.09MPa(c),且近乎整个熔池都可以发生超声空化效应。这是由于超声功率为1000 W时,初始声压较高,所以经过衰减后熔池处接收到的超声能量仍然较大。
不同超声功率下熔池内部声压云图:
(a)超声功率200 W;(b)超声功率600 W;(c)超 声功率1000 W
对熔池内部十个超声周期内平均声压的进行计算和分析,其曲线如图所示。从图中可以看出,当超声功率为 200W时,其内部声压变化趋于稳定后,负压为-0.3639MPa,正压为0.5541MPa(a),产生的空化泡可以在正压的作用下发生破裂,即 可以发生空化效应。结合其声压云图,判断该参数下熔池内接近二分之一的区域都发生 了超声空化效应。当超声功率为600W时,负压为-0.5705MPa,正压为0.7836MPa(b),符合超声空化效应发生条件,结合其声压云图,判断该参数下熔池内接近四分 之三的区域都可以发生超声空化效应。当超声功率继续增大到1000W时,负压为-1.0295MPa,正压为1.1737MPa(c),相比于超声功率200W时,其负压值提高了约 2.8 倍。综合以上分析,当超声功率为1000W时,熔池内部超声空化效应更易发生且近乎整个熔池范围内都存在超声空化效应。
不同超声功率下熔池内平均声压随时间变化曲线:
(a)超声功率200 W;(b)超声功率600 W;(c)超声功率1000 W
模拟了超声不同功率下熔池凝固过程的温度场的变化情况。熔池不同凝固时刻的温度场云图如图所示。当超声作用时间为0.01s时,此时超声作用时间较短,只有超声功率为1000 W时,温度存在明显变化。当超声作用时间提高到0.03s时,不同超声功率下温度场都出现明显变化。超声功率为200W时,熔池内低温区明显扩大,熔池内最低温度为1050K,超声功率增加到600W时,熔池内低温区进一步扩大,最低温度 降低为980K,超声功率增加到1000W时,熔池内低温区占比接近整个熔池的一半, 最低温度降低到849K。当超声作用时间提高到0.05s时,低温区进一步扩大且随着超 声功率的增加低温区面积不断增大,当超声功率为1000W时,熔池内最高温度为1560K,接近固-液转变温度,整个熔池温度大幅度降低。为了方便区分超声不同功率对固-液相变的影响,对熔池凝固过程固-液相变进行了模拟,如图所示。从图中可以看出, 当超声作用时间相同时,随着超声功率逐渐增大,固相所占面积越来越大。在超声作用 时间0.05s,超声功率为1000W时,几乎整个熔池内部全部由液相转变为固相。
不同超声功率下熔池凝固过程温度场:
(a)超声功率200 W;(b)超声功率600 W;(c) 超声功率1000 W
不同超声功率下熔池凝固过程固-液相变:
(a)超声功率200 W;(b)超声功率600 W; (c)超声功率1000 W
对不同超声功率下熔池凝固过程中不同时刻的最高温度进行了统计,如图 所 示。从图中可以看出,随着超声功率的增加,熔池内同一时刻最高温度不断降低。当超 声功率为200W时,熔池完全凝固需要0.11s,超声功率增大到600W时,完全凝固时 间降低为0.07s,超声功率为1000W时,完全凝固时间仅为0.05s,熔池完全凝固时间 大幅度降低。说明随着超声功率的增加,熔池冷却时间不断减小。综上所述,当超声功 率为1000W时,熔池的冷却速度最快,超声的声流效应最明显,对沉积层显微组织改 善效果最好。
不同超声功率下熔池内最高温度随时间变化曲线:
(a)超声功率200 W;(b)超声功率600 W;(c)超声功率1000 W

超声能场不同施加方式对熔池凝固过程的影响

超声不同施加方式对熔池内部声压场的影响

综合以上分析,结合上超声输出振幅20μm以及下超声功率1000W,对双超声耦合处理方式进行模拟,并比较了未经超声处理、上超声单独处理、下超声单独处理、双超声耦合处理四种超声施加方式对熔池凝固过程的影响。
首先对超声不同施加方式下熔池处实际振幅和声强进行了计算,如图所示。从图中可以看出,上超声处理时熔池处实际振幅和声强值相对较小。这是由于上超声处理 与熔池间隔为20mm,且是通过微锻造头与沉积层表面的接触来传播超声,接触面积较小,所以熔池接收到的超声总能量较小,因此上超声处理方式下熔池处实际振幅和声强值相对较小。下超声处理时熔池处实际振幅和声强值最大。这是由于本次模拟熔池凝固 过程选择是沉积层中间部位(L+S 区),下超声与熔池之间距离较小,且下超声是通过变幅杆与基板直接接触来传播超声,接触面积大,所以从熔池处接收到的超声总能量大,因此熔池处实际振幅和声强值比较大。双超声耦合处理时,熔池处实际振幅和声强值处在上超声和下超声处理之间。这是由于上超声和下超声共同作用时,超声波以正弦函数进行传播,可能上下超声存在部分抵消,导致双超声耦合处理下熔池处实际振幅和声强要比下超声单独处理时小。
超声不同施加方式下熔池处实际振幅和声强随时间变化曲线:
(a)实际振幅;(b)声强
对超声不同施加方式下熔池内部声压云图进行了绘制,如图所示。从图中可看出,当上超声单独作用时,最高声压值仅为0.52MPa,超声空化范围较小,只有熔池面积的五分之二可以发生超声空化效应。当下超声单独作用时,最高声压值为-1.09MPa, 是上超声单独作用时的2倍,且熔池内大部分区域可以发生超声空化效应。当双超声耦合作用时,最高声压值提高到了-1.73MPa,是上超声单独作用时的3.3倍,整个熔池范围内几乎都可以发生超声空化效应。说明双超声耦合作用下,超声空化范围进一步增大, 整个熔池都可以发生超声空化效应,进而使枝晶破碎更加彻底。
超声不同施加方式下熔池内部声压云图:
(a)上超声;(b)下超声;(c)双超声
对比了超声不同施加方式下熔池内平均声压随时间的变化,如图所示。从图中可以看出,在声压变化趋于稳定后,三种超声施加方式的负压值都大于超声空化的声压阈值,并且在正负压之间交替变化,从而不断重复空化气泡的产生和破裂,反映了超声对熔池的空化效应。
超声不同施加方式下熔池内平均声压随时间变化曲线

超声不同施加方式对熔池内部温度场的影响

超声不同施加方式下熔池凝固过程温度场云图如图所示。从图中可以看出,未经超声能场处理时,熔池内部的温度场分布与熔池的形状相类似,温度降低较慢,在冷却0.03s 后,熔池内最高温度为2460K,温度变化不大。在上超声处理后,温度降低速度增加,在冷却0.03s后,熔池内最高温度降低为2350K,相较未经超声处理熔池温度有了较大的降低,且低温区面积增大。在下超声处理0.03s后,熔池内最高温度降低为2390K,熔池内低温区面积明显扩大,最低温度降低到849K。在双超声耦合处理下,当冷却时间为0.02s 时,熔池温度就出现了极大地降低,此时熔池内最高温度为 1900K,且熔池内只有很小的区域温度高于1510K,当冷却时间增加到0.03s时,熔池内最 高温度为1300K,整个熔池温度迅速降低。
超声不同施加方式下熔池凝固过程温度场:
(a)无超声;(b)上超声;(c)下超声;(d) 双超声
对超声不同施加方式下熔池凝固过程中的固-液相变进行了分析,如图所示。 从图中可以看出,在冷却0.03s时,未经超声处理熔池和上超声处理下熔池基本还是以液相为主,下超声处理后熔池冷却速度加快,固相占到了熔池总面积的一半,双超声处理后,熔池完全凝固,全部转化为固相。
超声不同施加方式下熔池凝固过程固-液相变:
(a)无超声;(b)上超声;(c)下超声; (d)双超声
超声不同施加方式下熔池内最高温度随时间变化曲线如图所示。从图中可以看出,经过超声能场处理后,熔池完全凝固时间大大降低。未经超声处理熔池完全凝固时间为0.27s,上超声处理熔池完全凝固时间降低为0.06s,下超声处理熔池完全凝固时 间降低为0.05s,双超声耦合处理下熔池完全凝固时间降低为0.03s,说明双超声处理 后,熔池凝固速度有大幅度提高。
超声不同施加方式下熔池内最高温度随时间变化曲线:
(a)无超声;(b)上超声;(c)下 超声;(d)双超声
图为冷却时间0.03s时,熔池在超声不同施加方式下,从底部到顶部温度梯度和凝固速度变化曲线。从图中可以看出,经过超声能场处理后,熔池温度梯度明显降低。 上超声处理时,由于熔池接收到的能量较小,熔池内声流效应不明显,所以其温度梯度降低幅度较小。下超声处理和双超声耦合处理下,温度梯度大幅度降低。证明超声对熔 池的声流效应可以使熔池内温度场更加均匀,降低温度梯度。三种超声施加方式下,熔池凝固速度都有明显的提升,其中双超声耦合方式提升最大。低的温度梯度和高的凝固 速度有利于柱状晶向等轴晶转变,从而细化晶粒,提高沉积层力学性能。
超声不同施加方式下熔池从底部到顶部温度梯度和凝固速度变化曲线:
(a)温度梯度; (b)凝固速度

超声不同施加方式对凝固组织的影响

为了更直观的说明超声不同施加方式对凝固组织的影响,将熔池凝固过程中固液界 面处的温度梯度G和凝固速度R画在A517Q合金钢的凝固图中进行对比,如图所示。其中A517Q合金钢的凝固图是通过Hunt-Gäumann模型由公式绘制得到的。 从图中可以看出,未经超声处理样品,沉积层内以柱状枝晶为主,凝固组织一开始为柱 状枝晶,当凝固过程快结束时,存在部分柱状晶和等轴晶的混合组织。上超声处理时,熔池也是开始先生成柱状枝晶,但柱状枝晶只占沉积层的很小部分,后续生成大量等轴枝晶。下超声处理时,由于凝固速度的大幅度提高,沉积层凝固组织几乎全为柱状晶和等轴晶的混合组织。双超声耦合处理时,温度梯度降低,凝固速度大幅度提高,所以其最终凝固组织几乎全为等轴枝晶。这说明超声能场对熔池的声流效应可以明显降低熔池温度梯度,提高凝固速度,从而促进柱状晶向等轴晶转变。
超声不同施加方式下熔池凝固过程柱状晶向等轴晶转变对比

模拟与实验结果对比

模型有效性验证:
(a)无超声晶粒分布;(b)双超声晶粒分布;(c)凝固组织预测

超声能场对熔池凝固过程的影响机制

通过超声能场辅助下熔池凝固过程的模拟仿真,总结了超声能场对熔池凝固过程的影响机制,主要分为两个方面:空化效应和声流效应。
超声的空化效应对熔池的影响又分为引起形核率增加和引起枝晶破碎两种。首先当超声能量作用于液态熔池时,在极短的时间内(500μs),熔池内部声压会发生剧烈变化,正负压不断交替。当熔池内负压大于该材料的超声空化声压阈值时,就会生成空化泡,空化泡会随着负压的增加而膨胀,且会随着正压的增加而收缩。经过一个或多个压缩和膨胀循环后,空化泡最终会闭合破裂。空化泡的破裂会产生冲击波和强烈对流,使得熔池内局部区域产生高温高压,从而提高了过冷度,在高过冷度下熔池形核率增加,使得形核速率大于生长速率,所以生成了大量细小晶粒。其次由于空化泡膨胀、收缩、破裂的时间远小于熔池凝固时间,所以熔池内的枝晶会受到多次冲击波的影响,使得粗大的枝晶不断破碎,枝晶碎片又会作为新的形核点生成细小晶粒,抑制柱状枝晶的生长。且在这个过程中熔池中气泡更容易逸出,从而减少微孔等缺陷,提高沉积层的致密度。
超声的声流效应主要对熔池温度场产生影响,通过对熔池温度梯度和凝固速度的控 制,进而影响最终熔池凝固组织。超声能场作用于熔池时,由于超声波随着传播距离会逐渐衰减,进而形成声压梯度,这会加速熔池内的马兰戈尼流,使得熔池对流散热加快,熔池内温度场更加均匀,在降低熔池内的温度梯度的同时提高了凝固速度,根据Hunt Gäumann 模型,高凝固速度和低温度梯度有利于柱状晶向等轴晶转变,从而细化晶粒。 除此之外,熔池流动性的提高还可以使得沉积层显微组织更加均匀,降低了材料的各向异性,同时减少了未熔合缺陷,提高了沉积层的致密度。

结论

(1)在沉积层施加上超声处理后,超声能量可以通过沉积层传递到液态熔池中,并且随着超声输出振幅的增大,熔池处的实际振幅和声强不断增大,结合声压云图,超声空化范围也不断增大。当超声输出振幅为20μm时,熔池处实际振幅为5.26μm,声强 为580.33W/cm^2,远大于超声空化所需要的声强阈值(100 W/cm^2),超声空化范围接近熔池的三分之一。随着超声输出振幅的增大,熔池内最高温度也不断降低,冷却时间为 0.05s 时,在超声输出振幅20 μm的作用下,熔池内最高温度仅为1780K,熔池超过一半已经转化为固相,完全凝固时间为0.06s。
(2)在基板下方施加超声能场处理时,随着超声功率的提高,熔池处实际振幅和声强也不断增大,且比上超声处理时提高幅度更大,结合声压图,其超声空化范围也不断 增大。当超声功率为1000 W时,熔池处实际振幅高达10.46μm,声强为2622.95W/cm^2 远高于上超声处理时数值,且熔池内大部分范围都发生了超声空化效应。冷却时间为0.05s时,在超声功率1000W的作用下,熔池内最高温度为1560K,熔池基本已经完全凝固。
(3)经过双超声耦合处理后,熔池处实际振幅为8.95μm,声强为1825.57W/cm^2, 虽然其声强比下超声处理时小,但结合声压图,其超声空化范围比下超声处理时要大, 几乎整个熔池都可以发生超声空化效应。冷却时间为0.03s时,在双超声耦合作用下,熔池内最高温度为1300K,熔池已经完全凝固,熔池完全凝固时间大大降低。超声不同施加方式都促进了柱状晶向等轴晶转变,上超声处理对温度梯度的降低贡献更大,下超声处理对凝固速度的提高贡献更大,双超声耦合处理下,结合了上下超声的优点,所以其熔池凝固组织几乎全为等轴晶,晶粒得到明显的细化。

网络整理,仅限内部分享,禁止商用。
公众 号:机电君


来源:机电君
AbaqusComsol化学增材SYSWELDJMatProUM材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-12-26
最近编辑:1月前
ErNan.Chen🍃
硕士 | CAE工程师 即物而穷其理
获赞 32粉丝 49文章 66课程 0
点赞
收藏
作者推荐

Fluent如何优化汽车气动设计以提高性能

在深度学习与流体力学融合的背景下,科研边界不断拓展,创新成果层出不穷。从物理模型融合到复杂流动模拟,从数据驱动研究到流场智能分析,深度学习正以前所未有的力量重塑流体力学领域。目前在Nature和Science杂志上发表的深度学习驱动的流体力学方面的论文主要集中以下几个方面:1、流体力学方程的求解:利用深度学习模型来求解流体力学的基本方程,如纳维-斯托克斯方程。2、湍流模拟:应用深度学习技术来改进湍流模型,提高湍流模拟的准确性和效率。3、流场重建与超分辨率:使用深度学习算法对流场数据进行超分辨率重建,提升现有数据的分辨率,以更精确地模拟和分析流体流动。4、流动特征识别与分类:利用深度学习模型来识别和分类流体流动中的关键特征。5、流动控制与优化:应用深度学习进行流动控制策略的优化,以提高流体机械的性能。6、计算流体动力学(CFD)与机器学习的结合:将深度学习集成到传统的CFD软件中,以提高计算效率和精度。7、物理约束神经网络:开发满足物理守恒定律的神经网络模型。8、激波和边界层过渡:利用深度学习预测和分析流体中的激波以及边界层的过渡现象。9、实验数据与模拟数据的融合:使用深度学习来提高流体力学模型的预测能力。这些研究性成果不仅推动了流体力学领域的科学发展,也展示了深度学习技术在解决复杂流体动力学问题中的潜力。为促进科研人员、工程师及产业界人士对深度学习在流体力学领域应用的技术掌握,北京软研国际信息技术研究院特举办本次专题培训会议,承办方互动派(北京)教育科技有限公司,会议会务合作单位为北京中科万维智能科技有限公司,具体相关事宜通知如下:专题一(直播四天)基于机器学习深度学习驱动的流体力学流场重建技术与应用(详情内容点击上方名称查看)2024年11月16日-11月17日2024年11月23日-11月24日专题二(直播五天)基于OpenFOAM和深度学习驱动的流体力学计算与应用(详情内容点击上方名称查看)2024年10月26日-10月27日2024年11月01日-11月03日专题三(直播五天)基于Fluent和深度学习算法驱动的流体力学计算与应用(详情内容点击上方名称查看)2024年10月26日-10月27日2024年11月01日-11月03日适用人群流体力学相关领域的科研人员,力学、航空航天科学与工程、工业通用技术及设备、动力工程、船舶工业、建筑科学与工程、石油天然气工业、机械工业、汽车工业、环境科学与资源利用等领域的工程师,工业自动化、机器人、智能制造等相关行业从业者,跨领域研究人员。01培训讲师1 流场重建讲师国外某高校博士研究生, 研究方向集中于深度学习在流体动力学中的应用, 涵盖流场重建, 流场预测及流动控制等领域。过去三年内, 以第一作者及主要作者身份在 Journal of Fluid Mechanics (JFM), Physics of Fluids (POF), Physical Review Fluids (PRF), 以及 Nature Scientific Reports 等 SCI 期刊发表论文十余篇, 谷歌学术引用次数超过 300 次。2 OpenFOAM 讲师来自全球顶尖大学,具有丰富的流体力学工作经验,包括实验流体力学,计算流体动力学(CFD),近年来发表论文多篇。擅长领域:流体力学与人工智能的交叉科学,流场预测与重构,气动信息预测,基于深度强化学习的气动优化。3 Fluent 讲师来自国家“985工程”重点高校,主要从事物理和数据驱动的力学建模仿真研究,近年来发表SCI论文15篇,授权三项发明专利。擅长领域:计算流体力学、流体力学中的机器学习方法、数据驱动的计算力学、有限元方法等。02培训大纲基于机器学习深度学习驱动的流体力学流场重建技术与应用目录主要内容(*为重点内容)流体数值模拟及 Python编程和数据处理入门一、课程导论1、智能流体力学介绍2、机器学习驱动的流体力学流场重建技术介绍二、流体力学基础1、流体力学基础2、流体力学数据获得方法介绍(实验方法, CFD方法)3、*经典流场模型介绍(绕流,渠道流)4、*流体力学数据分析方法介绍三、OpenFOAM 数值模拟基础1、流体力学求解模型认知(RNAS, LES, DNS)2、OpenFOAM运行环境配制3、*OpenFOAM进行流体计算模拟的基本操作四、 Python 编程入门1、Python编程快速入门2、Python数据科学简介3、*Python流场数据后处理方法(流场云图, 参数统计曲线,概率图,能量谱等)课程实操:(所有示例流场数据和 Python 算法代码提供给学员)1、直接数值模拟(DNS)二维(绕流), 三维流动(渠道流)案例教学2、Python 对流场数据的读取与储存3、Python 绘制流场云图4、Python 计算与绘制常用流场统计结果图 机器学习基础一、深度学习基础1、*深度学习用于计算机视觉: 卷积神经网络(CNN)2、深度学习用于文本和序列: 长短记忆神经网络(LSTM)3、生成式神经网络: 生成式对抗神经网络(GAN)4、*常用的高级深度学习神经网络模型讲解二、强化学习基础1、强化学习快速入门2、高级深度强化案例介绍课程实操:(所有示例流场数据和 Python 算法代码提供给学员)1、在个人电脑上搭建深度学习Python 环境(Tensorflow在CPU以及GPU安装方法)2、使用 CNN 开发第一个深度学习算法(解决分类问题)3、基于深度学习算法开发第一个流场预测算法(数据驱动的方柱绕流流场预测生成)超分辨率问题和二维流场的三维重建问题一、超分辨率问题:1、*研究数据生成方法2、*超分辨率问题的常用的深度学习模型介绍3、*数据后处理方法二、二维流场的三维重建:1、*研究数据生成方法(三维流场)2、*二维流场的三维重建的常用的深度学习模型介绍3、*三维数据后处理方法课程实操:(所有示例流场数据和 Python 算法代码提供给学员)1、完成简单的超分辨率重构问题2、完成简单的三维重构问题3、三维流场数据后处理(结合 Python 代码和 Paraview) 流场去噪问题和流场参数重建问题一、流场去噪问题1、*研究数据生成方法2、*去噪常用的深度学习模型介绍(基于物理约束的强化学习方法以及自监督方法)二、流场参数重建问题1、*流场参数重建问题的常用的深度学习模型介绍2、*PIV 实验数据后处理方法课程实操:(所有示例流场数据和 Python 算法代码提供给学员)1、完成流场去噪问题2、完成流场参数重建问题 基于深度学习的高维插值方法前沿论文分享及SCI论文写作指导目标:以结构-声耦合模型 (Structure-Acoustic Coupling Models)中的传递函数预测为例讲解基于深度学习的高维插值方法一、 高维插值方法1、结构-声耦合模型问题介绍2、*研究数据生成方法(基于Comsol)3、*基于深度学习的高维插值方法二、论文分享以及 SCI 写作指导课程实操:(基于深度学习的高维插值方法应用) 答疑与互动课程期间帮助学员解决理论疑点、技术难点,并可一定程度地给予学员相关研究方向的科研以及论文写作指导基于OpenFOAM和深度学习驱动的流体力学计算与应用目录主要内容经典流体力学与OPENFOAM入门一、经典流体力学核心要点:1、回顾经典流体力学理论,掌握NS方程的基本求解方法和模型2、探索流体力学在工业领域的多元应用3、运用开源软件OpenFOAM进行流体计算模拟的基本操作4、流体力学求解模型认知(RNAS, LES)实操环节:1、OpenFOAM学习:2、掌握OpenFOAM后处理操作3、通过OpenFOAM获取流动信息4、OpenFOAM多种功能使用教程:网络生成,模拟设置5、基于OpenFOAM的矩形柱体LES模拟案例(数据与代码提供给学员)6、OpenFOAM模拟信息的后处理获取流场与压力信息(数据与代码提供给学员)计算流体动力学与人工智能二、机器学习基础与应用核心要点:1、了解Python语言的特征,熟悉常见的机器学习算法2、掌握使用python语言用于数据后处理3、了解计算流体动力学与AI的结合实操环节:1、基于Python语言的CFD数据后处理(数据与代码提供给学员)2、计算流体动力学与AI的结合案例讲解三、时空超分辨率技术核心知识点:1、了解时空超分辨率技术的基本原理与应用2、掌握人工智能技术在湍流时空超分辨率中的创新应用。3、深入理解深度学习与湍流超分辨率的耦合机制。实操环节:1、基于深度学习的流场时序超分辨率处理(数据与代码提供给学员)实验流体力学与人工智能四、实验流体力学核心知识点:1、掌握实验流体力学的基础知识,了解相关实验设备。2、了解机器学习技术在实验流体力学中的应用。3、掌握Python语言进行实验数据的后处理,增强数据处理能力。4、风洞试验实操环节:1、展示机器学习在实验流体力学领域的应用案例。2、运用Python语言处理实验数据(数据与代码提供给学员)五、人工智能与实验流体力学(流场部分)核心知识点:1、掌握实验流体力学数据处理的先进方法2、了解并掌握DNN、CNN等常见深度学习技术在流场重构与预测中的应用。实操环节:1、基于人工智能技术的流场预测与重构方法2、运用DNN技术进行流场预测(数据与代码提供给学员)六、人工智能与实验流体力学(压力部分)1、了解人工智能技术在压力预测预报的使用2、理解融合物理知识的神经网络在误差控制项的使用3、了解融合物理知识的神经网络(PINN)在流体力学的使用实操环节:1、基于人工智能技术的压力预测2、融合物理限制(Physical constrain) 神经网络(PCNN)的时序压力预测(数据与代码提供给学员)空气动力学与人工智能技术七、空气动力学与人工智能技术核心知识点:1、了解爬虫技术在网页数据获取的方法,掌握从开源网站获取信息的技术2、熟悉民航机翼空气动力学性能3、掌握基于MLP的气动性能预测方法实操环节:4、基于爬虫技术的网页数据获取5、基于深度学习的机翼气动性能预测实践案例:基于多层感知机(MLP)的民航超临界机翼气动性能预测(数据与代码提供给学员)深度强化学习学习在流体力学中的应用八、深度强化学习(Deep Reinforcement Learning, DRL)核心知识点:1、掌握深度强化学习的主要框架2、熟悉深度强化学习的常见算法及其应用场景。3、理解深度强化学习中动作空间与观察空间的定义与应用。4、理解代理模型(surrogate model)的概念与意义实操环节:5、构建融合物理限制(Physical constrain)的深度强化学习的环境空间(数据与代码提供给学员)九、深度强化学习的工程实践核心知识点:6、掌握定义离散动作空间/连续动作空间的方法,提升算法设计能力7、学习深度强化学习在工程领域的实际应用,增强解决复杂问题的能力实操环节:8、耦合代理模型的深度强化学习在民航飞机外形优化中的应用:运用深度强化学习进行离散动作空间/连续动作空间的优化(数据与代码提供给学员)课程互动与答疑回顾实践案例课程内容,巩固所学知识、通过答疑加深对知识点的理解与掌握。前沿文献的解读,如SORA技术、风乌技术等,了解人工智能技术在流体力学领域的最新进展,保持学术前沿性。部分案例展示: 基于Fluent和深度学习算法驱动的流体力学计算与应用目录主要内容机器学习与流体力学入门一、流体力学基础理论与编程实战1、流体力学的发展概述2、不可压缩流体力学的基本方程3、偏微分方程数值求解介绍4、傅里叶变换和流体的尺度分析5、伪谱法求解不可压缩流体力学方程案例实践:1、Matlab编程实现有限差分(案例数据与代码提供给学员)2、Python编程伪谱法求解NS方程(案例数据与代码提供给学员)二、Fluent简介与案例实战1、Fluent软件概述:软件功能和特点、Fluent在流体力学中的应用2、网格划分与计算流程:网格划分技术、Fluent计算流程和步骤3、基于Fluent软件的稳态与非稳态流体计算4、基于Fluent软件动网格技术的两相流求解5、Fluent仿真后处理案例实践:圆柱绕流、小球入水的Fluent求解流程(案例文件提供给学员)人工智能深度学习模型与流场超分辨技术三、人工智能基础理论与优化方法1、人工智能的基本概念2、最优化理论算法: a) 最优问题的定义 b) 优化算法介绍3、机器学习算法简介:支持向量机等机器学习算法4、深度学习的基本概念及实战案例实践:Python实现基础网络架构1、梯度下降算法的Python实现(案例数据代码提供给学员)2、二阶函数极值问题的求解(案例数据代码提供给学员)四、深度学习模型在流场超分辨中的应用1、超分辨的基本概念和应用场景2、经典超分辨算法a) 基于局部自适应对偶性先验的最优化方法 b) 超分辨算法的性能评估3、分别基于卷积神经网络(CNN)和生成对抗网络(GAN)的流场超分辨案例与实战案例实践:Python编程实现流场超分辨,不同模型超分辨的优势和劣势分析1、经典模型实现流体超分辨(案例数据代码提供给学员)2、深度学习模型实现流体的超分辨(案例数据代码提供给学员)五、深度学习模型的力学新范式及ODE求解实战1、深度学习模型的动力学解释2、残差神经网络(ResNet)与神经常微分方程(NeuralODE)3、Neural ODE与与流体力学方程求解4、循环神经网络(RNN)与流体动力学时序预测a) RNN的基本概念 b) RNN与隐式算法的对应关系c) 时间序列预测在流体动力学中的应用5、卷积神经网络(CNN)与流场特征提取a) CNN的基本原理 b) 卷积与微分算子的对应关系c) CNN如何用于流场图像分析,如涡识别案例实践:利用Neural ODE求解特定流体(案例数据代码提供给学员)深度学习模型在流体力学中的应用六、神经网络在湍流模拟中的应用1、物理信息神经网络(PINN)2、流动的拉格朗日结构提取与相互作用a) 基于图片的涡旋特征提取b) 基于图神经网络(GNN)的神经网络算法3、嵌入物理信息的神经网络 a) 基于几何对称性改造神经网络 b) 基于拉格朗日结构和几何对称性的神经网络案例实践:Python编程湍流的拉格朗日方法1、流体力学的拉格朗日算法(案例数据代码提供给学员)2、流体力学的拉格朗日神经网络(案例数据代码提供给学员)七、神经网络在空气动力学中的应用1、可压缩流体力学求解的数值方法和机器学习方法a) 可压缩流体力学的数值方法b) 神经网络在激波求解中的应用2、高精度格式在神经网络中的实现a) 高精度格式的主要思想和局限性b) 基于高精度格式的机器学习算法3、深度强化学习(Deep Reinforcement Learning)的理论与算法4、可压缩激波求解案例与编程实战案例实践:Python编程求解可压缩流体力学方程1、高精度格式求解可压缩流体力学方程(案例数据代码提供给学员)2、深度学习模型求解可压缩流体力学方程(案例数据代码提供给学员)流动可视化与新兴技术八、流动生成与后处理1、Tecplot可视化展示标量场、向量场等2、Houdini展示渲染高保真流场3、基于扩散模型(Diffusion Model)的流动生成4、动模态分解及流场预测案例实践:Matlab编程实现DMD(案例数据代码提供给学员)部分案例展示: 03培训特色 1流畅重建专题 1、前沿技术深度聚焦理论与实践结合:结合大量实战案例与项目演练,课程内容涵盖深度学习在流体力学中的最新应用,包括流场重建、超分辨率、三维重建、高维插值方法等。2、全方位技能提升:涵盖专业软件操作(OpenFOAM的流体数值模拟)、编程与数据处理(Python编程和数据处理)、深度学习(神经网络、CNN、LSTM和GAN等)、强化学习、流场数据分析(流场云图绘制、统计结果图计算等)、前沿技术应用、科研论文指导等,全方位提升您的流体力学计算与应用能力。3、个性化科研指导:根据学员的研究方向,提供一定程度的科研指导和论文写作建议。 2OpenFOAM 专题 1.前沿技术深度聚焦:结合大量实战案例与项目演练,聚焦人工智能技术在流体力学领域的最新研究进展。2.全方位技能提升:涵盖经典流体力学、机器学习、深度学习、实验流体力学、计算流体动力学、时空超分辨率、深度强化学习等核心知识,全方位提升您的流体力学计算与应用能力。3.专业优质资源:提供丰富案例数据与代码资源,确保学习效果与实践体验。 3Fluent 专题 1.全方位技能提升:涵盖先进的计算方法(如伪谱法、CNN、GAN、Neural ODE、PINN等)、软件工具应用(Fluent软件、Python编程)、深度学习流场超分辨率、神经网络在湍流模拟中的应用(物理信息神经网络(PINN)和基于图神经网络(GNN))、神经网络在空气动力学中的应用、流动生成与可视化技术,全方位提升您的流体力学计算与应用能力。2.专业优质资源:提供了多个经典案例实践机会,提供Python编程实现和案例数据代码的资源,确保学习效果与实践体验。3.新兴技术探讨:课程还包括了流体力学与深度学习融合的新兴技术,如基于扩散模型的流动生成、动模态分解及流场预测等核心知识的探讨,为学员提供了前沿技术的视野。来源:机电君

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈