凸輪運動Matlab仿真-Matlab課程設(shè)計

1、Matlab課程設(shè)計李俊機自091設(shè)計題目一：凸輪機構(gòu)設(shè)計已知輪廓為圓形的凸輪（圓的半徑為100mm、偏心距為20mm）,推桿與凸輪運動中心的距離20mm,滾子半徑為10mm,請利用matlab仿真出凸輪推桿的運動軌跡和運動特性（速度，加速度），并利用動畫演示出相關(guān)軌跡和運動特性。%總程序代碼clc;clf;clear;p=figure（position,1001001200600）;fori=1:360%畫圓形凸輪R=100;%圓形凸輪半徑A=0:0.006:2*pi;B=i*pi/180;e=20;%偏心距a=e*cos(B);b=e*sin(B);x=R*cos(A)+a;y=R*sin

2、(A)+b;subplot(1,2,1)plot(x,y,b,LineWidth,3);%填充fill(x,y,y)axis(-R-e,R+e,-R-e,R+e+100);set(gca,Xlim,-R-e,R+e)set(gca,Ylim,-R-e,R+e+100)axisequal;axismanual;axisoff;holdon;plot(a,b,og)plot(e,0,or)plot(0,0,or,LineWidth,3)%畫滾子gcx=0;%滾子中心X坐標(biāo)r=10;%滾子半徑gcy=sqrt(R+r)八2-a2)+b;%滾子中心Y坐標(biāo)gx=r*cos(A)+gcx;%滾子X坐標(biāo)gy

3、=r*sin(A)+gcy;%滾子Y坐標(biāo)plot(gx,gy,b,LineWidth,2);%畫其它部分plot(0a,0b,k,LineWidth,4)plot(33,170190,m,LineWidth,4)plot(-3-3,170190,m,LineWidth,4)%畫頂桿gc=120;dgx=00;dgy=gcygcy+gc;plot(dgx,dgy,LineWidth,4);holdoff%畫位移圖sx(i)=B;sy(i)=gcy;subplot(3,2,2)plot(sx,sy,b,LineWidth,3)title(位移線圖)gridonholdoff;%畫速度圖vx(i)=

4、B;vy(i)=20*cos(B)+(40*cos(B).*sin(B)./(121-4*cos(B)2)(1/2);subplot(3,2,4)plot(vx,vy,g,LineWidth,3)title(速度線圖)gridonholdoff;%畫加速度圖ay=(40*cos(B)-2)./(121-4*cos(B)/2)/(1/2)-20*sin(B)-(40*sin(B)2)/(1214*cos(B)/2)/(1/2)-(160*cos(B)/2.*sin(B)2)/(121-4*cos(B)2)(3/2);subplot2,6)plot(ax,ay,r,LineWidth,3),xlabel(B)titleC加速度線圖)gridonholdoff;M=getframe;end截圖附：通過求導(dǎo)求速度和加速度%求速度symsB;a=e*cos(B);b=e*sin(B);s=sqrt(R+r)2-a2)+b;v=diff(s)結(jié)果：v=20*cos(B)+(40*cos(B)*sin(B)/(121-4*cos(B)八2)八(1/2)%求加速度symsB;v=20*cos(B)+(40*cos(B)*sin(B)/(121-4*cos(B)八2)八(1/2);a=diff(v)結(jié)果：a=(40*cos(B廠2)/(121-4*cos(B2(1/2)-20*sin(B)