您现在的位置:e-works > 智造书屋 > 书籍列表 > 车辆工程仿真与分析 > 1.2 仿真计算

第一章 实例1 汽车锥齿轮差速器齿轮设计

第二节 1.2 仿真计算

根据表1-1的计算公式,可以编制M文件进行计算,MATLAB可以大大减少人工计算量。具体程序如下:
      z1=input('请输入行星齿轮齿数z1=');
     
      z2=input('请输入半轴齿轮齿数z2= ');

      if mod(2.*z2,4)~=0; %验证半轴齿数是否为4的倍数

            z2=input('半轴齿数应是4的倍数,请重新输入半轴齿轮齿数z2=' )
end

kd=2;   i1=4.48;   i0=4;   eta=0.9;   n=1;   Temax=121.5;

Tce=kd * Temax * i1 * i0 * eta / n;

G=874.8;   m2=1.36;   fai=0.8;   rr=0.37;   im=1;   etam=0.98;

Tcs=G * m2 *f ai * rr / im / etam;

T=[Tce,Tcs];

Td=min(T);%差速器计算转矩取Tce,Tcs中的小值

fprintf(1,'差速器计算转矩Td=%3 .3 fN·m\\ n',Td

disp'——差速器齿轮设计计算结果—— '%  注意,disp后应输入空格

Kb=2.8;

Rb=Kb * Td(1/3);

A0=0.985 * Rb;

F=0.28 * A0;

[m,gama1,gama2]=adm(z1,z2,Td);%调用求模数函数

m=round(m);

d1=m * z1;d2=m * z2;

hg=1.6 * m;%齿工作高

fprintf(1,'差速器齿轮工作高hg=%3.4f mm\ n',hg)

h=1.788 * m+0.051;%全齿高

fprintf(1,'差速器齿轮工作高h=%3 . 4f mm\ n',h)

h2=(0.430+0.370/(z2/z1). 2);%齿顶高

fprintf(1,'半轴齿轮齿顶高h2=%3 . 4f mm\ n',h2

h1=hg-h2; %齿顶高

fprintf(1,'行星齿轮齿顶高h1=%3 . 4f mm\ n',h1

h11=1.788 * m - h1; %齿根高;

fprintf(1,'行星齿轮齿根高h11=%34f mm\ n',h11

h22=1.788 * m - h2;

fprintf(1,'半轴齿轮齿根高h22=%3.4f mm \ n',h22

c=h - hg;%顶隙

fprintf(1,'顶隙c=%3 . 4f mm\ n',c)

delta1=atan(h11/A0);%齿根角

fprintf(1,'行星齿轮齿根角delta1=%3 . 4f度 \ n',delta1)

delta2=atan(h22/A0);

fprintf(1,'半轴齿轮齿根角delta2=%3 . 4f度 \ n',delta2)

gama01=gama1+delta1;%面锥角

fprintf(1,'行星齿轮面锥角gama01=%3 . 4f度 \ n',gama01

gama02=gama2+delta2;

fprintf(1,'半轴齿轮面锥角gama02=%3• 4f度 \ n',gama02

gamar1=gama1-delta1;%根锥角

fprintf(1,'行星齿轮根锥角gamar1=%3• 4f度 \ n',gamar1

gamar2=gama2-delta2;

fprintf(1,'半轴齿轮根锥角gamar2=%3 • 4f度 \ n,gamar2

d01=d1+2 * h1 * cos(gama1);%外圆直径

fprintf(1,'行星齿轮外圆直径d01=%3 • 4f mm \ n',d01

d02=d2+2 * h2 * cos(gama2);

fprintf(1,半轴齿轮外圆直径d02=%3 • 4f mm\\ n,d02

x01=d2/2h1 * sin(gama1);%节锥顶点至齿轮外缘距离

fprintf(1,'行星齿轮节锥顶点至齿轮外缘距离x01=%3 • 4f mm \ n',x01)

x02=d2/2h2 * sin(gama2);

fprintf(1,'半轴齿轮节锥顶点至齿轮外缘距离x02=%3 • 4f mm \ n',x02)

fprintf(1,'差速器齿轮模数m=%3 • 1f\  n',m)

%%%%%%%%%%

另外再编制求模数的函数文件adm•m。函数输入参数为z1、z2、Td。输出参数为m、gam、gama2

function[m,gama1,gama2]=adm(z1,z2,Td)

i=z1/z2;

gama1=atan(i);

gama2=atan(1/i);

Kb=2.8;

Rb=Kb * Td• (1/3);

A0=0.985 * Rb;

m=2 * A0 * sin(gama1). /z1;

运行第一个主程序:输入z1=10,z2=16得到下面的计算结果。

请输入行星齿轮齿数z1=10

请输入半轴齿轮齿数z2=16

差速器计算转矩Td=359*346N·m

——差速器齿轮设计计算结果——

差速器齿轮工作高hg=3.2000mm

差速器齿轮工作高h=3.6270mm

半轴齿轮齿顶高h2=0.5745mm

行星齿轮齿顶高h1=2.6255mm

行星齿轮齿根高h11=0.9505mm

半轴齿轮齿根高h22=3.0015mm

顶隙c=0.4270mm

行星齿轮齿根角delta1=0.0484度

半轴齿轮齿根角delta2=0.1519度

行星齿轮面锥角gama01=0.6070度

半轴齿轮面锥角gama02=1.1641度

行星齿轮根锥角gamar1=0.5102度

半轴齿轮根锥角gamar2=0.8603度

行星齿轮外圆直径d01=24.4528mm

半轴齿轮外圆直径d02=32.6090mm

行星齿轮节锥顶点至齿轮外缘距离x01=14.6085mm

半轴齿轮节锥顶点至齿轮外缘距离x02=15.5128mm

差速器齿轮模数m=2.0