基于FEPG的多圈管冻结温度场有限元计算软件的研发
|
摘要
在深井冻结法凿井工程中,多圈管冻结壁温度场的分析计算是冻结壁设计的基础。因此,研究掌握多圈管冻结条件下冻结壁温度场的计算方法,具有非常重要的理论和现实意义。
国际通用的有限元软件在温度场的数值模拟方面有着广泛的应用,但其价格高昂,对人工多圈管冻结法凿井工程中的一些实际问题应用起来并不方便,比如分期冻结方案的模拟问题、冻结管发生断裂后的温度场计算问题、考虑冻结管偏斜条件下的计算模型建立问题等。FEPG(Finite Element Program Generator)是一种有限元程序自动生成系统,其最大特点是可以根据虚功方程产生计算所需的源程序,用户可以直接修改内部参数,控制计算过程。本文即利用FEPG平台,开发多圈管冻结温度场的有限元计算程序。
论文对如何生成该源程序进行了深入细致的分析,对如何用FEPG的有限元(描述)语言编写控制方程文件及计算方法文件、各个文件中变量含义、子过程功能和文件接口做了详细说明,这对进一步利用FEPG程序做二次开发有重要意义。文中还对FEPG.GID及FEPG.GID的二次开发语言TCL/TK语言作了详细的介绍,并对FEPG.GID的二次开发接口作了详细说明,这对于今后的接口开发有重要的参考价值。
论文的研究成果为如何利用FEPG进行有限元计算程序开发提供了参考。论文最后对一个工程实例进行了模拟,模拟结果符合工程实际状况,可为实际工程应用提供参考。
In freezing method deep shaft sinking projects, analyzing and computing the temperature field of frozen wall is the base of design of the frozen wall. So, it has extremely important theoretical and realistic significance to study and master the computing method of the freeze wall temperature field under multi-freezing-tube cycles freezing condition.
The international general finite element software has the widespread application in the temperature field simulation aspect. But it is very expensive, and it is not convenient when it processing some actual problems in the shaft sinking projects using artificial multi-freezing-tube cycles freezing method, for example: how to simulate periodic freezing situation problem, how to create geometry model quickly and conveniently under deflective pipes? FEPG (Finite Element Program Generator) is a general purpose finite element program generator. For any kind of finite element modeling, FEPG can automatically generate complete source code based FEM (Finite Element Method) algorithm expressions. Users can control the process of calculation on modifying the internal parameters. The program special for artificial multi-freezing-tube cycles freezing temperature field is based on the source code which is generated by FEPG.
The generating procedure is described in detail, and this article also gives a detailed introduction to the method of programming of the control equation files and computing method files, deeply analysis to the program shows the meanings of parameters, function of subroutines and interface files. These works will be very helpful to the development of this program. The article gives a detailed introduction to FEPG.GID and the FEPG.GID's secondary development language—TCL/TK language, and have very important reference value to the interface development in future.
The findings in this paper provide reference for finite element program development using FEPG. In the last part of this paper, we simulate a project example, the simulating result is conformity with actual project situation, might provide reference for the actual project application.
国际通用的有限元软件在温度场的数值模拟方面有着广泛的应用,但其价格高昂,对人工多圈管冻结法凿井工程中的一些实际问题应用起来并不方便,比如分期冻结方案的模拟问题、冻结管发生断裂后的温度场计算问题、考虑冻结管偏斜条件下的计算模型建立问题等。FEPG(Finite Element Program Generator)是一种有限元程序自动生成系统,其最大特点是可以根据虚功方程产生计算所需的源程序,用户可以直接修改内部参数,控制计算过程。本文即利用FEPG平台,开发多圈管冻结温度场的有限元计算程序。
论文对如何生成该源程序进行了深入细致的分析,对如何用FEPG的有限元(描述)语言编写控制方程文件及计算方法文件、各个文件中变量含义、子过程功能和文件接口做了详细说明,这对进一步利用FEPG程序做二次开发有重要意义。文中还对FEPG.GID及FEPG.GID的二次开发语言TCL/TK语言作了详细的介绍,并对FEPG.GID的二次开发接口作了详细说明,这对于今后的接口开发有重要的参考价值。
论文的研究成果为如何利用FEPG进行有限元计算程序开发提供了参考。论文最后对一个工程实例进行了模拟,模拟结果符合工程实际状况,可为实际工程应用提供参考。
In freezing method deep shaft sinking projects, analyzing and computing the temperature field of frozen wall is the base of design of the frozen wall. So, it has extremely important theoretical and realistic significance to study and master the computing method of the freeze wall temperature field under multi-freezing-tube cycles freezing condition.
The international general finite element software has the widespread application in the temperature field simulation aspect. But it is very expensive, and it is not convenient when it processing some actual problems in the shaft sinking projects using artificial multi-freezing-tube cycles freezing method, for example: how to simulate periodic freezing situation problem, how to create geometry model quickly and conveniently under deflective pipes? FEPG (Finite Element Program Generator) is a general purpose finite element program generator. For any kind of finite element modeling, FEPG can automatically generate complete source code based FEM (Finite Element Method) algorithm expressions. Users can control the process of calculation on modifying the internal parameters. The program special for artificial multi-freezing-tube cycles freezing temperature field is based on the source code which is generated by FEPG.
The generating procedure is described in detail, and this article also gives a detailed introduction to the method of programming of the control equation files and computing method files, deeply analysis to the program shows the meanings of parameters, function of subroutines and interface files. These works will be very helpful to the development of this program. The article gives a detailed introduction to FEPG.GID and the FEPG.GID's secondary development language—TCL/TK language, and have very important reference value to the interface development in future.
The findings in this paper provide reference for finite element program development using FEPG. In the last part of this paper, we simulate a project example, the simulating result is conformity with actual project situation, might provide reference for the actual project application.
引文
[1]Brent B.Welch著,王道义,乔陶鹏等译.Tcl/Tk组合教程(第二版).电子工业出版社,2001.1
[2]Brent B.Welch著,崔凯译.Tcl/Tk编程权威指南.中国电力出版社,2002.6
[3]罗朝盛,余文芳,余平.Visual Basic6.0程序设计教程.北京:人民邮电出版社,2002.2
[4]徐谬.Visual Basic应用与开发案例教程.北京:清华大学出版社,2005.1
[5]曹青,邱李华,郭志强.Visual Basic程序设计教程.北京:机械工业出版社,2000.2
[6]赵经文,王宏钮编.有限元方法.哈尔滨工业大学出版社,1998:18 45
[7]王歇成.有限单元法.清华大学出版社,2005:131-137
[8]谭浩强,田淑清.Fortran语言-Fortran77结构化程序设计.清华大学出版社,2002:303-312
[9]飞箭软件公司.FEPG5.2初级教程,2002,12
[10]飞箭软件公司.FEPG5.2中级教程,2002,12
[11]飞箭软件公司.GID使用说明书,2002,12
[12]飞箭软件公司.基于FEPG的有限元方法,2002,12
[13]钱华山.在线生成有限元程序系统FEPG.计算机辅助设计与制造,2000年12期
[14]杨友卿.一种高效的有限元自动编程工具-FEPG系统.力学与实践,1998年02期
[15]徐学祖,王家澄,张立新.冻土物理学.北京:科学出版社,2001.4
[16]徐学祖,邓友生.冻土中水分迁移的实验研究.北京:科学出版社,1991
[17]徐学祖,王家澄等.温度梯度诱导薄膜水迁移的冻胀机理.科学通报,1997,42(9):576-579
[18]郭瑞平,霍雷声.冻结壁位移计算及其冻结施工优化设计.矿冶工程,1999,19(4):6-8
[19]汪仁和,王伟.冻结孔偏斜下冻结壁温度场的形成特征与分析.岩土工程学报,2003,25(6):658-661
[20]陈湘生.对深冻结井几个关键问题的探讨.煤炭科学技术,1999,27(1):36-38
[21]崔广心.厚表土层湿土结冰温度与冻结壁厚度确定的研究.中国矿业大学学报,1997,26(3):1-4
[22]汪仁和.冻结温度场的叠加计算与计算机方法.安徽理工大学学报[J],2003,23(1):25-29
[23]崔广心.深厚冲积层中冻结壁厚度研究[J].冰川冻土,1995,17(增刊):26-33
[24]Ma Wei.Review and prospect of the studies of ground freezing technology in China[J].Journal of Glaciology and Geocryology,2001,23(3):218-224
[25]付厚利.深厚表土层中冻结壁可靠性设计理论初探.建井技术,1999,20(4):32-34
[26]汪仁和,曹蓉斌.双排管冻结下冻结壁温度场形成特征的数值分析,2002,24(2):181-185
[27]崔广心.深土冻土力学-冻土力学发展的新领域.冰川冻土.1998,20(2):97-100
[28]胡省三,成玉琪.21世纪前期煤炭科技重点发展领域探讨.煤炭学报,2005,30(1):1-7
[29]崔广心.厚表土中地下工程特殊施工.陈明和,王建平等主编.矿山建设工程新进展(上册):中国矿业大学出版社,2005.8:1-3
[30]朱林楠.深井人工冻结壁温度场分析,冰川冻土,1981 4(3):1-4
[31]黄德发,邓文芳.冻结法凿井在河南煤矿立井施工中的应用.陈明和,王建平等主编.矿山建设工程新进展(上册):32-39.中国矿业大学出版社,2005.8
[32]佟凤健,朱德明.郭屯矿井建设若干技术问题研究.陈明和,王建平等主编.矿山建设工程新进展(上册):53-57.中国矿业大学出版社,2005.8
[33]卢相忠,宋世宏.丁集主副井冻结工程冻结壁分析.陈明和,王建平等主编.矿山建设工程新进展(上册):165-172.中国矿业大学出版社,2005.8
[34]刘明凯.丁集矿井深厚表土冻结技术.陈明和,王建平等主编.矿山建设工程新进展(上册):173-177.中国矿业大学出版社,2005.8
[35]李连华,王然芳.龙固矿副井深厚表土层冻结法施工.陈明和,王建平等主编.矿山建设工程新进展(上册):29-132.中国矿业大学出版社,2005.8,1
[36]李功洲,王长生.程村矿 主副井井筒深厚冲积层冻结法凿井技术研究.陈明和,王建平等主编.矿山建设工程新进展(上册):73-80.中国矿业大学出版社,2005.8
[37]林鸿苞,江新春,姚飞建.山东巨野矿区郭屯矿井特厚深表土复杂地层井简施工方法的选择[A].周兴旺,李明远.全国矿山建设学术会议论文选集(上册)[C].全国矿山建设学术会议论文选集(上册):94-99.中国矿业大学出版社,2004.8
[38]李述训.立井冻结法凿井工程中的热工计算.冰川冻土,1994,16(1)
[39]臧桂茂,谭杰.谈巨野煤田冻结法凿井技术问题.建井技术,2001(2)
[40]翁家杰.井巷特殊施工.北京:煤炭工业出版社.1991:4-5
[41]沈季良,崔云龙,王介峰.建井工程手册.北京:煤炭工业出版社.1986:60-317
[42]郭兰波等.竖井冻结壁温度场的有限元分析.中国矿业学院学报,1981,No.3
[43]洪伯潜.特殊凿井技术在我国的发展与前景.中国煤炭,2000,26(4):68-70
[44]李功洲等.深厚冲积层冻结凿井技术问题的探讨.周兴旺,李明远.全国矿山建设学术会议论文选集(上册)[C].全国矿山建设学术会议论文选集(上册):284-292:中国矿业大学出版社,2004.8
[45]苏立凡.对我国建设700m冻结立井的技术要素分析[A].王长生,苏立凡,张文.地层冻结工程技术和应用-中国地层冻结工程40年论文集[C].北京:煤炭工业出版社,1995
[46]Ma Wei(马巍).Review and Prospect ofthe Studies of Ground Freezing Technology in China(中国地层土冻结技术研究的回顾与展望).冰川冻土,2001,23(3):218-224
[47]苏立凡,徐兵壮.我国地层人工冻结工程技术的进展,应用和问题[A].周兴旺,张向东.全国矿山建设学术会议论文选集(上册)[C]:122-125.吉林:东北大学出版社,2003.8
[48]陈湘生.深冻结壁时空设计理论.岩土工程学报,1998,120(5):13-16
[49]虞相,王正廷,苏立凡.我国地层冻结技术的新发展-纪念冻结法凿井40周年[A].王长生,苏立凡,张文.地层冻结工程技术和应用-中国地层冻结工程40年论文集[C].北京:煤炭工业出版社,1995
[50]崔广心.关于冻结壁几个问题的探讨[A].王长生,苏立凡,张文.地层冻结工程技术和应用-中国地层冻结工程40年论文集[C].北京:煤炭工业出版社,1995
[51]姚志明,周世虎.立井特厚粘土层冻结施工技术.煤炭工程,2004(12):26-27
[52]李银桥等.深立井复杂冲积层冻结技术.全国矿山建设学术会议论文集,2004
[53]陈瑞杰,陈国栋,李述训等.人工地层冻结应用研究进展和展望.岩土工程学报.,2000,22(1):40-44
[54]陈文豹等.潘集矿区冻结壁平均温度及冻结孔布置圈径的探讨,煤炭学报,Nol,1982
[55]刘珣,梁鹏.城市地下工程中人工冻结法的防冻胀优化设计研究.武汉科技大学学报,2004,27(4):387-390
[56]陆卫国,郭圣澏,林稚华.广州地铁建设中冻结工法的应用分析.西部探矿工程.,2004(92):108-109
[57]郭晓江.冻结法在广州地铁二号线暗挖隧道中的应用.煤炭工程,2001(12):27-29
[58]周兴荣,王宗金等.冻结法在广州丫髻沙大桥桩基事故处理中的应用,2000.21(2):38-40
[59]程知言等.上海某冻结加固工程试验研究分析,地质与勘探,2002,38(3)
[60]程国栋,马巍.国际冻土工程研究进展-第五届冻十工程国际学术讨论会综述.冰川冻土,2003,25(3):303-308
[61]程国栋.冻土力学与工程的国际研究新进展.地球科学进展,2001,16(3):293-299
[62]朱正刚.寒区路基工程与多年冻土间相互作用问题研究进展.结构工程师,2004,20(6):43-48
[63]汪仁和,张瑞.多圈冻结管下冻土扩展速度的实测研究.陈明和,王建平等主编.矿山建设工程新进展(上册),2005.8:86-91
[64]周晓敏,张云利.竖井单圈布孔冻结壁形成平面温度场问题的数值试验研究.陈明和,王建平等主编.矿山建设工程新进展(上册),2005.8:104-110
[65]张开顺,李炳胜.井筒冻结壁径向位移规律计算研究.陈明和,王建平等主编.矿山建设工程新进展(上册),2005.8:116-120
[66]雷晓燕.岩土工程数值计算.中国铁道出版社.1999.11
[67]谢康和,周健.岩土工程有限元分析理论与应用.科学出版社.2003.5
[68]张铜生,张富德.简明有限元法及其工程应用.地震出版社.1990.8
[69]赵阳升等.冻结壁的温度场和弹塑性应力场的耦合分析.煤炭学报,1985.6
[70]朱林楠.深井人工冻结壁温度场分析.冻川冻土,1981.11
[71]张燕.预报冻结壁形成及其温度分布的有限元方法.煤炭学报,1979.2
[72]汤志斌等,冻结壁温度场的实测研究及在工程中的应用,地层冻结技术和应用--中国地层冻结工程40年论文集,1995.10
[73]赵阳升等.冻结壁的温度场和弹塑性应力场的耦合分析.煤炭学报,1985.6
[74]杨世铭编.传热学.人民教育出版社,1998.3
[75]俞昌铭编著.热传导及其数值分析.清华大学出版社,1982
[76]周晓敏,张云利.竖井单圈布孔冻结壁形成平面温度场问题的数值试验研究.矿山建设工程新进展(上册),2005.8.104-110
[77]李述训.人工冻结凿井中热学问题计算方法评价.世界冰川冻土,1982(3):23-32
[78]马思民,王金芝,李洪升.冻土温度场的参数辨识.冰川冻土.2002,24(3):299-303
[79]谢康和,周健.岩土工程有限元分析理论与应用.科学出版社.2003.5
[80]张铜生,张富德.简明有限元法及其工程应用.地震出版社.1990.8
[2]Brent B.Welch著,崔凯译.Tcl/Tk编程权威指南.中国电力出版社,2002.6
[3]罗朝盛,余文芳,余平.Visual Basic6.0程序设计教程.北京:人民邮电出版社,2002.2
[4]徐谬.Visual Basic应用与开发案例教程.北京:清华大学出版社,2005.1
[5]曹青,邱李华,郭志强.Visual Basic程序设计教程.北京:机械工业出版社,2000.2
[6]赵经文,王宏钮编.有限元方法.哈尔滨工业大学出版社,1998:18 45
[7]王歇成.有限单元法.清华大学出版社,2005:131-137
[8]谭浩强,田淑清.Fortran语言-Fortran77结构化程序设计.清华大学出版社,2002:303-312
[9]飞箭软件公司.FEPG5.2初级教程,2002,12
[10]飞箭软件公司.FEPG5.2中级教程,2002,12
[11]飞箭软件公司.GID使用说明书,2002,12
[12]飞箭软件公司.基于FEPG的有限元方法,2002,12
[13]钱华山.在线生成有限元程序系统FEPG.计算机辅助设计与制造,2000年12期
[14]杨友卿.一种高效的有限元自动编程工具-FEPG系统.力学与实践,1998年02期
[15]徐学祖,王家澄,张立新.冻土物理学.北京:科学出版社,2001.4
[16]徐学祖,邓友生.冻土中水分迁移的实验研究.北京:科学出版社,1991
[17]徐学祖,王家澄等.温度梯度诱导薄膜水迁移的冻胀机理.科学通报,1997,42(9):576-579
[18]郭瑞平,霍雷声.冻结壁位移计算及其冻结施工优化设计.矿冶工程,1999,19(4):6-8
[19]汪仁和,王伟.冻结孔偏斜下冻结壁温度场的形成特征与分析.岩土工程学报,2003,25(6):658-661
[20]陈湘生.对深冻结井几个关键问题的探讨.煤炭科学技术,1999,27(1):36-38
[21]崔广心.厚表土层湿土结冰温度与冻结壁厚度确定的研究.中国矿业大学学报,1997,26(3):1-4
[22]汪仁和.冻结温度场的叠加计算与计算机方法.安徽理工大学学报[J],2003,23(1):25-29
[23]崔广心.深厚冲积层中冻结壁厚度研究[J].冰川冻土,1995,17(增刊):26-33
[24]Ma Wei.Review and prospect of the studies of ground freezing technology in China[J].Journal of Glaciology and Geocryology,2001,23(3):218-224
[25]付厚利.深厚表土层中冻结壁可靠性设计理论初探.建井技术,1999,20(4):32-34
[26]汪仁和,曹蓉斌.双排管冻结下冻结壁温度场形成特征的数值分析,2002,24(2):181-185
[27]崔广心.深土冻土力学-冻土力学发展的新领域.冰川冻土.1998,20(2):97-100
[28]胡省三,成玉琪.21世纪前期煤炭科技重点发展领域探讨.煤炭学报,2005,30(1):1-7
[29]崔广心.厚表土中地下工程特殊施工.陈明和,王建平等主编.矿山建设工程新进展(上册):中国矿业大学出版社,2005.8:1-3
[30]朱林楠.深井人工冻结壁温度场分析,冰川冻土,1981 4(3):1-4
[31]黄德发,邓文芳.冻结法凿井在河南煤矿立井施工中的应用.陈明和,王建平等主编.矿山建设工程新进展(上册):32-39.中国矿业大学出版社,2005.8
[32]佟凤健,朱德明.郭屯矿井建设若干技术问题研究.陈明和,王建平等主编.矿山建设工程新进展(上册):53-57.中国矿业大学出版社,2005.8
[33]卢相忠,宋世宏.丁集主副井冻结工程冻结壁分析.陈明和,王建平等主编.矿山建设工程新进展(上册):165-172.中国矿业大学出版社,2005.8
[34]刘明凯.丁集矿井深厚表土冻结技术.陈明和,王建平等主编.矿山建设工程新进展(上册):173-177.中国矿业大学出版社,2005.8
[35]李连华,王然芳.龙固矿副井深厚表土层冻结法施工.陈明和,王建平等主编.矿山建设工程新进展(上册):29-132.中国矿业大学出版社,2005.8,1
[36]李功洲,王长生.程村矿 主副井井筒深厚冲积层冻结法凿井技术研究.陈明和,王建平等主编.矿山建设工程新进展(上册):73-80.中国矿业大学出版社,2005.8
[37]林鸿苞,江新春,姚飞建.山东巨野矿区郭屯矿井特厚深表土复杂地层井简施工方法的选择[A].周兴旺,李明远.全国矿山建设学术会议论文选集(上册)[C].全国矿山建设学术会议论文选集(上册):94-99.中国矿业大学出版社,2004.8
[38]李述训.立井冻结法凿井工程中的热工计算.冰川冻土,1994,16(1)
[39]臧桂茂,谭杰.谈巨野煤田冻结法凿井技术问题.建井技术,2001(2)
[40]翁家杰.井巷特殊施工.北京:煤炭工业出版社.1991:4-5
[41]沈季良,崔云龙,王介峰.建井工程手册.北京:煤炭工业出版社.1986:60-317
[42]郭兰波等.竖井冻结壁温度场的有限元分析.中国矿业学院学报,1981,No.3
[43]洪伯潜.特殊凿井技术在我国的发展与前景.中国煤炭,2000,26(4):68-70
[44]李功洲等.深厚冲积层冻结凿井技术问题的探讨.周兴旺,李明远.全国矿山建设学术会议论文选集(上册)[C].全国矿山建设学术会议论文选集(上册):284-292:中国矿业大学出版社,2004.8
[45]苏立凡.对我国建设700m冻结立井的技术要素分析[A].王长生,苏立凡,张文.地层冻结工程技术和应用-中国地层冻结工程40年论文集[C].北京:煤炭工业出版社,1995
[46]Ma Wei(马巍).Review and Prospect ofthe Studies of Ground Freezing Technology in China(中国地层土冻结技术研究的回顾与展望).冰川冻土,2001,23(3):218-224
[47]苏立凡,徐兵壮.我国地层人工冻结工程技术的进展,应用和问题[A].周兴旺,张向东.全国矿山建设学术会议论文选集(上册)[C]:122-125.吉林:东北大学出版社,2003.8
[48]陈湘生.深冻结壁时空设计理论.岩土工程学报,1998,120(5):13-16
[49]虞相,王正廷,苏立凡.我国地层冻结技术的新发展-纪念冻结法凿井40周年[A].王长生,苏立凡,张文.地层冻结工程技术和应用-中国地层冻结工程40年论文集[C].北京:煤炭工业出版社,1995
[50]崔广心.关于冻结壁几个问题的探讨[A].王长生,苏立凡,张文.地层冻结工程技术和应用-中国地层冻结工程40年论文集[C].北京:煤炭工业出版社,1995
[51]姚志明,周世虎.立井特厚粘土层冻结施工技术.煤炭工程,2004(12):26-27
[52]李银桥等.深立井复杂冲积层冻结技术.全国矿山建设学术会议论文集,2004
[53]陈瑞杰,陈国栋,李述训等.人工地层冻结应用研究进展和展望.岩土工程学报.,2000,22(1):40-44
[54]陈文豹等.潘集矿区冻结壁平均温度及冻结孔布置圈径的探讨,煤炭学报,Nol,1982
[55]刘珣,梁鹏.城市地下工程中人工冻结法的防冻胀优化设计研究.武汉科技大学学报,2004,27(4):387-390
[56]陆卫国,郭圣澏,林稚华.广州地铁建设中冻结工法的应用分析.西部探矿工程.,2004(92):108-109
[57]郭晓江.冻结法在广州地铁二号线暗挖隧道中的应用.煤炭工程,2001(12):27-29
[58]周兴荣,王宗金等.冻结法在广州丫髻沙大桥桩基事故处理中的应用,2000.21(2):38-40
[59]程知言等.上海某冻结加固工程试验研究分析,地质与勘探,2002,38(3)
[60]程国栋,马巍.国际冻土工程研究进展-第五届冻十工程国际学术讨论会综述.冰川冻土,2003,25(3):303-308
[61]程国栋.冻土力学与工程的国际研究新进展.地球科学进展,2001,16(3):293-299
[62]朱正刚.寒区路基工程与多年冻土间相互作用问题研究进展.结构工程师,2004,20(6):43-48
[63]汪仁和,张瑞.多圈冻结管下冻土扩展速度的实测研究.陈明和,王建平等主编.矿山建设工程新进展(上册),2005.8:86-91
[64]周晓敏,张云利.竖井单圈布孔冻结壁形成平面温度场问题的数值试验研究.陈明和,王建平等主编.矿山建设工程新进展(上册),2005.8:104-110
[65]张开顺,李炳胜.井筒冻结壁径向位移规律计算研究.陈明和,王建平等主编.矿山建设工程新进展(上册),2005.8:116-120
[66]雷晓燕.岩土工程数值计算.中国铁道出版社.1999.11
[67]谢康和,周健.岩土工程有限元分析理论与应用.科学出版社.2003.5
[68]张铜生,张富德.简明有限元法及其工程应用.地震出版社.1990.8
[69]赵阳升等.冻结壁的温度场和弹塑性应力场的耦合分析.煤炭学报,1985.6
[70]朱林楠.深井人工冻结壁温度场分析.冻川冻土,1981.11
[71]张燕.预报冻结壁形成及其温度分布的有限元方法.煤炭学报,1979.2
[72]汤志斌等,冻结壁温度场的实测研究及在工程中的应用,地层冻结技术和应用--中国地层冻结工程40年论文集,1995.10
[73]赵阳升等.冻结壁的温度场和弹塑性应力场的耦合分析.煤炭学报,1985.6
[74]杨世铭编.传热学.人民教育出版社,1998.3
[75]俞昌铭编著.热传导及其数值分析.清华大学出版社,1982
[76]周晓敏,张云利.竖井单圈布孔冻结壁形成平面温度场问题的数值试验研究.矿山建设工程新进展(上册),2005.8.104-110
[77]李述训.人工冻结凿井中热学问题计算方法评价.世界冰川冻土,1982(3):23-32
[78]马思民,王金芝,李洪升.冻土温度场的参数辨识.冰川冻土.2002,24(3):299-303
[79]谢康和,周健.岩土工程有限元分析理论与应用.科学出版社.2003.5
[80]张铜生,张富德.简明有限元法及其工程应用.地震出版社.1990.8