土木译文.docVIP

严寒地区隧道围岩冻融状况分析的导热与对流换热模型—东北走向. 由于大坂山地区隧道施工现场平均气温为负温的时间每年约长8个月，加之施工时间持续数年，围岩在施土过程中己经预冷，所以隧道开通运营后，洞内气体流动的形态主要由进出口的主导风速所确定，而受洞内围岩地温与洞外气温的温度压差的影响较小；冬季祁连山区盛行西北风，气流将从隧道出曰流向进口端，夏季虽然祁连山区盛行东偏南风，但考虑到洞口两端气压差、温度压差以及进出口地形等因素，洞内气流仍将由出口北端流向进口端.另外，由于现场年平均风速不大，可以认为洞内气体将以层流为主 基于以上基本情况，我们将隧道简化成圆筒，并认为气流、温度等关十隧道中心线轴对称，忽略气体温度的变化对其流速的影响，可有如下的方程: 其中t为时间，x为轴向坐标，r为径向坐标;U, V分别为轴向和径向速度，T为温度，P为有效压力(即空气压力与空气密度之比少，V为空气运动粘性系数，a为空气的导温系数，L为隧道长度，R为隧道的当量半径，D为时间长度, 分别为围岩的冻、融区域. ,分别为冻、融状态下的热传导系数，,分别为冻、融状态下的体积热容量， X=(x,r) , 为冻、融相变界面,To为岩石冻结临界温度(这里具体计算时取To=-0.10)，为水的相变潜热. 2 求解过程 由方程(1)知，围岩的温度的高低不影响气体的流动速度，所以我们可先解出速度，再解温度. 2.1 连续性方程和动量方程的求解 由于方程((1)的前3个方程不是相互独立的，通过将动量方程分别对x和r求导，经整理 化简，我们得到关于压力P的如下椭圆型方程: 于是，对方程(1)中的连续性方程和动量方程的求解，我们按如下步骤进行: (1)设定速度,; ( 2)将,代入方程并求解，得 (3)联立方程(1)的第一个和第二个方程，解得一组解,; (4)联立方程((1)的第一个和第三个方程，解得一组解,; (5)对((3) ,(4)得到的速度进行动量平均,得新的,返回(2) ; (6)按上述方法进行迭代，直到前后两次的速度值之差足够小.以,,作为本时段 的解,下一时段求解时以此作为迭代初值. 2. 2 能量方程的整体解法 如前所述，围岩与空气的温度场相互作用，壁面既是气体温度场的边界，又是固体温度场的边界，壁面的温度值难以确定，我们无法分别独立地求解隧道内的气体温度场和围岩温度场.为克服这一困难，我们利用在洞壁表面上，固体温度等于气体温度这一事实，把隧道内气体的温度和围岩内固体的温度放在一起求解，这样壁面温度将作为末知量被解出来.只是需要注意两点:解流体温度场时不考虑相变和解固体温度时没有对流项;在洞壁表面上方程系数的光滑化.另外，带相变的温度场的算法与文献[3]相同. 2. 3热参数及初边值的确定 热参数的确定方法: 用p=1013.25-0.1088H计算出海拔高度为H的隧道现场的大气 压强，再由计算出现场空气密度，其中T为现场大气的年平均绝对温度，G为空气的气体常数.记定压比热为,导热系数为，空气的动力粘性系数为.按 和 计算空气的导温系数和运动粘性系数.围岩的热物理参数则由现场采样测定. 初边值的确定方法:洞曰风速取为现场观测的各月平均风速.取卞导风进曰的相对有效 气压为0，主导风出口的气压则取为，这里k为隧道内的沿程阻力系数，L为隧道长度，d为隧道端面的当量直径，为进口端面轴向平均速度.进出口气温年变化规律由现场观测资料，用正弦曲线拟合，围岩内计算区域的边界按现场多年冻土下限和地热梯度确定出适当的温度值或温度梯度. 3 计算实例 按以上所述的模型及计算方法，我们对大兴安岭西罗奇2号隧道内气温随洞曰外气温变化的规律进行了模拟计算验证，所得结果与实测值[6]相比较,基本规律一致. 西罗奇2号隧道是位十东北嫩林线的一座非多年冻土单线铁路隧道，全长1160 m ,隧道 近西北一东南向，高洞口位于西北向，冬季隧道主导风向为西北风.洞口海拔高度约为700 m , 月平均最高风速约为3m/s,最低风速约为1.7m/s.根据现场观测资料，我们将进出口气温拟 合为年平均分别为-5和-6.4，年变化振幅分别为18.9和17.6的正弦曲线.隧道的当量直径为5.8 m,沿程阻力系数取为0.025.由于围岩的热物理参数对计算洞内气温的影响 远比洞口的风速、压力及气温的影响小得多，我们这里参考使用了大坂山隧道的资料. 图1给出了洞口及洞内年平均气温的计算值与观测值比较的情况，从进口到出口，两值之差都小于0.2. 图2给出了洞内 (距进出口l00m以上)月平均气温的计算值与观测值比较的情况，可以看出温度变化的基本规律完全一致，造成两值之差的主要原因是洞口气温年变化规律之正弦曲线的