响应面法与粒子群算法集成的激光熔覆工艺参数优化方法_胡言峰.pdf

年第 卷 月第 期机 械 科 学 与 技 术 ：收稿日期：基金项目：国家自然科学基金项目（）、重庆市自然科学基金联合基金重点项目（）、重庆市高校创新研究群体项目（）、重庆英才计划“包干制项目”（）及重庆英才计划（）作者简介：胡言峰（），硕士研究生，研究方向为激光增材再制造，通信作者：杜彦斌，教授，博士，胡言峰，杜彦斌，许磊，等响应面法与粒子群算法集成的激光熔覆工艺参数优化方法机械科学与技术，（）：响应面法与粒子群算法集成的激光熔覆工艺参数优化方法胡言峰，杜彦斌，许磊，周志杰，舒林森（重庆工商大学 制造装备机构设计与控制重庆市重点实验室，重庆；陕西理工大学 陕西省工业自动化重点实验室，陕西汉中）摘要：为制备高质量的熔覆层，保证再制造机械零部件再服役寿命，提出了一种响应面法与粒子群算法集成的激光熔覆工艺参数优化方法。该方法以熔覆层质量为优化目标，以工艺参数为优化变量，基于实验结果构建工艺参数与熔覆层质量间的响应面近似数学模型，使用粒子群算法对优化问题进行求解得到最优工艺参数组合。最后通过激光多道熔覆实验进行验证，结果表明该方法优化后的工艺参数能够有效改善熔覆层质量，进而节约实验成本、提高生产效率。关键词：激光熔覆；工艺参数；熔覆层质量；响应面法；粒子群算法中图分类号：；文献标志码：文章编号：（），（，；，）：，：；激光熔覆是一种在不改变基材芯部材质情况下，通过激光将合金粉末在基材表面熔融，并冷凝一层高性能（如高耐腐蚀性、耐磨性、抗氧化性等）的熔覆层，从而实现损伤修复及表面改性的绿色增材再制造技术，常用于损伤机械零部件再制造。激光熔覆是由工艺参数、自然环境、材料特征、设备性能等多因素耦合作用的复杂过程，所以熔覆层的机 械 科 学 与 技 术第 卷：质量随着因素的变化而变化，其中激光功率、扫描速度、送粉量、搭接率等工艺参数是决定熔覆层质量的关键因素，进而影响再制造零部件的再服役寿命，通过分析、调整、优选工艺参数，能够实现对熔覆层质量的有效控制。为此，针对制备高质量熔覆层的需求，通过优化得到最佳工艺参数组合对于保障并提高零部件再服役寿命具有重要意义。近年来，诸多高校、研究院主要通过统计学方法与智能算法两种方法对激光熔覆工艺参数优化展开了研究。等基于田口法实验结果，通过灰色关联分析得到最佳工艺参数组合。等基于正交实验结果，将 与 方法相结合，实现了工艺参数优化。蒋三生等通过极差分析法辅以激光熔覆实验结果分析工艺参数对熔覆层特征的影响，获得最优工艺参数组合。温海骏等结合全因子实验结果，构建了工艺参数与熔覆层综合质量间的反馈神经网络预测模型，并使用遗传算法实现寻优得到最佳的工艺参数。庞祎帆等基于正交实验结果建立工艺参数与沉积率间的遗传算法优化神经网络模型，同时采用遗传算法全局寻优得到最优工艺参数组合。等结合实验结果建立工艺参数与能量消耗、粉末利用率间的禁忌基因表达式编程模型，并使用非支配排序遗传算法得到最优工艺参数集。利用统计学方法实现工艺参数寻优的结果较可靠、操作简单，但存在一定的主观局限性，通过智能算法建模并寻优的方法系统性好、效率高，但模型对样本的需求量大，而激光熔覆是一个典型的小样本实验，因此模型的建立比较关键，会影响最终的优化结果。针对目前激光熔覆工艺参数优化存在的问题，本文引入响应面法（，）与粒子群算法（，），拟提出一种 与 集成的激光熔覆工艺参数优化方法。其中，响应面法能够利用较少的数据构建响应目标模型，同时还考虑了随机误差，改善了模型的精度，具有计算量小、且能够表征复杂的非线性关系的特点；粒子群算法是一种具有记忆能力的全局随机搜索算法，收敛速度快，全局搜索能力更强。该方法确定以熔覆层质量为优化目标，使用响应面法的中心复合实验设计法进行实验设计，根据实验方案进行激光熔覆实验，结合实验结果构建工艺参数与熔覆层质量目标间的响应面近似数学模型，并采用粒子群算法寻优获得最优工艺参数组合，从而有效地改善了熔覆层的质量。该方法可为生产实际中制备高质量熔覆层提供方法指导与决策参考。激光熔覆工艺参数优化模型建立 优化模型定义激光熔覆层质量是各种特征组成的综合体，决定了再制造零部件的再服役寿命，需要明确响应目标以表征熔覆层质量。熔覆层质量特征也随着工艺参数的变化而变化，但影响的工艺参数众多且相互耦合，因此将对其影响最大的工艺参数作为实验因素。由于表征熔覆层质量的响应目标有多个，各个子目标间可能会相互矛盾，即不能使子目标同时达到最优。因此需对各个子目标进行分析、比较、权衡使之成为一个单目标优化问题，数学表达式为：（）（）式中：为响应目标近似数学模型；为各个响应目标的对应权重，为响应目标的个数，；和 分别表示第 个实验因素的上下限，为实验因素的个数，。本文以表面平整度（）与稀释率（）表征熔覆层质量。表面平整度反映了多道熔覆表层的光滑程度，即实际可利用熔覆层的多少；稀释率是熔凝时基材渗入熔覆层导致成分的变化程度，也反映了熔覆层与基体结合质量，则将二者作为响应目标，计算公式为：（）（）式中：为熔覆层宽度；为熔覆层高度；为熔覆层面积；为熔池面积，如图 所示。图 激光多道熔覆层横截面示意图在激光熔覆实验中，激光功率（）、扫描速度（）、送粉量（）、搭接率（）对多道熔覆层的、影响最大，则将其作为实验因素。以表面平整度与稀释率作为工艺参数的质量响应目标，则加权使之成为一个最小单目标优化问题，数学表达式为：（）（）（）（）第 期 胡言峰，等：响应面法与粒子群算法集成的激光熔覆工艺参数优化方法：式中：（）为表面平整度近似数学模型；（）为稀释率近似数学模型，为表面平整度的最优值，愈趋近于 愈好，为稀释率的最优值，通常取.，、分别为表面平整度与稀释率对应权重，和 分别表示第 个实验因素的上下限，取。响应面近似数学模型构建由于工艺参数与熔覆层质量响应目标间是复杂的非线性关系，且激光熔覆实验是小样本实验，故需要保证通过实验结果回归拟合的模型能够近似准确的描述二者的关系。响应面法是由 和 提出的一种兼顾效率与精度的设计实验与构建模型的数理统计方法，其核心思想是按照恰当的实验设计方法设计的实验方案进行实验，以实验结果为基础利用简单的显式函数表达自变量（实验因素）与因变量（响应目标）间的复杂映射关系。为获得响应面近似数学模型，首先进行实验前准备，即确定实验所用的基材与熔覆材料。为避免相邻熔覆间的热影响，将大块基材切割成只进行一次熔覆实验的小样，并去除其表面氧化皮以及其它杂质，同时将熔覆用合金粉末干燥处理备用。确定实验因素（工艺参数）的取值范围。中心复合实验设计（，）是使用最多的响应面设计方法，通过 设计得到具体的激光熔覆实验方案，不仅能够减少实验次数，缩短实验周期，还能节省实验成本。以两水平析因设计点为基础，增加了保证模型精度的轴向极值点，同时为确保设计近似正交还增添了一定数量的中心点。以编码的形式编排，在 个实验因素的情况下，由中心点，析因点、轴向极值点 这 部分构成，实验时再转化为实际的操作值，实验次数 计算式为 （）式中：为析因设计点实验数；为轴向极值点实验数；为中心点处的重复实验数。在激光熔覆系统中按照实验方案进行实验。为防止熔覆过程中系统过热损坏全程打开冷却系统，同时通入保护气避免熔覆区域迅速被氧化；将激光发生器发出的激光作为热源，激光从激光头射出，辐照在基材上形成熔池，同时将从送粉器送出的合金粉末熔化；控制系统控制着激光头与工作台相对运动，最后在基材表面冷凝形成连续的熔覆层。对熔覆后的样品，通过切线割、打磨、抛光至镜面并腐蚀得到金相试样；使用三维超景深显微镜测量熔覆层的形貌特征；随后经分析处理得到响应目标数据。基于实验结果，选用多项式函数（如式（）所示）能够近似真实的表征实验因素与响应目标间的非线性关系，同时减少计算量、提升效率。以实验因素数据为输入，以响应目标数据为输出，在 中构建实验因素与响应目标间的近似数学模型，即（）（）式中：（）为响应目标函数；（，）表示第 个实验因素；，为待定各项回归系数；为误差项。构建的响应面近似数学模型直接影响着优化结果的正确性，需要对其进行验证分析。通过方差分析保证模型的显著性与可信度。一般地，模型的 值越大、.，失拟项的.，则认为回归拟合的响应面模型是显著的，能够较为真实反映实验因素与响应目标间的关系。残差是实验值与预测值的差值，基于残差分布情况分析响应面模型的适应性与准确性。引入决定系数、均方根误差 以及平均相对误差 评价响应面模型的精度，反映了模型的预测规律与真实规律的拟合程度，与 反映了模型预测值与实验值的差异程度，即 愈大和 与 愈小时，则回归拟合的响应面模型具有较高的精度，计算过程为：（）（）（）（）（）（）式中：为实验次数；为模型的计算预测值；为实验值；为均方根误差；为平均相对误差。为所有实验值的平均值。激光熔覆工艺参数优化问题求解粒子群算法是由 和 提出的一种全局随机并行优化算法。其核心思想是在解空间内，每个粒子在附近寻优得到个体最优解，然后通过种群信息共享，再群体迭代得到种群最优解，粒子的 与 影响粒子的速度，进而决定其搜索的方向和距离。粒子的速度与位置更新过程如图 所示。机 械 科 学 与 技 术第 卷：图 粒子更新过程计算式为：（）（）（）式中：，为种群的规模；，为迭代次数；，为粒子的维度；为第 次迭代粒子的飞行速度；为第 次迭代粒子的位置；，分别为个体与种群的学习因子；为惯性权重，调节搜索速度；，为随机常数。为有效地均衡粒子群算法的全局与局部的搜索能力，所以 不宜是一个固定的值，则采用惯性递减权重，即（）（）式中：为初始惯性权重；为最大迭代次数时的惯性权重。以式（）为适应度函数，将工艺因素的取值范围作为搜索范围，采用粒子群算法在搜索范围内迭代寻优，当达到最大迭代次数或者得到全局最优解时终止寻优，进而得到最优工艺参数组合。寻优过程如图 所示。图 基于粒子群算法的工艺参数优化求解流程图 实验验证与分析 实验设计与实验结果实验选用 钢为基材，粉末为熔覆材料，粉末形貌如图 所示，化学成分如表 所示。图 粉末的形貌表 钢与 各元素的质量分数 材料 钢 在 中根据 原理进行激光熔覆实验设计，工艺参数对应的编码水平数以及范围如表 所示，具体的实验方案如表 所示。表 激光熔覆工艺参数与编码水平数工艺因素水平（编码及真实值）（）（）.实验选用的激光熔覆系统如图）、图）所示，熔覆过程如图）所示。图 实验设备完成实验后制备金相试样，得到的 激光多道熔覆层横截面均无裂纹。选用 三维超景深显微镜测量多道熔覆层的、，测量结果如表 所示，按照式（）与式（）计算、，结果如表 所示。第 期 胡言峰，等：响应面法与粒子群算法集成的激光熔覆工艺参数优化方法：表 实验方案与实验结果序号 实验结果 响应面近似数学模型构建 表面平整度回归拟合模型在 中辅以表 中的工艺参数与表面平整度数据，利用式（）回归拟合表面平整度近似数学模型（）描述工艺参数与表面平整度间的映射关系为：机 械 科 学 与 技 术第 卷：（）（）构建表面平整度近似数学模型之后，对式（）进行方差分析，结果如表 所示。模型的 值为.，值小于.，表明由于噪音的影响而回归拟合出表面平整度近似数学模型的概率小于.，模型整体上是十分显著的，可信度高；失拟项的 值为.，相对于纯误差并不显著。所以表面平整度模型在回归区域内能够准确的描述工艺参数与表面平整度间的复杂关系。表 表面平整度模型方差分析结果方差来源 平方和 自由度均方 值 值模型 显著残差 失拟项 不显著纯误差 总和 图 图 中不同颜色的样本点代表不同工艺参数下的表面平整度，图 中所有样本点的残差正态概率近似呈一条直线分布，表明残差呈正态分布，表面平整度模型具有良好的适应性。图 中上下两条线代表残差的分布范围，所有表面平整度预测值的残差在分布范围内随机的分布在零线附近，表明表面平整度的残差之间无明显规律。图 中大部分样本点分布在直线（）上，小部分样本点分布在直线（）附近，表明表面平整度实验值与预测值十分吻合，模型拟合精度较高。综上，构建的表面平整度近似数学模型具有较高的适应性与准确性。图 表面平整度模型残差正态概率分布图图 表面平整度模型残差与预测值分布图图 表面平整度模型预测值与实际值分布图通过误差分析，可知决定系数 .，.，.，表明表面平整度近似数学模型的拟合程度高、预测精度高、稳定性好、实验值与预测值吻合。综上所述，获得的表面平整度近似数学模型非常显著，适应性好，可信度高，预测精度高，能够准确的描述工艺参数与表面平整度间的非线性关系。稀释率回归拟合模型以表 中的工艺参数与稀释率数据为基础，在 中利用式（）构建稀释率的近似数学模型（）表征工艺参数与稀释率间的非线性关系，即：（）（）构建稀释率近似数学模型之后，将式（）进行方差分析，结果见表。第 期 胡言峰，等：响应面法与粒子群算法集成的激光熔覆工艺参数优化方法：表 稀释率模型方差分析结果方差来源 平方和 自由度均方 值 值模型显著残差失拟项 不显著纯误差总和模型的.，.，表明仅有.的概率构建的稀释率近似模型会失效，意味着模型整体是显著的、可信度较高；模型失拟项的.，则相对于纯误差并不显著。因此，在回归区域内稀释率模型拟合程度高，能够通过工艺参数较为准确的预测其稀释率。图 图 中不同颜色的数据点代表不同工艺参数下的稀释率，图 中各个数据点的残差正态概率近似呈线性分布，表明所有数据点均是正常的，稀释率模型的适应性较好。图 稀释率模型残差正态概率分布图图 中上下两条线代表残差的分布范围，稀释率预测值的残差在分布范围内无规律的分布于零线附近，表明各个残差之间无明显联系，随机性好。图 稀释率模型残差与预测值分布图图 中小部分数据点分布在直线（）上，大部分数据点分布在直线（）周围，表明稀释率的实验值与预测值基本吻合，模型拟合精度较高。综上，构建的稀释率近似数学模型有较高的适应性与准确性。图 稀释率模型预测值与实际值分布图稀释率近似数学模型的决定系数.，.，.，表明模型的拟合程度较高，预测精度高。综上所述，回归拟合的稀释率近似数学模型的显著性高，适应性好，可信度高，拟合精度高，能够较为准确的反映工艺参数与稀释率间的非线性关系。工艺参数优化及结果取 .，.，即表面平整度为.，稀释率为.时，熔覆层质量最好。取、.，因此最终优化目标函数和约束条件为：（）（）（）|（）以式（）为适应度函数，以函数的最小值为寻优目标，在（）中以粒子群算法对激光功率、扫描速度、送粉量、搭接率这 个工艺参数进行优化，设置粒子群算法的初始参数如表 所示。表 粒子群算法的初始参数参数数值.粒子群算法在寻优过程中，仅需迭代 次便开始收敛，适应度值为，迭代 次时，适应度值为，趋近于，收敛精度最高，适机 械 科 学 与 技 术第 卷：应度曲线如图 所示，得出最优解.，.，.，.。即激光功率为 .、扫描速度为.、送粉量为.、搭接率为.。图 适应度曲线基于最优工艺参数的激光多道熔覆结果如下：表面平整度为.，即熔覆层表面较光滑，实际可利用熔覆层达到最多；稀释率为.，即熔覆层材料成分变化较小，且基材与熔覆层结合较好，不易脱落。结论）通过定义优化模型，以稀释率和表面平整度表征熔覆层质量并作为优化目标，并构建了响应面近似数学模型，最后采用粒子群算法对工艺参数优化问题求解，进而形成一种响应面法和粒子群算法集成的激光熔覆工艺参数优化方法。）辅以 钢为基体，合金粉末为熔覆材料的激光多道熔覆实验对该方法进行了验证分析，结果表明：所提出的工艺参数优化方法能够改善熔覆层质量，具有一定的实用性。）由于表征熔覆层质量的目标有多个，文中采用赋权法将其转化为单目标问题，存在一定的局限性，下一步研究将采用多目标优化算法求解工艺参数非劣解组合，从中挑选出相对最优解组合。参考文献，：徐国建，李春光，郭云强，等激光熔覆 混合粉末的熔覆层组织焊接学报，（）：，（）：（）邢如飞，许星元，黄双君，等激光沉积修复 钛合金微观组织及力学性能 材料工程，（）：，（）：（），：刘立君，刘大宇，王晓陆，等 钢激光熔覆陶瓷修复层的参数优化焊接学报，（）：，（）：（），：，：，（）：，（）：蒋三生，梁立帅，舒凤远 钢表面激光熔覆 基合金覆层工艺优化 材料导报，（）：，（）：（）温海骏，孟小玲，许向川，等基于神经网络和遗传算法的激光熔覆工艺参数多目标优化应用激光，（）：，（）：（）庞祎帆，傅戈雁，王明雨，等基于响应面法和遗传神经网络模型的高沉积率激光熔覆参数优化中国第 期 胡言峰，等：响应面法与粒子群算法集成的激光熔覆工艺参数优化方法：激光，（）：，（）：（），：，：唐忠，李文强，李彦一种区间不确定性参数的敏感度与可靠性分析方法计算机集成制造系统，（）：，（）：（）李聪波，杨青山，陈文倩，等面向响应性能的集成式电子液压制动系统执行器参数优化计算机集成制造系统，（）：，（）：（）朱刚贤，张安峰，李涤尘激光熔覆工艺参数对熔覆层表面平整度的影响 中国激光，（）：，（）：（）李伦翔，张德强，李金华，等激光熔覆镍基合金形貌优化及残余应力分析激光与光电子学进展，（）：，（）：（），：邓嵘，敬斌杰，张文汀响应面法在自激脉冲喷嘴结构优化中的应用机械科学与技术，（）：，（）：（）袁健基于响应面的汽车超高强度钢淬火参数优化金属热处理，（）：，（）：（）李莉，张赛，何强，等响应面法在试验设计与优化中的应用实验室研究与探索，（）：，（）：（）辛彬，李淑娟，李玉玺单晶硅电火花成形加工试验研究与工艺参数优化 兵工学报，（）：，（）：（），（）：王琪，张任，刘建，等基于粒子群算法的双层剪式液压升降台液压缸位置的多目标优化问题研究机械科学与技术，（）：，（）：（）单剑锋，杨雨粒子群优化的流形 模拟电路故障诊断机械科学与技术，（）：，（）：（）