一种有限元模型动力缩聚移频迭代法Ξ.pdfVIP

第 16 卷　第 2 期 应 用 力 学 学 报 . 16　 . 2 V o l N o 1999 年 6 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS Jun. 1999 一种有限元模型动力缩聚移频迭代法 瞿祖清　华宏星　傅志方 (上海交通大学　上海　200030) 摘　　要 提出了一种基于矩阵广义逆的有限元模型动力缩聚移频迭代方法。该方法首先 直接从原系统特征方程出发, 导出反映系统主、副自由度之间位移关系的动力缩聚矩 阵的控制方程, 然后给出了相应的迭代求解方法和收敛准则。为了减少求矩阵广义逆 的计算工作量, 本文给出了一种替代方法, 把对一个高阶满阵求逆转化为对一个同阶 高度稀疏矩阵求逆。与已有的动力缩聚迭代法相比, 本文提出的方法具有两个显著的 优点: 其一是迭代收敛速度高; 其二是通过改变移频量可达到提高收敛速度和使得降 阶后模态的动态性能在任意频段上逼近原系统。数值示例表明, 本文提出的方法是切 实可行的。 关键词: 有限元建模; 动力分析; 动力缩聚; 模型降阶; 广义逆 1　引　言 结构动力缩聚在结构模型的降阶、试验一分析模型的相关工作和振动控制等许多领域中 得到了广泛的应用, 因此引起了国内外一些学者的高度重视, 自从 [ 1 ] [ 2 ] 较完整地 Guyan Iron s 提出该方法以来, 已有大量的文献报道了这方面的研究。Suarez 和 Singh [ 3 ] 首先提出动力缩聚 法的迭代求解方法, 但该文是针对标准特征问题而言的, 对于广义特征问题则需要通过对质量 矩阵进行正交分解, 将其转化为标准特征问题, 然后才能进行动力缩聚。显然, 当质量矩阵为一 [ 4 ] 致矩阵时其计算量很大。为此, 作者对该方法进行了推广 , 使其能直接适用于广义特征问题, [ 5 ] [ 6, 7 ] [ 8 ] 并把该方法应用到浮筏减振装置降阶建模中 。 等 把 ’ 于 1989 年提 F risw ell O Callahan 出的 IR S 法进行了改进, 提出了迭代 IR S 法。尽管文[35 ]和文[6, 7 ]从两个不同角度提出了 动力缩聚迭代法, 但可以证明这两类方法在本质上是一致的。 以上各种方法都属于基于原系统的刚度矩阵和质量矩阵的动力缩聚迭代法, 它是近几年 才发展起来的一种有效的动力缩聚方法。与其它一些方法相比, 它具有缩聚精度高, 适用范围 广, 对主自由度的选择要求较低等优点。然而这些方法共同的缺点是: 迭代收敛速度低; 用这些 来稿日期: 修回日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 62 应 用 力 学 学 报 第 16 卷 方法得到的降阶模型只能在低频区逼近原系统的动态特性。 本文在已有方法基础上提出了一种基于矩阵广