混凝土水管冷却计算的一种复合算法.pdfVIP

北京力学会第18届学术年会论文集：计算力学 混凝土水管冷却计算的一种复合算法 邱永荣张国新 (中国水利水电科学研究院结构材料研究所(研)，100038) 摘要：提出一种新型的计算冷却水管效应的有限元与解析法结合的复合算法。该算法无需 通过剖分密集的网格来计算水管壁的温度梯度，也无需依靠反复的迭代来计算水温的沿程 变化，相比精细有限元算法，复合算法具有明显的速度优势。 关键词：温控计算，水管冷却，复合算法 ．1￡．—j- 一、丹{J吾 水管冷却计算一直是混凝土温度仿真计算中的难点，如果用热流耦合方法计算，需 要耗费大量计算时间。朱伯芳【lj提出过完整的精细计算算法，但该算法需要在水管周围 剖分密集的单元网格和反复的迭代，会使计算耗费大量的时间。国内一些学者针对该问 题也提出了各种算法，比如子结构方法12J，复合单元法pJ，多重有限元嘲格法14]等。但 这些方法都没有解决迭代计算的问题，计算时间仍然较慢。本文提出一种新型的计算冷 却水管效应的有限元与解析法结合的复合算法。该算法无需通过剖分密集的网格来计算 水管壁的温度梯度，也无需依靠反复的迭代来计算水温的沿程变化。 二、计算原理 2．1基本原理 复合算法的基本原理为：对于每步的温度场计算，都分两个子步计算，子步一，先 用有限法计算不考虑水管冷却的温度场，该子步的初始温度为上一步的温度场，考虑边 界散热和水化热放热影响；子步二， 初始温度为0，边界绝热，再按图1 模型逐个截面计算每个截面的等效 水温和温度场，截面的温度场计算可 根据解析公式l计算。最后该计算步 的温度场即为两子步温度场的叠加。 一．!L (1) re，，(，)=瓦，，(1一∑以(，．)P‰‘) ”1 图1水管计算模型 2．2算例 底板长30m，宽15m，厚lm，2009年9月15号浇筑。基础温度24℃，浇筑温度 15C，水管进水口水温15。C。水管布置按1．5×lm布置，布置在底板厚度方向中间。水 管为钢管，管内直径为2．0cm。通水流量1．5mJ／h，24小时换一次水流方向，通水时间 5d。混凝土和基础的热学参数见表l。 图2～图5显示，运用复合算法计算的温度结果与严格的精细算法相差较小，具有 较高的精度。表2显示，复合算法相比精细算法，计算效率有大幅提高。 Ⅵ．7I 1 “京力学套gg目学术年喜论￡采计算力学 盥2精细算法的Ⅺ度场分布 嘲3复音算法的温度场分布 ”～‘“”8”i。：乳”…5”……“’“…’一i：ii，…8 圈4水管进口淞鹰卅比 剧5涩凝土中心温瞳对比 重：盐墨塑兰型些 *』c§ #^＆ m∞n*¨自】㈧ ‘#* 0E口十M(Ⅻn) #日”# 1055∞I啪％ 12593 100 209 i☆H《{I*3Ⅲ】62937 【00 4895 ＆t￡^lE#AM，∞*《**nI■* Intel【R)co删口cPU6，∞@l86№，】006BlI自# 三、结论 复台算法计算结果与辅细算法基木‘致，误差较小。枉计算效唪方而．复合算洁仃 明显优势。如果该套算法应用于大型结构(如高拱坝)的温度计算中，可大幅¨省计算 时间。 参考文献 ：i??；i!‰；§；：i：：：!；镕i：i§籍i：勰嚣喘端。。。 ；凇黼鼎焉黠黼皴臻黼i蹴嚣麓黧i：慑j躲