陳明華，靳 良

(中交四航局第三工程有限公司，廣東 湛江 524005)

在水利水電工程建設(shè)中，水化熱導(dǎo)致的溫度變化不但可能導(dǎo)致裂縫的產(chǎn)生，對(duì)構(gòu)筑物的應(yīng)力狀態(tài)也有不可忽略的影響，有時(shí)溫度應(yīng)力在數(shù)值上可能超過其他荷載引起的應(yīng)力[1]。因此，在施工前進(jìn)行有限元計(jì)算，并針對(duì)開裂風(fēng)險(xiǎn)確定施工中裂縫控制措施，對(duì)控制裂縫的產(chǎn)生、保證工程質(zhì)量具有重要意義。曹周紅等[2]采用有限元分析，得到水化熱產(chǎn)生的溫度應(yīng)力導(dǎo)致船閘底板開裂的結(jié)論；徐可等[3]通過對(duì)船閘底板建立三維有限元模型進(jìn)行分析，認(rèn)為水化熱溫度場(chǎng)時(shí)效模式更能準(zhǔn)確地模擬工程實(shí)際；黃澤欽等[4]采用有限元仿真計(jì)算與實(shí)測(cè)數(shù)據(jù)對(duì)比的方法討論分析了水化熱溫度變化規(guī)律。但是在施工前對(duì)水化熱及其產(chǎn)生的應(yīng)力進(jìn)行有限元計(jì)算很難模擬準(zhǔn)確，其困難主要存在于軟件中的結(jié)構(gòu)放熱函數(shù)及結(jié)構(gòu)的強(qiáng)度發(fā)展函數(shù)的設(shè)置?？赏ㄟ^將軟件計(jì)算結(jié)果與試驗(yàn)節(jié)段實(shí)測(cè)結(jié)果進(jìn)行對(duì)比，通過調(diào)整這2個(gè)函數(shù)的設(shè)置使其與實(shí)測(cè)結(jié)果接近，從而適應(yīng)下一節(jié)段的施工分析。

1 工程概況

貴港二線船閘閘首、閘室及導(dǎo)航墻均屬于大體積混凝土。上閘首結(jié)構(gòu)單體方量為13 660 m3，下閘首結(jié)構(gòu)單體方量為18 094 m3，閘室單體方量為8 528 m3，產(chǎn)生水化熱較大，易受溫度收縮應(yīng)力和外部約束應(yīng)力等影響而生成裂縫。水泥使用魚峰牌中熱硅酸鹽水泥P·MH 42.5，水化熱較低，對(duì)混凝土的控裂有利，水泥的各項(xiàng)指標(biāo)如表1所示。

表1 水泥指標(biāo)

2 溫度監(jiān)測(cè)方法及結(jié)果

2.1 溫度監(jiān)測(cè)方法

2.1.1測(cè)溫計(jì)的布設(shè)

為探索主體混凝土澆筑后的水化熱溫升規(guī)律，進(jìn)而有針對(duì)性地采取合理的溫控措施，在船閘主體結(jié)構(gòu)中預(yù)埋溫度計(jì)進(jìn)行觀測(cè)。溫度計(jì)的布設(shè)是根據(jù)混凝土分層分塊澆筑圖，在混凝土分層中體積相對(duì)較大的典型倉面來布置。其中，上、下游導(dǎo)航墻各布設(shè)12套溫度計(jì)，閘室墻布設(shè)32套溫度計(jì)，上游導(dǎo)航墻A段的溫度計(jì)布設(shè)如圖1所示。

圖1 上游導(dǎo)航墻A段溫度計(jì)分布(高程：m； 尺寸：mm)

2.1.2溫度數(shù)據(jù)的記錄

混凝土澆筑完成后立即開始按設(shè)計(jì)規(guī)定時(shí)間間隔進(jìn)行溫度監(jiān)測(cè)及數(shù)據(jù)采集，混凝土單點(diǎn)溫度及混凝土內(nèi)外溫差達(dá)到峰值前每2 h監(jiān)測(cè)1次，峰值過后每4 h監(jiān)測(cè)1次。

2.2 溫度監(jiān)控結(jié)果

圖2給出導(dǎo)航墻在澆筑完成后的溫度變化時(shí)程曲線。構(gòu)筑物中心部位的溫度在72 h內(nèi)增長(zhǎng)了25 ℃，增長(zhǎng)速率為0.347 ℃h，下降到初始溫度用了384 h，下降速率為0.065 ℃h。

圖2 導(dǎo)航墻實(shí)測(cè)溫度

3 有限元仿真溫度計(jì)算結(jié)果及分析

從圖2可以看出，船閘內(nèi)部混凝土的溫度變化為高度非線性，不同體積的構(gòu)筑物、同一構(gòu)筑物不同部位的溫度變化規(guī)律都不一樣，實(shí)測(cè)構(gòu)筑物某一點(diǎn)的溫度變化規(guī)律無法推演到同一構(gòu)筑物的其他部位，更很難應(yīng)用到其他工程。若布置過多的溫度監(jiān)測(cè)點(diǎn)，通過統(tǒng)計(jì)來分析水化熱溫度的規(guī)律，也存在著成本過高、對(duì)工程影響較大等問題，因此采用有限元仿真計(jì)算來模擬分析大體積混凝土水化熱溫度規(guī)律從而指導(dǎo)工程施工很有必要。

3.1 有限元計(jì)算方法

采用有限元軟件Midas Civil對(duì)分階段澆筑的導(dǎo)航墻進(jìn)行水化熱分析。由于導(dǎo)航墻軸向具有對(duì)稱性，選取12進(jìn)行建模，如圖3所示。其中CS1為第1個(gè)施工階段(28.5～31.0 m層)，間隔7 d后進(jìn)行其上CS2為第2個(gè)施工階段(31.0～34.0 m層)。

圖3 導(dǎo)航墻計(jì)算模型

3.2 混凝土溫度計(jì)算結(jié)果

圖4為導(dǎo)航墻應(yīng)力計(jì)算實(shí)體單元。選取構(gòu)筑物表面4個(gè)節(jié)點(diǎn)進(jìn)行溫度應(yīng)力計(jì)算結(jié)果的分析，圖中每個(gè)單元格的長(zhǎng)寬均為0.5 m，選取分析的4個(gè)節(jié)點(diǎn)分布在距結(jié)構(gòu)物表面2 m內(nèi)，節(jié)點(diǎn)號(hào)從外到里分別為13 971、12 998、13 327、13 787。

圖4 導(dǎo)航墻應(yīng)力計(jì)算實(shí)體單元

3.3 溫度計(jì)算結(jié)果及適應(yīng)性調(diào)整

3.3.1溫度計(jì)算結(jié)果分析

圖5a)為施工前對(duì)導(dǎo)航墻進(jìn)行仿真計(jì)算的溫度變化時(shí)程曲線。計(jì)算結(jié)果顯示，混凝土在水泥水化過程中的最高溫度為63 ℃，處于結(jié)構(gòu)的內(nèi)部，約發(fā)生在澆筑完成后的48 h，小于實(shí)測(cè)的溫度峰值出現(xiàn)時(shí)間72 h。因此須減小放熱函數(shù)的斜率來達(dá)到推遲計(jì)算溫度峰值出現(xiàn)時(shí)間的目的，使得水化熱溫度分析結(jié)果接近實(shí)測(cè)結(jié)果，只有這樣，有限元分析的溫度應(yīng)力計(jì)算結(jié)果才具有參考意義。

圖5 導(dǎo)航墻溫度計(jì)算結(jié)果

3.3.2放熱函數(shù)調(diào)整

圖5b)為施工第1階段后，根據(jù)導(dǎo)航墻實(shí)測(cè)的溫度變化曲線調(diào)整放熱函數(shù)得到的導(dǎo)航墻溫度計(jì)算結(jié)果。從圖5b)可知，混凝土在水泥水化過程中的最高溫度為63 ℃，處于結(jié)構(gòu)的內(nèi)部，約發(fā)生在澆筑完成后的72 h，與實(shí)測(cè)值更為接近。

3.3.3強(qiáng)度發(fā)展函數(shù)調(diào)整

圖6為施工前根據(jù)規(guī)范[5]公式設(shè)置的強(qiáng)度發(fā)展函數(shù)，因?yàn)樵谑┕で盁o試驗(yàn)數(shù)據(jù)，相關(guān)參數(shù)選用規(guī)范建議值。施工前的應(yīng)力計(jì)算結(jié)果顯示，在澆筑完成后導(dǎo)航墻的拉應(yīng)力超過了容許抗拉強(qiáng)度，但經(jīng)拆模后觀察，導(dǎo)航墻的表面并未產(chǎn)生裂縫。除了放熱函數(shù)的設(shè)置偏離實(shí)際導(dǎo)致放熱過于迅猛的原因外，混凝土的強(qiáng)度發(fā)展函數(shù)設(shè)置過于平緩也可能是另外的原因，因此須對(duì)強(qiáng)度發(fā)展函數(shù)作進(jìn)一步調(diào)整。

圖6 (施工前分析)軟件設(shè)置強(qiáng)度發(fā)展函數(shù)

4 應(yīng)力計(jì)算結(jié)果及分析

4.1 結(jié)構(gòu)應(yīng)力計(jì)算結(jié)果

在施工第1階段后，根據(jù)導(dǎo)航墻實(shí)測(cè)溫度變化曲線對(duì)放熱函數(shù)和強(qiáng)度發(fā)展函數(shù)進(jìn)行調(diào)整后，計(jì)算得到的導(dǎo)航墻表面點(diǎn)溫度應(yīng)力如圖7所示。

圖7 導(dǎo)航墻表面點(diǎn)溫度應(yīng)力

4.2 結(jié)構(gòu)計(jì)算應(yīng)力

從圖7可以看出，混凝土的拉應(yīng)力與澆筑時(shí)間為高度非線性關(guān)系，但都呈現(xiàn)出先受壓再受拉的規(guī)律，混凝土表面處于受壓狀態(tài)的時(shí)間約為16 h，距表面0.5 m處的分析點(diǎn)處于受壓狀態(tài)的時(shí)間約為120 h，呈現(xiàn)出越接近結(jié)構(gòu)物中心的部位，其處于受壓狀態(tài)的時(shí)間越長(zhǎng)的規(guī)律。證明結(jié)構(gòu)在澆筑初期因溫度應(yīng)力而產(chǎn)生的張拉裂縫的深度很可能只在0～0.5 m的范圍。

4.3 結(jié)構(gòu)溫度裂縫產(chǎn)生的可能性

圖7給出結(jié)構(gòu)某點(diǎn)的應(yīng)力及容許抗拉強(qiáng)度的時(shí)程變化曲線。通過觀察某點(diǎn)應(yīng)力與容許抗拉強(qiáng)度曲線的接近程度來定性分析溫度裂縫產(chǎn)生的可能性，但是無法定量進(jìn)行分析。水化熱參數(shù)化分析指出，用溫度裂縫指數(shù)i值來預(yù)測(cè)是否發(fā)生裂縫[5]，如式(1)：

(1)

式中:i為裂縫指數(shù)，ft為混凝土抗拉強(qiáng)度，σt為發(fā)生的溫度應(yīng)力。

根據(jù)結(jié)構(gòu)對(duì)裂縫控制要求的不同對(duì)裂縫指數(shù)有不同的要求。韓國(guó)規(guī)范規(guī)定，當(dāng)需要防治裂縫發(fā)生時(shí)，裂縫指數(shù)應(yīng)在1.5以上；當(dāng)需要限制裂縫發(fā)生時(shí)，裂縫指數(shù)應(yīng)在1.2～1.5的范圍；當(dāng)需要限制有害裂縫發(fā)生時(shí)，裂縫指數(shù)應(yīng)在0.7～1.2的范圍。圖8為韓國(guó)規(guī)范中的溫度裂縫指數(shù)與裂縫發(fā)生幾率關(guān)系圖，該圖為韓國(guó)的實(shí)際工程統(tǒng)計(jì)結(jié)果。

圖8 溫度裂縫指數(shù)與裂縫發(fā)生幾率關(guān)系

4.3.1施工前的有限元分析結(jié)果

圖9a)為施工前有限元分析的不同結(jié)構(gòu)部位點(diǎn)裂縫指數(shù)隨澆筑完成時(shí)間的變化關(guān)系。由圖7得知，裂縫指數(shù)大于2時(shí)，溫度裂縫發(fā)生幾率幾乎為零，因此圖8略去了混凝土水化過程中裂縫指數(shù)大于3的數(shù)據(jù)。

從圖9b)可知，構(gòu)筑物內(nèi)部的裂縫指數(shù)幾乎在1.4以上，裂縫發(fā)生的幾率在10%以下。而構(gòu)筑物表面的裂縫指數(shù)在澆筑完成后48～96 h小于1.0，其裂縫發(fā)生的幾率在60%以上，但實(shí)際上，在導(dǎo)航墻相應(yīng)位置的表面并未發(fā)現(xiàn)裂縫。這是因?yàn)樵谑┕で叭狈?shí)測(cè)數(shù)據(jù)，計(jì)算參數(shù)采用規(guī)范推薦值，放熱函數(shù)及強(qiáng)度發(fā)展函數(shù)的相關(guān)設(shè)置偏離了實(shí)際情況，使得計(jì)算結(jié)果與實(shí)際有所偏離。

圖9 不同結(jié)構(gòu)部位點(diǎn)的裂縫指數(shù)

4.3.2施工第1階段后深化計(jì)算結(jié)果

在施工第1階段后，根據(jù)導(dǎo)航墻實(shí)測(cè)的溫度變化曲線對(duì)放熱函數(shù)和強(qiáng)度發(fā)展函數(shù)進(jìn)行調(diào)整后，計(jì)算得到的裂縫指數(shù)如圖9b)所示?？梢钥闯?，結(jié)構(gòu)的最小裂縫指數(shù)在結(jié)構(gòu)的外側(cè)，約發(fā)生在澆筑完成后36 h，imin=1.8(> 1.4)，滿足規(guī)范要求，分析結(jié)果也和結(jié)構(gòu)表面不存在裂縫的工程實(shí)際相符。

4.4 對(duì)相關(guān)規(guī)范的分析

4.4.1我國(guó)規(guī)范與韓國(guó)規(guī)范的對(duì)比

我國(guó)有關(guān)規(guī)范定義了抗裂安全系數(shù)Kf，Kf=εpEcσ，其中εpEc為混凝土的抗拉強(qiáng)度，σ為各種溫差所產(chǎn)生的溫度應(yīng)力之和。因此,韓國(guó)規(guī)范的裂縫指數(shù)i與我國(guó)規(guī)范的抗裂安全系數(shù)Kf在本質(zhì)上是一致的?！端\(yùn)工程大體積混凝土溫度裂縫控制技術(shù)規(guī)程》(簡(jiǎn)稱《水運(yùn)》)中規(guī)定K≥1.4[6]，《塊體基礎(chǔ)大體積混凝土施工技術(shù)規(guī)程》(簡(jiǎn)稱《塊體》)中規(guī)定K=1.15[7]，《混凝土拱壩設(shè)計(jì)規(guī)范》(簡(jiǎn)稱《拱壩》)中規(guī)定Kf取值為1.5～1.8[8]，《混凝土重力壩設(shè)計(jì)規(guī)范》(簡(jiǎn)稱《重力壩》)中規(guī)定Kf取值為1.5～2.0[9]，韓國(guó)規(guī)范中“需要防治裂縫發(fā)生時(shí)”所要求的裂縫指數(shù)i與我國(guó)相關(guān)規(guī)范規(guī)定相一致。通過圖8可知，抗裂安全系數(shù)Kf(溫度裂縫指數(shù))在1.5以上時(shí)，裂縫發(fā)生概率在5%以下，亦即有95%的保證率不發(fā)生溫度裂縫。

4.4.2我國(guó)不同行業(yè)規(guī)范之間的對(duì)比

《水運(yùn)》參考了《塊體》的溫度應(yīng)力計(jì)算公式，分為自約束拉應(yīng)力(表面拉應(yīng)力)與外約束拉應(yīng)力(內(nèi)部拉應(yīng)力)2種情況計(jì)算。而《重力壩》和《拱壩》一樣把溫差分為2部分：水化熱溫升與均勻溫差，其中水化熱溫升引起的溫度應(yīng)力計(jì)算采用的是基礎(chǔ)影響系數(shù)A，均勻溫差引起的溫度應(yīng)力計(jì)算采用的是基礎(chǔ)約束系數(shù)R。

《水運(yùn)》與《塊體》的外約束拉應(yīng)力計(jì)算中采用R值，因?yàn)镽>A，所以采用約束系數(shù)R進(jìn)行簡(jiǎn)化計(jì)算得到的溫度應(yīng)力也較大，因此《水運(yùn)》與《塊體》規(guī)定的抗裂安全系數(shù)比《重力壩》和《拱壩》的略小。在約束系數(shù)R的確定上，《重力壩》和《拱壩》的R值與基巖的彈性模量有關(guān)，在實(shí)際應(yīng)用中基巖的彈性模量大小一般較難確定，因此《水運(yùn)》參考了《塊體》中R值計(jì)算方法，計(jì)算結(jié)果也與工程實(shí)際吻合較好。

5 結(jié)論及建議

1)施工前應(yīng)用有限元軟件對(duì)大體積混凝土構(gòu)筑物進(jìn)行水化熱及溫度應(yīng)力分析，計(jì)算結(jié)果與實(shí)際工程監(jiān)測(cè)結(jié)果存在一定的誤差。在船閘的下一階段施工前，通過調(diào)整放熱函數(shù)和強(qiáng)度發(fā)展函數(shù)，采用經(jīng)過實(shí)踐驗(yàn)證的參數(shù)進(jìn)行分析，得到更符合實(shí)際情況的計(jì)算結(jié)果，為確定下一階段最佳的澆筑方法提供更準(zhǔn)確的參考依據(jù)。

2)通過施工前與施工第1階段導(dǎo)航墻后的有限元計(jì)算結(jié)果與實(shí)測(cè)結(jié)果的對(duì)比分析，介紹水化熱改進(jìn)計(jì)算的方法，該方法通過調(diào)整放熱函數(shù)和強(qiáng)度發(fā)展函數(shù)使水化熱的有限元計(jì)算結(jié)果更貼近工程實(shí)際情況。

3)通過裂縫指數(shù)i來評(píng)價(jià)大體積混凝土構(gòu)筑物溫度裂縫產(chǎn)生的可能性，施工前的分析結(jié)果表明，導(dǎo)航墻的表面產(chǎn)生溫度裂縫的可能性在60%以上，但實(shí)際上，在導(dǎo)航墻相應(yīng)位置的表面并未發(fā)現(xiàn)裂縫。表明在有限元計(jì)算相關(guān)參數(shù)的設(shè)置中可能存在偏差，在船閘施工第1階段后，通過調(diào)整放熱函數(shù)和強(qiáng)度發(fā)展函數(shù)進(jìn)一步計(jì)算，裂縫指數(shù)在1.8以上，裂縫產(chǎn)生的可能性在5%以下。

4)通過不同規(guī)范之間的對(duì)比，討論不同規(guī)范對(duì)大體積結(jié)構(gòu)溫度應(yīng)力計(jì)算規(guī)定的異同，分析不同規(guī)范對(duì)抗裂安全系數(shù)規(guī)定存在差別的原因以及相關(guān)計(jì)算參數(shù)取值的適用規(guī)范。

5)在本工程的導(dǎo)航墻中只埋設(shè)了溫度計(jì)，并未埋設(shè)應(yīng)力計(jì)和做回彈試驗(yàn)，無法得到混凝土實(shí)際的強(qiáng)度時(shí)程變化曲線來驗(yàn)證有限元的應(yīng)力計(jì)算結(jié)果，只能通過觀察結(jié)構(gòu)表面是否存在裂縫與軟件計(jì)算的溫度應(yīng)力進(jìn)行對(duì)比。建議下一步研究可據(jù)此做改進(jìn)。