机械创新实习报告

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