下载文档原格式
基于生死单元的激光熔覆温度场数值模拟摘要:计算了不同激光功率条件下粉末粒子到达基体前的温升情况,并将粉末粒子到达基体前的温度作为初始条件,采用生死单元法对单道和多道激光熔覆温度场进行了研究.利用熔池尺寸和形貌,验证了模型的可靠性.结果表明,粉末粒子温升和激光功率呈线性关系,单道熔覆层的温度变化呈一个锯齿状,升温过程近似呈直线上升,降温曲线近似呈双曲线的一支,而多道熔覆过程中,温度场呈后拖的偏椭圆状.节点上的热循环经过逐渐增大的峰值,峰值温度最终趋于稳态.0 序言激光熔覆按送粉工艺的不同可分为两类.即粉末预置法和同步送粉法.同步送粉法具有易实现自动化控制,激光能量吸收率高,无内部气孔,尤其熔覆金属陶瓷,可以显著提高熔覆层的抗开裂性能,使硬质陶瓷相可以在熔覆层内均匀分布等优点,具有广阔的应用空间.国内学者运用ANSYS对激光熔覆过程的温度场和应力场已经做了大量的研究工作[1-4].目前通过ANSYS 模拟激光熔覆温度场的研究并没有考虑激光束与粉末的交互作用,实际上激光束最先和粉末作用,激光束除了损失的能量以外,一部分被熔覆粉末吸收,另一部分则透过粉末被基体吸收.基体除了直接吸收激光束能量以外,还吸收粉末传递给基体的能量.因此很有必要在模拟之前搞清楚激光能量的分配情况,这样建立起来的模型与实际才更接近,模拟结果更有说服力.文中将粉末到达基体前吸收能量后的温升作为初始温度场加载给基体,同时采用ANSYS中生死单元技术模拟熔覆单元的生长过程,基体吸收的能量通过高斯体热源加载,模拟送粉激光熔覆的温度场分布情况.在此基础之上,模拟了多道次激光熔覆温度场,对多道激光熔覆温度场进行了研究.1 粉末到达基体前的温度粉末在遮挡激光的同时,也吸收了部分激光能量,从而使自身的温度升高.实际上,在不考虑等离子体影响(能量密度低于105W /cm2)情况下,粒子直接吸收激光辐射能,并放出辐射能.在空气中粉末颗粒也会由于空气对流散失能量,粒子之间也会相互加热这些能量在总能量中的比例很小,目前关于粉末颗粒温升的模型并不多见,而且在模型中需要建立假设条件[5].为了计算方便,在模型中假设:(1)粉末颗粒在气)粉射流中的体积分数很低,可以忽略激光的反射、折射、颗粒离子之间的相互加热和光束遮蔽等影响. (2)粉末颗粒是半径为rP的球体.由于粉末颗粒足够小,在能量计算时将其看成一个点,粒子的热导率为无限大,即认为粉末颗粒的温度是均匀一致的,在迎光面和背光面没有差异.(3)粉末颗粒只在迎光面吸收能量,但对外辐射则在整个球体表面发生. (4)粉末不吸收来自基体的反光.由以上假设,根据粒子的能量方程可以求出粉末粒子的温升.方程是非线性方程,使用Matlab软件采用迭代法求解.当激光功率P=2 kW时,在1 500~1 600 K范围内方程有解,于是初始值设为T=1 500 K,通过迭代求出方程的一实根为T=1 570 K.改变激光功率,得到不同激光功率下粉末离子到达基体时的温度,如图1所示,随着激光功率的增大,粉末粒子达到基体前的温升逐渐增高,且成近似线性关系增长.当激光功率P=1 500W时,计算出粉末粒子温升T=1 267 K,而Ni60粉末的熔点约为1 300 K.说明激光功率低于1 500W 时,粉末达到基体前不会熔化,这样粉末粒子就会以固体颗粒的形式和基体碰撞,飞溅严重,即使是部分粉末颗粒熔化形成熔池,基体熔深也得不到保障;当激光功率P=3. 0 kW时,T=2 111 K,这时粉末粒子的温度远大于熔点温度,会造成部分粉末颗粒的烧损.因此,在实际生产中应控制激光功率在一个合理的范围内图1 激光功率与粉末粒子温度关系曲线Fig11 Relation curve between laser power and tempera-ture rises of powder particles2 单道激光熔覆的模型在建立热源模型的过程中,粉末到达基体前与激光束发生交互作用,激光束对粉末的作用通过初始温度vT实现,即假设粉末飞行过程中吸收的有效能量全部用于升温,作为初始温度场施加给熔覆单元.根据Picasso的理论,基体吸收的能量可分为直接吸收和粉末热传递两部分,基体直接吸收的热量通过体热源形式实现.采用APDL语言编制热源的移动过程程序,实际操作中,始终以全局笛卡儿坐标系为求解坐标系,而载荷在局部坐标系下施加,这样通过熔覆方向上的坐标变换,实现全局坐标向局部坐标的转化,即z=Z-vt (1)式中:Z为全局笛卡尔坐标系下的坐标; z为局部笛卡尔坐标系下的坐标; v为光源移动速度; t为光源移动时间.采用ANSYS中生死单元技术模拟熔覆单元的生长过程,在计算开始时刻,所有熔覆层单元均设定为/死0单元.在随后每一步的计算中,首先判断所有/死0单元是否落入激光束的照射区域,如落入激光束照射区域,就将其激活,纳入计算模型中.2.1 热源模型文中粉末粒子的温升采用解析计算;基体直接吸收的激光束热源模型采用柱状高斯体热源,以此模拟熔覆层和基体中热量的三维分布和传导.体热密度表示为Q=Qmexp(-3r2/r2a)exp(-BZ) (2)r=(vt-x)2+y2(3)Qm=CPPabh(4)式中:Qm为加热斑点中心的最大体热流;B为激光体热密度沿厚度方向的衰减系数; r为某一深度任意一点(x,y)距体热流中心的距离; ra为激光的有效作用半径;a,b为高斯热源的长短轴;h为熔覆层深度;C为基体对激光的吸收系数.2.2 材料热物理性能参数激光熔覆过程就是加热)熔化)凝固)冷却的过程,其中包含着相变.相变热模型在数学上是一个强非线性问题,使计算发生困难.对于Q235基体材料,文中采用显热熔法,在糊状的熔化带内调整比热容c=$H/$T来近似计算[6],其中$H为熔化潜热,$T为熔化温度区间.粉末颗粒熔化潜热$HP=0. 26 J/mg,比热容c=4. 59@108J/kge,熔化温度区间为[1 230 K, 1 310 K].基体材料熔化潜热$HW=0. 273 9 J/mg. Ni60自熔合金的密度取常温下的值8. 378@103kg/m3.为了得到好的收敛解,激活牛顿)拉普森方法的线性搜索.3 单道激光熔覆温度场的模拟考虑到对称性,取工件的1/2进行分析.采用八节点六面体等参单元对工件进行离散.为保证计算精度,在熔覆层及其相邻部位对网格进行细化.采用生死单元技术编制程序,首先将熔覆层单元存入预先定义的数组中,将其全部/杀死0,然后建立局部坐标系,通过*DO循环实现光斑的移动,判断熔覆层单元是否落入激光光斑范围内,如果有则激活单元,如果没有则直接进入下一次循环.将粉末颗粒的温升作为初始载荷加载给被激活的熔覆单元.基体直接吸收的能量通过高斯体热源施加,用ANSYS自带的函数编辑器编写高斯函数.送粉速率va=71. 36 mg/s.模拟实现了熔覆层随着激光光斑的移动而逐渐生长的过程,图2为基体温度降温时任一温度场等值面图.图中最高温度区域边界,即基体对称面上距离基体上表面0. 50mm温度为1 653 K,而基体熔点温度为1 670 K,基体熔深即为0. 50 mm,这与实际检测结果0. 48 mm第5期赵洪运,等:基于生死单元的激光熔覆温度场数值模拟83 吻合较好,通过这种寻找熔点的方法确定熔宽为216mm,与实测结果2. 8 mm吻合较好.试验试样熔覆层金相组织形貌和数值模拟形貌对比如图3所示.图2基体温度场等值面图Fig. 2 Substrate iso-surface of temperature field图3试样金相组织形貌和数值模拟形貌对比Fig. 3 Metallography film and temperature field of moltenpool 根据以上分析,文中采用的模型是合理可行的,可以进行下一步的模拟分析.图4为熔覆层中序号为740的节点(x=0.000,y=0. 000 8 m,z=0. 021 m)的温度变化曲线.当时间t=6. 667 s时,该点开始受到激光光斑直接照射,温度迅速升高,当t=7. 0 s时激光光斑中心移至该点,该点温度升至最高点,当t=7. 333 s时,光斑移出该点,该点温度迅速降低.升温曲线近似呈直线,降温曲线近似呈双曲线的一支,整个曲线呈锯齿状.对曲线图求导数,就可以得到该点温度变化速率曲线.图4 温度变化曲线Fig14 Variation of temperature4 多道激光熔覆温度场模拟的实现4.1 多道搭接激光熔覆模型激光功率选取2. 0 kW,扫描速率为3. 0 mm/s,送粉速率取71. 36 mg/s,根据前文熔覆粉末到达基体时的温度的计算,初始温度为1 570 K.为了节省计算时间,只熔覆3道,搭接率取20%.这是因为多道熔覆数据文件较多,输出载荷步数超过1 000步.为了分析熔覆层上的热循环,在第一道顶点上取点A,在第二道和第三道顶点上取点B和点C,对不同熔覆层上的相同位置的点的热循环进行分析.图5为多道搭接取点示意图.图5 多道搭接顺序及取点示意图Fig15 O rder ofm ult-ipass laser cladding and taken pionts4.2 多道搭接激光熔覆模拟结果分析熔覆一道需要16. 667 s,激光器回程时间设置为1 s,即道间搭接时间间隔为1 s,熔覆完成需约52s,在空气中冷却.图6为激光器第一道回程过程恰好结束,即第二道结束时的温度场分布,虽然回程时间很短,但温度场却发生了很大变化,最高温度为861 K与最低温度630 K接近,温度场趋于平衡.图7为第三道熔覆过程中某时刻温度场分布,可见温度场呈后拖的偏椭圆状,即温度场椭圆不以光斑图6 第二道恰好结束时多道熔覆温度场云图Fig16 Contour of temperature field as second clad justfinished中心对称,而是偏向已形成熔覆层的一侧,这是多道温度场不同于单道熔覆的地方.多道熔覆之所以会形成偏椭圆形状,是因为前面先形成熔覆层的道次对后续熔覆产生了影响,先熔覆的熔覆层相对于未熔覆的区域来说有一个初始温差,使得温度场不能沿扫描线对称.图7 熔覆第三道时温度场分布云图Fig17 Contour of temperature field as cladding third clad按照图5所示的取点方法,分析各熔覆层上的热循环. A,B, C三点的温度随时间变化如图8所示.由图8知,三条曲线上最高峰值温度逐渐增高,这是因为A点是第一道熔覆层上的点, B, C 分别是第二道和第三道熔覆层上的点,前面道次的熔覆相当于给后续道次的熔覆起到预热作用.图8 A与B及C点温度变化曲线Fig18 Variation of temperature at points A, B and C5 结论(1)通过ANSYS生死单元技术实现了激光熔覆送粉过程的模拟计算.并通过试验证明该模型用于模拟送粉激光熔覆温度场是合理的,发现粒子温升和激光功率之间呈线性关系.(2)模拟发现单道熔覆过程熔覆层的温度变化呈一个锯齿状,升温过程近似呈直线上升,降温曲线近似呈双曲线的一支.(3)多道熔覆过程中温度场呈后拖的偏椭圆状,即温度场椭圆偏向已形成熔覆层的一侧.无论是基体还是熔覆层上的点,后一道上的最高温度均高于前一道次上的温度.但随着熔覆道次的增多,最高温度会趋于稳定.。
同轴送粉激光熔覆的粉末熔化过程建模与仿真
朱明;严瀚林;张宗智;杨骞;石玗;樊丁
【期刊名称】《焊接学报》
【年(卷),期】2024(45)5
【摘要】同轴送粉激光熔覆中,粉末与激光的交互作用会直接影响熔覆成形的精度与质量,红外摄像无法直接获得激光内粉末的熔化行为,因此通过分析粉末对热量的吸收情况,使用高速摄像系统采集粉末熔化的动态行为,并建立粉末熔化过程的动态解析模型,通过仿真分析激光功率对不同熔化阶段的影响,及粉末进入熔池时的温度特征.结果表明,高速摄像系统采集粉末在激光中的动态熔化行为中存在“固态→固液两相态→液态”3个典型熔化特征阶段,粉末熔化动态行为可以使用数学解析模型解析,而不同阶段的热物理行为具有粉末与激光热交互作用的动态解析模型,分析了激光功率、离焦量、载粉气流量对粉末熔化行为的影响,同时通过仿真分析不同激光功率对各个特征阶段持续时间的影响,预测粉末颗粒到达基材的温度分布,发现激光功率从100 W增加至1 500 W时,粉末进入熔池的温度呈非线性变化,温度从750℃增加至3 250℃.
【总页数】10页(P47-55)
【作者】朱明;严瀚林;张宗智;杨骞;石玗;樊丁
【作者单位】兰州理工大学;兰州职业技术学院
【正文语种】中文
【中图分类】TG456.7
【相关文献】
1.同轴送粉激光熔覆过程中粉光匹配影响因素及控制研究
2.基于UG的激光熔覆送粉器实体建模与运动仿真
3.激光熔覆中同轴送粉气体-粉末流数值模拟
4.送粉激光熔覆中送粉速率对激光束与粉末流相互作用的影响(英文)
5.同轴激光熔覆粉末的熔化行为表征与分析
因版权原因，仅展示原文概要，查看原文内容请购买。
同轴载气送粉激光熔覆粉末流参数研究同轴载气送粉激光熔覆技术是一种先进的制备方法，具有高效率、高质量、高可控性等优点，因此在加工领域得到了广泛的应用。
在使用该技术进行熔覆粉末流成形时，粉末流参数的选择对于熔覆质量的影响非常大。
本文将针对同轴载气送粉激光熔覆粉末流参数进行研究，从以下几个方面进行分析。
一、粉末流速度粉末流速度是影响熔覆质量的重要参数之一。
当粉末流速度过大时，会导致熔滴的飞溅和熔覆质量下降；当粉末流速度过小时，会导致熔池面积减小、熔覆质量下降。
因此，应根据熔覆材料的性质和设备的能力，选择合适的粉末流速度。
二、载气流量载气流量是指粉末颗粒在流动过程中所受到的空气流动的力量大小。
合适的载气流量可以保证熔覆粉末流形成正常，同时也可以将存在于熔合池中的杂质、气泡等物质排出，保证熔覆质量的提高。
过大或过小的载气流量都会影响熔覆质量，应选取合适的值。
三、喷嘴与底板的距离喷嘴和底板的距离也是影响熔覆质量的重要参数之一。
当喷嘴和底板的距离过大时，粉末流速度将变低，影响熔覆质量；当喷嘴和底板的距离过小时，容易导致熔滴飞溅，同样影响熔覆质量。
因此，喷嘴和底板的距离应根据熔覆材料的性质、熔覆质量要求和设备的能力，进行合理的选择。
四、粉末流角度和喷嘴半径粉末流角度和喷嘴半径也会影响熔覆质量。
当粉末流角度过大或喷嘴半径过小时，会导致粉末流速度过大，熔滴飞溅，熔覆质量下降；当粉末流角度过小或喷嘴半径过大时，会导致粉末流速度过小，熔覆层减薄，熔覆质量下降。
因此，粉末流角度和喷嘴半径也应进行合理的选择。
通过以上分析，我们可以得出同轴载气送粉激光熔覆的粉末流参数的研究，是保证熔覆质量和生产效率的必要措施。
因此，在使用该技术进行熔覆粉末流成形时，应根据熔覆材料的性质、熔覆质量要求和设备的能力，选择合适的粉末流速度、载气流量、喷嘴和底板的距离、粉末流角度和喷嘴半径。
光内送粉喷头空间变姿态粉气流场数值模拟姜付兵;石世宏;傅戈雁【摘要】光内送粉喷头是一种新型的激光熔覆同轴送粉喷头,具有空心光束、单根粉束指向性好、无干涉、光粉气一体同轴等优势.为了分析喷头倾斜角度和外层准直气流速对粉束汇聚特性的影响,采用FLUENT软件进行模拟计算,对喷头在空间变姿态作业过程中单喷嘴出口粉束的汇聚特性进行了理论分析和实验验证.结果表明,喷嘴出口粉末质量浓度分布沿x方向服从高斯分布,对称轴上沿-y方向质量浓度先基本保持恒定,后逐渐下降为0kg/m3;随喷头倾斜角度的增大,粉末汇聚长度h不断减小,h变化范围为11mm～17mm,偏移量△x大致呈增大趋势,△x变化范围为0.0mm ～ 0.4mm;随外层准直气速度增加,汇聚长度h呈轻微上升趋势,h上升幅度为1mm～2mm,而偏移量△x在不同倾斜角度范围内呈现出不同变化趋势,△x变化幅度为0.05mm～0.15mm;实验数据与模拟结果变化趋势基本吻合.所建模型和模拟结果对于选择和优化实验工艺参量具有参考价值.【期刊名称】《激光技术》【年(卷),期】2015(039)002【总页数】7页(P145-151)【关键词】激光技术;气固两相流;数值模拟;光内同轴送粉【作者】姜付兵;石世宏;傅戈雁【作者单位】苏州大学机电工程学院,苏州215021;苏州大学机电工程学院,苏州215021;苏州大学机电工程学院,苏州215021【正文语种】中文【中图分类】TG156.99;TN249Key words:laser technique; gas-solid two-phase flow; numerical simulation; inside-beam coaxial powder feeding*通讯联系人。
E-mail：*******************.cn金属零构件空间修复和3维自由成形在工程与国防领域具有很好的应用前景[1]，而光粉耦合工艺是实现激光直接成形的关键技术之一，对送粉喷头喷嘴出口流场特性进行数值模拟分析，对于优化喷嘴结构和工艺参量具有重要意义。
第29卷 第11期光 学 学 报V ol.29,No.112009年11月ACTA OPTICA SINICANovember,2009文章编号:0253-2239(2009)11-3114-07激光熔覆中同轴粉末流温度场的数值模拟杨洗陈 栗 丽 张 烨(天津工业大学激光技术研究所,天津300160)摘要 在激光同轴送粉熔覆中,由于激光与粉末流相互作用,粉末流整体温度分布直接影响激光熔覆的质量。
基于非预混燃烧模型,将激光相处理为连续性介质,粉末颗粒相看作离散相物质,建立了激光作用下粉末流的质量、动量和能量方程。
用Fluent 软件进行了不同激光功率和粉末流速度条件下粉末流整体温度场数值模拟,讨论了各种参数对温度场分布的影响。
为了验证该模型的准确度,利用CCD 比色测温方法测量了粉末流整体温度场分布。
结果表明,数值模拟与CCD 检测结果具有良好的一致性,数值模拟结果对激光熔覆具有指导意义。
关键词 激光熔覆;粉末流;温度场;数值模拟;Fluent 软件;CCD 相机;温度测量中图分类号 T N 249 文献标识码 A doi :10.3788/AOS 20092911.3114Nume rical Simulation of Temperature Fie ld of Coaxial Powder Flowin Lase r CladdingYang Xichen Li Li Zhang Ye(La ser Pr ocessing Cent er ,T ian jin P olytechnic Un iver sit y ,T ia n jin 300160,Chin a )Abstract I n la ser coaxia l cladding,laser casing qualities are directly affected by temperature field of powder flow for the interaction between laser and powers.According to the model of non -premixed c om bustion,regarding laser beam as c ontinuity medium pha se and powder a s disperse medium phase,the conservation equations of mass,m oment um and energy are established in laser and powder puter simulations of tem perature field in different para meters are finished using Fluent software.Some effects of laser c ladding parameters on temperat ure distribution such as laser power and powder flow velocity are discussed.Temperature field distribution in the powder flow is measured by CCD c amera.It is shown that simulation and experimental results agree well,numeric al simulation of temperature field in powder flow is important for laser cladding.Key wo rds la ser cladding;powder flow;temperature field;num eric al simulation;Fluent software;CCD camera;m ea surments temperature收稿日期:2009-07-20;收到修改稿日期:2009-09-28基金项目:国家自然科学基金(60478004)和天津市科技支撑计划重点项目(08ZCKF GX02300)资助课题。
送粉式激光熔覆裂纹机理的数值模拟蒙宣伊【摘要】在对送粉式激光熔覆应力产生机理进行深入分析的基础上,结合应力场数值模拟分析,结果表明熔覆层中心偏基体一侧是裂纹敏感区,该区域处于纵向拉伸、横向和厚向压缩的第三类主应变状态.【期刊名称】《湖南工程学院学报（自然科学版）》【年(卷),期】2010(020)001【总页数】4页(P39-42)【关键词】送粉式激光熔覆;数值模拟;应力场【作者】蒙宣伊【作者单位】湘电集团,风能有限公司,湘潭,411102【正文语种】中文【中图分类】TG404激光熔覆是一种新型的材料加工与表面改性技术,它可以在低成本材料上制成高性能表面,提高材料的综合性能,降低能源消耗,该技术引起了西方各国的高度关注[1-5]. Yunchang Fu,A.Loredo等建立了熔覆过程的解析模型[6],示意图如图1所示.图1 激光熔覆模型示意图目前已有模拟温度场和应力场的文献,但均没有考虑激光束与粉末的交互作用,因此,研究激光熔覆过程的应力产生机理,从而提出减少和消除残余应力的方法,对于激光熔覆技术的完善和推广具有重要实际意义.1 熔覆及凝固的有限元模型基体Q235力学性能参数见表1,熔覆层材料力学性能参数如表2所示.表1 基体Q235力学性能参数温度T/℃ 20 250 500 750 1000 1500 1700 2500弹性模量E/105 MPa 2.05 1.87 1.50 0.7 0.2 0.19 0.12 0.12屈服应力σs/MPa 220 - 130 40 25 - 1.0 0.1切向模量 Et/102 MPa 205 - 150 70 20 - 18 12泊松比μ 0.28 0.29 0.31 0.35 0.40 0.45 0.50 0.50线性膨胀系数at/10-6.0C 11.0 12.2 13.9 14.8 13.4 13.3 13.1 13.1表2 熔覆层材料力学性能参数温度T/℃ 20 100 200 300 400 600 1000 1200 2500弹性模量E/105 MPa 2.06- - - - - - - -切向模量Et/102 MPa 355 - 305 - - 299 268 245 230泊松比μ 0.25 - - - 0.28 - - 0.3 0.3线性膨胀系数at/10-6℃ 11.5 11.6 11.8 12.2 12.6 12.9 - 14.0 14.5本文中环境温度为293 K,工件无初始应力,采用对称模型,位移和对称约束见图2,这样让约束点尽量分散,且远离光斑加热区,避免计算产生偏差.本文在模拟激光熔覆的热应力场时,采用生死单元法实现熔覆层的生长过程,在熔覆开始前把熔覆单元“杀死”,并在每一步热应力计算时,将对应温度场的计算结果加载到应力场模型中.同时采用改变单元属性法,在开始计算前,先定义熔融金属的材料属性(熔融金属包括处于熔化状态的熔覆粉末和基体熔池金属)材料编号为3,使它们的屈服极限和弹性模量都很低,且不随温度变化,熔覆材料属性编号为2,基体材料编号为1.熔覆材料和基体的力学性能参数都随温度变化.图2 位移与对称约束在程序计算过程中,首先改变结构单元为热单元,在温度场后处理器中读取节点温度值,将超过熔点的熔覆层单元和基体单元分别存入温度表中.再转换热单元为结构单元,并改变这些单元材料属性.2 应力场模拟结果与讨论2.1 熔覆层中心应力应变分布图3 应力随时间变化曲线图图3 为熔覆层中心某点应力随熔覆时间变化曲线图,图中SZ表示纵向应力,SX表示横向应力,SY表示厚向应力.当时间t=6.667 s时,光斑照射到该点时,该点所在的单元被激活,并迅速达到熔点温度,应力场计算近似等于零;t=7.333 s,光斑移出该点,该点很快降至熔点温度,并开始凝固.材料冷却收缩,受到基材或自身温度不均匀的限制,处于拉应力状态.在塑性范围内,拉应力纵向最大,横向次之,厚度方向最小.拉应力的数值取决于材料的流动应力,随温度的下降,流动应力上升,拉应力表现为上升趋势,当材料进入弹性状态后,横向和厚度方向上的拉应力下降明显,纵向应力基本保持不变.图4 残余应力分布云图图3 (a)为冷却600 s后的纵向残余应力分布云图,由图看出,最大纵向残余应力σzmax=978 MPa,且靠近基体的熔覆层上,这是由于熔覆层材料的屈服强度、切变模量均高于基体,且熔覆层的温度梯度大,所以导致残余应力较大.图3(b)为冷却600 s后横向残余应力分布云图,横向残余应力最大值出现在熔覆层边界与基体的交界处,σxmax=387 MPa.图5 至图7分别为熔覆层中心点横向、厚向和纵向的应变变化情况,包括热应变(Thermal)、弹性应变(Elastic)、塑性应变(Plastic)和总应变(Total)总应变等于热应变、弹性应变、塑性应变三者之和.在时间t=6.667 s到t=7.333 s之间,该点在激光光斑内,处于熔化状态,因此出现应变突然增大的情况.当激光光斑移除该点后,该点温度迅速降低,应变都迅速减小;进入弹性状态后,各向总应变均趋于平缓.由图5至图6可以看出,各向应变中,热应变最大,且总大于零;x、y向塑性应变小于零,z向塑性应变大于零.无论弹性应变和塑性应变符号怎么变化,总应变(total)总大于零.说明,热应变是决定材料处于拉伸或压塑塑性变形的决定因素;最终的横向X方向总应变小于零;厚度Y方向的应变趋于零;纵向Z方向的应变最大.图7 纵向应变随时间变化曲线由于熔覆层心部组织冷却速率要低于周围组织,心部凝固收缩受到周围组织的拘束作用.在三向应力中,纵向应力要远大于其它两向应力,最终表现为纵向产生拉应变,其它两向产生压应变,熔覆层中心的主应变图为第三类主应变图[7],如图8所示.周围组织对其有类似于挤压和拉拔的作用.这进一步说明了实际生产中熔覆层总是产生横向裂纹这一事实,熔覆层横向裂纹如图9所示.图8 第三类主应变图图9 熔覆层横向裂纹2.2 熔覆层应力随厚度梯度分布图1 0和图11为冷却600 s后沿工件中心厚度方向从上至下纵向应力和横向应力变化曲线图.由图10看出,熔覆层(高1 mm)纵向应力从上至下逐渐增大,到熔池附近达到最大值,这说明熔池附近是裂纹高发区,容易产生横向裂纹.过熔池后纵向应力逐渐减小,且始终为拉应力,对称线上的拉应力由两侧压应力补偿.图11表明,熔覆层上横向应力沿厚度方向从上至下是由压变拉,到熔池附近拉应力达到最大.基体则是上下两侧受拉,中间受压.2.3 激光熔覆应力机理分析激光熔覆过程产生的内应力是典型的由于材料温度变化而形成的热应力.送粉式激光熔覆过程中,可以近似认为所有形成熔覆层的粉末在达到基体之前全部熔化,到达基体后将热量传递给基体,同时基体也直接受到激光束的照射,并在光斑内形成熔池.处于液态的金属流动应力较小,其受热膨胀由于受到周围材料限制所产生的不均匀的压缩塑性变形和压应力可以忽略不计.而与熔池相邻的高温区(未熔化)产生热膨胀,膨胀受到周围材料的限制,产生不均匀的压缩塑性变形和压应力;在冷却过程中,高温区一定程度上又被拉伸而卸载;同时,熔池的冷却凝固过程中,也受到周围材料的限制,难以自由收缩,产生拉应力与拉伸塑性变形;此外,由于金属相变后体积发生变化,也会产生相应的相变应力;熔覆粉末与基体材料热膨胀系数、弹性模量的差异,也会导致内应力的产生.3 结论(1)熔覆层中心区域一直受到拉伸应力,产生拉伸塑性变形;熔池附近的基体先受到压应力,产生压缩塑性变形,然后随着熔覆的进行和冷却再受到拉应力,产生拉伸塑性变形.(2)熔覆层纵向应力最大,横向应力次之,厚度方向应力最小.三向拉应力当材料在塑性状态均表现为上升趋势,当材料进入弹性状态后,横向和厚度方向上的拉应力下降明显,纵向应力基本保持不变.(3)熔覆层中心靠近基体一侧是裂纹敏感区,该区域处于纵向拉伸、横向和厚向压缩的第三类主应变状态,周围组织对其有类似于拉拔和挤压的作用.总之,送粉式激光熔覆智能控制与数值模拟技术方兴未艾,还有很多不确定因素需要探索,希望本文的研究工作能为该领域的发展起到抛砖引玉的作用.参考文献【相关文献】[1]王仲任,郭殿俭,王涛.塑性成形力学[M].哈尔滨工业大学出版社,1989:35-36.[2]las K.Tsirbas K.Salonitis G..Chryssolouris.Ananalytical Model of the Laser Clad Geometry[J].Int JAdv Manuf Technol,2007,32:34-41.[3]Edson Costa Santos,Masanari Shiomi,Kozo Osakada.Rapid Manufacturing of Metal components by laser forming[J].International Journal ofMachineTools&Manufacture,2006,46:1459-1468.[4]Jichang Liu,Lijun Li.Study on Cross-section Clad Profile in Coaxial Single-pass Cladding with a lowpower laser[J].Optics&Laser Technology,2005,37:478-482.[5]Kai Zhang,Weijun Liu,Xiaofeng Shang.Research on the Processing Experiments of Laser metal Deposition shaping[J].Optics&Laser Technology,2007,39:549-557.[6]U.deOliveira,V.Ocelik,J.Th.M.DeHosson.Analysis of Coaxial Laser Cladding Processing Congditions[J].Surface&Coatings Technology,2005,197:130-135.[7]Yunchang Fu,A.Loredo,B.Martin.A theoretical model for laser and powder particles interaction during laser cladding[J].Journal of Materials ProcessingTechnology,2002,128:106-112.。
同步送粉式激光熔覆过程温度场数值模拟孙福臻;曲文峰;杨立宁;邵东强【摘要】为测定同步送粉式激光熔覆过程温度场分布情况,需采用ANSYS软件对温度场进行数值模拟计算.建立了同步送粉式激光熔覆过程的数学模型与有限元模型,通过ANSYS软件得到了温度场数值模拟计算结果.针对激光熔覆过程中的拉应力消除与裂纹控制等问题,提出了双面熔覆的预热方法,并通过ANSYS数值模拟的方法确定了最佳翻转时间.通过模拟结果与实验结果对比了双面熔覆方法的优越性.双面熔覆方法可为同步送粉式激光熔覆的预热处理提供了行之有效的方法.%ANSYS software is needed for numerical simulation of temperature field to determination the temperature field distribution of synchronous powder feeding laser cladding. The mathematical model and finite element model of synchronous powder feeding laser cladding was established,and the temperature field numerical simulation result was obtained by ANSYS software. Aiming at the problems of tensile stress relieving and crack control in the process of laser cladding,two-sided cladding method was proposed,and the best turnover time was confirmed through the ANSYS numerical simulation method. The superiority was contrasted through the simulation results and the experimental results.Two-sided cladding method can provide an effective method for preheating treatment of synchronous powder feeding laser cladding.【期刊名称】《机械设计与制造》【年(卷),期】2017(000)010【总页数】3页(P126-128)【关键词】激光熔覆;温度场;双面熔覆方法;最佳翻转时间;数值模拟;ANSYS【作者】孙福臻;曲文峰;杨立宁;邵东强【作者单位】机械科学研究总院先进制造技术研究中心,北京 100083;先进成形技术与装备国家重点实验室,北京 100083;机械科学研究总院先进制造技术研究中心,北京 100083;先进成形技术与装备国家重点实验室,北京 100083;机械科学研究总院先进制造技术研究中心,北京 100083;先进成形技术与装备国家重点实验室,北京100083;烟台泰利汽车模具制造有限公司,山东烟台 264670【正文语种】中文【中图分类】TH16;TG456.7激光熔覆是一种重要的材料表面改性技术，它是以高能密度的激光为热源在基材表面熔覆一层熔覆材料，使之与基材实现冶金结合，在基材表面形成与基材具有完全不同成分和性能的合金层的表面改性方法[1]。
