第3期 胡向东，等：双持管冻结稳态温度场广义解析解 的取值为0.51.5.冻结管半径r。=0.054m:冻结 管表面温度T=一30℃：冻土边界温度T。=0℃ 数值模拟参数见表1及图3数值模拟模型及计算结 果见图4.数值解与解析解对比结果见图5，分 别为与第一排冻结管和第 二排冻结管相应的冻土韩 幕厚度：山1为第一排冻结管的间距，1。为第二排陈结 管的间距：L为第一排冻结管与第二排冻结管的排 间距，心为第二排冻结管的偏移量，参见图2. 5-1.0-05 表1数值模拟参数 ·面1数值解一面1解析解×面2数值解一面2解析 Tab.1 Simulatio a第1姐 组号 30 0.80 1.0 1.000.50 0.251.00 1.60 2 -10 00-000000 915-10-05005101520 面2数值解一面2解析 b第2组 第 面1 0…00… c第3组 d第4组 图3形象化的布管形式 Fig.3 Visualization of freezing pipe arrangement 915-10-0500510520 ·面1数 / 解析解面2数值 面2解析解 一面3解析 66 a网格划分 0 图单位. 图4第3组网格划分及温度云 Fig.4 mperature ficld /n 面1数值解 一面1解析 面2数值解一面2解析 第4 从上述4组数值模拟和解析解的结果对比知 道，在工程常用的参数情况下，解析公式具有较高精 图5解析解计算结果与数值计算结果对比 Fig.S on of analytical and 度，仅在冻结管处出现1℃左右的误差。 1994-2015 China Academic Joumal Electronic Publishing House.All rights http://www.cnki.ne 第３期 胡向东，等：双排管冻结稳态温度场广义解析解 的取值为０．５～１．５．冻结管半径ｒ０＝０．０５４ｍ；冻结 管表面温度Ｔｆ＝－３０ ℃；冻土边界温度 Ｔ０＝０ ℃． 数值模拟参数见表１及图３．数值模拟模型及计算结 果见图４．数值解与解析解对比结果见图５．ξ１，ξ２ 分 别为与第一排冻结管和第二排冻结管相应的冻土帷 幕厚度；ｌ１ 为第一排冻结管的间距，ｌ２ 为第二排冻结 管的间距；Ｌ 为 第 一排 冻 结 管 与 第 二 排 冻 结 管 的 排 间距，ｗ 为第二排冻结管的偏移量，参见图 ２． 表１ 数值模拟参数 Ｔａｂ．１ Ｓｉｍｕｌａｔｉｏｎｐａｒａｍｅｔｅｒｓ ｍ 组号 ξ１ ξ２ ｌ１ ｌ２ ｗ Ｌ  １ １．００ １．００ １．００ １．００ ０ １．００ ２ ０．８０ １．００ １．００ ０．５０ ０ １．００ ３ ０．８０ １．００ １．００ ０．５０ ０．２５ １．００ ４ １．６０ １．２０ ０．８０ １．６０ ０ １．６０ ａ第１组 ｂ第２组 ｃ第３组 ｄ第４组 图３ 形象化的布管形式 Ｆｉｇ．３ Ｖｉｓｕａｌｉｚａｔｉｏｎｏｆｆｒｅｅｚｉｎｇｐｉｐｅａｒｒａｎｇｅｍｅｎｔ ａ网格划分 ｂ温度云图（单位：℃） 图４ 第３组网格划分及温度云图 Ｆｉｇ．４ Ｍｅｓｈｉｎｇｆｉｇｕｒｅａｎｄｔｅｍｐｅｒａｔｕｒｅｆｉｅｌｄ ｃｏｌｏｒｍａｐｏｆＧｒｏｕｐ３ 从上述 ４ 组数值模拟和解析解的结果对比知 道，在工程常用的参数情况下，解析公式具有较高精 度，仅在冻结管处出现１ ℃左右的误差． ａ第１组 ｂ第２组 ｃ第３组 ｄ第４组 图５ 解析解计算结果与数值计算结果对比 Ｆｉｇ．５ Ｃｏｍｐａｒｉｓｏｎｏｆａｎａｌｙｔｉｃａｌａｎｄ ｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｓ ３９８