二维隧道裂隙岩体冻胀相场法仿真：从理论到完整实现
摘要
寒区隧道围岩中的裂隙水在低温条件下发生冻结相变，产生约9%的体积膨胀，由此引发的冻胀力常导致裂隙扩展乃至岩体渐进性破坏。本文基于相场断裂理论，构建了热-力耦合的裂隙岩体冻胀扩展数值模型，实现了温度场演化、冻胀力驱动开裂、新裂隙自发萌生以及多裂隙网络扩展的全过程仿真。模型采用等效热膨胀系数法施加裂隙冰相变膨胀，通过相场变量d∈[0,1]弥散化表征裂隙几何，自动处理裂隙萌生、扩展、分叉与合并等复杂拓扑变化。本文提供了完整的二维有限元实现代码（基于FEniCS/dolfinx平台），包含Python代码及详细的逐段解释，涵盖网格生成、温度场求解、力学平衡、相场演化、冻胀力加载及自适应时间步进等核心模块。最后通过典型算例验证了模型捕捉新裂隙萌生和多裂隙网络扩展的能力。
1 引言
1.1 研究背景
寒区隧道工程面临严峻的冻害问题。当围岩温度降至冰点以下时，赋存于裂隙中的水发生相变，体积膨胀约9%，在受限空间内产生巨大的冻胀力。当该力超过岩体抗拉强度时，裂隙即发生萌生与扩展。冻融循环的反复作用将进一步导致裂隙网络的渐进性发展，最终威胁隧道衬砌结构的安全性与耐久性。
1.2 现有方法的局限性
传统数值方法在模拟裂隙冻胀扩展时面临诸多挑战：
- 离散裂隙方法（如XFEM、DEM）需要显式追踪裂隙界面，难以处理复杂裂隙网络的萌生与拓扑变化；
