不动点迭代法在齿形系数和应力修正系数计算中的应用.docVIP

不动点迭代法在齿形系数和应力修正系数计算中的应用 　　摘要: 为了精确计算渐开线齿轮的齿形系数和应力修正系数,在计算过程中引入了不动点迭代法,并编制出计算齿形系数和应力修正系数的Matlab程序。最后,通过实例进行了验证。 　　Abstract: In order to calculate the involute gear’s tooth shape coefficient and stress modification coefficient accurately, fixed-point iteration method is applied to the calculation, and the matlab program which is used to calculate gear-tooth shape coefficient and stress modification coefficient is prepared. 　　关键词: 不动点迭代;齿形系数;应力修正系数;Matlab 　　Key words: fixed-point iteration;gear-tooth shape coefficient;stress modification coefficient;Matlab 　　中图分类号:TH123 文献标识码:A文章编号:1006-4311(2010)13-0147-02 　　 　　0前言 　　计算机科学中的一个基本要素是迭代。正如其名字所表示的含义,迭代是指重复执行一个计算过程,直到找到答案。迭代可以用来求解方程的根、线性和非线性方程组的解以及微分方程的解。函数g(x)的一个不动点是指一个实数P,满足P=g(P)。从图形上分析,函数y=g(x)的不动点是y=g(x)和y=x的交点。本文的主要目的是将不动点迭代法应用于渐开线齿轮齿形系数和应力修正系数的计算中。 　　1齿形系数和应力修正系数的计算[1][2]及迭代法的应用 　　1.1 齿条刀具参数 　　在图1中, 　　A――圆角中心 　　ρfp――圆角半径 　　I――圆角中心A与节线的距离 　　E――A与轮齿对称轴线的距离 　　通过计算[2], 　　ρ== 　　=≈0.38m 　　E=-htanα-ρcosα-(ρ-ρsinα)tanα 　　1.2 危险截面齿厚S的计算 　　确定危险截面齿厚S和弯曲力臂h的力学基础是弹性力学的光弹实验法,即H??ifer30°切线法。文献[1]基于共轭曲面原理中的包络法,文献[2]基于共轭曲面原理中的齿廓法线法,分别求出计算危险截面齿厚S和弯曲力臂h的方法。无论用何种方法都要求解非线性方程。下面仅就涉及到的非线性方程简述如下。 　　如图1所示,取被切轮齿的对称中线为Y轴,设B为刀具圆角中心A到Y轴的距离,B=-E 　　令齿轮逆时针转过φ角,使切削点S正好位于图示的30°切线与过渡曲线的切点处。PS为公法线。设SAP与Y轴夹角为θ, 　　由图2知, 　　AP= 　　又AP= 　　所以= 　　而φ=-θ,当量分度圆半径r= 　　所以θ=+- 　　令G=,H=--, 　　则θ=tanθ-H,即θ=g(θ)。 　　对此非线性方程用不动点迭代法求解。当θ很小时,tanθ≈θ,故可将初始值θ取为 　　θ= 　　由此建立迭代序列{θ},其序列模式为: 　　θ(初始值) 　　θ=g(θ) 　　θ=g(θ) 　　┇ 　　θ=g(θ) 　　θ=g(θ) 　　┇ 　　如果此无限序列趋向一个极限,则此极限就是求解的目标。不动点迭代法序列求解的matlab程序如下: 　　theta0=-H1/(1-2*G/zv1); 　　theta(1)=theta0; 　　 for k=2:max1 　　theta(k)=(2*G/zv1)*tan(theta(k-1))-H1; 　　 err=abs(theta(k)-theta(k-1)); 　　 relerr=err/(abs(theta(k))+eps); 　　thet=theta(k); 　　if (errtol)|(relerrtol) 　　break; 　　end 　　if k==max1 　　disp(maximum number of iterations exceeded) 　　end 　　 end 　　 theta1=thet; 　　k1=k 　　解出θ的同时输出迭代次数k,根据图2,就可求出危险截面齿厚与模数之比,结果为 　　=zsin-θ+- 　　1.3 载荷作用角α的计算 　　由图2(a)可知 　　α=α-γ 　　其中, 　　α=cos 　　r