基于人工边界方法的西藏旁多土石坝非线性动力分析

2021-07-03郭永刚

中国农村水利水电 2021年6期

吴悦，郭永刚，胡锦

（西藏农牧学院，西藏林芝860000）

0 引言

西藏地处高寒、高海拔、强震地区，土石坝的抗震安全成为关键问题。坝体结构的抗震研究多采用动力模型试验、数值模拟分析等方法［1-3］，由于坝体-库水-地基相互作用的模型试验技术尚未成熟，学者多采用数值模拟的方法来研究这一问题［4］。在土石坝非线性动力仿真分析中正确实现地震动的输入和地基辐射阻尼的模拟至关重要，人工边界的引入可以有效地解决地基无限域问题［5］，国内外研究者基于波动理论建立的不同类型人工边界主要分为全局人工边界和局部人工边界，无限元人工边界是全局人工边界的代表之一，无限元思想由Ungless［6］于1973年首先提出，后由Bettess、Beer 等［7-9］发展，许多学者研究经验［7-9］表明，无限元比边界元等其他数值方法求解无限域问题更具优势和实用性；局部人工边界包括黏性边界、透射边界、黏弹性边界等，黏弹性边界由Deeks 等［12］、刘晶波［5］等最先提出，是在黏性边界的基础上增加弹性元件以模拟远域地基的弹性恢复性能，具有良好高频和低频稳定性，应用广泛。目前土石坝动力响应分析中对于多种人工边界的对比研究较少，本文以西藏旁多水电站工程实际为分析对象，采用数值模拟的方法比较了固定边界、无限元边界、黏弹性边界条件下坝体结构动力响应结果，分析讨论不同的人工边界条件对模拟结果的影响，以期对数值分析中人工边界条件的选取提供依据。

1 人工边界条件

1.1 无限元人工边界

无限元单元的特点是坐标变换和位移模式采用不同的插值函数，其几何描述由一组映射函数实现，位移函数则采用与普通等参元相同的形函数逼近［13］。无限元单元具有方向性，其延伸方向是从近场到远场，趋于无穷远的方向上不能相交。如图1所示，节点1～4 与近场有限元节点相连，5～8 为无限元中间节点，其余4个节点在无穷远处。

节点编号顺序不同所对应的形函数和映射函数不同，其中形函数Ni为：

单向映射无限元映射函数Mi为：单元刚度矩阵按下式计算：

式中：ξ、η、ζ为局部坐标；［D］为弹性矩阵；［B］为应变矩阵，[B]=[B1，B2，B3，B4，B5，B6，B7，B8]。

无限元通过在边界上引入阻尼力吸收平面体波的辐射能量，假定无限元材料为线弹性，阻尼产生的应力为：

一般取系数为：

式中：Cp、Cs分别为纵波、横波波速。

ABAQUS 非线性有限元分析软件中不能对无限元单元直接定义，可先将无限元单元定义为某种单元类型（如C3D8R），随后在INP文件中修改单元类型为（如CIN3D8），从而实现对无限元单元的定义。

1.2 黏弹性人工边界

黏弹性边界是在黏性边界的基础上增加弹性元件以模拟远域地基的弹性恢复性能，具有良好高频和低频稳定性。在ABAQUS 非线性有限元分析软件中黏弹性边界条件下的地震动输入一般转化为人工边界节点上的等效节点力处理［14］。

式中：Ab为边界节点的影响面积；n为边界外法线方向余弦向量。

刘晶波等［15］将黏弹性人工边界节点的弹簧系数KB和阻尼系数CB取为：

切向边界：

法向边界：

式中：KBT和KBN分别为弹簧法向和切向刚度；R为波源到人工边界点的距离；CP和Cs分别为波和S波波速；G为介质剪切模量；ν为泊松比；ρ为密度；αT和αN分别为切向和法向人工边界参数。

根据谷音等［16］推荐的参数α的取值范围，文中取αT=0.67、αT=1.33。黏弹性边界模型示意图见图2。

2 工程概况及有限元模型

西藏旁多水电站实际工程位于海拔4 000.00 m 以上，大坝结构为碾压式沥青混凝土心墙砂砾石坝，大坝最大断面见图3。

本文的土石坝动力数值计算中，坝体左、右各取2 倍坝高；上、下游坡脚向上、下游各取3 倍坝高，根据旁多土石坝设计文件取覆盖层深度为132 和264 m 分别作为计算区域；并将地基边界分别做固定边界、无限元边界和黏弹性边界处理；计算工况为竣工期/蓄水位+地震动。三维模型有限元单元类型为C3D8，无限单元类型为CIN3D8。无限元和黏弹性边界条件下的三维模型网格划分见图4。

3 动力分析

动力分析采用迁安波作为地震动输入波，在模型底部边界施加双向地震波，在模型顺河向和横河向边界分别施加水平向地震波，地震波加速度时程采样间隔0.01 s，持时20 s。竖向边界施加的地震波加速度峰值取水平向的2/3，计算中输入的地震波幅值为0.1 g，如图5所示，并按比例放大得到幅值为0.3、0.4 g的地震波。

图6为覆盖层深度为132 m 时，黏弹性边界条件下的加速度时程曲线。通过与输入地震动时程曲线（图5）对比，可以看出坝体的加速度响应时程曲线与输入地震波时程曲线拟合度较高，且坝体加速度响应随着输入地震波峰值的增加而增大。

图7为覆盖层深度为132 m 时不同边界条件下的加速度峰值。随着输入地震波幅值的增加，坝体的加速度响应峰值也随之增加。0.1、0.3、0.4 g 地震波作用下，蓄水期竖向加速度峰值相对于竣工期竖向加速度峰值均有明显减小，固定边界条件下分别减小了6.1%、15.7%、13.4%；无限元边界条件下分别减小了10.0%、12.2%、12.7%；黏弹性边界条件下分别减小了12.5%、13.4%、12.8%；说明蓄水期库水荷载对竖向加速度的响应有明显的抑制作用。

图8为覆盖层深度为264 m 时不同输入地震波作用下的加速度峰值。可以看出在0.1和0.3 g地震波作用下库水荷载对加速度响应的影响较小，黏弹性边界条件下几乎没有影响，可以忽略不计，可见，较深的覆盖层对地震引起的动水压力有较强的吸收作用，从而较大程度上削弱了库水荷载对坝体加速度响应的影响。地震波幅值在0.4 g时，库水荷载使得无限元边界条件顺河向加速度峰值有明显增加，黏弹性边界条件下增加较小。

图9为覆盖层深度为132 m 时不同输入地震波作用下的顺河向加速度峰值放大系数，总体上固定边界条件下加速度放大系数最大，无限元边界次之，黏弹性边界最小，说明黏弹性边界条件可以更好地模拟地基的辐射阻尼效应。

图10和图11分别为不同覆盖层深度下加速度放大系数沿坝体高程的分布，覆盖层深度为132 m 时，固定边界、无限元边界、黏弹性边界条件下，放大系数的最小点分别位于坝高1/3处、1/6 处、坝底，逐渐向建基面靠近，说明黏弹性边界条件下的数值模拟结果更加精确。覆盖层深度为264 m 时，黏弹性边界与无限元边界的加速度响应结果基本一致，尤其是蓄水期，说明黏弹性边界条件的适用范围更广。