“单边不对称摆线-销齿圆弧型线”计算程序（matlab）

渴望

%% “单边不对称摆线-销齿圆弧型线”计算程序
%%  lishu  create 218-12-06 21:00
function c001
A=100;R=0.27*A;Z1=5;Z2=6;e=1;% 要输入的参数
%%
i1=Z1/Z2; k=i1+1; %齿数比
R1t=(Z1/(Z1+Z2))*A; %阳转子节圆
R2t=(Z2/(Z1+Z2))*A; %阴转子节圆
Alpha=0.035; % 齿顶保护角 α   Alpha=0.035
dd=sin(Alpha);
%% 曲线方程CD
b1=(R^2+R1t^2+2*R*R1t*cos(Alpha))^(1/2);%
Beta=asin(R*dd/b1) ; % β

RhoC=(R^2+R2t^2-2*R2t*R*cos(Alpha))^(1/2);%
RhoD=R2t-e;
qCDC=Beta-acos((A^2+b1^2-RhoC^2)/(2*A*b1));
qCDD=Beta-acos((A^2+b1^2-RhoD^2)/(2*A*b1));
qCD=linspace(qCDC,qCDD,200);

xCD=A*cos(i1*qCD)-b1*cos(Beta-k*qCD);
yCD=A*sin(i1*qCD)+b1*sin(Beta-k*qCD);
zCD=0*qCD;
cCD=[xCD',yCD',zCD'];
save cCD.txt cCD -ascii %保存为文本格式
%% 阳转子 曲线IJ
xD=A*cos(i1*qCDD)-b1*cos(Beta-k*qCDD);
yD=A*sin(i1*qCDD)+b1*sin(Beta-k*qCDD);
Theta2=atan(yD/xD);% θ
RhoI=b1;
Theta3=acos((R2t-e)/R2t);
RhoJ=(A^2+(R2t-e)^2-2*A*(R2t-e)*cos(Theta3))^(1/2);
qIJI=1/i1*(Theta2-acos((A^2+(R2t-e)^2-(RhoI)^2)/(2*A*(R2t-e))));
qIJJ=(1/i1)*(Theta2-acos((A^2+(R2t-e)^2-(RhoJ)^2)/(2*A*(R2t-e))));
qIJ=linspace(qIJI,qIJJ,200);
xIJ=A-(A*cos(qIJ)-(R2t-e)*cos(Theta2-k*qIJ));
yIJ=A*sin(qIJ)+(R2t-e)*sin(Theta2-k*qIJ);
zIJ=0*qIJ;
cIJ=[xIJ',yIJ',zIJ'];
save cIJ.txt cIJ -ascii %保存为文本格式
%% 阴转子 直线DE
RhoDE=linspace(R2t-e,R2t,200);
xDE=RhoDE*cos(Theta2);
yDE=RhoDE*sin(Theta2);
zDE=0*RhoDE;
cDE=[xDE',yDE',zDE'];
save cDE.txt cDE -ascii %保存为文本格式
%% 阳转子 曲线JK
qJKJ=(1/i1)*(Theta2-acos(k*(R2t-e)/A));%                                         
qJKK=(1/i1)*(Theta2-acos(k*R2t/A));
qJK=linspace(qJKJ,qJKK,200);
xJK=A-(A*cos(qJK)-cos(Theta2-k*qJK).*RhoDE);
yJK=   A*sin(qJK)+RhoDE.*sin(Theta2-k*qJK);
zJK=0*qJK;
cJK=[xJK',yJK',zJK'];
save cJK.txt cJK -ascii %保存为文本格式
%% ———————长边————————
%% 曲线 AB
RhoAB1=(R2t^2-R^2)^(1/2);
RhoAB2=R2t;
RhoAB=linspace(RhoAB1,RhoAB2,200);
ThetaAB=asin(R/R2t);
xAB=RhoAB*cos(ThetaAB);
yAB=-RhoAB*sin(ThetaAB);
zAB=0*RhoAB;
cAB=[xAB',yAB',zAB'];
save cAB.txt cAB -ascii
%% 曲线 GH
qGH=1/i1*(acos(k*RhoAB/A)-ThetaAB);
xGH=A-(-RhoAB.*cos(k*qGH+ThetaAB)+A*cos(qGH));
yGH=-RhoAB.*sin(k*qGH+ThetaAB)+A*sin(qGH);
zGH=0*qGH;
cGH=[xGH',yGH',zGH'];
save cGH.txt cGH -ascii
%% 曲线 BC
Alpha1=pi/2-ThetaAB;
tBC=linspace(-Alpha,Alpha1,200);
xBC=R2t-R*cos(tBC);
yBC=-R*sin(tBC);
zBC=0*tBC;
cBC=[xBC',yBC',zBC'];
save cBC.txt cBC -ascii
%%
plot(xCD,yCD,'m',xIJ,yIJ,'r',xDE,yDE,'co',xJK,yJK,'go',xAB,yAB,'rx',xGH,yGH,'go',xBC,yBC,'co')
axis square
xlabel('x轴');
ylabel('y轴');
title('单边不对称摆线-销齿圆弧 ')