摘要
建立了渗流作用下裂隙岩体双孔冻结的数值计算模型,采用COMSOL Multiphysics有限元分析软件,实现了裂隙岩体冻结过程中耦合及相变问题的模拟,分析了裂隙存在对岩体冻结温度场及渗流场的影响,研究了不同渗流速度下对应的冻结温度场的分布规律。结果表明:裂隙的存在会减薄上游冻结壁厚度,等温线在冻结壁外圈形成"心形"分布;渗流速度对冻结壁的发育及其交圈时间影响较大,渗流速度愈大,冻结壁发育愈慢,交圈时间愈长。在本文条件下,当渗流速度小于7 m/d时,流速与交圈时间基本呈线性关系;当渗流速度超过7 m/d时,交圈时间呈现急剧增长,达到10 m/d时,交圈时间趋于无限大,即交不上圈。