【Matlab】主减速器齿轮+行星齿轮/半轴齿轮设计计算校核代码

毕设代码

1.主减速器计算校核程序

clc
clear
close all
%%
z1=8;x1=z1;
z2=33;x2=z2;
x3=z1/z2;%传动比的倒数
F=39.33759;x4=F;
E=35;x5=E;
d2=251.9227;x6=d2;
rd=95.25;x7=rd;
x8=47.6589*pi/180;%弧度制
x9=tan(x8);
x10=1.2*x3;
x11=sin(acot(x10));
Rm2=0.5*(d2-F*x11);x12=Rm2;
x13=E*x11/Rm2;
x14=cos(asin(x13));
x15=x14+x9*x13;
x16=x3*Rm2;
Rm1=x15*x16;x17=Rm1;
TR=0.02*z1+1.06;x18=TR;%或TR=1.30;
x19=Rm2/x10+Rm1;
x20=0.0991;%%%%%%%%%%%%%%%%%%%%%不满足书133页要求需要重新试算x20-x65项,原公式为x20=E/x19，试凑法取x20=0.0991;
x21=(1+x20*x20).^0.5;
x22=x20/x21;
x23=asin(x22);
x24=(E-Rm1*x22)/Rm2;
x25=tan(asin(x24));
x26=x22/x25;
x27=cos(atan(x26));
x28=x24/x27;
x29=cos(asin(x28));
x30=(x15-x29)/x28;
x31=x28*(x9-x30);
x32=x3*x31;
x33=x24-x22*x32;
x34=tan(asin(x33));
x35=x22/x34;
x36=atan(x35);
x37=cos(x36);
x38=x33/x37;
x39=asin(x38);
x40=cos(x39);
x41=(x15+x31-x40)/x38;
x42=atan(x41);
x43=cos(x42);
x44=x42-x39;
x45=cos(x44);
x46=tan(x44);
x47=x22/x33;
x48=acot(x47);
x49=sin(x48);
x50=cos(x48);
x51=(Rm1+Rm2*x32)/x37;
x52=Rm2/x50;
x53=x51+x52;
x54=Rm2*x45/x49;
x55=x43*x51/x35;
x56=(x41*x55-x46*x54)/x53;
x57=-atan(-x56);
x58=cos(-x57);
x59=x41*x56/x51;
x60=x46*x56/x52;
x61=x54*x55;
x62=(x54-x55)/x61;
x63=x59+x60+x62;
x64=(x41-x46)/x63;
x65=x64/x58;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x66=x7/x65;
x67=x3*x50;%或x67=1-x3;
x68=x5/x34-x17*x35;%x68=x35*x37;
x69=x37+x40*x3*x50;
zm=x49*x51;x70=zm;
z=x12*x47-x70;x71=z;
Am=x12/x49;x72=Am;
A0=0.5*x6/x49;x73=A0;
x74=A0-Am;
hgm=3.7*x12*x45/x2;x75=hgm;%k=3.7;
x76=x12*x46/x7;
x77=x49/x45-x76;
x78=45/180*pi;%%%%%%%%%%%%%压力角，弧度制
x79=sin(x78);
x80=x78/2;
x81=cos(x80);
x82=tan(x80);
x83=x77/x82;
x84=10560*x83/x2;
x85=0.15;
x86=1.15-x85;
x87=hgm*x85;
x88=hgm*x86+0.05;
x89=(x84*x85)*pi/60/180;%按照双重收缩齿计算，单位为分，需转化为弧度制
x90=sin(x89);
x91=(x84-x89)*pi/60/180;%按照双重收缩齿计算，单位为分，需转化为弧度制
x92=sin(x91);
x93=x87+x74*x90;
x94=x88+x74*x92;
C=0.15*hgm+0.05;x95=C;
h=x93+x94;x96=h;
hg=h-C;x97=hg;
x98=x48+x89;
x99=sin(x98);
x100=cos(x98);
x101=x48-x91;
x102=sin(x101);
x103=cos(x101);
x104=cot(x101);
d02=x93*x50*2+x6;x105=d02;
x106=x70+x74*x50;
x02=x106-x93*x49;x107=x02;
x108=(x72*x90-x87)/x99;
x109=(x72*x92-x88)/x102;
z0=x71-x108;x110=z0;
zR=x71+x109;x111=zR;
x112=x12+x70*x104;
x113=x5/x112;
x114=cos(asin(x113));
x115=x113/x114;
x116=x103*x114;
x117=asin(x116);
x118=cos(x117);
x119=tan(x117);
x120=(x102*x111+C)/x103;
G0=(E*x113-x120)/x114;x121=G0;
x122=x38*x3*x50/x69;
x123=atan(x122);%x123=cos(28/60+9/3600);
x124=x39-x123;%x124=cos(x39-28/60+9/3600);
x125=x117-x36;%x125=cos(x117-x36);
x126=E/x112*(1-x3)-x35*x37;%x126=-E/x112*(1-x3)-x35*x37;
x127=cos(x123)/cos(x124);
x128=E/x34-Rm1*x35+x87*x35*x37;
x129=x118/cos(x125);
x130=x74*x127;
B0=x128+x130*x129+hgm*(E/x112*(1-x3)-x35*x37);x131=B0;%%%%%%%%%%%%%%%%%%%%%%%%%%
x132=x4*x127-x130;
Bi=x128-x132*x129+hgm*(-E/x112*(1-x3)-x35*x37);x133=Bi;
x134=G0+B0;
d01=x119*x134*2;x135=d01;
x136=x70*x100/x99+x12;
x137=E/x136;
x138=asin(x137);
x139=cos(x138);
x140=(x99*x110+x95)/x100;
GR=(E*x137-x140)/x139;x141=GR;
x142=x100*x139;
x143=asin(x142);
x144=cos(x143);
x145=tan(x143);
Bmin=0.1828887;x146=Bmin;
Bmax=0.2490302;x147=Bmax;
x148=x90+x92;
x149=x96-F*x148;
Ai=A0-F;x150=Ai;
%% %%%%%%%%%%%%%书129页式3-25
K=0.9;
RDX2=x72*sin(x44)+(2*K*K*x73*x73-x72*x72*(2-sin(x44)*sin(x44))).^0.5;

2.行星齿轮、半轴齿轮齿形计算、齿轮弯曲应力计算代码

%% %%%%%%行星齿轮、半轴齿轮齿形计算、齿轮弯曲应力计算
clc
clear
close all
%%
z1=12;
z2=18;
m=5;
F=15.98;
h=1.788*m+0.051;
hg=1.6*m;
alfa=22.5*pi/180;%压力角，单位是弧度
zjj=90;%轴交角，单位是度
d1=m*z1;d2=m*z2;
gama1=33.69*pi/180;gama2=56.31*pi/180;%节锥角，单位是弧度
A0=53.2655;
t=3.1416*m;
h21=(0.43-0.37/(z2/z1)/(z2/z1))*m;h11=hg-h21;%h11、h21是齿顶高
h111=1.788*m-h11;h211=1.788*m-h21;%h111、h211是齿根高
c=h-hg;
delta1=atan(h111/A0);delta2=atan(h211/A0);%齿根角,单位弧度
gama01=gama1+delta2;gama02=gama2+delta1;
gamar1=gama1-delta1;gamar2=gama2-delta2;
d01=d1+2*h11*cos(gama1);d02=d2+2*h21*cos(gama2);
X01=d2/2-h11*sin(gama1);X02=d1/2-h21*sin(gama2);
S2=t/2-(h11-h21)*tan(alfa)+0.05*m;S1=t-S2;
B=0.1684;
SX1=S1-S1.^3/6/d1.^2-B/2;SX2=S2-S2.^3/6/d2.^2-B/2;
hx1=h11+S1.^2*cos(gama1)/4/d1;hx2=h21+S2.^2*cos(gama2)/4/d2;
%% %%%%%%%%%%%%%%%差速器齿轮弯曲应力验算
Tje=5826.636;
Tjm=1230.568;
n=4;
Te=Tje*0.6/n;
Tm=Tjm*0.6/n;
K0=1;
Ks=(m/25.4).^(1/4);
Km=1;
Kv=1;
J=0.2672;
sigmaw_e=2000*Te*K0*Ks*Km/Kv/F/z2/m/m/J;%Tje对应弯曲应力，MPa
sigmaw_m=2000*Tm*K0*Ks*Km/Kv/F/z2/m/m/J;%Tjm对应弯曲应力，MPa