淬火过程应力场的计算机模拟-常州铁道高等职业技术学校.doc

淬火过程应力场的计算机模拟 史 东 丽， 朱 菊 香 （常州铁道高等职业技术学校电气工程系 江苏 常州 213011） 摘 要：利用有限元模拟软件对42Crmo试棒在淬火过程中的应力变化进行数值模拟，使得在淬火过程中的热应力和残余应力得以体现，并分析淬火过程中表面和心部的应力变化。 关键词：淬火；数值模拟；有限元分析 1﹒引言 淬火是一项能够明显提高金属工件机械性能的热处理工艺。但是淬火过程中热应力和组织应力综合作用常导致工件各部位塑性变形不一致并产生残余应力，对材料性能有很大的影响，严重时可导致工件变形甚至开裂 [1]。对于此类问题，通常“看不见，摸不着”，只能通过失效分析经验判断，缺乏科学的依据和预见性。淬火过程计算机数值模拟的出现，把过去靠经验进行生产变为科学地指导生产，对提高淬火质量、优化工艺、减少浪费具有现实意义。 2﹒淬火过程数值计算 2.1研究对象及试验条件 选用φ50×200 的试棒，材料选用42Crmo。试棒经从室温20℃加热到860℃保温1小时后，开始冷却，冷却介质为水，温度为20℃。考虑到试棒出炉到进入淬火介质之间的过程，设定出炉后在空气中冷却5s. 然后淬水直至冷却至室温。利用淬火过程中的冷却曲线反算求解得到表面换热系数[2] 表1表面换热系数 Surface Convection Coefficient 2.2 淬火过程数值计算 由于φ50的试棒为轴对称，所以选取1/2 截面进行网格划分，共200个单位，表面网格划分的细一些，心部则相对粗一些，可增加计算精度，并减少运算时间。并依次在表面下1mm、2mm、11mm、22mm选取4点，作为跟踪点，如图1。心部及表面各点的冷却曲线如图2。 图1 网格划分及跟踪点 图2 各跟踪点的冷却曲线 Chart 1: General mesh and track point Chart 2: The cooling curves of track point 3﹒计算结果 图3为试棒加热及淬冷过程中4个追踪点的应力曲线图。 图4为试棒1/2轴截面淬冷到最后的轴向残余应力等值线图。 图5为试棒1/2轴截面淬冷到30S(相变点)时的应力等值线图。 图6为试棒1/2轴截面淬冷到50S(应力转换)时的应力等值线图。 图3 为试棒加热及淬冷过程中 图4 为试棒1/2轴截面淬冷到 4个追踪点的应力曲线图 最后的轴向残余应力等值线图 Chart 3: The stress curves of 4 track point Chart 4: The residual stress isoline of during heat up and quench process 1/2 section specimen after quenched 图5 试棒1/2轴截面淬冷到 图6 为试棒1/2轴截面淬冷到 30S时的应力等值线图 50S时的应力等值线图 Chart 5: The stress isoline of 1/2 Chart 6: The stress isoline of 1/2 section after quench in 30 sec section after quench in 50 sec 4 结果分析及讨论 4.1根据表面及心部各追踪点在淬火过程中的应力曲线（图3），得知在淬火冷却时，试件纵截面Z轴方向上的表面位置，应力先呈拉应力后呈压应力(如图3，P1追踪点和P2追踪点的轨迹曲线)，而心部处应力先呈压应力后呈拉应力(如图3，P3追踪点和P4追踪点的轨迹曲线)。这是因为在试样淬火初期，试样表面收缩量大于心部，表面产生拉应力，中心产生压应力。表面和心部的内应力都随着淬火时间的延续而不断增大，经历一定时间后达到最大值。随后两部分的内应力都逐渐减小。随着淬火时间的延长，试样表面内应力变为压应力，试样中心的内应力变为拉应力，而且数值不断增大到一定值。 4.2 根据轴向残余应力等值线图，可得知试棒轴截面各个位置的残余应力（如图4），残余应力过高会显著影响产品的机械性能。通过有限元分析，使得通常不容易通过试验得到的残余应力数值通过计算机模拟展现出来。 4.3 通过图3可观察到，在淬火进行到30S时，（如图5）表面拉应力逐渐减小，心部压应力也逐渐减小，在淬火进行道50S时，（如图6）表面拉应力和心部压应力都达到最小值。这说明从淬火30S开始发生了相变，相变首先从边界开始，因此边界应力状态逐