齿轮有限元分析 - 图文

基于ANSYS的齿轮传动有限元分析和优化

摘要

ANSYS是随着电子计算机的发展而迅速发展起来的一种在计算数学，计算力学和计算工程科学领域最有效的通用有限元分析软件。它是融结构，热，流体，电磁，声学于一体的大型通用有限元商用分析软件。利用ANSYS有限元分析，可以对各种机械零件，构件进行应力，应变，变形，疲劳分析，并对某些复杂系统进行仿真，实现虚拟的设计，从而大大节省人力，财力和物力。由于其方便性、实用性和有效性，ANSYS软件在各个领域，特别是机械工程当中得到了广泛的应用。

齿轮是机械中常用的一种零件，其在工作的过程中会产生应力，应变和变形，为保证其正常工作需要对齿轮的轮齿和整体受力进行分析，保证其刚度和强度的要求。本论文采用ANSYS软件对齿轮进行静力学分析和优化实现对齿轮的虚拟设计。

齿轮是最重要的零件之一。它具有功率范围大，传动效率高，传动比正确，使用寿命长等特点，但从零件失效的情况来看，齿轮也是最容易出故障的零件之一。据统计，在各种机械故障中，齿轮失效就占故障总数的60%以上。其中轮齿的折断又是齿轮失效的主要原因之一。

齿轮啮合过程作为一种接触行为, 因涉及接触状态的改变而成为一个复杂的非线性问题。传统的齿轮理论分析是建立在弹性力学基础上的, 对于齿轮的接触强度计算均以两平行圆柱体对压的赫兹公式为基础,在计算过程中存在许多假设,不能准确反映齿轮啮合过程中的应力以及应变分布与变化。相对于理论分析,有限元法则具有直观、准确、快速方便等优点。

齿廓曲面是渐开线曲面，所以建模的难点和关键在于如何确定精确的渐开线。通过PDL命令流直接在ANSYS中创建标准直齿圆柱齿轮，学习应用ANSYS软件进行零件的几何建模和网格划分，并进行静力加载和求解，对求解的结果进行查看，分析和优化。

关键词：ANSYS；有限元；齿轮；CAE

Gear Transmission Of Finite Element Analysis And

Optimization

Abstract

ANSYS is along with the rapid development of electronic computers and developed a computational mathematics, computational mechanics and engineering science, the most effective general finite element analysis software. It is hot, the fluid, structure, electromagnetic, acoustics integration in the universal finite element analysis software for commercial. Using the ANSYS finite element analysis, all kinds of machine parts, can carry out stress, strain and structural deformation, fatigue analysis of some complex system, and the simulation, the design and realization of virtual human, to save money and material. Due to its convenience, practicability and validity, ANSYS software, especially in the field of mechanical engineering has been widely used.

Gear is commonly used in machinery, a part of the work in process of stress, strain and can produce deformation, so as to ensure the normal working of gear teeth and to overall analysis, ensure the stiffness and strength. This thesis of ANSYS software of gear static analysis and optimization of virtual design of gear.

Gear is one of the most important parts. It has big power range, high transmission efficiency and transmission ratio correctly, long using life, etc, but from the failure parts, gear is the most vulnerable parts of the fault. According to statistics, in all kinds of mechanical failure, gear failure is accounted for 60% of the total failure. One of the broken tooth gear is one of the main reasons.

Gear meshing process as a contact, because involves contact state changes a complex nonlinear problems. The traditional theory of gear analysis was based on the basis of elastic mechanics, the contact strength for gear with two parallel computation formula of the cylinder pressure, based in Hertz calculation process in many assumptions, was not accurate in reflecting gear meshing process of stress and strain distribution and change. Relative to the theoretical analysis, finite element method, the principle is convenient and fast accurate, etc. Involute tooth profile surface is curved, so the difficulties and modeling key lies in how to determine the precise involute. Through PDL coupler, single mode WDMS directly in order to create ANSYS flow standard spur gears, study on parts of ANSYS software, and the meshing geometric modeling and static load and the solving of solving the check, analysis and optimization.

Key words: ANSYS; Finite element; Gear; CAE

目录

1 绪论 .................................................................................................................................... - 1 -

1.1有限元概述 ................................................................................................................................. - 1 - 1.2选题背景 ..................................................................................................................................... - 3 - 1.3 本文主要工作 ............................................................................................................................ - 3 -

2 ANSYS准备工作 ................................................................................................................. - 4 -

2.1 ANSYS安装与启动 ..................................................................................................................... - 4 -

2.1.1 许可证服务器安装 ........................................................................................................ - 4 - 2.1.2 主程序安装 .................................................................................................................... - 5 - 2.1.3 启动许可证服务器 ........................................................................................................ - 7 - 2.1.4 ANSYS启动与配置 ......................................................................................................... - 8 - 2.2 设置工作目录 .......................................................................................................................... - 10 - 2.3 指定作业名与分析标题 .......................................................................................................... - 10 -

2.3.1 指定作业名 .................................................................................................................. - 10 - 2.3.2 指定分析标题 ...............................................................................................................- 11 - 2.4 定义图形界面过滤参数 ...................................................................................................- 11 - 2.5 ANSYS单位制 ................................................................................................................... - 12 - 2.6 选取和定义单元 .............................................................................................................. - 13 -

3 在ANSYS中建立齿轮分析模型 ...................................................................................... - 15 -

3.1 几何模型的建立 ...................................................................................................................... - 15 -

3.1.1 大小齿轮的具体基本参数和尺寸 .............................................................................. - 15 - 3.1.2 渐开线的生成原理 ...................................................................................................... - 16 - 3.1.3 创建渐开线曲线 .......................................................................................................... - 16 - 3.1.4 生成齿根过渡曲线 ...................................................................................................... - 18 - 3.1.5 生成完整齿廓线 .......................................................................................................... - 18 - 3.1.6 生成完整齿轮的面 ...................................................................................................... - 19 - 3.1.7 生成大齿轮 .................................................................................................................. - 20 - 3.1.8 生成两齿轮的啮合图 .................................................................................................. - 22 - 3.2 几何模型的网格划分 .............................................................................................................. - 22 -

3.2.1 定义单元属性 .............................................................................................................. - 23 - 3.2.2 定义网格生成控制并生成网格 .................................................................................. - 24 -

4 ANSYS静力加载与静力求解 ........................................................................................... - 27 -

4.1创建接触对 ............................................................................................................................... - 27 - 4.2 ANSYS施加边界条件和加载 ................................................................................................... - 29 - 4.3 ANSYS求解 ............................................................................................................................... - 31 -

I

5 求解结果的分析和优化 .................................................................................................. - 32 -

5.1 求解结果查看 .......................................................................................................................... - 32 - 5.2 结果分析及结论 ...................................................................................................................... - 34 - 5.3 对齿轮的优化 .......................................................................................................................... - 35 -

6 全文总结与展望 .............................................................................................................. - 36 -

6.1 全文总结 .................................................................................................................................. - 36 - 6.2 工作展望 .................................................................................................................................. - 36 -

参考文献 .............................................................................................................................. - 37 - 附录 ...................................................................................................................................... - 38 - 致 谢 .................................................................................................................................... - 47 -

II

绪论

1 绪论

1.1有限元概述

有限元是随着电子计算机的发展而迅速发展起来的一种现代计算方法。1960年，克拉夫（Clough）在他的一篇论文“平面分析的有限元法（The Finite Element Method in Plane Stress Analysis）”中最先引入了有限元（Finite Element）这一术语。这一方法是结构分析专家把杆件结构力学中的位移法推广到求解连续体介质力学问题（当时是解决飞机结构应力分析）而提出来的。这一方法的提出，引起广泛的关注，吸引了众多力学，数学方面的专家学者对此进行研究。

有限元法之所以能在1960年立刻获得成功，一是Clough从结构力学方法推导的刚度矩阵易于为广大工程师接受，而有限元法最初也被称为矩阵近似方法；二是在于这个方法所包含的大量数值运算，而这可以由新发展起来的数字计算机来完成。

在20世纪70，80年代，许多学者研究和推导出了许多精确，更高效的单元，在单元形状，单元节点和插值函的类型等方面都得到了长足的发展。20世纪70年代，等参元的提出为研发出新的单元开辟了新的途径，推动了有限元的发展。经过近几十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。

现在有限元方法的发展趋势是集成化、通用化、输入智能化和结构输出可视化。所谓集成化是一个有限元程序包往往包括了各种各样的单元（即单元库），并包括了许多材料的本构关系（即材料库），使用者可以根据需要选择和组合；通用化是一个通用程序同时又解决静力分析、动力分析、热传导、电场等各种问题的模块；输入智能化、图形化是计算机辅助输入，只要输入轮廓边界的关键点及计算所需节点数和单元类型，即可自动进行单元网格划分，并且其结果以图形方法表达出来。这样可以快捷，直观且易于发现错误而及时改正；输出结果可视化是计算所得的应力场、位移场、流态场等均可用多方位，多层次的图形或图像表示出来，非常直观，便于分析判断，有些学者称之为仿真或数值分析。

有限元法分析计算的思路和做法可归纳如下： （1） 物体离散化

将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用

- 1 -

基于ANSYS齿轮传动有限元分析及优化

有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。

（2） 单元特性分析 A、选择位移模式

在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。

当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法我们就将位移表示为坐标变量的简单函数。

B、分析单元的力学性质

根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。

C、计算等效节点力

物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上得力。

（3） 单元组集

利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程

（4） 求解未知节点位移

解有限元方程式得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。

通过上述分析，可以看出，有限单元法的基本思想是\一分一合\，分是为了就进行单元分析，合则为了对整体结构进行综合分析。 有限元的发展概况

1943年 courant在论文中取定义在三角形域上分片连续函数，利用最小势能原理研究St.Venant的扭转问题。

1960年 clough的平面弹性论文中用“有限元法”这个名称。

1965年 冯康发表了论文“基于变分原理的差分格式”，这篇论文是国际学术界承认我国独立发展有限元方法的主要依据。

1970年 随着计算机和软件的发展，有限元发展起来。

- 2 -

绪论

涉及的内容：有限元所依据的理论，单元的划分原则，形状函数的选取及协调性。 有限元法涉及：数值计算方法及其误差、收敛性和稳定性。 应用范围：固体力学、流体力学、热传导、电磁学、声学、生物力学

求解的情况：杆、梁、板、壳、块体等各类单元构成的弹性（线性和非线性）、弹塑性或塑性问题（包括静力和动力问题）。能求解各类场分布问题（流体场、温度场、电磁场等的稳态和瞬态问题），水流管路、电路、润滑、噪声以及固体、流体、温度相互作用的问题。 1.2选题背景

齿轮是最重要的零件之一。它具有功率范围大，传动效率高，传动比正确，使用寿命长等特点，但从零件失效的情况来看，齿轮也是最容易出故障的零件之一。据统计，在各种机械故障中，齿轮失效就占故障总数的60%以上。其中轮齿的折断又是齿轮失效的主要原因之一。

齿轮啮合过程作为一种接触行为, 因涉及接触状态的改变而成为一个复杂的非线性问题。传统的齿轮理论分析是建立在弹性力学基础上的, 对于齿轮的接触强度计算均以两平行圆柱体对压的赫兹公式为基础,在计算过程中存在许多假设,不能准确反映齿轮啮合过程中的应力以及应变分布与变化。相对于理论分析,有限元法则具有直观、准确、快速方便等优点。

ANSYS 是一个融结构、热、流体、电、磁、声学于一体的大型通用有限元软件。作为目前最流行的有限元软件之一,它具备功能强大、兼容性好、使用方便、计算速度快等优点,成为工程师们开发设计的首选,广泛应用于一般工业及科学研究领域。而在机械结构系统中,主要在于分析机械结构系统受到负载后产生的反应,如位移、应力、变形等,根据该反应判断是否符合设计要求。 1.3 本文主要工作

1. 利用有限元软件建立齿轮分析模型 2. 利用ANSYS对分析模型进行静力分析 3. 对有限元分析结果进行分析和优化

- 3 -

基于ANSYS齿轮传动有限元分析及优化

2 ANSYS准备工作

2.1 ANSYS安装与启动

ANSYS的许可证服务器（License Server）需要已经配置好的TCP/IP网络协议，后面

的安装过程将假定已经正确配置TCP/IP网络协议。下面介绍WindowsXP-32bit系统下ANSYS11.0的安装。ANSYS11.0的安装包含两部分安装：许可证服务器和主程序的安装和启动。

2.1.1 许可证服务器安装

首先对ANSYS 11.0许可证服务器的安装过程进行介绍和说明。

(1) 在ANSYS11.0安装光盘中找到LicenseSetup子文件夹，双击安装程序setup.exe，开

始安装许可证服务器。

(2) 在开始安装后会出现如图2.1和2.2所示的确认对话框。

图2.1 确认本机是否是服务器机器对话框

图2.2 确认是否有许可证文件对话框

(3) 单击【是】按钮，进入下一步安装，弹出如图2.3所示的许可证服务器安装进度条。

图2.3 许可证服务器安装进度条

(4) 许可证服务器安装进度条结束之后，会弹出一个【ANSYS FLEXlm file】对话框，如

图2.4所示。

- 4 -

ANSY准备工作

(5) 选择Browse for the location of an existing license file（指定许可证文件位置）选项，

然后单击【Next】按钮，打开Select the license file （选择许可证文件）对话框，如图2.5所示。

图2.4 浏览/创建许可证文件对话框 图2.5 选择许可证文件对话框

(6) 在该对话框中找到License.dat文件后单击【打开】按钮，继续安装。

(7) 安装完成后会弹出如图2.6所示的安装完成对话框，单击【Finish】按钮，完成许可

证服务器的安装。

图2.6 安装完成对话框

2.1.2 主程序安装

完成了许可证服务器安装后，接下来进行ANSYS 11.0主程序的安装。 (1) 在ANSYS 11.0安装光盘中找到双击安装程序AutoExec.exe，开始安装主程序。 (2) 安装程序启动后，弹出如图2.7所示的ANSYS 11.0的安装界面。

(3) 单击【Install Products】按钮，弹出如图2.8所示的模块安装选择界面。在界面的右

侧一栏有三项内容，第一项是安装ANSYS主程序；第二项是浏览光盘内容；第三项是浏览ANSYS官方网站。

- 5 -

基于ANSYS齿轮传动有限元分析及优化

图2.7 ANSYS安装界面图 图2.8 模块安装选择界面

选择安装ANSYS的常规模块，单击【ANSYS Products】，显示如图2.9所示的语言选择对话框。

图2.9 语言选择对话框

(4) 对话框中选择下拉菜单中的【English】，单击【OK】按钮，弹出如图2.10欢迎对话

框。

(5) 单击【Next】按钮，进入下一步，弹出如图2.11所示的ANSYS许可协议对话框。

图2.10 ANSYS欢迎对话框 图2.11 许可协议对话框

(6) 在许可协议对话框中选择AGREE后，单击【Next】按钮，弹出如图2.12所示的ANSYS

安装路径对话框。通过点击【Change】按钮，可以改变ANSYS默认的安装路径。 (7) 单击【Next】按钮，弹出如图2.13所示的安装配置特性对话框，通过点击复选框，

可以选择与ANSYS对接的各种CAD和CAE模块。

- 6 -

在ANSYS中建立齿轮分析模型

m=2.5 ！齿轮模数 z2=34 ！小齿轮齿数 pi=acos(-1) ！定义π

alfa2=20/180*pi ！定义齿轮的标准压力角 ha=1 ！定义齿顶高系数 c=0.25 ！定义顶隙系数 d2=m*z2 ！定义分度圆直径 alfa_ding=

acos(z2*cos(alfa2)/(2+z2)) ！定义齿顶压力角 sita_ding=

(tan(alfa_ding)-alfa_ding)*180/pi

sita=tan(alfa2)-alfa2 ！定义齿顶展角 *afun,rad ！角度转化为弧度 db2=m*z2*cos(alfa2) ！定义基圆直径 rb2=db2/2 ！定义基圆半径 df2=d2-2*(ha+c)*m ！定义齿根圆直径 rf2=df2/2 ！定义齿根圆半径 da2=d2+2*ha*m ！定义齿顶圆直径 ra2=da2/2 ！定义齿顶圆半径 zl=alfa_ding/50 ！定义压力角份数 *dim,alfa_12,array,50,1 ！定义50×1数组alfa_12 *dim,sita2,array,49,1 !定义49×1数组sita2 *dim,r2,array,49,1 !定义49×1数组r2 csys,1 ！当前坐标系为柱坐标 alfa_12(1,1)=0 ！数组数alfa_12(1,1)=0 /prep7 ！进入前处理器 k,1000,rb2,0 ！生成关键点1000 *do,j,1,49,1 ！循环生成其他关键点 alfa_12(j+1,1)=alfa_12(j,1)+zl

sita2(j,1)=(tan(alfa_12(j+1,1))-alfa_12(j+1,1))*180/pi r2(j,1)=rb2/cos(alfa_12(j+1,1)) /prep7

k,j+1000,r2(j,1),sita2(j,1)

*enddo ！循环结束生成49个关键点 k,1050,ra2,sita_ding ！生成齿顶圆上的关键点

- 17 -

基于ANSYS齿轮传动有限元分析及优化

bsplin,all ！用样条曲线连接所有点

在这个曲线中先创建51个点，再把所有的点连接起来就形成一条完整的渐开线曲线（图3.3、图3.4）。

图3.3 渐开线关键点 图3.4 渐开线

3.1.4 生成齿根过渡曲线

齿根过渡曲线原本比较复杂，不仅与齿轮工作的参数有关而且还与刀具齿顶形状等一系列参数有关。但在这里简化处理（参照参考文献[2]），对于齿根圆直径小于基圆直径的齿轮做渐开线的切线与基圆相交再利用倒角工具对其倒角，倒角的半径r=0.38*m。如果齿根圆的直径大于基圆直径直接倒角，倒角的半径r=0.38*m。生成的齿根过渡曲线如图3.5所示。

图3.5 齿根过渡曲线

3.1.5 生成完整齿廓线

通过镜像操作Main>Preprocessor>Modeling>Reflect>Lines对图3.5所示齿根过渡曲线进行镜像，在镜像操作中注意要将坐标系旋转θ+90/z2°见图3.7(见参考文献[1])，再用半径为Ra2的圆弧线连接齿顶两点生成完整的齿廓线，如图3.6所示。

图3.6 完整齿廓线

- 18 -

在ANSYS中建立齿轮分析模型

图3.7 镜像原理图

图中Rf—————齿根圆 R —————分度圆 Z —————齿数 ?k—————分度圆展角 3.1.6 生成完整齿轮的面

将生成的齿廓线进行整列，先改变当前坐标为柱坐标系，选择图3.6所示的完整齿廓线再通过阵列操作Main>Preprocessor>Modeling>Copy>Lines，弹出如图3.10所示对话框，在对话框中填入齿数和相应的角度如图3.11所示，用半径为rf2的曲线连接两齿间的间隙生成的图形如图3.8所示，连接所有的线生成面，生成完整的齿轮面如图3.9所示。

图3.8 齿轮线 图3.9 小齿轮面

先创建一个半径为20的圆，把线创建成面，再单击Main>Preprocessor>Modeling>

- 19 -

基于ANSYS齿轮传动有限元分析及优化

Operate>Subtract>Areas出现如图3.10所示的对话框，用图3.9所示的面减去半径为20的面最后生成的面如图3.12所示。

图3.10 布尔运算减 图3.11 整列线的对话框

图3.12 布尔运算后的齿轮

3.1.7 生成大齿轮

在生成大齿轮的渐开线前在笛卡尔坐标系下先将小齿轮面向右即正向平移227.5mm。

大齿轮的生成原理和小齿轮的生成原理一样，先生成渐开线，渐开线的APDL命令如下：

m=2.5 z1=68 pi=acos(-1) alfa1=20/180*pi ha=1 c=0.25 d1=m*z1

alfa_ding1=acos(z1*cos(alfa1)/(2+z1))

- 20 -

在ANSYS中建立齿轮分析模型

sita_ding1=(tan(alfa_ding1)-alfa_ding1)*180/pi *afun,rad

db1=m*z1*cos(alfa1) rb1=db1/2 df1=d1-2*(ha+c)*m rf1=df1/2 da1=d1+2*ha*m ra1=da1/2 zl1=alfa_ding1/50 *dim,alfa_11,array,50,1 *dim,sita1,array,49,1 *dim,r1,array,49,1 csys,1 alfa_11(1,1)=0 /prep7 k,10000,rb1,0 *do,i,1,49,1

alfa_11(i+1,1)=alfa_11(i,1)+zl1

sita1(i,1)=(tan(alfa_11(i+1,1))-alfa_11(i+1,1))*180/pi r1(i,1)=rb1/cos(alfa_11(i+1,1)) /prep7

k,i+10000,r1(i,1),sita1(i,1) bsplin,i+10000,i+9999 *enddo

k,10050,ra1,sita_ding1 bsplin,10050,10049

再将渐开线上的各线通过布尔运算连接起来操作如下Main>Pre processor>Modeling>Operate>Booleans>Add>Lines。大齿轮的齿根过渡曲线的生成方法与小齿轮的齿根过渡曲线的生成方法略有不同，因为大齿轮的齿根圆直径大于基圆直径，可先利用叠分操作再利用倒角操作对其倒角，倒角的半径r=0.38*m，如图3.13所示。其余创建方法与小齿轮同，生成的大齿轮的面如图3.14。

- 21 -

ANSYS静力加载与静力求解

4 ANSYS静力加载与静力求解

4.1创建接触对

利用ansys中的接触向导创建接触对，选择主菜单中的Preprocessor>Modeling>Cr eat>Contact Pair命令，弹出Contact Manager对话框，如图4.1所示。

图4.1 接触管理对话框

创建步骤如下： （1）单击对话框中的“的对话框，如图4.2所示。

（2）在对话框Target Type栏中选择Flexible如图4.2所示。

（3）选择目标面，在Target Surface下选择Areas，单击Pick Target…弹出拾取对话框。通过鼠标在大齿轮上选择目标面，然后单击OK按钮。

（4）单击向导对话框中的Next按钮，进入下一步，拾取接触面如图4.3所示。 ！注意：接触面不可是刚性面。

（5）在该对话框中，在Target Surface下选择Areas，单击Pick Contact…按钮，弹出图形拾取对话框，用鼠标在大齿轮上选择接触面。然后单击OK按钮。

（6）单击对话框中的Next按钮，进入接触设置对话框，如图4.4所示。 （7）单击对话框中的Optional settings…按钮，弹出的设置对话框如图4.5所示。 （8）完成设置，单击对话框中的OK按钮，回到接触向导对话框。

（9）单击对话框中的Create按钮，完成接触向导，创建好一个接触对，并出现如图4.6的信息提示，如果接触对正确单击Finish完成接触对的创建。

按照这个步骤创建好所有可能接触的所有接触对，见图4.7所示。

”创建接触向导（Contact Wizard）按钮，弹出接触向导

图4.2 接触向导 图4.3 拾取接触面对话框

- 27 -

基于ANSYS齿轮传动有限元分析及优化

图4.4 设置接触参数 图4.5 设置接触对属性对话框

图4.6 提示对话框

图4.7 三个接触对

回到接触管理对话框，可以看到所创建的所有接触对，如图4.8所示。

图4.8 所创建的接触对

在完成接触对创建后，还可以通过在接触管理对话框中选择该接触对，并单击属性按钮“

”进入接触属性修改对话框修改其属性，在这里重点介绍下对接触属性的设

- 28 -

ANSYS静力加载与静力求解

置，因为接触属性设置的正确与否直接影响到求解的收敛与正确，对三个接触对的接触属性的设置见图4.9。

图4.9 对接触属性的设置

选择拉格朗日法的原因是拉格朗日法容易得到良性条件，对接触刚度的敏感性较小。也可以通过设置实常数和关键字来实现，打开CONTA174 element type options菜单如图4.10设置这三个接触对的关键字。

图4.10 对关键字的设置

4.2 ANSYS施加边界条件和加载

单击应用菜单的Select>Entities在弹出的对话框中选择Lines/By num and pick选择小齿轮的内径的四条线，之后再单击应用菜单的Select>Entities在弹出的对话框选择Nodes/Attach to/Lines All,在单击Preprocessor>Modeling>Move/Modify>Rotate Node Cs

- 29 -

基于ANSYS齿轮传动有限元分析及优化

（注意：在做这一步的时候要将当前坐标系设置为坐标系12），则小齿轮的内径圆上的四条线全部转换为柱坐标系，此时X,Y分别代表R,θ。单击Preprocessor>Solution>Loads >Define Loads> Apply>Structural>Displacement> On Nodes如图4.9所示，单击Pick All在弹出的对话框中单击UX如图4.10，即约束小齿轮的径向。再次单击Preprocessor>Solution>Loads>Define Loads> Apply> Structural >Force/Moment > On Nodes如图4.11所示，单击Pick All在弹出的对话框中选择Fy，输入Fy的值为-200N如图4.12，则小齿轮的边界条件和载荷施加完毕如图4.13。同理，大齿轮的边界条件是约束其内径四条线上所有节点的所有自由度其操作与小齿轮施加边界条件的操作相似，完成大齿轮的边界条件后。单击应用菜单的Select>Everything约束的结果见图4.14。

图4.9 施加边界条件拾取对话框 图4.10 约束选择

图4.11施加载荷拾取对话框 图4.12 施加载荷

- 30 -

ANSYS静力加载与静力求解

图4.13 施加边界条件和载荷的小齿轮 图4.14施加边界条件和载荷的大小齿轮

4.3 ANSYS求解

在求解前先对求解控制进行设置单击Preprocessor>Solution>Analysis Type>Sol?n Control出现如图4.15所示的对话框，并按图4.15中的设置进行设置。完成设置后单击OK推出设置。

图4.15 求解控制设置

再单击Main Menu> Preprocessor>Solution> Solve>Current LS出现如图4.16的信息提示框，单击Ok开始求解。当出现“Solution is done”，求解完成。

图4.16 信息提示框

- 31 -

基于ANSYS齿轮传动有限元分析及优化

5 求解结果的分析和优化

5.1 求解结果查看

（1）变形图：单击Main Menu>Preprocessor>General Postproc>Plot Results>Def ormed Shape,弹出Plot Deformed Shape对话框，选择Def+undf edge，如图5.1，单击OK，变形图如图5.2

图5.1 Plot Deformed Shape对话框

图5.2 变形图

（2）节点位移图：

单击Main Menu>Preprocessor>General Postproc>Plot Results>Contour Plot> Nodal Solu，弹出Contour Nodal Solution Data对话框，如图5.3所示。选择DOF Solution/Displacement vector sum,单击OK，节点位移云图如图5.4所示。

- 32 -

求解结果的分析和优化

图5.3 Contour Nodal Solution Data对话框

图5.4 节点位移云图

（3）节点应力云图：单击Main Menu>Preprocessor>General Postproc>PlotResu lts>Contour Plot>Nodal Solu出现如图5.3所示的对话框，单击Stress>Von Miss stresss再点击OK,出现如图5.5、5.6所示的应力云图。

- 33 -

基于ANSYS齿轮传动有限元分析及优化

图5.5 节点应力云图

图5.6 啮合部位节点应力云图

5.2 结果分析及结论

从图5.4 节点位移云图中得出：从轮齿的节点位移情况看, 小轮齿距离啮合部位最

- 34 -

求解结果的分析和优化

远的部分节点位移最大，在大小齿轮的啮合部分的节点位移最小，但就大齿轮而言在与小齿轮的啮合部分的节点位移最大，距离啮合部分一段距离后没有节点位移。

从图5.5、5.6节点应力云图中得出：齿根应力分布特性是衡量齿轮传动性能的重要指标。啮合齿对中, 主动小齿轮的啮合轮齿齿根处Von Mises 应力分布，从动大齿轮的啮合轮齿齿根处Von Mises 应力分布如图5.6所示。从轮齿啮合状态下的应力和变形分布情况看, 齿轮轮齿的应力和变形主要分布啮合齿对上, 在与之相邻的轮齿和齿轮本体上, 则随着啮合点距离的增加而迅速减小。在啮合齿对上, 又以啮合点处的应力和变形最大。在啮合齿对的齿根处, 其应力值也较大。

利用ANSYS对齿轮接触问题进行分析可以得到较为精确的节点位移图和应力云图。

在接触问题的数值模型中，接触单元的选择、网格的划分、界面约束方法、摩擦模型、求解方法、接触刚度、初始接触条件、加载的方式都是影响计算结果的重要因素。 5.3 对齿轮的优化

从图5.4、5.6、5.7中看出齿轮最大应力集中在两齿啮合的部位和齿根部位，也就是说如果齿轮失效这些部位将最先失效，如何优化齿轮的关键就在于如何优化这些部位减少这些部位的应力增强其受力能力。

（1） 可采用更高强度的材料，但可能成本上升较大。

（2） 对齿轮的啮合部位和齿根部位进行加工处理增大它们的强度和韧性。 （3） 对啮合部位经常进行润滑处理减小摩擦。

（4） 采用更先进的技术加工更先进的齿形以改善齿轮啮合时的线接触。

- 35 -

基于ANSYS齿轮传动有限元分析及优化

6 全文总结与展望

6.1 全文总结

在完成本论文的过程中遇到了很多的问题，第一个问题是如何作出齿廓线，第二个

问题是如何划分网格，第三个是如何创建接触对，第四个是如何进行求解控制。

（1）在ANSYS中创建模型对一个初学ANSYS软件的人来说是有相当的难度的，我通过查阅资料自己学习和向老师请教最后编写了渐开线的APDL命令流，其余的工作大多利用界面操作完成。由于这是在ANSYS中直接建模，利用参数化建模技术，只要改变相应的齿轮参数就可以生成相应的齿轮模型，为齿轮啮合传动分析的建模，节约了大量的时间，提高了效率。

（2）划分网格的方法的有很多种，如何在这么多的方法中选择一个适用于我创建的齿轮模型费了很大的功夫，通过学习和查阅资料最后确定了利用扫略网格划分方式，在第一次对齿轮面划分网格的时候总划分失败，警告时有小角度。最后通过重新创建模型解决了这个问题。

（3）创建接触对相对比较简单，但是接触属性的设置至关重要，无论是算法、接触刚度、接触情况还是接触点的位置等等只要修改其中的一项得到的结果就可能千差万别，还可能导致不收敛。

（4）对于非线性的求解问题，加载是通过多载荷步方式逐步加载上的，对于每个载荷步，应以足够的载荷子步保证求解的准确。

虽然整个过程比较漫长和艰难，但是通过老师的帮助和自己的努力都一一克服。并通过这次的论文获益良多。 6.2 工作展望

尽管本论文在ANSYS中建模，同时对齿轮的变形和应力进行了分析，但由于在ANYS创建模型是从上到下建模方式，在创建齿廓线时只是相对精确。还有就是ANSYS的求解结果的正确性只能靠用户自己保证。

由于受到时间和本人知识的能力水平的限制，本文所做的工作还不够完整全面，任然有以下问题有待进一步的开展研究：

（1）本文只是对齿轮进行了接触分析，还可以对齿轮进行模态分析和瞬态动力分析，这一部分内容有待于作进一步的研究分析。

（2）本文只对齿轮进行静力分析，而且没考虑惯性、温度等，这一部分内容有待于作进一步的研究分析。

- 36 -

本文来源：https://www.bwwdw.com/article/w613.html