带滤频的逆迭代法借鉴.pdfVIP

带滤频的逆迭代法 结构动力分析时常需要求解体系的主振型和频率，当动力自由度较多时，需要采用迭代法进行电算分 析。由不考虑阻尼的自由振动方程可得： 1 2 [ K ]{ X } [ M ]{ X } （4-91 ） 这在数学上称为广义特征值问题，两边前乘 [ K ] 1 或 [ M ] 1 ，可以得到标准特征值问题： [ K ] 1 [ M ]{ X } 1 { X } （4-92 ） 2 或 1 2 [ M ] [ K ]{ X } { X } （4-93 ） 迭代法先通过假设一个初始的振动形状，然后通过逐步调整振动形状，直到前后两个振动形状充分接 近时为止，这样可以得到主振型，再由主振型获得对应的频率。采用（ 4-92 ）式进行迭代求解的方法称为 逆幂迭代（或逆迭代） ，采用（ 4-93 ）式进行迭代求解的方法称为幂迭代。逆幂迭代是求解特征向量和特征 值问题的很有效的方法。 1．计算过程 （1）求第一频率和主振型 在（ 4-92 ）式中，记 [ D] [ K ] 1[ M ] 、 = 1 ， [ D ] 称为动力矩阵，于是（ 4-92 ）式成为求 [ D ] 的特征 2 值问题： [ D] { X } = { X } （4-94 ） 当矩阵阶数较高时， [ K ] 1 不易求，可改用柔度矩阵 [ ] ，其系数 ij 可用单位力法求得，系数 ij 表示仅 在 i 层 （质点）作 用一单 位力时 在第 j 层 （质点 ）产生 的位 移，如 在第 n 层作 用单位 力， 则 1n = 1 , 2 n = 1 + 1 … nn = 1 + 1 +…+ 1 。 k k k k k k 1 1 2 1 2 n 先假定一初始振型向量，如全为 1 的向量 { X } ( 0) ={1 1 … 1} T ，按下式得到一个新的向量： { X } (1) = [ D] { X } ( 0) （4-95 ） 同时将 { X } (1) 标准化： { X } ( 1) = 1 { X } ( 1) （4-96 ） 标准化的方法可规定某个元素为