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

PAGEPAGE14理論力學(xué)項(xiàng)目專題研究題號(hào)：運(yùn)動(dòng)學(xué)1&動(dòng)力學(xué)學(xué)院： 機(jī)械工程學(xué)院專業(yè)：機(jī)械設(shè)計(jì)制造及其自動(dòng)班級(jí)： 機(jī)制171學(xué)號(hào)姓名： 顧繼偉日期： 2018-12-25目錄（目錄最好自動(dòng)生成，要有頁(yè)碼）運(yùn)動(dòng)學(xué)1：曲柄連桿機(jī)構(gòu)運(yùn)動(dòng)方程、速度、加速度滑塊BPP橢圓推導(dǎo)動(dòng)力學(xué)2：雙復(fù)擺三大定理建立運(yùn)動(dòng)微分方程動(dòng)靜法建立運(yùn)動(dòng)微分方程運(yùn)動(dòng)規(guī)律研究參考文獻(xiàn)附錄理論力學(xué)A運(yùn)動(dòng)學(xué)項(xiàng)目專題研究（三選一）一、曲柄連桿機(jī)構(gòu)1OAOAB=l，P點(diǎn)（AP=e）的運(yùn)動(dòng)方B的運(yùn)動(dòng)方程、速度和加速度。2、利用MATLAB軟件畫出滑塊B的位置、速度、加速度圖隨轉(zhuǎn)角的變化規(guī)律，02。自行給出所需參數(shù)。3、利用MATLAB軟件畫出P（在連桿上至少選擇4個(gè)等距離的點(diǎn)并進(jìn)行定性分析。自行給出所需參數(shù)。4、若用合成法計(jì)算連桿一點(diǎn)P的速度和加速度，不計(jì)算，給出解題思路。5、什么條件下？連桿上的點(diǎn)的運(yùn)動(dòng)軌跡是橢圓，請(qǐng)推導(dǎo)。運(yùn)動(dòng)方程、速度、加速度B位移、速度、加速度變化規(guī)律圖clear;clc;r=10;l=50;q=[0:0.005:2*pi];%%角度f(wàn)ai=asin(r.*sin(q)./l);clear;clc;r=10;l=50;q=[0:0.005:2*pi];%%角度f(wàn)ai=asin(r.*sin(q)./l);fore=1:4;%%P點(diǎn)位置fore=1:4;%%P點(diǎn)位置B點(diǎn)位移，速度，加速度變化圖B點(diǎn)位移，速度，加速度變化圖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點(diǎn)運(yùn)動(dòng)軌跡圖r=10,l=50時(shí)，e=1，2,3,4時(shí)的圖P點(diǎn)速度、加速度橢圓推導(dǎo)二、雙復(fù)擺雙擺桿在豎直平面內(nèi)擺動(dòng)，兩桿長(zhǎng)均為l，質(zhì)量均為m，初始靜止，利用動(dòng)力學(xué)普遍定理建立運(yùn)動(dòng)微分方程。利用動(dòng)靜法建立運(yùn)動(dòng)微分方程。MATLAB律，并定性判斷結(jié)果正確與否？三大定理建立運(yùn)動(dòng)微分方程動(dòng)靜法建立運(yùn)動(dòng)微分方程運(yùn)動(dòng)規(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(