【精选】精细积分的非线性动力学积分方程及其解法.pdfVIP

Vol. 22 No. 3 　　 　　　　　　 　 　　　　　　　　　 　 第22卷 第3期 固 体 力 学 学 报 2001年 　9 月 ACTA MECHANICA SOLIDA SINICA September 2001 精细积分的非线性动力学积分方程及其解法 吕和祥 　　于洪洁 　　裘春航 (大连理工大学工程力学系 ,大连 ,116023) 摘　要 　给出了非线性动力学积分方程的表达式 ,针对该方程提出了一个显式预测校正 的单步四阶精度的精细积分算法 ,适用于多自由度、强非线性 ,非保守系统. 算例表明该方法精 度高、计算量较少. 关键词 　精细积分 ,数值积分 ,非线性 ,积分方程 1 　引言 非线性动力学问题的计算研究 ,一直是解析法占主导地位 ,而直接的数值积分技术发展 相对滞后 ,直接的数值时程积分能够全面研究非线性动力行为 ,是一种强有力的计算手段. 哈密顿体系能本质地反映许多物理规律 ,提供了一个尝试利用这些规律来研究积分算 [4 ] 法的机会. 专著 将哈密顿混合变量方法引入弹性力学 ,成功地用新的、统一的方法论处理 弹性力学问题 ; 同时还提出了基于2N 类算法的精细积分方法[5 ] ,放弃了通常的差分格式 ,几 乎是计算机上的精确解 ,能够获得优于传统差分算法的数值结果[2～6 ] . 本文从哈密顿系统入手 ,继承精细积分方法高精度等特点 ,给出哈密顿体系下 ,非线性 动力学的积分方程表达式 ,在积分时采用了当前要计算的时间域内的未知函数的插值 ,因而 避免了如亚当姆斯 、哈明法等采用的外插的多步法 ,提出了一个显示单步预测校正四阶精 度的精细积分有效算法 ,能够自起步 ,适用于多自由度 、强非线性非保守动力学系统. 通过算 例说明 ,本算法精度高 、计算量小. 2 　结构的瞬态动力分析[ 2] 结构的动力方程 ( ) ( ) α ( ) β ( ) ( ) MX¨+ GX+ KX = r r , 　　X 0 = , 　　X = 　　 n 阶 1 0 ( ) ( ) 其中 X 是 n 维向量 , M 是质量矩阵 正定 、对称 , G 是阻尼矩阵 对称 、非负 和陀螺矩阵 ( ) ( ) ( ) α β 反对称 之和 , K 是刚度矩阵 对称 、非负 , r t 是随时间变化的载荷向量 , 和 是已知的 初始条件. 精细积分法是建立在下面的一阶常微分方程组 v= Hv + f , 　　v (0) = v0 (2) 的基础上 ,其中 v 是 2 n 维向量 , H 是由已知矩阵 M 、G 和 K 计算出来的 2 n ×2 n 的矩阵 , f 是由 r ( t) 构成的已知列向量. ( ) ( ) 为了把方程 1 化成方程 2 的形式 ,可借助哈密顿变换来实现. 令