Page 100 - 水利学报2021年第52卷第5期

地质，2012，39（1）：42-46 .
［ 33］ 刘晓燕，赵军，石成，等 . 土壤恒温层温度及深度研究［J］. 太阳能学报，2007（5）：494-498 .

Frost heave modeling with thermos-hydro-mechanical coupling for lined canal
in coldregions considering solar radiation and its application

1
1
JIANG Haoyuan ，WANG Zhengzhong 1，2 ，LIU Quanhong ，
1 1 1
WANG Yi ，TAN Zhixiang ，MENG Xiaodong
（1. a. Cold and Arid Regions Water Engineering Safety Research Center，b. Key Laboratory of Agricultural Soil and Water Engineering
in Arid and Semiarid Areas of Ministry of Education，Northwest A&F University，Yangling 712100，China；
2. State Key Laboratory of Frozen Soil Engineering，Northwest Institute of Eco-Environment and Resources，Lanzhou 730000，China）

Abstract： The sun shading caused by the canal orientation leads to the sunny-shady slope effect， which
brings about the asymmetrical distribution of moisture， heat and deformation of a canal in cold regions. In
order to explore the influence of solar radiation on canal frost heave damage，a shadow calculation method
of canal was proposed and thermal boundary balance equation was established based on radiation algorithm，
using HOTTEL clear sky model. Next， a coupled thermo-hydro-mechanical model was proposed based on
the frost heave theory of frozen soil， considering the interaction of “environment-frozen soil-engineering” .
Based on the on-site monitoring data of large-scale canal going from "E-W" in cold regions， the simula⁃
tion method is verified and the coupled effect of temperature field，moisture field and deformation field and
the asymmetry & asynchrony development are systematically analyzed. The simulation results show that solar
radiation and shady slope shading cause obvious differences in the thermal boundary between the shady and
sunny slopes. The difference in daily average temperature can reach 3.5℃ ， causing the shady slope to
freeze 15 days earlier than the sunny slope and 1 day later to melt with the maximum difference of 31 cm
in freezing depth. The ice content of the shady slope during the freezing period is significantly higher than
that of the sunny slope，while the ice melting on the sunny slope during the melting period is faster，and
the moisture field is distributed in strips. The maximum frost heave deformation of the slope slab and bot⁃
tom slab occurs at 1/4~1/3 slope and near the sunny slope respectively， and the development of frost
heave deformation between the shady and sunny slopes is obvious with the maximum difference of 4.0 cm.
The above characteristics of asymmetric and asynchrony freezing all arise from the difference of temperature
gradient and freezing rate caused by space effect of solar radiation and time effect of diurnal temperature
difference. Therefore，reasonable suggestions for the selection of the thermal boundary of the frost heave nu⁃
merical simulation are given. The results can provide a theoretical basis for the design of anti-frost heave，
review and repair of the large-scale canal in cold-arid regions.
Keywords： a cold-region canal； solar radiation； thermo-hydro-mechanics coupling of frozen soil； asymme⁃
try development；thermal boundary；freezing rate

（责任编辑：韩 昆）

— 602 —