齿轮渐开线齿廓的方程（斜齿轮）
function jiankaixianxiechilun(x)rb=140.95;%基圆半径
alpha=20*pi/180;%分度圆压力角
r=rb/cos(alpha);%分度圆半径
beta=5*pi/180;%分度圆螺旋角
p=r/tan(beta);%螺旋参数
b=50;%齿宽
Z=60;%工具齿轮齿数
mt=2*r/Z;%端面模数
mn=mt*cos(beta);%法向模数
ra=r+mn;%齿顶圆半径
rf=r-1.25*mn;%齿根圆弧半径
delta=4*pi/180;%齿厚控制参数
%%%%%%%%%%%%%%%左侧齿面%%%%%%%%%%%%
for zz=0:1:Z%齿数循环
dt=2*zz*pi/Z
for seta=-b/(2*p):0.05*pi/180:b/(2*p)%螺旋循环
x=[];y=[];z=[];
n=1 ;
for u=sqrt(rf^2-rb^2)/rb:1*pi/180:sqrt(ra^2-rb^2)/rb;%渐开线循环
x2=rb*sin(seta+u+dt)-rb*u*cos(seta+u+dt);
y2=rb*cos(seta+u+dt)+rb*u*sin(seta+u+dt);
z2=p*seta;
x(n)=x2;y(n)=y2;z(n)=z2;
n=n+1;
end;
plot3(x,y,z);
hold on;
end;
end
%%%%%%%%%%%%%%%右侧齿面%%%%%%%%%%
for zz=0:1:Z%齿数循环
dt=2*zz*pi/Z+delta
for seta=-b/(2*p):0.05*pi/180:b/(2*p)%螺旋循环
x=[];y=[];z=[];
n=1 ;
for u=-sqrt(rf^2-rb^2)/rb:-1*pi/180:-sqrt(ra^2-rb^2)/rb;%渐开线循环
x2=rb*sin(seta+u+dt)-rb*u*cos(seta+u+dt);
y2=rb*cos(seta+u+dt)+rb*u*sin(seta+u+dt);
z2=p*seta;
x(n)=x2;y(n)=y2;z(n)=z2;
n=n+1;
end;
plot3(x,y,z);
hold on;
end;
end
%画中心轴线
x=;
y=;
z=[-b/2 b/2];
plot3(x,y,z,'r');
hold on;
%画基圆
x=[];y=[];z=[];
n=1;
for s=0:5*pi/180:2*pi;
x(n)=rb*sin(s);
y(n)=rb*cos(s);
z(n)=0;
n=n+1;
end
plot3(x,y,z,'k');
hold on;
%画分度圆
x=[];y=[];z=[];
n=1;
for s=0:5*pi/180:2*pi;
x(n)=r*sin(s);
y(n)=r*cos(s);
z(n)=0;
n=n+1;
end
plot3(x,y,z,'r');
hold on;
axis equal
页:
[1]