【刚度计算】石川公式法齿轮时变啮合刚度计算【含Matlab源码 2680期】

一、⛄获取代码方式

获取代码方式1：
完整代码已上传我的资源：【刚度计算】基于matlab石川公式法齿轮时变啮合刚度计算【含Matlab源码 2680期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

clear;clc
z1=45;z2=90;%两齿轮的齿数
m=3;%模数
b=m8;%齿宽
ha=m;%齿顶高
c=0.25m;%顶隙
d1=mz1;d2=mz2;%分度圆直径
r1=mz1/2;r2=mz2/2;%分度圆半径
hf1=1.25m;%齿根高
hf2=1.25m;%齿根高
alpha=20pi/180;%分度圆压力角
invalpha=tan(alpha)-alpha;%压力角的渐开线函数
db1=d1cos(alpha);%基圆直径
db2=d2cos(alpha);%基圆直径
rb1=db1/2;%基圆半径
rb2=db2/2;%基圆半径
rf1=r1-hf1;%齿根圆半径
rf2=r2-hf2;%齿根圆半径
da1=d1+2ha;%齿顶圆直径
da2=d2+2ha;%齿顶圆直径
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=pim/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=r1cos(alpha)(tan(alpha_a1)-tan(alpha));%节点P到脱啮点B1的距离
PB2=r2cos(alpha)(tan(alpha_a2)-tan(alpha));%节点P到脱啮点B2的距离
B1B2=PB1+PB2;%实际啮合线
Pb=pimcos(alpha);%基圆齿距
Epsilona=B1B2/Pb;%重合度
N1B1=sqrt(ra1^2-rb12);
N1B2=N1B1-B1B2;
N2B2=sqrt(ra2^2-rb22);
N2B1=N2B2-B1B2;
rF2=sqrt(rb2^2+N2B12);%有效齿根圆半径
rF1=sqrt(rb1^2+N1B22);%有效齿根圆半径
alpha_F1=acos(rb1/rF1);
alpha_F2=acos(rb2/rF2);
sf1=2rF1sin(pi/2/z1+invalpha-tan(alpha_F1)+alpha_F1);
sf2=2rF2sin(pi/2/z2+invalpha-tan(alpha_F2)+alpha_F2);
h1=sqrt(ra1^2-(sk1/2)2)-sqrt(rf1^2-(sf1/2)2);
h2=sqrt(ra2^2-(sk2/2)2)-sqrt(rf2^2-(sf2/2)2);
hr1=sqrt(rF1^2-(sf1/2)2)-sqrt(rf1^2-(sf1/2)2);
hr2=sqrt(rF2^2-(sf2/2)2)-sqrt(rf2^2-(sf2/2)2);
hi1=(h1sf1-hr1sk1)/(sf1-sk1);
hi2=(h2sf2-hr2sk2)/(sf2-sk2);
N2C=N2B1+Pb;
B2C=B1B2-Pb;%双齿啮合区
CD=Pb-B2C;%单齿啮合区
N1C=N1B1-Pb;
Fn=1000;%外力
E=2e+008;%弹性模量
v=0.26;%泊松比
n=100;
step=B2C/n;%双齿啮合区的距离步长
step3=CD/n;%单齿啮合区的距离步长
nz1=3000;%齿轮1转速r/min
Tz=60/z1/nz1;%循环的周期
t1=B2C/PbTz;%双齿啮合的周期
t2=CD/PbTz;%单齿啮合的周期
step2=t1/n;%双齿啮合区的时间步长
step4=t2/n;%单齿啮合区的时间步长

cycle=5;%循环啮合次数
period=cycle*n;%啮合总时间

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]刘萍.如何理解空气质量分指数（IAQI）计算公式并速算[J].黑龙江环境通报. 2014,38(02)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除