演示文稿振动力学第五章

演示文稿振动力学第五章当前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页。当k>1时，