仙居抽水蓄能电站蜗壳混凝土施工仿真研究

2016-12-14 05:32何少云许其光王泽华蒋林华
关键词:仙居热传导蜗壳

何少云,许其光,王泽华,张 研,蒋林华

(1.浙江仙居抽水蓄能有限公司,浙江 仙居 317300;2.中国水利水电第十二工程局有限公司施工科学研究所,浙江 建德 311600;3.河海大学 力学与材料学院,江苏 南京 210098)

仙居抽水蓄能电站蜗壳混凝土施工仿真研究

何少云1,许其光2,王泽华3,张 研3,蒋林华3

(1.浙江仙居抽水蓄能有限公司,浙江 仙居 317300;2.中国水利水电第十二工程局有限公司施工科学研究所,浙江 建德 311600;3.河海大学 力学与材料学院,江苏 南京 210098)

浙江仙居抽水蓄能电站装配了大功率水轮发电机组,其蜗壳的大体积混凝土施工存在温度应力大并造成开裂的风险,影响机组正常运行。为探究蜗壳温度应力,对1#机组蜗壳混凝土进行了考虑施工过程的三维实时仿真计算,对关键位置进行了应力分析,提出施工期降低早期温度应力的建议措施。

仙居抽水蓄能电站;蜗壳;施工过程;三维仿真

1 研究背景

在大体积混凝土结构浇筑过程中,对温度以及温度应力的控制具有重要意义。温度变化对结构的应力具有很大的影响[1]。大体积混凝土浇筑完成后,在固结硬化过程中将会放出大量的水化热,进而导致混凝土自身温度的上升。由于混凝土的导热性能比较差,结构内部早期水化热产生较大的温度梯度,往往造成结构开裂,甚至出现贯穿性裂缝,影响工程安全[2]。

浙江仙居抽水蓄能电站位于浙江省仙居县湫山乡境内,枢纽工程主要由上水库、输水系统、地下厂房、地面开关站及下水库等建筑物组成。从2013年9月开始主体工程混凝土浇筑,至2014年6月已浇筑混凝土方量约8.5万m3。电站1#—4#机组将按次序先后浇筑,其蜗壳结构混凝土浇筑量大,1#机组首先浇筑且浇筑时间跨越夏、秋、冬三季,其早期温度变化幅度较大,引起的早期温度应力很可能会造成严重开裂事故。考虑到机组蜗壳的重要性,若其发生开裂影响整个电站的正常运营。本文应用基于等效时间理论的混凝土非线性热传导方程,借助GID有限元分析软件对仙居电站1#机组蜗壳大体积混凝土的浇筑过程开展三维实时仿真研究,结合理论和有限元手段,计算分析典型位置的应力特征,提出相应的温控措施。

2 理论基础

2.1 经典的热传导方程

2.1.1 热传导方程 将混凝土考虑成均匀、各向同性的固体,其内部或表面的某一点,在某一瞬时的温度Ti=f(x,y,z,t),考虑水泥的水化热作用,在混凝土内部释放放出水化热[3-4],则依据Fourier热传导理论,可得到下列三维不稳定热传导方程:

式中,D为导温系数,m2/h,D=λ/cρ,λ为导热系数(kJ/m·h·℃),c是比热容(k J/kg·℃),ρ为混凝土密度;θ为混凝土绝热温升(℃)。

2.1.2 定解条件 通过热传导方程,可建立物体的温度与空间、时间的关系,但满足热传导方程的解有无限多个,为了确定需要的温度场,还必须知道初始条件和边界条件。初始条件和边界条件合称为定解条件[5-6]。

初始条件:当t=0时,初始瞬时的温度分布为:

边界条件:确定需要的温度场所需的边界条件共有3种,分别为已知温度边界、绝热边界和热交换边界[7]。

(1)第一类边界条件(已知温度边界),已知温度边界上混凝土表面温度T是时间t的已知函数

式中,Tb为已知气温或水温,为给定的温度。

(2)第二类边界条件(绝热边界),若混凝土边界为绝热边界,与周围介质无热量交换,则有:

式中,n为边界的外法线方向。

(3)III第三类边界条件(热交换边界),第三类边界条件假定经过混凝土表面的热流量与混凝土表面温度T和气温Ta之差成正比,即:

式中,β为混凝土表面放热系数,k J/(m2·h·℃);Ta为外界气温。

当表面放热系数β趋于无穷大时,T=Ta,热交换边界条件转化为已知温度边界条件;当表面放热系数β=0,=0,热交换边界条件转化为绝热边界条件。

2.2 基于等效时间理论的热传导方程

2.2.1 等效时间的定义 水泥和水的水化反应是一种放热反应,一般情况下,只要有水和水泥,水化反应的速率就会随着温度的升高而加快,并且温度对化学反应速率的影响服从Arrhenius方程[8],可表示为:

式中,T为绝对温度;k为化学反应速率;R为气体常数(R=8.3144J/k-mole);E为与化学活动能有关的常数。

在1970年,Bazant教授根据Arrhenius方程提出了成熟函数用来计算相对于参考温度Tr的等效时间te:

式(7)可用来定量计算养护时间和温度对混凝土的影响,该式的离散形式为:

式中,Q定义为化学活动能与气体常数之商(Q=E/R);T为在时间间隔Δt内混凝土的平均温度。

应用气体常数R时,Tr和T需要采用绝对温度。美国ASTM规范建议,在缺乏试验资料的情况下,计算温度和时间对1型混凝土(Type 1)强度的影响时,可以采用Q=5000 K。

在研究温度对混凝土强度的影响时,Tank和Carino采用了如下表达式计算等效时间:

式中,T为养护温度;Tr为参考温度;Bt为温度敏感系数(℃-1)。

式(9)提供了一个表示等效时间te的一个更简便的形式。

2.2.2 等效绝热温升 研究表明,浇筑温度在混凝土水化热产生初期会对混凝土水化反应速率有较大的影响,但温度对最终绝热温升影响不大,主要是影响化学反应速率。不管浇筑温度如何,最高绝热温升由混凝土的给定而确定。等效时间是综合考虑混凝土的浇筑温度和时间对绝热温升的影响而得到表示在某一参考温度情况下产生的绝热温升所需要的时间,在其它不同温度下也能达到。对于给定的某种混凝土,最高绝热温升θu和等效时间之间存在着唯一关系,与浇筑温度无关,因此,

混凝土的绝热温升可以用等效时间的函数来唯一表示:

式中,θeq(te)为基于等效时间的绝热温升,称为等效绝热温升,℃;θu为混凝土最高绝热温升,℃;te为相对于参考温度的等效时间,d;m为某一个常数,d,表示达到最高绝热温升二分之一时所需的等效时间。

2.2.3 基于等效时间的热传导方程 如果假定混凝土在浇筑过程中满足能量守恒定律并且混凝土绝热温升可以用Arrhenius理论描述,则求解混凝土三维不稳定温度场的热传导方程为:

式中,x、y、z为直角坐标;T(x、y、z、t)为温度场;D为混凝土导温系数,D=λ/cρ,c为混凝土比热,ρ为混凝土质量密度。

2.3 材料热、力学参数计算公式 计算应力应用的是混凝土的弹性和线性徐变理论,即在同一时刻作用的应力所引起的各相同时刻的徐变值与应力的大小是成正比的。在任意时刻t,龄期为τ的混凝土作用的单位应力(σ=1)时,其总应变为:

式中,φ、α、b、C0、A0、B0、R0、C1、A1、B1、R1可根据试验资料确定。

3 蜗壳施工仿真

3.1 基本气候资料 浙江仙居气候属于亚热带季风区,年平均气温18.3℃,1月份平均气温5.6℃,7月份平均气温28.5℃。由于仙居水电站枢纽工程主厂房为地下厂房,水轮机组处于山洞中,并没有直接处于外部环境中,因此外界各月旬平均温度应取山洞内气温,如表1。

表1 山洞内各月平均气温 (单位:℃)

蜗壳结构地下厂房工程1#—4#机组段尾水管高程为93.8~101.95 m,蜗壳层高程为101.95~110.35 m,水轮机层高程为110.35~116.75m,中间层高程为116.75~122.95m,蜗壳层结构混凝土强度等级为C25F50。

3.2 施工方案 浙江仙居抽水蓄能电站1#发电机组蜗壳采用分期分批浇筑方式,分为两期共9批浇筑。一期混凝土为前8批,二期混凝土为第9批。详细施工进度如表2。

表2 施工进度

3.3 计算模型 计算模型为仙居水电站主厂房一号水轮机混凝土蜗壳及地基的简化,简化过程中忽略了其中的输水廊道和进人廊道的细节部分。根据圣维南原理,蜗壳地基长度尺寸大致取蜗壳基座的2倍,高度尺寸与长度尺寸大致相同。蜗壳底部海拔100.7 m,地基深30 m,长宽各45m。采用分期分批浇注方式,结构为典型的三维结构。前文的数学模型的数值实现可详见文献[1]。为了能够比较准确地反映出结构各部位的受力特点,采用了较密的剖分单元,用八结点等参单元作为结构的有限元单元形式。生成的模型包含22 705个六面体单元和25 056个节点。包括周围土体的蜗壳主体结构的有限元网格剖分如图1所示。

图1 蜗壳模型

3.4 材料热、力学属性参数 天然土层地基假设为一种材料,用I来表示;蜗壳混凝土假设为两种材料,分两期浇筑,分别用Ⅱ、Ⅲ表示。它们的主要力学参数如表3、表4及表5。

3.5 保温措施 蜗壳结构浇筑分为两期浇筑,一期混凝土8批,二期混凝土1批,共分9批,竖向模具为钢模,每批浇筑1 d后拆模,考虑实际施工过程,数值仿真时不考虑冷却水管,有入仓温控措施,表面无保温措施。拆模前后的表面放热系数分别为4.5和6.25。

表3 蜗壳材料常数

表4 弹性模量变化参数

表5 徐变参数

4 计算结果分析

4.1 应力计算结果 依据以上的理论研究和实际工程环境条件,开展了蜗壳混凝土的三维实时仿真计算。可以得到不同浇筑时刻的应力大小,典型时刻的整体应力分布如图2所示。

图2 典型时刻的蜗壳应力分布

对内部节点进一步分析,可得到不同批次混凝土在典型时刻的应力大小,列于表6。其中危险位置发生在第三批、第四批浇筑后4 d新浇筑层外表面;第九批浇筑后新老混凝土交界处外表面。其温度和应力时程曲线如图3和图4。现场对检测得到蜗壳混凝土的轴拉强度为3MPa。在标准养护条件下混凝土3 d强度能达到标准强度的35%、7 d能达到50%左右。在实际情况下混凝土同期所能达到的强度偏低。从表6应力计算结果可以看出,在整个浇筑过程中,蜗壳温度拉应力均未超过其轴拉强度。但是,第三批、第四批混凝土由于尺寸较大在浇筑完成4 d后就分别达到了1.12 MPa和0.99 MPa,超过早期混凝土抗拉强度,有可能产生早期温度裂缝。二期混凝土浇筑4 d后拉应力达到0.91 MPa,此处的拉应力是原来一期混凝土上的拉应力,不是新混凝土所受拉应力,因此并不会对新混凝土造成破坏。由于一期混凝土的约束作用,新混凝的拉应力也偏大,在浇筑第4 d即达到0.65 MPa,需要采取降低拉应力的措施。

表6 应力计算结果

图3 危险位置温度时程曲线

图4 危险位置应力时程曲线

4.2 温控建议与措施 通过对以上应力计算结果的分析可知,在当前的条件下,仙居电站1#发电机组蜗壳的浇筑过程中,第三、四、九批混凝土有可能出现早期温度裂缝,需要采取措施以降低早期温度应力,建议措施如下:(1)优化混凝土配合比。适当减少水泥用量,降低混凝土水化热,提高混凝土极限拉伸和降低混凝土热膨胀系数;(2)降低浇筑温度。入仓温度控制在23℃,在拌和前对骨料和石子进行预冷处理,冷却拌和水或掺加冰屑;(3)加强混凝土外表面保护。混凝土浇筑完毕后,保持混凝土表面湿润状态,禁止洒温度过低的水,防止内外温差引起裂缝;(4)混凝土浇筑应尽量避开高温时间。夏季浇筑混凝土时,尽量安排在早晨和夜间进行。

5 结语

通过对浙江仙居抽水蓄能电站1#发电机组蜗壳的施工仿真研究,可以得出,在浇筑过程中,第三、四、九批混凝土有可能出现早期温度裂缝,需要采取一定的温控措施降低早期温度应力。其余批次混凝土在当前施工条件下基本满足设计要求,能够保证施工期混凝土不开裂。

[1] 张研,韩林,蒋林华,等.大体积混凝土温度应力与裂缝控制[M].北京:科学出版社,2014.

[2] 张研,等.基于细观尺度的混凝土热膨胀性能研究[J].建筑材料学报,2011,14(3):310-316.

[3] 吕宏基,杨德福,等.大体积混凝土[M].北京:水利电力出版社,1990.

[4] 张子明,傅作新.龙滩碾压混凝土重力坝的仿真计算[J].红水河,1997,16(1):13-17,28.

[5] 龚召熊.水工混凝土的温控与防裂[M].北京:中国水利电力出版社,1999.

[6] Zhang ZM,Garga V K.State of temperature and thermal stress in mass concrete structures subjected to thermal shock[J].Dam Engineering,1996,7(4):21-27.

[7] 潘家铮.水工建筑物的温度控制[M].北京:水利电力出版社,1990.

[8] 张子明,宋智通,黄海燕.混凝土绝热温升和热传导方程的新理论[J].河海大学学报,2002,30(3):1-6.

Numerical analysis of the volute of the turbine generator sets at Xianju pumped storage power station

HE Shaoyun1,XU Qiguang2,WANG Zehua3,ZHANG Yan3,JIANG Linhua3
(1.Zhejiang Xianju Pumped-Storage Power Co.Ltd,Xianju 317300,China;2.Sinohydro Bureau 12 Co.,Ltd.Research on the construction of Science,Jiande 311600,China;3.College of Mechanics and Materials,Hohai University,Nanjing 210098,China)

Zhejiang Xianju pumped storage power station will be equiped with high-power turbine generator sets.The massive concrete of the volutes would cause high temperature stress and has the risk of thermal cracking.This would affect the normal operation of the whole station.To explore the temperature stress of the spiral case,for unit 1 volute concrete considering the construction process of three-dimensional real-time simulation calculation,the stress analysis was carried out on the key position,and recommended measures were put forward to reduce early thermal stress during construction.

Xianju pumped storage power station;volute;construction process;three-dimensional simulation

TU528

A

10.13244/j.cnki.jiwhr.2016.01.005

1672-3031(2016)01-0029-07

(责任编辑:祁 伟)

2015-09-13

国家自然科学基金项目(50808066,512781670);国家教育部基金项目(2009512812)

何少云(1975-),男,江西崇仁人,高级工程师,主要从事水利水电工程研究。E-mail:316381350@qq.com

许其光(1963-),男,浙江建德人,高级工程师,主要从事工程试验与研究。E-mail:xpg8889@126.com

猜你喜欢
仙居热传导蜗壳
仙居杨梅卖迪拜 每颗售价60元
焊接蜗壳泵模型开发
一类三维逆时热传导问题的数值求解
冬天摸金属为什么比摸木头感觉凉?
污水泵蜗壳断裂原因分析
杜光庭《題仙居觀》所咏地理位置考
热传导对5083 铝合金热压缩试验变形行为影响的有限元模拟研究
杨余龙
仙人居住的地方 浙江·仙居
高功率低燃油耗的新一代蜗壳式增压器