ANSYS計算大體積混凝土溫度場的關(guān)鍵技術(shù).

1、ANSYS計算大體積混凝土溫度場的關(guān)鍵技術(shù)王新剛1高洪生2聞寶聯(lián)3(1中交天津港灣工程研究院有限公司結(jié)構(gòu)所，中國天津,300222;Z中交一航局 五公司，河北秦皇島,066002;3天津市市政工程研究院，天津,300074摘要:為了能夠用ANSYS較為準確地計算大體積混凝土溫度場，本文結(jié)合工程 實例，對單元選擇、網(wǎng)格劃分以及各時變參數(shù)的技術(shù)處理進行了分析和探索。計算 結(jié)果表明，溫度場的計算值與實測值比較接近，在變化趨勢上也基本一致，同時也驗證 了本文所提出方法的可行性。關(guān)鍵詞:大體積混凝土 ；水化熱;溫度場；ANSYSkey tech no logies calculate the mass

2、con crete temp erature field with ANSYSW ANG Xin-ga ng1GAO Hon g-she ng2WEN Bao-lia n2(1. Tianjin Port Engineering Institute, Ltd. of CCCC-First Harbor Eng. Co., Ltd.,Tianjin 300222; 2. No.5 En gi neeri ng Com pany Ltd. Of CCCC First Harbor Engin eeri ngCompany Ltd., Qinhuangdao 0660023 Tianjin Muni

3、cipal administration Research in stitute, Tianjin 300074Abstract: In order to calculate the mass con crete temp erature field accurately withANSYS, the an alysis and the expIo rati on has carried on to the unit select, the grid divisio n as well as each time variable eleme nt's tech ni cal p roc

4、ess ing in this article un ifies the p roject exa mp le. The compu ted result in dicated that the temp erature field p redicted value and the actual value are quite close, the cha nge tendency is also quite con siste nt, it also con firmed that the method feasibility which prop osed in this article

5、at the same time.Keywords:Mass con crete; hydrati on heat; Temp erature field; A NSYS0引言1大體積混凝土廣泛應(yīng)用于船塢、船閘、大壩等工程中。大體積混凝土由于體積 厚大,導(dǎo)熱系數(shù)較低，容易生產(chǎn)溫度裂縫。因此在設(shè)計、施工以及監(jiān)理階段需要詳細 驗算水化熱引起的溫度應(yīng)力。但由于水泥水化過程中，系統(tǒng)的溫度、生熱率、熱流率、熱邊界條件等參數(shù)隨時間都有明顯變化。利用大型通用有限元程序 ANSYS對大體積混凝土施工過程進行仿真分析可以同時考慮各參數(shù)的隨時間變化，并利用強大的后處理功能形象地給 出溫度場分布情況。但在分析過程中

6、單元的選擇、網(wǎng)格劃分以及環(huán)境溫度、水化熱生熱率、表面放熱系數(shù)等時變參數(shù)如何在程序中實現(xiàn)是仿真分析的關(guān)鍵技術(shù)，直接影響到溫度場計算的準確性。1單元選擇及網(wǎng)格劃分ANSYS的實體單元SOLID701有8個節(jié)點，且每個節(jié)點上只有一個溫度自由 度,具有三作者簡介:王新剛（1973-,男，河北玉田人，碩士,中交天津港灣研究院工程有限公司結(jié)構(gòu)所工程師。研究方向：橋梁施工控制、大體積混凝土溫度應(yīng)力計算及裂縫處 理。通訊地址：天津市河西區(qū)大沽南路1002號郵編:300222E-mai I:wxg58聯(lián)系電話:135*個方向的熱傳導(dǎo)能力，并能實現(xiàn)勻速熱流的傳遞。該單 元可以用于三維靜態(tài)或瞬態(tài)的熱分析，同時此單元

7、也可以進行結(jié)構(gòu)分析。用 ANSYS計算大體積混凝土溫度場的目的是以此為基礎(chǔ)來計算溫度應(yīng)力。因此計算大體積混凝土三維溫度場時可選取三維實體熱單元SOLID70,該單元可以在前處理器通過“ETCHG,TTS命令進行單元轉(zhuǎn)換，原來的熱單元SOLID70將自動轉(zhuǎn)換為結(jié)構(gòu)單元SOLID45,以方便接下來的溫度應(yīng)力計算。劃分網(wǎng)格是建立有限元模型的一個重要環(huán)節(jié)，網(wǎng)格數(shù)量的多少將影響計算結(jié)果 的精度和計算規(guī)模的大小。一般來講，網(wǎng)格數(shù)量增加，計算精度會有所提高，但同時計 算規(guī)模也會增加，所以在確定網(wǎng)格數(shù)量時應(yīng)權(quán)衡兩個因數(shù)綜合考慮。圖1計算精度和計算時間隨網(wǎng)格數(shù)量變化圖圖1中的曲線1表示計算精度隨網(wǎng)格數(shù)量變化情況

8、；曲線2表示計算時間度隨 網(wǎng)格數(shù)量變化情況。從圖中可以看出網(wǎng)格較少時增加網(wǎng)格數(shù)量可以使計算精度明顯提高,而計算時間不會有大的增加。當網(wǎng)格數(shù)量增加到一定程度后，再繼續(xù)增加網(wǎng)格時計算精度提高甚微，而計算時間卻有大幅度增加。所以應(yīng)注意增加網(wǎng)格的經(jīng)濟 性。實際應(yīng)用時可以比較兩種網(wǎng)格劃分的計算結(jié)果 ，如果兩次計算結(jié)果相差較大，可 以繼續(xù)增加網(wǎng)格，相反則停止增加。計算大體積混凝土溫度場時，在計算機硬件允許 的條件下應(yīng)取盡可能多的網(wǎng)格數(shù)量。2. ANSYS計算大體積混凝土溫度場各參數(shù)的技術(shù)處理2.1.環(huán)境溫度氣溫資料可以從附近的氣象站或水文站取得。在ANSYS中環(huán)境溫度變化可采用正弦或余弦表達式來模擬：(?

9、 ?-+=012cos TTn a ama A T T (2-1 式中:a T 氣溫;am T 日平均氣溫，可取日最高氣溫與最低氣溫的平均值;a A 氣溫日變幅，取最高氣溫與最低氣溫差值的一半；T寸間,h ;0 T氣溫最高的時間,h 02.2.水泥水化熱的施加水泥水化熱是影響混凝土溫度場的關(guān)鍵因素，推薦采用朱伯芳院士的復(fù)合指數(shù) 式2，其表達式為:1(0batQ t Q -= (2-2式中:Q（t混凝土水化熱；OQ 最終水化熱;b 系數(shù)。0Q和系數(shù)a、b應(yīng)根據(jù)現(xiàn)場施工混凝土配合比通過試驗來確定，也可根據(jù)現(xiàn)場實測數(shù)據(jù)反推得到。粗略計算時可按表1取值計算。水泥水化熱常數(shù)表1水泥L毗柚（MJ/呂）(J

10、b普通汗酸盤水泥!第W3300.69O.SS525#3500. 360. 71普通訕戲鹽人壩水泥S25B0. T9O.TO礦渣磚酸鹽大壩水泥S25t3S50. 296 76在ANSYS中，不能簡單地將計算出來的水化熱值作為邊界條件來施加,ANSYS的水化熱是通過生熱率HGEN來施加的。生熱率就是單位時間內(nèi)混凝土的生熱量，即所產(chǎn)生的熱量對時間的導(dǎo)數(shù)，其表達 式為:dttdQWHGEN(=(2-3式中:HGEN 混凝土水化生熱速率 W/m 3; W C 單位體積混凝土水泥用 量,kg/m 3;由式2-2、2-3可得：batatQ W t Q W HGEN -=1O'(2-4其APDL語言實

11、現(xiàn)如下：TIMEJ !載荷步終止時間 ESEL,S,MAT,2 !選擇單元 FLST,2,9000,2,ORDE,2FITEM,2,17501 FITEM,2,-26500*SET,HEAT,Wc*Q0*a*J*(b-1*EX P(-a*J*b!計算水化熱BFE,P51X,HGEN,_Z4, HEAT, , ,!施加水化熱 2.3.邊界條件在大體積混凝土結(jié)構(gòu)中，廣泛采用的是第三類邊界條件3,即假定經(jīng)過混凝土表 面的熱流量與混凝土表面溫度T和混凝土周圍介質(zhì)溫度aT之差成正比，即：(a nT T T 二？ p 入(2-5式中P固體表面的放熱系數(shù)，也稱對流系數(shù),C m W 02/?。其值與風(fēng)速有密切

12、的關(guān)系，而且混凝土表面的粗糙程度對其影響特別大，如果混凝土表面有保溫層，則表面放熱系數(shù)主要取決于各保溫層的厚度 和導(dǎo)熱系數(shù)，可按下式計算：刀+=q（2-5式中,iS模板及各種保溫材料的厚度（m;i入模板及各種保溫材料的導(dǎo)熱系數(shù)（C m W O ?/;q P3體在空氣中的放熱系數(shù),可取23(C m W 02/?。其 APDL 語言實現(xiàn)如下：*SET,Ta,Tam+Aa*cos(3.14/12*(JtO!計算環(huán)境溫度*SET,Beta,1/(D1/L1+D2/L2+Bq!計算放熱系數(shù)FITEM,2,316FITEM,2,20197 !選擇表面節(jié)點SF, P51X,C0NV, Beta, Ta !施

13、加放熱系數(shù)24其他參數(shù) 澆筑溫度澆筑溫度的模擬是將新澆混凝土的初始溫度設(shè)置為澆筑溫度，然后在水化放熱和熱邊界條件的作用下，其溫度場開始發(fā)展變化。導(dǎo)熱系數(shù)混凝土的導(dǎo)熱系數(shù)是 反映熱量在混凝土內(nèi)傳導(dǎo)難易程度的一個系數(shù)。普通混凝土的導(dǎo)熱系數(shù)一般在8.3912.56kJ/(m h 之間，具體值可由混凝土的各組成成分的重量百分比按加權(quán)平均方法計算得出 比熱單位質(zhì)量的混凝土，溫度升高1C時所需吸收的熱量稱比熱，普通混凝土的比熱 一般在 0.841.05kJ/(kg C之間。密度混凝土密度是指單位體積混凝土的質(zhì)量，普通混凝土的密度大約在23002450kg/m 3具體值可由混凝土配合比計算得出。另外，有關(guān)地

14、基的一些熱學(xué)參數(shù)，如導(dǎo)熱系數(shù)、比熱、密度等可由該工程地基勘察報告中獲得。3.施工階段的模擬在ANSYS中要實現(xiàn)按大體積混凝土施工階段計算溫度場，就要用到ANSYS的單元生死”技術(shù)。所謂單元的 殺死”，ANSYS程序并不是將 殺死”的單元從模型中刪除，而是將其剛度矩陣乘以一個很小的因子，因子缺省值為1.0E-6。死單兀的單兀載荷將為0,從而不對載荷向量生效。與上面的過 程相似，如果單元的 生”并不是將其加到模型中，而是重新激活它們。在計算大體積混凝土溫度場時，首先在PREP7中生成所有單元，然后全部殺死，再根據(jù)施工進度在合適的時刻逐塊激活相應(yīng)的塊體或單元并施加邊界條件，同時以循環(huán)過程來實現(xiàn)順序澆

15、筑過程的仿真，相應(yīng) 的 APDL語言如下：*DO,J,1,nESEL,S,MAT,2 !選擇全部混凝土單元EKILL,ALL !殺死全部所選取的單元ESEL,S,根據(jù)施工進度選擇單元 EALIVE,ALL !激活所選取的單元SOLVE*ENDDO4.工程應(yīng)用實例4.1工程概況 山海關(guān)船舶重工有限責(zé)任公司造船塢、舾裝碼頭工程由中交一航局五公司承建。船塢主體中的水泵房由泵房主體、集水池及沉砂槽等組成，采用天然巖基上的鋼筋砼箱形結(jié)構(gòu)，平面尺寸24.0 X9.75m，由上至下分別為:面層、電機層、水泵層、流道層。泵房南北外墻厚1.6m，東墻厚1.5m，西墻厚1.15m底板厚2.0m。集水池和沉砂槽伸入

16、船塢9.0m。泵房砼總方量4704m 3砼強度等級為C30。水泵房C30混凝土配合比見表2。水泵房C30混凝土配合比（kg表2水泥 p.042,5水202611536X2181誠木劑UEA粉煤址曠粉6.34934.362.453每層施工完成后停工20天，拆模后覆蓋一層塑料薄膜和兩層 4mm厚土工布進行養(yǎng)護。NODAL SOLUTIOHTEMP(AVG)ft3T3wO3HH -13,703SKX -4斗日7$圖2水泵房分層施工示意圖4.2建立模型及確定計算參數(shù)長取按水泵房的實際尺寸建立 ANSYS有限元模型，采用自由網(wǎng)格劃分方法,單元邊0.1m。根據(jù)山海關(guān)地區(qū)氣候資料確定 T。由表1和表2,Q取

17、 330J/g,a取 0.69。式 2-4 中 Wc=202+34.3+0.3*(62.4+53=270.92kg。則有:TIMEJ !載荷步終止時間=18+8*COS(0.262*(J-14HGEN=61688.5*J*(-0.44*EXP(-0.69*J*0.56 =9.9kJ/(m h :C 其他計算參數(shù)取值如表3所列。其他計算參數(shù)取值表350-45 -AlDIL匚HOLAL SCLUTLOBJF-'i T ” T M riO-*1416 ?O 50 碉45 -特性.比熱(j/g r)0.950. 34密度（血/W）21261570導(dǎo)熱系數(shù)(kj/ir-h 'C)12.55g. 27澆筑溫啦（匸）1228daj2 (Mpa)3,01.0泊檸比0.180. 20綜合以上各節(jié)分析，并要據(jù)施工進度形成APDL文件,計算水泵房混凝土溫度場， 每個施工階段混凝土澆筑后的前 3天每半小時計算溫度場一次，47天每小時計算一次，以每天計算一次。圖3水泵房一層第48小時溫度場水泵房一層第48小時溫度場如圖3所示。一層中心某點溫度計算和實測值隨 時間變化如圖4所示。溫度計算值與實測值最大誤差為 6.6%最小誤差為-2.1%。圖4水泵房一層中心某點溫度水泵房底板上表面某點溫度計算和實測值隨時間變化如圖5所示。圖5水泵房一層上