有限元计算在多年冻土区混凝土灌注桩温度场分布中的应用.docxVIP

有限元计算在多年冻土区混凝土灌注桩温度场分布中的应用 黑龙江省混凝土及外加剂专家委员会 摘 要：针对在多年冻土地区建设青藏铁路时混凝土桥灌注桩水化放热引起周围冻土温度场变化这一实际工程问题，采用伽辽金法推导出带相变的瞬态温度场问题的有限元公式，在考虑混凝土作为放热边界的条件下综合考虑了气温变化、风速等多种因素，建立了多年冻土区混凝土桥灌注桩水化放热的传热模型，计算了由于混凝土水化放热引起的冻土温度场变化。结果表明，混凝土水化热在浇注后半年内对多年冻土的温度场影响很大，回冻时间（融化的冻土温度重新回到天然状态的时间）长达2年以上。而用粉煤灰和硅灰取代一定质量的水泥可以减少混凝土水化热对冻土热状况的影响。 关键词：多年冻土；水化热；温度场；有限元 我国北方各省和西部的青藏高原属于寒冷地区或冻土区。冻土是温度低于0℃且含有冰的土岩。而多年冻土在两年或两年以上都处于冻结状态，只有表层几米的土层处于夏融冬冻的状态[1]。在多年冻土地区的主要工程地质问题有融沉、冻胀等不良地质现象。西部开发以来，青藏铁路格尔木—拉萨段将穿越多年冻土地段约553 收稿日期：2004-07-06。 基金项目：国家自然科学基金资助项目； 哈尔滨工业大学校交叉学科基金项目（HIT. MD 2001.10）。 线路通过许多山脉和盆地，复杂的地形要求线路只能“以桥代路”。而混凝土桥灌注桩的使用却给环境带来了反作用，由于高强度、高耐久性等大体积混凝土在实际工程中的应用，混凝土在水化过程中放出大量的热，从而使周围冻土层不同程度的融化，造成建筑物的不均匀沉降。不仅如此，混凝土还是一种导热性很差的建筑材料。大体积Received date：2004-07-06。 HYPERLINK mailto:bahengjing2001@ bahengjing2001@1 传热模型的建立 1.1 控制微分方程 多年冻土的融化过程是伴随着相变的导热过程，这类问题的特点是控制方程是非稳态导热方程，区域内存在着一个随时间移动的两相界面，在该界面上放出或吸收潜热，因此称之为“移动边界问题”[2]。相界面把区域分为两部分，在固相区域DS的瞬态温度场热量平衡微分方程为： (1) 在液相区域DL的瞬态温度场热量平衡微分方程为： (2) 两个温度场在相变界面S(t)上耦合，即要满足温度连续条件和能量守恒条件： (3) (4) 式中，为物体的瞬态温度，℃；为过程进行的时间，s；和分别为冻土和融土的导热系数，W/(m·℃)；为密度，kg/m3；和分别为冻土和融土的定压比热，J/(kg·℃)；为内热源强度，W/m3。为冻结面的温度，℃；n为冻结面的法线方向；L为物质的相变潜热，W/m2。因此，伴随有相变过程的非稳态导热微分方程为： (5) 式中，fs为固相率，固相率的增加（或减小）与相变潜热的释放(或吸收)量成正比例。当fs＝0时处于液相；fs＝1时处于固相；0fs1则在两相共存相变区域内。假定相变发生在的温度范围内，那么在相变区内的导热系数为： (6) 比热为： (7) 本文利用伽辽金法解微分方程，得到的有限元计算公式为： (8) 式中，系数矩阵称为温度刚度矩阵；为非稳态变温矩阵；为未知温度值的列向量；称为等式右端项组成的列向量；下标t表示这些向量都取同一时间t时刻的值。 对于瞬态温度场问题利用向后差分得： (9) 这样在已知时刻的温度后，就可以求得时刻的温度。因此在给定初始时刻的温度值后就可以进行步进积分的计算了。 1.2 传热模型和边界、初始条件 根据多年冻土地区的环境条件和实际的工程条件，假定冻土是无限大的均匀物体，而按圆柱体浇注的混凝土为均质体，且热力学参数沿各个方向是一致的。因此可以抽象建立以下的传热模型。 图1 传热模型 Fig.1 model of heat transmission 图中的空白区域为冻土，阴影区域为混凝土。上边界AB为冻土周围的空气介质，属于第三类边界条件；右边界BCD为混凝土放热边界，可视为第二类边界条件；由于是轴对称问题，右边界DE为中心轴线，为绝热边界条件；下边界EF处于永久冻土层下限以下，常年温度保持不变，属于第一类边界条件；左边界AF可视为无限远处冻土的边缘，属于绝热边界条件；因此区域ABCDEF就是所研究的伴随有相变的非稳态温度场。 通常气温变化遵循三角函数规律： (10) 式中，为该地的大气温度，℃；A为气温变化的幅度；b是与当地环境相关的系数；t是时间，h。 经过胶结料水化放热试验得出其水化放热符合复合指数形式： (11) 式中，为龄期t时的累积水化热，KJ/㎏；为t趋近无穷大时的最终水化热，KJ/