Wilsonθ法直接积分的运动约束和计算扰动.pdfVIP

第22卷第4期 　 计 算 力 学 学 报 　 V o l. 2 2 , N o . 4 2005 年 8 月 Ch in ese Journa l of Com puta t iona l M echan ics A u gu st 2005 文章编号:(2005) W ilson - 法直接积分的运动约束和计算扰动 1 1 2 黄庆丰 , 　王全凤 , 　胡云昌 ( 1. 华侨大学 土木系, 福建 泉州 3620 11; 2. 天津大学 建工学院, 天津 300072) 摘　要:W ilson 法的积分过程一般不可能同时既符合计算假设规定的运动约束条件又满足动力平 衡方程, 时间步长内附加了一个计算扰动影响。由W ilson 法积分计算出的时间步长终点不平衡加 速度和系统的动力平衡方程, 本文导出了时间步长内计算扰动的确定方法, 并进一步采用① 同步计 算消除计算扰动效应和② 后续步计算消除计算扰动效应, 两种途径抵消其不利影响。算例指出, 本文 方法有效减少W ilson 法直接积分结果的误差和超越现象, 提高了计算稳定性。 关键词:W ilson 法; 运动约束; 动力平衡; 计算扰动; 抵消 中图分类号: TU 3 11. 3　　　文献标识码: A 动力平衡方程 1　引　言 M {x } + C {x } + K {x } = P ( ) ( ) t 1 动力方程求解的直接积分方法大致分为差分 式 中 M , C 和 K 分别为质量、阻尼和刚度矩阵, ( ) 型和积分型两类。差分型如 法、 P t 为动力荷载矩阵。给定初始位移{x 0 }、速度 W ilson N ew m ark 法和H oubo lt 法等存在的问题是仅能考虑系 {x 0 } 和动力P 0 , 求出初始加速度,W ilson 法用递 统在离散时刻的动力平衡, 而对于计算时间步长内 推算法求出其他离散时刻的系统运动参数。 各点的动力平衡则无法顾及, 容易引起计算误差和 第 i 个 = t 时间段的位移增量 ～ ～ [ 1 ] [2 ] {x , i } = K - 1 P , i (2) 计算不稳定 ; 积分型如精细积分法 、积分型直 , i 接积分法、样条积分法[3, 4 ] 等, 很好克服了这一问 ～ ～ 其中