matlab曲柄滑块机构的运动学仿真

《系统仿真与matlab》综合试题题目：曲柄滑块机构的运动学仿真编号： 21 难度系数： 姓名 班级 学号 联系方式 成绩 TOC\o"1-5"\h\z\o"CurrentDocument"《系统仿真与matlab》综合试题 1\o"CurrentDocument"一、 引言 3\o"CurrentDocument"二、 运动学分析 31、 实例题目 32、 运动分析 3\o"CurrentDocument"三、 MATLAB程序编写 5\o"CurrentDocument"四、 使用指南和实例仿真 8\o"CurrentDocument"五、 结语 10一、引言曲柄滑块机构是指用曲柄和滑块来实现转动和移动相互转换的平面连杆机构，也称曲柄连杆机构。曲柄滑块机构广泛应用于往复活塞式发动机、压缩机、冲床等的主机构中，把往复移动转换为不整周或整周的回转运动；压缩机、冲床以曲柄为主动件，把整周转动转换为往复移动。这里使用运动学知识，对其运动进行解析，并用MATLAB为其设计仿真模块。二、运动学分析1、实例题目对图示单缸四冲程发动机中常见的曲柄滑块机构进行运动学仿真。已知连杆长度：2=0.1m，r=0.4m，连杆的转速：①=0，①=0,设曲柄r2以匀速旋转，o=50r/s。初始条件：0=0=0。仿真以®为222322、运动分析建立封闭矢量方程：

r2+r2+r3=r19)将⑼式分解到X与y轴坐标上，得到：r2cos&2+r3cos03二rlr2sinO2+r3sinO3=0(10)可得：r1=r2cos92+r3cos9393=-arcsin(r2/r3)(ll)对(10)式对时间求导得：「 -r2«2sin92+r3«3sin93=v1r2«2cos92+r3«3cos93=0(12)将上式用矩阵形式表示，令：A=[r3sin93 1-r3cos93 0]X=[s3v1]B=[-r2«2sin92r2«2cos92]则(12)可表示为：AX=B。(13)从而可解出33与vl。三、MATLAB程序编写源代码如下：functionvarargout=zl(varargin)%ZlMATLABcodefor%Zl,byitself,createsanewZlorraisestheexisting%singleton*.%%H=ZlreturnsthehandletoanewZlorthehandleto% theexistingsingleton*.%% Zl('CALLBACK',hObject,eventData,handles,...)callsthelocal%functionnamedCALLBACKinwiththegiveninputarguments.%% Zl('Property','Value',...)createsanewZlorraisesthe%existingsingleton*.Startingfromtheleft,propertyvaluepairsare%appliedtotheGUIbeforezl_OpeningFcngetscalled.An% unrecognizedpropertynameorinvalidvaluemakespropertyapplication% stop.Allinputsarepassedtozl_OpeningFcnviavarargin.%%*SeeGUIOptionsonGUIDE'sToolsmenu.Choose"GUIallowsonlyone%instancetorun(singleton)".%%Seealso:GUIDE,GUIDATA,GUIHANDLES%Edittheabovetexttomodifytheresponsetohelpzl%LastModifiedbyGUIDE29-Dec-20l622:57:l3%Begininitializationcode-DONOTEDITgui_Singleton=l;gui_State=struct('gui_Name',mfilename,...'gui_Singleton',gui_Singleton,...'gui_OpeningFcn',zl_OpeningFcn,...'gui_OutputFcn',zl_OutputFcn,...'gui_LayoutFcn',[],...'gui_Callback',[]);ifnargin&&ischar(varargin{1})=str2func(varargin{1});endifnargout[varargout{1:nargout}]=gui_mainfcn(gui_State,varargin{:});elsegui_mainfcn(gui_State,varargin{:});end%Endinitializationcode-DONOTEDIT%---Executesjustbeforez1ismadevisible.functionz1_OpeningFcn(hObject,eventdata,handles,varargin)%Thisfunctionhasnooutputargs,seeOutputFcn.%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%varargincommandlineargumentstoz1(seeVARARGIN)%Choosedefaultcommandlineoutputforz1=hObject;%Updatehandlesstructureguidata(hObject,handles);axesmap1=imread('');imshow(map1)%UIWAITmakesz1waitforuserresponse(seeUIRESUME)%uiwait;%---Outputsfromthisfunctionarereturnedtothecommandline.functionvarargout=z1_OutputFcn(hObject,eventdata,handles)%varargoutcellarrayforreturningoutputargs(seeVARARGOUT);%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Getdefaultcommandlineoutputfromhandlesstructurevarargout{1}=;%---ExecutesonbuttonpressinpushbuttonRun.functionpushbuttonRun_Callback(hObject,eventdata,handles)%主要计算程序r2=;%单位mr3=;%单位momiga2=str2double(get,'String'));;%单位rad/sxll=l:500%单位msfori=1:500theta2(i)=i*omiga2/1000;theta3(i)=asin(-r2/r3*sin(theta2(i)));B=[-r2*omiga2*sin(theta2(i));r2*omiga2*cos(theta2(i))];A=[r3*sin(theta3(i))1;-r3*cos(theta3(i))0];X=inv(A)*B;omiga3(i)=X(1,1);v3(i)=X(2,1);endaxes%制表1plot(x11/1000,omiga3);xlabel('时间(t/s)')ylabel('连杆角速度®3(rad/s)')axes%制表2plot(x11/1000,v3);xlabel('时间(t/s)')ylabel('滑块速度vl(m/s)')%hObjecthandletopushbuttonRun(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)functioneditl_Callback(hObject,eventdata,handles)%hObjecthandletoeditl(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Hints:get(hObject,'String')returnscontentsofeditlastext% str2double(get(hObject,'String'))returnscontentsofeditlasadouble%---Executesduringobjectcreation,aftersettingallproperties.functioneditl_CreateFcn(hObject,eventdata,handles)%hObjecthandletoeditl(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesempty-handlesnotcreateduntilafterallCreateFcnscalled%Hint:editcontrolsusuallyhaveawhitebackgroundonWindows.%SeeISPCandCOMPUTER.ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');end%---ExecutesonbuttonpressinpushbuttonExit.functionpushbuttonExit_Callback(hObject,eventdata,handles)ss二questdlg('确认退出’，’退出信息窗口!','继续仿真!','退出仿真!','退出仿真!');switchsscase'退出仿真！'delete;end%hObject handletopushbuttonExit(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)四、使用指南和实例仿真进入MATLAB软件，打开并运行程序。初始界面：

对s2的值进行修改，修改为实例中的50*2*pi二，点击开始仿真按钮，得到仿真结果：THrl -亠一： - Lh*1曲柄滑块机构的

运动学仿真|uj2=da2/dt沁列朋皿t ui=dr1.-d1-'JSPHJ-号聖厨曲柄滑块机构的

运动学仿真|uj2=da2/dt沁列朋皿t ui=dr1.-d1-'JSPHJ-号聖厨Jd理r2=DImr3：Wm r 1uj2=314,16rac^i.m®点击退出仿真按钮，进入退出界面：自功优飞