激光毛化表面温度与应力的数值模拟与实验研究.pdfVIP

182 材料热处理学报 第29卷 表面非激光直接照射区： 数值计算过程中，采用简洁耦合法求解热应力。 3’7 首先，利用solid70单元计算脉冲激光作用的温度场。 (3) 一后；÷=^，(r一死)+盯。，￡(r一凭) u，‘ 在此基础上进行温度梯度，冷却速度的计算，并以温 其中，口为激光功率密度；^，为对流换热系数；盯。为 Stefan．Boltzman常数；疗为热交换面的外法线方向；￡ 法计算应力场。加载分为时间和温度两部分。采用 为材料的吸收率。 渐变的载荷变化方式。通过在脉冲激光作用前施加 1．2计算过程 一个加热半径较大，而能量密度较小且均匀的工况来 激光毛化示意图见图1。为简化计算，假设作用 模拟预热的作用”o。 于工件表面的脉冲激光为Gauss分布热源，光斑直径 O．2mm，激光功率400～1500w，脉冲宽度0．1。2IIls， 2 结果与讨论 x 初始温度20℃。工件尺寸15mm10mm×5mm。移 2．1激光毛化温度场分布 动速度10mm／s。对试样进行三维网格划分，单元总 图2为熔凝层不同位置点的温度一时间曲线，其 数1640，材料为Gcrl5轴承钢，密度78lokg／m3，泊松 中，z表示熔凝层不同位置点与工件表面的距离。可 比O．3，热导率35w／mK，比热450J／kg·K，热膨胀系数 见，工件表面与表面以下100肚m处的最高温度分别 1．1×lO。／℃，弹性模量208GPa。假设材料的物理性 质不随温度变化，固态和液态材料对激光的吸收系数 度，如图3(a)所示，随着激光功率的增加，熔凝层在 不变。通过选择不同的激光功率、脉冲宽度和预热温 并、，，、：三个方向的温度梯度G。、G，、G：都有所增加。 度计算试样表面及不同深度处的温度以及应力状况。 其中在深度方向的温度梯度最大，达到106—107℃／m 量级，大于另外两个方向的温度梯度一个数量级。特 p 童 三 芒 ＆ E 占 图l激光毛化示意图 Fig．1