基于有限元分析的特征值反问题求解的逆摄动方法.pdfVIP

第 卷第 期 计 算 力 学 学 报 18 1 Vol.18 No.1 年 月 2001 2 C~INESE JOURNAL OF COMPUTATIONAL MEC~ANICS F b a 2001 文章编号:1007-4708(2001D01-0048-07 基于有限元分析的特征值反问题 求解的逆摄动方法 苗玉彬 滕弘飞 丁金华 屈福政, , , (大连理工大学 机械工程学院 大连, 116024; 中科院现代制造CAD CAM/ 技术开放实验室 沈阳, 110003D 摘 要 本文研究特征值反问题的求解方法 根据广义特征值反问题理论和有限元法的特点 以转子: G , 系统平面梁单元有限元模型结构分析的特征值反问题求解为例,给出一种新的逆摄动方法 给出了G 本逆摄动法较完整的理论基础;给出了其逆摄动参数的显式计算公式及相应的取值方法 本逆摄动G 法也可推广到其他单元类型的有限元模型特征值反问题的求解G 关键词 特征值 反问题 摄动 算法 逆摄动: ; ; ; / 中图分类号: 122 文献标识码: T~ A 1 前 言 在机械或工程结构的设计中 经常要求结构满足指定的自然频率或避开指定的频率禁区 并且往往还要, , 求 结构的重量最小或变化最小G结构的自然频率或特征值通常是利用刚度阵和质量阵 它们是诸如单元 截面积 单元长度等一系列参数的函数 计算得到的 特征值反问题就是要确定这些参数来获得指定的特G 征值G 下面给出特征值反问题的数值定义 设 由下式求得 G A(xD : n A(xD = A + xA (1D 0 i i z = 1 其中 n 是 阶实对称阵 另设实数 为 的特征值 ,x R ,{A } n G (xD (xD A(xD G i 1 n 特征值的求解是根据广义特征方程 K/ = M/ , = 1,2, ,mj j j (2D 其中 为刚度阵 为质量阵 为特征向量 为特征值 下标 为特征值的阶数 :K