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

2、模塊。1、運動學分析1、實例題目對圖示單缸四沖程發(fā)動機中常見的曲柄滑塊機構(gòu)進行運動學仿 真。已知連桿長度：D 0.1m , r3 0.4m，連桿的轉(zhuǎn)速：22 ,33 ,設(shè)曲柄r2以勻速旋轉(zhuǎn)，2 50r/s。初始條件：23 0。仿真以2為輸入，計算3和A，仿真時間0.5 s。2、運動分析建立封閉矢量方程:r2+r3=r1(9)將(9)式分解到x與y軸坐標上，得到：r2cos 0 2+r3cos 0 3=r1r2sin 0 2+r3sin 0 3=0 (10)可得：r1=r2cos 0 2+r3cos 0 3.0 3-arcsin(r2r3)(11)對(10)式對時間求導得：-r2 w2sin 0

3、 2+3 w 3sin 0 3=v1r2 w2cos0 2+3w3cos0 3=0 (12)將上式用矩陣形式表示，令：A= r3sin 0 31-r3cos0 30X= w3v1B=-r2 w2sin 0 2r2 w2cos 0 2則(12)可表示為：AX二B(13)從而可解出w3與v1。、 MATLAB 程序編寫源代碼如下：function varargout = z1(varargin)% Z1 MATLAB code for z1.fig%Z1, by itself, creates a new Z1 or raises the existing% singleton*.%H = Z1

4、returns the handle to a new Z1 or the handle to%the existing singleton*.%Z1(CALLBACK,hObject,eventData,handles,.) calls the local%function named CALLBACK in Z1.M with the given input arguments.%Z1(Property,Value,.) creates a new Z1 or raises the%existing singleton*. Starting from the left, property

5、value pairs are%applied to the GUI before z1_OpeningFcn gets called. An%unrecognized property name or invalid value makes property application% stop. All inputs are passed to z1_OpeningFcn via varargin.%*See GUI Options on GUIDEs Tools menu. Choose GUI allows only one% instance to run (singleton).%

6、See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help z1 % Last Modified by GUIDE v2.5 29-Dec-2016 22:57:13 % Begin initialization code - DO NOT EDIT gui_Singleton = 1;gui_State = struct(gui_Name, m, .gui_Singleton, gui_Singleton, . gui_OpeningFcn, z1_OpeningFcn,

7、.gui_OutputFcn,gui_LayoutFcn,gui_Callback,z1_OutputFcn, . , .);if nargin & ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);end if nargoutvarargout1:nargout = gui_mainfcn(gui_State, varargin:);elsegui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes

8、just before z1 is made visible.function z1_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB % handlesstructure with handles and user data (see GUIDATA)% vara

9、rgin command line arguments to z1 (see VARARGIN)% Choose default command line output for z1 handles.output = hObject;% Update handles structure guidata(hObject, handles);axes(handles.axes3) map1=imread(1.bmp);imshow(map1)% UIWAIT makes z1 wait for user response (see UIRESUME)% uiwait(handles.figure1

10、);% - Outputs from this function are returned to the command line. function varargout = z1_OutputFcn(hObject, eventdata, handles)% varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB % handles st

11、ructure with handles and user data (see GUIDATA)% Get default command line output from handles structure varargout1 = handles.output;% - Executes on button press in pushbuttonRun.function pushbuttonRun_Callback(hObject, eventdata, handles) %主要計算程序r2=0.1;%單位 mr3=0.4;%單位 momiga2=str2double(get(handles

12、.edit1,String); % 單位 rad/s x11=1:500 %單位 msfor i=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);end axes(handles.axes1) %制表 1 plot(x11/1000,omiga

13、3);xlabel(時間(t/s)ylabel(連桿角速度3 3 (rad/s) axes(handles.axes2) %制表 2 plot(x11/1000,v3);xlabel( 時間( t/s) ) ylabel( 滑塊速度 v1( m/s ) )% hObject handle to pushbuttonRun (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)

14、function edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,String) returns contents of edit1 as text%str2double(get(hObject,

15、String) returns contents of edit1 as a double% - Executes during object creation, after setting all properties. function edit1_CreateFcn(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not cre

16、ated 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), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);end% - Executes on butt on press in pushbutt on Exit

17、.function pushbutt on Exit_Callback(hObject, even tdata, han dies)ss=questdlg(確認退出？,退出信息窗口！,繼續(xù)仿真！,退出仿真！,退出仿真！);switch sscase退出仿真！delete(ha ndles.figurel);end% hObject% even tdatahan dle to pushbutt on Exit (see GCBO) reserved - to be defi ned in a future version of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)四、使用指南和實例仿真進入MATLAB軟件，打開并運行程序初始界面：對 山的值進行修改，