齿轮系统的有限元分析

西京学院毕业设计（论文）

题目:齿轮系统的有限元分析

系（院）: 机电工程系

专业: 数控技术

班级: 数控0902

姓名: 方荣稳

学号: 05

指导老师: 李少海

日期: 2011年11月

摘要

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

利用有限元理论和数值分析方法, 对齿轮系统在加载和离心力共同作用下的变形和强度进行了分析, 研究了离心力对该系统的影响和动态响应。利用三维啮合弹塑性接触有限元方法对齿轮进行了接触强度分析, 并基于热弹耦合进行了轮齿的修形计算, 得到轮齿的理想修形曲线, 为齿轮动态设计提供了一种非常有效的方法。

将齿轮系统划分为传统系统和结构系统两部分, 通过轴承把两者耦合起来。采用有限元方法, 建立了实际单级齿轮减速器的有限元动力学模型, 在工作站上用I- DEA S 软件研究了该齿轮系统的固有特性, 所得结果既后映了系统的动力学性能, 又为齿轮系统的动态响应计算和分析奠定了基础。

关键词：齿轮；有限元法; 模态分析；接触; 修形；

目录

第一章绪论 (3)

有限元的概念 (3)

概述 (4)

第二章齿轮系统有限元模型的建立 (6)

第三章 I2DEA S 固有特性的计算方法 (8)

第四章齿轮系统有限元模态分析结果 (10)

结论 (12)

致谢 (14)

参考文献 (15)

第一章绪论

有限元的概念

有限元分析（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。用有限元法不仅能提高计算精度，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

有限元求解问题的基本步骤通常为：

第一步：问题及求解域定义：根据实际问题近似确定求解域的物理性质和几何区域。

第二步：求解域离散化：将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯上称为有限元网络划分。求解域的离散化是有限元法的核心技术之一。

第三步：确定状态变量及控制方法：一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示，为适合有限元求解，通常将微分方程化为等价的泛函形式。

第四步：单元推导：对单元构造一个适合的近似解，即推导有限单元的列式，其中包括选择合理的单元坐标系，建立单元函数，以某种方法给出单元各状态变量的离散关系，从而形成单元矩阵（结构力学中称刚度阵或柔度阵）。为保证问题求解的收敛性，单元形状应以规

则为好，内角避免出现钝角，避免出现畸形，因为畸形时不仅精度低，而且有缺秩的危险，将导致无法求解。

第五步：总装求解：将单元总装形成离散域的总矩阵方程（联合方程组），反映对近似求解域的离散域的要求，即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行，状态变量及其导数（可能的话）连续性建立在结点处。

第六步：联立方程组求解和结果解释：有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量，将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。简言之，有限元分析可分成三个阶段，前处理、处理和后处理。前处理是建立有限元模型，完成单元网格划分；后处理则是采集处理分析结果，使用户能简便提取信息，了解计算结果。

概述

齿轮系统是由齿轮、轴、轴承和箱体等组成的机械结构,在内部和外部激励下将发生机械振动。振动系统的固有特性, 一般包括固有频率和振型, 它是系统的动态特性之一, 对系统的动态响应、动载荷的产生与传递以及系统振动的形式等具有重要的影响。此外, 固有特性还是用振型叠加法求解系统响应的基础。然而, 在齿轮系统的设计阶段, 不能得到系统固有特性的实验数据, 只能通过理论计算得到

进行动力学分析的参数, 目前最好的方法是有限元动力分析方法。由于计算机软、硬件技术的发展, 在设计阶段计算结构的固有特性已成

可能。市面上有许多大型的商业化集成软件可供选择, 如M SCöNA STRAN ,MARC,AN SYS, I- DEA S 等。

齿轮传动的主要特点有高线速度带来较大的振动、噪声和动载荷; 齿轮的变形、应力较大。目前, 国内外学者对齿轮系统的振动进行了广泛的研究。Haruo对弹性轴上斜齿轮对的振动进行分析, 研究了轮齿啮合刚度的变化、啮合面阻尼及轴承阻尼对斜齿轮动态特性的影响。李润方等进行了轮齿啮合过程中应力应变的数值分析和轮齿修形研究。为提高齿轮传动的运行性能, 本文以GSC390F齿轮为对象, 用I- DEA S 软件求解分析了齿轮的离心力对齿轮的影响及其动态响应, 并用自主开发的齿轮三维接触有限元程序计算了该齿轮的接触应力, 用二维热弹接触有限元分析程序系统进行啮合轮齿的温度场及热弹

耦合分析。

第二章齿轮系统有限元模型的建立

齿轮系统分为传动系统(齿轮、传动轴) 和结构系统(主要是箱体) 两部分, 通过轴承把两者耦合起来。下面以渐开线单级圆柱齿轮减速器ZD- 10 为研究对象, 分别建立两个子系统和整个系统的有限元模型。

在建立齿轮轴的动力学模型时, 为减小模型的大小, 将齿轮轮

齿部分简化为分度圆柱, 将齿轮和轴的过盈配合联接看成刚性联接, 忽略键槽的影响。这样的简化对动力学研究来说, 误差很小。通过

I2DEA S 映射网格划分(M appedMesding) 方法, 选用八节点六面体

实体单元, 对齿轮轴部件进行了有限元网格划分。

齿轮箱体的结构比较复杂, 机体上分布有筋板、凸台、轴承孔和各螺栓联接孔等, 在建立有限元模型时, 对箱体实际结构进行了等

效简化处理。用自由网格划分(FreeMeshing)方法, 采用八节点四面

体实体单元, 对齿箱体进行了有元网络划分。

最后, 将高速轴和低速轴装入齿轮箱, 用I2DEA S 的A PPEND 命令将齿轮系统的各零部件集成为整个系统的有限元模型。其中, 每个轴承用四个弹簧模似, 弹簧刚度采用实验模态分析方法得到, 弹簧

刚度值如表所示。弹簧刚度的识别过程参见文献。

表齿轮轴轴承刚度识别结果