凸輪運(yùn)動(dòng)Matlab仿真-Matlab課程設(shè)計(jì)

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

2、=r*cos(a)+a;y =r*sin(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)axis equal;axis manual;axis off;hold on;plot(a,b,og)plot(e,0,or)plot(0,0,or,linewidth,3)%畫(huà)滾子gcx=0; %滾子中心 x 坐標(biāo)r=10; %滾子半徑gcy=sqrt(r+r)2-a2)+b; %滾子中

3、心 y 坐標(biāo)gx=r*cos(a)+gcx; %滾子 x 坐標(biāo)gy=r*sin(a)+gcy; %滾子 y 坐標(biāo)plot(gx,gy,b,linewidth,2);%畫(huà)其它部分plot(0 a,0 b,k,linewidth,4)plot(3 3,170 190,m,linewidth,4)plot(-3 -3,170 190,m,linewidth,4)%畫(huà)頂桿gc=120;dgx=0 0;dgy=gcy gcy+gc;plot(dgx,dgy,linewidth,4);hold off%畫(huà)位移圖sx(i)=b;sy(i)=gcy;subplot(3,2,2)plot(sx,sy,b,lin

4、ewidth,3)title(位移線圖)grid onhold off;%畫(huà)速度圖vx(i)=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(速度線圖)grid onhold off;%畫(huà)加速度圖ax(i)=b;ay(i)=(40*cos(b).2)./(121 - 4*cos(b).2).(1/2) - 20*sin(b) - (40*sin(b).2)/(121 -4*cos(b).2).(1/2) -(160*cos(

5、b).2.*sin(b).2)/(121 - 4*cos(b).2).(3/2);subplot(3,2,6)plot(ax,ay,r,linewidth,3),xlabel(b)title(加速度線圖)grid onhold off;m=getframe;end截圖附：通過(guò)求導(dǎo)求速度和加速度%求速度syms b;a =e*cos(b);b =e*sin(b);s=sqrt(r+r).2-a.2)+b;v=diff(s)結(jié)果：v =20*cos(b) + (40*cos(b)*sin(b)/(121 - 4*cos(b)2)(1/2)%求加速度syms b;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(b)2)(1/2) - 20*sin(b) - (4