有限元法基础-12动力学问题[文字可编辑].ppt

12. 动力学问题 有限元法基础 31 12. 动力学问题 有限元法基础 32 12. 动力学问题 12.4 特征值问题及其解法 系统的运动方程为 无阻尼自由振动退化为 设方程解的形式为 方程成为 有限元法基础 33 ( ) ( ) ( ) ( ) t t t t ? ? ? & & & Mq Cq Kq Q ( ) ( ) 0 t t ? ? & & Mq Kq ? ( ) i t t e ? ? q q 2 ? ( ) 0 i t e ? ? ? ? ? M K q 2 ? ( ) 0 ? ? ? K M q 广义特征值问题 2 ? ? ? ? 特征根 特征向量 q 12. 动力学问题 四种类型的解法： ? 直接矢量迭代法（幂法） ? 矩阵变换法 ? 多项式迭代求解法（行列式搜索法） ? 利用特征多项式的 Sturm 序列特性求解法 以及 i i i ? K φ = M φ ( ) I ( ) T T ? φ K φ = φ M φ = 对角阵 单位阵 ( ) det( ) 0 p ? ? ? ? ? K M ( ) det( ) p ? ? ? ? K M ( ) ( ) ( ) ( ) ( ) ( ) det( ), 1, 2, , 1 r r r r r i p r n ? ? ? ? ? ? K M L 12. 动力学问题 12.4.1 逆迭代法（幂法） 对方程 取近似解 按以下迭代格式求解 则序列 将收敛于相应的特征根 的特征 矢量。 ? K φ = M φ 1 x 1 ( 1, 2,3 ) k k k ? ? ? Kx Mx L 1 2 , , x x L 1 ? 12. 动力学问题 ? 因为对任一矢量可用特征矢量表示为 代入方程 按迭代方程有 若 ，当 时， 1 1 0 ? x M φ k ? ? 1 1 k ? ? x φ 1 1 { } { } n i i i x ? ? ? ? ? 1 1 1 1 { } { } { } n n i i i i i i i x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? M M K 1 2 1 1 1 { } { } { } n i i i i x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? K M 1 1 1 1 1 1 1 { } { } { } k k k n n k i i i i i i i i x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 2 n ? ? ? ? ? ? L 设 12. 动力学问题 ? 为了使 X i 不受计算的影响，常常需要归一化 ? 正迭代法的计算方案 迭代格式 若 ，当 时， ? 特征根 的近似解 1 1 1 1/2 1 ( ) k k k T k ? ? ? ? ? x x x M x 1 1 1 1 1/2 1 ( ) k k k k k T k ? ? ? ? ? ? ? Mx Kx x x x M x 1 0 T n ? x M φ k ? ? 1 k n ? ? x φ 1 1 1 1 1 ( ) k k k T k k ? ? ? ? ? ? ? x Kx x x M x 12. 动力学问题 12.4.2 变换法 ? 广义特征值问题化为标准特征值问题 有限元法中的质量矩阵 M 是对称正定的，则 故有 定义 得到 T ? M SS T ? ? ? ? K φ M φ SS φ ? 1 T ? ? ? ? S KS φ φ T ? φ S φ ? ? K φ φ 12. 动力学问题 ? 标准特征值问题 变换法中有 Jacobi 法、 Givens 法、 Householder ，其实质就 是通过一系列的变换矩阵，将 M 变换成单位矩阵，将 K 变换成对 角矩阵。 ? Jacobi 法 标准特征值问题的方程 设完成第 k 步变换成为 P k 是正交矩阵 ，即 ? ? K φ φ 1 T k k k k ? ? K P K P T k k ? P P I 12. 动力学问题 ? P k 矩阵的构造 12. 动力学问题 ? 特点 ? 在 时，矩阵 K 趋于对角阵 ? 由于只能做有限次变换，因此最后的矩阵是对角占 优 ? 变换后的矩阵总是对称的，可以减少计算次数 ? 在一次变换使非对角线为零元素，在下次变换中可 能成为非零，因此收敛缓慢 ? 需要结合一些其他策略提高计算效率 k ? ? 12. 动力学问题 12.4.3 子空间迭代法 子空间迭代法是求解大型特征值问题的低阶特征值有效方 法，它实际上是 Rayleigh-Ritz 法和同时逆迭代法的组合 。 ? 子空间迭