重力坝坝面动水压力的比例边界有限元分析.docxVIP

重力坝坝面动水压力的比例边界有限元分析 在地震的作用下，水库对重力坝的影响将形成对水库安全的不利动水压力，并合理确定水库面的动水压一直是国内外科学家的研究方向。从1933年的《拔格哥莱德》（1933年）到出版了一本关于计算垂直刚性水库面的动水压的经典文献后，许多科学家对动水压进行了深入研究。虽然韦加纳提出的动水压力增加的质量模型有一定的假设，但由于其简单、简单，它是许多国家抗强调整件方法所推荐的计算水库面的动水压的方法。后来的研究大多基于埃德尔菲德的想法，对动水压的计算进行了改进。ch浦a（1968）2利用势函数理论，导出了基于不同地震方向的垂直刚性水库运动压力解决方案，并考虑了水库运动的可压缩性。fenves和chuna（1983）3考虑到库水的可压缩性和库底边界的吸收，并使用有限方法求解库水的动水压。杜秀力和王进廷（1991）4回顾了动水压力的整个发展过程及其对水库的影响。 随着数值方法的不断发展进步, Wolf和Song ( 2000)［5］提出了一种新的半解析数值分析方法———比例边界有限元方法. 该方法在径向进行解析求解, 在环向仍采用有限元进行离散, 从而对所求解的问题可以达到降低一维的目的, 并且不需要边界元方法中的基本解, 可适用于各种复杂形状的结构求解. 更为重要的是, 比例边界有限元方法能够自动满足无穷远处的辐射边界条件, 可以方便应用于无限域问题. 林皋等 ( 2007)［6］、王翔等 ( 2011)［7 － 9］和李上明 ( 2011, 2013)［10 － 11］把比例边界有限元方法应用于大坝-库水相互作用中, 初步体现了比例边界有限元在该方面的优势. 本文基于比例边界有限元的方法, 对重力坝刚性坝面的动水压力进行求解, 并且考虑到了库水的可压缩性和库底岸坡边界的吸收作用, 在文献顺河向动水压力的比例边界有限元解法的基础上, 考虑了竖向地震荷载作用对动水压力的影响, 研究了动水压力沿库水上游方向的衰减规律. 1 考虑压缩性的坝面模型 对于重力坝, 假设库水为无粘小扰动理想液体, 可压缩性库水在求解域内的控制方程可用Helmholtz方程表示为 式中: ue0652为Laplace算子; p为动水压力; c为水中波速. 如若不考虑水的可压缩性, 其值将趋于无穷大. 库水的可压缩性一直以来就是研究的焦点, 当今已有学者把压缩性的研究扩展到气水孔隙混合流体中［13］. 本文后续的算例均考虑了库水的可压缩性. 坝面满足用力连续性边界条件, 即 库底和岸坡边界满足 式中: u为水平向地震动加速度; v为竖向地震动加速度; n为外法向; q表示库底和岸坡的吸收系数, 根据文献, 有 其中, α 表示库底和岸坡的反射系数, 当 α = 0时, 表示不考虑边界吸收作用, 即全反射情况, 此时库底相当于刚性, 完全没有淤沙的存在; 当 α = 1时, 则代表库底和岸坡对动水压力波的完全吸收作用. 库区表面边界条件 ( 忽略库区表面微幅重力波) 为 无穷远处边界条件为 公式 ( 6) 采用标准的Sommerfeld边界条件, 由于本文采用比例边界有限元方法, 因此, 无穷远处的边界条件是自动满足的, 可以不必引入方程进行计算, 具体的实施细节在下文中. 2 节点线性单元离散 采用比例边界有限元方法对库水进行离散, 由于库水可以理想化为沿着坝面上游方向延伸到无穷远处, 因此, 把比例边界有限元的相似中心O选在下游无穷远处即可以不离散库水表面和库底. 这里假定了坝前水深恒定不变, 从而可以使问题大为简化, 离散的区域只有大坝迎水面, 如图1所示. 用表示整体坐标, ξ, η 表示局部坐标, 其中, ξ 表示径向方向, 即水库的延伸方向; η 表示环向方向, 即垂直于库水延伸方向的坝面方向, 在此方向对求解边界进行离散, 单元离散方式如图2所示. 采用三节点线性单元对坝面进行离散. 这样整体坐标和局部坐标之间的转换可以表示为 式中, x、y表示坝面节点的坐标. 式 ( 7) 、 ( 8) 给出了库区求解域内中任意一个节点与同一深度的坝面离散节点之间的数学关系表达式. 显然, 在水平方向上, 二者仅相差一个延伸量; 在竖向则保持不变. 对重力坝迎水面用有限元方法进行离散, 即 式中, N为形函数. 从下文可以看出, 推导的比例边界有限元方程只和坝面节点坐标有关, 使问题达到降低一维的目的. 对于大坝迎水面一维线性单元, 采用三节点的二次形函数来表示, 即 这样就可以把整体坐标系转换到比例边界局部坐标系中, 其中Jacobian矩阵为 由式 ( 12) 可得到微分算子的局部坐标表达形式, 即 式中 故梯度算子可以表示为 单元面积为 式中, b1和b2分别为雅克比逆矩阵的第一列与第二列, 微分算子ue065和J及ds的计算都只和坝