《振动力学第五章.ppt》由会员分享，可在线阅读，更多相关《振动力学第五章.ppt（61页珍藏版）》请在金锄头文库上搜索。
1、第五章 自振频率和振型的实用计算,第一节 能量法求自振频率,根据能量守恒，在任何瞬时（忽略能量散失）,一，瑞利能量法,设图示系统中任一质点的运动方程为,振动速度,系统的动能,将振动速度代入得,动能的最大值发生在 时刻，即,若只考虑弯曲变形的影响，系统的应变能为,将运动方程代入得,当 时，应变能最大，即,使 ，即可得到,瑞利商,用外力做功的数值代替系统应变能的数值 图（b）系统上外力所做的总功为,将运动方程代入上式得,y(x,t)为静荷载（自重、F等）引起的位移，如自重等,当 时，应变能达到最大值，此时外 力所作的功亦为最大值，,这时系统的动能除了分布质量m(x)的动能外，还应 包括各集中质量
2、的动能，即,将振动速度代入得,当 时，动能达最大值,由 得到,例：如图（a） 所示均质等截面简支 梁。单位梁长的质量为 ，其抗弯 刚度EI为常数。若振型分别为图（b） 所示 （ 为梁中点的最大 挠度）和图（c）所示梁在自重作用下 的挠曲线。分别计算自振频率，并将 所得结果进行比较。,解：（1）振型为,从而得,自振频率,精确解,（2）取振型为梁在自重荷载上的挠曲线。图（c）所示为匀布 自重荷载作用下简支梁的静力挠曲线，即,最大动能,外力做功的最大值,因为,可以解得,此值与精确解相比较，偏大约2,例：计算重力坝沿水流方向 的自振频率时，可以取沿坝 轴线方向单位长度的坝体近 似地简化为图（a）所示的
3、变 截面悬臂梁。试用瑞利法计 算其自振频率。,解：选变截面悬臂梁在其自重作用下所引起的挠曲线作为 近似振型，如图（b）所示，即,从图（b）可以看出其分布质量为,最大动能和外力功的最大值为,根据,得到,14,例：等截面悬臂梁,端部有一集中质量,用瑞利法估计基频,解：,选择等截面悬臂梁在均布载荷下的静挠度曲线作为试函数：,选择端部集中质量作用下的静挠度曲线作为试函数：,因集中质量大于梁的分布质量，选用后一种试函数好,例.用能量法计算图示体系的基频.,解:,1.取自重引起的位移,精确解:,2.取直线,3.取常数,精确解:,二，李兹能量法,李兹给出了级数形式的近似振型,将上式代入瑞利商的表达式得,引进
4、下列记号为,所以,根据频率为极值的条件,得到,即,简化上式并将 代入得,或,上式为n个齐次线性方程，为了使方程组有非零解，必 须得到,上式展开后得到一个 的n次方程，该方程有n个根。对 于其中的每一个根 都可求得一组常数 ，因 此得到n个振型函数,求得的 就是所研究的系统前n个自振频率和振型 函数的近似解。,例：试用李兹法求图所示重力坝 的第一和第二阶自振频率。,解：为了使级数各项都满足位移边界条件，近似振型函数选为,假设经一次近似计算只取第一项，即,代入瑞利商的表达式得,若取级数前两项，即,将 代入相关式子计算出 ，这时 成为,展开系数行列式，并令其等于零,得频率方程：,解得,与精确解的相对
5、误差为0.6， 是较高一阶频率的近似值。,例：图所示等截面悬臂梁，用李兹法求自振频率。,解：选取两个函数：,这两个函数在梁的支承处满足固定端边界条件。于是， 近似振型函数可取为,求得 如下,于是，频率方程为,从上式可得到一个关于 的方程，方程的根为,这两个频率的精确值为,比较得，第二阶自振频率精读很差。,为了改善 得计算精读，采用以下四个函数：,求得结构的前四阶频率为,该结构第三阶和第四阶自振频率的精确值为,比较得， 的精读最高， 次之， 的精读最差。所 以说，为了得到精读较高的高阶频率，往往需要选取较多 的函数 。,31,例：等截面简支梁,梁中部有一集中质量 Ma，大小等于梁的质量,采用里兹
6、法，求：梁的模态函数近似解,选取无集中质量时的梁的模态函数作为里兹基函数：,解：,基函数满足自然边界条件（两端挠度和弯矩为零）,32,模态试函数：,若对第三阶固有频率的精度要求不高，取 n3,代入里兹法方程，求得系数：,33,模态试函数：,梁的模态函数近似解：,第二节 幂法计算自振频率和振型,一，最低阶频率和振型的计算,上式两边左乘以,首先假定一个任意的规准化振型 ，例如令其中第一 个自由度振幅值为1，即 ，亦即，假定规准化振型,上标（0）表示假设的初始形状，即零次迭代。,把这个假定的标准化振型代入 等号左边，经 过运算得 ，即,将 中第一个元素归一化为1后，得,式中,这就是频率和振型的第一次
7、近似值 。再把 代入 ，仿此继续循环迭代计算，直到经过连续迭代前后两 次的 给出相同或近乎相同的数值为止，这样得到的就是系 统的最低自振频率 及相应的振型 。,如果假定的形状 是一个真实的振型，则,因此，,那么， 中任何一对对应元素的比值都能得到相同的 频率 ，则,一般来说，经过一次迭代后的 和假定的 的形状是不一 样的。对于上式中的每一次位移坐标会得到不同的 值。,在这种情况下，,为了求出较好的频率近似值，通常采用以质量作为 加权系数的平均法，用 左乘以,当真实振型或是自重作用下的挠曲线都不能很快给出 时，习惯上总是把各质体的幅值假定为1，即取,例：如图所示三 层刚架，试用幂 法计算它的最低
8、 自振频率和振型。,解：该系统的劲度矩阵和质量矩阵分别为,因此，这个结构的柔度矩阵是,由此得,将假定的初始迭代振型 代入上式等号左边，得,将 代入，算得,将 代入，算得,将 代入，算得,因为前后两次迭代的振型基本相同，迭代至此停止，求得的 第一振型为 ，相应的自振频率为,精确解,如果按照 来求第一自振啤频率，则,可见，第一次迭代求得的频率精度较差。采用质量为权系数 平均，只需迭代一次就能得到较好的频率近似值,现在来证明上述迭代法求出频率和振型就是系统的最低自振 频率和相应的振型。,当给出假定的振型 后，逐次迭代可以作出如下的一系列 向量,对于开始所假定的振型 可表示为,将上式前乘以D,则,由于
9、 即 ，故当迭代次数 k充分大时， ，只要 时，则有,这就说明，在 迭代k次后，向量 与向量 仅相差一 常数倍数，或者说向量 收敛于向量 。由于每迭代 一次 都要归一化一次，所提出的因子就是 ，所以 迭代k次后，就收敛到系统的第一自振频率和对应的振型。,二，最高阶频率和振型的计算,用 左乘以,则基于上式的迭代计算结果将得到最高阶的自振频率和 相应的振型。为了论证这一点，令,按照前面同样的论证方法可以得到迭代k次的向量为,因为 ，当k充分大时， ，所以上面等式右端各项比最 后一项要小得多，略去前面（n-1）项，于是得到,这就说明，迭代k次后就收敛到系统的最高一阶振型。给出 第n阶自振频率的近似值
10、,或,例：图所示三层刚架， 试用幂法计算它的最高 自振频率和振型。,解：,按 ，有,假定初始振型 ，并设 ，代入 上式等号左边，得到,继续迭代计算，得,前后两次迭代振型已基本接近，迭代中止，得到第三阶振 型为,与前面所得第三阶振型一致，其相应的第三阶自振频率为,三，高阶频率和振兴的计算,假设振型,逐阶消频法：当要求第r+1阶振型 时，可以在假设振型 中清除掉所有前面r阶振型分量，逐步收敛到第r+1阶振型， 从而求出所需若干阶振型。,在上式等号两边前乘以 ,并利用振型的正交性得,从上式中解出,为了在假设的振型中清掉所有前面r阶振型分量，可取初始迭代 向量为,式中， 为r阶清型矩阵或淘汰矩阵；I为
11、主对角元素为1的对角 矩阵。,在实际迭代计算过程中，应该在每次迭代后都要重新清型。 也就是说，只是在求系统的第一阶振型时用矩阵D前乘，在 以后各阶振型的计算中，每次都要用清型后的矩阵来前乘。,经过清型后的各阶矩阵称为收缩矩阵，可表示为,收缩矩阵还可以写成递推公式的形式,对收缩矩阵作些说明。上式取r=1，则,上式两边右乘以,当k1时， 。 的特征值与特征向量和D的相同。,当k1时， 。 的 ，即第一阶特征值和对 应的特征向量给“收缩”掉了。,推至一般而言，收缩矩阵 的 ，而 的第r+1阶以上的特征值和各特征向量都和D的相同。,例：图所示三层 钢架，试用幂法 计算它的高阶自 振频率和振型。,解：计算第一阶频率和振型计算式的标准形式 为,用幂法迭代算得,清除第一阶振型的收缩矩阵 为,计算第二阶自振频率和振型计算式的标准形式 为,用幂法迭代算得,清除第一阶振型和第二阶振型的收缩矩阵 为,计算第三阶自振频率和振型计算式的标准形式 为,用幂法迭代算得,(与精确值基本 一致）,
