CN103926271A - 一种反演混凝土绝热温升的方法 - Google Patents

Abstract

一种获取混凝土绝热温升的方法，涉及一种使用温度传感器和温度采集模块测量采集温度，利用实测温度计算反演绝热温升的方法。本发明通过在混凝土同一剖面沿径向布设不同位置的温度传感器，测量不同时间不同测点混凝土的温度，通过差分计算温度场，利用实测温度与计算温度优化反演绝热温升。本发明并不着手于减少混凝土的散热，而是通过测量散热过程中的混凝土的温度反演绝热温升，经过数据分析后证明可以较好地反映出混凝土后期的放热过程。

Description

一种反演混凝土绝热温升的方法

技术领域

本发明涉及一种获取混凝土绝热温升的方法，特别涉及获取混凝土放热后期的绝热温升。

背景技术

绝热温升关系到大坝施工中的温度控制，是大体积混凝土的一项重要性能。目前一般只进行28d龄期的绝热温升试验，仪器的绝热室要达到绝热试验条件，即胶凝材料的水化热不散发到外界环境。试件表面是裸露的，绝热室温度要跟踪混凝土试样中心温度，相差不大于±0.1℃。随着龄期的增长，二者温度差值导致的热量累积损失随之变大，因而难以保证试验精度。

在该试验方法基础上进行改进，形成了一种新的试验方法，在混凝土试件外面包裹保温层，其他步骤依然按之前的试验方法进行，由于有了保温层，降低了等效散热系数，即使热量累计损失稍大一些，对混凝土温升的影响也小得多。保温层导热系数越小，厚度越厚，试验精度越高。该方法虽然一定程度上减小了误差，但无法完全阻止混凝土散热，在龄期较长时也不能够很好地反映混凝土后期放热过程，试验精度依旧不高。

发明内容

本发明的目的是提供一种反演混凝土绝热温升的方法，通过测量散热过程中的混凝土的温度反演绝热温升，使其可以更精确地获取绝热温升，并能较好地反映出混凝土后期的放热过程，更符合工程实际。

本发明的目的是通过如下技术方案实现的：

一种反演混凝土绝热温升的方法，其特征在于该方法包括如下步骤：

1）首先制作圆筒式空心模具，在空心模具下部铺设一层保温材料[2]，用来模拟径向一维散热，同时在模具四周分别铺设保温材料，厚度为下部保温材料的三分之一到五分之一，用来减缓散热速率；

2）浇筑混凝土至空心模具中间高度时，模具内沿径向布置并固定多个温度传感器[3]；

3）继续浇筑混凝土，振捣混凝土，振捣完成后，在模具上部铺设与下部相同厚度的保温材料，立即放入养护室中养护；

4）使用温度采集模块每隔固定时间Δτ采集混凝土内温度传感器的温度信号；

5）利用差分方法计算温度场，将模具内的混凝土分成多个薄圆环层，圆环圆心均为模具圆心，每层厚度设为h，试取出相邻三层中同一径向且相邻两点间距为h的i-1，i，i+1三点来分析，三点到圆心距离依次增大，由热传导方程及径向一维散热得到下式：

$T_{i,\mathrm{\τ}+\mathrm{\Δ\τ}}=(1-2b)T_{i,\mathrm{\τ}}+(b-\frac{\mathrm{bh}}{2r})T_{i-1,\mathrm{\τ}}+(b+\frac{\mathrm{bh}}{2r})T_{i+1,\mathrm{\τ}}+\mathrm{\Δ\θ}---\left(1\right)$

b=aΔτ/h²    （2）

式中，T_i,τ为第i点在时间τ的温度，T_i,τ+Δτ为第i点在时间τ+Δτ的温度，r为第i点到圆心距离，h为单层厚度，Δτ为时间步长，a为导温系数，Δθ为从τ到τ+Δτ时间段内的绝热温升值，b为中间变量；

6）由（1）式得到，第i点，时间τ+Δτ的温度可以由时间τ，i-1，i，i+1三点温度来表示，其中包含两个参数，Δθ和b，b≤1/2时（1）式计算结果稳定；

7）设温度传感器测温点数量为n个，n个点均在同一径向，点1～n到圆心距离依次增大，1，n两点温度为边界条件，根据（1）式，由初始时刻1～n点的温度以及时间Δτ，1，n两点的温度，得到时间Δτ，2～n-1点的计算温度T₂～T_n-1，设时间Δτ，2～n-1点的实测温度为t₂～t_n-1，进行优化计算求出Δθ和b使得（3）式最小：

$\underset{i=2}{\overset{n-1}{\mathrm{\Σ}}}{(T_i-t_i)}^2---\left(3\right)$

重复步骤7），由时间Δτ，1～n各点温度以及时间2Δτ，1，n两点的温度算得时间2Δτ，2～n-1点的计算温度，与实测温度比较，进行优化计算求出此时段的Δθ和b，以此类推，算出各个时段的温升。

本发明与现有技术相比，具有以下优点及突出性效果：①由于没有传统方法中的累积损失热量，本方法可以更好地反映出混凝土后期的放热过程。②本方法使用的公式简洁明了，便于编制相关程序进行计算。③传统方法在放热后期混凝土温度偏高，不符合工程实际。而本方法可以通过控制环境温度的方法，使得试件温度与工程中混凝土的温度变化相符合，这样所得到的绝热温升值更符合工程实际。

附图说明

图1为圆筒式空心模具纵剖面图。

图2为圆筒式空心模具横剖面图。

图3为圆筒式空心模具内点温度计算示意图。

其中：1-混凝土；2-保温材料；3-温度传感器，其中图2中从左到右分别是1～7号温度传感器。

具体实施方式

下面结合附图进一步说明本发明的实施方式：

本发明通过在混凝土同一剖面沿径向布设不同位置的温度传感器，测量不同时间不同测点混凝土的温度，再通过差分方法计算温度场，利用实测温度与计算温度差的平方和最小优化反演绝热温升。具体方式可按如下步骤实施：

1）首先制作圆筒式空心模具，两侧可拆卸便于搬运及拆模，中间高度处两侧留孔，在空心模具下部铺设一层保温材料[2]，用来模拟径向一维散热，同时在模具四周分别铺设保温材料，厚度为下部保温材料的三分之一到五分之一，用来减缓散热速率，如图1所示；

2）浇筑混凝土至空心模具中间高度时，将准备好的长方形塑料片通过模具两侧小孔放入模具中，长方形塑料片上按预设距离布设小孔，小孔尺寸与传感器测温头相同，将温度传感器[3]***长方形片预留的孔中以保证传感器位置沿径向按照一定间隔固定，传感器布设方式如图2所示；

3）继续浇筑混凝土，振捣混凝土，振捣完成后，在模具上部铺设与下部相同厚度的保温材料，立即放入养护室中养护；

4）使用温度采集模块每隔固定时间Δτ采集混凝土内温度传感器的温度信号；

5）差分方法计算温度场，柱坐标形式的热传导方程为：

$\frac{\∂T}{\∂\mathrm{\τ}}=a(\frac{{\∂}^{2}T}{{\∂r}^2}+\frac{1}{r}\frac{\∂T}{\∂r}+\frac{1}{r^2}\frac{{\∂}^{2}T}{{\∂\mathrm{\φ}}^2}+\frac{{\∂}^{2}T}{{\∂z}^2})+\frac{\∂\mathrm{\θ}}{\∂\mathrm{\τ}}$

式中，T为温度，a为导温系数，θ为绝热温升，τ为时间，r为径向，z为轴向，φ为角度方向。

由于本试验模拟径向一维散热，所以可以假定z、φ方向不传热，上式可以简化为：

$\frac{\∂T}{\∂\mathrm{\τ}}=a(\frac{{\∂}^{2}T}{{\∂r}^2}+\frac{1}{r}\frac{\∂T}{\∂r})+\frac{\∂\mathrm{\θ}}{\∂\mathrm{\τ}}$

将模具内的混凝土分成多个薄圆环层，圆环圆心均为模具圆心，每层厚度设为h，试取出相邻三层中同一径向且相邻两点间距为h的i-1，i，i+1三点来分析，三点到圆心距离依次增大，如图3所示。

根据差分原理，忽略截断误差，温度的各阶段偏导数可表示如下：

${\left(\frac{\∂T}{\∂r}\right)}_{i,\mathrm{\τ}}=\frac{1}{2h}(T_{i+1,\mathrm{\τ}}-T_{i-1,\mathrm{\τ}})$

${\left(\frac{\∂T}{\∂r}\right)}_{i+h/2,\mathrm{\τ}}=\frac{1}{h}(T_{i+1,\mathrm{\τ}}-T_{i,\mathrm{\τ}})$

${\left(\frac{\∂T}{\∂r}\right)}_{i-h/2,\mathrm{\τ}}=\frac{1}{h}(T_{i,\mathrm{\τ}}-T_{i-1,\mathrm{\τ}})$

$\begin{array}{c}{\left(\frac{{\∂}^{2}T}{{\∂r}^2}\right)}_{i,\mathrm{\τ}}=\frac{1}{h}[{\left(\frac{\∂T}{\∂r}\right)}_{i+h/2,\mathrm{\τ}}-{\left(\frac{\∂T}{\∂r}\right)}_{i-h/2,\mathrm{\τ}}]\\ =\frac{1}{h^2}(T_{i-1,\mathrm{\τ}}+T_{i+1,\mathrm{\τ}}-2T_{i,\mathrm{\τ}})\end{array}$

用向前差分计算及

${\left(\frac{\∂T}{\∂r}\right)}_{i,\mathrm{\τ}}=\frac{1}{\mathrm{\Δ\τ}}(T_{i+1,\mathrm{\τ}}-T_{i,\mathrm{\τ}})$

${\left(\frac{\∂\mathrm{\θ}}{\∂\mathrm{\τ}}\right)}_{i,\mathrm{\τ}}=\frac{\mathrm{\θ}(\mathrm{\τ}+\mathrm{\Δ\τ})-\mathrm{\θ}\left(\mathrm{\τ}\right)}{\mathrm{\Δ\τ}}=\frac{\mathrm{\Δ\θ}}{\mathrm{\Δ\τ}}$

把以上各式代到柱坐标一维热传导方程中，得到：

$T_{i,\mathrm{\τ}+\mathrm{\Δ\τ}}=(1-2b)T_{i,\mathrm{\τ}}+(b-\frac{\mathrm{bh}}{2r})T_{i-1,\mathrm{\τ}}+(b+\frac{\mathrm{bh}}{2r})T_{i+1,\mathrm{\τ}}+\mathrm{\Δ\θ}$

b=aΔτ/h²

式中，T_i,τ为第i点在时间τ的温度，T_i,τ+Δτ为第i点在时间τ+Δτ的温度，r为第i点到圆心距离，h为单层厚度，Δτ为时间步长，a为导温系数，Δθ为从τ到τ+Δτ时间段内的绝热温升值，b为中间变量；

6）由上式得到，第i点，时间τ+Δτ的温度可以由时间τ，i-1，i，i+1三点温度来表示，其中包含两个参数，Δθ和b，b≤1/2时上式计算结果稳定；

7）设传感器测温点数量为n个，n个点均在同一径向，点1～n到圆心距离依次增大，1，n两点温度为边界条件，根据上式，由初始时刻1～n点的温度以及时间Δτ，1，n两点的温度，得到时间Δτ，2～n-1点的计算温度T₂～T_n-1，设时间Δτ，2～n-1点的实测温度为t₂～t_n-1，进行优化计算求出Δθ和b使得下式最小：

$\underset{i=2}{\overset{n-1}{\mathrm{\Σ}}}{(T_i-t_i)}^2$

重复步骤7），由时间Δτ，1～n各点温度以及时间2Δτ，1，n两点的温度算得时间2Δτ，2～n-1点的计算温度，与实测温度比较，进行优化计算求出此时段的Δθ和b，以此类推，算出各个时段的温升。

实施例：

在具体的实施中，圆柱体内直径400mm，高600mm。保温材料采用橡塑海绵，上下部铺设80mm的保温材料，四周分别铺设20mm的保温材料。在中间高度剖面上共布设七个温度传感器，其中一个在圆心，沿径向布置，传感器间距均为50mm，如图2所示。每5分钟读取一次温度，取时间步长Δτ为5min，每层厚度h为50mm。

将七个测量位置分为两组，每组测点均在同一径向，分别为图2所示的1234号传感器测点和7654号传感器测点。取每组的最初的测量值为初始条件，第一组1、4号传感器测量值为边界条件，第二组7、4号传感器为边界条件。

先由初始值算到5分钟后的温度值，并与实测温度值作比较，进行优化计算求出导温系数与此时段的绝热温升，使得各个计算位置，即2、3、5、6号传感器位置计算温度与实测温度的差的平方和最小。再由5分钟后的温度值算到10分钟后的温度值，与实测温度值比较，进行优化计算求出此时段的导温系数与绝热温升。以此类推，算出各个时段的绝热温升，画出绝热温升曲线。

随后进行同配合比混凝土绝热室温度跟踪混凝土试样中心温度的传统方法进行绝热温升试验，画出绝热温升曲线，并与之前的计算绝热温升曲线进行比较。比较得出，在初始放热阶段两者比较接近，随着龄期的增长，前者趋于平缓，而后者仍保持着一定速率上升，两者的差距逐渐拉大，因而可以证明本方法可以体现出传统方法所无法体现出的混凝土后期放热过程。

Claims (1)

1.一种反演混凝土绝热温升的方法，其特征在于该方法包括如下步骤： 

1)首先制作圆筒式空心模具，在空心模具下部铺设一层保温材料[2]，用来模拟径向一维散热，同时在模具四周分别铺设保温材料，厚度为下部保温材料的三分之一到五分之一，用来减缓散热速率； 

2)浇筑混凝土至空心模具中间高度时，模具内沿径向布置并固定多个温度传感器[3]； 

3)继续浇筑混凝土，振捣混凝土，振捣完成后，在模具上部铺设与下部相同厚度的保温材料，立即放入养护室中养护； 

4)使用温度采集模块每隔固定时间Δτ采集混凝土内温度传感器的温度信号； 

5)利用差分方法计算温度场，将模具内的混凝土分成多个薄圆环层，圆环圆心均为模具圆心，每层厚度设为h，试取出相邻三层中同一径向且相邻两点间距为h的i-1，i，i+1三点来分析，三点到圆心距离依次增大，由热传导方程及径向一维散热得到下式： 

b＝aΔτ/h²    (2) 

式中，T_i,τ为第i点在时间τ的温度，T_i,τ+Δτ为第i点在时间τ+Δτ的温度，r为第i点到圆心距离，h为单层厚度，Δτ为时间步长，a为导温系数，Δθ为从τ到τ+Δτ时间段内的绝热温升值，b为中间变量； 

6)由(1)式得到，第i点，时间τ+Δτ的温度可以由时间τ，i-1，i，i+1三点温度来表示，其中包含两个参数，Δθ和b，b≤1/2时(1)式计算结果稳定； 

7)设温度传感器测温点数量为n个，n个点均在同一径向，点1～n到圆心距离依次增大，1，n两点温度为边界条件，根据(1)式，由初始时刻1～n点的温度以及时间Δτ，1，n两点的温度，得到时间Δτ，2～n-1点的计算温度T₂～T_n-1，设时间Δτ，2～n-1点的实测温度为t₂～t_n-1，进行优化计算求出Δθ和b使得(3)式最小： 

重复步骤7)，由时间Δτ，1～n各点温度以及时间2Δτ，1，n两点的温度算得时间2Δτ，2～n-1点的计算温度，与实测温度比较，进行优化计算求出此时段的Δθ和b，以此类推，算出各个时段的温升。