基于熱流量積分的混凝土溫控水管冷卻邊界模擬算法

朱振泱，劉敏芝，強(qiáng) 晟，相建方

（1. 中國(guó)水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100038；2. 鄞州區(qū)水利工程質(zhì)量與安全監(jiān)督站，寧波 315100； 3. 河海大學(xué)水利水電工程學(xué)院，南京 210098）

基于熱流量積分的混凝土溫控水管冷卻邊界模擬算法

朱振泱1，劉敏芝2，強(qiáng) 晟3，相建方1

（1. 中國(guó)水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100038；2. 鄞州區(qū)水利工程質(zhì)量與安全監(jiān)督站，寧波 315100； 3. 河海大學(xué)水利水電工程學(xué)院，南京 210098）

以往用冷卻水管離散模型迭代計(jì)算混凝土溫度場(chǎng)時(shí)，水管邊界被近似認(rèn)為是第三類邊界條件，此類邊界條件的參數(shù)獲取要通過試驗(yàn)并進(jìn)行反演，試驗(yàn)費(fèi)用較大且有時(shí)候可靠性不高。針對(duì)該問題，在熱量平衡條件的基礎(chǔ)上提出由與水管接觸混凝土的熱流量、水管的導(dǎo)熱系數(shù)、管壁厚度和水管內(nèi)壁溫度推算水管外壁溫度的新算法，并對(duì)原有迭代方法進(jìn)行改進(jìn)，解決迭代次數(shù)多和迭代可能無法收斂的問題。對(duì)比算例中，采用傳統(tǒng)算法，模型邊界處的誤差可達(dá)到1.67℃，而采用該新算法，誤差僅為0.3℃。應(yīng)用所提方法對(duì)某混凝土塊施工期混凝土溫度場(chǎng)進(jìn)行了仿真計(jì)算，計(jì)算值與實(shí)測(cè)值吻合較好，且迭代收斂速度較快，一般的迭代方法，需要迭代15次，而采用改進(jìn)的迭代方法，只需迭代7次即可以達(dá)到穩(wěn)定值。該算法能明確通水冷卻的邊界條件，節(jié)省試驗(yàn)費(fèi)用，提高計(jì)算效率，有較好的工程應(yīng)用價(jià)值。

溫度；模型；混凝土；迭代算法

0 引言

溫度應(yīng)力導(dǎo)致的裂縫為混凝土施工期最常見的裂縫，近年來很多學(xué)者對(duì)此進(jìn)行了持續(xù)深入的研究。冷卻水管是控制混凝土溫度應(yīng)力的有效措施，近年來廣泛應(yīng)用于工程實(shí)踐中[1-4]。

目前普遍認(rèn)為，金屬水管內(nèi)壁和外壁的溫差可以忽略，但塑料水管必須考慮內(nèi)外壁溫差的影響。塑料水管內(nèi)部的溫度場(chǎng)分布較為復(fù)雜，雖然有學(xué)者認(rèn)為塑料管內(nèi)溫度可以等效為線性分布，即塑料管內(nèi)的溫度梯度為常數(shù)[5-6]；但有關(guān)文獻(xiàn)和本文的研究成果均表明塑料水管內(nèi)部的溫度梯度非線性分布，不應(yīng)用常數(shù)表示[7]。以往含離散水管模型的大體積混凝土溫度場(chǎng)迭代計(jì)算中，一般將塑料水管邊界作為第三類邊界條件，在這種條件下溫度場(chǎng)計(jì)算精度與管壁放熱系數(shù)的取值有密切關(guān)系。不同的水管材質(zhì)、管壁厚度、水的流速、沿程水溫條件下，管壁放熱系數(shù)均有可能不同。由于管壁放熱系數(shù)對(duì)溫度場(chǎng)的影響較大，目前往往要通過現(xiàn)場(chǎng)試驗(yàn)反演計(jì)算獲得，耗費(fèi)的成本較大[8]。且目前用等效的放熱系數(shù)一個(gè)參數(shù)來涵蓋管壁導(dǎo)熱系數(shù)、管壁厚度和管壁表面溫度梯度等多個(gè)物理量綜合產(chǎn)生的效果，有時(shí)候可靠性不高。

目前大體積混凝土通水冷卻模型方面的研究較為豐富，許多專家都提出了各自的模型:等效算法[7]、子結(jié)構(gòu)算法[9-10]、埋置單元法[11]、模擬混凝土水管冷卻效應(yīng)的直接算法[12]、擴(kuò)展有限元法[13]。目前用于含水管大體積混凝土溫度場(chǎng)的算法中，如不考慮水管周圍混凝土的溫度場(chǎng)，則計(jì)算量小，如朱伯芳的等效算法在混凝土壩的溫控防裂中應(yīng)用廣泛[7]。離散水管模型的混凝土溫度場(chǎng)迭代算法是由朱伯芳院士提出，并由朱岳明教授進(jìn)行改進(jìn)并推廣，目前已經(jīng)廣泛應(yīng)用于各種薄壁結(jié)構(gòu)中，但在大壩中仍屬于初步的應(yīng)用[7,14-17]。該計(jì)算方法計(jì)算精度較高，且能很好地模擬工程現(xiàn)場(chǎng)的“蛇形水管”布置形式。但是節(jié)點(diǎn)數(shù)量龐大，計(jì)算速度慢，對(duì)計(jì)算機(jī)的性能要求較高，目前多用于諸如水閘等薄壁結(jié)構(gòu)中。由于混凝土通水冷卻時(shí)，塑料水管對(duì)混凝土薄壁結(jié)構(gòu)溫度場(chǎng)和應(yīng)力場(chǎng)影響較大，只有準(zhǔn)確考慮其影響才能對(duì)薄壁結(jié)構(gòu)溫度場(chǎng)和應(yīng)力場(chǎng)進(jìn)行可靠分析[18-19]。故針對(duì)塑料水管內(nèi)溫度和溫度梯度分布以及塑料管對(duì)混凝土冷卻效果的進(jìn)一步研究是十分有必要的。

本文根據(jù)熱平衡原理，推導(dǎo)了水管內(nèi)外壁的溫度關(guān)系式。由沿程水溫（沿水流方向水管內(nèi)壁的溫度）、與水管接觸混凝土的熱流量、管壁導(dǎo)熱系數(shù)及管壁厚度即可以推算出水管外壁溫度（水管與混凝土的接觸面）。由于水管和混凝土緊密接觸，可以認(rèn)為接觸面的熱阻很小，故在計(jì)算中將水管外壁的溫度作為第一類邊界條件。由于水管的導(dǎo)熱系數(shù)和水管的材質(zhì)有很大關(guān)系，同種材質(zhì)的水管，導(dǎo)熱系數(shù)相差不大，而有些水管在出廠時(shí)就標(biāo)明了導(dǎo)熱系數(shù)。管壁的厚度可以通過簡(jiǎn)單測(cè)量獲取。因此水管的導(dǎo)熱系數(shù)和管壁的厚度均可輕易獲取，而沿程水溫、與水管接觸混凝土的熱流量也可以在計(jì)算中獲取，因此該邊界條件處理方法可以避免將管壁作為第三類邊界條件而需要通過試驗(yàn)和反演獲取表面放熱系數(shù)的缺陷，提高了經(jīng)濟(jì)性及可靠性。

目前的最新研究成果表明，采用離散迭代算法時(shí)，即使在網(wǎng)格數(shù)量較少的情況下，如選擇最佳熱流量斷面進(jìn)行熱流量積分，依然能得到熱流量的精確解答[20-22]。以上研究成果表明，采用熱流量積分法計(jì)算塑料水管通水冷卻作用不需要增加額外網(wǎng)格節(jié)點(diǎn)的數(shù)量，是完全可行的。但應(yīng)用該邊界條件處理方法和傳統(tǒng)的水管冷卻離散迭代算法計(jì)算溫度場(chǎng)，會(huì)增加迭代次數(shù)，甚至出現(xiàn)迭代不收斂的現(xiàn)象，需要改進(jìn)?；谝陨戏治?，本文對(duì)混凝土溫控冷卻水管邊界模擬方法進(jìn)行深入研究，擬提出一種精確快速計(jì)算方法解決這些問題。

1 基本理論和方法

1.1 不穩(wěn)定溫度場(chǎng)的計(jì)算原理

在混凝土計(jì)算域R內(nèi)任何一點(diǎn)處，不穩(wěn)定溫度場(chǎng)T滿足以下熱傳導(dǎo)控制方程，

式中T為溫度，℃；α為導(dǎo)溫系數(shù)，m2/d；θ混凝土熱量釋量，℃；t和τ分別表示時(shí)間和齡期，d；R為計(jì)算域。

1.2 水管沿程水溫計(jì)算

根據(jù)熱傳導(dǎo)平衡條件，水管內(nèi)沿程水溫增量可以用以下方程表示[7,15]

式中cλ為混凝土的導(dǎo)熱系數(shù)，kJ/m·d·℃；cw為冷卻水比熱，kJ/kg·℃；ρw為冷卻水的密度，kg/m3；qw為冷卻水流量，m3/d；n為熱流量積分面的外法向向量。

1.3 熱流量積分的最佳斷面

本文所提的方法依賴于熱流量的計(jì)算，如果僅采用式（2）計(jì)算熱流量，在沒有選擇合適的積分?jǐn)嗝娴那闆r下，只有在網(wǎng)格十分密集的情況下，才能得到精確解，需要大量的計(jì)算節(jié)點(diǎn)。同樣，采用式（2）計(jì)算沿程水溫，需要十分密集的網(wǎng)格才能得到精確解，如果采用稀疏的網(wǎng)格計(jì)算沿程水溫則計(jì)算結(jié)果將明顯小于實(shí)際沿程水溫。為此，采用最佳熱流量積分法可解決該問題[20-21]，基本原理如下：

如將S1面定義為水管和混凝土單元A的交界面，S2面定義為混凝土單元A和另一個(gè)混凝土單元的交界面。將ζ=?1和ζ=1分別代表S1面和S2面；根據(jù)文獻(xiàn)[20-21]的研究結(jié)果，對(duì)于8節(jié)點(diǎn)的六面體單元，S1面的溫度梯度計(jì)算值小于真實(shí)值，而S2面的溫度梯度計(jì)算值大于真實(shí)值，只有S3面（ζ=0.4）的溫度梯度的計(jì)算值與真實(shí)值相符。因此，只有采用S3面（ζ=0.4）作為熱流量積分，才能得到真實(shí)的熱流量，為熱流量積分的最佳斷面。

2 含水管離散模型的混凝土溫度場(chǎng)迭代計(jì)算方法

2.1 水管內(nèi)壁和外壁溫差計(jì)算

圖1為水管垂直于水流方向截面的一部分，該部分的弧度為θ，內(nèi)壁弧長(zhǎng)為la，水管順?biāo)鞣较虻拈L(zhǎng)度為dn，水管的內(nèi)徑為ra，水管的外徑為rd，rx是某點(diǎn)到水管中心的距離。Na為該部分水管內(nèi)壁的溫度梯度，Nx是該部分水管內(nèi)某點(diǎn)（與水管中心的距離為rx）的溫度梯度，λ為水管材料的導(dǎo)熱系數(shù)，由于水和混凝土溫差較大，故可以認(rèn)為在水管的任意位置溫度梯度方向均與水管壁垂直，故在dt時(shí)間內(nèi)通過該段水管內(nèi)壁和距離水管中心長(zhǎng)度為rx截面（圖1中的X截面）的熱流量分別為：

圖1 一個(gè)典型水管垂直于水流方向截面Fig.1 A pipe section in direction perpendicular to water flow

同樣，由于水和混凝土溫差較大，可以認(rèn)為在水管的任意位置溫度梯度方向均與水管壁垂直，故dt時(shí)間內(nèi)流過該部分水管內(nèi)壁和圖1中流過X截面的熱量相等，故有。因此，，設(shè)水管內(nèi)壁的溫度為Tin，水管外壁的溫度為Tou。根據(jù)積分，則有，故

式中ra為水管的內(nèi)徑，m；rd為水管的外徑，m；Na為該水管微段水管內(nèi)壁的溫度梯度，℃/m；Tin和Tou分別為內(nèi)壁與外壁的溫度，℃。

計(jì)算時(shí)將水管分為m段，設(shè)第i段水管長(zhǎng)度為dn，單位時(shí)間dt內(nèi)由水管附近混凝土傳給該段水管的熱量為Q。根據(jù)熱流量守恒定則，Q與dt時(shí)間內(nèi)流經(jīng)水管內(nèi)壁的熱量相等，故有，故

式中λ為水管材料的導(dǎo)熱系數(shù)，kJ/m·d·℃；dn為計(jì)算微段水管長(zhǎng)度，m；Q為時(shí)間dt內(nèi)由水管附近混凝土傳給該段水管的熱量kJ/m；la為內(nèi)壁弧長(zhǎng)，m。

故水管內(nèi)壁和水管外壁的溫差為

2.2 沿水流方向水管外壁溫度分布

由于冷卻水的入口溫度已知，利用公式（2）在最佳的熱流量積分?jǐn)嗝嫔献鰺崃髁糠e分，對(duì)每根水管沿水流方向可以逐段推求出沿程水溫。計(jì)算時(shí)根據(jù)需要將水管分為m段，則第i段所在位置的沿程水溫Twi為式中進(jìn)口水溫為Tw0，℃；ΔTwj為沿程水管微段水溫的增量，℃。

水管內(nèi)壁直接與水接觸，因此可以認(rèn)為水管內(nèi)壁的溫度即為水溫。由公式（5）和公式（6）可知水管外壁的溫度Touwi為

式中ΔTi為水管內(nèi)壁和外壁的溫差，℃。

在式（2）中，水管沿程水溫的變化與溫度梯度有關(guān)，是一個(gè)邊界非線性問題，溫度場(chǎng)無法一步求出，必須采用迭代求解法逐步逼近真實(shí)解[7,15]。第一次迭代時(shí)可以假定整個(gè)冷卻水管沿程的內(nèi)壁和外壁的溫度均等于冷卻水的入口溫度，利用式（1）求出溫度場(chǎng)的近似解，再利用式（2）、式（6）和式（7）求出沿水流方向水管內(nèi)壁的溫度分布和沿水流方向水管外壁的溫度分布，重復(fù)以上過程直到獲得穩(wěn)定解。

當(dāng)管壁材料為鐵管時(shí)，沿程水溫十分接近鐵管外壁的溫度，采用上述方法計(jì)算大體積混凝土溫度場(chǎng)，一般迭代7次以內(nèi)就可以獲得穩(wěn)定解。而當(dāng)水管材料為塑料時(shí)，水管內(nèi)外壁溫差很大，甚至可達(dá)到5.0℃以上。沿程水溫變化對(duì)管外壁溫度的影響很大，只有精確求出水管沿程水溫才能獲得精確解。采用上述方法迭代，計(jì)算結(jié)果將迭代很多次以后才能收斂，甚至無法收斂。

2.3 溫度場(chǎng)的迭代算法改進(jìn)

由于第一次迭代前水管內(nèi)外壁的溫度均等于水管進(jìn)口溫度，小于穩(wěn)定解。因第一次迭代計(jì)算出的水管內(nèi)外壁的溫度梯度將大于穩(wěn)定解，故第一次迭代計(jì)算出的水管內(nèi)外壁的溫度也將高于穩(wěn)定解。而第二次迭代前水管內(nèi)外壁溫度即為第一次迭代后的水管內(nèi)外壁溫度，故第二次迭代計(jì)算出的水管內(nèi)外壁溫度和溫度梯度都將低于穩(wěn)定解。同理，第奇數(shù)次迭代計(jì)算出的水管內(nèi)外壁溫度和溫度梯度都將大于穩(wěn)定解，而第偶數(shù)次迭代計(jì)算出的水管內(nèi)外壁溫度和溫度梯度都將小于穩(wěn)定解。

設(shè)第（N-1）次迭代后的管外壁溫度為Tn1?，第N次迭代后的管外壁溫度為Tn，作為第（N+1）次迭代前的管外壁溫度，則計(jì)算結(jié)果更容易接近于穩(wěn)定解，收斂的速度也將提高。一般情況下，采用改進(jìn)后的迭代算法計(jì)算大體積混凝土溫度場(chǎng)時(shí)，迭代8次以內(nèi)誤差即可以控制在0.1℃以內(nèi)。

2.4 迭代收斂條件的數(shù)學(xué)依據(jù)

當(dāng)計(jì)算到第i步的第n次迭代結(jié)束時(shí)，根據(jù)式（3），以下關(guān)系式成立，則

式中Toun和Tin分別表示第n次迭代計(jì)算結(jié)束后某個(gè)水管微段的水管外壁和內(nèi)壁的溫度；Nn表示第n次迭代計(jì)算結(jié)束后得到的水管內(nèi)壁的溫度梯度。

以往的研究成果表明，傳統(tǒng)的離散迭代算法具有較好的收斂性[7,15]；此處只要證明采用該新方法不會(huì)增大迭代難度，即可證明該新方法也具有較好的迭代收斂性能。Tin雖然在迭代過程中也是個(gè)變量，由于該處只分析新方法對(duì)迭代的影響，故假設(shè)Tin是恒定的，那么

設(shè)Tour為真實(shí)的水管外壁溫度。如果將第n次迭代的計(jì)算結(jié)果作為第（n+1）次迭代的初始條件，那么，根據(jù)第2.3節(jié)的分析結(jié)果，當(dāng)當(dāng)。因此，由于當(dāng)時(shí)，是負(fù)數(shù)；而時(shí)，是正數(shù)。故如果足夠小，Tn+1總是向Tr逐漸逼近的。

3 算法對(duì)比分析

本節(jié)將傳統(tǒng)算法、新算法分別和理論解對(duì)比（將水管也作為實(shí)體單元剖分的有限元法），并分析新算法具有更高的精度。

傳統(tǒng)的模型用等效放熱系數(shù)模擬水管的邊界條件且未采用合理的熱流量積分。由于未考慮熱流量積分造成沿程水溫偏低，目前許多工程采用較小的表面放熱系數(shù)模擬水管邊界，如采用β=3 000 kJ/(m·d·℃)模擬（甚至更?。?，此時(shí)會(huì)造成計(jì)算上的誤差（例如進(jìn)出水口附近冷卻效果計(jì)算不準(zhǔn)確等）。

該算例為1.5 m×1.5 m×5 m的混凝土塊，中間埋置冷卻水管（如圖2所示）。澆筑后即開始冷卻，除水管邊界外，其余邊界條件均為絕熱邊界，水溫15℃?；炷两^熱溫升為，等效放熱系數(shù)為β=3 000 kJ/(m·d·℃)，水管的內(nèi)直徑和外直徑分別為28和32 mm。

圖2 有限單元網(wǎng)格Fig.2 Finite element mesh

通水冷卻10 d進(jìn)口附近切面的溫度分布如圖3所示。根據(jù)計(jì)算結(jié)果，離混凝土較遠(yuǎn)的水管絕熱邊界處，該算例下傳統(tǒng)算法的誤差可達(dá)1.67℃，而新算法的誤差僅為0.3℃。

圖3 澆筑10 d時(shí)斷面的溫度分布Fig.3 Temperature distribution on section 10 d after placement

4 工程算例

4.1 基本資料

某混凝土試驗(yàn)塊于2009年9月份施工，該混凝土塊長(zhǎng)27.0 m，寬15.0 m，高3.0 m。分兩層澆筑，每層高度為1.5 m，間歇時(shí)間為4 d。

澆筑溫度控制在18.0℃。距離地面0.75和2.25 m處分別布置兩根塑料冷卻水管，水管的水平間距1.5 m，各層混凝土開始澆筑后1.75 d通水。通水時(shí)間21.0 d，冷卻水初溫為8.5℃，通水12.0 d后水溫調(diào)整為15.0℃，通水流量均為35.0 m3/d。計(jì)算中，氣溫取實(shí)測(cè)值，澆筑塊表面覆蓋一層大壩保溫被。澆筑塊及地基的有限單元網(wǎng)格如圖4所示。特征點(diǎn)平面位置、水管水平布置和水流方向如圖5所示，水流方向始終未變。

圖4 有限單元網(wǎng)格Fig.4 Finite element mesh

圖5 水管平面布置和水流方向Fig.5 Layout of cooling pipe and water flow direction

混凝土的絕熱溫升曲線可以用以下公式表示

式中θ為絕熱溫升終值，℃，材料的其他熱學(xué)參數(shù)見表1。

表1 其他熱學(xué)參數(shù)Table 1 Other thermal parameters

4.2 計(jì)算結(jié)果分析

1）計(jì)算結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比

測(cè)點(diǎn)1和測(cè)點(diǎn)2位于水管的進(jìn)水口附近。測(cè)點(diǎn)1距離水管8 cm，測(cè)點(diǎn)2距離水管12 cm。受氣溫波動(dòng)的影響，進(jìn)口水溫很難保持穩(wěn)定，實(shí)際施工過程中冷卻水流量有時(shí)候也會(huì)出現(xiàn)小幅波動(dòng)，所以實(shí)測(cè)值的歷時(shí)曲線并不光滑?？傮w來看，計(jì)算值和實(shí)測(cè)值較為吻合。測(cè)點(diǎn)3和測(cè)點(diǎn)4位于水管的出口附近。測(cè)點(diǎn)3距離水管8 cm，測(cè)點(diǎn)4距離水管12 cm。由于流量較小且水管較長(zhǎng)，受沿程水溫（進(jìn)出口水溫差在10℃左右）影響，點(diǎn)3和點(diǎn)4的溫度要明顯高于點(diǎn)1和點(diǎn)2的溫度（見圖6），水溫的改變對(duì)點(diǎn)3和點(diǎn)4的溫度影響也較點(diǎn)1和點(diǎn)2小。故在實(shí)際施工中，定期改變通水的方向是很有必要的。

2）溫度場(chǎng)計(jì)算結(jié)果分析

圖7為沿水流方向水管內(nèi)壁和外壁的平均溫差分布。因進(jìn)口水溫較低，進(jìn)口處的管內(nèi)外壁溫差也大于出口處。

3）兩種不同迭代方法對(duì)比

采用一般的迭代方法，需要迭代15次，沿程水溫才能達(dá)到穩(wěn)定值（誤差小于0.08℃）。而采用筆者改進(jìn)的迭代方法，只需迭代7次即可以達(dá)到穩(wěn)定值。

圖6 計(jì)算值和實(shí)測(cè)值對(duì)比Fig.6 Comparison between calculation value and actual value

圖7 管內(nèi)外壁溫差沿水流方向分布Fig.7 Temperature difference distribution between inner and outer surface along pipe

5 結(jié)論

由于經(jīng)濟(jì)和施工方便等原因，塑料水管越來越廣泛地應(yīng)用于各類大體積混凝土溫控防裂中。在以往含離散水管模型的大體積混凝土溫度場(chǎng)迭代計(jì)算中，塑料水管邊界往往考慮成第三類邊界條件，在這種條件下溫度場(chǎng)計(jì)算的精度與水管管壁放熱系數(shù)有很大的關(guān)系。本文根據(jù)熱平衡原理，推導(dǎo)了水管內(nèi)外壁的溫度關(guān)系。由沿程水溫（沿水流方向水管內(nèi)壁的溫度）、與水管接觸混凝土的熱流量、水管的導(dǎo)熱系數(shù)及管壁厚度，即可推算出水管外壁溫度（水管和混凝土的接觸面），然后在計(jì)算中可將水管外壁的溫度作為第一類邊界條件。水管的導(dǎo)熱系數(shù)及管壁的厚度均可輕易獲取，而沿程水溫、與水管接觸混凝土的熱流量也可以在計(jì)算中獲取。從而避免了將管壁作為第三類邊界條件要通過試驗(yàn)和反演獲取管壁放熱系數(shù)以及可靠性較差的缺點(diǎn)。對(duì)比算例中，如采用傳統(tǒng)算法，模型邊界處的誤差達(dá)到1.67℃，而采用該新算法，誤差僅為0.3℃。工程算例表明：采用這種新的塑料水管邊界條件的處理方法，測(cè)點(diǎn)的溫度能與實(shí)測(cè)值吻合得較好。此外，筆者在原有迭代方法的基礎(chǔ)上還進(jìn)行了改進(jìn)，提高了收斂速度，節(jié)約了計(jì)算時(shí)間；本文的工程算例中，采用一般的迭代方法，需要迭代15次，而采用改進(jìn)的迭代方法，只需迭代7次即可以達(dá)到穩(wěn)定值。

[1] 強(qiáng)晟，鄭偉忠，張勇強(qiáng)，等. 基于改進(jìn)微粒群法和有限元法的混凝土溫控方案優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，30(16)：75－83. Qiang Sheng, Zheng Weizhong, Zhang Yongqiang, et al. Optimization of concrete temperature control measures based on improved particle swarm optimization and finite element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(16): 75－83. (in Chinese with English abstract)

[2] Zhu Zhenyang, Qiang Sheng, Chen Weimin. A model for temperature influence on concrete hydration exothermic rate (part one: theory and experiment)[J]. Journal of Wuhan University of Technology (Materials Science Edition), 2014, 29(3): 540－545.

[3] Wang Haibo, Qiang Sheng, Zhu Zhenyang. The whole process simulation and research on temperature field and thermal stress for large sluice pier[J]. Advanced Materials Research, 2011(163/164/165/166/167): 259－263.

[4] Zhu Zhenyang, Qiang Sheng, Liu Minzhi. Cracking mechanism of long concrete bedding cushion and prevention method[J]. Advanced Materials Research, 2011(163/164/165/166/167): 880－887.

[5] Chen Shenhong, Su Peifang, et al. Composite element algorithm for the thermal analysis of mass concrete simulation of cooling pipes[J]. International Journal of Numerical Methods for Heat & Fluid Flow, 2011, 21(3/4): 434－447.

[6] 蘇培芳，陳勝宏，田甜. 混凝土水管冷卻的復(fù)合單元算法[J].武漢理工大學(xué)學(xué)報(bào)，2010，32(24)：48－53. Su Peifang, Chen Shenghong, Tian Tian. Preliminary research on composite element algorithm for the concrete containing cooling pipes[J]. Journal of Wuhan University of Technology, 2010, 32(24): 48－53. (in Chinese with English abstract)

[7] 朱伯芳. 論混凝土壩的水管冷卻[J]. 水利學(xué)報(bào)，2010，41(5)：505－513. Zhu Bofang. On pipe cooling of concrete dams[J]. Journal of Hydraulic Engineering, 2010, 41(5): 505－513. (in Chinese with English abstract)

[8] 王振紅，張國(guó)新，劉毅，等. 混凝土水管冷卻試驗(yàn)與溫控參數(shù)的反分析[J]. 四川大學(xué)學(xué)報(bào)：工程科學(xué)版，2011，43(3)：56－60. Wang Zhenhong, Zhang Guoxin, Liu Yi, et al. Test and inverse analysis for temperature control parameters of concrete with cooling pipe[J]. Journal of Sichuan University (Engineering Science Edition), 2011, 43(3): 56－60. (in Chinese with English abstract)

[9] 劉寧，劉光廷. 水管冷卻效應(yīng)的有限元子結(jié)構(gòu)模擬技術(shù)[J].水利學(xué)報(bào)，1997(12)：43－49. Liu Ning, Liu Guangting. Sub-structural FEM for the thermal effect of cooling pipes in mass concrete structures[J]. Journal of Hydraulic Engineering, 1997(12): 43－49. (in Chinese with English abstract)

[10] 頡志強(qiáng)，強(qiáng)晟，許樸，等. 水管冷卻混凝土溫度場(chǎng)和應(yīng)力場(chǎng)計(jì)算的有限元子結(jié)構(gòu)法[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2011，27(5)：13－18. Xie Zhiqiang, Qiang Sheng, Xu Pu, et al. Finite element substructure method for calculation of pipe cooling concrete thermal field and stress field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(5): 13－18. (in Chinese with English abstract)

[11] 陳國(guó)榮，許文濤，楊昀，等. 含冷卻水管大體積混凝土溫度場(chǎng)計(jì)算的一種新方法[J]. 計(jì)算物理，2012，29(3)：411－416. Chen Guorong, Xu Wentao, Yang Yun, et al. Computation method for temperature field of mass concrete containing cooling water pipes[J]. Chinese Journal of Computaitonal Physics, 2012, 29(3): 411－416. (in Chinese with English abstract)

[12] 劉曉青，李同春，韓勃. 模擬混凝土水管冷卻效應(yīng)的直接算法[J]. 水利學(xué)報(bào)，2009，40(7)：892－896. Liu Xiaoqing, Li Tongchun, Han Bo. Direct algorithm for simulating cooling effect of water pipes in concrete[J]. Journal of Hydraulic Engineering, 2009, 40(7): 892－896. (in Chinese with English abstract)

[13] Zuo Zheng, Hu Yu, Li Qingbin, et al. An extended finite element method for pipe-embedded plane thermal analysis[J]. Finite Elements in Analysis and Design, 2015(102/103): 52－64.

[14] 強(qiáng)晟，朱岳明，丁兵勇，等. 基于冷卻水管離散算法的重力壩溫控研究[J]. 水利能源科學(xué)，2008，26(5)：93－95. Qiang Sheng, Zhu Yueming, Ding Bingyong, et al. Study on temperature control of gravity dam based on explicit algorithm for cooling pipe[J]. Water Resources and Power, 2008, 26(5): 93－95. (in Chinese with English abstract)

[15] 朱岳明，徐之青，賀金仁，等. 混凝土水管冷卻溫度場(chǎng)的計(jì)算方法[J]. 長(zhǎng)江科學(xué)院學(xué)報(bào)，2003，20(2)：19－22. Zhu Yueming, Xu Zhiqing, He Jinren, et al. A calculation method for solving temperature field of mass concrete with cooling pipes[J]. Journal of Yangtze River Scientific Research Institute, 2003, 20(2): 19－22. (in Chinese with English abstract)

[16] 許樸，朱岳明，賁能慧. 倒T型混凝土薄壁結(jié)構(gòu)施工期溫度裂縫控制研究[J]. 水利學(xué)報(bào)，2009，40(8)：969－975. Xu Pu, Zhu Yueming, Ben Nenghui. Study on thermal cracking control of inverted T-shaped concrete structures during construction[J]. Journal of Hydraulic Engineering, 2009, 40(8): 969－975. (in Chinese with English abstract)

[17] 王振紅，朱岳明，于書萍，等. 水閘閘墩施工期溫度場(chǎng)和應(yīng)力場(chǎng)仿真計(jì)算分析[J]. 天津大學(xué)學(xué)報(bào)，2008，41(4)：476－481. Wang Zhenhong, Zhu Yueming, Yu Shuping, et al. Simulation and analysis of temperature field and stress field of sluice pier concrete during construction[J]. Journal of Tianjin University, 2008, 41(4): 476－481. (in Chinese with English abstract)

[18] 朱振泱，強(qiáng)晟，鄭占強(qiáng)，等. 用遺傳算法確定考慮溫度歷程的混凝土水化放熱模型參數(shù)及試驗(yàn)驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(1)：86－92. Zhu Zhenyang, Qiang Sheng, Zheng Zhanqiang, et al. Determination of parameters for hydration exothermic model considering concrete temperature duration by genetic algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(1): 86－92. (in Chinese with English abstract)

[19] Hattel J H, Thorborg J. A numerical model for predicting the thermal mechanical conditions during hydration of early age concrete[J]. Applied Mathematical Modeling, 2003, 27(1): 1－26.

[20] Zhu Zhenyang, Qiang Sheng, Chen Weimin. A new method solving the temperature field of concrete around cooling pipes[J]. Comput. Concrete, 2013, 11(5): 441－462.

[21] 朱振泱，強(qiáng)晟，陳煒旻，等. 含水管混凝土溫度場(chǎng)的改進(jìn)離散迭代算法[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2014，22(3)：492－500. Zhu Zhenyang, Qiang Sheng, Chen Weimin, et al. Improvement on iteration method solving temperature field of mass concrete with pipes[J]. Journal of Basic Science and Engineering, 2014, 22(3): 492－500. (in Chinese with English abstract)

[22] 張軍，段亞輝. 混凝土冷卻水管的有限元沿程水溫改進(jìn)算法[J]. 華中科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2014(2)：56－58. Zhang Jun, Duan Yahui. FEM analysis of concrete with cooling pipes using an improved method on the calculation of water temperature along the pipes[J]. J. Huazhong Univ. of Sci. & Tech: Natural Science Edition, 2014(2): 56－58. (in Chinese with English abstract)

Algorithm to simulate concrete temperature control cooling pipe boundary based on heat flux integration

Zhu Zhenyang1, Liu Minzhi2, Qiang Sheng3, Xiang Jianfang1
(1. State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China; 2. Yinzhou District Quality and Safety Supervision Station of Water Conservancy Projects, Ningbo 315100, China; 3. College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China)

When using the explicit iterative method to solve the temperature of mass concrete with cooling pipes, it is generally considered that the inner and outer surface of metal pipes can be neglected but the temperature difference cannot be neglected when using the plastic pipes. And the plastic pipes are usually regarded as the third boundary condition. For the past researchers, the coefficient of this kind of boundary condition can be got by experiment and inversion, which is yet expensive and may also not be reliable sometimes. To solve the problem, on the base of heat balance condition, a new calculation method is brought forward. It is well known that concrete is a poor conductor of heat, and there is a large temperature difference between the concrete and the cooling water. So, in the shell of a cooling pipe and the concrete near it, it can be assumed that the heat flux is only discharged by cooling water in the pipe and the direction of the temperature gradient is perpendicular to the cooling pipe surface. So, the heat fluxes passing through any circle (take the center pipe as the center of those circles) in the shell of the plastic pipe are equal. Based on these basic principles, the temperature of cooling pipe outer surface can be obtained by the heat flux of the concrete around the pipe, the thermal conductivity, the thickness of the pipe and the temperature of pipe inner surface. When using the conventional iterative method to solve the temperature of the mass concrete with cooling pipes, the iterative method should be used for the unknown water temperature distributions along the cooling pipes. For this new method, the temperature distributions along the inner surface of the pipes is also unknown, so the iterative method should be also used. With this new method, when using the conventional iterative method, the convergence speed is relatively low, or even can not converge. To solve this problem, the iterative algorithm is also improved. When the iteration time is (N-1) and N separately, it is assumed that the corresponding calculated temperature on the outer surface of the cooling pipe is Tn-1and Tnrespectively. And then, when using 0.5(Tn-1+ Tn) as the initial calculation condition for the (N+1)thtime, the convergence of the iteration can be easily achieved. The convergence condition of the improved method was proved by mathematical deduction, and the deduction results showed that the convergence could be always reached in different engineering cases. A comparing numerical example was used to comparing the accuracy of the new method and the conventional explicit iterative method. In this comparing numerical example, the calculation results of the finite element method (FEM) considering the pipe as a part of mesh were considered as the theoretical solution. The calculation results showed that in the concrete near adiabatic boundary of the mesh in the comparing example, the temperature difference between the calculation result of conventional explicit iterative method and the theoretical solution was 1.67℃ , and the temperature difference between the calculation result of improved method and the theoretical solution was only 0.3℃ . So, the improved method can be more accurate than the conventional explicit iterative algorithm. Using these new achievements, the temperature field of a concrete block during construction period was simulated, and the calculation results and testing results were compared. The total number of the iteration times was 15 for the conventional iterative method and only 7 for the improved method in this engineering example. The results show that the calculation value is close to the actual value, and this algorithm has high convergence speed. So this method can be used in engineering projects to prevent mass concrete from cracking.

temperature; models; concrete; iterative algorithm

10.11975/j.issn.1002-6819.2016.09.012

TV315

A

1002-6819(2016)-09-0083-07

朱振泱，劉敏芝，強(qiáng) 晟，相建方. 基于熱流量積分的混凝土溫控水管冷卻邊界模擬算法[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(9)：83－89.

10.11975/j.issn.1002-6819.2016.09.012 http://www.tcsae.org

Zhu Zhenyang, Liu Minzhi, Qiang Sheng, Xiang Jianfang. Algorithm to simulate concrete temperature control cooling pipe boundary based on heat flux integration[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(9): 83－89. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2016.09.012 http://www.tcsae.org

2015-12-22

2016-03-14

國(guó)家自然科學(xué)基金資助項(xiàng)目（51409264，51509020，51209219）

朱振泱，男，福建三明人，中國(guó)水利水電科學(xué)研究院高級(jí)工程師，主要從事混凝土溫控防裂研究和大體積混凝土溫度裂縫擴(kuò)展研究。北京 中國(guó)水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室，100038。Email：1219921552@qq.com