Ti(C,N)基金属陶瓷干...件下瞬态温度场的有限元模拟_周黎明.pdf

第 卷 第 期 年 月 成都大学学报（自然科学版）（）文章编号：（）：收稿日期：基金项目：国家大学生创新创业训练计划项目（）；四川省大学生创新创业计划项目（）；四川省粉末冶金工程技术研究中心开放基金课题项目（）作者简介：周黎明（），男，博士，讲师，从事摩擦磨损及有限元仿真研究：通信作者：袁 容（），女，博士，副研究员，从事材料疲劳强度及有限元仿真研究：（，）基金属陶瓷干摩擦条件下瞬态温度场的有限元模拟周黎明，冉龙姣，葛巍东，杨向莙，袁 容，曾红志（成都大学 机械工程学院，四川 成都）摘 要：利用有限元数值模拟方法对（，）基金属陶瓷在摩擦过程中的瞬态温度场进行了分析，探讨了不同摩擦转速、载荷和环境温度等因素对（，）基金属陶瓷摩擦表面温度变化的影响 结果显示，（，）基金属陶瓷在摩擦过程中其表面最大温升随着转速和载荷的增加而增大，且转速对表面温升的影响大于载荷的作用，而环境温度对金属陶瓷表面温升的影响并不显著；在滑动过程中，（，）基金属陶瓷表面最高温度区域主要位于其与对磨小球的接触区；利用摩擦试验机对（，）基金属陶瓷表面摩擦温度变化进行了实验测量，结果与有限元分析相一致关键词：（，）基金属陶瓷；滑动摩擦；瞬态温度场；有限元分析中图分类号：文献标志码：引 言（，）基金属陶瓷是以（，）为基体，以 等过渡金属为粘结剂，并添加、和 等碳化物，采用粉末冶金方法制备的一种先进硬质材料，已在数控切削刀具、滑动轴承和喷嘴等耐磨领域得到了应用 等分别对干摩擦、水润滑和油润滑条件下（，）基金属陶瓷与 钢的对磨情况进行了探究，发现油润滑条件下其耐磨性最好，而干摩擦时其磨损状况最剧烈 等研究发现，原料粉末粒度越细，（，）基金属陶瓷耐磨损性能越好 虽然针对（，）基金属陶瓷的摩擦磨损研究已见诸多报道，但至今尚未形成系统研究，特别是摩擦温升对（，）基金属陶瓷摩擦性能影响方面的研究还较为缺乏 众所周知，由摩擦引起的表面温升变化对材料的磨损行为有着非常重要的影响 虽然有关摩擦温升问题的理论研究已有报道，但这些理论模型均存在一定的局限性且不能准确地反映出（，）基金属陶瓷材料在摩擦过程中温度的分布及变化状态 随着计算机技术的快速发展，有限元数值模拟方法在求解与分析摩擦过程的温升问题中得到了广泛应用 例如，等利用有限元数值模拟方法研究了齿轮传动过程中其表面瞬态温度场的分布，其分析结果在改善齿轮表面磨损方面提供了一定依据 等采用有限元数值模拟研究了螺旋锥齿轮啮合摩擦发热和瞬态热特性，揭示了螺旋锥齿轮在传动过程中的瞬态温度场分布，为螺旋锥齿轮副接触模式的改善提供了指导目前，有关（，）基金属陶瓷在摩擦过程中的瞬态温升问题研究较少，本研究利用有限元软件 对（，）基金属陶瓷在滑动过程中的瞬态温度场进行数值模拟，着重考查速度、载荷和环境温度等因素对该金属陶瓷瞬态摩擦温升的影响，并结合实际的摩擦试验验证，为后续（，）基金属陶瓷材料的设计及摩擦磨损机制分析提供一定的理论依据 模型与材料参数 有限元模型图 为球盘摩擦试验示意图 试验过程中，氧化铝（）小球作为对磨小球通过专用夹具保持周向固定，球上施加法向载荷 由于摩擦过程中对磨小球与试样盘属于高副接触，且基于节省模拟计算时间考虑，本研究对（，）基金属陶瓷的摩擦有限元模型进行了简化，如图 所示，对磨小球采用半球模型 简化后，有限元分析模型中试样盘外径为，内径为 ，厚度为 ，对磨球半径为 为了保证计算收敛精度，本研究对试样盘采用了渐进式网格划分，即从试样盘边缘至旋转轴心网格越来越细 网格划分后，金属陶瓷试样盘单元数为 个，节点数为 个，对磨半球单元数为 个，节点数 个 模拟过程中，试样盘和对磨小球的单元类型均采用六面体的 单元 此外，在有限元计算过程中假定条件如下：）材料各向同性，且其物性参数不受温度变化影响；）球盘接触面摩擦系数设定为 且保持不变；）不考虑材料表面的热辐射效应施加载荷旋转轴Al2O3小球试样盘图 球盘摩擦示意图图 有限元三维分析模型 材料物性参数为便于整个模拟过程的计算与分析，本研究假设试样材料组分均匀且完全致密（，）基金属陶瓷和 对磨小球的物理性能参数 见表表 （，）基金属陶瓷和 小球的物理性能参数材料杨氏模量 泊松比热膨胀系数（）热导率（）密度（）（，）基金属陶瓷 小球 结果与分析 转速对摩擦温升的影响假定初始环境温度为，施加载荷为 ，不同转速条件下金属陶瓷在摩擦过程中的瞬态温升变化情况如图 所示 在转速设定为 时，金属陶瓷接触表面摩擦温度随着时间逐渐上升，当 时，其表面摩擦温度大约为 ，即摩擦温升大约为 由图可知，随着转速的增加，金属陶瓷表面温升逐渐增大，当转速从 提高至 ，时的摩擦温升将由 跃升至，这表明转速对摩擦过程温度的变化具有显著的正相关影响 此外，从图中还可以看出，当转速由 提高至 时，虽然 内的温升也由 上升到了 ，但上升幅度明显减缓，这说明当转速提高到一定速度后，金属陶瓷的表面的温升变化趋势将会有所收敛，其原因可能是在较高的转速条件下，材料摩擦表面的对流换热也会增强，从而削弱了表面温升持续大幅升高的态势807060504030010 20 30 40 50 60时间/s温度/25 rad/s50 rad/s75 rad/s100 rad/s初始温度：25载荷：20 N图 转速对瞬态摩擦温升的影响 载荷对摩擦温升的影响假定初始环境温度为，转速为 ，图 显示的是不同载荷条件下金属陶瓷在摩擦过程中瞬态温升的变化情况 由图可知，随着施加载荷的增大，金属陶瓷接触表面温升也逐渐增大 当载荷分别为、和 时，金属陶瓷表面的最大温升 第 期 周黎明，等：（，）基金属陶瓷干摩擦条件下瞬态温度场的有限元模拟分别为 、和 分析可以发现，虽然载荷增大可以使金属陶瓷表面摩擦温度上升，但随着载荷的持续增加，其表面温升变化也将趋于减缓 根据赫兹接触理论可知，个物体在受压情况下实际接触表面积将随着载荷的增加而增大虽然增加的接触区域也会因摩擦产生热量而提高材料表面的温升，但同时也会增大对流换热的作用，故而金属陶瓷表面最高温度增加的幅度也会有所缓和504540353025010 20 30 40 50 60时间/s温度/10 N20 N30 N40 N初始温度：25转速：50 rad/s图 载荷对瞬态摩擦温升的影响 此外，结合图 和图 可以发现，虽然转速和载荷的增大都可以使金属陶瓷表面摩擦温升增大，但在相同的（载荷与转速的乘积）条件下，相较于载荷，转速对金属陶瓷表面温升的影响将更加显著例如，对比条件（转速 、载荷 ）和条件（转速 、载荷 ），种条件的 值一样，但表面最大温升分别是 和 初始环境温度对摩擦温升的影响假定转速为 ，施加载荷为 ，不同初始环境温度条件下金属陶瓷在摩擦过程中的瞬态温升变化情况如图 所示 可以看到，在 值特定的条件下，不同初始环境下金属陶瓷摩擦表面的最高温度均随时间延长而升高，且近似线性变化 当初始环境温度分别为、和 时，内金属陶瓷表面最大摩擦温升分别为 、和 ，即随着初始环境温度的升高，金属陶瓷在摩擦过程中的表面最大温升出现了略微的增大这可能是在较高的环境温度条件下，金属陶瓷摩擦表面的对流换热作用受到了部分抑制，从而增强了其表面的摩擦温升状况 为进一步分析初始环境温度的影响，本研究在此引入了温升变化率的概念，公式为，（）式中，为摩擦温升变化率，；为单位时间内的最大温升，；为单位时间，010 20 30 40 50 60时间/s温度/25100200400载荷：20 N转速：50 rad/s40032024016080图 初始环境温度对瞬态摩擦温升的影响 结合图 可知，不同初始环境温度条件下的金属陶瓷表面最大摩擦温升变化率分别为 、.、和 可以看到，其温升变化率相差不大，也就是说，初始环境温度的变化对金属陶瓷瞬态摩擦温升的影响并不显著 摩擦瞬态温度场的分布图 为初始环境温度 、转速 和载荷 时的（，）基金属陶瓷摩擦过程中瞬态温度场的分布云图 可以看出，摩擦过程中的瞬态最高温度位于（，）基金属陶瓷试样盘和 对磨小球接触区域，且沿着试样盘旋转方向（顺时针方向），金属陶瓷表面的摩擦温度逐渐减小 此外，由于自身的热传导作用，金属陶瓷试样盘表面温度也由摩擦区域沿半径方向以对称方式向四周递减TEMP（平均：75%）38.44338.36838.25438.12738.06237.98737.88737.75237.611图 金属陶瓷摩擦瞬态温度场分布云图 瞬态温升的实验测量为了验证有限元模拟的准确性，本研究利用奥地利 公司制造的 型摩擦磨损试验机对该金属陶瓷自润滑轴套进行了球盘摩擦磨损实验 实验采用的金属陶瓷试样盘尺寸为 ，对磨球为 的 小球 测试条件为，载荷 、转速 、初始温度 和成都大学学报（自然科学版）第 卷 摩擦半径 摩擦过程中的试样表面的温度变化由试验机自带的温度传感器测得图 展示了金属陶瓷试样盘表面温度随时间的变化情况 可以看出，金属陶瓷试样表面温度会随着摩擦时间增加而逐渐上升 此外，从图中可以发现，利用有限元模拟方法得到的对磨环表面温度结果略高于实验所测量的数值，这种差异的主要来源一方面是不可避免的测量误差，即在摩擦过程最高温度主要位于试样与对磨小球的接触区域，而试验机是无法对两者的直接接触区域进行测量；另一方面，为了便于有限元计算处理，模拟中对一些约束状态等条件进行了简化 虽然有限元结果与实测数值有所差异，但 种方法得到的摩擦温度变化规律基本一致，且 者温升数值基本近似，这表明有限元结果具有一定的可靠性 此外，相较于理论计算和实验测量等方法在预测或评估摩擦温度变化方面，有限元数值模拟对摩擦过程中的瞬态温度场分布情况可以进行更加直观地展现，这也是其独特的优势所在3936333027010 20 30 40 50 60时间/s温度/有限元结果实验测量值载荷：20 N；转速：50 rad/s；初始温度：25 图 实验测量值与有限元结果对比 结 论本研究利用有限元数值模拟方法研究了转速、载荷和初始环境温度等摩擦过程因素对（，）基金属陶瓷摩擦瞬态温度场的影响，主要结论如下：）（，）基金属陶瓷在摩擦过程中其表面最大温升随着转速和载荷的增加而增大，且相同的 条件下，转速对表面温升的影响要大于载荷的作用而在 条件不变的情况下，初始环境温度对表面最大温升的影响作用并不明显）（，）基金属陶瓷在摩擦过程中的最高温度区域主要位于试样盘与对磨小球的接触区，且试样盘表面温度沿着旋转方向和半径方向逐渐降低）摩擦试验测量结果显示，（，）基金属陶瓷试样盘摩擦表面温度变化趋势与有限元模拟结果相一致，从而验证了有限元模拟分析的可靠性参考文献：刘宁（，）基金属陶瓷材料 合肥：合肥工业大学出版社，（）：，（）：周黎明，吴明霞，廖富强，等 石墨（，）基金属陶瓷梯度自润滑材料的残余应力的有限元模拟 机械工程材料，（）：，（，），（）：，（，），（）：，：，（）：，（）：，：温诗铸，黄平 摩擦学原理 北京：清华大学出版社，（）：，（）：，：第 期 周黎明，等：（，）基金属陶瓷干摩擦条件下瞬态温度场的有限元模拟 ，：，（）：刘红兵，陶杰，张平则，等 功能梯度 涂层残余热应力分析 机械工程学报，（）：接触力学 徐秉业，罗学富，刘信声，等译 北京：高等教育出版社，（实习编辑：姚运秀）（，），（，）：（，）（），（，）（，），（，），（，）（，），：（，）；（上接第 页），：；成都大学学报（自然科学版）第 卷