机械原理作业-程序设计

1、.一、凸轮程序设计（matlab）a.位移程序% t表示转角，s表示位移t=0:0.01:pi/4; %升程阶段加速s=80*(t/pi*2).2;hold onplot(t/pi*180,s);t=pi/4:0.01:pi/2; %升程阶段减速s=40-80*4/pi.2*(pi/2-t).2;hold onplot(t/pi*180,s); t= pi/2:0.01:19*pi/18; %远休止阶段s=40;hold onplot(t/pi*180,s);t=19*pi/18:0.01:4*pi/3; %回程阶段s=40*(1-35*(18*t-19*pi)/(5*pi).4+84*(18

2、*t-19*pi)/(5*pi).5-70*(18*t-19*pi)/(5*pi).6+20*(18*t-19*pi)/(5*pi).7);hold onplot(t/pi*180,s);t=4*pi/3:0.01:2*pi; %近休止阶段s=0;hold onplot(t/pi*180,s);grid onhold offb速度程序（=1）% t表示转角，令1=1t=0:0.01:pi/4; %升程加速阶段v=160/(pi/2).2*t;hold onplot(t/pi*180,v/12.03211);t=pi/4:0.01:pi/2; %升程减速阶段v=160/(pi/2).2*(pi/

3、2-t);hold onplot(t/pi*180,v/12.03211);t= pi/2:0.01:19*pi/18; %远休止阶段v=0;hold onplot(t/pi*180,v/12.03211);t=19*pi/18:0.01:4*pi/3; %回程阶段v=-40*18/5/pi*(140*(18*t-19*pi)/(5*pi).3-420*(18*t-19*pi)/(5*pi).4+420*(18*t-19*pi)/(5*pi).5-140*(18*t-19*pi)/(5*pi).6);hold onplot(t/pi*180,v/12.03211);t=4*pi/3:0.01:

4、2*pi; %近休止阶段v=0;hold onplot(t/pi*180,v/2/12.03211);grid onhold offc加速度程序% t表示转角，令1=1t=0:0.001:pi/4; %升程加速阶段a=160/(pi/2).2;hold onplot(t/pi*180,a/(46.53378);t=pi/4:0.001:pi/2; %升程减速阶段a=-160/(pi/2).2;hold onplot(t/pi*180,a/(46.53378);t= pi/2:0.001:19*pi/18; %远休止阶段a=0;hold onplot(t/pi*180,a/(46.53378);

5、t=19*pi/18:0.001:4*pi/3; %回程阶段a=-40/pi/(5*pi/18).2*(420*(18*t-19*pi)/(5*pi).2-1680*(18*t-19*pi)/(5*pi).3+2100*(18*t-19*pi)/(5*pi).4-840*(18*t-19*pi)/(5*pi).5);hold onplot(t/pi*180,a/(46.53378);t=4*pi/3:0.001:2*pi; %近休止阶段a=0;hold onplot(t/pi*180,a/(46.53378);grid onhold offdds/d-s程序% t表示转角，x(横坐标)表示类速

6、度ds/d，y（纵坐标）表示位移st=0:0.001:pi/4; %升程加速阶段x=640/pi.2*t;y=2*40*(2*t/pi).2;hold onplot(x/10.02676158,y/4,'-r');t=pi/4:0.001:pi/2; %升程减速阶段x=640/pi.2*(pi/2-t);y=40-320*(pi/2-t).2/pi/pi;hold on plot(x/10.02676158,y/4,'-r');t=19*pi/18:0.001:4*pi/3; %回程阶段x=80/pi*(-140*(18*t-19*pi)/(5*pi).3+42

7、0*(18*t-19*pi)/(5*pi).4-420*(18*t-19*pi)/(5*pi).5+140*(18*t-19*pi)/(5*pi).6);y=40*(1-35*(18*t-19*pi)/(5*pi).4+84*(18*t-19*pi)/(5*pi).5-70*(18*t-19*pi)/(5*pi).6+20*(18*t-19*pi)/(5*pi).7);hold on plot(x/10.02676158,y/4,'-r');grid onhold offe求基圆半径和偏距程序x=-40:0.001:60; % 直线Dtdty=tand(60)*x-68.257

8、;hold on plot(x/10.02676158,y/4,'-r');x=-40:0.001:60; % 直线Dtdty=tand(150)*x-6.213;hold onplot(x/10.02676158,y/4,'-r');%t表示转角，x横坐标表示ds/dt,y纵坐标表示位移st=0:0.001:pi/4; %升程加速阶段x=640/pi.2*t;y=80*(2*t/pi).2;hold onplot(x/10.02676158,y/4,'-r');t=pi/4:0.001:pi/2; %升程减速阶段x=640/pi.2*(pi/2

9、-t);y=40-320*(pi/2-t).2/pi.2;hold onplot(x/10.02676158,y/4,'-r');t=19*pi/18:0.001:4*pi/3; %回程阶段x=144/pi.2*(-140*(18*t-19*pi)/(5*pi).3+420*(18*t-19*pi)/(5*pi).4-420*(18*t-19*pi)/(5*pi).5+140*(18*t-19*pi)/(5*pi).6);y=40*(1-35*(18*t-19*pi)/(5*pi).4+84*(18*t-19*pi)/(5*pi).5-70*(18*t-19*pi)/(5*pi

10、).6+20*(18*t-19*pi)/(5*pi).7);hold onplot(x/10.02676158,y/4,'-r');grid onhold offf确定滚子半径程序1.曲率半径%凸轮理论轮廓半径，t为转角，p为曲率半径,dxi为dx/d, dyi为dy/d,i=1,2,3,4,5,6h=40; %升程t0=pi/2; %升程角t01=5.*pi/18; %回程角ts=5.*pi/9; %远休止角ts1=2.*pi/3; %近休止角e=26.866; %偏距s0=(34.552-28.662).0.5;t=linspace(0,pi/4,1000); %升程加速阶

11、段s=80.*(2.*t/pi).2;dx1=cos(t).*640.*t/pi/pi-sin(t).*(40.*(2.*t/pi).2+43.7748)-cos(t).*26.866;dy1=sin(t).*640/pi.2.*t+cos(t).*( 40.*(2.*t/pi).2+ 43.7748)-26.866.*sin(t);p=sqrt(dx1.2+dy1.2); hold onplot(t,p);c=min(p);t=linspace(pi/4,pi/2,1000); %升程减速阶段s=40-320.*(pi/2-t).2/pi/pi;dx2=cos(t).*640/pi.2.*(

12、pi/2-t)-sin(t).*(40-320.*(pi/2-t).2/pi.2+43.7748)-26.866.* cos(t);dy2=sin(t).*640/pi.2.*(pi/2-t)+cos(t).*(40-320.*(pi/2-t).2/pi.2+43.7748)-28.866.* sin(t);hold onplot(t,p);t=linspace(pi/2,19.*pi/18,1000); % 远休止阶段 s=40;dx3 =-sin(t).*(40+43.7748)-28.866.*cos(t);dy3 =cos(t).*(40+43.7748)-28.866.*sin(t)

13、;p=sqrt(dx3.2+dy3.2); hold onplot(t,p);t=linspace(19.*pi/18,24.*pi/18,1000); % 回程阶段s=40.*(1-35.*(18.*t-19.*pi)/(5.*pi).4+84.*(18.*t-19.*pi)/(5.*pi).5-70.*(18.*t-19.*pi)/(5.*pi).6+20.*(18.*t-19.*pi)/(5.*pi).7);dx4=-cos(t).*144/pi.2.*(140.*(18.*t-19.*pi)/(5.*pi).3-420.*(18.*t-19.*pi)/(5.*pi).4+420.*(1

14、8.*t-19.*pi)/(5.*pi).5-140.*(18.*t-19.*pi)/(5.*pi).6)-sin(t).*( 40.*(1-35.*(18.*t-19.*pi)/(5.*pi).4+84.*(18.*t-19.*pi)/(5.*pi).5-70.*(18.*t-19.*pi)/(5.*pi).6+20.*(18.*t-19.*pi)/(5.*pi).7)+43.7748)-28.866.* cos(t);dy4=sin(t).*144/pi.2.*(140.*(18.*t-19.*pi)/(5.*pi).3-420.*(18.*t-19.*pi)/(5.*pi).4+420.

15、*(18.*t-19.*pi)/(5.*pi).5-140.*(18.*t-19.*pi)/(5.*pi).6)+cos(t).*(40.*(1-35.*(18.*t-19.*pi)/(5.*pi).4+84.*(18.*t-19.*pi)/(5.*pi).5-70.*(18.*t-19.*pi)/(5.*pi).6+20.*(18.*t-19.*pi)/(5.*pi).7)+43.7748)-28.866.*sin(t);p=sqrt(dx4.2+dy4.2);hold onplot(t,p);t=linspace(4.*pi/3,pi.*2,1000); %近休止阶段 s=0; dx6 =

16、-28.866.*cos(t)-sin(t).*43.7748;dy6 =43.7748.*cos(t -28.866.*sin(t); p=sqrt(dx6.2+dy6.2);hold onplot(t,p);grid onhold off2.压力角线t=linspace(0,pi/4,1000); % 升程加速阶段x=abs(640/pi/pi*t-28.866)./( 40*(2*t/pi).2+ 43.7748);y=atand(x);hold on plot(t,y);t=linspace(pi/4,pi/2,1000); %升程减速阶段x=abs(640/pi.2*(pi/2-t)

17、-28.866)./( 40-320*(pi/2-t).2/pi.2+43.7748);y=atand(x);hold onplot(t,y);t=linspace(pi/2,19*pi/18,1000); % 远休止阶段 x=28.866./( 40+43.7748);y=atand(x);hold onplot(t,y); t=linspace(19*pi/18,4*pi/3,1000); % 回程阶段x=abs(144/pi.2*(-140*(18*t-19*pi)/(5*pi).3+420*(18*t-19*pi)/(5*pi).4-420*(18*t-19*pi)/(5*pi).5+

18、140*(18*t-19*pi)/(5*pi).6)-28.866)./(40*(1-35*(18*t-19*pi)/(5*pi).4+84*(18*t-19*pi)/(5*pi).5-70*(18*t-19*pi)/(5*pi).6+20*(18*t-19*pi)/(5*pi).7)+43.7748);y=atand(x);hold onplot(t,y);t=linspace(4*pi/3,pi*2,1000); %近休止阶段 x=abs(0*t-28.866)/43.7748;y=atand(x);hold onplot(t,y);hold offgrid ong绘制凸轮理论轮廓线程序%

19、凸轮的理论轮廓，t为转角，x为横坐标，y为纵坐标s0=(34.55.2-26.866.2).0.5;e=26.866;%t表示转角 s表示位移t=0:0.001:pi/4; %升程加速阶段s=80*(2*t/pi).2;x=(s0+s).*cos(t)-e.*sin(t);y=(s0+s).*sin(t)+e.*cos(t);hold onplot(x,y);t=pi/4:0.001:pi/2; %升程减速阶段s=40-320*(pi/2-t).2/pi.2;x=(s0+s).*cos(t)-e.*sin(t);y=(s0+s).*sin(t)+e.*cos(t);hold onplot(x,

20、y);t=pi/2:0.001:19*pi/18; %远休止阶段s=40;x=(s0+s).*cos(t)-e.*sin(t);y=(s0+s).*sin(t)+e.*cos(t);hold onplot(x,y);t=19*pi/18:0.001:4*pi/3;%回程阶段s=40*(1-35*(18*t-19*pi)/(5*pi).4+84*(18*t-19*pi)/(5*pi).5-70*(18*t-19*pi)/(5*pi).6+20*(18*t-19*pi)/(5*pi).7);x=(s0+s).*cos(t)-e.*sin(t);y=(s0+s).*sin(t)+e.*cos(t);

21、hold onplot(x,y);t=4*pi/3:0.001:2*pi;%近休止阶段s=0;x=(s0+s).*cos(t)-e.*sin(t);y=(s0+s).*sin(t)+e.*cos(t);hold onplot(x,y);grid onhold ont=0:0.001:2*pi; %基圆x=34.55*cos(t);y=34.55*sin(t);hold onplot(x,y);t=0:0.001:2*pi; % 偏距圆x=26.866*cos(t);y=26.866*sin(t);hold onplot(x,y);grid onhold offh.绘制实际轮廓线%凸轮理论轮廓半

22、径，t为转角，p为曲率半径,dxi为dx/d, dyi为dy/d,i=1,2,3,4,5,6h=40; %升程t0=pi/2; %升程角t01=5.*pi/18; %回程角ts=5.*pi/9; %远休止角ts1=2.*pi/3; %近休止角e=26.866; %偏距s0=(34.552-28.662).0.5;rr=7.684; %滚子半径t=0:0.01:pi/4; %升程加速阶段s=80*(2*t/pi).2;x1=(s0+s).*cos(t)-e.*sin(t);y1=(s0+s).*sin(t)+e.*cos(t);dx1=cos(t).*640.*t/pi/pi-sin(t).*(

23、40.*(2.*t/pi).2+43.7748)-cos(t).*26.866;dy1=sin(t).*640/pi.2.*t+cos(t).*( 40.*(2.*t/pi).2+ 43.7748)-26.866.*sin(t);X1=x1-rr*dy1./(sqrt(dx1.2+dy1.2);Y1=y1+rr*dx1./(sqrt(dx1.2+dy1.2);hold onplot(x1,y1);plot(X1,Y1);t=pi/4:0.01:pi/2; %升程减速阶段s=40-320*(pi/2-t).2/pi.2;x2=(s0+s).*cos(t)-e.*sin(t);y2=(s0+s).

24、*sin(t)+e.*cos(t);dx2=cos(t).*640/pi.2.*(pi/2-t)-sin(t).*(40-320.*(pi/2-t).2/pi.2+43.7748)-26.866.* cos(t);dy2=sin(t).*640/pi.2.*(pi/2-t)+cos(t).*(40-320.*(pi/2-t).2/pi.2+43.7748)-28.866.* sin(t);X2=x2-rr*dy2./(sqrt(dx2.2+dy2.2);Y2=y2+rr*dx2./(sqrt(dx2.2+dy2.2);hold onplot(x2,y2);plot(X2,Y2);t=pi/2:

25、0.01:19*pi/18; % 远休止阶段s=40;x3=(s0+s).*cos(t)-e.*sin(t);y3=(s0+s).*sin(t)+e.*cos(t);dx3 =-sin(t).*(40+43.7748)-28.866.*cos(t);dy3 =cos(t).*(40+43.7748)-28.866.*sin(t);X3=x3-rr*dy3./(sqrt(dx3.2+dy3.2);Y3=y3+rr*dx3./(sqrt(dx3.2+dy3.2);hold onplot(x3,y3);plot(X3,Y3);t=19*pi/18:0.01:4*pi/3; % 回程阶段s=40.*(

26、1-35.*(18.*t-19.*pi)/(5.*pi).4+84.*(18.*t-19.*pi)/(5.*pi).5-70.*(18.*t-19.*pi)/(5.*pi).6+20.*(18.*t-19.*pi)/(5.*pi).7);x4=(s0+s).*cos(t)-e.*sin(t);y4=(s0+s).*sin(t)+e.*cos(t);dx4=-cos(t).*144/pi.2.*(140.*(18.*t-19.*pi)/(5.*pi).3-420.*(18.*t-19.*pi)/(5.*pi).4+420.*(18.*t-19.*pi)/(5.*pi).5-140.*(18.*t

27、-19.*pi)/(5.*pi).6)-sin(t).*( 40.*(1-35.*(18.*t-19.*pi)/(5.*pi).4+84.*(18.*t-19.*pi)/(5.*pi).5-70.*(18.*t-19.*pi)/(5.*pi).6+20.*(18.*t-19.*pi)/(5.*pi).7)+43.7748)-28.866.* cos(t);dy4=sin(t).*144/pi.2.*(140.*(18.*t-19.*pi)/(5.*pi).3-420.*(18.*t-19.*pi)/(5.*pi).4+420.*(18.*t-19.*pi)/(5.*pi).5-140.*(18

28、.*t-19.*pi)/(5.*pi).6)+cos(t).*(40.*(1-35.*(18.*t-19.*pi)/(5.*pi).4+84.*(18.*t-19.*pi)/(5.*pi).5-70.*(18.*t-19.*pi)/(5.*pi).6+20.*(18.*t-19.*pi)/(5.*pi).7)+43.7748)-28.866.*sin(t);X4=x4-rr*dy4./(sqrt(dx4.2+dy4.2);Y4=y4+rr*dx4./(sqrt(dx4.2+dy4.2); hold onplot(x4,y4);plot(X4,Y4);t=4*pi/3:0.01:2*pi; %近

29、休止阶段 s=0;x6=(s0+s).*cos(t)-e.*sin(t);y6=(s0+s).*sin(t)+e.*cos(t);dx6=-28.866.*cos(t)-sin(t).*43.7748;dy6=43.7748.*cos(t -28.866.*sin(t); X6=x6-rr*dy6./(sqrt(dx6.2+dy6.2); Y6=y6+rr*dx6./(sqrt(dx6.2+dy6.2); hold onplot(x6,y6);plot(X6,Y6);hold offt=0:0.001:2*pi; %基圆x=34.55*cos(t);y=34.55*sin(t);hold on

30、plot(x,y);t=0:0.001:2*pi; % 偏距圆x=26.866*cos(t);y=26.866*sin(t);hold onplot(x,y);grid onhold off二、连杆程序设计（VB） Public s4(3600) As Double Public v4(3600) As Double Public a4(3600) As Double Public s6(3600) As Double Public v6(3600) As Double Public a6(3600) As Double Public s8(3600) As Double Public v8(

31、3600) As Double Public a8(3600) As Double Public XG(3600) As Double Public YG(3600) As Double Public XB(3600) As Double Public YB(3600) As Double Public XC(3600) As Double Public YC(3600) As Double Public pi As Double Public pa As Double Sub main() Dim i As Long Dim f1(3600) As Double Dim RR1 As RR

32、Dim RR2 As RR Dim RRR1 As RRR Dim RRR2 As RRR Dim RRR3 As RRR Dim RPR1 As RPR Dim d1 As Double Dim d2 As Double Dim ss As Double Dim vXG(3600) As Double Dim vYG(3600) As Double Dim aXG(3600) As Double Dim aYG(3600) As Double Dim vG As Double Dim aG As Double Set RR1 = New RR Set RR2 = New RR Set RRR

33、1 = New RRR Set RRR2 = New RRR Set RRR3 = New RRR Set RPR1 = New RPR pi = 3.1415926 pa = pi / 180 For i = 0 To 3600 Step 1 RR1.f = i * pa / 10: RR1.delt = pi RR1.w = 10: RR1.e = 0 RR1.L = 100 RR1.XA = 542.3 RR1.YA = 41.57 RR1.vXA = 0 RR1.vYA = 0 RR1.aXA = 0: RR1.aYA = 0 RR1.cal XB(i) = RR1.XB: YB(i)

34、 = RR1.YB RRR1.XD = RR1.XB: RRR1.YD = RR1.YB RRR1.vXD= RR1.vXB: RRR1.vYD = RR1.vYB RRR1.aXD= RR1.aXB: RRR1.aYD = RR1.aYB RRR1.XB = 303: RRR1.YB = 200 RRR1.vXB = 0 RRR1.vYB = 0 RRR1.aXB = 0 RRR1.aYB = 0 RRR1.Lj = 273: RRR1.Li = 136 RRR1.M = -1 RRR1.cal XC(i) = RRR1.XC YC(i) = RRR1.YC RR2.XA = RRR1.XD

35、: RR2.YA = RRR1.yD RR2.vXA = RRR1.vXD: RR2.vYA = RRR1.vYD RR2.aXA = RRR1.aXD: RR2.aYA = RRR1.aYD RR2.f = RRR1.fj - 0.07157 * pi: RR2.delt = pi RR2.v= RRR1.vj: RR2.e = RRR1.aj RR2.L = 490 RR2.cal XG(i) = RR2.XB: YG(i) = RR2.YB RRR2.XB = RR2.XB: RRR2.YB = RR2.YB RRR2.vXB = RR2.vXB: RRR2.vYB = RR2.vYB RRR2.aXB = RR2.aXB: RRR2.aYB = RR2.aYB RRR2.XD = 249.8: RRR2.YD