非线性动力分析避免状态矩阵求逆的精细积分多步法.pdfVIP

振 动 与 冲 击 第27卷第4期 JOURNAL OF VIBRATION AND SHOCK 非线性动力分析避免状态矩阵求逆的精细积分多步法 王海波 ，余志武 ， 陈伯望 (1．中南大学 土木建筑学院，长沙 410075；2．湖南城市学院土木工程学院，湖南益阳 413000) 摘 要：将精细积分法和预估一校正Adams—Bashforth-Mouhon多步法相结合，提出了高精度的精细积分多步法，对 非线性动力状态方程进行求解，避免了对状态矩阵求逆。该方法与精确值和现有Adams多步法进行比较，数值计算结果 表明该方法是一种高精度、高效率和稳定性较好的方法。该方法可方便地进行不同阶次的积分运算。 关键词：非线性动力方程；多步法；预估一校正Adams—Bashforth—Moulton法；精细积分法；预估一校正法 中图分类号：TU311．3，0322 文献标识码：A 用时域直接积分法来计算结构系统的动力响应， ． 厂 ( ，t)变换为 =H- + ，t)+ t)。将精细积分 得到了越来越多学者的注意。钟万勰、林家浩等提出 法和预估一校正Adams—Bashforth—Mouhon多步法相结 的直接精细积分方法 J，用于求解线性定常结构动 合，提出了高精度的精细积分多步法，对非线性动力状 力系统，可以得到在数值上逼近于精确解的数值结果， 态方程进行求解，避免了对状态矩阵求逆。将本文的 为结构动力系统的时程分析开辟了崭新的途径。储德 方法与精确值和现有Adams多步法进行比较，数值计 文等 讨论了精细直接积分法中积分方法的选择问 算结果表明本文方法是一种高精度、高效率和稳定性 题，指出科茨积分、高斯积分是保持精细算法高精度的 较好的方法。该方法可方便地进行不同阶次的积分 较好算法。晋智斌等 利用精细积分方法在线性定常 运算。 系统方面计算精度优势，提出了一种分区域异步长显 式一精细混合积分方法以提高计算效率。 1 避免状态矩阵求逆的多步法 非线性动力学问题的计算研究，一直是解析法占 讨论非线性动力方程 主导地位，而直接的数值积分技术发展相对滞后。尽 = H· + ( ，t) (1) 管解析法本身也有了长足的进步，但仍局限于求解某 些特殊的问题，无法满足工程实际愈来愈多的更为复 将式(1)中的 ( ，t)分解为两项，其中一项为仅与时 杂的多自由度问题求解的需要，难以对非线性动力学 间相关的荷载向量r(t)，即转换为如下非线性动力 行为提供全面的认识。在非线性分析单步法方面，赵 方程 秋玲 j、裘春航等 提出了一种非线性系统线性化的 =H· + ，t)+ t) (2) 精细积分方法，但线性化的计算结果精度不是很高。 利用指数矩阵 可将其转化为如下同解积分方程 张洵安等 提出了一种精度稍高的线性化迭代计算方 tk+】 (f +1)=e砌‘ (f )+f eH(tk+l- ， )d 法。裘春航、吕和祥等 将非线性部分用t 时刻的J． J tk 次多项式来近似，然后借助于分段直接积分法，导出了 — +1