首页/文章/ 详情

中北大学李晓峰教授顶刊丨激光粉末床熔融成形Al-Cu-Mg-Si合金热裂行为的数值模拟研究

8月前浏览815

研究背景及目的


 Al-Cu系合金作为最早应用于工业的铸造铝合金,凭借其卓越的强度和耐热性脱颖而出,被广泛应用于国防、航空航天、汽车工业等领域。然而,Al-Cu合金在LPBF成形过程中,粉末床的快速熔化和凝固导致温度场变化剧烈,产生较大的温度梯度;由此产生的残余应力可能导致成形件出现变形、翘曲、开裂等缺陷,严重影响合金的成形质量和力学性能。本文将采用数值模拟的方法分析LPBF成形Al-Cu合金过程中热应力演变规律以及残余应力分布情况,从而为高性能、无裂纹Al-Cu合金的成形提供理论指导。            

           

论文亮点


(1)针对LPBF实际工况,建立了一种考虑粉末-固相转变、材料非线性特性、熔化潜热、熔池对流等因素的三维热-力耦合有限元模型;

(2)采用间接耦合的方式分析了LPBF成形过程中瞬态应力的分布及演化规律;

(3)通过组织观察以及残余应力测试,验证了三维热-力耦合分析模型的可靠性。


           

试验方法


本文借助ANSYS WORKBENCH 数值仿真平台建立了一种多层多轨三维热-力耦合有限元模型,并使用生死单元技术模拟分析了LPBF成形过程中Al-Cu-Mg-Si合金的温度场及应力场的变化规律。


           

结果


图1a为在210W、350mm/s下进行单层激光扫描时熔化层各轨道中心温度的热历史。由图可知,当激光束到达监测点时,温度急剧上升到熔点以上,然后迅速回落,形成一个非常尖锐的峰值形状。随着加工过程的进行,前一扫描轨道对后续未加工区域的预热作用逐渐稳定,温度变化幅度也趋于稳定。各轨道中心的温度变化率如图1b所示,随着轨道数量的增加,最大加热速率先增大后减小,最高达到10.27×10^6℃;最大冷却速率逐渐减小,最低为8.63×10^6℃,但差别不大。

图片            

Fig. 1 Transient thermal behavior of the Al-Cu-Mg-Si processed at P=210 W and V=350 mm/s: (a) Temperature history of P1–P5 and (b) temperature change rates of P1–P5.

图2和图3分别为激光功率和扫描速率对LPBF成形Al-Cu合金冷却后等效应力分布的影响。由图可知,在加工过程结束后,已成形区域开始发生收缩变形,其收缩受到周围低温区材料的限制,使成形区域区产生拉应力。不同成形条件下打印层的残余应力分布趋势基本一致,即打印层上表面的残余拉应力较大,而打印层边缘的残余拉应力较小。但是对于基板来说,除了与打印层界面接触的位置残余应力较大外,整体残余应力较小。当此处的残余应力一旦超过合金的屈服强度时,打印层极易发生翘曲,甚至开裂现象。

图片            

Fig. 2. Equivalent stresses after cooling at different laser powers with a scanning speed of 350 mm/s: (a) Residual stress distribution of different laser powers at room temperature; (b) Maximum equivalent stress variation with laser power; (c) Equivalent stress along the x-axis; (d) Equivalent stress along the y-axis, (e) Equivalent stress along the z-axis.


图片            

Fig. 3. Equivalent stresses after cooling at different scanning speeds with a laser power of 210 W: (a) Residual stress distribution for different scanning speeds at room temperature; (b) Maximum equivalent stress variation with scanning speed; (c) Equivalent stress along the x-axis; (d) Equivalent stress along the y-axis; (e) Equivalent stress along the z-axis.


           

结论


(1)熔池内的热累积导致温度随着扫描轨道数量的增加而升高。

(2)当激光功率从150 W提高到240 W时,熔池的最高温度、熔池寿命和冷却速度都有所提高;当扫描速度从250 mm/s增加到400 mm/s时,冷却速率从6.11×10^6 ℃/s提高至9.86×10^6 ℃/s,而熔池的最高温度和最低寿命均逐渐降低。

(3)LPBF成形后,Al-Cu合金中的残余应力主要表现为拉应力,集中在打印层的中心区域。提高扫描速度、增加激光功率均会提高等效残余应力。


           

前景与应用


本文建立了一种有效、可靠的可用于描述LPBF成形Al-Cu合金的三维热-力耦合有限元模型,通过该模型可快速分析成形过程中热应力演变规律以及残余应力分布情况,从而获得高性能、无裂纹Al-Cu合金的成形工艺窗口。

关注我们, 万物皆可3D打印

Yuling Niu, Li Zhang, Jianhong Wang, Jinfang Zhang, Xiaohui Yang, Bin Liu, Yuankui Cao, Hong Xu, Ruibin Duan, Xiaofeng Li. Cracking mechanis m of Al-Cu-Mg-Si alloy during laser powder bed fusion by numerical simulation. Additive Manufacturing Frontiers, Volume 4, Issue 3, 2025, 200213.

https://doi.org/10.1016/j.amf.2025.200213


来源:增材制造硕博联盟
ACTWorkbenchAdditive非线性航空航天汽车UM铸造裂纹理论材料单元技术试验ANSYS
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-12-11
最近编辑:8月前
增材制造博硕联盟
硕士 聚焦增材制造科研与工程应用,致...
获赞 129粉丝 96文章 625课程 0
点赞
收藏
作者推荐

华中科大顶刊《JMRT》丨激光增材制造新型TPMS结构:高换热性能和良好的能量吸收性能

航空航天需要轻量化、高强度、高效换热的组件。华中科技大学史玉生团队提出了一种结合材料成分和宏观结构设计的优化方法。通过激光粉末床融合(LPBF)技术制造了由氧化物弥散强化(ODS)CuCrZr合金组成的三周期最小表面(TPMS)点阵结构。 CuCrZr合金粉末本研究中使用的CuCrZr合金粉末是辽宁冠达增材制造有限公司通过气雾化工艺生产。粉末颗粒主要为球形或近球形,表面相对光滑。使用激光粒度分析仪分析了粒径分布,结果显示粒径范围为20.3 μm至55.1 μm,平均粒径为33.9 μm。这些特性——适中的粒径和高球形度——表明CuCrZr粉末非常适合用于激光粉末床融合(L-PBF)工艺。CuCrZr合金粉末的化学成分ElementCuCrZrOCNSContent (wt.%)Bal.0.990.0930.0350.0063<0.001<0.001本研究制备了四种具有不同Y₂O₃掺杂水平的CuCrZr-Y₂O₃复合材料:CuCrZr-0.5 wt.% Y₂O₃、CuCrZr-1.0 wt.% Y₂O₃、CuCrZr-1.5 wt.% Y₂O₃和CuCrZr-2.0 wt.% Y₂O₃。Y₂O₃,纯度≥99.99%,平均粒径为50 nm。CuCrZr合金粉末和Y₂O₃纳米粉末混合30小时。 TPMS点阵结构样品的尺寸设置为24 mm × 24 mm × 24 mm,由3×3×3个单元晶格组成,每个单元晶格尺寸为8.0 mm。L-PBF工艺参数LaserPower(W)ScanningVelocity(mm/s)LayerThickness(μm)HatchSpacing(mm)HatchAngle(°)VED(J/mm3)375500500.1267°125压缩有限元模拟CompositesDensity(g/cm3)Young’sModulus(GPa)YieldStrength(MPa)SpecificHeatCapacity(J/(kg·K))ThermalConductivity(W/(m·K))CuCrZr-8.920111.023200.8340.344114.258CuCrZr-0.5 wt.% Y2O38.901124.440214.2090.298123.519CuCrZr-1.0 wt.% Y2O38.881112.494174.2930.391102.736CuCrZr-1.5 wt.% Y2O38.861103.398111.4040.28987.296CuCrZr-2.0 wt.% Y2O38.94099.852105.1190.33074.163换热有限元模拟与试验ParameterValueSpecific heat of air at RM, Cpf1004 J/(kg·K)Thermal conductivity of air at RM, kf0.02607 W/(mK)Dynamic viscosity of air at RM, μf1.8477·10-5 Pa·sDensity of air at RM, ρf1.1305 kg/m3Applied heat flux, q70000 W/m2Inlet air temperature, Tin298 KReynolds number, Re1445.5-7227.5 材料成分设计如光学显微镜图像所示,所有样品均在不同程度上表现出微裂纹和孔隙。复合材料的成形质量随着Y₂O₃含量的增加先改善后恶化。 SEM分析显示,当Y₂O₃含量在1.0到2.0 wt.%范围内时,在高倍率图像中观察到显著的Y₂O₃纳米颗粒团聚。特别地,在掺杂量为2.0 wt.%的样品中,即使在低倍率SEM图像中也可见大量团聚的Y₂O₃,表明颗粒严重聚集。使用电子背散射衍射(EBSD)研究了通过L-PBF制造的CuCrZr和CuCrZr-0.5 wt.% Y₂O₃样品的晶粒取向和微应变分布。图(a)和(b)分别展示了两个样品垂直截面的反极图(IPF)。结果表明,CuCrZr样品在<101>//构建方向(BD)上表现出优选的晶体学织构,这可以归因于固化过程中定向热流对晶粒生长的影响。相比之下,CuCrZr-0.5 wt.% Y₂O₃样品显示出随机的晶粒取向,没有明显的织构。CuCrZr和CuCrZr-0.5 wt.% Y₂O₃样品的晶界(GB)图可以观察到,添加适量的Y₂O₃诱导形成了更多的LAGBs。LAGBs数量的增加可以通过位错阻碍和晶界强化等机制增强基体的强度。CuCrZr样品的平均晶粒尺寸为20.99 μm,而CuCrZr-0.5 wt.% Y₂O₃样品的平均晶粒尺寸减小到17.74 μm。Y₂O₃颗粒在低掺杂水平下通过提高成形质量来增强材料强度,而过量掺杂超过临界阈值则会导致颗粒团聚,显著降低机械性能。 CuCrZr和CuCrZr-0.5 wt% Y₂O₃样品的断口表面表现出韧性特征,表明韧性断裂。相比之下,CuCrZr-1.0 wt% Y₂O₃样品的横截面显示出类似“山川”的形态(用黄色虚线轮廓),这是脆性断裂的特征。 所有样品的热导率均随温度表现出明显的正相关关系。值得注意的是,随着Y₂O₃掺杂浓度的增加,热导率先提高后降低,CuCrZr-0.5 wt% Y₂O₃样品表现出最高的热导率。这种增强归因于最佳Y₂O₃含量能够改善复合材料的整体质量,减少冶金缺陷,从而促进更有效的热传导。相反,过量的Y₂O₃添加由于成形质量下降导致热导率降低。综合考虑机械和热性能,选择CuCrZr-0.5 wt% Y₂O₃复合材料用于后续TPMS点阵结构的制造。 点阵结构的成形质量 每个样品均表现出完整的表面结构,未观察到明显的缺陷或裂纹。I-WP点阵结构在正视图中表现出最严重的粉末粘附。换热性能努塞尔数(Nu)和摩擦系数(f)是评估点阵结构换热性能的关键无量纲参数,它们共同影响换热效率和能耗之间的平衡。Nu反映了流体-点阵界面处的对流换热强度,其中较高的Nu表示更好的换热效果。TPMS点阵结构,由于其复杂的孔隙几何形状,增强了湍流混合和边界层破坏,与传统散热结构相比具有更高的Nu。f衡量点阵的压降,较大的f表示更大的能量损失和流动阻力。TPMS点阵结构的高比表面积和曲折的流道增加了f,在流体流动过程中加剧了压降。 Nu随雷诺数显著增加,而f减小。这种行为的发生是因为较高的雷诺数对应较高的气流速度,这破坏了点阵表面上的热边界层,增强了换热效率。同时,速度的增加改善了边界层与外部流体之间的动量传递,减少了由流动分离引起的压力损失,从而降低了摩擦系数。总体而言,在此雷诺数范围内,基于Nu的换热性能排序为Diamond > Gyroid ≈ I-WP > Primitive。 点阵结构的比表面积顺序为:Diamond > I-WP > Gyroid > Primitive,TPMS点阵结构的换热性能与其比表面积密切相关。然而,本研究中的CFD模拟结果表明,换热性能不仅受比表面积影响,还受点阵结构内部流道特征的影响。 速度分布表明,Gyroid和Diamond点阵结构流道和孔隙内的空气速度显著高于入口和出口区域,且分布更为均匀。相比之下,Primitive和I-WP点阵结构内部流道中的空气速度极高,甚至超过了扩展入口、出口以及Gyroid和Diamond结构中的速度。然而,Primitive和I-WP点阵结构孔隙内的空气速度非常低,甚至低于扩展出口的速度。这些流动特性表明,Gyroid和Diamond点阵结构中相对均匀的气流确保了空气与换热表面之间的充分接触,增强了气流扰动并优化了热质传递,从而提高了对流换热性能。相反,Primitive和I-WP点阵结构中高度不均匀的速度分布导致局部高速区域快速移除热量,但低速区域热量积累,限制了整体换热效率。这解释了尽管I-WP点阵结构的比表面积大于Gyroid,但Gyroid表现出略高的努塞尔数。 压缩性能和能量吸收由弹性模量和峰值平台应力评估的压缩机械性能排序如下:Diamond > Gyroid > I-WP > Primitive。这些差异主要源于每种点阵结构特征所固有的不同变形和失效模式。 Gyroid点阵在整个压缩过程中倾斜表面角度显著减小(用白色虚线突出显示),样品表现出弯曲和屈曲变形特征。在10%应变时,Diamond点阵左下角板首先弯曲(白色虚线圆圈),随后相邻板逐渐倾斜;到20%时,左上角出现阶梯状失效,并伴有45°斜向剪切失效模式(白色虚线)。对于Primitive和I-WP点阵,在10%应变时圆形孔开始椭圆形化,到20%时中央层优先屈服(白色虚线框),且椭圆形变形越来越严重。在30%应变时,所有点阵均表现出不同程度的侧向“鼓胀”(黄色虚线),这是由于压头和底座端部附近变形受限,而中部区域向侧方扩展所致。Primitive和I-WP点阵显示出加剧的中央层失效(白色虚线框),而上层和下层与20%阶段相比相对稳定。在Diamond点阵中,左上角阶梯状失效进一步发展,同时右上角出现结构坍塌(白色箭头),表明即将发生连续层失效。随着压缩继续进行,点阵表面相互接触,表明即将致密化。在整个测试过程中,未在任何TPMS点阵中观察到明显断裂,表明其优异的韧性,并突出了它们在实际应用中作为有效能量吸收和缓冲组件的潜力。 在相同应变条件下,应力集中区的大小排序为Diamond > Gyroid > I-WP > Primitive。由于较大的应力集中区域与更大的承载能力和更高的应力水平相关,这些发现进一步证实了准静态压缩模拟的准确性。 在弹性变形阶段,所有四种晶格结构的能量吸收值都呈线性增加,并保持可比性。在进入应力平台和连续失效阶段后,能量吸收排名遵循Diamond>Gyroid>I-WP>Primitive。总体而言,能量吸收曲线的平滑性没有明显的波动,这表明晶格结构具有良好的塑性和稳定的能量吸收性能。结论1.随着Y2O3掺杂浓度的增加,复合材料的机械强度和导热系数先增加后降低。CuCrZr-0.5 wt.%Y2O3复合材料表现出最佳的整体性能,室温抗拉强度为309.92 MPa,导热系数为123.519 W/(m·K)。微观分析表明,机械性能的提高主要是由于四种强化机制(如沉淀强化)的协同作用和成形质量的提高,而导热性的提高主要归因于更好的材料成形性所实现的热传导途径的优化。2.四种TPMS晶格结构的体积分数偏差顺序为:Diamond>I-WP> Gyroid> Primitive。这种偏差主要归因于与更高比表面积相关的粉末粘附力增加。关于热交换性能,Diamond在相同的泵送流量下表现出最佳效果,这是由于其较大的比表面积和曲折的内部流动通道引起的湍流增强。在机械和能量吸收性能方面,Diamond也优于其他结构。这种优势主要是由于其剪切主导的破坏模式,与Primitive等结构中观察到的逐层破坏模式相比,它在承载能力和能量吸收方面具有显著优势。 来源:增材制造硕博联盟

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