机械原理课程设计matlab程序 凸轮机构

clear

h=70;w=2*pi*300/60; %行程h（单位mm），凸轮角转速（rad/s）d1=pi/2; d2=10/180*pi; d3=pi/2; d4=170/180*pi;%算出凸轮的推程角，远休止角，回程角，近休止角（弧度）

d=1:1:360; d0=d/180*pi; %定义向量

%推程

for i=1:45 %等加速运动

s(i)=2*h*d0(i)^2/d1^2;

v(i)=4*h*w*d0(i)/d1^2;

a(i)=4*h*w^2/d1^2;

end

for i=46:90 %等减速运动

s(i)=h-2*h*(d1-d0(i))^2/d1^2;

v(i)=4*h*w*(d1-d0(i))/d1^2;

a(i)=-4*h*w^2/d1^2;

end

for i=91:101

s(i)=s(i-1);

v(i)=0;

a(i)=0;

end

%回程

for i=102:146 %等加速运动

s(i)=h-2*h*d0(i-100)^2/d3^2;

v(i)=-4*h*w*d0(i-100)/d3^2;

a(i)=-4*h*w^2/d3^2;

end

for i=147:192 %等减速运动

s(i)=2*h*(d3-d0(i-100))^2/d3^2;

v(i)=-4*h*w*(d3-d0(i-100))/d3^2;

a(i)=4*h*w^2/d3^2;

end

for i=193:360

s(i)=s(i-1);

v(i)=0;

a(i)=0;

end

r0=61;e0=39;

s0=sqrt(r0^2-e0^2);

for i=1:1:360

x(i)=(s0+s(i))*sin(d0(i))+e0*cos(d0(i));

y(i)=(s0+s(i))*cos(d0(i))-e0*sin(d0(i));

end

subplot(2,2,1);plot(d,s);title('位移线图S(t)');xlabel('时间t s');ylabel('位移s mm')

subplot(2,2,2);plot(d,v);title('速度线图v(t)');xlabel('时间t s');ylabel('速度v mm/s')

subplot(2,2,3);plot(d,a);title('加速度线图a(t)');xlabel('时间t s');ylabel('加速度a mm/s^2')

subplot(2,2,4);plot(x,y);title('凸轮廓线');xlabel('mm');ylabel('mm');hold on

subplot(2,2,4);plot(0,0,'*');