冻结壁温度场水分场耦合分析研究.pdfVIP

冻结壁温度场水分场耦合分析研究 杨玉贵，高峰，刘冠男，李涛 中国矿业大学理学院，江苏徐州 (221008) E-mail ：yyg323@163.com 摘 要：人工冻结法施工技术是深厚表土矿井开挖的一种重要的技术，为了更好的应用这 种技术，本文利用ANSYS 二次开发功能，编写了考虑水分场对温度场参数影响的APDL 以及 UPFS 二次开发程序，通过该程序分析研究了冻结过程中温度场、水分场分布及扩展规律，计 算实例分析验证了该程序的正确性。 关键词：ANSYS；数值模拟；温度场；耦合 竖直方向上的传热比较弱，因此我们可以看 0 引言 作是平面问题。在冻结壁交圈之前，温度场 在矿井冻结法施工工程中，冻结壁的温 是依每一个冻结管为中心的近似轴对称问 度场分布一直是工程师最为关心的问题。因 [2] 题，考虑潜热释放的热传导方程如下 ： 为某时刻的冻结壁温度场分布决定着此时 2 ∂T ∂T 1∂T ∂θ C ( =+ )+L i 的冻结壁厚度，而冻结壁中某点的温度值又 ρ λ 2 ρ （1） ∂t ∂r r ∂r ∂T 决定着此处的冻结壁强度。冻结壁厚度和强 式中 ρ—粘土的密度； C —粘土的比热； 度是衡量冻结质量的重要指标，工程师往往 V λ—导热系数；L —潜热；θ —冰体积含量； 据此把握开挖时机、给出开挖方案。较一般 i 热传导问题不同的是，冻结会使土体中各组 冻土中的水流方程满足达西定律，并根 [3] 分的含量发生改变，水分在冻结时会逐渐凝 据连续性方程得到 ： 固成冰，水和冰的热容、导热系数等物性参 ∂ϕ c(θ) ∇[K θ ∇ϕ] ( ) 数是不一样的，因此冻土的整体物性参数也 ∂t 要随着温度的改变而改变； （2 ） 潜热并不只发生在一个狭窄的区域内，一般 式中，θ是体积含水率，K θ 是不饱和渗 ( ) 在负温区域内都会有潜热释放。归结起来， 水系数， c(