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

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

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

3、)(11)對(10)式對時間求導(dǎo)得：-r2 w2sin 0 2+3 w 3sin 0 3=v1r2 w2cos0 2-r3w3cos0 3=0 (12)將上式用矩陣形式表示，令：A= r3sin 0 31-r3cosO 30X= 口3v1B=-r2w2sin 0 2r2 w2cos 0 2則(12)可表水為：AX=R(13)從而可解出川3與丫1。三、MATLAB程序編寫源代碼如下：function varargout = z1(varargin)% Z1 MATLAB code for z1.fig%Z1, by itself, creates a new Z1 or raises the e

4、xisting% singleton*.% H = Z1 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 rais

5、es the% existing singleton*. Starting from the left, property 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 GUIDE's To

6、ols menu. Choose "GUI allows only one % instance to run (singleton)".% 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 EDITgui_Singleton = 1;gui_State = struct

7、('gui_Name',mfilename, .'gui_Singleton', gui_Singleton, .'gui_OpeningFcn', z1_OpeningFcn, .'gui_OutputFcn', z1_OutputFcn, .'gui_LayoutFcn', , .'gui_Callback', );if nargin && ischar(varargin1)gui_State.gui_Callback = str2func(varargin1);end if n

8、argoutvarargout1:nargout = gui_mainfcn(gui_State, varargin:);elsegui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before z1 is made visible.function z1_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hO

9、bject 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 line output for z1handles.output = hObject;% Update handles structureguid

10、ata(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 varargout = z1_OutputFcn(hObject, eventdata, handles) % varargo

11、ut 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 output from handles structure varargout1 = handles.output;% - E

12、xecutes on button press in pushbuttonRun.function pushbuttonRun_Callback(hObject, eventdata, handles)%主要計算程序r2=0.1;%單位 mr3=0.4;%單位 momiga2=str2double(get(handles.edit1,'String'); % 單位 rad/sx11=1:500 %單位msfor i=1:500theta2(i)=i*omiga2 /1000;theta3(i)=asin(-r2 /r3*sin(theta2(i);B=-r2*omiga2*si

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

14、;)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

15、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 c

16、reation, 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 created until after all CreateFcns called % Hint: edit controls usually have

17、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 button press in pushbuttonExit.function pushbuttonExit_C

18、allback(hObject, eventdata, handles)ss=questdlg('確認(rèn)退出？,退出信息窗口！,繼續(xù)仿真！,'退出仿真！,退出仿真！);switch sscase退出仿真！'delete(handles.figure1);end% hObject handle to pushbuttonExit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)四、使用指南和實例仿真進(jìn)入MATLA瞰件，打開并運行程序初始界面:對 應(yīng)的值進(jìn)行修改，修改為實例中的50*2*pi=314.16 ,點擊開始仿真按鈕，得到仿真結(jié)果:Tzi制作考：自動化148班前鞋飛曲柄滑塊機(jī)構(gòu)的運動學(xué)仿真曲柄滑塊札怕示意圖.2=妣什 w杷朋凱ft vl-drl/dljr2=0.1rrir3=0.4m0)2=314.16 rad，s仿真時間0 . 55汨薪真!,退出