MATLAB之理論力學課題作業(yè)

PAGEPAGE14理論力學項目專題研究題號：運動學1&動力學學院： 機械工程學院專業(yè)：機械設(shè)計制造及其自動班級： 機制171學號姓名： 顧繼偉日期： 2018-12-25目錄（目錄最好自動生成，要有頁碼）運動學1：曲柄連桿機構(gòu)運動方程、速度、加速度滑塊BPP橢圓推導(dǎo)動力學2：雙復(fù)擺三大定理建立運動微分方程動靜法建立運動微分方程運動規(guī)律研究參考文獻附錄理論力學A運動學項目專題研究（三選一）一、曲柄連桿機構(gòu)1OAOAB=l，P點（AP=e）的運動方B的運動方程、速度和加速度。2、利用MATLAB軟件畫出滑塊B的位置、速度、加速度圖隨轉(zhuǎn)角的變化規(guī)律，02。自行給出所需參數(shù)。3、利用MATLAB軟件畫出P（在連桿上至少選擇4個等距離的點并進行定性分析。自行給出所需參數(shù)。4、若用合成法計算連桿一點P的速度和加速度，不計算，給出解題思路。5、什么條件下？連桿上的點的運動軌跡是橢圓，請推導(dǎo)。運動方程、速度、加速度B位移、速度、加速度變化規(guī)律圖clear;clc;r=10;l=50;q=[0:0.005:2*pi];%%角度fai=asin(r.*sin(q)./l);clear;clc;r=10;l=50;q=[0:0.005:2*pi];%%角度fai=asin(r.*sin(q)./l);fore=1:4;%%P點位置fore=1:4;%%P點位置B點位移，速度，加速度變化圖B點位移，速度，加速度變化圖y=(l-e).*sin(fai);subplot(3,2,e);plot(x,y);holdon;end;y=(l-e).*sin(fai);subplot(3,2,e);plot(x,y);holdon;end;subplot(3,2,5);%%滑塊的數(shù)據(jù)plot(q,xb,'g');holdon;vb=diff(xb)/(q(2)-q(1));plot(q(1:end-1),vb,'r');holdon;ab=diff(vb)/(q(2)-q(1));subplot(3,2,5);%%滑塊的數(shù)據(jù)plot(q,xb,'g');holdon;vb=diff(xb)/(q(2)-q(1));plot(q(1:end-1),vb,'r');holdon;ab=diff(vb)/(q(2)-q(1));plot(q(1:end-2),ab,'b');P點運動軌跡圖r=10,l=50時，e=1，2,3,4時的圖P點速度、加速度橢圓推導(dǎo)二、雙復(fù)擺雙擺桿在豎直平面內(nèi)擺動，兩桿長均為l，質(zhì)量均為m，初始靜止，利用動力學普遍定理建立運動微分方程。利用動靜法建立運動微分方程。MATLAB律，并定性判斷結(jié)果正確與否？三大定理建立運動微分方程動靜法建立運動微分方程運動規(guī)律研究clearall;%????r1r2á???????ò?(???è)?í?é??±???°ú3?×′ì?1radμèóú57.?à??r1=0.6;r2=0.6;m1=1;m2=1;L1=1;L2=1;g=9.8;Da=inline(['[x(3);x(4);',...'inv([(m1+m2)*L1,m2*L2*cos(x(1)-x(2));',...'m1*L1*cos(x(1)-x(2)),m1*L2])*'...'[m2*L2*x(4)^2*sin(x(2)-x(1))-(m1+m2)*g*sin(x(1));',...'m2*L1*x(3)^2*sin(x(1)-x(2))-m2*g*sin(x(2))]]'],'t','x',...'flag','m1','m2','L1','L2','g');set(gcf,'DoubleBuffer','on');[t,x]=ode45(Da,[0,20],[r1,r2,0,0],[],m1,m2,L1,L2,g);axis([-(L1+L2),(L1+L2),-(L1+L2)*1.8,1]);axissquare;holdon;gh1=plot([0,L1*exp(i*(x(1)-pi/2))],'r-');set(gh1,'linewidth',2,'markersize',6,'marker','o');gh2=plot([L1*exp(i*(x(1)-pi/2)),L1*exp(i*(x(1)-pi/2))+L2*exp(i*(x(2)-pi/2))],'b-');set(gh2,'linewidth',2,'markersize',6,'marker','o');fork=2:size(x,1);C1=[0,L1*exp(i*(x(k,1)-pi/2))];C2=[L1*exp(i*(x(k,1)-pi/2)),L1*exp(i*(x(k,1)-pi/2))+L2*exp(i*(x(k,2)-pi/2))];','','','pause(0.1);endfigure;subplot(3,2,1);plot(t,x(:,1));title('t-\theta_1');xlabel('t');ylabel('\theta_1');subplot(3,2,2);plot(t,x(:,2));title('t-\theta_2');xlabel('t');ylabel('\theta_2');subplot(3,2,3);plot(t,x(:,3));title('t-\omega_1');xlabel('t');ylabel('\omega_1');subplot(3,2,4);plot(t,x(:,4));title('t-\omega_2');xlabel('t');ylabel('\omega_2');subplot(3,2,5);plot(x(:,1),x(:,3));title(