用边界单元法求解巷道围岩的散热量.pdf

第29卷第4期 中国矿业大学学报 C01．29No．4 of 2000年7月 JournalChina of University Jul．2000 MiningTechnology 文章编号；1000-1964(2000)04—040304 用边界单元法求解巷道围岩的散热量 秦跃平，党海政，刘爱明 (中国矿业大学资源开发工程系，jE京100083) 摘要：将阐岩温度场的导热微分方程厦其单值性条件转化为无因次形式，将拳解巷道壁面温度或 困岩散热量的问题归结为求解不稳定换热准数K。，利用二维不稳定导热方程的时间相关基夺 解建立边界积分方程，用边界单元法计算风流与躅岩不稳定挟热准数．在此基础上编常3了计算机 程序，根据计算结果绘制了巷道不稳定挟热准数醚毕渥准教Bi及傅立叶准数Fo的变化曲线．对 于实际问题，可先计算该问题对应的fb和西，再按本文方法确定jo，最后可计算出巷遭圃岩散 热量． 关键词：边界单元法；巷道固岩散热；无因次 727．2 中图分类号：TD 文献标识码：A 巷道围岩散热是通风巷道内的主要热源，它对 流在其平均温度下长期通风的不稳定换热．对于风 于矿井风温预测、制定降温措施和计算掘进工作面 流温度的周期性变化可在此基础上进行修正m． 的制冷量都是至关重要的．过去计算这一部分散热 量时，常常借助于有限差分法和有限单元法，但这 两种方法运用到半无限体时，势必要人为地假定外 部边界，使得数学模型与实际有一定的偏差{而且 在具体操作中，其求解规模和数据准备量也较大． 本文采用边界单元法，将区域内部的积分问题转化 到边界上，然后再以问题的基本解建立边界积分方 程，对时间和边界离散进行数值计算。 图1拱形巷道圃岩温度坜币惹图 Sketch isothermin rock F垃．1 map。f surrounding 1 巷道围岩温度场的导热微分方程 这样，对于常年通风巷道，其围岩灏度场可看 由于井下通风巷道散热过程的复杂性，作如下 成是二维的，其导热微分方程及边值条件为“]； 简化； 1)巷道周围岩体是均质的．由于计算的目的 n(嘉+骞)一塞，4l孬+霜』2磊’ u’㈣ 是确定围岩散热量而不是温度场分布，各岩层的热 g=一^差I一口(“一tO， (2) ml 物理性质差别叉不大，取一定范围内各岩层物理参 r=ro 数的平均值，对计算结果不会产生太大影响． t=tOr一。。． (3) 2)巷道截面为翻形．尽管巷谨的形状有梯形、 t=to，r=0， (4) 拱形、矩形等多种形式，但在均质各寓同性岩层中， 式中；t为岩体温度，℃；t。为原始岩体温度，℃；k 其等温线近似为圆形，只在巷道周围很小的范围内 为巷道壁面岩石温度，