HyperAI超神经

考虑动水压力的 Koyna 地震非线性动力响应分析

教程简介

Abaqus 简介

Abaqus 是一款功能强大的有限元分析 (FEA) 软件,广泛应用于工程模拟领域。它通过有限元方法对各种工程问题进行模拟和分析,能够处理从简单的线性问题到复杂的非线性问题。 Abaqus 最初于 1978 年发布,由 Hibbitt, Karlsson & Sorensen, Inc.(HKS) 开发,后更名为 ABAQUS 公司,并于 2005 年被达索系统 (Dassault Systèmes) 收购。

教程简介

该教程为 Abaqus 官方教程:混凝土重力坝的抗震分析。本例说明了混凝土损伤塑性材料模型在评估任意载荷作用下混凝土结构稳定性和损伤的典型应用。

教程将对 Koyna 大坝进行分析,该大坝于 1967 年 12 月 11 日遭受了里氏 6.5 级地震。之所以选择该问题,是因为它已被许多研究人员广泛分析过,包括 Chopra 和 Chakrabarti (1973) 、 Bhattacharjee 和 Léger (1993) 、 Ghrib 和 Tinawi (1995) 、 Cervera 等人 (1996) 以及 Lee 和 Fenves (1998) 。

基本参数介绍

koyna 大坝典型非溢流坝段的剖面几何形状如图 1 所示。该坝段高 103 米,底部宽 71 米。假设坝段的上游面直立。地震发生时水库的水深为 hw=91.75m 。根据研究人员的做法,考虑在平面应力条件下对非溢流坝段进行二维分析。

图 1 koyna 大坝典型非溢流坝段剖面

用于分析的有限元网格如图 2 所示。该网格由 760 个一阶、减缩积分、平面应力单元 (CPS4R) 组成。节点的定义以大坝左下角为中心的全局直角坐标系,其中垂直的 y 轴指向上方,水平的 x 轴指向下游方向。

图 2 有限元分析网格

在这个例子中,我们通过假设基础是刚性的来忽略坝体与基础之间的相互作用。坝体与水库之间的动态相互作用,由于地面运动的横向分量,可以通过 Westergaard 附加质量技术以简单形式进行建模。根据 Westergaard(1933 年)的理论,地震期间水对坝体施加的动水压力相当于一定体积的水与坝体一起来回移动,而水库的其余部分保持静止状态。上游面单位面积的附加质量可以用近似公式表示为,其中 y ≤ hw,ρw=1000kg/m3  是水的密度。在 Abaqus/Standard 分析中,附加质量方法是通过在用户子程序 UEL 中编码的简单 2 节点用户单元来实现的。由于地面运动的垂直分量产生的动水压力被认为很小,因此在所有模拟中都忽略不计。

koyna 地震期间记录到的地面加速度的横向和垂直分量如图 3 所示(单位为 g = 9.81 m/s²)。在地震激励之前,大坝受到其自重产生的重力荷载以及水库在上游面产生的静水压力的作用。

图 3.1 koyna 大坝实测地震动横向加速度时程
图 3.2 koyna 大坝实测地震动竖直向加速度时程

混凝土材料的力学行为通过混凝土损伤塑性本构模型来模拟。用于模拟的材料属性在表 1 和图 4 中给出。这些属性被认为是科伊纳大坝混凝土材料的代表性属性,并且基于之前研究人员的所用属性。在获得这些材料属性时,做了一些假设。特别值得注意的是混凝土拉伸行为的校准。抗拉强度估计为极限抗压强度 σcu = 24.1 MPa 的 10%,并乘以 1.2 的动态放大系数以考虑率效应;因此,σt0 = 2.9 MPa 。为了避免由于结构中缺乏钢筋而导致的网格敏感的不合理结果,拉伸破坏后的行为是通过指定应力/位移曲线而不是应力-应变曲线来给出的,如图 4.1 所示。这是通过后裂应力/位移曲线来实现的。

图 4.1 拉伸硬化

同样,拉伸损伤是以表格形式作为裂纹位移的函数来指定的,使用的是后裂损伤位移曲线。该曲线如图 4.2 所示。由于混凝土压缩破坏(压碎)引起的刚度退化损伤被假设为零。

图 4.2 拉伸损伤

工况 1:频域分析

混凝土坝频率分析在水利工程中具有重要意义,主要体现在以下几个方面:

一、抗震性能评估

高混凝土坝通常位于地震多发区,抗震性能是其安全运行的关键。通过频率分析,可以了解坝体在地震作用下的动力响应特性,评估其抗震能力。例如,频率分析可以帮助确定坝体在地震作用下的自振频率和振型,从而评估其抗震性能 3 。

二、动态特性分析

频率分析可以揭示坝体的动态特性,包括自振频率、振型和阻尼比等。这些动态特性对于理解坝体的动力响应和行为至关重要。例如,通过频率分析可以了解坝体在不同水位和荷载条件下的自振频率变化,从而评估其动态稳定性。

  • 下面图示为 koyna 的前 4 阶振型 
图 5.1 第 1 阶振型
图 5.2 第 2 阶振型
图 5.3 第 3 阶振型

图 5.4 第 4 阶振型
  • 下表为自振频率对比:  可以看到本次算例中的计算结果与 Chopra 计算保持一致,计算误差在可接受范围内。

可以看到本次算例中的计算结果与 Chopra 计算保持一致,计算误差在可接受范围内。

工况 3:地震响应分析-考虑动水压力

混凝土坝地震响应分析中,主要关注坝顶-坝锺的相对数据,如:相对位移数据、相对速度和相对加速数据。

*下面是坝顶-坝锺的顺河向相对位移时程曲线:

图 6 坝顶-坝锺的顺河向相对位移

*然后是坝顶-坝锺相对速度时程曲线:

图 7 坝顶-坝锺相对速度时程

*然后是坝顶-坝锺相对加速度时程曲线:

图 8 坝顶-坝锺相对加速度时程

在地震作用下,混凝土坝的材料处于循环拉压状态,混凝土会出现损伤,下面是不同时刻坝体的拉损伤云图。

图 9.1 1s 时刻坝体拉损伤
图 9.2 2s 时刻坝体拉损伤
图 9.3 3s 时刻坝体拉损伤
图 9.4 4s 时刻坝体拉损伤
图 9.5 5s 时刻坝体拉损伤
图 9.6 6s 时刻坝体拉损伤
图 9.7 7s 时刻坝体拉损伤
图 9.8 8s 时刻坝体拉损伤
图 9.9 9s 时刻坝体拉损伤
图 9.10 10s 时刻坝体拉损伤

可以看到,随着时间的增加,坝体混凝土的拉损伤逐渐累加,在地震的最后时刻,损伤达到最大,混凝土坝的拉损伤区域主要集中在坝锺和坝体上下游的折坡处。