结动振动摄动分析的新方法ξ.pdfVIP

第 15 卷　第 3 期 应 用 力 学 学 报 . 15　 . 3 V o l N o 1998 年 9 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS Sep. 1998 结动振动摄动分析的新方法 崔博文　沈允文　　　　　段晓革　　　　 ( ) ( ) 西北工业大学　西安　710072 　 西北航空公司飞机维修基地　西安　710082 摘　　要 提出了一种用于结构动力修改的设计灵敏度分析的新方法。它发展了N elson [ 1 ] [ 2 ] 陈 等人的方法, 较好地解决了结构修改量大时计算精度低之间的矛盾。数值算例表 明, 本文新方法计算精度高, 易于计算机实施。 关键词: 振动; 摄动分析; 特征值 1　引　言 复杂结构的振动分析通常可用有限元获得其固有频率和模态。由于结构共振要对结构进 行修改时, 传统的做法是对结构每进行一次修改则计算一次固有频率和模态, 直至满意为止。 这就需要惊人的计算机时。因此结构修改后的快速而有效的重分析技术是十分必要的。不仅 如此, 在结构动态设计中还需要有快速而有效的动态灵敏度分析的计算方法。 1976 年,N elson [ 1 ] 首次研究了实对称矩阵特征值特征向量导数计算问题。遗憾的是N el [ 2 ] son 在计算中忽略了特征等式中的二阶摄动项, 影响了其计算精度。陈 提出了矩阵摄动理 论, 通过一阶摄动或二阶摄动计算获得特征值和特征向量改变, 进而获得结构修改信息。但该 文同时指出, 当结构修改量超过 30% 时需要进行更高阶的摄动才可获得满意结果, 这必然会 加大计算工作量。因此有必要提出一种新方法, 以克服为提高计算精度而增大计算工作量的矛 盾。本文提出了一种振动摄动分析的新方法, 并通过算例检验了新方法的精度。 2　理论基础 本文研究的实对称矩阵的特征问题为 (K - M )X = 0 ( 1) i i T ( ) (2) 且满足如下的正交条件 X iM X i = 1　 i = I , …,N T ( ) ( ) X iM X j = 0　 i ≠j , i, j = 1, …,N 3 式中K 、M 分别为