杜 瑋，周超峰，鄧 旭，馬躍先

(1. 河南省農村水電及電氣化發(fā)展中心，鄭州 450003； 2. 河南省陸渾水庫管理局，河南 洛陽 471412；3. 河南鄭大水利科技有限公司，鄭州 450001； 4. 鄭州大學水利科學與工程學院，鄭州 450002)

0 引 言

大體積混凝土是現代工程中的重要結構形式，合理計算并確定混凝土的溫度場及溫度應力分布，是其溫控施工控制的重要環(huán)節(jié)[1,2]。目前，大體積混凝土溫度場及溫度應力的計算包括經驗公式法、差分算法以及有限單元法等[3,4]，但上述計算均需要依靠混凝土的溫度參數進行計算，而溫度參數的獲取主要通過經驗計算或試驗法獲得，由于現場施工環(huán)境、測量設備、施工工藝的不同，導致獲取的溫度度參數與實際工程仍然存在一定的誤差。

反演理論利用實測數據對計算參數進行修正，但針對混凝土溫度參數的反演研究相對較少，筆者曾提出混凝土溫度場、濕度場的差分計算模型[5,6]，并針對混凝土濕度場的濕度擴散系數提出差分反演模型[7]，但該模型的反演參數較少，模型較為簡單，而混凝土溫度參數較多，反演模型復雜，濕度場的反演模型和求解方法無法直接應用。

本文基于差分理論，構建混凝土溫度場的差分計算模型，利用實測數據建立混凝土溫度參數的差分反演模型，并針對反演模型，設計相應的遺傳算法模型，通過對混凝土溫度參數的反演計算，可以較為準確地獲取符合實際工程的溫度場參數，從而對大體積混凝土的溫控提供指導。

1 混凝土溫度場的差分計算模型

1.1 熱傳導方程

大體積混凝土的熱傳導方程可表達為[1]:

(1)

式中：c為混凝土比熱，kJ/(kg·℃)；λ為混凝土的導熱系數，kJ/(m·d·℃)；τ為時間,d；ρ為混凝土密度，kg/m3；θ為混凝土絕熱溫升；x,y,z為混凝土內部點的位置坐標；T為混凝土內部坐標點溫度，℃。

1.2 初始條件和邊界條件

對大體積混凝土的溫度場進行求解時，需要確定其初始條件和邊界條件。

初始條件：初始溫度為其初始瞬時溫度，即澆筑溫度，為常數。

T(x,y,z,0)=T0=const

(2)

邊界條件：根據混凝土的接觸介質不同，其邊界條件可以分為第一類邊界條件～第四類邊界條件，具體可參見文獻[1，5]；其中，第三類邊界條件使用較多，其假定混凝土表面的熱流量與其表面溫度T和氣溫Tf之差成正比。即：

(3)

式中：β為混凝土表面的放熱系數。

1.3 差分計算模型

對于上述構建的熱傳導方程，劃分差分單元，建立其差分計算模型，其中計算模型,如圖1所示，其中x方向格距為h，y方向格距為l，z方向格距為v，用中心差商代替微商。

溫度場的差分計算模型可表示為：

(4)

(5)

(6)

代入式(1)得：

T0,τ+Δτ=(1-2r1-2r2-2r3)T0,τ+r1(T1,τ+T2,τ)+

r2(T3,τ+T4,τ)+r3(T5,τ+T6,τ)+Δθ

(7)

式中：r1=aΔτ/h2，r2=aΔτ/l2，r3=aΔτ/v2，Δθ=θ(τ+Δτ)-θ(τ)，當r1+r2+r3<0.5時，計算結果穩(wěn)定。

第三類邊界條件的差分格式可利用牛頓后插公式，建立其差分格式為：

(8)

式中：T邊為混凝土邊界處溫度；T內1為混凝土內部距離邊界為1單位的節(jié)點的溫度值；T內2混凝土內部距離邊界為2單位的節(jié)點的溫度值，對于x方向，取s1=2βh/λ，對于y方向，取s2=2βl/λ，對于z方向，取s3=2βv/λ；Ta為外界介質溫度。

2 溫度場實例分析

2.1 計算參數

針對某混凝土工程，為準確獲取混凝土的溫度參數，在實驗室進行混凝土溫控試驗，選擇50 cm×50 cm×50 cm的混凝土試塊，在恒溫恒濕養(yǎng)護箱內保持恒溫(15 ℃)、恒濕(95%濕度)的工況下養(yǎng)護20 d，其中混凝土的各溫度參數通過經驗估算見表1。

表1 混凝土溫度參數Tab.1 Temperature parameters of concrete

混凝土的絕熱溫升可以表示為：

Q(τ)=Q0(1-e-aτb)

(9)

式中：Q(τ)為齡期為τ時的累積水化熱，kJ/kg；Q0為最終水化熱，kJ/kg；τ為齡期；a,b為固定系數。

其中，混凝土的絕熱溫升可根據下式計算：

(10)

式中：W為每方混凝土中的水泥含量，kg/m3；F為每方混凝土中摻合料含量，kg/m3；k為折減系數，對于粉煤灰，可取0.25。

根據經驗進行計算取值見表2。

表2 混凝土最終水化熱及系數取值Tab.2 Final hydration heat and coefficient of concrete

為對混凝土試塊的溫度變化進行測量，設置溫度傳感器進行測量，所述溫度傳感器布置位置見圖2。

2.2 溫度場計算

根據本文構建的差分模型，對上述試件進行溫度場差分計算，其中混凝土單元劃分為5 cm，時間劃分為0.002 d(即每隔2.88 min計算一次)，計算周期為20 d，分別計算各測點的混凝土溫度變化歷時曲線見圖3。

通過圖3可知，混凝土各測點的溫度先升高后降低，前期溫升主要是由于水泥水化熱產生的熱量大于空氣散熱，后期水泥水化熱產生的熱量小于空氣散熱，則溫度出現降低，并逐漸接近外界空氣氣溫。對于各測點而言，越接近空氣表面的測點，其受外界溫度影響越劇烈，其溫升越小。

為展現混凝土各點的溫度分布情況，分別計算2、20 d后，混凝土的溫度場分布云圖，計算結果見圖4、5。

由圖4、5分析可知，混凝土的中心點溫度最高，邊界點的溫度最低，且混凝土試件的四邊角溫度最低，主要是由于其同時受到兩個空氣擴散面的影響，其溫度受外界空氣溫度最為劇烈。

2.3 溫度誤差分析

通過采集各測點的溫度數據，每隔0.5 d提取一次數據，共提取20 d的40個溫度數據，通過經驗估算得到的溫度計算數據與實測數據進行對比，其誤差分析見圖6、7。

通過圖6、7分析可知，通過經驗估算得到的溫度計算數據與實測溫度數據存在較大的誤差，尤其前8 d的數據誤差較大，最大誤差接近1 ℃，而后期誤差相對較小，主要是前期溫升階段，受混凝土溫度參數的影響較大，如不合理確定溫度參數，會導致理論計算與實際產生較大的偏差。

3 基于協(xié)同遺傳進化的反演模型及求解

3.1 反演模型建立

為消除理論計算與實測數據的誤差，基于實測數據，建立反演模型，考慮誤差的正負變化，采用誤差值的累積平方和作為反演目標函數值，其反演模型可以表達為：

(11)

式中：W為計算值與實測值的誤差平方和；T為實測溫度值；T′為理論計算值；m為測點個數；n為各測點的溫度測量/計算的時段數；f為溫度計算的差分函數；a1,a2,…,as為各溫度參數，基于本例，可以選擇為混凝土的導熱系數、比熱系數、密度、表面放熱系數，絕熱溫升、經驗指數a、經驗指數b。

3.2 GA算法設計

上述反演模型為非線性多參數的復雜求解模型，無法利用理論計算求得其解析解，本文擬采用遺傳仿生算法進行求解，設計其算法流程見圖8。

其中，遺傳參數種群包括混凝土各溫度參數，適應度函數即為反演模型的 值，遺傳規(guī)則可制定為：適應度函數值小的個體向下一代遺傳進行，其中適應度函數值越小，其遺傳概率更大；變異規(guī)則可制定為：每代種群中的變異個數取種群個數的10%，以避免遺傳種群的局部收斂；進化代數作為遺傳終止判別條件，終止后，對種群的所有歷史解進行對比，尋找最小值，并輸出該最小值對應的各種群值。

3.3 遺傳反演求解

按照本文設計的遺傳算法，進行反演模型求解，其中各溫度參數種群規(guī)模設定為100，遺傳代數為40代，其中各溫度場參數的求解范圍可確定見表3。

表3 溫度參數的求解范圍表Tab.3 Solution range table of temperature parameters

以導熱系數和比熱系數為例，種群規(guī)模選擇確定其初始種群分布見圖9、10。

遺傳進化40代后的種群分布見圖11、12。

遺傳進化的適應度函數跟蹤見圖13。

4 反演結果及評價

4.1 反演結果

根據上述遺傳算法進行反演求解，40代后，適應度函數最小值為1.236，確定最優(yōu)種群的溫度參數見表4。

表4 反演后的溫度參數表Tab.4 Temperature parameter table after inversion

4.2 反演評價

利用反演后得到的參數進行溫度場差分計算，得到其各測點的溫度變化與實測值的歷時過程線及誤差分布見圖14、15。

其中分別對反演前后2 d及20 d的溫度誤差場的分布對比見圖16、17。

各測點的溫度誤差平方和比較見表5。

通過上述圖表可知，經過反演后的溫度參數經過差分計算，其計算溫度值與實際溫度值吻合度較好，各測點的溫度誤差平方和均下降75%以上，較大的消除了經驗參數帶來的誤差。

表5 各測點的溫度誤差平方和比較表Tab.5 Comparison table of sum of squared temperature error of measuring point

5 結 語

(1)通過溫度場差分模型，可以對混凝土溫度場進行仿真模擬，但溫度參數對溫度場的影響較大，僅依靠經驗確定溫度參數，誤差較大；

(2)利用實測溫度數據，可以建立混凝土溫度參數的反演計算模型，該模型可以較好地修正經驗參數，減少誤差；

(3)GA算法可以較好地求解非線性多參數的優(yōu)化問題，通過合理的適應度函數選擇、遺傳規(guī)則及變異規(guī)則制定，可以較好地求解混凝土溫度參數的差分-反演模型。

□