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

1、系统仿真与matlab综合试题目： 曲柄滑块机构的运动学仿真编号：21难度系数:姓名班级学号联系方式成绩系统仿真与matlab综合试题一、引言3二、运动学分析31、实例题目32、运动分析3三、MATLAB程厅编写5四、使用指南和实例仿真8五、结语102一、引言曲柄滑块机构是指用曲柄和滑块来实现转动和移动相互转换的平面连杆机构，也称曲柄连杆机构。曲柄滑块机构广泛应用于往复活塞式 发动机、压缩机、冲床等的主机构中，把往复移动转换为不整周或整 周的回转运动；压缩机、冲床以曲柄为主动件，把整周转动转换为往 复移动。这里使用运动学知识，对其运动进行解析，并用MAT L A B为其设计仿真模块。二、运动学

2、分析1、实例题目对图示单缸四冲程发动机中常见的曲柄滑块机构进行运动学仿”，连 杆的转速:，,真。己知连杆长度：，4mo.1 ?r 1333222设曲柄r以匀速 旋转，。初始条件：。仿真以为r/50723222So 0. 5，仿真时间输入, 计算和?T31.2、运动分析建立封闭矢量方程：3)(92+r3=rl r轴坐标上,得到：与y将式分解到xr2cos 0 2+r3cos 0 3=rlr2sin 0 2+r3sin 0 3=0 (10)可得：rl=r2cos 0 2+r3cos 0 3(11)0 3=-arcsin(r2/r3)对(10)式对时间求导得：-r2 3 2sin 0 2+r3 &l

3、t;*)3sin 0 3=vlr2 3 2cos ° 2+ r3 3 3cos 0 3=0 (12)将上式用矩阵形式表示，令：A= r3sin 0 31-r3cos 0 303 3 X=vlB=-r2 3 2sin 0 2r2 a 2cos 0 2则(12)可表示为：AX= B。(13)o V133从而可解出与 4三、MATLAB程序编写源代码如下：function varargout = zl(varargin)% Z1 MATLAB code for zl.fig% Zl, by itself creates a new Z1 or raises the existing% si

4、ngleton*. %H = Z1 returns the handle to a new Z1 or the handle to%the existing singleton*.%ZlCCALLBACKhObject.eventDataandles,.) calls the local% function named CALLBACK in Zl.M with the given input arguments.% Zl('Property,/,Value,/.) creates a new Z1 or raises the% existing singleton*. Startin

5、g from the left, property value pairs are% applied to the GUI before zl_OpeningFcn gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to zl_OpeningFcn via varargin.% *See GUI Options on GUIDE'S Tools menu. Choose GUI allows only o

6、ne% instance to run (singleton).% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help zl% Last Modified by GUIDE v2.5 29-Dec-2016 22:57:13% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', m,.'gui-Singleton',

7、 gui-Singleton,.,gui-OpeningFcn1, zl_OpeningFcnz.'gui-OutputFcn1, (©zl-OutputFcn,.,gui-LayoutFcn1,.'gui_Callback1, );if nargin && ischar(vararginl)gui_State.gui_Callback = str2func(vararginl);endif nargout5va ra rgo ut 1: n a rgo ut = gui_mainfcn(gui_State, varargin:);elsegui_ma

8、infcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before zl is made visible.function zl-OpeningFcnfhObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure % eventdata reserved - to be defined in a future

9、 version of MATLAB % handles structure with handles and user data (see GUIDATA)% varargin command line arguments to zl (see VARARGIN)% Choose default command line output for zl handles.output = hObject;% Update handles structureguidata(hObject, handles);axes(handles.axes3)mapl=imread('Lbmp')

10、;imshow(mapl)% UIWAIT makes zl wait for user response (see UIRESUME) % uiwait(handles.figurel);% Outputs from this function are returned to the command line, function varargout = zl-OutputFcnfhObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject hand

11、le to figure% eve nt data reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structure varargoutl = handles.output;% Executes on button press in pushbuttonRun.function pushbuttonRun-Callback

12、fhObject, eventdata, handles)%主要计算程序 r2=0;单位m r3=0.4;%单位m 6 omiga2=str2double(get(handles.editl/String,); %单位 rad/s xll=l:500 %单位 msfori=l: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) l;-r3*cos(theta3(i) 0;X=inv(

13、A)*B;omiga3(i)=X(l,l);v3(i)=X(2,l);endaxes(handles.axesl) %制表 1 plot(xll/1000zomiga3);xlabel(时间(t/s) 1)丫1m0|(，连杆角速度3 3 (rad/s) ,)axes(handles.axes2) %制表 2 plot(xll/1000zv3);xlabel(时间(t/s)")ylabel(,滑块速度 vl (m/s),)% hObject handle to pushbuttonRun (see GCBO)% eventdata reserved - to be defined in

14、 a future version of MATLAB% handles structure with handles and user data (see GUI DATA) function editl-CallbackfhObject, eventdata, handles)% hObject handle to editl (see GCBO)% eve nt data reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUI

15、 DATA) % Hints: getfhObject/String1) returns contents of editl as textstr2double(get(hObject/String,) returns contents of editl as a double% - Executes during object creation, after setting all properties.function editl_CreateFcn(hObject, eventdata, handles)% hObject handle to editl (see GCBO)% even

16、tdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject/BackgroundColor,)/ getfO/defaultUicontro

17、lBackgroundColor1) sehObject/BackgroundColorYwhite');7end% - Executes on button press in pushbuttonExit.function pushbuttonExit_Callback(hObjectz eventdata, handles)ss=questd©确认退出？ ?退出信息窗口！ 7继续仿真！ ?退出仿真！ ?退出仿真！ switch sscase，退出仿真! 1delete(handles.figurel);end% hObject handle to pushbuttonExit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles struc