某水库上游围堰的流固耦合分析研究.pdfVIP

某水库上游围堰的流固耦合分析研究 景 亮 张青峰 (宝鸡市冯家山水库管理局 陕西 宝鸡 721300) 摘 要 水利工程中，坝体渗流场与应力场的分析越来越受到关注，但是渗流和应力是相互关联密不可 分的，单独分析应力场或渗流场会带来很大的误差。本文选用某水库上游围堰土石坝坝体的典 型剖面，采用Geo—Studio软件中的渗流分析模块 SEE 和应力应变分析模块 SIGMAW／联合，实 现流固耦合即应力场与渗流场的耦合分析 ，以期较真实地反应围堰的渗流与应力情况，为工程 设计提供参考。 关键词 耦合；应力；渗流；土石坝 E—H 中图分类号：1、，551．3 文献标识码：A 在水电工程中，坝体的应力与渗流情况 时，有限元划分网格中，在每个节点创建 限元平衡方程 ，其方程如下： H 是工程中最为关注的要素，它直接影响着水 了三个方程 ，两个是平衡方程 (位移)，一 利工程安全运行与效益的发挥。坝体结构的 个渗流连续方程(孔隙水压力)。同时求解 ∑[B】D【】B[]△【8]+∑[B】D【】{m】N 应力场和渗流场是相互影响的，单一的考 三个方程就能给出位移和孔隙水压力 的 虑渗流场和应力场并不能真实客观的反映 变化。 (△uw}_∑F 4() 实际渗流和应力状态。将应力场与渗流场 耦合分析中渗流方程可以用 Darcy方 [KI=[BIXDI[B】 在一个模型中同时进行计算分析是当前的 程表示如下： 发展方向，这更能突现坝区应力分布、坝体 + + ：0 (1) m }：f导，1，1，0。)) 坝基变形 以及孔隙水压力的变化，从而更 ax 。v 0t 确切地反映坝体应力渗流的真实情况。 式中，[B卜—一梯度矩阵 (也叫做应变矩 式中，k，kr——分别表示 x和Y方向 Geo—Studio软件具有全耦合分析功能， 阵)； 的渗流； 即在求解土体应力应变的时候也能同时求 [D卜—_{{E水本构矩阵； u——渗流流速； K【卜一 刚度矩阵； 解土体材料中孔隙水压力的变化情况，实 现应力渗流的耦合求解。软件中，应力应变 ^y——单位水的重量； I【廿—-耦合矩阵； 通过应力应变分析模块 SIGMM W模块求 0——体积含水量； △【8]——增量位移矢量； 解，孔隙水压力的变化通过渗流求解模块 t-_一 时间。 (△u }——增量孔隙水压力矢量。 SEEP／W求解。本文应用 Geo—Studio软件 对于弹性材料其体积含水量可用 1984