摘要:碾壓混凝土壩的非穩(wěn)定溫度場(chǎng)是進(jìn)行仿真應(yīng)力計(jì)算和設(shè)計(jì)的基礎(chǔ)。根據(jù)理論分析和數(shù)值計(jì)算,分析各熱學(xué)參數(shù)和非熱學(xué)參數(shù)對(duì)非穩(wěn)定溫度場(chǎng)的影響。通過(guò)對(duì)石門子工程(在建)進(jìn)行工程實(shí)際監(jiān)測(cè),取得了一些高程測(cè)點(diǎn)實(shí)測(cè)溫度變化過(guò)程,并與仿真計(jì)算結(jié)果相比較。為模擬混凝土中粉煤灰后期放熱,提出雙e 曲線模型,使仿真計(jì)算結(jié)果更符合實(shí)際。按所提出的方法進(jìn)行非穩(wěn)定溫度場(chǎng)的計(jì)算預(yù)測(cè),從而對(duì)大壩后期澆筑作指導(dǎo),具有較好的工程應(yīng)用價(jià)值。
關(guān)鍵詞:碾壓混凝土壩;非穩(wěn)定溫度場(chǎng);工程實(shí)測(cè);
塔西河水利樞紐工程位于新疆維吾爾自治區(qū)昌吉州瑪納斯縣西南塔西河中游“石門”峽谷段,峽谷長(zhǎng)約350 m ,河谷呈“U ”型,谷底寬80 m,布置109m 高的碾壓混凝土拱壩。
壩址地區(qū)冬季寒冷,月平均氣溫在0 ℃ 以下長(zhǎng)達(dá)5 個(gè)月,夏季氣溫最高達(dá)38 ℃ ;工地日變幅常達(dá)
15～ 19 ℃ 。這不僅縮短了砼的良好澆筑時(shí)間,而且還給混凝土溫控帶來(lái)較大困難。
面對(duì)惡劣的自然條件,在快速施工方面,用低熱混凝土連續(xù)大倉(cāng)面澆筑兼用保溫與降溫措施;在結(jié)
構(gòu)方面,在壩肩設(shè)人工短縫、拱冠設(shè)中縫以釋放拱壩溫度應(yīng)力,并在拱壩上游早期降溫區(qū)設(shè)“鉸點(diǎn)”,使拱壩在拱向溫度收縮的同時(shí)能保證平面拱向整體作用[1,2]。為提前蓄水提供安全保障。
1 參數(shù)分析
影響壩體溫度場(chǎng)的因素有很多,其中主要有:混凝土澆筑層厚度、邊界氣溫及澆筑間隙期內(nèi)部溫控措施、混凝土入倉(cāng)溫度、混凝土絕熱溫升θ0 及混凝土放熱系數(shù)m 等,見(jiàn)表1[3 ]。各部分參數(shù)作用耦合在一起。
1. 1 混凝土厚度的影響
溫度實(shí)測(cè)表明,邊界的影響在1m 范圍之內(nèi)較顯著。如果澆筑間歇期小于6～ 10 d,超過(guò)2 m 深度處的混凝土溫度受外界影響很小,所以單從溫度場(chǎng)來(lái)看,在澆筑過(guò)程中,將每次累積層厚度定為小于2m,對(duì)散熱有利。為了減少層間間隔,加快施工進(jìn)度及考慮到模板工藝,石門子工程采用3m 澆筑層。
1. 2 邊界條件
邊界條件變化有兩種,散熱系數(shù)的變化以及邊界氣溫值的變化。第1 種變化影響較大,可以說(shuō)是改 變了熱交換條件,如固氣交換變成固液交換等;第2種情況較普遍,如周圍空氣溫度上升、下降等。
1. 2. 1 散熱系數(shù)變化的影響
本文比較的兩種散熱條件,一種是用河水養(yǎng)護(hù)(比初始邊界氣溫低8 ℃ ,初始邊界氣溫15 ℃ ),一
種是靠邊界空氣散熱(散熱系數(shù)β 取0. 8,初始邊界氣溫15 ℃ ,氣溫降溫幅度為- 1 ℃ / d)。計(jì)算條件:入倉(cāng)溫度15 ℃ ,水化熱溫升17 ℃ ,混凝土放熱系數(shù)0. 16。計(jì)算結(jié)果為:12 d 后,在1 m 深處,河水冷卻條件下,混凝土溫度要多降低4 ℃ 左右。
1. 2. 2 氣溫值變化
石門子工程所在山區(qū),氣溫日變幅大。圖1為± 8,9 ℃ 時(shí)混凝土溫度場(chǎng)的變化圖,可以看出日變
幅對(duì)后期混凝土溫度場(chǎng)影響較小,1 m 深處混凝土只有0. 2 ℃ 的變幅,即使在0. 5 m 深處,也只有0.5℃ 左右的變化。不管是白天,還是夜晚對(duì)最后結(jié)果影響不大。但由于夜間周圍氣溫較低,對(duì)骨料散熱有利;再加上在運(yùn)輸路上及入倉(cāng)后的混凝土不受太陽(yáng)直射,所以其入倉(cāng)溫度比白天要低,也就降低了混凝土澆筑塊溫度。所以在實(shí)際工程施工中,一般在夜間澆筑混凝土。
1. 3 入倉(cāng)溫度
其他條件不變,入倉(cāng)溫度提高+1℃,6d后在1m 深處,仍維持升高+0.8℃;12d 后,仍有+0.6℃ 。而在周圍氣溫變化不大的情況下,通過(guò)澆筑面散熱,6d才散熱1.4℃ 。也就是說(shuō),只要入倉(cāng)溫度高2 ℃ ,則6d 散熱的效果就會(huì)被抵消。石門子工程屬于快速施工間歇期較短的情況,因此必須注意混凝土入倉(cāng)溫度的控制。
1. 4 混凝土的絕熱溫升
由于混凝土絕熱溫升曲線常用表達(dá)式為θ=θ0 (1- e- m t),絕熱溫升θ0 的影響是成比例關(guān)系。若放熱系數(shù)不變,取0. 16,邊界為固氣熱交換,則絕熱溫升增加+ 1 ℃ ,在混凝土3 m 深處,混凝土溫度值增加+ 1 ℃ 左右;1m深處6d為+ 0.5℃ ,12d 為+0.60℃ 。
通常水泥絕熱溫升e 曲線擬合,曲線前期值上升較快,后期(一個(gè)月左右)持平。而在工程中用的粉
煤灰后期水化發(fā)熱,擬合采用雙e 曲線更好。圖2是我國(guó)南方某一工程混凝土絕熱溫升測(cè)試結(jié)果及曲線 模擬,可以看出采用雙e 曲線更能符合實(shí)際情況[4]。
1. 5 放熱系數(shù)
為了研究混凝土放熱系數(shù)對(duì)混凝土溫度的影響,在計(jì)算中,放熱系數(shù)m 分別取0. 16 和0. 32。其他計(jì)算條件相同,即邊界散熱系數(shù)取0.8,澆筑塊入倉(cāng)溫度0 ℃ ,絕熱溫升17 ℃ 。計(jì)算表明混凝土的放 熱系數(shù)大時(shí),早期放熱快,所以在一定間歇期內(nèi)散出熱量也多,12 d 后,1 m 深處,放熱系數(shù)0. 32 要比放熱系數(shù)0. 16 散熱多,達(dá)1. 3 ℃ 。
1. 6 層間間隙時(shí)間的影響
邊界氣溫升溫、降溫速率假設(shè)為±1℃/d。表2中所列溫度下降值,均是相對(duì)于混凝土在絕熱溫升 條件下正常溫升的差值。
由表2 可以看出,如果外界日平均氣溫上升,又無(wú)冷卻措施,6d 后澆筑第2層和12d后澆筑第2層,兩者相差不大,不到1 ℃ 左右。所以對(duì)這種情況間歇期的降溫效果不大。
但如果日平均氣溫不變或降低,12d后要比6d要散熱效果好得多,即以時(shí)間代價(jià)換取混凝土溫度值的降低。
如果又想溫度下降快,又不想時(shí)間擱得太長(zhǎng),就只能用河水冷卻或其他冷卻方法。為減少層間缺陷,簡(jiǎn)化工藝,加高模板,加快施工進(jìn)度,石門子工程采用3 m 澆筑砼,間歇期間散熱采用低溫河水浸泡方式。
2 三維仿真數(shù)值計(jì)算
圖3是石門子碾壓混凝土拱壩1295高程溫度計(jì)埋設(shè)位置圖。仿真計(jì)算所用邊界溫度值、混凝土入倉(cāng)溫度以及材料熱學(xué)參數(shù)均是實(shí)測(cè)資料。因?yàn)樵趯?shí)測(cè)數(shù)據(jù)時(shí)間段內(nèi),大壩沒(méi)有蓄水,所以計(jì)算中不考慮大壩滲流的影響。
根據(jù)石門子拱壩實(shí)際開(kāi)挖與澆筑的體形,共計(jì)離散為3 352 個(gè)六面體或五面體(過(guò)渡)等參單元,
4580個(gè)節(jié)點(diǎn),計(jì)算得到施工過(guò)程仿真非穩(wěn)定溫度場(chǎng)的累計(jì)變化過(guò)程[5,6]。
2. 1 邊界條件
施工期從下往上逐層澆注混凝土?xí)r,上下游邊界暴露于大氣中。拱壩上游面基坑有上游滲水,所以 邊界條件為第一類邊界條件。由于水是上游天山雪水,所以水溫度較低(7 ℃ )且變化值不大。拱壩頂部暴露面隨澆筑高度上升,依靠自然冷卻的碾壓混凝土塊的溫度隨氣溫及輻射熱的變化而變化。在早期壩澆筑過(guò)程中,為了施工方便,用土堆了上壩道路,如圖4 所示。土對(duì)混凝土散熱產(chǎn)生了一定負(fù)面影響,但在嚴(yán)冬中可起一定保護(hù)作用。
2. 2 混凝土入倉(cāng)溫度
加高混凝土,實(shí)際上可以說(shuō)是一個(gè)新過(guò)程的開(kāi)始。對(duì)老混凝土來(lái)說(shuō),只是改變了邊界條件,其初始
條件為:老混凝土的當(dāng)時(shí)溫度作為初溫,還有剩余的沒(méi)有散發(fā)的水化熱,也就是在水化熱方程中時(shí)間參數(shù)要減去一個(gè)初值;對(duì)于新混凝土,其初始溫度便是澆筑溫度。所以在有限元計(jì)算時(shí),用單元溫度, 即用單元離散混凝土的溫度[3,4]。
2. 3 混凝土熱學(xué)參數(shù)
由于抗?jié)B性、抗凍耐久性等要求不同,石門子拱壩各部位混凝土配合比也不盡相同。大壩內(nèi)部三級(jí)配混凝土含水泥60 kg,粉煤灰110kg,θ0 為17 ℃ , m 為0. 16。外部抗?jié)B二級(jí)配混凝土含水泥90 kg, 粉煤灰110 kg,θ0 為23 ℃ ,m 為0. 32。其中熱學(xué)參數(shù)由本課題組研制的混凝土絕熱溫升測(cè)試儀(同溫絕熱)測(cè)出。該儀器可按照混凝土試驗(yàn)規(guī)程,進(jìn)行長(zhǎng)時(shí)間混凝土絕熱溫升曲線跟蹤測(cè)試。
3 計(jì)算預(yù)測(cè)值與工程實(shí)測(cè)值的比較
計(jì)算和實(shí)測(cè)的溫度變化規(guī)律見(jiàn)圖4(a)～ (f)。
1)計(jì)算預(yù)測(cè)與實(shí)測(cè)的溫度過(guò)程線的曲線走勢(shì)和值大小基本一致。走勢(shì)均有兩個(gè)高峰期,第一個(gè)出
現(xiàn)在上面混凝土覆蓋之前,第二次出現(xiàn)在覆蓋之后。對(duì)于覆蓋之后,離上下游散熱面較近的測(cè)點(diǎn)1 和測(cè) 點(diǎn)5 的第二個(gè)高峰期的溫度比第一次要低,而對(duì)于中間散熱條件較差的點(diǎn),情況正好相反,第二次高峰
期溫度要比第一次高。測(cè)點(diǎn)2由于離廊道較近,所以其后期溫度沒(méi)有3,4測(cè)點(diǎn)上升的高;5號(hào)測(cè)點(diǎn)因?yàn)?后期有土覆蓋,所以其溫度變化形式與1號(hào)測(cè)點(diǎn)也不一致。
4 結(jié)論
本文經(jīng)過(guò)理論分析和數(shù)值計(jì)算,分析了各熱學(xué)參數(shù)和非熱學(xué)參數(shù)對(duì)非穩(wěn)定溫度場(chǎng)的影響。對(duì)比實(shí)
測(cè)和仿真計(jì)算結(jié)果,說(shuō)明決定計(jì)算預(yù)測(cè)可靠性的關(guān)鍵是熱學(xué)參數(shù)(入倉(cāng)溫度、絕熱溫升和邊界溫度等)
的正確取值。除需獲得工地實(shí)際氣溫變化情況和施工進(jìn)度,還要正確測(cè)試和模擬大壩混凝土的絕熱溫升曲線。文中為了模擬碾壓混凝土中粉煤灰后期放熱,采用雙e 曲線模型,取得較好效果。在大壩設(shè)計(jì) 階段,能獲得準(zhǔn)確的仿真溫度場(chǎng),對(duì)研究壩體應(yīng)力分布,調(diào)整施工計(jì)劃及采取溫控措施要求,都是很有實(shí)
際工程價(jià)值的。
參考文獻(xiàn)
[1] Dunstan M R H. Recent developments in RCC dams [J]. Int J H yd ropow
er & Dams,1999,6(1):40 45.
[2] 張國(guó)新,劉光廷,劉志輝. 整體碾壓混凝土拱壩工藝及溫度場(chǎng)仿真計(jì)算[J]. 清華大學(xué)學(xué)報(bào),1996,36(1):1 7.
[3] 劉光廷,邱德隆. 含橫縫碾壓混凝土拱壩的變形和應(yīng)力重分布[J]. 清華大學(xué)學(xué)報(bào),1996,36(1):20 26.
[4] 張富德,薛慧,陳鳳歧,等. 碾壓混凝土壩非穩(wěn)定溫度場(chǎng)計(jì)算預(yù)測(cè)與工程實(shí)踐的比較[J]. 清華大學(xué)學(xué)報(bào),1998, 38(1):75 78.
[5] 劉光廷,麥家煊,張國(guó)新. 溪柄碾壓混凝土薄拱壩的研究[J]. 水力發(fā)電學(xué)報(bào),1997,57(2):19 28.
[6] ZHEN G G R,Y AN G S C. Solution of steady periodic heat conduction by the
FEM [J]. N umerical H eat T ransf er, 1989,15:525 534.
