这个是我在网上的代码,画的是这个东西,希望对你有用,不过这个是斜齿轮的,你可以进行一点改动!
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;