高水头平面闸门流激振动数值模拟.pdfVIP

现代水利水电工程抗震防灾研究与进展 高水头平面闸门流激振动的数值模拟 刘权1张燎军1关超年1’2 (1．河海大学水利水电学院，江苏南京210098 2．中国水电顾问集团中南勘测设计研究院，湖南长沙410014) 摘要：高水头下平面闸门在启闭过程中受水流脉动压力作用可能引起强烈振动而容易发 生破坏。本文以某水电站导流底孔中工作闸门为例，采用势流体流固耦合理论，利用 ADINA有限元软件建立水体一闸门流固耦舍模型。通过输入水力学模型试验实测脉动 压力时程曲线，计算分析了闸门局部开启时的流激振动特性，研究了闸门振动位移、应 力、速度、加速度等分布规律。讨论了闸门相关振动控制标准，据此进行了闸门的振动 安全评价。研究成果可以为高水头平面闸门振动设计分析提供参考。 关键词：平面闸门；高水头流激振动；流固耦合；数值模拟；ADINA 0引言 某水电站大坝坝顶高程为1767m，最大坝高158m。按照计划，深埋在山体的两条施工期导流 洞将被改造成一条主底孔和一条次底孔，以备在必要时放空水库。其中在距离主底孔进水口约 500m的长隧洞内部布置了一扇工作闸门和一扇事故闸门，皆为平面型式。按设计要求，工作闸门 必须能够在运行水头高达130m的动水中启闭。当闸门局部开启泄水时，高水头有压流水流条件 复杂，水流流速可达到40m／s，可能引起闸门产生较大的振动。另外，与国内一般平面闸门的设 计不同，该平面闸门挡水面板位于下游侧，这就意味着其梁格系统完全浸泡在水中，在启闭时上 游复杂的水流条件引起的脉动压力可能会加剧闸门的振动响应。根据对国内外许多平面闸门破坏 事故分析u’，平面闸门多是在启闭过程中发生破坏，并往往伴随着剧烈的振动，容易诱发共振比1。 随着高库大坝日益增多，对于这类闸门在超高水头长隧洞有压流作用下产生的流激振动问题口1的 研究意义已愈加重要。 本文基于ADINA有限元软件针对该平面钢闸门建立三维水体一闸门流固耦合数值仿真模型， 结合水力学模型试验成果，研究此闸门在局部开启时过水的流激振动特性，并对闸门进行振动安 全评价，为高水头闸门设计提供科学依据。 I流固耦合理论 流固耦合有限元分析理论是解决深孔平面闸门流激振动特性的有效方法。 根据文献[4]，在势流体中，连续方程近似采用标准Galerkin方法获得，并可以改写为以下 形式： (1) 艿‘=ff妥拿五却．∥V≯．V彩py+f-p五．hands 多＼、On ／ ； 式中，万E是质量流量率变量，S是流体域y的边界，，z是边界面S向内的法线，p是密度，h 是特定的焓(定义为h=I咖／p)，P是压力。如果没有边界条件施加在边界面s上，则暗示着 边界条件pi：·n=O，即没有流体流过这个边界。所以pil·刀叫做自然边界条件，并利用这个边界 条件把结构的运动施加到流体域内，同样无穷流水域也可以使用这个边界条件来模拟。 方程(1)通过线性化可以得到只考虑势流体单元对系统矩阵的贡献的有限元运动方程，如 考虑结构项的贡献，可得到整个系统的流固耦合方程为陆1； [警M口rr煳A≯+[肾怕％一阱 201 现代水利水电工程抗震防灾研究与进展 J‰+K ％ ]～1f‰1f毛 ] 心’ 【- 置凡， l△≯j 一(盔0F+(j0．，)s)J【口j 【j_+(，’F)sj 式中：未知势矢量矽的增量写为△≯，未知位移矢量H的增量写为△Ⅳ；互，是外力矢量，B和 (B)s分别来自方程(1)的体积积分项和面积积分项；朋o、％、墨F，分别为流体的质量 矩阵、阻尼矩阵、刚度矩阵；朋基、Ck、如、‰分别为结构的质量矩阵、阻尼矩阵、刚 度矩阵和荷载矢量；cUF和％代表流体与固体的耦合阻尼矩阵，k、％代表流体