matlab曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)學(xué)仿真設(shè)計(jì)

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

2、。二、運(yùn)動(dòng)學(xué)分析、實(shí)例題目對(duì)圖示單缸四沖程發(fā)動(dòng)機(jī)中常見的曲柄滑塊機(jī)構(gòu)進(jìn)行運(yùn)動(dòng)學(xué)仿真。已知連桿長(zhǎng)度：，連桿的轉(zhuǎn)速：，設(shè)曲柄r2以勻速旋轉(zhuǎn)，。初始條件：。仿真以為輸入，計(jì)算和，仿真時(shí)間0.5s。、運(yùn)動(dòng)分析建立封閉矢量方程：r2+r3=r1（9）將(9)式分解到x與y軸坐標(biāo)上，得到：r2cos2+r3cos3=r1r2sin2+r3sin3=0(10)可得：r1=r2cos2+r3cos33=-arcsin(r2/r3)(11)對(duì)(10)式對(duì)時(shí)間求導(dǎo)得：-r22sin2+r33sin3=v1r22cos2+ r33cos3=0(12)將上式用矩陣形式表示，令：A= r3sin3 1-r3cos3 0

3、X= 3 v1B=-r22sin2r22cos2則(12)可表示為：AX=B。(13)從而可解出3與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 returns the handle to a new Z1 or the handle to% the existing singleton*.% Z1(CALLBACK,hObj

4、ect,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 value pairs are% applied to the GUI before z1_OpeningFcn gets called. An% unrecognized

5、 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).% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to

6、 help z1% 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, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, z1_OpeningFcn, . gui_OutputFcn, z1_OutputFcn, . gui_LayoutFcn, , . gui_Callback, );if nargin

7、& ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);endif nargout varargout1:nargout = gui_mainfcn(gui_State, varargin:);else gui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before z1 is made visible.function z1_OpeningFcn(hObject, eventdata

8、, 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% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to z1 (see VARARGIN)% Choose default command

9、 line output for z1handles.output = hObject;% Update handles structureguidata(hObject, handles);axes(handles.axes3)map1=imread(1.bmp);imshow(map1)% UIWAIT makes z1 wait for user response (see UIRESUME)% uiwait(handles.figure1);% - Outputs from this function are returned to the command line.function

10、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 structure with handles and user data (see GUIDATA)% Get default command line ou

11、tput from handles structurevarargout1 = handles.output;% - Executes on button press in pushbuttonRun.function pushbuttonRun_Callback(hObject, eventdata, handles)%主要計(jì)算程序r2=0.1;%單位mr3=0.4;%單位momiga2=str2double(get(handles.edit1,String); %單位rad/sx11=1:500 %單位msfor i=1:500 theta2(i)=i*omiga2/1000; theta

12、3(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) %制表1plot(x11/1000,omiga3); xlabel(時(shí)間（t/s）) ylabel(連桿角速度3（rad/s）) axes(handles.axes2) %制表2plot(x11/1000,v3)

13、; xlabel(時(shí)間（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)function edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eve

14、ntdata 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,String) returns contents of edit1 as a double% - Executes during object creation, after sett

15、ing 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 created until after all CreateFcns called% Hint: edit controls usually have a white background o

16、n Windows.% See ISPC and COMPUTER.if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);end% - Executes on button press in pushbuttonExit.function pushbuttonExit_Callback(hObject, eventdata, handles)ss=questdlg(確認(rèn)退出？,退出信息窗口！,繼續(xù)仿真！,退出仿真！,退出仿真！);switch ss case 退出仿真！ delete(handles.figure1);end% hObject handle to pushbuttonExit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structu