第四章实用计算方法精要.pptVIP

* 例：设 求X0 * 例：试用子空间迭代法求图示刚架的前两阶频率和振型。 解：由前例可知： 1）、选取X0 * 2）、求 由： 3）、求 ( ) 4）、求缩减后的矩阵 此矩阵方程的特征方程是： 将 代入，解此方程，得 求特征向量A1 将 回代，令a11=1,得： * a21=1.0000114 同理， 令a12=1,得：a22=-2.5 故得： 5）、求振型X1 按第一个元素归一化，得： 本例精确解为： 一次迭代精度已接近精确解 * 练习：试用子空间迭代法图示三个自由度的体系前两阶频率和振型。 EI=常数，m1=m2=m3=m，a=L/4 m1 m2 m3 a a a a EI * §4.2 矩阵特征值问题及解法 则其下三角 i=1,2, …,n行的元素为: * §4.2 矩阵特征值问题及解法 2、标准特征值问题解法 特征值问题: 一是求解它的全部特征值问题，即所有的特征值和对应的特征向量； 另一是求解它的部分特征值问题，即部分（通常是最小或最大的一部分）特征值和对应的特征向量。 * 在结构动力学中，往往矩阵的阶数都很高，有时不可能，而且也没必要求解全部特征值和特征向量。 求解方法分为两大类: 一类是变换法,另一类是向量迭代法。 雅可比法 幂法 子空间迭代法 §4.2 矩阵特征值问题及解法 1）雅可比法 雅可比法是求解实对称矩阵全部特征值和特征向量的简单有效方法。 * [A]为n×n阶的实对称矩阵。 根据线性代数理论，任何一个 n×n实对称矩阵[A] ，可通过一个的正交矩阵[S]，经相似变换化为一对角矩阵，即 矩阵 [D]的n个对角元素就是[A] 的n 个特征值，而[S] 的第i 列，就是[D]中第i个元素所对应的特征向量。 §4.2 矩阵特征值问题及解法 正交矩阵[S]实际上是一个坐标变换矩阵 * 功能: 实现坐标系的旋转，使旋转后的n个坐标轴的方向与矩阵[A]的n个互相正交的特征向量的方向一致。 雅可比法的基本思想就是通过多次坐标系统的旋转来逐渐实现所想达到的最终旋转 ，而其中每一次分步旋转则是通过正交矩阵 所实现。 正交矩阵[S]不容易寻找。 §4.2 矩阵特征值问题及解法 * 幂迭代法是一种渐进解法，可用来计算频率及相应的振型，且随着迭代轮次的增加，可以得到任意精度的解。 一、求基本频率和基本振型 2）幂迭代法 设A为多自由度体系的振型列向量，则由频率方程： 令 H= -----称为动力矩阵，代入上式，得： （1） §4.2 矩阵特征值问题及解法 * （1） 若已知A，则 即可求出，要确定振型A，可先取任一振型，如取 代人（1）右边，即可求得新的向量A1，然后将A1按首项归一化，得到新的 ， 再将 代人（1）右边，又可得到A2，再归一化得到 如此继续下去得到A3、 …，An、 当 时，就停止迭代，此时可把 作为基本振型，将其代人（1）即可求得 * 例：图示为三个自由度的体系，求其基本频率和振型。 EI=常数，m1=m2=m3=m，a=L/4 m1 m2 m3 a a a a EI 解：采用图乘法可知 设 代人（1） * 再把 代人（1） 再把 代人（1）得 第四次迭代 此时 故取 作为基本振型，代人（1）得： 据此有： * 二、迭代法求最高频率和振型 如果把迭代法应用于刚度法建立的动平衡方程，则可得到最高频率和振型。 由 两式左乘 得： 而 令 代人（a）得： 与（1）式比较，当采用同样的迭代法计算时，显然得到 是最小的，因而 是最高频率，对应振型也为最高。 * 例：试用迭代法求图示刚架的最高频率和振型。 解：可求得： 于是 * 由迭代公式 取 至此迭代结果已收敛，相应的频率： * 三、求较高阶频率和振型 前述方法一求出的是第一振型，为什么？ 这是因为第一振型分量的影响随着迭代次数的增加而愈来愈大，故最后迭代结果将收敛于第一振型，说明如下： 由式（1） 即： 设以任一初始向量 去迭代，该向量实际上包含各振型分量，可用各振型分量的线性组合表示： 故迭代k次后为： 当k很大时， （2） * 如选择初始向量时，能设法消除第一振型的影响，使 那么迭代结果自然就会收敛于第二振型