模拟混凝土水管冷却效应的直接算法 Direct algorithm for simulating cooling effect of water pipes in concrete.pdfVIP

水 利 学 报 2009年7月 SHUILIXUEBAO 第加卷第7期 文章编号：0559．9350(2009)07．0892．05 模拟混凝土水管冷却效应的直接算法 刘晓青，李同春，韩勃 (河海大学水利水电学院，江苏南京210098) 摘要：本文提出了模拟混凝土水管冷却效应的直接算法，该法根据混凝土和冷却水管间热量交换平衡条件，建立同 时以混凝土内部温度和冷却水管结点温度为未知量的有限元直接求解方程，克服了常规迭代解法需计算水管周边 混凝土温度梯度的缺点。经典型算例验证，该方法简单可行，可应用于实际工程中冷却水管的直接模拟。 关键词：混凝土结构；冷却水管；直接算法 中图分类号：TV315 文献标识码：A 1研究背景 水管冷却是混凝土温控的重要措施…。在混凝土大坝的温度场仿真计算中冷却水管温度效应的模 拟总体上可以分为三大类：第一类是文献[2—5]的等效算法；第二类是文献[6]提出的水管冷却理论解 与有限元结合的计算方法；第三类是将冷却水管直接模拟的计算方法。虽然现阶段仍以文献[2—5]提 出的等效算法为主，但随着计算机性能的不断发展，直接对冷却水管进行模拟的方法越来越得到重视。 在对水管进行直接模拟过程中，对于冷却水管水温的计算，文献[5]提出的迭代算法已被大多数研 究者采用一‘17|，其基本原理是利用混凝土与冷却水管交界面上结点的温度相同条件，假定水管沿程各结 点水温，根据有限元方程求解混凝土内各结点温度。在此基础上求解水管表面各点的温度梯度，然后进 行表面积分求解混凝土流入水管中的热量，进而由进VI水温推求沿程各点水温。上述步骤重复进行，直 至收敛。为了提高温度梯度的计算精度，文献[5]提出了如下两个方法：(1)在孔口附近采用密集的计算 网格；(2)孔口附近的结点沿半径方向布置。 本文根据混凝土和冷却水管间热量交换平衡条件，建立同时以混凝土内部温度和冷却水管结点温 度为未知量的有限元直接求解方程，利用该法无需计算水管周边的温度梯度，以期简化冷却水管模拟中 的计算分析过程。 2混凝土不稳定温度场求解的有限元方法 在区域R内，混凝土温度满足热传导方程㈨ T 。 2 ，一V xOTO磐+雾+磊Oe一一00(a’，2。az21＼_+T 口、 Or Or)=。 并满足以下4种可能边界条件之一： (1)表面温度已知边界条件：T(r)=f(r)5 收稿日期：2008．05．27 基金项目：国家自然科学基金资助项目；国家科技支撑计划(2006BAB04A03) 作者简介：刘晓青(1965一)，女，四川成都人，博士，副教授，主要从事有限元在水工结构工程中的实用化研究 E·mail：Ixqhhu@163．coin _。——892·-—— 万方数据 (2)表面热流量是时间的已知函数：一A磊01=厂(r)； (3)与空气接触的边界条件：一A丽OT=卢(r—L)； (4)两种不同固体接触边界条件：耻¨。等咄等； 其中：T为温度；r为时间；0为绝热温升；卢为表面放热系数；A为导热系数；口为导温系数；L为 气温。 对式(1)运用Galerkin方法进行空间离散，可得到有限元方程 [H]…“s]{筹)=T1{Q，…Qo} (2) ’ — 式中：[日]=』nVⅣ’VJn Ndl2；㈥=』n』＼r7扣溉协}-A．『，∥O-dTdF；旧}=．f。1口Nt／,Jn 口 Jr d Jn口 1薯帆其中Ⅳdr 为温度场插值的形函数；0为混凝土温度场计算区域；r为与该区域相对应的外边界。 在考虑冷却水效果问题中主要采用两类边界r。和乃，前者