东北大学学报:自然科学版 ›› 2020, Vol. 41 ›› Issue (5): 635-641.DOI: 10.12068/j.issn.1005-3026.2020.05.005

• 资源与土木工程 • 上一篇    下一篇

寒区隧道围岩水热耦合数值分析

张泽, 王述红, 杨天娇, 张雨浓   

  1. (东北大学 资源与土木工程学院, 辽宁 沈阳110819)
  • 收稿日期:2019-04-10 修回日期:2019-04-10 出版日期:2020-05-15 发布日期:2020-05-15
  • 通讯作者: 张泽
  • 作者简介:张泽(1995-),男,山西大同人,东北大学博士研究生; 王述红(1969-),男,江苏泰州人,东北大学教授,博士生导师.
  • 基金资助:
    国家自然科学基金资助项目(51474050,U1602232); 中央高校基本科研业务费专项资金资助项目(N17010829); 辽宁省自然科学基金资助项目(20170540304,20170520341); 硅酸盐建筑材料国家重点实验室(武汉理工大学)开放基金资助项目(SYSJJ2017-08); 中建股份科技研发项目(CSCEC-2016-Z-20-8).

Numerical Analysis on Hydro Thermal Coupling of Surrounding Rocks in Cold Region Tunnels

ZHANG Ze, WANG Shu-hong, YANG Tian-jiao, ZHANG Yu-nong   

  1. School of Resources & Civil Engineering, Northeastern University, Shenyang 110819, China.
  • Received:2019-04-10 Revised:2019-04-10 Online:2020-05-15 Published:2020-05-15
  • Contact: WANG Shu-hong
  • About author:-
  • Supported by:
    -

摘要: 寒区隧道土体中的水分迁移和相变是冻害问题的主要诱因.基于混合物理论,建立考虑水分迁移和水冰相变的联合求解微分方程对水热耦合问题进行求解,并使用 COMSOL Multi-physics 软件进行模块开发,实现渗流-温度耦合数值模拟,进而将模拟结果与土柱冻结实验的结果进行对比,证明水热耦合模型是正确的.最后以西藏自治区米林隧道为例,对温度场、水分场模拟分析并对是否考虑水分迁移的温度场进行对比.结果表明:随着时间的增加,隧道顶部边界气温由-0.82℃降低到-9℃,隧道内部边界温度由-0.74℃ 降低到-11.11℃,并在3月份隧道温度回升;含冰量峰值出现在1月份,在3月份含冰量开始下降.同时,未考虑水分迁移的温度场中热传导速度较快,证明相变潜热对隧道中温度场的分布影响远大于液态水靠重力迁移造成的热对流传热.研究成果直观反映富水寒区隧道的冻害发生过程,具有一定的参考价值.

关键词: 寒区隧道, 水冰相变, 水分迁移, 水热耦合, COMSOL软件

Abstract: The mechanism of freezing damage of tunnels in the cold regions is complicated. These freezing damages are mostly related to water migration and phase transition in the soil. Based on the theory of mixture, a joint differential equation for the hydrothermal coupling problem of cold region tunnels considering water migration and water ice phase transition was established. The numerical simulation of temperature field and water field coupling was carried out by further development of COMSOL Multi-physics software, and the numerical simulation results were compared with the results of soil column freezing experiment to verify the effectiveness of the hydrothermal coupled numerical simulation model. Finally, taking the Milin Tunnel in Tibet Autonomous Region as an example, the temperature field and water field were simulated and analyzed, and the temperature field considering water migration was compared. The results show that with the increase of time, the temperature at the top boundary of the tunnel reduces from -0.82℃ to -9℃, the internal boundary temperature of the tunnel decreases from -0.74℃ to -11.11℃, and the tunnel temperature rises in March; the peak of ice content appears in January, and the ice content begins to decline in March. At the same time, the heat conduction velocity in the temperature field without considering moisture migration is faster, which proves that the influence of latent heat of phase transition on the distribution of temperature field in the tunnel is greater than that of liquid convection caused by gravity migration. The research results can reflect the occurrence process of freezing damage in the rich water tunnel in cold area, which has certain reference value.

Key words: tunnel in the cold regions, water ice phase change, moisture transfer, hydro thermal coupling, COMSOL software

中图分类号: