地震波斜入射对高拱坝地震反应的影响.pdfVIP

. . 第 卷 第 期 水 力 发 电 学 报 ., 年 月 地震波斜入射对高拱坝地震反应的影响 徐海滨 ,杜修力 ,赵 密 ,王进廷 .北京工业大学城市与工程安全减灾教育部重点实验室,北京 ; .清华大学水沙科学与水利水电工程国家重点实验室,北京 摘 要:高坝地震反应分析通常假定地震波垂直入射,在显式有限元法结合黏弹性人工边界的时域波动分析方法的基 础上,建立了三维平面地震波倾斜入射的输人方法。采用所建立的斜入射输入方法,研究了地震波斜入射对小湾拱坝地 震反应的影响。分析结果表明:与垂直入射相比,地震波倾斜入射可能显著增加对拱坝抗震安全性起控制作用部位位 移、速度、加速度反应的频谱幅值放大系数,应该在拱坝地震反应分析中考虑地震波斜入射的影响。 关键词:水工结构;地震动斜入射;时域显式有限元;高拱坝 中图分类号:文献标识码:,, , . ,, ; . ,,: ,., ? ? . ,. , ,. . : ; ;; 引言 我国水电资源丰富的西南和西北地区已建和拟建米级高拱坝 余座,该地区属于高地震烈度区,因而 高拱坝抗震研究极其重要。数值仿真是研究高拱坝抗震问题的一种可行和有效手段。杜修力等 采用多次透 射公式模拟无限域能量辐射,采用可以考虑行波效应和非线性因素的显式有限元方法模拟结构和有限地基,建立 了将高坝一地基开放系统作为整体在时域内进行分析的波动求解方法。杜修力等 进一步应用建立的时域整体 分析方法研究了库水可压缩性和坝缝开合对拱坝地震反应的影响。为了避免多次透射公式的高频和低频数值失 稳问题,杜修力等又提出了一种黏弹性人工边界条件 ,并引入拱坝地震反应分析,建立了黏弹性边界结合显式 有限元法的时域整体分析方法 。 收稿日期: . 基金项目:国家自然科学基金 , ;北京市优秀博士论文指导教师科技项目 。 作者简介:徐海滨 一 ,男,博士研究生, .;杜修力 一 ,男,长江学者特聘教授,..水 力 发 电 学 报 在时域整体分析方法的基础上,合理确定地震动输入形式并且建立相应的输入方法是进行高拱坝地震反应 分析的一个关键问题 。对于工程场址附近的深源地震,入射地震波通常假定为垂直向上传播的平面体波。事 实上,当距离工程场址适度远的浅源地震发生时,地震波并不是垂直向上入射的。对于高拱坝等大体积结构,地 震波斜入射引起的地面运动的非一致变化对结构地震反应有很大影响 ,工程师们已经认识到考虑地震波入射 角度不确定性的必要性。李小军 研究了小角度斜入射时,局部介质非均匀非线性场地反应。吴兆营 研究了 倾斜入射条件下土石坝的地震反应。张如林等 利用有限元软件进行了不同角度入射条件下 波作用的时域 数值模拟计算。上述斜入射研究工作均为二维问题。 本文采用黏弹性边界结合显式有限元法的时域整体分析方法,将二维平面地震波斜入射输入方法拓展至三 维情况,并将其应用于小湾拱坝的地震反应分析,研究地震波斜入射对小湾拱坝地震反应的影响。 平面地震波三维斜入射的实现方法 . 基于黏弹性边界的地震波输入方法 本文采用显式有限元法结合黏弹性边界的时域整体分析方法,人工边界面 结点 方向的集中质量有限元总 波场运动方程为 ∑∑。衄 ∑∑“巧式中, 为结点 的集中质量; 触,和 分别为结点 方向 对于结点 方向 的刚度和阻尼系数; 和五 分别 为结点 方向.『的位移和速度; 为结点 方向 的加速度; 为在结点 方向 处截去的无限远场对有限近场的 作用应力; ,为人工边界面上结点 的影响面积。对于三维问题,,即下标 , : , , ,分别相应于直角坐 标 , , 。将人工边界处的总波场分解为内行场 上标 表示 和外行场 上标 表示 ,人工边界结点 方向 的 总位移和总作用应力可以分别写为 “ ¨“ 。 外行场采用黏弹性人工边界条件模拟,人工边界结点 方向 的应力运动关系可以写为 一 “ “ ? 五“ 式中, 和 为结点 方向 的人工边界参数,参照文献 ? 选取。 将式的外行场位移及其时间导数代人式,进一步将式 代人式,最终将式代入式,整理 得到考虑无限域辐射阻尼和地震波输入条件下人工边界结点的集中质量有限元运动方程 Ⅱ ∑∑埘 五 ∑∑ 占““ “五 : 』式中,,≠ 。 与荷载为零的有限元方程相比,式等号左边增加的两项是为了模拟无限域辐射阻尼而施加的黏弹性人 工边界条件,物理上相当于在结点 方向 施加一个另一端固定的并联弹簧一阻尼器单元,如图 所示。 式 等号右边表示地震荷载,即通过人射波得到的人工边界结点 方向 处的内行场所对应的等效荷载, 其中前两项表示产生内行场反应所需抵抗人工边界物理元件的结点力,第三项表示产生内行场反应所需抵抗近 场介质的结点力。 . 由斜入射平面波计算人工边界内行场的基本思路 上述基于