有效地模拟非饱和多孔介质中的固液气耦合效应在许多领域都具有重要意义。由于三相多孔介质全耦合方程的强非线性特性，必须采用有效的数值求解方法，如有限元法和有效的迭代算法。本文基于自适应松弛Picard方法，建立了一种分析多孔介质中固液气耦合作用的有效有限元方法。将耦合模型和有限元分析程序集成到计算机程序PorousH2M中，通过数值模拟与实验基准的比较，验证了所提方法的有效性。结果表明，自适应松弛Picard方法在效率和鲁棒性方面都优于传统方法，特别是在时间步长较大的情况下。与Picard格式相比，后者虽然能较好地避免非物理条件下的塑性变形现象。该方法为分析多孔介质中多相、多场耦合问题提供了重要的参考。版权所有©2014 John Wiley&Sons，Ltd。
