多自由度体系的振动PPT精选文档

1、1第三章第三章 多多 自自 由由 度度体系的振动体系的振动2第三章 多自由度体系的振动 3.1 两个自由度体系的自由振动 3.2 多自由度体系的自由振动 3.3 主振型的正交性和正则坐标 3.4 多自由度体系的强迫振动3 很多结构的振动问题不能按单自由度体系计算，如多层房屋的侧向振动，不等高排架的振动，柔性较大的高耸的结构在地震作用下的振动等，都应按多自由度体系（multi-degree of freedom system）计算。3.1 两个自由度体系的自由振动 对具有无限个自由度的弹性结构，精确地处理其振动问题：有时是非常困难的，在某些情况下也并不必要。在某些特定条件下可对问题作一些简化假定

2、，使一个无限自由度体系离散为有限多个自由度体系，使原来的问题变得容易求解，能获得原结构体系的主要属性和特征。 针对两个自由体系，介绍三种常用的方法： 平衡力系法、刚度法、柔度法4v1、平衡力系法 如图，两集中质量 和 通过三个弹簧 、 和 相互联结，在任意一时刻它们偏离其平衡位置的水平位移分别为 和 。 1m2m1k2k3k)(1ty)(2ty3.1 两个自由度体系的自由振动23122221221111)()(ykyykymyykykym 根据两质量块的平衡条件，可以得到：5表示成矩阵形式：0)(0)(23212222212111ykkykymykykkym 0ykym 式中：1212 , ,

3、 ,ttyy yyy y 整理：32222121 ,00kkkkkkkmmm3.1 两个自由度体系的自由振动2个自由度体系的自由振动写成一般形式：0021222112112122211211yykkkkyymmmm 60021222112112122211211yykkkkyymmmm 对于图中结构体系，有22112322221112112222111,0,kkkkkkkkkmmmmmm3.1 两个自由度体系的自由振动7 假设两个质点为简谐振动，上式的解设为：)sin()()sin()(2211tytytyty 位移振幅 和 ，以及频率 和相位角 均为待定参数。1y2y002122211211

4、2122211211yykkkkyymmmm 3.1 两个自由度体系的自由振动8 1）、在振动过程中，两个质点具有相同的频率 和相同的相位角 。常数2121)()(yytyty)sin()()sin()(2211tytytyty2）、在振动过程中，两个质点的位移在数值上随时间而变化，但两者的比值始终保持不变：3.1 两个自由度体系的自由振动900222121222212111121ykykmyykykmy)sin(t0021222112112122211211yykkkkyymmmm )sin()()sin()(2211tytytyty21111122221 12222()0()0km yk

5、yk ykm y3.1 两个自由度体系的自由振动齐次方程有非零解的条件为其系数行列式等于零，即:0222221121211mkkkmkd 该式是固有频率应满足的条件，称为频率方程或特征方程。（eigen equation or characteristic equation）利用这个方程可计算固有频率102121122211222211122211122121mmkkkkmkmkmkmk展开上式，求得 的两个根为: 2正实根，仅依赖于结构体系的物理性质，即质量和弹簧刚度。0222221121211mkkkmkd3.1 两个自由度体系的自由振动具有两个自由度的体系共有两个自振频率， 12 表示其

6、中最小的圆频率，称为第一圆频率或基本圆频率（fundamental frequency）； 称为第二圆频率。11 比值所确定的振动形式就是与第一圆频率 相应的振型，称为第一振型或基本振型（fundamental mode）21/yy11211111220221 12222()0()0dkm yk yk ykm y 分析频率各自对应的振型(1)112(1)221111ykykm 3.1 两个自由度体系的自由振动1221112)2(2)2(1mkkyy 和 表示第二振型中质点1和2的振幅。 )2(1y(2)2y1m2m下标与质量 和 相对应，上标表示模态号码。)1(2)1(1/yy)2(2)2(1

7、/yy由于模态方程是齐次的，所以 及 只有相对关系。12振型计算公式频率计算公式频率方程)sin()()sin()(2211tytytyty002221212221211111ykykymykykym.振型方程0)(0)(2222212121211211ymkykykymk为了得到y1、y2的非零解，应使系数行列式=00222221121211mkkkmkdf展开是2的二次方程，解得2 两个根为：2121122211222211122211122, 12121mmkkkkmkmkmkmk可以证明这两个根都是正根。与2相应的第二振型：12211122212mkkyy因为d=0，两个振型方程式线性

8、相关的，不能求出振幅的值， 只能求出其比值 求与1相应的第一振型：12111122111mkkyy 132121122211222211122211122121mmkkkkmkmkmkmk212112222211122211122121mmkkmkmkmkmk2 的两个根均为实根；212221121121yykkkkrr2122211211212121yykkkkyyrryy02211ryryyykijjiij矩阵k为正定矩阵的充分必要条件是：它的行列式的顺序主子式全部大于零。故矩阵k为正定矩阵。k11k22-k12k2102 的两个根均为正根；14与2相应的第二振型：212211122212

9、mkkyyf求与1相应的第一振型：112111122111mkkyy2个自由度体系能够按某个主振型自由振动的条件是：初始位移和初始速度应当与此主振型相对应。在这种特定的初始条件下出现的振动，在数学上称为微分方程组的特解，其线性组合即一般解。3.1 两个自由度体系的自由振动15方程的全解:0)(0)(23212222212111ykkykymykykkym )sin()sin()()sin()sin()(22)2(2211)1(21222)2(1211)1(111tyatyatytyatyaty1a122a其中， 、 、 和 由初始条件确定。3.1 两个自由度体系的自由振动一般情况下，体系的自由

10、振动不是主振动，而是两种不同频率及其振型的组合振动:16例 试分析图示结构体系的固有频率和振型。已知： 。kkkkmmm32121,0202212211kykyymkykyym 解：体系的运动方程为：0)(0)(23212222212111ykkykymykykkym 3.1 两个自由度体系的自由振动170202212211kykyymkykyym )sin()()sin()(2211tytytyty体系的运动方程为：设方程的解为：0222122yymkkkmk3.1 两个自由度体系的自由振动上式有非零解的条件是系数行列式为零：02222mkkkmk2/122/113,mkmk展开行列式，可以

11、求得18 第一模态（振型）为两个质量一起振动，无相对位移，中间一个弹簧不起作用，只有第一和第三个弹簧起作用，其结果类似于质量为2m、弹簧系数为2k的单自由度体系的振动；以横坐标表示系统的静平衡位置，纵坐标表示各点的振幅，体系的主振型图。第二模态（振型）为两个质量作相反振动，中间一个弹簧的中点始终不动。计算振型：11221)1(2)1(1)1(mkkyyy11222)2(2)2(1)2(mkkyyy3.1 两个自由度体系的自由振动19 1、主振动:结构体系以某一阶固有频率按其相应的主振型作振动，称为体系的主振动。2、各点同时经过静平衡位置，并同时到达最大偏移位置，以确定的频率和振型作简谐振动。3

12、.1 两个自由度体系的自由振动3、自由度体系自振频率的个数= 其自由度数，自振频率由特征方程求出。4、每个自振频率相应一个主振型。主振型是多自由度体系能够按单自由度体系振动时所具有的特定形式。5、自振频率和主振型是体系本身的固有特性。总结：20y1(t)r2r1y2(t)y1(t)y2(t)r2r122ym.002221212221211111ykykymykykym.0, 0222111rymrym.11ym.2、刚度法3.1 两个自由度体系的自由振动22212122121111ykykrykykr如图，具有两个集中质量的结构体系，两个自由度21y1(t)y2(t)r2r1乘 y1(t)k1

13、1k21乘 y2(t)k12k2211fr1=k11y1+k12y2 r2=k21y1+k22y2kij表示使j点产生单位位移(其它点位移=0)时,在i点需施加的力(称为刚度系数).3.1 两个自由度体系的自由振动22例：m2m1k2k1质量集中在楼层上m1、m2 ，层间侧移刚度为k1、k2k21k111解：求刚度系数： k11=k1+k2 , k21=k2 ,k22k121k22=k2 , k12=k20222221121211mkkkmkd0)(222221221kmkmkk1)当m1=m2=m,k1=k2=kmkmk61803. 225322mkmk38197. 025321()()km

14、kmk02222mkmk61803. 161803. 021 代入频率方程：2121122211222211122211122, 12121mmkkkkmkmkmkmk+231)当m1=m2=m,k11=2k，k12=mkmk61803. 225322mkmk38197. 025321求振型：618. 1138197. 02kkk12k12111mk2111yy1第一主振型：y21=1.618y11=1第一主振型618. 0161803. 22kkk12k12211mk2212yy2第二主振型：y22=0.618y11=1第二主振型240)(222221221kmkmkk 2)当m1=nm2

15、, k1=nk2k11=（1+n）k2，k12=k20)() 1(22222222kmknmkn求频率：求振型：如n=90时1101121yy191222yy当上部质量和刚度很小时，顶部位移很大。（鞭梢效应）2192112yy第一振型：第二振型：特征方程：222214122112mknnn+4121)4121() 1() 1(222212121112nnnnknmknkmkyy+2121122211222211122211122, 12121mmkkkkmkmkmkmk+25y1(t)y2(t)建立振动微分方程：（建立位移协调方程） m1、m2的位移y1(t)、 y2(t)应等于体系在当时惯性

16、力作用下所产生的静力位移。11ym.22ym.1212( ),( )mtm y ty1121p1=11222p2=1 3、柔度法3.1 两个自由度体系的自由振动26222221112122211111)()()()()()(tymtymtytymtymty.22ym.y1(t)y2(t)11ym.1121p1=11222p2=13.1 两个自由度体系的自由振动 是结构体系的柔度系数(flexibility coefficient)，即体系在点j承受单位力时，在点i产生的位移。 ij27)sin()()sin()(2211tytytyty)sin(t222221112122211111)()(y

17、mymtyymymty 设解的形式为:3.1 两个自由度体系的自由振动设各质点按相同频率和初相角作简谐振动。2222221112212222111121)()()()(ymymyymymy),(21yy),(222112ymym 主振型的位移幅值 就是结构体系在此主振型惯性力幅值 作用下所引起的静力位移。280)1(0)1(222221121221212111ymymymym 令系数行列式等于零，可得到 和 的非零解，即:1y2y01122221212122111mmmmd2222221112212222111121)()()()(ymymyymymy用柔度系数表示的频率方程或特征方程。3.1

18、 两个自由度体系的自由振动292)(4)()(2)(4)()(212112221122221112221112212112221122221112221111mmmmmmmmmmmm求得固有圆频率的两个值为:1112212101122221212122111mmmmd解出 的两个根3.1 两个自由度体系的自由振动体系频率的数目总等于其自由度数目30体系的第一阶主振型:21111212)1(2)1(11mmyy0)1(0)1(222221121221212111ymymymym13.1 两个自由度体系的自由振动2体系的第二阶主振型:22111212)2(2)2(11mmyy31例:求简支梁的自振

19、 频率和主振型。l/3l/3l/3解：1）求柔度系数 p=1 p=132l32leil243432211eil4867321122)(4)()(2121122211222211122211112mmmmmm2)(4)2(222122112111112mmmmm121112mm121112eimlmm31211148615eimlmm3121124861322311221,69. 51mleimlei求得频率：求得主振型：1111112122111mmyy1121112122212mmyymm32例: 求简支梁的自振 频率和主振型。l/3l/3l/3mml/3另解：如果结构本身和质量分布都是对称

20、的，则主振型不是对称就是反对称。故可取半边结构计算 ：1对称情况：eil1625311311169. 51mleiml/91反对称情况：eil4863223222221mleim33例：求图示体系对称振动情况下的频率。mmmeieiei3m3m3m3mm/2m 1210.51m110.8750.252m 1102m01m33eimm5 . 4:,11011相乘eimmmm125. 1,2112021012相乘，相乘或eimm6875. 1:,22022相乘ei5 . 411ei125. 12112ei6875. 1222)(4)()(2121122211222211122211112mmmmm

21、mmeimei943. 01,596. 012211)125. 16875. 15 . 4(214)6875. 125 . 4()6875. 125 . 4(22212eimeimeim/125. 1/8125. 2123412/825. 2/5 . 2/125. 111112122111eimeimeimmmyy11/125. 1/5 . 2/125. 121112122212eimeimeimmmyy2111 yij为正时为正时表示质量表示质量mi的的运动方向与计运动方向与计算柔度系数时算柔度系数时置于其上的单置于其上的单位力方向相同，位力方向相同，为负时，表示为负时，表示与单位力方向与单

22、位力方向相反。相反。35例例3.3.求图示体系的频率、振型求图示体系的频率、振型解解: :eil31134令令21111m02/18/34/31032/9)2/1)(1 (1637. 0336. 1213231140. 2;749. 0mleimleimleimeil1y2y12ym222ym1y2y11211122211221212112yymym2222212212yymym02)1 (211121yy0)2(2112211121yylleil3211221eil3223136例例3.3.求图示体系的频率、振型求图示体系的频率、振型解解: :令令21111m1221212112yymym2

23、222212212yymym02)1 (211121yy02/18/34/31032/9)2/1)(1 (1637. 0336. 1213231140. 2;749. 0mleimleimleimeil1y2y12ym222ym1y2y1121112221ll23. 214/312111yy897. 014/322212yy 1897. 0;123. 221yy37mleimeil1y2y12ym222ym1y2y1121112221ll例例3.3.求图示体系的频率、振型求图示体系的频率、振型解解: :令令21111m1221212112yymym2222212212yymym02)1 (21

24、1121yy23. 214/312111yy897. 014/322212yy 1897. 0;123. 221yy123. 2 1y1897. 0 2y380.5a例例4. 4. 试求图示梁的自振频率和主振型，梁的试求图示梁的自振频率和主振型，梁的eiei已知。已知。12aaamm解：（解：（1 1）计算频率）计算频率1a1m12meiaeiaeia6,4,322321123113231203. 3967. 0maeimaei（2 2）振型）振型61. 31277. 0122122111yyyy10.27713.61第一振型第一振型第二振型第二振型39例5 求图示体系对称振动情况下的频率。m

25、mmaaaaaaeieiei解： 因为结构和质量分布均匀对称，其振型也是对称和反对称的，分别取半边结构计算。对称结构m2=m/2m1=maaaeiei40 以求对称振型为例说明 中系数的求解。首先求出半边结构在集中质量上分别作用有单位集中力产生的弯矩图。aaa1320a1740aaaa710a310a1(a) m1图(b) m2 图 为了求柔度系数，可以在另外的静定基本结构上加单位力并作弯矩图。(c)图1m图2m(d)1aaa2a1aaaa41总结：在2个自由体系自由振动问题中4）体系的自振频率和主振型是体系本身的固有性质。自振频率只与体系本身的刚度系数和质量分布有关，与外荷无关。3）每个自振

26、频率有其相应的主振型；2）振动频率个数与自由度个数一致，自振频率可通过特征方程计算；1）主要问题是确定体系自振频率及其相应主振型3.1 两个自由度体系的自由振动42y1yiynri动平衡方程：riy1yiynri 应满足刚度方程),.,2 , 1(.2211niykykykrniniiikij是结构的刚度系数，使点j产生单位位移（其它点位移为零）时在点i所需施加的力。iiym.),.,2 , 1(0nirymiii. 3.2 多自由度体系的自由振动430.0.0.2211222212122121211111nnnnnnnnnnnykykykymykykykymykykykym.),.,2 ,

27、1(.2211niykykykrniniii或： ),.,2 , 1(0nirymiii.0ykym.0.00.212122221112112121nnnnnnnnnyyykkkkkkkkkyyymmm. 多自由度体系自由振动方程为: 3.2 多自由度体系的自由振动44得振幅方程：2 sin()yyt 3.2 多自由度体系的自由振动0ykym km线性二阶常微分方程组，矩阵 和 是时不变和对称的方阵。 )sin(tyy设方程的解的形式为: y 是位移幅值向量:tnyyyy,210)(2ymk45得频率方程： 3.2 多自由度体系的自由振动 齐次方程，有非零解的条件为系数行列式等于零，即:02m

28、k n个自由度的结构体系，上式展开后得到关于圆频率 的n次代数方程:2), 2 , 1(niitn,21求方程的n个根 ，得到体系的n个自振频率 。1最小的频率 称为基本频率。把全部自振频率按照从小到大的顺序进行排列而成的向量，称为频率向量。46 tniiiiyyyy,21)(0)()(2iiymk), 2 , 1()(niyi线性齐次方程组，如果 是方程组的解：tniiiiyyyy,21)(计算相应的主振型向量: 令： 3.2 多自由度体系的自由振动i对于第i个频率 ，有: i)(iy 为与频率 相应的主振型向量。tniiiicycycycy,21)(也是方程组的解（这里c为任一常数）。47

29、 1)从特征问题中解得的振型的幅值是任意的， 2)只有主振型的形状是唯一的。为了使主振型具有确定值，可以通过以下几种振型标准化（mode normalization）的方法来进行处理。 3.2 多自由度体系的自由振动1）、指定某元素为1 规定主振型向量的某个元素为1，例如取第一元素的值为1，这样以这个元素为标准就可确定其他元素的大小。2）、指定最大元素为1 取每一振型向量的最大值为1，即可确定其他元素的大小。48例： 质量集中在楼层上， 层间侧移刚度如图。k11=4k/3解：1）求刚度系数：m2mmk3k5kk21=-k/3k31=0k12=-k/3k22=8k/15k32=-k/51k13=0k23=-k/5k33=k/5 刚度矩阵k和质量矩阵m：100010002330385052015mmkk1149215,03303850522015kmk其中解得：mk0862. 021mk4453. 022mk8685. 023mk2936. 01mk6673. 02mk9319. 032）求频率：代入频率方程： k2 m0322422252250展开得：1231.293,6.680,13.027503）求主振型：振型方程：（k2 m）y0的后两式：（令y3i=1）0)3(303)8(5221iiiiiyyy（a）122025058300331i