基于刚性冰模型,应用有限差分法离散方程组,对于饱和颗粒土开放系统的冻结过程进行了一维数值模拟,给出了冻结缘内参数的分布.计算过程中修正了刚性冰模型中分凝冰产生条件,提出以冻结缘内冰压与载荷的关系作为分凝冰产生的判据,计算得到的冻胀量与实验室冻胀实验的测量数据较为吻合.将计算结果与现场实测资料进行了比较,温度场、含水量分布等结果在分布规律方面与现场实测资料相符.