吉林大学学报(工学版) ›› 2026, Vol. 56 ›› Issue (9): 2395-2407.doi: 10.13229/j.cnki.jdxbgxb.20250212

• 交通运输工程·土木工程 • 上一篇    

寒区路基水热耦合理论模型及变形规律

何永明(),张鹏飞,王锦扬,佟鑫,梁晨   

  1. 东北林业大学 土木与交通学院,哈尔滨 150040
  • 收稿日期:2025-02-20 出版日期:2026-09-01 发布日期:2026-09-07
  • 作者简介:何永明(1979-),男,副教授,博士.研究方向:超高速公路. E-mail: hymjob@nefu.edu.cn
  • 基金资助:
    长沙理工大学极端环境绿色长寿道路工程全国重点实验室开放基金项目(kfj230105);黑龙江省自然科学基金项目(LH2023E011)

Theoretical model of water-heat coupling and deformation law of roadbed in cold area

Yong-ming HE(),Peng-fei ZHANG,Jin-yang WANG,Xin TONG,Chen LIANG   

  1. School of Civil and Transportation Engineering,Northeast Forestry University,Harbin 150040,China
  • Received:2025-02-20 Online:2026-09-01 Published:2026-09-07

摘要:

为探究寒区路基在水?热耦合作用下的变形机理,以鹤哈高速(G1111)苔青至伊春段(62.44 km)为对象,构建了二维水热迁移耦合模型,结合现场监测与数值模拟方法,揭示了冻融循环过程中路基温度场、水分场与变形场的相互作用规律。基于相变传热与非等温水分迁移理论,采用Abaqus软件进行仿真分析,在路基横断面上自左侧路肩至右侧路肩按1.25 m间距均匀布设9个综合监测点,量化路基表层冻胀与融沉阶段的变形特征。结果表明:冻结期模拟路段路基表层最大冻胀位移可达8.48 mm,横截面变形呈抛物线分布;春季融沉速率(0.12 mm/d)显著高于冻胀期(0.07 mm/d),且11月至次年1月的冻胀速率较1~2月高18.7%;竖向位移极值与地温梯度呈负相关。研究成果为寒区道路抗冻设计与病害防控提供了理论依据。

关键词: 水?热耦合, 变形规律, 数值模拟, 冻土路基, 有限元

Abstract:

To investigate the deformation mechanisms of subgrade in cold regions under hydro-thermal coupling effects, this study focuses on the TQ-Yichun section (62.44 km) of the G1111 He-ha Expressway. A two-dimensional coupled hydro-thermal migration model was developed to analyze the interactions among temperature fields, moisture fields, and deformation fields during freeze-thaw cycles. Based on the theory of phase change heat transfer and non-isothermal moisture migration, Abaqus software was used for simulation and analysis, and nine integrated monitoring points were deployed along the roadway cross-section at a spacing of 1.25 m to quantify the deformation characteristics of the surface layer of the roadway in the freezing and thawing stages. Key findings include: The maximum frost heave displacement at the subgrade surface during freezing periods reached 8.48 mm, with cross-sectional deformation exhibiting a parabolic distribution. The thaw settlement rate in spring (March-May, 0.12 mm/day) was significantly higher than the frost heave rate during freezing periods (0.07 mm/day), with frost heave rates between November-January surpassing those of January-February by 18.7%. Vertical deformation extremes showed a negative correlation with geothermal gradients. By introducing dynamic ice-water phase-change correction coefficients, the study improves simulation accuracy and proposes a predictive model for subgrade deformation. These results provide critical insights for optimizing frost-resistant designs in cold-region road engineering, particularly in determining interfacial heat transfer coefficients and insulation layer thickness. The research establishes a theoretical foundation for mitigating subgrade deformation and controlling freeze-thaw-induced pavement deterioration in permafrost environments.

Key words: hydro-thermal coupling, deformation law, numerical simulation, permafrost roadbed, finite element

中图分类号: 

  • U416.1

图1

计算区域结构"

图2

计算区域网格划分"

表1

土体的热物理参数"

土壤类型土体状态

导热系数k/

[W·(m·℃)-1

比热容C/

[J·(kg·℃)-1

亚黏土冻土1.7471 864
未冻土1.2111 574
砂砾土冻土2.5861 087
未冻土1.9761 516

表2

土体的水分运动参数"

土壤参数单位亚黏土砂砾土
导水率kucm/s1×10-52×10-6
扩散率Ducm2/s11×10-5

表3

伊春市2023年气候数据"

月份

月平均

气温

/℃

月最高

气温

/℃

月最低

气温

/℃

日太阳辐射量/[(MJ·m-2·d-1)]有效日照时间/(h·d-1
1-22.242.30-41.305.424.41
2-13.645.80-29.508.416.96
3-1.2518.40-19.6012.327.21
45.3720.90-5.8014.245.92
513.4928.90-2.8017.698.22
618.4534.104.6015.736.89
721.2230.4014.1012.354.87
820.2831.2011.3010.904.76
914.5229.400.1011.497.24
106.1820.90-8.208.706.59
11-9.783.90-26.105.083.31
12-19.503.50-35.204.292.97

表4

计算模型各边界处风速表 (m/s)"

月份阴坡风速阳坡风速路基风速地基风速
12.422.672.842.72
22.783.013.083.06
33.293.373.523.44
43.844.034.204.16
53.783.903.983.94
63.083.113.213.15
72.602.672.812.73
82.542.612.732.68
92.792.883.012.92
103.213.273.443.38
113.083.193.363.27
122.632.722.832.79

图3

融化-冻结阶段路基温度场分布"

图4

冻结-融化阶段路基温度场分布"

图5

一年内A点沿深度方向的温度趋势图"

图6

冻土路基水分场分布"

图7

路基各月份纵向水分饱和度分布图"

图8

11月份的冻土路基位移场分布"

图9

路基表面竖向位移图"

图10

路基阴阳坡不均匀沉降值"

表5

路面平均竖向相对位移极值"

路面宽

度/m

路面平均竖向相对变形量/mm
最大值最小值极差
0.004.800.853.95
1.252.23-0.883.11
2.500.86-1.642.50
3.751.19-1.352.54
5.001.37-1.763.13
6.251.23-1.272.50
7.501.15-1.352.50
8.753.02-0.743.76
10.008.480.787.70

图11

道路横截面位移场极值分布概率图"

表6

道路横截面竖向位移场模型反演"

月份线性拟合抛物线拟合
公式R2公式R2
1y=0.223d+3.4680.070 3y=-2.388d2+0.261d+7.2760.842 4
2y=0.219d+3.3430.064 4y=-2.487d2+0.271d+7.2890.851 4
3y=0.214d+2.8310.113 1y=-1.761d2+0.198d+5.7110.887 4
4y=0.152d+1.9030.082 2y=-1.272d2+0.142d+3.9790.658 5
5y=0.091d+3.2190.026 5y=-1.689d2+0.178d+5.8160.833 9
6y=0.061d+3.6290.009 9y=-1.921d2+0.198d+6.5190.860 6
7y=0.067d+4.0460.012 4y=-1.865d2+0.193d+6.8630.836 1
8y=-0.032d+3.1500.099 8y=-0.372d2+0.034d+3.6450.990 4
9y=-0.118d+3.4060.151 7y=-1.087d2+0.097d+4.8190.969 9
10y=0.084d+3.0140.113 2y=-0.365d2+0.045d+3.6680.374 9
11y=0.181d+2.3280.130 3y=-1.162d2+0.134d+4.2860.709 9
12y=0.231d+2.5960.107 3y=-1.819d2+0.205d+5.5860.783 6

图12

地表温度模型验证"

图13

路基中心沉降模型验证"

[1] 陈佩佩, 尚智, 王熳祺, 等. 寒区高速铁路冻土路基的水热耦合响应特性分析[J]. 水文地质工程地质, 2024, 52(3): 152-162.
Chen Pei-pei, Shang Zhi, Wang Firm-qi, et al. Characteristic analysis of hydrothermal coupling response of frozen soil roadbed for high-speed railway in cold region[J]. Hydrogeology Engineering Geology, 2024, 52(3): 152-162.
[2] 彭勇. 路堑边坡防护设计方案[J]. 中国科技信息, 2024(14): 93-95.
Peng Yong. Design scheme of road graben slope protection[J]. China Science and Technology Information, 2024(14): 93-95.
[3] 于景飞, 张文月, 李明浩. 季节性冻土区铁路路基冻融变形规律[J]. 长江科学院院报, 2020, 37(5):133-138.
Yu Jing-fei, Zhang Wen-yue, Li Ming-hao. Freezing and thawing deformation law of railway roadbed in seasonal permafrost area[J]. Journal of Changjiang Academy of Sciences, 2020, 37(5): 133-138.
[4] 李晓芳, 官少龙, 滕伟玲. 基于高山草甸区季冻土的冻融变形特性影响因素研究[J]. 公路,2018, 63(4): 47-51.
Li Xiao-fang, Guan Shao-long, Teng Wei-ling. Research on factors affecting freeze-thaw deformation characteristics based on seasonal permafrost in alpine meadow area[J]. Highway, 2018, 63(4): 47-51.
[5] 李万志, 马海玲, 庞昕玮, 等. 气候变暖背景下青海三江源区季节冻土冻融特征研究[J]. 冰川冻土, 2023, 45(4): 1233-1241.
Li Wan-zhi, Ma Hai-ling, Pang Xin-wei, et al. Study on permafrost and thawing characteristics of seasonally frozen soil in the San-jiang-yuan area of Qinghai under the background of climate warming[J]. Glacial Permafrost, 2023, 45(4): 1233-1241.
[6] 赫金良, 李晓宁, 李青龙, 等. 基于优化SVM的粗粒土路基冻胀特性影响因素敏感性分析[J]. 科学技术与工程, 2024, 24(22): 9577-9586.
He Jin-liang, Li Xiao-ning, Li Qing-long, et al. Sensitivity analysis of factors influencing the freezing and expansion characteristics of coarse-grained soil subgrade based on optimised SVM[J]. Science, Technology and Engineering, 2024, 24(22): 9577-9586.
[7] 张玉芝, 杨威, 李锐, 等. 季冻区重载铁路路基冻胀融沉监测及规律分析[J]. 中国铁道科学, 2022, 43(5): 1-10.
Zhang Yu-zhi, Yang Wei, Li Rui, et al. Monitoring and law analysis of freezing and thawing subsidence of heavy railway roadbase in monsoon frozen area[J]. China Railway Science, 2022, 43(5): 1-10.
[8] 甘海龙, 唐先习, 李小博. 季冻区高速铁路路基冻融变形特征[J]. 科学技术与工程, 2022, 22(9): 3760-3767.
Gan Hai-long, Tang Xian-xi, Li Xiao-bo. Freezing and thawing deformation characteristics of high-speed railway roadbase in monsoon frozen area[J]. Science Technology and Engineering, 2022, 22(9): 3760-3767.
[9] Lu B, Zhao W, Li S, et al. Study on seasonal permafrost roadbed deformation based on water-heat coupling characteristics[J]. Buildings, 2024, 14(9): No.2710.
[10] Su A, Hai M, Wang M, et al. Analytical study on water and heat coupling process of black soil roadbed slope in seasonal frozen soil region[J]. Sustainability, 2024, 16(19): No.8427.
[11] 张哲, 付伟, 张军辉, 等. 循环荷载下冻融路基黏土长期塑性行为[J]. 吉林大学学报: 工学版, 2023, 53(6): 1790-1798.
Zhang Zhe, Fu Wei, Zhang Jun-hui, et al. Long-term plastic behavior of frozen-thawed roadbed clay under cyclic loading[J]. Journal of Jilin University(Engineering and Technology Edition), 2023, 53(6): 1790-1798.
[12] 侯雅, 李伟平, 左金清. 土壤冻融过程中的水热参数化方案研究进展[J]. 高原气象, 2024, 44(1): 1-15.
Hou Ya, Li Wei-ping, Zuo Jin-qing. Progress of hydrothermal parameterisation scheme during soil freezing and thawing[J]. Plateau Meteorology, 2024, 44(1): 1-15.
[13] 赵亮, 景立平, 单振东. 冻土冻胀模型研究现状与进展[J]. 自然灾害学报, 2020, 29(4): 43-52.
Zhao Liang, Jing Li-ping, Shan Zhen-dong. Current status and progress of research on permafrost frost expansion modelling[J]. Journal of Natural Hazards, 2020, 29(4): 43-52.
[14] 宋存牛. 冻融过程中土体水热力耦合作用理论和模型研究进展[J]. 冰川冻土, 2010, 32(5): 982-988.
Song Cun-niu. Progress in theory and modelling of soil hydrothermal coupling action during freeze-thaw process[J]. Glacial Permafrost, 2010, 32(5): 982-988.
[15] Zhan Y, Lu Z, Yao H. Numerical analysis of thermo-hydro-mechanical coupling of diversion tunnels in a seasonally frozen region[J]. Journal of Cold Regions Engineering, 2020, 34(3): No.04020018.
[16] 王自力, 李金峰, 滕继东, 等. 基于未冻水膜压力判据的冻土水热力耦合冻胀模型研究[J]. 岩土工程学报, 2023, 45(5): 997-1007.
Wang Zi-li, Li Jin-feng, Teng Ji-dong, et al. Study on hydrothermal coupled freezing model of frozen soil based on unfrozen water film pressure criterion[J]. Journal of Geotechnical Engineering, 2023, 45(5): 997-1007.
[17] 李长雨, 马桂霞, 郝光, 等. 季节冻土地区路基冷阻层温度场效应[J]. 吉林大学学报: 地球科学版, 2018, 48(4): 1174-1181.
Li Chang-yu, Ma Gui-xia, Hao Guang, et al. Temperature field effect of cold resistance layer of roadbed in seasonal frozen soil area[J]. Journal of Jilin University(Earth Science Edition), 2018, 48(4): 1174-1181.
[18] 苗祺, 牛富俊, 林战举, 等. 季节冻土区高铁路基冻胀研究进展及展望[J]. 冰川冻土, 2019, 41(3):669-679.
Miao Qi, Niu Fu-jun, Lin Zhan-ju, et al. Progress and prospect of research on freezing and expansion of high-speed railway roadbed in seasonal frozen soil zone[J]. Glacial Permafrost, 2019, 41(3): 669-679.
[19] 毛雪松, 李宁, 王秉纲, 等. 多年冻土路基水-热-力耦合理论模型及数值模拟[J]. 长安大学学报: 自然科学版, 2006(4): 16-19, 62.
Mao Xue-song, Li Ning, Wang Bing-gang, et al. Theoretical model and numerical simulation of water-heat-force coupling of multi-year frozen soil roadbed[J]. Journal of Chang'an University (Natural Science Edition), 2006(4): 16-19, 62.
[20] 赵刚. 寒区高等级公路路堑边坡春季浅层滑塌机理研究[D]. 哈尔滨: 哈尔滨工业大学交通科学与工程学院, 2010.
Zhao Gang. Study on the mechanism of shallow landslide of high-grade highway cutting slope in cold region in spring[D]. Harbin: School of Transportation Science and Engineering, Harbin Institute of Technology, 2010.
[21] 毛雪松, 王秉纲, 胡长顺, 等. 冻土路基水热迁移问题的理论模型及数值模拟[J]. 中外公路, 2006(1): 23-26.
Mao Xue-song, Wang Bing-gang, Hu Chang-shun, et al. Theoretical modelling and numerical simulation of water-heat migration problem in permafrost roadbeds[J]. Chinese and Foreign Highway, 2006(1): 23-26.
[22] 龚强, 晁华, 朱玲, 等. 东北地区地温和冻结深度时空特征的细化分析[J]. 冰川冻土, 2021, 43(6): 1782-1793.
Gong Qiang, Chao Hua, Zhu Ling, et al. Refinement analysis of spatial and temporal characteristics of ground temperature and freezing depth in Northeast China[J]. Glacial Permafrost, 2021, 43(6): 1782-1793.
[23] 杨科. 多年冻土路基温度场变化及数值模拟[J]. 中国科技信息, 2022(14): 81-83.
Yang Ke. Variation and numerical simulation of permafrost subgrade temperature field[J]. China Science and Technology Information, 2022(14): 81-83.
[24] 柳波. 多年冻土地区公路路基不稳定变形研究[D]. 西安: 长安大学公路学院, 2015.
Liu Bo. Study on unstable deformation of highway subgrade in permafrost area[D]. Xi'an: School of Highway, Chang-an University, 2015.
[25] 汪海年, 侯仲杰, 杨玉杰. 多年冻土路基对流换热试验方法[J]. 公路交通科技, 2009, 26(3): 43-48.
Wang Hai-nian, Hou Zhong-jie, Yang Yu-jie. Test method of convective heat transfer in roadbeds with multiyear frozen soil[J]. Highway and Transportation Science and Technology, 2009, 26(3): 43-48.
[26] 李方昊, 范昊明, 石昊, 等. 东北地区极端大风时空分布特征及其对风蚀量的影响[J]. 应用生态学报, 2024, 35(1): 87-94.
Li Fang-hao, Fan Hao-ming, Shi Hao, et al. Temporal and spatial distribution characteristics of extreme high winds in Northeast China and their effects on wind erosion[J]. Journal of Applied Ecology, 2024, 35(1): 87-94.
[27] 毛云程, 李国玉, 张青龙, 等. 季节冻土区黄土路基水分与温度变化规律研究[J]. 冰川冻土, 2014, 36(4): 1011-1016.
Mao Yun-cheng, Li Guo-yu, Zhang Qing-long, et al. Study on the variation of moisture and temperature of loess subgrade in seasonal permafrost area[J]. Journal of Glaciology and Geocryology, 2014, 36(4): 1011-1016.
[28] 徐泉涌. 某公路岛状多年冻土路基处理[J]. 低温建筑技术, 2005(3): 91-92.
Xu Quan-yong. Treatment of island-shaped perennial frozen soil roadbed for a highway[J]. Low Temperature Construction Technology, 2005(3): 91-92.
[29] 纪伦, 张红菊, 郭宏斌, 等. 北方寒冷地区温度场作用下公路路表形变规律分析[J]. 哈尔滨工业大学学报, 2024, 56(9): 113-123.
Ji Lun, Zhang Hong-ju, Guo Hong-bin, et al. Analysis of surface deformation of highway road surface under the action of temperature field in cold areas of northern China[J]. Journal of Harbin Institute of Technology, 2024, 56(9): 113-123.
[1] 王国强,马超洋,郭志远,王帅,左国棚,刘世博. 液压挖掘机履带链轨节疲劳分析[J]. 吉林大学学报(工学版), 2026, 56(9): 2292-2301.
[2] 赵勇,杨佳乐,金文明,郑祺峰. 减速器壳体轴承座胀断加工断裂本质及参数[J]. 吉林大学学报(工学版), 2026, 56(3): 681-688.
[3] 鲁亮,颜浩天. 基于耐震时程法的耗能自复位铰节点RC框架抗震性能[J]. 吉林大学学报(工学版), 2026, 56(3): 700-710.
[4] 李义,刘轶,姚卫国. 汽车内饰件胶接过程中产品质量的提升[J]. 吉林大学学报(工学版), 2025, 55(9): 2926-2934.
[5] 于征磊,张超磊,陈立新,胡平,徐涛,郭滨恺. 高速公路声屏障板仿生结构设计及声学力学性能[J]. 吉林大学学报(工学版), 2025, 55(9): 3079-3088.
[6] 李欢,刘千喜,张长鑫,张健. 大功率超声波焊接纯铜的动态摩擦及超声软化过程[J]. 吉林大学学报(工学版), 2025, 55(8): 2548-2554.
[7] 范亮,曾文,文强,赵富裕,徐英铭. 集束群钉装配式钢-混组合梁桥自振特性与车桥耦合分析[J]. 吉林大学学报(工学版), 2025, 55(7): 2354-2364.
[8] 姜歌东,王昊,荆亚彬. 接触热阻对高速滚珠丝杠副温升特性的影响[J]. 吉林大学学报(工学版), 2025, 55(6): 1915-1922.
[9] 赵勇,金文明,郑祺峰,寇淑清. 裂解槽深度对减速器壳体轴承座裂解性能的影响[J]. 吉林大学学报(工学版), 2025, 55(5): 1552-1558.
[10] 卢荡,王晓凡,吴海东. TWEEL轮胎接地压力均布特性分析[J]. 吉林大学学报(工学版), 2025, 55(3): 811-819.
[11] 周正峰,唐虎城,欧信旺. 基于黏聚区模型的沥青混凝土SCB试验模拟分析[J]. 吉林大学学报(工学版), 2025, 55(12): 3955-3963.
[12] 马天忠,杨嘉俊,王正振,陈璋佳,郭保文. 浸水状态下黄土地区长短桩基础承载特性试验[J]. 吉林大学学报(工学版), 2025, 55(12): 3942-3954.
[13] 吴金星,马宇翔,肖嘉邦,徐耀,李松歌. H型翅片椭圆管壁面磨损特性的数值模拟[J]. 吉林大学学报(工学版), 2024, 54(9): 2495-2501.
[14] 刘化民,杨舒涵,李义,梁策,韩奇钢. 推力杆球铰仿生表面改进及有限元分析[J]. 吉林大学学报(工学版), 2024, 54(9): 2733-2740.
[15] 梁策,李敏,李义,梁继才,韩奇钢. 轿车前轴摇臂衬套仿生柔性接触表面摩擦特性数值模拟[J]. 吉林大学学报(工学版), 2024, 54(8): 2181-2186.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!