《带滤频的逆迭代法》由会员分享，可在线阅读，更多相关《带滤频的逆迭代法（4页珍藏版）》请在金锄头文库上搜索。
1、带滤频的逆迭代法结构动力分析时常需要求解体系的主振型和频率，当动力自由度较多时，需要采用迭代法进行电算分析。由不考虑阻尼的自由振动方程可得：（4-91）12XMK这在数学上称为广义特征值问题，两边前乘 或 ，可以得到标准特征值问题：1K1M（4-92 ）121或（4-93 ）1XK迭代法先通过假设一个初始的振动形状，然后通过逐步调整振动形状，直到前后两个振动形状充分接近时为止，这样可以得到主振型，再由主振型获得对应的频率。采用（4-92）式进行迭代求解的方法称为逆幂迭代（或逆迭代） ，采用（4-93）式进行迭代求解的方法称为幂迭代。逆幂迭代是求解特征向量和特征值问题的很有效的方法。1计算过程（
2、1）求第一频率和主振型在（4-92）式中，记 、 = ， 称为动力矩阵，于是（4-92 ）式成为求 的特征1MKD21DD值问题：= （4-94 ）DX当矩阵阶数较高时， 不易求，可改用柔度矩阵 ，其系数 可用单位力法求得，系数 表示仅1 ijij在 i 层（质点）作用一单位力时在第 j 层（质点）产生的位移，如在第 n 层作用单位力，则 = , =n1kn2+ = + + 。1k2n1k2nk先假定一初始振型向量，如全为 1 的向量 =1 1 1T，按下式得到一个新的向量：)0(X= )1(XD)0(X（4-95 ）同时将 标准化：)1(= )1(X)1(X（4-96 ）标准化的方法可规定某
3、个元素为 1（如顶层为 1） ，其余取与该元素的比值；或按下式进行：= )1(X5.0)1()1()(XMXT（4-97）再以 为主振型的第二次近似值，前乘动力矩阵 ，得到：)1( D)2(2)1()2( XXDX （4-98）重复上述过程，到第 k 次迭代： （4-99 ）)()1()( kkk XD比较 与 ，如果二者很接近，则迭代结束，得到第一主振型)(kX)1(k 1X和第一频率)(kXk11（2）求高阶频率和振型采用迭代法求解高阶振型和频率在于在迭代过程中消除前面低阶的振型成分，这个过程称为滤型或滤频。任一向量 可按主振型展开：Y，两边同乘以 ：niiXa1 MXTj（4-100）n
4、k kTjTjM1由振型正交性原理得到：（4-101）*jTjj Ya式中， 为广义质量。jTjj X从 中滤去第一主振型成分，余下的振型向量为：（4-102 ）)/(/ *11*111)2( YMXIMYYY TT令 ，称为一阶滤型矩阵， （4-102）式成为：*1 /MIQT（4-103）)2(此时， 中已不包含一阶主振型分量。)2(将 代入（4-94）式右端，得到：)(Y（4-104）1 21)2()2(2 YDQD式中的 就是迭代求解第二阶振型所用的动力矩阵：2（4-105 ）*12*112 / MXMXTT 类似地，可得到 阶滤频后的动力矩阵：s（4-106）*21 ssTssD迭代
5、过程为： 设初始向量 =1 1 1T 或其它值)0(X 求广义质量 sss XM* 求第 阶滤频后的动力矩阵1s*21 ssTss MXD)(1)0(1ssDX21221 ss X )(1)1(11 kskkssks 当 与 很接近时，迭代结束，得到第 阶振型 和第 阶频率)(X)(k 1s1sX)(1ksks12.分析示例三层刚架结构的 13 层的集中质量分别为 360、270 和 180t，侧移刚度分别为 294103、19610 3 和98103kN/m。求结构的振动频率和振型。质量矩阵为： 1802736M刚度矩阵为： 904964K柔度矩阵为： （m/kN ）617.85.33 （1
6、）求第一阶频率和振型 1MD306.29.4.10设初始向量 =1 1 1T)0(X第 1 次迭代： = =6.886 =)0(1)1( XD31086.57420.17349)1(1X第 2 次迭代： = =5.539 =)2(1X1)1( 39.623 )2(1第 3 次迭代： = =5.294 =)3(11)2(4.586310310.5)3(1X第 4 次迭代： = =5.238 = /)4(1X1D)3(.33.649)4(1第 5 次迭代： = =5.227 =)5(11)4(27.59130310.2)5(15X如果振型精度为 ，则迭代可以结束，第一振型 ，第一频率 。40 1X
7、.64983.27.501（2）求第二阶频率和振型设初始向量 =1 1 1T)0(2X求广义质量 56.32 11*1 XMXMT求第 2 阶滤频后的动力矩阵：-*1212 DT 36.29.4.5018356.21073.492= x10-31649.038073 85.09.56.347第 1 次迭代： = =-0.535 =)0(2)1(2XD.13310.429)1(2X第 2 次迭代： = =1.196 =)1(22)2( X96.35310310.5279)(2X第 6 次迭代： = =1.144 =)5(2)6(2XD14.690783310.682)6(26X第二阶振型 ，第二
8、阶频率 。第 9 次迭代得到振型为 ，频2X0.138.4.2 0.1578率为 29.57。（3）求第三阶频率和振型设初始向量 =1 1 1T)0(3X求广义质量 62.4522*2 XM求第 3 阶滤频后的动力矩阵：= -2.566*10-6*2223 DT 8.09.1.037310.165821605.=18076 2.4. 66573第 1 次迭代： = =-0.034 =)0(3)1(3XD4.130310.6827)1(3X第 2 次迭代： = =-0.98 =)1(3)2(398.67233.)2(32第 3 次迭代： = =0.508 =)2(3)3(X50.31310.874)3(3X第 4 次迭代： = =0.53 =)3()4(3D.9633.2)4(34第 5 次迭代： = =0.53 =)4(3)5(3X50.183310.87)5(3X第三阶振型 ，第三阶频率 ，周期 s3X0.18724.314.023T
