大体积混凝土温度场及温度应力的有限元分析

天津城市建设学院学报第18卷第4期2012年12月JournalofTianjinInstituteofUrbanConstructionVol.18No.4Dec.2012收稿日期：2012-06-15；修订日期：2012-07-03作者简介：杨杰（1987—），女，河北邯郸人，天津城市建设学院硕士生.土木工程$大体积混凝土温度场及温度应力的有限元分析杨杰，毛毳，侯霞，王汪阳(天津城市建设学院土木工程学院，天津300384)摘要：结合工程实例，采用ANSYS有限元软件对大体积混凝土浇筑及冷却过程中的温度场及温度应力进行了模拟，通过APDL语言编写程序控制了模拟过程.分析结果与工程计算及实际测试结果相近，验证了所建模型与选取参数的正确性，保证了运用ANSYS研究大体积混凝土水化热影响参数的可行性．分析发现：混凝土温度峰值与最大拉应力与浇筑温度呈正相关；采用低热量水泥时，可降低温度峰值、推迟峰值出现时间，并降低结构温度应力；环境温度变化时，温差成为关键因素，应根据实际情况制定季节性施工方案．关键词：大体积混凝土；有限元方法；温度场；温度应力中图分类号：TU37文献标志码：A文章编号：1006-6853(2012)04-0270-05大体积混凝土浇筑后，在内部水化热和外部环境温度的影响下，其内部温度与外界气温相差很大，极易使混凝土产生温度裂缝，影响结构的安全与正常使用．因此，有必要对大体积混凝土浇筑和养护过程中的温度场及温度应力进行分析，为控制温度裂缝提供依据．笔者采用有限元方法，运用ANSYS热分析模块，结合华中科技大学协和医院外科医疗大楼智能化超高层建筑主楼筏板浇筑实例[1]，建立有限元模型，确定模型材料参数，对主楼筏板浇筑过程进行了模拟，并将温度场计算结果作为荷载施加于模型，计算了混凝土温度应力．在此基础上，对不同水泥标号、不同浇筑温度及不同环境因素对大体积混凝土温度场及应力场的影响进行了分析．1大体积混凝土计算模型的建立1.1模型及材料参数以文献[1]中的工程实例为基础建立有限元模型，尺寸为80,m×36.1,m×2.5,m，混凝土下层为地基，尺寸为120,m×80,m×5,m．主楼筏板大体积混凝土选用亚东P.O42.5级水泥，28,d抗压强度为51.0,MPa.设计配合比如表1所示．表1主楼筏板大体积混凝土配合比3kgm−⋅材料用量水180水泥365砂688石1080II级粉煤灰80UEA29外加剂9.8根据文献[1]提供的混凝土配合比以及施工期间的具体情况，选取模型材料参数如表2所示．表2主要材料参数选取选取参数结构部位导热系数/()1kJmhK−⋅⋅比热/()1kJkgK−⋅密度/3kgm−⋅混凝土8.3880.972431.8基础9.270.836819701.2单元选择及有限元模型的建立模型采用SOLID70单元用于三维瞬态传热分析，该单元具有“单元生死”功能，可根据工序实现浇筑混凝土的动态模拟[2]．此外，该单元可转化为SOLID65结构单元，方便后续的应力分析[3]．依据结构和荷载的对称性，采用1/4模型进行模拟．上层混凝土与下层基础连接处采用共节点耦合温天津城市建设学院学报杨杰等：大体积混凝土温度场及温度应力的有限元分析·271·度自由度来模拟触面上温度协同变化．所建有限元模型如图1所示．图1有限元计算模型1.3初始条件的施加主楼筏板混凝土浇筑时间在11月中旬．选取12,℃作为地基温度初始条件施加于其节点上；混凝土入模温度控制在30,℃左右，以30,℃作为其温度初始条件施加于节点．1.4边界条件的施加(1)在大体积混凝土浇筑过程模拟中，与空气接触的表面采用第三类边界条件[4]，假设通过混凝土表面的热流量与混凝土表面温度T和大气温度Tα的差值成正比，即()TTTnαλβ∂−=−∂(1)式中：β为固体表面放热系数，()12kJmhK−⋅⋅，可根据风速、混凝土表面粗糙程度取值[5]．(2)地基四周与底面假设为绝热边界条件．(3)混凝土对称面采用绝热边界条件．1.5水化热的施加在ANSYS中，水泥水化热是通过生热率HGEN来施加.生热率即单位时间内混凝土的生热量，所以需要对水化热函数求导[4]()()01ebatQtQ−=−(2)()10HGENebbatQtQabt−−′==(3)式中：0Q为混凝土最终水化热，由混凝土配合比计算得出；a、b是和水泥品种相关的水化热常数，查表可知，a＝0.69，b＝0.56．2模拟结果及分析2.1温度场模拟结果主楼筏板大体积混凝土浇筑温度场，12,d后温度场如图2所示．图2主楼筏板浇筑12d后的温度场混凝土中心沿板厚度方向4个测温点的位置如图3所示.图3沿混凝土竖向测温点布置计算温度与实测温度的对比曲线如图4所示．（a）a测温点计算值与实测值对比（b）b测温点计算值与实测值对比（c）c测温点计算值与实测值对比·272·天津城市建设学院学报2012年第18卷第4期（d）d测温点计算值与实测值对比图4计算温度与实测温度对比曲线由图4可以看出：①最高温度出现在c测温点，靠近中间层，为65.6587,℃，与实测最高温度70,℃很接近，表明采用此方法可以较好地模拟混凝土温度的峰值；②模拟的温度变化趋势与实测相近，即在浇筑初期温度急剧上升，温度达到峰值后在养护期间缓慢下降．结果表明，利用ANSYS模拟温度场是可行的，模拟结果可以用来模拟温度应力场，并为工程实际提供参考．2.2温度应力的模拟考虑到混凝土材料的非线性，使用VonMises屈服准则，并采用多线性随动强化选项确定混凝土的非线性本构关系．将混凝土养护28,d时的温度场计算结果作为荷载施加于计算模型，计算温度应力，结果如图5所示．图5混凝土养护28d时第一主应力分布从计算结果可以看出，混凝土中心为压应力，表面为拉应力，最大拉应力为0.700,MPa，与文献[1]计算结果0.874,MPa相近且小于C40混凝土的抗拉设计值．承台中心表面和承台与地基接触边缘拉应力最大．前者因承台中心水化热温度最高，引起混凝土膨胀，使其表面产生较大的拉应力；后者因地基与混凝土温差较大(地基温度为12,℃)，产生较大的拉应力；这表明温差是影响温度应力的关键因素．在此边缘选取四个距中心不同距离的路径，路径上的节点沿板厚方向的第一主应力大小如图6所示．图6距中心不同距离沿板厚方向第一主应力的变化从图6可以看出，在距中心10,m处拉应力最大，比较危险，是分析的重点．3部分影响因素的分析3.1不同水泥品种的影响水泥品种是影响混凝土水化热的重要因素．选取了四种不同标号的水泥，根据表3[6]提供的水化热常数，在其他条件不变的情况下，计算温度场和温度应力．读取浇筑后每天的最高温度，不同水泥标号的最高温度历时曲线如图7所示．表3不同水泥品种的水化热常数水泥品种10/kJkgQ−⋅ab普通硅酸盐水泥425号3300.690.56普通硅酸盐水泥525号3500.360.74普通硅酸盐大坝水泥525号2700.790.70矿渣硅酸盐大坝水泥425号2850.290.76图7使用不同标号水泥时结构最高温度历时曲线从图7可以看出，水泥标号不同，不仅能影响浇筑及养护时的最高温度，而且会改变温度峰值出现的时间．与普通硅酸盐425号水泥相比，普通硅酸盐大坝水泥525号温度峰值较小，相同养护条件下温度下降也较快，28,d时能达到比较低的温度；普通硅酸盐水泥525号温度峰值最高，但推迟了峰值出现的时间，养护期间温度下降较慢；而矿渣硅酸盐大坝水泥天津城市建设学院学报杨杰等：大体积混凝土温度场及温度应力的有限元分析·273·425号温度峰值最低，反应期间水化热上升速度也较慢，峰值出现时间也比较靠后，但养护28,d后与普通硅酸盐425号水泥达到的温度相同．由前面可知，长边边缘距混凝土中心10,m处出现较大的拉应力，比较危险．选取不同水泥品种时，此处沿板厚方向的由表面到底部的路径，查看其第一主应力，结果如图8所示.图8长边边缘距混凝土中心10,m处拉应力变化由图8可知：使用普通硅酸盐水泥525号时，此处拉应力最大；普通硅酸盐水泥425号和矿渣硅酸盐大坝水泥425号次之，且两者比较接近；普通硅酸盐大坝水泥525号最小．说明温度水化热上升缓慢，峰值较小时，结构拉应力相应也较小．3.2浇筑温度的影响通过冷水搅拌、加冰搅拌、加快浇筑速度等方式，都可以控制混凝土的浇筑温度．在环境温度为22.7,℃，相同养护条件下，模拟不同浇筑温度时混凝土温度场，并计算温度应力．不同浇筑温度时最高温度历时曲线如图9所示，最大拉应力随浇筑温度的变化如图10所示．图9浇筑温度不同时结构最高温度历时曲线由图9和图10可以看出：混凝土的温度与浇筑温度呈正相关关系，峰值升高幅度比初温升高的幅度略大；相同养护条件下，温度下降速度大致相同；混凝土结构的最大拉应力与浇筑温度的关系几乎为线性关系，说明混凝土入模时的温度是影响温度应力的关键因素之一，施工时应特别注意控制混凝土的初温．图10结构最大拉应力随浇筑温度变化曲线3.3环境温度的影响环境温度不同，混凝土结构表面与外界的温差就不同，直接影响混凝土的温度场和温度应力．取浇筑温度为20,℃，其他条件均相同，模拟不同环境温度下混凝土温度场和温度应力，得出不同环境温度下最高温度历时曲线，如图11所示．图11环境温度不同时结构最高温度历时曲线由图11可知：环境温度越高，混凝土温度峰值也越高；环境温度的影响主要体现在养护阶段，浇筑期间的影响不大，养护期间温度下降速度大致相同；随着环境温度的降低，混凝土养护28,d后达到的温度也越低，且环境温度越低时，最终温度降低的幅度也越小．将28,d时的温度场结果作为温度荷载施加于混凝土结构，计算温度应力，最大拉应力随环境温度变化如图12所示．从图12可以看出，随着环境温度的升高，混凝土结构最大拉应力呈现先下降后上升的趋势．说明此时混凝土表面与外界的温差成为影响温度应力的主要因素．温差越大，最大拉应力越大，越有可能出现温度裂缝，且低温比高温时拉应力大，表明严寒时期·274·天津城市建设学院学报2012年第18卷第4期更要注意控制内外温差．但并不是外界温度与浇筑温度(20,℃)相同时最大拉应力最小，而是外界温度比浇筑温度略低时，结构最大拉应力最小．图12结构最大拉应力随环境温度变化曲线4结论首先结合工程实例，利用ANSYS模拟了大体积混凝土浇筑过程，又分别模拟了三种因素影响下的大体积混凝土温度场和温度应力，分析结构最高温度和最大拉应力，得到以下结论：(1)使用不同标号的水泥会改变大体积混凝土温度场的峰值，且峰值出现的时间也会不同．其中矿渣硅酸盐大坝水泥425号在浇筑初期放热较慢，且峰值较小，结构出现的拉应力也较小，在考虑降低水化热和温度应力时应优先选用；(2)混凝土浇筑温度越高，达到的温度峰值越高，且最大拉应力与浇筑温度几乎成线性关系．混凝土温度峰值和温度应力的大小与入模温度有直接关系，应尽量采取措施控制浇筑时的初始温度；(3)环境温度越高，混凝土温度峰值越高，但最大拉应力与环境温度并非正相关．此时，混凝土与外界温差成为关键因素，应根据工程实际情况，制定不同季节的不同施工及养护方案，采取措施降低内外温差，防止出现较大的温度应力和应变．参考文献：［1］王玉岚，蒋沧如，苏俊，等.大体积混凝土的温控检测实例[J].新型建筑材料，2009，6：25-27.［2］张朝晖.ANSYS8.0热分析教程与实例解析[M].北京：中国铁道出版社，2005．［3］王国强.实用工程数值模拟技术及其在ANSYS上的实践[M].西安：西北工业大学出版社，1999．［4］王新刚，高洪生，闻宝联.ANSYS计算大体积混凝土温度场的关键技术[J].中国港湾建设，2009，2(1)：41-44.［5］司政，李守义，陈培培，等.基于ANSYS的大体积混凝土温度场计算程序开发[J].长江科学院院报，2011，28(9)：53-56.［6］朱伯芳.大体积混凝土温度应力与温度控制[M].北京：中国电力出版社，1998．FiniteElementAnalysisofMassConcreteTempera