多孔介质中化学-热-水力-力学耦合分析和混合元方法.pdf

摘要 摘 要 近年来，非饱和多孔介质的热一水力一力学祸合行为的研究引起了工程界的广泛的 注意。人们己经投入了很大的精力来进行多孔介质瞬态THM祸合行为的定量研究。 本论文致力于这个课题的三方面的工作:1)建立了求解多孔介质中化学一热一水力 一力学祸合行为的数学模型;2)发展了这个数学模型初值、边值问题的混合元方法的数 值解法;3)建立多孔介质中的化学一热一水力一力学((CTHM)祸合本构模型。 这个数学模型由一组偏微分方程组成:包括孔隙水、孔隙气的质量守恒方程:混溶 于孔隙流体中的化学污染物的质量守恒方程;热量守恒方程和多孔介质混合物的总体动 量守恒方程。为了模拟混溶污染物在非饱和多孔介质中的传输，模型中包含了污染物传 输的六种控制机制:对流、分子弥散、机械逸散、吸附、降解、不动水效应。 采用有限元法求解数学模型的初值、边值问题。考虑到化学污染物质量守恒方程在 物理上的对流一扩散的特性和数学上的非自伴随和双曲线的特性，采用交错算法进行求 解，即:对化学污染物的质量守恒方程和其余的控制方程分别进行离散和求解。 化学污染物的质量守恒方程采用依赖于时间的对流一扩散方程的隐式特征线 Galerkin方法进行离散和求解，其中污染物的浓度作为基本未知量。该方法的基本思想 是对浓度在时域中进行针对物质粒子(Lagrangain)而不是空间点((Eulerain)的离散。 非饱和多孔介质中热一水力一力学((THM)祸合行为的有限元模拟通常要求解 u一P，一Pa-T混合形式的方程。混合方程的基本未知数u,P,P,,T分别为固相位移、 孔隙水压力、孔隙气压力和温度。位移场、水压力和气压力场及温度场分属于不同的函 数空间，因此应该对他们采用不同的有限元插值近似，以进行满足自然边界条件的祸合 控制方程弱形式的空间离散。 流体和固体力学u一P混合形式的有限元公式研究表明，u一P有限元插值的函数空 间必须满足Babuska-Brezzi条件或更简明的Zienkiewiez-Taylor分片试验。这些限制条 件判定了应用方便的位移(或速度)和压力采用等阶插值的低阶形函数的单元在不可压缩 条件下压力场会产生虚假的数值振荡。 与u一P形式的混合公式相比，非饱和多孔介质的u-P，一Pa-T公式更加复杂，求 解以u,P,P,,T为基本未知数的整体半离散方程、高斯积分点上祸合的非线性本构方程 迭代及计算一致性切线模量矩阵所需要的时间要长很多。因此就更加迫切需要发展低阶 高精度粗网格的有限元方法避免虚假的数值振荡，同时在满足数值计算精度基础上又能 降低计算费用。 本论文基于一点积分混合应变元在固体中的成功应用，提出了饱和、非饱和多孔介 质中THM祸合问题的稳定的一点积分混合应变元。本文提出的偏微分方程组的Galerkin 弱形式基于固体力学中的胡海昌·Washizu三变量广义变分原理，将之推广到由固体骨架 和不混溶的两种孔隙流体组成的三相非饱和多孔介质与温度场的相互作用的祸合场分 析中。 摘要 基于多孔介质连续理论基础上的热一水力一力学祸合行为的u,P.,P.,T形式偏微 分方程组的Galerkin有限元法弱形式，导出了混合元公式;在推导单元公式的过程中考 虑了固体骨架的材料和几何非线性。 人们己经认识到孔隙水溶液中特定化学物质会对土的水力一力学性质产生不利的 影喻。污染物的存在引起了孔隙水的电一化学性质变化，导致了土体微观结构和构造的 变化，产生导致土体体积、孔隙度和水力渗透性变化的化学应变。深入的了解化学对土 的性质的影响是工程土障、隧道开挖及钻井等岩土工程的设计和稳定性分析所必须的。 本论文在相关试验结果和化学一力学本构模型及热一水力一力学(THM)本构模型 的基础上，发展了一个多孔介质中的化学一热一水力一力学((CTHM)合本构模型。本 文的本构模型基于非饱和多孔介质CAP模型，考虑了热一化学对非饱和多孔介质的水 力一力学性质的影响。在本文的本构模型中引入化学软化函数以模拟孔隙水中污染物对 多孔介质的力学性质的影响。本文的CTHM本构模型描述了在多孔介质的化学一热一 水力一力学的祸合行为，在应力一吸力一温度四维空间中包含了五重屈服面。 对非饱和多孔介质的非关联塑性的多重屈服面的的屈服准则，针对本文提出的混合 元发