机械原理内燃机课程设计

电算的源程序（MATLAB）

1.滑块的位移源程序

x1=0:0.05:4*pi; %x1--原动件的角度变量

b=150; %b--原动件的长度

c=200; %c--连杆的长度

k=b/c;

l=150*cos(x1)+sqrt(200^2-(150*sin(x1).^2)) %l--连杆的位移plot(x1,l);

title('滑块的位移图像')

xlabel('\it角度','FontSize',8)

ylabel('位移大小','FontSize',8)

2.滑块的速度源程序

x1=0:0.05:4*pi; %x1--原动件的角度变量

b=150; %b--原动件的长度

c=200; %c--连杆的长度

k=b/c;

w=(2*pi*650)/60; %w--原动件的角速度x2=asin(-k*sin(x1)); %x2--连杆的角度

v=b*w*sin(x1-x2)./cos(x2);

plot(x1,v);

title('滑块的速度图像')

xlabel('\it角度','FontSize',8)

ylabel('速度大小','FontSize',8)

3.滑块的加速度源程序：

x1=0:0.05:4*pi; %x1--原动件的角度变量

b=150; %b--原动件的长度

c=200; %c--连杆的长度

k=b/c;

w=(2*pi*650)/60; %w--原动件的角速度x2=asin(-k*sin(x1)); %x2--连杆的角度

a=b*(w^2)*(cos(x1-x2)./cos(x2)+k*cos(x1).^2./cos(x2).^3); plot(x1,a);

title('滑块的加速度图像')

xlabel('\it角度','FontSize',8)

ylabel('加速度大小','FontSize',8)

4.从动件滚子的位移源程序：

e=5; %凸轮的偏心距

r=35; %凸轮的基圆半径

u=0:0.001:360;

s0=sqrt(r^2-e^2);

s1=(4*(1-cos(pi*u/50))).*(0<=u & u<50); %凸轮推程运动阶段

s2=8.*(u>=50&u<55); %远休止阶段

s3=4*(1+cos(pi*(u-55)/50)).*(u>=55&u<105); %凸轮回程阶段

s4=0.*(u>=105&u<=360); %近休止阶段

s=s1+s2+s3+s4;

plot(u,s);

axis([0,360,0,15]);

title('滚子的运动线图');

xlabel('角度');

ylabel('位移');

5.凸轮的轮廓线源程序：

e=5; %凸轮的偏心距

r=35; %凸轮的基圆半径r1=5; %滚子的半径

s0=sqrt(r^2-e^2); %区间分段

x1=50;

x2=55;

x3=105;

x4=360;

w=pi*650/60; %凸轮的角速度u=0:0.05:360;

s1=(4*(1-cos(pi*u/50))).*(0<=u & u

v1=pi*8.*sin(pi*u/50).*(0<=u & u

dx1=((v1-e).*sin(u*pi/180)+(s0+s1).*cos(u*pi/180)).*(0<=u & u

dy1=((v1-e).*cos(u*pi/180)-(s0+s1).*sin(u*pi/180)).*(0<=u & u

X1=((s0+s1).*sin(u*pi/180)+e*cos(u*pi/180)).*(0<=u &

u

Y1=((s0+s1).*cos(u*pi/180)-e*sin(u*pi/180)).*(0<=u &

u

XS1=X1-r1*dy1./sqrt(dx1.^2+dy1.^2);

YS1=Y1+r1*dx1./sqrt(dx1.^2+dy1.^2);

s2=8; %远休止阶段

v2=0;

dx2=((v2-e).*sin(u*pi/180)+(s0+s2).*cos(u*pi/180)).*(x1<= u & u

dy2=((v2-e).*cos(u*pi/180)-(s0+s2).*sin(u*pi/180)).*(x1<= u & u

X2=((s0+s2)*sin(u*pi/180)+e.*cos(u*pi/180)).*(x1<=u &

u

Y2=((s0+s2)*cos(u*pi/180)-e.*sin(u*pi/180)).*(x1<=u &

u

XS2=X2-r1*dy2./sqrt(dx2.^2+dy2.^2);

YS2=Y2+r1*dx2./sqrt(dx2.^2+dy2.^2);

s3=4*(1+cos(pi*(u-x2)/50)).*(x2<=u&u

v3=-pi*8.*sin(pi*(u-x2)/50).*(x2<=u&u

dx3=((v3-e).*sin(u*pi/180)+(s0+s3).*cos(u*pi/180)).*(x2<= u&u

dy3=((v3-e).*cos(u*pi/180)-(s0+s3).*sin(u*pi/180)).*(x2<= u&u

X3=((s0+s3).*sin(u*pi/180)+e.*cos(u*pi/180)).*(x2<=u&u

Y3=((s0+s3).*cos(u*pi/180)-e.*sin(u*pi/180)).*(x2<=u&u

XS3=X3-r1*dy3./sqrt(dx3.^2+dy3.^2);

YS3=Y3+r1*dx3./sqrt(dx3.^2+dy3.^2);

%近休止阶段

v4=0;

s4=0;

dx4=((v4-e).*sin(u*pi/180)+(s0+s4).*cos(u*pi/180)).*(x3<= u&u

dy4=((v4-e).*cos(u*pi/180)-(s0+s4).*sin(u*pi/180)).*(x3<= u&u

X4=(s0*sin(u*pi/180)+e.*cos(u*pi/180)).*(x3<=u & u<=x4); Y4=(s0*cos(u*pi/180)-e.*sin(u*pi/180)).*(x3<=u & u<=x4);

XS4=X4-r1*dy4./sqrt(dx4.^2+dy4.^2);

YS4=Y4+r1*dx4./sqrt(dx4.^2+dy4.^2);