带有小齿齿廓的啮合齿轮matlab代码clear all clc 带有小齿的齿廓程序 。 N=1000; T=02piN2pi; r=50;基圆半径为50 h=5;齿高 z1=6;齿数 z2=z1; d=100;中心距 ri_1=dz1z1+z2;节圆半径 h1=08;小齿齿高 m=6;小齿齿数 z_1=z1m; z_2=z_1; l=h1sinz1mT;极坐标系下的齿廓 u=hsinz1T+ r

clear all clc

% 定义参数 N=1000; T=0:2pi/N:2pi; r=50; % 基圆半径为50 h=5; % 齿高 z1=6; % 齿数 z2=z1; d=100; % 中心距 ri_1=dz1/(z1+z2); % 节圆半径 h1=0.8; % 小齿齿高 m=6; % 小齿齿数 z_1=z1m; z_2=z_1; l=h1sin(z1m*T); % 极坐标系下的齿廓

% 绘制齿廓 u=hsin(z1T)+ r; c= (hsin(z1T)+r)+h1sin(z1m*T); axis equal plot(T,l,T,u); axis equal plot(T,c); x=c.*cos(T); y=c.*sin(T); [t,r]=cart2pol(x,y); t=unwrap(t); axis equal plot(x,y,'b-'); hold on xi_1=ri_1.*cos(T); yi_1=ri_1.*sin(T); axis equal plot(xi_1,yi_1,'m'); hold on

% 计算法向量斜率 dxdt=((z1hcos(z1T))+z1mh1cos(z1mT)).cos(T)+((hsin(z1T)+r)+h1sin(z1mT)).-sin(T); dydt=((z1hcos(z1T))+z1mh1cos(z1mT)).sin(T)+((hsin(z1T)+r)+h1sin(z1m*T)).cos(T); tt=atan2(dydt,dxdt); tt_1=unwrap(tt,1.8pi); tn=tt_1-pi/2;

% 求齿廓和节圆的交点 [xi,yi]=polyxpoly(xi_1,yi_1,x,y);

% 将交点旋转到节点时 xc=xicos(T)-yisin(T); yc=xisin(T)+yicos(T);

% 绘制啮合点 for i=1:length(T) [x0(i),y0(i)]=polyxpoly(c,sin(tn(i))*xc+cos(tn(i))*yc,x,y); end

% 按照传动比关系求共轭点 z_1=6; z_2=8; d_1=dz_2/(z_1+z_2); ri_2=d_1z_2/(z_1+z_2); theta=2piz_2/(z_1+z_2); x1=ri_2cos(T-theta); y1=ri_2sin(T-theta);

% 绘制共轭点 plot(x0,y0,'ro'); plot(x1,y1,'k-'); axis equal title('带有小齿齿廓的啮合齿轮')