斜齿轮变形及啮合刚度仿真剖析.docVIP

斜齿轮变形及啮合刚度仿真剖析 　　［摘 要］齿轮是大型成套工业设备中重要的基础部件之一，正确分析斜齿轮的啮合刚度是对齿轮进行合理的修形设计的基础，因而对提高整套设备的总体性能十分重要。本文在广泛吸收国内外齿轮刚度和修形理论及实验研究成果的基础上，采用材料力学方法、矩阵理论、空间啮合理论、数值分析方法，对大型斜齿轮的受载变形和啮合刚度进行了理论研究，并进行了利用Matlab软件进行仿真计算，编制了求解斜齿轮啮合线上各点的变形与啮合刚度的计算程序，并探讨了齿轮设计参数对轮齿啮合刚度的影响，为工程实际齿轮设计时提供参考。 　　［关键词］斜齿轮 变形 啮合刚度 Matlab 　　［中图分类号］TB ［文献标识码］A ［文章编号］1009－5489(2009)06－0098－02 　　 　　近20年来，采用有限元法、边界元法等现代方法，计算轮齿变形的研究报道频频出现。与此同时，轮齿变形的测量方法也相应得到发展，由机械杠杆式加载测量到激光全息法测量和激光散斑法测量，但这些研究还仅处于对轮齿变形的静态测量阶段。到目前为止，对轮齿变形的动态测量研究尚在起步阶段。到目前为止，已有多种计算轮齿变形和刚度的方法和公式，但由于所采用的分析手段不同，有时计算结果差别较大。 　　轮齿受载弹性变形的计算方法大体上可归为解析法、实验法、数值法。本文基于解析法中材料力学方法，这种方法是把轮齿简化为弹性基础的悬臂梁，然后用材料力学原理及求解方法进行处理，再加上轮齿对之间的局部接触变形和相对弹性变形基础的倾斜变形。 　　一、坐标系的建立 　　斜齿圆柱齿轮的啮合是一个空间啮合问题。为了计算轮齿变形，本文首先建立坐标系如图1所示。 　　图1 啮合齿轮坐标系 　　图中??O??1?Ш酮?O??2?Х直鹞?齿轮1和齿轮2的转动中心，中心距为??a=??O??1O??2???АＭ贾歇?x??1y??1z??1?Ш酮?x??2y??2z??2??是分别与齿轮1和齿轮2固联的旋转坐标系。??x??1O??1z??1?Ш酮?x??2O??2z??2??位于齿轮端面，且使齿轮上各点的??z??1和z??2?У闹捣钦?，??y??1和y??2?Х直鹞?平面??x??1O??1z??1和x??2O??2z??2??上一对啮合轮齿的对称中心线。??x??0O??0y??0??是以??O??1??为原点的固定坐标系。其中，??y??0轴与连心线O??1O??2?б恢拢???x??0轴与??O??1O??2???Т怪保???z??0轴z??1与轴?е睾稀＊?xyz??是齿面啮合接触线上的坐标系，??y?е嵊虢哟ハ咧睾希?且总是指向??z??0?Ъ跎俚姆较颉Ｗ?标系的原点在端面齿廓与接触线的交点。如图2－1中的??I?У恪Ｕ庋?加载就仅在??x?е岬母悍较颍?而其他两个坐标分量为零。 　　由微分几何可得到坐标系??x??1y??1z??1到x??0y??0z??0转换关系??为： 　　??x??0y??0z??01=??cos??θ??1－??sin??θ??100 　　??sin??θ??1??cos??θ??10000100001?Ax??1y??1z??11?? 　　坐标系??x??2y??2z??2?У姜?x??0y??0z??0?ё?换关系为： 　　??x??0y??0z??01=－??cos??θ??2－??sin??θ??200 　　??sin??θ??2－??cos??θ??20Ax??2y??2z??21?? 　　由固定坐标系??(O??1－x??0y??0z??0)向总体坐标系(I－xyz)?У淖?换矩阵为： 　　??xyz1=a????11??a????12??????13??????14??a???┆?21a????22??????23??????24??a????31??a????32??????33??????34??0001?Ax??0y??0z??01?? 　　二、轮齿变形计算 　　1963年，K?Hayashi首先提出了一种求解斜齿轮啮合接触线载荷分布和变形状态的积分方程法。设齿轮宽度为??b??(mm)、齿轮所受的齿面法向载荷为??F??n??(N)、轮齿的总变形量为?Е莫?(μm)，变形指向轮齿中心为正。那么轮齿单位齿宽上的载荷?Е鬲?(N/m)与变形有如下关系： 　　?Е?=∑Ni=1∫K????bi??(x,ε,y)p??i(ε,y)dε+∑Ni=1∫K????ci??(x,ε,y)p??i(ε,y)dε?? 　　式中，??N?А?同时参与啮合的齿对数； 　　??p??i(ε,y)?А?轮齿啮合接触线上的分布载荷； 　　??K????bi??(x,ε,y)?А?轮齿弯曲、剪切变形的综