一个高效率的有限元矩阵摄动法程序.pdfVIP

⑦ 维普资讯 弓一j＆ 第 l7卷 第 2期 太 原 重 型 机 械 学 院 学 报 17 No．2 L996年 6月 JOURNALOFTAIYUAN HEAVY MACHINERY INSTITUTE J％in L996 一 个高效率的有限元矩阵摄动法程序 侯新录 A、 摘 要 本文介绍 了有限元矩阵摄动法的编程原理 ，研完了提高计算速度 的措 关键词 … 力eX~-#t ．堡堕堡堑鎏；，广义特征值；结构动力 l 。J可 中图分类号 0327 0 引言 在各种工程 问题 的分析和求解 中，经常遇到参数小变化 的情形 ，例如结构 的小修改，制 造误差 ，迭代设计 ，灵敏度分析和弱耦 台系统等等 结构参数小变化 引起结构振动特性变化 的问题，对工程结构 的优化设计有重要意义．在工程结构设计 中，为 了取得最优设计方案， 需要多次修改方案 ，而常规作法是每修改一次方案就需要求解一次广义特征值 问题 ，对大型 复杂结构进行迭代振动分析是非常麻烦和昂贵的，因此，需要寻求一种快速的灵敏度分析和 重分析的计算方法 ．矩阵摄动法就是这种快速灵敏度分析和动力重分析 的有力工具 正因为 这样，它引起了国内外 力学工作者 的浓厚兴趣 ，并 已发表 了一系列论文_1 ]．本文将通常 的矩阵摄动法改造成更为有效的有限元形式，用 c语言编写 了一个有 限元矩 阵摄动法程序 ， 可在配备 DOS3．3以上版本的各种微机上运行 本文是这一工作的总结 介绍 了有限元矩阵 摄动法的算法原理和提高程序运行速度的技术措旋，并用算例说明了程序的效率 1 矩阵摄动法的有限元公式 结构振动系统的固有频率和固有振型，是该系统的刚度矩阵和质量矩阵的函数．若结构参 数 (几何尺寸．截面特性，材料特性，质量分布，边界条件等)发生小变化，则刚度矩阵K和质量 矩阵 M 发生小变化，即 K = K0十 ￡Kj． M = Mo+eM 1 (1) 其 中 e是一个小参数．与 e=0对应的系统称为原系统 K0和 M0是原系统 的刚度矩阵和质 收犒 日期 1995— 12— 12 维普资讯 太 原 重 型 机 械 学 院 学 报 996丘 量矩阵 sKl和 sM1代表两者的变化 当 e很小时系统 的特征值和特征 向量都只有小变化，并各有展开式 ： ^J0+ ^l￡+ ，2￡ + ’· (2) 竹 = o+ l￡ 许2￡+… (3) 其 中t。和 许0是原系统 的特征值和特征向量 ，^，l和 竹1是系统特征值和特征向量的一阶摄动 ， ^2和 2是系统特征值和特征 向量 的二阶摄动 在实际计算 中大多只计算到一阶摄动，偶尔计 算二阶摄动，但很少计算更高阶的摄动 ． 将 (2)、(3)代入固有振动的基本方程 (K 一^M) -一O (4) 和归一条件 M吼 ：1 tsl 得 到 ((K0+eK )一 (M0+eM 1))许 =u (6) (Mo+eM 1)竹 ： 1