1、多自由度系统的振动、响应和求解4.1 多自由度系统的动力学方程多自由度系统的动力学方程4.2 无阻尼系统的自由振动与模态无阻尼系统的自由振动与模态4.3 振动分析的模态叠加法振动分析的模态叠加法4.4 无阻尼系统的响应无阻尼系统的响应4.5 阻尼系统及其求解阻尼系统及其求解4.6 模态问题的一些特殊情况模态问题的一些特殊情况4.7 传递矩阵法传递矩阵法第四章 多自由度系统的振动大部分实际系统都是多自由度系统，其中的一类，系统本身为近似的集中参数系统，可以简化为多自由度系统，另一类是将分布参数系统通过一定的建模方法简化得到的。本章只学习线性多自由度系统的分析方法和基本规律，解决问题的基本方法是模
2、态叠加法，就是将n自由度系统分解成 n 个单自由度系统，每个单自由度系统对应于原系统的一种特定的振动形态（即模态），将各个单自由度系统的振动叠加便得到原系统的振动。因此，本章的学习重点是要理解和掌握模态的求解和使用。4.1 多自由度系统的动力学方程我们先来考察多自由度线性系统动能和势能的数学结构。对n个自由度的完整系统，设系统的广义坐标为qi , i = 1, 2, , n；系统的势能为V(q1, q2, , qn)，用Taylor 级数展开nijinjqqjiniiqqiqqnndqdqqqVdqqVqqqVqqqViiiiii11212121)(21)(),(),(000将广义坐标原点和零
3、势位均取在系统的静平衡位置，得nijinjqjinqqqqVqqqVi110221)(21),(4.1)njiqqVkqqkqqqVninjqjiijnijinjijni, 2, 1,21),(11021121其中量，得将上式只保留到二阶小因为偏导数的结果与求导次序无关，因此 kij = kji 。对于线性系统，(4.2)式是精确的；对于非线性系统，上式近似成立，可用来研究系统在平衡位置附近的微振动、或用来逐步逼近非线性系统。设系统受到定常约束，则动能为广义速度的二次齐函数(4.2)(4.3)nijinjjkikNkkniiikNkkkNkkqqqrqrmqqrmrmT111121212121
4、21nnijnnijTnTTmMkKqqqqqMqTKqqV,21,2121其中jijkikNkkijnijinjijnijinjjkikNkkmqrqrmmqqmqqqrqrm 1111112121其中(4.2)、(4.4)式可写成矩阵形式(4.4)(4.5)(4.6)矩阵K 称为刚度矩阵，它是一个对称正定或半正定矩阵；矩阵M 称为质量矩阵，它是一个对称正定矩阵。根据Lagrange方程，系统的振动微分方程为0M qKqQQM qKq其中为非保守力形成的广义力列阵对自由振动系统，振动方程为(4.7)(4.8)下面解释一下刚度矩阵的物理意义。由（4.7）式，系统的静力平衡方程为iijjmjjn
5、jjjQkQkjmqqjKkQqkQKq,)(0, 1,1时当列的第为矩阵可写成(4.10)上式的物理意义是：刚度矩阵的第j列是使qj产生单位变形所需的广义力列阵。可以应用这一物理含义直接列写系统的刚度矩阵，称为刚度影响系数法。(4.9)刚度矩阵的逆矩阵F 称为柔度矩阵qfjmQQQfFQQKqjmjnjjj时当, 0, 111柔度矩阵的第j列是Qj为单位力而使系统产生的广义坐标变形列阵。也可以应用这一物理含义直接列写系统的柔度矩阵，称为柔度影响系数法。 同样，我们也可以用影响系数法来列写质量矩阵，方法是：令第 j 个广义坐标的加速度为(4.11)其余广义坐标的加速度为 0 ，为此而需要在各个
6、广义坐标方向上施加的广义力向量就是质量矩阵的第 j 列。1jq 对于直梁，经常用几个位置的挠度作为广义坐标，来近似描述直梁的振动。这时，采用影响系数法，建立梁的柔度矩阵更方便的，因而需要用到简单边界条件下梁的挠度公式。 简支梁在横向集中力作用下的挠度公式为xPabllxaxxblaxblEIlPbfaxbxlEIlPbxf,)()(60),(63223222 例4.1 写出图示梁的柔度矩阵，梁的抗弯刚度为EI。如果将梁的质量按分段区间均分到区间的两个端点，写出梁的质量矩阵，设梁单位长度的质量为 rl。4l4l4l4l1m2m3m2y1y3y例4.1图ijibjbl图(a) 解：根据材料力学，参
7、见图a，在梁上 j 点作用单位横向载荷，将在 i 点产生挠度 fij ：222(),6ijijijijjibb lbbfffEIl 现在，分别在三个挠度坐标（本题的广义坐标） y1、 y2、 y3 方向作用单位力，每次得到一个挠度向量，将它们写成矩阵就是柔度矩阵。比如在y1方向作用单位力，按照前面的挠度公式，可得三个坐标方向的挠度为：222311( /4)(3 /4)( /4)(3 /4) 96768llllllfEIlEI222321( /4)( /2)( /4)( /2) 116768llllllfEIlEI222331( /4)( /4)( /4)( /4) 76768llllllfEI
8、lEI因此，柔度矩阵的第一列为3112131,9,11,7768TTlfffEI类似可算出柔度矩阵的第二、第三列。柔度矩阵为11121332122233132339117 11 16117687119ffflFfffEIfff系统的动能为1122211223312322330011(), 002200myTm ym ym yy yymymy。所以，系统的质量矩阵为1230010000010400001lmlMmmr其中1224llmmmr 例4.2 建立图示三级摆的线性自由振动方程。设1m1232m3m1l2l3l例4.2图 解：我们用Lagrange方程来建立振动方程。 各质点的速度为123
9、123,mmmmllll11111(cossin)vlij211221122(coscos)(sinsin) vlij3112233112233(coscoscos)(sinsinsin) vlij所以系统的动能为32221112222123121223233 1311(cos)(sin) 21324cos()22cos()2cos()kkiiiikiiTmlml 系统的势能为1236(3cos2coscos)Vmglmgl（注意：为了得到线性振动方程，能量表达式必须保留到二阶微量）微振动时， 为小量，将以上能量保留到二阶小量，得,ii222212312233 1121232313242223
10、211 , 2212111Tmlml 2221221123231(32)23001 , 0202001Vmglmgl 代人Lagrange方程0,1,2,3iiidTTVidt1122233321300022102001110010mlmgl 得系统的振动方程为 例4.3 如图所示结构，刚性矩形板由三根长度均为L的无重弹性支柱支撑，支柱与板和地面刚结，每根支柱抵抗端点位移产生的弯曲刚度为12EI/L3，支柱的扭转刚度不计。如图所示，取板的广义坐标为A、B和E三点在水平面内的位移，即广义坐标列向量为 q = v1, v2, v3T。用刚度影响系数法求刚度矩阵。例4.3图2Dv 2Av 11v A
11、BCDE瞬心位移图受力图AkvABCDEDkv1Q2Q3Q1231,0vvv时板的位移和受力图图(a) 解：（1）求刚度矩阵第一列参见图 a，可得板的力平衡方程：31120()202()02DAADADQkvkvQ LkvkvLQQkvkv；其中 312EIkL解得1234 ,2 ,0QkQkQ 因此，刚度矩阵第一列为123, , 4, 2, 0TTQQQk2131,0vvv时板的位移和受力图图(b)位移图1Dv 2Bv 21v 1Av ABCDE瞬心1Ev 受力图AkvEkv1Q2Q3QCBADEDkv （2）求刚度矩阵第二列参见图 b，可得板的力平衡方程：31120()00ADADEQkv
12、kvQ LkvkvLQQkv；其中 312EIkL因此，刚度矩阵第二列为123, , 2, 3, 0TTQQQk解得1232 ,3 ,0QkQkQ 3121,0vvv时板的位移和受力图图(c)受力图Akv1Q2Q3QCBADEDkv3kv位移图1Dv 1Av ABCD31v E板平动 （3）求刚度矩阵第三列参见图 c，可得板的力平衡方程：331120()00ADDAQkvkvkvQ LkvkvLQQ；其中 312EIkL因此，刚度矩阵第三列为123, , 0, 0, 3TTQQQk解得1230,0,3QQQk 综合以上结果，得系统的刚度矩阵为342012230 ,003EIKkkL4.2 无阻
13、尼系统的自由振动与模态1. 自由振动与模态的产生设系统自由振动方程为0M xKx(4.12)其中 x 为 n 维列向量。设方程的解为nRtx),sin(4.13)代入（4.12），得0|0)(22MKMK(4.14)方程（4.14）为方程（4.12）的特征（本征）方程，因此线性系统自由振动的求解转化为相应特征值问题的求解。由方程（4.14）可以求出 n 个特征值22221,n因为M、K 矩阵为实对称正定和半正定矩阵，根据线性代数理论，将保证以上所有特征值大于或等于零。因此，常数2222211,nn对于自由振动解（4.13）是有意义的，它们代表了自由振动的固有频率。由此有结论，n自由度系统有 n
14、 个固有频率。进一步，对应于每个特征值（每个固有频率），由方程（4.14）可求出一个特征向量，有)()2()1 (,n只要所有特征值都是单根，所有特征向量将是线性独立的。这样，我们得到了n个特征对niiii, 2, 1;,)(每一个特征对叫做一个模态（mode），i 称为模态频率（或固有频率）， i 称为模态振型（modal shape）；经常也将振型称为模态。因为以上模态都是实数，因此称为实模态实模态。对于自由振动，n个模态就是n个线性无关的自由振动解（但模态还有更大的意义，下面逐渐介绍）。由此得线性系统自由振动通解为( )100sin(),0(0),(0)iniiiiiiiiiixatat
15、xxxx常数 、由初始条件确定：时：(4.15)2. 模态的代数性质每个模态应满足特征方程（4.14）：0)(, 0)(0)(, 0)()(2)()(2)()()()(2)(2jjTiiiTjTiTjjjiiMKMKMKMK，得、两个方程分别左乘两式相减并考虑到K、M的对称性，得(4.16)2)()()()()()()()()()(22)16. 4(; 0; 0,0)(iiTiiTiiTjiTjjiiTjjiMKjijiKjiMM得时，由当进而得当(4.17)(4.18)(4.17)、(4.18)表示模态的一个最重要的性质，称为振型关于质量矩阵和刚度矩阵的正交性，简称振型的正交性或模态的正交性
16、。模态的其它性质是特征值问题共有的。第二个重要性质是：任一振型乘以一个不为零的常数仍为相应模态的振型。第三个重要性质是：线性系统（4.12）经相似变换不改变系统的特征值。(4.18)(4.19)4.3 振动分析的模态叠加法应用模态的性质，现在推导多自由度系统振动分析的分解、叠加方法。因为在所有特征值为单根的条件下，所有振型是线性无关的，因此我们取坐标变换称为模态矩阵其中nTPnPPPnnnPRxxxxRxx,21)()2()1 (代入无阻尼振动方程 PPTTTTPPPPPPPMxKxQMxKxQM xK xQ 得两边同乘，得应用模态正交性，得(4.20)(4.23)M xKxQ(4.21)(4
17、.22)11diag(,),diag(,),TPPPnTPPPnMMmmKKkk TPnPTPiTiPiiTiPiQQQQKkMm,1)()()()(模态刚度模态质量这样，原方程已分解为（4.23）的n个独立的单自由度振动方程，这种过程称为解耦，解耦后的坐标 xP 称为模态坐标或主坐标。求出主坐标的解xP(t) 后，由坐标变换式（4.20）得到原广义坐标的解，这个变换的本质就是线性叠加。到此，利用模态对线性系统振动进行分解、叠加的完整分析方法已建立，称为模态叠加法。 对方程（4.23）再作变换，可写成标准形式。令1,1,1diag21PnPPPNPPmmmMxMx其中(4.24)PTPPNnP
18、PPNNNPTNPPTNPPTMMQMQMKMQxxQMxMKMxMMMPPPP ,diag)23. 4(22221其中变为方程标准方程（4.26）在多自由度系统理论中称为正则形式或简正形式(normal form)，坐标 xN 称为正则坐标或简正坐标。(4.26)(4.25)即由广义坐标方程（4.21）到简正坐标方程（4.26），可写成一个总的坐标变换，称为正则（简正）变换(1)(2)( )12(1)(2)( ), PPNNNnNPPPPnnNNNxxM xxMmmm 其中 正则模态矩阵(4.27)下面我们回过头去考察模态叠加法的数学和力学本质。将变换过程重写于下：)()()()()(eQx
19、KxMdQxKxMcQxKxMbxxaQKxxMPPPPPTPTPTTPPP 应用模态正交性，得，得两边同乘得取变换原广义坐标方程因为方程(a)的解归根到底是一个 n 维向量，它可用任意一组n维向量基线性表示，并且向量基选定后，这种表示是一一对应的，因为变换(b)是非奇异变换，因此变换(b)使方程（c）与方程（a）等价（两个方程的解相互满足）。只要方程（c）成立，方程（d）必成立，反过来，若方程（d）成立，则因此方程（c）与方程（d），即方程（e）等价。最终方程（a）与方程（e）等价。00)(QxKxMQxKxMPPPPT 以上从数学上证明了模态叠加法的正确性。同时看到，在方程（c）两边同乘任
20、何非奇异矩阵都是可行的，但是同乘 T ，不但可以利用模态的正交性来对方程解耦，而且符合动力学基本原理。因为经（a）式变换后2121111,222111,222nTTTPPPiPiinTTTPPPiPiiTTTPTx MxxMxm xVx KxxKxk xWx QxQ动能势能广义力的虚功因此，模态叠加法不但数学上是正确的，而且完全满足动力学基本原理；同时还将系统的动能、势能变换为各模态能量的简单叠加，这一结果有实际应用价值。Lagrange ( )TTTPPMxKxQd 由方程得变换后系统的动力学方程为此式正是上面式4.4 无阻尼系统的强迫振动响应1. 简谐激励响应本节应用模态叠加法和其它方法，
21、对多自由度无阻尼系统在不同激励下的响应，给出明确的分析计算方法。设系统简谐激励方程为00010,i tTnM xKxF eFFF(4.28)设稳态解为ntiCXXex,代入（4.28）得01202)()(FMKXFXMK所以稳态响应为titieFHeFMKx0012)()(4.29)122)()()(MKHMK称为系统的阻抗矩阵或动刚度矩阵，称为系统的频率响应矩阵或动柔度矩阵。阻抗矩阵和频响矩阵的物理含义是不同激励频率处激励力幅（列阵）与响应幅值（列阵）之间的比例系数矩阵。实际中可以用实验测试建立这两个矩阵。以上方法求响应直观、简单，但对高阶系统，实际应用却有很大的困难，原因是频响矩阵的获得需
22、要求矩阵的逆，高阶矩阵求逆不但比较困难，而且精度很难保证；还有，在设计系统时，一般要求知道频响矩阵在一段频带上的分布，需要逐点求出频响矩阵，计算量太大、可信度难以保证。下面用模态叠加法来求简谐响应。对方程（4.28）作模态坐标变换，得0)(10)()(1)()(1,),(diag),(diagFFFFFFKkkkKKMmmmMMeFxKxMTjPjTPnPTPjTjPjPnPTPjTjPjPnPTPtiPPPPP (4.30)下面用模态叠加法来求简谐响应。对方程（4.28）作模态坐标变换，得解（4.30）得jjtinjjPjTjjtinjjPjTjjnjPjjPtijPjTjtijPjTjti
23、PjPjPjPjseFskeFmxxxeskFemFemkFx其中;)1 ()()1 ()(102)()(1022)()(1)(20)(220)(2(4.31)比较（4.31）与（4.29）可得，频响函数矩阵的模态展开式为njjPjTjjnjjPjTjjmskH122)()(12)()()()1 ()(4.32)如果上式中的振型取为正则（简正）振型，则有njjjTjNjNsH122)()()1 ()(4.33)显然，当sk = 1时，H()的非零元素，系统振幅 ，系统共振，因此 n 自由度系统有n 个共振频率。当简谐激励频率在某个共振频率 k 附近时，系统的响应将主要由（4.31）式右端的第
24、k 项确定，即有( )( )02,(1)kk Ti tkPkkxF eks当时(4.34)任意周期激励的响应，只要将周期激励展开成Fourier级数，再将级数中每项简谐激励的响应叠加即可。例4.4 求图示系统的稳态受迫振动。解：系统的动力学方程为TTFFkkkkkkkKmmmMxxxxtFKxxM, 0 , 0,0202,diag,sin03210其中 可求出模态为例4.4图555. 0247. 11,802. 0445. 01,247. 2802. 11247. 3,555. 1,198. 0)3()2()1(232221振型固有频率mkmkmk振型可表示成下图：11.8022.2470.4
25、4510.80211.2470.555tFtFmmmmmmMmjtFmxxxxTjPjPPPjTjPjPjPjPjjPjPsin)(863. 2,841. 1,296. 9,3 , 2 , 1),(1,0)(321)()(2)3()2()1 ( 得解耦方程应用模态变换tmFmtFxjPjTjjPjPjPjsin)()()(220)(22解得tFHHHtFHtFmxxxnjjPjTjjnjPjjPsin)()()(sin)(sin)(33231301022)()(1)(所以132222222322222233222210.2420.4360.194( )0.198/1.555/3.247/10.
26、4360.1940.242( )0.198/1.555/3.247/10.5430.3500.108( )0.198/1.555/3.247BBBBBBBBHkHkHk代入数据可算得22/B2113132232333333/( )sin( )sin( )sin( )sin( )sin( )sinBk mFxHFtHtkFxHFtHtkFxHFtHtk其中因此各坐标分量的响应为可以画出图形如下数，实际上就是振幅放大系，上式中的)()()(332313HHHH13()H23()H33()21)(B22)(B23)(B2)(B2. 瞬态响应任意激励的振动方程为1( ),( )( ),( )TnM x
27、KxF tF tF tF t(4.35)两边作Fourier变换，得所以221( )( )( )( )( )( )( ) ( )1( )( ) ( )2i tM XK XFXKMFXHFx tHFed(4.36)(4.37)0( )( ) ()( ) ()tx thF tdhF td(4.38)由卷积定理直接得时域表达式，当 时，有 0( )0tF t。产生的响应为阵，其各元素称为系统的脉冲响应矩逆变换。的为其中)()()()()()()(21)(Fourier)()(thtxttFththdeHthHthijijijti根据频响函数矩阵H()的模态展开式（4.32）或（4.33），可将频率响
28、应（4.36）写成(4.39)njjPjTjjmFFHX122)()()()()()()(查Fourier变换表得到22Fourier1sinjjjt变换njtjjPjTjjdtFmtx10)()()(sin)(4.40)(4.40)式可由另一条路径得到。对(4.35)式作模态变换得)()(,)(,),()()(),(diag),(diag)()(1)()(1)()(1tFtFtFtFtFtFKkkkKKMmmmMMtFxKxMxxTjPjTPnPTPjTjPjPnPTPjTjPjPnPTPPPPPPP 得(4.41)对(4.42)中的各个方程应用Duhamel 积分，得tjjPjTjtPjj
29、jPjPjdtFmdtFmtx0)(0)(sin)(sin1)(4.42)由(4.41)就得到(4.40)。4.5 阻尼系统及其求解1.阻尼的结构、系统动力学方程 实际的阻尼机理都是很复杂的，对其建立精确的模型是很难的。当系统的阻尼比较弱时，处理方法还是将其等效为粘性阻尼。因此在多自由度系统中，一般将阻尼产生的广义力假设为系统广义速度的线性函数qCQniqcQnjjijdi即, 2 , 1,1(4.43)将q 改记为x，得系统动力学方程( )M xCxKxF t(4.44)矩阵C一般为对称正定或半正定矩阵。2. 响应的求解方法（1）模态阻尼方法这种方法认为矩阵C 可用实模态矩阵近似解耦，即,d
30、iag1PnPPTccCC因此通过模态变换 x = xP 方程（4.40）变换成n个独立的单自由度方程njtFxkxcxmPjPjPjPjPjPjPj , 2 , 1),(4.45)(4.46)njtBxxxPjPjjPjjjPj , 2 , 1),(22化成标准形式(4.47)求出xP后，由模态变换返回到系统的原坐标，就得到系统的响应。在实际应用中，这种方法对阻尼有两种处理方法：(1) 不去建立阻尼系数矩阵C ，而是凭经验或实验指定各阶模态的阻尼比j。(2) 认为阻尼为比例阻尼，即假设阻尼系数矩阵C为bKaMC(4.48)通过实验确定比例系数a、b。例4.4图F(t)例4.5 对如图阻尼系统
31、，(1) 建立阻尼矩阵，考察是否能表示成比例阻尼。(2) 如果有k1 = k2 = k3 = k， m1 = m2 = m3 = m， c1 = c2 = c3 = c，F(t) = F cos t， 求系统的稳态振动。解：阻尼力产生的虚功为1 1122121332321212212 123233232333122112322332333() ()() ()()()()0,0dWc xxc xxxxc xxxxcc xc xxc xcc xc xxc xc xxcccxxxxccccxccx 所以广义力12212233233300dcccxQccccxCxccx 阵为系统动力学方程为tFtFK
32、xxCxMcos)(0 333322221333322221321032100,00,diag, 0 , 0,cccccccccCkkkkkkkkkKmmmMFFxxxxTT其中bKaMC如果阻尼能表示成比例阻尼，则必须有等价于时，以上方程满足。当33221133223333232221211, 0)()(kckckcbacbkcbkcbkamcckkbamcckkbam因此，本题的阻尼当满足上述条件时，可以表示为比例阻尼。如果有k1 = k2 = k3 = k ， c1 = c2 = c3 = c ，则满足上述条件，系统为比例阻尼，可通过模态变换解耦。取模态变换,4.4Pxx 振型矩阵 已由
33、例求出tFtFmkmkmkmcmcmcmmmmmmjtFmxxxTjPjPPPPjPjPjjPjjPjcos)(247. 3,555. 1,198. 0247. 3,555. 1,198. 0863. 2,841. 1,296. 93 , 2 , 1),(10)(2322212113212 得解耦方程)cos()()(tan)()()(22220)(221)(22220)(220)(jjjPjTjPjjjjtijjPjTjtijjPjTjPjtmFxemFeimFxj对应于余弦激励的解为解为对应于指数简谐激励的3122220)()()cos()()(jjjjPjTjjPtmFxx所以(2) 复
34、模态方法当阻尼系数矩阵C不能解耦时，以上方法不能采用。)(tFKxxCxM (4.49)设M、K、C 均为实对称正定矩阵。求解这一方程的基本思想仍然是采用解耦方法，先来求出齐次方程的特征对0KxxCxM 设 x = el t 对应的特征方程为02llKCM(4.50)由此可解出2n个特征根和特征向量niii2, 2 , 1,l与实模态不同，现在 li 可以是实根或复根。如果li 是实根，则一定是负数，因为我们假定了阻尼矩阵C为正定，这时对应于衰减自由运动。如果li 是复根，它一定具有负实部，且由于方程的系数都是实的，所以复特征根一定是共轭成对地出现，进一步复特征向量也是共轭成对地出现。每一对共
35、轭复根对应于特定频率和衰减率的一种衰减自由振动。系统的复特征值和复特征向量称为系统的复模态。现在的复振型矩阵为nnnC2221,因此它的各列不独立，进而不能直接用它来对方程（4.49）解耦，但可用方程扩阶的方法解决。将方程（4.49）改写为(4.51)(4.52)0)(00,0)(0)(,)(lKMKMKMKCMMMtFtFxxytFyKyM特征值问题为都是对称矩阵。对应的和显然，其中niiiiii2, 2 , 1,ll(4.54)(4.53)容易验证，(4.51)式的复模态经组合将是(4.54)的特征问题解再来看特征向量i 的正交性。由(4.54)有0, 0jjjiiiKMKMll(4.55
36、)0,0)00,0,TTjiTTTTijijijijijTijjiijTjiTjiTTjjPjjjPjPjjPjMKMKMMijKijMmKkkmllll lll 以上两式分别左乘、得两式相减得 （因此，只要，就有正交关系进而而因此我们可以用复特征向量系i , i = 1, 2, , 2n对方程(4.53)解耦，解耦变换为,;,21221nndiagzyll其中，(4.59)(4.56)(4.57)(4.58)njtFmzztFtFtFtFtFkkKKmmMMnjtFzkzmtFzKzMPjPjjjjTjPjTTPnPPTPnPPTPPjjPjjPjPPP2 , 2 , 1),(1)()(),
37、()()(,diag,diag2 , 2 , 1),()()(2121l或写成其中程解耦为由复模态的正交性，方(4.60)方程(4.53)变为)(tFzKzMTTT(4.61)方程(4.61)对应于零初始条件的解，可用积分变换法求解。方程(4.61)两边作Fourier变换，得Fourier1( )( )( ),1( )1( )(/ )1Fourier(/ )jjjjPjPjPjjPjjtji ZZFimFZmiieiillll变换查变换表知道：由卷积定理得ttTjPjttPjPjtPjPjjdeFmdeFmetFmtzjjj0)()(0)()(1)(1)(1)(lll(4.62)(4.63)
38、其中算符 “ ” 表示卷积。根据解耦变换式(4.59)，有zxzxxzzy( )( )22()( )011( )jjj TnnttjjjjPjxzzFedml (4.64)(4.65)到此，我们已将线性阻尼方程(4.49)的任意激励解表达成复模态叠加的卷积积分形式。当系统受到简谐激励时ntiRFeFtF00,)(由(4.65)式，系统的响应为(4.66)njttijPjTjjnjtitjPjTjjnjtitPjTjjnjttiPjTjjjjjjjjeeimFeeimFdeemFdeeFmx210)()(210)(0)()(210)(0)()(210)(0)()()()()(llllllll最后
39、来求系统对应于初始扰动的响应。设系统的初始条件为：当 t = 0时，有0000(0),(0);,nxxxxxxR(4.67)(4.68)0()0(,)0(100yzxxy于是0001001111,0,)0()0(CxxMMxMxxCMMMMyMzMMMMTTPTTPTPTPPT因此有故有，考虑到njCxxMMxmzjPjTjj2 , 2 , 1,)0(000)(l可写成njzzjjj2 , 2 , 1, 0l写成系统的自由振动方程可( )000(0)( )(0),1,2,2jjj TttjjjPjzz tzeMxMxCx emjnlll对应于初始条件的解为(4.69)2( )1( )( )20
40、001(4.64)( )jnjjjjj TntjjPjxxzz tMxMxCx emll根据式返回到原坐标 ，有(4.70)4.6 模态问题的一些特殊情况 1. 等固有频率（重特征值）的情形当特征值均为单根时，n 维系统有 n 个线性无关特征向量，它们构成 n 维线性系统的完备解耦基。当特征方程出现重根时，对应的特征向量不能唯一确定。比如，设特征方程有二重根 l1 = l2 = lr ，对应的特征向量为 (1) 与 (2) ，所以的特征向量。对应于的任意线性组合仍然是、即（因此rrrrbaMbaKMKMKllll)2() 1 ()2() 1 ()2() 1 ()2()2() 1 () 1 ()
41、(),现在的问题是：当 r 重根对应于 r 个线性无关的特征向量时，如何选定这 r 个线性无关的特征向量？回答是应按正交性条件来选取，即所选定的特征向量必须满足M、K的正交条件。举例说明如下：如图4.1系统，质量和刚度矩阵为,401001010100111113,1, 1, 1, 1 diag4321mkmkkKmMllll，容易求得特征值为：容易确定对应于 l1 、 l4 的特征向量为：TT 1, 1, 1, 3, 1, 1, 1, 1 )4()1 (m图4.1231123423141234234340211120100001000010000;01,00,1kmkKMmll对应于的特征向量
42、由下式求解：即得、 、 中有两个数可自由选取。令 、 为自由变量，它们的两组独立的取值对应于二重根的两个独立的特征向量，不妨取和，因此有TTTTTTTTTaaMaaKMmkKMKM2, 1, 1, 0 1,21,21, 021012,01,00 1, 0, 1, 0,0, 1, 1, 0)3()3()3()2()3()2()3()2()3()2()2()3()2()3()2()3()2()3()2()3()2()3()2()3()2()3()2(，规范化得得由，设程如下：的一个即可。正交化过个方程中，所以只需满足上述两因为，使得、线性组合的两组、正交化，这就是要求出、进行、再将48, 6, 2
43、, 0diag12, 6, 2, 4diag1201111111113001kKmMTT正交性检验：型矩阵为现在可以写出系统的振可见结果完全正确。2. 固有频率随系统参数的变化研究固有频率随参数的变化规律，目的是寻找系统的动态设计和改进振动系统的方法。设系统的特征值问题为MKDDDMKjjjj2)(2, 00)(假定质量、刚度矩阵随参数 s 变化，上式对 s 求偏导数( )( )( )( )( )( )( )( )( )( )( )( )2( )00(4.71)00(2)0jjjjjjj Tj Tjj Tjjj Tjj Tjjjjj TjjjDDssDDssDDDsKMMsss左乘得由式(4.
44、71)(2)()()()(211jjTjjjjTjsMsKsM假定振型已经正则化，即(4.72)这就是固有频率随参数的变化规律，它实际上是固有频率的参数灵敏度。如果系统有多个参数可以调整，通过计算和比较灵敏度的大小，可知道调整哪些参数更有效。例4.6 图示为发动机的扭振模型。各个 J 为转动惯量，各个k 为轴的扭转刚度；具体值见下表。考察二阶固有频率对各个J和k的灵敏度。例4.6图下标 jJj / kgm2kj / 105 Nm/rad10.101.0020.100.5030.500.2540.301.0051.10解：系统的质量、刚度矩阵为44443333222211115432100000
45、0000000,diagkkkkkkkkkkkkkkkkKJJJJJM求得系统的固有频率和振型矩阵为0002. 03620. 01128. 05377. 06901. 00052. 06438. 12777. 02792. 06901. 01164. 02698. 07905. 09014. 06901. 04946. 21324. 04416. 10977. 16901. 09259. 12668. 01038. 21479. 16901. 0s/rad102953. 2,5037. 0,3148. 0,0437. 0, 062521：根据灵敏度公式(4.72)可得，5 , 1,4 , 1,
46、21)2()2(22)2()2(22imMmikKkiTiiTi下标 irad/s/kgm2rad/s/Nm/rad12.751026.02 10-622.52 1029.22 10-531.70 1023.33 10-341.63 101.06 10-456.05 10iJ2ik2灵敏度的计算结果如下表：3. 约束对固有频率的影响多自由度系统加上约束后，固有频率会发生相应的变化，现在来定性考察这一问题。设系统已经解耦成正则坐标形式：现在假设对正则坐标集加上一个线性约束阶单位矩阵为，1,01111111nEzzycccycEzzczzbTnTnTnjjjnniii,diag, 01, 1, 0
47、2212niiizznizz 即(4.73)(4.74)原方程(4.73)变成0,ycEcEycEcETT ,diag, 021212nTnTyccyyccy 对应的特征值方程为222222()()0( )det()()0TnTnEccEcc频率方程为这是一个关于 2 的n 1次多项式求根问题，它有n 1个根，为了判断这n 1个根的值，我们来考察(i)，i = 1, , n。222111211 21122221222 12211222111 11 21()det()()0()TnnnnnnnnEcccccccc ccc cccccc)()()()()()()()()(:)()(000000)(
48、0)(2212222212221221222121211222222322212222122123212221112121212212121221nnnnnnnnnnnnnnnnnccccccccc而同理可得上式第一行元素乘以(ci / c1)再加到第i行上去，可得假设n 个固有频率已经安递增次序排列n21那么以上各个(i)的正负号刚好是交替出现的。这表明 i 与 i+1 之间有() 0 的一个根，因此约束后系统的n 1个特征根（固有频率）镶嵌在原系统的n 个固有频率之间。如图6.3所示。(1)(2).()123n-1n约束系统的 n 1个固有频率图4.2(n-1)(n)4.7 传递矩阵法 1
49、. 轴的扭转振动传递矩阵法是链式结构动力学建模的一种离散集结方法，可以用来计算模态振动。多盘扭振系统的建模方法如下：将图5.1的整体系统分割成结构形式相同的两端单元，如图5.2。将盘视为刚性质量元件、轴段视为柔性无质量元件。无质量轴段的扭转运动取决于轴段两端的转角和扭矩，因此取各元件的转角和扭矩作为状态变量TTX,2101RLiiRLiiiii tiRPLiiiPiiTTJeXS XSJ 设盘作频率为 的模态振动，即，点传则得称为递矩阵。于是，第 i 个盘的动力学方程为（参见图5.3）对于第 i 个轴段（参见图5.3） ，平衡与变形方程为111(),1101RRiiLLRiPiiiiiiiLF
50、RiiiFiiTTG ITkklXS XkS写成状态变量的矩场阵形式为称为传递矩阵。于是得到第 i 单元两端状态的传递关系（从第 i 1 个盘右侧到第 i 个盘右侧）RiiRiFiPiRiXSXSSX1122111iPFiiiiiikSS SJJk单元传其中为递矩阵。对于有n个盘的轴系，将各单元的传递关系连乘就得到整个系统两端状态的传递关系：RRRnRnRRnnRnTSSSSTSXXSSSX00222112110011)()()()(即再根据轴系两端的边界条件，比如两端自由的情况，有就是系统的固有频率。满足以上第二个方程的，有0)()(00210110SSTTnRRn2. 梁的横向振动将梁结构
