《高等传热学》课程教学资源（文献资料）环形单圈管冻结稳态温度场一般解析解.pdf_大学文库

第46卷第6期中南大学学报（自然科学版） .46No.6 2015年6月 Journal of Central South University (Science and Technology) June 2015 D0:10.11817i.issn.1672-7207.2015.06. 环形单圈管冻结稳态温度场一般解析解胡向东3，韩延广12 (1同济大学岩土及地下工程教育部重点实验室，上海，200092： 2.同济大学地下建筑与工程系，上海，200092) 摘要：隧道联络通道和矿山竖井冻结法施工经常采用环形布置冻结管。对于环形布管冻结温度场的解，目前仅有冻结圈内完全冻实的单圈管冻结温度场解析解，但现实冻结工程中，内部通常是未冻实的。首先应用保角变换将环形单圈管冻结模型变换为易于求解的特殊的直线型挂布的模型，再结合调和方程的边界条件可分离的性质，将日分去雀挂鬼实指杨的最后藏单要管快表森实物鞋华气析华。一九子玫值模8 对解析解的验证结果表明：在冻土幕充分交圈后，果基木致，得到的解析解具有高的精度。同时本文解析解也适用于冻结圈内完全冻实的温度场。关键词：人工地层冻结法：环形单圆冻结管：温度场：解析解：保角变换：调和方程中图分类号：TK124 文献标志码：A 文章编号：1672-7207201506-2342-08 General analytical solution to steady-state temperature field of single-circle-pipe freezing HU Xiangdong',HAN Yanguang (1.Key Laboratory of Geotechnical and Underground of the Ministry of Education Tongji University.Shanghai 200092.China: 2.Department of Geotechnical Engineering Tongji University.Shanghai00092.China) :Single-circ-ipenthd sha present,there is only one analytical solution which is aimed at the situation that the ground inside the freezing-pipe circl s frozen solid.However,a more usual situation is that the area enclosed inside the freezing-pipe circle is not frozen completely.For this situation,an analytical solution to the steady-state temperature field of single-circle-pipe freezing was proposed.The solution was derived by converting the sing-r-pipefreezing model toa special single-row-pipe freezing model using conformal mapping.and then decomposing the single-row-pipe freezine model into the single-row-pipe freezing model and a linear temperature field model according to the separability of boundary conditions for harmonic qtions.Comprison of the analytical solution with the numerical thermal analysis shows that the Key words:artificial gr mapping:harmonic equation 自然科学基资助项日(5117836,51478340：浙江省自然科学基金资助项目(LZ13E08002XPot(51178336 nce foun .du.c 1994-2016 China Academic Joumal Electronic Publishing House.All rights reserved.http://www.cnki.ne

第 46 卷第 6 期中南大学学报(自然科学版) Vol.46 No.6 2015 年 6 月 Journal of Central South University (Science and Technology) June 2015 DOI: 10.11817/j.issn.1672-7207.2015.06. 环形单圈管冻结稳态温度场一般解析解胡向东 1, 2，韩延广 1, 2 (1. 同济大学岩土及地下工程教育部重点实验室，上海，200092； 2. 同济大学地下建筑与工程系，上海，200092) 摘要：隧道联络通道和矿山竖井冻结法施工经常采用环形布置冻结管。对于环形布管冻结温度场的解，目前仅有冻结圈内完全冻实的单圈管冻结温度场解析解，但现实冻结工程中，内部通常是未冻实的。首先应用保角变换将环形单圈管冻结模型变换为易于求解的特殊的直线型排布的模型，再结合调和方程的边界条件可分离的性质，将问题分解为单排问题和线性温度场问题，最后完成单圈管冻结圈内未冻实的温度场解析解。热力学数值模拟计算对解析解的验证结果表明：在冻土帷幕充分交圈后，数值模拟结果和解析解结果基本一致，得到的解析解具有较高的精度。同时本文解析解也适用于冻结圈内完全冻实的温度场。关键词：人工地层冻结法；环形单圈冻结管；温度场；解析解；保角变换；调和方程中图分类号：TK124 文献标志码：A 文章编号：1672−7207(2015)06−2342−08 General analytical solution to steady-state temperature field of single-circle-pipe freezing HU Xiangdong1, 2, HAN Yanguang1, 2 (1. Key Laboratory of Geotechnical and Underground Engineering of the Ministry of Education, Tongji University, Shanghai 200092, China; 2. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China) Abstract: Single-circle-pipe freezing method is often applied in cross passage construction and mine shaft sinking. At present, there is only one analytical solution which is aimed at the situation that the ground inside the freezing-pipe circle is frozen solid. However, a more usual situation is that the area enclosed inside the freezing-pipe circle is not frozen completely. For this situation, an analytical solution to the steady-state temperature field of single-circle-pipe freezing was proposed. The solution was derived by converting the single-circle-pipe freezing model to a special single-row-pipe freezing model using conformal mapping, and then decomposing the single-row-pipe freezing model into the single-row-pipe freezing model and a linear temperature field model according to the separability of boundary conditions for harmonic equations. Comparison of the analytical solution with the numerical thermal analysis shows that the analytical solution is precise enough. Furthermore, the solution can also be employed for the temperature field of the frozen-solid situation. Key words: artificial ground freezing; single-circle-pipe freezing; temperature field; analytical solution; conformal mapping; harmonic equations 收稿日期：2014−04−13；修回日期：2014−08−20 基金项目(Foundation item)：国家自然科学基金资助项目(51178336，51478340)；浙江省自然科学基金资助项目(LZ13E080002)(Project (51178336, 51478340) supported by the National Natural Science Foundation of China; Project (LZ13E080002) supported by the Natural Science Foundation of Zhejiang Province, China) 通信作者：胡向东，副教授，博士生导师，从事隧道及地下工程研究；E-mail：anton.geotech@tongji.edu.cn

第6物胡向东。等：环形单网管冻结拉态温度场一般解析解 2343 人工地层冻结技术已发展成为一种成熟工法，在和双拜管冻结巴霍尔金公式的准确性进行了分析叫，井矿工程、隧道工程、地基临时加固、地下水污染招结果表明，在冻结的中后期(>0.7，其中：为冻制、废弃物掩埋等领域已有广泛应用。在施工过程中唯幕厚度，1为冻结管间距)，温度场任何一点的计算堂据连土推幕的厚度和力学性质等参数是非常重要误差不超过1℃。因此在冻结的中后期，人工冻结温的，而这些参数均依赖于族结温度场的分布。因此，温度场理论是冻结法理论的基础。目前冻结温度场的度场将可以近似当作稳态温度场。则9=0：8。计算方法主要有解析法模拟法以及数值分析方法。 0,方程1)简化为：而解析解由于其理论性强，始终是研究温度场的重要部分。到目前为止，各国专家、学者们已经取得了 (2) 系列冻结管规则布置形式下的经典稳态温度场解析解。例如单管冻结温度场公式山，两管~五管直线布其极坐标下的表达式为：冻结温度场公式，单排管冻结温度场公式以及双管冻结温度场公式等，并且沿用至今。作者对这些 80,180180 ar2 r ar 0 公式进行了完善与应用性研究5，并完成了三排管冻结温度场解析解的推导0-15，对于环形布置冻结管的情形。完成了该结国内冻实的单图管冻结温度场解折单管的数学模型及解析解解，对于内部未冻实的情祝尚无解析解答文州者基于保角变换和调和方程边界条件可分离性，单管在无限大地层形成的温度场的问题如图1所环形单圈管冻结圈内未冻实的稳态温度场解析解进行示。了求解。 1二维热传导方程冻士边界在土体中，热传递主要有3种形式：传导、对流和辐射。因为对流和辐射对温度场的影响相对于传导可以忽路不计，所以仅仅老虑传导形成的温度场陶东结管根据Fourier热传导定律和热量守恒律，可以得到二年热传导方程的形式如下： 1) 图1平西内管的温度场问其中：日为温度场的分布：k为介质的热传导系数： Fig.1 Model of single-pipe freezing 为密度：c为此热容：，x,)为单位时间单位体利系统吸收或者放出的热量。其数学模型如下：人工地层冻结过程是一个瞬态热传导过程，是个具有相变并且相变界面移动的问题。对于采用稳态温度场沂似避本温度场的活用性问顺。学术果和丁程界 6,)=,；冻结管处条件普遍接受的观点是：由于人工地层冻结是个发展相 ξ.)=0：冻土边界条件缓慢的过程，尤其是在冻结的后期，其温度场非常拉通过分离变量法可，求得满足上述调和方程的通近稳态导热温度场，在人工地层冻结法的温度场理论解为：与工程中，对此状态可按稳态导热近似求解人工冻结 0(r.B)=co++ 温度场。世界上最流行和实用的人工冻结温彦场解析解有前苏联，乳、美国和日本公式，他们均为稳 a.cos m+-.sin m4产+"）(5) 温度场公式。对此问题，作者曾针对文献3)]的单排将边界条件进行Fourier级数展开，仍为常数，即 1002016 Academic Jou al Eleetronie Publishing House.All right http://www.cnki.ne

第 6 期胡向东，等：环形单圈管冻结稳态温度场一般解析解 2343 人工地层冻结技术已发展成为一种成熟工法，在井矿工程、隧道工程、地基临时加固、地下水污染控制、废弃物掩埋等领域已有广泛应用。在施工过程中，掌握冻土帷幕的厚度和力学性质等参数是非常重要的，而这些参数均依赖于冻结温度场的分布。因此，温度场理论是冻结法理论的基础。目前冻结温度场的计算方法主要有解析法、模拟法以及数值分析方法。而解析解由于其理论性强，始终是研究温度场的重要部分。到目前为止，各国专家、学者们已经取得了一系列冻结管规则布置形式下的经典稳态温度场解析解。例如单管冻结温度场公式[1]，两管~五管直线布置冻结温度场公式[2]，单排管冻结温度场公式以及双排管冻结温度场公式[3−4]等，并且沿用至今。作者对这些公式进行了完善与应用性研究[5−9]，并完成了三排管冻结温度场解析解的推导[10−15]；对于环形布置冻结管的情形，完成了冻结圈内冻实的单圈管冻结温度场解析解[16]，对于内部未冻实的情况尚无解析解答。本文作者基于保角变换和调和方程边界条件可分离性[17]，对环形单圈管冻结圈内未冻实的稳态温度场解析解进行了求解。 1 二维热传导方程在土体中，热传递主要有 3 种形式：传导、对流和辐射。因为对流和辐射对温度场的影响相对于传导可以忽略不计，所以仅仅考虑传导形成的温度场[18]。根据 Fourier 热传导定律和热量守恒律，可以得到二维热传导方程的形式[18]如下：   2 2 2 2 k g txy , , tc c x y                  (1) 其中：θ 为温度场的分布；k 为介质的热传导系数；ρ 为密度；c 为比热容；g(t，x，y)为单位时间单位体积系统吸收或者放出的热量。人工地层冻结过程是一个瞬态热传导过程，是个具有相变并且相变界面移动的问题。对于采用稳态温度场近似瞬态温度场的适用性问题，学术界和工程界普遍接受的观点是：由于人工地层冻结是个发展相对缓慢的过程，尤其是在冻结的后期，其温度场非常接近稳态导热温度场，在人工地层冻结法的温度场理论与工程中，对此状态可按稳态导热近似求解人工冻结温度场。世界上最流行和实用的人工冻结温度场解析解有前苏联[1, 3]、美国[4]和日本公式[2]，他们均为稳态温度场公式。对此问题，作者曾针对文献[3]的单排管和双排管冻结巴霍尔金公式的准确性进行了分析[20]，结果表明，在冻结的中后期(ξ/l＞0. 7，其中 ξ 为冻土帷幕厚度，l 为冻结管间距)，温度场任何一点的计算误差不超过 1 ℃。因此在冻结的中后期，人工冻结温度场将可以近似当作稳态温度场。则 0 t    ；g(t，x， y)=0，方程(1)简化为： 2 2 2 2 0 x y         (2) 其极坐标下的表达式为： 2 2 2 22 1 1 0 r r r r            (3) 2 单管的数学模型及解析解单管在无限大地层形成的温度场的问题如图 1 所示。图 1 平面内单管的温度场问题 Fig. 1 Model of single-pipe freezing 其数学模型如下： 2 2 2 22 0 0 1 1 0 ( , ) ; ( , ) ; f r r r r r                        冻结管处条件冻土边界条件 (4) 通过分离变量法[17]，求得满足上述调和方程的通解为： 10 20  (, ) r   cc r   ln 1 2 1 ( cos sin )( ) m m m m mm m a m b m cr c r         (5) 将边界条件进行 Fourier 级数展开，仍为常数，即

2344 中南大学学报（自然料学板算45卷与0无关，则a0 b-0 示为Z-R©（共中R为极径，B为极角）.而将平面表再将冻结管和冻土边界条件代入式(5)可以得到：示为=+其中为平面中的横坐标，为其纵坐标)。将2个平面分别带入变换函数可以得到： 0=8, + (6) utiv-B+ilnR R (9) 则物平面问题经过此种变换，其像平面问题如图上式即为特鲁巴克单管冻结温度场 3所示。单圈管冻结圈内未冻实的温度场外边界 R. 31保角变换的应用在冻土帷幕交圈后，冻土边界呈波浪形状。所以冻结圈内未冻实的环形单园管冻结问题如图2所示。内边界图3平面温定场问题外边界 Fig.Model of temperature field nplane 限据保角变换，=：=h=h冬内边界。 3.2 单排管冻结的巴霍尔金问题主而 .01 单排管冻结时，在冻土帷幕交阁后，冻土帷幕边界为波浪形状，则其温度场问题如图4所示。图2单国音冻结图内未冻实的问题求解模型冻土边界 Fig.Modelofsi e-pipe freezing 其数学模型可以表达为 ,月第)管 =C ar 8+2红 0：冻结管处条付图4平面内单排管的温度场问题 Fig.4 Model ofsingle-row-pipe freezing o(习-风：体士内边界条作其数学模型为 ©（-：漆士外边界条件 a20a20 =0 3a2 共中：n为冻结管的数量，j从0到1取整数值，均)=：冻土边界条 (10 由于直接求解上述问题较为困难，现引入保角变 U1,6)=0,:冻结管处条换，将其变换为较易解决的问题。变换函数为利用调和方程的边界条件可分离性可，将族土边 (8) 界条件和冻结管处条件均讲行分解。将原问题分解为其中Z平面为原平面（即物平面）》 (平面为像平面。根无限个单管冻结问题，则针对第管可以建立如下的据复变函数可以表示一个平面中的点，因此Z平面表数学模型： 994-2016 China Academic Joumal Electronic Publishing House.All rights reserved.http://www.enki.ne

2344 中南大学学报(自然科学版) 第 46 卷与 θ 无关，则 am=0；bm=0。再将冻结管和冻土边界条件代入式(5)可以得到： 0 0 0 0 ln ln ln ln f r r r T r r        (6) 上式即为特鲁巴克单管冻结温度场。 3 单圈管冻结圈内未冻实的温度场 3.1 保角变换的应用在冻土帷幕交圈后，冻土边界呈波浪形状。所以冻结圈内未冻实的环形单圈管冻结问题如图 2 所示。图 2 单圈管冻结圈内未冻实的问题求解模型 Fig. 2 Model of single-circle-pipe freezing 其数学模型可以表达为 2 2 2 22 2 0 1 0 3 0 1 1 0 2π , ; 2π , ; 2π , ; f r r r r R Rj n R j n R j n                                                    冻结管处条件冻土内边界条件冻土外边界条件 (7) 其中：n 为冻结管的数量，j 从 0 到 n−1 取整数值。由于直接求解上述问题较为困难，现引入保角变换，将其变换为较易解决的问题。变换函数为 2 ln Z i R   (8) 其中 Z 平面为原平面(即物平面)；ζ 平面为像平面。根据复变函数可以表示一个平面中的点，因此 Z 平面表示为 Z=Re iβ (其中 R 为极径，β 为极角)。而将 ζ 平面表示为 ζ=u+iv(其中 u 为 ζ 平面中的横坐标，v 为其纵坐标)。将 2 个平面分别带入变换函数可以得到： 2 i i ln R u v R   (9) 则物平面问题经过此种变换，其像平面问题如图 3 所示。图 3 像平面温度场问题 Fig. 3 Model of temperature field in ζ plane 根据保角变换， u   ； 2 ln R v R  ； 2 1 1 ln R R   ； 3 2 2 ln R R   ； 2π l n  ； 0 w 2 R R R  。 3.2 单排管冻结的巴霍尔金问题单排管冻结时，在冻土帷幕交圈后，冻土帷幕边界为波浪形状，则其温度场问题如图 4 所示。图 4 平面内单排管的温度场问题 Fig. 4 Model of single-row-pipe freezing 其数学模型为 2 2 2 2 0 0 0 ( , ) ; ( , ) ; f x y jl jl r                       冻土边界条件冻结管处条件 (10) 利用调和方程的边界条件可分离性[17]，将冻土边界条件和冻结管处条件均进行分解，将原问题分解为无限个单管冻结问题，则针对第 j 管可以建立如下的数学模型：

第 6 期胡向东，等：环形单圈管冻结稳态温度场一般解析解 2345 1 2 1 2 2 2 2 2 0 0 0 0 0 0 0 (, ) ; ( ,) ( 2, ) (,) ; ( ,) ( 2, ) j j j j j j j j j jf j jf j jf x y jl j jl l jl l jl r j jl l r jl l r                                                                   第管冻土边界条件第管冻结管处条件 (11) 选取适当的 1 j0  ， 2 j0  ，…及 1 jf  ， 2 jf  ，…，可以保证问题(11)为单管冻结的温度场模型。则 0 f 0 0 0 ln ln ln ln j j jj j r r r r r        (12) 其中： 2 2 ( ) j r y x jl   。根据冻土边界和冻结管条件的周期性，则 θjf；θj0 为常数。即每根管温度场问题相同。单排管冻结温度场解为无限个单管冻结温度场解的叠加：   j     (13) 根据贝塞特求和公式[21]： 2 2 2π 2π ln[ ( ) ] ln ch cos y x y x jl l l             (14) 整理式(13)得： 1 2π 2π ln ch cos 2 y x C D l l          (15) 0 0 0 0 ln ln ln ln jf j C r r r              ； 0 0 0 ln ln j jf D r r       将冻土边界条件代入式(15)，得： 0 1 2π 2π ln ch cos 2 jl C D l l           (16) 再将冻结管处条件代入式(15)，得： 0 f 1 2 2π π ln ch cos 2 r jl C D l l          (17) 因为 ch2x=2sh2 x+1 ，工程中 0r l ＜＜ 1 则 0 0 π π sh r r l l  ；冻结中后期 π 1 l   ，则 π π 1 sh e 2 l l    ， 2π 2π 1 ch 1 e 2 l l     ，化简并求解式(16)和(17)得： 0 1 2π ln 2 2 C D l            ； 0 0 2π π ln f D r l l       将其代入式(15)，则 0 0 0 π 2π π ln f A r l l l                  (18) 其中： 1 2π 2π ln 2 ch cos 2 y x A l l            。上式与巴霍尔金单排管冻结温度场的解析解[3] 相同。 3.3 像平面温度场问题的解对于图 3 表示的像平面温度场问题，其数学模型为： 2 2 2 2 1 0 2 0 w 0 ( , ) ; ( , ) ; ( , ) ; f u v jl jl jl R                          冻土内边界条件冻土外边界条件冻结管处条件 (19) 根据边界条件可以分离性[17]，上述问题可以分解为求解以下 2 个问题，然后再将解加起来即可得到原问题的解，即 θ=θ1+θ2。 2 2 1 1 2 2 1 10 12 0 1 w 0 ( , ) ; ( , ) ; ( , ) ; a f b u v jl jl jl R                              冻土内边界条件冻土外边界条件冻结管处条件 (20) 2 2 2 2 2 2 2 1 2 2 2 w 0 ( , ) 0; ( , ) ; ( , ) ; a b u v jl jl jl R                          冻土内边界条件冻土外边界条件冻结管处条件 (21) 通过选取适当的 θa 可以保证式(20)为单排巴霍尔金问题。由冻土内边界和冻结管条件，根据单排巴霍尔金解(18)，其解可以表达为： 0 1 1 2π π ln f b Rw l l           1 0 1 2π 2π π ln2 ch cos 2 v u l ll                  (22)

2346 中南大学学报(自然科学版) 第 46 卷则将冻土外边界条件代入式(22)得到： 0 2 1 w 1 1 2π π ln2 ch 1 2π π 2 ln f b a R l l l l                         由 2 2π l  ＞＞1 则 2 2π 2 2π ch 1 0.5 e l l     ，化简后得到： 0 2 1 w 1 π ( ) 2π π ln f b a R l l l             (23) 而对于式(21)，由于冻结管半径很小，在不引起大的误差的前提下，为了便于计算，将冻结管处条件 2 (, ) w b  jl R   简化为 2 ( ,0) b  jl   。选取适当的 Tb可以保证其解为线性温度场。则 2 1 1 2 ( ) a v         (24) 将冻结管处条件代入(24)得到 1 1 2 a b        (25) 联立式(23)和(25)可以解出 11 2 1 2 π b A l           (26) 1 2 π ( ) a A l       (27) 其中： f 0 w 1 2 1 2 2π π 2 ln A R l l            。将式(26)代入式(22)得： 1 10 1 2π 2π π ln2 ch cos 2 A v u l ll                    (28) 将式(27)代入式(24)得： 11 21 2 121 π A v l                  (29) 12 1 2 0 12 21 π π 2 A A v l l                      (30) 其中： 1 2π 2π ln 2 ch cos 2 v u A l l               。 3.4 物平面温度场的解析解根据保角变换，u=β； 2 ln R v R  ； 2 1 1 ln R R   ； 3 2 2 ln R R   ； 2π l n  ； 0 w 2 R R R  。全部带入式(30)，得到单圈管冻结圈内未冻实的温度场解析解： 2 2 3 2 1 2 1 3 0 3 3 2 1 1 ln ln ln ln ln 2ln M R R R n n R R R R R M R R R R R             (31) 其中， 0 2 3 0 1 2 2 3 1 ln ln ln ln f M R R n nR R R R R R       ，   2 2 1 ln 2cos 2 n n R R M n R R                       。 3.5 冻土帷幕形状特征理论上讲，冻土帷幕边界应随冻结管位置呈波浪形，现考察其特征。按命题，冻土边界各处温度相等，故外边界主面(β=0°的面)与外边界界面(β=±180°/n 的面)位置处的温度之差t=0(见图 2)。由式(31)，得 2 3 2 2 2 3 13 3 3 3 3 2 1 2 3 2 ln 1 = ln + ln 2 ln 2 n n M n n R R R n R R RR R t R R R R R R R                                      (32) 在冻结中后期， R R 3 2  ； R R 3 2   ，因工程中 n 较大，则 2 3 0 n R R        ， 2 3 0 n R R         ； 3 2 n R R       ＞＞2， 3 2 n R R       ＞＞2。所以式(32)化简为 2 3 2 2 13 3 3 3 3 1 2 ln ln + ln 0 ln n n R R n R RR R R R R R R               (33) 则可得到 2 3 1 3 2 ln ln 0 R R n R R   ，即 3 3 1 R R   。从上述判定过程可以看出，在满足冻结中后期(帷幕厚度较大)和冻结管数较多(工程常用布孔参数)的情况下，即满足化简过程的参数条件时，主面厚度和界面厚度基本一致，则冻土边界形状可以当作环形，其波浪形可以忽略。这一结论不适用于不能构成环形布孔的过少冻结管的情况

第6期胡向东，等：环形单网管冻结稳态温度场一般解析解 2347 3.6解析解准确性检验此时式(31)简化为解析解推导过程中采用了一定的简化处理，故有必要对解的准确性进行验证。因为数值模拟能够较好 8 R 的反应人工冻结温度场的分布，所以用ANSYS热力 = -(34 学数值模拟来检验解析解的结果。根据一般的工程实 - 际，选取参数。环半径R2一般在2.5-8.0m:管间距/ 般取值为05-12m :帷幕厚度与管间距比位般为已有的冻结圈内部完全冻实的温度场解析解© 0.5~1.5:而且内侧冻土惟幕与外侧栋土准幕比值5/3 一般介于1和0.55：0.45之间。冻结管半径0.054m: 冻结管表面温度=-30℃：冻土边界温度a-0℃。 (瓮-a -(35) 本次数值模拟共分4组。每组的参数列于表1。 )(-a 表1单管冻结数值模拟参教 Table 1 Simulation parameters of single-cirele-pipe freezing Group R/m Ry/m R3/m n I/m 14 2534119 151.05 2 34119 30 6.0 7.0 8.0 1 401.10 46.07.08.01800.55 上述模拟参数属工程常用参数范围内，满足3.5 节的参数要求，故采用环形边界进行数值模拟。又由于问题的周期性，采取包含一个冻结管的圆心角为 360m的扇形进行模拟。其中主面为00°的面：界面为=士180m的面。第1组模拟的网格划分图和温度场云图如图5和6所示。各组模拟对比结果如图7 所示通过上述对比结果可以知道，数值解和解析解基本一致。说明在工程中常用的参数范围内，本文求得的单圈冻结圈内未冻实的温度场的解析解具有较高的准确性。不仅说明了采用圆弧面边界模型讲行数伯模拟是可行的，同时也说明了在工程常用布孔参数范围图5第1组数值模椒模型及网格刻分内并在冻结的中后期，波浪形边界可以看作环形边界 Fig.Simulation model and meshing figure of first group 的结论的正确性。 37与冻结圈内冻实的温度场解析解对比对于式(31)，若考虑R,趋向于0，即极限时，冻结图内处于完全冻实状态，则 in-h飞 R2 In R Rs -InR -0雾男86西 =n-nR1 图6第1组数值模拟结米温度云图 Fig.6 Temperature field color map of first group 1994-2016Chim al Electronic Publishing House.All right erved. www.cnki.ne

第 6 期胡向东，等：环形单圈管冻结稳态温度场一般解析解 2347 3.6 解析解准确性检验解析解推导过程中采用了一定的简化处理，故有必要对解的准确性进行验证。因为数值模拟能够较好的反应人工冻结温度场的分布，所以用 ANSYS 热力学数值模拟来检验解析解的结果。根据一般的工程实际，选取参数。环半径 R2一般在 2.5~8.0 m；管间距 l 一般取值为 0.5~1.2 m；帷幕厚度与管间距比值一般为 0.5~1.5；而且内侧冻土帷幕与外侧冻土帷幕比值 ξ1/ξ2 一般介于 1 和 0.55:0.45 之间。冻结管半径 0.054 m；冻结管表面温度 θf=−30 ℃；冻土边界温度 θ0=0 ℃。本次数值模拟共分 4 组。每组的参数列于表 1。表 1 单圈管冻结数值模拟参数 Table 1 Simulation parameters of single-circle-pipe freezing Group R1/m R2/m R3/m ξ1/ξ2 n l/m 1 1.4 2.5 3.4 11/9 15 1.05 2 1.4 2.5 3.4 11/9 30 0.52 3 6.0 7.0 8.0 1 40 1.10 4 6.0 7.0 8.0 1 80 0.55 上述模拟参数属工程常用参数范围内，满足 3.5 节的参数要求，故采用环形边界进行数值模拟。又由于问题的周期性，采取包含一个冻结管的圆心角为 360°/n 的扇形进行模拟。其中主面为 θ=0°的面；界面为 θ=±180°/n 的面。第 1 组模拟的网格划分图和温度场云图如图 5 和 6 所示。各组模拟对比结果如图 7 所示。通过上述对比结果可以知道，数值解和解析解基本一致。说明在工程中常用的参数范围内，本文求得的单圈冻结圈内未冻实的温度场的解析解具有较高的准确性。不仅说明了采用圆弧面边界模型进行数值模拟是可行的，同时也说明了在工程常用布孔参数范围内并在冻结的中后期，波浪形边界可以看作环形边界的结论的正确性。 3.7 与冻结圈内冻实的温度场解析解对比对于式(31)，若考虑 R1 趋向于 0，即极限时，冻结圈内处于完全冻实状态，则 1 1 2 3 12 2 1 3 3 0 0 3 23 1 2 1 ln ln ln ln lim lim ln ln ln ln ln R R R R RR R R R R R RR R R R       ， 1 1 2 2 2 2 1 13 3 0 0 3 3 1 1 ln ln ln lim lim 1 ln ln ln R R R R R RR R R R R R       此时式(31)简化为： 2 3 2 3 2 1 2 0 2 ln 2ln 2cos n f n n n n R RR R R R n nR R R R                                  (34) 已有的冻结圈内部完全冻实的温度场解析解[16] 为： 2 3 2 3 2 1 2 0 2 2cos ln 2ln 2cos n f n n n n R n RR R R R n nR R R R                                    (35) 图 5 第 1 组数值模拟模型及网格划分 Fig. 5 Simulation model and meshing figure of first group 图 6 第 1 组数值模拟结果温度云图 Fig. 6 Temperature field color map of first group

第6期胡向东。等：环形单网管冻结检态温度场一般解析解 2349 wall by 报 Pysics.20101333-138 HU Xia ons.GUO Wang.ZHANG 胡向东，黄峰。白楠，考电土层冻结温度时人工冻结温度场模型中图矿业大学学报（自然科学，2008,3 inn infinite region Jou of ShanghaiJ 550-555 2012.47911367.137 HU Xianedone.HUANG Feng.BAl Nan.Models of artificia 14胡向东，汪洋。三排管冻结温度场的势函数叠加法解析解 sidering soil freezing point 岩石力学与工程学报，2012,3151071-1080. China University of Mining Technology(Natr HU Xiangdong.WANG Yang.Analytical solution of 8 cience1200837M4r550-555 three-row-piped frozen temperature field by means of [刀胡向东，白楠.余锋.单择管冻结温度场Trupak和Bakholdin superposition of potentia Chinese Journal of Rock 公式的适用性同济大学学报（自然科学版，2008,367 s and En ing2012.31(511071-1080 906-910. [Ⅱ]胡向东，任辉。3排管冻结梯形-抛物弓叠合等效温度场模型 HU Xi edong.BAl Nan.YU Feng.Analysis of Trupak and 和平均温度U几.煤炭学报.2014.39列1上78-83」 Bakholdin formulas for temperature field of single-row-pipe HU Xiangdong,REN Hui.A traperoidal-parabolic arch rozen soil walllJ]Journal of Tongji University (Natural ture field model and its averape Science,)2008,36(7作906-910， temperature formula for triple-row-pipe freezing Joual of 8 胡向东.直线形单挂管冻士幕平均温度计算方法冰川 China Coal So iety,2014.3%178-83 冻土.2010.32(4：778-78 6胡向东，陈锦，汪洋，等环形单圈管冻结稳态温度场解析解 HU Xiangdong A calculation for the straigh .岩士力学2013 soil wallJ]Joumal of Glaciology and CHEN Jin,WANG Yang.et al.Analytical 胡向东 ,赵 2010.3 直线双排管结平均温度的等 371k28 上之西 2012.37028-2 ng.QIN Te nM]s s520035 ZHANG LO 1060.1 r adiabatic/CDMA2013.Qindao.2013: 9jLM.Heat conduct n[M].Berlin:Springer Berlin 7.360 [11]HU Xiangdong.ZHANG Luoyu Analytical solution to 2训胡向东，赵俊杰人工冻结温度场巴霍尔金模型准确性研究 rature ficld of two freczing pipes with n地下空间与T积学报20101上g6-10 different temperatures Joural of Shanghai Jiaotong HU Xianedone ZHAO Juniie.Research on precision of University (Science)2013.1816 706-711 Bakholdin model for temperature field of artificial ground 2]胡向东，郭旺.张洛喻.无限大区域内少量冻结管稳态温度 freezing[J].Chinese Journal of Underground Space and 场解析解几煤炭学报.2013.3811：1933.1960 ing2010.6196-101. Hu Xiansdong.Guo Wang.Zhang Luovu.Analytical solution of 1211 Hurst W.Reservoir engi g and conformal mapping of oil steady state temperature ficld of a few freezing pipes in infin andgs ficld[M].Tuls Oklahoma:The Petroleum Publishing region].Journal of Chins Coal Society,2013, Company.,1979:115-118. 38M11上1953.-1960 (编辑陈爱华) 1994-2016 China Academic Joumal Eleetronic Publishing House.All rights reserved. http://www.cnki.net

第 6 期胡向东，等：环形单圈管冻结稳态温度场一般解析解 2349 frozen soil wall by equivalent trapezoid method[C]//AIP Conference Proceedings. New York: American Institute of Physics, 2010: 1333−1338. [6] 胡向东, 黄峰, 白楠. 考虑土层冻结温度时人工冻结温度场模型[J]. 中国矿业大学学报(自然科学版), 2008, 37(4): 550−555. HU Xiangdong, HUANG Feng, BAI Nan. Models of artificial frozen temperature field considering soil freezing point[J]. Journal of China University of Mining & Technology (Natural Science), 2008, 37(4): 550−555. [7] 胡向东, 白楠, 余锋. 单排管冻结温度场 Тrupak 和 Bakholdin 公式的适用性[J]. 同济大学学报(自然科学版), 2008, 36(7): 906−910. HU Xiangdong, BAI Nan, YU Feng. Analysis of Trupak and Bakholdin formulas for temperature field of single-row-pipe frozen soil wall[J]. Journal of Tongji University (Natural Science), 2008, 36(7): 906−910. [8] 胡向东. 直线形单排管冻土帷幕平均温度计算方法[J]. 冰川冻土, 2010, 32(4): 778−785. HU Xiangdong. Average temperature calculation for the straight single-row-pipe frozen soil wall[J]. Journal of Glaciology and Geocryology, 2010, 32(4): 778−785. [9] 胡向东, 赵飞, 佘思源, 等. 直线双排管冻结壁平均温度的等效抛物弓形模型[J]. 煤炭学报, 2012, 37(1): 28−32. HU Xiangdong, ZHAO Fei, SHE Siyuan, et al. Equivalent parabolic arch method of average temperature calculation for straight double-row-pipe frozen soil wall[J]. Journal of China Coal Society, 2012, 37(1): 28−32. [10] HU Xiangdong, ZHANG Luoyu. Analytical solution to steady-state temperature field of one and two freezing pipes near linear adiabatic boundary[C]//ICDMA 2013. Qindao, 2013: 257-260. [11] HU Xiangdong, ZHANG Luoyu. Analytical solution to steady-state temperature field of two freezing pipes with different temperatures[J]. Journal of Shanghai Jiaotong University (Science), 2013, 18(6): 706-711. [12] 胡向东, 郭旺, 张洛瑜. 无限大区域内少量冻结管稳态温度场解析解[J]. 煤炭学报, 2013, 38(11): 1953-1960. Hu Xiangdong, Guo Wang, Zhang Luoyu. Analytical solution of steady state temperature field of a few freezing pipes in infinite region[J]. Journal of China Coal Society, 2013, 38(11):1953-1960. [13] 胡向东, 郭旺, 张洛瑜. 无限大区域内四管冻结的稳态温度场解析解[J]. 上海交通大学学报, 2012, 47(9): 1367-1371. HU Xiangdong, GUO Wang, ZHANG Luoyu. Analytical solution of steady state temperature field of four freezing pipes in an infinite region[J]. Journal of Shanghai Jiaotong University, 2012, 47(9): 1367-1371. [14] 胡向东, 汪洋. 三排管冻结温度场的势函数叠加法解析解[J]. 岩石力学与工程学报, 2012, 31(5): 1071−1080. HU Xiangdong, WANG Yang. Analytical solution of three-row-piped frozen temperature field by means of superposition of potential[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(5): 1071−1080. [15] 胡向东, 任辉. 3 排管冻结梯形-抛物弓叠合等效温度场模型和平均温度[J].煤炭学报, 2014, 39(1):78-83. HU Xiangdong, REN Hui. A trapezoidal-parabolic arch superimposed equivalent temperature field model and its average temperature formula for triple-row-pipe freezing[J]. Journal of China Coal Society, 2014, 39(1): 78-83. [16] 胡向东, 陈锦, 汪洋, 等. 环形单圈管冻结稳态温度场解析解 [J]. 岩土力学,2013, 34(3):874−880. HU Xiangdong, CHEN Jin, WANG Yang, et al. Analytical solution to steady-state temperature field of single-circle-pipe freezing[J]. Chinese Journal of Rock and Soil Mechanics, 2013, 34(3): 874−880. [17] 陈恕行, 秦铁虎, 周忆. 数学物理方程[M]. 上海: 复旦大学出版社, 2003: 35. CHEN Shuxing, QIN Tiehu, ZHOU Yi. Mathematical physics equation[M]. Shanghai: Fudan University Press, 2003: 35. [18] Carslaw H S, Jaeger J C. The conduction of heat in solids[M]. London: Oxford University Press, 1959: 1. [19] Jiji L M. Heat conduction[M]. Berlin: Springer Berlin Heidelberg, 2009: 2−9. [20] 胡向东, 赵俊杰. 人工冻结温度场巴霍尔金模型准确性研究 [J]. 地下空间与工程学报, 2010, 6(1): 96−101. HU Xiangdong, ZHAO Junjie. Research on precision of Bakholdin model for temperature field of artificial ground freezing[J]. Chinese Journal of Underground Space and Engineering, 2010, 6(1): 96−101. [21] Hurst W. Reservoir engineering and conformal mapping of oil and gas field[M]. Tulsa. Oklahoma: The Petroleum Publishing Company, 1979: 115−118. (编辑陈爱华)