固体表面激光加工熔池特性的数值分析.pdfVIP

应用基础与工程科学学报 第 9卷 223期 Vol. 9 ,No. 223 2001年 9月 September 2001 JOURNAL OF BASIC SCIENCE AND ENGINEERING 文章编号 :1005209302001202? 0320215207 中图分类号 :TK124 文献标识码 :A 固体表面激光加工熔池特性的数值分析 杨立新 , 彭晓峰 清华大学热能工程系 ,北京 100084 摘要 :通过引入液相分数标量描述界面熔融区 ,能量方程中将潜热转换为热焓源 项 ,而动量方程则引入‘Darcy’ 源项 ,由此建立统一的液和固相控制方程 ,实现在 固定网格下求解具有移动界面的熔化和凝固过程.应用该模型详细分析不锈钢 表面激光加工过程热传导、 Marangoni 和热浮力三种换热机制对熔化过程和熔池 形状的影响. 关键词 :熔化 ,激光 ,Marangoni流 ,数值分析 高能束流 High Energy Density Beam 通常缩写为 HEDB加工技术是利用高能量密 度的束流 激光束、 电子束、 等离子 作为热源 ,对材料或构件进行特种加工或表面处理的 技术 ,包括焊接、 切割、 打孔、 喷涂、 表面改性和刻蚀、 精细加工等. HEDB 经历了近半个世 纪的研究和开发 ,在 20世纪 80年代呈现出加速发展的趋势.在我国‘九五’经济发展关键 技术规划中 ,该技术也入选作为支柱产业之一的制造技术的重点发展方向. HEDB 加工技 术的基础研究涉及到传热、 传质、 流动动力、 材料、 计算机等多学科 ,呈现出多学科协同和 交叉的特点.目前国内在这方面的基础研究比较薄弱 ,往往为得到某种加工工艺理想的加 工结果 ,需要多次地实验才能得到一组适合的工艺参数.显然 ,研究高能束作用下的熔体 流动特征 ,如溶池形状、 温度梯度、 局部冷却速率和凝固的晶体组织等 ,是一个更好地理解 和掌握液体相变和凝固组织形成规律的基础.快速、 准确的分析方法会对工艺参数的选择 和控制将起到良好的理论指导作用. [1 ,2 ] 熔池的形成和发展是典型的固、 液相变传热问题.从 20世纪 70年代后期 ,Sparrow [3 ] 等人开始考虑液相区中对流对相变过程的影响.彭晓峰等人 应用小扰动分析方法 ,对激 光加热过程初始熔池的形成和 Marangoni 流在熔池形成发展过程中的影响给出了机理性 的解释.近年来 ,随着计算机和数值方法的发展 ,更多的研究者致力于通过数值模型来模 [4 ,5 ] 拟 HEDB 加工过程的溶池细节特性 .这些数值方法可以划分为固定网格法和移动网格 法两大类.前者物理意义清晰 ,但存在一些缺点 ,如 :流动区初始溶池的确定问题 ,不同的 [6 ] 初始溶池形状的假设将决定不同的数值模拟结果 ;对合金熔化时 Mushy区处理的困难 ; 巨大的网格重新划分的工作量等.后者虽然不存在这些问题 ,但在 Mushy区模型建立以及 计算收敛性上加大了数值模拟的难度.然而由于缺乏相关的 HEDB 加工实验结果的验证 , 使这些数值模型很难在实际工艺中得到广泛应用. 收稿日期 :2001203205 ;修订日期 :2001206208 ? 基金项目 :国家自然科学基金No.清华大学博士论文基金资助项目 作者简介 :杨立新1969~ ,男 ,博士生1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd All rights reserved.216 应用基础与工程科学学报 Vol. 9 本文通过引入液相分数标量描述界面熔融区 ,能量方程中将潜热转换为热焓源项 ,而 动量方程则引入‘Darcy’ 源项 ,由此建立统一的液和固相控制方程 ,实现在固定网格下求 解具有移动界面的熔化和凝固过程.应用该模型对 AISI 304 不锈钢材料激光加工过程进 行详细分析 ,考察热传导、 Marangoni 和热浮力三种换热机制对熔化过程和熔池形状的影 响. 1 数学物理模型 本文的分析模型如图 1所示 ,为长 10mm ,宽 8mm ,高 2mm的 AISI 304 不锈钢板 ,钢板 上面作用激光束.在高能量的激光作用下金 属开始熔化形成熔池 ,金属熔池的熔体对流 驱动力主要来自两种不同的机制 :表面张力 梯度引起的 Marangoni 流和熔池水平温度梯 度决定的自然对流.熔化过程实质上是一个 移动界面的相变问题 ,这里将导出可以在固 定网格下统一求解的数值模型.数值模型中 所做的主要假设包括 : 1 流动为不可压层 图 1 分析模型 [7 ,8 ] 流 ; 2 导热系数、 液体密度、 表面张力和 Fig.