第三章 多自由度系統(tǒng)振動6.19

1、第三章 多自由度系統(tǒng)振動多自由度系統(tǒng)和單自由度系統(tǒng)的振動特性是有區(qū)別的。單自由度系統(tǒng)受初始擾動后，按系統(tǒng)的固有頻率作簡諧振動。多自由度系統(tǒng)有多個固有頻率，當(dāng)系統(tǒng)按某一個固有頻率作自由振動時，各獨立坐標(biāo)在振動過程中相互關(guān)系是固定的，這個關(guān)系叫振幅比，也叫作主振型或模態(tài)。主振型是多自由度系統(tǒng)以及彈性體振動的重要特征。多自由度系統(tǒng)的振動方程是多個二階微分方程組，這些方程一般是耦合的。多自由度振動的求解有兩種方法：直接積分法和振型疊加法。直接積分法可直接根據(jù)微分方程求出響應(yīng)，涉及的概念不多且有應(yīng)用軟件，本章不做介紹。振形疊加法要先求出系統(tǒng)的固有頻率和振型，在此基礎(chǔ)用疊加法求響應(yīng)，物理概念清楚、并且是模

2、態(tài)分析與參數(shù)識別的理論基礎(chǔ)。因此本章將先用較多的篇幅介紹多自由度系統(tǒng)的固有振動特性、振型疊加法和傳遞函數(shù)。3.1 振動微分方程雖然一些多自由度系統(tǒng)數(shù)目較多，有些相當(dāng)復(fù)雜，但建立多自由度系統(tǒng)振動微分方程并沒有新理論和方法，都是動力學(xué)基本理論和方法，本節(jié)只通過例題介紹多自由度系統(tǒng)振動微分方程基本形式。例一 試建立圖3-1所示3自由度系統(tǒng)的運動微分方程。三個質(zhì)量只作水平方向的運動，并分別受到激振力，和的作用，質(zhì)量塊的質(zhì)量分別為，和，彈簧剛度分別為， 和，阻尼分別為， 和。圖3-1 3自由度系統(tǒng)解：分別用三個獨立坐標(biāo)，和描述三個質(zhì)量塊的運動，坐標(biāo)原點分別取在，和的靜平衡位置。質(zhì)量塊的速度分別為，和，加

3、速度分別為，和。每個質(zhì)量塊的受力圖如3-2（a、b、c）所示，則由受力圖根據(jù)牛頓第二定律，得系統(tǒng)的運動方程為： 圖3-2 (a) 圖3-2(b)圖3-2(c)或上述方程組可以用矩陣表示為：三個二階微分方程是耦合的，這是因為矩陣中有非零的非對角元素。若質(zhì)量、剛度和阻尼矩陣都是對角矩陣，則三個微分方程是獨立的，相當(dāng)于三個獨立的單自由度系統(tǒng)，其求解變?yōu)槿齻€單自由度系統(tǒng)求解。質(zhì)量矩陣中出現(xiàn)耦合項稱為慣性耦合，剛度矩陣中出現(xiàn)耦合項稱為彈性耦合。此例中沒有慣性耦合，因為質(zhì)量矩是對角的。但一般情況下質(zhì)量矩陣并不是對角的，所以一般情況下多自由度系統(tǒng)既有彈性耦合、也有慣性耦合。下面我們通過一個例子來說明質(zhì)量矩陣

4、不是對角的情況。例二 寫出圖3-3所示系統(tǒng)振動微分方程系統(tǒng)中均質(zhì)剛性桿AB的質(zhì)量為m，轉(zhuǎn)動慣量為，前后兩端分別用剛度為和的兩個彈簧由承于地面上，桿全為長。圖3-3若用桿兩端的豎向位移、來描述剛桿的運動狀態(tài)，則受力圖如圖3.4所示， 圖3.4顯然、質(zhì)心處的加速度為，根據(jù)牛頓第二定律，在豎直方向有：桿的轉(zhuǎn)動加速度為（順時針為正），對C點應(yīng)用動力矩定理：整理并寫成矩陣形式有：質(zhì)量矩陣并不是對角的。當(dāng)然，此例中若選質(zhì)心的平動及繞質(zhì)心的轉(zhuǎn)動來描述運動，質(zhì)量矩陣將是對角的。一般地，對n自由度系統(tǒng)，振動微分方程為：寫成矩陣形式有： (3.1)根據(jù)分析力學(xué)，具有定常約束的系統(tǒng)的動能T與勢能U可寫為下列二次型

5、(3.2)對于穩(wěn)定平衡的振動系統(tǒng),系統(tǒng)的動能T總是大于零的(除非系統(tǒng)是靜止的)，所以質(zhì)量矩陣一般是正定的。同樣，系統(tǒng)的勢能U也總大于零，所以剛度矩陣也是正定的。此外，系統(tǒng)的動能和勢能不會因為表達形式不同而改變，對式(3.2)轉(zhuǎn)置，比較可知，剛度矩陣和質(zhì)量矩陣必須是對稱矩陣，因而有： (3.3)3.2無阻尼自由振動一、固有頻率和振形本節(jié)主要目的是通過無阻尼自由振動系統(tǒng)來介紹多自由系統(tǒng)的固有頻率和振型，它們是多自由振動系統(tǒng)的重要特征。在無阻尼情況下，系統(tǒng)的自由振動微分方程可以表達為： (3.4)在單自由度系統(tǒng)中,我們得到無阻尼自由振動解為正弦函數(shù)或余弦函數(shù),不失一般性。對于多自由度系統(tǒng)振動解可設(shè)為

6、： (3.5)列向量和均為待定復(fù)常數(shù)。若系統(tǒng)是振動的，則解必為實數(shù)。將式(3.5)代入(3.4)，得到下列代數(shù)齊次方程組： (3.6)上面的方程組存在非零解的充分必要條件是系數(shù)行列式為零，即： (3.7)式(3.7)為系統(tǒng)的特征方程，具體寫出為： (3.8)上式左端的行列式展開后是關(guān)于的n次代數(shù)多項式： (3.9)稱為特征多項式，由式(3.8)或（3.9）可解出n個稱為特征值或特征根，將其按升序排列為：顯然特征值僅取決于系統(tǒng)本身的剛度和質(zhì)量參數(shù)。這n個特征值在大多數(shù)情況下互不相等且不為零，重根的零根說明系統(tǒng)有剛體運動。有零根和情況本書不再討論，有興趣的讀者可參考相關(guān)的線性代數(shù)和振動理論書籍。在

7、求得特征值后把某一個代回式(3.6)，可求對應(yīng)的列向量。由于式(3.6)的系數(shù)矩陣不滿秩，在沒有重根和零根情況下只有(n-1)個是獨立的，故只能求出列向量中各元素、的比例關(guān)系。我們?nèi)サ羝渲胁华毩⒌哪骋皇?例如最后一式)，并將剩下的n-1個方程式中某一相同的項(如項)移到等式右邊，可得代數(shù)方程組： (3.10)解上面的方程，可得到用表達的解、，顯然都與的值成比例。我們可將這些比例常數(shù)用表示，并補充，可得列向量，則有： (3.11)列向量是確定的常數(shù)，反映列向量中各數(shù)的比例關(guān)系，叫作特征向量。同比例放大或減小特征向量并不改變其比例關(guān)系，所以應(yīng)用時常根據(jù)需要來放大或減小特征向量。不失一般性，我們可在

8、式（3.11）中用待定復(fù)常數(shù)取代，式（3.11）可寫為： (3.12)這樣，當(dāng)成比例變化時，有相應(yīng)的變化，對應(yīng)不同的特征值，可得到不同的特征向量。對應(yīng)于n個特征值可得n個特征向量 ，且每一個特征向量都滿足式（3.6）。對于一個振動系統(tǒng)，特征值就是系統(tǒng)的固有頻率，特征值相對應(yīng)的特征向量就是系統(tǒng)的振形。顯然，對應(yīng)于n個固有頻率可得n個振形。我們將在后面論述。二、無阻尼自由振動的解顯然，將及代入式(3.5)，可得n組滿足方程（3.4）的解，將這些解相加，可得多自由度系統(tǒng)自由振動的一般解為： (3.13)其中2n個待定常數(shù)由系統(tǒng)運動的初始位移和初始速度確定。如果系統(tǒng)在某一特殊的初始條件下，使得待定常數(shù)

9、中只有0，則式(3.13)所表示的系統(tǒng)運動方程只保留第k項: (3.14)參見前一章，多自由度系統(tǒng)振動一般解的方程可表達為： (3.15)這時整個系統(tǒng)按圓頻率、振幅比作同步簡諧運動。振幅分別為，振幅之間都保持固定不變的比值。因此特征向量完全確定了系統(tǒng)按固有頻率振動時的形態(tài)，所以特征向量就是按相應(yīng)固有頻率振動時的振型向量，對應(yīng)的特征向量稱為它的第階主振型或主模態(tài)，相應(yīng)的振動叫主振動。在振動過程中，一般還會產(chǎn)生其它階主振動。對于一個n自由度系統(tǒng)，一般可以找到n個固有頻率，以及相應(yīng)的n個主振型。我們把各階主振型組成的矩陣叫做振型矩陣： (3.16)例三在圖3-5所示的三自度系統(tǒng)中，設(shè)，,，求系統(tǒng)的固

10、有頻率、振型。圖3-5 解：分別用3個獨立坐標(biāo)、和描述三個質(zhì)量塊的水平運動，可寫出系統(tǒng)的質(zhì)量矩陣和剛度矩陣。 ； (1)系統(tǒng)自由振動微分方程為： (2)令 代入上式整理可得特征值方程為： (3)展開整理后，有： (4)可求得的三個根為：所以系統(tǒng)的固有頻率分別為： 不妨設(shè)所對應(yīng)的列向量為，且=那么解特征向量的代數(shù)方程為：將=代人(上)式中二式有，且：對應(yīng)的列向量為=解得:一階振型,見圖3.6（a）同理，將=代人(上)式中第式有：解得: 二階振型，見圖3.6（b）將=代人(上)式中二式有：解得: 三階振型，見圖3.6（c）其各階陣形如圖所示：圖3.6（a）圖3.6（b）圖3.6（c）3.3振型的正

11、交性在前面我們分析了無阻尼系統(tǒng)自由振動的一般性質(zhì)，指出一個n自由度的系統(tǒng)具有n個固有頻率及n個主振型 (j=1,2,n)，本節(jié)主要討論主振型之間關(guān)于質(zhì)量矩陣和剛度矩陣的正交性。設(shè)、分別是相應(yīng)于固有頻率、的主振型，由式（3.6）有： (3.17)同理可得： (3.18)將式（3.17）和式（3.18）分別左乘和得： (3.19) (3.20)將式（3.19）轉(zhuǎn)置，考慮到質(zhì)量和剛度矩陣的對稱性有： (3.21)將式(3.21)減去式(3.20)得： (3.22)當(dāng)不等于時，若無重根，,則有： （ij） (3.23)將式（3.23）回代到式 (3.20)得到： （ij） (3.24)上兩式(3.23

12、)、式(3.24)表示固有頻率不相等的兩個振型之間，存在著對質(zhì)量矩陣和和剛度矩陣的正交性。當(dāng)i等于j時，可記為： (3.25) (3.26)和分別叫作第階振動的主質(zhì)量和主剛度，都是正實數(shù)，顯然： (3.27)推廣到整個振型矩陣有： (3.28)其中為對角陣，稱為主質(zhì)量矩陣，即： (3.29) 也是對角陣，稱為主剛度矩陣，即： (3.30)在近似計算或?qū)嶒灧治鰰r，給出相應(yīng)的振型函數(shù)后，就可得到系統(tǒng)固有頻率的近似值，其精度與振型函數(shù)的選取有關(guān)。如果主質(zhì)量為1，則主剛度就是固有頻率的平方（式3.27），使表達式更簡潔。使主質(zhì)量等于1的振型叫正則振型。由于主振型列陣只表示系統(tǒng)作主振動時各幅值的比例關(guān)系

13、，同除主質(zhì)量平方根并不影響這個比例關(guān)系。各主振型正則化的過程就是求出主振型后，再求出主質(zhì)量，然后將主振型除以主質(zhì)量的平方根。若記為正則振型，要使主質(zhì)量為1，只需令。顯然，振形正則化后的主質(zhì)量矩陣是單位矩陣?！纠摹壳罄鞠到y(tǒng)的主剛度矩陣、主質(zhì)量矩陣、正則振型矩陣以及振形正則化后的主剛度矩陣、主質(zhì)量矩陣。解：例三中已得出質(zhì)量矩陣、剛度矩陣和各階振型，振形矩陣為：從而主質(zhì)量矩陣主剛度矩陣由可得各正則振型列陣如下于是得正則振型矩陣主剛度矩陣3.4振型的疊加法在無阻尼自由振動中，我們得到： (3.31)從式中看出，每一個獨立坐標(biāo)的振動量是各階振型的疊加，這啟發(fā)我們用各階主振動來描述系統(tǒng)的振動，這種

14、方法叫振型疊加法。振型疊加法求解多自由度系統(tǒng)就是將按獨立坐標(biāo)振動的問題轉(zhuǎn)變?yōu)榘粗髡裥驼駝拥膯栴}。用獨立坐標(biāo)（）來描述多自由度系統(tǒng)的一般振動微分方程為： (3.32) 在求得系統(tǒng)主振型后，設(shè)描述各階主振型的振動量的獨立坐標(biāo)為,則各獨立坐標(biāo)的振動是各主階振動量的疊加，可設(shè)： (3.33)為按第主振型振動的坐標(biāo)，叫作主坐標(biāo)。在數(shù)學(xué)上，這相當(dāng)于做了一個變換。經(jīng)過這個變換后，相應(yīng)的初始條件為： (3.34) (3.35)將其式(3.33)代入式(3.32)，并前乘振型矩陣的轉(zhuǎn)置得： (3.36) 將式（3.25），（3.26）代入上式（3.31），并記, (3.37a) (3.37b)可得： (3.38

15、)雖然主質(zhì)量矩陣與主剛度矩陣是對角陣，矩陣一般并非對角陣，所以式(3.38) 是一組通過速度項相互耦合的微分方程。如果也是一個對角矩陣，則式(3.38)是不耦合的，是n個獨立的二階微分方程，相當(dāng)于n個單自由度系統(tǒng)，使求解大為簡化。當(dāng)阻尼較小，略去的非對角線元素組成的各阻尼項，即令的非對角線元素的值為零可使方程解耦；另外一種方法是假設(shè)阻尼矩陣是質(zhì)量矩陣和剛度矩陣的線性組合的情況,使阻尼矩陣對角化。即： (3.39) 解耦后式(3.38)為： (3.40)轉(zhuǎn)化為n個單自由度系統(tǒng)，可根據(jù)初始條件求解，得。 最后，將方程代入坐標(biāo)變換式(3.33)，即可得到廣義坐標(biāo)下的響應(yīng)，即 (3.41)對于大阻尼情

16、況，引入復(fù)模態(tài)理論可使方程解耦，本書不做討論。 例五：用振型疊加法求解上例題例三初始速度及初始位移為下系統(tǒng)的響應(yīng)。解：由上例可知：無阻尼自由振動微分方程為： (1)振型函數(shù)為： (2)主質(zhì)量矩陣：(3)主剛度矩陣：(4) 將代入(1)，并左乘可知： (5)展開得到： (6)解(6)式得到： (7)其中： 根據(jù)初始條件，由可知，主坐標(biāo)就滿足的初始條件，。對振形矩陣求逆得：從而，對式進行求導(dǎo)得到：將初始條件，代入式和，得到待定系數(shù)的解：同理，得到： 代入式，得到：代入中，得到：3.5多自由度傳遞函數(shù)多自由度振動微分方程為： (3.42)這是多輸入多輸出問題，按傳遞函數(shù)的定義，每一個輸入對每一個輸出

17、都有一個傳遞關(guān)系，若輸入，其余均為零，則可得j點輸入時對各輸出點的傳遞函數(shù)，對全部輸入可得傳遞函數(shù)矩陣 (3.43)事實上對于復(fù)數(shù)形式激勵，設(shè)系統(tǒng)的穩(wěn)態(tài)響應(yīng)的復(fù)數(shù)形式 ，代入微分方程 (3.42)得： (3.44) 其中，為傳遞函數(shù)， (3.45)顯然，由于系統(tǒng)參數(shù)是對稱矩陣，傳遞函數(shù)矩陣也是對稱的。若輸入為單位脈沖，其余均為零，則可得j點輸入時對各輸出點的單位脈沖響應(yīng)，對全部單位脈沖激勵可得單位脈沖響應(yīng)矩陣 (3.46)利用上式作傅氏變換，可得到傳遞函數(shù)矩陣。為了更好的分析傳遞函數(shù)的性質(zhì),設(shè)，代入式（3.38）解耦后得到： (3.47)當(dāng)激勵為時，主坐標(biāo)的穩(wěn)態(tài)響應(yīng)有修改解為： (3.48)激

18、勵,代入上式（3.28）得到： (3.49)激勵對主坐標(biāo)的傳遞函數(shù)為： (3.50)若變換中振型矩陣式為正則振型矩陣，則上式為： (3.51)因為，所以代入得每一項 (3.51)下標(biāo)表示第點激勵，第點響應(yīng)之間的傳遞關(guān)系。顯然，傳遞函數(shù)矩陣是對稱的，是項多項式的和，其中每一項相都當(dāng)于一個單自由度系統(tǒng)的傳遞函數(shù)。上式分離后可得實頻和虛頻關(guān)系和，從而其幅頻和相頻分別如下： 幅值： (3.52)相位： (3.53)傳遞函數(shù)反映了系統(tǒng)參數(shù)與固有頻率、振形之間的關(guān)系，是振動測試和分析的理論基礎(chǔ)。 例六 在圖3-7所示的三自度系統(tǒng)中，設(shè)，求傳遞函數(shù)，并圖示其實頻、虛頻、幅頻和相頻特性。圖3.7解：第三節(jié)已經(jīng)求出系統(tǒng)的固有頻率為；正則振形為： 將已知條件和上述結(jié)果代入式（3.51）可得：我把式（3.51）重新推導(dǎo)了一下，發(fā)現(xiàn)原來求和下標(biāo)錯了，已