基于ANSYS平台冻土路基温度场数值模拟.docVIP

基于ANSYS平台冻土路基温度场数值模拟 　　摘要：应用有限元软件ANSYS热分析模块对冻土路基温度场进行了数值模拟。考虑全球气候变暖的影响，青藏高原未来50年平均气温上升2.6℃的气候背景条件下，分析了多年冻土区路基温度场在工程运营使用中温度场的响应趋势，以及路堤结构的修筑和气候变暖分别对地层温度场带来的影响，指出了路堤下部融化盘产生的原因。 　　关键词：冻土；路基；温度场；数值模拟 　　1 引言 　　冻土是气候变化的产物。全球气候变暖的背景将会对青藏高原冻土的发展演化带来长期而又显著的影响。各类工程设施，尤其是穿越冻土区的线性工程，如公路、铁路及管线工程等，将会直接受此影响，从而带来诸多工程问题。另一方面，冻土是热敏感性土体，各类工程的修筑不可避免地改变土体天然水热状况，即工程活动引起局部冻土发生变化，诱发工程问题。因此，气候变化和工程活动是冻土区工程问题最直接的诱因，而其核心则是土体的热状况。 　　2 控制方程 　　在自然气候的影响下，路基土体发生冻融交替变化，因此冻土路基温度场问题是带相变的导热问题，控制方程是非稳态导热方程，问题域内存在着随时间移动的两相界面，界面上可吸收或释放相变潜热。在数学上这是强非线性问题，计算中存在困难[1]。问题的控制方程为二维带相变的非稳态导热微分方程： 　　在非相变区内：，则上式就转变为通常的导热微分方程式。其中为固相率，它是无因次量，固相率的增加（或减小）与相变潜热的释放（或吸收）量成正比。的含义为： 　　式中：、分别为融化及冻结温度；为相变区内节点温度（℃）；为该节点的固相率；为土冻结或融化相变潜热（）；导热系数（）和比热（）；为材料的密度（）；为材料的内热源强度（）。 　　3 计算模型 　　计算中取路基宽度为10.0，路堤高度为5.0，路基坡度为1：1.5，地基计算宽度为87.4，有坡脚各向外延伸30，路基计算深度为天然地表以下30。计算区域有四层介质组成：Ⅰ区为碎石道碴层，Ⅱ区为路堤填土层，Ⅲ区为粉质粘土层，Ⅳ区为强风化泥岩层，它们的物理参数见表-1。 　　4.1 上边界条件 　　参照有关资料[2，3]，考虑全球气候变暖的影响，取青藏高原未来50年平均气温上升2.6℃，设初始年平均气温为-4.0℃，由于路堤模型上表面温度变化不仅与环境空气有关，而且受到太阳辐射等复杂因素的综合作用，故根据附面层理论及相关资料[4，5]，对各计算模型的热边界条件设定如下： 　　6 结果分析 　　图（6）和图（7）分别是路堤修筑完成后第2年和第50年四个时间节点（10月15日、1月15日、4月15日和7月15日）的温度分布情况。从两幅图可以清楚地看出：热能和冷能在温度场云图中可以形象地看作封闭的场线构成的圆圈，热能对应中心凸起的系列圆圈，冷能对应中心下凹的系列圆圈。随着季节变更，这些圆圈由路堤顶部形成，逐渐向地层底部扩散。图（6）和图（7）对照来看，第50年的正等温线明显上移，表明地层升温，与图（5）反映的温度趋势相类似。 　　图（6） 路堤修筑完成后第2年四个时间节点的温度分布： 　　（a） 10月15日； （b） 1月15日； （c） 4月15日； （d） 7月15日 　　从局部来看，路堤修筑完成后第2年7月15日（如图（6d）所示）的冻土上限（0℃等温线）与初始位置（大约天然地表下5m处）相比几无变化，只是路堤下方的上限略有抬升，但是-3℃等温线明显下移，甚至在路堤下部消失。第50年的7月15日（如图（7d）所示）的温度分布图可以看出：路堤以下部最低温度已经升至-0.2℃，冻土上限与初始位置相比出现明显下移，最低位置在地表以下约14m处，且在路堤下存在一个巨大的融化盘，最大深度约24m。 　　从图（6）和图（7）个时间节点的温度分布对比还可以发现：气候变暖使冻土天然上限下移，但路堤的修筑使其下部上限下移得更为严重，反映这种路堤结构加速了冻土退化。 　　从图（7b）可以看出路堤内巨大的融化核上边界为y=-2m，下边界达到y=-23m。因此在气候变暖和路堤修筑的双重影响下，路地下21m范围的多年冻土发生融化，此时路堤及其下部已有很大的跨年度融化核存在，这将导致路堤发生大面积融沉，对路基稳定非常不利。 　　图（7） 路堤修筑完成后第50年四个时间节点的温度分布： 　　（a） 10月15日； （b） 1月15日； （c） 4月15日； （d） 7月15日 　　7 结论 　　（1）应用有限元软件ANSYS热分析模块对冻土路基温度场进行了数值模拟。得出了青藏高原未来50年平均气温上升2.6℃的气候条件下路基温度场的变化规律，直观看出能量通过路堤输入进地层，并直观模拟了冻土上限的下移及融化盘的形成等常见现象。 　　（2）通过比较分析，得出了气候变暖和路堤修筑对