CN107085370B - 混凝土仓面小气候自适应控制方法 - Google Patents
混凝土仓面小气候自适应控制方法 Download PDFInfo
- Publication number
- CN107085370B CN107085370B CN201710259197.4A CN201710259197A CN107085370B CN 107085370 B CN107085370 B CN 107085370B CN 201710259197 A CN201710259197 A CN 201710259197A CN 107085370 B CN107085370 B CN 107085370B
- Authority
- CN
- China
- Prior art keywords
- concrete
- temperature
- hours
- coefficient
- old
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000005855 radiation Effects 0.000 claims abstract description 15
- 239000007921 spray Substances 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 35
- 239000010410 layer Substances 0.000 claims description 30
- 230000036571 hydration Effects 0.000 claims description 17
- 238000006703 hydration reaction Methods 0.000 claims description 17
- 238000005266 casting Methods 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 239000011229 interlayer Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims 6
- 238000005507 spraying Methods 0.000 abstract description 41
- 230000007246 mechanism Effects 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- On-Site Construction Work That Accompanies The Preparation And Application Of Concrete (AREA)
Abstract
本发明公开一种混凝土仓面小气候自适应控制方法,利用数据采集子系统采集的混凝土仓面的温度、风速、太阳辐射等参数及喷雾机的水温参数等,根据设定的目标浇筑温度,确定混凝土仓面所需的目标喷雾量,根据计算得到的目标喷雾量自适应地调整喷雾机的功率,以使喷雾机的喷雾量达到目标喷雾量,进而保证混凝土仓面可达到目标浇筑温度;在现有的混凝土仓面小气候控制系统基础上,能够根据实时采集的混凝土仓面的各项参数,自适应地调节喷雾量,实现混凝土仓面温度的准确控制调节。
Description
技术领域
本发明涉及一种混凝土仓面小气候自适应控制方法,属于水利水电工程技术领域。
背景技术
控制浇筑温度是混凝土温控防裂的重要措施之一,控制混凝土仓面环境温度是控制浇筑温度的重要手段。
中国专利CN201520908963.1公开了一种混凝土仓面小气候控制系统,其包括喷雾机构、数据采集子系统、仓面气候控制子系统,数据采集子系统采集的混凝土仓面的气候数据(温湿度数据、风速数据、太阳辐射数据等)传输至仓面气候控制子系统,仓面气候控制子系统根据该气候数据及预设的气候阈值或气候阈值范围控制喷雾机构的喷雾方式、喷雾范围及喷雾水压。该系统能够实现根据混凝土仓面的各参数实时控制喷雾机构动作,实现混凝土仓面的温度控制。其中的喷雾机构形式多种多样,喷雾机即是常用的喷雾装置,通过调节不同的功率可实现不同的喷雾量,达到不同的降温效果。现有的喷雾机已经具有自动调节功率和喷雾量的功能,但是,在混凝土仓面温度控制过程中,尚没有根据采集的混凝土仓面的温度等参数,自适应调节喷雾机功率、喷雾量的方法。
发明内容
鉴于上述原因,本发明的目的在于提供一种混凝土仓面小气候自适应控制方法,可根据混凝土入仓温度、仓面温度、风速、太阳辐射等参数及喷雾机水温参数自适应地调节喷雾量,使得混凝土仓面的浇筑温度达到目标浇筑温度,实现混凝土仓面温度的准确控制调节。
为实现上述目的,本发明采用以下技术方案:
一种混凝土仓面小气候自适应控制方法,由数据采集子系统采集混凝土仓面的温度、风速、太阳辐射、喷雾机水温参数,包括:
确定喷雾机的喷雾量为:
其中,Tw为喷雾机水温,vw为混凝土仓面的风速,kf和ke为系数,通过试验或现场反分析确定;Tay为外界环境温度,其计算公式为:
Tay=Tat+ΔTa (2)
其中,Tat外界气温;ΔTa为太阳辐射热引起的日平均环境温度增量,计算公式为;
其中,S为混凝土仓面的太阳辐射;
Tas为达到设计的目标浇筑温度所对应的目标混凝土仓面环境温度,其计算公式为:
其中,T1为入仓温度,φ1为平仓作用影响系数,φ2为铺筑层间歇影响系数,为水化放热引起的温度回升;Tpg为去除老铺筑层混凝土热传导影响因素后的浇筑温度,其计算方法为:
Tpg=Tps-ΔTP (5)
其中,Tps为目标浇筑温度,ΔTP为老铺筑层混凝土热传导引起的浇筑温度修正项。
所述平仓作用影响系数φ1的计算方法为,
φ1=kt (6)
其中,t为混凝土入仓后到平仓前所经历的时间,单位为分钟;k为经验系数;小型工程在缺乏资料情况下,取k=0.003(1℃/min);对于大中型坝体,需要通过工程类比或实测方法获得经验系数;
所述铺筑层间歇影响系数φ2的计算方法为:
其中:为导热系数和表面放热系数的比值,取值范围为超过上下限范围时按照对应的上下限值考虑;c为比热,取值范围为0.6~1.2,单位为kJ/(kg·℃),超过上下限范围时按照对应的上下限值取值。
老铺筑层混凝土热传导引起的浇筑温度的修正项ΔTP的计算方法为:
其中,为老混凝土热传导影响系数,TPo为老铺筑层混凝土的实测浇筑温度。
对于水化放热引起的温度回升缺少试验数据但具有绝热温升拟合公式时,依据拟合公式确定水化放热引起的温度回升量,包括:
采用指数形式拟合时,
采用双曲线形式拟合时,
其中,θ0为绝热温升终值;a和b、n均为常数,根据绝热温升曲线性质决定;
等效水化放热时间为:
其中,Tc为绝热温升试验块的初始温度;Δτ为铺筑薄层间歇时间。
当铺筑层厚度为0.3m时,新老混凝土温差为1℃时,老混凝土热传导影响系数取值为:
比热/0.9 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.667 | 0.117 | 0.161 | 0.200 | 0.234 | 0.263 | 0.291 |
0.889 | 0.082 | 0.119 | 0.153 | 0.183 | 0.211 | 0.235 |
1.000 | 0.070 | 0.103 | 0.135 | 0.164 | 0.190 | 0.214 |
1.111 | 0.059 | 0.091 | 0.120 | 0.147 | 0.173 | 0.196 |
1.333 | 0.044 | 0.070 | 0.096 | 0.121 | 0.144 | 0.166 |
其中,比热的单位为kJ/(kg·℃)。
当铺筑层厚度为0.4m时,新老混凝土温差为1℃时,老混凝土热传导影响系数取值为:
比热/0.9 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.667 | 0.036 | 0.060 | 0.084 | 0.108 | 0.132 | 0.155 |
0.889 | 0.020 | 0.036 | 0.054 | 0.072 | 0.091 | 0.109 |
1.000 | 0.015 | 0.028 | 0.044 | 0.060 | 0.076 | 0.093 |
1.111 | 0.012 | 0.022 | 0.036 | 0.050 | 0.065 | 0.080 |
1.333 | 0.007 | 0.015 | 0.025 | 0.036 | 0.048 | 0.060 |
当铺筑层厚度为0.5m时,新老混凝土温差为1℃时,老混凝土热传导影响系数取值为:
比热/0.9 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.667 | 0.011 | 0.021 | 0.035 | 0.052 | 0.070 | 0.088 |
0.889 | 0.005 | 0.010 | 0.018 | 0.028 | 0.039 | 0.052 |
1.000 | 0.003 | 0.007 | 0.013 | 0.020 | 0.030 | 0.040 |
1.111 | 0.002 | 0.005 | 0.010 | 0.016 | 0.023 | 0.031 |
1.333 | 0.001 | 0.003 | 0.006 | 0.009 | 0.014 | 0.020 |
本发明的优点是:
本发明的混凝土仓面小气候自适应控制方法,利用数据采集子系统采集的混凝土仓面的温度、风速、太阳辐射等参数及喷雾机的水温参数等参数,根据设定的目标浇筑温度,确定混凝土仓面所需的目标喷雾量,根据计算得到的目标喷雾量自适应地调整喷雾机的功率,以使喷雾机的喷雾量达到目标喷雾量,进而保证混凝土仓面可达到目标浇筑温度,在现有的混凝土仓面小气候控制系统基础上,能够根据实时采集的混凝土仓面的各参数,自适应地调节喷雾量,达到准确控制调节混凝土仓面温度的目的。
附图说明
图1是不同导热系数和表面放热系数的比值条件下单位温度差引起的混凝土温升的公式计算值和有限元计算验证值的对比结果示意图。
图2是图1所示对比结果的误差分析图。
图3是不同比热条件下单位温度差引起的混凝土温升的公式计算值和有限元计算验证值的对比结果示意图。
图4是图3所示对比结果的误差分析图。
图5是本发明于一具体实施例中喷雾机的部署方式图。
具体实施方式
以下结合附图和实施例对本发明作进一步详细的说明。
本发明公开的混凝土仓面小气候自适应控制方法,包括:
1、浇筑期间水化热温升。
等效水化放热时间为:
其中,Tc为绝热温升试验块的初始温度;T1为入仓温度;为等效水化放热时间,即绝热温升试验进行的时间;Aτ为铺筑薄层间歇时间。
水化放热引起的温度回升为最好由试验数据直接得到。缺少试验数据但拥有绝热温升拟合公式时候,可依据拟合公式确定水化放热引起的温度回升量。
采用指数形式拟合时,水化放热引起的温度回升为:
其中,θ0为绝热温升终值;a和b均为常数,根据绝热温升曲线性质决定。
采用双曲线形式拟合时,水化放热引起的温度回升为:
其中,θ0为绝热温升终值;n为常数,根据绝热温升曲线性质所决定。
使用式(1)-式(3)时应特别注意铺筑薄层间歇时间的单位,避免出现因时间单位而引起的错误。
2、入仓温度、外界环境温度与浇筑温度的关系。
混凝土浇筑期间,外界环境温度引起的浇筑温度增量可采用迭代方法计算。
计算初值按下式考虑:
其中,Ta为环境温度;为水化放热引起的温度回升;T1为入仓温度;φ1为平仓作用影响系数,φ2为铺筑层间歇影响系数,计算公式分别为:
φ1=kt (5)
其中,t为混凝土入仓后到平仓前所经历的时间,单位为分钟;k为经验系数。小型工程缺乏资料情况下,取k=0.003(1℃/min);对于大中型坝体,需要通过工程类比或实测方法获得经验系数。
其中,为导热系数和表面放热系数的比值,取值范围为超过上下限范围时按照对应的上下限值考虑;c为比热,取值范围为0.6~1.2,单位为kJ/(kg·℃),超过上下限范围时按照对应的上下限值考虑。
浇筑温度通过迭代计算得到,迭代递推计算公式为:
当n≥3时,即进行三次以上迭代可获得精确的浇筑温度计算值。
混凝土浇筑过程中,下层混凝土的热传导将影响新浇筑混凝土的浇筑温度,新老混凝土温差为1℃时,老混凝土热传导影响系数取值参见表1~表3。铺筑层厚度小于0.3m时,应专题研究可行性;铺筑层厚度在0.3m~0.5m之间时,热传导修正项取值按插值取值;铺筑层厚度大于0.5m时,按照铺筑层厚度为0.5m对应的热传导修正项取值确定。
最终,得到混凝土浇筑温度为:
Tp=Tp估n+ΔTP (8)
其中,Tp为混凝土浇筑温度,单位为℃;ΔTP为老混凝土热传导引起的浇筑温度的修正项,单位为℃;Tp估n为第n次的估算值,一般取3次估算值。
由老混凝土热传导引起的浇筑温度修正项ΔTP的计算方法为:
表1当铺筑层厚度为0.3m时,老混凝土热传导影响系数值
比热/0.9 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.667 | 0.117 | 0.161 | 0.200 | 0.234 | 0.263 | 0.291 |
0.889 | 0.082 | 0.119 | 0.153 | 0.183 | 0.211 | 0.235 |
1.000 | 0.070 | 0.103 | 0.135 | 0.164 | 0.190 | 0.214 |
1.111 | 0.059 | 0.091 | 0.120 | 0.147 | 0.173 | 0.196 |
1.333 | 0.044 | 0.070 | 0.096 | 0.121 | 0.144 | 0.166 |
表2当铺筑层厚度为0.4m时,老混凝土热传导影响系数值
表3当铺筑层厚度为0.5m时,老混凝土热传导影响系数值
比热/0.9 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.667 | 0.011 | 0.021 | 0.035 | 0.052 | 0.070 | 0.088 |
0.889 | 0.005 | 0.010 | 0.018 | 0.028 | 0.039 | 0.052 |
1.000 | 0.003 | 0.007 | 0.013 | 0.020 | 0.030 | 0.040 |
1.111 | 0.002 | 0.005 | 0.010 | 0.016 | 0.023 | 0.031 |
1.333 | 0.001 | 0.003 | 0.006 | 0.009 | 0.014 | 0.020 |
3、对上述预测混凝土浇筑温度的方法进行验证。
式(4)中,铺筑层间歇影响系数φ2的取值和验证。
设:
φ2=φ21φ22Δτ (10)
其中,φ21为表面热交换影响系数,φ22为内部热传导的影响系数。
φ21为表面热交换影响系数,其与导热系数和表面放热系数的比值有关,故设:
其中,a和b为待定系数。
φ22为内部热传导影响系数,其与导温系数相关,考虑到混凝土密度相差不大,故设:
其中,d和e为待定系数。
1)表面热交换影响系数系数φ21确定
根据朱伯芳院士的《大体积混凝土温度应力与温度控制》的理论解,可采用以下计算条件确定φ21和φ22:
浇筑温度按0℃考虑,外界温度按10℃考虑。混凝土浇筑模型顶面散热,其余面绝热。混凝土浇筑模型高度为0.5m,导热系数为164kJ/m·d·℃,比热为0.9kJ/kg·℃,计算表面放热系数为100-1200kJ/m2·d·℃情况下,单位温度差引起的混凝土温升。将系数d的取值定为182.2,则φ22=1,根据计算结果即得到比热0.9kJ/kg·℃且导热系数取值为164kJ/m·d·℃时,表面放热系数变化情况下单位温度差引起的单位时间混凝土温升,并由此得到系数φ21。表4为在浇筑间歇期Δτ分别为3~8个小时情况下,φ21φ22值(φ22=1)。根据计算结果,取值在0.164到1.64期间,φ21的取值和浇筑间歇的时间无关且和取值密切相关。
表4浇筑间歇期Δτ分别为3~8个小时情况下,φ21φ22值(φ22=1)
λ/β | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.164 | 0.053 | 0.053 | 0.052 | 0.050 | 0.049 | 0.047 |
0.205 | 0.045 | 0.045 | 0.044 | 0.043 | 0.042 | 0.041 |
0.273 | 0.036 | 0.036 | 0.035 | 0.035 | 0.034 | 0.033 |
0.410 | 0.025 | 0.025 | 0.025 | 0.025 | 0.024 | 0.024 |
0.547 | 0.020 | 0.020 | 0.020 | 0.019 | 0.019 | 0.019 |
0.820 | 0.013 | 0.014 | 0.014 | 0.014 | 0.013 | 0.013 |
1.093 | 0.010 | 0.010 | 0.010 | 0.010 | 0.010 | 0.010 |
1.640 | 0.007 | 0.007 | 0.007 | 0.007 | 0.007 | 0.007 |
当间歇时间为5小时、比热为0.9kJ/kg·℃、导热系数取值164kJ/m·d·℃、φ22=1,且铺筑层间歇影响系数取值为时,有限元和公式拟合计算结果及相应的误差分析结果见图1和图2。公式拟合值和有限元计算值十分吻合。
2)内部热传导的影响系数φ22确定
计算条件:浇筑温度按0℃考虑,外界温度按10℃考虑。混凝土浇筑模型顶面散热,其余面绝热;导热系数为164kJ/m·d·℃,表面放热系数为600kJ/m2·d·℃,计算比热为0.6-1.2kJ/kg·℃情况下单位温度差引起的混凝土温升,并由此得到φ22。
表5和表6分别为导热系数为164kJ/m·d·℃,表面放热系数为600kJ/m2·d·℃时,不同比热情况下φ21φ22和φ22的取值。根据计算结果可知,φ22值和浇筑间歇的时间无关,仅和比热相关。
当间歇时间为5小时、导热系数为164kJ/m·d·℃、表面放热系数为600kJ/m2·d·℃且时,有限元和公式拟合计算结果及相应的误差分析结果见图3和图4。
表5表面放热系数为600kJ/m2·d·℃时,不同比热情况下的φ21φ22值
比热 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.8 | 0.040 | 0.040 | 0.039 | 0.038 | 0.037 | 0.036 |
0.9 | 0.036 | 0.036 | 0.035 | 0.035 | 0.034 | 0.033 |
1.0 | 0.032 | 0.032 | 0.032 | 0.032 | 0.031 | 0.030 |
1.2 | 0.026 | 0.027 | 0.027 | 0.027 | 0.027 | 0.026 |
表6表面放热系数为600kJ/m2·d·℃时,不同比热情况下φ22值
比热 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.8 | 1.12 | 1.12 | 1.12 | 1.09 | 1.09 | 1.08 |
0.9 | 1.00 | 1.00 | 1.01 | 0.99 | 0.99 | 0.99 |
1.0 | 0.88 | 0.90 | 0.93 | 0.91 | 0.91 | 0.92 |
1.2 | 0.72 | 0.75 | 0.77 | 0.77 | 0.79 | 0.80 |
3)φ2适用范围的验证
根据式(10)式(11)和式(12),可以确定φ2的取值:
本发明采取了不同的混凝土材料参数验证式(13)的正确性,表7和表8分别为表面放热系数为900kJ/m2·d·℃、导热系数为164kJ/m·d·℃,不同的比热情况下计算得到的φ21φ22值和φ22值;表9和表10分别为表面放热系数为300kJ/m2·d·℃、导热系数为164kJ/m·d·℃,不同的比热情况下,计算得到的φ21φ22值和φ22值。根据计算结果,φ21φ22值和φ22值均与铺筑层浇筑间歇时间无关,且均能较好的利用式(13)拟合。
表7导热系数为164kJ/m·d·℃,不同比热情况下φ21φ22值
比热 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.8 | 0.055 | 0.055 | 0.053 | 0.051 | 0.049 | 0.048 |
0.9 | 0.049 | 0.049 | 0.048 | 0.047 | 0.045 | 0.044 |
1.0 | 0.044 | 0.045 | 0.044 | 0.043 | 0.042 | 0.041 |
1.2 | 0.036 | 0.037 | 0.037 | 0.037 | 0.036 | 0.035 |
表8导热系数为164kJ/m·d·℃,不同比热情况下φ22值
比热 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
0.8 | 1.13 | 1.11 | 1.11 | 1.09 | 1.10 | 1.08 |
0.9 | 1.01 | 1.01 | 1.01 | 1.00 | 1.01 | 1.00 |
1.0 | 0.90 | 0.91 | 0.92 | 0.92 | 0.93 | 0.92 |
1.2 | 0.74 | 0.76 | 0.78 | 0.79 | 0.80 | 0.80 |
表9导热系数为164kJ/m·d·℃,不同比热情况下φ21φ22值
比热 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
800 | 0.022 | 0.022 | 0.022 | 0.021 | 0.021 | 0.020 |
900 | 0.020 | 0.020 | 0.020 | 0.019 | 0.019 | 0.019 |
1000 | 0.017 | 0.018 | 0.018 | 0.018 | 0.017 | 0.017 |
1200 | 0.014 | 0.015 | 0.015 | 0.015 | 0.015 | 0.015 |
表10导热系数为164kJ/m·d·℃,不同比热情况下φ22值
比热 | 3小时 | 4小时 | 5小时 | 6小时 | 7小时 | 8小时 |
800 | 1.10 | 1.11 | 1.10 | 1.13 | 1.10 | 1.07 |
900 | 0.98 | 0.99 | 0.99 | 1.02 | 1.00 | 0.98 |
1000 | 0.87 | 0.89 | 0.90 | 0.93 | 0.92 | 0.90 |
1200 | 0.70 | 0.73 | 0.75 | 0.79 | 0.78 | 0.77 |
4、由老混凝土热传导引起的浇筑温度的修正项ΔTP的取值方法。
计算条件:浇筑温度按0℃考虑,外界温度按10℃考虑。混凝土浇筑模型底面散热,其余面绝热;导热系数为164kJ/m·d·℃,表面放热系数为100000kJ/m2·d·℃,计算比热为0.6-1.2kJ/kg·℃情况下,单位温度差引起的混凝土温升。计算结果见表1~表3,分别对应铺筑层厚度分别为0.3m、0.4m、0.5m情况下,老混凝土热传导引起的老混凝土热传导影响系数值。
5、确定喷雾量
上述式(7)~式(9)是根据已知的入仓温度、混凝土仓面环境温度得到预测的混凝土浇筑温度。那么,在设计院提供的目标浇筑温度(已知定值)的基础上,根据式(7)~式(9)进行反推,可得到为达到该目标浇筑温度所对应的目标混凝土仓面环境温度,进一步的,根据所要达到的目标混凝土仓面温度,确定喷雾机所需提供的喷雾量。具体为:
根据目标浇筑温度,得到达到该目标浇筑温度所需的目标混凝土仓面环境温度为:
其中,T1为入仓温度,可由人工实时监测获得,φ1为平仓作用影响系数,φ2为铺筑层间歇影响系数,为水化放热引起的温度回升;Tpg为去除老铺筑层热传导影响因素后的浇筑温度,其计算方法为:
Tpg=Tps-ΔTP (15)
其中,Tps为给定的目标浇筑温度,ΔTP为老混凝土热传导引起的浇筑温度修正项,其根据式(19)确定;即在目标混凝土浇筑温度基础上,去除老混凝土热传导影响的混凝土浇筑温度增量。
根据式(14)计算得到目标混凝土仓面环境温度,确定喷雾机的喷雾量为:
其中,Tw为水温,其由喷雾机中装设的温度传感器采集的温度信号确定;vw为风速,其由混凝土仓面上方安装的风速传感器采集的风速信号确定;kf和ke为系数,通过试验或现场反分析确定;Tay为外界环境温度,其根据式(17)确定:
Tay=Tat+ΔTa (17)
其中,Tat为气温,其由混凝土仓面上方安装的温度传感器采集的温度信号确定;ΔTa为太阳辐射热引起的日平均环境温度增量,单位为℃,计算公式为;
S为太阳辐射的实测值,单位为kJ/(m2·h),其由混凝土仓面上方安装的太阳辐射传感器采集的太阳辐射信号确定。
本发明中,老铺筑层是指混凝土已浇筑结束的铺筑层,新铺筑层是指在已浇筑结束的老铺筑层的表面正在浇筑的混凝土层。由老铺筑层混凝土热传导引起的浇筑温度的修正项ΔTP的计算方法为:
其中,为老混凝土热传导影响系数,TPo为老铺筑层混凝土的实测浇筑温度。
即,本发明是利用现有的数据采集子系统采集混凝土仓面的各项参数(包括温度、风速、太阳辐射等)及喷雾机中的水温参数等参数,根据设定的目标浇筑温度,利用式(16)确定混凝土仓面所需的目标喷雾量,根据计算得到的目标喷雾量调整喷雾机的功率,使得喷雾机的喷雾量达到目标喷雾量,进而使得混凝土仓面达到目标浇筑温度。
于一具体实施例中,乌东德电站坝体沿河流左右岸的宽度一般为20m,大坝一般修建在河谷中,往往有大风的作用。如图5所示,在河流两岸部署数据采集子系统及喷雾机,喷雾机的间距为20m,喷雾机距离上游和下游坝体表面10m左右;为使喷雾量均匀,喷雾机的水平转角设置在180°以上,垂直转角设置为向下10°以上、向上60°;根据设计的目标浇筑温度及计算得到的目标喷雾量,调整喷雾机的功率,使得喷雾机的喷雾量达到目标喷雾量,保证混凝土仓面的浇筑温度达到目标浇筑温度。
以上所述是本发明的较佳实施例及其所运用的技术原理,对于本领域的技术人员来说,在不背离本发明的精神和范围的情况下,任何基于本发明技术方案基础上的等效变换、简单替换等显而易见的改变,均属于本发明保护范围之内。
Claims (7)
1.混凝土仓面小气候自适应控制方法,由数据采集子系统采集混凝土仓面的温度、风速、太阳辐射、喷雾机水温参数,其特征在于,
确定喷雾机的喷雾量为:
其中,Tw为喷雾机水温,vw为混凝土仓面的风速,kf和ke为系数,通过试验或现场反分析确定;Tay为外界环境温度,其计算公式为:
Tay=Tat+ΔTa (2)
其中,Tat外界气温;ΔTa为太阳辐射热引起的日平均环境温度增量,计算公式为;
其中,S为混凝土仓面的太阳辐射;
Tas为达到设计的目标浇筑温度所对应的目标混凝土仓面环境温度,其计算公式为:
其中,T1为入仓温度,φ1为平仓作用影响系数,φ2为铺筑层间歇影响系数,为水化放热引起的温度回升;Tpg为去除老铺筑层混凝土热传导影响因素后的浇筑温度,其计算方法为:
Tpg=Tps-ΔTP (5)
其中,Tps为目标浇筑温度,ΔTP为老铺筑层混凝土热传导引起的浇筑温度修正项。
2.根据权利要求1所述的混凝土仓面小气候自适应控制方法,其特征在于,所述平仓作用影响系数φ1的计算方法为,
φ1=kt (6)
其中,t为混凝土入仓后到平仓前所经历的时间,单位为分钟;k为经验系数;小型工程在缺乏资料情况下,取k=0.003℃/min;对于大中型坝体,需要通过工程类比或实测方法获得经验系数;
所述铺筑层间歇影响系数φ2的计算方法为:
其中:为导热系数和表面放热系数的比值,取值范围为超过上下限范围时按照对应的上下限值考虑;c为比热,取值范围为0.6~1.2,单位为kJ/kg·℃,超过上下限范围时按照对应的上下限值取值;Δτ为铺筑薄层间歇时间。
3.根据权利要求1所述的混凝土仓面小气候自适应控制方法,其特征在于,老铺筑层混凝土热传导引起的浇筑温度的修正项ΔTP的计算方法为:
其中,为老混凝土热传导影响系数,为老铺筑层混凝土的实测浇筑温度。
4.根据权利要求1所述的混凝土仓面小气候自适应控制方法,其特征在于,对于水化放热引起的温度回升缺少试验数据但具有绝热温升拟合公式时,依据拟合公式确定水化放热引起的温度回升量,包括:
采用指数形式拟合时,
采用双曲线形式拟合时,
其中,θ0为绝热温升终值;a和b、n均为常数,根据绝热温升曲线性质决定;
等效水化放热时间为:
其中,Tc为绝热温升试验块的初始温度;Δτ为铺筑薄层间歇时间。
5.根据权利要求3所述的混凝土仓面小气候自适应控制方法,其特征在于,当铺筑层厚度为0.3m时,新老混凝土温差为1℃时,老混凝土热传导影响系数取值为:
其中,比热的单位为kJ/kg·℃。
6.根据权利要求3所述的混凝土仓面小气候自适应控制方法,其特征在于,当铺筑层厚度为0.4m时,新老混凝土温差为1℃时,老混凝土热传导影响系数取值为:
其中,比热的单位为kJ/kg·℃。
7.根据权利要求3所述的混凝土仓面小气候自适应控制方法,其特征在于,当铺筑层厚度为0.5m时,新老混凝土温差为1℃时,老混凝土热传导影响系数取值为:
其中,比热的单位为kJ/kg·℃。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710259197.4A CN107085370B (zh) | 2017-04-19 | 2017-04-19 | 混凝土仓面小气候自适应控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710259197.4A CN107085370B (zh) | 2017-04-19 | 2017-04-19 | 混凝土仓面小气候自适应控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107085370A CN107085370A (zh) | 2017-08-22 |
CN107085370B true CN107085370B (zh) | 2019-07-23 |
Family
ID=59611666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710259197.4A Expired - Fee Related CN107085370B (zh) | 2017-04-19 | 2017-04-19 | 混凝土仓面小气候自适应控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107085370B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109164852B (zh) * | 2018-07-27 | 2020-11-27 | 同济大学 | 一种大体积混凝土自适应温度与应力控制方法 |
CN109092587B (zh) * | 2018-08-29 | 2020-05-19 | 中国水利水电科学研究院 | 控制用于混凝土仓面降温的喷雾机的方法 |
CN113886926B (zh) * | 2021-10-14 | 2023-02-17 | 中国水利水电科学研究院 | 高压水道引水岔管温控防裂方法 |
CN113898346B (zh) * | 2021-10-14 | 2023-11-03 | 中国水利水电科学研究院 | 盾构工作井底板温控防裂方法 |
CN114277799A (zh) * | 2021-10-14 | 2022-04-05 | 中国水利水电科学研究院 | 抽水蓄能电站地下厂房蜗壳外围混凝土结构温控防裂方法 |
CN114319248B (zh) * | 2022-01-05 | 2024-01-23 | 中国水利水电第一工程局有限公司 | 一种混凝土坝仓面养护方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011065238A (ja) * | 2009-09-15 | 2011-03-31 | Ecogold Co Ltd | 構造物内温度計算方法及び当該計算方法を実施するプログラム |
CN103605888A (zh) * | 2013-11-14 | 2014-02-26 | 中国水利水电科学研究院 | 一种基于实时监测数据的大体积混凝土温度过程预测方法 |
CN105138048B (zh) * | 2015-08-05 | 2017-06-30 | 中国水利水电科学研究院 | 混凝土浇筑仓面小环境温湿度智能控制装置及方法 |
CN205091640U (zh) * | 2015-11-13 | 2016-03-16 | 中国水利水电科学研究院 | 混凝土仓面小气候控制系统 |
-
2017
- 2017-04-19 CN CN201710259197.4A patent/CN107085370B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107085370A (zh) | 2017-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107085370B (zh) | 混凝土仓面小气候自适应控制方法 | |
CN106844989B (zh) | 最高环境温度浇筑时段混凝土浇筑温度的快速计算方法 | |
Alavi et al. | Prediction of asphalt pavement temperature profile with finite control volume method | |
CN105353620B (zh) | 一种基于地表水深信息的地面灌溉控制方法 | |
CN105842278B (zh) | 一种室内测量混凝土表面放热系数的方法 | |
CN105326386B (zh) | 水温沸点判断方法和装置以及电热水壶 | |
CN104112065B (zh) | 热水器和热水器剩余洗浴时间的计算方法 | |
JP6646391B2 (ja) | コンクリートの養生制御システム | |
US10429827B2 (en) | Method of dynamically calculating refractive index to determine the thickness of roofing materials | |
CN110377981B (zh) | 一种数字化的施工现场混凝土预测式防裂方法 | |
CN109871618B (zh) | 一种在土壤水热耦合模拟中确定底层温度边界的方法 | |
CN116109136A (zh) | 大体积混凝土结构开裂风险评估与预警方法 | |
CN105178605A (zh) | 大体积混凝土智能通水系统 | |
CN104360984A (zh) | 一种基于两点边值测量下冷却塔填料热力性能计算方法 | |
CN111267276A (zh) | 固化参数优化方法、装置、薄膜固化系统及计算机设备 | |
CN209144034U (zh) | 一种混凝土冬季养护装置 | |
CN104090604A (zh) | 热处理设备的温度补偿方法、温度控制方法及系统 | |
CN205197768U (zh) | 水温沸点判断装置和电热水壶 | |
CN105133470B (zh) | 摊铺机的熨平板加热系统、加热方法及摊铺机 | |
CN202854601U (zh) | 通水智能温度控制试验系统 | |
CN204027954U (zh) | 径流泥沙含量实时测量装置 | |
Falk et al. | Investigation of a simple approach to predict rainscreen wall ventilation rates for hygrothermal simulation purposes | |
US20120261402A1 (en) | Apparatus for Remotely Measuring Surface Temperature Using Embedded Components | |
CA2891456C (en) | Dynamically calculated refractive index for determining the thickness of roofing materials | |
WO2019086899A2 (en) | Method and system for determining flow rate in a pipe |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190723 |