單自由度系統(tǒng)的振動及matlab分析

1、精選優(yōu)質文檔-傾情為你奉上單自由度系統(tǒng)的振動及matlab分析 摘要：以彈簧質量系統(tǒng)為力學模型,研究單自由度系統(tǒng)的特性有著非常普遍的實際意義。根據(jù)單自由度振動系統(tǒng)數(shù)學模型,利用Matlab軟件設計了單自由度振動系統(tǒng)的數(shù)學仿真實驗。通過實驗可以得到單自由度振動方程的數(shù)值關鍵字：有阻尼自由振動、有阻尼自由振動、matlab正文：無阻尼自由振動:如圖所示的單自由度振動系統(tǒng)可以用如下微分方程描述：圖1 (1-1)令 ，方程的通解為 (1-2)式(1-2)表示了圖示(1)中質量m的位置隨時間而變化的函數(shù)關系，反映了振動的形式與特點，稱為振動函數(shù)。式(1-2)中，a、b為積分常數(shù)，它決定于振動的初始條件。

2、如假定t=0時，質量塊的位移 x=x0，其速度 ,則即 （1-3）或寫成 （1-4），其中A為振幅，為振動圓頻率， 為相位角，（赫茲）稱為固有頻率。固有頻率與外界給予的初始條件無關，它是系統(tǒng)本身所具有的一種重要特性。有阻尼自由振動圖2圖1所示的自由振動中，由于系統(tǒng)的能量守恒，如果振動一旦發(fā)生，它就會持久的，等幅的一直進行下去。但是，實際上所遇到的自由振動都是逐漸衰減而至最終停止的，即系統(tǒng)存在阻尼。阻尼有相對運動表面的摩擦力，液體與氣體的介質阻力，電磁阻力以及材料變形時的內阻力等。圖2所示為考慮了阻尼的單自由度振動系統(tǒng)模型。其運動微分方程為 (2-1)令，則 （2-2）其通解為 （2-3）式中c

3、1、c2為積分常數(shù)，由振動初始條件確定。令 ，稱為相對阻尼系數(shù)或阻尼率。則式 （2-3）可寫為 （2-4）由此可以討論阻尼對系統(tǒng)的自由振動產生的影響。一、當時，稱為弱阻尼狀態(tài)此時，為虛數(shù)，式 (2-4) 變?yōu)?（2-5）利用歐拉公式，式(2-5)可寫為 （2-6）括號內為兩個簡諧振動相加,即式 (2-5) 可寫為 （2-7）由式(2-7)可以看出，弱阻尼自由振動具有如下幾種特性:它是一個簡諧振動，振動的頻率為，這是為無阻尼時系統(tǒng)的固有頻率。一般情況下，常在0.1左右，因此對固有頻率的影響不大，即認為 。2. 振動的振幅為，其中、皆為定值。所以振幅隨時間變化的規(guī)律是一條指數(shù)遞減曲線(圖3)。二、

4、當時，稱為強阻尼狀態(tài)此時，式(1.2.2-4)可寫成圖3 （2-8）由于，故式(2-8 )中二項指數(shù)皆為實數(shù)。又因為，故二項之指數(shù)皆為負值，所以，式（2-8）所表示的是一根指數(shù)遞減曲線。這表示系統(tǒng)將不再產生前面所述的振動，而是產生一按指數(shù)規(guī)律衰減的曲線。三、當時，稱為臨界阻尼狀態(tài) 由于 ，則有 (2-9) 這里 cc 為臨界阻尼狀態(tài)下的阻尼系數(shù)，稱為臨界阻尼系數(shù)。顯然它是系統(tǒng)本身所具有的特性之一。由及，有。也就是說，相對阻尼系數(shù)（阻尼率）反映了系統(tǒng)的實際阻尼與臨界阻尼的關系。在臨界阻尼狀態(tài)下，有 (2-10)其中。顯然，在這種狀態(tài)下不能形成振動。有阻尼自由振動響應計算與 MATLAB實現(xiàn)根據(jù)式

5、(2-7)、(2-8)、(2-10) 編寫的程序如下：function VTB1(m,c,k,x0,v0,tf)%VTB1用來計算單自由度有阻尼自由振動系統(tǒng)的響應%VTB1繪出單自由度有阻尼自由振動系統(tǒng)的響應圖%m為質量;c為阻尼;k為剛度;x0為初始位移;v0為初始速度;tf為仿真時間%VTB1(zeta,w,x0,v0,tf)繪出單自由度有阻尼自由振動系統(tǒng)的響應圖%zeta為阻尼系數(shù);n為固有頻率%程序中z為阻尼系數(shù);A為振動幅度;phi為初相位clc%該循環(huán)確定輸入方式是VTB1(m,c,k,x0,v0,tf),還是%VTB1(zeta,w,x0,v0,tf)if nargin=5 z=

6、m;wn=c;tf=v0;v0=x0;x0=k;m=1;c=2*z*w;k=w2;endwn=sqrt(k/m);%固有頻率z=c/2/m/wn;wd=wn*sqrt(1-z2);fprintf('固有頻率為%.3g.rad/s.n',wn);fprintf('阻尼系數(shù)%.3g.n',z);fprintf('有阻尼的固有頻率%.3g.n',wd);t=0:tf/1000:tf;if z<1 A=sqrt(v0+z*wn*x0)2+(x0*wd)2)/wd2); phi=atan2(x0*wd,v0+z*wn*x0); x=A*exp(-z*

7、wn*t).*sin(wd*t+phi); fprintf('A=%.3gn',A); fprintf('phi=%.3gn',phi);elseif z=1 a1=x0; a2=v0+wn*x0; fprintf('a1=%.3gn',a1); fprintf('a2=%.3gn',a2); x=(a1+a2*t).*exp(-wn*t);else a1=(-v0+(-z+sqrt(z2-1)*wn*x0*wn*x0)/2/wn/sqrt(z2-1); a2=(v0+(z+sqrt(z2-1)*wn*x0)/2/wn/sqrt(

8、z2-1);圖4 fprintf('a1=%.3gn',a1); fprint('a2=%.3gn',a2); x=exp(-wn*t).*(a1*exp(wn*sqrt(z2-1)*t) +a2*exp(wn*sqrt(z2- 1)*t);endplot(t,x),gridxlabel('時間(s)')ylabel('位移')title('位移相對時間的關系')運行該程序時，只需要給出相應參數(shù)，例如>>VTB1(1,0.05,1,1,1,100)則顯示固有頻率為1（rad/s）,阻尼系數(shù)0.03，幅值為A=1.43相位角為phi=0.773。其響應曲線如圖 (4) 所示。程序中if語句就是判斷大小的，即判斷是弱阻尼狀態(tài)、強阻尼狀態(tài)還是臨界阻尼狀態(tài)。如果運行>>VTB1(1,2,1,0.1,1,20)，則顯示固有頻率為1（rad/s），阻尼系數(shù)1。其響應曲線如圖 (5) 所示。如果要想求出振動的速度(=xd)和加速度(=xdd)，只要式(2-1)、式(