高温高含冰量冻土路基流变特性数值分析 Numerical Analysis of Rheological Character of Embankment with Warm and Ice-Rich Frozen Soil.pdfVIP

第29卷，第5期 中 国 铁 道 科 学 V01．29No．5 200 8年9月 CHINARAII。WAYSCIENCE September，2008 文章编号：1001—4632(2008)05—0013—07 含冰量冻土路基流变特性数值分析 许健，牛富俊，林战举 (中国科学院寒区旱区环境与工程研究所冻土工程国家重点实验室，甘肃兰州 730000) 摘 要：推导冻土相变温度场的数学模型及有限元计算公式，在此基础上分别利用热弹塑性和热弹塑性一蠕 变的本构关系，考虑融土和冻土的应力和变形，给出相应的有限元计算公式。以青藏铁路北麓河试验段为计算 模型，计算冻土路基的温度场和温度影响下的应力、变形场，并且与实测温度和路基变形数据进行对比分析。计 算和分析结果表明：冻土路基下多年冻土层温度随着路基使用年限的增加而逐步升高，形成高温冻土层；冻土路 基变形随着时间的增加而增加，最后趋于稳定，这表明冻土的蠕变为衰减型蠕变；路基修筑后3～5年的变形量 占变形的绝大部分，此后路基变形进入一个相对稳定的阶段，可满足铁路的正常运营。 关键词：高温高含冰量冻土；流变；热弹塑性一蠕变；有限元法；冻土路基 中图分类号：U213．14文献标识码：A 关于冻土的力学特性，一些专家学者进行了研 一A婴一h(T—L) (2) 究[1。4]。但是，对冻土尤其是高温高含冰量冻土流 变[5’6]力学特性的分析，尚未建立严密的分析理论式中：L为环境温度；q为热流密度边界条件；甩 和设计方法。为此，本文采用弹塑性一蠕变的分离 为边界面的外法线方向；h为对流换热系数。 初始条件为 型本构关系模拟高温高含冰量冻土的流变特性，同 时考虑与温度的耦合作用，对青藏铁路北麓河沉降 Tl，。一To (3) 构造温度的插值函数： 试验段高温高含冰量冻土路基(以下简称“冻土路 基”)进行计算分析，以便更有效地对高温冻土区 T—ENi(z，y)L(￡) (4) 的工程稳定性做出预报。 式中：Ni为单元型函数；L为单元各节点的温度 基本理论及有限元方程 值。 应用Galerkin法‘8|，可求得二维相变温度场的 1．1 热量控制微分方程及有限元方程 有限元方程，即 考虑介质的热传导、冰水相变而忽略热对流及 EK-]{T}“M]f可aT)一{F) (5) 其他作用，并认为未冻水含量仅是温度的函数时， 其中 利用显热容法E73得到冻土路基二维相变温度场控制 方程，即 CNiNjm ㈣ 坞一∑j-n ㈣ ON,差+ c等=去(A罢)+寺(A哥) K。一∑．f盯A 式中：A，C分别为等效导热系数和等效热容。 固定边界条件为 可O (8) ‘ T；L，一A婴：q， dn