学系统工程.PDF

1998年3月 重庆大学学报 (自然科学版) Vo1．31。№ ．2 第21卷第2期 Journal of Chongqing University(Natural Science Edition) Mar．1998 动力学方程Newmark方法的 Neumann级数解及其收敛性。 岁c=) (重庆大学系统工程及应用敦学系，重庆， I第一作者6萝，男，教授)； 庆，4 0 0 0．414 者5 萝，男，教授 岁，岁 重庆大学系统工程萄 摘 要 讨论了莉力学方程的iNknmnark-Neunmnn级数解及其收敛性，给出了级数解的 收敛性条件，以及近似解和误差的估计武。 关键词 ． 塾 堂立堡 Newmrk方法 竺竺磐墼墼；选墼挂 中国图书资料分类法分类号 0302；(Y29 0引 言 ／、／， 见眯坛 在考虑结构系统的动力学问题时，经过离散化可以得到以下形式的二阶动力学方程 [^ { +【c]{壬)+[g3{对={同：})) (1) 这里，[加 为广义质量阵，[网 为广义刚度阵，[c]为广义阻尼阵。式(1)是一个具有常阵的 二阶线性常檄分方程组，由于系数矩阵的阶数较高，甩古典方法(如乘幂法dacobi法)无法 计算求解；甩提取前若干个最小特征值的方法(如子空间迭代或Iznms法)也不太理想，以 及当F(￡)为任意非循环荷载时求解也较困难。因此，实际上目前仍限于按时段逐次迭代推 进，但若时段取大就不精确，取得太小又使计算量惊人增大。文Fz3在Newmark方法的基础 上应用Neuman级数求矩阵的逆，给出求增量的级数解。但是，该文并未对级数解的收敛性及 误差分析进行讨论。众所周知，线性代数方程的计算机求解及级数的计算相对成熟，也较容 易处理，问题就在如何保证计算精确。 对于动力学方程(1)，一般求得的都是近似解，因而对其收敛性或误差分析是很重要的 研究内容。文[1]对Ne~ rk方法的稳定性和精度分析已作了介绍和讨论，因此笔者只对应 用№unmnn级数求得的级数解的收敛性和误差分析进行讨论。 1 级数解及其收敛性 对方程(1)，由Newraark方法，在一个时间步长It，￡+At]内，有 斟 △I一△￡+ 五 一面 ’△￡一 立一at。 (2) 一 面。’△ 一ai 南一面 孟 ·收文日期 1996—12—06 重庆大学学报 (ill然科学暇1 1998盎 其中 西 )一 口At 西。)一l／a(At) 口{ )= 口一1 一 1／ At ”一 ( 口一2) 一壶一1 且当 ≥0．5，a一0．25(O．5+曲。时，N Ⅵn日rk方法是无条件稳定的伽。 考虑t+At时刻的动力学方程 [ { + +EO]{立+△I+[明 {斟 ={F) (3) 将式(2)代入(3)，有 [j {斟 △f一{p) (4) 或写成增量形式 [足]{△耐 △I=IF) + (5) 式中 [j ；[明 +面1’[c]+西 EM]