资源描述：
《齿轮接触有限元分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第20卷第2期 计算力学学报 Vol.20,No.22003年4月ChineseJournalofComputationalMechanicsApril2003文章编号:100724708(2003)0220189206齿轮接触有限元分析杨生华(煤炭科学研究总院上海分院,上海200030)摘 要:通过接触仿真分析研究了通用接触单元在轮齿变形和接触应力计算中的应用。建立了一对齿轮接触仿真分析的模型,并使用新的接触单元法计算了轮齿变形和接触应力,与赫兹理论比较,同时也计算了摩擦力对接触应力的影响。计算分析了单元离散、几何、边界范围与加载或约束处理方式的误差,建立了一个计算轮齿变形和接
2、触应力的标准,说明了新的接触单元法的精确性、有效性和可靠性。关键词:接触单元;轮齿变形;接触应力;计算标准;仿真分析中图分类号:TP391 文献标识码:A明了新的接触单元法的精确性、有效性和可靠性。1 引 言建立了一个计算轮齿变形和接触应力的标准或基计算接触非线性问题有许多方法,例如罚函数准,给力学研究和机械设计人员一个参考。法、拉格朗日乘子法等,其中罚函数法由于其经济2 通用接触单元的赫兹计算和方便而得到广泛使用。过去使用点2点接触单元,求解接触问题,对于象齿轮类接触,模型构造很麻为了检验通用接触单元的有效性和精确性,赫烦,计算结果精度和准确性很难保证。随着计算机兹计算验证
3、是必要的。两无限长圆柱有限元计算网和有限元法的发展,新的接触单元法产生精确的几格模型如图1所示。结构单元是具有附加形状函数何模型,自动划分网格,自适应求解。新的单元计算的四节点等参单元(一次单元)。图中接触处网格边精度更高,更有效,功能更强大。其中接触单元能非长为二十分之一接触半宽,该模型节点为7444,单常有效地求解接触非线性问题,新的通用接触单元元为7280(其中接触单元为80个点2面单元)。计(包括点2面和面2面单元)特别适合于计算齿轮接算参数和结果如表1所示,理论结果按公式(1)2触问题。在微机上能实现齿轮接触仿真分析,大大(4)计算[1]。计算结果表明:有限元计算结果和
4、理论地促进了齿轮CAE的形成和发展。计算结果一致,两圆柱变形计算误差仅分别为0.轮齿变形的有限元分析20世纪70年代已开08%和0.045%。注意到公式(2)、(4)是按赫兹接始,但仅仅计算挠曲变形。接触变形和接触应力的触半无限空间推导的公式,因而是理论近似的(变有限元分析在20世纪90年代才真正开始。总之,形误差为1.7%、0.6%,应力误差为0.6%、0.过去的计算是基于试验的计算方法,计算方法是简4%),在接触点不远处一点的变形和应力与有限元化的、近似的,不够精确更不够可靠;没有使用有限计算结果基本一致,有限元计算结果略大于公式元法研究轮齿接触变形和应力,并说明与赫兹变形[
5、1](2)和(4)与理论一致。和应力之间的差别,没有分析计算误差,没有考虑齿轮本体变形对轮齿变形的影响,没有计算摩擦力对接触应力的影响。文中使用ANSYS大型通用有限元分析软件,在个人计算机上建立齿轮接触仿真分析模型。通过两圆柱赫兹接触变形和应力验证其有效性和精度,分析计算了一对直齿轮的轮齿变形和接触应力,说收稿日期:2001204228;修改稿收到日期:20022062241基金项目:上海自然科学基金资助项目1作者简介:杨生华(19632),男,硕士生,工程师1©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreser
6、ved.190计算力学学报 第20卷 表1 两无限长圆柱接触分析Tab.1Twocylindercontactanalysis半径距离理论计算有限元计算参数d(m)S2d2R(mm)d(mm)DH(m)DHmax(Nömm)DH(m)DH(m)Smax(Nömm)圆柱138.0737.1986.1513.898229.5786.1463.966230.85圆柱263.6706.0056.7103.701229.5786.7073.722230.28 注:弹性模量:E521=E2=2.06×10Nömm,泊松比:M1=M2=0.3,kn=10E12表2 有摩擦接触应力分析(单
7、位:Nömm)2Tab.2Contactstressanalysiswithfriction(Unit:Nömm)摩擦系数00(0.001)1(10)20.1(0.001)1(42)20.2(0.05)1(93)20.2(0.001)1(14)2计算应力赫兹应力S3333max(误差)误差Smax(误差)误差Smax(误差)误差Smax(误差)误差圆柱1229.578231.07(3.88)1.0%234.58(4.04)2.2%245.19(3.32)6.8%245.55(0.
