哈工大机械原理大作业凸轮教材.docx

所选为第4小题。 推程 远休止程 回程 近休止程 Matlab 绘制图像 取=10rad/s 程序代码： x=0:0.01:2*pi; for m=1:length(x); if x(m)=pi/4 s(m)=120*(2*x(m)/pi).^2; v(m)=9600/(pi*pi).*x(m); a(m)=96000/(pi*pi); elseif x(m)=pi/2 s(m)=60-480/(pi*pi).*(pi/2-x(m)).^2; v(m)=9600/(pi*pi).*(pi/2-x(m)); a(m)=-96000/(pi*pi); elseif x(m)=pi s(m)=60; v(m)=0; a(m)=0; elseif x(m)=13/9*pi s(m)=60.*(1-9.*(x(m)-pi)/(4*pi)+sin(9/2.*(x(m)-pi))/(2*pi)); v(m)=-1350/pi*(1-cos(9/2.*(x(m)-pi))); a(m)=-12000*pi/(4/9*pi)^2*sin(9/2*(x(m)-pi)); else s(m)=0; v(m)=0; a(m)=0; end end d=v/10; hold on figure(1) plot(x,s) xlabel(\phi),ylabel(s),title(s-\phi) figure(2) plot(x,v) xlabel(\phi),ylabel(v),title(v-\phi) figure(3) plot(x,a) xlabel(\phi),ylabel(a),title(a-\phi) figure(4) plot(d,s) xlabel(ds/d\phi),ylabel(s),title(ds/d\phi-s) hold off 绘制的图像： 1、 凸轮的图像 确定凸轮基圆和半径 x=0:0.01:2*pi; for m=1:length(x); if x(m)=pi/4 s(m)=120*(2*x(m)/pi).^2; v(m)=9600/(pi*pi).*x(m); a(m)=96000/(pi*pi); elseif x(m)=pi/2 s(m)=60-480/(pi*pi).*(pi/2-x(m)).^2; v(m)=9600/(pi*pi).*(pi/2-x(m)); a(m)=-96000/(pi*pi); elseif x(m)=pi s(m)=60; v(m)=0; a(m)=0; elseif x(m)=13/9*pi s(m)=60.*(1-9.*(x(m)-pi)/(4*pi)+sin(9/2.*(x(m)-pi))/(2*pi)); v(m)=-1350/pi*(1-cos(9/2.*(x(m)-pi))); a(m)=-12000*pi/(4/9*pi)^2*sin(9/2*(x(m)-pi)); else s(m)=0; v(m)=0; a(m)=0; end if a(m)~=0 K=v(m)*10./a(m); k1=tan(30*pi/180); k2=-tan(70*pi/180); if abs(K-k1)0.005 m1=m; m1 end if abs(K-k2)0.5 m2=m; m2 end end end d=v/10; y1=k1*(d-v(m1)/10)+s(m1); y2=k2*(d-v(m2)/10)+s(m2); hold on plot (d,s) plot (d,y1) plot (d,y2) hold off 由图可以选 所以 凸轮理论轮廓绘制 取 x=0:0.01:2*pi; for m=1:length(x); if x(m)=pi/4 s(m)=120*(2*x(m)/pi).^2; p(m)=(100+s(m)).*cos(x(m))+50.*sin(x(m)); q(m)=(100+s(m)).*sin(x(m))-50.*cos(x(m)); pp(m)=960/(pi*pi).*x(m).*cos(x(m))-(1