频响函数计算高精度级数展开法.pdfVIP

第 15 卷第 2 期 计 算 力 学 学 报 V o l. 15 N o. 2 1998 年 5 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS M ay 1998 频响函数计算的高精度级数展开法 瞿祖清　　傅志方 (上海交通大学振动、冲击与噪声国家重点实验室, 上海, 200030) 摘　要　把作者近期提出的用于频响函数计算的级数展开法进行了推广, 使其计算精度进 一步提高, 并且能适用于自由自由系统。数值示例表明, 本文提出的方法是有效的。 关键词　结构动力分析; 频响函数; 模态截断; 模态加速级数展开法 分类号　 325 O 1　引　　言 在振动控制、结构动力分析、参数识别和故障诊断等领域中频响函数的计算是一个重要环 节, 往往也是工作量所在。如何快速精确地获得频响函数, 在工程应用中显得越来越重要。如 果用线性方程组理论直接求解, 则会因系数矩阵为频率的函数而使计算量极其庞大。于是经典 模态展开法应运而生, 当只需计算部分自由度上的频响函数时, 该方法的优点尤其突出。然而, 在非完备模态集的条件下, 该方法给出的频响函数截断误差较大。文[1 ]提出的改进模态法和 [ 2 ] 作者近期提出的级数展开法 能使截断误差大大减小。本文在此基础上, 导出了一种高精度级 数展开法, 该方法的计算工作量与级数展开法相当, 但计算精度可进一步提高, 并且它可适用 于自由自由结构。 2　理论基础 ( ) 设 n 自由度无阻尼系统 其结果可很方便地推广到有阻尼情况 在傅氏域中的运动方程 为[ 3 ] ( 2 ) ( ) ( ) ( ) K - M X = F 1 其中　K 和M 分别为该系统的刚度矩阵和质量矩阵; X () 和F () 分别为位移列向量和载荷 列向量; 为激励频率。由上式可得 X ( ) ( ) ( ) ( ) = F 2 式中H () 即为频响函数矩阵, 定义为 ( ) ( 2 ) - 1 ( ) H = - M 3 ( ) ( ) 2 ( ) 在式 2 和 3 中假定 ≠ i = 1, 2, …, n , 为系统的特征值, 当该条件不满足时可用 i i M oo reP en ro se 广义逆求解。 假定该系统的正则振型为 = [ … ],