石川法，啮合冲击载荷程序

MATLAB程序 ，利用石川法计算齿轮啮合冲击载荷。

clear;
z1=19;z2=45;%两齿轮的齿数
m=1.75;
b=m*z1*0.6;
ha=m;
c=0.25*m;
d1=m*z1;d2=m*z2;
r1=m*z1/2;r2=m*z2/2;%分度圆半径
hf1=1.25*m;
hf2=1.25*m;
alpha=20*pi/180;%分度圆压力角
invalpha=tan（alpha）-alpha;
db1=d1*cos（alpha）;
db2=d2*cos（alpha）;
rb1=db1/2;
rb2=db2/2;%基圆半径
rf1=r1-hf1;
rf2=r2-hf2;%齿根圆半径
da1=d1+2*ha;
da2=d2+2*ha;
ra1=da1/2;
ra2=da2/2;%齿顶圆半径
alpha_a1=acos（rb1/ra1）;
alpha_a2=acos（rb2/ra2）;%齿顶圆压力角
alpha_f1=acos（rb1/rf1）;
alpha_f2=acos（rb2/rf2）;%齿根圆压力角
s=pi*m/2;%分度圆弧齿厚
e=s;%分度圆齿槽宽
sk1=ra1*（s/r1+2*（（tan（alpha）-alpha）-（tan（alpha_a1）-alpha_a1）））;
sk2=ra2*（s/r2+2*（（tan（alpha）-alpha）-（tan（alpha_a2）-alpha_a2）））;%齿顶圆齿厚
PB1=r1*cos（alpha）*（tan（alpha_a1）-tan（alpha））;
PB2=r2*cos（alpha）*（tan（alpha_a2）-tan（alpha））;
B1B2=PB1+PB2;
Pb=pi*m*cos（alpha）;%基圆齿距
Epsilona=B1B2/Pb;
N1B1=sqrt（ra

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       8222  2013-10-28 16:33  shichuan.m

----------- ---------  ---------- -----  ----

                 8222                    1