大体积混凝土温度应力场变化分析

第３０卷第１期　２０１３年３月　土木工程与管理学报　Ｖｏｌ＿３０　Ｎｏ．１　Ｍａｒ．２０１３　Ｊｏｕｒｎａｌ　ｏｆ　Ｃｉｖｉｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　ａｎｄ　Ｍａｎａｇｅｍｅｎｔ　大体积混凝土温度应力场变化分析　谭广柱　，　刘书贤　，　张弛　，　麻凤海　（１．辽宁工程技术大学２．辽宁工程技术大学产业开发处，辽宁土木工程学院，辽宁阜新阜新１２３０００；　１２３０００；　３．本溪钢铁集团建设有限责任公司，辽宁４．大连大学摘建筑工程学院，辽宁本溪１　１７０００；　大连１１６６２２；）　要：大体积混凝土水泥水化放热产生的大量热量，由于边界条件的存在和，温度场和应力场的存在非　常容易破坏混凝土，所以大体积混凝土温度场和应力场是其施工质量的关键因素。通过对大体积混凝土温度　场和应力场的现场测试，基于现场实际的测量方案，建立正确的数值分析计算模型，根据大体积混凝土的现场　力学性能非线性增长的特性，分析养护过程中大体积混凝土的温度应力场的分布、变化过程，通过对比现场测　试和数值计算的结果分析可知：需要实时监测大体积｝昆凝土的温度场和现场气象状况，采取合理的养护方案，　保证水泥水化放热能够及时散出，为有效控制大体积混凝土灾害温度应力场的产生和温度裂缝提供理论支持，　其结果可为类似的大体积混凝土工程提供借鉴参考。　关键词：大体积混凝土；温度场；温度应力；有限元；数值分析　中图分类号：ＴＵ３５２．１　文献标识码：Ａ　文章编号：２０９５－０９８５（２０１３）０１－００２０－０５　Ａｎａｌｙｓｉｓ　ｏｎ　Ｔｅｍｐｅｒａｔｕｒｅ　Ｓｔｒｅｓｓ　Ｆｉｅｌｄ　Ｃｈａｎｇｅ　ｏｆ　Ｍａｓｓ　Ｃｏｎｃｒｅｔｅ　ＴＡＮ　Ｇｕａｎｇ．ｚｈｕ　，ＬＩＵ　Ｓｈｕ—ｘｉａｎ　，ＺＨＡＮＧ　Ｃｈｉ　，ＭＡ　Ｆｅｎｇ—ｈａｉ　（１．Ｉｎｄｕｓｔｒｙ　Ｄｅｖｅｌｏｐｍｅｎｔ　Ｄｉｖｉｓｉｏｎ，Ｌｉａｏｎｉｎｇ　Ｔｅｃｈｎｉｃａｌ　Ｕｎｉｖｅｒｓｉｔｙ，Ｆｕｘｉｎ　１２３０００，Ｃｈｉｎａ；　２．Ｃｏｌｌｅｇｅ　ｏｆ　Ｃｉｖｉｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｌｉａｏｎｉｎｇ　Ｔｅｃｈｎｉｃａｌ　Ｕｎｉｖｅｒｓｉｔｙ，Ｆｕｘｉｎ　１　２３０００，Ｃｈｉｎａ；　３．Ｂｅｎｘｉ　Ｉｒｏｎ　ａｎｄ　Ｓｔｅｅｌ　Ｇｒｏｕｐ　Ｃｏｎｓｔｒｕｃｔｉｏｎ　Ｃｏ　Ｌｔｄ，Ｂｅｎｘｉ　１　１　７０００，Ｃｈｉｎａ；　４．Ｃｏｌｌｅｇｅ　ｏｆ　Ｃｉｖｉｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　ａｎｄ　Ａｒｃｈｉｔｅｃｔｕｒｅ，Ｄａｌｉａｎ　Ｕｎｉｖｅｒｓｉｔｙ，Ｄａｌｉａｎ　１　１　６６２２，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ａｓ　ｔｈｅ　ｃｅｍｅｎｔ　ｈｙｄｒａｔｉｏｎ　ｒｅａｃｔｉｏｎ　ｇｉｖｅｓ　ｏｆｆ　ａ　ｌｏｔ　ｏｆ　ｈｅａｔ　ｉｎ　ｃａｓｔｉｎｇ　ｍａｉｎｔｅｎａｎｃｅ　ｐｒｏｃｅｓｓ　ｏｆ　ｍａｓｓ　ｃｏｎｃｒｅｔｅ，ａｎｄ　ｂｏｕｎｄａｒｙ　ｃｏｎｄｉｔｉｏｎｓ　ｏｆ　ｔｈｅ　ｃｏｎｓｔｒａｉｎｔｓ　ｌｅｄ　ｔｏ　ｔｈｅ　ｅｘｉｓｔｅｎｃｅ　ｏｆ　ｔｅｍｐｅｒａｔｕｒｅ　ｓｔｒｅｓｓ，　ｔｅｍｐｅｒａｔｕｒｅ　ｆｉｅｌｄ　ｉｓ　ｔｈｅ　ｃｏｒｅ　ａｎｄ　ｋｅｙ　ｆａｃｔｏｒ　ｏｆ　ｍａｓｓ　ｃｏｎｃｒｅｔｅ　ｉｎ　ｔｈｅ　ｃｏｎｓｔｒｕｃｔｉｏｎ　ｑｕａｌｉｔｙ　ｃｏｎｔｒｏ１．Ｂａｓｅｄ　ｏｎ　ｔｈｅ　ａｃｔｕａｌ　ｍｅａｓｕｒｅｍｅｎｔ　ｓｙｓｔｅｍ，ｔｈｅ　ｃｏｒｒｅｃｔ　ｎｕｍｅｒｉｃａｌ　ｃａｌｃｕｌａｔｉｏｎ　ｍｏｄｅｌ　ｉｓ　ｅｓｔａｂｌｉｓｈｅｄ．Ａｃｃｏｒｄｉｎｇ　ｔｏ　ｔｈｅ　ｒｅａｌ　ｃｏｎｓｔｒｕｃｔｉｏｎ　ｐｒｏｃｅｓｓ　ｏｆ　ｔｅｍｐｅｒａｔｕｒｅ　ｍｅａｓｕｒｉｎｇ　ａｎｄ　ｔｅｍｐｅｒａｔｕｒｅ　ｃｏｎｔｒｏｌ，ａ　ｍｏｒｅ　ｒｅａｓｏｎａｂｌｅ　ｉｆｎｉｔｅ　ｅｌｅｍｅｎｔ　ａｎａｌｙｓｉｓ　ｍｏｄｅｌ　ｉｓ　ｅｓｔａｂｌｉｓｈｅｄ．Ｔｈｒｏｕｇｈ　ｔｈｅ　ｃｏｎｓｉｄｅｒａｔｉｏｎ　ｏｆ　ｔｈｅ　ａｃｔｕａｌ　ｍｅｃｈａｎｉｃｓ　ｐｅｒｆｏｒｍａｎｃｅ　ｏｆ　ｃｏｎｃｒｅｔｅ　ｎｏｎｌｉｎｅａｒ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｇｒｏｗｔｈ，ｔｈｅ　ｐｒｏｃｅｓｓ　ｏｆ　ｍａｓｓ　ｃｏｎｃｒｅｔｅ　ｃｏｎｓｔｒｕｃｔｉｏｎ　ｉｎ　ｔｈｅ　ｃｈａｎｇｅ　ｏｆ　ｔｅｍｐｅｒａｔｕｒｅ　ｉｓ　ａｎａｌｙｚｅｄ．Ｔｈｅ　ｔｅｍｐｅｒａｔｕｒｅ　ｆｉｅｌｄ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ａｎｄ　ｓｔｒｅｓｓ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｓｉｔｕａｔｉｏｎ，ａｓ　ｗｅｌｌ　ａｓ　ｔｈｅ　ｔｅｍｐｅｒａｔｕｒｅ　ｆｉｅｌｄ　ｏｆ　ｍａｓｓ　ｃｏｎｃｒｅｔｅ　ａｎｄ　ｔｈｅ　ｗｅａｔｈｅｒ　ｃｏｎｄｉｔｉｏｎ　ｎｅｅｄ　ｒｅａｌ—ｔｉｍｅ　ｍｏｎｉｔｏｒｉｎｇ，ａｎｄ　ｔａｋｅ　ｒｅａｓｏｎａｂｌｅ　ｍａｉｎｔｅｎａｎｃｅ　ｐｌａｎ　ｉｎ　ｏｒｄｅｒ　ｔｏ　ｅｎｓｕｒｅ　ｔｈａｔ　ｃｅｍｅｎｔ　ｈｙｄｒａｔｉｏｎ　ｈｅａｔ　ｃａｎ　ｐｒｏｍｐｔ　ｓｃａｔｔｅｒｅｄ　ｏｕｔ．Ｔｈｅ　ｎｕｍｅｒｉｃａｌ　ａｎａｌｙｓｉｓ　ｒｅｓｕｌｔｓ　ｍａｋｅ　ａ　ｇｏｏｄ　ａｇｒｅｅｍｅｎｔ　ｗｉｔｈ　ｔｈｅ　ｍｅａｓｕｒｅｄ　ｄａｔａ，　ｔｈｅ　ｒｅｓｕｌｔｓ　ｃａｎ　ｂｅ　ｔｈｅ　ｒｅｆｅｒｅｎｃｅ　ｆｏｒ　ｏｔｈｅｒ　ｓｉｍｉｌａｒ　ｐｒｏｊｅｃｔｓ．　Ｋｅｙ　ｗｏｒｄｓ：ｍａｓｓ　ｃｏｎｃｒｅｔｅ；ｔｅｍｐｅｒａｔｕｒｅ　ｆｉｅｌｄ；ｔｅｍｐｅｒａｔｕｒｅ　ｓｔｒｅｓｓ；ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ；ｎｕｍｅｒｉｃａｌ　ａｎａｌｙｓｉｓ　热造成的，大体积混凝土的温度裂缝主要是由水泥水化　为了提高大体积混凝土结构的施工质　量　于　，必须控制大体积混凝土施工和养护过程中由　水泥水化放热所产生的温度变化和温度应力。　收稿日期：２０１２－０５－１８修回日期：２０１２－０６－０５　作者简介：谭广柱（１９７４一），男，辽宁海城人，硕士，工程师，研究方向为｝昆凝土结构理论（Ｅｍａｉｌ：ｌｎｔｕｌｉｕｓｈｕｘｉａｎ＠１６３．ｃｏｎ１）　基金项目：国家自然科学基金（５１１７４０３８）；辽宁省教育厅科研项目（Ｌ２０１２１ｌ１）；中国煤炭工业协会科学技术研究指导性计划项目　（ＭＴＫＪ２０１２．３１９）；辽宁工程技术大学市场调研信息立项基金（ＳＣＤＹ２０１２０１０）　第１期　谭广柱等：大体积混凝土温度应力场变化分析　・２１・　由于大体积混凝土结构体系自身的复杂性和所处　环境的多样性，大体积混凝土的开裂问题必须考　虑其养护过程中温度应力的影响　。　Ａ　＝一Ａ　（４）　式中，Ａ为导热系数，ｋＪ／（ｍ・ｈ・ｏＣ）。　１．２数学模型　对于大体积混凝土温度应力场的研究分析，　美国主要是利用有限元时问过程的分析方法，日　本则主要是利用约束系数矩阵法；在法国和英国　ＡＢＡＱＵＳ、ＡＮＳＹＳ等有限元分析软件则应用比较　广泛。美国的威尔逊教授　是最早利用有限元　时间过程分析法来分析混凝土温度场的；日本的　考虑到现场施工中的大体积混凝土工程的实　际边界条件的影响因素较多，需要依据实际工程　的详细情况来确定边界条件。基于四维温度场理　论的大体积混凝土的数学力学模型为：　专家学者　不仅考虑了混凝土徐变应力场的计　ｆ　＋　１（０ａ一０　ｒＯＴ）＝０　（５）　算，而且对温度应力场也进行了深入研究；我国朱　｛　（　）：　）　（初始条件）　伯芳院士　编制了我国第一个用于计算混凝土　温度徐变应力的有限元程序，并将其应用于三门　峡重力坝温度应力分析中；刘光廷　等将断裂力　学的研究成果融入功能强大的仿真程序中，应用　“人工短缝”成功地解决了溪柄碾压混凝土薄拱　坝两岸的温度拉应力问题。　为了更进一步的分析研究在养护过程中大体　积混凝土裂缝控制方法，本文基于大体积混凝土　四维温度场与温度应力场理论，通过现场实测，利　用有限元模拟分析具体工程的大体积混凝土四维　温度场，分析数值模拟结果与现场实测结果的温　度场变化趋势及吻合度，以验证数值模拟结果的　可靠性和准确性，为大体积混凝土工程施工方案　制定提供理论依据。　１　大体积混凝土四维温度场理论　１．１基本概念　“四维温度场”理论是指在三维空间（　，Ｙ，　）　的基础上考虑时间一维ｔ的温度场，此时为不稳　定的三维瞬态温度场：　Ｔ＝Ｔ（　，Ｙ，Ｚ，ｔ）　（１）　等温面是指在某一时刻，温度场内具有相同　的温度值各点所构成的一个曲面；温度梯度是指　在温度场内的某点的最大升温速率的单位矢量，　可表示为：　ＡＴ＝ｎｏ　ａＴ＝　ａ凡Ｏｘ＋　＋　（２）　式中，　为单位矢量，沿等温面的法向的温度增　加方向。在单位时间内的单位面积所通过的热量　为热流密度（Ｗ／ｍ　），Ｒ［Ｊ　ｑ＝÷　。　根据热传导定律，热流密度和温度梯度之间　的函数关系可表示为：　ｑ＝一Ａ　Ｖ　Ｔ＝一ｎｏａ　０—１　（３）　【根据具体情况给出的边界条件　本文需要对某综合楼大体积混凝土基础的温　度场在空间域和时间域中的实际分布进行分析研　究，参考文献［５～１０］‘‘四维温度场”理论中的实　际计算模型的相关知识，建立该大体积混凝土基　础四维温度场的理论计算模型如下：　ｆＩｌ　Ｖ２ｌ　Ｔ＋　１（ａ　ｏ一　ＯＴ）＝０　＝Ｔｏ（Ｐ）　｛　．　”　（６）　ｌ　ｌ＝　（Ｐ，ｔ）　Ｉ　．　１．３　控制大体积混凝土温度裂缝的方法　大体积混凝土的温度裂缝在施工养护过程中　分为升温过程和降温过程：　大体积混凝土处于升温阶段时：　ｏｒｒ　一％＋　０＋＋丽　　丽△　（７）　大体积混凝土处于降温阶段：　０．０５ｏ－　．一　＋　丽　＝。。＋　０＋　＿＿　瓦　差　△　（８）　定义温度应力：　（ｔ）＝　（ｔ）＋　Ｎ（ｔ）　（９）　其中：ｏｒ　（ｔ）为外界约束条件下，大体积混凝　土整体变形所产生的应力，　（ｔ）的产生主要是　由于在外界约束存在的条件下，大体积混凝土结　构内部由水化放热所导致的最高温度　降至周　・２２・　土木工程与管理学报　２０１３年　围环境温度Ｔ　（ｔ）产生的收缩应力所引起的；　（　）为内部约束条件下，大体积混凝土表面质点与　内部质点所产生的应力，应力　（ｔ）是在内部约　束条件下，由大体积混凝土结构的中心与其表面　温度之差△　（ｔ）产生的变形不协调所产生的。　温度差值△　为：　ＡＴ（ｔ）＝［　＋Ｔ（ｔ）一　（ｔ）一Ｔｄ（ｔ）］　ｌ　一　（ｔ）ｌｌ￡＝ｔｏ　　（１０）　式中，ｔ　为大体积混凝土结构内部温度降至环周　围境温度时的时间；Ｔｄ（￡）为大体积混凝土非收缩　应力所产生的收缩当量温度；　（ｔ）为大体积混凝　土结构所处的环境温度值；　要控制大体积混凝土结构的温度裂缝，需要　满足以下条件：　Ｒｆ（ｔ）＞　Ｎ（ｔ）＋ｏｒ：（　）＋ｏｒ　（ｔ）＋　ｏ（￡）　一　（ｔ）　（１１）　其中，Ｒ，（ｆ）为大体积混凝土结构在养护过程处于　ｔ龄期时的抗拉强度；　：（ｔ）为外界约束条件下，大　体积混凝土结构降温产生的应力；ｒｏ。（ｆ）为其它环　境因素约束所产生的应力；ｏｒ　（ｔ）为混凝土结构　由于徐变所产生的应力。　２　工程实例　２．１工程概况及计算模型　某圆形基础，外部直径３５　１Ｔｌ，内部直径１０　１ＴＩ，　厚度４　ｍ，属于典型的大体积混凝土，板面标高为　一５．７　ｍ，板底标高为一７．７　ｍ，设置二条后浇式膨　胀加强带，后浇式膨胀加强带宽度为２一　一一加加　　～一．　　ｉｎ，将底板　分成三块。底板后浇式膨胀加强带以外混凝土为　补偿收缩混凝土，强度等级为Ｃ３０，抗渗等级Ｐ６；　后浇式膨胀加强带混凝土强度等级为Ｃ３５，抗渗　等级为Ｐ６，并掺加膨胀剂，在地下室外墙浇筑前　完成，分析模型如图１所示。　图１有限元分析模型　在有限元分析计算过程中，考虑到能量守恒　原理，需要利用数值分析软件ＡＮＳＹＳ热分析来完　成，选取具有三个方向的热传导能力Ｓｏｌｉｄ７０三维　热实体单元，它可以实现匀速热流的传递…　。　该筏板基础大体积混凝土的热力学性能以及物理　力学性能参数如表１。　表１计算参数　名　称　取值　砼材料的导热系数Ａ［ｋＪ／（ｍ・ｈ・ｏＣ）］　土的导热系数Ａ［ｋＪ／（ｍ・ｈ・ｏＣ）］　砼木模热交换系数　。［１（Ｊ／（ｎｌ・ｈ・℃）］　养护过程热交换系数　ｋＪ／（ｍ・ｈ・℃）］　大体积砼初始温度　（℃）　砼周围空气温度　（℃）　线性膨胀系数（ｏｃ）　泊松比　２．２现场测试数据分析　为了保证数值分析的可靠性、准确性和对比　性，首先对工程施工中大体积混凝土的温度变化　和应力场变形进行分析。　为了保证测量方案具有较好的代表性，准确　掌握大体积混凝土结构的温度场变化过程，将该　基础分为上、中、下三组多个温度测试点，以保证　能够同时准确测试不同深度部位的混凝土内部温　度变化过程；具体布置见图２。　／一　＼　）　图２测温点布置　通过现场监测获得的大体积混凝土各个部位　的温度变化曲线如图３～８所示。　图３养护过程中基础Ａ处的温度变化　通过分析大体积混凝土基础的温度变化曲线　和温度应力变化可以看出：在养护的过程中，随着　大体积混凝土水泥水化热反应的结束及混凝土热　量的不断散失，大体积混凝土基础逐渐完成了由　升温阶段向降温阶段的过渡。　浇筑初期大体积混凝土的温度是呈线性变化　上升的，随着养护时间的不断增加，大体积混凝土　第１期　谭广柱等：大体积混凝土温度应力场变化分析　・２３・　赠　∞　如　如　加　ｍ　Ｏ　图４养护过程中基础Ｂ处的温度变化　羹　５０　－一】．　一－ｔ．－一　举器　景　凄　图５养护过程中基础Ｃ处的温度变化　１．０×１　０６　８．０Ｘ１　０５　∞６．０×１　４　０×１０５　２Ｏ×１０　０．０　．２．０×１０　图６养护过程中基础Ａ处的温度应力变化　。＇、、—≥　Ｉ＇４　如　如：　养护过程中基础Ｂ处的温度应力变化　。～　．　。６　１。８　２　２　２。６：　图８养护过程中基础Ｃ处的温度应力变化　基础的中心部分温度迅速增加，并且温度峰值出　现在厚度相对较大的部位。在短时间内，大体积　混凝土温度增加了４０％，此时由于结构边界条件　的约束，在不利组合条件下将可能超过混凝土的　抗拉强度而导致温度裂缝的出现。基于此，在工　程现场进行施工时，应该对大体积混凝土的温度　场和现场气象状况进行连续的实时监测，采取适　合于当地实际情况的养生方式，在养护期间需要　根据天气的变化状况，定期对　昆凝土结构洒水，保　证混凝土内部相对湿度的基本稳定。　２．３大体积混凝土温度场有限元数值分析　为了能够显示大体积混凝土在养护过程中，　其温度应力场的分布变化过程，利用有限元分析　软件ＡＮＳＹＳ对其养护过程进行分析，如图９所　示。　通过分析图ｌＯ～１２可知，在整个养护的过程　中，大体积混凝土的放热升温过程一般是在３～４　天内基本完成，浇筑养护的过程中期温度变化的　升温速率一般可以达到０．４℃／ｈ左右；当温度达　到峰值后一般能够维持１８　ｈ左右，之后的养护过　程中，大体积混凝土开始进入温度下降状态，一般　持续１８天左右，之后的温度场的分布比较均匀，　降温基本趋于稳定。在后期养护过程中，高温集　中区始终存在于大体积混凝土基础的内部，降温　速度比较缓慢的也同样出现在厚度相对较大的部　位。养护实际到达２０天以后，大体积混凝土的温　度场基本上是均匀分布的，并且逐渐趋于稳定。　从空间角度分析发现由于大体积混凝土水化　热效应的显著性导致其应力分布呈现出不均匀性　的特点，这主要是由于大体积混凝土不同部位与　周围环境进行热交换的差异所引起的，基础　上表面与环境的空气发生热交换，热量散失较好，　故其温度曲线明显低于中心和底部温度曲线；而　基础的下表面和中心由于与外界接触的面积　较小，加之边界约束条件的，发生热交换不够　充分，故而温度较高。　通过对比数值分析结果与现场监测的数据可　知：大体积混凝土的温度变化曲线虽然发展变化　趋势吻合较好，但是数值分析的曲线明显比现场　测试的数据曲线要光滑，这主要是由于在有限元　分析中未考虑施工工艺、养护条件的偶然性变化　以及弹性模量、边界条件、热力学公式等参数的选　择过于理想化，导致数值分析与现场实测曲线不　一样，说明ＡＮＳＹＳ模拟分析大体积混凝土温度场　变化趋势是基本可行的，但是仍存在一定的误差。　・４４・　土木工程与管理学报　２０１３正　晏，刘新荣．深基坑开挖对临近建筑　］　］　］　］　ｔｕｒｅｓ　ｄｉｓｐｌａｃｅｍｅｎｔ　ａｎｄ　ｄｅｓｉｇｎ［Ｊ］．Ｇｅｏｔｅｃｈｎｉｑｕｅ，　１９９２，４２（４）：５４１—７６７．　［７］薛］　莲，傅］　物的影响研究［Ｊ］．地下空间与工程学报，２００８，４　Ｉ二　（５）：８４７—８５１．　［３］　Ｏｕ　Ｃ　Ｙ，Ｌａｉ　Ｃ　Ｈ．Ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ　ａｎａｌｙｓｉｓ　ｏｆ　ｄｅｅｐ　ｅｘ—　ｃａｖａｔｉｏｎ　ｉｎ　ｌａｙｅｒｅｄ　ｓａｎｄｙ　ａｎｄ　ｃｌａｙｅｙ　ｓｏｉｌ　ｄｅｐｏｓｉｔｓ［Ｊ］．　Ｃａｎａｄｉａｎ　Ｇｅｏｔｅｃｈｎｉｃａｌ　Ｊｏｕｒｎａｌ，１９９４，３　１（２）：２０４—　２１４．　［８］　李萍，邓小鹏，相建华，等．基坑开挖与支护模　拟的位移迭代法［Ｊ］．岩土力学，２００５，２６（１１）：　１８１５—１８ｌ８．　［４］　刘国彬，王卫东．基坑工程手册（第二版）［Ｍ］．北　京：中国建筑工业出版社，２００９．　［５］　ＪＧＪ　１２００９，建筑基坑支护技术规程［Ｓ］．　［９］　武亚军，栾茂田，任汉锋．深基坑支护设计增量法　的理论分析及其应用［Ｊ］．工业建筑，２００４，３４　（９）：１－４．　［６］　王浩然，王卫东，徐中华．基坑开挖对临近建筑物　影响的三维有限元分析［Ｊ］．地下空间与工程学　报，２００９，５（ｓ２）：１５１２—１５１７．　［１Ｏ］张克恭，刘松玉．土力学［Ｍ］．北京：中国建筑工　业出版社，２００１．　“＋“＋”＋”＋”—・卜“—＿卜”—卜“－４－“—＋－　（上接第２４页）　参考文献　分区异步长算法［Ｊ］．水利发电，１９９６，（１）：３８－４３．　［７］刘光廷，谢树南，辉，等．碾压？昆凝土拱坝设人　朱伯芳．大体积｝昆凝土温度应力与温度控制［Ｍ］．　北京：中国电力出版社，１９９９．　Ｍａｌｋａｗｌ　Ａｂｄａｌｌａｈ　Ｉ　Ｈｕｓｅｉｎ，Ｍｕｔａｓｈｅｒ　Ｓａａｄ　Ａ，Ｑｉｕ　Ｔｏｎｙ　Ｊ．Ｔｈｅｒｍａｌ—ｓｔｒｕｃｔｕｒａｌ　ｍｏｄｅｌｉｎｇ　ａｎｄ　ｔｅｍｐｅｒａｔｕｒｅ　ｃｏｎｔｒｏｌ　工短缝的应力释放及止裂作用［Ｊ］．水利学报，　２００２，（６）：９－１４．　［８］　李守巨，刘迎曦．基于模糊理论的混凝土热力学参　数识别方法［Ｊ］．岩土力学，２００４，２５（４）：５７０—　５７３．　ｏｆ　ｒｏｌｌｅｒ　ｃｏｍｐａｃｔｅｄ　ｃｏｎｃｒｅｔｅ　ｇｒａｖｉｔｙ　ｄａｍ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｐｅｒｆｏｒｍａｎｃｅ　ｏｆ　Ｃｏｎｓｔｕｃｔｅｄ　Ｆａｃｉｌｒｉｔｉｅｓ，２００３，１７　［９］ＧＢ　５０４９６—２００９，大体积｝昆凝土施工规范［Ｓ］．　［１０］刘文燕，耿耀明．热工参数对混凝土结构温度场影　响研究［Ｊ］．混凝土与水泥制品，２００５，（１）：１１—　１５．　（４）：１７７—１８７．　Ｋａｗａｇｕｅｈｉ　Ｔｏｈｒｕ，Ｎａｋａｎｅ　Ｓｕｎａｏ．Ｉｎｖｅｓｔｉｇａｔｉｏｎｓ　ｏｎ　ｄｅ－　ｔｅｒｍｉｎｉｎｇ　ｔｈｅｒｍａｌ　ｓｔｒｅｓｓ　ｉｎ　ｍａｓｓｉｖｅ　ｃｏｎｃｒｅｔｅ　ｓｔｒｕｃｔｕｒｅｓ　［Ｊ］．ＡＣＩ　Ｍａｔｅｒｉａｌｓ　Ｊｏｕｒｎａｌ，１９９６，９３（１）：９６—１０１．　Ｋｉｍ　Ｊａｎｇ—Ｈｏ　Ｊａｙ，Ｊｅｏｎ　Ｓａｎｇ—Ｅｕｎ，Ｋｉｍ　Ｊｉｎ—Ｋｅｕｎ．Ｄｅ—　ｖｅｌｏｐｍｅｎｔ　ｏｆ　ｎｅｗ　ｄｅｖｉｃｅ　ｆｏｒ　ｍｅａｓｕｒｉｎｇ　ｔｈｅｒｍａｌ　ｓｔｒｅｓｓ　［１１］刘［１２］周睫，陈兵．大体积混凝土水化热温度场数值　模拟［Ｊ］．混凝土与水泥制品，２０１０，１０（５）：１５—２７．　莹，王双喜．塑料连栋温室覆盖层的能量平衡　模型的理论研究［Ｊ］．中国农学通报，２０１０，２６　（２０）：４０２—４０５．　［Ｊ］．Ｃｅｍｅｎｔ　ａｎｄ　Ｃｏｎｃｒｅｔｅ　Ｒｅｓｅａｒｃｈ，２００２，３２（１０）：　１６５１—１６５４．　朱伯芳．多层混凝土结构仿真分析的并层算法［Ｊ］．　水力发电学报，１９９４，（３）：２１—３０．　［１３］石亦平，周玉蓉．ＡＢＡＱＵＳ有限元分析实例详解　［Ｍ］．北京：机械工业出版社，２００６．　朱伯芳，许平．混凝土坝仿真计算的并层算法和