基于计算流体力学和岩土力学,考虑渗流场和应力场的相互作用,推导介质的平衡方程和饱和孔隙流体的连续性方程,建立流–固耦合模型。结合Galerkin法,对平衡方程进行有限元离散,同时,采用向后差分格式和二阶精度Warming-Beam差分格式,对连续性方程进行有限差分离散,从而对模型进行数值求解。最后,以一富水双向隧道为例,进行施工数值仿真,得到不同施工阶段隧道围岩的孔隙水压力、位移场、应力场的分布及其支护锚杆的受力特性,并将计算结果与相关现场监测得到的一般性规律进行对比,为类似隧道工程的设计和施工提供参考。