matlab曲柄滑塊機構(gòu)的運動學(xué)仿真

1、系統(tǒng)仿真與matlab綜合試題目： 曲柄滑塊機構(gòu)的運動學(xué)仿真編號：21難度系數(shù):姓名班級學(xué)號聯(lián)系方式成績系統(tǒng)仿真與matlab綜合試題一、引言3二、運動學(xué)分析31、實例題目32、運動分析3三、MATLAB程廳編寫5四、使用指南和實例仿真8五、結(jié)語102一、引言曲柄滑塊機構(gòu)是指用曲柄和滑塊來實現(xiàn)轉(zhuǎn)動和移動相互轉(zhuǎn)換的平面連桿機構(gòu)，也稱曲柄連桿機構(gòu)。曲柄滑塊機構(gòu)廣泛應(yīng)用于往復(fù)活塞式 發(fā)動機、壓縮機、沖床等的主機構(gòu)中，把往復(fù)移動轉(zhuǎn)換為不整周或整 周的回轉(zhuǎn)運動；壓縮機、沖床以曲柄為主動件，把整周轉(zhuǎn)動轉(zhuǎn)換為往 復(fù)移動。這里使用運動學(xué)知識，對其運動進行解析，并用MAT L A B為其設(shè)計仿真模塊。二、運動學(xué)

2、分析1、實例題目對圖示單缸四沖程發(fā)動機中常見的曲柄滑塊機構(gòu)進行運動學(xué)仿”，連 桿的轉(zhuǎn)速:，,真。己知連桿長度：，4mo.1 ?r 1333222設(shè)曲柄r以勻速 旋轉(zhuǎn)，。初始條件：。仿真以為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)式對時間求導(dǎo)得：-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繼續(xù)仿真！ ?退出仿真！ ?退出仿真！ 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