胶凝砂砾石化石坝抗震性能试验研究.docxVIP

胶凝砂砾石化石坝抗震性能试验研究 该胶凝砂砾坝是在总结传统板堆坝和混凝土瓣效应的基础上发展起来的。具有良好的安全、可靠性、经济合理性。被认为是抗旱性好的水库。它被认为是抗旱性好的水库。它特别适用于地震多发区，具有很高的竞争力。该坝型是由堆石料、砂砾料、胶凝材料等填筑而成, 其坝料的物理力学性质极其复杂, 已有研究成果对该筑坝材料的动力力学性质和变形特征认识尚存在不足, 从而阻碍了该坝型的实际应用。目前, 关于胶凝砂砾石坝的动力响应研究主要采用剪切楔法, 将坝体假定为底部嵌固在基岩上的变截面线弹性梁[2~4], 其分析结果可在一定程度上反映该坝型的动力响应特性, 但由于剪切楔法只能计算坝体各水平截面的动力反应平均值, 不能给出水平截面上不同部位的动力反应值, 且胶凝砂砾石坝坝料具有明显的非线性力学特性, 将其假定为线弹性材料不能客观地反映大坝的应力应变特征。因此, 该方法不能很好地满足该坝型的实际应用。鉴此, 本文针对目前胶凝砂砾石坝动力响应研究中存在的不足, 采用数值方法, 基于动力三轴试验提出的胶凝砂砾石料动力本构关系, 对胶凝砂砾石坝进行地震反应计算, 研究该坝型的动力响应特征, 旨在为工程应用提供参考依据。 1 橡胶凝胶固结石块的动态本构关系 本文分别拟合了胶凝含量为40、60、80kg/m3三种胶凝砂砾石料 (简称C40、C60、C80) 对应该模型的试验参数。 1.1 材料的动力材料n 胶凝砂砾石料的非线性动弹性模量为: 式中, Ed为动弹性模量;a、k、n均为试验拟合参数;pa为标准大气压;σ3为单元最小主应力;εd为动应变。 三种胶凝含量的胶凝砂砾石料对应于式 (1) 的动力材料参数见表1。其中Ra2, Rk2,n分别为参数a和k、n的拟合曲线回归系数。 1.2 胶凝砂的阻尼比参数 胶凝砂砾石料的非线性阻尼比为: 三种胶凝含量的胶凝砂砾石料对应于式 (2) 的阻尼比参数见表2 (其中R2为参数d和λmax的拟合曲线回归系数) 。 2 baqus二次开发功能 本文根据胶凝砂砾石料动力本构关系模型, 借助ABAQUS的二次开发功能, 编写了胶凝砂砾石料的动力本构关系UMAT子程序, 对该坝型进行动力计算分析, 可较准确地反映大坝的实际地震响应特征。 2.1 采样时间和加速度 按照规范选取2条人工地震波A1、A2, 见图1。这2条地震波震动时程的基准期为100年, 超越概率为2%, 加速度峰值为2.88m/s2, 时程采样间隔0.02s, 持续时间为40s。计算中顺河向和竖直向输入的地震波分别为A1和A2, 并根据规范, 将竖直向的加速度峰值调整为水平向的2/3。 2.2 数值模型的建立 目前已建的胶凝砂砾石坝, 坝料胶凝含量多为50~70kg/m3, 断面形状为上下游对称的梯形, 上下游坝坡坡比多采用1.0∶0.6~1.0∶0.8。为研究胶凝砂砾石坝的地震响应特征, 本文同时考虑胶凝含量和坝坡变化等因素的影响, 建立均质胶凝砂砾石坝有限元模型见图2。模型坝高100m, 正常蓄水位90m, 坝顶宽10m, 坝体断面为上、下游对称的梯形断面, 上、下游坝坡坡比分别选取1.0∶0.5、1.0∶0.6、1.0∶0.7、1.0∶0.8、1.0∶0.9, 筑坝材料分别采用胶凝含量为40、60、80kg/m3的胶凝砂砾石料, 材料参数采用本次试验结果数据。为简化计算, 不考虑上游面的混凝土防渗面板, 而将水压力直接作用于上游坝面上, 采用ABAQUS的单元移除和激活功能模拟大坝的分级填筑, 大坝模型分20级施工。地震动力荷载按顺河向和竖直向加速度时程曲线输入, 每0.02s为一个荷载步, 有限元模型的坐标方向以顺河向为x轴, 向下游为正;以竖直向为y轴, 向上为正。 地震动水压力采用Westergaard提出的库水模型理论计算, 上游坝面倾斜的堆石坝库水的附加质量为: 式中, φ为上游坝坡与水平面的夹角;H0为节点i所在断面从水面到库底的水深;yi为水面到节点i的水深;Ai为节点i的控制面积。 由式 (3) 计算胶凝砂砾石坝模型坝面上所有节点, 即可得到库水的附加质量矩阵。针对动水压力的处理方法, 基于式 (3) , 利用ABAQUS中的集中质量单元模拟库水的附加质量, 可得到较精确的结果。 3 计算与分析 3.1 坝体自振周期与地震卓越周期相结合 坝体的基频和自振周期是分析坝体动力反应的重要因素之一, 坝体的动力响应特性和地震波频谱特性密切相关, 若预先了解大坝的自振周期, 尽量避免坝体自振周期与地震卓越周期重叠, 可减小地震损伤。 3.1.1 自振周期和频率 以C60为例, 计算上、下游坝坡坡比分别为1.0∶0.5、1.0∶0.6、1.0∶0.7、1.0∶0.8、1.0∶0.9的胶凝砂砾石坝模型的自振周期和频率, 结果见表3