面向黏弹性边界的拱坝地震反应分析方法.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 基于黏弹性边界的拱坝地震反应分析方法 杜修力 赵密 北京工业大学城市与工程安全减灾省部共建教育部重点实验室 北京 摘要 无限地基的能量辐射效应是影响拱坝地震反应的一个重要因素 本文提出将一种新的黏弹性人工边界结合 显式有限元的时域波动求解方法用于拱坝 地基系统的地震反应分析 地震动输入通过分别将底边界入射位移时 程和侧边界自由场时程对应的解析应力时程作为边界荷载实现 该方法具有较高的计算精度 且不存在高 低频 计算失稳问题 文中以小湾拱坝地震反应分析为例 比较了透射边界 黏弹性边界两种计算方法的计算结果 论证 了本文建议方法的有效性 关键词 拱坝 地基系统 时域地震反应 黏弹性边界 透射边界 中图分类号 文献标识码 研究背景 由于拱坝与无限地基间存在不能忽视的能量交换 因此 地震作用下的拱坝 地基系统是一能量开 放系统 对于能量开放系统的动力反应分析 人工边界问题是将难以在计算机上实现的无限模型转换 为便于在计算机上实现的有限模型的关键 各种常用的人工边界处理方法参见张楚汉 等 的 评述 文献 将无限域的远场波动位移解假定为无衰减的平面波位移解与无限域格林函数解的远场 近似位移解的组合 建立了一种新的黏弹性边界 算例表明 该方法对于近场波源引起的波散射问题模 拟具有不低于二阶透射边界 的计算精度 拱坝地震反应分析方法研究中 基于有限元法开发了最早的拱坝地震反应分析程序 其将地基模拟为带有固定边界的无质量弹性体 没有考虑无限地基能量辐射影响 等 和 等 在频域内用有限元 边界元混合法研究了拱坝地震反应 等 也提出了拱坝 地基动力相互作用分析的有限元 边界元 无限元耦合模型 这些频域求解方法采用的是全局人工边 界条件模拟无限地基能量辐射效应 存在不便反映结构和近场地基的非线性特征且计算效率低的问题 为提高计算效率 等 和 等 采用了近似的阻抗曲线来计算地基施加给拱坝上的作用力 这无疑会影响计算精度 实际上 拱坝 地基系统求解规模巨大 且往往包含基岩介质的节理裂隙 断 层和软弱夹层 坝体中缝界面开合 坝体混凝土材料损伤等非线性因素 为了解决这一问题 杜修力 等 将局部人工边界中的一种位移型人工边界 透射边界引入了拱坝地震反应分析 为抑制透射边 界引起的高频失稳现象 等 在文献 的分析模型中引进了与介质应变成正比的阻尼 在一 定条件下实现了数值模型的稳定分析 但这一措施并未在理论上圆满解决透射边界引起的高频失稳问 题 事实上 在使用有限元与透射边界的结合法计算拱坝 地基系统的地震反应时 为获得直观的稳定 数值解 一般需要多次反复试算 这给工程应用带来了不便 也使得可能存在的误判引起计算精度损失 为克服位移型人工边界方法存在的这一固有缺陷 本文提出将文献 建立的应力型人工边界方法 收稿日期 基金项目 国家自然科学基金资助项目 作者简介 杜修力 四川人 教授 博士 主要从事水工结构抗震研究 一种新的黏弹性边界用于拱坝 地基开放系统的地震动力反应分析 并针对位于 度地震设防区 坝高 的小湾拱坝进行地震动力反应分析 并与透射边界方法的分析结果进行对比 有限元运动方程 地震作用下的人工边界运动方程 拱坝 地基系统内结点有限元运动方程见文献 无外源作 用的黏弹性人工边界结点运动方程见文献 本文进一步给出外源地震作用下黏弹性人工边界结点 的运动方程及等效地震荷载的计算方法 人工边界设置为 盒子 型截断形式 地震波以平面波形式自底边界竖直入射 这种边界截断和地震 动入射形式是多数结构 地基地震反应分析通常采用的 与所有的人