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

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

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

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

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

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

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

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

8、式（3.11）中用待定復常數(shù)取代，式（3.11）可寫為： (3.12)這樣，當成比例變化時，有相應的變化，對應不同的特征值，可得到不同的特征向量。對應于n個特征值可得n個特征向量 ，且每一個特征向量都滿足式（3.6）。對于一個振動系統(tǒng)，特征值就是系統(tǒng)的固有頻率，特征值相對應的特征向量就是系統(tǒ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)，所以特征向量就是按相應固有頻率振動時的振型向量，對應的特征向量稱為它的第階主振型或主模態(tài)，相應的振動叫主振動。在振動過程中，一般還會產(chǎn)生其它階主振動。對于一個n自由度系統(tǒng)，一般可以找到n個固有頻率，以及相應的n個主振型。我們把各階主振型組成的矩陣叫做振型矩陣： (3.16)例三在圖3-5所示的三自度系統(tǒng)中，設，,，求系統(tǒng)的固

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

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

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

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

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

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

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

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

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