机械创新实习报告

目录目录 第一章 绪 论.- 2 - 1.1 齿轮的淬火处理.- 2 - 1.2 有限元方法.- 2 - 1.3 齿轮有限元仿真的现状.- 2 - 1.4 课题意义.- 3 - 1.5 研究内容.- 3 - 第二章 圆柱齿轮的几何特征和模型建立 .- 4 - 2.1 圆柱齿轮的几何特征.- 4 - 2.2 圆柱齿轮的三维模型建立.- 5 - 第三章 圆柱齿轮渗碳淬火过程热力学模型 .- 5 - 3.1 相关材料热力学参数.- 6 - 3.1.1 热学性能参数.- 6 - 3.1.2 力学性能参数.- 6 - 3.2 建立齿轮模型.- 6 - 3.3 淬火仿真方案及其相关参数.- 7 - 第四章 圆柱齿轮淬火过程的计算机仿真 .- 8 - 4.1 温度场分布的仿真.- 8 - 4.1.1 瞬态传热的数学模型.- 8 - 4.1.2 齿轮淬火温度场的仿真.- 8 - 4.1.3 温度场仿真.- 9 - 4.2 淬火应力分布的仿真.- 11 - 第五章 结论 .- 12 - 5.1 结论.- 12 - 参考文献.- 13 - 精选文库 - 2 - 第一章第一章 绪绪 论论 1.11.1 齿轮的淬火处理齿轮的淬火处理 齿轮传动是机器中最常见的一种机械传动，是传递机器动力和运动的一种 主要形式，是各种机械产品的重要基础零部件。它与带、链、摩擦、液压等机 械传动相比，具有功率范围大、传动效率高、圆周速度高、传动比准确、使用 寿命长、结构尺寸小等一系列特点。因此，它已成为许多机械产品不可缺少的 传动部件也是机器中所占比重最大的传动形式。齿轮的设计与制造水平将直接 影响到机械产品的性能和质量。由于齿轮在工业发展中的突出地位，致使齿轮 被公认为工业化的一种象征。 随着我国国民经济的快速发展，作为国家基础产业的机械设备也相应朝着 大型化方向发展。作为机械设备的重要元件齿轮传动装置的传递功率也越来 越大，产品规格尺寸也越来越大。为了满足机械设备传递动力的需要，同时减 少齿轮传动的结构尺寸和重量、提高齿轮的承载能力以及克服热处理技术和加 工机床带来的限制，最大化提升齿轮装置的综合性能指标，要求应用于机械中 的齿轮等机械零部件具有高强度、高可靠性。因此，为提高机械零部件的材料 强度，大多数采用各种热处理及表面处理等方法。目前，常常是通过淬火实施 表面硬化处理，以取代传统的齿轮调质处理。 1.21.2 有限元方法有限元方法 ANSYS是一个功能十分强大的有限元分析软件，不仅适用于常规工程问题的 静态或动态有限元分析，还能在诸如流体力学、热力学、温度场、电磁场等方 面进行有限元计算。 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算物体内部各 点的温度、并导出其它热物理参数。运用ANSYS软件可以进行热传导、热对流、 热辐射、相变、热应力及接触热阻等热问题的分析求解。 ANSYS不仅能解决纯粹的热分析问题，还能解决与热相关的其它诸多问题， 如热-应力分析、热-电分析、热-磁分析等。一般称这类涉及两个或多个物 理场相互作用的问题为耦合场分析。ANSYS提供了两种分析耦合场的方法：直接 耦合法和间接耦合法。 1.31.3 齿轮有限元仿真的现状齿轮有限元仿真的现状 传统的齿轮淬火都是利用表面硬化处理，或进行钢的调质处理，其分析过 程大都是通过实际测量，或经验进行对结果的处理和对比，并没有成熟的理论 和数据分析说明，存在很大的偶然性，而且由于实验条件的不稳定性，结果有 精选文库 - 3 - 时会出现很大的偏差，所以必须通过大量的实验，才能得出相对稳定的结果， 耗时又耗力。从有限元仿真分析的理论建立开始，经过很短的时间，便被各领 域进行了参数化分析和改进，目前，有限元分析仿真已经发展成为了一套成熟 的理论，应用到机械、建筑、电子等各个领域中，取得了很好的研究成果。 目前，在国外，通过对齿轮进行有限元仿真分析来获得齿轮进行热处理时 的各项参数和数据，已成为一种主流的求解参数的方法，特别是针对 ANSYS 分 析软件，不仅可以得到详细的数据和参数，而且可以看到清晰的仿真过程，便 于我们 对整个过程进行分析控制，所以在国外普及化程度很高；在国内，在齿 轮加工过程中，大都依然沿用传统的制作工艺，又加上对 ANSYS 软件的接触较 晚，普及率相对较低，总体上还在一种学习的过程中，相信随着国内外技术的 交流和进步，齿轮有限元分析会很快运用到国内的生产技术中，取得很大的发 展。 1.1.4 4 课题意义课题意义 采用齿面硬化处理制成硬齿面齿轮，是提高齿轮强度及承载能力的有效途 径，也是齿轮传动的主要发展趋势，目前大多数齿轮制造业发达国家已普遍采 用了硬齿面齿轮。在实现硬齿面的各种热处理工艺（渗碳淬火、氮化、表面淬 火等）中，渗碳淬火工艺虽然比较复杂，但在传递相同功率（扭矩）的情况下， 齿轮的减速器体积最小、重量最轻，整机价格最低，是生产应用中最主要的工 艺方法。为得到硬齿面齿轮，各国至今仍然采用机械加工（或塑性成形）-渗 碳-热处理的传统工艺存在如下问题： 齿轮表面渗碳层厚度不一致，不利于齿轮综合性能的提高。渗碳处理是在 齿轮切削或塑性成形后进行，由于没有确立渗碳层控制技术，齿根与齿面在同 一渗碳氛围中进行渗碳处理，齿面、齿顶及齿根的渗碳层浓度、梯度、厚度大 致相同。然而由于齿面与齿根的工作特性不同，要求齿面渗厚、齿根渗碳层薄。 为了保证齿面的耐磨性，使含碳量达到或超过0.18%。渗碳层深度与模数之比 0.120.125，而有效硬化层与齿轮模数之比为0.1时，齿曲疲劳强度最高；因此 为了保证齿根弯曲强度，则渗碳层的含碳量和渗碳层深度就满足齿面接触疲劳 强度的要求。在实际生产中，很难二者兼顾，达到理想状态。为了防止轮齿硬 化层的剥落，有效硬化层深度应不小于最大剪应力深度的1.5倍。 1.1.5 5 研究内容研究内容 1）选择汽车圆柱齿轮为研究对象，获取其工艺尺寸参数，为建立数学模型 提供参考； 2）根据其尺寸参数，利用SolidWorks软件进行三维建模，再把模型导入 精选文库 - 4 - ANSYS软件中，建立热力学参考模型； 3）查阅必要的资料，获取齿轮材质20CrNi2MoA在201000的热、力学参 数，为后继的热、力学分析提供依据； 4) 对圆柱齿轮实体在淬火过程中的模型进行计算机仿真,简要地对大模数 圆柱齿轮淬火变形进行分析。 第二章第二章 圆柱齿轮的几何特征和圆柱齿轮的几何特征和模型建立模型建立 2.2.1 1 圆柱齿轮的几何特征圆柱齿轮的几何特征 本课题选用 齿轮的几何特性 如下: 图图2-12-1圆柱齿轮示意图圆柱齿轮示意图 表表 2-12-1：圆柱齿轮参数及技术要求：圆柱齿轮参数及技术要求 技术要求 精选文库 - 5 - 1，齿面淬火处理，有效硬化层深度 5-5.6mm；齿面硬度 HRC57-64；齿心硬 度 HRC35-40； 2，机械性能 1080Mpa；785Mpa；s8%；35%；AX47Jcm2； b s 3，淬火后对齿部进行喷丸处理； 4，进行探伤检查，齿部进行磁粉探伤，不允许存在线性磁痕显示，整体进 行超声波探伤，内部不允许有白点、裂痕，同时内部质量应符合 GBT5000.15-1988 标准规定的级要求； 2.22.2 圆柱齿轮的三维模型建立圆柱齿轮的三维模型建立 首先，根据齿轮的各项数据和参数，利用 SolidWorks 软件对齿轮进行三维 模型的建造，图 2-2 即为该齿轮的三维模型。 图图 2-22-2 用用 SolidWorksSolidWorks 创建的齿轮三维模型图创建的齿轮三维模型图 由于齿轮各轮齿参数相同，热处理的边界条件相同，生产中对大模数齿轮 一般采用逐步淬火的方法，即对齿轮单齿进行逐一淬火，为了能够准确和方便 地建立几何模型，可选取齿轮的 1/54 部分，即单齿进行模型建立，这样可以简 化计算模型，节省内存空间，提高仿真速度，其仿真结果与实际完整齿轮基本 一致。图 2-3 即为用 SolidWorks 创建的单齿的三维模型。 图图 2-32-3 用用 SolidWorksSolidWorks 建立的齿单齿的三维模型建立的齿单齿的三维模型 生成模型后，将当前模型另存为*.x_t 格式，留作之后导入 ANSYS 软件。 精选文库 - 6 - 第三章第三章 圆柱齿轮渗碳淬火过程热力学模型圆柱齿轮渗碳淬火过程热力学模型 20CrNi2MoA低碳合金钢是大模数齿轮首选的渗碳钢之一，长期以来具有稳 定的化学成分、成熟的冶炼工艺过程。其热力学性能参数已作为标准列入国家 渗碳钢材料手册。 3.13.1 相关材料热力学参数相关材料热力学参数 3.1.13.1.1 热学性能参数热学性能参数 圆柱齿轮材料选用 20CrNi2MoA，当利用 ANSYS 软件对其淬火过程进行温度 场仿真时，需要提供材料密度、比热容、热传导率等参数，其对应的参数数值 如表 3-1【8】： 表表 3-13-1 20CrNi2MoA20CrNi2MoA 热学性能参数热学性能参数 密度（Kg/m3） 7850 比热容 Cp【J/(Kg K) 】 20 460 温度 t() 01002004006008001000 热传导率 【W（m K） 】 14.716.61820.823.526.328.2 3.1.23.1.2 力学性能参数力学性能参数 同样的，当利用 ANSYS 软件对齿轮的淬火过程进行应力场仿真时，也需要 提供在不同温度下的弹性模量、屈服强度、切变模量、泊松比、线膨胀系数等 力学性能参数，见表 3-2： 表表 3-23-2 20CrNi2MoA20CrNi2MoA 力学性能参数力学性能参数 温度 t () 弹性模量 E （Pa） 屈服强度 （Pa） 切变模量 G （Pa） 泊松比线膨胀系数 1 201.93e110.2e91.93e10 5001.5e110.933e91.5e10 10000.7e110.435e90.7e10 15000.1e110.07e90.1e10 0.291.78e-5 精选文库 - 7 - 20000.01e110.007e90.01e10 3.23.2 建立齿轮模型建立齿轮模型 利用 ANSYS 建立分析模型最常用的方法有两种：一种是利用 ANSYS 的绘图 功能，生成分析模型，再进行仿真模拟；另外一种，就是利用 ANSYS 与其它三 维绘图软件，如 AutoCAD，SolidWorks，Pro/E 等强大的数据传输能力，本课题 就是利用 SolidWorks 软件进行建模，然后导入 ANSYS 软件中，具体步骤为：运 行 ANSYS 软件， 【File】【Import】【PARA】,弹出窗口中查找到之前保存 的*.x_t 文件,选中， 【OK】 ，导入完成，现在看到的是线框，接着【Ploctrls】 【Style】【Solid Model Facets】,下拉框中选择【Normal Faceting】 【OK】,接着鼠标右键，选择【Replot】即可看到实体，如图 3-1： 图图 3-13-1 ANSYSANSYS 建模分析图建模分析图 3.33.3 淬火仿真方案及其相关参数淬火仿真方案及其相关参数 依赖有限元方法可以仿真淬火齿轮内部温度变化过程，这样可以得到齿轮 内部的仿真温度场，然后计算、仿真齿轮的淬火变形情况。但是，首先需要确 定初始的温度边界条件和对流换热系数。 齿轮齿面在油性介质中淬火，其介质主要在上下深度方向流动（既齿宽方 向），其流速一般取0.7-1.5m/s的中间值，将粗车齿轮放入电热炉中加热到 780-820，接着通过行车将工件立即置入室温淬火油井中，冷却2小时后取出， 根据长期经验数据取淬火油平均温度30。 数据：淬火油为AN32淬火油，使用温度为30，淬火前齿轮加热温度800； AN32淬火油：黏度=32X10-6（Pas）,热导率=0.1276W（mK）,密 m m 度=882.5m3。 mc 我们根据研究需要将淬火分为2个阶段：淬火变形期和淬火冷却期。本文使 精选文库 - 8 - 用经典传热理论进行相似计算获得淬火变形期不同阶段（我们根据研究需要取 0-15秒）的平均对流换热系数红hl,随着淬火冷却的进行,齿轮和淬火油温差的减 少, hl开始衰减【12】。齿轮的变形主要发生在淬火变形期,这期间齿轮和淬火油 温差最大,对流剧烈(且可能存在相变时，流体吸收或放出汽化潜热，其局部对 流换热系数有可能大于hl)，因此齿轮内部将产生较大的内应力，导致齿轮的不 规律的变形。我们把有限元方法引入齿轮淬火的应力应变分析中，结合最优化 方法中的试探法，逐步掌握齿轮淬火变形趋势，为生产加工提供数据参考和改 良依据。 第四章第四章 圆柱齿轮淬火过程的计算机仿真圆柱齿轮淬火过程的计算机仿真 4.14.1 温度场分布的仿真温度场分布的仿真 4.1.14.1.1 瞬态传热的数学模型瞬态传热的数学模型 当齿轮淬火时，会发生组织转变（奥氏体可能分解为铁素体、珠光体、贝 氏体、马氏体）同时释放潜热，虽然释放潜热不像熔化或凝固时潜热那么大， 但也是不可忽略的因素，潜热释放将使得整个仿真过程波动和非线性变化，为 研究带来难度。根据Fourier定律，运用能量守恒原则，对于物体温度随时间而 变（非稳态问题） 4.1.24.1.2 齿轮淬火温度场的仿真齿轮淬火温度场的仿真 ANSYS能够很好地完成瞬态热分析，瞬态 传热分析计算一个系统地随时间变化的温度场 及其他热参数，在工程上一般用其计算温度场 并作为热载荷进行应力分析。瞬态传热分析中 的载荷是随时间变化的，必须将载荷-时间曲 线分为载荷步，载荷-时间曲线中的每一个拐 点为一个载荷步，如图所示对于每个载荷步， 必须定义载荷值及时问值，同时必须选择载荷 步为简便或阶越。ANSYS瞬态热分析基本步骤包括构建模型、施加载荷、求解与 后处理。步骤1 建模的基本流程如图4-1 精选文库 - 9 - 图图4-14-1 ANSYSANSYS建模流程图建模流程图 图图4-24-2 热力学模型热力学模型的网格划分图的网格划分图 根据以上步骤，设置工件名为chilundanchi，定义材料密度 ,比热，热传导率取800时， 3 /7820mKgKKgJCp/460 ，为得到较精确的仿真结果，划分网格密度为2，经以上步骤，KmW/ 3 . 26 得出模型的网格划分图如图4-2 步骤2 载荷求解的基本流程： 图图4-34-3 ANSYSANSYS仿真载荷求解流程图仿真载荷求解流程图 根据步骤，定义分析类型为模型三维实体的温度场仿真，设置初始温度为 800，淬火介质温度为30，载荷设置为30秒，每载荷步为1秒。 步骤3 后处理的基本流程： 图图4-44-4 ANSYSANSYS仿真数据仿真数据 后处理流程图后处理流程图 由于齿轮渗碳后淬火的特殊性，整个瞬态过程的不稳定结果，使得我们无 法给出确切的即时温度场，所有的仿真都是对可能持续时间的估计和推断，也 就是说，我们只能仿真出淬火温度场出现的趋势和范围，而无法给出某点的即 时温度。 精选文库 - 10 - 通常淬火是为了强化齿面硬度、增加耐磨性，同时增加齿轮内部应力，产 生应变。淬火增效强化的时间通常在开始淬火的几十秒内，但由于淬火冷却的 实际情况比较复杂，淬火增效强化的时间范围难以把握，为了研究需要，我们 选取05秒、010秒、015秒为可能出现的淬火增效强化范围，逐一进行温 度场仿真。淬火增效进行之后的温度场变化对齿轮齿面性能影响较小，产生的 内应力及应变小于前者，基本上属于典型的热胀冷缩类型，故不作为研究对象。 以下仅列出05秒，010秒和015秒的温度场、应力场仿真分析。 4.1.34.1.3 温度场仿真温度场仿真 齿廓和内圈边界受淬火油作用，从被加热的780820骤冷至与淬火油边 界的平衡温度（不等于平均温度），考虑到液体流动速度，冷却时间和平衡温 度需要根据计算过程设定，轮齿 两侧结合边界属于对称边界。这 样即可建立相应的温度场计算模 型，通过23次迭 代计算，求出轮齿与淬火油边界 的平衡温度，然后进行最后的温 度场计算。图4-5、4-6、4-7分 别为圆柱齿轮5秒,10秒和15秒淬火温度场。随着时间延长，冷却层深入轮齿内 部。 图图4-54-5 圆柱齿轮圆柱齿轮5 5秒淬火温度场秒淬火温度场 图图4-64-6 圆柱齿轮圆柱齿轮1010秒淬火温度场秒淬火温度场 精选文库 - 11 - 图图4-74-7 圆柱齿轮圆柱齿轮1515秒淬火温度场秒淬火温度场 由仿真图可以看出齿心与齿面的温差变化，在5秒时，T1=53；在10秒 时，T2=79；在15秒时，T3=87。随着温度的推移，齿心与齿面的温差 逐渐增大，其热应力也随之增大，所以在淬火刚开始的这个时间段内，是累积 热应力的主要过程。 由于在淬火的过程中，齿面与齿心的温差较大，就整个齿轮而言，尺寸较 小的轮齿部分淬火过程中受作用最为明显，会引起不同程度的尺寸和形状误差， 所以在齿轮的加工过程中，要根据仿真的结果，分析其综合变形量，应用到实 际的生产加工中，为齿轮的下一步工艺加工，取一个合适的叫加工余量。 通过对温度场仿真结果进行分析，我们可以清晰地得出齿轮淬火过程中各 时间段的温度变化，各部分间的温差范围，推算出齿轮的各部分的淬火变形方 式和形变量，同样，根据对仿真结果个时间段的对比分析，可以较快捷准确地 确定齿轮淬火过程中的最佳时间，为齿轮淬火分析提供理论依据。 4.24.2 淬火应力分布的仿真淬火应力分布的仿真 由于齿轮的不同部分之间的热膨胀系数不匹配，在淬火的急速冷却过程中， 彼此的收缩不一致，从而导致热应力的产生，同时也达到了齿面增效强化的目 的。热应力问题实际又是两个物理场之间的相互作用，故属于耦合场问题。 ANSYS提供了两种分析热应力的方法：直接法和间接法。 分析 1 分析 2 结果 1 结果 2 图图4-84-8 间接法分析数据流程图间接法分析数据流程图 精选文库 - 12 - 使用间接法的步骤基本上是先进行热分析，重新进入前处理，然后设置结 构分析，读入热分析的节点温度后再 设置参考温度，最后进行求解和后处 理。我们采用后者，先进行热分析获 得淬火温度场分布，然后将求得的节 点温度作为体载荷加到结构应力分析 中。淬火过程产生不均匀的应力分布， 齿廓应力变化明显，应力集中在齿根 部位，图4-9、4- 10、4-11分别为5秒、10秒和15秒轮齿增效强化的应力分布。 图图4-94-9 圆柱齿轮圆柱齿轮5 5秒淬火热应力秒淬火热应力 图图4-104-10 圆柱齿轮圆柱齿轮1010秒淬火热应力秒淬火热应力 图图4-114-11 圆柱齿轮圆柱齿轮1515秒淬秒淬火热应力火热应力 精选文库 - 13 - 通过对图4-9进行分析，可以得到以下几点结论： 1）齿廓应力变化明显，应力集中在齿根部位，且齿根部位应力最大，极大 值为