结构动力分析中Kunar方法的改进和收敛性.pdfVIP

⑤ ， 维普资讯 口 第 l7卷 第 3期 太 原 重 型 机 械 学 院 学 报 Vol 17 No．3 1996年 9月 JOURNALOFTAIYUAN HEAVYMACHINERY INSTITUTE Sep 1996 结构动力分析中 Kuhar方法的改进和收敛性 侯新录 (太原 I业大 学，太原 020024) 摘 要 本文介绍 了结构动力分析 中基于 Kuhar方法 的动力缩减迭代法 ，应 用矩 阵摄动理论证明了其收敛性 ，算例说 明了方法的有效性 葶中键图词分类1．竺号兰0童 ；，缩减矩 阵摄动法 327 ／／ 0 引言 在大型复杂结构动力分析有 限元模型 中，其 自由度数量少则 几十几百 、多则 几万几十 万 从计算机费用和计算机 内存考虑 ，经常希望减少求解 的规模 为此 目的，Guvan一“和 Irons 提 出了缩减方 法 ，这个 方法 至今 仍被广 泛应 用 于大型结构动力分析程 序 ，如 NAS． TRAN Guyan缩减采用了仅考虑静力效应 的位移变换 ，略去了被消去的副 自由度 的惯 陛效 果 由于 主 副 自由度仅 存 在 静 力耦 合 ，因此 ，在 高频 段 内可 能产 生 大 的 误 差 为此 ． Kuhar 等人作了改进 ，提 出动力变换 的方法 ．这是 一个很好 的思路 ，但作 者未对 方法 的收 敛性 、参数选取等问题作理论探讨 ，此外 ，作者 只限于讨论 由模 态综 台法导出的一类特征值 问题 ，其副 自由度 自然地选为远离关心频段 的模态座标 ．对于有 限元 法导出的大型广义特征 值 问题如何执行 该法 ，作者未予考虑 ．本 文给 出了 Kuhar方 法在一般有 限元动力分析 中的 迭代格式 ，并以矩阵摄动法为工具对这一方 法作 了理论探讨 ． 1 两种常用的缩减降价方法 1．1 Guyan方法 对结构作有限元动力分析 ，其广义特征值 问题 为 Kx = 6o2MX 收藕 日期 1996—03—25 维普资讯 第 17卷第 3期 候新录 ：结构动力分析 中Kuhar方法的改进和收敛性 255 · · (2) 这里 X 是 副 自由度 ．在考 虑它的位 移时，认为是按准静态的方式 (即 w 2≈ 0)，由主 自由度 X1 (3) (4) (5) 将 (4)代人 (1)，两端左乘 丁，就得 降阶 的广义特征值 问题 (K 一ⅢM )X．=0 (6) 其中 K = TKT，M = 丁7~MT (7) 在将 K，M 分别缩减 为K ，M 的过程 中，并未将副 自由度 的动能忽 略掉 ，而是假定 了一种静