可压缩流体渗流流固耦合问题数值模拟的研究’.pdfVIP

维普资讯 第 21卷 第 5期 地 F 空 间 Vo1．21 No．5 2001年 12月 UNDERGROUND SPACE Dec．200l 文章编号 ：1001—811X(2001)一05—0386—06 可压缩流体渗流 的流固耦合 问题 数值模拟研究’ 薛 强， 梁 冰 ， 李宏艳 (辽宁工程技术大学力学与工程科学系，辽宁阜新 123000) 摘 要 ：可压缩流体在多孔介质 中的流动过程是流体流动和固体 变形相互 耦合作用的结果 ，研 究多孔介质 中流 固耦合渗流规律对于完善流 固耦合理论 以及 I 程应用有着重要意义 。本文在考虑 j流体 的可压缩性和 固体变形影响 的前提下 ，基于 流 固耦合力学理论 ，建立 j地下水渗流 问题 的流 固耦合数学模型 ，采用 了编制 FSC— SP程序 的对数学模型进行数值求解 。通过数值计算 ，比较 了考虑耦合效应 与非耦合 效应之 间的存在较大差异 ，这对于实际I程 中流 固耦合 问题 的数值模拟有着重要 的 理论意义和指导意义 。 关键词 ：可压缩流体 ；流 固耦合 ；渗流问题 }数值模拟 中圉分类号 ：TU452，TB12 文献标识码 ：A 引 言 目前，水文地质工程、环境工程、石油工程 以及采矿工程等众多领域常常涉及到可压缩流 体渗流 问题[1-8]。可压缩流体在多孔介质中的渗流是流体流动和固体变形相互作用的结果 。当 多孔介质变形时，其 内部孔隙空间也随之改变 ，引起储存在其中的流体 的体积发生改变 ，而孔 隙流体压力作为一种 内部载荷直接参与孔 隙的改变和多孔介质的变形 ，因此地下水在岩土中 的渗流 问题一方面是岩土的变形会引起流体渗透性能的改变 ，从而使流体 内的孔隙压力改变 ； 另一方面，流体的孔隙压力改变不仅直接改变岩土的应力状态 ，而且直接引起岩土的物理力学 性质的改变 ，两者是相互影响相互作用的。而在 以往的数值模拟过程 中均采用了纯渗流力学理 论进行研究，不考虑流体可压缩性和 固体骨架变形耦合作用及其影响。近几年来，随着流固耦 合理论和数值计算技术的发展 ，人们对可压缩流体渗流过程 中的力学机制有了更深入的了解 和认识 ：进行流固耦合分析是很有必要的。本文在考虑了流体的可压缩性和固体变形影响的前 提下，建立了可压缩流体渗流 问题的流 固耦合数学模型 ，这对进一步准确地模拟多孔介质中流 体的流动有着重要的意义。 · 收藕 日期 t2001—05—30 怍着简开 ：薛强(1976一)，男，博士 ，主要从事地下水污染及渗流方面的研究。 维普资讯 2001年第 5期 薛强等 ：可压缩流体渗流的流固耦舍 问题数值模拟研究 387 1 基 本理论 流俸在多孔介质中渗流问题属于流固耦合问题 ，不仅流体质点要在多孔介质中发生渗流， 而且岩土由于载荷作用要发生变形 ，岩土质点要发生位移而产生运动。由于以往渗流运动方程 中没有考虑流固耦合效应，因此必须建立流固耦合渗流运动方程来描述地下水中的渗流问题 ， 该方程包括应力场方程和渗流场方程。 建立方程前首先 引入几个假设 (1)岩土是各向同性的线弹性地质材料 及变形的多孔介质 。 (2)固体骨架的变形遵从 Terzaghi有效应力规律 ：％= 十户 (1) 式中q为总应力张量，矗为有效应力张量，P为孔隙水压力， 为Kronecker符号 ^ (3)土体内水渗流遵守达西定律，即 ql暑 (2) 式中q为水渗流