齿轮齿根应力和轮齿变形的三维整轮仿真 - 上海超级计算中心.PDF

高性能计算应用 43 齿轮齿根应力和轮齿变形的三维整轮仿真  杨生华 中国煤炭科工集团上海研究院 上海200030  李丽君 上海超级计算中心 上海201203  李根国 上海超级计算中心 上海201203 摘要： 本文使用高性能并行计算机和ANSYS有限元软件，计算和仿真一个真实齿轮整体模型的 三维齿根应力和轮齿变形，显示了二维模型和三维等效模型的计算误差和可靠性。为齿根应力 和轮齿变形提供了一个计算标准，为齿根应力、轮齿变形和齿轮强度研究，齿轮动态、优化 和可靠性设计，齿轮CAE提供新的基础，也为齿轮力学研究提供了方向。另一方面，也测试了 ANSYS软件的有限元并行计算能力，计算是在曙光4000A超级计算机上进行，模型有659万自由 度，5个工况。 关键词：齿根应力，轮齿变形，并行计算，计算标准，CAE 1. 前言 算齿轮齿根应力和轮齿变形。本文对一个采煤机齿 齿轮齿根应力和轮齿变形的计算是齿轮设计 轮进行了计算，计算是在曙光4000A超级计算机上进 和研究的基本问题。过去由于计算机的限制，许多 行，计算模型的自由度为659万，工况5个，计算时 研究和分析人员计算齿根应力和轮齿变形时所用的 模型为二维或三维等效轮齿模型，由于计算机的限 位CPU，内存为8个GB。 制，计算模型存在误差，计算模型的可靠性低，因 此计算结果很难比较，特别是与试验结果不好比 2. 齿轮齿根应力和轮齿变形的计算 较，造成很多重复试验。所以有必要建立一个齿根 表1为计算齿轮和加工刀具的参数，刀具为凸头 应力和轮齿变形计算的标准，给齿轮设计和研究人 滚刀，齿顶高系数为1，齿轮配合的轴直径（齿轮内 员和齿轮CAE研究人员参考。使用高性能计算，现 圆直径）为90毫米，表中凸台高度仅为加工后挖根 在能够对齿轮三维整轮模型进行仿真分析，精确计 凸台高度。 表1 计算齿根应力齿轮参数 模数 齿数 压力角 变位系数 顶圆半径 齿长 刀具圆角 凸台高度 内圆半径 m Z α X ra B ρao Pro rb 8 35 20° 0.865 153.08 100 0.4m 0.2m 45 2.1 齿轮三维计算模型 应力集中，应该选择三维二次等参单元，ANSYS软 载荷均匀或线性分布作用在齿顶，作用载荷大 件中使用单元SOLID95表示，齿轮整轮网格如图2所 小设定为F = 100N/mm，不同工况总载荷大小一样， 示。整轮网格划分分三部分，载荷作用轮齿扇、相 齿轮内圆固定约束。齿轮弹性模量E = 2.06e5 MPa泊 邻轮齿扇（左、右）和其余轮齿扇，每个轮齿扇分 松比ν=0.3。 轮齿和轮缘两部分。载荷作用轮齿扇的轮齿部分由 2.1.1 网格模型 于齿根网格密度高，因此其它部分网格密度高也相 图1是计算齿轮实体几何模型，由于齿轮齿根有 应提高，轮缘部分周向网格密度相对其它轮缘部分 44 《高性能计算发展与应用》 2012年第一期 总第三十八期 最高。由于应力主要集中在载荷作用轮齿扇局部范 围内，因此与它间隔轮齿扇可以使用较稀疏的网 格，相邻轮齿扇的网格是载荷作用轮齿扇和其余间 隔轮齿扇之间的过渡网格。其余轮齿扇网格都一 样，可以重复复制。所有轮缘部分径向网格密度一