第3 章(II)多自由度系统的振动ppt课件.ppt
第3章 多自由度系统的振动李映辉西南交通大学2015.09,2022年11月19日,振动力学,2,2022年11月19日,中国力学学会学术大会2005,2,2022年11月19日,2,声 明,本课件可供教师教学和学生学习中免费使用。不可用于任何商业目的。本课件的部分内容参阅了上海交通大学陈国平教授和太原科技大学杨建伟教授的课件，作者在此向二位教授表示衷心感谢。如该课件无意中损害了二位教授利益，作者在此致歉。本课件以高淑英、沈火明编著的振动力学（中国铁道出版社，2011年）的前四章为基础编写。感谢研究生蒋宝坤、王金梅在文字录入方面的工作,2022年11月19日,振动力学,3,教学内容,多自由度系统的振动,2022年11月19日,振动力学,3,教学内容,两自由度系统的振动多自由度系统的振动多自由度系统固有特性的近似解法,2022年11月19日,振动力学,4,教学内容,多自由度系统的振动/多自由度系统固有特性的近似解法,2022年11月19日,振动力学,4,多自由度系统固有特性的近似解法邓柯莱法（Dunkerley）瑞雷法（Rayleigh)李兹法（Ritz)传递矩阵法矩阵迭代法,2022年11月19日,振动力学,5,多自由度系统的振动/多自由度系统固有特性的近似解法/邓柯莱法,1 邓柯莱法（Dunkerley） 该方法由邓柯莱在用实验确定多圆盘轴横向振动频率时提出，作为系统基频估算公式。 设n自由度系统质量阵、柔度阵为自由振动方程为,2022年11月19日,振动力学,6,多自由度系统的振动/多自由度系统固有特性的近似解法/邓柯莱法,特征方程为 式中=1/2。展开设上式根1=1/12 ，2=1/22 ，, n=1/n2 ,则（3.101）可表为,2022年11月19日,振动力学,7,多自由度系统的振动/多自由度系统固有特性的近似解法/邓柯莱法,展开得比较（3.101）和（3.103）得即因12、 3、n,1/2、1/ 3、1/n较小，得式中ii=1/kii，则,2022年11月19日,振动力学,8,多自由度系统的振动/多自由度系统固有特性的近似解法/邓柯莱法,故上式即为邓柯莱公式，ii是系统在质量mii单独作用下（其他质量为零）系统的固有频率。 因（3.105）左边舍去了一些正项，由（3.105）计算的1/12值比实际大，12实际值小。,2022年11月19日,振动力学,9,多自由度系统的振动/多自由度系统固有特性的近似解法/邓柯莱法,【例3.10】图3.15为等截面简支梁。其有3集中质量是m1、 m2、 m3，梁弯曲刚度为EI，质量不计。用邓柯莱法计算系统第一阶固有频率的近似值。已知： m1= m3 =m, m2 =2m。【解】由材料力学知，简支梁在单位下的挠曲线公式为a、b分别为力作用点到左右端的距离。,2022年11月19日,振动力学,10,多自由度系统的振动/多自由度系统固有特性的近似解法/邓柯莱法,求得柔度影响系数为由（3.107）得求得的1值比精确值小2.5%。,2022年11月19日,振动力学,11,2022年11月19日,振动力学,11,教学内容,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,2022年11月19日,振动力学,11,多自由度系统固有特性的近似解法邓柯莱法（Dunkerley）瑞雷法（Rayleigh)李兹法（Ritz)传递矩阵法矩阵迭代法,2022年11月19日,振动力学,12,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,2 瑞雷法（Rayleigh) 多自由度系统的动能T与势能U的表达式为系统作某一阶主振动时代入（3.108）和(3.109)得系统在作第i阶主振动时，最大动能Tmax与最大势能Umax,2022年11月19日,振动力学,13,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,由机械能守恒定律， Tmax = Umax 得在（3.115）中A(i)代入假设振型A，结果以R1表示，则称为瑞雷商。这种计算系统固有频率的方法称为瑞雷法. 因很难选A(i)接近高阶主振型，通常不用瑞雷法求高阶固有频率，只用它求低阶固有频率。,2022年11月19日,振动力学,14,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,取接近一阶主振型的假设振型A代入（3.115），则瑞雷商为一阶固有频率平方近似值。证明如下： 如假设振型A不是主振型，将其用正则振型线性表示 有瑞雷商为,2022年11月19日,振动力学,15,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,若A接近于一阶主振型A（1），则C2/C11, C3/C11, Cn/C11, 则由（3.118）有一般以静变形作假设振型，可得相当准确的结果。如选取A有困难，可任选一A。与动力矩阵D(=M)相乘，得B1=DA，然后以B1或与其成比例的B1作A(1)的近似振型，再按（3.116)计算R1，可得12好的近似。,2022年11月19日,振动力学,16,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,瑞雷法也可用于以柔度阵建立振动方程的情况，这时系统势能U等于外力的功，即在振动过程中，只有惯性力作用，即因x为得,2022年11月19日,振动力学,17,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,势能、动能最大值由Tmax= Umax ，得当A为第i阶主振型，由（3.122）得第i阶固有频率的平方值i2。在（3.122）中代入假设振型A，结果用R2表示，则有上式称为第二瑞雷商。,2022年11月19日,振动力学,18,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,注意：用（3.115）或（3.123）计算的12总比精确值12大。因任选一A，即对系统增加了约束，提高了系统刚度，使频率增大。【例3.11】用瑞雷法求例3.10中一阶固有频率的近似值。【解】由例3.10系统质量阵M和柔度阵为,2022年11月19日,振动力学,19,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,三点处静挠度为最大势能为各点最大速度为y1、y2、y3，最大动能为,2022年11月19日,振动力学,20,多自由度系统的振动/多自由度系统固有特性的近似解法/瑞雷法,由（3.124）得此结果比真实值略高，误差为0.02%。,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,2022年11月19日,振动力学,21,作业,第97页3.19,多自由度系统的振动/多自由度系统的振动,2022年11月19日,振动力学,22,2022年11月19日,振动力学,22,2022年11月19日,振动力学,22,教学内容,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,2022年11月19日,振动力学,22,多自由度系统固有特性的近似解法邓柯莱法（Dunkerley）瑞雷法（Rayleigh)李兹法（Ritz)传递矩阵法矩阵迭代法,2022年11月19日,振动力学,23,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,3 李兹法（Ritz) 瑞雷法理论上可求系统的各阶固有频率，但实际上难以用于求高阶固有频率。 李兹法对瑞雷商进行了改进，采用其极值形式，能找到较精确的低阶和高阶振型，不仅可求出更精确的基频，还可计算高阶频率和振型，故李兹法也称为瑞雷李兹法。 李兹法需先假定若干振型，并进行线性组合，用瑞雷法计算前几阶固有频率。若系统自由度很大，矩阵阶数很高，其存储量大，运算速度慢。,2022年11月19日,振动力学,24,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,如希望有s阶频率与振型为准确值，需假设n1个振型(2sn1n) ，矩阵阶数大为降低，故李兹法是一种缩减系统自由度的近似解法。介绍如下： 任取n1个线性无关的特征向量1 、 2 、 n1 ，用其线性组合作为新的假设振型A,即式中C1 、 C2 、 Cn1为待定常数，将（3.125）表为矩阵形式其中,2022年11月19日,振动力学,25,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,将A=C代入（3.116）得由瑞雷商极值性质，可得待定常数Cj，即令将这n1个方程表为矩阵形式其中 分别为n1n1的广义刚度阵和广义质量阵.,2022年11月19日,振动力学,26,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,（3.130）为特征值问题，因阶数n1远小于系统自由度数n，求解简便。 由（3.130）得s个特征值R1 、R2 、Rn1为系统最低的n1个固有频率，C(1)、 C(2) 、 C(n1) 为对应的n1个主振型【例3.12】图示弹簧质量系统,用李兹法求其前三阶固有频率和主振型。,2022年11月19日,振动力学,27,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,【解】取广义坐标x1 、x2 、 x3、x4 ，系统质量阵、刚度阵为设则广义刚度阵为,2022年11月19日,振动力学,28,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,广义质量阵为代入（3.130）得解得,2022年11月19日,振动力学,29,多自由度系统的振动/多自由度系统固有特性的近似解法/李兹法,系统最低二阶固有频率的近似值为主振型的近似值为同样可以求出另两阶频率和振型。,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,2022年11月19日,振动力学,30,作业,第98页3.22,多自由度系统的振动/多自由度系统的振动,2022年11月19日,振动力学,31,2022年11月19日,振动力学,31,2022年11月19日,振动力学,31,教学内容,多自由度系统的振动/多自由度系统固有特性的近似解法,2022年11月19日,振动力学,31,多自由度系统固有特性的近似解法邓柯莱法（Dunkerley）瑞雷法（Rayleigh)李兹法（Ritz)传递矩阵法矩阵迭代法,2022年11月19日,振动力学,32,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,4.传递矩阵法 质量阵、刚度阵形成后，前述方法有广泛应用。 对一环连一环，呈链状的工程结构（如发动机曲轴、连续梁等），可采用另一计算方法传递矩阵法。 该法只需对低阶次的传递矩阵进行乘法运算，并计算其特征值，节省计算工作量。 由界面上的力和位移组成列矢量，称为状态矢量。联系相邻单元间状态矢量的矩阵，称传递矩阵。传递矩阵把状态矢量从一个位置转换到另一位置，因此称为传递矩阵法，又称变换矩阵法。,2022年11月19日,振动力学,33,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,（1）梁上有集中质量的横向振动系统 连续梁或汽轮机的发动机转子可简化为无质量的梁附加若干集中质量的横向振动系统，如图3.17（a）。图 3.17,2022年11月19日,振动力学,34,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,设第i个单元集中质量mi，梁长li，抗弯刚度EIi。梁段及集中质量受力如图3.17（c）、（d）。各截面挠度y、截面转角、剪力Q及弯矩M约定为正值。任一截面状态向量有4个分量，即广义位移y与及广义力Q与M，表示为,2022年11月19日,振动力学,35,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,由图3.17（d）的力平衡条件有图 3.17,2022年11月19日,振动力学,36,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,设第i段梁上距左端x处截面的弯矩、剪力、转角、挠度分别为Mi(x)、Qi(x) 、i(x) 及yi(x) ，则上式中令x=li,得,2022年11月19日,振动力学,37,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,表为矩阵形式简写成式中Hi f称为场传递矩阵。 由图3.17（c），集中质量两边的挠度、转角、弯矩、剪力满足,2022年11月19日,振动力学,38,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,当系统以频率振动时，（3.142）为（3.139）-（3.143）表示矩阵形式,2022年11月19日,振动力学,39,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,HiP称为点传递矩阵。 由（3.138）、（3.145）得Zi-1R到ZiR的传递关系为其中Hi称为第i单元的传递矩阵,2022年11月19日,振动力学,40,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,由此得到记总传递矩阵为则从最左端与最右端间的传递关系为,2022年11月19日,振动力学,41,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,具体形式为H中各元素依赖于，表示为 两端边界条件已知，可以得到梁的固有频率。,2022年11月19日,振动力学,42,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,（2）轴盘扭转振动系统 对图3.18（a）示链状轴盘扭振系统，其典型单元包括无质量轴段和刚性圆盘。 设第i单元内轴段扭转刚度ki，长度li ，圆盘转动惯量Ji ，轴段及圆盘受力如图3.18（b）-（c）。 各截面转角、扭矩M约定为正值。截面状态向量为不计轴段转动惯量，由图3.18知两边扭矩相等,2022年11月19日,振动力学,43,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,轴段两边转角满足：,2022年11月19日,振动力学,44,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,（3.154）与（3.155）写成矩阵形式即表示从状态向量Zi-1R到ZiL的传递关系，Hif称为场传递矩阵。 由图3.18（c）知圆盘两边转角相等，即圆盘振动方程为,2022年11月19日,振动力学,45,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,当轴盘系统以频率振动时，有 ，代入得（3.158）与（3.157）的矩阵形式得表示从状态向量ZiL到ZiR的传递关系，HiP称为点传递矩阵。,2022年11月19日,振动力学,46,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,由（3.156）、（3.160）得从ZiR到Zi-1R的传递关系Hi称为第i单元的传递矩阵。等于Hi是频率的函数。通过各单元的传递矩阵，可建立结构最左端与最右端的状态向量之间的传递关系。,2022年11月19日,振动力学,47,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,【例3.13】3圆盘扭振系统。各圆盘的转动惯量J1=4.9kg.m2，J1=4.9kg.m2 ，J3=19.6kg.m2 ，轴段扭转刚度k2=98103N.m， k3=196103N.m .用传递矩阵法求各阶固有频率和主振型。【解】由传递矩阵法，系统最左端与最右端状态向量之间传递关系为即,2022年11月19日,振动力学,48,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,边界条件为M1L=M3R=0。对左端边界条件，在（a）中M1L =0，得因1L任意，是固有频率，M3R=0，代入（b），有上式即为频率方程。设最左端状态向量为,2022年11月19日,振动力学,49,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,Z1R=H1PZ1L, Z2R=H1fZ1R, Z3R=H3Z2R的具体形式为可得该扭转系统的固有频率和模态（注意刚体运动）,2022年11月19日,振动力学,50,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,对应主振型为,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,2022年11月19日,振动力学,51,作业,第98页3.24,多自由度系统的振动/多自由度系统固有特性的近似解法/传递矩阵法,2022年11月19日,振动力学,52,2022年11月19日,振动力学,52,2022年11月19日,振动力学,52,2022年11月19日,振动力学,52,教学内容,多自由度系统的振动/多自由度系统固有特性的近似解法,2022年11月19日,振动力学,52,多自由度系统固有特性的近似解法邓柯莱法（Dunkerley）瑞雷法（Rayleigh)李兹法（Ritz)传递矩阵法矩阵迭代法,2022年11月19日,振动力学,53,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,5.矩阵迭代法（1）求一阶固有频率和振型 对无阻尼多自由度振动系统，其固有频率及主振型可由下式求出：上式也可写为引入 D=M 和 =1/2，则（3.164）可写为D称为系统的动力矩阵。,2022年11月19日,振动力学,54,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,用矩阵迭代法是计算固有频率和主振型步骤如下：（1）任假设一初始振型A；（2）按下格式计算位形列阵序列 Am，m=1,2,；当m足够大时，位形列阵趋于第一主振型，即 且,2022年11月19日,振动力学,55,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,式中（ai）m为列阵Am=a1 a2 aimT的第i个元素。证明如下： 因任意初始振型A0可表为主振型的线性组合，即第一次迭代得同样,2022年11月19日,振动力学,56,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,继续迭代有当迭代次数较大时，(2/ 1 )m， (3/ 1 )m ， (n/ 1 )m 均小于1，除了第一项外，其他各项均可忽略.# 影响迭代收敛速度因素：（1）系统本身性质，即2/ 1的大小；（2）初始列阵A0，越接近一阶主振型，即C2/ C1 、 C3/ C1 、等越小，收敛越快。 迭代收敛速度主要取决于2/ 1的比值。,2022年11月19日,振动力学,57,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,(2)求高阶固有频率和振型 矩阵迭代法可求出全部固有频率及主振型。下面说明求各阶固有频率及主振型的清除矩阵迭代法。 设第一阶固有频率1和主振型A(1)已求得,将A(1)对质量阵正则化得AN(1) 。为求二频率，构造如下动力矩阵：D2称为清除矩阵（清除第一振型的动力矩阵）。 用上述迭代步骤，任取初始振型A0，用D2替代原来D，则迭代将收敛于2及A(2) 。,2022年11月19日,振动力学,58,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,证明如下：由于第一次迭代注意到得到,2022年11月19日,振动力学,59,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,可见，第一次迭代不含第一主振型，迭代将收敛于第二主振型和二阶固有频率。 类似在求得1、 2和AN(1)、 AN(2)后，可构造三阶固有频率清除矩阵D3：迭代可得AN(3)和3。,2022年11月19日,振动力学,60,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,【例3.14】用矩阵迭代法解例3.10中系统的基频和主振型，【解】已知得动力矩阵式中,2022年11月19日,振动力学,61,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,设初始振型A0=1 1 1T，第一次迭代得：第二次迭代继续迭代有,2022年11月19日,振动力学,62,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,因，A3 A4，迭代停止，得第一阶主振型和基频为,2022年11月19日,振动力学,63,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,【例3.15】求例3.14的第二阶固有频率和主振型。【解】将上例中A(1)=1.000 1.4277 1.000T对M正则化，得由（3.168），有,2022年11月19日,振动力学,64,2022年11月19日,振动力学,64,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,2022年11月19日,振动力学,65,多自由度系统的振动/多自由度系统固有特性的近似解法/矩阵迭代法,因二阶主振型将有一节点（？），取初始列阵为按A m= D2 Am-1 （m=1,2,）迭代，迭代10次的结果为得即主振型为,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,2022年11月19日,振动力学,66,作业,第98页3.24,多自由度系统的振动/多自由度系统的振动,
