1. 引言
季冻区路基每年都发生季节性融化和冻结,由此会产生一系列特殊的工程地质问题。季冻区路基冻胀会引起道路变形、裂缝;暖季冻土融化,导致道路出现翻浆冒泥、路基下沉、路肩滑蹋和路面沉陷等病害。路基温度场与环境因素、路基结构型式、路基填筑材料及路基湿度状况相关,影响因素较为复杂。为了了解季冻区路基的冻胀机理及发展规律、控制冻土路基工程中的冻害、确保冻土区道路工程的稳定性、耐久性,必须深入研究路基冻融过程中温度场的变化规律,为季冻区路基的冻害处理提供理论依据。
刘志强分析了影响冻土路基温度场的主要因素,建立了冻土路基模拟试验装置 [1] 。郑仕建对特殊环境下路面、路基的温度场进行了模拟研究,得出了不同参数条件下的温度场分布规律 [2] 。毛雪松建立了水泥混凝土路面结构的路基温度场模型,分别在第一类边界条件(附面层理论)和第二、第三类边界条件下对路基温度场进行数值模拟 [3] 。令锋运用传热学理论和渗流理论导出了路堤侧向有地表积水渗入情形下,多年冻土区路基温度场渗流场藕合作用的控制微分方程,对路基温度场的变化状态进行了数值分析 [4] 。葛建军采用数值模拟方法分析了青藏铁路多年冻土区路基保温护道的效果 [5] 。邹泽雄结合实际监测数据,利用ANSYS软件对路基温度场进行有限元数值模拟 [6] 。汪海年对青藏公路五道梁地区路基温度场进行了有限元分析,路线走向对冻土路基温度场的对称性有着重要影响 [7] 。汤涛运用ABAQUS及其二次开发平台,建立了多年冻土地区路基温度场有限元分析模型,运用该模型对低等级公路窄幅路基和高等级公路宽幅路基温度场进行对比分析 [8] 。白青波建立了一维路面路基结构温度场模拟模型,并通过路基温度场反演分析,验证了数值模型的正确性 [9] 。岳祖润分析了天然地表及路基不同位置的地温分布规律,建立温度场的仿真模型,研究温度沿深度方向的变化规律 [10] 。王威娜结合季节冻土地区实际气候条件,基于第II类与第III类混合热边界条件,模拟季节冻土地区路基温度场变化规律,分析阴阳坡面温度差异 [11] 。汪双杰针对典型工况条件下冻土路基内温度场的分布特征及其随时间的变化规律,分析了路面宽度、风速、路基高度和走向等因素对冻土路基温度场的影响规律 [12] 。韩铠屹运用有限元软件模拟路基冬季温度场,分析路基温度场在冬季各月份的分布情况及不同年份下路基温度场的变化 [13] 。许顺国运用有限元分析,对多年冻土区草皮护坡路基温度场进行了数值模拟 [14] 。张尧禹分析了多年冻土区路基温度场与外界环境的响应趋势,得出了路基的修筑对天然冻土上线的影响程度 [15] 。
本文结合黑龙江省北安至古城高速公路所在地的气候特点,建立了路基温度场分析模型,分析了路基温度场的分布规律及时变特性,为路基冻胀病害的产生提供了理论依据。
2. 寒区路基温度场数值模型建立
2.1. 有限元模型
黑龙江省北安至古城高速公路,双向四车道,路基宽度为26 m,路基边坡坡率为1:1.5。建立的二维温度场分析模型尺寸如图1所示。设置4个路径:竖向路径1 (SLJ1),竖向路径2 (SLJ2),横向路径1 (HLJ1),横向路径2 (HLJ2)。在考虑当地冻深的条件下,选取路基计算深度为3 m。
采用ANSYS软件进行瞬态温度场分析,单元类型为SOLID42,划分单元数量为6870个。网格划分结果如图2所示。
研究表明,土基底面的温度不论是年温度变化或是日温度变化均较稳定。所以,在进行有限元热分析时土基底部温度视为常温,取值大小根据路面厚度及当时气温决定。瞬态热分析中,需给定路表近地空气温度值和对流系数以及模型底面的地层温度。
2.2. 温度场分析参数
参数选取如下:路基土的密度为1800 kg/m3,焓值取值为980 J。传热系数取值如表1所示,比热取值如表2所示。
2.3. 气象资料
项目所在地为克东地区,该地区全年各月的温度数据如图3所示。全年温度荷载加载时为瞬态分析,主要分析关键点的温度随时间变化的时变规律特性,分为12个荷载子步,在ANSYS软件中由表格数据的方式加载。
Figure 1. Two-dimensional temperature field analysis model
图1. 二维温度场分析模型
表1. 传热系数
表2. 比热
Figure 3. The temperature of the whole year
图3. 全年各月温度
3. 路基温度场分布规律
3.1. 温度场竖向分布规律
一月份为项目所在地的全年最低温度,选取一月份每天的最高温度、最低温度及平均温度分别进行数值模拟,得到路基内部的温度分布规律。根据模拟结果,得到路基内部竖向路径1、竖向路径2 (SLJ1、SLJ2)在三个不同工况(高温、低温、平均温度)条件下温度的分布如图4、图5所示。
从图4可以看出,在竖向路径1位置,随着深度的增加,路基内部的温度升高,在路表处,路基的温度接近环境温度。在不同的环境温度下,路基深度较浅时温度差别较大,随着路基深度的增加,路基温度的差值逐渐减小。同时从图4也可以看出,竖向路径1处的大地冻深大致为2.5 m。
从图5可以看出,在竖向路径2位置,温度沿深度的分布与路径1的分布规律完全相同。与路径1的主要区别是路径2处的大地冻深约为2.0 m。即路径2处的冻深明显小于路径1处的冻深。这主要是由于路径的位置不同造成的。由图1可知,路径1接近路基两侧边坡的位置,路基内部的温度受到环境温度的影响更大,路基顶面及边坡坡面都与路基存在着热交换;而竖向路径2离边坡较远,其热交换主要来自路基顶面,所以其冻深要小于路径1的冻深。综合以上分析,在不考虑边坡的影响情况下,当地的路基冻深应当为2.0 m,这与当地的实测大地冻深基本一致。说明数值模拟结果较为准确。
3.2. 温度场横向分布规律
通过路基温度场的数值模拟,得到一月份路基横向路径1 (HLJ1)和横向路径2 (HLJ2)在三个不同工况(高温、低温、平均温度)条件下温度的分布状况如图6及图7所示。
通过图6及图7可以看出,路基横向的温度分布与路基竖向的温度分布相类似,从边坡至路基中心温度逐渐升高。环境温度对路基横向的影响范围大致为2.0 m,即在距离边坡表面2.0 m以内,路基温度随横向距离的增加而升高,当距离边坡表面的距离超过2.0 m以后,路基内部的横向温度接近一致。环境温度对路基横向温度的影响范围与大地冻深基本一致。
(a) 1月10日 (b) 1月20日(c) 1月30日
Figure 4. The temperature distribution of vertical path 1
图4. 竖向路径1的温度分布
(a) 1月10日 (b) 1月20日(c) 1月30日
Figure 5. The temperature distribution of vertical path 2
图5. 竖向路径2的温度分布
(a) 1月10日 (b) 1月30日
Figure 6. The temperature distribution of crosswise path 1
图6. 横向路径1的温度分布
(a) 1月10日 (b) 1月20日(c) 1月30日
Figure 7. The temperature distribution of crosswise path 2
图7. 横向路径2的温度分布
(a) A点 (b) B点(c) C点
Figure 8. The time curve of temperature
图8. 温度的时变曲线
(a) 最高温度 (b) 最低温度(c) 平均温度
Figure 9. The comparison of ABC point temperature change
图9. A/B/C温度变化对比
4. 温度场时变特性分析
根据项目所在地一年内气温的变化情况,分别计算环境最高温度、最低温度、平均温度三个工况下路基内部的温度场变化情况。图1所示的A、B、C三个点的温度随时间变化如图8所示(A点接近于路表,B,C点处于路基内部)。
按全年环境最高温度、最低温度、平均温度三个工况计算结果得到路基内部A、B、C三点的温度对比如图9所示。
由图8可以看出,路基内部不同位置处,温度随时间的变化总体趋势基本相同,并且可以看出路基内部不同位置处,全年温度变化趋势与环境温度变化趋势一致。由于A、B、C三点所处位置的不同,导致不同点之间的温度变化存在一定的差别。
由图9可以看出,点A处的温度对环境温度的变化更为敏感,而点B与点C的温度变化基本一致,且其对环境温度的变化的敏感程度要小于A点。
由于A点处于路表位置,所在位置的温度受环境温度影响最大,而B、C点处于路表以下,受环境温度变化的影响相对较小。如果按全年的最低气温来计算,A点处的温度与环境温度的变化基本一致,每年大致有6个月的时间温度处于零度以下(10月至第二年的4月),而B、C点全年的温度都处在零度以下。如果按全年的最高气温来计算,A点处的温度每年大致有3个月的时间处于零度下,而B、C点的温度全年有5个月处于零度以下。路基内部处于多年冻结或季节性冻结状态,路基容易出现冻胀和融沉病害。
5. 研究结论
结合ANSYS软件建立了路基温度场的分析模型,分析了路基内部温度场的分布规律及时变特性。
路基内部竖向温度场分布规律为随深度的增加温度升高,一月份最冷时的路基冻深约为2.0 m;路基横向的温度场分布特征为,环境温度对边坡横向的影响范围约为2.0 m,与路基冻深基本一致。
路基内部不同位置处全年的温度变化与环境温度变化趋势基本一致。路表处的温度与环境温度变化基本相同,路基内部温度变化受环境温度变化影响较小,路基深处土体处于多年冻结或季节性冻结状态。