无限域波动问题的局部非协调网格法.pdf

No．1 第19卷第1期 振 动 工 程 学 报 V01．19 2006年3月 JournalofVibration Mar．2006 Engineering 无限域波动问题的局部非协调网格法 贺向丽，李同春 (河海大学水利水电工程学院，江苏南京210098) 摘要：基于远置人工边界的原理采用了局部非协调网格法，即以协调点集的位移、速度、加速度变量作为基本变量· 非协调点集的变量作为从变量，通过线性插值建立从变量与基本变量之间的线性关系，从而导出以基本变量作为 未知量的总体线性方程组。并通过算例验证了本文方法的精度。本文方法不仅大大减少了计算量，提高了求解效 率，而且避开了大角度入射和人工边界稳定性的问题，还可以用于其他人工边界的检验工作。 关键词：无限域；波动；远置人工边界；局部非协调网格；插值 中图分类号：TV312文献标识码：A 文章编号：1004—4523(2006)01—0124-04 构的影响，而对无限域地基本身动力响应的精度没 引 言 有过高要求。 针对以上问题，本文以人工边界早期的处理方 广’r 无限域地基的波动问题是高坝抗震的关键问 法为基础，即取人工边界至结构的距离L≥～-6-，并 题。当用有限离散模型模拟无限地基时，将在人工截 将李同春等提出的局部非协调网格法的思想应用于 取的边界上发生波的反射，从而引起波的振荡，导致 无限域波动问题口州，对结构采用较密的网格。随着 模拟失真。解决这个问题最有效的办法就是引人人 计算区至结构的距离增大，网格逐渐变粗，从而大大 工边界条件。为此诸多学者提出了多种人工边界以 减少了计算量，提高了求解效率。并通过算例进行了 解决有限截取模型边界上波的反射问题[1]。 比较，验证了本文方法的精度。 人工边界早期的处理方法是取人工边界至结构 ^7r 的距离L足够大[2]，即取L≥等，T为计算时间，c 厶 1局部非协调网格法 为无限域介质的最大波速，则结构反应在o≤f≤T 的时段内将不受边界条件的影响，边界条件可以任 1．1 网格模型 意给定。虽然这种远置人工边界的方法计算量大，效 率低，尤其对于大型的三维结构将对计算机的内存 常规的有限元分析一般采用全协调网格，即研 提出更高的要求，但是它至今仍在各种建议的人工 究域内的任意一节点相连单元的集合包含以该节点 边界条件的检验工作中，用于提供精确地消除了边 为中心的微小空间全区域，如图1。 界影响的数值解。为了减少计算量，提高求解效率， 局部非协调网格的特点是局部区域内与某节点 近30年来，国内外的学者陆续提出多种人工边界条 相连单元的集合不能覆盖以该节点为中心的微小空 件。其中比较有代表性的有粘性边界、旁轴边界、叠 间全区域，如图2。其