哈工大機械原理大作業(yè)凸輪

1、所選為第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(

2、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; endendd=v/10;hold onfigure(1)plot(x,s)xlabel('phi'),ylabel('s'),title('s-ph

3、i')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/dphi'),ylabel('s'),title('ds/dphi-s')hold off繪制的圖像：1、凸輪的圖像確定凸輪基圓和半徑x=0:0.01:

4、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

5、(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 endendd=v/1

6、0;y1=k1*(d-v(m1)/10)+s(m1);y2=k2*(d-v(m2)/10)+s(m2);hold onplot (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

7、)-(100+s(m).*sin(x(m)+50.*cos(x(m); qq(m)=960/(pi*pi).*x(m).*sin(x(m)+(100+s(m).*cos(x(m)+50.*sin(x(m);elseif x(m)<=pi/2 s(m)=60-480/(pi*pi).*(pi/2-x(m).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)=480/(pi*pi).*(pi-2*x(m).*cos(x(m)-(100+s(m).*sin(x(m)+50.

8、*cos(x(m); qq(m)=480/(pi*pi).*(pi-2.*x(m).*sin(x(m)+(100+s(m).*cos(x(m)+50.*sin(x(m);elseif x(m)<=pi s(m)=60; p(m)=160.*cos(x(m)+50.*sin(x(m); q(m)=160.*sin(x(m)-50.*cos(x(m); pp(m)=-160.*sin(x(m)+50.*cos(x(m); qq(m)=160.*cos(x(m)+50.*sin(x(m);elseif x(m)<=13/9*pi s(m)=60.*(1-9.*(x(m)-pi)/(4*p

9、i)+sin(9/2.*(x(m)-pi)/(2*pi); 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)=-135/pi.*(1-cos(9/2.*(x(m)-pi).*cos(x(m)-(100+s(m).*sin(x(m)+50*cos(x(m); qq(m)=-135/pi.*(1-cos(9/2.*(x(m)-pi).*sin(x(m)+(100+s(m).*cos(x(m)+50*sin(x(m);else s(m)=0; p(m)=100*cos(x(m)+50*sin(x(m); q(m)=100*sin(x(m)-50*cos(x(m); pp(m)=-100*sin(x(m)+50*cos(x(m); qq(m)=100*cos(x(m)+50*sin(x(m); end X(m)=p(m)-20*qq(m)./sqrt(pp(m).*pp(m)+qq(m).*qq(m); Y(m)=q(m)+2