寒区圆形截面隧道温度场的解析解.pdfVIP

第 23 卷 　第 2 期 冰 　川 　冻 　土 Vol. 23 　No. 2 2 0 0 1 年 6 月 　　　　　　　　 　　JOURNAL OF GLACIOLOGY AND GEOCRYOLOGY　　　　　　 　　　 Jun. 2 0 0 1 ( ) 文章编号 2001 寒区圆形截面隧道温度场的解析解 1 1 1 1 2 赖远明 , 　喻文兵 , 　吴紫汪 , 　何 　平 , 　张孟喜 ( 1. 中国科学院 寒区旱区环境与工程研究所冻土工程国家重点实验室 , 甘肃 兰州 　730000 ; 2. 兰州铁道学院 土木工程系 , 甘肃 兰州 　730000) 摘 　要 : 根据冻土地区的实际情况对圆形隧道正冻区和未冻区的热传导方程进行简化 , 应用无量纲量 和摄动技术对简化方程进行求解 ,给出了圆形隧道冻结过程的近似解析解. 通过和有限元数值解结果比 较 , 发现该近似解析解具有较高的精度 , 能满足工程实际的精度要求. 该解既可用于工程估算也可用来 检验数值计算的计算程序及结果. 关键词 : 寒区; 隧道 ; 温度场 ; 近似解析解 中图分类号 : TU445 文献标识码 : A 微分方程 , 应用伽辽金法导出了这一问题的有限元 1 　引言 计算公式 , 以大坂山隧道为例考察了渗流和保温层 　　对于冻土研究中遇到的大多数有相变的热传导 对隧道围岩冻深的影响规律. 到 目前为止还没有圆 问题只能求近似解和数值解 , 对于高维问题 , 即使 形隧道冻结问题的解析解. 本文应用摄动技术导出 用近似方法也很难给出解析解 , 个别问题尽管可以 圆形隧道冻结问题的近似解析解 , 该解既可用于工 用近似方法使问题简化 , 但最终结果的获得还需要 程估算也可用来检验数值计算的计算机程序. 应用数值方法. 所以, 数值方法便成为求解有相变 2 　数学模型及其求解过程 热传导问题的非常重要的方法. 尽管如此 , 近似方 法能够给出系统中各物理量和时间及空间变量之间 　　对于图 1 所示圆形隧道 , 假设大气的温度为 明确的数值关系 , 近似解析解也能用来检查数值解 Ta , 隧道内径为 a , 冻结锋面半径为 R ( t ) , 正 和计算机程序的可靠性. 此外 , 在工程实际中 , 技 冻和未冻围岩的热物理学参数保持不变 , 热量传布 术人员偏爱于简便易行又满足工程精度的简单解析 是通过热传导进行的 , 在正冻和未冻状态中的热传 解 , 以便省去繁琐的数值计算前的数据准备和输入 导方程为 : 在冻结区 Ω 内有 : 工