浅层黄土温度场数值分析.pdfVIP

维普资讯 第 39卷 第 4期 西 安 建 筑 科 技 大 学 学 报(自然科学版) Vo1．39 NO．4 2007年 8月 J．XianUniv．ofArch． Tech．(NaturalScienceEdition) Aug． 2007 浅层黄土温度场数值分析 王铁行 ，刘 自成，岳彩坤 (西安建筑科技大学土木工程学院，陕西 西安 710055) 摘 要 ：基于模拟辐射 、蒸发、风速等 自然条件 的浅层黄土非稳态温度场计算结果和野外裸土场地 的实测结 果，分析 了土体温度场随气候的演变过程 ，探讨 了不同月份土体温度场随深度 的变化规律 ，揭示 出土体所处深 度越小，土体温度随季节 的变化幅度越大，土体年最高温度越大，最低温度越小．进一步探讨 了土体月平均温 度最大年变化随深度增加而减小的规律，并得到了不同深度土体月平均温度最大年变化幅度的计算公式． 关键词 ：黄土 ；温度场 ；气候 ；数值分析 中图分类号 ：U416．1 文献标识码 ：A 文章编号 ：1006—7930(2007)04-0463-05 黄土高原是我国西部大开发和环境保护建设的战略要地．黄土高原公路 、铁路 、市政 、水利等领域的 工程活动 日益活跃，环境工程建设事业亦蓬勃发展，这些工程项 目基本上修筑于浅层黄土之上．由于浅 层土体大面积暴露于空间，受到诸如辐射、蒸发、降水、边坡坡率等外界因素的影响，浅层黄土温度是变 化的，土体温度变化引起水分迁移使含水量重分布，黄土因含水量增大产生湿陷，含水量增大使黄土强 度降低口]，冻融使土体 自身体积发生变化，常常导致一系列病害的发生．因此，为了保证浅层黄土工程 的安全，也为黄土高原山川I绿化事业和农业生产提供技术支持，有必要对浅层黄土温度场随气候的变化 过程进行研究．此方面研究文献尚少 ，基于此 ，以西安地 区为例 ，模拟西安气候的年变化过程，研究该地 区浅层黄土温度场随气候的变化规律，这对于确定浅层黄土温度场以及进行地下贮藏、保温设施的设计 都是有意义的． 1 黄土温度场计算模型及参数 辐射 、蒸发 、湿度 、风速等气候因素随时间是变化的，受此影响，浅层黄土温度场属非稳态相变温度 场 ，其基本方程如下 ： ([K]+ ){T}一 {P)+ ){T (1) 式中：[K]为温度刚度矩阵；[N]为非稳态变温矩阵；{T}为温度值的列向量 ；At为时间步长；{P}为合 成列阵，下标 t表示时间． 采用等参四边形单元 ，刚度矩阵和变温矩阵参数确定如下 ]： 一 肌 (D。等+D1等 s (f，” (2) 一 IIlDC J，JHHd 4-Aas (z，一i，，是，) (3) 式中 J、是、 为单元节点；为导热系数；D』c 为容积比热容；lJl、D 、Dz、H都是积分基点坐标 (， )的函数 ，并与节点坐标有关；口为对流换热系数；为换热边界长；A对非第三类边界取 0，对第三类边 界，Z等于rt时A取 1／3，Z不等于 ”时A 取 1／6． 列阵 {P}由四项组成，按下式进行计算-8]： 收稿日期 ：2006-1卜28 修改稿 日期 ：2007-07-05 基金项 目：国家 自然科学基金资助项 目；陕西省 自然科学基金资助项 目(E208) 作者简介 ：王铁行(1968一)，男，陕西富平人，博士，教授 ，主要从事黄土、冻土工程理论和实践方法等方面研究 维普资讯 464 西 安 建 筑 科 技 大 学 学 报 (自然科学版) 第 39卷 {P)= {P )+ {P )+ {P。)+ {P ) (4) 式中：{