机械原理matlab编程

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

机械原理（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

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

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

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;