泄水陡槽三维流固耦合问题的解析-有限元法.docVIP

泄水陡槽三维流固耦合问题的解析-有限元法* 谢省宗、吴一红、方咏来**、薛小香 [提要] 本文采用解析-有限元联合解法,对泄水陡槽(或渡槽)三维流固耦合问题进行了研究，求出了在平行水流作用下矩形断面泄水陡槽结构振动时的三维速度势φ 和动水压力p的解析解,及与结构有限元分析相对应的解析解;同时提出了陡槽结构流固耦合振动时,动水压力的附加质量矩阵、附加阻尼矩阵和附加刚度矩阵，流固耦合时，流体对结构自振特性的影响化为这些附加矩阵的影响。 [关键词] 泄水陡槽,流固耦合,解析-有限元联合解法,附加质量矩阵、附加阻尼矩阵、附加刚度矩阵。 1.前 言 高速水流作用下泄水陡槽（或渡槽）简单结构的二维耦合问题，可以有解析解[1][2]。对于三维复杂结构流固耦合问题的分析，目前大多采用 纯数值解法，其中尤以有限元法为代表，它对流体和结构都采用有限元离散。有限元法虽然适用范围广，但当用它解流体为无限域或半无限域时，必须引进人为的外部有限边界进行截断，这样就会带来一定的误差; 特别是当流体含有流速场时,用有限元法求解流固耦合问题时,难度很大。由于流体域的边界相对较为简单,有时其流场可以用解析法求解,而结构型式一般较为复杂.通常只能用有限元法求解,则流固耦合问题适合用本文所提出的解析-有限元联合解法，其特点是对结构部分采用有限元法，对流体部分采用解析法,也就是说，在解流体方程时，通过严格的数学方法，求得作用在流固耦合面上的动水压力，只是耦合面的边界条件用的是有限元离散近似；而在解固体方程时,采用的是离散型的有限元方程,但作用在固体边界面上的流体荷载却是通过严格的数学方法求得的。解析-有限元法对于复杂结构流固耦合的特征值问题或动力响应问题，在低频和高频范围内都能得到足够精确的结果。 考虑到我们所研究的是高速水流作用下泄水陡槽(或渡槽)的三维流固耦合问题,其流场两端延伸到无穷远，断面为一矩形规则断面，泄水陡槽的自振频率较低，且要求论证在高速水流作用下的动力稳定性,所以适于用解析-有限元法进行分析。 2.流体运动的基本方程及其解析解 研究如图1所示的泄水陡槽三维流固耦合问题，陡槽为一由底板和两个侧墙(分别为侧墙1和侧墙2) 组成的矩形规则断面,槽内水深h,水流流速U平行于底板,取一个长为L的分段进行分析,建立的坐标系如图1所示。 图1 陡槽――水流示意图 *本文为水利电力重点科技项目《溢流混凝土面板堆石坝关键技术开发》（编号：92-SZ04）的部分内容，项目曾获2004年新疆维吾尔自治区科技进步二等奖。 **现在浙江省水利勘测设计研究院 结构作微幅振动时,流体运动考虑为理想不可压缩无旋运动,其扰动速度势 φ满足Laplace方程和相应的边界条件: (1) 自由表面处： （1a） 底板处： （1b） 侧墙1处： (1c) 侧墙2处： (1d) 两端无穷远处： (1e) 其中w(x,y)为底板在z方向的振动位移,v1(x,z)、v2（x，z）别为侧墙1、侧墙2在y方向的振动位移, ω为振动圆频率。因各物理量都以简谐方式()随时间t而变化，所以所有的分析都只考虑它们的幅值。 上述方程和边界条件都是线性的，所以 可看作底板、侧墙1和侧墙2分别振动时产生的速度势、、和的叠加，即： (2) 它们分别满足Laplace方程和各自的边界条件（表1）。 表1 、和所满足的边界条件 自 由 面 处： 0 0 0 底 板 处： (0≤x≤L) 0 0 侧 墙1 处： 0 (0≤x≤L) 0 侧 墙2 处： 0 0 (0≤z≤L) 无穷远处: 0 0 0 0 0 0 先用一次无限付氏变换和一次分离变量法,再加付氏逆变换得到它们的解析解: （3） （3a） （4）