演示文稿振动力学第五章

演示文稿振动力学第五章现在是1页\一共有59页\编辑于星期三（优选）振动力学第五章现在是2页\一共有59页\编辑于星期三振动速度系统的动能将振动速度代入得动能的最大值发生在时刻，即现在是3页\一共有59页\编辑于星期三若只考虑弯曲变形的影响，系统的应变能为将运动方程代入得当时，应变能最大，即使，即可得到瑞利商现在是4页\一共有59页\编辑于星期三用外力做功的数值代替系统应变能的数值图（b）系统上外力所做的总功为将运动方程代入上式得y(x,t)为静荷载（自重、F等）引起的位移，如自重等现在是5页\一共有59页\编辑于星期三

当时，应变能达到最大值，此时外力所作的功亦为最大值，

这时系统的动能除了分布质量m(x)的动能外，还应包括各集中质量的动能，即将振动速度代入得现在是6页\一共有59页\编辑于星期三当时，动能达最大值由得到现在是7页\一共有59页\编辑于星期三例：如图（a）所示均质等截面简支梁。单位梁长的质量为，其抗弯刚度EI为常数。若振型分别为图（b）所示（为梁中点的最大挠度）和图（c）所示梁在自重作用下的挠曲线。分别计算自振频率，并将所得结果进行比较。现在是8页\一共有59页\编辑于星期三解：（1）振型为从而得自振频率精确解现在是9页\一共有59页\编辑于星期三（2）取振型为梁在自重荷载上的挠曲线。图（c）所示为匀布自重荷载作用下简支梁的静力挠曲线，即最大动能外力做功的最大值现在是10页\一共有59页\编辑于星期三因为,可以解得此值与精确解相比较，偏大约2％现在是11页\一共有59页\编辑于星期三例：计算重力坝沿水流方向的自振频率时，可以取沿坝轴线方向单位长度的坝体近似地简化为图（a）所示的变截面悬臂梁。试用瑞利法计算其自振频率。现在是12页\一共有59页\编辑于星期三解：选变截面悬臂梁在其自重作用下所引起的挠曲线作为近似振型，如图（b）所示，即从图（b）可以看出其分布质量为最大动能和外力功的最大值为现在是13页\一共有59页\编辑于星期三根据得到现在是14页\一共有59页\编辑于星期三例：等截面悬臂梁端部有一集中质量用瑞利法估计基频解：选择等截面悬臂梁在均布载荷下的静挠度曲线作为试函数：选择端部集中质量作用下的静挠度曲线作为试函数：因集中质量大于梁的分布质量，选用后一种试函数好现在是15页\一共有59页\编辑于星期三例.用能量法计算图示体系的基频.mmm321解:1.取自重引起的位移mgmgmg精确解:现在是16页\一共有59页\编辑于星期三2.取直线mmm321mgmgmg3.取常数精确解:现在是17页\一共有59页\编辑于星期三二，李兹能量法李兹给出了级数形式的近似振型将上式代入瑞利商的表达式得现在是18页\一共有59页\编辑于星期三引进下列记号为所以根据频率为极值的条件现在是19页\一共有59页\编辑于星期三得到即简化上式并将代入得或现在是20页\一共有59页\编辑于星期三

上式为n个齐次线性方程，为了使方程组有非零解，必须得到

上式展开后得到一个的n次方程，该方程有n个根。对于其中的每一个根都可求得一组常数，因此得到n个振型函数

求得的就是所研究的系统前n个自振频率和振型函数的近似解。现在是21页\一共有59页\编辑于星期三例：试用李兹法求图所示重力坝的第一和第二阶自振频率。现在是22页\一共有59页\编辑于星期三解：为了使级数各项都满足位移边界条件，近似振型函数选为假设经一次近似计算只取第一项，即代入瑞利商的表达式得现在是23页\一共有59页\编辑于星期三若取级数前两项，即

将代入相关式子计算出，这时成为展开系数行列式，并令其等于零,得频率方程：现在是24页\一共有59页\编辑于星期三解得与精确解的相对误差为0.6％，是较高一阶频率的近似值。现在是25页\一共有59页\编辑于星期三例：图所示等截面悬臂梁，用李兹法求自振频率。现在是26页\一共有59页\编辑于星期三解：选取两个函数：

这两个函数在梁的支承处满足固定端边界条件。于是，近似振型函数可取为求得如下现在是27页\一共有59页\编辑于星期三现在是28页\一共有59页\编辑于星期三于是，频率方程为从上式可得到一个关于的方程，方程的根为现在是29页\一共有59页\编辑于星期三这两个频率的精确值为比较得，第二阶自振频率精读很差。为了改善得计算精读，采用以下四个函数：现在是30页\一共有59页\编辑于星期三求得结构的前四阶频率为该结构第三阶和第四阶自振频率的精确值为

比较得，的精读最高，次之，的精读最差。所以说，为了得到精读较高的高阶频率，往往需要选取较多的函数。现在是31页\一共有59页\编辑于星期三例：等截面简支梁梁中部有一集中质量Ma，大小等于梁的质量采用里兹法，求：梁的模态函数近似解Ma选取无集中质量时的梁的模态函数作为里兹基函数：解：基函数满足自然边界条件（两端挠度和弯矩为零）现在是32页\一共有59页\编辑于星期三模态试函数：若对第三阶固有频率的精度要求不高，取n＝3代入里兹法方程，求得系数：现在是33页\一共有59页\编辑于星期三模态试函数：梁的模态函数近似解：现在是34页\一共有59页\编辑于星期三

第二节幂法计算自振频率和振型一，最低阶频率和振型的计算上式两边左乘以

首先假定一个任意的规准化振型，例如令其中第一个自由度振幅值为1，即，亦即，假定规准化振型上标（0）表示假设的初始形状，即零次迭代。现在是35页\一共有59页\编辑于星期三

把这个假定的标准化振型代入等号左边，经过运算得，即将中第一个元素归一化为1后，得式中

这就是频率和振型的第一次近似值。再把代入，仿此继续循环迭代计算，直到经过连续迭代前后两次的给出相同或近乎相同的数值为止，这样得到的就是系统的最低自振频率及相应的振型。现在是36页\一共有59页\编辑于星期三如果假定的形状是一个真实的振型，则因此，

那么，中任何一对对应元素的比值都能得到相同的频率，则

一般来说，经过一次迭代后的和假定的的形状是不一样的。对于上式中的每一次位移坐标会得到不同的值。现在是37页\一共有59页\编辑于星期三在这种情况下，

为了求出较好的频率近似值，通常采用以质量作为加权系数的平均法，用左乘以

当真实振型或是自重作用下的挠曲线都不能很快给出时，习惯上总是把各质体的幅值假定为1，即取现在是38页\一共有59页\编辑于星期三例：如图所示三层刚架，试用幂法计算它的最低自振频率和振型。现在是39页\一共有59页\编辑于星期三解：该系统的劲度矩阵和质量矩阵分别为因此，这个结构的柔度矩阵是现在是40页\一共有59页\编辑于星期三由此得将假定的初始迭代振型代入上式等号左边，得现在是41页\一共有59页\编辑于星期三将代入，算得将代入，算得现在是42页\一共有59页\编辑于星期三将代入，算得

因为前后两次迭代的振型基本相同，迭代至此停止，求得的第一振型为，相应的自振频率为精确解现在是43页\一共有59页\编辑于星期三如果按照来求第一自振啤频率，则

可见，第一次迭代求得的频率精度较差。采用质量为权系数平均，只需迭代一次就能得到较好的频率近似值现在是44页\一共有59页\编辑于星期三

现在来证明上述迭代法求出频率和振型就是系统的最低自振频率和相应的振型。

当给出假定的振型后，逐次迭代可以作出如下的一系列向量现在是45页\一共有59页\编辑于星期三对于开始所假定的振型可表示为将上式前乘以D,则现在是46页\一共有59页\编辑于星期三

由于即，故当迭代次数k充分大时，，只要时，则有

这就说明，在迭代k次后，向量与向量仅相差一常数倍数，或者说向量收敛于向量。由于每迭代一次都要归一化一次，所提出的因子就是，所以迭代k次后，就收敛到系统的第一自振频率和对应的振型。现在是47页\一共有59页\编辑于星期三二，最高阶频率和振型的计算用左乘以

则基于上式的迭代计算结果将得到最高阶的自振频率和相应的振型。为了论证这一点，令按照前面同样的论证方法可以得到迭代k次的向量为现在是48页\一共有59页\编辑于星期三

因为，当k充分大时，，所以上面等式右端各项比最后一项要小得多，略去前面（n-1）项，于是得到

这就说明，迭代k次后就收敛到系统的最高一阶振型。给出第n阶自振频率的近似值或现在是49页\一共有59页\编辑于星期三例：图所示三层刚架，试用幂法计算它的最高自振频率和振型。现在是50页\一共有59页\编辑于星期三解：按，有

假定初始振型，并设，代入上式等号左边，得到现在是51页\一共有59页\编辑于星期三继续迭代计算，得现在是52页\一共有59页\编辑于星期三

前后两次迭代振型已基本接近，迭代中止，得到第三阶振型为与前面所得第三阶振型一致，其相应的第三阶自振频率为现在是53页\一共有59页\编辑于星期三三，高阶频率和振兴的计算假设振型

逐阶消频法：当要求第r+1阶振型时，可以在假设振型中清除掉所有前面r阶振型分量，逐步收敛到第r+1阶振型，从而求出所需若干阶振型。在上式等号两边前乘以,并利用振型的正交性得现在是54页\一共有59页\编辑于星期三从上式中解出

为了在假设的振型中清掉所有前面r阶振型分量，可取初始迭代向量为

式中，为r阶清型矩阵或淘汰矩阵；I为主对角元素为1的对角矩阵。现在是55页\一共有59页\编辑于星期三

在实际迭代计算过程中，应该在每次迭代后都要重新清型。也就是说，只是在求系统的第一阶振型时用矩阵D前乘，在以后各阶振型的计算中，每次都要用清型后的矩阵来前乘。经过清型后的各阶矩阵称为收缩矩阵，可表示为现在是56页\一共有59页\编辑于星期三收缩矩阵还可以写成递推公式的形式对收缩矩阵作些说明。上式取r=1，则上式两边右乘以现在是57页\一共有59页\编辑于星期三当