第六节-两个自由度体系的自由振动PPT课件

1、2021/7/24111-6 11-6 两个自由度体系的自由振动两个自由度体系的自由振动 在实际工程中，很多问题都是简化为多自由度体系来计算。多在实际工程中，很多问题都是简化为多自由度体系来计算。多自由度体系在强迫振动时的动力反应与体系的动力特性自由度体系在强迫振动时的动力反应与体系的动力特性( (自振自振频率、主振型等频率、主振型等) )有密切关系。有密切关系。分析多自由度体系自由振动的分析多自由度体系自由振动的目的是确定体系的自振频率和相应的主振型。目的是确定体系的自振频率和相应的主振型。 在单自由度体系自由振动的分析中已经看到，阻尼对自振频率在单自由度体系自由振动的分析中已经看到，阻尼对

2、自振频率的影响很小，在多自由度体系中也是如此。另外，在分析多自的影响很小，在多自由度体系中也是如此。另外，在分析多自由度体系强迫振动的动力反应时，常要用到不考虑阻尼情况下由度体系强迫振动的动力反应时，常要用到不考虑阻尼情况下体系的主振型。所以在分析多自由度体系的自由振动时，体系的主振型。所以在分析多自由度体系的自由振动时，不考不考虑阻尼的影响。虑阻尼的影响。 本节先讨论两个自由度体系，下节再推广到本节先讨论两个自由度体系，下节再推广到n个自由度体系。个自由度体系。与单自由度体系一样，两个自由度体系建立运动方程也有柔度与单自由度体系一样，两个自由度体系建立运动方程也有柔度法和刚度法，分别讨论如下

3、。法和刚度法，分别讨论如下。 2021/7/242一、柔度法一、柔度法 1.运动方程的建立运动方程的建立2021/7/243按动静法，将惯性力按动静法，将惯性力 和和 分别作用在质点分别作用在质点m1 1和和m2 2上上( (图图11-11-34b)34b)，则质点位移，则质点位移y1 1( (t) )、y2 2( (t) )应等于这两个惯性力共同作用所产生的静应等于这两个惯性力共同作用所产生的静力位移。根据叠加原理可得力位移。根据叠加原理可得 )(11tym )(22tym 222221112122211111)()()()()()(tymtymtytymtymty 或或0)()()(0)(

4、)()(222221121122121111tytymtymtytymtym 式中式中 ，的物理意义如图，的物理意义如图11-34c11-34c、d d所示，它们是结构所示，它们是结构的柔度系数。根据位移互等定理，的柔度系数。根据位移互等定理， .22122111、21122021/7/2442.频率和振型的计算频率和振型的计算0)()()(0)()()(222221121122121111tytymtymtytymtym 注意到单自由度体系的自由振动为简谐振动，假定微分方程组的特解为两注意到单自由度体系的自由振动为简谐振动，假定微分方程组的特解为两个质点作同频率、同相位的简谐振动，个质点作同

5、频率、同相位的简谐振动，)sin()()sin()(2211tAtytAty)sin()()sin()(222211tAtytAty 即即 二阶导数为二阶导数为将上面式子代入运动方程，消去公因子将上面式子代入运动方程，消去公因子sin(t+ )，经整理后得，经整理后得 0)1(0)1(222221112221212111AmAmAmAm2021/7/2450)1(0)1(222221112221212111AmAmAmAm右式是以质点振幅右式是以质点振幅A1和和A2为未知量的齐次线性为未知量的齐次线性代数方程组。其中零解对应于无振动的情况，代数方程组。其中零解对应于无振动的情况，不是所要求的解

6、答。为使方程组具有非零解，不是所要求的解答。为使方程组具有非零解，则其系数行列式必须等于零，即则其系数行列式必须等于零，即 01122221212122111mmmmD上式称为频率方程，用它可求出体系的自振频率上式称为频率方程，用它可求出体系的自振频率。 令令 = ，并将上式展开得，并将上式展开得 210)()(2122211212221112mmmm由此可解出由此可解出 的两个正实根的两个正实根 (大值大值)和和 (小值小值)如下：如下： 214)()(212122122221112221112, 1mmmmmm2021/7/2464)()(212122122221112221112, 1m

7、mmmmm111221于是求得频率的两个值为于是求得频率的两个值为 两个自由度体系有两个自振频率，其中较小的一个用两个自由度体系有两个自振频率，其中较小的一个用1表示，称为第一频率表示，称为第一频率或基本频率；另一个用或基本频率；另一个用2表示，称为第二频率。相应的两个自振周期分别为表示，称为第二频率。相应的两个自振周期分别为: : 221122TT0)1(0)1(222221112221212111AmAmAmAm将将1、2分别代入下式可求相应的分别代入下式可求相应的A A1 1和和A A2 22021/7/247当当=1或或 =2使方程组使方程组(11-46)的系数行列式等于零，因此它的两

8、个方程的系数行列式等于零，因此它的两个方程不是独立的，只能由其中的任一方程求出不是独立的，只能由其中的任一方程求出A1与与A2的比值。的比值。当当=1时时, ,此时此时A1用用 表示，表示，A2 用用 表示，则由式表示，则由式(11-46)的第一式得的第一式得 )1(1A)1(2A12121112)1(1)1(21mmAA相应地，得到质点位移相应地，得到质点位移y1(t)、y2(t)的一个特解的一个特解 )sin()()sin()(11)1(2211)1(11tAtytAty由此可知由此可知 。它表明：在自由振动过程中，两质点位。它表明：在自由振动过程中，两质点位移的比值保持为常数移的比值保持

9、为常数 ，也就是说在，也就是说在任何时刻体系的振动都保持同一形任何时刻体系的振动都保持同一形状。这种相对位移保持不变的振动形式称为主振型，简称振型。状。这种相对位移保持不变的振动形式称为主振型，简称振型。 1)1 (1)1 (212)()(AAtyty12021/7/248同理，对于同理，对于 的情况，有的情况，有 2)sin()()sin()(22)2(2222)2(11tAtytAty221211122)2(1)2(2121)()(mmAAtyty当体系按当体系按 振动时，质点位移振动时，质点位移y1(t)与与y2(t)之比为之比为1: ，称为，称为第一振型或基本第一振型或基本振型振型。当

10、体系按。当体系按 振动时，质点位移振动时，质点位移y1(t)与与y2(t)之比为之比为1: ，称为，称为第二振第二振型型。当多自由度体系按某个主振型作自由振动时，由于振动形式不变，只需。当多自由度体系按某个主振型作自由振动时，由于振动形式不变，只需一个几何坐标即能确定全部质点的位置，因此它一个几何坐标即能确定全部质点的位置，因此它实际上如同一个单自由度体实际上如同一个单自由度体系那样在振动。系那样在振动。 11222021/7/249要使体系按其某一主振型作简谐自由振动，只有在特定的初要使体系按其某一主振型作简谐自由振动，只有在特定的初始条件下才能出现。例如，对应于第一振型，应有始条件下才能出

11、现。例如，对应于第一振型，应有 111)1(111)1(21211)1(11)1(212coscos)0()0(sinsin)0()0(AAvvAAyy这表明只有当质点这表明只有当质点2的初位移和初速度均分别为质点的初位移和初速度均分别为质点1的初位移和初速度的的初位移和初速度的 倍倍时，体系才会按第一振型作自由振动。这种在特定初始条件下出现的运动形式在时，体系才会按第一振型作自由振动。这种在特定初始条件下出现的运动形式在数学上称为微分方程组的特解。数学上称为微分方程组的特解。 13.运动方程的通解运动方程的通解将两个特解进行线性组合就得到通解：将两个特解进行线性组合就得到通解： 2)2(1)

12、2(21)1(1)1(222)2(211)1(2222)2(111)1(11,)sin()sin()()sin()sin()(AAAAtAtAtytAtAty2021/7/24102)2(1)2(21)1(1)1(222)2(211)1(2222)2(111)1(11,)sin()sin()()sin()sin()(AAAAtAtAtytAtAty式中有四个独立的待定常数，它们可由四个初始条件来确定。所式中有四个独立的待定常数，它们可由四个初始条件来确定。所以，给定任意四个初始条件后，即可完全确定体系的自由振动。以，给定任意四个初始条件后，即可完全确定体系的自由振动。由上式可知，在一般初始条件

13、下，质点的位移是由具有不同频率由上式可知，在一般初始条件下，质点的位移是由具有不同频率的简谐振动叠加而成的，它的简谐振动叠加而成的，它不再是简谐振动不再是简谐振动。不同质点的位移的。不同质点的位移的比值也不再是常数，而是随时间变化。比值也不再是常数，而是随时间变化。 需要指出，需要指出，体系能否按某个主振型作自由振动由初始条件决定体系能否按某个主振型作自由振动由初始条件决定，但由式但由式(11-47)(11-50)可以看出，可以看出，体系的自振频率和主振型则完体系的自振频率和主振型则完全取决于体系的质量和柔度系数，而与初始条件无关。全取决于体系的质量和柔度系数，而与初始条件无关。 2021/7

14、/2411例例11-9 试求图试求图a所示等截面简支梁的自振频率和主振型。所示等截面简支梁的自振频率和主振型。 2021/7/2412解：解：(1)求柔度系数求柔度系数体系有两个自由度。作体系有两个自由度。作 图，如图图，如图b、c所示。由图乘法求得柔度系数所示。由图乘法求得柔度系数21MM 、EIlEIl4867,243432112322112021/7/2413(2)求自振频率求自振频率将柔度系数及将柔度系数及m1=m2=m代入式代入式(11-48)求得求得 EImlmEImlm486)(,48615)(312112312111于是得到两个自振频率于是得到两个自振频率 3322331104

15、5.224861692. 5154861mlEImlEImlEImlEI2021/7/2414(3)求主振型求主振型由式由式(11-49)求得第一振型为求得第一振型为 11)(121112112121111)1(1)1(2mmmmmAA由式由式(11-50)求得第二振型为求得第二振型为 11)(121112112121112)2(1)2(2mmmmmAA这表明体系按第一频率振动时，两质点保持同向且相等的位移，其振型是对这表明体系按第一频率振动时，两质点保持同向且相等的位移，其振型是对称的称的 ;按第二频率振动时，两质点的位移是等值而反向的，振型为反对称形按第二频率振动时，两质点的位移是等值而反

16、向的，振型为反对称形状状 2021/7/2415由此例可以看出，由此例可以看出，若结构本身及质量分布都是对称的，则其主振型不是对称的若结构本身及质量分布都是对称的，则其主振型不是对称的便是反对称的便是反对称的。因此，在求自振频率和振型时，可以分别就对称和反对称两种。因此，在求自振频率和振型时，可以分别就对称和反对称两种情况取情况取半结构半结构来进行计算。例如图来进行计算。例如图11-36a所示的对称刚架，质量分布也是对称所示的对称刚架，质量分布也是对称的。其一个主振型必是对称的，另一个主振型必是反对称的的。其一个主振型必是对称的，另一个主振型必是反对称的(图图11-36b、c)。分。分别取半结

17、构别取半结构(图图11-36d、e)，从而变成两个单自由度体系来求自振频率。，从而变成两个单自由度体系来求自振频率。 2021/7/2416例例11-10 试求图试求图11-37a所示刚架的自振频率和主振型。所示刚架的自振频率和主振型。EI=常数。常数。 解：这是两个自由度体系。设质点解：这是两个自由度体系。设质点的水平位移为的水平位移为y1(t)，向右为正；质，向右为正；质点的竖向位移为点的竖向位移为y2(t)，向下为正。，向下为正。则沿则沿y1(t)方向的质量方向的质量m1=2m，沿，沿y2 2(t)方向的质量方向的质量m2=m (1)求柔度系数求柔度系数 作作 图，如图图，如图11-37

18、b、c所示。由图乘法，可得所示。由图乘法，可得 21MM 、EIlEIlEIl2,34,3321123223112021/7/2417(2)(2)求自振频率求自振频率 将柔度系数及将柔度系数及m m1 1=2=2m m、m m2 2= =m m代入式代入式(11-48)(11-48)求得求得 EImlEImlEImlEImlEImlEIml333233312183. 0622,7817. 1622两个自振频率为两个自振频率为 3223111403. 21,7492. 01mlEImlEI(3)求主振型求主振型由式由式(11-49)求得第一振型为求得第一振型为 123.222 37817.133

19、32121111)1(1)1(2mEIlmEIlEImlmmAA2021/7/2418由式由式(11-50)求得第二振型为求得第二振型为 190. 022 32183. 03332121112)2(1)2(2mEIlmEIlEImlmmAA两个主振型的形两个主振型的形状如图状如图11-37d、e所示所示 2021/7/2419二二. . 刚度法刚度法1.1.运动方程的建立运动方程的建立刚度法是列动力平衡方程来建立运动刚度法是列动力平衡方程来建立运动方程，列动力平衡方程时，可以取质方程，列动力平衡方程时，可以取质点为隔离体用平衡条件建立运动方程，点为隔离体用平衡条件建立运动方程，也可以不将质点分

20、离，而按第八章位也可以不将质点分离，而按第八章位移法的步骤来处理。即先在移法的步骤来处理。即先在m m1 1、m m2 2处处沿位移方向加入附加链杆阻止质点的沿位移方向加入附加链杆阻止质点的位移位移( (图图11-38b)11-38b)，然后施加质点的惯，然后施加质点的惯性力性力 ，这时两链杆，这时两链杆的反力分别为的反力分别为 ；其次，；其次，令两链杆发生与两质点实际位置相同令两链杆发生与两质点实际位置相同的位移的位移y y1 1( (t t) )、y y2 2( (t t)()(图图11-38c)11-38c)，此，此时两链杆上所需施加的力为时两链杆上所需施加的力为S S1 1( (t t

21、) )、S S2 2( (t t) )： )()(2211tymtym 及)()(2211tymtym 、)()()()()()(22212122121111tyktyktStyktyktS式中式中k11、k21、k12、k22的物理意的物理意义如图义如图11-38d、e所示，它们是所示，它们是结结构的刚度系数构的刚度系数。 2021/7/2420图图11-38b、c两种情况的叠加，两种情况的叠加，使体系恢复到图使体系恢复到图11-38a的实际运的实际运动位置而处于瞬时平衡，附加链动位置而处于瞬时平衡，附加链杆上的总反力应等于零杆上的总反力应等于零(本来没本来没有附加链杆有附加链杆)。即。即

22、0)()(0)()(22221111tStymRtStymR 将将 代入上式得两个自由代入上式得两个自由度体系自由振动的运动方程为度体系自由振动的运动方程为0)()()(0)()()(2221212221211111tyktyktymtyktyktym )(),(21ttss2021/7/2421仍设其特解为以下形式：仍设其特解为以下形式： )sin()()sin()(2211tAtytAty将上式代入运动方程式，消去公因子将上式代入运动方程式，消去公因子sin( t+ )后，经整理得后，经整理得 0)(0)(2222212121212111AmkAkAkAmk上式是以质点振幅上式是以质点振幅

23、A1、A2为未知量的齐次线性代数方程组，称为振型方程。有为未知量的齐次线性代数方程组，称为振型方程。有非零解的充要条件是方程组的系数行列式等于零，即非零解的充要条件是方程组的系数行列式等于零，即 0222221122111mkkkmkD上式可用来确定频率值，称为频率方程。将其展开得上式可用来确定频率值，称为频率方程。将其展开得 2021/7/24220)()()(212221121222112221kkkmkmkmm由此解得由此解得 的两个根为的两个根为 22121222221112221112214)()(21mmkmkmkmkmk，由式由式(11-53)的第一式可求得两个主振型为的第一式可

24、求得两个主振型为 21211221)2(1)2(211211211)1(1)1(2kkmAAkkmAA2021/7/2423例例11-11 试求图试求图a所示体所示体系的自振频率和主振型。系的自振频率和主振型。已知梁的已知梁的EI=常数，质点常数，质点m1=m，m2=2m，弹，弹簧的刚度系数簧的刚度系数 。 33lEIk 解：设质点解：设质点m1、m2的竖向的竖向位移分别为位移分别为y1 1(t)、y2 2(t),向下向下为正为正 (1)求刚度系数求刚度系数 分别令质点分别令质点m1、m2沿振动正沿振动正方向发生单位位移，作方向发生单位位移，作 图图 如图如图b、c所示。取质点为隔离所示。取质

25、点为隔离体，利用平衡条件求得刚度系体，利用平衡条件求得刚度系数：数： 21MM 、2021/7/2424kklEIk323311kklEIk23322 k12=k21=- -k (2)求自振频率求自振频率 将刚度系数及将刚度系数及m1=m、m2=2m代入式代入式(11-55)，求得，求得 mkmkmk)262(262221，两个自振频率两个自振频率 32311103. 37958. 15250. 18805. 0mlEImkmlEImk，(3)求主振型求主振型 由式由式(11-56)求得两个主振型为求得两个主振型为 12247. 0,12248. 2)2(1)2(2)1(1)1(2AAAA20

26、21/7/2425作振型图如图作振型图如图11-39d、e所示所示 :2021/7/242611-7 11-7 一般多自由度体系的自由振动一般多自由度体系的自由振动 本节讨论本节讨论n个自由度体系的自由振动，并采用矩阵表示形式。个自由度体系的自由振动，并采用矩阵表示形式。一、柔度法一、柔度法 1.运动方程的建立运动方程的建立的图的图11-40a所示为所示为n个自由度体系，个自由度体系，在自由振动的任一时刻在自由振动的任一时刻t，质点，质点mi的位的位移为移为yi(t)，此位移可看作是由各质点，此位移可看作是由各质点的惯性力的惯性力 共共同作用产生的静力位移。采用图同作用产生的静力位移。采用图1

27、1-40b所示的柔度系数，仿照两个自由所示的柔度系数，仿照两个自由度体系的情形，可写出度体系的情形，可写出n个自由度体个自由度体系自由振动系自由振动n个位移方程个位移方程 ), 2 , 1)()(nitymtIiii nnnnnnnnnnnnntymtymtymtytymtymtymtytymtymtymty)()()()()()()()()()()()(22211122222211121122211111 2021/7/2427nnnnnnnnnnnnntymtymtymtytymtymtymtytymtymtymty)()()()()()()()()()()()(2221112222221

28、1121122211111 上式用矩阵形式表示如下上式用矩阵形式表示如下 000)()()()()()(212121 tytytytytytymmmnnnnnnnnn212222111211上式可简写为上式可简写为 0yyM 式中，式中， 和和M分别是柔度矩阵和质量矩阵：分别是柔度矩阵和质量矩阵： 2021/7/2428nnnnnn212222111211nmmmM21由位移互等定理，由位移互等定理， ，故，故 是对称矩阵。对于质点体系，是对称矩阵。对于质点体系，M是是对角矩阵。对角矩阵。 jiij2.频率和振型的计算频率和振型的计算设其特解为设其特解为 yi(t)=Aisin( t+ ) (i=1,2,n) 即各质点的振动频率和相位相同。上式写成矩阵形式即各质点的振动频率和相位相同。上式写成矩阵形式 )sin(tAy式中式中 为振幅列阵，即为振幅列阵，即 AnAAAA212021/7/2429将特解代入运动方程将特解代入运动方程,并令并令 得得21 0)AIM (上式称为上式称为振型方程振型方程，它是位移幅值的齐次线性代数方程组，欲使具有非零解，它是位移幅值的齐次线性代数方程组，欲使具有非零解，则方程组的系数行列式必须等于零，即则方程组的系数行列式必须等于零，即 DIM=0上上式式就是就是n个自由度体系的频率方程