拱坝地震输入模型中黏弹性边界与黏性边界比较.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 拱坝地震输入模型中黏弹性边界与黏性边界的比较 刘云贺 张伯艳 陈厚群 西安理工大学水利水电学院 陕西西安 中国水利水电科学研究院工程抗震研究中心 北京 摘要 本文首先介绍了动力分析中无限地基数值模拟采用的各种方法 包括无穷元方法 边界元方法 无穷边界元 方法 无限边界处理有叠加边界 黏性边界 黏弹性边界和人工透射边界 介绍了拱坝地震自由场输入模型 通 过算例和实际工程计算 验证了黏性边界和黏弹性边界均有较好的吸能效果 且两种边界条件下 所得计算结果几 乎一致 黏性边界或黏弹性边界均可作为拱坝地震自由场输入模型中的吸能边界 关键词 拱坝 黏性边界 黏弹性边界 自由场 中图分类号 文献标识码 研究背景 在静力分析中 地基对结构的作用可通过取足够大的地基范围来模拟其影响 但在动力分析中 由 于地震波向无限远域传播 因此 增大地基范围也不能有效地模拟地基与结构的动力相互作用 必须采 用人工边界吸收截断边界上的外传波 对无限地基的数值模拟是吸引众多研究者的一个热点问题 主要方法有无穷元方法 边界元方法 无穷边界元方法 在无限边界处理上有叠加边界 黏性 边界 黏弹性边界 和人工透射边界 无穷元是外部区域向无限远处延伸的辐射状窄条 在形函 数中包含反映位移衰减特征的衰减函数 衰减函数与所模拟的波的波数有关 常在频域分析中使用 边 界元在近几十年来发展较迅速 其优点是能自动满足无限远处的幅射条件 只需对边界积分方程进行数 值积分 精度较高 但形成的刚度阵往往是满阵 叠加边界的基本思想是将透射边界分解为两个已知的 非透射边界 分别求出两个边界问题的解 再通过叠加获得满足透射边界条件的解 其缺点是难以处理 多次反射问题 上述方法可分为两类 一类与频率求解相适应 如无穷元 边界元等 另一类与时域求解 相对应 如黏性边界 人工透射边界等 它们的共同特点是模拟地震波向远域的传播 对于工程实际应 用而言 边界条件应满足两点要求 一是有一阶以上的精度 二是程序容易实现 有较高的计算效率 由于时域显式有限元法在动力分析特别是非线性动力分析上的优越性 大型结构的抗震分析大多 使用时域显式有限元方法 二次透射人工边界具有二阶精度 但透射边界的编程较复杂 稳定性较差 且要外加很多节点 使求解规模过大 比较费机时 基于此 本文重点关注黏性边界和黏弹性边界的应 用 黏性边界具有一阶精度 概念清楚 实现起来简单 其缺点是仅考虑了对散射波能量的吸收 不能模 拟半无限地基的弹性恢复能力 文献 认为 黏弹性边界在计算精度和稳定性两方面都比黏性边界优 越 但其只针对简单的算例进行了计算 其结论是否适合实际工程问题有待进一步验证 本文将对黏性 和黏弹性边界条件作深入的分析和比较 为拱坝地基系统有限元分析中合理选用吸能边界条件提供依 据 收稿日期 基金项目 国家自然科学基金资助项目 作者简介 刘云贺 男 辽宁人 博士 副教授 从事结构抗震研究 拱坝地震自由场输入模型中的黏滞性边界问题 对于计入无限地基影响的拱坝地震反应分析 已建立了许多种计算模型 其中拱坝地震自由场输入 模型 完全能满足工程实际应用的需要 图 是基于黏弹性边界条件的地震自由场输入模型 其一端与内域有限元相连 另一端连接远域地 基 很显然 当黏弹性边界完全吸收有限元区域的外传散射波时 与远域地基相连的一端承受的是地震 自由场运动 这样 地震的输入问题转化为在黏弹性边界的弹簧 阻尼器上作用自由场运动的问题 对 于层状地基的自由场计算已有较成熟的计算方法 与远域地基相连的黏弹性边界承受 的是地震自由场运动 可转化为作用在结构地基系统上的边界力 图 坝的地震自由场输入模型 设黏弹性边界外侧的边界节点 处自由场位移向量为 自由场速度向量为　 自由场应力张量 为 构成黏弹性边界的弹簧系数为 阻尼系数为 则作用在弹簧阻尼器上的力矢量为