【正文】 18 對方程 ([K]－ ?2[M]){XN}={0}兩邊左乘 {XN(i)}T 可得到 ( ) ( ) 2{ } [ ] { }i T iN N iX K X ?? 注意： 這里的 {XN(i)}均為正規(guī)化后的振型 ， 而不是求解的原始主振型 {X(i)} 。 振型的基準(zhǔn)化和標(biāo)準(zhǔn)化 第 6章 多自由度系統(tǒng)的振動 17 2. 振型的標(biāo)準(zhǔn)化 另外一種確定振型各元素數(shù)值的方法是 ， 以某個限制條件來確定振型中的常數(shù)因子 。通常假設(shè)振型的某個元素為 1， 則其它元素就可以表示為此元素的倍數(shù) ， 這種方法或過程就是振型的基準(zhǔn)化 。 如在 A1格 “ 插入 ” “名稱 ” “定義 ” … w (2)輸入公式 。 多自由度系統(tǒng)的運(yùn)動微分方程式 第 6章 多自由度系統(tǒng)的振動 13 解：代入數(shù)值得 代入 |[K]?2[M]|=0得： 1 0 0[ ] 0 1 00 0 1M???????????2 1 0[ ] 1 2 10 1 1K?????? ? ??????6 4 25 6 1 0? ? ?? ? ? ? 理論求解很困難 ， 一般通過試算或利用工具軟件 ， 如 Excel、 MATLAB、 Mathematica等 。 2 ( )( [ ] [ ] ) { } { 0 } , ( 1 , 2 )iiK M X i n?? ? ?第 6章 多自由度系統(tǒng)的振動 12 例 T610中：設(shè) m1＝ m3＝ 1， m2＝ 2， r＝ 1， k1＝ k2＝ k3＝ 1。 顯然 : 和兩自由度一樣 ， 由上式只能求出振幅的比值 ， 而不能確定各振幅大小 。 由此可求出 n個特征根 （ 固有頻率 ） ?2。 作業(yè)： T69， 611 多自由度系統(tǒng)的運(yùn)動微分方程式 m1 m2 m3 m4 m5 k1 k2 k3 k4 k5 k8 k9 k10 k6 k7 第 6章 多自由度系統(tǒng)的振動 9 無阻尼自由振動的運(yùn)動方程為 主 振型 方程式 特征值和特征向量 [ ] { } [ ] { } { 0 }M x K x?? 利用兩自由度系統(tǒng)的分析結(jié)果 ， 假設(shè)方程解的形式為 ? ?{ } s i n ( )x X t????第 6章 多自由度系統(tǒng)的振動 10 代入振動方程得到廣義特征值問題： [K]?2[M]稱為 特征矩陣 。 4. [M]、 [K] 、 [C]均為對稱陣 。 2. [K]或 [C]中的主對角線元素 kii或 cii為連接在質(zhì)量 mi上所有彈簧剛度或阻尼系數(shù)的代數(shù)和 。 多自由度系統(tǒng)的運(yùn)動微分方程式 1 1 1 1 2 1()m x k x k x r?? ? ? ?22 2 1 3 31 ( ) ( )2 m r k x r r k x r r? ? ?? ? ? ?3 3 3 3()m x k x r?? ? ?解法 2：計算動勢能 。 多自由度系統(tǒng)的運(yùn)動微分方程式 多自由度系統(tǒng)的運(yùn)動微分方程式 第 6章 多自由度系統(tǒng)的振動 4 【 T610】 求系統(tǒng)的微振動微分方程 。 多自由度振動系統(tǒng)的很多概念和研究方法在兩自由度系統(tǒng)中已經(jīng)討論 。第 6章 多自由度系統(tǒng)的振動 1 第 6章 多自由度系統(tǒng)的振動 第 6章 多自由度系統(tǒng)的振動 2 多自由度系統(tǒng)指的是可以用有限個自由度描述的振動系統(tǒng) 。 一般來說 ， 一個 n自由度的振動系統(tǒng) ， 其廣義位移可以用 n個獨(dú)立坐標(biāo)來描述 ， 其運(yùn)動規(guī)律通?？捎?n個二階常微分方程來確定 。 第 6章 多自由度系統(tǒng)的振動 3 建立振動系統(tǒng)運(yùn)動微分方程的方法 ，包括一般的動力學(xué)方法 、 影響系數(shù)法 （ 剛度影響系數(shù)和柔度影響系數(shù) ） 、 拉格朗日方程和能量方法等 。 多自由度系統(tǒng)的運(yùn)動微分方程式 第 6章 多自由度系統(tǒng)的振動 5 解法 1：用動力學(xué)方法 。 2 2 2 21 1 3 3 21 1 1 12 2 2 2T m x m x m r ?? ? ?第 6章 多自由度系統(tǒng)的振動 6 多自由度系統(tǒng)的運(yùn)動微分方程式 1223001[ ] 0 0200mM m rm?????????????1 2 222 2 3 3330[ ] ( )0k k k rK k r k k r k rk r k??????? ? ? ??? ???2 2 21 1 2 1 3 31 1 1( ) ( )2 2 2V k x k x r k x r??? ? ? ? ?第 6章 多自由度系統(tǒng)的振動 7 對于有分支結(jié)構(gòu)的 mkc振動系統(tǒng) ， 可以用直觀目測方法直接形成振動系統(tǒng)的 [M]、 [K] 、 [C]: 1. [M]為各個質(zhì)量形成的對角陣 。 3. [K]或 [C]中的非對角線元素 kij或 cij為連接在質(zhì)量 mi 和 mj之間所有彈簧剛度或阻尼系數(shù)的串并聯(lián)等效剛度或阻尼 ， 且為負(fù)值 。 多自由度系統(tǒng)的運(yùn)動微分方程式 第 6章 多自由度系統(tǒng)的振動 8 例：用直觀目測方法直接形成標(biāo)準(zhǔn) mkc自由振動系統(tǒng)的 [M]、[K]。 要使 上 式有解 ，必須使其系數(shù)行列式為零 ： 2( [ ] [ ] ) { } { 0 }K M X???2[ ] [ ] 0KM??? 上式稱為 頻率方程 或 特征方程 。 第 6章 多自由度系統(tǒng)的振動 11 將每個特征根 ?i（ 固有頻率 ） 代入廣義特征值問題 ([K]－ ?2[M]){X}={0}, 可得到相應(yīng)的非零向量 {X(i)}, 稱為 特征矢量 ,或稱 特征向量 、固有振型 、 固有向量 、 模態(tài)向量 等 。 固有頻率和特征向量只決定于系統(tǒng)本身的物理特性 ， 而與外部激勵和初始條件無關(guān) ，它們都是系統(tǒng)的固有屬性 。 求固有頻率和振型 。 第 6章 多自由度系統(tǒng)的振動 14 利用 Excel計算固有頻率步驟： (1)定義變量 。 如在 A2格輸入 =w^35*w^2+6*w1 (3)“工具 ” “單變量求解 ” … (只能求第一固有頻率 ) (4)高階特征值的求解要用 “ 工具 ” “規(guī)劃求解 ” … 固有頻率為： 2 2 21 2 30 .1 9 8 , 1 .5 5 5 , 3 .2 4 7? ? ?? ? ?第 6章 多自由度系統(tǒng)的振動 15 分別代入 ([K]?2[M]){X}=0得： ( 1 )1{ } 1 .8 0 22 .2 4 7X????? ??????( 2 )1{ } 0 . 4 4 50 . 8 0 2X????? ???????( 3 )1{ } 1 . 2 4 70 . 5 5 5X????????????作業(yè)： T613 第 6章 多自由度系統(tǒng)的振動 16 由于固有振型 {X(i)} 只是振幅的比例關(guān)系 ，各階振型均有一個未確定的常數(shù)比例因子 。 一般假設(shè)振型的第一個元素為 1。 通常規(guī)定 {XN(i)}滿足條件 ( ) ( ){ } [ ] { } 1i T iNNX M X ? 滿足這個限制條件的振型 {XN(i)}稱為 標(biāo)準(zhǔn)化 （ 或正規(guī)化 、 歸一化 ） 的振型 。 第 6章 多自由度系統(tǒng)的振動 19 3. 標(biāo)準(zhǔn)化振型與主振型的關(guān)系 將主振型 {X(i)}進(jìn)行如下運(yùn)算： ( ) ( ){ } [ ] { }i T i iX M X M? Mi稱為廣義質(zhì)量 （ 主質(zhì)量 、 模態(tài)質(zhì)量 ） 。 例如給出 t＝ 0時的位移向量 {x0}和速度向量 {v0} ， 則得到含有 2n個方程的方程組 ()1{ } { } ( sin c os )nii i i iix X A t B t??????()01()01{ } { } si n{ } { } c osniiiinii i iix c Xv c X??????? ?????????()01()01{ } { }{ } { }niiiniiiix B Xv A X????? ?????????或 第 6章 多自由度系統(tǒng)的振動 22 【 T626】 圖示系統(tǒng)中 , m1＝ m2＝ m3＝ m, k1＝ k2＝ k3＝ k, 設(shè)初始位移為 1, 初始速度為 0, 求初始激勵的自由振動響應(yīng) 。 所以： ( ) ( ){ } [ ] { } 0j T iX M X ? （ i≠j） 主 振型的正交性 第 6章 多自由度系統(tǒng)的振動 30 這里的 Mi和 Ki是兩個實常數(shù) ， 分別稱為系統(tǒng)的 主質(zhì)量 和 主剛度 （ 或稱 模態(tài)質(zhì)量 和模態(tài)剛度 ） 。 即 [ ] [ ] [ ] [ ]00TpiM Q M Q M?????? ??????第 6章 多自由度系統(tǒng)的振動 33 [ ] [ ] [ ] [ ]00TpiK Q K Q K?????? ?????? 由主質(zhì)量矩陣 [Mp]和主剛度矩陣 [Kp]可得到