机械原理matlab编程

机械原理（Matlab绘图）
专业班级
学号
姓名：
图一：导杆的角位移曲线和刨刀的位移曲线
x1=linspace(0,2*pi,100);
l1=100;
l3=360;
l4=200;
l6=200；
y=336;
w1=2*pi;
x3=atan((l6+l1*sin(x1))./(l1*cos(x1)));
for i = 1 : 100;
if x3(i)<0
x3(i)=x3(i)+pi;
end
end
S3=l1*cos(x1)./cos(x3);
V23=(-1)*w1*l1*sin(x1-x3);
w3=w1*l1*cos(x1-x3)./S3;
a23=w3.^2.*S3-w1.^2*l1*cos(x1-x3);
A3=(w1.^2.*l1.*sin(x3-x1)-2.*w3.*V23)./S3;
x4= atan((y-l3*sin(x3))./l4);
Se=l3*cos(x3)+l4*cos(x4);
w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4));
Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4);
A4=(x3.^2.*l3.*sin(w3)+x4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,180*x3./pi,180*x1./pi,Se/1000);
xlabel('x1/(°)'),ylabel('x3(°)');
title('导杆的角位移曲线和刨刀的位移曲线');
text(150,108,'x3');
text(200,60,'Se');
grid on;
plotedit on;
图二：导杆的角速度曲线和刨刀的速度曲线x1=linspace(0,2*pi,100);
l1=100;
l3=360;
l4=200;
l6=200;
y=336;
w1=2*pi;
x3=atan((l6+l1*sin(x1))./(l1*cos(x1)));
for i = 1 : 100;
if x3(i)<0
x3(i)=x3(i)+pi;
end
end
S3=l1*cos(x1)./cos(x3);
V23=(-1)*w1*l1*sin(x1-x3);
w3=w1*l1*cos(x1-x3)./S3;
a23=w3.^2.*S3-w1.^2*l1*cos(x1-x3);
A3=(w1.^2.*l1.*sin(x3-x1)-2.*w3.*V23)./S3;
x4= atan((y-l3*sin(x3))./l4);
Se=l3*cos(x3)+l4*cos(x4);
w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4));
Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4);
A4=(x3.^2.*l3.*sin(w3)+x4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,w3,180*x1./pi,Ve/1000);
xlabel('x1/(°)'),ylabel('w3(rad/s)');
title('导杆的角速度曲线和刨刀的速度曲线');
grid on;
text(100,-5,'Ve');
text(100,1.8,'w3');
axis([0 400 -7 4]);
plotedit on;
图三：导杆的角加速度曲线和刨刀的加速度曲线
x1=linspace(0,2*pi,100);
l1=100;
l3=360;
l4=200;
l6=200;
y=336;
w1=2*pi;
x3=atan((l6+l1*sin(x1))./(l1*cos(x1)));
for i = 1 : 100;
if x3(i)<0
x3(i)=x3(i)+pi;
end
end
S3=l1*cos(x1)./cos(x3);
V23=(-1)*w1*l1*sin(x1-x3);
w3=w1*l1*cos(x1-x3)./S3;
a23=w3.^2.*S3-w1.^2*l1*cos(x1-x3);
A3=(w1.^2.*l1.*sin(x3-x1)-2.*w3.*V23)./S3;
x4= atan((y-l3*sin(x3))./l4);
Se=l3*cos(x3)+l4*cos(x4);
w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4));
Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4);
A4=(x3.^2.*l3.*sin(w3)+x4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,A3,180*x1./pi,ae/1000);
xlabel('x1/(°)'),ylabel('A3(rad/s^2)'),grid;
title('导杆的角加速度曲线和刨刀的加速度曲线');
text(250,25,'ae');
text(300,50,'A3');
plotedit on;
图四：凸轮从动件位移曲线h=20;
F1=2*pi/3;
F2=pi/2;
i=1;
for x=0:0.01:2*pi;
if x>=0&x<=2*pi/3;
S(i)=h*(1-cos(pi*x./F1))./2;
elseif x>=2*pi/3&x<=pi;
S(i)=h;
elseif x>=pi&x<=5*pi/4;
S(i)=h-2*h*(x-pi).^2./(F2.^2);
elseif x>=5*pi/4 &x<=3*pi/2;
S(i)= 2*h*(F2-x+pi).^2./(F2.^2); elseif x>=3*pi/2 &x<=2*pi;
S(i)=0;
end
i=i+1;
x=0:0.01:2*pi;
plot(x,S);
title('从动件位移曲线');
axis([0 7 0 25]);
xlabel('x/(rad)');
ylabel('S/(mm)');
grid on;
图五：凸轮的理论轮廓曲线和实际轮廓曲线
e=10;
Rt=10;
Rb=25;
w=2*pi;
F1=2*pi/3;
F2=pi/2;
i=1;
for f=0:0.01:2*pi;
if f>=0&f<=2*pi/3;
S(i)=h*(1-cos(pi*f./F1))./2;
elseif f>=2*pi/3&f<=pi;
S(i)=h;
elseif f>=pi&f<=5*pi/4;
S(i)=h-2*h*(f-pi).^2./(F2.^2);
elseif f>=5*pi/4 &f<=3*pi/2;
S(i)= 2*h*(F2-f+pi).^2./(F2.^2);
elseif f>=3*pi/2 &f<=2*pi;
S(i)=0;
end
i=i+1;
end
f=0:0.01:2*pi;
s0=sqrt((Rb.^2)-(e.^2));
x =e*cos(f)+(s0+S).*sin(f);
y =(s0+S).*cos(f)-e*sin(f);
d1=e*(-1).*sin(f).*w+(s0+S).*cos(f).*w;
d2=(s0+S).*sin(f).*(-1).*w-e*cos(f).*w;
x1=x+Rt*d2./sqrt((d1.^2)+(d2.^2));
y1=y-Rt*d1./sqrt((d1.^2)+(d2.^2));
plot(x,y,'k-',x1,y1,'k-.');
title('凸轮的理论轮廓曲线和实际轮廓曲线'); text(-10,-30,'实际轮廓曲线');
text(-10,-40,'理论轮廓曲线');
axis([-40 50 -55 40 ]);
hold on;
plot(0,0);
text(0,-3,'O');。