计算结构动力响应的状态方程直接积分法ξ.pdfVIP

第 17 卷第 1 期 计 算 力 学 学 报 V o l. 17 N o. 1 2000 年 2 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS February 2000 文章编号:(2000) 计算结构动力响应的状态方程直接积分法 1 2 张森文 , 　曹开彬 ( 1 暨南大学应用力学研究所, 广州 510632; 2 中国计算机报, 北京 100038) 摘　要: 利用钟万勰等发展的指数矩阵精细算法, 提出了状态方程直接积分法。它能适 用于确定情形各种激励作用下系统的动力响应分析; 与分段等效线性化方法相结合, 也 可用于某些非线性系统的响应计算。算例表明, 状态方程直接法具有精度高、不受时间 步长的严格限制等特点。 关键词: 振动; 动力响应; 状态方程; 直接积分 中图分类号: O 332　　　文献标识码:A 1　引　　言 对一般动力系统 ¨ M x + Cx + Kx = f ( ) x = c , 　x = c ( ) t , 　 0 1 0 2 1 　　它在时域内的数值求解, 已有很多方法, 如杜哈美积分法, 加权残值法, 中心差分法, 呼伯 [ 1 ] 特法, 纽马克方法, 威尔逊 法 等。但这些方法有个共同的缺点, 即对时间步长很敏感, 且精 度不是很高。 近来, 钟万勰、林家浩等提出了精细时程积分法[ 2, 3 ] , 具有很高的精度。在他们的研究中, 给 A t ( ) 出了指数矩阵e 的精细数值算法, 很有推广价值。对上述动力系统, 钟万勰将激励f t 在每个 [ 2 ] 时间步长内假设为线性的, 从而推导出了一个计算系统响应的迭代公式 ; 而林家浩则将激励 [ 3 ] 展开为简谐激励的和, 从而发展了另一计算系统响应的迭代公式 。另外, 他们所推导的迭代 公式都限于在线性系统中。 本文在上述成果的基础上, 将其作了推广, 发展了状态方程直接积分法( ISEM ) 本法不 D 。 仅精度高, 而且对激励不作基本的假定, 因而可用于各种激励作用下系统的响应分析; 与分段 线性化方法相结合, 也能适用于一些非线性系统。 2　 的导出及其在确定性线性系统中的应用 D ISEM 一般而言, 确定性线性系统的动力方程可以写成 ¨