王立成，吳 迪，鮑玖文，梁永欽

（大連理工大學(xué) 海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024）

1 研究背景

早齡期混凝土是指28 d齡期內(nèi)的混凝土［1］。非荷載因素造成的早齡期混凝土開(kāi)裂是影響鋼筋混凝土結(jié)構(gòu)耐久性的重要因素［2］。造成其開(kāi)裂的原因主要是因?yàn)槟z凝材料的水化會(huì)引起混凝土內(nèi)部溫度發(fā)生變化，并在內(nèi)外約束下產(chǎn)生溫度應(yīng)力，進(jìn)而引起混凝土開(kāi)裂［3］。因此，研究早齡期混凝土內(nèi)部溫度場(chǎng)的分布規(guī)律，可為混凝土結(jié)構(gòu)的早期抗裂性能評(píng)估提供有利依據(jù)。

近年來(lái)國(guó)內(nèi)外研究者對(duì)早齡期混凝土溫度場(chǎng)分布的研究，主要從試驗(yàn)研究、理論分析和數(shù)值模擬等方面展開(kāi)。Kim等［4］試驗(yàn)研究了在正常條件下溫度、濕度和骨料等對(duì)混凝土導(dǎo)熱系數(shù)的影響，并提出了考慮以上各因素影響的混凝土導(dǎo)熱系數(shù)模型。Wang等［5］基于微觀水化模型，研究了摻加粉煤灰和礦粉的早齡期混凝土溫度分布規(guī)律。理論方面，王甲春等［6］通過(guò)建立混凝土材料的水化放熱模型，引入了Arrhenius公式描述早齡期混凝土放熱速率受溫度的影響規(guī)律，并分析了混凝土在早齡期時(shí)溫度應(yīng)力隨齡期的發(fā)展規(guī)律。張君等［7］建立了基于絕熱溫升的早齡期混凝土溫度場(chǎng)計(jì)算模型，通過(guò)等效齡期法修正了溫度對(duì)水泥水化及其放熱量的影響，并通過(guò)試驗(yàn)驗(yàn)證了模型的合理性。在數(shù)值模擬方面，崔溦等［8］從理論上考慮了早齡期混凝土熱學(xué)參數(shù)隨水化度的變化規(guī)律，更加準(zhǔn)確地模擬了早齡期混凝土的溫度場(chǎng)。Klemczak等［9］研究表明，早齡期混凝土的濕熱擴(kuò)散具有明顯的耦合特性。但是，以上研究工作通常基于混凝土的宏觀層次，認(rèn)為混凝土是單一均質(zhì)材料，未考慮混凝土的非均質(zhì)特性。實(shí)際混凝土各相材料的性能存在差異，且水泥是水化過(guò)程中熱量的主要來(lái)源，忽視混凝土各相材料性能差異而開(kāi)展的數(shù)值仿真分析，將不能準(zhǔn)確描述混凝土的實(shí)際溫度場(chǎng)及熱量分布，進(jìn)而不能有效地解決混凝土的早期開(kāi)裂問(wèn)題。

在細(xì)觀層次上，混凝土可以看作是由粗骨料、砂漿及二者界面過(guò)渡區(qū)（ITZ）組成的三相復(fù)合材料［10］。本文考慮混凝土各相組分熱學(xué)性能的差異，建立早齡期混凝土的熱傳導(dǎo)方程并進(jìn)行數(shù)值離散，進(jìn)而開(kāi)展早齡期混凝土溫度場(chǎng)分布的細(xì)觀數(shù)值仿真分析，最后分析骨料預(yù)冷和邊界條件對(duì)溫度場(chǎng)分布及變化規(guī)律的影響。

2 早齡期混凝土的熱傳輸理論

2.1 熱傳導(dǎo)控制方程 混凝土熱量傳遞主要有3種方式：熱對(duì)流、熱傳導(dǎo)和熱輻射，其中熱傳導(dǎo)是影響混凝土內(nèi)部溫度場(chǎng)分布的主要機(jī)制?；炷恋囊痪S（設(shè)傳熱方向?yàn)閤軸方向）非穩(wěn)態(tài)熱傳導(dǎo)控制方程可表示為［11］：

式中：T（x，t）為經(jīng)歷時(shí)間t后在位置x處的溫度值；D為混凝土導(dǎo)溫系數(shù)，D=λ/（ρc），其中，λ為混凝土導(dǎo)熱系數(shù)，ρ為混凝土密度，c為混凝土比熱容；Q為混凝土水化引起的熱源，即絕熱溫升速率，與水泥的水化度有關(guān)。

2.2 等效齡期與絕熱溫升 混凝土中的水泥在水化過(guò)程中釋放熱量，引起自身溫度升高，而水化反應(yīng)速率隨溫度的升高而加快。通常，溫度對(duì)水化反應(yīng)速率的影響可通過(guò)Arrhenius方程描述［12］：

式中：k（T）為與溫度T有關(guān)的化學(xué)反應(yīng)速率；η為指前因子；Ea為混凝土的活化能；R為氣體常數(shù)。Freiesleben等［13］基于Arrhenius方程提出了混凝土的等效齡期成熟度函數(shù)，其表達(dá)式為：

式中：te為相對(duì)于參考溫度的混凝土等效齡期；Tr為參考溫度。

一般而言，混凝土水化度指某一時(shí)刻混凝土的水化反應(yīng)程度，通常與等效齡期te密切相關(guān)，所以Knudsen［14］在試驗(yàn)基礎(chǔ)上提出了基于等效齡期的雙曲線(xiàn)式絕熱溫升計(jì)算模型，其表達(dá)式為：

式中：α（te）為基于等效齡期te的水化度；C為水化度曲線(xiàn)形狀參數(shù)；θ（α）為基于等效齡期的混凝土絕熱溫升；θu為最終絕熱溫升。

早齡期混凝土的水化反應(yīng)是一個(gè)動(dòng)態(tài)過(guò)程，溫度（熱）與水化度（化學(xué)）之間是互相影響的，混凝土水化會(huì)引起自身溫度的升高，溫度對(duì)其水化反應(yīng)速率產(chǎn)生影響。本文在熱傳導(dǎo)控制方程中引入的水化熱源Q，既考慮了溫度對(duì)其反應(yīng)速率的影響，又可影響早齡期混凝土內(nèi)溫度場(chǎng)分布，這充分體現(xiàn)了水泥水化過(guò)程中的“熱-化學(xué)”耦合作用。

3 早齡期混凝土熱量傳輸?shù)募?xì)觀數(shù)值分析

3.1 細(xì)觀格構(gòu)網(wǎng)絡(luò)模型 在細(xì)觀層次上，混凝土材料具有明顯的非均質(zhì)特性，為考慮混凝土骨料分布的隨機(jī)性和各相組分不同的傳輸性能，可將混凝土看作由粗骨料、砂漿和界面過(guò)渡區(qū)（ITZ）組成的三相復(fù)合材料。根據(jù)利用已知點(diǎn)的集合將平面或者空間劃分成凸多邊形或凸多面體的Voronoi法進(jìn)行單元?jiǎng)澐?，則整個(gè)混凝土區(qū)域被劃分成有限個(gè)Voronoi單元（如圖1）。連接Voronoi單元的中心與其各邊界的中點(diǎn)，形成Voronoi單元內(nèi)部的格構(gòu)單元；連接Voronoi單元邊界的中點(diǎn)和角點(diǎn)，形成邊界格構(gòu)單元。根據(jù)該劃分原則，混凝土細(xì)觀模型將形成兩類(lèi)內(nèi)部格構(gòu)單元（骨料和砂漿格構(gòu)單元）和三類(lèi)邊界格構(gòu)單元（砂漿-砂漿、骨料-骨料、砂漿-骨料）。

對(duì)于混凝土在一定面積上的熱傳導(dǎo)問(wèn)題，可用格構(gòu)單元來(lái)代替，從而將二維的熱傳導(dǎo)問(wèn)題簡(jiǎn)化為一維問(wèn)題。各相組分的傳輸性能與其格構(gòu)單元的截面面積和傳輸系數(shù)有關(guān)，規(guī)定如下：內(nèi)部格構(gòu)單元的傳輸系數(shù)為常數(shù)，同時(shí)假定混凝土厚度為單位長(zhǎng)度，因此該格構(gòu)單元截面面積數(shù)值上可取為Voronoi單元邊長(zhǎng)的一半，例如，格構(gòu)單元ij的截面面積為mn長(zhǎng)度的一半；砂漿-骨料邊界格構(gòu)單元的傳輸系數(shù)為砂漿格構(gòu)單元的10倍［15］，截面面積取決于界面層的厚度，通常取為20μm［15］；砂漿-砂漿和骨料-骨料邊界格構(gòu)單元的傳輸系數(shù)為常數(shù)，截面面積為0。

細(xì)觀格構(gòu)網(wǎng)絡(luò)模型的計(jì)算效率與Voronoi單元大小及時(shí)間步長(zhǎng)有關(guān)，在確保計(jì)算精度的條件下，通過(guò)增大Voronoi單元或時(shí)間步長(zhǎng)，均可提高計(jì)算效率。

圖1 混凝土細(xì)觀格構(gòu)網(wǎng)絡(luò)模型

3.2 數(shù)值求解及模型驗(yàn)證 混凝土結(jié)構(gòu)暴露在空氣中，與周?chē)h(huán)境進(jìn)行熱交換，屬于第三類(lèi)邊界條件（對(duì)流邊界條件）［16］，該條件認(rèn)為混凝土的表面熱流量與表面溫度T和大氣氣溫Ta之差成正比，其表達(dá)式為：

式中：β為放熱系數(shù)，kJ/（m2·h·℃）。

對(duì)于圖1所示長(zhǎng)度為L(zhǎng)的一維格構(gòu)單元ij，采用Galerkin法和Crank-Nicholson差分法，以Δt為時(shí)間步長(zhǎng)，對(duì)熱傳導(dǎo)控制方程進(jìn)行數(shù)值求解。

假定格構(gòu)單元節(jié)點(diǎn)的溫度分別為T(mén)i、Tj，加權(quán)余函數(shù)為Φi=（xj-x）/L和Φj=（x-xi）/L，則溫度T的任意函數(shù)f（T）沿格構(gòu)單元ij的加權(quán)余函數(shù)可表示為f（T）=Φif（Ti）+Φjf（Tj），對(duì)式（1）通過(guò)積分并整理，可確定格構(gòu)單元ij的有限元方程和差分離散形式如下：

式中：A為格構(gòu)單元的截面面積；w為與維數(shù)有關(guān)的參數(shù)，對(duì)于平面問(wèn)題，w取為2［17］；Q為水化引起的熱源，對(duì)于骨料和邊界格構(gòu)單元，Q=0；對(duì)于砂漿格構(gòu)單元，在時(shí)間段［t，t+Δt］內(nèi)的絕熱溫升速率為：

為了驗(yàn)證所建立的細(xì)觀格構(gòu)網(wǎng)絡(luò)模型在早齡期混凝土熱傳導(dǎo)問(wèn)題中應(yīng)用的可行性，結(jié)合張子明等［18］開(kāi)展的養(yǎng)護(hù)溫度對(duì)混凝土絕熱溫升影響的試驗(yàn)研究結(jié)果，開(kāi)展早齡期混凝土絕熱升溫的數(shù)值計(jì)算。文獻(xiàn)［18］的絕熱溫升試驗(yàn)采用3種不同的初始溫度（4.4℃、23.3℃和40℃），在砂漿拌合1 h內(nèi)，將試件置于水池中。采用溫度控制精度在±0.05℃的3 kW加熱器來(lái)調(diào)整控制水的溫度，借助熱電偶監(jiān)測(cè)試件的溫度，并通過(guò)多通道數(shù)據(jù)記錄器每分鐘測(cè)量一次溫度并輸出每天的平均溫度。采用圖2所示的二維細(xì)觀幾何模型，骨料含量為40%并隨機(jī)分布，水灰比0.45，骨料和砂漿單元的計(jì)算參數(shù)與試驗(yàn)條件保持一致，界面過(guò)渡區(qū)的比熱容、密度與砂漿單元一致，導(dǎo)熱系數(shù)λ=20 W/（m·K），主要參數(shù)的取值情況見(jiàn)表1。

圖2 混凝土二維計(jì)算模型

表1 混凝土絕熱溫升計(jì)算參數(shù)

圖3 不同初始溫度下早齡期混凝土絕熱溫升曲線(xiàn)

采用細(xì)觀格構(gòu)網(wǎng)絡(luò)模型，分別計(jì)算得到了試件中心A處混凝土在不同初始溫度下考慮與不考慮“熱-化學(xué)”耦合的早齡期混凝土絕熱溫升曲線(xiàn)，與試驗(yàn)結(jié)果的對(duì)比如圖3所示。由圖3可知，在不同初始溫度下，不考慮“熱-化學(xué)”耦合的數(shù)值模擬計(jì)算結(jié)果明顯低于試驗(yàn)值，這是因?yàn)闇囟壬邥?huì)導(dǎo)致水化速率加快，進(jìn)而導(dǎo)致絕熱溫升速率加快。當(dāng)考慮“熱-化學(xué)”耦合時(shí)，數(shù)值模擬計(jì)算結(jié)果與試驗(yàn)值吻合較好，后期溫度變化略高于試驗(yàn)值。產(chǎn)生偏差的原因，一是數(shù)值計(jì)算中采用經(jīng)驗(yàn)性的導(dǎo)溫系數(shù)，可能與試驗(yàn)中混凝土試件的導(dǎo)溫系數(shù)存在一定偏差；二是文獻(xiàn)［18］的絕熱試驗(yàn)條件不可能達(dá)到理論上的完全絕熱，從而使得試驗(yàn)中存在一定的熱量散失，造成試驗(yàn)結(jié)果低于理論計(jì)算值。從整體上來(lái)看，混凝土絕熱溫升數(shù)值計(jì)算結(jié)果與試驗(yàn)值變化趨勢(shì)基本一致，這說(shuō)明通過(guò)細(xì)觀格構(gòu)網(wǎng)絡(luò)模型可實(shí)現(xiàn)考慮“熱-化學(xué)”耦合作用的早齡期混凝土內(nèi)部熱傳導(dǎo)問(wèn)題的數(shù)值仿真分析。

3.3 骨料預(yù)冷效果分析 骨料預(yù)冷是大體積水工混凝土夏季溫度控制的常用措施，一般將骨料預(yù)冷至0～10℃，通過(guò)低溫骨料吸收水泥水化釋放的熱量，進(jìn)而達(dá)到降低混凝土內(nèi)部溫度的效果。分別采用圖4所示的單骨料混凝土和多骨料混凝土的細(xì)觀格構(gòu)網(wǎng)絡(luò)模型，針對(duì)混凝土的骨料預(yù)冷效果開(kāi)展相對(duì)應(yīng)的細(xì)觀數(shù)值仿真分析。假定混凝土水灰比為0.45，骨料含量約為27%，砂漿初始溫度為40℃，骨料初始溫度為5℃，邊界絕熱，骨料導(dǎo)熱系數(shù)為3 W/（m·K），密度為1500 kg/m3，比熱容為1000 J/（kg·K）。僅有砂漿水化釋放熱量，其他計(jì)算參數(shù)如表1所示，界面過(guò)渡區(qū)的導(dǎo)熱系數(shù)λ=20 W/（m·K），比熱容、密度與砂漿單元一致。

圖5為單骨料混凝土中心A處溫度歷程曲線(xiàn)以及溫度沿水平向的分布曲線(xiàn)和溫度分布云圖。從圖5可知，在澆筑初期，骨料溫度急劇上升至30℃，砂漿溫度逐漸降低，直至骨料與砂漿溫度一致，此時(shí)砂漿溫度大約降低10℃（如圖5（b）），這主要是由砂漿與預(yù)冷骨料間的熱傳導(dǎo)造成的；此后，隨著砂漿水化，混凝土溫度逐漸上升，28 d后溫度約為50℃。圖5（c）給出了水化1、2和5 min后混凝土內(nèi)溫度分布云圖，說(shuō)明建立的細(xì)觀格構(gòu)網(wǎng)絡(luò)可直觀地看出各位置處溫度的變化情況。

圖4 混凝土二維計(jì)算模型（單元數(shù)量：Voronoi941；格構(gòu)10 704）

多骨料混凝土的溫度分布曲線(xiàn)和云圖如圖6所示。從圖6可以看出，在澆筑10 min后，混凝土內(nèi)部溫度趨于一致，砂漿溫度降低了10℃。對(duì)比圖5（b）可知，在骨料體積含量相同的條件下，混凝土整體降溫效果與骨料粒徑無(wú)關(guān)。由圖5、圖6可知，骨料預(yù)冷條件下的混凝土在澆筑初期內(nèi)部溫度分布不均。若采用宏觀數(shù)值模擬理論，假定混凝土為均質(zhì)材料，則無(wú)法真實(shí)反應(yīng)混凝土內(nèi)部的溫度分布情況。在細(xì)觀數(shù)值模擬方法中，假定骨料和砂漿具有不同的物理性能，可真實(shí)模擬骨料預(yù)冷條件下的混凝土溫度場(chǎng)。

3.4 散熱邊界影響分析 采用圖7所示的混凝土細(xì)觀格構(gòu)網(wǎng)絡(luò)模型（500mm×150mm），針對(duì)散熱邊界影響下的溫度場(chǎng)分布開(kāi)展細(xì)觀仿真分析?；炷了冶葹?.45，骨料含量約為14%，混凝土初始溫度26℃，周?chē)h(huán)境溫度為20℃，底部為散熱邊界，放熱系數(shù)為18 kJ/（m2·h·℃），骨料和界面過(guò)渡區(qū)的計(jì)算參數(shù)與3.3節(jié)一致，僅有砂漿水化釋放熱量，其他計(jì)算參數(shù)如表1所示。

圖6 多骨料混凝土水平向溫度分布曲線(xiàn)

圖7 混凝土計(jì)算模型（單位：mm）

圖8 溫度發(fā)展歷程曲線(xiàn)

圖9 散熱邊界下溫度分布云圖

散熱邊界下混凝土溫度發(fā)展歷程曲線(xiàn)和溫度分布云圖分別如圖8和圖9所示。由圖可知，混凝土內(nèi)部各點(diǎn)溫度在澆筑初期迅速上升，約在第2 d達(dá)到溫度峰值，最高溫度約為34℃，隨后溫度逐步下降，約在第28 d達(dá)到穩(wěn)定，即與外界溫度相等。另外，混凝土澆筑后自底部至3/4高度范圍內(nèi)的溫度相差不大，這主要是由于混凝土通過(guò)散熱邊界與外界進(jìn)行熱交換導(dǎo)致的，一般來(lái)說(shuō)，散熱邊界越大，其影響的高度范圍也越大。

4 結(jié)論

（1）基于混凝土的細(xì)觀結(jié)構(gòu)，提出了適用于混凝土內(nèi)部熱傳導(dǎo)的細(xì)觀格構(gòu)網(wǎng)絡(luò)模型，建立了考慮“熱-化學(xué)”耦合的一維非穩(wěn)態(tài)熱傳導(dǎo)控制方程并進(jìn)行數(shù)值求解。通過(guò)與文獻(xiàn)中的試驗(yàn)結(jié)果對(duì)比，驗(yàn)證了細(xì)觀格構(gòu)網(wǎng)絡(luò)模型在考慮“熱-化學(xué)”耦合作用的早齡期混凝土內(nèi)部熱傳導(dǎo)問(wèn)題中應(yīng)用的可行性。（2）開(kāi)展了骨料預(yù)冷和散熱邊界對(duì)早齡期混凝土溫度場(chǎng)分布影響的數(shù)值模擬研究，結(jié)果表明，預(yù)冷骨料可在短期內(nèi)迅速降低混凝土的整體溫度，且在骨料含量相同的條件下，降溫效果與骨料粒徑無(wú)關(guān)；散熱邊界會(huì)影響混凝土溫度場(chǎng)分布，距離散熱邊界越近，峰值溫度越低。

致謝：感謝日本北海道大學(xué)的上田多門(mén)教授在混凝土細(xì)觀數(shù)值方法方面給予的幫助和指導(dǎo)。

參 考 文 獻(xiàn)：

［1］ 陳軍.早齡期混凝土水化進(jìn)程及宏觀與細(xì)微觀性能相關(guān)性研究［D］.杭州：浙江大學(xué)，2014.

［2］ 王鐵夢(mèng).工程結(jié)構(gòu)裂縫控制［M］.北京：中國(guó)建筑工業(yè)出版社，1997.

［3］ KOVLER K，LANGE D A，STANG H.Early age concrete—properties and performance［J］.Cement and Con?crete Composites，2004，26（5）：413-415.

［4］ KIM K H，JEON S E，KIM J K，et al.An experimental study on thermal conductivity of concrete［J］.Cement and Concrete Research，2003，33（3）：363-371.

［5］ WANG X Y，CHO H K，LEE H S.Prediction of temperature distribution in concrete incorporating fly ash or slag using a hydration model［J］.Composites Part B Engineering，2011，42（1）：27-40.

［6］ 王甲春，閻培渝.早齡期混凝土結(jié)構(gòu)的溫度應(yīng)力分析［J］.東南大學(xué)學(xué)報(bào)，2005，35（S1）：15-18.

［7］ 張君，祁錕，侯東偉.基于絕熱溫升試驗(yàn)的早齡期混凝土溫度場(chǎng)的計(jì)算［J］.工程力學(xué)，2009，26（8）：155-160.

［8］ 崔溦，吳甲一，宋慧芳.考慮水化度對(duì)熱學(xué)參數(shù)影響的早期混凝土溫度場(chǎng)分析［J］.東南大學(xué)學(xué)報(bào)，2015，45（4）：792-798.

［9］ KLEMCAZK B.Prediction of coupled heat and moisture transfer in early-age massive concrete structures［J］.Nu?merical Heat Transfer Part A：Applications，2011，60（3）：212-233.

［10］ 王立成，鮑玖文.干濕交替下開(kāi)裂混凝土中水分傳輸?shù)募?xì)觀數(shù)值分析［J］.水利學(xué)報(bào)，2016，47（8）：986-995.

［11］ 何志霞，王謙，袁建平.熱流體數(shù)值計(jì)算方法與應(yīng)用［M］.北京：機(jī)械工業(yè)出版社，2013.

［12］ 田野，金賢玉，金南國(guó).基于水泥水化動(dòng)力學(xué)和等效齡期法的混凝土溫度開(kāi)裂分析［J］.水利學(xué)報(bào)，2012，43（S1）：179-186.

［13］ HANSEN P F，PEDERSEN E J.Maturity computer for controlled curing and hardening of concrete［J］.Nordisk Betong，1977，1（19）：19-34.

［14］ KNUDSEN T.Proceedings of the Engineering Foundation Conference on Characterization and Performance Predic?tion of Cement and Concrete［C］.New Hampshire：United Engineering Trustees，1982.

［15］ OH B H，JANG S Y.Prediction of diffusivity of concrete based on simple analytic equations［J］.Cement and Con?crete Research，2004，34（3）：463-480.

［16］ 王寧.考慮性態(tài)變化的早期混凝土多場(chǎng)耦合分析及其應(yīng)用［D］.天津：天津大學(xué)，2013.

［17］ NAKAMURA H，SRISOROS W，YASHIRO R，et al.Time-dependent structural analysis considering mass trans?fer to evaluate deterioration process of RC structures［J］.Journal of Advanced Concrete Technology，2006，4（1）：147-158.

［18］ 張子明，宋智通，黃海燕.混凝土絕熱溫升和熱傳導(dǎo)方程的新理論［J］.河海大學(xué)學(xué)報(bào)，2002，30（3）：1-6.