有缝拱坝地震自由场输入模型及其工程应用研究.pdfVIP

水电 2006 国际研讨会 有缝拱坝地震自由场输入模型及其工程应用 张伯艳1 2 3 陈厚群 李德玉 1.中国水利水电科学研究院，北京车公庄西路20 号，100044，zhangby@ 2. 中国水利水电科学研究院，北京车公庄西路20 号，100044，chenhq@ 3. 中国水利水电科学研究院，北京车公庄西路20 号，100044，lideyu@ 摘要：为了反映坝与库水以及岩石基础的动力相互作用影响，无限地基幅射阻尼对外传波 的吸能作用，拱坝横缝在地震荷载作用下的开合运动，以及由此引起的坝体应力重分布， 本文提出了以LDDA 理论为基础的有缝拱坝地震自由场输入模型。在该模型中，层状地基 的自由场位移、速度和应力被转化为有限元边界上的作用力荷载。在求解动接触力时，为 了提高收敛精度和收敛速度，本文提出了一个新的动接触力的迭代算法，该算法的松驰因 子是依据块体接触的物理特性而非按照数学理论选用的。作为本文模型的实际应用，对拟 建于西南地区的大岗山拱坝，分别计算分析了其整体坝和分缝坝的地震响应。针对这一具 体工程，得出了以下结论：无限地基幅射阻尼可有效地降低大坝的动拱梁应力，对整体坝 而言，动拱应力降幅达 42%~52% ，动梁应力降幅达 53%~61% ；对分缝坝而言，大坝横缝 在设定地震作用下，会产生开合运动，进而引起坝体应力的重分布，横缝的张开，使拱应 力进一步减小，而梁应力有所增加。 关键词：拱坝；自由场；地震输入模型；粘弹性边界 1 前言 建立合理的拱坝地震输入模型是正确评估拱坝抗震能力的基础，为此，人们进行了不断 的探索和改进。Clough 等[1]描述了拱坝地震输入的四种模型，它们是：1）标准基底输入模型， 这一模型的特点是将坝与较大范围的可变形基础划分成有限元网格，假设基础底部是刚性的， 直接将地面运动加速度作用于刚性基底。显然，这一模型不能反应地表自由场运动与深部基 岩运动的差异。2 ）无质量地基模型，是对标准基底输入模型的改进，众所周知的 ADAP 程 序[2]使用的就是无质量地基模型。由于地基假设为无质量的，消除了地震波在岩石基础内的 传播的影响，其结果是坝与基础交接面上受到完全相同的地震动。3 ）反演输入模型，是将地 表地震动反演至深部基岩，其反演结果作用于标准基底输入模型中。4 ）河谷自由场输入模型， 是先求出二维河谷的地震自由场反应，然后将这一反应作为输入，求解三维坝——河谷地基 系统。 已有的研究表明：1)河谷地形对强震有较大的影响[3] ；2)用数值分析方法计算弹性波在无 限介质中的传播时，在截断边界上应计入无限介质的吸能效应[4,5] ，如：引入粘滞边界或人工 透射界；3)拱坝横缝在地震作用下会张开，从而导致坝的刚度减少，延长坝的自振周期，并 基金项目：国家自然科学基金90510017） 310 水电 2006 国际研讨会 [6] 释放拉应力和导致应力重分布 。因此，更为合理的拱坝地震输入，应能反映坝与库水以及 岩石基础的动力相互作用影响，无限地基幅射阻尼对外传波的吸能作用，拱坝横缝在地震荷 载作用下的开合运动，以及由此引起的坝体应力重分布。 对坝与近域地基进行有限元网格剖分，引入吸能边界条件，将拱坝横缝作为接触非线性 问题处理似乎是一个不错的选择。最近提出的人工透射边界[5]与时域显式有限元相结合的非 [7] [8] 线性地震波动分析方法 ，采用动接触力模型 模拟拱坝横缝，能考虑上述各种复杂因素对拱 坝地震反应的影响。但是，动接触力模型对接触力未进行迭代计算，因此，接触力的计算精 度可能不高。此外，人工透射边界的引入，极大地增加了有限元节点的数目，使拱坝地基系 [9] [10] 统的求解极费机时 。而以LDDA