一种大体积混凝土水管冷却温度场的模拟方法
【专利摘要】本发明公开了一种大体积混凝土水管冷却温度场的模拟方法，其特征是，包括以下步骤：步骤一，划分混凝土结构与地基网格；步骤二，划分水管网格；步骤三，创建复合单元；步骤四，定义材料属性；步骤五，创建水管编号；步骤六，创建冷却水信息文件；步骤七，计算浇筑时地基温度场；步骤八，编写自定义单元子程序计算温度场：将地基温度场作为初始温度场，计算得到各浇筑阶段的温度场。本发明能够有效地模拟冷却水管对大体积混凝土结构的整体冷却效果，反应水管周围的温度场梯度变化；划分网格时不需考虑水管布置，求解混凝土节点温度的同时可以得到水管的节点温度，计算时不需要添加额外的迭代步来求解水管温度，具有很高的前处理效率和计算效率。
【专利说明】
-种大体积混凝±水管冷却溫度场的模拟方法
技术领域
[0001] 本发明设及一种溫度场模拟方法，具体设及一种基于ABA卵S的大体积混凝±水管 冷却溫度场的模拟方法，属于建筑结构分析技术领域。
【背景技术】
[0002] 溫度应力对大体积混凝±早期的安全性有重要影响，是大体积混凝±结构早期裂 缝产生的主要原因之一。在大体积混凝±结构中布置冷却水管，能够有效地降低混凝±结 构溫升，从而降低溫度应力。工程计算中，一般使用有限元法模拟混凝±结构早期的溫度应 力场。但是，应用有限元法准确模拟包含冷却水管的大体积混凝±结构的溫度应力场是一 项很复杂的工作，有诸多困难。首先，大体积混凝±结构模型复杂，一般采用分期诱筑，有数 十个至数百个诱筑期;通水冷却时，每隔一到两天就要改变通水方向，运成倍增加了分析步 的个数。此外，在冷却水管附近，溫度梯度变化剧烈，容易在水管周围产生较大的拉应力，为 保证混凝±结构的安全，有必要对水管周围的溫度场进行准确的模拟。最后，在有限元前处 理过程中，大体积混凝±结构模型的复杂，有限元网格划分工作量大，而在划分混凝±网格 的同时需要考虑蛇形布置的冷却水管，运大大增加了前处理工作的工作量与难度。
[0003] 目前模拟含有冷却水管的大体积混凝±结构溫度场的方法总体上可W分为不考 虑水管周围溫度梯度的等效算法和考虑水管周围溫度梯度的精细算法。朱伯芳院±在《大 体积混凝±非金属水管冷却的降溫计算K《水利水电技术》，1997,6:28,朱伯芳）中提出了 冷却水管的等效算法，将冷却水管看作热汇，不考虑水管周围的溫度梯度变化，推导出了同 时考虑水管冷却效果和混凝±绝热溫升的混凝±等效溫升公式。由于该方法求得的溫度场 与实际溫度场在整体上差别不大，同时计算十分简便，因此得到了广泛的应用。但是由于该 方法未考虑水管周围的溫度梯度，不能准确计算水管周围的溫度应力。
[0004] 精细算法根据水管的布置，将水管作为线单元或管状的实体单元，与混凝±单元 共用节点，在模拟水管对混凝±整体冷却效果的同时，能够准确地模拟水管周围的溫度变 化，为计算水管周围的溫度应力提供依据。在运些算法中，水管单元作为实体网格，需要与 混凝±网格共节点。因此在划分网格时必须要考虑水管的布置，加大了前处理的工作量。
[0005] 为减少使用精细算法计算冷却水管溫度场的工作量，陈胜宏等人将复合单元应用 于大体积混凝±的溫度场计算当中，把水管单元看作混凝±单元中的虚单元，在划分混凝 ±网格时不需要考虑水管的布置，减少了前处理的工作量。计算时，初始水管溫度为假设水 溫，并且水管溫度和混凝±溫度不能同时求解，因此每一个分析步，都需要迭代来完成溫度 场的求解。运样不仅成倍增加了计算时长，而且由于增加了额外的迭代计算，一般商用软件 并没有提供修改迭代步的接口，因此运种复合单元法很难在商用有限元软件中实现，运也 就意味着如果使用复合单元法求解冷却水管溫度场，需要耗费大量的时间编写一整套有限 元程序，因此很难在推广到工程应用中。
【发明内容】
[0006] 本发明的目的在于克服现有技术中的不足，提供了一种大体积混凝±水管冷却溫 度场的模拟方法，解决了现有技术求解过程中前处理耗时多和计算效率低的技术问题。
[0007] 为解决上述技术问题，本发明提供了一种大体积混凝±水管冷却溫度场的模拟方 法，其特征是，包括W下步骤：
[0008] 步骤一，划分混凝±结构与地基网格:将混凝±与地基实体模型划分成混凝±单 元和地基单元，单元类型为二次六面体单元，每个单元包含20个节点和27个积分点，并对节 点进行编号；
[0009] 步骤二，划分水管网格:水管被混凝±单元轮廓切分为若干段，每一段记为一个水 管辅助单元，每个单元包含2节点和2个积分点，并对节点进行编号，包含水管辅助单元的混 凝±单元记为混凝±辅助单元；
[0010] 步骤=，创建复合单元:将每对水管辅助单元和混凝±辅助单元创建新的复合单 元，并对复合单元的节点进行编号；
[0011] 步骤四，定义材料属性:定义混凝±单元、地基单元的材料属性，水管辅助单元材 料属性为0，混凝上辅助单元的材料属性为混凝上属性；
[0012] 步骤五，创建水管编号:将属于同一根水管的水管辅助单元创建单元集合，集合名 称为水管编号；
[0013] 步骤六，创建冷却水信息文件:将所有诱筑阶段的冷却水信息按诱筑阶段划分来 创建文件；
[0014] 步骤屯，计算诱筑时地基的溫度场:将气溫作为地基表面的第=类边界条件，将当 地年平均气溫作为地基的初始溫度，计算得到诱筑时地基的溫度场；
[0015] 步骤八，创建自定义单元子程序计算各阶段的溫度场:将地基溫度场作为初始溫 度场，按照诱筑顺序逐层激活混凝±辅助单元、水管辅助单元、混凝±单元和复合单元，计 算得到各诱筑阶段的溫度场。
[0016]进一步的，在所述步骤=中，复合单元为用户自定义的单元，包含22个节点和29个 积分点，其中1~20节点编号对应混凝±辅助单元的节点，21及22节点编号对应水管辅助单 元的节点。
[0017] 进一步的，在所述步骤四中，材料属性包括材料的导热系数、比热容和密度。
[0018] 进一步的，在所述步骤六中，冷却水信息包括水管编号、该条信息的起始时间、是 否通水、通水方向是否逆向、水管内径、流速、水管厚度、水管导热系数、水管入口水溫。
[0019] 进一步的，在所述步骤八中，自定义单元子程序中复合单元的雅克比矩阵AMATRX 和残余向量R服的计算公式为：
[0020]
[0021]
[0022]
[0023]
[0024]
[0025]
[0026]
[0027]
[002引
[0029] 其中，TB为计算时间段末的节点溫度，调用自定义单元子程序时作为参数传入。JI 为圆周率，d为水管直径，h为水管的散热系数，1为水管长度，0^。为二次六面体单元单元高 斯积分的第i个积分点的权重，N。为混凝±子单元形函数，鸣为线性杆单元高斯积分的第i 个积分点的权重，化0为水管入口在混凝上子单元局部坐标下的形函数矩阵，化1为水管出口 在混凝±子单元局部坐标下的形函数矩阵，V为冷却水流速，Pw为水密度，Cw为水的比热容， 4为系数，用于调整求解瞬态溫度场收敛准则中水管溫度的容许误差。
[0030] 与现有技术相比，本发明所达到的有益效果是:本发明能够有效地模拟冷却水管 对大体积混凝±结构的整体冷却效果，并且可W反应水管周围的溫度场梯度变化;划分网 格时不需考虑水管布置，求解混凝±节点溫度的同时可W得到水管的节点溫度，计算时不 需要添加额外的迭代步来求解水管溫度，具有很高的前处理效率和计算效率。
【附图说明】
[0031 ]图1是本发明复合单元模型的结构示意图；
[0032] 图2是本发明实施例重力巧的几何模型图；
[0033] 图3是本发明实施例第二层诱筑块诱筑两天后溫度场云图；
[0034] 图4是本发明实施例第二层诱筑块诱筑两天后水管溫度变化；
[0035] 图5是本发明实施例第一年冬天溫度场云图。
【具体实施方式】
[0036] 下面结合附图对本发明作进一步描述。W下实施例仅用于更加清楚地说明本发明 的技术方案，而不能W此来限制本发明的保护范围。
[0037] 本发明提供的一种大体积混凝±水管冷却溫度场的模拟方法，利用大型非线性有 限元软件ABA卵S对溫度场进行模拟计算，在S维建模软件中建立S维的混凝±结构与地基 的实体模型，通过将实体模型划分为有限元，并定义复合单元，建立同时求解混凝±节点溫 度与水管节点溫度的有限元支配方程，能够有效地模拟冷却水管对大体积混凝±结构的整 体冷却效果;划分网格时不需要考虑水管布置，同时求解混凝±溫度和水管溫度，不需要增 加额外的迭代步求解水管沿程水溫，减小了计算量，大大地提高了计算效率。
[0038] 如图1所示，本发明的一种大体积混凝±水管冷却溫度场的模拟方法，其特征是， 包括W下步骤：
[0039] 步骤一，划分混凝±与地基网格:根据设计图建立结构和地基的=维实体模型，将 模型划分成二次六面体单元，每个二次六面体单元包含20个节点和27个积分点，并对节点 进行编号；
[0040]步骤二，划分水管网格:水管被混凝±单元轮廓切分为若干段，每一段记为一个水 管辅助单元，每个单元包含2个节点和2个积分点，并对节点进行编号，包含水管辅助单元的 混凝±单元记为混凝±辅助单元；
[0041 ]步骤=，创建复合单元:保留辅助单元，将每对水管辅助单元和混凝±辅助单元创 建新的复合单元，复合单元的单元类型为用户自定义单元UEL
[0042] 复合单元由一个混凝±辅助单元和一个水管辅助单元组成，其单元控制方程反应 了水管对混凝±的降溫作用和因水管从混凝±中吸收热量造成的沿程水溫升高。混凝±辅 助单元包含20个节点，27个积分点，水管辅助单元包含2个节点，2个积分点，复合单元共有 22个节点和29个积分点，节点编号时，1~20节点编号对应混凝±辅助单元的节点，21及22 节点编号对应水管辅助单元的节点。由于辅助单元与复合单元公用节点，在辅助单元上施 加面边界条件等同于在复合单元上时间面边界条件，显示辅助单元的节点溫度也就是显示 复合单元的节点溫度。复合单元在ABA卵S中的单元类型为用户自定义单元，而ABAQUS中不 支持对用户自定义单元时间面边界条件、显示节点溫度，因此借助辅助单元来完成创建表 面集合施加面边界条件、显示节点溫度的功能；
[0043] 步骤四，定义材料属性:定义混凝±单元、地基单元的材料属性，水管辅助单元仅 用于施加面边界条件、显示节点溫度，因此其材料属性为0,而混凝±辅助单元还用于计算 单元的比热容矩阵和热传导矩阵，因此材料属性为混凝±属性；
[0044] 材料属性包括材料的导热系数、比热容和密度；
[0045] 步骤五，创建水管编号:将属于同一根水管的水管辅助单元创建单元集合，集合名 称为水管编号；
[0046] 步骤六，创建冷却水信息文件:将所有诱筑阶段的冷却水信息按诱筑阶段划分来 创建文件；
[0047] 冷却水信息包括水管编号、该条信息的起始时间、是否通水、通水方向是否逆向、 水管内径、流速、水管厚度、水管导热系数和水管入口水溫；
[0048] 步骤屯，根据当地的天气预报，将气溫作为地基表面的第S类边界条件，在ABA卵S 的"interaction"模块中使用生死单元（"model change")功能杀死除地基W外的其他单 元，在"step"模块中创建瞬态溫度场分析步，将当地年平均气溫作为地基的初始溫度，计算 诱筑前两年至诱筑前一天的计算地基溫度场，所得的最后一天溫度场即诱筑时地基的溫度 场;其中计算地基溫度场的前推时间一般为诱筑前1~2年，为了能够更精准推算出诱筑时 地基的初始溫度场，运里选择2年；
[0049] 步骤八，创建自定义单元子程序计算各个诱筑阶段的溫度场:将地基溫度场作为 初始溫度场，在ABA卵S的"interaction"模块中使用生死单元（"model change")功能，按照 施工过程中诱筑每个诱筑层的日期，逐步激活每个诱筑层所包含的混凝±辅助单元、水管 辅助单元混凝±单元和复合单元，在"job"模块创建分析作业，通过"write i吨Ut"功能输 出作业的inp文件。
[0050] 利用ABAQUS的用户自定义单元子程序功能编写FORTRAN程序，定义复合单元。在 ABA卵S中使用牛顿法迭代求解瞬态溫度场，求解过程中需要计算牛顿迭代法的雅克比矩阵 AMATRX和残余向量列阵畑S。复合单元中，AMATRX和畑S的计算公式为：
[0化1 ]
[0化2]
[0化3]
[0化4]
[0化5]
[0化6]
[0化7]
[0化引 [0化9]
[0060]其中，TB为计算时间段末的节点溫度，为圆周率，d为水管直径，h为水管的散热系 数，1为水管长度，鸣。为二次六面体单元高斯积分的第i个积分点的权重，N。为混凝±单元的 形函数，4为线性杆单元高斯积分的第i个积分点的权重，化〇为水管入口在混凝±子单元 局部坐标下的形函数矩阵，化1为水管出口在混凝上子单元局部坐标下的形函数矩阵，V为冷 却水流速，Pw为水密度，Cw为水的比热容，&为系数，用于调整求解瞬态溫度场收敛准则中水 管溫度的容许误差。
[0061 ]使用修改后的inp文件创建ABAQUS分析作业，将复合单元子程序路径添加到"user subroutine fiIe"输入框，提交分析作业，得到每个诱筑阶段的溫度场。
[0062]本发明能够有效地模拟冷却水管对大体积混凝±结构的整体冷却效果，并且可W 反应水管周围的溫度场梯度变化;划分网格时不需考虑水管布置，水管节点溫度为有限元 支配方程的基本未知量，求解混凝±节点溫度的同时可W得到水管的节点溫度，计算时不 需要添加额外的迭代步来求解水管溫度，具有很高的前处理效率和计算效率。
[00创实施例一
[0064] W某重力巧为例，已知两年中，该巧段从基础约束区向上诱筑至泄洪道，分86次诱 筑，每个诱筑块高Im左右，总诱筑高度为70m，计算该巧段上块、中块两年内的溫度场。计算 施工期溫度场需要考虑诱筑次序，水管冷却等因素。该巧段宽21米，上块、中块诱筑块长度 为25m、38m。巧段的几何模型如图2。混凝±和地基的热学参数通过室内试验获得，混凝±的 绝热溫升18.2°C，其余参数如表1所示。水管为铁管，在诱筑层面上蛇形布置，通水后每两天 变换一次通水方向。初期通水自诱筑后12小时开始。第一年入冬之前11月至12月，对当年秋 天诱筑的导流底孔部分进行中后期通水冷却。混凝±诱筑块的初始溫度为混凝±的入仓溫 度，冬季入仓溫度较低，在1 (TC左右，夏季时入仓溫度较高，为14 °C左右。
[0(?日]表1重力巧材料参数 [0066]
[0067]采用本发明的方法进行计算溫度场的具体步骤如下：
[006引第一步，划分混凝±网格;使用S维建模软件AutoCAD建立巧体和地基的S维实体 模型，然后按照诱筑块分割巧体，再划分混凝上单元和地基单元，不同的诱筑块相邻的单元 共用节点;模型共有1，446,801个节点，340,902个单元。
[0069] 第二步，划分水管网格。由第一根水管开始，根据水管的布置方式计算水管各个拐 点的坐标，建立水管的=维模型。按照水管与混凝上单元轮廓的相交关系，将水管划分为水 管辅助单元。
[0070] 第S步，创建复合单元。将水管辅助单元与包含水管辅助单元对应的混凝±辅助 单元组合构成新的复合单元。复合单元的单元类型为用户自定义单元UEL。复合单元共22个 节点，29个积分点，在ABA卵S的i吨文件中添加如下两行语句来定义肥L的属性：
[0071 ] 水 user element ,variables = 27, coordinates = 3 ,unsymm, nodes = 22, type = ul, properties = 3
[0072] 11
[0073] 第一行中，*user element为关键词，表示该段代码用于定义用户自定义单元。其 他变量为该类型的用户自定义单元的参数，variable = 27表示单元需要在计算过程中存储 27个临时变量，coordinates = 3表示单元用于S维模型的计算，unsymm表示单元的热传导 矩阵和比热容矩阵为非对称矩阵，nodes = 22表示单元有22个节点，type = ul表示该用户自 定义单元的名称为"ur，propedies = 3表示单元需要定义=个材料参数。第二行中11表示 该单元用于计算瞬态溫度场；
[0074] 第四步，定义材料属性。按照表1给混凝±单元、地基单元、混凝±辅助单元定义材 料的导热系数、比热容、密度。并编写ABA卵S子程序肥TVAL来定义混凝±的绝热溫升。
[0075] 第五步，定义水管编号。将每一根水管的水管辅助单元创建单元集合，在inp文件 中通过初始场变量命令*initial conditions定义水管编号。
[0076] 第六步，定义冷却水信息。文件中每一行为某一诱筑阶段的冷却水信息，共9个值， 分别为:水管编号、该条信息的起始时间、是否通水、通水方向是否逆向、水管内径、流速、水 管厚度、水管导热系数、水管入口水溫。
[0077] 第屯步，计算地基溫度场。使用ABAQUS中"interaction"模块中的生死单元 ("model change")功能杀死除地基单元外的其他单元。计算诱筑前两年内地基溫度场，地 表溫度根据拟合的气溫曲线定义。
[0078] 第八步，编写用户自定义单元子程序肥L，计算诱筑期溫度场。将地基溫度场作为 初始溫度场，按照诱筑顺序逐层激活诱筑层所包含的混凝上辅助单元、水管辅助单元、混凝 ±单元和复合单元，计算诱筑期溫度场。
[0079] 附图3至5给出了该方法的计算结果云图。其中，图3为第二层诱筑层诱筑两天后的 溫度场云图，图4为当时的水管沿程水溫，最高溫度位于第一层混凝±下表面，达到了20.71 °C，第二层混凝±下表面由于有冷却水管的作用，溫度得到了明显的控制。同时，前两天冷 却水入口位于诱筑块左侧，左侧水管附近的溫度(8.57度）明显低于右侧水管附近的溫度 (约为15至17度）。图5给出了第一年冬天(12月27日）的溫度场云图，图中可W看出冷却水管 周围溫度较低，为8 °C左右，离水管越远，冷却效果逐渐减弱。
[0080] 本发明的模拟方法所得的结果符合混凝±溫度场的变化规律，能够有效地模拟冷 却水管对大体积混凝±结构的整体冷却效果，并且可W反应水管周围的溫度场梯度变化； 划分网格时不需考虑水管布置，计算时不需要迭代求解水管溫度，具有很高的前处理效率 和计算效率。
[0081] W上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明技术原理的前提下，还可W做出若干改进和变型，运些改进和变型 也应视为本发明的保护范围。
【主权项】
1. 一种大体积混凝土水管冷却温度场的模拟方法，其特征是，包括以下步骤： 步骤一，划分混凝土结构与地基网格:将混凝土与地基实体模型划分成混凝土单元和 地基单元，单元类型为二次六面体单元，每个单元包含20个节点和27个积分点，并对节点进 行编号； 步骤二，划分水管网格:水管被混凝土单元轮廓切分为若干段，每一段记为一个水管辅 助单元，每个单元包含2节点和2个积分点，并对节点进行编号，包含水管辅助单元的混凝土 单元记为混凝土辅助单元； 步骤三，创建复合单元:将每对水管辅助单元和混凝土辅助单元创建新的复合单元，并 对复合单元的节点进行编号； 步骤四，定义材料属性:定义混凝土单元、地基单元的材料属性，水管辅助单元材料属 性为0，混凝土辅助单元的材料属性为混凝土属性； 步骤五，创建水管编号:将属于同一根水管的水管辅助单元创建单元集合，集合名称为 水管编号； 步骤六，创建冷却水信息文件:将所有浇筑阶段的冷却水信息按浇筑阶段划分来创建 文件； 步骤七，计算浇筑时地基的温度场:将气温作为地基表面的第三类边界条件，将当地年 平均气温作为地基的初始温度，计算得到浇筑时地基的温度场； 步骤八，创建自定义单元子程序计算各阶段的温度场:将地基温度场作为初始温度场， 按照浇筑顺序逐层激活混凝土辅助单元、水管辅助单元、混凝土单元和复合单元，计算得到 各浇筑阶段的温度场。2. 根据权利要求1所述的一种大体积混凝土水管冷却温度场的模拟方法，其特征是，在 所述步骤三中，复合单元为用户自定义的单元，包含22个节点和29个积分点，其中1~20节 点编号对应混凝土辅助单元的节点，21及22节点编号对应水管辅助单元的节点。3. 根据权利要求1所述的一种大体积混凝土水管冷却温度场的模拟方法，其特征是，在 所述步骤四中，材料属性包括材料的导热系数、比热容和密度。4. 根据权利要求1所述的一种大体积混凝土水管冷却温度场的模拟方法，其特征是，在 所述步骤六中，冷却水信息包括水管编号、该条信息的起始时间、是否通水、通水方向是否 逆向、水管内径、流速、水管厚度、水管导热系数、水管入口水温。5. 根据权利要求1所述的一种大体积混凝土水管冷却温度场的模拟方法，其特征是，在 所述步骤八中，自定义单元子程序中复合单元的雅克比矩阵AMATRX和残余向量RHS的计算其中，Γ为计算时间段末的节点温度，调用自定义单元子程序时作为参数传入。31为圆周 率，d为水管直径，h为水管的散热系数，1为水管长度，为二次六面体单元高斯积分的第i 个积分点的权重，N。为混凝土子单元形函数，4为线性杆单元高斯积分的第i个积分点的权 重，N2Q为水管入口在混凝土子单元局部坐标下的形函数矩阵，N 21为水管出口在混凝土子单 元局部坐标下的形函数矩阵，V为冷却水流速，Pw为水密度，(^为水的比热容，Φ为系数，用于 调整求解瞬态温度场收敛准则中水管温度的容许误差。
【文档编号】G06F17/50GK106021720SQ201610334937
【公开日】2016年10月12日
【申请日】2016年5月19日
【发明人】董义佳, 苏超, 徐超, 任志明
【申请人】河海大学