满库状态和强震作用下戈兰滩重力坝的损伤研究论文.docxVIP

满库状态和强震作用下戈兰滩重力坝的损伤研究论文 0 引言 重力坝是挡水部分的重要部分，大量的震害调查发现，重力坝地震作用下在大坝不弱部位产生裂缝损伤，受损部分的混凝土进入受损，开裂后退出工作会产生应力释放，对大坝的挡水功能造成一定的影响。如在四川汶川地震中[1],多数重力坝都产生混凝土开裂等损伤。在研究挡水重力坝的动力响应时，库水与重力坝接触面由于地震会产生共振作用，1933 年，Westerguard[2]在对直立坝面的重力坝提出了动水压力的计算公式。徐汉忠[3]应用边界元法将韦氏公式加以扩充用以计算挡水面倾斜时的动水压力，对正常水位下的动水压力计算有较好的效果。王铭明[4]则针对不同坝高考虑坝体弹性修正系数、库底吸收修正项对 Westerguard 公式进行改进，得出在不同高度的压力下，上游坝面的流固效应在坝体下部会有所减小。刘依松[5]对重力坝的动力分析中采用二维无质量地基模拟，得出坝体反映峰值滞后地震加速度峰值，动力响应与地震加速度变化规律相似。程冬[6]采用二维模型模拟丰满重力坝的动力分析中指出重力坝上部转折点为损伤薄弱处。谯雯[7]采用谱分析法研究口水重力坝的模态分析中得出对于目前建成的挡水重力坝在震级不强的地震作用下，其抗震性能较好，混凝土一直处于线弹性阶段。目前，在坝工领域运用最为广泛的数值计算是有限单元法[8],然而有限元方法是一种基于连续介质力学的方法，认为岩体是连续的。本文利用 ansys 分析软件，建立挡水重力坝三维立体模型，考虑地基质量、材料损伤引起的`刚度退化和应力释放，对戈兰滩重力坝在满库状态的强震作用下进行动力响应研究和损伤区域研究。 1 工程实例 1.1 工程概况 戈兰滩水电站位于云南省李仙江流域，挡水重力坝坝长466m,顺流方向长度为 95m,高度 113m,正常蓄水位 103m,结构设计抗震烈度等级为 8 度。图 1 是挡水重力坝的结构示意图。 1.2 工程有限元模型 1.2.1 有限元计算模型 坝体选用了 solid65 单元，该单元是专为混凝上和岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元， 用于各种三维实体模型；坝基选用 solid45 模拟，计算时假定坝体和坝基的材料是均匀的，两者之间紧密连续在一起；基岩则采用线弹性本构模型。有限元模型如图 2,计算模型中坝体建立13680 个单元，坝基建立 48918 个单元。 1.2.2 荷载类型 （1）动水压力荷载。在地震作用下，库水与坝上游面有耦合作用，产生流固耦合效应[9].在计算上游面的动水压力时，由于上游坝面的挡水面是倾斜的，为了准确地表达水压力，采用徐汉忠[3]利用边界元法将 Westergaard[2]公式扩充以后的动水压力计算公式来表达： 式中：V咬n为法向加速度；hw为水深；y′为从坝基算起到计算处的高度，β 为上游坝面的倾角。 （2）地震加速度荷载。对于地震荷载，采用人工合成地震波。地震波时程采用由规范普推演出的规范普地震波加速度时程曲线。按照规范规定，大坝动力分析应考虑垂直坝轴线方向的水平地震作用和竖向地震作用，其中加速度峰值 0.2g,在计算中竖向地震动和水平地震动同时输入，地震时程曲线见图 3,图 4. 2 地震动力分析 2.1 动位移分析 图 5~ 图 8 是整个地震过程中挡水重力坝在 X 方向（顺河）和 Y 方向（竖直）最大位移变形图。从分析结果可以看出，在地震作用下，挡水重力坝的位移变形主要表现在 X 方向（顺河），坝顶跨中位移最大值 30.7mm,转折点位移最大值11.8mm,坝肩与地基连接处最大位移 1.4mm;在 Y 方向（竖向）由于结构自身自重的原因，在一定程度上限制了自身的震动，坝顶跨中最大位移为 7.5mm,转折点最大位移6.6mm,坝肩与地基连接处最大位移 0.67mm. 通过软件对坝顶，坝踵转折处坝肩与地基连接处的位移随时间的变化，得出挡水重力坝关键部位的位移时程变化曲线如图 9,图 10 所示。由于重力坝在转折点上部的结构的刚度较小，下部结构刚度较大，坝顶动位移响应比坝踵动位移响应大，坝踵比坝肩动位移响应大，位移迟滞效应较明显，位移最大值出现在 3.8s.Y 方向个节点的位移曲线重合度较高，表明各个部位在 Y 方向（竖向）的迟滞效应不明显，动力响应较小，最大值出现在 4.1s,由于自身自重的原因，位移峰值比X 方向（顺河）动力响应滞后。 2.2 动应力分析 图 11~ 图 14 是挡水重力坝在地震动力和满库动水压力作用下叠加后的主应力图。从分析的结果来看，在地震作用下，由于坝顶震