基于ANSYS的斜齿轮淬火温度场和应力场仿真分析

1、20122012 届本科毕业论文（设计）届本科毕业论文（设计） 论文题目：论文题目：基于基于 ANSYSANSYS 的斜齿轮淬火温度场和应力的斜齿轮淬火温度场和应力 场仿真分析场仿真分析 学生姓名：学生姓名： 所在院系：所在院系： 机电学院机电学院 所学专业：所学专业： 机械设计制造及其自动化机械设计制造及其自动化 导师姓名：导师姓名： 完成时间：完成时间： 20122012 年年 5 5 月月 1010 日日 基于基于 ANSYSANSYS 的斜齿轮淬火温度场和应力场的仿真分析的斜齿轮淬火温度场和应力场的仿真分析 摘摘 要要 齿轮传动是现代机器中最常见的一种机械传动，是传递机器动力和运动的

2、一种主要形式，是各种机械产品的重要基础零部件。在齿轮传动过程中，需要有 较好的硬度和耐磨性，这可以通过齿轮的淬火来完成，但在淬火过程中，齿轮的 尺寸和形状都会出现不同程度的改变，特别是大型淬火齿轮的变形，我们需要进 行很好地控制，但这种研究，需要受到条件的限制，较高的成本和较长的实验时 间，通过计算机仿真，可以大大提高研究效率，降低实验费用，还可以仿真在不 同工艺参数下的实验，获得合理可行的最佳的工艺参数。 本课题利用了 ANSYS 有限元仿真分析软件，对齿轮淬火过程的温度场和应 力场进行了仿真，得出了在不同时间段内淬火过程的温度和应力分布图，并进行 了数据的分析对比，得出了影响淬火效果的各因

3、素之间的关系，为之后对齿轮淬 火变形的研究分析提供了数值依据。 关关键词键词：斜斜齿圆齿圆柱柱齿轮齿轮，淬火，有限元，淬火，有限元，计计算机仿真算机仿真 Temperature Field and Stress Field in the Quenching Progress of Helical Gear Based on ANSYS Abstract Gear transmission is one of the most common mechanical transmission in modern machines ,its a main form of transferring po

4、wer and sport of machine ,and its the important foundation of spare parts in all kinds of mechanical products .In the process of gear transmission ,it needs good hardness and wear resistance .This can be completed ,through the quenching ,but during the quenching ,the size and shape will change ,espe

5、cially the quenching of big gear .We need to control it perfectly ,but this study needs the limit of the condition ,more money and longer time .Through the computer we an save much time and money ,and it can get the most suitable results. This subject uses ANSYS software to simulate the quenching pr

6、ogress of gear of the temperature field and stress field ,which concludes that the different period quenching progress temperature and stress distribution and the analysis of the data contrast, and the influence of quenching effect that the relationship between the various factors , As for the gear

7、after quenching deformation of numerical analysis provides the basis. Keywords :Helical-Spur Gear，Quenching，Finite Element，Computer Simulation 目目录录 第一章 绪 论.- 1 - 1.1 齿轮的淬火处理.- 1 - 1.2 有限元方法.- 1 - 1.3 齿轮有限元仿真的现状.- 2 - 1.4 课题意义.- 2 - 1.5 研究内容.- 3 - 第二章 斜齿圆柱齿轮的几何特征和模型建立.- 3 - 2.1 斜齿圆柱齿轮的几何特征.- 3 - 2.2

8、斜齿圆柱齿轮的三维模型建立.- 5 - 第三章 斜齿圆柱齿轮渗碳淬火过程热力学模型.- 6 - 3.1 相关材料热力学参数.- 6 - 3.1.1 热学性能参数.- 6 - 3.1.2 力学性能参数.- 7 - 3.2 建立齿轮模型.- 7 - 3.3 淬火仿真方案及其相关参数.- 8 - 第四章 斜齿圆柱齿轮淬火过程的计算机仿真.- 10 - 4.1 温度场分布的仿真.- 10 - 4.1.1 瞬态传热的数学模型.- 10 - 4.1.2 齿轮淬火温度场的仿真.- 11 - 4.1.3 温度场仿真.- 13 - 4.2 淬火应力分布的仿真.- 15 - 第五章 全文总结.- 18 - 5.1

9、 结论.- 18 - 5.2 研究不足及展望.- 18 - 致 谢.- 20 - 参考文献.- 21 - - 0 - 第一章第一章 绪绪 论论 1.1 齿轮齿轮的淬火的淬火处处理理 齿轮传动是近代机器中最常见的一种机械传动，是传递机器动力和运动的 一种主要形式，是各种机械产品的重要基础零部件。它与带、链、摩擦、液压等机 械传动相比，具有功率范围大、传动效率高、圆周速度高、传动比准确、使用寿命 长、结构尺寸小等一系列特点。因此，它已成为许多机械产品不可缺少的传动部 件也是机器中所占比重最大的传动形式。齿轮的设计与制造水平将直接影响到 机械产品的性能和质量。由于齿轮在工业发展中的突出地位，致使齿轮

10、被公认为 工业化的一种象征。 随着我国国民经济的快速发展，作为国家基础产业的机械设备也相应朝着 大型化方向发展。作为机械设备的重要元件齿轮传动装置的传递功率也越来 越大，产品规格尺寸也越来越大。为了满足机械设备传递动力的需要，同时减少 齿轮传动的结构尺寸和重量、提高齿轮的承载能力以及克服热处理技术和加工 机床带来的限制，最大化提升齿轮装置的综合性能指标，要求应用于机械中的齿 轮等机械零部件具有高强度、高可靠性。因此，为提高机械零部件的材料强度， 大多数采用各种热处理及表面处理等方法。目前，常常是通过淬火实施表面硬化 处理，以取代传统的齿轮调质处理【1】【2】。 1.2 有限元方法有限元方法 A

11、NSYS是一个功能十分强大的有限元分析软件，不仅适用于常规工程问题 的静态或动态有限元分析，还能在诸如流体力学、热力学、温度场、电磁场等方 面进行有限元计算。 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算物体内部 各点的温度、并导出其它热物理参数。运用ANSYS软件可以进行热传导、热对流、 热辐射、相变、热应力及接触热阻等热问题的分析求解。 ANSYS不仅能解决纯粹的热分析问题，还能解决与热相关的其它诸多问题， 如热-应力分析、热-电分析、热-磁分析等。一般称这类涉及两个或多个物理场 相互作用的问题为耦合场分析。ANSYS提供了两种分析耦合场的方法：直接耦合 法和间接耦合法。

12、直接耦合法的耦合单元包含与相应的耦合场相关的所有必须的自由度，采 用耦合单元仅仅通过一次求解就能得出耦合场的分析结果。这种方法实际上是 - 1 - 通过计算包含所有必须项的单元矩阵或单元载荷向量来实现的。间接耦合法又 称序贯耦合法，主要通过第一次场分析的结果作为第二次场分析的载荷来实现 两种场的耦合。热应力耦合分析就是将热分析得到的节点温度作为载荷施加 到后继的应力分析中来实现耦合【3】【4】。 1.3 齿轮齿轮有限元仿真的有限元仿真的现现状状 传统的齿轮淬火都是利用表面硬化处理，或进行钢的调质处理，其分析过程 大都是通过实际测量，或经验进行对结果的处理和对比，并没有成熟的理论和数 据分析说明

13、，存在很大的偶然性，而且由于实验条件的不稳定性，结果有时会出 现很大的偏差，所以必须通过大量的实验，才能得出相对稳定的结果，耗时又耗 力。 从有限元仿真分析的理论建立开始，经过很短的时间，便被各领域进行了参 数化分析和改进，目前，有限元分析仿真已经发展成为了一套成熟的理论，应用 到机械、建筑、电子等各个领域中，取得了很好的研究成果。 目前，在国外，通过对齿轮进行有限元仿真分析来获得齿轮进行热处理时的 各项参数和数据，已成为一种主流的求解参数的方法，特别是针对 ANSYS 分析 软件，不仅可以得到详细的数据和参数，而且可以看到清晰的仿真过程，便于我 们 对整个过程进行分析控制，所以在国外普及化程

14、度很高；在国内，在齿轮加工 过程中，大都依然沿用传统的制作工艺，又加上对 ANSYS 软件的接触较晚，普 及率相对较低，总体上还在一种学习的过程中，相信随着国内外技术的交流和进 步，齿轮有限元分析会很快运用到国内的生产技术中，取得很大的发展。 1.4 课题课题意意义义 采用齿面硬化处理制成硬齿面齿轮，是提高齿轮强度及承载能力的有效途 径，也是齿轮传动的主要发展趋势，目前大多数齿轮制造业发达国家已普遍采用 了硬齿面齿轮。在实现硬齿面的各种热处理工艺（渗碳淬火、氮化、表面淬火等） 中，渗碳淬火工艺虽然比较复杂，但在传递相同功率（扭矩）的情况下，齿轮的减 速器体积最小、重量最轻，整机价格最低，是生产

15、应用中最主要的工艺方法。为 得到硬齿面齿轮，各国至今仍然采用机械加工（或塑性成形）-渗碳-热处理的传 统工艺存在如下问题： 齿轮表面渗碳层厚度不一致，不利于齿轮综合性能的提高。渗碳处理是在齿 轮切削或塑性成形后进行，由于没有确立渗碳层控制技术，齿根与齿面在同一渗 - 2 - 碳氛围中进行渗碳处理，齿面、齿顶及齿根的渗碳层浓度、梯度、厚度大致相同。 然而由于齿面与齿根的工作特性不同，要求齿面渗厚、齿根渗碳层薄。为了保证 齿面的耐磨性，使含碳量达到或超过0.18%。渗碳层深度与模数之比0.120.125， 而有效硬化层与齿轮模数之比为0.1时，齿曲疲劳强度最高；因此为了保证齿根弯 曲强度，则渗碳层

16、的含碳量和渗碳层深度就满足齿面接触疲劳强度的要求。在实 际生产中，很难二者兼顾，达到理想状态。为了防止轮齿硬化层的剥落，有效硬 化层深度应不小于最大剪应力深度的1.5倍。出于安全性等考虑，一般规定齿廓中 点处的有效硬化层深度不小于节点处最大剪应力深度的2倍。正常工艺条件下， 齿廓的硬化层深度从齿顶到齿根平缓过渡（减小），齿根处的渗层深,其淬火变形 量难以计算，使得必须提前留有充分的加工余量，保证淬火后、磨削前有足够的 磨削量，实际生产中的磨削量过多、磨削难度较大，直接影响到生产成本。故有 必要建立轮齿的数学分析模型，计算热处理过程中的渗碳层厚度和加工余量问 题，以此为生产工艺的参数设定提供参考

17、【5】。 1.5 研究内容研究内容 1）选择某减速器设备上斜齿圆柱齿轮为研究对象，获取其工艺尺寸参数，为 建立数学模型提供参考； 2）根据其尺寸参数，利用SolidWorks软件进行三维建模，再把模型导入 ANSYS软件中，建立热力学参考模型； 3）查阅必要的资料，获取齿轮材质20CrNi2MoA在201000的热、力学参 数，为后继的热、力学分析提供依据； 4) 对斜齿圆柱齿轮实体在淬火过程中的模型进行计算机仿真,简要地对大模 数斜齿圆柱齿轮淬火变形进行分析。 第二章第二章 斜斜齿圆齿圆柱柱齿轮齿轮的几何特征和的几何特征和模型建立模型建立 2.1 斜斜齿圆齿圆柱柱齿轮齿轮的几何特征的几何特征

18、 本课题选用齿轮的几何特性如下: - 3 - 图图2-1斜斜齿圆齿圆柱柱齿轮齿轮示意示意图图 表表 2-1：斜：斜齿圆齿圆柱柱齿轮齿轮参数及技参数及技术术要求要求 渐开线圆柱齿轮特性表 序号 项目名称 代号 数值 1 齿数 Z87 2模数Mn28 3 齿形角 an20 4 基本 参数 齿高系数 Ha/Hf1/1.4 5 分度圆柱螺旋角 8.5 6 螺旋线方向 LR左 7 变位系数 X-0.224 8 精度等级 GB10095-886KL 9 齿距累积公差 Fp0.16 - 4 - 10 齿形公差 ff0.042 11 齿距极限偏差 fpt0.036 12 齿向公差 F0.028 13 公法线跨

19、越齿数 K10 14 公法线理论长度 WK815.805 15 公法线平均长度极限偏差 Ews/EWI-0.42-0.51 技术要求 1，齿面淬火处理，有效硬化层深度 5-5.6mm；齿面硬度 HRC57-64；齿心硬度 HRC35-40； 2，机械性能1080Mpa；785Mpa；s8%；35%；AX47Jcm2； b s 3，淬火后对齿部进行喷丸处理； 4，进行探伤检查，齿部进行磁粉探伤，不允许存在线性磁痕显示，整体进行 超声波探伤，内部不允许有白点、裂痕，同时内部质量应符合 GBT5000.15- 1988 标准规定的级要求； 5，齿形两端棱边及齿顶沿齿长棱边倒角 245； 6，棱边倒角

20、 245 2.2 斜斜齿圆齿圆柱柱齿轮齿轮的三的三维维模型建立模型建立 首先，根据齿轮的各项数据和参数，利用 SolidWorks 软件对齿轮进行三维模 型的建造【6】，图 2-2 即为该齿轮的三维模型。 - 5 - 图图 2-2 用用 SolidWorks 创创建的建的齿轮齿轮三三维维模型模型图图 由于齿轮各轮齿参数相同，热处理的边界条件相同，生产中对大模数齿轮一 般采用逐步淬火的方法，即对齿轮单齿进行逐一淬火，为了能够准确和方便地建 立几何模型，可选取齿轮的 1/87 部分，即单齿进行模型建立，这样可以简化计算 模型，节省内存空间，提高仿真速度，其仿真结果与实际完整齿轮基本一致【7】。 图

21、 2-3 即为用 SolidWorks 创建的单齿的三维模型。 - 6 - 图图 2-3 用用 SolidWorks 建立的建立的齿单齿齿单齿的三的三维维模型模型图图 生成模型后，将当前模型另存为*.x_t 格式，留作之后导入 ANSYS 软件。 第三章第三章 斜斜齿圆齿圆柱柱齿轮齿轮渗碳淬火渗碳淬火过过程程热热力学模型力学模型 20CrNi2MoA低碳合金钢是大模数齿轮首选的渗碳钢之一，长期以来具有 稳定的化学成分、成熟的冶炼工艺过程。其热力学性能参数已作为标准列入国家 渗碳钢材料手册。 3.1 相关材料相关材料热热力学参数力学参数 3.1.1 热热学性能参数学性能参数 斜齿圆柱齿轮材料选用

22、 20CrNi2MoA，当利用 ANSYS 软件对其淬火过程进 行温度场仿真时，需要提供材料密度、比热容、热传导率等参数，其对应的参数 数值如表 3-1【8】： - 7 - 表表 3-1 20CrNi2MoA 热热学性能参数学性能参数 密度 （Kg/m3）7820 比热容 Cp【J/(Kg K) 】 20 460 温度 t()01002004006008001000 热传导率 【W（m K ）】 14.716.61820.823.526.328.2 3.1.2 力学性能参数力学性能参数 同样的，当利用 ANSYS 软件对齿轮的淬火过程进行应力场仿真时，也需要 提供在不同温度下的弹性模量、屈服强

23、度、切变模量、泊松比、线膨胀系数等力 学性能参数，见表 3-2： 表表 3-2 20CrNi2MoA 力学性能参数力学性能参数 温度 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 - 8 - 3.2 建立建立齿轮齿轮模型模型 利用 ANSYS 建立分析模型最常用的方法有两种：一种是利

24、用 ANSYS 的绘 图功能，生成分析模型，再进行仿真模拟；另外一种，就是利用 ANSYS 与其它三 维绘图软件，如 AutoCAD，SolidWorks，Pro/E 等强大的数据传输能力，本课题就 是利用 SolidWorks 软件进行建模，然后导入 ANSYS 软件中，具体步骤为：运行 ANSYS 软件， 【File】【Import】【PARA】,弹出窗口中查找到之前保存的*.x_t 文 件,选中， 【OK】，导入完成，现在看到的是线框，接着【Ploctrls】【Style】【Solid Model Facets】,下拉框中选择【Normal Faceting】【OK】,接着鼠标右键，选择

25、 【Replot】即可看到实体【9】，如图 3-1： 图图 3-1 ANSYS 建模分析建模分析图图 3.3 淬火仿真方案及其相关参数淬火仿真方案及其相关参数【 【10】【 】【11】 】 依赖有限元方法可以仿真淬火齿轮内部温度变化过程，这样可以得到齿轮 内部的仿真温度场，然后计算、仿真齿轮的淬火变形情况。但是，首先需要确定 初始的温度边界条件和对流换热系数。 齿轮齿面在油性介质中淬火，其介质主要在上下深度方向流动（既齿宽方向） - 9 - ，其流速一般取0.7-1.5m/s的中间值，将粗车齿轮放入电热炉中加热到780-820， 接着通过行车将工件立即置入室温淬火油井中，冷却2小时后取出，根据

26、长期经 验数据取淬火油平均温度30。 数据：淬火油为AN32淬火油，使用温度为30，淬火前齿轮加热温度800； AN32淬火油：黏度=32X10-6（Pas）,热导率=0.1276W（mK）,密度 m m =882.5m3。 mc 根据理论公式： （3-1） 3/18 . 0 Pr23500Re037 . 0 mlmm Nu m m Nu l h （3-2） 式中：tm-特征温度； x-特征尺寸； l-板长。 淬火油与温度高于其相应压力饱和温度的壁面接触时传热属于对流传热， 罗森诺通过大量实验数据的整理，最后得出下列计算式： = （3-3） r ttc swpt )( n l vl wl r

27、q CPr )( 33 . 0 1 通过上式可以计算出壁面热流密度q,其中： cpt-饱和液体的比定压热容，J（kgK）; tw-壁面温度， ts-饱和液体温度，； r汽化潜热，JKg； q壁面热流温度，Wm2; -饱和液体动力黏度，Pas; l -液气表面张力，Nm; -饱和液体密度，Kgm3; l -相应压力下饱和蒸汽密度，Kgm3; v n指数，对于水，n=1,对于其他液体，n=1.7; Cwr取决于高温表面和液体组合情况的经验系数。 根据以上理论公式，可以得到Num=1.76105,ho=3.68104W/(m2K)。 我们根据研究需要将淬火分为2个阶段：淬火变形期和淬火冷却期。本文使

28、用经 - 10 - 典传热理论进行相似计算获得淬火变形期不同阶段（我们根据研究需要取0-15秒） 的平均对流换热系数红hl,随着淬火冷却的进行,齿轮和淬火油温差的减少, hl开 始衰减【12】。齿轮的变形主要发生在淬火变形期,这期间齿轮和淬火油温差最大,对 流剧烈(且可能存在相变时，流体吸收或放出汽化潜热，其局部对流换热系数有 可能大于hl)，因此齿轮内部将产生较大的内应力，导致齿轮的不规律的变形。我 们把有限元方法引入齿轮淬火的应力应变分析中，结合最优化方法中的试探法， 逐步掌握齿轮淬火变形趋势，为生产加工提供数据参考和改良依据。 第四章第四章 斜斜齿圆齿圆柱柱齿轮齿轮淬火淬火过过程的程的计

29、计算机仿真算机仿真 4.1 温度温度场场分布的仿真分布的仿真 4.1.1 瞬瞬态传热态传热的数学模型的数学模型 当齿轮淬火时，会发生组织转变（奥氏体可能分解为铁素体、珠光体、贝氏体、 马氏体）同时释放潜热，虽然释放潜热不像熔化或凝固时潜热那么大，但也是不 可忽略的因素，潜热释放将使得整个仿真过程波动和非线性变化，为研究带来难 度。根据Fourier定律，运用能量守恒原则，对于物体温度随时间而变（非稳态问题） 并且内部具有相交潜热qv时的Fourier导热微分方程,可推导出直角坐标系下的传 热方程: +qv= （4-1） 2 2 2 2 2 2 z T y T x T t T cp 式中:-热传

30、导率; -材料密度; Cp-材料的定压比热; T-任意位置的温度; t-时间, x,y,z-直角坐标。 温度场随时间发生变化的传热过程为非稳态传热。这类传热又分为周期性传 热和非周期性传热两种。周期性传热过程中，导热物体内部温度以一定的规律 随时问周期性地变化。而在非周期性传热过程中，物体内部温度随时问不断升高 - 11 - 或降低,并经历相当时间后逐渐趋于周围介质的温度而最终海达到平衡，这类传 热过程又称为瞬态传热【13】。对于瞬态传热过程，系统的温度、热流率、热边界条 件以及内能都将随时问发生变化，根据（4-1）瞬态传热可以用以下公式表示： （4-2） QTKTC k 式中： K-传导矩阵

31、，包含导热系数，对流系数及辐射率和形状系数： C-比热容矩阵； Tk-温度对时间的导数； Q-节点热流率向量，包含热生成。 4.1.2 齿轮齿轮淬火温度淬火温度场场的仿真的仿真【 【14】 】 ANSYS能够很好地完成瞬态热分析，瞬态传热分析计算一个系统地随时间 变化的温度场及其他热参数，在工程上一般用其计算温度场并作为热载荷进行 应力分析。瞬态传热分析中的载荷是随时间变化的，必须将载荷-时间曲线分为 载荷步，载荷-时间曲线中的每一个拐点为一个载荷步，如图所示对于每个载荷 步，必须定义载荷值及时问值，同时必须选择载荷步为简便或阶越。ANSYS瞬态 热分析基本步骤包括构建模型、施加载荷、求解与后

32、处理。 步骤1 建模的基本流程如图4-1： 图图4-1 ANSYS建模流程建模流程图图 确定工件名，工作 标题，单位 定义单元类型并 设置选项 定义单元 实常数 划分网格 导入几何 模型 定义材料热性能参数（热传 导率，密度，比热） - 12 - 根据以上步骤，设置工件名为chilundanchi，定义材料密度, 3 /7820mKg 比热，热传导率取800时，为得到较精KKgJCp/460KmW/ 3 . 26 确的仿真结果，划分网格密度为2，经以上步骤，得出模型的网格划分图如图4-2： 图图4-2 热热力学模型的网格划分力学模型的网格划分图图 步骤2 载荷求解的基本流程： 定义分析类型 定

33、义均匀的温度场，设置初始温度 - 13 - 图图4-3 ANSYS仿真仿真载载荷求解流程荷求解流程图图 根据步骤，定义分析类型为模型三维实体的温度场仿真，设置初始温度为 800，淬火介质温度为30，载荷设置为30秒，每载荷步为1秒。 步骤3 后处理的基本流程： 基本数据 导出数据 列表显示 图图4-4 ANSYS仿真数据后仿真数据后处处理流程理流程图图 由于齿轮渗碳后淬火的特殊性，整个瞬态过程的不稳定结果，使得我们无法 给出确切的即时温度场，所有的仿真都是对可能持续时间的估计和推断，也就是 说，我们只能仿真出淬火温度场出现的趋势和范围，而无法给出某点的即时温度。 通常淬火是为了强化齿面硬度、增

34、加耐磨性，同时增加齿轮内部应力，产生 应变。淬火增效强化的时间通常在开始淬火的几十秒内，但由于淬火冷却的实际 情况比较复杂，淬火增效强化的时间范围难以把握，为了研究需要，我们选取 05秒、010秒、015秒为可能出现的淬火增效强化范围，逐一进行温度场仿 真。淬火增效进行之后的温度场变化对齿轮齿面性能影响较小，产生的内应力及 应变小于前者，基本上属于典型的热胀冷缩类型，故不作为研究对象。以下仅列 出05秒，010秒和015秒的温度场、应力场仿真分析。 设置载荷步的选择 - 14 - 4.1.3 温度温度场场仿真仿真 齿廓和内圈边界受淬火油作用，从被加热的780820骤冷至与淬火油边界 的平衡温度

35、（不等于平均温度），考虑到液体流动速度，冷却时间和平衡温度需要 根据计算过程设定，轮齿两侧结合边界属于对称边界。这样即可建立相应的温度 场计算模型，通过23次迭代计算，求出轮齿与淬火油边界的平衡温度，然后进 行最后的温度场计算。图4-5、4-6、4-7分别为斜齿圆柱齿轮5秒,10秒和15秒淬火 温度场。随着时间延长，冷却层深入轮齿内部。 图图4-5 斜斜齿圆齿圆柱柱齿轮齿轮5秒淬火温度秒淬火温度场场 - 15 - 图图4-6 斜斜齿圆齿圆柱柱齿轮齿轮10秒淬火温度秒淬火温度场场 图图4-7 斜斜齿圆齿圆柱柱齿轮齿轮15秒淬火温度秒淬火温度场场 由仿真图可以看出齿心与齿面的温差变化，在5秒时，T

36、1=53；在10秒时， - 16 - T2=79；在15秒时，T3=87。随着温度的推移，齿心与齿面的温差逐渐增大， 其热应力也随之增大，所以在淬火刚开始的这个时间段内，是累积热应力的主要 过程。 由于在淬火的过程中，齿面与齿心的温差较大，就整个齿轮而言，尺寸较小 的轮齿部分淬火过程中受作用最为明显，会引起不同程度的尺寸和形状误差，所 以在齿轮的加工过程中，要根据仿真的结果，分析其综合变形量，应用到实际的 生产加工中，为齿轮的下一步工艺加工，取一个合适的叫加工余量。 通过对温度场仿真结果进行分析，我们可以清晰地得出齿轮淬火过程中各 时间段的温度变化，各部分间的温差范围，推算出齿轮的各部分的淬火

37、变形方式 和形变量，同样，根据对仿真结果个时间段的对比分析，可以较快捷准确地确定 齿轮淬火过程中的最佳时间，为齿轮淬火分析提供理论依据。 4.2 淬火淬火应应力分布的仿真力分布的仿真 由于齿轮的不同部分之间的热膨胀系数不匹配，在淬火的急速冷却过程中， 彼此的收缩不一致，从而导致热应力的产生，同时也达到了齿面增效强化的目的。 热应力问题实际又是两个物理场之间的相互作用，故属于耦合场问题。ANSYS提 供了两种分析热应力的方法：直接法和间接法【15】。 分析 1 分析 2 结果 1 结果 2 图图4-8 间间接法分析数据流程接法分析数据流程图图 使用间接法的步骤基本上是先进行热分析，重新进入前处理

38、，然后设置结构 分析，读入热分析的节点温度后再设置参考温度，最后进行求解和后处理。我们 采用后者，先进行热分析获得淬火温度场分布，然后将求得的节点温度作为体载 荷加到结构应力分析中。淬火过程产生不均匀的应力分布，齿廓应力变化明显， 应力集中在齿根部位，图4-9、4-10、4-11分别为5秒、10秒和15秒轮齿增效强化的 应力分布。 - 17 - 图图4-9 斜斜齿圆齿圆柱柱齿轮齿轮5秒淬火秒淬火热应热应力力 - 18 - 图图4-10 斜斜齿圆齿圆柱柱齿轮齿轮10秒淬火秒淬火热应热应力力 图图4-11 斜斜齿圆齿圆柱柱齿轮齿轮15秒淬火秒淬火热应热应力力 通过对图4-9进行分析，可以得到以下几

39、点结论： 1）齿廓应力变化明显，应力集中在齿根部位，且齿根部位应力最大，极大值 为530.4Mpa,齿顶最小,极小值为2.7Mpa; 2）两侧齿面应力较芯部大，形成狭窄齿面张力的集中区域，且由齿根向齿 顶衰减，变化幅度较大； 3）内圈应力集中明显，内圈到齿顶应力变化平稳； 4）其它时间段的轮齿增效强化情况与图4-9类似，随着时间段拉长，应力向 轮齿内部扩散。 第五章第五章 全文全文总结总结 5.1 结论结论 本课题采用了 ANSYS 有限元分析计算软件，对斜齿圆柱齿轮进行温度场 应力场的耦合、温度场、应力分布和应变的仿真计算。在此期间，我根据齿轮参 - 19 - 数，用 SolidWorks

40、三维绘图软件绘制出了斜齿圆柱齿轮的三维造型图，并把数据 导入 ANSYS 软件，建立了斜齿圆柱齿轮淬火过程热力学模型，接着用计算机仿 真了斜齿圆柱齿轮淬火的过程，并简要分析了大模数斜齿圆柱齿轮淬火的变形。 首先根据参数，用 SolidWorks 三维绘图软件进行了斜齿圆柱齿轮的三维造 型，并利用 SolidWorks 软件和 ANSYS 软件的的数据传输功能，将齿轮的三维造 型图导入了 ANSYS 软件，并采用了 ANSYS 有限元分析计算软件，完成了对斜 齿圆柱齿轮淬火过程中温度场和应力场的有限元仿真分析，并通过对仿真结果 的分析，得出了淬火过程中的各项数据。 通过此次对大模数齿轮的淬火过程进行计算机仿真和分析，所揭示的规律 和结果，对以后的深入研究都有一定的参考价值和意义。 5.2 研究不足及展望研究不足及展望 由于自身对 ANSYS 软件的接触时间不长，尚达不到很熟练的运用，所以仿 真出来的结果会存在偏差，加之在进行齿轮淬火过程中温度场和应力场的 ANSYS 有限元仿真过程中，所选的有关参数数据和计算方法，以及所建模型，都 只考虑了主要因素的作用，忽略了次要因素的影响，同时，也没有考虑