无限地基辐射阻尼对重力坝地震响应的影响.pdf

无限地基辐射阻尼对重力坝地震响应的影响 李志全，杜成斌，洪永文 (河海大学土木工程学院工程力学系，江苏 南京 210098) 摘要：大坝地震响应计算是混凝土大坝抗震安全评估的关键，其中无限地基辐射阻尼对坝体地震响应的 影响是不可忽视的。本文在大型商业软件ABAQUS二次开发的基础上，对结构与地基的动力相互作用进 行研究，采用有质量粘弹性边界来处理地基中辐射阻尼，并与无质量固定边界进行比较。结合在建中的 金安桥混凝土重力坝，对其厂房坝段的地震响应进行了数值分析，从位移、应力和动能等角度对两种边 界的响应进行了比较，得到了一些有价值的结论。 关键词：辐射阻尼；重力坝；地震响应；坝- 基岩相互作用；人工边界；ABAQUS 1 引 言 二十世纪七十年代以来，与土－结构动力相互作用问题相关的动力基础振动、地 震波散射等近场波动问题的数值模拟方法得到了很大的发展[1]。在实际工程中，大体积 建筑物和地基的相互作用，不仅影响结构的动力特性，而且，结构振动的反馈又影响 地震动输入。更为重要的是结构振动的能量将沿地基向远处逸散，形成所谓的辐射阻 尼效应[2]。无限地基辐射阻尼对结构动力相互作用有重要影响已达成共识，但不同的研 究成果有相当差异，目前仍处于探讨阶段，尚难得出公认的方法和结论。通常有粘性 边界、无穷边界元和人工透射边界等处理方法。粘性边界概念清楚、简单方便，但粘 性边界存在精度不高，且低频失稳等问题。无穷边界元采用无穷元（或半无穷元）， 用无限地基模拟远场无限地基的辐射阻尼效应，在空间或时间上是非局部的，在频域 中确定地基动力刚度时要求解巨型方程组，且对地基介质非均匀性的适应差。廖振鹏 院士提出的人工透射边界的处理方法在波动分析中得到了广泛的应用，但存在高频振 荡失稳和漂移失稳问题。可见，目前工程界还没有形成一个统一有效的关于边界处理 方法。 本文在大型商业软件ABAQUS二次开发的基础上，对结构和地基的动力相互作用 进行研究，研究采用粘弹性边界条件探讨无限地基辐射阻尼对重力坝地震响应的影 响，结合在建的金安桥重力坝进行分析研究，以期能得到一些启发和认识。 1 2 无限地基数值模型的建立及地震波的输入 在考虑重力坝地基结构动力问题时，应注意到天然地基是无限延伸的，地震能量 将沿地基向无穷远处逸散，形成所谓的辐射阻尼效应。在数值计算中如何正确有效地 模拟无限地基的能量辐射作用，是决定能否更精确的模拟结构地震响应的关键技术之 一 [3] 。本文采用二维平面有限元模型，如图1。由于受计算规模的限制，截取的地基只 能为有限范围，利用在边界设置弹簧和阻尼器单元的方法来模拟无限地基的辐射阻尼 效应。采用ABAQUS软件对重力坝厂房进行建模，单元采用四节点等参元，单元数和 节点数分别为18625和19021。 图1 重力坝厂房坝段有限元模型 在地基－结构动力相互作用有限元分析中，采用地基底部均匀输入水平方向和竖 直方向地震波的方法来实现对地震的模拟。本文采用的地震为依据《水工建筑物抗震 设计规范》（DL5073-2000 ）关于重力坝抗震设计提供的有关反应谱参数生成的地震波 [4] 。 表1 坝址设计地震波峰值加速度及反应谱参数的选取 Amax(g) T (sec) T (sec) β γ 备注 1 2 max 0.400 0.1 0.2 2.0 0.9 采用三角级数法生成人工地震波，对于给定的功率谱密度函数S (ω) ，按式（1）