有缝重力坝-地基系统地震反应分析方法.pdfVIP

有缝重力坝一地基系统地震反应分析方法 涂劲 杜修力 陈厚群 t中雷水利木电科牛研究院工程抗震研究中心，北京looD“) 提要 本文针对目前重力坝地震反应分析中存在的一些问题，提出了一种考虑无限地基的辐射阻尼作 用．考虑坝俸分缝时的接触非线性的静动组合的地震反应分析方法，并根据这一方法的思想绾制出有限 元分析程序，通过实际工程的分析计算验证其有效可行． 关键词重力坝。人工透射边界，动接触力法，静动组合 1引言 重力坝地震反应分析通常是采用无质量地基模型的有限元法和悬臂梁法，这种方法通 常只能考虑结构为线性的情况，而且忽视了无限地基的能量辐射效应．近年来在重力坝和 拱坝坝体系统的分析计算中出现了多种试图考虑无限地基作用的方法，金峰等““采用的 有限元一边界元一无限边界元耦合模型，能够充分考虑近场土的非均匀性和远场辐射条 件，但这是在频域内的耦合模型，无法考虑坝体及地基中存在缝隙及软弱夹层等非线性的 情况；赵崇斌等…采用无限元与有限元耦合，使得有限元的离散区域缩小到工程中感兴 趣的最小范围，用无限元法进行动力计算目前主要局限于频域内的求解．由于这些方法都 涉及到解析解的求驭，而解析解一般都是在频域内获得，近年来一些学者致力于将以上频 域内的计算模型转化到时域内的工作，如Ⅱ3mⅡ4目，其要点是将频域内表示无限地基 作用的地基动阻抗SI∞)用f国的多项式进行拟合，将其转化为一个由等效集总参数组成 的、与频率无关的力学体暴—附加于边界结点的弹簧、质量、阻尼器，从而建立起考虑无 限地基作用的运动方程，采用一般的逐步积分方法即可求解．但是应当注意到，这一方法 存在以下明显的缺陷：1．由频域向时域转化进行曲线拟合的过程中，拟合点的选取是人 为的．目前采用的是用抛物线模拟动阻抗实部，用直线模拟虚部的方法，因此两个拟合点 国．，珊，的选取直接影响到所得集总参数的取值．2．对于坝基地震动输入问题，若需要考 虑河谷地形对地震波的散射作用，则要求解地震波散射场，散射场的求解只能在频域内进 行，即先将半无限空间自由场{《lf)}进行傅氏变换，得到{W(∞l}，通过求解得到散射场 的处理采用的仍是频域内的动阻抗概念，无法处理地基内存在非线性的情况． 为克服上述问题，本文采用时域显式有限元与人工透射边界相结合的方法考虑无限地 基的辐射阻尼作用，求孵重力坝一地基系统的地震反应．对于坝中的纵缝则采用动接触力 模型来模拟缝面的动态开合情况．同对，由于注意到在地震作用之前，由于库水、淤沙和 潦劲．1973年9月l日出生．博士 ·374· ‘工程力学)增刊 1999年 坝体自重作用，坝体一地基系统已处于受力状态，缝界面的接触状态已经改变，不再适合 采用直接叠加的方式计算结构最终的静动组合反应，因此本文提出一种静动连续计算的组 合方法，来进行有缝重力坝一地基系统的地震反应分析。 2显式有限元一人工透射边界方法 显式有限元结合人工透射边界这一解耦的时空离散近场波动数值模拟方法，不仅考虑 了地基中地震波动能量向无穷远处的逸散，而且能处理地震动输入的不均匀性及基岩介质 非均匀、非线性等情况．同时，由于其“解耦”特征，处理简单，精度可控，且无需求解 大的方程组，这对于求解自由度数目极大和介质呈非线性的问题优势极为突出，因此我们 将它用于重力坝一地基系统的地震反应分析。 显式有限元一人工透射边界方法首先将结构和感兴趣的部分地基切割出来进行离散 化，称为近场计算区，在近场计算区的地基外附加一定厚度的线性介质构成人工边界区， 人工边界区的外部结点称为人工边界结点，人工边界以内的所有结点称为内点． 2．1显式有限元的内点计算方法 采用集中质量、集中阻尼的有限元模型，用中心差分法得到内点解耦的显式积分格式： (寿m，+赤。，Ⅳt”“=R卜Ft‘+孛mtuj一(寿mt一壶c，Ⅳj“’ (2一1) 2．2人工边界结点位移计算方法 2．2．1散射场引起位移的计算 散射波是外传波，由不同入射角度和传播速度的若干子波组成，假定人工边界区介质 是线性的，则外传波可表为n个外传子渡沿人工边界外法向的视传播渡的线性组合。 毗蹦)=弘●一晶，0 现取一常值人工波速c。近似替代c，／cosp；，于是有 ue，y，f+△f)z∑ulG—c。△f，y，f)=ue—c。△r，y，f) 式(2