基于MATLAB曲柄滑塊機(jī)構(gòu)運(yùn)動(dòng)仿真報(bào)告

1、*計(jì)算機(jī)仿真技術(shù)matlab報(bào)告*曲柄滑塊機(jī)構(gòu)目錄一、基于GUI的曲柄滑塊機(jī)構(gòu)運(yùn)動(dòng)仿真二、基于simulink的曲柄滑塊機(jī)構(gòu)運(yùn)動(dòng)仿真曲柄滑塊機(jī)構(gòu)1.題目描述題目：對(duì)如圖1所示的曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)過(guò)程進(jìn)行仿真，并用動(dòng)畫(huà)的方式顯示曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)過(guò)程，位移曲線、速度曲線和加速度曲線。圖中，AB長(zhǎng)R2，BC長(zhǎng)R3mm，A點(diǎn)為坐標(biāo)原點(diǎn)。 圖1 曲柄滑塊機(jī)構(gòu)示意圖2.實(shí)現(xiàn)方法利用GUI界面設(shè)計(jì)來(lái)對(duì)曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)過(guò)程進(jìn)行仿真，并用動(dòng)畫(huà)的方式顯示曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)過(guò)程。3.界面設(shè)計(jì)1. Gui 設(shè)計(jì)1）新建GUI：菜單-新建-gui，并保存為test52）界面設(shè)計(jì)：拖拽左側(cè)圖標(biāo)到繪圖區(qū)，創(chuàng)建GUI界面

2、拖拽左側(cè)圖標(biāo)值繪圖區(qū)設(shè)置如下的按鈕最終的仿真界面如圖所示3）代碼添加：進(jìn)入代碼界面4.代碼編程%模型求解a1=str2double(get(handles.edit1,'String');a2=str2double(get(handles.edit2,'String');a3=str2double(get(handles.edit3,'String');a4=str2double(get(handles.edit4,'String');a5=str2double(get(handles.edit5,'String'

3、); a=a1*(1-cos(a4*a5)+0.25*(a1/a2)*(1-cos(2*a4*a5); set(handles.edit6,'String',a); a0=(a4*a1)*(sin(a4*a5)+0.5*(a1/a2)*sin(2*a4*a5); set(handles.edit7,'String',a0); a6=(a4*a4*a1)*(cos(a4*a5)+(a1/a2)*cos(a4*a5); set(handles.edit8,'String',a6);%繪制位移、速度、加速度曲線axes(handles.axes3);r

4、1=str2double(get(handles.edit1,'String');r2=str2double(get(handles.edit2,'String');omiga1=str2double(get(handles.edit4,'String'); x11=1:720;for i=1:720 x1(i)=i*pi/180; %sin(x2(i)=r1/r2*sin(x1(i); x2(i)=asin(-r1/r2*sin(x1(i); x22(i)=x2(i)*180/pi; r3(i)=r1*cos(x1(i)+r2*cos(x2(i

5、); B=-r1*omiga1*sin(x1(i);r1*omiga1*cos(x1(i); A=r2*sin(x2(i) 1;-r2*cos(x2(i) 0; X=inv(A)*B; omiga2(i)=X(1,1); v3(i)=X(2,1);endplot(x11/60,0.5*r1*sin(x1);xlabel('ʱ¼äÖá t/sec')ylabel('Á¬¸ËÖÊÐÄÔÚYÖá&

6、#201;ϵÄλÖÃ/mm')figure(2)plot(x11/60,r3);title('λÒÆÏßͼ')grid onhold off;xlabel('ʱ¼ät/sec')ylabel('»¬¿éλÒÆ r3/mm')figure(3)plot(x11/6

7、0,omiga2);title('Á¬¸Ë½ÇËÙ¶È')grid onhold off;xlabel('ʱ¼ä t/sec')ylabel('Á¬¸Ë½ÇËÙ¶È omiga2/rad/sec')figure(4)plot(x11/60,v3*pi/180);title('»¬&

8、#191;éËÙ¶È')grid onhold off;xlabel('ʱ¼ä t/sec')ylabel('»¬¿éËÙ¶È v3/mm/sec')%繪制表格axes(handles.axes3);grid onaxes(handles.axes1);grid on%制作動(dòng)畫(huà)axes(handles.axes1);hf=figure('name','Ç

9、;ú±ú»¬¿é»ú¹¹');set(hf,'color','r');hold onaxis(-6,6,-4,4);grid onaxis('off');xa0=-5;%»îÈû×󶥵ã×ø±êxa1=-2.5;%»îÈûÓÒ&#

10、182;¥µã×ø±êxb0=-2.5;%Á¬¸Ë×󶥵ã×ø±êxb1=2.2;%Á¬¸ËÓÒ¶¥µã×ø±ê x3=3.5;%תÂÖ×ø±êy3=0;%

11、5;ªÂÖ×ø±êx4=xb1;%ÉèÖÃÁ¬¸ËÍ·µÄ³õʼλÖúá×ø±êy4=0;%ÉèÖÃÁ¬¸ËÍ·µÄ³õ&#

12、202;¼Î»ÖÃ×Ý×ø±êx5=xa1;y5=0;x6=x3;%ÉèÖÃÁ¬Öá³õʼºá×ø±êy6=0;%ÉèÖÃÁ¬Öá³õʼ×Ý×ø&

13、#177;êa=0.7;b=0.7c=0.7 a1=line(xa0;xa1,0;0,'color','b','linestyle','-','linewidth',40); %ÉèÖûîÈûa3=line(x3,y3,'color',0.5 0.6 0.3,'linestyle','.','markersize',300);%ÉèÖ

14、;ÃתÂÖa2=line(xb0;xb1,0;0,'color','black','linewidth',10);%ÉèÖÃÁ¬¸Ëa5=line(x5,y5,'color','black','linestyle','.','markersize',40);%ÉèÖÃÁ¬

15、4;Ë»îÈûÁ¬½ÓÍ·a4=line(x4,y4,'color','black','linestyle','.','markersize',50);%ÉèÖÃÁ¬¸ËÁ¬½ÓÍ·a6=line(xb1;x3,0;0,'color','blac

16、k','linestyle','-','linewidth',10);a7=line(x3,0,'color','black','linestyle','.','markersize',50);%ÉèÖÃÔ˶¯ÖÐÐÄa8=line(-5.1;-0.2,0.7;0.7,'color','y','lin

17、estyle','-','linewidth',5);%ÉèÖÃÆû¸×±Úa9=line(-5.1;-0.2,-0.72;-0.72,'color','y','linestyle','-','linewidth',5);%ÉèÖÃÆû¸×±Úa10=line(-5.1;-5.1

18、,-0.8;0.75,'color','y','linestyle','-','linewidth',5);%ÉèÖÃÆû¸×±Úa11=fill(-5,-5,-5,-5,0.61,0.61,-0.61,-0.61,a,b,c);%ÉèÖÃÆû¸×ÆøÌå len1=4.8;%Á¬

19、;¸Ë³¤len2=2.5;%»îÈû³¤r=1.3;%Ô˶¯°ë¾¶dt=0.015*pi;t=0; while 1 t=t+dt;if t>2*pi t=0;end lena1=sqrt(len1)2-(r*sin(t)2);%Á¬¸ËÔÚÔ˶¯¹ý³ÌÖ

20、кáÖáÉϵÄÓÐЧ³¤¶È rr1=r*cos(t);%°ë¾¶ÔÚÔ˶¯¹ý³ÌÖкáÖáÉϵÄÓÐЧ³¤

21、82;È xaa1=x3-sqrt(len12-(sin(t)*r)2)-(r*cos(t);%»îÈûÔÚÔ˶¯¹ý³ÌÖеÄÓÒ¶¥µã×ø±êλÖà xaa0=xaa1-2.5;%»îÈûÔÚ

22、12;˶¯¹ý³ÌÖеÄ×󶥵ã×ø±êλÖà x55=x3-cos(t)*r;%Á¬¸ËÔÚÔ˶¯¹ý³ÌÖкá×ø±

23、34;λÖà y55=y3-sin(t)*r;%Á¬¸ËÔÚÔ˶¯¹ý³ÌÖÐ×Ý×ø±êλÖà set(a4,'xdata',x55,'ydata',y55);%ÉèÖÃÁ¬¸Ë¶¥µãÔ˶¯ set(a1,'xdata',xaa1-2.5;xaa1,'ydata',0;