结构动力方程的更新精细积分方法.pdfVIP

第36卷第2期 力 学 学 报 V01．36．No．2 2004年3月 ACTAMECHANICAS1NICA Mar，．2004 结构动力方程的更新精细积分方法 。湖南大牡差塑，长沙410082，。中国建筑科学研究盆嚣。京100013，俐弓J月(湖南大学土木工程学院，长沙 ) (中国建筑科学研究院抗震所，北京 )～。， 。J 摘要将高斯秘分方法与精细积分方法中的指数矩阵运算技巧结合起来，建立了精细积分法的更新形式及计算 挝程，对该更新精细积分方法的稳定性进行了论汪与探讨．在实施精细积分过程中不必进行矩阵求逆，整个积 分方法的精度取决于所选高斯积分点的数量．这种方法理论上可实现任意高精度，计算效率较高，其稳定性条 件极易满足数值例题也显示了这种方法的有效性． 关键词结构动力学，时程分析，指数矩阵运算，高斯积分，稳定性分析 最近，Zhoag[7J提出了1种具有任意高精度的精 细积分方法(PTSIM)，这种方法对线性定常结构动 自从20世纪70年代初以来，结构动力反应计力系统，能够得到在数值上逼近于精确解的结果，但 算的数值方法已经发展．这些数值方法的建立主要 是对于非齐次动力方程涉及到矩阵求逆的困难，计 基于两个假设，即： (1)不是刻意要求在任何时刻 算精度取决于非齐次项的拟合精度等问题，这些问 满足结构动力方程，而是在离散的时间间隔点上满 题的存在极大地限制了精细积分方法的广泛应用． 足结构动力方程；(2)在时间间隔At内，加速度、 本文将精细积分方法的基本原理与高斯积分方 速度、位移反应按某种假定的规律变化，直接积分方 法结合起来，建立了1个新的精细积分格式，即更 法，如差分方法(中心差分法，Houbolt法)，线性加新精细积分方法(RPTSIM)．新的积分格式只需进行 速度法及改进的线性加速度法(Wilson0，Newmark指数矩阵运算，避免了矩阵求逆问题，无需对非齐 次项进行数学拟台，整个积分格式的计算精度取决 ／3，Hilber法)【lJ已广泛应用于工程实践．然而这些方 法存在明显的弊端．首先，这些方法一般只具有一 于高斯积分点的数量，从理论上说，该算法可达任 阶精度或二阶精度，精度不高，尤其在高频阶段更 意高精度，通过算例，比较了本文的精细积分方法 是精度较低；另外，这些方法均具有算法阻尼，这种 与传统精细积分方法的精度及计算效率． 算法阻尼难以有效地控制，从而使计算结果产生较 1更新精细积分方法的基本原理与公式 大误差． 离散结构模型的动力方程为 为r克服上述方法的不足，一些高精度的数值 积分方法已被提出． Golley[2j将时间间隔△t内的 M爻+c文+KX=F(t、 高斯点作为配点，得到了1种三阶精度的有条件稳 定格式． Rift和Batch[3J根据哈密尔顿变分原理， X(O)=Xo，x(o)=Xo 应用3次内插函数，构造了1个暇阶精度的有条件 式中，M，c，K分别表示质量矩阵、阻尼矩阵和刚 稳定格式．Argyris[4J等应用Hermitian内插函数， 由配点法，得到了1种四阶精度的无条件稳定的积 度矩阵，F(t)为给定的外力列向量．作变换n令 分格式． Kujawski和Gallagher[5J应用广义最小二 f