一种新的求解坝面动水压力的半解析方法.pdfVIP

振动与冲击 No·32009 VIBRATIONANDSHOCK V01．28 OF 第28卷第3期 JOURNAL 一种新的求解坝面动水压力的半解析方法 杜建国1，林皋2，谢清粮1 (1．总参工程兵科研三所，洛阳471023；2．大连理工大学土木水利学院，大连116024) 摘 要：基于比例边界有限元理论，推导了综合考虑库水可压缩件和库底边界吸收的坝面动水压力方程，提出了 一种新的求解坝面动水压力的半解析方法，并通过箅例验证了这种方法的精度和效率。结果表明，与传统方法相比，这种 方法具有精度高、计算工作餐小的优点。 关键词：动水压力；比例边界有限元法；半解析法；水的可压缩性；库底边界吸收 中图分类号：TV312 文献标识码：A 大坝在地震作用下，库区的地面运动将会使坝面 (2)只需对求解域边界上进行有限元离散，使问 承受很大的动水压力，所以在研究大坝的地震响应时， 题空间维数降低一维，从而节约了一定的计算量； 不得不考虑坝面动水压力对大坝的影响。早在上个世 (3)由于无限远处边界条件已经得到满足，从而 不用基本解，使得可以解决由于得不到基本解而无法 纪30年代，Westergaard…就研究了刚性直立坝面的地 震动水压力问题，并提出了附加质量的概念，至今仍被 利用边界无法求解的问题； 工程界所接受和采纳。然而，坝面动水压力问题是一 (4)基本理论是以有限元法为基础的，容易被科 个非常复杂的问题，大量有关研究表明：库水可压缩 研及工程应用人员接受和掌握，并且可以方便地与有 性、坝体特性、淤砂层特性、地基特性以及库水一坝体 限元法结合来解决问题。 相互作用、库水一地基相互作用、库水一库底淤积砂层 鉴于以上原因，本文中将比例边界有限元理论应 相互作用、淤砂层一坝体相互作用、淤砂层一地基相互 用于动水压力的求解，详细推导了综合考虑水的可压 作用和坝体一地基相互作用等因素都会对坝面动水压 缩性和库底边界吸收的坝面动水压力比例边界有限元 力产生影响心]。其中库水的可压缩性(包括坝体与库 控制方程，实现了对库水的模拟，从而提出了一种基于 水的相互作用)、淤砂层和地基的特性(包括库水与淤 比例边界有限元理论的求解坝面动水压力的半解析方 砂层和地基的相互作用)是对库水动水压力计算产生 法，并通过算例对这种方法进行了验证。 影响的两个重要因素。 1坝面动水压力问题的完整描述 目前，对坝面动水压力的计算多采用有限元法p’41 或边界元法”’6J。采用有限元法计算时，为了满足远处 如图1所示，以坝前的库水为求解域，整体Carte— 无反射的边界条件即Sommerfeid辐射条件，一般要引 人人工边界，这就使得计算量增加。如果不引人人工 小扰动的理想流体。同时，为了把本文的重点放在方 边界，库水的范围则应大于一3～5倍坝高以上，在这样 法的推导上，假定水面与库水边界平行，对于复杂的库 大的范嗣里进行库水单元离散化，形成关于动水压力 底和岸坡形状的模拟，将另外研究讨论。 的有限元代数方程组，其计算工作量是相当可观的。 Z 采用边界元法进行求解，由于要求基本解，其计算程序 n S 。 的复杂性往往大大超过有限元法。同时这种方法得到 一一，