matlab二自由度系統(tǒng)振動(dòng)講解學(xué)習(xí)

1、Good is good, but better carries it.精益求精，善益求善。matlab二自由度系統(tǒng)振動(dòng)-利用Adams和Matlab對(duì)二自由度系統(tǒng)振動(dòng)進(jìn)行仿真與分析一、實(shí)驗(yàn)思想Adams是一種可以對(duì)一些典型運(yùn)動(dòng)進(jìn)行高效仿真的軟件，本實(shí)驗(yàn)是利用Adams對(duì)二自由度系統(tǒng)振動(dòng)進(jìn)行仿真及分析，再和理論公式對(duì)比，并用另外一種常見(jiàn)的仿真軟件Matlab的仿真結(jié)果進(jìn)行對(duì)比，觀察兩者的差異，分析軟件仿真產(chǎn)生差異的原因，加深對(duì)二自由度系統(tǒng)振動(dòng)的理解。二、二自由度系統(tǒng)振動(dòng)分析固有頻率取決于系統(tǒng)本身物理性質(zhì)，而與初始條件無(wú)關(guān)。對(duì)于二自由度的振動(dòng)系統(tǒng)是有兩種頻率的簡(jiǎn)諧波組成的復(fù)合運(yùn)動(dòng)，這兩個(gè)頻率都是

2、系統(tǒng)的固有頻率。主振型是當(dāng)系統(tǒng)按固有頻率作自由振動(dòng)時(shí)，稱為主振動(dòng)。系統(tǒng)作主振動(dòng)時(shí)，任何瞬時(shí)各個(gè)運(yùn)動(dòng)坐標(biāo)之間具有一定的相對(duì)比值，即整個(gè)系統(tǒng)具有確定的振動(dòng)形態(tài)，稱為主振型。強(qiáng)迫振動(dòng)是振動(dòng)系統(tǒng)在周期性的外力作用下，其所發(fā)生的振動(dòng)稱為強(qiáng)迫振動(dòng)，這個(gè)周期性的外力稱為驅(qū)動(dòng)力。三、二自由度系統(tǒng)自由振動(dòng)1.建立二自由度系統(tǒng)振動(dòng)模型1）創(chuàng)建底座：先生成一個(gè)尺寸合適的長(zhǎng)方體基體，再使用addtopart指令創(chuàng)建底座的側(cè)壁。2）使用newpart指令分別創(chuàng)建兩個(gè)滑塊，創(chuàng)建滑塊時(shí)應(yīng)注意滑塊與滑塊、滑塊與側(cè)壁之間的尺寸適當(dāng)。3）彈簧連接：分別用彈簧鏈接滑塊、側(cè)壁的中心點(diǎn)。彈簧生成后，依次選中彈簧，在modify選項(xiàng)中的

3、stiffnessanddamping下拉菜單中將dampingcoefficient設(shè)置成nodamping，即彈簧無(wú)阻尼。添加約束：底座和地面固定，滑塊和底座用滑動(dòng)副連接。彈簧剛度分別改為1、1、2(newton/mm)滑塊質(zhì)量分別為1.02.0滑塊與機(jī)體滑動(dòng)副的阻尼改為1.0E-007模型展示運(yùn)動(dòng)仿真結(jié)果設(shè)置x10=12經(jīng)過(guò)Adams運(yùn)算后，滑塊1、2運(yùn)動(dòng)狀態(tài)如圖所示：matlab驗(yàn)證程序：k1=1000;k2=1000;k3=2000;m1=1;m2=2;a=(k1+k2)/m1;b=k2/m1;c=k2/m2;d=(k2+k3)/m2;x1x2=dsolve(D2x1+2000*x1

4、-1000*x2=0,2*D2x2-1000*x1+3000*x2=0,x1(0)=0.012,x2(0)=0,Dx1(0)=0,Dx2(0)=0,t)t1=0:0.01:2;x1=subs(x1,t,t1);x2=subs(x2,t,t1);figureplot(t1,x1,-);title(系統(tǒng)響應(yīng)x(1)曲線);xlabel(時(shí)間/s);ylabel(位移/m);figureplot(t1,x2,-);title(系統(tǒng)響應(yīng)x（2）曲線);xlabel(時(shí)間/s);ylabel(位移/m);計(jì)算結(jié)果：5.結(jié)果分析存在差異的原因是Adams仿真中并沒(méi)有完全忽略摩擦力，而Matlab計(jì)算時(shí)沒(méi)有

5、考慮摩擦，故存在差異，但是在允許范圍內(nèi)。綜上所述，利用兩種軟件得出的結(jié)果輸出比較接近，可認(rèn)為仿真結(jié)果正確。二自由度系統(tǒng)受迫振動(dòng)將機(jī)體與地面的鎖改為滑動(dòng)副，同時(shí)將滑塊移動(dòng)副初始狀態(tài)設(shè)為0，即可2.運(yùn)動(dòng)結(jié)果仿真將底座和地面的滑動(dòng)副上輸入不同運(yùn)動(dòng)方程x=sin(w*t)。當(dāng)w=10時(shí)，得到滑塊1的一運(yùn)動(dòng)曲線；當(dāng)w為其固有頻率時(shí)，得到另一曲線。曲線如圖所示：參數(shù)設(shè)置仿真結(jié)果將w改為固有頻率31.6仿真結(jié)果3.實(shí)驗(yàn)結(jié)果檢驗(yàn)和單自由度系統(tǒng)一樣，二自由度系統(tǒng)在受到持續(xù)的激振力作用下就會(huì)產(chǎn)生強(qiáng)迫振動(dòng)，而且在一定條件下也會(huì)產(chǎn)生共振。共振是指一物理系統(tǒng)在必須特定頻率下，相比其他頻率以更大的振幅做振動(dòng)的情形；這些特

6、定頻率稱之為共振頻率。利用Matlab仿真得幅頻特性曲線，它表明系統(tǒng)位移對(duì)頻率的響應(yīng)特性。程序如下k1=1000;k2=1000;k3=2000;m1=1;m2=2;w=0:1:100;a=(k1+k2)/m1;b=k2/m1;c=k2/m2;d=(k2+k3)/m2;q1=1/m1;q2=1/m2;dw=w.4-(a+d)*(w.2)+a*d-b*c;B1=abs(d-w.2)*q1+b*q2)./dw);B2=abs(c*q1+(a-w.2)*q2)./dw);figureholdonplot(w,B1,g-);title(幅頻特性曲線);xlabel(./s);ylabel(B1、B2/m);plot(w,B2,r-);legend(B1,B2,1);圖像4.實(shí)驗(yàn)結(jié)果分析當(dāng)w=10時(shí)，相應(yīng)的質(zhì)量塊1幅值較小，僅有略微的振動(dòng)，但是當(dāng)w=31.6時(shí)，即共振時(shí)，其最終的幅值很大。由于在運(yùn)用Adams進(jìn)行