齿轮接触有限元分析.pdfVIP

第 卷第 期 ， 20 2 计 算 力 学 学 报 Vol . 20 No . 2 年 月 Chinese Journal of Computational Mechanics April 2003 2003 4 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号： （ ） 1007-4708 2003 02-0189-06 齿轮接触有限元分析 杨生华 （煤炭科学研究总院上海分院，上海 200030 ） 摘 要：通过接触仿真分析研究了通用接触单元在轮齿变形和接触应力计算中的应用。建立了一对齿轮接触仿 真分析的模型，并使用新的接触单元法计算了轮齿变形和接触应力，与赫兹理论比较，同时也计算了摩擦力对接 触应力的影响。计算分析了单元离散、几何、边界范围与加载或约束处理方式的误差，建立了一个计算轮齿变形 和接触应力的标准，说明了新的接触单元法的精确性、有效性和可靠性。 关键词：接触单元；轮齿变形；接触应力；计算标准；仿真分析 中图分类号： 文献标识码： TP391 A 建立了一个计算轮齿变形和接触应力的标准或基 1 引 言 准，给力学研究和机械设计人员一个参考。 计算接触非线性问题有许多方法，例如罚函数 2 通用接触单元的赫兹计算 法、拉格朗日乘子法等，其中罚函数法由于其经济 和方便而得到广泛使用。过去使用点-点接触单 为了检验通用接触单元的有效性和精确性，赫 元，求解接触问题，对于象齿轮类接触，模型构造很 兹计算验证是必要的。两无限长圆柱有限元计算 麻烦，计算结果精度和准确性很难保证。随着计算 网格模型如图 所示。结构单元是具有附加形状 1 机和有限元法的发展，新的接触单元法产生精确的 函数的四节点等参单元（一次单元）。图中接触处 几何模型，自动划分网格，自适应求解。新的单元 网格边长为二十分之一接触半宽，该模型节点为 计算精度更高，更有效，功能更强大。其中接触单 ，单元为 （其中接触单元为 个点面单 7444 7280 80 - 元能非常有效地求解接触非线性问题，新的通用接 元）。计算参数和结果如表 所示，理论结果按公 1 触单元（包括点面和面 面单元）特别适合于计算 - - ［］ 1 式（）（）计算 。计算结果表明：有限元计算结