考虑辐射阻尼拱坝随机振动分析.pdfVIP

水 利 学 报 2005 年1 月 SHUILI XUEBAO 第1 期 文章编号: 0559-9350(2005)01-0110-06 考虑辐射阻尼的拱坝随机振动分析 吴健，金峰，徐艳杰 (清华大学 水利水电工程系，北京 100084) 摘要：本文基于有限元-边界元-无限边界元耦合的时域模型，提出了结合随机振动虚拟激励法的拱坝-地基动力相 互作用分析模型，并用蒙特卡罗法对模型进行了验证。对某拱坝的计算结果表明，辐射阻尼对拱坝随机动力响应 有显著影响，结果同时表明，地震动的随机性对于拱坝动力响应的影响明显，用有限的几条确定性地震动时程输 入不能全面反映地震动的随机性。 关键词：随机振动；无限地基；辐射阻尼；虚拟激励法；拱坝-地基动力相互作用 中图分类号：TV642.4 文献标识码：A [1～3] 拱坝和岩石地基的动力相互作用是涉及工程抗震安全的重要研究课题，一些学者 的研究表明，无 限地基的辐射阻尼不应该被忽略。Lysmer和Richart[4]首先提出了用三参数模型来模拟半无限地基；Wolf[5] 发展了一种系统化的方法，使用分式展开建立连续集中参数来模拟无限地基，得到近似的地基阻抗函数； 张楚汉和金峰[6]提出了FE-BE-IBE时域模型，将与频率相关的地基阻抗函数转换为与频率无关的质量-弹簧- 阻尼三参数模型，然后耦合到时域内的拱坝动力计算程序中，该模型采用自由场地震输入，能够考虑无限 地基辐射阻尼以及地震动输入的空间差异对拱坝动力响应的影响；杜修力等[7]也采用显式有限元-外推人工 边界方法分析了拱坝与地基的动力相互作用，可以考虑地基的辐射阻尼效应、坝基不均匀地震动影响等。 [8，9] 另一方面，地震动是不确定的随机过程，国内外众多学者 ， 对于结构随机振动响应做了大量的研究，但大都运算量过大， 难于应用于工程实际。陈健云[10]等研究了拱坝-地基动力体系 的随机动力响应，计算中需要位移的功率谱密度函数，对功率 谱的形式要求较高，计算过程比较烦琐，为了减少工作量需要 忽略各振型之间的相关性，给结果带来一定的误差；林家浩[11， 12]等从计算力学的角度提出虚拟激励法，由于其形式简单，便 于编程，在结构随机响应计算中得到了较为广泛的应用；梁爱 虎[13]等利用虚拟激励法计算了拱坝的随机地震响应，并给出了 可靠度的结果，但该模型中没有考虑到无限地基辐射阻尼对拱 坝动力响应的影响。本文采用FE-BE-IBE理论模型结合虚拟激 励法，对拱坝的随机动力响应做了分析。 图1 拱坝-地基-库水系统离散示意 1 FE-BE-IBE数值模型 如图1，拱坝坝体用有限元离散，地基用边界元和无限边界元划分网格。采用子结构方法分析系统的 相互作用，把地基与坝体分开考虑，便于确定地基的阻抗函数。 收稿日期：2003-10-27 作者简介：吴健(1978-)，男，安徽宣城人，博士生，从事结构抗震分析方面的研究。 1 水 利 学 报 2005 年1 月 SHUILI XUEBAO 第1 期 将地基阻抗函数转换成质量-弹簧-阻尼三参数模型，得到地基的有限元动力方程： M cc uc +Ccc uc +K cc uc F c