水管冷却混凝土温度场计算三维p型有限元法.pdfVIP

第12@ 第2期 中国水利水 电科学研究院学报 Vo1．12No．2 2014年6月 JournalofChinaInstituteofWaterResourcesandHydropowerResearch June,2014 文章编号：1672—3031(2014)02—0122—08 水管冷却混凝土温度场计算的三维P型有限元法 王海波，王仁坤 (中国水电顾问集团 成都勘测设计研究院有限公司，四川 成都 610072) 摘要 ：引入阶谱与 自适应技术 ，将P型有限元法应用于水管冷却混凝土温度场 的仿真计算 ，通过对单元升阶谱 ， 以扩大容许函数空间来提高数值解精度 ，并编制 了P型 自适应有限元计算程序。针对同一计算模型，分别采用不 同网格密度的常规有限元法和P型有限元法进行了水管冷却混凝土温度场的计算分析和对 比研究。结果表明，P 型有限元法的计算结果在网格稀疏的情况下仍然能够达到较高的精度，特征点温度时程线和温度场等值线与精细 算法值吻合得较好 ，而常规有限元法则误差相对较大。同时 ，采用P型有限元法进行水管冷却混凝土温度场模 拟 ，可以在保证计算精度的前提下减少单元数量 ，从而降低前处理的难度 、工作量 以及计算的时间成本 。 关键词 ：P型有限元法 ；水管冷却；温度场 ；混凝土 中图分类号 ：TV315 文献标识码 ：A doi：10．13244~．cnki．jiwhr．2014．02．002 1 研究背景 水管冷却是混凝土坝施工中控制温度 、防止裂缝的重要方法之一 。自20世纪30年代在美 国胡佛 拱坝施工中首次成功应用 以来 ，已在大体积混凝土施工 中得到广泛应用 ，但水管冷却效果的计算一 直是混凝土温度仿真计算 的难点。美国垦务局利用分离变量法得到了无热源平面问题的严格解答和 空间问题的近似解答…，由于没有考虑内热源，该解答只适用于混凝土二期水管冷却温度场的计算 。 为了解决该难题 ，朱伯芳院士 提出了混凝土水管冷却效果的等效算法 ，该方法将冷却水管看成热 汇 ，从能量等效的角度出发 ，在平均意义上反映水管的冷却效果 ，不考虑坝体 内水管 的具体布置型 式 ，因此无需划分密集的有限元网格就可 以计算考虑水管冷却 的混凝土温度场 ，使冷却水管问题得 到极大的简化 ，有效地避免 了计算 网格过密 、计算工作量过大的问题 ，在计算机软硬件条件仍不具 备的情况下，该方法不失为一种有效的办法 ，更符合工程实际应用的需要。然而 ，水管冷却效果其 实是一个 由内向外逐渐扩散 的过程 ，等效算法不能真实反映沿程水温变化对混凝土温度 的影响，对 水管周边混凝土的温度梯度无法模拟 ，无法反映局部 区域 的应力变化过程。 水管离散算法 在计算 中考虑了密集而走 向复杂的水管，能够精确地模拟水管在混凝土坝中的 具体位置及布置型式 ，可以较好地体现水管周边混凝土温度 由近及远逐渐扩散的过程 ，以及水管周 围区域中温度梯度所导致的应力变化差异 。但为了保证计算精度 ，在水管周围必须划分密集的有限 元网格 ，对于混凝土高坝 ，浇筑层数多、施工工期长 ，温度场的仿真分析需要很长的计算时间，计 算效率低 。因此 ，精细算法 目前只能用于个别特别重要的、需要反映局部区域温度和应力变化过程 的工程结构 ，其扩大计算能力需从硬件进步和算法改进方面进行研究。针对此问题 ，子结构技术 、 复合单元算法 等引入到水管离散算法中后 ，在一定程度上降低了算法对计算机的要求 ，但是需要使 用节点凝聚技术或一些必要的假定，其有限元格式 比较复杂。 有限单元法作为一种数值模拟方法，对于一个给定的问题，从逼近真实解的途径分类 ，可以分 为三种 ：(1)h型 。即传统的有限元方法 ，通过减小单元尺寸 ，增加单元个数和节点个数 ，来提高 有限元解 的精度 ；(2)P型 。不改变网格 ，通过增加基 函数多项式 的阶数 ，扩大容许函数空间，从 收稿 日期 ：2013—11-叭 基金项 目：国家 自然科学基金资助项 目 作者简介 ：王海波(1984一)，男，湖北人，博士，主要从事混凝土温控防裂及试验研究。E-mail：whblzp@gmail．com 一 122 —