水管冷却混凝土温度场的精确计算与应用.docVIP

水管冷却混凝土温度场的精确计算与应用 第23卷第2期 2005年4月 水电能源科学 WaterReSOLIrcesandPower V01.23No.2 文章编号:1000—7709(2005)02—0073.04 水管冷却混凝土温度场的精确计算与应用 林志祥朱岳明 (河海大学水利水电工程学院,江苏南京210098) 摘要:简要介绍了基于水管与混凝土之闻热量交换的能量守恒原理,水管冷却混凝土温度场的精确算法及有 限元数值计算方法.算例和试验验证,该迭代方法具有很高的计算效率和精度.应用该精确算法为几个大坝 仿真分析,均取得了很好的效果,证明该算法具有良好的工程应用前景. 关键词:混凝土温度场;水管冷却;有限元;仿真计算;迭代求解;温控防裂 中图分类号:TV315文献标志码:A 水管冷却技术是大体积混凝土温控防裂最有 效和最常用的措施之一.然而,由于混凝土温度 场自身的非线性,内部数目众多的冷却水管中开 放式流动水体和混凝土的热交换以及浇筑层面与 表面的散热,使整个混凝土体的温度场变得非常 复杂,成为极难准确求解的复杂的三维不稳定温 度场问题.文献[1]曾提出把冷却水管看成负热 源,建立考虑水管冷却效果的等效热传导方程,在 平均意义上考虑水管冷却的效果.文献[2,3]提 出一种从平面扩充列空间的近似解法,奠定了水 管问题精确求解的基础.文献[4]在该方法基础 上,基于水管与混凝土之间热量交换的能量守恒 原理提出了理论上正确严密的新算法,并应用到 数值计算中提高计算精度.针对实际工程中大量 的蛇形冷却水管,通过水管截面的形心坐标和水 管半径方便地实现沿程搜索水管的边界单元,再 进行曲面积分,就可以严格地模拟水管的实际走 向来计算管内沿程水温的增量,避免以往裁弯取 直所引起的计算误差. 1水管冷却温度场的精确计算 1.1温度场的基本理论 众所周知,根据变分原理,三维非稳定温度场 的求解可取如下泛函: IcT=『『jR .{丢[()+()+(aT)]+ 1(aT一霎)T)dxdydz+『』(一T.)T一 (1) 式中,R为计算域;T为温度;a为导温系数;为绝 热温升;为导热系数;卢为表面散热系数;T.为环 境温度;t为时间;r为龄期;I1.为第三类边界域. 由泛函的驻值条件I/T一0和时间域的向 后差分,可得温度场求解的递推方程组: (H+)计-一Tn+F计-一0(2) 由已知上一时刻的T,根据式(2)计算可得 下一时刻的T+.. 1.2三维水管冷却有限元算法 水管冷却的解析解只适用于具有简单边界的 温度问题,对于具有复杂边界的水管冷却温度问 题,一般采用有限元数值解法].如图1所示, 设A为水管周围的混凝土单元,在单元A范围 内水温的增量为△T,则单位时间内水流在单元 A内所吸收的热量为 △Q一clDqAT+ClDVAT(3) 式中,c为水的比热;ID为水的密度;q为单位时 间内水的流量;V为单元A范围内水管内静止 水的体积;△T为单位时间内静止水上升的温度. 第2相比第1项小得多,其热量不到第1项的万 分之五,可忽略不计.在单元A内,混凝土向水 流放出的热量为 收稿日期:2005-01-17 基金项目:水利部科技创新重点项目资助(SCXC2003一10);河海大学院士学科发展基金资(HHUYS003) 作者简介:林志祥(1978一),男,博士研究生,研究方向为水工结构工程温控防裂,E-mail:lzx7808t~163.tom ? 74?水电能源科学2005年 圄…图1冷却水管及单元划分示意图 △Q=一f1ds(4) JSO,(,,f 式中,积分是沿着水管外缘进行的.由热量平衡, △Q一△Q,得到单元A内水温的增量为 一一 j(筹)㈣ 因此,经过单元A后的水温为 T. 斗1一T.i-F△T.f(6) 各点温度在计算前未知,因而式(5)中的温度 梯度(aT~an)也是未知的,须迭代求解.先假定 水管沿程各点的水温都等于进口水温,加上其他 的初始条件和边界条件,由式(2)求出混凝土的整 体温度场,按式(5)和(6)求出水管沿程各点的第 一 次近似水温T;然后以T.in..作为管沿程各点 的初始水温,其他的初始条件和边界条件不变,求 出混凝土的整体温度场,同样按式(5)和(6)求出 水管沿程各点的近似水温T;重复上述过程,比 较前后两次的水温差值,直到满足收敛准则: 1TlaxIT一TI≤e(e为一指定的小数)(7) 求解非稳定温度场时,第_『+1时段的初始水 温可采用第_『时段迭代后的水温,这样将减少迭 代步数,大大节省计算工作量. 上述算法是我们新近的研究成果,利用其进 行有水管冷却作用的混凝土试块的室内温升试验 和淮河人海水道二河新泄洪闸现场试验,对其准 确性和计算效果进行了验证.结果证明,该算法 不但理论上严密,而且迭代计算的