基于ANSYS汽车发动机前端齿轮的接触应力分析毕业设计论文.doc

本科毕业设计（论文）论文题目：基于ANSYS汽车发动机前端齿轮的接触应力分析毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作 者 签 名： 日 期： 指导教师签名： 日期： 使用授权说明本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名： 日 期： 学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名： 日期： 年 月 日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权 大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期： 年 月 日导师签名： 日期： 年 月 日注 意 事 项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词 5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订摘要在汽车行业中，有限元法广泛应用于各大汽车总成，包括车架、车身、车桥、离合器、轮胎等零部件以及驾驶室噪声分析，大大提高了汽车设计水平，大幅度提高了汽车的设计水平，减少了试验次数，缩短了开发周期，正在成为设计计算的强有力工具。利用ANSYS进行有限元计算，分析了渐开线斜齿轮啮合最不利加载线位置,并在此位置对齿轮施加载荷，提高齿轮强度分析计算的精度，为齿轮的设计提供更准确的数据。同时，研究不同的过渡圆角对齿根最大接触应力的影响，可以为提高轮齿弯曲强度的表面硬化热处理方法提供定量依据。本文基于CAD/CAE,对发动机前段齿轮做了一下几方面工作:1利用Pro/e软件强大参数建模方法建立渐开线齿轮的三维实体模型。在PRO/E 4.0中，通过绘制渐开线的方法创建渐开线齿轮相对来说比较复杂。2发动机齿轮接触应力的有限元计算 将Pro/e中的一对齿轮以IGES格式文件导入到ANSYS中生成有限元模型，并设置合理的接触对，对齿轮的模型及其边界条件进行适当的简化，计算出齿轮在运动过程中几个极限的接触点，再施加约束和载荷，分析这几个极限情况下，齿轮啮合的接触应力变化。3. 发动机齿轮接触应力的理论值计算 根据机械设计理论，用赫兹应力公式计算出接触应力的理论值，和有限元计算的结果进行比较分析。发现赫兹应力公式的计算结果要比有限元大一点。实际上，滑动的大小和方向、摩擦系数、润滑状态等，都会影响实际接触应力的大小，而赫兹应力公式没有考虑到这些因素。通过以上三方面的工作，为齿轮简单的设计和分析提供一套有效的方法；和传统的计算方法相比，显示出有限元分析方法快速有效、降低成本的优点。关键词：发动机前端齿轮，有限元，接触应力，ANSYSABSTRACTIn the automotive industry, the finite element method is widely used in major automobile assembly, comprising a frame, body, axle, clutch, tires and other parts as well as the cab interior noise analysis, greatly improving the automobile design level, greatly improves the cars design level, reduce test times, shorten the development cycle, is a design calculation of the powerful tool. The use of ANSYS finite element calculation, analysis of involute helical gear meshing the most unfavorable loading line position, and in this position on the wheel applied load, improve the analysis of gear strength calculation accuracy, gear design to provide more accurate data. At the same time, study the different transition fillet on the tooth root of maximum contact stress effect, can improve the tooth bending strength surface hardening heat treatment method to provide quantitative basis.Based on CAD/CAE software, the following aspects are done for the engine front-end gear transmission system: The build of 3d model for engine gearBy the pro/e 4.0 , the 3d model of gear is established .In the pro/e 4.0 , the method of creating standard gear by drawing involute ways is quite complex. The calculation of bending stress by FEM.Import the three dimensional finite element model of involute gear to theSoftware of ansys , and then make a simplification to the model and boundary conditions . Calculate several points of gear contact in the athletic and constraint the gear . Load resistance in the gear and analysis the change of bending stress of gear engagement during the limit condition. The theoretical calculation of bending stress. According to the theory of mechanical design , use the national standard method to calculate the bending stress of the theoretical value. Then compare with the finite element calculation result.Through the work of above three aspects ,a set of effective method for gear design and analysis is provided . Compared with the traditional calculation method, the finite element analysis method has the rapid and effective advantages.Keywords: engine gear , finite element analysis, bending stress, theoretical calculation目录第一章 绪论1第一节 引言1第二节 齿轮接触应力研究历史和现状1第三节 齿轮CAD/CAE的研究现状1第四节 有限元方法概述1第五节 ANSYS的主要功能2第二章 发动机齿轮的三维建模4第一节 几何建模技术4第二节 齿轮三维模型建立5一、 创建渐开线5二、 创建出一个齿廓8三、 创建齿根过渡曲线9四、 计算齿宽并对2D草图进行拉伸9五、 对轮齿进行阵列从而创建出完整的齿轮11第三章 齿轮静态接触应力有限元计算14第一节 齿轮传动的失效形式及设计准则14一、 失效形式14二、 设计准则16第二节 静力学分析基础16一、 结构静力学分析简介16二、 AWE静力学分析基本过程17第三节 三维齿轮模型接触应力分析17一、 模型的简化17二、 导入模型到ANSYS并添加材料信息18三、 网格划分19四、 设置边界条件21五、 分析求解24六、 求解并观察结果26第四节 分析结果概述27一、 总变形云图27二、 应力奇异27三、 齿根圆角半径大小对齿根处应力的影响28第五章 结论与展望38第一节 结论38第二节 展望38致谢39参考文献39第一章 绪论1.1 引言齿轮机构是应用较为广泛的一种传动机构，齿轮齿条传动作为齿轮传动的特殊机构，可以实现直线运动与圆周运动之间转换，由于齿轮齿条传动是线接触机构，其接触应力的大小对于齿轮齿条的使用平稳性和使用寿命有着重要影响。齿轮是传动系统中承受载荷和传动动力的主要零部件，也是最容易出故障的零件之一。据统计，在各种机械故障中，齿轮失效就占总数的6 0%以上，其齿面损坏又是齿轮失效的主要原因之一。因此，工程中需要发大量工作对齿面强度及其应力进行分析。齿轮啮合过程作为一种接触行为, 因涉及接触状态的改变而成为一个复杂的非线性问题。传统的齿轮理论分析是建立在弹性力学基础上的, 对于齿轮的接触强度计算均以两平行圆柱体对压的赫兹公式为基础,在计算过程中存在许多假设,不能准确反映齿轮啮合过程中的应力以及应变分布与变化。相对于理论分析,有限元法则具有直观、准确、快速方便等优点。ANSYS是随着电子计算机的发展而迅速发展起来的一种在计算数学，计算力学和计算工程科学领域最有效的通用有限元分析软件。它是融结构，热，流体，电磁，声学于一体的大型通用有限元商用分析软件。利用ANSYS有限元分析，可以对各种机械零件，构件进行应力，应变，变形，疲劳分析，并对某些复杂系统进行仿真，实现虚拟的设计，从而大大节省人力，财力和物力。由于其方便性、实用性和有效性，ANSYS软件在各个领域，特别是机械工程当中得到了广泛的应用。1.2 课题研究背景随着计算技术的迅速发展与广泛应用，以有限元法为代表的数值计算方法为齿轮应力和变形分析提供了一种方便、可靠的研究方法目前齿轮工程中实用的数值解法主要有三种：有限差分法(FDM)、边界元法（BEM)和有限元法(FEM)。在数值计算方法中最引人注目的是有限元法。有限元法用于齿根应力分析大约起始于二十世纪六十年代末、七十年代初，此后迅速发展，国外不少研究人员如Chabert、Wilcox、户部、Chang、Bibel等都进行过这方面的研究工作因此，在用有限元方法对直齿轮的齿根应力进行分析时，都把它简化为力学中的平面应变问题。为了进行齿面接触强度计算，分析齿面失效和润滑状态，必须分析齿面的接触应力。经典的齿面接触应力计算公式是建立在弹性力学基础上，而对于齿轮的接触强度计算均以两平行圆柱体对压的赫兹公式为基础。由于齿轮副啮合齿面的几何形状十分复杂，采用上面的方法准确计算轮齿应力和载荷分配等问题非常困难甚至无法实现。随着计算机的普及，齿轮接触问题的数值解法获得了越来越广泛的应用。数值解法可以求解复杂的齿面接触问题，但不能给出一般性的函数关系。在工程应用上数值解法具有很大的实用价值，己经取得了很多重要成果，例如有限元法、边界元法、有限差分法以及与数值方法相配合的各种变分法、实变函数法、泛函分析法等。在所有这些方法中，有限元法的应用最为广泛，可以求解边界条件、几何形状和载荷方式复杂的工程接触问题。轮副在工作时，在内部和外部激励下将发生机械振动。振动系统的固有特性，一般包括固有频率和主振型，它是系统的动态特性之一，同时也可以作为其它动力学分析的起点，对系统的动态响应、动载荷的产生与传递以及系统振动的形式等都具有重要的影响。在进行结构设计时，使激振力的频率与系统的固有频率错开，可以有效的避免共振的发生。然而，在齿轮的设计阶段，往往很难得到齿轮固有特性的实验数据，只能通过理论计算得到进行动力学分析的参数，目前最好的方法是有限元分析法。1.3 PRO/E概述Pro/Engineer操作软件是美国参数技术公司（PTC)旗下的CAD/CAM/CAE一体化的三维软件。Pro/Engineer软件以参数化著称，是参数化技术的最早应用者，在目前的三维造型软件领域中占有着重要地位，Pro/Engineer作为当今世界机械CAD/CAE/CAM领域的新标准而得到业界的认可和推广。是现今主流的CAD/CAM/CAE软件之一，特别是在国内产品设计领域占据重要位置。Pro/E第一个提出了参数化设计的概念，并且采用了单一数据库来解决特征的相关性问题。另外，它采用模块化方式，用户可以根据自身的需要进行选择，而不必安装所有模块。Pro/E的基于特征方式，能够将设计至生产全过程集成到一起，实现并行工程设计。它不但可以应用于工作站，而且也可以应用到单机上。Pro/E采用了模块方式，可以分别进行草图绘制、零件制作、装配设计、钣金设计、加工处理等，保证用户可以按照自己的需要进行选择使用。 1 参数化设计相对于产品而言，我们可以把它看成几何模型，而无论多么复杂的几何模型，都可以分解成有限数量的构成特征，而每一种构成特征，都可以用有限的参数完全约束，这就是参数化的基本概念。 2 基于特征建模 Pro/E是基于特征的实体模型化系统，工程设计人员采用具有智能特性的基于特征的功能去生成模型，如腔、壳、倒角及圆角，您可以随意勾画草图，轻易改变模型。这一功能特性给工程设计者提供了在设计上从未有过的简易和灵活。 3 单一数据库（全相关） Pro/Engineer是建立在统一基层上的数据库上，不像一些传统的CAD/CAM系统建立在多个数据库上。所谓单一数据库，就是工程中的资料全部来自一个库，使得每一个独立用户在为一件产品造型而工作，不管他是哪一个部门的。换言之，在整个设计过程的任何一处发生改动，亦可以前后反应在整个设计过程的相关环节上。例如，一旦工程详图有改变，NC（数控）工具路径也会自动更新；组装工程图如有任何变动，也完全同样反应在整个三维模型上。这种独特的数据结构与工程设计的完整的结合，使得一件产品的设计结合起来。这一优点，使得设计更优化，成品质量更高，产品能更好地推向市场，价格也更便宜。1.4 ANSYS概述ANSYS有限元软件包是一个多用途的有限元法计算机设计程序，目前，有限元法从它最初应用的固体力学领域，已经推广到温度场、流体场、电磁场、声场等其他连续介质领域。在固体力学领域，有限元法不仅可以用于线性静力分析，也可以用于动态分析，还可以用于非线性、热应力、接触、蠕变、断裂、加工模拟、碰撞模拟等特殊问题的研究。软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。前处理模块前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型。ANSYS的前处理模块主要有两部分内容：实体建模和网格划分。分析计算模块分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力。后处理模块后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。软件提供了200种以上的单元类型，用来模拟工程中的各种结构和材料。ANSYS软件提供的分析类型如下：结构静力分析用来求解外载荷引起的位移、应力和力。静力分析很适合求解惯性和阻尼对结构的影响并不显著的问题。ANSYS程序中的静力分析不仅可以进行线性分析，而且也可以进行非线性分析，如塑性、蠕变、膨胀、大变形、大应变及接触分析。结构动力学分析结构动力学分析用来求解随时间变化的载荷对结构或部件的影响。与静力分析不同，动力分析要考虑随时间变化的力载荷以及它对阻尼和惯性的影响。ANSYS可进行的结构动力学分析类型包括：瞬态动力学分析、模态分析、谐波响应分析及随机振动响应分析。结构非线性分析结构非线性导致结构或部件的响应随外载荷不成比例变化。ANSYS程序可求解静态和瞬态非线性问题，包括材料非线性、几何非线性和单元非线性三种。动力学分析结构动力学分析研究结构在动载荷作用的响应（如位移、应力、加速度等得时间历程），以确定结构的承载能力的动力特性等。ANSYS程序可以分析大型三维柔体运动。当运动的积累影响起主要作用时，可使用这些功能分析复杂结构在空间中的运动特性，并确定结构中由此产生的应力、应变和变形。热分析程序可处理热传递的三种基本类型：传导、对流和辐射。热传递的三种类型均可进行稳态和瞬态、线性和非线性分析。热分析还具有可以模拟材料固化和熔解过程的相变分析能力以及模拟热与结构应力之间的热结构耦合分析能力。电磁场分析主要用于电磁场问题的分析，如电感、电容、磁通量密度、涡流、电场分布、磁力线分布、力、运动效应、电路和能量损失等。还可用于螺线管、调节器、发电机、变换器、磁体、加速器、电解槽及无损检测装置等的设计和分析领域。流体动力学分析ANSYS流体单元能进行流体动力学分析，分析类型可以为瞬态或稳态。分析结果可以是每个节点的压力和通过每个单元的流率。并且可以利用后处理功能产生压力、流率和温度分布的图形显示。另外，还可以使用三维表面效应单元和热流管单元模拟结构的流体绕流并包括对流换热效应。声场分析ANSYS把声学归为流体，程序的声学功能用来研究在含有流体的介质中声波的传播，或分析浸在流体中的固体结构的动态特性。这些功能可用来确定音响话筒的频率响应，研究音乐大厅的声场强度分布，或预测水对振动船体的阻尼效应。压电分析压电效应分析是一种结构电场耦合分析，给压电材料加电压会产生位移，反之使压电材料振动则产生电压，一个典型的压电分析的应用是压力换能器。ANSYS压电分析用于分析二维或三维结构对AC（交流）、DC（直流）或任意随时间变化的电流或机械载荷的响应。这种分析类型可用于换热器、振荡器、谐振器、麦克风等部件及其它电子设备的结构动态性能分析。可进行四种类型的分析：静态分析、模态分析、谐波响应分析、瞬态响应分析。1.4 ANSYS的接触能力接触问题是一种高度非线性行为，需要较大的计算资源，分为两种基本类型：刚体柔体的接触，柔体柔体的接触，在刚体柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体柔体的接触，许多金属成形问题归为此类接触，另一类，柔体柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。ANSYS支持三种接触方式：点点，点面，面面，每种接触方式使用的接触单元只适用于某一类问题。为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元。在齿轮的有限元分析问题上的接触单元为面面的接触单元，故只介绍这一类的接触单元。 ANSYS支持刚体柔体的面面的接触单元，刚性面被当作“目标”面，分别用Targe169和Targe170来模拟2D和3D的“目标”面，柔性体的表面被当作“接触”面，用Conta171,Conta172,Conta173,Conta174来模拟。一个目标单元和一个接单元叫作一个“接触对”程序通过一个共享的实常号来识别“接触对”，为了建立一个“接触对”给目标单元和接触单元指定相同的实常的号。 与点面接触单元相比，面面接触单元有以下几项优点： 1.支持低阶和高阶单元； 2.支持有大滑动和摩擦的大变形，协调刚度阵计算，单元提法不对称刚度阵 的选项； 3.提供工程目的采用的更好的接触结果，例如法向压力和摩擦应力； 4.没有刚体表面形状的限制，刚体表面的光滑性不是必须允许有自然的或网格离散引起的表面不连续； 5.与点面接触单元比，需要较多的接触单元，因而造成需要较小的磁盘空间和CPU时间； 6.允许多种建模控制，例如：绑定接触、渐变初始渗透、目标面自动移动到补始接触 、平移接触面（老虎梁和单元的厚度）、支持死活单元。 使用这些单元，能模拟直线（面）和曲线（面），通常用简单的几何形状例如圆、抛物线、球、圆锥、圆柱采模拟曲面，更复杂的刚体形状能使用特殊的前处理技巧来建模。 接触单元被约束住，不能穿透目标面，然而目标单元可以穿透接触面。指定目标面和接触面时应考虑以下准则： 1、当凸面和平面或凹面接触时，应指定平面或凹面作为目标面； 2、如果一个面上的网格较细，而另一个面上的网格较粗，应指定细网格所在面作为接触面，粗网格所在面作为目标面； 3、如果两个面的刚度不同时，应当较硬的面为目标面，较软的面作为接触面； 4、如果一个面上的基础单元是高阶单元，而另一个面上的基础单元是低阶单元，应将基础单元为高阶单元的面作为接触单元； 5、如果两个面的大小明显不同，应将大面作为目标面。第二章 发动机齿轮的三维建模2.1 大小齿轮的基本参数和尺寸大齿轮材料泊松比齿宽模数齿数分度圆压力角齿顶高系数顶隙系数20CrMnTi=0.340433=2010.25小齿轮材料泊松比齿宽模数齿数分度圆压力角齿顶高系数顶隙系数20CrMnTi=0.350421=2010.252.2 建立齿轮底板三维模型1. 在Pro/e中绘制如图2-1所示底板草图图2-1.齿轮底板二维草图2.采用“拉伸”命令建立如图2-2所示齿轮底板三维零件图图2-2.齿轮底板三维零件图2.3 建立齿轮固定轴三维模型在Pro/e中分别绘制如图固定轴零件图，两个轴的半径分别与底板的两个孔半径相同2.4 建立大齿轮三维模型1. 定义齿轮参数选择菜单“工具”，点击“参数”命令，此时系统弹出“参数”窗口。分别添加如图2-3所示4个参数，其中m为模数，z为齿数，b为齿宽，pa为压力角。图2-3.定义新参数2. 创建旋转特征选择菜单“工具”，点击“关系”命令，打开“关系”窗口。输入草图尺寸和设置的参数的关系，如图2-4所示。并采用“旋转”命令建立如图2-5所示的大齿轮的旋转特征在对话框中输入关系式为：Sd0=bSd2=0.25*bSd3=0.5*m*z+m图2-4.设置关系式图2-5.创建的旋转特征3. 草绘曲线，画出四个圆选择菜单“工具”，点击“关系”命令，打开“关系”窗口。输入草图尺寸和设置的参数的关系，如图2-6所示。在对话框中输入关系式为：Sd0=m*(z=2)Sd1=m*zSd2=m*z*cos(pa)Sd3=m*z-2.5*mDB=sd2图2-6.绘制圆4. 创建渐开线（1） 单击（插入基准曲线）按钮，弹出“曲线选项”菜单（2） 在“曲线选项”菜单中，选择“从方程”，点击“完成”，如图2-7所示。图2-7.“曲线选项”菜单（3） 在模型树中选择PRT_CSYS_DEF基准坐标系，弹出如图所示“设置坐标类型”菜单，点击“笛卡尔”，如图2-8所示图2-8.“设置坐标类型”菜单（4） 弹出记事本编辑器，输入函数方程，如图2-9所示，保存记事本文件，并退出图2-9.定义渐开线方程（5） 完成渐开线的创建，如图2-10所示。图2-10.完成一条渐开线5. 镜像渐开线，如图2-11所示。图2-11.镜像结果6. 绘制齿槽草图，如图2-12所示。图2-12.绘制渐开线草图7. 拉伸并阵列齿槽，如图2-13所示。图2-13.拉伸阵列齿槽8. 完成齿轮的三维建模，如图2-14所示。图2-14.完成的大齿轮2.5 建立小齿轮三维模型按照大齿轮建模的步骤建立小齿轮三维模型，如图2-15所示。图2-15.完成的小齿轮2.6 装配组合体1. 装配齿轮底板和两根固定轴，如图2-16所示。图2-16.装配底板和轴2. 装配齿轮，如图2-17所示。图2-17.装配齿轮3. 设置两个齿轮正好啮合，如图2-18所示。图2-18.齿轮啮合4. 导出IGES打开菜单“文件”，点击“保存副本”，设置类型为IGES，名称为jwx1_asm，如图2-19所示。图2-19.保存副本IGES第三章 齿轮静态接触应力有限元计算3.1 Pro/e模型导入ANSYS（1）选择Utility MenuFileImportIGES命令，弹出import IGES file对话框，单击import IGES file对话框上的【OK】按钮。单击【Browse】按钮，打开jwx1_asm.igs所在的文件夹，选中文件名为jwx1_asm.igs的文件，单击import IGES file对话框上的【OK】按钮。（2）选择Utility MenuPoltVolumes命令，ANSYS显示窗口将显示如图3-1所示齿轮组件模型显示。图3-1.齿轮组件模型显示3.2 设置计算类型（1）选择Main MenuPrepferences命令，出现Prepferences for GUI Filering对话框，选择Structural，单击【OK】按钮。如图3-2所示。图3-2.设置计算类型3.3 定义单元类型（1）选择Main MenuPreprocessorElementAddEditDelete命令，出现Element Types对话框，单击【Add】按钮，出现Library of Element Types对话框。（2）在Library of Element Types对话框复选框中选择Structural Solid ，Tet 20node 187（选择用于金属材料分析的三维六面体等参元SOLID187,该单元具有10个节点，适用于带中间节点的复杂四面体结构），如图3-3所示“单元类型列表”对话框，单击【OK】按钮，关闭该对话框。图3-3.“单元类型列表”对话框3.3 定义材料性能参数（1）选择Main MenuPreprocessorMaterial PropsMaterial Models命令，出现Define Material Model Behavior对话框。（2）在Material Models Available一栏中依次单击Structural、Linear、Elastic、Isotropic选项，出现Linear Isotropic Propeties for Material Number 1对话框，在EX输入栏中输入2.11E11，在PRXY输入栏中输入0.3查阅资料可知，汽车发动机前端齿轮采用合金钢制造，在本文中，齿轮以20CrMnTi为材料，其它各零件均采用45号钢制造,故由汽车发动机前端齿轮的材料的弹性模量与泊松比分别为211GPa、0.3。如图3-4所示“输入材料弹性模量和泊松比”对话框，单击【OK】按钮，关闭该对话框。图3-4. “输入材料弹性模量和泊松比”对话框3.4 网格划分 (1) 选择Main MenuPreprocessorMeshingMesh ToolGlobalSize命令，弹出Global Size对话框，在SIZE Element egde length输入栏中输入8，设置单元尺寸为8，如图3-5所示“定义单元尺寸”对话框，单击【OK】按钮，关闭该对话框。(2) 单击【Mesh】按钮，弹出Mesh Volumes拾取菜单，在显示窗口中选中大齿轮模型，单击【OK】按钮，关闭该对话框，对大齿轮进行网格划分。图3-5. “定义单元尺寸”对话框(3) 重复大齿轮网格划分步骤，设置单元尺寸为5，ANSYS显示窗口将显示如图3-6所示两个啮合齿轮模型网格划分结果显示。图3-6. 两个啮合齿轮模型网格划分结果显示对于给定的材料的单元类型，可对模型进行划分单元，并可以根据实际计算需要，设置网格划分大小。采用ANSYS智能网格划分技术，根据技算的需要，单元划分得比较小，两个齿轮共划分了25426个单元，43817个节点。3.6 创建接触对（1）选择Main MenuPreprocessorModelingCreateContact Pair命令，弹出Contact Manager对话框。如图3-7所示图3-7.“Contact Manager”对话框（2）创建第一个接触面单击左上角图标（Contact Wizard），弹出Contact Wizard对话框，单击对话框中的“Areas”单选按钮，指定接触目标表面为面，如图3-8所示“添加接触”对话框，然后单击【Pick Target】按钮。选择具体目标面，弹出Select Area for Target拾取菜单，用鼠标在屏幕上拾取齿轮啮合的第一个接触面，如图3-9所示，单击【OK】按钮，将其关闭。图3-8.选择第一个接触面第一个面对话框图3-9.选择第一个接触面的显示（2）创建第二个接触面这时，Contact Wizard对话框中【next】按钮将被激活，单击【next】按钮进入下一步，将弹出选中接触面的对话框。单击对话框中的“Areas”单选按钮，指定接触目标表面为面，如图4-7所示“添加接触”对话框，然后单击【Pick Contact】按钮，如图3-10所示。选择具体目标面，弹出Select Area for Target拾取菜单，如图3-11所示，用鼠标在屏幕上单击拾取齿轮啮合的第二个接触面，单击【OK】按钮，将其关闭。这时，Contact Wizard对话框中【next】按钮将被激活，单击【next】按钮进入下一步，对接触面属性进行设置。 图3-10.选择第二个接触面第一个面对话框 图3-11.选择第二个接触面的显示在对话框中单击“Include initial penetration”按钮，将其选中，使分析包括初始渗透，单击Material ID下拉框中的1，指定接触材料属性为定义的一号材料。在Coeffieient of Friction中输入0.2，如图3-12所示“设置接触对属性”对话框。单击【