喷涂层残余热应力的有限元分析研究.pdf

喷涂层残余热应力的有限元分析 大连海事大学 严志军 高阳 f内容摘要}采用有限元方法对等离子喷涂材料在制备过程中产生的残余热应力进行了分 析。详细论述了基体和涂层的热膨胀系数、加工温度以及基体表面粗糙度对应力大小和分 布的影响。分析结论对实现等离子喷涂工艺参数优化具有指导意义。 关键词:等离子喷涉残余热应专有限元捧 等离子喷涂技术以其流焰温度高(可以高达13000K)，可以熔化象陶瓷类的高熔点材料 的特点，日益得到广泛的应用。但是在等离子喷涂材料冷却后，由于涂层内部应力作用、 使用不当等因素，可能造成涂层部分脱落，其特点是在局部点产麻点状涂层脱落。根据脱 落特点可以推断产生脱落的原因是由于在涂层中产生了不均匀的残余热应力，因此有必要 对喷涂层的热残余应力进行分析。 本文通过采用有限元((FiniteElementMethod)对等离子喷涂涂层在制备过程中产生的残 余热应力进行了线弹性分析。详细论述了涂层材料、基体表面形貌、喷涂温度等因素对涂 层应力大小及分布的影响，其结论对实现等离子喷涂工艺参数的优化控制具有重要意义。 1残余热应力计算模型 喷涂后的冷却过程是一个十分复杂的物理过程，在这个过程中，温度场、组织场和应 力场都在连续不断的变化，三者互相联系、互相影响。为了分析方便，本分析作如下假设: 忽略组织场的变化而把分析限于温度场和应力场的变化;喷涂的基体材料的厚度无穷大， 基体表面具有一定的粗糙度Ra:涂层厚度一定，涂层材料的热胀系数a=小于基体热胀系数 at口 根据以上假设选取一定长度和厚度的基体和其表面的涂层进行二维有限元网格划分， 如图所示。由于估计在涂层内部的应力分布情况比较复杂，所以划分单元时增大了涂层的 单元密度。 图1有限元网格划分 1.1温度场分析 无内热源二维温度场传热微分方程为: aT a(.aT) a .‘OT) PO二，=:一}k二一!十:甲}k-:一1 of ax}oxj 卿l卿少 式中p为密度，c为比热容，k为导热系数。 边界条件为对流换热。初始条件为两种情况考虑:(1)涂层的初始温度Tz。大于基体的 初始温度Tio;(2)涂层与基体的初始温度相等并同时冷却。 LZ应力场分析 应力场分析采用线弹性模型，应力变化需考虑温度产生的热应变影响和弹性应变的影 响。 弹性应变公式为: - 1 介 G Z dal= 式中el为弹性应变张量，G为剪切模量，Ql为应力张量· 热应变公式为: 衅=艺mkczkdT Mk和ak分别为不同组织的体积分数(k=1.2依次代表基体材料和 式中:r1h为热应变张It, 涂层材料)n热膨胀系数 总应变公式为: d--,j=de,十介derlh 式中eT}为总的应变张it, 2 热应力计算与分析 本文以有限元分析软件 Algor为平台进行涂层材料应 力分析计算。 首先假设基体(钢)的热膨 胀系数a.为6.5e-6,涂层的热 膨胀系数。2为1e-6;基体的初 图2 涂层应力分布 始温度为2000C,涂层的初始温度为3000C,同时冷却至室温(200C)时，在表面凸峰处的 VanMise应力分布如图2所示. 由图2可见.由于基体和涂层的热膨胀系数不同，初始温度不同，因而在冷却过程中 荃体和涂层之间由于热变形影响而产生残余热应力，热应力分布的基本规律是:最大应力 出现在基体和涂层之间偏向于基体一侧，如果应力较小，则涂层和基体之间产生徽观裂纹 并开裂的可能性较小，反之则增大:受表面形貌影响，基体和涂层之间对应表面凹谷处的 应力小干对应凸峰处的应力。 通过对最大主应力分布的进一步分析表明，由于涂层热膨胀系数小于基体的热膨胀系 数，喷涂后